湖泊科学   2023, Vol. 35 Issue (6): 2111-2122.  DOI: 10.18307/2023.0641
0

研究论文——流域水文与水资源安全

引用本文 [复制中英文]

梁新歌, 王涵, 赵爽, 宋春桥, 21世纪以来泛北极湖泊水位变化时空特征及原因探讨. 湖泊科学, 2023, 35(6): 2111-2122. DOI: 10.18307/2023.0641
[复制中文]
Liang Xinge, Wang Han, Zhao Shuang, Song Chunqiao. Spatial-temporal characteristics and causes of pan-Arctic lake water level changes since the 21st century. Journal of Lake Sciences, 2023, 35(6): 2111-2122. DOI: 10.18307/2023.0641
[复制英文]

基金项目

国家重点研发计划(2022YFF0711603)、中国科学院战略性先导科技专项(A类)(XDA23100102)和国家自然科学基金项目(41971403, 41801321)联合资助

通信作者

宋春桥, E-mail: cqsong@niglas.ac.cn

文章历史

2022-11-07 收稿
2023-03-14 收修改稿

码上扫一扫

21世纪以来泛北极湖泊水位变化时空特征及原因探讨
梁新歌1,3 , 王涵2 , 赵爽1,3 , 宋春桥3     
(1: 河南理工大学测绘与国土信息工程学院, 焦作 454000)
(2: 中国矿业大学环境与测绘学院, 徐州 221000)
(3: 中国科学院南京地理与湖泊研究所流域地理学重点实验室, 南京 210008)
摘要:在全球气候变暖和极端气候事件增加的背景下, 流域水文循环过程受到的影响越来越强烈, 导致湖泊水位变化表现出复杂的时空特征。而泛北极地区是地球上湖泊数量与面积分布最为集中的区域之一, 该地区湖泊对气候变化响应非常敏感。因此, 了解这些湖泊近期水文变化特征十分必要。本研究共搜集了36个泛北极大型湖泊(>500 km2)基于遥感或站点观测的近20年水位数据, 分析其时空变化特征。本文使用线性回归模型来估算湖泊水位的变化趋势, 进而利用皮尔逊相关分析了其主要水文影响变量和大气环流机制, 并运用Mann-Kendall突变检验法探讨了水位突变的原因。结果表明, 泛北极湖泊的水位整体上呈现不同程度上升(平均速率为0.013 m/a), 有23个(64%)湖泊的水位呈上升趋势; 研究湖泊中有10个通过90%统计显著性检验。其中, 水位上升速率最大的湖泊是位于哈萨克斯坦的腾吉兹湖, 上升速率为0.078 m/a。泛北极湖泊水位的波动主要与径流有关, 有19个(53%)湖泊的水位波动与径流的增加更为相关; 相比而言, 位于亚洲的极地湖泊水位的上升与流域蒸发的降低显著相关, 尤其是库苏古尔湖。从区域大气环流影响来看, 泛北极湖泊水位变化主要与厄尔尼诺-南方涛动有关, 其次是北极涛动和北大西洋涛动。本研究有助于加深对泛北极湖泊近20年水位变化规律及气候影响特征的科学理解。
关键词水位    泛北极    湖泊    气候变化    大气环流    
Spatial-temporal characteristics and causes of pan-Arctic lake water level changes since the 21st century
Liang Xinge1,3 , Wang Han2 , Zhao Shuang1,3 , Song Chunqiao3     
(1: College of Surveying and Geotechnical Engineering, Henan Polytechnic University, Jiaozuo 454000, P.R. China)
(2: School of Environment and Spatial Informatics, China University of Mining and Technology, Xuzhou 221000, P.R. China)
(3: Key Laboratory of Watershed Geographic Sciences, Nanjing Institute of Geography and Limnology, Chinese Academy of Sciences, Nanjing 210008, P.R. China)
Abstract: In the context of climate warming and intensified extreme climate events, the hydrological cycle and processes are increasingly affected, leading to certain spatial and temporal changes in hydrological elements. The pan-Arctic region is one of the most concentrated areas of lake distribution in both count and area on the Earth, and the lakes in this region are very sensitive to climate change. It is essential to understand the characteristics of recent hydrological changes of these lakes. In this study, the water level data for a total of 36 pan-Arctic large lakes (>500 km2) were collected over the past 20 years based on remote sensing or station observations, and their spatio-temporal variation characteristics were analyzed. We used a linear regression model to estimate the changing trend of lake water level, then used the Pearson correlation to analyze its main hydrological variables and atmospheric circulation factors. The Mann-Kendall mutation test was applied to explore the possible reasons for the abrupt changes in lake water level. The results indicated that the water levels of the pan-Arctic lakes showed different degrees of increase overall (mean rate at 0.013 m/a), with 23 lakes (64%) showing an increase in water levels. Only 10 lakes were featured by the linear trends at the 90% statistical significance. Lake Tengiz in Kazakhstan ranked the fastest among these rising lakes, with a change rate of 0.078 m/a. The changes in water level of these pan-Arctic lakes were mainly related to runoff variations, with 19 lakes (53%) showing correlations with increases in runoff. In comparison, increased water levels of the lakes located in Asia were significantly correlated with decreased evaporation, especially Lake Khovsgol. The pan-Arctic lake level changes were mainly associated with El Nino Southern Oscillation, followed by Arctic Oscillation and North Atlantic Oscillation. This study is expected to deepen our scientific understanding on the water level variation patterns and climate impact on the pan-Arctic lakes over the last 20 years.
Keywords: Water level    pan-Arctic    lake    climate change    atmospheric circulation    

