冰川冻土, 2022, 44(1): 133-146 doi: 10.7522/j.issn.1000-0240.2022.0025

冰冻圈与全球变化

三江源地区1980—2019年积雪时空动态特征及其对气候变化的响应

陈龙飞,1,2, 张万昌,1, 高会然1,2

1.中国科学院 空天信息创新研究院 数字地球重点实验室,北京 100094

2.中国科学院大学,北京 100049

Spatiotemporal dynamic characteristics of snow cover from 1980 to 2019 in the Three-River-Source region and its response to climate change

CHEN Longfei,1,2, ZHANG Wanchang,1, GAO Huiran1,2

1.Key Laboratory of Digital Earth Science,Aerospace Information Research Institute,Chinese Academy of Sciences,Beijing 100094,China

2.University of Chinese Academy of Sciences,Beijing 100049,China

通讯作者: 张万昌,研究员,主要从事冰冻圈及水资源量的变化研究. E-mail: zhangwc@radi.ac.cn

收稿日期: 2021-08-21   修回日期: 2022-02-03  

基金资助: 青海省重点研发与转化计划项目.  2020-SF-C37

Received: 2021-08-21   Revised: 2022-02-03  

作者简介 About authors

陈龙飞,硕士研究生,主要从事积雪、冰川及水资源量的变化研究.E-mail:chenlongfei19@mails.ucas.ac.cn , E-mail:chenlongfei19@mails.ucas.ac.cn

摘要

三江源地区气象站点稀疏,依靠地面台站数据难以反映地面真实积雪情况。利用卫星遥感数据引入重心模型分析了三江源地区1980—2019年4个积雪参数(积雪日数、积雪深度、积雪初日和积雪终日)的时空动态特征,利用Mann-Kendall检验和Sen斜率估计分析了积雪和气候因子的变化趋势,并探究积雪对气候变化的响应。结果表明:1980—2019年三江源地区呈现积雪日数和积雪深度减少、积雪初日推迟、积雪终日提前的变化趋势,而该区域同期的气温和降水量则呈现上升趋势;4个积雪参数重心均呈现出东移趋势,而同期气温重心则呈现西移趋势,气温重心位置西移速率分别是积雪日数和积雪深度重心位置东移速率的6倍和2倍。这表明该区域4个积雪参数以及气候因子的变化趋势具有较强的空间异质性,西部气温升高速率大于东部,导致西部积雪日数和积雪深度减少速率同样大于东部,从而导致气温重心西移而积雪参数重心东移。澜沧江源区积雪日数减少、积雪深度减少、积雪初日推迟以及积雪终日提前的速率最大,其次是长江源区和黄河源区。进一步的相关性分析表明,三江源地区年平均气温的升高是导致积雪日数和积雪深度减少、积雪初日推迟、积雪终日提前的主要影响因子,积雪日数对气温升高响应最敏感,其次是积雪深度、初日和终日;而年降水量与4个积雪参数的相关性均不显著。研究可为三江源地区水资源和生态环境保护提供基础资料和理论依据。

关键词: 三江源地区 ; 积雪 ; 重心模型 ; 气候变化

Abstract

Meteorological stations in the Three-River-Source region is sparse, and only analyzing ground station data fails to reflect true snow cover on ground. In this study, based on remote sensing data, we applied gravity center model to analyze spatiotemporal dynamic characteristics of four snow cover parameters including snow cover days, snow depth, snow cover onset date and snow cover end date from 1980 to 2019. Mann-Kendall test and Sen’s slope estimation were used to analyze trends of the four snow cover parameters and climate factors. We finally explored its response to variations in annual mean air temperature and annual precipitation in the Three-River-Source region. The results indicated that snow cover days and snow depth showed downward trends, and snow cover onset date delayed and snow cover end date advanced from 1980 to 2019 in the Three-River-Source region, while air temperature and precipitation showed upward trends. Gravity center of the four snow cover parameters showed eastward trends, while gravity center of air temperature showed a westward trend. The westward speed of gravity center of air temperature was 6 times and 2 times the eastward speed of snow cover days and snow depth, respectively. These indicated that trends of the four snow cover parameters and climate factors in this region showed strong spatial heterogeneity. The increase rate of air temperature in the west was greater than that in the east, which leads to a greater decrease rate of snow cover days and snow depth in the west than in the east. It was responsible for the westward of gravity center of air temperature, and the eastward of gravity center of snow cover parameters. The Lancang River Source region had the highest rate of decrease in snow cover days and snow depth, delay of snow cover onset date and advance of snow cover end date, followed by the Yangtze River Source region and the Yellow River Source region. Further correlation analysis indicated that the increase in air temperature was primarily responsible for the decrease in snow cover days and snow depth, delay of snow cover onset date and advance of snow cover end date. Snow cover days was most sensitive to air temperature rise, followed by snow depth, snow cover onset date and end date. The correlation between annual precipitation and the four snow cover parameters was not significant. This research can provide basic data and theoretical basis for the protection of water resources and ecological environment in the Three-River-Source region.

Keywords: Three-River-Source region ; snow cover ; gravity center model ; climate change

