板栗树种在中国水土流失区的分布及其环境因子

唐兴港, 袁颖丹, 张 星, 张金池

(1.南方现代林业协同创新中心 江苏省水土保持与生态修复重点实验室南京林业大学 林学院, 江苏 南京 210037; 2.扬州大学 园艺与植物保护学院, 江苏 扬州 225009)

改革开放以来中国的快速发展立足于资源与环境的消耗，这使得当前环境与发展的矛盾逐渐成为新的瓶颈[1]。作为生态环境决定因素和生产生活物质基础的水土资源面临着严峻的挑战[2]。据统计，中国2/3以上的国土面积是山地丘陵，现有水土流失面积2.74×106km2，每年新增治理水土流失面积约5.50×104km2[3-4]。水土资源的流失影响生态安全和经济发展，威胁人民生命财产安全[5]。在2015年8 224起地质灾害中，滑坡和泥石流分别占68%和6%，因灾害造成的直接经济损失达24.9亿元[6]。因此，防治水土流失是确保社会稳定发展的重要前提。林业作为中国生态建设领域的主体能够起到保护生态环境和维持生态平衡的作用，植树造林增加森林覆盖率是防治水土流失的关键。树木根系通过物理和生物作用来提高土体的抗剪强度，同时地表的枯落物层也可以起到增加下渗和削弱洪峰的作用[7]。植被的叶片可以拦截降雨，减少雨水降落时对土壤表面的冲刷和侵蚀，有效提高土壤的稳定性，进而起到保护水土资源的作用[8]。由于地理空间的不同，水土流失区的地形地貌、生态环境和土壤性质等因素不尽相同，这使得水土保持林的树种选择及配置成为保护水土资源的重点和难点。随着地理信息技术以及生态环境与统计模型的发展，生态位模型可以通过最大程度地分离，量化和比较研究区域的气候和空间环境条件，确定环境因子和物种分布之间的稳定关系[9]。根据最大信息熵的状态参数可以估计物种的适宜分布范围，这将为不同地理空间下水土保持树种的选择提供理论依据，同时有利于制定适宜的植物引种和森林经营策略。目前，BIOCLIM，GMPGIS，CLIMEX，GARP和MaxEnt等[10-13]模型在动植物潜在分布区的预测方面应用颇多且具有广泛的前景。在这些模型中，MaxEnt模型是一种在JAVA和最大熵理论的基础上发展起来的用来评估物种潜在适生分布的工具，其可以模拟气候对森林植被潜在分布的影响并合理预测在不同环境条件下物种可能出现的地区[14]。在物种潜在分布的预测方面，其较强的预测能力和相对较高的预测精度是MaxEnt模型相较于其他模型的优点，同时可以支持多种不同的变量类型，灵活性强，结果容易解释，且在样本数量较少的情况下依然可以获得令人满意的预测效果[15]。

板栗树种(Castaneamollissima)属于被子植物门(Angiospermae)，最早记载于《诗经》，栽培历史悠久，以落叶乔木为主，少数为灌木[16]。其具有较好的对温度和水分的适应性，在平均温度8～15 ℃，降水量650～750 mm的地区都适合栽培[17]。板栗树种多分布在沙砾质壤土为主的山地丘陵地区，在土壤贫瘠、坡度较大的边坡也能正常生长。将板栗林建设与水土保持的工程措施相结合，在增加地表覆盖、改良土壤和减少水土流失的同时，还能获得经济效益，提高人民生活水平。为此，本研究采用MaxEnt模型对板栗树种在全国范围内的空间分布格局进行预测并确定其主要影响因子。通过数据库和文献检索确定261个板栗树种分布点和40个环境变量数据，结合ArcGIS平台进行建模和分析，科学地预测当前板栗树种的空间分布现状及其生态特征。对比中国水土流失区的划分情况，为板栗树种在经济林建设中的应用提供理论依据。这有助于板栗树种的适地种植和资源保护，同时也为实现落后地区的水土保持和经济可持续发展提供参考。

1 材料与方法

1.1 地理分布数据

从全球生物多样性数据库(GBIF，http:∥www.gbif.org/)，中国数字植物标本馆(CVH，http:∥www.cvh.org.cn/)，中国在线植物志(eFlora，http:∥www.eflora.cn)和中国知识基础设施工程(CNKI，http:∥www.cnki.net/)获取1950—2017年板栗树种在中国的分布点记录。剔除地理信息缺失和坐标经纬度重复的样本信息，同时核查样本点的植物名录。对于经纬度缺失的有效记录，可以利用经纬度坐标拾取网站(http:∥www.gpsspg.com/)进行查询，确定分布点对应的经纬度信息。另外，为避免分布点的地理自相关，板栗树种分布点在同一个像元内仅保留1个，将获得的261个板栗树种分布点保存为(.csv)文件供后续使用。最后，基于ArcGIS 10.5软件和MaxEnt模型的预测结果来绘制板栗树种的适生区。ArcGIS 10.5是由美国环境系统研究所公司开发的用于数据格式转换和重新分类的综合地理信息系统平台，中国行政区划比例尺矢量图由国家基础地理信息系统提供(http:∥mail.nsdi.gov.cn/)。板栗树种在中国的分布点如图1所示。

