土地利用/土地覆盖变化(LUCC)是全球气候变化和全球环境变化研究中的重要内容, 是人类社会经济活动行为与自然生态过程交互衔接的纽带^[1-2]。土地利用变化的同时也引起生态系统服务功能和结构的变化, 从而带来生态系统服务价值的改变^[3]。生态系统服务是通过生态系统的功能直接或间接得到的产品和服务, 是由自然资本的能流、物流、信息流构成的生态系统服务和非自然资本结合在一起所产生的人类福利^[4-6]。1997年Costanza等^[5]在Nature上发表的关于全球生态系统服务价值评估的文章, 生态系统服务开始成为生态学的热点研究。国内外学者先后对Costanza提出的评估方法进行相关的研究, Dawson N.等发现由于生态系统服务理论未对社会复杂性进行考虑, 使得相关研究不能为实际社会中的生态管理提供有效的指导方案^[7]。Hassan M.认为生态系统本身存在异质性, 不同地区使用同一套生态系统服务价值系数的做法, 其结果的准确性将会受到质疑^[8]。因此, 我国学者谢高地等^[9]在Costanza理论基础上针对中国的实际情况对其进行系数修正, 得到了中国陆地生态系统服务价值当量表。当量因子法与常见的条件价值法、影子工程法、市场机会法和资产价值法等相比具有使用简便、数据需求少、结果可比性高、评估全面等优点, 该方法被较多研究人员用于区域的生态系统服务价值评估^[10-15]。

赣南地区是典型的南方山地丘陵区, 属于赣江流域重要的生态屏障区, 该区林业资源丰富, 森林覆盖率高, 其丰富的林业资源和保护良好的原生生态系统对区域气候调节和生态平衡至关重要。生态系统服务价值与社会经济的发展状况及区域的自然地理要素分布密切相关, 生态系统服务空间格局是人类针对不同区域特征实施生态系统服务功能保护和利用的重要依据。因此本文从区域尺度上对生态系统服务价值格局的时空变化进行量化和描述, 并且利用SPSS和GeoDa软件对生态系统服务价值时空变化的驱动因素进行多元回归分析和空间相关性分析, 进一步为科学开展丘陵生态系统服务空间格局研究、生态系统的管理和资源的合理有效配置、生态环境保护、制定社会经济发展政策提供重要依据, 对促进该区域生态环境改善和生态环境的可持续发展, 人与自然的和谐发展具有重要意义。

1 研究区概况

赣南地区位于江西省南部, 赣江上游, 包括赣州市全区的18个县市, 全区总面积3934700 hm², 占江西总面积的23.6%, 为江西省最大的行政区, 地处24°29′—27°09′N, 113°54′—116°38′E之间(图 1)。全区平均海拔高度在300—500 m之间, 山地、丘陵占总面积的80.98%, 是典型的山地丘陵区。该区多年平均气温19.3℃, 多年平均降水量1568.75 mm, 属于亚热带山区湿润季风气候, 水热条件好, 气候条件优越。其林业资源丰富, 森林覆盖率达到74.54%, 是我国商品林基地和重点开发的林区之一。该区丰富的林地资源和保护良好的原生生态系统, 对全省的生态系统服务的价值有较大的贡献, 但其耕地面积小, 耕地后备资源不足, 且该区农业生产条件较差, 经济技术水平落后, 城镇化水平较低。故对该区域的生态系统服务价值的空间格局进行研究, 有助于优化该区的产业分布, 在不破坏生态环境的前提下, 实现该区域的生态系统服务最优策略, 从而促进当地经济的可持续发展。

2 研究方法 2.1 数据来源及处理

本研究采用1990年、1995年、2000年、2005年、2010年和2015年6个时段的Landsat TM/ETM遥感影像, 结合研究区域1:25万地形图、土壤和植被专题图进行分析。利用面向对象的遥感分类方法, 结合区域土地类型特征, 将赣南地区土地利用类型分为两级(图 2), 一级用地分为6类：耕地、林地、草地、水体、建设用地和未利用地, 二级分类中把林地分为针叶林、阔叶林、针阔混交林和灌木林, 耕地分为水田和旱地, 六期遥感影像解译精度均达到80%以上。

2.2 土地利用动态度

土地利用动态度是评价土地利用类型数量变化、土地利用程度变化及土地利用变化的区域差异的指标^[16]。本文计算了赣南地区的单一土地动态度和综合土地动态度, 以分析1990—2015年赣南地区的土地利用动态。

综合土地利用类型动态度公式:

