生态可持续发展已成为备受各界关注的焦点^[1-2], 其一方面强调社会和经济的发展必须要尊重地球的生物物理极限^[3], 另一方面强调生态系统是不可分割的有机整体^[4-5]。故从整体的观点研究生态可持续性是必然趋势。但当前关于可持续性的研究多是从社会、经济、环境3个子系统中分别遴选评价指标, 并计算得出综合指数来反映某区域的可持续状态^[6]。为了与生态可持续性强调的生态系统有机性和整体性相适应, 很多学者倡导转变研究思路, 以生态系统为研究对象, 从它的基本特征出发——生态系统是一个开放的、具有自组织能力的生命体(即“生命体”观), 能够不断优化自身的结构和功能, 朝着有序的方向发展^[7-8]。

随着研究思路的转变, 熵这一重要分析方法受到关注。将熵用于可持续性问题的探索已有近三十年的历史。Addiscott^[9]在1995年初步阐述了熵视角下生态系统的可持续性过程, 并提出以“最小熵产生”作为衡量系统是否可持续的概念模型；Zurlini等^[10]认为生态可持续性可以通过衡量生态系统的功能或过程指标动态演变的有序性来表征；Nielsen等^[11]更是指出从熵的视角出发, 通常可以更好地对传统科学提出的基本问题(如可持续性、复原力、完整性和生态系统健康等)的因果关系进行解释。熵被应用于生态可持续性研究主要有以下3个方面的原因：一是在生态学领域, 将生态系统看作有活力的生命体, “生命赖负熵为生”的观念对生态系统研究影响深远^[12]；二是生态可持续性本质上是生态系统的结构和功能不退化的状态, 而熵被广泛用于度量系统的退化(或无序性)程度^[13-14], 反映复杂系统的自组织能力^{[11, 15]}；三是熵具有反应系统复杂性的整体性的特性^[16], 复杂性度量可以有效捕获生态系统状态的系统级变化^[17], 具有“透过树木以见森林”的优势^[5]。

熵在生态可持续性研究中的具体应用分为两大类——生态系统格局和生态系统动态。

关于生态系统格局, 如何定量区分和刻画不同的景观格局是生态学研究的重点之一^[18]。许多研究指出, 生态系统格局的空间分布及其演变规律遵循熵原理^[19-20]。熵度量可以对景观格局的均匀性和随机性进行分类, 是复杂性度量的良好替代^[17]。相对于传统景观指数理论, 基于熵方法表征生态系统格局在空间分布的有序或无序程度, 可以更好地深入生态过程的机理层面^[21-22]。如何选择合适的熵模型, 以及如何让熵方法在有效表达空间格局分布特征的同时具有生态学含义是当前学者们研究的重点和难点^[22]。

关于生态系统动态, 对长时间序列数据记录的生态系统变化过程进行客观的表达与分析, 可以反映生态系统的健康状况^[23]。受自然和人为双重耦合因素的影响, 生态参数的变化呈现出非稳定和非线性特点^[24]。同时, 生态系统是有自组织能力的生命体, 具有自调节、自适应的能力, 其变化除了对外界条件改变的响应, 还是其具有生命力的表现, 应同时考虑其自身的属性特征^[25]。与时间序列分割法、主成分分析法、傅里叶变换等常用方法相比, 熵方法从“生命体”的角度出发, 可以更好地表达生态系统复杂的变化信息, 揭示生态系统的自组织和弹性能力等自身属性特征^[26-27]。目前基于熵方法分析生态系统动态特征已得到较为广泛应用^{[23, 28]}, 部分熵方法存在的先验函数模型选择困难以及计算结果的生态学意义难以解释等问题, 使其应用范围受到一定约束。

此外, 已有基于熵方法的生态可持续性研究多从生态系统格局或动态分途讨论生态可持续性, 对二者之间关联的分析较少。但生态系统是一有机整体, 格局和动态之间存在着密切的关联, 将二者结合起来分析, 格局(空间)和动态(时间)相辅相成, 无疑可以提供更加客观的可持续性分析结果。

