基于CDOM吸收特性的东北地区水库DOC浓度遥感估算*

赖锋发，杨 倩，宋开山，尚盈辛

(1：吉林建筑大学测绘与勘查工程学院，长春 130018) (2：中国科学院东北地理与农业生态研究所，长春 130102)

有色溶解性有机物(CDOM)又称黄色物质，是溶解性有机质(DOM)中有色的部分，占溶解有机碳(DOC)的10%～90%[1]，CDOM是水体中主要的吸光物质，其光学特性比较单一，在水中主要表现为吸收特性，其散射可以忽略不计。CDOM从紫外到蓝波段的吸收非常强烈，在阻止紫外线对水生生态系统不利影响的同时，还改变着水下光强的分布，影响水生生物群落的结构和功能[2]。此外CDOM组成成分复杂，主要由富里酸、腐殖质、芳烃聚合物等物质组成，通过陆源和内源两个途径输入到水体中，陆源通常是陆生植被降解后形成的腐殖质，内源通常由水体中藻类和水生植被降解所形成[3]。CDOM的来源、性质以及空间分布在不同水体中存在着显著差异，CDOM光谱吸收特性能有效表征其物质组成和来源[4-6]，并在一定程度上指示水体的污染状况。

DOC是内陆水体最主要的碳组分之一，也是DOM的主要成分，在水环境生态系统中扮演着多种角色，对全球碳循环过程及气候变化都具有至关重要的作用[2]。正确认识河流溶解有机碳的行为特征，了解内陆水体中有机碳的储量和规律，对区域碳循环和碳估算具有重要意义。近年来，自然水体中的CDOM已成为研究热点，国内外学者在大洋Ⅰ类水体[7]、河口近岸水体[8]以及内陆自然水体[9]中均发现CDOM与DOC有良好的关系，并利用MODIS、SeaWiFS和Landsat系列卫星构建了CDOM吸收系数与遥感反射率反演模型，根据DOC与CDOM吸收的关系完成对内陆水体DOC浓度进行估算。

水库作为极其重要的Ⅱ类水体具有复杂的光学特性，其中CDOM的浓度和组成差异直接导致了水库之间水下光场性质和结构的差异性，对于整个内陆水体碳循环系统具有重要贡献。目前我国内陆水体CDOM和DOC的相关研究在南方地区典型湖库已相当系统，东北地区水库CDOM光学特性的研究近年来也在不断深入。以往研究表明东北地区内陆水体的DOC浓度普遍高于南方地区内陆水体，在我国内陆水体碳循环作用中占有重要位置[10]。随着东北老工业基地城市的快速发展，大量人类污染物和氮磷营养盐排入水库导致水体生态系统结构和功能发生变化，因此东北地区水库CDOM的光学特性研究亟待解决。利用CDOM作为水体生态环境指示因子、水体碳循环示踪物对表征水库水生生态系统有着重要的作用，这对于构建东北地区水库水质参数反演监测理论模型以及水资源管理具有重要的国家战略意义。本研究以位于松花江上的东北地区三大水库中的丰满水库和白山水库为研究对象，其水系一脉相承，两水库区域内具有显著人类活动及地表径流输入特征。通过解析两个水库的CDOM吸收与DOC浓度特征，利用Landsat系列遥感数据进行两水库的CDOM年际时空反演研究，结合CDOM吸收特征与DOC浓度的相关性，对水库DOC浓度进行时空演变分析，并解析DOC浓度变化的驱动因素，以期进一步揭示CDOM在东北地区内陆水库型湖泊中的光学特性变化，为水质参数反演监测提供理论依据。

1 数据与方法

1.1 研究区概况

