榆林市生态系统固碳价值对林地系统结构的响应变化_张凡.pdf

DOI10.12401/j.nwg.2023103 榆林市生态系统固碳价值对林地系统结构的响应变化 张凡1 ，田旭荣1 ，薛彤1 ，林永航1 ，齐佳敏1 ，张鹏飞2 ，张茂省1 （1. 西安交通大学人居环境与建筑工程学院，陕西 西安 712000；2. 榆林市自然资源和规划局，陕西 榆林 719000） 摘 要将生态资源纳入经济体系统筹管理是环境可持续发展的必然要求。对于生态脆弱地区， 加强环境压力变化下的资源管理尤为重要。笔者从系统结构的角度理解生态系统固碳功能变化 特征，在榆林地区结合景观结构指数对固碳功能进行空间分析。2000～2020年榆林市生态系统 固碳服务价值持续增加，2000～2010年增长较快，2010～2020年增速减缓。在空间分布方面，单 位面积固碳价值从东南向西北递减，东南部的增长速率也更快。工业活动是形成榆林市生态系 统固碳价值西北、东南分划的主要因素，东南部区县受生态系统结构变化的影响更强。研究结 果表明，在生态脆弱区实施生态建设和修复工程时，应综合考虑生态系统的结构性以提高地区 生态系统潜在的功能性。 关键词生态系统服务；固碳价值；系统结构；生态脆弱区 中图分类号P69；F205文献标志码A文章编号1009-6248（2023）03-0051-07 Response of Ecosystem Carbon Fixation Value to Forest System Structure in Yulin City ZHANG Fan1，TIAN Xurong1，XUE Tong1，LIN Yonghang1，QI Jiamin1，ZHANG Pengfei2，ZHANG Maosheng1 （1. School of Human Settlements and Civil Engineering, Xi’an Jiaotong University, Xi’an 712000, Shaanxi, China；2. Yulin Bureau of Natural Resources and Planning, Yulin 719000, Shaanxi, China） AbstractIt is an inevitable requirement of environmentally sustainable development to bring ecological re- sources into the overall management of the economic system. For ecologically fragile areas, it is particularly im- portant to strengthen resource management under changing environmental pressures. This paper attempts to un- derstand the change characteristics of carbon fixation function of ecosystem from the perspective of system structure, and carries out spatial analysis of carbon fixation function in Yulin area combined with landscape structure index. From 2000 to 2020, the carbon sequestration service value of Yulin’s ecosystem continued to in- crease, and the growth rate was fast in 2000～2010 and slowed down in 2010～2020. In terms of spatial distri- bution, the carbon sequestration value per unit area s decreases from southeast to northwest, and the growth rate in southeast is also faster. Industrial activities are the main factor to form the northwest and southeast division of carbon sequestration value of Yulin, and the southeast counties are more affected by the change of ecosystem structure. The results of the study suggest that we should comprehensively consider the structure of the ecosys- tem to improve the potential functionality of the regional ecosystem when implementing ecological construction and restoration projects in the ecological fragile areas in the future. 收稿日期2023-02-16；修回日期2023-05-23；责任编辑吕鹏瑞 基金项目国家自然科学基金面上项目（42173019）资助。 作者简介张凡（1978−），女，副教授，博士，从事生态评价研究工作。E−mailzhangfanxjtu.edu.cn。 第 56 卷 第 3 期西 北 地 质Vol. 56 No. 3 2023 年 （总 229期 ） NORTHWESTERN GEOLOGY 2023（Sum229） Keywordsecosystem services；carbon fixation value；system structure；ecologically fragile area 实现“双碳”目标是国家根据世界发展新形势 作出的重大战略决策。如何着眼地区经济实际和资 源环境禀赋，统筹产业结构与发展布局，推动高质量 发展，是地方政府助力“双碳”目标的重要举措。生 态系统是地球系统正常运转的基础。2030年前碳达 峰行动方案对生态系统碳汇能力巩固提升提出明确 要求。对区域生态资源赋存状态与变化进行准确评 估，是构建“双碳”目标下国土空间开发与生态环境 保护新格局的基础。基于自然的资本属性，根据生态 系统及其固碳功能服务衡量经济价值（Costanza et al.， 1997）的区域生态系统价值评估，不仅可以为从国家 层面探索自然资源资产负债表提供技术支持，也是地 区可持续发展目标导向政策优化调整的基础。 国内很多学者已对中国生态系统固碳服务价值 开展评估研究。早期研究表明，中国陆地生态系统效 益价值约为GDP的1.25倍（陈仲新等，2000），固碳服 务价值约为773.5109～3 619.2109元（欧阳志云等， 1999）。后续研究进一步发现，生态系统的结构、状态 和动态可以影响系统的功能和稳定性（Elmqvist et al.， 2003；Cortina et al.，2006），从而影响生态系统固碳服 务价值，因而不同生态系统类型间差异很大（谢高地 等，2015；潘勇军等，2018）。然而，当前研究对于地区 生态系统结构变化的影响考虑较少，而这却是生态脆 弱与退化地区进行生态修复和国土空间规划布局的 基础。基于地区生态系统结构进行固碳功能评估，对 于脆弱生态的修复措施和经济建设的政策调整有重 要意义。 黄河流域的生态保护和高质量发展是重大国家 战略。榆林地处黄河中游地区，丰富的资源储量使榆 林具有明显的区位优势。然而，经济发展过程中的环 境效应问题也逐渐成为限制榆林市经济与生态可持 续发展的主要影响因素。生态系统的脆弱性加剧了 环境保护难度，在当前全球变化背景下如何改善生态 环境的同时增加应对环境压力和干扰的能力是榆林 地区高质量发展的必然要求。笔者基于高分辨率遥 感影像对榆林市生态系统固碳价值进行评估，分析环 境变化和人类活动的影响，结合景观结构指数对固碳 功能进行空间分析，尝试从系统的角度理解榆林市生 态系统固碳功能变化特征，以期为榆林市生态建设和 高质量发展提供科学依据。 1 研究方法与数据 1.1 研究区概况 榆林市位于黄河中游西岸，地理坐标为 E 10714′51″～11114′31″，N3648′58″～3935′7″，包 含1市2区9县，总面积为42 920 km2。全市地势西北 高、东西低，古长城以北的风沙草滩区占总面积的 42，风蚀沙化较严重，后备资源开发潜力大。南部黄 土高原丘陵沟壑区占58，地形破碎、水土流失严重。 因为地处农牧交错区，榆林在相当长的时间一直以农 牧业为主，从20世纪末开始发展轻重工业，到2006年 开始以能化产业为主导，榆林经济进入了快速发展阶 段。富集的能源矿产资源开发引发了对土地、植被资 源的破坏，生态环境质量评价较差，局部地区甚至仍 在恶化（陕西省生态环境状况公报，2021）。 1.2 数据来源 遥感数据采用MOD17A3HGF数据（https//ladsweb. modaps.eosdis.nasa.gov/），图像空间分辨率为500 m 500 m，数据选用2000、2005、2010、2015和2020年。 土地利用数据来源于资源环境科学与数据中心 （https//www.resdc.cn/），空间分辨率为1 km1 km，包 括耕地、林地、草地、水域、建设用地和未利用土地6 个一级分类和25个二级分类。 1.3 研究方法 1.3.1 固碳服务功能价值 MOD17A3HGF数据是基于BIOME-BGC模型通 过Terra卫星遥感参数计算的全球植被净初级生产力 年际遥感监测产品。将MOD17A3HGF数据转化为 NPP值，利用Modis Reprojection Tool对数据进行投影 转换，生成1 km1 km的栅格数据。将生成数据与土 地利用遥感监测数据和榆林市行政区划进行重采样、 重分类、栅格计算等处理，得到2000～2020年榆林市 6大类生产性土地的NPP统计数据。基于净生态系 统生产力法对固碳功能量进行估算（国家市场监督管 理总局，2020），采用市场价值法对固碳价值进行换算。 CO2吸收价值根据2022年中国碳交易市场价格取58 元/t。区县单位面积固碳价值通过各区县固碳价值量 的总和与区划总面积换算得到，各区县面积加权平均 52 西 北 地 质 NORTHWESTERN GEOLOGY 2023 年 得到全市单位面积固碳价值。 1.3.2 林地系统综合指数 基于遥感数据通过Fragstats对景观格局进行计 算。分别在景观和类型水平下选取斑块面积（CA）、 斑块所占景观面积比例（PLAND）、斑块数量（NP）、 斑块密度（PD）、最大斑块占景观面积比（LPI）、景观 形状指数（LSI）、平均斑块面积（AREA_MN）、蔓延度 指数（CONTAG）、斑块内聚力指数（COHESION）、聚 集度指数（AI）等指标，对各类土地利用类型的景观指 数进行计算。对各区县森林景观计算结果标准化处 理后通过熵权法确定指标权重，计算得到各区县林地 综合指数。对各年份分别计算，得到各区县林地综合 指数趋势变化。 1.3.3 主成分分析 基于气象条件、土地利用及工业GDP占比对各 区县固碳价值进行PCA分析（Principal Component Analysis），分析环境因素、生态系统结构和人类活动 对固碳价值的影响。采矿是榆林地区最主要的重工 业产业，由于缺乏直接的采矿数据，以工业GDP占比 代表地区工业活动。数据分析使用R-software （version- 3.5.3）。 2 分析结果 榆林市2000年陆地生态系统固碳价值为5 423 万元， 2020年达到11 449万元。西部的榆阳区、神木 市、靖边县和定边县固碳价值总量最高，平均占全市 10以上。东南部的绥德县、米脂县、佳县、清涧县 和子洲县固碳价值量占比较低，吴堡县最低仅约占 1。2020年榆林市固碳价值总量比2000年增加了 111，但各区县间增幅不同。西部的榆阳区、神木市、 定边县增幅低于全市，分别增长了89和96，靖边 县略高增加了114。东南部增幅高于全市水平，增 幅基本高于120，清涧县增加了175。 从土地利用类型看，榆林全市范围以草地和耕地 为主，面积占比82。2000年至2020年间全市草地 和耕地固碳价值量分别增加了116和101。林地 在榆林地区占比较小，然而到2020年固碳价值量增 加了176。建设用地增长速率最快固碳价值变幅最 大。全市范围未利用土地持续缓慢减少，固碳价值量 增幅最小。 3 讨论 3.1 榆林市固碳价值的空间分布与影响因素 榆林市固碳价值总量的空间分布西北高、东南低。 因为价值总量估算是基于行政区划面积，对区划面积 较小的区县会低估，笔者基于单位面积固碳价值对各 区县生态系统固碳功能进行对比。榆林市单位面积 固碳价值量多年平均21.25 元/ha，东南部区县单位面 积固碳价值水平高于全市，多年平均为21.96～24.98 元/ha，西北部区县基本低于全市平均水平（图1）。文 中估算结果低于已报道的陕北榆林地区结果（李晶等， 2011；封建民等，2018），可能与数据分辨率以及采用 的市场交易价格低于造林成本价格有关。榆林市单 位面积固碳价值从东南向西北递减的空间分布趋势， 与陕北地区NPP从南向西北递减的趋势一致（李登科 等，2022），显示出区域间土地利用类型的结构差异。 对12区县固碳价值的PCA运算显示东南部区县 与西北部区县在PC1（贡献率为44.1～57.5）上表 榆阳区 神木市 府谷县 横山区 靖边县定边县绥德县 米脂县 佳县 吴堡县 清涧县 子洲县 107 108 109 110 111 37 38 39 40 纬度 经度 10 13 16 19 22 25 28 31 34 单位面积固碳价值b 榆阳区 神木市 府谷县 横山区 靖边县 定边县绥德县 米脂县 佳县 吴堡县 清涧县 子洲县 107 108 109 110 111 37 38 39 40 纬度 经度 90 290 490 690 890 1 090 1 290 1 490 固碳价值a a. 多年平均固碳价值量（104元）；b. 多年平均单位面积固碳价值量（元/ ha） 图1 榆林市生态系统固碳服务价值空间分布 Fig. 1 Spatial distribution of value of carbon sequestration service in Yulin ecosystem 第 3 期 张凡等榆林市生态系统固碳价值对林地系统结构的响应变化 53 现出显著差异（图2a），西北部区县正向分布，东南部 区县负向分布。20年间工业GDP占比对PC1的贡献 率始终较高，建筑用地的相关性也逐渐增加，表明榆 林市生态系统固碳价值主要依据工业活动形成东南 部区县和西北部区县的分划，这与榆林市矿产资源集 中分布在市域西北部吻合。西北部区县矿产资源富 集，采矿活动等重工业发达，对地表的破坏可以直接 影响生态系统固碳服务功能。神木市是榆林典型的 采矿工业区，煤矿资源丰富，工业GDP占比在20年间 从86增加到98，对PC1的负荷也从2.09增加到 3.18，表明该区域生态系统固碳价值与工业、建设活 动的相关性在逐渐增加。榆林市西北部区县如榆阳 区、府谷县也都表现出受工业建设活动影响增加的相 似趋势（图2b）。气象条件对PC2轴的贡献率较高。 可能受退耕还林还草政策导致的土地利用结构调整 影响，各区县没有表现出显著差异。 −4 −3 −2 −1 0 1 2 3 4−4 −3 −2 −1 0 1 2 3 4 PC 2 PC1 2 000 2 005 2 01 0 2 02 0 a −8 −6 −4 −2 0 2 4 6 8 −4 −2 0 2 4 P/T 工业GDP 占比 耕地 林地草地 建筑用地 未利用地 榆阳 神木 府谷 横山 靖边 定边 绥德 米脂 佳县 吴堡 清涧 子洲 −0.6 − 0.4 − 0.2 0.0 0.2 0.4 0.6 −1 .0 −0.5 0.0 0.5 1 .0 PC 2 21 .6 PC1 53 .8 2 02 0b a.各区县2000、2005、2015与2020年的PCA排序，实心点为西北部区县分组， 空心点为东南部区县分组；b.各区县2020年主分量分布 图2 榆林市12区县的PCA排序 Fig. 2 The PCA distribution figure of 12 counties of Yulin 榆林市地貌类型南北不同，西北部是毛乌素沙漠 南缘风沙草滩区，南部是黄土高原丘陵沟壑区。全市 土地利用类型基本以草地和耕地为主，贡献了全市约 83的固碳价值。尽管耕草地比例在各区县不同，两 种土地利用类型的固碳价值总量并不存在显著差异。 东南部区县的林地面积比例较高，因而林地固碳贡献 在东南部区县相应增加。西北部靖边县、榆阳区、神 木市、府谷县、横山县、定边县3～4的林地固碳 贡献，在东南部区县增加到10以上，绥德县和清涧 县可以达到17，提升了东南部区县的单位面积固碳 价值。单位面积固碳价值与林地比例的显著正相关 （图3），以及与建设用地及工业活动的负相关表明，在 采矿、工业建设等人类活动干扰会损害地区生态系统 固碳功能，生态系统结构尤其是林地比例对固碳价值 的提高有重要影响。 全市范围看生态系统固碳价值的增长变化也表 现出明显的空间差异。全市单位面积固碳价值在 2000～2020年间持续增长，然而区域间增长速率不同 （图4a）。2020年西北部6区县固碳价值总量增加了4 357 万元，平均增长率了107，平均增长速率为0.67。东 南部6区县虽然行政区划面积小，固碳价值总量仅增 加1 668万元，然而区域内林地系统的面积比例、平均 斑块面积更高（图4b），增长速率也显著高于西北部区 县（p＜0.05）。林地系统较高的比例和增长率使东南 部区县单位面积固碳价值平均增加了147， 0.955的 增速显著高于西北部区县（p＜0.01）。 榆林地处中国干旱半干旱地区，黄土高原与沙漠 过渡地带，脆弱的生态条件导致植被生长敏感易变。 全球变化背景下，榆林地区也表现出气候暖干化和水 资源减少趋势（杨文峰等，2002；王德丽等，2011），多 坡地的地貌结构可能进一步加深地表干旱。矿产开 发对土地资源的破坏如果不加以维护管理，容易受到 人类活动干扰而退化，从而降低生态系统的稳定和功 能。林地综合结构指数显示榆林市东南部区县相对 较高（图4a）。这种差异不仅仅因为林地总面积比例 在东南部区县相对较高，而是区域内林地类型密度、 54 西 北 地 质 NORTHWESTERN GEOLOGY 2023 年 大小、形状、临近度、连通性等的综合表现，表明森 林生态系统的结构性在东南部区县相对较好，合理的 生态修复和管理措施必然能够更好的发挥系统功能 性。西北部区县林地系统结构性相对较差，不仅受地 理条件限制，也与地方产业有关。榆林市作为中国能 源化工基地，大量的矿产资源集中分布在西北部区县， 矿山开发挤占了资源并造成景观破碎化，不利于植被 发展。神木市林地类型比例、平均面积大小、连通性 在全市基本处于最低值，且在不断恶化，林地综合指 数较低并具有较大波动。矿山开发对林地生态系统 结构的破坏不仅影响了当前的固碳功能，而且可能会 对系统的稳定性及水土保持等生态系统服务产生影 响。由这种潜在影响得知，在易发生土壤侵蚀或地质 灾害的黄土高原地区（孙萍萍等，2022；仲原等，2022） 的生态建设，应尽可能加强系统的整体性、结构性， 协调好地区土地利用规划，避免景观的破碎化，加强 水资源利用成效，从而实现地区长期的“双碳” 目标。 工业 GDP 占比 年均温 年降水 未利用地比例 建设用地比例 水域比例 草地比例 林地比例 耕地比例 固碳价值1.0 0.8 −0.7 −0.6 0.6 0.6 −0.6 1.0 −0.7 −0.7 −0.7 −0.8 1.0 −0.6 0.6 −0.6 1.0 0.8 0.8 −0.6 0.8 1.0 1.0 −0.7 −0.6 0.9 1.0 −0.6 0.6 1.0 0.5 −0.7 1.0 −0.6 1.0 −1.0 −0.8 −0.6 −0.4 −0.2 0.0 0.2 0.4 0.6 0.8 1.0固碳价值 耕地比例 林地比例 草地比例 建设用地比例 未利用地比例 年降水 年均温 工业 GDP 占比 水域比例 图中数字为相关系数（p˂0.05）；代表无显著相关；固碳价值为单位面积固碳价值；土地类型比例为面积比例 图3 榆林市固碳价值影响因素相关性分析结果 Fig. 3 Correlation analysis results of influencing factors of carbon sequestration value in Yulin 107 108 109 110 111 37 38 39 107 108 109 110 111 37 38 39 综合指数 固碳价值 增长速率 0.5 0.8 1.1 1.4 0.51.0a PLAND AREA_MN 20 40 60 80 100 120 10 20b 经度 纬度 经度 纬度 榆阳区 神木市 府谷县 横山区 靖边县定边县绥德县 米脂县 佳县 吴堡县 清涧县 子洲县 榆阳区 神木市 府谷县 横山区 靖边县 定边县绥德县 米脂县 佳县 吴堡县 清涧县 子洲县 a.圆代表榆林市森林类型综合指数，颜色代表单位面积固碳价值增长速率； b.圆代表森林景观类型面积比例，颜色代表森林类型平均斑块面积 图4 榆林市生态系统固碳价值增长速率与森林系统景观结构综合指数 Fig. 4 The growth rate of carbon sequestration value in Yulin and forest ecosystem landscape structure comprehensive index 第 3 期 张凡等榆林市生态系统固碳价值对林地系统结构的响应变化 55 3.2 榆林固碳价值的年代变化 2000～2020年榆林市固碳价值持续稳定增长，然 而2000～2010年增长速率显著高于2010～2020年 （图5）。定边县和绥德县位于榆林市南部，两县都没 有煤矿开采，工业GDP占比不高，2020年分别仅占全 县GDP总值的10和15。两县均是以耕地为主的 传统农业县，建筑用地很少，土地利用类型的差异主 要是林地和未利用地（图5b）。林地较多的绥德县单 位面积固碳价值更高（图5b），增长率也更高（图5a）。 为了在更小尺度进行对比，分别在两县基于相似土地 利用结构选取了10 km 10 km样方，绥德县样方不仅 固碳价值提高了27.5，20年间固碳价值的增长率也 提高了92。 2000 2005 2010 2015 20200 5 10 15 20 25 30 时间 a 固碳价值 元 /ha 绥德县 定边县 全市平均 a 耕地 林地 草地 水域 建筑用地 未利用土地 绥德县定边县 b a.榆林市生态系统单位面积固碳价值随时间变化，虚线为12区县 平均值，阴影为标准偏差，红线为绥德县，蓝线为定边县；b.定 边县与绥德县土地利用结构 图5 榆林市生态系统单位面积固碳价值动态变化 Fig. 5 The dynamic change of value of carbon sequestration of Yulin 1999年，陕西在全国率先开展了退耕还林还草试 点，榆林市重点实施了建设工程，并从2006年起陆续 出台了一系列土地利用规划，加强土地整治与生态环 境建设。全市耕地面积在2000～2010年下降较快 （p＜0.01），2010～2020年下降速率减缓。20年间草 地和林地整体呈上升趋势，然而草地面积在2010～ 2020年出现下降，只有林地始终保持增加。耕地和草 地面积变化的年代差异达到显著水平，林地增长达到 极显著水平（p＜0.01）。土地整治政策的实施改善了 榆林市的生态环境，自然植被的逐渐恢复有效促进了 生态系统固碳功能的提高。2000～2010年榆林市林 地固碳价值增加337.66万元，2010～2020年增加 116.13万元。相应各区县固碳价值总量增加量和增 长速率均表现出极显著年代差异（p＜0.01）。西北部 区县受采矿等工业建设活动影响，固碳价值增长速率 普遍较低。加强全域统筹规划，因地制宜提高生态系 统结构的协调性，对榆林市西北部地区的生态修复与 建设有重要意义。 4 结论 （1）2000～2020年榆林市生态系统固碳服务价值 持续增长，2000年代固碳价值年增长480万元/a， 2010年代增速减缓年增长120万元/a。全市单位面积 固碳价值空间分布呈现从东南向西北递减的趋势，东 南部区县单位面积固碳价值增长速率也更快。 （2）采矿工业活动是形成榆林市西北、东南各区 县分划的主要因素，东南部工业活动少、林地占比较 高，固碳价值主要受生态系统结构性变化的影响。 （3）榆林地区长期的资源开发导致地表景观的破 碎化，降低了系统的稳定性。长期来看，未来实施退 耕还林还草与生态修复建设工程时，应尽可能扩大单 个系统类型面积、避免分散，综合考虑生态系统的整 体性，以改善林地系统的结构性与稳定性，实现提质 增效。后续可进一步对生态系统服务价值总量进行 核算，深入分析系统结构变化影响下的系统功能变化。 参考文献（References） 陈仲新，张新时. 中国生态系统效益的价值[J]. 科学通报，2000， 45（1）17–22. CHEN Zhongxin, ZHANG Xinshi. The value of ecosystem benefits in China[J]. Chinese Science Bulletin，2000，45（1）17–22. 封建民，文琦，郭玲霞. 风沙过渡区土地利用变化对生态系统服 务价值的影响以榆林市为例[J]. 水土保持研究，2018， 25（4）304–308. FENG Jianmin, WEN Qi, GUO Lingxia. Effects of land use change on the ecosystem service value in ecotone of wind and sand[J]. Research of Soil and Water Conservation， 2018， 25（4） 304–308. 国家市场监督管理总局, 国家标准化管理委员会. 生态系统评 估、生态系统生产总值GEP核算技术规范[S]. 国家市场 监督管理总局, 国家标准化管理委员会, 2020. 李登科，王钊. 气候变化和人类活动对陕西省植被NPP影响的 定量分析[J]. 生态环境学报，2022，31（6）1071–1079. LI Dengke, WANG Zhao. Quantitative analysis of the impact of cli- mate change and human activities on vegetation NPP in Shaanxi Province[J]. Ecology and Environmental Sciences，2022，31（6） 56 西 北 地 质 NORTHWESTERN GEOLOGY 2023 年 1071–1079. 李晶，任志远. 基于GIS的陕北黄土高原土地生态系统固碳释 氧价值评价[J]. 中国农业科学，2011，44（14）2943–2950. LI Jing, REN Zhiyuan. Research on the values of CO2 fixation and O2 release by landuse ecosystem in Loess Plateau in Northern Shaanxi province[J]. Scientia Agricultura Sinica，2011，44（14） 2943–2950. 欧阳志云，王效科，苗鸿. 中国陆地生态系统服务功能及其生态 经济价值的初步研究[J]. 生态学报，1999，19（5）607–613. OUYANG Zhiyun, WANG Xiaoke, MIAO Hong. A primary study on Chinese terrestrial ecosystem services and their ecological- economic values[J]. Acta Ecologica Sinica， 1999， 19（5） 607–613. 潘勇军，王兵，牛香. 让自然资本成为主流，用生态GDP核算美 丽中国. 温带林业研究[J]. 温带林业研究，2018，1（3） 10–18. PAN Yongjun, WANG Bin, NIU Xiang. Mainstreaming the econom- ics of natural capital and accounting the eco-GDP of beautiful China[J]. Journal of Temperate Forestry Research，2018，1（3） 10–18. 陕西省生态环境厅. 2021年陕西省生态环境状况公报[R]. 陕西 省生态环境厅, 2022. 孙萍萍，张茂省，贾俊，等. 中国西部黄土区地质灾害调查研究 进展[J]. 西北地质，2022，55（3）96–107. SUN Pingping, ZHANG Maosheng, JIA Jun, et al. Geo-hazards re- search and investigation in the Loess Regions of Western China[J]. Northwestern Geology，2022，55（3）96–107. 王德丽，殷淑燕，王海燕，等. 近50年陕北地区的气候变化[J]. 干旱区研究，2011，28（2）262–267. WANG Deli, YIN Shuyan, WANG Haiyan, et al. Analysis on cli- mate change in North Shaanxi Province[J]. Arid Zone Research， 2011，28（2）262–267. 谢高地，张彩霞，张雷明，等. 基于单位面积价值当量因子的生 态系统服务价值化方法改进[J]. 自然资源学报，2015， 30（8）1243–1254. XIE Gaodi, ZHANG Caixia, ZHANG Leiming, et al. Improvement of the evaluation method for ecosystem service value based on per unit area[J]. Journal of Natural Resources， 2015， 30（8） 1243–1254. 杨文峰，李星敏，鲁渊平. 陕西省近50年气候变化及其对水资 源的影响[J]. 陕西气象，2002，61–6. YANG Wenfeng, LI Xingmin, LU Yuanping. YANG Wenfeng, LI Xingmin, LU Yuanping[J]. Journal of Shaanxi Meteorology， 2002，61–6. 仲原，王春梅，庞国伟，等. 黄土高原浅沟空间分异特征[J]. 生 态学报，2022，42（5）1798–1810. ZHONG Yuan, WANG Chunmei, PANG Guowei, et al. Spatial vari- ation of ephemeral gully in the Loess Plateau[J]. Acta Ecolo- gica Sinica，2022，42（5）1798–1810. Cortina J, Maestre FT, Vallejo R, et al. Ecosystem structure, function, and restoration success are the related[J]. Journal for Nature Conservation，2006，14152–160. Costanza R, dArge R, de Groot R, et al. The value of the world’s ecosystem services and natural capital[J]. Nature，1997，387 253–260. Elmqvist T, Folke C, Nystrm M, et al. Response diversity, ecosys- tem change, and resilience[J]. The Ecological Society of Amer- ica，2003，1（9）488–494. 第 3 期 张凡等榆林市生态系统固碳价值对林地系统结构的响应变化 57