黄渤海悬浮颗粒物密度的空间分布及其对固有光学特性的影响*

苏校平，孙德勇，丘仲锋，王胜强，路 颖，吴晨颖，叶之翩，岳小媛

(南京信息工程大学海洋科学学院，江苏南京 210044)

黄渤海悬浮颗粒物密度的空间分布及其对固有光学特性的影响*

苏校平，孙德勇**，丘仲锋，王胜强，路 颖，吴晨颖，叶之翩，岳小媛

(南京信息工程大学海洋科学学院，江苏南京 210044)

悬浮颗粒物密度 固有光学特性 时空分布 黄渤海

0 引言

【研究意义】水色受悬浮颗粒物的影响最为显著,这是我国近海水体光学特性的普遍特点之一[1-2]，悬浮颗粒物在很大程度上控制着海水固有光学特性的变化，因此它在沿海水域中起着十分重要的作用[3]。除此之外，悬浮颗粒物还影响水体的光场变化，进而影响海洋的初级生产力[4]。悬浮颗粒物密度，是描述悬浮颗粒物的一个重要物理参数。了解悬浮颗粒密度的空间分布特征，对悬浮颗粒动力学过程、海洋环境控制因素以及海洋水色遥感机理与应用的深入理解具有重要意义。在悬浮颗粒物含量较高的黄渤海海域，悬浮颗粒物密度对水体生物和化学过程都起着极其重要的作用，其分布以及对固有光学特性的影响机制研究也对海洋水色遥感研究具有重要意义[5]。【前人研究进展】国外学者在悬浮颗粒物与固有光学特性的关系上已经做了大量工作。Neukermans等[6]开展了在光学性质复杂的海域中悬浮颗粒物密度对衰减系数和后向散射系数的影响研究，得到了悬浮颗粒物密度与固有光学特性的理论公式，并将悬浮颗粒物密度与平均粒径结合研究，得到了比衰减系数变化几乎全部受这两个量影响的结论；Boss等[7]研究证实了悬浮颗粒物浓度与固有光学特性之间的密切关系；Babin等[4]对于近岸和大洋水体悬浮颗粒物浓度与光散射特性的关系进行了研究探讨。Bowers等[8]证明矿物质占多数的陆架浅海中散射系数的变化主要受密度影响。同时，国内学者已经开展了关于悬浮颗粒物浓度以及水体固有光学特性的研究。刘炜等[1]利用2003年4月黄东海实验数据对中国黄东海海区的悬浮颗粒物的散射特性进行了研究；李铜基等[2]对中国近海水体的后向散射系数进行了相关研究；虞兰兰等[9]基于2006年6月以及2007年1月实测数据，对中国黄、东海的悬浮颗粒物浓度及其分布进行了研究；孙德勇等[10]主要针对内陆二类水体固有光学特性及其季节性差异进行了研究。【本研究切入点】前人多以实验室测量和理论运算为主研究悬浮颗粒物对固有光学特性的影响，而对于现场测量数据的应用以及悬浮颗粒物密度对固有光学特性的影响研究比较少。【拟解决的关键问题】利用现场测量仪器实测黄渤海水体固有光学特性(包括衰减系数、散射系数、后向散射系数、比衰减系数、比散射系数、比后向散射系数、衰减效率、散射效率和后向散射效率)，以及悬浮颗粒物浓度，同时结合其他相关参量间接计算悬浮颗粒物的密度，得到其空间分布特征并分析悬浮颗粒物密度对固有光学特性的影响。

1 材料与方法

1.1 观测站位

研究区域为黄海和渤海，研究数据集由2014年11月和2015年8月两个航次的测量数据组成。航次站位分布分别如图1所示。

图1 航次调查站位(a:2014年11月航次;b:2015年8月航次)

Fig.1 Locations of the sampling station(a:In November 2014 voyage;b:In August 2015 voyage)

1.2 数据测量及处理

总体积浓度VC、粒子数浓度NC均由LISST-100 X (Type C)现场激光粒度仪测量得到， LISST-100 X (Type C)根据激光的衍射原理，运用Mie散射理论，可区分32个呈对数分布的粒级并可以给出每个粒级的体积浓度[VC]i(i=1～32)。仪器得到的数据为不同粒径粒子的体积浓度，将前32列数据求和得到悬浮颗粒物总体积浓度VC。将悬浮颗粒物近似看做球体，则单位体积内所含的粒子个数由下式得：

