宮麗嬌，黃 玨*，吳 銘

(1.山東科技大學(xué) 測繪與空間信息學(xué)院，山東 青島 266590；2.浙江大學(xué) 海洋學(xué)院，浙江 舟山 316021)

0 引言

水體固有光學(xué)特性(Inherent Optical Properties，IOPs)僅受水體自身組分影響，不隨光照條件變化而變化，是生物光學(xué)模型的重要輸入?yún)?shù)，對水色遙感的研究具有重要意義。水體總吸收系數(shù)(a)是水體固有光學(xué)特性的組成成分之一，由浮游植物吸收系數(shù)(aph)、有色溶解和碎屑顆粒物吸收系數(shù)(adg)以及純水吸收系數(shù)(aw)組成。了解總吸收系數(shù)的各組成成分，有助于更好地認識生態(tài)系統(tǒng)與生物地球化學(xué)相互作用的問題，值得深入研究[1]。aph是浮游植物生物量的指標，可納入初級生產(chǎn)模型[2-3]；由于光譜相似，非浮游植物顆粒和有色溶解有機物的吸收系數(shù)通常相加并記為有色溶解和碎屑顆粒物吸收系數(shù)。有色溶解和碎屑顆粒物光化學(xué)活性顯著，在地球化學(xué)和碳循環(huán)以及海洋中的其他過程扮演著重要的角色[4-6]。

國內(nèi)外學(xué)者對不同地區(qū)水體光學(xué)特性反演模型進行了諸多的研究，其中，半分析模型是應(yīng)用最廣泛的模型[7-10]。Hoge和Lyon[11-13]開發(fā)了線性矩陣求逆模型(The Linear Matrix Inversion Algorithm，LMI)，后經(jīng)Boss和Roesler[14]修改；Garver和Siege[15]開發(fā)了半分析海洋水色模型(Garver-Siegel-Maritorena，GSM)，Maritorena等[16]開發(fā)了模擬退火程序，改善GSM模型性能；Lee等[6]基于輻射傳輸理論，開發(fā)了Quasi-Analytical Algorithm(QAA)來反演水體固有光學(xué)特性，并通過解析函數(shù)評估模型反演固有光學(xué)特性的可行性[17]。3種反演模型得到了廣泛的驗證與應(yīng)用[1,5,15-16,18-23]。

隨著全球氣候變暖，北極生態(tài)系統(tǒng)和生物地球化學(xué)循環(huán)正在不斷變化[24-25]，研究北冰洋IOPs的需求日益增加。國際海洋水色協(xié)調(diào)小組[15]在2006年利用實測和模擬光學(xué)數(shù)據(jù)集對各種反演模型的性能進行了比較，但是目前的反演模型通常是基于低緯度的海洋數(shù)據(jù)集開發(fā)出來的，北極地區(qū)氣候嚴寒，北冰洋常年被冰雪覆蓋，僅在6—10月部分海冰融化可通航，惡劣的地理條件使得北冰洋的水色遙感研究困難重重。現(xiàn)有研究證明在大陸架區(qū)域內(nèi)的北冰洋，生物光學(xué)特性通常很復(fù)雜，與低緯度地區(qū)觀察到的關(guān)系大不相同[26-32]。

與楚科奇海和波弗特海的實測數(shù)據(jù)相比，Wang等[27]的研究表明GSM模型反演的adg(443)誤差較大；Matsuoka等[30]根據(jù)北冰洋的IOPs優(yōu)化GSM模型，基于實測數(shù)據(jù)的波譜形狀，將數(shù)據(jù)集分為2類，優(yōu)化模型反演的adg(443)的誤差分別為0.35和0.50，但是基于有限數(shù)據(jù)分析得到的結(jié)論具有偶然性。Zheng等[31]發(fā)現(xiàn)，QAA_V5模型在低緯度地區(qū)和北極地區(qū)，均低估了浮游植物的吸收，特別是在北極地區(qū)，模型反演的aph出現(xiàn)大量負值。最新形式的QAA_V6模型，根據(jù)遙感反射率Rrs(670)將水體分為2類并選用不同的參考波段進行計算，可提高模型在復(fù)雜水體中的性能，但目前還未應(yīng)用于北冰洋。

