膠州灣及青島近海水中顆粒物散射特性及其分類

劉曉燕，楊 倩*，陳樹果，胡連波，宗芳伊，劉巧君

1. 齊魯工業大學(山東省科學院)，山東省科學院海洋儀器儀表研究所，山東 青島 2660002. 中國海洋大學信息科學與工程學院，山東 青島 266000

引 言

水中粒子的散射系數和后向散射系數是衡量粒子屬性的兩個重要參數，利用水體的散射特性可以研究海水粒子的組成成分，這是近岸海洋光學研究的一個重要問題。研究表明粒子散射和后向散射系數的光譜模型都服從冪函數[1]，宋慶君等[2]利用中國黃、東海的現場測量數據，得到該海區的水體總后向散射系數bb與水體總散射系數b的關系為乘冪關系。Van de Hulst[3]的研究顯示粒子的散射特性取決于粒子的大小、折射率、組成及形狀。因此，利用多個測量的光學參數可以推斷懸浮粒子的組成成分[4],海中懸浮粒子后向散射比bbp/bp可用于繪制近岸區域不同類型粒子的分布圖[5]。Zhang等[6]利用中國黃、東海的夏季現場光學實驗數據，計算粒子折射率信息分析了該海域的粒子組成。該研究顯示黃東海的粒子后向散射比光譜，即粒子后向系數與散射系數的比值bbp(λ)/bp(λ)，在442～676 nm不隨波長變化。McKee等[7]利用英國沿岸水體的現場測量數據，通過回歸分析表明，對于礦物粒子主導的水體，粒子后向散射比具有光譜依賴性的。

本研究利用2014年8月和2015年10月在膠州灣及青島近海進行的兩次海洋光學現場實驗數據，包括海水總吸收系數a和衰減系數c，以及粒子后向散射系數bbp。膠州灣是伸入內陸的半封閉性海灣，屬于典型的二類渾濁水體,光學性質較復雜,而且其富營養化問題突出,對其水質進行實時有效監測意義重大。分析了膠州灣及青島近海的水中顆粒物的散射特性，并應用其散射特性研究了該海域的粒子組成及空間分布情況，最后研究了散射特性與懸浮顆粒物(SPM)濃度之間的關系。這是改進膠州灣及青島近海水色衛星數據生物光學反演算法、實現衛星遙感監測水質的重要基礎。

1 數 據

圖1是2014年8月和2015年10月在膠州灣及青島近海進行現場光學實驗的站點位置圖。其中在膠州灣內設有9個站位(JZ1—JZ9)，在青島近海設有7個站位(QD1—QD7)。其中，2015年航次不包含JZ8站位。

圖1 膠州灣及青島近海實驗站位圖Fig.1 The map of Jiaozhou Bay and Qingdao coastalexperiment stations

兩次實驗使用高光譜吸收衰減儀ac-s測量吸收系數剖面a(λ,z)和衰減系數c(λ,z)剖面，波長范圍為400～760 nm,共85個波段。實驗前對ac-s儀器進行空氣和純水校準。ac-s直接獲得水體吸收系數a(λ)和衰減系數c(λ),為了得到更為精確的吸收系數和衰減系數,需對其作溫度、鹽度校正和散射校正[8-9]

amts-nw(λ)=am-nw(λ)-[Ψt(t-tr)+Ψs(s-sr)]

(1)

anw(λ)=amts-nw(λ)-amts-nw(λref)

(2)

式(1)是溫鹽校正，am-nw是測量的水體總吸收系數減去純水的吸收值，Ψt和Ψs是溫度和鹽度的校正系數，t和s分別是現場實測溫度和鹽度，tr和sr分別是校正的參考溫度和鹽度，amts-nw是溫鹽校正后的吸收系數； 溫鹽校正后做散射校正[式(2)],λref是散射校正的參考波長，選擇750 nm，anw是經溫鹽和散射校正后的吸收系數。對衰減系數數據做同樣的校正處理，進而可以得到顆粒物散射系數bp(λ)=cnw(λ)-anw(λ)。

