基于地形梯度的河北省怀来县山地生态系统空间分布及变化特征分析

许 策, 李 超, 张天柱, 束继年

(1.中国科学院 地理科学与资源研究所, 北京 100101； 2.中国科学院 生态环境研究中心,北京 100085; 3.中国科学院大学; 4.中国农业大学 资源与环境学院, 北京 100193)

基于地形梯度的河北省怀来县山地生态系统空间分布及变化特征分析

许 策1,2,3, 李 超4, 张天柱4, 束继年2,3

(1.中国科学院 地理科学与资源研究所, 北京 100101； 2.中国科学院 生态环境研究中心,北京 100085; 3.中国科学院大学; 4.中国农业大学 资源与环境学院, 北京 100193)

[目的] 从地形梯度角度,对区域生态系统空间分布及变化特征进行分析,提出基于地形梯度的区域生态系统保护对策建议,为制定山区生态规划及实现山区生态系统的持续、健康发展提供参考。[方法] 以遥感解译获得的1993,2003,2013年生态系统类型空间分布数据为基础,选取高程级别、地形起伏度和坡度变率为地形梯度表征指标,采用分布指数,对区域农田、水域与湿地、森林、草地、聚落和荒漠6类生态系统的地形梯度差异分布情况及优势位进行分析,并结合1993—2013年各生态系统的服务价值变化特征,探讨区域生态系统的保护对策。[结果] 农田、水域与湿地、聚落、草地生态系统集中于低地形梯度区;森林集中于高地形梯度区;荒漠生态系统分布特征随地形梯度差异不明显;在低地形梯度区,农田、水域与湿地的生态系统服务价值降低,草地、聚落生态系统服务价值有明显的增加趋势;在高地形梯度区,虽然森林生态系统服务价值变化较为剧烈,但总体呈现出持续增加的趋势。[结论] 在制定针对山地的区域生态规划时,应考虑对位于优势位上的生态系统实行特殊保护,禁止被侵占为生态系统服务功能较弱的生态系统服务类型,以保证生态系统能够持续、健康地提供生态服务。

地形梯度; 生态系统; 空间分布; 优势位

生态系统服务作为当今社会关注的焦点以及生态学、地理学等学科的研究前沿和热点，引起了众多学者与研究组织的广泛关注。近些年来，随着社会经济的快速发展以及土地利用方式的转变，区域生态系统格局也发生了相应的改变，使生态系统的稳定性、可持续发展受到了严峻挑战。与平原地区相比，山区的生态环境更加脆弱，其社会经济与生态环境保护的和谐发展显得尤其重要。山区生态系统格局的改变会使社会的生态安全风险问题更加凸显，生态系统服务的改变甚至丧失将会使人类面临着的生存与发展的多重挑战[1]。

地形是区域生态系统中的重要组成因素，其变化将会影响生态系统内部的物质流动和能量转化。区域地形特征的差异也会引起土地利用方式的空间差异，进而影响生态系统的空间分布格局[2]。地形梯度逐渐被广泛用于土地利用格局变化方面的相关研究[3]。

生态系统的内部结构及分布格局反映了不同生态系统的自身特征以及各类生态系统之间的空间结构关系，决定了生态系统服务功能的整体状况及其空间差异,也为人类因地制宜对不同区域的生态系统服务功能采取不同的保护和利用措施提供了重要依据。近年来，国内外众多学者对生态系统的空间分布格局及其结构变化开展了一系列相关研究[4-5]。但从地形梯度角度，对区域生态系统格局及空间结构变化情况展开分析的研究较为罕见。

本文拟从高程、地形起伏度、坡度变率3个方面，揭示研究区生态系统在地形梯度上的分布特征及不同地形梯度上的区域生态系统服务价值变化，并提出基于地形梯度的区域生态系统保护对策建议，以期为山区生态规划的制定及实现生态系统的持续、健康发展提供参考。

1 研究区概况

怀来县坐标地处115°16′—115°58′E，40°04′—40°35′N，位于张家口市东南部，紧邻北京市昌平区和延庆区，距离北京市中心约100 km，土地面积约1 973.18 km2。由于地处燕山山脉西北侧，永定河上游，属于中温带半干旱区，光照充足，降水稀少，生态环境较为敏感、脆弱。怀来县境内主要地貌类型为河川平原、丘陵和山地等，其南部和北部多为山地，占全县面积的绝大部分，中部夹有山间盆地，海拔大致在320～1 954 m。怀来县作为北京及其周边地区的绿色屏障和生态缓冲区，其生态环境状况对北京地区的生态安全问题有不可忽视的影响。