本文从熵的视角出发, 以延河流域为研究区, 尝试提出将格局和动态相结合的时空信息熵方法。延河流域是典型的黄土丘陵沟壑区, 其生态环境问题尤其是退耕还林还草工程实施过程中的植被恢复可持续性问题一直吸引着众多学者的关注^[29-31]。目前关于延河流域生态可持续性研究, 主要侧重具体地类的变化特征(如耕地向林、草地的转化), 以及由此带来的生态效益(如各地类所提供的生态服务功能), 而鲜有从整体的视角分析该流域生态系统格局和动态的变化特征。以时空信息熵方法定量分析该流域的生态可持续性, 希望为当地及黄土高原其他类似流域的生态保育和管理工作提供辅助决策支持。

1 研究区概况与数据来源 1.1 研究区概况

延河是黄河的一级支流, 全长约286.9 km, 发源于榆林市靖边县天赐湾乡, 自西北向东南依次流经靖边县、志丹县、安塞区、宝塔区和延长县, 最后注入黄河。延河流域位于陕北黄土高原腹地, 地理位置36°23′—37°17′N, 108°45′—110°28′E, 总面积约7687 km², 地势起伏较大, 呈现出西北高、东南低的特点, 平均海拔约1214 m。在地质构造上, 该区域黄土地貌沟壑纵横, 主要是第四纪的黄土岩系, 土质疏松, 抗冲刷、抗侵蚀能力差, 极易被分散和搬运。该流域属半干旱大陆性季风气候区, 年平均气温约6.2—12.4℃。年平均降雨量约为520 mm, 降水季节性分配不均匀, 6—9月份降雨量占全年总量的75%以上, 且多以暴雨形式出现。流域内植被多样性较为丰富, 植被分布从南向北呈现出梯度分布特征, 大致分为三大植被带分区：森林带、森林草原过渡带和草原带。延河流域位置及2015年土地利用情况见图 1。

过去受不合理的人类活动和自然因素的影响, 延河流域的生态系统明显退化, 天然植被锐减、水土流失等现象严重^[32]。自1999年起, 为了遏制该地区生态系统的退化, 国家开展了一系列生态建设政策, 退耕还林还草工程在延河流域大范围开展^{[29, 32-33]}。然而延河流域(乃至整个黄土高原)在该工程实施过程中, 由于所引进的植被类型的耗水量较高, 且未充分考虑当地土壤水分的承载力, 导致土壤水分过度消耗、土壤水分的涵养功能退化, 由此引发土壤干层, 反过来又影响到植被恢复的可持续性^[29-30]。

1.2 数据来源

本研究所使用的数据包括：

