辽河口盐地碱蓬时空动态遥感监测及其识别机理研究

doi:10.6046/zrzyyg.2023293

辽河口盐地碱蓬时空动态遥感监测及其识别机理研究

李钰彬^,¹^,², 王宗明², 赵传朋^,², 贾明明², 任春颖², 毛德华², 于皓¹

1.吉林建筑大学测绘与勘查工程学院,长春 130118

2.中国科学院东北地理与农业生态研究所,长春 130102

Remote sensing-based monitoring and identification mechanisms of the spatiotemporal dynamics of Suaeda salsa in the Liaohe estuary, China

LI Yubin^,¹^,², WANG Zongming², ZHAO Chuanpeng^,², JIA Mingming², REN Chunying², MAO Dehua², YU Hao¹

1. School of Surveying, Mapping and Exploration Engineering, Jilin Jianzhu University, Changchun 130118, China

2. Key Laboratory of Wetland Ecology and Environment, Northeast Institute of Geography and Agroecology, Chinese Academy of Sciences, Changchun 130102, China

通讯作者: 赵传朋(1991-),男,副研究员,研究方向为滨海湿地遥感。Email:zhaochuanpeng@iga.ac.cn。

责任编辑: 陈昊旻

收稿日期: 2023-09-22 修回日期: 2024-03-11

基金资助:

国家重点研发计划青年科学家项目“地上植被生物量广域精细多模观测技术与装备”(2022YFF1302000)
国家自然科学基金青年基金项目“基于样本扩增方法与多源卫星影像的无瓣海桑扩散进程监测方法研究”(42201422)
“基于Sentinel光学和雷达遥感影像的泥炭沼泽识别方法研究”(42101399)
第五批吉林省青年科技人才托举工程“泥炭沼泽的信息提取研究”(QT202101)

Received: 2023-09-22 Revised: 2024-03-11

作者简介 About authors

李钰彬(1999-),男,硕士研究生,研究方向为滨海湿地遥感。Email: 18636183172@163.com。

摘要

作为全球面积最大的红海滩景观,监测辽河口盐地碱蓬的时空动态变化对揭示其“退养还湿”等保护措施成效具有重要意义。目前,卫星遥感技术已广泛应用于包括盐地碱蓬在内的滨海植被识别与制图,但现有分类方法依赖于难以解释的黑箱模型,忽视了对识别机理的探究,制约了相关方法的改进和发展。可解释人工智能的发展为黑箱算法的解析指出了新的方向。考虑到构成随机森林的决策规则具有可解释性,本研究构建了一套从已训练随机森林模型中抽取最优决策规则的新方法,最终重构得到识别盐地碱蓬的最优决策规则,即B3/B4<0.90且B5/B3≥1.46,数据整体精度优于90%; 以2017—2022年的Sentinel-2影像为数据源,实现了对辽河口盐地碱蓬的逐年动态提取,并结合质心迁移法,分析了“退养还湿”工程实施后盐地碱蓬时空变化,揭示了该区域盐地碱蓬呈现快速恢复的现状。

关键词： 退养还湿; 碱蓬; 可解释人工智能; 随机森林; 决策规则

Abstract

The Liaohe estuary of China boasts the largest red beach landscape in the world. Monitoring the spatiotemporal dynamics of Suaeda salsa in this region is of great significance for revealing the performance of conservation measures such as returning aquaculture to wetlands. Currently, satellite remote sensing technology has been widely applied to the mapping and identification of coastal vegetation including Suaeda salsa. However, existing classification methods rely on black-box models, which are difficult to interpret, while overlooking exploring identification mechanisms. This has hindered the improvement and development of related methods. Fortunately, the advancement in explainable artificial intelligence (XAI) has provided new directions for analyzing the black-box models. Considering that the decision rules in random forests are interpretable, this study developed a new method to extract the optimal decision rules from trained random forest models. Using this method, this study ultimately reconstructed the optimal decision rules used to identify Suaeda salsa, i.e., B3/B4<0.90 & B5/B3≥1.46, with an overall data accuracy exceeding 90%. Using annual Sentinel-2 images from 2017 to 2022 as a data source, the study successfully extracted the annual dynamics of Suaeda salsa in the Liaohe Estuary. Accordingly, by combining the centroid migration method, this study analyzed the spatiotemporal changes in the Suaeda salsa following the implementation of returning aquaculture to wetlands, revealing the current status that the Suaeda salsa in this region is undergoing rapid restoration.