PDF (10964KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

陈龙飞, 张万昌, 高会然. 三江源地区1980—2019年积雪时空动态特征及其对气候变化的响应[J]. 冰川冻土, 2022, 44(1): 133-146 doi:10.7522/j.issn.1000-0240.2022.0025

CHEN Longfei, ZHANG Wanchang, GAO Huiran. Spatiotemporal dynamic characteristics of snow cover from 1980 to 2019 in the Three-River-Source region and its response to climate change[J]. Journal of Glaciology and Geocryology, 2022, 44(1): 133-146 doi:10.7522/j.issn.1000-0240.2022.0025

0 引言

三江源地区是长江、黄河、澜沧江三条大河的发源地,水资源量丰富,是中国淡水资源的重要补给地。同时,也是高海拔地区中生物多样性最丰富的地区之一。2016年8月,习近平总书记在青海考察时指出,“中华水塔”是国家的生命之源,保护好三江源,对中华民族发展至关重要1。三江源地区的生态系统服务功能具有重要且不可替代的地位,对中国的生态状况及国民经济发展起着不可或缺的作用2-4

积雪是冰冻圈的重要组成部分,也是对气候变化反应最为敏感的环境要素之一,通常作为检验与监测全球变化的一个重要指标5-7。中国是中、低纬度地区冰川和积雪最为丰富的国家,冬季积雪大约相当于7.4×1010 m3水量,而青藏高原是我国三大积雪区之一,同时也是对全球气候变化反应最为敏感的区域之一8-9。研究表明,在全球变暖的背景下,青藏高原地区呈现积雪日数、雪水当量减少,积雪初日推迟而积雪终日提前的变化趋势10-12。目前绝大多数积雪遥感相关研究主要针对整个青藏高原地区,对于三江源地区积雪时空变化的研究较少,而利用遥感手段研究该区域积雪特征目前尚未有报道。开展三江源地区积雪遥感研究探究积雪时空特征,可为该区域水资源和生态环境保护提供科学资料以及理论依据。

获取可靠的积雪数据是研究积雪时空动态特征的基础,积雪数据的获取手段目前主要有地面台站观测以及卫星遥感反演。地面台站观测可以直接获取积雪参数,是一种非常有效的观测手段13。刘晓娇等14利用8个气象站点记录的逐日积雪深度数据研究了黄河源区1978—2016年积雪变化特征,结果表明黄河源区呈现出积雪初日推迟、终日提前和积雪日数减少的趋势;许显花等15利用2个气象站点资料分析了青海黄南南部区域1960—2015年的积雪变化,结果表明该区域积雪日数呈增加趋势,但趋势并不显著。这些结果虽然能够反映研究区积雪总体变化特征,但单点观测结果只能代表台站附近的情况,而对于广大未设置台站的区域,仍然无法获取数据。卫星遥感技术具有资料获取速度快、覆盖范围广、受地面条件限制少、信息量丰富等特点,被广泛应用于积雪变化研究。杨志刚等16利用MOD10A1积雪产品分析了2000—2014年青藏高原积雪面积和覆盖率的时空分布和变化特点;王晓茹等17利用逐日MODIS积雪产品分析了亚洲高山区2001—2016年融雪末期雪线高度的时空变化特征;Wang等6利用AMSR-E和MODIS积雪产品研究了青藏高原2003—2010年积雪日数和积雪面积的变化趋势;Dai等18利用AMSR-E和MODIS数据对青藏高原积雪深度进行了估算。这些研究表明,应用积雪遥感产品是目前进行大尺度研究优先选择的数据源,它能够全面地反映地面积雪状况,即便在站点稀疏的高海拔地区如青藏高原,同样能够快速获取可靠的积雪资料。三江源地区地处青藏高原东部,整体海拔较高,其西部区域气象站点稀疏,依靠地面台站数据几乎无法反映地面真实积雪情况。因此,本研究采用卫星遥感数据来研究三江源地区的积雪动态特征。

目前,大多数研究主要探究积雪特征变量的时间变化特征。许显花等15分析了青海黄南南部区域1960—2015年的积雪变化趋势,但未对积雪的空间变化进行探究;除多等19分析了青藏高原1981—2010年积雪日数的变化趋势,发现青藏高原地区年积雪日数出现显著减少趋势,冬季减幅最为明显,但对于该区域积雪的空间变化仍然只是定性分析。重心模型能够用于衡量某种属性在区域总体的分布状况,其变化趋势则可以揭示该属性整体的时空演变格局20。目前重心模型已被广泛应用于人口、经济、气象、环境、植被等要素的时空演变格局研究20-25。鉴于此,本研究引入重心模型对三江源地区积雪的空间变化情况进行定量分析。

三江源地区近几十年来积雪及气候因子具有怎样的时空动态特征,该区域积雪变化和气温之间存在怎样的关系是本研究重点解决的问题。针对这两个问题,利用逐日无云积雪面积和逐日积雪深度遥感产品,引入重心模型对三江源地区1980—2019年4个积雪参数(积雪日数、积雪深度、积雪初日和积雪终日)以及气温和降水量的空间演变格局进行研究,并利用Mann-Kendall(M-K)检验和Sen斜率估计分析积雪参数和气候因子的变化趋势,进而探究积雪对气候变化的响应,并为三江源地区水资源和生态环境保护提供依据。

1 数据与方法

1.1 研究区概况

三江源地区位于中国西部,是青藏高原的组成部分(图1)。其地理坐标为89°24′~102°27′ E,31°32′~37°08′ N,东西长1 156 km,南北宽592 km,总面积达36.9万km2。该区域整体海拔较高,且海拔呈现自西向东逐渐降低的特点。因此该区域气象站点主要分布在东部海拔较低的黄河源区,而西部区域气象站点较少。该区域雪山、冰川分布广阔,水资源丰富,是长江、黄河、澜沧江三条大河的发源地。本文根据集水区将研究区划分为长江源区、澜沧江源区、黄河源区。

图1

图1   研究区地形及气象站点分布

Fig. 1   Topography and distribution of meteorological stations over the study area


1.2 数据来源

中国逐日无云积雪面积数据集26-28来源于国家冰川冻土沙漠科学数据中心(http://www.ncdc.ac.cn),时间和空间分辨率分别为1 d和5 km,时间跨度为1980—2019年。该数据是由美国国家海洋与大气管理局(NOAA)发布的AVHRR-CDR反射率产品经过处理得到的。处理步骤为:(1)以Landsat TM数据为真值,利用多级决策树积雪判别算法来训练AVHRR-CDR表面反射率数据,从而获取初级产品;(2)利用隐马尔科夫算法、多源数据融合方法对该初级产品进行去云处理,最终获得逐日无云积雪面积产品。

中国逐日积雪深度长时间序列数据集29来源于国家青藏高原科学数据中心(http://data.tpdc.ac.cn/zh-hans/),时间和空间分辨率分别为1 d和25 km,时间跨度为1980—2019年。该数据集是由被动微波数据反演得到的,原始数据来自美国国家雪冰数据中心(NSIDC)发布的SMMR(1978—1987年),SSM/I(1987—2007年)和SSMI/S(2008—2019年)逐日被动微波亮温数据,经过4个步骤的处理最终获得逐日积雪深度数据30-34:(1)不同传感器的交叉定标;(2)被动微波遥感雪深反演算法系数的确定;(3)识别积雪的被动微波遥感分类树算法;(4)积雪反演。目前已有多项研究应用该数据集来研究积雪的变化特征,表明该数据集能够有效地反映积雪的实际情况35-37

三江源地区1980—2019年逐日、逐月、逐年气象站点数据来源于国家气象科学数据中心(http://data.cma.cn/)。中国1980—2015年年平均气温和年降水量插值数据集来源于资源环境科学与数据中心(https://www.resdc.cn/Default.aspx),原始数据的空间分辨率为1 km。将年平均气温和年降水量数据分别重采样为5 km和25 km的空间分辨率,用于后续相关系数的计算。

1.3 研究方法

本研究利用重心模型来探究4个积雪参数以及气候因子的时空变化特征。重心是物理学中的概念,指物体各部分所受重力产生合力的作用点38。重心的计算方法见式(1)~(2)。

X=i=1nXiVi/i=1nVi
Y=i=1nYiVi/i=1nVi

式中:XY为物理量Vxy方向的距离权重和总量之比,反映了V在空间上最大值的所在位置,即重心;n为研究区域所包含的栅格个数;XiYi 分别为第i个栅格中心的地理坐标;Vi 为该栅格上的物理量;(XY)为物理量重心的地理坐标。随着时间的推移,积雪参数和气候因子的空间分布会发生变化,从而导致其重心的位置也发生变化39。因此,重心的移动距离和方向能够反映地理要素的空间变化趋势和轨迹,从而揭示地理要素的时空演变格局。

本研究利用Mann-Kendall(M-K)检验法40-42来分析4个积雪参数以及气候因子变化趋势的显著性。M-K检验法是世界气象组织推荐的一种非参数统计检验方法。在双侧检验下,给定显著性水平α时,当时间序列的标准化方差Z>Z1-α/2时,时间序列存在显著趋势。Z>1.96以及Z>2.576分别表示在显著性水平α=0.05(即P<0.05)以及α=0.01(即P<0.01)下,时间序列存在显著变化趋势。

Sen斜率估计方法43-44用于计算4个积雪参数以及气候因子的变化趋势大小,其计算方法为

b=medianXj-Xij-i, j>i

式中:b为时间序列的趋势大小;ji为年份;XjXi 分别为第j和第i年的变量值;median( )为取中位数的函数。计算所得结果为正值(负值)时,表示该时间序列具有上升(下降)趋势,趋势大小为b,计算结果为0则表明该时间序列不存在趋势。

本研究利用积雪日数、积雪深度、积雪初日和积雪终日4个关键积雪参数来探究三江源地区积雪的时空动态特征。本文定义9月1日至次年8月31日为一个水文年5945,积雪日数定义为一个水文年内像元被积雪覆盖的日数之和。为了避免短时间降雪事件的影响,积雪初日定义为如果水文年内像元首次出现积雪且往后连续14天均被积雪覆盖,则该像元首次出现积雪的日期即为积雪初日。相反,积雪终日定义为如果水文年内像元最后一次出现积雪且往前连续14天均被积雪覆盖,则该像元最后一次出现积雪的日期即为积雪终日46-48

2 结果与讨论

2.1 积雪和气候因子的时空变化特征

2.1.1 时间变化

三江源地区1980—2019年积雪日数整体呈现减少趋势,倾向率为-2.18 d·a-1图2(a)],这表明该区域每经过10年,积雪日数将会减少约22 d。1980—1988年该区域积雪日数波动较大,1988—1990年积雪日数减少幅度较大,由139.72 d减少至38.31 d,而后1990—1999年则呈现缓慢上升趋势。1999—2001年又出现快速减少趋势,由73.34 d减少至16.68 d;而2001—2019年积雪日数均保持在较低水平,变化范围为7.66~41.34 d。三江源地区1980—2019年积雪日数具有明显的季节性特征。积雪主要分布在春季和冬季,积雪日数最大和最小的月份分别为1月和8月,1—8月积雪日数逐渐减小,而8—12月则逐渐增大[图2(c)~(d)]。分析积雪深度可知,三江源地区1980—2019年的平均积雪深度整体呈现减少趋势,倾向率为-0.22 mm·a-1,这表明每经过10年,该区域的平均积雪深度减少2.2 mm[图2(a)]。1985年该区域平均积雪深度最大,达62.30 mm,随后1986年则减少为12.19 mm。1986—1999年呈现缓慢上升趋势,1999—2001年则又出现快速减少趋势,由22.96 mm减少至7.32 mm;而2001—2019年积雪深度均保持在较低水平,变化范围为3.75~26.82 mm。积雪深度同样具有强烈的季节性特征。积雪深度较大的月份主要分布在10月至次年5月,10—12月积雪深度逐渐增大,而12月至次年5月积雪深度逐渐减小,6—9月则积雪深度几乎为0[图2(c)~(d)]。三江源地区积雪深度与积雪日数具有相似的变化特征,尤其是1999年之后,两者的变化趋势几乎一致,表明本文所采用的积雪覆盖和积雪深度数据集具有良好的一致性,能够有效地反映三江源地区真实的积雪动态情况。结合积雪日数及深度变化特征可知,该区域2001—2019年积雪日数和深度均明显小于1980—2000年的积雪日数和深度。

图2

图2   三江源地区1980—2019年积雪日数、积雪深度、积雪初日和积雪终日时间序列

Fig. 2   Time series of snow cover days, snow depth, snow cover onset date and end date in the Three-River-Source region from 1980 to 2019: yearly variation of snow cover days and snow depth (a), yearly variation of snow cover onset date and end date (b), seasonal variation of snow cover days and snow depth (c) and monthly variation of snow cover days and snow depth (d)


三江源地区积雪初日呈现缓慢的上升趋势,倾向率为0.35 d·a-1图2(b)],这表明该区域每经过10年,积雪初日推迟约3.5 d。积雪初日最早的年份为2009年,10月21日开始出现连续积雪覆盖;而积雪初日最晚的为2018年,直到次年1月27日才出现连续积雪覆盖。而积雪终日则呈现下降趋势,倾向率为-1.22 d·a-1图2(b)],这表明该区域每经过10年,积雪终日提前约12 d。2009年积雪终日时间最早,为11月16日;而2018年最晚,为次年3月28日。由此可见,三江源地区积雪初日逐渐推迟,而积雪终日逐渐提前,且积雪终日的提前速率是积雪初日推迟速率的3.5倍,这也是造成该区域积雪日数快速减少的原因。

通过分析三江源地区同期气温和降水量可知,无论月、季还是年尺度变化,三江源地区的气温和降水量均呈现出相似的变化趋势(图3)。三江源地区月平均气温在-15~15 ℃之间,每年7月的月平均气温最高,而1月最低;月降水量则为0~160 mm,最大和最小值同样分别出现在7月和1月。该区域年平均气温和年降水量均呈现上升的趋势,二者倾向率分别为0.05 ℃·a-1和1.87 mm·a-1图3(b)]。这表明每经过10年,该区域的年平均气温升高0.5 ℃,而年降水量则会增加18.7 mm。

