基于卫星遥感影像的藻型湖泊水体溶解二氧化碳估算方法与流程

1.本发明涉及环境工程与图像处理技术领域，具体涉及一种基于卫星遥感影像的藻型湖泊水体溶解co2估算方法。


背景技术：

2.水是生命之源，而湖泊是地球上最重要的淡水资源之一、是湖泊流域地区经济可持续发展和人们赖以生存的重要基础。同时，湖泊虽然面积远小于海洋，但其生态系统生产力很高，且与陆地生态系统物质、能量和信息交换强烈，是全球碳(c)循环的重要组成部分目前，我国内陆湖泊面临的一个主要问题是水体的富营养化，其重要特征是藻类，特别是蓝藻大量繁殖，异常生长，极易堆积，形成水华。一方面，藻类大量繁殖使得湖泊初级生产力大幅提高，消耗更多的co2；另一方面，导致湖泊富营养化的外源输入营养物质会加强有机质的分解产生更多的co2，同时堆积的藻类腐烂降解也会产生co2，进而抵消初级生产力消耗的co2。因此，在富营养化的藻型湖泊中，水体的溶解co2浓度存在着极高的时空分异。
3.目前，传统研究主要通过野外采样手段，在不同区域进行大规模采样分析，不仅费时费力，而且由于湖泊水体高度的时空异质性，有限的采样点数据无法代表全湖水体；或者有限的采样频次无法代表全年的湖泊水体。因此，目前关于湖泊碳排放的估算差异高达近10倍，而这些不确定性导致湖泊等内陆水体与大陆和全球尺度碳循环的界限和联系并不清晰。因此，如何准确评估湖泊co2排放量及其在全球碳循环中的作用是一个目前亟需解决的一个关键问题。
4.卫星遥感快捷、大范围和周期性的特点为湖泊水体碳通量研究提供了更好的研究手段。特别是随着拥有适合湖泊水体水色波段的卫星，如美国国家航空航天局aquamodis(全球1 天可以重返一次)、欧空局sentinel-3olci(目前已经发射2颗，全球2天可以重返1次)和我国高分系列卫星(gf-5有高光谱传感器)等已经发射。因此，开展基于卫星遥感的湖泊水体溶解co2浓度，对准确估算湖泊水气界面co2交换通量具有重要的研究意义和科学价值。


技术实现要素：

