CC BY
34
Спутниковая гидрофизика
УДК 504(262.5)
Н.А. Тимофеев, А.В. Юровский
Климатический мониторинг атмосферных осадков и радиации на Черном море по спутниковым данным
Исследуются связи между радиацией и атмосферными осадками на поверхности Черного моря, с одной стороны, и параметрами, их определяющими, с другой. Показано, что по температуре поверхности воды и облачности с учетом ее оптической плотности можно с высокой точностью восстанавливать поля месячных и годовых величин атмосферных осадков, радиационного баланса и поглощенной в воде суммарной солнечной радиации. Исходная информация для этих целей определяется на основе данных по спектральным и интегральным потокам уходящих в космос длинноволновых излучений, измеряемых с ИСЗ.
В ХХ в. на Черном море наблюдались линейные тренды, характеризующие возрастание атмосферных осадков на 6 см и уменьшение радиационного баланса на 80 МДж-м-2 соответственно за 68 и 72 года. Температурный режим поверхности моря оказался стабильным, что обусловлено противоположным влиянием двух факторов: атмосферных осадков и глобального парникового эффекта.
В процессе спутникового мониторинга в XXI в. необходимо на основе предложенной модели сформировать банк данных с целью более точного описания климата и многолетней изменчивости полей облачности, температуры поверхности воды, радиации и атмосферных осадков.
Введение. Изучение влагооборота и радиационного режима в системе океан (море) - атмосфера (СОА) принадлежит к числу приоритетных задач глобальной климатологии [1]. В отличие от речного и подземного стоков пресных вод с суши, а также испарения, которые вычисляются на основе многолетней метеорологической и гидрологической информации, количество атмосферных осадков AO, выпавших на открытые акватории морей и океанов, оценивалось по данным, полученным с маяков, островных и прибрежных станций, с учетом их экстраполяции и коррекции по данным эпизодических судовых наблюдений. Для Черного моря разнородные ряды AO с 27 отечественных и 10 зарубежных станций [2] откорректированы поправками для приведения дождемерных наблюдений к осадкомерным и на недоучет осадкомерами ветрового выдувания, смачивания и испарения [3, 4]. При построении карт AO использовались косвенные оценки, полученные в результате среднеарифметического или средневзвешенного осреднения значений, экстраполированных по данным с выборочной сети береговых станций, и планиметрирования изогиет с учетом судовых наблюдений. Кроме того, количество атмосферных осадков определялось как остаточный член уравнения водного баланса [2].
Ввиду фрагментарности данных по осадкам в задачах диагноза климата используются средние месячные и годовые величины, полученные в 1923 -1940, 1945 - 1985, 1986 - 1998 и 1999 - 2002 гг. [2]. Последние два периода характеризуются соответственно повышенным (61,6 см год-1) и понижен© Н.А. Тимофеев, А.В. Юровский, 2009
68
ным (40,0 смтод-1) уровнями осадкообразований. Средняя для акватории Черного моря климатическая величина AO составила к 1998 г. слой воды в 55 - 57 смтод- с хорошо выраженным годовым ходом месячных величин. В интервале с абсолютными минимумом (38,2 см-год-1) и максимумом (84,8 см-год-1) осадков выделялись колебания с амплитудами больше 12 -13 сммес-1. Несмотря на недостатки указанных методик, полученные оценки AO признаются близкими к реальным величинам [2, 5].
Главный недостаток упомянутых способов - их непригодность для мониторинга осадков. В этом отношении перспективен метод [6 - 8], основанный на использовании спутниковых наблюдений, непрерывных по пространству и времени. Из космоса обнаруживаются области экстремумов, зоны ослабления или усиления осадков и радиации, их изменения во времени. Для расчетов AO и радиации используются три вида информации: сведения о температуре поверхности воды (ТПВ) t [9], разности температур вода - воздух (t -t) и об эффективной облачности ЭО = n f , численно равной произведению общего количества облаков (0 < n < 1) на их условную оптическую плотность (0 < f < 1). Температуры воды и воздуха характеризуют энергетику процессов осадкообразования и теплообмена океана (моря) с атмосферой [6, 10]; ЭО - регулятор термодинамического состояния СОА при облачно-радиационных взаимодействиях [8], а AO являются источником доступной потенциальной энергии в атмосфере за счет скрытой теплоты конденсации. Рассмотрим вышеуказанные методы.
Атмосферные осадки. Среди физических механизмов, определяющих формирование поля AO на поверхности Черного моря на сезонных масштабах, можно выделить следующие [10].
1. Фоновые поля, связанные с крупномасштабной циркуляцией атмосферы, безусловно определяемые с ИСЗ. Летом над морем, особенно в северозападной части, формируются преимущественно антициклоны с дивергенцией воздушных потоков у поверхности; зимой ситуация противоположная -циклоническая завихренность с конвергенцией. Результаты реанализов NCEP/NCAR по скорости ветра на стандартных горизонтах 10 м, 850, 500 и 200 мбар, относящихся к периоду 1958 - 1997 гг., действительно показали сезонную изменчивость циркуляционных процессов, проявляющуюся в смене циклонической завихренности ветра в холодную половину года на антициклоническую в теплый период.
2. Сезонные эффекты термической и динамической природы, обусловленные наличием в регионе высоких гор.
3. Муссонные эффекты, определяемые потоками тепла и влаги, формируемыми над морем и окружающей его сушей.
Таким образом, пространственно-временные распределения AO на поверхности Черного моря связаны как с особенностями атмосферной циркуляции, так и с орографией побережья. Главной причиной выпадения осадков является циклоническая активность в регионе. Преобладающее число атлантических и средиземноморских циклонов проходит с запада на восток над
южными районами моря. В течение всего года южные, в особенности юго-восточные, районы отличаются наибольшим количеством осадков. Это обусловлено наличием на пути движения влагонесущих потоков воздуха горных хребтов Кавказа, способствующих задержанию и накоплению влаги в атмосфере, при этом существенно повышается количество и оптическая плотность облаков. Под влиянием орографии побережья на Южном берегу Крыма также выпадает несколько большее количество осадков по сравнению с соседними районами. Меньше всего осадков выпадает в центральных районах моря, особенно мало - в его северо-западной части. В течение года около 40% осадков выпадает зимой, меньше всего (порядка 15%) - в конце весны и начале лета, когда под влиянием стационарного Азорского барического максимума атлантические антициклоны, чередой перемещаясь с запада на восток, стабилизируются в северо-западной части региона. Здесь в узкой прибрежной полосе наблюдается сравнительно небольшое увеличение осадков, обусловленное термической конвекцией.
Основная формула для расчета атмосферных осадков, выпавших на свободную ото льда поверхность морей и океанов, записывается в виде [6 - 8]:
АО = а ЭО + Ь ЭО2 + с ЭО3 + й ЭО4 + ЭО Ф(Д). (1)
Коэффициенты а , Ь, с, й для Северного и Южного полушарий в отдельности задаются нелинейными комбинациями от ТПВ, а функция Ф(Д) описывает влияние на формирование поля осадков апвеллингов, морских и атмосферных течений, обуславливающих пространственно-временную изменчивость температур поверхности воды I и воздуха t' в приводном слое атмосферы с учетом их высокой взаимной коррелируемости. Для климатических расчетов величин осадков, выпадающих на поверхности Черного моря (при размерности скорости выпадения осадков - ммсут-1) функция Ф(Д) ~ 0,3 ^ - . Здесь перепад температур вода - воздух ^ - имеет годовой ход, при среднем значении Дt = 1,1 °С [11].
В работах [6 - 8] на основе многолетней спутниковой информации построены и опубликованы климатические карты атмосферных осадков над океанами с пространственным разрешением 500x500 км, по времени - месяц и год, погрешность ~ 5 - 10 мм мес-1. Исследована пространственно-временная изменчивость АО, в частности, при Эль-Ниньо [7]. Эти события в 1982 - 1983 и 1997 - 1998 гг. были самыми мощными в XX в. Аномалии осадков в экваториальной зоне Тихого океана достигали ± 350 - 400 мммес-1.
В горизонтально-однородной системе подстилающая поверхность - атмосфера среднегодовые климатические величины эффективной облачности ЭОп можно найти по формуле глобального парникового эффекта [12]:
0,27+0,0Ш = 0,516-0,060ЭОп -0,456ЭОП2 . (2)
Откуда имеем
70
ЭОп = -0,066+ д/0,5436- 0,0263161
(3)
При климатической температуре 1 = 15 °С эффективная облачность ЭОП = 0,320. Результирующее влияние ЭОП на глобальный климат является охлаждающим. При переходе от условий безоблачной атмосферы ЭОП = 0, когда 1 = 20,5°С, к состоянию сплошного облачного покрова предельно высокой оптической плотности ЭОП = 1 глобальная температура понизится относительно современного уровня на 43 °С и составит 20,5 - 43 = - 22,5°С - аналог эпох оледенения, неоднократно происходивших в истории Земли. Парниковый эффект в этом случае перестает оказывать влияние на климат.
При расчетах ЭОпнад Черным морем по формуле (3) использовались данные по ТПВ за 1986 - 2002 гг. с шагом по пространству 18 км, полученные с помощью радиометров высокого разрешения, установленных на орбитальных спутниках NOЛЛ-7, -9, -11, -14 и -16 [9]. Среднее по акватории моря за указанный период времени годовое значение 1 = 14,80°С в точности соответствует аналогичной величине, полученной обобщением данных судовых наблюдений за 1928 - 1940, 1945 - 1985 гг. [11]. Температурный режим поверхности Черного моря оказался относительно стабильным.
На рис. 1 представлены климатические карты среднегодовых значений ЭОП, обусловленных парниковым эффектом (3), величин ЭО, рассчитанных в соответствии с уравнением (1) по ТПВ, Д1 = 1,1 °С и атмосферным осадкам, а также их разности ДЭО = ЭО - ЭОи . На рис.2 показаны климатические карты среднегодовых значений атмосферных осадков АО, АО и их разности ДАО = АО - АО; величины АО заимствованы из справочного издания [2]. Разности ДЭО = ЭО - ЭОп и ДАО = АО - АО характеризуют влияние атмосферной циркуляции и орографии побережья в регионе на пространственное распределение облачности и осадков. В табл. 1 годовые характеристики температур, облачности и осадков выделены по различным временным интервалам. По этим данным построен график (рис. 3) зависимости разностей осадков (АО - АО) от разностей облачности (ЭО - ЭО )• С учетом данных по 20 выделенным временным подмассивам климатического периода 1923 - 2002 гг., с разрывом в 1941 - 1944 гг., имеем почти универсальное соотношение (АО - АО)/(ЭO - ЭО) ~ 388 см-год"1, по которому средние, максимальные и минимальные разности ДАО определяются с ошибками не более 5%. Как видно из рис.1 и 2, пространственное распределение нулевых разностей ДАО и ДЭО соответствует положению климатической изогиеты АО « 55 - 57 см-год-1, проходящей между 42 и 43° с.ш., а восточнее 37° в.д. она резко смещается к северу, по направлению к Керченскому проливу, отделяя области с преобладанием значительных осадков на юге и, в особенности, в восточной части моря. Здесь под влиянием циклонов и гор Кавказа в атмосфере формируются области высоких значений ЭО и АО •
71
Р и с. 1. Климатические карты среднегодовых значений эффективной облачности ЭО , ЭО и их разности АЭО = ЭО - ЭО
Р и с. 2. Климатические карты среднегодовых значений атмосферных осадков АО , АО и их разности ДАО = АО - АО, см год-1 (штриховой линией на верхней карте показано пространственное положение годовой климатической изотермы 14,6 - 14,8°С, совпадающее с положением изогиеты АО ~ 56 см год-1)
73
Т а б л и ц а 1
Годовые данные по ТПВ (¿,°С), атмосферным осадкам (АО, см) и эффективной облачности (ЭО), осредненные по акватории Черного моря за различные временные интервалы
Годы г ^шах ^шт АО ЭО эоп АОп
1923 - 1940 1945 - 1985 14,80 16,4 13,4 56,0 0,322 0,326 54,0
1986 14,75 16,0 13,2 50,9 0,318 0,328 54,7
1987 13,54 15,4 12,6 64,7 0,350 0,367 71,6
1988 14,34 16,0 13,2 84,8 0,399 0,342 60,6
1989 14,69 16,2 13,0 53,8 0,325 0,330 55,6
1990 14,71 16,2 13,2 55,8 0,331 0,330 55,5
1991 14,54 16,6 13,8 60,7 0,343 0,336 58,1
1992 14,12 15,5 13,2 65,7 0,354 0,350 64,1
1993 14,00 15,2 12,6 38,2 0,278 0,353 65,4
1994 15,49 16,6 13,8 51,2 0,321 0,304 45,1
1995 14,69 16,4 13,6 76,0 0,380 0,330 55,6
1996 14,88 16,8 14,0 59,6 0,340 0,324 53,1
1997 14,49 17,0 13,2 78,8 0,386 0,338 58,9
1998 15,36 16,8 13,4 61,6 0,348 0,310 47,4
1999 - 2002 15,46 17,4 14,3 40,0* 0,290* 0,306 45,8
1923 - 1998** 14,75 16,4 13,4 57,1 0,344 0,328 54,7
1986 - 1998 14,60 16,4 13,4 61,6 0,345 0,334 57,4
1986 - 2002 14,80 16,4 13,4 56,0 0,332 0,326 54,0
1923 - 2002** 14,80 16,4 13,4 56,0 0,332 0,326 54,0
П р и м е ч а н и е: * - восстановленные величины, ** - отсутствие наблюдений в 1941 - 1944 гг.
74
с <
I
с с
с
-30
-0.12 -0.08 -0.04 0 0.04 0.08
ЭО-ЭОп
Р и с. 3. Соотношение между разностями (AO - АО ) и (ЭО - ЭОп ), осредненными по акватории Черного моря за различные временные интервалы (данные табл. 1)
Формула (3) в комбинации с выражением (1) обеспечивают корректное вычисление годовых климатических норм атмосферных осадков, обусловленных полем ТПВ, за 1986 - 2002 гг. [9]. Пространственные положения изогиеты AO « 56 см-год-1 и климатической изотермы 14,6 - 14,8°С совпадают (см. рис. 2).
С учетом влияния гор Кавказа и особенностей атмосферной циркуляции в регионе для расчетов по известным значениям ЭО и ТПВ как климатических, так и текущих месячных и годовых величин AO и их пространственных распределений целесообразно использовать соотношение (1). Текущие месячные и годовые значения ТПВ в период 1986 - 2002 гг. изменяются по акватории Черного моря относительно их климатических величин в сравнительно узких интервалах. По формуле (1) на основе данных о фактических осадках [2], климатических значений ТПВ [9] и At [11] найдем ряды эффективной облачности, а затем получим соотношения вида АО(ЭО) и ЭО(AO) . Из табл. 2 следует, что «переходы» ЭО ^ АО характеризуются среднеквад-ратическими погрешностями аппроксимаций: для осадков aAO = 0,05 -0,55 см мес-1 и 0,52 - 1,01 см год-1; соответственно аЭО = 0,006 - 0,024 и 0,003 для отдельных месяцев и лет. Погрешности aAO более чем на порядок меньше величин a'AO (рис. 4), характеризующих природные межгодовые вариации атмосферных осадков.
Р и с. 4. Межгодовые колебания месячных и годовых сумм атмосферных осадков на Черном море в период 1986 - 1998 гг. (сг'АО - среднеквадратические величины)
По данным наблюдений (рис. 5) временной тренд годовых сумм осадков
ЛЛО( И) =
1 И
- У Л01 - 51 N У
г=1
относительно начального уровня ~ 51 смгод-1 в 1923 - 1926 гг. составил к 1998 г., т.е. за 68 лет, приращенный слой воды в 6 см со скоростью выпадения осадков 0,088 смгод-1. Таким образом, на Черном море имела место линейная тенденция роста атмосферных осадков с 51 смгод-1 в 1923 - 1926 гг. до 57 смгод-1 к 1998 г.
ААО(М), см/(М лет)
1920 1930 1940 1950 1960 1970 1980 1990 2000
Р и с. 5. Временные тренды атмосферных осадков на Черном море (получены скользящим осреднением по одному году; N - годы наблюдений)
76
Т а б л и ц а 2
Параметры регрессионных зависимостей между месячными и годовыми величинами атмосферных осадков ( АО, см) и эффективной облачности ( ЭО ) для акватории Черного моря при заданных разностях температур вода - воздух [11] для соответствующих интервалов ТПВ [9], наблюдавшихся с ИСЗ
1986 - 1998 гг. АО = а' ЭО + Ъ' ЭО2 + с' ЭО аАО а'лО ЭО = а" АО + Ъ" АО 2 + с" АО3 аЭО Интервал ТПВ (°С) Д/ = =/ - /' (°С)
а' Ъ с' а" Ъ" с"
Январь -7,6 82,4 -55,9 0,26 4,0 0,124 -0,0135 0,00060 0,013 7,1 - 9,9 2,7
Февраль -0,5 42,0 -14,7 0,30 3,6 0,163 -0,0241 0,00137 0,009 5,8 - 7,5 1,7
Март -3,9 46,4 -12,6 0,36 3,1 0,163 -0,0219 0,00114 0,024 6,0 - 7,5 0,3
Апрель -3,3 48,9 -14,0 0,21 2,6 0,184 -0,0344 0,00244 0,016 7,3 - 9,7 -0,3
Май -6,0 64,0 -17,8 0,05 1,3 0,204 -0,0505 0,00477 0,006 11,2 - 15,1 0,1
Июнь 18,6 -90,8 221,6 0,33 1,5 0,166 -0,0312 0,00232 0,016 18,5 - 21,9 0
Июль 18,0 -54,0 160,0 0,41 1,7 0,095 -0,0092 0,00019 0,019 21,9 - 24,1 -0,6
Август 17,9 -17,1 96,8 0,50 2,6 0,090 -0,0137 0,00099 0,019 22,8 - 25,5 -0,3
Сентябрь 8,5 1,6 69,9 0,55 3,5 0,102 -0,0090 0,00032 0,019 20,3 - 23,8 1,3
Октябрь 4,9 13,4 47,2 0,32 3,7 0,116 -0,0117 0,00044 0,017 16,8 - 22,0 2,0
Ноябрь -6,5 81,9 -44,8 0,10 4,5 0,109 -0,0105 0,00041 0,009 11,7 - 15,9 2,7
Декабрь -8,8 94,3 -66,2 0,35 4,7 0,102 -0,0088 0,00033 0,011 8,1 - - 12,6 3,1
Год
АО = 77 ЭО +17 ЭО2 + 793 ЭОъ (для АО < АОКл )
АО = -57 ЭО + 68 ЭО2 (для АО > АОкл )
0,76 12,6
ЭО = 9,24 • 10-3 АО - 6,9110 5 х х АО2 +1,95 -10 7 АО3 (для АО < АОкл )
ЭО = 0,042 + 40,0018+0.0015АО (для АО > АОкл )
0,003 13,4 - 16,4
1,1
0233-7584. Мор. гидрофиз. журн., 2009, № 1
77
Радиационный режим. Радиационный баланс на поверхности моря (МДж-м 2-сут согласно [8], записывается в виде
Яп =у (1 - ЭО) [ПСР - (6,8 - 0,31 1 + 0,01712 - 0,0003713 + 0,5 А1)], (4)
а поглощенная в воде суммарная солнечная радиация (ПСР ) в условиях «нормальной» безоблачной атмосферы и при фактическом расстоянии Земли до Солнца - в виде
пср = о = о — sin h
sin h + 0,04
(5)
[35-0,385 р + (0,01644р-0,2348) 8 при8> 0;! I 35 - 0,385р + 0,58 при8< 0 J ,
где sin h - среднее от восхода до захода Солнца значение синуса его высоты; р и 8 - географическая широта и склонение Солнца в градусах.
Характеристикой оптического состояния безоблачной атмосферы является свободный от «эффекта Форбса» параметр
у = 0,9 - 0,4 sin h + (0,1 + 0,4 sin h) [S (h)/SH (p, h)],
который вычисляется по общему влагосодержанию атмосферы p и интенсивности прямой солнечной радиации S(h) . Если S(h) = SM (p, h), то параметр у , равный 1, характеризует оптическое состояние «нормальной» океанической атмосферы. Величина QH (МДж-м^-сут1) затабулирована в работах [8, 13] с учетом p и спектральной оптической толщины атмосферы аЛ = 0,05 - 0,10 (Л = 0,66 мкм), типичной для акваторий, не граничащих с большими пустынями Земли. Кроме регионов с сильно замутненной атмосферой, где у значительно меньше 1, над остальной существенно большей частью Мирового океана, включая акваторию Черного моря, в среднем за год у = 0,99 [11]. Под влиянием сезонности выносов пыли из пустынных районов Африки на Европу параметр у для Черного моря имеет следующий годовой ход: от значений 1,022 в декабре до 0,966 в августе [11].
Эффективная облачность и ТПВ в годовом ходе коррелируют между собой. Формулы (4) и (5) непосредственно предназначены для расчетов месячных величин радиационного баланса, а 12
Rn,год ^ ^и,мес .
1
На рис. 6 показана карта климатических годовых сумм радиационного баланса. Необходимые для расчетов значения эффективной облачности (рис. 1) определены по атмосферным осадкам (рис. 2); данные по ТПВ получены из спутниковых наблюдений [9], а значения Al взяты из табл. 2. В среднем для всей площади моря Rn = 3200 МДж-м2-год1 или 8,8 МДж-м^-сут1, что несколько меньше оценки 9,1 МДж-м-2-сут-1, полученной в [11] при осреднении за 1923 - 1985 гг. с учетом разрыва в наблюдениях 1941 - 1944 гг. Вблизи Кавказского побережья величины Rn на годовой карте из работы [11]
представляются существенно завышенными, что является следствием использования в этих расчетах данных только по общему количеству облаков п без учета их оптической плотности / .
Р и с. 6. Климатическая карта радиационного баланса на поверхности Черного моря К, МДж-м 2тод 1
На фоне временного тренда годовых величин радиационного баланса порядка -80 МДж-м-2 за 72 года (1923 - 1940, 1945 - 1998 гг.) имеют место значительные его межгодовые вариации 220 МДж-м-2-год-1 (рис. 7).
Р и с. 7. Межгодовая изменчивость средних месячных (точки) и годовых (ромбы) величин радиационного баланса Черного моря
Климатические величины атмосферных осадков и аномалий радиационного баланса
АЯп = К - К = К - К (1 - ЭО)=к ЭО
(6)
(Л0 соответствует условию ЭО = 0), определенные по формулам (4), (5) и полученные по картам на рис. 2, 6 в узлах одноградусной координатной сетки, тесно связаны между собой (рис. 8):
АО = 1,27-10 ДЯ„ + 5,91 • 10 + 4,6110^ Щ, для АО < АОкл I АО = -1,27 АЯ + 0,291-10-4 АЯ2 для АО > АО ) '
при этом 5АО = 4,9 смгод 1
Р и с. 8. Соотношения между годовыми величинами атмосферных осадков и соответствующими аномалиями радиационного баланса ДЯп = Я0 ЭО на Черном море
Осредненная по акватории Черного моря и по времени за периоды 1923 - 1940, 1945 - 1998 гг. аномалия ДКп = 1630 МДж-м-2-год-1. Этой величине, согласно формулам (7), соответствует современное климатическое значение атмосферных осадков в 57 см год-1. Тенденция увеличения АО с 51 до 57 см год-1 за 68 лет (рис. 5) согласуется с временным трендом радиационного баланса ДКп (К) « -80 МДж-м-2 за временной интервал в 72 года (рис. 7).
При максимуме эффективной облачности ЭО = 0,65 аномалиям ДЯп = 1630^°65 =
3160 МДж-м-2-год-1 соответствуют осадки « 250 см год-1, выпадающие в районе г. Батуми [2]. Это немного меньше аналогичных величин, которые наблюдаются в тропической зоне Мирового океана в сезон дождей [6 - 8].
Таким образом, если из наблюдений с ИСЗ при известных климатических значениях Д^ [11] имеем надежные данные по ТПВ и эффективной облачности, можно будет с высокой точностью проводить на Черном море мониторинг согласованных между собой месячных и годовых величин атмосферных осадков АО и радиации (Яп и Qn).
Эффективная облачность и ТПВ. Величины ЭО определяются следующим образом [10]:
ЭР= К ~ = (22,6 + 0,361 -0.00512 - 0,0001213) -^ (8)
К - 11,2 + 0,261 •
Интегральные по спектру Я = 3,8 - 50 мкм среднесуточные полусферические потоки ^ уходящей в космос длинноволновой радиации (УДР) при произвольных условиях облачности измеряются широкоугольными приборами спутников или вычисляются по интенсивностям излучений Д (©) в спектральной полосе Я = 10,3 - 11,3 мкм. Второй метод основан на существовании тесных природных связей (коэффициенты корреляции 0,96 - 0,99) интенсивностей излучений Д (©) , /2 (©), /3 (©) в спектральных полосах Я = 10,3 - 11,3 мкм, Я = 8 - 12 мкм, Я = 3 - 30 мкм между собой, а также со среднесуточными полусферическими потоками ^ [13 - 20]. Максимальные ^ и минимальные потоки УДР, соответствующие условиям ЭО = 0 и ЭО = 1, параметризованы по ТПВ с учетом данных по температуре воздуха на уровне тропопаузы в связи с их высокой взаимной коррелируемостью [8].
Пусть над океаном (морем) для точки, соответствующей углу визирования © (град), в участке спектра Я = 10,3 - 11,3 мкм получена (без фильтрации облачности) радиационная температура ^ (©), которой, согласно формуле Планка, соответствует интенсивность длинноволнового излучения [113]
0 81
/1<©>- 13322 7 (9)
ехр--1
273,1 + ^(©)
размерностью кВт-м-2 в телесном угле 1 ср, оконтуривающем на поверхности водоема площадку в один пиксел. Под ^(©) понимается температура поверхности абсолютно черного тела, при визировании которой на выходе измерительного комплекса в спектральном интервале Я = 10,3 - 11,3 мкм создается такой же сигнал, как и от наблюдаемого пиксела. Последовательно переходим от /1 (©) к интенсивностям длинноволновых излучений /2, 13 в подспутниковой (© = 0) точке соответственно для спектральных интервалов Я2 = 8 - 12 мкм, Я = 3 - 30 мкм и далее к среднесуточному полусферическому потоку ¥сп , МДж-м2-сутЛ Согласно [8, 10, 13 - 20],
Л!2(© = 0) = /а =(3,66-12-10-7 ©3)ц©-2-10-8 ©3, (10) 8-12 2 1 -11110-8 ©3 , ( )
/3_3 0(© = 0)= /3 =(0,716 /2 + 0,0024)0,616 - 0,0122, (11)
^ =370/3 -1590/3 + 4980/1, (12)
^ =240/3 + 2,8 , (13)
^ = 443/3 - 2948/3 + 9970/1. (14)
Варианты (12) - (14) в пределах — 0,5 МДж-м 2-сут равнозначны по точности. По трем значени-
ям средняя величина
Fc„ = 351/3 -1513/32 + 4983/З3 + 0,93. (15)
Суточный ход УДР над морями и океанами, по данным с ИСЗ «Метеор» [13], не обнаруживается. Этот вывод подтвержден наблюдениями с ИСЗ «ШтЪш-7» [21], который находился на солнечно-синхронной орбите и пересекал экватор в 12 ч дня и ночи по местному времени. В области Мирового океана 60°с.ш. - 60°ю.ш. в течение 5 лет отношение (Рси )день/(Рси )ночь= 1,012 с разбросом ± 0,0077.
Алгоритмы различных модификаций восстановления ТПВ по наблюдениям с ИСЗ с учетом фильтрации облачности подробно проанализированы в статьях [8, 22 - 24]. В широкомасштабных судово-спутниковых экспериментах [20], проведенных в Атлантике в 1986 - 1991 гг., при пространственной изменчивости ТПВ от 6 до 30°С, влагосодержании о = 5 - 75 мм [25] среднеквадратиче-ские ошибки восстановления £ составили 0,2 - 0,3°С на углах сканирования 0 - 55°. При этом наблюдаемые аэрозольные оптические толщины атмосферы для длины волны X = 0,66 мкм изменялись от значений 0,05 - 0,10, близких к океанической «норме» [13], до предельно высокой величины 0,76 в «море мрака» [26], где имели место систематические выносы в атмосферу большого количества пыли из районов Сахары [17]. Ошибки восстановления ТПВ оказались независимыми от угла ©, влагосо-держания и условий аэрозольной замутненности атмосферы. Аналогичные результаты получены в судово-спутниковых экспериментах 1989 - 1993 гг. на Черном и Средиземном морях [27].
В настоящей работе данные по ТПВ Черного моря получены, как отмечалось, с орбитальных спутников №АА-7, -9, -11, -14, -16 за 1986 - 2002 гг. [9].
Температурный режим на поверхности Черного моря в течение последних 70 - 80 лет ХХ в. был относительно стабильным: £ср = 14,8°С. Некоторое понижение ТПВ, обусловленное возрастанием
эффективной облачности и осадков, компенсировано глобальным парниковым эффектом [12]. Однако за короткие интервалы времени ТПВ на Черном море может заметно меняться под влиянием региональной изменчивости в климатической системе. Так, например, в 11 -летнем цикле солнечной активности [9] в результате извержения вулкана Пинатубо на Филиппинах в 1991 г. [12] - самого мощного явления, зарегистрированного со времени начала спутниковых наблюдений, - на Черном море отмечено к 1993 г. понижение температуры на 0,7°С, а затем через 1 - 2 года произошло возвращение к климатической норме. После событий Эль-Ниньо в 1982 - 1983 гг. и 1997 - 1998 гг. [6 -8], которые были самыми сильными в ХХ в., на Черном море, наоборот, наблюдалось повышение ТПВ на 1 °С. По данным спутниковых наблюдений (табл. 1), самым теплым и сухим был период с 1999 г. по 2002 г.
Способ определения локальных текущих месячных величин Д = £ — £ на основе спутниковой информации по ТПВ и ЭО описан в работе [10].
Основные результаты и заключение. Алгоритмы, учитывающие угловую структуру поля УДР, параметризации Fc, Fmn соответственно при ЭО = 0 и ЭО = 1 [8], позволяют вычислять 0 < ЭО < 1 в любой географической точке морских акваторий. При этом используются среднемесячные величины Fcn и ТПВ. Информация, полученная с орбитальных спутников NOAA, ERS и других с шагом по пространству 18 - 50 км, содержит сведения о внутри- и межгодовых вариациях полей УДР, ТПВ и эффективной облачности в глобальном и региональном масштабах [8, 9]. Показано, что по этим данным для акватории Черного моря могут восстанавливаться с высокой точностью согласованные меж-
ду собой поля месячных и годовых величин атмосферных осадков, радиационного баланса и поглощенной в воде суммарной солнечной радиации.
В ХХ в. на Черном море имели место согласованные между собой линейные тренды, характеризующие возрастание атмосферных осадков на 6 см за 68 лет и уменьшение радиационного баланса на 80 МДж-м-2 за 72 года.
Температурный режим поверхности моря оказался стабильным, что обусловлено противоположным влиянием двух факторов: роста атмосферных осадков и увеличения глобального парникового эффекта [12].
Следующий этап - проведение на Черном море в XXI в. спутникового мониторинга с целью формирования более полного электронного банка данных, а также создание атласа климатических величин и характеристик изменчивости вышеупомянутых параметров.
СПИСОК ЛИТЕРАТУРЫ
1. Мировой водный баланс и водные ресурсы Земли. - Л.: Гидрометеоиздат, 1974. - 638 с.
2. Гидрометеорология и гидрохимия морей СССР. T.IV. Черное море. - СПб.: Гидрометеоиздат, 1991. - Вып. 1. - 429 с.
3. Методические указания. Устранение неоднородности между рядами дождемерных и осадкомерных наблюдений. - Л.: ГГО, 1964. - 31 с.
4. Методические указания. Введение поправок на недоучет осадков осадкомерами. - Л.: ГГО, 1966. - 42 с.
5. Горячкин Ю.Н., Иванов В.А. Уровень Черного моря: прошлое, настоящее, будущее. - Севастополь: МГИ НАН Украины, 2006. - 210 с.
6. Тимофеев Н.А., Юровский А.В. Атмосферные осадки над океанами: климат и изменчивость. Спутниковая информация. - Севастополь, 2002. - 40 с. - (Препринт / НАН Украины. МГИ).
7. Тимофеев Н.А., Юровский А.В. Реконструкция атмосферных осадков на океанах на основе спутниковой информации // Исследование Земли из космоса. - 2003. - № 3. - С.51 - 62.
8. Тимофеев Н.А., Юровский А.В. Радиационные, тепло- и водобалансовые режимы океанов. Климат и изменчивость. -Севастополь: МГИ НАН Украины, 2004. - 255 с.
9. Бабий М.В., Букатов А.Е., Станичный С.В. Температура поверхности Черного моря по спутниковым данным 1986 -2002 гг. - Севастополь: МГИ НАН Украины, 2005. - 263 с.
10. Тимофеев Н.А., Юровский А.В. О климатическом спутниковом мониторинге внешнего теплового баланса океанов и Черного моря // Морской гидрофизический журнал. - 2007. - № 6. - С 59 - 75.
11. Ефимов В.В., Тимофеев Н.А. Теплобалансовые исследования Черного и Азовского морей. - Обнинск: ВНИИГМИ-МЦД, 1990.- 237 с.
12. Тимофеев Н.А. Парниковый эффект атмосферы и его влияние на климат Земли (спутниковая информация) // Морской гидрофизический журнал. - 2006. - № 6. - С. 13 - 28.
13. Тимофеев Н.А. Радиационный режим океанов. - Киев: Наук. думка, 1983.- 247 с.
14. Тимофеев Н.А. Определение интегрального по спектру полусферического потока уходящей в космос длинноволновой радиации с помощью ИСЗ // Исследование океана дистанционными методами. - Севастополь: МГИ АН УССР, 1981. -С. 145 - 156.
15. Тимофеев Н.А., Ермоленко А.И., Шутова Е.Н. О корреляции интегральной и спектральных интенсивностей уходящего в космос длинноволнового излучения // Там же. - С. 157 - 166.
16. Тимофеев Н.А. Соотношения между интенсивностями уходящей в космос радиации в различных участках ИК-спектра // Морской гидрофизический журнал. - 1988. - № 1. - С. 41 - 46.
17. Тимофеев Н.А., Каба М.Л. Взаимодействие атмосферы и океана // Тропическая Атлантика. Регион Гвинеи. - Киев: Наук. думка, 1988. - 1. - Гл. 1 - 2. - С. 13 - 47.
18. Тимофеев Н.А., Иванчик М.В., Севостьянов А.И. Восстановление потоков радиации на поверхности океана по данным спутников NOAA // Морской гидрофизический журнал. - 1990. - № 3. - С. 30 - 36.
19. Тимофеев Н.А., Иванчик М.В., Севостьянов А.И. Солнечная радиация и радиационный баланс на поверхности океана по наблюдениям с искусственных спутников // Исследование Земли из космоса. - 1991. - № 2. - С. 15 - 25.
20. Тимофеев Н.А. Методы и результаты спутникового мониторинга гидрометеорологического режима Атлантического океана // Гидрофизика Тропической Атлантики.- Киев: Наук. думка, 1993. - С. 51 - 74.
21. Kyle HL., Ardanuy P.E., Hurley E.J. The Nimbus-7 Earth Radiation Budget data set. - Washington: NASA, 1985. - 248 p.
22. Тимофеев Н.А., Иванчик М.В., Севостьянов А.И. и др. Спутниковый мониторинг температуры поверхности океана по данным ИСЗ NOAA // Морской гидрофизический журнал.. - 1990. - № 5. - С. 39 - 48.
23. Тимофеев Н.А., ИванчикМ.В., Севостьянов А.И. Спектральный метод восстановления температуры поверхности океана по данным ИСЗ NOAA // Исследование Земли из космоса. - 1991. - № 3. - С. 82 - 88.
24. Тимофеев Н.А. Анализ алгоритмов восстановления температуры поверхности океана на основе наблюдений с оперативных спутников NOAA и ERS-1 // Там же. - 1992. - № 3. - С. 35 - 41.
25. Тимофеев Н.А., Юровский А.В. Определение температуры, влажности воздуха, скорости ветра, потоков явного и скрытого тепла на поверхности океана по наблюдениям с ИСЗ // Морской гидрофизический журнал. - 2007. - № 2. - С. 56 -70.
26. Фетт В. Атмосферная пыль. - М.: Иностранная литература, 1961. - 336 с.
27. Гришин Г.А., Тимофеев Н.А., Калинин Е.И. и др. Спутниковый мониторинг средиземных морей в оптическом диапазоне спектра // Исследование Земли из космоса. - 1991. - № 6. - С. 45 - 54.
Морской гидрофизический институт НАН Украины, Материал поступил
Севастополь в редакцию 07.09.07
ABSTRACT Relations between the radiation and the atmospheric precipitation on the Black Sea surface, on the one hand, and the parameters determining it, on the other hand, are investigated. It is shown that data on sea surface temperature and cloudine ss (with regard to its optical density) permit to reconstruct the fields of monthly and annual values of atmospheric precipitation, radiation balance and total solar radiation absorbed in water. Initial information for this purpose is obtained from the satellite-derived data on spectral and integral fluxes of the outgoing long wave radiation.
In the XX century the Black Sea was noted for the linear trends characterizing increase of atmospheric precipitation by 6 cm and decrease of radiation balance by 80 MJ-m"2 during 68 and 72 years, respectively. The temperature regime of the sea surface is found stable due to opposite influence of two factors: atmospheric precipitation and greenhouse effect.
Satellite monitoring in the XXI century provides an opportunity to create a databank based on the proposed model with the purpose to describe climate and long-term variability of the fields of cloudiness, sea surface temperature, radiation and atmospheric precipitation more correctly.