图3

图3   三江源地区1980—2019年气温和降水量时间序列

Fig. 3   Time series of air temperature and precipitation in the Three-River-Source region from 1980 to 2019: monthly variation (a), yearly variation (b) and seasonal variation (c)


由此可见,三江源地区积雪、气温和降水量均具有强烈的季节性特征,1980—2019年三江源地区积雪日数和积雪深度呈现减少趋势,而同期气温和降水量则呈现上升趋势。气温升高会导致部分降雪转变为降雨,导致该区域降雪时间推迟。同时气温升高会加快积雪消融速度,从而导致积雪消融完毕时间逐渐提前。因此,随着年平均气温的升高,三江源地区积雪日数和积雪深度逐渐减少,积雪初日逐渐推迟,而积雪终日则逐渐提前。

2.1.2 空间变化

就空间变化特征而言,三江源地区积雪日数和积雪深度具有相似的空间分布特征(图4)。1980—1998年积雪日数和深度较大的区域分布较广阔,主要分布在该区域的西部和中南部地区,且积雪区域主要呈现面状分布;而2004年之后,积雪区域只零散的分布在三江源地区海拔较高的区域。随着时间的推移,积雪初日的图像颜色逐渐加深,而积雪终日则逐渐变浅(图5),表明三江源地区积雪初日逐渐推迟,积雪终日则逐渐提前。此外,积雪初日和终日的分布面积正在逐渐减小,1980—1998年主要呈现大面积面状分布;2004年之后则仅零散地分布在三江源地区海拔较高的区域,表明三江源地区能够保持连续14天积雪覆盖的区域正在快速减少。由此可见,三江源地区积雪区域面积正在逐渐缩小,由大面积的面状分布逐渐变成零散的点状分布,且积雪初日逐渐推迟,而积雪终日逐渐提前。

图4

图4   三江源地区1980—2019年部分年份积雪日数和积雪深度空间分布

Fig. 4   Spatial distributions of snow cover days (a) and snow depth (b) in some years from 1980 to 2019 in the Three-River-Source region


