气候、径流及河漫滩湿地淹没状态耦合关系研究——以海拉尔河下游段为例

董 熙,陈志泊,胡春明,刘 平

1 北京林业大学信息学院, 北京 100083

2 国家林业和草原局智能信息处理工程技术研究中心, 北京 100083

3 中国科学院生态环境研究中心城市与区域生态国家重点实验室, 北京 100085

河漫滩湿地是水文过程的空间集成介质,径流量影响着地表河漫滩湿地资源,径流不仅是陆地水循环与水量平衡的重要组成部分,也最有价值和最易开发利用的自然资源之一[1]。同时,河水淹没频率反应了水体淹没时间的长短和次数,被认为是影响湿地植被生态系统的最重要的一种水文因素[2]。因此综合分析径流、水体淹没频率以及其影响因素,对改善区域水资源和河流管理具有非常重要的意义[3]。

国内外学者对径流变化特征、演变规律及其驱动因素开展了大量的研究工作。丁艳峰等通过研究发现,气候变化以及日益频繁的人类活动对黄河流域水资源影响严重,造成黄河径流量减少,产生断流[4]。穆兴民等研究表明,黄河天然径流量年际振荡较大,变差系数在0.19—0.20之间波动,其径流量年际变化具有阶段性特征[5]。胡弟弟等在北极开展了勒拿河流域气候变化及其对径流的影响的研究,发现径流受温度和降雨影响较大,与气候因素呈正相关[6]。吴锦奎等在河西走廊石羊河、黑河流域的研究表明,河流径流的变化特征与多年冻土活动层增厚导致的流域退水减慢有关[7];王晓颖等基于气候变化和人类活动对白河流域径流变化影响的定量研究,研究发现人类活动对白河径流有抑制作用[8]。上述前人学者在气候对径流的驱动以及径流的变化规律等方面进行一系列的研究,但是以往对径流的研究主要集中在气候对径流的影响以及径流变化的驱动因素等方面,而对气候与径流以及河流面积变化的相互关系鲜有研究。

近年来,由于全球河漫滩湿地不断退化,国内外学者也开展了河漫滩湿地的相关研究。由于径流量的变化进而引起河流淹没区域以及湿地面积的变化,这一系列变化对生态系统和环境健康产生了巨大威胁[9—12]。例如,胡春明以二卡湿地为为研究区,开展了生态阈值研究,研究发现中、高径流期湿地面积比例可近似作为保护区湿地面积的适宜值和理想值[13]；苟芳珍在苏干湖湿地开展了河流径流补给对漫滩湿地植被生物量空间格局的影响研究[14]；李理在淇河流域进行了径流供给、调节对生态系统服务功能的研究[15]；在国外,如在恒河,印度河,阿木河谷和锡尔河等一些大型流域,发生了周期性的断流情况[10]。由于流域周期性的断流使周围的湿地生态系统遭到破坏,造成生物多样性减少[9—12]。近年来遥感(RS)技术已成为监测水体信息的重要工具[8—10]。遥感技术的发展使得湖泊、河流监测技术从传统的水文站发展到更加全面的遥感监测[16—19]。Dong等人基于遥感的方法监测湿地情况,以Water Performance Frequency(WPF)为基础,提出淹没保证率模型Water Inundate Guarantee Rate (WIGR),来分析周期性淹没对于河漫滩湿地健康的重要性[20—21]。Sarda同样以遥感监测为技术手段,结合Water Performance Frequency(WPF)和深度数据分析了湿地的水文稳定性[22]。Feng等人指出长江下游由于水利枢纽建设而导致的流量减少了17%—27%,影响流域淹没状态,并提出建议如果不尽快采取适当的管理策略,将导致水资源短缺[23]。在上述学者河漫滩湿地淹没状态研究工作中,往往集中于遥感尺度,没有考虑河流径流对湿地生态系统健康的影响,并且结合水文数据的研究较少。

针对上述研究存在的问题,本研究将遥感数据与河流水文数据相结合,综合考虑气候与径流对淹没状态的影响,通过建立气候-径流-淹没状态的关系,填补在寒区河流漫滩湿地淹没状态对气候响应的研究空白。本研究选取位于海拉尔河下游段作为研究对象,基于淹没频率的模型,通过对气候变化的特征、趋势、规律的分析,从径流变化的角度了解气候变化在寒区河流径流变化中的作用,分析气候对径流产生的影响以及对河漫滩湿地淹没状态的影响,分析探讨三者的关系,以期为河流管护决策者提供理论参考和方法支持。

1 研究区与数据

