Результаты производственных испытаний технологии краткосрочных численных ледовых прогнозов для юго-западной части Охотского моря презентация

Охотского моря

В период 2008-2010 г.г. в рамках раздела 5.2 «Разработка, испытание и внедрение моделей, методов и технологий гидрометеорологических, в том числе ледовых прогнозов и расчетов различной заблаговременности для арктических и замерзающих морей России» ЦНТП Росгидромета «Научные исследования и разработки в области гидрометеорологии и мониторинга окружающей среды» в ФГБУ «ААНИИ» была разработана модель краткосрочного прогноза ледовых условий юго-западной части Охотского моря в зимний период. Тогда же, в 2009-2010 г.г., были выполнены серии авторских испытаний данной методики, показавшие удовлетворительное качество работы модели. Производственные испытания состоялись в январе-марте 2013 г. Испытания проводились на базе ФГБУ «ААНИИ».

С.В. Клячкин, З.М. Гудкович, Р.Б. Гузенко, Р.И. Май ФГБУ "ААНИИ"

время, Fa, Fh - потоки тепла через поверхность воды и через нижнюю границу однородного слоя соответственно, Mwa - изменение солености на границе вода-воздух (при отсутствии льда) или вода-лед (при наличии льда), Qsh - поток солей через нижнюю границу однородного слоя соответственно, cw - теплоемкость воды, ρw - плотность воды, hc - толщина квази-однородного слоя, μ - коэффициент горизонтального турбулентного обмена, Δ - оператор Лапласа.

Прогноз теплового состояния моря основан на уравнениях баланса тепла и солей, проинтегрированных по вертикали в пределах от поверхности до глубины верхнего квази-однородного слоя. При отсутствии льда происходит изменение солености, связанное с испарением. При наличии льда температура подледного слоя воды равна температуре замерзания и зависит только от солености. На нижней границе верхнего квази-однородного слоя происходит турбулентный подток тепла и солености из нижележащих слоев воды. Также учитывается горизонтальный перенос тепла и солености как адвективного, так и турбулентного характера.

Краткое описание модели

уравнения движения, гидростатики и неразрывности.
Модель позволяет рассчитывать баротропную составляющую поля течений, бароклинная составляющая задается как поле постоянных течений. Баротропная мода воспроизводится по проинтегрированным по вертикали уравнениям движения и неразрывности от поверхности моря до глубины квази-однородного слоя.

где Uh Uw - горизонтальные и вертикальные компоненты скорости течения соответственно; ρ - плотность, Р - давление; F - турбулентный (вертикальный и горизонтальный) обмен количеством движения, - вектор градиента.
Учет приливных явлений осуществляется путем в задании в каждой ячейке сетки гармонических постоянных главных приливных волн (M2, S2, N2, K2, K1, O1, P1, Q1) и вычисления уровня и составляющих течений.

где ht, Ut, Vt – уровень моря и компоненты приливных течений, f – редукционный множитель, А, Au, Av – амплитуды гармоник, q – угловая скорость, Vо+u – астрономическая часть фазы, g°, gu, gv – фазы гармоник, t - время.
Числовые значения амплитуд и фаз получаются с помощью модели [Padman, Erofeeva, 2004]

модели

где Но - начальная толщина льда, hs - толщина снега, λs - теплопроводность снега, Ta, Ts - температура воздуха и снежно-ледяной поверхности соответственно.
Толщина снега определяется по толщине льда с помощью эмпирических соотношений.

Температура верхней поверхности снега оценивается по уравнению теплового баланса.

где Θ - температура замерзания воды, Fa - поток тепла на границе снег - воздух.

где H – толщина льда, λл - теплопроводность льда, Tл - температура льда, k - удельная теплота плавления льда, ρл – плотность льда. Qw - поток тепла от воды к нижней поверхности льда.
В период нарастания льда уравнение теплового баланса решается относительно толщины льда:

движения в нестационарной постановке.
Тангенциальное напряжение на верхней поверхности льда определяется по скорости приземного ветра с помощью известного квадратичного закона.
Тангенциальное напряжение на нижней границе ледяного покрова определяется относительной скоростью дрейфа льда.
Реология ледяного покрова основана на вязко-пластической параметризации сил внутреннего взаимодействия. Внутренние напряжения σξξ, σχχ, соответствующие главным осям тензора скоростей деформаций, рассчитываются по формуле [Аппель, Гудкович, 1992]:

