光衰減系數(shù)與懸浮顆粒物濃度的關(guān)系＊

王魯寧，魏 皓，趙 亮

（1.中國海洋大學(xué)海洋環(huán)境學(xué)院，山東 青島266100；2.天津科技大學(xué)海洋科學(xué)與工程學(xué)院，天津300457）

海水的光學(xué)性質(zhì)對(duì)于研究海洋生態(tài)系統(tǒng)的結(jié)構(gòu)與功能尤其是初級(jí)生產(chǎn)力有重要意義，而海水中的懸浮顆粒物對(duì)于光輻射的影響在國際上也一直是1個(gè)重要的研究課題，例如 HEBBLE（the High Energy Benthic Boundary Layer Experiment）[1]和 PROVESS（the Processes of Vertical Exchange in Shelf Seas）[2－3]都將海水懸浮顆粒物動(dòng)力學(xué)作為研究重點(diǎn)。

海洋中的浮游植物用來進(jìn)行光合作用的那部分太陽輻射我們稱為光合作用有效光強(qiáng)PAR（photosynthetically available radiation），其波段為400～700nm，對(duì)PAR的度量有2種計(jì)量系統(tǒng)[4]：一種是能量系統(tǒng)，即光合有效輻照度，單位是w/m2；另一種是量子系統(tǒng)，即光合有效量子通量密度，單位是 mol·m－2·s－1）。本文中單位采用后者。在海洋生態(tài)模式中PAR是1個(gè)重要的強(qiáng)迫因子，海表的PAR可以從衛(wèi)星遙感資料得到，但PAR隨深度的垂直分布通常是利用PAR在水中的指數(shù)衰減規(guī)律得到[5－6]：

式中：k為衰減系數(shù)；I0為海面處的PAR值；z為海面以下深度（z＜0）?？梢?，PAR在海水中衰減系數(shù)k是海洋生態(tài)模型的關(guān)鍵參數(shù)。造成PAR在海水中衰減的物理過程主要為反射與吸收[7]，而海水中能夠反射、吸收PAR的組分可以分為純海水、可溶物質(zhì)與懸浮顆粒物[8－10]，三者以懸浮顆粒物影響最大，特別是在中國近海這樣典型的高濁度海域，懸浮顆粒物的濃度與衰減系數(shù)之間的關(guān)系也成為研究重點(diǎn)。Bader[11]認(rèn)為不同粒徑的懸浮顆粒物在水中的分布滿足關(guān)系式：

式中：r為粒徑累積斜率（cumulative slope）；ND為單位體積海水中大于直徑為D的粒徑個(gè)數(shù)。在懸浮顆粒物分布穩(wěn)定的情況下，懸浮物的體積濃度與光的衰減系數(shù)間存在線性相關(guān)性。對(duì)于分布性質(zhì)穩(wěn)定的懸浮物，累積斜率的標(biāo)準(zhǔn)值為3[12－13]。在此基礎(chǔ)上，很多工作相繼展開，McCave[13]與 Peterson[14]研究得到，累積斜率的值在3.0±0.3之間時(shí)，懸浮物體積濃度與衰減系數(shù)之間存在較好的線性關(guān)系：P＝－203.4＋496.0k[13]與P＝－265.1＋569.5k[14]。然而兩者工作僅限于對(duì)長(zhǎng)波段（λ＝660nm）光衰減的研究，且所用懸浮顆粒物粒徑范圍較窄，基本分布在1～32μm之間。本文利用測(cè)得的同一深度剖面上PAR與懸浮顆粒物體積濃度，建立了PAR的衰減系數(shù)與懸浮顆粒物總濃度之間的相關(guān)關(guān)系，并進(jìn)一步將32種粒徑大小不一的懸浮顆粒物分為3種粒級(jí)，建立起了各粒級(jí)與衰減系數(shù)之間的相關(guān)關(guān)系。

1 數(shù)據(jù)資料與觀測(cè)方法

