低丘缓坡建设开发对水源涵养功能的影响

　　摘要 运用InVEST模型并结合实地调研结果，对云南10个典型低丘缓坡建设项目区开发前后生态系统水源涵养功能进行定量测算，比较分析不同类型项目区水源涵养功能的差异。结果表明，开发后项目区平均产水量和平均产水深度均低于开发前，城市建设区的产水总量高于工业区，下降程度也高于工业区;大部分项目区水源涵养总量均低于开发前，仅光华项目区水源涵养总量高于开发前，城市建设区的水源涵养总量高于工业区，工业区水源涵养量下降程度高于城市建设区;导致项目区开发前后水源涵养量变化的原因有土地利用结构、植被、地形、径流、降水、土壤等自然和人为因素，其中受土地利用类型影响较大;低丘缓坡开发导致区域生态环境受到干扰和破坏，但不同的开发模式和功能定位对其影响存在差异。因此在低丘缓坡开发过程中，应采取工程、监管、法律和经济等综合措施，实现低丘缓坡开发建设与生态保护的协调。

　　关键词 低丘缓坡;建设开发;水源涵养;InVEST模型

　　生态系统的水源涵养功能主要体现在改善水文状况、调节区域水分循环等方面[1]。生态系统中的植被如森林、灌丛、草地对水源涵养功能发挥着巨大的作用，其通过林冠层截留、枯枝落叶层持水和土壤层蓄水来调节地表径流和补充地下水，在防治旱涝灾害、净化水质等方面具有至关重要的作用[2]。生态系统服务和权衡的综合评估模型(integrated valuation of ecosystem services and trade-offs，InVEST)从水文角度以中小尺度流域单元为对象，反映不同气候状况、土地利用、地形起伏、土壤性质下的产水量大小，功能模型以地图形式表达了流域的产水能力，间接地对区域的自然价值赋予了适当的衡量标准，而数据与参数的适宜性成为模型结果是否可信的关键[3]。近年来国内学者将InVEST模型广泛运用于北京山区、陕北黄土高原、秦岭地区、横断山区、若尔盖高原等山地[4-8]，均取得较好的应用效果。国外学者将InVEST模型广泛应用于美国、英国等国家和地区的产水量计算中[9-11]。

　　目前，国内关于生态系统服务的研究报道较多，但对整个云南低丘缓坡生态系统服务的研究较少且时序较早。陈武强等[12]研究昆明西华湿地公园水分调节量和价值，赵元藩等[13]估算玉溪市森林生态系统涵养水源量和调节水量价值，杨芳等[14]评估玉龙县天保工程森林生态系统水源涵养能力，都是从宏观上估算了区域的水源涵养能力以及为保护环境创造的价值。云南省山地面积占全省土地面积的90%以上，地质环境复杂，生态环境多样。山地环境具有空间异质性和垂直分异性，海拔、地形起伏和土地覆被在一定的空间尺度上镶嵌复合。水循环受土地利用、气候、地形、土壤、植被等自然和人为因素的复杂影响，如何量化这种水源涵养能力，仍是一个悬而未决的问题[15]。鉴于上述原因，笔者针对气候、地形、土壤、土地利用对山地水循环的影响，采用InVEST产水量模型与地形指数、土壤饱和导水率、流速系数结合的方法，估算研究区水源涵养能力，旨在揭示低丘缓坡综合开发前后水源涵养功能的时空变化与影响因素，为低丘缓坡合理开发和水资源保护提供科学依据。

　　1 资料与方法

　　1.1 研究区概况 该研究在选取研究区时，主要考虑的是功能定位、地理位置、气候类型和开发情况4个因素，具体理由如下：①功能定位，依据功能定位，项目区大致可分为2种类型，一类是以加工、能源、仓储物流、高新技术为主的工业区，且具有特色产业;另一类是集旅游、康养、居住等功能的城市建设区。②地理位置，该研究根据地理区位划分，从滇中选择6个，滇西选择2个，滇南选择2个项目区。③气候类型，云南省大部分地区属于亚热带季风气候，滇南部分地区属于热带季风气候。为了避免气候类型对结果产生较大影响，该研究所选项目区8个属于亚热带季风气候，景洪的光华和景大项目区属于热带季风气候。④开发情况，对项目区进行实地调研，选择开发程度较高的项目区进行研究。

　　研究区位于云南省昆明市、曲靖市、玉溪市、大理白族自治州和西双版纳傣族自治州(表1和图1)。研究区面积在550.98～2 158.40 hm 研究区内土壤类型多样，有水稻土、暗黄棕壤、黄棕壤、红壤、山原红壤、赤红壤、砖红壤、褐红土。植被类型丰富，主要为阔叶林、针叶林、针阔混交林、阔叶灌丛、灌草丛和山地草甸。