图5

图5   三江源地区1980—2019年部分年份积雪初日和积雪终日空间分布

Fig. 5   Spatial distributions of snow cover onset date (a) and end date (b) in some years from 1980 to 2019 in the Three-River-Source region


为了定量化衡量积雪的空间变化情况,本文分别计算了三江源地区积雪日数、深度、初日和终日1980—2019年逐年重心坐标并绘制其时间序列图[图6(a)~(d)]。结果表明,近40年来,三江源地区积雪日数重心在93.58°~96.54° E和33.87°~34.28° N之间变化,积雪深度重心在94.20°~97.62° E和33.29°~34.32° N之间变化。积雪日数、深度、初日和终日的重心位置在经度方向均呈现上升趋势,倾向率分别为0.0112° E·a-1、0.0249° E·a-1、0.0104° E·a-1和0.0123° E·a-1,而在纬度方向位置变化较小。这表明三江源地区的积雪日数、深度、初日和终日的分布区域逐渐向东移动,且每经过10年,其重心位置分别向东移动约0.112°(约10 km)、0.249°(约23 km)、0.104°(约10 km)和0.123°(约11 km)。

图6

图6   三江源地区4个积雪参数(1980—2019年)和年平均气温、年降水量(1980—2015年)的重心位置变化

Fig. 6   Changes in gravity center of the four snow cover parameters from 1980 to 2019, and annual mean air temperature and annual precipitation from 1980 to 2015 in the Three-River-Source region: snow cover days (a), snow depth (b), snow cover onset date (c), snow cover end date (d), annual mean air temperature (e) and annual precipitation (f)


采用相同的方法分别对研究区1980—2015年的年平均气温和年降水量数据进行处理[图6(e)~(f)],结果表明,1980—2015年三江源地区年平均气温的重心位置在经度方向呈现下降趋势,倾向率为-0.0651° E·a-1,在纬度方向变化趋势较小;而年降水量的重心位置在经度和纬度方向几乎无变化。由此可见,每经过10年,该区域年平均气温的重心位置向西移动约0.651°(约60 km),速率分别是积雪日数和积雪深度重心位置向东移动速率的6倍和2倍。

由此可见,三江源地区1980—2019年积雪日数、深度、初日和终日的重心位置呈现出东移的变化趋势,而同时期的年平均气温重心呈现出西移的变化趋势。重心位置的变化实际上是由地理要素在不同区域之间的变化存在差异而导致的。对于4个积雪参数而言,在三江源地区积雪逐渐减少的变化趋势下,该区域西部地区的积雪减少速率大于东部区域,从而导致4个积雪参数重心逐渐东移。对于年平均气温而言,在三江源地区气温逐渐升高的变化趋势下,该区域西部的气温升高速率大于东部区域,从而导致年平均气温重心逐渐西移。由此可见,三江源地区积雪参数和气候因子的变化趋势均呈现出较强的空间异质性。西部气温升高速率大于东部导致西部积雪日数和深度减少速率同样大于东部,从而导致气温重心西移而积雪参数重心东移。

2.2 积雪和气候因子的变化趋势

上文结果表明三江源地区积雪呈现减少趋势,而年平均气温和年降水量则呈现上升趋势,但其变化趋势均呈现较强的空间异质性。为了定量探究积雪日数、积雪深度、积雪初日、积雪终日、年平均气温和年降水量的变化趋势在不同区域的差异情况,本文联合M-K检验和Sen斜率估计方法从栅格层面分别检验了6个变量变化趋势的显著性以及趋势大小。

表1可知,澜沧江源区的积雪日数和积雪深度均呈现减少趋势,且减少趋势最大,分别为 -2.44 d·a-1和-0.32 mm·a-1;其次是长江源区,分别为-1.79 d·a-1和-0.17 mm·a-1;而对于黄河源区,其积雪日数呈现减少趋势(-1.36 d·a-1),而积雪深度却呈现微弱的上升趋势(0.03 mm·a-1),但绝大多数区域并未通过显著性检验[图7(b)]。对于积雪日数变化趋势的空间分布而言[图7(a)],三江源地区绝大多数区域的积雪日数都呈现出显著减少趋势,其中减少速率较大的区域主要分布在长江源区东部和南部、澜沧江源区以及黄河源区的中部和南部,倾向率达-5 d·a-1;而减少速率较小的区域则主要分布在黄河源区的东北部,倾向率约为-1 d·a-1。积雪深度的变化趋势与积雪日数呈现出相似的空间分布,大部分区域同样呈现出显著减少趋势,减少速率较大的区域主要分布在澜沧江源区北部和长江源区西南部,倾向率达-0.8 mm·a-1;而在澜沧江源区南端以及黄河源区的东南部部分区域,积雪深度呈现出轻微上升趋势,但该趋势并不显著。

表1   三江源地区4个积雪参数(1980—2019年)和年平均气温、年降水量(1980—2015年)的倾向率

Table 1  Tendency of the four snow cover parameters from 1980 to 2019, and annual mean air temperature and annual precipitation from 1980 to 2015 in the Three-River-Source region