注：该图基于中华人民共和国自然资源部标准地图服务网站下载的审图号GS(2019)1823的标准地图，边界无修改。下同。

1.2 环境变量

为预测板栗树种在中国的适生区分布，本研究选取40个环境变量数据，包括19个气候变量，20个土壤变量和1个地形变量。生物气候变量(Bio1—19)是从WorldClim数据库(http:∥www.worldclim.org)下载的，版本为2.1(1970—2000年气候数据，2020年发布)，坐标系为WGS84，空间分辨率为30″(约1 km2)[18]。气候变量反映了温度和降水的特点及其季节性变化特征[19]。土壤变量和地形变量来自世界土壤数据库(http:∥www.fao.org/,HWSDV.1.2)，空间分辨率为30″，地形和土壤变量主要反映了陆地海拔高程和表层土壤(0—30 cm)的理化性质等特征。40个环境变量的多重共线性会影响模型预测的精度和准确性，所以对上述变量的相关性采用Pearson相关分析法进行检验。若两个变量的相关系数大于0.8说明两者具有较强的相关性，应删除两个相关变量中的一个来确保模型模拟的准确性[20]。气候变量、土壤和地形变量中各变量间的相关性如图2所示，根据相关性检验最终确定25个环境变量用于模型建模(表1)。

图2 影响板栗树种分布的40个环境变量(气候、土壤)相关性

表1 用于板栗树种分布模型预测的环境变量

1.3 模型预测及评估

基于板栗树种分布点和环境变量，MaxEnt模型V3.4.1被用于适生区分布的分析和预测。将25个环境变量和261个板栗树种分布点加载到MaxEnt模型中进行建模和分析。参数设置：将数据按照3∶1的比例分配给训练集(training data)和测试集(testing data)，即25%的数据用于测试模型的预测能力[21]。通过R语言对特征组合(FC)和正则化系数(RM)进行优化，其他参数则保留默认值。模型进行500次迭代运算，预测结果以ASCII文件类型输出[22]。以假阳性率(FPR)和真阳性率(TPR)为横、纵坐标的ROC曲线，其曲线面积为AUC值[23]。ROC曲线下的面积是评估模型可靠性的重要方法，通常，AUC值介于0.8～0.9时，认为预测结果很准确，在0.9到1.0范围内时则说明预测结果非常优秀[24-25]。

1.4 适生等级划分和主导环境因子分析

利用ArcGIS 10.5软件对板栗树种预测结果的(.asc)文件进行重分类，参照IPCC报告有关评估可能性的划分方法，p≥0.6为高度适生区，0.4≤p<0.6为中度适生区，0.2≤p<0.4为低度适生区，利用ArcGIS 10.5软件统计不同等级适生区的分布面积。根据Jackknife检验得到每个环境变量的贡献率和置换重要值，以此确定影响板栗树种当前分布的主要环境变量。最后将板栗树种不同等级适生区分布与中国水土流失重点防治区进行比较，确定板栗树种在中国水土流失区的适种范围。

2 结 果

2.1 模型准确性

MaxEnt模型预测结果的准确性通过ROC曲线和坐标轴所包围区域的面积(AUC)来确定。在0～1的范围内，AUC值与模型的预测精度成正相关，即AUC值越高，模型结果的可靠性越高[26]。模型在匹配板栗树种的发生记录和预测适生区方面表现良好，训练集与测试集的AUC平均值分别为0.901和0.869，表明模型对板栗树种适生区的预测达到“极准确”与“很准确”的水平，可以很好地用于模拟和预测板栗树种的空间分布，预测结果可信度高。

2.2 板栗树种潜在地理分布

