删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

中国生态过渡带分布的空间识别及情景模拟

本站小编 Free考研考试/2021-12-29

范泽孟,1,2,31.中国科学院地理科学与资源研究所 资源与环境信息系统国家重点实验室,北京100101
2.中国科学院大学资源与环境学院,北京100049
3.江苏省地理信息资源开发与利用协同创新中心,南京 210023

Spatial identification and scenario simulation of ecotone distribution in China

FAN Zemeng,1,2,31. State Key Laboratory of Resources and Environmental Information System, Institute of Geographic Sciences and Natural Resources Research, CAS, Beijing 100101, China
2. College of Resources and Environment, University of Chinese Academy of Sciences, Beijing 100049, China
3. Jiangsu Center for Collaborative Innovation in Geographical Information Resource Development and Application, Nanjing 210023, China

收稿日期:2020-10-19修回日期:2021-03-9网络出版日期:2021-03-25
基金资助:国家重点研发计划.2018YFC0507202
国家重点研发计划.2017YFA0603702
国家自然科学基金项目.41971358
国家自然科学基金项目.41930647
中国科学院先导专项.XDA20030203
资源与环境信息系统国家重点实验室自主部署创新研究计划


Received:2020-10-19Revised:2021-03-9Online:2021-03-25
Fund supported: National Key R&D Program of China.2018YFC0507202
National Key R&D Program of China.2017YFA0603702
National Natural Science Foundation of China.41971358
National Natural Science Foundation of China.41930647
Strategic Priority Research Program of the Chinese Academy of Sciences.XDA20030203
Innovation Research Project of State Key Laboratory of Resources and Environment Information System, CAS

作者简介 About authors
范泽孟(1977-), 男, 云南镇雄人, 博士, 副研究员, 中国地理学会会员(S110007694M), 研究方向为气候变化与生态系统响应、生态模型与系统模拟、土地覆盖情景预测及生物多样性模拟。E-mail: fanzm@lreis.ac.cn



摘要
在全球变化及其生态环境效应研究中,如何对生态过渡带的空间分布格局及变化情景进行空间定量识别和模拟分析,对揭示气候变化和人类活动对全球变化的响应及反馈具有指示性意义。在对HLZ模型进行修正和拓展的基础上,建立了生态过渡带类型的空间识别方法。并基于1981—2010年的全国782个气候观测站点数据,在实现全国生态过渡带类型及分布的空间识别基础上,结合3种气候情景数据CMIP5 RCP 2.6、RCP 4.5和RCP 8.5,实现了T0(1981—2010年)、T1(2011—2040年)、T2(2041—2070年)和T3(2071—2100年)4个时段内全国生态过渡带的空间分布格局及其未来情景模拟。另外,引入平均中心空间分析模型,对全国生态过渡带平均中心的时空偏移趋势进行了定量分析。结果显示:在T0~T3时段内,全国共出现41种生态过渡带类型,约占全国陆地面积的18%;冷温带草原/湿润森林与暖温带干旱森林过渡带(564238.5 km2)、冷温带湿润森林与暖温带干旱/湿润森林过渡带(566549.75 m2)、北方湿润/潮湿森林与冷温带湿润森林过渡带(525750.25 km2)是最主要的3种生态过渡带类型。面积占到全国生态过渡带总面积的35%;2010—2100年期间的冷温带荒漠灌丛与暖温带荒漠灌丛/有刺草原过渡带的增加速度最快,在3种情景RCP 2.6、RCP 4.5和RCP 8.5下,其面积将分别增加3604.2 km2/10a、10063.1 km2/10a和17242 km2/10a;寒冷型生态过渡带类型总体上呈向暖湿型过渡带类型增加的趋势;北方潮湿森林与冷温带湿润/潮湿森林过渡带的平均中心偏移幅度最大,在4个时段内整体向东北方向偏移,其偏移幅度将超过150 km。另外,随着气温的逐渐上升和降水量的增加,中国北方的生态过渡带整体呈向北偏移趋势,南方生态过渡带则逐渐减少且平均中心呈现逐渐向高海拔地区退缩的趋势,气候变化对青藏高原区生态过渡带时空格局的影响日益显著。
关键词: 生态过渡带;平均中心;空间识别模型;情景模拟;中国

Abstract
Explicitly identifying the ecotone distribution and scenario change is of important significance to understand the response of ecosystem to climatic change. In this paper, a spatial identification method was developed to analyze the ecotone distribution in terms of the improved Holdridge life zone (iHLZ) model. Based on the climatic observation data of 782 weather stations of China in the T0 (1981-2010) period, and the climatic scenario data of IPCC CMIP5 RCP2.6, RCP4.5 and RCP8.5 in T1 (2011-2040), T2 (2041-2070) and T3 (2071-2100), the ecotones distribution and scenarios in China were simulated in the four periods. Moreover, a spatial shift trend model of mean center was introduced to quantitatively calculate the shift direction and distance of each ecotone type during the periods from T0 to T3. The simulated results show that there are 41 ecotone types in China, accounting for 18% of the total land area of China. The ecotones of cold temperate grassland / humid forest and warm temperate arid forest (564238.5 km2), cold temperate humid forest and warm temperate arid / humid forest (566549.75 km2), and northern humid / humid forest and cold temperate humid forest (525750.25 km2) are the main ecotone types, accounting for 35% of the total area of ecotones in China. Between T0 and T3, the area in the ecotone of cold temperate desert shrub and warm temperate desert shrub / thorn steppe will increase at a rate of 4% per decade, which is up by 3604.2, 10063.1 and 17242 km2 per decade under the RCP2.6, RCP4.5 and RCP8.5 scenarios, respectively. The cold ecotones will transform to the warm humid ecotones in the future. The average shift distance of mean center in the ecotone of north wet forest and cold temperate desert shrub / thorn grassland will be generally larger than that of other ecotones, whose mean center will move to the northeast, and the shift distance will be more than 150 km between T0 and T3. In addition, with a gradual increase of temperature and precipitation, the ecotones in northern China will show a shifting northward trend, while the ecotones in southern China will decrease gradually, and their mean center move to the high-altitude areas. The effects of climate change on ecotones will show an increasing trend in China, especially in the Qinghai-Tibet Plateau.
Keywords:ecotone;mean center;spatial identification method;scenario simulation;China


PDF (6829KB)元数据多维度评价相关文章导出EndNote|Ris|Bibtex收藏本文
本文引用格式
范泽孟. 中国生态过渡带分布的空间识别及情景模拟. 地理学报[J], 2021, 76(3): 626-644 doi:10.11821/dlxb202103010
FAN Zemeng. Spatial identification and scenario simulation of ecotone distribution in China. Acta Geographica Sinice[J], 2021, 76(3): 626-644 doi:10.11821/dlxb202103010


1 引言

自1905年生态过渡带的概念被Clements[1]提出后,其定义和内涵经历了很长时期的讨论,于1988年在巴黎召开的第七届环境问题科学委员会(SCOPE)上被明确定义为:在时间和空间上特征具有不确定性,且相互影响的2个相邻生态系统之间的生态地带[2]。而且会议认为生态过渡带的概念把生态系统界面理论以及非稳定的脆弱性特征结合起来,可以作为识别全球变化的基本指标,呼吁国际生态学界对生态过渡带开展研究[2]。随着生态过渡带概念定义的明确和研究热度的增加[3,4,5],国际上将生态过渡带的研究逐渐延伸到100 m~10 km的空间尺度上,认为生态过渡带在空间格局水平上,可定义为是某一空间范围内的优势植被类型之间,在气候变化和人类活动驱动下发生快速交替演变的生态过渡区域[6,7]

基于上述定义,可认为生态过渡带不仅具有时间和空间2个维度的不确定性特征[8],生态过渡带内还通常叠加有两种或两种以上的生态系统类型,而不同生态系统间容易相互渗透、相互联系和相互作用。尤其是在一定的空间格局和空间尺度上,在气候变化和人类活动驱动下,生态过渡带内的优势生态系统之间会发生快速交替演变[9],进而引起生态过渡带时空分布的演替变化[5, 10]。譬如,生态过渡带的植被类型分布变化及其对气候变化的响应速度均快于相邻的非生态过渡带区域[11,12];气候干旱和变暖对生态过渡带植物物种消长影响的风险要高于非过渡带区域[13,14,15];生态过渡带的土地覆盖变化强度及其景观多样性均高于非过渡带[16,17];降雨量对生态过渡带常绿植物冠层的影响会引起的植物组成结构变化幅度高于非过渡带[18];生态过渡带是植物生长及其光合作用对气候变暖响应敏感性较高的区域[19]

一方面,现有研究均表明生态过渡带是对气候变化和人类活动最为敏感的区域[20,21,22],不仅生态过渡带区域的生物物种易于扩散和传播,而且生态过渡带生物多样性更容易受到气候变化的影响而发生改变[23]。开展生态过渡带的生态系统时空变化研究比非过渡带更需要关注[24],这将有助于分析生态过渡带物种传播模式,有效构建物种类群优先保护框架[25],进而降低气候变化对生物多样性影响的风险系数[26]。另一方面,目前虽在针对某一类生态过渡带区域内的土地利用/覆盖变化[17, 27-30]、植被结构与分布[6]、生物多样性[31,32,33]、生态系统评估[34]、遥感监测[35]等方面开展系列研究,但这些研究主要是开展某一类型的生态过渡带区域内的生态系统对气候变化和人类活动响应规律,而在生态过渡带空间分布范围和边界的定量识别领域,则一直主要是处于生态过渡带边界确定的讨论及概念模型构建,两种生态系统类型过渡带以及某一区域的生态过渡带边界识别的研究阶段[36]。譬如,基于景观格局过渡与边界概念关系的生态过渡带层次概念模型[37];运用线性和“S”型回归方法构建植被—气候关系方程识别美国明尼苏达州草原和森林生态系统边界[38];基于Alpha和Gamma多样性分析澳大利亚荒漠沙丘区域的三叶草地和金合欢灌丛之间的过渡和突变边界[39];依据食草动物空间活动特征区分陆地与水生生态系统的过渡带边界[40]