本文所用數(shù)據(jù)資料來源于4次出海觀測(cè)：2005年3月自然科學(xué)重大基金項(xiàng)目“上層海洋與低空大氣”（SOLAS）黃海第一航次觀測(cè)中的大面站與2個(gè)連續(xù)站的觀測(cè)；國家重點(diǎn)基礎(chǔ)研究發(fā)展規(guī)劃項(xiàng)目“我國近海生態(tài)系統(tǒng)食物產(chǎn)出的關(guān)鍵過程及其可持續(xù)機(jī)理”第二課題“近海重要界面物質(zhì)交換的關(guān)鍵過程”的3個(gè)航程：2005年12月14～15日膠州灣定點(diǎn)連續(xù)觀測(cè)；2006年5月4～5日桑溝灣定點(diǎn)連續(xù)觀測(cè)；2007年4月黃海春季水華調(diào)查。在每個(gè)測(cè)站利用加拿大RBR公司生產(chǎn)的XR－620CTD所集成的球形光學(xué)探頭觀測(cè)得到PAR的剖面分布，同時(shí)將現(xiàn)場(chǎng)激光粒度儀LISST－100與CTD一同施放，觀測(cè)得到與PAR同一剖面上1.36～230μm范圍內(nèi)32種粒徑的懸浮顆粒物體積濃度（μL/L）的分布。圖1為這4次出海觀測(cè)站位的分布。

圖1 4個(gè)航次觀測(cè)站位圖Fig.1 The site locations in four investigations

1.1 2005年3月大面與定點(diǎn)連續(xù)觀測(cè)

SOLAS黃海第一航次觀測(cè)區(qū)域在（32.50°N～39.50°N，120.50°E～124.50°E）（見圖1a）。在17個(gè)大面站與2個(gè)25h連續(xù)觀測(cè)站A1、A2處進(jìn)行了PAR與懸浮物濃度的剖面觀測(cè)。A1站觀測(cè)時(shí)間（3月26日15：00～27日15：00）；A2站觀測(cè)時(shí)間（3月20日15：00～27日15：00）。此航次得到PAR與懸浮物濃度的數(shù)據(jù)共有21組。

1.2 2005年12月定點(diǎn)連續(xù)觀測(cè)

2005年12月14 ～15日利用中國海洋大學(xué)東方紅2號(hào)考察船于青島外海（36.04°N，120.32°E）處（見圖1b）進(jìn)行了1個(gè)周日定點(diǎn)連續(xù)觀測(cè)，觀測(cè)間隔1h，共得到10組數(shù)據(jù)。

1.3 2006年5月定點(diǎn)連續(xù)觀測(cè)

2006年5月4 ～5日于榮成桑溝灣內(nèi)（37.13°N，122.55°E）站點(diǎn)（見圖1c）做了1次定點(diǎn)連續(xù)觀測(cè)，觀測(cè)間隔1h，得到15組數(shù)據(jù)。

1.4 2007年4月大面與定點(diǎn)連續(xù)觀測(cè)

2007年4月份利用北斗號(hào)考察船在黃海進(jìn)行了為期25天的觀測(cè)（見圖1d），此次在大面觀測(cè)其間做了3次定點(diǎn)連續(xù)觀測(cè)，共得到35組數(shù)據(jù)。

以上4次觀測(cè)均采用船基吊放方式剖面觀測(cè)，XR－620CTD與 LISST－100同時(shí)施放，其中 XR－620CTD的數(shù)據(jù)采樣頻率為6Hz（水華航次為0.3Hz），LISST－100的數(shù)據(jù)采樣頻率為1Hz。在施放儀器時(shí)，已盡量保證了儀器保持在船向陽側(cè)，以便接受充足的外部光照。

2 數(shù)據(jù)處理與分析

2.1 光衰減系數(shù)的計(jì)算

首先對(duì)PAR的原始觀測(cè)數(shù)據(jù)做物理剔除，得到各測(cè)站PAR隨深度的剖面變化，此處分別選取SOLAS航次中連續(xù)站A1的3月26日23時(shí)與27日11時(shí)2次觀測(cè)PAR的剖面圖（見圖2）為例。通過對(duì)比發(fā)現(xiàn)，23時(shí)測(cè)得的PAR隨深度變化不大，主要分布在3～5個(gè)單位之間，一般認(rèn)為是儀器背景值；而11時(shí)測(cè)得的PAR隨深度變化較為規(guī)律，基本呈指數(shù)分布，可見水體中PAR的值首先取決于外部光學(xué)環(huán)境（日照條件）。因此在隨后分析中同測(cè)站的PAR與懸浮顆粒物濃度都只采用日照充分的白天時(shí)段數(shù)據(jù)，篩選后得到的數(shù)據(jù)對(duì)共有77組。

圖2 SOLAS航次定點(diǎn)觀測(cè)站A1站在26日23時(shí)與27日11時(shí)PAR的剖面分布Fig.2 The profiles of PAR at the mooring site A1 （SOLAS）at 23：00（26th，March）and 11：00（27th，March）separately

