Spatio-temporal variation in landscape leakiness and the influencing factors in the riparian vegetation buffer zone of the Beijiang River
XUShanshan1,2,, ZHAOQinghe1,2,, CAOZihao1,2, ZHANGYifan1,2 1. Key Laboratory of Geospatial Technology for the Middle and Lower Yellow River Regions, Ministry of Education, Kaifeng 475004, China2. College of Environment and Planning, Henan University, Kaifeng 475004, China 通讯作者:通讯作者:赵清贺,E-mail:zhaoqinghe@henu.edu.cn 收稿日期:2017-01-9 修回日期:2017-09-28 网络出版日期:2018-06-25 版权声明:2018《资源科学》编辑部《资源科学》编辑部 基金资助:国家自然科学基金项目(41301197)河南省高等学校重点科研项目(18A170004)河南省高校科技创新团队支持计划(16IRTSTHN012) 作者简介: -->作者简介:徐珊珊,女,河南平顶山人,博士生,主要从事流域景观格局与生态过程研究。E-mail:xus199@163.com
关键词:景观渗透性;植被和地形;水土保持;河岸植被缓冲带;时空演变;影响因素;北江干流 Abstract Landscape leakiness describes the extent to which landscapes have lost their capacity to regulate rainwater and soil nutrients, which can provide an approach to assess landscape function. In this study, taking the riparian vegetation buffer zone of the Beijiang River as a case study, the leakiness index (LI), an indicator for reflecting landscape leakiness, was applied to evaluate the soil and water conservation function of the riparian vegetation buffer zone. Based on the DEM and vegetation data derived from remote sensing images from 1995-2015, we explained the calculation process of LI and verified the applicability of LI in the study area using the USLE model. We analyzed spatial-temporal variation in landscape leakiness and associated factors. The results showed that LI presented better applicability in the riparian vegetation buffer zone of the Beijiang River according to a significantly positive correlation between LI and soil erosion modulus. Temporal and spatial variation in landscape leakiness was obvious and generally showed an increasing tendency during 1995-2015. While landscape leakiness in several typical areas was higher with fluctuated change, such as Shaoguan, Nanxiong and Qingyuan. Land use changes, especially forest and construction land were major factors affecting landscape leakiness variation. Specifically, LI presented inverse variation with changes in the proportion of forest land area, and the same tendency with construction land. Results from this study can provide a reference for the applicability of LI and the evaluation of water and soil conservation function of riparian vegetation buffer zone in subtropical areas.
Keywords:leakiness index;vegetation and terrain;soil and water conservation;riparian vegetation buffer zone;spatio-temporal variation;influencing factors;Beijiang River -->0 PDF (6957KB)元数据多维度评价相关文章收藏文章 本文引用格式导出EndNoteRisBibtex收藏本文--> 徐珊珊, 赵清贺, 曹梓豪, 张祎帆. 北江干流河岸植被缓冲带景观渗透性时空变化及其影响因素[J]. 资源科学, 2018, 40(6): 1267-1276 https://doi.org/10.18402/resci.2018.06.16 XUShanshan, ZHAOQinghe, CAOZihao, ZHANGYifan. Spatio-temporal variation in landscape leakiness and the influencing factors in the riparian vegetation buffer zone of the Beijiang River[J]. RESOURCES SCIENCE, 2018, 40(6): 1267-1276 https://doi.org/10.18402/resci.2018.06.16
1 引言
水土流失引起的环境问题已严重威胁到区域生态平衡,给人们的生产和生活带来严重危害,成为世界关注的焦点[1, 2]。河岸植被缓冲带作为流域陆生生态系统和水生生态系统的生态过渡带,通过阻挡地表径流、降低径流流速、吸收养分等,在拦截过滤悬浮颗粒物、阻控非点源污染、保持水土等方面发挥着重要作用[3]。但是,近年来在城市化和现代农业发展的影响下,大量的河岸植被缓冲带被破坏并急剧减少,造成河岸带水土保持功能退化,严重威胁流域生态系统功能的完整性[4]。因此,对河岸植被缓冲带水土保持状况进行研究将有助于河岸带的保护和管理,从而确保其生态功能的可持续发挥。 针对水土保持研究,以往常用的方法有水文法、遥感解译法、模型估算法、径流小区法等[5],此类方法为水土保持的定性和定量研究奠定了基础。关于河岸植被缓冲带水土保持的研究多以RS、GIS、GPS技术和相关模型的结合为主[6],如SWAT模型(Soil and Water Assessment Tool Model)、REMM模型(Riparian Ecosystem Management Model)和“源-汇”景观模型(“Source-Sink” Landscape Model)等[7]。同时,亦有不少研究通过小区试验研究河岸植被缓冲带对径流量、悬浮颗粒物以及N和P的阻控能力[8]。近年来,随着相关研究的深入,一些新方法和技术被逐渐应用到水土保持研究领域中,如地形测量、核素示踪、沉积泥沙反演、格局与过程耦合指数等,在改进研究方法的同时,也提高了研究的时效性和准确性[5, 9]。其中,格局与过程耦合指数,如连通性指数(Connectivity Index,CI)、径流长度指数(Flowlength,FL)和景观渗透性指数(Leakiness Index,LI)等[10,11,12],通过分析径流在景观单元之间运动的连续性、建立景观连接度与水土流失之间的关系,能有效地将景观格局与过程相结合从而反映景观的功能特征,而备受关注。在此类耦合指数中,景观渗透性指数(LI)能够定量化地评价景观水土保持功能,具有可靠、操作简便、易于使用等优点[13],已在干旱半干旱景观中(如澳大利亚牧场)得到应用和验证[11, 14]。但是,该指数在亚热带尤其是河岸植被缓冲带的应用还鲜有报道。 基于此,本文以北江干流河岸植被缓冲带为例,结合USLE(Universal Soil Loss Equation)模型,在验证景观渗透性指数适用性的基础上,对研究区景观渗透性的时空动态进行分析,并根据土地利用状况分析其动态的影响因素,以期为北江干流河岸植被缓冲带水土保持评价提供一定的方法参考。
2 研究区概况、数据来源与研究方法
2.1 研究区概况
北江是珠江流域的第二大水系(23°0′N—25°41′N, 111°52′E—114°41′E),是广东省四大河流之一。上游浈水发源于江西信丰县石碣大茅山,在广东韶关市与发源于湖南临武县的武水汇合后称为北江,本研究范围为北江干流10 km河岸植被缓冲带,区域内共涉及14个县市,其分布状况如图1所示。该区属亚热带季风气候,夏季高温多雨,冬季低温少雨,年均气温在18℃~21℃,最高气温出现在7、8月份[15];年均降雨量为1600~1800 mm,降雨多集中在汛期(4—9月),占全年降雨量的70%以上[16]。北江年平均径流量为482亿 m3,水资源丰富。北江流域地势北高南低,地貌以山地和丘陵为主,地形陡峭,水土流失严重。植被类型主要为马尾松、桉树林、竹林以及稀树草坡[17]。研究区土壤以红壤为主,养分充足且保水能力强,有利于植物和农作物的生长与发育,导致河岸被严重开垦。另外,北江流域人口增长快且增幅高[18],流域内森林砍伐、基础建设、开山造田等人类活动频繁,土壤侵蚀加剧[19],水土流失是其面临的主要生态问题之一。 显示原图|下载原图ZIP|生成PPT 图1研究区示意 -->Figure 1The study area -->
2.2 研究方法
2.2.1 景观渗透性指数计算方法 水土流失与景观格局关系密切,景观格局可通过植被、地形等景观要素影响水土流失过程,而水土流失过程可通过土壤养分、水分等资源的再分配影响景观格局[20]。因此,建立景观格局与水土流失过程之间的相互联系是深入研究水土保持的有效方法。景观渗透性指数(LI)是在DLI(Directional Leakiness Index)和CDLI(Cover-based Directional Leakiness Index)基础上加以改进而来[11, 21, 22],相比于后者,其在丘陵、山地等复杂地形中的适用性较强并得到相关研究的证实[14]。该指数通过地表径流在景观要素(植被和地形)之间的连通性,表征景观水土保持能力。植被斑块和洼地,能够有效减缓地表径流、提高土壤抗蚀性、增加入渗,从而达到保持水土的功效,因此在LI计算过程中作为“汇”景观;而裸土斑块和坡度较大的坡面,对径流阻力小、增加水流流速,从而加剧土壤侵蚀,因此被作为“源”景观。“源-汇”景观之间的连通性是影响水土流失的重要因素,在径流从源区到达目标区域的路径上,LI充分考虑了“源-汇”景观单元在水流过程中的作用,最终达到定量评价景观水土保持功能的效果。 径流从源区到达目标区域的路径上,植被斑块作为抑制侵蚀的景观要素,影响径流连通性,进而影响水土流失。植被盖度是表征植被斑块水土保持潜力的有效指标,因此,植被盖度是景观渗透性指数计算过程中的重要参数之一。在众多植被指数中,归一化植被指数(NDVI)被认为是监测植被生长状态及植被盖度的最佳指示因子并得以广泛应用[23]。本研究中植被盖度通过NDVI计算得到,公式如下: (1) 式中IR和R分别表示植被在遥感影像近红外波段和红光波段上的反射率。在NDVI基础上,植被盖度的计算方法为: (2) 式中C为单个像元植被盖度;NDVIsoil和NDVIveg分别为最小和最大NDVI值,即无植被覆盖或裸土像元的NDVI值和纯植被覆盖的像元的NDVI值。 基于提取的植被盖度,结合DEM数据,根据公式(3)计算景观渗透性指数: (3) 式中Lcal为沿水流方向由逐步累积汇流算法计算的景观渗透值;Lmin和Lmax分别为最小(植被覆盖度为100%时的Lcal)和最大渗透值(植被覆盖度为0时的Lcal);k为LI与植被覆盖之间关系的衰减曲线的斜率,根据Ludwig等的研究,并结合北江干流实际情况,本研究中k取值为5[11, 13]。LI取值范围值在[0,1],其值越小,景观水土保持功能越强;其值越大,景观水土保持功能越弱。Lcal采用多流向累积算法计算。多流向算法的设计及计算多以DEM为数据源,在计算过程中为避免伪地形造成的误差,如由生成DEM时的内插误差、低质量的数据输入等原因造成的伪洼地对径流路径的影响,在数据处理时,需对DEM进行填洼处理,图2a为填洼后的DEM数据,在此基础上提取每个栅格的水流流向(图2b)。Lcal的算法如下: 显示原图|下载原图ZIP|生成PPT 图2基于DEM的水流流向矩阵示意 -->Figure 2The diagram of flow direction matrix based on DEM -->
如图3所示,1995—2015年间,研究区景观渗透性指数呈现先上升后下降的趋势,并且变化明显,值分别为0.09,0.12,0.21,0.24和0.23,最低值出现在1995年,最高值出现在2010年。 显示原图|下载原图ZIP|生成PPT 图31995—2015期间北江河岸植被缓冲带景观渗透性 -->Figure 3Landscape leakiness of the Beijiang riparian vegetation buffer zone from 1995 to 2015 -->
受流域内水源涵养林的建设以及退耕还林政策等因素影响,区域内大于25°的坡耕地基本上实现了退耕还林[19],区域植被覆盖度增加,水土保持功能增强,使1995年和2000年的渗透性指数值偏低。另外,随城市化进程的加快,北江流域作为广东省重要的金属矿资源分布区,如韶关段有“有色金属之乡”的美誉,凡口铅锌矿等均为中国著名的大型金属矿山[28]。但由于矿产资源的过度开采、旅游业的开发以及水利水电等工程的建设等人为干扰的作用下,大面积的植被和表土受损,水土流失严重[17],致使2005—2015年间景观水土保持功能下降,渗透性指数有所升高。根据20年间景观渗透性指数的明显变化可知该指数在研究区具有一定的敏感性。 为了一步验证LI在研究区的适用性,本研究运用USLE模型对该指数进行验证。如图4所示,1995—2015年间,研究区侵蚀模数与LI变化趋势相似,并呈现显著正相关关系(P<0.05),表明LI具有指示河岸植被缓冲带水土流失量的潜力,亦能有效评价河岸植被缓冲带的水土保持状况,在研究区具备一定的适用性。 显示原图|下载原图ZIP|生成PPT 图41995—2015期间北江河岸植被缓冲景观渗透性与土壤侵蚀模数关系 -->Figure 4The relationship between landscape leakiness and erosion modulus of the Beijiang riparian vegetation buffer zone from 1995 to 2015 -->
3.2 研究区景观渗透性指数时空分异特征
基于验证结果,为进一步分析北江河岸植被缓冲带景观渗透性指数的时空分异状况,本研究按照行政边界对研究区进行空间划分。根据研究区5个时期景观渗透性指数的变化,将LI分为6个等级,各分区景观渗透性指数的空间分布如图5所示。 显示原图|下载原图ZIP|生成PPT 图51995—2015年北江河岸植被缓冲带景观渗透性指数的时空分异示意 -->Figure 5Spatio-temporal variation of landscape leakiness index in the Beijiang riparian vegetation buffer zone from 1995 to 2015 -->
3.3.1 1995—2015年研究区土地利用变化分析 表1为不同时期研究区土地利用变化。1995—2015年,农田、水体和未利用地面积变化不大,建设用地和林地的面积变化相对较大,其中建设用地由最低2000年的335.82 km2增加至2015年的510.48 km2,林地的最低值为1995年的4629.69 km2,最高值为2005年的4817.77 km2。20年间,林地和建设用地的面积变化比率较大,分别为3.58% 和18.29%。基于林地和建设用地的大幅度变化,同时考虑这两种土地类型对水土流失作用,本文重点分析林地和建设用地面积的变化对时空动态较大的分区的景观渗透性指数的影响。 Table 1 表1 表11995—2015年北江河岸植被缓冲带土地利用类型面积变化 Table 1Area changes in land use types in the Beijiang riparian vegetation buffer zone from 1995 to 2015 (km2, %)
1995年
2000年
2005年
2010年
2015年
变化百分比
农田
2 320.28
2 418.20
2 336.81
2 326.42
2 332.68
0.53
水体
548.54
544.74
566.36
566.36
556.22
1.40
林地
4 629.69
4 810.02
4 817.77
4 765.00
4 795.51
3.58
建设用地
431.56
335.82
410.30
475.79
510.48
18.29
未利用地
17.87
13.14
15.38
25.15
18.46
3.30
新窗口打开 3.3.2 林地和建设用地面积的变化对渗透性指数的影响 各分区的林地面积和建设用地面积占该分区总面积的比例如图6所示,1995—2015年间,位于上游的南雄、始兴、韶关的林地面积占比的变化趋势相似,呈现较大的波动状态,均以2000年面积占比最小,分别为41%、51%和39%;南雄和韶关均以2015年面积占比最大,分别为51%和49%,始兴2005年面积占比最大,为68%。韶关作为粤北地区的政治、经济、交通、文化中心,人口增长快,工业基础雄厚,旅游业发达,经济的发展加速了建设用地的扩张,导致韶关的渗透性值增高,且明显高于其他地区。位于下游的清远林地面积占比呈现先上升后下降的变化,四会、三水、花都变化趋势相似,且与上游的南雄、始兴和韶关相反。林地面积占比的增加或减少直接影响着植被盖度的增减,进而影响景观渗透性的变化。对比各个时期景观渗透性指数的变化可知,7个典型分区景观渗透性值的变化与林地面积占比的变化呈相反的趋势。 显示原图|下载原图ZIP|生成PPT 图61995—2015年典型区域林地面积和建设用地面积占比与景观渗透性指数的变化 -->Figure 6Variations in landscape leakiness index and the proportion of forest land and construction land in the typical areas from 1995 to 2015 -->
1995—2015年间,上游的南雄和韶关建设用地面积占比波动性变化相似。下游的清远、四会、三水、花都呈现相同的波动性变化,且最高值均出现在2015年,这与下游地区地势平坦、水源丰富、有利于人口的居住、产业的发展有关。各个时期建设用地面积的变化表明,7个典型分区建设用地面积占比逐渐增大,景观渗透性与其呈相似的变化,表明建设用地占比的增加可通过降低植被盖度、增加径流产生与汇集等对景观渗透性产生影响。 相关分析表明(表2),林地面积比例与渗透性指数呈现负相关关系,建设用地面积比例与渗透性指数呈现正相关关系。其中除始兴和清远外,其他各区的林地面积比例与渗透性指数均显著负相关,这可能与始兴位于流域上游、受降雨侵蚀等影响较大有关,而清远位于流域下游,人为干扰严重,林地面积年际波动较大致使其与渗透性指数的相关性相对较低。区域林地面积的增加,可有效减少水土流失。除韶关、四会和三水外,其他各区的建设用地面积占比与渗透性指数均显著正相关,表明建设用地面积的增加可引起更严重的水土流失。 Table 2 表2 表21995—2015年典型区域林地面积和建设用地面积占比与景观渗透性指数的相关性 Table 2Correlation between landscape leakiness index and the proportion of forest land and construction land in the typical areas from 1995 to 2015
景观渗透性指数虽是评价景观水土保持功能的有效方法,但鲜见在热带亚热带地区的应用与分析。本文在验证景观渗透性指数在亚热带的北江干流河岸植被缓冲带的适用性的基础上,进一步分析了其时空分异与影响因素。为该指数未来在河岸植被缓冲带的应用提供了一定理论基础的同时,也为北江干流河岸植被缓冲带的管理与相关水土流失治理的目标区域提供了参考。主要结论如下: (1)1995—2015年间,研究区景观渗透性指数变化较明显,呈逐渐增加的趋势。相同时期的侵蚀模数与景观渗透性指数表现出显著的正相关(P<0.05),表明景观渗透性指数在研究区具有一定的适用性。 (2)在研究区不同分区景观渗透性时空分异的分析中,多数地区变化并不明显。但韶关、南雄、清远等地区呈现波动性变化,是景观渗透性指数时空差异较大的典型分区。 (3)1995—2015年间,整体上典型分区内林地面积所占比例的变化与景观渗透性指数呈现负相关关系,而建设用地与景观渗透性指数呈现正相关关系。林地和建设用地均能通过增强或减弱地表拦截径流及泥沙的能力进而影响景观的水土保持状况,表明林地和建设用地面积的变化是影响景观渗透性指数变化的重要因素。 The authors have declared that no competing interests exist.
[Wang XK, Ouyang ZY, XiaoH, et al. Distribution and division of sensitivity to water-caused soil loss in China [J]. Acta Ecologica Sinica, 2001, 21(1): 14-19. ] [本文引用: 1]
[2]
PimentelD.Soil erosion: a food and environmental threat [J]. Environment, Development and Sustainability, 2006, 8(1): 119-137. [本文引用: 1]
[Rao LY, Cui JG.Research advances on the eco-hydrological functions of riparian buffer [J]. Science of Soil and Water Conservation, 2008, 6(4): 121-128. ] [本文引用: 1]
[4]
BastieneN, KirstukaJ.Riparian Buffer Strips-Foreign Experience and Readjustment to Lithuanian Conditions [C]. Vilnius: Environmental Engineering The 8th International Conference, 2011. [本文引用: 1]
[TangQ, Bao YH, He XB, et al. New soil erosion monitoring methods and techniques [J]. Science of Soil and Water Conservation, 2011, 9(2): 11-18. ] [本文引用: 2]
[ZengY.Construction Technology of Vegetation Buffer Zones Based on "Source-Sink" Model: A Case Study of Kaixian in Chongqing[D] . Wuhan: Huazhong Agricultural University, 2014. ] [本文引用: 1]
[8]
SrivastavaP, SinghS.Conservation of soil, water and nutrients in surface runoff using riparian plant species [J]. Journal of Environmental Biology, 2012, 33(1): 43-49. [本文引用: 1]
[9]
LiuY, Fu BJ, Lü YH, et al. Linking vegetation cover patterns to hydrological responses using two process-based pattern indices at the plot scale [J]. Science China Earth Science, 2013, 56(11): 1888-1898. [本文引用: 1]
[10]
BorselliL, CassiP, TorriD.Prolegomena to sediment and flow connectivity in the landscape: a GIS and field numerical assessment [J]. Catena, 2008, 75(3): 268-277. [本文引用: 1]
[11]
Ludwig JA, Bastin GN, Chewings VH, et al. Leakiness: a new index for monitoring the health of arid and semiarid landscapes using remotely sensed vegetation cover and elevation data [J]. Ecological Indicators, 2007, 7(2): 442-454. [本文引用: 5]
[12]
Mayor ÁG, BautistaS, Small EE, et al. Measurement of the connectivity of runoff source areas as determined by vegetation pattern and topography: a tool for assessing potential water and soil losses in drylands [J]. Water Resources Research, 2008, 44(10): 2183-2188. [本文引用: 1]
[LiuY, Wu BF, ZengY, et al. Assessment of the impacts of soil erosion on water environment based on the integration of soil erosion process and landscape pattern [J]. Chinese Journal of Applied Ecology, 2013, 24(9): 2581-2589. ] [本文引用: 4]
[14]
BastinG, AbbottB, KarfsB, et al. Landscape Leakiness in the Grazed Rangelands of the Burdekin [C]. Queensland: Australian Rangeland Society, 15th Biennial Conference, 2008. [本文引用: 2]
[Gao CY, Zhao QH, LiuQ.Fractal characteristic of soil particle size under different vegetation types in riparian zone of the main stream of Beijiang River [J]. Research of Soil and Water Conservation, 2016, 23(3): 37-42. ] [本文引用: 1]
[LuoL, Wu GC, Wu JN.Analysis on the climatic characteristics of rainfall in flood season of Beijiang River Basin from 1965 to 2010 [J]. Guangdong Meteorology, 2012, 34(3): 13-15. ] [本文引用: 1]
[Wang JG, Jin PW, Liu CQ, et al. Analysis on soil and water losses in upper reaches of Beijiang River basin based on Remote Sensing [J]. Pearl River, 2012, 34(3): 13-15. ] [本文引用: 2]
[LiY, Zhang PF.Analysis of change characteristics and variability of runoff in Beijiang river basin under influence of human activity [J]. Journal of Water Resources and Water Engineering, 2014, 25(2): 61-65. ] [本文引用: 1]
[Zou ZX, Xu XM.Problems and solutions in the sustainable utilization of water resources in Guangdong Beijiang River Basin [J]. Heilongjiang Science and Technology Information, 2008, 11(7): 119-119. ] [本文引用: 2]
[Fu BJ, Xu YD, Lv YH.Scale characteristics and coupled research of landscape pattern and soil and water loss [J]. Advances in Earth Science, 2010, 25(7): 673-681. ] [本文引用: 1]
[21]
Ludwig JA, Eager RW, Bastin GN, et al. A leakiness index for assessing landscape function using remote sensing [J]. Landscape Ecology, 2002, 17(2): 157-171. [本文引用: 1]
[22]
Ludwig JA, Eager RW, Liedloff AC, et al. A new landscape leakiness index based on remotely sensed ground-cover data [J]. Ecological Indicators, 2006, 6(2): 327-336. [本文引用: 1]
[Wei HB, LiR, Yang QK.Research advances of vegetation effect on soil and water conservation in China [J]. Acta Phytoecologica Sinica,2002, 26(4): 489-496. ] [本文引用: 1]
[Cai CF, Ding SW.Study of applying USLE and geographical nformation system IDRISI to predict soil erosion in small watershed [J]. Journal of Soil and Water Conservation, 2000, 14(2): 19-24. ] [本文引用: 1]
[Zhou HR, Xiao DN, Zhou KF.Corridor effect of landscape spatial pattern change in arid area-taking the corridor area of the middle and lower reaches of Tarim River as an example [J]. Chinese Science Bulletin, 2006, 51(S1): 66-72. ] [本文引用: 1]
[ZhouT, Peng SL, Ren WT.Influence of landscape pattern changes on the restoration of stream in Dongjiang River riparian buffer [J]. Acta Ecologica Sinica, 2009, 29(1): 231-239. ] [本文引用: 1]
[National Earth System Science Data Sharing Infrastructure, National Science & Technology Infrastructure of China. The 1: 4000 000 Digital Soil Map of China (the Chinese Soil Taxonomy (CST) [EB/OL]. (2017-01-02)[2017-01-09]. ]URL [本文引用: 1]
[Wu HJ, Yuan SF, Guo LL.Environment-economy analysis on mineral resources development in North River valley [J]. China Environmental Protection Industry, 2013, (5): 63-66. ] [本文引用: 1]