V(D)=N′(D)(πD3/6)，

(1)

其中，V(D)为32个不同粒级粒子各自的体积浓度，D为对应的每个粒级粒径大小。

N(D)=N′(D)/ΔD，

(2)

△D为每个粒级的尺寸范围，N(D)为每单位尺寸、单位体积浓度所含的粒子个数(粒子数·m-3·μm-1)。具体原理及处理步骤可见仪器说明[11]。

粒子光束衰减系数cp和散射系数bp均由AC-S水体光吸收-衰减系数测量仪测得。实际操作中需要对AC-S测得的衰减系数和吸收系数的原始数据进行温度、盐度和散射校正，最终得到的水体散射系数bnw是校正后的水体衰减系数cnw与吸收系数anw的差，即

bnw=cnw-anw，

(3)

设纯水的散射系数为bw，则悬浮颗粒物的散射系数bp为

bp=bnw-bw。

(4)

后向散射系数bbp的测量仪器为六通道后向散射仪(HS6)。HS6在412～852 nm之间有6个波段配置，通过测量出140°处的体散射函数U(θ)然后计算得到后向散射系数[12]。HS6使用时的数据校正主要为sigma校正。假设后向散射系数为bbp，没有校正过的是bpu，校正系数为e,则校正方法如下：

bbp=e·bpu。

(5)

TSM是运用LISST-100 X (Type C)结合GF/F滤膜测定得到。通过Neukermans等[6]的研究，得出ρa与TSM和VC的关系见式(6):

(6)

(7)

(8)

(9)

平均衰减效率Qce是由式(10)(11)推导而来，其中，AC为总粒子截面积浓度，是由LISST-100 X(Type C)粒度仪测出32个粒级的截面积浓度后求和得到的总的截面积浓度。

(10)

(11)

其中，粒径D由LISST-100 X (Type C)粒度仪测得。由式(10)(11)推导出式(12)(文献[12])。

(12)

同理，由式(10)(13)推导出式(14)，计算得到平均散射效率Qbe：

(13)

(14)

由式(10)和式(15)推导出式(16)，计算得到平均后向散射效率Qbbe

(15)

(16)

2 结果与分析

2.1 实测数据处理

表1 2014年11月航次光学参量统计结果

Table 1 Statistical parameters of data collected in the Yellow Sea and Bohai Sea in November 2014

参量Parameters最小值Minimum最大值Maximum平均值Mean标准差SD变异系数CV(%)TSM(mg/L)0．04447．685．8999．514161．3ρa(kg/L)0．0060．8750．2220．15569．82VC(μL/L)3．31151．225．5232．28126．5cp(m-1)12．1520．8116．312．44615．00bp(m-1)0．09518．072．2213．741168．4bbp(m-1)0．0010．6630．0570．116203．5c∗p(m2/g)0．402428．119．2653．03275．3b∗p(m2/g)0．006141．54．90419．22391．9b∗bp(m2/g)4E040．1360．0090．017188．9Qce1．503131．547．8538．8981．27Qbe0．01494．908．11515．65192．9Qbbe0．0020．0520．0250．01560．00

表2 2015年8月航次光学参量统计结果

Table 2 Statistical parameters of data collected in the Yellow Sea and Bohai Sea in August 2015

参量Parameters最小值Minimum最大值Maximum平均值Mean标准差SD变异系数CV(%)TSM(mg/L)0．77331．415．5785．692102．1ρa(kg/L)0．0310．9510．2510．17870．88VC(μL/L)7．469118．125．7822．687．64cp(m-1)10．0526．6715．112．62917．4bp(m-1)0．1168．3310．6961．13162．3bbp(m-1)0．0010．3260．0160．045285．3c∗p(m2/g)0．55519．35．064．02279．49b∗p(m2/g)0．0321．230．1570．171108．9b∗bp(m2/g)0．0010．0120．0020．002119．5Qce0．251．9141．0670．45342．47Qbe0．0050．1560．0370．03183．44Qbbe0．0010．0030．0010．001127．1

图2直观分析了两个航次cp、bp、bbp数值分布情况。2014年11月航次cp数值分布以18左右最多(图2a)；数值上bp比较小，bbp最小，但二者分布皆比较均匀(图2b和2c)。2015年8月航次cp数值分布在15左右较多，整体分布比较规律，无值处很少(图2d)；该航次bp仍然比较小但大于bbp，二者分布比较规律，缺值区域较少(图2e和2f)。

