李晓婧, 白艳萍, 李 萌, 马金辉

(兰州大学 资源环境学院， 甘肃 兰州 730000)

水资源是我国西北内陆地区社会经济发展和生态平衡的关键制约因素[1]。在气候变化及人类活动影响不断加剧的背景下，从区域到全球尺度监测陆地总水储量及各要素的变化对人类生存和可持续发展具有重要意义[5]。重力反演与气候试验卫星GRACE(gravity recovery and climate experiment)广泛地被应用与中长空间尺度流域的水储量变化研究[2]。地球各系统间的质量重新分布和迁移引起地球重力场的变化，在年或季较短时间尺度上，主要由地球表层大气、海洋、陆地水质量交换引起。通过精确观测地球重力场随时间的变化，可以反演地球水质量变化的情况[3]。目前，利用重力卫星重建地球重力场取得了丰富的成果[4-9]，是传统陆地水文观测技术的有效补充[5]。河西走廊位于甘肃省西北部，行政上属于甘肃省张掖、酒泉、武威、金昌、嘉峪关5市。研究区气候干旱，属于大陆性气候区，年降水量仅为40～400 mm，年蒸发量却高达1 500～3 000 mm。该地区地势南高北低，南部祁连山地的冰川储量丰富[10]，是河西地区降水较多的地区，也是河西走廊石羊河、黑河和疏勒河3大内陆河的发源地。河西走廊干旱缺水，生态环境极为脆弱，是制约河西走廊地区社会经济可持续发展的主要瓶颈[11]。加之人类活动对水资源的过度开发利用，造成了天然绿洲退缩、地下水位下降、内陆河断流、沙漠化扩大等严重的生态问题[12]。本文对河西走廊水资源和生态环境的时空变化规律、以及时空关联关系进行分析，旨在为干旱区水资源合理开发及实现生产、生活和生态可持续发展提供科学依据。

1 研究资料及研究方法

1.1 数据来源

国际上应用解算过的GRACE 卫星数据来反演全球以及区域的陆地水储量的变化，通行的做法[13]有：一种是应用已解算的月重力场模型数据直接求解；另一种是利用卫星轨道数据推求区域陆地水储量变化，即Mascon方法。本文GRACE数据采用空间研究中心(CSR)提供的2003—2015年的RL05 GRACE MASCON解决方案(http:∥www.csr.utexas.edu/grace)，空间分辨率为0.5°，时间分辨率为1个月，栅格值表示等效水高。由于GRACE部分月度数据缺失，采用缺失月份多年平均值表示该月的等效水高。降水采用NASA网站提供的2003—2015年TRMM 3B43数据(https:∥pmm.nasa.gov)，NDVI(归一化差分植被指数)数据来源于美国国家航空航天局(NASA)提供的MOD13产品。2000，2005，2010，2015年的土地利用/覆盖类型数据由中国科学院资源环境科学数据中心(http:∥www.resdc.cn)提供。

1.2 研究方法

1.2.1 趋势分析 采用一元线性回归趋势分析法模拟每个栅格的变化趋势，计算公式为：

式中：i——时间段序号；n——监测时间段累计次数； rasteri——第i时间段栅格图层。

1.2.2 相关性分析 计算两个时间序列x(t)和y(t)相关系数通常用以下公式：

2 研究结果

2.1 水储量时空变化

利用一元线性回归分析法，计算研究区每个像元2003—2015年水储量变化趋势，结果如图1所示。由图1可知，研究区祁连山地区水储量呈增长趋势，等效水柱每年最大增长3.1 mm，沿着径流方向水储量增长趋势愈弱，直至为减少趋势，等效水柱每年最大减少1.9 mm。

图1 河西走廊GRACE等效水高变化趋势空间分布

通过区域统计得到河西走廊及河西走廊内3大流域GRACE等效水高月度变化折线图(图2)。结果显示研究区水储量变化在月尺度上表现出明显的周期性，波动范围基本保持在-20～20 mm之间，周期为12个月，表示研究区整体及3大流域水储量变化存在显著的年周期；疏勒河流域与黑河流域水储量变化在年内波动变化的最高值和最低值出现月份一致。经过12次滑动平均及线性回归，可观察到研究区及3大流域水储量变化均呈下降趋势。通过观察可以发现疏勒河流域2003—2009年水储量变化趋势平稳，2009年后水储量在波动中下降，2011—2015年水储量年内变化波动范围与极大值逐年增大；黑河流域2003—2010年水储量稳定下降，2010—2015年水储量在波动中下降且年内变化波动范围较大，年内变化的极大值呈增长趋势；石羊河流域水储量呈逐年下降趋势，年内变化波动范围较稳定。研究区3大流域中石羊河流域水储量降低速度最快，约为0.146 3 mm/月，黑河流域水储量降低速度较快，约为0.134 7 mm/月，疏勒河流域水储量下降速度较慢，约为0.085 mm/月，研究区整体水储量呈较快的下降趋势，约为0.109 2 mm/月。

