Известия РАН. Физика атмосферы и океана, 2023, T. 59, № 2, стр. 200-216

Вариации приземных концентраций и общего содержания СО2 и СН4 в центральной части Европейской территории России

К. Н. Вишератин a***, Е. Л. Баранова a, Г. И. Бугрим a, В. Н. Иванов a, Е. И. Краснопеева a, Д. Г. Сахибгареев a, В. П. Устинов a, А. В. Шилкин ab

a ФГБУ “НПО “Тайфун”,
249038 Обнинск, ул. Победы, 4, Россия

b Институт физики атмосферы им. А.М. Обухова РАН
119017 Москва, Пыжевский пер., 3, Россия

* E-mail: kvisher@rpatyphoon.ru
** E-mail: kvisher@mail.ru

Поступила в редакцию 30.08.2022
После доработки 14.11.2022
Принята к публикации 11.01.2023

Полный текст (PDF)

Аннотация

Представлены результаты измерений методом ИК-спектроскопии приземных концентраций и общего содержания углекислого газа и метана в атмосфере центральной части ЕТР на ст. Обнинск. Приведено описание комплекса аппаратуры MР-32/МГС для мониторинга газового состава атмосферы. Проанализированы межгодовые и сезонные вариации и тренды приземных концентраций СО2 и CH4 за 1998–2021 гг. и общего содержания за 2015–2021 гг. По результатам кросс-корреляционного вейвлетного анализа фаза годовых вариаций средних по высоте концентраций CH4 опережает фазу приземных вариаций на 2–3 месяца, а вариации средних по высоте концентраций СО2 отстают по фазе от годовых вариаций приземных значений на 1–2 месяца. Минимальные приземные концентрации метана в мае-августе совпадают со средними концентрациями в толще атмосферы, а для углекислого газа в этот же период они ниже на 20–30 ppm. В зимние месяцы минимальные приземные концентрации метана и углекислого газа выше средних по высоте на 70–150 ppb и на 6–15 ppm соответственно. Результаты измерений сопоставлены со спутниковыми данными GOSAT и данными наземных Европейских станций.

Ключевые слова: углекислый газ, метан, приземные концентрации, средние по высоте концентрации, вейвлетный анализ, сезонные вариации, тренды

1. ВВЕДЕНИЕ

Углекислый газ (СО2) и метан (СН4) являются основными атмосферными газами, влияющими на изменение климата. Рост глобальной концентрации углекислого газа и метана в 1750–2019 гг. составил 47–49% и 156–162% соответственно, при этом вклад в усиление парникового эффекта за счет СО2 и СН4 оценивается в 60–80% [Кондратьев, Крапивин, 2004; IPCC, 2021; WMO, 2021]. Сжигание ископаемого топлива (уголь, природный газ и нефть) является основным (более 70%) источником углекислого газа. Другие источники – сжигание деревьев и других биологических материалов, дыхание животных и человека, а также результат определенных химических реакций (например, при производстве цемента) [IPCC, 2021]. Двуокись углерода удаляется из атмосферы в основном при поглощении растениями в рамках биологического цикла углерода. Помимо нагрева тропосферы, увеличение содержания CO2 приводит к радиационному охлаждению стратосферы [Кондратьев, Крапивин, 2004]. Метан является вторым по важности антропогенным парниковым газом после двуокиси углерода. По оценкам, приведенным в [IPCC, 2021; WMO, 2021], антропогенные источники метана (примерно 62%) связаны с животноводством и другими методами ведения сельского хозяйства и распадом органических отходов на свалках бытовых отходов. Около 30% выделяется при добыче и транспортировке угля, природного газа и нефти. Заболоченные земли являются основным (более 80%) естественным источником метана. Метан имеет значительно более короткое время жизни, чем CO2 и, следовательно, стабилизация или сокращение выбросов CH4 в течение нескольких десятилетий приводит к стабилизации или уменьшению его концентрации в атмосфере. Метан также важен для химии тропосферы и стратосферы, реакции CH4 с гидроксильным радикалом ОН являются основным стоком метана [Кондратьев, Крапивин, 2004]. Следует заметить, что по последним оценкам NOAA в 2020–2021 гг. был зафиксирован рекордный глобальный рост CO2 и CH4 [NOAA, 2021].

Одной из актуальных задач современности является развитие сети регулярных наблюдений за содержанием парниковых газов. В настоящее время мониторинг приземных концентраций СО2 и СН4 осуществляется на сети фоновых и региональных станций. На ряде станций проводится периодический забор проб воздуха во фляги и дальнейший анализ на стационарных установках с использованием прецизионных калибровочных смесей. Наряду с методом забора проб используются также методы непрерывных измерений различными газоанализаторами. Данные измерений проходят предварительную фильтрацию и помещаются в открытый доступ с задержкой 1–2 года [Dlugokencky et al., 2021; WDCGG, 2021; ESRL, 2021]. Важную информацию о вариациях общего содержания и средней по высоте концентрации парниковых газов позволяют получить интегральные методы, основанные на регистрации солнечного излучения. Систематические измерения общего содержания ряда атмосферных газов приборами высокого спектрального разрешения (в основном Bruker IFS-125 HR) в настоящее время проводятся на станциях NDACC (Network for the Detection of Atmospheric Composition Change) и TCCON (Total Carbon Column Observation Network) [Ostler et al., 2014, Buschmann et al., 2016]. В последние годы для наблюдений за общим содержанием CO2 и CH4 на сети станций COCCON (COllaborative Carbon Column Observing Network) начали применяться портативные Фурье-спектрометры Bruker EM27/SUN [Frey et al., 2019]. В России спектроскопические измерения общего содержания ряда атмосферных газов были начаты в 1960–1980 гг. в Главной геофизической обсерватории (ГГО) [Шашков и др., 1985], Санкт-Петербургском государственном университете [Мироненков и др., 1996], Институте физики атмосферы (Москва) [Дианов–Клоков, Юрганов,1982], Институте оптики атмосферы (Томск) [Кабанов и Сакерин, 1995].

Полевой спектрометрический комплекс и его модификации [Арефьев и др., 1978; Арефьев и др., 1995] с 1980 г. применялись в Институте экспериментальной метеорологии (Обнинск) для измерений общего содержания атмосферных газов в Киргизии (ст. Иссык-Куль), Антарктиде (ст. Новолазаревская) и др. местах [Арефьев и др., 2010]. В 1998 г. на станции Обнинск были начаты регулярные измерения концентрации диоксида углерода и метана в приземном воздухе [Baranov et al., 2000]. С 2015 г. в НПО “Тайфун” введен в эксплуатацию газоаналитический комплекс MP-32/МГС, и начаты регулярные измерения метана, углекислого газа, закиси азота (N2O) и окиси углерода (CO) в пробах приземного воздуха, а также общее содержание в столбе атмосферы перечисленных компонентов, включая дополнительно водяной пар (H2O) и озон (O3) [Visheratin et al., 2021].

В настоящей статье представлено описание комплекса аппаратуры для мониторинга газового состава атмосферы MP-32/МГС, а также результаты измерений приземных концентраций диоксида углерода и метана в 1998–2021 гг. Более подробно, с использованием кросс-корреляционного вейвлетного анализа и других методов, рассмотрены сезонная и долговременная изменчивость и взаимосвязь между приземными и средними в толще атмосферы концентрациями углекислого газа и метана в период с 2015 по 2021 гг.

2. АППАРАТУРА И МЕТОДИКА ИЗМЕРЕНИЙ

2.1. Спектрометрический комплекс MР-32/МГС