鉴于此,如何构建能够适用于大尺度生态过渡带类型及其分布范围的空间识别方法,定量揭示生态过渡带类型空间分布范围、时空变化特征及其类型间差异,是厘清各类生态过渡带及其内部的各种生态系统类型之间对气候变化和人类活动响应的时空差异特征亟需研究的热点问题。因此,本文旨在对HLZ(Holdridge Life Zone)模型进行拓展,建立生态过渡带类型判别规则,进而构建生态过渡带类型的空间分析模型,实现全国生态过渡带类型及其分布范围的空间识别。并在此基础上,结合平均中心的时空偏移分析模型,运用全国气候观测数据和IPCC CMIP5的不同气候情景模式数据,对全国生态过渡带类型及空间分布及未来变化趋势的情景模拟。

2 数据与方法

2.1 数据收集与处理

用于生态过渡带类型空间识别及模拟分析的气候数据,主要包括观测数据和模式模拟的情景数据。其中,气候观测数据来源于全国752个气象观测站点的月观测数据(1981—2010年)。气候情景数据采用IPCC数据分布中心的CMIP5情景中代表温室气体排放的高中低3种情景[41],即:RCP 2.6(低排放情景)、RCP 4.5(中间排放情景)、RCP 8.5(高排放情景)(https://esgf-node.llnl.gov/projects/cmip5/)的2011—2040年(T1)、2041—2070年(T2)和2071—2100年(T3)3个时段的气候情景数据。全国1 km×1 km的DEM采用SRTM数据来自于美国航空航天局(NASA)的数据网站(http://srtm.csi.cgiar.org)。

如何对气候站点观测数据进行空间插值和对气候模式情景数据进行空间降尺度,进而获得高精度的气候要素空间分布数据,对生态过渡带空间分布及情景变化模拟结果具有重要影响,甚至直接关系到生态过渡带类型定量识别的准确度。空间插值和降尺度模型主要包括反距离加权模型(IDW)、三角网模型(TIN)、克里金模型(Kriging)和样条插值模型(Spline),以及基于微分几何学曲面理论发展的高精度曲面建模(HASM,High Accuracy Surface Modeling)方法[42]。HASM方法能够克服IDW、TIN、Kriging、Spline等方法存在的理论缺陷,大幅度提高气候数据的空间插值和降尺度精度[43,44]。因此,文中采用高精度曲面建模(HASM)方法,结合经纬度及高程数据,实现气候观测数据的空间插值和未来气候情景数据的降尺度处理,分别获得全国1 km×1 km空间分辨率的T0(1981—2010年)、T1(2011—2040年)、T3(2041—2070年)和T3(2071—2100年)4个时段的年平均生物温度、年降水量和潜在蒸散比率数据。

2.2 生态过渡带类型的空间识别模型

平均生物气温、降水及潜在蒸散比率等关键气候要素空间分布的连续性特征[45],决定两种甚至多种生态系统类型在地球表层空间上并不是独立分布,而是在这些生态系统的空间分布边界地带形成一定范围的重叠区域。在这些区域内,平均生物气温、降水及潜在蒸散比率等关键气候要素的时空变化将直接影响到各种生态系统类型边界的拓展和萎缩。这些多种生态系统类型重叠区域则称之为生态过渡带(Ecotone)。在模型构建的过程中,基于GIS空间分析方法,在对HLZ模型参数格式、运行模式进行改进的基础上[46,47,48,49],将其分类体系中平均生物温度、平均降水和潜在蒸散比率刻度线相交而成的正三角形定义为生态过渡带(原HLZ模型中相邻正六边形的相交地带),从而建立生态过渡带类型的理论判别体系(图1)。

图1

新窗口打开|下载原图ZIP|生成PPT
图1生态过渡带空间识别机理及模型框架

Fig. 1Mechanism and scheme for identifying ecotones



具体计算方法为:根据重新定义的生态过渡带判别体系,计算每一类生态过渡带的年均生物温度的边界阈值( MAB0i(℃))、年均降水量的边界阈值( TAP0i(mm))和潜在蒸散比率的边界阈值( PER0i),进而构建定量识别所有生态过渡带类型边界阈值的判定标准集合,整个集合中包括49类生态过渡带类型的判别标准值(表1)。譬如,如果某一生态过渡带的MABTAPPER 3个关键气候要素的值同时满足MAB > 0.375 ℃、TAP < 125 mm和PET < 0.25,那么该生态过渡带则被识别为风荒漠与冰雪过渡带。生态过渡带类型的空间定量识别理论计算公式可表达如下:

MABx,y,t=1365j=1365TEM(j,x,y,t)
TAPx,y,t=j=1365P(j,x,y,t)
PERx,y,t=58.93MAB(x,y,t)TAP(x,y,t)
ifMABx,y,tMAB0iTAPx,y,tTAP0iPERx,y,tPER0iEcotonex,y,t=i1,2,3,,49;elseEcotonex,y,t=0
式中, MABx,y,tTAPx,y,tPERx,y,t分别为t时段每一个栅格 x,y的年平均生物温度、年平均降水量和潜在蒸散比率; Ecotonex,y,t表示实现分类后的每种生态过渡带的类型值; MAB0iTAP0iPER0i分别代表第i种生态过渡带类型的判别标准值。

Tab. 1
表1
表1生态过渡带类型的判别标准
Tab. 1Identification criteria of ecotones
编码BTZ类型MAB(℃)TAP(mm)PER
1风荒漠与冰雪过渡带> 0.375< 125< 0.25
2风荒漠—冰缘与冰雪过渡带< 0.75> 125> 0.25
3风荒漠—冻荒漠与冰缘过渡带> 0.75< 125< 0.50
4冰缘与冰雪过渡带> 0.75< 250< 0.25
5冻荒漠—冰缘与高寒草原过渡带< 1.50> 125> 0.50
6冰缘—冰雪与高寒草甸过渡带< 1.50> 250> 0.25
7冻荒漠—高寒荒漠/草原过渡带> 1.50< 125< 1.00
8冰缘—高寒草原/草甸过渡带> 1.50< 250< 0.50
9冰雪—高寒草甸与苔原过渡带> 1.50< 500< 0.25
10高寒荒漠/草原与干旱灌丛过渡带< 3.00> 125> 1.00
11高寒草原/草甸与北方湿润森林过渡带< 3.00> 250> 0.50
12高寒草甸/苔原与北方潮湿森林过渡带< 3.00> 500> 0.25
13高寒荒漠与北方荒漠/干旱灌丛过渡带> 3.00< 125< 2.00
14高寒草原与北方干旱灌丛/湿润森林过渡带> 3.00< 250< 1.00
15高寒草甸与北方湿润/潮湿森林过渡带> 3.00< 500< 0.50
16北方潮湿/雨林与苔原过渡带> 3.00< 1000< 0.25
17北方荒漠/灌丛与冷温带荒漠灌丛过渡带< 6.00> 125> 2.00
18北方干旱灌丛/湿润森林与冷温带草原过渡带< 6.00> 250> 1.00
19北方湿润/潮湿森林与冷温带湿润森林过渡带< 6.00> 500> 0.50
20北方潮湿/雨林与冷温带潮湿森林过渡带< 6.00> 1000> 0.25
21北方荒漠与冷温带荒漠/灌丛过渡带> 6.00< 125< 4.00
22北方干旱灌丛与冷温带荒漠灌丛/草原过渡带> 6.00< 250< 2.00
23北方湿润森林与冷温带草原/湿润森林过渡带> 6.00< 500< 1.00
24北方潮湿森林与冷温带湿润/潮湿森林过渡带> 6.00< 1000< 0.50
25北方雨林与冷温带潮湿森林/雨林过渡带> 6.00< 2000< 0.25
26冷温带荒漠/灌丛与暖温带荒漠灌丛过渡带< 12.00> 125> 4.00
27冷温带荒漠灌丛/草原与暖温带有刺草原过渡带< 12.00> 250> 2.00
28冷温带草原/湿润森林与暖温带干旱森林过渡带< 12.00> 500> 1.00
29冷温带湿润/潮湿森林与暖温带湿润森林过渡带< 12.00> 1000> 0.50
30冷温带潮湿森林/雨林与暖温带潮湿森林过渡带< 12.00> 2000> 0.25
31冷温带荒漠与暖温带荒漠/灌丛过渡带> 12.00< 125< 8.00
32冷温带荒漠灌丛与暖温带荒漠灌丛/有刺草原过渡带> 12.00< 250< 4.00
33冷温带草原与暖温带有刺草原/干旱森林过渡带> 12.00< 500< 2.00
34冷温带湿润森林与暖温带干旱/湿润森林过渡带> 12.00< 1000< 1.00
35冷温带潮湿森林与暖温带湿润/潮湿森林过渡带> 12.00< 2000< 0.50
36冷温带雨林与暖温带潮湿森林/雨林过渡带> 12.00< 4000< 0.25
37亚热带荒漠/灌丛与热带荒漠灌丛过渡带< 24.00> 125> 8.00
38亚热带灌丛/有刺草原与热带有刺疏林过渡带< 24.00> 250> 4.00
39亚热带有刺疏林/干旱森林与热带很干森林过渡带< 24.00> 500> 2.00
40亚热带干旱/湿润森林与热带干旱森林过渡带< 24.00> 1000> 1.00
41亚热带湿润/潮湿森林与热带湿润森林过渡带< 24.00> 2000> 0.50
42亚热带潮湿森林/雨林与热带潮湿森林过渡带< 24.00> 4000> 0.25
43亚热带荒漠与热带荒漠/灌丛过渡带> 24.00< 125< 16.00
44亚热带荒漠灌丛与热带荒漠灌丛/有刺疏林过渡带> 24.00< 250< 8.00
45亚热带有刺疏林与热带有刺疏林/很干森林过渡带> 24.00< 500< 4.00
46亚热带干旱森林与热带很干/干旱森林过渡带> 24.00< 1000< 2.00
47亚热带湿润森林与热带干旱/湿润森林过渡带> 24.00< 2000< 1.00
48亚热带潮湿森林与热带湿润/潮湿森林过渡带> 24.00< 4000< 0.50
49亚热带雨林与热带潮湿/雨林过渡带> 24.00< 8000< 0.25

新窗口打开|下载CSV

在对生态过渡带类型进行空识别的过程中,当空间栅格单元 x,y处的 MABx,y,tTAPx,y,tPERx,y,t不满足任何一类生态过渡带的判别标准时,将空间栅格单元 x,y处的生态过渡带类型值赋为0,表示栅格单元 x,y不属于生态过渡带;如果栅格单元 x,y处的 MABx,y,tTAPx,y,tPERx,y,t满足49种生态过渡带类型中任何一类的判别标准时,则将空间栅格单元处的 x,y生态过渡带类型值赋为相应的生态过渡带类型编码。在该模型及所有的定量判别过程,均在栅格层次上利用程序设计和算法编程实现,运用识别标准和判别规则,将其与生态过渡带类型的识别标准集合进行对比,从而对每一个栅格单元进行判断识别并赋值,直到完成所有栅格单元的判断识别为止。最终获得研究区所有生态过渡带类型及其空间分布格局的定量识别。

2.3 生态过渡带平均中心偏移的时空分析模型

在计算生态过渡带空间分布平均中心的过程中,引入中心模型并赋予每个输入因子相应的含义[23,24,25],构建了生态过渡带平均中心偏移的时空分析模型,其理论计算公式可表达为:

xjt=i=1IjtsijtXijtSjt
yjt=i=1IjtsijtYijtSjt
式中:t为时间变量;Ij(t)为 t时段内第 j种生态过渡带类型的斑块数;sij(t)为 t时段内第 j种生态过渡带类型的第 i个斑块的面积;sj(t)为t时段内第 j种生态过渡带类型每个斑块的面积;(Xij(t), Yij(t))为 t时段内第 j种生态过渡带类型的第 i个斑块的几何中心的经纬度坐标;(Xj(t), Yj(t))为 t时段内第 j种生态过渡带类型空间分布的平均中心的经纬度坐标。

j种生态过渡带平均中心的偏移距离和方向可分别为:

dj=xj(t+1)-xj(t)2+yj(t+1)-yj(t)2
θj=arctgyj(t+1)-yj(t)xj(t+1)-xj(t)
其中,dj为第 j种生态过渡带平均中心从 tt+1时段的偏移距离;(xj(t), yj(t))和(xj(t+1), yj(t+1))分别代表第 j种生态过渡带 t时段和 t+1时段的平均中心坐标;θj为第j种生态过渡带平均中心从 tt+1时段的移动方向,当345°< θj ≤ 15°、75°< θj ≤ 105°、165°< θj ≤195°和225°< θj ≤ 285°时,可近似认为分别代表第 j种生态过渡带的平均中心从 t时段到 t+1时段分别向东、北、西和南偏移,而当5°< θj ≤ 75°、105°< θj ≤ 165°、195°< θj ≤ 255°和285°< θj ≤ 345°时,则分别代表第 j种生态过渡带的平均中心从 t时段到 t+1时段向东北、西北、西南和东南方向偏移。

3 结果分析

3.1 生态过渡带空间分布特征

IPCC RCP 2.6、RCP 4.5、RCP 8.5气候情景驱动下的T0~T3时段的生态过渡带未来情景模拟结果显示,全国总共将出现41种生态过渡带类型(图2~图4)。其中,东部地区的生态过渡带类型呈东北向西南方向的条带状分布,西部地区的生态过渡带类型分布零散且空间异质性高,生态过渡带区域主要分布在山区,而平原、盆地内部分布较少。在所有生态过渡带类型中,冷温带草原/湿润森林与暖温带干旱森林过渡带、冷温带湿润森林与暖温带干旱/湿润森林过渡带、北方湿润/潮湿森林与冷温带湿润森林过渡带是主要的过渡带类型,3种生态过渡带类型面积占所有过渡带面积的1/3。冷温带草原/湿润森林与暖温带干旱森林过渡带主要分布在长白山、吕梁山和太行山南部地区以及渭河北岸山地丘陵地区。冷温带湿润森林与暖温带干旱/湿润森林过渡带主要分布在长江中下游沿岸山区及云贵高原西部。北方湿润/潮湿森林与冷温带湿润森林过渡带主要分布在大兴安岭北部、横断山地区及四川盆地西北部的山区。风荒漠与冰雪过渡带、风荒漠—冰缘与冰雪过渡带、风荒漠—冻荒漠与冰缘过渡带、冰缘与冰雪过渡带、冻荒漠—冰缘与高寒草原过渡带、冻荒漠—高寒荒漠/草原过渡带、高寒荒漠/草原与干旱灌丛过渡带、高寒荒漠与北方荒漠/干旱灌丛过渡带、北方潮湿/雨林与苔原过渡带、北方荒漠/灌丛与冷温带荒漠灌丛过渡带、北方潮湿/雨林与冷温带潮湿森林过渡带、北方雨林与冷温带潮湿森林/雨林过渡带、冷温带潮湿森林/雨林与暖温带潮湿森林过渡带、亚热带荒漠/灌丛与热带荒漠灌丛过渡带、亚热带湿润/潮湿森林与热带湿润森林过渡带的面积之和仅占全国生态过渡带面积的2%左右,均主要分布在在青藏高原地区。

图2

新窗口打开|下载原图ZIP|生成PPT
图2基于RCP 2.6情景的中国生态过渡带空间分布

注:基于自然资源部标准地图服务网站GS(2016)2923号标准地图制作,底图无修改,图3~图5同。
Fig. 2Ecotones distribution in China under RCP2.6 scenario



图3

新窗口打开|下载原图ZIP|生成PPT
图3基于RCP 4.5情景的中国生态过渡带空间分布

Fig. 3Ecotones distribution in China under RCP 4.5 scenario



图4

新窗口打开|下载原图ZIP|生成PPT
图4基于RCP 8.5情景的中国生态过渡带空间分布

Fig. 4Ecotones distribution in China under RCP 8.5 scenario



3.2 生态过渡带分布面积及格局变化

对RCP 2.6、RCP 4.5、RCP 8.5气候情景驱动下中国生态过渡带的模拟结果(图2~图4)进行统计分析表明(表2),在2011—2100年间不同气候情景下的全国生态过渡带将呈现出系列变化趋势。

Tab. 2
表2
表2不同情景下的生态过渡带面积(km2)
Tab. 2Area of ecotones under the three scenarios of RCP 2.6, RCP 4.5 and RCP 8.5 (km2)
编码T0RCP 2.6RCP 4.5RCP 8.5
T1T2T3T1T2T3T1T2T3
1018112833201619124
2033131918635923460358448
301126144998232106366
4413325231929225222491362125523111251177
56917331621106620431487444193723521
6672071016554175630996827052324908523281271
704965698549450338462330
8217101183385751043611382630659471098661761778
95094261529541995764256865398313505058659301233972
1062268108756046027505765135346888189951
114281372193626996318277481518973368674401298729700
1251665607987472676468562648119491293578368651164255
1336982938256613013330281089130975689
1433304439703535028835445762691222438396532619420468
1553039104570121523125735105814129326145756109206147273124791
1620394101614141271415414
171448216502134371379417491119401106715377100013115
1831289451094636643487399865005253911402085719832904
1952387135375135381167223152001164383160644151677152096182743
2055370000212101514
211319422047180901488720233977957391724445322311
2263353421354234647085451894469046499474164557928390
231484361089811027199681310481591829778951063368694858513
2421484200644952412190922125104170215809054
25186000000000
261300801451201364291437181437691057731021511509629261023052
2749453107082126616131613116791144363139982120822144384119965
28206067247976236233227930215464268572267315215277230023200860
2937587590878222878204138630711333151645
3017433232391134130
311749110764911933710509992396114786145654104675152720101876
323540210993259939582234607021610417126351115379175955
339832519043298743235827524503844688428507102668167583
34221686226182188416212098239735163580154365225159132460132566
35152146769387840466482209713626260129116
37000000000884
3813000000000
3962000000000
4065438400143949436285482988900410961766368237548515481
412516636433491005323045937029917498
4718348396975569754925414248586213461042064157828317434

新窗口打开|下载CSV

在RCP 2.6情景下的T0~T3时段内,北方潮湿/雨林与冷温带潮湿森林过渡带、北方雨林与冷温带潮湿森林/雨林过渡带、亚热带荒漠/灌丛与热带荒漠灌丛过渡带、亚热带灌丛/有刺草原与热带有刺疏林过渡带、亚热带有刺疏林/干旱森林与热带很干森林过渡带将有可能消失。T0~T3时段内,北方湿润/潮湿森林与冷温带湿润森林过渡带、冷温带荒漠与暖温带荒漠/灌丛过渡带、冷温带荒漠灌丛/草原与暖温带有刺草原等过渡带的面积增幅最大,将分别增加11.4万km2、8.7万km2和8.2万km2;冷温带草原与暖温带有刺草原/干旱森林过渡带、冰缘—冰雪与高寒草甸过渡带、北方湿润森林与冷温带草原/湿润森林过渡带的面积减少最多,将分别平均减少6.5万km2、6.1万km2和5.1万km2;冻荒漠—冰缘与高寒草原过渡带、冷温带荒漠灌丛与暖温带荒漠灌丛/有刺草原过渡带增长速度最快,分别以144%/10a和101%/10a的速度增长;北方潮湿/雨林与苔原过渡带、亚热带湿润/潮湿森林与热带湿润森林过渡带减少速度最快,分别以9.9%/10a和9.6%/10a的速度减少。

在RCP 4.5情景下的T0~T3时段内,北方雨林与冷温带潮湿森林/雨林过渡带、亚热带荒漠/灌丛与热带荒漠灌丛过渡带、亚热带灌丛/有刺草原与热带有刺疏林过渡带、亚热带有刺疏林/干旱森林与热带很干森林过渡带将可能完全消失。冷温带荒漠与暖温带荒漠/灌丛过渡带、亚热带湿润森林与热带干旱/湿润森林过渡带、北方湿润/潮湿森林与冷温带湿润森林过渡带相均呈快速增长趋势,在T0~T3时段内将分别增长12.8万km2、11.6万km2和10.8万km2。北方湿润森林与冷温带草原/湿润森林过渡带、冷温带湿润森林与暖温带干旱/湿润森林过渡带、冰缘—冰雪与高寒草甸过渡带将呈严重缩减趋势,在T0~T3时段内将分别减少7万km2、6.7万km2和6.4万km2。冷温带荒漠灌丛与暖温带荒漠灌丛/有刺草原过渡带、冷温带荒漠与暖温带荒漠/灌丛过渡带增长速度最快,在T0~T3时段内将分别以284%/10a和73%/10a的速度增长,而北方潮湿/雨林与苔原过渡带、北方潮湿/雨林与冷温带潮湿森林过渡带则均以9.9%/10a的速度减少。

在RCP 8.5情景下的T0~T3时段内,北方雨林与冷温带潮湿森林/雨林过渡带、亚热带灌丛/有刺草原与热带有刺疏林过渡带、亚热带有刺疏林/干旱森林与热带很干森林过渡带将可能消失。亚热带干旱/湿润森林与热带干旱森林过渡带、亚热带湿润森林与热带干旱/湿润森林过渡带、冷温带荒漠灌丛与暖温带荒漠灌丛/有刺草原过渡带呈快速扩张趋势,T0~T3时段分别扩张45万km2、29.9万km2和17.2万km2。冷温带荒漠/灌丛与暖温带荒漠灌丛过渡带、北方湿润森林与冷温带草原/湿润森林过渡带、冷温带湿润森林与暖温带干旱/湿润森林过渡带缩减面积最大,T0~T3时段内分别减少10.7万km2、9万km2和8.9万km2。冷温带荒漠灌丛与暖温带荒漠灌丛/有刺草原过渡带、亚热带湿润森林与热带干旱/湿润森林过渡带增长速度最快,T0~T3将分别以每10年487%和163%的速度增长。冷温带潮湿森林与暖温带湿润/潮湿森林过渡带、高寒荒漠与北方荒漠/干旱灌丛过渡带、北方潮湿/雨林与冷温带潮湿森林过渡带的缩减速度最快,T0~T3时段内将分别以每10年9.9%的速度减少。

3.3 生态过渡带平均中心的时空偏移

运用生态过渡带平均中心时空偏移模型对基于RCP 2.6、RCP 4.5和RCP 8.5三种情景的生态过渡带平均中心的时空偏移进行求算(表3~5,表中的空白均表示某一种生态过渡带类型在相邻两个时段内的某一个时段没有出现),并对平均中心偏移距离大于50 km的生态过渡带类型的平均中心进行时空变化趋势制图(图4),发现RCP 2.6情景驱动下的生态过渡带平均中心偏移距离大于50 km的生态过渡带有15种类型,RCP 4.5情景驱动下的生态过渡带有27种类型,RCP 8.5情景驱动下的生态过渡带有31种类型。

Tab. 3
表3
表3RCP 2.6情景下生态过渡带平均中心的时空偏移(km)
Tab. 3Shift trends of mean center in ecotones under RCP 2.6 scenario (km)
编码T0~T1T1~T2T2~T3
偏移距离偏移方向偏移距离偏移方向偏移距离偏移方向
1575.89东南158.44西
236.46364.35东南
3524.49西北440.09东南
4313.37西45.65西北33.72西
5163.57东北204.43西北80.98东南
675.23111.08西北34.18西南
762.35183.79西北
8344.76西9.82东北34.16
965.0166.15西3.16东北
10114.76西北102.64西北80.54西北
11133.9990.37西北25.27
12265.35东北75.64西10.24
13221.6西33.91西北23.62西
14226.46西南72.05西北153.18西北
15171.44东南110.42西36.23东南
162827.93东北251.73西南251.45东北
17289.66西南134.86西218.66西
18447.7564.87西南127.8
19119.67东南150.92西南197.31
2186.68东南56.18西南97.37西南
22981.11西南167.08西南30.74东北
23533.87西南66.24西309.09东北
242268.19东北843.93西南1032.88东北
2654.23西南18.27西9.99东南
27101.5东北38.75西北167.87
28216.51东北101.74东北72.23东北
29160.22西20.78西46.98
30766.46东南13.120
31259.19东南95.16东南44.37西北
32839.21415.56东南112.89西北
33294.23西南48.46西北21.06东南
34220.07西南52.519.88西北
35378.4835.47东南9.45东南
40488.65西北48.12东北381.42
41641.45东北216.49东北169.64东北
4729.03西北61.22.36

新窗口打开|下载CSV

Tab. 4
表4
表4RCP 4.5情景下生态过渡带平均中心的时空偏移(km)
Tab. 4Shift trends of mean center in ecotones under RCP 4.5 scenario (km)
编码T0~T1T1~T2T2~T3
偏移距离偏移方向偏移距离偏移方向偏移距离偏移方向
1293.08123.63
236.08东北279.36东南
3473.97西北494.86东南
4320.13西52.11西北57.99西
5170.54203.39西北70.2东南
687.43237.81西161.8西北
7120.74西北104.33
8321.52西28.5112.16西北
964.65115.89西北71.13西北
10170.76西北119.2西北33.71西
11131.38114.57西北101.63西北
12268.45东北149.39西95.55西
13226.55西28.1354.75西
14217.42西南248.67西175.08西北
15213.48东南120.03西北119.29西
162827.64东北2225.08西南60.36
17384.24西372.57西72.11西南
18325.97西南234.8252.36
19174.34432.43西南308.11西南
200.28东北
2169.1377.76西南412.93西南
22919.58西南332.35西南150.7西南
23445.81西南260.71西南155.94西南
242650.17东北179.94西南1561.57西南
2663.85西82.43西57.92西南
27306.92146.99西北140.7西北
28320.71东北191.07东北24.79西北
29153.5西北214.07西103.68西
30768.54东南19.0824.07
31235.48东南346.3767.52东南
32743.09664.54东南20.01
33274.56西南99.58西北89.56西北
34322.29西南137.42东北76.99东北
35384.46105.23东南94.5东南
40556.84西253.32东北108.94东北
411017.64东北58.91东南7.72东北
4734.71西北136.54东北74.47

新窗口打开|下载CSV

Tab. 5
表5
表5RCP 8.5情景下生态过渡带平均中心的时空偏移(km)
Tab. 5Shift trends of mean center in ecotones under RCP 8.5 scenario (km)
编码T0~T1T1~T2T2~T3
偏移距离偏移方向偏移距离偏移方向偏移距离偏移方向
1122.27东南333.54西
2930.5东南285.74东南
3297.21东南81.38
4326.75西136.14西北64.26
5173.46246.25东南840.32东南
686.9445.58西北133.04西北
7295.09东南
8333.53西147.15西北78.14西
967.49185.05西北150.35西北
10189.19西北76.72西74.34
11137.95239.1西北265.65西北
12280.61东北257.8西290.64西
13224.04西79.49西南50.82西北
14259.6西南386.82西北113.5西北
15196.98东南202.81西北299.93西北
162827.64东北2282.75西南206.86
17479.18西346.49西南49.96西
18326.78西南349.86428.34西南
19231.84727.43西南656.35西
201.12
21100.831140.32西南89.01西北
22956.6西南508.54西南719.47西南
23389.8西南791.86西南1331.83西南
242649.15东北744.07西南1838.76西南
2664.76西206.18西427.11西
27289.31东北278.89西北378.4西
28364.48东北312.33东北176.63
29153.47西北267.16西138.46
30767.04东南39.75
31260.29东南388.81189.04
32720.01631东南250.86西北
33280.88西南215.3西北137.05
34349.58西南165.45953.92东北
35401.18232.65东南202.82东南
40510.16西418.01东北128.65东北
4113.76西北
4736.96西北233.08东北102.5

新窗口打开|下载CSV

其中,① 高寒荒漠/草原与干旱灌丛过渡带、高寒草原/草甸与北方湿润森林过渡带、高寒草原与北方干旱灌丛/湿润森林过渡带等生态过渡带的平均中心整体上呈向西北偏移趋势。揭示出随着平均气温的逐渐上升,随着主要位于青藏高原和天山区域的冰川消融的增加,相应的高寒生态过渡带将逐渐演替成其邻域内的非过渡带类型。② 北方干旱灌丛/湿润森林与冷温带草原过渡带、北方湿润/潮湿森林与冷温带湿润森林过渡带、北方荒漠与冷温带荒漠/灌丛过渡带、北方干旱灌丛与冷温带荒漠灌丛/草原过渡带、北方湿润森林与冷温带草原/湿润森林过渡带等生态过渡带的平均中心整体上呈向西南偏移趋势。③ 冰雪、高寒草甸与苔原过渡带、亚热带湿润森林与热带干旱/湿润森林过渡带等生态系统类型的平均中心整体上呈向北偏移趋势。

3.4 3种情景下的生态过渡带时空变化对比分析

RCP 2.6、RCP 4.5和RCP 8.5的3种气候变化情景下的全国生态过渡带时空分布及变化趋势模拟表明:随着平均生物温度的不断上升和和降水的不断增加,全国生态过渡带类型空间格局及其平均中心在2011—2100年间,将会发生一系列的时空变化趋势,且不同的气候变化情景对生态过渡带类型的变化将产生不同的驱动效应。其中,① 冰缘与冰雪过渡带、冰缘—冰雪与高寒草甸过渡带、冰缘—高寒草原/草甸过渡带、冷温带潮湿森林与暖温带湿润/潮湿森林过渡带的面积,在RCP 4.5和RCP 8.5气候变化情景下均将呈持续减少趋势,而亚热带湿润森林与热带干旱/湿润森林过渡带的面积将呈持续增加趋势;② 北方湿润/潮湿森林与冷温带湿润森林过渡带的面积在RCP 2.6和RCP 8.5气候变化情景下将持呈续增加趋势,而RCP 4.5情景下也仅是在T2~T3时段有呈减少趋势;③ 高寒草甸/苔原与北方潮湿森林过渡带、高寒草甸与北方湿润/潮湿森林过渡带的面积在RCP 2.6和RCP 4.5气候变化情景下将呈持续增加趋势,而RCP 8.5情景下也仅是在T2~T3时段有呈减少趋势;④ 2种情景下的高寒荒漠与北方荒漠/干旱灌丛过渡带、北方湿润森林与冷温带草原/湿润森林过渡带的面积均将呈持续减少趋势而冷温带荒漠灌丛与暖温带荒漠灌丛/有刺草原过渡带的面积将呈持续增加趋势。

总之,在RCP 2.6、RCP 4.5和RCP 8.5的3种气候变化驱动的未来情景下,高寒冰雪等高寒过渡带类型约占总生态过渡带总面积的20%。虽有个别高寒生态过渡带的面积呈增长趋势,如高寒草甸/苔原与北方潮湿森林过渡带、高寒草甸与北方湿润/潮湿森林过渡带,但大部分高寒生态过渡带呈减少趋势,甚至部分过渡带面临消失,如风荒漠与冰雪过渡带、风荒漠、冻荒漠与冰缘过渡带、冻荒漠、高寒荒漠/草原过渡带的面积均不超过50 km2。其中,冻荒漠、高寒荒漠/草原过渡带将在RCP 8.5情景的T3时段消失。这一研究显示,未来的全球气候变暖,将会使冰雪消融增加,导致土壤水分增加,进而引起冻荒漠、高寒荒漠/草原过渡带朝其邻域的非过渡带类型发生演替。

图5

新窗口打开|下载原图ZIP|生成PPT
图5不同情景下的中国生态过渡带平均中心偏移趋势

Fig. 5Shift trends of mean center in ecotones under the three scenarios of RCP 2.6, RCP 4.5 and RCP 8.5



另外,在3种气候变化情景的驱动下,高寒草原与北方干旱灌丛/湿润森林过渡带、北方荒漠/灌丛与冷温带荒漠灌丛过渡带、北方干旱灌丛/湿润森林与冷温带草原过渡带、北方湿润/潮湿森林与冷温带湿润森林过渡带、北方干旱灌丛与冷温带荒漠灌丛/草原过渡带、北方湿润森林与冷温带草原/湿润森林过渡带,冷温带潮湿森林/雨林与暖温带潮湿森林过渡带、冷温带荒漠灌丛与暖温带荒漠灌丛/有刺草原过渡带、亚热带干旱/湿润森林与热带干旱森林过渡带等生态过渡带类型的平均中心时空偏移幅度整体上大于其它生态过渡带类型的平均中心偏移幅度。以上生态过渡带的平均中心在T0~T3时段内,每个时段的平均偏移幅度均超过150 km。这表明在同样的气候变化条件下,以上生态过渡带类型对气候相关因子的灵敏性高于其它的生态过渡带类型。

4 讨论

如何获取高精度的年平均生物温度、年平均降水和潜在蒸散比率等关键模型参数的空间分布数据,直接关系到生态过渡带空间识别模型模拟结果的质量[16]。文中结合经纬度及高程数据,采用HASM方法,对年平均生物温度、年平均降水和潜在蒸散比率三个关键模型参数的现状数据和未来情景进行空间插值和降尺度模拟,能够有效保证生态过渡带模型参数的精度和质量。

基于HLZ模型修正和拓展构建的生态过渡带类型的空间识别模型,能够通过年平均生物温度、年平均降水和潜在蒸散比率3个关键模型参数对生态过渡带类型及空间分布进行定量刻画和求算,有效弥补了现有研究仅通过景观格局过渡与边界概念关系确定生态过渡带层的次概念模型[37]、运用线性和“S”型回归方法构建植被—气候线性关系模型[38],以及仅局限于局地区域的边界识别或者是两种生态系统类型的突变边界识别[40]的研究缺陷。另外,生态过渡带空间识别模型能够结合不同的未来气候变化情景数据对生态过渡带类型未来变化的时空情景进行模拟分析。

运用RCP 2.6、RCP 4.5和RCP 8.5三种未来气候变化情景数据,2011—2100年间不同情景下全国生态过渡带的时空变化模拟显示,在未来80年内,随着平均气温的逐渐增加,高寒型生态过渡带的面积总体上呈现减少趋势,而且主要位于青藏高原的寒冷型生态过渡带的平均中心呈向西北偏移趋势。这一研究能够有效阐释随着气温的逐渐上升,青藏高原和天山区域的冰雪将呈逐渐向高海拔退缩的趋势,而且由于冰雪消融水导致相应区域的土壤湿度的增加,使原有的生态过渡带类型演替为邻域范围内的其他非过渡带生态类型[48,49]。另外,由于干旱区域的气温升高,蒸腾强度上升,致使干旱型的生态过渡带普遍呈向西南偏移趋势[50]。总体上看,寒冷型生态过渡带和干旱型过渡带在同样的气候变化驱动下,对气候相关因子的灵敏性普遍高于其他的生态过渡带类型。

另外,生态过渡带作为对气候变化和人类活动最为敏感的区域,区域内的植被结构及其分布、土地利用覆盖不仅在同等气候变化条件下发生的变化高于其他的非过渡带区域[20,21,22],而且在人类活动干扰下也较非过渡带容易发生变化[16]。全国的农牧交错带大部分位于生态过渡带区域,其耕地和草地之间的变化强度和幅度要高于其他的区域。尤其是处于生态过渡带的农田控制水源涵养区域,在农牧生产、林果种植及生态旅游等人类活动的作用下,植被结构和土地覆盖格局的变化要高于其他的区域[51]。该模型方法主要阐释了在不同气候变化情景驱动下的全国生态过渡带类型空间分布的时空变化趋势,有效反应了气候变化对不同生态过渡带类型时空演替的驱动效应,揭示不同生态过渡带类型对气候变化驱动效应的敏感性问题。但目前主要侧重于气候变化对生态过渡带的驱动效应分析。在后续的研究过程中,将深入探讨人类活动作用对生态过渡带时空格局变化的驱动效应,综合分析自然气候要素与人类活动双重耦合驱动下的全国生态过渡带时空变化规律,并对其变化幅度和强度进行归因分析。

5 结论

基于HLZ模型修正和拓展构建的生态过渡带类型空间识别模型,能够有效实现全国生态过渡带类型及空间分布格局的定量识别。引入平均中心空间分析模型,能够实现全国生态过渡带平均中心的偏移方向和距离的时空刻画和表征。冷温带草原/湿润森林与暖温带干旱森林过渡带、冷温带湿润森林与暖温带干旱/湿润森林过渡带、北方湿润/潮湿森林与冷温带湿润森林过渡带是全国分布最广的三种生态过渡带类型,在三种气候变化情景下的平均面积分别为564238.5 km2、566549.75 km2和525750.25 km2,其总面积约占全国生态过渡带总面积的1/3。在RCP 2.6、RCP 4.5和RCP 8.5未来气候变化情景的驱动下,冷温带荒漠灌丛与暖温带荒漠灌丛/有刺草原过渡带的面积增加最快,将分别增加3604.2 km2/10a、10063.1 km2/10a和17242 km2/10a。高寒荒漠与北方荒漠/干旱灌丛过渡带、冰缘—冰雪与高寒草甸过渡带、冷温带潮湿森林与暖温带湿润/潮湿森林过渡带等生态过渡带类型,在未来气候变化驱动下的面积减少总体上高于其他的生态过渡带类型。

参考文献 原文顺序
文献年度倒序
文中引用次数倒序
被引期刊影响因子

Clements F C. Research Methods in Ecology. Lincoln: University Publishing Corporation, 1905.
[本文引用: 1]

Holland M M. SCOPE/MAB technical consultations on landscape boundaries: Report on a ACOPE /MAB Workshop on Ecotones
Biology International, 1988,17:47-106.

[本文引用: 2]

Castri F, Hansen A J, Holland M M. A new look at ecotones: Emerging international projects on landscape boundaries
Biology International, 1988,17(Special Issue):1-167.

[本文引用: 1]

Niu Wenyuan. The discriminatory index with regard to the weakness, overlapness, and breadth of ecotone
Acta Ecological Sinica, 1989,9(2):97-105.

[本文引用: 1]

[ 牛文元. 生态环境脆弱带ECOTONE的基础判定
生态学报, 1989,9(2):97-105.]

[本文引用: 1]

Ma Shijun. A Perspective of Modern Ecology
Beijing:Science Press, 1990.

[本文引用: 2]

[ 马世骏. 现代生态学透视
北京:科学出版社, 1990.]

[本文引用: 2]

Bestelmeyer B T, Wiens J A. Local and regional-scale responses of ant diversity to a semiarid biome transition zone
Ecography, 2001,24(4):381-392.

DOI:10.1111/eco.2001.24.issue-4URL [本文引用: 2]

Gosz J R. Ecotone hierarchies
Ecological Applications, 1993,3(3):369-376.

DOI:10.2307/1941905URLPMID:27759251 [本文引用: 1]
Ecological phenomena are evident over a broad spectrum of space and time scales. Ecotones, being defined as zones of transition between adjacent ecological systems, also must occur over an equally broad spectrum of space and time scales. Scale-dependent constraints influence ecological phenomena and resulting zones of transition; however, in traditional ecotone studies little treatment has been given to these influences. This paper addresses aspects of the ecotone concept that relate to the strength of interactions between adjacent systems for a hierarchy of ecotones in a biome transition area in central New Mexico on the Sevilleta National Wildlife Refuge. Zones of transition occur at plant, population, patch, landscape, and biome levels in the hierarchy suggested here. Constraints differ across this hierarchy, primarily because of the different scales at which these constraints exert their influences. The basic strategy to understand these cross-scale influences must be to perform studies at different scales, and a hierarchical approach identifies those scales. This also is important for identifying the appropriate technologies that focus at the scales where transition zones between ecological systems/phenomena are expressed. A broad array of technologies are available for integrating the pattern-process relationships that occur across the many scales in ecological systems.

Ferro I, Morrone J J. Biogeographical transition zones: A search for conceptual synthesis
Biological Journal of the Linnean Society, 2014,113(1):1-12.

DOI:10.1111/bij.12333URL [本文引用: 1]
We revise concepts, definitions and examples of biogeographical transition zones to help develop a conceptual framework and differentiate them from other transitions that occur in the geographical space. A biogeographical transition zone is defined as a geographical area of overlap, with a gradient of replacement and partial segregation between biotic components (sets of taxa that share a similar geographical distribution as a product of a common history). It is an area where physical features, environmental conditions and ecological factors allow the mixture and the co-occurrence of two or more biotic components, but also constrain their distribution further into one another. The biogeographical affinities of the taxa are the most fundamental information to consider to analyse biogeographical transition zones accurately. By plotting the frequency of different distribution patterns on maps, gradual changes in their relative contribution to a given area can be perceived. Thus, the most heterogeneous places in terms of distributional patterns can be defined numerically on strictly geographical grounds. Biogeographical transition zones can occur at every hierarchical level of a biogeographical regionalization as long as different biotic components come into contact geographically. Ecological boundaries or ecotones when characterized only by differences in dominance of some species, life forms or sets of characteristic species (not endemic) are not biogeographical transition zones as defined herein. (C) 2014 The Linnean Society of London,

Wuyts B, Champneys A R, House J I. Amazonian forest-savanna biostability and human impact
Nature Communication, 2017,8(1):15519. DOI: 10.1038/ncomms15519

DOI:10.1038/ncomms15519URL [本文引用: 1]

Peteet D. Sensitivity and rapidity of vegetational response to abrupt climate change
PNAS, 2000,97(4):1359-1361.

DOI:10.1073/pnas.97.4.1359URLPMID:10677467 [本文引用: 1]
Rapid climate change characterizes numerous terrestrial sediment records during and since the last glaciation. Vegetational response is best expressed in terrestrial records near ecotones, where sensitivity to climate change is greatest, and response times are as short as decades.

Breshears D D, Cobb N S, Rich P M, et al. Regional vegetation die-off in response to global-change-type drought
PNAS, 2005,102(42):15144-15148.

URLPMID:16217022 [本文引用: 1]

Williamsa A P, Allen C D, Millar C I, et al. Forest responses to increasing aridity and warmth in the southwestern United States
PNAS, 2010,107(50):21289-21294.

DOI:10.1073/pnas.0914211107URLPMID:21149715 [本文引用: 1]
In recent decades, intense droughts, insect outbreaks, and wildfires have led to decreasing tree growth and increasing mortality in many temperate forests. We compared annual tree-ring width data from 1,097 populations in the coterminous United States to climate data and evaluated site-specific tree responses to climate variations throughout the 20th century. For each population, we developed a climate-driven growth equation by using climate records to predict annual ring widths. Forests within the southwestern United States appear particularly sensitive to drought and warmth. We input 21st century climate projections to the equations to predict growth responses. Our results suggest that if temperature and aridity rise as they are projected to, southwestern trees will experience substantially reduced growth during this century. As tree growth declines, mortality rates may increase at many sites. Increases in wildfires and bark-beetle outbreaks in the most recent decade are likely related to extreme drought and high temperatures during this period. Using satellite imagery and aerial survey data, we conservatively calculate that approximately 2.7% of southwestern forest and woodland area experienced substantial mortality due to wildfires from 1984 to 2006, and approximately 7.6% experienced mortality associated with bark beetles from 1997 to 2008. We estimate that up to approximately 18% of southwestern forest area (excluding woodlands) experienced mortality due to bark beetles or wildfire during this period. Expected climatic changes will alter future forest productivity, disturbance regimes, and species ranges throughout the Southwest. Emerging knowledge of these impending transitions informs efforts to adaptively manage southwestern forests.

Dolezal J, Dvorsky M, Kopecky M, et al. Vegetation dynamics at the upper elevational limit of vascular plants in Himalaya
Scientific Reports, 2016,6:24881. DOI: 10.1038/srep24881.

DOI:10.1038/srep24881URLPMID:27143226 [本文引用: 1]
A rapid warming in Himalayas is predicted to increase plant upper distributional limits, vegetation cover and abundance of species adapted to warmer climate. We explored these predictions in NW Himalayas, by revisiting uppermost plant populations after ten years (2003-2013), detailed monitoring of vegetation changes in permanent plots (2009-2012), and age analysis of plants growing from 5500 to 6150 m. Plant traits and microclimate variables were recorded to explain observed vegetation changes. The elevation limits of several species shifted up to 6150 m, about 150 vertical meters above the limit of continuous plant distribution. The plant age analysis corroborated the hypothesis of warming-driven uphill migration. However, the impact of warming interacts with increasing precipitation and physical disturbance. The extreme summer snowfall event in 2010 is likely responsible for substantial decrease in plant cover in both alpine and subnival vegetation and compositional shift towards species preferring wetter habitats. Simultaneous increase in summer temperature and precipitation caused rapid snow melt and, coupled with frequent night frosts, generated multiple freeze-thaw cycles detrimental to subnival plants. Our results suggest that plant species responses to ongoing climate change will not be unidirectional upward range shifts but rather multi-dimensional, species-specific and spatially variable.

Rich R L, Stefanski A, Montgomery R A, et al. Design and performance of combined infrared canopy and belowground warming in the B4WarmED (Boreal Forest Warming at an Ecotone in Danger) experiment
Global Change Biology. 2015,21(6):2334-2348.

DOI:10.1111/gcb.12855URLPMID:25640748 [本文引用: 1]

Feddema J J, Oleson K W, Bonan G B, et al. The importance of land-cover change in simulating future climates
Science, 2005,310(5754):1674-1678.

DOI:10.1126/science.1118160URLPMID:16339443 [本文引用: 1]
Adding the effects of changes in land cover to the A2 and B1 transient climate simulations described in the Special Report on Emissions Scenarios (SRES) by the Intergovernmental Panel on Climate Change leads to significantly different regional climates in 2100 as compared with climates resulting from atmospheric SRES forcings alone. Agricultural expansion in the A2 scenario results in significant additional warming over the Amazon and cooling of the upper air column and nearby oceans. These and other influences on the Hadley and monsoon circulations affect extratropical climates. Agricultural expansion in the mid-latitudes produces cooling and decreases in the mean daily temperature range over many areas. The A2 scenario results in more significant change, often of opposite sign, than does the B1 scenario.

Fan Z M, Li J, Yue T X. Land-cover changes of biome transition zones in Loess Plateau of China
Ecological Modelling, 2013,252:129-140.

DOI:10.1016/j.ecolmodel.2012.07.039URL [本文引用: 3]
The Holdridge life zone (HLZ) model has been improved to help classify the biome transition zone (BTZ) in China's Loess Plateau. A positive and negative transformation index of land-cover (PNTIL) was developed to quantitatively evaluate the land-cover changes in every type of BTZ. Three bioclimatic datasets, with a spatial resolution of 1 km x 1 km, were used to classify the BTZ type in Loess Plateau. These include the mean annual biotemperature (MAB), average total annual precipitation (TAP) and potential evapotranspiration ratio (PER). In 1985, 1995 and 2005 land cover data was used to analyze the changes within BTZs. The results show that there are 14 BTZ types, which account for 25.21% of the total land-cover area in Loess Plateau. From 1985 to 2005, cultivated land decreased 0.93% per decade; on average wetland and water areas, woodland and grassland increased 3.47%, 0.24% and 0.06% respectively per decade. During this period the total rate of whole BTZ land-cover transformation decreased from 28.53% to 21.91%. Overall the total positive and negative transformed areas of land cover in BTZs displayed a decreasing trend. Moreover, the results indicate that the transition zones may have exhibited a greater change and landscape diversity than the adjacent biomes in Loess Plateau from 1985 to 2005. (C) 2012 Elsevier B.V.

Guan K Y, Pan M, Li H B, et al. Photosynthetic seasonality of global tropical forests constrained by hydroclimate
Nature Geoscience, 2015,8(4):284-289.

DOI:10.1038/ngeo2382URL [本文引用: 2]

Reich P B, Sendall K M, Rice K, et al. Geographic range predicts photosynthetic and growth response to warming in co-occurring tree species
Nature Climate Change, 2015,5(2):148-152.

DOI:10.1038/nclimate2497URL [本文引用: 1]

Gao Hongwen. Advancement of theoretical research in ecotone
Chinese Journal of Ecology, 1994,13(1):32-38.

[本文引用: 1]

[ 高洪文. 生态交错带(Ecotone)理论研究进展
生态学杂志, 1994,13(1):32-38.]

[本文引用: 1]

Smith T B, Wayne R K, Girman D J, et al. A role for ecotones in generating rainforest biodiversity
Science, 1997,276(20):1855-1857.

DOI:10.1126/science.276.5320.1855URL [本文引用: 2]

Mayle F E, Burbridge R, Killeen T J. Millennial-scale dynamics of southern Amazonian rain forests
Science, 2000,290(22):2291-2294.

DOI:10.1126/science.290.5500.2291URL [本文引用: 2]

Chapin II I F S, Zavaleta E S, Eviner V T, et al. Consequences of changing biodiversity
Nature, 2000,405(6783):234-242.

DOI:10.1038/35012241URLPMID:10821284 [本文引用: 2]
Human alteration of the global environment has triggered the sixth major extinction event in the history of life and caused widespread changes in the global distribution of organisms. These changes in biodiversity alter ecosystem processes and change the resilience of ecosystems to environmental change. This has profound consequences for services that humans derive from ecosystems. The large ecological and societal consequences of changing biodiversity should be minimized to preserve options for future solutions to global environmental problems.

Pauli H, Gottfried M, Dullinger S, et al. Recent plant diversity changes on Europe's mountain summits
Science, 2012,336(20):353-355.

DOI:10.1126/science.1219033URL [本文引用: 2]

Marshall C J, Liebherr J K. Cladistic biogeography of the Mexican transition zone
Journal of Biogeography, 2000,27(1):203-216.

DOI:10.1046/j.1365-2699.2000.00388.xURL [本文引用: 2]

Loehle C. Disequilibrium and relaxation times for species responses to climate change
Ecological Modelling, 2018,384:23-29.

DOI:10.1016/j.ecolmodel.2018.06.004URL [本文引用: 2]

Wang X Y, Li Y, Chen Y P, et al. Temporal and spatial variation of extreme temperatures in an agropastoral ecotone of northern China from 1960 to 2016
Scientific Reports, 2018,8:8787. DOI: 10.1038/s41598-018-27066-0.

DOI:10.1038/s41598-018-27066-0URLPMID:29884865 [本文引用: 1]
The agro-pastoral ecotone of northern China is one of the areas most sensitive to global temperature change. To analyze the temporal and spatial trends of extreme temperature events in this area, we calculated the values of 16 extreme-temperature indices from 1960 to 2016 based on data from 45 national meteorological stations. We found that the coldest-temperature indices decreased significantly and the warmest-temperature indices increased significantly. The warming of night temperatures contributed more than warming of day temperatures to the overall warming trend. In addition, the warm-temperature indices appeared to be increasing since the late 1980s and early 1990s. Overall, though the four extremal indices showed an increasing trend, the rate of change in the minimum temperature was greater than that of the maximum temperature; thus, the minimum temperature contributed most strongly to the overall temperature increases. The growing season is being prolonged in higher-elevation areas, but vegetation maturation in lower-elevation areas has been accelerated by the high temperatures, potentially leading to a shorter growing season at low altitudes. However, the impacts of land-use changes caused by human activities on the temperature increases will require additional study.

Li S J, Sun Z G, Tan M H. Changing patterns in farming-pastoral ecotones in China between 1990 and 2010
Ecological Indicators, 2018,89:110-117.

DOI:10.1016/j.ecolind.2018.01.067URL [本文引用: 1]

Fan Z M, Zhang X, Li J, et al. Land-cover changes of national nature reserves in China. Journal of
Geographical Sciences. 2013,23(2):258-270.



Fan Z M, Li J, Yue T X, et al. Scenarios of land cover in karst area of southwestern China
Environmental Earth Sciences, 2015,74(8):6407-6420.

DOI:10.1007/s12665-015-4223-zURL

Rich R L, Stefanski A, Montgomery R A, et al. Design and performance of combined infrared canopy and belowground warming in the B4WarmED (Boreal Forest Warming at an Ecotone in Danger) experiment
Global Change Biology, 2015,21(6):2334-2348.

DOI:10.1111/gcb.12855URLPMID:25640748 [本文引用: 1]

Fan Zemeng, Zhang Xuan, Li Jing, et al. Transition trends of land-cover in national nature reserves of China
Acta Geographica Sinica, 2012,67(12):1623-1633.

DOI:10.11821/xb201212004URL [本文引用: 1]
A large number of nature reserves have been established in China aiming to prevent ecosystem degradation, protect natural habitats and conserve the biodiversity within the habitats. By the end of 2010, a total of 2588 nature reserves has been established in China and their total area was 149.44 million hectares, covering over 15% of China's total land area. As the primary driver of biodiversity change, land-cover change has direct effect on ecosystem structures and functions. Thus, a quantitative analysis of changes in the land-cover of nature reserves is a critical step for evaluating the effectiveness and improving the management policies of nature reserves. In terms of the ecosystem characteristics and its major protected objects, 180 National Nature Reserves (NNRs) are chosen and classified into 7 types in this paper. A Land-cover Transform Direction Index (LTDI) is developed on the basis of the contribution of each land-cover type to maintaining the ecosystem stability. In Northeast China, North China, East China, South China, Central-southern China, Northwest and Southwest China, LTDI is used to calculate the transition trend of land-cover in the core zone, buffer zone and experimental zone of each NNRs type during the period from the late 1980s to 2005. The results show that the mean transition rate of all selected NNRs types has become lower during the two decades. The land-cover transform rate of Southwest China was the largest, while that of East China was the smallest among the six regions. The mean positive and negative transform rates of land-cover in all core zones decreased by 0.69% and 0.16% respectively. The landscape pattern of land-cover in the core zones was more stable than that in the buffer zones and the experimental zones. The land-cover transformed rate of NNRs was less than that of Non-NNRs in general. Furthermore, the ecological diversity and patch connectivity of land-cover in the whole selected area increased generally during the period 1995-2005. In summary, the land-cover of NNRs in China has a beneficial change trend after the NNRs were established, especially during the period from 1995 to 2005.
[ 范泽孟, 张轩, 李婧, . 国家级自然保护区土地覆盖类型转换趋势
地理学报, 2012,67(12):1623-1633.]

[本文引用: 1]

Fan Z M, Bai R Y, Yue T X. Spatio-temporal distribution of vascular plant species abundance on Qinghai-Tibet Plateau
Journal of Geographical Sciences, 2019,29(10):1625-1636.

DOI:10.1007/s11442-019-1667-1URL [本文引用: 1]

Caneva G, Bartoli F, Savo V, et al. Combining statistical tools and ecological assessments in the study of biodeterioration patterns of stone temples in Angkor (Cambodia)
Scientific Reports, 2016,6:32601. DOI: 10.1038/srep32601.

DOI:10.1038/srep32601URLPMID:27597658 [本文引用: 1]
Biodeterioration is a major problem for the conservation of cultural heritage materials. We provide a new and original approach to analyzing changes in patterns of colonization (Biodeterioration patterns, BPs) by biological agents responsible for the deterioration of outdoor stone materials. Here we analyzed BPs of four Khmer temples in Angkor (Cambodia) exposed to variable environmental conditions, using qualitative ecological assessments and statistical approaches. The statistical analyses supported the findings obtained with the qualitative approach. Both approaches provided additional information not otherwise available using one single method. Our results indicate that studies on biodeterioration can benefit from integrating diverse methods so that conservation efforts might become more precise and effective.

Xiao Qianguang, Chen Weiying, Du Peng, et al. Monitoring the ecological transect in East Asia monsoon region by meteorological satellite remote sensing
Acta Botanica Sinica, 1997,39(9):826-830.

[本文引用: 1]

[ 肖乾广, 陈维英, 杜鹏, . 用气象卫星对东亚季风区的生态过渡带的遥感监测研究
植物生态学报, 1997,39(9):826-830.]

[本文引用: 1]

van Vuuren D P, Edmonds J, Kainuma M, et al. The representative concentration pathways: An overview
Climatic Change, 2011,109(1/2):5-31.

DOI:10.1007/s10584-011-0148-zURL [本文引用: 1]

Yarrow M M, Salthe S N. Ecological boundaries in the context of hierarchy theory
BioSystems, 2008,92(3):233-244.

DOI:10.1016/j.biosystems.2008.03.001URLPMID:18448239 [本文引用: 1]
Ecological boundaries have been described as being multiscalar or hierarchical entities. However, the concept of the ecological boundary has not been explicitly examined in the context of hierarchy theory. We explore how ecological boundaries might be envisioned as constituents of scalar hierarchical systems. Boundaries may be represented by the surfaces of constituents or as constituents themselves. Where surfaces would correspond to abrupt transition zones, boundary systems might be quite varied depending on hierarchical context. We conclude that hierarchy theory is compatible with a functional vision of ecological boundaries where functions can be largely represented as the processing or filtering of ecological signals. Furthermore, we postulate that emergent ecological boundaries that arise on a new hierarchical level may contribute to the overconnectedness of mature ecosystems. Nevertheless, a thermodynamic approach to the emergence and development of boundary systems does indicate that in many situations, ecological boundaries would persist in time by contributing to the energy production of higher hierarchical levels.

Yarrow M M, Marín V H. Toward conceptual cohesiveness: A historical analysis of the theory and utility of ecological boundaries and transition zones
Ecosystems, 2007,10(3):462-476.

DOI:10.1007/s10021-007-9036-9URL [本文引用: 2]
Ecological transition zones are increasingly recognized as systems that play a critical role in controlling or modifying flows of organisms, materials, and energy across landscapes. Many concepts describing transitional areas have been proposed over the years, such as the prevalent and durable ecotone concept. Confusion among ecologists and land managers about transition zone concepts and the isolation of studies that use only one transition concept can hinder unified progress in understanding these key systems. Currently, a movement toward conceptual synthesis under the umbrella concept of ‘ecological boundary’ is underway. Here we examine the history and theoretical baggage of the ecotone, riparian zone, and several other concepts. Subsequently, we present a conceptual cluster analysis, which facilitates a better understanding of the similarities and differences between boundary and transition concepts. We emphasize the hierarchical nature of these concepts: higher-level synthetic concepts can be used in the development of theory, whereas lower-level concepts allow more specificity and the formulation of operational definitions. Finally, we look briefly at the utility and future use of boundary and transition zone concepts.]]>