(1) 土地利用数据：生态系统格局多以土地利用/覆盖的形式表征, 土地利用的组成及其空间分布特征会直接影响生态系统功能^[34-35], 与生态可持续性密切相关。本文采用来源于中科院数据中心(http://www.resdc.cn)的三期(2000、2010、2015年)土地利用数据, 其空间分辨率为100 m, 包括15种土地利用类型, 数据集总体精度超过80%。

(2) 归一化植被指数数据：归一化植被指数(NDVI)是反映植被覆盖水平的指示因子, 且植被是陆地生态系统的重要组成部分, 可以有效反映生态环境状况^[23]。因此, 时间序列NDVI数据可以作为研究生态系统动态演变的关键状态变量。本文采用来源于NASA数据中心(https://ladsweb.nascom.nasa.gov/data)的MODIS植被指数产品(MOD13Q1 V6), 时间范围为2000—2018年, 共计432景数据, 空间分辨率为250 m, 时间分辨率为16天。经过投影、格式转换、数据重采样、裁剪、对缺失像元的数据修复等环节, 得到延河流域NDVI月均值数据集, 共计226组。进一步利用最大值合成法, 得到2000—2018年NDVI年极大值数据。

2 研究方法

将流域生态系统看成是一个开放的、具有自组织能力的生命体, 其可持续性可以通过刻画系统功能和过程指标动态演变的有序性来表征^[10]。利用空间信息熵和时间信息熵, 研究生态系统格局以及生态系统动态演变的特征, 并将两者结合分析生态系统可持续性, 初步建立结合“格局”和“动态”的时空信息熵方法。同时, 需要指出的是, 新方法的目的是探索生态可持续性分析的新视角, 可视为对已有研究方法的补充而非替代。

2.1 空间信息熵

玻尔兹曼熵和信息熵是研究生态系统格局常用的两种熵方法。基于玻尔兹曼熵的方法在实际应用中存在随着景观维度和类型数目的增加, 计算量剧增的问题, 且其景观组分和理想气体之间是否存在热力学联系的前提假设尚不明确。与玻尔兹曼熵相比, 信息熵由于其自身的“数学理论”特性, 在景观格局分析中有着更广泛的应用^[22]。然而由于信息熵的研究对象是一维的, 而空间是二维或多维的, 故信息熵难以反映实际空间格局的分布特征, 即存在“维度不匹配”问题。

空间信息熵(H_s)是一种可以有效解决该问题的景观格局分析方法^[22], 其基于信息熵, 结合地理学第一定律^[36], 以考虑公共边界长和质心距离的邻近度作为定量刻画空间属性影响的关键因子, 使得H_s可以同时包含格局的组成信息及结构信息, 更为客观地表征格局的空间异质性和无序性。H_s的计算式如下：

(1)

(2)

(3)

式中, n为土地利用类别的数量, p_i为土地利用类别i所占面积比例。L_i/d_i被定义为邻近度, 其中L_i表示不同类别之间的公共边界长, 它等于与类别i相邻的其他类别的所有边界长之和(公式(2))。d_i表示不同类别间的距离, 即类别i与其他不同类别质心间的距离之和(公式(3)), 本文所指距离为欧式距离, 为避免当不同类别的质心重合时, d_i为空值产生噪声效应, 需将此时的d_i设为相对较小的常量, 本文取为栅格长度的一半。由上式可知, H_s的取值范围为[0, +∞), 对于给定的类别数目和比例的不同空间格局, H_s会随着邻近度的变化而变化。

关于空间信息熵的生态学意义, 对于同一研究区域, 不同时期生态系统格局的熵值降低, 反映系统生态复杂度提升, 朝着有序的方向发展, 生态可持续状况好转^[22]。但需要补充的是, 由于信息熵具有概率属性^[37], 故熵值的下降并不绝对意味着系统状态的好转。为避免此种情况, 在实际研究中, 可选择具有一定时间间隔的生态系统格局数据, 将生态系统的自组织能力考虑进来。

2.2 时间信息熵

信息熵从生态系统为“生命体”的角度出发, 揭示生态系统动态演变的机制以及系统的组织能力^[11]。时间信息熵(temporal information entropy, H_t)^[27]是信息熵在反映生态参数在时间维动态演变相对不确定性的延伸, 其基于观测数据直接计算熵值, 克服了以往信息熵方法需要选择合适先验模型的局限性^[32], 适用范围更广。H_t的计算式如下：

(4)

其中,

式中, n为所选时间序列的项数, y₁ ≤ y₂≤ … ≤ y_n是像元不同时期的观测值x₁, …, x_n按照从小到大的顺序排列后得到的。m表示时间频率因子, 为不超过n/2的正整数, 该值的选择与所研究的时间范围相关。Δ为缩放系数, 使计算结果恒为正值, 一般取为时间序列数据的基本变化单元。

时间信息熵的生态学意义需要视具体情况而定, 需考虑生态系统所处的宏观状态, 如稳定态、生长期或恢复期、退化期等。并且由于H_t反映的是生态弹性增强的程度, 而非绝对值, 因此需要兼顾不同的景观格局(特别是不同的植被类型)。H_t的生态学含义通常与系统的弹性能力结合起来讨论, 弹性反映的是系统能够维持自身稳定性的能力^[38], 可以反应生态系统的可持续性。在稳定态的条件下, 生态参数的H_t值越小, 反映其动态变化越稳定、有序, 相应的生态弹性能力越强, 系统越可持续。而在生长期或恢复期条件下, 对于同一类型的植被, H_t值越大, 生态弹性越强, 生态系统越有序, 生态可持续性越强。

2.3 时空信息熵

基于时间信息熵和空间信息熵, 将“格局”和“动态”两个维度初步结合起来, 建立定量分析生态可持续性的时空信息熵方法, 生态系统格局的熵值(即空间信息熵)用于反映生态可持续性的整体状态, 提供关于系统所处宏观状态的客观信息, 为生态系统动态特征的分析提供参考依据。同时, 将时间序列生态参数动态演变的熵值结果(即时间信息熵)与土地利用格局(特别是不同植被类型)结合起来, 理解计算结果的生态学含义, 判断生态可持续性。

时空信息熵方法具体步骤如下：

步骤一：利用多期生态系统格局数据计算各时期的空间信息熵, 根据空间信息熵的变化情况, 分析生态系统格局的有序(或无序)程度, 判断生态可持续性的整体情况和生态系统所处的宏观状态(见2.1节)。本文以土地利用数据表征生态系统格局, 考虑到信息熵的概率属性以及生态系统的自组织能力, 可能需要一定的时间才可以反映出退耕还林还草工程对延河流域生态系统的影响, 因此选用2000、2010、2015年三期土地利用数据计算空间信息熵；

步骤二：选取合适的生态参数, 基于该参数的长时间序列数据, 计算时间信息熵, 分析生态系统动态演变状况(见2.2节)。本文选取NDVI数据作为研究生态系统动态演变的关键状态变量。由于本文的目的是探究生态系统在一段时间范围内的总体变化特征, 侧重对年际间动态特点的讨论。因此, 忽略了年内物候因素, 选取2000—2018年NDVI年极大值数据作为时间信息熵值计算的输入数据；

步骤三：由于不同生态系统宏观状态下, 时间信息熵所反映的生态学含义不同：在稳定态的条件下, H_t与生态可持续呈负相关, 即熵值越大, 生态可持续性越弱；而在生长期或恢复期条件下, H_t与生态可持续呈正相关。因此, 将步骤一得到的生态系统宏观状态作为前提条件, 根据步骤二所得时间信息熵值大小判断研究区各像元的生态可持续性。

3 结果与分析 3.1 延河流域生态格局分析

选取延河流域2000、2010、2015年三期土地利用数据进行基于空间信息熵的生态系统格局分析, 按前述方法计算的空间信息熵值结果如图 2所示。由图 2可知, 三期土地利用数据的H_s连续降低, 反映生态系统格局的有序性提升, 生态复杂度增强, 生态可持续性整体状态好转。据此, 判断出延河流域处于生长期或恢复期。

结合延河流域空间信息熵的具体计算过(表 1)程, 进一步分析影响熵值变化的主要因素。根据表 1中2000、2010、2015年的H_si值, 可知11种土地利用类型中, 耕地、中覆盖度草地、低覆盖度草地、灌木林4种类型对H_s计算结果的贡献最多。从各地类H_si值的变化可以看出, 耕地面积及其邻近度的改变对H_s的影响最大, 中覆盖度草地的影响次之。

3.2 延河流域生态动态特征分析

基于延河流域2000—2018年NDVI年极大值数据, 按前述方法计算得到NDVI时间信息熵结果(图 3)。图 3中掩模区为非植被区域, 包括居民用地、水域、其他建设用地和未利用地。由图 3可知, H_t呈现空间异质性, 总体分布特征为：研究区南部、西南部的熵值最低, 北部和中部的熵值较高, 西北、东南两端的熵值呈交错分布的特点。

由于NDVI的大小与植被类型有关, 结合延河流域的土地利用类型图, 比较植被覆盖区域中不同植被类型的H_t值情况, 结果如图 4所示。由图 4可知, 有林地、灌木林、高覆盖度草地熵值较低；耕地、中低覆盖度草地、其他林地熵值较高, 且此四类之间的彼此差异较小。

3.3 延河流域生态可持续性分析

根据3.1节, 延河流域生态可持续性整体状态好转, 处于生长期或恢复期。对应于这一宏观状态, 时间信息熵值越大, 意味着生态弹性能力越强, 生态系统越有序, 可持续性越强。因此, 根据H_t计算结果判断生态可持续性, 利用自然间断点法, 将H_t结果划分为5个等级：弱(9.04—10.21)、较弱(10.21—10.77)、中等(10.77—11.13)、较强(11.13—11.49)、强(11.49—12.41), 以此表征生态可持续性等级, 各等级的空间分布及各自所占比例如图 5所示。

由图 5可知, 延河流域生态可持续性以“强”和“较强”为主(61%), 且分布范围广(尤其是中部和北部), 表明流域的生态弹性能力总体增强, 生态可持续性较好；“中等”类型则散乱地分布在其他类型周围；“弱”、“较弱”占比最低, 集中在流域的西北和南部两端, 表明该区域的生态弹性能力尚未增强或增强幅度较少, 生态可持续性需进一步提高。

3.4 结果验证与建议

关于生态可持续性的研究通常难以进行直接的验证, 但仍可从以下两个方面作进一步地讨论：

(1) 自2000年起, 一系列生态建设工程在延河流域得到有效开展, 多数研究结果表明生态系统质量、生态系统服务功能得到提升^{[29, 32-33]}。由于流域生态系统是一有机整体, 生态系统质量和生态系统服务功能的提升, 也反映出生态系统结构和功能的改善, 即生态可持续性状况好转。与本文得出的延河流域生态系统整体处于生长期或恢复期, 朝着有序的方向发展的结论相符合。

(2) 在2013年延河流域遭遇了百年一遇的大暴雨, 此次暴雨无论降雨强度、笼罩范围都远超其他年份记录, 但并未造成以往常见的洪水与水土流失现象^[39]。从生态系统自身的属性特征出发, 这得益于生态系统的复杂度和弹性能力的增强, 有效地抵抗了外界的干扰, 这与时间信息熵及生态可持续性分析结果相吻合。

生态可持续性的研究, 可为延河流域生态保护和治理提出建议。在退耕还林还草诸项措施中, 耕地、中覆盖度草地的面积及邻近度的改变是空间信息熵值改变的主要影响因素。并且耕地、中覆盖度草地的时间信息熵值较高, 表明其生态弹性能力更强。因此, 建议今后重点采取“耕地-中覆盖度草地”这一模式。同时, 二者的分布应相对集中, 使其邻近度降低(如增大质心距离, 减少公共边界长), 有助于降低空间信息熵值。对于时间信息熵的计算结果中生态弹性无明显增强的区域, 应继续开展生态保护和修复措施。由于生态弹性与植被类型有关, 在开展生态保护和修复时, 还应兼顾不同植被类型自身的特点。基于生态可持续性分区图(图 5), 提出以下两项具体措施(图 6)：(1)对于等级为“弱”和“较弱”的区域, 应加强生态修复。根据植被类型, 修复措施又可分为人工调节和自然恢复, 前者针对耕地、经济林地, 后者包括天然林、灌木林及各类草地。(2)对于等级为“中等”的区域, 应加强生态保护。

4 结论

面向生态可持续性研究, 本文视生态系统为具有自组织能力的生命体分析其生态可持续性分析, 是一种新的尝试, 或可以提供研究生态可持续性的新视角。该视角具有交叉学科的特点, 初步将生态系统格局和动态结合起来, 暂称之为时空信息熵方法, 其中表征格局有序性的空间信息熵结果反映生态可持续的整体状态, 提供关于系统所处宏观状态的客观信息, 为动态特征分析提供参考；时间信息熵与土地利用格局结合起来, 分析生态系统动态演变的有序性, 可望更综合诠释研究区域生态学意义。

由于已有关于延河流域研究鲜有从整体的视角分析生态系统格局的变化特征, 本文进一步研究延河流域生态系统格局的总体状态及其变化特征是非常有必要的。对于延河流域在退耕还林还草工程实施过程中的植被恢复可持续性问题, 本文从生命体观出发, 通过植被覆盖连续动态变化的有序程度, 分析当地植被恢复的效果, 可为其今后的生态修复和治理工作提供参考。此外, 延河流域地处陕北黄土高原腹地, 是具有代表性的黄土丘陵沟壑区。黄土高原其他类似流域可同样被视为生命体, 通过选取符合当地流域特点的生态监测参数, 时空信息熵方法可望为其生态保护和修复提供借鉴和参考。

无论是时空信息熵概念, 其表达式及结合遥感数据的具体计算, 还是其尺度效应问题, 都需要进一步系统化的研究。毕竟, 生态可持续性问题本身是非常复杂的, 但基于生命体观的熵的研究视角, 是值得继续的有益尝试。未来研究需要从生命体观及时空尺度出发, 探讨流域生态可持续问题的时空关联性。