2 数据来源与研究方法

2.1 数据获取与处理

地理空间数据云网站下载获得1993，200与2013年空间分辨率为30 m的怀来县Landsat TM/ETM+影像以及空间分辨率为30 m的DEM图，考虑研究区的气候特点，并保证影像的清晰，因此，影像选择在8—9月。以生态系统分类方法、土地遥感调查分类标准为基础，结合研究区生态系统分布特征及《中国5 a间隔陆地生态系统空间分布数据集》，将怀来县生态系统类型划分为森林生态系统、农田生态系统、草地生态系统、聚落生态系统、水体与湿地生态系统、荒漠生态系统6个生态系统(表1)。将Landsat影像经大气校正、辐射定标、投影变换、掩膜提取等预处理后，采用监督分类并结合人机交互的方法完成遥感影像的生态系统解译分类，并利用实地随机选取若干样点进行实际测量的数据来检验解译精度，最后利用混淆矩阵法进行精读评价，结果表明1993，2003与2013年3期遥感影像的解译分类精度分别为83.83%，83.21%和87.69%。运用ArcGIS软件对DEM影像进行空间分析得出地形起伏度、高程、坡度变率、地形位指数分布图。

2.2 研究方法

2.2.1 地形梯度等级划分 综合考虑地形因子的可定量性和对生态系统空间分布及变化的影响程度，选取高程(altitude)、地形起伏度(landform relief)、坡度变率(slope change rate)分析地形对生态系统类型空间分布及变化的影响[6]。其中，高程级别数据可运用ArcGIS对DEM图进行重分类获得[7]；地形起伏度数据可运用ArcGIS栅格邻域分析工具计算获得[8]；地面坡度变率数据可依据坡度计算原理，利用ArcGIS软件在坡度分布图的基础上再次求取坡度。

表1 怀来县生态系统类型划分与《全国遥感监测土地利用覆盖分类体系》对照

结合研究区地形特点，本文对高程、地形起伏度、坡度变率因子进行分级(表2)。采用等间隔重分类法[9]对地形位指数(0.69～3.03)分为10级，并将这10个地形位级别分为低、中、高3个梯度。运用ArcGIS空间分析功能，对各地形因子进行重分类，统一输出分辨率30 m的栅格数据(本文将栅格计算中地形因子划分的等级视为地形分区)。

表2 怀来县高程、地形起伏度、坡度变率分级及面积比例

2.2.2 地形梯度表征指标 为消除面积差异产生的量纲影响，采用分布指数(distribution index， DI)反映研究区各类生态系统在不同地形梯度的分布特征，公式如下[10-11]：

(1)

式中：P——分布指数；e——地形因子，分别为地形起伏度、坡度变率、坡向和地形位；Sie——e地形因子特定等级下的i类生态系统的面积(hm2)；Si——i类生态系统的面积(hm2)；Se——整个研究区e地形因子特定等级下的全部生态系统总面积(hm2)；S——研究区总面积(hm2)。分布指数变化越小，表明某类生态系统分布与标准分布的偏离越小，其对地形差异的适宜性越大；反之，则表明某类生态系统的分布对地形条件的选择较为苛刻[12]。当p=1时，表示某类生态系统在某种地形上的比重等于该类生态系统总面积在研究区的比重；当p>1时，表明某类生态系统在该地形上的比重大于该类生态系统总面积在研究区的比重，故将p>1的区间设定为该生态系统的优势位。

2.2.3 生态系统服务价值测算 参考谢高地等[13]研究成果，单位当量因子的经济量为全国平均粮食单产市场价格的1/7，中国2005年单个生态系统价值当量的经济价值为707.92元/hm2，结合2006年《河北省统计年鉴》的相关数据，计算得出2005年河北省平均粮食产量为4 683.35 kg/hm2。2005年全国平均粮食产量为5 896.50 kg/hm2，据此对全国生态系统服务价值当量价值进行系数修正，并通过2005与2013年间CPI指数修正物价水平间的差异，确定怀来县2013年单个生态当量的价值为732.95元/hm2。由于陆地生态系统单位面积生态系统服务价值系数=单位当量价值×当量，因此借助中国陆地生态系统服务价值当量价值表[14]，计算得出单位面积各生态系统服务功能价值系数(VC)，结合生态系统空间分布数据，可得出研究区生态系统服务价值量。具体公式为：

