基于GRACE数据的中国水储量变化特征分析
陈坤1,2,3, 蒋卫国2,*, 何福红3,*, 陈征2
1. 中国环境科学研究院,北京 100012
2. 北京师范大学地理科学学部环境演变与自然灾害教育部重点实验室, 北京 100875
3. 鲁东大学资源与环境工程学院,山东 烟台 264025
*通信作者简介:蒋卫国(1976- ),男,湖南衡阳人,教授,研究方向为洪水灾害及湿地生态遥感监测与评价研究。E-mail:jiangweiguo@bnu.edu.cn;何福红(1978—2017),男,湖南衡阳人,副教授,研究方向为土壤侵蚀与3S技术应用。

第一作者简介:陈坤(1989- ),男,山东威海人,硕士研究生,主要从事水资源遥感方面的研究。E-mail:chenkun0697@126.com

摘要

论文以中国2003—2015年GRACE水储量数据、GLDAS-Noah地表蒸散发数据、地面实测降水数据为数据源,基于Theil-Sen Median 趋势分析和Mann-Kendall检验法,分析了中国水储量、降水量以及蒸散发量的变化趋势及相关关系,并利用经验正交函数(EOF)分析了中国水储量变化时空分布特征,结果表明:1) 2006年以后中国整体水储量年际变化剧烈,且降水和蒸散发可以识别水储量异常时段内的特征信息;2) 2003—2015年全国各地水储量、降雨量和蒸散发量变化空间分布差异明显;3) 全国共有18.6%的地区气候因素对水储量变化的方差解释量超过50%,主要位于中国东北、东南沿海、四川盆地、青海高原以及新疆西北部地区;4) EOF分解的GRACE水储量各主成分信息可以表征研究时段内中国地区整体和局部的水储量时空变化特征。

关键词: GRACE; 水储量; Theil-Sen Median 趋势分析; Mann-Kendall检验法; EOF
中图分类号:TV213 文献标志码:A 文章编号:1000-3037(2018)02-0275-12
Tempoal and Spatial Variations of Water Storage Changes in China with GRACE Data
CHEN Kun1,2,3, JIANG Wei-guo2, HE Fu-hong3, CHEN Zheng2
1. Chinese Research Academy of Environmental, Beijing 100012, China
2. Key Laboratory of Environmental Change and Natural Disaster of Ministry of Education, Faculty of Geographical Science, Beijing Normal University, Beijing 100875, China
3. School of Resources and Environmental Engineering, Ludong University, Yantai 264025, China
Abstract

In this paper, we collected Gravity Recovery and Climate Experiment (GRACE) Tellus products data, GLDAS-Noah surface evapotranspiration data and ground observation rainfall data from 2003 to 2015 in China. The trends of water storage, rainfall and evapotranspiration and the correlations between them were analyzed with Theil-Sen Median trend analysis and Mann-Kendall test. Furthermore, empirical orthogonal function (EOF) decomposition was used to analyze the temporal and spatial distribution characteristics of total water storage in China. The results show that: 1) The annual variation of water reserves in China became more severe after 2006, and the characteristics of abnormal water storage periods can be characterized by precipitation and evapotranspiration. 2) From 2003 to 2015, the changes of water storage, rainfall and evapotranspiration had significant spatial variations. 3) In the past 13 years, climatic factors could explain more than 50% variation of water storage change in almost 18.6% of total area of China, mainly distributing in the Northeast China, the southeast coast area, the Sichuan Basin, the Qinghai Plateau and the northwest of Xinjiang. 4) The principal component information of GRACE water storage decomposed by EOF can characterize the temporal and spatial characteristics of the global and local water storage in China in the study period.

Key words: GRACE; water storage; Theil-Sen Median trend analysis; Mann-Kendall test; EOF

随着气候变化影响的加剧和社会经济的快速发展, 中国水资源短缺和时空分布不均等问题十分突出[1]。因此, 准确地估算中国水储量的时空变化, 分析水储量与气象因素的响应规律, 对中国水资源宏观的开发利用和管理具有重要的参考意义。长期以来, 全球及区域大尺度陆地水储量研究主要是以水文过程模型(如Global Land Data Assimilation System, GLDAS)模拟为主[2], 但是受观测资料和模型参数等因素的影响, 模拟结果不确定性较大。GRACE重力卫星(Gravity Recovery and Climate Experiment)的发射成功, 其采用的近极低轨道星载GPS跟踪、非保守力加速仪以及高精度星间距离测量等新技术, 显著提高了重力场观测的精度和时空分辨率, 从而为大、中尺度区域陆地水储量变化研究提供了新的途径[3]