本文基于北冰洋3次巡航取樣的數(shù)據(jù)集，比較LMI、GSM和QAA_V6模型在高緯度水域反演IOPs的性能，并對誤差原因進行分析，結(jié)合模型特性，進行區(qū)域參數(shù)優(yōu)化以適用于北冰洋，為北冰洋光學(xué)遙感提供更全面的參考。

1 數(shù)據(jù)與方法

1.1 研究區(qū)域

北冰洋是以北極點為中心的水域，其海岸線十分曲折，形成了許多淺而寬的邊緣海，其中，波弗特海面積為47.60萬平方千米，楚科奇海為59.50萬平方千米，均在70°N以上，氣溫低，海面幾乎全年被冰雪覆蓋。楚科奇海通過白令海峽與營養(yǎng)豐富的太平洋連通，具有較高的生產(chǎn)力；波弗特海受麥肯齊河影響，其養(yǎng)分利用率低，但是在阿蒙森灣入口處形成的冰穴延長了浮游植物的生長季節(jié)，因而在波弗特海東部觀察到相對較高的初級生產(chǎn)力[11]。

2009年8月，MALINA探險隊在波弗特海巡航取樣，2010年6月和2011年7月，作為美國國家航空航天局項目“ICESCAPE”的一部分，在楚科奇海進行巡航取樣。3次采樣共83個站點，除去幾個觀測異常的站點，本文使用了77個站點的數(shù)據(jù)，站位分布如圖1所示。使用的數(shù)據(jù)集同時測定了水體固有光學(xué)特性和表觀光學(xué)特性，是迄今為止在北冰洋進行的較為全面的光學(xué)測量結(jié)果，這對全面評估模型的性能尤為重要。

圖1 北極地區(qū)測量點位分布Fig.1 Distribution of measurement points in the arctic

1.2 數(shù)據(jù)獲取

1.2.1 遙感反射率測量

Rrs(λ)被定義為上行輻照度Lu(λ，z= 0+)≡Lw(λ)和下行輻照度Ed(λ，z= 0+)≡Es(λ)的比值，其中，Lw(λ)為離水輻亮度，Lw(λ)為水面的入射輻照度。在ICESCAPE航行中，使用剖面反射率輻射儀(Biosperical)測量了包括412，443，490，510，555和665 nm在內(nèi)的18個光譜波段的垂直輻射剖面。在MALINA航行中，使用便攜式水體光學(xué)剖面測量系統(tǒng)C-OPS(Biosperical Instruments Inc.)測量了包括412，443，490，510，555和670 nm在內(nèi)的18個光譜波段。開闊海域在波長為555 nm或者更長處有峰值，而沿海水域則在412～490 nm處有峰值，在這2種數(shù)據(jù)中，各隨機選取4個光譜數(shù)據(jù)繪圖，如圖2所示。

圖2 遙感反射率光譜Fig.2 Spectra of remote-sensing reflectance

1.2.2 吸收系數(shù)和散射系數(shù)測量

吸收系數(shù)ap和ad是通過測量放置在Perkin Elmer Lambda 35紫外/可見分光光度計積分球內(nèi)的樣品過濾器的光譜吸光度獲得的，詳細步驟參考文獻[33-34]。將樣品過濾器放在樣品光束檢測器的一側(cè)，在0°，90°和180°測量3次，得到ap和ad三次掃描的平均值，ag通過UltraPath儀器測量，其中aph=ap-ad，adg=ag+ad。

水體的后向散射系數(shù)bb是用垂直剖面模式放置的原位儀器測量的。在2次ICESCAPE巡航中，bb由2套Hydract-6測量儀器分別進行測量。

1.3 實驗方法

Hoge和Lyon[11-13]開發(fā)的LMI模型利用412，443和555 nm三個波長的遙感反射率，構(gòu)建一個由aph，adg和bbp組成的線性矩陣方程，導(dǎo)出IOPs。LMI模型通過IOPs的共同變化來限制誤差，光譜斜率S的變化范圍在0.01～0.02，冪律指數(shù)Y的變化范圍在0～2，適用于全局優(yōu)化。

Garver和Siegel[15]基于遙感反射率與a和bbp之間的二次關(guān)系開發(fā)了GSM模型，引入了葉綠素a比吸收系數(shù)a*ph(λ)，針對全球海洋葉綠素濃度反演，確定了GSM模型的幾個重要參數(shù)：S=0.020 6，Y=1.037 0，在443～555 nm處，a*ph(λ)分別為0.006 65，0.055 82，0.020 55，0.019 10和0.010 15；Maritorena等[30]采用全局最小化技術(shù)，改進Garver和Siegel的模擬退火模型。GSM模型需要滿足3個條件才能運行：0