在對(duì)所得數(shù)據(jù)作平滑處理后，利用PAR指數(shù)衰減規(guī)律對(duì)PAR值做擬合，擬合時(shí)所需數(shù)據(jù)集的深度區(qū)間選取基本遵循以下原則：取表層1m以下的濃度值，以降低船體對(duì)PAR測(cè)量的影響；另外PAR在水中的衰減主要集中在真光層內(nèi)，真光層以外的深度PAR值很小可認(rèn)為是光場(chǎng)背景值。海洋生態(tài)學(xué)上定義光強(qiáng)僅為表層光強(qiáng)1%的深度為真光層，海洋中的浮游植物光合作用主要集中在這一深度層內(nèi)，同時(shí)也是海洋浮游生物活動(dòng)最為活躍的區(qū)域[15]。根據(jù)定義與各站PAR隨深度的分布，可以求得各站真光層深度hk（如果hk超過測(cè)站水深hmax，則認(rèn)為hk＝hmax），圖3所示所有測(cè)站真光層深度hk與站點(diǎn)水深hmax具有較好的線性相關(guān)（R2＝0.85）。

對(duì)各站真光層內(nèi)PAR值做擬合后得到PAR隨深度指數(shù)衰減曲線以及衰減系數(shù)k。仍然以A1站27日11時(shí)為例，圖2中虛線為該測(cè)站PAR的擬合曲線，衰減指數(shù)為0.3/m。

圖3 所有測(cè)站真光層深度hk與測(cè)站水深hmax的關(guān)系Fig.3 The relation of the euphotic layer depth with the maximum depth form all situs

作者將膠州灣與桑溝灣定點(diǎn)觀測(cè)作為在近岸海區(qū)的觀測(cè)，SOLAS航次以及水華航次觀測(cè)作為外海區(qū)觀測(cè)。對(duì)比近岸與外海PAR衰減系數(shù)的分布（見圖4）可以看到：近岸海區(qū)PAR衰減系數(shù)較大，分布在0.5～1之間，說明近岸海區(qū)PAR衰減的較快；外海區(qū)PAR衰減系數(shù)較小，在0.1～0.4之間，只有個(gè)別觀測(cè)有較大值。

圖4 近岸海區(qū)與外海區(qū)PAR衰減系數(shù)相對(duì)頻數(shù)的分布Fig.4 The distributions of frequency of attenuation coefficient in coastal and open waters separately

2.2懸浮顆粒物體積濃度分布特征

LISST－100測(cè)得的懸浮顆粒物濃度數(shù)據(jù)的初步處理同PAR類似。根據(jù)Soulsby[16]提出的溫氏分級(jí)表，本文進(jìn)一步將32種粒徑分為3種粒級(jí)：黏土（d≤12μm）；粉砂（d＝12～62μm ）；細(xì)砂 （d≥62μm），d為粒徑大小。

圖5為SOLAS航次A1站3月27日11時(shí)測(cè)得的3種粒徑體積濃度的剖面，其中同一深度層懸浮顆粒物體積濃度細(xì)砂＞粉砂＞黏土，濃度主要取決于粒徑大小。黏土與粉砂體積濃度隨深度的變化很小，整個(gè)剖面分布比較均勻，而細(xì)砂剖面分布比較雜亂且體積濃度隨深度增大。根據(jù)Bader[11]提出的累計(jì)斜率的定義，經(jīng)過擬合可分別求得各站每一深度層上所有懸浮顆粒物的累積斜率值r，對(duì)每一站的r在真光層內(nèi)做深度平均，可得到每個(gè)站點(diǎn)所有懸浮顆粒物累積斜率值的平均值（見圖6）?？梢姡姓军c(diǎn)的基本成正態(tài)分布，且主要分布在2.4～3.0之間，平均值為2.71，與 Mc－Cave[13]與Peterson[14]等得到的值較為接近，說明觀測(cè)站點(diǎn)懸浮顆粒物的分布基本穩(wěn)定，可推測(cè)懸浮物顆粒濃度與光衰減系數(shù)之間存在良好的線性相關(guān)性。同時(shí)，通過對(duì)比近岸與外海海區(qū)3種粒級(jí)懸浮顆粒物累計(jì)斜率的分布（見圖7）可以發(fā)現(xiàn)：近岸黏土的累計(jì)斜率平均值分布在2～2.5之間，粉砂則在2.5～3.5之間，細(xì)砂的值則主要分布在2.5～4之間；外海3種粒級(jí)懸浮物的累計(jì)斜率平均值基本呈正態(tài)分布，其中黏土r值主要分布在2～3之間，而粉砂、細(xì)砂的r值變化范圍較大，分布在2～5之間。通過對(duì)比發(fā)現(xiàn)，不管在近岸還是外海區(qū)粉砂類懸浮物的分布最穩(wěn)定，細(xì)砂類分布最不穩(wěn)定，而外海區(qū)黏土類懸浮物的分布相對(duì)于近岸要更穩(wěn)定。