5.本发明的目的在于提供一种基于卫星遥感影像的藻型湖泊水体溶解co2估算方法，利用现有公开的湖泊水环境参数卫星遥感产品，实现准确地估算目标水域的水体溶解co2浓度及其空间变化状况，为区域碳排放和环境保护提供重要的科学依据与技术支撑。
6.为达成上述目的，本发明所采用的技术方案如下：
7.一种基于卫星遥感影像的藻型湖泊水体溶解二氧化碳估算方法，包括：
8.获取目标水域叶绿素a浓度chla、水温lst、光合有效辐射par、光合有效辐射漫射衰减系数k
d
(par)遥感产品数据；
9.采集目标水域表层水体样本，获取实测溶解co2浓度；
10.以获取的chla、lst、par、k
d
(par)数据为自变量，co2浓度为因变量建立回归方程；
11.基于所述回归方程结合待测湖泊遥感数据估算藻型湖泊水体溶解二氧化碳。
12.进一步的，所述遥感产品数据基于modis-aqua卫星数据反演获取。
13.进一步的，对获取的遥感产品数据进行水华像元和水生植被像元剔除处理；通过浮藻指数fai剔除水华像元，并通过掩膜处理剔除包含水生植被的像元。
14.进一步的，通过ph-碱度计算法获取溶解co2浓度。
15.进一步的，以
①
水体样本采样时间与遥感产品卫星观测时间时差不超过5小时；
②
采样点坐标中心像元周围3
×
3像元窗口中的一半像素具有有效数据，且有效数据的变异系数小于 10％为准则进行星地同步数据匹配，获得样点-像元配对数据集；基于样点-像元配对数据集中的数据建立回归方程。
16.进一步的，根据co2浓度从小到大的顺序将配对数据集排列，每间隔3个样本挑出一个样本作为验证样本，建立模型训练数据集和模型验证数据集，基于模型训练数据集建立回归方程后采用模型验证数据集进行数据验证。
17.进一步的，对lst、k
d
(par)产品数据进行对数化处理，建立chla、ln(lst)、ln(k
d
(par)、 par和ln(cco2)的回归方程，确定回归方程的系数。
18.进一步的，建立的回归方程如下：
[0019][0020]
式中，x1＝chla，x2＝ln(lst)，x3＝ln(k
d
(par))，x4＝par，a
0-a
14
为方程系数。
[0021]
进一步的，所述方法还包括，计算全湖或子区域的水体溶解co2浓度的平均值、标准差等统计指标，以作为长时序时空变化分析的数据支撑。
[0022]
进一步的，所述方法还包括，生成水体溶解co2浓度影像，获得其空间分布情况。
[0023]
本发明的算法原理如下：
[0024]
理论上，水体中溶解co2浓度主要受溶解度泵和生物泵驱动。溶解度泵取决于风速、温度、盐度和碳酸的离解。co2在水体中的溶解度与温度和盐度成反比。生物泵主要在表层水中通过浮游植物的光合作用与生物的呼吸分解作用共同影响溶解co2。浮游植物的光合作用从水体吸收溶解的co2，导致溶解co2降低；浮游植物与细菌等微生物的呼吸作用释放co2到水体，导致溶解co2升高。在富营养化的藻型湖泊水体中，浮游植物在水体中的物质能量循环中占主导地位，同时浮游植物还是富营养化程度的直接指示指标，间接反映了外源输入污染物，是影响水体中溶解co2的重要因素。
[0025]
叶绿素a作为浮游植物共有的色素，其浓度大小可以直接表征浮游植物生物量的多少。并且，叶绿素a浓度(chla)是水色三要素之一，是常规水色遥感产品。此外，对水体中各种生物化学反应都有重要影响的水温(lst)也可以通过遥感手段反演。光合有效辐射(par) 以及对应波段的水体漫衰减系数(k
d
(par))分别影响了水体接受到用于生化反应的能量强弱以及水下光场中能量的衰减快慢，同时作用于藻类光合作用与有机质的光解，因此本发明中也利用两者的遥感产品来解释水体中溶解co2的变化。综上所述，本发明基于卫星遥感的chla、 lst、k
d
(par)、par产品，通过多元二次多元式回归算法建立了这些自变量与水体溶解co2浓度的定量关系，实现了对水体溶解co2浓度的遥感估算。
[0026]
应当理解，前述构思以及在下面更加详细地描述的额外构思的所有组合只要在这样的构思不相互矛盾的情况下都可以被视为本公开的发明主题的一部分。另外，所要求保
(collection6)lst，原始分辨率为1km日间产品，重采样为250m分辨率与其他产品匹配； k
d
(par)产品使用huang et al.(2017)提出的半分析算法反演；par产品使用nasa obpg (https://oceandata.sci.gsfc.nasa.gov/)网站提供的modis-aqua卫星反演得到的日平均par 的l2级产品，原始分辨率为1km日间产品，重采样为250m分辨率与其他产品匹配。
[0041]
之后，以卫星chla、lst、k
d
(par)、par产品数据为基础进行太湖水体溶解co2浓度的估算具体处理流程如下：
[0042]
1)对遥感获取到的chla、lst、k
d
(par)、par产品数据进行预处理。由于在蓝藻水华区域无法反演有意义的chla与k
d
(par)参数，因此采用hu et al.(2010)提出的浮藻指数(fai) 来剔除影像数据中的水华像元(fai>-0.004)。此外，由于太湖东部区域存在大量水生植被，不适用本发明所述的方法，故进行掩膜处理剔除这部分数据。
[0043]
2)对lst、k
d
(par)产品数据进行对数化处理；
[0044]
3)以
①
采样时间与卫星观测时间时差不超过5小时，
②
采样点坐标中心像元周围3
×
3 像元窗口中的一半像素具有有效数据，且这些有效数据的变异系数小于10％为准则将太湖站样点的水体co2浓度与预处理后的遥感影像栅格数据像元进行匹配，获得样点-像元配对数据集。共获得样点-像元配对80对，根据co2浓度从小到大的顺序将配对数据集排列，每间隔3个样本挑出一个样本作为验证样本，得到配对数据数为54的模型训练数据集和配对数据数为26的模型验证数据集。
[0045]
4)通过模型训练数据集样本进行回归计算以确定方程系数，回归方程如下：
[0046][0047]
式中，x1＝chla，x2＝ln(lst)，x3＝ln(k
d
(par))，x4＝par，a
0-a
14
为方程系数，a0＝4.736， a1＝-0.0791，a2＝-0.412，a3＝-0.275，a4＝0.00383，a5＝0.00121，a6＝0.194，a7＝0.232，a8＝0.00106， a9＝-0.00434，a
10
＝0.00769，a
11
＝0.000781，a
12
＝-0.108，a
13
＝-0.0334，a
14
＝0.00675。
[0048]
利用训练数据集样本与验证数据集样本对模型的预测精度进行评估，结果如图1所示。模型的预测结果同野外实测co2浓度结果的均方根误差(rmse)与无偏差百分比差异(upd) 均较小，且散点均匀的分布在1：1线两侧，说明模型对于1-100μmol l-1
接近2个数量级范围内的co2浓度具有很好的预测能力。
[0049]
5)将完成预处理后的遥感影像栅格数据通过经纬网格关联起来，并以像元为单位将chla、 lst、k
d
(par)、par数据代入回归方程估算对应像元的水体溶解co2浓度(cco2，μmol l-1
)。
[0050]
6)计算太湖全湖以及五个子区域(西北湖区、梅梁湾湖区、贡湖湾湖区、湖心区、西南湖区)的水体溶解co2浓度的平均值、标准差等统计指标，以作为长时序时空变化分析的数据支撑。
[0051]
最终，湖泊水体溶解co2浓度影像数据后处理：包括生成水体溶解co2浓度影像，获得其空间分布情况。
[0052]
利用上述算法，对太湖水体溶解co2浓度进行自动计算，发现对于不同子湖区水域、不同水体特征的影像，算法都发挥了较好的作用，估算结果具体如图2-4所示，图5进一步对太湖五个子区域长时序日平均co2浓度变化曲线统计显示，以分析水体溶解co2的时空
分异规律。
[0053]
通过上述方法即可建立基于遥感影像的湖泊水体溶解co2估算模型，估算精度较高；据此，卫星遥感监测将在大型藻型湖泊水体碳循环及温室气体排放研究中具有较好的前景。本发明能够全面获得250m空间分辨率的湖泊水体溶解co2浓度空间分布，对评估湖泊水-气界面co2源汇功能具有重要借鉴和指导意义，可以为区域碳排放和环境保护提供重要的科学依据与技术支撑。
[0054]
虽然本发明已以较佳实施例揭露如上，然其并非用以限定本发明。本发明所属技术领域中具有通常知识者，在不脱离本发明的精神和范围内，当可作各种的更动与润饰。因此，本发明的保护范围当视权利要求书所界定者为准。
[0055]
参考文献
[0056]
hu,c.m.,lee,z.p.,ma,r.h.,yu,k.,li,d.q.,&shang,s.l.(2010).moderate resolution imaging spectroradiometer (modis)observations of cyanobacteria blooms in taihu lake,china.journal of geophysical research-oceans,115, 20
[0057]
huang,c.c.,yao,l.,huang,t.,zhang,m.l.,zhu,a.x.,&yang,h.(2017).wind and rainfall regulation of the diffuse attenuation coefficient in large,shallow lakes from long-term modis observations using a semianalytical model. journal of geophysical research:atmospheres,122,6748-6763
[0058]
shi,k.,zhang,y.,zhou,y.,liu,x.,zhu,g.,qin,b.,&gao,g.(2017).long-termmodis observations of cyanobacterial dynamics in lake taihu:responses to nutrient enrichment and meteorological factors.scientific reports, 7,40326
[0059]
xiao,q.,xu,x.,duan,h.,qi,t.,qin,b.,lee,x.,hu,z.,wang,w.,xiao,w.,&zhang,m.(2020).eutrophic lake taihu as a significant co2 source during 2000
–
2015.water research,170,115331。