19942018年松嫩平原西部月亮泡水淹区生态变化及其对水淹频次的响应*

李晓东，宋开山，闫守刚

(1：滨州学院，山东省黄河三角洲生态环境重点实验室，滨州 256603)(2：中国科学院东北地理与农业生态应用研究所,中国科学院湿地生态与环境重点实验室，长春 130102)(3：枣庄学院旅游与资源环境学院，枣庄 277160)

湖泊水淹频次是指一定周期内单元水淹的次数.湖泊水域范围的时空变化特征决定了湖泊湿地的景观结构与格局.因此，湖泊水域淹没范围与淹没频次是影响湖泊湿地生态系统结构与功能的重要因子[1-3].湖泊水文过程与湖泊生态系统相互作用和影响，不同程度的水文过程对湖泊湿地生态系统的作用并不一致.而一定程度合适的湖泊淹没频率可以促进周围湿地生态系统的多样性[4].另外，不同类型的湖泊也有着不同的水文生态过程.相关研究表明，半干旱地区浅水型湖泊的时空变化相对复杂，枯水期与丰水期的湖泊面积往往相差甚远[5].湖泊和水库的库存水是人类活动、农业灌溉和社会发展所必须的根本保障[6].因此，湖泊水淹区的生态变化是湖泊生态系统安全与健康的重要保障[7-9].

近年来，在气候变化和人类活动的共同影响背景下，国内湖泊的生态状况发生了显著变化，严重影响了湖泊水文过程[10-12].同样，湖泊及其周围地表的蒸散量和地下水、地表水流量的交互作用对维持与恢复环境的生态意义更为明显[13-14].监测湖泊水文时空变化，评价湖泊水淹区域的生态变化及其生态效应具有重要意义[15].Casanova和Brock[16]通过实验证明了洪水淹没的持续时间和频率调节着湖泊及其周边的生态结构与功能.Amir等统计分析了湖泊水淹面积变化对局地气候因子的影响，发现世界第二大咸水湖——伊朗尔米亚湖水淹面积的降低减缓了当地大气流动的速度；在暖干的气候条件下，湖泊对局地气候的影响程度随着湖泊水淹面积的缩减而降低[17].湖泊水淹频次与其周边生态因子的响应关系是近期的研究热点，目前研究把主要关注方向均放在了时序变化上[18-21].

随着遥感技术的发展，多类型传感器提供的多源遥感数据可以实现不同尺度数据的融合与互用，定量衡量不同尺度下地表生态评价及其变化成为可能[22-24].不同水文模型的构建也是湖泊生态变化检测在全球遥感技术实现过程中的一个有益探索，物理模型[25]、随机模型[26]和系统动态模型[24]等更多数学模型被构建用于湖泊生态变化检测.徐涵秋[27]提出了用于衡量地表生态状况的生态指数(remote sensing base ecological index，RSEI).RSEI综合指数通过主成分分析方法集成了绿度、湿度、干度和温度等生态指标来评价整个研究区域的生态环境水平.地表生态变化检测研究主要是基于单一时相影像对某区域地表生态状况进行评估.如何利用年内多时相数据和综合多维生态指数共同表征地表的状态变化，从而提高地表(包括湖泊生态)变化信息的提取精度，是另一个主要关注的问题.

本文以月亮泡湖泊湿地为研究对象，以其在水淹区域内地表生态的时空变化特征为重点研究内容.1994-2018年月亮泡的水淹区域的研究涉及以下内容:(1)完成研究区湖泊湿地遥感分类方案的构建，计算提取月亮泡湖泊湿地生态变化的(1995-2006、2006-2016年)年际变化范围及其方向和实际年际变化类型；(2)完成月亮泡湖泊湿地生态变化的(1995、2006、2016年)年内波动等级和主要生态因子的生长季累积量的计算提取；(3)定量衡量湖泊湿地生产能力等相关生态参数，进一步确定湖泊年际水淹频次和湖泊湿地生态变化的时空相关关系.

1 研究区概况与数据准备

1.1 研究区概况

