CC BY
9DOI: 10.17516/1999-494X-0312 УДК 621.317.335.3:551.578.465
Methods for Assessing the State of Snow-Ice Cover
Viktor G. Mashkov*, Vladimir A. Malyshev and Pavel A. Fedyunin
Military Education and Research Centre of Military-Air Forces
«Military-Air Academy Named After Professor N. E. Zhukovsky and Yu. A. Gagarin»
Voronezh, Russian Federation
Received 13.03.2021, received in revised form 15.04.2021, accepted 21.05.2021
Abstract. The development new methods for measuring the electrophysical and geometric parameters multilayer dielectric medium in order to identify their layers, as well as to detect inhomogeneities in them, is an urgent scientific task. Methods for assessing the state the snow-ice cover are proposed, based on the identification layers of the snow-ice cover by the dielectric permittivity obtained indirectly as a result inclined sounding of the underlying surface by an electromagnetic wave with vertical and horizontal polarization. Under normal sensing the underlying surface, in the form snow-ice cover, by radio waves with vertical and horizontal polarization, the obtained values the reflection coefficient from the boundary layers the received signals are identical. Oblique sensing by a vertically polarized wave in the range from 25 to 90 degrees leads to a sharp dip in the values the reflection coefficient from the boundary layers at certain angles, so for dry snow - 47...55°, dry firn - 55 ... 58° and dry ice - 58.61°. These angles characterize the dielectric properties the layers (permittivity), which are explained by the complete polarization the reflected wave in a plane perpendicular to the plane incidence, since at this moment the vibrations the layer's electrons along the direction their movement do not lead to the emission electromagnetic waves in the direction the reflected wave (absent). The given estimated values the complex relative permittivity snow-ice cover as two component medium air and ice at a temperature from minus one to minus forty degrees, or three component medium air, ice and water at zero temperature are a sample for establishing the identity the characteristic properties layers. The purpose this article is to develop methods for assessing the state snow-ice cover used in determining the possibility performing a safe landing a helicopter - type aircraft on an unprepared site with snow-ice cover based on the identification the obtained characteristics snow-ice cover layers based on the results radar sensing with calculated data.
Keywords: state snow-ice cover, underlying surface, permittivity, layer identification, oblique sensing.
Citation: Mashkov V.G., Malyshev V.A., Fedyunin P.A. Methods for assessing the state of snow-ice cover, J. Sib. Fed. Univ. Eng. & Technol., 2021, 14(3), 316-331. DOI: 10.17516/1999-494X-0312
© Siberian Federal University. All rights reserved
This work is licensed under a Creative Commons Attribution-Non Commercial 4.0 International License (CC BY-NC 4.0). Corresponding author E-mail address: [email protected]
Методы оценки
состояния снежно-ледяного покрова
В. Г. Машков, В. А. Малышев, П. А. Федюнин
Военный учебно-научный центр ВВС «Военно-воздушная академия имени профессора Н. Е. Жуковского и Ю. А. Гагарина» Российская Федерация, Воронеж
Аннотация. Предложены методы оценки состояния снежно-ледяного покрова, основанные на идентификации слоев снежно-ледяного покрова по диэлектрической проницаемости, полученной косвенным путем в результате наклонного зондирования подстилающей поверхности электромагнитной волной с вертикальной и горизонтальной поляризацией. При нормальном зондировании подстилающей поверхности в виде снежно-ледяного покрова радиоволнами с вертикальной и горизонтальной поляризацией полученные значения коэффициента отражения от границ раздела слоев принятых сигналов идентичны. Наклонное зондирование вертикально поляризованной волной в пределах от 25 до 90° приводит к резкому провалу значений коэффициента отражения от границ раздела слоев при определенных углах: так, для сухого снега - 47...55°, сухого фирна - 55...58° и сухого льда - 58...61° Эти углы характеризуют диэлектрические свойства слоев (диэлектрическую проницаемость), объясняемые полной поляризацией отраженной волны в плоскости, перпендикулярной плоскости падения, поскольку в этот момент колебания электронов слоя вдоль направления своего движения не приводят к излучению электромагнитных волн в направлении отраженной волны (отсутствует). Приведенные оценочные значения комплексной относительной диэлектрической проницаемости снежно-ледяного покрова как двухкомпонентной среды воздуха и льда при температуре от минус 1 до минус 40° либо трехкомпонентной среды воздуха, льда и воды при нулевой температуре являются образцом для установления тождественности характеризующих свойств слоев.
Цель данной статьи - разработка методов оценки состояния снежно-ледяного покрова, используемых в определении возможности выполнения безопасной посадки воздушного судна вертолетного типа на неподготовленную площадку со снежно-ледяным покровом, на основе идентификации полученных характеристик слоев снежно-ледяного покрова по результатам радиолокационного зондирования с расчетными данными.
Ключевые слова: состояние снежно-ледяного покрова, подстилающая поверхность, диэлектрическая проницаемость, идентификация слоев, наклонное зондирование.
Цитирование: Машков, В. Г. Методы оценки состояния снежно-ледяного покрова / В. Г. Машков, В. А. Малышев, П. А. Федюнин // Журн. Сиб. федер. ун-та. Техника и технологии, 2021, 14(3). С. 316-331. DOI: 10.17516/1999-494Х-0312
Introduction
Remote sensing the underlying surface the landing site a helicopter-type aircraft (HTA), in the form snow-ice cover, which is a multi-layer structure, to assess its condition, is almost the only justified method, from the point view performing a safe landing on unprepared (unequipped) sites. The task assessing the state snow-ice cover is reduced, first, to the identification its layers (snow, firn, ice), since the number layers snow-ice cover formed in the natural environment will be different, which is due to their formation during snowfalls, thaws, cold spells, etc., determined by a complex hydrometeorological conditions. Secondly, to determine the depth snow and the thickness the ice cover, which is a key factor
in making a decision to land a HTA on an unprepared (unequipped) snow-ice platform. Landing on the snow-covered pond with snow depth above the permissible or thickness the ice cover below the permissible in accordance with the type HTA, may lead to breaking through the snow, ice, or roll HTA day and night in simple and adverse weather conditions (fog, haze, rain, snow, dust, or smoke the atmosphere), and also in conditions raised snow rotating propeller [1, 2].
Method for estimating the state snow-ice cover by Brewster angle
Identification snow-ice cover layers. Measurement complex relative permittivity materials with losses in the microwave range [3-5], based on measuring the dependence the coefficient of Fresnel reflection a flat sample the material from the angle incidence within the range 6 = 40°.. .90° when the electric field the wave parallel to the plane incidence, then the graph the dependence the coefficient of Fresnel reflection from the angle incidence determine the Brewster angle 6B and the relative complex dielectric permittivity s the sample material or portion the surface the earth is calculated by the formula
s = (tg6e)2-
Since the snow-ice cover is a multi-layer structure with different densities, water content shares, and its own structures, it is important to determine the boundaries the air-snow, snow-ice, and ice-water layers. It is not always possible to obtain peaks echo signals from the boundary snow-i ce layers as a result the contrast the complex relative permittivity the layers, since if the snow is inhomogeneous, compacting to the bottom so that at the s-i boundary es ~ s;, then there is no reflection from the s-i boundary at all [6]. But, as a rule, the snow cover has layers with various complex relative permittivity, increasing in depth, which is due to their compaction during snowfall, thawing, cooling, etc., which allows you to receive echoes from the boundary layers snow-ice cover. Moreover, for example, for freshwater ice (at X = 3 cm), the amplitude the echo signal from the i-w border will be greater than from the s-i border, this is observed when the thickness the ice cover is less than 1.6 m [6].
On the territory the Russian Federation, the thickness the ice cover on rivers and lakes varies from 0.1 m in the southern regions to 3 m in the polar zone Siberia and up to 4 m on the ice sections rivers in Eastern Siberia. In most cases, and especially in highly salty sea ice, the thickness does not exceed 2 m [6].
The dependence the reflection coefficient of the Fresnel with vertical polarization on the complex relative permittivity the medium er and the sensing angle within the limits 6 = 40°.. .90° has a pronounced dip in values (Fig. 1a) compared to the dependence the reflection coefficient of the Fresnel with horizontal polarization Rhh (Fig. 1b). For vertical and horizontal polarization (vv, hh - the first index the probing, the second received radio signal).
The reflection coefficient from the boundary snow-ice layers is determined by [4, 7] the sensing angle and the relative permittivity the corresponding layer:
1) w _
1,2 "
where 6i - sensing angle;
nw _ 2,3 "
, where 62 = arcsin (sin 6J^s^ j;
^2 COS 6 '„1 (^2 "^rl1 (sin 0
[^2 COS 6 <sJ (^2 "^rl1 (sin 0
cos [si3-£i2( Sin6>2)2) 0
U3COS#2+^ (£r3-£r2(sin#2)2) 0
a) b)
Fig. 1. Dependence the Fresnel reflection coefficient: a - with vertical polarization Rvv; b - with horizontal polarization Rhh on the complex relative permittivity the medium and the sensing angle within 0 = 0° .. 90°
-,-m: COSé/-
£d+l COS 6> +Jsn (snl] - £u (sin Qt )2
, where di = arcsin (sin oj Jc
(1)
The total reflection coefficient (1) from the snow-ice cover without taking into account multiple reflections between the layer boundaries is determined by the recurrent formula [4]:
r:
v(hh) _
R
(2)
.... " i + RJJ^R^ exp(-jAxhJI
where hvi+1 - is the depth (/ + l)-layer; X - is the wavelength the probing signal; ]VJ(hh> = o, k /, k i + 1.
The speed propagation an electromagnetic wave (EW) in a snow-ice cover is lower than the speed propagation in the air, depending on its density, water content, and structure. For dry snow Vds = 278,1...212,7 mVs, dry firn Vdf = 212,7...189,0 m/^s and dry ice Yd, = 189,0...167,9 m/^s, vary very markedly depending on the proportion water content, since for pure water Vpw = 32,74.41,97 m/^s, the preferred orientation and shape ice and air inclusions in snow [2], so it is necessary to take into account the wavelength in the layer.
Given the change in the wavelength in the layers, in formula (2) we obtain:
An h
where y] = -j r2
(/ + 1)-layer's; Anll =
- relative depth
It is known that if the probing signal falls on the boundary two dielectrics at an angle 0i equal to the Brewster angle 0B, then the reflected 01 and refracted 02 signals are perpendicular to 01 + 02 = 90°.
In this case, the reflected signal will be completely polarized in the plane perpendicular to the plane incidence (absent), and the degree polarization the refracted beam will be maximum. According to Snellius law (3), for a wave with the angle incidence 61 and the angle 62 refraction at the boundary two dielectrics, the equality is valid for the medium under study (with ^ = 1: snow, firn, ice):
sK (3)
where n2, n1 - is the refractive index the first and second layers the medium; n21 - the refractive index the boundary the first and second layers; e1 and e2 - the relative permittivity the first and second layers, respectively.
The Brewster angle is defined:
(
dy, = arctg(//2 j) = arctg
The Brewster angle for the z-th layer boundary, at which there is no reflected signal with vertical polarization from the z-th layer boundary, will be determined
6b 1, z+1 = arctg(nmj), where nz+11 = ni+1 / n1 - the refractive index the z-th layer boundary, then
L | = arctg
I "i
= arctgI ^ I = arctg
(4)
the sensing angle the first layer the underlying surface 6B 1, z+1 ^ 6B 1, k. The refractive z-th layer boundary increases in depth with increasing snow cover density.
According to the proposed method, probing the radio signal with vertical polarization the controlled area the snow-ice cover within the limits 6 = 40°. ..90°, allows us to determine the dependence the Fresnel reflection coefficients R™+1, (',/ + l)-boundary layers on the angle incidence the radio wave.
Fig. 2 shows the Fresnel reflection coefficients RJJ+1 for layers with complex relative permittivity: er2 = 1,3 - j0,0008, sr3 = 1,8 - j0,0008, er4 = 2,3 - j0,0008, sr5 = 2,8 - j0,0008, er6 = 3,1 - j0,0008, er7 = 74 - j, depending on the sensing angle within 6 = 40°.. .90°. The increase in electrical conductivity for seawater (er7 = 74 - j) is due to the presence salinity Ssw = 35 g/kg.
A sequential increase in the sensing angle, starting from 6 = 40°, will lead to a sequential decrease (disappearance) each subsequent peak the echo signal from the subsequent layer boundary at the appropriate 6, since the density the snow-ice cover increases as the depth increases, ranging from pr = 10 kg/m3 - for freshly fallen snow to pr = 917 kg/m3 - for dry ice without air inclusions.
Thus, the intervals the snow-ice cover densities will correspond to the intervals of Brewster angles. So for dry snow pds = 100.500 kg/m3 (s'ds = 1,162.1,984) - 6Bds = 47°.55°, dry firn pdf = 500...700 kg/m3 (e'df = 1,984.2,51) - 6Bdf = 55°.58°, dry ice pdi = 700..917 kg/m3 (sdi = 2,51.3,179) - 6 Bdi = 58°.. .61°. As the moisture content increases, the values of Brewster angles shift to the values for water. For sea water ssw = 74 salinity is Ssw = 35 g/kg, and for pure (meltwater) water epw = 87 with - 6 Bsw = 84°.
Fig. 2. Fresnel reflection coefficients RJ^+i for layers with complex relative permittivity: er2 = 1,3 - j0,0008, er3 = 1,8 - j0,0008, er4 = 2,3 - j0,0008, er5 = 2,8 - j0,0008, er6 = 3,1 - j0,0008, er7 = 74 - j depending on the sensing angle within 0 = 40°... 90°
The graphs (Fig. 2) show the absence reflected signals from the boundary layers at 0, for medium with complex relative permittivity: er2 = 1,3 - j0,0008, er3 = 1,8 - j0,0008, er4 = 2,3 - j0,0008, er5 = 2,8 - j0,0008, er6 = 3,1 - j0,0008, er7 = 74 - j depending on the sensing angle within the range 0 = 40° ..90°, respectively, equal to the Brewster angles 0BM+1: 0B1,2 = 49°, 0B2,3 = 53°, 0B3,4 = 57°, 0B4,5 = 59°, 0B5,6 = 60°, 0B6,7 = 83°.
The solution the inverse problem the reconstruction the layers. The essence solving the inverse problem estimating the dielectric permittivity individual layers snow-ice cover based on the determination of Brewster angles is as follows. Let the object sensing the underlying surface is a snow-ice cover of k layers irradiated by a linear frequency modulated (LFM) signal at an angle of 40 to 90 degrees by an electromagnetic wave with vertical polarization.
Brewster angles are determined using the minimum measured Fresnel reflection coefficient for i,i + 1 the boundary between layers snow and ice cover 0B!-,!+i. At the same time, the Brewster angles 0Bi,i+1 are determined sequentially, which allows us to determine the permittivity each subsequent layer from the formula (4):
eri+1 = (tg0B1,i+1)2. (5)
The methodological error in determining the dielectric permittivity layers (Fig. 2) at Brewster angles 0Bi,i+1 at specified design values eri+1 (for example, obtained 0B1,2 = 49°, 0B2,3 = 53°, 0B34 = 57°, 0B4,5 = 59°, 0B5,6 = 60°, 0B6,7 = 83°) constructed according to formula (1) and determined according to formula (5), that is, their assessment without taking into account the impact of noise is no more than 3 %, which indicates the validity using formula (5).
To estimate errors in measuring the permittivity layers corresponding to the actual practical situation, an additive Gaussian noise was added to the values of Brewster angles 0Bii+1
0Bi,i+1exp = 0Bi,i+1 + n(f),
where n(f^) - is the noise affecting £ measurement, due to both the error the measurement complex and other experimental errors during the experiment, and the reliability the selected electrodynamic model the real physical situation.
With this in mind, formula (5) will take the form
Srz+1exp = (tg6B1,/+1)2 + n(f).
Fig. 3 shows the dependence the measured values the relative permittivity sn+1exp the k layers (Fig. 3a) and the error their estimation (Fig. 3b) when the values RMS the noise level increase from 3.8 to 4.8 in increments of 0.1 for 100 implementations each, where 1 - s^exp, 2 - s^exp, 3 - s^xp, 4 - s^xp,
5 sr6exp-
The analysis dependencies in Fig. 3 shows that the value errors in estimating the measured values the relative permittivity sn+1exp of k layers when the values RMS the noise level increase from 3.8 to 4.8 in increments of 0.1 for 100 implementations each with a confidence probability of 0.95 does not exceed 10 %, which indicates the validity using this method.
The obtained values permittivity layers sn+1 (z + 1)-layers are identified with the sample values dielectric permeability layers eVrA and the estimated state the snow-ice condition sn+1 = sVrA: ''snow cover'', ''firn'', ''ice cover'' or ''water''.
At negative temperatures T = -1.-40 °C, the real part the complex relative permittivity layers er S - snow, sf - firn, si - ice) with the density layers pr = 100... 917 kg/m3 (p = 100... 500 kg/m3; Pf = 500.700 kg/m3; pi = 700. 917 kg/m3) do not depend on f = 1.10 GHz, but only on T in small limits. The real part the sample values relative permittivity wet medium ev^ (for the General formula a three-component medium): s'vsa - snow; svfA - firn;
For example, T = -1. -40 °C for snow, as a two-component medium air and ice, the actual part the permittivity s'vsa is between air evaA = 1 and ice sviA = 3,15 (for dry dense ice without air inclusions with Pi = 917 kg/m3 - sviA = 3,20±0,02), for firn (densely Packed, granular and partially recrystallized, usually multi-year snow, i. e. an intermediate stage between snow and glacier ice), the relative permittivity approaches the values for ice. For water under a snow-ice cover at T = 0 °C frequencies f = 2.8 GHz,
1.5
0-..+,..........o------O.....0-—-0.....o-..
3.8 3.9 4.0 4.1 4.2 4.3 4.4 4.5 4.6 4.7 4.8 a 3.8 3.9 4.0 4.1 4.2 4.3 4.4 4.5 4.6 4.7 4.8 c
a) b)
Fig. 3. Dependence: a - measured values the relative permittivity sn+1exp of k layers; b - errors in their estimation when the values RMS the noise level increase from 3.8 to 4.8 in increments of 0.1 for 100 implementations each
there is a smooth decrease in the actual part the dielectric constant for meltwater epwA = 83,84.51,03, for sea water with a salinity of Ssw = 35 g/kg - e'swA = 74,97.48,42 determined by the Debaev model [2].
The values the real part the dielectric constant the medium is calculated using the formula Of G. Lueng [2, 10-12] for such medium with spherical air and ice inclusions
^(^(Vi-lj + l)3, (6)
where = pr / p; - the volume content ice; pr - the density the dry medium (dry snow pds, dry firn pdf, dry ice pdi); p; = 917 kg/m3 - the density dry ice without air inclusions; ej - the actual part the dielectric permittivity ice.
At temperature T = 0 °C, a significant contribution to the permittivity will be made by humidity (the proportion water content in the layer) for a two-component ice-water mixture with water-filled pores
+ (7)
where Pw - is the total percentage water content; ej and sw - the relative permittivity ice and water, respectively. General formula for a three-component medium [2, 10-12] consisting ice with water and air inclusions including special cases (6) and (7)
=(V? (i))3, (8)
where Pwa - is the total percentage water and air content.
Fig. 4 shows the dependence the real part the complex relative permittivity the medium (snow, firn, ice) on the density the medium pr (6) and the proportion water content (7) Pw = 0.0,1 for f = 2 GHz (Fig. 4a), for f = 8 GHz (Fig. 4b) when Tr = 0 °C.
Thus, based on the dielectric constant obtained indirectly as a result inclined sensing the underlying EW surface with vertical polarization, it is possible to restore the parameters the layers snow-ice cover (density pr, percentage water content Pw).
a) b)
Fig. 4. Dependence the real part the complex relative permittivity the medium s'vr (snow, firn, ice) on the density the medium pr (6) and the proportion water content Pw = 0.. .0,1 (7) for: a - f = 2 GHz; b -f = 8 GHz when Tr = 0 °C
The depth snow-ice cover layers is measured as follows.
The geometrical parameters the layers the snow-ice cover. Probing a controlled section snow-ice cover with a linear frequency modulated (LFM) signal [6, 13] that implements the frequency principle determining the depth layers snow-ice cover, based on the selection the beat frequency (difference signal) obtained by multiplying the received and probing (reference) signals. Generated in the transmitting device LFM signal with frequency
ftx(t) = f0 + atm, V 0 < tm < Tm, (9)
where f = 2 GHz - is the starting frequency, a = 600 GHz/s - the rate change frequency (LFM slope), tm - the time within a single period modulation chirp signal (fast time), a Tm = 10 ms - the period modulation (chirp LFM signal) and receiving echo signals by a receiver with a frequency
frx(t) = f + a(tm - t), V t < tm < Tm, (10)
collected according to the classical scheme. At the same time, it is necessary to note the requirement for the formation a LFM signal in the transmitting device: a constant signal amplitude at the output and a high linearity the frequency-time dependence.
The received echo signal from the boundary snow-ice layers at a distance r has a time delay t defined by the expression
t = 2r / Vrt, (11)
where Vrz - is the speed propagation the EW in the z-layer [2].
The frequency component the beat fb is directly related to the delay the echo signal (11), which is the difference between the expressions (9) and (10)
fb = fx -frx = aT = 2rB / Vrz Tm, (12)
where B - the spectral width (deviation) LFM signal.
The depth the snow-ice layers is determined by the resulting distance difference from formula (12), which pass the probing signals according to the formula
r = fb Vrz Tm / 2B
under normal probing to the underlying surface.
The total time delay, without taking into account the multi-layer snow-ice cover, Td to the boundary the ''ice cover - water'' Tiw medium is generally defined
2h 2hs 2h
c K K
where Tas - time delay to the boundary the ''air - snow cover'' environment; Ts, t - time delay in snow and ice cover, respectively; ha - distance from AHT to the ''air - snow'' medium boundary; hs - snow depth; hi - the depth the ice cover; Vs, Vi - the speed propagation electromagnetic effects in the snow and ice, respectively.
To improve the accuracy measuring the depth snow and thickness ice cover, when restoring the structure the underlying surface, it is necessary to take into account the propagation velocity EW [2] to each layer the medium (11).
The signal delay to the ''snow cover - ice cover'' medium boundary is determined Tis = Tas + ts.
The resolution the depth the FM-CW radar, when using an ultra-wide-band LFM signal that overlaps the C.S (3.75.15 cm) bands, with a frequency from 2 GHz to 8 GHz (B = 6 GHz) is about 4 cm [1, 2, 13], which is a good indicator when determining the possibility performing a safe landing a helicopter-type aircraft on an unprepared area with snow and ice cover.
Thus, the probability identifying the components the structure the snow-ice cover increases, thereby increasing the level safety landing a helicopter - type aircraft on an unprepared site with snow-ice cover, transporting (delivering) cargo and objects along a reservoir with snow-ice cover, by increasing the depth resolution, which when using an ultra-wide-band LFM signal with a frequency from 2 GHz to 8 GHz (B = 6 GHz) is about 4 cm. In this case, the error in determining the permittivity the layers (Fig. 2) at Brewster angles 6 BM+1 at the specified calculated values eri+1 and according to the formula eri+1 = (tg6Bi,i+1)2 in other words the methodological error is no more than 3 %. The value the error in estimating the measured values the relative permittivity eri+1exp the k layers when the values the noise level RMS increase a from 3.8 to 4.8 in increments of 0.1 for 100 implementations each with a confidence probability of 0.95 does not exceed 10 %, which indicates the validity using this method.
Method for remote identification the state snow-ice cover by the ratio of Fresnel reflection coefficients
Identification snow-ice cover layers. In [14] based on the radiation controlled area the sea surface microwave radio waves on an inclined polarization in the interval of 6 = 25°.. .75°, registering the scattered back signal, defined Bragg mechanism scattering, both on the vertical and horizontal polarizations, calculate the polarization ratio normalized backscattering cross sections measured at vertical avv0 and horizontal ahh0 polarizations and calculated relative permittivity the medium below the boundary the atmosphere-ocean.
The total coefficient reflection the Fresnel from the snow and ice cover without taking into account multiple reflections between the layer boundaries during oblique sounding by flat microwave radio waves with vertical and horizontal polarization is determined by the recurrent formula [4, 15] (2).
For horizontal polarization (hh - the first index the probing, second received radio signal) [15]:
cos °< - J(s<-+1 -s' (sm )2)
= 7-/ ■ 2 V (13)
I ^ cos ei + ^sM - s, (sm 9, f J 1
The dependence the reflection coefficient the Fresnel with vertical polarization calculated by the formula (1) on the complex relative permittivity the medium er and the sensing angle within the limits 6 = 0°.90° has a pronounced dip in values (Fig. 1a), called the Brewster angle, compared with the dependence the reflection coefficient the Fresnel with horizontal polarization Rhh calculated by the formula (13) (Fig. 1b).
When restoring the structure, the snow-ice cover according to this method, the EW [2, 11] propagation velocity is taken into account in the same way as the above method for estimating the state of the snow-ice cover by the Brewster angle.
The power the reflected signal from the snow and ice cover is determined by the formula [13]:
ptxx2g2 ch
jx-| | , (14)
r (4;:) (2/z)
where Ptx - power the transmitting device; X - wavelength; G - antenna gain; h - height in HTA.
On both polarizations, the radio signal (14) at a particular time has the same parameters Ptx, X, G, h and since this dependence is leveled when calculating the ratio the Fresnel reflection coefficients vertically and horizontally polarized signals, the ratio will be determined by three parameters: srm, 6 and fx. Expression (14) does not take into account the roughness layers snow-ice cover, since the radio wave is not mirrored, but is scattered at different angles, but when calculating the ratio, the Fresnel reflection coefficients vertically and horizontally polarized signals, this dependence is also leveled.
According to the proposed method, probing the radio signal simultaneously with vertical and horizontal polarization the controlled area the snow-ice cover within the limits 6 = 25°.. .45°, allows us to determine the dependence the reflection coefficients the Fresnel with vertical lV']'+] and horizontal
polarization, (i,i + 1)-boundary layers on angle incidence radio wave.
Relative permittivity snow-ice layers is determined from the ratio of Fresnel reflection coefficients by power Pim (Fig. 5) signals with vertical (Fig. la) obtained by formulas (1), (13) and |/^h|2 horizontal (Fig. 1b) obtained by formulas (1), (13) polarizations:
I Dhhl2
(15)
I ™ I
where |/^h|2, |i?™|2 - the reflection coefficients of the Fresnel on power (reflectivity) measured on the horizontal (hh and vertical (vv) polarizations, respectively (the first index is the polarization the probing, the second - the received radio signal); m - the number peaks the echo signal (the boundaries the layers snow-ice cover with different relative permittivity).
0
/ /
-C =1,2-10,0008= - r2 J ds dry snow;
--er3=l.9-j0.0008^ e^ - dry snow; ' / /
----£ =2,5-¡0,0008= rAJ df dry fim; 4 / #
-£ =3,240,0008= £- r5 J di dry ice; / # /
rfi J SW SW / /
jT 1 ✓ / ✓ / * ** ^
5
25
30
35
40
0°
Fig. 5. Ratios of Fresnel reflection coefficients by power Prm signals with vertical obtained by formulas (2), (4) and horizontal obtained by formulas (3), (4) polarizations: 1 - er2 = 1,2 - j0,0008, 2 - sr3 = 1,9 - j0,0008, 3 -sr4 = 2,5 - j0,0008, 4 - sr6 = 3,2 - j0,0008, 5 - sr7 = 74 - j, depending on the sensing angle within 6 = 25°.. .45°
The maximum angle oblique sensing 6 = 25°.. .45° the controlled section snow-ice cover is determined by the lower Brewster angle 6 = 46° for the layer snow-ice cover with the lowest relative permittivity dry snow - erds = 1,07 - j0,0008.
Fig. 5 shows the ratio of Fresnel reflection coefficients by power Prm for layers with complex relative permittivity: 1 - er2 = 1,2 - j0,0008, 2 - er3 = 1,9 - j0,0008, 3 - er4 = 2,5 - j0,0008, 4 - er6 = 3,2 - j0,0008, 5 - er7 = 74 - j, depending on the sensing angle within 6 = 25°.. .45°.
Successive determination the permittivity each subsequent layer up to erm, where m = 1.n - is the number layers snow-ice cover, is carried out by the formula:
, , 4. s = 1 +--n-sm(u)
r m / ,-\2 V /
(l-V^W)
which corresponds to the graphs in Fig. 5 the dependency (15).
The interval snow-ice cover densities will correspond to the intervals the Fresnel reflection coefficient ratios by power Prm. For example, when 6 = 46° for dry snow pds = 100. 500 kg/m3 (eds = 1,162.1,984) -Prm = 5,6915.3,3266, dry firn pdf = 500.700 kg/m3 (edf = 1,984.2,51) - Prm = 3,3255.2,8311, dry ice pdi = 700.917 kg/m3 (edi = 2,51.3,179) - Prm = 2,8311.2,4753. With an increase in the moisture content the values the Fresnel reflection coefficient ratios by power Prm decrease, aiming for values for water. For sea water esw = 74 (sea water) salinity Ssw = 35 g/kg this Prm = 1,1923, a for pure (thawed) water with epw = 87 - Prm = 1,1760.
To identify the layers snow-ice cover, the identity the obtained values the dielectric permittivity the layers erm with the specified calculated (model) values of the dielectric permittivity the layers evrA is established by the condition erm = evrA: ''snow cover'', "firn", ''ice cover'' or "water" similar to the above method.
Thus, the methodological error in determining the dielectric permittivity layers (Fig. 5) based on the ratios Fresnel reflection coefficients Prm for the given values erm and according to formula (16) is no more than 1,5 %, which indicates the validity of using this method.
Method for remote determination the state snow-ice cover by backscattering polarization relations
Identification snow-ice cover layers. The normalized cross section the backscattering avv(hh)0 in the case propagation the resonant components the wave field over a flat surface can be represented by the formula [14]:
,
where vv(hh) - vertical and horizontal polarization, the first index corresponds to the polarization the radiated signal, the second taken; krk - wave number radio waves; 5 (K]irk j - the range snow, ice or water surface, according to the relevant wave vector KBik the resonant components.
On both polarizations, a signal with the same wave net krk is proportional to the roughness the snow, ice, and water surface, respectively, and since this dependence is leveled when calculating the polarization ratio vertically and horizontally polarized signals, the polarization ratio will be determined by three parameters: erk, 6, and f. Knowing the polarization relations of the signals Prk, the sensing angle 6, and the frequency f iteratively determine the permittivity each subsequent layer erk+1.
,
(16)
The relative permittivity of snow-ice layers is determined from the polarization relations
backscattering prk (Fig. 6) signals with vertical crr7° (Fig. 7a) and horizontal erf0 (Fig. 7b) polarizations:
_wO
p _ ^rk _\ rg v ' ™ (17)
rrk - hho ~ I .. - . 1?. '
where crj0, erf0 - normalized cross section backscattering measured on the vertical (vv) and horizontal (hh) polarizations, respectively (the first index is the polarization the probing signal, the second - the received signal); k - the number peaks of the echo signal (the boundary between layers of snow-ice cover with different permittivity corresponding to the number of layers of snow-ice cover); |£Ik)f, |G^ {0, £-rA)| - geometric coefficients that depend on the type polarization, the angle incidence radio waves at the boundary layers, and the electrophysical (dielectric permittivity) parameters layers snow-ice cover a reservoir [16],
G- (0, ert) = cos4 0^ ^L^V-- - -J, (18)
.(«.-1)1 £lk (l + sin 0^ — sin
£rk cos0 + -sin2 <9J
G?(0, £[k) = cos40--^^ _ n2; (I9)
sxk cos 0 + yje!k - sin2 6
6 = 25° . .75° - the angle sensing the controlled area the snow-ice cover the reservoir determined by the resonant (Bragg) scattering mechanism.
From the formulas (17)-(19), the complex relative permittivity snow-ice layers is determined:
(20)
1 + sin #
The interval snow-ice cover densities will correspond to the intervals polarizing relations Prk backscattering signals. So for example, when 6 = 65° for dry snow pds = 100...500 kg/m3
Fig. 6. Polarization relations Prk backscattering signals with vertical avv°, obtained by formulas (17), (18) and horizontal ahh0, obtained by formulas (17), (19) polarizations: 1 - Pr2 = 1,8613, 2 - Pr3 = 9,0212, 3 - Pr4 = 12,6012, 4 - Pr5 = 23,2799, depending on the sensing angle within 0 = 25° . .75°
0.9 I 0.9
a) b)
Fig. 7. Dependence the normalized cross-section backscattering: a - with vertical polarization avv0; b - with horizontal polarization ahh0 on the complex relative permittivity the medium er and the sensing angle within 6 = 0°...90°
(£& = 1,162... 1,984) - Prk = 1,6772...7,7967, for dry firn pdf = 500.700 kg/m3 (e'df = 1,984.2,51) -Prk = 7,7967.14,0648, for dry ice pdm = 700.917 kg/m3 (e'di = 2,51.3,179) - Prk = 14,0648.24,6891. As the moisture content increases значения, the values the polarization ratios Prk the backscattering signals increase significantly, tending to the values for water. For sea water e'sw = 74 (sea water) salinity is Ssw = 35 g/kg this Prk = 17949, а for pure (meltwater) water with = 87 - PA = 24854.
The obtained values of the dielectric permittivity the layers erk are compared with the calculated (model) values the dielectric permittivity the layers evrA. By establishing the identity which condition erk = evrA determines the state the snow-ice cover the reservoir: ''snow cover'', ''firn'', ''ice cover'' or ''water''.
Methodological error in determining the permittivity of layers (Fig. 6) based on the polarization relations Plk the normalized cross-sections backscattering <x™°- <y,f0 at the specified values slk and according to the formula (20) is no more 1 %, which indicates the legality using this method.
Conclusion
Thus, the probability of correct identification of the constituent elements of the snow and ice cover structure increases, thereby increasing the level of safety of landing a helicopter-type aircraft on an unprepared site with snow and ice cover, by increasing the depth resolution, which when using ultrawideband LFM signal with a frequency from 2 GHz to 8 GHz (B = 6 GHz) is about 4 cm.
The error in determining the dielectric permittivity layers (Fig. 2) at Brewster angles 6Bi,M for the specified calculated values eri+1 and according to formula (5), i. e. the methodological error is no more than 3 %. The value of the error in estimating the measured values of the relative permittivity of eri+1exp k layers with increasing values of the RMS noise level a from 3.8 to 4.8 in increments of 0.1 for 100 implementations of each with a confidence probability of 0.95 does not exceed 10 %.
Thus, the methodological error in determining the dielectric permittivity layers (Fig. 5) based on the ratios Fresnel reflection coefficients Prm for the given values erm and according to formula (16) is no more than 1,5 %.
Methodological error in determining the permittivity of layers (Fig. 6) based on the polarization relations Plk the normalized cross-sections backscattering er"°, £7rf0 at the specified values slk and according to the formula (20) is no more 1 %.
The results obtained allow us to draw a conclusion about the applicability the proposed methods for assessing the state snow-ice cover used in determining the possibility performing a safe landing a helicopter-type aircraft on an unprepared site with snow-ice cover based on the identification the obtained characteristics snow-ice cover layers based on the results radar sensing with calculated data.
References
[1] Машков В. Г., Малышев В. А. Способ выбора площадки для посадки воздушного судна вертолетного типа. Пат. 2707275 Российская Федерация, МПК G01S13/94, заявитель и патентообладатель ВУНЦ ВВС «ВВА» (г. Воронеж), 2019100117, заявл. 09.01.2019, опубл. 26.11.2019 [Mashkov V. G., Malyshev V. A. Method selecting a landing site for a helicopter - type aircraft. Pat. 2707275 Russian Federation, IPC G01S13/94, applicant and patent holder VUNTS VVS «VVA» (Voronezh), 2019100117, declared 09.01.2019, publ. 26.11.2019 (in Russian)]
[2] Малышев В. А., Машков В. Г. Скорость распространения электромагнитной волны в снежно-ледяной подстилающей поверхности. Радиотехника. М.: Радиотехника, 2020, 3 (5), 40-54. DOI 10.18127/j00338486-202003(05)-05 [Malyshev V. A., Mashkov V. G. Electromagnetic wave propagation speed in the snow-ice underlying surface. Radio Engineering. M., Radio Engineering, 2020, 3 (5), 40-54. DOI 10.18127/j00338486-2003(05)-05 (in Russian)]
[3] Харвей А. Ф. Техника сверхвысоких частот. М.: Советское радио, 1965, 41-42 с. [Harvey A. F. Technique of ultrahigh frequencies. Moscow, Sovetskoe radio, 1965, 41-42 р. (in Russian)]
[4] Шостак А. С., Загоскин В. В., Лукьянов С. П., Карауш А. С. О возможности определения диэлектрической проницаемости верхних слоев подстилающих сред по измеренным коэффициентам отражения при наклонном зондировании плоскими волнами вертикальной и горизонтальной поляризации в СВЧ диапазоне. Журнал радиоэлектроники, 1999, 11 [Shostak A. S., Zagoskin V. V., Lukyanov S. P., Karaush A. S. On the possibility determining the dielectric permittivity the upper layers underlying medium by measured reflection coefficients during oblique sounding by plane waves vertical and horizontal polarization in the microwave range. Journal of Radioelectronics, 1999, 11 (in Russian)]
[5] Валеев Г. Г. Способ измерения относительной комплексной диэлектрической проницаемости материала с потерями в СВЧ диапазоне. Пат. 2613810 Российская Федерация, МПК G01R27/00 (2006.01), заявитель и патентообладатель Валеев Г. Г. 2015142390, заявл. 06.10.2015, опубл. 21.03.2017 [Valeev G. G. Method for measuring the relative complex permittivity of a material with losses in the microwave range. Pat. 2015142390 Russian Federation, 2613810, IPC G01R27/00 (2006.01), applicant and patent holder Valeev G. G. 2015142390, declared 06.10.2015, publ. 21.03.2017 (in Russian)]
[6] Финкельштейн М. И., Лазарев Э. И., Чижов А. Н. Радиолокационные аэроледомерные съемки рек, озер, водохранилищ. Под ред. М. И. Финкельштейн. Ленинград: Гидрометеоиздат, 1984, 112 с. [Finkelstein M. I., Lazarev E. I., Chizhov A. N. Radar aeroledomeric surveys of rivers, lakes, reservoirs. Ed. by M. I. Finkelstein. Leningrad, Hydrometeoizdat, 1984, 112 р. (in Russian)]
[7] Канарейкин Д. Б., Павлов Н. Ф., Потехин В. А. Поляризация радиолокационных сигналов. Под ред. В. Е. Дулевича. М.: Сов. радио, 1966. 440 с. [Kanareykin D. B., Pavlov N. F., Potekhin V. A. Polarization of radar signals. Ed. by V. E. Dulevich. M., Sov. radio, 1966. 440 p. (in Russian)]
[8] Гринев А. Ю., Темченко В. С., Багно Д. В. Радары подповерхностного зондирования. Мониторинг и диагностика сред и объектов. Монография. М.: Радиотехника, 2013. 74-85 c. [Grinev A. Yu., Temchenko V. S., Bagno D. V. Radars of subsurface sounding. Monitoring and diagnostics among objects. Monograph. Moscow: Radio Engineering, 2013. 74-85 p. (in Russian)]
[9] Поздняк С. И., Мелитицкий В. А. Введение в статистическую теорию поляризации радиоволн. М.: Сов. радио, 1974. 479 с. [Pozdnyak S. I., Melitinsky V. A. Introduction to the statistical theory of radio wave polarization. M.: Sov. radio, 1974. 479 p. (in Russian)]
[10] Мачерет Ю. Я. Оценка содержания воды в ледниках по гиперболическим отражениям. Материалы гляциологических исследований. Институт географии РАН, 2000, 89, 3-10. [Macheret Yu. Ya. Assessment water content in glaciers by hyperbolic reflections. Materials glaciologicalresearch. Institute geography the Russian Academy Sciences, 2000, 89, 3-10 (in Russian)]
[11] Котляков В. М., Мачерет Ю. Я., Сосновский А. В., Глазовский А. Ф. Скорость распространения радиоволн в сухом и влажном снежном покрове. Лед и снег. М.: Институт географии РАН, 2017, 57(1), 45-56. [Kotlyakov V. M., Macheret Yu. Ya., Sosnovsky A. V., Glazovsky A. F. Speed radio wave propagation in dry and wet snow cover. Ice and snow. Moscow, Institute geography the Russian Academy of Sciences, 2017, 57(1), 45-56 (in Russian)]
[12] Фролов А. Д., Мачерет Ю. Я. Оценка содержания воды в субполярных и теплых ледниках по данным измерений скорости распространения радиоволн. МГИ, 1998, 84, 148-154 [Frolov A. D., Macheret Yu. Ya. Estimation of water content in Subpolar and warm glaciers based on radio wave propagation velocity measurements. MGI, 1998, 84, 148-154 (in Russian)]
[13] Sudarsan Krishnan B. E. Modeling and simulation analysis of an FMCW radar for measuring snow thickness. Electronics and communication engineering. University of Madras, 2000, 93.
[14] Запевалов А. С. Способ дистанционного определения относительной диэлектрической проницаемости среды под границей атмосфера-океан. Пат. 2623668 Российская Федерация, МПК G01N27/06 (2006.01), G01R27/26 (2006.01), заявитель и патентообладатель МГИ РАН. 2015156757, заявл. 28.12.2015, опубл. 28.06.2017 [Zapevalov A. S. Method for remote determination of the relative permittivity of the medium under the atmosphere-ocean boundary. Pat. 2623668 Russian Federation, IPC G01N27/06 (2006.01), G01R27/26 (2006.01), applicant and patent holder of ICAI RAS. 2015156757, application 28.12.2015, publ. 28.06.2017 (in Russian)]
[15] О возможности определения диэлектрической проницаемости верхних слоев подстилающих сред по измеренным коэффициентам отражения при наклонном зондировании плоскими волнами вертикальной и горизонтальной поляризации в СВЧ диапазоне [Электронный ресурс] - Режим доступа: http://jre.cplire.ru/jre/nov99/4/text.html - Заглавие с экрана. [On the possibility determining the dielectric permittivity the upper layers underlying medium from the measured reflection coefficients for inclined sounding by plane waves vertical and horizontal polarization in the microwave range [Electronic resource] - Access: http://jre.cplire.ru/jre/nov99/4/text.html]
[16] Valenzuela G. Theories for the interaction of electromagnetic and ocean waves. A. Review. Boundary Layer Meteorology, 1978, 13(1-4), 61-85.
