为了探讨气候变化下沙河流域地表径流演变趋势以及其对径流的影响,本文基于1960~2020年沙河流域逐日气象水文数据,采用统计分析法对比分析了流域降水、温度和径流的演变趋势性和突变性,并结合沙河流域1980、1990、2000、2010、2020年五期土地利用数据,建立了多时段沙河流域SWAT (Soil and Water Assessment Tool)模型。根据未来趋势设置不同气候变化情景,分析了沙河流域径流对气候变化的响应。研究结果表明:1) SWAT模型在沙河流域径流模拟中具有较好的适用性,率定期和验证期的相关系数(R2)均在0.6以上,Nash-Sutcliffe纳什效率系数(Ens)均在0.5以上,相对误差(Re)均在±25%以内。2) 过去61年沙河流域温度呈显著上升趋势,降水量下降不明显,而实测径流量呈显著下降趋势。3) 不同气候变化情景下的模拟表明,研究区径流与降水呈正相关,与温度呈负相关,且径流对降水变化的敏感程度高于径流对温度变化的敏感程度。 To explore the evolution trends of surface runoff in the Shahe River Basin under climate change and its impacts, we used daily meteorological and hydrological data from 1960 to 2020 in the Shahe River Basin. We employed a statistical analysis method to compare and analyze the trend and abrupt changes in precipitation, temperature, and runoff in the basin. Additionally, we established a multi-period SWAT (Soil and Water Assessment Tool) model for the Shahe River Basin, considering land use data for five periods in 1980, 1990, 2000, 2010, and 2020. We set different climate change scenarios according to the future trends to analyze the response of runoff in the Shahe River Basin to climate change. The research results indicate: 1) The SWAT model demonstrates good applicability in simulating runoff in the Shahe River Basin, with correlation coefficients (R²) above 0.6 during the calibration and validation periods, Nash-Sutcliffe efficiency coefficients (Ens) above 0.5, and relative errors (Re) within ±25%. 2) Over the past 61 years, the temperature in the Shahe River Basin showed a significant upward trend, with a less pronounced decrease in precipitation, while observed runoff exhibited a significant decreasing trend. 3) Simulations under different climate change scenarios revealed that runoff in the study area is positively correlated with precipitation and negatively correlated with temperature. Runoff is more sensitive to changes in precipitation than to changes in temperature.
为了探讨气候变化下沙河流域地表径流演变趋势以及其对径流的影响,本文基于1960~2020年沙河流域逐日气象水文数据,采用统计分析法对比分析了流域降水、温度和径流的演变趋势性和突变性,并结合沙河流域1980、1990、2000、2010、2020年五期土地利用数据,建立了多时段沙河流域SWAT (Soil and Water Assessment Tool)模型。根据未来趋势设置不同气候变化情景,分析了沙河流域径流对气候变化的响应。研究结果表明:1) SWAT模型在沙河流域径流模拟中具有较好的适用性,率定期和验证期的相关系数(R2)均在0.6以上,Nash-Sutcliffe纳什效率系数(Ens)均在0.5以上,相对误差(Re)均在±25%以内。2) 过去61年沙河流域温度呈显著上升趋势,降水量下降不明显,而实测径流量呈显著下降趋势。3) 不同气候变化情景下的模拟表明,研究区径流与降水呈正相关,与温度呈负相关,且径流对降水变化的敏感程度高于径流对温度变化的敏感程度。
气候变化,SWAT模型,沙河流域,径流
Xian Zhao1, Jianwei Sun2, Zhenghu Lyu3, Guanghui Liu1, Shuai Ren4, Ruiting Liu1
1School of Water Conservancy and Hydroelectric Power, Hebei University of Engineering, Handan Hebei
2Handan Water Conservation Promotion Center, Handan Water Conservancy Bureau, Handan Hebei
3Xingtai Hydrological Survey Research Center of Hebei Province, Xingtai Hebei
4School of Earth Science and Engineering, Hebei University of Engineering, Handan Hebei
Received: Oct. 23rd, 2023; accepted: Jan. 8th, 2024; published: Jan. 18th, 2024
To explore the evolution trends of surface runoff in the Shahe River Basin under climate change and its impacts, we used daily meteorological and hydrological data from 1960 to 2020 in the Shahe River Basin. We employed a statistical analysis method to compare and analyze the trend and abrupt changes in precipitation, temperature, and runoff in the basin. Additionally, we established a multi-period SWAT (Soil and Water Assessment Tool) model for the Shahe River Basin, considering land use data for five periods in 1980, 1990, 2000, 2010, and 2020. We set different climate change scenarios according to the future trends to analyze the response of runoff in the Shahe River Basin to climate change. The research results indicate: 1) The SWAT model demonstrates good applicability in simulating runoff in the Shahe River Basin, with correlation coefficients (R²) above 0.6 during the calibration and validation periods, Nash-Sutcliffe efficiency coefficients (Ens) above 0.5, and relative errors (Re) within ±25%. 2) Over the past 61 years, the temperature in the Shahe River Basin showed a significant upward trend, with a less pronounced decrease in precipitation, while observed runoff exhibited a significant decreasing trend. 3) Simulations under different climate change scenarios revealed that runoff in the study area is positively correlated with precipitation and negatively correlated with temperature. Runoff is more sensitive to changes in precipitation than to changes in temperature.
Keywords:Climate Change, SWAT Model, Shahe River Basin, Runoff
Copyright © 2024 by author(s) and beplay安卓登录
This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).
http://creativecommons.org/licenses/by/4.0/
随着社会的快速发展,温室气体大量排放,全球平均温度明显上升。2021年,IPCC在发布的第六次气候变化评估报告中表明,全球的平均温度在未来20年内将再升高大约1.5℃左右 [
近年来,众多学者针对气候变化对国内外不同流域径流的响应开始进行大量研究 [
沙河为海河流域子牙河水系滏阳河支流,主要是邢台市行洪的主要河道,担负起城市补水以及生态恢复功能的作用。但目前分析沙河流域气候变化对径流影响的研究成果十分稀少,仅是在海河流域相关研究中作为一个子流域进行了分析,研究成果远不能满足该流域水资源管理和开发利用需求。因此,本文以沙河流域为研究区,通过统计方法对沙河流域的降水、温度和径流演变趋势和特征进行分析,建立了不同时期沙河流域SWAT水文模型,在此基础上采用任意情景假设法模拟不同气候变化情景对径流的影响,这将为该地区水资源的合理规划和配置提供重要的科学依据,从而对流域的生态环境保护具有极其重要的意义。
沙河发源于河北省内丘县白鹿角川的山丘区地带,沙河又名大沙河,全长166 km,流域面积2280 km2,其中在邢台境内全长53 km,流域面积1804 km2,是邢台市境内流经距离最长、流域面积最大、流经县市最多、流量最大的河流。本文以邢台境内的沙河为研究对象,研究流域面积为1681.37 km2。该流域属于暖温带半湿润大陆性季风型气候,季节交替明显,年平均温度12.8℃,年平均降水量571.4 mm,降水量年际变化大,而且年内分配极不均匀,汛期(6~9月)降水量约占全年的75%以上。沙河流域地理位置如图1所示。
图1. 沙河流域地理位置
构建SWAT模型所需的基础数据如表1所示。主要包括数字高程DEM图、五期土地利用数据、土壤数据以及气象水文数据,其气象数据来源于6个雨量站(日降水量)和2个气象站(日平均风速、最高、最低温度、相对湿度和日照时数);水文数据来源于朱庄水库水文站(逐月实测径流)。此外,运行SWAT模型之前,需要构建DEM数据库、土地利用数据库、土壤数据库以及气象数据库,并且对空间数据进行相同的投影坐标处理,选取统一UTM中WGS_1984_UTM_Zone_49N的投影坐标。
数据类型 | 分辨率 | 数据格式 | 数据年份 | 数据来源 |
---|---|---|---|---|
数字高程图 | 30 m | GRID | 2020 | 地理空间数据云 |
土地利用图 | 30 m | GRID | 1980、1990、2000、2010和2020 | 中国科学院资源环境中心 |
土壤类型图 | 1 km | GRID | - | 世界土壤数据库(HWSD) |
气象数据 | 日 | TXT | 1960~2020 | 中国气象站 |
水文数据 | 月 | TXT | 1960~2020 | 水文年鉴 |
表1. SWAT模型基础数据来源
SWAT模型是由美国农业部农业研究中心(USDA-ARS)开发的分布式水文模型,其能够模拟复杂流域内不同的土地利用类型、土壤类型和土地管理措施下的各水文要素的长期变化 [
利用数字高程DEM数据库、土地利用数据库、土壤数据库以及气象数据库,分别对五期土地利用情况进行构建沙河流域SWAT模型,划分为30个子流域。根据径流突变以及五期土地利用数据,将1960~2020年气象水文数据划分为1960~1978年、1979~1990年、1991~2000年、2001~2010年、2011~2020年五个时段,为了减小误差,都设置一年预热期。利用SWAT-CUP软件对朱庄水库水文站1960~2020年逐月实测径流数据对模型的参数进行率定和验证 [
对模型参数进行率定和验证之前,首先对模型的参数进行参数敏感性分析,然后选用敏感性最高的参数。参考其他文献 [
E n s = 1 − ∑ i = 1 n ( Q 0 i − Q m i ) 2 ∑ i = 1 n ( Q 0 i − Q 0 ¯ ) 2 (2-1)
R 2 = [ ∑ i = 1 n ( Q 0 i − Q 0 ¯ ) ( Q m i − Q m ¯ ) ] 2 ∑ i = 1 n ( Q 0 i − Q 0 ¯ ) 2 ∑ i = 1 n ( Q m i − Q m ¯ ) 2 (2-2)
R e = Q m − Q 0 Q m × 100 % (2-3)
在月尺度模拟径流情况下,当率定期与验证期的数据要同时满足纳什效率系数(Ens)大于0.5、相关系数(R2)大于0.6和相对误差(Re)小于25%时,说明满足模型精度要求,表明模型在研究流域具有适用性。
本文基于沙河流域长时间序列1960~2020年降水、温度和径流数据,采用滑动平均法 [
将基于未来可能的气候变化范围进行研究,这些范围包括温度增加或减少1℃、2℃、3℃,以及降水增加或减少10%、20%等情况。将通过数量化方法来分析单一气候要素对径流的影响,或者对不同气候要素进行不同的组合,以研究在这些不同组合情景下径流的影响过程。
根据IPCC第六次评估报告的结果 [
情景SC1:采用1980年土地利用数据,保持沙河流域1960~2020年降水不变,将各气象站逐日温度均升高0℃、1℃、2℃、3℃和降低1℃、2℃、3℃的7种方案。
情景SC2:采用1980年土地利用数据,保持沙河流域1960~2020年温度不变,将各雨量站逐日降水均增加0%、10%、20%和减少10%、20%℃的5种方案。
情景SC3:采用1980年土地利用数据,将各气象站逐日温度均升高1℃、2℃、3℃和降低1℃、2℃、3℃,各雨量站逐日降水均增加10%、20%和减少10%、20%℃,温度和降水两两综合组合情景下的23种方案。
对朱庄水库水文站所在子流域出口的不同时段月实测径流量与模拟径流量进行对比来看(图2),SWAT模型模拟沙河流域的径流与实测径流拟合较好,率定期和验证期的相关系数(R2)均在0.6以上,纳什效率系数(Ens)均在0.5以上,相对误差(Re)均在±25%以内,均满足三个评价指标,表明了SWAT模型在沙河流域具有良好的适用性。
对沙河流域6个雨量站:朱庄水库站、渡口站、蝉房站、路罗站、侯家庄站以及槲树滩站1960~2020年年降水量数据进行线性趋势分析,可以看出年降水量极大值出现在1963年,多年降水量最小值在朱庄
图2. 五个时段率定期和验证期月径流实测值和模拟值
水库站的1987年的29.6 mm。之后根据线性回归检验和M-K非参数检验法进行显著性程度可知,六个雨量站线性倾向率整体呈现不显著减少,说明降水在这61年里表现出不显著减少趋势,除了槲树滩站倾斜率呈现正向,其他几个站都呈现负向,倾向率不同,表现的降低速率有所不同,沙河流域降水量趋势变化特征见图3(a)~图3(f)。采用M-K突变检验和累积距平法对降水数据进行突变检测(图4、图5),发现朱庄水库站、渡口站、蝉房站雨量站的降水时间序列突变年为1978年,路罗站和侯家庄站突变年为1974年,槲树滩站突变年为1968年。
图3. 沙河流域不同站降水、温度年际趋势变化特征图
沙河、和顺气象站1960~2020平均温度趋势变化特征见图3(g)、图3(h),不同气象站线性倾向率不同,呈现出来的增加速率也有所不同,沙河站温度以0.1℃/10a的速率呈显著增加,和顺站温度以0.2℃/10a的速率呈显著增加,表明沙河流域1960~2020年总体气温呈现显著上升趋势,这与全球气候变化趋势一致 [
图4. 沙河流域不同站点降水、温度序列M-K突变点检验
图5. 沙河流域不同站点降水序列累积距平
沙河流域在1960~2020年的年平均径流量为17.32亿m3,最大年径流量为190.86亿m3,出现的年份为1963年,最小年径流量为2011年,径流量为2.99亿m3。朱庄水库水文站径流量趋势变化见图6(a),根据M-K非参数检验和线性回归检验(表2)可知,沙河流域以−5.0亿m3/10a的速率呈显著减少,说明近61年沙河流域径流量呈现显著下降趋势。首先对朱庄水库水文站实测年径流序列进行M-K突变检验,结果如图6(b),UF和UB两条曲线在0.05显著性临界值直线之间有一个交点,说明径流量突变点发生在1978年,为了进一步检验证明,采用了累积距平法来判断突变点(图6(c)),由图可知,在1978年径流量累积距平值达到最高值之后,累积距平值的曲线开始呈现下降趋势,结合M-K突变检验和累积距平值的结果,表明1978年是朱庄水库水文站径流量时间序列发生的突变年份。当流域发生显著人类活动后,水文系列的自然过程会受到破坏,从而导致趋势性变化或跳跃性突变,故选择不受人类活动影响天然状态下1960~1978年为基准期,干扰期为1979~2020年。
图6. 沙河流域径流量变化趋势、M-K突变检验、累积距平值
气象水文 | 站点 | 倾向率(mm/10a或 ℃/10a或亿m³/10a) | 变化趋势 | 线性回归 | Mann-Kendall | 显著性 | ||
---|---|---|---|---|---|---|---|---|
检验值|T| | T(0.05/2) | 统计量|Zc| | Zc(0.05) | |||||
降水 | 朱庄水库站 | −21.6 | 下降 | 1.23 | 1.64 | 0.752 | 1.96 | 不显著 |
渡口站 | −3.7 | 下降 | 0.24 | 1.64 | 0.285 | 1.96 | 不显著 | |
蝉房站 | −25.1 | 下降 | 1.33 | 1.64 | 0.888 | 1.96 | 不显著 | |
路罗站 | −12.1 | 下降 | 0.72 | 1.64 | 0.148 | 1.96 | 不显著 | |
侯家庄站 | −21.2 | 下降 | 1.33 | 1.64 | 0.752 | 1.96 | 不显著 | |
槲树滩站 | 4.6 | 上升 | 0.25 | 1.64 | 1.458 | 1.96 | 不显著 | |
温度 | 沙河站 | 0.1 | 上升 | 3.08 | 1.64 | 3.007 | 1.96 | 显著 |
和顺站 | 0.2 | 上升 | 8.35 | 1.64 | 6.389 | 1.96 | 显著 | |
径流 | 朱庄水库水文站 | −5.0 | 下降 | 2.58 | 1.64 | 3.496 | 1.96 | 显著 |
表2. 标准试验系统结果数据
为了进一步研究流域径流对气候变化的响应,设定了35种的降水和温度两个主要因素不同的气候变化情景,气候变化中的降水和温度的数据来源于雨量站和气象站来进行情景设计,采用率定好的SWAT水文模型,并假设保持土地利用不变的条件下,分析不同气候变化情景对径流的影响。本文采用1980年土地利用类型图,根据1960~1978年已率定好的SWAT模型,模拟1960~2020年天然径流,并与不同气候变化情景组合依次输入已率定SWAT水文模型模拟径流进行对比分析,结果见表3。
温度变化(℃) | 降水变化(%) | ||||
---|---|---|---|---|---|
−20% | −10% | 0% | 10% | 20% | |
−3 | −28.64 | −5.99 | 26.73 | 58.71 | 89.90 |
−2 | −36.64 | −9.83 | 17.28 | 48.77 | 79.64 |
−1 | −43.86 | −18.34 | 7.97 | 38.77 | 69.29 |
0 | −49.56 | −25.42 | - | 30.14 | 60.22 |
1 | −53.68 | −30.90 | −6.33 | 23.03 | 52.61 |
2 | −56.92 | −35.56 | −12.14 | 16.21 | 44.96 |
3 | −59.86 | −39.80 | −17.53 | 9.70 | 37.68 |
表3. 沙河流域不同气候变化情景下径流变化率
1) 温度对径流的影响
由图7的情景SC1模拟结果可知,从趋势方面进行分析,当降水保持不变时,随着温度从1℃、2℃、3℃增加,模拟径流逐渐呈现减少趋势,模拟径流的减少率为−6.33%、−12.14%、−17.53;当温度从−1℃、−2℃、−3℃逐渐降低,模拟径流逐渐呈增加趋势,模拟径流的增加率为7.97%、17.28%、26.73%。影响结果可以看出,径流随着温度减少而增加,表明温度与径流呈负相关关系,与第三章数理统计分析温度趋势一致。当温度减少和温度增加的幅度相同时,发现温度减少对径流的影响程度较大,而且随着温度的降低,径流的增加比例也随之增大,可见温度下降对径流影响的敏感程度高于温度上升的敏感程度。总体可以看出,当降水不变时,温度对径流整体的影响程度不大。
图7. 温度变化情景下径流相对变化率
2) 降水对径流的影响
由图8的情景SC2模拟结果可知,从趋势方面进行分析,当温度保持不变时,随着降水从−10%、−20%减少,模拟径流逐渐呈现减少趋势,模拟径流的减少率为−25.42%、−49.56%;当降水从10%、20%逐渐增加,模拟径流逐渐呈增加趋势,模拟径流的增加率为30.14%、60.22%。影响结果可以看出,径流随着降水增加而增加,表明降水与径流呈正相关关系,与第三章数理统计分析降水趋势一致。当降水减少和降水增加的幅度相等时,观察到降水增加对径流的影响程度较大,而且随着降水量的增加,径流的增加比例也相应增大,观察表明,对于径流影响,降水增加的敏感程度高于降水减少。总体而言,径流对降水变化的响应相对于温度变化的响应更加敏感,当温度变化在±3℃范围内时,径流的变化率比较小,仅在26.73%以内。然而,当降水发生±20%的变化时,径流的变化率可以达到60.22%以内。这表明,降水对于径流的影响更大,并且径流对降水的变化更加敏感。
图8. 降水变化情景下径流相对变化率
3) 温度和降水变化组合对径流的影响
由图9的情景SC3模拟结果可知,当降水和温度同时变化情景下,当降水减少20%时,温度从1℃、2℃、3℃逐渐增加,模拟径流的减少率为−53.68%、−56.92%、−59.86%。从趋势方面进行分析,模拟径流呈减少趋势,当降水减少20%时,温度从−1℃、−2℃、−3℃逐渐降低,模拟径流的减少率为−43.86%、−36.64%、−28.64%,温度降低和降水减少对径流起相互抑制作用,模拟径流逐渐呈现减少趋势,进一步说明降水相对温度对径流的影响程度大。当降水减少20%时,温度增加比温度减少对径流的影响程度大,可见温度增加对径流影响的敏感程度高于温度减少的敏感程度。当降水减少10%时,与降水减少20%趋势基本一致,可见温度增加对径流影响的敏感程度高于温度减少的敏感程度。当降水增加20%时,温度从1℃、2℃、3℃逐渐增加,模拟径流的增加率为52.61%、44.96%、37.68%,可以看出,温度增加和降水增加对径流起相互抑制作用,模拟径流逐渐呈现增加趋势,当降水增加20%时,温度从−1℃、−2℃、−3℃逐渐减少,模拟径流的增加率为69.29%、79.64%、89.90%,模拟径流呈增加趋势。进一步说明降水相对温度对径流的影响程度大。当降水增加20%时,径流随着温度减少而增加,温度减少比温度增加对径流的影响程度大,可见温度减少对径流影响的敏感程度高于温度增加的敏感程度。降水增加10%与降水增加20%趋势基本一致,可见温度降低对径流影响的敏感程度高于温度增加的敏感程度。表明降水是导致径流变化的直接因素,而温度是导致径流变化的次要因素。
图9. 温度和降水组合变化情景下径流相对变化率
文章基于1960~2020年沙河流域的气象水文数据,采用数理统计方法分析了降水、温度和径流的变化规律,并建立了多时段沙河流域SWAT模型,分析了不同气候情景设置下径流变化,得到以下结论:
1) 沙河流域1960~2020年温度呈现显著上升趋势,径流呈显著减少趋势,降水整体呈不显著下降趋势,仅有槲树滩站降水量呈不显著上升趋势。
2) 构建了不同土地利用情景下沙河流域SWAT模型,对不同时段10个较敏感参数进行了率定和验证,相关系数(R2)、纳什效率系数(Ens)以及相对误差(Re)均满足评价指标要求,表明所构建的SWAT模型适用于沙河流域。
3) 模拟了三种气候变化情景对径流的影响,表明沙河流域径流与降水呈正相关,径流对降水变化的敏感程度高于径流对温度变化的敏感程度。
本文考虑到不同时期下垫面的水文物理特性有所不同,前期的下垫面参数无法充分反映后期的水文物理特性。因此,为了更好地了解径流的变化,对径流突变前和突变后进行多时段调参,建立了基于多时段参数的SWAT模型,多时段参数的水文模型能够显著提升模拟精度,具有良好的效果,与前人采用两三期土地利用数据相比较,本文采用每十年五期土地利用数据构建模型,使模拟结果更加精确,能更好地模拟沙河流域径流过程。
赵 娴,孙建伟,吕正虎,刘光辉,任 帅,刘芮廷. 基于SWAT模型的沙河流域气候变化对径流的影响The Impact of Climate Change on Runoff in the Shahe River Basin Based on the SWAT Model[J]. 地球科学前沿, 2024, 14(01): 12-29. https://doi.org/10.12677/AG.2024.141002
https://doi.org/10.1111/j.1752-1688.1998.tb05961.x
https://doi.org/10.1016/S0959-3780(99)00017-5
https://doi.org/10.1016/j.jhydrol.2006.09.014
https://doi.org/10.1016/j.jhydrol.2011.04.038
https://doi.org/10.2166/nh.2010.009
https://doi.org/10.3390/hydrology1010020
https://doi.org/10.1002/hyp.7876
https://doi.org/10.1080/07900629849231
https://doi.org/10.1002/clen.201100652
https://doi.org/10.1016/j.quaint.2014.03.048
https://doi.org/10.1002/hyp.10307