Danz N P, Frelich L E, Reich P B, et al. Do vegetation boundaries display smooth or abrupt spatial transitions along environmental gradients? Evidence from the prairie-forest biome boundary of historic Minnesota, USA
Journal of Vegetation Science, 2013,24(6):1129-1140.

DOI:10.1111/jvs.12028URL [本文引用: 2]

Nicholas A M M, Franklin D C, Bowman D M J S. Floristic uniformity across abrupt boundaries between Triodia hummock grassland and Acacia shrubland on an Australian desert sandplain
Journal of Arid Environments, 2011,75(11):1090-1096.

DOI:10.1016/j.jaridenv.2011.06.016URL [本文引用: 1]

Sarneel J M, Huig N, Veen G F, et al. Herbivores enforce sharp boundaries between terrestrial and aquatic ecosystems
Ecosystems, 2014,17(8):1426-1438.

DOI:10.1007/s10021-014-9805-1URL [本文引用: 2]
The transitions between ecosystems (ecotones) are often biodiversity hotspots, but we know little about the forces that shape them. Today, often sharp boundaries with low diversity are found between terrestrial and aquatic ecosystems. This has been attributed to environmental factors that hamper succession. However, ecosystem properties are often controlled by both bottom-up and top-down forces, but their relative importance in shaping riparian boundaries is not known. We hypothesize that (1) herbivores may enforce sharp transitions between terrestrial and aquatic ecosystems by inhibiting emergent vegetation expansion and reducing the width of the transition zone and (2) the vegetation expansion, diversity, and species turnover are related to abiotic factors in the absence of herbivores, but not in their presence. We tested these hypotheses in 50 paired grazed and ungrazed plots spread over ten wetlands, during two years. Excluding grazers increased vegetation expansion, cover, biomass, and species richness. In ungrazed plots, vegetation cover was negatively related to water depth, whereas plant species richness was negatively related to the vegetation N:P ratio. The presence of (mainly aquatic) herbivores overruled the effect of water depth on vegetation cover increase but did not interact with vegetation N:P ratio. Increased local extinction in the presence of herbivores explained the negative effect of herbivores on species richness, as local colonization rates were unaffected by grazing. We conclude that (aquatic) herbivores can strongly inhibit expansion of the riparian vegetation and reduce vegetation diversity over a range of environmental conditions. Consequently, herbivores enforce sharp boundaries between terrestrial and aquatic ecosystems.