Lee等[6]開發(fā)的QAA模型，基于輻射傳輸方程的理論分析和數(shù)值模擬，導(dǎo)出遙感反射率與IOPs之間的關(guān)系。QAA模型主要分為2部分：第1部分模型反演得到總吸收系數(shù)a和后向散射系數(shù)bbp；第2部分是根據(jù)解析方程將a分解為aph和adg。QAA模型已更新了好幾個版本，本次研究采用最新版本的QAA_V6，步驟如下所示。

將水面遙感反射率Rrs(λ)轉(zhuǎn)換為水面下的遙感反射率rrs(λ)[6]：

rrs(λ)=Rrs(λ)/[0.52+1.7Rrs(λ)]。

(1)

Gordon等[9]通過模擬得到了遙感反射率與固有光學(xué)特性的關(guān)系：

rrs(λ)=g0u(λ)+g1u(λ)，

(2)

u(λ)=bb/(a+bb)，

(3)

式中，a表示水體總吸收系數(shù)；bb表示后向散射系數(shù)；g0和g1為模型參數(shù)，分別取值0.089和0.124 5。結(jié)合式(1)和式(2)，u也可以表示為：

(4)

QAA_V6第一部分，選定合適的參考波段：

若Rrs(670)<0.001 5sr-1，選取555 nm作為參考波段，a(λ0)為：

(5)

a(λ0)=a(555)=aw(555)+10-1.146-1.1366χ-0.469 χ2,

(6)

式中，aw(555)是555 nm處純海水的吸收系數(shù)。

若Rrs(670)≥0.001 5sr-1，選取670 nm作為參考波段，a(λ0)為：

(7)

式中，aw(670)是670 nm處純海水的吸收系數(shù)。

bb(λ)可表示為[6]：

bb=bbw(λ)+bbp(λ)，

(8)

(9)

式中，bbw(λ)和bbp(λ)為純海水和懸浮顆粒物的后向散射系數(shù)。

當a(λ0)，u(λ0)和bbw(λ0)已知，bbp(λ0)為：

(10)

光譜斜率η可表示為：

(11)

據(jù)此，可以計算其他波長的bbp(λ)。

將bb(λ)帶入式(3)，可計算其他波段的a(λ)：

(12)

QAA_V6第二部分，據(jù)解析方程分解a為：

a(λ)=aph(λ)+adg(λ)+aw(λ)，

(13)

adg可表示為：

adg(λ)=ag(443)exp[-S(λ-443)]，

(14)

式中，光譜斜率S表示為：

(15)

ag(443)可表示為：

(16)

式中，ζ和ξ為：

(17)

ξ=exp[S(442.5-415.5)]。

(18)

據(jù)此，可以計算其他波長的adg和aph，即：

adg(λ)=ag(443)exp[-S(λ-443)]，

(19)

aph(λ)=a(λ)-adg(λ)-aw(λ)。

(20)

1.4 反演結(jié)果評估

根據(jù)確定系數(shù)(R2)、均方根誤差(Root Mean Square Error，RMSE)和平均相對誤差(Mean Relative Error，MRE)來評估模型性能：

(21)

(22)

(23)

2 實驗結(jié)果分析

本文共77個數(shù)據(jù)，但是研究發(fā)現(xiàn)，LMI模型反演的3個吸收系數(shù)均有負值，約占了總數(shù)據(jù)的1/7；GSM模型受3個條件的約束，有2個無效值；QAA_V6模型反演的aph在5個波段出現(xiàn)多個負值，與Zheng等[31]用QAA_V5模型反演得到的結(jié)論相似，本文基于有效數(shù)據(jù)值對3個模型進行評估。

采用對數(shù)坐標，如圖3所示。

圖3 412，443，490，510和555 nm處LMI,GSM和QAA模型反演和實測a的散點圖