圖5 SOLAS航次A1站27日11時(shí)3種粒級(jí)懸浮顆粒物濃度的剖面分布Fig.5 The profiles of concentration of suspended particles of 3classes at mooring site A1（SOLAS），at 11：00on 27th，March of 2005

圖6 全部站點(diǎn)真光層內(nèi)所有粒徑懸浮顆粒物平均累積斜率值相對(duì)頻數(shù)的分布Fig.6 The distribution of frequency of cumulative slope of all suspended particles in euphotic layers of all sites

圖7 近岸（膠州灣、桑溝灣定點(diǎn)測(cè)站）與外海（2005年3月、2007年4月黃海站點(diǎn)）海區(qū)3種懸浮顆粒物粒級(jí)累積斜率平均值的分布Fig.7 The distributions of cumulative slope of suspended particle of 3classes in coastal and open waters separately

3 結(jié)果與討論

3.1 PAR衰減系數(shù)與懸浮顆粒物濃度的關(guān)系

對(duì)各站點(diǎn)懸浮顆粒物濃度特征值的選取較為關(guān)鍵，Karen[2]利用12L采水瓶同光學(xué)儀器同層取水，選取深度層在上表層（大約3m左右）的濃度為測(cè)站特征值；Bowers[17]則每站基本采3層水樣（近表層、中層或躍層、近底層），取表層水樣的濃度作為所測(cè)站點(diǎn)的懸浮顆粒物濃度特征值。同2.1節(jié)中處理PAR的原則相同，本文選取對(duì)真光層深度內(nèi)懸浮顆粒物濃度做積分平均后的濃度值Pk作為測(cè)站的特征值。

本文選取桑溝灣觀測(cè)與2007年4月黃海航次觀測(cè)數(shù)據(jù)作分析，得到PAR衰減系數(shù)K與真光層內(nèi)所有粒徑懸浮顆粒物濃度Pk（mL/L）之間的相關(guān)關(guān)系（式3.1），圖8所示所有測(cè)站PAR衰減系數(shù)與懸浮顆粒物濃度存在較好的線性相關(guān)關(guān)系（R2＝0.80），但分別在近岸區(qū)（桑溝灣觀測(cè)）與外海區(qū)（黃海航次觀測(cè)）測(cè)站兩者相關(guān)性較較差（R2＝0.67；R2＝0.43），同時(shí)近岸海區(qū)的懸浮顆粒物濃度普遍高于外海海區(qū)，這主要與近岸海區(qū)水動(dòng)力環(huán)境以及陸源輸入有關(guān)。

圖8 RAR衰減系數(shù)與真光層內(nèi)所有粒徑懸浮顆粒物濃度的關(guān)系Fig.8 The relationship of PAR attenuation coefficient with concentrations of all suspended particles in euphotic layer

2.2節(jié)中本文將測(cè)得的32種不同粒徑的懸浮顆粒物分為3種粒級(jí)（黏土、粉砂、細(xì)砂），這3種粒級(jí)對(duì)PAR衰減的影響由于濃度、組成成分以及分布性質(zhì)的不同而不同，需要對(duì)這3種粒級(jí)進(jìn)行對(duì)比分析。圖9給出了真光層內(nèi)3種粒級(jí)懸浮顆粒物的平均濃度分別與PAR衰減系數(shù)K之間的相關(guān)關(guān)系。經(jīng)過對(duì)比不難發(fā)現(xiàn)，黏土、粉砂濃度P1、P2與PAR衰減系數(shù)具有較好的相關(guān)性（R2＝0.87；R2＝0.87），而細(xì)砂濃度P3與衰減系數(shù)的相關(guān)性較差（R2＝0.59）。

