Понятие климата. По определению [4, 7] «...климат - это статистический ансамбль состояний, проходимых системой атмосфера-океан-суша за периоды времени несколько десятилетий». В этом определении содержится два недостатка: в системе генерируется не один, а множество взаимосвязанных ансамблей разных пространственно-временных масштабов колебаний и противоречиво понятие «климатические характеристики». С одной стороны утверждается, что климатическими характеристиками считаются дескриптивные статистики («нормы»), полученные за стандартный период 30 лет [3]. Следовательно, характер изменения климата за 100 лет невозможно надежно оценить, поскольку на этом интервале мы имеем всего лишь три независимых значения. Однако далее говорится, что «...важнейшим параметром состояния климата является именно аномалия температуры у поверхности планеты, осредненная по всему Земному шару, или полушарию, или региону». Поэтому оценками изменения климата являются изменения аномалий. Причем, аномалии должны быть определены за период менее 30 лет, т.к. в противном случае они тождественно будут равны нулю. На практике изменения климата оценивается по значениям среднегодовых аномалий, усредненных по полушарию. При этом возникает еще одно противоречие: локальный климат есть проявление глобального климата [7]. Но осреднение по полушарию означает, что все локальные климаты тождественны между собой. А если средние аномалии вычисляются независимо для каждого региона, то это означает, что глобальный климат есть сумма локальных климатов, что противоречит определению климата в [7].
Альтернативное определение. Дадим другое понятие климата: климат - это глобальное макросостояние атмосферы. Измерение климата означает определение параметров макросостояния (климатических характеристик) по данным метеорологических наблюдений. Климатическая система представляет неразрывное целое. Атмосфера - это связующее звено всех внутренних подсистем климатической системы. В колебаниях ее характеристик отражаются изменения всей системы. Макросостояние можно характеризовать глобальными пространственными функциями текущих значений параметров годового хода основных процессов в атмосфере - мезомасштабных и иных флуктуаций. Характеристика мезомасштабных колебаний есть среднеквадратическое отклонение от функции суточного хода за данные сутки. Параметрами суточных колебаний являются наименьшее значение за сутки и амплитуда суточного хода. Каждый из этих параметров можно характеризовать за каждый год среднегодовым значением и амплитудой годового хода. Отклонения мезомасштабных флуктуаций от их годового хода обусловлены процессами синоптического масштаба за год. Эти отклонения также можно характеризовать среднеквадратическими отклонениями.
Для определения климатических характеристик (параметров глобальных полей) можно использовать системный синтез. Основные принципы системного синтеза для климатической системы, сформулированы в [5]. Они позволяют практически решать задачу реконструкции аттрактора состояния системы по эмпирическим данным [2].
Три основные предпосылки системного синтеза
Макросостояние системы при заданных внешних условиях является инвариантом. Повторение внешних условий порождает множество состояний, подобных инвариантной структуре [1,4].
Для описания макросостояния знание физических механизмов взаимодействия внутренних подсистем не является обязательным условием.
Каждому макросостоянию соответствует вполне определенная функция распределения микросостояний. Это соответствие можно найти эмпирически.
Поясним кратко сущность системного синтеза для определения климатических характеристик. Пусть ujk(φ,λ,t) - значения температуры воздуха на станции с координатами φ, λ за j-тый день k-того года, t - время в течение суток. Внешние условия для этого дня в году за весь период наблюдений будем считать постоянными. Тогда ujk(t) можно представить в виде:
ujk(φ,λ,t) = y1jk(φ,λ) + y2jk(φ,λ)f(φ,λ,t) + εjk(φ,λ,t), (1)
где f(φ,λ,t) - инвариант (эталон) суточного хода температуры воздуха.
Из (1) f(φ,λ,t) заменим многолетним суточным ходом за j-тый день u¯ j (φ,λ,t), откуда:
ujk(φ,λ,t)] = y1jk(φ,λ) + y2jk(φ,λ) (φ,λ,t) + εjk(φ,λ,t). (2)
Определим параметры y1jk(φ,λ), y2jk(φ,λ) методом наименьших квадратов и найдем за каждый день среднеквадратические отклонения y3jk(φ,λ) значений εjk(φ,λ,t). Изменение каждого параметра в течение года можно представить в виде (2), где в качестве u используются средние многолетние сезонные изменения параметра. В итоге получим шесть новых параметров, не зависящих от номера дня в году j. Далее поле каждого из шести параметров представим в виде линейной зависимости от его же среднего многолетнего поля. Коэффициенты этих зависимостей можно найти эмпирически. Они не зависят от пространственных координат, т.е. являются глобальными характеристиками сезонных изменений определенного процесса. Межгодовые изменения этих параметров являются полной характеристикой колебаний климата. Обозначим вектор параметров, характеризующий годовой ход как p¯, а вектор параметров, характеризующий флуктуации от годового хода как q¯. Изменение этих параметров происходит взаимообусловлено, т.е. можно записать
Ф(p¯ ,q¯ ) =ε¯ , (4)
где ε¯- неизвестный вектор погрешностей параметров p¯,q¯ .
Выражение (4) можно называть уравнением состояния климатической системы, связывающее ее микро- и макросостояния. Оно позволяет оптимально согласовать значения искомых параметров, минимизируя вектор ε¯.
Эмпирическим обоснованием альтернативного подхода является адекватность линейной зависимости реализаций полей и процессов от их средних многолетних значений. Это может быть лишь в том случае, если правильно выбран инвариант и выполняется условие подобия (правило «золотой пропорции»). Адекватность линейных зависимостей легко проверяется визуально. Дополнительным доказательством обоснования альтернативного подхода является зависимость инвариантов только от внешних условий. Покажем это на примере зависимости среднего многолетнего годового хода температуры воздуха, осредненной по полушарию y¯(ξ) от склонения Солнца δ(ξ), ξ - время в течение года.
Представим y¯(ξ) и δ(ξ) в виде сумм:
δ(ξ) = δ1(ξ) + δ2(ξ), (5)
y¯(ξ) =y¯ 1(ξ) +y¯2(ξ) , (6)
где y¯1(ξ), δ1(ξ) - гармонические функции времени, а y¯2(ξ), δ2(ξ) - отклонения от гармонических функций (от равномерного движения Земли по окружности).
Для построения эмпирической модели был использован архив "National Climatic Data Center. Data Set 9300 (DSI-9300). October 18, 2004. Asheville, NC 28801-5001 USA", содержащий минимальные и максимальные значения температуры воздуха за сутки. Из архива были выбраны данные по 435 станциям, примерно равномерно расположенных на Евро-Азиатском и Северо-Американском континентах за период 1948-1999 годы, не содержащий существенных пропусков. Анализ показал, что составляющая y¯1(ξ) зависит от δ1(ξ) с запаздыванием 30 дней. (Для сравнения: гармонические составляющие средней амплитуды годового хода на полушарии запаздывает на 22 дня, а гармоническая составляющая среднеквадратического отклонения от годового хода всего на 17 дней.)
Методом регрессионного анализа получено следующее выражение для y¯2(ξ):
y¯2(ξ) = с0 + с1x1 + c2x3 + c3x1x3 + (d0 + d1x2 + d2x3 +d3x1x2 + d4x1x3)x4 + ε, (7)
где: x1 = δ2(ξi), x2 = Δδ2(ξi), x3 = Δ2δ2(ξi), x4 = Δδ1(ξi), c0 - c3, d0 - d4 - параметры.
Коэффициент детерминации модели (7) равен 90.7%, все параметры статистически значимы не менее чем при 95% уровне вероятности. Результаты аппроксимации показывают, что модель (7) можно признать адекватной эмпирическим данным. Отметим, что аналогичные модели для амплитуды годового хода и среднеквадратического отклонения имеют коэффициенты детерминации близкие к 100%. Изменения климатических характеристик (параметров подобия глобальных полей), получаемых с помощью системного синтеза имеют две составляющие: «медленные» и «быстрые» колебания. Климатические характеристики, получаемые системным синтезом изменяются взаимосвязано (как и «положено» параметрам макросостояния системы). При этом взаимосвязь для медленных и быстрых изменений необходимо может быть, очевидно, разной, т.к. быстрые и медленные колебания отражают влияние на атмосферу разных компонентов климатической системы. Регрессионный анализ медленных изменений позволяет для всех параметров получить модели с коэффициентами детерминации 97% и более. Для короткопериодных колебаний коэффициенты детерминации заметно меньше - от 81% и выше. Это вполне понятно, т.к. погрешности климатических характеристик, обусловленные ограниченным объемом данных в пространстве и во времени, в основном, сосредоточены в короткопериодной составляющей.
Сравнение альтернативного и традиционного способа определения климатических характеристик. Основное принципиальное отличие двух подходов состоит в том, что в классической климатологии многолетнее поле никак не используется для оценки текущего климата. В системном синтезе, напротив, многолетнее поле используется как эталон для оценки климата за текущий год.
Применение климатических характеристик, полученных системным синтезом. Зная климатические параметры и многолетние поля, легко получить «нормы» на каждый текущий год, если в этом есть необходимость для решения прикладных задач. Высокая степень взаимосвязи «медленных» и «быстрых» колебаний климатических параметров может служить основанием для разработки эмпирических методов их прогнозирования. Это можно сделать, например, на основе теоремы Такенса [2] с помощью модели динамической авторегрессии, в которой случайная компонента может быть равной нулю (в отличии от статистической авторегрессии). При этом, модельные значения должны, очевидно, удовлетворять тем же соотношениям подобия полей, что и результаты наблюдений. Особенно важно, чтобы модельные поля были подобны эмпирическим инвариантам (эталонам) и траектория движения точки в фазовом пространстве переменных (4) была бы продолжением траектории, полученной по эмпирическим данным. В этом случае можно надеяться на то, что можно получить не только надежные значения климатических характеристик, но и детальный долгосрочный прогноз погоды.
Резюме. Любая теория должна опираться на эмпирические факты. Изложенный подход можно рассматривать как измерительно-вычислительный комплекс получения климатических характеристик, учитывающий современные достижения современной науки, в частности, теорию нелинейной динамики и самоорганизации. Надежные измерения климатических характеристик позволит плодотворно развивать и теоретические, и эмпирические методы долгосрочного прогнозирования погоды и климата.
Список литературы
- Колесников А.А. Прикладная синергетика: основы системного синтеза. - Таганрог: ТТИ ЮФУ, 2007. - 384 с.
- Малинецкий Г.Г., Потапов А.Б. Современные проблемы нелинейной динамики. М.: Едиториал УРСС. - 2002. - 360 с.
- Результаты исследований изменений климата для стратегий устойчивого развития Российской федерации. - М.: ООО ВИВА-СТАР. - 2005. - 212 с.
- Сороко Э.М. Золотые сечения, процессы самоорганизации и эволюции систем. - М.: КомКнига. - 2009. - 264 с.
- Степаненко С.Р., Воронцов А.А. Принципы системного синтеза эмпирических законов на примере моделирования климата. // Проблемы безопасности и чрезвычайных ситуаций. - 2009. - № 5. - с. 65-74.
- Физические основы климата и его моделирования - Л.: Гидрометеоиздат. - 1977. - 271 с.