1. K = K0 H (3C - 2) при выполнении условий:

2. K = 0 при невыполнении хотя бы одного из этих условий
где W – скорость дрейфа льда, С – сплоченность льда, К0 = 1010 кг с-1 м-2

Если напряжения в ледяном покрове превышают некоторый предел, то это приводит к торошению. Этот предел оценивается как устойчивость ледяной пластины, лежащей на упругом основании, по отношению к продольному изгибу.
В противном случае поле дрейфа льда корректируется таким образом, чтобы не происходило накопление «лишнего» объема льда

мили), общая размерность матрицы 81×212 ячеек

характеризуется пространственными координатами в декартовой системе, толщиной, торосистостью, разрушенностью и скоростью.

Термическое изменение толщины рассчитывается отдельно для каждого маркера.

Расчет дрейфа льда основан на лагранжевом методе. Силы внутреннего взаимодействия, возникающие при дрейфе, определяются в соответствии с эйлеровским подходом и считаются одинаковыми в пределах ячейки.

Выторашиваемые маркеры исключаются из дальнейшего расчета, а их суммарная толщина равномерно распределяется среди оставшихся маркеров данной ячейки как прибавка толщины за счет торосов.

соответствии с международным кодом «СИГРИД-3»

Цифровой файл с фактическим данными о состоянии ледяного покрова в каждой ячейке расчетной сетки.

температуры воздуха в бинарном формате GRIB (ЕЦСПП). Временная дискретность 12 часов, пространственный шаг 2,5о.

С помощью специальной процедуры эти поля интерполируются в регулярную сетку модели.

в формате ГИС ArcView составлялись в ФГБУ «ААНИИ» 2 раза в неделю (по вторникам и пятницам). В эти же дни составлялись ледовые прогнозы заблаговременностью соответственно 3 и 4 суток.

Первый прогноз был составлен 29 января, последний – 22 марта.

Общее количество прогнозов – 16.

Порядок проведения испытаний

обозначены на фактических (исходных и контрольных) ледовых картах, а именно по общей и частной сплоченности.

Для оценки прогностического метода применяются две характеристики: оправдываемость и эффективность.
Оправдываемость прогноза – отношение количества ячеек сетки, в которых прогноз оказался успешным, к общему количеству ячеек сетки.
Эффективность прогноза – разность между обеспеченностью модельного прогноза и инерционного прогноза (инерционный прогноз – это неизменное состояние ледяного покрова в течение срока заблаговременности).

Предельно допустимая ошибка прогноза определяется стандартными градациями сплоченности льда (0, 1-3, 4-6, 7-8, 9-10 баллов).

Чтобы метод был признан удовлетворительным, необходимо, чтобы обеспеченность была достаточно высокой (формальных пороговых значений нет, но традиционно принимается величина 70-75 %), а эффективность – положительной. Если хотя бы одно из этих условий не выполняется, то прогностический метод следует признать неприемлемым.

Каждая последующая фактическая ледовая карта использовалась для оценки достоверности предыдущего прогноза.

Порядок оценки достоверности прогнозов

достоверности прогнозов

сплоченности

Оправдываемость прогнозов частной сплоченности составила в среднем 86,5% при эффективности +3,1%.

Систематическая ошибка не выявлена.

Результаты оценок достоверности прогнозов

на 15 февраля 2013 года

на 05 марта 2013 года

запись результатов) отработала без единого сбоя.

- Прогностическая методика не содержит заметной систематической ошибки.

- Порядка 80% ошибок прогнозов группируются в пределах ±1 балл, повторяемость грубых ошибок прогнозов общей сплоченности не превышает 1-2%, частной сплоченности - примерно 3-5%.

- Наиболее удачными были прогнозы частной сплоченности 1-летнего среднего льда (оправдываемость – более 92%, эффективность – более 3,5%), наименее – прогнозы молодого льда (оправдываемость – около 81%, эффективность 1,6%).

- Основные недостатки методики – (1) расчет горизонтальных потоков тепла в прикромочной области, который в наибольшей степени влияет на эволюцию молодых льдов, и (2) прогноз поля течений.

- В целом, уровень оправдываемости и эффективности прогнозов оказался достаточно удовлетворительным и близким к традиционным значениям, полученным на других регионах. Это дает основание заключить, что методика в целом адекватно отражает природный процесс и может быть представлена на рассмотрение ЦМКП Росгидромета.