中国生态农业学报(中英文)  2021, Vol. 29 Issue (6): 1030-1041  DOI: 10.13930/j.cnki.cjea.200806
0

引用本文 

石嘉丽, 张晓龙, 李红军, 沈彦军. 河北坝上地区绿水时空变化及其驱动因素研究[J]. 中国生态农业学报(中英文), 2021, 29(6): 1030-1041. DOI: 10.13930/j.cnki.cjea.200806
SHI J L, ZHANG X L, LI H J, SHEN Y J. Spatial-temporal changes in green water and its driving factors in the Bashang area of Hebei Province[J]. Chinese Journal of Eco-Agriculture, 2021, 29(6): 1030-1041. DOI: 10.13930/j.cnki.cjea.200806

基金项目

国家自然科学基金项目(41807177,42001037,41971262)资助

通信作者

沈彦军, 主要从事流域生态水文模拟与水环境管理方面的研究。E-mail: shenyanjun@sjziam.ac.cn

作者简介

石嘉丽, 主要从事生态水文模拟方面的研究。E-mail: shijiali19@mails.ucas.ac.cn

文章历史

收稿日期:2020-10-10
接受日期:2020-12-29
河北坝上地区绿水时空变化及其驱动因素研究*
石嘉丽1,2, 张晓龙1, 李红军1, 沈彦军1,2     
1. 中国科学院遗传与发育生物学研究所农业资源研究中心/中国科学院农业水资源重点实验室/河北省节水农业重点实验室 石家庄 050022;
2. 中国科学院大学 北京 100049
摘要:河北坝上半干旱农牧交错带是首都西北重要的生态屏障。绿水在维持坝上半干旱生态系统稳定、保障京津冀地区生态环境安全中发挥着重要作用,但绿水资源的时空分布及其驱动因素还缺乏定量理解。本文基于GLASS蒸散产品,定量分析了2001—2015年河北坝上(康保县、沽源县、尚义县、张北县)县域尺度的绿水时空变化特征,通过线性变化系数、相关系数、贡献率方法评估了绿水变化对气温、降水、净辐射(Rn)、总第一性生产力(GPP)以及土地利用的响应规律,以期为坝上地区“首都水源涵养功能区和生态环境支撑区”建设、水资源高效利用提供科学依据。研究结果表明:1)2001—2015年坝上4县的绿水量总体呈下降趋势,多年平均绿水量为371.11 mm;绿水量季节变化明显,表现为夏季>春季>秋季>冬季。2)绿水量由西北至东南递增,其中沽源县绿水量最大,康保县绿水量最小;研究期间21.2%的区域绿水量呈增加趋势,78.8%的区域呈减少趋势。3)绿水与降水、GPP呈正相关关系,与气温和净辐射呈负相关关系,各影响因子对绿水变化的贡献率排序为GPP >气温>降水> Rn,GPP的贡献率高达51%。4)研究区内各土地利用类型绿水量表现为林地>草地>耕地>建设用地。绿水量变化趋势受土地利用变化影响显著,土地利用更直接地影响绿水的空间分布。
关键词绿水    蒸散    气候因素    总第一性生产力(GPP)    土地利用    河北坝上    
Spatial-temporal changes in green water and its driving factors in the Bashang area of Hebei Province*
SHI Jiali1,2, ZHANG Xiaolong1, LI Hongjun1, SHEN Yanjun1,2     
1. Center for Agricultural Resources Research, Institute of Genetics and Developmental Biology, Chinese Academy of Sciences/Key Laboratory of Agricultural Water Resources, Chinese Academy of Sciences/Hebei Key Laboratory of Water-saving Agriculture, Shijiazhuang 050022, China;
2. University of Chinese Academy of Sciences, Beijing 100049, China
Abstract: The Bashang area is a semi-arid farming-pastoral zone in Hebei Province and represents an important ecological barrier northwest of the capital. Green water plays an important role in maintaining the stability of the semi-arid ecosystem in Bashang and ensuring the safety of the ecological environment in the Beijing-Tianjin-Hebei region. However, a quantitative understanding of the temporal and spatial distribution of green water in this area and underlying driving factors is lacking. This study used Global Land Surface Satellite (GLASS) evapotranspiration products to quantitatively analyze the spatio-temporal changes in green water and evaluated the effect on changes in green water of temperature, precipitation, net radiation (Rn), gross primary productivity (GPP), and land use at the county level (Kangbao, Guyuan, Shangyi, and Zhangbei Counties) from 2001 to 2015. These results provide a scientific basis for the construction of the "Capital Water Conservation Functional Zone and Ecological Environment Support Zone" and more efficient use of the water resources in the Bashang area. This study used the linear trend method to assess the changing trends of green water, and based on the pixel and correlation coefficient methods to analyze the driving factors of green water variation. Multiple regression analysis was used to calculate the contribution rate of each driving factor. The results showed that: 1) the amount of green water in the four counties showed an overall insignificant downward trend from 2001 to 2015. The maximum value was in 2003 (415.34 mm), and the minimum value was in 2009 (322.35 mm). The annual average green water amount was 371.11 mm. The seasonal variation of green water amount was obvious. Summer had the highest amount of green water, followed by spring, autumn, and winter. 2) The amount of green water increased from northwest to southeast regions, most of which showed a decreasing trend. Guyuan County had the highest amount of green water, and Kangbao County had the least amount. Across the regions, 21.2% showed an increasing trend in the volume changes of green water, and 78.8% of the regions showed a decreasing trend from 2001 to 2015. 3) From 2001 to 2015, the precipitation, temperature, and GPP showed an overall increasing trend, whereas Rn showed an insignificant downward trend. Green water was positively correlated with precipitation and GPP and negatively correlated with temperature and Rn. The contribution rate of each influencing factor was in the following order: GPP > temperature > precipitation > Rn. The GPP contribution rate was as high as 51%. The above-mentioned factors primarily affected temporal changes in green water. 4) The highest amount of green water was in the forest land, followed by grassland, cultivated land, construction land, and unused land. There were no significant differences in the amount of green water among land-use types, but the trend of green water changes was significantly affected by land-use changes. Land use directly affected the spatial distribution of green water.
Keywords: Green water    Evapotranspiration    Climate factors    Gross Primary Productivity (GPP)    Land use    Bashang area of Hebei Province    