圖9 三種粒級(jí)體積濃度分別與PAR衰減系數(shù)的相關(guān)性比較Fig.9 The relationships of PAR attenuation coefficient with concentration of suspended particles of 3classes separately

由式3.1得到，當(dāng)海水中懸浮顆粒物濃度Pk的值為零，波長(zhǎng)范圍在400～700nm的PAR衰減系數(shù)K的值為0.14m－1，這是由純海水、黃色物質(zhì)以及1.36～230μm范圍以外的懸浮物造成的衰減，明顯高于楊生光[18]在黃海千里崖海域測(cè)得的純海水消光系數(shù)值（0.065 4m－1）。另外，Spinrad[1]和 Peterson[14]分別求得的660nm單一波段的光在純水中的衰減系數(shù)值為0.4與0.47m－1，遠(yuǎn)高于PAR在純水中的衰減，反映了不同波段的光在水中衰減的差異。

由式3.2、3.3、3.4得到，海水中粒級(jí)較小的黏土（≤12μm）與粉砂（12～62μm）類懸浮顆粒物濃度與PAR衰減系數(shù)間呈較好的線性相關(guān)，而粒級(jí)較大的細(xì)砂（＞62μm）類懸浮顆粒則與光衰減相關(guān)較差。Pak[19]認(rèn)為光的衰減主要是由小于20μm的顆粒組分控制的，Siegel[20]也認(rèn)為粒徑小于5μm的顆粒濃度分布能更好地反映出光衰減系數(shù)的分布，都突出了細(xì)顆粒物質(zhì)對(duì)于光衰減的重要性。

3.2 結(jié)果的驗(yàn)證

本文利用SOLAS走航測(cè)站與膠州灣定點(diǎn)測(cè)站數(shù)據(jù)分別對(duì)關(guān)系式3.1、3.2、3.3、3.4做精度檢驗(yàn)，用判定系數(shù)R2衡量回歸方程的擬合度，計(jì)算公式為：

式中：yi、Yi分別為數(shù)據(jù)實(shí)測(cè)值與預(yù)報(bào)值。

表1所示為SOLAS航次與膠州灣定點(diǎn)觀測(cè)數(shù)據(jù)與各回歸方程（式3.1、3.2、3.3、3.4）之間的判定系數(shù)，其中近岸膠州灣數(shù)據(jù)與各方程之間的判定系數(shù)值較高，均在0.50左右，說明觀測(cè)數(shù)據(jù)中的誤差有50%可以用回歸方程中的線性關(guān)系解釋。SOLAS航次數(shù)據(jù)與式3.1、3.4的擬合度較好，而與式3.2、3.3擬合度相對(duì)較差，判定系數(shù)分別為0.30、0.39，即只有30%、39%的誤差用回歸方程解釋。

表1 SOLAS航次與膠州灣定點(diǎn)觀測(cè)數(shù)據(jù)與各回歸方程之間的判定系數(shù)R2值Table 1 The determination coefficients（R2）of observations in SOLAS and Jiaozhou Bay to the results form regression equations

4 結(jié)論

本文通過對(duì)2個(gè)航次與2次定點(diǎn)觀測(cè)得到的PAR與懸浮顆粒物體積濃度數(shù)據(jù)進(jìn)行分析，得到以下結(jié)論：

（1）近岸海區(qū)懸浮顆粒物濃度普遍要高于外海區(qū)，PAR的衰減系數(shù)也大于外海區(qū)。不同粒徑的懸浮顆粒物在不同海區(qū)的分布特征不同：近岸與外海海區(qū)的粉砂類懸浮物累計(jì)斜率平均值在3左右，分布較穩(wěn)定；外海區(qū)黏土類懸浮物分布較外海海區(qū)穩(wěn)定；細(xì)砂類懸浮物在近岸與外海的穩(wěn)定性較差。

（2）PAR衰減系數(shù)K與32種粒徑懸浮顆粒物總體積濃度Pk（mL/L）間呈較好的線性相關(guān)：K＝0.92Pk＋0.14。對(duì)于不同粒級(jí)的懸浮顆粒物，PAR衰減系數(shù)與分布較穩(wěn)定的黏土、粉砂類懸浮物濃度線性相關(guān)，與分布較不穩(wěn)定的細(xì)砂類懸浮物濃度的相關(guān)性較差，可見PAR的衰減特性從粒徑較小的懸浮顆粒物分布上得到更好的反映。