Keywords： returning aquaculture to wetlands; Suaeda salsa; explainable artificial intelligence (XAI); random forests; decision rule

PDF (5252KB) 元数据多维度评价相关文章导出 EndNote| Ris| Bibtex 收藏本文

本文引用格式

李钰彬, 王宗明, 赵传朋, 贾明明, 任春颖, 毛德华, 于皓. 辽河口盐地碱蓬时空动态遥感监测及其识别机理研究[J]. 自然资源遥感, 2025, 37(1): 195-203 doi:10.6046/zrzyyg.2023293

LI Yubin, WANG Zongming, ZHAO Chuanpeng, JIA Mingming, REN Chunying, MAO Dehua, YU Hao. Remote sensing-based monitoring and identification mechanisms of the spatiotemporal dynamics of Suaeda salsa in the Liaohe estuary, China[J]. Remote Sensing for Land & Resources, 2025, 37(1): 195-203 doi:10.6046/zrzyyg.2023293

0 引言

辽河口国家级自然保护区是我国沿海最大的河口湿地,其中的碱蓬形成了全球面积最大的“红海滩”景观^[1-2],具有稳定海岸线、净化污染等生态价值和旅游观光等经济价值^[3⇓⇓-6]。近30 a来,受海岸开发等人为因素的影响,该区域盐地碱蓬群落退化,覆盖面积大幅减少,影响了辽河口湿地生态系统的稳定性。2018年起,该区域全面启动“退养还湿”工程。为了监测这一工程的实施成效,需要对碱蓬的时空动态变化进行高效监测。

卫星遥感技术具有成本低、时效性强和观测范围广等优势,已逐渐应用于碱蓬群落的动态监测研究^[7-8]。依据时间序列卫星数据的使用方式,现有研究主要分为2类。其一是获取像元尺度的土地覆被类型变化轨迹,例如,彭剑伟^[9] 基于CCDC算法对1986—2018年共803景Landsat卫星影像进行分析研究,得到辽河口地区碱蓬群落的逐年分布数据,但该方法易受不同时相、不同卫星姿态及不同传感器等外界因素影响。其二为直接获取不同时段土地覆被分类结果,该方法通常结合人工智能算法以提高分类结果的准确性。Song等^[10] 基于Sentinel-2卫星影像与随机森林算法获取了辽河口湿地景观演化,但用户精度偏低; 王文硕^[11] 通过归一化植被指数(normalized difference vegetation index,NDVI)时间序列表征辽河口典型植被物候信息,并基于支持向量机(support vector machine,SVM)分类方法得到包含碱蓬在内的各土地覆被分布信息,但升维导致模型难以被解释。

深度学习模型Alex-Net能改善高光谱遥感影像在盐地碱蓬等湿地植被分类过程中出现的错分和漏分现象^[12],但海量参数也导致模型难以被高效解释。基于人工智能算法构建的黑箱模型出现决策偏见等问题,推动了可解释人工智能的发展^[13]。可解释性成为解决当前黑箱模型信任危机的关键,即人们可直接理解抽象复杂模型内部决策过程的能力^[14-15]。李营等^[16] 基于不同地物类型的反射率信息,提出了盐地碱蓬植被指数(suaeda salsa vegetation index,SSVI),在通过多次实验阈值调节后,实现了对局部地区翅碱蓬分布信息的获取,但这类方法往往依赖于大量实验,构建和改进难度较大。在当前大数据时代,通过挖掘已训练机器学习模型,从数据中发现知识,或可解决上述问题。根据可解释人工智能,决策规则是一种可被人们理解的知识表达形式。目前广泛使用的随机森林算法内部集成有大量决策规则,其具备训练效率高、稳健性强、准确度高等优势^[17-18]。因此,如何根据已训练随机森林模型内部的众多决策规则重构特征数量少、准确度高的最优决策规则^[19],是解释盐地碱蓬遥感识别机理的核心问题。