1.1 研究区概况

研究区位置见图1,海拉尔河位于中国内蒙古自治区呼伦贝尔市境内,发源于牙克石市乌尔其汉镇境内大兴安岭吉鲁契那山西麓,是额尔古纳河的上游河段,呈东西流向。干流全长622km。该地区属于大陆性草原气候,位于中高纬度温带半干旱区。冬季严寒漫长,春季干旱多大风,夏季温凉短促,秋季降温急剧、霜冻早,年光照充足,春季降水少,空气干燥,是全年相对湿度最小的季节,平均在47%左右；夏季降水多,空气中水汽含量丰富,相对湿度达65%—73%；秋季气温下降,相对湿度略高于春季,在58%—66%之间。

图1 海拉尔河位置

海拉尔河为补给型河流,径流年内分配不均,主要集中在汛期。根据资料统计,多年平均年径流量为36.95亿m3,6—9月径流量占年径流量的56.7%—71.4%；尤其是主汛期7—8月的径流量,占年径流量的32.3%—41.9%；枯水期5月来水较少,仅占年径流量的12.1%—15.6%。海拉尔河洪水分春季融冰洪水和夏季暴雨洪水。融冰洪水一般出现在4—5月,洪水特点为峰高量小,历时较短,约一个月左右；夏季洪水一般出现在7—8月,洪水特点为峰平量大,历时长,一般历时达60—70d,洪峰流量出现的次数每年1—2次。

1.2 数据来源与预处理

1.2.1数据来源

从美国国家航空航天局(NASA)获取本研究所用Landsat 5 TM/8 OLI TIRS 1T级图像,影像覆盖整个研究区,行列号为123026,影像获取时间1989—2018年年间(5- 10月),所有可用影像云量低于2%,共计45期。

