由于洪泽湖水体和大气浑浊，常用的水色大气校正算法难以获取高精度的水体遥感反射率数据，因此，本文基于去除瑞利散射校正后的反射率（Rrc, Rayleigh-Corrected Reflectance），利用红光与短波红外的差值构建了适用于洪泽湖悬浮物的高精度算法（R2=0.64, RMS = 39.4%），该算法能够有效地应用于MODIS、VIIRS等传感器，在10-80mg/L的悬浮物范围具有较好的估算精度。

船只在Landsat数据的短波红外波段伤具有较高反射率，明显区别于水体，因此基于Landsat数据的短波红外波段，结合真彩色合成影像，对洪泽湖的采砂船只进行识别提取，获取了2002-2017年的洪泽湖采砂船数量。

VIIRS的夜间灯光数据能够获取地球表面的光信号，夜晚，洪泽湖面出现的灯光物体为非法采砂船聚集造成，基于梯度阈值分割的方法，对洪泽湖的采砂船聚集的灯光区域进行提取，获取采砂船灯光2012-2017年的采砂船像元分布。

3.1 非法采砂活动的影响

利用2012-2015年洪泽湖111个星地同步观测数据集，建立了基于MODIS Aqua数据的较高精度的水体悬浮物浓度估算模型(RMSE = 7.64mg/L，RMS = 46.9%）。将该算法成功应用到了2002-2015年1602景MODIS Aqua无云影像上，获取了洪泽湖的悬浮物14年来的时空分布。洪泽湖悬浮物浓度在2002-2011年全湖悬浮物浓度较低且变动较小，2012年之后湖体悬浮物浓度显著上升；同时，洪泽湖水体悬浮物浓度存在显著的年内季节差异，2002-2011年悬浮物浓度在夏季最高，春季最低；但2012-2015则展现出完全不同的季节变化，在秋、冬季高，春、夏季低。实际上2012年前后两个阶段，洪泽湖地区降雨和风速等气象数据并没有出现较大变化，这种年际和季节上的显著变化是由于2012年开始的非法采砂活动造成的。初步估算全湖悬浮物上升了21.54%，在12.5米常水位下，采砂活动每年约向洪泽湖贡献150.48亿吨的泥沙。

3.2 采砂活动真的停止了吗？

洪泽湖非法采砂严重，中央和当地政府出台了各种政策法规，采取了多种措施，严厉打击非法采砂活动。利用2002年以来Landsat系列数据，发现2012年2月第一观测到采砂船，2016年达到最多（月均658艘），2017年3月后禁止采砂船采砂，移除船只的采砂设备，并统一注册并集中停放到韩桥、界集水域。但是，很多采砂船昼伏夜出，是否真正停止偷采，采砂船究竟何时真正停止，政策的效力如何？作者利用VIIRS DNB灯光数据，很好地监测了夜晚非法采砂船只，结果与Landsat观测结果趋势一致，且由于其观测频率高（每天1次），更准确的展示了采砂船的时空分布特征。同时，发展了一种利用白天和夜晚采砂船数据的采砂强度评估方法，可以有效区分采砂活动对水体的扰动强度。最后，科学评估了政府不同阶段政策对打击采砂船的效果，发现2015年8月前，只对白天采砂有所控制，且只能维持两个月；之后到2016年8月，对白天和晚上都有作用，偷采强度显著下降；但直到2017年“利剑2017”行动，非法采砂全面停止。2002年以来MODIS观测结果表明，洪泽湖水体悬浮物浓度明显呈现3个阶段，与采砂分布时间分布一致：2002-2011年偏低，2012-2017年显著升高，2017年3月份后浓度开始下降，已经接近2012年之前水平。由此，说明在2017年3月以后洪泽湖的非法采砂活动已经全面停止。

3.3 利用VIIRS延续MODIS在未来的观测

利用MODIS获取了洪泽湖近16年的悬浮物数据集，但是MODIS传感器因远超设计年限（Terra, 1999-; Aqua, 2002-），面临随时停止运行的可能；因此，2012年NASA新发射了VIIRS传感器，拟接替MODIS工作；但是两者数据的一致性，特别是在内陆湖泊水环境监测的一致性有待评价。利用2012-2015年大量MODIS和VIIRS同步数据（±1 hour），对比了卫星的Rrc反射率的差异，基于之前构建的悬浮物反演算法，评价了它们在水体悬浮物浓度估算的一致性。MODIS和VIIRS的Rrc数据在红光和绿光波段具有很好的一致性（UPD在10%左右）。利用之前构建的洪泽湖的悬浮物浓度区域经验模型，分别应用于MODIS和VIIRS，都具有较高的估算精度(R2>0.70,UPD<15%)；并且，二者在悬浮物浓度估算具有很好的一致性，对于全湖80%的像元差异都在20%以下。这表明VIIRS能够延续MODIS在内陆水体悬浮物浓度的观测，未来可以提供更长时间序列的水体悬浮物卫星产品数据。

洪泽湖日、夜采砂活动

灯光卫星数据监测全国湖泊晚间采砂活动