黄晓翰,娄全胜, 2,王 孟,苏 文,姜广甲, 2, 4*
(1. 国家海洋局南海环境监测中心,广东 广州 510300;2. 自然资源部南沙珊瑚礁生态系统野外科学观测研究站,广东 广州 510300;3. 国家海洋局南海规划与环境研究院,广东 广州 510310;4. 广东省海洋资源与近岸工程重点实验室,广东 广州 510310)
南海是太平洋最大、最深的边缘海,面积约为350×104 km2,海底地形复杂、岛礁密布[1]。南海近岸海域常年存在上升流,有河口、海湾、深海区等各种复杂生态系统并存,并受淡水、咸水、咸淡水、潮滩湿地和沙洲湿地等生态系统的综合影响[2-3],囊括了光学I类和II类水体,光学特性十分复杂[4]。其中,珠江是南海近岸海域最大的输入河流,平均每年有3 260亿m3的径流量及其携带物进入南海,加之南海近岸复杂的水动力条件,影响了该海域水体的水质环境。
透明度是评价水质环境的重要参数之一。水下光场分布影响了水下生物、物理和化学动力变化过程,例如,浮游植物利用水下光照进行光合作用将无机碳转化为有机碳储存于水体中,而溶解有机物在光照条件下进行降解释放无机碳。同时,水下光场分布影响着水体的热传输过程,进而影响了水动力变化[5-6]。透明度是确定光照在水下传输界面和评估水体营养水平的重要因子,在无法进行水下辐照度测量的情况下,透明度是评价水下光场分布的唯一指标[7]。
与透明度的传统调查方式相比,卫星遥感技术可以高频次、大尺度、长时间序列的监测海洋环境动力变化特征。VIIRS(Visible Infrared Imaging Radiometer Suite)发射于2011年10月28日,是新一代的近极地水色卫星传感器,其水色波段与MODIS(Moderate Resolution Imaging Spectroradiometer)具有相似的波段配置,但具有较大幅宽(3 000 km)和较高空间分辨率(750 m)[8-9]。目前,针对VIIRS卫星遥感数据,相关研究主要集中在数据本身的校正方法[9-10]以及将VIIRS卫星遥感数据的水色产品与EOS MODIS数据水色产品进行比对、验证[11-14],而开展海洋生态环境监测的技术研究相对较少[15-16],并且基于VIIRS数据对南海近岸海域透明度评价方面还未见报道。
利用遥感技术评价水体透明度的方法主要包括经验法(empirical method)和半分析方法(semi-analytical method)[7,17-21]。本文以涵盖光学I类水体和II类水体的南海近岸海域为研究区,利用珠江口、大亚湾和粤西海域的现场生物-光学数据,构建适于南海近岸海域水体光学特征的透明度评价经验算法,并应用于2013—2016年每天的VIIRS卫星遥感影像,以揭示该海域透明度时空变化特征,为新一代水色卫星数据产品验证提供案例,并为海洋生态环境遥感监测业务化运行提供方法论支撑。
分别于2014年5月11日—18日和2015年4月27日—5月15日在南海近岸海域布设样点采集表层水样,测量水体的遥感反射光谱、水体透明度、风速风向、水深等现场数据,同时记录水环境状况,共获取85个有效样点(图1)。2015年10月在粤西海域获取13个样点的现场生物-光学数据,用于南海近岸海域透明度遥感评价模型的精度验证(图2a)。同时,于2015年3月—8月在南海近岸海域卫星准同步采集样点(图2b),用于VIIRS卫星遥感影像预处理结果的验证。
图2 (a)实测验证样点分布和(b)VIIRS卫星遥感数据准同步样点分布
采用水面之上测量法测定水体遥感反射光谱[22]。所用光谱仪为美国ASD公司生产的FieldSpec Pro Dual VNIR野外双通道光谱仪,可在350~1 050 nm波段范围内连续测量,测定水体、灰板以及天空光光谱,通过计算得到水体遥感反射比。测量时天气晴朗无云,风速小于3 m/s,测量时间为9:00—15:00。
在避开太阳直射反射、忽略或避开水面泡沫的情况下,水体总光谱信号(Lsw)可近似看作水体离水辐亮度(Lw)和天空光辐亮度(Lsky)之和,即:
Lsw=Lw+rLsky
(1)
式中,r为气-水界面对天空光的反射率,与风速、太阳高度角和观测几何等有关,取值为0.026~0.028[23]。
由于Lw在遥感中容易测定,因此,遥感反射比(Rrs, rs-1)可以定义为水体离水辐亮度与恰好位于水体表面的向下辐照度(Ed(0+))的比值。则:
Rrs=Lw/Ed(0+)
(2)
Ed(0+)可表示为:
Ed(0+)=Lpπ/ρp
(3)
其中,ρp为标准白板的反射率,校正值为30%;Lp为折合到100%标准板的信号。
综上,可得水体遥感反射比:
(4)
根据VIIRS卫星遥感的波段响应,将实测光谱重采样为模拟的VIIRS数据,重采样规则与文献[24]保持一致。
利用塞氏盘深度(Secchi Disk Depth)测量法测量水体透明度(SD)。使用的测量工具为20 cm直径的黑白塞氏盘。测量时,将塞氏盘水平放入待测水体中,直至其到达“可见”与“不可见”的深度临界值,该深度临界值即为测量点的塞氏盘深度[25]。
在NOAA官网(Climate Services and Monitoring Division, CLASS)下载2013—2016年每天VIIRS L1B(Sensor Data Record, SDR)数据产品(空间分辨率为750 m),共计1 460景。在SeaDAS 7.4软件(OBGP, NASA)中进行数据预处理,包括:几何校正、云掩膜处理、水体提取、大气校正等。其中,大气校正采用近红外迭代法(Iterative NIR correction-BFW10)[15,26]。
数据统计、多元线性拟合等数据处理均在Origin 8.5统计软件中完成。
采用观测值与模拟值之间的相对均方根误差(relative root mean square error,rRMSE)来评价估算模型的优劣[27]。
(5)
南海近岸海域透明度在0.3~30.0 m之间变化,均值为6.1 ± 7.7 m,具有较大的空间变异系数(1.3),说明南海近岸海域透明度空间差异较大。该海域40%以上的透明度小于2.0 m(图3),主要分布在珠江口海域(0.8 ± 0.4 m),但变异系数较小(0.5),说明珠江口海域透明度空间差异不明显。大亚湾海域具有较高的透明度,变化范围为1.2~5.5 m,均值为2.4 ± 1.3 m。在离岸较远海域透明度远大于珠江口和大亚湾,均值为13.5 ± 8.7 m,变化范围为3.0~30.0 m。
图3 南海近岸海域透明度直方图统计
2.2.1 透明度对水体遥感反射光谱的影响
透明度与水体中的光学活性物质密切相关[7],包括叶绿素a、悬浮物和有色溶解有机物(Chromphoric Dissolved Organic Matter, CDOM)[28],因此间接影响了水体的离水辐亮度或者遥感反射比。选取透明度梯度为1(0.9)、5(5.3)、10(9.0)、15(15.0)和25(24.0) m对应的遥感反射比以分析不同透明度对VIIRS波段遥感反射光谱的影响。
图4给出透明度对水体遥感反射信号的影响。不同透明度对遥感反射比具有显著影响:总体上,透明度较低(0.9 m)的遥感反射比较高,而透明度较高(24.0 m)的遥感反射比较低(Rrs(λ)<0.01 sr-1)。特别是Rrs(555),随着透明度的增加而逐渐降低(图4a)。为分析透明度对水体遥感反射比光谱形状的影响,将不同透明度对应的遥感反射光谱进行归一化,即每条光谱的每个波段的遥感反射比除以该光谱中遥感反射比的最大值(图4b、c)。发现当透明度≤5.3 m时,在555 nm处出现明显反射峰且遥感反射比最大,表现为412—555 nm波段范围内Rrs(λ)随波长增加逐渐增大,在555—865 nm波段范围内,Rrs(λ)随波长增加逐渐减小(图4b);同时,透明度较高水体的归一化Rrs(λ)高于透明度较低水体。而透明度>5.3 m时,拐点出现在488 nm处,且归一化Rrs(555)随着透明度升高而降低,而透明度为9.0 m时,555 nm出现较小反射峰(图4c)。
图4 透明度对水体遥感反射光谱的影响:(a)原始遥感反射比;(b、c)归一化遥感反射比
2.2.2 透明度的遥感评价模型
通过分析透明度对水体遥感反射光谱的影响,发现随着透明度的增加,遥感反射光谱在488和555 nm处变化敏感(图4)。因此,构建基于VIIRS的波段比值算法以评价南海近岸海域透明度的时空分布特征,即评价(Rrs(488)/Rrs(555))与透明度之间的相关关系。
两者在南海近岸海域具有较好的正相关关系,随着透明度逐渐增大,Rrs(488)/Rrs(555)逐渐增大(变化范围为0.5~3.5,均值1.1 ± 0.6),但当透明度>10.0 m时,两者的相关性有所降低。经过拟合,发现对数形式结果最好(图5)。基于VIIRS波段的南海海域透明度遥感评价模型表达式为:
图5 透明度与Rrs(488)/Rrs(555)的相关关系
SD = 15.1× ln[Rrs(488)/Rrs(555)] + 8.2 (R2= 0.77,N= 85,p< 0.01)
(6)
利用VIIRS传感器波段的比值构建南海近岸海域透明度的遥感评价算法,并应用于预处理后的VIIRS卫星遥感影像,以获取该海域透明度的时空分布特征。
为分析VIIRS大气校正效果,将准同步的现场遥感反射比与预处理后的VIIRS卫星遥感数据进行对比,发现两者形状相近(图6a),且不同波长的遥感反射比的绝对值相差很小(图6b),可用于评价南海海域透明度的动态变化特征。
图6 (a)准同步采样点实测遥感反射光谱与遥感影像反射光谱的对比;(b)特定波长处实测与遥感影像遥感反射比的对比:412、445、488、555、672、746和865 nm
将公式(6)应用于2015年10月在粤西海域的现场数据,以验证南海近岸海域透明遥感评价算法的稳定性(图7)。结果表明,模型计算结果与实测值具有较好的一致性,均匀分布在1∶1线周围,两者具有较好的相关关系(R2= 0.84,rRMSE= 19.3%,N= 13,p< 0.01),说明该模型具有较好的稳定性,可用于南海近岸海域透明度的业务化监测。
图7 实测透明度与估算透明度之间的关系
将该算法应用于2013—2016年南海近岸海域的VIIRS卫星遥感影像,获取该海域的月均产品,从整体上揭示了该海域透明度的时空变化规律。基于不同时相的VIIRS卫星遥感影像获取了南海近岸海域透明度的时空分布状况(图8),发现该海域的透明度在0.0~ 40.0 m范围内变化。整体上看,不同年度5月份透明度的空间变化显著,但离岸较远海域透明度高于离岸较近海域,特别是沿岸水体透明度较低,大部分海域透明度小于12.0 m。受水动力影响,珠江口东部沿岸透明度高于西部。在海南附近海域,不同年度透明度变化较大。
图8 2013—2016年5月南海近岸海域透明度的空间分布(白色代表有云区域)
经验模型受限于训练数据集的选择,并且与水体光学特性密切相关[29]。但遥感反射光谱是对水体及周围环境信息的综合响应,基于水体生物光学模型的半分析算法仍与样点个数、代表性以及水体环境具有较大的相关性,限制了其在不同水体中的推广应用。针对南海近岸海域,样点分布于悬浮泥沙和CDOM占主导地位的珠江口海域[30]、相对封闭且初级生产力较高的大亚湾海域[31]以及属于光学I类水体的远离近岸水体,同时,使用不同年度的样点数据对模型进行了验证(R2= 0.84,rRMSE= 19.3%),样点代表性较强,说明该模型在南海近岸海域有较强的适用性,为海洋生态环境的业务化监测提供了方法论基础。
针对不同透明度大小对应的遥感反射光谱以及VIIRS传感器的波段设置(图4),488和555 nm处遥感反射比的斜率对不同透明度具有较强的区分性。CDOM在紫外和可见光波段吸收强烈,悬浮泥沙在555 nm附近具有较强的散射特征,浮游植物色素在490 nm处吸收强烈,而叶绿素a在665 nm处具有显著的吸收峰。水体中的悬浮泥沙、浮游植物色素和CDOM是影响透明度的主要因素[7],因此,考虑488和555 nm处的遥感反射比与透明度的相关关系具有一定的物理意义。
从图5可以看出,当透明度>10.0 m时,模型敏感度有所降低,变化平缓,但对透明度仍有较好的区分能力(图7)。而仅考虑透明度<10.0 m的样点,透明度与Rrs(488)/Rrs(555)的相关关系有所提高(图9,R2=0.82),模型敏感度较高。将模型应用于VIIRS影像数据,南海近岸海域透明度<40.0 m(图8),且大部分海域透明度<32.0 m,说明该模型的适用性较强,稳定性较高。
图9 透明度(<10 m)与Rrs(488)/Rrs(555)的相关关系
为了进一步评价VIIRS卫星遥感影像产品的优劣,将2013年3月7日VIIRS和MODIS透明度产品进行对比。由于MODIS传感器中心波长为555 nm的波段(545—565 nm)空间分辨率为500 m,而中心波长为488 nm的波段(483—493 nm)空间分辨率为1 000 m,故采用重采样为1 000 m空间分辨率的555 nm波段参与透明度评价模型运算。同时,VIIRS和MODIS Aqua的成像时间相近,因此选择VIIRS与MODIS Aqua卫星遥感影像数据透明产品进行比较。
基于VIIRS和MODIS卫星遥感影像数据的透明度产品均能反映出南海近岸海域透明度的空间分布特征(图10),总体上表现为离岸较近海域透明度较低,反之较高,且均可表征该海域水动力状况。两者成像时间差别较小(VIIRS:5:45 UTC,MODIS Aqua:5:17 UTC),可从空间上进行比较。在珠江口以东海域,两者透明度表现出相似的空间分布特征,但在细节上VIIRS透明度产品刻画的较为细致,主要与其较高的空间分辨率有关。而在珠江口以东海域,特别是湛江海域附近,两者透明度绝对值有一定差异,基于VIIRS计算得到的透明度较高,而MODIS较低,可能与两者传感器的性能[9]以及大气校正参数输入有关[15]。为保证VIIRS传感器在轨校准和数据质量,卫星上搭载了一系列在轨校准器(onboard calibrators, OBCs),并携带漫反射板和稳定器以减少MODIS存在的非均匀光照带来的误差[9]。对于大气校正方法——近红外迭代法(Iterative NIR correction-BFW10)更适于清澈至中等浑浊水体[15,26],并将该方法应用于VIIRS卫星遥感影像,大气校正效果较好。对于MODIS来说,大气校正输入的大气参数与VIIRS获取的大气参数有一定差别,导致两者透明度产品在珠江口以西海域有一定差异。
图10 VIIRS和MODIS透明度产品对比
总之,作为新一代在轨的水色传感器VIIRS,具有较高的性能和空间分辨率,将在未来水色遥感水体生态参量的定量评价中发挥重要作用,并将推动水色遥感的业务化应用。
将南海近岸海域透明度遥感评价模型应用于2013—2016年每月VIIRS卫星遥感影像(2016年4月整个海域云覆盖,无影像),发现南海近岸海域透明空间差异显著(图11),透明度在0.0~40.0 m之间变化,总体表现为离岸较近海域透明度低于离岸较远海域,西部海域低于东部。特别是珠江口海域,透明度均小于8.0 m。南海近岸海域透明度季节变化显著,总体上夏季最高,其次是春季和秋季,冬季最低(图11)。薛宇欢等[32]研究了中国近海海水透明度分布特征和季节变化,发现南海近岸海域水体透明度变化范围为0.2~32.0 m,且冬季最低,夏季最高,与本文研究结果一致。
图11 南海近岸海域透明度空间和季节变化
图12给出南海近岸海域透明度均值年际变化,最大值出现在2015年8月,透明度均值为21.2 m。最小值出现在2016年5月(2.6 m),可能由于该月获取的卫星遥感影像主要集中在离岸较近海域。从不同年度进行比较,2015年透明度均值最大,为13.7 m,其次2013年(12.0 m),2014年和2015年透明度均值相差不大,分别为10.8 m和9.7 m。
图12 南海近岸海域透明度均值年际变化(2013—2016) (其中,不同颜色的横线代表透明度年均值)
(1)南海近岸海域透明度变化范围为0.3~30.0 m(均值为6.1 ± 7.7 m),且该海域40%以上的透明度在2.0 m以下,而离岸较远海域透明度较大,均值为13.5 ± 8.7 m。
(2)利用涵盖光学I类和II类水体的生物-光学数据,参照VIIRS卫星遥感传感器波段设置,构建了南海近岸海域透明度遥感评价模型,该模型基于VIIRS的波段比值(Rrs(488)/Rrs(555))和透明度之间的定量关系,精度较高(R2= 0.77,N= 85,p< 0.01)。该模型在透明度>10.0 m时敏感性有所降低,但对透明度(10.0~30.0 m)仍有较好的区分能力。VIIRS具有较高的空间分辨率和性能,对透明度的细节刻画能力优于MODIS。
(3)将南海近岸海域透明度遥感评价模型应用于预处理后的VIIRS卫星遥感影像,发现南海近岸海域透明度在0.0~40.0 m之间变化,总体表现为离岸较近海域透明度低于离岸较远海域,西部海域低于东部。