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

    東海及黃海南部海域走航海-氣CO2通量季節(jié)性變化*

    2016-12-13 01:18:55胡登輝高郭平翟惟東張春玲
    海洋科學(xué)進(jìn)展 2016年4期
    關(guān)鍵詞:海表黃海東海

    胡登輝,高郭平*,翟惟東,張春玲

    (1.上海海洋大學(xué) 海洋科學(xué)學(xué)院,上海 201306;2.大洋漁業(yè)資源可持續(xù)開發(fā)省部共建教育部重點(diǎn)實(shí)驗(yàn),上海 201306;3.國家海洋環(huán)境監(jiān)測(cè)中心,遼寧 大連 116023)

    ?

    東海及黃海南部海域走航海-氣CO2通量季節(jié)性變化*

    胡登輝1,2,高郭平1,2*,翟惟東3,張春玲1,2

    (1.上海海洋大學(xué) 海洋科學(xué)學(xué)院,上海 201306;2.大洋漁業(yè)資源可持續(xù)開發(fā)省部共建教育部重點(diǎn)實(shí)驗(yàn),上海 201306;3.國家海洋環(huán)境監(jiān)測(cè)中心,遼寧 大連 116023)

    根據(jù)2003—2008年東海及黃海南部海域多個(gè)航次現(xiàn)場(chǎng)觀測(cè)獲得的海表溫度、鹽度及海水表層pCO2觀測(cè)數(shù)據(jù),分析該海域海水表層pCO2及海-氣CO2通量的季節(jié)變化特征,探討了海-氣界面CO2轉(zhuǎn)移與海表溫度、鹽度分布之間的聯(lián)系。結(jié)果表明:該海域海水表層pCO2及海-氣CO2通量具有顯著的季節(jié)性差異。近海區(qū)域,春季受海表溫度上升、生物作用加強(qiáng)的影響黃海南部、東海近岸區(qū)及陸架中部、東海南部表現(xiàn)為大氣CO2匯,其海-氣界面季平均通量分別為(-7.77±6.59),(-11.08±8.99),(-2.94±6.78) mmol·m-2·d-1。夏季黃海南部區(qū)域表現(xiàn)為大氣CO2源(2.99±6.09) mmol·m-2·d-1,與該海域的下層海水涌升有關(guān),東海中部陸架區(qū)及東海南部近岸區(qū)由于淡水輸入,形成躍層阻礙水體混合,再加上光合作用增強(qiáng)等的綜合作用為大氣CO2匯,通量為(-4.81±8.92),(-0.75±12.14) mmol·m-2·d-1。秋季北風(fēng)逐漸增強(qiáng)水體混合加劇,向冬季格局轉(zhuǎn)變,底層富含CO2的海水上涌,致使海表pCO2升高,整個(gè)海區(qū)表現(xiàn)為大氣CO2源。在年際變化上,春季碳匯呈減弱趨勢(shì),而秋季碳源則逐漸增強(qiáng)。

    海-氣CO2通量; 二氧化碳分壓; 季節(jié)性變化; 影響因素

    全球氣候變化背景下,海洋與大氣之間氣體交換通量,特別是對(duì)二氧化碳通量的監(jiān)測(cè)、估算,對(duì)深刻理解生物地球化學(xué)循環(huán)以及全球氣候變遷具有重大意義[1]。但迄今為止,全球海洋碳循環(huán)研究仍存在諸多不確定性,其中之一就是陸架邊緣海碳通量及其控制因素不確定性問題,并逐漸被人們所重視[2]。東海及黃海南部是西北太平洋以陸架區(qū)為特征的邊緣海,具有較為復(fù)雜的物理、化學(xué)和生物過程,長江攜帶的大量淡水和陸源無機(jī)物注入該海域,并顯著影響該海域的物理、化學(xué)、生物過程[3]。由于長江沖淡水、沿岸流、臺(tái)灣暖流的相互作用,黑潮的持續(xù)入侵以及東亞季風(fēng)等的影響,使該海域呈現(xiàn)生產(chǎn)力高、環(huán)流復(fù)雜、陸海作用強(qiáng)烈的特點(diǎn),是全球物質(zhì)循環(huán)“源-匯”效應(yīng)的過渡區(qū)。

    東海及南海海域CO2源/匯分布及其通量得到較多的關(guān)注和研究。國外,Cai等[4]和Borges等[5]對(duì)全球邊緣海海-氣CO2通量進(jìn)行了研究,認(rèn)為全球中、高緯度的陸架邊緣??傮w表現(xiàn)為大氣CO2匯,而低緯則表現(xiàn)為大氣CO2源。Tsunogai等[6]基于PN斷面提出大陸架碳泵機(jī)制—冬季高密度水團(tuán)對(duì)CO2的輸運(yùn)。Shim等在長江口北部海水表層pCO2研究認(rèn)為,在近岸區(qū)主要受海水混合及上升流影響,而外海則受對(duì)馬暖流影響較顯著[7]。Oliveira等通過數(shù)值模型得出淡水輸入、表層流等因素對(duì)塔霍河河口區(qū)域內(nèi)冬季pCO2分布有很大影響[8]。以上在對(duì)近岸區(qū)及深水區(qū)海水表層pCO2影響因素直接的差異可以得出,近岸區(qū)主要由生物作用和陸源因素的影響;外海,除上升流區(qū)外,海水混合對(duì)海水表層pCO2起重要影響[9]。國內(nèi),在東海及黃海海域,張龍軍、宋金明、翟惟東等對(duì)東海、黃海海域內(nèi)化學(xué)及生物作用對(duì)近海海-氣pCO2的影響機(jī)制也進(jìn)行了比較全面的研究,得出夏季受生物和長江沖淡水影響較大,冬季則受海水混合影響較顯著[10-14]。由于近海海域生物化學(xué)作用強(qiáng)烈,其對(duì)海-氣CO2通量也有著決定性影響。翟惟東等通過對(duì)2005—2008年長江口河口及外海海水表層pCO2同溶解氧間年際變化得出,該海域較暖季節(jié)(4—10月)在初級(jí)生產(chǎn)力較強(qiáng)區(qū)域出現(xiàn)pCO2降低,溶解氧增強(qiáng),為大氣CO2匯。其他月份則主要受海水混合影響,為大氣CO2源。并且長江口外海-氣CO2年平均通量約為長江口北部2倍,為(-1.9±1.3) mol·m-2·a-1[15]。在此基礎(chǔ)上,對(duì)黃海南部東海北部碳酸鹽與海水表層pCO2間相關(guān)關(guān)系進(jìn)行研究,并提出碳酸鹽體系能夠減緩CO2達(dá)到平衡的過程[16]。以上研究充分說明陸架邊緣海區(qū)CO2源、匯問題的復(fù)雜性。但目前大多聚集于較小區(qū)域的單一季節(jié)性研究,無法從大的范圍內(nèi)系統(tǒng)地研究多種因素對(duì)海-氣CO2通量的變化及其影響。

    本研究在多年東海及黃海南部海域走航海-氣CO2觀測(cè)資料的基礎(chǔ)上,系統(tǒng)地研究該海域CO2源/匯分布問題,并對(duì)該海域走航斷面上海-氣CO2通量的季節(jié)性變化及其控制因素進(jìn)行分析,以及其通量年際變化進(jìn)行初步探討。為進(jìn)一步研究陸架邊緣海在全球碳循環(huán)系統(tǒng)中對(duì)CO2吸收、轉(zhuǎn)移的物理、生物和化學(xué)機(jī)制提供基礎(chǔ)。

    1 資料來源及計(jì)算方法

    1.1 研究區(qū)域劃分及觀測(cè)數(shù)據(jù)

    整個(gè)調(diào)查區(qū)域覆蓋60 m水深附近的大部分東海海域(約70%)及部分黃海海域(約50%)。由于陸架邊緣海的CO2體系受到陸地、海洋等多種生物地球化學(xué)過程的共同作用,具有顯著的區(qū)域性差異,如北部的黃海暖流(主要在冬季)、蘇北沿岸流(主要在夏季),中部的長江沖淡水以及南部的沿岸流和臺(tái)灣暖流等。故根據(jù)水文及氣候特點(diǎn),將研究海域分為黃海南部海域(Ⅰ區(qū),范圍:120°00′~ 124°00′E , 32°30′~35°00′N)、東海近岸及中部陸架區(qū)域(Ⅱ區(qū),范圍:120°00′~124°00′E , 30°00′~32°30′N)、東海南部海域(Ⅲ區(qū),范圍:120°00′~124°00′E, 26°00′~30°00′N)三個(gè)局部區(qū)域(圖1),進(jìn)而分別計(jì)算各個(gè)區(qū)域的海-氣CO2通量,并分析其季節(jié)變化特征。

    圖1 研究海區(qū)走航觀測(cè)航線和調(diào)查站位Fig.1 Cruise routes and survey stations

    文中所使用的觀測(cè)數(shù)據(jù)由廈門大學(xué)以及國家海洋環(huán)境中心提供,包括2003—2008年的多個(gè)航次收集的走航調(diào)查數(shù)據(jù)(圖1a和表1)以及2006-09—10的站位調(diào)查數(shù)據(jù)(圖1b)。調(diào)查航次觀測(cè)數(shù)據(jù)包括海表溫度、鹽度、海水表層CO2摩爾分?jǐn)?shù)(xCO2)、大氣CO2摩爾分?jǐn)?shù)(xCO2)、海表大氣壓、風(fēng)速等多個(gè)觀測(cè)要素。但數(shù)據(jù)主要集中在春季、夏季和秋季。海水表層pCO2由干空氣下海水表層CO2摩爾分?jǐn)?shù)、海表大氣壓、水汽平衡器內(nèi)飽和水汽壓計(jì)算獲得[17]:

    pCO2(干空氣下)=(P-VP(H2O,s/w))×xCO2,

    (1)

    式中,P為海表大氣壓;VP(H2O,s/w)為水汽平衡器內(nèi)飽和水汽壓,由水汽平衡器內(nèi)溫度和鹽度計(jì)算獲得。并且,通過((αpCO2/αT)/pCO2=0.042 3 ℃-1)[18]對(duì)水汽平衡器內(nèi)海水表層pCO2校正到現(xiàn)場(chǎng)實(shí)測(cè)溫度。

    表1 調(diào)查數(shù)據(jù)分布情況Table 1 Observational data

    1.2 海-氣CO2通量計(jì)算

    目前海-氣CO2通量是表征單位面積、時(shí)間的海-氣CO2轉(zhuǎn)移速率,一般通過如下塊體公式[20]進(jìn)行計(jì)算。

    FCO2=K×L×ΔpCO2,

    (2)

    式中,F(xiàn)CO2為海-氣CO2通量(mmol·m-2·d-1);K為水-氣氣體交換輸送率(mmol·m-2·d-1·μatm-1);L為CO2在與溫度、鹽度有關(guān)的溶解度(mmol·L-1·μatm-1);ΔpCO2為海-氣CO2偏壓差,常用uatm為單位,由ΔpCO2=pCO2water-pCO2air計(jì)算獲得。由于水-氣交換速度較難獲得,目前使用的K值大多采用Wanninkhof 1992年提出的水-氣氣體交換輸送率計(jì)算模式[21]:

    (3)

    式中,U10是離海面10 m高處的風(fēng)速(m/s);Sc是一定溫度下CO2的Schmidt常數(shù);660是指在20 ℃下的Schmidt常數(shù)。

    1.3 補(bǔ)充缺失數(shù)據(jù)

    由于觀測(cè)數(shù)據(jù)中部分航次風(fēng)速數(shù)據(jù)缺失,本文使用QuikSCAT衛(wèi)星風(fēng)速計(jì)在相應(yīng)位置、時(shí)刻的風(fēng)速數(shù)據(jù)對(duì)缺失數(shù)據(jù)進(jìn)行補(bǔ)充。為驗(yàn)證QuikSCAT風(fēng)速計(jì)數(shù)據(jù)精度,以船載氣象站獲得2008-04的航測(cè)風(fēng)速數(shù)據(jù)(時(shí)間分辨率為1 min)同QuikSCAT風(fēng)速數(shù)據(jù)進(jìn)行對(duì)比,結(jié)果見圖2。同實(shí)測(cè)數(shù)據(jù)相比,QuikSCAT日分辨率數(shù)據(jù)雖然在時(shí)間上存在一定的延遲,但較好地保留了風(fēng)速大小的變化趨勢(shì)。而月平均數(shù)據(jù)則完全過濾掉中、高風(fēng)速數(shù)據(jù),對(duì)計(jì)算水-氣交換速度結(jié)果的精確性具有重要影響[22],因此,選用日分辨率數(shù)據(jù)對(duì)缺失數(shù)據(jù)進(jìn)行補(bǔ)充。并且使用最鄰近差值法對(duì)大氣pCO2air進(jìn)行差值,計(jì)算水-氣氣體交換輸送率以及海-氣CO2通量。在此基礎(chǔ)上探討東海北部及黃海南部海域海-氣CO2通量的季節(jié)性變化及其控制因素。

    圖2 2008-04風(fēng)速對(duì)比分析Fig.2 Comparison of wind speed in April, 2008

    2 結(jié)果分析

    2.1 海-氣CO2通量的區(qū)域分布特征

    從碳循環(huán)角度來說,東海及黃海海域的CO2體系受控于一系列物理和生物地球化學(xué)過程,CO2源/匯分布呈現(xiàn)顯著的時(shí)空差異性[23]。根據(jù)2003—2008年間11個(gè)航次海-氣CO2通量結(jié)果平面分布如圖3所示,從圖中可以看出,整個(gè)研究區(qū)域內(nèi)海-氣CO2通量分布存在著較大的不均勻性,CO2通量變化較大(-29.41~79.15 mmol·m-2·d-1),這也從側(cè)面反映了邊緣海區(qū)的復(fù)雜性。冬季(1月份),觀測(cè)數(shù)據(jù)集中在Ⅱ區(qū)區(qū)域內(nèi),在長江口南支出現(xiàn)海-氣CO2通量最高值(79.15 mmol·m-2·d-1),表現(xiàn)為大氣CO2較強(qiáng)源,并且隨離岸距離增加而逐漸降低,從122°12′E處開始由大氣CO2源轉(zhuǎn)變?yōu)榇髿釩O2匯(圖4)。但整個(gè)區(qū)域均值為0.22 mmol·m-2·d-1,為大氣CO2弱源。這與李巖等[24]對(duì)河口海域的咸、淡水混合形成的渾濁帶的變化相似,并且相較于李巖等的研究,航線長度更長、變化也更為明顯。

    春季(3—5月),偏北風(fēng)逐漸減弱,溫度逐漸上升,Ⅱ區(qū)區(qū)域內(nèi)鹽度變化較大,整個(gè)海域海-氣CO2通量分布以長江口為界呈現(xiàn)南高北低的分布特征。同一季節(jié)不同區(qū)域海-氣CO2通量差異性較大。3-5月黃海南部海-氣CO2通量變化最為明顯,2008-03和2008-04在(121°30′E,34°30′N)區(qū)域分別出現(xiàn)通量最小、最大值(-23.15,2.41 mmol·m-2·d-1),該區(qū)域的大氣CO2匯逐漸減弱。并且東海南部由北向南大氣CO2匯強(qiáng)度也呈現(xiàn)逐漸減弱的趨勢(shì)。到5月份,東海南部到26°48′N以南由大氣CO2匯區(qū)轉(zhuǎn)變?yōu)榇髿釩O2源區(qū)(圖4)。東海中部外海區(qū)則表現(xiàn)為大氣CO2匯,并且在30°00′N附近出現(xiàn)最強(qiáng)匯。

    夏季(7—8月),整個(gè)海域海-氣CO2通量分布同春季相似,只是黃海南部CO2源進(jìn)一步加強(qiáng),并且在33°30′N附近呈現(xiàn)為大氣CO2源。在東海海域,除近岸河口區(qū)為源區(qū)外,其他大部分仍為大氣CO2匯。并且近岸河口區(qū)隨離岸距離越來越大,海-氣CO2源逐漸減弱,在122°36′E以東由CO2源轉(zhuǎn)化為CO2匯區(qū)(圖4),說明河口區(qū)CO2源面積隨淡水輸入量的增加逐漸擴(kuò)大。

    進(jìn)入秋季(09—11月)后,整個(gè)海域海-氣CO2源區(qū)面積達(dá)全年最大,南部CO2通量高于北部。并且,不同區(qū)域內(nèi)CO2源/匯變化較大。從緯向平均變化曲線可以看出,除2007-11外,其他月份僅在中部約28°30′~31°30′N范圍內(nèi)表現(xiàn)為大氣CO2匯,其他大部分區(qū)域?yàn)镃O2弱源(小于10 mmol·m-2·d-1),并且在2006年表現(xiàn)最為明顯。東西方向上變化較小,整體上表現(xiàn)為大氣CO2源。

    圖3 各航次海-氣CO2通量分布Fig.3 Distribution of air-sea CO2 flux

    圖4 CO2通量經(jīng)向、緯向平均結(jié)果Fig.4 Meridional and zonal mean of CO2 flux

    2.2 海-氣CO2通量變化控制因素分析

    從東海及黃海南部海域海-氣CO2通量分布上看,該海域海-氣CO2源、匯分布及其控制因素呈現(xiàn)明顯的季節(jié)變化,并且同溫度、季風(fēng)、生物、淡水輸入等因素的季節(jié)差異密切相關(guān)[25]。其中,溫度是從熱力學(xué)角度影響海水中CO2分布的一個(gè)重要因素,它與鹽度、壓力是表征海水物理過程、化學(xué)過程的基本參數(shù)[26]。并且,溫度的分布與變化,除取決于海區(qū)的熱量平衡外,還與地理環(huán)境、海流強(qiáng)弱、氣象條件等因素變化有關(guān)。而海表鹽度則能夠反映淡水輸入對(duì)陸架邊緣海的影響。因此,可以通過分析實(shí)測(cè)溫度(圖5)、鹽度同海-氣CO2通量間關(guān)系,再結(jié)合海水表層pCO2分布特征(圖6),分析該區(qū)域CO2源/匯格局的控制因素。

    圖5 各航次海水表層溫度分布Fig.5 Distribution of sea surface temperature

    圖6 各航次海水表層pCO2分布Fig.6 Distribution of sea surface pCO2

    長江口南支冬季海-氣CO2通量變化趨勢(shì)同觀測(cè)到的海水表層pCO2相同,并且同溫度、鹽度變化呈負(fù)相關(guān),相關(guān)系數(shù)分別為-0.55,-0.86,而李巖[24]得出的pCO2同鹽度的相關(guān)系數(shù)為-0.92,說明河口區(qū)域受陸源無機(jī)物輸入及咸、淡水混合影響比較顯著。在通量最大值區(qū)域出現(xiàn)溫、鹽最小值(6.23 ℃,0.47)。這是由于河口區(qū)冬季受富含無機(jī)物的低溫、低鹽水輸入以及海水混合的影響。淡水的輸入帶來了大量的無機(jī)物,導(dǎo)致海水透明度降低,光合作用減弱,CO2吸收減弱。而淺水區(qū)強(qiáng)烈的海水混合使得富含CO2的底層水上涌,加上冬季低溫抑制浮游植物對(duì)碳的同化,導(dǎo)致海水表層pCO2上升,成為大氣CO2源。隨著離岸距離越遠(yuǎn),海水溫度、透明度逐漸升高,光合作用逐漸增強(qiáng),混合減弱,使表層海水中CO2含量逐漸降低,致使河口區(qū)由大氣CO2源逐漸向大氣CO2匯過渡。這些都充分說明,在河口區(qū)冬季源/匯格局主要受淡水輸入和海水混合作用影響,以往研究也表明該區(qū)域內(nèi)受潮致混合影響比較明顯[27]。

    隨著春季向夏季的過渡,整個(gè)區(qū)域內(nèi)溫度逐漸上升,南北差異減小。而鹽度受中部海區(qū)長江沖淡水逐漸加強(qiáng)的影響逐漸降低。從局部區(qū)域看,黃海南部(Ⅰ區(qū))3-5月海表溫度平均上升了約6 ℃,而鹽度變化則較小。并且海水表層pCO2同海表溫度間呈現(xiàn)正相關(guān)(0.52 mmol·m-2·d-1)關(guān)系。說明山東半島沿岸南下低溫水水勢(shì)逐漸減弱,導(dǎo)致海水表層增溫明顯,光合作用越來越顯著,從而促進(jìn)生物對(duì)CO2的吸收,影響該區(qū)域源/匯分布格局[28]。最為明顯的是在黃海南部(121°30′E,34°30′N)附近3月份由大氣CO2較強(qiáng)匯,到5月份轉(zhuǎn)換為大氣CO2源區(qū)。根據(jù)目前研究,該區(qū)域除受沿岸流及生物作用的影響外還與底層富含CO2海水涌升作用加強(qiáng)有關(guān)[29-30]。東海近岸及中部陸架區(qū)(Ⅱ區(qū))溫、鹽變化較大,明顯受淡水輸入量增加的影響,海水穩(wěn)定性增加,混合作用減弱,光合作用占主導(dǎo)作用,因而整個(gè)區(qū)域表現(xiàn)CO2匯區(qū),并呈增強(qiáng)趨勢(shì)。東海南部海域(Ⅲ區(qū))由于受黑潮、臺(tái)灣暖流影響溫度上升,適宜的海水溫度再加上沿岸流營養(yǎng)物的輸入促進(jìn)了光合作用,海表pCO2降低,為CO2匯區(qū)。但舟山群島到臺(tái)州列島海域出現(xiàn)CO2強(qiáng)匯區(qū)(-56.24 mmol·m-2·d-1),與朱連磊等觀測(cè)結(jié)果差異較大[31],這可能與沿岸淡水輸入的年際波動(dòng)及海水混合有關(guān),但具體原因需進(jìn)一步研究。

    夏季,整個(gè)海區(qū)溫度、鹽度分布同春季相似,只是溫度升至全年最高,南北溫度分布較均勻,地區(qū)差異性小。黃海南部夏季由于上升流作用,底層富含CO2海水上涌,導(dǎo)致海表pCO2上升,出現(xiàn)pCO2最高值(813.98 μatm),表明該區(qū)域物流作用明顯強(qiáng)于生物作用。因而夏季黃海南部表現(xiàn)為CO2弱源(2.99 mmol·m-2·d-1),這與江春波等[29]提出的南黃海主要受底層水涌升以及長江沖淡水影響而導(dǎo)致夏季呈現(xiàn)為CO2弱源的結(jié)果較為一致的,但他的研究位置更靠近涌升區(qū),從而他觀測(cè)到的結(jié)果也高于本次計(jì)算的結(jié)果。東海中部及南部陸架區(qū)受沿岸流和長江沖淡水的影響,大量懸浮物及營養(yǎng)物的輸入促進(jìn)了浮游植物的生長、繁殖,吸收大量的碳。并且,密度較低的淡水浮于海水上層,使水體層化加劇,阻礙水體混合。底部生物呼吸產(chǎn)生的CO2無法轉(zhuǎn)移到海表,使得表層pCO2降低,呈現(xiàn)為CO2匯。通過對(duì)比2003-08與2006-01在東海近岸區(qū)CO2通量分布情況可以看出,徑流量大的月份CO2源區(qū)面積越大,說明徑流量的大小對(duì)近岸河口區(qū)CO2源/匯格局轉(zhuǎn)換具有重要影響。但局部如舟山群島(122°00′E,30°30′N)附近由于夏季西南風(fēng)的影響產(chǎn)生上升流導(dǎo)致該區(qū)域呈現(xiàn)為CO2源(38.24 mmol·m-2·d-1)[32]。

    受偏北季風(fēng)加強(qiáng),冷空氣南下的影響,導(dǎo)致表層水溫逐漸下降,南北溫差增大,整個(gè)區(qū)域秋季總體上表現(xiàn)為大氣CO2源區(qū)。在黃海南部黃海沿岸流和黃海暖流,以及東海海域臺(tái)灣暖流北上與長江沖淡水交匯,水體發(fā)生較強(qiáng)混合。溫度較高的黃海暖流以及臺(tái)灣暖流使區(qū)域內(nèi)CO2溶解度降低,使黃海南部以及東海南部總體上秋季呈現(xiàn)為碳源。不僅如此,根據(jù)2006-09(P1,P2,P3斷面)和2006-10(PN,S1斷面)五個(gè)斷面34個(gè)站位觀測(cè)結(jié)果計(jì)算該區(qū)域秋季混合層深度(MLD)、障礙層厚度等水文信息發(fā)現(xiàn),由于秋季季風(fēng)加強(qiáng),障礙層厚度較薄,平均為1~2 m,因此在季風(fēng)加強(qiáng)的情況下較容易使海水混合達(dá)到底部。并且從混合層深度與海表pCO2以及海氣CO2通量關(guān)系(圖7a)可以看出。秋季混合層深度達(dá)底部,致使夏季形成躍層被打破,海水底部由春季、夏季生物呼吸作用而產(chǎn)生的富含CO2海水上涌,海水表層pCO2升高,圖7a也證明這一觀點(diǎn),即混合層達(dá)底部的海表pCO2普遍較高。不同區(qū)域內(nèi),黃海南部9月混合層深度較淺(平均約15 m),為大氣CO2匯,CO2通量隨著混合層深度的加深而逐漸升高。從S1,P2,PN斷面可以看出(圖7b),中部受淡水輸入影響,海水穩(wěn)定性較強(qiáng),混合較弱,2006-10為CO2匯,并且隨混合層深度加深CO2匯逐漸增強(qiáng),這可能與該區(qū)域浮游植物、海水化學(xué)性質(zhì)以及該區(qū)域不同來源水團(tuán)CO2含量不同等有關(guān),這也需要對(duì)底層海水的來源及變化做進(jìn)一步研究[33]。東海南部9月份受混合層深度最為明顯,且由于混合層深度較深,表現(xiàn)為CO2源。由此可以認(rèn)為,水體垂直混合是引起東海及黃海南部海域秋季源/匯變化的一個(gè)重要因素。

    圖7 混合層深度(MLD)與海表pCO2、FCO2關(guān)系Fig.7 Relationship between mixed layer depth and pCO2 and FCO2

    為進(jìn)一步分析東海及黃海南部海域海-氣CO2通量季節(jié)性變化,對(duì)3個(gè)區(qū)域不同季節(jié)海表pCO2及CO2通量進(jìn)行統(tǒng)計(jì)(表2)。整個(gè)海域的3個(gè)局部區(qū)域,除秋季和黃海南部夏季為大氣CO2源外,其他區(qū)域在春季和夏季均表現(xiàn)為大氣CO2匯。黃海南部主要受長江沖淡水(夏季最強(qiáng))、黃海暖流(冬季最強(qiáng))、黃海沿岸流(冬季)以及底層水涌升的影響,海水表層pCO2波動(dòng)較大,總體上,春季為大氣CO2匯,夏季、冬季為大氣CO2源。東海近岸區(qū)則由于大量富含無機(jī)物淡水的輸入,再加上近岸海水混合阻礙了光合作用,導(dǎo)致該區(qū)域呈現(xiàn)為大氣CO2源,并且源區(qū)面積同淡水輸入量呈正比。另外,由于生物作用及溫躍層的影響,東海部陸架區(qū),冬季、春季、夏季表現(xiàn)為大氣CO2匯。秋季則由于海水混合導(dǎo)致底層富集的CO2上涌,造成表層pCO2升高,表現(xiàn)為CO2源。東海南部則受臺(tái)灣暖流和沿岸流季節(jié)性變化影響,該海域CO2通量分布也呈現(xiàn)出季節(jié)性變化。

    從不同月份海-氣CO2通量變化上可以看出,夏季海-氣CO2通量變化比較平穩(wěn),而春季和秋季CO2通量波動(dòng)較大。從春季到夏季,CO2匯逐漸增強(qiáng);從夏季到秋季由于海水混合的影響,整個(gè)海域逐漸由大氣CO2匯轉(zhuǎn)化為CO2源,并且在9月份達(dá)到最高,這可能是由于春季、夏季CO2在底層積累量的變化有關(guān)??梢灶A(yù)見,冬季該海域大部分區(qū)域?qū)⒊尸F(xiàn)為弱源,這也與冬季海水混合加強(qiáng)、生物作用減弱、溫度降低表層pCO2溶解度增大等因素有關(guān)。這些均與翟惟東[9-10]以及張龍軍等[14]的研究結(jié)果相近,但是從年際變化上看,同一區(qū)域不同年份,春季CO2匯呈減弱趨勢(shì),秋季CO2匯呈逐漸減弱甚至變?yōu)樵吹内厔?shì)。這可能是陸源污染物使水體富營養(yǎng)化、海表升溫等因素導(dǎo)致。

    表2 東海及黃海南部海域CO2通量Table 2 CO2 flux in the East China Sea and the South Yellow Sea

    需要引起重視的是,陸架邊緣海復(fù)雜的物理、生物和化學(xué)過程,使得不同年份觀測(cè)的水文、氣象要素間存在著較大差異,這也必然導(dǎo)致不同年份的海-氣CO2通量有一定的波動(dòng)。因此需要我們得到覆蓋范圍更大、精確度更高、時(shí)間尺度更長的CO2參數(shù)來深入研究控制碳循環(huán)的關(guān)鍵生物地球化學(xué)過程[34]。

    3 結(jié) 論

    依據(jù)東海及黃海南部海域海水表層pCO2走航觀測(cè)結(jié)果,并結(jié)合同步獲取的溫度、鹽度等參數(shù),對(duì)該海域海-氣CO2通量影響因素進(jìn)行分析。研究結(jié)果表明,近岸河口區(qū)受陸地徑流量季節(jié)性變化和咸淡水混合的影響,總體上表現(xiàn)為大氣CO2源,并且隨徑流量增大,CO2源面積也逐漸擴(kuò)大。春季黃海南部、東海中部及南部海域則由于海水混合較弱,適宜的海水溫度再加上沿岸流營養(yǎng)物的輸入促使光合作用增強(qiáng),因而整個(gè)海域表現(xiàn)為大氣CO2匯,分別為(-7.77±6.59),(-11.08±8.99),(-2.94±6.78) mmol·m-2·d-1。夏季除黃海南部底層富含CO2海水上涌為大氣CO2源,其他區(qū)域均為大氣CO2匯,這主要是由于淡水的輸入使海水穩(wěn)定性增強(qiáng),生物活動(dòng)產(chǎn)生的CO2因躍層的作用在海底富集,海水表層pCO2較低,大部分海區(qū)為大氣CO2匯。并且通過走航和站位數(shù)據(jù)的分析,秋季由于海表溫度降低、季風(fēng)加強(qiáng)等因素導(dǎo)致該區(qū)域總體上呈現(xiàn)為大氣CO2源,平均通量為(1.31±5.87) mmol·m-2·d-1。而水體垂直混合是引起長江口及其鄰近海域秋季源/匯變化的一個(gè)重要因素。

    由于陸架邊緣海物理、化學(xué)、生物過程較為復(fù)雜,致使陸架邊緣碳通量及其控制因素出現(xiàn)不確定性問題。而用于碳收支估算的直接觀測(cè)數(shù)據(jù)往往非常有限。僅僅通過時(shí)空覆蓋度都很局限的斷面調(diào)查來研究近海海域CO2海-氣通量的季節(jié)分布,并進(jìn)而估計(jì)其年度凈通量,可能存在偏歧。因此需要在對(duì)區(qū)域源/匯格局及其影響因素初步認(rèn)識(shí)的基礎(chǔ)上,進(jìn)一步研究物理、化學(xué)、生物過程對(duì)海-氣CO2通量的影響。把握該海區(qū)CO2海-氣通量調(diào)控機(jī)制,為未來實(shí)現(xiàn)通過數(shù)值模擬以及遙感等技術(shù)進(jìn)行大面、精確、實(shí)時(shí)CO2觀測(cè)、模擬提供可能,也為今后研究陸架邊緣海在全球氣候變化中的作用提供研究基礎(chǔ)。

    致謝:廈門大學(xué)及國家海洋監(jiān)測(cè)中心提供CO2相關(guān)數(shù)據(jù),美國遙感系統(tǒng)(RSS)提供部分風(fēng)速數(shù)據(jù)下載,國家海洋局第二研究所許建平研究員對(duì)本文提出了寶貴意見!

    [1] GUO X H. A comparative study on the carbonate system between the Pearl and Mississippi River estuaries[D].Xiamen:Xiamen University, 2009. 郭香會(huì). 珠江與密西西比河口碳酸鹽系統(tǒng)的比較研究[D].廈門: 廈門大學(xué), 2009.

    [2] FENNEL K. The role of continental shelves in nitrogen and carbon cycling: Northwestern North Atlantic case study[J].Ocean Science Discussions, 2010, 7(1):177-205.

    [3] LIU Z Y. The changing behavior carbon form estuarine dynamics in Yellow River Estuary[D].Qingdao:Ocean University of China, 2011. 劉志媛. 黃河河口過程中碳的行為變化[D].青島:中國海洋大學(xué), 2011.

    [4] CAI W J, DAI M, WANG Y. Air-sea exchange of carbon dioxide in ocean margins: A province-based synthesis[J].Geophysical Research Letters, 2006, 33(12):347-366.

    [5] BORGES A V, DELILLE B, FRANKIGNOULLE M. Budgeting sinks and sources of CO2in the coastal ocean: Diversity of ecosystems counts[J].Geophysical Research Letters, 2005, 32(14):301-320.

    [6] TSUNOGAI S, WATANABE S, SATO T. Is there a "continental shelf pump" for the absorption of atmospheric CO2?[J].Tellus(Series B), 1999, 51(3): 701-712.

    [7] SHIM J H, KIM D, KANG Y C, et al. Seasonal variations inpCO2and its controlling factors in surface seawater of the northern East China Sea[J].Continental Shelf Research, 2007, 27(20): 2623-2636.

    [8] OLIVEIRA A P, MATEUS M D, CABECADAS G, et al. Water-air CO2fluxes in the Tagus estuary plume(Portugal) during two distinct winter episodes[J].Carbon Balance & Management, 2015, 10(1):1-15.

    [9] ZHAI W D, DAI M, CAI W J. Coupling of surfacepCO2and dissolved oxygen in the northern South China Sea: impacts of contrasting coastal processes[J].Biogeosciences, 2009, 6(11): 2589-2598.

    [10] ZHAI W D, DAI M, CAI W J, et al. The partial pressure of carbon dioxide and air sea fluxes in the northern South China Sea in spring, summer and autumn[J].Marine Chemistry, 2005, 96(1): 87-97.

    [11] ZHANG L, XUE L, SONG M, et al. Distribution of the surface partial pressure of CO2in the southern Yellow Sea and its controls[J].Continental Shelf Research, 2010, 30(3): 293-304.

    [12] TAN Y, ZHANG L J, WANG F, et al. Summer surface waterpCO2and CO2flux at air-sea interface in western part of the east China sea[J].Oceanologia et Limnologia Sinica, 2004, 35(3): 239-245. 譚燕, 張龍軍, 王凡, 等. 夏季東海西部表層海水中的pCO2及海-氣界面通量[J].海洋與湖沼, 2004, 35(3): 239-245.

    [13] QU B X, SONG J M, YUAN H M, et al. Advances of seasonal variations and controlling factors of the sea-air CO2flux in the East China Sea[J].Advances in Earth Science, 2013, 28(7): 783-793. 曲寶曉, 宋金明, 袁華茂, 等. 東海海—?dú)饨缑娑趸纪康募竟?jié)變化與控制因素研究進(jìn)展[J].地球科學(xué)進(jìn)展, 2013, 28(7): 783-793.

    [14] ZHANG L J, WANG B Y, ZHANG J.pCO2in the Surface Water of the East China Sea in Winter and Summer[J].Periodical of Ocean University of Qingdao, 1999(Suppl.1): 149-153. 張龍軍, 王彬宇, 張經(jīng). 東海冬、夏兩季表層海水的二氧化碳分壓[J].青島海洋大學(xué)學(xué)報(bào)(自然科學(xué)版), 1999(增1): 149-153.

    [15] ZHAI W D, DAI M. On the seasonal variation of air-sea CO2fluxes in the outer Changjiang(Yangtze River) Estuary, East China Sea[J].Marine Chemistry, 2009, 117(1): 2-10.

    [16] ZHAI W D, CHEN J F, JIN H Y, et al. Spring carbonate chemistry dynamics of surface waters in the northern East China Sea: Water mixing, biological uptake of CO2, and chemical buffering capacity[J].Journal of Geophysical Research: Oceans, 2014, 119(9): 5638-5653.

    [17] DICKSON A G, GOYET C. Handbook of methods for the analysis of the various parameters of the carbon dioxide system in sea water. Version 2[R].United States:Oak Ridge National Lab., TN, 1994.

    [18] TAKAHASHI T, OLAFSSON J, GODDARD J G, et al. Seasonal variation of CO2and nutrients in the high-latitude surface oceans: A comparative study[J].Global Biogeochemical Cycles, 1993, 7(4): 843-878.

    [19] CHEN C T A, ZHAI W, DAI M. Riverine input and air-sea CO2exchanges near the Changjiang (Yangtze River) Estuary: status quo and implication on possible future changes in metabolic status[J].Continental Shelf Research, 2008, 28(12): 1476-1482.

    [20] FRANKIGNOULLE M, BOURGE I, WOLLAST R. Atmospheric CO2fluxes in a highly polluted estuary (the Scheldt)[J].Limnology and Oceanography, 1996, 41(2): 365-369.

    [21] WANNINKHOF R. Relationship between wind speed and gas exchange over the ocean[J].Journal of Geophysical Research: Oceans (1978-2012), 1992, 97(C5): 7373-7382.

    [22] ZHAI W D, DAI M H, CHEN B S, et al. Seasonal variations of sea air CO2fluxes in the largest tropical marginal sea(South China Sea) based on multiple-year underway measurements[J].Biogeosciences, 2013, 10(11): 7775-7791.

    [23] LI N. Dissolved organic and inorganic carbon nearby Changjiang Estuary[D].Qingdao: Ocean University of China, 2011. 李寧. 長江口鄰近海域的溶解有機(jī)碳和溶解無機(jī)碳系統(tǒng)[D].青島:中國海洋大學(xué), 2011.

    [24] LI Y, ZHANG L J, SU Z, et al.pCO2and its influencing factors at the initial mixing stage of freshwater and saltwater in the Changjiang River estuary[J].Periodical of Ocean University of China, 2006, 36(2): 295-298. 李巖, 張龍軍, 蘇征,等. 長江口淡水端淡、鹽水混合表層pCO2的急劇變化及其影響機(jī)制[J].中國海洋大學(xué)學(xué)報(bào)(自然科學(xué)版), 2006, 36(2): 295-298.

    [25] BAI Y, CAI W J, HE X, et al. A mechanistic semi-analytical method for remotely sensing sea surfacepCO2in river-dominated coastal oceans: A case study from the East China Sea[J].Journal of Geophysical Research: Oceans, 2015, 120(3): 2331-2349.

    [26] MILLERO F J. The marine inorganic carbon cycle[J].Chemical Reviews, 2007, 107(2): 308-341.

    [27] LU X G, QIAO F L, XIA C S, et al. The Yangtze River Estuary and Zhejiang coastal upwelling in summer tide mechanism[J].Science in China, 2007, 37(1): 133-144. 呂新剛, 喬方利, 夏長水,等. 長江口外及浙江沿岸夏季上升流的潮生機(jī)制[J].中國科學(xué), 2007,37(1): 133-144.

    [28] XUE L, XUE M, ZHANG L, et al. Surface partial pressure of CO2and air-sea exchange in the northern Yellow Sea[J].Journal of Marine Systems, 2012, 105-108(12): 194-206.

    [29] JIANG C B, ZHANG L J, WANG F. A study ofpCO2in the surface water of the southern Yellow Sea in summer II: The respective contribution of the upwelling and the Yangtze River diluted water to the air-sea of CO2flux[J].Periodical of Ocean University of China, 2006, 36(Suppl.1): 147-152. 江春波, 張龍軍, 王峰. 南黃海夏季海水pCO2研究Ⅱ——下層海水涌升和長江沖淡水對(duì)海-氣界面 CO2通量的貢獻(xiàn)[J].中國海洋大學(xué)學(xué)報(bào)(自然科學(xué)版), 2006, 36(增1): 147-152.

    [30] SONG M Q, ZHANG L J, JIANG C B. The distribution of partial pressure of CO2and its maintaining mechanism: vertical mixing versus biological activity in the Yellow Sea in early spring[J].Periodical of Ocean University of China, 2007 (Suppl.1): 67-72. 宋美芹, 張龍軍, 江春波. 初春(3月)黃海水體垂直混合和生物活動(dòng)對(duì)水-氣界面pCO2分布的控制作用[J].中國海洋大學(xué)學(xué)報(bào)(自然科學(xué)版), 2007(增1): 67-72.

    [31] ZHU L L. Dissolved inorganic carbon and air-sea carbon fluxes in the min-west of the East China Sea[D].Qingdao: Graduate University of Chinese Academy of Sciences, 2011. 朱連磊. 東海中西部海域水體無機(jī)碳與海氣界面碳通量[D].青島: 中國科學(xué)院研究生院, 2011.

    [32] PAN Y T, SHA W Y. Numerical study on the summer coastal upwelling off Fujian and Zhejiang[J].Marine Science Bulletin, 2004, 23(3): 1-11. 潘玉萍, 沙文鈺. 夏季閩浙沿岸上升流的數(shù)值研究[J].海洋通報(bào), 2004, 23(3): 1-11.

    [33] WANNINKHOF R, BARBERO L, BYME R, et al. Ocean acidification along the Gulf Coast and east coast of the USA[J].Continental Shelf Research, 2015, 88:54-71.

    [34] SONG J M. Carbon sources and sinks in oceans[J].Marine Environmental Science, 2003, 22(2): 75-80. 宋金明. 海洋碳的源與匯[J].海洋環(huán)境科學(xué), 2003, 22(2): 75-80.

    Received: December 7, 2015

    Seasonal Variability of Shipboard Air-Sea CO2Flux in the East China Sea and South Yellow Sea

    HU Deng-hui1,2,GAO Guo-ping1,2,ZHAI Wei-dong3,ZHANG Chun-ling1,2

    (1.CollegeofMarineSciences,ShanghaiOceanUniversity, Shanghai 201306, China; 2.KeyLaboratoryofSustainableExploitationofOceanicFisheriesResources,MinistryofEducation, Shanghai 201306, China; 3.NationalMarineEnvironmentalMonitoringCenter, Dalian 116023, China)

    Seasonal variability of sea surfacepCO2and air-sea CO2flux and relationship among air-sea CO2transfer, sea surface temperature and salinity are analyzed using field observation data from 2003 to 2008 in Yangtze estuary and adjacent sea area. Results show that the sea surfacepCO2and air-sea CO2have significant seasonal variations. The South Yellow Sea, the offshore area and middle continental shelf of the East China Sea, the South East China Sea all demonstrate atmospheric CO2convergence in spring due to the rising sea surface temperature and strengthening biological action, and their seasonal averaged air-sea CO2flux is(-7.77±6.59),(-11.08±8.99) and(-2.94±6.78)(mmol·m-2·d-1), respectively. While, in summer, the South Yellow Sea is a CO2source(2.99±6.09)(mmol·m-2·d-1) due to weaker upwelling, and the offshore area and middle continental shelf of the East China sea(-4.81±8.92)(mmol·m-2·d-1), the South East China Sea(-0.75±12.14)(mmol·m-2·d-1) are atmospheric CO2convergence zones, because thermocline is formed when fresh water input blocks water mixing, coupled with enhancement of photosynthesis. In autumn, due to enhanced northerly winds and increasing water mixing as well as upwelling of CO2-rich sea water which increases sea surfacepCO2, the entire sea area becomes an atmospheric CO2source. In summary, carbon convergence presents a weakening trend in spring but a gradual increasing trend in autumn.

    air-sea CO2flux;pCO2; seasonal variability; factors

    2015-12-07

    上海市科委重點(diǎn)支撐項(xiàng)目——典型海域海-氣CO2通量分析系統(tǒng)研制及應(yīng)用示范(13230502000)

    胡登輝(1990-),男,湖北黃岡人,碩士研究生,主要從事海-氣CO2交換方面研究.E-mail:hudh_shou@163.com

    *通信作者:高郭平(1972-),男,江蘇泰州人,博士,教授,主要從事海洋觀測(cè)、海洋動(dòng)力學(xué)方面研究.E-mail:gpgao@shou.edu.cn

    P734

    A

    1671-6647(2016)04-0474-12

    10.3969/j.issn.1671-6647.2016.04.003

    猜你喜歡
    海表黃海東海
    東方濕地 黃海明珠
    基于無人機(jī)的海表環(huán)境智能監(jiān)測(cè)系統(tǒng)設(shè)計(jì)與應(yīng)用
    2016與1998年春季北大西洋海表溫度異常的差異及成因
    融合海表溫度產(chǎn)品在渤黃東海的對(duì)比分析及初步驗(yàn)證
    黃海簡介
    東海采風(fēng)作品選
    中華詩詞(2019年2期)2019-08-27 01:59:16
    太陽總輻照度對(duì)熱帶中太平洋海表溫度年代際變化的可能影響
    三角恒等變換去哪兒了
    精衛(wèi)鳴,歌東海
    李東海的詩
    詩潮(2014年7期)2014-02-28 14:11:30
    91久久精品电影网| 免费黄色在线免费观看| 国产有黄有色有爽视频| 久久综合国产亚洲精品| 国产探花在线观看一区二区| 久久精品久久久久久噜噜老黄| 能在线免费观看的黄片| 肉色欧美久久久久久久蜜桃 | 国产真实伦视频高清在线观看| 欧美精品一区二区大全| 亚洲,欧美,日韩| 久热久热在线精品观看| 色播亚洲综合网| 99久久中文字幕三级久久日本| av专区在线播放| 舔av片在线| 人人妻人人看人人澡| or卡值多少钱| 国产一区二区三区综合在线观看 | 国产精品久久久久久久电影| 精品久久久精品久久久| 亚洲欧美中文字幕日韩二区| 亚洲欧美一区二区三区国产| 国产精品蜜桃在线观看| 久久久a久久爽久久v久久| 日韩视频在线欧美| 欧美激情国产日韩精品一区| 国产高清不卡午夜福利| 日韩人妻高清精品专区| 亚洲精品亚洲一区二区| 久久99精品国语久久久| 婷婷色综合大香蕉| 亚洲av免费在线观看| 91久久精品国产一区二区三区| 男女边吃奶边做爰视频| 69人妻影院| 亚洲欧美中文字幕日韩二区| 国产亚洲一区二区精品| 99热这里只有是精品在线观看| 亚洲美女搞黄在线观看| 成人一区二区视频在线观看| 国产成人精品福利久久| 男女啪啪激烈高潮av片| 自拍偷自拍亚洲精品老妇| 亚洲一级一片aⅴ在线观看| 高清av免费在线| 久久综合国产亚洲精品| 亚洲天堂国产精品一区在线| 建设人人有责人人尽责人人享有的 | 国产片特级美女逼逼视频| 国产精品一二三区在线看| 日本一本二区三区精品| 高清视频免费观看一区二区 | 中文乱码字字幕精品一区二区三区 | 九草在线视频观看| av网站免费在线观看视频 | 日韩欧美精品免费久久| 国产精品嫩草影院av在线观看| 精品一区二区三区视频在线| 在线 av 中文字幕| 99久国产av精品| 一级爰片在线观看| 国产欧美日韩精品一区二区| 国产探花极品一区二区| 三级男女做爰猛烈吃奶摸视频| 国产黄色免费在线视频| www.色视频.com| 国产免费视频播放在线视频 | 亚州av有码| 午夜福利在线观看吧| 熟妇人妻不卡中文字幕| 国产老妇女一区| 日日撸夜夜添| 中文乱码字字幕精品一区二区三区 | 最新中文字幕久久久久| 免费观看无遮挡的男女| 中文在线观看免费www的网站| 国产高清不卡午夜福利| 九色成人免费人妻av| 男人和女人高潮做爰伦理| 亚洲在线自拍视频| 成人毛片a级毛片在线播放| 国产成人精品婷婷| 国产中年淑女户外野战色| a级一级毛片免费在线观看| 亚洲精品aⅴ在线观看| 精品午夜福利在线看| 亚洲四区av| 亚洲人成网站高清观看| 成人毛片60女人毛片免费| 亚洲自偷自拍三级| 少妇人妻精品综合一区二区| 欧美日本视频| 看非洲黑人一级黄片| 啦啦啦中文免费视频观看日本| 亚洲一区高清亚洲精品| 免费无遮挡裸体视频| 少妇被粗大猛烈的视频| 一区二区三区免费毛片| 国产高清三级在线| 日韩一本色道免费dvd| 国产黄a三级三级三级人| 日韩亚洲欧美综合| 亚洲在线自拍视频| 1000部很黄的大片| 国产精品女同一区二区软件| 麻豆av噜噜一区二区三区| 人妻一区二区av| av免费观看日本| 亚洲高清免费不卡视频| 18禁在线无遮挡免费观看视频| 熟女人妻精品中文字幕| 国产乱来视频区| 久久久成人免费电影| 亚洲欧美成人综合另类久久久| 精品人妻视频免费看| 亚洲精品色激情综合| 淫秽高清视频在线观看| 久久97久久精品| 卡戴珊不雅视频在线播放| 一区二区三区乱码不卡18| 男人爽女人下面视频在线观看| av黄色大香蕉| 午夜福利高清视频| 综合色av麻豆| 秋霞在线观看毛片| 亚洲成色77777| 高清日韩中文字幕在线| 水蜜桃什么品种好| 日韩成人伦理影院| 国产精品av视频在线免费观看| 免费观看无遮挡的男女| 九九爱精品视频在线观看| 亚洲经典国产精华液单| 国产 亚洲一区二区三区 | 91久久精品国产一区二区三区| 99热这里只有是精品在线观看| 中文字幕av在线有码专区| 亚洲av中文av极速乱| 91久久精品国产一区二区三区| 国产亚洲一区二区精品| 国产男女超爽视频在线观看| 国产一区亚洲一区在线观看| 亚洲欧美中文字幕日韩二区| 亚洲精品一二三| 久久99蜜桃精品久久| 国产精品女同一区二区软件| 九九爱精品视频在线观看| 亚洲国产高清在线一区二区三| 亚洲18禁久久av| 国产精品一区www在线观看| 久久久久久久久中文| 国产精品国产三级专区第一集| 听说在线观看完整版免费高清| 啦啦啦啦在线视频资源| 国产高清不卡午夜福利| av一本久久久久| 国产亚洲一区二区精品| 精品国产露脸久久av麻豆 | 国产精品伦人一区二区| 白带黄色成豆腐渣| 午夜激情久久久久久久| 国产单亲对白刺激| 女人久久www免费人成看片| 日韩欧美三级三区| 亚洲一区高清亚洲精品| 全区人妻精品视频| 久久精品夜色国产| 国产精品1区2区在线观看.| 精品午夜福利在线看| 久久人人爽人人片av| 人妻一区二区av| 草草在线视频免费看| 五月天丁香电影| 18+在线观看网站| 日本午夜av视频| 欧美xxxx黑人xx丫x性爽| 精品一区二区三卡| 人人妻人人澡欧美一区二区| videossex国产| 亚洲三级黄色毛片| 91av网一区二区| videossex国产| 成人欧美大片| 色综合亚洲欧美另类图片| 在线免费十八禁| 亚洲,欧美,日韩| 国产成人免费观看mmmm| 国产视频内射| 亚洲国产高清在线一区二区三| a级毛片免费高清观看在线播放| 超碰97精品在线观看| 亚洲成色77777| 麻豆国产97在线/欧美| 亚洲成人一二三区av| 欧美区成人在线视频| 床上黄色一级片| 久久99热6这里只有精品| 国产精品一区二区三区四区免费观看| 亚洲内射少妇av| 久久精品久久久久久久性| h日本视频在线播放| 777米奇影视久久| 国产精品国产三级国产专区5o| 国产黄片美女视频| 99久久精品热视频| 男人和女人高潮做爰伦理| 国产激情偷乱视频一区二区| 欧美xxⅹ黑人| 亚洲av电影在线观看一区二区三区 | 免费看日本二区| av福利片在线观看| 麻豆国产97在线/欧美| 欧美丝袜亚洲另类| 久久久久久久久久人人人人人人| 久久精品综合一区二区三区| 免费看日本二区| 免费观看在线日韩| 日韩av免费高清视频| 最近2019中文字幕mv第一页| 国产精品.久久久| 欧美精品国产亚洲| 狠狠精品人妻久久久久久综合| 不卡视频在线观看欧美| 麻豆久久精品国产亚洲av| 天堂网av新在线| av在线亚洲专区| 日韩精品青青久久久久久| 亚洲自偷自拍三级| 国产精品嫩草影院av在线观看| 老师上课跳d突然被开到最大视频| 少妇熟女aⅴ在线视频| 日本色播在线视频| 啦啦啦中文免费视频观看日本| 免费看日本二区| www.av在线官网国产| 3wmmmm亚洲av在线观看| 日韩av在线大香蕉| 直男gayav资源| av天堂中文字幕网| 亚洲一级一片aⅴ在线观看| 2021少妇久久久久久久久久久| 精品久久久久久成人av| 国产精品人妻久久久久久| 肉色欧美久久久久久久蜜桃 | 自拍偷自拍亚洲精品老妇| 成人漫画全彩无遮挡| 精品一区二区三区视频在线| 少妇人妻精品综合一区二区| 人人妻人人看人人澡| 日本黄大片高清| 国产精品av视频在线免费观看| 日日摸夜夜添夜夜爱| 18禁动态无遮挡网站| 九草在线视频观看| 日韩中字成人| 国产精品女同一区二区软件| 日韩 亚洲 欧美在线| 丝袜喷水一区| 成人亚洲欧美一区二区av| 亚洲怡红院男人天堂| 欧美zozozo另类| 蜜桃亚洲精品一区二区三区| 天天一区二区日本电影三级| 午夜福利高清视频| 777米奇影视久久| 婷婷色综合www| 国产色爽女视频免费观看| 夜夜看夜夜爽夜夜摸| 亚洲av免费在线观看| 日本欧美国产在线视频| 天堂网av新在线| .国产精品久久| 汤姆久久久久久久影院中文字幕 | 一边亲一边摸免费视频| 久久久久精品性色| 亚洲丝袜综合中文字幕| freevideosex欧美| 日韩亚洲欧美综合| 欧美最新免费一区二区三区| 一个人免费在线观看电影| 国产色婷婷99| 国产不卡一卡二| 亚洲精品第二区| 欧美区成人在线视频| 最近视频中文字幕2019在线8| 天堂影院成人在线观看| 又大又黄又爽视频免费| 51国产日韩欧美| 欧美区成人在线视频| 视频中文字幕在线观看| 日韩成人伦理影院| 国产 一区精品| 日韩亚洲欧美综合| 国产黄片美女视频| 午夜激情久久久久久久| 水蜜桃什么品种好| 又爽又黄无遮挡网站| 久久人人爽人人片av| av在线亚洲专区| 汤姆久久久久久久影院中文字幕 | 97人妻精品一区二区三区麻豆| 搡女人真爽免费视频火全软件| 91久久精品国产一区二区成人| 中文字幕人妻熟人妻熟丝袜美| 国产亚洲精品av在线| 久久精品国产自在天天线| 久久久久性生活片| 亚洲精品aⅴ在线观看| 大片免费播放器 马上看| 国产精品不卡视频一区二区| 色哟哟·www| 日本一本二区三区精品| 日韩大片免费观看网站| 99久久人妻综合| 亚洲成人精品中文字幕电影| 能在线免费观看的黄片| 欧美精品一区二区大全| 国产淫片久久久久久久久| av福利片在线观看| 欧美成人a在线观看| 国产伦在线观看视频一区| 欧美精品国产亚洲| 国产综合懂色| 色综合亚洲欧美另类图片| 久久久久久久大尺度免费视频| 男人舔奶头视频| 老司机影院成人| 成人亚洲精品一区在线观看 | 日韩欧美精品v在线| 国产黄色视频一区二区在线观看| 久久久久久久久久久免费av| 国产精品人妻久久久影院| 国产片特级美女逼逼视频| 国产成人一区二区在线| 菩萨蛮人人尽说江南好唐韦庄| av国产久精品久网站免费入址| 夜夜爽夜夜爽视频| 日本午夜av视频| 好男人在线观看高清免费视频| 国产伦一二天堂av在线观看| 高清日韩中文字幕在线| videossex国产| 秋霞伦理黄片| 久久精品夜夜夜夜夜久久蜜豆| 高清欧美精品videossex| 午夜免费男女啪啪视频观看| 在线观看一区二区三区| 日本欧美国产在线视频| 国产在线男女| 日韩伦理黄色片| 日本一二三区视频观看| 午夜视频国产福利| 哪个播放器可以免费观看大片| 成人无遮挡网站| 一级毛片aaaaaa免费看小| 国产精品国产三级专区第一集| 一二三四中文在线观看免费高清| 观看免费一级毛片| 中文字幕av成人在线电影| 最近视频中文字幕2019在线8| 欧美丝袜亚洲另类| 大香蕉久久网| 在线免费观看的www视频| a级毛色黄片| 天天躁夜夜躁狠狠久久av| 久久国内精品自在自线图片| 久久久久久久久久人人人人人人| 好男人视频免费观看在线| 亚洲最大成人中文| 91在线精品国自产拍蜜月| 中文资源天堂在线| 我要看日韩黄色一级片| 国产伦理片在线播放av一区| 乱人视频在线观看| 午夜免费观看性视频| 日韩三级伦理在线观看| 乱码一卡2卡4卡精品| 草草在线视频免费看| 亚洲精品国产av蜜桃| 亚洲成人精品中文字幕电影| 婷婷色av中文字幕| 深夜a级毛片| 麻豆乱淫一区二区| 国产精品国产三级专区第一集| 国产麻豆成人av免费视频| 中文天堂在线官网| 亚洲精品成人av观看孕妇| 天天躁夜夜躁狠狠久久av| 国产人妻一区二区三区在| 国产一级毛片七仙女欲春2| 亚洲va在线va天堂va国产| 在线a可以看的网站| av国产久精品久网站免费入址| 国产精品1区2区在线观看.| 卡戴珊不雅视频在线播放| 十八禁网站网址无遮挡 | 99久久精品一区二区三区| 亚洲熟妇中文字幕五十中出| 国产探花在线观看一区二区| 国产一区二区亚洲精品在线观看| 久久久欧美国产精品| 97在线视频观看| 美女被艹到高潮喷水动态| 2018国产大陆天天弄谢| 成人一区二区视频在线观看| 欧美另类一区| 免费av观看视频| 亚州av有码| 国产日韩欧美在线精品| 欧美精品国产亚洲| 欧美变态另类bdsm刘玥| 99九九线精品视频在线观看视频| 国产毛片a区久久久久| 91精品伊人久久大香线蕉| 亚洲色图av天堂| 91精品伊人久久大香线蕉| 99久久九九国产精品国产免费| 91精品伊人久久大香线蕉| 亚洲色图av天堂| 97热精品久久久久久| 中文字幕久久专区| 我的老师免费观看完整版| 国产熟女欧美一区二区| 一级a做视频免费观看| 日韩一区二区三区影片| 精品国产一区二区三区久久久樱花 | 免费播放大片免费观看视频在线观看| 久久99热6这里只有精品| 国产爱豆传媒在线观看| 日本一二三区视频观看| 日韩人妻高清精品专区| 99久国产av精品| 亚洲国产av新网站| 欧美97在线视频| 一级毛片电影观看| 免费看光身美女| 一区二区三区免费毛片| 亚洲欧美一区二区三区黑人 | 亚洲av二区三区四区| 国产视频首页在线观看| 亚洲人成网站在线播| 国产精品综合久久久久久久免费| 69人妻影院| 国产伦理片在线播放av一区| 黑人高潮一二区| 亚洲图色成人| 少妇被粗大猛烈的视频| 久久草成人影院| 深夜a级毛片| 干丝袜人妻中文字幕| 亚洲av国产av综合av卡| av女优亚洲男人天堂| 成年女人看的毛片在线观看| 一级毛片我不卡| 亚洲成色77777| 国精品久久久久久国模美| 伦精品一区二区三区| 淫秽高清视频在线观看| 午夜免费观看性视频| 床上黄色一级片| 亚洲人成网站高清观看| 亚洲va在线va天堂va国产| 国产伦在线观看视频一区| 日韩三级伦理在线观看| 欧美 日韩 精品 国产| 午夜福利高清视频| 少妇猛男粗大的猛烈进出视频 | 99久久九九国产精品国产免费| 可以在线观看毛片的网站| 看黄色毛片网站| 中文欧美无线码| 精品久久国产蜜桃| 在线 av 中文字幕| 简卡轻食公司| 国产永久视频网站| 亚洲精品日本国产第一区| 日韩av在线免费看完整版不卡| 免费在线观看成人毛片| 亚洲欧美中文字幕日韩二区| 少妇的逼好多水| 熟女人妻精品中文字幕| 国产成人精品一,二区| 1000部很黄的大片| 五月玫瑰六月丁香| 亚洲成人久久爱视频| 精品酒店卫生间| 国产高清有码在线观看视频| 免费无遮挡裸体视频| 一级毛片久久久久久久久女| 少妇的逼水好多| 免费看美女性在线毛片视频| 18+在线观看网站| 校园人妻丝袜中文字幕| 欧美人与善性xxx| 久久97久久精品| 十八禁网站网址无遮挡 | 亚洲av中文av极速乱| 日本欧美国产在线视频| 国产精品福利在线免费观看| 中文字幕免费在线视频6| 亚洲成人精品中文字幕电影| 国产精品蜜桃在线观看| 免费观看a级毛片全部| 精品人妻一区二区三区麻豆| 七月丁香在线播放| 赤兔流量卡办理| 亚洲内射少妇av| 高清视频免费观看一区二区 | 亚洲国产最新在线播放| 欧美区成人在线视频| 春色校园在线视频观看| 成人毛片60女人毛片免费| 精品欧美国产一区二区三| 国产亚洲5aaaaa淫片| 蜜桃久久精品国产亚洲av| www.色视频.com| 在线观看人妻少妇| 日韩强制内射视频| 中国美白少妇内射xxxbb| 久久韩国三级中文字幕| 纵有疾风起免费观看全集完整版 | 中文资源天堂在线| 亚洲精品亚洲一区二区| 久久精品久久精品一区二区三区| 国产成年人精品一区二区| 久久99热这里只频精品6学生| 成年免费大片在线观看| 麻豆成人av视频| 欧美日韩一区二区视频在线观看视频在线 | 亚洲美女搞黄在线观看| 九九在线视频观看精品| 九九久久精品国产亚洲av麻豆| 国产美女午夜福利| 草草在线视频免费看| 欧美3d第一页| av免费观看日本| 又黄又爽又刺激的免费视频.| or卡值多少钱| 伦精品一区二区三区| 春色校园在线视频观看| 国产午夜精品久久久久久一区二区三区| 国产亚洲一区二区精品| 高清日韩中文字幕在线| 亚洲人成网站高清观看| 久久精品久久精品一区二区三区| 国内精品美女久久久久久| 国产淫片久久久久久久久| 久久精品国产自在天天线| 狠狠精品人妻久久久久久综合| 欧美精品国产亚洲| 日本熟妇午夜| 在线观看免费高清a一片| 国产爱豆传媒在线观看| av女优亚洲男人天堂| 精华霜和精华液先用哪个| 精品国内亚洲2022精品成人| 国产成人精品婷婷| 少妇熟女aⅴ在线视频| 一区二区三区四区激情视频| 99久久九九国产精品国产免费| 久久久精品欧美日韩精品| 欧美日本视频| 91久久精品国产一区二区成人| 免费av不卡在线播放| 日本欧美国产在线视频| 久久精品国产自在天天线| 成人午夜精彩视频在线观看| 国产亚洲午夜精品一区二区久久 | 亚洲经典国产精华液单| 成人av在线播放网站| 国产精品久久久久久久电影| 日本黄色片子视频| 国产成年人精品一区二区| 午夜视频国产福利| 成人无遮挡网站| 汤姆久久久久久久影院中文字幕 | 日韩一区二区视频免费看| 久久久成人免费电影| 十八禁国产超污无遮挡网站| 色综合站精品国产| 男人狂女人下面高潮的视频| 久久久a久久爽久久v久久| 亚洲精华国产精华液的使用体验| 亚洲成人精品中文字幕电影| 日韩欧美精品免费久久| 99re6热这里在线精品视频| 国产欧美另类精品又又久久亚洲欧美| 国产一区有黄有色的免费视频 | 国产免费一级a男人的天堂| 欧美另类一区| 久久久a久久爽久久v久久| 国产黄片美女视频| 嫩草影院精品99| 日韩欧美一区视频在线观看 | 国产精品一区二区三区四区免费观看| 日韩一本色道免费dvd| 七月丁香在线播放| 人妻少妇偷人精品九色| 久久鲁丝午夜福利片| 亚洲怡红院男人天堂| 人妻系列 视频| 人人妻人人澡人人爽人人夜夜 | 91久久精品国产一区二区三区| 十八禁国产超污无遮挡网站| 国产 一区精品| 国产中年淑女户外野战色| 两个人的视频大全免费| 97热精品久久久久久| 网址你懂的国产日韩在线| 欧美性猛交╳xxx乱大交人| 欧美+日韩+精品| 国产精品女同一区二区软件| 久久这里有精品视频免费| 有码 亚洲区|