Yue T X. Surface Modelling:High Accuracy and High Speed Methods. New York: CRC Press, 2011.
[本文引用: 1]

Yue T X, Zhao N, Fan Z M, et al. CMIP5 downscaling and its uncertainty in China
Global and Planetary Change, 2016,146(1):30-37.

DOI:10.1016/j.gloplacha.2016.09.003URL [本文引用: 1]

Yue T X, Zhao N, Liu Y, et al. A fundamental theorem for eco-environmental surface modelling and its applications
Science China Earth Sciences, 2020,63(8):1092-1112.

DOI:10.1007/s11430-019-9594-3URL [本文引用: 1]

Yue T X, Fan Z M, Liu J Y. Changes of major terrestrial ecosystems in China since 1960
Global and Planetary Change, 2005,48(4):287-302.

DOI:10.1016/j.gloplacha.2005.03.001URL [本文引用: 1]

Biermann F. "Earth system governance" as a crosscutting theme of global change research
Global Environmental Change, 2007,17(3/4):326-337.

DOI:10.1016/j.gloenvcha.2006.11.010URL [本文引用: 1]

Mayle F E, Beerling D J, Gosling W D, et al. Responses of Amazonian ecosystems to climatic and atmospheric carbon dioxide changes since the last glacial maximum
Philosophical Transactions of the Royal Society B: Biological Sciences, 2004,359(1443):499-514.