2014年8月航次實驗應用后向散射儀BB9測得了隨水深變化的后向散射系數bbp(λ,z)剖面數據，波長為412,440,488,510,532,595,660,676和715 nm。使用Scale factor(由Wetlabs 公司提供) 校正由散射引起的衰減部分,得到117°的水體總體散射函數β(λ)，然后對體積散射函數β(λ)進行吸收校正[10]

β(λ)=βraw(λ)exp(0.031a(λ))

(3)

其中βraw(λ)和β(λ)分別是吸收校正前、后的總體散射函數; 是吸收糾正前的總體散射函數;a是相應波段的校正后的吸收系數，由ac-s獲得。由校正后得到的總體散射函數,減掉純水的體散射函數,可獲得顆粒物的體散射函數βb(λ)，根據下列公式可得到顆粒物后向散射系數bbp(λ)

bbp(λ)=2πXβb(λ)

(4)

其中X=1.1。

2015年10月航次實驗應用HS6測得了隨水深變化的剖面后向散射系數bbp(λ,z)數據，波長為410,442,488,532,550和640 nm。HS6測量的是水體140°方向上的體散射函數，可以計算得到后向散射系數。但由于水體的衰減作用，使得后向散射系數被低估，因此必須對測量的數據進行校正，即SIGMA校正?？衫胊c-s的吸收系數與散射系數數據進行SIGMA校正,校正過程如式(5)[11]

bbp=2πxbb[β(140)-βw(140)]

β=σ(Kbb)βu

σ(Kbb)=k1exp(kexpKbb)

Kbb=a+0.4b

(5)

其中： 2πxbb=6.79；r=4.32；λ0=525 nm；βw(λ0)=0.000 167；k1和kexp是與儀器相關的常數，可在定標文件中找到；βu是測量未校正的總體散射函數。

本研究中，計算各站位表層平均散射系數bp(λ)和后向散射系數bbp(λ)數據時，根據水深和剖面特性取深度前1～5 m內的平均值。

2 粒子散射特性

圖2是2014年和2015年在膠州灣及青島近海測量的(a)散射系數、(b)后向散射系數的光譜曲線圖，其中紅色表示膠州灣內站位，藍色表示青島近海站位?？紤]到BB9后向散射儀在較短波段(<532 nm)的校正問題，圖2僅給出2015年HS6測量的后向散射系數光譜曲線。散射系數bp值隨波長的增大而減小，且膠州灣內站位的散射系數bp光譜斜率(m=0.56)較青島近海站位大(m=0.44)?？偟膩碚f，膠州灣內粒子的散射系數bp和后向散射系數bbp值普遍高于青島近海，且膠州灣內越靠近岸邊的站位，其值越大。

Morel[1]的懸浮粒子散射系數光譜模型bp(λ)被廣泛應用如式(6)

bp(λ)=bp(λ0)[λ/λ0]-m

(6)

利用膠州灣及青島近?，F場站位的bp數據擬合式(6)中的系數m，得出m的范圍在0.21～0.81之間，m的均值為0.50。

懸浮粒子后向散射系數光譜模型bbp(λ)與bp(λ)類似[1]，為

bbp(λ)=bbp(λ0)[λ/λ0]-n

(7)

利用膠州灣及青島近?，F場站位532～715 nm的bbp數據擬合式(7)中的系數n，得出n的范圍在0.36～0.95之間，n的均值為0.70。

圖2 2014年和2015年現場實驗的散射系數(a)與2015年現場實驗的后向散射系數(b)光譜曲線圖

海洋水色遙感中常用的是粒子后向散射系數模型，但現場數據測量多以測量水體的散射系數為主(可由吸收系數和衰減系數計算得到)，如果知道了后向散射系數與散射系數之間的關系，那么將增加大量的后向散射系數數據。利用2014年和2015年的膠州灣及青島近海的現場光學實驗數據，對水體后向散射系數和散射系數之間的關系進行了研究與分析。如圖3所示，以532 nm為例，對后向散射系數與散射系數進行了線性回歸(黑色曲線)和冪指數回歸(紅色曲線)兩種回歸形式。結果發現兩種回歸結果的相關性都很好，冪指數回歸形式略優于線性回歸形式。值得注意的是該bbp導出bp的關系模型[式(8)]僅適應于膠州灣及青島近海海域，對于具有不同顆粒物特性的水體不具有普適性。