海拉尔河1989—2018年逐日流量数据来自水文勘测局,数据记录完整；研究区1989—2018年逐月的降雨数据,通过中国气象数据网(http://data.cma.cn)下载得到,水文数据用于后续计算水体淹没保证率模型,气象数据用于分析计算海拉尔地区气候变化趋势。坝后水文站位于海拉尔河中段,所以该水文站所测径流信息能代表研究区水文状况。

1.2.2数据预处理

本研究的数据预处理包括影像的辐射定标、大气校正、几何校正、正射校正、影像裁剪等。本研究使用ENVI 5.3中的辐射度校准工具(Radiomotric Calibration)计算辐射度,将辐射文件输入到FLAASH大气校准模块中。有研究表明[23]使用多光谱数据计算光谱指数时,可以使用简单的快速大气校正或暗像元法,为减少由大气散射和吸收引起的辐射误差,故本研究ENVI 5.3中FLAASH大气校正法。所有反射率产品的值都在0到1之间,并对获得的所有反射产品进行几何校正,并且误差均小于1,完成对所有遥感数据的预处理后,通过波段计算工具计算水体指数,并对数据进行裁剪。

2 研究方法

2.1 水体面积的确定

各种遥感指数可用于水体提取以及从卫星图像进行湿地测绘。诸如NDWI[19],MNDWI[24], AWEIsh[20]等。在本研究中,利用MNDWI[24]水体指数绘制研究区水体地图,并提取水体面积。已有研究表明,MNDWI指数方法提取水体精度更高,同时,错、漏分误差也较低,且在狭长状水体的提取上效果较好[25]。为了验证湿地图,本研究于2019、2020、2021年3次实地调研海拉尔干流及支流,共采集水体样点100个, 并结合相对较高分辨率的Google地球历史影像(50个点)验证水体湿地提取精度,并且计算Kappa系数。

(1)

式中,Green,NIR,MIR是指Landsat影像的绿、近红外和中红外波段,MNDWI值大于0表示水体。

本研究对研究区域的水体提取准确性进行了定量评估,基于误差矩阵的整体准确性和Kappa系数被用来评估水体提取的准确性,Kappa系数越接近1则表明准确性越高[25],Kappa系数可计算为：

(2)

式中,P是参考数据的像素总数；Pii类别中正确像素的总数；Pi+是从分类数据中得出的类别的像素总数；P+i是从参考数据中得出的类别的像素总数；n是类别总数。

2.2 水体淹没频率

利用水体淹没频率Water Performance Frequency(WPF)来[21—22]表示海拉尔河下游段及两岸湿地的淹没状态,用以反映径流对淹没状态的影响,其数学表达式为：

(3)

式中,(xi,yi)表示影像中某一像素,T表示该像素在研究的所有遥感影像中被淹没的期数,Z表示研究的所有遥感影像总期数。

用Matlab编程计算逐个像素的WPF,基于WPF值将各像素的淹没情况分为：极低(0—20%)、低(21%—40%)、中(41%—60%)、高(61%—80%)、极高(81%—100%)五类。

2.3 径流变化对气候的响应

以坝后水文站逐日径流资料为基础,研究海拉尔河下游段年径流总量年际变化特征,并计算径流量与气温、降雨的距平、累计距平直观反映各要素的变化趋势,并分析气候与径流以及淹没面积的关系。

3 研究结果

3.1 水体面积提取与精度验证

研究使用MNDWI指数提取水体面积,并整理了对应日期的流量数据。在42期遥感影像对应的42个日均流量数据中(图2),超200m3/s大流量数据共计14个,其中：5月产生的大流量数据6个,占42.85%；6—10月大流量数据8个,占57.15%。水体面积提取结果显示(表1)：研究区水体面积平均值45.31km2；最大值出现在1995年6月28日为107.15km2,对应日均流量381m3/s；最小值出现在1998年5月3日,对应日均流量49.9m3/s。

图2 海拉尔河坝后站45个日均流量数据曲线

Kappa系数验证结果表明(表1)：MNDWI水体提取精度大多在85%以上,MNDVI所描绘的湿地地图可以基本反映真实地面情况,数据可用。

表1 提取水体精度、kappa系数与流域面积

3.2 水体淹没频率

图3代表研究区的淹没频率图,WPF高于80%区域主要为海拉尔河主河道,WPF低于20%区域则主要分布在海拉尔河两岸湿地。

不同淹没频率区域的面积统计见表2：淹没频率区域总面积为194.48km2；极低WPF(0—20%)区域面积132.38km2,占比达到68.07%；而低(21%—40%)、中(41%—60%)、高(61%—80%)、极高(81%—100%)WPF区域占比分别为12.71%、6.01%、4.45%和8.76%。

表2 不同水体淹没频率区域的面积统计

3.3 海拉尔气候变化

3.3.1年降雨量变化

1989—2018年期间(图4),海拉尔气象站多年平均降水量356.45mm,其中2017年降水量最小为215.4mm,2013年降水量最大为619.1mm,最大降雨量与最小降雨量相差403.7mm,表明海拉尔地区年际间降雨量差异较大。

海拉尔地区年降水量整体表现为下降趋势,年均减少11.64mm。1989年开始到2018年底,降水距平正负不断交替出现,1998年以前累计距平持续大于0,表明该时期降水相对较多(图4),1998年后累计距平持续小于0,表明降雨量低于均值。

图4 海拉尔年际降雨变化、距平、累计距平

3.3.2年际气温变化

如图5所示,海拉尔地区气温年际变化幅度较大,海拉尔气象站多年平均气温为-0.2℃,其中2012年平均温度最低为- 2.15℃,2007年平均温度最高为1.26℃,最高均温与最低均温相差3.4℃,年均气温整体呈下降趋势,年均变化0.01℃。

由图5可知,2009年之前距平值几乎均为正,2009年后2015年距平温度为正,研究期内正负距平值交替出现,累计距平值出现三次上升三次下降的趋势。

图5 海拉尔年际气温变化、距平、累计距平

3.4 径流年际变化

多年均径流量32.96×108m3,年径流量最大值出现在 2013年为87.48×108m3,年径流量最小值出现在 2017年为10.23×108m3,最大值与最小值相差77.25×108m3,年径流量减少,年均减少2.05×108m3。

坝后站实测年径流量变幅较大,研究时间段内呈“s”形变化(图6),在1989—2018研究时间段内,坝后站所监测的区域内年径流呈下降趋势。由累计距平曲线先上升、后下降的变化趋势可以看出年径流量不断减少(图6)。

图6 坝后站年径流量变化趋势、距平、累计距平

4 讨论

4.1 气候变化与人类活动对径流的影响

表3显示了坝后站所测年径流与年均气温和年降水量的相关与偏相关系数。由表3可知,年降雨量、年均气温与坝后站所测径流量的相关关系。年均降雨量、年均气温与径流呈现显著相关与负相关关系；在年降雨量作为控制变量的条件下,气温与径流的偏相关系数为-0.095,说明二者存在线性影响但影响有限,在年均气温作为控制变量的条件向,降雨与径流偏相关系数为0.805,明二者存在线性影响且影响显著。

表3 年度径流与年降雨、气温相关性与偏相关分析

见图7,坝后站所监测的径流从1989—2018年整体下降,气象站所测降雨和气温下降。1989—1993年降雨数据下降,1989—1993年的径流明显下降；1995—1998年降雨有明显的上升过程,导致这个时间段内的年度径流也有一个明显的增加过程；2011—2013年降雨有明显的上升过程,径流也有明显上升的过程,随后降雨下降,径流下降。

图7 海拉尔河年度径流与气温、降雨相关变化趋势

见图7,温度对径流的影响是综合且复杂的[26—27]。如2001—2005年间,温度变化与径流变化表现出一致性,即：气温升高,径流增加,造成这一现象的主要原因是：海拉尔河属于东北冻土区,在这一时间段内气温升高,使得冰雪消融和土壤冻结推后,特别是进入 21世纪以来,暖冬的出现导致地面升温明显,冻土退化释放地下水补给河流,同时冻土活动层厚度增加,季节性冻结减少,土壤冻融的变化加强了水分向地下的迁移,使得河流径流增多,对气候变暖的响应显著,在这一时间段内,虽然气温升高,对河流的蒸散有一定的作用,但是在这期间降雨量缓慢增加,直接补给了河流径流,所以该时间段内,河流径流增加明显[28]。另一方面温度与径流在一定程度上表现出了相反的变化趋势,如2010—2013年间,气温下降,径流增加,造成这一现象的可能原因是：随着年度温度的下降导致冻土的增加,冻土冷储构成冷控系统,导致土壤含水量增大,增加下渗增加基流,一定程度上增加了河流径流量,同时,温度降低河流蒸散减少,上述原因总体对河流产流有增强作用[28—30],与陆胤昊等人对海拉尔河径流的结果相同[31]。

刘腊山等人的研究表明[26—31],在高寒地区,气温对径流的影响在高寒主要表现在2个方面：1)影响积雪的消融；2)影响河流的蒸发量,积雪消融产生的融水是高寒地区河流径流的主要补给来源之一。近年来,研究区气温变化明显,呈现“S”增减趋势(图7),随着气温的变化,海拉尔河研究区冰雪融水量发生变化,进一步引起河流径流量的变化,进而导致淹没区域的变化。