“绿水”的概念最早由瑞典水文学家Falkenmark[1]提出, 是指在热力作用和植物根系吸水作用下由降水转化而来的水资源[2], 包括绿水流和绿水蓄量。绿水流主要为地表实际蒸散发, 满足了世界83%的雨养农业用水需求[3], 是植物生长的必要条件。作为水量平衡和热量平衡的纽带, 绿水流的时空变化深刻地影响着生态系统的协调均衡。绿水蓄量主要为土壤水分变化量[4], 直接影响土壤水分传输、土壤蒸发、植物根系吸水等过程。一般情况下, 绿水蓄量的年际变化可忽略不计。开展绿水时空动态变化研究可以为理解植被耗水、水文过程变化及水资源合理利用提供新的思路。

绿水水量计算方法主要包括水量平衡法和水文模型模拟法[5-6]。水量平衡法是通过计算区域降水量与径流量(包括地表径流与地下径流)之差得到区域绿水量[5-7], 因其原理易于理解, 可操作性强而广受青睐。水文模型法对水文过程的反映详尽, 但对驱动数据要求较高, 机理相对复杂, 在实际应用中存在一定局限性。应用较广的水文模型有LPJ (Lund-Potsdam-Jena)[6]、SWAT (Soil and Water Assessment Tool)[8]、HIMS-VIH (Hydro-Informatic Modeling System-Vegetation Impacts on Hydrology)[9]、MIKE SHE[10]等模型。众多学者已在多种尺度下开展了绿水资源的研究, 如从小流域尺度[5]到区域[8]再到全球尺度[6], 日尺度[10]到年际尺度[8]再到年代尺度[7], 这些研究从不同角度揭示了绿水多时空尺度的变化过程。在气候变化和人类活动对水循环影响日益深刻的背景下, 探究绿水变化的驱动因素成为绿水研究的热点之一。Zang等[11]定量衡量了土地利用变化对黑河流域蓝绿水转化的影响, 发现土地利用变化使2.06亿m3的蓝水转化为绿水; 赵安周等[12]探讨了气候变化和人类活动对黄河支流渭河流域蓝绿水资源的影响, 发现人类活动对绿水量变化影响较小, 气候变化是导致研究区绿水量变化的主要原因; 冯畅等[13]评估了变化环境对涟水流域蓝绿水的影响, 发现在时间变化特征上气候变化使绿水流呈一定的上升趋势而绿水蓄量呈持续下降趋势, 在空间分布上绿水资源与土地利用空间分布格局表现出较显著的相关性。这些研究对了解当地绿水变化机理具有重要贡献, 但普遍对蒸散发过程考虑不足, 气象站点驱动数据在从站点到区域的尺度转换中也存在一定不确定性。

不断发展的遥感技术为绿水研究提供了更加便利的条件, 可以避免数据资料不足、尺度转化造成的问题, 已成为当前研究绿水流的主要手段。其中MODIS (Moderate Resolution Imaging Spectroradiometer)蒸散产品在全球各区域得到了广泛应用[14-17], 由北京师范大学开发的GLASS (Global Land Surface Satellite)蒸散产品[18]在MODIS数据集的基础上, 利用贝叶斯方法改进的PM、PT-JPL、MS-PT及半经验彭曼算法计算得到, 以通量站观测数据为参考值确定每种算法的权重, 其数据精度得到较大提升, 具有良好的应用前景。

河北张家口坝上地区为典型的农牧交错地带, 属于内流区, 肩负着构建京津冀绿色生态屏障的任务, 具有防风固沙和涵养水源等生态功能[19]。但该地区生态环境脆弱, 处理好水资源开发与环境保护的关系, 对实现区域经济、社会、环境协调发展和保障京津冀地区生态环境安全具有重要意义。气候变化背景下, 以退耕还林、还草为主体的国家生态建设工程极大地改变了下垫面状况, 进而直接影响该地区绿水的时空演变过程。作为水热平衡的重要组分, 绿水资源的时空变化规律及其驱动因素尚不明确, 开展该研究是科学进行区域水资源管理及合理规划的前提, 也是“首都水源涵养功能区和生态环境支撑区”建设的迫切需求。因此, 本研究选择河北张家口坝上地区为研究区, 基于广泛使用的蒸散产品、气象再分析数据产品、土地利用、地表净辐射等遥感产品, 分析该地区绿水时空变化特征及其主要驱动因素, 以期为坝上地区“两区”建设、水资源高效利用提供科学依据。

1 研究区概况与数据来源 1.1 研究区概况

“坝上”由地形陡然升高而形成, 位于华北平原向内蒙古高原过渡的地区。张家口坝上地区包括张北县、沽源县、康保县和尚义县, 位于113.81°~ 116.06°E, 40.73°~42.17°N, 海拔多在1200~1800 m, 南高北低(图 1)。该区域为东亚大陆性季风气候, 是半干旱区向干旱区的过渡带, 年平均气温−0.3~3.5 ℃, 年均降水量400 mm左右[20], 集中在6—8月。特殊的地形和气候条件使该地区生态系统相对脆弱, 自我修复能力差。主要土地利用类型为耕地和草地, 产流微弱, 绿水是降水量消耗的主要途径。近20年来高耗水作物种植规模不断扩大, 过度的灌溉活动使地表水断流、湿地湖泊萎缩, 水资源形势愈加严峻。

图 1 张家口市坝上4县研究区地理位置 Fig. 1 Location of the study area of the four counties of Zhangjiakou in Bashang area
1.2 数据来源与处理