同理，两个航次Qce、Qbe、Qbbe的数值分布如图4所示：2014年11月航次Qce分布比较均匀，以小于10为主(图4a)；数值上Qbe较小，1附近较多且分布比较均匀(图4b)；Qbbe数值上最小，分布规律而且均匀(图4c)。2015年8月航次的Qce分布比较均匀，数值以1.5左右为主(图4d)；Qbe比较小，多数分布在0.05以内，缺值区域较多(图4e)；Qbbe数值上最小，多数分布在0.000 16左右，缺值区域较少，分布较均匀(图4f)。

为了直观分析比较悬浮颗粒物密度变化，图5给出了两个航次密度的频数分布，由图5可见，两个航次ρa整体分布都比较均匀且分布趋势相似，缺值区域较少。两个航次的数值分布区间大致相同，2014年11月航次在0.2附近的样本数比较密集，而2015年8月航次中落在0.1左右的样本数较多。

图2 2014年11月(a～c)和2015年8月(d～f)航次cp、bp、bbp的频数分布

图3 2014年11月(a～c)和2015年8月(d～f)航次、、的频数分布

图4 2014年11月(a～c) 和2015年8月(d～f)航次Qce、Qbe、Qbbe的频数分布

图5 2014年11月(a)和2015年8月(b)航次ρa的频数分布

Fig.5 Frequency distribution of ρain November 2014(a) and August 2015(b)

2.2 悬浮颗粒物密度的空间分布

对于2014年11月航次，表层在渤海大部和黄海北部以及中部悬浮颗粒物密度比较小，黄河口及山东半岛附近增大，黄海南部密度从中间向两侧降低，两侧出现极低值(图6a)；中层悬浮颗粒物密度除渤海和黄海中部有小片高值外，大部分海域密度比表层要低，最大值出现在黄海中部(图6b)；底层在渤海东部、大清河附近以及山东半岛南部密度较大，其他海域较小，底层最大值出现在黄海中部，围绕最大值向四周密度降低，至黄海南部最低(图6c)。对于2015年8月航次，表层极大值出现在大清河附近，除了黄河口附近，整个渤海密度都比较大，黄海密度较低，以黄海南部为中心向两侧密度增加(图6d)；中层渤海密度明显小于表层，大清河附近海域以及渤海东北海域密度较大，整个中层密度最大值出现在黄海东北部和东南部小片海域，山东半岛附近及黄海中部偏西部分海域密度较大，其他海域密度均较小(图6e)，由于数据缺失，只能得到底层黄海中部和南部海域密度分布；底层整体趋势为以黄海中部的极大值为中心向四周密度降低，盐城东侧出现极小值(图6f)。综上可知，两个航次悬浮颗粒物密度的空间分布比较复杂，但整体趋势显示中层密度要小于表层和底层。在各层密度分布图中，黄河口附近均有密度偏小值区域，山东半岛附近有密度偏大值区域，随着深度增加该大值区域有南移并增大的趋势。

图6 2014年11月航次(a～c)和2015年8月航次(d～f)表、中、底3个层次密度分布

Fig.6 Density distribution of surface,middle and bottom in November 2014 (a～c) and August 2015 (d～f)

2.3 密度对固有光学特性影响分析

2.3.1 对cp、bp、bbp的影响

分别利用两个航次数据做回归分析，结果如图7所示，各量之间回归方程及决定系数见表3。2014年11月航次有7%的cp受ρa影响，而2015年8月航次中有1%的cp受ρa影响，两年的汇总量中，有约2%的cp受ρa影响(图7a)。2014年11月航次有约6%的bp受ρa影响，2015年8月航次有约4%的bp受ρa影响，两年的汇总量中，有约2%的bp受ρa影响(图7b)。2014年11月航次，有11%的bbp受ρa影响，2015年8月航次中有约7%的bbp受ρa影响，两年的汇总量中，有约7%的bbp受ρa影响(图7c)。由此可知，整体上两个航次的悬浮颗粒物密度与cp、bp、bbp的相关性都较小。

图7 两个航次以及两个航次汇总的ρa与 cp、bp、bbp拟合

Fig.7 Density of two voyages and the summary of two voyages fitted to cp,bpand bbp