白山水库和丰满水库位于吉林省东南部地区，是松花江南源上建立的2个梯级水库。其水系一脉相承，是东北地区排名前三的大型水库，年内气温变化较大，封冻期长，降水量由东南山区向西北平原递减，来水相对集中在春季和夏季两个汛期。白山水库位于吉林省桦甸与靖宇两县交界处，于1992年正式完工，全河段位于长白山山脉中，层峦迭嶂，河谷狭长弯曲，水深能达80 m以上。丰满水库位于吉林市东部，是1937年丰满水电站建成后形成的人工湖泊，平均水深30～40 m，水库总库容107.8亿m3，库区水量丰富，落差集中，以发电防洪为主，兼具旅游开发功能，流域内森林覆盖率由东南向西北呈递减趋势，林区随着海拔的降低转变为适宜耕种的农田。

1.2 野外采样

本研究于2015年5月、2018年10月、2019年6月对白山水库采样，2015年8月、2018年10月、2020年10月对丰满水库采样，采集丰满水库33个表层水样，白山水库76个表层水样，5次共采集109个有效采样点(图1)。现场利用YSI多参数水质仪(YSI EXO2)测定温度、pH、溶解氧DO、盐度、浊度、叶绿素等。在每个采样点用手持式GPS定位仪记录样点的位置和海拔高度。水样采集水面表层以下约0.2 m的水体，使用预清洁的PET塑料瓶收集约3 L水样，当天使用混合纤维滤膜(半岛0.45 μm)过滤水样，将初步过滤的水样放入车载冰箱冷冻保存，带回实验室进行测试。

图1 研究区采样点分布 (a：丰满水库；b：白山水库)Fig.1 Samplingsites distribution of study area (a: Fengman Reservoir, b: Baishan Reservoir)

1.3 实验室DOC及CDOM参数测定及计算

DOC浓度是由初步过滤水样通过岛津总有机碳分析仪(TOC-VCPN)测试得出，CDOM吸收系数的测试是将初步过滤水样使用聚碳酸酯膜(Whatman 0.2 μm)再次过滤得到用于CDOM测试的水样，以超纯水作为参照，利用岛津紫外-可见光分光光度计(UV-6200)测量CDOM滤液在200～800 nm处的吸光度，利用公式1得到CDOM的吸收系数：

aCDOM(λ)=2.303A(λ)/L

(1)

式中，aCDOM(λ)表示CDOM的吸收系数(m-1)，L为光路路径1 cm，本文CDOM浓度使用355 nm处吸收系数来表示。

SUVA254(比紫外吸光度)表示CDOM在波长254 nm处的吸收系数与DOC浓度的比值，用来表征单位DOC浓度的CDOM对光的吸收能力，其大小反映水体中有机质的芳香性和腐殖化程度，SUVA254越高其芳香性和腐殖化程度越高[11]。

M值表示aCDOM(250)/aCDOM(365)的比值，用来表征CDOM分子量的大小，M值与CDOM相对分子质量的大小呈反比，M值越大CDOM分子量越小，腐殖酸的比例越低[12]。

CDOM吸收光谱特性通常用公式(2)表示，采用最小二乘法对275～400 nm之间波段进行非线性拟合得到不同光谱斜率值。光谱斜率S275-295表示CDOM在275～295 nm波段处吸收光谱的斜率，S350-400表示CDOM在350～400 nm波段吸收光谱的斜率，SR表示光谱斜率比，为S275-295/S350-400(λ0为440 nm)：

aCDOM(λ)=aCDOM(λ0)exp [S(λ0-λ)]

(2)

吸收光谱斜率S值的大小与DOC的来源和组成有关，当DOC以陆源即流域内土壤腐殖质淋溶析出为主时S值偏小，当DOC以内源即水体中生物自身分解为主时S值偏大。S275-295能表征富里酸和腐殖酸的比例差异，斜率S越高表明其富里酸含量越高，斜率S越低表明腐殖质含量越高[12]。

1.4 遥感影像获取与模型构建