本文以辽河口盐地碱蓬为研究对象,Sentinel-2影像为数据源,构建了一套从已训练随机森林模型中抽取最优决策规则以提取盐地碱蓬的可解释分类方法。在此基础上,获取2017—2022年辽河口逐年盐地碱蓬分布数据,揭示了“退养还湿”工程实施后盐地碱蓬的时空动态,为保护区的精细化管理、可持续利用和海岸生态文明建设提供了数据支撑和科学依据。

1 研究区概况及其数据源

1.1 研究区概况

辽河口国家级自然保护区(121°30'~122°00'E, 40°50'~41°20'N),位于渤海北部(图1),地势平坦低洼。该区域属温带季风气候,春季干旱少雨,降雨集中在6—8月^[10]。受河水径流量、泥沙淤积等因素影响,该地区形成了植被类型较为单一的河口三角洲湿地生态系统,其作为珍稀水禽栖息地的同时,也是重要的全球候鸟迁移通道^[20]。研究区河网纵横,芦苇碱蓬等植被分布广泛,其中碱蓬呈现条带状、块状聚集生长,对潮汐变化、人类活动干扰等高度敏感。

图1

新窗口打开| 下载原图ZIP| 生成PPT

图1 研究区地理位置及样点分布

Fig.1 Location of the study area and layout of the sampling points

1.2 数据来源

Sentinel-2遥感影像广泛用于土地覆被制图,如盐沼^[21]、滩涂^[22]等。Sentinel-2影像来自2颗极地轨道卫星(Sentinel-2A与Sentinel-2B),具有优于5 d的重访周期,10 m,20 m与60 m的空间分辨率,以及覆盖可见光、近红外至短波红外的13个光谱波段。本研究使用了Level-1C与Level-2A反射率产品,应用三次卷积重采样方法将空间分辨率转换为10 m。

结合目视解译和实地调查方式,获取样本共580个,其包含的碱蓬与非碱蓬样本数量相等。训练样本从样本库中随机选择,个数为180,剩余样本作为测试集以检验模型性能。其中,非碱蓬样本分为2组: 与碱蓬相似的地物,即滩涂、盐沼、芦苇和草本; 与碱蓬相似度较低的地物,即不透水面、养殖池和水体。

2 研究方法

2.1 分类特征集构建

辽河口盐地碱蓬于每年4—5月长出地面,其绿色特征与湿地植被如芦苇、蒲草等相似,难以利用遥感影像进行区分; 至6月,盐地碱蓬进入生长期,颜色变为浅红; 8—9月,盐地碱蓬进入花期,呈现深红色,分类特征明显; 11月后,盐地碱蓬进入枯萎阶段^[1,10]。在综合考虑盐地碱蓬不同生长阶段的影像特征后,选择9月作为盐地碱蓬的最佳分类时段,并基于Google Earth Engine (GEE)平台筛选该时段内受云及其阴影影响较小的相对低潮影像,得到2017—2022年覆盖研究区的高质量遥感影像共6景(表1)。根据已选定的影像构建了包含12个影像波段、27个波段比值特征以及13个指数特征共计52个特征的特征集,具体计算公式见表2。

表1 选中的2017—2022年卫星影像

Tab.1 Selected satellite images during 2017 and 2022

时间	产品类型	云量/%
2017-09-29	Level-1C	0
2018-09-04	Level-1C	1.08
2019-09-14	Level-2A	0.51
2020-09-18	Level-2A	0.01
2021-09-23	Level-2A	0
2022-09-18	Level-2A	0

新窗口打开| 下载CSV

表2 Sentinel-2影像特征介绍

Tab.2 Features derived from the Sentinel-2 Images