松嫩平原腹地的月亮泡湖泊湿地(45°35′～45°50′N，123°37′～124°9′E)位于吉林省西部洮儿河下游，嫩江西侧，界于吉林省大安、镇赉等县市之间.月亮泡湖泊湿地主要包括月亮泡和新荒泡.月亮泡湖底平坦，为游泥及细沙底质.平水期水体面积达181 km2，最大面积达205.7 km2.常年平均水深3～4 m，最深5 m.结冰期为每年10月-次年3月，最大冰厚1.13 m，多年年平均冰厚0.89 m(图1).

图1 研究区位置图

月亮泡为天然淡水湖泊，对洮儿河、嫩江有拦蓄调洪作用.1976年修筑围堤、排水节制闸等湖区水利配套工程，成为吉林西部最大调蓄水利枢纽.2000年，月亮泡湿地作为东北地区重要湿地被列入中国重要湿地名录[28].近年来，吉林西部的“河湖连通”工程实现了月亮泡湖区的水资源调配，进一步改变了湖泊生态系统的分布格局[29-30].同时，1994-2018年极端气象事件对吉林西部月亮泡湖泊湿地的生态功能和结构产生了巨大影响，特别是1998和2013年2次大规模洪水交替发生.因此，月亮泡湖泊湿地既有平原内部湖泊的变化共性，又有自己的独特性.

1.2 数据预处理