地区积雪日数/(d·a-1积雪深度/(mm·a-1积雪初日/(d·a-1积雪终日/(d·a-1年平均气温/(℃·a-1年降水量/(mm·a-1
长江源区-1.79-0.170.25-0.430.061.17
黄河源区-1.360.030.24-0.310.05-0.06
澜沧江源区-2.44-0.321.18-1.080.051.00

新窗口打开| 下载CSV


图7

图7   三江源地区4个积雪参数(1980—2019年)和年平均气温、年降水量(1980—2015年)趋势显著性及倾向率的空间分布

Fig. 7   Spatial distributions of trend significance and tendency of the four snow cover parameters from 1980 to 2019, and annual mean air temperature and annual precipitation from 1980 to 2015 in the Three-River-Source region: snow cover days (a), snow depth (b), snow cover onset date (c), snow cover end date (d), annual mean air temperature (e) and annual precipitation (f)


三个源区均呈现出积雪初日推迟、积雪终日提前的趋势(表1)。澜沧江源区的变化趋势最大,积雪初日推迟、积雪终日提前速率分别为1.18 d·a-1和-1.08 d·a-1;其次是长江源区,分别为0.25 d·a-1和-0.43 d·a-1;最后是黄河源区,分别为0.24 d·a-1和-0.31 d·a-1。积雪初日和终日的变化趋势具有相似的空间分布,倾向率较大的区域主要分布在长江源区东部、澜沧江源区北部以及黄河源区西部,积雪初日和终日的倾向率分别达到2 d·a-1和-3 d·a-1。这表明每经过10年,该区域积雪初日推迟20 d,同时积雪终日提前30 d。

对于气候要素而言,长江源区、澜沧江源区和黄河源区年均气温倾向率分别为0.06 ℃·a-1、0.05 ℃·a-1和0.05 ℃·a-1。除了黄河源区南端小部分区域外,整个三江源地区的年均气温均呈现出显著上升趋势,且其倾向率自东南向西北逐渐增大。黄河源区东南部区域的倾向率为0~0.02 ℃·a-1,而长江源区南部和西北部的倾向率则超过0.06 ℃·a-1。长江源区和澜沧江源区同期的年降水量呈现上升趋势,倾向率分别为1.17 mm·a-1和1.00 mm·a-1;而黄河源区则呈现下降趋势,倾向率为-0.06 mm·a-1,但几乎整个三江源地区年降水量的变化趋势均未通过显著性检验。

由此可见,三江源地区绝大部分区域1980—2019年积雪日数和积雪深度呈现显著减少的变化趋势,而同期该区域气候则朝暖湿方向变化。胡豪然等49利用38个站点资料研究了青藏高原东部地区1961—2010年冬季积雪的时空变化特征;姜琪等50利用110个气象站数据研究了青藏高原地区1961—2014年积雪时空特征;车涛等51利用AVHRR以及被动微波数据研究了青藏高原地区1980—2016年的积雪变化特征,均发现积雪日数或积雪深度呈现减少的变化趋势。本研究还发现黄河源区积雪深度呈现微弱的上升趋势(0.03 mm·a-1),这与刘晓娇等14利用8个气象站点数据研究黄河源区1978—2016年积雪深度时所得结果一致,其上升趋势同样为0.03 mm·a-1。由此表明,本研究所采用的积雪面积数据集以及积雪深度数据集能够很好地反映三江源地区积雪日数和积雪深度的真实情况。此外,三江源地区积雪参数的变化与该区域气候要素的变化具有一定的相关关系。

2.3 积雪对气候变化的响应

上文研究结果表明三江源地区4个积雪参数的变化与该区域气候要素的变化具有一定的相关关系。为了定量探究4个积雪参数对气候变化的响应,本文分别计算了三江源地区1980—2015年积雪日数、积雪深度、积雪初日、积雪终日与年平均气温和年降水量的相关系数。

表2可知,整个三江源地区积雪日数与气温呈现负相关关系,相关系数达-0.52,通过了P<0.001的显著性检验。其中长江源区、黄河源区、澜沧江源区的负相关系数依次为-0.47、-0.57、 -0.65,均通过了P<0.01的显著性检验。就相关系数空间分布而言,除极少数栅格外,几乎整个三江源地区积雪日数与气温都呈现正相关关系[图8(a)],其中相关系数较大的区域主要分布在长江源区西部和东部、澜沧江源区北部以及黄河源区的中部和西部,相关系数超过-0.6,部分地区甚至超过-0.8。整个三江源地区积雪深度与气温同样呈现负相关关系,相关系数为-0.40,长江源区、黄河源区、澜沧江源区的相关系数分别为-0.46、-0.34、-0.34,均通过P<0.05的显著性检验。但相比于积雪日数与气温的相关系数,积雪深度的相关系数相对较小。积雪深度与气温的相关系数较大的区域主要分布在长江源区、澜沧江源区北部以及黄河源区西部,系数值超过-0.4[图8(b)]。澜沧江源区南部以及黄河源区东部区域虽出现少量正相关系数,但均未通过显著性检验。由此可见,三江源地区积雪日数和深度随着该区域年平均气温的升高在逐渐减少,且积雪日数对于年平均气温升高的响应更为敏感。

表2   三江源地区1980—2015年4个积雪参数与年平均气温、年降水量的相关系数均值

Table 2  Mean correlation coefficients between the four snow cover parameters and annual mean air temperature, annual precipitation from 1980 to 2015 in the Three-River-Source region

相关参数三江源地区长江源区黄河源区澜沧江源区
积雪日数-气温-0.52***-0.47**-0.57***-0.65***
积雪深度-气温-0.40*-0.46**-0.34*-0.34*
积雪初日-气温0.34*0.39*0.260.43**
积雪终日-气温-0.33*-0.38*-0.25-0.41**
积雪日数-降水量0.080.080.100.01
积雪深度-降水量0.190.220.150.19
积雪初日-降水量-0.06-0.07-0.090.04
积雪终日-降水量0.050.050.07-0.04

注:*、**、***分别表示通过了P<0.05、P<0.01、P<0.001的显著性检验。在本研究中,样本数量为36时,相关系数绝对值分别大于0.32、0.41、0.51。

新窗口打开| 下载CSV


图8

图8   三江源地区1980—2015年4个积雪参数与年平均气温的相关系数空间分布

Fig. 8   Spatial distributions of correlation coefficients between the four snow cover parameters and annual mean air temperature from 1980 to 2015 in the Three-River-Source region: snow cover days and air temperature (a), snow depth and air temperature (b), snow cover onset date and air temperature (c) and snow cover end date and air temperature (d)


整个三江源地区积雪初日与气温呈现显著的正相关关系(P<0.05),其中相关性最强的为澜沧江源区(0.43),通过P<0.01的显著性检验;其次为长江源区(0.39),通过P<0.05的显著性检验;而黄河源区相关系数最小(0.26),且相关性不显著。积雪初日与气温的正相关系数几乎遍布整个研究区域[图8(c)],其中相关性较强的区域主要分布在长江源区东部、澜沧江源区北部和黄河源区西部,相关系数超过0.4;而黄河源区东部则相关系数较小,甚至出现负相关系数,但相关性并不显著。而三江源地区积雪终日与气温的相关系数分布情况与积雪初日相似,但积雪终日与气温呈现负相关关系。由此可见,除黄河源区东部区域外,三江源地区积雪初日和终日与气温分别呈现显著正相关关系和显著负相关关系。即随着年平均气温的升高,该区域积雪初日逐渐推迟,而积雪终日则逐渐提前。

4个积雪参数与年降水量的相关性均不显著,因此本研究未进一步从栅格层面分析4个积雪参数与年降水量的相关系数空间分布情况。由于年降水量包含降雪量,理论上,当年平均气温恒定的情况下,年降水量的增加会导致降雪量也有所增加,从而使积雪日数和积雪深度增大。然而,由于三江源地区的年平均气温逐渐升高,导致降雪的温度条件更加难以满足,部分降雪转变为降雨,从而导致该区域降雪时间推迟52。此外气温升高会加快积雪消融速度,从而导致积雪消融完毕的时间逐渐提前;同时积雪消融期由于温度不足达到降雪条件,此时降雨量的增加也会加速积雪的消融9。因此,年平均气温的升高减弱了4个积雪参数与年降水量的相关关系,从而导致降水对积雪参数虽产生一定的影响但在数理统计上并不显著。谭秋阳等53研究了1979—2017年雅鲁藏布江流域积雪深度与气候因子的关系;白淑英等54研究了近30年来西藏地区积雪深度与气候因子的关系;刘金平等9研究了2000—2014年雅鲁藏布江流域积雪对气候变化的响应,研究结果均表明降水量与积雪参数的相关性不显著。因此,三江源地区积雪日数、深度、初日和终日的变化主要受年平均气温的影响。积雪日数对气温升高响应最敏感,其次是积雪深度、积雪初日和积雪终日;而年降水量与4个积雪参数的相关性均不显著。

3 结论

本文基于AVHRR-CDR逐日无云积雪面积数据集和SMMR、SSM/I、SSMI/S逐日被动微波亮温数据集引入重心模型分析了三江源地区1980—2019年积雪日数、积雪深度、积雪初日和积雪终日的时空动态特征,并探究了4个积雪参数对研究区年平均气温和年降水量变化的响应。研究结论如下:

(1)在时间变化特征上,1980—2019年三江源地区积雪日数和积雪深度均呈现减少的变化趋势,减少速率分别为-2.18 d·a-1和-0.22 mm·a-1;积雪初日逐渐推迟而积雪终日逐渐提前,变化速率分别为0.35 d·a-1和-1.22 d·a-1,积雪终日提前速率是积雪初日推迟速率的3.5倍。该区域同期的气温和降水量则呈现上升趋势,上升速率分别为0.05 ℃·a-1和1.87 mm·a-1。在空间变化特征上,4个积雪参数的重心呈现出东移的变化趋势,东移速率分别为 1.0 km·a-1、2.3 km·a-1、1.0 km·a-1、1.1 km·a-1。同期的年降水量重心几乎无变化,而年平均气温的重心则呈现出西移的变化趋势,西移速率为6.0 km·a-1,分别是积雪日数和积雪深度重心位置东移速率的6倍和2倍,表明三江源地区4个积雪参数和气候因子的变化趋势均呈现出较强的空间异质性。西部气温升高速率大于东部导致西部积雪日数和深度减少速率同样大于东部,从而导致气温重心西移而积雪参数重心东移。

(2)澜沧江源区的积雪日数减少、积雪深度减少、积雪初日推迟以及积雪终日提前速率最大,分别为-2.44 d·a-1、-0.32 mm·a-1、1.18 d·a-1和-1.08 d·a-1,其次是长江源区和黄河源区。除黄河源区部分区域的积雪深度呈现微弱的不显著上升趋势外,三江源地区积雪日数、积雪深度逐渐减少,积雪初日逐渐推迟而积雪终日逐渐提前。对于同期年平均气温和年降水量而言,除黄河源区的年降水量呈现微弱的下降趋势外,其余均呈现上升趋势。因此,三江源地区积雪参数的变化与该区域气候要素的变化具有一定的相关关系。

(3)进一步的相关性分析表明,整个三江源地区积雪日数、积雪深度、积雪终日与年平均气温均呈现显著负向相关关系(P<0.05),而积雪初日与气温则呈现显著正向相关关系(P<0.05)。其中相关系性较强的区域主要分布在长江源区东部,澜沧江源区北部以及黄河源区西部。三江源地区4个积雪参数与年降水量的相关性均未通过显著性检验。由此可见,三江源地区气温的升高是导致该区域积雪日数和积雪深度减少、积雪初日推迟而积雪终日提前的主要影响因子。积雪日数对气温升高响应最敏感,其次是积雪深度、积雪初日和积雪终日;而降水量的变化对于该区域积雪参数的变化影响较小。

因此,随着全球气候变暖的持续,三江源地区的积雪日数和积雪深度持续减少,同时积雪初日持续推迟而积雪终日持续提前,这可能导致积雪在春季来临之前完全融化。冬季积雪融水是植物生长季开始时土壤水分的主要来源,同时也是初春时期河流径流的主要来源。积雪的提前融化会影响该区域植物的生长,这对于生态环境原本就脆弱的三江源地区而言是不利的。春季河流径流量的减少也会影响下游地区农作物的生长,对区域粮食安全造成威胁。

参考文献

Xuan Wei.

Protection of the Three-River-Source region, how to act in culture

[J]. China Water Resources, 201617): 1-6.

[本文引用: 1]

轩玮.

三江源保护, 文化如何行动

[J]. 中国水利, 201617): 1-6.

[本文引用: 1]

Sun QinglingLi BaolinLi Feiet al.

Review on the estimation of net primary productivity of vegetation in the Three-River Headwater region, China

[J]. Acta Geographica Sinica, 2016719): 1596-1612.

[本文引用: 1]

孙庆龄李宝林李飞.

三江源植被净初级生产力估算研究进展

[J]. 地理学报, 2016719): 1596-1612.

[本文引用: 1]

Li SisiYao ZhijunWang Ruiet al.

Dryness/wetness pattern over the Three-River Headwater region: variation characteristic, causes, and drought risks

[J]. International Journal of Climatology, 2020407): 3550-3566.

Wang SuiziFan JiangwenLi Yuzheet al.

Dynamic response of water retention to grazing activity on grassland over the Three River Headwaters region

[J]. Agriculture, Ecosystems & Environment, 2019286106662.

[本文引用: 1]

Wang Hai’eLi ShengchenZhang Qingmeiet al.

Changes of the days with snow cover on the Qinghai Plateau during 1961—2013

[J]. Journal of Glaciology and Geocryology, 2016385): 1219-1226.