(2)

式中：ESV——生态系统年服务总价值(元/a)；Ai——研究区第i种生态类型的分布面积(i=1，2，…，5，为生态系统类型)(hm2)； VCi——第i种生态系统服务价值当量〔元/(hm2·a)〕。

根据研究区的特点，以上述Costanza的计算模型为基础，建立区域生物量差异性系数S，适当修正干扰生态系统服务价值的某些因素[15-16]，更改后的计算模型如下：

(3)

式中：fi——研究区内i类生态系统的区域生物量差异性系数。下同。

生态系统服务功能的大小与其生物量密切相关，生物量越大，生态服务功能越强。分别从《中国1∶100万植被类型图》和全国NPP数据中提取研究区生态系统类型和生物量数据，叠加得到研究区和全国生态系统平均生物量数据，利用区域生物量差异系数对各生态类型的生态服务价值进行空间修正，假定生态系统服务功能大小与生物量大小成线性关系，用区域生物量差异性系数来定量表达生物量大小对生态服务功能的影响，公式如下：

(4)

式中：Ci——研究区i类生态系统的平均生物量；C——中国一级 类生态系统类型平均生物量。

3 结果与分析

3.1 基于地形梯度的生态系统空间分布特征分析

运用ArcGIS的raster calculator功能，将通过遥感解译的各年份生态系统类型空间分布图与DEM高程分级图进行叠加，得出各年份生态系统类型的高程梯度差异分布情况(表3)。从生态系统类型分布上看，随着高程级别的上升，水域与湿地生态系统的面积和分布指数呈下降趋势，且变化幅度最大，在300～500 m高程级别上的分布指数均>5，是绝对优势的生态系统类型，而高程上升到500 m以上时，其分布指数降低至<0.4，说明水域与湿地生态系统类型对高程的适应性较差。

随高程级别的上升，农田、聚落、草地及荒漠生态系统的面积呈现先上升后下降的变化趋势，其中，农田和草地生态系统的优势位分布在300～500和500～700 m高程级别上；聚落生态系统除1993年的优势位在300～500 m高程级别上外，其余2个年份的优势位均在500～700 m高程级别上，说明在500～700 m高程级别上的人类活动最为活跃；荒漠生态系统的优势位分布在500～1 100 m高程级别上，而且在700～900 m高程级别上达到最优，以未利用地为主；森林生态系统作为怀来县分布面积最大的生态系统，主要分布于高程>900 m的区域，而且在900～1 100和1 100 m高程以上的优势位显著，该区域应作为生态环境重点保护区域。

表3 怀来县1993-2013年各生态系统的高程梯度差异

地形起伏度是表征地貌类型的重要指标，运用ArcGIS的Raster Calculator功能，将通过遥感解译的各年份生态系统类型空间分布图与地形起伏度图进行叠加，得出各年份生态系统类型的地形起伏度梯度差异分布情况(表4)。从各个生态系统的分布来看，不同生态系统类型在不同地形起伏度上的分布表现出明显的梯度变化特征。地形起伏度在30 m以下和30～70 m的区域是水域和湿地生态系统的优势位，地形起伏度超过200 m的区域很少再有水域和湿地的分布；森林生态系统的面积以及分布指数与地形起伏度级别呈现明显的正相关趋势，其优势位分布于70～200和200 m以上的起伏度级别上；农田面积及其分布指数与地形起伏度的级别的呈现明显的负相关趋势，其优势位在30 m以下的起伏度级别上，在优势位范围内3个不同时点的分布指数分别为2.27，2.27，2.30，均高于同期其他5类生态系统的分布指数，这表明农田生态系统在0～30 m地形起伏度级别上具有相对较强的稳定性和适应性；聚落生态系统和和草地生态系统的面积和分布指数均与地形起伏度级别呈负相关，其中聚落生态系统的优势位分布于0～30 m起伏度级别上，表明人类的生产建设活动大多都在此起伏度范围内；草地生态系统的优势位分布于0～30和30～70 m起伏度级别上；比较特殊的是荒漠生态系统，因为荒漠生态系统除了包括沙地、戈壁、盐碱地等还包括裸地和裸岩，所以其在30～70，70～200及200 m以上的起伏度级别区域均有分布。

