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

    乳山灣鄰近海域有機(jī)碳的分布與底界面過程

    2017-04-11 14:39:06臧家業(yè)趙晨英謝琳萍王以斌張愛軍冉祥濱國家海洋局第一海洋研究所海洋生態(tài)研究中心山東青島26606國家海洋局第一海洋研究所海洋地質(zhì)與地球物理研究室山東青島26606
    中國環(huán)境科學(xué) 2017年3期
    關(guān)鍵詞:乳山表層沉積物

    臧家業(yè),趙晨英,劉 軍,2,謝琳萍,王以斌,張愛軍,冉祥濱*(.國家海洋局第一海洋研究所,海洋生態(tài)研究中心,山東 青島 26606;2.國家海洋局第一海洋研究所,海洋地質(zhì)與地球物理研究室,山東 青島 26606)

    乳山灣鄰近海域有機(jī)碳的分布與底界面過程

    臧家業(yè)1,趙晨英1,劉 軍1,2,謝琳萍1,王以斌1,張愛軍1,冉祥濱1*(1.國家海洋局第一海洋研究所,海洋生態(tài)研究中心,山東 青島 266061;2.國家海洋局第一海洋研究所,海洋地質(zhì)與地球物理研究室,山東 青島 266061)

    根據(jù)2009年和2014年夏季在乳山灣口及鄰近海域的綜合調(diào)查結(jié)果,分析了該海域夏季有機(jī)碳的時(shí)空分布、底界面過程與影響因素.結(jié)果表明,2009年夏季乳山灣近海水體溶解有機(jī)碳(DOC)含量介于0.70~3.19mg/L之間,平均值為1.80mg/L;DOC 的平均值在8月最高,7月與 9月次之,6月最低;2014年 8月份 DOC的變化范圍為 1.79~15.2mg/L,高于 2009年同期水平,顆粒有機(jī)碳(POC)的變化范圍為0.04~1.33mg/L;水體有機(jī)碳的分布受陸源輸入、海洋初級生產(chǎn)以及潮汐的顯著影響.研究區(qū)域夏季顆粒有機(jī)碳(POC)的沉降通量為(25±0.8)g/m2,約占初級生產(chǎn)固碳量的66%;沉積物上層(0~4cm)間隙水中DOC的濃度是沉積物上覆水的8~9倍,DOC在沉積物—水界面存在向上覆水釋放的現(xiàn)象;乳山灣灣口DOC交換通量為14.4—97g/(m2·a),占水體存量的1.1%~13.4%.人類活動一定程度上影響了乳山灣及其近海有機(jī)碳的構(gòu)成與循環(huán)收支過程,是區(qū)域環(huán)境變化的重要驅(qū)動因子之一.有超過50%的有機(jī)碳會隨潮流輸送到外海,顯示潮流在有機(jī)碳輸送中巨大作用;沉積物-水界面 DOC的交換會影響底界面有機(jī)碳的收支與循環(huán)過程,有機(jī)碳的收支表明研究海域底界面有機(jī)碳的降解所產(chǎn)生的溶解有機(jī)碳是水體 DOC的重要來源,最終保存在沉積物中的碳其埋藏量約占初級生產(chǎn)的 13%;相對較低的溶解氧水平可能會增加DOC的交換通量,影響碳在陸架邊緣海的埋藏.

    乳山灣;有機(jī)碳;通量;初級生產(chǎn);底界面過程

    海洋是地球上最大的碳儲庫,碳儲量為 38Tt (T=1012),為大氣碳儲庫的50倍,是陸地碳儲庫的20倍[1],因此在調(diào)節(jié)大氣碳含量和全球氣候變化中的作用關(guān)鍵.碳在海洋中的存在形式主要是溶解無機(jī)碳、溶解有機(jī)碳、顆粒有機(jī)碳和碳酸鹽等

    [1].海水中溶解有機(jī)碳(DOC)作為水生態(tài)系統(tǒng)中最大的一類有機(jī)碳庫,約占海洋有機(jī)碳的80%~95%[2],是海水中主要的有機(jī)碳形式;同時(shí), DOC也是表征水體中有機(jī)物含量、生物活動水平、初級生產(chǎn)力和有機(jī)質(zhì)污染[3]的重要參數(shù)之一,在海洋碳循環(huán)和海洋生態(tài)系統(tǒng)功能變化以及環(huán)境評價(jià)的研究中具有重要的價(jià)值.

    海洋碳循環(huán)分別受到外部和內(nèi)部過程的綜合影響;外部過程多指陸-海、海-氣交換等界面過程;內(nèi)部過程則指碳在海洋環(huán)流和海洋生態(tài)系統(tǒng)作用的驅(qū)動下進(jìn)行的遷移與轉(zhuǎn)化等過程[4-6].實(shí)際上,底界面過程在區(qū)域碳循環(huán)中的作用不容忽視,有機(jī)碳(主要是顆粒有機(jī)碳)在底界面經(jīng)歷反復(fù)的吸附-解吸、沉積-再懸浮和礦化降解過程,最終僅有一小部分惰性有機(jī)碳得以長期保存[4-5].陸架區(qū)的碳循環(huán)是全球碳收支的重要一環(huán),同時(shí)還是全球物質(zhì)源-匯角色的敏感區(qū).海灣及其鄰近海域作為陸海相互作用最集中的區(qū)域,在區(qū)域碳的生物地球化學(xué)循環(huán)研究中具有重要的地位[5-6].

    目前國內(nèi)關(guān)于海洋有機(jī)碳的研究主要集中在有機(jī)碳的分布及影響因素等方面,且大多集中在渤、黃海[6-9]和東海[10]等開闊的近海區(qū)域,對于典型海灣的研究較少,且缺少底界面過程的研究,導(dǎo)致無法準(zhǔn)確量化區(qū)域有機(jī)碳循環(huán)的過程和人類活動對碳循環(huán)的影響.乳山灣及鄰近海域是中國北方重要的經(jīng)濟(jì)貝類和魚類養(yǎng)殖區(qū),其海洋環(huán)境狀況深受日益加強(qiáng)的人類活動的影響,特別是近年來在該海域發(fā)現(xiàn)的富營養(yǎng)化[11]、綠潮[12]、重金屬污染[13]和季節(jié)性低氧等現(xiàn)象[14-15]可能影響該區(qū)域碳循環(huán)及環(huán)境變化.本研究通過2009年夏季和2014年8月對乳山灣外鄰近海域的綜合調(diào)查,分析了乳山灣外鄰近海域水體和沉積物中有機(jī)碳的分布特征,并結(jié)合水文、生物等同步觀測數(shù)據(jù),探討了有機(jī)碳與水文、化學(xué)及生物參數(shù)之間的關(guān)系,旨在識別典型海灣(近岸養(yǎng)殖區(qū))影響有機(jī)碳分布和循環(huán)過程的主要因素,并為揭示人類活動影響下養(yǎng)殖區(qū)有機(jī)碳的生物地球化學(xué)過程提供依據(jù).

    1 材料與方法

    1.1 采樣站位

    圖1 乳山灣口及其鄰近海域采樣站位Fig.1 Sampling stations in the coastal area of Rushan Bay (虛線為收支計(jì)算的邊界)

    于2009年6~9月(分別為6月14~17、7月15~19、8月15~23與9月15~17)在山東威海乳山灣口及其鄰近海域進(jìn)行了逐月共 4個(gè)航次的綜合調(diào)查,觀測期間調(diào)查水域的平均水深約為20m,用于分析有機(jī)碳的分布與影響因素;采樣站位如圖1所示,其中A、B、C、D、E和F站位位于灣口,其它位于灣口外鄰近水域.另于 2014年5月、7月和8月大潮期在灣口鄰近海域選取了 6個(gè)代表性站位采集了懸浮顆粒物和表層沉積物(圖1),沉積物采樣站位的分布力求覆蓋整個(gè)調(diào)查區(qū)域,并考慮了陸源輸入和潮汐等因素的影響;同時(shí)于2014年7月和8月在C1、C2和C5站位現(xiàn)場進(jìn)行了初級生產(chǎn)力的測定;于2014年8月在C1和C2站采集了沉積物柱狀樣,并在C2站位進(jìn)行了26h的連續(xù)觀測;2014年獲取的數(shù)據(jù)用于量化碳循環(huán)的主要過程.

    1.2 樣品采集與分析

    1.2.1 水樣采集 調(diào)查中使用溫鹽深采水器(JFE Co. Ltd,AAQ122,Japan)采集水樣,同時(shí)獲取水體溫度、鹽度以及透明度等水文參數(shù);采樣層次為表層(0.5m)、5m、10m和近底層,并根據(jù)水深增加或減少層次.水樣采集后,先用虹吸管取一定已知體積的水樣于碘量瓶中,加入固定劑后利用溫科勒法進(jìn)行DO的分析.再取一定體積水樣用孔徑 0.70μm玻璃纖維濾膜(GF/F) (預(yù)先用1:1000HCl溶液浸泡24h,并以Milli-Q水洗至中性,450℃馬弗爐中灼燒備用)過濾,濾液用于溶解有機(jī)碳(DOC)的分析,濾膜用于顆粒有機(jī)碳(POC)和顆??偟?PN)的分析.另取一定體積水樣用孔徑 0.45μm醋酸纖維膜收集(預(yù)先用 1:1000HCl浸泡24h,并以Milli-Q水洗至中性,烘箱內(nèi)45℃烘干),濾膜避光-20℃冷凍保存用于葉綠素 a (Chl-a)的分析.最后在現(xiàn)場分裝一定體積的未過濾水樣,一份避光冷凍保存用于總氮(TN)和總磷(TP)的測定,另一份避光保存用于COD的分析.

    1.2.2 沉積物采集 沉積物柱狀樣使用重力取樣器(KC-Denmark Ltd., Denmark)獲取.柱狀樣采集后,立刻用一次性注射器采集上覆水,現(xiàn)場過濾后,-20℃冷凍保存,用于測定 DOC含量;上覆水采集后,柱狀樣立刻放入充氮?dú)獾氖痔紫鋬?nèi)分割取樣,取樣間隔為1~2cm;其中深度為0~8cm的樣品,取一部分分割后的沉積物離心過濾制取間隙水(4000r/min,10min),-20℃冷凍保存用于DOC的分析;另取一部分沉積物于密封袋中冷凍保存,用于總有機(jī)碳(TOC)、總氮(TPN)和穩(wěn)定碳同位素的分析.

    1.2.3 樣品分析方法 DO和 COD等均按照《海洋監(jiān)測規(guī)范》(GB 17378-2007)中方法進(jìn)行分析處理.TN和TP采用過硫酸鉀(K2S2O8)氧化法(GB17378-2007)測定,回收率為90%,相對標(biāo)準(zhǔn)偏差<10% (n=5).Chl-a樣品用 9:1丙酮提取后,使用熒光分光光度計(jì)(儀電科儀,970CRT,上海)測定(GB 17378-2007).DOC采用總有機(jī)碳分析儀(日本島津公司,TOC-CCPH型)測定,分析誤差為±0.005mg/L,相對標(biāo)準(zhǔn)偏差<2% (n=6).POC和PN采用元素分析儀測定(德國EURO公司,EA 3000型);測定前,將濾膜置于烘箱內(nèi)(45℃)烘72h以上,恒重后稱重;之后用5‰的HCl溶液洗掉無機(jī)態(tài)的碳再用 Milli-Q水清洗掉濾膜上殘留的鹽酸,烘干后稱重(45℃,72h),上機(jī)測試;POC和PN的分析誤差分別為±0.01mg/L和0.003mg/L,相對標(biāo)準(zhǔn)偏差均小于10% (n=6)[9].

    沉積物樣品冷凍干燥后研磨過200目篩,充分混合后取0.3g左右加入5mL HCl溶液(1+1)浸泡6h,再加入新的HCl溶液直至不再有氣泡冒出,反應(yīng)完全后離心去掉酸液(2000r/min, 5min), Milli-Q水清洗數(shù)遍至中性,低溫烘干至恒重后稱重;稱取 10mg左右處理后的樣品使用元素分析儀—連續(xù)流同位素比值質(zhì)譜聯(lián)用儀(Vario Micro Cube-IsoPrime 100,Elementar,德國)進(jìn)行沉積物中TOC、TPN和穩(wěn)定碳同位素分析[9,16].

    沉積物利用酸洗前后重量差校正得到沉積物樣品TOC和TPN的百分含量.TOC和TPN的分析誤差分別為±0.02%和±0.005%,相對標(biāo)準(zhǔn)偏差均<10% (n=6)[15].穩(wěn)定碳同位素值校正到Vienna Peedee Belemnite (V-PDB)標(biāo)準(zhǔn),有機(jī)碳穩(wěn)定同位素豐度(δ13C)的計(jì)算方法為:

    式中:R樣品和 R標(biāo)準(zhǔn)分別是樣品和標(biāo)準(zhǔn)的13C/12C值.δ13C分析精度優(yōu)于±0.1‰ (n=10).

    1.3 初級生產(chǎn)力與相應(yīng)的DOC的生產(chǎn)量

    初級生產(chǎn)力固碳量的測定采用葉綠素法,具體計(jì)算如下[17]:

    式中:FP為真光層浮游植物初級生產(chǎn)固碳量,g C/(m2·a);Ps為表層水體中浮游植物的潛在生產(chǎn)力,mg C/(m3·h);Z為真光層深度,m;Z=7.5m,為透明度的3倍;t為每日光照時(shí)間,h,取8h.式(2)中浮游植物的潛在生產(chǎn)力即浮游植物光合作用速率的測定采用黑白瓶法,現(xiàn)場使用CTD采水器采集表層水樣,一部分水樣用于測定DO和Chl-a的濃度,另一部分水樣分裝于4個(gè)1L的透明磨口瓶中,其中 2瓶用鋁箔紙包裹(黑瓶);將黑白瓶放回到原采樣深度培養(yǎng),8h后取出,立即測定其DO含量;根據(jù)DO的變化、葉綠素a的濃度以及培養(yǎng)時(shí)間等,求出調(diào)查區(qū)域的平均同化系數(shù)(Q),之后根據(jù)葉綠素濃度和同化系數(shù)求出浮游植物的潛在生產(chǎn)力.

    式中:Chl-a為葉綠素a的濃度,mg/m3;Q為同化系數(shù),mg C/[(mg Chl-a)·h];h為光照時(shí)間,h,取8h;Ol和Od分別為白瓶和黑瓶的DO濃度,mg/L.

    1.4 界面通量計(jì)算

    1.4.1 有機(jī)碳沉積通量 乳山灣有機(jī)碳沉積通量根據(jù)公式(5)和(6)進(jìn)行計(jì)算[15]:

    式中:FS為有機(jī)碳沉積通量,g/(m2·a);C0為表層沉積物 TOC含量,%;MAR為沉積物累積速率,g/(m2·a).

    1.4.2 底界面溶解有機(jī)碳交換通量 DOC在沉積物—水界面的交換主要考慮分子擴(kuò)散、生物擾動和生物灌溉等[18-20],分子擴(kuò)散、生物擾動引起的交換通量根據(jù)以下公式計(jì)算[18-19,21]:

    式中:FE為DOC交換通量,g/(m2·a),負(fù)號表示擴(kuò)散是由沉積物向水體方向進(jìn)行的;φ為沉積物孔隙率;DD為交換系數(shù);DM為沉積物中DOC的分子擴(kuò)散系數(shù); DB為生物擾動系數(shù);?C/?z為沉積物—水界面處DOC的濃度梯度;D0為無限稀釋溶液中溶質(zhì)的分子擴(kuò)散系數(shù)[19];m 為經(jīng)驗(yàn)系數(shù)(φ≤0.7, m=2;φ>0.7,m=2.5~3.0)[22];r為沉積速率.FG為生物灌溉作用引起的 DOC交換通量,本文利用Hong等[20]在九龍江口基于224Ra/228Th不平衡得到的224Ra的分子擴(kuò)散通量(FM)和生物灌溉通量(FM)的比值(1:2.5~1:21,夏季河口和入??谡疚?來估算生物灌溉作用對DOC在沉積物—水界面交換通量的貢獻(xiàn).

    1.4.3 區(qū)域邊界有機(jī)碳交換通量 以 121.4°~121.7°E和36.4°~36.7°N為邊界(圖1),進(jìn)行系統(tǒng)內(nèi)部有機(jī)碳的收支計(jì)算(圖2).利用美國普林斯頓大學(xué)海洋模式(POM)建立了的變網(wǎng)格的潮汐潮流數(shù)值模式[23],模擬研究海區(qū)的分潮,并提取主要分潮的調(diào)和常數(shù).以潮汐潮流模式為基礎(chǔ),在POM的計(jì)算框架下,利用T_tide預(yù)報(bào)程序提取研究區(qū)域流場(包括流速和流向的數(shù)據(jù)),結(jié)合水深和收支系統(tǒng)內(nèi)外水體有機(jī)碳數(shù)據(jù),計(jì)算水平輸送作用對區(qū)域有機(jī)碳的輸入和輸出的貢獻(xiàn).

    圖2 平流作用下各界面水體有機(jī)碳收支計(jì)算框架Fig.2 Schematic of organic carbon budget across the boundaries of the study areaFi為各界面有機(jī)碳的凈通量,由水量和對應(yīng)的濃度計(jì)算得到;i=1~4

    1.4.4 內(nèi)源貢獻(xiàn)率 通過沉積物—水界面DOC的交換通量,評估沉積物作為內(nèi)源對水體DOC的貢獻(xiàn)[15],內(nèi)源釋放的貢獻(xiàn)比例按如下公式計(jì)算:

    式中:LDOC為內(nèi)源DOC的貢獻(xiàn)率,%;FE為DOC在底界面的交換通量;CDOC為水體 2009年夏季DOC的平均濃度,為 1.80mg/L;twr為水體存留時(shí)間(11.6d,根據(jù) ADCP所測的水體平均流速計(jì)算得到[15]);D為平均水深,為19.5m.

    1.5 物源分析

    陸源輸入對顆粒有機(jī)碳貢獻(xiàn)(WT,%)通過端元模型[24]計(jì)算得到:

    式中:δ13CM(-19.5‰)和 δ13CT(-27‰)分別為海源和陸源有機(jī)碳的δ13C值[24],在中國東部陸架海區(qū)有機(jī)碳物源的研究中廣泛應(yīng)用[6,16].

    2 結(jié)果

    2.1 溫度、鹽度、COD、DO、總氮、總磷和葉綠素a

    如表1所示,6月份溫度較低,7月份、8月份和 9月份的溫度并無明顯差異,4個(gè)月份溫度均表現(xiàn)為表層顯著高于底層,并且出現(xiàn)了層化現(xiàn)象.整個(gè)調(diào)查區(qū)域夏季的鹽度穩(wěn)定,沒有明顯的表底層差異.TN和TP均具有顯著的時(shí)空分布差異,6月和9月TN和TP濃度相對較高,7-8月相對較低,并在 9月份底層均顯著高于表層水體(P<0.05).COD和葉綠素a的濃度在7月份、8月份和9月份較高,6月份最低.4個(gè)月份DO濃度均表現(xiàn)為表層高于底層,底層水體DO低值現(xiàn)象在8月份最為突出.

    表1 2009年夏季乳山灣外鄰近海域溫度、鹽度、COD、DO、總氮、總磷和葉綠素a的變化范圍Table 1 Temperature, salinity, chemical oxygen demand, total and chlorophyll in the water column in the adjacent area of Rushan Bay in summer 2009

    2.2 水體溶解有機(jī)碳

    2.2.1 季節(jié)變化 乳山灣鄰近海域DOC含量6月份最低,8月份含量最高(表3),顯著高于其它月份(P<0.01,n=71),7月份和9月份DOC含量沒有顯著差異.各月中,6月份表層DOC含量高于底層(P<0.01,n=23),7月份和8月份表、底層DOC含量相差不大,而9月份底層DOC含量較表層高(P<0.01,n=22).

    表2 乳山灣外鄰近海域2009年夏季DOC的濃度與范圍(mg/L)Table 2 Results of DOC in the adjacent area of Rushan Bay in summer 2009 (mg/L)

    2.2.2 表層 DOC分布 調(diào)查海域外側(cè) DOC的含量較低,中部海域DOC含量較高,表層水體在C1和D6站位出現(xiàn)2個(gè)高值中心.6月份與7月份相似,表層高值區(qū)主要出現(xiàn)在乳山灣口及灣口東部的近岸海域,整體呈由近岸向離岸DOC濃度遞減的變化趨勢.8月份表層中部海域DOC濃度較高,DOC高值水舌由西南部向中部海域擴(kuò)展.9月份水體表層調(diào)查海域DOC的分布呈垂直與岸線的帶狀分布,中間 C斷面DOC含量較高.可見乳山灣鄰近海域DOC的分布的共同特征是灣口 DOC濃度高于離岸水域(圖3).

    圖3 乳山灣外鄰近海域2009年夏季表層 DOC分布(mg/L)Fig.3 Horizontal distribution of DOC in the adjacent area of Rushan Bay in summer 2009 (mg/L)

    2.2.3 定點(diǎn)觀測 2014年8月份C2站26h定點(diǎn)觀測中,表層和底層 DOC的含量分別介于2.02~ 15.2mg/L和1.79~10.4mg/L之間,平均值為5.63, 4.48mg/L.定點(diǎn)觀測中表層DOC顯著高于底層(P<0.05);表層DOC含量隨潮汐表現(xiàn)出規(guī)律性的變化,在高潮位濃度較高,而在低潮位含量普遍較低(圖4a).

    圖4 乳山灣外2014年8月份定點(diǎn)觀測(C2)中有機(jī)碳隨時(shí)間的變化Fig.4 Time-dependent variation of DOC and POC at C2station in the coastal area of Rushan Bay in August 2014

    2.3 顆粒有機(jī)碳氮

    2014年夏季調(diào)查海域 POC的變化范圍為0.04~1.33mg/L,平均濃度為(0.50±0.33)mg/L.PN的變化范圍為 0.01~0.27mg/L,平均濃度為(0.12± 0.08)mg/L.乳山灣外POC和PN表、底層分布特征一致,與 DOC一致呈現(xiàn)近岸相對較高,離岸相對較低的分布趨勢;高值區(qū)位于灣口附近水域,低值區(qū)位于灣外深水區(qū)域.

    26h定點(diǎn)觀測中,表層和底層POC的含量分別介于0.52~1.33mg/L和0.04~0.92mg/L之間,平均值為 0.82mg/L,0.57mg/L,顯著高于底層(P<0.05);表層POC也隨潮汐表現(xiàn)出規(guī)律性的變化,在高潮位濃度較高,而在低潮位含量普遍較低(圖4b).

    2.4 沉積物中有機(jī)碳氮

    2.4.1 間隙水中溶解有機(jī)碳 間隙水中 DOC的剖面分布如圖 5所示.C1和 C2站位間隙水DOC的變化范圍分別為6.55~24.2mg/L和3.05~25.6mg/L,平均值分別為 14.0mg/L和 18.1mg/L.上層(0~4cm)間隙水中 DOC的含量是沉積物上覆水的8~9倍,上覆水—沉積物界面DOC明顯的濃度梯度表明調(diào)查站位存在DOC從沉積物向上覆海水?dāng)U散的現(xiàn)象.

    2.4.2 沉積物中總有機(jī)碳和總氮 研究區(qū)域表層沉積物中TOC變化范圍為0.40%~0.88%,平均值為(0.60±0.16)%;TPN的變化范圍為 0.10%~0.17%,平均值為(0.14±0.12)%.柱狀沉積物中TOC的含量在沉積物上層(0~15cm)的變化較大,15cm以深相對穩(wěn)定(圖5).

    圖5 沉積物柱狀樣中TOC和間隙水中DOC的剖面變化Fig.5 Vertical distribution of DOC in pore water and TOC content in cores

    2.5 初級生產(chǎn)

    浮游植物光合作用速率的測定中,白瓶中DO7.29~8.90mg/L,黑瓶DO介于6.59~7.71mg/L, Chl-a的濃度介于 1.0~2.79μg/L,同化系數(shù)為0.86~85.2mgC/[(mg Chl-a)·h].

    2.6 顆粒和沉積有機(jī)碳物源分析

    圖6 乳山灣外懸浮顆粒物(SPM)和沉積物中C/N和δ13C的分配及來源Fig.6 Distinctive source combinations of atomic C/N ratios and organic δ13C values of SPM and sediment in the coastal area of Rushan Bay

    2014年夏季水體懸浮顆粒物和沉積物中的有機(jī)碳和氮的比值介于2.15~19.2之間,平均值為5.56(C/N,物質(zhì)的量比).δ13C 介于-24.63‰~-21.79‰之間,平均值為-23.11‰;乳山灣外懸浮顆粒物(SPM)和沉積物中C/N和δ13C的分配及來源(圖6)顯示該區(qū)域有機(jī)碳分別由陸源和海源有機(jī)碳組成;通過端元模型計(jì)算陸源輸入的有機(jī)碳比例介于31%~68%,平均值為47%±10%;海源輸入的有機(jī)碳所占比例介于 32%~69%之間,平均值為53%±10%.

    2.7 沉積物底界面過程

    根據(jù)公式(6)~公式(9)及現(xiàn)場觀測(表3),乳山灣C1和C2站沉積物—水界面DOC的交換通量分別為14.4~76.5g/(m2·a)和18.3~97.0g/(m2·a),其貢獻(xiàn)量為水體 DOC儲量(依據(jù)公式(10)計(jì)算)的 1.1%~13.4%.C1和C2站TOC沉積通量分別為102,97g/ (m2·a),約占初級生產(chǎn)固碳量的68%和64%.

    表3 沉積物底界面通量計(jì)算相關(guān)參數(shù)Table 3 Parameters for calculation of organic carbon flux at the sediment-water interface

    3 討論

    3.1 海洋水環(huán)境變化對溶解有機(jī)碳分布和水平的影響

    2009年夏季調(diào)查海域 DOC濃度在 0.70~3.19mg/L之間,平均值為1.80mg/L.這與2007年ST03區(qū)塊水體環(huán)境調(diào)查與研究(908-01-ST03)中重疊區(qū)域的平均值(1.87mg/L)無顯著性差異.不過,乳山灣鄰近海域不同時(shí)期的 DOC含量均顯著高于南黃海外部海域[8-9],同時(shí)灣口處的DOC含量(2.15mg/L)也明顯高于灣外的DOC含量(1.75mg/L; P<0.05).

    對于近岸而言,表層和次表層 DOC含量較高很大程度上是因?yàn)楦叩年懺摧斎牒宛B(yǎng)殖活動造成的.乳山灣近岸人口密集,人類活動產(chǎn)生的大量有機(jī)物和灣內(nèi)養(yǎng)殖活動排放入海的有機(jī)質(zhì)可能使近岸水體的 DOC含量較高.本次調(diào)查結(jié)果中,灣口的COD、Chl-a、TN和TP平均含量分別為(1.24mg/L、3.91μg/L、13.1μmol/L 和0.83μmol/L),而灣外其平均含量則分別為(1.13mg/L、1.67μg/L、12.0μmol/L和0.74μmol/L);均具有灣口高于灣外的特征,表明人類活動對DOC的分布具有一定的影響.在 2009年的調(diào)查中,乳山灣口 COD (0.71~1.77mg/L)明顯高于歷史同期(20世紀(jì)90年代的平均值為0.59mg/L[25]),這也從一個(gè)側(cè)面表明乳山灣養(yǎng)殖業(yè)及工農(nóng)業(yè)的快速發(fā)展導(dǎo)致排海有機(jī)物的量顯著增加,導(dǎo)致調(diào)查區(qū)域DOC含量進(jìn)一步升高.

    從表4可以看出,DOC與其他參數(shù)均存在一定的相關(guān)性,夏季乳山灣外鄰近海域水體中DOC含量與鹽度呈負(fù)相關(guān),表明 DOC的分布即受到陸源輸入的影響,還與外部海水的交換有關(guān).8月份乳山灣低值現(xiàn)象較為突出;水體DOC和DO呈顯著負(fù)相關(guān)關(guān)系,表明次氧化條件有利于顆粒態(tài)及表層沉積物中的有機(jī)質(zhì)降解,生成的 DOC釋放到水體中造成的8月份含量最高.

    調(diào)查區(qū)域水體DOC含量與COD、TN和TP濃度也有一定的正相關(guān)關(guān)系,尤其是8月份,水體中大量的營養(yǎng)元素滿足浮游植物的生長需求,較其它月份普遍較高的浮游植物光合作用(以Chl-a濃度表示)產(chǎn)生了相對較多的DOC,這也使得水體中的DOC含量在8月份明顯升高.由此可見,較高的初級生產(chǎn)也是影響 DOC含量和分布的重要因素.本研究中,DOC的高值往往與葉綠素 a的高值對應(yīng),這也間接表明浮游植物生長對DOC含量的貢獻(xiàn).

    表4 乳山灣外鄰近海域2009年夏季DOC與環(huán)境因子的相關(guān)關(guān)系Table 4 Correlations between DOC and other properties of Rushan Bay in summer 2009

    潮汐也會對表層DOC與POC的分布起到一定的作用,根據(jù)2014年8月份C2站位的分析發(fā)現(xiàn),高潮位時(shí)水體有機(jī)碳濃度較高,而低潮位時(shí)較低(圖 4).一方面,潮汐作用可能使得表、底層高含量的DOC和POC相互混合,從而造成高潮期水體有機(jī)碳的含量的升高.另一方面,潮汐作用促進(jìn)了近岸水體和離岸水體的混合,有利于將近岸高濃度的有機(jī)碳向離岸輸送,這種因潮汐而產(chǎn)生的稀釋作用間接造成了低潮期有機(jī)碳含量較低;連續(xù)站觀測中,表層水體DOC和POC濃度均與鹽度呈顯著負(fù)相關(guān)(P<0.05),表明潮汐作用的確影響了水體 DOC的時(shí)空分布.2個(gè)年份(2009與2014)DOC濃度的顯著差異很大程度上也是由于潮汐作用引起的,但 2009年的數(shù)據(jù)為大面調(diào)查所獲得的,難以體現(xiàn)潮汐因素對 DOC分布的影響.

    3.2 生物生長與溶解有機(jī)碳含量變化

    8月份,水體中DOC不僅與Chl-a呈顯著正相關(guān)(P<0.01),與 TN和 TP也呈顯著正相關(guān)(P<0.01),一方面反映了初級生產(chǎn)是 DOC的主要來源之一,浮游植物量的增加使得水體中 DOC含量升高;另一方面也說明了浮游植物通過初級生產(chǎn)將水體中無機(jī)態(tài)的氮、磷等營養(yǎng)元素轉(zhuǎn)化成有機(jī)物.懸浮顆粒物中有機(jī)質(zhì)的C/N表明有機(jī)碳主要來源于海洋浮游藻類[25],進(jìn)一步表明水體有機(jī)碳(POC和DOC)主要受初級生產(chǎn)的控制,這與碳穩(wěn)定同位素的分析結(jié)果相一致.定點(diǎn)觀測中表層DOC和 POC顯著相關(guān)的線性關(guān)系(R2=0.32, P<0.05)也表明DOC與POC大致受共同因素的影響.

    本研究中也發(fā)現(xiàn)在夏季表層水體中的葉綠素含量高于底層,平均值達(dá)到1.73μg/L,而底層為1.56μg/L,這和 DOC的表底層分布差異相一致;此外,8月份表層葉綠素相對其他月份較高,平均值達(dá)到2.13μg/L,同時(shí)DOC含量也高于其它月份,說明了海洋自生初級生產(chǎn)是水體DOC的重要貢獻(xiàn)者.根據(jù)2014年黑白瓶培養(yǎng)實(shí)驗(yàn)結(jié)果計(jì)算得到初級生產(chǎn)力介于 26.4~311gC/(m2·a)之間,平均值為 151gC/(m2·a).已有報(bào)道表明浮游植物光合作用產(chǎn)生的DOC為固碳總量的5%~50%[26].乳山灣外浮游植物初級生產(chǎn)力固碳量介于 26.4~311gC/(m2·a)之間,根據(jù)上述比例估算可得到浮游植物光合作用產(chǎn)生 DOC的量約為 7.5~75.3gC/(m2·a).初級生產(chǎn)力的變化會使其產(chǎn)品結(jié)構(gòu)發(fā)生變化,最終導(dǎo)致溶解有機(jī)碳和顆粒有機(jī)碳的比例發(fā)生變化,從而影響有機(jī)碳的循環(huán)和保存.

    3.3 與其他海域比較

    與其他海區(qū)相比(表5),本次調(diào)查乳山灣外部海域 DOC的濃度與青島近岸DOC濃度相近(1.70mg/L)[27-28],較渤海(2.19mg/L)低[8];明顯高于東海、南海和大洋等開闊的海區(qū).乳山灣養(yǎng)殖排放大量的有機(jī)物使得DOC的濃度整體偏高于上述海域.另外,相對東海、南海和大洋高的富營養(yǎng)化水平也產(chǎn)生了較高的初級生產(chǎn)力,這也使得乳山灣外浮游植物所固定下來的有機(jī)物碳含量較高,DOC含量也隨之升高.

    表5 中國近海典型區(qū)域以及大洋水體中DOC質(zhì)量濃度Table 5 Concentrations of DOC in the oceans and typical offshore areas in China

    此外,夏季乳山灣鄰近海域 DO濃度(3.21~7.16mg/L)[29]顯著低于南黃海外海海域[30-31],除人類活動以及養(yǎng)殖業(yè)對近海DOC的分布的影響外,季節(jié)性低氧所產(chǎn)生的氧化還原環(huán)境的更替還有利于顆粒態(tài)和表層沉積物有機(jī)碳的降解,這一定程度上提高了乳山灣外水體DOC的濃度.

    3.4 有機(jī)碳底界面過程與收支

    3.4.1 底界面過程 海洋沉積物中的碳一部分是以有機(jī)碳形式存在,使得沉積物成為一個(gè)海洋碳循環(huán)中重要的匯.同時(shí),水-沉積物界面還存在著釋放、溶解和再懸浮等一系列物理化學(xué)作用,使得沉積物間隙水中的 DOC含量較高,進(jìn)而又通過擴(kuò)散作用進(jìn)入水體,成為水體的碳源.DOC在底界面的交換通量和 TOC的沉積通量表明,乳山灣沉積物總體具有碳匯這一性質(zhì),但同時(shí)又是水體中DOC的源.調(diào)查區(qū)域沉積物-水界面向上覆水?dāng)U散的DOC的量占水體 DOC存量的1.1%~13.4%,同時(shí)占上層水體浮游植物初級生產(chǎn)固碳量的9.5%~64%,說明沉積物對水體DOC具有較大的影響.

    從表6中可以看出,研究區(qū)域內(nèi)DOC的交換通量低于亞馬遜河口等生物活動劇烈的區(qū)域

    [43-45](如紅樹林、珊瑚礁和沼澤地等),顯著高于一些陸地水域、近岸河口區(qū)和大洋等開放水域.有機(jī)碳沉降通量和溶解有機(jī)碳交換通量表明研究海域底界面有機(jī)碳的降解所產(chǎn)生的溶解有機(jī)碳是水體DOC的重要來源.

    表6 不同海洋系統(tǒng)沉積物-水界面的溶解有機(jī)碳(DOC)的交換通量對比[g/(m2·a)]Table 6 Summary of DOC fluxes at the sedimentinterface for different marine systems [g/(m2·a)]

    3.4.2 有機(jī)碳的收支 根據(jù)2009年和2014年的調(diào)查資料和上述的討論,基于有機(jī)碳的物質(zhì)平衡模式,對夏季乳山灣鄰近海域水體和沉積物中有機(jī)碳的收支進(jìn)行了研究.以研究區(qū)域的水體和沉積物為核心,區(qū)域水深為 8~30m,面積為1800km2;根據(jù)表7中各界面流速和界面面積,計(jì)算研究區(qū)域水收支,水收支計(jì)算誤差為1%;再根據(jù)各界面水收支和有機(jī)碳濃度計(jì)算平流作用對區(qū)域有機(jī)碳的水平輸送的貢獻(xiàn).水體初級生產(chǎn)在夏季產(chǎn)生的POC和DOC的總量約為68Gg,底界面POC的沉積通量(FS)約為45Gg,DOC的交換通量為 24Gg,有機(jī)碳的埋藏?zé)o機(jī)碳的產(chǎn)出為21Gg.

    根據(jù)穩(wěn)定碳同位素?cái)?shù)據(jù)可知該區(qū)域陸源有機(jī)碳比例為47%,海源有機(jī)碳的比例為53%.為便于計(jì)算,將這個(gè)比例應(yīng)用到了輸入、輸出和沉積過程,這在相對狹小的區(qū)域具有一定的合理性.其它過程,如有機(jī)碳在水體和沉積物中的降解是根據(jù)物質(zhì)平衡計(jì)算得到.有機(jī)碳的收支結(jié)果顯示潮汐作用和沿岸流對研究區(qū)域有機(jī)碳的輸入(96Gg)和輸出(88Gg)的貢獻(xiàn)較大,是研究區(qū)域有機(jī)碳分布和向外海輸送的關(guān)鍵過程.從圖7還可以看出,水體有機(jī)碳的濃度受初級生產(chǎn)的影響,同時(shí)受底界面 DOC的釋放的影響;初級生產(chǎn)產(chǎn)生的有機(jī)碳向沉積物中埋藏的量約占初級生產(chǎn)的 35%;若最終埋藏的有機(jī)碳比例也符合該海域海陸源有機(jī)碳的比例的話,研究區(qū)域內(nèi)初級生產(chǎn)所固定的有機(jī)碳的最終埋藏效率約為 16%,最終保存在沉積物中的碳約占初級生產(chǎn)和外部輸入的 13%,這與桑溝灣海源有機(jī)碳的埋藏效率相近(20%)[46],但遠(yuǎn)高于世界大洋的平均值(0.3%)[47],可見近岸海域全球有機(jī)碳的主要沉積區(qū)域,高初級生產(chǎn)和沉積速率以及穩(wěn)定陸源有機(jī)碳的輸入是導(dǎo)致近岸有機(jī)碳得以有效保存的主要原因.另外,輸入到該區(qū)域的陸源有機(jī)碳有 19%在水柱內(nèi)被降解,大約53%的部分隨潮流向外海擴(kuò)散;同樣,海源有機(jī)碳中的大部分有機(jī)碳(54%)也隨潮流輸送到外海,表明潮流對物質(zhì)輸送的巨大作用.

    表7 收支區(qū)域邊界相關(guān)參數(shù)與有機(jī)碳交換通量Table 7 Parameters for fluxes calculation of organic carbon across the boundaries of the study area

    研究區(qū)域內(nèi)水量的輸入/輸出比為1.01,水量總體平衡,表明收支計(jì)算的合理性,有機(jī)碳的水平輸送主要是由濃度梯度導(dǎo)致的,且應(yīng)主要來自于乳山灣內(nèi)高濃度有機(jī)碳的貢獻(xiàn).盡管在收支模型中,潮流所攜帶的有機(jī)碳在量級上與初級生產(chǎn)固碳量大致相當(dāng)(圖7),但研究區(qū)域?qū)﹃懺从袡C(jī)碳的滯留和轉(zhuǎn)化作用較為有限(27%).值得說明的是潮汐和沿岸流所攜帶的有機(jī)碳中的一部分會在潮流等的作用下往復(fù)多次進(jìn)入研究區(qū)域,從而在表觀上顯著提高了潮流對有機(jī)碳的輸送作用,即潮流輸送的有機(jī)碳量顯著高于初級生產(chǎn)的固碳量.實(shí)際上,研究區(qū)域及鄰近水域初級生產(chǎn)的水平較為一致,潮流的搬運(yùn)作用變得相對有限,控制該區(qū)域有機(jī)物質(zhì)循環(huán)和埋藏的主要因素應(yīng)該還是水柱內(nèi)的初級生產(chǎn)和隨后向沉積物的埋藏作用.由于地下水對水體有機(jī)碳的貢獻(xiàn)研究較小,且貢獻(xiàn)不大[20],上述收支計(jì)算沒有考慮地下水對有機(jī)碳的貢獻(xiàn).

    水體自生有機(jī)碳的再生循環(huán)和降解(FR)主要為顆粒有機(jī)碳在沉降過程中的早期礦化降解,一部分為轉(zhuǎn)化為 DOC,另一部分則轉(zhuǎn)化為無機(jī)碳

    [36].沉積物中有機(jī)碳惰性部分得以埋藏保存,較活性部分經(jīng)歷早期礦化降解,其中部分有機(jī)碳降解為 DOC,并有一部分向上釋放到上覆水體(FE),另一部分則通過早期成巖作用以碳酸鹽形態(tài)保存下來,其余部分則可能生成甲烷與二氧化碳等.王娟等[48]研究表明夏季乳山灣沉積物中甲烷釋放的潛在速率平均值為 8g/(m2·a),接近于本研究中 DOC的釋放通量,可見甲烷等溫室氣體的釋放也是控制該海域沉積物中有機(jī)碳的保存的重要過程.

    值得一提的是,目前的收支模式考慮了初級生產(chǎn)參與下的碳循環(huán)過程,同時(shí)關(guān)注了陸源有機(jī)碳的貢獻(xiàn)和海水水平運(yùn)動的影響.不過,近年來日益加重的富營養(yǎng)化趨勢和環(huán)境變化(如低氧現(xiàn)象)可能會改變區(qū)域初級生產(chǎn)的水平,從而影響碳的循環(huán).受人類活動的影響,養(yǎng)殖海域有機(jī)碳的生物地球化學(xué)過程應(yīng)該有別于其它環(huán)境系統(tǒng),如碳的埋藏量與甲烷的釋放速率都比開放水域要高,同時(shí)也可能顯著高于歷史同期,而系統(tǒng)研究養(yǎng)殖區(qū)有機(jī)碳的收支對于深入研究全球碳循環(huán)具有重要的意義.沉積物-水界面DOC的交換會關(guān)系到沉積物中有機(jī)碳的收支與循環(huán)過程,相對較低的溶解氧水平(如本區(qū)域溶解氧水平較低,屬于新發(fā)現(xiàn)的近岸低氧區(qū)之一[14])還可能會增加 DOC的交換通量[49],減少碳的埋藏和有機(jī)碳的礦化速率[50].因此應(yīng)該關(guān)注低氧水域碳的循環(huán)過程,包 括溫室氣體的排放.

    圖7 夏季乳山灣鄰近海域水體和沉積物中有機(jī)碳的收支過程 (Gg,以碳計(jì))Fig.7 Organic carbon budget in the water and sediment in the coastal area of Rushan Bay in summerFP:初級生產(chǎn);FR-M:系統(tǒng)內(nèi)部海源有機(jī)碳的再生循環(huán)和降解; FS:有機(jī)碳沉積量,FS-M:海源有機(jī)碳沉積量(FS-M= FS×0.53),FS-T:陸源有機(jī)碳沉積量(FS-T= FS×0.47);FE:沉積物—水界面DOC交換通量;FB:有機(jī)碳的埋藏和無機(jī)碳的產(chǎn)出量;Fin:平流作用下有機(jī)碳的輸入(Fin= F1+ F3),其中Fin-M:為海源有機(jī)碳[Fin-M:=(Fin+ FP)×0.53- FP],Fin-T:為陸源有機(jī)碳(Fin-T=Fin- Fin-M);Fout:平流作用下有機(jī)碳的輸出(Fout= F2+ F4),其中Fout-M: 為海源有機(jī)碳(Fout-M= Fout×0.53),Fout-T:為陸源有機(jī)碳(Fout-T= Fout×0.47);FSGD:地下水的輸入,FC-T:系統(tǒng)外部凈輸入的陸源有機(jī)碳的循環(huán)過程

    4 結(jié)論

    4.1 2009年夏季乳山灣外鄰近海域DOC的含量范圍為 0.70~3.19mg/L,平均值為 1.80mg/L. DOC的分布呈現(xiàn)由近岸向外部海域逐漸降低的變化趨勢,灣口處DOC含量較高,明顯高于灣外;剖面分布上,在近岸區(qū)域隨著水深的增加 DOC的含量逐漸降低,而在外部海域 DOC的含量底層稍高于表層;時(shí)間變化上,8月份的DOC含量較其它月份高.2014年8月份DOC周日變化范圍為1.79~15.2mg/L,POC的周日變化范圍為0.04~1.33mg/L;表層水體有機(jī)碳和潮汐具有相似變化規(guī)律.

    4.2 調(diào)查區(qū)域 DOC的含量與分布受陸源輸入的影響,同時(shí)還受水體中浮游植物生長的影響.水體有機(jī)碳還受到潮汐的影響.

    4.3 研究區(qū)域顆粒有機(jī)碳沉降通量為(25± 0.8)g/(m2?a),約占初級生產(chǎn)固碳量的66%.DOC在沉積物—水界面存在向上覆水釋放的現(xiàn)象,DOC交換通量為14.4~97g/(m2·a),占水體存量的1.1%~13.4%.有機(jī)碳的收支表明研究海域底界面有機(jī)碳的降解所產(chǎn)生的溶解有機(jī)碳是水體DOC的重要來源,沉積物有機(jī)碳的埋藏和早期成巖過程是沉積物顆粒有機(jī)碳保存的重要途徑.

    4.4 最終保存在沉積物中的碳約占初級生產(chǎn)的13%.輸入到該區(qū)域的陸源有機(jī)碳有 19%在水柱內(nèi)被降解,超過 50%的海源和陸源有機(jī)碳隨潮流輸送到外海.沉積物-水界面有機(jī)碳的降解礦化過程會影響有機(jī)碳的循環(huán)過程,相對較低的溶解氧水平可能會增加 DOC的交換通量,從而減少碳的埋藏.

    [1] Falkowski P, Scholes R J, Boyle E, et al. The global carbon cycle: A test of our knowledge of earth as a system [J]. Science, 2000, 290(5490):291-296.

    [2] Duursma E K, Dawson R. Marine organic chemistry [M]. Elsevier Oceanography Series, 31, Elsevier Scientific Publishing Company, 1981.415-437.

    [3] Collins K J, Le P J, Williams B. An automated photochemical method for the determination of dissolved organic carbon in sea and estuarine waters [J]. Marine Chemistry, 1977,5(2):123-141.

    [4] Canuel E A, Hardison A K. Sources, Ages, and Alteration of Organic Matter in Estuaries [J]. Annual Review of Marine Science, 2016,8:409-434.

    [5] Bauer J E, Cai W J, Raymond P A, et al. The changing carbon cycle of the coastal ocean [J]. Nature, 2013,504(7478):61-70.

    [6] Tao S, Eglinton T I, Montlu?on D B, et al. Diverse origins and pre-depositional histories of organic matter in contemporary Chinese marginal sea sediments [J]. Geochimica et Cosmochimica Acta, 2016,191:70-88.

    [7] Chen Y, Yang G P, Wu G W, et al. Concentration and characterization of dissolved organic matter in the surface microlayer and subsurface water of the Bohai Sea, China [J]. Continental Shelf Research, 2013,52:97-107.

    [8] 李鴻妹,石曉勇,商容寧,等.秋季黃渤海溶解有機(jī)碳的分布特征及影響因素 [J]. 海洋環(huán)境科學(xué), 2013,32(2):161-164.

    [9] 劉 軍,于志剛,臧家業(yè),等.黃渤海有機(jī)碳的分布特征及收支評估研究 [J]. 地球科學(xué)進(jìn)展, 2015,30(5):564-578.

    [10] Ogawa H, Usui T, Koike I. Distribution of dissolved organic carbon in the East China Sea [J]. Deep Sea Research Part : ⅡTopical Studies in Oceanography, 2003,50(2):353-366.

    [11] 辛福言,陳碧鵑,曲克明,等.乳山灣表層海水 COD與氮,磷營養(yǎng)鹽的分布及其營養(yǎng)狀況 [J]. 海洋水產(chǎn)研究, 2004,25(5):52-56.

    [12] Liu D, Keesing J K, Dong Z, et al. Recurrence of the world’s largest green-tide in 2009in Yellow Sea, China: Porphyra yezoensis aquaculture rafts confirmed as nursery for macroalgal blooms [J]. Marine pollution bulletin, 2010,60(9):1423-1432.

    [13] 崔 毅,辛福言,馬紹賽,等.乳山灣沉積物重金屬污染及其生態(tài)危害評價(jià) [J]. 中國水產(chǎn)科學(xué), 2005,12(1):83-90.

    [14] 冉祥濱,臧家業(yè),韋欽勝,等.乳山灣口及其鄰近海域溶解氧分布特征及影響因素研究 [J]. 海洋學(xué)報(bào), 2011,33(4):173-180.

    [15] Liu J, Zang J Y, Zhao C Y, et al. Phosphorus speciation, transformation, and preservation in the coastal area of Rushan Bay [J]. Science of the Total Environment, 2016,565:258-270.

    [16] 高立蒙,姚 鵬,王金鵬,等.渤海表層沉積物中有機(jī)碳的分布和來源 [J]. 海洋學(xué)報(bào), 2016,38(6):8-20.

    [17] Behrenfeld M J, Falkowski P G. Photosynthetic rates derived from satellite-based chlorophyll concentration [J]. Limnology and oceanography, 1997,42(1):1-20.

    [18] Boudreau B P. Is burial velocity a master parameter for bioturbation? [J]. Geochimica et Cosmochimica Acta, 1994,58(4): 1243-1249.

    [19] Burdige D J, Berelson W M, Coale K H, et al. Fluxes of dissolved organic carbon from California continental margin sediments [J]. Geochimica et Cosmochimica Acta, 1999,63(10):1507-1515.

    [20] Hong Q Q, Cai P H, Shi X M, et al. Solute transport into the Jiulong River estuary via pore water exchange and submarine groundwater discharge: New insights from224Ra/228Th disequilibrium [J]. Geochimica et Cosmochimica Acta, 2016. Doi:http://dx.doi.org/10.1016/j.gca. 2016.11.002.

    [21] Lettmann K A, Riedinger N, Ramlau R, et al. Estimation of biogeochemical rates from concentration profiles: A novel inverse method [J]. Estuarine, Coastal and Shelf Science, 2012,100:26-37.

    [22] Ullman W J, Aller R C. Diffusion coefficients in nearshore marine sediments [J]. Limnology and Oceanography, 1982, 27(3):552-556.

    [23] 孫寶楠,尹訓(xùn)強(qiáng),連 展.乳山灣污染物擴(kuò)散模型數(shù)值實(shí)驗(yàn)與結(jié)果分析 [A]. 第十三屆全國水動力學(xué)學(xué)術(shù)會議暨第二十六屆全國水動力學(xué)研討會文集, 2014:601-608.

    [24] Wu J P, Calvert S E, Wong C S. Carbon and nitrogen isotope ratios in sedimenting particulate organic matter at an upwelling site off Vancouver Island [J]. Estuarine, Coastal and Shelf Science, 1999,48(2):193-203.

    [25] 馬紹賽,趙 俊,周詩賚,等.乳山灣水化學(xué)環(huán)境 [J]. 海洋水產(chǎn)研究, 1996,17(1):63-70.

    [26] Handson R B , Ducklow H W, Field J G, et al .The Changing Ocean Carbon Cycle: A midterm of the JGOFS [M]. London: Cambridge University Press, 2000.

    [27] 譚麗菊,張 哲,梁成菊,等.青島鄰近海域海水中有機(jī)碳的分布特征 [J]. 地球科學(xué)進(jìn)展, 2011,26(4):426-432.

    [28] 梁成菊,王江濤,譚麗菊.青島近海夏季溶解有機(jī)碳的分布 [J].中國海洋大學(xué)學(xué)報(bào):自然科學(xué)版, 2009(S1):217-221.

    [29] 冉祥濱,臧家業(yè),韋欽勝,等.乳山灣鄰近海域低氧現(xiàn)象及成因淺析 [J]. 海洋科學(xué)進(jìn)展, 2012,30(3):347-356.

    [30] 韋欽勝,魏修華,戰(zhàn) 閏,等.夏季南黃海西部溶解氧的分布特征及其影響因素分析 [J]. 海洋環(huán)境科學(xué), 2010,29(6):808-814.

    [31] 辛 明,馬德毅,王保棟,等.黃海溶解氧的平面分布特征及其季節(jié)變化 [J]. 中國海洋大學(xué)學(xué)報(bào)(自然科學(xué)版), 2013,7:011.

    [32] 梁成菊.青島近海有機(jī)碳的分布特征及影響因素 [D]. 中國海洋大學(xué), 2008.

    [33] 夏 斌,陳碧鵑,崔 毅,等.夏季桑溝灣養(yǎng)殖水域有機(jī)碳的平面分布特征及其來源分析 [J]. 漁業(yè)科學(xué)進(jìn)展, 2013(1):44-49.

    [34] 賀志鵬,宋金明,張乃星.南黃海溶解有機(jī)碳的生物地球化學(xué)特征分析 [J]. 海洋科學(xué)進(jìn)展, 2006,24(4):477-488.

    [35] 謝琳萍,王宗靈,王保棟,等.春季南黃海溶解有機(jī)碳的分布特征及其受控因素 [J]. 海洋環(huán)境科學(xué), 2010,29(5):636-640.

    [36] Hung J J, Wang S M, Chen Y L. Biogeochemical controls on distributions and fluxes of dissolved and particulate organic carbon in the Northern South China Sea [J]. Deep Sea Research Part II: Topical Studies in Oceanography, 2007,54(14):1486-1503.

    [37] Thomas C, Cauwet G, Minster J F. Dissolved organic carbon in the equatorial Atlantic Ocean [J]. Marine Chemistry, 1995,49(2): 155-169.

    [38] Wang D, Henrichs S M, Guo L. Distributions of nutrients, dissolved organic carbon and carbohydrates in the western Arctic Ocean [J]. Continental Shelf Research, 2006,26(14):1654-1667.

    [39] 倪建宇.赤道東北太平洋沉積物間隙水中溶解有機(jī)碳的分布特征 [J]. 海洋學(xué)報(bào), 2007,29(1):155-160.

    [40] Dadi T, V?lkner C, Koschorreck M. A sediment core incubation method to measure the flux of dissolved organic carbon between sediment and water [J]. Journal of Soils and Sediments, 2015, 15(12):2350-2358.

    [41] Martin W R, McCorkle D C. Dissolved organic carbon concentrations in marine pore waters determined by hightemperature oxidation [J]. Limnology and oceanography, 1993, 38(7):1464-1479.

    [42] Burdige D J, Alperin M J, Homstead J, et al. The role of benthic fluxes of dissolved organic carbon in oceanic and sedimentary carbon cycling [J]. Geophysical Research Letters, 1992,19(18): 1851-1854.

    [43] Dittmar T, Hertkorn N, Kattner G, et al. Mangroves, a major source of dissolved organic carbon to the oceans [J]. Global Biogeochemical Cycles, 2006,20,GB1012,doi:10.1029/2005GB 002570.

    [44] Nakajima R, Yoshida T, Azman B A R, et al. In situ release of coral mucus by Acropora and its influence on the heterotrophic bacteria [J]. Aquatic ecology, 2009,43(4):815-823.

    [45] Hassen M B. Nutrient and particulate matter annual budget. Interactions in a western European coastal salt marsh [J]. Comptes Rendus de l'Académie des Sciences-Series IIA-Earth and Planetary Science, 2000,330(6):429-435.

    [46] 劉 賽,楊 茜,楊 庶,等.桑溝灣養(yǎng)殖海域沉積物中碳埋藏通量的長期記錄 [J]. 海洋學(xué)報(bào), 2014,36(8):30-38.

    [47] Muller-Karger F E,Varela R,Thunell R,et al. The importance of continental margins in the global carbon cycle [J]. Geophysical Research Letters, 2005,32,L01602,doi:10.1029/2004GL021346.

    [48] 王 娟,曲克明,徐 勇,等.乳山灣養(yǎng)殖區(qū)甲烷釋放速率及其分布特征 [J]. 海洋水產(chǎn)研究, 2008,6:019.

    [49] Peter S, Isidorova A, Sobek S. Diffusion of dissolved organic carbon enhances carbon loss from anoxic lake sediment [J]. Journal of Geophysical Research: Biogeosciences, 2016, 2016JG003425.

    [50] 古小治,張啟超,孫淑雲(yún),等.富氧-缺氧過程對氧氣分布及交換過程影響 [J]. 中國環(huán)境科學(xué), 2015,35(5):1495-1501.

    致謝:本研究中POM潮汐潮流數(shù)值模式的結(jié)果由國家海洋局第一海洋研究所海洋環(huán)境科學(xué)與數(shù)值模擬國家海洋局重點(diǎn)實(shí)驗(yàn)室楊永增研究員和尹訓(xùn)強(qiáng)副研究員協(xié)助完成,在此表示感謝.

    Characteristics and benthic processes of organic carbon in the adjacent area of Rushan Bay.

    ZANG Jia-ye1, ZHAO Chen-ying1, LIU Jun1,2, XIE Lin-ping1, WANG Yi-bin1, ZHANG Ai-jun1, RAN Xiang-bin1*(1.Research Center for Marine Ecology, First Institute of Oceanography, State Oceanic Administration, Qingdao 266061, China;2.Laboratory of Marine Geology and Geophysics, First Institute of Oceanography, State Oceanic Administration, Qingdao 266061, China). China Environmental Science, 2017,37(3):1089~1102

    Estuaries, lagoons and bays are considered as the most important ocean margins due to their important influence in land-sea interaction under human activities, playing key roles in the regional carbon biogeochemical cycles. Based on measurements obtained from comprehensive survey in the adjacent area of Rushan Bay in summer of 2009 and 2014, organic carbon (OC) distributions, its benthic processes and influence factors were determined. The results showed that the concentration of dissolved organic carbon (DOC) in summer 2009 ranges from 0.70 to 3.19mg/L, with an average concentration of 1.80mg/L. Concentration of DOC was high in August, and low in June. Time-dependent variation of DOC and particulate organic carbon (POC) were ranged in1.79~15.2mg/L and 0.04~1.33mg/L, respectively. DOC patterns were mainly affected by human activities, primary production and tide. Estimated POC sedimentation rate was (25±0.8) g/m2in summer, accounting for 66% of the carbon fixation by phytoplankton related to primary production. DOC concentrations in the pore water of the upper sediment (0~4cm) were 8~9times larger than those in the overlying water, and benthic DOC flux at the sediment-water interface ranged from 14.4 to 97g/(m2a), accounting for 1.1%~13.4% of DOC in the water DOC pool. Anthropogenic activities were the key factors that influence the OC distribution and cycling in the coastal areas of Rushan Bay. The OC budget shows that regeneration of sedimentary TOC contributes to a large share of diffusive DOC flux, indicating that the sediment is a significant source of DOC to water column in the study area. More than 50% of OC can be transported by tides and currents. The net burial of OC in sediments represents 13% of the carbonpool that yield by the primary production in summer. Low-oxygen level may enhance benthic diffusive of DOC and other diagenetic processes, which would play a key role in controlling OC conservation in the sediment.

    Rushan Bay;organic carbon;flux;primary production;benthic process

    X55,P72

    A

    1000-6923(2017)03-1089-14

    臧家業(yè)(1962- ),男,山東青島人,研究員,碩士,主要從事海洋環(huán)境科學(xué)研究,發(fā)表論文30余篇.

    2016-07-20

    中央級公益性科研院所基本科研業(yè)務(wù)費(fèi)專項(xiàng)(GY0216Q03);國家自然科學(xué)基金(41376093,41106072)

    * 責(zé)任作者, 副研究員, rxb@fio.org.cn

    猜你喜歡
    乳山表層沉積物
    文化經(jīng)濟(jì)融合視域下乳山民間文學(xué)的傳承與發(fā)展研究
    晚更新世以來南黃海陸架沉積物源分析
    山東乳山:草莓小鎮(zhèn)的“莓”好生活
    渤海油田某FPSO污水艙沉積物的分散處理
    海洋石油(2021年3期)2021-11-05 07:43:12
    半潛式平臺表層卡套管處理與認(rèn)識
    海洋石油(2021年3期)2021-11-05 07:43:10
    水體表層沉積物對磷的吸收及釋放研究進(jìn)展
    乳山長牡蠣(Crassostrea gigas)的抗性基因表達(dá)和生存環(huán)境的季節(jié)差異
    討論用ICP-AES測定土壤和沉積物時(shí)鈦對鈷的干擾
    氬弧熔覆原位合成Ti(C,N)-WC增強(qiáng)鎳基表層復(fù)合材料的研究
    焊接(2015年6期)2015-07-18 11:02:25
    超聲波光整強(qiáng)化40Cr表層顯微硬度研究
    免费观看a级毛片全部| 老司机在亚洲福利影院| 国产在线观看jvid| 18禁国产床啪视频网站| 51午夜福利影视在线观看| 久9热在线精品视频| 中文字幕最新亚洲高清| 麻豆乱淫一区二区| 国产精品久久久人人做人人爽| 性色av乱码一区二区三区2| 免费在线观看亚洲国产| 一级毛片高清免费大全| 人人妻,人人澡人人爽秒播| 亚洲熟妇中文字幕五十中出 | 国产精品久久视频播放| 色尼玛亚洲综合影院| 久久精品亚洲熟妇少妇任你| 精品一区二区三卡| 精品福利永久在线观看| 欧美日韩亚洲国产一区二区在线观看 | 日韩免费高清中文字幕av| 精品国产超薄肉色丝袜足j| 啦啦啦免费观看视频1| 国产主播在线观看一区二区| 99国产精品免费福利视频| 叶爱在线成人免费视频播放| 满18在线观看网站| cao死你这个sao货| 欧美日韩成人在线一区二区| 欧美日韩福利视频一区二区| 国产成人一区二区三区免费视频网站| 亚洲第一欧美日韩一区二区三区| 国产成人欧美在线观看 | 国产一区在线观看成人免费| 久热这里只有精品99| 中文字幕人妻丝袜一区二区| 一个人免费在线观看的高清视频| 国产成人精品无人区| 香蕉国产在线看| 欧美黄色淫秽网站| 丁香六月欧美| 亚洲欧美日韩另类电影网站| 丝瓜视频免费看黄片| 久久精品熟女亚洲av麻豆精品| 亚洲第一青青草原| 在线免费观看的www视频| 亚洲av日韩精品久久久久久密| 国产成人精品在线电影| 亚洲专区中文字幕在线| 国产精品香港三级国产av潘金莲| 欧美日韩乱码在线| 精品国产乱码久久久久久男人| 视频在线观看一区二区三区| 日韩一卡2卡3卡4卡2021年| 在线播放国产精品三级| 日日夜夜操网爽| 日韩免费高清中文字幕av| 亚洲五月色婷婷综合| 老熟女久久久| 日韩欧美一区二区三区在线观看 | 丝袜人妻中文字幕| 很黄的视频免费| 亚洲欧洲精品一区二区精品久久久| 精品高清国产在线一区| 久久国产精品影院| 亚洲精品在线观看二区| 国产不卡av网站在线观看| 色94色欧美一区二区| 亚洲精品自拍成人| 日韩一卡2卡3卡4卡2021年| 人人澡人人妻人| 国产精品电影一区二区三区 | 人成视频在线观看免费观看| 欧美人与性动交α欧美精品济南到| 成年女人毛片免费观看观看9 | av线在线观看网站| 国产蜜桃级精品一区二区三区 | 一级毛片高清免费大全| 亚洲专区国产一区二区| 午夜免费成人在线视频| 欧美老熟妇乱子伦牲交| 亚洲一码二码三码区别大吗| 久久九九热精品免费| 天堂中文最新版在线下载| 另类亚洲欧美激情| 黄片小视频在线播放| 国产欧美日韩精品亚洲av| 俄罗斯特黄特色一大片| 亚洲国产欧美一区二区综合| 又大又爽又粗| 热99re8久久精品国产| 婷婷丁香在线五月| 欧美最黄视频在线播放免费 | av片东京热男人的天堂| svipshipincom国产片| 一级毛片高清免费大全| 国产精品一区二区精品视频观看| 亚洲国产欧美日韩在线播放| 飞空精品影院首页| 啦啦啦免费观看视频1| 日本精品一区二区三区蜜桃| 精品国产乱码久久久久久男人| 国产精品成人在线| 老司机午夜十八禁免费视频| 国产精品 欧美亚洲| 欧美日韩国产mv在线观看视频| 国产精品乱码一区二三区的特点 | 欧美精品av麻豆av| 黄色a级毛片大全视频| netflix在线观看网站| 一级a爱视频在线免费观看| 大型黄色视频在线免费观看| 久久精品aⅴ一区二区三区四区| 五月开心婷婷网| 电影成人av| 精品一区二区三区四区五区乱码| 国产精品欧美亚洲77777| 少妇猛男粗大的猛烈进出视频| 一级黄色大片毛片| 欧美av亚洲av综合av国产av| 99热只有精品国产| 亚洲欧美一区二区三区黑人| 麻豆国产av国片精品| 亚洲avbb在线观看| 丰满迷人的少妇在线观看| 国产免费男女视频| 午夜视频精品福利| 国产av一区二区精品久久| 狠狠婷婷综合久久久久久88av| 后天国语完整版免费观看| 欧美日韩亚洲国产一区二区在线观看 | 久久久久视频综合| 淫妇啪啪啪对白视频| 又黄又爽又免费观看的视频| 日韩免费av在线播放| 久久精品熟女亚洲av麻豆精品| 捣出白浆h1v1| 大型av网站在线播放| 如日韩欧美国产精品一区二区三区| 欧美另类亚洲清纯唯美| 搡老熟女国产l中国老女人| 亚洲三区欧美一区| 亚洲国产看品久久| 午夜福利在线免费观看网站| 中亚洲国语对白在线视频| 精品无人区乱码1区二区| 怎么达到女性高潮| 又紧又爽又黄一区二区| 涩涩av久久男人的天堂| 日韩三级视频一区二区三区| 精品福利永久在线观看| 女人精品久久久久毛片| 免费一级毛片在线播放高清视频 | 精品国产一区二区久久| www.熟女人妻精品国产| 欧美 亚洲 国产 日韩一| 999精品在线视频| av视频免费观看在线观看| 亚洲av电影在线进入| 国产精华一区二区三区| 99久久国产精品久久久| 高潮久久久久久久久久久不卡| 国产乱人伦免费视频| 亚洲精品在线观看二区| 18禁国产床啪视频网站| 满18在线观看网站| 在线av久久热| 免费在线观看日本一区| 亚洲精品成人av观看孕妇| 在线视频色国产色| 亚洲熟妇熟女久久| 欧美日韩成人在线一区二区| 成人三级做爰电影| 午夜91福利影院| 999久久久国产精品视频| 一本一本久久a久久精品综合妖精| 亚洲国产欧美一区二区综合| 欧美丝袜亚洲另类 | 欧洲精品卡2卡3卡4卡5卡区| 女人爽到高潮嗷嗷叫在线视频| 飞空精品影院首页| 亚洲成人国产一区在线观看| 天天操日日干夜夜撸| 精品久久久久久电影网| 99国产精品99久久久久| 999久久久国产精品视频| 国产精品免费一区二区三区在线 | 精品欧美一区二区三区在线| 午夜福利影视在线免费观看| 免费看a级黄色片| 久久久水蜜桃国产精品网| 久久精品亚洲av国产电影网| 国产在线精品亚洲第一网站| 激情视频va一区二区三区| 在线观看免费视频网站a站| cao死你这个sao货| 自拍欧美九色日韩亚洲蝌蚪91| 一区二区三区国产精品乱码| 18禁观看日本| 99在线人妻在线中文字幕 | 亚洲精品在线美女| 国产欧美日韩一区二区三| 亚洲黑人精品在线| 人妻久久中文字幕网| 亚洲va日本ⅴa欧美va伊人久久| 涩涩av久久男人的天堂| 免费观看人在逋| 免费不卡黄色视频| 91在线观看av| 久久久久久人人人人人| 在线天堂中文资源库| 桃红色精品国产亚洲av| 正在播放国产对白刺激| 成年版毛片免费区| 成年动漫av网址| 狠狠婷婷综合久久久久久88av| 国产欧美日韩综合在线一区二区| 久久人人爽av亚洲精品天堂| 精品一区二区三区视频在线观看免费 | 黄频高清免费视频| 色婷婷av一区二区三区视频| 亚洲自偷自拍图片 自拍| 亚洲精品久久成人aⅴ小说| 这个男人来自地球电影免费观看| 丰满迷人的少妇在线观看| 午夜精品国产一区二区电影| 久久国产精品大桥未久av| 欧美成人免费av一区二区三区 | 国产欧美日韩精品亚洲av| 一区在线观看完整版| 男男h啪啪无遮挡| 无遮挡黄片免费观看| 亚洲精品中文字幕一二三四区| 久久草成人影院| 亚洲va日本ⅴa欧美va伊人久久| av在线播放免费不卡| 叶爱在线成人免费视频播放| 国产成人欧美在线观看 | 欧美日韩亚洲综合一区二区三区_| bbb黄色大片| 91在线观看av| 国产男靠女视频免费网站| 人人妻人人澡人人爽人人夜夜| 日本wwww免费看| 在线观看免费高清a一片| 精品国产亚洲在线| 亚洲久久久国产精品| 制服诱惑二区| 夜夜躁狠狠躁天天躁| 老汉色∧v一级毛片| 久久久国产精品麻豆| 中文字幕色久视频| 老熟妇乱子伦视频在线观看| 丁香欧美五月| 亚洲av第一区精品v没综合| 成人三级做爰电影| 国产精品1区2区在线观看. | 成人手机av| 免费观看精品视频网站| 天天躁夜夜躁狠狠躁躁| 色94色欧美一区二区| 国产区一区二久久| 国产免费男女视频| 99re在线观看精品视频| 男女下面插进去视频免费观看| 免费久久久久久久精品成人欧美视频| 水蜜桃什么品种好| 人妻一区二区av| 操出白浆在线播放| 一本综合久久免费| 亚洲国产看品久久| 亚洲五月婷婷丁香| 美女福利国产在线| 嫩草影视91久久| 18禁美女被吸乳视频| 国产精品免费大片| 侵犯人妻中文字幕一二三四区| 亚洲精品美女久久av网站| 久久这里只有精品19| 搡老熟女国产l中国老女人| 精品一区二区三卡| 最新的欧美精品一区二区| 欧美最黄视频在线播放免费 | 色老头精品视频在线观看| 久久亚洲精品不卡| 久久九九热精品免费| 亚洲国产欧美网| 免费黄频网站在线观看国产| 国产精品 国内视频| 欧美精品av麻豆av| 成人特级黄色片久久久久久久| 国产精品.久久久| 大码成人一级视频| 国产精品1区2区在线观看. | 日本wwww免费看| 又黄又爽又免费观看的视频| 国产不卡av网站在线观看| 麻豆国产av国片精品| 亚洲一区二区三区不卡视频| 国产淫语在线视频| 丝袜人妻中文字幕| xxx96com| 韩国精品一区二区三区| 精品免费久久久久久久清纯 | 精品一品国产午夜福利视频| 嫁个100分男人电影在线观看| 黄色毛片三级朝国网站| tube8黄色片| 电影成人av| 亚洲欧美日韩高清在线视频| av线在线观看网站| 交换朋友夫妻互换小说| 亚洲在线自拍视频| 三上悠亚av全集在线观看| 91av网站免费观看| av一本久久久久| 亚洲五月色婷婷综合| www.999成人在线观看| 如日韩欧美国产精品一区二区三区| 亚洲午夜理论影院| 村上凉子中文字幕在线| 视频在线观看一区二区三区| 国精品久久久久久国模美| 美女午夜性视频免费| 久久精品亚洲av国产电影网| 亚洲专区国产一区二区| 欧美日韩瑟瑟在线播放| 丰满的人妻完整版| 99精品久久久久人妻精品| 超碰成人久久| 啦啦啦 在线观看视频| 国产成人欧美在线观看 | 最近最新免费中文字幕在线| 亚洲精品久久午夜乱码| 免费少妇av软件| 麻豆国产av国片精品| 久久精品亚洲熟妇少妇任你| 国产又爽黄色视频| 欧美日韩成人在线一区二区| 在线国产一区二区在线| av在线播放免费不卡| 男人的好看免费观看在线视频 | 高潮久久久久久久久久久不卡| 国产一区二区三区视频了| 日韩视频一区二区在线观看| 精品久久蜜臀av无| 精品国产亚洲在线| 91av网站免费观看| 丰满迷人的少妇在线观看| 国产在线一区二区三区精| 夜夜爽天天搞| 精品少妇久久久久久888优播| 国产日韩欧美亚洲二区| 最新的欧美精品一区二区| 热re99久久国产66热| 99re在线观看精品视频| 亚洲精华国产精华精| 日日夜夜操网爽| 精品一区二区三卡| 麻豆国产av国片精品| 亚洲一区中文字幕在线| 成年人免费黄色播放视频| 欧美日韩亚洲国产一区二区在线观看 | 麻豆成人av在线观看| 99国产精品免费福利视频| av片东京热男人的天堂| 又紧又爽又黄一区二区| 日韩熟女老妇一区二区性免费视频| 久久久国产一区二区| 丝袜在线中文字幕| 久久99一区二区三区| 99精品在免费线老司机午夜| 一级a爱片免费观看的视频| 香蕉国产在线看| 欧美日韩瑟瑟在线播放| 国产一区二区三区综合在线观看| 久久久久国产一级毛片高清牌| 亚洲精品国产区一区二| 狠狠狠狠99中文字幕| 亚洲精品久久午夜乱码| 午夜福利乱码中文字幕| 视频区图区小说| 在线观看免费视频日本深夜| 99精国产麻豆久久婷婷| 国产高清国产精品国产三级| 自线自在国产av| 一级,二级,三级黄色视频| 午夜91福利影院| 黄片小视频在线播放| 狂野欧美激情性xxxx| 国产一区二区三区视频了| 男女之事视频高清在线观看| 人人妻人人添人人爽欧美一区卜| 免费黄频网站在线观看国产| 成人手机av| 免费日韩欧美在线观看| 欧美黑人欧美精品刺激| 老司机福利观看| 嫩草影视91久久| 久久九九热精品免费| 老熟妇乱子伦视频在线观看| 国产亚洲精品久久久久5区| 亚洲在线自拍视频| 精品一区二区三卡| 一进一出抽搐动态| 欧美 日韩 精品 国产| 日韩大码丰满熟妇| 亚洲专区字幕在线| 五月开心婷婷网| 久久久久久久久免费视频了| 天天影视国产精品| 别揉我奶头~嗯~啊~动态视频| 亚洲欧美精品综合一区二区三区| 国产精品av久久久久免费| 热re99久久国产66热| 亚洲av美国av| 欧美另类亚洲清纯唯美| 国产高清视频在线播放一区| 久久天躁狠狠躁夜夜2o2o| 不卡av一区二区三区| 国产欧美日韩精品亚洲av| 91大片在线观看| 国产淫语在线视频| 免费不卡黄色视频| 三上悠亚av全集在线观看| 日韩大码丰满熟妇| 不卡av一区二区三区| 亚洲性夜色夜夜综合| 亚洲精品av麻豆狂野| 人人妻人人爽人人添夜夜欢视频| www.精华液| videos熟女内射| x7x7x7水蜜桃| 每晚都被弄得嗷嗷叫到高潮| 亚洲aⅴ乱码一区二区在线播放 | 免费在线观看影片大全网站| 久久久久久久久久久久大奶| 成人av一区二区三区在线看| 久久亚洲精品不卡| 国产区一区二久久| 中文亚洲av片在线观看爽 | 韩国av一区二区三区四区| 久久久精品免费免费高清| 日本vs欧美在线观看视频| 欧美 亚洲 国产 日韩一| 成年人免费黄色播放视频| 丰满迷人的少妇在线观看| tube8黄色片| 一级作爱视频免费观看| 亚洲专区国产一区二区| 18禁美女被吸乳视频| 99久久99久久久精品蜜桃| 欧美黑人欧美精品刺激| 免费av中文字幕在线| 日本欧美视频一区| 一进一出抽搐动态| 麻豆成人av在线观看| 国产伦人伦偷精品视频| 91精品国产国语对白视频| 日韩三级视频一区二区三区| 久久ye,这里只有精品| 久久国产精品男人的天堂亚洲| 女警被强在线播放| 999久久久精品免费观看国产| a级毛片在线看网站| 美国免费a级毛片| 欧美激情极品国产一区二区三区| 人妻久久中文字幕网| 伊人久久大香线蕉亚洲五| 波多野结衣一区麻豆| 黄色成人免费大全| 亚洲欧美一区二区三区黑人| 很黄的视频免费| 午夜视频精品福利| 成年人黄色毛片网站| 一边摸一边做爽爽视频免费| 男女下面插进去视频免费观看| 国产在视频线精品| 午夜福利免费观看在线| 最近最新中文字幕大全电影3 | av线在线观看网站| 一区在线观看完整版| 久久久精品区二区三区| 色精品久久人妻99蜜桃| 成年人午夜在线观看视频| 亚洲国产精品一区二区三区在线| 欧美日韩国产mv在线观看视频| 极品教师在线免费播放| 久久国产精品人妻蜜桃| 国产高清激情床上av| av网站在线播放免费| 日本欧美视频一区| 欧美黄色淫秽网站| 99国产极品粉嫩在线观看| 国产色视频综合| 久久香蕉精品热| 久久性视频一级片| 国产精品免费视频内射| 久久久久久久精品吃奶| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲免费av在线视频| 国产精品影院久久| 国产麻豆69| 国产日韩一区二区三区精品不卡| 日韩精品免费视频一区二区三区| 精品久久蜜臀av无| 丁香欧美五月| 操美女的视频在线观看| 婷婷精品国产亚洲av在线 | ponron亚洲| 久久久国产成人免费| 丝袜在线中文字幕| 欧美亚洲 丝袜 人妻 在线| 一二三四在线观看免费中文在| 满18在线观看网站| www日本在线高清视频| 精品国产乱码久久久久久男人| 美国免费a级毛片| 亚洲性夜色夜夜综合| 精品卡一卡二卡四卡免费| 大香蕉久久网| 真人做人爱边吃奶动态| av有码第一页| 色综合婷婷激情| 成人免费观看视频高清| 欧美最黄视频在线播放免费 | 一区福利在线观看| 国产精品欧美亚洲77777| 国产免费av片在线观看野外av| 午夜成年电影在线免费观看| 国产主播在线观看一区二区| 熟女少妇亚洲综合色aaa.| 国产不卡一卡二| 美女扒开内裤让男人捅视频| 国产又色又爽无遮挡免费看| 亚洲avbb在线观看| av在线播放免费不卡| 久久久国产成人精品二区 | 久久亚洲精品不卡| 国产精品 欧美亚洲| 久久久精品免费免费高清| 欧美国产精品一级二级三级| 久热这里只有精品99| 男女高潮啪啪啪动态图| 另类亚洲欧美激情| 国产人伦9x9x在线观看| 欧美黑人欧美精品刺激| 中文字幕av电影在线播放| videosex国产| 一本大道久久a久久精品| cao死你这个sao货| 一级a爱视频在线免费观看| 十八禁人妻一区二区| 黄色 视频免费看| 人人妻人人澡人人看| 国产男靠女视频免费网站| 欧美一级毛片孕妇| 亚洲精品久久成人aⅴ小说| 高清在线国产一区| 99精品欧美一区二区三区四区| 丁香六月欧美| 黑丝袜美女国产一区| 欧美 亚洲 国产 日韩一| 国产精品久久久久久人妻精品电影| 国产精品1区2区在线观看. | 国精品久久久久久国模美| 成人亚洲精品一区在线观看| 高清欧美精品videossex| 人成视频在线观看免费观看| 国产精品免费视频内射| 成年人黄色毛片网站| 纯流量卡能插随身wifi吗| 亚洲欧美激情综合另类| 免费黄频网站在线观看国产| 免费久久久久久久精品成人欧美视频| 国产精品免费一区二区三区在线 | 久久久久久久国产电影| 国产成人av教育| tube8黄色片| 精品午夜福利视频在线观看一区| 国产蜜桃级精品一区二区三区 | 亚洲成av片中文字幕在线观看| 91精品国产国语对白视频| 人人妻人人澡人人爽人人夜夜| netflix在线观看网站| 久久中文字幕人妻熟女| 亚洲va日本ⅴa欧美va伊人久久| 一级片免费观看大全| 国产精品久久久久久精品古装| 久久这里只有精品19| 国产不卡一卡二| 最近最新中文字幕大全电影3 | 国产午夜精品久久久久久| 男男h啪啪无遮挡| 久久99一区二区三区| 女同久久另类99精品国产91| 在线观看日韩欧美| 亚洲国产精品sss在线观看 | 99热网站在线观看| 国产又爽黄色视频| 久久影院123| 久久精品人人爽人人爽视色| 91在线观看av| 亚洲一码二码三码区别大吗| 国产精品亚洲一级av第二区| 69精品国产乱码久久久| 日韩制服丝袜自拍偷拍| 婷婷丁香在线五月| 男女床上黄色一级片免费看| 亚洲av成人一区二区三| 男女之事视频高清在线观看| 国产亚洲精品久久久久久毛片 | 999久久久国产精品视频| 少妇的丰满在线观看| 国产精品国产av在线观看|