　　1.2 InVEST产水量模型

　　InVEST模型基于水量平衡原理，通过研究区的降水、地表蒸发、植物蒸腾、土壤深度、根系深度、地形等数据计算每个栅格的产水量。充分考虑不同土地利用类型下土壤渗透的空间差异和地形等因素对径流的影响，定量估算不同土地利用类型的产水量，结果更为科学。同时，对北京山区、黄土高原和横断山区等山区产水或水源涵养评估，也体现了该模型具备参数调整灵活、空间表达性强等優点。

　　Yjx=1-AETxjPx×Px(1)

　　式中，Yjx为第j类土地利用/覆被上栅格单元x的年产水量(mm);Px为栅格单元x的年平均降水量(mm);AETxj为第j类土地利用/覆被类型上栅格单元x的实际年平均蒸散量(mm)，通过公式(2)计算。

　　AETxjPx=1+ωxRxj1+ωxRxj+1/Rxj(2)

　　式中，Rxj为土地利用类型j上栅格单元x的干燥指数，无量纲，表示潜在蒸发量与降水量的比值，由公式(3)计算;

　　k为蒸散系数，由植被叶面积指数LAI计算获得，见公式(4);ET0为潜在蒸散量(mm/d)，计算公式如式(5)。

　　ωx为修正植被年可利用水量与降水量的比值，无量纲，由公式(6)计算。

　　Rxj=k×ET0Px(3)

　　k=min( LAI3)(4)

　　ET0=0.001 3×0.408×RA×(Tavg+17)×(TD-0.012 3P)0.76(5)

　　ωx=ZAWCxPx(6)

　　式中，RA为太阳大气顶层辐射[MJ/(m2·d)];Tavg是日最高温均值和日最低温均值的平均值(℃);TD是日最高温均值和日最低温均值的差值(℃)。太阳大气顶层辐射用气象站太阳平均总辐射除以50%计算获得[16]。Z为zhang系数;AWCx为植被有效可利用水，由土壤深度和理化性质决定，计算公式如下：

　　AWCx=min(maxSDx，RDx)×PAWCx(7)

　　PAWCx=54.509-0.132Sand-0.003(Sand)2-0.055Silt-0.006(Silt)2-0.738Clay+0.007(Clay)2-2.688OM+0.501(OM)2(8)

　　式中，maxSDx为最大土壤深度;RDx为根系深度;PAWCx为植被可利用水含量;Sand为土壤砂粒含量(%);Silt为土壤粉粒含量(%);Clay为土壤黏粒含量(%);OM为土壤有机质含量(%)。

　　1.3 水源涵养模型 该研究选取研究区均为山地，必须考虑地形对水源涵养功能的影响，因此在产水量的基础上，再利用地形指数、土壤饱和导水率和流速系数对产水量进行修正，得到研究区的水源涵养量。

　　WR=min( 249V)×min( 0.9×TI3)×min( Ks300)×Yjx(9)

　　TI=lg(DASD×PS)(10)

　　lnKs=20.62-0.96×ln(C)-0.66×ln(S)-0.46×ln(OC)-8.43×BD(11)

　　式中，WR为水源涵养量(mm);TI为地形指数，无量纲;Ks为土壤饱和导水率(cm/d)，该研究选用土壤传递函数间接计算得出[17]。V为流速系数，数据无量纲;Yjx为产水量。DA为集水区栅格数量，无量纲;SD为土壤深度(mm);PS为百分比坡度。C为土壤黏粒含量(%);S为土壤砂粒含量(%);OC为土壤有机碳含量(%);BD为土壤容重(g/cm3)。

　　1.4 水源涵养功能影响因素 为定量分析影响因素与不同水源涵养功能等级之间的关系，利用地图代数原理，利用公式(12)对研究区的土地利用类型、坡度与水源涵养功能之间的相关变化信息进行提取。

　　Cij=Rij×100+Sij(12)

　　式中，Cij为各个影响因素图层与水源涵养功能图层进行空间叠加的综合图层，Rij为各个影响因素图层，Sij为水源涵养功能图层。通过对2个不同的专题图层进行空间叠加运算，计算结果既能在空间上反映不同专题图层的等级变化，同时又能得到不同图层之间的相互转换关系。

　　1.5 数据来源及处理

　　1.5.1 产水量数据。