特征类型	特征
光谱特征	B2,B3,B4,B5,B6,B7,B8,B8A,B11,B12,VV,VH
波段比值特征组	B2/B4,B3/B4,B3/B8,B4/B2,B4/B3,B4/B5,B4/B8,B5/B3,B5/B4,B6/B3,B6/B5,B7/B3,B7/B4,B8/B2,B8/B3,B8/B4,B8/B5,B8/B11,B8/B12,B8A/B5,B11/B8,B11/B12,B12/B4,B12/B8,B12/B11,VV/VH,VH/VV
指数特征组^①	归一化植被指数(normalized difference vegetation index, NDVI)	$N D V I = \frac{B 8 - B 4}{B 8 + B 4}$
	增强植被指数(enhanced vegetation index, EVI)	$E V I = \frac{B 8 - B 4}{B 8 + 6 B 4 - 7.5 B 2 + 1}$
	植被衰减指数(plant senescence reflectance index, PSRI)	$P S R I = \frac{B 4 - B 2}{B 6}$
	归一化差异红色变异指数(normalized difference red edge, NDRE)	$N D R E = \frac{B 8 - R E}{B 8 + R E}$
	地表水分指数(land surface water index, LSWI)	$L S W I = \frac{B 8 - S W I R}{B 8 + S W I R}$
	归一化水分指数(normalized difference water index, NDWI)	$N D W I = \frac{B 3 - B 8}{B 3 + B 8}$
	修正的归一化差异水体指数(modified normalized difference water index, mNDWI)	$m N D W I = \frac{B 3 - S W I R}{B 3 + S W I R}$
	植被近红外反射指数(near-infrared reflectance of vegetation, NIRv)	$N I R v = B 8 (\frac{B 8 - B 4}{B 8 + B 4})$

①SWIR表示B11和B12这2个波段; RE表示B5, B6, B7, B8A这4个波段。

新窗口打开| 下载CSV

2.2 随机森林模型训练与最优决策规则重建

基于构建的训练样本和特征集,应用R语言的randomforest包进行随机森林模型训练。随机森林是一种集成随机构建的决策树的强分类器,其通过随机选择特征子集和样本子集构建决策树,并通过投票综合这些决策树的结果以提高分类精度^[23]。本研究中,随机森林模型决策树数量设置为500,构建决策树每个结点的备选分类特征数目为7。

已训练的随机森林模型内部决策规则需要重构才能表征投票过程。考虑到最简单的决策规则的由特征、关系和阈值构成,因此可遍历随机森林模型的各决策规则,并按照特征与阈值的关系分组,将决策规则离散化为特征频数和相应特征的阈值分布。参照遗传算法,按照特征频数随机初始化种群,得到一系列特征组合及其对应的阈值区间,然后依据遗传算法求解各个特征组合对应的最优阈值。其中,遗传算法参数保持默认设置,即种群个体数目为200,最大迭代次数为100,变异率为0.01。各个特征对应阈值的最大和最小值通过获取的阈值分布估算得到。遗传算法的适应度函数采用Jaccard指数。Jaccard指数可以衡量2个集合之间相似性,公式为A∩B/A∪B,其中A和B分别表示盐地碱蓬的分类结果与参考分布^[21]。Jaccard指数值为1时表示分类结果与参考分布完全相同,为0则表示分类结果与参考分布完全不同。对于得到的一系列解集,按照Jaccard指数值进行排序,优选特征数量少、准确度高的特征及其阈值作为潜在解。通过局部分类结果筛选选定的潜在解,获取经过重构的决策规则。

综上,从已训练的随机森林模型中提取决策规则识别盐地碱蓬工作流程如图2所示: ①观测最佳时相确定; ②分类特征集构建; ③随机森林模型训练; ④遍历模型内部决策规则; ⑤决策规则重构。

图2

新窗口打开| 下载原图ZIP| 生成PPT

图2 从已训练随机森林模型中提取决策规则识别盐地碱蓬工作流程

Fig.2 Workflow of extracting decision rules from trained RF models to derive Suaeda Salsa map

2.3 决策规则的迁移应用

为获取研究区内不同年份的盐地碱蓬分布数据,对已构建的决策规则进行阈值优化是必要的。这是因为不同年份遥感影像成像时间不同,气候条件、卫星姿态等差异会导致同一地物的反射率发生变化。因此,在前述的遗传算法基础上,根据土地覆被类型未发生变化的样本,对决策规则阈值进行优化。在此基础上,通过目视校验方法进一步改进盐地碱蓬识别结果,提高分类结果的可用性。

