全 文 ：未来气候变化对中国天然橡胶种植
气候适宜区的影响∗
刘少军１，２　 周广胜３∗∗　 房世波３　 张京红２
（ １海南省气象科学研究所， 海口 ５７０２０３； ２海南省南海气象防灾减灾重点实验室， 海口 ５７０２０３； ３中国气象科学研究院， 北京
１０００８１）
摘　 要　 全球气候变暖将严重影响中国天然橡胶种植的气候适宜区分布．根据影响中国橡胶
种植的 ５个主导气候因子，即最冷月平均温度、极端最低温度平均值、月平均温度≥１８ ℃月
份、年平均气温和年平均降水量，基于最大熵 ＭａｘＥｎｔ模型，利用 １９８１—２０１０ 年全国气候数据
和 ＲＣＰ４．５情景的气候预估，分析了 １９８１—２０１０、２０４１—２０６０、２０６１—２０８０ 年中国天然橡胶种
植的气候适宜区变化．结果表明： 随着未来气候变化，２０４１—２０６０ 和 ２０６１—２０８０ 年中国天然
橡胶的种植气候适宜区范围总体呈北扩趋势，对橡胶树北移有利． ２０４１—２０６０、２０６１—２０８０ 年
中国天然橡胶气候适宜区总面积较 １９８１—２０１０ 年呈增长趋势，高适宜区和中适宜区的面积
均有增加趋势，而低适宜区面积呈减少趋势．局部区域气候适宜性发生明显变化：云南的橡胶
主产区的适宜区总面积减少，其中，云南省的景洪、勐腊等地将由现在的高适宜区转变为中适
宜区，海南岛及广东雷州半岛的橡胶种植高适宜区面积明显增加，在台湾岛出现了新的橡胶
种植低适宜区等．
关键词　 气候变化； 天然橡胶； 气候适宜性
文章编号　 １００１－９３３２（２０１５）０７－２０８３－０８　 中图分类号　 Ｓ１６２．５　 文献标识码　 Ａ
Ｅｆｆｅｃｔｓ ｏｆ ｆｕｔｕｒｅ ｃｌｉｍａｔｅ ｃｈａｎｇｅ ｏｎ ｃｌｉｍａｔｉｃ ｓｕｉｔａｂｉｌｉｔｙ ｏｆ ｒｕｂｂｅｒ ｐｌａｎｔａｔｉｏｎ ｉｎ Ｃｈｉｎａ． ＬＩＵ
Ｓｈａｏ⁃ｊｕｎ１，２， ＺＨＯＵ Ｇｕａｎｇ⁃ｓｈｅｎｇ３， ＦＡＮＧ Ｓｈｉ⁃ｂｏ３， ＺＨＡＮＧ Ｊｉｎｇ⁃ｈｏｎｇ２ （１Ｈａｉｎａｎ Ｉｎｓｔｉｔｕｔｅ ｏｆ Ｍｅｔｅｏ⁃
ｒｏｌｏｇｉｃａｌ Ｓｃｉｅｎｃｅ， Ｈａｉｋｏｕ ５７０２０３， Ｃｈｉｎａ； ２Ｈａｉｎａｎ Ｐｒｏｖｉｎｃｅ Ｋｅｙ Ｌａｂｏｒａｔｏｒｙ ｏｆ Ｓｏｕｔｈ Ｃｈｉｎａ Ｓｅａ
Ｍｅｔｅｏｒｏｌｏｇｉｃａｌ Ｄｉｓａｓｔｅｒ Ｐｒｅｖｅｎｔｉｏｎ ａｎｄ Ｍｉｔｉｇａｔｉｏｎ， Ｈａｉｋｏｕ ５７０２０３， Ｃｈｉｎａ； ３Ｃｈｉｎｅｓｅ Ａｃａｄｅｍｙ ｏｆ
Ｍｅｔｅｏｒｏｌｏｇｉｃａｌ Ｓｃｉｅｎｃｅｓ， Ｂｅｉｊｉｎｇ １０００８１， Ｃｈｉｎａ） ． ⁃Ｃｈｉｎ． Ｊ． Ａｐｐｌ． Ｅｃｏｌ．， ２０１５， ２６（７）： ２０８３－２０９０．
Ａｂｓｔｒａｃｔ： Ｇｌｏｂａｌ ｗａｒｍｉｎｇ ｍａｙ ｓｅｒｉｏｕｓｌｙ ａｆｆｅｃｔ ｔｈｅ ｃｌｉｍａｔｉｃ ｓｕｉｔａｂｉｌｉｔｙ ｄｉｓｔｒｉｂｕｔｉｏｎ ｏｆ ｒｕｂｂｅｒ ｐｌａｎｔａ⁃
ｔｉｏｎ ｉｎ Ｃｈｉｎａ． Ｆｉｖｅ ｍａｉｎ ｃｌｉｍａｔｅ ｆａｃｔｏｒｓ ａｆｆｅｃｔｉｎｇ ｒｕｂｂｅｒ ｐｌａｎｔｉｎｇ ｗｅｒｅ ｍｅａｎ ｔｅｍｐｅｒａｔｕｒｅ ｏｆ ｔｈｅ ｃｏｌ⁃
ｄｅｓｔ ｍｏｎｔｈ， ｍｅａｎ ｅｘｔｒｅｍｅｌｙ ｍｉｎｉｍｕｍ ｔｅｍｐｅｒａｔｕｒｅ， ｔｈｅ ｎｕｍｂｅｒ ｏｆ ｍｏｎｔｈｌｙ ｍｅａｎ ｔｅｍｐｅｒａｔｕｒｅ
≥１８ ℃， ａｎｎｕａｌ ｍｅａｎ ｔｅｍｐｅｒａｔｕｒｅ ａｎｄ ａｎｎｕａｌ ｍｅａｎ ｐｒｅｃｉｐｉｔａｔｉｏｎ． Ｃｌｉｍａｔｉｃ ｓｕｉｔａｂｉｌｉｔｙ ａｒｅａｓ ｏｆ ｒｕｂ⁃
ｂｅｒ ｐｌａｎｔａｔｉｏｎ ｉｎ １９８１－２０１０， ２０４１－２０６０， ２０６１－２０８０ ｗｅｒｅ ａｎａｌｙｚｅｄ ｂｙ ｔｈｅ ｍａｘｉｍｕｍ ｅｎｔｒｏｐｙ
ｍｏｄｅｌ ｂａｓｅｄ ｏｎ ｔｈｅ ｆｉｖｅ ｍａｉｎ ｃｌｉｍａｔｅ ｆａｃｔｏｒｓ ａｎｄ ｔｈｅ ｃｌｉｍａｔｅ ｄａｔａ ｏｆ １９８１－２０１０ ａｎｄ ＲＣＰ４．５ ｓｃｅ⁃
ｎａｒｉｏ ｄａｔａ． Ｔｈｅ ｒｅｓｕｌｔｓ ｓｈｏｗｅｄ ｔｈａｔ ｕｎｄｅｒ ｔｈｅ ｂａｃｋｇｒｏｕｎｄ ｏｆ ｔｈｅ ｆｕｔｕｒｅ ｃｌｉｍａｔｅ ｃｈａｎｇｅ， ｔｈｅ ｃｌｉｍａｔｉｃ
ｓｕｉｔａｂｉｌｉｔｙ ａｒｅａ ｏｆ ｒｕｂｂｅｒ ｐｌａｎｔａｔｉｏｎ ｗｏｕｌｄ ｈａｖｅ ａ ｔｒｅｎｄ ｏｆ ｅｘｐａｎｓｉｏｎ ｔｏ ｔｈｅ ｎｏｒｔｈ ｉｎ ２０４１－２０６０，
２０６１－２０８０． Ｔｈｅ ｃｌｉｍａｔｉｃ ｓｕｉｔａｂｉｌｉｔｙ ａｒｅａｓ ｏｆ ｒｕｂｂｅｒ ｐｌａｎｔａｔｉｏｎ ｉｎ ２０４１－２０６０ ａｎｄ ２０６１－２０８０ ｉｎ⁃
ｃｒｅａｓｅｄ ｍｏｒｅ ｏｂｖｉｏｕｓｌｙ ｔｈａｎ ｉｎ １９８１－２０１０． Ｔｈｅ ｓｕｉｔａｂｌｅ ａｒｅａ ａｎｄ ｏｐｔｉｍｕｍ ａｒｅａ ｗｏｕｌｄ ｉｎｃｒｅａｓｅ，
ｗｈｉｌｅ ｔｈｅ ｌｅｓｓ ｓｕｉｔａｂｌｅ ａｒｅａ ｗｏｕｌｄ ｄｅｃｒｅａｓｅ． Ｔｈｅ ｃｌｉｍａｔｉｃ ｓｕｉｔａｂｉｌｉｔｙ ｍｉｇｈｔ ｃｈａｎｇｅ ｉｎ ｓｏｍｅ ａｒｅａｓ，
ｓｕｃｈ ａｓ ｔｈｅ ｔｏｔａｌ ｓｕｉｔａｂｌｅ ａｒｅａ ｗｏｕｌｄ ｄｅｃｒｅａｓｅ ｉｎ Ｙｕｎｎａｎ Ｐｒｏｖｉｎｃｅ， ａｎｄ ｔｈｅ ｓｕｉｔａｂｉｌｉｔｙ ｇｒａｄｅ ｉｎ ｂｏｔｈ
Ｊｉｎｇｈｏｎｇ ａｎｄ Ｍｅｎｇｎａ ｗｏｕｌｄ ｃｈａｎｇｅ ｆｒｏｍ ｏｐｔｉｍｕｍ ａｒｅａ ｔｏ ｓｕｉｔａｂｌｅ ａｒｅａ． Ｈｏｗｅｖｅｒ， ｔｈｅ ｏｐｔｉｍｕｍ ａｒｅａ
ｏｆ ｒｕｂｂｅｒ ｐｌａｎｔａｔｉｏｎ ｗｏｕｌｄ ｉｎｃｒｅａｓｅ ｓｉｇｎｉｆｉｃａｎｔｌｙ ｉｎ Ｈａｉｎａｎ Ｉｓｌａｎｄ ａｎｄ Ｌｅｉｚｈｏｕ Ｐｅｎｉｎｓｕｌａ ｏｆ Ｇｕａｎｇ⁃
ｄｏｎｇ Ｐｒｏｖｉｎｃｅ， ａｎｄ ａ ｎｅｗ ｌｅｓｓ ｓｕｉｔａｂｌｅ ａｒｅａ ｏｆ ｒｕｂｂｅｒ ｐｌａｎｔｉｎｇ ｗｏｕｌｄ ａｐｐｅａｒ ｉｎ Ｔａｉｗａｎ Ｉｓｌａｎｄ ｄｕｅ
ｔｏ ｔｈｅ ｃｌｉｍａｔｅ ｃｈａｎｇｅ．
Ｋｅｙ ｗｏｒｄｓ： ｃｌｉｍａｔｅ ｃｈａｎｇｅ； ｒｕｂｂｅｒ； ｃｌｉｍａｔｉｃ ｓｕｉｔａｂｉｌｉｔｙ．
∗国家自然科学基金项目（４１４６５００５，４１１７５０９６，４１２６５００７）、中国气象局气象关键技术集成与应用项目（ＣＭＡＧＪ２０１３Ｍ３８）和海南省自然科学基
金项目（２０１５４１７２，４０９００５）资助．
∗∗通讯作者． Ｅ⁃ｍａｉｌ： ｇｓｚｈｏｕ＠ ｉｂｃａｓ．ａｃ．ｃｎ
２０１４⁃０８⁃２８收稿，２０１５⁃０３⁃０９接受．
应 用 生 态 学 报　 ２０１５年 ７月　 第 ２６卷　 第 ７期　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　
Ｃｈｉｎｅｓｅ Ｊｏｕｒｎａｌ ｏｆ Ａｐｐｌｉｅｄ Ｅｃｏｌｏｇｙ， Ｊｕｌ． ２０１５， ２６（７）： ２０８３－２０９０
　 　 气候变化已是不争的事实，根据政府间气候变
化专门委员会（ＩＰＣＣ）第五次评估报告，气候变化已
经对海洋和陆地的人类和自然系统产生了广泛影
响［１］ ．气候条件与作物的生产关系密切相关，气候要
素的改变必然导致作物种植格局的变化［２］ ．中国自
１９０４年在云南引种橡胶树以来，至今已有百年橡胶
种植历史，但橡胶产业真正开始蓬勃发展是在 ２０ 世
纪 ５０年代之后．截至 ２０１１年，中国植胶面积已增至
１０７０５．３５ ｋｍ２，年总产值近 ２００ 亿元［３］ ．橡胶树原产
于亚马逊河流域，属典型的热带雨林树种，中国属于
非传统植胶区，天然橡胶种植的气候适宜区域非常
有限．目前中国天然橡胶种植区主要分布在海南、云
南、广东、广西、福建等五省．橡胶对温度、降水、光
照、风速条件均有严格要求．研究表明：在全球气候
变暖背景下，中国热带作物安全种植北界发生了明
显北移，且北移速率呈加快趋势［４－５］ ．因此，与气候条
件密切相关的天然橡胶生产也必将受到未来气候变
化的严重影响．ＩＰＣＣ第五次评估报告（ＡＲ５）给出了
典型浓度路径（ｒｅｐｒｅｓｅｎｔａｔｉｖｅ ｃｏｎｃｅｎｔｒａｔｉｏｎ ｐａｔｈｗａｙｓ，
ＲＣＰｓ）下的气候情景，包括 ＲＣＰ８．５、ＲＣＰ６、ＲＣＰ４．５
及 ＲＣＰ２．６［６－７］ ．其中，ＲＣＰ８．５ 为高端路径，ＲＣＰ６ 和
ＲＣＰ４．５ 为中间稳定路径，ＲＣＰ２．６ 为低端路径，但
ＲＣＰ４．５ 的优先性大于 ＲＣＰ６［８］和 ＲＣＰ２．６［９］ ．目前，
国内外研究者已经开展了大量的基于未来气候情景
的作物影响评价研究，包括作物产量影响评
估［１０－１４］、东北春玉米种植布局［２，１５－１６］、水稻高温热
害［１７－１８］、冬小麦生产潜力［１９－２０］、马铃薯单产趋
势［２１］、柑桔气候适宜性［２２］等．在研究未来气候变化
对物种生境变化的影响时，不同研究者采用的方法
并不相同，包括气候图技术［２３］、农业气候相似距方
法［２４－２７］、模糊综合评判法［２８］、回归模型［２９］、广义线
性模型［３０－３１］、广义加法模型［３２］、分类回归树模
型［３３－３４］、ＣＬＩＭＥＸ 生物种群生长模型［３５－３６］、最大熵
模型［３７－３８］、ＢＩＯＣＬＩＭ 模型［３９－４０］等．不同的研究者也
采用不同的指标组合，包括极端最低温度出现
≤０ ℃概率、阴雨大于 ２０ ｄ内平均温度出现≤１０ ℃
概率、日平均气温≥１５ ℃的活动积温、月平均气温
≥１８ ℃的月份、年平均降水量、年平均风速、≥１０
级风出现的概率等，开展了天然橡胶种植的气候适
宜性区划研究［４１－４４］ ．但迄今关于未来气候变化情景
下中国天然橡胶种植的气候适宜区变化研究还未见
报道，制约着气候变化背景下中国天然橡胶的区域
布局和科学规划决策的制定．
本研究试图利用全球气候模式（ＢＣＣ⁃ＣＳＭ１⁃１）
给出的 ＲＣＰ４． ５［４５］ 排放情景下 ２０４１—２０６０ 和
２０６１—２０８０ 年的气候预估资料以及基准时段
（１９８１—２０１０ 年）的气候资料，结合最大熵模型
（ＭａｘＥｎｔ），构建中国橡胶种植分布与气候因子的关
系模型，确定影响中国橡胶种植的主导气候因子，探
讨中国天然橡胶种植的气候适宜性及其对未来气候
变化的响应，以期为中国天然橡胶应对气候变化和
种植布局调整提供决策依据．
１　 数据来源与研究方法
１􀆰 １　 数据来源
１９８１—２０１０年气候标准年的中国地面气象逐
日数据来源于中国气象局国家气象信息中心，包括
温度、降水、风速、辐射等要素．１９８１—２０１０年影响中
国的台风路径数据源于中国台风网（ｗｗｗ． ｔｙｐｈｏｏｎ．
ｇｏｖ．ｃｎ） “ ＣＭＡ⁃ＳＴＩ 热带气旋最佳路径数据集”．
ＩＰＣＣ第五次评估报告 （ＡＲ５）基于全球气候模式
（ＢＣＣ⁃ＣＳＭ１⁃１） 给出的 ＲＣＰ４． ５ 气候情景数据
（２０４１—２０６０、２０６１—２０８０ 年）来源于 ｈｔｔｐ： ／ ／ ｗｏｒｌｄ⁃
ｃｌｉｍ．ｏｒｇ ／ ｃｍｉｐ５．中国橡胶林地理分布数据（１９５４—
２００９年）来自中国数字植物标本馆、国内外公开发
表的相关论文和实地考查数据等，共 １０６ 条．中国国
界、省界和县界行政区划图来源于国家基础地理信
息网站提供的 １ ∶ ４００ 万基础地理信息数据
（ｈｔｔｐ： ／ ／ ｎｇｃｃ．ｓｂｓｍ．ｇｏｖ．ｃｎ ／ ） ．
１􀆰 ２　 研究方法
最大熵原理是在已知部分知识前提下关于未知
分布最合理的推断，即符合已知知识最不确定或最
随机的推断［４６］ ．ＭａｘＥｎｔ 模型是一种基于最大熵原理
的物种分布模型，主要应用于物种的潜在生境预
测［４７］ ．具体算法和原理参考文献［３７－３８］．天然橡胶
种植的适宜区预测采用最大熵 ＭａｘＥｎｔ 模型 ３．３．３ｋ
版实现 （ ｈｔｔｐ： ／ ／ ｗｗｗ． ｃｓ． ｐｒｉｎｃｅｔｏｎ． ｅｄｕ ／ ～ ｓｃｈａｐｉｒｅ ／
ｍａｘｅｎｔ ／ ），该模型在预测物种分布方面具有很强的
优势［４８－４９］，具体步骤如下：
１）评价因子的选择．基于已有的天然橡胶种植
区划［４１－４４］，结合天然橡胶生产区域的气候特征，筛
选出 ９个可能影响橡胶种植分布的潜在气候因子：
年平均降水量、最冷月平均温度、最暖月平均温度、
极端最低温度平均值、年辐射量、年平均温度、月平
均气温≥１８ ℃的月份、台风影响该区域概率、年平
均风速．
２）数据格式的转换．根据 ＭａｘＥｎｔ模型运行的需
求，将年平均降水量、最冷月平均温度、最暖月平均
４８０２ 应　 用　 生　 态　 学　 报　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 ２６卷
温度、极端最低温度平均值、年辐射量、年平均温度、
月平均气温≥１８ ℃的月份、台风影响该区域概率、
年平均风速等 ９ 个因子转换为 ＡＳＣＩＩ 文件，坐标系
为 ＷＧＳ⁃８４，作为环境变量输入到最大熵模型；将
１０６个橡胶种植分布信息点数据按经度和纬度顺序
储存成 ｃｓｖ 格式的文件，作为训练样本输入到最大
熵模型．
３）主导因子的确定．基于最大熵模型的 Ｊａｃｋ⁃
ｋｎｉｆｅ模块评价可以得到，９ 个潜在气候因子的重要
性排序为：最冷月平均温度＞极端最低温度平均值＞
月平均温度≥１８ ℃月份＞年平均气温＞年平均降水
量＞台风影响概率＞最暖月平均温度＞年辐射量＞年
平均风速．９ 个可能影响气候因子对橡胶林存在的
贡献率分别为：最冷月平均温度 １５．９％，极端最低温
度平均值 ６９．２％，月平均温度≥１８ ℃月份 ２．９％，年
平均气温 １．０％，年平均降水量 ４．４％，台风影响概率
０．３％，最暖月平均温度 １．７％，年辐射量 ０．５％，年平
均风速 ４．０％．根据最大熵模型计算出 ９ 个潜在气候
因子对橡胶种植区分布的重要性排序和贡献率大
小，确定影响橡胶种植的主导气候因子为最冷月平
均温度、极端最低温度平均值、月平均温度≥１８ ℃
月份、年平均气温、年平均降水量，５ 个因子的累积
贡献率为 ９３．５％，具体方法参考文献［５０］．
４）气候适宜性分区．将最冷月平均温度、极端最
低温度平均值、月平均温度≥１８ ℃月份、年平均气
温、年平均降水量等 ５个因子作为环境因子，输入最
大熵模型，最大熵模型预测结果即是天然橡胶在待
预测地区的存在概率．根据概率的大小，划分中国橡
胶种植的气候适宜区．
本研究的适宜性分区标准结合不同气候保证率
来划分，取概率＜５０％保证率为不适宜区、５０％～８０％
保证率为低适宜区、８０％ ～ ９０％保证率为中适宜区、
≥９０％保证率为高适宜区．考虑到气候资源不同保
证率（ａ）以及影响中国橡胶林种植分布的 ５ 个主导
气候因子，则某地可安全种植橡胶树的概率拟为适
宜气候条件下的 ａ５，因此，０．５５ ＝ ０．０３，０􀆰 ８５ ＝ ０．３２，
０􀆰 ９５ ＝ ０．６０．根据概率的大小，取 Ｐ ＝ ０􀆰 ０３ 为分界线，
确定天然橡胶的气候适宜区的范围，并建立了气候
适宜性分区标准：Ｐ≥０．０３为天然橡胶可种植区，Ｐ＜
０．０３ 为天然橡胶种植不适宜区，０．０３≤Ｐ＜０．３２ 为低
适宜区，０．３２≤Ｐ＜０．６０ 为中适宜区，Ｐ≥０．６０ 为高适
宜区．
１􀆰 ３　 评价结果精度分析
最大熵模型的精度检验采用受试者工作特征曲
线（ｒｅｃｅｉｖｅｒ ｏｐｅｒａｔｉｎｇ ｃｈａｒａｃｔｅｒｉｓｔｉｃ ｃｕｒｖｅ，ＲＯＣ）与横
坐标围成的面积（即 ＡＵＣ 值）来评价模型预测结果
的精准度，ＡＵＣ的大小作为模型预测准确度的衡量
指标，其取值范围为［０，１］，值越大表示模型判断力
越强［５１］ ．ＡＵＣ值取 ０．５ ～ ０．６ 为失败，０．６ ～ ０．７ 为较
差，０．７～０．８ 为一般，０．８ ～ ０．９ 为好，０．９ ～ １．０ 为非常
好［５２］ ．
１􀆰 ４　 模型适用性检验
为检验最大熵模型在预测天然橡胶种植分布区
的适用性，以 １９８１—２０１０ 年的气候数据为基础，随
机选取 ７５％的天然橡胶分布点数据用于构建模型，
剩下 ２５％的天然橡胶分布点用于模型的验证．通过
最大熵模型和 ５个可能的气候因子构建天然橡胶种
植分布⁃气候关系模型，模型运算结果的训练集和验
证集 ＡＵＣ 值分别为 ０．９９４ 和 ０．９８９，表明所构建模
型的预测精度达到“非常好”的标准，可以用于预测
天然橡胶种植区范围．
２　 结果与分析
２􀆰 １　 气候变化背景下橡胶种植北界演变特征
将 １９８１—２０１０、２０４１—２０６０、２０６１—２０８０ 年的
数据分别代入 ＭａｘＥｎｔ模型，计算橡胶在待预测地区
的存在概率．模型运行结果显示，３ 个时间段数据计
算的 ＡＵＣ 值分别为 ０．９９３、０．９９１、０．９９１，表明 Ｍａｘ⁃
Ｅｎｔ模型可以用于橡胶种植空间分布的预测．
根据最大熵模型给出的存在概率，分别取
５０％、８０％气候资源保证率下的概率作为分界．在考
虑 ５０％气候保证率条件下，中国天然橡胶种植的北
界在 １９８１—２０１０、２０４１—２０６０、２０６１—２０８０ 年的趋
势基本一致，其中，橡胶种植的北界在 ２０４１—２０６０、
２０６１—２０８０年较 １９８１—２０１０ 年略有北扩，同时台
湾岛出现橡胶种植气候适宜区（注：台湾岛出现的
适宜区界线是环形状，并非单纯北界）（图 １ａ）．在考
虑 ８０％气候保证率条件下，中国天然橡胶种植的北
界发生了明显变化．在云南省范围内，橡胶种植的北
界并未出现北扩趋势，而是向南逐渐缩小适宜区域，
１９８１—２０１０ 年橡胶种植的北界明显大于 ２０４１—
２０６０、２０６１—２０８０年橡胶种植北界．而在广西、广东、
福建，橡胶种植的北界保持北扩趋势，１９８１—２０１０
年橡胶种植的北界明显小于 ２０４１—２０６０、２０６１—
２０８０ 年橡胶种植北界，２０４１—２０６０ 与 ２０６１—２０８０
年橡胶种植北界基本一致（图 １ｂ）．
２􀆰 ２　 天然橡胶种植气候适宜区的变化
根据天然橡胶种植气候适宜区的划分标准，得
５８０２７期　 　 　 　 　 　 　 　 　 　 　 刘少军等： 未来气候变化对中国天然橡胶种植气候适宜区的影响　 　 　 　 　 　
图 １　 ５０％（ａ）和 ８０％（ｂ）气候保证率预测橡胶种植适宜区
北界
Ｆｉｇ．１　 Ｎｏｒｔｈｅｒｎ ｂｏｕｎｄａｒｉｅｓ ｏｆ ｒｕｂｂｅｒ ｐｌａｎｔａｔｉｏｎ ｗｉｔｈ ５０％ （ａ）
ａｎｄ ８０％ （ｂ） ｇｕａｒａｎｔｅｅ ｒａｔｅ ｉｎ Ｃｈｉｎａ．
到 １９８１—２０１０、２０４１—２０６０、２０６１—２０８０ 年中国天
然橡胶种植气候适宜区的分布图．从图 ２ 可以看出，
１９８１—２０１０年中国天然橡胶种植气候适宜区分布
在海南、广东、福建、广西、云南 ５ 省；低适宜区主要
分布在云南盈江—潞西—永德—双江—景谷—元
江—元阳—河口以南区域，广西低适宜区分布在靖
西—田东—马山—来宾—藤县—岑溪以南区域，广
东的低适宜区分布在德庆—四会—三水—惠东—海
丰—丰顺等以南区域，福建的低适宜区分布位于平
和—漳州—福清以南；中适宜区主要分布在海南的
北部、广东的雷州半岛、广西的南部；气候高适宜区
主要分布在海南省的儋州、乐东、白沙、保亭，云南省
的景洪、勐腊等地；低适宜区以北的其他地区为天然
橡胶种植区的气候不适宜区．
２０４１—２０６０年中国天然橡胶种植气候适宜区
分布较 １９８１—２０１０年的计算结果发生了明显变化，
其中，天然橡胶种植的低适宜区除分布在广东、福
建、广西、云南外，在台湾岛的西部和东部的沿海区
域也出现了天然橡胶种植的低适宜区；气候高适宜
区只出现在海南岛、广东省的雷州、徐闻等地；云南
种植橡胶的适宜区面积整体减少，而且云南的景洪
和勐腊等主产区，由高适宜区变为了中适宜区．
２０６１—２０８１ 年中国天然橡胶种植适宜区分布与
２０４１—２０６０年中国天然橡胶种植适宜区分布趋势
基本一致，但天然橡胶种植的高适宜区较２０４１—
图 ２　 １９８１—２０１０（ａ）、２０４１—２０６０（ｂ）和 ２０６１—２０８０（ｃ）年
橡胶种植气候适宜区
Ｆｉｇ． ２ 　 Ｃｌｉｍａｔｉｃ ｓｕｉｔａｂｉｌｉｔｙ ａｒｅａ ｏｆ ｒｕｂｂｅｒ ｐｌａｎｔａｔｉｏｎ ｉｎ Ｃｈｉｎａ
ｄｕｒｉｎｇ １９８１－２０１０ （ａ）， ２０４１－２０６０ （ｂ） ａｎｄ ２０６１－２０８０ （ｃ）．
Ⅰ： 不适宜区 Ｕｎｓｕｉｔａｂｌｅ ａｒｅａ； Ⅱ： 低适宜区 Ｌｅｓｓ ｓｕｉｔａｂｌｅ ａｒｅａ； Ⅲ：
中适宜区 Ｍｉｄｄｌｅ ｓｕｉｔａｂｌｅ ａｒｅａ； Ⅳ： 高适宜区 Ｏｐｔｉｍｕｍ ａｒｅａ．
２０６０年预测结果略有减少，而且在分布上整体向东
南偏移．
从面积上看，２０４１—２０６０、２０６１—２０８０ 年中国
天然橡胶气候适宜区总面积较 １９８１—２０１０ 年呈增
长趋势（表 １）．１９８１—２０１０年，中国天然橡胶的气候
适宜性种植总面积为 ３９．６１×１０４ ｋｍ２，占全国陆地面
积的 ４．１％，其中，高适宜区面积为 ２．９６×１０４ ｋｍ２，中
适宜区面积为 １０． ０９ × １０４ ｋｍ２，低适宜区面积为
２６􀆰 ５６×１０４ ｋｍ２ ．２０４１—２０６０ 年，中国天然橡胶的适
宜性种植总面积为 ４５．５１×１０４ ｋｍ２，占全国陆地面积
的 ４􀆰 ７％，其中，高适宜区面积为 ４．０４×１０４ ｋｍ２，中适
宜区面积为 １８．８７×１０４ ｋｍ２，低适宜区面积为 ２２．０６×
１０４ ｋｍ２ ．２０６１—２０８０ 年，中国天然橡胶的适宜性种
植总面积为 ４６． ０２ × １０４ ｋｍ２，占全国陆地面积的
４􀆰 ８％，其中，高适宜区面积为 ３．９０×１０４ ｋｍ２，中适宜
区面积为 １７．７３×１０４ ｋｍ２，低适宜区面积为 ２４．３９×
１０４ ｋｍ２ ．
６８０２ 应　 用　 生　 态　 学　 报　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 ２６卷
表 １　 不同时期天然橡胶种植气候适宜区面积及比例
Ｔａｂｌｅ １　 Ｐｏｔｅｎｔｉａｌ ｄｉｓｔｒｉｂｕｔｉｏｎ ｏｆ ｃｌｉｍａｔｉｃ ｓｕｉｔａｂｌｅ ａｒｅａｓ ａｎｄ ｒａｔｉｏ ｆｏｒ ｒｕｂｂｅｒ ｐｌａｎｔａｔｉｏｎ ｉｎ ｄｉｆｆｅｒｅｎｔ ｐｅｒｉｏｄｓ
适宜性分级
Ｇｒａｄｅ ｏｆ ｃｌｉｍａｔｉｃ
ｓｕｉｔａｂｌｉｌｉｔｙ
１９８１—２０１０
面积
Ａｒｅａ （１０４ ｋｍ２）
比例
Ｒａｔｉｏ （％）
２０４１—２０６０
面积
Ａｒｅａ （１０４ ｋｍ２）
比例
Ｒａｔｉｏ （％）
２０６１—２０８０
面积
Ａｒｅａ （１０４ ｋｍ２）
比例
Ｒａｔｉｏ （％）
不适宜区 Ｕｎｓｕｉｔａｂｌｅ ａｒｅａ ９１９．９４ ９５．９ ９１４．０４ ９５．３ ９１３．５３ ９５．２
低适宜区 Ｌｅｓｓ ｓｕｉｔａｂｌｅ ａｒｅａ ２６．５６ ２．８ ２２．６０ ２．４ ２４．３９ ２．５
中适宜区 Ｍｉｄｄｌｅ ｓｕｉｔａｂｌｅ ａｒｅａ １０．０９ １．１ １８．８７ ２．０ １７．７３ １．９
高适宜区 Ｏｐｔｉｍｕｍ ａｒｅａ ２．９６ ０．３ ４．０４ ０．４ ３．９０ ０．４
３　 讨　 　 论
针对气候变化背景下中国橡胶生产布局及其应
对气候变化政策制定的需求，基于已有研究成果从
年尺度筛选出的影响橡胶种植区分布的潜在气候因
子，结合橡胶种植地理分布信息，利用最大熵模型和
ＡｒｃＧＩＳ 空间分析技术，以及 １９８１—２０１０ 年气候数
据和 ＲＣＰ４． ５ 情景的气候预估数据 （２０４１—２０６０、
２０６１—２０８０年），研究种植区橡胶林分布的气候适
宜性的变化规律．
ＲＣＰ４．５ 情景下，２０４１—２０６０、２０６１—２０８０ 年中
国天然橡胶的总气候适宜区面积（低适宜、中适宜、
高适宜）较 １９８１—２０１０ 年种植适宜面积有所增加，
界限有北扩趋势，主要原因在于 ＲＣＰ４． ５ 情景下
２０４１—２０６０、２０６１—２０８０ 年的气候预估与 １９８１—
２０１０年数据相比，年平均气温和年平均降水量在橡
胶种植区域内整体呈增加趋势，对橡胶树北移有利，
但最冷月平均温度、极端最低温度平均值、月平均温
度≥１８ ℃月份呈减小趋势，对橡胶树北移不利．最
大熵模型充分考虑了各影响因子的相互作用，能够
更加客观地反映橡胶种植的潜在分布．由于未来气
候变化的影响，未来中国橡胶种植的气候适宜区有
出现不适宜的可能，如云南的橡胶主产区面临气候
变化的负影响，适宜区总面积将减少，其中，云南省
的景洪、勐腊等地将由现在的高适宜区转变为中适
宜区．导致云南橡胶适宜区减少的主要因素为最冷
月平均温度、极端最低温度平均值、月平均温度
≥１８ ℃月份呈减小趋势，其影响程度大于年平均气
温和年平均降水量增加所产生的作用．但一些地区
的橡胶种植气候适宜性将增加，如海南岛及广东雷
州半岛的橡胶种植高适宜区面积明显增加，主要原
因在于 ＲＣＰ４．５ 情景下预估的最冷月平均温度、极
端最低温度平均值、月平均温度≥１８ ℃月份、年平
均气温和年平均降水量等因子在空间上存在差异，
导致低适宜区、中适宜区、高适宜区在空间分布上与
１９８１—２０１０年中国橡胶气候适宜区相比发生了较
大变化．
总体而言，在未来气候情景下，中国天然橡胶种
植适宜区总面积总体呈增加趋势，高适宜区和中适
宜区的面积均有增加趋势，而低适宜区面积呈减少
趋势．中国天然橡胶种植气候适宜区非常有限，目前
仅海南、云南、广西、广东、福建等 ５省存在种植气候
适宜区．根据预测结果，中国橡胶种植气候适宜区的
总面积在 ３９．６１×１０４ ～４６．０２×１０４ ｋｍ２，占全国陆地面
积的 ４．２％～４．８％．橡胶是国民经济建设不可缺少的
重要战略物资，因此利用预测的结果，在有限的橡胶
种植气候适宜区开展天然橡胶的布局和规划具有重
要的现实意义．
根据中国农林气候区划协作组提出的中国橡胶
种植区［４２］、农业部热带作物区划办公室编制的中国
橡胶种植适宜区［５３－５４］、王利溥［４３］、王菱［４１］的研究成
果及目前中国橡胶主产区实际种植范围［海南省植
胶现 状 图 （ ４９００ ｋｍ２ ）、 云 南 省 植 胶 现 状 图
（４９１３ ｋｍ２）、广东省植胶现状图（４１３ ｋｍ２）］ ［５３］，对
比分析 １９８１—２０１０ 年中国天然橡胶的种植区及等
级划分可以发现，利用最大熵模型获取的结果基本
符合当前中国天然橡胶的种植现状，能客观反映目
前橡胶种植的实际情况，从一定程度上说明该方法
是可行的．
气候变化对不同区域的影响存在差异．虽然气
候变化对中国天然橡胶种植气候适宜区的影响是利
大于弊，但对局部区域的橡胶种植可能产生负作用，
如在云南的橡胶主产区景洪、勐腊等地，气候适宜性
将有所降低，因此，需要积极应对气候变化的影响．
当然，影响不同品种橡胶种植分布的主导气候因子
和气候阈值并不相同，同时橡胶种植适宜区的北移
一定程度上会增加寒害的风险，因此，需要根据实际
情况针对不同品种开展橡胶种植区分布的气候适宜
性研究，以取得更为准确的种植分布区信息．需要说
明的是，影响橡胶种植的因素不单是气候因子，还需
要考虑气象灾害的风险、土壤类型、品种的差异、经
营管理和栽培技术的差异等，今后将对此深入研究．
７８０２７期　 　 　 　 　 　 　 　 　 　 　 刘少军等： 未来气候变化对中国天然橡胶种植气候适宜区的影响　 　 　 　 　 　
由于 ＲＣＰ４．５数据是模式输出的对未来气候变
化情景的模拟，是 ＩＰＣＣ 在全球气候变化预测中被
广泛认可的数据集．以目前全球模式的预测水平，必
定与未来真实的情况有偏差，误差必客观存在，难免
会对本研究结果产生一定误差．但迄今暂无更好的
解决方法．
致谢　 感谢中国气象局国家气候中心的段居琦博士和尹宜
舟博士在模型和数据处理方面给予的大力帮助．
参考文献
［１］　 Ｊｉａｎｇ Ｔ （姜　 彤）， Ｌｉ Ｘ⁃Ｃ （李秀仓）， Ｃａｏ Ｑ⁃Ｃ （巢清
尘）， ｅｔ ａｌ． Ｈｉｇｈｌｉｇｈｔｓ ａｎｄ ｕｎｄｅｒｓｔａｎｄｉｎｇ ｏｆ ｃｌｉｍａｔｅ
ｃｈａｎｇｅ ２０１４： Ｉｍｐａｃｔｓ， ａｄａｐｔａｔｉｏｎ， ａｎｄ ｖｕｌｎｅｒａｂｉｌｉｔｙ．
Ａｄｖａｎｃｅｓ ｉｎ Ｃｌｉｍａｔｅ Ｃｈａｎｇｅ Ｒｅｓｅａｒｃｈ （气候变化研究
进展）， ２０１４， １０（３）： １５７－１６６ （ｉｎ Ｃｈｉｎｅｓｅ）
［２］　 Ｈｕ Ｙ⁃Ｌ （胡亚南）， Ｌｉｕ Ｙ⁃Ｊ （刘颖杰）． Ｐｌａｎｔｉｎｇ ｄｉｓｔｒｉ⁃
ｂｕｔｉｏｎ ｏｆ ｓｐｒｉｎｇ ｍａｉｚｅ ａｎｄ ｉｔｓ ｐｒｏｄｕｃｔｉｖｉｔｙ ｕｎｄｅｒ ＲＣＰ４．５
ｓｃｅｎａｒｉｏ ｉｎ Ｎｏｒｔｈｅａｓｔ Ｃｈｉｎａ ｉｎ ２０１１－２０５０． Ｓｃｉｅｎｔｉａ Ａｇ⁃
ｒｉｃｕｌｔｕｒａ Ｓｉｎｉｃａ （中国农业科学）， ２０１３， ４６ （ １５）：
３１０５－３１１４ （ｉｎ Ｃｈｉｎｅｓｅ）
［３］　 Ｑｉ Ｄ⁃Ｌ （祁栋灵）， Ｗａｎｇ Ｘ⁃Ｑ （王秀全）， Ｚｈａｎｇ Ｚ⁃Ｙ
（张志扬）， ｅｔ ａｌ． Ｃｕｒｒｅｎｔ ｓｉｔｕａｔｉｏｎ ｏｆ Ｃｈｉｎｅｓｅ ｎａｔｕｒａｌ
ｒｕｂｂｅｒ ｉｎｄｕｓｔｒｙ ａｎｄ ｄｅｖｅｌｏｐｍｅｎｔ ｓｕｇｇｅｓｔｉｏｎｓ． Ｃｈｉｎｅｓｅ
Ｊｏｕｒｎａｌ ｏｆ Ｔｒｏｐｉｃａｌ Ａｇｒｉｃｕｌｔｕｒｅ （热带农业科学 ），
２０１３， ３３（２）： ７９－８７ （ｉｎ Ｃｈｉｎｅｓｅ）
［４］　 Ｌｉ Ｙ （李　 勇）， Ｙａｎｇ Ｘ⁃Ｇ （杨晓光）， Ｗａｎｇ Ｗ⁃Ｆ （王
文峰）， ｅｔ ａｌ． Ｔｈｅ ｐｏｓｓｉｂｌｅ ｅｆｆｅｃｔｓ ｏｆ ｇｌｏｂａｌ ｗａｒｍｉｎｇ ｏｎ
ｃｒｏｐｐｉｎｇ ｓｙｓｔｅｍｓ ｉｎ Ｃｈｉｎａ． Ⅴ． Ｔｈｅ ｐｏｓｓｉｂｌｅ ｅｆｆｅｃｔｓ ｏｆ
ｃｌｉｍａｔｅ ｗａｒｍｉｎｇ ｏｎ ｇｅｏｇｒａｐｈｉｃａｌ ｓｈｉｆｔ ｉｎ ｓａｆｅ ｎｏｒｔｈｅｒｎ
ｌｉｍｉｔ ｏｆ ｔｒｏｐｉｃａｌ ｃｒｏｐｓ ａｎｄ ｔｈｅ ｒｉｓｋ ａｎａｌｙｓｉｓ ｏｆ ｃｏｌｄ
ｄａｍａｇｅ ｉｎ Ｃｈｉｎａ． Ｓｃｉｅｎｔｉａ Ａｇｒｉｃｕｌｔｕｒａ Ｓｉｎｉｃａ （中国农业
科学）， ２０１０， ４３（１２）： ２４７７－２４８４ （ｉｎ Ｃｈｉｎｅｓｅ）
［５］　 Ｚｈａｏ Ｊ （赵　 锦）， Ｙａｎｇ Ｘ⁃Ｇ （杨晓光）， Ｌｉｕ Ｚ⁃Ｊ （刘
志娟）， ｅｔ ａｌ． Ｔｈｅ ｐｏｓｓｉｂｌｅ ｅｆｆｅｃｔ ｏｆ ｇｌｏｂａｌ ｃｌｉｍａｔｅ ｃｈａｎ⁃
ｇｅｓ ｏｎ ｃｒｏｐｐｉｎｇ ｓｙｓｔｅｍｓ ｂｏｕｎｄａｒｙ ｉｎ Ｃｈｉｎａ． Ⅱ． Ｔｈｅ
ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ ｏｆ ｃｌｉｍａｔｉｃ ｖａｒｉａｂｌｅｓ ａｎｄ ｔｈｅ ｐｏｓｓｉｂｌｅ
ｅｆｆｅｃｔ ｏｎ ｎｏｒｔｈｅｒｎ ｌｉｍｉｔｓ ｏｆ ｃｒｏｐｐｉｎｇ ｓｙｓｔｅｍｓ ｉｎ Ｓｏｕｔｈ
Ｃｈｉｎａ． Ｓｃｉｅｎｔｉａ Ａｇｒｉｃｕｌｔｕｒａ Ｓｉｎｉｃａ （中国农业科学），
２０１０， ４３（９）： １８６０－１８６７ （ｉｎ Ｃｈｉｎｅｓｅ）
［６］　 Ｍｏｓｓ Ｒ， Ｅｄｍｏｎｄｓ Ｊ， Ｈｉｂｂａｒｄ Ｋ， ｅｔ ａｌ． Ｔｈｅ ｎｅｘｔ ｇｅｎｅ⁃
ｒａｔｉｏｎ ｏｆ ｓｃｅｎａｒｉｏｓ ｆｏｒ ｃｌｉｍａｔｅ ｃｈａｎｇｅ ｒｅｓｅａｒｃｈ ａｎｄ ａｓ⁃
ｓｅｓｓｍｅｎｔ． Ｎａｔｕｒｅ， ２００９， ４６３： ７４７－７５６
［７］　 Ｗａｎｇ Ｓ⁃Ｗ （王绍武）， Ｌｕｏ Ｙ （罗 　 勇）， Ｚｈａｏ Ｚ⁃Ｃ
（赵宗慈）， ｅｔ ａｌ． Ｎｅｗ ｇｅｎｅｒａｔｉｏｎ ｏｆ ｓｃｅｎａｒｉｏｓ ｏｆ ｇｒｅｅｎ⁃
ｈｏｕｓｅ ｇａｓ ｅｍｉｓｓｉｏｎ． Ａｄｖａｎｃｅｓ ｉｎ Ｃｌｉｍａｔｅ Ｃｈａｎｇｅ Ｒｅ⁃
ｓｅａｒｃｈ （气候变化研究进展）， ２０１２， ８（４）： ３０５－３０７
（ｉｎ Ｃｈｉｎｅｓｅ）
［８］　 Ｃｈｅｎ Ｍ⁃Ｐ （陈敏鹏）， Ｌｉｎ Ｅ⁃Ｄ （林而达）． Ｇｌｏｂａｌ
ｇｒｅｅｎｈｏｕｓｅ ｇａｓ ｅｍｉｓｓｉｏｎ ｍｉｔｉｇａｔｉｏｎ ｕｎｄｅｒ ｒｅｐｒｅｓｅｎｔａｔｉｖｅ
ｃｏｎｃｅｎｔｒａｔｉｏｎ ｐａｔｈｗａｙｓ ｓｃｅｎａｒｉｏｓ ａｎｄ ｃｈａｌｌｅｎｇｅｓ ｔｏ Ｃｈｉ⁃
ｎａ． Ａｄｖａｎｃｅｓ ｉｎ Ｃｌｉｍａｔｅ Ｃｈａｎｇｅ Ｒｅｓｅａｒｃｈ （气候变化研
究进展）， ２０１０， ６（６）： ４３６－４４２ （ｉｎ Ｃｈｉｎｅｓｅ）
［９］　 Ｚｈａｎ Ｍ⁃Ｊ （占明锦）， Ｙｉｎ Ｊ⁃Ｍ （殷剑敏）， Ｋｏｎｇ Ｐ （孔
萍）， ｅｔ ａｌ． Ｐｒｅｄｉｃｔｉｏｎ ｏｎ ｔｈｅ ｐｏｓｓｉｂｌｅ ｃｌｉｍａｔｅ ｃｈａｎｇｅ ｏｆ
Ｐｏｙａｎｇ Ｌａｋｅ ｂａｓｉｎ ｉｎ ｔｈｅ ｆｕｔｕｒｅ ５０ ｙｅａｒｓ， ｕｎｄｅｒ ＲＣＰ
ｓｃｅｎａｒｉｏ． Ｓｃｉｅｎｃｅ Ｔｅｃｈｎｏｌｏｇｙ ａｎｄ Ｅｎｇｉｎｅｅｒｉｎｇ（科学技
术与工程）， ２０１３， １３（３４）： １０１０７－ １０１１５ （ ｉｎ Ｃｈｉ⁃
ｎｅｓｅ）
［１０］　 Ｎｋｕｌｕｍｏ Ｚ， Ｏｌｉｖｉｅｒ Ｃ， Ｓｅｐｏ Ｈ． Ｃｒｏｐ ｒｅｓｐｏｎｓｅ ｔｏ ｃｌｉ⁃
ｍａｔｅ ｃｈａｎｇｅ ｉｎ ｓｏｕｔｈｅｒｎ Ａｆｒｉｃａ： Ａ ｃｏｍｐｒｅｈｅｎｓｉｖｅ ｒｅ⁃
ｖｉｅｗ． Ｇｌｏｂａｌ ａｎｄ Ｐｌａｎｅｔａｒｙ Ｃｈａｎｇｅ， ２０１３， １１１： １１８－
１２６
［１１］　 Ｂｕｒｋｅ ＭＢ， Ｌｏｂｅｌｌ ＤＢ， Ｇｕａｒｉｎｏ Ｌ． Ｓｈｉｆｔｓ ｉｎ Ａｆｒｉｃａｎ ｃｒｏｐ
ｃｌｉｍａｔｅｓ ｂｙ ２０５０， ａｎｄ ｔｈｅ ｉｍｐｌｉｃａｔｉｏｎｓ ｆｏｒ ｃｒｏｐ ｉｍｐｒｏｖｅ⁃
ｍｅｎｔ ａｎｄ ｇｅｎｅｔｉｃ ｒｅｓｏｕｒｃｅｓ ｃｏｎｓｅｒｖａｔｉｏｎ． Ｇｌｏｂａｌ Ｅｎｖｉｒｏｎ⁃
ｍｅｎｔａｌ Ｃｈａｎｇｅ， ２００９， １９： ３１７－３２５
［１２］　 Ｓｍｉｔｈ ＷＮ， Ｇｒａｎｔ ＢＢ， Ｄｅｓｊａｒｄｉｎｓ ＲＬ， ｅｔ ａｌ． Ａｓｓｅｓｓｉｎｇ
ｔｈｅ ｅｆｆｅｃｔｓ ｏｆ ｃｌｉｍａｔｅ ｃｈａｎｇｅ ｏｎ ｃｒｏｐ ｐｒｏｄｕｃｔｉｏｎ ａｎｄ
ＧＨＧ ｅｍｉｓｓｉｏｎｓ ｉｎ Ｃａｎａｄａ． Ａｇｒｉｃｕｌｔｕｒｅ， Ｅｃｏｓｙｓｔｅｍｓ ａｎｄ
Ｅｎｖｉｒｏｎｍｅｎｔ， ２０１３， １７９： １３９－１５０
［１３］　 Ｓｕｐｉｔａ ＩＣ， ｖａｎ Ｄｉｅｐｅｎ ＡＪ， ｄｅ Ｗｉｔ Ｗ， ｅｔ ａｌ． Ａｓｓｅｓｓｉｎｇ
ｃｌｉｍａｔｅ ｃｈａｎｇｅ ｅｆｆｅｃｔｓ ｏｎ Ｅｕｒｏｐｅａｎ ｃｒｏｐ ｙｉｅｌｄｓ ｕｓｉｎｇ ｔｈｅ
ｃｒｏｐ ｇｒｏｗｔｈ ｍｏｎｉｔｏｒｉｎｇ ｓｙｓｔｅｍ ａｎｄ ａ ｗｅａｔｈｅｒ ｇｅｎｅｒａｔｏｒ．
Ａｇｒｉｃｕｌｔｕｒａｌ ａｎｄ Ｆｏｒｅｓｔ Ｍｅｔｅｏｒｏｌｏｇｙ， ２０１２， １６４： ９６ －
１１１
［１４］　 Ｚｏｍｅｒ ＲＪ， Ｔｒａｂｕｃｃｏ Ａ， Ｗａｎｇ ＭＣ， ｅｔ ａｌ． Ｅｎｖｉｒｏｎｍｅｎｔａｌ
ｓｔｒａｔｉｆｉｃａｔｉｏｎ ｔｏ ｍｏｄｅｌ ｃｈａｎｇｅ ｉｍｐａｃｔｓ ｏｎ ｂｉｏｄｉｖｅｒｓｉｔｙ ａｎｄ
ｒｕｂｂｅｒ ｐｒｏｄｕｃｔｉｏｎ ｉｎ Ｘｉｓｈｕａｎｇｂａｎｎａ， Ｙｕｎｎａｎ， Ｃｈｉｎａ．
Ｂｉｏｌｏｇｉｃａｌ Ｃｏｎｓｅｒｖａｔｉｏｎ， ２０１４， １７０： ２６４－２７３
［１５］　 Ｌｉｕ Ｚ⁃Ｊ （刘志娟）， Ｙａｎｇ Ｘ⁃Ｇ （杨晓光），Ｗａｎｇ Ｗ⁃Ｆ
（王文峰）， ｅｔ ａｌ． Ｔｈｅ ｐｏｓｓｉｂｌｅ ｅｆｆｅｃｔｓ ｏｆ ｇｌｏｂａｌ ｗａｒｍｉｎｇ
ｏｎ ｃｒｏｐｐｉｎｇ ｓｙｓｔｅｍｓ ｉｎ Ｃｈｉｎａ． Ⅳ． Ｔｈｅ ｐｏｓｓｉｂｌｅ ｉｍｐａｃｔ
ｏｆ ｆｕｔｕｒｅ ｃｌｉｍａｔｉｃ ｗａｒｍｉｎｇ ｏｎ ｔｈｅ ｎｏｒｔｈｅｒｎ ｌｉｍｉｔｓ ｏｆ
ｓｐｒｉｎｇ ｍａｉｚｅ ｉｎ ｔｈｒｅｅ ｐｒｏｖｉｎｃｅｓ ｏｆ Ｎｏｒｔｈｅａｓｔ Ｃｈｉｎａ． Ｓｃｉ⁃
ｅｎｔｉａ Ａｇｒｉｃｕｌｔｕｒａ Ｓｉｎｉｃａ （中国农业科学）， ２０１０， ４３
（１１）： ２２８０－２２９１ （ｉｎ Ｃｈｉｎｅｓｅ）
［１６］　 Ｇｏｕ Ｓ⁃Ｗ （苟诗薇）， Ｚｈａｎｇ Ｙ⁃Ｘ （张颖娴）， Ｘｕ Ｙ⁃Ｌ
（许吟隆）． Ａｎａｌｙｓｉｓ ｏｆ ｃｌｉｍａｔｅ ｒｅｓｏｕｒｃｅ ｃｈａｎｇｅｓ ｄｕｒｉｎｇ
ｍａｉｚｅ ｇｒｏｗｔｈ ｐｅｒｉｏｄ ｉｎ Ｎｉｎｇｘｉａ ｕｎｄｅｒ ＳＲＥＳ Ａ１Ｂ
ｓｃｅｎａｒｉｏ． Ｃｈｉｎｅｓｅ Ｊｏｕｒｎａｌ ｏｆ Ｅｃｏ⁃Ａｇｒｉｃｕｌｔｕｒｅ （中国生态
农业学报）， ２０１２， ２０（１０）： １３９４－１４０３ （ｉｎ Ｃｈｉｎｅｓｅ）
［１７］　 Ｗａｎｇ Ｌ⁃Ｘ （王连喜）， Ｒｅｎ Ｊ⁃Ｑ （任景全）， Ｌｉ Ｑ （李
琪）． Ｓｉｍｕｌａｔｉｏｎ ｏｆ ｔｈｅ ｈｅａｔ ｉｎｊｕｒｙ ｏｎ ｒｉｃｅ ｐｒｏｄｕｃｔｉｏｎ ｉｎ
Ｊｉａｎｇｓｕ Ｐｒｏｖｉｎｃｅ ｕｎｄｅｒ ｔｈｅ ｃｌｉｍａｔｅ ｃｈａｎｇｅ ｓｃｅｎａｒｉｏ． Ⅱ．
Ａｄａｐｔａｂｉｌｉｔｙ ａｎａｌｙｓｉｓ ｏｆ ｔｈｅ ｒｉｃｅ ｔｏ ｈｅａｔ ｉｎｊｕｒｙ ｆｒｏｍ ｂｏｏ⁃
ｔｉｎｇ ｔｏ ｈｅａｄｉｎｇ ｓｔａｇｅ． Ｃｈｉｎｅｓｅ Ｊｏｕｒｎａｌ ｏｆ Ａｇｒｏｍｅｔｅｏｒｏｌｏｇｙ
（中国农业气象）， ２０１４， ３５（２）： ２０６－２１３ （ ｉｎ Ｃｈｉ⁃
ｎｅｓｅ）
８８０２ 应　 用　 生　 态　 学　 报　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 ２６卷
［１８］　 Ｌｉ Ｑ （李　 琪）， Ｒｅｎ Ｊ⁃Ｑ （任景全）， Ｗａｎｇ Ｌ⁃Ｘ （王
连喜）． Ｓｉｍｕｌａｔｉｏｎ ｏｆ ｔｈｅ ｈｅａｔ ｉｎｊｕｒｙ ｏｎ ｒｉｃｅ ｐｒｏｄｕｃｔｉｏｎ
ｉｎ Ｊｉａｎｇｓｕ Ｐｒｏｖｉｎｃｅ ｕｎｄｅｒ ｔｈｅ ｃｌｉｍａｔｅ ｃｈａｎｇｅ ｓｃｅｎａｒｉｏ．
Ⅰ． Ｉｍｐａｃｔ ａｓｓｅｓｓｍｅｎｔ ｏｆ ｔｈｅ ｈｅａｔ ｉｎｊｕｒｙ ｏｎ ｒｉｃｅ ｙｉｅｌｄ
ｆｒｏｍ ｂｏｏｔｉｎｇ ｔｏ ｈｅａｄｉｎｇ ｓｔａｇｅ． Ｃｈｉｎｅｓｅ Ｊｏｕｒｎａｌ ｏｆ Ａｇｒｏ⁃
ｍｅｔｅｏｒｏｌｏｇｙ （中国农业气象）， ２０１４， ３５（１）： ９１－９６
（ｉｎ Ｃｈｉｎｅｓｅ）
［１９］　 Ｐａｎｇ Ｙ⁃Ｍ （庞艳梅）， Ｃｈｅｎ Ｃ （陈 　 超）， Ｐａｎ Ｘ⁃Ｂ
（潘学标）． Ｃｈａｎｇｅ ｏｆ ｃｌｉｍａｔｉｃ ｒｅｓｏｕｒｃｅｓ ａｎｄ ｐｏｔｅｎｔｉａｌ
ｐｒｏｄｕｃｔｉｖｉｔｙ ｄｕｒｉｎｇ ｔｈｅ ｅｎｔｉｒｅ ｇｒｏｗｔｈ ｐｅｒｉｏｄ ｏｆ ｗｉｎｔｅｒ
ｗｈｅａｔ ｕｎｄｅｒ ｆｕｔｕｒｅ ｃｌｉｍａｔｅ ｓｃｅｎａｒｉｏ ｉｎ Ｓｉｃｈｕａｎ ｂａｓｉｎ．
Ｃｈｉｎｅｓｅ Ｊｏｕｒｎａｌ ｏｆ Ａｇｒｏｍｅｔｅｏｒｏｌｏｇｙ （中国农业气象），
２０１４， ３５（１）： １－９ （ｉｎ Ｃｈｉｎｅｓｅ）
［２０］　 Ｚｈａｏ Ｊ⁃Ｆ （赵俊芳）， Ｇｕｏ Ｊ⁃Ｐ （郭建平）， Ｗｕ Ｄ⁃Ｒ （邬
定荣）， ｅｔ ａｌ． Ｃｌｉｍａｔｉｃ ｐｏｔｅｎｔｉａｌ ｐｒｏｄｕｃｔｉｖｉｔｙ ｏｆ ｗｉｎｔｅｒ
ｗｈｅａｔ ａｎｄ ｓｕｍｍｅｒ ｍａｉｚｅ ｉｎ Ｈｕａｎｇｈｕａｉｈａｉ Ｐｌａｉｎ ｉｎ
２０１１－２０５０． Ｃｈｉｎｅｓｅ Ｊｏｕｒｎａｌ ｏｆ Ａｐｐｌｉｅｄ Ｅｃｏｌｏｇｙ （应用
生态学报）， ２０１１， ２２（１２）： ３１８９－３１９５ （ｉｎ Ｃｈｉｎｅｓｅ）
［２１］　 Ｇａｏ Ｔ （高　 涛）， Ｙａｎｇ Ｚ⁃Ｌ （杨泽龙）， Ｗｅｉ Ｙ⁃Ｒ （魏
玉荣）， ｅｔ ａｌ． Ｐｒｅ⁃ｅｓｔｉｍａｔｉｎｇ ｔｈｅ ｙｉｅｌｄｓ ｏｆ ｔｈｒｅｅ ｍａｊｏｒ
ｃｒｏｐｓ ｉｎ Ｉｎｎｅｒ Ｍｏｎｇｏｌｉａ ｕｎｄｅｒ ｔｈｅ ｎｕｍｅｒｉｃａｌ ｓｉｍｕｌａｔｉｏｎ
ｓｃｅｎａｒｉｏｓ ｆｏｒ ｔｈｅ ｆｕｔｕｒｅ ３０ ｙｅａｒｓ． Ｃｈｉｎｅｓｅ Ｊｏｕｒｎａｌ ｏｆ
Ａｇｒｏｍｅｔｅｏｒｏｌｏｇｙ （中国农业气象）， ２０１３， ３４ （ ６）：
６８５－６９５ （ｉｎ Ｃｈｉｎｅｓｅ）
［２２］　 Ｄｕ Ｙ⁃Ｄ （杜尧东）， Ｄｕａｎ Ｈ⁃Ｌ （段海来）， Ｔａｎｇ Ｌ⁃Ｓ
（唐力生）． Ａｄａｐｔａｂｉｌｉｔｙ ｏｆ ｃｉｔｒｕｓ ｉｎ ｓｕｂｔｒｏｐｉｃｓ ｏｆ Ｃｈｉｎａ
ｕｎｄｅｒ ｆｕｔｕｒｅ ｃｌｉｍａｔｉｃ ｓｃｅｎａｒｉｏ． Ｃｈｉｎｅｓｅ Ｊｏｕｒｎａｌ ｏｆ Ｅｃｏｌｏ⁃
ｇｙ （生态学杂志）， ２０１０， ２９（５）： ８３３－８３９ （ ｉｎ Ｃｈｉ⁃
ｎｅｓｅ）
［２３］　 Ｃｏｏｋ ＷＣ． Ｎｏｔｅｓ ｏｎ ｐｒｅｄｉｃｔｉｎｇ ｔｈｅ ｐｒｏｂａｂｌｅ ｆｕｔｕｒｅ ｄｉｓｔｒｉ⁃
ｂｕｔｉｏｎ ｏｆ ｉｎｔｒｏｄｕｃｅｄ ｉｎｓｅｃｔｓ． Ｅｃｏｌｏｇｙ， １９３１， １２： ２４５－
２４７
［２４］　 Ｗｅｉ Ｓ⁃Ｘ （魏淑秋）， Ｌｉｕ Ｇ⁃Ｌ （刘桂莲）． Ｓｔｕｄｙ ｏｆ Ｂｉｏ⁃
ｌｏｇｉｃ Ｃｌｉｍａｔｅ Ｓｉｍｉｌａｒｉｔｙ ｂｅｔｗｅｅｎ Ｃｈｉｎａ ａｎｄ ｔｈｅ Ｗｏｒｌｄ．
Ｂｅｉｊｉｎｇ： Ｏｃｅａｎ Ｐｒｅｓｓ， １９９４ （ｉｎ Ｃｈｉｎｅｓｅ）
［２５］　 Ｗｅｉ Ｓ⁃Ｘ （魏淑秋）． Ｉｎｔｒｏｄｕｃｔｉｏｎ ｏｆ ａｇｒｉｃｕｌｔｕｒａｌ ｃｌｉｍａｔｅ
ｓｉｍｉｌａｒ ｄｉａｇｎｏｓｉｓ ａｎａｌｙｓｉｓ． Ｊｏｕｒｎａｌ ｏｆ Ｂｅｉｊｉｎｇ Ａｇｒｉｃｕｌｔｕｒａｌ
Ｕｎｉｖｅｒｓｉｔｙ （北京农业大学学报）， １９８４， １０（４）： ４２７－
４２８ （ｉｎ Ｃｈｉｎｅｓｅ）
［２６］　 Ｊｉａｎｇ Ｑ （蒋 　 青）． Ｓｔｕｄｙ ｏｎ ｔｈｅ ｓｕｉｔａｂｌｅ ｄｉｓｔｒｉｂｕｔｉｏｎ
ａｒｅａ ｏｆ Ｊｏｈｎｓｏｎｇｒａｓｓ ｉｎ ｏｕｒ ｃｏｕｎｔｒｙ ｂｙ ａｇｒｉｃｕｌｔｕｒａｌ ｃｌｉ⁃
ｍａｔｅ ａｎａｌｏｇｉｅｓ ｄｉｓｔａｎｃｅ． Ｐｌａｎｔ Ｑｕａｒａｎｔｉｎｅ （植物检
疫）， １９９４， ８（５）： ２５７－２６２ （ｉｎ Ｃｈｉｎｅｓｅ）
［２７］　 Ｌｉｕ Ｈ⁃Ｘ （刘红霞）， Ｗｅｎ Ｊ⁃Ｂ （温俊宝）， Ｌｕｏ Ｙ⁃Ｑ
（骆有庆）， ｅｔ ａｌ． Ｒｅｖｉｅｗ ｏｎ ｐｅｓｔ ｒｉｓｋ ａｎａｌｙｓｉｓ ｏｆ ｆｏｒｅｓｔ
ｐｌａｎｔｓ． Ｊｏｕｒｎａｌ ｏｆ Ｂｅｉｊｉｎｇ Ｆｏｒｅｓｔｒｙ Ｕｎｉｖｅｒｓｉｔｙ（北京林业
大学学报）， ２００１， ２３（６）： ４６－５１ （ｉｎ Ｃｈｉｎｅｓｅ）
［２８］　 Ｆａｎ Ｊ⁃Ａ （范京安）． Ｓｔｕｄｉｅｓ ｏｎ ｔｈｅ ａｓｓｅｓｓｍｅｎｔ ｏｆ ｐｏｔｅｎ⁃
ｔｉａｌ ｅｓｔａｂｌｉｓｈｍｅｎｔ ｏｆ ｏｒｉｅｎｔａｌ ｆｒｕｉｔ （Ｂａｃｔｒｏｃｅｒａ ｄｏｒｓａｌｉｓ）
ｉｎ Ｃｈｉｎａ ｗｉｔｈ ｓｙｎｔｈｅｔｉｃ ｊｕｄｇｍｅｎｔ ｏｎ ｆｕｚｚｙ ｍａｔｈｅｍａｔｉｃｓ．
Ｐｌａｎｔ Ｑｕａｒａｎｔｉｎｅ （植物检疫）， １９９８， １２（２）： ７６－８０
（ｉｎ Ｃｈｉｎｅｓｅ）
［２９］ 　 Ｓｃｈａｄｔ Ｓ， Ｒｅｖｉｌｌａ Ｅ， Ｗｉｅｇａｎｄ Ｔ， ｅｔ ａｌ． Ａｓｓｅｓｓｉｎｇ ｔｈｅ
ｓｕｉｔａｂｉｌｉｔｙ ｏｆ ｃｅｎｔｒａｌ Ｅｕｒｏｐｅａｎ ｌａｎｄｓｃａｐｅｓ ｆｏｒ ｔｈｅ ｒｅｉｎｔｒｏ⁃
ｄｕｃｔｉｏｎ ｏｆ Ｅｕｒａｓｉａｎ ｌｙｎｘ． Ｊｏｕｒｎａｌ ｏｆ Ａｐｐｌｉｅｄ Ｅｃｏｌｏｇｙ，
２００２， ３９： １８９－２０３
［３０］　 Ａｕｓｔｉｎ ＭＰ， Ｍｅｙｅｒｓ ＪＡ． Ｃｕｒｒｅｎｔ ａｐｐｒｏａｃｈｅｓ ｔｏ ｍｏｄｅｌｉｎｇ
ｔｈｅ ｅｎｖｉｒｏｎｍｅｎｔａｌ ｎｉｃｈｅ ｏｆ ｅｕｃａｌｙｐｔｓ： Ｉｍｐｌｉｃａｔｉｏｎｓ ｆｏｒ
ｍａｎａｇｅｍｅｎｔ ｏｆ ｆｏｒｅｓｔ ｂｉｏｄｉｖｅｒｓｉｔｙ． Ｆｏｒｅｓｔ Ｅｃｏｌｏｇｙ ａｎｄ
Ｍａｎａｇｅｍｅｎｔ， １９９６， ８５： ９５－１０６
［３１］　 Ｂｒｉｔｏ ＪＣ， Ｃｒｅｓｐｏ ＥＧ， Ｐａｕｌｏ ＯＳ． Ｍｏｄｅｌｉｎｇ ｗｉｌｄｌｉｆｅ ｄｉｓ⁃
ｔｒｉｂｕｔｉｏｎｓ： Ｌｏｇｉｓｔｉｃ ｍｕｌｔｉｐｌｅ ｒｅｇｒｅｓｓｉｏｎ ｖｓ ｏｖｅｒｌａｐ ａｎａｌｙ⁃
ｓｉｓ． Ｅｃｏｇｒａｐｈｙ， １９９９， ２２： ２５１－２６０
［３２］　 Ｆｒｅｓｃｉｎｏ ＴＳ， Ｅｄｗａｒｄｓ ＴＣ， Ｍｏｉｓｅｎ ＧＧ． Ｍｏｄｅｌｉｎｇ ｓｐａ⁃
ｔｉａｌｌｙ ｅｘｐｌｉｃｉｔ ｆｏｒｅｓｔ ｓｔｒｕｃｔｕｒａｌ ａｔｔｒｉｂｕｔｅｓ ｕｓｉｎｇ Ｇｅｎｅｒａ⁃
ｌｉｚｅｄ Ａｄｄｉｔｉｖｅ Ｍｏｄｅｌｓ． Ｊｏｕｒｎａｌ ｏｆ Ｖｅｇｅｔａｔｉｏｎ Ｓｃｉｅｎｃｅ，
２００１， １２： １５－２６
［３３］　 Ｖａｙｓｓｉｒｅｓ ＭＰ， Ｐｌａｎｔ ＲＥ， Ａｌｌｅｎ⁃Ｄｉａｚ ＢＨ． Ｃｌａｓｓｉｆｉｃａｔｉｏｎ
ｔｒｅｅｓ： Ａｎ ａｌｔｅｒｎａｔｉｖｅ ｎｏｎ⁃ｐａｒａｍｅｔｒｉｃ ａｐｐｒｏａｃｈ ｆｏｒ ｐｒｅ⁃
ｄｉｃｔｉｎｇ ｓｐｅｃｉｅｓ ｄｉｓｔｒｉｂｕｔｉｏｎｓ． Ｊｏｕｒｎａｌ ｏｆ Ｖｅｇｅｔａｔｉｏｎ Ｓｃｉ⁃
ｅｎｃｅ， ２０００， １１： ６７９－６９４
［３４］　 Ｃａｏ Ｍ⁃Ｃ （曹铭昌）， Ｚｈｏｕ Ｇ⁃Ｓ （周广胜）， Ｗａｎｇ Ｅ⁃Ｓ
（翁恩生）． Ａｐｐｌｉｃａｔｉｏｎ ａｎｄ ｃｏｍｐａｒｉｓｏｎ ｏｆ ｇｅｎｅｒａｌｉｚｅｄ
ｍｏｄｅｌｓ ａｎｄ ｃｌａｓｓｉｆｉｃａｔｉｏｎ ａｎｄ ｒｅｇｒｅｓｓｉｏｎ ｔｒｅｅ ｉｎ ｓｉｍｕｌａ⁃
ｔｉｎｇ ｔｒｅｅ ｓｐｅｃｉｅｓ ｄｉｓｔｒｉｂｕｔｉｏｎ． Ａｃｔａ Ｅｃｏｌｏｇｉｃａ Ｓｉｎｉｃａ （生
态学报）， ２００５， ２５（８）： ２０３１－２０４０ （ｉｎ Ｃｈｉｎｅｓｅ）
［３５］　 Ｓｕｔｈｅｒｓｔ ＲＷ， Ｍａｙｗａｌｄ ＧＦ． Ａ ｃｏｍｐｕｔｅｒｉｚｅｄ ｓｙｓｔｅｍ ｆｏｒ
ｍａｔｃｈｉｎｇ ｃｌｉｍａｔｅｓ ｉｎ ｅｃｏｌｏｇｙ． Ａｇｒｉｃｕｌｔｕｒｅ， Ｅｃｏｓｙｓｔｅｍｓ ＆
Ｅｎｖｉｒｏｎｍｅｎｔ， １９８５， １３： ２８１－２９９
［３６］　 Ｗｏｍｅｒ ＳＰ． Ｅｃｏｃｌｉｍａｔｉｃ ａｓｓｅｓｓｍｅｎｔ ｏｆ ｐｏｔｅｎｔｉａｌ ｅｓｔａｂ⁃
ｌｉｓｈｍｅｎｔ ｏｆ ｅｘｏｔｉｃ ｐｅｓｔｓ． Ｊｏｕｒｎａｌ ｏｆ Ｅｃｏｎｏｍｉｃ Ｅｎｔｏｍｏｌｏ⁃
ｇｙ， １９８８， ８１： ９７３－９８３
［３７］　 Ｓｔｅｖｅｎ ＪＰ， Ｒｏｂｅｒｔ ＰＡ， Ｒｏｂｅｒｔ ＥＳ． Ｍａｘｉｍｕｍ ｅｎｔｒｏｐｙ
ｍｏｄｅｌｉｎｇ ｏｆ ｓｐｅｃｉｅｓ ｇｅｏｇｒａｐｈｉｃ ｄｉｓｔｒｉｂｕｔｉｏｎ． Ｅｃｏｌｏｇｉｃａｌ
Ｍｏｄｅｌｌｉｎｇ， ２００６， １９０： ２３１－２５９
［３８］　 Ｓｔｅｖｅｎ ＪＰ， Ｍｉｒｏｓｌａｖ Ｄ． Ｍｏｄｅｌｉｎｇ ｏｆ ｓｐｅｃｉｅｓ ｄｉｓｔｒｉｂｕｔｉｏｎｓ
ｗｉｔｈ Ｍａｘｅｎｔ： Ｎｅｗ ｅｘｔｅｎｓｉｏｎｓ ａｎｄ ａ ｃｏｍｐｒｅｈｅｎｓｉｖｅ
ｅｖａｌｕａｔｉｏｎ． Ｅｃｏｇｒａｐｈｙ， ２００８， ３１： １６１－１７５
［３９］　 Ｆｉｓｃｈｅｒ Ｊ， Ｌｉｎｄｅｎｍａｙｅｒ ＤＢ， Ｎｉｘ ＨＡ， ｅｔ ａｌ． Ｃｌｉｍａｔｅ
ａｎｄ ａｎｉｍａｌ ｄｉｓｔｒｉｂｕｔｉｏｎ： Ａ ｃｌｉｍａｔｉｃ ａｎａｌｙｓｉｓ ｏｆ ｔｈｅ Ａｕｓ⁃
ｔｒａｌｉａｎ ｍａｒｓｕｐｉａｌ Ｔｒｉｃｈｏｓｕｒｕｓ ｃａｎｉｎｕｓ． Ｊｏｕｒｎａｌ ｏｆ Ｂｉｏｇｅ⁃
ｏｇｒａｐｈｙ， ２００１， ２８： ２９３－３０４
［４０］　 Ｓｈａｏ Ｈ （邵　 慧）， Ｔｉａｎ Ｊ⁃Ｑ （田佳倩）， Ｇｕｏ Ｋ （郭　
柯）， ｅｔ ａｌ． Ｅｆｆｅｃｔｓ ｏｆ ｓａｍｐｌｅ ｓｉｚｅ ａｎｄ ｓｐｅｃｉｅｓ ｔｒａｉｔｓ ｏｎ
ｐｅｒｆｏｒｍａｎｃｅ ｏｆ ＢＩＯＣＬＩＭ ｉｎ ｐｒｅｄｉｃｔｉｎｇ ｇｅｏｇｒａｐｈｉｃａｌ ｄｉｓ⁃
ｔｒｉｂｕｔｉｏｎ ｏｆ ｔｒｅｅ ｓｐｅｃｉｅｓ： Ａ ｃａｓｅ ｓｔｕｄｙ ｗｉｔｈ １２ ｄｅｃｉｄｕｏｕｓ
Ｑｕｅｒｃｕｓ ｓｐｅｃｉｅｓ ｉｎｄｉｇｅｎｏｕｓ ｔｏ Ｃｈｉｎａ． Ｃｈｉｎｅｓｅ Ｊｏｕｒｎａｌ ｏｆ
Ｐｌａｎｔ Ｅｃｏｌｏｇｙ （植物生态学报）， ２００９， ３３（５）： ８７０－
８７７ （ｉｎ Ｃｈｉｎｅｓｅ）
［４１］ 　 Ｗａｎｇ Ｌ （王 　 菱）． Ｇｅｏｇｒａｐｈｉｃａｌ ｅｎｖｉｒｏｎｍｅｎｔ ｅｖａｌｕａ⁃
９８０２７期　 　 　 　 　 　 　 　 　 　 　 刘少军等： 未来气候变化对中国天然橡胶种植气候适宜区的影响　 　 　 　 　 　
ｔｉｏｎ ｏｆ ｎｏｒｔｈｅｒｎ ｇｒｏｗｔｈ ｂｏｕｎｄａｒｙ ｏｆ ｒｕｂｂｅｒ ｐｌａｎｔａｔｉｏｎ ｉｎ
Ｃｈｉｎａ． Ｊｏｕｒｎａｌ ｏｆ Ｎａｔｕｒａｌ Ｒｅｓｏｕｒｃｅｓ （自然资源学报），
１９８７， １１（２）： ５４－６１ （ｉｎ Ｃｈｉｎｅｓｅ）
［４２］　 Ｃｏｏｐｅｒａｔｉｖｅ Ｇｒｏｕｐ ｏｆ Ａｇｒｉｃｕｌｔｕｒｅ， Ｆｏｒｅｓｔｒｙ ａｎｄ Ｃｒｏｐ Ｒｅ⁃
ｇｉｏｎａｌｉｚａｔｉｏｎ ｉｎ Ｃｈｉｎａ （中国农林作物气候区划协作
组）． Ｃｌｉｍａｔｉｃ Ｒｅｇｉｏｎａｌｉｚａｔｉｏｎ ｏｆ Ａｇｒｉｃｕｌｔｕｒｅ， Ｆｏｒｅｓｔｒｙ
ａｎｄ Ｃｒｏｐ ｉｎ Ｃｈｉｎａ． Ｂｅｉｊｉｎｇ： Ｃｈｉｎａ Ｍｅｔｅｏｒｏｌｏｇｉｃａｌ Ｐｒｅｓｓ，
１９８７ （ｉｎ Ｃｈｉｎｅｓｅ）
［４３］　 Ｗａｎｇ Ｌ⁃Ｆ （王利溥）． Ｒｕｂｂｅｒ Ｍｅｔｅｏｒｏｌｏｇｙ． Ｂｅｉｊｉｎｇ：
Ｃｈｉｎａ Ｍｅｔｅｏｒｏｌｏｇｉｃａｌ Ｐｒｅｓｓ， １９８９ （ｉｎ Ｃｈｉｎｅｓｅ）
［４４］　 Ｚｈａｏ’ ａｎ Ｃｏｕｎｔｙ Ｒｕｂｂｅｒ Ｓｔａｔｉｏｎ Ｄｉｖｉｓｉｏｎ Ｇｒｏｕｐ （诏安
县橡胶站区划组）． Ｒｅｐｏｒｔ ｏｆ Ｚｈａｏ’ ａｎ Ｃｏｕｎｔｙ ｒｕｂｂｅｒ
ｐｒｏｄｕｃｔｉｏｎ ａｎｄ ｄｉｖｉｓｉｏｎ． Ｆｕｊｉａｎ Ｓｃｉｅｎｃｅ ＆ Ｔｅｃｈｎｏｌｏｇｙ ｏｆ
Ｔｒｏｐｉｃａｌ Ｃｒｏｐｓ （福建热作科技）， １９８５， １０（３）： １－９
（ｉｎ Ｃｈｉｎｅｓｅ）
［４５］　 Ｔｈｏｍｓｏｎ ＡＭ， Ｃａｌｖｉｎ ＫＶ， Ｓｍｉｔｈ ＳＪ， ｅｔ ａｌ． ＲＣＰ４．５： Ａ
ｐａｔｈｗａｙ ｆｏｒ ｓｔａｂｉｌｉｚａｔｉｏｎ ｏｆ ｒａｄｉａｔｉｖｅ ｆｏｒｃｉｎｇ ｂｙ ２１００．
Ｃｌｉｍａｔｉｃ Ｃｈａｎｇｅ， ２０１１， １０９： ７７－９４
［４６］　 Ｗａｎｇ Ｙ⁃Ｙ （王言英）， Ｈａｏ Ｌ （赫　 亮）． Ｍａｘｉｍｕｍ ｅｎ⁃
ｔｒｏｐｙ ｍｅｔｈｏｄ ａｎｄ ｉｔｓ ａｐｐｌｉｃａｔｉｏｎ ｏｎ ｓｐｅｃｔｒａｌ ａｎａｌｙｓｉｓ．
Ｊｏｕｒｎａｌ ｏｆ Ｄａｌｉａｎ Ｍａｒｉｔｉｍｅ Ｕｎｉｖｅｒｓｉｔｙ （大连海事大学
学报）， ２０１１， ３７（１）： １－４ （ｉｎ Ｃｈｉｎｅｓｅ）
［４７］　 Ｓｕｎ Ｙ （孙　 瑜）， Ｓｈｉ Ｍ⁃Ｃ （史明昌）， Ｐｅｎｇ Ｈ （彭　
欢）， ｅｔ ａｌ． Ｆｏｒｅｓｔ ｌｉｇｈｔｉｎｇ ｆｉｒｅ ｆｏｒｅｃａｓｔｉｎｇ ｆｏｒ Ｄａｘｉｎｇ’
ａｎｌｉｎｇ Ｍｏｕｎｔａｉｎｓ ｂａｓｅｄ ｏｎ ＭＡＸＥＮＴ ｍｏｄｅｌ． Ｃｈｉｎｅｓｅ
Ｊｏｕｒｎａｌ ｏｆ Ａｐｐｌｉｅｄ Ｅｃｏｌｏｇｙ （应用生态学报）， ２０１４， ２５
（４）： １１００－１１０６ （ｉｎ Ｃｈｉｎｅｓｅ）
［４８］　 Ｓｕｎ Ｊ⁃Ｓ （孙敬松）， Ｚｈｏｕ Ｇ⁃Ｓ （周广胜）． Ｉｎｔｅｒ⁃ｄｅｃａｄａｌ
ｖａｒｉａｂｉｌｉｔｙ ｏｆ ｗｉｎｔｅｒ ｗｈｅａｔ ｐｌａｎｔｉｎｇ ｚｏｎｅ ｉｎ Ｃｈｉｎａ ｄｕｒｉｎｇ
１９６１ ｔｏ ２０１０ ｓｉｍｕｌａｔｅｄ ｂｙ ｍａｘｉｍｕｍ ｅｎｔｒｏｐｙ （ＭａｘＥｎｔ）．
Ｃｈｉｎｅｓｅ Ｊｏｕｒｎａｌ ｏｆ Ａｇｒｏｍｅｔｅｏｒｏｌｏｇｙ （中国农业气象），
２０１２， ３３（４）： ４８１－４８７ （ｉｎ Ｃｈｉｎｅｓｅ）
［４９］　 Ｐｅｔｅｒｓｏｎ ＡＴ， Ｐａｐｅｓ Ｍ， Ｅａｔｏｎ Ｍ． Ｔｒａｎｓｆｅｒａｂｉｌｉｔｙ ａｎｄ
ｍｏｄｅｌ ｅｖａｌｕａｔｉｏｎ ｉｎ ｅｃｏｌｏｇｉｃａｌ ｎｉｃｈｅ ｍｏｄｅｌｉｎｇ： Ａ ｃｏｍ⁃
ｐａｒｉｓｏｎ ｏｆ ＧＡＲＰ ａｎｄ ＭａｘＥｎ． Ｅｃｏｇｒａｐｈｙ， ２００７， ３０：
５５０－５６０
［５０］　 Ｄｕａｎ Ｊ⁃Ｑ （段居琦）， Ｚｈｏｕ Ｇ⁃Ｓ （周广胜）． Ｃｌｉｍａｔｉｃ
ｓｕｉｔａｂｉｌｉｔｙ ｏｆ ｄｏｕｂｌｅ ｒｉｃｅ ｐｌａｎｔｉｎｇ ｒｅｇｉｏｎｓ ｉｎ Ｃｈｉｎａ． Ｓｃｉ⁃
ｅｎｔｉａ Ａｇｒｉｃｕｌｔｕｒａ Ｓｉｎｉｃａ （中国农业科学）， ２０１２， ４５
（２）： ２１８－２２７ （ｉｎ Ｃｈｉｎｅｓｅ）
［５１］　 Ｗａｎｇ Ｙ⁃Ｓ （王运生）， Ｘｉｅ Ｂ⁃Ｙ （谢丙炎）， Ｗａｎ Ｆ⁃Ｈ
（万方浩）， ｅｔ ａｌ． Ａｐｐｌｉｃａｔｉｏｎ ｏｆ ＲＯＣ ｃｕｒｖｅ ａｎａｌｙｓｉｓ ｉｎ
ｅｖａｌｕａｔｉｎｇ ｔｈｅ ｐｅｒｆｏｒｍａｎｃｅ ｏｆ ａｌｉｅｎ ｓｐｅｃｉｅｓ’ ｐｏｔｅｎｔｉａｌ
ｄｉｓｔｒｉｂｕｔｉｏｎ ｍｏｄｅｌｓ． Ｂｉｏｄｉｖｅｒｓｉｔｙ Ｓｃｉｅｎｃｅ （生物多样
性）， ２００７， １５（４）： ３６５－３７２ （ｉｎ Ｃｈｉｎｅｓｅ）
［５２］　 Ｃｈｅ Ｌ （车　 乐）， Ｃａｏ Ｂ （曹 　 博）， Ｂａｉ Ｃ⁃Ｋ （白成
科）， ｅｔ ａｌ． Ｐｒｅｄｉｃｔｉｖｅ ｄｉｓｔｒｉｂｕｔｉｏｎ ａｎｄ ｈａｂｉｔａｔ ｓｕｉｔａｂｉｌｉｔｙ
ａｓｓｅｓｓｍｅｎｔ ｏｆ Ｎｏｔｈｏｌｉｒｉｏｎ ｂｕｌｂｕｌｉｆｅｒｕｍ ｂａｓｅｄ ｏｎ ＭａｘＥｎｔ
ａｎｄ ＡｒｃＧＩＳ． Ｃｈｉｎｅｓｅ Ｊｏｕｒｎａｌ ｏｆ Ｅｃｏｌｏｇｙ （生态学杂
志）， ２０１４， ３３（６）： １－６ （ｉｎ Ｃｈｉｎｅｓｅ）
［５３］　 Ｚｈｅｎｇ Ｗ⁃Ｒ （郑文荣）． Ｔｈｅ Ｄｅｖｅｌｏｐｍｅｎｔ ｏｆ Ｃｈｉｎａ’ ｓ
Ｒｕｂｂｅｒ ａｎｄ Ｌａｔｅｘ Ｔｒｅｎｄ ［ ＥＢ ／ ＯＬ ］． （ ２０１１⁃０８⁃１６ ）
［ ２０１４⁃０６⁃３０］． ｈｔｔｐ： ／ ／ ｗｗｗ． ｄｏｃｉｎ． ｃｏｍ ／ ｐ⁃２４５９４４８６９．
ｈｔｍｌ （ｉｎ Ｃｈｉｎｅｓｅ）
［５４］　 Ｏｆｆｉｃｅ ｏｆ Ｒｅｇｉｏｎａｌ Ｐｌａｎｎｉｎｇ ｏｆ Ｔｒｏｐｉｃａｌ Ｃｒｏｐｓ ｏｆ Ｐａｓｔｏｒａｌ
Ｆｉｓｈｅｒｙ Ｍｉｎｉｓｔｒｙ （农牧渔业部热带作物区划办公室）．
Ｒｅｇｉｏｎａｌ Ｐｌａｎｎｉｎｇ ｏｆ Ｔｒｏｐｉｃａｌ Ｃｒｏｐｓ ｉｎ Ｃｈｉｎａ． Ｇｕａｎｇ⁃
ｚｈｏｕ： Ｇｕａｎｇｄｏｎｇ Ｓｃｉｅｎｃｅ Ｐｒｅｓｓ， １９８９ （ｉｎ Ｃｈｉｎｅｓｅ）
作者简介　 刘少军，男，１９８０ 年生，博士，副研究员． 主要从
事应用气象研究． Ｅ⁃ｍａｉｌ： ｃｄｕｔｌｓｊ＠ １６３．ｃｏｍ
责任编辑　 杨　 弘
０９０２ 应　 用　 生　 态　 学　 报　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 ２６卷