白洋淀—大清河流域位于京津冀腹地，是太行山区与渤海湾水资源连通的廊道。自2017年中共中央、国务院决定在河北省设立国家级新区——雄安新区以来，新区的建设与城市发展对水资源的需求量显著增加，流域水资源短缺与不合理利用、水土流失问题亟待解决。白洋淀作为北方的典型湿地，是雄安新区重要的水源地和生态屏障。白洋淀—大清河流域上游山区森林的水源涵养功能是该区重要的生态系统服务，该区水源涵养的空间分布状况与自然条件对生态系统服务的影响受到广泛关注。《全国重要生态系统保护和修复重大工程总体规划(2021—2035年)》将白洋淀—大清河流域的水源涵养林与水土保持工程列入“十四五”建设重要内容。森林植被有利于改良流域土壤，拦截降水、调节坡面径流，消减洪峰，延缓地下径流汇流速度，减少区域水土流失量。定量评估白洋淀—大清河流域产水和水源涵养功能，探明流域森林重点建设和优先保护区，通过在水源涵养能力提升潜力较大的地块优化水源涵养林林分结构配置，精准替换高耗水植被、提升森林覆盖率等方式，提高白洋淀—大清河流域水源涵养能力，解决因植被耗水量过大、生态环境脆弱导致的京津冀地区城市发展水量短缺、环境恶化等问题，以保证雄安新区水资源的持续供应。结合该区生态与生产生活共融、“以水定城”“生态护城”^[1]的发展理念，探寻人口密集区经济—社会—生态可持续发展的新出路，以在京津冀地区建设一个高标准的绿色城区。

应用InVEST模型可评估多种生态功能，能够利用情景化模拟分析预测变化^[2]，分析生态系统服务能力的空间分布，在变化监测、使用成本和开源性等诸多方面有显著优势^[3]。陈龙等^[4]进行澜沧江流域的水源涵养功能研究，为澜沧江水域水文状况的评估与有效管理提供依据；刘宥延等^[5]定量分析西安市丘陵区的水源涵养分布格局与冷热点空间特征，帮助黄土高原地区的相关部门对自然资源进行管理和规划。

笔者定量评估白洋淀—大清河流域水源涵养量，进行水源涵养功能的等级划分和空间可视化表达，结合自然条件和土地利用现状，探明重点建设和优先保护区，筛选水源涵养提升潜力较大的区域，从而为制定京津冀区域水源涵养保育林规划和生态修复决策等提供参考。

白洋淀—大清河流域处于海河流域中部(E 113°39′~117°34′，N 38°10′~40°102′)，是海河五大水系之一^[6-7]，总面积为4万3 060 km²，山区约为1万8 659 km²，丘陵平原区约2万4 401 km²，各为河系面积的43.33%和56.67%。河系内地形西北高东南低^[7]。该区多为温带大陆性季风气候，属半干旱地区，季节特征显著^[8]；年平均气温7.6~13.1 ℃；无霜期129.0~199.2 d；年日照时间2 518.7~2 924.8 h；常年降雨500~600 mm且集中在汛期(7—9月)，年内时空分布不均。研究区土壤以棕壤土和褐土为主，山区一般土层有机质含量低，质地疏松，土壤抗侵蚀能力较弱^[9]。研究区的地理位置如图 1所示。

图 1(Fig. 1)

图 1 研究区位置分布图 Fig. 1 Location map of the study area

本研究所需数据主要为年平均降雨量、年平均潜在蒸散、土壤深度、植被可利用水含量、土地覆被类型、流域边界、Zhang系数^[10]、地形指数、土壤饱和导水率和流速系数。