GRACE卫星的重力观测提供了地球重力场变化信息, 在季节性或年际较短时间尺度上地球重力场变化主要原因是地球表层大气压力、海底压力和陆地水储量变化, 在陆地区域重力场变化主要是由水储量变化造成[4]。在此基础上, Swenson等利用GRACE数据提取陆地水储量变化数据并对结果进行验证, 研究表明在研究区域大于4× 105km2时, GRACE数据估算的水储量变化的精度可达1 cm等效水柱高度[5]。此后, 很多学者利用GRACE数据对全球各大流域和典型区域进行了水储量变化研究, 并验证了GRACE数据适用于大中尺度水储量变化研究的可靠性[6, 7, 8, 9, 10]

在陆地水循环过程中, 降水、蒸散发是影响水储量变化的重要因素[11]。目前, 很多学者利用GRACE数据分析水储量与气象因素之间的关系, Penatti等利用TRMM数据、MOD16数据和GRACE数据分析了潘塔纳湿地水储量与降水和蒸散发之间的相关性, 结果表明该地区水储量变化与降水和蒸散发都有较高的相关性[12]。许民等利用GRACE数据和站点降水数据, 分析了2003— 2010年天山地区水储量变化与月降水量年内和年际变化趋势, 结果表明两者在变化趋势上保持很好的一致性[13]。严家宝等利用GRACE数据和TRMM降水数据对比分析了中国各大流域陆地水储量与降水年内和年际变化规律[14]。以上研究主要从区域尺度分析了降水、蒸散发与水储量变化的相关关系, 然而很少有人从全国角度出发, 分析中国蒸散发和降水共同作用对水储量变化的影响。

因此, 本文利用2003年1月到2015年12月的GRACE陆地水储量数据、GLDAS- Noah的地表蒸散发数据以及地面实测降水数据, 对中国近13 a陆地水储量、降水量和蒸散发量变化进行趋势分析和相关分析, 最后采用经验正交函数(EOF)分解方法分析了中国水储量时空变化特征, 研究结果将有助于我们宏观了解中国陆地水储量与降雨和蒸散发的关系, 并为我国水资源的合理分配与管理提供依据。

1 研究数据
1.1 GRACE重力卫星数据

GRACE重力卫星计划是由美国国家航空航天局(National Aeronautics and Space Administration, NASA)和德国宇航中心(Deutsches Zentrum fü r Luft- und Raumfahrt, DLR)联合开发, 并于2002年3月发射成功。GRACE卫星的主要任务是对地球重力场的时空变化进行精准测量, 获取的时变重力场信息可以等效转换为地球表面的质量变化。

本文采用的陆地水储量变化数据是NASA提供的GRACE Tellus陆地栅格数据产品, 数据时间序列为2003年1月至2015年12月共142期数据, 产品主要包括CSR、JPL和GFZ三个版本, 本文参考前人研究, 采用了CSR机构的产品数据[15, 16]。产品数据反演过程中对重力场球谐系数截断至60阶, 并经过去条带滤波和窗口宽为300 km的高斯滤波处理[17]。最终产品为全球陆地区域1.0° × 1.0° 栅格大小的月等效水柱高度变化数据, 单位为cm。由于卫星传感器问题, 2003— 2015年期间共丢失14期数据, 当数据缺失时, 对缺失月份为中心的5个月窗口用四次多项插值填补。

1.2 降水数据