图2 研究区水储量变化

利用2002—2015年河西走廊整体及其3大流域月尺度水储量变化值多年平均值，分析研究区年内水储量变化情况(图3)。疏勒河流域水储量在4月份进入累积状态，7月份达到最高值，然后开始逐渐下降，12月降至最低值。黑河流域水储量在1至7月缓慢上升，7月份进入累积状态，同时达到全年最高值，然后缓慢下降，在10月下降至亏损状态。石羊河流域水储量1至9月在波动中上升，7月进入累积状态，9月达到全年最高，11月下降至亏损状态。疏勒河流域水储量年内波动幅度最大，石羊河流域水储量年内变化相对于疏勒河流域及黑河流域呈现一定的滞后效应，晚2个月达到全年最高值。

2.2 水储量变化与降水的关系

通过区域统计得到2003—2015年河西走廊及3个流域分区的历年年降水量，结果显示研究区及3个流域分区的年降水总体趋势平稳，线性趋势斜率接近于0，黑河流域年降水年际变化波动范围较大。利用2003—2015年河西走廊整体及其3大流域月尺度降水数据，分别计算各月累年均值，结果如图4所示。研究区及其3大流域在冬季降水量最少，春季降水开始增多，夏季降水量最大，秋季降水量开始减少。

图3 研究区各月等效水高累积多年平均值

图4 研究区年降水变化趋势和各月降水累积多年平均值

通过基于像元通过对多年GRACE等效水高与降水进行相关性分析(图5)，河西走廊南部的祁连山地年水储量与年降水为正相关，研究区祁连山地以北年水储量与年降水为负相关，是由于河西走廊3大流域水源补给主要来源为祁连山冰雪融水，研究区流域中下游干旱少雨，潜在蒸散发量大于降水量；另一方面，由甘肃省水资源公报可知2003—2015年研究区中游及下游的农田有效灌溉面积呈增加趋势，水资源的补给和消耗除了受降水量影响外，还深受地表径流和地下水的补给、径流和排泄以及人类生产、生活和生态用水等因素影响。因此，研究区除祁连山地区以外，水储量的变化受降水影响较小。

图5 河西走廊年水储量变化与年降水相关系数分布

利用区域统计得到的研究区2003—2015年水储量变化值与降水量，进行月尺度与年尺度相关性分析(表1)。月尺度下，研究区整体的水储量变化与降水正相关性最好，相关系数为0.62，且置信水平为0.01，其次为疏勒河流域、黑河流域，石羊河流域的水储量变化与降水在月尺度下相关系数最低。年尺度下，研究区整体及3大流域分区的水储量变化与降水均不存在相关性。

2.3 水储量变化与生态环境时空耦合

河西走廊地区在“移民政策”实施以后，农业生产大大增加了研究区内水资源的消耗，进一步加剧了研究区内荒漠化的状况[14]。由中国科学院资源环境科学数据中心(http:∥www.resdc.cn)提供2000，2005，2010和2015年的土地利用/覆盖类型栅格数据，通过区域统计得到河西走廊三大流域2000，2005，2010以及2015年土地利用/覆盖类型面积比重(表2)。疏勒河流域土地利用/覆盖类型主要为未利用地(约79%)和草地(约18%)，林地和草地无明显变化，耕地和建筑用地增加明显，冰川与未利用地面积比重下降，水域面积比重轻微上升。黑河流域土地利用/覆盖类型主要为未利用地(约57%)、草地(约24%)与耕地(约10%)，耕地与建筑用地面积比重明显上升，水域面积比重先下降后上升，冰川面积比重在2015年出现上升。石羊河流域土地利用/覆盖类型主要为未利用地(约44%)、草地(约29%)与耕地(约18%)，耕地与建筑用地面积比重明显上升，水域面积比重轻微上升，草地面积比重出现下降。疏勒河流域冰川呈负增长趋势，但结合前文可知疏勒河流域南部的祁连山地水储量呈正增长趋势，原因可能为冻土退化，随着冻土层的隔水作用减小，活动层加大和加厚，使地下水库库容增大，更多的地表水入渗变成地下水，造成地下水水库的储水量增加[15]。

表1 水储量变化与降水的年度和月度相关性

注：**表示显著性水平为0.01； *表示显著性水平为0.05。下同。

植被是干旱区衡量生态环境质量高低及演变规律最有效的指示器[16]。通过区域统计得到2003—2015年河西走廊及三个流域分区的历年年均NDVI，2009年前NDVI年均值趋势较平稳，增长速度较慢，2009年以后增长速度较快(图6)。通过NDVI各月累年均值可知，石羊河变化幅度最高，其次为黑河流域，疏勒河流域NDVI波动变化较小，年内一直保持在一个比较低的水平，3个流域分区NDVI最大值出现在7月份，11月份回落到较低水平。研究区整体及3个流域分区植被生长的物候特征比较相似，植被在春季(3月至5月)进入生长期，夏季(6月至8月)达到旺盛期，秋季(9月至11月)开始衰退，冬季(12月至次年2月)进入休眠期。

