• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    夏季長(zhǎng)江口水體可溶性有機(jī)質(zhì)組成的空間分布特征及其控制因素

    2021-09-15 01:54:32王安月范代讀孫永革
    地球化學(xué) 2021年3期
    關(guān)鍵詞:長(zhǎng)江口站位鹽度

    王安月, 張 凱, 何 丁, 范代讀, 孫永革*

    夏季長(zhǎng)江口水體可溶性有機(jī)質(zhì)組成的空間分布特征及其控制因素

    王安月1, 張 凱1, 何 丁1, 范代讀2, 孫永革1*

    (1. 浙江大學(xué)地球科學(xué)學(xué)院 有機(jī)地球化學(xué)研究組, 浙江 杭州 310027; 2. 同濟(jì)大學(xué) 海洋地質(zhì)國家重點(diǎn)實(shí)驗(yàn)室, 上海 200092)

    本研究通過測(cè)定有色溶解有機(jī)物(CDOM)的吸收光譜、熒光可溶性有機(jī)質(zhì)(FDOM)的激發(fā)-發(fā)射-矩陣三維熒光光譜(excitation-emission-matrix spectra, EEMs)和穩(wěn)定碳同位素組成(13C), 系統(tǒng)探討了長(zhǎng)江口夏季水體可溶性有機(jī)質(zhì)(DOM)的組成、來源、空間分布及河口混合行為等。研究結(jié)果表明, 可溶性有機(jī)碳(DOC)濃度整體呈現(xiàn)由陸到海逐漸降低的趨勢(shì), 表征有色溶解有機(jī)物含量水平的吸收系數(shù)(355)與鹽度呈負(fù)相關(guān)關(guān)系, 指示可溶性有機(jī)質(zhì)中熒光組分在河口的分布主要受稀釋作用調(diào)控。利用EEMs并結(jié)合平行因子分析(PARAFAC)鑒定出代表陸源有機(jī)質(zhì)的類腐殖質(zhì)的熒光組分C2和C3, 以及代表原地生產(chǎn)力的類蛋白質(zhì)熒光組分C1和C4。由近岸到外海, 表征海洋藻類生產(chǎn)力的C1組分在水體熒光有機(jī)質(zhì)中所占比例增高, 陸源信號(hào)則呈現(xiàn)逐漸降低趨勢(shì), 與此相對(duì)應(yīng), 外海水可溶性有機(jī)質(zhì)具有高的光譜斜率275-295/350-400比值和重碳同位素組成?;邴}度.可溶性有機(jī)碳濃度、鹽度.13C值的河流-海洋雙端元混合模型, 發(fā)現(xiàn)長(zhǎng)江口最大混濁帶存在可溶性有機(jī)碳的移除過程, 且表現(xiàn)為以物理稀釋作用為主的缺失性非保守混合行為。研究成果揭示了水動(dòng)力條件是控制長(zhǎng)江口-東海陸架系統(tǒng)可溶性有機(jī)質(zhì)組成和空間分布的關(guān)鍵因素。

    可溶性有機(jī)質(zhì); 碳同位素; 混合行為; 三維熒光光譜-平行因子分析; 長(zhǎng)江口

    0 引 言

    可溶性有機(jī)碳(DOC)是海洋中最大的有機(jī)碳庫, 與大氣中的CO2含量相當(dāng)(相當(dāng)于0.75×1018g 碳)。據(jù)估算, 每年約有0.25×1015g可溶性有機(jī)碳通過河流進(jìn)入到海洋[1], 約占流入海洋總碳量(有機(jī)碳和無機(jī)碳共0.9×1015g)的30%, 更是占到流入海洋有機(jī)碳的60%之多[2]??扇苄杂袡C(jī)碳從近海至遠(yuǎn)洋, 濃度發(fā)生急劇下降, 每年只有約0.1 Gt的可溶性有機(jī)碳在深海中流通, 說明有接近一半的可溶性有機(jī)碳從河口進(jìn)入海洋的過程中發(fā)生了改變[3]。然而, 至今對(duì)于水體可溶性有機(jī)質(zhì)(DOM)從河口至大洋的轉(zhuǎn)化降解形式以及最終歸縮的了解仍十分有限, 尤其是陸源有機(jī)質(zhì)的輸入和遷移轉(zhuǎn)化途徑, 已經(jīng)成為了當(dāng)前海洋碳循環(huán)研究的熱點(diǎn)問題之一。

    河口作為陸地生態(tài)系統(tǒng)和海洋生態(tài)系統(tǒng)的交匯處, 水動(dòng)力條件、鹽度梯度變化較大, 具有水體存留時(shí)間短、溶解態(tài)物質(zhì)周轉(zhuǎn)速率快和生物地球化學(xué)過程復(fù)雜且劇烈等特點(diǎn), 是開展可溶性有機(jī)碳來源和歸宿研究最復(fù)雜的水體之一, 其中可溶性有機(jī)碳在水體中是否遵循保守行為一直存有爭(zhēng)議。理論上, 當(dāng)只有淡水與海水兩端元混合時(shí), 可溶性有機(jī)碳濃度與鹽度之間存在線性關(guān)系, 表現(xiàn)出保守性(不存在添加與去除過程)[4-5]。但越來越多的研究揭示, 可溶性有機(jī)碳在河口存在降解消耗、絮凝沉降、以及與顆粒有機(jī)碳(POC)的相互轉(zhuǎn)化等行為, 使得可溶性有機(jī)碳濃度偏離于混合保守曲線, 通常表現(xiàn)出移除或者添加的行為[6-11]。álvarez-Salgado.[12]在對(duì)Humber河口的研究中發(fā)現(xiàn), 可溶性有機(jī)碳只有在鹽度大于5的區(qū)域表現(xiàn)出良好的保守混合行為, 而當(dāng)鹽度小于5時(shí), 可溶性有機(jī)碳表現(xiàn)出隨鹽度上升而快速的移除過程。張向上[9]在研究黃河口可溶性有機(jī)碳的混合行為時(shí)發(fā)現(xiàn), 在鹽度大于5的區(qū)域同樣有較好的保守性混合行為, 海水稀釋占了主導(dǎo)地位, 與前人不同的是, 在鹽度小于5的區(qū)域, 由于顆粒有機(jī)碳向可溶性有機(jī)碳的轉(zhuǎn)化, 使得可溶性有機(jī)碳濃度高于保守混合線, 表現(xiàn)出增高的趨勢(shì)。Cauwet[13]對(duì)Yugoslavia Krka河口、Callahan.[14]對(duì)珠江口的研究發(fā)現(xiàn), 可溶性有機(jī)碳僅在一個(gè)小范圍的鹽度梯度內(nèi)與鹽度呈現(xiàn)出較好的相關(guān)性, 若以整個(gè)河口水體為研究對(duì)象, 則沒有明顯的保守行為。Medeiros.[11]在觀測(cè)亞馬遜河-西大西洋體系中可溶性有機(jī)碳的動(dòng)態(tài)變化過程中, 發(fā)現(xiàn)大量可溶性有機(jī)碳在河口內(nèi)濃度低于保守混合線, 存在大量去除過程。由此可見, 河口區(qū)內(nèi)部的物理、化學(xué)和生物過程對(duì)可溶性有機(jī)碳濃度的控制可能遠(yuǎn)比外部混合重要。因此, 多因素作用下可溶性有機(jī)碳化學(xué)組成的動(dòng)態(tài)過程成為揭示河口區(qū)可溶性有機(jī)碳濃度動(dòng)態(tài)過程控制因素的重要手段。

    長(zhǎng)江作為中國第一大河流, 其入海徑流量占東海所有河流徑流輸入的90%~95%[15]。長(zhǎng)江口區(qū)域的物質(zhì)輸送受到潮汐活動(dòng)、陸源徑流量以及東海復(fù)雜水系等多重因素的影響, 其中最大渾濁帶的存在, 為可溶性有機(jī)碳物理、化學(xué)和生物反應(yīng)提供了理想的場(chǎng)所, 是河口“過濾器”最突出的表現(xiàn)[16]。Ogawa.[17]沿長(zhǎng)江口沖淡水方向?qū)扇苄杂袡C(jī)碳濃度的空間分布及其控制因素進(jìn)行了系統(tǒng)研究, 結(jié)果顯示, 可溶性有機(jī)碳濃度介于65~75 μmol/L之間, 主要受控于水文條件的變化而非葉綠素a (a)表征的生物量。對(duì)可溶性有機(jī)質(zhì)的光譜研究揭示, 有色溶解有機(jī)物(CDOM)的組成主要受外源輸入、保守性混合以及相態(tài)轉(zhuǎn)化影響。有色溶解有機(jī)物主要熒光組分均表現(xiàn)為近岸高遠(yuǎn)岸低的特征, 熒光組分濃度與葉綠素a以及鹽度的相關(guān)性指示浮游植物對(duì)有色溶解有機(jī)物組成的直接影響較小, 夏季長(zhǎng)江的徑流輸入是該海域有色溶解有機(jī)物的主要來源。同時(shí), 長(zhǎng)江口外區(qū)域表層有色溶解有機(jī)物濃度顯著受潮汐活動(dòng)影響[18-20]??傮w上, 生物來源的有色溶解有機(jī)物比例隨著鹽度的增加而增加[21], 而類腐殖質(zhì)組分濃度隨著鹽度的增加而線性減小, 說明類腐殖質(zhì)組分主要受到咸淡水物理混合的影響[22]。Guo.[23]應(yīng)用保守混合模型提出長(zhǎng)江口的有色溶解有機(jī)物遵守保守混合行為, 物理稀釋作用是影響可溶性有機(jī)碳濃度分布的關(guān)鍵。然而, 平行因子分析(PARAFAC)揭示河口區(qū)表現(xiàn)出非保守性混合行為, 在最大渾濁帶存在有色溶解有機(jī)物添加過程[24]。最大渾濁帶區(qū)域熒光組分濃度的增加也被認(rèn)為可能與沉積物再懸浮作用有關(guān)[22]。另一方面, 長(zhǎng)江口最大混濁帶混合行為還存在時(shí)間和空間上的差異性, 長(zhǎng)江口北港在洪水季表現(xiàn)為保守性混合, 在枯水季則存在顯著的添加, 而難降解的可溶性有機(jī)質(zhì)主要通過南港入??趨R入東海, 黃浦江匯入南港后帶入了大量的生物可利用性較高的可溶性有機(jī)質(zhì), 受人類活動(dòng)影響大[21,25]。由此可見, 長(zhǎng)江口最大渾濁帶對(duì)于可溶性有機(jī)質(zhì)的組成和空間展布起到直接調(diào)控作用, 在最大渾濁帶很可能存在一個(gè)顯著的可溶性有機(jī)碳添加過程。由于基于可溶性有機(jī)碳與鹽度間的線性關(guān)系只能代表輸入與消耗之間的平衡, 不能代表絕對(duì)的保守行為[8]。因此, 至今對(duì)于最大渾濁帶的調(diào)控作用的了解仍十分有限, 關(guān)鍵在于可溶性有機(jī)質(zhì)來源的動(dòng)態(tài)過程難以很好地刻畫。本研究擬利用紫外吸收光譜、激發(fā)-發(fā)射-矩陣三維熒光光譜和平行因子分析(EEMs- PARAFAC)等技術(shù), 并結(jié)合可溶性有機(jī)碳碳同位素組成和保守混合模型, 試圖探討: (1) 夏季長(zhǎng)江口可溶性有機(jī)質(zhì)的組成、來源、空間分布特征及控制因素; (2) 最大渾濁帶對(duì)于陸源物質(zhì)的過濾作用。以期幫助了解可溶性有機(jī)質(zhì)在河口的內(nèi)部循環(huán)機(jī)制, 及其對(duì)于河口區(qū)碳循環(huán)、儲(chǔ)存的影響。

    1 材料與方法

    1.1 樣品的采集與預(yù)處理

    2017年7月24日至30日, 利用同濟(jì)大學(xué)海洋地質(zhì)國家重點(diǎn)實(shí)驗(yàn)室組織的YES航次獲得YE斷面18個(gè)站位表層水樣。該斷面西起121°E徐六涇附近, 東至124.5°E, 南北跨度大, 覆蓋水域廣(圖1)。每個(gè)站位取表層水樣(約3 m深), 用孔徑0.7 μm的GF/F膜(預(yù)先用馬弗爐在450 ℃時(shí)灼燒5 h)和0.2 μm濾膜(millipore)過濾并收集于Nalgene瓶中(提前用pH=2的鹽酸清洗并烘干), 于?20 ℃冷凍保存。

    1.2 可溶性有機(jī)碳濃度的測(cè)定

    水樣帶回實(shí)驗(yàn)室后, 解凍至室溫。一部分用于可溶性有機(jī)碳濃度測(cè)定, 采用高溫催化氧化法, 使用的儀器為日本島津總TOC-L有機(jī)碳分析儀。

    1.3 可溶性有機(jī)質(zhì)富集和穩(wěn)定碳同位素分析

    按照Dittmar.[26]的流程, 采用PPL柱進(jìn)行固相萃取富集可溶性有機(jī)質(zhì), 并用甲醇進(jìn)行洗脫, 洗脫液置于安培瓶中封口, 置于?20 ℃冰箱中保存。

    圖1 長(zhǎng)江口2017年7月采樣站位圖

    可溶性有機(jī)質(zhì)穩(wěn)定碳同位素測(cè)定在連接有Thermo Fisher Flash 2000元素分析儀的Thermo Fisher MAT-253同位素質(zhì)譜儀上完成。標(biāo)準(zhǔn)物質(zhì)尿素監(jiān)控儀器狀態(tài), 每6次進(jìn)樣測(cè)試1次尿素。分析精度小于±0.2‰, VPDB標(biāo)準(zhǔn)。

    1.4 三維熒光光譜測(cè)定

    采用三維熒光光譜儀(Aqualog)對(duì)過濾水樣進(jìn)行三維熒光光譜掃描。光譜激發(fā)波長(zhǎng)范圍為240~600 nm, 間隔3 nm; 發(fā)射波長(zhǎng)為250~550 nm, 間隔3 nm。以Milli-Q水為空白參比, 對(duì)水樣進(jìn)行連續(xù)掃描, 獲得各水樣的光吸收曲線。

    1.5 吸收光譜測(cè)定

    采用紫外-可見光分光光度計(jì)(島津UV-2450, 日本)對(duì)過濾水樣進(jìn)行紫外-可見光吸收光譜分析。掃描波長(zhǎng)范圍為200~700 nm, 步長(zhǎng)1 nm, 光程1 cm, 以Milli-Q水作參比。采用式(1)將吸光度換算為吸收系數(shù)()。

    式中,A為波長(zhǎng)處的吸光度,0為680~700nm的平均吸光度,為實(shí)驗(yàn)用比色皿的寬度。用于校正可溶性有機(jī)質(zhì)過濾液中膠體顆粒引起的光散射效應(yīng)。

    有色溶解有機(jī)物的光吸收特性表現(xiàn)為在紫外可見光區(qū)最高, 到紅外光譜區(qū)降為零, 光譜吸收近似呈指數(shù)衰減的規(guī)律[27], 可用式(2)進(jìn)行表示。

    式中,(0)為有色溶解有機(jī)物的吸收系數(shù);為波長(zhǎng);0為參照波長(zhǎng);為指數(shù)函數(shù)曲線斜率參數(shù), 表示吸收系數(shù)隨波長(zhǎng)降低的程度。

    本研究采用355 nm處的紫外吸收系數(shù)(355)表征水體中有色溶解有機(jī)物含量[10,23,28]。

    光譜斜率比275-295/350-400是將指數(shù)方程(2)經(jīng)自然對(duì)數(shù)轉(zhuǎn)化后, 對(duì)275~295 nm或350~400 nm波長(zhǎng)范圍內(nèi)的吸收系數(shù)進(jìn)行線性擬合而得到。通過275-295和350-400的比值可以獲得斜率比(S), 該值可用于有色溶解有機(jī)物來源識(shí)別。一般認(rèn)為,S越大, 可溶性有機(jī)質(zhì)相對(duì)分子質(zhì)量越小, 海源有色溶解有機(jī)物的成分越大[29]。

    1.6 熒光光譜-平行因子分析

    PARAFAC是一種基于交互最小二乘法獲得數(shù)據(jù)殘差平方和最小的算法, 是一種解析三維熒光光譜常用的方法。PARAFAC利用模型來確定組分的峰位置, 使分析過程簡(jiǎn)便快捷, 并且提高了分析的靈敏度。利用其分析解的唯一性, 可以成功區(qū)分化學(xué)結(jié)構(gòu)相似的熒光組分[30-31]。本研究采用PARAFAC模型對(duì)有色溶解有機(jī)物的三維熒光光譜進(jìn)行分解, 以識(shí)別長(zhǎng)江口水體有色溶解有機(jī)物的熒光組分特征。利用Matlab軟件對(duì)獲得的三維熒光光譜圖進(jìn)行模擬識(shí)別, 共得到4種組分, 其中有2種是類腐殖質(zhì)組分, 2種是類蛋白質(zhì)組分, 通過折半分析法驗(yàn)證結(jié)果的可靠性。

    2 結(jié)果與討論

    2.1 長(zhǎng)江口環(huán)境梯度分區(qū)及其對(duì)可溶性有機(jī)碳濃度的控制作用

    一般認(rèn)為, 長(zhǎng)江沖淡水鹽度值主要介于5~31之間, 主體鹽度值為18~28[32]。夏季受臺(tái)灣暖流影響, 沖淡水會(huì)在122°30′~123°50′E附近轉(zhuǎn)向東北濟(jì)州島方向[33]。本航次水體鹽度實(shí)測(cè)結(jié)果顯示, 沿采樣斷面由西向東表現(xiàn)出鹽度逐漸增加的過程, 其中16站位以東可能受臺(tái)灣暖流的影響, 形成的上升流使得鹽度快速遞增(圖2)。站位15以西, 即河口至123°E范圍內(nèi), 屬于長(zhǎng)江沖淡水區(qū)域。站位14至17區(qū)間同時(shí)受到了長(zhǎng)江沖淡水和臺(tái)灣暖流的影響, 鹽度梯度較小, 位于水舌邊緣(圖1)。由水體表底層濁度分布圖可以發(fā)現(xiàn)在站位9至13范圍內(nèi)濁度值出現(xiàn)高值, 指示“最大渾濁帶”的存在(圖2)。為了便于討論不同水文條件對(duì)可溶性有機(jī)質(zhì)混合行為的影響, 將研究區(qū)域分為4部分(圖1)。河流區(qū): 站位1至8, 鹽度接近于0; 最大渾濁帶: 站位9至13, 以表底層水體高濁度為特征, 光照弱, 不利于初級(jí)生產(chǎn)者的生長(zhǎng); 沖淡水邊緣: 122°30′E與123°30′E之間, 包含站位14至17; 外海區(qū): 站位18以東的廣闊海域。

    由圖2可見, 可溶性有機(jī)碳濃度整體呈現(xiàn)由陸到海逐漸降低的趨勢(shì)。光譜研究進(jìn)一步揭示, 表征有色溶解有機(jī)物含量水平的吸收系數(shù)(355)與鹽度呈負(fù)相關(guān)(圖3a), 與大多數(shù)河口研究結(jié)果相一致[25,34-36], 說明可溶性有機(jī)質(zhì)中熒光組分在河口的分布主要受稀釋作用調(diào)控[19,37,38]。8站位有色溶解有機(jī)物出現(xiàn)異常高值, 可能是由于8站位位于沖淡水與渾濁帶的交匯處, 水動(dòng)力較強(qiáng), 強(qiáng)烈的再懸浮作用可以使海底表層沉積物釋放含有高濃度類腐殖質(zhì)熒光組分的間隙水, 增加水體中類腐殖質(zhì)含量[39]。(355)與可溶性有機(jī)碳濃度同樣具有顯著的正相關(guān)關(guān)系(<0.01)。由于同等可溶性有機(jī)碳濃度下, 渾濁帶水體(包括邊緣站位8)具有更高的(355)值, 因此, 在不考慮最大渾濁帶的情況下, 發(fā)現(xiàn)可溶性有機(jī)碳與(355)線性擬合程度更高(圖3b)。渾濁帶區(qū)水體具有更高(355)值的一個(gè)主要原因是渾濁帶強(qiáng)的混合作用使海底表層沉積物再懸浮釋放有色溶解有機(jī)物至表層水中, 從而造成水體有色溶解有機(jī)物濃度的升高[40-41]。另一方面, 有研究揭示, 在渾濁帶區(qū)顆粒態(tài)有機(jī)物(POM)可以通過解吸作用釋放可溶性有機(jī)質(zhì)[39],同時(shí), 顆粒態(tài)有機(jī)物和絮凝體的光降解作用也可以產(chǎn)生有色溶解有機(jī)物[42-44], 最終導(dǎo)致渾濁帶水體有色溶解有機(jī)物含量水平異常升高。

    圖2 YE斷面表層鹽度(a)、表層可溶性有機(jī)碳濃度(b)、表底層濁度(c、d)平面分布圖

    圖3 吸收系數(shù)a(355)與鹽度(a)、可溶性有機(jī)碳濃度(b)的關(guān)系圖(站位8未計(jì)入統(tǒng)計(jì))

    2.2 水體可溶有機(jī)質(zhì)有色溶解有機(jī)物的熒光組分特征、空間展布及其生源指示

    利用PARAFAC鑒定出長(zhǎng)江口表層水樣品熒光可溶性有機(jī)質(zhì)(FDOM)的4個(gè)熒光組分C1、C2、C3和C4(圖4, 表1)。通過與文獻(xiàn)比對(duì), 并根據(jù)長(zhǎng)江口水體實(shí)際背景, 確定C1為類色氨酸, 主要代表原地生產(chǎn)力, 該熒光信號(hào)在海洋樣品中往往占主導(dǎo)地位[45-47]; C2和C3代表陸源類腐殖質(zhì)生源, 其中C2主要與高營養(yǎng)水體或廢水有關(guān)[45,48,49]。而C3的來源則比較廣泛, 在生物活動(dòng)較強(qiáng)的環(huán)境中較普遍[50], 如酵母細(xì)菌能夠產(chǎn)生與它相似的組分[51], 廢水中也發(fā)現(xiàn)過它的存在[52]; C4為類絡(luò)氨酸, 其來源廣泛, 可以是生物活動(dòng), 陸源輸入, 也可以是人類活動(dòng)產(chǎn)生[49]。長(zhǎng)江口位于繁華的上海市, 大量的工業(yè)、生活廢水都有可能成為長(zhǎng)江口水體有機(jī)質(zhì)的潛在源, 盡管有嚴(yán)格的管理存在。已經(jīng)有研究人類活動(dòng)所產(chǎn)生的類蛋白質(zhì)組分是長(zhǎng)江口水體陸源可溶性有機(jī)碳的重要組成部分之一[28]。

    圖4 PARAFAC鑒別出的長(zhǎng)江口4個(gè)熒光有機(jī)質(zhì)熒光組分

    表1 PARAFAC鑒別出的長(zhǎng)江口4個(gè)熒光有機(jī)質(zhì)熒光組分及其特征

    注: T-類蛋白質(zhì)(色氨酸)峰; A-類腐殖質(zhì)峰; M-類腐殖質(zhì)峰; B-類蛋白質(zhì)峰(絡(luò)氨酸)。

    由圖5a可見, 由河口, 經(jīng)最大渾濁帶至外海, 表層水體熒光有機(jī)質(zhì)熒光組分呈現(xiàn)規(guī)律性變化。代表原地生產(chǎn)力的熒光組分C1絕對(duì)含量表現(xiàn)為先穩(wěn)定-快速上升-快速下降-再穩(wěn)定的過程??焖偕仙涂焖傧陆颠^程主要與河口最大渾濁帶的水動(dòng)力過程有關(guān)。河流入海處由于強(qiáng)的水動(dòng)力條件, 雖然富含營養(yǎng)鹽, 但不利于藻類等水生生物繁育, 原地生產(chǎn)力基本穩(wěn)定[53-56]。遠(yuǎn)離河口又未至最大渾濁帶區(qū)域內(nèi)水動(dòng)力減弱(站位5、7、8和9), 河流帶來豐富的營養(yǎng)物質(zhì), 促使初級(jí)生產(chǎn)力快速提升, 在物理-生物化學(xué)作用下釋放大量的可溶性有機(jī)質(zhì)[57]。至最大渾濁帶, 由于光的限制降低了初級(jí)生產(chǎn)力, 使得補(bǔ)給速率小于移除速率[58-60], 其濃度出現(xiàn)斷崖式下降。到Y(jié)E-17站位以東水體環(huán)境再次處于平衡, 生產(chǎn)力漸穩(wěn)定。這種由于水動(dòng)力條件改變引起的生產(chǎn)力動(dòng)態(tài)過程較好的在熒光組分C1上得到了體現(xiàn), 同時(shí)也體現(xiàn)在可溶性有機(jī)碳的碳同位素組成上。與站位1、2和4相比, 站位5、7和8雖然鹽度均接近于0, 但由于藻類貢獻(xiàn)增加, 可溶性有機(jī)碳顯著富集13C(圖6), 這是由于與陸源有機(jī)質(zhì)相比, 海洋藻類生源有機(jī)質(zhì)往往富集13C所致[3,61-63]。代表陸源類腐殖質(zhì)生源的C2和C3組分則在河口與最大渾濁帶的緩沖區(qū)稍有富集, 然后呈現(xiàn)急劇下降趨勢(shì), 至外海達(dá)到一個(gè)相對(duì)穩(wěn)定的含量水平, 這也證實(shí)了海水稀釋的調(diào)控作用。熒光組分C4沿整個(gè)斷面并未發(fā)生顯著變化, 可能來自于難降解的海洋自生源[48-49], 這一點(diǎn)在熒光組分相對(duì)強(qiáng)度中得到了更好的證實(shí)(圖5b)。與絕對(duì)含量不同, 各組分相對(duì)含量沿?cái)嗝嬷镣夂? 典型地表現(xiàn)為表征陸源類腐殖質(zhì)生源的C2和C3組分相對(duì)含量減少, 而表征海洋藻類生產(chǎn)力的C1組分及其難降解組分的C4, 尤其以C4組分顯著富集為特征, C1組分只是略有富集。這說明在由陸到海的過程中, C4組分和表征海洋藻類生產(chǎn)力的C1組分在有色溶解有機(jī)物中占的比例越來越多, 而陸源信號(hào)則呈現(xiàn)出逐漸降低的趨勢(shì),13C值也同樣呈現(xiàn)出富13C的海源信號(hào)(圖6)。

    圖5 熒光組分絕對(duì)強(qiáng)度(a)和相對(duì)強(qiáng)度(b)站位分布圖

    圖6 固相萃取所得可溶性有機(jī)質(zhì)δ13C和可溶性有機(jī)碳濃度隨站位分布圖

    通過光譜斜率275-295和350-400的比值獲得的斜率比S通常用來區(qū)分有色溶解有機(jī)物的來源。研究表明, 陸源有色溶解有機(jī)物平均相對(duì)分子質(zhì)量較大, 所以S值較低, 而海源有色溶解有機(jī)物擁有較低的平均相對(duì)分子質(zhì)量,S值較高[29]。Stedmon.[64]研究揭示, 在保守混合情況下, 斜率比S與鹽度存在指數(shù)關(guān)系。由圖7所示,S在長(zhǎng)江口的分布從近岸到遠(yuǎn)岸, 沿著保守混合趨勢(shì)逐漸增大(圖7), 說明近岸陸源輸入對(duì)于S的值具有調(diào)控作用, 遠(yuǎn)岸陸源輸入勢(shì)力減弱, 微生物降解能力增強(qiáng), 陸源帶來的腐殖化程度較高的有機(jī)質(zhì)減少, 引起S值增大。值得一提的是, 河口至最大渾濁帶區(qū)域S值存在脫離于保守混合線的顯著增高過程, 而后回落至水舌邊緣再次快速升高, 與熒光組分C1的動(dòng)態(tài)過程相一致, 再次證實(shí)了藻類生產(chǎn)力沿?cái)嗝娴膭?dòng)態(tài)過程。相應(yīng)地, 可溶性有機(jī)質(zhì)碳同位素值由陸到??傮w呈富13C趨勢(shì)(圖6), 說明海洋自生源產(chǎn)生的可溶性有機(jī)碳越來越占主導(dǎo)地位。

    2.3 水體可溶有機(jī)質(zhì)在長(zhǎng)江口的混合行為

    普遍認(rèn)為, 在不考慮光降解、溶解態(tài)與顆粒態(tài)相互轉(zhuǎn)化、其他影響可溶性有機(jī)質(zhì)性質(zhì)和濃度變化的情況下, 河流和海水兩端元混合后水體可溶性有機(jī)碳與鹽度呈線性關(guān)系, 屬于保守性行為[5]。本研究中1站位鹽度為0.12, 可以作為河流端元, 19站位鹽度為33.63, 是此次采樣站位獲得的最高鹽度, 作為海水端元值, 根據(jù)每個(gè)站位的鹽度值計(jì)算出咸、淡水的混合比例[4], 可以生成可溶性有機(jī)碳保守混合理論線(圖8)。實(shí)際樣品可溶性有機(jī)碳濃度若高于該線, 表示存在可溶性有機(jī)碳的添加過程, 反之, 若低于該線, 則反映可溶性有機(jī)碳的移除過程。由圖8可見, 在最大渾濁帶核心區(qū)(9、10和11), 可溶性有機(jī)碳濃度明顯低于理論線, 可能指示在兩端元的混合過程中, 部分可溶性有機(jī)碳被選擇性移除。相對(duì)而言, 其他區(qū)域水體可溶性有機(jī)碳保守性比較強(qiáng)。雖然可溶性有機(jī)碳與鹽度的關(guān)系, 可以大致說明可溶性有機(jī)碳的行為保守與否, 但是這種關(guān)系只能說明可溶性有機(jī)碳量的變化, 無法從源上區(qū)分可溶性有機(jī)碳的變化。因此, 本研究依據(jù)同位素質(zhì)量守恒的原理, 選用Raymond.[65]提出的雙端元保守混合模型, 通過碳同位素值的變化來辨別可溶性有機(jī)碳的來源。13C保守混合值的計(jì)算公式如下。

    圖7 斜率比SR與鹽度的關(guān)系

    式中,代表河流混合比例; (1?)代表海水混合比例;R代表河流可溶性有機(jī)碳濃度(以1站位可溶性有機(jī)碳濃度代替);M代表海水可溶性有機(jī)碳濃度(以19站位可溶性有機(jī)碳濃度代替);mix代表河水和海水不同比例混合后碳同位素的計(jì)算結(jié)果,r代表河流末端成員同位素組成;m代表海洋末端成員同位素組成;代表保守混合時(shí)可溶性有機(jī)碳的理論值。

    由圖8可見, 雖然在河流區(qū)鹽度變化微弱, 但是由于光-微生物的降解作用[11], 可溶性有機(jī)碳的濃度有所降低, 并在河流區(qū)和最大渾濁帶的分界處, 由于咸淡水的混合稀釋作用, 出現(xiàn)斷崖式的降低趨勢(shì)。同時(shí), 由于水動(dòng)力條件控制的生源分布情況, 初級(jí)生產(chǎn)力的貢獻(xiàn)使得13C在5、7和8站位富集(上文已討論)。及至最大渾濁帶, 隨著鹽度的增加,13C值逐漸升高, 海源可溶性有機(jī)碳逐漸占據(jù)主導(dǎo)地位。值得注意的是, 在最大渾濁帶實(shí)際碳同位素值與理論曲線吻合度較高, 說明在最大渾濁帶區(qū), 可溶性有機(jī)碳濃度的下降以物理過程為主, 由海水稀釋和可溶性有機(jī)碳絮凝沉降過程共同支配。

    圖8 可溶性有機(jī)碳濃度和固相萃取有機(jī)質(zhì)δ13C隨鹽度分布圖

    3 結(jié) 論

    (1) 水動(dòng)力條件是控制長(zhǎng)江口-東海陸架系統(tǒng)可溶性有機(jī)碳分布的主要因素之一。在河口之外的海域, 可溶性有機(jī)碳濃度隨著采樣站位的離岸距離呈現(xiàn)明顯降低的趨勢(shì), 而碳同位素值表現(xiàn)為富13C趨勢(shì), 指示陸源可溶性有機(jī)碳相對(duì)貢獻(xiàn)減少, 海源可溶性有機(jī)碳貢獻(xiàn)增加,S值從陸到海逐漸升高, 反映可溶性有機(jī)質(zhì)生源由陸源向海源藻類的演變。

    (2) 在最大渾濁帶區(qū), 河流-海水兩端元橫向混合劇烈, 可溶性有機(jī)碳濃度明顯低于理論值, 說明最大渾濁帶對(duì)可溶性有機(jī)碳有去除作用。從熒光組分的分布特征來看, 陸源類腐殖質(zhì)組分C2、C3和類蛋白質(zhì)組分C1在該區(qū)域均有減小的現(xiàn)象, 并與鹽度呈現(xiàn)較好的保守關(guān)系, 說明陸源有機(jī)質(zhì)的減少以海水稀釋的物理過程為主, 并得到碳同位素?cái)?shù)據(jù)的支持。

    (3) 長(zhǎng)江口水體可溶性有機(jī)質(zhì)性質(zhì)主要受控于水體鹽度、有機(jī)質(zhì)來源及其所經(jīng)歷的降解程度, 河口區(qū)、最大渾濁帶、沖淡水邊緣帶和海域等4個(gè)區(qū)域水體內(nèi)可溶性有機(jī)質(zhì)性質(zhì)存在陸源和海源的差異, 并在最大渾濁帶的過濾作用下, 可溶性有機(jī)碳濃度出現(xiàn)明顯的降低趨勢(shì), 說明水動(dòng)力條件是控制長(zhǎng)江口-東海陸架系統(tǒng)可溶性有機(jī)質(zhì)組成和空間分布的重要控制因素。

    本次YES航次得到了同濟(jì)大學(xué)海洋地質(zhì)國家重點(diǎn)實(shí)驗(yàn)室的經(jīng)費(fèi)支持, 在此表示感謝。同濟(jì)大學(xué)的李超博士和南方科技大學(xué)的李朋輝博士分別在采樣過程中以及后期三維熒光光譜數(shù)據(jù)處理中提供了幫助和指導(dǎo), 在此一并誠謝!匿名專家和胡建芳研究員為此文的發(fā)表提出了非常寶貴的意見, 在此表示感謝!

    [1] Hedges J I. Global biogeochemical cycles: Progress and problems[J]. Mar Chem, 1992, 39(1): 67-93.

    [2] Bade D L, Carpenter S R, Cole J J, Pace M L, Kritzberg E, Van de Bogert M C, Cory R M, McKnight D M. Sources and fates of dissolved organic carbon in lakes as determined by whole-lake carbon isotope additions[J]. Biogeochemistry, 2007, 84(2): 115-129.

    [3] Hedges J I. What happens to terrestrial organic matter in the ocean? [J]. Org Geochem, 1997, 27(5/6): 195-212.

    [4] Liss P S. Conservative and non-conservative behavior of dissolved constituents during estuarine mixing[M]//Burton J D, Liss P S. Estuarine Chemistry. London: Academic Press, 1976: 93-130.

    [5] Wen L, Shiller A, Santschi P H, Gill G. Trace element behavior in Gulf of Mexico estuaries[M]//Bianchi T S, Pennock J R, Twilley R R. Biogeochemistry of Gulf of Mexico Estuaries. New York: Wiley, 1999: 303-346.

    [6] 王江濤, 于志剛, 張經(jīng). 鴨綠江口溶解有機(jī)碳的研究[J]. 中國海洋大學(xué)學(xué)報(bào)(自然科學(xué)版), 1998, 40(3): 471-475.

    Wang Jiang-tao, Yu Zhi-gang, Zhang Jing. Study on dissolved organic carbon in the Estuary of the Yalujiang River[J]. Period Ocean Univ China (Nat Sci Ed), 1998, 40(3): 471-475 (in Chinese with English abstract).

    [7] Benner R, Opsahl S. Molecular indicators of the sources and transformations of dissolved organic matter in the Mississippi river plume[J]. Org Geochem, 2001, 32(4): 597-611.

    [8] Abril G, Nogueira M, Etcheber H, Cabecadas G, Lemaire E, Brogueira M J. Behaviour of organic carbon in nine contrastingeuropean estuaries[J]. Estuarine Coastal Shelf Sci, 2002, 54(2): 241-262.

    [9] 張向上. 黃河口有機(jī)碳的時(shí)空分布及影響因素研究[D]. 青島: 中國海洋大學(xué), 2004.

    Zhang Xiang-shang. Temporal and spatial distribution of organic carbon and its influence aspects In Huanghe Estuary[D]. Qingdao: Ocean University of China, 2004 (in Chinese with English abstract).

    [10] Yang F, Huang Q H, Li J H, Zhu X M. Chracterization of chromophoric dissolved organic matter in the Yangtze Estuary by absoption and fluorescence spectroscopy[J]. J Environ Sci Sustain Soc, 2007, 1(1): 55-60.

    [11] Medeiros P M, Seidel M, Ward N D, Carpenter E J, Gomes H R, Niggemann J, Krusche A V, Yager P L, Dittmar T. Fate of the Amazon River dissolved organic matter in the tropical Atlantic Ocean[J]. Global Biogeochem Cycle, 2015, 29(5): 677-690.

    [12] álvarez-Salgado X A, Miller A E J. Dissolved organic carbon in a large macrotidal estuary (the Humber, UK): Behaviour during estuarine mixing[J]. Mar Pollut Bull, 1999, 37(3/7): 216-224.

    [13] Cauwet G. Carbon inputs and biogeochemical processes at the halocline in a stratified estuary: Krka River, Yugoslavia[J]. Mar Chem, 1991, 32(2/4): 269-283.

    [14] Callahan J, Dai M H, Chen R F, Li X L, Lu Z M, Huang W. Distribution of dissolved organic matter in the Pearl River Estuary, China[J]. Mar Chem, 2004, 89(1/4): 211-224.

    [15] Zhang J, Liu S M, Ren J L, Wu Y, Zhang G L. Nutrient gradients from the eutrophic Changjiang (Yangtze river) Estuary to the oligotrophic Kuroshio waters and re-evaluation of budgets for the East China Sea Shelf[J]. Prog Oceanogr, 2007, 74(4): 449-478.

    [16] 沈煥庭, 賀松林, 潘定安, 李九發(fā). 長(zhǎng)江河口最大渾濁帶研究[J]. 地理學(xué)報(bào), 1992, 47(5): 472-479. Shen Huan-ting, He Song-lin, Pan Ding-an, Li Jiu-fa. A study of turbidity maximum in the Changjiang Estuary[J]. Geogr J, 1992, 47(5): 472-479 (in Chinese with English abstract).

    [17] Ogawa H, Usui T, Koike I. Distribution of dissolved organic carbon in the East China Sea[J]. Deep Sea Res Pt II Topic Stud Oceanogr, 2003, 50(2): 353-366.

    [18] 朱偉健, 沈芳, 洪官林. 長(zhǎng)江口及鄰近海域有色溶解有機(jī)物(CDOM)的光學(xué)特性[J]. 環(huán)境科學(xué), 2010, 31(10): 2292- 2298.

    Zhu Wei-jian, Shen Fang, Hong Guan-lin. Optical characteristics of colored dissolved organic material (CDOM) in Yangtze Estuary[J]. Environ Sci, 2010, 31(10): 2292-2298 (in Chinese with English abstract).

    [19] 范冠南, 毛志華, 陳鵬, 王天愚, 張琳. 長(zhǎng)江口及其鄰近海域CDOM光譜吸收特性分析[J]. 海洋學(xué)研究, 2013, 31(1): 53-58.

    Fan Guan-nan, Mao Zhi-hua, Chen Peng, Wang Tian-yu, Zhang lin. Optical characteristics of colored dissolved organic matter (CDOM) in Changjiang River Estuary[J]. J Mar Sci, 2013, 31(1): 53-58 (in Chinese with English abstract).

    [20] 呂麗莎, 趙衛(wèi)紅, 苗輝. 三維熒光結(jié)合平行因子分析在東海溶解有機(jī)物研究中的應(yīng)用[J]. 光譜學(xué)與光譜分析, 2013, 33(3): 653-658.

    Lü Li-sha, Zhao Wei-hong, Miao hui. Application of excitation- emission matrix spectrum combined with parallel factor analysis in dissolved organic matter in East China Sea[J]. Spect Spect Anal, 2013, 33(3): 653-658 (in Chinese with English abstract).

    [21] Sun Q Y, Wang C, Wang P F, Hou J, Ao Y H. Absorption and fluorescence characteristics of chromophoric dissolved organic matter in the Yangtze Estuary[J]. Environ Sci Pollut Res Int, 2014, 21(5): 3460-3473.

    [22] Li P H, Chen L, Zhang W, Huang Q H. Spatiotemporal distribution, sources, and photobleaching imprint of dissolved organic matter in the Yangtze Estuary and its adjacent sea using fluorescence and parallel factor analysis[J]. Plos One, 2015, 10(6): e0130852.

    [23] Guo W D, Stedmon C A, Han Y C, Wu F, Yu X X, Hu M H. The conservative and non-conservative behavior of chromophoric dissolved organic matter in Chinese estuarine waters[J]. Mar Chem, 2007, 107(3): 357-366.

    [24] 徐亞宏, 姚鵬, 蘇榮國, 趙彬. 長(zhǎng)江口鹽度梯度下有色溶解有機(jī)物的分布、來源與季節(jié)變化[J]. 海洋學(xué)報(bào), 2017, 39(6): 21-32.

    Xu Ya-hong, Yao Peng, Su Rong-guo, Zhao Bin. Distribution, sources, and seasonal variation of colored dissolved organic matter along salinity gradients in the Changjiang Estuary[J]. Acta Oceanol Sinica, 2017, 39(6): 21-32 (in Chinese with English abstract).

    [25] Guo W D, Yang L Y, Zhai W D, Chen W Z, Osburn C L, Huang X, Li Y. Runoff‐mediated seasonal oscillation in the dynamics of dissolved organic matter in different branches of a large bifurcated estuary—The Changjiang Estuary[J]. J Geophys Res Biogeosci, 2014, 119(5): 776-793.

    [26] Dittmar T, Koch B P, Hertkorn N, Kattner G. A simple and efficient method for the solid‐phase extraction of dissolved organic matter (SPE-DOM) from seawater[J]. Limnol Oceanogr Method, 2008, 6(6): 230-235.

    [27] Dong Q Q, Li P H, Huang Q H, Abdelhafez A A, Chen L. Occurrence, polarity and bioavailability of dissolved organic matter in the Huangpu River, China[J]. J Environ Sci, 2014, 26(9): 1843-1850.

    [28] Twardowski M S, Boss E, Sullivan J M, Donaghay P L. Modeling the spectral shape of absorption by chromophoric dissolved organic matter[J]. Mar Chem, 2004, 89(1/4): 69-88.

    [29] Helms J R, Stubbins A, Ritchie J D, Minor E C. Absorption spectral slopes and slope ratios as indicators of molecular weight, source, and photobleaching of chromophoric dissolved organic matter[J]. Limnol Oceanogr, 2008, 53(3): 955-969.

    [30] Bai Y, Su R G, Shi X Y. Assessing the dynamics of chromophoric dissolved organic matter in the southern Yellow Sea by excitation-emission matrix fluorescence and parallel factor analysis (EEM-PARAFAC) [J]. Continent Shelf Res, 2014, 88: 103-116.

    [31] Stedmon C A, Markager S, Bro R. Tracing dissolved organic matter in aquatic environments using a new approach to fluorescence spectroscopy[J]. Mar Chem, 2003, 82(3/4): 239-254.

    [32] 黃昌春, 李云梅, 王橋, 施坤, 金鑫, 王彥飛. 基于三維熒光和平行因子分析法的太湖水體CDOM組分光學(xué)特征[J]. 湖泊科學(xué), 2010, 22(3): 375-382.

    Huang Chang-chun, Li Yun-mei, Wang Qiao, Shi Kun, Jin Xin, Wang Yan-fei. Components optical property of CDOM in Lake Taihu based on three-dimensional excitation emission matrix fluorescence[J]. J Lake Sci, 2010, 22(3): 375-382 (in Chinese with English abstract).

    [33] 鄭琰明, 康建成, 鐘巍, 周尚哲, 王芳, 徐韌, 孫瑞文, 吳濤. 2004年夏初、秋末長(zhǎng)江口外海區(qū)沖淡水及羽狀鋒的鹽度特征[J]. 海洋與湖沼, 2009, 40(3): 275-282.

    Zheng Yan-ming, Kang Jian-cheng, Zhong Wei, Zhou Shang-zhe, Wang Fang, Xu Ren, Sun Rui-wen, Wu Tao. Salinity distribution of diluted water and plume front of early summer & later autumn in 2004 off Changjiang (Yangtze) Estuary[J]. Oceanol Limnol Sinica, 2009, 40(3): 275-282 (in Chinese with English abstract).

    [34] 王凱敏, 熊學(xué)軍, 郭炳火, 華鋒, 江志輝. 2006-2007年長(zhǎng)江沖淡水的擴(kuò)展形態(tài)及季節(jié)變化[J]. 海岸工程, 2012, 31(1): 46-54.

    Wang Kai-min, Xiong Xue-jun, Guo Bing-huo, Hua Feng, Jiang Zhi-hui. The extension form and seasonal variation of the Changjiang diluted water during 2006-2007[J]. Coast Eng, 2012, 31(1): 46-54 (in Chinese with English abstract).

    [35] Hong H S, Wu J Y, Shang S L, Hu C M. Absorption and fluorescence of chromophoric dissolved organic matter in the Pearl River Estuary, South China[J]. Mar Chem, 2005, 97(1/2): 78-89.

    [36] Vecchio R D, Blough N V. Spatial and seasonal distribution of chromophoric dissolved organic matter and dissolved organic carbon in the Middle Atlantic Bight[J]. Mar Chem, 2004, 89(1/4): 169-187.

    [37] Rochelle-Newall E J, Fisher T R. Chromophoric dissolved organic matter and dissolved organic carbon in Chesapeake Bay[J]. Mar Chem, 2002, 77(1): 23-41.

    [38] 潘慧慧, 姚鵬, 趙彬, 孟佳, 李棟, 王金鵬. 基于水淘選分級(jí)的長(zhǎng)江口最大渾濁帶附近顆粒有機(jī)碳的來源、分布和保存[J]. 海洋學(xué)報(bào), 2015, 37(4): 1-15.

    Pan Hui-hui, Yao Peng, Zhao Bin, Meng Jia, Li Dong, Wang Jin-peng. Sources, distribution and preservation of size-fractionated particulate organic carbon in the turbidity maximum zone of the Changjiang Estuary based on water elutriation[J]. Acta Oceanol Sinia, 2015, 37(4): 1-15 (in Chinese with English abstract).

    [39] Markager S, Vincent W F. Spectral light attenuation and the absorption of UV and blue light in natural waters[J]. Limnol Oceanogr, 2000, 45(3): 642-650.

    [40] Yang L Y, Guo W D, Hong H S, Wang G Z. Non-conservative behaviors of chromophoric dissolved organic matter in a turbid estuary: Roles of multiple biogeochemical processes[J]. Estuarine Coastal Shelf Sci, 2013, 133(4): 285-292.

    [41] Bodineau L, Thoumelin G, Béghin V, Wartel M. Particulate organic matter composition in the Estuarine Turbidity Maxima (ETM) of the Seine River estuary[J]. Hydrobiologia, 1998, 373: 281-295.

    [42] Wang Y, Zhang D, Shen Z Y, Feng C H, Chen J. Revealing sources and distribution changes of dissolved organic matter (DOM) in pore water of sediment from the Yangtze estuary[J]. Plos One, 2013, 8(10): e76633.

    [43] Kieber R J, Whitehead R F, Reid S N, Willey J D, Seaton P J. Chromophoric dissolved organic matter (CDOM) in rainwater, southeastern North Carolina, USA[J]. J Atmos Chem, 2006, 54(1): 21-41.

    [44] Mayer L M, Schick L L, Skorko K, Boss E. Photodissolution of particulate organic matter from sediments[J]. Limnol Oceanogr, 2006, 51(2): 1064-1071.

    [45] Helms J R, Mao J D, Stubbins A, Schmidt-Rohr K, Spencer R G M, Hernes P J, Mopper K. Loss of optical and molecular indicators of terrigenous dissolved organic matter during long-term photobleaching[J]. Aquat Sci, 2014, 76(3): 353-373.

    [46] Yu H R, Liang H, Qu F S, Han Z S, Shao S L, Chang H Q, Li G B. Impact of dataset diversity on accuracy and sensitivity of parallel factor analysis model of dissolved organic matter fluorescence excitation-emission matrix[J]. Sci Rep, 2015, 5: 10207.

    [47] Kothawala D N, Stedmon C A, Müller R A, Weyhenmeyer G A, K?hler S J, Tranvik L J. Controls of dissolved organic matter quality: Evidence from a large-scale boreal lake survey[J]. Global Chang Biol, 2014, 20(4): 1101-1114.

    [48] Osburn C L, Stedmon C A. Linking the chemical and optical properties of dissolved organic matter in the Baltic-North Sea transition zone to differentiate three allochthonous inputs[J]. Mar Chem, 2011, 126(1/4): 281-294.

    [49] Murphy K R, Hambly A, Singh S, Henderson R K, Baker A, Stuetz R, Khan S J. Organic matter fluorescence in municipal water recycling schemes: Toward a unified PARAFAC model[J]. Environ Sci Technol, 2011, 45(7): 2909-2916.

    [50] J?rgensen L, Stedmon C A, Kragh T, Markager S, Middelboe M, S?ndergaard M. Global trends in the fluorescence characteristics and distribution of marine dissolved organic matter[J]. Mar Chem, 2011, 126(1/4): 139-148.

    [51] Nimptsch J, Woelfl S, Osorio S, Valenzuela J, Ebersbach P, von Tuempling W, Palma R, Encina F, Figueroa D, Kamjunke N, Graeber D. Tracing dissolved organic matter (DOM) from land-based aquaculture systems in North Patagonian streams[J]. Sci Total Environ, 2015, 537: 129-138.

    [52] Marose S, Lindemann C, Scheper T. Two-dimensional fluorescence spectroscopy: A new tool for on-line bioprocess monitoring[J]. Biotechnol Prog, 2010, 14(1): 63-74.

    [53] Stedmon C A, Markager S. Resolving the variability in dissolved organic matter fluorescence in a temperate estuary and its catchment using PARAFAC analysis[J]. Limnol Oceanogr, 2005, 50(2): 686-697.

    [54] 顧新根, 袁騏, 沈煥庭, 周月琴. 長(zhǎng)江口最大渾濁帶浮游植物的生態(tài)研究[J]. 中國水產(chǎn)科學(xué), 1995, 2(1): 16-27.

    Gu Xin-gen, Yuan Qi, Shen Huan-ting, Zhou Yue-qin. The ecological study on phytoplankton in maximum turbid zone of Changjiang Estuary[J]. J Fish Sci China, 1995, 2(1): 16-27(in Chinese with English abstract).

    [55] Ha K, Cho E A, Kim H W, Joo G J. Microcystis bloom formation in the lower Nakdong River, South Korea: Importance of hydrodynamics and nutrient loading[J]. Mar Freshwater Res, 1999, 50(1): 89-94.

    [56] Mitrovic S M, Chessman B C, Davie A, Avery E L, Ryan N. Development of blooms ofandspp. (Bacillariophyceae) in a shallow river and estimation of effective suppression flows[J]. Hydrobiologia, 2008, 596(1): 173-185.

    [57] 蔡卓平, 段璐洋, 肖群, 段舜山. 風(fēng)浪擾動(dòng)促進(jìn)中肋骨條藻和球形棕囊藻增殖的模擬研究[J]. 海洋環(huán)境科學(xué), 2011, 30(4): 473-476.

    Cai Zhuo-ping, Duan Lu-yang, Xiao Qun, Duan Shun-shan. Simulation study on promoted proliferation of Skeletonema costatum and Phaeocystis globose disturbed by wave[J]. Mar Environ Sci, 2011, 30(4): 473-476(in Chinese with English abstract).

    [58] Biddanda B, Benner R. Carbon, nitrogen, and carbohydrate fluxes during the production of particulate and dissolved organic matter by marine phytoplankton[J]. Limnol Oceanogr, 1997, 42(3): 506-518.

    [59] 沈志良. 長(zhǎng)江口海區(qū)理化環(huán)境對(duì)初級(jí)生產(chǎn)力的影響[J]. 海洋湖沼通報(bào), 1993, 15(1): 47-51.

    Shen Zhi-liang. The effects of the physic-chemical environment on the primary productivity in the Yangtze River Estuary[J]. Transact Oceanol Limnol, 1993, 15(1): 47-51(in Chinese with English abstract).

    [60] Wienke S M, Cloern J E. The phytoplankton component of seston in San Francisco Bay[J]. Netherland J Sea Res, 1987, 21(1): 25-33.

    [61] 黃小平, 黃良民. 河口最大渾濁帶浮游植物生態(tài)動(dòng)力過程研究進(jìn)展[J]. 生態(tài)學(xué)報(bào), 2002, 22(9): 1527-1533.

    Huang Xiao-ping, Huang Liang-min. Progress in researches on dynamical processes of phytoplankton ecology in maximum turbidity zone of estuary[J]. Acta Ecol Sinica, 2002, 22(9): 1527-1533(in Chinese with English abstract).

    [62] Fry B, Sherr E B.13C measurements as indicators of carbon flow in marine and freshwater ecosystems[M]. Rundel P W, Ehleringer J R, Nagy K A. Stable Isotopes in Ecological Research. Berlin: Springer, 1984: 196-229.

    [63] Lamb A L, Wilson G P, Leng M J. A review of coastal palaeoclimate and relative sea-level reconstructions using13C and C/N ratios in organic material[J]. Earth Sci Rev, 2006, 75(1/4): 29-57.

    [64] Stedmon C A, Markager S. Behaviour of the optical properties of coloured dissolved organic matter under conservative mixing[J]. Estuarine Coastal Shelf Sci, 2003, 57(5/6): 973- 979.

    [65] Raymond P A, Bauer J E. Use of14C and13C natural abundances for evaluating riverine, estuarine, and coastal DOC and POC sources and cycling: A review and synthesis[J]. Org Geochem, 2001, 32(4): 469-485.

    Spatial distribution and controlling factors of dissolved organic matter composition in the Yangtze Estuary in summer

    WANG An-yue1, ZHANG Kai1, HE Ding1, FAN Dai-du2and SUN Yong-ge1*

    1. Organic Geochemistry Unit, School of Earth Sciences, Zhejiang University, Hangzhou 310027, China; 2. State Key Laboratory of Marine Geology, Tongji University, Shanghai 200092, China

    In this study, the composition, source, spatial distribution, and mixing behavior of dissolved organic matter (DOM) in the Yangtze Estuary in summer were systematically determined using the measured absorption spectra of chromophoric dissolved organic matter (CDOM), Excitation-Emission-Matrix spectra (EEMs) of fluorescence dissolved organic matter (FDOM) and stable carbon isotope compositions. The results show that the dissolved organic carbon (DOC) concentration decreases gradually from land to sea. The absorption coefficient(355), which indicates the concentration of CDOM, is negatively correlated with salinity, indicating that the distribution of CDOM is primarily controlled by dilution. EEMs and parallel factor analysis (PARAFAC) were used to identify the two fluorescent components C2and C3representing terrestrial organic matter and two protein-like fluorescent components C1and C4representing in-situ productivity. From nearshore to offshore, the proportion of C1, which represents the productivity of marine algae, increased in fluorescent organic matter, whereas the land-source signal gradually decreased. Correspondingly, the DOM in offshore waters exhibited steep spectral slope S275-295and S350-400ratios and heavy carbon isotopes. Based on the river-ocean two-terminal mixing model of salinity vs DOC concentration and salinity.13C value, we determined that the DOC removal process occurs in the maximum turbidity zone of the Yangtze Estuary, which is characterized by a lack of non-conservative mixing behavior dominated by physical dilution. The results reveal that hydrodynamic conditions are the key factors controlling the composition and spatial distribution of DOM in the Yangtze Estuary-East China Sea shelf system.

    dissolved organic matter; carbon isotope; mixing behavior; EEMs-PARAFAC; Yangtze Estuary

    P593; P66

    A

    0379-1726(2021)03-0317-12

    10.19700/j.0379-1726.2021.03.009

    2019-10-11;

    2019-11-30;

    2019-12-02

    國家自然科學(xué)基金(41973070, 41773098)

    王安月(1994-), 女, 碩士研究生, 地球化學(xué)專業(yè)。E-mail: 1104681767@qq.com

    Corresponding author): SUN Yong-ge, E-mail: ygsun@zju.edu.cn; Tel: +86-571-87951336

    猜你喜歡
    長(zhǎng)江口站位鹽度
    提高政治站位 對(duì)標(biāo)國內(nèi)一流
    建黨百年說“站位”
    水上消防(2021年3期)2021-08-21 03:12:00
    提升站位講政治 創(chuàng)新?lián)?dāng)爭(zhēng)出彩
    長(zhǎng)江口鳳鱭繁殖群體的年齡結(jié)構(gòu)和生長(zhǎng)特性
    長(zhǎng)江口及其鄰近水域仔稚魚種類組成及分布特征
    鹽度和pH對(duì)細(xì)角螺耗氧率和排氨率的影響
    鹽度脅迫對(duì)入侵生物福壽螺的急性毒性效應(yīng)
    適用于高鹽度和致密巖層驅(qū)油的表面活性劑
    風(fēng)場(chǎng)對(duì)夏季長(zhǎng)江口外東北海域低鹽水團(tuán)的作用
    春、夏季長(zhǎng)江口及鄰近海域溶解甲烷的分布與釋放通量
    免费少妇av软件| 99香蕉大伊视频| 久久久国产一区二区| 国产免费视频播放在线视频| 亚洲精品456在线播放app| 十八禁高潮呻吟视频| 日韩,欧美,国产一区二区三区| 精品人妻偷拍中文字幕| 日韩三级伦理在线观看| 国语对白做爰xxxⅹ性视频网站| 欧美日韩亚洲高清精品| 免费看光身美女| 亚洲精品,欧美精品| 我的女老师完整版在线观看| av不卡在线播放| 97在线人人人人妻| 性高湖久久久久久久久免费观看| 日本色播在线视频| 在线观看美女被高潮喷水网站| 亚洲精品第二区| 日韩精品免费视频一区二区三区 | 色视频在线一区二区三区| 少妇 在线观看| 丝袜喷水一区| 亚洲第一区二区三区不卡| av片东京热男人的天堂| 日本欧美国产在线视频| 国产色婷婷99| 男女免费视频国产| 香蕉国产在线看| 日韩制服骚丝袜av| 桃花免费在线播放| av又黄又爽大尺度在线免费看| 久久久国产一区二区| 亚洲欧美日韩卡通动漫| 色婷婷久久久亚洲欧美| 欧美 日韩 精品 国产| 青春草视频在线免费观看| 欧美老熟妇乱子伦牲交| 插逼视频在线观看| 国产黄色免费在线视频| 亚洲国产看品久久| 午夜福利视频在线观看免费| 日日啪夜夜爽| 精品亚洲乱码少妇综合久久| 久久久久久人妻| 精品一区二区三卡| 人妻一区二区av| 免费在线观看完整版高清| a级片在线免费高清观看视频| 蜜桃在线观看..| 久久精品aⅴ一区二区三区四区 | 亚洲,一卡二卡三卡| 色视频在线一区二区三区| 美国免费a级毛片| 99久久精品国产国产毛片| 女性被躁到高潮视频| 青春草国产在线视频| av线在线观看网站| 男女无遮挡免费网站观看| 熟女av电影| 最近中文字幕2019免费版| 亚洲精品视频女| 免费在线观看黄色视频的| 夫妻午夜视频| 卡戴珊不雅视频在线播放| a级毛片黄视频| 免费看av在线观看网站| 国产亚洲av片在线观看秒播厂| 国产精品人妻久久久久久| 波野结衣二区三区在线| 寂寞人妻少妇视频99o| 日韩中文字幕视频在线看片| 51国产日韩欧美| 少妇被粗大的猛进出69影院 | 人人妻人人澡人人看| 欧美国产精品一级二级三级| 久久国内精品自在自线图片| 国产成人午夜福利电影在线观看| 久久99蜜桃精品久久| 久久人人爽av亚洲精品天堂| av线在线观看网站| 母亲3免费完整高清在线观看 | 国精品久久久久久国模美| 亚洲精品国产av蜜桃| 少妇高潮的动态图| 黑人高潮一二区| 色婷婷久久久亚洲欧美| 丰满饥渴人妻一区二区三| 亚洲欧美中文字幕日韩二区| 国产麻豆69| 色吧在线观看| 九色成人免费人妻av| 欧美日韩视频高清一区二区三区二| 欧美人与性动交α欧美软件 | 免费女性裸体啪啪无遮挡网站| 男人操女人黄网站| 国产精品三级大全| 一级毛片我不卡| 亚洲精品中文字幕在线视频| 亚洲三级黄色毛片| 亚洲成色77777| 大香蕉久久成人网| 国产男人的电影天堂91| 啦啦啦视频在线资源免费观看| av不卡在线播放| 久久人人97超碰香蕉20202| 午夜免费男女啪啪视频观看| 成年美女黄网站色视频大全免费| 欧美人与性动交α欧美精品济南到 | 久久久久久久久久久久大奶| 日本爱情动作片www.在线观看| 中文字幕免费在线视频6| 99精国产麻豆久久婷婷| 久久ye,这里只有精品| 国产一级毛片在线| 肉色欧美久久久久久久蜜桃| 久久人人爽人人爽人人片va| 女人被躁到高潮嗷嗷叫费观| 精品国产露脸久久av麻豆| 日韩中文字幕视频在线看片| 亚洲精品456在线播放app| 国产成人精品在线电影| 国产精品国产三级国产av玫瑰| 五月天丁香电影| 18+在线观看网站| 成人黄色视频免费在线看| 91精品伊人久久大香线蕉| 国产亚洲av片在线观看秒播厂| 天天操日日干夜夜撸| 搡老乐熟女国产| 国产一级毛片在线| 满18在线观看网站| 啦啦啦视频在线资源免费观看| 91成人精品电影| 男女下面插进去视频免费观看 | 女性生殖器流出的白浆| 欧美国产精品一级二级三级| 天天躁夜夜躁狠狠躁躁| 久久免费观看电影| 一二三四中文在线观看免费高清| 青春草视频在线免费观看| 一区二区三区精品91| 热re99久久国产66热| av在线老鸭窝| 亚洲精品乱码久久久久久按摩| 亚洲精品久久午夜乱码| 成年女人在线观看亚洲视频| 一区在线观看完整版| 亚洲高清免费不卡视频| 18禁国产床啪视频网站| 最近最新中文字幕免费大全7| 亚洲国产精品成人久久小说| 少妇高潮的动态图| 夫妻午夜视频| 久久久久久人人人人人| 男女免费视频国产| 久久综合国产亚洲精品| 日韩视频在线欧美| 国产欧美日韩综合在线一区二区| 曰老女人黄片| 看免费成人av毛片| 少妇的逼好多水| 日本av免费视频播放| 国产免费视频播放在线视频| 欧美变态另类bdsm刘玥| www.色视频.com| 美女国产高潮福利片在线看| 国产白丝娇喘喷水9色精品| 一本—道久久a久久精品蜜桃钙片| 最近2019中文字幕mv第一页| 久久精品人人爽人人爽视色| 午夜激情久久久久久久| 欧美成人精品欧美一级黄| 看十八女毛片水多多多| 精品午夜福利在线看| 少妇的逼好多水| 国产成人精品久久久久久| 亚洲第一区二区三区不卡| 日韩av在线免费看完整版不卡| 国产精品人妻久久久影院| 性色avwww在线观看| 亚洲精品国产av蜜桃| 不卡视频在线观看欧美| 亚洲美女黄色视频免费看| 欧美少妇被猛烈插入视频| 女性被躁到高潮视频| 久久国产精品男人的天堂亚洲 | 丰满乱子伦码专区| 国产精品久久久久久精品电影小说| tube8黄色片| 免费女性裸体啪啪无遮挡网站| 国产一区亚洲一区在线观看| 婷婷成人精品国产| 成人毛片a级毛片在线播放| 久久97久久精品| 美女中出高潮动态图| 国产熟女午夜一区二区三区| 99久久精品国产国产毛片| 欧美老熟妇乱子伦牲交| 午夜激情久久久久久久| 亚洲国产精品国产精品| 亚洲美女黄色视频免费看| 一本大道久久a久久精品| 亚洲国产精品999| 美女国产视频在线观看| 国产精品久久久久久精品古装| 全区人妻精品视频| 中文字幕av电影在线播放| 日韩视频在线欧美| 亚洲av.av天堂| 欧美成人午夜精品| 亚洲av日韩在线播放| 青青草视频在线视频观看| 丝袜喷水一区| 看非洲黑人一级黄片| 亚洲欧美中文字幕日韩二区| 全区人妻精品视频| 国产成人a∨麻豆精品| 人妻 亚洲 视频| 日韩不卡一区二区三区视频在线| 99热这里只有是精品在线观看| 久久精品人人爽人人爽视色| 大片电影免费在线观看免费| 99久国产av精品国产电影| 国产成人免费无遮挡视频| 人人澡人人妻人| 一二三四在线观看免费中文在 | 在线天堂最新版资源| 各种免费的搞黄视频| a 毛片基地| 国产一区有黄有色的免费视频| 国产成人精品久久久久久| 少妇的丰满在线观看| 亚洲 欧美一区二区三区| 美女视频免费永久观看网站| 五月开心婷婷网| 精品一区二区三区四区五区乱码 | 中文字幕精品免费在线观看视频 | xxx大片免费视频| 在现免费观看毛片| 男女无遮挡免费网站观看| 国产永久视频网站| 国国产精品蜜臀av免费| 欧美xxxx性猛交bbbb| 亚洲国产日韩一区二区| 国产亚洲午夜精品一区二区久久| 国产高清三级在线| 日本wwww免费看| 人妻系列 视频| 在线天堂最新版资源| xxxhd国产人妻xxx| 少妇的逼水好多| 久久久久久伊人网av| 在线观看www视频免费| 国产免费福利视频在线观看| videossex国产| 亚洲在久久综合| 国产亚洲午夜精品一区二区久久| 久久国产精品男人的天堂亚洲 | 亚洲国产看品久久| 久久久久视频综合| 91成人精品电影| av福利片在线| 亚洲av综合色区一区| 国产成人精品无人区| 好男人视频免费观看在线| 国产毛片在线视频| 制服人妻中文乱码| 在线观看免费视频网站a站| 日韩制服骚丝袜av| 九九爱精品视频在线观看| 最新的欧美精品一区二区| 春色校园在线视频观看| 老司机影院毛片| 人人澡人人妻人| 久热这里只有精品99| 日日撸夜夜添| 亚洲美女搞黄在线观看| 国产欧美日韩综合在线一区二区| 国产精品久久久久久久久免| 只有这里有精品99| 老司机影院成人| 亚洲人成网站在线观看播放| 免费看av在线观看网站| 久久久久久人妻| 久久精品人人爽人人爽视色| 大码成人一级视频| 少妇人妻精品综合一区二区| 十八禁网站网址无遮挡| 成人国产av品久久久| videos熟女内射| 永久免费av网站大全| 丰满迷人的少妇在线观看| 久久精品国产a三级三级三级| 一级毛片黄色毛片免费观看视频| 七月丁香在线播放| 黄色毛片三级朝国网站| 欧美成人午夜免费资源| 免费高清在线观看日韩| 亚洲精品,欧美精品| 免费观看无遮挡的男女| 国产精品久久久久久精品古装| 久久女婷五月综合色啪小说| 日韩av免费高清视频| av一本久久久久| 少妇熟女欧美另类| 久久久久久久久久成人| 欧美日韩精品成人综合77777| 国产1区2区3区精品| 99久久人妻综合| 插逼视频在线观看| 亚洲国产精品专区欧美| 一本大道久久a久久精品| 80岁老熟妇乱子伦牲交| 欧美日韩精品成人综合77777| 国产精品久久久久成人av| a级片在线免费高清观看视频| 哪个播放器可以免费观看大片| 精品卡一卡二卡四卡免费| videos熟女内射| 国产免费福利视频在线观看| 亚洲欧洲精品一区二区精品久久久 | 亚洲精品第二区| 精品酒店卫生间| 欧美精品国产亚洲| 99热这里只有是精品在线观看| 成人亚洲精品一区在线观看| 精品午夜福利在线看| 欧美丝袜亚洲另类| 在线天堂最新版资源| 中文字幕最新亚洲高清| 国产深夜福利视频在线观看| 交换朋友夫妻互换小说| 日本爱情动作片www.在线观看| 欧美精品高潮呻吟av久久| 如日韩欧美国产精品一区二区三区| 狂野欧美激情性bbbbbb| 日韩一本色道免费dvd| 曰老女人黄片| av网站免费在线观看视频| 五月天丁香电影| 亚洲第一区二区三区不卡| 一本久久精品| 男女免费视频国产| 欧美日韩av久久| 边亲边吃奶的免费视频| 亚洲欧美一区二区三区国产| 视频在线观看一区二区三区| 熟女人妻精品中文字幕| 亚洲精品乱码久久久久久按摩| 97人妻天天添夜夜摸| 人成视频在线观看免费观看| 最近2019中文字幕mv第一页| 蜜臀久久99精品久久宅男| 99热全是精品| 视频在线观看一区二区三区| 久久99热6这里只有精品| 亚洲一区二区三区欧美精品| 国产精品一二三区在线看| 国产极品天堂在线| 麻豆精品久久久久久蜜桃| 婷婷色麻豆天堂久久| 51国产日韩欧美| 午夜日本视频在线| 99热全是精品| 波多野结衣一区麻豆| 亚洲精品456在线播放app| 国产 一区精品| 黄色 视频免费看| 国产精品无大码| 亚洲成人一二三区av| 黄片无遮挡物在线观看| 女的被弄到高潮叫床怎么办| 日韩精品有码人妻一区| 国产日韩欧美亚洲二区| 少妇的丰满在线观看| 在线观看一区二区三区激情| 少妇的逼好多水| 久热这里只有精品99| 丰满乱子伦码专区| 伦精品一区二区三区| 校园人妻丝袜中文字幕| videosex国产| 午夜91福利影院| 一级毛片黄色毛片免费观看视频| 成年人午夜在线观看视频| 满18在线观看网站| 亚洲精品第二区| av在线观看视频网站免费| 性色av一级| 两性夫妻黄色片 | 久久久亚洲精品成人影院| 少妇 在线观看| 久久99精品国语久久久| 女的被弄到高潮叫床怎么办| 精品一区二区免费观看| 久久国产精品男人的天堂亚洲 | 国产在线视频一区二区| 99久久综合免费| 大陆偷拍与自拍| 国语对白做爰xxxⅹ性视频网站| av.在线天堂| 成年人午夜在线观看视频| 精品人妻偷拍中文字幕| 性色av一级| 婷婷色综合大香蕉| 国产乱人偷精品视频| 亚洲国产精品专区欧美| 男人添女人高潮全过程视频| 中文字幕最新亚洲高清| 看免费成人av毛片| 人妻人人澡人人爽人人| 欧美精品人与动牲交sv欧美| 久久久精品区二区三区| 国产极品粉嫩免费观看在线| 国产深夜福利视频在线观看| 五月玫瑰六月丁香| 精品视频人人做人人爽| 狂野欧美激情性xxxx在线观看| a级片在线免费高清观看视频| 综合色丁香网| 国产成人一区二区在线| 黄色毛片三级朝国网站| 中文字幕人妻丝袜制服| 丝袜脚勾引网站| 国产极品天堂在线| 精品99又大又爽又粗少妇毛片| 国产成人91sexporn| 亚洲第一av免费看| 久久ye,这里只有精品| 久久精品国产亚洲av涩爱| 亚洲欧美日韩卡通动漫| 欧美日韩一区二区视频在线观看视频在线| 国产 精品1| freevideosex欧美| 狠狠精品人妻久久久久久综合| 菩萨蛮人人尽说江南好唐韦庄| tube8黄色片| 亚洲欧洲日产国产| 国产成人午夜福利电影在线观看| 久久97久久精品| 最近中文字幕2019免费版| 在线 av 中文字幕| av播播在线观看一区| 成年av动漫网址| 欧美精品av麻豆av| 日韩成人伦理影院| 热re99久久国产66热| 久久人妻熟女aⅴ| 亚洲人与动物交配视频| 最黄视频免费看| 极品少妇高潮喷水抽搐| 国产亚洲午夜精品一区二区久久| 十八禁网站网址无遮挡| 国产在视频线精品| 日韩一本色道免费dvd| 国产精品久久久久久精品古装| 一区二区av电影网| 免费人成在线观看视频色| 国产高清国产精品国产三级| 免费大片黄手机在线观看| 一区二区三区乱码不卡18| 在线看a的网站| 在线免费观看不下载黄p国产| 亚洲五月色婷婷综合| 只有这里有精品99| 国产在线一区二区三区精| 精品午夜福利在线看| kizo精华| 91精品国产国语对白视频| 国产av码专区亚洲av| 亚洲三级黄色毛片| 9色porny在线观看| 色网站视频免费| 久久久久网色| 欧美精品一区二区免费开放| 极品少妇高潮喷水抽搐| 久久久精品94久久精品| 亚洲中文av在线| 99久久精品国产国产毛片| 亚洲av国产av综合av卡| 高清毛片免费看| 啦啦啦在线观看免费高清www| 尾随美女入室| 日韩欧美精品免费久久| 婷婷色麻豆天堂久久| 天天躁夜夜躁狠狠躁躁| 少妇的逼好多水| 色5月婷婷丁香| 91精品伊人久久大香线蕉| 国产福利在线免费观看视频| 婷婷色av中文字幕| 七月丁香在线播放| 亚洲欧美成人精品一区二区| 18禁在线无遮挡免费观看视频| 久久精品人人爽人人爽视色| 91精品伊人久久大香线蕉| 成人免费观看视频高清| 卡戴珊不雅视频在线播放| 宅男免费午夜| 狂野欧美激情性xxxx在线观看| 秋霞伦理黄片| 亚洲内射少妇av| 亚洲第一av免费看| 亚洲欧洲国产日韩| 亚洲精品国产av蜜桃| 国产极品粉嫩免费观看在线| 久久久久久久大尺度免费视频| 美女主播在线视频| 女的被弄到高潮叫床怎么办| 中文字幕人妻熟女乱码| 18禁观看日本| 久久精品国产自在天天线| 国产爽快片一区二区三区| 国产极品粉嫩免费观看在线| 欧美成人午夜免费资源| 少妇人妻久久综合中文| av黄色大香蕉| 亚洲成国产人片在线观看| 亚洲第一av免费看| 又粗又硬又长又爽又黄的视频| 精品一区二区三卡| 色哟哟·www| 国产又色又爽无遮挡免| 9色porny在线观看| 亚洲久久久国产精品| 国产免费一区二区三区四区乱码| 午夜日本视频在线| 精品人妻一区二区三区麻豆| 国产日韩一区二区三区精品不卡| 熟女电影av网| 亚洲成人一二三区av| 国产一区亚洲一区在线观看| 精品人妻熟女毛片av久久网站| 99久国产av精品国产电影| 亚洲熟女精品中文字幕| 精品国产一区二区三区四区第35| 最近中文字幕高清免费大全6| 男女免费视频国产| 久久久久久人人人人人| 一级a做视频免费观看| 看免费成人av毛片| 午夜激情av网站| 精品国产乱码久久久久久小说| 亚洲精品久久成人aⅴ小说| 91精品国产国语对白视频| 久久久亚洲精品成人影院| 久久精品国产亚洲av天美| 免费av中文字幕在线| 亚洲av电影在线进入| 一边摸一边做爽爽视频免费| 十分钟在线观看高清视频www| 免费在线观看完整版高清| 中文字幕精品免费在线观看视频 | 亚洲图色成人| 亚洲国产精品一区三区| 肉色欧美久久久久久久蜜桃| 久热久热在线精品观看| 国产精品三级大全| 成人亚洲欧美一区二区av| 久久99热6这里只有精品| 爱豆传媒免费全集在线观看| 国产在线视频一区二区| 人妻少妇偷人精品九色| 国产精品秋霞免费鲁丝片| 亚洲欧美成人精品一区二区| 又黄又粗又硬又大视频| 婷婷色综合大香蕉| 少妇被粗大的猛进出69影院 | 免费大片黄手机在线观看| 国产精品一区www在线观看| 国产一区二区在线观看av| 亚洲,欧美精品.| 亚洲av免费高清在线观看| 精品国产乱码久久久久久小说| 成人午夜精彩视频在线观看| 精品久久久精品久久久| 波多野结衣一区麻豆| 秋霞伦理黄片| 日韩欧美精品免费久久| 男人操女人黄网站| 两个人免费观看高清视频| 最后的刺客免费高清国语| 欧美日韩视频精品一区| 久热这里只有精品99| 一区二区三区乱码不卡18| 波多野结衣一区麻豆| 久热这里只有精品99| 丰满少妇做爰视频| 日本色播在线视频| av视频免费观看在线观看| 日韩免费高清中文字幕av| 美女主播在线视频| 国产又色又爽无遮挡免| av不卡在线播放| 亚洲欧美成人综合另类久久久| 亚洲,欧美精品.| 精品一区二区三区视频在线| 菩萨蛮人人尽说江南好唐韦庄| 美女主播在线视频| 一区二区日韩欧美中文字幕 | 在线观看一区二区三区激情| 男女国产视频网站| 在线天堂中文资源库| 国产1区2区3区精品| 亚洲中文av在线| 亚洲欧美清纯卡通| 成人毛片60女人毛片免费| 建设人人有责人人尽责人人享有的| 国产精品偷伦视频观看了| 中文字幕亚洲精品专区| 国产精品人妻久久久影院| 人人妻人人爽人人添夜夜欢视频| 精品久久久精品久久久|