2.4 精度评价

使用混淆矩阵评价盐地碱蓬识别决策规则的精度,包括总体精度 (overall accuracy,OA) 、生产者精度 (producer accuracy,PA) 和用户精度 (user accuracy,UA) 。验证样本为2019年采集,包括数目相等的碱蓬样本与非碱蓬样本共计400个,对于其余年份的验证样本,利用现有样本及高分辨率的谷歌地球影像进行采集且确保每年样点数量保持一致。各评价指标公式为:

(1)

O A = \frac{\sum_{i = 1}^{n} p_{i i}}{N}

(2)

P A = \frac{p_{i i}}{p_{+ i}}

(3)

U A = \frac{p_{i i}}{p_{i +}}

式中: n为类别数; N为验证样本点数; p_ii为第 i 行和第i列中被正确分类的样本数; p₊_i为第i列的总和; p_i₊为第i行的总和。

3 结果与分析

3.1 盐地碱蓬遥感识别决策规则的构建与应用

针对2019年盐地碱蓬构建的随机森林模型,重构得到其遥感识别规则为B3/B4<0.90 & B5/B3≥1.46。依据验证数据,该决策规则OA达到94.0%,碱蓬类别的UA为89.6%,PA为99.5%(表3)。重构的决策规则在2017—2022年盐地碱蓬提取结果见表4。对比不同年份下盐地碱蓬提取精度,本研究重构的决策规则在2020-09-18,2021-09-23及2022-09-18影像中提取效果最好,均达到了95%以上的精度。同时,盐地碱蓬提取时相跨度达6 a,而盐地碱蓬提取效果未出现较大差距,这表明重构后的决策规则对时间敏感性不强,其对盐地碱蓬的提取具有较好的适用性。

表3 混淆矩阵与精度分析

Tab.3 Classification confusion matrix and precision analysis

类别	碱蓬	非碱蓬
碱蓬	199	23
非碱蓬	1	177
UA/%	89.6	99.4
PA/%	99.5	89.0
OA/%	94.0

新窗口打开| 下载CSV

表4 不同年份对应决策规则的分类精度结果

Tab.4 Classification accuracy results of decision rules in different years

时间	阈值信息	OA/%
2017-09-29	B3/B4<0.98 & B5/B3≥1.16	86.0
2018-09-04	B3/B4<1.0 & B5/B3≥1.20	92.5
2020-09-18	B3/B4<0.88 & B5/B3≥1.46	96.8
2021-09-23	B3/B4<0.83 & B5/B3≥1.46	97.3
2022-09-18	B3/B4<0.96 & B5/B3≥1.26	98.5

新窗口打开| 下载CSV

应用迁移的决策规则,得到2017—2022年辽河口盐地碱蓬时空分布见图3。该地区盐地碱蓬主要分布于辽河口两岸,呈条带状分布; “退养还湿”工程实施后,盐地碱蓬面积先衰减后再迅速扩张。根据该地区盐地碱蓬面积统计结果(表5),2017—2019年间,该地区从1 154.69 hm²衰减至658.15 hm²。此后,盐地碱蓬斑块迅速向退养区域扩张,面积从2019年的658.15 hm²快速恢复至4 221.75 hm²,表明“退养还湿”工程在盐地碱蓬恢复过程中起到了积极作用。

图3

新窗口打开| 下载原图ZIP| 生成PPT

图3 2017—2022年辽河口盐地碱蓬时空变化

Fig.3 Temporal and spatial distribution changes of Suaeda Salsa in Liaohe estuary from 2017 to 2022

表5 2017—2022年辽河口盐地盐地碱蓬面积

Tab.5 Area information of Suaeda salsa in Liaohe estuary from 2017 to 2022 (hm²)

年份	2017年	2018年	2019年	2020年	2021年	2022年
面积	1 154.69	732.43	658.15	1 380.20	3 048.38	4 221.75

新窗口打开| 下载CSV

3.2 最优决策规则的可靠性分析