表4 怀来县1990-2010年各生态系统的地形起伏度梯度差异

运用ArcGIS的Raster Calculator功能，将通过遥感解译的各年份生态系统类型空间分布图与坡度变率图进行叠加，得出各年份生态系统类型的坡度变率梯度差异分布情况(表5)。

从生态系统类型分布上看，随坡度变率的增加，各生态系统呈现出明显的坡度变率梯度分布特征，农田、聚落和草地的优势位均分布在0°～2°和2°～6°两个坡度变率级别上，森林和荒漠的优势位则分布在大于6°的坡度变率级别上，水域与湿地的优势位主要分布在0°～2°和2°～6°坡度变率级别上，但在2003年之后，6°～15°上也成为其分布的优势位。

3.2 基于地形梯度的生态系统服务价值时空变化分析

以遥感解译获得的生态系统空间分布图为基础，运用生态系统服务价值计算公式，分别计算1993，2003，2013年怀来县生态系统服务价值，计算得到1993—2013年怀来县生态系统服务价值变化，分别与DEM高程分级、地形起伏度、坡度变率叠加，汇总得出基于地形梯度的生态系统服务价值变化情况(图1)。

表5 怀来县1993-2013年各生态系统的坡度变率梯度差异

图1a为1993—2013年不同高程级别上的生态系统服务价值变化情况。农田、水域与湿地生态系统服务价值在各高程级别上均表现为持续减少，其中农田、水域与湿地分别在500～700和300～500 m高程级别上的减少量最为明显，20 a间分别减少了2 941.42万元和1.42亿元；森林、草地、聚落生态系统服务价值在各高程级别上均表现为持续增加，其中森林、草地、聚落分别在1 100 m以上、300～500和500～700 m高程级别上的增加量最为明显，20 a间分别增加了3 955.07,1 920.62和496.74万元；荒漠生态系统服务价值与其他5个生态系统相比，在各高程级别上的变化规律不太明显。

图1b为1993—2013年不同坡度变率级别上的生态系统服务价值变化情况。农田、草地、水域与湿地、聚落生态系统服务价值量在0°～6°坡度变率较小的区域变化较为剧烈，其中，农田、水域与湿地在0°～6°坡度变率级别上的价值量持续减少，20 a间分别减少了3 613.02,1.38亿元；草地、聚落在0°～6°坡度变率级别上的价值量持续增加，20 a间分别增加了3 575.84,556.17万元。森林生态系统服务价值量在6°以上坡度变率较大的区域变化较为剧烈，20 a间共增加了6 184.50万元。从怀来县整体来看，在0°～6°坡度变率级别区域的生态系统服务价值量呈减少趋势，20 a间共减少了1.26亿元；而在6°以上坡度变率级别区域的生态系统服务价值量呈增加趋势，20 a间共增加了5 950.45万元。草地在300～500 m高程级别上的增加量，水域与湿地生态系统的价值量持续减少，20 a间减少了1.42亿元；草地生态系统的价值量增加明显，20 a间增加了1 920.62万元。在500～700和700～900 m高程级别上，生态系统服务价值变化情况较为相似，其中农田生态系统的价值量表现为持续减少，20 a间共减少了3 805.60万元；森林、草地生态系统的价值量持续增加，20 a间共增加了2 692.85万元。

图1c为1993—2013年不同地形起伏度级别上的生态系统服务价值变化情况。在0～30 m低起伏度级别区域，各生态系统的价值量均有所增加，20 a间共增加467.82万元；在30 m以上起伏度级别区域，各生态系统的价值量变化产生分异，农田、水域与湿地、荒漠的价值量持续减少，20 a间分别减少了4 406.81万元、1.50亿元、279.64万元；森林、草地、聚落的价值量持续增加，20 a间分别增加了6 622.22,7 456.80和643.11万元。从各生态系统的价值量变化来看，森林生态系统在200 m以上起伏度级别上的价值增加量最为明显、农田、水域与湿地在30～70 m起伏度级别上的价值减少量最为明显。

图1 基于地形梯度的1993-2013年怀来县生态系统服务价值空间变化

3.3基于地形梯度的区域生态系统保护对策建议