圖3比較了LMI、GSM和QAA_V6三個模型在412，443，490，510和555 nm五個波段處反演a的性能。3個模型反演的a基本散落在1∶1的擬合線附近，其中，QAA模型反演結(jié)果最為理想，如表1所示，在412～555 nm處，R2均大于0.80，RMSE在0.13～0.63之間波動，MRE在0.33～0.68之間波動，但是在555 nm處，圖3(e)高值區(qū)域，明顯高估；GSM模型嚴重低估a(443)，而且在a(412)處，誤差較大；LMI模型反演的a在含量較高區(qū)域出現(xiàn)負值，且相關(guān)性隨波長的增加而降低，在555 nm處相關(guān)性急劇下滑，其RMSE和MRE分別為0.03，0.13。

QAA模型反演的adg與實測數(shù)據(jù)相關(guān)性最好，如圖4所示。

表1 反演a(λ)，adg(λ)和aph(λ)的統(tǒng)計

圖4 412，443，490，510和555 nm處LMI，GSM和QAA模型反演和實測adg的散點圖

在adg含量高的區(qū)域，反演結(jié)果均高估，且對比GSM和LMI，RMSE和MRE精度較低，特別是在412 nm處，RMSE高達1.03；GSM模型導(dǎo)出的adg表現(xiàn)也較好，但均表現(xiàn)為低估，在555 nm處，其MRE為0.76；LMI反演負值出現(xiàn)在含量較高的區(qū)域，除掉10個無效值后，LMI模型反演的adg，其R2均大于0.70，誤差波動也較少。

相比于adg，3個模型反演的aph誤差都較大。QAA模型反演的aph，其R2均低于0.60，且MRE誤差大，555 nm處高達19.52；如圖5所示，GSM模型反演的aph基本圍繞1∶1擬合線附近，分布較為分散，R2均小于0.40，但是精度較穩(wěn)定；LMI模型反演的aph，其R2在0.30左右，RMSE在0.05～0.16間波動，MRE在1.5～3.2間波動，反演結(jié)果均高于實測數(shù)據(jù)且誤差較大。從表1可以看出，3個模型均在555 nm處的反演結(jié)果較差。

圖5 412，443，490，510和555 nm處LMI，GSM和QAA模型反演和實測aph的散點圖

綜合來看，QAA模型反演結(jié)果最為理想，GSM模型次之，LMI模型在北冰洋的反演結(jié)果較不理想，而且3個模型在反演a，adg和aph時，均在555 nm處得到較差的反演結(jié)果。

3 討論與改進

3.1 討論

GSM模型是應(yīng)用于全球大尺度的反演模型，但是應(yīng)用于高緯度地區(qū)時并不準確。分析數(shù)據(jù)發(fā)現(xiàn)，從開闊海域到沿海海域，adg(443)的值變化很大，范圍分別是0.02～0.22 m-1，0.49～1.31 m-1，約占總吸收系數(shù)的77+26%，占據(jù)了主導(dǎo)地位，與低緯度地區(qū)的情況大不同，這與Matsuoka等[32]的研究相一致。模型反演的a和adg在各個波段整體表現(xiàn)為在高值區(qū)域高估，在低值區(qū)域低估，尤其低估a(443)，但是反演的aph基本分散1∶1擬合線。而且不同于前人在低緯度水域的研究[16]，GSM模型沒有顯著低估aph(412)?？梢钥闯?，模型反演的a和adg的變化趨勢相似，而adg的反演精度依賴于參數(shù)S和Y的精度。由于顆粒的散射和高CDOM吸收，光在開闊海域的穿透深度很淺，但在沿海區(qū)域很深。因此，需要為海洋和沿海水域分別設(shè)置S和Y，Matsuoka[30]等研究了波弗特海南部海域的光譜特征，發(fā)現(xiàn)開闊海域在波長為555 nm或者更長處有峰值，而沿海水域則在412～490 nm處有峰值，如圖2所示，依次將77個數(shù)據(jù)劃分為2簇，對海洋和沿海水域分別設(shè)置S和Y參數(shù)，進行靈敏度分析，優(yōu)化GSM模型。依次將本次77個數(shù)據(jù)代入，發(fā)現(xiàn)反演的aph的R2明顯增加，在412 nm處，R2提高為0.56；在443 nm處，R2提高為0.63；在490 nm處，R2提高為0.63；在510 nm處，R2提高為0.58；在555 nm處，R2提高為0.53，盡管R2有所提高，但是MRE誤差較大，在1.06～1.96區(qū)間波動；在412～555 nm范圍內(nèi)，模型反演的adg，其R2穩(wěn)定在0.60左右；但是反演的a(443)明顯改善，其R2由0.62提高為0.96，RMSE由0.25降為0.08，MRE由0.58降為0.19。改進的GSM模型在應(yīng)用于北冰洋時，盡管性能有所提高，但是整體的精度仍不高，可能是受到數(shù)據(jù)集的限制。