Общий вид комплекса аппаратуры MР-32/МГС показан на рис. 1. Комплекс предназначен для мониторинга вариаций малых газовых составляющих (МГС) атмосферы и используется в двух режимах:

Рис. 1.

Блок-схема комплекса МР-32/МГС и его общий вид. 1 – спектрометр ФСМ 2203, 2 – кювета КГ40, 3 – солнцеследящая система ССС-КГМА-С, 4 – управляющий компьютер, 5 – система откачки и заполнения кюветы, 6 – глобар, 7 – приемники излучения, 8 – перемещаемые зеркала.

– режим регистрации спектров пропускания пробами воздуха в оптической многоходовой кювете;

– режим регистрации спектров поглощения солнечного излучения атмосферой.

Для перехода от режима измерений концентраций газов в приземном слое к измерениям общего содержания газов в столбе атмосферы изменяется положение подвижных зеркал (8), расположенных внутри Фурье-спектрометра ФСМ 2203 (разработка ООО “Инфраспек”, [Infraspek, 2021]). Спектральный диапазон канала регистрации прошедшего через кювету излучения 1700–5000 см–1, канала регистрации солнечного излучения 800–7700 см–1, максимальная разность хода в интерферометре 41 мм, аподизация Нортон-Бир слабая. Спектральная ширина инструментального контура на уровне 50% высоты 0.12–0.15 см–1, предел допустимой погрешности измерений интенсивности солнечного излучения и излучения искусственного источника (глобар) не более 3.0%.

Многоходовая оптическая кювета (КГ40) [Infraspek, 2021] обеспечивает получение поглощающих слоев газа толщиной до 40 м. Технологическая система осуществляет вакуумирование и заполнение кюветы воздухом от пробоотборника, расположенного на высоте 8 м от земной поверхности. В случае измерений общего содержания, солнцеследящая система ССС-КГМА-С, расположенная на крыше здания, сканирует положение Солнца и направляет во входную аппертуру спектрометра поток солнечного излучения диаметром 80 мм. Скорость пространственного сканирования по азимуту и углу места в режиме слежения от 0.25 до 2.5 угловых минут/с. Абсолютная погрешность точного наведения на Солнце и последующего слежения за Солнцем не более 2 угловых минут. Одно измерение занимает 2.5 минуты, в течение которых записывается одна спектрограмма солнечного излучения.

2.2. Методика определения приземных концентраций

Методика измерений приземных концентраций углекислого газа и метана в период с 1998 по 2014 гг. приведена в работах [Baranov et al., 2000; Арефьев и др., 2004; Kashin, 2013]. Концентрации СН4 и CO2 в пробах воздуха определялись по ИК-спектрам поглощения, полученным с помощью Фурье-спектрометров Perkin-Elmer 1720 или Bruker IFS-113V. Разработанная и изготовленная в НПО “Тайфун” кювета обеспечивала поглощающий слой 30 м при базовой длине 1 м. Для определения концентраций СО2 использовался спектральный интервал 2220–2280, для СН4 – интервал 2940–3110 см–1. При систематической погрешности менее 2% случайная погрешность составляла 30 ppb для СН4 и 2 ppm для СО2.

Методика измерений приземных концентраций углекислого газа и метана комплексом МР-32/МГС в основном аналогична использованной ранее [Baranov et al., 2000; Арефьев и др., 2004;Kashin, 2013]. Пропускание рассчитывается как отношение спектра интенсивности излучения, прошедшего через кювету с поглощающим слоем 40 м, к усредненному спектру интенсивности излучения, прошедшего через вакуумированную кювету до и после забора пробы. Остаточное давление в вакуумированной кювете составляет менее 0.5 мб. Одновременно регистрируется температура внутри и вне кюветы. Один цикл измерений составляет около 20 минут.

Разработанные ранее программы адаптированы применительно к спектральному разрешению и шагу спектрограммы ФСМ 2203, уточнены параметры аппаратной функции, сделаны оценки реального спектрального разрешения и качества интерферограмм. Рабочие спектральные интервалы выбраны на основе анализа полос поглощения СО2 и СН4. Основные критерии – степень перекрывания линий поглощения определяемых и мешающих газов при максимально возможных вариациях концентраций измеряемых газов и длине оптического пути 40 м. Для измерений СО2 выбраны микроокна 2217.3–2234.5 и 2242.2–2249.1 см–1, а для измерений СН4 использовались микроокна 2946.93–2949.28; 2956.81–2959.52; 2967.48–2970.07; 3014.37–3020.15; 3037.69–3039.38; 3054.44–3058.84; 3075.59–3078.79; 3085.00–3087.83; 3103.62–3105.49 см–1. Концентрация основного мешающего газа (Н2О) определяется в интервалах 2965–2967 и 3024–3028 см–1 и используется в качестве первого приближения при расчете концентраций СО2 и СН4. Стандартное отклонение измеренного спектра от рассчитанного по параметрам спектральных линий HITRAN 2012 [Rothman, 2013] минимизировалось попеременным варьированием концентраций СО2, СН4 и водяного пара. Величины их концентраций, полученные при наилучшем согласии экспериментального и расчетного спектров, принимались за результат измерения.

Для оценки систематических погрешностей проведен ряд калибровочных измерений концентраций газов в многоходовой кювете. Для этого использовались баллоны со стандартными смесями ФГУП ВНИИМ им. Менделеева с погрешностью отношения смеси СО2 3–4%, СH4 4–5%. В конце 2018 г. была проведена поверка этих баллонов в Лаборатории мониторинга малых газовых компонент атмосферы ГГО им А.И. Воейкова, что позволило снизить погрешность калибровочных смесей до 0.02%. В 2021 г. были проведены совместные измерения СО2 и СН4 в приземном слое атмосферы приборами МР-32/МГС и Picarro G2311-f. Принцип действия Picarro G2311-f основан на методе CRDS (cavity ring-down spectroscopy) с использованием полупроводникового ИК лазера и газовой кюветы с эффективной длиной пути лучей более километра [Chen et al., 2010]. Результаты сопоставления в качестве примера приведены на рис. 2 для периода с 9 по 12 августа 2021 г. Расхождение между совпадающими по времени измерениями МР-32/МГС и Picarro G2311-f менее 40 ppb для СН4 и менее 2 ppm для СО2.

Рис. 2.

Вариации приземных концентраций углекислого газа и метана 9–12 августа 2021 г. 1 – Picarro G2311-f, 2 – МР‑32/МГС.

По результатам измерений калибровочных смесей систематическая погрешность определения СН4 составляет менее 2%, а СО2 менее 1%. Основной вклад в систематическую погрешность вносит температурный дрейф. Случайная погрешность оценивалась при повторных анализах одной и той же калибровочной смеси и равна 1.5 и 0.5% для СН4 и СО2 соответственно.

2.3. Методика определения средних по высоте концентраций

Предварительный анализ спектров солнечного излучения показал, что узкие спектральные интервалы в средней инфракрасной области, используемые на станциях NDACC [Ostler et al., 2014; Buschmann et al., 2016], не вполне подходят для определения общего содержания газов прибором МР-32/МГС. Это связано с тем, что спектральное разрешение спектрометра ФСМ 2203 значительно ниже, чем у приборов Bruker IFS-125. Поэтому поглощение в центрах линий оказывается меньшим, а из-за большей полуширины линий происходит перекрывание с соседними спектральными линиями. По мере накопления данных измерений МР-32/МГС в 2016–2018 гг. были проведены сравнительные расчеты общего содержания СН4 и СО2 в различных спектральных интервалах. При выборе узких спектральных интервалов (микроокон) общей проблемой является перекрывание линий поглощения определяемых и мешающих газов, а также влияние на структуру спектров водяного пара, содержание которого существенно меняется с сезоном. В настоящей работе приведены результаты вычислений средних по высоте концентраций CH4 и CO2 для широких спектральных интервалов 5897–6145 и 6173–6390 см–1, используемых на станциях TCCON [Ostler et al., 2014, Buschmann et al., 2016] и C-OCCON [Frey et al., 2019]. Эти интервалы являются рабочими также для Фурье-спектрометра TANSO-FTS, входящего в состав GOSAT (Greenhouse gases Observing SATellite) на борту японского спутника “Ibuki” [GOSAT, 2021].