（3）光衰減系數(shù)在生態(tài)模型中是一重要的因子，能夠影響浮游植物葉綠素a以及初級(jí)生產(chǎn)力的分布[4]。然而，由于光衰減系數(shù)很難直接觀測(cè)到，以往的生態(tài)模型在光衰減系數(shù)取值上一般采用經(jīng)驗(yàn)估算得到，忽略了光衰減系數(shù)的空間分布。本文建立起海水中光衰減系數(shù)與懸浮顆粒物濃度的關(guān)系，可以利用較易觀測(cè)的懸浮物濃度計(jì)算光衰減系數(shù)，同時(shí)，也能反映出水動(dòng)力環(huán)境下懸浮物濃度空間分布的變化對(duì)光衰減系數(shù)空間分布的影響，為生態(tài)模型提供更準(zhǔn)確的光衰減系數(shù)值。

[1]Spinrad R W，Zaneveld J R V，Kitchen J C.A study of the optical characteristics of the suspended partickes in the benthic nepheloid layer of the scotian rise [J].J Geophys Res，1983，88（C12）：7641－7645.

[2]Karen W A，Andrew L，Paul T.Phytoplankton，sediment and optical observations in Netherlands coastal water in spring[J].J Sea Res，2002，47：303－315.

[3]McCandliss R R，Jones S E，Hearn M，et al.Dynamics of suspended particles in coastal waters during a spring Bloom [J].J Sea Res，2002，47：285－302.

[4]曹文熙，楊躍忠.海洋光合有效輻射分布的計(jì)算模式 [J].熱帶海洋學(xué)報(bào)，2002，21（3）：47－54.

[5]Kirk J T O.Light and Photosynthesis in Aquatic Ecosystem [M].New York：Cambridge University Press，1983.

[6]Morel A.Bio－optical properties of oceanic waters：A reappraisal[J].J Geophys Res，2001，106（C4）：7163－7180.

[7]McKee D，Cunningham A，Slater J，et al.Inherent and apparent optical properties in coastal waters：A study of the Clyde Sea in early summer[J].Estuarine，Coastal and Shelf Science，2003，56：369－376.

[8]Babin M，Stramski D，F(xiàn)errari G M，et al.Variations in the light absorption coefficients of phytoplankton，nonalgal particles and dissolved organic matter in coastal waters around Europe [J].J Geophys Res，2003，108（C7）：3211.

[9]Bracchini L，Dattilo A M，F(xiàn)alcucci M，et al.Spatial and temporal variations of the inherent and apparent optical properties in a shallow coastal lake[J].J Photoche Photobio B：Biology，2005，80：161－177.

[10]Figueroa F L.Bio－optical characteristics of gerlache and bransfield strait waters during an Antarctic Summer Cruise[J].Deep－Sea Research II，2002，49：675－691.

[11]Bader J.The hyperbolic distribution of particle sizes[J].J Geophys Res，1970，75：2822－2830.

[12]Baker E T，Lavelle J W.The effect of particle size on the light attenuation coefficient of natural suspensions[J].J Geophys Res，1984，89（C5）：8199－8203.

[13]McCave I N.Particulate size spectra，behavior and origin of nepheloid layers over the Nova Scotian Continental Rise[J].J Geophys Res，1983，88（C12）：7647－7666.

[14]Peterson R L.A Study of suspended particulate matter：arctic O－cean and Northern Oregon Continental Shelf[D].Oregon：Oregon State Univ，1978.

[15]陳敏，黃奕普.U不平衡法在真光層顆粒動(dòng)力學(xué)研究中的應(yīng)用[J].地球科學(xué)進(jìn)展，1999，14（4）：365－370.

[16]Soulsby R.Dynamics of marine sands[M].London：Thomas Telford，1997：25－32.

[17]Bowers D G，Harker G E L，Smith P S D，et al.Optical properties of a region of freshwater influence [J].Estuarine，Coastal and Shelf Science，2000，50：717－726.

[18]楊生光.千里崖海域真光層光衰減 [J].海洋與湖沼，1992，23（3）：245－251.

[19]Pak H，Kiewfer D A，Kitchen J C.Meridionalvariations in the concentration of chlorophy11and microparticles in the North Pacific Ocena[J].Deep－Sea Research，1988，35：1151－1171.

[20]Siegel D A，Iturriaga R，Bidigare R R，et al.Meridional variations of the springtime phytoplankton community in the Sargasso Sea[J].Journal of Marine Research，1990，48：379－412.