表3ρa分别与cp、bp和bbp的回归方程及决定系数

Table 3 The regression equation and correlation factor ofρawithcp,bpandbbprespectively

所在图表Figure数据来源Dataset方程类型Equationtype回归方程Regression决定系数R2DeterminationcoefficientR2P值PvalueFig．7a2015年8月航次InAugust2015voyage指数Exponentialy=15．18e-0．07x0．010．66>0．05Fig．7a2014年11月航次InNovember2014voyage对数Logarithmicy=0．97ln(x)+17．930．070．12>0．05Fig．7a两个航次汇总Inbothvoyages对数Logarithmicy=0．46ln(x)+16．420．020．13>0．05Fig．7b2015年8月航次InAugust2015voyage指数Exponentialy=0．34e0．91x0．040．42>0．05Fig．7b2014年11月航次InNovember2014voyage乘幂Powery=0．35x-0．410．060．16>0．05Fig．7b两个航次汇总Inbothvoyages对数Logarithmicy=-0．49ln(x)+0．610．020．23>0．05Fig．7c2015年8月航次InAugust2015voyage乘幂Powery=0．01x0．430．110．48>0．05Fig．7c2014年11月航次InNovember2014voyage乘幂Powery=0．04x0．680．070．11>0．05Fig．7c两个航次汇总Inbothvoyages乘幂Powery=0．02x0．510．070．19>0．05

2.3.3 对Qce，Qbe,Qbbe的影响

利用得到的数据做散点图并做回归分析，结果如图9所示，各量之间回归方程及决定系数见表5。2014年11月航次有6%的Qce变化受ρa影响，而2015年8月航次中只有4%，两年的汇总量中有6%的Qce变化受ρa影响(图9a)； 2014年11月航次中有10%的Qbe变化受ρa影响，而2015年8月航次有19%，两年的汇总量中有10%的Qbe变化受ρa影响(图9b)； 2014年11月航次有16%的Qbbe变化受ρa影响，而2015年8月航次中有24%，两年的汇总量中有1%的Qbbe变化受ρa影响(图9c)。综上所述，ρa与Qbbe的相关性比其他两个量大很多，2014年11月航次ρa与Qbbe的相关性要比与Qce和Qbe两个量大，但是仍然不及ρa与比固有光学量的相关性。

3 讨论

渤海的大清河和滦河附近密度值比较大，原因可能是滦河向渤海注入的大量泥沙产生堆积，使得近岸水深较浅，较浅的海水受风、波浪和流的影响使得水体底层运动比较活跃，易形成再悬浮，造成滦河附近底层密度高于表层。黄河口附近的密度偏小区域可能是由于黄河河口属于弱潮陆相河口，其突出特点有二：一是水少沙多、沙粗、洪枯悬殊、洪峰陡涨陡落，二是潮差相对小、海潮流速也小、海洋动力弱等[13]。由于较弱的海洋动力以及较粗的泥沙导致再悬浮不活跃，因此该区域悬浮颗粒物密度较小。山东半岛附近随着深度增加有南移趋势的密度大值区域可能是由黄海沿岸流引起的。黄海沿岸流水色浑浊，带有大量悬浮颗粒物而且具有冬强夏弱的特点。本研究采用的两个航次数据均来自秋冬季节，北风占主导地位，向南流动的黄海沿岸流增强，因此呈现上述现象及趋势。

所在图表Figure数据来源Dataset方程类型Equationtype回归方程Regression决定系数R2DeterminationcoefficientR2P值PvalueFig．8a2015年8月航次InAugust2015voyage指数Exponentialy=8．08e-3．05x0．515．53E10<0．05Fig．8a2014年11月航次InNovember2014voyage乘幂Powery=0．96x-1．070．510．005<0．05Fig．8a两个航次汇总Inbothvoyages乘幂Powery=0．97x-0．940．500．001<0．05Fig．8b2015年8月航次InAugust2015voyage指数Exponentialy=0．18e-1．98x0．210．007<0．05Fig．8b2014年11月航次InNovember2014voyage乘幂Powery=0．02x-1．460．330．018<0．05Fig．8b两个航次汇总Inbothvoyages乘幂Powery=0．03x-1．080．260．019<0．05Fig．8c2015年8月航次InAugust2015voyage指数Exponentialy=0．01e-1．27x0．070．396>0．05Fig．8c2014年11月航次InNovember2014voyage乘幂Powery=0．01x-0．380．100．052>0．05Fig．8c两个航次汇总Inbothvoyages乘幂Powery=0．01x-0．420．090．031<0．05