1993—2013年，怀来县生态系统服务价值总量呈减少趋势，20 a间共减少6 631.07万元。生态系统服务价值持续降低，说明当前生态系统的空间分布格局不利于当地生态文明建设。高地形梯度区，由于人为干扰较少，加之近年来荒山治理、封山育林、退耕还林等政策的实施，使得生态系统服务价值量有所增加，说明高地形梯度区的生态系统运行良好。因此，造成怀来县生态系统恶化的主要原因，在于低地形梯度区域生态系统的不合理利用。从生态系统服务价值变化特征来看，低地形梯度区由于人为活动较为剧烈，建设用地的不断扩张，导致农田、水域与湿地、森林、草地等生态系统服务价值量较高的区域不断被侵占，是造成低地形梯度区生态系统服务价值减少的根本原因。例如在0°～6°低坡度变率级别上，1993—2013年农田、水域与湿地面积共减少了约10 828.12 hm2，但与此同时，聚落面积增加了5 316.05 hm2，这直接导致生态系统服务价值降低了约1.68亿元。

从前文研究可以看出，不同生态系统均有其对应的地形梯度优势位(表6)，笔者认为，实现生态系统持续健康发展的必要条件，是要根据各生态系统的地形梯度优势位，对区域生态系统进行科学、合理空间布局，制定生态规划，以确保生态系统的稳定性。根据生态规划的定义：“生态规划是以生态学原理为指导，以区域生态系统整体优化为基本目标，应用环境科学、系统科学、生态学等多学科的知识来组织、协调生态系统内部各种生态关系，是一种确定资源开发利用和保护的生态适宜性，促进人与生态环境协调、持续发展的规划方法。”因此，生态规划必须建立在自然综合体不做重大变化的前提下进行，故不能将区域各生态系统进行重新布局，应以保护现有适宜于生态系统健康发展的生态类型为规划思路。建议山地在制定区域生态规划时，注意对位于地形梯度优势位上的各生态类型实行特殊保护，禁止被侵占为生态系统服务功能较弱的生态类型，只有这样才能保证生态系统持续、健康地提供服务。

表6 怀来县各生态系统的地形梯度优势位特征

4 结 论

(1) 1993—2013年怀来县生态系统空间分布随地形梯度差异呈现明显的层级性。农田、水域与湿地、聚落、草地生态系统集中于低地形梯度区；森林集中于高地形梯度区；荒漠生态系统分布特征随地形梯度差异不明显。

(2) 1993—2013年怀来县生态系统服务价值变化随地形梯度差异呈现明显的空间分异特征。在低地形梯度区，农田、水域与湿地的生态系统服务价值降低，草地、聚落生态系统服务价值有明显的增加趋势；在高地形梯度区，虽然森林生态系统服务价值变化较为剧烈，但总体呈现出持续增加的趋势。

(3) 在山区的生态规划过程中，对位于优势位的生态系统要考虑实施实施特殊保护，防止其被侵占为生态服务功能较弱的生态类型，以保证生态系统稳定、持续、健康地提供服务。

[1] 王根绪,邓伟,杨燕,等.山地生态学的研究进展、重点领域与趋势[J].山地学报,2011,29(2):129-140.

[2] 龚文峰,袁力,范文义.基于地形梯度的哈尔滨市土地利用格局变化分析[J].农业工程学报,2013,29(2):250-259.

[3] 郭洪峰,许月卿,吴艳芳.基于地形梯度的土地利用格局与时空变化分析:以北京市平谷区为例[J].经济地理,2013,33(1):160-166.

[4] 徐新良,刘纪远,邵全琴,等.30 a来青海三江源生态系统格局和空间结构动态变化[J].地理研究,2008,27(4):829-838.

[5] 朱瑜馨,赵军,曹静.祁连山山地生态系统稳定性评估模型[J].干旱区研究,2002,38(4):33-37.

[6] 温晓金,杨海娟,刘焱序.基于地形因子的km尺度景观生态廊道布局研究[J].地理科学进展,2013,32(2):298-307.

[7] 田瑞云,王玉宽,傅斌,等.基于DEM的地形单元多样性指数及其算法[J].地理科学进展,2013,32(1):121-129.

[8] 孙丕苓,许月卿,王数.环京津贫困带土地利用变化的地形梯度效应分析[J].农业工程学报,2014,30(14):277-288.

[9] 李升发,董玉祥.广州市土地利用垂直分层模式及其演变过程[J].自然资源学报,2013,28(6):1019-1029.