本研究从已训练随机森林模型中重构得到的决策规则为B3/B4<0.90 & B5/B3≥1.46,其不仅能够有效获取盐地碱蓬的分布信息,还可帮助人们理解盐地碱蓬的遥感识别机理。重构的决策规则信息来自绿(B3)、红(B4)以及红边(B5)波段,其中B3/B4<0.90对盐地碱蓬、沿岸滩涂及泥沙含量较大的水体均有积极响应(图4(c)),而B5/B3≥1.46对包含盐地碱蓬、芦苇在内的众多植被较为敏感(图4(d))。

图4

新窗口打开| 下载原图ZIP| 生成PPT

图4 决策规则内部各特征分类结果

Fig.4 Classification results of each elements of the decision rule

通过案例区域对比本研究重构的决策规则和SSVI的有效性。根据文献,利用SSVI提取2019-09-23 GF-1 WFV影像中的盐地碱蓬,指数阈值为0.5,提取结果见图4(f)。依据验证数据(表6),该提取结果OA为83.3%,碱蓬类别UA为96.5%,PA为69.5%。与本研究构建的决策规则相比,SSVI表现较差,尤其是在盐地碱蓬-芦苇混合生长区域。二者对光滩、如果表面及陆地植被等土地覆被类型存在一定错分问题,但SSVI漏分问题较为严重。

表6 SSVI提取盐地碱蓬结果的混淆矩阵

Tab.6 Confusion matrix of extraction results of Suaeda Salsa using SSVI

类别	碱蓬	非碱蓬
碱蓬	139	5
非碱蓬	61	195
UA/%	96.5	75.9
PA/%	69.5	97.5
OA/%	83.3

新窗口打开| 下载CSV

SSVI在多期影像中的盐地碱蓬提取精度见表7。对比表4可以看出,重构的决策规则对盐地碱蓬提取效果更好,除2017-09-29外,其在不同时期影像的盐地碱蓬提取精度均在90%以上且变化幅度较小。这说明本研究重构的决策规则其时间敏感性不强,具有较好的适用性。综合来看,本研究重构的决策规则对盐地碱蓬的提取效果优于SSVI。

表7 2017—2022年SSVI提取盐地碱蓬的阈值与精度

Tab.7 Threshold and accuracy of SSVI extraction Suaeda Salsa from 2017 to 2022

时间	阈值信息	OA/%
2017-09-20	SSVI > 4	73.0
2018-09-27	SSVI > 0.3	84.0
2020-09-30	SSVI > 1	71.5
2021-09-01	SSVI > 3	66.8
2022-09-01	SSVI > 5	79.5

新窗口打开| 下载CSV

3.3 辽河口盐地碱蓬的时空变化分析

应用质心分析^[24]进一步讨论2017—2022年辽河口盐地碱蓬的时空变化特征。根据图5可知,2017—2022年的盐地碱蓬质心坐标发生了明显迁移,其质心迁移轨迹与盐地碱蓬年际面积变化趋势相吻合。在盐地碱蓬缓慢衰退阶段(2017—2019年),其质心先向东南方向偏移,后向西方向迁移,该质心迁移轨迹表明辽河口北部盐地碱蓬在此阶段萎缩严重,泥沙淤积可能是造成此现象的主要原因;

图5

新窗口打开| 下载原图ZIP| 生成PPT

图5 2017—2022年辽河口国家自然保护区碱蓬的质心迁移情况

Fig.5 Centroids shifting of Suaeda Salsa in Liaohe estuary from 2017 to 2022

在盐地碱蓬快速增长阶段(2019—2022年),其质心首先向西北偏移,然后又折向东南,最后再向西南方向迁移。结合表4,盐地碱蓬向退养区域的扩张是导致其质心向西南迁移的主要原因,表明辽河口西岸养殖池设施的拆除为盐地碱蓬提供了更多生存空间。

3.4 盐地碱蓬、芦苇与养殖池的转化案例分析

