适地适树是造林的基本原则, 由于过去的数据获取条件和分析技术有限, 人们只能从较大的区域尺度进行立地类型划分, 尤其在气象条件方面区域尺度划分更为粗放, 导致人工林多年以后会出现部分区域林地低效或退化现象。为了解决这一问题, 人们试图通过研究植物生长的复杂环境来科学分析其适宜的区域, 如Mark J.Kimsey等使用地理加权回归模型, 在区域尺度上对美国北部爱德华州地区的冷杉进行立地指数计算, 并进行分析评价, 得到相关结果^[1]。Bourque等利用DEM、温度、降水量、潜在蒸发量、土壤酸碱性等数据, 对牙买加西南海岸的五种生物质能源树种进行了生境适宜性分析及产碳特性分析, 得到了良好的分析结果^[2]。巩垠熙等在旺业甸林场, 以落叶松为研究目标, 利用遥感和地面调查数据, 引入反向传播人工神经网络(Back Propagation Artificial Neural Network, BPANN)模型, 在小班尺度上评价落叶松立地质量^[3]。高蓓、卫海燕等利用GIS和最大熵模型根据气候指标、土壤指标、地形指标来确定100×100m精度上秦岭冷杉的潜在地理分布^[4]。李程利用地理加权回归(Geographically Weighted Regression, GWR)的方法, 结合辽宁红松落界数据, 得到了辽宁省红松适宜性结果和立地指数空间分布格局^[5], 等等。从近年的研究成果来看, 尽管植物适宜性空间单元度逐渐缩小, 对生产决策有一定指导性, 但到具体地块上仍需对适宜性进行精准研判。

人工林低效或退化的原因有多种, 树种不适宜也是其中原因之一。若能精准的分析具体地块上树种适宜程度, 对人工林低效类型和成因分析以及林地改造具有一定的实际应用价值。

本文以湖北省种植面积较大的马尾松为例, 将空间划分为30m的格网, 通过空间拟合、地形分析等模型计算获得每个格网的气象和地形要素, 从地块尺度上排除不适宜马尾松生长的区域, 再结合地面调查单元数据利用最大熵模型计算马尾松最适宜的分布空间, 得到湖北省马尾松适宜性精准空间分布, 为退化人工林生态功能提升改造提供理论支撑和技术参考。

1 研究区概况

湖北省位于108°21′42″E—116°07′50″E、29°01′53″N—33°6′47″N之间, 地势呈三面高, 中间低的状势, 西部整体而言地势较高, 西、北、东部均有山脉环绕, 中南部为平原, 总面积18.59万km²。该研究区位于亚热带季风气候区, 夏季高温多雨, 冬季寒冷干燥。省内林地面积793.89万hm², 占土地总面积的42.71%, 主要分布在山区和丘陵地区, 平原地区分布较少, 林种也以马尾松(Pinus massoniana)、杉木(Cunninghamia lanceolata)、柏树(Platycladus orientalis(L.)Francoptmxjjkmsc)等树种为主, 研究表明, 现有马尾松面积为114.23万hm², 占总林地面积的13.04%, 并且省内大部分马尾松林存在退化现象。

2 数据来源与研究方法 2.1 数据来源

