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

    陸架海碳循環(huán)模式研究現狀與進展*

    2014-10-08 12:49:32羅曉凡
    海洋科學進展 2014年2期
    關鍵詞:陸架碳循環(huán)營養(yǎng)鹽

    羅曉凡,魏 皓

    (天津科技大學 海洋科學與工程學院,天津 300457)

    碳元素(主要是CO2)在大氣圈、水圈、巖石圈、生物圈之間轉移和交換的過程稱為碳循環(huán)。由于人類活動的影響,全球海面上大氣CO2的年均濃度已由1980年的338.80μatm升高至2012年的392.53μatm,并保持持續(xù)攀升的狀態(tài)(http://www.esrl.noaa.gov/gmd/ccgg/trends/global.html)。在1958—2008年間,平均約43%的人為CO2停留在大氣中,其余被陸地和海洋吸收[1]。據最新觀測數據估算全球海洋每年吸收人為CO22.0±1.0PgC[2],一定程度上緩解了大氣CO2濃度的上升速度以及人為驅動的氣候變化幅度[3],但這同時也擾動了海洋碳酸鹽體系,導致了海水pH的下降趨勢,由此引發(fā)的海洋環(huán)境問題對海洋生態(tài)系統(tǒng)產生了潛在的影響[4-5]。海洋作為全球碳循環(huán)重要的組成部分,在保持其自身環(huán)境健康的前提下,能夠在吸收大氣CO2、減緩氣候變暖中起多大的作用,未來又將扮演什么樣的角色,這已然成為當前研究的熱點問題。

    陸架海是陸地與大洋的連接帶,受人類活動影響最為顯著,并且具有復雜的潮波系統(tǒng)、環(huán)流結構、季節(jié)性河流輸入、較高的生物生產力等特點,各種過程變化劇烈,空間梯度大。陸架海的面積雖僅占全球海洋總面積的8%,其初級生產卻是全球總初級生產的19%~28%[6],因此陸架海生態(tài)系統(tǒng)的固碳作用對全球碳收支平衡有不可忽略的貢獻,尤其是“碳失匯”的問題被提出后,陸架海對大氣CO2的源/匯作用成為平衡全球CO2收支不確定性因素之一,引起熱議[7-11]。Takahashi[2]給出的全球氣候態(tài)年均凈海-氣 CO2交換的數據中,變化迅速且復雜的近岸水體數據相對匱乏。IPCC-AR4所使用的氣候-碳循環(huán)耦合模式比較計劃的多模式數據(Coupled Climate-Carbon Cycle Model Intercomparison Project,C4MIP)也沒有包含時空變化劇烈的陸架海過程,給氣候預測帶來了不確定性[12]。而為了更準確地理解海洋在大氣CO2濃度不斷上升情景下的作用與響應,陸架海及近岸水體碳通量的量化研究得到越來越多的關注[7-11]。

    本研究綜述了近年來對陸架海CO2源/匯格局的爭論,指出通過建立3維-物理-生物-碳化學耦合數值模式量化陸架海碳循環(huán)關鍵過程的科學發(fā)展趨勢?;陉懠芎L佳h(huán)模式的研究現狀,從物理場、生物泵、碳化學系統(tǒng)等方面總結了建模中應考慮的關鍵控制過程,提出了構建東中國陸架海碳循環(huán)模式的基本思路以及擬解決的科學問題。

    1 陸架海對大氣CO2的源/匯作用

    陸架海對大氣CO2的源/匯作用是碳循環(huán)研究中最基礎也是最重要的部分,它的強弱反映陸架海對氣候變化的調節(jié)能力,目前還處于熱烈的爭論之中。部分學者認為全球陸架海是向大氣釋放CO2的源[13-15],其根據是河流輸入的有機碳通量大于海洋有機碳埋藏速率,剩余部分很有可能被分解為CO2返回大氣。Tsunogai等[16]在1999年基于東海觀測結果提出了“陸架泵”的概念,認為如果以東中國海吸收大氣CO2的效率去估算,全球陸架海每年可吸收1PgC。Yool等[17]采用GCM模型量化“陸架泵”作用探究其在全球陸架海的潛在強度并指出“陸架泵”效率在不同海區(qū)變化不一,外推估算全球陸架海每年可吸收0.6PgC,但該模型粗糙的分辨率對研究陸架海問題存在很大局限性與不確定性。Thomas等[18]與Bozec等[19]使用觀測數據對北海的碳通量及碳輸出的研究證實并支持了“陸架泵”假說。Bates等[20]指出陸架碳泵作用使得楚科奇海陸架成為大氣CO2較強的碳匯海區(qū)。然而,陸架海是多過程共同作用的海區(qū),其物理、生物、化學機制不盡相同,因此Cai等[7]根據陸架海的流場、形態(tài)、緯度等進行分類研究并指出中高緯陸架海區(qū)是大氣CO2匯,吸收CO2約0.33PgC yr-1;低緯陸架海區(qū)是大氣CO2源,釋放CO2約0.11PgC yr-1。Chen等[10]則搜集了大量前人關于陸架海、河口、鹽沼濕地及紅樹林等海-氣CO2通量的論文,給出一個認識:陸架海是大氣CO2匯,近岸生態(tài)系統(tǒng)是大氣CO2源。但Xue等[21]指出地處中緯度的開闊陸架?!宵S海中部缺少“陸架泵”需要的物理環(huán)境進行有效的碳輸出,其年均尺度上體現為大氣CO2凈源,這與早期Cai等[7]和Chen等[10]的結論有出入。這些爭論主要源于陸架海表層海水CO2分壓(pCO2)強烈的時間變化與巨大的空間分布不均勻性,同時,一次調查數據不能絕對代表一個月份或是一個季節(jié)的特征,這便是實測數據的局限性所在,而“陸架泵”假說是否成立,還需結合數值模式進行驗證。此外,Chen等[11]總結了全球106個河口的海-氣CO2交換,指出大多數河流與河口都是大氣CO2源,但是大河輸送的營養(yǎng)鹽在河流羽流區(qū)被用來進行光合作用,一定程度上使水體呈現CO2不飽和狀態(tài),因此,大河的羽流區(qū),甚至是距離河口幾百公里遠的海區(qū)都可能成為CO2匯。Tseng等[22]認為長江沖淡水流量影響著CO2的吸收,三峽大壩建成后,長江沖淡水流量減小,可能導致CO2吸收能力的下降,甚至會使東海陸架由大氣CO2匯轉為源,這個結論有待驗證。Dai等[23]根據對南海的研究提出了大洋主控型邊緣海碳循環(huán)概念,認為來自大洋深部的外源溶解無機碳通過上升流和躍密混合被輸送至真光層后,如果不被同樣來源的營養(yǎng)鹽所支持的有機碳生產完全消耗,其過剩部分就會溢出至大氣中,而在河流主控型邊緣海,陸源有機物與營養(yǎng)鹽輸入含量的競爭調控該海區(qū)生物的新陳代謝作用,從而影響CO2通量??梢姡煌懠芎^(qū)水文環(huán)境、動力過程存在差別,不同生態(tài)系統(tǒng)生物泵運作不同,陸架海源/匯格局的分布各具特色。碳酸鹽系統(tǒng)實測資料的匱乏及其時空局限性,限制了人們對海洋碳循環(huán)整體的進一步認識,這使得海洋碳循環(huán)模式快速發(fā)展,同時也促進了人們對海洋碳循環(huán)關鍵過程進行量化的研究。

    2 陸架海碳循環(huán)模式的發(fā)展

    相較于大洋,陸架海系統(tǒng)的復雜性使得陸架海碳循環(huán)模式的研究更為困難,目前,國際上對陸架海碳循環(huán)模式的報道不多,我國對東中國陸架海3維-物理-生物-碳循環(huán)模式的研究尚在起步階段,陸架海碳循環(huán)模式的建立已成為各國氣候預測模型急需攻克的難題。

    自20世紀60年代起,大洋碳循環(huán)模式經歷了由簡單到復雜的發(fā)展。從早期的箱式收支模式到三維海洋無機碳循環(huán)模式,再到基于海洋動力學耦合生物地球化學過程的三維碳循環(huán)模式,直至今天的氣候-碳循環(huán)耦合模式。然而直至20世紀90年代,陸架海碳循環(huán)模式才開始發(fā)展。CO2在海洋中的轉移和歸宿,即海洋吸收、轉移大氣CO2的能力及CO2在海洋中的循環(huán)機制均集中體現為其物理、化學和生物過程[24]。這也是陸架海碳循環(huán)模式的三大模塊,三者的耦合完成了溶解無機碳(DIC)、溶解有機碳(DOC)、顆粒有機碳(POC)等在陸架海的遷移轉化。

    本研究在Hofmann等[25]對陸架海碳循環(huán)模式綜述的基礎上,列出近年來有關模擬陸架海碳收支的文章(表1),雖然不同海區(qū)模式所含變量以及邊界條件不同,但模式基本包括了陸源輸入的淡水、營養(yǎng)鹽、碳,海-氣交換的CO2及大氣沉降的營養(yǎng)鹽,沉積物-水界面交換通量,生物作用下營養(yǎng)鹽及碳的遷移轉化,跨陸架輸運及陸架海與外海的交換。同時,基于地球物理流體力學原理所建立的海洋環(huán)流模式提供了更為真實,結構更復雜的海洋物理場,使得對陸架海碳循環(huán)的模擬也更為準確。

    表1 部分研究陸架邊緣海海-氣交換通量及碳輸出的模式Table 1 Sea-Air CO2flux(SACF)and carbon transportation models developed for the continental shelf seas

    3 陸架海碳循環(huán)模式中的關鍵過程

    3.1 關鍵物理過程

    現有的數值模式對海洋碳循環(huán)的模擬都存在缺點和不足之處,Doney等[38]對海洋碳循環(huán)模式比較計劃——階段2(the phase 2of the Ocean Carbon-cycle Model Intercomparison Project,OCMIP-2)的13個全球碳循環(huán)模式的水動力場以及子午圈翻轉流進行比較指出,這些缺點與不足主要與物理環(huán)流有關,碳循環(huán)模式結果在生物地球化學過程的差異有40±25%都是來自海洋環(huán)流模式[38]。生源要素必然隨海水運動而遷移,除了由海洋各個界面輸入、輸出外,海洋環(huán)流是決定各種物質時空分布的主要因素[39]。因此,陸架海碳循環(huán)模式對海洋動力模式有較高的要求,同時也促進了海洋模式向高分辨率、多過程耦合發(fā)展。陸架海多上升、下降流系統(tǒng),其垂直方向輸送的營養(yǎng)鹽以及碳通量也影響著表層pCO2的分布[27]。在南海,上升流是導致越南沿岸與呂宋島西北海區(qū)pCO2空間分布差異的因素之一[30],南海北部陸架海區(qū)碳酸鹽系統(tǒng)的動態(tài)格局也主要由河流羽流和近岸上升流影響下的物理和生物強迫控制[40]。此外,在有利于下降流形成的風驅動下,沿岸流和底層離岸輸運有機碳也是陸架海碳循環(huán)模型中需要考慮的物理過程。

    海洋的暖化降低了海水CO2溶解度,可能使海洋碳存儲減小,從而對大氣CO2濃度起正反饋作用[41],加劇全球氣候變化。因此除環(huán)流輸運外,對海水溫度長期變化的模擬、抓住區(qū)域海溫度對氣候變化的響應特征是非常重要的。Chai等[30]模式結果表明水溫對海水pCO2季節(jié)變化的影響比生物活動更強,是控制南海pCO2的主導因素。已有研究表明東中國陸架海增溫更加顯著[42-43],由此而帶來的對碳循環(huán)的影響還需科學評估。

    在歐洲北海,顯著不同的南北兩個生物地球化學省的海-氣CO2交換是溫度與生物作用影響平衡的結果,而這很大程度上取決于水體的層化以及在層化影響下生物作用的差異[31]。夏末,北海南部水體混合均勻,呈非自養(yǎng)型,向大氣凈釋放CO2,而在出現季節(jié)性分層的北海北部,表層水體呈自養(yǎng)型,該區(qū)是大氣CO2的匯[19]??梢?,陸架海碳循環(huán)模式對海洋模式中混合的參數化方案提出了更精細的要求,以能夠反演時空快速變化的層化/混合、鋒面、河流羽流結構等。

    3.2 生物泵

    生物泵作用是陸架海碳循環(huán)中重要的環(huán)節(jié),碳從真光層至深層的傳遞要依靠生物的作用。在陸架海,大約一半的有機碳通過生物泵作用轉移到深海[25]。生物過程驅動了海洋對CO2的吸收,Kühn等[32]和Lorkowski等[36]對北海碳收支的模擬中指出,如果模式不包含生物過程,北海北部將會和北海南部一樣表現為大氣CO2的源,這與實際調查結果不符。在羅斯海,海冰消融的海區(qū)因較高的生物初級生產而成為大氣CO2的強匯區(qū)[29]。因此陸架海碳循環(huán)模式耦合低營養(yǎng)級生態(tài)模式是當前碳循環(huán)模式發(fā)展主流,包含詳細的生物過程與營養(yǎng)鹽循環(huán)的碳循環(huán)模式可以更準確的對碳化學系統(tǒng)收支進行估算。Borges等[44]指出在具有高初級生產力的近岸環(huán)境,碳化學系統(tǒng)對富營養(yǎng)化的響應要強于對酸化的響應。此外,底棲生物過程影響著底層營養(yǎng)鹽的循環(huán)以及堿度的變化,模式中對碳化學系統(tǒng)估算的不確定性有一部分是來自底棲作用[4]。在經典生物泵中,微生物通常被認為是分解者,而研究表明,細菌還可以把活性DOC高效地轉化為惰性DOC,形成巨大的惰性DOC碳庫,進入的碳循環(huán)的慢性循環(huán)(千年尺度)中,起到儲碳作用,這一過程就是“微生物碳泵”[45]。目前,雖然對微生物碳泵有了理論的認識,但將這一過程進行模式化,還是面臨很大的挑戰(zhàn)。復雜的生物過程不僅增加了模式的難度,同時也增加了誤差在模式中的傳遞。此外,碳循環(huán)與營養(yǎng)鹽循環(huán)都有其獨特的控制過程,模式中應考慮將兩者解耦,這樣也可以使浮游植物的生長不受碳而只受營養(yǎng)鹽的限制[4,36]。

    3.3 碳化學系統(tǒng)

    3.3.1 界面碳通量

    在現有日趨完善的生態(tài)動力學模型的發(fā)展基礎之上耦合碳化學系統(tǒng)是當前陸架海碳循環(huán)模式的發(fā)展方向,因此,對碳化學系統(tǒng)的整體認知是模式研究者建模中重要的環(huán)節(jié)。陸架海碳化學系統(tǒng)參與了海-氣界面、陸-海界面、沉積物-水界面以及水體內部的化學反應,調控著水體pH,并影響其他化學反應,對形成和維持海洋生命活動有重要作用。在海洋中,碳的存在形式主要為DIC、DOC、POC,且約98%的碳是以溶解無機碳的形式存在[12]。

    評估陸架海在全球氣候變化中的作用,海-氣CO2交換是最基礎的一環(huán),也是對碳酸鹽系統(tǒng)的平衡起重要作用的界面過程。海水碳酸鹽系統(tǒng)是由DIC、總堿度(Talk)、pH值、pCO2等通過基本的熱力學關系相聯(lián)系,包括了碳酸、硼酸、磷酸、硅酸、水等弱酸的電離平衡[46]。在模式中,為了減小其復雜度,通常只考慮碳酸-硼酸-水體系。目前,模式估算海-氣CO2通量(SACF)多使用公式(1):

    其中,fice為海冰覆蓋率;為本文章上述所指的pCO2,即海水表層CO2分壓;為海面上大氣CO2分壓;α為CO2的溶解度,模式中通常使用Weiss[47]給出的與溫鹽相關的經驗公式;k為液相傳輸速率,模式中通常使用Wanninkhof[48]給出的與海面10m風速以及溫度相關的經驗公式。盡管使用上述方法存在不確定性,主要體現在氣體交換系數(k﹒α),但使用該方法估算SACF是被普遍接受的[46]。海-氣界面是碳循環(huán)與人類活動影響直接相關的關鍵界面,Omstedt等[49]對波羅的海pCO2長期變化的模擬表明,工業(yè)革命前波羅的海表現為凈源作用,而工業(yè)革命后,pCO2季節(jié)變化增強,年際變異上源/匯作用都有出現;同時對pH的長期模擬指出,在工業(yè)化之前pH維持一個穩(wěn)定狀態(tài),而工業(yè)化后,隨著大氣CO2濃度的增加,pH有小幅下降。在地中海,1 980s其表現為大氣CO2的凈源,每年釋放1.5×1012mol C,而到了2 000s,大氣CO2濃度較1 980s增加了30μatm,使得大氣向海洋的CO2通量增加,地中海表現出弱匯作用,平均每年吸收0.2×1012mol C[50]。Lu等[51]對南海SEATS站位建立1維模型的模擬結果表明,海-氣CO2交換調控pCO2的季節(jié)變化,甚至超過SST的影響,且因其為寡營養(yǎng)鹽海區(qū),生物過程的貢獻較??;在日變化上,pCO2的調控因子與局地物理和生物環(huán)境有關。

    與大洋不同,陸架海直接受人類擾動的另一個界面就是陸-海界面,主要表現為河流輸入帶來的影響。河流碳通量在特定時間和地區(qū)是單向輸送的,主要包括DIC、DOC、POC、PIC(顆粒無機碳)的輸送[52]。2004—2008年,長江平均每年向東海輸送DIC為1.56×1013g[53]。Wang等[54]研究表明2009年,黃海和長江向海輸送的DOC分別為3.20×1010g和1.58×1012g,POC分別為3.89×1011g和1.52×1012g,可能對河口及其近岸的生物化學循環(huán)和生態(tài)系統(tǒng)產生重要影響。Raymond等[55]指出當考慮近岸區(qū)域DIC和Talk的季節(jié)、年際、長期變化趨勢時,河流輸入的準確性是重要的。但在模式中,河流碳輸入資料的匱乏,時間變化的不確定性,模式分辨率不足等都可能使得河口低鹽區(qū)的模擬不甚理想[4,25],另外河口過程、濕地等增加了問題的復雜性,陸-海界面過程是陸架海碳循環(huán)模式研究的一個難點。由于河流帶來流域植被代謝的碳、營養(yǎng)物質以及農耕區(qū)的化肥等,因此,將陸架海碳循環(huán)模型與陸地碳循環(huán)模型、水文模型銜接,獲得連續(xù)的河流邊界條件是河流主控型陸架海碳循環(huán)模型的一個發(fā)展方向。此外,海底地下水(Submarine Groundwater Discharge-SGD)也是陸地向海洋輸送水量和物質的重要通道,被認為是近海營養(yǎng)鹽、碳等溶解物質的源[56]。全球尺度上,SGD碳通量輸送粗略估算約為0.23PgC yr-1,相當于全球河流碳輸送的30%[57]。在南海北部,SGD雖占珠江徑流量的12%~21%,但其攜帶的DIC可能高達(153~347)×109mol yr-1,相當于珠江DIC輸入量的23%~53%[58]。可見,SGD對近岸碳收支平衡的潛在影響。然而現今對SGD碳通量的研究還不全面,在陸架海碳循環(huán)模式中對SGD的定量化是比較困難的,SGD數值模式的發(fā)展雖然為研究SGD開拓了視野,但目前可以用來對模式進行有效驗證的實測數據及經驗數據非常少,由此而帶來的不確定性限制了模式的廣泛應用與發(fā)展[59]。

    沉積物-水界面交換體現在有機碳沉積與再懸浮以及有機碳在微生物作用下的再礦化作用,在動力較強、海底有機污染嚴重的陸架海,這個界面尤其重要。大洋沉積通過對有機碳和CaCO3的埋藏可以作為永久碳庫,然而在東中國陸架海,海底是非河流影響區(qū)域重要的營養(yǎng)鹽來源,強烈的潮流、風場的擾動使得顆粒物再懸浮嚴重,加之活躍的底棲生物作用,有機碳的保存效率大幅降低[60-62]。由于該界面精細過程觀測難以實施,僅有零星底界面交換實驗數據,而交換又與動力過程密切相關,目前還沒有一個完善的底界面交換模型,這個過程的參數化是陸架碳循環(huán)模型中另一難點,需要突破。

    陸架海-外海界面碳通量和營養(yǎng)鹽通量受開邊界條件的影響。Kühn等[32]指出使用來自二氧化碳信息處理中心(Carbon Dioxide Information Analysis Center,CDIAC)的氣候態(tài)DIC、Talk數據作為北海的開邊界條件會高估北海DIC濃度。Artioli等[4]指出完善初始條件與邊界條件,尤其是河流輸入以及與波羅的海通量交換,可以顯著減小對北海營養(yǎng)鹽模擬的偏差。在中國東海,高溫高鹽的黑潮水入侵為陸架海輸送了可觀的營養(yǎng)鹽,其中硝酸鹽輸送入東海的平均值為170.8±41.6kmol s-1,約為磷酸鹽輸送量的13.6倍[63],這對于陸架海營養(yǎng)鹽循環(huán)、生物生產以及整個碳循環(huán)系統(tǒng)都有重要影響。

    3.3.2 碳的遷移轉化

    在碳的遷移轉化過程中(圖1),生物作用主導有機碳與無機碳的相互轉化。而DOM的光降解主要表現為有色溶解有機物(CDOM)的光降解,陸源CDOM入海后通過光降解作用在近海海域發(fā)生轉化與消除,潛在地影響海洋生態(tài)系統(tǒng)及海洋碳循環(huán)[64-66]。賀志鵬等[67]指出黃海表層至20m層DOC濃度呈上升趨勢可能就是表層DOC發(fā)生光降解作用的直接證據。Yang等[65]估算得黃東海每年約5.4~10.5TgC的DOC發(fā)生光礦化,約占初級生產的2.9%~5.7%。DIC以碳酸鹽形式移出水體主要表現為CaCO3的產出,包括珊瑚、顆石藻外殼、某些浮游動物的碳酸鹽外殼和骨架等。全球海洋碳酸鹽的生產和輸出以顆粒無機碳的形式創(chuàng)造了較高的碳輸出效率(0.4~1.8GtC yr-1)[68]。Cao等[69]指出CaCO3的形成與溶解調節(jié)著海洋 DIC系統(tǒng)并最終調控海洋對大氣CO2的緩解能力,在南海,近岸淺水海區(qū)CaCO3的溶解對外海上層水體中累積的Ca2+和總堿度有貢獻。近岸水體酸化很可能會減少碳酸鹽的生產,加速碳酸鹽基底的溶解,從而打破碳化學系統(tǒng)的平衡,不利于鈣化生物的生存。因此,對于碳酸鹽沉積底質的陸架海以及珊瑚生態(tài)系統(tǒng),在模式中應酌情考慮碳酸鹽的生成與溶解過程。對于POC的源匯項,除受陸源與生物作用影響外,還受到沉降及再懸浮過程的影響。海洋對大氣CO2的調節(jié)作用重點在于碳的垂直遷移并最終被埋藏于海底的過程。郭志剛等[70]等指出有機碳從水體向海底沉積物轉移并保存的條件要滿足4點:1)有充足的POC向海底轉移;2)有一定的沉積速率;3)適于細顆粒沉積的弱水環(huán)境;4)相對還原的海底環(huán)境。然而陸架海動力調整變化快,沉降并保存在沉積物中的POC在受到動力擾動時,便會發(fā)生再懸浮,使POC沉降受限。據統(tǒng)計,大陸架及大陸坡的沉積物中,有40%~85%的有機碳要發(fā)生再懸?。?1]。長江口及鄰近海域海-氣交換的CO2僅有13%以顆粒有機碳的形式垂直轉移而形成表層沉積物的有機碳,其余87%沒有被埋藏的碳通過生物泵或動力過程進入水中再循環(huán)[72]。

    圖1 陸架海碳的遷移轉化示意圖Fig.1 A sketch map showing the migration and transformation of carbon in the continental shelf seas

    4 東中國陸架海碳循環(huán)模式的構建思想

    東中國海地處東亞季風區(qū),渤海、黃海的全部以及東海的2/3都屬于寬闊的陸架海,受到長江、黃河等陸地徑流的強烈影響以及臺灣暖流、黑潮入侵等帶來的強烈作用。目前,基于水動力模型的營養(yǎng)鹽-低營養(yǎng)級生態(tài)動力學模型在東中國海已取得較好的應用[73],在此基礎上耦合碳化學系統(tǒng)是碳循環(huán)模型發(fā)展的基本方向。東中國陸架海碳循環(huán)模型的基本構想是在生態(tài)動力學模型中增加DIC作為主要的狀態(tài)變量,同時增加Talk、pCO2、SACF作為診斷變量,其中為減小模式復雜度,可以考慮將Talk與鹽度建立關系,pCO2通過碳酸鹽體系平衡計算[74],SACF的計算前文已給出(詳見3.3.1節(jié))。相較于大洋模型,陸架海碳循環(huán)模型顯著的特點在于陸-海界面和沉積物-水界面碳通量的參數化。目前對河流碳通量動態(tài)變化規(guī)律的研究并不細致,但河流徑流量主導其無機碳和有機碳入海通量變化[75-76],由此可根據徑流量以及DIC、DOC、POC在不同年代入海通量的年際特征值推算碳通量長時間序列作為河流邊界條件。沉積物-水界面所表現的固碳及去碳作用,在年或年際尺度上處于動態(tài)平衡狀態(tài),季節(jié)尺度上有變化,這在現代陸架環(huán)境下主要與生物過程相關[70]。這一界面過程目前還沒有精細的數值模型可供采用,結合水動力過程的沉積物-水界面物質交換、再懸浮過程的模型量化,應該是東中國陸架海碳循環(huán)模型的一個主要特色。

    東中國陸架海3維-物理-生物-碳循環(huán)模型的建立,可進一步認識東中國海表層pCO2、海-氣碳通量、向陸架外碳輸出效率的季節(jié)與年際變化;定量評估東中國海源/匯動態(tài)格局以及生物泵、“陸架泵”的貢獻;河流徑流量變化所致的物理、生物效應對河口及其近岸碳通量的影響;沉積-水界面交換通量動力學參數化在碳循環(huán)中的作用;并預測未來陸架海碳循環(huán)系統(tǒng)的演變規(guī)律,為全球碳循環(huán)模型提供可嵌入的東中國陸架海碳循環(huán)模塊,這對于改善全球碳循環(huán)模型,減少其不確定性,準確評估東亞海區(qū)在氣候變化中的作用都是非常必要的。

    (References):

    [1] LEQUéRéC,RAUPACH M R,CANADELL J G,et al.Trends in the sources and sinks of carbon dioxide[J].Nature Geoscience,2009,2(12):831-836.

    [2] TAKAHASHI T,SUTHERLAND S C,WANNINKHOF R,et al.Climatological mean and decadal change in surface ocean pCO2,and net sea-air CO2flux over the global oceans[J].Deep-Sea Research II,2009,56*8-10):554-577.

    [3] FRIEDLINGSTEIN P,COX P,BETTS R,et al.Climate-carbon cycle feedback analysis:results from the C4MIP model intercomparison[J].Journal of Climate,2006,19(14):3337-3353.

    [4] ARTIOLI Y,BLACKFORD J C,BUTENSCH?N M,et al.The carbonate system in the North Sea:Sensitivity and model validation[J].Journal of Marine Systems,2012,102-104:1-13.

    [5] MCGRATH T,KIVIMāE C,TANHUA T,et al.Inorganic carbon and pH levels in the Rockall Trough 1991-2010[J].Deep-Sea Research I,2012,68:79-91.

    [6] LONGHURST A.Seasonal cycles of pelagic production and consumption[J].Progress in Oceanography,1995,36:77-167.

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

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

    [9] GUO X H,DAI M H,ZHAI W D,et al.CO2flux and seasonal variability in a large subtropical estuarine system,the Pearl River Estuary,China[J].Journal of Geophysical Research,2009,114(G3):G03013.

    [10] CHEN C T A,BORGES A V.Reconciling opposing views on carbon cycling in the coastal ocean:Continental shelves as sinks and nearshore ecosystems as sources of atmospheric CO2[J].Deep-Sea Research II,2009,56(8-10):578-590.

    [11] CHEN C T A,HUANG T H,FU Y H,et al.Strong sources of CO2in upper estuaries become sinks of CO2in large river plumes[J].Current Opinion in Environmental Sustainability,2012,4:179-185.

    [12] IPCC.Summary for Policymakers of the Synthesis Report of the IPCC Fourth Assessment Report[M].Cambridge,UK:Cambridge University Press,2007.

    [13] VER M B,MACKENZIE F T,LERMAN A.Biogeochemical responses of the carbon cycle to natural and human perturbations:Past,present,and future[J].American Journal of Science,1999,299:762-801.

    [14] MACKENZIE F T,VER L M,LERMAN A.Coastal-zone biogeochemical dynamics under global warming[J].International Geology Review,2000,42(3):193-206.

    [15] FASHAM M J R,BALINO B M ,BOWLES C,et al.A new vision of ocean biogeochemistry after a decade of the Joint Global Ocean Flux Study(JGOFS)[J].Ambio(Special Issue),2001,10:4-30.

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

    [17] YOOL A,FASHAM M J R.An examination of the"continental shelf pump"in an open ocean general circulation model[J].Global Biogeochemical cycles,2001,15(4):831-844.

    [18] THOMAS H,BOZEC Y,ELKALAY K,et al.Enhanced open ocean storage of CO2from shelf sea pumping[J].Science,2004,304(5673):1005-1008.

    [19] BOZEC Y,THOMAS H,ELKALAY K,et al.The continental shelf pump for CO2in the North Sea-evidence from summer observation[J].Marine Chemistry,2005,93:131-147.

    [20] BATES N R.Air-sea CO2fluxes and the continental shelf pump of carbon in the Chukchi Sea adjacent to the Arctic Ocean[J].Journal of Geophysical Research,2006,111(C10):C10013.

    [21] XUE L,ZHANG L J,CAI W J,et al.Air-sea CO2fluxes in the southern Yellow Sea:An examination of the continental shelf pump hypothesis[J].Continental Shelf Research,2011,31(18):1904-1914.

    [22] TSENG C M,LIU K K,GONG G C,et al.CO2uptake in the East China Sea relying on Changjiang runoff is prone to change[J].Geophysical Research Letters,2011,38(24):L24609.

    [23] DAI M H,CAO Z M,GUO X H,et al.Why are some marginal seas sources of atmospheric CO2?[J].Geophysical Research Letters,2013,40(10):2154-2158.

    [24] XU Y F,PU Y F,ZHAO L.Advances in the studies of ocean carbon-cycle model[J].Advances in Earth Science,2005,20(10):1106-1115.徐永福,浦一芬,趙亮.海洋碳循環(huán)模式的研究進展[J].地球科學進展,2005,20(10):1106-1115.

    [25] HOFMANN E E,CAHILL B,FENNEL K,et al.Modeling the Dynamics of Continental Shelf Carbon[J].Annual Review of Marine Science,2011,3:93-122.

    [26] WALSH J J,DIETERLE D A.CO2cycling in the coastal ocean.I-A numerical analysis of the southeastern Bering Sea with applications to the Chukchi Sea and the northern Gulf of Mexico[J].Progress in Oceanography,1994,34(4):335-392.

    [27] IANSON D,ALLEN S E.A two-dimensional nitrogen and carbon flux model in a coastal upwelling region[J].Global Biogeochemical Cycles,2002,16(1):1-16.

    [28] FENNEL K,WILKIN J,PREVIDI M,et al.Denitrification effects on air-sea CO2flux in the coastal ocean:Simulations for the northwest North Atlantic[J].Geophysical Research Letters,2008,35(24):L24608.

    [29] ARRIGO K R,DIJKEN G V,LONG M.Coastal Southern Ocean:A strong anthropogenic CO2sink[J].Geophysical Research Letters,2008,35(1):L21602.

    [30] CHAI F,LIU G M,XUE H J.Seasonal and Interannual Variability of Carbon Cycle in South China Sea:A Three-Dimensional Physical-Biogeochemical Modeling Study[J].Journal of Oceanography,2009,65:703-720.

    [31] PROWE A E F,THOMAS H,PāTSCH J,et al.Mechanisms controlling the air-sea CO2flux in the North Sea[J].Continental Shelf Research,2009,29:1801-1808.

    [32] KüHN W,PāTSCH J,THOMAS H,et al.Nitrogen and carbon cycling in the North Sea and exchange with the North Atlantic-A model study,Part II:Carbon budget and fluxes[J].Continental Shelf Research,2010,30:1701-1716.

    [33] DRUON J N,MANNINO A,SIGNORINI S,et al.Modeling the dynamics and export of dissolved organic matter in the Northeastern U.S.continental shelf[J].Coastal and Shelf Science,2010,88(4):488-507.

    [34] KUZNETSOV I,NEUMANN T.Simulation of carbon dynamics in the Baltic Sea with a 3Dmodel[J].Journal of Marine Systems,2013,111-112:167-174.

    [35] W?HLSTR?M I,OMSTEDT A,BJ?RK G,et al.Modelling the CO2dynamics in the Laptev Sea,Arctic Ocean:Part I[J].Journal of Marine Systems,2012,102-104:29-38.

    [36] LORKOWSKI I,PāTSCH J,MOLL A,et al.Interannual variability of carbon fluxes in the North Sea from 1970to 2006-Competing effects of abiotic and biotic drivers on the gas-exchange of CO2[J].Estuarine,Coastal and Shelf Science,2012,100:38-57.

    [37] ?RTHUN M,BELLERBY R G J,OMAR A M,et al.Spatiotemporal variability of air-sea CO2fluxes in the Barents Sea,as determined from empirical relationships and modeled hydrography[J].Journal of Marine Systems,2012,98-99:40-50.

    [38] DONEY S C,LINDSAY K,CALDEIRA K,et al.Evaluating global ocean carbon models:The importance of realistic physics[J].Global Biogeochemical Cycles,2004,18(3):GB3017.

    [39] SUN W X,JIANG W S,LI L.Numerical Models of Hydrodynamics in the Coastal Environment[M].Beijing:Science Press,2004.孫文心,江文勝,李磊.近海環(huán)境流體動力學數值模型[M].北京:科學出版社,2004.

    [40] CAO Z M,DAI M H,ZHENG N,et al.Dynamics of the carbonate system in a large continental shelf system under the influence of both a river plume and coastal upwelling[J].Journal of Geophysical Research,2011,116(G2):G02010.

    [41] DONEY S C,TILBROOK B,ROY S,et al.Surface-ocean CO2variability and vulnerability[J].Deep-Sea Research II,2009,56(8-10):504-511.

    [42] WEI H,SHI J,LU Y Y,et al.Interannual and long-term hydrographic changes in the Yellow Sea during 1977—1998[J].Deep-Sea Research II,2010,57(11-12):1025-1034.

    [43] LI J X,WEI H,LUO X F.Spatial and temporal inter-annual variations of sea surface temperature in the continental shelf of the East China Sea[J].Oceanologia et Limnologia Sinica,2012,43(3):624-634.李家星,魏皓,羅曉凡.東海陸架表層水溫年際變化時空特征分析[J].海洋與湖沼,2012,43(3):624-634.

    [44] BORGES A V,GYPENS N.Carbonate chemistry in the coastal zone responds more strongly to eutrophication than to ocean acidification[J].Limnology & Oceanography,2010,55(1):346-353.

    [45] JIAO N Z.Carbon fixation and sequestration in the ocean,with special reference to the microbial carbon pump[J].Scientia Sinica Terrae,2012,42(10):1473-1486.焦念志.海洋固碳與儲碳-并論微生物在其中的重要作用[J].中國科學:地球科學,2012,42(10):1473-1486.

    [46] XU Y F,ZHAO L,PU Y F,et al.Uncertainties in the estimate of the air-sea exchange flux of carbon dioxide[J].Earth Science Frontiers,2004,11(2):565-571.徐永福,趙亮,浦一芬,等.二氧化碳海氣交換通量估計的不確定性[J].地學前緣,2004,11(2):565-571.

    [47] WEISS R F.Carbon dioxide in water and seawater:the solubility of a non-ideal gas[J].Marine Chemistry,1974,2:203-215.

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

    [49] OMSTEDT A,GUSTAFSSON E,WESSLANDER K.Modelling the uptake and release of carbon dioxide in the Baltic Sea surface water[J].Continental Shelf Research,2009,29(7):870-885.

    [50] TAILLANDIER V,ORTENZIO F D,ANTOINE D.Carbon fluxes in the mixed layer of the Mediterranean Sea in the1980sand the2000s[J].Deep-Sea Research I,2012,65:73-84.

    [51] LU Z M,GAN J P,DAI M H.Modeling seasonal and diurnal pCO2variations in the northern South China Sea[J].Journal of Marine Systems,2012,92(1):30-41.

    [52] GAO Q Z,SHEN C D.Riverine carbon flux and continental erosion[J].Advances in Earth Science,1998,13(4):369-375.高全洲,沈承德.河流碳通量與陸地侵蝕研究[J].地球科學進展,1998,13(4):369-375.

    [53] LI D.The Study on the Hydro-Chemical Characteristics and the Flux to the sea about the Rivers in the East of China[D].Shanghai:East China Normal University,2009.李丹.中國東部若干入海河流水化學特征與入海通量研究[D].上海:華東師范大學,2009.

    [54] WANG X C,MA H Q,LI R H,et al.Seasonal fluxes and source variation of organic carbon transported by two major Chinese Rivers:The Yellow River and Changjiang(Yangtze)River[J].Global Biogeochemical Cycles,2012,26:GB2025.

    [55] RAYMOND P A,OH N H,TUMER R E,et al.Anthropogenically enhanced fluxes of water and carbon from the Mississippi River[J].Nature,2008,451:449-452.

    [56] MOORE W S.The role of submarine groundwater discharge in coastal biogeochemistry[J].Journal of Geochemical Exploration,2006,88(1-3):389-393.

    [57] ZHANG J,MANDAL A K.Linkages between submarine groundwater systems and the environment[J].Current Opinion in Environmental Sustainability,2012,4(2):219-226.

    [58] LIU Q,DAI M H,CHEN W,et al.How significant is submarine groundwater discharge and its associated dissolved inorganic carbon in a river-dominated shelf system?[J].Biogeosciences,2012,9:1777-1795.

    [59] SANTOS I R,EYRE B D,HUETTEL M.The driving forces of porewater and groundwater flow in permeable coastal sediments:A review[J].Estuarine,Coastal and Shelf Science,2012,98:1-15.

    [60] KIM G,RYU J W,YANG H S,et al.Submarine groundwater discharge(SGD)into the Yellow Sea revealed by 228Ra and 226Ra isotopes:implications for global silicate fluxes[J].Earth and Planetary Science Letters,2005,237(1-2):156-166.

    [61] YUAN Y,WEI H,ZHAO L,et al.Observations of sediment resuspension and settling off the mouth of JiaozhouBay,Yellow Sea[J].Continental Shelf Research,2008,28(19):2630-2643.

    [62] HUNG J J,LIN P L,LIU K K.Dissolved and particulate organic carbon in the southern East China Sea[J].Continental shelf Research,2000,20(4-5):545-569.

    [63] GUO X Y,ZHU X H,WU Q S,et al.The Kuroshio nutrient stream and its temporal variation in the East China Sea[J].Journal of Geophysical Research,2012,117(C1):C01026.

    [64] GUO W D,CHENG Y Y,YU X X,et al.An overview of the Photochemistry of marine chromophoric dissolved organic matter[J].Marine Science Bulletin,2008,27(3):107-114.郭衛(wèi)東,程遠月,余翔翔,等.海洋有色溶解有機物的光化學研究進展[J].海洋通報,2008,27(3):107-114.

    [65] YANG G P,REN C Y,LU X L,et al.Distribution,flux,and photoproduction of carbon monoxide in the East China Sea and Yellow Sea in spring[J].Journal of geophysical research,2011,116(C2):C02001.

    [66] XIE H X,AUBRY C,BéLANGER S,et al.The dynamics of absorption coefficients of CDOM and particles in the St.Lawrence estuarine system:Biogeochemical and physical implications[J].Marine Chemistry,2012,128-129:44-56.

    [67] HE Z P,SONG J M,ZHANG N X.Analysis of biogeochemical characteristics of dissolved organic carbon in the southern Yellow Sea[J].Advances in Marine Science,2006,24(4):477-488.賀志鵬,宋金明,張乃星.南黃海溶解有機物的生物地球化學特征分析[J].海洋科學進展,2006,24(4):477-488.

    [68] BERELSON,W M,BALCH W M,NAJJAR R,et al.Relating estimates of CaCO3production,export,and dissolution in the water column to measurements of CaCO3rain into sediment traps and dissolution on the sea floor:A revised global carbonate budget[J].Global Biogeochemical Cycles,2007,21:GB1024.

    [69] CAO Z M,DAI M H.Shallow-depth CaCO3dissolution:Evidence from excess calcium in the South China Sea and its export to the Pacific Ocean[J].Global Biogeochemical Cycles,2011,25(2):GB2019.

    [70] GUO Z G,YANG Z S.Vertical transfer processes for the particulate carbon in the East China Sea[J].Oceanologia et Limnologia Sinica,1997,28(6):659-664.郭志剛,楊作升.東海顆粒碳的垂直轉移過程[J].海洋與湖沼,1997,28(6):659-664.

    [71] ZHANG N X,SONG J M,HE Z P.Biogeochemical mechanism of particulate organic carbon(POC)variations in sea waters[J].Acta Ecologica Sinica,2006,26(7):2328-2339.張乃星,宋金明,賀志鵬.海水顆粒有機碳(POC)變化的生物地球化學機制[J].生態(tài)學報,2006,26(7):2328-2339.

    [72] SONG J M.Biogeochemistry of China Seas[M].Jinan:Shandong Science & Technology Press,2004.宋金明.中國近海生物地球化學[M].濟南:山東科技出版社,2004.

    [73] ZHAO L.GUO X Y.Influence of cross-shelf water transport on nutrients and phytoplankton in the East China Sea:a model study[J].Ocean Science,2011,7(1):27-43.

    [74] KANTHA L H.A general ecosystem model for applications to primary productivity and carbon cycle studies in the global oceans[J].O-cean Modelling,2004,6(3-4):285-334.

    [75] DAI M H,YIN Z Q,MENG F F,et al.Spatial distribution of riverine DOC inputs to the ocean:an updated global synthesis[J].Current Opinion in Environmental Sustainability,2012,4:170-178.

    [76] LIU X C,SHEN H T,HUANG Q H.Concentration variation and flux estimation of dissolved inorganic nutrient from the changjiang river into its estuary[J].Oceanologia et Limnologia Sinica,2002,33(3):332-340.劉新成,沈煥庭,黃清輝.長江入河口區(qū)生源要素的濃度變化及通量估算[J].海洋與湖沼,2002,33(3):332-340.

    猜你喜歡
    陸架碳循環(huán)營養(yǎng)鹽
    瓊東南盆地鶯歌海組陸架邊緣軌跡演化及主控因素
    沉積學報(2023年1期)2023-04-29 18:57:49
    食物網與碳循環(huán)
    晚更新世以來南黃海陸架沉積物源分析
    海洋通報(2022年2期)2022-06-30 06:07:04
    涼水河子河營養(yǎng)鹽滯留能力評估
    層序地層格架內白云凹陷北坡珠江組沉積相分析
    南京城市系統(tǒng)碳循環(huán)與碳平衡分析
    多措并舉構筑綠色低碳循環(huán)發(fā)展的制造業(yè)體系
    瓊東海域冬季、夏季營養(yǎng)鹽結構特征及其對浮游植物生長的影響
    2012年冬季南海西北部營養(yǎng)鹽分布及結構特征
    南海北部陸架表層沉積物重礦物分布特征及物源意義
    91九色精品人成在线观看| 午夜福利成人在线免费观看| 99riav亚洲国产免费| 最近最新中文字幕大全免费视频| 亚洲精品美女久久av网站| 男女下面进入的视频免费午夜 | 波多野结衣av一区二区av| 丝袜在线中文字幕| 精品国产国语对白av| 夜夜躁狠狠躁天天躁| 国产av一区在线观看免费| www国产在线视频色| 最新在线观看一区二区三区| 国产av又大| 久久精品91无色码中文字幕| 一级片免费观看大全| 91av网站免费观看| 日韩三级视频一区二区三区| 亚洲精品美女久久久久99蜜臀| 午夜免费成人在线视频| 一区二区三区高清视频在线| 精品少妇一区二区三区视频日本电影| 国产黄片美女视频| 午夜激情av网站| 国产成+人综合+亚洲专区| 黄色毛片三级朝国网站| 一级毛片女人18水好多| 久久国产精品人妻蜜桃| 精品久久久久久成人av| 国产精品 欧美亚洲| 亚洲国产精品久久男人天堂| 久久久久精品国产欧美久久久| 亚洲男人的天堂狠狠| 特大巨黑吊av在线直播 | 国产一区二区三区视频了| 久久香蕉精品热| 久久精品91蜜桃| 精品国产超薄肉色丝袜足j| 中文亚洲av片在线观看爽| 国产精品1区2区在线观看.| 国产乱人伦免费视频| 国产精品98久久久久久宅男小说| 国产成+人综合+亚洲专区| 国产亚洲精品一区二区www| 色综合婷婷激情| 亚洲国产日韩欧美精品在线观看 | 国产亚洲精品av在线| 久久国产亚洲av麻豆专区| 看免费av毛片| 精品国内亚洲2022精品成人| 亚洲 欧美一区二区三区| 日本一本二区三区精品| 欧美久久黑人一区二区| 成人亚洲精品av一区二区| 国产av一区在线观看免费| 成人18禁在线播放| 黑人欧美特级aaaaaa片| 啦啦啦韩国在线观看视频| 国产精品一区二区三区四区久久 | 亚洲成人国产一区在线观看| 欧美午夜高清在线| 一区福利在线观看| 两个人免费观看高清视频| 1024香蕉在线观看| 女性生殖器流出的白浆| 亚洲久久久国产精品| 欧美激情极品国产一区二区三区| 露出奶头的视频| 成在线人永久免费视频| 久久精品成人免费网站| 亚洲五月婷婷丁香| 三级毛片av免费| 精品无人区乱码1区二区| 一进一出抽搐动态| 97超级碰碰碰精品色视频在线观看| 精品熟女少妇八av免费久了| 一级作爱视频免费观看| 午夜激情福利司机影院| 99久久精品国产亚洲精品| aaaaa片日本免费| 特大巨黑吊av在线直播 | 天天躁狠狠躁夜夜躁狠狠躁| 啦啦啦免费观看视频1| 丰满的人妻完整版| 国产欧美日韩精品亚洲av| 久久久久久大精品| 午夜亚洲福利在线播放| 欧美成人一区二区免费高清观看 | 99久久久亚洲精品蜜臀av| 亚洲av成人不卡在线观看播放网| 午夜福利18| 两个人免费观看高清视频| 国产一卡二卡三卡精品| 波多野结衣高清作品| 日韩视频一区二区在线观看| 国产黄色小视频在线观看| 国产成人影院久久av| 人人妻人人澡人人看| 十分钟在线观看高清视频www| 国产精品久久电影中文字幕| 欧美一级a爱片免费观看看 | 亚洲国产日韩欧美精品在线观看 | 午夜视频精品福利| av天堂在线播放| 日韩 欧美 亚洲 中文字幕| 丰满的人妻完整版| 国产欧美日韩精品亚洲av| 少妇裸体淫交视频免费看高清 | 欧美午夜高清在线| 精品久久久久久久毛片微露脸| а√天堂www在线а√下载| 黄片大片在线免费观看| 很黄的视频免费| 欧美中文综合在线视频| 国产精品 欧美亚洲| 脱女人内裤的视频| 久久精品国产综合久久久| 一二三四在线观看免费中文在| 身体一侧抽搐| 国产精品久久久久久亚洲av鲁大| 国产精品久久久av美女十八| 久久婷婷人人爽人人干人人爱| 天堂√8在线中文| 中文字幕另类日韩欧美亚洲嫩草| 一个人观看的视频www高清免费观看 | 无遮挡黄片免费观看| 国产三级在线视频| 最好的美女福利视频网| 一区二区日韩欧美中文字幕| 岛国视频午夜一区免费看| 长腿黑丝高跟| 制服丝袜大香蕉在线| 性色av乱码一区二区三区2| 久久亚洲真实| 最近最新免费中文字幕在线| 欧美乱妇无乱码| 女人被狂操c到高潮| 亚洲人成网站高清观看| 韩国av一区二区三区四区| 日韩精品中文字幕看吧| 曰老女人黄片| 99久久久亚洲精品蜜臀av| 午夜福利欧美成人| 精品乱码久久久久久99久播| 俄罗斯特黄特色一大片| 久久亚洲真实| 亚洲aⅴ乱码一区二区在线播放 | 男人的好看免费观看在线视频 | 亚洲aⅴ乱码一区二区在线播放 | 欧美成狂野欧美在线观看| 免费在线观看亚洲国产| 国产精品永久免费网站| 午夜福利在线在线| 成人亚洲精品一区在线观看| 精品欧美一区二区三区在线| 免费在线观看影片大全网站| 满18在线观看网站| 免费一级毛片在线播放高清视频| 日本 av在线| 久久九九热精品免费| 99久久99久久久精品蜜桃| 无遮挡黄片免费观看| 18美女黄网站色大片免费观看| 高潮久久久久久久久久久不卡| av中文乱码字幕在线| 色播亚洲综合网| 精品国产乱子伦一区二区三区| 最好的美女福利视频网| 亚洲成av人片免费观看| 亚洲最大成人中文| 波多野结衣高清作品| 99热6这里只有精品| 老熟妇乱子伦视频在线观看| 亚洲无线在线观看| 午夜久久久久精精品| 国内少妇人妻偷人精品xxx网站 | а√天堂www在线а√下载| 制服人妻中文乱码| 国产男靠女视频免费网站| 色精品久久人妻99蜜桃| 中文资源天堂在线| 亚洲久久久国产精品| 夜夜夜夜夜久久久久| 国产精品免费视频内射| 国产精品精品国产色婷婷| 国产久久久一区二区三区| ponron亚洲| 久久久久久久精品吃奶| 一级毛片女人18水好多| 国产免费av片在线观看野外av| 日本一区二区免费在线视频| 在线观看免费视频日本深夜| 亚洲欧美一区二区三区黑人| 亚洲av中文字字幕乱码综合 | 中文字幕久久专区| 草草在线视频免费看| 亚洲成人久久爱视频| 亚洲国产精品合色在线| 精品久久久久久久毛片微露脸| 欧美国产日韩亚洲一区| 成人国产一区最新在线观看| 在线观看免费日韩欧美大片| 日韩精品中文字幕看吧| 黑丝袜美女国产一区| 精品久久久久久久久久久久久 | 精品欧美国产一区二区三| 久久久久免费精品人妻一区二区 | 欧美成人免费av一区二区三区| 999久久久精品免费观看国产| 日韩欧美免费精品| 国产激情欧美一区二区| 久久精品国产清高在天天线| 老司机福利观看| 亚洲国产欧美日韩在线播放| 一级片免费观看大全| 成人精品一区二区免费| 国产精品一区二区三区四区久久 | 日韩精品青青久久久久久| 午夜福利在线观看吧| 人人妻人人看人人澡| 他把我摸到了高潮在线观看| 欧美一级毛片孕妇| 精品人妻1区二区| 丝袜人妻中文字幕| 美女国产高潮福利片在线看| 亚洲狠狠婷婷综合久久图片| www.熟女人妻精品国产| 天天躁狠狠躁夜夜躁狠狠躁| 99久久国产精品久久久| 身体一侧抽搐| 亚洲国产欧美一区二区综合| 91九色精品人成在线观看| 在线免费观看的www视频| 欧美丝袜亚洲另类 | 亚洲国产精品合色在线| 极品教师在线免费播放| 日日摸夜夜添夜夜添小说| 国产蜜桃级精品一区二区三区| 国内精品久久久久久久电影| 欧美一区二区精品小视频在线| 老鸭窝网址在线观看| 波多野结衣巨乳人妻| 91麻豆精品激情在线观看国产| 亚洲免费av在线视频| 日本黄色视频三级网站网址| a级毛片在线看网站| 在线观看舔阴道视频| 亚洲精品色激情综合| 婷婷丁香在线五月| 日韩欧美国产一区二区入口| 中文字幕人妻熟女乱码| 在线观看www视频免费| 国产高清有码在线观看视频 | 久久久久国产精品人妻aⅴ院| 国产精品日韩av在线免费观看| 亚洲第一电影网av| 日韩 欧美 亚洲 中文字幕| 三级毛片av免费| 日本黄色视频三级网站网址| 视频在线观看一区二区三区| 老司机在亚洲福利影院| 日日干狠狠操夜夜爽| 国产激情偷乱视频一区二区| 亚洲av成人av| 国产精品一区二区免费欧美| 黄色女人牲交| 亚洲国产中文字幕在线视频| 巨乳人妻的诱惑在线观看| 在线观看免费午夜福利视频| 日韩中文字幕欧美一区二区| 国产精品久久视频播放| 午夜激情av网站| 黑丝袜美女国产一区| 国内精品久久久久精免费| 两个人免费观看高清视频| 国产精品亚洲av一区麻豆| 一区二区三区高清视频在线| 精品欧美国产一区二区三| 黄色成人免费大全| 日韩av在线大香蕉| 日韩精品中文字幕看吧| а√天堂www在线а√下载| 男男h啪啪无遮挡| 久久精品国产清高在天天线| 国产亚洲欧美精品永久| 亚洲国产精品sss在线观看| bbb黄色大片| 精品一区二区三区视频在线观看免费| 一边摸一边抽搐一进一小说| 一级a爱视频在线免费观看| 久久人妻福利社区极品人妻图片| 国产亚洲精品av在线| 99在线人妻在线中文字幕| 在线播放国产精品三级| 国产av一区二区精品久久| 久久草成人影院| 波多野结衣高清无吗| 最近在线观看免费完整版| 又紧又爽又黄一区二区| 麻豆av在线久日| 欧美日韩瑟瑟在线播放| 亚洲欧美一区二区三区黑人| 韩国av一区二区三区四区| 久久香蕉国产精品| 日日爽夜夜爽网站| 国产三级黄色录像| 欧美最黄视频在线播放免费| 久久久水蜜桃国产精品网| 最近最新中文字幕大全电影3 | 巨乳人妻的诱惑在线观看| xxxwww97欧美| 2021天堂中文幕一二区在线观 | 男女下面进入的视频免费午夜 | 精品电影一区二区在线| 可以在线观看的亚洲视频| 国产一区二区三区视频了| 又大又爽又粗| 国产乱人伦免费视频| 亚洲七黄色美女视频| av电影中文网址| 露出奶头的视频| 久久 成人 亚洲| 香蕉av资源在线| 国产亚洲欧美精品永久| 国产不卡一卡二| 国产成人欧美在线观看| 久久国产乱子伦精品免费另类| bbb黄色大片| 午夜视频精品福利| 18禁黄网站禁片午夜丰满| 亚洲国产欧美日韩在线播放| 一级作爱视频免费观看| 18禁观看日本| 人人妻人人澡人人看| 欧美一区二区精品小视频在线| 久久性视频一级片| 久久精品91蜜桃| 18禁裸乳无遮挡免费网站照片 | 日本免费a在线| 亚洲中文日韩欧美视频| 午夜福利免费观看在线| 在线观看舔阴道视频| 露出奶头的视频| 成人精品一区二区免费| 黑丝袜美女国产一区| 国产成人精品久久二区二区免费| 久久久久久大精品| 久久久久久久久中文| 啦啦啦 在线观看视频| 国内精品久久久久精免费| 岛国视频午夜一区免费看| 久久国产精品人妻蜜桃| 亚洲熟妇中文字幕五十中出| 日日夜夜操网爽| 亚洲熟妇熟女久久| 中文字幕高清在线视频| 欧美激情极品国产一区二区三区| 一区二区三区国产精品乱码| 女性生殖器流出的白浆| 91九色精品人成在线观看| 国产精品自产拍在线观看55亚洲| 一本大道久久a久久精品| 听说在线观看完整版免费高清| 亚洲av成人一区二区三| 亚洲av日韩精品久久久久久密| 熟女少妇亚洲综合色aaa.| 亚洲精品在线观看二区| 欧美中文日本在线观看视频| 国产爱豆传媒在线观看 | 国产激情久久老熟女| 久久狼人影院| 久久久久精品国产欧美久久久| 丰满的人妻完整版| 亚洲人成网站高清观看| 黄色视频不卡| 亚洲午夜理论影院| 性色av乱码一区二区三区2| 亚洲va日本ⅴa欧美va伊人久久| 亚洲av第一区精品v没综合| 少妇熟女aⅴ在线视频| 欧美成人性av电影在线观看| 午夜激情av网站| 国产精品久久视频播放| 精品国产亚洲在线| 麻豆成人午夜福利视频| 国产精品久久久久久亚洲av鲁大| 久久久国产成人精品二区| 欧美av亚洲av综合av国产av| 亚洲国产中文字幕在线视频| 久久中文看片网| 50天的宝宝边吃奶边哭怎么回事| 亚洲精品国产精品久久久不卡| 真人做人爱边吃奶动态| av在线天堂中文字幕| 中文字幕精品免费在线观看视频| 亚洲黑人精品在线| 成人特级黄色片久久久久久久| 青草久久国产| 亚洲 国产 在线| 日韩一卡2卡3卡4卡2021年| 91成人精品电影| 久久草成人影院| 国产又色又爽无遮挡免费看| 亚洲天堂国产精品一区在线| 99在线人妻在线中文字幕| 免费一级毛片在线播放高清视频| 99国产综合亚洲精品| 老司机深夜福利视频在线观看| 满18在线观看网站| 国产视频一区二区在线看| 午夜免费激情av| 久久婷婷成人综合色麻豆| 可以在线观看的亚洲视频| 亚洲精品在线观看二区| 色综合欧美亚洲国产小说| av欧美777| 亚洲avbb在线观看| 亚洲美女黄片视频| 欧美日韩福利视频一区二区| 美国免费a级毛片| 无限看片的www在线观看| 亚洲精品美女久久av网站| 最新在线观看一区二区三区| 最好的美女福利视频网| 国产精品爽爽va在线观看网站 | 黑人巨大精品欧美一区二区mp4| 亚洲性夜色夜夜综合| 国产国语露脸激情在线看| 成人永久免费在线观看视频| 欧美激情 高清一区二区三区| 精品午夜福利视频在线观看一区| 99久久久亚洲精品蜜臀av| 在线观看午夜福利视频| 少妇 在线观看| 久久精品91无色码中文字幕| 国产黄片美女视频| 亚洲欧美一区二区三区黑人| 久久精品亚洲精品国产色婷小说| 亚洲九九香蕉| 中文字幕高清在线视频| 欧美激情极品国产一区二区三区| 亚洲 欧美一区二区三区| 国产欧美日韩精品亚洲av| 九色国产91popny在线| 久久99热这里只有精品18| 亚洲国产中文字幕在线视频| 欧美在线黄色| 国产精品亚洲av一区麻豆| 国产精品影院久久| 久99久视频精品免费| 国产一区二区三区视频了| 18禁裸乳无遮挡免费网站照片 | av天堂在线播放| 国产麻豆成人av免费视频| 久久婷婷成人综合色麻豆| 黄片大片在线免费观看| 一级a爱视频在线免费观看| cao死你这个sao货| 欧美国产精品va在线观看不卡| 一个人免费在线观看的高清视频| 成人一区二区视频在线观看| 在线免费观看的www视频| 在线播放国产精品三级| 男人舔女人下体高潮全视频| 国产成人欧美在线观看| 婷婷精品国产亚洲av在线| 一进一出抽搐动态| 99久久久亚洲精品蜜臀av| 成人手机av| 亚洲熟女毛片儿| 免费看十八禁软件| 黄色成人免费大全| 国产精品久久视频播放| 久久香蕉精品热| 香蕉国产在线看| 亚洲精品国产区一区二| 午夜免费激情av| 欧美成人午夜精品| 色精品久久人妻99蜜桃| 老司机在亚洲福利影院| 国产成人精品无人区| 在线十欧美十亚洲十日本专区| 免费观看人在逋| 一本综合久久免费| 国产精品精品国产色婷婷| 99久久99久久久精品蜜桃| 久久久久亚洲av毛片大全| 国产真实乱freesex| 亚洲精品中文字幕在线视频| 国产亚洲av嫩草精品影院| 欧美成人午夜精品| 国产成人影院久久av| 欧美三级亚洲精品| 色精品久久人妻99蜜桃| 一边摸一边做爽爽视频免费| 欧美成人一区二区免费高清观看 | 欧美丝袜亚洲另类 | 久久九九热精品免费| av电影中文网址| 丝袜美腿诱惑在线| 国产亚洲精品一区二区www| 久久久精品欧美日韩精品| 日日爽夜夜爽网站| 国产午夜精品久久久久久| 成人一区二区视频在线观看| 国产av不卡久久| 精品国内亚洲2022精品成人| 黄色视频不卡| 日日摸夜夜添夜夜添小说| 国产亚洲欧美在线一区二区| 国产1区2区3区精品| 久久久久久国产a免费观看| 一区二区三区国产精品乱码| 欧美性猛交黑人性爽| 亚洲av片天天在线观看| 视频区欧美日本亚洲| 波多野结衣巨乳人妻| 午夜福利免费观看在线| 亚洲av成人一区二区三| 俄罗斯特黄特色一大片| av电影中文网址| 日日干狠狠操夜夜爽| 老司机在亚洲福利影院| 最近最新免费中文字幕在线| 在线观看66精品国产| 欧美乱妇无乱码| 天天躁狠狠躁夜夜躁狠狠躁| 我的亚洲天堂| 日本a在线网址| 亚洲精品中文字幕在线视频| 亚洲欧美精品综合久久99| 波多野结衣高清作品| 亚洲va日本ⅴa欧美va伊人久久| 午夜精品在线福利| 18禁黄网站禁片免费观看直播| 亚洲天堂国产精品一区在线| 亚洲九九香蕉| 国产99白浆流出| 亚洲成国产人片在线观看| 国产精品,欧美在线| 亚洲av成人一区二区三| 十八禁人妻一区二区| 国产三级在线视频| 人人澡人人妻人| 一夜夜www| 亚洲精品中文字幕在线视频| 人妻久久中文字幕网| 国产欧美日韩精品亚洲av| 久久久久国产精品人妻aⅴ院| 中文字幕人妻熟女乱码| 草草在线视频免费看| 色播在线永久视频| 男人的好看免费观看在线视频 | 精品免费久久久久久久清纯| 动漫黄色视频在线观看| 婷婷丁香在线五月| 国产三级黄色录像| 在线观看66精品国产| 国产视频内射| 欧美黑人欧美精品刺激| 不卡av一区二区三区| 99精品久久久久人妻精品| 久热这里只有精品99| 99久久99久久久精品蜜桃| 桃红色精品国产亚洲av| 美女免费视频网站| 日日摸夜夜添夜夜添小说| 亚洲精华国产精华精| 人人妻,人人澡人人爽秒播| 最近最新中文字幕大全电影3 | 人人妻人人看人人澡| 这个男人来自地球电影免费观看| 色在线成人网| 欧美最黄视频在线播放免费| 欧美黑人巨大hd| av中文乱码字幕在线| 91老司机精品| 亚洲黑人精品在线| 午夜免费成人在线视频| 精品一区二区三区av网在线观看| 两个人视频免费观看高清| 男人舔女人的私密视频| 级片在线观看| 日韩欧美免费精品| 少妇粗大呻吟视频| 黄色丝袜av网址大全| 制服丝袜大香蕉在线| 99久久久亚洲精品蜜臀av| 欧美国产日韩亚洲一区| 免费看美女性在线毛片视频| 每晚都被弄得嗷嗷叫到高潮| 国产精品久久电影中文字幕| 国产精品自产拍在线观看55亚洲| 欧美日韩黄片免| 熟女少妇亚洲综合色aaa.| 久久久精品欧美日韩精品| 999久久久精品免费观看国产| 亚洲avbb在线观看| 黄色视频,在线免费观看| 国产成人欧美| 黄色丝袜av网址大全| 黄频高清免费视频| 亚洲 国产 在线| 久久热在线av| 色播亚洲综合网| 2021天堂中文幕一二区在线观 | 亚洲精品中文字幕一二三四区| 日韩免费av在线播放| 欧美大码av| 国产成人精品无人区| av天堂在线播放| 脱女人内裤的视频| 丝袜人妻中文字幕| 美女午夜性视频免费| 国产亚洲精品综合一区在线观看 |