-
气候变化与植物物种地理分布之间的相互关系,不仅是传统的植物生态学和生物地理学的重点研究领域,也是目前全球生态学的研究热点之一。从植物生长环境的角度出发,植物种的地理分布受到土壤、气候、地形等生态因子的影响,其中气候是决定植物种分布的最主要因素,同时植物地理分布区的动态变化也是全球气候变化的最直接反应[1]。政府间气候变化专门委员会(IPCC)第5次评估报告(AR5)指出,自1980年以来,全球平均气温呈上升趋势,至2012年为止,全球平均气温升高(0.85 ± 0.20)℃,并且近60 a升高了0.72 ℃。根据政府间气候变化专门委员会第5次报告的预测,全球平均气温将升高0.30~4.80 ℃[2]。全球气温升高必然改变生态环境,从而对植物种的地理分布格局造成影响。蛛网萼Platycrater arguta为绣球花科Hydrangeaceae蛛网萼属Platycrater的落叶灌木,第三纪孑遗植物,分布于中国浙江南部、江西东北部和福建北部,日本九州、四国等地也有,是东亚植物区系中典型的中国-日本间断分布物种[3],属于国家Ⅱ级保护植物。蛛网萼夏季开花,具二型花,花型优美,具有良好的观赏价值;叶片代为茶饮用,有较好的开发前景。根据多年的野外调查,蛛网萼分布区域极其狭窄,对生境要求极其苛刻,必须要在荫蔽度较高且常年有水分的石壁上才能生长。物种分布模型(SDMs),又称生态位模型,是根据物种的分布信息和相应的环境变量信息,通过特定的算法计算出两者之间的关系,从而将这种关系应用于所研究的地理区域,并对物种可能的生态位进行估计[4-5]。目前,物种分布模型已经被广泛用于外来生物的入侵机制研究[6-7],珍稀濒危植物保护[8-10],生物多样性监测[11-13]等领域。DIVA-GIS是国际马铃薯中心(CIP)开发的一款结合地理信息系统和物种分布数据处理的软件,兼具绘图功能,操作简便,其整合的BIOCLIM模型和DOMAIN模型已经被广泛用于预测物种的潜在分布区[14-15]。本研究利用该软件分析蛛网萼的分布格局,并预测其在气候变化条件下可能的分布区,丰富蛛网萼的生物地理学研究内容,为其保护工作提供依据。
-
蛛网萼实际地理分布点通过查阅标本信息、相关文献和实地调查获得。标本信息通过中国数字植物标本馆(CVH,http://www.cvh.ac.cn/),教学标本资源共享平台(MNH,http://mnh.scu.edu.cn/)和全球生物多样性信息网络(GBIF,https://www.gbif.org/)筛选。同时查阅相关文献[16-17],将每份标本的地理位置尽量精确到乡镇及以下区域,排除重复记录和模糊记录,去除引种记录,利用百度地图拾取坐标系统(http://api.map.baidu.com/lbsapi/getpoint/index.html)析出每份标本的精确坐标,并结合Google Earth确认,最终获得蛛网萼80个精确的有效实际地理分布点。由于用于预测的空间插值气候数据精度为2.5',所以为了消除误差,在做分布区预测时每个2.5'×2.5'的网格只取1条有效记录。
-
气候数据来源于世界气候数据库(WorldClim,http://www.worldclim.org/)。目前的气候数据是根据1950-2000年世界各地的气象站观测信息,利用空间插值法生成的全球气候数据,空间分辨率为2.5';未来气候数据是由美国国家大气研究中心(NCAR)研发的全球大气环流谱模式CCM3模型,在二氧化碳倍增情境下模拟的2100年气候数据[18],空间分辨率为2.5'。数据内包含的19个气候环境变量用作预测物种分布的主要变量,包括年平均气温(BIO1),昼夜温差月平均值(BIO2),等温性(BIO3),温度季节性变化(BIO4),最暖月最高气温(BIO5),最冷月最低气温(BIO6),气温年较差(BIO7),最湿季平均气温(BIO8),最干季平均气温(BIO9),最暖季平均气温(BIO10),最冷季平均气温(BIO11),年降水量(BIO12),最湿月降水量(BIO13),最干月降水量(BIO14),降水量季节性变化(BIO15),最湿季降水量(BIO16),最干季降水量(BIO17),最暖季降水量(BIO18),最冷季降水量(BIO19)。世界各国地图数据来源于DIVA-GIS官方网站(http://www.diva-gis.org/)。
-
将蛛网萼的精确坐标导入DIVA-GIS 7.5.0软件包,1个分布点代表 1个蛛网萼群落,利用Analysis-Point to Grid-Richness功能,统计每个1°×1°网格内出现的蛛网萼群落的数量。利用其内置的Design功能绘制多度分析图。
-
利用DIVA-GIS软件包的Data-Extract Values by Points-From Climate Data功能,提取蛛网萼分布点的19个气候因子数据,将所有80个分布点按纬度分为3个样本较为均匀的梯度,分别为33.50°~36.00°,28.00°~33.50°和26.50°~28.00°,对其气候数据进行方差分析;将气候数据导入PAST 3.16(http://folk.uio.no/ohammer/past/)软件,通过Multivariate-Ordination-Principal components(PCA)功能对其进行主成分分析,得出影响蛛网萼分布的主要气候因子。
-
BIOCLIM模型利用已知的物种分布信息和环境信息,确定物种分布所需的环境变化范围,只要在这个范围的区域就认为是适合物种存在的[19]。DOMAIN模型则是利用Gower算法,生成分布点和分布点之间的相似矩阵,通过矩阵计算目标点上的环境变量的适宜程度,从而确定物种的分布范围[20]。BIOCLIM模型和DOMAIN模型的操作通过DIVA-GIS软件实现。将蛛网萼精确坐标处理成CSV格式,导入DIVA-GIS,然后在Tools-Options选项卡中将气候数据导入,在Modeling-Bioclim/Domain-Precdict选项卡中实现物种分布预测。
-
受试者工作特征曲线(ROC,receiver operating characteristic curve)被广泛用于各种物种分布模型的精度检验[21-22],ROC曲线下的面积即为AUC值(area under curve),AUC值的取值范围为[0.5~1.0],值越大代表模型预测结果越可信,因其不受阈值影响,是目前公认的评价物种分布模型预测结果的最优方法之一[23]。具体操作方法为:在Modeling-Evaluation模块中,将随机取样获得的物种分布点的75%作为训练子集(sample size training data),剩余的25%作为测试子集(sample size test data),训练子集用于训练模型,测试子集用于验证模型,重复3次(表 1)。
表 1 蛛网萼精确分布点
Table 1. Distribution records of Platycrater arguta
序号 分布点 纬度/(°) 经度/(°) 1 中国浙江 30.32 119.43 2 中国浙江 29.48 119.53 3 中国江西 28.63 117.68 4 中国江西 28.58 117.64 5 中国浙江 28.42 121.26 6 中国浙江 28.37 121.11 7 中国浙江 28.37 121.07 8 中国浙江 28.37 121.10 9 中国浙江 28.37 121.10 10 中国浙江 28.36 121.06 11 中国浙江 28.20 118.91 12 中国浙江 28.18 118.91 13 中国江西 28.10 118.12 14 中国浙江 28.06 119.77 15 中国浙江 28.03 121.00 16 中国浙江 28.03 119.51 17 中国福建 28.00 118.17 18 中国江西 27.99 117.04 19 中国浙江 27.94 119.18 20 中国江西 27.94 117.38 21 中国江西 27.94 117.68 22 中国福建 27.94 117.85 23 中国江西 27.93 117.37 24 中国浙江 27.93 119.22 25 中国江西 27.91 117.35 26 中国江西 27.91 117.26 27 中国江西 27.90 117.33 28 中国浙江 27.89 119.19 29 中国浙江 27.85 119.15 30 中国江西 27.80 117.22 31 中国江西 27.80 117.50 32 中国福建 27.76 117.58 33 中国浙江 27.76 119.21 34 中国江西 27.75 117.19 35 中国福建 27.73 117.72 36 中国福建 27.73 117.78 37 中国浙江 27.69 119.14 38 中国福建 27.63 117.90 39 中国福建 27.45 118.02 40 中国浙江 27.05 119.15 41 中国江西 26.51 114.17 42 日本九州 32.82 131.37 43 日本九州 32.57 131.10 44 日本九州 31.43 131.27 45 日本四国 32.72 131.18 46 日本四国 33.55 133.52 47 日本四国 33.58 133.17 48 日本本州 33.72 135.98 49 日本本州 34.98 137.78 50 日本本州 35.10 137.82 51 日本九州 32.11 131.40 52 日本九州 35.43 136.78 53 日本本州 33.65 135.67 54 日本四国 33.67 133.83 55 日本本州 36.00 138.00 56 日本九州 32.10 130.88 57 日本九州 31.72 131.37 58 日本九州 33.20 129.95 59 日本九州 31.74 131.37 60 日本本州 35.65 139.27 61 日本本州 35.67 139.32 62 日本本州 35.02 137.84 63 日本本州 33.63 135.94 64 日本九州 33.35 131.19 65 日本本州 34.41 135.80 66 日本四国 33.56 133.12 67 日本四国 33.52 133.06 68 日本四国 33.17 132.63 69 日本四国 33.53 133.21 70 日本四国 33.80 133.97 71 日本四国 33.73 134.04 72 日本四国 33.48 133.08 73 日本四国 33.65 133.24 74 日本四国 33.53 133.06 75 日本四国 33.81 133.44 76 日本四国 33.51 133.24 77 日本四国 33.74 133.94 78 日本四国 33.70 133.36 79 日本四国 33.45 133.21 80 日本四国 33.81 133.97 -
多度分析结果表明:蛛网萼在中国的分布点主要集中在江西、浙江、福建3个省份,在江西与福建交界的武夷山脉一带分布尤为集中;而在日本,九州、本州、四国均有分布,其中在四国高知一带分布较为集中。结合当地的地理情况,可知蛛网萼的生长环境是需要一定海拔的。
-
由表 2可知:各气候因子之间的差异性不大,可知气候因子对蛛网萼分布的影响与纬度的关系不大;而主成分分析表明,第1主成分的方差贡献率为95%以上,可以代表蛛网萼分布点的主要气候变量信息(表 3)。从表 4可知:第1主成分中年降水量(BIO12),最湿季降水量(BIO16)和最暖季降水量(BIO18)的载荷值较高,可将其定义为“水分因子”。由此可见,水分是对蛛网萼分布影响较为显著的气候限制因子。
表 2 蛛网萼分布的气候指标
Table 2. Climatic indexes of Platycrater arguta
气候指标 33.50°~36.00°N 28.00°~33.50°N 26.50°~28.00°N 平均值±标准差 最小值 最大值 平均值±标准差 最小值 最大值 平均值±标准差 最小值 最大值 BIO1 13.19 ± 2.21 b 8.62 16.78 15.07 ± 2.45 a 8.92 17.51 15.28 ± 1.73 a 11.72 18.25 BIO2 8.15 ± 1.05 a 6.75 10.65 7.82 ± 1.01 6.36 9.63 7.64 ± 0.49 a 6.41 8.73 BIO3 27.14 ± 1.94 a 23.84 29.92 26.8 ± 2.78 a 23.29 32.00 27.25 ± 1.20 a 23.73 29.88 BIO4 795.69 ± 37.84 a 723.60 914.83 780.63 ± 38.96 a 694.42 877.20 739.66 ± 30.86 b 674.42 799.19 BIO5 28.03 ± 2.35 a 22.70 32.10 29.37 ± 2.62 a 22.50 33.20 28.85 ± 2.17 a 24.40 32.50 BIO6 -1.89 ± 2.54 b -8.20 2.80 0.23 ± 2.55 a -5.70 3.70 0.81 ± 1.45 a -2.00 4.70 BIO7 29.92 ± 1.93 a 27.00 35.60 29.14 ± 1.43 a 27.30 32.50 28.04 ± 1.11 b 26.20 30.00 BIO8 21.25 ± 2.25 ab 17.87 25.63 22.28 ± 2.63 a 17.45 26.07 19.83 ± 1.98 b 17.67 27.18 BIO9 4.09 ± 2.66 b -1.68 7.93 7.43 ± 2.71 a 1.25 12.90 8.51 ± 2.85 a 4.68 14.55 BIO10 23.02 ± 2.05 b 18.55 26.23 24.54 ± 2.36 a 18.48 27.33 24.03 ± 1.97 ab 19.90 27.42 BIO11 3.45 ± 2.53 b -1.68 7.93 5.34 ± 2.62 a -1.13 8.97 5.86 ± 1.48 a 2.92 8.97 BIO12 2 207.15 ± 402.63 a 1 420.00 3 254.00 2 159.81 ± 553.20 a 1 453.00 3 424.00 1 979.33 ± 136.29 a 1 611.00 2 131.00 BIO13 325.08 ± 57.21 a 203.00 442.00 355.15 ± 114.36 a 222.00 623.00 326.89 ± 29.50 a 221.00 353.00 BIO14 60.77 ± 15.04 a 43.00 101.00 59.96 ± 15.48 a 38.00 98.00 48.96 ± 3.38 b 40.00 52.00 BIO15 50.19 ± 3.38 b 42.71 53.39 55.34 ± 5.22 a 48.43 64.00 57.18 ± 4.82 a 49.59 66.77 BIO16 859.19 ± 137.88 a 562.00 1 156.00 925.19 ± 301.14 a 568.00 1 622.00 876.33 ± 94.39 a 587.00 958.00 BIO17 221.77 ± 42.88 a 142.00 327.00 216.48 ± 44.17 a 163.00 331.00 186.59 ± 9.36 b 163.00 197.00 BIO18 801.96 ± 134.29 a 550.00 1 137.00 803.89 ± 253.23 a 486.00 1 384.00 686.52 ± 63.26 b 533.00 798.00 BIO19 222.69 ± 42.86 a 142.00 327.00 236.41 ± 42.39 a 185.00 348.00 240.22 ± 19.22 a 184.00 259.00 说明:同行两两之间若有1个相同字母则代表差异性不显著(P>0.05),反之则代表差异显著(P<0.05) 表 3 前3个主成分方差贡献
Table 3. Contribution of variance of first three principal compenent
主成分 特征值 方差贡献率/% 累积方差贡献率/% PC1 238 016.00 95.414 0 95.414 0 PC2 6 891.52 2.762 6 98.176 6 PC3 2 676.09 1.072 8 99.249 4 表 4 第1主成分19个气候因子的得分系数
Table 4. Loadings of first principal compenent relative to 19 environmental variable
气候因子 PC1 BIO1 -0.001 630 0 BIO2 -0.000 110 0 BIO3 0.000 444 0 BIO4 -0.022 710 0 BIO5 -0.002 220 0 BIO6 -0.001 390 0 BIO7 -0.000 830 0 BIO8 -0.000 180 0 BIO9 -0.002 300 0 BIO10 -0.001 760 0 BIO11 -0.001 310 0 BIO12 0.837 050 0 BIO13 0.147 850 0 BIO14 0.023 403 0 BIO15 0.002 295 0 BIO16 0.380 960 0 BIO17 0.072 615 0 BIO18 0.349 890 0 BIO19 0.059 784 0 -
BIOCLIM模型对蛛网萼当下分布区的预测结果表明:目前适宜蛛网萼生存的区域除了已知的武夷山一带、雁荡山一带,在桂林、吉安和黄山一带也有较高的分布可能性。日本的尾鹫一带、长门一带、四国整体以及九州北部均有较高的潜在分布可能性。韩国济州也有极高的分布可能。DOMAIN模型对当下的预测表明蛛网萼的潜在最适区分布在建德和长沙等地,而在日本的潜在最适分布区则与目前已知的分布区域无甚差异。
在二氧化碳倍增的未来气候背景下,BIOCLIM模型对蛛网萼的预测显示其最适分布区逐步向东南方向偏移,主要在云和、丽水一带,并且日本的最适分布区整体向东迁移,集中在四万十市和和歌山;DOMAIN模型则表示未来蛛网萼的分布区与当下相比有一定的缩小。
-
以1-特异度为横坐标,灵敏度为纵坐标绘制ROC曲线,结果如图 1所示。BIOCLIM模型的ROC检验结果显示其AUC值为0.967,DOMAIN模型的AUC值则为0.997,均显著高于随机测试AUC(0.500),说明这2种模型预测的结果都是较为可信的。
-
蛛网萼繁育系统为兼性异交类型与专性异交类型且需要传粉者[24],其繁殖策略是产生大量的小种子用来占据更多适宜的生活位。但张丽芳等[25]的研究发现:蛛网萼种子的发芽条件极其苛刻,导致其发芽率低下。同时我们在野外观察发现蛛网萼的生境要求非常严格,必须在常年有水的石壁附近才能生长,且其种群自我更新绝大部分是以根蘖的形式,极难见到实生苗。这导致了蛛网萼种群遗传多样性的缺失,加之人为等因素使其生境破碎化,蛛网萼种群很难进行大面积的扩张,这也是其种群数量有限的原因之一。生境要求的苛刻,加之自我更新能力有限,推测是导致其如今极其狭窄的分布格局的部分原因。
一般来说,气候因素是决定物种大尺度空间分布的主要因素[26],气候的变化往往会影响物种的分布格局。朱弘等[27]对浙闽樱桃Cerasus schneideriana分布相关的气候因子做主成分分析,得出年降水量和夏季降水量是影响其分布的主要气候因子。王娟等[13]研究了5种锦鸡儿属Caragana的植物分布与15个气候指标之间的关系,得出年降水量、湿润指数等对其中2种锦鸡儿属植物的地理分布有重要影响。以上结果均表明,水分对植物的地理分布有不可忽视的作用。本研究通过差异性分析和主成分分析法分析了与蛛网萼相关的19个气候因子,结果显示:19个气候因子之间的差异并不是很大,表明纬度对于蛛网萼的地理分布影响不大;在足以代表 19个气候变量信息的第1主成分中,年降水量的载荷值最高,结合实地对其生境的调查,可知该气候变量对蛛网萼的地理分布影响较大。可推测蛛网萼分布格局与水分的关系较为密切。
基于当下气候数据的预测结果表明:蛛网萼的分布需要一定海拔,除了已知的分布带以外,中国广西湖南交界、湖南江西交界、湖南西北部以及韩国济州一带都是适生区。目前,蛛网萼的标本信息都较陈旧,缺乏分布信息的更新,而模型预测结果显示江西湖南交界的北部以及湖南广西交界处存在蛛网萼的可能性较大,可作为野外调查的重点区域。由于算法的不同[19-20],2种模型预测出的结果不尽相同,但高度适生区都显示出破碎化,反映出蛛网萼的生存现状已经极其恶劣。BIOCLIM模型中对于蛛网萼中低度适生区的预测集中在赣湘交界和赣浙皖交界区,因为算法的原因,DOMAIN模型对于蛛网萼中低度适生区的预测较为连续,且DOMAIN模型的结果基本覆盖了BIOCLIM模型的结果,两者的预测结果较为接近;2种结果均表明蛛网萼的中低度适生区范围较大,一定程度上反应了这些地区的生境有成为蛛网萼生活区的潜力。陶晓瑜[17]对蛛网萼的群体遗传结构研究表明,蛛网萼的种群之间有一定程度的分化,在保护其遗传多样性的前提下,就地保护的措施是较为妥当的。
在未来气候背景的预测结果中,2种模型均显示蛛网萼的适生区有明显的缩小,且高度适生区有向东部高海拔地区迁移的趋势,这与曹福祥等[28]的报道中提到的全球变暖环境下物种向高海拔地区迁移的结论一致。对于中低度适生区的预测,2种模型均显示有向西迁移的趋势。邵慧等[29]研究表明:BIOCLIM模型用于狭域分布型的物种预测具有较高的模拟准确度;而张琴等[30]对天麻Gastrodia elata的研究表明:对于广布型物种的预测,DOMAIN模型比BIOCLIM模型表现更为出色。尽管两者对不同分布型物种的预测结果不一致,但在本研究中,两者对于蛛网萼的分布区迁移趋势预测基本一致。根据本研究结果,建议对蛛网萼可能的分布区域开展专项调查,以完善其地理分布与资源信息;在做好已有资源保护的同时,加强对浙江东部和福建北部潜在分布地区的生态保护,为其分布扩散提供有利条件。
Modeling the geographic distribution of Platycrater arguta
-
摘要: 为探究国家Ⅱ级野生保护植物蛛网萼Platycrater arguta的地理分布格局,了解它在气候变化下地理迁移的趋势,运用DIVA-GIS软件及其整合的BIOCLIM和DOMAIN模型,结合经过筛选的80个精确分布点信息,对其潜在分布区进行预测并分析影响其分布的主要气候因子。2种模型的预测结果均表明:当下蛛网萼的潜在分布区集中在中国浙赣闽交界处和日本四国岛,且未来气候下蛛网萼的最适生活区有明显的缩小,低度适生区整体向西偏移,并有向高海拔迁移的趋势。受试者工作特征曲线(ROC曲线)检验表明:2种模型的ROC曲线下的面积(AUC值)均达到0.9以上,说明模型预测结果可信。对19个气候变量的主成分分析表明:第1主成分方差贡献率为95%以上,其中年降水量(BIO12),最湿季降水量(BIO16)和最暖季降水量(BIO18)载荷值较高,是影响蛛网萼分布的主要气候因子,表明水分是影响蛛网萼分布的主要因素。建议基于其预测的分布区域开展相应的调查与保护工作。Abstract: To find the geographical distribution pattern of the National Level Ⅱ wild protected plant Platycrater arguta and to understand its geographical migration patterns with climate change, DIVA-GIS software and its integrated BIOCLIM and DOMAIN models were used in combination with 80 precision distribution points. Distribution point information was used to forecast the potential distribution area and analysis included the Receiver Operating Characteristic (ROC) curve test for the main climatic factors affecting distribution and a principal component analysis. The precise distribution information of P. arguta was imported into DIVA-GIS software, and the potential distribution area of P. arguta was predicted through the modelling panel. Prediction results from the two models both indicated that current potential distribution areas of P. arguta were concentrated at the junction of Zhejiang, Jiangxi, and Fujian Provinces in China as well as Shikoku Island in Japan. The optimum living area for P. arguta was greatly reduced with future climate estimates, and the low-grade habitat area shifted overall to the west and to higher altitudes. The ROC curve test showed that the area under the curve (AUC) values of both models reached more than 0.9 indicating that model prediction results were credible. The principal component analysis of 19 climatic variables showed that annual precipitation (BIO12), the most humid seasonal precipitation (BIO16), and the warmest seasonal precipitation (BIO18) were the major climatic factors that affected distribution of Platycrater arguta indicating that water was the main factor. Thus, research and protection strategies should be carried out based on predicted distribution areas.
-
Key words:
- botany /
- Platycrater arguta /
- potential geographic distribution /
- DIVA-GIS /
- BIOCLIM model /
- DOMAIN model
-
林分密度是评定林分内林木间拥挤程度的指标,林分密度不仅影响木材的产量和质量[1],还会影响林分稳定性[2]。在影响林分生长的几个因素中,林分密度是营林工作中能够有效控制的因子。同时,在建立林分生长收获模型时,林分密度也是主要的考虑因素之一。评定林分密度的常用指标中,Reineke的林分密度指数、Beekhuis的相对植距,因为测定容易、与年龄或立地条件关系不密切、有实际经营意义等,近年来被越来越多的人所采用[3-6]。但是,Reineke的林分密度指数、Beekhuis的相对植距存在过度依赖现存数学模型,且对林分要求苛刻。目前,人工神经网络(ANN)在森林生长模型建模方面已有少量应用。谢鹏芳等[7]对马尾松Pinus massoniana直径分布模型进行了研究,认为神经网络建模具有较强的适应性。廖世涛[8]对杉木Cunninghamia lanceolata人工林密度效应进行研究,认为采用人工神经网络建立的密度效应模型科学合理,预测精度优于以往传统模型。马翔宇等[9]以白桦Betula platyphlla人工林为研究对象,建立了单木生长神经网络模型。金星姬等[10]以大兴安岭落叶松Larix gemelinii为研究对象建立了神经网络全林分生长模型。陈晨等[11]对伏牛山栓皮栎Quercus variabilis天然次生林地位指数进行了研究,认为神经网络建模技术具有过度拟合数据的缺点,使用时必须定性定量相结合。由于人工神经网络具有无限逼近非线性系统的能力,且拟合速度快,不依赖于现存的函数[12-13],是目前解决复杂对象建模的最佳手段。林分的生长过程是一个复杂的非线性过程,适合采用人工神经网络来建模。本研究以薄山林场马尾松人工林为研究对象,采用人工神经网络建模技术建立密度指数ANN模型,并与Reineke的林分密度指数进行比较,建立密度指数模型,旨在为提高木材产量和质量,实现马尾松人工林科学抚育经营提供参考。
1. 研究地概况
研究区为河南省驻马店薄山林场,薄山林场为伏牛山东段余脉浅山丘陵区,地理坐标为32°36′~32°44′N,113°50′~113°58′E。山势平缓,坡度为15°~30°,海拔为150~250 m,最高的山峰狠头寨高为430.5 m。土壤为黄棕壤,厚度为30~70 cm,pH 5~6,为弱酸性土壤。该林场位于中国南北气候过渡带,水资源充足,土壤肥沃湿润,物种资源丰富。该地区在气候地理区划中属亚热带向暖温带过渡区,年极端最高气温36.7°,年极端最低气温-11.3°,年均气温为15.5°,年降水多分布在6-8月,年降水量为808~1 206 mm,雨季明显。林场土地总面积为6 017 hm2,林业用地面积5 212 hm2,马尾松人工林和栎类Quercus为该林场的主要树种。
2. 研究方法
2.1 材料收集
于2015年6-8月进行外业调查,并在2016年6-8月做了补充调查。马尾松立地类型有下坡阴坡、下坡阳坡、中坡阴坡、中坡阳坡、上坡阴坡、上坡阳坡等6个。依据立地类型设置不同的标准地。共设147块标准地,标准地面积为20 m × 30 m。标准地主要林分调查因子特征:年龄4~40 a,平均胸径4.2~26.4 cm,平均高3.4~22.7 m。由年龄范围可知,标准地林分年龄组涵盖了幼龄林、中龄林、近熟林、成熟林等。由立地类型、龄组划分、平均高分布等情况可知,所选标准地具有代表性,可以用来建立模型。标准地数据见表 1。
表 1 林分密度统计结果Table 1. Statistical results of stand density径阶/cm 平均胸径/cm 林分密度/(株·hm-2) 疏密度 样地数 4 4.25 9 900 0.9 1 6 6.23 6 285 1.0 10 8 7.85 4 680 1.1 3 10 10.34 3 450 1.0 3 12 12.04 2 520 1.0 16 14 13.87 1 860 1.1 31 16 15.75 1 590 1.0 25 18 17.50 1 275 1.1 25 20 19.36 1 125 1.0 24 22 21.53 1 050 1.1 6 24 23.70 870 0.9 2 26 25.64 600 1.1 1 2.2 数据可行性分析
Reineke的林分密度指数的建立是假设在完满立木度(疏密度为1.0的林分),未经间伐的同龄林中,林分单位面积株数与林分平均胸径遵从幂函数关系。本研究所用标准地的马尾松为未经间伐的人工同龄纯林,林分的郁闭度采用样点法测定,依据郁闭度与疏密度的关系确定疏密度。由表 1可以看出:各组的平均疏密度近似为1.0,因此可认为所用数据符合林分密度指数的完满立木度概念。
2.3 研究方法
2.3.1 林分密度指数法
林分密度指数(ISDI)是Reineke在1933年提出的,指林分在标准平均胸径处应具有的单位面积株数。Reineke所采用标准平均胸径为25.4 cm,中国林业生产中,平均胸径通常取20 cm。本研究采用拐点法确定标准平均胸径D0[13]。N=ISDI(Dg/D0)b;ISDI=N(D0/Dg)b。其中:N为单位面积株数,即林分密度(株·hm-2);Dg为林分平均胸径(cm);D0为标准平均胸径(cm);b为模型参数。
2.3.2 人工神经网络法
在Matlab软件中,应用nntool工具箱,以平均胸径Dg为输入变量,以单位面积株数N为输出变量,建立林分密度指数的ANN模型,中间层为隐层,隐层神经元个数为待定的s个。神经元个数的确定采用由少到多逐步训练的方式确定。图 1中,ws11(s=1, 2, 3…)为输入层节点Dg至隐层神经元的权值,bs1(s=1, 2, 3…)为隐层神经元的阈值,①为输入值为1的输入节点,w1s2(s=1, 2, 3…)为隐层至输出层节点N的权值,b2为输出层阈值,“S”形符号表示神经元的传递函数为对数sigmoid函数,即Matlab软件中的的logsig函数。“/”符号表示神经元的传递函数为线性函数,即Matlab软件中的purelin函数。
3. 结果与分析
3.1 林分密度指数参数求解
依据表 1数据,经回归分析可得幂函数模型参数,N=677 9Dg-1.516 3。最大密度线斜率b=-1.516 3,与Reineke研究的斜率-1.605很接近。图 2可知:在林分发育初期,林分密度随着林分平均胸径的增大而减小,到后期,随着平均胸径的增大,林分密度逐渐趋于平缓或略微下降。平均胸径14 cm处为拐点,因此马尾松人工林标准平均胸径D0为14 cm。除第1个点外,实测值与拟合值拟合效果很好,拟合精度为92.11%。t检验的P常数项为3.78 × 10-14,线性项为7.87 × 10-10,均小于0.05,说明在回归方程中常数项与线性项均是显著的。将b=-1.516 3,D0=14 cm代入林分密度指数计算式:ISDI=N(14/Dg)-1.516 3。
3.2 林分密度指数ANN拟合结果
以平均胸径为输入变量,以单位面积株数为输出变量,依表 1数据为训练样本,对图 1所建ANN模型,在Matlab系统里,经反复训练,同时结合林分密度随林分平均直径的变化规律,即随着林分平均直径的增大,林分密度曲线呈反“S”型,并逐渐趋于平稳,最终确定薄山林场马尾松人工林密度指数模型为1:2:1的ANN模型,模型拟合精度为92.57%,均方误差为0.001 469 7,所得ANN模型共3层,包括输入层1个节点Dg,隐层有2个神经元,输出层有1个节点N。将该模型命名为network2,network2里面包含了所有的模型参数。模型结构为:
$ N=\text { purelin }\left(-0.095 \;403 D_{\mathrm{gl}}-4.706\;9 D_{\mathrm{g} 2}+3.767\;2\right); $
(1) $ D_{\mathrm{g}}=\log \operatorname{sig}\left(7.917\;2 D_{\mathrm{g}}-7.379\;2\right); $
(2) $ D_{\mathrm{g} 2}=\log \operatorname{sig}\left(2.574\;5 D_{\mathrm{g}}+3.267\;6\right)。 $
(3) 式(1)~(3)中:Dg1、Dg2分别为隐层2个神经元的输出值。也可将模型简单写为:
$ N=\operatorname{sim}\left(\text { network } 2, D_{\mathrm{g}}\right)。 $
(4) 拟合效果如图 3,林分密度随平均胸径的增大逐渐减小,达14 cm之后,林分密度逐渐稳定,这与Reineke的林分密度指数曲线变化趋势是一致的。在平均胸径达6 cm之后,ANN能很好地拟合林分密度随林分平均胸径的变化趋势,但在林分发育初期,却未能更好地拟合林分的生长发育过程,这可能与数据的处理有关。若能增加幼龄林组的数据,可能会进一步提升拟合的精度。在式(4)中,当Dg为林分标准平均胸径D0时,公式可写为:
$ I_{\mathrm{SDI}}=\operatorname{sim}\left(\text { network } 2, D_{0}\right)。 $
(5) 由式(4)、式(5)可得林分密度指数表达式:
$ I_{\mathrm{SDI}}=N \operatorname{sim}\left(\text { network } 2, D_{0}\right) / \operatorname{sim}\left(\text { network } 2, D_{\mathrm{g}}\right); $
(6) $ N=I_{\mathrm{SDI}} \operatorname{sim}\left(\text { network } 2, D_{\mathrm{g}}\right) / \operatorname{sim}\left(\text { network } 2, D_{0}\right)。 $
(7) 将D0=14 cm代入式(6)和式(7)可得:
$ I_{\mathrm{SDI}}=125.287\;6 \;{N} / \mathrm{sim}\left(\text { network } 2, D_{\mathrm{g}}\right); $
(8) $ N=0.007\;98 I_{\mathrm{SDI}} \operatorname{sim}\left(\text { network } 2, D_{\mathrm{g}}\right)。 $
(9) 依据式(9)可绘制林分密度指数曲线簇(图 4)。当林分平均胸径相同时,随着林分密度指数的增加,林分密度也在增加。当林分密度相同时,随着林分密度指数的增加,林分的平均胸径在增大。
4. 结论与讨论
本研究数据来自不同立地类型、不同龄级、未经间伐的同龄林标准地,经数据可行性分析,所用标准地具有代表性,运用标准地数据建立的Reineke林分密度指数模型和人工神经网络密度指数模型能够准确反映林分密度的变化情况。采用147块标准地的林分密度—平均胸径数据,以Reineke的林分密度指数曲线为数学模型,拟合林分株数密度与平均胸径的关系曲线。经过回归分析,求得最大密度线斜率b为-1.516 3,模型精度为92.11%,t检验结果显著。用拐点法求得马尾松标准平均胸径为14 cm,林分密度与平均胸径的关系为N=6 779Dg-1.516 3。采用147块标准地的林分密度—平均胸径数据,以林分平均胸径为输入向量,以林分密度为输出向量,构建了人工神经网络密度指数模型,结构为1:2:1。人工神经网络密度指数模型简写为N=sim(network2, Dg),模型拟合精度为92.57%,均方误差为0.001 469 7。模型精度及均方误差检验结果表明:采用标准地数据所建立的人工神经网络林分密度指数模型精度高,能够客观反映该地区林分密度变化情况。
人工神经网络建模技术与Reineke的林分密度指数模型拟合精度均达92%以上,都有较好的拟合效果。与Reineke的林分密度指数曲线相比,人工神经网络密度指数模型具有不依赖于现存数学函数,建模简单、拟合精度高等优势,在幼龄林组数据较少的情况下,人工神经网络具有更好的拟合效果。Reineke林分密度指数模型要求林分为未经间伐的同龄林且平均疏密度近似为1.0,相对于Reineke的林分密度指数,人工神经网络对林分要求不高,适用范围更宽。2种方法在拟合株数密度随林分平均胸径的变化趋势时,幼龄林组拟合效果都不理想,这与幼龄林组数据数量有关,后期研究需增加幼龄林的数量,进一步提高林分密度指数精度。
-
表 1 蛛网萼精确分布点
Table 1. Distribution records of Platycrater arguta
序号 分布点 纬度/(°) 经度/(°) 1 中国浙江 30.32 119.43 2 中国浙江 29.48 119.53 3 中国江西 28.63 117.68 4 中国江西 28.58 117.64 5 中国浙江 28.42 121.26 6 中国浙江 28.37 121.11 7 中国浙江 28.37 121.07 8 中国浙江 28.37 121.10 9 中国浙江 28.37 121.10 10 中国浙江 28.36 121.06 11 中国浙江 28.20 118.91 12 中国浙江 28.18 118.91 13 中国江西 28.10 118.12 14 中国浙江 28.06 119.77 15 中国浙江 28.03 121.00 16 中国浙江 28.03 119.51 17 中国福建 28.00 118.17 18 中国江西 27.99 117.04 19 中国浙江 27.94 119.18 20 中国江西 27.94 117.38 21 中国江西 27.94 117.68 22 中国福建 27.94 117.85 23 中国江西 27.93 117.37 24 中国浙江 27.93 119.22 25 中国江西 27.91 117.35 26 中国江西 27.91 117.26 27 中国江西 27.90 117.33 28 中国浙江 27.89 119.19 29 中国浙江 27.85 119.15 30 中国江西 27.80 117.22 31 中国江西 27.80 117.50 32 中国福建 27.76 117.58 33 中国浙江 27.76 119.21 34 中国江西 27.75 117.19 35 中国福建 27.73 117.72 36 中国福建 27.73 117.78 37 中国浙江 27.69 119.14 38 中国福建 27.63 117.90 39 中国福建 27.45 118.02 40 中国浙江 27.05 119.15 41 中国江西 26.51 114.17 42 日本九州 32.82 131.37 43 日本九州 32.57 131.10 44 日本九州 31.43 131.27 45 日本四国 32.72 131.18 46 日本四国 33.55 133.52 47 日本四国 33.58 133.17 48 日本本州 33.72 135.98 49 日本本州 34.98 137.78 50 日本本州 35.10 137.82 51 日本九州 32.11 131.40 52 日本九州 35.43 136.78 53 日本本州 33.65 135.67 54 日本四国 33.67 133.83 55 日本本州 36.00 138.00 56 日本九州 32.10 130.88 57 日本九州 31.72 131.37 58 日本九州 33.20 129.95 59 日本九州 31.74 131.37 60 日本本州 35.65 139.27 61 日本本州 35.67 139.32 62 日本本州 35.02 137.84 63 日本本州 33.63 135.94 64 日本九州 33.35 131.19 65 日本本州 34.41 135.80 66 日本四国 33.56 133.12 67 日本四国 33.52 133.06 68 日本四国 33.17 132.63 69 日本四国 33.53 133.21 70 日本四国 33.80 133.97 71 日本四国 33.73 134.04 72 日本四国 33.48 133.08 73 日本四国 33.65 133.24 74 日本四国 33.53 133.06 75 日本四国 33.81 133.44 76 日本四国 33.51 133.24 77 日本四国 33.74 133.94 78 日本四国 33.70 133.36 79 日本四国 33.45 133.21 80 日本四国 33.81 133.97 表 2 蛛网萼分布的气候指标
Table 2. Climatic indexes of Platycrater arguta
气候指标 33.50°~36.00°N 28.00°~33.50°N 26.50°~28.00°N 平均值±标准差 最小值 最大值 平均值±标准差 最小值 最大值 平均值±标准差 最小值 最大值 BIO1 13.19 ± 2.21 b 8.62 16.78 15.07 ± 2.45 a 8.92 17.51 15.28 ± 1.73 a 11.72 18.25 BIO2 8.15 ± 1.05 a 6.75 10.65 7.82 ± 1.01 6.36 9.63 7.64 ± 0.49 a 6.41 8.73 BIO3 27.14 ± 1.94 a 23.84 29.92 26.8 ± 2.78 a 23.29 32.00 27.25 ± 1.20 a 23.73 29.88 BIO4 795.69 ± 37.84 a 723.60 914.83 780.63 ± 38.96 a 694.42 877.20 739.66 ± 30.86 b 674.42 799.19 BIO5 28.03 ± 2.35 a 22.70 32.10 29.37 ± 2.62 a 22.50 33.20 28.85 ± 2.17 a 24.40 32.50 BIO6 -1.89 ± 2.54 b -8.20 2.80 0.23 ± 2.55 a -5.70 3.70 0.81 ± 1.45 a -2.00 4.70 BIO7 29.92 ± 1.93 a 27.00 35.60 29.14 ± 1.43 a 27.30 32.50 28.04 ± 1.11 b 26.20 30.00 BIO8 21.25 ± 2.25 ab 17.87 25.63 22.28 ± 2.63 a 17.45 26.07 19.83 ± 1.98 b 17.67 27.18 BIO9 4.09 ± 2.66 b -1.68 7.93 7.43 ± 2.71 a 1.25 12.90 8.51 ± 2.85 a 4.68 14.55 BIO10 23.02 ± 2.05 b 18.55 26.23 24.54 ± 2.36 a 18.48 27.33 24.03 ± 1.97 ab 19.90 27.42 BIO11 3.45 ± 2.53 b -1.68 7.93 5.34 ± 2.62 a -1.13 8.97 5.86 ± 1.48 a 2.92 8.97 BIO12 2 207.15 ± 402.63 a 1 420.00 3 254.00 2 159.81 ± 553.20 a 1 453.00 3 424.00 1 979.33 ± 136.29 a 1 611.00 2 131.00 BIO13 325.08 ± 57.21 a 203.00 442.00 355.15 ± 114.36 a 222.00 623.00 326.89 ± 29.50 a 221.00 353.00 BIO14 60.77 ± 15.04 a 43.00 101.00 59.96 ± 15.48 a 38.00 98.00 48.96 ± 3.38 b 40.00 52.00 BIO15 50.19 ± 3.38 b 42.71 53.39 55.34 ± 5.22 a 48.43 64.00 57.18 ± 4.82 a 49.59 66.77 BIO16 859.19 ± 137.88 a 562.00 1 156.00 925.19 ± 301.14 a 568.00 1 622.00 876.33 ± 94.39 a 587.00 958.00 BIO17 221.77 ± 42.88 a 142.00 327.00 216.48 ± 44.17 a 163.00 331.00 186.59 ± 9.36 b 163.00 197.00 BIO18 801.96 ± 134.29 a 550.00 1 137.00 803.89 ± 253.23 a 486.00 1 384.00 686.52 ± 63.26 b 533.00 798.00 BIO19 222.69 ± 42.86 a 142.00 327.00 236.41 ± 42.39 a 185.00 348.00 240.22 ± 19.22 a 184.00 259.00 说明:同行两两之间若有1个相同字母则代表差异性不显著(P>0.05),反之则代表差异显著(P<0.05) 表 3 前3个主成分方差贡献
Table 3. Contribution of variance of first three principal compenent
主成分 特征值 方差贡献率/% 累积方差贡献率/% PC1 238 016.00 95.414 0 95.414 0 PC2 6 891.52 2.762 6 98.176 6 PC3 2 676.09 1.072 8 99.249 4 表 4 第1主成分19个气候因子的得分系数
Table 4. Loadings of first principal compenent relative to 19 environmental variable
气候因子 PC1 BIO1 -0.001 630 0 BIO2 -0.000 110 0 BIO3 0.000 444 0 BIO4 -0.022 710 0 BIO5 -0.002 220 0 BIO6 -0.001 390 0 BIO7 -0.000 830 0 BIO8 -0.000 180 0 BIO9 -0.002 300 0 BIO10 -0.001 760 0 BIO11 -0.001 310 0 BIO12 0.837 050 0 BIO13 0.147 850 0 BIO14 0.023 403 0 BIO15 0.002 295 0 BIO16 0.380 960 0 BIO17 0.072 615 0 BIO18 0.349 890 0 BIO19 0.059 784 0 -
[1] 倪健, 宋永昌.中国亚热带常绿阔叶林优势种及常见种分布与气候的相关分析[J].植物生态学报, 1997, 21(2):115-129. NI Jian, SONG Yongchang. Relationship between climate and distribution of main species of subtropical evergreen broad-leaved forest in China[J]. Chin J Plant Ecol, 1997, 21(2):115-129. [2] IPCC. IPCC, 2013, Technical summary[G]//STOCKER T, QIN D, PLATTNER G, et al. IPCC, Climate Change 2013: the Physical Science Basis, Contribution of Working Group Ⅰ to the Fifth Assessment Report of the Intergovernmental Panel on Climate Chang. Cambridge: Cambridge University Press, 2013: 710-719. [3] 祁新帅.东亚特有间断分布植物蛛网萼属的生物地理学研究[D].杭州: 浙江大学, 2013. QI Xinshuai. Biogeographic Studies on the Genus Platycrater (Hydrangeaceae) Endemic to East China and South Japan[D]. Hangzhou: Zhejiang University, 2013. [4] 李国庆, 刘长成, 刘玉国, 等.物种分布模型理论研究进展[J].生态学报, 2013, 33(16):4827-4835. LI Guoqing, LIU Changcheng, LIU Yuguo, et al. Advances in theoretical issues of species distribution models[J]. Acta Ecol Sin, 2013, 33(16):4827-4835. [5] 许仲林, 彭焕华, 彭守璋.物种分布模型的发展及评价方法[J].生态学报, 2015, 35(2):557-567. XU Zhonglin, PENG Huanhua, PENG Shouzhang. The development and evaluation of species distribution models[J]. Acta Ecol Sin, 2015, 35(2):557-567. [6] 陆建忠, 翁恩生, 吴晓雯, 等.加拿大一枝黄花在中国的潜在入侵区预测[J].植物分类学报, 2007, 45(5):670-674. LU Jianzhong, WENG Ensheng, WU Xiaowen, et al. Potential distribution of Colidago canadensis in China[J]. Acta Phytotaxon Sin, 2007, 45(5):670-674. [7] 张海涛, 罗渡, 牟希东, 等.应用多个生态位模型预测福寿螺在中国的潜在适生区[J].应用生态学报, 2016, 27(4):1227-1284. ZHANG Haitao, LUO Du, MU Xidong, et al. Predicting the potential suitable distribution area of the apple snail Pomacea canaliculata in China based on multiple ecological niche models[J]. Chin J Appl Ecol, 2016, 27(4):1277-1284. [8] 龚维, 夏青, 陈红锋, 等.珍稀濒危植物伯乐树的潜在适生区预测[J].华南农业大学学报, 2015, 36(4):98-104. GONG Wei, XIA Qing, CHEN Hongfeng, et al. Prediction of potential distributions of Bretschneidera sinensis, an rare and endangered plant species in China[J]. J South China Agric Univ, 2015, 36(4):98-104. [9] 高蓓, 卫海燕, 郭彦龙, 等.应用GIS和最大熵模型分析秦岭冷杉潜在地理分布[J].生态学杂志, 2015, 34(3):843-852. GAO Bei, WEI Haiyan, GUO Yanlong, et al. Using GIS and MaxEnt to analyze the potential distribution of Abies chensiensis[J]. Chin J Ecol, 2015, 34(3):843-852. [10] 马松梅, 聂迎彬, 段霞, 等.蒙古扁桃植物的潜在地理分布及居群保护优先性[J].生态学报, 2015, 35(9):2960-2966. MA Songmei, NIE Yingbin, DUAN Xia, et al. The potential distribution and population protection priority of Amygdalus mongolica[J]. Acta Ecol Sin, 2015, 35(9):2960-2966. [11] 张兴旺, 李垚, 方炎明.麻栎在中国的地理分布及潜在分布区预测[J].西北植物学报, 2014, 34(8):1685-1692. ZHANG Xingwang, LI Yao, FANG Yanming. Geographical distribution and prediction of potential ranges of Quercus acutissima in China[J]. Acta Bot Boreal-Occident Sin, 2014, 34(8):1685-1692. [12] 王刚, 孟艳琼, 孙婧, 等.麻花杜鹃的地理分布模拟[J].西北农林科技大学学报(自然科学版), 2013, 41(5):173-177. WANG Gang, MENG Yanqiong, SUN Jing, et al. Modeling the geographic distribution of Rhododendron maculiferum[J]. J Northwest A & F Univ Nat Sci Ed, 2013, 41(5):173-177. [13] 王娟, 倪健.中国北方温带地区5种锦鸡儿植物的分布模拟[J].植物生态学报, 2009, 33(1):12-24. WANG Juan, NI Jian. Modeling the distribution of five Caragana species in temperate northern China[J]. Chin J Plant Ecol, 2009, 33(1):12-24. [14] 徐晓婷, 杨永, 王利松.白豆杉的地理分布及潜在分布区估计[J].植物生态学报, 2008, 32(5):1134-1145. XU Xiaoting, YANG Yong, WANG Lisong. Geographic distribution and potential distribution estimation of Pseudotaxus chienii[J]. Chin J Plant Ecol, 2008, 32(5):1134-1145. [15] 李垚, 张兴旺, 方炎明.气候变暖对中国栓皮栎地理分布格局影响的预测[J].应用生态学报, 2014, 25(12):3381-3389. LI Yao, ZHANG Xingwang, FANG Yanming. Predicting the impact of global warming on the geographical distribution pattern of Quercus variabilis in China[J]. Chin J Appl Ecol, 2014, 25(12):3381-3389. [16] 乐新贵, 谭策铭, 李林海, 等.贵溪市阳际峰自然保护区新发现成片蛛网萼群落[J].江西林业科技, 2008(4):42-43. LE Xingui, TAN Ceming, LI Linhai, et al. New discovery on flaky community of Platycrater arguta of Yangjifeng Nature Reserve in Guixi City, Jiangxi, China[J]. J Jiangxi For Sci Technol, 2008(4):42-43. [17] 陶晓瑜.东亚特有濒危植物蛛网萼的遗传多样性与亲缘地理学研究[D].杭州: 浙江大学, 2008. TAO Xiaoyu. Genetic Diversity and Phylogeography of Genus Platycrater (Hydrangeaceae) Endemic to East China and South Japan[D]. Hangzhou: Zhejiang University, 2008. [18] GOVINDASAMY B, DUFFY P B, COQUARD J. High-resolution simulations of global climate, part 2:effects of increased greenhouse cases[J]. Clim Dynam, 2003, 21(5/6):391-404. [19] BUSBY J R. BIOCLIM:a bioclimate analysis and prediction system[J]. Plant Protect Quart, 1991, 6(1):64-68. [20] CARPENTER G, GILLISON A N, WINTER J. DOMAIN:a flexible modelling procedure for mapping potential distributions of plants and animals[J]. Biodiv Conserv, 1993, 2(6):667-680. [21] 曹向锋, 钱国良, 胡白石, 等.采用生态位模型预测黄顶菊在中国的潜在适生区[J].应用生态学报, 2010, 21(12):3063-3069. CAO Xiangfeng, QIAN Guoliang, HU Baishi, et al. Prediction of potential suitable distribution area of Flaveria bidentis in China based on niche models[J]. Chin J Appl Ecol, 2010, 21(12):3063-3069. [22] RAES N, ter STEEGE H. A null-model for significance testing of presence-only species distribution models[J]. Ecography, 2007, 30(5):727-736. [23] CHEN Pingfu, WILEY E O, MCNYSET K M. Ecological niche modeling as a predictive tool:silver and bighead carps in North America[J]. Biol Invas, 2007, 9(1):43-51. [24] 张丽芳, 裘利洪.蛛网萼开花物候、花部特征及繁育系统研究[J].广西植物, 2017, 37(10):1301-1311. ZHANG Lifang, QIU Lihong. Flowering phenology, floral traits and breeding system of Platycrater arguta[J]. Guihaia, 2017, 37(10):1301-1311. [25] 张丽芳, 林昌勇, 俞群, 等.珍稀濒危植物蛛网萼的种子形态及萌发特性[J].江西农业大学学报, 2015, 37(3):497-503. ZHANG Lifang, LIN Changyong, YU Qun, et al. The morphological structure and germination characters of seed of rare and endangered species Platycrater arguta[J]. Acta Agric Univ Jiangxi. 2015, 37(3):497-503. [26] 朱耿平, 刘国卿, 卜文俊, 等.生态位模型的基本原理及其在生物多样性保护中的应用[J].生物多样性, 2013, 21(1):90-98. ZHU Gengping, LIU Guoqing, BU Wenjun, et al. Ecological niche modeling and its applications in biodiversity conservation[J]. Biodivers Sci, 2013, 21(1):90-98. [27] 朱弘, 尤禄祥, 李涌福, 等.浙闽樱桃地理分布模拟及气候限制因子分析[J].热带亚热带植物学报, 2017, 25(4):315-322. ZHU Hong, YOU Luxiang, LI Yongfu, et al. Modeling the geographical distribution pattern and climatic limited factors of Cerasus schneideriana[J]. J Trop SubtropBot, 2017, 25(4):315-322. [28] 曹福祥, 徐庆军, 曹受金, 等.全球变暖对物种分布的影响研究进展[J].中南林业科技大学学报, 2008, 28(6):86-89. CAO Fuxiang, XU Qingjun, CAO Shoujin, et al. Advances of global warming impact on sepcies distribution[J]. J Cent South Univ For Technol, 2008, 28(6):86-89. [29] 邵慧, 田佳倩, 郭柯, 等.样本容量和物种特征对BIOCLIM模型模拟物种分布准确度的影响:以12个中国特有落叶栎树种为例[J].植物生态学报, 2009, 33(5):870-877. SHAO Hui, TIAN Jiaqian, GUO Ke, et al. Effects of sample size and species traits on performance of BIOCLIM in predicting geographical distribution of tree species:a case study with 12 deciduous Quercus species indigenous to China[J]. Chin J Plant Ecol, 2009, 33(5):870-877. [30] 张琴, 张东方, 吴明丽, 等.基于生态位模型预测天麻全球潜在适生区[J].植物生态学报, 2017, 41(7):770-778. ZHANG Qin, ZHANG Dongfang, WU Mingli, et al. Predicting the global areas for potential distribution of Gastrodia elata based on ecological niche models[J]. Chin J Plant Ecol, 2017, 41(7):770-778. -
-
链接本文:
https://zlxb.zafu.edu.cn/article/doi/10.11833/j.issn.2095-0756.2019.02.005