研究用数据分别为中国资源卫星应用中心(http://www.cresda.com/CN/index.shtml)提供的环境1A和1B的CCD数字产品和美国地质调查局(USGS)地球资源观测与科学中心(EROS)提供的Landsat TM和 OLI图像(http://earthexplorer.usgs.gov/).其中，1994-2009年卫星影像数据为Landsat TM5卫星数字产品；2014-2018年数据为Landsat 8 OLI_TIRS 卫星数字产品；2013年数据为中国的环境卫星HJ-1A和1B的CCD数字产品.在HJ-1A卫星和HJ-1B卫星上装载的两台CCD相机为4波段的推扫式成像仪，星下点的空间分辨率为30 m，用于环境与灾害的监测预报.为避免单一时相影像对年际分类结果的影响，利用5年内的清晰无云TM和OLI影像(<10%的云量)组建基准年(重点观测年份).根据月亮泡年际水淹没区域的实际卫星数据收集情况，1995、2006和2016年集中了研究区地表生长盛期的关键月份(4-10月)，被确定为月亮泡水淹区生态变化检测的基准年.

卫星影像数据的处理涉及辐射定标、大气校正和几何精校正.辐射定标将原始影像的像素值(digital number，DN)转换为传感器处反射率(top of atmosphere reflectance, TOA)，大气校正算法为FLAASH；为了减少不同年份的影像反射差异，相对辐射校正又称辐射归一化，被用于不同卫星传感器数据的配准.不同时相影像之间的配准采用二次多项式和最邻近象元法，配准的均方根误差小于0.5个象元.不同于数据预处理过程中的图像辐射校正，处理的目的是减小两期卫星影像数据中辐射波动的影响[31].遥感影像的相对辐射校正主要应用于时序影像数据集合的构建，消除因大气和光照条件的不同而引起的两期影像数据的辐射差异[32].

1.3 地表生态特征指数及湖泊年际水淹频次计算

为了提高地面物体的变化检测精度，修正的归一化水体指数(MNDWI)结合归一化植被指数(NDVI)、归一化土壤指数(NDSI)被用于计算提取研究区域内土地覆被类型的变化范围，并进一步确定其变化类型.

1.3.1NDVI的计算NDVI计算公式为[33]：

NDVI=(NIR-Red)/(NIR+Red)

(1)

式中,Red代表红光波段，而NIR为近红外波段.

1.3.2MNDWI及湖泊年际水淹频次的计算MNDWI计算公式为[34]：

MNDWI=(Green-SWIR1)/(Green+SWIR1)

(2)

式中,Green表示绿光波段，SWIR1表示Landsat卫星传感器的短波红外光谱波段.月亮泡年际水淹频次的获取过程：① 计算修正的水体归一化指数(MNDWI)；②确定水体划分阈值，累计二值化水体指数，并作为湖泊年际水淹频次；③根据参与计算的24期(1994-2018)年际数据，恒定水体的年际水淹频次为24，季节性水体的年际水淹频次为1～23.

1.3.3NDSI的计算NDSI也被称为归一化差值裸土指数，被用于裸露土地信息的提取[35].

NDSI=(SWIR1-NIR)/(SWIR1+NIR)

(3)

2 方法与验证

2.1 变化检测方案

变化检测方法的实现过程如图2.从生态学角度，本次研究与其他方法的主要区别在于(1)基于MNDWI、NDVI和NDSI等主要生态指标，计算生态指数的年内最大差值和年际动态变化率，用于定量衡量研究区的变化程度；(2)综合图像变换与像素点值运算，提取年内最大差值和年际动态变化率的第一主成分，标准化后的第一主成分分别用于定量衡量年内和年际地表类型的变化范围.

图2 生态变化检测研究方案(MNDWIs、NDVIs、NDSIs为生态指数的生长季(4-10月)月值数据集合)

2.1.1 年内变化检测 (1)最大差异值

最大差异值(max-difference)是遥感监测区域年内变化的最大波动幅度.年内变化检测研究实际收集的时序数据为j=4、5、6、7、8、9、10月(东北地区生长季的范围为4-10月)遥感影像数据，所以松嫩平原地区年内变化实际为年内生长季的月动态变化检测.生态指数(eco-index，EI)包含了研究区地表水体指数、植被指数和土壤指数等[36].各月生态指数进行最大差值计算的公式如下：

Max-Difference=Max(EIi,4,EIi,5,EIi,6,…,EIi,9,EIi,10)-Min(EIi,4,EIi,5,EIi,6,…,EIi,9,EIi,10)

(4)

式中，Max-Difference为研究区地表年内生态变化的最大差异值；EIi,j为第i个地表生态指数(MNDWI、NDVI和NDSI)第j月的值.研究区的土壤和水体信息的最大差异计算同最大NDVI差异值计算.

(2)显著性分级

综合生态指数(RSEI)是由植被、水体和土壤等指标通过主成分分析得到的一个集成数值，用于评价监测区域的生态质量.通常，RSEI被转换成85分制或百分制，进一步等间隔分割成5个值域梯度，分别对应差、较差、中等、良和优5个生态等级[37].

2.1.2 年际变化检测 (1)动态变化率

根据遥感数据的实际识别能力和陆表的生态状况，动态变化率实际上是两个样本中影像数值增加或减少的定量度量.最终计算公式如下：

DynamicRatio=(EI2-EI1)/Mean(EI1,EI2)

(5)

式中，EI1和EI2分别表示第1、2年地表生态指数(MNDWI、NDVI和NDSI)的月极大值.两幅图像的平均值(Mean(EI1,EI2))来表示同一位置的生态指数(MNDWI、NDVI和NDSI)的恒定状态，并使用两期数据集的差分来表示生态环境的变化程度.最后，将陆表的动态比率法描述为算法公式(5).

(2)显著性分级

通常，陆地覆被类型的未变化区域具有接近于零的像素值，但变化区域的像素值通常大于(或小于)其他区域[38].当差分结果大于(或小于)给定阈值T(-T)时，该情况发生区域被确定为变化区域[39].常规计算方法只关注差分结果，由此引出的问题是差值较大的点未必发生了变化，差值较小的点未必没有变化.动态比率算法(DRM)反映的是前后期陆表生态指数差异与其生态均值的对比结果.综合水体、植被和土壤等信息的生态指数，动态变化率的正、负值分别代表综合生态指数的正变化方向(改善)与负变化方向(衰退)；年际动态比率第一主成分的显著性分级方法与年内最大差值的显著性分级方法相同.

2.2 陆表类型的遥感提取方案

基于重建的NDVI、NDSI与MNDWI数据序列集合，分别提取研究区的植被信息、土地利用信息和水体信息，进一步计算并获取研究区的植被、土地利用与水体的变化信息(包括全年极大值、变异系数、生长季累积量的全年占比)，最终计算提取的变化信息构建适合于研究区湿地分类提取的多维特征数据集[40].分类方案如下：NDVI极大值+NDVI变异系数+NDVI生长季累积量的全年占比+NDSI极大值+NDSI变异系数+MNDWI极大值+MNDWI变异系数+面向对象+SVM分类算法，实现平台分别为：ERDAS 9、eCognition 9等遥感软件.ERDAS软件完成相关比值的运算提取，同时，eCognition软件通过面向对象+SVM分类算法实现研究区湿地类型的遥感解译提取，得到湿地分类结果.

2.3 验证点与遥感提取精度

为了确保实验区内地表覆被类型的代表性，已有的遥感解释分类结果——全球表面覆盖率遥感数据产品(GlobeLand30)和 Google Earth 高分辨率图像作为采样点上地表类型数据集的参考依据.同时，109个2016年9月的野外采样点和217个基于Google Earth高清卫星影像的人工解译样点和2006、2010年东北地区土地利用现状统计资料，在遥感影像图上选出各类具有代表性的样区.把选取的样本随机分为2部分：70%用于分类，30%用于精度评价[41-42].在1994-2018共24年的分类结果中，7个年份的遥感解译精度在90%以上，比较真实反映了实际土地覆被类型分布(表1)，用于研究区湖泊湿地主要覆被类型的变化检测.

表1 研究区1994-2016年的土地覆被类型分类精度

2.4 湖泊年际水淹频次与生态因子的相关性研究

为确定月亮泡的年际水淹频次是否影响湖滨湿地生态状态，研究区采样方案如下：首先，利用叠加函数计算湖泊水体的多年年内波动频次(图3a)与生长季MNDWI、NDVI和NDSI累积量；其次，利用eCognition 软件平台的面向对象分割算法，将研究区内的非湿地类型和恒定水体从矢量数据(SHP)中剔除.最后，通过区域性统计，得到四级缓冲距离的提取结果(图3b)，以验证湿地生态状态与湖泊年际水淹频次的时空相关性.

图3 月亮泡湖泊年际水淹范围的多年水淹频次(a)及四级缓冲区(b)

3 结果分析

3.1 月亮泡的主要湿地类型和年际水淹范围

3.1.1 主要湿地类型的长时间序列变化趋势 通过研究区主要湿地类型的面积变化来分析自然过程和人类活动强度对月亮泡湖泊湿地的影响.研究区湿地类型分为水田、恒定水体、沼泽湿地3类(图4).利用线性回归分析方程对3种湿地类型的区域变化特征进行了统计分析.统计结果表明，水田面积持续性增加最明显，趋势值为11.746 km2/a(函数方程y=11.746x-2.909,R2=0.817).水田是该区域人类活动的主要形式，成为湖泊湿地的最为持续的干扰因子.研究区水体分析结果未通过显著性检验(P>0.05)，恒定水体主要变化特征为周期性增减，总体面积稳定.此外，研究区沼泽的函数方程为y=-20.663x+314.9(R2=0.736)，表示沼泽的区域变化呈持续性下降趋势，趋势值为-20.663 km2/a.

图4 主要湿地类型的变化特征

3.1.2 湖泊年际水淹没区域 1994-2018年，湖泊的年际水淹没区域主要分布在月亮泡北侧和西侧尾闾，年际水淹频次的均值分别为18(波动范围12～24)、15(波动范围13～23).湖泊北侧年际水淹没区域的面积为58.36 km2，月亮泡西侧尾闾的年际水淹没面积为25.07 km2.综上，月亮泡年际水淹没区域总面积为83.43 km2，占整个湖泊影响范围(恒定水体面积与季节性水体面积)的56.85%.月亮泡是一个典型的内陆潮汐性湖泊，湖泊北侧的年际水淹没区域面积更大，年际水淹频次更高(图5).月亮泡湖泊湿地的范围划分依据，沿嫩江河堤和省、地级公路为主要边界圈定本次研究范围(图5).A与B区域的分割来自区域内主要人工建筑物(道路和桥梁)和实际年际水淹频次.

图5 1994-2018年月亮泡的水淹没范围及其水淹频次(A、B分别为月亮泡北侧和西侧尾闾年际水淹没区域，下同)

3.2 湖泊水淹区内湿地生态的年际变化与类型

3.2.1 湖泊湿地生态的年际变化 (1)变化等级与范围

为了便于分析地表生态变化方向和显著性等级，动态比率的第一主成分被标准化成[-1，1]的值域范围，进一步将其分隔成5个生态等级：[-1，-0.6)、[-0.6，-0.2)、[-0.2，0.2]、(0.2，0.6]、(0.6，1]，分别代表显著负变化方向、偏负方向、稳定、偏正方向、显著正变化方向.最终，[-1，-0.6)和(0.6，1]作为两个极显著变化值域，被确定为年际动态检测的变化范围.

1995、2006和2016年作为主要年际变化检测时段，研究区地表主要生态指数(MNDWI、NDVI、NDSI)年际变化的动态比率计算第一主成分，用于定量衡量该区域内的年际变化方向.研究区地表生态的年际变化方向的结果表明(图6)：1995-2006年，月亮泡北侧年际水淹区为该检测时段湿地生态的主要正向变化(改善)区域.2006-2016年，月亮泡西侧尾闾的年际水淹区为该检测时段湿地生态的主要正向变化(改善)区域，同时月亮泡北侧年际水淹区的地表生态变化相对稳定.

图6 1995-2006(a)和2006-2016年(b)月亮泡湿地生态的年际变化

(2)面积占比

1995-2006年，湖泊北侧水淹没区域内湿地生态的年际变化以改善为主要特点，面积占比为39.96%(A区域面积)；湖泊西侧尾闾湿地生态改善区域占B区域面积的19.56%.2006-2016年，月亮泡水淹没区域湿地生态的年际变化相对稳定，北侧湖泊湿地生态的年际变化改善区域面积占比为1.51%(A区域面积)；西侧尾闾湿地生态的年际改善面积占B区域面积的38%.综上，月亮泡西侧尾闾的地表生态年际变化以正向变化(改善)为主要特点，同时湖泊北侧的地表生态年际变化趋于稳定(表2).

表2 月亮泡湖泊湿地生态的年际变化面积占比

3.2.2 主要年际变化类型 月亮泡西侧尾闾沼泽湿地在水淹没区域B内的面积占比从1995年的80.29%下降到2016年的23.69%.同时，恒定水体面积占比从1995年的15.1%增加到2016年的74.87%；北侧沼泽湿地在水淹没区域A内的面积占比从1995年的92.94%减少到2016年的25.39%.同时，恒定水体面积从1995年的15.1%增加到2016年的74.87%(图7).综上，在主要观测时段内，恒定水体面积的增加与沼泽湿地面积的减少是月亮泡水淹区的主要变化.

图7 月亮泡北侧(a)和西侧尾闾(b)水淹区内土地植被覆盖类型的面积占比

3.3 湖泊水淹区内湿地生态的年内变化

3.3.1 湖泊湿地生态的年内波动级别 为了方便进行显著性分析，最终该结果的第一主成分被正规化到[0，1]区间，年内的最大变化差异划分为：[0，0.2)、[0.2，0.4)、[0.4，0.6]、(0.6，0.8]、(0.8，1]，分别代表极不显著、不显著、中等、显著、极显著5个生态等级.

从1995、2006、2016标准年的生态波动级别分析表明，2006年，月亮泡北侧湿地的年内生态波动的极显著等级(L-5)面积占比为88.44%；2016年，月亮泡西侧尾闾湿地的一般变化等级(L-3)为年内生态波动等级的最大面积占比(56.8%)(图8).综上，月亮泡北侧湿地的年内生态波动级别比湖泊尾闾更为显著.

图8 月亮泡湿地北侧(a)和西侧尾闾(b)年内生态波动等级的面积占比

3.3.2 湖泊湿地主要生态因子的生长季累积量 研究区地表生态因子的累积量变化表明，2016年月亮泡尾闾水体指数的年内累积量增加最为显著，为3.487±1.503.同时，湖泊西侧尾闾地表植被指数的年内累积量的衰减也较北侧显著.在主要观测时段内，月亮泡湖泊尾闾具有较低的年际水淹频次，年内生长季水体信息(增加)与植被累积量(减少)的时序变化更为显著.综上，月亮泡湖泊尾闾湿地水进土退、水进植被退的时序变化特点更为显著(表3).

表3 月亮泡水淹区内主要生态因子的生长季累积量

3.4 生态因子与湖泊年际水淹频次的空间相关性

选择Spearman秩相关分析来定量衡量NDVI生长季累积量与月亮泡水淹频次的多年平均值的空间相关性.结果表明，在距离月亮泡湖泊4.5 km的缓冲区上，地表具有相对稳定的湖滨湿地生态环境，湖泊年际水淹频次最低，平均值为0.441；地表植被累积量最高(1.587±0.529)(表4、5).在1.5～3.0 km范围内，湖泊年际水淹频次与生长季MNDWI、NDVI和NDSI累积量的空间相关性相对显著，相关系数分别为0.719、-0.683和-0.561；湿地植被累积量较高(1.581±0.421).在缓冲区0～1.5 km范围内，湖泊年际水淹频次与生长季MNDWI、NDVI和NDSI累积量的空间相关性最高(0.878、-0.890和-0.756)，通过了0.01显著双侧检验，但是湿地植被累积量最低(1.148±0.519).湖泊年际水淹频次影响着湖泊湿地MNDWI、NDVI和NDSI累积量，最终导致湖滨湿地生态变化.综上，月亮泡湖泊年际水淹频次在一定时间和空间上影响着水淹区内的植被生产力(负相关关系)，通过年内生态波动，进一步影响地表生态年际变化.

表4 月亮泡地表 NDVI与多年平均水淹没次数

表5 月亮泡水淹区生态因子与多年平均水淹频次的空间相关性*

3.5 气候因子与湖泊年际水淹没面积的时序相关性

月亮泡湖泊湿地变化既有平原内部湖泊的变化共性，又有自己的独特性.月亮泡湖泊湿地位于东北平原腹地，地形结构相对封闭，接受外界供水较少.洮儿河和霍林河只有在汛期可以维持一定的水量，很难成为河流中下游沿岸湖泊湿地的有效补给资源.同时，统计分析结果表明，以4年作为降水的累积周期，月亮泡湖泊面积与降水年累积量的相关系数最高，R2为0.973(图9).研究区4年的降水累积量直接影响月亮泡恒定水体面积，进而影响湖泊年际水淹频次，是湖泊恒定水体面积和年际水淹频次的主要影响因素.

图9 气候因子(降水/气温)与恒定水体的相关性分析

4 讨论

基于陆地资源卫星与环境卫星作为数据源，利用动态比率与最大差值算法分别定量衡量湿地生态的1995、2006和2016标准年年际/年内变化，进一步分析吉林西部月亮泡湖泊湿地1994-2018年年际水淹频次及其对地表生态环境的影响.近年来，月亮泡年际水淹频次下降、湖面扩展、沼泽湿地逐年萎缩，进而湖泊湿地生态因子中植被生产能力在短时期内上升.

4.1 月亮泡湖泊湿地水淹区的气候变化

吉林省是内陆湖泊湿地资源的重要分布省份.湖泊湿地主要集中在吉林中西部的干旱、半干旱地区，湖泊湿地的空间变化具有较强的干旱、半干旱地域性特征[43-44].同时，实际研究工作也表明，吉林省西部湿地总面积一直呈稳步减少趋势.Mao等指出，在1990-2010年中国最常见的自然湿地中，沼泽湿地是呈现下降趋势的主要湿地类型.但是，湖泊湿地(恒定水体)经历了3.9%的面积扩展[45].2012-2017年，近地表湿度的函数趋势方程为y=-1.5781x+72.017(R2=0.8294)，近地表大气湿度降低趋势显著(图10).生长季气候因子的变化不同于气候因子的年际变化，因此生长季的主要气候因子与内陆湖泊水淹范围及其水淹频次的相关关系值得进一步研究.

图10 2003-2017年研究区生长季(5-10月)空气湿度和降水/蒸发量差异

4.2 月亮泡湖泊水淹区的土地覆被类型转换

统计分析月亮泡湖滨湿地生态变化(改善/衰退)区域内的主要土地覆被转换类型(表6).首先，湖泊恒定水体面积扩展是1995-2006和2006-2016年月亮泡北侧湿地主要改善类型，分别为31.093、28.429 km2.其次，2006-2016年，月亮泡西侧尾闾湿地表现为大规模沼泽化，各种地物类型共计转出成为沼泽湿地的面积为10.4 km2.月亮泡湖泊湿地类型与其他地表植被覆被类型间的转换是本区域内生态变化的主要形式.研究区地表生态衰退的主要形式是湖泊湿地开垦为耕地.

表6 月亮泡水淹区生态变化与主要土地覆被类型转换(km2)

近年，大量的农业用水在一定程度上减少了湖泊水淹区域，降低了湖泊水淹频次，最终促进了当地湖泊的沼泽化.周静等[46]的研究表明，湖泊水位下降对湖泊湿地中芦苇分布面积的扩张起促进作用.湖泊水淹区域减少，最终水淹频次下降，地表植被在短时间内得到一定程度的生物量累积.另外，月亮泡是松嫩平原内部的半天然内陆湖泊湿地，具有占比最大的湖泊水淹区域.月亮泡是莫莫格湿地的重要组成部分，湖泊淹水与盐度交互作用对水生生物群落的株高及其生物量有着显著影响[47].刘舒[48]发现2000-2016年吉林西部湿地植被净初级生产力指数(NPP)不断增加，湿地草本植物夏季累积量的上升趋势明显.

4.3 月亮泡水淹区的土地盐碱化

松嫩平原是世界上主要的苏打盐碱地分布区之一.土壤中微量的钠盐会显著增加土壤的酸碱度，从而导致土壤盐渍化[49].研究区土地盐渍化会直接改变水从土壤深层向地表运移的过程，表层积水逐渐转化为碱性土地，碱性土地具有相对稳定的面积和无生物量产出等特点[50-51].1994-2006年，月亮泡湖泊湿地由于土地盐碱化导致的生态退化主要发生在湖泊西侧尾闾湿地，为8.94 km2；2006-2016年，月亮泡湿地由于土地盐碱化导致的生态退化主要发生在月亮泡湖泊北侧湿地，为2.33 km2(表7).盐碱化土地的扩大是月亮泡湖泊湿地土壤裸露程度增加的主要形式，也是研究区湖滨湿地生态稳定的主要威胁[52].

表7 月亮泡水淹区生态变化与土地盐碱化(km2)

5 结论

基于遥感卫星数据的月值序列集合，综合水体、植被和土壤指数(MNDWI、NDVI、NDSI)对研究区地表生态变化进行评价，进一步检测和分析湿地生态因子对湖泊年际水淹频次的响应，最终综合时空变化分析结果表明：

1)月亮泡是一个典型的内陆潮汐性湖泊.1994-2018年湖泊北侧和西侧尾闾存在2个湖泊水淹区，湖泊北侧年际水淹频次明显大于湖泊西侧尾闾的水淹频次.从年际变化类型和生态因子累积量的统计结果表明，月亮泡年际水淹区内具有明显的水进土退、水进植被退化的生态变化特点，其中湖泊北侧的变化更为显著.

2)月亮泡湖泊湿地生态变化存在显著的空间区域性，西侧尾闾的地表生态年际变化以改善主要特点，湖泊北侧的地表生态年际变化趋于稳定.同时，湖泊湿地生态变化存在显著的季节性，湖泊尾闾湿地的年内生态波动属于中等级别，湖泊北侧湿地的年内生态波动更为显著.

3)研究区4年降水量与恒定水体呈显著正相关关系，即降水直接影响恒定水体面积的变化，间接影响湿地生态功能.同时，湿地生态对恒定水体的响应范围约为1.5～3.0 km.综上，研究区湿地的变化机制表现为气候-水-湿地的相关性.因此，合理科学调配嫩江的入湖水量和维持研究区自然湿地水位是月亮泡湖泊湿地恢复的必要条件；月亮泡北侧与西侧自然沼泽湿地的恢复是该区域生态保护的关键.