由于白山水库和丰满水库较长河段都在山谷中蜿蜒，水面宽度大多不超过500 m，Landsat系列卫星能提供30 m的空间分辨率和较好的时间序列支持，因此本研究使用Landsat系列卫星进行反演。遥感影像数据利用Google earth engine(GEE)平台提供的地表反射率数据集“USGS Landsat 5/7/8 Surface Reflectance Tier 1”。地表反射率是将已经通过大气校正的地面反射辐射量与入射量之比，反映地表对太阳辐射吸收和反射的能力，反射率越大，地表吸收的太阳辐射就越小。其中“Landsat5/7 Surface Reflectance Tier 1”数据集是Landsat 5 TM和landsat7 ETM+的原始数据利用LEDAPS程序根据6S辐射传输模型进行处理得到地表反射率产品，“Landsat 8 Surface Reflectance Tier 1”使用LaSRC程序进行大气校正后生成。其中Landsat 5 TM拥有7个波段，包含3个可见光、1个近红外、2个短波红外和1个热红外波段，Landsat 7 ETM+增加一个15 m分辨率的全色波段并将热红外分辨率提高至30 m，Landsat 8 OLI在Landsat 7 ETM+的基础上新增蓝色波段(0.433～0.453 μm)和用于云检测的短波红外波段(1.36～1.39 μm)。

通过GEE平台挑选出采样当天或者相近时间无云影像，提取出采样点位置处的遥感反射率，从白山水库和丰满水库分别随机抽取三分之二的样本进行模型建立，剩余三分之一的样本用于模型精度验证，通过回归分析将相关系数较高的波段或波段组合作为自变量进行拟合分析，建立CDOM遥感反演模型。基于构建的遥感CDOM波段组合模型，估算2000-2020年白山水库和丰满水库时空变化情况(除去2012年无合适遥感影像)，利用GEE筛选出无云非结冰期184景Landsat系列卫星影像，影像分布如图2所示，为减少因不同月份可利用的遥感影像数量分布不均的误差，在计算时空分布时，首先计算非冰期春、夏、秋季影像CDOM的均值，再以非冰期不同季节平均值计算CDOM年均值。

图2 遥感影像时间分布(a：丰满水库；b：白山水库)Fig.2 Time distribution of remote sensing image (a: Fengman Reservoir, b: Baishan Reservoir)

1.5 其他数据