Вместе с развитием дистанционных методов определения вертикальных профилей и общего содержания атмосферных газов получили развитие методы решения обратных задач, реализованные в таких программах как SFIT, GFIT, PROFIT9 [Тимофеев, 2016]. Для расчетов общего содержания газов нами использовался свободно распространяемый пакет программ SFIT4 V0.9.4.4 [SFIT4]. В качестве атласа спектральных линий был использован атлас HITRAN-2012 [Rothman et al., 2013], априорные среднеширотные профили температуры, давления и профилей газов заимствовались из [WACCAM, 2013]. В используемой версии SFIT4 при решении обратной задачи для определения общего содержания углекислого газа или метана восстанавливался вертикальный профиль этих газов (48 слоев, 0.02–120 км) и в дальнейшем проводилось суммирование. Итеративная процедура приближения экспериментальных и расчетных значений основана на алгоритме Левенберга–Марквардта. Переход от общего содержания газа в столбе атмосферы (осГАЗ, молекул/см2) к его средней по высоте мольной концентрации (хГАЗ, ppm, ppb), осуществлялась с учетом общего содержания водяного пара, т.е. данные измерений пересчитывались для “сухого” воздуха. Средняя по высоте мольная концентрация газа хГАЗ определяется делением общего содержания газа осГАЗ на общее содержание сухого воздуха Рсв [Buschmann et al., 2016]:

(1)
${\text{хГАЗ}} = {\text{осГАЗ/}}Рсв.$

Величина Рсв находится конвертированием приземного давления Р в общее содержание столба воздуха при известной величине общего содержания водяного пара:

(2)
$Р{\text{св}} = Р/\left( {{{m}_{{air}}}g} \right) - \left( {{{m}_{{{{{\text{H}}}_{{\text{2}}}}{\text{O}}}}}/{{m}_{{{\text{air}}}}}} \right){\text{ос}}{{{\text{H}}}_{{\text{2}}}}{\text{O}},$
где Р – приземное давление в Па,

mair – молекулярная масса сухого воздуха (28.96 × × 10–3 Na кг/молекула),

${{m}_{{{{{\text{H}}}_{{\text{2}}}}{\text{O}}}}}$ – молекулярная масса водяного пара (18 × × 10–3 Na кг/молекула),

осH2O – общее содержание в столбе атмосферы водяного пара (молекул/см2),

Na – число Авогадро (6.022 × 1023 молекул/моль),

g – ускорение силы тяжести в месте измерений.

Для вычислений общего содержания водяного пара был выбран интервал 4680–4710 см–1, содержащий одиночную линию Н2О с волновым числом 4699 см–1. В расчетах осН2О учитывалось поглощение в “крыльях” линий N2O и CH4, а положение центра линии поглощения водяного пара использовалось также для предварительной оценки возможного смещения шкалы волновых линий спектрометра. При определении общего содержания CO2 в анализ включались мешающие газы CH4 и Н2О, а при вычислениях CH4 мешающими газами являлись CO2 и Н2О. Для примера, на рис. 3а приведен дневной ход CO2, CH4 и Н2О для дневных безоблачных условий в течение 7 июля 2021 г., а на рис. 3б показаны измеренные солнечные спектры и разность между измеренными и рассчитанными спектрами в полдень. На рис. 3а для всех газов заметен тренд, исключение которого дает значения стандартных отклонений 1.1 ppm, 3.9 ppb и 0.03 г/см2 для хСО2, хСН4 и осН2О соответственно. Полученные оценки характеризуют случайную погрешность среднедневных величин. Аналогичные оценки для других сезонов из-за погодных условий основаны на более коротких сериях измерений и в среднем дают величины погрешности среднедневных значений для углекислого газа 0.7% (~3 ppm) и метана 1.6% (~30 ppb).

Рис. 3.

Результаты измерений 7 июля 2021 г. средних по высоте концентраций хСО2, хСН4 и общего содержания осН2О. (а) – внутридневной ход, приведены средние за день (mean) значения и их стандартное отклонение (std), (б) – атмосферные солнечные спектры и разность между измеренными и рассчитанными спектрами. Приведены также стандартные отклонения разности измеренных и рассчитанных спектров, отношение сигнал/шум (SNR) и число степеней свободы (DOF).

В результате работы программы SFIT для каждого измерения вычисляются отношение сигнал/шум (SNR) и число степеней свободы (DOF). Для спектров, представленных на рис. 3б, значения SNR и DOF для CO2 и СН4 примерно соответствуют средним значениям в летний период. Отношение сигнал/шум испытывает значительную изменчивость и растет в зимнее время, что связано с ростом оптического пути при больших зенитных углах. Величины (1/SNR)2 обычно используются для оценки случайной погрешности (измерительного шума) единичного измерения [Takeda et al., 2021]. Для значений SNR из рис. 3б расчет случайной ошибки единичного измерения дает величины 1.1% для CO2 и 1.9% для СН4. Систематические ошибки, связанные с неточностями спектроскопических параметров, составляют для CO2 2.9% (~12 ppm) и для СН4 6.6% (~125 ppb). Для расчетов систематических ошибок были взяты максимальные значения погрешностей интенсивностей и полуширин линий из [Rothman et al., 2013].

3. РЕЗУЛЬТАТЫ И ОБСУЖДЕНИЕ

3.1. Концентрации СН4 и СО2 в приземном слое в 1998–2021 гг.

Для расположенной в центре ЕТР ст. “Обнинск” общей тенденцией является рост парниковых газов в приземном слое с сезонными минимумами в летний период и максимумами – в зимний [Арефьев и др. 2014; Арефьев и др. 2015]. В течение суток максимумы концентрации метана и углекислого газа наблюдаются обычно в утренние часы. Эффект накопления усиливается в те дни, когда в ночное время суток при отсутствии конвективных потоков воздуха возникает инверсионный температурный слой, что приводит к повышенным концентрациям метана и углекислого газа и суточному ходу с минимумом в дневное время (см. в качестве примера рис. 2). В период наблюдений иногда регистрировались аномально высокие концентрации СН4, обычно им сопутствовали также повышенные концентрации СО2 и других измеряемых компонентов (СО, N2O), что указывает на перенос загрязненных масс воздуха от удаленных источников. Например, в 2021 г. значительные концентрации метана (2500–4000 ppb), углекислого газа (470–500 ppm) и окиси углерода (более 500 ppb) отмечались 15–19 января. Расчеты обратных траекторий с помощью модели [HYSPLIT, 2021] показали, что воздушные массы в этот период пришли преимущественно с северо-восточного направления, возможно, из района Москвы.

