Forest biomass characteristics and influencing factors in different restoration stages in the Daxing’anling forest region of Inner Mongolia, northern China
-
摘要:目的针对内蒙古大兴安岭林区中火烧、采伐干扰后退化森林生态系统恢复演替进程,研究其林分生物量特征,确定林分、地形和气候因子对林分生物量的影响,建立林分生物量通用模型,进行森林碳汇功能研究。方法利用野外实测的胸径和树高推算生物量。采用方差分析确定林分、地形和气候因子以及恢复时间、干扰类型等对林分生物量有无影响,使用一般线性模型建立林分生物量的多元线性模型,比较退化森林生态系统在不同生境下林分生物量特征。结果(1)不同恢复阶段的林分生物量差异明显,森林各器官的生物量分配有所不同,生物量最大比重多出现在4 ~ 8 cm、14 cm和18 ~ 26 cm径阶。(2)方差分析结果表明,恢复时间、干扰类型、林分优势树种、海拔对林分生物量恢复有显著性影响,而气候因子或气候因子的2个主成分对林分生物量的影响都不显著,2个林分生物量模型的确定系数在0.85以上,可以用来预测退化森林生态系统中的森林生物量。结论不同恢复阶段森林生物量有较大差异,森林乔木各器官生物量分配、径阶分配受到干扰类型和恢复时间影响。林分生物量与时间和海拔成正比,干扰方式与林型会影响林分生物量,研究结果可为大兴安岭林区森林生物量预测以及森林碳库研究提供借鉴。Abstract:ObjectiveTo study the succession process of degraded forest ecosystems recovered from the interference of fire and logging in the Inner Mongolia Daxing’anling forest area of northern China, a general stand biomass model was established to analyze the forest biomass characteristics and figure out the effects of stand, topography and climate factors on stand biomass.MethodBiomass was calculated from the published biomass equations based on corresponding diameter at breast (DBH) and tree height (H) measured in the field. Variance analysis was used to determine the significance of stand, topographic and climatic factors, restoration time and interference type to the stand biomass. General linear model method was used to establish multivariate linear model of stand forest biomass and the forest biomass characteristics of degraded forest ecosystem were compared among different habitats.Result(1) The forest biomass in different recovery stages of forest ecosystems was significantly different. The forest biomass organic allocation was also significantly different in a way that the maximum proportion of biomass mostly occurred in 4−8 cm, 14 cm and 18−26 cm diameter class. (2) Variance analysis showed that restoration time, interference type, dominant tree species and altitude had significant effects on stand biomass, while the two main climatic factors and climatic factors had no effect on stand biomass. The coefficient of determination for two stand biomass models was above 0.85 and could be used to predict forest biomass in degraded forest ecosystem.ConclusionThere are large differences in forest biomass in varied degraded stages of forest ecosystems. The forest biomass organic and diameter allocation are affected by interference type and restoration time. The stand biomass is proportional to restoration time and altitude. The interference type and forest type will affect the stand biomass. The research results can provide reference for forest biomass prediction and forest carbon pool research in Daxing’anling forest area.
-
Keywords:
- burned land /
- cutover land /
- restoration time /
- biomass /
- influencing factor
-
我国城镇化进程带来经济增长与生活便利的同时,也给城市生态空间带来了挑战。一些大型城市与城市群面临生态资源保护不足、价值转化受阻、生态环境发展滞后,以及区域间生态服务功能差异等问题[1]。为此,国家提出了新型城镇化战略,旨在探索“城镇高质量发展”的城市转型发展新路径。在这个过程中,国土空间规划扮演着关键角色。随着 “五级三类”国土空间规划体系的逐步完善,跨行政边界、区域尺度的生态空间体系构建变得尤为重要。这不仅有助于实现“多规合一”,也是推进生态环境保护和新型城镇化发展的重要支撑。城市群与都市圈同时具有城市空间的结构特征和国土空间的尺度规模[2]。都市圈作为中国城镇化的重要空间单元,其空间结构和组织方式对整个国土空间的规划实践和政策制定具有重大影响,是探索和实现城镇高质量发展的重要一环[3]。为实现经济、社会和生态环境的整体优质发展,都市圈需要提升城市发展质量和环境宜居度。因此,有必要从系统耦合协调的视角,对都市圈尺度国土空间规划进行深入研究,提出有效的对策建议。以上海大都市圈为代表的现代化城市群,进行城镇化与生态环境耦合协调发展的研究,有助于探索都市圈跨区域尺度的高质量城市建设和生态空间体系构建的发展路径。
城镇化与生态环境协调优质发展的关键是促进两个系统达到各自功能和整体功能最优,即二者在时间、功能和发展速度上实现良性耦合与协同完善[4]。正视城镇化对生态安全的威胁与城市发展的失衡问题,学者们愈加关注城市生态环境对可持续发展的影响,城镇化与生态环境间的交互耦合研究成为学术前沿,主要包括以下3个方面。一是关注城镇化对生态环境的影响,即由城镇化造成的城市生态环境风险升级。如Alberti等[5]综述了不同城市模式产生不同生态效应的实证证据;岳邦瑞等[6]提出土地覆被变化影响生态系统服务的形成与供给。二是揭示生态环境对城市发展的制约作用,多侧重城市资源禀赋及生态环境对城镇化的响应研究。如徐海龙等[7]基于3种生态安全格局,进行了城市用地空间扩展的多情景模拟;Wu等[8]指出水资源短缺是制约部分高度城镇化城市的普遍因素。三是探索城镇化与生态环境间的耦合作用关系。理论研究方面, Kijima等[9]提出城市与生态环境间存在环境库兹涅茨EKC倒U型曲线及非线性关系。实证研究方面,Shafiei等[10]证实了城镇化和生态环境间存在U型曲线变化关系,邵明等[3]强调了自然和社会经济结构的相互作用对城市群发展的限制。在研究方法上,学者常采用多方法交叉叠合的方式,包括耦合协调模型[7]、系统动力学模型[11]和动态计量分析[12]等。在研究尺度上,涉及城市、城市群或流域等多尺度研究,同时中国国土空间规划逐渐重视跨行政边界的大尺度规划,更多学者从城市群等区域尺度开展研究。如针对珠三角地区[13]、长江流域地区[14]等多地域进行的耦合协调关系的实证研究。
综上所述,关于城镇化与生态环境关系的研究思路已由单向影响机制转为耦合协调的交互探索;研究方法由基于现状评价的静态定性分析转为多时空维度的定量动态评价;研究对象与尺度更注重多尺度、多要素及多时空维度的交互作用研究。但现有研究仍存在不足。首先,长三角城市群作为研究对象备受关注,而上海大都市圈作为其一体化发展格局的重要部分却受到忽视。其次,以都市圈为对象的研究中,城市尺度的研究较多,而区县尺度的考察分析缺乏。基于此,本文以上海大都市圈“1 + 8”城市群为研究区域,构建城镇化系统与生态环境系统综合评价指标体系,基于2011—2020年面板数据,聚焦城市和区县两个空间尺度,运用耦合协调度模型定量分析两系统时空演变特征,旨在揭示2011—2020年间上海大都市圈城镇化与生态环境协调发展的基本规律,以期为都市圈未来高质量健康发展提供参考。
1. 研究区概况与数据来源
1.1 研究对象
以上海大都市圈为研究对象,范围包括上海、无锡、常州、苏州、南通、宁波、湖州、嘉兴、舟山在内的“1 + 8”市域行政区域。该都市圈地处长三角核心区,陆域总面积达5.6万km2,以长三角1/6的陆域面积承载了长三角1/3的人口和约1/2的经济总量,呈现空间集约高效的城镇化空间形态。随着长三角一体化的国家战略地位提升,上海大都市圈作为承接我国经济高质量发展的最重要载体之一,是国家新型城镇化的主体区[15],对城镇化与生态环境的协同优质发展提出了更高的要求。
1.2 数据来源
研究基于上海大都市圈“1 + 8”城市群2011—2020年的面板数据,数据来源于《中国环境年鉴》、《中国城市统计年鉴》、《中国城市建设年鉴》及各城市及辖区的相关年鉴。鉴于数据的完整性和权威性,对个别缺失及异常数据采用插值法、类比法等方法进行数据补全与修正。
2. 研究方法
2.1 指标体系构建
研究以中国学术期刊出版总库(CNKI)为文献来源,检索了2011—2022年间的核心期刊和CSSCI文献。通过“关键词 = 耦合协调”AND“关键词 = 城镇/城市化”AND“关键词 = 生态环境”的检索式,共找到相关文献185篇,剔除综述和无数据文献后,最终得到163篇相关文献。
统计发现,城镇化指标系统通常涵盖“人口–社会–经济–空间”4个维度,包括人口城镇化、经济城镇化、空间城镇化和社会城镇化4方面;生态环境指标系统可以从“规模–压力–治理–效率”4个维度与城镇化指标系统对标,囊括生态规模、环境水平、环境治理、资源效率4个方面。考虑到各地区数据统计口径及类型的不同,基于数据的均衡性和可得性,兼顾指标体系数据量的平衡,研究最终在所梳理的92个指标中为两系统各选取了研究对象在10年内均覆盖的9项基础指标作为指标因子,以反映各研究单元城镇化与生态环境系统的综合发展水平,构建上海大都市圈城镇化和生态环境综合评价指标体系(表1)。其中,人口城镇化包括市域人口密度和每万人拥有登记失业人数,旨在反映人口聚集度和就业机会水平,而每万人拥有登记失业人数作为衡量人口流动和就业机会的指标,间接反映城市对劳动力的吸引容纳能力,以评估人口城镇化的程度;社会城镇化包括每万人拥有医院床位数、每万人拥有中小学生人数、科教支出占一般公共预算支出比重,旨在衡量地区公共服务状况和科教水平;经济城镇化包括人均固定资产投资额、人均地区生产总值、第三产业增加值占地区生产总值比重,旨在从经济可持续性角度衡量城市经济活力和产业结构的发展;空间城镇化包括人均拥有居住面积,旨在对城市的空间布局和资源利用进行评估。生态规模包括建成区绿化覆盖率、人均公园绿地面积,旨在反映城市绿化水平和生态规模,为城市生态平衡提供依据;环境水平包括二氧化硫年日均浓度、可吸入颗粒物年日均浓度、年空气污染指数达优良的百分率,旨在反映城市空气质量状况和主要污染物来源;环境治理包括人均生活垃圾清运量、人均道路清扫面积、单位市域面积环保投入资金,旨在衡量城市卫生状况、清洁度和环境管理水平;资源效率包括市域单位面积GDP,旨在衡量城市空间集约化利用与自然资源的经济潜力水平,以反映城市土地资源需求、生态足迹以及对自然资源的利用效率。指标设计采用比率或每万人为单位的标准化处理方式,将各地区的指标数据转换为相对统一的尺度,以消除指标上的绝对数值差异;指标权重采用熵值法和CRITIC权重法2种客观评价方法进行比较和验证后求取平均值确定综合权重,以全面反映各子系统指标数据的特征和规律。
表 1 上海大都市圈城镇化与生态环境系统指标体系及指标权重Table 1. Index system and index weights of urbanization and ecological environment system in Shanghai metropolitan area系统
System子系统
Subsystem具体指标
Specific index属性
Attribute熵值法
Entropy
methodCritic权重
Critic weight综合权重
Comprehe-nsive
weight参考来源
Reference城镇化系统
Urbanization
system人口城镇化
Population urbanization市域人口密度/(人·km−2)
Urban population density/(person·km−2)+ 0.306 9 0.121 7 0.214 3 [16] 每万人拥有登记失业人数
Registered unemployment per 10 000 people− 0.032 1 0.109 6 0.070 9 [17] 社会城镇化
Social urbanization每万人拥有医院床位数
Number of hospital beds per 10 000 people+ 0.144 7 0.087 2 0.116 0 [17] 每万人拥有中小学生人数
Number of primary and secondary school students
per 10 000 people+ 0.105 0 0.109 1 0.107 0 [16] 科教支出占一般公共预算支出比重
Proportion of expenditure on science and education in the general public budget/%+ 0.032 4 0.099 5 0.066 0 [18] 经济城镇化
Economic urbanization人均固定资产投资额/元
Per capita investment in fixed assets/CNY+ 0.119 9 0.080 7 0.100 3 [18] 人均地区生产总值/元
Per capita regional GDP/CNY+ 0.091 8 0.068 5 0.080 2 [16] 第三产业增加值占地区生产总值比重
Proportion of added value of the tertiary industry to regional GDP/%+ 0.104 6 0.1376 0.1210 [19] 空间城镇化
Spatial urbanization人均拥有居住面积 Per capita living space/m2 + 0.062 6 0.186 1 0.124 3 [18] 生态环境系统
Eco-environment system生态规模
Ecological scale建成区绿化覆盖率
Green coverage rate of built-up area/m2+ 0.044 0 0.176 8 0.110 4 [19] 人均公园绿地面积 Per capita green park area/m2 + 0.098 4 0.182 5 0.140 4 [19] 环境水平
Environmental level二氧化硫年日均浓度
Annual average daily concentration of sulfur
dioxide /(μg·m−3)− 0.013 5 0.080 1 0.077 3 [20] 可吸入颗粒物年日均浓度
Annual average daily concentration of inhalable
particulate matter/(μg·m−3)− 0.031 3 0.074 7 0.092 0 [20] 年空气污染指数达优良的百分率
Percentage of excellent annual air pollution index/%+ 0.035 0 0.091 7 0.207 7 [16] 环境治理
Environmental governance人均生活垃圾清运量
Per capita domestic garbage removal volume/t+ 0.055 2 0.099 4 0.077 3 [20] 人均道路清扫面积
Per capita road clearing area/m2+ 0.082 4 0.101 6 0.092 0 [21] 单位市域面积环保投入资金 /(万元·km−2)
Environmental protection investment funds per unit city area/(104 CNY·km−2)+ 0.315 2 0.100 2 0.207 7 [22] 效率水平
Efficiency level市域单位面积GDP/(108元·km−2)
GDP per unit area of a city/(108 CNY·km−2)+ 0.325 0 0.093 0 0.209 0 [18] 2.2 综合评价模型
基于综合计算的权重和无量纲化后的指标值,根据主成分分析原理建立城镇化和生态环境综合评价模型。
城镇化综合评价模型:Ui=∑nj=1(X′ijWij) (1) 生态环境综合评价模型:Ei=∑nj=1(X′ijWij) (2) 式中:Ui、Ei为第i个地区城镇化、生态环境综合指数,Wij为第i个地区第j项指标的权重;X′ij为第i地区第j项指标数值标准化后的值。
2.3 耦合协调度关系模型
耦合度作为物理学概念,描述两个或多个系统间关系的程度,即互相依赖的程度。耦合度值介于0到1之间,数值越大表示系统间更趋于良性共振耦合和有序发展状态,其具体模型如式(3);在城镇化与生态环境两系统分析时,模型可表示如式(4)所示。
Cn={u1⋅u2⋅⋅⋅un[∏ui+uj]}1/n (3) C={UiEi(Ui+Ei)2}1/2 (4) 式中:Cn为物理学中的耦合度,n为耦合度模型子系统个数;C为城镇化与生态环境之间的耦合度;u1,u2,···,un为各系统的综合指数;ui为各系统的序参量;uj为具体的评估指标;Ui为第i个地区的城镇化综合指数,Ei为第i个地区的生态环境综合指数。
耦合度虽可以反映两系统间的相互作用强度,但难以衡量二者的协调水平。因此在耦合度模型基础上,进一步构建耦合协调度(D)模型。
D=√CT,T=αUi+βEi (5) 式中:T为综合发展指数,α、β是城镇化和生态坏境在系统中的权重,由于二者同等重要,取α=β=0.5。
参考相关学者研究[13],根据耦合协调度D值大小,将耦合协调类型划分为3个类型和5个亚类型。根据Ui和Ei的差异,比较两系统间的相对发展状态,划分为3种耦合协调特征,由此共得到15种子类型(表2)。
表 2 城镇化与生态环境耦合协调类型划分标准Table 2. Classification criteria for the coupling and coordination of urbanization and ecological environment类型
Category耦合协调度
Coupling and
coordination
degree (D)亚类型
Subcategory
系统比较
System
comparison子类型
Subtype协调发展
Coordinated
development[0.8 ~ 1.0) 优质协调
High-quality coordination (Ⅰ)Ui−Ei>0.1 优质协调-城镇化受阻型
High-quality coordination-hindered urbanization typeⅠ-1 Ei−Ui>0.1 优质协调-生态环境受阻型
High-quality coordination-ecological environment obstruction typeⅠ-2 0⩽∣Ei−Ui∣⩽0.1 优质协调-同步发展型
High-quality coordination-synchronous development typeⅠ-3 [0.7 ~ 0.8) 优化协调
Good coordination (Ⅱ)Ui−Ei>0.1 优化协调-城镇化受阻型
Good coordination-hindered urbanization typeⅡ-1 Ei−Ui>0.1 优化协调-生态环境受阻型
Good coordination-ecological environment obstruction typeⅡ-2 0⩽∣Ei−Ui∣⩽0.1 优化协调-同步发展型
Good coordination-synchronous development typeⅡ-3 转型发展
Transformational development[0.5 ~ 0.7) 基本协调
Basic coordination (Ⅲ)Ui−Ei>0.1 基本协调-城镇化受阻型
Basic coordination-hindered urbanization typeⅢ-1 Ei−Ui>0.1 基本协调-生态环境受阻型
Basic coordination-ecological environment obstruction typeⅢ-2 0⩽∣Ei−Ui∣⩽0.1 基本协调-同步发展型
Basic coordination-synchronous development typeⅢ-3 不协调发展
Uncoordinated
development[0.3 ~ 0.5) 轻度失调
Mild imbalance (Ⅳ)Ui−Ei>0.1 轻度失调-城镇化滞后型
Mild imbalance-hindered urbanization typeⅣ-1 Ei−Ui>0.1 轻度失调-生态环境受阻型
Mild imbalance-ecological environment obstruction typeⅣ-2 0⩽∣Ei−Ui∣⩽0.1 轻度失调-同步发展型
Mild imbalance-synchronous development typeⅣ-3 (0.0 ~ 0.3) 严重失调
Severe imbalance (Ⅴ)Ui−Ei>0.1 严重失调-城镇化受阻型
Severe imbalance-hindered urbanization typeⅤ-1 Ei−Ui>0.1 严重失调-生态环境受阻型
Severe imbalance-ecological environment obstruction typeⅤ-2 0⩽∣Ei−Ui∣⩽0.1 严重失调-同步发展型
Severe imbalance-synchronous development typeⅤ-3 注:Ui为第i个地区的城镇化综合指数,Ei为第i个地区的生态环境综合指数。Notes: Ui, comprehensive index of urbanization in the ith region. Ei, comprehensive index of ecological environment in the ith region. 3. 结果与分析
3.1 城镇化综合水平时空演变特征
区域层面,上海大都市圈城市群城镇化质量整体呈上升趋势,城镇化综合指数均值从2011年0.313上升至2020年0.412(图1)。城市层面, 9市城镇化综合水平波动式上升,各城市间差异明显,在2020年呈现以上海市中心城区、苏州–南通、常州为核心的“多核”空间结构。上海市辖区层面,10年间16个辖区城镇化综合水平缓慢提高,中心城区和郊区分异明显,呈两极分化局面(图2)。其中,中心城区各辖区城镇化综合指数均值由0.449提升至0.562,城镇化水平较高;郊区各辖区城镇化综合指数均值由0.195提升至0.260,仍处于低水平状态。
2011—2020年间,上海都市圈内9市城镇化综合水平在各阶段发展水平分化较大,大体呈现:上海中心城区 > 都市圈内9市 > 上海市郊区。具体来看, 9市增长步调存在差异:2011年嘉兴(0.346)、南通(0.344)、无锡 (0.332) 、湖州(0.328) 4市的城镇化综合指数略高,位于上海大都市圈(0.313)之上;随着城市发展的加快,至2016年,苏州、常州等市城镇化水平明显提高,各市差异趋小 ;至2020年底,苏州(0.442)、无锡(0.426)、常州(0.450)、南通市(0.457)4市的城镇化综合指数已经超过了整个都市圈的平均水平 (0.412)。上海中心城区与郊区的差异显著: 2011—2016年二者增幅极缓,一是中心城区本身城镇化水平较高,增幅有限;二是上海郊区开发迟缓,致整体发展低下,其中,崇明区城镇化综合指数介于[0.157,0.163],城镇化水平最低,且2011—2014年间呈现下降趋势,主要原因是人均居住面积的减少。2016—2020年,随着增幅的逐渐提升,中心城区各辖区城镇化综合指数均值由0.486增至0.562。其中,黄浦区始终处于高值区间,普陀、静安由于政府科教支出资金比例的缩小在个别年份出现小幅降低,其他4区保持同步发展状态;郊区各辖区城镇化综合指数均值由0.217增至0.260,始终处于低城镇化水平状态。
3.2 生态环境综合水平时空演变特征
区域层面,上海大都市圈城市群生态环境水平整体呈上升趋势,生态环境综合指数均值从2011年0.256上升至2020年0.376(图3)。城市层面, 9市生态环境综合水平差异逐渐缩小,至2020年呈现各地多片区协同发展的格局(图4)。上海市辖区层面,中心城区和郊区生态环境综合水平大致相同,但郊区由于单位面积环保投入资金较中心城区偏低,与中心城区生态环境指数的差距逐年增加。
2011—2020 年间,上海都市圈内9 市生态环境综合水平差距逐步缩小。具体来看,2011年,舟山(0.342)、苏州(0.291)、无锡(0.288)、湖州(0.285) 4市的生态环境综合水平领先,高于都市圈平均水平(0.256),南通市(0.203)最低;至2016年,湖州(0.353)、舟山(0.370)两市生态环境综合水平居首位,南通市因人均公园绿地面积等指标的提高,生态环境综合水平显著大幅提升,位于都市圈平均水平(0.335)之上;至2020 年,各城市间差异进一步缩小,仅苏州(0.369)、宁波(0.355)、常州 (0.351)、嘉兴(0.324) 4市较都市圈平均水平(0.376)略低。上海市辖区层面,10年间各辖区生态环境综合水平波动较大:2011—2016年间,除崇明、嘉定、虹口、杨浦区保持稳步增长外,其余辖区在不同年份均出现短暂下降,其中,普陀和金山区由于空气质量指标的下降而变化最大;2016—2020年,徐汇、静安、虹口等区增长较为稳健,其余辖区在个别年份生态环境综合水平均出现下降,分析原因大都由于市域单位面积环保投入资金及人均道路清扫面积的变化引起。
3.3 耦合协调度时空演变特征
3.3.1 耦合协调等级演变特征分析
时间维度上,10年间,上海大都市圈城镇化与生态环境耦合协调度(D)不断提升,由 0.502提升至0.700,耦合协调类型由“转型发展状态”转为“协调发展状态” (图5)。2011—2016年间,城市层面,除南通和无锡市在2014年由于生态环境综合指数下降及城镇化综合指数提升,耦合程度小幅下降外,各市均呈稳步增长态势。其中,南通市增幅最大,D值由0.447增至0.605,由轻度失调转为基本协调。苏州、无锡、湖州、舟山市D值均高于9市平均水平,始终维持在基本协调阶段。至2016年,都市圈内9市已全部步入基本协调亚类型,耦合协调类型以“转型发展状态”为主。上海市辖区层面,各区耦合协调度波动式提升,崇明区波动最为明显,5年中D值呈先下降后提升的趋势,原因是2011—2014年间社会和空间城镇化指数的下降,导致城镇化发展相对滞后。至2016年底,嘉定、金山、松江、青浦、奉贤、崇明区处于轻度失调状态,其他辖区均为基本协调阶段;黄浦、长宁区最优,分别到达优质协调和优化协调阶段。2016—2020年间,城市层面宁波、嘉兴、湖州3市保持在基本协调阶段,处于转型发展状态,其余城市均步入协调发展状态。上海市内,除长宁区在2019年由于生态环境综合指数值下降导致D值同步下降外,其余各辖区均稳步增长。至2020年16个辖区全部脱离“不协调发展状态”,其中黄浦、长宁区于2020年达到优质协调阶段。
空间维度上,上海大都市圈城镇化与生态环境的耦合协调度呈现出以上海市区、苏州–无锡–南通、舟山3个片区为中心辐射同步发展的格局(图6)。上海郊区耦合协调度提升最为缓慢;苏州–无锡、舟山两个片区周边城市耦合协调度水平提升明显,各市间存在良性同步发展趋势。
3.3.2 耦合协调类型识别及演变路径分析
依据表2的耦合协调度等级与类型划分标准,按照城镇化与生态环境综合指数的差异,划分各地区耦合协调类型,得到城镇化与生态环境耦合协调类型演变图(图7)。总体来看,两系统的耦合协调相对发展状态呈现由“同步发展为主,兼城镇化受阻”向“同步发展为主,兼生态环境受阻” 转变的演变态势。
2011—2014年,大部分地区处于同步发展型(62.5%),25%的地区处于城镇化受阻型,仅有12.5%地区属于生态环境受阻类型。说明此期间大部分地区生态环境和城市发展的步调一致,但部分地区由于前期大规模的城镇化建设,对生态环境发展的关注不高,形成城镇化建设与生态环境发展的错位。2014—2017年,同步发展型(62.5%)仍占较大比重,城镇化受阻型(20.8%)占比小幅下降,生态环境受阻型(16.7%)占比略有上升,说明在此期间各地区对生态环境的建设已初见成效,城镇化发展受阻的问题日益解决。2017—2020年,同步发展型(54.2%)小幅下降,城镇化受阻型(20.8%)占比未变,生态环境受阻型(25%)进一步上升。可见随着各地区环境保护政策的落实,使生态环境建设的速度更快,但两系统的同步发展状态逐渐被打破,耦合协调类型开始转向生态环境受阻型(表3)。
表 3 2011—2020年上海大都市圈及上海市各区城镇化与生态环境耦合协调度测算结果与耦合协调子类型划分Table 3. Calculation results of the coupling coordination between urbanization and ecological environment in the Shanghai metropolitan area and various districts from 2011 to 2020 and division of coupling coordination subtypes年份
Year上海大
都市圈
Shanghai
metropolitan
area苏州市
Suzhou
City无锡市
Wuxi
City常州市
Changzhou
City南通市
Nantong
City宁波市
Ningbo
City嘉兴市
Jiaxing
City湖州市
Huzhou
City舟山市
Zhoushan
City上海市
Shanghai
City黄浦区
Huangpu
District徐汇区
Xuhui
District长宁区
Changning
District2011 0.502 0.540 0.570 0.462 0.447 0.454 0.469 0.563 0.567 0.443 0.743 0.501 0.604 Ⅲ-2 Ⅲ-3 Ⅲ-3 Ⅳ-3 Ⅳ-1 Ⅳ-3 Ⅳ-1 Ⅲ-3 Ⅲ-3 Ⅳ-2 Ⅱ-1 Ⅲ-1 Ⅲ-1 2012 0.536 0.577 0.577 0.478 0.525 0.502 0.507 0.581 0.596 0.486 0.816 0.533 0.603 Ⅲ-2 Ⅲ-3 Ⅲ-3 Ⅳ-3 Ⅲ-1 Ⅲ-3 Ⅲ-1 Ⅲ-3 Ⅲ-3 Ⅳ-2 Ⅰ-1 Ⅲ-1 Ⅲ-1 2013 0.560 0.615 0.585 0.509 0.595 0.515 0.533 0.587 0.608 0.498 0.845 0.536 0.606 Ⅲ-2 Ⅲ-3 Ⅲ-3 Ⅲ-3 Ⅲ-3 Ⅲ-3 Ⅲ-1 Ⅲ-3 Ⅲ-3 Ⅳ-2 Ⅰ-1 Ⅲ-1 Ⅲ-1 2014 0.583 0.676 0.581 0.563 0.587 0.563 0.562 0.598 0.602 0.522 0.863 0.559 0.646 Ⅲ-2 Ⅲ-3 Ⅲ-3 Ⅲ-1 Ⅲ-3 Ⅲ-3 Ⅲ-3 Ⅲ-3 Ⅲ-3 Ⅲ-2 Ⅰ-1 Ⅲ-1 Ⅲ-1 2015 0.601 0.683 0.613 0.593 0.605 0.583 0.565 0.619 0.622 0.533 0.866 0.577 0.661 Ⅲ-2 Ⅲ-3 Ⅲ-3 Ⅲ-1 Ⅲ-1 Ⅲ-3 Ⅲ-3 Ⅲ-3 Ⅲ-3 Ⅲ-2 Ⅰ-1 Ⅲ-1 Ⅲ-1 2016 0.631 0.675 0.633 0.639 0.670 0.610 0.597 0.637 0.649 0.571 0.924 0.616 0.701 Ⅲ-2 Ⅲ-3 Ⅲ-3 Ⅲ-3 Ⅲ-3 Ⅲ-3 Ⅲ-3 Ⅲ-3 Ⅲ-3 Ⅲ-2 Ⅰ-1 Ⅲ-1 Ⅱ-1 2017 0.642 0.693 0.652 0.668 0.675 0.623 0.579 0.647 0.654 0.587 0.942 0.649 0.75 Ⅲ-2 Ⅲ-3 Ⅲ-3 Ⅲ-1 Ⅲ-3 Ⅲ-3 Ⅲ-3 Ⅲ-3 Ⅲ-3 Ⅲ-2 Ⅰ-1 Ⅲ-1 Ⅱ-3 2018 0.662 0.711 0.663 0.656 0.695 0.642 0.602 0.671 0.695 0.624 0.944 0.675 0.819 Ⅲ-2 Ⅱ-3 Ⅲ-3 Ⅲ-1 Ⅲ-3 Ⅲ-3 Ⅲ-3 Ⅲ-3 Ⅲ-3 Ⅲ-2 Ⅰ-1 Ⅲ-1 Ⅰ-1 2019 0.685 0.716 0.712 0.67 0.721 0.668 0.613 0.684 0.712 0.667 0.982 0.705 0.857 Ⅲ-2 Ⅱ-1 Ⅱ-3 Ⅲ-1 Ⅱ-3 Ⅲ-3 Ⅲ-3 Ⅲ-3 Ⅱ-3 Ⅲ-2 Ⅰ-1 Ⅱ-1 Ⅰ-1 2020 0.700 0.717 0.726 0.705 0.730 0.675 0.632 0.694 0.727 0.687 0.995 0.735 0.812 Ⅱ-2 Ⅱ-3 Ⅱ-3 Ⅱ-3 Ⅱ-3 Ⅲ-3 Ⅲ-3 Ⅲ-3 Ⅱ-3 Ⅲ-2 Ⅰ-1 Ⅱ-1 Ⅰ-1 年份
Year静安区
Jing’an
District普陀区
Putuo
District虹口区
Hongkou
District杨浦区
Yangpu
District浦东新区
Pudong
New Area闵行区
Minhang
District宝山区
Baoshan
District嘉定区
Jiading
District金山区
Jinshan
District松江区
Songjiang
District青浦区
Qingpu
District奉贤区
Fengxian
District崇明区
Chongming
District2011 0.618 0.452 0.287 0.437 0.546 0.444 0.471 0.312 0.395 0.343 0.307 0.349 0.286 Ⅲ-1 Ⅳ-1 Ⅴ-1 Ⅳ-1 Ⅲ-3 Ⅳ-3 Ⅳ-3 Ⅳ-3 Ⅳ-3 Ⅳ-3 Ⅳ-3 Ⅳ-3 Ⅴ-3 2012 0.642 0.530 0.409 0.476 0.577 0.513 0.497 0.349 0.435 0.364 0.374 0.405 0.259 Ⅲ-1 Ⅲ-1 Ⅳ-1 Ⅳ-1 Ⅲ-3 Ⅲ-3 Ⅳ-3 Ⅳ-3 Ⅳ-3 Ⅳ-3 Ⅳ-3 Ⅳ-3 Ⅴ-3 2013 0.617 0.580 0.447 0.502 0.617 0.515 0.537 0.370 0.393 0.396 0.368 0.379 0.277 Ⅲ-1 Ⅲ-1 Ⅳ-1 Ⅲ-1 Ⅲ-3 Ⅲ-3 Ⅲ-3 Ⅳ-3 Ⅳ-3 Ⅳ-3 Ⅳ-3 Ⅳ-3 Ⅴ-2 2014 0.655 0.645 0.48 0.547 0.615 0.533 0.582 0.4 0.432 0.382 0.394 0.402 0.223 Ⅲ-1 Ⅲ-3 Ⅳ-1 Ⅲ-1 Ⅲ-3 Ⅲ-3 Ⅲ-3 Ⅳ-3 Ⅳ-3 Ⅳ-3 Ⅳ-3 Ⅳ-3 Ⅴ-2 2015 0.662 0.613 0.528 0.55 0.626 0.519 0.582 0.395 0.417 0.417 0.395 0.392 0.343 Ⅲ-1 Ⅲ-1 Ⅲ-1 Ⅲ-1 Ⅲ-3 Ⅲ-3 Ⅲ-3 Ⅳ-3 Ⅳ-3 Ⅳ-3 Ⅳ-3 Ⅳ-3 Ⅳ-2 2016 0.688 0.674 0.596 0.603 0.651 0.588 0.589 0.456 0.442 0.418 0.414 0.409 0.362 Ⅲ-1 Ⅲ-3 Ⅲ-1 Ⅲ-1 Ⅲ-3 Ⅲ-3 Ⅲ-3 Ⅳ-3 Ⅳ-2 Ⅳ-3 Ⅳ-3 Ⅳ-3 Ⅳ-2 2017 0.699 0.675 0.601 0.611 0.666 0.572 0.616 0.494 0.436 0.436 0.425 0.418 0.384 Ⅲ-1 Ⅲ-3 Ⅲ-1 Ⅲ-1 Ⅲ-3 Ⅲ-3 Ⅲ-2 Ⅳ-2 Ⅳ-2 Ⅳ-3 Ⅳ-3 Ⅳ-3 Ⅳ-2 2018 0.696 0.730 0.637 0.654 0.708 0.607 0.639 0.509 0.511 0.502 0.475 0.464 0.418 Ⅲ-1 Ⅱ-3 Ⅲ-1 Ⅲ-1 Ⅱ-2 Ⅲ-3 Ⅲ-2 Ⅲ-2 Ⅲ-3 Ⅲ-3 Ⅳ-3 Ⅳ-3 Ⅳ-2 2019 0.739 0.809 0.671 0.691 0.736 0.671 0.653 0.56 0.52 0.554 0.529 0.492 0.505 Ⅱ-1 Ⅰ-1 Ⅲ-1 Ⅲ-1 Ⅱ-2 Ⅲ-3 Ⅲ-3 Ⅲ-2 Ⅲ-2 Ⅲ-2 Ⅲ-2 Ⅳ-2 Ⅲ-2 2020 0.781 0.856 0.747 0.735 0.724 0.678 0.659 0.545 0.538 0.562 0.552 0.565 0.514 Ⅱ-1 Ⅰ-1 Ⅱ-1 Ⅱ-1 Ⅱ-3 Ⅲ-3 Ⅲ-3 Ⅲ-2 Ⅲ-2 Ⅲ-2 Ⅲ-2 Ⅲ-2 Ⅲ-2 4. 结论与讨论
4.1 结 论
本研究立足上海大都市圈的国土空间,基于文献爬梳,在城镇化系统和生态环境系统中各遴选9项基础指标作为指标因子,构建城镇化与生态环境综合评价体系,涵盖“人口–社会–经济–空间”和“规模–压力–治理–效率”8个维度。利用2011—2020年间上海大都市圈内9个城市和上海16区的面板数据,运用耦合协调度模型,从城市和上海市辖区两个层面,在时间和空间两个维度,分析对比上海大都市圈城镇化与生态环境耦合协调发展的时空分异特征,揭示上海大都市圈城镇化与生态环境耦合协调关系及其发展趋势。
(1)2011—2020年,上海大都市圈城镇化–生态环境耦合协调度稳步上升,协调发展等级呈现“基本协调—优化协调”的演进态势。各地区耦合协调发展等级同步增长,但发展水平与步调差异较大。城市层面,常州、南通、苏州、无锡和舟山市协调水平较其余4市更高,其中常州、南通两市的增速最快,嘉兴、湖州两市呈低水平缓慢发展。辖区层面,上海市部分中心城区起步较高且增速较快,如黄浦、长宁、普陀区分别在2012、2018、2019年率先进入“优质协调阶段”;松江、金山、嘉定区等郊区起步较低,在2018—2020年才相继脱离“失调阶段”。
(2)上海大都市圈城镇化–生态环境的动态耦合关系在空间分布上并不均衡,呈现出以上海中心城区、苏州–无锡–南通、舟山3个片区为中心辐射同步发展的格局特征。城市层面,嘉兴、宁波耦合协调度水平较低,南通、无锡、舟山等城市则处于较高水平。上海市辖区层面,各辖区之间差异较大,中心城区耦合协调度水平整体优于郊区,呈两极分化格局。
(3)2011—2020年,上海大都市圈城镇化与生态环境综合水平总体趋于向良性耦合协调发展的状态,但存在城镇化发展与生态环境建设错位的现象。通过两系统相对发展水平的比较,发现二者耦合协调特征的演变趋势呈现由 “同步发展为主,兼城镇化受阻”向“同步发展为主,兼生态环境受阻”转变。
4.2 讨 论
(1)都市圈的高质量发展要求城镇化与生态环境在多维度上达到良性共振的耦合发展状态,根据上海大都市圈内各地区城镇化与生态环境协调水平分析,应在严控生态环境保护与城市发展的底线约束的基础上,针对性地提出平衡生态环境保护与城镇化发展的国土空间规划策略:一是对于金山、奉贤、崇明区等生态环境受阻地区,应秉承生态产业化的国土空间规划策略,关注生态产品价值的稳态转换,发挥生态价值资源禀赋。二是对于徐汇、静安、虹口区等城镇化受阻地区,应秉承生态基础化的国土空间规划策略,改善地区发展战略取向,将生态环境保护置于基础性地位,培育绿色转型发展的新业态模式。三是对于无锡市、闵行、松江区等同步发展型地区,应推进生态现代化的国土空间规划策略,在“生态优先,绿色发展”的规划原则基础上进一步实现生态现代化。
(2)上海大都市圈作为长三角地区高质量一体化发展的功能载体,其发展模式具有跨尺度、多类型的国土空间规划特点,应立足国土空间尺度进行统筹规划,以系统观念提出跨区域尺度的高质量规划发展路径。一是增强上海中心城区等内部核心地区对整个都市圈范围内的辐射带动功能,打破行政边界与自然条件形成的城市壁垒,畅通内部高质量发展渠道,缩小发展差异,推动各地区高质量协同发展。二是深化区域一体的生态建设保护机制和城市协同发展机制。城镇化与生态环境耦合系统是人地关系地域系统的重要一环[23],人地系统的耦合协调性又是生态安全的重要表现形式[24]。作为实现生态安全的基本载体,都市圈尺度下的生态空间由城市绿地转变为山水林田湖草居的复杂生命共同体[25]。通过多层级—多尺度的生态空间体系构建,建立生态安全格局总体框架,统筹整合山水林田湖草等生态要素和重点生态空间,可以在推进生态环境的综合治理的基础上,同步促进经济社会的发展。
(3)城镇化与生态环境间存在着复杂的耦合协调机制,通过对上海大都市圈城镇化与生态环境交互耦合类型及演变路径的研究,探索两系统的耦合协调发展规律,可以为都市圈未来可持续协调发展提供参考。本研究以国土空间规划和新型城镇化建设为目标,通过耦合协调度模型对上海大都市圈进行城镇化系统与生态环境系统的耦合协调关系研究,不仅能够推进城市建设和生态环境保护的可持续发展,也能为上海大都市圈高质量发展的落地实施提供理论依据。此外,本研究存在一定局限性。对城镇化与生态环境两系统内部要素的交互耦合作用讨论较少,在今后的相关研究中可基于不同类别的要素与城镇化过程的关联度分析,对关键影响要素与驱动机制进一步探索。
-
图 2 不同恢复阶段的森林林分生物量分配格局
LF表示轻度火烧;MF表示中度火烧;SF表示重度火烧;C表示采伐。数字12、14、20、22、29、31表示恢复时间(年)。下同。LF, mild fire; MF, moderate fire; SF, severe fire; C, cutover. Number 12, 14, 20, 22, 29, 31 represent 12, 14, 20, 22, 29, 31 years restoration time, respectively. The same below.
Figure 2. Stand biomass distribution pattern of different restoration stages
表 1 不同阶段的森林林分生物量 (t·hm− 2)
Table 1 Stand biomass of different restoration stages
(t·ha− 1) 恢复时间/年
Restoration
time/year干扰类型
Interference
type树干生物量
Stem
biomass树皮生物量
Bark
biomass树枝生物量
Branch
biomass树叶生物量
Leaf
biomass地上生物量
Aboveground
biomass地下生物量
Root
biomass总生物量
Total
biomass12 轻度火烧
Mild fire78.54 11.36 23.44 5.46 118.80 41.63 160.43 中度火烧
Moderate fire23.70 3.45 7.72 1.90 36.77 13.56 50.33 重度火烧
Severe fire1.39 0.38 0.74 0.26 2.77 1.09 3.86 22 轻度火烧
Mild fire68.92 9.39 18.87 4.19 101.37 34.93 136.30 中度火烧
Moderate fire57.45 7.60 15.06 3.31 83.42 28.81 112.23 重度火烧
Severe fire54.82 7.34 14.65 3.28 80.09 27.80 107.89 31 轻度火烧
Mild fire42.72 8.70 22.81 4.26 78.49 24.09 102.58 中度火烧
Moderate fire30.71 6.48 13.56 3.28 54.03 17.64 71.67 重度火烧
Severe fire28.69 6.00 10.77 3.08 48.54 16.17 64.71 14 采伐 Cutover 1.76 0.51 0.97 0.34 3.58 1.44 5.02 20 采伐 Cutover 52.56 10.23 15.42 4.69 82.90 21.27 104.17 29 采伐 Cutover 26.45 5.65 9.79 2.90 44.79 14.38 59.17 表 2 不同阶段森林林分生物量的径阶分配
Table 2 Forest biomass proportion in different diameter classes of varied restoration stages
% 径阶
Diameter
class/cm12年
12 years22年
22 years31年
31 years14年
14 years20年
20 years29年
29 yearsLF MF SF LF MF SF LF MF SF C C C 2 — 0.48 6.1 — — — 0.09 0.52 0.75 9.04 0.36 1.58 4 0.14 0.48 21.77 0.08 — — 0.36 1.22 5.49 43.98 4.18 9.21 6 0.52 — 47.28 0.22 0.06 0.29 0.29 1.67 16.63 36.8 11.22 21.53 8 1.91 1.33 24.85 0.12 — 0.46 0.93 4.73 26.18 10.19 19.13 32.99 10 5.25 3.07 — 0.58 0.7 1.16 1.78 13.81 17.46 — 28.35 9.11 12 6.23 5 — 1.03 — 2.12 1.91 14.67 10.08 — 22.94 2.59 14 6.08 19.66 — 2.91 2.08 5.97 5.03 21.89 3.78 — 8.35 2.35 16 11.47 11.91 — 4.77 7.31 13.04 8.97 20.45 3.01 — 4.41 — 18 8.15 29.65 — 18.25 9.31 7.25 3.68 9.57 3.44 — 1.07 5.82 20 17.93 9.11 — 17.97 15.23 18.03 14.59 — 5.33 — — 5.09 22 10.38 19.31 — 16.07 10.86 8.77 11.79 8.75 3.34 — — 3.31 24 8.81 — — 22.55 16.43 12.86 12.62 2.72 — — — 6.42 26 8.76 — — 3.95 21.37 11.51 8.6 — 4.5 — — — 28 4.36 — — 9.28 11.94 — 7.12 — — — — — 30 5.27 — — 2.21 — 7.69 — — — — — — 32 — — — — — 3 4.53 — — — — — 34 — — — — — — — — — — — — 36 — — — — — 7.85 — — — — — — 38 4.69 — — — — — — — — — — — 40 — — — — — — 7.22 — — — — — 50 — — — — — — 10.5 — — — — — 表 3 林分生物量与各因子的方差分析
Table 3 Variance analysis result of forest biomass with other factors
影响因子 Influencing factor 平方和 Sum Sq. 自由度 df 均方 Mean Sq. F 值 F value 显著性 Sig. 恢复时间 Restoration time (T) 520.9 1 520.9 8.6 0.007 海拔 Altitude (AL) 887.7 1 887.7 14.7 0.001 干扰类型 Interference type (I) 2 920.0 3 973.3 16.1 < 0.001 优势树种 Dominant species (S) 454.7 2 227.3 3.8 0.036 变量 Variable (Nˉd2) 7 545.2 1 7 545.2 124.8 < 0.001 残差 Residual 1 633.0 27 60.5 模型 Model 48 002.0 8 6 000.2 99.2 < 0.001 -
[1] Tang S Z, Liu S R. Conservation and sustainability of natural forests in China[J]. Review of China Agricultural Science and Technology, 2000, 2(1): 42−46.
[2] 夏自谦, 滕秀玲. 世界森林资源现状及前景展望[J]. 北京林业大学学报(社会科学版), 2003, 2(3):24−28. doi: 10.3969/j.issn.1671-6116.2003.03.006 Xia Z Q, Teng X L. Status quo and prospect of world forest resources[J]. Journal of Beijing Forestry University (Social Sciences), 2003, 2(3): 24−28. doi: 10.3969/j.issn.1671-6116.2003.03.006
[3] Dobson A P, Bradshaw A D, Baker A J M. Hopes for the future: restoration ecology and conservation biology[J]. Science, 1997, 277: 515−522. doi: 10.1126/science.277.5325.515
[4] 李俊清. 森林生态学[M]. 北京: 高等教育出版, 2010. Li J Q. Forest ecology[M]. Beijing: Higher Education Press, 2010.
[5] 徐化成. 中国大兴安岭森林[M]. 北京: 科学出版社, 1998. Xu H C. Daxing’anling forest, China[M]. Beijing: Science Press, 1998.
[6] 郑焕能, 温广玉. 林火灾变阈值[J]. 火灾科学, 1999, 8(3):1−5. Zheng H N, Wen G Y. Forest fire catastrophe and its threshold[J]. Fire Safety Science, 1999, 8(3): 1−5.
[7] Houghton R A. Converting terrestrial ecosystems from sources to sinks of carbon[J]. Arnbio, 1996, 25(4): 267−272.
[8] Hornbeck J W, Smith C T, Martin C W, et al. Effects of intensive harvesting on nutrient capitals of three forest types in New England[J]. Forest Ecology and Management, 1990, 30(1−4): 55−64. doi: 10.1016/0378-1127(90)90126-V
[9] 张厚华, 傅德志, 孙谷畴. 森林植被恢复重建的理论基础[J]. 北京林业大学学报, 2004, 26(1):97−99. doi: 10.3321/j.issn:1000-1522.2004.01.020 Zhang H H, Fu D Z, Sun G C. Theoretical fundamentals of restoration and reconstruction of degraded forest vegetation[J]. Journal of Beijing Forestry University, 2004, 26(1): 97−99. doi: 10.3321/j.issn:1000-1522.2004.01.020
[10] Amiro B. Net primary productivity following forest fire for Canadian ecoregions[J]. Canadian Journal of Forest Research, 2000, 30(30): 939−947.
[11] Wang C, Gower S T, Wang Y, et al. The influence of fire on carbon distribution and net primary production of boreal Larix gmelinii forests in northeastern China[J]. Global Change Biology, 2010, 7(6): 719−730.
[12] 丛燕, 魏荣华, 袁强. 黑河地区阔叶林和白桦落叶松林火烧迹地生物量动态变化的研究[J]. 林业科技, 2011, 36(2):17−19. doi: 10.3969/j.issn.1673-3290.2011.02.03 Cong Y, Wei R H, Yuan Q. Research of the biomass dynamics in burned areas of broad-leaved forest and birch and larch forest in Heihe Area[J]. Forestry Science & Technology, 2011, 36(2): 17−19. doi: 10.3969/j.issn.1673-3290.2011.02.03
[13] 解伏菊, 肖笃宁, 李秀珍, 等. 大兴安岭北坡火烧迹地森林景观恢复及其影响因子: 以郁闭度指标为例[J]. 应用生态学报, 2005, 16(9):1711−1718. doi: 10.3321/j.issn:1001-9332.2005.09.025 Xie F J, Xiao D N, Li X Z, et al. Forest landscape restoration and its affecting factors in burned area of northern Great Xing’an Mountains: taking forest coverage as an example[J]. Chinese Journal of Applied Ecology, 2005, 16(9): 1711−1718. doi: 10.3321/j.issn:1001-9332.2005.09.025
[14] 王晓莉, 常禹, 陈宏伟, 等. 黑龙江省大兴安岭主要森林生态系统生物量分配特征[J]. 生态学杂志, 2014, 33(6):1437−1444. Wang X L, Chang Y, Chen H, et al. Biomass allocation characteristics of the main forest ecosystems in the Great Xing’an Mountains, Heilongjiang Province[J]. Chinese Journal of Ecology, 2014, 33(6): 1437−1444.
[15] 孙伟迪. 乌尔旗汗林业局森林碳储量分布特征研究[J]. 内蒙古林业调查设计, 2017, 40(6):96−100, 102. Sun W D. Study on distribution characteristics of forest carbon storage in Wuerqihan Forestry Bureau[J]. Inner Mongolia Forestry Survey and Design, 2017, 40(6): 96−100, 102.
[16] 孟宪宇. 测树学[M]. 北京: 中国林业出版社, 2006. Meng X Y. Forest mensuration[M]. Beijing: China Forestry Publishing House, 2006.
[17] Wang T, Wang G, Innes J L, et al. ClimateAP: an application for dynamic local downscaling of historical and future climate data in Asia Pacific[J]. Frontiers of Agricultural Science and Engineering, 2017, 4(4): 448−458. doi: 10.15302/J-FASE-2017172
[18] Weiskittel A R, Crookston N L, Radtke P J. Linking climate, gross primary productivity, and site index across forests of the western United States[J]. Canadian Journal of Forest Research, 2011, 41(8): 1710−1721. doi: 10.1139/x11-086
[19] 国家林业局. 国家森林资源连续清查技术规定[S]. 北京: 国家林业局, 2014. State Forestry Administration. National forest resources continuous inventory technical regulations[S]. Beijing: State Forestry Administration, 2014.
[20] 国家林业局. LY/T 2655—2016立木生物量模型及碳计量参数: 云杉[S]. 北京: 中国标准出版社, 2017. State Forestry Administration. LY/T 2655—2016, Tree biomass models and related parameters to carbon accounting for Picea asperata[S]. Beijing: China Forestry Publishing House, 2017.
[21] 国家林业局. LY/T 2654—2016立木生物量模型及碳计量参数: 落叶松[S]. 北京: 中国标准出版社, 2017. State Forestry Administration. LY/T 2654—2016, Tree biomass models and related parameters to carbon accounting for Larix gmelinii[S]. Beijing: China Forestry Publishing House, 2017.
[22] 国家林业局. LY/T 2659—2016, 立木生物量模型及碳计量参数—桦树[S]. 北京: 中国标准出版社, 2017. State Forestry Administration. LY/T 2659—2016, Tree biomass models and related parameters to carbon accounting for Betula platyphylla[S]. Beijing: China Forestry Publishing House, 2017.
[23] 李海奎, 雷渊才. 中国森林植被生物量和碳储量评估[M]. 北京: 中国林业出版社, 2010. Li H K, Lei Y C. Estimation and evaluation of forest biomass carbon storage in China[M]. Beijing: China Forestry Publishing House, 2010.
[24] 唐守正, 郎奎建, 李海奎. 统计和生物数学模型计算: ForStat教程[M]. 北京: 科学出版社, 2009. Tang S Z, Lang K J, Li H K.Statistics and biomathematical models[M]. Beijing: Science Press, 2009.
[25] Kutner M H, Nachtsheim C J, Neter J, et al. Applied linear statistical models[M]. Boston: McGraw-Hill Irwin, 2005.
[26] Farrar D E, Glauber R R. Multicollinearity in regression analysis: the problem revisited[J]. Review of Economics and Statistics, 1967, 49(1): 92−107. doi: 10.2307/1937887
[27] Mcintyre S H, Montgomery D B, Srinivasan V, et al. Evaluating the statistical significance of models developed by stepwise regression[J]. Journal of Marketing Research, 1983, 20(1): 1−11. doi: 10.1177/002224378302000101
[28] Cornell J. Classical and modern regression with applications[J]. Technometrics, 1987, 29(3): 2.
[29] 曾锋, 张金池. 重金属在森林生态系统中的迁移规律研究进展[J]. 世界林业研究, 2001, 14(2):16−22. Zeng F, Zhang J C. Advances in the research on movement pattern of heavy metal ions in forest ecosystem[J]. World Forestry Research, 2001, 14(2): 16−22.
[30] 范兆飞, 徐化成, 于汝元. 大兴安岭北部兴安落叶松种群年龄结构及其与自然干扰关系的研究[J]. 林业科学, 1992, 28(1):2−11. Fan Z F, Xu H C, Yu R Y. A study on the species group age structure of Larix Gmelini population and its relation to disturbance on the North Daxing’anling Mountains[J]. Scientia Silvae Sinicae, 1992, 28(1): 2−11.
[31] 李卓凡. 兴安落叶松林生物量与碳储量的研究[D]. 呼和浩特: 内蒙古农业大学, 2013. Li Z F. Biomass and carbon storage of the Larix gmelinii forest’s research[D]. Hohhot: Inner Mongolia Agricultural University, 2013.
[32] 王丽红, 辛颖, 邹梦玲, 等. 大兴安岭火烧迹地植被恢复中植物多样性与生物量分配格局[J]. 北京林业大学学报, 2015, 37(12):41−47. Wang L H, Xin Y, Zou M L, et al. Plants diversity and biomass distribution of vegetation restoration in burned area of Great Xing’an Mountains[J]. Journal of Beijing Forestry University, 2015, 37(12): 41−47.
[33] Brown I F, Martinelli L A, Thomas W W, et al. Uncertainty in the biomass of Amazonian forests: an example from Rondonia, Brazil[J]. Forest Ecology and Management, 1995, 75(1−3): 175−189.
[34] Condés S, Roberedo F G. An empirical mixed model to quantify climate influence on the growth of Pinus halepeensis Mill. stands in southeastern Spain[J]. Forest Ecology and Management, 2012, 284(15): 59−68.
[35] 王轶夫. 基于神经网络的森林生物量估测模型研究[D]. 北京: 北京林业大学, 2013. Wang Y F. The study of forest biomass estimation model based on neural network[D]. Beijing: Beijing Forestry University, 2013.
[36] Taylor K D. Fire history of a sequoia-mixed conifer forest[J]. Ecology, 1979, 60(1): 129−142. doi: 10.2307/1936475
[37] Reed R A, Finley M E, Romme W H, et al. Above-ground net primary production and leaf area index in early postfire vegetation in Yellowstone National Park[J]. Ecosystems, 1999, 2(1): 88−94. doi: 10.1007/s100219900061
[38] 张玉红, 覃炳醒, 孙铭隆, 等. 林火对大兴安岭典型林型林下植被与土壤的影响[J]. 北京林业大学学报, 2012, 34(2):7−13. Zhang Y H, Qin B X, Sun M L, et al. Impact of forest fire on understory vegetation and soil in typical forest types of Daxing’an Mountains, northeastern China[J]. Journal of Beijing Forestry University, 2012, 34(2): 7−13.
[39] 郑丽凤, 周新年, 李丹, 等. 森林采伐对闽北天然次生林碳储量及其动态的影响[J]. 安全与环境学报, 2013, 13(6):162−167. Zheng L F, Zhou X N, Li D, et al. Impact of harvesting intensity of the natural secondary forestry on the carbon stocks reduction in northern Fujian[J]. Journal of Safety and Environment, 2013, 13(6): 162−167.
[40] 胡海清, 罗碧珍, 魏书精, 等. 小兴安岭7种典型林型林分生物量碳密度与固碳能力[J]. 植物生态学报, 2015, 39(2):140−158. doi: 10.17521/cjpe.2015.0014 Hu H Q, Luo B Z, Wei S J, et al. Biomass carbon density and carbon sequestration capacity in seven typical forest types of the Xiaoxing ’an Mountains, China[J]. Chinese Journal of Plant Ecology, 2015, 39(2): 140−158. doi: 10.17521/cjpe.2015.0014
[41] He Y J, Qin L, LI Z Y, et al. Carbon storage capacity of monoculture and mixed-species plantation in subtropical China[J]. Forest Ecology and Management, 2013, 295(1): 193−198.
[42] Fan H B, Liu W F, Wu J P, et al. Ecosystem carbon pools in mixed stands of hardwood species and Masson pine[J]. Journal of Tropical Forest Science, 2013, 25(2): 154−165.
[43] 李江, 陈宏伟, 冯弦. 云南热区几种阔叶人工林C储量的研究[J]. 广西植物, 2003, 23(4):294−298. Li J, Chen H W, Feng X. Carbon stock and rate of carbon sequestration assessment of hardwood plantations in tropical Yunnan, China[J]. Guihaia, 2003, 23(4): 294−298.