Влияние глобального изменения климата на питание и ресурсы подземных вод

Содержание

Слайд 2

Содержание Почему меняется климат Прогнозы изменения климата – последние результаты Используемая

Содержание
Почему меняется климат
Прогнозы изменения климата – последние результаты
Используемая модель для

оценки влияния изменения климата на питание
Результаты моделирования на примере бассейна р. Жиздры ( в пределах Калужской области)
Выводы
Слайд 3

Прогнозы изменения климата – последние результаты - наблюдения (по данными: IPCC

Прогнозы изменения климата – последние результаты - наблюдения (по данными: IPCC Fourth

Assessment Report, 2007 Изменение климата, . Обобщенный доклад, доклад 2007 г.
Слайд 4

Радиационный баланс Средняя глобальная приповерхностная температура T зависит от баланса поступающих

Радиационный баланс

Средняя глобальная приповерхностная температура T зависит от баланса поступающих

потоков тепла Qin и уходящих потоков Qout
Поток Qin состоит из падающего на Землю потока тепла от солнца ( и эндогенного потока
Поток Qout это переизлучение поглощённой солнечной энергии в космическое пространство
Слайд 5

Радиационный баланс и катастрофические извержения вулканов В случае, если входящие и

Радиационный баланс и катастрофические извержения вулканов

В случае, если входящие и

выходящие потоки сбалансированы, глобальная температура остается постоянной
При изменении соотношения потоков меняется средняя температура
Каждое извержение вулканов приводит к выбросу в атмосферу аэрозольных частиц, которые поднимаются в стратосферу, уменьшая ее прозрачность для солнечной радиации, уменьшая тем самым поток тепла, приходящий в приземную атмосферу
Наблюдения показывают, что этот эффект заметен в течении года или немногим более
Слайд 6

Радиационный баланс и ядерная зима, в 80-х годах прошлого века в

Радиационный баланс и ядерная зима,

в 80-х годах прошлого века

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

Радиационный баланс и парниковый эффект Парниковый эффект – различная прозрачность атмосферы

Радиационный баланс и парниковый эффект

Парниковый эффект – различная прозрачность атмосферы для

входного и выходных потоков тепла
На диапазон длин волн 400—​1500 нм (видимый свет и ближний инфракрасный диапазон) приходится 75 % энергии солнечного излучения, большинство атмосферных газов не поглощают в этом диапазоне
Нагретая поверхность планеты и атмосфера излучают в дальнем инфракрасном диапазоне: так, 75 % теплового излучения приходится на диапазон 7,8—28 мкм
Атмосфера, содержащая газы, поглощающие в этой области спектра (т. н. парниковые газы. парниковые газы — H2OO, CO2, CH4) существенно непрозрачна для такого излучения, направленного от её поверхности в космическое пространство, Вследствие такой непрозрачности атмосфера становится хорошим теплоизолятором
Слайд 8

Радиационный баланс, парниковый эффект и загрязнение атмосферы Загрязнение атмосферы за счет

Радиационный баланс, парниковый эффект и загрязнение атмосферы

Загрязнение атмосферы за счет выбросов

в нее аэрозольных частиц от промышленности, транспорта и пр. приводит к эффекту, который называется потускнение атмосферы – он снижает ее прозрачность для входящего потока тепла от солнца
Парниковый эффект – снижает ее прозрачность для потока переизлучения
Компенсируют ли они друг друга?
По данным IPCC:
Потускнение – снижение доступной солнечной радиации на 0.7 Вт/кв.м
Рост концентраций парниковых газов уменьшает выходящий за пределы атмосферы поток переизлучения на 2.6 Вт/кв.м
Слайд 9

Изменение климата за последние 650 тыс лет Климат в течение тысячелетий

Изменение климата за последние 650 тыс лет

Климат в течение тысячелетий

меняется в связи с изменением потока тепла от Солнца за счет циклов Миланковича
Быстрый рост температуры в последние годы не может быть объяснен этими циклами
Слайд 10

Циклы Миланковича и климат Циклы Миланковича описывают периодически возникающие отклонения от

Циклы Миланковича и климат

Циклы Миланковича описывают периодически возникающие отклонения от средней

достигающего Земли количества солнечной радиации за большой период времени в пределах от 5 до 10 процентов. Причиной этих отклонений от средней интенсивности солнечного излучения на Земле являются три эффекта:
Прецессия: поворот земной оси с периодом около 25 750 лет, в результате которого меняется сезонная амплитуда интенсивности солнечного потока на северном и южном полушариях Земли;
Нутация: долгопериодические (так называемые вековые) колебания угла наклона земной оси к плоскости её орбиты с периодом около 41 000 лет;
долгопериодические колебания эксцентриситета орбиты Земли с периодом около 93 000 лет
Рост температуры в последние годы не объясняется циклами Миланковича
Слайд 11

Прогнозы изменения климата – последние результаты: факторы Наиболее значимыми среди факторов

Прогнозы изменения климата – последние результаты: факторы

Наиболее значимыми среди факторов

являются:
увеличение концентраций парниковых газов аэрозолей (находящихся в воздухе микроскопических частиц и капель) в атмосфере – парниковый эффект
колебания солнечной активности – поток солнечного тепла
оба фактора могут изменить радиационный баланс Земли и, соответственно, климат.
Измерение потока солнечного тепла при помощи спутников, показывает что его интенсивность практически не меняется за весь период наблюдений (порядка 50 лет), составляя 1367 Вт/кв.м
Слайд 12

Глобальные циркуляционные модели (ГЦМ) Для анализа процессов переноса тепла, влаги, воздушных

Глобальные циркуляционные модели (ГЦМ)

Для анализа процессов переноса тепла, влаги, воздушных масс

в атмосфере, тепла и масс воды в океане используются Глобальные циркуляционные модели
ГЦМ решает в планетарном масштабе на конечно-разностной сетке уравнения динамики атмосферы и океана, позволяя рассчитать температуру атмосферы и океана, выпадение осадков, направление ветров и течений в узлах конечно разностной сетки имеющих размер порядка 1-го градуса
Если ввести в ГЦМ динамику прозрачности атмосферы за счет эмиссии парниковых газов, то при помощи такой модели можно сделать прогноз – как во времени будет меняться температура и интенсивность выпадения осадков в различных регионах, а также в среднем для всей Земли
Слайд 13

Прогнозы изменения климата – последние результаты: эпигнозное моделирование при помощи ГЦМ

Прогнозы изменения климата – последние результаты: эпигнозное моделирование при помощи ГЦМ


Наблюдаемые (черная кривая) и моделируемые глобальные среднегодовые аномалии температуры (°С) относительно средних значений наблюдений за период с 1900 г. по 1930 г.
Контрольное моделирование -зеленая кривая без учета парникового эффекта
три независимых глобальных модели с использованием одного и того же фактора воздействия парникового газа плюс аэрозольные воздействия

Слайд 14

Прогнозы изменения климата – последние результаты: прогнозное моделирование Результаты простой модели:

Прогнозы изменения климата – последние результаты: прогнозное моделирование

Результаты простой модели: историческая

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

Возможное изменение суммарного стока с водосбора при наиболее вероятном сценарии изменения

Возможное изменение суммарного стока с водосбора при наиболее вероятном сценарии изменения

климата для территории Московского Артезианского Бассейна

Глобальные циркуляционные модели предсказывают, что наиболее вероятен на Европейской части Росси рост температуры на примерно 4 градуса за сто лет
Изменение осадков пересказываются существенно хуже. По сценарию B1 предсказывается их небольшой рост (5-15%)
Как изменится сток рек (поверхностный и подземный)?

Слайд 16

Возможное изменение суммарного стока с водосбора при наиболее вероятном сценарии изменения

Возможное изменение суммарного стока с водосбора при наиболее вероятном сценарии изменения

климата для территории МАБ

Для оценки изменения суммарного стока при изменении средней температуры и осадков использовалась зависимость Fu:
здесь E0 – испаряемость мм/год, E – суммарное испарение, P – годовая сумма осадков
Сток = Осадки - Испарение

Слайд 17

Возможное изменение суммарного стока с водосбора при наиболее вероятном сценарии изменения климата для территории МАБ

Возможное изменение суммарного стока с водосбора при наиболее вероятном сценарии изменения

климата для территории МАБ
Слайд 18

Возможное изменение суммарного стока с водосбора при наиболее вероятном сценарии изменения

Возможное изменение суммарного стока с водосбора при наиболее вероятном сценарии изменения

климата для территории МАБ -резюме

Глобальные циркуляционные модели предсказывают, что наиболее вероятен на Европейской части Росси рост температуры на примерно 4 градуса за сто лет
Изменение осадков пересказываются существенно хуже. По сценарию B1 пересказывается их небольшой рост (5-15%)
Следовательно существенного изменения общего не следует ожидать – увеличение испаряемости примерно на 80-100 мм возможно компенсируется увеличением осадков
Как изменится соотношение поверхностной и подземной части стока???

Слайд 19

Подземные воды и ландшафты – как изменение условий на поверхности земли скажется на подземных водах?

Подземные воды и ландшафты – как изменение условий на поверхности земли

скажется на подземных водах?
Слайд 20

Модель для оценки влияния изменения климата Модель состоит из трех расчетных

Модель для оценки влияния изменения климата

Модель состоит из трех расчетных блоков.


Первый блок – трансформация осадков на поверхности земли
Второй блок – влагоперенос в зоне аэрации
Третий блок – геофильтрация в зоне полного насыщения
Питание подземных вод в данной модели – расход потока через нижнюю границу зоны аэрации – интерфейс между вторым и третьим блоками
Слайд 21

Модель для оценки влияния изменения климата – первый и второй блоки

Модель для оценки влияния изменения климата – первый и второй блоки

Моделируемые

процессы, определяющие баланс влаги на поверхности земли в зоне аэрации:
Выпадение осадков, задержка растительностью, накопление снега, испарение с поверхности листьев, впитывание в почву
Формирование поверхностного стока с талой и мерзлой почвы
Эвапотранспирация влаги из корнеобитаемого слоя
Вертикальный влагоперенос в зоне аэрации от поверхности земли до ее нижней границы зоны
Слайд 22

Первый блок: Трансформация осадков на поверхности земли: моделируемые процессы Выпадение осадков,

Первый блок: Трансформация осадков на поверхности земли: моделируемые процессы

Выпадение осадков, их

накопление в виде снега, таяние, задержкой растительностью, формирование стока с учетом состояния почвы: талая или мерзлая
Расчет потенциальной эвапотранспирации,
Расчет испарения с поверхности листьев растительности, на которых задержались осадки
Результат – суточный объем влаги, который впитывается в почву с поверхности
Слайд 23

Процессы на поверхности земли – пример модели трансформации влаги Трансформация влаги

Процессы на поверхности земли – пример модели трансформации влаги

Трансформация влаги зимой:

Базовые переменные: Vs – слой снега, Vl – слой воды, Hs – мощность снега, Ts – температура снега
Слайд 24

Трансформация осадков на поверхности земли (входные и выходные данные) Входные данные

Трансформация осадков на поверхности земли (входные и выходные данные)

Входные данные ряды

суточных:
осадков O(t),
Минимальной и максимальной температуры T(t),
солнечной радиации R(t)
динамики индекса развития поверхности листьев LAI(t) растительности
Выходные данные ряды суточных объемов:
Влаги, впитывающихся в почву с поверхности V(t)
Потенциальной эвапотранспираци Ep(t)
Поверхностного стока R(t)
Слайд 25

Второй Блок: Влагоперенос в зоне аэрации (модель Hydrus 1-d) Одномерная вертикальная

Второй Блок: Влагоперенос в зоне аэрации (модель Hydrus 1-d)

Одномерная вертикальная

насыщенно-ненасыщенная фильтрация от поверхности земли до нижней границы, расположенной :
глубже зоны колебания уровня грунтовых вод при их близком залегании
На глубине зоны гравитационного стекания влаги (dH/dz=-1) при глубоком залегании уровня
В зоне аэрации выделяется корнеобитаемая зона из которой происходит отбор воды корнями растений
Слайд 26

Влагоперенос в зоне аэрации (входные и выходные данные) Входные данные ряды

Влагоперенос в зоне аэрации (входные и выходные данные)

Входные данные ряды

суточных объемов:
Влаги, впитывающихся в почву с поверхности V(t)
Потенциальной эвапотранспираци Ep(t) за вычетом испарившейся с поверхности листьев влаги
Выходные данные ряды суточных объемов:
Суммарного испарения из почвы E(t)
Стока через нижнюю границы модели – то есть питания подземных вод W(t)
Слайд 27

Третий Блок: региональная нестационарная модель потока подземных вод (MODFLOW) Модель потока

Третий Блок: региональная нестационарная модель потока подземных вод (MODFLOW)

Модель потока подземных

вод в его естественных границах
Входные данные: переменное во времени питание подземных вод, рассчитанное в предыдущем блоке (для каждого узла, или серии узлов, объединенных в зону с одинаковыми условиями питания)
Выходные – изменение уровней подземных вод и их разгрузки в реки
Слайд 28

Как использовать эту модель для оценки влияния изменения климата? Промоделировать процесс

Как использовать эту модель для оценки влияния изменения климата?

Промоделировать процесс при

стационарных входных рядах:
осадков O(t),
температур T(t), с
солнечной радиации R(t)
То же самое для измененных с учетом прогноза изменения климата рядах
Разница – влияние изменения климата на питание
Слайд 29

Как использовать эту модель для оценки влияния изменения климата? (продолжение) Для

Как использовать эту модель для оценки влияния изменения климата? (продолжение)

Для моделирования

входных рядов использовался генератор метеоданных:
LARS-WG: A Stochastic Weather Generator for Use in Climate Impact Studies (Semenov M. 2002)
LARS-WG позволяет:
Генерировать стационарную серию суточных метеоданных любой длины, “обучившись” на наблюденных метеоданных за несколько десятков лет
Вводить в генерируемые ряды сценарии изменения любого из метеорядов: осадков, мин. суточных температур, max. суточных температур, радиации
Слайд 30

Как использовать эту модель для оценки влияния изменения климата? (продолжение) Сценарии

Как использовать эту модель для оценки влияния изменения климата? (продолжение)

Сценарии вводятся

на месячном уровне, то есть задается, как изменяются среднемесячные параметры генерируемых рядов относительно исходных параметров рядов на которых проводилось “обучение” генератора
[USE OF A STOCHASTIC WEATHER GENERATOR IN THE DEVELOPMENT OF CLIMATE CHANGE SCENARIOS, by M. Semenov, E Barrow, Climatic Change 35: 397–414, 1997.c1997 Kluwer Academic Publishers. Printed in the Netherlands.]
Слайд 31

Применение модели к юго-западной части Московского Артезианского Бассейна (МАБ) – территория Калужской области

Применение модели к юго-западной части Московского Артезианского Бассейна (МАБ) – территория

Калужской области
Слайд 32

Моделирование : наблюденные тренды температуры по метеостанции г. Сухиничи за последние 30 лет

Моделирование : наблюденные тренды температуры по метеостанции г. Сухиничи за последние

30 лет
Слайд 33

Моделирование на примере бассейна Жиздры Калужской области Более чем 40 летний

Моделирование на примере бассейна Жиздры Калужской области

Более чем 40 летний ряд

суточных наблюдений по метеостанции Сухиничи применялся для “обучения”LARS WG моделированию стационарного ряда метеоданных в р-не г. Сухиничи
В качестве прогнозного сценария изменения климата использовалось следующее:
За будущие сто лет среднегодовая температура вырастет на 4.5 градусов
Летние температуры вырастут не более чем на 2 градуса, а основной рост температур будет в холодный период года
Солнечная радиация не изменится
Осадки увеличатся в зимний период примерно на 15% и сократятся в летний период примерно на 5%
Слайд 34

Моделирование на примере Калужской области: Базовый и измененный сценарии климата

Моделирование на примере Калужской области: Базовый и измененный сценарии климата

Слайд 35

Результаты моделирования питания Для расчета питания моделировалась зона аэрации, представленная в

Результаты моделирования питания

Для расчета питания моделировалась зона аэрации, представленная в основном

песчанистым разрезом (что характерно для флювиогляциальных отложений, развитых в районе г.Сухиничи)
Полагалось, что уровень грунтовых вод расположен достаточно глубоко на глубине 5м, поэтому на нижней границе зоны аэрации использовалось граничное условие постоянного напора
Слайд 36

Осредненные входные ряды Модельная температура Расчетная испаряемость

Осредненные входные ряды

Модельная температура

Расчетная испаряемость

Слайд 37

Результаты моделирования Высота снежного покрова Интенсивность снеготаяния

Результаты моделирования

Высота снежного покрова

Интенсивность снеготаяния

Слайд 38

Предвесенний запас влаги в снегу

Предвесенний запас влаги в снегу

Слайд 39

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

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

Слайд 40

Результаты моделирования: среднемноголетнее питание

Результаты моделирования: среднемноголетнее питание

Слайд 41

Результаты моделирования: временная изменчивость питания

Результаты моделирования: временная изменчивость питания

Слайд 42

Третий блок: региональная нестационарная модель потока подземных вод бассейна Жиздры Модель

Третий блок: региональная нестационарная модель потока подземных вод бассейна Жиздры

Модель бассейна

р. Жиздры вырезана из региональной модели юго-западной части МАБ (С.О. Гриневский, 2006)
Модель имеет размеры 60*50 км и включает 12 основных пластов и разделяющих толщ, представленных в регионе
Питание, разгрузка в реки, водоотбор – основные статьи баланса
Испарение подземных вод при близком УГВ не учитывается
Слайд 43

Региональная нестационарная модель потока подземных вод бассейна Жиздры – два сценария

Региональная нестационарная модель потока подземных вод бассейна Жиздры – два сценария

изменения среднегодового питания в 21 веке
Слайд 44

Питание подземных вод Среднемноголетняя карта питания подземных вод Питание подземных вод

Питание подземных вод

Среднемноголетняя карта питания подземных вод
Питание подземных вод каждый год

менялось в соответствии с питанием полученным на модели влагопереноса для
А)стационарного климата
Б) нестационарного климата
w*(t)=wsr(x,y)*K(t)
K(t)=Wmod(t)/WsrMod
Слайд 45

Результаты моделирования 100 летней динамики уровней подземных вод на региональной модели:

Результаты моделирования 100 летней динамики уровней подземных вод на региональной модели:

куст в долине р. Жиздры

Стационарный климат

Нестационарный климат

Слайд 46

Результаты моделирования 100 летней динамики уровней подземных вод на региональной модели:

Результаты моделирования 100 летней динамики уровней подземных вод на региональной модели:

куст в районе водозабора г. Сухиничи

Стационарный климат

Нестационарный климат

Слайд 47

Выводы Моделирование показало существенное влияние изменения климата на питание подземных вод,

Выводы

Моделирование показало существенное влияние изменения климата на питание подземных вод, приводящее

к его увеличению в исследуемом бассейне и возможному подъему уровней подземных вод
Основным механизмом определяющим увеличения питания в рассмотренной модели является сокращение поверхностного талого стока и, соответственно:
Инфильтрация в период зимних оттепелей