QAA_V6模型先根據(jù)經(jīng)驗公式反演獲得a(λ0)和bb(λ0)，再根據(jù)半分析模型推算得到adg(λ)和aph(λ)。分析可知，在a和adg時，相關(guān)性較好，但是精度有待改進；在反演aph時，結(jié)果并不理想。這是因為經(jīng)驗公式通常只適用于光學(xué)特性與模型開發(fā)中使用的光學(xué)特性相似的水域，因此驗證模型在北冰洋的適用性，對模型進行大量的評估，并根據(jù)結(jié)果優(yōu)化模型性能是必要的。

3.2 QAA_V6模型改進

本文從77個站點中隨機選取23個站點，對QAA_V6進行改進，再用其余站點進行驗證。

QAA_V6模型是針對開闊海域和沿海海域提出的，根據(jù)遙感反射率的不同，選取不同的參考波段。圖3中的a(555)明顯高估，特別是高吸收區(qū)域。Lee等[6]發(fā)現(xiàn)在高吸收的渾濁水域，使用大于555 nm的參考波段精確度更高，但是，高緯度地區(qū)不同于低緯度地區(qū)，氣候嚴寒，極晝極夜現(xiàn)象使得人類與植物的活動痕跡較少，Zheng等[31]基于QAA_V5反演的北冰洋的a(555)與實測數(shù)據(jù)吻合良好，在670 nm處的相關(guān)性卻較差。根據(jù)張亭祿等水質(zhì)分類，77個站點中一類水體約占71%，而且實測a(555)>0.24 m-1的站點約占8%，水質(zhì)清澈，再加上Rrs(670)遠小于Rrs(555)，過小的數(shù)值使得實測數(shù)據(jù)的檢測較為困難，而實測數(shù)據(jù)的準確度是優(yōu)化模型的關(guān)鍵。QAA_V6優(yōu)化模型根據(jù)Austin和Petzold的方法，選取對水質(zhì)敏感的a(443)i來估算a(555)，參數(shù)根據(jù)吸收系數(shù)進行調(diào)整[6]：

(24)

412，443，490，510和555 nm處QAA_V6和改進模型反演和實測a的散點圖如圖6所示。

由圖6可以看出，改進后的a(555)結(jié)果明顯改善，特別是在高值區(qū)域，R2由0.81提高為0.91，RMSE由0.20降為0.03，MRE由0.68降為0.13。延續(xù)QAA_V6模型的分類方法，將數(shù)據(jù)分為2簇，選取不同的參考波段。在ICESCAPE中，可以根據(jù)實測的后向散射系數(shù)直接擬合η值，共14個數(shù)據(jù)值，反演的bb(443)相關(guān)性由0.66到0.84；MALINA數(shù)據(jù)缺少實測后向散射系數(shù)，將實測的總吸收系數(shù)代入到QAA_V6模型，讓η在0～3之間迭代，直到反演與實測的a在5個波段中的平均相對誤差小于0.20，迭代停止，得到參考η值。當Rrs(670)<0.001 5 sr-1，η與rrs(443)/rrs(555)的R2由0.47提高為0.57；當Rrs(670)≥0.001 5 sr-1，η與rrs(443)/rrs(555)的R2提也提高為0.57。根據(jù)a(555)反演的其他4個波段的總吸收系數(shù)也緊緊圍繞在1∶1擬合線附近。

QAA_V6模型將a(λ)劃分為adg(λ)和aph(λ)依賴參數(shù)ζ和ξ。通過擬合光譜比率(rrs(443)/rrs(555))與實測數(shù)據(jù)aph(412)/aph(443)求得ζ，adg(412)/adg(443)求取S，進而得到ξ：

(25)

ξ=exp[S(442.5-415.5)]。

(26)