实际蒸散发(ET)数据与总第一性生产力(GPP)数据来自于北京师范大学开发完成的GLASS数据集[18]。气温、降雨数据来源于中国科学院青藏高原研究所开发的中国区域地面气象要素驱动数据集[21-22]。净辐射(Rn)数据来源于张晓龙[23]估算的1982—2017年中国大陆区域逐日Rn。土壤水分数据来自中国科学院青藏高原研究所开发的基于微波数据同化的中国土壤水分数据集[24-27]。土地利用数据来源于中国科学院资源环境与数据中心1∶100万土地利用数据产品。数据详细信息见表 1

表 1 研究所用数据类型、来源及时空分辨率 Table 1 Data source and spatio-temporal resolution used in the study

ET数据为8 d数据, 对其进行格式、投影转化, 影像裁切处理, 再将8 d数据线性插值得到每日数据, 累加求和得到多年月均数据及年尺度数据。年均降雨量、Rn及气温由日值数据累加及求均值得到, 最后将以上数据重采样统一为0.01°。将2000年、2010年、2015年3期土地利用数据重分类为耕地、林地、草地、建设用地、水体、未利用地6类, 并结合野外实测数据对错分类别进行修正, 将栅格格式转化为矢量后分别对绿水量进行提取。本研究除土壤水分数据以外, 其他所有研究数据研究时段均为2001—2015年, 受数据获取限制, 土壤水分数据研究时段仅为2002—2011年。

2 研究方法 2.1 绿水蓄量计算方法

绿水蓄量即土壤水分含量变化量, 计算公式如下[28]:

$ \Delta W = h \times \Delta \theta \times 10/100 $ (1)

式中: $\Delta W$为土壤水分含量变化量即绿水蓄量(mm); h为土层厚度(cm); Δθ为土壤体积含水量变化量(%)。土壤水分含量变化量分为0~5 cm、5~20 cm、20~100 cm 3层计算。

2.2 贡献率计算方法

采用多元回归分析法计算各因子对绿水变化的贡献率。为衡量各因子对绿水量相对影响程度, 需对研究数据进行标准化处理, 使之介于0~1。对标准化后的数据进行多元回归分析, 计算公式如下[29]:

${Y_{\rm{s}}} = \alpha {X_{1{\rm{s}}}} + b{X_{2{\rm{s}}}} + c{X_{3{\rm{s}}}} + \cdots $ (2)
$ {\eta _1} = |a|/\left( {|a| + |b| + |c| + \cdots } \right) $ (3)

式中: Ys为因变量变化标准值, X1sX2sX3s$ \cdots $为自变量标准化值, abc为各因子标准化后回归系数, η1X1变化对于Y变化的相对贡献率。

2.3 评估指标

用最小二乘法基于像元拟合得到线性变化趋势线, 趋势线的斜率即为线性变化系数θslope, 以此评估绿水量及各影响因子的变化趋势。在IDL的环境实现计算, 计算公式为[30]:

${\theta _{{\rm{slope}}}} = \frac{{n \times \sum\limits_{i = 1}^n {\left( {i \times {Y_i}} \right) - \sum\limits_{i = 1}^n {i\sum\limits_{i = 1}^n {{Y_i}} } } }}{{n \times \sum\limits_{i = 1}^n {{i^2} - {{\left( {\sum\limits_{i = 1}^n i } \right)}^2}} }}$ (4)

式中: θslope为线性变化系数, i为从2001年开始的第i年, n为研究年数15, Yi为第i年的变量值。利用绿水量及各影响因子年值序列和时间序列的相关关系来判断年际间变化的趋势, θslope < 0表示变化趋势为减少, θslope > 0表示变化趋势为增加。

利用皮尔逊相关系数衡量绿水与降水、温度、Rn的变化相关性, 计算公式为[31]:

$r = \frac{{\sum \left( {{X_i} - \frac{{\sum {X_i}}}{n}} \right)\left( {{Y_i} - \frac{{\sum {Y_i}}}{n}} \right)}}{{\sqrt {\sum {{\left( {{X_i} - \frac{{\sum {X_i}}}{n}} \right)}^2}} \sqrt {\sum {{\left( {{Y_i} - \frac{{\sum {Y_i}}}{n}} \right)}^2}} }}$ (5)

式中: rXY的相关系数, Xi为第i年的影响因子值, Yi为第i年的绿水量(mm), N为研究年数15。r取值在−1到1之间: r > 0, 表示两要素正相关; r < 0, 表示两要素负相关。r的绝对值越大, 表明二者的相关性越强。

3 结果与分析 3.1 绿水量时间变化特征

本项研究时间尺度较短, 土壤水分变化量微乎其微, 故年均绿水量计算时仅将绿水流即ET纳入考虑, 在多年月均绿水量计算时则同时考虑了多年月均绿水流和绿水蓄量。坝上4县绿水量的年际变化趋势如图 2所示, 在2001—2015年间, 坝上4县绿水量的变化总体呈不显著(P=0.47)的减小趋势, 多年平均值为371.11 mm, 其中2003—2006年、2008年、2010年、2012年、2013年共8个年份高于平均水平, 其余7个年份绿水量低于平均水平, 在2003—2007年间连续下降, 2007—2015年间绿水量波动起伏较大。波动最为突出的年份为2008年和2009年, 相对变化率分别为16.5%和−16.5%。最大值出现在2003年, 绿水量为415.34 mm; 最小值出现在2009年, 绿水量为322.35 mm, 年最大绿水量与最小绿水量差值为92.99 mm。

图 2 2001—2015年张家口市坝上4县绿水量年际变化趋势图 Fig. 2 Annual variation of green water of the four counties of Zhangjiakou in Bashang area during 2001−2015

