文章信息
- 董李勤, 章光新, 张昆
- DONG Liqin, ZHANG Guangxin, ZHANG Kun
- 嫩江流域湿地生态需水量分析与预估
- Analysis and prediction of wetland ecological water requirements in the Nenjiang Basin
- 生态学报, 2015, 35(18): 6165-6172
- Acta Ecologica Sinica, 2015, 35(18): 6165-6172
- http://dx.doi.org/10.5846/stxb201401100077
-
文章历史
- 收稿日期:2014-01-10
- 网络出版日期:2014-11-19
2. 西南林业大学, 昆明 650224
2. Southwest Forest University, Kunming 650224, China
湿地是发育于水陆环境过渡地带的生态系统,其特定的水文条件是湿地生态系统形成及其功能维持的重要因子[1, 2]。气候变化通过改变全球水文循环现状引起水资源在时空上的重新分布,并对气温、降水与干旱洪涝等极端水文事件发生频率和强度造成直接影响,进而改变湿地蒸散发、径流、水深、水文周期等水文条件和水质状况,最终影响到整个生态系统结构和功能,致使湿地生态系统演替规律及发展趋势发生变化[1, 2, 3]。
嫩江流域属于全球气温变暖幅度最大的北方高纬地区,是全国七大水系之一松花江的北源,也是我国湿地主要集中分布区[4]。分布在流域内的自然湿地,不仅能调节径流、卸载河川径流,同时具有调蓄、净化、涵养水的功能,在维系流域水安全和生态安全,调节区域小气候、水资源补给与生物多样性保护等方面发挥着极其重要的作用,也是阻止松嫩平原西部荒漠化向中东部扩展的重要生态屏障。然而,近百年气候变化带来的气温升高、水资源时空上的重新分布等问题,以及人类对水土资源不合理的开发利用,引起下游湿地入流量减少和蒸散量增加,改变湿地水文特征,造成湿地水资源短缺、旱化现象,进而导致湿地面积萎缩和功能下降[5]。
湿地作为一个有机体,对气候变化和人类活动具有一定的自我调节能力。因此,维持湿地健康所需的水分不是在一个特定的点上,而是在一定范围内变化的,变化的范围就构成了生态系统水分需求的阈值区间。从目前研究情况来看,对湿地生态需水量的研究主要是在年内和年际两个时间尺度上进行估算,年内分为生长期与非生长期,年际分为平水年和枯水年[6]。而对气候变化影响下的湿地生态需水量的变化研究极少[7],更缺乏对嫩江流域整体湿地生态需水量的计算与预测分析。因此,研究气候变化背景下嫩江流域湿地生态需水量的变化状况以及未来气候变化情景模式下湿地生态需水量的变化趋势对湿地保护与湿地水资源管理意义重大,也是未来湿地生态水文的重要研究领域。
1 研究区和研究方法 1.1 研究区概况嫩江流域位于我国东北地区中西部,地跨东北三省,自北向南流经黑龙江省的黑河市、大兴安岭地区、嫩江县、齐齐哈尔市、大庆市,内蒙古的呼伦贝尔盟和吉林省的白城市等县市区。流域全长1370 km,流域面积29.7万 km2,属于温度大陆性的季风气候。近55年嫩江流域多年平均温度约为2.1 ℃,多年平均降水量为454.9 mm,其中夏季降水量占全年降水量的70%—80%,因此夏季易发生洪涝灾害。
嫩江流域北部、西部和南部3面地势较高,为隔水边界,东南部地势低平,形成广阔的松嫩平原,为排泄区。自西向东,嫩江流域大致可分为大兴安岭山区、山区丘陵过渡地带、山前倾斜平原区与中部低平原区。其中大兴安岭山区海拔为1000—1400 m,自北向南分布在流域西侧,西部山前倾斜平原为大兴安岭东麓的山前地带,主要由扇形台地构成,地面高程为143—153 m,水系密布并在低平原过渡地带形成沼泽湿地。低平原地区地势低平开阔,洼地、湖泡星罗棋布,盐碱地、沼泽湿地发育,是各类湿地集中分布区(图 1)。
![]() |
图 1 嫩江流域湿地分布图 Fig.1 Distribution of the wetlands in Nenjiang Basin |
本文的研究应用嫩江流域内国家气象台站水文气象监测资料(包括最高、最低气温、风速、水汽压、相对湿地、日照时数、净辐射、降水等)计算参照作物的潜在蒸散发量;此外,由于嫩江流域内沼泽湿地植被分布以芦苇为主,因此选取芦苇植被为代表并结合FAO(Food and Agriculture Organization)提供的单系数法 Kc (The single crop coefficient)确定沼泽湿地的实际蒸散发量,各时期作物参数指标如表 1所示。
沼泽湿地 Marsh wetland | 时间段/月-日 Period of time/month-date | 芦苇 Kc Reed Kc | 水面 Kc Surface of the water Kc |
Kcini(The initial crop coefficient) | 04-01—05-31 | 1.00 | 0.65 |
Kcmid(The crop coefficient in the mid) | 06-01—08-31 | 1.20 | 0.65 |
Kcend(The crop coefficient in the end) | 09-01—09-30 | 1.00 | 0.65 |
由于FAO推荐的作物系数是针对美国东南部地区,其气候温暖潮湿,而嫩江流域气候干燥,因此根据前人实验结果[8, 9]调整沼泽系数,Kcini、Kcmid和Kcend分别调整为1.48、1.12和1.28,水面 Kc 调整为0.94。利用ArcGIS空间分析软件,将嫩江流域土壤类型数据与湿地遥感数据叠加处理,得出沼泽湿地范围内各土壤类型的分布面积以及土壤厚度等参数。发现嫩江流域内分布着暗棕壤、黑钙土、栗钙土、草甸土、沼泽土、潮土、水稻土等约47种亚类土壤类型,其土壤厚度从34 cm到150 cm不等,平均土壤厚度为102.7 cm,平均土壤容重为1.14 g/cm3,平均田间持水量为40.4%。
湿地生态需水量预测采用的气象数据来自CMIP5(the Fifth phase of the Couple Model Intercomparison Project)全球气候模式的模拟结果,是由包括BCC-CSM1-1、BNU-ESM、CCSM4、CNRM-CM5在内的共21个全球气候模式的模拟结果,经过差值计算统一降尺度到统一分辨率下,并利用简单平均方法进行多模式集合制作而成。CMIP5全球气候模式包括1901—2005年Historical和2006—2100年RCP2.6(Representative Concentration Pathways)、RCP4.5和RCP8.5排放情景下的月平均资料,本文将选取3个排放情景下2006—2100年的模拟结果进行气候变化趋势的分析,并选取2030年、2050年和2100年的模拟结果进行湿地生态需水量的预测分析。
1.2.2 方法广义上讲,湿地生态需水量是指湿地为维持自身发展过程和保护生物多样性所需要的水量,包括植物需水量、土壤需水量和野生生物栖息地需水量3部分;狭义而言,湿地生态需水量是指湿地每年用于生态消耗而需要补充的水量,即补充湿地生态系统蒸散需要的水量。文中湿地生态需水量的计算用的是广义上的概念,其中湿地植物、土壤和野生生物栖息地需水量的内涵如下:
(1) 湿地植物需水量
植被蒸散发量大小与植物种类、组成结构、植被面积及盖度等因素有关,在正常生育状况下(即水分充足),常常采用彭曼公式计算植物潜在蒸散发,再根据植被类型用 Kc 单系数法进行校正。
估算大尺度区域或流域湿地植物需水量时,常常将湿地植被面积与蒸散发量的乘积看作植物需水量,公式为:
式中, dWp 为植物需水量; A(t) 为湿地植被面积; ETm 为蒸散发量;t 为时间。其中,计算蒸散发用的是FAO提供的彭曼公式。
(2) 湿地土壤需水量
湿地土壤需水量与湿地植物生长及其需水量密切相关。不同类型的湿地土壤,其田间持水量、土壤容重、土壤厚度等特性均不相同,需水量也会出现差异。根据研究的需要,通常按照田间持水量或用饱和持水量参数进行计算,从理论公式上可表达为:
式中, Q 为土壤需水量; α 为田间持水量或者饱和持水量百分比;γ 为土壤容重; Ht 为湿地土壤厚度; At 为湿地土壤面积。
(3) 野生生物栖息地需水量
在计算大区域或者流域内野生生物栖息地需水量时,由于湿地分布广而又分散,各类指标不可能一一测出,因此,通常根据栖息地水面面积百分比和水深进行计算。嫩江流域内沼泽湿地面积广而分布比较分散,因而计算野生生物栖息地需水量时,以各重点湿地保护区内的野生生物栖息地需水量代表整个嫩江流域。
依据以上分析,确定通过水面面积百分比和水深要素计算野生生物栖息地需水量的公式为:
式中, Wp 为野生生物栖息地需水量; A(t) 为湿地面积;C 为水面面积百分比; H(t) 为水深;t为时间。
目前为止,湿地生态需水量用到的预测方法有传统预测方法、专家评价方法和情景分析方法等。情景分析预测法是建立在对研究对象未来状态或趋势进行多种可能性推断的一种预测方法,即假设事物的发展变化存在很大不确定性,难以精确预测,因此给出未来发展的多种可能性,情景分析得出的是定性、定量相结合的结论,可以有效地降低预测风险。本文在充分考虑嫩江流域湿地气象因子历史变化趋势的前提下,选取情景分析方法,进行湿地生态需水量的预测分析。
2 嫩江流域湿地生态需水量的计算嫩江流域内沼泽植被分布面积较广,因此首先对不同降水频率下流域内各月沼泽湿地植被需水量进行计算,并根据上述公式对嫩江流域内各种土壤类型需水量以及野生生物栖息地需水量进行计算,得到不同降水频率下嫩江流域沼泽湿地适宜生态需水量,在丰水年(25%降水频率)、平水年(50%降水频率)和枯水年(75%降水频率)湿地适宜生态需水量分别为70.284亿m3、118.696亿m3、169.343亿m3。由于嫩江流域存在多尺度丰枯水变化周期,并且降水与径流两者演变周期基本一致[10],在枯水年,降水量的减少,导致湿地生态需水量的增加。并且从不同降水频率下嫩江流域沼泽湿地适宜生态需水量年内变化过程曲线(图 2)可以看出:嫩江流域沼泽湿地生态需水量年内季节变化比较大,特别是在4—9月份的生长季,枯水年湿地生态需水量呈现较大增幅。具体表现为:5、6月份湿地生态需水量最大,占全年湿地生态需水量的30%以上,主要原因在于春季干燥,降水量少而蒸发量大;7月份湿地生态需水量急剧下降,是由于嫩江流域进入丰水季节,降水量增多,补充湿地需水。从各季节湿地生态需水量变化趋势来看,春季和夏季的湿地生态需水量最大,秋季次之。
![]() |
图 2 不同降水频率下嫩江流域湿地适宜生态需水量年内过程 Fig.2 Annual distribution of optimum ecological water requirements in the Nenjiang Basin under different precipitation frequency |
在枯水年、平水年和丰水年,月平均湿地适宜生态需水量的最大值和最小值如表 2所示:其中最大值出现在生长季,占全年的70%—80%,而最小值出现在冬季,仅占全年的3%左右。从以上分析可以看出,不同降水频率下嫩江流域沼泽湿地生态需水量的变化过程符合季节气候水文条件和湿地生态规律。
不同降水频率/时期 different precipitation frequency/period | 生长季(—9月份) The growing season (April—September) | 冬季(12—2月份) Winter (December—February) |
75%降水频率 75% precipitation frequency | 136.470 | 5.717 |
50%降水频率 50% precipitation frequency | 91.875 | 4.246 |
25%降水频率 25% precipitation frequency | 53.235 | 1.910 |
为了预估未来全球和区域气候变化情况,必须事先提供未来温室气候和气溶胶等的排放情况,即排放情景。排放情景通常是根据一系列因子假设得到的,包括经济发展、人口增长、技术进步、全球化、环境条件、公平原则等。对应未来可能出现的不同社会经济发展状况,要制作不同的排放情景。
此前IPCC曾先后发展了两套温室气体和气溶胶排放情景,及1992的IS92排放情景和2000年的SRES排放情景,但目前来看,早已需要更新和补充[11]。2011年Climatic Change出版专刊,介绍新一代排放情景,并对4种情景分别作了详细的分析[12, 13]。新一代情景叫做“典型浓度目标(Representative Concentration Pathways)”,简称RCPs。4种情景分别为RCP8.5、RCP6.0、RCP4.5和RCP2.6。其中前3个情景模式大致同2000年方案中的SRES A2、A1B和B1相对应。本文选取RCP8.5情景、RCP4.5情景和RCP2.6情景分析未来2030年、2050年和2100年气候变化的情况,以及未来各种情景下流域湿地生态需水量的变化趋势。RCP的3种情景的简单情况如表 3所示。
情景 Scenarios | 具体描述 Detailed description |
RCP8.5 | 辐射强迫上升到8.5 W/m2,2100年之后CO2的当量浓度达到1370×10-6 |
RCP4.5 | 辐射强迫稳定在4.5 W/m2,2100年之后CO2的当量浓度达到650×10-6 |
RCP2.6 | 辐射强迫在2100年之前就达到峰值,到2100年下降至2.6 W/m2,CO2的当量浓度峰值约为490×10-6 |
对2006—2100年嫩江流域范围内RCP2.6、RCP4.5和RCP8.5排放情景下的降水、均温、最高温、最低温变化趋势(图 3)进行分析发现:
![]() |
图 3 嫩江流域1906—2100年年降水量、年均温、年均最低温和年均最高温变化趋势 Fig.3 Trend of annual precipitation,annual temperature,annual minimum temperature and annual maximum temperature in Nenjiang Basin from 1906—2100 |
RCP2.6、RCP4.5和RCP8.5等3种排放情景下的年降水量在2006—2010年差异很大,从图 3中可以看出,RCP2.6情景下的降水量趋势稳定,无明显变化;RCP4.5情景下的降水量呈现先增加、后趋于稳定的变化趋势;RCP8.5情景下的降水量呈持续显著增加趋势。对RCP2.6、RCP4.5和RCP8.5 3种排放情景下嫩江流域年均气温、最高温和最低温进行分析(图 3)发现,2006—2030年3种情景下年均温、最低温和最高温的变化趋势基本一致。2030—2100年,RCP2.6情景下的年均温度增长缓慢并趋于平稳,RCP4.5情景下的年均温度稳定增长,RCP8.5情景下的年均温度增长幅度最大(图 3)。
对比图 3,3种排放情景以年均最低温的增温幅度最大,为0.665 ℃/10a,其次为年均温度,为0.629 ℃/10a,年均最高温度增温幅度最小,为0.611 ℃/10a。年均最高温和最低温之间差值的缩小,导致昼夜温差变小[14],会使一天中的热量循环周期发生变化,使热量累积达到峰值的时间和大小都发生改变,导致某些作物的最佳生长周期的推迟或提前。
3.3 3种不同排放情景下湿地生态需水量变化趋势通过对湿地生态需水量的计算可以看出,湿地生态需水量中以植物需水量所占比重最大,并且受气候变化影响也较大。因此,对3种排放情景下湿地生态需水量的计算只计算植被需水量部分,并根据相应的最高温、最低温和降水变化预测未来3种排放情景下嫩江流域湿地生态需水量的变化趋势(图 4)。
![]() |
图 4 嫩江流域不同排放情景下湿地生态需水量预测 Fig.4 Ecological water requirements prediction under different emission scenarios in the Nenjiang Basin |
对比3种排放情景下嫩江流域的湿地生态需水量可以发现,RCP2.6情景、RCP4.5情景和RCP8.5情景3种排放情景下湿地生态需水量在2030—2100年期间的变化趋势并不一致。到2030年,嫩江流域湿地生态需水量为RCP2.6情景<RCP4.5情景<RCP8.5情景,具体3种排放情景下湿地生态需水量年变化趋势分析如下:
(1) RCP2.6情景下湿地生态需水量呈现先增加后减少的趋势,主要是由于该情景下提倡应用生物质能、恢复森林,是最低端排放情景,其中21世纪后半叶能源应用为负排放,因此到从2050年到2100年,其湿地生态需水量呈现明显减少趋势;
(2) RCP4.5情景下湿地生态需水量在2030—2050年变化不大,主要是由于在此期间,降水量增加幅度明显高于温度,对湿地生态需水有一定的补充作用;2050—2100年期间,湿地生态需水量大幅增加,主要是由于期间降水量呈现减少并趋于稳定变化趋势,但温度却持续增加,使该情景下需水量达到最大值;
(3) RCP8.5情景下湿地生态需水量在2030—2050年期间呈现减少趋势,是由于在此期间降水量增加幅度远远大于温度的增加幅度,增加的降水量不仅可以抵消由于温度升高带来的蒸散发的增加,还可以对湿地有一定的水资源补给作用;但2050—2100年期间,湿地生态需水量有大幅增加,主要是由于在最高温室气体排放情景下,即使降水量有增加趋势,但温度的持续快速增长,以及人口增加、能源改善缓慢带来的影响,都必然会导致湿地生态需水量的增加。
选择3种情景下2030年月湿地生态需水量分析其变化趋势(图 5),可以看出,植物生长季(4—9月份)湿地生态需水量占全年湿地生态需水量的比例达到80%左右,说明植被需水量是湿地生态需水量最主要的组成部分;一年之内以5月份需水量最大,是由于5月份降水不多,但植被开始生长,需水量较大,符合植物生长规律;湿地需水量在5月份达到顶峰后,在7月份迅速降低,是由于7月份降水量最大,补充湿地水资源需求;冬季需水量最小。
![]() |
图 5 不同情景下2030年嫩江流域湿地适宜生态需水量年内过程预测 Fig.5 Predication of annual distribution of optimum ecological water requirements in the Nenjiang Basin under different emission scenarios in 2030 |
(1) 关于嫩江流域内湿地生态需水量以往的研究都是关于扎龙湿地、向海湿地等具体湿地保护区进行需水量计算,缺乏对整个嫩江流域整体沼泽湿地生态需水量的认识和了解。本研究湿地生态需水量的计算是为蓄水量和耗水量之和,不包含降水量,并同时考虑降水频率的影响。研究表明嫩江流域不同降水频率下的湿地生态需水量分别为:丰水年(25%降水频率)70.284亿m3,平水年(50%降水频率)118.696亿m3,枯水年(75%降水频率)169.343亿m3;在年内变化的趋势表现为春夏较大、秋季次之,冬季最小。研究结果也反映了湿地生态需水量与气候条件的相关性。
(2) 选择CMIP全球气候模式,在RCP2.6、RCP4.5和RCP8.5 3种排放情景下,预测2030年、2050年和2100年嫩江流域湿地生态需水量的变化趋势,分析结果表明:RCP2.6情景下需水量呈先增加后减少的趋势;RCP4.5和RCP8.5情景下需水量整体呈增加趋势,到2100年分别达到147.337亿m3和132.659亿m3,其变化受到最高气温、最低气温和降水量变化的共同影响。
在气候变暖、降水量及格局发生变化的现状下,嫩江流域水资源供需矛盾将会更加尖锐。而如何协调水资源需求间的矛盾,维持湿地生态系统健康稳定,实现水资源持续发展,将是未来研究关注的重点。另外,湿地生态需水量是一个动态的值,随着生态保护目标的改变或者经济社会的发展其值也必然不同。
[1] | Brock T C M, Vierssan W V. Climatic change and hydrophte-dominated communities in inland wetland ecosystem. Wetland Ecology and Management, l992, 2(1):37-49. |
[2] | Nielsen D L, Brock M A. Modified water regime and salinity as a consequence of climate change:prospects for wetlands of Southern Australia. Climatic Change, 2009, 95(3/4):523-533. |
[3] | Finlayson C M, Davis J A, Gell P A, Kingsford R T, Parton K A. The status of wetlands and the predicted effects of global climate change:the situation in Australia. Aquatic Sciences, 2013, 75(1):73-93. |
[4] | 刘兴土.东北湿地.北京:科学出版社, 2005. |
[5] | 董李勤, 章光新.嫩江流域沼泽湿地景观变化及其水文驱动因素分析.水科学进展, 2013, 24(2):177-183. |
[6] | 李九一, 李丽娟, 姜德娟, 李俊伟.沼泽湿地生态储水量及生态需水量计算方法探讨.地理学报, 2006, 61(3):289-296. |
[7] | Wang H, Xu S G, Sun L S. Effects of climatic change on evapotranspiration in Zhalong Wetland, Northeast China. Chinese Geographical Science, 2006, 16(3):265-269. |
[8] | 唐藴, 王浩, 严登华.向海自然保护区湿地生态需水研究.资源科学, 2005, 27(5):101-106. |
[9] | 王昊, 许士国, 孙砳石.扎龙湿地芦苇沼泽蒸散耗水预测.生态学报, 2006, 26(5):1352-1358. |
[10] | 徐东霞, 章光新, 冯夏清.嫩江流域径流量多时间尺度特征分析.资源科学, 2009, 31(9):1592-1598. |
[11] | IPCC. Climate Change 2007:The Physical Science Basis Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change. Solomon S, Qin D, Manning M, Chen Z, Marquis M, Averyt K B, Tignor M, Miller H L. Cambridge:Cambridge University Press, 2007:996-996. |
[12] | Van Vuuren D P, Edmonds J A, Kainuma M, Riahi K, Weyant J. A special issue on the RCPs. Climatic Change, 2011, 109(1/2):1-4. |
[13] | Van Vuuren D P, Edmonds J A, Kainuma M, Riahi K, Thomson A, Hibbard K, Hurtt G C, Kram T, Krey V, Lamarque J F, Masui T, Meinshausen M, Nakicenovic N, Smith S J, Rose S K. The representative concentration pathways:an overview. Climatic Change, 2011, 109(1/2):5-31. |
[14] | Easterling D R, Horton B, Jones P D, Peterson T C, Karl T R, Parker D E, Salinger M J, Razuvayev V, Plummer N, Jamason P, Folland C K. Maximum and minimum temperature trends for the globe. Science, 1997, 277(5324):364-367. |