本文采用的降水数据为中国地面气候资料日值数据集(V3.0)中地面实测数据, 来源于中国气象科学数据共享服务网(http://cdc.nmic.cn/home.do), 该数据集共包含了839个地面气象站数据, 其中国家级基准、基本气象站824个, 其他级别测站15个, 日降水数据的比例因子为0.1, 单位为mm, 本文选取2003年1月到2015年12月的降水日值数据, 通过月累加求和和反距离加权插值处理, 最终得到空间分辨率为1.0° × 1.0° 的全国月降水空间分布数据。

1.3 蒸散发数据

本文采用的2003— 2015年地表蒸散发数据来源于GLDAS-Noah陆面模式的Version 2.1产品, 该产品的时间分辨率为月, 空间分辨率为1° × 1° , 根据研究需要提取了主要的输出变量蒸散发量, 为了与GRACE数据匹配, 将其处理成单位为cm的等效水柱高度, 且对2003年1月至2015年12月的数据做距平处理。

2 分析方法
2.1 经验正交函数

经验正交函数分析方法(Empirical Orthogonal Function, EOF)也称主成分分析或特征向量分析。它主要是利用数据集的方差贡献来识别和强化信号的相似性, 通过将随时间变化的地球物理场分解为正交的时间函数和空间函数部分[18], 一个连续的时空随机场F(x, t)的分解可表示为:

Fx, t=k=1Kcktukx(1)

式中:向量(x, t)表示时间为t、空间位置为x的时空位置; K是正交时空随机场中的模态数, 这些模态由正交空间函数 ukx表示, 即EOFk; 时间关联函数 ck(t)F(x, t)在 ukx上的投影, 又称经验正交函数的膨胀系数, 即PCk。为了解月水储量的时空变化模式, 在EOF分解之前将各栅格的月水储量时间序列数据进行标准化处理。

2.2 Mann-Kendall(M-K)检验法

在水文、气象时间序列趋势分析中, 非参数检验方法比参数检验方法更加适用于非正态分布的数据检验, 非参数M-K检验法常广泛用于降雨、气温和径流等水文气象要素的趋势分析及变异检验, 其结构简单, 计算方便, 且时间序列无需遵从一定的分布, 也不受少数异常值的干扰[19]

定义一组时间序列X=(x1, x2, …, xn), 其标准化的检验统计量公式为:

Z=(S-1VarS(S> 0)0(S=0)(S+1)VarS(S< 0)(2)

式中:

S=i=1n-1j=i+1nsgnxj-xi(3)

sgnxj-xi=1(xj-xi> 0)0(xj-xi=0)-1(xj-xi< 0)(4)

VarS=nn-12n+518(5)

式中:xixj分别表示像元ij时刻的像元值; n表示时间序列的长度; sgn是符号函数; 统计量Z的取值范围为( -+)。在给定显著性水平α 下, 当 Z> u1-α2时, 表示研究序列在α 水平上存在显著的变化。一般取α =0.05, 即本文判断在0.05置信水平上水储量时间序列变化趋势的显著性。

2.3 Theil-Sen Median 趋势分析法

Theil-Sen Median趋势分析是一种稳健的非参数统计的趋势计算方法, 可以减少数据异常值的影响[20]。Theil-Sen Median趋势计算n(n– 1)/2个数据组合的斜率的中位数, 其计算公式为:

SV=MedianVj-Vij-i(2003i< j2015(6)

式中:V代表空间上某点的水储量或降水量时间序列; ij代表时间序列上两个不同的年份。当 SV> 0时, 表示变量呈现增长趋势, 反之, 则表示变量 SV呈现减少趋势。

3 结果与分析
3.1 中国水储量、降水量、蒸散发量时间变化分析

从2003— 2015年全国降水、蒸散发和水储量变化过程可以看出 [图1(a)], GRACE水储量变化与降水量和蒸散发量变化具有明显的一致性, 并表现为明显的季节性变化; GRACE水储量变化滞后于降水量和蒸散发量变化, 且降水量变化的幅度大于水储量和蒸散发量变化的幅度, 其中降水量、蒸散发量和GRACE水储量变化(距平值)幅度分别为12、6.3、10.5 cm, 这在一定程度上反映了降水需要经过地表径流、蒸发、渗透等复杂的水文过程使得水分重新分配, 最终形成陆地水储量[21]

图1 2003— 2015年中国陆地水储量(GRACE)、降水量、蒸散发量变化Fig. 1 Inter-annual changes of GRACE water storage, precipitation and evapotranspiration of China from 2003 to 2015

为了消除水储量、降水量和蒸散发量年内变化的影响, 本文对原时间序列数据做了13点滑动平均, 得到研究时段内中国水储量(GRACE)、降水量和蒸散发量变化图, 如图1(b)所示。降水量、蒸散发量、GRACE水储量三者的年际变化存在明显差异, 降水量虽然与水储量有一定差异, 但是降水量变化能够很好地识别出水储量变化过程中的主要峰谷特征; 从整体上看蒸发量在研究时段内呈上升趋势, 在一定程度上可以解释水储量下降; 2003— 2005年水储量、降水量和蒸散发量的年际变化相对平稳, 2006年以后中国水储量年际波动剧烈, 2006、2007、2009、2011以及2013年以后水储量出现明显下降, 而且多年处于偏低状态, 仅在2010和2012年间回升到平均状态以上, 降水量和蒸散发量的年际变化幅度虽然不大, 但是在水储量年际波动异常年份也能表现出一致的变化趋势。由EM-DAT数据库(① 比利时鲁汶大学灾害流行病学研究中心的紧急灾害数据库(Emergency Events Database, EM-DAT), 该数据在国际灾害管理和研究领域被广泛应用, 该数据库收录了自1900年以来17 000多例大灾害发生和影响数据。)提供的灾害信息可知(图2), 2006年之后全国干旱和洪涝事件频发, 其中2006年、2009— 2011年和2014— 2015年连续干旱可能是直接导致水储量常年偏低的主要原因, 而2011年下半年至2012年水储量的迅速回升可能与降水量增加和各地区连续的洪水事件有关。

图2 2003— 2015年中国地区洪水和干旱发生次数Fig. 2 The number of floods and droughts in China from 2003 to 2015

3.2 中国水储量、降水量、蒸散发量空间变化分析

为了分析中国陆地水储量变化与降水、蒸散发的季节性响应的空间分布规律, 本文计算了研究时段内多年平均的月降水量、月蒸散发量以及月水储量, 并提取了各栅格最大值、最小值出现的月份。由图3(a)、(b)、(d)、(e)可知, 全国大部分地区降水量和蒸散发量最大值主要发生在夏季, 受大陆性季风气候作用, 南方降水量最大值对应的月份要早于北方, 而在高原冰川地区的蒸散发量会受积雪升华影响, 最终导致最大值出现在春季; 从降水量和蒸散发量最小值对应的月份来看, 全国大部分地区降水量和蒸散发量出现在冬季。由图3(c)、(f)可知, 全国各地区水储量最大值和最小值对应的月份具有明显差异, 其中全国最大值出现的月份以长江为分界线由南向北分层明显, 长江以南最大值主要出现在夏季, 长江以北大部分地区出现在秋季, 而华北和东北以及新疆大部分地区分别出现在冬季和春季; 从全国各地区水储量最小值对应的月份来看, 长江以南和黄河上游地区分别出现在冬季和春季, 东北、新疆、内蒙古以北地区和华北地区分别出现在秋季和夏季。综合图3可以看出, 我国大部分地区水储量最大最小值出现的月份要滞后于降水和蒸散发1~2个月, 而对于华北、东北以及新疆大部分地区水储量最大最小值对应的月份明显晚于其他地区, 主要是因为这些地区的水储量变化受人类活动影响和水汽来源不同造成的。前人研究发现, 中国北方地区的水库管理、调水以及运煤等人类活动造成的质量增益在一定程度上可以补偿水储量亏损导致的质量减少[22], 最终导致北方地区水储量最大最小值对应的月份具有延迟性; 东北和新疆地区水储量最大值出现于早春, 主要是因为冬季大部分液态水以固态形式保存, 加上蒸发量和径流量较小, 导致这些地区水储量最大值出现在春汛之前。新疆地区水储量最小值主要出现在秋季, 主要是因为夏季蒸发量较大导致, 而东北地区由于经历春汛和夏汛期导致水储量最小值出现在秋季。吐鲁番盆地由于地形、地势特点形成了与邻近地区截然不同的气候特征, 最终导致该区域内水储量最大最小值对应的月份与周边地区差异明显。

图3 月多年平均降水量、蒸散发量和水储量最大、最小值对应月份空间分布Fig. 3 The months corresponding to the maximum and the minimum monthly average precipitation, evaporation and GRACE water storage

本文利用Theil-Sen Median趋势分析法和M-K检验方法对全国降水量、蒸发量和水储量(GRACE)变化趋势及显著性空间分布特征进行分析。从图4(a)看, 在研究时段内华北平原中部和南部、西藏中部以及横断山脉地区年降水量呈减少趋势, 而中国东北、陕北高原北部、四川盆地、长江中下游平原以及华南部分地区的降水量呈增加趋势, 部分地区降水量年际变化超过1 cm/a, 且通过了95%置信区间的M-K检验; 西北干旱区由于降水量少, 降水量没有明显的变化趋势。图4(b)是蒸散发量年际变化趋势及显著性的空间分布, 我国大部分地区蒸散发呈增加趋势, 其中东北、黄淮平原东部、四川盆地以及华南沿海地区蒸散发量显著增加; 我国西部干旱地区虽然实际蒸散发较小, 但是局部地区仍具有显著变化趋势, 例如新疆西部虽然蒸散发量年增加率不足0.4 cm/a, 但是仍然通过了95%置信区间的M-K检验。从图4(c)可知, 研究时段内我国有6个地区水储量变化显著, 其中东北边境、藏北高原、长江中下游及华南沿海地区水储量整体呈显著上升趋势, 华北地区、西南诸河流域和新疆吐鲁番盆地周边地区水储量呈显著下降趋势, 局部地区虽然水储量整体变化趋势不明显, 但是研究时段内该区域内水储量变化显著。对比图4(a)、(b)、(c)可知, 我国东部和南部地区降水和蒸发量与水储量变化趋势大体一致, 局部地区变化趋势显著, 说明降水和蒸散发在一定程度上影响了该地区的水储量变化, 而西北地区水储量变化与降水和蒸散发的变化趋势截然不同, 这是由于该地区水储量变化主要是受冰川融雪和内流河补给的影响, 近年来的气候变暖导致冰川融化、冻土水分释放增加了高原及周边地区水储量的补给, 加上可可西里和三江源地区生态保护政策的实施促使了水储量增加[23, 24]

图4 2003— 2015年降水量、蒸发量、水储量(GRACE)年变化趋势空间分布Fig. 4 Spatial distribution of annual changes of precipitation, evaporation and GRACE water storage in China from 2003 to 2015

3.3 中国水储量变化对气象因素敏感性分析

本文采用偏相关和复相关分析了2003— 2015年中国陆地水储量变化对气象因素(降水和蒸散发)的敏感性, 即水储量变化与气象因素之间的相关关系。图5(a)、(b)为中国陆地水储量变化与气象因素(降水和蒸散发)的偏相关显著性空间分布, 本文将显著性水平低于0.05的区域认为两者具有显著的线性相关性。从图5(a)看, 全国大部分地区水储量与降水量变化呈正相关, 其中显著正相关区域主要集中在东南沿海、西南中部、黑龙江北部以及新疆西部, 共占全国面积的15.9%, 而显著负相关区域主要在内蒙古高原东北部, 仅占全国面积的1.6%。从图5(b)看, 中国水储量变化与蒸发量变化呈负相关的地区主要分布在中国东部、青藏高原北部以及新疆天山附近, 其中显著负相关区域主要分布于我国东北和青藏高原中部, 共占全国面积的12.5%, 而显著正相关地区较少, 零星分布于山西、陕西、四川、西藏以及新疆等地区。图5(c)为中国陆地水储量变化与气象因素(降水和蒸散发)复相关决定系数的空间分布, 研究时段内全国共有18.6%的地区气候因素对水储量变化的方差解释量超过50%, 这些地区主要集中在中国东北、东南沿海、四川盆地、青海高原以及新疆西北部, 而华北、西北以及西藏北部等地区的方差解释量均小于30%, 说明这些地区水储量变化受气象因素影响较小。

图5 2003— 2015年陆地水储量与降水量和蒸散发量变化的敏感性空间分布Fig. 5 Sensitivity of GRACE water storage to precipitation and evapotranspiration

3.4 中国水储量变化的EOF分析

为了分析水储量变化的空间分布特征, 本文采用2003— 2015年中国GRACE水储量变化数据进行EOF分解, 分解后的第一、二、三特征向量分别占总方差的40.6%、20.1%和11.5%, 前3个分量共占72.2%的比重, 各特征向量的空间分布和时变特征如图6所示。

图6 EOF第一、二、三模态空间分布和时间系数变化Fig. 6 The first, second and third characteristic vector and corresponding time series of GRACE water storage changes

1)第一主成分

图6(a)为中国GRACE水储量变化经EOF分解得到的第一特征向量和时间系数(方差贡献值为40.6%), 表征了全国水储量时空变化的整体趋势。从空间分布来看, EOF1是以强正值为主的西南诸河流域为中心向四周呈递减趋势。在时间系数上, 2003— 2015年整体变化趋势不明显, 年内表现为正负交替, 说明水储量变化具有季节性变化规律; 其中2006、2011和2013年时间系数变化明显减弱, 说明在2006、2011、2013年全国水储量变化减弱, 对比图1的水储量变化曲线, 也可以发现与之对应的规律。

2)第二主成分

图6(b)为中国GRACE水储量变化经EOF分解得到的第二特征向量和时间系数(方差贡献值为20.1%)。从EOF2可以明显看出中国主要有6个水储量变化异常区, 其中华南、可可西里和三江平原地区为强正值区, 西南诸河流域、华北平原和新疆东部吐鲁番盆地等地区为强负值区。在时间系数上, 2003— 2015年呈现明显上升趋势, 表明近13 a华南、可可西里和三江平原地区水储量呈增加趋势, 而西南诸河流域、华北平原和新疆东部吐鲁番盆地等地区水储量呈下降趋势。

3)第三主成分