由上文图2、表1大流量数据出现时间可知,超200m3/s流量数据共14个, 5月产生的大流量数据6个,由于研究区5月几乎没有降雨(雪),主要依靠冰雪融化补给径流,所以在该时段内,温度变化起到主要影响径流的作用,6—10月大流量数据8个,该时间段内研究区降雨充沛,降雨直接作用于径流,所以6—10月海拉尔河径流的变化由降雨占主导影响[28]。

本研究以1989、2018年部分海拉尔河地区统计年鉴为反应人类活动因素的指标,见表4,首先当地区生产总值增加,且增速较快,相比前一时间段增长40倍有余,人口、工业生产总值、粮食产量、耕地面积均有明显增长,不同反应人生活动因素的指标的增长,导致对海拉尔河水资源的过度利用与开发,直接导致径流减少,在研究时限内,海拉尔地区工农业得到空前发展,基础设施也不断完善,人类活动通过改变下垫面条件进而影响流域水循环过程,其对流域径流变化的影响因素众多,如土地利用变化、水资源开发、水利工程等均有明显的影响。水资源的不合理利用,造成流域水资源开发利强度加剧,在一定程度上改变了径流状态。综上,径流的下降是自然因素及人类活动影响综合作用的结果,由于研究区地广人稀,自然因素对径流的影响占主导因素,同时,将自然因素与人类活动因素对径流的影响的具体量化将是下一步的研究重点。

表4 海拉尔1989年、2018年统计年鉴

4.2 径流对淹没状态的影响

如图8所示,本文将45期数据的提取面积与对应日期的流量(表1)的数据进行了拟合,研究发现,海拉尔河的当日遥感影像所提取的水体面积和当日的流量数据呈高度拟合,拟合系数R2分别达到0.8525,又对该数据进行了皮尔逊相关性分析,海拉尔河的当日遥感影像所提取的河流面积和当日的流量数据相关系数是0.923**(表5)。

图8 海拉尔河水体面积与对应日期流量的拟合关系

表5 研究时限内不同时间的流域面积与月均径流量

当流量数据增大时,对应日期遥感影像所提取的湿地面积越大,流量与淹没面积关系均表现出同样趋势,随着流量的增大,研究区及两侧区域的淹没面积也相应增加,其原因在于：在没有人为因素干扰的情况下,研究区河段为平原河段,流量较大时,洪水的漫滩效应较明显,使得坝后站流量与所提取的淹没面积表现出较好的一致性,由此可见二者呈现显性正相关,并有良好的线性关系。

