富营养化藻型水体甲烷扩散排放的卫星遥感监测方法

1.本发明属于遥感环境监测技术领域，具体涉及富营养化湖泊藻型水体甲烷扩散排放的卫星遥感监测方法。


背景技术：

2.甲烷是重要的大气温室气体，对气候变暖影响显著。目前大气甲烷浓度呈现快速增长趋势，但导致大气浓度变异原因还不明了，因此亟需开展甲烷不同排放源的研究。富营养化水体对气候变化敏感，初级生产力高，已被发现是大气甲烷重要的自然排放源。扩散是富营养化水体甲烷重要排放方式。富营养化水体地球化学循环过程较为复杂，基于传统人工调查等方法很难捕捉甲烷扩散排放的时空变化规律，也难以应用于区域尺度水体甲烷扩散排放监测，致使甲烷扩散排放量估算存在较大的不确定性。卫星遥感具有大面积同步连续观测优势的遥感技术可以克服传统观测方法的局限，为富营养化水体甲烷扩散排放研究提供关键技术支撑。


技术实现要素：

3.本发明的目的在于克服现有技术的不足，结合遥感影像，提供一种富营养化藻型水体甲扩散排放的卫星遥感监测方法，实现富营养化水体甲烷扩散排放的大面积、连续、同步观测。
4.本发明的上述技术目的通过如下技术方案实现：一种富营养化藻型水体甲烷扩散排放的卫星遥感监测方法，包括：获取研究区域卫星遥感影像数据、光合有效辐射数据、地表温度数据和实测甲烷扩散通量数据；利用卫星遥感影像数据建立漫衰减系数kd(490)和叶绿素a浓度的估算模型；以漫衰减系数kd(490)、叶绿素a浓度、光合有效辐射数据、地表温度数据和/或其不同数学变换形式的组合为自变量，实测甲烷扩散通量数据为因变量，利用随机森林算法建立预测模型并进行精度验证，选取使模型精度最佳的自变量建立甲烷扩散排放预测模型；利用所述甲烷扩散排放预测模型监测富营养化藻型水体甲烷扩散排放。
5.作为一种优选的实施方式，所述利用卫星遥感影像数据建立漫衰减系数kd(490)和叶绿素a浓度的估算模型的方法为：对卫星遥感影像数据进行预处理，分别获取瑞利散射校正后的遥感反射率r
rc
数据及精确大气校正后的遥感反射率r
rs
数据；利用r
rs
数据结合半分析算法建立漫衰减系数kd(490)估算模型；利用r
rc
数据结合实测叶绿素a浓度数据建立叶绿素a浓度估算模型。
6.作为一种优选的实施方式，利用645nm和859nm波段对应的遥感反射率r
rc
值建立叶绿素a的估算模型。
7.作为一种优选的实施方式，所述卫星遥感影像数据选用modis/aqua的l1a级数据，
光合有效辐射数据和地表温度数据选用modis对应的光合有效辐射par产品和地表温度产品数据。
8.作为一种优选的实施方式，将具有不同分辨率的遥感反射率与遥感产品数据重采样到相同的空间分辨率，并剔除云覆盖较高和太阳耀斑的影像，仅保留无云或者云覆盖低的影像用于计算。
9.作为一种优选的实施方式，所述不同数学变换形式的组合指对所述漫衰减系数kd(490)、叶绿素a浓度、光合有效辐射数据、地表温度数据四个自变量参数中的多个进行不同数学形式的组合，形成新的自变量参数。
10.作为一种优选的实施方式，所述不同数学变换形式的组合指对所述漫衰减系数kd(490)、叶绿素a浓度、光合有效辐射数据、地表温度数据四个自变量参数中的两个进行不同数学形式的组合，形成新的自变量参数。
11.作为一种优选的实施方式，所述不同数学变换形式的组合包括对所述四个自变量参数中的多个以乘积形式，或乘积后取对数的形式组合形成的新的自变量参数。
12.作为一种优选的实施方式，利用随机森林算法建立预测模型后，利用独立样本验证模型精度。
13.作为一种优选的实施方式，所述方法还包括，利用长时序的遥感数据结合所述甲烷扩散排放预测模型获取长时序的水体甲烷扩散排放时空分布结果。
14.本发明的方法利用遥感数据结合机器学习，通过卫星影像数据实现对目标水域大面积、同步、连续的水体甲烷排放监测，获取甲烷扩散排放的空间分布及其在不同时间尺度（天-月-年）分布，有效解决了水体甲烷排放时空异质性强的问题，减少富营养化水体甲烷排放的观测成本，并通过输入参数的研究及对其交互作用考虑，提高了预测模型的准确性，从而提高水体甲烷排放估算精度。
附图说明
15.图1为本发明方法的流程图。
16.图2为估算甲烷扩散通量(fm)的验证图。
17.图3为本发明遥感监测的2002-2020年太湖甲烷扩散通量(fm)的时间序列图。
具体实施方式
18.下面结合附图和实施例对本发明的技术方案做进一步说明。
19.选取太湖为研究区域。太湖(119.55-120.34e， 30.55-31.32n)水域面积为2338 km2，是我国第三大淡水湖，长期藻华富营养化问题困扰，蓝藻水华频繁暴发，是典型的富营养化藻型湖泊。
20.从nasa海洋水色网站(https://oceandata.sci.gsfc.nasa.gov/)下载太湖富营养化水域2002年7月-2020年12月modis/aqua的l1a级数据供遥感反射率的计算使用，并下载modis/aqua的水面日平均光合有效辐射par产品和地表温度产品(myd11a1)，其中日平均光合有效辐射par是模型的输入变量之一，地表温度作为表层水温（lst）变量来估算甲烷扩散通量(fm)。
21.影像数据预处理，主要包括辐射校正、几何校正和研究区裁剪，其中所述辐射校正
指辐射定标和大气校正，经seadas7.5处理得到经过瑞利散射等去除后的遥感反射率r
rc
，另外采用6sv模型对原始数据进行精确大气校正以获取遥感反射率r
rs
。为保证具有相同的空间分辨率，并支持下一步的计算，将具有不同分辨率的modis反射率与产品数据重采样到相同的空间分辨率(250m)。另外，剔除云覆盖较高和太阳耀斑的影像，仅保留无云或者云覆盖低的影像。
22.结合2012-2017年太湖星地同步数据，实测数据是甲烷扩散通量(fm)，基于同步卫星影像，反演/计算与实测甲烷扩散通量(fm)同步的kd(490)和chl-a数据，其中漫衰减系数kd(490)的遥感提取主要是根据huang等(2017)开发的半分析算法得到（huang,c.c.,yao,l.,huang,t.,zhang,m.l.,zhu,a.x.,&yang,h..(2017).windandrainfallregulationofthediffuseattenuationcoefficientinlarge,shallowlakesfromlong-termmodisobservationsusingasemianalyticalmodel.journalofgeophysicalresearch:atmospheres,122(13):6748-6763.），如下：（1）（2）（3）（3）（4）（4）（5）其中，b
bp
(645)为645nm的颗粒后向散射系数，rrs(645)与rrs(531)分别为645nm及531nm波段的遥感反射率，a(645)为吸收系数，bb(645)为后向散射系数，u(645)为中间变量，r
rs
(645)为恰在水面以下的遥感反射率。
23.进一步的，叶绿素a浓度遥感提取使用shi等(2017)等针对富营养湖区开发的算法（shi,k,,zhang,y.l,,zhou,y.q.,liu,x.h.,zhu,g.w.,qin,b.q.,&gao,g..(2017).long-termmodisobservationsofcyanobacterialdynamicsinlaketaihu:responsestonutrientenrichmentandmeteorologicalfactors.scientificreports,7:40326），该算法是利用实测叶绿素a数据和645nm、859nm波长的遥感反射率建立的拟合模型：（6）其中，r
rc
(λ)是经过瑞利校正λ波长处的反射率，则r
rc
(645)、rrc(859)分别为modis数据中645nm及859nm波段r
rc
值。
24.结合上述步骤，以遥感反演的水体叶绿素a浓度(chl-a)、表层温度(lst)、漫衰减系数(kd)、光合有效辐射产品（par）以及不同参数之间的乘积或乘积后取对数的组合为自变量，实测甲烷扩散通量为因变量，结合随机森林算法(参数配置为：bootstrap=true,max_features=0.1,min_samples_leaf=2,min_samples_split=4,n_estimators=100)，
以chl-a、lst、kd、par、chl-a
×
lst、log
10
(chl-a
×
lst)、chl-a
×
kd、lst
×
kd、log
10
(lst
×
kd)、kd×
par为最终输入变量，进行拟合回归，构建甲烷扩散排放的预测模型并采用独立样本数据验证预测模型的精度，如图2所示；利用长时序的遥感影像数据，计算得到2002年7月到2020年12月日尺度太湖甲烷扩散排放通量，如图3所示。