其中降水、温度、日照时间等数据均来自中国气象数据网1981—2019年数据表。多年降水数据用流域内外共74个站点数据及GPM卫星降水数据做补充，利用线性回归残差克里金插值进行精细化插值计算^[11]得到。年平均潜在蒸散量使用Modified-Hargreaves法计算并精细化插值。土壤质地数据来源为中科院南土壤所提供的1∶100万土壤数据，使用SPAWHydrology模型中的土壤属性计算，求得PAWC和饱和导水率。土地利用数据和DEM数据均来源于地理空间数据云，土地利用数据通过2019年Landsat-8遥感影像解译与山区40个县林业资源二类调查数据^[12]结合获得。流域边界运用ArcGIS水文分析得到；Zhang系数是经验常数，白洋淀—大清河流域属于半干旱温带大陆性季风气候区，季节特征显著，夏季降水集中，因此Zhang系数选取3~12的不同值模拟，随着Zhang系数值的增大，模型产水量降低。根据研究区年降水次数^[13]，同时结合海河流域近50年(1957—2000年)平均河川径流量观测数据分析对照，当Zhang系数值为8时，模型模拟产水量与多年观测数据最接近。流速系数参照InVEST模型手册^[14-15]进行设定。

笔者基于InVEST模型Water Yeild模块^[16]计算流域产水量，数据获取方便，参数设定灵活，综合不同植被覆盖条件下的蒸散发、不同土壤质地的土壤渗透特性、地形、地表粗糙度等因素对地表径流产生的影响，能较好地考虑到影响产水的主要因素，表达流域产水量的分布状况。

模型原理基于Budyko水热耦合平衡假设^{[2, 17-18]}计算产水量，结合流速系数、土壤饱和导水率(mm/d)、地形指数求出径流在区域上存留时间^[19-20]，求出水源涵养量。模型相关参数如图 2所示。

图 2(Fig. 2)

图 2 模型相关参数分布图 Fig. 2 Distribution of model related parameters

计算结果显示，白洋淀—大清河流域产水总量为35.8亿m³，产水量介于-0.14~444.96 mm之间，平均产水量为83.91 mm。由图 3可知，研究区产水空间分布规律表现为自西北太行山区向中部浅山丘陵区缓慢增加，再自中部浅山丘陵区往东南平原区逐渐减少。产水高值位于中部地区，该地雨量充沛且裸岩和灌木林占地面积广阔，发生降雨事件后极容易在裸露岩石上产生径流，由于地形不平坦，人类活动影响较平原区少，故产水量高；中值区位于东南部平原地区，该地地势平坦，以耕地和建设用地为主，人类活动影响最大，产水量高处多为城市建设用地，降水后易形成地表径流；低值区位于西北部低山区，常年降水少，该地地势起伏大且林木繁多，基本上不受限于人类活动。

图 3(Fig. 3)

图 3 白洋淀—大清河流域产水量空间分布图 Fig. 3 Spatial distribution of water yield in Baiyangdian-Daqinghe River Basin

计算结果显示，白洋淀-大清河流域水源涵养量达到1.7亿m³，深度介于0~114.94 mm之间。由图 4可知，白洋淀大清河流域水源涵养能力的空间的差异性表现为自西北太行山区向中部浅山丘陵区缓慢增加，再自中部浅山丘陵区往东南平原区逐渐减少，变化趋势与产水的变化规律基本一致。结合地统计分析方法统计分析流域的各类土地利用的水源涵养(表 1)可以看出，涵养能力由高到低排序为灌木>乔木>平原林地>耕地>草地>建设用地>裸土>裸岩。山区灌木林地和乔木林地的水源涵养功能最强，其次是平原区林地。

图 4(Fig. 4)

图 4 白洋淀—大清河流域涵养量空间分布图 Fig. 4 Spatial distribution of water conservation in Baiyangdian-Daqinghe River Basin

表 1(Tab. 1)

表 1 白洋淀—大清河流域不同地类涵养量 Tab. 1 Conservation scale of different land types in Baiyangdian-Daqinghe River Basin 土地利用类型 Landuse type 流域面积 Watershed area/km2 总涵养量 Total conservation/106 m3 占总涵养量比例 Proportion in total conservation/% 单位面积涵养深度 Conservation depth per unit area/mm 乔灌混交林Mixed forest 7 104.67 31.93 18.61 4.49 建设用地Construction land 8 900.62 15.34 8.94 1.72 耕地Farmland 12 924.01 51.62 30.09 3.99 乔木Arbor 5 471.24 39.62 23.09 7.24 裸土Bare soil 242.19 0.38 0.22 1.58 灌木Shrub 3 153.07 25.72 14.99 8.16 裸岩Bare rock 3 688.89 5.75 3.35 1.56 草地Grassland 203.00 0.54 0.31 2.65