Динамика многолетней изменчивости средних годовых концентраций СО2 и СН4 и усредненный за 1998–2021 гг. сезонный ход приведены на рис. 4 с оценкой стандартных отклонений результатов измерений. Небольшой отрицательный тренд СН4, наблюдавшийся в 1998–2007 гг., в дальнейшем сменился положительным трендом (рис. 4а). Рост концентрации метана после 2007 г., хотя и имеет региональные особенности [Ракитин и др. 2017; Rakitin, et al., 2020], но носит глобальный характер [WMO, 2021]. Причины, приведшие к изменению тренда метана, в ряде работ объясняются влиянием глобального потепления на рост эмиссии метана на Арктическом континентальном шельфе, см. например, [Dlugokencky et al., 2011]. Некоторые авторы связывают начавшийся в 2007 г. рост метана c увеличением темпов добычи сланцевого газа в Северной Америке (с 2005 по 2015 гг. производство сланцевого газа выросло в 14 раз, подробнее см. [Howarth, 2019]).

Рис. 4.

Вариации средних годовых значений СН4 (а), СО2 (в) в приземном слое атмосферы и средний сезонный ход (б, г). 1 – Обнинск, 2 – Хегихатсал.

Концентрация углекислого газа за время наблюдений также значительно возросла (рис. 4в). Увеличение концентрации СО2 на ст. Обнинск в процентном отношении (0.65% в год) в полтора раза превосходит наблюдавшийся с 2007 г. рост СН4 (0.42% в год). Для примера, на рис. 4 представлены также данные станции мониторинга парниковых газов Венгерской метеорологической службы HUN [ESRL, 2021]. Станция расположена в западной части Среднедунайской низменности вблизи селения Хегихатсал (46.57° N, 16.38 °E) [Haszpra et al., 2019]. Выбор этой станции для сопоставления обусловлен тем, что Альпы на западе и Карпаты на востоке определяют для этой станции преимущественно северное и южное направление ветров при небольшом влиянии местных сельскохозяйственных источников. Также, как и для ст. Обнинск, максимальные значения концентрации метана и углекислого газа на станции Хегихатсал приходятся на зимние месяцы, а минимальные на летние (рис. 4б, 4г). Амплитуда сезонного хода метана на станции Хегихатсал сравнима с амплитудой вариаций на ст. Обнинск, хотя абсолютные значения ниже в среднем на 70 ppb. Для углекислого газа сезонные вариации на станциях Обнинск и Хегихатсал практически совпадают, а тренды, амплитуда колебаний и среднее значение концентрации СО2 на этих станциях показывают достаточно хорошее согласие (табл. 1).

Таблица 1.  

Коэффициенты линейного тренда, амплитуды колебаний и средние значения (в скобках) приземных концентраций метана и углекислого газа

Станция/годы Тренд СН4 Тренд СО2 Амплитуда (Среднее)
ppb/год ppm/год СН4, ppb СО2, ppm
1998–2020 2007–2020 1998–2020 1998–2020 1998–2020
Обнинск 2.8 ± 1.1 8.3 ± 2.6 2.6 ± 0.6 49 (1997) 11 (397)
Хегихатсал 4.6 ± 1.0 7.5 ± 2.1 2.2 ± 0.5 43 (1940) 13 (393)

В табл. 1 дополнительно приведены тренды метана для периода роста с 2007 г. Для сравнения, глобально усредненные тренды в период с 2007 по 2020 гг., рассчитанные по данным [ESRL, 2021], составляют 7.4 ± 2.1 ppb/год и 2.0 ± 0.4 ppm/год для метана и углекислого газа соответственно.

3.2. Средние по высоте концентрации метана и углекислого газа в 2015–2021 гг.

3.2.1. Сопоставление со спутниковыми данными GOSAT

GOSAT (Greenhouse Gases Observing Satellite) является первым спутником, предназначенным для наблюдения за парниковыми газами из космоса. С 2009 г. GOSAT проводит глобальные наблюдения за СО2, СН4 и Н2О с пространственным разрешением ~10 км2, временем пересечения экватора около 13:00 по местному времени и 3-дневным циклом повторения. Фурье-спектрометр TANSO-FTS на борту GOSAT измеряет коротковолновое инфракрасное излучение (SWIR), отраженное от поверхности Земли в полосах поглощения кислорода (12950–13200 см–1), углекислого газа (6180–6380 см–1) и метана (5900–6150 см–1) со спектральным разрешением 0.2 см–1. Точность единичного измерения менее 3.5 ppm для хСО2 и 15 ppb для хСН4 [Yoshida et al., 2011; Kivimäki et al., 2019].

Первичные данные GOSAT (версия V02.хх) содержатся в [GOSAT, 2021]. Для формирования временных рядов хСО2 и хСН4 над наземными пунктами наблюдений обычно проводят пространственное усреднение данных GOSAT для ячеек 6° × 6° и более [Kivimäki et al., 2019, Oshchepkov et al., 2013, Gavrilov et al., 2014]. Нами были проанализированы 3 варианта усреднения для ячеек 3° × 3°, 5° × 5°, 10° × 10°. При размерах ячеек меньших 5° × 5° число совпадающих значений в интервале ±1 час относительно времени пролета спутника становится малым, и из сравнения могут выпадать целые сезоны. Сопоставление с данными GOSAT для области 5° × 5°, приведено на рис. 5.

Рис. 5.

Сопоставление наземных измерений хСО2 (а) и хСН4 (б) на ст. Обнинск (55.11N, 36.60E) прибором МР-32/МГС и данных GOSAT, усредненных для ячейки 5° × 5°. На рисунке приведено число совпадающих измерений (N), коэффициент корреляции (R2), отношение средних значений GOSAT/МР-32 и стандартное отклонение (в скобках) для 95% доверительного интервала.

Сравнение данных наземных и спутниковых измерений xСН4 и xСО2 показало, что результаты измерений на ст. Обнинск систематически завышены. Отношения средних значений GOSAT/МР-32 для ячеек 5° × 5° и 10° × 10° в пределах погрешности совпадают и в среднем составляют 0.98 ± 0.01 и 0.97 ± 0.01 для xСО2 и xСН4 соответственно. Для компенсации систематического расхождения и привязки данных МР-32/МГС к шкале GOSAT, нами были введены поправочные коэффициенты, равные 0.98 для углекислого газа и 0.97 для метана. Следует отметить, что введение систематических поправок не сказывается на величинах трендов и амплитуде годового хода, также как и на результатах спектрального или кросс-корреляционного анализа.

3.2.2. Сопоставление с данными измерений на станциях TCCON

Среднемесячные значения xСН4 и xСО2, полученные с помощью прибора МР-32/МГС в 2015–2021 гг., были сопоставлены со среднемесячными данными станций TCCON. Для сравнения были использованы данные европейских станций Карлсруэ, Бремен, Орлеан и Париж, расположенных ниже 200 м над уровнем моря [TCCON, 2021]. Поскольку не для всех станций имеются данные за 2021 г., параметры сопоставляемых рядов – линейный тренд, амплитуда вариаций и среднее значение приведены в табл. 2 за одинаковый период с января 2016 по декабрь 2020 г.

Таблица 2.  

Коэффициенты линейного тренда, амплитуды колебаний и средние значения (в скобках) концентраций xСН4 и xСО2 за период 2016–2020 гг.

Станция Тренд xСН4 Тренд xСО2 Амплитуда (Среднее)
ppb/год ppm/год xСН4, ppb xСО2, ppm
Обнинск
55.11N 36.60E
12.5 ± 1.9 2.3 ± 0.5 11.5 (1843) 5.6 (409)
Карлсруэ
49.10N 8.44E
10.2 ± 1.5 2.2 ± 0.3 8.8 (1842) 3.3 (408)
Бремен
53.10N 8.85E
10.7 ± 1.6 2.4 ± 0.4 8.0 (1845) 2.5 (408)
Орлеан
47.97N 2.11E
9.8 ± 1.5 2.2 ± 0.4 10.8 (1838) 3.2 (407)
Париж
48.85N 2.36E
10.0 ± 1.5 2.0 ± 0.4 8.9 (1843) 3.5 (408)
GOSAT (Обнинск, 5° х 5°). 10.5 ± 1.4 2.3 ± 0.3 8.0 (1846) 3.6 (406)