对多年月均绿水量进行空间统计得到多年月均绿水量分布图(图 3), 绿水量的多年月均值介于2.5~121.8 mm, 月份之间差异较大, 呈现周期性的单峰变化趋势, 峰值出现在7、8月。多年月均绿水量的季节性差异很强, 春季(3—5月)绿水总量介于70.64~121.54 mm, 夏季(6—8月)为141.15~ 324.94 mm; 秋季(9—11月)绿水总量为50.88~ 96.56 mm, 冬季(12—2月)仅为14.42~27.52 mm, 整体表现为夏季 > 春季 > 秋季 > 冬季, 坝上4县绿水量的变化基本与水热同期。

图 3 2001—2015年张家口市坝上4县多年月均绿水量分布图 Fig. 3 Average monthly green water of the four counties of Zhangjiakou in Bashang area during 2001−2015
3.2 绿水量空间变化特征

张家口坝上4县2001—2015年平均绿水量为279.02~540.50 mm, 且具有显著的空间异质性(图 4a), 由西北至东南呈现很明显的递增规律。其中沽源县的绿水量最高, 多年平均值为442.5 mm, 该地降水量较大, 植被覆盖度高且多为林地, 绿水量相应较大; 康保县的绿水量最低, 多年平均值为343.7 mm, 该地降水量低且空气干燥, 绿水量较小。像元尺度绿水量变化趋势介于−10.60~2.59 mm∙a−1(图 4b), 仅有0.2%的区域变化显著(P < 0.05)。总体上呈减小趋势, 其中有21.2%的区域绿水量变化呈增加趋势, 有78.8%的区域呈减少趋势。

图 4 2001—2015年张家口市坝上4县年均绿水量(a)及绿水变化趋势(b)图 Fig. 4 Annual average (a) and linear trend (b) of green water of the four counties of Zhangjiakou in Bashang area during 2001−2015
3.3 绿水量变化驱动因素分析 3.3.1 气象因子时空变化特征分析

降水是蒸散发水分的主要来源, 降水量直接影响ET。2001—2015年坝上4县的平均降水量为348.84~435.38 mm, 年均降水量空间上由西北向东南逐渐增加(图 5a), 与多年平均绿水量的空间分布有明显的递变一致性。其中康保县的年均降雨量最小, 沽源县最大。2001—2015年降水量整体呈不显著的增加趋势(图 5d), 研究区东部降水增加较为缓慢, 西部增加较快, 高值区集中分布于尚义县和张北县西部。气温是反映气候变化的另一关键因子, 也是影响绿水量的另一重要气象要素。坝上4县2001—2015年的年均温为0.85~7.67 ℃ (图 5b), 由东北向西南逐渐减小, 尚义县气温最低, 沽源县气温最高。气温的年际变化幅度较小(图 5e), θslope的绝对值最大为0.11, 整体呈不显著的增加趋势(P=0.62)。其中有68.5%的区域呈增加趋势, 31.5%呈减小趋势。研究区Rn的多年平均值为2712.30~3097.22 MJ∙m−2, 随海拔升高由北向南不断增加(图 5c)。Rn年际波动较大, 如图 5f所示, 变化趋势自西向东不断增加, 在西部呈减少趋势, 在东部呈增加趋势。其中有68.9%的区域呈减小趋势, 31.1%的区域呈增加趋势, 区域整体上呈不显著的下降趋势(P=0.65)。

图 5 2001—2015年张家口市坝上4县各气象因子年均值及其变化趋势分布图 Fig. 5 Annual average and linear trends of various weather influencing factors of the four counties of Zhangjiakou in Bashang area during 2001−2015
3.3.2 GPP、土地利用时空变化分析

水循环、碳循环是陆地表层耦合的两个基本生物物理过程[32]。GPP与绿水分别作为水、碳循环的重要产物, 二者的变化密切相关。如图 6a所示, GPP的多年平均值为140.29 g(C)∙m−2, 2001—2015年间变化幅度较大, 最小值为61.16 g(C)∙m−2, 最大值为255.37 g(C)∙m−2。空间上由西北至东南逐渐增加, 高值区集中在沽源县。变化趋势(图 6b)也由西北向东南逐渐增加, 整体呈不显著的增加趋势(P=0.42)。其中31.3%呈减少趋势, 68.7%呈增加趋势, 增加最快的区域多位于沽源县中部和尚义县南部。尚义县南部年均GPP较低但增长速度较快, 这与该区域林地面积扩大有密切关系。

图 6 2001—2015年张家口市坝上4县年均总第一性生产力(a)及其变化趋势(b)分布图 Fig. 6 Annual average (a) and linear trend (b) of gross primary productivity (GPP) of the four counties of Zhangjiakou in Bashang area during 2001−2015

土地利用变化直观地反映了人类活动对下垫面的作用, 是影响绿水变化的重要原因。因土地利用在短时间内变化较小且土地利用产品为多年合成产品, 故此处基于2000年、2010年、2015年土地利用数据进行分析, 并将原始产品土地利用类型合并为耕地、林地、草地、水体、建设用地和未利用地6类, 对原始数据中错分的水体部分基于野外调查数据进行更正(图 7)。研究区内不同土地利用类型面积为耕地 > 草地 > 林地 > 未利用地 > 建设用地 > 水体(表 2), 2001—2015年耕地、草地、水体面积有所缩减, 林地面积连续增加, 建设用地面积不断扩张, 未利用地面积未发生较大变化。林地的增长面积小于耕地、草地、水体的缩减面积, 耕地、草地被建设用地侵占, 湖泊干涸演变为未利用地, 这在一定程度上使得ET总体上呈减少趋势。

图 7 2001年、2010年、2015年张家口市坝上4县土地利用类型图 Fig. 7 Spatial distribution of land use types of the four counties of Zhangjiakou in Bashang area in 2001, 2010 and 2015
表 2 2001年、2010年、2015年张家口市坝上4县各土地利用类型面积占比 Table 2 Proportions of different land use types of the four countiesof Zhangjiakou in Bashang area in 2001, 2010 and 2015  
3.3.3 绿水量对主要驱动因子的响应关系