图6(c)为中国GRACE水储量变化经EOF分解得到的第三特征向量和时间系数(方差贡献值为11.5%)。从EOF3可以看出, 强正值区位于中国东南部和新疆吐鲁番盆地周边区, 而强负值区主要位于云南西南部地区。从时间系数上看, 2006年下半年至2008年呈下降趋势, 表明云南西南部地区水储量呈增加趋势, 东南和新疆地区水储量呈减少趋势; 2009年下半年至2010年上半年时间系数显著上升, 表明云南地区水储量呈减少趋势, 而东南沿海地区水储量增加, 资料显示该时段内云南境内发生了一次干旱事件, 同时湖南、福建东南部地区在该时段内共发生了两次流域性洪水, 因此这种时间系数的异常可能与极端旱涝事件有关; 2011至2015年时间系数变化呈上升趋势, 表明东南和新疆地区水储量增加, 云南等负值区水储量减少。

4 结论

本文利用GRACE水储量变化数据、GLDAS地表蒸散发数据和地面实测降水数据, 对中国2003— 2015年水储量、降水量、蒸散发量变化进行趋势分析和相关分析, 并通过EOF分解对中国水储量进行时空变化分析, 本研究的主要结论如下:

1)GRACE水储量、降水量、蒸散发量变化在时间序列上具有一致性, 呈现明显的季节性变化特征, 但是水储量的季节性变化要滞后于降水和蒸散发变化, 同时GRACE水储量和蒸散发量的变化幅度小于降水量变化幅度。水储量、降水量、蒸散发量年际变化结果表明, 2006年以后水储量年际变化剧烈, 水储量长时期处于减少状态, 同时降水量和蒸散发量在相应时段内变化显著, 而且EM-DAT数据显示该时段内洪涝灾害频发。