Из табл. 2 следует, что результаты измерений на ст. Обнинск за анализируемый период достаточно хорошо согласуются с наземными наблюдениями на станциях TCCON и дополнительно приведенными в табл. 2 данными GOSAT. Коэффициенты корреляции между измерениями на ст. Обнинск и станциях TCCON превышают 0.84 для метана и 0.73 для углекислого газа при значимости, превышающей уровень 2σ.

Среднемесячные временные ряды xСН4 и xСО2 за 2015–2021 гг. по данным ст. Обнинск и GOSAT, а также по данным ст. Карлсруэ, для которой результаты измерений за 2021 г. были опубликованы в [Hase et al., 2017], приведены на рис. 6а, 6в, а средний сезонный ход за время наблюдений – на рис. 6б, 6г. Экстремумы сезонного хода xСН4 и xСО2 по данным ст. Обнинск, GOSAT и ст. Карлсруэ совпадают, при этом амплитуда сезонных вариаций на ст. Обнинск немного выше. В 2021 г. концентрации углекислого газа и метана на ст. Обнинск выше данных GOSA-T и ст. Карлсруэ, особенно это заметно для метана в октябре–ноябре 2021 г. В осенне-зимний период измерения немногочисленны и зачастую проводятся при больших зенитных углах Солнца. В частности, в октябре 2021 г. повышенные значения углекислого газа и метана были, вероятно, обусловлены резкой сменой погодных условий. Если с 6 по 8 октября давление составляло 764–765 мм рт. ст., влажность 40–50%, то 11 октября, при сохранении ясной погоды, давление упало до 740 мм рт. ст., влажность повысилась до 97%, а температура понизилась на 10°С. Вероятно, реальные профили метеопараметров в этот период существенно отличались от априорных, что могло привести к завышенным расчетным значениям xСН4 и xСО2. На данном этапе исследований результаты измерений в 2021 г. являются предварительными, более детальный анализ будет темой дальнейших работ.

Рис. 6.

Вариации среднемесячных значений xСН4 (а), xСО2 (в) в толще атмосферы и средний сезонный ход (б, г) за 2015-2021 гг. 1 – Обнинск, 2 – Карлсруэ, 3 – GOSAT (Обнинск, усреднение 5° × 5°).

3.3. Корреляционные связи и фазовые соотношения между приземными и средними в толще концентрациями

Измерения концентраций метана и углекислого газа в приземном слое и в толще атмосферы на ст. Обнинск проводятся практически одновременно. Это позволяет оценить взаимосвязи между приземными и средними в толще концентрациями и провести сравнительный анализ сезонной и долговременной изменчивости, а также сделать оценки возможного влияния местных источников или переноса из других регионов. Для сокращения мы в дальнейшем будем использовать обозначения gСН4 и gСО2 для приземных, а хСН4 и хСО2 для средних по высоте мольных концентраций.

Рассмотрим динамику корреляционных связей и фазовые соотношения gСН4 и gСО2, хСН4 и хСО2 более подробно с помощью кросс-корреляционного вейвлетного анализа [Grinsted et al., 2004]. Для этого предварительно все ряды были стандартизованы и удален линейный тренд. Кросс-корреляционный вейвлетный анализ между gCH4 и gСО2, хСН4 и xСО2 (вейвлетограммы не приводятся) показал, что для приземных концентраций годовые вариации метана и углекислого газа происходят в фазе, а для средних по высоте концентраций почти в противофазе, со сдвигом фазы 4–5 мес. Этот сдвиг фаз нашел отражение и в корреляционных связях между gCH4 и хСН4, gСО2 и xСО2 (рис. 7). Значимые (доверительный интервал более 95%) корреляции характерны только для годовой гармоники. Высокий уровень значимости связан с тем, что по данным спектрального анализа (в статье не приводятся) амплитуда годовой гармоники более чем в два раза превосходит амплитуды других колебаний.

Рис. 7.

Кросс-корреляционная вейвлетограмма среднемесячных значений gCH4 и хСН4 (а), gСО2 и xСО2 (б) за 2015–2021 гг., ст. Обнинск. Цветовая шкала в отн. ед. Значения больше 2.5 отн. ед. соответствуют доверительному интервалу 65%, более 3 отн. ед. – доверительному интервалу 95% (выделен жирной линией). Направления стрелок вправо соответствуют синхронности колебаний.

Фаза годовых вариаций xCH4 опережает фазу годовых вариаций gCH4 на 2–3 месяца (стрелки на рис. 7а направлены вверх или вправо и вверх). Аналогичный анализ, за более короткий интервал с 2017 по 2020 гг., проведенный для приземных [Kubistin et al., 2021] и средних по высоте [Hase et al., 2017] концентраций на ст. Карлсруэ, также показал небольшое опережение фазы годовых вариаций хСН4 относительно gCH4 в период 2017–2018 гг. Следует отметить, что для периодов колебаний больше 28 мес. вариации метана в приземном слое и в тропосфере происходят в основном в одной фазе (при “отбеливании” рядов, т.е. исключении сезонного хода, корреляции рядов хСН4 и gCH4 с периодами более 28 мес. становятся значимыми). Это, возможно, указывает на ведущую роль наземных источников метана во внутригодовых и годовых вариациях и сравнительно быстрое, в течение 2–3 лет, затухание короткопериодных приземных пульсаций в тропосфере. При этом фазовые соотношения между приземными и средними по высоте концентрациями метана для годовых вариаций и колебаний с периодами более 2.5 лет стабильны для всего анализируемого периода.

В случае углекислого газа годовые вариации xСО2 отстают по фазе на 1–2 месяца от годовых вариаций gСО2 (стрелки на рис. 7б направлены вправо и вниз). Аналогичные результаты получены при кросс-корреляционном вейвлетном анализе данных ст. Карлсруэ. Возможно, стабильные фазовые соотношения имеются также в области квазидвухлетних колебаний, для которых, судя по рис. 7б, вариации gСО2 и хСО2 противофазны. В целом, представляется заслуживающим внимания анализ фазовых соотношений приземных и средних по высоте концентраций СН4 и СО2 на основе более длительных измерений.

3.4. Отношение концентраций СН4/CО2 в приземном слое и в толще атмосферы

Одним из методов оценки антропогенных эмиссий метана и углекислого газа и их взаимосвязи является анализ изменчивости отношения СН4/СО2 (например [Tohjima et al., 2014; Wong et al., 2015; Ghasemifard et al., 2019] и цитируемая литература). При этом исходят из предположения, что источники и стоки наблюдаемых в данном месте вариаций CO2 и CH4 могут быть достаточно удалены, а соотношение CH4/CO2 используется для фильтрации наблюдаемых данных, на которые, возможно, влияют также местные источники – стоки. Корреляция внутри и межгодовых вариаций gСН4 и gСО2 отмечалась ранее в ряде работ, например в [Yoshida et al., 2011] показано, что для различных сезонов коэффициенты регрессии СН4/СО2 для Европейского региона варьируют от 4.7 до 15.6.

