湖泊遥感中溶解性有机碳的研究综述

摘 要：从研究文献现状、DOC光学特性、反演模型开发、卫星传感器及数据源等方向展开讨论。梳理了过去30多年来的发展历程及趋势，对未来DOC浓度遥感反演算法需要解决的问题和难题进行总结。

关键词：遥感；湖泊；溶解性有机碳；水质反演

湖泊生态系统中的溶解性有机碳（DOC）在全球碳循环中起着至关重要的作用[1]，但随着DOC浓度的增加，会导致湖泊有机酸增加，也会改变水下光场，降低水体透明度，从而影响水生生物和水生生态系统。因此，DOC监测对深入研究湖泊碳循环过程和估算湖泊碳收支以及研究DOC累积对湖泊生态系统的影响有着重要意义。

传统DOC监测采用TOC仪和分光光度计检测方法[2]，这些方法通常都是采用人工布点、收集水样、带回实验室进行仪器检测。这种方法往往会耗时费力，且数据量很少。随着水体富营养化和水环境问题日益突出，水体监测任务日益频繁，这种传统的方法不再适应现代社会的监测技术需求。卫星遥感具有快速、大尺度检测的技术优势，受到广泛青睐。湖面遥感是指在长距离、不直接接触目标表面的前提下，传感器或地测方法可以透过表面所进行反射、散射或传播的电磁波数据[3]，并对其资料加以获取、分析、处理研究、再利用，进而掌握并得到湖泊特征信息的一种技术[1，4]。近年来，遥感技术在水质监测中的应用日渐广泛[5]。湖泊DOC遥感监测也得到快速发展。因此本文从文献现状、DOC光学特性、反演模型开发、卫星传感器及数据源等方面，总结了湖泊DOC遥感的研究现状，这对于湖泊DOC遥感技术发展有着重要意义。

1 文献计量分析

中国知网（CNKI）作为本研究的文献数据库。以“溶解性有机碳”、“DOC”和“湖泊遥感”为主题词。1989—2022年期间国内外研究成果发表总计422篇，除去与主题不相关的文献，确定剩余398篇文献符合研究课题。其中，溶解性有机碳219篇，详细讲述溶解性有机碳成分的13篇，光学特性的5篇，遥感揭示时空变化31篇，遥感揭示DOC和其他生物关系128篇。其中中国科学院南京地理与湖泊研究所共发表该研究方向文献24篇，是目前国内研究机构中发表文献数最多的。

过去这30多年中，溶解性有机碳遥感反演研究不断深入，一开始研究范围多围绕于土壤，随着研究范围的逐渐扩大，湖泊水质遥感监测的需要，湖泊遥感中的溶解性有机碳浓度监测进一步发展并深入应用。从溶解性有机碳遥感反演的发文数量上来看，2010年至今是一个持续增加的趋势（图1），这种趋势可以分为三个阶段，一是2010—2013年，初次深入研究及稳步上升阶段，该阶段对其研究发展迅速；二是2013—2020年，这期间不断发生波动，有疫情等多方面因素影响其研究进展；最后一阶段是2020年至今，随着2020年双碳目标的提出，该研究方向开始呈现稳定发展并持续上升阶段。从这三个发展阶段看来，湖泊溶解性有机碳的遥感反演日趋成熟，同时也希望这可以为改善湖泊水生态环境提供部分技术参考和理论研究[6]。

2 DOC与CDOM的光学特性

湖泊在地球碳循环中释放和储存着大量的碳，每年湖泊碳排放和碳储存与地球海洋的数量相当， DOC是湖泊中最大的还原碳库[7]，它可以为湖泊中的浮游动植物、藻类、微生物等提供能量和碳，所以溶解性有机碳是湖泊碳循环中的关键部分[1，4]。

因为DOC是非光学活性的参数，因此利用具有光学活性的CDOM进行遥感反演是目前常用且有效的方法。虽然CDOM并不能代表整个DOC碳库，但许多研究均发现在某些海岸带水体、河口及内陆水域中，CDOM的光吸收与DOC浓度有很强的相关性[8]。目前，CDOM领域的研究主要是根据CDOM的特性来分析复杂内陆水域中溶解性有机物的相关特征，其光学特性主要体现在光吸收和荧光性上。

2.1 吸收特性