对MaxEnt模型的预测结果进行重分类，将板栗树种适生区划分为高、中和低三级，其具体分布范围如图3所示。板栗树种潜在适生区面积总计2.92×106km2，约占国土总面积的30.46%。高度适生区面积为4.62×105km2，占其适生区总面积的15.8%，主要分布在四川省和云南省的东北部，湖北省、湖南省和江西省的大部分地区，陕西省、河南省、安徽省、浙江省和山东省的部分地区，尤其在四川省的东北部分布最为集中。中度适生区面积为1.78×106km2，占其适生区总面积的60.9%，主要分布在中国800 mm降水量线以东地区，包括华北平原、长江中下游平原、云贵高原和四川盆地等地区。低度适生区面积为6.82×105km2，占其适生区总面积的23.3%，主要分布在中国400 mm降水量线以东的地区，包括陕西省，山西省、河北省和辽宁省的大部分地区。

图3 基于MaxEnt模型的中国板栗树种潜在适生区分布

2.3 板栗树种在水土流失地区的适用性

水利部《全国水土保持规划(2015—2030年)》中划分的国家级水土流失重点预防和重点治理区如图4所示。其中国家级水土流失重点预防区主要分布在大小兴安岭、呼伦贝尔、长白山、燕山、阴山北麓、武陵山以及大江大河流域等地区。国家级水土流失重点治理区主要包括东北漫川漫岗、大兴安岭东麓、西辽河和大凌河中上游等地区。水土流失区普遍存在于高原山地和大江大河流域，在降雨集中但植被覆盖率不高的山地丘陵地区尤为严重。中国长江流域由于横跨三级阶梯地形复杂，同时受到降雨和植被破坏等因素的影响，水土流失现象也十分突出。根据MaxEnt模型的预测结果，板栗树种高度适生区主要集中在长江流域，整个长江流域范围内均适合板栗树种的生长。在中国水土流失重点预防区中，桐柏山大别山地区、金沙江岷江及三江并流区、丹江口库区、武陵山、新安江流域、湘资沅上游、东江上中游、黄泛平原适合引种板栗树种作为预防水土流失的经济树种。在中国水土流失重点治理区中，沂蒙山泰山地区、西南诸河高山峡谷地区、金沙江下游、嘉陵江和沱江中下游、三峡库区、湘资沅上游、乌江赤水河上中游、滇黔桂岩溶石漠区和粤闽赣红壤区适合引种板栗树种作为水土保持的经济树种。根据板栗树种的潜在适生区分布，青藏地区、西北地区和东北地区不太适宜板栗树种的引种。

图4 国家级水土流失重点预防和重点治理区

2.4 影响板栗树种分布的主导环境变量

在板栗树种适生区分布预测中，Jackknife检验被用于评估不同变量的重要性，进而研究环境因子对板栗树种分布的影响。在仅单一变量被使用的情况下，将25个环境变量按照训练增益由高到低排列，前十位的环境变量依次为年均降水量(Bio12)、年平均温度(Bio1)、最暖季降水量(Bio18)、平均日温差(Bio2)、最干月降水量(Bio14)、温度季节变化方差(Bio4)、最湿季均温(Bio8)、季节降水量变异系数(Bio15)、极端最高温(Bio5)和海拔(Elev)。环境变量按照测试增益由高到低排列，前十位的环境变量依次为年均降水量(Bio12)、年平均温度(Bio1)、最暖季降水量(Bio18)、平均日温差(Bio2)、最干月降水量(Bio14)、温度季节变化方差(Bio4)、最湿季均温(Bio8)、季节降水量变异系数(Bio15)、极端最高温(Bio5)和表层土壤可交换钠盐含量(T_ESP)。贡献率和置换重要值也是评估环境变量影响程度的重要指标。由表2可以看出，年均降水量(Bio12)、年平均温度(Bio1)、表层土壤酸碱度(T_Ph_H2O)、平均日温差(Bio2)和温度季节变化方差(Bio4)这5个环境变量的贡献率位居前列，累积贡献率在83%以上，环境变量的响应曲线见图5。

图5 影响板栗树种分布的主导环境变量的响应曲线

表2 影响板栗树种分布的前10环境变量建模贡献率和置换重要性

3 讨 论

最大熵理论最早提出于1957年，MaxEnt模型便是基于这一理论和JAVA语言发展而来的，目前已成为最常用的物种分布模型[27]。在评估当前气候环境下动植物潜在分布的同时还能够合理预测未来气候变化条件下物种可能出现的适宜分布区，以此判断环境对物种分布区迁移的影响[28-30]。板栗树种的中、高度适生区主要分布在中国800 mm降水量线以东地区，其中高度适生区主要集中在长江流域，在中国的整个南方地区具有重要的应用潜力。对其他物种的预测发现中国西南地区的秦巴山、横断山脉、云贵高原和四川盆地等是漆树的主要适生区，这为漆树的合理栽培提供了依据[31]。赤桉主要分布在东南沿海丘陵、南岭山地和云贵高原西部，未来气候变化将使得赤桉的适生环境进一步扩展[32]。但同样的气候变化背景下，未来胡杨不同等级的潜在适生区面积均有不同程度的缩小，且整体上呈现向高海拔区域迁移的趋势[33]。与此同时，MaxEnt模型还在物种入侵、病虫害防治和其他动物的适生区研究中取得了良好的预测效果[34-36]。值得注意的是，在应用模型预测时由于参数设置缺乏准确的评价指标，参数设置对模型过度拟合的影响尚需进一步研究。同时概率指数模型在预测结果中的应用可能导致预测值超出研究区域气候条件的范围。因此，当外推转移到另一个研究区域或未来(过去)的气候条件时，需要特别注意。

