УДК 551.582 (571.56) М. С. Васильев, С. В. Николашкин, Р. Р. Каримов
СРАВНЕНИЕ ПРИЗЕМНОЙ ТЕМПЕРАТУРЫ ВОЗДУХА В ЯКУТИИ ПО ДАННЫМ РЕАНАЛИЗА И НАЗЕМНЫХ НАБЛЮДЕНИЙ
Приводятся результаты сравнительного анализа приземной температуры воздуха на территории Якутии по данным реанализа и наземных наблюдений. Использованы данные метеорологических станций ФГБУ «ВНИИГМИ-МЦД» за период с января 1979 г. по декабрь 2012 г., находящихся в северной, центральной, южной, западной и восточной части Якутии, а также данные климатических архивов ERA-Interim и NCEP/NCAR, которые основаны на регулярных метеорологических наблюдениях, аэрологической и спутниковой информации. Проведено сопоставление среднегодовых и среднесезонных (зима, лето) значений приземной температуры воздуха и рассчитаны коэффициенты корреляции, показано преимущество одного типа реанализа над другим. Для визуализации приземной температуры воздуха на территории Якутии по данным реанализа с широтно-долготным разрешением 1 x 1 градус использовано программное обеспечение «CMPAS» и геоинформационный продукт ArcView 3.2. По данным реанализа построены различные пространственно-временные карты приземной температуры воздуха на высоте 2 м от поверхности земли. По результатам сравнительного анализа сделаны выводы о том, что данные реанализа ERA-Interim и NCEP/NCAR в целом достаточно хорошо отображают изменчивость температуры воздуха в Якутии относительно наземных наблюдений, могут быть дополнены и применены в различных областях исследований.
Ключевые слова: приземная температура воздуха, среднегодовые и среднесезонные значения, корреляционный анализ, реанализ, European Interim ReAnalysis, National Centers Environmental Prediction / National Center for Atmospheric Research, наземные наблюдения, Cloud Map Processing Script, ArcView 3.2, Якутия.
M. S. Vasiliev, S. V. Nikolashkin, R. R. Karimov
Comparison of Surface Air Temperatures Variations in Yakutia Derived from Reanalysis and Ground Measurements
The results of a comparative analysis of surface air temperature variations in Yakutia derived from reanalysis data and ground-based observations are presented. The meteorological data (1979-2012) from the RIHMI-WDC stations located throughout Yakutia and ERA-Interim along with NCEP/NCAR reanalysis data were used. Reanalysis data are based on regular meteorological observations, including upper-air satellite measurements. A comparison of annual and seasonal (winter, summer) variations of surface air temperature is made. Analysis of correlation coefficients showed that one type of reanalysis has some advantages over another. "CMPAS" and ArcView 3.2 software were used to visualize surface air temperature data derived from reanalysis with spatial resolution of 1 x 1 degree. Maps of spatial distribution of surface air temperature at different heights (from 2 m to the ground surface) were plotted using reanalysis data. Comparative analysis shows that ERA-Interim and NCEP/NCAR reanalysis data
ВАСИЛЬЕВ Михаил Семенович - м. н. с. лаборатории оптики атмосферы Института космофизических исследований и аэрономии им. Ю. Г. Шафера СО РАН.
E-mail: [email protected]
VASILIEV Mikhail Semyonovich - Junior Scientific Researcher of the Laboratory of Atmospheric Optics, the Institute of Cosmophysical Research and Aeronomy named after Yu. G. Shafer, the Siberian Branch of Russian Academy of Sciences.
E-mail: [email protected]
НИКОЛАШКИН Семен Викторович - к. ф.-м. н., зав. лабораторией оптики атмосферы Института космофизических исследований и аэрономии им. Ю. Г. Шафера СО РАН.
E-mail: [email protected]
NIKOLASHKIN Semyon Viktorovich - Candidate of Physical-Mathematical Sciences, Head of the Laboratory of Atmospheric
Optics, the Institute of Cosmophysical Research and Aeronomy named after Yu. G. Shafer, the Siberian Branch of Russian Academy of Sciences.
E-mail: [email protected]
КАРИМОВ Рустам Рамильевич - к. ф.-м. н., с. н. с., лаборатории радиоизлучений ионосферы и магнитосферы Института космофизических исследований и аэрономии им. Ю. Г. Шафера СО РАН.
E-mail: [email protected]
KARIMOV Rustam Ramilyevich - Candidate of Physical-Mathematical Sciences, Senior Scientific Researcher of the Laboratory of Radio Emissions of Ionosphere, the Institute of Cosmophysical Research and Aeronomy named after Yu. G. Shafer, the Siberian Branch of Russian Academy of Sciences.
E-mail: [email protected]
generally well reflect the variability of temperatures in Yakutia in comparison with ground-based observations and can be used in various research fields.
Key words: air temperatures, annual and seasonal mean values, correlation analysis, reanalysis, European Interim ReAnalysis, National Centers Environmental Prediction / National Center for Atmospheric Research, ground-based observations, Cloud Map Processing Script, ArcView 3.2, Yakutia.
Введение
Для оценки происходящих на обширной территории Якутии климатических изменений и связанных с ними изменений окружающей среды необходима подробная климатическая информация. Однако в связи с недостаточным количеством сети метеостанций (~100) отличаются пропуски в наблюдениях и существует сложность получения имеющихся данных. В настоящее время назрела необходимость регионализации существующих математических моделей глобального климата, увеличение сети метеостанций и организация комплексных инструментальных наблюдений большого числа одновременно измеряемых метеорологических параметров.
В последнее десятилетие все большее число исследований основывается на доступных климатических архивах, в которых первичные наблюдения подвергнуты обработке и приведены к узлам регулярной сетки. При построении таких архивов используются различные методы пространственного обобщения, от пространственной интерполяции до сложного модельного реанализа [1-6]. Данные реанализа представляют собой результаты расчетов глобальных моделей циркуляции атмосферы при реальных начальных и граничных условиях, основанных на регулярных метеорологических наблюдениях, спутниковой и аэрологической информации. Вопрос, насколько хорошо реанализ описывает изменчивость климата, в значительной степени остается открытым, что приводит к необходимости сопоставления таких рядов данных с данными, полученными экспериментальным путем [7-8].
Безусловно, одним из основных метеорологических параметров является приземная температура воздуха (ПТВ, оС), который, по мнению многих научных сообществ, характеризует крупномасштабные аномалии климатической системы и энергетического баланса Земли [9].
Можно отметить ряд работ М. К. Гавриловой, посвященных исследованиям пространственно-временных вариаций ПТВ на территории Якутии с использованием данных метеорологических станций [10-13], а ПТВ на территории России и за рубежом изучались, например, в работах [14-19].
Целью настоящего исследования является проведение сравнительного анализа ПТВ в Якутии по данным реанализа и наземных наблюдений за период с января 1979 г. по декабрь 2012 г.
Используемые материалы
В работе были использованы среднемесячные значения ПТВ метеорологических станций на территории Якутии (ФГБУ «ВНИИГМИ-МЦД»), а также реанализов ERA-Interim и NCEP/NCAR.
Реанализ ERA-Interim (European Interim ReAnalysis) использует модель ECMWF (European Centre for Medium-Range Weather Forecasts) CY29R1 с 6-часовым шагом по времени, пространственным разрешением 50 км для 60 слоев в атмосфере, а также четырехмерным усвоением данных [2]. Анализировались данные ПТВ на высоте 2 м с наименьшим широтно-долготным разрешением 0,125 x 0,125 градусов.
Реанализ NCEP/NCAR (National Centers Environmental Prediction / National Center for Atmospheric Research) - это объединение двух национальных центров США: экологического прогнозирования и атмосферных исследований. NCEP/NCAR использует прогностическую модель NCEP T62L28 с пространственным разрешением около 210 км (на экваторе), с 28 слоями в атмосфере и 6-часовым разрешением по времени [4]. Данные ПТВ на высоте 2 м представлены на гауссовой сетке с шагом ~1,9о по широте и 1,875° по долготе.
Данные реанализа ERA-Interim (с января 1979 г. по настоящее время - www.ecmwf.int) и NCEP/NCAR (с января 1948 г. по настоящее время - www.esrl.noaa.gov) находятся в свободном доступе.
Методика пространственного сопоставления данных
Пространственное сопоставление данных реанализа было проведено с данными метеорологических станций в разных частях Якутии: север - ст. Тикси (7P,35'N, 128^55^); центр - ст. Якутск (62о,05^, 129^45^); юг - ст. Алдан (58о,37^, 125^22^); запад -ст. Шелагонцы (66оД5^, 114оД7,E); восток - ст. Усть-Мома (66о,27^, 143^14^).
Вокруг метеостанций были выбраны «ячейки» с соответствующим широтно-долготным разрешением данных реанализа. Пример метода пространственного сопоставления данных реанализа и ст. Якутск показан на рис. 1.
Сравнение среднегодовых значений ПТВ
В табл. 1 представлены среднегодовые значения ПТВ за период 1979-2012 гг., а также их средне-квадратические отклонения в межгодовой изменчивости.
Среднегодовые значения ПТВ в Якутии за исследуемый период по данным наземных наблюдений находятся в диапазоне от -5.5 до -13.2, а по данным
Рис. 1. Метод пространственного сопоставления данных реанализа и наземных наблюдений на примере вокруг ст. Якутск
реанализа ERA-Interim и NCEP/NCAR - от -6.8 до -14.0 и от -7.6 до -13.1 градуса соответственно. Среднегодовые значения ПТВ по данным ERA-Interim близки к данным наземных наблюдений в центральной, южной и восточной частей Якутии относительно данных NCEP/NCAR. В северной и западной частях Якутии наблюдается схожесть данных NCEP/NCAR с данными метеорологических станций.
На рис. 2 для примера представлен график годового хода среднемесячных значений ПТВ по данным реанализа и данным ст. Якутск за период 1979-2012 гг.
Видно (рис. 2), что вариации ПТВ по данным ERA-Interim (штрихпунктирная линия с квадратным маркером) практически легли на кривую вариаций ПТВ по данным ст. Якутск (сплошная линия). Таким образом, исходя из табл. 1, можно считать, что данные ERA-Interim на большей части Якутии в целом наиболее отчетливее
отображают среднегодовую изменчивость ПТВ в Якутии за период 1979-2012 гг., однако, если рассматривать различные периоды усреднения, картина преобладания одного реанализа над другим может измениться.
Корреляционный анализ среднегодовых значений ПТВ
В работе нами были рассчитаны коэффициенты корреляции Пирсона среднегодовых значений ПТВ между данными наземных наблюдений и реанализа за исследуемый период с уровнем значимости 99 % при r = 0.45 (табл. 2).
А крит v '
Из табл. 2 видно, что коэффициенты корреляции между данными реанализа и наземных наблюдений находятся в достаточно хорошей согласованности с преобладанием данных ERA-Interim.
Необходимо отметить, что при использовании и сопоставлении данных реанализа (среднегодовых,
Таблица 1
Среднегодовые значения ПТВ (в скобках указано стандартное квадратичное отклонение)
Данные ст. Тикси ст. Якутск ст. Алдан ст. Шелагонцы ст. Усть-Мома
Наземные наблюдения -12.8 (1.5) -8.9 (1.2) -5.5 (1.0) -12.6 (1.3) -13.2 (1.2)
ERA-Interim -14.0 (1.6) -8.9 (1.1) -6.8 (1.0) -10.4 (1.3) -12.9 (1.0)
NCEP/NCAR -13.1 (1.2) -9.7 (0.8) -7.6 (1.0) -10.9 (1.1) -12.8 (1.0)
Центральная часть Якутии
Месяцы
-ст. Якутск о ЕЯА-Ыепш —а— ЩРЩК
Рис. 2. Годовой ход среднемесячных значений ПТВ по данным реанализа и ст. Якутск за период 1979-2012 гг.
среднесезонных и среднесуточных вариаций) на территории Якутии с данными наземных наблюдений есть возможность расчета и ввода дополнительных поправочных коэффициентов, исходя из интересующих пользователей определенных промежутков времени, что приводит к получению более близкой картины, отображающей изменчивость ПТВ относительно данных метеорологических станций. Например, временной ход ПТВ ст. Тикси и ЕКА-Шепт за период 1979-2012 гг. практически одинаков (Я = 0.916), но амплитуды вариаций несколько отличаются. Зависимость ПТВ для ст. Тикси и ЕЯЛ-Мент имеет следующий вид: Т Т = 0.916 х ТЕ_... . . Используя
ст. Тикси ЕкА-тгепт ^
аналогичные уравнения регрессии, можно скорректировать данные реанализа относительно данных метеостанций.
Сравнение среднесезонных значений ПТВ
В табл. 3 представлены среднесезонные
значения ПТВ в зимний (декабрь-февраль) и летний (июнь-август) периоды 1979-2012 гг.
Среднесезонные значения ПТВ (1979-2012) в зимний период на территории Якутии по данным наземных наблюдений находятся в диапазоне от -25 оС до -41 оС, а по данным реанализа ЕИА-Щепт и NCEP/NCAR от -28 до -37 и от -27 до -33 градусов по Цельсию соответственно. В летний период значения по данным наземных наблюдений находятся в диапазоне от 6 до 17 оС, а по данным реанализа ЕИА-Щепт и NCEP/NCAR от 8 до 17 и от 8 до 15 градусов Цельсия соответственно.
Исходя из табл. 3, в зимний и летний период преимущественно наблюдается схожесть данных ЕИА-Щепт и наземных наблюдений в западной, центральной и восточной частях Якутии. Данные NCEP/NCAR в зимний период наиболее отчетливо отображают среднесезонную изменчивость ПТВ в
Таблица 2
Корреляционный анализ среднегодовых значений ПТВ
Корреляция (Я) ст. Тикси ст. Якутск ст. Алдан ст. Шелагонцы ст. Усть- Мома
Наземные наблюдения & ЕЯЛ-Шеит 0.916 0.914 0.899 0.924 0.848
Наземные наблюдения & ЫСЕР/ЫСАК 0.878 0.669 0.872 0.879 0.810
Таблица 3 Среднесезонные значения ПТВ
Данные ст. Тикси ст. Якутск ст. Алдан ст. Шелагонцы ст. Усть-Мома
Зима/Лето
Наземные наблюдения -29/6 -37/17 -25/15 -38/12 -41/14
ЕЯЛ-Шеит -34/8 -37/17 -28/14 -33/12 -36/10
ЫСЕР/ЫСАЯ -31/8 -33/15 -27/14 -31/12 -32/10
Рис. 3. Карта-схема преобладания одного реанализа над другим в зимний и летний периоды 1979-2012 гг. Знак «~» означает возможность одновременного (равнозначного) использования данных реанализа
северной и южной частях Якутии относительно данных ERA-Interim. В летний период на западе Якутии данные ERA-Interim и NCEP/NCAR одинаково схожи с данными наземных наблюдений и т. д.
На рис. 3 показано преобладание (сходство данных с наземными наблюдениями) одного реанализа над другим в зимний и летний периоды 1979-2012 гг. Визуализация данных реанализа ПТВ в Якутии В данной работе для визуализации данных реанализа было принято решение использовать программное обеспечение (ПО) «Cloud Map Processing
Script (CMPAS)» [20], так как оно изначально создавалось как универсальное (с возможностью модернизации в зависимости от тех или иных требований пользователя). «CMPAS» использовано для сортировки данных реанализа ERA-Interim с широтно-долготным разрешением 1 x 1 градус в «удобочитаемый» вариант, доступный для визуализации данных геоинформационным продуктом ArcView 3.2.
На рис. 4 представлен вид фрагмента файла данных ERA-Interim, полученных с помощью обработки ПО «CMPAS».
result.Ш Ь.юнног
ИДИ
Файл Правка Формат Вид Справка
"Ion "tat id order val ue
66 ы 1 1 13
67 62 1 2 12
67 61 1 3 13
66 61 1 4 11
6 С 62 1 5 14
67 62 2 1 10
68 62 2 2 12
68 61 2 3 11
67 61 2 4 13
67 62 2 5 13
68 62 3 1 12
Рис. 4. Фрагмент файла данных ERA-Interim, полученных на выходе ПО «CMPAS». В верхней строке: lon - географическая долгота; lat - географическая широта; id - идентификационный номер; value - значение ПТВ (оС)
На рис. 5 представлен возможный пример отображения пространственного распределения ПТВ в Якутии по данным ЕЯА-Шепт в среде геоинформационного продукта ArcView 3.2, у которого существует возможность наложения элементов топоосновы и гидрографии.
Стоит отметить (рис. 5), что при использовании ПО «СМРА8» и геоинформационного продукта ArcView 3.2 для визуализации данных реанализа, пользователям необходимо провести расчеты стандартных отклонений, исходя из которых осуществить выбор количества градаций абсолютных значений ПТВ.
Заключение
В работе показано, что данные климатических архивов ПТВ (ЕЯА-Шепт, ЫСЕР/МСАЯ) на высоте 2 метров в разной степени характеризуют пространственно-временную изменчивость среднегодовых и среднесезонных (зима, лето) вариаций за период 1979-2012 гг. относительно данных наземных наблюдений. Однако они в целом достаточно хорошо отображают изменчивость ПТВ на территории Якутии.
Данные ЕЯА-Шепт и ЫСЕР/МСАЯ могут быть использованы (с учетом уравнений регрессии) для восстановления архива недостающей информации наземных наблюдений.
Результаты, полученные в работе, могут быть дополнены и применены в ФГБУ «Якутское управление по гидрометеорологии и мониторингу окружающей среды», а также в сельскохозяйственной деятельности
и при построении модели прогноза изменяющегося климата на территории Якутии.
Л и т е р а т у р а
1. Uppala S. M., Kallberg P. M., Simmons A. J., Andrae U. et. al. The ERA-40 re-analysis // Quart. J. R. Met. Soc. 2005. - V. 131, № 612. - P. 2961-3012.
2. ECMWF Newsletter No. 110 - Winter 2006/07. / Riddaway B. (ed.). Reading: European Centre for Medium-Range Weather Forecast. 2007. - 53 p.
3. Berrisford, P., et. al. The ERA-Interim Archive. ERA Report Series. 1. Technical Report. European Centre for MediumRange Weather Forecasts, Shinfield Park, Reading, 2009. - Р. 16.
4. Kistler R., Kalnay E., Collins W., Saha S. et. al. The NCEP-NCAR 50-Year Reanalysis: Monthly Means CD-ROM and Documentation // Bull. Amer. Met. Soc. 2001. - V. 82, № 2. - P. 247-267.
5. Kanamitsu M., Ebisuzaki W., Woollen J., Yang S. K. et. al. NCEP-DOE AMIP-II Reanalysis (R-2) // Bull. Amer. Met. Soc. 2002. - V. 83, № 11. - P. 1631-1643.
6. Onogi K., Tsuitsui J., Koide H., Sakamoto M. et. al. The JRA-25 Reanalysis // J. Met. Soc. Japan. 2007. - V. 85, № 4. -P. 369-432.
7. Жильцова Е. Л., Анисимов О. А. О точности воспроизведения температуры и осадков на территории России глобальными климатическими архивами // Метеорология и гидрология. - 2009. - № 10. - C. 79-89.
8. Золотов С. Ю. и др. Сравнение данных реанализа NCEP/NCAR профилей температуры почвы с данными
измерений сети станций на территории Западной Сибири // Криосфера Земли. - 2011. - Т. XV, № 2. - C. 14-20.
9. Lambert S. and Boer G., 2001. CMIP1 evaluation and intercomparison of coupled climate models. Climate Dynamics. -Vol. 17. - Р. 83-106.
10. Гаврилова М. К. Климат Центральной Якутии. -Якутск: Якутское кн. изд-во, 1962. - 63 c.
11. Гаврилова М. К. Климат Центральной Якутии. Изд. 2-е, перераб. и доп. - Якутск: Якутское кн. изд-во, 1973. - 119 c.
12. Гаврилова М. К. Климат холодных регионов Земли: Учебное пособие. - Якутск: ИМЗ СО РАН, 1998. - 206 c.
13. Гаврилова М. К. Изменения климата (температуры воздуха и осадки) на территории Якутии и возможности его влияния на сельское хозяйство // Наука и образование. -Якутск: Изд-во ЯНЦ СО РАН, 2009. - № 3. - C. 48-54.
14. Зукерт Н. В., Замолодчиков Д. Г. Изменения температуры воздуха и осадков в тундровой зоне России // Метеорология и гидрология. - 1997. - № 8. - С. 45-52.
15. Груза Г. В., Ранькова Э. Я. Колебания и измерения климата на территории России // Известия РАН. Физика атмосферы и океана. 2003. - Т. 39, № 2 - C. 166-185.
16. Гаврилова М. К. Изменение температуры воздуха на территории Азиатской России в XIX-XX веках // Наука и образование. - Якутск: Изд-во ЯНЦ СО РАН, 2008. - № 4 (52). - С. 40-44.
17. Smith T. M., Peterson T. C., Lawrimore J. H., and Reynolds R. W. New surface temperature analyses for climate monitoring. Geophys. Res. Lett., 2005. - Vol. 32, L14712, doi:10.1029/2005GL023402.
18. Vose R. S., Easterling D. R., and Gleason B. Maximum and minimum temperature trends for the globe: An update through 2004. Geophys. Res. Lett., - 2005. Vol. 32, L23822, doi:10.1029/2004GL024379.
19. Brohan P., Kennedy J. J., Harris I., Tett S. F. B., and Jones P. D. Uncertainty estimates in regional and global observed temperature changes: A new dataset from 1850. J. Geophys. Res., 2006. - Vol. 111, D12106, doi:10.1029/2005JD006548.
20. Соловьев В. С., Васильев М. С. Свидетельство о государственной регистрации программы для ЭВМ №20126103912 «Cloud Map Processing Script (CMPAS)», выдан Федеральной службой по интеллектуальной собственности, патентам и товарным знакам (Роспатент) 10 января 2012 г. Правообладатель: Учреждение Российской академии наук Институт космофизических исследований и аэрономии им. Ю. Г. Шафера СО РАН.
R e f e r e n c e s
1. Uppala S. M., Kallberg P. M., Simmons A. J., Andrae U. et. al. The ERA-40 re-analysis // Quart. J. R. Met. Soc. - 2005. -V. 131, № 612. - P. 2961-3012.
2. ECMWF Newsletter No. 110 - Winter 2006/07. / Riddaway B. (ed.). Reading: European Centre for Medium-Range Weather Forecast. 2007. - 53 p.
3. Berrisford, P., et. al. (2009). The ERA-Interim Archive. ERA Report Series. 1. Technical Report. European Centre for Medium-Range Weather Forecasts, Shinfield Park, Reading. - P. 16.
4. Kistler R., Kalnay E., Collins W., Saha S. et. al. The
NCEP-NCAR 50-Year Reanalysis: Monthly Means CD-ROM and Documentation // Bull. Amer. Met. Soc. - 2001. - V. 82, № 2. - P. 247-267.
5. Kanamitsu M., Ebisuzaki W., Woollen J., Yang S. K. et. al. NCEP-DOE AMIP-II Reanalysis (R-2) // Bull. Amer. Met. Soc. 2002. - V. 83, № 11. - P. 1631-1643.
6. Onogi K., Tsuitsui J., Koide H., Sakamoto M. et. al. The JRA-25 Reanalysis // J. Met. Soc. Japan. 2007. - V. 85, № 4. - P. 369-432.
7. Zhil'cova E. L., Anisimov O. A. O tochnosti vosproizvedenija temperatury i osadkov na territorii Rossii global'nymi klimaticheskimi arhivami // Meteorologija i gidrologija, 2009. - № 10. - S. 79-89.
8. Zolotov S. Ju. i dr. Sravnenie dannyh reanaliza NCEP/ NCAR profilej temperatury pochvy s dannymi izmerenij seti stancij na territorii Zapadnoj Sibiri // Kriosfera Zemli, 2011. -T. XV, № 2. - S. 14-20.
9. Lambert S. and Boer G., 2001. CMIP1 evaluation and intercomparison of coupled climate models. Climate Dynamics. -Vol. 17. - P. 83-106.
10. Gavrilova M. K. Klimat Central'noj Jakutii. Jakutsk: Jakutskoe kn. izd-vo, 1962. - 63 s.
11. Gavrilova M. K. Klimat Central'noj Jakutii. Izd. 2-e, pererab. i dop. - Jakutsk: Jakutskoe kn. izd-vo, 1973. - 119 s.
12. Gavrilova M. K. Klimat holodnyh regionov Zemli: Uchebnoe posobie. - Jakutsk: IMZ SO RAN, 1998. - 206 s.
13. Gavrilova M.K. Izmenenija klimata (temperatury vozduha i osadki) na territorii Jakutii i vozmozhnosti ego vlijanija na sel'skoe hozjajstvo // Nauka i obrazovanie. - Jakutsk: Izd-vo JaNC SO RAN, 2009. - № 3. - S. 48-54.
14. Zukert N. V., Zamolodchikov D. G. Izmenenija temperatury vozduha I osadkov v tundrovoj zone Rossii // Meteorologija i gidrologija, 1997. - № 8. - S. 45-52.
15. Gruza G. V., Ran'kova Je. Ja., 2003. Kolebanija i izmerenija klimata na territorii Rossii. Izvestija RAN. Fizika atmosfery i okeana. - T. 39, № 2. - S. 166-185.
16. Gavrilova M. K. Izmenenie temperatury vozduha na territorii Aziatskoj Rossii v XIX-XX vekah // Nauka i obrazovanie. Jakutsk: Izd-vo JaNC SO RAN, 2008. - № 4 (52). - S. 40-44.
17. Smith T. M., Peterson T. C., Lawrimore J. H., and Reynolds R. W. New surface temperature analyses for climate monitoring. Geophys. Res. Lett., 2005. - Vol. 32, L14712, doi:10.1029/2005GL023402.
18. Vose R. S., Easterling D. R., and Gleason B. Maximum and minimum temperature trends for the globe: An update through 2004. Geophys. Res. Lett., 2005. - Vol. 32, L23822, doi:10.1029/2004GL024379.
19. Brohan P., Kennedy J. J., Harris I., Tett S. F. B., and Jones P. D. Uncertainty estimates in regional and global observed temperature changes: A new dataset from 1850. J. Geophys. Res., 2006. - Vol. 111, D12106, doi:10.1029/2005JD006548.
20. Solov'ev V. S., Vasil'ev M. S. Svidetel'stvo o gosudarstvennoj registracii programmy dlija JeVM №20126103912 «Cloud Map Processing Script (CMPAS)», vydan Federal'noi sluzhboj po intellektual'noj sobstvennosti, patentam i tovarnym znakam (Rospatent) 10 janvarja 2012 g. Pravoobladatel': Uchrezhdenie Rossijskoj akademii nauk Institut kosmofizicheskih issledovanij i ajeronomii im. Ju. G. Shafera SO RAN.