[10] 龚文峰,袁力,范文义.基于地形梯度的哈尔滨市土地利用格局变化分析[J].农业工程学报,2013,29(2):250-259.

[11] Mottet A, Ladet S, Coqué N, et al. Agricultural land-use change and its drivers in mountain landscapes: A case study in the Pyrenees[J]. Agriculture, Ecosystems & Environment, 2006,114(2):296-310.

[12] 郭洪峰,许月卿,吴艳芳.基于地形梯度的土地利用格局与时空变化分析:以北京市平谷区为例[J].经济地理,2013,33(1):160-166.

[13] 谢高地,鲁春霞,冷允法,等.青藏高原生态资产的价值评估[J].自然资源学报,2003,18(2):189-196.

[14] 李超,杜哲,陈亚恒,等.环京津地区土地生态服务价值时空分异特征[J].土壤通报,2015,46(1):42-47.

[15] 谢高地,甄霖,鲁春霞,等.一个基于专家知识的生态系统服务价值化方法[J].自然资源学报,2008,23(5):911-919.

[16] 石晓丽,王卫.生态系统功能价值综合评估方法与应用:以河北省康保县为例[J]生态学报,2008,28(8):3998-4006.

DistributionandVariationAnalysisofMountainEcosystemBasedonTopographicGradientinHuailaiCountyofHebeiProvince

XU Ce1,2,3, LI Chao4, ZHANG Tianzhu4, SHU Jinian2,3

(1.InstituteofGeographicSciencesandNaturalResourcesResearch,ChineseAcademyofScience,Beijing100101,China; 2.ResearchCenterforEco-EnvironmentalSciences,ChineseAcademyofSciences,Beijing100085,China; 3.UniversityofChineseAcademyofSciences,Beijing100049,China; 4.CollegeofResourcesandEnvironmentalSciences,ChinaAgriculturalUniversity,Beijing100193,China)

[Objective] To provide relevant suggestions for regional ecosystem protection and sustainable development based on topographic gradient via analyzing the spatial distribution of regional ecological system and the variation characteristic of it. [Methods] Topographic gradient index was referred to elevation level, topographic relief and slope variation rate. Distribution index was used to designate topographic gradient distribution and dominant position of cultivated land, waters and wetlands, forestland, grassland, settlement and desert and was analyzed using the data of ecosystem-type-based space distributions in the years of 1993, 2003 and 2013 which was acquired by remote sensing interpretation. Finally, the variations of each ecosystem service item in the years of 1993—2013 were analyzed to provide relevant suggestion for ecological protection. [Results] Cultivated land, waters and wetlands, settlement and grassland mainly distributed on low terrain gradient area; forestland distributed on high terrain gradient area; the distribution characteristic of desert ecosystem was not significantly influenced by gradient. In low terrain gradient area, the values of ecosystem service of cultivated land, waters and wetlands decreased while the values of grassland and settlement showed increasing trends. In high terrain gradient area, the value of ecosystem service of forestland increased as a whole but changed dramatically. [Conclusion] To make plan for regional ecological system of mountain land, the protection for the ecosystem at dominant position should be taken into consideration to prevent it from being replaced by those ecosystem service types having weak ecosystem service function, and to keep it sustainable and healthy in providing ecosystem service.

topographicgradient;ecosystemsystem;spacedistribution;dominantposition

A

1000-288X(2017)05-0198-07

X321

文献参数： 许策, 李超, 张天柱, 等.基于地形梯度的河北省怀来县山地生态系统空间分布及变化特征分析[J].水土保持通报，2017,37(5)：198-204.

10.13961/j.cnki.stbctb.2017.05.033； Xu Ce, Li Chao, Zhang Tianzhu, et al. Distribution and variation analysis of mountain ecosystem based on topographic gradient in Huailai County of Hebei Province[J]. Bulletin of Soil and Water Conservation, 2017,37(5)：198-204.DOI:10.13961/j.cnki.stbctb.2017.05.033

2017-01-26

2017-03-31

国土资源部公益性行业科研专项“环渤海盐碱土地景观整治与植物修复技术集成”(201311060)

许策(1992—),男(汉族),河北省保定市人,硕士研究生,研究方向为生态系统服务。E-mail:417716152@qq.com。

束继年(1966—),男(汉族),安徽省六安市人,博士，研究员,博士生导师,主要从事生态环境治理研究。E-mail:Jshu@rcees.ac.cn。