CC BY
36
GEODYNAMICS & TECTONOPHYSICS
PUBLISHED BY THE INSTITUTE OF THE EARTH'S CRUST SIBERIAN BRANCH OF RUSSIAN ACADEMY OF SCIENCES
ISSN 2078-502X
2016 VOLUME 7 ISSUE 4 PAGES 545-553
http://dx.doi.org/10.5800/GT-2016-7-4-0222
The determination of atmospheric water content from
meteorological and GPS data
O. F. Lukhneva1, M. G. Dembelov2, A. V. Lukhnev1
1 Institute of the Earth's Crust, Siberian Branch of RAS, Irkutsk, Russia
2 Institute of Physical Materials Science, Siberian Branch of RAS, Ulan-Ude, Russia
Abstract: The Global Positioning System (GPS) based on satellites and the networks of dual frequency receivers are actively used for geodetic and geophysical applications, as well as for studying the ionosphere and troposphere. The atmospheric water content is in the focus of research as a key parameter for determining of the accuracy of weather forecasting and hydrological monitoring. The precision of atmospheric water content calculations depends on the accuracy of determination of the delays of signals propagating from GPS satellites to ground-based GPS receivers when geodynamic measurements are conducted. This paper describes a technique that allows us to estimate the integrated water vapor (IWV) in the atmosphere from measurements of GPS satellite signal delays.
We consider remote sensing of the lower atmosphere by GPS measurements to detect the water vapor content in the conventional vertical column to the top level of the troposphere (up to 12 km above the Earth's surface). In studies of the propagation of signals from GPS satellites to ground receivers, the atmospheric water vapor is taken into account as a 'wet' component (ZWD) of the zenith tropospheric delay (ZTD). ZTD is the sum of ZHD (hydrostatic or 'dry' delay) and ZWD ('wet' delay). ZWD values can be converted with a very high confidence in integrated water vapor (IWV) values for each installed GPS receiver.
Key words: GPS; remote sensing; meteorological data; troposphere; tropospheric zenith delay; moisture content of troposphere; precipitable water
Recommended by K.G. Levi
For citation: Lukhneva O.F., Dembelov M.G., Lukhnev A.V. 2016. The determination of atmospheric water content from meteorological and GPS data. Geodynamics & Tectonophysics 7 (4), 545-553. doi:10.5800/ GT-2016-7-4-0222.
Определение атмосферного влагосодержания по
метеорологическим и GPS-данным
О. Ф. Лухнева1, М. Г. Дембелов2, А. В. Лухнев1
1 Институт земной коры СО РАН, Иркутск, Россия
2 Институт физического материаловедения СО РАН, Улан-Удэ, Россия
Аннотация: Система спутникового позиционирования GPS с использованием сетей двухчастотных приемников активно применяется не только для решения задач геодинамики, но и для исследования ионосферы и тропосферы. Особый интерес представляет оценка атмосферного влагосодержания, так как это один из ведущих параметров определения точности прогнозов погоды и гидрологического мониторинга. Точность оценки влагосодержания определяет точность оценки задержки GPS-сигнала при геодинамических измере-
ниях. В работе описывается методика, позволяющая оценивать значение интегрального влагосодержания атмосферы по измеряемым фазовым задержкам сигнала спутников GPS.
Рассматривается дистанционное зондирование нижней части атмосферы посредством GPS-измерений с целью определения содержания водяного пара в условном вертикальном столбе до уровня верхней части тропосферы (до 12 км над поверхностью Земли). Атмосферный водяной пар учитывается в процессе распространения сигналов от GPS-спутников до наземных приемников в виде «влажной» компоненты полной тропосферной задержки (ZWD). Полная тропосферная зенитная задержка (ZTD) является суммой «сухой», или гидростатической (ZHD), и «влажной» (ZWD) компонент. Данные по ZWD могут быть преобразованы с очень высокой достоверностью в данные по суммарному водяному пару (IWV) над каждым установленным GPS-приемником.
Ключевые слова: GPS; дистанционное зондирование; метеорологические данные; тропосфера; тропосферная зенитная задержка; влагосодержание тропосферы; осаждаемая вода
1. Introduction
Permanent GPS stations in Irkutsk, IRKT and IRKM, are the base stations in the Baikal geodynamic network. High-precision geodetic measurements were conducted on the IRKT - ULAZ profile (from Irkutsk to Ulan-Ude) from 1999 to 2014 [Sankov et al., 2014]. In addition to the geodynamic parameters of crustal movements and deformations, the total tropospheric zenith delay data was collected. Changes in the refractive index with height in the lower neutral atmosphere cause tropospheric refraction or deviation of the rays of GPS signals. The index of radio wave refraction, N depends on temperature, pressure and aqueous tension. The zenith tropospheric delay (ZTD) is one of the most significant corrections [Davis et al., 1985] in high-precision geodetic calculations from GPS measurements [Lukhnev et al., 2013; Dembelov et al., 2015]. Then, using long term observations, it is also possible to study its changes in time (e.g., Emardson et al., 1998; Nilsson, Elgered, 2008; Wang, Zhang, 2009), which constitutes the beginning of its use in climate research. Water vapor is the most variable of all parameters of the troposphere and thus plays an important role in atmospheric processes [Khutorova et al., 2012; Kunitsyn et al., 2015]. The distribution of water vapor are closely related to the distribution of clouds and rain [Eminov, Magerramov, 2012]. The 'wet' component (ZWD) is practically proportional to the moisture content of clouds overlying a GPS station and can be estimated from GPS data on ZTD. Thus, GPS measurements can be useful for continuous remote sensing of the atmosphere. The emerging denser GPS networks, established mainly for studies of geodynamics [Sankov et al., 2014] and the ionosphere, offer the possibility of using GPS data for efficient weather forecasting [Dembelov et al., 2015].
GPS antennas on monitoring sites work with the cut-off angle of 10° relative to the horizon. Therefore, a satellite monitoring area is a cone with the base radius
of 68 km, and the tropopause border altitude is 12 km. The tropospheric delay varies from about 6 to 8 ns in time (1.9-2.4 m), depending on weather conditions and location. With decreasing elevation angle a this parameter increases as the cosecant of the angle, so a delay at L1=1575.42 MHz with angle a of 20° can range from 30 to 36 cycles (5.71-6.85 m).
A real tropospheric delay is estimated as a zenith direction in the form of cosec a to consider ZTD. Surface weather data for Irkutsk are taken from www.rp5.ru. Meteorological vertical sounding data for Angarsk are provided by the Irkutsk Department of Hydrometeo-rology and Environmental Monitoring.
GPS site (IRKM) is located on a distance of 6 km from the Irkutsk meteorological station (WMO ID: 30710) and 43 km from the Angarsk meteorological station (WMO ID: 30715). The error of vertical meteorological data records in Angarsk for IRKM station does not exceed 5 %.
2. Tropospheric refraction index and zenith tropospheric delay
A tropospheric delay depends on refraction factor, n. Radio physics widely uses the refraction index from [Smith, Weintraub, 1953]:
N = (n-l)-106 =^p + ^e = NT + Ne, (1)
where Nt is the refraction index of dry air, which depends on changes in air temperature and pressure; Ne is the refraction index of water vapour; ki=77.6 is the first refraction constant, K/mbar (1 mbar=102 kg/m-s2); k2=3.73-105 is the second refraction constant, K2/mbar; T is absolute temperature, K; p is atmospheric pressure, mbar; e is aqueous tension, mbar. For average weather conditions in Irkutsk, T=288 K, p=960 mbar, e=19 mbar in summer, and T=255 K, p=975 mbar, e=1 mbar in winter.
N, Nt
Months, 2015
Fig. 1. Refraction indices N, Nt and Ne by months in 2015 (Irkutsk).
Рис. 1. Индексы рефракции N, Nt и Ne по месяцам 2015 г. для г. Иркутска.
Figure 1 shows the annual variations of parameters N, Nt and Ne in Irkutsk in 2015. In general, 'wet' component Ne makes a much smaller contribution to the determination of the refractive index, especially in the cold time of the year. The refractive index is predetermined by 'dry' component Nt. In the warm season, from the middle of April to the middle of October, the variations of the refractive index N correlate well with the changes in the refractive index of water vapor Ne (K=0.95).
The total tropospheric zenith delay (ZTD) is the sum of 'dry' (hydrostatic) and 'wet' components, i.e. ZHD and ZWD. By definition, the total delay of the signal path from the satellite to the GPS receiving antenna (ZTD) is equal to the difference of the real signal path in the atmosphere and the geometric distance:, ZTD = f . n(h)dh— f, dh, where h is
JAtmosphere v y J Vacuum
height variable.
Using equation (1), components of the zenith delay can be expressed by integrating over the vertical profile of the corresponding refractive index values, NT(h) and Ne(h):
ZHD = 10~6 Ç NT (h)dh, ZWD = 10"6/0°°Л?е (h)dh.
(2) (3)
In practice, the integration is limited for water vapor to the upper part of the troposphere, i.e. about 12 km, and for dry air may be continued in the stratosphere.
3. Vertical profiles of tropospheric refraction and formulas for ZTD
Given the atmosphere in the state of hydrostatic equilibrium, components ZHD and ZWD (depending on the height) are well modeled with reference to pressure near the ground data, air humidity and temperature. The vertical profile of the total refraction is well approximated by an exponential formula [Moshkov, Pozhidaev, 2014]:
N(h) = N0exp(-/3h),
(4)
where N0 is the refraction index at the Earth's surface; P, km-1 is the rate of decrease of the refraction index with changes in height h, km. Parameter P in formula (4) is estimated from N measurements at different heights by the least squares method. For the Angarsk station, parameter P is 0.125 in summer, 0.132 in winter, and 0.13 in spring and autumn. Figure 2 shows curves of altitude profiles obtained from weather sounding data and approximations using formula (4) for Angarsk on 03 January 2014 and 14 July 2014. Seasonal values of parameters N0 and P in formula (4) allow obtaining reasonably reliable data on refraction index N(h) through the year. For a specified observation point, values on parameter P can be tabulated by seasons of the year, which can allow with a sufficient accuracy to calculate height profiles of N(h) without using the weather radio-sounding data.
Formula (4) for the vertical profile is also represented as the sum of the 'dry' and 'wet' components:
h, m a)
Angarsk, 03 January 2014
35000
— Exponential model • Sounding data
30000
25000
1
Ï
5 é
\
h, m b)
Angarsk, 14 July 2014
20000
— Exponential model • Sounding data
I Fig. 2. Comparative height profiles of the refractive index based on sounding data and approximations over the Angarsk station in winter (a) and summer (b).
I Рис. 2. Сравнительные высотные профили индекса рефракции по данным зондирования и их аппроксимации над пунктом наблюдения Ангарск для зимних (a) и летних условий (b).
N(h) = NTOexp(-ßDh) + Ne0exp(-ßwh),
(5)
where Nt0 and Ne0 are the near-surface refraction indices for dry air and water vapour, respectively; Pd and Pw are rates of decrease of the refraction index for dry air and water vapour, respectively. Values of Pd and Pw depend on the location of the monitoring station on the Earth. For the latitudes of the Baikal region, Pd is about 5/hD, and pW is about 0.5 [Dembelov et al., 2016]. In this formula, hD = 40136 + 148.72 • t is the effective height of the 'dry' (hydrostatic) component, m; t is surface temperature, °C [Bevis et al., 1992].
Representing ZTD according to formulas (2) and (3), given the components of the refraction indices according to formula (5), and performing the integration, the total tropospheric delay is given as the following sum:
ZTD = ZHD + ZWD = 10~6NT-^+10~6Ne-^. (6)
pd ßw
Formula (6) relates the refraction indices for dry air and water vapor with corresponding ZHD and ZWD. Hydrostatic delay ZHD can be estimated from meteorological data, using the well-known model described in [Saastamoinen, 1972], which, in its turn, was derived by transforming the integral (2) into the form convenient for use:
ZHD =
0.002277p
f(q>,hs) .
(7)
In this formula, hs) = 1 — 0.00266 cos 2^ — —0.00028hs where 9 is the geographical latitude of the station's location (degree); hs is the receiver's height above the seal level (km). For all the monitoring stations in the Baikal region, value f (ty,hs) in formula (7) deviates from 1 by no more than 0.08 %. Thus, in formula (7), this value can be ignored. Actually, in formula (7), only the atmospheric pressure should be
taken into account, which greatly simplifies the calculations of ZWD from GPS measurements by subtracting ZHD from ZTD.
4. Integrated water vapor in the troposphere
In formula (3), parameter ZWD is given with account of Ne as follows:
ZWD = 10~6k2 C-^dh.
z J0 T2
(8)
Weighted 'mean temperature' Tm over the reception station is introduced as the following ratio [Bevis et al., 1992]:
Ti _ Jo т
rr
1 rn e
L —^dh
Jo r2
(9)
With account of ratio (9), formula (8) can be given as follows:
ZWD = 10~6 —
Tm J0 T
(10)
For normal atmospheric conditions, a ratio can be derived for the specific gas constant for water vapor Rw [Bevis et al., 1992]:
eV mT
— = =R
PwvT W'
(11)
where pwv is the density of water vapor, kg/m3. With account of ratio (11) and the equation of the integrated vapor, IWV = JQro pwv dh, formula (10) can be given as follows:
ZWD = 10~6 —RwIWV.
(12)
Integrated water vapour (IWV) in the atmosphere over the measurement station is estimated as the vertically integrated mass of water vapor per unit area (kg/m2). The height of the precipitable water column (PW, m) corresponding to IWV can be calculated as PW = IWV/p, where p is the density of liquid water, kg/m3.
Fig. 3 shows the summer and winter high-altitude temperature profiles over Angarsk according to the vertical sounding data. The profiles have different heights (up to 23 km in winter, and up to 32 km in summer). Measurements taken at different heights may have errors and mistakes due to weather conditions, such as strong wind. Besides, in real conditions, it is quite difficult to accurately determine water vapor pressure (e) at different heights. Nonetheless, based on
350001
30000
25000
20000
15000
10000
5000
h, m
Angarsk
о
О
о «•
0
о*
° 03 January 2014 • 14 July 2014
A
-80
-60
-40~ -20 0
Temperature, °С
20
40
IFig. 3. High-altitude temperature profiles over Angarsk in winter and summer.
I Рис. 3. Высотные температурные профили над г. Ангарском для зимних и летних условий.
the available data, the relationship between the variations of parameter Tm and surface temperature is established.
Based on the vertical sounding data up to 35 km, average temperature Tm is calculated from formula (9). Figure 4 shows the linear regression relationships between surface temperature T and average temperature Tm for Angarsk in 2014 and 2015. Thus, the linear regression for Angarsk is Tm=70.27+0.73T in 2014, and Tm=78.3+0.7T in 2015. It is shown that the 'average temperature' up to the height of the lower stratosphere above the measurement station can be estimated from the surface temperature of air with the use of the regression equation.
Having processed the primary GPS data series [Lukhnev et al., 2013] by GAMIT software package precision, in addition to high-precision geodetic data, we obtained the atmospheric components in the form of the total zenith tropospheric delay (ZTD). Figure 5, a, shows the curve of annual ZTD variations for IRKM station. According to formula (7) from [Saastamoinen, 1972], taking into account the surface atmospheric pressure, the annual ZHD variations are estimated. From the difference of ZHD and ZTD values obtained by processing the GPS measurements, the annual ZWD variations are determined. Using formula (10), given the current temperature dependence of parameter Tm
310 300 290 280 270
^ 260
250 240
230 — 230
310 300 290 280
^ 270 ^ 260 250 240
230
230
Angarsk, 2014
a)
^"YCo о о
240 250 260 270
Angarsk, 2015
280 290
300
Ь)
•__•• •
240
250
260
270
Г,К
280
290
300
310
310
Fig. 4. Parameters T and Tm, estimated from the weather sounding data for Angarsk. The data for 2014 (a) and 2015 (b)
Рис. 4. Соотношение п 2014 г. (a) и 2015 г. (b).
Рис. 4. Соотношение параметров Т и Тга, выявленных по данным метеорологических зондирований в г. Ангарск в
for Angarsk in 2015, the annual IWV variations over IRKM stations are estimated (Fig. 5, b). It should be noted that similar to the case for the refraction index and its components for water vapor, the cross-correlation coefficient of ZTD and IWV functions is quite high (0.92) in the warm period, from April to October.
5. Moisture content and precipitation
Using equation (12), it is possible to estimate an integrated water vapor (IWV) value over a permanent
GPS station. During warm months, the measured ZWD value is practically proportional to the current humidity, and meteorologists can thus use the GPS network for remote sensing of the atmosphere.
In this paper, we consider variations in the humidity content of the troposphere in the warm season of 2015 (from April to October) over IRKM station (Irkutsk). According to the annual ZTD data obtained by processing the GPS data series (see Fig. 5, a), we estimate the 'wet' component of ZWD and the precipitable water (PW) levels that correspond to IWV. Figure 6 shows the humidity content of the troposphere on the scale of PW and the levels of precipitation accumulated during the
Months, 2015
Fig. 5. Annual ZTD (a) and IWV (b) variations according to GPS measurements. The data for IRKM station (Irkutsk).
I Рис. 5. Графики годового хода полной тропосферной зенитной задержки (а) и суммарного водяного пара (b) по данным GPS-измерений для пункта IRKM (Иркутск).
50
IV V VI VII VIII IX X
Months, 2015
IFig. 6. Variations of the humidity content (PW) and the surface temperature in Irkutsk in the period from April to October 2015. Columns show the levels of accumulated precipitation during six hours.
Рис. 6. Вариации влагосодержания (PW) и приземной температуры в г. Иркутске за период апрель-октябрь 2015 г. Столбиками показаны уровни осадков, накопленных за 6 часов.
six hours of observation in Irkutsk. It should be noted that a high level of PW does not always correlate with precipitation due to the insufficient degree of vapour saturation at the given moment. An increase in the humidity content may be associated with precipitation, as a rule, in case of a sharp drop of air temperature and the presence of updrafts of air masses [Marchenko et al., 2012]. In Figure 6, precipitation correlated with the increase in the humidity content and a sharp decrease in surface temperature.
6. Conclusion
Based on the exponential model, the integrated refraction N(h) is calculated to the altitude of 35 km, and the value of the zenith tropospheric delay (ZTD) is predicted on the basis of the analysis and processing of the surface meteorological data. In this paper, we have described the technique that allows estimating the integrated water vapour (IWV) from ZTD that is obtained by processing the primary GPS data series
using GAMIT software package. The 'wet' component is the difference between ZTD and the 'dry' (hydrostatic) component calculated according to [Saastamoinen, 1972].
The linear regressions between surface temperature (T) and weighted mean temperature (Tm) are established for the Angarsk station. Using dependences of parameter Tm from surface temperatures, it is possible to more accurately estimate the moisture content of the troposphere over the observation point.
Global Positioning Systems (GPS/GLONASS) provide for more efficient hydrological monitoring of the atmosphere and better forecasting of atmospheric precipitation.
7. Acknowledgments
The study was financially supported by the Russian Fundation for Basic Research and the Ministry of Education and Research of the Buryatia Republic (Research Project No.15-45-04449-r_Siberia_a).
8. References
Bevis M., Businger S., Herring T., Rocken C., Anthes V., Ware R., 1992. GPS meteorology: remote sensing of atmospheric water vapor using the Global Positioning System. Journal of Geophysical Research: Atmospheres 97 (D14), 1578715801. http://dx.doi.org/10.1029/92JD01517.
Davis J.L., Herring T.A., Shapiro 1.1., Rogers A.E.E., Elgered G., 1985. Geodesy by radio interferometry: effects of atmospheric modeling errors on estimates of baseline length. Radio science 20 (6), 1593-1607. http://dx.doi.org/10. 1029/RS020i006p01593.
Dembelov M.G., Bashkuev Yu.B., Lukhnev A.V., Lukhneva O.F., San'kov V.A., 2016. The moisture content in the troposphere in the Baikal region from GPS measurements. Zhurnal Radioelektroniki (Journal of Radio Electronics) (3), 1-16 (in Russian) [Дембелов М.Г., Башкуев Ю.Б., Лухнев А.В., Лухнева О.Ф., Саньков В.А. Влагосодержание тропосферы в Байкальском регионе по данным GPS измерений // Журнал радиоэлектроники. 2016. № 3. С. 1-16]. Available from: http://jre.cplire.ru/jre/mar16/10/text.pdf.
DembelovM.G., Bashkuev Yu.B., LukhnevA.V., Lukhneva O.F., San'kov V.A., 2015. Diagnostics of atmospheric water vapor content according to GPS measurements. Atmospheric and Oceanic Optics 28 (4), 291-296. http://dx.doi.org/ 10.1134/S1024856015040053.
Emardson T.R., Elgered G., Johannson J.M., 1998. Three months of continuous monitoring of atmospheric water vapor with a network of Global Positioning System receivers. Journal of Geophysical Research: Atmospheres 103 (D2), 1807-1820. http://dx.doi.org/10.1029/97JD03015.
Eminov R.A., Magerramov E.I., 2012. To the issue of GPS measurements of the total amount of precipitated water. Zhurnal Radioelektroniki (Journal of Radio Electronics) (7) (in Russian) [Эминов Р.А., Магеррамов Э.И. К вопросу о GPS измерениях общего количестве осаждаемой воды // Журнал радиоэлектроники. 2012. № 7. C. 1-7]. Available from: http://jre.cplire.ru/iso/jul12/8/text.pdf.
Khutorova O.G., Kalinnikov V.V., Kurbangaliev T.R., 2012. Variations in the atmospheric integrated water vapor from phase measurements made with receivers of satellite navigation systems. Atmospheric and Oceanic Optics 25 (6), 429-433. http://dx.doi.org/10.1134/S1024856012060073.
Kunitsyn V.E., Nesterov I.A., Tereshin N.A., 2015. Analysis of atmospheric moisture content according to GPS receivers data. Zhurnal Radioelektroniki (Journal of Radio Electronics) (6), 1-12 (in Russian) [Куницын В.Е., Нестеров И.А., Терешин Н.А. Анализ влагосодержания атмосферы по данным приемников GPS // Журнал радиоэлектроники. 2015. № 6. C. 1-12]. Available from: http://jre.cplire.ru/jre/jun15/12/text.pdf.
Lukhnev A.V., San'kov V.A., Miroshnichenko A.I., Ashurkov S.V., Byzov L.M., San'kov A.V., Bashkuev Yu.B., Dembelov M.G., Calais E., 2013. GPS-measurements of recent crustal deformation in the junction zone of the rift segments in the central Baikal rift system. Russian Geology and Geophysics 54 (11), 1417-1426. http://dx.doi.org/10.1016/j.rgg. 2013.10.010.
Marchenko O.Yu., Mordvinov V.I., Antokhin P.N., 2012. The study of long-term variability and conditions for the formation of precipitation in the Selenga river basin. Atmospheric and Oceanic Optics 25 (12), 1084-1090 (in Russian)
[Марченко О.Ю., Мордвинов В.И., Антохин П.Н. Исследование долговременной изменчивости и условий формирования атмосферных осадков в бассейне реки Селенга // Оптика атмосферы и океана. 2012. Т. 25. № 12. С. 1084-1090].
Moshkov A.V., Pozhidaev V.N., 2014. Features of radio wave refraction in the near-polar regions. Journal of Communications Technology and Electronics 59 (11), 1107-1111. http://dx.doi.org/10.1134/S1064226914110163.
Nilsson T., Elgered G., 2008. Long-term trends in the atmospheric water vapor content estimated from ground-based GPS data. Journal of Geophysical Research: Atmospheres 113 (D19), D19101. http://dx.doi.org/10.1029/2008JD 010110.
Saastamoinen J., 1972. Atmospheric correction for the troposphere and stratosphere in radio ranging satellites. In: S.W. Henriksen, A. Mancini, B.H. Chovitz (Eds.), The use of artificial satellites for geodesy. AGU Geophysical Monograph Series, vol. 15, p. 247-251. http://dx.doi.org/10.1029/GM015p0247.
Sankov V.A., Lukhnev A.V., Miroshnitchenko A.I., Dobrynina A.A., Ashurkov S.V., Byzov L.M., Dembelov M.G., Calais E., Déverchère J., 2014. Contemporary horizontal movements and seismicity of the south Baikal Basin (Baikal rift system). Izvestiya, Physics of the Solid Earth 50 (6), 785-794. http://dx.doi.org/10.1134/S106935131406007X.
Smith E.K., Weintraub S., 1953. The constants in the equation for atmospheric refractive index at radio frequencies. Proceedings of the IRE 41 (8), 1035-1037. http://dx.doi.org/10.1109/JRPR0C.1953.274297.
Wang J., Zhang L., 2009. Climate applications of a global, 2-hourly atmospheric precipitable water dataset derived from IGS tropospheric products. Journal of Geodesy 83 (3), 209-217. http://dx.doi.org/10.1007/s00190-008-0238-5.
Lukhneva, Olga F., Candidate of Geology and Mineralogy, Researcher Institute of the Earth's Crust, Siberian Branch of RAS 128 Lermontov street, Irkutsk 664033, Russia Tel.: (3952)429534; e-mail: olgal@crust.irk.ru
Лухнева Ольга Федоровна, канд. геол.-мин. наук, н.с. Институт земной коры СО РАН 664033, Иркутск, ул. Лермонтова, 128, Россия Тел.: (3952)429534; e-mail: olgal@crust.irk.ru
Dembelov, Mikhail G., Candidate of Physics and Mathematics, Senior Researcher Institute of Physical Materials Science, Siberian Branch of RAS 6 Sakhyanova street, Ulan-Ude 670047, Russia Tel.: (3012)433210; e-mail: mdembelov@mail.ru
Дембелов Михаил Георгиевич, канд. физ.-мат. наук, с.н.с. Институт физического материаловедения СО РАН 670047, Улан-Удэ, ул. Сахьяновой, 6, Россия Тел. (3012) 433210; e-mail: mdembelov@mail.ru
Lukhnev, Andrei V., Candidate of Geology and Mineralogy, Senior Researcher
Institute of the Earth's Crust, Siberian Branch of RAS
128 Lermontov street, Irkutsk 664033, Russia
Tel.: +7(3952)429534; e-mail: loukhnev@crust.irk.ru
Лухнев Андрей Викторович, канд. геол.-мин. наук, с.н.с. Институт земной коры СО РАН 664033, Иркутск, ул. Лермонтова, 128, Россия Тел.: (3952)429534; e-mail: loukhnev@crust.irk.ru