　　水源涵养功能受到多种因素的影响，其计算过程必然关系到多种因素的信息提取与集成利用。利用InVEST模型进行产水量评估需要的输入数据包括土地利用/覆被图、年平均降水量图、年平均潜在蒸散量图、土壤根系深度图、植物可利用水含量(PAWC)、集水区以及生物物理参数表。

　　(1)年平均降水量。降水数据来源于中国气象数据网(https：//data.cma.cn/)的中国地面降水月值0.5°×0.5°格点数据集，为对应研究区开发时段，也为了避免单年数据的低代表性，该研究选取2000—2020年的降水数据，运用Python 3.8.3提取项目区周围格点的月降水量，在ArcGIS中以地形作为协变量进行空间插值获得研究区多年平均降水量栅格图层。

　　(2)年平均潜在蒸散量。太阳辐射数据来源于WorldClim(https：//www.worldclim.org/)，气温数据来源于中国气象数据网(https：//data.cma.cn/)的中国地面累年值日值数据集(1981—2010年)的累年平均日最高气温和日最低气温。估算潜在蒸散量的方法有很多，该研究使用InVEST模型推荐的Modified-Hargreaves。在ArcGIS中进行克里金插值、重采样和裁剪等处理，并运用栅格计算器对太阳辐射、平均日最高气温和最低气温计算得到多年平均潜在蒸散量栅格图。

　　(3)土地利用类型。土地利用类型由资源三号卫星遥感影像解译得到。根据数据分析要求将土地利用类型划分为7个一级地类、14个二级地类，分别是耕地(水田、旱地)、林地(有林地、灌木林地、其他林地)、园地、草地、水域(水库坑塘、滩涂)、建设用地(农村居民点用地、城镇用地、其他建设用地)、未利用地(裸地、荒草地)。同时，利用Google Earth 高分辨率影像(2020年)，对土地利用解译结果和植被覆盖情况进行野外选点验证并访谈当地居民和工作人员，对土地利用数据进行现场复核。

　　(4)土壤深度。土壤数据来源于北京大学城市与环境学院地理数据平台(http：//geodata.pku.edu.cn)的中国土壤数据集。对全国土壤数据进行空间插值后得到云南省土壤根深图，再运用项目区范围数据对云南省土壤根深图进行裁剪得到研究区的土壤根深栅格数据。

　　(5)植被可利用水含量。结合中国土壤数据集中土壤质地的划分，将砂粒、粉粒、黏粒、有机质等质量分数导入土壤有效含水量的经验公式中，在ArcGIS中利用字段算计器求出PAWC，由空间分析和转换工具处理，得到植被可利用含水量图。

　　(6)Z系数。Z系数是表征降水季节性特征的一个常数，其值在1～30，降水主要集中在冬季时，其值越大，降水主要集中在夏季或季节分布较均匀时，其值越小;但在降水量相等的地区，降水事件越频繁，Z值越大。产水量的大小受到Z系数影响，在热带地区Z值取4;在季风气候区Z值取1;在温带海洋性气候区Z值取9[18];该研究区均属于季风气候，故Z值取1。

　　(7)InVEST模型参数。蒸散系数和流速系数数据是参考相关文献与结合InVEST用户指南推荐使用的参考数据以及根据研究区地表植被覆盖实际情况确定。LULC_veg表示采用何种实际蒸散量公式的信息，含有植被土地利用/覆被賦值为 其他赋值为0(表2)。

　　1.5.2 水源涵养数据来源及处理。

　　(1)地形指数。通过项目区资料获得2012年等高线、高程点、项目区范围，通过ArcGIS 10.2创建TIN后转为栅格，得到开发前的DEM;运用2020年资源三号卫星影像立体像对，采用ENVI 5.3提取DEM。通过ArcGIS 10.2的水文分析工具进行填洼、流向、流量等计算步骤，获得流量栅格，根据项目区情况确定阈值，利用Con条件函数筛选主要径流，经过盆域分析得到流域图。再采用面转栅格工具得到流域单元集水区栅格的数量;通过 ArcGIS中3D分析提取DEM中的百分比坡度，此时会得到坡度为0的栅格，为了不影响后面的计算，将百分比坡度为0的点替换成0.000 0 由于该值远小于地形指数的计算位数，误差可忽略不计。最后综合百分比坡度、流域单元集水区栅格数量以及土壤深度可得到流域的地形指数。

期刊VIP网，您身边的高端学术顾问

文章名称： 低丘缓坡建设开发对水源涵养功能的影响

文章地址： http://www.qikanvip.com/nengyuan/60646.html