Регрессии среднемесячных значений gСН4/gСО2 и средних по высоте отношений концентраций хСН4/хСО2 представлены на рис. 8. По данным ст. Обнинск наблюдается более широкий диапазон вариаций gСН4/gСО2 сравнительно с данными ст. Хегихатсал, Венгрия. При этом обращает на себя внимание совпадение коэффициентов регрессии. Это, вероятно, связано с тем, что большое время жизни метана и углекислого газа приводит к эффективному перемешиванию в пограничном слое атмосферы и наблюдаемые на ст. Хегихатсал и Обнинск вариации приземных концентраций возможно обусловлены эмиссией удаленного общего источника.

Рис. 8.

Регрессия среднемесячных значений СН4 на СО2. 1 – Обнинск, приземные концентрации, 2 – Хегихатсал, приземные концентрации, 3 – Обнинск, средние по высоте концентрации, 4 – GOSAT (Обнинск, усреднение 5° × 5°), средние по высоте концентрации. На рисунке приведены также числовые значения коэффициентов регрессии (доверительная вероятность более 95%).

Коэффициенты регрессии для средних по высоте концентраций на ст. Обнинск и по данным GOSAT меньше на 6%. Причины различия могут быть связаны как с относительным ростом тропосферного углекислого газа, так и, что более вероятно, с относительным уменьшением тропосферного метана, и требуют дальнейшего изучения. Учитывая выводы работы [Howarth, 2019] о масштабной интенсификации добычи сланцевого газа в США и Канаде, как о возможной причине роста концентрации метана в период после 2007 г., представляет интерес анализ более обширного материала, охватывающего средние широты всего Северного полушария, с включением в анализ пространственно-временных вариаций отношений СО2/CO, СО2/N2O, СН4/CO и СН4/N2O.

3.5. Сезонные вариации в толще атмосферы и вблизи земной поверхности