维持河流中的生态流量至关重要,河流流量的稳定有利于该河流的生态环境稳定[21—23]。各淹没频率下的淹没区域,受影响程度不同,大于80%极高淹没频率区域集中在主河道附近,河流径流影响最小,0—20%极低频率区域分布河岸两侧的河漫滩湿地当中,该区域最易受径流变化影响,由于径流与淹没面积呈现正相关,则径流量的损失会导致横向水体的交换减少,减少了淹没面积的大幅度波动,导致对河漫滩湿地水文状况的产生负面影响,进而影响淹没频率0—20%极低淹没频率区域；相反的是,当径流量增加,会增加河水出槽次数,极低频率淹没区域淹没次数增加,增大了淹没面积的波动,同时,保证各频率下淹没区域面积稳定,对河漫滩湿地生态系统的发育起到促进作用。

4.3 气候-径流-淹没状态的相关关系

水文情势与气候之间存在显着相关性,水文情势的变化通常认为是反映气候变化的常用指标,因为气候变化可以通过改变水文过程事件的幅度、持续时间、频率、时间、可预测性等,进而改变河流系统的水文周期[32—34]。

河道径流量的影响因素包括气温、降水等因素,地质条件、下垫面状况等自然因素以及人口、资源开发、项目建设等社会经济因素[35—36],但是,在自然状态下的河流,降雨以及融雪起到对河流的径流量进行调节的主要作用,在海拉尔河流域地区,气温升高或者降低会导致的河流的蒸发加大以及冻土退化,这也是影响径流年内分配的重要因素之一。1989年以来,该地区气温呈缓慢下降趋势,但下降趋势缓慢,由于温度的降低,导致海拉尔河地区冰封时间增长,会抑制春季径流的产生,导致春夏秋季的径流的延迟,影响河流对河漫滩的湿地水份供应[37—40]。湿地面积随时间的推移而发生的变化,由于受研究区降雨、融雪以及周边地区的城市发展,导致河流径流的减少,河流的水份供给不足会减少漫滩与河流间的横向连接度,进而减少河水漫出河槽次数,减少河漫滩湿地淹没频率[41—42],从而对河岸的河岸两侧极低淹没频率区域的生态系统服务功能产生抑制作用。

因此,可以得出结论,气候变化对径流产生影响,由于气温下降,降雨的减少,致使河道天然径流过程,趋于稳定、均匀的水文过程改变水沙情势,进而影响海拉尔河的淹没状态,河水出槽次数减少,不同河流淹没频率都会受到影响,导致河流及两岸廊道生物栖息环境的多样性下降,对水生生物、湿地生态系统的健康产生影响。

5 结论

本文以Landsat TM/OIL遥感数据的光谱指数为数据基础,根据中国气象数据中心及海拉尔河坝后水文站提供的 1989—2018 年间的气温、降雨和径流数据,分析了海拉尔河下游淹没状态对气候变化的响应。结果表明：

(1)海拉尔河地区的气候经历了变干变冷的过程,温度、降雨量下降,研究期内海拉尔站监测降雨量多年平均356.45mm,年均减少11.64mm,多年平均温度-0.2℃,年均变化0.01℃；

(2)海拉尔河年径流呈减少趋势,坝后监测站所测段海拉尔河径流量年均减少2.05×108m3,在自然状态下导致海拉尔河径流减少的主要原因是:一方面海拉尔河地区属于东北冻土区,随着温度下降导致冻土的增加,春季融雪补给径流减少,从温度层面上对海拉尔河产生影响,另一方面夏秋季降雨直接补给径流,随着降雨减少,径流减少；

(3)径流与淹没面积呈现正相关,径流量的损失会导致横向水体的交换减少,减少淹没面积的大幅度波动,导致对河漫滩湿地水文状况的产生负面影响,各淹没频率下的淹没区域,受影响程度不同,大于80%极高淹没频率区域受河流径流影响最小,0—20%极低频率区域受径流变化影响最大。

(4)气候—径流—淹没状态存在相关关系。海拉尔河温度下降,导致海拉尔河地区冰封时间增长,抑制春季径流的产生,降雨减少,导致夏秋季的径流延迟,影响河流对河漫滩的湿地水份供应,河流的水份供给不足会减少漫滩与河流间的横向连接度,进而减少河水漫出河槽次数,减少河漫滩湿地淹没频率,影响河漫滩湿地淹没状态,从而对河岸的河岸两侧20%以下极低淹没频率区域的生态系统服务功能产生抑制作用。