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

    張掖濕地甲烷通量動(dòng)態(tài)特征及其影響因子

    2017-10-13 12:28:14蔣國(guó)慶徐自為劉紹民
    生態(tài)學(xué)報(bào) 2017年17期
    關(guān)鍵詞:張掖土壤溫度通量

    張 強(qiáng),蔣國(guó)慶,孫 睿,*,徐自為,劉紹民

    1 遙感科學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室(北京師范大學(xué),中國(guó)科學(xué)院遙感應(yīng)用研究所),北京 100875 2 北京師范大學(xué)地理學(xué)與遙感科學(xué)學(xué)院,北京 100875 3 環(huán)境遙感與數(shù)字城市北京市重點(diǎn)實(shí)驗(yàn)室,北京 100875

    張掖濕地甲烷通量動(dòng)態(tài)特征及其影響因子

    張 強(qiáng)1,2,3,蔣國(guó)慶1,2,3,孫 睿1,2,3,*,徐自為1,2,劉紹民1,2

    1 遙感科學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室(北京師范大學(xué),中國(guó)科學(xué)院遙感應(yīng)用研究所),北京 100875 2 北京師范大學(xué)地理學(xué)與遙感科學(xué)學(xué)院,北京 100875 3 環(huán)境遙感與數(shù)字城市北京市重點(diǎn)實(shí)驗(yàn)室,北京 100875

    于2012年7月—2014年6月對(duì)地處干旱區(qū)的張掖濕地甲烷(CH4)通量進(jìn)行觀測(cè),分析其CH4通量的變化特征及其影響因子。結(jié)果表明:CH4通量的日變化趨勢(shì)總體表現(xiàn)為白天大于夜間;不同季節(jié)CH4通量排放特征差異明顯,夏季最大,春秋次之,冬季最??;CH4通量日總量與空氣溫度、土壤溫度之間指數(shù)相關(guān)關(guān)系顯著,其中4 cm處土壤溫度與之相關(guān)性最強(qiáng);1—6月摩擦風(fēng)速(U*)與CH4通量顯著正相關(guān);結(jié)合CO2通量觀測(cè)數(shù)據(jù),研究時(shí)段張掖濕地凈碳吸收量為495.92 g C m-2a-1,為明顯碳匯。

    濕地;渦動(dòng)相關(guān);CH4通量;溫度;摩擦風(fēng)速

    Abstract: Analyzing the characteristics and variations of CH4flux in wetland would be very helpful in order to understand the importance of the carbon cycle in the terrestrial ecosystem. In the arid area, the Zhangye wetland features complex meteorological conditions, which have introduced obvious uncertainties into the carbon fluxes. However, there are few research studies on the carbon flux of wetland in the arid areas. It is helpful to strengthen the study of wetland ecological system in the arid areas. In order to study the dynamics of carbon flux in a wetland ecosystem, CO2flux and CH4flux were measured using an eddy covariance technique over a phragmites-dominated wetland in Zhangye, Gansu Province (from June in 2012 to August in 2014). In the present study, the CH4flux variations and its main driving factors were analyzed based on the eddy covariance observations. The raw 10 Hz data were processed into 30 min mean values based on data quality control (spike detection, coordinate rotation, frequency response correction, etc.) and gap-filling (linear interpolation and Look-up Table method). The result showed the CH4flux reaching its maximum in summer (0.34 g CH4m-2d-1), then spring and autumn, while getting the minimum values in winter (0.02 g CH4m-2d-1). Diurnal variation of CH4flux indicated that the daytime values were higher than that of nighttime. There was an exponential relationship between CH4flux and air temperature and soil temperature, and it was mostly correlated with the soil temperature in 4-cm depth, where the relationship wasy=0.07614e0.03402xandR2is 0.3560. CH4fluxes also exhibited a strong relationship with friction velocity (R2=0.91 and 0.94, respectively, from January to March and from April to June). The enhanced turbulence could increase the ventilation and exchange of CH4flux. CH4flux significantly influenced the total carbon budget. By considering the absorption and emission of CO2and CH4fluxes, the net carbon uptake was 495.92 g C m-2a-1. The carbon sink in Zhangye wetland was apparent during the study period.

    KeyWords: wetland; eddy covariance; CH4flux; temperature; friction velocity

    目前,國(guó)內(nèi)外學(xué)者已從多個(gè)方面對(duì)濕地CH4變化特性和影響因素進(jìn)行過(guò)研究。不同環(huán)境和植被類型,濕地CH4通量有明顯差異。黃國(guó)宏等[5]用封閉式箱法測(cè)定的遼河三角洲蘆葦濕地CH4觀測(cè)結(jié)果表明,其排放有明顯的季節(jié)變化規(guī)律,淹水前土壤為CH4匯,淹水期間有大量的CH4排放,排水后CH4排放明顯減少。Kang等[6]對(duì)英國(guó)北威爾士酸性泥炭沼澤和森林沼澤CH4的測(cè)定,發(fā)現(xiàn)不同類型濕地CH4排放存在差異:酸性泥炭沼澤CH4排放高峰在春季和秋季,而在森林沼澤,其釋放峰值出現(xiàn)在10月和11月。Rinne[7]等利用渦動(dòng)相關(guān)技術(shù)觀測(cè)北方沼澤濕地CH4排放,發(fā)現(xiàn)無(wú)明顯的日變化規(guī)律,非雪期CH4排放量占全年的91%。Long[8]等分析了加拿大北部CH4的日變化和季節(jié)變化規(guī)律,CH4排放白天明顯大于夜間,且5月從0附近開(kāi)始增大,7月底達(dá)到峰值,8月CH4排放開(kāi)始變小。

    張掖濕地地處河西走廊中部、黑河沖積扇形成的三角洲之上,是由河流、草本沼澤、濕草甸等天然濕地,以及人工湖、池塘、溝渠等人工濕地為主體構(gòu)成的復(fù)合濕地生態(tài)系統(tǒng),濕地類型多樣,原生態(tài)特征突出[9]。目前關(guān)于干旱區(qū)濕地CH4通量的研究很少,對(duì)張掖濕地CH4通量的觀測(cè)研究有助于深入分析干旱區(qū)綠洲生態(tài)系統(tǒng)碳循環(huán)特征及變化規(guī)律,并有利于干旱區(qū)綠洲生態(tài)系統(tǒng)的保護(hù)。本文基于2012年7月至2014年6月張掖濕地生態(tài)系統(tǒng)野外觀測(cè)資料,擬解決的關(guān)鍵問(wèn)題如下:(1)分析CH4通量變化特征及影響因素,為研究干旱區(qū)CH4通量研究提供科學(xué)依據(jù);(2)結(jié)合CO2與CH4的吸收和排放,準(zhǔn)確評(píng)價(jià)干旱區(qū)濕地生態(tài)系統(tǒng)碳源匯功能。

    1 材料與方法

    1.1 研究區(qū)概況

    1.1.1 地理位置

    張掖市位于黑河流域中游,其濕地共有天然濕地、人工濕地2個(gè)大類4個(gè)類型13個(gè)類別,總面積為210420.42 hm2,占全市國(guó)土面積的5.02%。其中天然濕地有河流濕地、湖泊濕地和沼澤濕地3個(gè)類型,9個(gè)類別,面積為199709.97 hm2,占全市濕地總面積的94.90%;人工濕地有水產(chǎn)池塘、灌溉地、蓄水區(qū)、鹽田4個(gè)類型,面積為10710.45 hm2,占全市濕地總面積的5.10%[9]。通量觀測(cè)站位于甘肅省張掖市國(guó)家濕地公園內(nèi)。張掖國(guó)家濕地公園位于張掖城區(qū)北郊(圖1),總面積4108 hm2,其中濕地面積1733 hm2,多以蘆葦濕地為主。

    圖1 張掖濕地通量觀測(cè)站地理位置(100.44640°E, 38.97514°N)Fig.1 Position of eddy covariance tower of Zhangye wetland (100.44640°E, 38.97514°N)紅色十字標(biāo)記處為通量塔位置

    1.1.2 氣候條件及植被長(zhǎng)勢(shì)分析

    張掖地區(qū)氣候?qū)倜黠@的溫帶大陸性氣候,其顯著特點(diǎn)是:降水稀少且集中,年內(nèi)降水分布不均,年際變化較大;蒸發(fā)強(qiáng)烈,年平均蒸發(fā)量2047 mm。多大風(fēng)、風(fēng)沙天氣,風(fēng)向以北風(fēng)、西風(fēng)為主;年平均氣溫6℃,7月份氣溫最高,1月份氣溫最低。

    經(jīng)地面調(diào)查,濕地水位在4月和10月最高,4月水位高是由于冰雪融化,10月則是附近農(nóng)田停止灌溉,耗水量降低,使?jié)竦厮簧?。此?由于政策原因,在“均水”期水位也較高,即黑河上游全線禁止灌溉,直接供給內(nèi)蒙古額濟(jì)納旗。一年內(nèi)“均水”3次,分別在6、7和9月,每次10—15 d。

    葉面積指數(shù)(LAI)是植被冠層結(jié)構(gòu)的重要參數(shù)之一,對(duì)于了解植被生長(zhǎng)發(fā)育過(guò)程、植被長(zhǎng)勢(shì)等具有重要意義。本研究選擇觀測(cè)期間MODIS LAI產(chǎn)品(http://modis.gsfc.nasa.gov/)對(duì)張掖濕地植被長(zhǎng)勢(shì)進(jìn)行分析(圖2)。生長(zhǎng)季(5—10月)LAI均表現(xiàn)出先增大后減小的趨勢(shì),非生長(zhǎng)季LAI較小。不同年份LAI峰值有差異, LAI峰值差異體現(xiàn)了植被長(zhǎng)勢(shì)的差異。

    圖2 研究時(shí)段MODIS LAI變化趨勢(shì)Fig.2 MODIS LAI during the study period in Zhangye wetland

    1.2 研究數(shù)據(jù)

    張掖濕地通量觀測(cè)站地理坐標(biāo)為100.44640°E,38.97514°N,海拔1460.00 m,是“黑河流域生態(tài)-水文過(guò)程綜合遙感觀測(cè)試驗(yàn):水文氣象要素與多尺度蒸散發(fā)觀測(cè)試驗(yàn)”(Heihe Watershed Allied Telemetry Experiment Research-the Multi-Scale Observation Experiment on Evapotranspiration over heterogeneous land surfaces, HiWATER-MUSOEXE)[10- 12]觀測(cè)站點(diǎn)之一,自2012年6月以來(lái),采用渦動(dòng)相關(guān)系統(tǒng)對(duì)濕地生態(tài)系統(tǒng)碳通量進(jìn)行連續(xù)測(cè)量。

    渦動(dòng)相關(guān)儀架設(shè)高度為5.2 m。渦動(dòng)相關(guān)系統(tǒng)主要由三維超聲風(fēng)速儀(CSAT3, Campbell Scientific, USA)和開(kāi)路CH4/H2O 紅外氣體分析儀(Li- 7500A, Licor Inc., USA)組成,原始數(shù)據(jù)采樣頻率為10 Hz,每30 min輸出一組通量平均值以及10 Hz原始數(shù)據(jù)。超聲風(fēng)速儀朝向正北,與CH4/H2O紅外氣體分析儀間距為25 cm。通量塔上配有風(fēng)、溫、濕梯度儀以及四分量輻射分析儀,另外,對(duì)塔下不同深度的土壤溫度也進(jìn)行測(cè)定。

    1.3 數(shù)據(jù)處理

    1.3.1 數(shù)據(jù)質(zhì)量控制

    在對(duì)CH4通量數(shù)據(jù)處理時(shí),基本處理步驟參考已較成熟的CO2通量處理方法[8]。采用Eddypro軟件(http://www.licor.com/env/products/eddy_covariance/software.html)對(duì)渦動(dòng)相關(guān)原始10 Hz數(shù)據(jù)進(jìn)行處理,主要步驟包括:野點(diǎn)值剔除、延遲時(shí)間校正、坐標(biāo)旋轉(zhuǎn)(二次坐標(biāo)旋轉(zhuǎn))、頻率響應(yīng)修正、角度訂正、超聲虛溫和密度(WPL)修正等[13]。由于大氣穩(wěn)定度、天氣原因以及儀器的物理限制等影響,Eddypro軟件輸出的30 min通量值仍然存在一些野點(diǎn),需要進(jìn)一步對(duì)通量數(shù)據(jù)進(jìn)行質(zhì)量控制,剔除如下數(shù)據(jù)[14]:(1)降雨前后1 h數(shù)據(jù);(2)超出儀器測(cè)量量程或者合理范圍的數(shù)據(jù);(3)夜間摩擦風(fēng)速(U*)低于臨界值的數(shù)據(jù);(4)CH4通量中的負(fù)值[15]。共剔除約60%的CH4通量無(wú)效數(shù)據(jù)。

    1.3.2 數(shù)據(jù)插補(bǔ)

    CH4通量數(shù)據(jù)插補(bǔ)時(shí),簡(jiǎn)單建立有效數(shù)據(jù)與溫度或輻射的相關(guān)關(guān)系均不理想,考慮到影響CH4通量排放的因素及現(xiàn)有數(shù)據(jù),最終利用土壤溫度和輻射建立查找表分別對(duì)不同階段的CH4通量數(shù)據(jù)進(jìn)行插補(bǔ)。在建立查找表時(shí),由于CH4通量的量級(jí)較小,微小的溫度或輻射變化均會(huì)對(duì)插補(bǔ)結(jié)果產(chǎn)生較大影響。因此,查找表中土壤溫度和輻射的劃分等級(jí)相對(duì)于CO2通量更為精細(xì)。利用查找表仍未插補(bǔ)的數(shù)據(jù)用臨近值的線性平均值來(lái)代替。

    2 結(jié)果

    2.1濕地CH4通量動(dòng)態(tài)特征

    本研究選擇2013年4個(gè)不同時(shí)期CH4通量日均值分析其日變化規(guī)律:(1)3月1日至3月10日,春季,生長(zhǎng)季初期;(2)7月21日至7月30日,夏季,生長(zhǎng)旺季;(3)11月11日至11月20日,秋季,生長(zhǎng)季末期;(4)1月21日至1月30日,冬季,非生長(zhǎng)季。

    從圖3可以看出,CH4通量范圍在0—0.5 μmol m-2s-1之間,整體排放量級(jí)明顯小于CO2通量。日變化趨勢(shì)總體表現(xiàn)為白天通量大于夜間通量,且峰值集中在正午前后。不同季節(jié)CH4通量排放特征差異顯著。CH4通量變化范圍從冬季的0.01 μmol m-2s-1到夏季0.1—0.4 μmol m-2s-1,其余兩時(shí)段CH4通量均在0—0.1 μmol m-2s-1內(nèi)波動(dòng)。

    圖3 2013年張掖濕地不同季節(jié)CH4通量日變化規(guī)律Fig.3 Daily variations in CH4 flux in different seasons of Zhangye wetland

    圖4可以看出,濕地CH4通量有明顯季節(jié)變化規(guī)律。冬季最小,春秋次之,夏季通量值最大。2012年7月至2014年6月平均CH4排放量為37.96 g CH4m-2a-1。

    圖4 CH4季節(jié)變化規(guī)律及4 cm處土壤溫度變化趨勢(shì)Fig.4 Seasonal variations of CH4 flux and soil temperature (4 cm)

    2.2CH4通量的影響因子分析

    2.2.1CH4通量對(duì)溫度的響應(yīng)

    本文分析CH4通量日總量與空氣溫度、土壤溫度之間的相關(guān)關(guān)系。CH4通量日總量與空氣溫度之間指數(shù)關(guān)系顯著,擬合關(guān)系式為y=0.07614e0.03402x,R2=0.3560。對(duì)比不同深度土壤溫度與CH4通量日總量相關(guān)關(guān)系發(fā)現(xiàn)(表1),4 cm處土壤溫度與之相關(guān)性最強(qiáng),即相對(duì)于其它深度的土壤溫度,4 cm處土壤溫度更能代表產(chǎn)甲烷菌活性的最佳溫度條件。

    表1 CH4通量與土壤溫度擬合關(guān)系相關(guān)參數(shù)

    2.2.2CH4通量對(duì)U*的響應(yīng)

    圖5 CH4平均值與U*的相關(guān)關(guān)系Fig.5 Related relationship between CH4 flux and U*

    為了分析CH4通量與大氣湍流之間可能的相關(guān)關(guān)系,采用Long等[8]在2010年提出的方法進(jìn)行分析。將2013年1月至11月的所有通量分為4個(gè)不同時(shí)間段:1—3月、4—6月、7—9月及10—11月。每個(gè)時(shí)間段數(shù)據(jù)按照U*排序,并將數(shù)據(jù)等分為10組,分別計(jì)算每組U*平均值、標(biāo)準(zhǔn)差及CH4通量平均值(圖5)。在1—3月和4—6月期間,CH4通量與U*相關(guān)性較強(qiáng),擬合關(guān)系式分別為y=0.087x+0.019,R2=0.91和y=0.170x+0.074,R2=0.94。7—9月、10—11月,CH4通量與U*相關(guān)性較弱,擬合關(guān)系式分別為y=0.151x+0.113,R2=0.77和y=0.107x+0.080,R2=0.83。

    2.3 濕地碳源匯分析

    本研究選擇2012年7月至2014年6月CH4通量數(shù)據(jù),結(jié)合CO2通量觀測(cè)數(shù)據(jù)[16],以月為步長(zhǎng)分析CO2和CH4的吸收、排放關(guān)系。在圖6中,碳的凈吸收或排放為CO2和CH4的總碳排放與碳吸收之和,正值表示碳排放,負(fù)值表示碳吸收??梢钥闯?3、4、5月由于溫度增加,生態(tài)系統(tǒng)呼吸大幅增加,而受蘆葦生長(zhǎng)的影響,光合作用固定的碳仍較小,加之CH4排放,濕地整體表現(xiàn)為碳源。除3、4、5月外,其他月份碳吸收量均大于排放量,為明顯碳匯。

    研究時(shí)段年均碳吸收量為1287.30 g C m-2a-1,排放量為791.38 g C m-2a-1,即濕地凈吸收碳495.92 g C m-2a-1,為碳匯。

    圖6 2012年7月至2014年6月張掖濕地CO2和CH4月總量變化趨勢(shì)Fig.6 Seasonal variations of CO2 and CH4 flux of Zhangye wetland during the study period

    3 討論

    3.1CH4通量變化特征分析

    3.1.1CH4通量日變化特征

    大量觀測(cè)資料表明,天然濕地CH4通量日變化大致有3種形式:日間極大值型、夜間極大值型和隨機(jī)型[17]。張掖濕地CH4通量日變化特征屬于日間極大值型,不同季節(jié)日變化規(guī)律的差異主要是來(lái)自不同時(shí)期CH4排放和傳輸控制因素的差別。春季和冬季,溫度是制約CH4排放的主要因素;夏秋季由于天氣因素、植被生長(zhǎng)及水位變化等各種因素的綜合影響,使CH4排放的日變化更加復(fù)雜。

    3.1.2CH4通量季節(jié)變化特征

    冬季由于溫度較低,產(chǎn)甲烷菌活性降低,CH4通量接近于零(0.02 g CH4m-2d-1)。春季隨著溫度升高,產(chǎn)甲烷菌活性增強(qiáng),CH4通量開(kāi)始增大。同時(shí)3月底4月初有一 明顯的CH4排放峰值。分析4 cm處土壤溫度數(shù)據(jù)可知(圖4),該時(shí)間段溫度開(kāi)始回升至零度以上,結(jié)合CH4產(chǎn)生的過(guò)程,可推測(cè)是由于冬季結(jié)冰,造成CH4傳輸不暢,產(chǎn)生的CH4在冰面下累積,到春季冰面消融時(shí),阻礙消失,累積的CH4排放到大氣中[18-19]。2013年春季CH4通量值比2014年高的原因主要是由于濕地公園通量站附近開(kāi)挖河道,一方面增加了水體覆蓋的面積,擴(kuò)大了CH4產(chǎn)生的環(huán)境;另一方面挖取土壤有利于深層CH4的排放。夏季由于植被生長(zhǎng)旺盛,有足夠的碳有機(jī)質(zhì)底物,并且溫度較高,產(chǎn)甲烷菌的活性增大,因此CH4通量值最大(0.34 g CH4m-2d-1)。秋季溫度降低,CH4通量開(kāi)始變小。但在9、10月份出現(xiàn)了CH4通量峰值,Mastepanov等[19]研究也發(fā)現(xiàn)秋季的排放峰值,并稱之為“a late-autumn shoulder”。在關(guān)于張掖濕地水位的介紹中提到,10月份由于農(nóng)田不再需要灌溉,濕地水位會(huì)升高。這有可能使得之前裸露地表被水覆蓋,造成更大面積的厭氧環(huán)境,加之夏季植被生物量增加,土壤有機(jī)質(zhì)豐富,為CH4的產(chǎn)生提供了充足的底物,因此有更多的CH4氣體產(chǎn)生。關(guān)于這一解釋有待進(jìn)一步探討。

    有關(guān)CH4通量季節(jié)變化的研究眾多,不同地區(qū)變化規(guī)律有一定差別。Singh等[20]發(fā)現(xiàn)印度北部濕地CH4排放最大值出現(xiàn)在夏季,其次是雨季,冬季最小。Khalila等[21]研究發(fā)現(xiàn)北半球大氣中CH4出現(xiàn)兩個(gè)峰值,分別在春季(4月)和秋季(10月)。Song等[22]在中國(guó)三江平原的淡水濕地發(fā)現(xiàn)CH4排放從5月開(kāi)始增加,到7月達(dá)到最大值((30.5±23.5) mg C m-2h-1),8月中旬開(kāi)始下降。張掖濕地CH4排放高值出現(xiàn)在夏季(植被生長(zhǎng)旺盛)、秋季(水位升高)以及春季(冰雪融化)。不同濕地類型CH4排放的差別與影響因子差異有關(guān)。表2可以看出不同濕地類型CH4年排放量差別較大,在3—33 g CH4m-2a-1之間,但均小于張掖濕地的CH4排放量。

    表2 不同濕地類型的CH4年排放量比較

    3.2CH4通量影響因子分析

    3.2.1 溫度

    溫度作為影響濕地CH4排放的主要?dú)夂蛞蜃?對(duì)CH4的產(chǎn)生、氧化以及傳輸過(guò)程都會(huì)產(chǎn)生影響[27],CH4通量的時(shí)間動(dòng)態(tài)(日變化和季節(jié)變化)與之也有一定的相關(guān)性[1]。產(chǎn)甲烷菌活性的最佳溫度在30—40℃,低溫會(huì)導(dǎo)致產(chǎn)甲烷菌及參與甲烷發(fā)酵過(guò)程的其他微生物活性降低從而減少CH4通量[1]。

    大量研究認(rèn)為CH4通量與空氣溫度或土壤溫度之間有明顯的指數(shù)相關(guān)關(guān)系。Herbst等[28]發(fā)現(xiàn)日CH4通量與20 cm處土壤溫度存在指數(shù)關(guān)系。Hanis等[26]認(rèn)為CH4通量與空氣溫度和50 cm處土壤溫度之間呈指數(shù)相關(guān)關(guān)系。Song等[22]則發(fā)現(xiàn)CH4通量與5 cm處土壤溫度之間的相關(guān)關(guān)系最為明顯。

    3.2.2 摩擦風(fēng)速

    近地表的湍流與垂直風(fēng)速密切相關(guān),諸多研究都表明CH4通量與大氣湍流之間有一定的相關(guān)關(guān)系[29-31]。CH4通量依賴近地表的湍流可解釋為水面和大氣之間氣體傳輸?shù)臄U(kuò)散和湍流。湍流和風(fēng)速的增大可能會(huì)產(chǎn)生一個(gè)薄的層流邊界層,導(dǎo)致在富集CH4的土壤層與湍流邊界層之間有更高的CH4濃度梯度,從而加強(qiáng)了CH4通量的擴(kuò)散[32]。此外,湍流的增強(qiáng)可以增加在平靜期儲(chǔ)存的CH4的通氣和交換[33]。

    3.3 碳源匯分析

    受氣候變化等因素的影響,濕地生態(tài)系統(tǒng)碳吸收與碳排放之間的平衡常常被打破。濕地能否繼續(xù)作為碳匯,或即將、已經(jīng)成為碳源問(wèn)題受到廣泛關(guān)注[34]。B?ckstrand等[35]發(fā)現(xiàn)全球變暖使得北方高緯地區(qū)富含有機(jī)碳的永久凍土開(kāi)始融化、分解,由之前的碳匯變?yōu)樘荚?。Koehler等[36]對(duì)愛(ài)爾蘭北方泥炭地進(jìn)行了連續(xù)6a的觀測(cè),發(fā)現(xiàn)其中有兩年CH4和溶解性有機(jī)碳(DOC)的損失之和超過(guò)了CO2吸收量。Waddington等[37]發(fā)現(xiàn)在夏季較暖或干旱的年份,北方泥炭地為大氣碳源。

    濕地生態(tài)系統(tǒng)碳源匯的年際差異,主要來(lái)自于氣候條件的年際變化。另外,CO2、CH4對(duì)溫度的響應(yīng)及濕地面積的減少,可能導(dǎo)致厭氧與好氧環(huán)境相互轉(zhuǎn)換,這都可能造成張掖濕地碳源匯功能發(fā)生轉(zhuǎn)變[38]。

    4 結(jié)論

    本研究通過(guò)對(duì)張掖濕地CH4通量分析,得出以下結(jié)論:

    (1)CH4通量的日變化趨勢(shì)總體表現(xiàn)為白天大于夜間,在春季和冬季其變化趨勢(shì)均呈倒“U”型分布,且峰值集中在正午前后;不同季節(jié)CH4通量排放特征差異明顯。夏季最大,春秋次之,冬季最?。?/p>

    (2)CH4通量日總量與空氣溫度之間指數(shù)相關(guān)關(guān)系顯著。不同深度土壤溫度中,4 cm處土壤溫度與之相關(guān)性最強(qiáng),因此4 cm處土壤溫度更能代表產(chǎn)甲烷菌活性的溫度條件。1—6月U*與CH4通量顯著正相關(guān);

    (3)2012年7月至2014年6月年張掖濕地平均年CH4排放量為37.96 g CH4m-2a-1。結(jié)合CO2通量觀測(cè)數(shù)據(jù),研究時(shí)段張掖濕地凈吸收碳495.92 g C m-2a-1為碳匯。

    致謝:本文數(shù)據(jù)來(lái)自HiWATER試驗(yàn),感謝HiWATER試驗(yàn)全體參加人員。

    [1] 陳槐, 周舜, 吳寧, 王艷芬, 羅鵬, 石福孫. 濕地甲烷的產(chǎn)生、氧化及排放通量研究進(jìn)展. 應(yīng)用與環(huán)境生物學(xué)報(bào), 2006, 12(5): 726- 733.

    [2] 張堅(jiān)超, 徐鐿欽, 陸雅海. 陸地生態(tài)系統(tǒng)甲烷產(chǎn)生和氧化過(guò)程的微生物機(jī)理. 生態(tài)學(xué)報(bào), 2015, 35(20): 6592- 6603.

    [3] 沙晨燕, 譚娟, 王卿, 王敏. 不同類型河濱濕地甲烷和二氧化碳排放初步研究. 生態(tài)環(huán)境學(xué)報(bào), 2015, 24(7): 1182- 1190.

    [4] 欒軍偉, 崔麗娟, 宋洪濤, 王義飛. 國(guó)外濕地生態(tài)系統(tǒng)碳循環(huán)研究進(jìn)展. 濕地科學(xué), 2012, 10(2): 235- 242.

    [5] 黃國(guó)宏, 肖篤寧, 李玉祥, 陳冠雄, 楊玉成, 趙長(zhǎng)偉. 蘆葦濕地溫室氣體甲烷(CH4)排放研究. 生態(tài)學(xué)報(bào), 2001, 21(9): 1494- 1497.

    [6] Kang H, Freeman C. The influence of hydrochemistry on methane emissions from two contrasting northern wetlands. Water, Air, and Soil Pollution, 2002, 141(1/4): 263- 272.

    [7] Rinne J, Riutta T, Pihlatie M, Aurela M, Haapanala S, Tuovinen J P, Tuittila E S, Vesala T. Annual cycle of methane emission from a boreal fen measured by the eddy covariance technique. Tellus Series B-Chemical and Physical Meteorology, 2007, 59(3): 449- 457.

    [8] Long K D, Flanagan L B, Cai T B. Diurnal and seasonal variation in methane emissions in a northern Canadian peatland measured by eddy covariance. Global Change Biology, 2010, 16(9): 2420- 2435.

    [9] 牛贇, 劉賢德, 張宏斌, 孟好軍. 黑河流域中上游濕地生態(tài)功能評(píng)價(jià). 濕地科學(xué), 2007, 5(3): 215- 220.

    [10] Li X, Cheng G D, Liu S M, Xiao Q, Ma M G, Jin R, Che T, Liu Q H, Wang W Z, Qi Y, Wen J G, Li H Y, Zhu G F, Guo J W, Ran Y H, Wang S G, Zhu Z L, Zhou J, Hu X L, Xu Z W. Heihe watershed allied telemetry experimental research (HiWATER): scientific objectives and experimental design. Bulletin of the American Meteorological Society, 2013, 94(8): 1145- 1160.

    [11] 李新, 劉紹民, 馬明國(guó), 肖青, 柳欽火, 晉銳, 車濤, 王維真, 祁元, 李弘毅, 朱高峰, 郭建文, 冉有華, 聞建光, 王樹(shù)罘. 黑河流域生態(tài)——水文過(guò)程綜合遙感觀測(cè)聯(lián)合試驗(yàn)總體設(shè)計(jì). 地球科學(xué)進(jìn)展, 2012, 27(5): 481- 498.

    [12] Liu S M, Xu Z W, Wang W Z, Jia Z Z, Zhu M J, Bai J, Wang J M. A comparison of eddy-covariance and large aperture scintillometer measurements with respect to the energy balance closure problem. Hydrologyand Earth System Sciences, 2011, 15: 1291- 1306.

    [13] 于貴瑞, 孫曉敏. 陸地生態(tài)系統(tǒng)通量觀測(cè)的原理與方法. 北京: 高等教育出版社, 2006: 212- 223.

    [14] 李霞, 孫睿, 李遠(yuǎn), 王修信, 謝東輝, 嚴(yán)曉丹, 朱啟疆. 北京海淀公園綠地二氧化碳通量. 生態(tài)學(xué)報(bào), 2010, 30(24): 6715- 6725.

    [15] Wille C, Kutzbach L, Sachs T, Wagner D, Pfeiffer E. Methane emission from Siberian arctic polygonal tundra: eddy covariance measurements and modeling. Global Change Biology, 2008, 14(6): 1395- 1408.

    [16] Zhang Q, Sun R, Jiang G Q, Xu Z W, Liu S M. Carbon and energy flux from aPhragmitesaustraliswetland in Zhangye oasis-desert area, China. Agricultural and Forest Meteorology, 2016, doi:10.1016/j.agrformet.2016.02.019.

    [17] 郝慶菊, 王躍思, 江長(zhǎng)勝, 王長(zhǎng)科, 王明星. 濕地甲烷排放研究若干問(wèn)題的探討. 生態(tài)學(xué)雜志, 2005, 24(2): 170- 175.

    [18] Juutinen S, Alm J, Larmola T, Huttunen J T, Morero M, Saarnio S, Martikainen P J, Silvola J. Methane (CH4) release from littoral wetlands of boreal lakes during an extended flooding period. Global Change Biology, 2003, 9(3): 413- 424.

    [19] Mastepanov M, Sigsgaard C, Dlugokencky E J, Houweling S, Str?m L, Tamstorf M P, Christensen T R. Large tundra methane burst during onset of freezing. Nature, 2008, 456(7222): 628- 630.

    [20] Singh S N, Kulshreshtha K, Agnihotri S. Seasonal dynamics of methane emission from wetlands. Chemosphere-Global Change Science, 2000, 2(1): 39- 46.

    [21] Khalil M A K, Rasmussen R A. Sources, sinks, and seasonal cycles of atmospheric methane. Journal of Geophysical Research: Oceans, 1983, 88(C9): 5131- 5144.

    [22] Song C C, Sun L, Huang Y, Wang Y S, Wan Z M. Carbon exchange in a freshwater marsh in the Sanjiang Plain, northeastern China. Agricultural and Forest Meteorology, 2011, 151(8): 1131- 1138.

    [23] Nyk?nen H, Heikkinen J E P, Pirinen L, Tiilikainen K, Martikainen P J. Annual CO2exchange and CH4fluxes on a subarctic palsa mire during climatically different years. Global Biogeochemical Cycles, 2003, 17(1): 1018.

    [24] Zhou L, Zhou G S, Jia Q Y. Annual cycle of CO2exchange over a reed (Phragmitesaustralis) wetland in Northeast China. Aquatic Botany, 2009, 91(2): 91- 98.

    [25] Jackowicz-Korczyński M, Christensen T R, B?ckstrand K, Crill P, Friborg T, Mastepanov M, Str?m L. Annual cycle of methane emission from a subarctic peatland. Journal of Geophysical Research: Biogeosciences, 2010, 115(G2), doi: 10.1029/2008JG000913.

    [26] Hanis K L, Tenuta M, Amiro B D, Papakyriakou T N. Seasonal dynamics of methane emissions from a subarctic fen in the Hudson Bay Lowlands. Biogeosciences, 2013, 10(7): 4465- 4479.

    [27] Cao M K, Gregson K, Marshall S. Global methane emission from wetlands and its sensitivity to climate change. Atmospheric Environment, 1998, 32(19): 3293- 3299.

    [28] Herbst M, Friborg T, Ringgaard R, Soegaard H. Interpreting the variations in atmospheric methane fluxes observed above a restored wetland. Agricultural and Forest Meteorology, 2011, 151(7): 841- 853.

    [29] Beringer J, Livesley S J, Randle J, Hutley L B. Carbon dioxide fluxes dominate the greenhouse gas exchanges of a seasonal wetland in the wet-dry tropics of northern Australia. Agricultural and Forest Meteorology, 2013, 182- 183: 239- 247.

    [30] 仝川, 鄂焱, 廖稷, 姚順, 王維奇, 黃佳芳, 張林海, 楊紅玉, 曾從盛. 閩江河口潮汐沼澤濕地CO2排放通量特征. 環(huán)境科學(xué)學(xué)報(bào), 2011, 31(12): 2830- 2840.

    [31] Fan S M, Wofsy S C, Bakwin P S, Jacob D J, Anderson S M, Kebabian P L, Mcmanus J B, Kolb C E, Fitzjarrald D R. Micrometeorological measurements of CH4and CO2exchange between the atmosphere and subarctic tundra. Journal of Geophysical Research: Atmospheres, 1992, 97(D15): 16627- 16643.

    [32] Sachs T, Wille C, Boike J, Kutzbach L. Environmental controls on ecosystem-scale CH4emission from polygonal tundra in the Lena River Delta, Siberia. Journal of Geophysical Research: Biogeosciences, 2008, 113(G3), doi: 10.1029/2007JG000505

    [33] Hargreaves K J, Fowler D, Pitcairn C E R, Aurela M. Annual methane emission from Finnish mires estimated from eddy covariance campaign measurements. Theoretical and Applied Climatology, 2001, 70(1/4): 203- 213.

    [34] 楊永興. 國(guó)際濕地科學(xué)研究的主要特點(diǎn)、進(jìn)展與展望. 地理科學(xué)進(jìn)展, 2002, 21(2): 111- 120.

    [35] B?ckstrand K, Crill P M, Jackowicz-Korczyski M, Mastepanov M, Christensen T R, Bastviken D. Annual carbon gas budget for a subarctic peatland, Northern Sweden. Biogeosciences, 2010, 7(1): 95- 108.

    [36] Koehler A K, Sottocornola M, Kiely G. How strong is the current carbon sequestration of an Atlantic blanket bog? Global Change Biology, 2011, 17(1): 309- 319.

    [37] Waddington J M, Roulet N T. Carbon balance of a boreal patterned peatland. Global Change Biology, 2000, 6(1): 87- 97.

    [38] IPCC. Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge, United Kingdom and New York, USA, Cambridge University Press, 2013.

    CH4fluxvariationsandmainfactorsfromareedwetlandoasis-desertareainZhangye,China

    ZHANG Qiang1,2,3, JIANG Guoqing1,2,3, SUN Rui1,2,3,*, XU Ziwei1,2, LIU Shaomin1,2

    1StateKeyLaboratoryofRemoteSensingScience,JointlySponsoredbyBeijingNormalUniversityandtheInstituteofRemoteSensingApplications,ChineseAcademyofSciences,Beijing100875,China2SchoolofGeographyandRemoteSensingSciences,BeijingNormalUniversity,Beijing100875,China3BeijingKeyLaboratoryforRemoteSensingofEnvironmentandDigitalCities,Beijing100875,China

    國(guó)家自然科學(xué)基金(41471349);國(guó)家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目(2016YFB0501502);國(guó)家科技支撐計(jì)劃項(xiàng)目(2013BAC03B02);中央高校基本科研業(yè)務(wù)費(fèi)專項(xiàng)(2014kJJCA02)

    2016- 06- 13; < class="emphasis_bold">網(wǎng)絡(luò)出版日期

    日期:2017- 04- 24

    10.5846/stxb201606131130

    *通訊作者Corresponding author.E-mail: sunrui@bnu.edu.cn

    張強(qiáng),蔣國(guó)慶,孫睿,徐自為,劉紹民.張掖濕地甲烷通量動(dòng)態(tài)特征及其影響因子.生態(tài)學(xué)報(bào),2017,37(17):5681- 5690.

    Zhang Q, Jiang G Q, Sun R, Xu Z W, Liu S M.CH4flux variations and main factors from a reed wetland oasis-desert area in Zhangye, China.Acta Ecologica Sinica,2017,37(17):5681- 5690.

    猜你喜歡
    張掖土壤溫度通量
    冬小麥田N2O通量研究
    情暖張掖大地 讓愛(ài)不再孤單
    到張掖看黑河
    文學(xué)港(2019年5期)2019-05-24 14:19:42
    遼東山區(qū)3種人工林土壤呼吸對(duì)土壤溫度和土壤水分的響應(yīng)
    西夏時(shí)期的張掖
    西夏學(xué)(2016年2期)2016-10-26 02:20:26
    大美張掖
    絲綢之路(2016年13期)2016-08-10 07:58:08
    管群間歇散熱的土壤溫度響應(yīng)與恢復(fù)特性
    土壤濕度和土壤溫度模擬中的參數(shù)敏感性分析和優(yōu)化
    緩釋型固體二氧化氯的制備及其釋放通量的影響因素
    春、夏季長(zhǎng)江口及鄰近海域溶解甲烷的分布與釋放通量
    色哟哟·www| 中文字幕人妻熟女乱码| 叶爱在线成人免费视频播放| 亚洲男人天堂网一区| 国产免费又黄又爽又色| 欧美成人午夜免费资源| 18禁动态无遮挡网站| 99久久精品国产国产毛片| 天美传媒精品一区二区| 欧美日韩亚洲国产一区二区在线观看 | 99久久综合免费| 制服人妻中文乱码| 性少妇av在线| 日韩制服骚丝袜av| 日本色播在线视频| 男女边摸边吃奶| 肉色欧美久久久久久久蜜桃| 久久精品人人爽人人爽视色| 午夜福利一区二区在线看| 一级爰片在线观看| 熟妇人妻不卡中文字幕| 韩国高清视频一区二区三区| 精品久久久精品久久久| 免费观看性生交大片5| 亚洲精品国产色婷婷电影| xxx大片免费视频| 亚洲综合色网址| 搡老乐熟女国产| 久久久久久伊人网av| 一级片免费观看大全| 亚洲欧美精品自产自拍| 中文天堂在线官网| 天堂中文最新版在线下载| 久久国产精品大桥未久av| 女性生殖器流出的白浆| 伦理电影大哥的女人| 女人久久www免费人成看片| 久久av网站| 久久毛片免费看一区二区三区| 久久精品熟女亚洲av麻豆精品| 国产xxxxx性猛交| 亚洲av在线观看美女高潮| 精品亚洲乱码少妇综合久久| av线在线观看网站| av在线app专区| 精品亚洲乱码少妇综合久久| 国产欧美日韩一区二区三区在线| 最近的中文字幕免费完整| 9热在线视频观看99| 午夜日韩欧美国产| 亚洲综合精品二区| 精品久久久精品久久久| videossex国产| 好男人视频免费观看在线| 99久久综合免费| 一级毛片我不卡| 日日爽夜夜爽网站| 久久韩国三级中文字幕| 亚洲精品国产av成人精品| 亚洲欧美中文字幕日韩二区| 97精品久久久久久久久久精品| 欧美日韩精品网址| 久久久久精品性色| a 毛片基地| 97精品久久久久久久久久精品| 久久久精品区二区三区| 9热在线视频观看99| 黄色视频在线播放观看不卡| 在线免费观看不下载黄p国产| 亚洲国产精品成人久久小说| 欧美成人午夜精品| 精品少妇一区二区三区视频日本电影 | 中文字幕亚洲精品专区| 纯流量卡能插随身wifi吗| 涩涩av久久男人的天堂| 精品国产一区二区三区四区第35| 1024香蕉在线观看| 中文字幕人妻丝袜一区二区 | 国产免费现黄频在线看| 色网站视频免费| 日韩人妻精品一区2区三区| 美女主播在线视频| 亚洲av成人精品一二三区| 日韩不卡一区二区三区视频在线| 国产一区亚洲一区在线观看| 色哟哟·www| 91在线精品国自产拍蜜月| 国产在视频线精品| 亚洲国产精品成人久久小说| 建设人人有责人人尽责人人享有的| 热re99久久精品国产66热6| 国产精品久久久久久久久免| 伊人久久大香线蕉亚洲五| 国产亚洲精品第一综合不卡| 亚洲三级黄色毛片| 久久久欧美国产精品| 九草在线视频观看| 久久人妻熟女aⅴ| 日韩av在线免费看完整版不卡| 黄片无遮挡物在线观看| 亚洲美女视频黄频| 免费不卡的大黄色大毛片视频在线观看| 欧美另类一区| 久久午夜综合久久蜜桃| 久久韩国三级中文字幕| 热99国产精品久久久久久7| 精品国产超薄肉色丝袜足j| 亚洲美女搞黄在线观看| 一个人免费看片子| 人妻 亚洲 视频| 欧美人与善性xxx| 欧美日韩综合久久久久久| 天天操日日干夜夜撸| 日本色播在线视频| 一本色道久久久久久精品综合| 亚洲精品国产色婷婷电影| 大片免费播放器 马上看| 亚洲国产精品成人久久小说| 婷婷成人精品国产| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 久久久久精品人妻al黑| 人人妻人人爽人人添夜夜欢视频| 亚洲成色77777| 满18在线观看网站| 波多野结衣一区麻豆| 日韩大片免费观看网站| 午夜久久久在线观看| 久久久久久久久免费视频了| 亚洲国产最新在线播放| 在线免费观看不下载黄p国产| 免费看av在线观看网站| 熟妇人妻不卡中文字幕| kizo精华| 国产男女超爽视频在线观看| 亚洲精品日韩在线中文字幕| 久久这里只有精品19| 亚洲精品国产一区二区精华液| 亚洲精品久久午夜乱码| 咕卡用的链子| 国产乱来视频区| 亚洲av在线观看美女高潮| 午夜福利乱码中文字幕| a级片在线免费高清观看视频| 国产男女超爽视频在线观看| 一级片免费观看大全| 国产麻豆69| 久久人人爽av亚洲精品天堂| 性色avwww在线观看| 国产精品国产三级国产专区5o| 自线自在国产av| 精品国产一区二区三区四区第35| 90打野战视频偷拍视频| 2022亚洲国产成人精品| 亚洲国产av影院在线观看| 有码 亚洲区| 在线观看免费高清a一片| 丝袜脚勾引网站| 超碰97精品在线观看| 精品亚洲成a人片在线观看| 久久毛片免费看一区二区三区| 欧美日本中文国产一区发布| 免费观看av网站的网址| 亚洲精品自拍成人| 国产精品一国产av| 日日撸夜夜添| 中文字幕av电影在线播放| 9热在线视频观看99| 少妇精品久久久久久久| 欧美黄色片欧美黄色片| 久久精品国产亚洲av高清一级| 国产一区二区三区av在线| av又黄又爽大尺度在线免费看| 日韩成人av中文字幕在线观看| www日本在线高清视频| 久久久久久免费高清国产稀缺| 女人高潮潮喷娇喘18禁视频| 赤兔流量卡办理| 在线天堂最新版资源| 欧美日韩成人在线一区二区| 亚洲国产av影院在线观看| 青春草亚洲视频在线观看| 一本大道久久a久久精品| 亚洲国产成人一精品久久久| 久热这里只有精品99| 99精国产麻豆久久婷婷| 国产精品无大码| 欧美少妇被猛烈插入视频| 欧美国产精品一级二级三级| 国产欧美日韩综合在线一区二区| 老女人水多毛片| 一区二区三区激情视频| 国产免费又黄又爽又色| 乱人伦中国视频| 中文字幕精品免费在线观看视频| 欧美精品一区二区免费开放| 国产乱人偷精品视频| 久久国产精品男人的天堂亚洲| 欧美激情 高清一区二区三区| 在线看a的网站| 一区二区av电影网| 丰满迷人的少妇在线观看| 黄色视频在线播放观看不卡| 日本91视频免费播放| 日本黄色日本黄色录像| 亚洲av欧美aⅴ国产| 亚洲综合色惰| 日日摸夜夜添夜夜爱| 国产精品99久久99久久久不卡 | 欧美xxⅹ黑人| 久久国产精品大桥未久av| 天天躁狠狠躁夜夜躁狠狠躁| www.自偷自拍.com| 伊人久久大香线蕉亚洲五| 亚洲伊人色综图| 乱人伦中国视频| 丁香六月天网| 波多野结衣av一区二区av| 人妻 亚洲 视频| 人人妻人人爽人人添夜夜欢视频| 美女高潮到喷水免费观看| 久久久久久久久久久免费av| 国产av码专区亚洲av| 国产精品二区激情视频| 99国产精品免费福利视频| 亚洲中文av在线| 免费高清在线观看日韩| 天堂俺去俺来也www色官网| 下体分泌物呈黄色| 亚洲国产最新在线播放| av在线播放精品| 9191精品国产免费久久| 一本久久精品| 亚洲国产精品成人久久小说| 日韩中文字幕欧美一区二区 | 亚洲精品一区蜜桃| www.自偷自拍.com| 女人高潮潮喷娇喘18禁视频| 巨乳人妻的诱惑在线观看| 精品99又大又爽又粗少妇毛片| 天天躁夜夜躁狠狠躁躁| 巨乳人妻的诱惑在线观看| 少妇被粗大的猛进出69影院| 黄片小视频在线播放| 黄色一级大片看看| 亚洲欧美一区二区三区黑人 | 99香蕉大伊视频| 日韩不卡一区二区三区视频在线| 日韩一卡2卡3卡4卡2021年| 一级毛片 在线播放| 国产亚洲欧美精品永久| 美女xxoo啪啪120秒动态图| 美女福利国产在线| av一本久久久久| 久久影院123| 中文字幕av电影在线播放| 老司机影院成人| 9热在线视频观看99| 亚洲欧美日韩另类电影网站| 丰满饥渴人妻一区二区三| 97人妻天天添夜夜摸| 最近手机中文字幕大全| 久久久久久久精品精品| 亚洲av.av天堂| 欧美日韩成人在线一区二区| 色网站视频免费| 久久婷婷青草| 久久精品国产a三级三级三级| 欧美国产精品va在线观看不卡| 亚洲,一卡二卡三卡| 日本猛色少妇xxxxx猛交久久| 亚洲国产av新网站| 狠狠婷婷综合久久久久久88av| 18禁国产床啪视频网站| av在线老鸭窝| 亚洲国产成人一精品久久久| 丝袜脚勾引网站| 久久韩国三级中文字幕| 欧美少妇被猛烈插入视频| 久久99蜜桃精品久久| 国产精品av久久久久免费| 国产精品不卡视频一区二区| 赤兔流量卡办理| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 精品人妻在线不人妻| 国产毛片在线视频| 在线天堂中文资源库| 性高湖久久久久久久久免费观看| 国产男女超爽视频在线观看| 成人毛片60女人毛片免费| 男人操女人黄网站| 亚洲av福利一区| 丝袜美足系列| 亚洲欧美日韩另类电影网站| 自线自在国产av| 少妇被粗大的猛进出69影院| 午夜福利网站1000一区二区三区| freevideosex欧美| 亚洲色图 男人天堂 中文字幕| 成人国产av品久久久| 色视频在线一区二区三区| 一本—道久久a久久精品蜜桃钙片| 亚洲精品国产色婷婷电影| xxxhd国产人妻xxx| 不卡av一区二区三区| 韩国高清视频一区二区三区| 老熟女久久久| 少妇的逼水好多| 在线天堂中文资源库| 熟女av电影| 久久女婷五月综合色啪小说| a级片在线免费高清观看视频| 欧美最新免费一区二区三区| 天天躁夜夜躁狠狠久久av| 最黄视频免费看| 日韩av免费高清视频| 成人18禁高潮啪啪吃奶动态图| videossex国产| 99热全是精品| 狠狠精品人妻久久久久久综合| 亚洲少妇的诱惑av| 成人国产麻豆网| 精品国产一区二区久久| 日韩av免费高清视频| 青青草视频在线视频观看| 中文字幕精品免费在线观看视频| 亚洲欧美精品综合一区二区三区 | 王馨瑶露胸无遮挡在线观看| √禁漫天堂资源中文www| 97在线视频观看| 国产亚洲最大av| 国产免费一区二区三区四区乱码| a级毛片黄视频| 男人舔女人的私密视频| 美女中出高潮动态图| 不卡视频在线观看欧美| 成人亚洲精品一区在线观看| 欧美变态另类bdsm刘玥| 99国产精品免费福利视频| 久久人妻熟女aⅴ| 国产精品一国产av| 欧美中文综合在线视频| 一本久久精品| 国产精品久久久久久精品古装| 纯流量卡能插随身wifi吗| 久久久久久久国产电影| 国产乱人偷精品视频| 18禁国产床啪视频网站| 久久亚洲国产成人精品v| 女人被躁到高潮嗷嗷叫费观| 亚洲三区欧美一区| 国产成人aa在线观看| 欧美 亚洲 国产 日韩一| 91精品伊人久久大香线蕉| 9191精品国产免费久久| 亚洲婷婷狠狠爱综合网| 久久青草综合色| 日本爱情动作片www.在线观看| 成人国产麻豆网| 观看av在线不卡| 你懂的网址亚洲精品在线观看| 色网站视频免费| 亚洲国产精品一区二区三区在线| 国产福利在线免费观看视频| 国产黄色免费在线视频| 亚洲精品乱久久久久久| 99九九在线精品视频| 18禁观看日本| 宅男免费午夜| 天天躁夜夜躁狠狠久久av| 午夜激情久久久久久久| 咕卡用的链子| 国产亚洲精品第一综合不卡| 精品国产乱码久久久久久男人| 日本wwww免费看| 午夜免费鲁丝| 性色av一级| 亚洲在久久综合| 青草久久国产| av一本久久久久| 美女国产高潮福利片在线看| 在线观看人妻少妇| 性高湖久久久久久久久免费观看| 国产精品无大码| 人人妻人人添人人爽欧美一区卜| 少妇 在线观看| 免费少妇av软件| 国产精品.久久久| 久久精品亚洲av国产电影网| 国产男人的电影天堂91| 一本久久精品| 国产日韩欧美亚洲二区| 国产精品一区二区在线观看99| 国产免费现黄频在线看| 亚洲精品久久久久久婷婷小说| 日本猛色少妇xxxxx猛交久久| 国产精品99久久99久久久不卡 | 天美传媒精品一区二区| 国产精品久久久久成人av| 国产日韩欧美亚洲二区| 丝袜美腿诱惑在线| 亚洲欧美成人综合另类久久久| 亚洲国产欧美网| 丝袜美足系列| 亚洲精品中文字幕在线视频| 秋霞伦理黄片| 一区二区av电影网| 热re99久久国产66热| 日韩制服丝袜自拍偷拍| 九草在线视频观看| 精品人妻偷拍中文字幕| 亚洲成国产人片在线观看| 热re99久久国产66热| 亚洲美女搞黄在线观看| 美女高潮到喷水免费观看| 日本猛色少妇xxxxx猛交久久| 97人妻天天添夜夜摸| 宅男免费午夜| 大码成人一级视频| 下体分泌物呈黄色| 男女无遮挡免费网站观看| 欧美最新免费一区二区三区| 男女边吃奶边做爰视频| 免费黄色在线免费观看| 日韩在线高清观看一区二区三区| 日韩av不卡免费在线播放| 国产精品人妻久久久影院| 我要看黄色一级片免费的| 高清在线视频一区二区三区| 伦理电影大哥的女人| 搡女人真爽免费视频火全软件| 丰满少妇做爰视频| 亚洲欧美色中文字幕在线| 寂寞人妻少妇视频99o| 各种免费的搞黄视频| 在现免费观看毛片| 久久久a久久爽久久v久久| av一本久久久久| 丝袜脚勾引网站| 男女无遮挡免费网站观看| 精品人妻偷拍中文字幕| 中国国产av一级| 18禁动态无遮挡网站| www.av在线官网国产| 又大又黄又爽视频免费| 好男人视频免费观看在线| 在线观看人妻少妇| 亚洲精品自拍成人| 少妇猛男粗大的猛烈进出视频| 欧美成人午夜免费资源| 黄色一级大片看看| 亚洲人成电影观看| 伦理电影免费视频| 人妻系列 视频| 久久国产精品男人的天堂亚洲| 妹子高潮喷水视频| 亚洲第一青青草原| 国产白丝娇喘喷水9色精品| 一区二区av电影网| 我的亚洲天堂| 久久久久久久大尺度免费视频| 亚洲久久久国产精品| 人成视频在线观看免费观看| 亚洲一级一片aⅴ在线观看| 美国免费a级毛片| 亚洲国产看品久久| 成人毛片60女人毛片免费| 日韩av免费高清视频| 亚洲五月色婷婷综合| 宅男免费午夜| 日韩av不卡免费在线播放| 高清av免费在线| 一二三四在线观看免费中文在| 日韩欧美一区视频在线观看| 男女边摸边吃奶| 国产成人免费无遮挡视频| 女性被躁到高潮视频| 天堂俺去俺来也www色官网| 国产探花极品一区二区| 99热国产这里只有精品6| 欧美日韩国产mv在线观看视频| 亚洲伊人色综图| 亚洲第一区二区三区不卡| 一区在线观看完整版| 色吧在线观看| 国产片特级美女逼逼视频| 尾随美女入室| 欧美人与性动交α欧美精品济南到 | 99国产精品免费福利视频| 国产欧美日韩综合在线一区二区| 精品少妇内射三级| 成人二区视频| 黄片小视频在线播放| 欧美日韩视频精品一区| 黄片播放在线免费| 亚洲国产欧美日韩在线播放| 国产成人午夜福利电影在线观看| 欧美亚洲 丝袜 人妻 在线| 亚洲,一卡二卡三卡| 亚洲av在线观看美女高潮| 黄片小视频在线播放| 中文欧美无线码| 免费观看无遮挡的男女| av在线app专区| 成人毛片a级毛片在线播放| 9色porny在线观看| 三上悠亚av全集在线观看| 最近中文字幕2019免费版| 日日撸夜夜添| www日本在线高清视频| 免费看不卡的av| 少妇被粗大猛烈的视频| 国产97色在线日韩免费| 夜夜骑夜夜射夜夜干| 日韩精品免费视频一区二区三区| 黑丝袜美女国产一区| 午夜激情久久久久久久| 丰满饥渴人妻一区二区三| 高清视频免费观看一区二区| 久久精品国产自在天天线| 高清不卡的av网站| 久久精品国产a三级三级三级| av.在线天堂| 在线观看国产h片| 18禁观看日本| 欧美人与善性xxx| 极品人妻少妇av视频| 欧美日韩精品网址| 亚洲国产日韩一区二区| 桃花免费在线播放| 亚洲欧美一区二区三区国产| 久久久久国产网址| 日韩av不卡免费在线播放| 黑人欧美特级aaaaaa片| 女人被躁到高潮嗷嗷叫费观| 90打野战视频偷拍视频| 午夜激情久久久久久久| 免费黄色在线免费观看| 国产成人一区二区在线| 最黄视频免费看| 搡老乐熟女国产| 久久精品国产a三级三级三级| 丝袜美腿诱惑在线| 黑人猛操日本美女一级片| 飞空精品影院首页| 婷婷色综合大香蕉| 老司机亚洲免费影院| 欧美xxⅹ黑人| 一区二区三区四区激情视频| 亚洲第一av免费看| 国产成人免费观看mmmm| 啦啦啦啦在线视频资源| 日本免费在线观看一区| 免费在线观看黄色视频的| 欧美成人午夜免费资源| 最近最新中文字幕大全免费视频 | 热99国产精品久久久久久7| 午夜福利在线观看免费完整高清在| 色视频在线一区二区三区| 免费黄频网站在线观看国产| 观看av在线不卡| 777久久人妻少妇嫩草av网站| 国产不卡av网站在线观看| 国产男女内射视频| 女人被躁到高潮嗷嗷叫费观| 大话2 男鬼变身卡| 欧美日韩av久久| 亚洲成色77777| 国产精品 欧美亚洲| 国产欧美日韩一区二区三区在线| 精品少妇黑人巨大在线播放| 成人国产av品久久久| 考比视频在线观看| 色网站视频免费| 日韩三级伦理在线观看| 精品一品国产午夜福利视频| 色婷婷av一区二区三区视频| 激情视频va一区二区三区| 久久精品国产亚洲av涩爱| 午夜av观看不卡| 中国三级夫妇交换| 夜夜骑夜夜射夜夜干| 国产日韩欧美亚洲二区| 国产黄频视频在线观看| 麻豆精品久久久久久蜜桃| 免费女性裸体啪啪无遮挡网站| 久久综合国产亚洲精品| av线在线观看网站| 亚洲精品日韩在线中文字幕| 亚洲av综合色区一区| 欧美精品一区二区大全| 国产1区2区3区精品| 欧美变态另类bdsm刘玥| 亚洲国产精品一区三区| 亚洲国产看品久久| 人人妻人人澡人人看| 少妇人妻精品综合一区二区| 最新的欧美精品一区二区| av电影中文网址| 欧美激情极品国产一区二区三区| 国产激情久久老熟女| 一区二区日韩欧美中文字幕| 中国三级夫妇交换| 三上悠亚av全集在线观看| 这个男人来自地球电影免费观看 | 久久狼人影院| 日本午夜av视频| 国产黄色视频一区二区在线观看| 美女脱内裤让男人舔精品视频| 天堂8中文在线网| 国产色婷婷99| 亚洲一区二区三区欧美精品| 午夜免费鲁丝| 80岁老熟妇乱子伦牲交| 五月伊人婷婷丁香| 国产老妇伦熟女老妇高清|