2017—2022年期间,辽河口盐地碱蓬生长状况经历了显著变化。在2017—2019年间,研究区东部(图6)临近海岸线生长的盐地碱蓬逐渐退化消亡(图6(a)—(c))。然而,在2020—2021年间,泥沙淤积使得盐地碱蓬逐渐恢复,原本的破碎斑块逐渐扩张。泥沙沉积形成的新栖息地由盐地碱蓬首先占据并建立优势群落,随后芦苇逐渐侵入(图6(d)—(e))。截至2022年,该地区盐地碱蓬恢复良好,但存在盐地碱蓬被芦苇取代的问题,而部分区域仅残留小斑块盐地碱蓬(图6(f))。

图6

新窗口打开| 下载原图ZIP| 生成PPT

图6 辽河口东岸盐地碱蓬转化示意图

Fig.6 Schematic diagram of succession of Suaeda salsa in the east bank of Liaohe estuary

在研究区西部(图7),2017—2019年间,受养殖等人类活动干扰,养殖池边缘生长的盐地碱蓬面积萎缩(图7(a)—(c)); 2020—2022年间,盐地碱蓬面积回升并随着退养而侵入废弃的养殖池(图7(d)—(f))。

图7-1

新窗口打开| 下载原图ZIP| 生成PPT

图7-1 辽河口西岸盐地碱蓬转化示意图

Fig.7 Schematic diagram of succession of Suaeda salsa in the west bank of Liaohe estuary

图7-2

新窗口打开| 下载原图ZIP| 生成PPT

图7-2 辽河口西岸盐地碱蓬转化示意图

Fig.7-2 Schematic diagram of succession of Suaeda salsa in the west bank of Liaohe estuary

以上2个案例揭示了盐地碱蓬的细节变化,说明“退养还湿”工程的成效显著,但也应注意芦苇扩张对盐地碱蓬的影响。盐地碱蓬作为迁徙水鸟的重要休息觅食地^[25],土地覆被类型转化可能会影响其栖息和觅食行为^[26],因此,通过对不同植被类型与候鸟行为之间的关系研究或可进一步揭示影响该地区盐地碱蓬转化的原因。

4 结论

本文以Sentinel-2影像为数据源,基于随机森林模型与遗传算法,重构了一套可自我解释的盐地碱蓬决策规则(B3/B4<0.90 & B5/B3≥1.46),后利用该分类决策规则高精度快速提取了2017—2022年辽河口盐地碱蓬分布信息,并得出以下结论:

1)本研究重构得到的盐地碱蓬决策规则对2019年辽河口国家自然保护区的盐地碱蓬OA达到94.0%,UA达到89.6%,PA达到99.5%。在2017—2019年间,该地区盐地碱蓬从1 154.69 hm²衰减至658.15 hm²,而在“退养还湿”工程实施后,盐地碱蓬向退养区域快速扩张,其面积也于2022年恢复至4 221.75 hm²。

2)相较于现有的SSVI指数,本研究重构得到的决策规则对盐地碱蓬的提取完整度更高,尤其是在缓解提取小型盐地碱蓬斑块时的漏分现象方面。

3)“退养还湿”工程的实施在辽河口国家自然保护区西部创造了适宜盐地碱蓬生长的大面积裸滩,这是该地区盐地碱蓬质心位置总体向西南方向迁移的主要原因。

4)对辽河口国家自然保护区盐地碱蓬的转化情况分析可知,盐地碱蓬可在浅滩快速生长扩张,但泥沙堆积导致的地形抬升会迫使其被芦苇演替。同时,盐地碱蓬在退养区域的快速扩张表明“退养还湿”工程成效显著。

参考文献

原文顺序

文献年度倒序

文中引用次数倒序

被引期刊影响因子

[1]

温广玥.

1997—2018年辽河口翅碱蓬生物群落时空变化特征研究[D]. 北京: 中国地质大学(北京), 2020.

[本文引用: 2]

Wen

G Y

Temporal and spatial variation of Suaeda salsa community in Liaohe River from 1997 to 2018[D]. Beijing: China University of Geosciences, 2020.

[本文引用: 2]

[2]

, Jin

, Chen

et al.

Areal extent,species composition,and spatial distribution of coastal saltmarshes in China

[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2021,14:7085-7094.