在溶解有机碳库中，有一组部分光学活性物质，并且可以强烈吸收紫外光（UV-A和UV-B），这些物质在波段中呈现出明显的黄色且在范围内吸收最不明显，于是将这类物质称为“黄色物质”。20世纪90年代后，国际上研究者便将黄色物质改称为“有色可溶性有机物（chromophoric dissolvable organic matter，CDOM）” [9]。CDOM的吸收主要发生在紫色和蓝色光区域，CDOM在可见的蓝色光区域的吸收与叶绿素a和总悬浮物（TSM）的吸收重叠，这对评估内陆水域的生物生产力有重要影响[9-10]。CDOM吸收紫外辐射光后，导致光漂白或是光降解反应， CDOM的彩色组分可以光解为无色组分，因此随着波长的增加，吸收作用会迅速减小。CDOM的光学活性物质主要是通过紫外-可见吸收来描述，CDOM在200～800 nm的波长范围内呈指数型， 200～800 nm之间呈指数级下降。通常选择a（440），a（355）或者a（350）吸收系数来表示短波长下CDOM的浓度，紫外线吸收波长越短，则CDOM的吸收系数与水中的DOC浓度之间的关系越好，一般来说，水中的DOC浓度与CDOM吸收系数之间可以建立一定的模型关系，从而估算出湖泊DOC的浓度[8，10]。

2.2 荧光特征

内陆水体的DOC浓度可以用来表示溶解性有机物（DOM）的浓度，而简单的特定波长吸收系数aCDOM可以表示DOM内部组分，为了区分水体中DOC的组分，通常可以测定水中CDOM的三位荧光光谱，从而分析组分来源于组成[11]。

利用CDOM的荧光特性：短波长的光作为激发波长，CDOM吸收较长的波长并发出荧光，利用PARAFAC（平行因子分析）对荧光光谱进行三维分析，能够客观定量和定性地分析CDOM含量和其他成分特性的变化。分析三维荧光法（TEE）方法可以解析水环境中CDOM的来源、降解、组分特性等[12]。 这种方法被广泛用于识别水生环境中CDOM的来源荧光成分和其他相关信息。CDOM的荧光部分主要是来自陆地输入物和腐烂的植物残体的蛋白质和腐殖酸部分，其中腐殖酸主要来自自由形成的细菌和水体中的微生物活动，但也来自工业、农业和生活废水。土壤类型、降雨量、径流量、植被类型、人为活动等都会对CDOM浓度及组成产生一定影响。现在，国际上许多研究表明，CDOM的荧光成分是通过三维荧光来描述的，分析海洋水域和内陆水域中CDOM的三位荧光光谱经常使用PARAFAC方法。PARAFAC方法提取3～13种荧光成分，同时还可以分析这些荧光成分在水体中的动态变化，进而得出它们的分布模式，将水质参数与CDOM吸收系数联系起来，并为利用水色遥感技术估计内陆水域的CDOM浓度提供了一个基础。

3 溶解性有机碳的反演模型

遥感技术可以获取湖泊中各种物质的光谱特征，所以可以进行与水体要素密切相关的要素进行反演，例如溶解性有机碳、透明度、浊度等[12-13]。

目前DOC的遥感反演有两种方法，一是直接反演法，根据实际测得的DOC浓度值与遥感波段反射率直接建立经验关系；二是是间接反演法，先通过遥感波段反射率反演CDOM浓度，再根据DOC浓度与CDOM浓度之间的关系计算DOC浓度。利用遥感数据进行湖泊水体要素反演的主要有三种模型：经验模型、半经验和半分析模型、人工智能模型[11，14]。

经验模型又叫统计模型，根据大量遥感数据（水表层以上或以下的反射率）和其相应现场同步测量的水体要素数值，使用统计学的方法建立遥感数据和实验数据之间的关系。这种模型的优点是算法更简单，需要的参数更少，当理论模型的参数很重、很复杂、更难获得时，这是最可取的。其缺点是空间和时间上的限制，低代表性，低稳定性以及重复利用性也不足。

随着DOC光学特性的了解和遥感机理研究的逐渐深入，经验模型逐渐放弃，具有机理模型的半经验算法被认可。半经验模型，也称为统计物理模型，是指在分析水体特征的光谱特性的基础上，在测量的水体特征数据和同时从遥感获得的光谱特征波段或组合之间建立统计关系。这种类型的模型的优点是，水体的光谱特征和统计分析的结合有一定的物理基础。其缺点是，模型的准确性必须很高，方法的应用很复杂且模型的通用性很低。

