• <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)江口及鄰近海域溶解甲烷的分布與釋放通量
    一级毛片我不卡| 亚洲中文av在线| 国产福利在线免费观看视频| 最近最新中文字幕免费大全7| 午夜免费鲁丝| 如何舔出高潮| 男女床上黄色一级片免费看| 大香蕉久久成人网| 日本av手机在线免费观看| 国产亚洲av高清不卡| 无遮挡黄片免费观看| 成人亚洲精品一区在线观看| 欧美人与性动交α欧美精品济南到| 超色免费av| 国产亚洲欧美精品永久| 免费日韩欧美在线观看| 日韩 欧美 亚洲 中文字幕| 国产av码专区亚洲av| 别揉我奶头~嗯~啊~动态视频 | 亚洲av欧美aⅴ国产| 99久久精品国产亚洲精品| 精品少妇久久久久久888优播| 国产精品一二三区在线看| 80岁老熟妇乱子伦牲交| 国产精品国产三级专区第一集| 精品午夜福利在线看| 婷婷色综合www| 国产亚洲精品第一综合不卡| 在线看a的网站| 少妇的丰满在线观看| 国产免费一区二区三区四区乱码| 国产深夜福利视频在线观看| 色网站视频免费| 欧美精品人与动牲交sv欧美| 夫妻午夜视频| www.精华液| 成人黄色视频免费在线看| 免费日韩欧美在线观看| 国产在视频线精品| 国产在视频线精品| 久久人人爽人人片av| 伊人久久大香线蕉亚洲五| 久久精品熟女亚洲av麻豆精品| 99国产精品免费福利视频| 99国产综合亚洲精品| 亚洲av国产av综合av卡| bbb黄色大片| 国产亚洲午夜精品一区二区久久| 亚洲色图综合在线观看| a 毛片基地| 亚洲专区中文字幕在线 | 久久精品人人爽人人爽视色| 色精品久久人妻99蜜桃| 最近最新中文字幕免费大全7| 在线观看国产h片| 97精品久久久久久久久久精品| 人妻 亚洲 视频| 亚洲av欧美aⅴ国产| www.精华液| 国产精品一区二区精品视频观看| av视频免费观看在线观看| 日韩人妻精品一区2区三区| 少妇的丰满在线观看| 国产精品人妻久久久影院| 女性生殖器流出的白浆| 亚洲图色成人| 丝瓜视频免费看黄片| avwww免费| 亚洲在久久综合| 最近中文字幕2019免费版| 国产精品久久久久久人妻精品电影 | 少妇 在线观看| 韩国高清视频一区二区三区| 精品亚洲成国产av| 波多野结衣av一区二区av| 中文字幕最新亚洲高清| 久久热在线av| 中文字幕av电影在线播放| 无限看片的www在线观看| 男女无遮挡免费网站观看| 亚洲三区欧美一区| 中文精品一卡2卡3卡4更新| 色婷婷av一区二区三区视频| 十分钟在线观看高清视频www| av网站在线播放免费| 一边亲一边摸免费视频| 亚洲综合精品二区| 免费黄色在线免费观看| 久热爱精品视频在线9| 国产无遮挡羞羞视频在线观看| 成人手机av| 天天操日日干夜夜撸| 精品国产一区二区久久| 国产一区二区 视频在线| 一区二区三区乱码不卡18| 在线观看www视频免费| 青春草视频在线免费观看| 久久久国产精品麻豆| 欧美亚洲 丝袜 人妻 在线| 18禁裸乳无遮挡动漫免费视频| 综合色丁香网| 亚洲欧美成人综合另类久久久| 大香蕉久久成人网| 欧美久久黑人一区二区| 黄网站色视频无遮挡免费观看| 宅男免费午夜| 亚洲精品自拍成人| 国产精品欧美亚洲77777| 美女国产高潮福利片在线看| 国产又爽黄色视频| 夫妻性生交免费视频一级片| 又大又黄又爽视频免费| 天美传媒精品一区二区| 性色av一级| 国产一卡二卡三卡精品 | 日韩av在线免费看完整版不卡| 99热国产这里只有精品6| 一本色道久久久久久精品综合| 久久这里只有精品19| e午夜精品久久久久久久| 一本色道久久久久久精品综合| 王馨瑶露胸无遮挡在线观看| 亚洲专区中文字幕在线 | 亚洲欧洲精品一区二区精品久久久 | 久久精品亚洲熟妇少妇任你| 精品亚洲成a人片在线观看| 看十八女毛片水多多多| 在线看a的网站| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品久久久久成人av| 搡老岳熟女国产| 天天操日日干夜夜撸| 欧美中文综合在线视频| 老司机在亚洲福利影院| 青春草视频在线免费观看| 99国产综合亚洲精品| 欧美日韩视频精品一区| 18禁裸乳无遮挡动漫免费视频| 大片免费播放器 马上看| 午夜精品国产一区二区电影| 日韩,欧美,国产一区二区三区| 欧美日韩视频精品一区| 亚洲欧洲精品一区二区精品久久久 | 久久精品国产综合久久久| 亚洲图色成人| 最近中文字幕2019免费版| 亚洲国产av新网站| 黄色怎么调成土黄色| 亚洲欧美激情在线| 亚洲欧美成人精品一区二区| 大香蕉久久成人网| 男女边吃奶边做爰视频| 一本大道久久a久久精品| 欧美日韩亚洲高清精品| 亚洲精华国产精华液的使用体验| 叶爱在线成人免费视频播放| 亚洲天堂av无毛| 亚洲av欧美aⅴ国产| 一本色道久久久久久精品综合| 国产精品熟女久久久久浪| 天天影视国产精品| 国产精品久久久久久人妻精品电影 | 国产一区亚洲一区在线观看| 女人爽到高潮嗷嗷叫在线视频| 自线自在国产av| 久久女婷五月综合色啪小说| 亚洲伊人久久精品综合| 色播在线永久视频| 日本91视频免费播放| 你懂的网址亚洲精品在线观看| 国产一区亚洲一区在线观看| 2021少妇久久久久久久久久久| 精品亚洲成国产av| 亚洲一码二码三码区别大吗| 中文字幕色久视频| 午夜激情av网站| 国产淫语在线视频| 亚洲精品,欧美精品| 日日爽夜夜爽网站| 999久久久国产精品视频| 天天添夜夜摸| 七月丁香在线播放| 美女福利国产在线| 久久久久久免费高清国产稀缺| 另类亚洲欧美激情| 在线观看免费日韩欧美大片| 国产99久久九九免费精品| 一级黄片播放器| 一区二区日韩欧美中文字幕| 最近中文字幕2019免费版| 久久久久视频综合| 久久影院123| 女性被躁到高潮视频| 欧美成人精品欧美一级黄| 久久亚洲国产成人精品v| 欧美亚洲 丝袜 人妻 在线| 久久精品久久精品一区二区三区| 悠悠久久av| 午夜福利,免费看| 视频在线观看一区二区三区| 国产精品女同一区二区软件| 亚洲成人一二三区av| av有码第一页| 精品国产一区二区久久| 啦啦啦中文免费视频观看日本| 啦啦啦啦在线视频资源| 亚洲国产精品999| 久久99精品国语久久久| 国产视频首页在线观看| 精品少妇一区二区三区视频日本电影 | 天天躁夜夜躁狠狠躁躁| 久久久国产精品麻豆| 成人国产麻豆网| 热re99久久国产66热| 亚洲国产精品国产精品| 99精国产麻豆久久婷婷| 国产精品人妻久久久影院| 老司机亚洲免费影院| 青春草视频在线免费观看| 制服丝袜香蕉在线| 久久韩国三级中文字幕| 免费黄频网站在线观看国产| 中文字幕人妻熟女乱码| 日韩一区二区三区影片| 欧美精品一区二区大全| 午夜91福利影院| 宅男免费午夜| 国产成人啪精品午夜网站| 18禁观看日本| 国产黄频视频在线观看| 黄色视频在线播放观看不卡| 美女脱内裤让男人舔精品视频| 久久综合国产亚洲精品| 久久精品熟女亚洲av麻豆精品| 99热国产这里只有精品6| 国产黄色视频一区二区在线观看| 欧美日韩一区二区视频在线观看视频在线| 高清不卡的av网站| 尾随美女入室| 国产亚洲欧美精品永久| 无遮挡黄片免费观看| 亚洲色图 男人天堂 中文字幕| 老司机深夜福利视频在线观看 | 亚洲精品久久成人aⅴ小说| 成人国产麻豆网| 韩国av在线不卡| 中文字幕另类日韩欧美亚洲嫩草| 午夜福利视频精品| 蜜桃在线观看..| 中文字幕色久视频| 高清在线视频一区二区三区| 最黄视频免费看| 国产黄色视频一区二区在线观看| 99国产综合亚洲精品| 青青草视频在线视频观看| 欧美在线黄色| 一本—道久久a久久精品蜜桃钙片| 国产一级毛片在线| 国产精品国产三级专区第一集| 欧美 亚洲 国产 日韩一| av有码第一页| 久久久久久久久免费视频了| 久久婷婷青草| av又黄又爽大尺度在线免费看| 亚洲国产成人一精品久久久| 国产精品一国产av| 在线免费观看不下载黄p国产| 欧美日韩一区二区视频在线观看视频在线| 99久久99久久久精品蜜桃| 97在线人人人人妻| 亚洲av国产av综合av卡| 国产精品99久久99久久久不卡 | 制服人妻中文乱码| h视频一区二区三区| 日韩一卡2卡3卡4卡2021年| 久久精品久久久久久久性| 最近最新中文字幕大全免费视频 | 97在线人人人人妻| 亚洲美女视频黄频| 狠狠精品人妻久久久久久综合| 亚洲人成电影观看| 女性被躁到高潮视频| 夫妻性生交免费视频一级片| 男男h啪啪无遮挡| 宅男免费午夜| 日本猛色少妇xxxxx猛交久久| 精品亚洲乱码少妇综合久久| 国产极品粉嫩免费观看在线| 在线亚洲精品国产二区图片欧美| 国产熟女午夜一区二区三区| 另类精品久久| 丝袜人妻中文字幕| 亚洲欧美一区二区三区国产| 中文天堂在线官网| 国产精品嫩草影院av在线观看| 久久久国产精品麻豆| 青青草视频在线视频观看| 午夜免费观看性视频| 日韩大片免费观看网站| videos熟女内射| 精品国产乱码久久久久久小说| 久久精品久久久久久久性| 国产精品 国内视频| 亚洲四区av| 天天添夜夜摸| 午夜91福利影院| 国产一区二区激情短视频 | 国产精品三级大全| 日韩一区二区视频免费看| 国产精品亚洲av一区麻豆 | 2018国产大陆天天弄谢| av线在线观看网站| 男的添女的下面高潮视频| 国产极品粉嫩免费观看在线| 亚洲天堂av无毛| 人人澡人人妻人| 99热全是精品| 19禁男女啪啪无遮挡网站| 成人国产av品久久久| 美国免费a级毛片| 两个人免费观看高清视频| 中文字幕人妻熟女乱码| 2021少妇久久久久久久久久久| 亚洲精品一二三| 人人妻人人澡人人看| 9热在线视频观看99| 啦啦啦 在线观看视频| 可以免费在线观看a视频的电影网站 | 精品人妻在线不人妻| 欧美日韩福利视频一区二区| 婷婷成人精品国产| 制服诱惑二区| 国产精品蜜桃在线观看| 69精品国产乱码久久久| 狠狠精品人妻久久久久久综合| 国产精品一国产av| 2021少妇久久久久久久久久久| 99国产精品免费福利视频| 国产1区2区3区精品| 我的亚洲天堂| 中文字幕制服av| 肉色欧美久久久久久久蜜桃| 久久ye,这里只有精品| 自线自在国产av| 久久久久久久久久久久大奶| 亚洲,一卡二卡三卡| 巨乳人妻的诱惑在线观看| 欧美精品一区二区免费开放| 高清欧美精品videossex| 人妻一区二区av| 一本—道久久a久久精品蜜桃钙片| 青春草国产在线视频| videosex国产| www.精华液| 日韩制服丝袜自拍偷拍| 久久99精品国语久久久| 女人爽到高潮嗷嗷叫在线视频| 人妻一区二区av| 可以免费在线观看a视频的电影网站 | 99re6热这里在线精品视频| 久久久久精品国产欧美久久久 | 一级黄片播放器| 亚洲精品视频女| 狂野欧美激情性bbbbbb| 国产精品国产三级专区第一集| 韩国av在线不卡| 欧美亚洲 丝袜 人妻 在线| 一边摸一边做爽爽视频免费| 亚洲精品自拍成人| 久久久精品区二区三区| 99香蕉大伊视频| 蜜桃国产av成人99| 国产国语露脸激情在线看| 国产欧美日韩一区二区三区在线| 久久久久精品久久久久真实原创| 国产乱人偷精品视频| 色婷婷av一区二区三区视频| 国产日韩欧美在线精品| 国产激情久久老熟女| 午夜福利乱码中文字幕| 久久影院123| 黄片小视频在线播放| 亚洲国产精品一区三区| 啦啦啦 在线观看视频| 久久 成人 亚洲| 成人免费观看视频高清| 一区二区三区乱码不卡18| 最近中文字幕2019免费版| 午夜福利影视在线免费观看| 欧美在线一区亚洲| 麻豆乱淫一区二区| 婷婷色综合大香蕉| 国产乱人偷精品视频| 最新在线观看一区二区三区 | 亚洲第一青青草原| xxxhd国产人妻xxx| 国产成人精品在线电影| 不卡av一区二区三区| 欧美在线一区亚洲| 高清欧美精品videossex| 捣出白浆h1v1| 欧美老熟妇乱子伦牲交| 国产熟女午夜一区二区三区| 成年人午夜在线观看视频| 亚洲,欧美,日韩| 国产成人欧美在线观看 | 伦理电影免费视频| 国产精品麻豆人妻色哟哟久久| 美女主播在线视频| 老司机在亚洲福利影院| 人体艺术视频欧美日本| 久久青草综合色| 另类亚洲欧美激情| av又黄又爽大尺度在线免费看| 韩国精品一区二区三区| 成人手机av| 美女国产高潮福利片在线看| 大陆偷拍与自拍| 亚洲色图综合在线观看| 欧美老熟妇乱子伦牲交| 伦理电影免费视频| 日韩视频在线欧美| 亚洲综合色网址| 日韩精品免费视频一区二区三区| 国产日韩一区二区三区精品不卡| 国产一区二区激情短视频 | av免费观看日本| 天堂俺去俺来也www色官网| 午夜福利网站1000一区二区三区| 曰老女人黄片| 国产精品无大码| 久久久久久久大尺度免费视频| 日韩不卡一区二区三区视频在线| 少妇精品久久久久久久| 精品一区在线观看国产| 日韩,欧美,国产一区二区三区| 午夜日韩欧美国产| 男女下面插进去视频免费观看| 一级黄片播放器| 亚洲天堂av无毛| 99久久人妻综合| 爱豆传媒免费全集在线观看| 午夜福利一区二区在线看| 成人18禁高潮啪啪吃奶动态图| 在线天堂中文资源库| 欧美日韩国产mv在线观看视频| 制服丝袜香蕉在线| 青春草国产在线视频| 看非洲黑人一级黄片| videos熟女内射| 日本猛色少妇xxxxx猛交久久| 黄片小视频在线播放| 一二三四中文在线观看免费高清| 欧美xxⅹ黑人| 建设人人有责人人尽责人人享有的| 欧美人与性动交α欧美精品济南到| 人人妻,人人澡人人爽秒播 | 免费观看性生交大片5| 国产高清不卡午夜福利| 亚洲国产最新在线播放| 亚洲七黄色美女视频| 国产精品一区二区精品视频观看| 亚洲精品aⅴ在线观看| 无限看片的www在线观看| 国产不卡av网站在线观看| 女人久久www免费人成看片| 成人三级做爰电影| 黑人巨大精品欧美一区二区蜜桃| 欧美精品一区二区大全| 人人妻人人添人人爽欧美一区卜| 久热这里只有精品99| 国产成人欧美在线观看 | 久久人人爽av亚洲精品天堂| 青春草视频在线免费观看| 极品人妻少妇av视频| 久久久精品区二区三区| 交换朋友夫妻互换小说| 人妻人人澡人人爽人人| 大香蕉久久网| 成人漫画全彩无遮挡| 日韩精品免费视频一区二区三区| 日韩大片免费观看网站| 韩国精品一区二区三区| 欧美人与性动交α欧美精品济南到| 一边摸一边做爽爽视频免费| av网站免费在线观看视频| 美女国产高潮福利片在线看| 亚洲色图 男人天堂 中文字幕| av在线app专区| 丰满乱子伦码专区| 久久精品人人爽人人爽视色| 黄色视频在线播放观看不卡| 日韩伦理黄色片| 国产在线一区二区三区精| 日韩精品有码人妻一区| 中文字幕高清在线视频| 亚洲精品美女久久久久99蜜臀 | 国产黄频视频在线观看| 精品少妇一区二区三区视频日本电影 | 亚洲国产欧美一区二区综合| 男人爽女人下面视频在线观看| 天堂中文最新版在线下载| 婷婷色综合www| 你懂的网址亚洲精品在线观看| 日韩免费高清中文字幕av| 久久女婷五月综合色啪小说| 成人18禁高潮啪啪吃奶动态图| 波多野结衣av一区二区av| 国产黄频视频在线观看| 不卡av一区二区三区| 国产一区二区激情短视频 | 一级片免费观看大全| 亚洲国产精品成人久久小说| 久久国产精品大桥未久av| 中文字幕高清在线视频| www.av在线官网国产| 久久青草综合色| 高清视频免费观看一区二区| 搡老乐熟女国产| 亚洲精品久久成人aⅴ小说| 18禁裸乳无遮挡动漫免费视频| 国产日韩欧美亚洲二区| 亚洲一卡2卡3卡4卡5卡精品中文| 国产乱来视频区| 视频在线观看一区二区三区| 久久久久久免费高清国产稀缺| 亚洲欧洲国产日韩| 波多野结衣av一区二区av| 久久久国产一区二区| 男人添女人高潮全过程视频| 国产精品女同一区二区软件| 欧美日韩一区二区视频在线观看视频在线| 久久久精品区二区三区| 可以免费在线观看a视频的电影网站 | 大片电影免费在线观看免费| 啦啦啦 在线观看视频| 美女国产高潮福利片在线看| 日本一区二区免费在线视频| 午夜免费观看性视频| 美女福利国产在线| 视频在线观看一区二区三区| 精品人妻一区二区三区麻豆| 9色porny在线观看| 亚洲欧美精品自产自拍| 女人被躁到高潮嗷嗷叫费观| www.精华液| 成人亚洲精品一区在线观看| 亚洲精品成人av观看孕妇| 操美女的视频在线观看| 亚洲男人天堂网一区| 超碰成人久久| 欧美日韩一级在线毛片| 日韩欧美一区视频在线观看| 9191精品国产免费久久| 亚洲精品,欧美精品| 国产一区有黄有色的免费视频| 日韩成人av中文字幕在线观看| 亚洲国产精品一区二区三区在线| 中文字幕制服av| 国产精品欧美亚洲77777| 午夜福利免费观看在线| 亚洲三区欧美一区| 精品久久蜜臀av无| 大香蕉久久网| 少妇人妻精品综合一区二区| 国产成人免费观看mmmm| 亚洲在久久综合| 99久久综合免费| 大陆偷拍与自拍| 热99国产精品久久久久久7| 国产av码专区亚洲av| 日本欧美国产在线视频| 天天操日日干夜夜撸| 一边摸一边抽搐一进一出视频| 在线天堂最新版资源| a级毛片黄视频| 国产黄色免费在线视频| 妹子高潮喷水视频| 看非洲黑人一级黄片| 久久久久久久大尺度免费视频| 少妇人妻 视频| 蜜桃国产av成人99| 在线观看免费午夜福利视频| 纯流量卡能插随身wifi吗| 成人国语在线视频| 考比视频在线观看| 久久久久久久久免费视频了| 亚洲成色77777| www日本在线高清视频| 亚洲av中文av极速乱| 美女视频免费永久观看网站| 日韩制服骚丝袜av| 免费观看av网站的网址| 欧美 日韩 精品 国产| 丰满乱子伦码专区| 中文字幕另类日韩欧美亚洲嫩草| 欧美成人精品欧美一级黄| 欧美少妇被猛烈插入视频| 美女脱内裤让男人舔精品视频| 蜜桃国产av成人99| 欧美少妇被猛烈插入视频| 欧美激情高清一区二区三区 | 爱豆传媒免费全集在线观看| 国产一区二区激情短视频 | 欧美精品亚洲一区二区| 女人精品久久久久毛片| 亚洲第一青青草原| 国产成人a∨麻豆精品| 成人毛片60女人毛片免费| 99热全是精品| 精品国产乱码久久久久久小说| 亚洲av在线观看美女高潮| 2018国产大陆天天弄谢| 国产成人一区二区在线| 成人影院久久|