式中, LU_i监测起始时第i类土地利用类型的面积；ΔLU_i-j为监测时段内第i类土地利用类型转为非i类土地利用类型面积的绝对值；T为监测时段长度。当T的时段设定为年时, LC的值就是该研究区土地利用年变化率。

2.3 生态系统服务价值评估

本文参照近年来关于生态系统价值的研究^[17-22], 结合谢高地等^[9]提出的价值当量换算方法确定赣南地区生态系统服务价值当量表。谢高地等将生态系统生产的净利润看作该生态系统所能提供的生产价值, 将单位面积农田生态系统粮食生产的净利润当作1个标准当量因子的生态系统服务价值量。确定中国2010年1个生态系统价值当量的经济价值为3406.50元/hm², 根据《2010年赣州市统计年鉴》(赣州市统计局, 2011)资料, 赣州市单位面积粮食产量为5345 kg/hm², 而同期全国单位面积粮食产量为4974 kg/hm², 以此对各生态系统服务当量价值进行系数修订, 得出研究区的修正系数为1.075, 即赣南地区的一个生态系统服务价值当量为3528.79元/hm²。根据修正后的价值当量和各用地类型的面积得到赣南地区生态系统服务价值。生态系统服务功能价值的评估公式为：

式中, ESV表示生态系统的总服务价值；A_i表示第i种土地利用类型的面积(hm²)；VC_i表示第i种土地利用类型的生态系统服务价值系数；ESV_f表示生态系统第f项服务功能的价值；VC_fi表示第i种土地利用类型的第f项生态系统服务价值系数。

2.4 驱动因素分析

生态系统服务价值时空差异的驱动因素包括自然因素和人为因素两大类^[23]。根据相关研究^[10-11]可知, 生态系统服务价值在短时间内主要受到人为因素的干扰。本文参考其他学者已有研究并结合数据可获取性, 选取人口、经济和旅游3个指标的21个人为因素在SPSS 19中进行主成分分析及相关性分析, 筛选出对区域生态系统服务价值时空变化差异影响较大的驱动因素进行多元线性回归分析, 得到线性回归方程, 从而得到生态系统服务价值时空变化的主要驱动因素。

2.5 空间自相关分析

根据驱动因素分析的结果选取人均GDP、人口密度和城市化率3个指标进行空间相关性分析。利用GeoDa模型中的Bivariate Local Moran′s Ⅰ分析赣南地区ESV与人口密度、人均GDP和城镇化率相关性, 绘制BiLISA集聚图, 把集聚区域划分为：不显著、高-高型、低-低型、低-高型和高-低型^[24-25]。

3 结果与分析 3.1 赣南地区土地利用类型空间变化特征

1990—2015年林地一直是赣南地区的主要土地利用类型, 占总面积的比例分别为75.248%, 75.266%, 75.175%, 74.969%, 75.055%, 74.542%, 处于生态系统的主导地位。耕地为面积比重第二高的土地利用类型, 主要分布在赣南地区东北部的宁都县, 于都县, 兴国县及西北部的南康市和信丰等县。赣南地区的草地、水体、建设用地和未利用地所占面积比重较少。

根据土地利用动态评估结果(表 1)可知, 1990—2015年, 赣南地区土地利用综合动态度呈上升趋势, 但整体变化幅度较小。2010—2015年, 土地利用动态度达到最大值, 2005—2010年最小。这说明2010—2015年土地利用最为活跃, 而2005—2010年土地利用最不活跃。林地在1990—2005年间一直呈现下降的趋势, 在2005—2010年间有所上升, 2010年之后又呈现下降的趋势。与林地变化趋势恰好相反, 耕地面积呈增加趋势, 由1990年的670800 hm², 上升到2015年的680900 hm², 净增加11100 hm², 尤其是大余县、寻乌县、赣县的耕地面积呈现迅速上升趋势。在2005年之前, 赣南大多数地区的耕地面积都在增加, 只有少部分呈现下降趋。2005年开始, 除大余县之外, 大部分县都较好地执行国家退耕还林政策, 耕地面积基本保持下降或者平衡状态。草地面积在2010年以前持续减少, 在2010—2015年期间出现上升, 整体呈略微上升的趋势, 面积的变化也表明了2010年以后退耕还林还草政策的作用更加显现。2000年后赣南地区开始步入快速城镇化发展期, 随着经济发展, 大量农民外出务工, 部分耕作条件差的农田被撂荒变为草地, 这也是导致耕地面积减少的部分原因。赣南地区水体的面积基本保持稳定状态。建设用地主要分布在市区和县城, 一直呈快速增加之势, 其面积由1990年的33000 hm²快速增加到2015年的47900 hm², 净增14900 hm², 为1990年的1.5倍。