半分析模型是重要的分析模型。分析模型使用水的辐射传递模型来确定水体的光谱反射值和水体成分之间的关系，然后通过水体的反射光谱值反向推导计算得到反估水体各成分的含量。利用光学遥感数据反演水体成分，只提供了分析模型反演中的部分水体固有光学量，而反演过程中的大部分水体参数需要引入经验关系，所以在分析模型的推导过程中采用经验方程来引入，得到半分析模型[9]。半分析模型也是随着高光谱技术的不断发展而发展，因为半分析模型，一方面需要水体辐射传输方程中的参数比较多，必要使用具有多波段的高光谱遥感数据才能满足；另一方面因为高光谱数据可以提取水体中的特征波段，从而加强了水体要素反演的精度和稳定性。半分析模型是湖泊水体遥感反演模型日后发展的方向，随着大气校正模型的不断发展和大量水体固有光学数据的积累，半分析模型将在湖泊水体中发挥更显著的作用和优势。

除了上述所提到的经验模型、半经验/半分析模型以外，还有基于遗传算法、神经网络、支持向量机等人工智能模型。人工智能模型的优势在于它可以使用多种人工智能算法来模拟水体辐射传输的过程，可以达到快速反演多个水体要素。下面将国内外直接反演法与间接反演法做了总结（表1）。

3.1 DOC直接反演法

Arenz等[15]（1996）率先利用遥感技术建立了遥感反射率与DOC浓度之间的经验关系，并开发了一种利用波段匹配来反演DOC浓度的算法。陈楚平等[16]（2001）在珠江口进行了观测实验，得到了遥感反射率与实测DOC浓度的关系，最后确定了波段的最佳组合，构建了lg（DOC）=1.241 9 lg（R670/R412）-0.2614的模型用于该流域的DOC浓度反演。张运林[17]等（2005）通过研究太湖水体中DOC浓间的关系，得出：lg（DOC）=0.654（±0.012）lg[R670/R530] +1.007（±0.086）。Kust等[18]（ 2015） 使用MERIS数据，将CDOM、DOC 和 TOC等水体各个碳组分遥感反演模型。Cherukuru等[19]（2016）在澳大利亚莫尔顿湾基于野外数据，开发了DOC浓度相关经验算法，从R412/R488的比值中得出DOC浓度，建立了模型：lg（DOC）=282.01exp（-1.72（R412/R488）。

3.2 CDOM-DOC关系间接反演法

1993年，Frank E.H.等提出了根据荧光测量结果反演CDOM的光吸收活性，建立CDOM的荧光发射和DOC浓度之间的适当关系，从而快速获取DOC浓度的设想。Mannino等[20]（2008）取得了美国大西洋中部海湾的MODIS和SeaWiFS数据，同时利用DOC与CDOM吸收系数之间的经验关系建立了第一个CDOM吸收系数与遥感反射率的经验模型，从而可以反演该区域的DOC浓度。Fang Cao等[21]利用MERIS卫星数据，在墨西哥湾北部、切萨皮克湾河口和大西洋中部沿海水域进行算法研究，利用多光谱遥感反射率检索不同波长CDOM的光吸收活性，最终建立了CDOM吸收斜率与DOC浓度的关系式：DOC=aCDOM（300）/（exp（-15.05-33.95×S275-295）+exp（-1.502-104.3×S275-295）。起初大部分研究在欧洲海岸带、北美海湾等地区，还有少部分关于东太平洋的研究，因此国内对CDOM与DOC关系认识还不够充分，中国更是鲜少有对内陆水体的研究。白雁等[22]在我国黄东海海区进行试点研究，利用SeaWiFS数据资料，比较了SIO算法、Lee算法、GSM01算法的遥感反演结果，最终建立了SIO半分析遥感反演算法，即DOC=90.91×aCDOM（355）+84.5。虽然该算法需要进一步完善与验证，但总体来说可以较好地表述黄东海区域的趋势变化，也是国内较为成功的一项研究案例。

4 卫星传感器及数据源

随着技术的发展以及对湖泊DOC浓度动态监测的需求，监测的重要信息多源自于遥感数据，高光谱遥感技术的发展和高分辨率陆地卫星的应用也促进了湖泊DOC浓度的发展，同时可以将遥感反演算法成功运用到卫星影像上，实现湖泊DOC浓度在时空上的遥感估算，也为水质监测提供技术支撑[21]。遥感体系需要由传感器、遥感平台、地面站等构成，传感器可以远距离测量地物环境辐射。湖泊DOC浓度监测数据源可分为微波遥感数据和可见光数据[23]。