DOI:10.1098/rstb.2003.1434URL [本文引用: 1]

Yue T X, Fan Z M, Liu J Y, et al. Scenarios of major terrestrial ecosystems in China
Ecological Modelling, 2006,199(3):363-376.

DOI:10.1016/j.ecolmodel.2006.05.026URL [本文引用: 1]

Fan Z M, Fan B, Yue T X. Terrestrial ecosystem scenarios and their response to climate change in
Eurasia. Science China-Earth Sciences, 2019,62(10):1607-1618.

[本文引用: 2]

Fan Z M, Bai R Y, Yue T X. Spatio-temporal distribution of vascular plant species abundance on Qinghai-Tibet Plateau
Journal of Geographical Sciences, 2019,29(10):1625-1636.

DOI:10.1007/s11442-019-1667-1URL [本文引用: 2]

Fan Z M, Fan B. Shifts of the mean center of potential vegetation ecosystems under future climate change in Eurasia
Forest, 2019,10(873). DOI: 10.3390/f10100873.

[本文引用: 1]

Shi Wenjiao, Liu Yiting, Shi Xiaoli. Quantitative methods for detecting the impacts of climate change on the fluctuation of farming-pastoral ecotone boundaries in northern China
Acta Geographica Sinica, 2017,72(3):407-419.