近年来,气候变化及其影响受到各界学者的广泛关注。气候变化加剧水循环的过程,驱动蒸发、降水、径流等水文要素发生变化,增强极端事件的发生频率,改变区域水量收支平衡[1-2]。2021年政府间气候变化专门委员会(IPCC)发布的第六次评估报告(AR6)指出,21世纪初全球气温将会升高1.5℃[3]。而北极是地球上变暖速度最快的地区,假设在最极端的情况下,北极年平均气温预计在本世纪末将会升高10℃[3-4]。温度升高会导致北极永久冻土的退化、冰川消融,影响该地区的水文过程和能量平衡[5]。因此,在全球气候变暖的背景下,针对(泛)北极湖泊的研究具有重要的意义。

湖泊作为地球水圈的重要组成部分,储存着地球表面最主要的液态淡水,对全球水循环具有重要影响[6]。湖泊为人类和生态系统提供了重要的水源,在维持生物多样性、人类农业生产和各种工业设施运营中发挥着重要作用[7]。水位是湖泊水量平衡的重要指标之一,对人类干预和气候变化影响非常敏感,水位的波动对湖泊生态系统也产生了显著的影响[8]。湖泊水位的极端波动可能会导致洪旱灾害、水资源短缺和湖泊水质变化等问题,对生态系统、渔业和旅游业都会产生潜在的影响[9]。北极湖泊数量众多,分布范围广泛,是北极地表景观的重要组成部分,为其栖居的动物和人类提供着重要的水源和食物[10]。有研究表明,近几十年来,在全球气候变化影响下,北极冻土带消退,水体流失,导致北极一些湖泊发生显著变化[11]。北极湖泊受气候影响非常显著,开展该区域湖泊水文变化以及驱动机制的研究,对科学认识北极响应全球气候变化特征具有重要意义。

目前,基于遥感手段开展湖泊水位变化监测的技术已经较为成熟,相关技术也在泛北极湖泊上得到了广泛地应用,但是系统地调查该地区湖泊水位变化特征及影响湖泊水位的主要因素(例如,气候和大气环流等)仍存在不足。本研究以泛北极(50°N以北)36个大型湖泊(面积>500 km2)为研究对象,分析湖泊21世纪以来水位的变化特征,结合水文气象变量和大气环流因子,分析湖泊水位与水文和大气环流的潜在联系。本研究旨在促进对泛北极湖泊水文变化的认识,为今后科学认识该地区湖泊变化对气候和大气环流的响应关系提供参考。

1 研究区概况

北极地区的气候终年寒冷,极昼与极夜交替出现[12]。该地区以寒带气候和温带大陆性气候为主,对气候变化响应非常敏感。北极的水资源在不同季节之间有着巨大的变化,春季水资源充裕,而冬季水资源稀少[13]。北极年平均气温接近0℃,夏季最高气温接近10℃[14]。大多数北极湖泊全年或季节性地被冰覆盖,破冰期主要发生在每年的5—7月,受人类活动的影响较小[15-16]。近几十年来,随着北极气温的升高,大部分高纬度湖泊都有延缓冰期和提前破冰的趋势[17]。本研究综合考虑湖泊的大小、分布地理位置和水位数据的可获取性,选择北纬50°以北的泛北极36个面积大于500 km2的湖泊为研究对象(图 1),分别为:(1)阿特托林湖;(2)大熊湖;(3)霍塔湖;(4)大奴湖;(5)布法罗湖;(6)克莱尔湖;(7)阿萨巴斯卡湖;(8)诺纳科湖;(9)小奴湖;(10)多雷湖;(11)马尼托巴湖;(12)多芬湖;(13)温尼伯戈西斯湖;(14)温尼伯湖;(15)穆斯湖;(16)锡皮韦斯克湖;(17)南印第安湖;(18)赖恩迪尔湖;(19)卡斯巴湖;(20)纽埃尔延湖;(21)希尼克湖;(22)图利马卢湖;(23)贝克湖;(24)阿马朱瓦克湖;(25)维纳恩湖;(26)韦特恩湖;(27)佩普西湖;(28)伊尔门湖;(29)拉多加湖;(30)奥涅加湖;(31)伊纳里湖;(32)钱尼湖;(33)贝加尔湖;(34)滕吉兹湖;(35)乌布苏湖;(36)库苏古尔湖。湖泊的具体信息见附表Ⅰ。本研究的湖泊对象占泛北极大型湖泊(>500 km2)数量的32%,面积和水量占比约为70%和92%,该统计信息来自于数据集HydroLAKES[18]