表 1 白洋淀—大清河流域不同地类涵养量 Tab. 1 Conservation scale of different land types in Baiyangdian-Daqinghe River Basin

中西部丘陵低山区主要以灌木林和乔木林为主。该区雨量充沛，降雨发生后，较缓的地势使径流流动速度放缓，增加区域内的松散土壤水分渗入的时间，地表植被根系也能固持土壤。灌木、乔木的涵养能力强，但乔木林植被丰茂，对水分的消耗也大。因此灌木林地的水源涵养能力要高于乔木林地。东南部平原区主要为耕地和建设用地，林地面积较少，且平原区林地多为经济林或人工纯林，区域的降水量明显低于山区降水，平原林地的水源涵养能力低于山区林地。耕地包括山区坡耕地和平原旱地水浇地，但考虑到山区坡耕地面积较少且分散，与平原区耕地归并为一类。因为耕地作物生长周期短，根系不如灌木和乔木发达，短期耗水量大，对降水的截留能力弱且缺少枯枝落叶层辅助吸收水分，因此水源涵养能力要低于灌木、乔木。研究区内草地主要分布于西北部中山区，地势较高，坡度相对陡峭，土壤难以附着，植被较少，不利于拦截降水，并且当地降水量相对较少，涵养量就低于林地和耕地。

根据《国家生态功能红线划定技术指南》中生态系统服务重要性分级方法，在ArcGIS软件中对研究区水源涵养量处理，按照几何间隔原理划分为5级(图 5)。研究区水源涵养功能受到气候条件、地形地貌和土地利用条件的综合影响，通过对水源涵养量重要性评估，可以在规划决策时突出潜力区，为后续优先保护提供参考。中部丘陵区雨量充沛，受到人类活动影响小且裸土裸岩和灌木林占地面积广阔，可通过优化植被种类与结构，增加灌木面积，增强小流域的涵养能力。平原区多为耕地和建设用地，地势平坦，受到人类活动影响大，退耕还林还草成本高，易对当地经济社会发展产生不利影响。综合考虑流域的自然因素与社会经济因素，在水源涵养潜力较大的小流域，进行开荒造林，优化植被，增强小流域的涵养能力，助力雄安新区的规划建设。

图 5(Fig. 5)

图 5 白洋淀—大清河流域水源涵养重要性及潜力区分布图 Fig. 5 Distribution map of important areas and potential area of water conservation in Baiyangdian-Daqing River Basin

1) 白洋淀—大清河流域产水总和为35.8亿m³，水源涵养总量为1.7亿m³，平均产水量是83.91 mm。大清河流域产水量及水源涵养量的空间布局具有显著差异性。

2) 土地利用类型不同的区域的水源涵养能力存在差异，由高到低排序为灌木>乔木>平原林地>耕地>草地>建设用地>裸土>裸岩。中部丘陵区雨量充沛，受到人类活动影响小且裸土占地面积广阔，后续水源涵养提升潜力较大。

白洋淀—大清河流域是雄安新区的水源地，新区的生态城市建设的关键一环。要通过流域水源涵养功能的综合评估，优化流域水资源配置来增强水资源可持续利用，提高生态系统的服务效率，加快推进清洁流域治理和山水林田湖草人综合流域建设。后续研究需要综合考虑海河流域的上下游流域水资源的关系及地下水与地表水的差别，采取其他方法减少人类直接消耗水量的影响；增加模型对中国自然生态环境的适应性，改进模型参数，结合区域的实际径流观测数据作进一步验证分析等；另外应重视流域内的植被，可通过模拟不同植被结构下水源涵养功能的强弱，得到最优状态下的林木结构和植被物种，积极推进水源涵养林工程更好更快发展。