[本文引用: 2]

王海娥李生辰张青梅.

青海高原1961—2013年积雪日数变化特征分析

[J]. 冰川冻土, 2016385): 1219-1226.

[本文引用: 2]

Wang WeiHuang XiaodongDeng Jieet al.

Spatio-temporal change of snow cover and its response to climate over the Tibetan Plateau based on an improved daily cloud-free snow cover product

[J]. Remote Sensing, 201571): 169-194.

[本文引用: 1]

Wang HuiWang ShengliYu Xingjieet al.

Spatial-temporal variation of snow cover in Xinjiang based on surface observation from 1961 to 2017

[J]. Journal of Glaciology and Geocryology, 2020421): 72-80.

[本文引用: 1]

王慧王胜利余行杰.

1961—2017年基于地面观测的新疆积雪时空变化研究

[J]. 冰川冻土, 2020421): 72-80.

[本文引用: 1]

Qin DaheXiao CundeDing Yongjianet al.

Progress on cryospheric studies by international and Chinese communities and perspectives

[J]. Journal of Applied Meteorological Science, 2006176): 649-656.

[本文引用: 1]

秦大河效存德丁永建.

国际冰冻圈研究动态和我国冰冻圈研究的现状与展望

[J]. 应用气象学报, 2006176): 649-656.

[本文引用: 1]

Liu JinpingZhang WanchangDeng Caiet al.

Spatiotemporal variations of snow cover over Yarlung Zangbo River basin in Tibet from 2000 to 2014 and its response to key climate factors

[J]. Journal of Glaciology and Geocryology, 2018404): 643-654.

[本文引用: 4]

刘金平张万昌邓财.

2000—2014年西藏雅鲁藏布江流域积雪时空变化分析及对气候的响应研究

[J]. 冰川冻土, 2018404): 643-654.

[本文引用: 4]

Xu WenfangMa LijuanMa Minnaet al.

Spatial-temporal variability of snow cover and depth in the Qinghai-Tibetan Plateau

[J]. Journal of Climate, 2017304): 1521-1533.

[本文引用: 1]

Sun YanhuaHuang XiaodongWang Weiet al.

Spatio-temporal changes of snow cover and snow water equivalent in the Tibetan Plateau during 2003—2010

[J]. Journal of Glaciology and Geocryology, 2014366): 1337-1344.

孙燕华黄晓东王玮.

2003—2010年青藏高原积雪及雪水当量的时空变化

[J]. 冰川冻土, 2014366): 1337-1344.