水热条件的变化极大地影响着蒸散, 结合绿水量年内变化与水热同期的现象, 本文选取了降水、温度、Rn 3个气象因子探讨绿水量对气候变化的响应, 以降水表征水分条件、以气温和Rn表征热力条件。同时将GPP纳入主要驱动因子考虑范畴, 以GPP反映研究区内绿色植物蒸散能力。

基于IDL分别计算各因子与绿水量的相关系数, 在像元尺度上分别分析绿水和降水、温度、Rn、GPP的相关性。2001—2015年各影响因子与绿水的相关系数分布如图 8所示。坝上4县降水量与绿水量相关系数(图 8a)空间分布不均, 介于−0.04到0.68之间, 显著相关区域仅占13.9% (P < 0.05)。其中呈正相关关系的区域占研究区面积的98%以上, 尤其以康保县西北部及尚义县南部的小部分区域最为明显, 低值区多位于张北县西部和沽源县东南部。研究区内温度与绿水整体上呈负相关关系(图 8b), 相关系数最大值为0.18, 最小值为−0.76, 28.1%的区域呈显著相关(P < 0.05)。有大于99%的区域呈负相关, 仅有极少区域呈正相关。相关性较强的区域集中于沽源县。Rn与绿水相关性较弱, 仅有0.02%的像元呈显著相关(P < 0.05), 相关系数介于−0.62到0.43之间(图 8c), 总体上与绿水量呈现负相关的关系, 只存在15.9%的正相关区域。负相关相关性最强的区域分布于研究区南部边界, 正相关区域集中分布于康保县西部、尚义县北部以及沽源县东北部的小部分区域。如图 8d所示, GPP与绿水的正相关性很强, 平均相关系数为0.63, 其中75.0%的区域呈显著正相关(P < 0.05), 相关系数最大值高达0.96, 只有极小部分区域呈负相关关系。空间分布自西向东递减, 其中尚义县的GPP与绿水的正相关关系最为明显, 平均相关系数达0.74。

图 8 张家口市坝上4县各影响因子与绿水相关系数分布图 Fig. 8 Spatial distribution of correlation coefficients between influencing factors and green water of the four counties of Zhangjiakou in Bashang area

研究区内绿水变化受多种因素共同作用, 故以标准化后的气温、降水、Rn和GPP为自变量, 以绿水量为因变量进行多元回归分析, 定量计算各影响因子对绿水变化的贡献。结果表明绿水量与气温、Rn呈负相关关系, 与降水、GPP呈正相关关系, 与相关分析结果一致。对绿水的影响程度为GPP > 气温 > 降水 > Rn, 其中GPP对绿水量变化的贡献率最大, 为45%, 气温的贡献率为28%, 降水的贡献率为16%, Rn对绿水量变化的贡献率最小(11%)。

3.3.4 绿水量与土地利用的关系

提取2001年、2010年、2015年各土地利用类型实际绿水量, 结果如图 9所示。各年份绿水量均表现为林地 > 草地 > 耕地 > 建设用地。未利用地在研究区内多为干涸湖泊、沼泽等地类, 且分布较为分散, 受来自周围林地、水体等高绿水量区水汽的影响, 绿水量也较高。绿水量高值区的土地利用类型多为草地, 草地的植被覆盖度较高且该区域降水量较大, ET相应较大; 林地集中分布于绿水量的增长区, 这可能由植树造林工程和总初级生产力增加造成。

图 9 2001年、2010年、2015年张家口市坝上4县各土地利用类型绿水量 Fig. 9 Green water of different land use types of the four counties of Zhangjiakou in Bashang area in 2001, 2010 and 2015
4 讨论

坝上地区绿水量整体呈减小趋势, 年内变化呈现较强的季节性, 基本与区域水热同期。吴桂平等[33]在鄱阳湖流域、杨秀芹等[34]在淮河流域对实际蒸散发的同期研究中也得到了相同的季节变化特征。绿水的时间变化主要受气温、降水、净辐射等水热条件的影响。坝上4县降水、温度、GPP整体上呈不显著增加趋势, Rn呈不显著下降趋势, 但变化趋势存在显著的空间差异, 这与其他学者的研究结果一致[23, 35-37]。研究发现绿水与降水、GPP呈正相关关系, 与气温、Rn呈负相关关系。坝上4县地处干旱与半干旱区, 夏季贡献了主要的蒸散量, 夏季平均气温约18 ℃, 降水量仅为200 mm左右, 热量充足而降水不足, 温度升高会使植被受水分胁迫, 降水较气温对该地区绿水变化的限制性更强, 这可能是造成绿水与降水呈正相关而与温度呈负相关的原因。绿水与降水呈正相关关系得到了广泛认同[38-39], 温度与绿水的负相关关系在其他流域也有类似表现。如在干旱区石羊河流域, Xiong等[40]研究发现地表温度与ET呈负相关关系, 水分限制是造成这种温度与绿水量显著负相关关系的主要因素, 这和本研究结论大体一致。Goroshi等[41]在印度干旱半干旱区、张明明[42]在中国干旱半干旱区的坝上地区、马梓策等[38]在锡林河流域、杨宇娜等[43]在吉林泰周边等的研究中也得到绿水与气温呈负相关的结论。此外, 坝上4县绿水量与GPP呈现显著的正相关关系, 气温的变化与GPP的年际波动均相反(图 10), 可能也是气温与绿水量负相关关系形成的重要原因。

图 10 2001—2015年张家口市坝上4县总第一性生产力与温度变化图 Fig. 10 Annual variations of gross primary productivity and temperature of the four counties of Zhangjiakou in Bashang area from 2001 to 2015

