• <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)最大殺手
    啪啪无遮挡十八禁网站| 噜噜噜噜噜久久久久久91| 欧美zozozo另类| e午夜精品久久久久久久| 精品一区二区三区av网在线观看| 国产伦一二天堂av在线观看| 伊人久久大香线蕉亚洲五| 亚洲国产欧洲综合997久久,| 亚洲欧美日韩高清专用| 日韩欧美在线二视频| 国产亚洲精品综合一区在线观看| 国产三级在线视频| 超碰av人人做人人爽久久 | 一卡2卡三卡四卡精品乱码亚洲| 精品久久久久久久人妻蜜臀av| 麻豆久久精品国产亚洲av| 精品国产三级普通话版| 久久天躁狠狠躁夜夜2o2o| 精品日产1卡2卡| 国产亚洲精品av在线| 日韩人妻高清精品专区| 麻豆国产97在线/欧美| 成年人黄色毛片网站| av黄色大香蕉| 亚洲 欧美 日韩 在线 免费| 国产欧美日韩精品亚洲av| tocl精华| 国产激情偷乱视频一区二区| 人妻夜夜爽99麻豆av| 国产精品免费一区二区三区在线| 久久久久国内视频| 欧美中文综合在线视频| 网址你懂的国产日韩在线| 亚洲国产精品sss在线观看| 听说在线观看完整版免费高清| 国产成人aa在线观看| 一级毛片女人18水好多| 12—13女人毛片做爰片一| 国产成年人精品一区二区| 啪啪无遮挡十八禁网站| 国产 一区 欧美 日韩| 变态另类成人亚洲欧美熟女| 亚洲av免费在线观看| www.www免费av| 丰满人妻熟妇乱又伦精品不卡| 老熟妇乱子伦视频在线观看| 亚洲激情在线av| 18美女黄网站色大片免费观看| 国产精品久久电影中文字幕| 国产视频一区二区在线看| 国产精品久久久人人做人人爽| 一个人观看的视频www高清免费观看| 成熟少妇高潮喷水视频| 欧美日韩综合久久久久久 | av天堂中文字幕网| 99热只有精品国产| 毛片女人毛片| 美女被艹到高潮喷水动态| 两个人的视频大全免费| 黄色视频,在线免费观看| 国产亚洲欧美98| 特级一级黄色大片| 久久亚洲精品不卡| 国产精品 国内视频| 亚洲男人的天堂狠狠| 国产精品国产高清国产av| 日韩 欧美 亚洲 中文字幕| 99国产精品一区二区蜜桃av| 欧美成人免费av一区二区三区| 伊人久久精品亚洲午夜| 久久精品国产自在天天线| av视频在线观看入口| 色综合欧美亚洲国产小说| 757午夜福利合集在线观看| 少妇丰满av| 母亲3免费完整高清在线观看| 校园春色视频在线观看| av天堂中文字幕网| 少妇的逼好多水| 亚洲国产色片| netflix在线观看网站| av天堂在线播放| 久久久精品欧美日韩精品| 国内精品久久久久久久电影| 久久久久免费精品人妻一区二区| 搡女人真爽免费视频火全软件 | 搞女人的毛片| 国产真人三级小视频在线观看| 狠狠狠狠99中文字幕| 欧美成人免费av一区二区三区| 欧美不卡视频在线免费观看| 亚洲av免费在线观看| 久9热在线精品视频| 国产伦精品一区二区三区四那| 大型黄色视频在线免费观看| 18禁美女被吸乳视频| 欧美日韩亚洲国产一区二区在线观看| 麻豆成人av在线观看| 18禁国产床啪视频网站| 麻豆国产av国片精品| 丁香六月欧美| 成人三级黄色视频| 国产成人影院久久av| 国产高清三级在线| 一区二区三区高清视频在线| 在线观看美女被高潮喷水网站 | 看免费av毛片| 亚洲精品成人久久久久久| 久久亚洲真实| 国产av麻豆久久久久久久| 日本五十路高清| 免费观看人在逋| 欧美极品一区二区三区四区| 成人特级黄色片久久久久久久| 国产高清视频在线观看网站| 亚洲欧美激情综合另类| 手机成人av网站| 久久精品国产清高在天天线| 精品无人区乱码1区二区| 亚洲一区二区三区不卡视频| 丰满人妻一区二区三区视频av | 熟女电影av网| 99在线人妻在线中文字幕| 老司机午夜十八禁免费视频| 亚洲精品一卡2卡三卡4卡5卡| 国产精品一及| 国产日本99.免费观看| 久久精品国产亚洲av涩爱 | xxxwww97欧美| 久久久久国产精品人妻aⅴ院| 国产伦精品一区二区三区视频9 | 婷婷精品国产亚洲av| 国产91精品成人一区二区三区| 亚洲成av人片免费观看| 特级一级黄色大片| 日本精品一区二区三区蜜桃| 国产精品美女特级片免费视频播放器| 欧美绝顶高潮抽搐喷水| 丰满人妻一区二区三区视频av | 精品国产亚洲在线| 亚洲性夜色夜夜综合| 婷婷丁香在线五月| 亚洲成人免费电影在线观看| 99国产综合亚洲精品| 亚洲avbb在线观看| 99久久久亚洲精品蜜臀av| 欧美一区二区亚洲| 国产综合懂色| 2021天堂中文幕一二区在线观| 亚洲aⅴ乱码一区二区在线播放| 99久久成人亚洲精品观看| 身体一侧抽搐| 亚洲黑人精品在线| 真实男女啪啪啪动态图| 亚洲专区中文字幕在线| 亚洲aⅴ乱码一区二区在线播放| 在线播放国产精品三级| 国产精品综合久久久久久久免费| 啪啪无遮挡十八禁网站| 看黄色毛片网站| 日韩成人在线观看一区二区三区| 久久精品91无色码中文字幕| 欧美日本视频| 精品乱码久久久久久99久播| 国产视频内射| 小说图片视频综合网站| 99国产综合亚洲精品| 日本在线视频免费播放| 国产男靠女视频免费网站| 69av精品久久久久久| 观看免费一级毛片| 97人妻精品一区二区三区麻豆| 18美女黄网站色大片免费观看| 国产欧美日韩精品一区二区| 热99在线观看视频| 在线免费观看不下载黄p国产 | 国产亚洲精品av在线| 热99在线观看视频| 夜夜夜夜夜久久久久| 欧美最新免费一区二区三区 | 亚洲国产精品999在线| 免费搜索国产男女视频| www.熟女人妻精品国产| 国产精品久久久久久久电影 | 国产乱人伦免费视频| 亚洲熟妇中文字幕五十中出| 亚洲人与动物交配视频| 国产综合懂色| 91在线观看av| 级片在线观看| 麻豆一二三区av精品| 国产乱人视频| 高清日韩中文字幕在线| 国产黄色小视频在线观看| 在线观看日韩欧美| 法律面前人人平等表现在哪些方面| 国产色婷婷99| 久久香蕉精品热| 黄色女人牲交| 久久6这里有精品| 久久久国产成人精品二区| 看黄色毛片网站| 一进一出抽搐gif免费好疼| 久久中文看片网| 一个人观看的视频www高清免费观看| 欧美黄色淫秽网站| 午夜视频国产福利| 三级国产精品欧美在线观看| 久久婷婷人人爽人人干人人爱| 亚洲欧美一区二区三区黑人| 91麻豆av在线| 亚洲专区国产一区二区| 天美传媒精品一区二区| 国产探花在线观看一区二区| 99久久精品热视频| 很黄的视频免费| 国产老妇女一区| 亚洲性夜色夜夜综合| 久久人妻av系列| 亚洲美女黄片视频| 熟女人妻精品中文字幕| 亚洲电影在线观看av| 超碰av人人做人人爽久久 | 一区二区三区高清视频在线| 午夜免费激情av| 欧美日韩国产亚洲二区| 琪琪午夜伦伦电影理论片6080| 蜜桃久久精品国产亚洲av| 大型黄色视频在线免费观看| 久久久色成人| 变态另类成人亚洲欧美熟女| 欧美性感艳星| 成人欧美大片| 黄片小视频在线播放| 蜜桃亚洲精品一区二区三区| 久久精品人妻少妇| 日韩亚洲欧美综合| 亚洲成人免费电影在线观看| 亚洲五月天丁香| 久久精品国产综合久久久| 在线免费观看的www视频| 欧美中文日本在线观看视频| 成人国产一区最新在线观看| 免费看十八禁软件| 亚洲欧美日韩东京热| 老汉色av国产亚洲站长工具| 久久人妻av系列| 啦啦啦免费观看视频1| 国产真实乱freesex| 波野结衣二区三区在线 | 怎么达到女性高潮| 国产精品一及| 日本在线视频免费播放| 村上凉子中文字幕在线| 在线看三级毛片| 欧美日韩黄片免| 色播亚洲综合网| 人人妻人人看人人澡| 一进一出好大好爽视频| 日本三级黄在线观看| 国产野战对白在线观看| 国产成人av教育| 全区人妻精品视频| 国模一区二区三区四区视频| 网址你懂的国产日韩在线| 黄色视频,在线免费观看| 国产色婷婷99| 深夜精品福利| 天堂影院成人在线观看| 亚洲片人在线观看| 蜜桃亚洲精品一区二区三区| 村上凉子中文字幕在线| 亚洲欧美激情综合另类| 亚洲av成人不卡在线观看播放网| 国产淫片久久久久久久久 | 夜夜爽天天搞| 亚洲欧美日韩东京热| 九九热线精品视视频播放| 久久人妻av系列| 亚洲片人在线观看| 欧美日韩亚洲国产一区二区在线观看| 成人三级黄色视频| 性欧美人与动物交配| 91字幕亚洲| 禁无遮挡网站| 一区二区三区激情视频| 最后的刺客免费高清国语| 免费看光身美女| 综合色av麻豆| 免费一级毛片在线播放高清视频| 中文字幕人妻丝袜一区二区| 精品久久久久久久久久免费视频| 久久久久久久久大av| 色综合欧美亚洲国产小说| av中文乱码字幕在线| 1024手机看黄色片| 老汉色∧v一级毛片| 老鸭窝网址在线观看| 国产亚洲精品一区二区www| 国产精品三级大全| av专区在线播放| 免费电影在线观看免费观看| 高清在线国产一区| 搡女人真爽免费视频火全软件 | 亚洲精品在线观看二区| 黄色片一级片一级黄色片| 亚洲 欧美 日韩 在线 免费| 亚洲五月婷婷丁香| 九九在线视频观看精品| 十八禁人妻一区二区| h日本视频在线播放| 精品免费久久久久久久清纯| 国产免费男女视频| 精品免费久久久久久久清纯| 伊人久久精品亚洲午夜| 亚洲七黄色美女视频| 亚洲国产日韩欧美精品在线观看 | 久久久久国产精品人妻aⅴ院| 国内毛片毛片毛片毛片毛片| 国产熟女xx| 综合色av麻豆| netflix在线观看网站| 亚洲精品国产精品久久久不卡| 精品久久久久久久久久免费视频| 亚洲18禁久久av| 中文在线观看免费www的网站| 精品国产三级普通话版| 久久伊人香网站| 午夜福利欧美成人| 欧美日韩亚洲国产一区二区在线观看| 最新中文字幕久久久久| 亚洲avbb在线观看| 草草在线视频免费看| 99久久综合精品五月天人人| 亚洲精品成人久久久久久| 亚洲一区高清亚洲精品| 精品久久久久久,| 在线看三级毛片| 国产精品综合久久久久久久免费| 亚洲成人久久爱视频| 日韩欧美 国产精品| 最后的刺客免费高清国语| 亚洲精品一卡2卡三卡4卡5卡| 欧美在线一区亚洲| 亚洲av五月六月丁香网| 亚洲五月天丁香| 老司机在亚洲福利影院| 日本 欧美在线| 国产乱人伦免费视频| 亚洲国产精品成人综合色| 免费看光身美女| 国产伦人伦偷精品视频| 国产精品乱码一区二三区的特点| 最近最新中文字幕大全免费视频| 亚洲国产精品合色在线| 老鸭窝网址在线观看| 亚洲色图av天堂| 99国产精品一区二区蜜桃av| 亚洲精品一卡2卡三卡4卡5卡| 夜夜爽天天搞| 在线观看免费视频日本深夜| 免费观看的影片在线观看| 全区人妻精品视频| 真人一进一出gif抽搐免费| 亚洲电影在线观看av| 毛片女人毛片| 日韩欧美 国产精品| 麻豆成人午夜福利视频| 国产高清视频在线播放一区| 非洲黑人性xxxx精品又粗又长| 亚洲专区中文字幕在线| 97超级碰碰碰精品色视频在线观看| 日韩精品青青久久久久久| 一进一出抽搐动态| 久久欧美精品欧美久久欧美| 美女高潮的动态| 日韩欧美精品v在线| 欧美3d第一页| 午夜福利欧美成人| 在线播放国产精品三级| 男女床上黄色一级片免费看| 人妻丰满熟妇av一区二区三区| 五月玫瑰六月丁香| 国产精品电影一区二区三区| 亚洲狠狠婷婷综合久久图片| 国产精品99久久久久久久久| 色播亚洲综合网| 黄色成人免费大全| 久久久久九九精品影院| 国产一区二区三区视频了| 久9热在线精品视频| 午夜福利高清视频| 在线观看免费午夜福利视频| 久久欧美精品欧美久久欧美| 免费看日本二区| 成人三级黄色视频| 国产久久久一区二区三区| 搡女人真爽免费视频火全软件 | 少妇丰满av| 亚洲人成网站在线播| 久久精品夜夜夜夜夜久久蜜豆| 国产精品永久免费网站| 岛国视频午夜一区免费看| 国产色婷婷99| 国产亚洲欧美在线一区二区| 亚洲一区二区三区不卡视频| 亚洲国产精品成人综合色| 黄色成人免费大全| tocl精华| av天堂在线播放| www国产在线视频色| 一a级毛片在线观看| 国产久久久一区二区三区| 久久精品国产清高在天天线| 欧美最黄视频在线播放免费| av在线天堂中文字幕| 91在线观看av| 香蕉丝袜av| 国产毛片a区久久久久| 欧美成人a在线观看| 欧美不卡视频在线免费观看| 男人的好看免费观看在线视频| 精品国产超薄肉色丝袜足j| 亚洲国产色片| 久久久久亚洲av毛片大全| 最新在线观看一区二区三区| 国产不卡一卡二| 高清在线国产一区| 国产国拍精品亚洲av在线观看 | 一二三四社区在线视频社区8| 看片在线看免费视频| 偷拍熟女少妇极品色| 亚洲欧美日韩东京热| 日本 欧美在线| 亚洲 欧美 日韩 在线 免费| 51国产日韩欧美| 女人高潮潮喷娇喘18禁视频| 国产精品久久电影中文字幕| 国产高清视频在线观看网站| 久久精品国产综合久久久| 91在线精品国自产拍蜜月 | 在线观看日韩欧美| 色在线成人网| 亚洲国产精品久久男人天堂| 亚洲av成人av| 国产伦一二天堂av在线观看| 午夜福利在线观看吧| 精品一区二区三区视频在线观看免费| 热99re8久久精品国产| 少妇人妻一区二区三区视频| 久久伊人香网站| 国产精品一区二区三区四区久久| 久久精品国产综合久久久| 亚洲av熟女| 成人特级黄色片久久久久久久| 99精品久久久久人妻精品| 国产三级在线视频| 丰满人妻熟妇乱又伦精品不卡| 亚洲精品色激情综合| 色综合婷婷激情| av视频在线观看入口| 国产精品久久久久久久电影 | 3wmmmm亚洲av在线观看| 丰满人妻熟妇乱又伦精品不卡| 亚洲av成人av| 波多野结衣高清无吗| 久久精品国产清高在天天线| 精品久久久久久久人妻蜜臀av| 欧美性感艳星| 国产精品一区二区三区四区免费观看 | 国产麻豆成人av免费视频| 色在线成人网| 午夜久久久久精精品| 国产三级黄色录像| 欧美日韩瑟瑟在线播放| 成年人黄色毛片网站| 日本五十路高清| 亚洲精品乱码久久久v下载方式 | 熟妇人妻久久中文字幕3abv| 欧美xxxx黑人xx丫x性爽| 日韩欧美 国产精品| 久久久国产成人精品二区| 他把我摸到了高潮在线观看| 窝窝影院91人妻| 欧美一区二区亚洲| 欧美日韩综合久久久久久 | 国产精品自产拍在线观看55亚洲| 搞女人的毛片| 国产亚洲av嫩草精品影院| 99国产极品粉嫩在线观看| 欧美中文综合在线视频| 久久香蕉国产精品| 色老头精品视频在线观看| 久久午夜亚洲精品久久| 无遮挡黄片免费观看| 国产高清videossex| www.999成人在线观看| 91在线观看av| 日韩高清综合在线| 波多野结衣高清无吗| 亚洲国产中文字幕在线视频| 有码 亚洲区| 一个人看视频在线观看www免费 | 亚洲成人精品中文字幕电影| 最近最新中文字幕大全电影3| 国产亚洲欧美在线一区二区| 白带黄色成豆腐渣| 国产精品美女特级片免费视频播放器| 色噜噜av男人的天堂激情| 精品久久久久久久久久免费视频| 国产精品女同一区二区软件 | 又黄又爽又免费观看的视频| 麻豆国产av国片精品| 91av网一区二区| 国产成人av激情在线播放| 叶爱在线成人免费视频播放| 精品不卡国产一区二区三区| 在线看三级毛片| 他把我摸到了高潮在线观看| netflix在线观看网站| 一区二区三区免费毛片| ponron亚洲| 成人亚洲精品av一区二区| 中文字幕人妻丝袜一区二区| 免费看光身美女| 精品欧美国产一区二区三| 美女大奶头视频| 国产色爽女视频免费观看| 极品教师在线免费播放| 久久香蕉精品热| av片东京热男人的天堂| 97超视频在线观看视频| 国产精品亚洲美女久久久| 美女cb高潮喷水在线观看| 村上凉子中文字幕在线| 精品久久久久久久久久免费视频| 性色av乱码一区二区三区2| 亚洲中文字幕一区二区三区有码在线看| 女人被狂操c到高潮| 亚洲精品色激情综合| 国产精品99久久久久久久久| 久久午夜亚洲精品久久| 精华霜和精华液先用哪个| 国产色婷婷99| 亚洲国产精品成人综合色| 亚洲专区中文字幕在线| 高清在线国产一区| 亚洲性夜色夜夜综合| 国产精品久久久久久久电影 | 宅男免费午夜| 特级一级黄色大片| 亚洲欧美日韩高清专用| 少妇丰满av| 成人精品一区二区免费| 日韩欧美精品免费久久 | 午夜亚洲福利在线播放| 我的老师免费观看完整版| 欧美日韩亚洲国产一区二区在线观看| 国产一区二区亚洲精品在线观看| 波野结衣二区三区在线 | 每晚都被弄得嗷嗷叫到高潮| 精品久久久久久成人av| 一个人免费在线观看的高清视频| 全区人妻精品视频| 国产v大片淫在线免费观看| 亚洲内射少妇av| 老熟妇乱子伦视频在线观看| 国产中年淑女户外野战色| 久久中文看片网| 99国产精品一区二区三区| 成人av在线播放网站| 国产高清激情床上av| 亚洲成av人片在线播放无| 国产精华一区二区三区| 美女大奶头视频| 国产激情欧美一区二区| 啪啪无遮挡十八禁网站| 精品人妻1区二区| 国产成人a区在线观看| 成人午夜高清在线视频| 97碰自拍视频| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 日本a在线网址| 一进一出抽搐gif免费好疼| 亚洲第一电影网av| 深爱激情五月婷婷| 精品国产亚洲在线| 成人亚洲精品av一区二区| 午夜福利免费观看在线| 亚洲电影在线观看av| 成人亚洲精品av一区二区| 露出奶头的视频| 亚洲成人免费电影在线观看| 特大巨黑吊av在线直播| 欧美性猛交黑人性爽| 亚洲成人免费电影在线观看| 精品午夜福利视频在线观看一区| 国产一区二区在线av高清观看| 变态另类成人亚洲欧美熟女| 国产成人欧美在线观看| 床上黄色一级片| 在线观看免费视频日本深夜| 国产v大片淫在线免费观看| 久久久久久久午夜电影| 在线天堂最新版资源| 中出人妻视频一区二区| 亚洲精品国产精品久久久不卡| 99视频精品全部免费 在线| 久久精品国产99精品国产亚洲性色| 长腿黑丝高跟| 小说图片视频综合网站| 亚洲欧美精品综合久久99| or卡值多少钱| 免费在线观看影片大全网站| 级片在线观看|