Shen S S PYao RuzhenNgo Jet al.

Characteristics of the Tibetan Plateau snow cover variations based on daily data during 1997—2011

[J]. Theoretical and Applied Climatology, 20151203): 445-453.

[本文引用: 1]

Ma LijuanQin Dahe.

Spatial-temporal characteristics of observed key parameters for snow cover in China during

19572009[J]. Journal of Glaciology and Geocryology, 2012341): 1-11.

[本文引用: 1]

马丽娟秦大河.

1957—2009年中国台站观测的关键积雪参数时空变化特征

[J]. 冰川冻土, 2012341): 1-11.

[本文引用: 1]

Liu XiaojiaoChen RenshengLiu Junfenget al.

Variation of snow cover and its influence on spring runoff in the source region of Yellow River

[J]. Plateau Meteorology, 2020392): 226-233.

[本文引用: 2]

刘晓娇陈仁升刘俊峰.

黄河源区积雪变化特征及其对春季径流的影响

[J]. 高原气象, 2020392): 226-233.

[本文引用: 2]

Xu XianhuaLi YanlinLiu Yihuaet al.

Analysis on the variation of snow cover in the south of Huangnan in the last 56 years

[J]. Plateau and Mountain Meteorology Research, 2016364): 65-70.

[本文引用: 2]

许显花李延林刘义花.

黄南南部近56年积雪变化分析研究

[J]. 高原山地气象研究, 2016364): 65-70.

[本文引用: 2]

Yang ZhigangDawaChu Duo.

Spatiotemporal variations of snow cover on the Tibetan Plateau over the last 15 years

[J]. Remote Sensing Technology and Application, 2017321): 27-36.

[本文引用: 1]

杨志刚达娃除多.

近15 a青藏高原积雪覆盖时空变化分析

[J]. 遥感技术与应用, 2017321): 27-36.

[本文引用: 1]

Wang XiaoruTang ZhiguangWang Jianet al.

Monitoring of snowline altitude at the end of melting season in High Mountain Asia based on MODIS snow cover products

[J]. Acta Geographica Sinica, 2020753): 470-484.

[本文引用: 1]

王晓茹唐志光王建.

基于MODIS积雪产品的高亚洲融雪末期雪线高度遥感监测

[J]. 地理学报, 2020753): 470-484.

[本文引用: 1]

Dai LiyunChe TaoXie Hongjieet al.

Estimation of snow depth over the Qinghai-Tibetan Plateau based on AMSR-E and MODIS data

[J]. Remote Sensing, 20181012): 1989.

[本文引用: 1]

Chu DuoYang YongJiancan Luobuet al.

The variations of snow cover days over the Tibetan Plateau during 1981—2010

[J]. Journal of Glaciology and Geocryology, 2015376): 1461-1472.

[本文引用: 1]

除多杨勇罗布坚参.

1981—2010年青藏高原积雪日数时空变化特征分析

[J]. 冰川冻土, 2015376): 1461-1472.

[本文引用: 1]

Zhang GuojunHuang WanlingZhou Chunshanet al.

Spatio-temporal characteristics of demographic distribution in China from the perspective of urban agglomeration

[J]. Acta Geographica Sinica, 2018738): 1513-1525.

[本文引用: 2]

张国俊黄婉玲周春山.

城市群视角下中国人口分布演变特征

[J]. 地理学报, 2018738): 1513-1525.

[本文引用: 2]

Yan ZihongTang MenggeWei Lanet al.

Population distribution in northeastern Sichuan: spatio-temporal evolution characteristics based on GIS

[J]. Journal of Agriculture, 2020106): 94-100.

晏自红唐梦鸽魏兰.

基于GIS的川东北地区人口分布时空演变特征

[J]. 农学学报, 2020106): 94-100.

Chen SujingLi LijuanLi Jiuyiet al.

Analysis of the temporal and spatial variation characteristics of precipitation in the Lancang River basin over the past 55 years

[J]. Journal of Geo-Information Science, 2017193): 365-373.

陈素景李丽娟李九一.

近55年来澜沧江流域降水时空变化特征分析

[J]. 地球信息科学学报, 2017193): 365-373.

Wang HaijunZhang BinLiu Yaolinet al.

Urban expansion patterns and their driving forces based on the center of gravity-GTWR model: a case study of the Beijing-Tianjin-Hebei urban agglomeration

[J]. Journal of Geographical Sciences, 2020302): 297-318.

Wang ZhenChen ShutingCui Canet al.

Industry relocation or emission relocation? Visualizing and decomposing the dislocation between China’s economy and carbon emissions

[J]. Journal of Cleaner Production, 20192081109-1119.

Wang JianguoZhang Fei.

Spatial-temporal pattern and gravity center change of fractional vegetation cover in Xinjiang, China from 2000 to 2019

[J]. Transactions of the Chinese Society of Agricultural Engineering, 20203620): 188-194.

[本文引用: 1]

王建国张飞.

2000—2019年新疆植被覆盖度时空格局及重心变化分析

[J]. 农业工程学报, 20203620): 188-194.

[本文引用: 1]

Hao XiaohuaHuang GuanghuiChe Taoet al.

The NIEER AVHRR snow cover extent product over China: a long-term daily snow record for regional climate research

[J]. Earth System Science Data, 20211310): 4711-4726.

[本文引用: 1]

Wang XuanHao XiaohuaWang Jianet al.

Accuracy evaluation of long time series AVHRR snow cover area products in China

[J]. Remote Sensing Technology and Application, 2018336): 994-1003.

王轩郝晓华王建.

中国地区AVHRR长时间序列积雪范围产品精度评估

[J]. 遥感技术与应用, 2018336): 994-1003.

Hao XiaohuaJi WenzhengSun Xinglianget al.

Daily 5-km gap-free AVHRR snow cover extent product over China

[DS/OL]. Lanzhou: National Cryosphere Desert Data Center, 20202022-01-25]. .

[本文引用: 1]

郝晓华纪文政孙兴亮.

中国AVHRR逐日无云5 km积雪面积产品数据集

[DS/OL]. 兰州: 国家冰川冻土沙漠科学数据中心, 20202022-01-25]. .

[本文引用: 1]

Che TaoDai Liyun.

Long-term series of daily snow depth dataset in China (1979—2020)

[DS/OL]. Beijing: National Tibetan Plateau Data Center, 20202022-01-25] [. [车涛, 戴礼云. 中国雪深长时间序列数据集(19792020)[DS/OL]. 北京: 国家青藏高原科学数据中心, 2020 [2022-01-25]. ]

[本文引用: 1]

Che TaoLi XinJin Ruiet al.

Snow depth derived from passive microwave remote-sensing data in China

[J]. Annals of Glaciology, 2008491): 145-154.

[本文引用: 1]

Dai LiyunChe Tao.

Cross-platform calibration of SMMR, SSM/I and AMSR-E passive microwave brightness temperature

[C]// Sixth International Symposium on Digital Earth: data processing and applications. BeijingThe International Society for Digital Earth20107841103.

Dai LiyunChe TaoWang Jianet al.

Snow depth and snow water equivalent estimation from AMSR-E data based on a priori snow characteristics in Xinjiang, China

