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

    北極高緯度苔原CO2凈通量的空間變異規(guī)律及影響因素

    2021-03-05 02:08:08陳清清包韜朱仁斌徐華
    極地研究 2021年4期
    關(guān)鍵詞:苔原海鳥(niǎo)北極

    陳清清包韜朱仁斌 徐華

    研究論文

    北極高緯度苔原CO2凈通量的空間變異規(guī)律及影響因素

    陳清清1,2包韜3朱仁斌2徐華4

    (1宜春學(xué)院生命科學(xué)與資源環(huán)境學(xué)院, 江西 宜春 336000;2中國(guó)科學(xué)技術(shù)大學(xué)地球和空間科學(xué)學(xué)院極地環(huán)境與全球變化安徽省重點(diǎn)實(shí)驗(yàn)室, 安徽 合肥 230026;3中國(guó)科學(xué)院東亞區(qū)域氣候-環(huán)境重點(diǎn)實(shí)驗(yàn)室, 中國(guó)科學(xué)院大氣物理研究所, 北京 100029;4中國(guó)科學(xué)院南京土壤研究所土壤與農(nóng)業(yè)可持續(xù)發(fā)展國(guó)家重點(diǎn)實(shí)驗(yàn)室, 江蘇 南京 210091)

    北極夏季無(wú)冰區(qū)苔原CO2的釋放對(duì)全球碳循環(huán)起著重要的作用。在第4次北極科學(xué)考察(2008年7月26日—8月5日)期間, 采用密閉箱法在北極新奧爾松(Ny-?lesund)黃河站區(qū)附近的鳥(niǎo)類保護(hù)區(qū)苔原(Tundra in Seabird Sanctuary, TSB)、普通海灘苔原(Tundra in Non-seabird Colony, TNS)和兩者間的過(guò)渡苔原(Tundra in Transition Zone, TTR)監(jiān)測(cè)CO2凈通量(Net CO2Fluxes, NEE)空間變異規(guī)律及其影響因素。結(jié)果表明: 鳥(niǎo)類保護(hù)區(qū)苔原TSB的平均 NEE為(–39.0±6.0)mg·m–2·h–1, 為CO2吸收匯, 且對(duì)CO2的吸收隨著植被覆蓋度和海鳥(niǎo)活動(dòng)強(qiáng)度的增強(qiáng)而顯著增強(qiáng); 普通海灘苔原TNS和過(guò)渡苔原TTR的平均NEE分別為 (12.0± 6.3)mg·m–2·h–1和(40.5±29.3)mg·m–2·h–1, 均為CO2凈排放源。普通海灘苔原TNS的CO2排放強(qiáng)度隨土壤水分的減少和地勢(shì)增高而增強(qiáng); 過(guò)渡苔原TTR中高地苔原區(qū)的平均NEE為(106.4±23.1)mg·m–2·h–1, 為CO2強(qiáng)排放源; 而泥炭苔原區(qū)的平均NEE為(–58.3±9.5)mg·m–2·h–1, 為CO2強(qiáng)吸收匯。苔原土壤的理化性質(zhì)對(duì)苔原NEE有重要影響, 鳥(niǎo)類保護(hù)區(qū)苔原TSB和普通海灘苔原TNS的NEE與土壤含水率呈顯著負(fù)相關(guān)關(guān)系(= –0.44,=0.003)。鳥(niǎo)類保護(hù)區(qū)苔原TSB的NEE與土壤溫度弱正相關(guān)(=0.32,=0.06), 與NH4+-N(<0.05)和NO3–-N(<0.05)含量均顯著負(fù)相關(guān)。在全球變暖的驅(qū)使下, 不同地形地貌和海鳥(niǎo)活動(dòng)特征的北極苔原區(qū)域CO2排放的復(fù)雜性將顯著增強(qiáng), 對(duì)全球碳循環(huán)產(chǎn)生較大的影響, 本研究也將為此提供科研資料。

    CO2凈通量北極新奧爾松高緯度苔原

    0 引言

    北極和亞北極生態(tài)系統(tǒng)的有機(jī)碳主要儲(chǔ)存在上層凍土中, 且儲(chǔ)量巨大, 約占全球土壤有機(jī)碳的50%[1-3]。氣候變暖導(dǎo)致凍土融化, 表層土壤微生物活動(dòng)增強(qiáng), 北極和亞北極苔原生態(tài)系統(tǒng)具有強(qiáng)大的CO2排放潛力[4], 但從已有研究來(lái)看, 北極和亞北極是CO2的排放源還是吸收匯尚不明確。亞北極的俄羅斯?jié)駶?rùn)苔原區(qū)域的CO2凈通量(NEE)波動(dòng)范圍為–79~–41g·m–2·a–1, 為CO2吸收匯[5]。北極泥炭地和沼澤的年平均CO2凈通量(NEE)波動(dòng)范圍為–170~–458 mg·m–2·h–1, 平均值為–218 kg·ha–1·a–1, 為CO2強(qiáng)吸收匯[6]。而干旱苔原區(qū)域的NEE監(jiān)測(cè)結(jié)果為275 kg·ha–1·a–1, 為CO2凈排放源[7]。除南極外, 北極高緯度苔原是海鳥(niǎo)主要棲息地, 格陵蘭和斯瓦爾巴群島聚居大量的海鳥(niǎo)[8-9]。海鳥(niǎo)活動(dòng)是極地苔原外源碳(C)、氮(N)、磷(P)、硫(S)等營(yíng)養(yǎng)物質(zhì)輸入的主要途徑, 并能夠在長(zhǎng)時(shí)間尺度下改變影響苔原土壤和植被的發(fā)育[10], 促進(jìn)海鳥(niǎo)聚居區(qū)及附屬區(qū)域生態(tài)系統(tǒng)的發(fā)育[11-12], 最終影響溫室氣體的排放。海洋動(dòng)物聚居區(qū)的溫室氣體通量比普通苔原區(qū)的顯著增高, 尤其是南極企鵝聚居區(qū)和海豹聚居區(qū)[13-16], 對(duì)北極高緯度苔原區(qū)域的CO2排放及其影響的研究較少。Wüthrich等[17]和Lloyd[18]研究發(fā)現(xiàn)艾登布克塔(Eidembukta)和新奧爾松(Ny-?lesund)苔原在1995年為CO2凈吸收匯, 而在1996年轉(zhuǎn)變?yōu)閮襞欧旁础?/p>

    苔原CO2的排放取決于光合作用和呼吸作用強(qiáng)度的差異, 是由苔原植被變化特征和土壤理化性質(zhì)等決定的[19-20], 受到土壤溫度、土壤水分、植被群落和動(dòng)物干擾的影響[19-21]。海洋動(dòng)物活動(dòng)強(qiáng)度的區(qū)域差異性也會(huì)造成極地苔原系統(tǒng)土壤和植被的發(fā)育的差異; 對(duì)于大部分不受海洋動(dòng)物糞便影響的區(qū)域, 地形地貌造成的土壤和植被發(fā)育也是CO2排放的重要因素。全球變暖的趨勢(shì)下, 海鳥(niǎo)活動(dòng)增強(qiáng)、夏季無(wú)冰區(qū)的擴(kuò)張和生長(zhǎng)季延長(zhǎng)都能增加苔原匯的強(qiáng)度[22-24], 而鳥(niǎo)糞營(yíng)養(yǎng)增加, 土壤呼吸作用增強(qiáng)等則會(huì)增加苔原CO2的排放[25-26], 因此, 影響北極高緯度苔原生態(tài)系統(tǒng)CO2排放的因素更多也更復(fù)雜, 但尚未開(kāi)展系統(tǒng)科學(xué)的研究。

    本文在第4次北極科學(xué)考察(2008年7月26日—8月5日)期間, 選取了北極高緯度的斯瓦爾巴德群島新奧爾松(Ny-?lesund)的鳥(niǎo)類保護(hù)區(qū)苔原(TSB)、普通海灘苔原(TNS)和兩者之間的過(guò)渡苔原(TTR)為研究區(qū)域, 觀測(cè)北極高緯度苔原不同空間尺度CO2凈通量(NEE)變化規(guī)律并分析其影響因素。本文的研究目標(biāo)是: (1)鳥(niǎo)類保護(hù)區(qū)苔原TSB、普通海灘苔原TNS和過(guò)渡苔原TTR的CO2凈通量(NEE)的空間變化規(guī)律; (2)苔原土壤理化性質(zhì)的差異; (3)苔原土壤CO2凈通量(NEE)的影響因素。

    1 研究區(qū)域與研究方法

    1.1 研究區(qū)域概況

    本文的研究區(qū)域?yàn)楸睒O挪威群島西部的斯瓦爾巴德半島的新奧爾松(Ny-?lesund)的苔原區(qū)域, (圖1a)。北大西洋的高爾夫暖流流向斯瓦爾巴德群島西部, 其支流西斯匹次卑爾根環(huán)流流向北極最北部的開(kāi)闊海域[24,27], 與氣旋共同作用下, 斯瓦爾巴德半島西部呈海洋性氣候, 夏天較涼爽, 冬天較溫和, 年平均氣溫–5.8 ℃, 比同緯度的格陵蘭高6~7 ℃[28]。全年以東南風(fēng)為主導(dǎo)風(fēng)向, 1~2 m·s–1風(fēng)速的出現(xiàn)頻率最高[29]。年平均降水量約400 mm, 且往內(nèi)陸逐漸減少, 新奧爾松降水量為每年200~300 mm。植物種類較少, 以北部極地苔原植物區(qū)系和極地沙漠植物區(qū)系為主, 包括、、、和等[30]。野生動(dòng)物主要包括北極狐、北極馴鹿、北極鵝及30余種海鳥(niǎo)[31-32]。新奧爾松典型的無(wú)冰期是從6月持續(xù)到9月[33]。6月初凍土開(kāi)始融化, 9月初最大融化深度為1.6~2.0 m, 活躍層厚度也達(dá)到最大, 并持續(xù)到9月中旬[34]。夏季有80%的時(shí)間為多云天氣, 導(dǎo)致太陽(yáng)輻照強(qiáng)度較低(http://www.yr.no/place/ Norway/Svalbard/)。新奧爾松區(qū)域整體為淤泥黏性土壤和沙質(zhì)發(fā)育土壤覆蓋, 有機(jī)質(zhì)含量較低。采樣期間平均溫度為(2.1±0.9) ℃, 平均降水量大約是400 mm。

    本文在第4次北極科學(xué)考察期間, 在3個(gè)苔原區(qū)域監(jiān)測(cè)CO2凈通量(NEE)的空間變異規(guī)律(圖1b)。

    區(qū)域Ⅰ為鳥(niǎo)類保護(hù)區(qū)苔原(TSB), 位于新奧爾松(Ny-?lesund)的西北部。海鳥(niǎo)活動(dòng)的重點(diǎn)區(qū)域面積為0.5 km2, 海拔為5~20 m。鳥(niǎo)類保護(hù)區(qū)苔原的海鳥(niǎo)種類較多, 主要品種有: razor-billed auk ()Arctic duck ()kittiwake()Arctic Tern()fulmar()和 albatross()。苔原植被覆蓋區(qū)域面積>98%, 裸露的土壤和石頭區(qū)域面積約占2 %。主要植被類型是苔蘚(和)、地衣()和莎草(和)。土壤有機(jī)層厚度約為15~20 cm。

    區(qū)域Ⅱ?yàn)槠胀ê┨υ?TNS), 位于新奧爾松(Ny-?lesund)的東南。此苔原基本沒(méi)有海鳥(niǎo)活動(dòng), 但有北極兔和馴鹿活動(dòng)蹤跡。普通海灘苔原的 48%為苔蘚或者地衣覆蓋, 20 %生長(zhǎng)著維管植物, 28 %是裸露的土壤, 4 %是石頭。土壤有機(jī)層厚約為10~20 cm。

    圖1 北極斯瓦爾巴(Svalbard)群島、新奧爾松(Ny-?lesund)地區(qū)的地理位置、研究區(qū)域分布圖和鳥(niǎo)類保護(hù)區(qū)采樣點(diǎn)分布圖.a)北極斯瓦爾巴(Svalbard)群島和新奧爾松(Ny-?lesund)地區(qū)的地理位置; b)研究區(qū)域分布圖(Ⅰ為鳥(niǎo)類保護(hù)區(qū)苔原TSB, Ⅱ?yàn)檎緟^(qū)附近的普通海灘苔原TNS, Ⅲ為過(guò)渡苔原TTR); c)鳥(niǎo)類保護(hù)區(qū)采樣點(diǎn)分布圖

    Fig.1.The location of Ny-?lesund, Svalbard in Arctic and the distribution map for the study during the net CO2fluxes observation period.a) the location of Ny-?lesund, Svalbard; b) the study area (Ⅰ-TSB, Ⅱ-TNS, Ⅲ-TTR); c) observation sites in the bird sanctuary

    區(qū)域Ⅲ為過(guò)渡苔原區(qū)域(TTR), 位于黃河站區(qū)鳥(niǎo)類保護(hù)區(qū)苔原(TSB)和普通海灘苔原(TNS)之間的過(guò)渡區(qū)域, 離黃河站很近, 土壤情況較復(fù)雜, 地表植被稀疏。

    1.2 野外CO2凈通量(NEE)的觀測(cè)點(diǎn)實(shí)驗(yàn)設(shè)計(jì)

    1.2.1 鳥(niǎo)類保護(hù)區(qū)苔原(TSB)

    在鳥(niǎo)類保護(hù)區(qū)苔原(TSB)采用靜態(tài)箱法監(jiān)測(cè)CO2凈通量(NEE)的空間規(guī)律。根據(jù)植被覆蓋度、水文條件以及海鳥(niǎo)活動(dòng)等的差異將鳥(niǎo)類保護(hù)區(qū)苔原TSB分成3類觀測(cè)區(qū): 植被覆蓋度最高、強(qiáng)烈的海鳥(niǎo)活動(dòng)區(qū)HB(High Seabird-activity Area)觀測(cè)點(diǎn)HB1~HB9, 植被覆蓋度較高、中等強(qiáng)度的海鳥(niǎo)活動(dòng)區(qū)MB (Medium Seabird-activity Area)觀測(cè)點(diǎn)MB1~MB6, 以及植被覆蓋度低、海鳥(niǎo)活動(dòng)邊緣區(qū)LB (Low Seabird-activity Area)測(cè)點(diǎn)LB1~LB4。共隨機(jī)設(shè)置19組采樣點(diǎn), 探討苔原TSB的CO2凈通量(NEE)的空間變化。這19組采樣點(diǎn)覆蓋了苔原TSB的典型海鳥(niǎo)活動(dòng)和植被變化區(qū), 且每組內(nèi)設(shè)置兩個(gè)采樣點(diǎn)互為對(duì)照。

    為使鳥(niǎo)類保護(hù)區(qū)苔原的CO2凈通量(NEE)具有可對(duì)比性, 綜合考慮植被情況、水文條件和海鳥(niǎo)活動(dòng)情況等, 本研究也選取了普通海灘苔原(TNS)和過(guò)渡苔原(TTR)為研究區(qū)域并進(jìn)行NEE空間變化規(guī)律觀測(cè)。

    1.2.2 非海鳥(niǎo)聚居地的普通海灘苔原(TNS)

    依據(jù)地形地貌特征, 在普通海灘苔原(TNS)區(qū)域從地勢(shì)低的濕潤(rùn)苔原區(qū)到地勢(shì)高的干旱苔原區(qū)隨機(jī)設(shè)置5大類43組觀測(cè)點(diǎn), 采用靜態(tài)箱法, 觀測(cè)苔原TNS的CO2凈通量(NEE)的空間變化規(guī)律。5大類觀測(cè)點(diǎn)的設(shè)置為: 地勢(shì)低的濕地苔原觀測(cè)區(qū)Wet Sites(SW)的8組觀測(cè)點(diǎn)(SW1~SW8)、地勢(shì)略高的半干旱苔原觀測(cè)區(qū)Mesic Sites(SM)的8組觀測(cè)點(diǎn)(SM1~SM8)、地勢(shì)更高的干旱斜坡苔原觀測(cè)區(qū)Dry Slope Sites(SD)的8組觀測(cè)點(diǎn)(SD1~ SD8)、地勢(shì)高的高地苔原觀測(cè)區(qū)Upland Sites(UD)的9組觀測(cè)點(diǎn)(UD1~UD9)和地勢(shì)最高的頂部苔原觀測(cè)區(qū)Higher Upland Sites(UL)的10組觀測(cè)點(diǎn)(UL1~UL10), 囊括了普通海灘苔原TNS的代表性區(qū)域, 每組觀測(cè)點(diǎn)設(shè)置兩個(gè)重復(fù)。

    1.2.3 過(guò)渡苔原(TTR)

    過(guò)渡苔原(TTR)在黃河站附近, 充分考慮了苔原地勢(shì)等表觀特征, 設(shè)置高地苔原HT觀測(cè)區(qū)3組采樣點(diǎn)(HT1~HT3)和泥炭地苔原NM觀測(cè)區(qū)2組采樣點(diǎn)(NM1和NM2), 采用靜態(tài)箱法, 觀測(cè)苔原TTR的CO2凈通量(NEE)的空間變化規(guī)律。

    1.3 CO2凈通量(NEE)的測(cè)定

    本文采用靜態(tài)箱法監(jiān)測(cè)北極高緯度苔原生態(tài)系統(tǒng)的CO2凈通量(NEE)[35]。CO2采樣系統(tǒng)即靜態(tài)箱系統(tǒng)由一個(gè)一端開(kāi)口的丙烯酸樹(shù)脂的透明箱體和一個(gè)底座構(gòu)成。本次采樣共用兩組規(guī)格的采樣箱, 一組箱體外徑的長(zhǎng)×寬×高是50cm×50cm× 50 cm, 另一組箱體外徑的長(zhǎng)×寬×高為50cm× 50cm×30 cm。為取樣后不改變箱內(nèi)氣壓, 不影響氣體擴(kuò)散, 箱子露出土壤的高度至少為25 cm[36]。將PVC底座插入土壤約5 cm, 放置透明箱子后, 向底座注水密封, 密封面積0.25 m2。采集CO2氣體, 將兩通針的一端插入透明箱上方, 另一端連接帶橡膠塞的真空瓶, 每隔15 min采一次氣體樣品, 共三次。測(cè)定順序不會(huì)影響采樣結(jié)果, 同步記錄箱內(nèi)溫度和采樣時(shí)的氣象條件。真空瓶體積為17.5 mL, 采樣前抽成高真空。在第4次北極科學(xué)考察期間, 采用靜態(tài)箱法, 兩組不同尺寸采樣箱進(jìn)行流動(dòng)采樣[37], 觀測(cè)鳥(niǎo)類保護(hù)區(qū)苔原(TSB)、普通海灘苔原(TNS)和過(guò)渡苔原(TTR)的CO2凈通量(NEE)的空間變化規(guī)律。

    1.4 CO2濃度測(cè)定和凈通量計(jì)算

    氣體樣品CO2濃度在中國(guó)科學(xué)院南京土壤研究所土壤與農(nóng)業(yè)可持續(xù)發(fā)展國(guó)家重點(diǎn)實(shí)驗(yàn)室用氣相色譜分析儀(Shimadzu GC-14B)檢測(cè)。檢測(cè)條件: 檢測(cè)器為熱導(dǎo)檢測(cè)器(TCD), 色譜柱為80/100目Chromosorb 102填充柱, 進(jìn)樣器、檢測(cè)器和填充柱的溫度分別為100 ℃、60 ℃和60 ℃, 高純H2作載氣, 載氣流速為80 mL·min–1[38]。標(biāo)準(zhǔn)氣體的濃度為641 ppm, 由南京特種氣體廠生產(chǎn), 并經(jīng)日本國(guó)立農(nóng)業(yè)食品產(chǎn)業(yè)技術(shù)綜合研究機(jī)構(gòu)農(nóng)業(yè)環(huán)境技術(shù)研究所(原國(guó)立農(nóng)業(yè)環(huán)境技術(shù)研究所)校正, 平均值為(640.9±1.43) ppm (=15)。

    CO2凈通量(NEE)的計(jì)算公式為=ρ×× (d/ d)×273/(273+)。式中:為CO2凈通量(NEE), 單位時(shí)間單位面積交換CO2的質(zhì)量, 單位為mg·m–2·h–1; ρ為標(biāo)準(zhǔn)狀態(tài)下CO2的密度1.964 kg·m–3;為箱子有效高度(高出地表的高度), 單位為m; d/d為濃度隨時(shí)間的變化率, 單位為ppm·h–1; t為箱內(nèi)氣溫, 單位為℃。文中CO2凈通量(NEE)的正值表示排放CO2即為CO2的排放源, 負(fù)值表示苔原吸收CO2即為CO2的吸收匯。

    1.5 氣象資料、土壤采樣及其理化性質(zhì)檢測(cè)

    采樣氣體樣品時(shí), 同時(shí)記錄新奧爾松地區(qū)的氣象數(shù)據(jù)。在CO2凈通量觀測(cè)點(diǎn)插入溫度計(jì)到地表以及地表以下5、10和15 cm, 檢測(cè)0 cm、–5 cm、–10 cm和–15 cm的土壤溫度。在對(duì)應(yīng)氣體樣品點(diǎn)旁邊采集表層0~20 cm處的土壤作為土壤樣品, 冷凍保存至實(shí)驗(yàn)室進(jìn)行檢測(cè)。土壤含水率(Wc)是將新鮮土壤樣品在105 ℃條件下烘至恒重, 根據(jù)公式Wc=(水分重量/干土重) ×100%計(jì)算。pH測(cè)定方法: 用天平稱取4 g室溫下風(fēng)干的土壤樣品, 放入小燒杯中, 按照質(zhì)量體積比 1∶2.5, 量取10 mL超純水倒入小燒杯, 加磁力振子在磁力攪拌器上高速旋轉(zhuǎn)2 min后, 靜置至上層液澄清, 測(cè)量上清液pH??傆袡C(jī)碳(TOC)用重鉻酸鉀體積法測(cè)定, 測(cè)量方法的分析誤差為2.5%。土壤NH4+-N 和 NO3–-N含量分別用靛酚藍(lán)比色法和格里斯分光光度法測(cè)定[39]。土壤總磷(TP)含量采用鉬酸銨分光光度法檢測(cè)。土壤C、N、S、H含量是將土壤在室溫下晾干, 然后用德國(guó)Vario·EL元素分析儀(Elementar Vario EL, Hanau, Germany)測(cè)定, 相對(duì)誤差小于1 %[30]。

    1.6 統(tǒng)計(jì)分析

    實(shí)驗(yàn)室內(nèi)主要是利用微軟Excel 2016、 SPSS以及OriginPro7.5計(jì)算、分析數(shù)據(jù), 繪制圖和表。微軟Excel 2016計(jì)算CO2凈通量(NEE), 使用標(biāo)準(zhǔn)誤差()來(lái)衡量NEE的不確定性。使用OriginPro7.5繪制NEE空間變化圖。用OriginPro7.5和SPSS分析不同觀測(cè)點(diǎn)和不同苔原NEE的顯著性差異, 并分析NEE與土壤化學(xué)性質(zhì)(土壤溫度、土壤含水率、土壤碳氮比(C/N), 及土壤中的TC、TN、NH4+-N、NO3–-N和TP含量)的相關(guān)性及繪圖。

    2 結(jié)果與討論

    2.1 苔原土壤理化性質(zhì)

    苔原TSB和苔原TNS的土壤理化性質(zhì)如表1所示。植被覆蓋度高、海鳥(niǎo)活動(dòng)強(qiáng)烈區(qū)HB的平均TC和TN含量(9.1 %和1.1 %)最高。P和S是海鳥(niǎo)糞便的典型的基本元素[40], 海鳥(niǎo)會(huì)增加苔原土壤C、N、P和S等營(yíng)養(yǎng)物質(zhì)的輸入, 會(huì)顯著地增加苔原土壤TN、NH4+-N和NO3–-N含量[14,41]。因此, 在鳥(niǎo)類保護(hù)區(qū)苔原(TSB)和普通海灘苔原(TNS)的所有觀測(cè)區(qū)中, 植被覆蓋度高、海鳥(niǎo)活動(dòng)強(qiáng)烈的觀測(cè)區(qū)HB土壤的TP和TS含量極高, 比植被覆蓋和海鳥(niǎo)活動(dòng)中等強(qiáng)度觀測(cè)區(qū)MB及植被覆蓋度低和海鳥(niǎo)活動(dòng)邊緣區(qū)觀測(cè)區(qū)LB以及苔原TNS土壤高一個(gè)數(shù)量級(jí); NH4+-N和NO3–-N的平均含量比觀測(cè)區(qū)MB、LB和苔原TNS土壤高1到2個(gè)數(shù)量級(jí); 而C/N值則低于這些區(qū)域。觀測(cè)區(qū)HB和MB土壤的pH值略微高于觀測(cè)區(qū)LB。苔原TSB的pH值與苔原TNS的相近, 苔原TNS土壤的pH從地勢(shì)低的濕潤(rùn)苔原到地勢(shì)高的高地干旱苔原呈減小趨勢(shì)。由于海鳥(niǎo)糞便增加了苔原土壤有機(jī)物的輸入以及海鳥(niǎo)踩踏等海鳥(niǎo)活動(dòng)的干擾, 所以苔原TSB土壤的TN、NH4+-N、NO3–-N、TC、TP、TS和土壤含水率(Wc)高于苔原TNS土壤, 其中NH4+-N含量約為苔原TNS土壤的3倍, NO3–-N含量比苔原TNS土壤高約2個(gè)數(shù)量級(jí), TS含量比苔原TNS土壤高約1個(gè)數(shù)量級(jí)。

    表1 新奧爾松苔原土壤的理化性質(zhì)

    續(xù)表

    2.2 鳥(niǎo)類保護(hù)區(qū)苔原(TSB)CO2凈通量(NEE)的空間變化

    鳥(niǎo)類保護(hù)區(qū)苔原(TSB)的CO2凈通量(NEE)的空間結(jié)果如表2和圖2a所示。在植被覆蓋度最高、海鳥(niǎo)活動(dòng)最強(qiáng)烈觀測(cè)區(qū)HB(HB1~HB9), 采樣點(diǎn)的所有NEE數(shù)據(jù)中約有35 %均小于零, 變化范圍為–188.0~84.5 mg·m–2·h–1, 平均值為(–76.4±10.6)mg·m–2·h–1, 強(qiáng)烈吸收CO2。在植被覆蓋度較高、海鳥(niǎo)活動(dòng)強(qiáng)度中等的觀測(cè)區(qū)MB(MB1~ MB6), 超過(guò)65 %的NEE數(shù)據(jù)為負(fù)值, 變化范圍為–145.3~86.9 mg·m–2·h–1, 均值為(–29.9±8.3) mg·m–2·h–1, 為CO2吸收匯, 但吸收強(qiáng)度顯著(≤0.05)低于觀測(cè)區(qū)HB。在植被覆蓋度低、海鳥(niǎo)活動(dòng)邊緣觀測(cè)區(qū)LB(LB1~LB4), 超過(guò)80 %的NEE數(shù)據(jù)為負(fù)值, 變化范圍–98.7~175.9 mg·m–2·h–1, 平均NEE為(–1.7±7.4)mg·m–2·h–1, 弱吸收CO2, 吸收強(qiáng)度顯著(≤0.05)低于觀測(cè)區(qū)HB和MB。

    表2 鳥(niǎo)類保護(hù)區(qū)苔原(TSB)、普通海灘苔原(TNS)和過(guò)渡苔原(TTR)CO2凈通量(NEE)在夏季(2008年7月26—8月5日)的空間變化

    注: 不同字母a, b, c表示具有顯著性差異(≤0.05); 不同字母A, B,C表示在≤0.1水平上存在差異。

    圖2 新奧爾松地區(qū)研究區(qū)域CO2凈通量(NEE)的空間變化.Ⅰ為鳥(niǎo)類保護(hù)區(qū)苔原(TSB), Ⅱ?yàn)檎緟^(qū)附近的普通海灘苔原(TNS), Ⅲ為過(guò)渡苔原(TTR)

    Fig.2.Spatial variation of net CO2fluxes (NEE) in the study area of Ny-?lesund.Ⅰ-TSB, Ⅱ-TNS, Ⅲ-TTR

    就苔原TSB整體而言, NEE的波動(dòng)范圍為–188.0~175.9 mg·m–2·h–1, 均值為(–39.0±6.0) mg·m–2·h–1, 為CO2強(qiáng)吸收匯。海鳥(niǎo)活動(dòng)也增加了苔原TSB的外源營(yíng)養(yǎng)物質(zhì)(C、N、P等)的輸入, 土壤營(yíng)養(yǎng)豐富, 促進(jìn)了苔蘚的發(fā)育, 植被覆蓋度大, CO2吸收強(qiáng)烈。同時(shí)由于土壤營(yíng)養(yǎng)物質(zhì)和海鳥(niǎo)活動(dòng)的差異, 苔原TSB的植被分布也呈現(xiàn)明顯的差異, 導(dǎo)致NEE的空間變異顯著(為173.0%)[35]。

    2.3 普通海灘苔原(TNS)CO2凈通量(NEE)的空間變化

    本研究設(shè)置了觀測(cè)點(diǎn)地勢(shì)低的濕地苔原觀測(cè)區(qū)SW對(duì)應(yīng)點(diǎn)位SW1~SW8、地勢(shì)略高的半干旱苔原觀測(cè)區(qū)SM對(duì)應(yīng)點(diǎn)位SM1~SM8地勢(shì)更高的干旱斜坡苔原觀測(cè)區(qū)SD對(duì)應(yīng)點(diǎn)位SD1~SD8、地勢(shì)高的高地苔原觀測(cè)區(qū)UD對(duì)應(yīng)點(diǎn)位UD1~UD9和地勢(shì)最高的頂部苔原觀測(cè)區(qū)UL對(duì)應(yīng)點(diǎn)位UL1~UL10, 觀測(cè)苔原TNS的CO2凈通量(NEE)空間變化情況, 結(jié)果如表2和圖2b所示。地勢(shì)低的濕地苔原觀測(cè)區(qū)SW的所有NEE數(shù)據(jù)中約有60 %小于零, 整體波動(dòng)范圍為–142.2~39.2 mg·m–2·h–1, 均值為(–23.8±13.4) mg·m–2·h–1, 為CO2吸收匯。地勢(shì)略高的半干旱苔原觀測(cè)區(qū)SM的所有NEE數(shù)據(jù)中小于零的超過(guò)75 %, 變化范圍在–113.2~ 117.9 mg·m–2·h–1, 平均值為(–15.8±13.9) mg·m–2·h–1,為CO2吸收匯, 但其吸收強(qiáng)度在0.1水平上弱于觀測(cè)區(qū)SW。地勢(shì)更高的干旱斜坡苔原觀測(cè)區(qū)SD的所有NEE數(shù)據(jù)中約有70 %小于零, 整體波動(dòng)范圍為–58.4~79.6 mg·m–2·h–1, 平均NEE為(–0.2±9.8) mg·m–2·h–1, 為CO2弱吸收匯, 吸收強(qiáng)度在0.1水平上弱于觀測(cè)區(qū)SW。地勢(shì)高的高地苔原觀測(cè)區(qū)UD的所有NEE數(shù)據(jù)中約有25 %小于零, 整體變化范圍為–64.2~161.8 mg·m–2·h–1, 平均NEE為(43.5±10.7) mg·m–2·h–1, 為CO2強(qiáng)排放源, 在0.1水平上與觀測(cè)區(qū)SW和SD具有差異性。地勢(shì)最高的頂部苔原觀測(cè)區(qū)UL所有NEE數(shù)據(jù)中小于零的少于25%, 變化范圍在–25.4~194.7 mgCO2·m–2·h–1之間, 平均NEE為44.2±12.4 mg·m–2·h–1, 為CO2強(qiáng)排放源, 在0.1水平上與觀測(cè)區(qū)SW和SD具有差異性, 在0.05水平上與觀測(cè)區(qū)UD具有顯著差異性。

    綜合來(lái)說(shuō), 苔原TNS的NEE為(12.0± 5.7) mg·m–2·h–1, 整體表現(xiàn)為CO2排放源;為438.1 %, 苔原TNS的CO2排放具有很大的空間變異性。苔原TNS的NEEs的空間監(jiān)測(cè)結(jié)果表明, 地勢(shì)低的濕地苔原觀測(cè)區(qū)SW和地勢(shì)略高半干旱苔原觀測(cè)區(qū)SM均為CO2強(qiáng)吸收匯, 地勢(shì)更高的干旱斜坡苔原觀測(cè)區(qū)SD為CO2弱吸收匯, 地勢(shì)高的高地苔原觀測(cè)區(qū)UD和地勢(shì)最高的頂部苔原觀測(cè)區(qū)UL為CO2強(qiáng)排放源(圖2b); 不同苔原觀測(cè)區(qū)CO2的排放隨土壤水量的增加而降低, 隨地勢(shì)的升高而升高。苔原TNS受海鳥(niǎo)活動(dòng)影響極為微小, 以砂石和早期土壤發(fā)育特征為主, 地形和地貌的差別導(dǎo)致了植被分布的差異, 最終導(dǎo)致了不同苔原觀測(cè)區(qū)NEE的空間差異[42]。

    2.4 過(guò)渡苔原(TTR)CO2凈通量(NEE)的空間變化

    在過(guò)渡苔原(TTR)設(shè)置的高地苔原觀測(cè)區(qū)HT(HT1~HT3)和泥炭地觀測(cè)區(qū)NM(NM1~ NM2), 其CO2凈通量(NEE)結(jié)果如表2和圖2c所示。高地苔原觀測(cè)區(qū)HT的所有NEE數(shù)據(jù)均高于30 mg·m–2·h–1, 變化范圍為36.7~165.4 mg·m–2·h–1之間, 平均值為(106.4±23.1) mg·m–2·h–1, 是苔原TSB、苔原TNS和苔原TTR所有觀測(cè)區(qū)中的最強(qiáng)CO2排放源。而泥炭地觀測(cè)區(qū)NM的所有NEE數(shù)據(jù)均低于–40 mg·m–2·h–1, 波動(dòng)范圍為–71.9~ –44.1 mg·m–2·h–1, 平均NEE為(–58.3±9.5) mg·m–2·h–1, 是苔原TSB、苔原TNS和苔原TTR所有觀測(cè)區(qū)中的強(qiáng)吸收匯。

    苔原TTR的NEE為(40.5±29.3) mg·m–2·h–1, 為CO2強(qiáng)排放源;為228.4 %, 苔原TTR的CO2排放具有很大的空間變異性。苔原TTR受海鳥(niǎo)活動(dòng)影響微弱, 地形、地貌特征和植被分布情況是影響苔原CO2排放的主要因素[42]。

    2.5 新奧爾松(Ny–?lesund)苔原CO2凈通量(NEE)的分析與討論

    總體而言, 苔原TSB的整體NEE為(–19.1±9.6) mg·m–2·h–1, 為CO2吸收匯; 苔原TNS的平均NEE為(12.0±5.7) mg·m–2·h–1, 為CO2排放源; 苔原TTR中高地苔原的平均NEE為(106.4±23.1) mg·m–2·h–1, 為CO2強(qiáng)排放源, 而泥炭地的平均NEE為 (–58.3± 9.5) mg·m–2·h–1, 為CO2強(qiáng)吸收匯。苔原TSB、苔原TNS和苔原TTR這3類不同苔原NEE的空間差異原因?yàn)楹xB(niǎo)活動(dòng)強(qiáng)度、植被覆蓋程度和地勢(shì)的差異[42-43]。苔原TSB的海鳥(niǎo)活動(dòng)密集, 為此區(qū)域輸入了C、N、P和S等營(yíng)養(yǎng)物質(zhì), 促進(jìn)了苔原植被的發(fā)育, 強(qiáng)烈吸收CO2[35, 43]。苔原TNS和苔原TTR的海鳥(niǎo)活動(dòng)較少, 影響微弱。其中, 地勢(shì)較低的苔原區(qū)域, 水量豐沛, 植被發(fā)育較好, 為CO2吸收匯; 地勢(shì)較高的區(qū)域, 土壤含水率較低, 植被稀疏, 對(duì)CO2的排放大于吸收, 為CO2排放源[42]。

    此后, 不同學(xué)者在新奧爾松(Ny-?lesund)不同苔原區(qū)域也開(kāi)展了溫室氣體的相關(guān)研究。2009年, 陳清清等[44]再次在新奧爾松苔原的鳥(niǎo)類保護(hù)區(qū)苔原(TSB)開(kāi)展了CO2凈通量的相關(guān)研究。發(fā)現(xiàn)海鳥(niǎo)活動(dòng)強(qiáng)烈區(qū)的CO2凈通量(NEE)為(–88.2±8.8) mg·m–2·h–1, 海鳥(niǎo)活動(dòng)中等的苔原區(qū)NEE為(14.6±9.1) mg·m–2·h–1, 海鳥(niǎo)活動(dòng)邊緣區(qū)的NEE為(42.2±13.3) mg·m–2·h–1, 苔原TSB的整體平均值(–5.1±8.3) mg·m–2·h–1, 為CO2吸收匯, 與本文研究結(jié)果一致。同時(shí), 陳清清等[44]也在與本研究的普通海灘苔原類似的苔原區(qū)域開(kāi)展了CO2凈通量(NEE)研究, 其平均NEE為(12.1±6.2) mg·m–2·h–1, 為CO2排放源, 與本文研究中結(jié)果極為相近。2014年和2015年, 李芳芳[45]也在普通海灘苔原TNS開(kāi)展了監(jiān)測(cè)CO2凈通量(NEE)的變化規(guī)律相關(guān)研究, 2014年NEE的波動(dòng)范圍為–131.5~ 73.2 mg·m–2·h–1, 平均值為(7.92±7.3) mg·m–2·h–1,為CO2排放源; 2015年NEE 的波動(dòng)范圍為–90.9~ 30.4 mg·m–2·h–1, 平均值為(8.9±5.1) mg·m–2·h–1, 為CO2排放源。誤差范圍內(nèi)與本文研究結(jié)果基本相同。

    將本研究結(jié)果與其他極地區(qū)域苔原研究結(jié)果進(jìn)行相比。鳥(niǎo)類保護(hù)區(qū)苔原TSB的觀測(cè)區(qū)中, 植被覆蓋度最高、強(qiáng)烈的海鳥(niǎo)活動(dòng)區(qū)HB和植被覆蓋度較高、中等強(qiáng)度的海鳥(niǎo)活動(dòng)區(qū)MB的CO2凈通量(NEE)處于低緯度北極的阿拉斯加和西伯利亞苔原區(qū)域已報(bào)告的數(shù)據(jù)范圍–63.8~–38.5 mg·m–2·h–1內(nèi)[22,46-47]。觀測(cè)區(qū)HB的NEE與南極海洋動(dòng)物聚集區(qū)苔原如企鵝和海豹聚居區(qū)(NEE為–37.2~ 5.2 mg·m–2·h–1)的結(jié)果相近, 且均為CO2強(qiáng)吸收匯[48]。鳥(niǎo)類保護(hù)區(qū)苔原TSB的整體NEE均處于格陵蘭不同植被區(qū)域的平均NEE通量范圍(–200~30 mg·m–2·h–1)[23, 49]。

    2.6 土壤理化性質(zhì)對(duì)新奧爾松(Ny-?lesund)苔原CO2凈通量(NEE)的影響

    CO2凈通量受氣溫、土壤溫度、光照強(qiáng)度、降水、植被覆蓋程度、植被類型等因素的綜合影響[40, 50-52]。土壤基質(zhì)含量決定了土壤CO2產(chǎn)生和同化的速率, 控制了CO2凈通量[53-57]。海鳥(niǎo)糞便能增加苔原土壤C、N、P等營(yíng)養(yǎng)物質(zhì)輸入, 增強(qiáng)植被的發(fā)育, 加強(qiáng)了苔原碳匯作用。同時(shí), 海鳥(niǎo)活動(dòng)增加了苔原植被N輸入, 能夠促進(jìn)該區(qū)域?qū)O2的吸收, 增強(qiáng)其碳匯功能[56-57]。鳥(niǎo)類保護(hù)區(qū)苔原TSB的CO2凈通量(NEE)與TN、TP、TC含量和C/N值有一定的相關(guān)性, 但比較微弱; 而與土壤NH4+-N(<0.05)和NO3–-N(<0.05)含量顯著負(fù)相關(guān)(圖3)。苔原TNS的NEE與土壤TC、TN、TP、NH4+-N、NO3–-N含量及C/N沒(méi)有明顯的相關(guān)關(guān)系。土壤的水分條件會(huì)影響土壤物質(zhì)的形態(tài)和轉(zhuǎn)化率及微生物活性等[50, 57]。土壤水分孔隙率(WFPS)≥60 %, 高于厭氧呼吸作用范圍, 降水和土壤水分會(huì)增加CO2同化作用增強(qiáng), CO2排放減弱或吸收增強(qiáng)[50, 57]。土壤水分含量較低, 呼吸作用可能超過(guò)光合作用, CO2同化作用增強(qiáng), 增加CO2的排放[50, 57], 因此水分含量的高低對(duì)土壤CO2的排放起著重要作用。鳥(niǎo)類保護(hù)區(qū)苔原(TSB)和普通海灘苔原(TNS)的CO2凈通量(NEE)與土壤含水率呈顯著負(fù)相關(guān)關(guān)系(<0.01) (圖3)。

    土壤溫度是影響生態(tài)系統(tǒng)CO2排放的重要因素[35]。植被覆蓋度高、海鳥(niǎo)活動(dòng)強(qiáng)烈的觀測(cè)區(qū)HB的CO2凈通量(NEE) 與0~15 cm平均土壤溫度顯著正相關(guān)(=0.32,=0.06)(圖4)。在普通海灘苔原(TNS)區(qū)域, 地勢(shì)低的濕地苔原觀測(cè)區(qū)SW的CO2凈通量(NEE)與0~15 cm土壤平均溫度呈顯著極負(fù)相關(guān)關(guān)系(≤0.01), 地勢(shì)略高的半干旱苔原觀測(cè)區(qū)SM的CO2凈通量(NEE)與0~15 cm土壤平均溫度呈極顯著正相關(guān)(≤0.01), 地勢(shì)更高的干旱苔原觀測(cè)區(qū)SD的CO2凈通量(NEE)與0~15 cm土壤平均溫度呈顯著正相關(guān)(≤0.05), 而與地勢(shì)高的高地苔原觀測(cè)區(qū)UD的NEE與0~15 cm土壤平均溫度也具有弱的負(fù)相關(guān)關(guān)系(=–0.41,= 0.09)。

    3 本研究的不確定性分析

    極地生長(zhǎng)季溫室氣體的排放, 尤其高密度海鳥(niǎo)活動(dòng)的7—8月份, 對(duì)全年溫室氣體排放及碳循環(huán)起著重要的作用。本研究選取了新奧爾松夏季無(wú)冰區(qū)的鳥(niǎo)類保護(hù)區(qū)苔原(TSB)、普通海灘苔原(TNS)和過(guò)渡苔原(TTR), 監(jiān)測(cè)CO2凈通量(NEE)的空間變化規(guī)律, 整個(gè)觀測(cè)期持續(xù)10天左右, 本研究的不確定性分析也聚焦于此。本研究?jī)H關(guān)注北極高緯度新奧爾松苔原的生長(zhǎng)期短暫的夏季CO2凈通量的空間規(guī)律, 采用密閉箱法監(jiān)測(cè)體系, 觀測(cè)時(shí)間均在9:00—11:00(當(dāng)?shù)貢r(shí)間), 與Oechel 等[58]以及Lafleur 和Humphreys[24]的監(jiān)測(cè)方法體系一致。本研究結(jié)果與其他在北極苔原7月和8月的研究結(jié)果具有可比較性, 如前述的格陵蘭、阿拉斯加、西伯利亞和南極苔原區(qū)域[22-23, 46-48, 59-61]。

    圖3 鳥(niǎo)類保護(hù)區(qū)苔原(TSB, n=16)和普通海灘苔原(TNS, n=29)NEE與土壤理化學(xué)性質(zhì)(土壤中NH4+-N、NO3–-N、TN、TC、TP、C∶N及土壤含水率)的相關(guān)性

    Fig.3.The correlation between NEE and soil chemical properties (NH4+-N, NO3–-N, TN, TC, TP, C∶N ratio and soil moisture) at the different observation sites in seabird sanctuary tundra TSB (=16) and beach tundra TNS (=29)

    4 結(jié)論

    1.鳥(niǎo)類保護(hù)區(qū)苔原TSB整體為CO2吸收匯, 植被覆蓋度高、海鳥(niǎo)活動(dòng)強(qiáng)烈的觀測(cè)區(qū)HB對(duì)CO2的吸收顯著高于植被覆蓋度和海鳥(niǎo)活動(dòng)中等觀測(cè)區(qū)MB, 觀測(cè)區(qū)HB和MB的NEE均顯著低于植被覆蓋度低、海鳥(niǎo)活動(dòng)邊緣觀測(cè)區(qū)LB。苔原植被發(fā)育越好, 海鳥(niǎo)活動(dòng)越密集, 對(duì)CO2吸收越強(qiáng)烈。

    2.非海鳥(niǎo)聚居區(qū)的普通海灘苔原TNS和過(guò)渡苔原TTR均為CO2排放源。苔原TNS中, 地勢(shì)越低越濕潤(rùn)的苔原區(qū)域?qū)O2吸收越強(qiáng), 如從地勢(shì)略高的半干旱苔原觀測(cè)區(qū)SM到地勢(shì)低的濕地苔原觀測(cè)區(qū)SW; 地勢(shì)越高越干旱的苔原觀測(cè)區(qū)的CO2排放越強(qiáng)烈, 如從地勢(shì)高的高地苔原觀測(cè)區(qū)UD到地勢(shì)最高的頂部苔原觀測(cè)區(qū)UL。苔原TTR中泥炭區(qū)為強(qiáng)吸收匯, 而高地苔原區(qū)域?yàn)閺?qiáng)排放源。

    圖4 鳥(niǎo)類保護(hù)區(qū)苔原(TSB)和普通海灘苔原(TNS)的NEE與土壤溫度的相關(guān)性.*表示P≤0.05, **表示P≤0.01顯著相關(guān)

    Fig.4.The correlation between NEE and soil temperature at the different observation sites in tundra TSB and tundra TNS.Significant correlations were indicated by *(≤0.05) and ** (≤0.01)

    3.土壤理化性質(zhì)會(huì)對(duì)新奧爾松(Ny-?lesund)苔原CO2凈通量(NEE)產(chǎn)生重要影響。鳥(niǎo)類保護(hù)區(qū)苔原TSB和普通海灘苔原TNS的NEE與土壤水分顯著負(fù)相關(guān)(<0.01); 鳥(niǎo)類保護(hù)區(qū)苔原TSB的NEE與土壤溫度顯著正相關(guān), 與NH4+- N(<0.05)、NO3–-N(<0.05)含量顯著負(fù)相關(guān)。

    致謝感謝賈楠博士在第4次北極科學(xué)考察中采集溫室氣體和土壤樣品, 感謝梁有慶高級(jí)工程師在實(shí)驗(yàn)上的指導(dǎo)。

    1 VIRKKALA A M, VIRTANEN T, LEHTONEN A, et al.The current state of CO2flux chamber studies in the Arctic tundra[J].Progress in Physical Geography: Earth and Environment, 2018, 42(2): 162-184.

    2 HUGELIUS G, STRAUSS J, ZUBRZYCKI S, et al.Estimated stocks of circumpolar permafrost carbon with quantified uncertainty ranges and identified data gaps[J].Biogeosciences, 2014, 11(23): 6573-6593.

    3 TARNOCAI C, CANADELL J G, SCHUUR E A G, et al.Soil organic carbon pools in the northern circumpolar permafrost region[J].Global Biogeochemical Cycles, 2009, 23(2): GB2023.DOI: 10.1029/2008GB003327.

    4 WATSON R T, NOBLE I R, BOLIN B, et al.Land use, land-use change and forestry.A special report of the international panel on climate change (IPCC)[EB/OL].(2000).https://www.ipcc.ch/report/land-use-land-use-change-and-forestry/.

    5 MARUSHCHAK M E, KIEPE I, BIASI C, et al.Carbon dioxide balance of subarctic tundra from plot to regional scales[J].Biogeosciences, 2013, 10(1): 437-452.

    6 BLAIS A M, LORRAIN S, TREMBLAY A.Greenhouse gas fluxes (CO2, CH4and N2O) in forests and wetlands of boreal, temperate and tropical regions[M]//TREMBLAY A, VARFALVY L, ROEHM C.et al.Greenhouse gas emissions–Fluxes and processes.Berlin, Heidelberg: Springer, 2005: 87-127.

    7 HEIKKINEN J E P, ELSAKOV V, MARTIKAINEN P J.Carbon dioxide and methane dynamics and annual carbon balance in tundra wetland in NE Europe, Russia[J].Global Biogeochemical Cycles, 2002, 16(4): 62-1-62-9.DOI:10.1029/2002GB001930.

    8 BLACKALL T D, WILSON L J, THEOBALD M R, et al.Ammonia emissions from seabird colonies[J].Geophysical Research Letters, 2007, 34(10): L10801.

    9 RIDDICK S N, DRAGOSITS U, BLACKALL T D, et al.The global distribution of ammonia emissions from seabird colonies[J].Atmospheric Environment, 2012, 55: 319-327.

    10 KRUPA S V.Effects of atmospheric ammonia (NH3) on terrestrial vegetation: A review[J].Environmental Pollution, 2003, 124(2): 179-221.

    11 LINDEBOOM H J.The nitrogen pathway in a penguin rookery[J].Ecology, 1984, 65(1): 269-277.

    12 MIZUTANI H, KABAYA Y, WADA E.Ammonia volatilization and high15N/14N ratio in a penguin rookery in Antarctica[J].Geochemical Journal, 1985, 19(6): 323-327.

    13 SUN L G, ZHU R B, XIE Z Q, et al.Emissions of nitrous oxide and methane from Antarctic Tundra: Role of penguin dropping deposition[J].Atmospheric Environment, 2002, 36(31): 4977-4982.DOI: 10.1016/S1352-2310(02)00340-0.

    14 PARK J H, DAY T A, STRAUSS S, et al.Biogeochemical pools and fluxes of carbon and nitrogen in a maritime tundra near penguin colonies along the Antarctic Peninsula[J].Polar Biology, 2007, 30(2): 199-207.

    15 ZHU R B, LIU Y S, XU H, et al.Methane emissions from three sea animal colonies in the maritime Antarctic[J].Atmospheric Environment, 2008, 42(6): 1197-1205.

    16 ZHU R B, LIU Y S, XU H, et al.Nitrous oxide emissions from sea animal colonies in the maritime Antarctic[J].Geophysical Research Letters, 2008, 35(9): L09807.DOI:10.1029/2007GL032541.

    17 WüTHRICH C, M?LLER I, THANNHEISER D.CO2-fluxes in different plant communities of a high-Arctic tundra watershed (Western Spitsbergen)[J].Journal of Vegetation Science, 1999, 10(3): 413-420.

    18 LLOYD C R.The measurement and modelling of the carbon dioxide exchange at a high Arctic site in Svalbard[J].Global Change Biology, 2001, 7(4): 405-426.

    19 CANNONE N, BINELLI G, WORLAND M R, et al.CO2fluxes among different vegetation types during the growing season in Marguerite Bay (Antarctic Peninsula)[J].Geoderma, 2012, 189-190: 595-605.DOI:10.1016/j.geoderma.2012.06.026.

    20 ZHU R B, BAO T, WANG Q, et al.Summertime CO2fluxes and ecosystem respiration from marine animal colony tundra in maritime Antarctica[J].Atmospheric Environment, 2014, 98: 190-201.

    21 LIU Y S, ZHANG W Y, ZHU R B, et al.Potential methane and nitrous oxide production and respiration rates from penguin and seal colony tundra soils during freezing–thawing cycles under different water contents in coastal Antarctica[J].Advances in Polar Science, 2017, 28(1): 61-74.DOI: 10.13679/j.advps.2017.1.00061.

    22 HARAZONO Y, MANO M, MIYATA A, et al.Inter-annual carbon dioxide uptake of a wet sedge tundra ecosystem in the Arctic[J].Tellus B: Chemical and Physical Meteorology, 2003, 55(2): 215-231.

    23 GROENDAHL L, FRIBORG T, SOEGAARD H.Temperature and snow-melt controls on interannual variability in carbon exchange in the high Arctic[J].Theoretical and Applied Climatology, 2007, 88(1/2): 111-125.

    24 LAFLEUR P M, HUMPHREYS E R.Spring warming and carbon dioxide exchange over low Arctic tundra in central Canada[J].Global Change Biology, 2008, 14(4): 740-756.

    25 OECHEL W C, HASTINGS S J, VOURLRTIS G, et al.Recent change of Arctic tundra ecosystems from a net carbon dioxide sink to a source[J].Nature, 1993, 361(6412): 520-523.

    26 OECHEL W C, VOURLITIS G L, HASTINGS S J, et al.Acclimation of ecosystem CO2exchange in the Alaskan Arctic in response to decadal climate warming[J].Nature, 2000, 406(6799): 978-981.

    27 BIRKS H J B, JONES V J, ROSE N L.Recent environmental change and atmospheric contamination on Svalbard as recorded in lake sediments–An introduction[J].Journal of Paleolimnology, 2004, 31(4): 403-410.

    28 SCHWERDTFEGER W.Weather and climate of the Antarctic[M].Amsterdam: Elsevier, 1984.

    29 鄧海濱, 陸龍驊, 卞林根.北極苔原新奧爾松地區(qū)的地表輻射特征[J].極地研究, 2006, 18(4): 254-264.

    30 劉雅淑.南極苔原溫室氣體通量時(shí)空變化規(guī)律與產(chǎn)生機(jī)理研究[D].合肥: 中國(guó)科學(xué)技術(shù)大學(xué),2011.

    31 袁林喜, 龍楠燁, 謝周清, 等.北極新奧爾松地區(qū)現(xiàn)代污染源及其指示植物研究[J].極地研究, 2006, 18(1): 9-20.

    32 鄧海濱, 陸龍驊, 卞林根.北極苔原Ny-?lesund地區(qū)短期氣候特征[J].極地研究, 2005, 17(1): 32-44.

    33 WINTHER J G, GODTLIEBSEN F, GERLAND S, et al.Surface albedo in Ny-?lesund, Svalbard: Variability and trends during 1981-1997[J].Global and Planetary Change, 2002, 32(2/3): 127-139.

    34 WESTERMANN S, WOLLSCHL?GER U, BOIKE J.Monitoring of active layer dynamics at a permafrost site on Svalbard using multi-channel ground-penetrating radar[J].The Cryosphere, 2010, 4(4): 475-487.

    35 朱仁斌, 孫立廣, 謝周清, 等.南極菲爾德斯半島植被微區(qū)CO2濃度的監(jiān)測(cè)[J].環(huán)境科學(xué), 2001, 22(4): 6-10.

    36 HUTCHINSON G L, MOSIER A R.Improved soil cover method for field measurement of nitrous oxide fluxes[J].Soil Science Society of America Journal, 1981, 45(2): 311-316.

    37 VAN HUISSTEDEN J, MAXIMOV T C, DOLMAN A J.High methane flux from an Arctic floodplain (Indigirka Lowlands, eastern Siberia)[J].Journal of Geophysical Research: Biogeosciences, 2005, 110(G2): G02002.DOI:10.1029/2005JG000010.

    38 郎漫, 蔡祖聰.百菌清對(duì)土壤氧化亞氮和二氧化碳排放的影響[J].應(yīng)用生態(tài)學(xué)報(bào), 2008, 19(12): 2745-2750.

    39 KEENEY D R, NELSON D W.Nitrogen-inorganic forms[M]//PAGE A L.Methods of soil analysis: Part 2 Chemical and microbiological properties.2nd ed.New Jersey: Wiley,1983: 643–698.

    40 SUN L, XIE Z, ZHAO J.A 3, 000-year record of penguin populations[J].Nature, 2000, 407(6806): 858.

    41 ZHU R B, SUN J J, LIU Y S, et al.Potential ammonia emissions from penguin guano, ornithogenic soils and seal colony soils in coastal Antarctica: Effects of freezing-thawing cycles and selected environmental variables[J].Antarctic Science, 2011, 23(1): 78-92.

    42 MOURA P A, FRANCELINO M R, SCHAEFER C E G R, et al.Distribution and characterization of soils and landform relationships in Byers Peninsula, Livingston Island, Maritime Antarctica[J].Geomorphology, 2012, 155-156: 45-54.

    43 RAICH J W, TUFEKCIOGUL A.Vegetation and soil respiration: Correlations and controls[J].Biogeochemistry, 2000, 48(1): 71-90.

    44 陳清清, 朱仁斌, 丁瑋, 等.海鳥(niǎo)活動(dòng)對(duì)北極新奧爾松苔原-大氣CO2交換的影響[J].極地研究, 2012, 24(3): 254-265.

    45 李芳芳.高緯度北極苔原斷面溫室氣體通量變異規(guī)律研究[D].合肥: 中國(guó)科學(xué)技術(shù)大學(xué), 2017.

    46 VOURLITIS G L, OECHEL W C.Eddy covariance measurements of CO2and energy fluxes of an Alaskan tussock tundra ecosystem[J].Ecology, 1999, 80(2): 686-701.

    47 VAN DER MOLEN M K, VAN HUISSTEDEN J, PARMENTIER F J W, et al.The growing season greenhouse gas balance of a continental tundra site in the Indigirka Lowlands, NE Siberia[J].Biogeosciences, 2007, 4(6): 985-1003.

    48 ZHU R B, MA D W, XU H.Summertime N2O, CH4and CO2exchanges from a tundra marsh and an upland tundra in maritime Antarctica[J].Atmospheric Environment, 2014, 83: 269-281.

    49 CHRISTENSEN T R, FRIBORG T, SOMMERKORN M, et al.Trace gas exchange in a high-Arctic valley: 1.Variationsin CO2and CH4flux between tundra vegetation types[J].Global Biogeochemical Cycles, 2000, 14(3): 701-713.

    50 賈志軍, 宋長(zhǎng)春.濕地生態(tài)系統(tǒng)CO2凈交換、水汽通量及二者關(guān)系淺析[J].生態(tài)與農(nóng)村環(huán)境學(xué)報(bào), 2006, 22(2): 75-79.

    51 宋長(zhǎng)春, 閻百興, 王躍思, 等.三江平原沼澤濕地CO2和CH4通量及影響因子[J].科學(xué)通報(bào), 2003, 48(23): 2473-2477.

    52 李海防, 夏漢平, 熊燕梅, 等.土壤溫室氣體產(chǎn)生與排放影響因素研究進(jìn)展[J].生態(tài)環(huán)境, 2007, 16(6): 1781-1788.

    53 SUBKE J A, REICHSTEIN M, TENHUNEN J D.Explaining temporal variation in soil CO2efflux in a mature spruce forest in Southern Germany[J].Soil Biology and Biochemistry, 2003, 35(11): 1467-1483.

    54 MICHAELIS L, MENTEN M L, GOODY R S, et al.The kinetics of invertase action[J].Biochemistry, 1913, 49(2): 333–369.

    55 STITT M, KRAPP A.The interaction between elevated carbon dioxide and nitrogen nutrition: The physiological and molecular background[J].Plant, Cell & Environment, 1999, 22(6): 583-621.DOI: 10.1046/j.1365-3040.1999.00386.x.

    56 MAGNANI F, MENCUCCINI M, BORGHETTI M, et al.The human footprint in the carbon cycle of temperate and boreal forests[J].Nature, 2007, 447(7146): 849-851.

    57 陳全勝, 李凌浩, 韓興國(guó), 等.典型溫帶草原群落土壤呼吸溫度敏感性與土壤水分的關(guān)系[J].生態(tài)學(xué)報(bào), 2004, 24(4): 831-836.

    58 OECHEL W C, VOURLITIS G, HASTINGS S J.Cold season CO2emission from Arctic soils[J].Global Biogeochemical Cycles, 1997, 11(2): 163-172.

    59 RHEW R C, TEH Y A, ABEL T.Methyl halide and methane fluxes in the northern Alaskan coastal tundra[J].Journal of Geophysical Research: Biogeosciences, 2007, 112(G2): G02009.DOI:10.1029/2006JG000314.

    60 VON FISCHER J C, RHEW R C, AMES G M, et al.Vegetation height and other controls of spatial variability in methane emissions from the Arctic coastal tundra at Barrow, Alaska[J].Journal of Geophysical Research: Biogeosciences, 2010, 115(G4): G00I03.DOI:10.1029/2009JG001283.

    61 SOEGAARD H, NORDSTROEM C, FRIBORG T, et al.Trace gas exchange in a high-Arctic valley: 3.Integrating and scaling CO2fluxes from canopy to landscape using flux data, footprint modeling, and remote sensing[J].Global Biogeochemical Cycles, 2000, 14(3): 725-744.

    The spatial variation of net CO2fluxes (NEE) and their influence factors in tundra in high Arctic

    Chen Qingqing1,2, Bao Tao3, Zhu Renbin2, Xu Hua4

    (1College of Life Science, Resources and Environment, Yichun University, Yichun 336000, China;2Anhui Province Key Laboratory of Polar Environment and Global Change, School of Earth and Space Science, University of Science and Technology of China, Hefei 230026, China;3Key Laboratory of Regional Climate-Environment for Temperate East Asia, Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing 100029, China;4State Key Laboratory of Soil and Sustainable Agriculture, Institute of Soil Science, Chinese Academy of Sciences, Nanjing 210091, China)

    The release of CO2from summer ice-free tundra of the high Arctic plays an important role in the global carbon cycle.These net ecosystem exchanges (NEEs) and their influence factors from Seabird Sanctuary Tundra (TSB), Tundra in Non-Seabird Colony (TNS) and Tundra in Transition Zone (TTR) were measured using the closed chamber method during the Fourth Arctic Scientific Exploration (26 July–5 August 2008) in the High Arctic.The spatial variations of NEE show that TSB was a sink for CO2, with an average NEE of (?39.0±6.0) mg·m?2·h?1.Typically, absorption of CO2in areas with high vegetation cover and strong seabird usage was greater than in areas with lower vegetation cover and seabird usage.In contrast, TNS and TTR were emission sources for CO2, with average NEEs of (12.0±6.3) mg·m?2·h?1and (40.5±29.3) mg·m?2·h?1, respectively.Highland TTR ((106.4±23.1) mg·m?2·h?1) was the strongest emission source, while peatland TTR [(?58.3±9.5) mg·m?2·h?1] was a strong sink.Spatial variation of NEEs was linked to differences in vegetation cover and hydrological conditions related to seabird activity in all tundra types.Physical and chemical properties of tundra soils also affected NEE values, showing negative correlations with soil moisture (=?0.44,=0.003) in TSB and TNS, but a positive correlation with soil temperature (=0.32,=0.06) and negative correlations with NH4+-N (<0.05) and NO3–-N (<0.05) concentrations in TSB.

    CO2,net CO2fluxes, Arctic, Ny-?lesund, high Arctic tundra

    2020年11月收到來(lái)稿, 2021年1月收到修改稿

    國(guó)家重點(diǎn)研發(fā)計(jì)劃重點(diǎn)專項(xiàng)(2020YFA0608500)、國(guó)家自然科學(xué)基金(41776190,41976220)資助

    陳清清, 女, 1984年生。博士, 主要從事元素生物地球化學(xué)循環(huán)研究。E-mail: chenqq@nju.edu.cn

    朱仁斌, E-mail: zhurb@ustc.edu.cn

    10.13679/j.jdyj.20200073

    猜你喜歡
    苔原海鳥(niǎo)北極
    生活在西伯利亞苔原上的人們
    照相機(jī)(2022年12期)2022-02-09 09:13:48
    In Russia's far north,legends and lives are frozen in time
    北極有個(gè)“放屁湖”
    北極兔乖乖,唱起歌來(lái)
    海鳥(niǎo)遇上惡房東
    極地苔原
    北極
    心聲歌刊(2017年4期)2017-09-20 11:43:48
    去北極嘍!
    幼兒園(2016年22期)2017-03-01 13:57:52
    苔原地區(qū)
    小布老虎(2016年12期)2016-12-01 05:47:09
    塑料垃圾成海鳥(niǎo)最大殺手
    欧美黑人欧美精品刺激| 黄片无遮挡物在线观看| 日本vs欧美在线观看视频| 国产精品女同一区二区软件| 国产无遮挡羞羞视频在线观看| 国产精品香港三级国产av潘金莲 | 天堂中文最新版在线下载| av免费观看日本| 一级毛片电影观看| 老司机在亚洲福利影院| 男女无遮挡免费网站观看| 欧美日韩综合久久久久久| 久久精品久久久久久噜噜老黄| 天堂8中文在线网| 老司机影院成人| 丁香六月天网| 久久精品人人爽人人爽视色| 少妇人妻久久综合中文| 亚洲国产欧美网| 国产精品蜜桃在线观看| 久久女婷五月综合色啪小说| 桃花免费在线播放| av不卡在线播放| 十八禁高潮呻吟视频| 看免费成人av毛片| 精品一区二区三区四区五区乱码 | 91精品国产国语对白视频| 999精品在线视频| 亚洲精华国产精华液的使用体验| 国产精品偷伦视频观看了| 激情五月婷婷亚洲| 99热全是精品| 免费在线观看完整版高清| 亚洲精品自拍成人| 黑人欧美特级aaaaaa片| 久久久久久人人人人人| 国产成人精品福利久久| 精品一区二区三区四区五区乱码 | 曰老女人黄片| 另类亚洲欧美激情| 亚洲天堂av无毛| 国产成人啪精品午夜网站| 桃花免费在线播放| 水蜜桃什么品种好| 天美传媒精品一区二区| 国产亚洲最大av| 欧美中文综合在线视频| 成人国语在线视频| 欧美少妇被猛烈插入视频| 国产激情久久老熟女| 精品少妇久久久久久888优播| 亚洲成人手机| 久久97久久精品| 好男人视频免费观看在线| 美女国产高潮福利片在线看| 韩国高清视频一区二区三区| 亚洲图色成人| 日韩欧美一区视频在线观看| 宅男免费午夜| 久久免费观看电影| 天美传媒精品一区二区| 成年av动漫网址| videosex国产| 一本—道久久a久久精品蜜桃钙片| 午夜91福利影院| 国产乱人偷精品视频| 成人亚洲精品一区在线观看| 性少妇av在线| 亚洲四区av| 在线观看一区二区三区激情| 色精品久久人妻99蜜桃| 综合色丁香网| 18禁观看日本| 久久人妻熟女aⅴ| 中文字幕色久视频| 国产精品免费视频内射| 国产在线视频一区二区| 久久这里只有精品19| 精品国产国语对白av| 日日啪夜夜爽| 男女下面插进去视频免费观看| av网站在线播放免费| 亚洲国产精品成人久久小说| h视频一区二区三区| 亚洲少妇的诱惑av| 中文欧美无线码| 大香蕉久久成人网| 新久久久久国产一级毛片| 性色av一级| 亚洲伊人色综图| 亚洲精品久久成人aⅴ小说| 亚洲视频免费观看视频| 少妇被粗大的猛进出69影院| 色网站视频免费| 超色免费av| 一本色道久久久久久精品综合| 久久久久视频综合| 午夜福利视频精品| 国产免费又黄又爽又色| 中文欧美无线码| 丁香六月欧美| 国产成人系列免费观看| 免费人妻精品一区二区三区视频| 大片免费播放器 马上看| 午夜福利视频在线观看免费| 久久天躁狠狠躁夜夜2o2o | 女人爽到高潮嗷嗷叫在线视频| 午夜福利视频在线观看免费| 汤姆久久久久久久影院中文字幕| 在线观看免费高清a一片| 亚洲精品一区蜜桃| 亚洲欧美精品自产自拍| 中国国产av一级| 欧美亚洲日本最大视频资源| 欧美日韩亚洲综合一区二区三区_| 国产精品人妻久久久影院| a 毛片基地| 久久av网站| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲精品国产av成人精品| 欧美97在线视频| 又大又爽又粗| 午夜福利视频精品| 国产成人系列免费观看| 丝袜人妻中文字幕| 99久久综合免费| 天天躁狠狠躁夜夜躁狠狠躁| 只有这里有精品99| 99九九在线精品视频| 日韩电影二区| 十八禁人妻一区二区| 久久久久久久精品精品| 免费黄频网站在线观看国产| 国产精品国产三级专区第一集| 精品少妇久久久久久888优播| 亚洲精品一区蜜桃| 汤姆久久久久久久影院中文字幕| 亚洲色图综合在线观看| 无遮挡黄片免费观看| 欧美精品一区二区大全| 久久97久久精品| 日本91视频免费播放| 伊人亚洲综合成人网| 在线观看国产h片| 中文字幕亚洲精品专区| 国产成人精品久久二区二区91 | 精品国产一区二区三区四区第35| 午夜av观看不卡| 好男人视频免费观看在线| 美女中出高潮动态图| 国产精品国产三级专区第一集| 黄色毛片三级朝国网站| 在线观看免费视频网站a站| 欧美日韩成人在线一区二区| 国语对白做爰xxxⅹ性视频网站| 男女无遮挡免费网站观看| 国产精品成人在线| av电影中文网址| 少妇被粗大的猛进出69影院| 99久久精品国产亚洲精品| 一边摸一边抽搐一进一出视频| 国产伦理片在线播放av一区| 女人精品久久久久毛片| 天堂俺去俺来也www色官网| 国产精品.久久久| 亚洲成国产人片在线观看| 亚洲欧美精品综合一区二区三区| 成人亚洲欧美一区二区av| 多毛熟女@视频| 黄色 视频免费看| 欧美激情高清一区二区三区 | 国产一区二区激情短视频 | 操美女的视频在线观看| 国产精品99久久99久久久不卡 | 欧美 亚洲 国产 日韩一| 国产不卡av网站在线观看| 国产精品一国产av| 国产精品久久久久久精品电影小说| 人人妻人人爽人人添夜夜欢视频| 狂野欧美激情性xxxx| 亚洲三区欧美一区| 欧美xxⅹ黑人| 精品久久蜜臀av无| 色综合欧美亚洲国产小说| 岛国毛片在线播放| 亚洲av日韩精品久久久久久密 | 欧美黑人精品巨大| 男女边摸边吃奶| 国产激情久久老熟女| 亚洲精品第二区| 伊人久久大香线蕉亚洲五| 亚洲男人天堂网一区| 国产成人精品久久二区二区91 | 97在线人人人人妻| av免费观看日本| www日本在线高清视频| 国产成人一区二区在线| 精品午夜福利在线看| 久久久精品区二区三区| 成人毛片60女人毛片免费| 人人妻人人添人人爽欧美一区卜| 日韩一区二区视频免费看| 国产熟女午夜一区二区三区| 中文乱码字字幕精品一区二区三区| 亚洲精华国产精华液的使用体验| 亚洲av男天堂| 国产精品久久久久久久久免| 婷婷色麻豆天堂久久| 午夜福利视频精品| 1024香蕉在线观看| 中国三级夫妇交换| 国产精品蜜桃在线观看| 免费高清在线观看视频在线观看| 亚洲自偷自拍图片 自拍| 国产成人一区二区在线| 国产亚洲午夜精品一区二区久久| 精品亚洲成国产av| 免费高清在线观看视频在线观看| 亚洲色图综合在线观看| 国产精品久久久久成人av| av网站免费在线观看视频| 欧美精品一区二区大全| 精品久久久久久电影网| 国产一区二区激情短视频 | 黄色毛片三级朝国网站| 亚洲精品国产av蜜桃| 在线天堂中文资源库| 亚洲熟女精品中文字幕| 看免费av毛片| 下体分泌物呈黄色| 久久毛片免费看一区二区三区| 51午夜福利影视在线观看| www.精华液| 丰满饥渴人妻一区二区三| 久久国产精品男人的天堂亚洲| 一级,二级,三级黄色视频| 久久午夜综合久久蜜桃| 黄网站色视频无遮挡免费观看| 人妻一区二区av| 久久人妻熟女aⅴ| 亚洲av男天堂| 色播在线永久视频| 丁香六月天网| 国产男女超爽视频在线观看| 最新的欧美精品一区二区| 大片电影免费在线观看免费| 在线观看免费高清a一片| 国产精品久久久久久精品古装| 午夜福利,免费看| 亚洲av国产av综合av卡| 中文精品一卡2卡3卡4更新| 99久久99久久久精品蜜桃| 熟妇人妻不卡中文字幕| 免费高清在线观看日韩| 99久久精品国产亚洲精品| 国产在视频线精品| 一区二区日韩欧美中文字幕| 亚洲精品aⅴ在线观看| 卡戴珊不雅视频在线播放| 久久久欧美国产精品| 久久久久视频综合| 精品亚洲成国产av| 精品一区二区三区四区五区乱码 | 中文字幕制服av| 国产黄色免费在线视频| 久久久久精品性色| 亚洲欧美一区二区三区黑人| 99久久综合免费| 国产精品二区激情视频| 国产精品一区二区在线不卡| 天天躁狠狠躁夜夜躁狠狠躁| av片东京热男人的天堂| 韩国精品一区二区三区| 国产 一区精品| 少妇人妻精品综合一区二区| 51午夜福利影视在线观看| 在线观看免费高清a一片| 亚洲国产精品999| 国产男女内射视频| 欧美97在线视频| 免费不卡黄色视频| av线在线观看网站| 一区二区av电影网| 999久久久国产精品视频| 亚洲人成电影观看| 国产日韩欧美视频二区| 亚洲激情五月婷婷啪啪| 亚洲欧美一区二区三区国产| 中文天堂在线官网| 丝袜脚勾引网站| 美女脱内裤让男人舔精品视频| 天美传媒精品一区二区| 精品少妇黑人巨大在线播放| 人人澡人人妻人| 秋霞伦理黄片| 高清av免费在线| 丝袜喷水一区| 亚洲av成人精品一二三区| 婷婷成人精品国产| 街头女战士在线观看网站| 国产精品一区二区精品视频观看| 精品福利永久在线观看| 一边摸一边抽搐一进一出视频| 人体艺术视频欧美日本| 王馨瑶露胸无遮挡在线观看| 亚洲中文av在线| 观看美女的网站| 赤兔流量卡办理| 精品午夜福利在线看| 女人高潮潮喷娇喘18禁视频| 精品一区二区免费观看| 国产欧美日韩一区二区三区在线| 一区在线观看完整版| 建设人人有责人人尽责人人享有的| 成人毛片60女人毛片免费| 无限看片的www在线观看| 五月开心婷婷网| 我的亚洲天堂| 国产精品二区激情视频| 精品少妇久久久久久888优播| 熟女av电影| 国产精品三级大全| av视频免费观看在线观看| 大码成人一级视频| 在线观看一区二区三区激情| 国产成人欧美在线观看 | 美女脱内裤让男人舔精品视频| 91精品国产国语对白视频| 国产一区二区三区综合在线观看| 国产高清国产精品国产三级| 80岁老熟妇乱子伦牲交| 黑人猛操日本美女一级片| 国产免费又黄又爽又色| 一边亲一边摸免费视频| 午夜福利乱码中文字幕| 久久久精品94久久精品| 久久精品熟女亚洲av麻豆精品| 精品午夜福利在线看| 少妇精品久久久久久久| xxx大片免费视频| 国产亚洲av片在线观看秒播厂| 精品一区二区三区四区五区乱码 | h视频一区二区三区| 又粗又硬又长又爽又黄的视频| 国产成人系列免费观看| av网站在线播放免费| 日本欧美视频一区| 啦啦啦在线观看免费高清www| 日本91视频免费播放| 国产 一区精品| 久久 成人 亚洲| 日韩av免费高清视频| 在线 av 中文字幕| 黄片小视频在线播放| 天天躁夜夜躁狠狠久久av| 人人妻人人澡人人看| 国语对白做爰xxxⅹ性视频网站| 91成人精品电影| 超色免费av| 精品少妇内射三级| 亚洲美女搞黄在线观看| 亚洲av电影在线进入| 成人国语在线视频| 热99国产精品久久久久久7| 国产淫语在线视频| 91国产中文字幕| 亚洲国产精品一区三区| 日韩中文字幕欧美一区二区 | 91精品三级在线观看| 久久久久视频综合| 性色av一级| 国产av码专区亚洲av| 丰满迷人的少妇在线观看| 午夜精品国产一区二区电影| 免费高清在线观看日韩| 曰老女人黄片| 精品卡一卡二卡四卡免费| 欧美97在线视频| 成人毛片60女人毛片免费| 亚洲色图综合在线观看| 久久精品熟女亚洲av麻豆精品| 亚洲精品久久久久久婷婷小说| 亚洲av电影在线观看一区二区三区| av在线app专区| 亚洲精品日韩在线中文字幕| 在线观看免费高清a一片| 日韩一区二区视频免费看| 亚洲国产av影院在线观看| 麻豆乱淫一区二区| 欧美激情高清一区二区三区 | 精品少妇久久久久久888优播| 天堂俺去俺来也www色官网| 男女下面插进去视频免费观看| 国产免费一区二区三区四区乱码| 另类亚洲欧美激情| 美女扒开内裤让男人捅视频| 一个人免费看片子| 99精品久久久久人妻精品| 丝袜脚勾引网站| 亚洲七黄色美女视频| 卡戴珊不雅视频在线播放| 秋霞伦理黄片| 看非洲黑人一级黄片| netflix在线观看网站| 日韩电影二区| 中文欧美无线码| 亚洲精品国产av成人精品| 欧美激情极品国产一区二区三区| 中国三级夫妇交换| 一级a爱视频在线免费观看| 国产成人精品在线电影| 又粗又硬又长又爽又黄的视频| 狠狠精品人妻久久久久久综合| 成年美女黄网站色视频大全免费| xxxhd国产人妻xxx| 亚洲成人免费av在线播放| 99精品久久久久人妻精品| 欧美日韩av久久| 国产女主播在线喷水免费视频网站| 国产精品久久久久久精品电影小说| 日韩av在线免费看完整版不卡| 天堂俺去俺来也www色官网| 观看av在线不卡| 丝瓜视频免费看黄片| 亚洲一级一片aⅴ在线观看| 美女中出高潮动态图| 久久久久精品国产欧美久久久 | 亚洲精品中文字幕在线视频| 国产极品粉嫩免费观看在线| 欧美av亚洲av综合av国产av | 九草在线视频观看| 亚洲综合精品二区| 交换朋友夫妻互换小说| 最近2019中文字幕mv第一页| av.在线天堂| 视频区图区小说| 黑人欧美特级aaaaaa片| 亚洲国产成人一精品久久久| 香蕉丝袜av| 国产日韩一区二区三区精品不卡| 在线观看人妻少妇| 69精品国产乱码久久久| 国产 精品1| 精品国产露脸久久av麻豆| 久久韩国三级中文字幕| 亚洲国产精品999| 男女高潮啪啪啪动态图| 国产免费又黄又爽又色| 日本vs欧美在线观看视频| 欧美成人精品欧美一级黄| 亚洲精品视频女| 9热在线视频观看99| 亚洲第一区二区三区不卡| 日本一区二区免费在线视频| 丝瓜视频免费看黄片| 中文字幕色久视频| 18禁国产床啪视频网站| 中文字幕高清在线视频| 永久免费av网站大全| 久久人人爽人人片av| 亚洲熟女精品中文字幕| 丝袜喷水一区| 久久影院123| 性高湖久久久久久久久免费观看| 老司机亚洲免费影院| 婷婷色麻豆天堂久久| 国产成人欧美| 国产97色在线日韩免费| 久久天堂一区二区三区四区| 毛片一级片免费看久久久久| 欧美在线一区亚洲| 亚洲av综合色区一区| 日本wwww免费看| 亚洲国产看品久久| 亚洲精品久久午夜乱码| 丝袜美腿诱惑在线| 一边摸一边做爽爽视频免费| 成人国产麻豆网| 日本欧美视频一区| 91老司机精品| 最近中文字幕2019免费版| 国产精品一区二区精品视频观看| 99香蕉大伊视频| 国产成人系列免费观看| 亚洲情色 制服丝袜| 免费高清在线观看日韩| 国产精品免费视频内射| 国产国语露脸激情在线看| 午夜影院在线不卡| 青青草视频在线视频观看| 欧美日韩综合久久久久久| 国产精品99久久99久久久不卡 | 久久久久久免费高清国产稀缺| 亚洲国产av新网站| 一区二区日韩欧美中文字幕| 日日摸夜夜添夜夜爱| 大话2 男鬼变身卡| 婷婷成人精品国产| 国产一区亚洲一区在线观看| 国产精品二区激情视频| 中文字幕精品免费在线观看视频| 亚洲av国产av综合av卡| 国产在线一区二区三区精| 一边亲一边摸免费视频| 欧美激情极品国产一区二区三区| 国产片特级美女逼逼视频| 国产 一区精品| 欧美av亚洲av综合av国产av | 久久国产精品大桥未久av| 一边摸一边抽搐一进一出视频| 男人操女人黄网站| 一级爰片在线观看| 日本爱情动作片www.在线观看| 欧美日韩亚洲综合一区二区三区_| 一级黄片播放器| 视频区图区小说| 在线观看免费视频网站a站| 亚洲,一卡二卡三卡| 国产免费福利视频在线观看| 免费看av在线观看网站| 成年女人毛片免费观看观看9 | 精品国产乱码久久久久久男人| 秋霞在线观看毛片| av片东京热男人的天堂| 久久亚洲国产成人精品v| 日韩欧美精品免费久久| 久久99精品国语久久久| 无遮挡黄片免费观看| 丰满迷人的少妇在线观看| 国产精品av久久久久免费| 亚洲精品第二区| 日本vs欧美在线观看视频| 国产又色又爽无遮挡免| 亚洲国产最新在线播放| 亚洲精品,欧美精品| 免费黄频网站在线观看国产| 91老司机精品| 在线看a的网站| 久久久久久人人人人人| 亚洲欧洲日产国产| 国产精品.久久久| 国产黄色视频一区二区在线观看| 亚洲国产av影院在线观看| 热99久久久久精品小说推荐| 精品少妇内射三级| 香蕉丝袜av| 国产精品久久久久久精品电影小说| 国产高清不卡午夜福利| 一区二区三区精品91| 色婷婷久久久亚洲欧美| 国产福利在线免费观看视频| 日本欧美国产在线视频| 亚洲七黄色美女视频| 在线免费观看不下载黄p国产| 天天操日日干夜夜撸| 欧美日韩亚洲高清精品| 国产亚洲午夜精品一区二区久久| 亚洲人成77777在线视频| 麻豆乱淫一区二区| 狂野欧美激情性xxxx| 捣出白浆h1v1| 一边摸一边抽搐一进一出视频| bbb黄色大片| 午夜免费男女啪啪视频观看| 欧美乱码精品一区二区三区| 日本av免费视频播放| 亚洲国产看品久久| 91精品伊人久久大香线蕉| 这个男人来自地球电影免费观看 | 亚洲美女视频黄频| 国产一卡二卡三卡精品 | 久久久久久久久久久免费av| 亚洲欧美日韩另类电影网站| 一本—道久久a久久精品蜜桃钙片| 亚洲精品美女久久久久99蜜臀 | 国产不卡av网站在线观看| 亚洲欧美一区二区三区久久| 久久性视频一级片| 青春草视频在线免费观看| 亚洲成人手机| 国产精品偷伦视频观看了| 国产xxxxx性猛交| 久久韩国三级中文字幕| 亚洲人成电影观看| 国产免费又黄又爽又色| 精品少妇内射三级| 街头女战士在线观看网站| 久久狼人影院| 黄色怎么调成土黄色| 男的添女的下面高潮视频| 亚洲少妇的诱惑av| xxxhd国产人妻xxx| 久久天躁狠狠躁夜夜2o2o | 高清黄色对白视频在线免费看| 免费黄网站久久成人精品| 久久影院123| 国产成人精品久久二区二区91 | av免费观看日本| 欧美变态另类bdsm刘玥| 我要看黄色一级片免费的| 美女主播在线视频| 国产不卡av网站在线观看| 国产精品嫩草影院av在线观看| 一区二区av电影网| 欧美人与善性xxx| 丝袜脚勾引网站| 亚洲免费av在线视频| 久久鲁丝午夜福利片| 精品少妇一区二区三区视频日本电影 | 又粗又硬又长又爽又黄的视频| 久久精品人人爽人人爽视色| 国产老妇伦熟女老妇高清| 久久99精品国语久久久| 国产深夜福利视频在线观看| 亚洲欧美成人精品一区二区|