表2 河西走廊各流域2000-2015年土地利用/覆盖类型面积比重

注：标“↑”符号表示面积比重增加明显。

图6 研究区年均NDVI变化趋势和各月NDVI累积多年平均值

统计得到的研究区2003—2015年水储量变化值与NDVI，进行月尺度与年尺度相关性分析(表3)。月尺度下，研究区整体的水储量变化与降水相关性较好，相关系数为0.49，置信水平为0.01，其次为黑河流域，疏勒河流域与石羊河流域的水储量变化与NDVI在月尺度下相关系数最低。年尺度下，研究区整体及石羊河流域与黑河流域的水储量变化与NDVI均存在显著的高度负相关，置信水平为0.01，疏勒河流域水储量变化与NDVI存在较显著的低度负相关。

表3 水储量变化与NDVI的年度和月度相关性

通过对空间尺度进一步细化，深入研究水储量与植被之间的耦合关系。对研究区像元尺度的水储量年际变化量与NDVI年均值做相关分析(图7)，发现主要灌区水储量年际变化量与植被指数年均值呈现高度负相关，主要原因为研究时间段内研究区耕地面积的持续增加，大量地消耗了该区域的地下水资源的同时增加了地面水分蒸发量，导致水资源流失，最终表现为水储量年际变化量与植被指数呈明显负相关的特征。疏勒河流域南部山区呈正相关，原因为该区高寒草甸生长的主要限制因素为温度[17]，全球气候变暖有利于高寒草甸植被增长[18]，进而蓄水功能得到增加。

图7 水储量年际变化与NDVI相关性空间分布

3 讨论与结论

河西走廊及3个流域分区的水储量变化周期为12个月，年尺度上水储量均呈现下降趋势，其中石羊河流域与黑河流域水储量下降较快，疏勒河流域下降较慢。由于生态建设工程的实施[19]，研究区水储量下降速度趋于放缓。从多年月平均尺度看，夏季水储量为积累期，冬季为亏损期，由于河西走廊3个流域分区的自然条件与社会经济发展状况等方面有所差异[20]，水储量年内变化规律存在差异。

河西走廊及3个流域分区2003—2015年历年降水量趋势平稳，疏勒河流域与黑河流域的降水与水储量变化的月度累年均值趋势一致，由于灌溉用水需求量大，石羊河流域水储量变化月度累年均值最大值出现的月份比降水月度累年均值的最大值晚一个月。月尺度上降水量与水储量变化呈现较好的正相关性，年尺度上降水量与水储量变化不存在相关性。河川径流是干旱区内陆河流域的主要水资源[21]，由于研究区3个流域分区的中下游水资源的补给与消耗过程比较复杂，例如上游来水对中下游水储量的影响，因此基于像元的水储量变化与降水量的相关性不显著。

2003—2015年在流域降水量保持相对稳定的背景下，研究区及流域分区的水储量均不断下降，而植被指数则在持续上升，在流域尺度和整体上均呈现高度负相关。河西走廊3个流域分区中石羊河流域NDVI累年均值最大，并且水储量月度累年均值出现的最大值的月份比NDVI月度累年均值最大值出现的月份滞后2个月，说明植被影响了水储量的累积。自然条件与人类活动共同作用于水循环系统，使水资源在人工绿洲和天然绿洲、流域上中下游之间重新进行分配。河西走廊人工绿洲规模的扩大主要是人类改造自然绿洲或未利用土地获得的[21]，并且人工绿洲的演变主要跟人类的经济活动有关。从像元尺度上看，灌区内水储量年际变化量与植被指数年均值呈高度负相关，南部祁连山地区水储量年际变化量与植被指数呈正相关。随着研究区耕地面积的不断扩大，农业生产对水资源的消耗量也大幅增加，导致该区水储量逐年减少，因此研究区主要灌区水储量年际变化与植被指数年均值表现为显著的负相关关系。

综上，2003—2015年研究区水储量变化历年呈下降趋势，并在流域尺度及像元尺度表现出明显的空间异质性和尺度效应。研究时间段内研究区降水保持相对稳定，没有明显的变化趋势，证明降水不是水储量发生下降的主要原因。从土地利用/覆被类型变化和2003—2015年历年NDVI均值变化来看，河西走廊人工植被呈增加趋势。研究区水储量的下降受人类活动影响较大，主要表现在耕地面积和建筑面积比率的增长，从而增加了农业生产及人类生活耗水量，并直观地表现为水储量年际变化与植被指数年均值显著负相关，石羊河流域表现最为明显。由于GRACE数据空间分辨率较低，文章仅从宏观角度对研究区水储量变化进行了分析，后期可结合径流、土壤含水率与地下水实测数据开展进一步研究。