气象数据来源于中国气象数据网站(http://data.cma.cn/site/index.html)的1981—2010年中国地面累年值月值数据, 包含103个气象站点数据(其中湖北省内67个, 相邻省外36个), 并通过整理得到1981—2010年年均温、年均昼夜温差、最热月最高温、最冷月最低温、年温日较差、最小气温日较差、气温标准差、年温范围、等温性这9个气温要素、4个降雨量要素含年降雨量、最湿季降雨量、最干月降雨量、降雨相对标准差和气压、湿度2个要素。

地形数据来源于地理空间数据云(http://www.gscloud.cn)的30m分辨率的数字高程数据(DEM)。利用DEM计算坡度、坡向、高程数据。土壤数据来源于寒区旱区科学数据中心(http://westdc.westgis.ac.cn)的1:100万土壤数据集^[6], 表 1为模型所需的环境数据。

马尾松样点数据主要是按照中山、低山等地形特点在湖北省范围内均匀选取实验区进行调查采样获得, 选取以马尾松为优势种的不同高程、坡度、坡向、坡位、土壤类型、龄组的样地进行调查, 采用每木检尺方法记录树高、胸径、郁闭度、材积量以及经纬度等信息, 分布位置如图 1所示。在实验区随机调查共得到294个20×20m马尾松样方数据, 主要包括本人外业调查的60个样方以及项目组共享的234个样方数据, 内业对数据进行整理, 统一到WGS-84坐标系, 转成CSV格式用于模型计算。

2.2 空间数据处理方法

本文主要采用多元线性回归克里金模型(MLRK)对湖北省的气象空间要素进行精准插值, 得到各气象数据的空间分布图。多元回归克里金插值法是线性回归法和普通克里金插值法的共同模型, 它采用的是最小二乘法(OLS)估计的线性回归, 是优化改进的克里金插值方法, 其预测结果分为两部分, 回归拟合的确定部分；以及普通克里金(OK)对回归残差的插值结果部分。

使用该模型可将每个预测点的确定值和残差部分相加, 从而得到待预测点的预测值, 本文采用此插值方法, 对气温、降雨等气象因子进行插值时, 考虑高程、太阳辐射、坡度、坡向等因素对目标变量的影响, 在R语言环境下处理, 采用逐步回归方法选择关联性较大的辅助变量, 具体形式如下：

式中, PreT表示预测变量, DEM-std表示标准化后的数字高程, slope-std表示标准化后的坡度, aspect变量中包含各个坡向aspectN、aspectNW、aspectW、aspectSW、aspectS、aspectSE、aspectE、aspectNE及平地aspectPT^[7]。

为保证边界处数据精度, 本文在湖北省边界基准外扩30km的基础上, 将15个气象因子插值成30m栅格数据, 再将插值结果中省界外数据裁剪掉得到湖北省的气象空间因子。15个因子差值结果交叉验证值R²均大于90%。图 2为18个空间因子(前三个为地形因子, 后为MLRK模型插值得出的气象因子)：

2.3 空间适宜性分析方法

马尾松属阳性树种, 不耐庇荫, 喜光、喜温。适生于年均温13—22℃, 年降水量800—1800mm, 绝对最低温度到-13℃的区域。根系发达, 主根明显, 有根菌。对土壤要求不严格, 喜微酸性土壤, 但怕水涝, 不耐盐碱, 在石砾土、沙质土、粘土、山脊和阳坡的冲刷薄地上, 以及陡峭的石山岩缝里都能生长。一般在长江下游高程600—700m以下, 中游约1200m以上, 上游约1500m以下均有分布^[8]。本研究依据该树种的生态适宜性, 将各地理因子统一坐标系后分层处理, 分三步分析其适宜性在空间格网上的精准分布。首先利用限制性地理要素进行空间叠加排除不适宜的区域, 然后在可能适宜的区域内利用最大熵模型计算每个格网的适宜指数, 最后按照适宜指数进行适宜性分类。

2.3.1 马尾松可能适宜空间单元筛选

根据马尾松生物生态学习性, 分析出马尾松适宜的生长环境条件, 本文主要是利用限制性因子通过空间叠加, 将不适宜空间单元进行排除。通过对前人研究成果分析, 马尾松限制性因子包括年均温、最冷月最低温、年降水量、土壤类型。高程虽然不是独立的因子, 但在适宜性中有明显的指示作用, 也作为限制因子之一考虑。

SR(Suitable Raster)指马尾松的适宜空间单元, 当SR=1时表示该格网位置理论上适宜马尾松生长, 当SR=0时表示该格网所在位置不适宜。Tem1表示年均温, 当Tem1 < 13℃或Tem1>22℃时, 取“0”, 当13℃≤Tem1≤22℃时, 取“1”；Pre1表示年降水量, 当Pre1 < 800mm或Pre1>1800mm时, 取“0”, 当800mm≤Pre1≤1800mm时, 取“1”；DEM表示高程因子, 当DEM>1500m时, 取“0”, 当DEM≤1500m时, 取“1”；Soil表示土壤要素, 当Soil为水体、沼泽土、草甸土时, 取“0”, 当Soil其它类型土壤时, 取“1”。SR计算结果如图 3所示。

2.3.2 马尾松适宜性空间计算方法

物种分布模型(Species Distribution Models.SDMs), 也称生态位模型^[9], 即使用环境和地理信息解释观测到的物种出现信息关联, 并对目标物种分布进行估计和预测。其中的最大熵(MAXENT)模型^[10], 由于在物种现实生境模拟、生态环境因子筛选等方面作用显著, 能够更好的预测物种的潜在分布区域, 近年来被林业及生态学者广泛应用。本研究在2.3.1计算的基础上利用最大熵(MAXENT)模型, 计算SR=1的格网上马尾松适宜性空间分布, 进一步分析马尾松在不同精准空间上的适宜程度。

最大熵模型的思想是在满足已知条件的情况下, 不对未知情况做任何有偏的假设, 这样得到的模型预测风险最小。也就是概率相等, 这实际上就是说条件熵最大。其公式^[11]为：

式中, X∈{x₁, x₂, …, x_n }, Y∈{y₁, y₂, …, y_n}为离散变量, 在计算中, 模型根据物种“出现点”(样方数据)的环境特征变量运算出目标物种分布的约束条件, 并探寻此约束条件下最大熵的可能分布, 预测目标物种在研究地的生境分布及适宜性^[12]。在树种潜在分布的预测研究中, 同一生态系统中的树种和树种所依附的生境, 通过数理统计, 计算系统具有最大熵时的状态参数确定物种和环境之间的稳定关系, 估计物种的分布^[13]。

3 结果与分析 3.1 MAXENT模型检测与精度验证

最大熵模型提供的刀切法(Jackknife)检验, 能对环境因子贡献率和重要性进行分析, 并利用受试者工作曲线ROC(Receiver operator characteristic curve)评价模型模拟结果的精准度, ROC曲线下的面积AUC(Area Under Curve)越大代表模型模拟效果越好^[14]。

为验证模型的精准度, 本文采用ROC曲线来验证适宜性结果, 如图 4结果显示：训练集(Training data)与测试集(Test data)这两个数据集的AUC值分别为0.826和0.795, 均大于75%, 高于随机预测值, 这表明模型精度较为良好。

3.2 环境变量对适宜性的贡献分析

为获得影响湖北省马尾松生长的关键因素, 如图 5本文采用刀切法(Jackknife)计算了各项因子的贡献率, 结果显示：从贡献百分比来看, 影响湖北省马尾松适生性的主要因子为高程、降雨相对标准差、最湿季降雨量、坡度、年均温, 分别达到42%、19%、10.4%、8.7%和6.3%；气压、年均昼夜温差、气温标准差、年温范围、等温性等因子增益值较小, 贡献率均小于0.30%；在排列重要性上, 高程、年均温、降雨相对标准差排列靠前, 而年均昼夜温差、年温范围、等温性排序靠后。可见, 高程、降雨相对标准差、最湿季降雨量、坡度、年均温为影响湖北省马尾松适宜性分布的主要因子。

年降水量、最冷月最低温对马尾松适宜性分布影响较小的潜在原因可能是其在最大熵计算之前已经将其作为限制性因子进行了不适宜性扣除, 所以在因子贡献分析上表现不突出, 此外, 湖北省处于亚热带季风性气候区, 夏季高温多雨, 雨热同期, 在一定高程范围内, 马尾松可能对于降雨量的需求较大。

3.3 马尾松适生性状况

由最大熵模型(MAXENT)得到精确到地块尺度上的马尾松适宜性分布, 并按照等距法重分类为不适宜区、弱适宜区、中适宜区和强适宜区。如图 6、7所示。

空间分布显示, 湖北全省强适宜区面积为1.28万km², 中适宜区为3.37万km², 弱适宜区为4.53万km², 不适宜区有9.41万km², 中等及以上适生区面积共计4.65万km², 占湖北省总面积的25.01%。其中中、强适宜区主要分布在湖北省东北部及中北部地区, 即黄冈市、孝感市地区、随州市、十堰市的部分地区等；弱适宜区及不适宜区主要分布在中南部及西部地区, 即荆州、仙桃市、潜江市和天门市等地。从县域空间上统计强适宜面积和比率(如表 2)：马尾松的中等及以上适生区主要分布在湖北省东北部及中北部的蕲春县、英山县、罗田县、麻城市、红安县、大悟县、广水市、随州市、浠水县、团风县等地区。其中麻城市适宜性面积最大, 达到0.203万km², 罗田县适宜性区域占比最大, 占全县的83.35%。

从立地条件上看, 土壤因子对马尾松适宜性分布影响较大, 本文2.3.1部分也予以考虑。表 3主要基于地形因子和气象因子对马尾松适宜区进行划分：中等及以上适宜区主要分布在高程低于1200m、坡度在40°以下、年均温13—18℃的区域。

中、强适宜区主要分布在东北部及中北部等地区, 原因是湖北省西部海拔较高, 部分地区甚至高于3000m, 因而相对西部地区, 东部地区更适宜马尾松生长, 其次根据图 2中降雨相对标准差(Pre4)和最湿季降雨量(Pre2)数据显示, 相比于中南部及西部地区, 东部及中北部地区的水分条件较好, 更利于马尾松的生长, 这与马尾松生物生态学习性情况相吻合, 表明在排除限制性因子的基础上采用最大熵模型分析得到马尾松适宜性分布的方法是合理的。

4 讨论

过去调查物种资源适宜性分布的方法一般是通过野外实地调研确定分布范围、数量、种类, 该方法在小尺度区域得到的结果比较准确可靠, 但传统的方法地面调查耗时、耗力、耗费资金, 在地形复杂区耗资更大, 全面的现场调查难度很大。随着地理信息系统、遥感等技术的发展, 将新技术与物种适宜性分析相结合进行评价, 利用最大熵模型(MAXENT)的预测结果使精准分析植物适宜性分布成为可能^[15]。

本研究从马尾松生态适宜性出发来分析其适宜的空间分布, 综合分析得出高程、降雨相对标准差、最湿季降雨量、坡度、年均温对马尾松适宜性分布的影响较大, 分别达到42%、29%、10.4%、8.7%和6.3%, 相关研究也表明, 对马尾松生长影响最显著的气候因子主要是年均温和降雨量, 特别是最湿季降雨量对马尾松的生长有较大的影响^[16-17]；高程、坡度对马尾松的空间分布有一定影响, 杨俊松等对不同坡向30年生的马尾松进行生长指标测定, 结果显示中山区域, 马尾松在阳坡的生长状况要优于阴坡^[18], 而本文坡向的因子贡献率仅占2.7%, 与其产生差异的原因可能是湖北省低山丘陵区占总面积的64%^[19], 加之前期限制因子空间筛选时扣除了大于1500m的山区, 削弱了坡向的影响程度。因此, 在实际应用中, 中山地带坡陡山高地区较多, 阴、阳坡环境差异大, 需要考虑马尾松在阴坡适宜性较差的问题；此外, 土壤也是影响马尾松生长适宜性的重要因素, 土壤类型、土壤质地等直接影响马尾松苗木生长质量^[20-22], 这与本文2.3.1空间单元筛选将土壤要素作为限制因子相吻合。

本文着重从自然环境出发研究马尾松生态适宜性空间分布, 并未扣除城市建筑、水体等不能造林地, 也未考虑林木经营管理和经济因素对马尾松的影响。在经营管理中, 整地方式、种子质量、苗木品质、造林密度、间伐透光等对马尾松生长都会有较大影响^[23-24], 微环境改变也对适宜性有一定影响；从经济角度来看, 交通条件、造林难易程度、管护投入等因素会直接影响造林和经营管理成本, 间接影响马尾松在具体地块上是否适宜；另外, 湖北省建有三峡大坝、葛洲坝、丹江口水库、南水北调渠系等大型水利工程, 水面增加改变了周围微环境, 但影响范围有限。上述情况未在模型中作为因子进行分析, 因此, 在马尾松造林或退化林地改造中, 应结合实际情况全面考虑。