图9 两个航次以及两个航次汇总的ρa与Qce、Qbe、Qbbe散点图

Fig.9 Density of two voyages and the summary of two voyages fitted to Qce,Qbeand Qbbe

表5ρa分别与Qce、Qbe及Qbbe的回归方程及决定系数

Table 5 The regression equation and correlation factor ofρawithQce,QbeandQbberespectively

所在图表Figure数据来源Dataset方程类型Equationtype回归方程Regression决定系数R2DeterminationcoefficientR2P值PvalueFig．9a2015年8月航次InAugust2015voyage乘幂Powery=1．18x0．130．040．152>0．05Fig．9a2014年11月航次InNovember2014voyage指数Exponentialy=42．59e-1．94x0．060．198>0．05Fig．9a两个航次汇总Inbothvoyages乘幂Powery=15．05x-0．340．060．116>0．05Fig．9b2015年8月航次InAugust2015voyage乘幂Powery=0．55x0．410．190．045<0．05Fig．9b2014年11月航次InNovember2014voyage乘幂Powery=0．34x-0．730．100．219>0．05Fig．9b两个航次汇总Inbothvoyages乘幂Powery=0．08x-0．470．100．160>0．05Fig．9c2015年8月航次InAugust2015voyage乘幂Powery=0．01x0．570．240．009<0．05Fig．9c2014年11月航次InNovember2014voyage乘幂Powery=0．04x0．360．160．003<0．05Fig．9c两个航次汇总Inbothvoyages对数Logarithmicy=0．01ln(x)+0．020．010．383>0．05

4 结论

在研究黄渤海的悬浮颗粒物密度对固有光学特性的影响时，本文所研究的散射系数和后向散射系数均选取波段640 nm作为参考波长，并未对其它波段进行研究，下一步可对不同波长处的影响进行分析，量化其光谱差异性。同时，在空间分布上，本研究仅选取两个航次数据，接下来可选取长时间序列数据对空间分布进行定性研究。

致谢

感谢中国海洋大学信息科学与工程学院张亭禄老师的支持，感谢陈树果博士在野外数据观测过程中给予的帮助。

[1] 刘炜,李铜基,朱建华,等.黄东海海区总悬浮物散射特性研究[J].海洋技术,2007,26(2):42-46. LIU Y,LI T J,ZHU J H,et al.Study of scattering properties of total suspended matter in the Yellow Sea and East China Sea[J].Ocean Technology,2007,26(2):42-46.

[2] 李铜基,陈清莲,杨安安,等.黄东海春季水体后向散射系数的经验模型研究[J].海洋技术,2004,23(3):10-14. LI T J,CHEN Q L,YANG A A,et al.Study of the spring empirical model for backscattering coefficient in Yellow Sea and East China Sea[J].Ocean Technology,2004,23(3):10-14.

[3] ASTORECA R,DOXARAN D,RUDDICK K,et al.Influence of suspended particle concentration,composition and size on the variability of inherent optical properties of the Southern North Sea[J].Continental Shelf Research,2012,35:117-128.

[4] BABIN M,MOREL A,FOURNIER-SICRE V,et al. Light scattering properties of marine particles in coastal and open ocean waters as related to the particle mass concentration[J].Limnology and Oceanography,2003,48(2):843-859.

[5] 陈亚慧,丘仲锋,孙德勇,等.黄渤海悬浮颗粒物粒径的遥感反演研究[J].光学学报,2015,35(9):50-59. CHEN Y H,QIU Z F,SUN D Y,et al.Remote sensing of suspended particle size in Yellow Sea and Bohai Sea[J].Acta Optica Sinica,2015,35(9):50-59.

[6] NEUKERMANS G,LOISEL H,MÉRIAUX X,et al.In situ variability of mass-specific beam attenuation and backscattering of marine particles with respect to particle size,density,and composition[J].Limnology and Oceanography,2012,57(1):124-144.

[7] BOSS E,TAYLOR L,GILBERT S,et al.Comparison of inherent optical properties as a surrogate for particulate matter concentration in Coastal waters[J].Limnology and Oceanography,2009,7(11):803-810.