3.2 赣南地区生态系统服务价值时空变化 3.2.1 生态系统服务价值时间维度的差异性分析

本研究以5年为一个时间尺度, 估算了基于土地利用二级分类的1990—2015年赣南地区的生态系统服务价值量(表 2)。1990—2015年赣南地区的生态系统服务价值量分别为2526.48、2519.66、2516.93、2508.86、2510.63、2506.06亿元, 整体呈现下降的趋势。赣南地区最突出的生态系统服务为调节服务, 占总服务价值的68.50%, 其次依次为支持服务, 供给服务和文化服务。

调节服务主要包括气体调节、气候调节、净化环境和水文调节4个方面的内容。赣南地区气体调节、气候调节和净化环境服务功能的变化基本一致, 呈下降趋势, 在1990—1995年及2005—2010年有略微的回升, 其他阶段均为下降趋势。这3种调节功能主要受耕地、林地和水体的影响较大, 由于赣南地区林地的面积比重很大, 所以调节功能变化与林地的面积变化基本保持一致。水文调节整体呈现出先下降后上升再下降的过程, 该功能受水体的影响较大, 水文调节价值的变化与水体的面积变化一致, 均呈现出下降, 上升再平稳下降的趋势。

赣南地区第二大服务功能为支持服务, 该功能包括土壤保持功能、生物多样性功能和维持养分循环功能, 土壤保持功能和生物多样性功能的变化趋势表现一致, 在2005年之前呈现下降趋势, 之后略有回升, 2010年后再次下降。维持养分循环功能的价值则呈现波动下降状态, 尤其是2010年后下降迅速, 原因是该功能除了受草地和林地的影响之外, 还受到了耕地的较大影响。耕地在2010年前稳步上升, 耕地的上升弥补了林地和草地变化带来的影响, 这也是维持养分循环功能在2010年之前的价值变化相比另外两个功能的价值变化较为平稳的主要原因, 而2010年后林地和耕地都呈现较大的下降, 所以维持养分功能下降最为剧烈。

供给服务和文化服务比重较小且在研究期内变化不大。供给服务中的食物生产功能主要受到耕地的影响, 故赣南地区食物生产功能价值的变化与耕地面积变化保持一致, 呈现出缓慢增加, 快速增加, 缓慢减少和快速减少四个阶段。原料生产受林地、草地和旱地的影响都较大, 该区的原料供给功能价值变化与林地变化呈现较为同步的趋势, 主要还是由于林地面积比重大, 对其价值变化的引导起到了主要的作用。水资源供给受到水体和水田的影响较大, 水体对于水资源供给有很大的贡献, 而水田在很大程度上削弱了水资源服务的价值。赣南地区供给服务的主要贡献者是食物生产和原料生产, 水资源供给价值较小, 其变化对总供给服务的影响几乎为零, 故其变化主要受到林地、耕地和草地3个用地类型综合影响, 变化原因较为复杂, 整体呈减少的趋势。文化服务价值受林地和水体的影响较大, 基本处于下降, 这与林地的变化趋势较为一致。

整体上, 赣南地区的生态系统服务价值呈现下降的趋势, 在1990—2015年期间总价值减少, 总价值变化与调节服务的变化趋势一致。

3.2.2 生态系统服务价值空间维度的差异性分析

生态系统服务价值总量随着时间的变化呈减少趋势, 但总量的变化并不能反映空间上的变化, 本文采用单位面积(hm²)生态系统服务价值来分析各县生态系统服务价值(ESV)的空间变化(表 3, 图 3), 分析结果可知：崇义县的生态系统服务价值最高, 1990—2015年的生态系统服务价值均排在第一, 超过70000元/hm², 赣州市市区的生态系统服务价值最低, 1990—2015年均排在最后, 且处于持续下降状态。把生态系统服务价值变化分为Ⅰ型(增长型)和Ⅱ型(衰退型), 赣南地区18个县中有4个Ⅰ型和14个Ⅱ型, 在14个Ⅱ型的县中有8个县的生态系统服务价值下降的值均超过1000元/hm²。从图 3可知, 增长型主要分布在赣南地区的东部多个县市和西北部的崇义县, 其他地区生态系统服务价值处于下降的状态, 其中西南地区生态系统服务价值下降最为明显。

3.3 赣南地区生态系统服务价值时空差异驱动因素分析

通过相关性分析和主成分分析后选取人口、经济及旅游3个指标的10个因素进行多元线性回归分析, 得到生态系统服务价值与各因素之间的线性回归方程, 该线性方程的拟合度大于0.99, 拟合效果极好。