图 1 研究区湖泊分布 Fig.1 Map of the distribution of lakes in the study area
2 数据与方法 2.1 数据来源 2.1.1 湖泊数据

湖泊边界数据来自Messager等[18]在2016年发布的全球湖泊数据集HydroLAKES(https://www.hydrosheds.org/products/hydrolakes),此数据包含面积大于0.1 km2约143万个天然湖泊和人造水库的矢量边界。此外,该数据还提供了湖泊的面积、水量、经纬度等属性信息供人参考,本研究所选的湖泊信息均来自HydroLAKES数据集。湖泊流域数据来源于HydroBASINS数据集(https://www.hydrosheds.org/products/hydrobasins),它表示一系列包含子流域边界的全球矢量化多边形图层,包括100万个单独的子流域边界,并提供12个分层的子流域边界的矢量化数据[19]。本研究选择能包含湖泊的最小等级的流域作为此湖泊的流域。

2.1.2 水位数据

本研究收集的水位数据,主要来源于水文站、测高卫星和文献,具体信息见附表Ⅰ。本研究使用的水位数据都是根据日或月平均水位计算平均值得到的年水位。为了获得较长的水位时间序列,部分湖泊整合了多源的水位数据,具体见附表Ⅰ

水文站数据源主要包括两个水文信息网站:加拿大水资源调查局(WSC)和瑞典气象水文研究所(SMHI)。WSC提供的是加拿大水文站的实测水位数据,它包含大量水文站点实时和历史的水位数据,还可以通过各个水文站点来获得湖泊或者水库的日平均值、月平均值、年最大值、最小值和瞬时值的水位[20],数据可从https://wateroffice.ec.gc.ca/站点获取。SMHI包含大约200个水文站点的数据,它可以提供水位、径流、雪密度、雪深等多种数据产品。水位产品主要包括维纳恩湖和韦特恩湖这两个湖泊的历史和实时的月度和年度的水位最大值、最小值和平均值[21],可在http://vattenwebb.smhi.se/查看。

测高卫星数据主要来自HYDROWEB和DAHITI,水位信息经过严格的质量控制和精度评价。HYDROWEB(https://hydroweb.theia-land.fr/)提供的水位数据是基于Topex/Poseidon、Jason-1、Jason-2、ENVISAT和GFO这5颗卫星数据整合得到的,它可以免费提供从多卫星高度测量推导出的大约150个湖泊和水库的月水位变化[22]。DAHITI由德国理工大学(DGFI-TUM)开发,利用多任务卫星测高数据来获取内陆水域的水位时间序列,可在http://dahiti.dgfi.tum.de.免费获取[23]

2.1.3 ERA5-Land气象数据

ERA5-Land来源于欧洲中期天气预报中心(European Centre for Medium-Range Weather Forecasts,ECMWF),它是一个高分辨率(9 km)的再分析数据集,提供了从1950年1月到现在的50个陆地上的水和能量循环变量的全球0.1°×0.1°网格化月度数据,可从https://cds.climate.copernicus.eu/cdsapp#!/dataset/reanalysis-era5-land-monthly-means?tab=overview上获取[24]。本文使用的是2000—2020年的月度数据,选择了7个与湖泊水位有关的变量,包括总蒸发量(E)、潜在蒸发量(PE)、雪覆盖(SC)、径流(R)、雪水当量(SDWE)、总降水(TP)和温度(T),并以湖泊的流域边界为掩膜提取湖泊的变量值。

2.1.4 大气环流指数数据

参考以往学者报道的对极地气候影响比较显著的大气环流异常情势[25-28],本研究选取北极涛动(AO)、厄尔尼诺-南方涛动(ENSO)和北大西洋涛动(NAO)3个大气环流因子来分析湖泊水位变化与大气环流之间可能存在的联系。北极涛动是北半球温带大气变率的主要模式,主导着20°N以北的大气压力、风、温度和降水的大规模变化[29]。厄尔尼诺-南方涛动是由热带太平洋与大气之间复杂相互作用引起的气候准周期性变化[30]。北大西洋涛动是决定北半球中高纬度的主要大气环流类型,与北半球中高纬度的降水、气温等有着紧密的联系[31-32]。各个大气环流因子都是以指数的形式表示,此数据来源于美国国家气象局气候预报中心(https://psl.noaa.gov/gcos_wgsp/Timeseries/)。

2.2 方法 2.2.1 水位数据同化分析

一元线性回归方程是利用最小二乘函数对一个自变量和因变量之间关系的一种回归分析,公式如下:

$ \hat{y}=a_0+a_1 \hat{x} $ (1)

式中,$\hat{y}$表示最终水位值,$\hat{x}$表示测高水位值,a0a1为待求系数。对水位数据同化分析时,使用具有相同年份水位的水文站水位和测高水位数据,根据这两列值得到一元线性回归方程。随后将测高数据的年份水位值作为自变量输入到模型中,得到完整的水位时间序列。

2.2.2 皮尔逊相关分析

皮尔逊相关系数(R)表示的是两组变量之间相关性,它还可以表示这种相关性的方向和强度[33-34]R取值范围在-1到+1之间,R的绝对值越接近1,说明两个变量之间的相关性越强。如果R>0,就说明两个变量的关系为正相关;反之表明两个变量的关系为负相关。R的计算公式如下:

$ R=\frac{\sum\limits_{i=1}^n\left(x_i-\bar{x}\right)\left(y_i-\bar{y}\right)}{\sqrt{\sum\limits_{i=1}^n\left(x_i-\bar{x}\right)^2 \sum\limits_{i=1}^n\left(y_i-\bar{y}\right)^2}} $ (2)

式中,xi表示变量x的观测值,x表示变量x的平均值,yi表示变量y的观测值,y表示变量y的平均值。

本文通过对湖泊水位与水文变量和大气环流因子进行相关性分析,确定它们之间的潜在联系,通过相关系数的大小判断它们之间的相关性强弱。

2.2.3 Mann-Kendall突变检验

目前,许多研究人员使用Mann-Kendall方法检测时间序列的突变问题[35-36]。本文主要使用Mann-Kendall方法检测水位、水文气象变量和大气环流因子时间序列的突变。

对具有n个样本大小的时间序列x构造一个秩序列:

$ s_k=\sum\limits_{i=1}^k r_i, r_i=\left\{\begin{array}{l} 1, x_i>x_j(j=1, 2, \cdots, i) \\ 0, \text { else } \end{array}\right. $ (3)

假设时间序列x是随机且相互独立,定义以下统计量:

$ U F_k=\frac{\left[s_k-E\left(s_k\right)\right]}{\sqrt{\operatorname{Var}\left(s_k\right)}}(k=1, 2, \cdots, n) $ (4)

式中,UF1=0,E(sk)和Var(sk)分别是sk的平均值和方差。当序列x1x2,…,xnn个变量相互独立并且具有相同的连续分布时,E(sk)和Var(sk)可以使用以下公式计算:

$ E\left(s_k\right)=\frac{n(n-1)}{4} $ (5)
$ \operatorname{Var}\left(s_k\right)=\frac{n(n-1)(2 n+5)}{72} $ (6)

式中,UFk遵循标准正态分布,且给定的显著性水平a为0.05。

根据时间序列x的倒序,即xnxn-1,…,x1,重复上述过程,并使UBk=-UFk(k=n, n-1, …, 1)且UB1=0,则可以得到UBk(k=1, 2, …, n)。如果UFkUBk两条曲线相交,并且交点在两条临界线之间,则该交点即为突变开始的时间。如果UFkUBk两条曲线没有交点或者交点超出临界线,则说明突变趋势不显著。本研究通过湖泊水位与水文气象变量或大气环流因子在同一时间发生突变确定为该气象变量或该大气环流因子与该湖泊的突变有关。

3 结果与分析 3.1 湖泊水位变化

36个泛北极湖泊水位在21世纪初的20年期间变化情况如图 2所示,图 3展示了泛北极湖泊水位变化速率的空间分布。其中,滕吉兹湖(0.078 m/a)水位上升最快,平均振幅(4.218 m)也最大;图利马卢湖(-0.016 m/a)和库苏古尔湖(-0.015 m/a)的水位下降得最快。泛北极36个湖泊水位总体上呈上升趋势(平均速率为0.013 m/a),其中23个湖泊的水位以平均速率为0.024 m/a的趋势上升,13个湖泊水位以平均速率为-0.007 m/a的趋势下降。但是只有10个湖泊(如阿特托林湖、伊纳里湖、钱尼湖等)的水位在90%统计显著性水平下呈现显著上升或者下降趋势。水位上升湖泊总面积为126484.250 km2,水位下降湖泊总面积为79326.040 km2,水位上升湖泊的面积占总面积的61.46%,此面积统计数据来自HYRROLAKES数据集[18]。泛北极湖泊的年水位平均波动幅度为1.066 m,有20个(56%)的湖泊水位振幅小于1 m,水位振幅较大的几个湖泊分别为滕吉兹湖(4.218 m)、锡皮韦斯克湖(2.485 m)和阿萨巴斯卡湖(2.286 m)。

图 2 36个湖泊水位时间序列 Fig.2 Time series of water level of 36 lakes
图 3 湖泊年际水位变化速率(黑框表示该湖泊的水位没有通过90%统计显著性水平) Fig.3 Interannual water level variation rate map of lakes (The black box indicates that the water level of the lake did not pass the 90% significance level)

在位于北美洲的24个湖泊中,有16个湖泊的水位呈上升趋势,8个湖泊的水位下降,而水位下降的湖泊多集中于60°N以北的地区。北美洲湖泊水位的平均振幅为1.043 m,大部分湖泊(例如湖泊7~22)的水位振幅均超过1 m。位于欧洲的9个湖泊中,其中有6个湖泊的水位呈上升趋势,3个湖泊水位下降。在50°N以北的亚洲区域只研究了3个湖泊,其中滕吉兹湖水位上升速率远超于其他两个湖泊的下降速率。欧亚大陆极区12个湖泊水位的平均振幅为1.113 m,但其中只有4个湖泊的水位振幅超过1 m,其他几个湖泊的振幅均小于1 m。

3.2 湖泊水位变化与水文气象变量的关系

湖泊月水位变化与对应月份水文变量的相关性如图 4所示。泛北极湖泊水位变化主要受到流域径流的影响,有19个湖泊的水位与径流最为相关。其中,有9个湖泊水位与径流表现为强相关(相关系数大于0.6),而多芬湖的水位与径流的相关性最为显著(相关系数大于0.8)。

图 4 湖泊水位与水文变量的关系(P < 0.01) Fig.4 Relationship between lake water levels and hydrological variables (P < 0.01)

位于北美洲和欧洲的湖泊水位变化主要与径流有关,相较而言,位于亚洲的极地湖泊水位变化较大程度受到温度和降水的影响。其中,阿特托林湖的水位与其流域内的蒸发、雪水当量和雪覆盖相关性极强,与蒸发的相关系数更是高达0.90。布法罗湖、克莱尔湖和阿萨巴斯卡湖的水位均与降水和温度的相关性更强,图利马卢湖、贝克湖和阿马朱瓦克湖与雪水当量、雪覆盖和温度表现得相关性更强。近年来,由于大气变暖,一些高纬度地区的冰川正在快速融化,永久冻土融化,尤其是加拿大和俄罗斯北极地区的径流在稳步增加,从而导致湖泊水位上升[37-39]

3.3 湖泊水位变化与大气环流指数的关系

近年来,许多研究发现AO、ENSO和NAO对北半球中高纬度地区的气候变化有着显著的影响[29-31]。为了探讨湖泊水位年际波动与大气环流的关系,本研究计算了每个大气环流指数年值序列的3年滑动平均值与水位时间序列的相关性,分析结果如图 5所示。泛北极湖泊水位变化主要与厄尔尼诺-南方涛动有关,其次是北极涛动和北大西洋涛动。其中,11个湖泊的水位与厄尔尼诺-南方涛动中等程度相关(相关系数大于0.4),7个湖泊的水位与北极涛动中等程度相关,6个湖泊的水位与北大西洋涛动中等程度相关。有研究发现,最近随着北极涛动的变化,北极的海平面压力下降,夏季海冰面积逐渐减少,湖泊的水位也跟着发生变化[40]

图 5 湖泊水位与大气环流指数的关系(P < 0.01) Fig.5 Correlations between lake water levels and atmospheric circulation index (P < 0.01)

北美洲湖泊水位波动与北大西洋涛动和北极涛动较为相关。其中,诺纳科湖水位与北极涛动表现为强相关,与北大西洋涛动和厄尔尼诺-南方涛动表现为中等程度相关;穆斯湖水位与北大西洋涛动表现为强相关,与北极涛动表现为中等程度相关;赖恩迪尔湖水位与厄尔尼诺-南方涛动的相关性强。欧亚大陆湖泊水位变化主要受北极涛动和厄尔尼诺-南方涛动影响:9个湖泊与北极涛动负相关,8个湖泊与厄尔尼诺-南方涛动负相关,4个湖泊与北大西洋涛动负相关。其中,佩普西湖的水位与北极涛动具有强相关,与厄尔尼诺-南方涛动和北大西洋涛动中等程度相关,奥涅加湖的水位与厄尔尼诺-南方涛动表现为强相关。

4 讨论

为探讨水文气象条件和大气环流对湖泊水位时间序列突变的潜在影响,研究分别对湖泊水位、水文气象变量和大气环流指数时间序列进行了突变检验,初步分析了引起水位发生突变的原因。值得注意的是,在所有湖泊中,穆斯湖水位的突变并不明显,没有探测到水位的突变年份,表 1展示了湖泊年水位的突变检验结果。

表 1 研究区36个湖泊年水位突变年份 Tab. 1 Year of abrupt change in annual water level of 36 lakes in the study area

图 6展示的是湖泊水位突变与水文气象变量的关系。泛北极湖泊水位突变主要受到潜在蒸发量、总蒸发量、温度和降水的影响。其中,有8个湖泊水位突变与潜在蒸发量有关,有6个湖泊水位突变受到总蒸发量和温度的影响。其中,多雷湖、大奴湖、卡斯巴湖、库苏古尔湖、纽埃尔延湖、赖恩迪尔湖和南印第安湖水位的突变并没有受到气象要素的影响;比如,位于北美洲的几个湖泊(阿马朱瓦克湖、贝克湖、图利马卢湖等)水位的突变与流域内蒸发量和温度更相关。而位于欧亚大陆的湖泊(如佩普西湖、韦特恩湖、维纳恩湖等)水位的突变主要受到潜在蒸发量、温度和降水的影响。而位于60°N以北的霍塔湖和诺纳科湖主要受到了雪水当量的影响。

图 6 湖泊水位突变与水文气象变量的关系:(a)总蒸发量、(b)潜在蒸发量、(c)径流、(d)雪水当量、(e)温度、(f)总降水 Fig.6 Relationship between abrupt changes in lake level and hydrometeorological variables: (a) total evaporation, (b) potential evaporation, (c) runoff, (d) snow depth water equivalent, (e) temperature, (f) total precipitation

图 7展示的是湖泊水位突变与大气环流的可能联系。北极涛动和北大西洋涛动均是在2014年发生的突变,而厄尔尼诺-南方涛动在多年(如2006年、2015年、2017—2019年)发生了突变。泛北极湖泊水位突变主要与厄尔尼诺-南方涛动有关,其次是北极涛动和北大西洋涛动。值得注意的是,阿特托林湖、大熊湖、库苏古尔湖等这几个湖泊并没有受到大气环流的显著影响。

图 7 湖泊水位突变与大气环流指数的关系:(a)北极涛动和北大西洋涛动,(b)厄尔尼诺-南方涛动 Fig.7 The relationship between abrupt changes in lake levels and atmospheric circulation index: (a) Arctic Oscillation and North Atlantic Oscillation, (b) El Nino Southern Oscillation

ENSO主要影响了中纬度北太平洋中部以及热带东太平洋的蒸发量,导致水位发生变化。在拉尼娜事件中,水汽在北太平洋中部异常辐散(辐合),大气的水汽含量减小(增大),造成蒸发量的增大(减小)。而在热带东太平洋由于降水量的增大(减小),导致蒸发量的增大(减小)[41]。在北大西洋上空,AO的型态与NAO高度相关,关系最为密切[42]。在北半球的冬季,AO/NAO会影响欧洲的气温及降水,从而使湖泊的水位发生变化。当AO/NAO处于正相位时,北大西洋地区高纬的气旋性环流异常加强,将暖湿空气带到北欧大陆,使欧洲气温异常偏高,降水偏多[43]

5 结论

本研究利用多源水位数据,包括从水文站、测高卫星以及文献中获取的水位资料,分析了泛北极36个湖泊近20年水位的变化特征,进而分析其变化规律与水文气象因素和大气环流的潜在联系,还进一步探讨了湖泊水位突变特征及影响因素。本研究调查的泛北极大型湖泊在21世纪以来水位整体上呈现上升趋势,水位上升湖泊的数量和面积均高于水位下降湖泊,但是只有10个湖泊通过了90%统计显著性检验。泛北极湖泊的年际水位变化振幅平均为1.066 m,水位振幅最大的湖泊为腾吉兹湖(4.218 m)。

径流是影响泛北极湖泊水位发生变化的主要因子,有19个湖泊的水位与径流最为相关。北美洲和欧洲的湖泊水位主要与径流更为相关,而亚洲的湖泊水位的上升则更大程度受到蒸发增加的影响。泛北极湖泊水位变化主要受到了厄尔尼诺-南方涛动的影响,其次是北极涛动和北大西洋涛动。位于北美洲的湖泊水位波动与北大西洋涛动和北极涛动较为相关,而位于欧亚大陆湖泊水位变化主要受北极涛动和厄尔尼诺-南方涛动主导。泛北极湖泊水位的突变主要与流域蒸发、温度和降水有关,其中位于北美洲湖泊水位的突变主要与蒸发和温度有关,而位于欧亚大陆湖泊水位的突变主要受到潜在蒸发量、降水和温度的影响。

6 附录

附表Ⅰ见电子版(DOI: 10.18307/2023.0641)。

附表Ⅰ 研究区36个湖泊的具体信息 Appendix Ⅰ Detailed information on the 36 lakes in the study area
7 参考文献

[1]
郭靖. 气候变化对流域水循环和水资源影响的研究. 武汉: 武汉大学, 2010.
[2]
Yin JB, Guo SL, Yang Y et al. Projection of droughts and their socioeconomic exposures based on terrestrial water storage anomaly over China. Science China Earth Sciences, 2022, 52(10): 2061-2076. [尹家波, 郭生练, 杨妍等. 基于陆地水储量异常预估中国干旱及其社会经济暴露度. 中国科学: 地球科学, 2022, 52(10): 2061-2076. DOI:10.1360/SSTe-2021-0276]
[3]
Masson-Delmotte V, Zhai P, Pirani A et al. Climate change 2021: The physical science basis, 2021: 2. DOI: 10.1017/9781009157896.
[4]
Hollister KV, Thomas EK, Raynolds MK et al. Aquatic and terrestrial plant contributions to sedimentary plant waxes in a modern Arctic Lake setting. Journal of Geophysical Research: Biogeosciences, 2022, 127(8): e2022JG006903. DOI:10.1029/2022JG006903
[5]
In't Zandt MH, Liebner S, Welte CU. Roles of thermokarst lakes in a warming world. Trends in Microbiology, 2020, 28(9): 769-779. DOI:10.1016/j.tim.2020.04.002
[6]
Sheng Y, Smith L, Li J et al. A Pan-Arctic assessment of high-latitude lake change ~25 years apart. AGU Fall Meeting Abstracts, 2011, GC31C-01.
[7]
Song K, Wang Q, Liu G et al. A unified model for high resolution mapping of global lake (>1 ha) clarity using Landsat imagery data. Science of the Total Environment, 2022, 810: 151188. DOI:10.1016/j.scitotenv.2021.151188
[8]
Zohary T, Ostrovsky I. Ecological impacts of excessive water level fluctuations in stratified freshwater lakes. Inland Waters, 2011, 1(1): 47-59. DOI:10.5268/iw-1.1.406
[9]
Leira M, Cantonati M. Effects of water-level fluctuations on lakes: An annotated bibliography. In: Ecological Effects of Water-Level Fluctuations in Lakes. Dordrecht: Springer Netherlands, 2008: 171-184. DOI: 10.1007/978-1-4020-9192-6_16.
[10]
Ayala-Borda P, Lovejoy C, Power M et al. Evidence of eutrophication in Arctic lakes. Arctic Science, 2021, 7(4): 859-871. DOI:10.1139/as-2020-0033
[11]
Finger-Higgens R. Diminishing Arctic lakes. Nature Climate Change, 2022, 12(9): 782-783. DOI:10.1038/s41558-022-01466-7
[12]
Hassol SJ. Arctic climate impact assessment. Cambridge: Cambridge University Press, 2004.
[13]
Stuefer SL, Kane DL. Snow retention for increased water supply of shallow Arctic lakes. Cold Regions Science and Technology, 2016, 123: 32-43. DOI:10.1016/j.coldregions.2015.11.011
[14]
Thompson R, Ventura M, Camarero L. On the climate and weather of mountain and sub-Arctic lakes in Europe and their susceptibility to future climate change. Freshwater Biology, 2009, 54(12): 2433-2451. DOI:10.1111/j.1365-2427.2009.02236.x
[15]
Liu KS, Yao TD, Pearce DA et al. Bacteria in the lakes of the Tibetan Plateau and polar regions. Science of the Total Environment, 2021, 754: 142248. DOI:10.1016/j.scitotenv.2020.142248
[16]
Schütte UME, Cadieux SB, Hemmerich C et al. Unanticipated geochemical and microbial community structure under seasonal ice cover in a dilute, dimictic Arctic lake. Frontiers in Microbiology, 2016, 7: 1035. DOI:10.3389/fmicb.2016.01035
[17]
Lei RB, Leppäranta M, Cheng B et al. Changes in ice-season characteristics of a European Arctic lake from 1964 to 2008. Climatic Change, 2012, 115(3): 725-739. DOI:10.1007/s10584-012-0489-2
[18]
Messager ML, Lehner B, Grill G et al. Estimating the volume and age of water stored in global lakes using a geo-statistical approach. Nature Communications, 2016, 7: 13603. DOI:10.1038/ncomms13603
[19]
Lehner B, Grill G. Global river hydrography and network routing: Baseline data and new approaches to study the world's large river systems. Hydrological Processes, 2013, 27(15): 2171-2186. DOI:10.1002/hyp.9740
[20]
Pelletier PM. A review of techniques used by Canada and other northern countries for measurement and computation of streamflow under ice conditions. Hydrology Research, 1990, 21(4/5): 317-340. DOI:10.2166/nh.1990.0023
[21]
Bergström S. Principles and confidence in hydrological modelling. Hydrology Research, 1991, 22(2): 123-136. DOI:10.2166/nh.1991.0009
[22]
Crétaux JF, Arsen A, Calmant S et al. SOLS: A lake database to monitor in the near real time water level and storage variations from remote sensing data. Advances in Space Research, 2011, 47(9): 1497-1507. DOI:10.1016/j.asr.2011.01.004
[23]
Schwatke C, Dettmering D, Bosch W et al. DAHITI—an innovative approach for estimating water level time series over inland waters using multi-mission satellite altimetry. Hydrology and Earth System Sciences, 2015, 19(10): 4345-4364. DOI:10.5194/hess-19-4345-2015
[24]
Muñoz-Sabater J, Dutra E, Agustí-Panareda A et al. ERA5-Land: A state-of-the-art global reanalysis dataset for land applications. Earth System Science Data, 2021, 13(9): 4349-4383. DOI:10.5194/essd-13-4349-2021
[25]
Lopez LS, Hewitt BA, Sharma S. Reaching a breaking point: How is climate change influencing the timing of ice breakup in lakes across the Northern Hemisphere?. Limnology and Oceanography, 2019, 64(6): 2621-2631. DOI:10.1002/lno.11239
[26]
Clancy R, Bitz C, Blanchard-Wrigglesworth E. The influence of ENSO on Arctic Sea ice in large ensembles and observations. Journal of Climate, 2021, 1-50. DOI:10.1175/jcli-d-20-0958.1
[27]
Zhang XD, Sorteberg A, Zhang J et al. Recent radical shifts of atmospheric circulations and rapid changes in Arctic climate system. Geophysical Research Letters, 2008, 35(22): L22701. DOI:10.1029/2008GL035607
[28]
Overland JE, Wang M. The Arctic climate paradox: The recent decrease of the Arctic oscillation. Geophysical Research Letters, 2005, 32(6): L06701. DOI:10.1029/2004GL021752
[29]
Armitage TWK, Bacon S, Kwok R. Arctic Sea level and surface circulation response to the Arctic oscillation. Geophysical Research Letters, 2018, 45(13): 6576-6584. DOI:10.1029/2018GL078386
[30]
Gloersen P. Modulation of hemispheric sea-ice cover by ENSO events. Nature, 1995, 373(6514): 503-506. DOI:10.1038/373503a0
[31]
Wrzesiński D, Ptak M, Plewa K. Effect of the north Atlantic oscillation on water level fluctuations in lakes of northern Poland. Geographia Polonica, 2018, 91(2): 243-259. DOI:10.7163/gpol.0119
[32]
Gong DY, Guo D, Gao YQ et al. Boreal winter Arctic oscillation as an indicator of summer SST anomalies over the western tropical Indian Ocean. Climate Dynamics, 2017, 48(7): 2471-2488. DOI:10.1007/s00382-016-3216-2
[33]
Barzegar R, Aalami MT, Adamowski J. Coupling a hybrid CNN-LSTM deep learning model with a boundary corrected maximal overlap discrete wavelet transform for multiscale lake water level forecasting. Journal of Hydrology, 2021, 598: 126196. DOI:10.1016/j.jhydrol.2021.126196
[34]
德吉玉珍. 羌塘高原湖泊面积变化及其对气候变化的响应[学位论文]. 南京: 南京信息工程大学, 2022.
[35]
Xing LT, Huang LX, Chi GY et al. A dynamic study of a Karst spring based on wavelet analysis and the Mann-Kendall trend test. Water, 2018, 10(6): 698. DOI:10.3390/w10060698
[36]
Zhang AL, Gao RZ, Liu YX et al. Identification on hydrometeorology mutation characteristics and ecological evolution pattern of the plateau inland river basin—Taken Xilin River and Balager River of Inner Mongolia for instance. China Environmental Science, 2019, 39(12): 5254-5263. [张阿龙, 高瑞忠, 刘廷玺等. 高原内陆河流域气候水文突变与生态演变规律——以内蒙古锡林河和巴拉格尔河流域为例. 中国环境科学, 2019, 39(12): 5254-5263.]
[37]
Bliss A, Hock R, Radić V. Global response of glacier runoff to twenty-first century climate change. Journal of Geophysical Research: Earth Surface, 2014, 119(4): 717-730. DOI:10.1002/2013jf002931
[38]
Jansson P, Hock R, Schneider T. The concept of glacier storage: A review. Journal of Hydrology, 2003, 282(1/2/3/4): 116-129. DOI:10.1016/s0022-1694(03)00258-0
[39]
Colby GA, Ruuskanen MO, St Pierre KA et al. Warming climate is reducing the diversity of dominant microbes in the largest high Arctic lake. Frontiers in Microbiology, 2020, 11: 561194. DOI:10.3389/fmicb.2020.561194
[40]
Holland MM. The north Atlantic oscillation-Arctic oscillation in the CCSM2 and its influence on Arctic climate variability. Journal of Climate, 2003, 16(16): 2767-2781. DOI:10.1175/1520-0442(2003)0162767:tnaooi>2.0.co;2
[41]
Han ZX, Su T, Zhi R et al. Effects of atmospheric water budget changes under different Pacific interdecadal oscillations and ENSO phases on Pacific evaporation in winter in the Northern Hemisphere. Chinese Journal of Atmospheric Sciences, 2017, 41(6): 1316-1331. [韩子轩, 苏涛, 支蓉等. 不同太平洋年代际振荡和ENSO位相下大气水分收支变化对北半球冬季太平洋蒸发量的影响. 大气科学, 2017, 41(6): 1316-1331. DOI:10.3878/j.issn.1006-9895.1702.16257]
[42]
Wanner H, Brönnimann S, Casty C et al. North Atlantic oscillation—Concepts and studies. Surveys in Geophysics, 2001, 22(4): 321-381. DOI:10.1023/A:1014217317898
[43]
泮琬楠. 北半球中高纬前后冬气温转向与AO/NAO位相转变的联系[学位论文]. 南京: 南京信息工程大学, 2021.