通过对板栗树种潜在适生区分布和环境因子的关系进行研究，结果发现年均降水量(Bio12)、年平均温度(Bio1)、表层土壤酸碱度(T_Ph_H2O)、平均日温差(Bio2)和温度季节变化方差(Bio4)是影响板栗树种适生区分布的最主要环境变量，其累计贡献率在83%以上。年平均温度、年均降水量和温度季节变化方差分别指征热量指标、水分指标以及气候的稳定性。在57%～64%的田间持水量下，板栗树种具有良好的生长状态，其适宜栽培在多雨湿润的气候条件下。在一定的范围内，降雨量与板栗树种产量成正比，但板栗树种在土壤水分过高时其根系容易被泡胀，造成板栗树种的减产甚至枯萎，说明其耐涝性较差[37]。酸碱度pH值为4.6～7.5和含盐量小于0.2%的土壤最适宜板栗树种的栽植，当土壤含盐、碱量较高时会引起渗透压的改变，造成根毛细胞脱水，最后导致植物死亡，在需水怕涝的同时，酸性土壤更加适合板栗树种的生长[38]。热量影响植物的生长速度和光合效率，板栗树种生长环境的温度高低是决定其丰欠的关键因素，因为冬季低温会造成栗树冻害，影响其来年的生长[39]。属于北亚热带季风气候的浙江省中部和东部地区年均温15～17 ℃，板栗树种生长良好，其产量占全省的60%以上[40]。温差有利于植物养分的累积，促使植物更好更快的生长。在一定范围内适当扩大温差能够提高植株内源激素，通过提高植物光合作用增加植物干物质的生成量[41]。相反，过高的温差会抑制植物生长以及营养物质的合成与积累。

水土流失是中国水土资源面临的严重威胁，地形和气候条件的异质性使得全国范围内水土流失的发生程度不同。在水土流失重点预防区，林草覆盖度较大，但依然存在潜在的水土流失风险，对此应采取保护措施来防止水土流失的发生和发展。而重点治理区内由于水土流失严重，需要开展大规模的集中治理。板栗树种作为兼具水土保持生态功能与经济价值为一身的树种，其在水土流失区的应用十分广泛。研究表明坡耕地改种板栗树种后明显减少了年均地表径流量、土壤侵蚀量和养分流失量[42]。其引种到长汀县河田地区后，在防止水土流失的同时实现了丰产和早熟[43]。此外，美国山核桃、香樟、槐树和杨树也都是优秀的兼具水土保持和经济效益的树种，同样具有重要的应用潜力[44]。当前对水土保持树种的研究主要集中在生长状况和生理特征、土壤改良和造林技术等方面，缺乏对不同树种在当前气候条件下的适生分布研究[45-47]。这使得水土流失区的树种选择以经验为主，缺乏必要的理论依据。本研究以板栗树种为研究对象，研究其适生区及主要影响因子，同时结合中国水土流失重点防治区的分布来确定适宜引种板栗树种的地区，为板栗树种在水土流失区的应用和引种栽培提供了理论基础。

4 结 论

(1) 板栗树种潜在适生区面积总计为2.92×106km2，约占国土总面积的30.46%，主要分布在中国400 mm降水量线以东的地区，高度适生区在四川省东北部分布最为集中。受热量指标、水分指标和土壤酸碱度的影响，中国西北地区、青藏地区和东北地区不太适合板栗树种的引种和培育。

(2) 年均降水量、年平均温度、表层土壤酸碱度、平均日温差和温度季节变化方差这5个环境变量对板栗树种适生区分布具有较大的影响，贡献率总和超过83%。板栗树种喜水怕涝，适合在酸性土壤中生长，在一定范围内增加降雨量和温差有利于板栗树种的生长。

(3) 板栗树种的适生区分布基本覆盖中国整个南方地区和大部分的北方地区，对比中国水土流失重点预防区，南方地区和华北地区的水土流失地区都可以考虑选择板栗树种作为主要的经济林树种。