气温、降水等因素在时间尺度上影响绿水, 土地利用则更直接地影响绿水的空间分布。2000年以来, 国家大力实施退耕还林还草等生态建设工程, 研究区内植被覆盖变化明显。各土地利用类型绿水量多年平均值差异较小, 但绿水量变化趋势受土地利用变化影响显著。研究区内绿水增加区集中于尚义县南部和沽源县, 2001—2015年间两地林地面积增加显著, 绿水增加区与林地面积增加区相对一致; 绿水量剧烈减少区(θslope≤−5)位于张北县中部, 安固里淖干涸可能是导致绿水量剧烈减少的重要原因。安塞[44]在京津冀地表蒸散发的研究中发现了ET的空间分布与土地利用类型分布非常一致, 与植被覆盖度有良好的相关性; 李发鹏等[45]在黄河三角洲源区的研究中强调了陆面蒸散发的空间分布主要受人类活动的影响。这与本研究得出土地利用影响绿水空间分布的结果相对一致。

本研究基于遥感数据产品进行绿水的时空特征及其影响分析, 遥感数据产品虽在许多学者的研究中得到验证且具有一定区域适用性, 但仍存在一定不确定性。其中, 绿水数据分辨率为0.01°, 土地利用数据为30 m, 二者空间分辨率的差异可能是造成各土地利用类型年均绿水量差异不明显的重要原因。此外, 绿水变化还受其他气候因子(如风速、日照时数等)作用的影响, 且各气候因子之间也存在复杂的相互作用, 多因素水热耦合对绿水时空演变的影响机理还有待进一步研究。

5 结论

基于多源数据产品分析了2001—2015年张家口坝上康保县、沽源县、尚义县、张北县4县绿水量的时空变化特征, 并利用相关系数、线性变化系数、贡献率评估了不同影响因子的时空特征及其对绿水变化的影响。得到以下结论:

1) 坝上4县的绿水量总体呈不显著下降趋势, 多年平均值为371.37 mm。绿水量存在明显的季节分异, 整体表现为夏季 > 春季 > 秋季 > 冬季。

2) 坝上4县多年平均绿水量由西北至东南递增, 沽源县的绿水量最高, 康保县的绿水量最低。空间上, 78.8%的区域绿水量呈减少趋势, 21.2%的区域绿水量呈增加趋势且主要集中在东南部和西南部。

3) 坝上4县绿水量整体上与降水、GPP呈正相关关系, 与气温、Rn总体呈负相关关系。其中, GPP对绿水变化的贡献率最高, 依次为GPP > 气温 > 降水 > Rn

4) 研究区内各土地利用类型绿水量表现为林地 > 草地 > 耕地 > 建设用地。土地利用变化对绿水空间分布影响较大。