式中, Y：生态系统服务价值(亿元)；X₁：农业人口数(人)；X₂：农业总产值(万元)；X₃：接待入境旅游人数(万人)；X₄：城市化率(%)；X₅：人均GDP(元)。

通过多元回归分析可知, 赣南地区生态系统服务价值的整体变化受到以上5个因素的影响。本文中城市化率由城镇人口数占总人口数的比重求得, 根据回归方程的相关性系数可知, 城市化率越高, 生态系统服务价值越低。这是因为城镇人口的增加会导致区域的经济和生态结构变化, 从而导致生态系统服务价值发生变化。

3.4 赣南地区社会经济因素与生态系统服务价值空间相关性分析

根据多元回归分析结果结合数据的可获取性, 选取赣南各县市表征人口和经济的3个因素与生态系统服务价值变化率进行空间相关性分析。局部分析结果(图 4)显示：人均GDP、人口密度和城市化率与ESV的局部空间相关性具有高度的一致性, 主要集中分布在赣南地区东北部的会昌、于都和瑞金3个县市。空间相关性分析的显著性水平达到95%以上, 说明3个因素与ESV变化率的局部空间相关的显著性良好。这3个县市的ESV都呈现出较高的增长率, 其中于都县和瑞金市的人均GDP变化率和城市化变化率低, 人口密度变化较大, 故分别在局部形成低-高型和高-高型集聚。会昌县因为高ESV变化率和高城市化率变化率, 在城市化率变化率与ESV变化率的集聚上表现出与另外两个县不同的集聚模式。

4 讨论

本文在谢高地等^[9]于2015年更新的中国陆地生态系统服务价值当量表的基础上结合赣南地区的实际情况进行系数修正, 得到赣南地区生态系统服务价值当量表, 该法使用简单, 数据需求少, 对于区域的生态系统服务价值评估特别适用。近年来关于江西及赣江区域较多研究^{[23, 26-28]}均使用该法进行区域的生态系统服务价值评估, 因此运用该方法对赣南地区的生态系统服务价值进行评估方法合理, 得到结果较为准确及可靠, 且能够较好的反应出赣南地区生态系统服务价值的时空动态变化。生态系统服务价值时空变化的原因较多, 根据选取的指标进行相关性分析和主成分分析可以剔除大部分对其影响较小且相关性不够显著的因素, 最终选取10个因素进行多元回归分析, 得到生态系统服务价值的多元回归方程, 拟合度高, 驱动因素分析较为准确。并且文章在此基础上对关键驱动因素进行空间相关性分析, 能够更清楚赣南地区生态系统服务价值变化受人为因素影响较大的空间分布, 也有利于该区的生态系统管理。

从分析结果可知, 赣南地区的林地比例高, 但其在25年间的处于下降趋势, 呈现与全国同期的林地变化相反的趋势, 这与赣南地区原本就极高的森林覆盖率有极大的关联。但高的森林覆盖率并未赋予该区较高的生态系统服务价值, 这与其地形特点有极大的关系。山地丘陵地形限制了其农业的发展, 因此提高农业耕作技术和改善农业种植环境对于该区的生态系统服务价值提高极为重要。

5 结论

本文利用江西赣南地区6期遥感分类数据, 研究了25年间该区土地利用类型结构的变化, 并对生态用地所提供的生态系统服务功能的时空演变进行评估, 并分析了人均GDP变化率、人口密度变化率及城市化率变化率与ESV变化率的空间相关性。主要结论如下：

(1) 赣南地区的主要用地类型为林地, 根据土地利用动态度结果可知，1990—2015年, 赣南地区土地利用综合动态度呈上升趋势, 但整体变化幅度较小, 其中建设用地, 草地和耕地的单一动态度比其他四类用地单一动态度更大, 说明这3种类型的用地变化最为明显。

(2) 赣南地区调节服务价值为总价值的主要贡献者, 占总价值的68.49%, 因此赣南地区对于临近城市乃至全省起到了很好的调节作用。但在研究期内生态系统服务价值整体呈现下降趋势。

(3) 通过多元回归分析得出, 影响赣南地区生态系统服务价值时空差异的主要人为因素包括人均GDP、城市化率、农业总产值、农业人口数和接待入境旅游人数。通过GeoDa空间相关分析可知, 社会经济因素与ESV变化率的空间相关性主要集中的赣南东北地区三县, 在空间分布上具有极高的显著性。

(4) 林地对于赣南地区生态系统服务价值的维持起到关键性作用, 未来发展应注重保持现有的林地优势, 积极开展生态旅游观光等提升当地的经济效益。同时要提高农业生产技术, 以改善其山地丘陵地势所带来的农业生产薄弱的劣势。