2)从中国月水储量、月降水量和月蒸散发量最大、最小值对应月份空间分布来看, 全国降水量和蒸散发量最大值出现在夏季, 最小值出现在冬季; 由于各地区水汽来源和人类活动影响程度的不同, 导致各地区水储量最大最小值出现的月份具有显著差异。从全国水储量、降雨量和蒸散发量的年际变化趋势及显著性空间分布看, 我国各地区水储量变化差异明显, 且全国大部分地区具有显著性变化, 而降水量和蒸散发的显著变化区域主要集中在中国东部地区。

3)从全国水储量变化与气象因素(降水、蒸散发)的偏相关显著性空间分布来看, 降水量与水储量变化呈显著正相关地区主要在华南, 蒸散发量与水储量变化呈显著正相关的区域主要分布在东北和青海高原西部地区。全国水储量变化与气象因素的复相关的决定系数空间分布表明, 全国共有18.6%的地区气象因素对水储量变化的方差解释量超过50%, 主要集中在中国东北、东南沿海、四川盆地、青海高原以及新疆西北部。

4)通过对中国GRACE水储量和降水量变化进行EOF分解, 得到水储量变化的第一、二、三特征向量和时间系数, 前3个特征向量的方差贡献之和达72.2%。第一特征向量的贡献值为40.6%, 因此EOF1中空间分布和时间系数代表整个中国地区水储量的时空变化特征; 第二特征向量的贡献值为20.1%, 结合EOF2的空间变化和时间系数可知, 华南、三江平原和三江源地区水储量呈增加趋势, 西南诸河流域、华北平原和新疆吐鲁番盆地周边地区水储量呈减少趋势; 第三特征向量的贡献值仅为11.5%, 由EOF3的空间分布和时间系数可知, 2010年的时间系数异常波动与云南和东南地区的干旱和洪水有关。

