Применение параметризации эмпирическими ортогональными функциями в задаче восстановления термической структуры тропосферы по радиометрическим данным

В работе описан математический аппарат метода оптимальной оценки для восстановления профилей температуры тропосферы (0-10 км) по радиометрическим данным с параметризацией искомых профилей в виде разложения по эмпирическим ортогональным функциям, получаемым путём сингулярного анализа ковариационной матрицы данных радиозондовых измерений. Показано, что в рамках такого подхода можно использовать сетку высокого разрешения по высоте при сравнительно небольшой размерности обратной задачи. Это позволяет корректно анализировать статистику радиозондовых измерений и оперативно восстанавливать профили температуры без привлечения больших численных мощностей. Проведено тестовое восстановление температуры по большому ансамблю модельных зашумлённых спектров собственного излучения атмосферы на основе статистики радиозондовых измерений температуры над Нижним Новгородом. Установлено, что среднеквадратичная ошибка восстановления температуры по спектрам в диапазонах частот 50-55 ; 55-59 и 50-59 ГГц варьируется в диапазонах 0,5-2,7 К (на высотах 0-10 км), 0,3-1,0 К (на высотах 0-2 км) и 0,3-2,8 К (на высотах 0-10 км) соответственно. При этом оптимальное количество эмпирических ортогональных функций (размерность задачи) составляет 6-8 в случае спектров в диапазоне частот 50-55 ГГц, 8 в случае 55-59 ГГц и 10 в случае 50-59 ГГц.

На английском языке

Application of Empirical Orthogonal Functions Parameterization in the Problem of Retrieval of the Tropospheric Thermal Structure by Radiometric Data

Belikovich M.V.

Kulikov M.Yu.

Ryskin V.G.

Shvetsov A.A.

Krasilnikov A.A.

Skalyga N.K.

Serov E.A.

Feigin A.M.

We describe the mathematical tools of the optimal estimation method for retrieval of tropospheric temperature profiles (0-10 km) from radiometric data with the desired profiles parameterized in the form of expansion over empirical orthogonal functions obtained by a singular analysis of the covariance matrix of the radiosonde measurement data. It is shown that within the framework of such an approach one can use a high-resolution altitude grid for a relatively small dimension of the inverse problem. This permits one to properly analyze the statistics of radiosonde measurements and, at the same time, retrieve the temperature profiles without using a large computing power. We have conducted test retrievals of temperature using a large ensemble of simulated noisy spectra of atmospheric self-radiation based on the statistics of radiosonde temperature measurements above Nizhny Novgorod. It is found that the r.m.s. error of temperature retrieval from spectra in frequency ranges 50-55 GHz, 55-59 GHz, and 50-59 GHz varies within the ranges 0.5-2.7 K (at altitudes of 0-10 km), 0.3-1 K (at altitudes of 0-2 km), and 0.3-2.8 K (at altitudes of 0-10 km), respectively. The optimal number of empirical orthogonal functions (dimension of the problem) is 6-8 in the case of the 50-55 GHz spectra, 8 in the case of the 55-59 GHz spectra, and 10 in the case of the 50-59 GHz spectra.