412，443，490，510和555 nm處QAA_V6和改進模型反演和實測adg的散點圖如圖7所示。優(yōu)化模型反演的adg(443)改善明顯，尤其是在高值區(qū)域，基于adg(443)反演的adg(λ)與實測數(shù)據(jù)的相關(guān)性也取得較為滿意的結(jié)果，R2均大于0.85。

圖6 412，443，490，510和555 nm處QAA_V6和改進模型反演和實測a的散點圖

圖7 412，443，490，510和555 nm處QAA_V6和改進模型反演和實測adg的散點圖

QAA_V6模型反演的aph出現(xiàn)了大量負值，約占總數(shù)據(jù)的20%～43%，在統(tǒng)計分析中會被刪除。本文對誤差源進行敏感性分析，按照控制變量的原則，對參數(shù)進行驗證，發(fā)現(xiàn)參數(shù)η是QAA模型反演出現(xiàn)負值的主要原因，與Zheng等[31]的結(jié)果相一致，模型反演的bb(λ)一般會高估短于551 nm波長的后向散射和吸收系數(shù)，并可能在北冰洋的反演中產(chǎn)生負值。優(yōu)化后的aph負值個數(shù)約占總數(shù)據(jù)的18%～35%。412，443，490，510和555 nm處QAA_V6和改進模型反演和實測aph的散點圖如圖8所示。由圖8可以看出，aph(555)整體下降，更接近于實測數(shù)據(jù)，且在412～510 nm處，反演結(jié)果均有所改善，R2提高為0.60左右，但是在555 nm處，R2為0.29，仍不理想。優(yōu)化模型反演的a和adg與實測數(shù)據(jù)較為吻合，但是aph和實測數(shù)據(jù)之間產(chǎn)生了較大的偏差，很大原因是受氣候影響，在北冰洋中，有色溶解和碎屑顆粒物對總吸收系數(shù)的貢獻更大。

圖8 412，443，490，510和555 nm處QAA_V6和改進模型反演和實測aph的散點圖

QAA_V6模型優(yōu)化前后3個吸收系數(shù)的評價指標如圖9所示。由圖9可以看出，adg(λ)的R2有所下降，但整體保持在0.85以上，aph(λ)的R2明顯提高；除了adg(555)的MRE稍增大，RMSE和MRE均取得了令人滿意的結(jié)果。其中，adg(412)，adg(443)，aph(555)和a(555)的誤差大幅下降。

3個模型反演的3個吸收系數(shù)中，浮游植物吸收系數(shù)與實測數(shù)據(jù)差異最大，很大可能是由于模型的性能。早先的研究表明，局部生物光學(xué)模型應(yīng)優(yōu)先于全局模型[36]，但是優(yōu)化參數(shù)，提高模型的整體性能并不簡單，因為模型是高度非線性的，有大量的潛在的聯(lián)系。此外，400～520 nm是浮游植物中多種輔助色素的吸收譜段，導(dǎo)致浮游植物吸收系數(shù)與葉綠素a的相關(guān)性降低[35]。輸入數(shù)據(jù)Rrs(λ)的準確性對模型反演的結(jié)果至關(guān)重要，但是，實測數(shù)據(jù)通常包含由各種實驗和環(huán)境因素引起的大量不確定性，這也是模型反演產(chǎn)生誤差的一個不可避免的來源。

4 結(jié)束語

基于北冰洋的數(shù)據(jù)，測試LMI，GSM和QAA_V6三個反演模型在北冰洋的適用性。研究發(fā)現(xiàn)：① QAA_V6模型反演結(jié)果最好，GSM次之，LMI模型的反演結(jié)果出現(xiàn)多個負值，可能不適合北冰洋的反演；② 3個模型均在555 nm處得到較差的反演結(jié)果；③ 根據(jù)北冰洋特有的光學(xué)特性優(yōu)化的GSM模型，其反演結(jié)果仍不如QAA_V6模型；④ 用555 nm代替670 nm作為參考波段，QAA_V6優(yōu)化模型反演的a與實測數(shù)據(jù)的相關(guān)性更高。

精確的參數(shù)是模型優(yōu)化成功的關(guān)鍵。本文建立了適用于北冰洋的QAA_V6優(yōu)化模型，為水域反演提供了更全面的參考，但水域的光學(xué)性質(zhì)變化很大，受到實測點數(shù)和采樣環(huán)境的限制，未來還需要更多的北冰洋的實測數(shù)據(jù)來提高模型的性能。