bbp(532)=0.015bp(532)1.173

(8)

圖3 532 nm處后向散射系數與散射系數的關系模型

3 粒子折射率及分類

利用水體的散射特性可以分析水中粒子的組成。利用測量的粒子散射系數和后向散射系數，基于Twardowski[5]定義粒子折射率的公式[式(9)]研究了膠州灣及青島近海的粒子組成。

2.295 0r2+2.311 3r4]

(9)

式(9)中的r是粒子衰減系數的光譜斜率，可以用ac-s測量的去除純水貢獻的衰減系數光譜cnw(λ)擬合獲得每個站位的r；bbp可以從BB9和HS6的測量得到。代入式(9)，計算每個站位的粒子折射率np，結果如圖4所示?？梢钥闯?，膠州灣及青島近海的粒子折射率np值分布在1.097～1.197范圍內。說明該海域水中顆粒物主要以無機礦物粒子為主。2014年夏季和2015年秋季，灣口站位(JZ1，QD1和QD7)的折射率變化不大。灣內北部站位(JZ5，JZ6，JZ7和JZ9)的折

圖4 膠州灣及青島近海各站位點的粒子折射率np橫坐標表示站位，從左往右站位點順序依次為JZ1—9，QD1—7

射率相對較低，相比其他站位含有比例稍高的有機粒子。如果膠州灣內以JZ5，JZ6，JZ7站位為斷面，以南水域的np值高于以北水域。這可能與環灣陸地上人類的活動狀態有關。

圖5是依次選取JZ9，JZ6，JZ3，JZ1，QD7，QD6和QD5站位為一不規則斷面(見圖1黑色實線)，計算其粒子折射率np隨站位及剖面的分布圖?？梢钥闯鲈谀z州灣灣口站位JZ1和QD7處，粒子折射率np大于灣內與灣外的站位。灣內從近岸(JZ9)到灣口(JZ1)，粒子折射率逐漸增大，灣內粒子無機成分更多； 青島近海離岸越遠(QD7—5)，粒子折射率逐漸減小，有機粒子的成分增加。JZ6站位為靠近膠州灣大橋的站位點，粒子折射率明顯小于大橋兩邊的站位。2014年，QD6和QD5站位在7 m以下出現無機粒子再懸浮增大的趨勢。

圖5 2014年和2015年粒子折射率剖面分布圖Fig.5 The vertical distribution of particlerefractive index in 2014 and 2015

4 懸浮物后向散射系數與濃度的關系

基于散射特性的粒子組成分霄表明膠州灣及青島近海的粒子組成主要是無機顆粒物。因此，利用2015年海洋光學現場實驗中測量的懸浮物濃度SPM和后向散射系數數據，研究了后向散射系數與SPM濃度的關系，如圖6所示，為基于海洋水色遙感數據的SPM濃度反演工作提供依據?？梢钥闯?，532 nm后向散射系數與懸浮物濃度SPM之間的線性擬合結果與冪指數擬合結果的相關性都較好(相關系數分別為0.74和0.85)，但冪指數擬合的結果優于線性擬合結果。因此，在膠州灣及青島近海海域bbp(532)與懸浮物濃度SPM的冪指數關系為

bbp(532)=0.012[SPM]1.229

(10)

圖6 532 nm的后向散射系數與SPM濃度的關系

5 結 論

利用2014年8月和2015年10月在膠州灣及青島近海進行的兩次海洋光學現場實驗測量的海水吸收系數a和衰減系數c，以及后向散射系數bb，分析了該海域懸浮粒子的散射特性。膠州灣內的散射系數和后向散射系數值普遍高于青島近海； 膠州灣及青島近海海域的粒子后向散射系數和散射系數之間存在冪指數的回歸關系。利用粒子折射率信息分析了該海域的粒子組成及其空間分布情況，該海域的粒子折射率在1.097～1.197之間，粒子組成以無機礦物粒子為主。分析了后向散射系數與懸浮物濃度SPM的關系，冪指數回歸結果優于線性回歸結果，但是SPM的實測數據量較少，后續工作將增加現場實驗數據繼續進行研究。