DOI:10.11821/dlxb201703004URL [本文引用: 1]
The quantitative analysis of the effect of climate change on the fluctuation of farming-pastoral ecotone (FPE) boundary in northern China is a current focus in the field of response to climate change in ecological vulnerable regions. Previous studies have given profound descriptions about the effects of climate change on the boundary shifts of the ecotones qualitatively, but lacked quantitative analysis of the contribution of climate change on both spatial and temporal scales. Here, climate data from national meteorological stations and land use data interpreted from remote sensing images of the farming-pastoral ecotone (FPE) in northern China since 1970 were used to describe boundaries of the FPE based on both climate and land use in the 1970s, 1980s, 1990s and 2000s. Detection in horizontal and vertical directions method (FishNet) and Digital Shoreline Analysis System method (DSAS) were applied to detect the spatial pattern of the FPE boundaries and to examine how much of the boundaries shifts can be explained by climate change in different periods. The results showed that the spatial pattern of the FPE boundaries and contributions of climate varied in different regions and periods. The FPE boundaries moved slightly in the northwest part of the FPE and violently in the northeast part. The shift of climate and land use boundaries in northwest segment of the Greater Hinggan Mountains showed the most highly coupling relationship and the contribution rate of climate change reached 10.7%-44.4% in east-west direction and 4.7%-55.9% in north-south direction based on FishNet method and 1.1%-16.8% based on DSAS method from the 1970s to the 2000s. Most of the detections based on the two methods had consistent results. Moreover, DSAS method was better than FishNet method in overall accuracy and suitable for the precise detection in small range, while FishNet method was more appropriate for intuitive, rapid and low-precision analysis. Our findings highlight the importance of different adaption measures to climate change in the FPE of northern China in different regions and periods.
[ 史文娇, 刘奕婷, 石晓丽. 气候变化对北方农牧交错带界线变迁影响的定量探测方法研究
地理学报, 2017,72(3):407-419.]

[本文引用: 1]

相关话题/生态 空间 过渡带 数据 生物