[J]. Remote Sensing of Environment, 20121271): 14-29.

Dai LiyunChe TaoDing Yongjian.

Inter-calibrating SMMR, SSM/I and SSMI/S data to improve the consistency of snow-depth products in China

[J]. Remote Sensing, 201576): 7212-7230.

Dai LiyunChe TaoDing Yongjianet al.

Evaluation of snow cover and snow depth on the Qinghai-Tibetan Plateau derived from passive microwave remote sensing

[J]. The Cryosphere, 2017114): 1933-1948.

[本文引用: 1]

Wang XuejiaoWang XiaomingLiu Shiweiet al.

Snow cover loss compounding the future economic vulnerability of western China

[J]. Science of the Total Environment, 2021755143025.

[本文引用: 1]

Luo JiangxinShihua Wang Tinget al.

Variation characteristics of snow cover and its influence on soil water and heat transfer on the Qinghai-Xizang Plateau

[J]. Plateau Meteorology, 2020396): 1144-1154.

罗江鑫吕世华王婷.

青藏高原积雪变化特征及其对土壤水热传输的影响

[J]. 高原气象, 2020396): 1144-1154.

Bi ZheruiSachulaLiu Guixiang.

Remote sensing analysis of temporal and spatial changes of snow depth in Inner Mongolia from 1982 to 2015

[J]. Grassland and Prataculture, 2019314): 12-18.

[本文引用: 1]

毕哲睿萨楚拉刘桂香.

1982—2015年内蒙古雪深时空变化遥感分析

[J]. 草原与草业, 2019314): 12-18.

[本文引用: 1]

Chen ShutingGuo BingYang Feiet al.

Spatial and temporal patterns of NPP and its response to climate change in the Qinghai-Tibet Plateau from 2000 to 2015

[J]. Journal of Natural Resources, 20203510): 2511-2527.

[本文引用: 1]

陈舒婷郭兵杨飞.

2000—2015年青藏高原植被NPP时空变化格局及其对气候变化的响应

[J]. 自然资源学报, 20203510): 2511-2527.

[本文引用: 1]

Cheng ChaoTong ShaoyuPeng Haiyinget al.

Ecological carrying capacity of water resources in the central Yunnan urban agglomeration area

[J]. Resources Science, 2016388): 1561-1571.

[本文引用: 1]

程超童绍玉彭海英.

滇中城市群水资源生态承载力的平衡性研究

[J]. 资源科学, 2016388): 1561-1571.

[本文引用: 1]

Mann H B.

Nonparametric tests against trend

[J]. Econometrica: Journal of the Econometric Society, 1945133): 245-259.

[本文引用: 1]

Kendall M G. Rank correlation methods[M]. LondonCharles Griffin & Co. Ltd.1975.

Gilbert R O. Statistical methods for environmental pollution monitoring[M]. Hoboken, New Jersey, USAJohn Wiley & Sons1987.

[本文引用: 1]

Theil H.

A rank-invariant method of linear and polynomial regression analysis

[J]. Proceedings of the Koninklijke Nederlandse Akademie van Wetenschappen, 1992122): 345-381.

[本文引用: 1]

Sen P K.

Estimates of the regression coefficient based on Kendall’s tau

[J]. Journal of the American Statistical Association, 196863324): 1379-1389.

[本文引用: 1]

Dietz A JKuenzer CConrad C.

Snow-cover variability in Central Asia between 2000 and 2011 derived from improved MODIS daily snow-cover products

[J]. International Journal of Remote Sensing, 20133411): 3879-3902.

[本文引用: 1]

Li ChenhaoSachulaLiu Guixianget al.

Spatiotemporal changes of snow cover and its response to climate changes in the Mongolian Plateau from 2000 to 2017

[J]. Chinese Journal of Grassland, 2020422): 95-104.

[本文引用: 1]

李晨昊萨楚拉刘桂香.

2000~2017年蒙古高原积雪时空变化及其对气候响应研究

[J]. 中国草地学报, 2020422): 95-104.

[本文引用: 1]

Wang ZengyanChe Tao.

Spatiotemporal distribution of snow cover in arid regions in China

[J]. Arid Zone Research, 2012293): 464-471.

王增艳车涛.

2002—2009年中国干旱区积雪时空分布特征

[J]. 干旱区研究, 2012293): 464-471.

Gao YangXie HongjieYao Tandong.

Developing snow cover parameters maps from MODIS, AMSR-E, and blended snow products

[J]. Photogrammetric Engineering & Remote Sensing, 201177351-361.

[本文引用: 1]

Hu HaoranLiang Ling.

Spatial and temporal variations of winter snow over east of Qinghai-Tibet Plateau in the last 50 years

[J]. Acta Geographica Sinica, 20136811): 1493-1503.

[本文引用: 1]

胡豪然梁玲.

近50年青藏高原东部冬季积雪的时空变化特征

[J]. 地理学报, 20136811): 1493-1503.

[本文引用: 1]

Jiang QiLuo SiqiongWen Xiaohanget al.

Spatial-temporal characteristics of snow and influence factors in the Qinghai-Tibetan Plateau from 1961 to 2014

[J]. Plateau Meteorology, 2020391): 24-36.

[本文引用: 1]

姜琪罗斯琼文小航.

1961—2014年青藏高原积雪时空特征及其影响因子

[J]. 高原气象, 2020391): 24-36.

[本文引用: 1]

Che TaoHao XiaohuaDai Liyunet al.

Snow cover variation and its impacts over the Qinghai-Tibet Plateau

[J]. Bulletin of Chinese Academy of Sciences, 20193411): 1247-1253.

[本文引用: 1]

车涛郝晓华戴礼云.

青藏高原积雪变化及其影响

[J]. 中国科学院院刊, 20193411): 1247-1253.

[本文引用: 1]

Yao YananXu JunrongGuo Lingpenget al.

Relationships between temperature and precipitation forms based on the data from meteorological stations: a case study in the Kaidu River basin

[J]. Arid Zone Research, 2015321): 94-101.

[本文引用: 1]

姚亚楠徐俊荣郭玲鹏.

气温与降水形态关系研究: 以开都河流域为例

[J]. 干旱区研究, 2015321): 94-101.

[本文引用: 1]

Tan QiuyangCheng LeiXu Zongxueet al.

Spatiotemporal distribution of snow cover depth and its driving factors in the Yarlung Zangbo River basin, 1979—2017

[J]. Journal of Glaciology and Geocryology, 2021434): 1049-1059.

[本文引用: 1]

谭秋阳程磊徐宗学.

1979—2017年雅鲁藏布江流域雪深时空分布特征及其影响因素分析

[J]. 冰川冻土, 2021434): 1049-1059.

[本文引用: 1]

Bai ShuyingShi JianqiaoShen Weishouet al.

Spatial-temporal variation of snow depth in Tibet and its response to climatic change in the past 30 years

[J]. Remote Sensing for Land & Resources, 2014261): 144-151.

[本文引用: 1]

白淑英史建桥沈渭寿.

近30年西藏雪深时空变化及其对气候变化的响应

[J]. 国土资源遥感, 2014261): 144-151.

[本文引用: 1]

/