本文的研究结果虽然在一定程度上反映了中国各地区水储量变化对气象因素敏感性差异, 为水资源的管理和分配提供了依据, 但是受数据限制, 本文仅从气象因子(降水和蒸散发)的角度分析了中国水储量变化, 而区域水储量变化是降水、气温、径流、人类活动等多种因素共同作用的结果, 因此未来研究会综合考虑气象、人类活动、径流量等因素对中国水储量变化的影响。

The authors have declared that no competing interests exist.

参考文献
[1] 刘佳骏, 董锁成, 李泽红. 中国水资源承载力综合评价研究[J]. 自然资源学报, 2011, 26(2): 258-269.
[LIU J J, DONG S C, LI Z H. Comprehensive evaluation of China’s water resources carry capacity. Journal of Natural Resources, 2011, 26(2): 258-269. ] [本文引用:1]
[2] RODELL M, HOUSER P R, JAMBOR U, et al. The global land data assimilation system[J]. Bulletin of the American Meteorological Society, 2004, 85(3): 381-394. [本文引用:1]
[3] HU X G, CHEN J, ZHOU Y, et al. Seasonal water storage change of the Yangtze River Basin detected by GRACE[J]. Science in China: Series D Earth Sciences, 2006, 49(5): 483-491. [本文引用:1]
[4] 段建宾, 钟敏, 闫昊明, . 利用重力卫星观测资料解算中国大陆水储量变化[J]. 大地测量与地球动力学, 2007, 27(3): 68-71.
[DUAN J B, ZHONG M, YAN H M, et al. Recovery of land water storage variations in Chinese mainland by use of GRACE data. Journal of Geodesy and Geodynamics, 2007, 27(3): 68-71. ] [本文引用:1]
[5] SWENSON S, WAHR J, MILLY P C D. Estimated accuracies of regional water storge variations inferred from the gravity recovery and climate experiment (GRACE)[J]. Water Resources Research, 2003, 39(8): 1223. [本文引用:1]
[6] LINAGE C D, FAMIGLIETTI J S, RANDERSON J T. Statistical prediction of terrestrial water storage changes in the Amazon Basin using tropical Pacific and North Atlantic sea surface temperature anomalies[J]. Hydrology and Earth System Sciences, 2014, 18(6): 2089-2102. [本文引用:1]
[7] CHEN J L, LI J, ZHANG Z Z, et al. Long-term groundwater variations in Northwest India from satellite gravity measurements[J]. Global and Planetary Change, 2014, 116: 130-138. [本文引用:1]
[8] AWANGE J L, KHANDU, SCHUMACHER M, et al. Exploring hydro-meteorological drought patterns over the Greater Horn of Africa (1979-2014) using remote sensing and reanalysis products[J]. Advances in Water Resources, 2016, 94: 45-59. [本文引用:1]
[9] ABELEN S, SEITZ F, ABARCA-DEL-RIO R, et al. Droughts and floods in the La Plata basin in soil moisture data and GRACE[J]. Remote Sensing, 2015, 7(6): 7324-7349. [本文引用:1]
[10] 杨元德, 鄂栋臣, 晁定波, . GRACE估算陆地水储量季节和年际变化[J]. 地球物理学报, 2009, 52(12): 2987-2992.
[YANG Y D, E D C, CHAO D B, et al. Seasonal and inter-annual change in land water storage from GRACE. Chinese Journal of Geophysics, 2009, 52(12): 2987-2992. ] [本文引用:1]
[11] 金晓龙, 邵华, 张弛, . GPM卫星降水数据在天山山区的适用性分析[J]. 自然资源学报, 2016, 31(12): 2074-2085.
[JIN X L, SHAO H, ZHANG C, et al. The applicability evaluation of three satellite products in Tianshan Mountains. Journal of Natural Resources, 2016, 31(12): 2074-2085. ] [本文引用:1]
[12] PENATTI N C, ALMEIDA T I R D, FERREIRA L G, et al. Satellite-based hydrological dynamics of the world’s largest continuous wetland [J]. Remote Sensing of Environment, 2015, 170(1): 1-13. [本文引用:1]
[13] 许民, 叶柏生, 赵求东, . 利用GRACE重力卫星监测新疆天山山区水储量时空变化[J]. 干旱区研究, 2013, 30(3): 404-411.
[XU M, YE B S, ZHAO Q D, et al. Spatiotemporal change of water reserves in the Tianshan Mountains, Xinjiang based on GRACE. Arid Zone Research, 2013, 30(3): 404-411. ] [本文引用:1]
[14] 严家宝, 贾绍凤, 吕爱锋, . 近十年中国陆地水储量变化及其时空分布规律[J]. 南水北调与水利科技, 2016, 14(4): 21-28.
[YAN J B, JIA S F, LÜ A F, et al. Spatial- temporal variation characteristics of China terrestrial water storage in the last ten year. South-to-North Water Transfers and Water Science & Technology, 2016, 14(4): 21-28. ] [本文引用:1]
[15] 曹艳萍, 南卓铜, 胡兴林. 利用GRACE重力卫星数据反演黑河流域地下水变化[J]. 冰川冻土, 2012, 34(3): 680-689.
[CAO Y P, NAN Z T, HU X L. Changes of groundwater storage in the Heihe River Basin derived from GRACE gravity satellite data. Journal of Glaciology and Geocryology, 2012, 34(3): 680-689. ] [本文引用:1]
[16] MO X, WU J J, WANG Q, et al. Variations in water storage in China over recent decades from GRACE observations and GLDAS[J]. Natural Hazards and Earth System Sciences, 2016, 16(2): 469-482. [本文引用:1]
[17] SWENSON S, WAHR J. Post-processing removal of correlated errors in GRACE data[J]. Geophysical Research Letters, 2006, 33(8): 1223. doi: DOI:10.1029/2005GL025285. [本文引用:1]
[18] 邹正波, 李辉, 吴云龙, . 日本M_W9. 0地震震区及其周缘2002—2015年卫星重力变化时空特征[J]. 地震学报, 2016, 38(3): 417-428.
[ZOU Z B, LI H, WU Y L, et al. Spatial and temporal characteristics of long-term satellite gravity change in the epicenter of Mw9. 0 Japan earthquake and its surrounding region. Acta Seismologica Sinica, 2016, 38(3): 417-428. ] [本文引用:1]
[19] 刘兆飞, 王翊晨, 姚治君, . 太湖流域降水、气温与径流变化趋势及周期分析[J]. 自然资源学报, 2011, 26(9): 1575-1584.
[LIU Z F, WANG Y C, YAO Z J, et al. Trend and periodicity of precipitation, air temperature and runoff in the Taihu Lake Basin. Journal of Natural Resources, 2011, 26(9): 1575-1584. ] [本文引用:1]
[20] 袁丽华, 蒋卫国, 申文明, . 2000—2010年黄河流域植被覆盖的时空变化[J]. 生态学报, 2013, 33(24): 7798-7806.
[YUAN L H, JIANG W G, SHEN W M, et al. The spatio-temporal variations of vegetation cover in the Yellow River Basin from 2000 to 2010. Acta Ecologica Sinica, 2013, 33(24): 7798-7806. ] [本文引用:1]
[21] 魏浩翰, 何秀凤, 郑加柱. 联合重力卫星数据及气候水文模型监测西南地区陆地水储量变化[J]. 河海大学学报(自然科学版), 2013, 41(6): 488-492.
[WEI H H, HE X F, ZHENG J Z. Terrestrial water storage variations in Southwest China revealed by gravity mission and hydrologic and climate model. Journal of Hohai University (Natural Sciences Edition), 2013, 41(6): 488-492. ] [本文引用:1]
[22] TANG Q H, ZHANG X J, TANG Y. Anthropogenic impacts on mass change in North China[J]. Geophysical Research Letters, 2013, 40(15): 3924-3928. [本文引用:1]
[23] YAO X J, LIU S Y, LONG L, et al. Spatial-temporal characteristics of lake area variations in Hoh Xil region from 1970 to 2011[J]. Journal of Geographical Sciences, 2014, 24(4): 689-702. [本文引用:1]
[24] XIANG L W, WANG H S, STEFFEN H, et al. Groundwater storage changes in the Tibetan Plateau and adjacent areas revealed from GRACE satellite gravity data[J]. Earth and Planetary Science Letters, 2016, 449: 228-239. [本文引用:1]