CC BY
60
ДОКЛАДЫ АКАДЕМИИ НАУК РЕСПУБЛИКИ ТАДЖИКИСТАН __________________________________2009, том 52, №4______________________________
ЭКОЛОГИЯ
УДК 989.5: 201
Б.Р.Нахшиниев, Т.Сато ЧИСЛЕННОЕ МОДЕЛИРОВАНИЕ РАСПРОСТРАНЕНИЯ ПЫЛИ С ВЫСОХШЕЙ АКВАТОРИИ АРАЛЬСКОГО МОРЯ
(Представлено академиком АН Республики Таджикистан И.А.Абдусалямовым 18.08.2009 г.)
Аральское море оценено как одна из девяти опасных экологических зон, где здоровье и благосостояние населения резко ухудшилось в результате неблагоприятных экологических изменений [1]. Вследствие увеличения площади орошаемых земель и строительства гидроэлектростанций наряду с увеличением населения в регионе Аральское море деградировало по масштабам в самое большое внутриконтинентальное соленое озеро (с 10 г/л до 40 г/л) в мире. Уровень воды в море упал более чем на 15 м, объем воды уменьшился с 50 до 75%, сократилась поверхность моря, а площадь высохшего морского дна в настоящее время составляет более 33 000 км [2]. Хотя эта экологическая катастрофа называется “Проблемой Аральского моря”, кризис Арала включает более широкий спектр проблем и не ограничивается лишь самим озером. Большая область периферии моря также серьёзно пострадала и определяется исключительно как “Экологическая зона бедствия” [3].
Каждый год песчаные бури поднимают в атмосферу миллионы тонн солей и песков с высохшего морского дна и распространяют их на расстояние более чем сотни километров. Согласно данным Singer и др. [4], пески и пыль могут быть подняты на высоту более 4 км, разносятся по всей территории региона, а затем оседают в виде осадков и пыли.
Несмотря на эти факты, отсутствует комплексное изучение распространения и переноса пыли с высохшей акватории моря, распределения, сезонности, ее количества и направления переноса. По данному вопросу опубликован ряд статей [5-10]. Однако результаты, представленные в этих статьях, получены главным образом экспериментальным путём, а не современным методом вычислительного моделирования.
В настоящей статье представлены результаты численного моделирования распространения пыли с высохшей акватории Аральского моря, включая количество и направления переноса пыли. Сначала описывается система и обсуждаются результаты моделирования на фоне спутниковых наблюдений. Работа была выполнена на базе Лаборатории по моделированию окружающей среды Токийского университета.
Моделирование поля ветра
В странах Центральной Азии существуют национальные метеорологические станции, которые систематически проводят метеорологические наблюдения и мониторинг качества атмосферного воздуха. Несмотря на то, что станции обеспечивают наиболее реальную
Приблизительное местоположение
формацию в пределах своих территорий, они все еще имеют большой разрыв в местах, особенно там, где необходимо выяснить концентрацию загрязнения. В связи с этим в наших следованиях в качестве потенциального источника атмосферных данных были выбраны данные NCAR (Национальный центр атмосферных исследований США), записанные в 00:00, 06:00, 12:00, 18:00 ч с разрешением 2.50 широты и 2.50 долготы.
Другим наиболее важным фактором, с учетом которого можно достичь успеха в числовом моделировании атмосферного поля ветра над трехмерной топографией, является определение топографии модели как граничное условие в процессе вычисления. С этой целью мы использовали сеточный метод (рис.1) с координационной системой, приспособленной к рельефу местности. Модель состоит из 143 и 110 сеток в горизонтальном (Х и У) направлении с разрешением 15 км и 30 слоев в вертикальном ^) направлении с высотой 4 м в самом низком слое и до 2000 м в верхнем слое.
Обычно NCAR-атмосферные данные доступны в координационной системе с разрешением 2.50 широты и 2.50 долготы. Однако атмосферные данные в этой модели потребовались в координатной сети с пространственным разрешением 15x15 км. Поэтому данные с географической координационной системы сначала были преобразованы на декартовую координационную систему и методом бикубического сплайна интерполированы на более мелкую сетку.
В вертикальном направлении рассматривался линейный логарифмический профиль ветра для того, чтобы использовать интерполированные данные и эффект шероховатости на интенсивность и направление ветра. Эти величины также зависят от атмосферной стабильности (нейтральная, стабильная и неустойчивая) поверхностного атмосферного слоя. Выше поверхностного слоя (Иповерх.), проведена линейная интерполяция с использованием гео-строфического ветра. Логарифмический профиль ветра в вертикальном направлении дан по формуле
Рис.1. Моделирование рельефа местности района Аральского моря с использованием сеточного использованием сеточного метода. В центре - Аральское море.
/ \ и*
и{ г) = —
к
108--^ (у)
V
где и* - скорость поверхностного трения, к - безразмерный коэффициент Вон-Кармана, -высота шероховатости, г - высота, где определяется скорость ветра, и Ь - длина Обухова.
Значение Х^т зависит от стабильности атмосферы.
Чтобы поддержать теорию гладкого изменения направления ветра выше поверхностного слоя, линейная интерполяция с геострофическим ветром иГ выражена следующими формулами [11]:
и (г) = р(г)и(кРШЙСЧ) + \-р{г)Ъ
Г п ^ РПВЕСЧ. — 2 — Ьрсс ’
где р(£) находится по формулам
z-h
РПВЕСЧ
h^-h
К -К
p(z) = l-
V hrcc 'hРПВЕСЧ. J \ hrcc hРПВЕСЧ J
Здесь высота поверхностного слоя определена по формуле [12]:
/2 = И /10
,1 РПВЕСЧ *1ГСС 1 5
где (ИГРс) - высота атмосферного граничного слоя и выбрана так, что интенсивность и направление ветра всех слоев той высоты одинаковая. В заключение модель предполагает, что и (г) =иГ , если z > Игрс, и и(г) = и (Иповерх), если г < Иповерх. .
Расчет вертикальной скорости ветра
Вертикальная скорость ветра - обычно не измеряемая переменая, и ее расчет до настоящего времени считается одной из самых трудных проблем для метеорологов. Однако, если данные поля ветра доступны на основе координатной сетки, расчет кинематической вертикальной скорости ветра можно получить путем интеграции уравнения непрерывности:
Кдл ду
dz.
Расчет гео строфичес к ой скорости ветра
Геострофическая скорость ветра - это функция, которая зависит от широты, градиента атмосферного давления и плотности воздуха, которая не измеряется регулярно. Однако существуют эквивалентные способы расчета геострафической скорости ветра. По Holton [13], геостро-фический ветер выражен в более подходящей форме, которая имеет то преимущество, что не учитывает плотность воздуха. Используя геострофическое равновесие, геострофический ветер можно определить по формуле:
f0=-«az
/ AN
где / - параметр Кориолиса, А/ - разница в высоте между соседними контурными линиями, АЛ - расстояние между контурными линиями.
Моделирование переноса пыли
С точки зрения управления качеством воздуха, наибольшая угроза человеческому здоровью представляет концентрация частиц пыли меньше 10 /ш. Однако, с целью получения наиболее важных особенностей пыльной бури, верхний предел размера частицы пыли в нашем исследовании равен 100¡мп в диаметре. Во время вычисления движения частиц размеры диаметра частиц выбраны случайно на основе функции Гауссовского распределения вероятности со средним отклонением 50 и среднеквадратичным отклонением 60 и 10 микрон, соответственно.
N - количество частиц, выпущенных из источника каждого временного шага. Предполагается представить вертикальный поток пыли, испущенный с поверхности земли в данном временном шаге. В этой модели рассматривается только сухое гравитационное осаждение. Хотя известно, что мокрое осаждение является важным процессом, существует большая неопределенность в разработанных системах. Учитывая эти недостатки и ограничения в существующей системе, в данном исследовании не был использован процесс мокрогбисшкдентареноса предсказывает концентрацию пыли, решая следующую систему уравнения в Лагранжевой структуре:
г 1
* =-УрЧР + шр%-^рг8С0иг |иг|
тп +■ р 2
V 2 У
ёХр ’
<и Пр
где хр обозначает местоположение частицы, тр - масса частиц, р; - плотность воздуха, V -объем частиц, и - скорость перемещения частиц, VР - градиент давления, g - сила тяжести, 5 - поперечное сечение частиц в направлении и г, Св - коэффициент аэродинамического сопротивления, и г - относительная скорость частиц к воздуху.
Значение коэффициента аэродинамического сопротивления (Св (Яер) ) зависит от режима течения воздуха вокруг частиц. Соотношение Сс (Яе ) аппроксимативно было установлено, используя следующее функции:
к к
С о (Яе )= + —2— + кг.
Число Рейнольдса частиц пыли определено по формуле:
Кєр=-
и„ -пЮг
где Vап. - вязкость воздуха и Б - диаметр частиц.
Результаты моделирования
Значительное усилие было предпринято, чтобы разработать модель, которая могла бы предсказать особенности переноса пыли, испускаемой с высохшей акватории Аральского моря. Модель полностью программирована на компьютерном языке ФОРТРАН и была запущена в чение одного года (январь-декабрь 2003 г.) с временным шагом 3 ч. Мы не претендуем, что наша
модель совершенно точно
Рис.2. Пыльная буря в районе Аральского моря: а) - модель переноса пыли, построенная б) с использованием спутникового снимка NASA.
зывает перенос пыли, поскольку вершенствованная модель переноса пыли все еще не разработана из-за существования сложных
денных физических процессов во время ветровой эрозии. Учитывая описанные сложности, результаты этого исследования очевидны: тема достаточно аккуратно рует перенос пыли с высохшего дна Аральского моря.
Моделируемые переносы ли достаточно хорошо согласуются при сравнении со спутниковыми изображениями песчаной бури 18 апреля 2003 г (рис.2), опубликованных в интернет-сайте NASA [14]. Как видно, формы и направления переноса частиц пыли совпадают с моделированным полем
странения частиц песка и пыли.
Как видно из моделируемой поверхностной температуры
(рис.3), главной характерной особенностью поля ветра является сильный северо-восточный ветер, связанный с сильным прохладным воздушным потоком в течение почти целого дня 18
апреля 2003 г. Местоположение фронтальной системы с острыми температурными градиен-ми можно легко установить с помощью рисунка.
За моделируемый период 2003 г. модель сумела обнаружить количество и основные правления песчаных бурь с высохшего дна Аральского моря. Установлено, что основными транспортными направлениями пыли, с частотой более 12, были Ю, Ю-ЮЗ, ЮЗ и З-ЮЗ. Перенос пыли на направления В, В-ЮВ, ЮВ, Ю-ЮВ был обнаружен 6 раз, а на С, С-СВ, СВ, В-СВ - 4 раза.
Перенос пыли на направления З, З-СЗ, СЗ и С-СЗ
установлен только 2-3 раза.
Институт зоологии и паразитологии им. Е.Н.Павловского Поступило 18 августа 2008 г. АН Республики Таджикистан,
Лаборатория по моделированию окружающей среды Токийского университета, Япония
ЛИТЕРАТУРА
1. UNEP/CBD/SBSTTA. - Habitat Destruction in the World’s Waters: Progress report by the Global International Waters Assessment. Montreal, 2003, p.5.
2. IFAS. - International Fund for saving the Aral Sea, 2002 (http://www.ec-ifas.org).
3. Micklin, P., Williams W.D. - NATO ASI Series 2: Env., Berlin, 1996, v.12, pp.5-6.
4. Singer A., Zobeck T. et al. - J. of Arid Env., 2003, №54, рр.705-728.
5. Denisov Yu.M. - J. Prob. Des. Dev., 1998, №3, pp.21-25.
6. Esenov P. - J. Prob. Des. Dev., 1998, №34, pp.70-75.
7. Tolkacheva G.A., Kovalevskaya Yu.I. et al. - J. Prob. Des. Dev., 1998, №3, pp.55-59.
8. Kamalov S. - J. Prob. Des. Dev., 1998, №3, pp.25-35.
9. Razakov R.M., Kosnazarov K.A. - NATO ASI Series 2: Env., Berlin, 1996, v.12, pp.17-22.
10. Semenov O.E. - J. Prob. Des. Dev., 1998, №3, pp.37-44.
11. Montero G., Montenegro R., Escobar J.M. - J. of Win. Eng. & Ind. Aer., 1998, v.74, pp. 249-261.
12. Holton J. R. - An Introduction to Dynamic Meteorology. - New York, Acad. Press, 2004, p.74.
13. Gillette D., Ono D., Richmond K. - J. Geophys. Res., 2004, p.109.
14. http://visibleearth.nasa.gov/ view_rec.php?id=5252
271.42 272.87 274.51 276.06 277.60 279.15 260.68 282.23 283.76 265.32 266.87
Рис.3. Моделируемая поверхностная температура и характерные особенности поля ветра в районе Аральского моря 18.04.2003 г.
Б.Р.Нахшиниев, Т.Сато ТАХ,РСОЗИИ ПА^НШАВИИ ЧАНГУ FYEOF АЗ САТ^И ХУШКГАРДИДАИ БА^РИ АРАЛ
Дар мак;ола натичах,ои тахрсозии ракдмии пахдшавии чангу губор аз сатхи хушкгардидаи бахри Арал, бо назардошти микдор ва самтх,ои харакати чанг мух,окима гардидааст.
B.R.Nakhshiniev, T.Sato NUMERICAL SIMULATION OF WINDBLOWN DUST EMISSION FROM DRIED-UP BOTTOM OF ARAL SEA
This article describes the results of the modeling of windblown dust emission from dried-up bottom of Aral Sea including its quantities and transport directions.