本文使用的DEM数据为ASTER GDEM V3数字高程数据集(https://search.earthdata.nasa.gov/)，该数据集空间分辨率为1弧秒(约30 m)，参考大地水准面为世界大地测量系统(WGS84)和1996年地球重力模型(EGM96)大地水准面投影。土地利用数据来自于GlobeLand30 V2020版(http://globeland30.org/)，气象数据来自国家气象科学数据中心(http://data.cma.cn/)中国地面气候标准值月值数据集，收集其中降水、温度信息；归一化植被指数(NDVI)来自中分辨率成像光谱仪(MODIS)数据集MODIS/Terra Vegetation Indices Monthly L3 Global 1km SIN Grid V061，分辨率为1 km的月均值数据；人口、化肥使用量、工业废水、生活废水等统计数据来自吉林省统计局公布的统计年鉴。

1.6 统计分析和随机森林模型

本文采用IBM SPSS Statistics 26软件进行统计分析，包括方差分析(ANOVA)用于两个及两个以上样本均数差别的显著性检验，分别用**(P<0.01)和*(P<0.05)表示，Pearson相关性分析、Spearman相关性分析、拟合分析、均方根误差(RMSE)和平均相对误差(MRE)来分析估计值和实测值是否一致，其中RMSE和MRE的计算公式为:

(3)

(4)

将气象数据和NDVI月均数据计算为年均值，结合水体所在行政区统计年鉴的人口、污水排放量、化肥使用量等多种类型驱动因子，通过Spearman相关性分析计算单个因子与CDOM的相关性，筛选出相关性显著的因子(P<0.05)作为自变量，将年均CDOM值作为因变量，基于多元回归分析方法[13]，利用SPSS软件对丰满水库和白山水库的自然和人为因素进行回归分析，最后根据标准化回归系数和回归方程确定不同因子贡献度。

利用R语言Trend包中的mk.test函数实现Mann-Kenddall(M-K)检验，M-K趋势检验法是一种非参数检验，它不需要数据服从特定的分布，是一种常用的趋势检验方法，其中Uf的值大于零表示序列呈现上升趋势，反之呈下降趋势，当超过显著水平临界线时表明上升或下降趋势明显。

随机森林模型由Breiman等[14]基于bagging算法构建，近年来也逐步被应用于水质参数反演[15-16]，本文利用R语言中randomForest包实现随机森林模型，其中最重要的参数为mtry和ntree，mtry表示指定节点中用于二叉树的最佳变量个数，ntree表示随机森林所包含的决策树数目。将Landsat 8 B1～B7这7个波段的地表反射率作为输入自变量，aCDOM(355)作为因变量，随机挑选三分之二采样点用于模型构建，三分之一用于检验。

2 结果

2.1 CDOM吸收系数特征分析

CDOM吸收系数从紫外到可见光的范围内大致呈指数递减(图3)，以355 nm处吸收系数表征CDOM浓度，其中丰满水库CDOM浓度均值为(4.73±1.37) m-1，白山水库为(7.18±2.75) m-1，总体上白山水库CDOM浓度高于丰满水库，两个水库样点除数值上有一定差距外，吸收特性一致表现为紫外波段吸收强烈，在275～295 nm左右呈现出一个吸收肩，这与施坤等[17]的结果一致，丰满水库和白山水库的吸收肩斜率S275-295变化范围(均值)分别为0.0157～0.0234 (0.018±0.002)和0.0154～0.0206 (0.017±0.0015) nm-1，丰满水库的均值略高于白山水库，反映丰满水库CDOM中富里酸占比高于白山水库。

图3 CDOM吸收系数曲线(a：丰满水库；b：白山水库)Fig.3 CDOM absorption coefficient curve(a: Fengman Reservoir, b: Baishan Reservoir)

从图4可以看出，白山水库与丰满水库的SUVA254均值分别为(6.63±0.73)、(5.21±0.59) L/(mg·m)，白山水库和丰满水库的S275-295均值分别为0.017±0.001和0.018±0.002，M值均值分别为6.31±0.52和7.29±1.0，CDOM各项吸收参数显示白山水库腐殖酸含量略高于丰满水库，两个水库CDOM组分和来源具有差异性(P<0.05)。

图4 白山水库和丰满水库CDOM吸收参数M值、SUVA254和S275-295Fig.4 CDOM absorption parameters M, SUVA254 and S275-295 values of Baishan Reservoir and Fengman Reservoir

2.2 CDOM建模与验证

通过CDOM吸收曲线可以看出CDOM在短波段更为敏感，因此选择CDOM常用建模波段蓝波段(430～450 nm)、绿波段(530～590 nm)、红波段(640～670 nm)，利用波段比值的经验方法进行建模估算CDOM。如表1所示，单波段反演aCDOM(355)与蓝波段的相关性最好(验证R2=0.32，RMSE=2.21 m-1，MRE=36.47%)，其次是绿波段(验证R2=0.19，RMSE=2.43 m-1，MRE=39.93%)，最后通过对比R2、RMSE选择最优模型——红/蓝，绿波段模型效果最好，验证R2=0.86，RMSE=0.97 m-1，MRE=15.18%。

表1 波段组合Tab.1 Band combination

2.3 基于CDOM的DOC浓度估算

CDOM与DOC具有较好的一致性，不仅体现在海洋、河口和海岸带，还体现在内陆湖泊河流中。本研究将实测CDOM和DOC数据分为三分之二用于建模(N=73)，三分之一用于验证(N=35)，通过回归分析得到DOC=0.478aCDOM(355)+2.515，对比丰满水库和白山水库DOC与aCDOM(355)呈现良好线性关系，其R2分别为0.79和0.92，与先前诸多研究表现一致，证实CDOM反演DOC的可能性。

随机将实测DOC数据三分之二用于建模，三分之一用于验证，即丰满水库22个点用于训练模型，10个点用于精度验证，白山水库51个点用于训练模型，25个点用于精度验证。通过分析丰满水库、白山水库采样点数据中CDOM与DOC的相关关系，得出aCDOM(355)与DOC呈显著正相关(图5a)，验证结果如图5b所示，方程斜率接近1，散点均匀分散于两侧，R2=0.78，RMSE=0.65 m-1，MSE=10.9%。

图5 CDOM吸收系数与DOC浓度模型构建和验证Fig.5 Construction and verification of CDOM absorption coefficient and DOC concentration model

2.4 CDOM与DOC的时空演变

利用GEE平台筛选出2000-2020年无云非结冰期184景Landsat系列卫星影像，基于CDOM波段组合模型反演CDOM值，通过CDOM与DOC较好的一致性，得到2000-2020年DOC空间分布平均值(图6，7)。白山水库年均DOC浓度范围(均值)为5.05～8.6(5.9±0.97)m-1，而丰满水库年均浓度(均值)为4.84～6.97(5.35±0.6)m-1，从整体上看白山水库和丰满水库DOC浓度一直上下波动但变化幅度不大，从时间上看年均最高值均为2017年，最低值均为2011年，白山水库在绝大多数年份均值高于丰满水库(除2009年)。从空间上看，高值区域的分布大致符合东部和南部高，西北部较低的情况，丰满水库干流和支流汇入的河段即水库的南部和东部区域DOC浓度相对较高，汇入水库之后DOC浓度都有一定程度的减小，白山水库DOC高浓度区域主要出现在上游段头道松花江的汇入段以及东部二道松花江汇入段，丰满水库DOC高浓度区主要出现在南部支流辉发河的汇入段以及东部蛟河支流的汇入段。

图6 2000-2020年白山水库DOC浓度时空分布Fig.6 Spatio-temporal distribution of DOC concentration in Baishan Reservoir from 2000 to 2020

图7 2000-2020年丰满水库DOC浓度时空分布Fig.7 Spatio-temporal distribution of DOC concentration in Fengman Reservoir from 2000 to 2020

通过M-K检验分析DOC年均值变化可以看出，UF均在显著区间内波动，白山水库和丰满水库的变化基本一致，呈现变化不显著，在2010年之前，变化呈现上下波动，2010年之后均表现出上升趋势(图8)。

图8 白山水库(a)与丰满水库(b)DOC年均变化的M-K检验(橙色虚线表示Z值(P=0.05))Fig.8 M-K test of average annual variation of DOC in Baishan Reservoir (a) and Fengman Reservoir (b)

2.5 驱动力分析结果

所有的驱动因素数据时间段均选择2000-2020年，分别计算出水库年均CDOM与气象和人类活动因子Spearman相关系数和显著性检验P值(表2)，通过构建多元线性回归模型，评价驱动因子变化对水库CDOM的贡献情况。白山水库年均CDOM变化回归方程由6个驱动因子组成，可解释其年均变化的83.1%(R2=0.831)，贡献率排前3位的因素是降水、生活污水排放量和磷肥的使用量，分别贡献54%、9.07%和8.73%，丰满水库年均CDOM变化回归方程由6个驱动因子组成，可解释其年均变化的79.4%(R2=0.794)，贡献率排前3位的是降水、氮肥使用量和人口，分别贡献60%、12%、7.2%。其中对白山水库和丰满水库CDOM影响最大的因子均为降水，分别占贡献率的54%和65%(图9)。年际变化研究中选取的影像以秋季居多，秋季白山地区正逢锋面过境以及受台风北上影响雨水极为丰富，白山水库和丰满水库两岸地形极为复杂，沟叉港湾众多，周边农田环绕，因此受降水和农田施肥的影响较为密切。

表2 水库CDOM与气候和人为因素的相关性Tab.2 Correlation between CDOM of reservoir and climate and human factors

图9 各驱动因子对白山水库和 丰满水库CDOM的贡献率 (*表示P<0.01，**表示P<0.05)Fig.9 Contribution rate of driving factors to CDOM of Baishan Reservoir and Fengman Reservoir (*indicates P<0.01, **indicates P<0.05)

3 讨论

3.1 CDOM吸收特征分析

CDOM各项吸收参数多年均值均显示白山水库的腐殖酸含量略大于丰满水库，即白山水库CDOM的分子量略大于丰满水库，并表现出良好的一致性。SUVA254与CDOM中芳香族成分有强相关性，更高的SUVA254反映芳香族含量更高(表3)，对比国内其他河流湖泊，白山水库和丰满水库中SUVA254值相对较高，即CDOM芳香族的含量较高且白山水库高于丰满水库。研究发现[18]，河流或者受河流影响较大的河口及沿岸区域，其S值通常较小，在0.011～0.018 nm-1之间，而白山水库和丰满水库的S275-295均值分别为(0.017±0.0015)和(0.018±0.002) nm-1，反映出水库受河流输入的影响较大，白山水库与丰满水库中aCDOM(355)的吸收系数与斜率S275-295呈现相反的关系，较低的S值对应较高的CDOM吸收系数，这同殷燕等在新安江水库[19]以及段洪涛等在太湖[20]的研究一致。

表3 相关文献中湖库CDOM吸收参数对比Tab.3 Comparison of CDOM absorption parameters in lakes and reservoirs in related literature

从年际变化反演和驱动力分析的结果看，CDOM浓度的变化主要受到人类活动和气候变化方面的影响，在各支流处汇入时浓度变化尤为明显。辉发河口与蛟河汇入口水质较差，在反演的年际变化中可以清晰地看出，支流汇入处CDOM浓度较大，支流给丰满水库库区带来大量污染物，但是由于丰满水库中下游段湖水较深、湖面逐渐变宽、水流速度逐渐放缓，在水体的自净作用下，污染物逐渐被降解稀释，因此越靠近下游，丰满大坝水质越好，CDOM浓度越低，而白山水库河道相对狭小，从支流汇入的高浓度CDOM逐渐被降解稀释的效果不如丰满水库明显。

CDOM各项吸收参数均表明白山水库和丰满水库CDOM偏向于陆源大分子腐殖酸。通常陆源有机物主要来源于土壤表层中生物的降解，通过降水的方式进入库区[21]，根据驱动力分析可以看出降水和污水排放量贡献率超过2/3，降水将水库周边林地、草地和农田等生产降解出的部分有机物冲刷进入库区，而随着城镇化的进程，白山水库和丰满水库周围常伴有城镇污水进入[22]，污水排放分别占白山水库和丰满水库CDOM贡献率的16.52%和8.6%，未处理或者处理不完全的污水注入库区导致CDOM浓度增加。

3.2 CDOM模型讨论

通过波段比值将遥感影像与地面信息联系起来是很常见的，例如用于叶绿素反演[27]、水体悬浮颗粒物反演[28]和监督分类[29]等，诸多学者研究表明[30-31]通过波段比值可以校正大气的一些影响，增强所需的信息或者减弱其他信息的干扰。本文采用红、绿、蓝3个波段建立波段比值模型，得到良好的效果(R2=0.86，RMSE=0.97，MRE=15.1%)，由于CDOM的光学吸收特性随着波长增加呈现指数型递减，理论上短波的CDOM浓度更加敏感，蓝波段是反演CDOM的理想波段，但是蓝波段大气校正结果通常误差较大，在内陆水体中易受其他水体组分影响，而CDOM在红波段的响应可以忽略不计，因此红波段的加入起到归一化的作用，减少大气校正误差的影响提升反演精度，绿波段的加入一定程度上能够减弱浮游植物和非藻类颗粒物在蓝波段的吸收对CDOM反演造成的影响。

利用遥感影像对CDOM进行估测不仅仅只有波段比值方法，还有机器学习算法[16]、QAA算法[32]和矩阵反演法[33]等，但这些方法同样具有不足之处，比如计算复杂、计算量大或者模型构建复杂需要更多实测固有光学特性参数进行模型构建，而波段比值法进行反演建模的优点是算法简单、便于快速计算实现，能够一定程度上消除大气影响，在面对获取的遥感影像具有薄雾干扰时，波段比值模型表现出更好的抗干扰能力，缺点在于得出的反演模型只能用做特定的研究区域，普适性较差。本文尝试使用机器学习中随机森林算法对白山水库和丰满水库CDOM进行估算，在两水库采样点中随机挑选72个采样点作为训练集用于模型构建，36个采样点作为验证集用于模型检验，将Landsat 8 B1～B7这7个波段的地表反射率作为输入自变量，aCDOM(355)作为因变量，采用R语言的随机森林包实现反演。随机森林在验证精度上和波段比值模型相似(表4)，虽然该算法的相关性较高，但由于样品数量较少，会对随机森林算法的准确性产生影响，因此不具有可行性。在今后研究中将进一步进行补充实验采样，尝试在大尺度地区进行机器学习算法模型的构建，最终构建具有东北地区普适性的精确CDOM遥感反演模型和DOC估算模型。

表4 经验模型与随机森林模型对比Tab.4 Comparison between empirical model and stochastic forest model

3.3 DOC估算的不确定性

我国内陆水体中DOC的浓度受区域地理条件、人类活动以及季节的影响较大，赵莹[34]、尚盈辛等[35]研究表明各类内陆水体在不同的地域、季节以及人类活动的影响下，DOC与CDOM的相关性具有很大的不确定性，因此需要在不同地区建立DOC与CDOM的相关性模型，为进一步精确估算DOC浓度提供理论数据支撑。

由于DOC并没有明显的光谱特性，因此DOC的遥感反演通常是依靠间接的形式进行估算[36]，利用CDOM与DOC的强相关性估算DOC避免不了误差的累计传播。同时，DOC和CDOM的来源不一定是相同的，尽管CDOM与DOC呈现明显线性关系，但是根据来源的不同线性关系强弱会发生变化[37]，在河流交汇口，水体的垂直和水平交换速率较大，水体各项物质浓度不稳定，CDOM比DOC更容易发生光降解[38]，这些复杂的原因增加了DOC反演的不确定性。

白山水库和丰满水库位于我国东北地区，四季分明，冬季严寒漫长，夏季雨热同期，每年从12月上旬开始结冰到次年3月份陆续融化，在这期间由于冰雪覆盖的影响，导致卫星遥感影像很难对CDOM进行很好地估算，Song等[39]对东北松嫩平原40个湖泊进行研究，发现同一点的冰样CDOM浓度远小于水样，因此在估算的年均情况中去除了冬季的影响，容易导致年均值有一定程度的高估。在今后的研究中会进一步扩大东北地区水库的采样和研究范围，选取不同季节、不同类型的水库进行CDOM与DOC的相关性研究以及CDOM和DOC的遥感反演模型构建，为今后构建东北地区具有普适性的大尺度区域DOC估算和相关部门提供水库水质监测应用提供数据支撑。

4 结论

本文通过野外实测数据和Landsat系列卫星验证了通过遥感影像反演东北地区典型内陆水库——丰满水库和白山水库无冰期CDOM变化的能力(R2=0.86，RMSE=0.97 m-1)，反演2000-2020年(除去2012年无遥感影像)aCDOM(355)年均值，体现出Landsat系列卫星利用波段比值法反演丰满水库和白山水库CDOM的可靠性，根据CDOM与DOC的强相关性反演出水库DOC浓度年均值(R2=0.78，RMSE=0.65)，证明了遥感影像对DOC浓度变化监测的能力，增加了对东北典型水库碳循环的了解。

从反演结果来看，白山水库和丰满水库CDOM与DOC呈现出空间差异和季节差异。从空间上看，白山水库与丰满水库的高值区都分布在东部和南部支流的汇入处，越靠近下游大坝附近值越低。通过分析CDOM的吸收特性，发现白山水库的腐殖酸含量略高于丰满水库。总体来看，白山水库和丰满水库CDOM更倾向于陆源大分子的腐殖酸，降水、化肥使用量和污水排放量对白山水库和丰满水库CDOM时空变化起到主要影响作用。因此，在未来工作中将利用三维荧光技术结合平行因子以及分子结构分析等方法对东北典型水库CDOM空间分布特征及其来源做进一步遥感监测建模和时空演变分析。