参考文献
[1]
FALKENMARK M. Land-water Linkages: A Synopsis. Land and Water Integration and River Basin Management[R]. Rome: Landand Water Bulletin, 1995: 15-16
[2]
刘昌明. 蓝水与绿水的确切定义及其转化机制问题[J]. 南水北调与水利科技, 2020, 18(1): 25.
LIU C M. Precise definition of blue water and green water and their mutual transformation[J]. South-to-North Water Transfers and Water Science & Technology, 2020, 18(1): 25.
[3]
ALEXANDRATOS N, WILEY J, CHICHESTES S. World Agriculture: Towards 2010[R]. Rome: FAO, 1995
[4]
ROCKSTRÖM J, FALKENMARK M, KARLBERG L, et al. Future water availability for global food production: The potential of green water for increasing resilience to global change[J]. Water Resources Research, 2009, 45(7): W00A12.
[5]
许炯心. 黄河中游绿水系数变化及其生态环境意义[J]. 生态学报, 2015, 35(22): 7298-7307.
XU J X. Increasing trend of green water coefficient in the middle Yellow River Basin and the eco-environmental implications[J]. Acta Ecologica Sinica, 2015, 35(22): 7298-7307.
[6]
GERTEN D, HOFF H, BONDEAU A, et al. Contemporary "green" water flows: Simulations with a dynamic global vegetation and water balance model[J]. Physics and Chemistry of the Earth, Parts A/B/C, 2005, 30(6/7): 334-338.
[7]
夏露, 宋孝玉, 符娜, 等. 黄土高塬沟壑区绿水对土地利用和气候变化的响应研究——以南小河沟流域为例[J]. 水利学报, 2017, 48(6): 678-688.
XIA L, SONG X Y, FU N, et al. Impacts of land use change and climate variation on green water in the Loess Plateau Gully Region-A case study of Nanxiaohegou basin[J]. Journal of Hydraulic Engineering, 2017, 48(6): 678-688.
[8]
ABBASPOUR K C, FARAMARZI M, GHASEMI S S, et al. Assessing the impact of climate change on water resources in Iran[J]. Water Resources Research, 2009, 45(10): W10434.
[9]
刘昌明, 李艳忠, 刘小莽, 等. 黄河中游植被变化对水量转化的影响分析[J]. 人民黄河, 2016, 38(10): 7-12.
LIU C M, LI Y Z, LIU X M, et al. Impact of vegetation change on water transformation in the middle Yellow River[J]. Yellow River, 2016, 38(10): 7-12. DOI:10.3969/j.issn.1000-1379.2016.10.002
[10]
夏露. 基于绿水理论的砚瓦川流域生态水文过程对变化环境的响应[D]. 西安: 西安理工大学, 2019
XIA L. Responses of eco-hydrological process to changing environment in Yanwachuan basin based on green water theory[D]. Xi'an: Xi'an University of Technology, 2019
[11]
ZANG C, LIU J, JIANG L, et al. Impacts of human activities and climate variability on green and blue water flows in the Heihe River Basin in Northwest China[J]. Hydrology and Earth System Sciences Discussions, 2013, 10(7): 9477-9504.
[12]
赵安周, 赵玉玲, 刘宪锋, 等. 气候变化和人类活动对渭河流域蓝水绿水影响研究[J]. 地理科学, 2016, 36(4): 571-579.
ZHAO A Z, ZHAO Y L, LIU X F, et al. Impact of human activities and climate variability on green and blue water resources in the Weihe River basin of northwest China[J]. Scientia Geographica Sinica, 2016, 36(4): 571-579.
[13]
冯畅, 毛德华, 周慧, 等. 气候变化对涟水流域蓝水绿水资源的影响[J]. 长江流域资源与环境, 2017, 26(10): 1525-1537.
FENG C, MAO D H, ZHOU H, et al. Impacts of climate change on blue and green water resources in the Lianshui River basin[J]. Resources and Environment in the Yangtze Basin, 2017, 26(10): 1525-1537.
[14]
MOREIRA A A, RUHOFF A L, ROBERTI D R, et al. Assessment of terrestrial water balance using remote sensing data in South America[J]. Journal of Hydrology, 2019, 575: 131-147. DOI:10.1016/j.jhydrol.2019.05.021
[15]
HE M Z, KIMBALL J S, YI Y H, et al. Satellite data-driven modeling of field scale evapotranspiration in croplands using the MOD16 algorithm framework[J]. Remote Sensing of Environment, 2019, 230: 111201. DOI:10.1016/j.rse.2019.05.020
[16]
贺添, 邵全琴. 基于MOD16产品的我国2001-2010年蒸散发时空格局变化分析[J]. 地球信息科学学报, 2014, 16(6): 979-988.
HE T, SHAO Q Q. Spatial-temporal variation of terrestrial evapotranspiration in China from 2001 to 2010 using MOD16 products[J]. Journal of Geo-Information Science, 2014, 16(6): 979-988.
[17]
TANG R L, SHAO K, LI Z L, et al. Multiscale validation of the 8-day MOD16 evapotranspiration product using flux data collected in China[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2015, 8(4): 1478-1486. DOI:10.1109/JSTARS.2015.2420105
[18]
梁顺林, 程洁, 贾坤, 等. 陆表定量遥感反演方法的发展新动态[J]. 遥感学报, 2016, 20(5): 875-898.
LIANG S L, CHENG J, JIA K, et al. Recent progress in land surface quantitative remote sensing[J]. Journal of Remote Sensing, 2016, 20(5): 875-898.
[19]
王彦芳, 裴宏伟. 1980-2015年河北坝上地区生态环境状况评价与对策研究[J]. 生态经济, 2018, 34(1): 186-190.
WANG Y F, PEI H W. Evaluation on eco-environment quality of Bashang region in Hebei Province from 1980 to 2015[J]. Ecological Economy, 2018, 34(1): 186-190.
[20]
李志祥, 田明中, 武法东, 等. 河北坝上地区生态环境评价[J]. 地理与地理信息科学, 2005, 21(2): 91-93.
LI Z X, TIAN M Z, WU F D, et al. Ecological environment evaluation on Bashang District in Hebei Province[J]. Geography and Geo-Information Science, 2005, 21(2): 91-93. DOI:10.3969/j.issn.1672-0504.2005.02.022
[21]
HE J, YANG K, TANG W J, et al. The first high-resolution meteorological forcing dataset for land process studies over China[J]. Scientific Data, 2020, 7(1): 25. DOI:10.1038/s41597-020-0369-y
[22]
ANG K, HE J, TANG W J, et al. On downward shortwave and longwave radiations over high altitude regions: Observation and modeling in the Tibetan Plateau[J]. Agricultural and Forest Meteorology, 2010, 150(1): 38-46. DOI:10.1016/j.agrformet.2009.08.004
[23]
张晓龙. 基于遥感和通量观测的实际蒸散发机理研究与模型应用[D]. 西安: 西安理工大学, 2019
Zhang X L. Research on the mechanism and model of actual evapotranspiration from satellite and flux observations[D]. Xi'an: Xi'an University of Technology, 2019
[24]
YANG K, CHEN Y Y, HE J, et al. Development of a daily soil moisture product for the period of 2002-2011 in Chinese mainland[J]. Science China Earth Sciences, 2020, 63(8): 1113-1125. DOI:10.1007/s11430-019-9588-5
[25]
YANG K, ZHU L, CHEN Y Y, et al. Land surface model calibration through microwave data assimilation for improving soil moisture simulations[J]. Journal of Hydrology, 2016, 533: 266-276. DOI:10.1016/j.jhydrol.2015.12.018
[26]
YANG K, KOIKE T, KAIHOTSU I, et al. Validation of a dual-pass microwave land data assimilation system for estimating surface soil moisture in semiarid regions[J]. Journal of Hydrometeorology, 2009, 10(3): 780-793. DOI:10.1175/2008JHM1065.1
[27]
YANG K, WATANABE T, KOIKE T, et al. Auto-calibration system developed to assimilate AMSR-E data into a land surface model for estimating soil moisture and the surface energy budget[J]. Journal of the Meteorological Society of Japan, 2007, 85A: 229-242. DOI:10.2151/jmsj.85A.229
[28]
王琥. 华北农牧交错带地区人工草地水分利用与水分消耗研究[D]. 北京: 中国农业科学院, 2007: 24
WANG H. Study on the utilization and consumption of water for sward in semiarid agro-pastoral transitional zone of North China[D]. Beijing: Chinese Academy of Agricultural Sciences, 2007: 24
[29]
张守红, 刘苏峡, 莫兴国, 等. 阿克苏河流域气候变化对潜在蒸散量影响分析[J]. 地理学报, 2010, 65(11): 1363-1370.
ZHANG S H, LIU S X, MO X G, et al. Assessing the impact of climate change on reference evapotranspiration in Aksu River Basin[J]. Acta Geographica Sinica, 2010, 65(11): 1363-1370. DOI:10.11821/xb201011006
[30]
王强, 张廷斌, 易桂花, 等. 横断山区2004-2014年植被NPP时空变化及其驱动因子[J]. 生态学报, 2017, 37(9): 3084-3095.
WANG Q, ZHANG T B, YI G H, et al. Tempo-spatial variations and driving factors analysis of net primary productivity in the Hengduan mountain area from 2004 to 2014[J]. Acta Ecologica Sinica, 2017, 37(9): 3084-3095.
[31]
马磊, 闫浩文, 何毅, 等. 2001-2015年喜马拉雅南麓地区植被变化遥感监测[J]. 干旱区地理, 2017, 40(2): 405-414.
MA L, YAN H W, HE Y, et al. Vegetation changes in south Himalayas areas based on remote sensing monitoring during 2001-2015[J]. Arid Land Geography, 2017, 40(2): 405-414.
[32]
于贵瑞, 王秋凤, 于振良. 陆地生态系统水-碳耦合循环与过程管理研究[J]. 地球科学进展, 2004, 19(5): 831-839.
YU G R, WANG Q F, YU Z L. Study on the coupling cycle of water-carbon and process management in terrestrial ecosystem[J]. Advance in Earth Sciences, 2004, 19(5): 831-839. DOI:10.3321/j.issn:1001-8166.2004.05.022
[33]
吴桂平, 刘元波, 赵晓松, 等. 基于MOD16产品的鄱阳湖流域地表蒸散量时空分布特征[J]. 地理研究, 2013, 32(4): 617-627.
WU G P, LIU Y B, ZHAO X S, et al. Spatio-temporal variations of evapotranspiration in Poyang Lake Basin using MOD16 products[J]. Geographical Research, 2013, 32(4): 617-627.
[34]
杨秀芹, 王磊, 王凯. 基于MOD16产品的淮河流域实际蒸散发时空分布[J]. 冰川冻土, 2015, 37(5): 1343-1352.
YANG X Q, WANG L, WANG K. Spatio-temporal distribution of terrestrial evapotranspiration in Huaihe River basin based on MOD16 ET data[J]. Journal of Glaciology and Geocryology, 2015, 37(5): 1343-1352.
[35]
孙雷刚, 王绍强, 徐全洪, 等. 京津冀地区生态系统生产力估算及时空格局分析[J]. 地理与地理信息科学, 2019, 35(5): 109-117.
SUN L G, WANG S Q, XU Q H, et al. Ecosystem productivity and its spatio-temporal pattern in Beijing-Tianjin-Hebei region[J]. Geography and Geo-Information Science, 2019, 35(5): 109-117.
[36]
杨佳, 郝桂珍, 张婧, 等. 1960-2016年冀西北地区降水时空变化特征研究[J]. 水资源与水工程学报, 2019, 30(5): 117-123.
YANG J, HAO G Z, ZHANG J, et al. The spatial and temporal variation characteristics of precipitation in Northwest Hebei Province during 1960-2016[J]. Journal of Water Resources and Water Engineering, 2019, 30(5): 117-123.
[37]
李元菲, 张兰霞, 曹永强, 等. 河北省潜在蒸散量时空变化特征及气候影响因素分析[J]. 南水北调与水利科技, 2019, 17(3): 67-78.
LI Y F, ZHANG L X, CAO Y Q, et al. Spatiotemporal variations of potential evapotranspiration and its climate influencing factors in Hebei Province[J]. South-to-North Water Transfers and Water Science & Technology, 2019, 17(3): 67-78.
[38]
马梓策, 于红博, 张巧凤. 2000-2017年锡林河流域地表蒸散量的时空特征及其影响因素[J]. 中国农村水利水电, 2020(3): 18-24.
MA Z C, YU H B, ZHANG Q F. Spatio-temporal characteristics of evapotranspiration and its influencing factors in Xilin river basin from 2000 to 2017[J]. China Rural Water and Hydropower, 2020(3): 18-24.
[39]
黄葵, 卢毅敏, 魏征, 等. 土地利用和气候变化对海河流域蒸散发时空变化的影响[J]. 地球信息科学学报, 2019, 21(12): 1888-1902.
HUANG K, LU Y M, WEI Z, et al. Effects of land use and climate change on spatiotemporal changes of evapotranspiration in Haihe River basin[J]. Journal of Geo-Information Science, 2019, 21(12): 1888-1902. DOI:10.12082/dqxxkx.2019.190269
[40]
XIONG Y J, ZHAO S H, YIN J, et al. Effects of evapotranspiration on regional land surface temperature in an arid oasis based on thermal remote sensing[J]. IEEE Geoscience and Remote Sensing Letters, 2016, 13(12): 1885-1889. DOI:10.1109/LGRS.2016.2616409
[41]
GOROSHI S, PRADHAN R, SINGH R P, et al. Trend analysis of evapotranspiration over India: Observed from long-term satellite measurements[J]. Journal of Earth System Science, 2017, 126(8): 1-21. DOI:10.1007/s12040-017-0891-2
[42]
张明明. 2000-2015年中国干旱半干旱区蒸散发时空变化及其影响因素分析[D]. 西安: 长安大学, 2019: 28
ZHANG M M. Analysis of the temporal and spatial variation of evapotranspiration and its driving factors in arid and semi-arid region of China from 2000 to 2015[D]. Xi'an: Chang'an University, 2019: 28
[43]
杨宇娜, 汪季, 张成福, 等. 吉兰泰及周边地区蒸散发的时空变化规律[J]. 灌溉排水学报, 2019, 38(S2): 30-36.
YANG Y N, WANG J, ZHANG C F, et al. Spatial and temporal variations of evapotranspiration in Jilantai and its surrounding areas[J]. Journal of Irrigation and Drainage, 2019, 38(S2): 30-36.
[44]
安塞. 京津冀地区土地利用/覆被与地表蒸散发动态变化及关系研究[D]. 石家庄: 河北科技大学, 2018: 47-56
AN S. Research on dynamic change and relationship between land use/cover and evapotranspiration in Beijing-Tianjin-Hebei region[D]. Shijiazhuang: Hebei University of Science and Technology, 2018: 47-56
[45]
李发鹏, 徐宗学, 李景玉. 基于MODIS数据的黄河三角洲区域蒸散发量时空分布特征[J]. 农业工程学报, 2009, 25(2): 113-120.
LI F P, XU Z X, LI J Y. Characteristics of the spatial and temporal distribution for regional evapotranspiration in the Yellow River Delta based on MODIS data[J]. Transactions of the Chinese Society of Agricultural Engineering, 2009, 25(2): 113-120.