摘要：【目的】油松Pinustabulaeformis作为白洋淀大清河流域上游山区主要造林树种，研究其精准适宜性空间分布对流域内水源涵养林的建设与提升具有生产指导作用。【方法】依据植物生态学、植物地理学原理，采用空间限制性因子叠加及最大熵模型(MaxEnt)，以白洋淀大清河流域油松为研究对象，基于野外190个样方数据，结合24个地理环境变量(高程、坡度、坡向、气象、岩石类型、土壤类型等)，其中气象因子通过多元线性回归克里金插值方法(multiplelinearregressionKriging,MLRK)得到，对白洋淀大清河流域油松的适宜性空间分布进行了研究。【结果】油松适宜性预测模型效果达到准确水平(受试者工作曲线ROC=0.869)，积温、年均最大日降水量、高程、年均气温标准差为流域油松适宜性主导地理环境变量，油松高适宜区主要分布在灵丘县、涞源县、涞水县、易县、房山区等县市。【结论】白洋淀大清河流域油松的适宜性条件为：积温2800~3800℃、海拔500~1200m、年均日最大降水量70~110mm、年均气温标准差0.31~0.42℃，按照该适宜性条件进行油松林营造，可充分发挥其有益机能。图8表2参30
关键词：适宜性;多元线性回归克里金(MLRK);空间叠加;最大熵模型;油松
白洋淀及大清河水系是雄安新区水资源安全保障体系的重要部分，但气候干旱、上游土壤流失严重、森林水源涵养功能低下，导致白洋淀干淀、上游断流，难以为新区发展提供保障[1]。土地或岩石的裸露、造林树种的不适宜以及森林空间结构不良等是森林水源涵养功能低下的主要原因[2−3]。早期植树造林时，主要根据区划及经验进行判断，导致油松等部分防护林违背适宜性原则，致使其逐渐衰退为低效林，未达预期效果[4−5]。张春霞等[6]以海拔、坡度、坡向为主要立地因子对陕西黄陵油松Pinustabulaeformis进行立地质量评价，赖文豪等[7]通过地形、土壤、气象等立地因子分析了油松的适宜性情况。可见，已有研究主要集中在小流域等小尺度范围内，难以为当地区域防护林的营造提供支撑，且各立地因子精度较低，如仅用个别站点数据直接代表整个区域，或克里金插值(Kriging)的气象因子站点数据较少，不够精准化，进而导致适宜性分析结果无法贴合生产实际。本研究基于植物生态学特性，考虑了地形、土壤地质、气象等因子，且将各因子落实到30m分辨率，采用了空间限制性因子叠加和最大熵模型(MaxEnt)，对较大尺度下油松的精准适宜性空间分布进行了研究，以期为白洋淀大清河流域造林规划提供科学依据。
1研究区概况与研究方法
1.1研究区概况
大清河系地处海河流域中部(38°10′~40°10′N，113°39′~117°34′E)，西起太行山，东临渤海湾，北临永定河，南至子牙河。河系跨山西、河北、北京、天津4省市，总面积43060km2，其中山区18659km2，丘陵平原24401km2，分别占河系总面积的43.33%和56.67%。大清河流域地处温带半干旱大陆性季风气候区，多年平均气温为6.58~13.26℃，多年平均降水量为365~626mm，年内降水量极为不均，多年平均积温为858~4664℃，多年平均无霜期为73~217d。研究区内优势植物有华北落叶松Larixprincipisrupprechtii、油松Pinustabulaeformis、侧柏Platycladusorientalis、杨树Populusspp.、刺槐Robiniapseudoacacia、臭椿Ailanthusaltissima、绣线菊Spiraeasalicifolia、荆条Vitexnegundovar.heterophylla、胡枝子Lespedezabicolor等。
1.2数据来源
本研究数据分为地形、土壤地质、气象等3类，其中：地形数据来源于地理空间数据云(http://www.gscloud.cn/);土壤地质数据来源于中国寒区旱区科学数据中心(http://westdc.westgis.ac.cn)和全国地质资料馆(http://ngac.org.cn/Map/List);气象数据为站点数据，来源于中国气象数据网(http://data.cma.cn/)各数据集。站点数量为79个，包含流域内36个，流域外43个，各地理环境变量如表1。
在白洋淀大清河流域上游山区，以油松为优势种，随机选取不同高程、坡度、坡向、坡位、土壤类型、龄组的样方进行调查，高程按照中山、低山区间结合太行山植物垂直地带性进行考虑[8]，记录树高、胸径、郁闭度、植被覆盖度、林下灌草种类及盖度和地理位置等信息。在样方调查时，山脉、高大植物会遮挡全球定位系统(GPS)卫星信息接收[9]，故地理位置信息通过手持GPS仪结合卫星地图上比对位置并标注后自动获取，其精确度要高于手持GPS仪直接获取的位置数据，样方分布如图1所示。根据植物学和植物地理学等原理，在对油松可能适宜性空间单元筛选时，考虑人为干预较少的样方和海拔100m以上地区，因此，本研究调查样方主要分布于流域西部山区，共计190个30m×30m油松样方数据(图1)。
相关期刊推荐：《浙江农林大学学报》主要报道林学基础学科、森林培育学、森林经理学、经济林学、林业工程、森林保护学、林木遗传育种学、森林生物学、生态学、生物技术、园林学和园艺学等学科的学术论文、问题讨论和研究简报，适当刊登与农林相关的其他学科的稿件。
在进行气象数据插值时，将高程、坡度、坡向数据标准化后作为辅助变量，利用逐步回归法筛选后进行回归拟合，计算多元线性回归在气象站点位置的残差，再利用普通克里金法对回归残差进行插值，将每个预测点的确定部分和残差插值结果部分相加得到各点的预测值[12]。本研究在插值时考虑到研究区面积较大且跨越不同气候区，故将研究区按照气候区划分为2个区，即半湿润大区与半干旱大区，分别扩宽30km进行插值。为提升流域边缘结果的准确性，站点除选取流域内的36个站点外，另选取流域周围43个站点纳入插值过程，插值完成后再进行拼接得到连续的气象要素空间分布数据，分辨率为30m，根据精度检验，各气象数据MLRK插值结果交叉验证值R2均大于0.9，插值结果达到准确水平。
1.4精准适宜性研究方法
1.4.1油松可能适宜空间单元筛选根据油松的生态习性，利用限制性因子叠加[13]，得到油松生长的可能适宜单元。油松是中国特有树种，主产华北、西北;生于海拔100~2600m，为研究区中、低山区重要的造林树种。喜光，适干冷气候，产区极端低温−25℃，年降水量400~750mm，耐干旱，也能生长在年降水量250mm的部分山区;在酸性、中性、钙质土壤中均能生长，极耐干瘠，石质荒山亦能生长，不耐水涝和盐碱[14−15]。
2结果与分析
2.1最大熵模型精度分析
经MaxEnt软件计算，得到受试者工作曲线(ROC)图，该图反映模型灵敏度，可评价预测模型的准确程度，ROC曲线下面积AUC值越高，评价模型越准确[20]。由图5可见：基于MaxEnt模型和地理环境变量构建的油松适宜性预测模型的训练数据集AUC为0.869，大于0.8，达到很准确水平;检验数据集AUC为0.817，大于0.8，达到很准确水平。综合模型准确性，说明油松适宜性预测模型有效性较高，可用于白洋淀大清河流域油松精准适宜性空间分布研究。
刀切法(Jackknife)[21−22]可预测各地理环境变量对油松适宜性空间分布影响的得分情况，可确定各个环境变量的重要程度。从图6可见：多年年平均积温、多年年平均最大日降水量、高程、多年年平均气温标准差依次是对油松预测贡献率最高的地理环境变量，为流域内油松精准适宜性预测提供了丰富信息。此外，当忽略岩石类型变量时，会最大程度上降低地理环境变量总得分，因此，它具有最多其他变量中没有的信息，也需要综合考虑。
通过MaxEnt模型得到各地理环境变量与发生概率之间的反馈曲线，该曲线可反映单个地理环境变量值与发生概率之间的响应关系。从图7A可以看出：随着多年年平均积温的升高，分布概率呈先增大后减小的趋势，以分布概率0.5为界，积温适宜区间为2800~3800℃。由图7B看出：随着高程升高，分布概率先增大后减小，以分布概率0.5为界，高程适宜区间为500~1200m。由图7C看出：随着多年年平均最大日降水量升高，分布概率先增大后减小，以分布概率0.5为界，年平均最大日降水量适宜区间为70~110mm。由图7D看出：随着多年年平均气温标准差升高，分布概率先增大后减小，以分布概率0.5为界，年平均气温标准差适宜区间为0.31~0.42℃。
2.2适宜性结果分析
参考李培琳[23]、麻亚鸿[24]的适宜性分级标准，确定概率在0~0.05为不适宜区，0.05~0.30为低适宜区，0.30~0.50为中适宜区，0.50~1.00为高适宜区。白洋淀大清河流域油松高适宜区面积为1924km2，中适宜区为3477km2，低适宜区为7621km2，不适宜区为30038km2，中、高适宜区面积达到流域总面积的12.54%，中、高适宜区主要分布在大清河流域上游山区的西北部，在灵丘县、涞源县、涞水县、易县、房山区、阜平县等区县。油松精准适宜性空间分布及其细节放大图如图8，分辨率为30m。表2为油松高适宜区面积大于20km2的县，可见，高适宜区主要集中在涞源县、灵丘县、易县、涞水县、清苑县、涿鹿县，其中涞源县面积最大，为462.60km2，占流域总面积的1%，这是由于这些县符合油松适宜条件的区域较多，在利用油松林进行绿化时可着重考虑。
3讨论
油松是华北地区主要造林树种，由于早期造林时，没有考虑油松适宜性或仅根据气候区划种植了大面积的人工油松林，使得油松生长缓慢，或衰退现象明显[25−26]。近年来，针对油松适宜性的研究仅停留在小尺度，且数据和研究成果精度不高，难以为油松防护林营造提供科学支撑。
本研究从油松生态习性的角度，得出影响油松适宜分布的主导地理环境变量有积温、最大日降水量、高程、年平均气温标准差，贡献率均大于20%，高适宜区分布在积温2800~3800℃、高程500~1200m、年平均最大日降水量70~110mm、年平均气温标准差0.31~0.42℃的区域内。张春霞等[6]认为：油松中低海拔(900~1100m)立地质量为中等以上水平，与本研究高程最大值接近，最小值有差异，是由于其调查样方集中于900m以上。魏安琪[27]认为：坡向是影响油松适宜性最大的因子，但在本研究中坡向不是主导因子，这是由于该研究仅考虑了地形因子未考虑气象因子，且本研究分析后认为坡向通过影响区域局部水热条件而对气象因子产生影响，这在气候因子中有所体现。
在油松林管护过程中，针对不适宜导致的油松低效林，可间伐及间种其他适宜植物，通过逐渐演替由新植物种代替油松成为优势种;其次，在油松防护林营造时，需根据上述条件判断该地块是否适宜油松生长。本研究由于研究区范围较大，仅考虑了土壤类型对油松生长的影响，没有考虑土壤理化性质[28]、土壤厚度[21]等，在防护林建设中需综合考虑。此外，在以地理环境为角度进行油松适宜性研究时，未考虑林业经营管理、经济因素、其他人为因素等。在防护林建设中，整地方式、苗木品质、森林采伐等经营管理技术均对油松有较大影响;从经济角度看，造林成本、管护投入等因素也会造成一定影响[29−30]，也需综合考虑。——论文作者：郭虹扬1,2，史明昌1,2，杨建英1,2，陈春阳1,2
