CC BY
60
Вычислительные технологии
Том 12, № 5, 2007
ПРОГНОЗИРОВАНИЕ РАСПРОСТРАНЕНИЯ ДЫМОВЫХ АЭРОЗОЛЕЙ НА ТЕРРИТОРИИ
СИБИРИ
О. А. ДУБРОВСКАЯ, Е. Г. КЛИМОВА Институт вычислительных технологий СО РАН, Новосибирск, Россия e-mail: olga@ict.nsc.ru, klimova@ict.nsc.ru
Gas components emitted at fires may propagate far away from the source of fire due to the atmospheric circulation. Therefore one of the important problems is a description of the spatial and temporal distribution of the smoke aerosol using actual meteorological information and the data on fires. For an estimation of the territory covered by the fire smoke, a technique for restoration of the values of concentration of the gas components is proposed which utilizes both meteorological and fire data. Pollutant concentrations in the given region were estimated using the reconstruction of the backward trajectories. Data assimilation technique was used in the restoration process. Numerical experiments on estimation of the gas component concentrations were compared to the actual trajectories of smoke plumes obtained by NOAA satellites.
Введение
Газовые составляющие атмосферы, такие как СО, СО2, СН4, называемые парниковыми газами, существенно влияют на характер атмосферных процессов и экологическую обстановку в исследуемом регионе [1]. Одним из источников этих газов являются лесные пожары в бореальных лесах на территории Западной и Восточной Сибири, Якутии, Дальнего Востока.
Так как в Сибири сосредоточены большие запасы наземного углерода, то увеличение площадей пожарищ, длительности пожарного сезона и силы пожаров приводит к тому, что освобождаются значительные объемы углерода [2, 3]. В годы с экстремальной пожарной активностью общая эмиссия углерода может быть на 37-41 % больше, чем в годы с нормальной пожарной активностью, — по причине увеличения сгорания органических веществ в почве. Средние оценки для стандартного сценария динамики выбросов рассматриваемых парниковых газов составляют: С02 (555-1031 тераграмм (Тг)), СО (43-80), СН4 (2.4-4.5 Тг). Эти оценки представляют 10, 15, 19% соответственно от глобальных оценок, сделанных по всему Земному шару для пожаров растительности [2].
Газовые составляющие, выделяющиеся при пожарах, распространяются далеко за их пределы благодаря атмосферной циркуляции. Поэтому одна из важных задач — описать распространение дымового аэрозоля в пространстве и во времени, привлекая для этого фактическую метеорологическую информацию и данные о пожарах.
© Институт вычислительных технологий Сибирского отделения Российской академии наук, 2007.
В настоящей работе описывается методика восстановления значений концентрации газовых составляющих на примере СО2 в узлах сетки заданного региона по данным о лесных пожарах и по метеорологическим материалам международного архива "Ре-анализ". Концентрации СО2 в заданном регионе оцениваются с помощью построения обратных траекторий [4]. При этом данные о ветре могут быть заданы на более крупной сетке, чем сетка, на которой вычисляется концентрация. Для восстановления значений концентрации вдоль траекторий применялась процедура усвоения данных [5]. В работе использовалась информация о состоянии атмосферы с 10 по 20 августа 2002 года, что соответствует пику пожарной активности на территории Республики Саха (Якутия), который был зафиксирован спутником МОЛЛ [6]. На рис. 1 представлена карта-схема выжженной территории за пожарный сезон с 1 марта по 1 ноября 2002 года. На ней представлена карта Сибири и Дальнего Востока. Черным цветом изображены территории, которые были уничтожены огнем в результате пожаров. Видно, что массовые лесные пожары наибольшей мощности прошли на территории Якутии. На рис. 2 показана динамика пожарной активности за этот же пожарный сезон на азиатской части России.
Рис. 1. Картосхема распределения пожаров за пожарный сезон 2002 года
Рис. 2. Динамика пожарной активности 2002 года (с 1 марта по 1 ноября) на азиатской части России; шкала слева показывает число пожаров, шкала справа — выгоревшую площадь
1. Оценка концентрации СО2 по данным о сгоревшей биомассе
Использование спутниковых данных о выгоревшей территории позволяет сделать более точные оценки эмиссий углерода от лесных пожаров. Количество горючего материала, в действительности сгоревшего в пожаре, будет зависеть от запаса и состояния топлива, его типа, климатических и метеорологических факторов, а также интенсивности пожара. Модели расхода горючих материалов в общем случае включают эмпирические коэффициенты для прогнозирования относительных объемов тления и горения при различных условиях горения.
Подходы к оценке общих эмиссий углерода (С4) и эмиссий конкретных шлейфовых газов (Е3) основаны на известных моделях, предложенных в работе [7]:
С = А • В • /о • в, (1)
Es = С • в/3, (2)
где А — площадь пожарищ, га; В — средний запас биомассы, т/га; /о — доля углерода в биомассе; в — масштабирующий коэффициент для доли биомассы, которая сгорела (связан с типом и интенсивностью пожара, а также с видом лесного горючего материала); е/ — коэффициент эмиссии для данного вида шлейфового газа в (обычно выражается как доля данного вещества, сгоревшего во время пожара, по отношению к сухому топливу).
В настоящей работе расчеты оценки эмиссии СО2 по формулам (1) и (2) были сделаны на примере якутских пожаров 2002 года с использованием данных о площади, пройденной огнем, со следующими эмпирическими коэффициентами: средний запас горючего материала рассматриваемого региона — 8.0 т углерода на 1 га, доля углерода / = 1, масштабирующий коэффициент в = 1, коэффициент эмиссии для СО2 — 0.45 [8].
2. Вычисление распространения аэрозоля с помощью обратных траекторий
Для оценки территории, охваченной дымовым аэрозолем, разработана методика, основанная на вычислении обратных траекторий [4]. В расчетах используются данные архива "Реанализ" Европейского центра среднесрочных прогнозов погоды (ECMWF).
Для расчета траекторий привлекаются данные о ветре, температуре и высоте на 23 изобарических поверхностях с шагом по времени 6 ч.
Уравнения траекторий записываются в виде системы обыкновенных дифференциальных уравнений:
dA 1 , . , ,
— =-u[t, A, ^p,p), (3)
dt a cos t
= av(t, A, p,p), (4)
dt a
= w{t,A,t,p), (5)
где A(t), t(t),p(t) — текущие координаты перемещающейся частицы (широта, долгота, давление); u, v — горизонтальные составляющие скорости ветра в сферической системе координат; w — вертикальная скорость в p-системе координат; a — радиус Земли; w рассчитывается из уравнения притока тепла [9]:
и \ \ gv
RT (la - Y)
dT u dT v dT dt a cos t dA a dt
Система уравнений (3)-(5) решается на временном интервале т. е. в обрат-
ном времени, поэтому в качестве начальных данных задаются координаты частицы N NN 4-М
(А1 ,рМ) в момент времени г .
Для нахождения решения системы (3)-(5) используется вычислительная схема, предложенная в работе [9]. На каждом временном шаге находится значение координаты р из уравнения (5) с использованием простейшей схемы односторонних разностей:
Рп+1 = Рп + 'ш(А,(р,р,Ь)пДЬ, (6)
где ДГг — шаг по времени.
Затем для нахождения горизонтальных координат решается система уравнений (3)-(4) методом Рунге—Кутты четвертого порядка точности [10]. Так как достоверность данных о поле относительно мала по сравнению с полями и и V, использование схемы, более точной по сравнению со схемой односторонних разностей, для нахождения рп+1 не представляется оправданным.
3. Фоновое значение СО2
Естественные источники образования углекислого газа — растения, океан, вулканические выбросы, лесные пожары. В первом приближении можно считать СО2 в атмосфере равномерно перемешанным, однако для умеренных широт был выявлен максимум
содержания СО2 (340-350 млн-1) на высотах 25-30 км. В начале эпохи бурного индустриального развития фоновое содержание у поверхности земли составляло всего лишь 280 млн-1. В настоящее время оно равно 334-356 млн-1 [3]. Скорость роста концентрации СО2 в год достигает 0.5%, или 0.7 млн-1. Данный газ имеет самый большой относительный вклад в антропогенный парниковый эффект (60 %). Мы рассматриваем поступление этой примеси в атмосферу в результате пиролиза древесины и подстилающей поверхности при горении бореальных лесов.
4. Усвоение данных при восстановлении концентрации дымовых аэрозолей
В работе рассматривается задача о восстановлении значений объемной концентрации примеси в дискретных пространственных узлах (Лi, <fi,p), расположенных на выбранной изобарической поверхности p = const, где (Лг,рг) — узлы пространственной прямоугольной широтно-долготной сетки.
Под усвоением данных принято понимать совместный учет математической модели и данных наблюдений для описания пространственно-временной картины распределения исследуемых полей. Системы усвоения данных разрабатываются в настоящее время в различных областях: метеорологии, океанологии, охране окружающей среды [11, 12]. В настоящей работе в качестве данных наблюдений берутся значения эмиссии СО2 в зонах пожара, вычисленные по изложенной выше, в п. 1, методике. Методика усвоения данных частично заимствована из работы [5].
Усвоение данных проводится на основе цикла прогноз—анализ, при этом под анализом данных понимается интерполяция наблюдений в заданную точку. Наиболее популярным методом анализа является в настоящее время метод оптимальной интерполяции [13], однако для его реализации требуется точное знание среднеквадратических ошибок прогноза и наблюдений. Такой информацией о данных, получаемых при наблюдениях за лесными пожарами, мы не располагаем. По этой причине анализ осуществляется с помощью метода последовательных коррекций [14].
Предполагается, что значение этой концентрации в случае отсутствия данных наблюдений остается неизменным вдоль траектории частицы [4]. Для i-го сеточного узла широтно-долготной сетки производится вычисление координат траектории частицы (Л(£), p(t),p(t)) в обратном времени. Процедура предназначена для определения значения концентрации qi в i-м сеточном узле в момент времени tN. Обозначим через qn (n = 0,..., N) значение объемной концентрации углекислого газа в i-м сеточном узле в момент времени tn. В начальный момент времени tN зададим значение концентрации qN = g. Шаг "прогноза" вдоль траектории для одного временного шага At можно формально записать следующим образом:
qn = qn+1- (7)
В настоящей работе мы используем следующие индексные обозначения: f — шаг прогноза, a — шаг анализа.
Пусть данные о сгоревшей биомассе имеются в K точках земной поверхности с координатами {(Лк(Р), <fk(tj)), k = 1,...,K} в моменты времени t1,...,tJ. Процедура усвоения данных производится следующим образом: от момента времени tj+1 до tj
считается "прогноз" по формуле (7). В момент Ь7, в котором присутствуют данные наблюдений, производится "анализ", т. е. оценка значений концентрации в г-м сеточном узле по вычисленному на этот момент значению концентрации и данным наблюдений. Полученное значение "анализа" является начальным для вычисления следующего шага "прогноза" (7).
Значение концентрации qla, соответствующей точке траектории (А7 , р7) в момент времени Ь7, в котором есть информация о пожарах, определяется по значениям данных наблюдений, известным в окрестности этой точки радиусом К. В момент времени Ь производится уточнение значения концентрации с помощью метода последовательных коррекций [14]:
К1
Е ajkAqd
= ^ + -' (8)
Е — к= 1
где К1 — количество точек наблюдений, попавших в окрестность точки (А7, ^ ,р7) ради-
_0 5 _ , , ,
усом К; а^к = е ' (вт)2; г7к — расстояние между точкой траектории (А(Ь),^(Ь3),р(Ь)) и к-й точкой наблюдений; Вт — радиус корреляции; Д^ = ДЬ ■ Пк. Величина пк определяется по данным наблюдений о выделении газовых составляющих в единицу времени. Для каждой точки наблюдений величина Пк сохраняет свое постоянное значение в течение 7-10 сут. Значение Д^ считается постоянным до высоты 1.5 км. Для высот, больших 1.5 км, изменение Дд{ по вертикали вычисляется по формуле
rjk qm\2
'J
__„лт,,,«,,,,,,,,,,. Л J
■ ■ _ (Н1-Н2)2
д^(р0 = Д^(Р2)е 2Ь2 ,
где р1, р2 — давление; Н1, к2 — высоты на уровнях р1, р2; Ь — характерный масштаб по вертикали [15]. В проведенных экспериментах бралось Вт = 100 км, т = 0, К = 1000 км, Ь = 1.5 км.
5. Численные эксперименты
Цель численных экспериментов состояла в том, чтобы получить пространственное распределение концентрации СО2 и сравнить со спутниковыми изображениями, на которых можно проследить динамику переноса дымовых шлейфов. При сравнении результатов расчетов с реальными данными была использована программа приема изображений облачности с космических аппаратов NOAA в формате Analog Picture Transmission (APT) с последующей тематической обработкой.
Восстановление концентрации проводилось для региона 63.5° — 80° с.ш. и 70° — 172.5° в.д. в узлах сетки 2.5° х 2.5°. Расчеты соответствуют интервалу времени с 10 по 20 августа 2002 года. В численных экспериментах данные о ветре и температуре линейно интерполировались по времени на шести уровнях. Эти уровни соответствуют высотам с давлением 1000, 925, 850, 775, 700, 600 мбар. Фоновая концентрация СО2 задавалась равной 340 млн-1 на этих уровнях. Информация о пожарах включает в себя дату (число, месяц, год), время (час, минута, секунда), пространственные координаты пожара, площадь горения. При расчете эмиссии использовались данные о сгоревшей биомассе при массовых лесных пожарах на территории Республики Саха (Якутия) [6].
Длительность пожарного сезона в численных экспериментах бралась 90 дней с шагом по времени 1 ч. Выгоревшая площадь, зафиксированная по данным спутника МОЛЛ в 2002 году при пожарах в Якутии, составила около 5 млн га, а восстановленные значения газовых примесей по данным о сгоревшей биомассе — СО2 = 28.6 • 109 кг, СО= 2.89 • 109, СН4 = 0.286 • 109 кг.
На рис. 3 представлен спутниковый снимок в АРТ-формате, на котором замкнутой линией очерчена зона задымления и пожаров. На ней отмечены контрольные точки, координаты которых приведены в таблице. На рис. 4 представлены результаты одного из численных экспериментов. Результаты расчетов (рис. 4) и космический снимок (см. рис. 3) соответствуют дате 14 августа 2002 года. Все пожары на рис. 4 отмечены
Координаты контрольных точек на рис. 3
N С.ш. В.д.
1 60 127
2 63 120
3 68 121
4 69 129
5 63 132
6 61 129
Рис. 3. Снимок пожаров и зоны задымления в Якутии 14.08.2002
в виде черных точек. Изолиниями очерчены зоны распространения эмиссии двуокиси углерода и представлено ее количественное распределение. Более детальный анализ показывает, что территория с самым высоким уровнем концентрации на рис. 4 имеет схожие очертания с зоной пожаров и дымового шлейфа на спутниковом изображении. На рис. 5 представлены результаты численного эксперимента для района 63.5° — 73.5° с.ш. и 90° — 152.5° в.д. Процедура позволяет оценивать концентрацию примеси в различ-
Рис. 4. Результаты расчетов распространения примеси и ее концентрации для региона 63.5° 80° с.ш. и 70° - 172.5° в.д.
Рис. 5. Результаты расчетов распространения примеси и ее концентрации для региона 63.5° 73.5° с.ш. и 90° - 152.5° в.д.
ных заданных регионах. На рис. 4 и 5 представлено распределение концентрации для различных областей. Из этих рисунков видно, что изолинии в области локализации пожаров совпадают.
Заключение
Таким образом, представленная методика восстановления распространения объемной концентрации примеси при использовании реальных метеорологических данных и данных о лесных пожарах позволяет описывать распространение газовых составляющих при лесных пожарах за счет атмосферной циркуляции. Результаты численных экспериментов, полученные при использовании данной методики, хорошо согласуются с фактическим распределением дымовых шлейфов по данным спутников серии NOAA.
Выражаем благодарность за предоставленные данные о пожарах Отделу приема и обработки космической информации ВНИИ ГОЧС г. Красноярска.
Список литературы
[1] Парниковый эффект, изменение климата и экосистемы / Под редакцией Б. Болина, Б.Р. Дееса, Дж. Ягера, Р. Уоррика. Л.: Гидрометеоиздат, 1989. С. 33-234.
[2] SojA A.J., Cofer W.R., Shugart H.H. et al. Estimating fire emissions and disparities in boreal Siberia (1998 through 2002) // J. of Geophys. Res. 2004. Vol. 109, D14S06, doi:10.1029/2004JD004570.
[3] МицЕль А.А., Фирсов К.М., Фомин Б.А. Перенос оптического излучения в молекулярной атмосфере. Томск, 2001. C. 67-80.
[4] Sutton R.T., O'Neill A., Taylor F.M. High-resolution stratospheric tracer fields estimated from satellite observations using lagrangian trajectory calculations //J. Atmos. Sci. 1994. Vol. 51. P. 2995-3005.
[5] Khattatov B.V., Gill J.C., Lyjak L.V. et al. Assimilation of photo chemically active species and a case analysis of UARS data // J. of Geophys. Res. 1999. Vol. 104. P. 18715-18737.
[6] Sukhinin A.I., Ivanov V.V., Ponomarev E.I. et al. The 2002 Fire season in the Asian part of the Russian Federation: A view from space (2003) // Intern. Forest Fire News. 2003. N 28. Р. 18-28.
[7] Seiler W., Crutzen P.J. Estimates of gross and net fluxes of carbon between the biosphere and atmosphere // Clim. Change. 1980. Vol. 2. P. 207-247.
[8] Kasiscke E.S., Hyer E.J., Novelli P.C. et al. Influences of boreal fire emissions on Northern Hemisphere atmospheric carbon and carbon monoxide // Global Biogeochem. Cycles. Vol. 19, GB1012, doi:10.1029/2004GB002300.
[9] Белоусов С.Л., Юсупов Ю.И. Расчет трехмерных траекторий воздушных частиц // Метеорология и гидрология. 1991. № 12. С. 41-48.
[10] КАлиткин Н.Н. Численные методы М.: Наука, 1978. С. 237-260.
[11] Климова Е.Г., КилАновА Н.В. Численные эксперименты по оценке эмиссии метана на основе системы усвоения данных о пассивной примеси в атмосфере Северного полушария // Оптика атмосферы и океана. 2006. Т. 11. С. 961-964.
[12] КилАновА Н.В., Климова Е.Г. Численные эксперименты по оценке систематической ошибки модели в задаче усвоения данных о концентрации пассивной примеси // Вычисл. технологии. 2006. Т. 11, № 5. С. 32-40.
[13] ГАндин Л.С. Объективный анализ метеорологических полей Л.: Гидрометеоиздат, 1963. 287 с.
[14] Bratseth A.M. Statistical interpolation by means of successive correction // Tellus. 1986. N 38A. P. 439-447.
[15] Владимиров А.М., Ляхин Ю.И., Матвеев Л.Т., Орлов В.Г. Охрана окружающей среды. Л.: Гидрометеоиздат, 1991. C. 87-116.
Поступила в редакцию 16.05.07 2007 г., в переработанном виде — 2.07.07 2007 г.