Средние сезонные разности концентраций в приземном слое и в толще атмосферы за период с 2015 по 2021 гг. на ст. Обнинск приведены на рис. 9. Наибольшие превышения приземных концентраций gСН4 и gСО2 над концентрациями в толще атмосферы наблюдаются в зимний период. В летний период для углекислого газа концентрации примерно совпадают, а для метана характерной чертой сезонного хода является постоянное превышение приземных концентраций над средними концентрациями в толще атмосферы. Ранее, при анализе данных ст. Обнинск, для определения фоновой составляющей приземных концентраций использовался метод фильтрации данных с выделением минимальных внутридневных величин [Арефьев и др., 2014; Арефьев и др., 2015]. Аналогичный подход был использован в [Макарова и др., 2015] при анализе вариаций метана на станции Петергоф (СПбГУ). На рис. 9 квадратами показаны сезонные разности концентраций, полученные для минимальных величин. Обращает на себя внимание почти параллельный сдвиг разности gСН4 – хСН4 и gСО2 – хСО2 при переходе от средних к минимальным приземным концентрациям. Для метана минимальные приземные и средние по высоте концентрации в мае-августе совпадают в пределах 30 ppb, а в остальные месяцы превышение приземных концентраций сохраняется и составляет 70–150 ppb. Качественно это может означать, что для ст. Обнинск в холодное время года возрастает вклад местных источников метана, а в летнее время местные источники метана незначительны и регистрируются концентрации метана, определяемые средним переносом. Учитывая, что для ст. Обнинск основные направления ветров западные и юго-западные (до 40%) и южные (до 15%) [https://world-weather.ru], то воздушные массы в течение года в основном приходят из центральной и южной Европы.

Рис. 9.

Сезонный ход разности среднемесячных значений gСН4 – хСН4 (а) и gСО2 – хСО2 (б). 1 – усредненные с 2015 по 2021 гг. среднемесячные значения, 2 – среднемесячные значения для минимальных величин.

В случае углекислого газа сезонный ход средних и минимальных значений разности gСО2 – хСО2 качественно такой же, как и для метана, однако, количественные показатели различаются. Превышение минимальных приземных концентраций над концентрациями в толще атмосферы (6–15 ppm) сохраняется только для декабря, января и февраля, а с мая по август хСО2 больше gСО2 на 20–30 ppm. Этот результат согласуется с анализом сезонных вариаций СО2 в работе [Olsen, Randerson, 2004], где было показано, что для средних широт Северного полушария характерно падение амплитуды сезонных колебаний СО2 с высотой, причем летом средние по высоте концентрации превышают приземные, а зимой и весной, в период роста эмиссии за счет сжигания топлива, хСО2 ниже приземных концентраций.

4. ЗАКЛЮЧЕНИЕ

Результаты измерений в приземном слое в период 1998–2021 гг. показали значительный рост концентраций углекислого газа и метана. Небольшой отрицательный тренд СН4, наблюдавшийся в 1998–2007 гг., в дальнейшем сменился положительным трендом. За период наблюдений концентрация метана выросла на 124 ppb, а углекислого газа на 62 ppm. Скорость роста концентрации СО2 на ст. Обнинск в процентном отношении (0.65% в год) в полтора раза превосходит наблюдавшийся с 2007 г. рост СН4 (0.42% в год). Тренды, амплитуда колебаний и сезонный ход приземной концентрации метана и углекислого газа на ст. Обнинск хорошо согласуются с данными ст. Хегихатсал, Венгрия.

Сопоставление средних по высоте концентраций метана хСН4 и углекислого газа xСО2 со спутниковыми данными GOSAT показало, что для наземных наблюдений в холодное время года наблюдается небольшое систематическое превышение. Для привязки к данным GOSAT были введены поправочные коэффициенты, равные 0.97 для метана и 0.98 для углекислого газа. Тренды и средние значения по данным ст. Обнинск за сопоставимый период с 2016 по 2020 гг. достаточно хорошо согласуются с наземными наблюдениями на европейских станциях Бремен, Карлсруэ, Орлеан и Париж, при этом для ст. Обнинск несколько выше амплитуда годовых вариаций. Коэффициенты корреляции между измерениями на ст. Обнинск и на станциях TCCON превышают 0.84 и 0.73 для метана и углекислого газа соответственно.

Практически одновременные измерения прибором МР-32/МГС концентраций в приземном слое (gСН4 и gСО2) и в толще атмосферы (xСН4 и xСО2) позволили провести сравнительный анализ сезонной и долговременной изменчивости и оценить взаимосвязи между приземными и средними в толще концентрациями метана и углекислого газа. По результатам кросс-корреляционного вейвлетного анализа фаза годовых вариаций средних по высоте концентраций xCH4 опережает фазу годовых вариаций приземных значений gCH4 на 2–3 месяца, а вариации xСО2 отстают по фазе от годовых вариаций gСО2 на 1–2 месяца. Годовые вариации gСН4 и gСО2 проходят синхронно, а вариации хСН4 и хСО2 примерно в противофазе.

Для сезонного хода среднедневных приземных концентраций метана и углекислого газа характерно их превышение над средними концентрациями в толще атмосферы. В то же время, минимальные приземные концентрации углекислого газа в мае-августе ниже средних концентраций в толще атмосферы на 20–30 ppm, а превышение (6–15 ppm) сохраняется только для декабря, января и февраля. В период с мая по август минимальные приземные концентрации метана совпадают со средними концентрациями в толще атмосферы, в остальные месяцы их превышение составляет 70–150 ppb. Качественно это может означать, что для ст. Обнинск в летнее время местные источники метана и углекислого газа незначительны и регистрируются концентрации, определяемые адвективным переносом при преобладании западных и юго-западных ветров. В холодное время года, вероятно, возрастает вклад местных источников за счет роста эмиссии при сжигании топлива.

Список литературы

  1. Арефьев В.Н., Дианов–Клоков В.И., Малков И.П. Полевой спектральный комплекс для исследований содержания в атмосфере загрязняющих газов // Труды ИЭМ. М.: Гидрометеоиздат, 1978. Вып. 8(81). С. 73–78.

  2. Арефьев В.Н., Вишератин К.Н., Кашин Ф.В., Устинов В.П. Аппаратура для измерений спектроскопическим методом интегрального содержания газов в атмосфере // Труды ИЭМ. 1995. Вып. 25(160). С. 119–125.

  3. Арефьев В.Н., Баранов Ю.И., Баранова Е.Л., Бугрим Г.И., Каменоградский Н.Е., Кашин Ф.В. Изменчивость содержания метана в приземном слое и в толще атмосферы // Изв. РАН. Физика атмосферы и океана. 2004. Т. 40. № 3. С. 403–409.

  4. Арефьев В.Н., Кашин Ф.В., Акименко Р.М., Баранов Ю.И., Бугрим Г.И., Вишератин К.Н., Кальсин А.В., Каменоградский Н.Е., Сизов Н.И., Устинов В.П., Упэнэк Л.Б. Исследования в области атмосферной спектроскопии // “Проблемы гидрологии и мониторинга окружающей среды”. Сб. науч. тр., Обнинск, 2010. С. 85–104.

  5. Арефьев В.Н., Каменоградский Н.Е., Кашин Ф.В., Шилкин А.В. Фоновая составляющая концентрации двуокиси углерода в приземном воздухе. (Станция мониторинга “Обнинск”) // Изв. РАН. Физика атмосферы и океана. 2014. Т. 50. № 6. С. 655–662.

  6. Арефьев В.Н., Акименко Р.М., Кашин Ф.В., Упэнэк Л.Б. Фоновая составляющая концентрации метана в приземном воздухе. (Станция мониторинга “Обнинск”) // Изв. РАН. Физика атмосферы и океана. 2015. Т. 51. № 6. С. 1–9.

  7. Дианов–Клоков В.И., Юрганов Л.Н. Измерения интегрального содержания примесей СО, СН4 и N2О в атмосфере // Изв. АН СССР. ФАО. 1982. Т. 18. № 7. С. 1159–1167.

  8. Кабанов Д.М., Сакерин С.М. Результаты исследований общего влагосодержания атмосферы методом оптической гигрометрии. Ч. 1. Aнализ методики и результатов калибровки // Оптика атмосферы и океана. 1995. Т. 8. № 06. С. 852–860.

  9. Кондратьев К.Я., Крапивин В.Ф. Моделирование глобального круговорота углерода. М.: Физматлит, 2004. 336 с.

  10. Макарова М.В., Кирнер О., Тимофеев Ю.М., Поберовский А.В., Имхасин Х.Х., Осипов С.И., Макаров Б.К. Годовой ход и долговременный тренд атмосферного метана в районе СанктПетербурга // Изв. РАН. Физика атмосферы и океана. 2015. Т. 51. № 4. С. 493–501.

  11. Мироненков А.В., Поберовский А.В., Тимофеев Ю.М. Методика интерпретации инфракрасных спектров прямой солнечной радиации для определения общего содержания атмосферных газов // Изв. РАН, Физика атмосферы и океана. 1996. Т. 32. № 2. С. 207–215.

  12. Ракитин В.С., Еланский Н.Ф., Панкратова Н.В., Скороход А.И., Джола А.В., Штабкин Ю.А., Ван Пусай, Ван Ген Чен, Васильева А.В., Макарова М.В., Гречко Е.И. Исследование трендов общего содержания CO и CH4 над Евразией на основе анализа наземных и орбитальных спектроскопических измерений // Оптика атмосферы и океана. 2017. Т. 30. № 06. С. 449–456. https://doi.org/10.15372/AOO20170601

  13. Тимофеев Ю.М. Исследование атмосферы Земли методом прозрачности. СПб.: Наука, 2016. 367 с.

  14. Шашков А.А., Броунштейн А.М., Жуков А.В., Привалов В.И., Фабер Е.В. Спектроскопические методики определения интегрального содержания СО, N2О и СН4 в вертикальном столбе атмосферы // Труды ГГО им. А.И. Воейкова. 1985. Вып. 496. С. 7–23.

  15. Baranov Yu.I., Baranova E.L., Bougrim G.I., Kashin F.V. Temporal variability of methane, carbon oxide and dioxide and dinitrogen oxide in surface air // Proc. SPIE. 2000. V. 4341. P. 438–442.

  16. Buschmann M., Deutscher N.M., Sherlock V., Palm M., Warneke T., Notholt J. Retrieval of xCO2 from ground-based mid-infrared (NDACC) solar absorption spectra and comparison to TCCON // Atmos. Meas. Tech. 2016. V. 9. P. 577–585. https://doi.org/10.5194/amt-9-577-2016

  17. Chen H., Winderlich J., Gerbig C., Hoefer A., Rella C.W., Crosson E.R., Van Pelt A.D., Steinbach J., Kolle O., Beck V., Daub B.C., Gottlieb E.W., Chow V.Y., Santoni G.W., Wofsy S.C. High-accuracy continuous airborne measurements of greenhouse gases (CO2 and CH4) using the cavity ring-down spectroscopy (CRDS) technique // Atmos. Meas. Tech. 2010. V. 3. P. 375–386.

  18. Dlugokencky E.J., Nisbet E.G., Fisher R., Lowry D. Global atmospheric methane: Budget, changes and dangers // Philos. Trans. R. Soc. A. 2011. V. 369. P. 2058–2072. https://doi.org/10.1098/rsta.2010.0341

  19. Dlugokencky E.J., Mund J.W., Crotwell A.M., Crotwell M.J., Thoning K.W. Atmospheric Carbon Dioxide Dry Air Mole Fractions from the NOAA GML Carbon Cycle Cooperative Global Air Sampling Network, 1968–2020, Version: 2021-07-30. https://doi.org/10.15138/wkgj-f215

  20. ESRL, 2021. Earth System Research Laboratory, Global Monitoring Division, Carbon Cycle Greenhouse Gases, http://www.esrl.noaa.gov/.

  21. Frey M., Sha M.K., Hase F., Kiel M., Blumenstock T., Harig R., Surawicz G., Deutscher N.M., Shiomi K., Franklin J.E., Bösch H., Chen J., Grutter M., Ohyama H., Sun Y., Butz, A., Mengistu Tsidu G., Ene D., Wunch D., Cao Z., Garcia O., Ramonet M., Vogel F., Orphal J. Building the COllaborative Carbon Column Observing Network (COCCON): long-term stability and ensemble performance of the EM27/SUN Fourier transform spectrometer // Atmos. Meas. Tech. 2019. V. 12. P. 1513–1530. https://doi.org/10.5194/amt-12-1513-2019

  22. Gavrilov N.M., Makarova M.V., Poberovskii A.V., Timofeyev Yu.M. Comparisons of CH4 ground-based FTIR measurements near Saint Petersburg with GOSA-T observations // Atmos. Meas. Tech. 2014. V. 7. P. 1003–1010.

  23. Ghasemifard H., Vogel F.R., Yuan Y., Luepke M., Chen J., Ries L., Leuchner M., Schunk C., Noreen Vardag S., Menzel A. Pollution Events at the High-Altitude Mountain Site Zugspitze-Schneefernerhaus (2670 m a.s.l.), Germany // Atmosphere. 2019. V. 10(6). P. 330. https://doi.org/10.3390/atmos10060330

  24. GOSAT, 2021. Greenhouse Gases Observing Satellite, https://data2.gosat.nies.go.jp; https://www.gosat.nies. go.jp/en/about_2_observe.html

  25. Grinsted A., Moore J.C., Jevrejeva S. Application of the cross wavelet transform and wavelet coherence to geophysical time series // Nonlin. Processes Geophys. 2004. № 11. P. 561–566. https://doi.org/10.5194/npg-11-561-2004

  26. Hase F., Blumenstock T., Dohe S., Groß J., Kiel M. TCCON data from Karlsruhe, Germany, Release GGG2014R1. 2017. https://doi.org/10.14291/tccon.ggg2020.karlsruhe01.R0

  27. Haszpra L., Ferenczi Z., Barcza Z. Estimation of greenhouse gas emission factors based on observed covariance of CO2, CH4, N2O and CO mole fractions // Environ Sci. Eur. 2019. P. 31–95. https://doi.org/10.1186/s12302-019-0277-y

  28. Howarth R.W. Ideas and perspectives: is shale gas a major driver of recent increase in global atmospheric methane? // Biogeosciences. 2019. V. 16. P. 3033–3046. https://doi.org/10.5194/bg-16-3033-2019

  29. HYSPLIT, 2021. The Hybrid Single-Particle Lagrangian Integrated Trajectory model, https://www.ready.noaa. gov/hypub-bin/trajtype.pl.

  30. Infraspek. 2021. http://www.infraspek.ru/produktsiya/ spektrometryi/fsm-2203-2/; http://www.infraspek.ru/ produktsiya/kyuvetyi/kyuvetyi-gazovyie/.

  31. IPCC, 2021. Global Carbon and other Biogeochemical Cycles and Feedbacks. In: Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge University Press. 2021.

  32. Kashin F. Variations of CO2 Mixing Ratios in the Air near the Ground in the European Territory of Russia // J. Environmental Science and Engineering A. 2013. V. 2. № 9. P. 531–536.

  33. Kivimäki E., Lindqvist H., Hakkarainen J., Laine M., Sussmann R., Tsuruta A., Detmers R., Deutscher N.M., Dlugokencky E.J., Hase F., Hasekamp O., Kivi R., Morino I., Notholt J., Pollard D.F., Roehl C., Schneider M., Sha M.K., Velazco V.A., Warneke T., Wunch D., Yoshida Y., Tamminen J. Evaluation and Analysis of the Seasonal Cycle and Variability of the Trend from GOSAT Methane Retrievals // Remote Sens. 2019. V. 11(7). P. 882. https://doi.org/10.3390/rs11070882

  34. Kubistin D., Plaß-Dülmer C., Arnold S., Lindauer M., Müller–Williams J. // ICOS RI, 2021. ICOS ATC CO Release, Karlsruhe (30.0 m), 2019-01-31–2021-01-31. https://hdl.handle.net/11676/7l7HG_GiOgS_u3TEuXB6HkFO

  35. NOAA, 2021. https://www.noaa.gov/news-release/i-ncrease-in-atmospheric-methane-set-another-record-during-2021.

  36. Olsen S.C., Randerson J.T. Differences between surface and column atmospheric CO2 and implications for carbon cycle research // J. Geophys. Res. 2004. V. 109. D02301. https://doi.org/10.1029/2003JD003968

  37. Oshchepkov S., Bril A., Yokota T., Yoshida Y., Blumenstock T., Deutscher N.M., Dohe S., Macatangay R., Morino I., Notholt J., Rettinger M., Petri C., Schneider M., Sussman R., Uchino O., Velazco V., Wunch D., Belikov D. Simultaneous retrieval of atmospheric CO2 and light path modification from space-based spectroscopic observations of greenhouse gases: Methodology and application to GOSAT measurements over TCCON sites // Appl. Opt. 2013. V. 52. P. 1339–1350.

  38. Ostler A., Sussmann R., Rettinger M., Deutscher N.M., Dohe S., Hase F., Jones N., Palm M., Sinnhuber B.-M. Multistation intercomparison of column–averaged methane from NDACC and TCCON: impact of dynamical variability // Atmos. Meas. Tech. 2014. V. 7. P. 4081–4101.

  39. Rakitin V.S., Skorokhod A.I., Pankratova N.V., Shtabkin Yu.A., Rakitina A.V., Wang G., Vasilieva A.V., Makarova M.V., Wang P. Recent changes of atmospheric composition in background and urban Eurasian regions in XXI-th century // IOP Conf. Series: Earth and Environmental Science. 2020. T. 606. № 1. P. 012 048. https://doi.org/10.1088/1755-1315/606/1/012048

  40. Rothman L.S., Gordon I.E., Babikov I.E., Barbe A., Benner C.D. et al. The HITRAN 2012 molecular spectroscopic database // J. Quant. Spectrosc. Radiat. Transfer. 2013. V. 130. P. 4–50.

  41. SFIT4. Version V0.9.4.4. The University Corporation for Atmospheric Research, https://wiki.ucar.edu/di-splay/ sfit4.

  42. Takeda M., Nakajima H., Murata I., Nagahama T., Morino I., Toon G.C., Weiss R.F., Mühle J., Krummel P.B., Fraser P.J., Wang H.-J. First ground-based Fourier transform infrared (FTIR) spectrometer observations of HFC-23 at Rikubetsu, Japan, and Syowa Station, Antarctica // Atmos. Meas. Tech. 2021. V. 14. P. 5955–5976.https://doi.org/10.5194/amt-14-5955-2021

  43. TCCON, 2021. https://tccondata.org/; https://tccon-wiki. caltech.edu/Sites.

  44. Tohjima Y., Kubo M., Minejima C., Mukai H., Tanimoto H., Ganshin A., Maksyutov S., Katsumata K., Machida T., Kita K. Temporal changes in the emissions of CH4 and CO from China estimated from CH4/CO2 and CO/CO2 correlations observed at Hateruma Island // Atmos. Chem. Phys. 2014. V. 14. P. 1663–1677.

  45. Visheratin K.N., Baranova E.L., Bugrim G.I., Ivanov V.N., Krasnopeeva E.I., Sakhibgareev D.G., Ustinov V.P., Shilkin A.V., Baranov Yu.I., Kashin F.V. MR-32/MGS setup for monitoring gas composition of atmosphere // XXVII International Symposium “Atmospheric and Ocean Optics. Atmospheric Physics”. 2021 july 05–09. Moscow. 2021. https://symp.iao.ru/files/symp/aoo/ 27/presentation_13663.pdf.

  46. WACCAM, 2013. Whole Atmosphere Community Climate Model, https://www2.acom.ucar. edu/gcm/waccm; ftp://acd.ucar.edu/user/jamesw/IRWG/2013/.

  47. WDCGG, 2021 World Data Centre for Greenhouse Gases, https://gaw.kishou.go.jp.

  48. WMO, 2021. The State of Greenhouse Gases in the Atmosphere Based on Global Observations through 2020 // Greenhouse Gas Bulletin. No. 17, 25 October 2021.

  49. Wong K.W., Fu D., Pongetti T.J., Newman S., Kort E.A., Duren R., Hsu Y.-K., Miller C.E., Yung Y.L., Sander S.P. Mapping CH4: CO2 ratios in Los Angeles with CLARS-FTS from Mount Wilson, California // Atmos. Chem. Phys. 2015. V. 15. P. 241–252. https://doi.org/10.5194/acp-15-241-2015

  50. Yoshida Y., Ota Y., Eguchi N., Kikuchi N., Nobuta K., Tran H., Morino I., Yokota T. Retrieval algorithm for CO2 and CH4 column abundances from short-wavelength infrared spectral observations by the Greenhouse gases observing satellite // Atmos. Meas. Tech. 2011. V. 4. P.717–734.https://doi.org/10.5194/amt-4-717-2011

Дополнительные материалы отсутствуют.