[8] BOWERS D G,BRAITHWAITE K M,NIMMOSMIT- H W A M,et al.Light scattering by particles suspended in the sea:The role of particle size and density[J].Continental Shelf Research,2009,29(14):1748-1755.

[9] 虞兰兰,江文胜.黄、东海悬浮细颗粒物浓度和粒径分布变化研究[J].海洋与湖沼,2011,42(4):474-481. YU L L,JIANG W S.Seasonal variations in the distributions of suspended fine particulate matter in the Yellow Sea and the East China Sea[J].Oceanologia et Limnologia Sinica,2011,42(4):474-481.

[10] 孙德勇,李云梅,王桥,等.内陆湖泊水体固有光学特性的典型季节差异[J].应用生态学报,2008,19(5):1117-1124. SUN D Y,LI Y M,WANG Q,et al.Differences of inherent optical properties of inland lake water body in typical seasons[J].Chinese Journal of Applied Ecology,2008,19(5):1117-1124.

[11] 中华人民共和国国家质量监督检验检疫总局,中国国家标准化管理委员会.海洋监测规范 第4部分:海水分析:GB 17378.4-2007[S].北京:中国标准出版社,2008. General Administration of Quality Supervision,Inspection and Quarantine of the People’s Republic of China,Standardization Administration of the People’s Republic of China.The Specification for Marine Monitoring-Part 4:Seawater Analysis:GB 17378.4-2007[S].Beijing:China Standards Press,2008.

[12] BOSS E,PEGAU W S,LEE M,et al.Particulate backscattering ratio at LEO 15 and its use to study particle composition and distribution[J].Journal of Geophysical Research:Oceans,2004,109(C1):C01014.

[13] 梁若皓,樊彦国.黄河口演变分析研究[J].灌溉排水学报,2007,26(S1):224-225. LIANG R H,FAN Y G.Analysis on the evolution of Yellow River Estuary[J].Journal of Irrigation and Drainage,2007,26(S1):224-225.

(责任编辑：陆 雁)

Spatial Distribution of Suspended Particle Density in the Yellow Sea and Bohai Sea and Its Effect on Inherent Optical Properties

SU Xiaoping，SUN Deyong，QIU Zhongfeng，WANG Shengqiang，LU Ying, WU Chenying，YE Zhipian，YUE Xiaoyuan

(School of Marine Science，Nanjing University of Information Science &Technology，Nanjing，Jiangsu，210044，China)

【Objective】To study the spatial distribution characteristic of suspended particle density in the Yellow Sea and Bohai Sea and its impact on the inherent optical properties (IOPs) of the water bodies.【Methods】The in situ data in November 2014 and August 2015 were collected to get the suspended particle density and draw the map of spatial distribution.Fitting analysis was conducted between the data of suspended particle density and those of inherent optical properties.【Results】In the surface,middle and bottom layers,the suspended particle density in the areas around the Yellow River Estuary was relatively low,while that around the Shandong Peninsula was relatively high.And the areas with high suspended particle density tended to move southwardly with the depth increasing.The P value of the fitting between the suspended particle density and the mass-specific inherent optical properties was less than 0.05,but the other IOPs were more than 0.05.【Conclusion】The effect of the suspended particle density is greater than 51% of the variability on the mass-specific beam attenuation coefficient,and is small on the other IOPs.

suspended particle density，inherent optical properties，spatial distribution，the Yellow Sea and Bohai Sea

2016-08-04

苏校平(1993-)，女，硕士研究生，主要从事海洋水色遥感研究。

*国家自然科学基金项目(41276186，41576172，41506200)，江苏省自然科学基金项目(BK20151526，BK20150914)，江苏省高校自然科学基金项目(15KJB170015)，“全球变化与海气相互作用”专项(GASI03030101)，江苏省“青蓝工程”优秀青年骨干教师项目和南京信息工程大学大学生实践创新训练计划项目(201510300073，201610300035，201610300064)资助。

http://www.cnki.net/kcms/detail/45.1206.G3.20170106.1039.003.html

P73

A

1005-9164(2016)06-0538-10

**通信作者：孙德勇(1984-)男，副教授，主要从事水光学、水环境遥感等方面的研究，E-mail：sundeyong1984@163.com。

广西科学Guangxi Sciences 2016,23(6):538～547

网络优先数字出版时间：2017-01-06 【DOI】10.13656/j.cnki.gxkx.20170106.001