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

    黃土高原小麥田土壤呼吸季節(jié)和年際變化

    2013-12-20 11:27:08周小平王效科張紅星龐軍柱任玉芬侯培強(qiáng)劉文兆
    生態(tài)學(xué)報(bào) 2013年23期
    關(guān)鍵詞:年際土壤溫度土壤水分

    周小平,王效科,張紅星,* ,龐軍柱,任玉芬,侯培強(qiáng),劉文兆

    (1.中國科學(xué)院生態(tài)環(huán)境研究中心,城市與區(qū)域生態(tài)重點(diǎn)實(shí)驗(yàn)室,北京 100085;2.青海省環(huán)境科學(xué)研究設(shè)計(jì)院,西寧 810007;3.西北農(nóng)林科技大學(xué)林學(xué)院,楊陵 712100)

    過去100年間,地球表面的溫度增加了0.6℃,預(yù)計(jì)到下個(gè)世紀(jì),氣溫將會增加1.4—5.8℃,甚至更多[1]。空氣中CO2濃度在很大程度上控制著全球氣候變化[1],土壤作為僅次于海洋的全球第二大有機(jī)碳庫,其碳存儲量為1300—2000 Pg C[2],是陸地植被碳庫的2—3倍,是大氣碳庫的2倍多[3]。土壤呼吸是生態(tài)系統(tǒng)CO2交換的重要組成部分[4-5],是土壤碳輸出的主要途徑,每年約有50—75 Pg C通過土壤呼吸排放到大氣中,約占全球總排放的5%—25%,土壤呼吸的變化可能會導(dǎo)致大氣CO2濃度的較大改變,進(jìn)而影響全球氣候變化[6]。因此,土壤呼吸作為碳循環(huán)的核心問題而引起了國內(nèi)外的極大關(guān)注和廣泛研究[6-8]。

    土壤呼吸的季節(jié)變化和年際變化主要受環(huán)境因素和植物生長的影響。研究土壤呼吸的季節(jié)變化和年際變化以及調(diào)控因素,是準(zhǔn)確預(yù)測未來氣候變化條件下土壤呼吸變化的關(guān)鍵。許多研究表明,土壤溫度是決定土壤呼吸季節(jié)變化和年際變化的最重要因素[9-11],溫度升高會加速有機(jī)碳的分解,進(jìn)而促進(jìn)土壤呼吸的增加[12-16]。很多模型都是基于土壤溫度與土壤呼吸之間的指數(shù)關(guān)系來模擬土壤呼吸,但也有很多研究的結(jié)果有所不同。Grace和Rayment指出短期的增溫實(shí)驗(yàn)可能提高土壤呼吸,但是這個(gè)結(jié)果并不能代表溫度長期升高情況下土壤呼吸的變化趨勢[17]。Giardina和Ryan研究表明,僅僅增加土壤溫度并不能刺激森林土壤的C的降解[18]。也有很多研究表明,土壤呼吸的季節(jié)變化和年際變化由土壤溫度和土壤水分共同調(diào)控[19-25]。但Epron等研究表明,與土壤溫度相比,土壤水分對土壤呼吸季節(jié)變化的影響更為重要[1]。

    以往對土壤呼吸的年際變化主要通過年平均值[26]、變異系數(shù)[27]和變化范圍等來表示。但為提高全球碳循環(huán)預(yù)測的準(zhǔn)確性,僅僅研究土壤呼吸年際變化是不夠的,還需要對土壤呼吸年際變化的成因進(jìn)行分析。Hui等[28]用斜率同質(zhì)性模型將土壤呼吸年際變化的成因分為4種:(1)氣候因子年際變化對土壤呼吸的直接影響;(2)氣候因子季節(jié)變化對土壤呼吸的直接影響;(3)氣候因子通過影響生物生態(tài)過程而間接對土壤呼吸的調(diào)控(FC);(4)隨機(jī)誤差。Hui對森林生態(tài)系統(tǒng)呼吸年際變化的分析表明,功能變化對土壤呼吸的影響要大于環(huán)境因子年際變化對土壤呼吸年際變化的直接影響[28]。

    農(nóng)田是我國最重要的生態(tài)系統(tǒng)之一,我國農(nóng)田生態(tài)系統(tǒng)土壤呼吸連續(xù)多年原位監(jiān)測的研究較少,缺少對農(nóng)田生態(tài)系統(tǒng)土壤呼吸年際變化及其成因的分析。因此,本研究對黃土高原小麥田生態(tài)系統(tǒng)CO2交換進(jìn)行了4a的連續(xù)觀測,分析了土壤呼吸的日變化、季節(jié)變化和年際變化特征,分析了主要環(huán)境因子(土壤溫度、空氣溫度、土壤濕度、降水量)對土壤呼吸的影響,并對土壤呼吸年際變化的成因進(jìn)行了分解,旨在為準(zhǔn)確預(yù)測未來氣候變化下的土壤呼吸變化提供數(shù)據(jù)支持和理論依據(jù)。

    1 材料與方法

    1.1 研究區(qū)概況

    本研究場地布設(shè)于中國科學(xué)院長武農(nóng)業(yè)生態(tài)試驗(yàn)站,該站位于黃土高原中南部的陜西甘肅交界處,長武縣東王村(北緯35°12',東經(jīng)107°40'),海拔1200 m,地下水位60 m以下。該地區(qū)農(nóng)業(yè)生產(chǎn)全部依靠自然降水,屬典型的旱作農(nóng)業(yè)區(qū),作物為一年一熟,主要有冬小麥(Triticum aestivum L.)、春玉米(Zea mays L.)、糜子(Panicum miliaceum L.)等。該地區(qū)屬暖溫帶半濕潤大陸型季風(fēng)氣候,1957—2009年年平均降水量為580.5 mm,其中休閑期降水量占年降水量的54%,降水分布不均,年際間差異大,干旱頻繁。年均氣溫9.1℃,無霜期171d。土壤屬黑壚土,土壤有機(jī)質(zhì)含量約3%,pH值為8.4,粘粒(<0.002 mm)含量為24%,土質(zhì)均勻疏松。

    1.2 研究材料及管理方法

    供試冬小麥(Triticum aestivum L.)品種為長旱58,是黃土高原典型的小麥品種。冬小麥9月中旬播種,翌年4月中旬人工鋤草,6月中旬用聯(lián)合收割機(jī)收獲。播種時(shí)同時(shí)施肥,施肥品種主要為尿素和磷肥,施肥標(biāo)準(zhǔn)為:尿素(含氮量為46%)300 kg/hm2,過磷酸鈣750 kg/hm2。本研究區(qū)為雨養(yǎng)農(nóng)業(yè)區(qū),無灌溉措施。小麥?zhǔn)斋@后機(jī)械翻耕,深度為30 cm。

    本研究時(shí)間范圍為2005年7月1日—2009年6月30日,經(jīng)歷4個(gè)小麥生長期。在年際變化研究時(shí),每年的時(shí)間劃分為:2005年7月1日—2006年6月30日(以下簡寫為2005—2006年),2006年7月1日—2007年6月30日(以下簡寫為2006—2007年),2007年7月1日—2008年6月30日(以下簡寫為2007—2008年),2008年7月1日—2009年6月30日(以下簡寫為2008—2009年)。在研究期間,小麥田管理措施見表1。

    表1 黃土高原冬小麥管理措施Table 1 Schedule of management activities in a winter wheat ecosystem in the Loess Plateau,China

    1.3 CO2通量的測定方法

    本研究采用多通道通量箱法對CO2通量進(jìn)行測定[29]。多通道全自動通量箱系統(tǒng)由兩個(gè)子系統(tǒng)組成:CO2濃度分析和記錄子系統(tǒng)及全自動通量箱子系統(tǒng)(圖1)。兩個(gè)子系統(tǒng)都在一臺可編程邏輯控制器(PLC,Programmable Logic Controller)控制下工作。全自動通量箱子系統(tǒng)由8個(gè)全自動通量箱和1臺空氣壓縮機(jī)組成。全自動通量箱為長、寬、高均為50 cm的正方體,框架由鋁合金制成。箱體四壁采用透明的亞克力板(透光率98%),并通過雙面密封膠條和鋁合金框架粘接起來用螺釘固定。箱蓋采用較厚的亞克力板(透光率98%),與箱體間用不銹鋼合頁連接。為了防止箱蓋熱脹冷縮變形,箱蓋外表面固定有高強(qiáng)度U型鋁合金板。箱體與箱蓋間鑲有高密度密封條以提高氣密性。箱蓋中部位置和箱體的一側(cè)中部位置之間固定一個(gè)適當(dāng)長度的氣缸。當(dāng)氣缸完全伸展時(shí),箱蓋完全打開,保持與外界很好的氣體交換;當(dāng)氣缸回縮后,箱蓋完全關(guān)閉,進(jìn)行通量的測定。氣缸往返動作由PLC控制的電磁閥控制,由空氣壓縮機(jī)的壓縮空氣驅(qū)動。在氣缸缸體上固定著1個(gè)風(fēng)扇,用于通量測定時(shí)進(jìn)行箱內(nèi)氣體的混合。通量箱蓋關(guān)閉時(shí),風(fēng)扇開始運(yùn)轉(zhuǎn)以充分?jǐn)嚢柰肯渲械臍怏w。箱體的一側(cè)靠上位置安裝有出氣管,對側(cè)的靠下位置安裝進(jìn)氣管,分別與CO2分析子系統(tǒng)相連。在通量箱的頂部接1根長1.5 m的塑料管(4 mm×6 mm)以平衡通量箱內(nèi)外氣壓差。

    CO2濃度分析和記錄子系統(tǒng)由紅外 CO2分析儀(LI-820,LI-COR Inc.,Lincoln,Nebraska)、緩沖管、干燥器、過濾器、流量計(jì)、氣泵、多通閥和數(shù)據(jù)采集器組成。氣泵將通量箱中的氣樣抽出,并經(jīng)過多通閥、緩沖管、干燥器、過濾器和流量計(jì)后,進(jìn)入紅外CO2分析儀,測定通量箱中的空氣CO2濃度,分析后的氣體經(jīng)過另外一個(gè)多通閥流回通量箱。每個(gè)箱子測定時(shí)間為3 min[30],CO2分析儀連接在數(shù)據(jù)采集器上(CR10,Campbell Scientific Inc.,Logan,UT,USA)。數(shù)據(jù)采集器設(shè)定為每10s記錄1次紅外CO2分析儀的讀數(shù)。

    圖1 多通道全自動通量箱系統(tǒng)結(jié)構(gòu)示意圖Fig.1 Multi-channel automated chamber system

    用透明通量箱系統(tǒng)來測定土壤CO2通量時(shí),在小麥田行間放置通量箱,小麥出苗后,去除箱子內(nèi)部的小麥,保證通量箱內(nèi)沒有小麥植株,箱體內(nèi)部的氣體交換僅為土壤呼吸。

    本研究以通量箱為重復(fù),共設(shè)4個(gè)重復(fù)。

    1.4 環(huán)境因子的測定方法

    降水量和大氣壓數(shù)據(jù)來源于中國科學(xué)院長武農(nóng)業(yè)生態(tài)試驗(yàn)站氣象站。土壤溫度(5cm)和空氣溫度采用T型熱電偶(T型)測定;土壤水分采用EC-H2O傳感器測定(測定土壤水分為0—30cm土壤的體積含水量),T型熱電偶和EC-H2O傳感器均連接在數(shù)據(jù)采集器CR10x上進(jìn)行數(shù)據(jù)記錄和保存。數(shù)據(jù)采集頻率為1次/10min。

    1.5 數(shù)據(jù)處理方法

    1.5.1 原始數(shù)據(jù)分箱及CO2通量計(jì)算

    將數(shù)據(jù)采集器里的數(shù)據(jù)下載到計(jì)算機(jī)中,用C++語言編寫計(jì)算機(jī)程序,將整塊數(shù)據(jù)分到各對應(yīng)的通量箱,然后轉(zhuǎn)化為Excel形式,在此形式下,利用下式計(jì)算CO2通量:

    式中,A是單位面積上單位時(shí)間內(nèi)CO2釋放量(μmol·m-2·s-1),c是CO2摩爾濃度(μmol/mol);t是時(shí)間(s);v是通量箱體積(m3);S是通量箱底面積(m2);P是大氣壓(kPa);R是氣體常數(shù)(8.3×10-3m3·kPa·mol-1·k-1);T是通量箱內(nèi)氣體溫度(K)。方程(1)中d c/d t是通量箱中氣體CO2濃度變化率,即為所測得的一組CO2濃度及其相應(yīng)的時(shí)間回歸所得直線方程的斜率。該直線方程的R2<0.85時(shí),舍棄擬合所得的斜率值。

    1.5.2 土壤呼吸IAV的成因分析

    為研究引起土壤呼吸年際變化的因素,采用HOS模型對土壤呼吸年際變化的成因進(jìn)行分析[28]。將引起土壤呼吸年際變化的原因分為4種:(1)環(huán)境因子對土壤呼吸的間接影響(功能變化),(2)環(huán)境因子年際變化對土壤呼吸的直接影響,(3)環(huán)境因子季節(jié)變化對土壤呼吸的直接影響,(4)隨機(jī)誤差。

    HOS模型是一種基于回歸分析的模型,該模型用來檢驗(yàn)環(huán)境因子和年之間的交互作用對土壤呼吸的影響:

    式中,Yij表示土壤呼吸,i表示第 i年,i=1,2,…,y(本研究中,y=4);j表示第 j天,j=1,2,…,n;Yij表示第 i年第j天的土壤呼吸,k表示第k個(gè)因素,k=1,2,…,m(本研究中,m=4,是土壤溫度、空氣溫度、土壤濕度和降水量)。

    其假設(shè)為:各環(huán)境因素與土壤呼吸的關(guān)系斜率在各年際間均相同。即不存在功能變化(bik=0)。用F檢驗(yàn)對每年使用不同的斜率與所有年使用相同的斜率相比較,是否能顯著提高方程(2)的擬合度。如果原假設(shè)H0被接受,則表明斜率的年際間差異并不顯著,認(rèn)為功能變化不存在,方程2可以簡化為多元線性回歸模型:

    如果原假設(shè)H0被拒絕,則表明斜率有顯著的年際間差異,認(rèn)為功能變化存在,方程2則被簡化為多元線性回歸模型:

    由方程3和方程4模擬出來的因變量的差值即為功能變化。

    通過用統(tǒng)計(jì)學(xué)的方法將因變量變異的平方和進(jìn)行分解,從而實(shí)現(xiàn)土壤呼吸年際變化的成因分解:

    式中,SST為土壤呼吸的變異平方和,SSf為環(huán)境因素對土壤呼吸的間接影響(功能變化)所能解釋的土壤呼吸的變異平方和,SSic為由環(huán)境因素年際變化引起的因變量的變異平方和,SSsc為由環(huán)境因素的季節(jié)變化引起的因變量的變異平方和,SSe為由系統(tǒng)誤差引起的因變量的變異平方和。通過對SST的分解,從而得出不同成因?qū)ν寥篮粑觌H變化的貢獻(xiàn)率。如果HOS模型的檢驗(yàn)結(jié)果表明FC并不顯著,則將土壤呼吸的變異分解為3部分,環(huán)境因素的年際變化,環(huán)境因素的季節(jié)變化,以及隨機(jī)誤差。

    用SPSS 16.0進(jìn)行所有的統(tǒng)計(jì)分析。用sigmaplot 10.0進(jìn)行圖形繪制。

    2 結(jié)果與分析

    2.1 土壤呼吸和各環(huán)境因素的變化

    2.1.1 土壤呼吸的日變化

    除冬天外,土壤呼吸日變化呈顯著的單峰曲線(圖2)。土壤呼吸日變化的主要驅(qū)動因素是土壤溫度的日變化。夜間(0:00—5:00)土壤溫度的變化較大,但土壤呼吸夜間的變化幅度較小。從7:00左右開始,土壤呼吸隨著土壤溫度的增加而增加,在12:00左右達(dá)到最大值,不同季節(jié)土壤溫度和土壤呼吸達(dá)到最大值的時(shí)間不同。冬季土壤溫度在0℃以下,雖然土壤溫度的變化幅度較大,但土壤呼吸變化幅度很小。

    2.1.2 土壤呼吸和各環(huán)境因素的季節(jié)變化

    土壤溫度、空氣溫度均表現(xiàn)出明顯的季節(jié)變化,由于季風(fēng)型氣候的影響,降水量和土壤含水量也表現(xiàn)出明顯的季節(jié)變化(圖3)。溫度和水分的季節(jié)變化趨勢較為一致,6月份至9月份是全年中溫度最高的季節(jié),也是降水量集中分布的季節(jié),占全年降水量的60%。土壤含水量的季節(jié)變化趨勢與降水量季節(jié)變化趨勢一致。

    各環(huán)境因子的年平均值年際間變化較大(表2)。年平均空氣溫度的變化幅度為11.0—9.5℃,年平均土壤溫度的變化幅度為11.7—10.5℃,年平均土壤濕度的變化幅度為28.0%—19.2%,年總降水量變化幅度為568.3—495.9 mm。

    表2 每年土壤呼吸和環(huán)境因子年平均值Table 2 Mean annual values of soil temperature,air temperature,precipitation and soil water content

    圖2 土壤呼吸和土壤溫度(10cm)的典型日變化Fig.2 The diurnal variations of R soil and T soil

    土壤呼吸表現(xiàn)出明顯的逐日變化和季節(jié)變化(圖4)。土壤呼吸在冬季較低,最小值大概為0.25 μmol·m-2·s-1,從冬季到7月下旬逐漸增加,7月下旬達(dá)到最大值(17.92μmol·m-2·s-1),土壤呼吸年際間變化較大(表2),其變化幅度為815.72—980.12 g C·m-2·a-1。土壤呼吸極大值出現(xiàn)在小麥?zhǔn)崭詈?,耕地?dǎo)致了土壤呼吸的大量增加[31].

    2.2 土壤呼吸與環(huán)境因子的關(guān)系

    對土壤呼吸與環(huán)境因子的關(guān)系進(jìn)行了回歸分析(圖5)。土壤溫度、空氣溫度以及土壤濕度都對土壤呼吸有顯著的影響(P<0.0001),其中土壤溫度對土壤呼吸的影響最大(r2=0.40),土壤含水量與土壤呼吸之間有較弱的正相關(guān)關(guān)系(r2=0.12),降水量與土壤呼吸之間的相關(guān)關(guān)系不顯著。土壤溫度和空氣溫度與土壤呼吸之間的關(guān)系也可以用指數(shù)關(guān)系來描述,回歸系數(shù)略高于線性回歸。

    多元逐步回歸分析結(jié)果表明,用土壤濕度和土壤溫度可以很好的擬合土壤呼吸。逐步回歸分析得到的最優(yōu)擬合方程為:Rsoil=1.761+0.119Tsoil- 1.30SWC,R2=0.48。用通徑分析對土壤溫度和土壤濕度對土壤呼吸的影響進(jìn)行了分析,土壤溫度對土壤呼吸的通徑系數(shù)為0.73,SWC對土壤呼吸通徑系數(shù)為0.07,表明土壤溫度對土壤呼吸的影響要大于土壤濕度的影響。盡管逐步回歸模型可以很好的模擬土壤呼吸,但仍然有很52%的土壤呼吸的變異是由于其他因素引起的(圖6)。

    2.3 土壤呼吸年際變化成因分析

    土壤呼吸年際間變化較大,其變化幅度為815.72—980.12 g C·m-2·a-1。2005—2006年土壤呼吸最低(815.72 g C·m-2·a-1),2007—2008年土壤呼吸最高(980.12 g C·m-2·a-1)。

    逐步回歸分析中,土壤濕度和土壤水分都顯著影響土壤呼吸。用HOS模型對Tsoil和SWC對土壤呼吸年際變化的影響進(jìn)行分析。結(jié)果表明,Tsoil和Tsoil×year顯著影響土壤呼吸,而SWC和SWC×year與土壤呼吸的關(guān)系不顯著,最終HOS模型中,Tsoil和Tsoil×year兩個(gè)變量進(jìn)入了模型(表3)。表明Tsoil對土壤呼吸的影響存在顯著的年際間差異,F(xiàn)C存在。用4a同一斜率的Tsoil-Rsoil直線關(guān)系模擬的土壤呼吸的模擬結(jié)果和每年不同斜率Tsoil-Rsoil直線關(guān)系模擬的土壤呼吸相比較發(fā)現(xiàn),每年不同斜率模擬提高了模擬效果,R2從0.44增加到0.56(圖 5)。平方和的分解結(jié)果表明,SSf、SSic、SSsc和 SSe對土壤呼吸變異的貢獻(xiàn)分別為 10.6%、4.5%、58.4%、和26.5%。

    3 討論

    在日尺度上,土壤呼吸主要受溫度變化的調(diào)控。分析發(fā)現(xiàn),用指數(shù)關(guān)系可以很好的描述土壤呼吸與土壤溫度日變化的關(guān)系。以2008年7月份為例(圖7,每小時(shí)的土壤呼吸和土壤溫度為2008年7月所有晴天該小時(shí)值的平均值),用指數(shù)關(guān)系回歸土壤溫度和土壤呼吸,其R2為0.84,但同一土壤溫度下,上午土壤呼吸要高于下午(圖7)。

    圖4 土壤呼吸的季節(jié)變化Fig.4 Seasonal variation of soil respiration(R soil)

    圖5 土壤呼吸與空氣溫度、土壤溫度和土壤含水量的關(guān)系Fig.5 Simple relationships between soil respiration(R soil)and air temperature(T air),soil temperature(T soil)and soil water content

    圖6 土壤呼吸的實(shí)測值與模擬值比較Fig.6 Comparisons of model estimated by the multiple regression model

    表3 土壤呼吸的HOS模型的ANOVA結(jié)果Table 3 An ANOVA of a homogeneity-of-slopes model of soil respiration(R soil)

    圖7 一天內(nèi)土壤溫度和土壤呼吸小時(shí)值的關(guān)系Fig.7 Relationships between hourly soil respiration and soil temperature and the relationships that before and after peak point was described with different slopes

    韓廣軒等指出,環(huán)境因子(水分、溫度等)、作物生物學(xué)特性和農(nóng)業(yè)管理活動是造成中國農(nóng)田生態(tài)系統(tǒng)土壤呼吸作用時(shí)間變異的主要因素[32]。土壤溫度和土壤水分是土壤呼吸季節(jié)變化的最主要環(huán)境因子,但不同的生態(tài)系統(tǒng)中,研究結(jié)果有所不同。許多研究中,土壤呼吸的季節(jié)變化都是用土壤溫度的變化來解釋[11-13],也有一些土壤呼吸季節(jié)變化是用土壤溫度和土壤水分的共同變化來解釋的[14-19]。在熱帶森林中,土壤水分對土壤呼吸的季節(jié)變化的影響比土壤溫度更大[2,33]。本研究中,土壤溫度的變化范圍較大(33℃),土壤溫度和土壤水分的季節(jié)變化是土壤呼吸季節(jié)變化的最重要原因,與土壤呼吸呈顯著的正相關(guān)關(guān)系。土壤呼吸和土壤水分的關(guān)系通常用直線關(guān)系來描述[24,26,34],本研究中,土壤水分和土壤呼吸的關(guān)系較弱。降水規(guī)律、年總降水量以及土壤結(jié)構(gòu)等的差異可能是導(dǎo)致這種不同結(jié)果的原因[33]。在土壤含水量較高或土壤粘度較大的地區(qū),土壤水分和土壤呼吸之間可能是單峰的關(guān)系,在土壤水分低于限制值時(shí)以及土壤排水性好的地區(qū),因?yàn)閰捬跻种频那闆r極少可能發(fā)生,土壤呼吸隨著土壤水分的增加而增加,土壤含水量與土壤水分之間可能會是指數(shù)或者直線關(guān)系[33,35]。但當(dāng)土壤水分含量高于限制值時(shí),由于土壤處于厭氧條件而抑制了土壤呼吸[23,25]。本研究區(qū)土壤含水量較低,因此采用直線關(guān)系來描述土壤呼吸和土壤含水量之間的關(guān)系。本研究中,土壤溫度對土壤呼吸的影響要大于土壤水分,通徑分析的結(jié)果表明,土壤溫度對土壤呼吸的通徑系數(shù)為0.73,SWC對土壤呼吸通徑系數(shù)0.06。指數(shù)關(guān)系對土壤呼吸和土壤溫度的關(guān)系的描述要優(yōu)于線性關(guān)系。

    土壤呼吸年際變化是普遍存在的[35],影響土壤呼吸年際變化的主要因素是環(huán)境因子。不同生態(tài)系統(tǒng)土壤呼吸年際變化的影響因子不同,Gaumont-Guay等的研究結(jié)果表明,干旱對土壤呼吸年際變化的貢獻(xiàn)大于土壤溫度對土壤呼吸年際變化的貢獻(xiàn)[26]。本研究結(jié)果有所不同,在本研究中土壤溫度是引起土壤呼吸年變量的最重要因素。環(huán)境因子對土壤呼吸的影響包括直接影響和間接影響。間接影響是指環(huán)境因子通過影響生物生態(tài)過程而對土壤呼吸的調(diào)控,即功能變化,比如,溫度變化引起土壤微生物群落和微生物活性的變化[36,37],進(jìn)而影響到土壤呼吸的變化。功能變化的存在,使得每一年氣候因子與土壤呼吸之間的關(guān)系有了很大的差異。很多研究發(fā)現(xiàn)了這一現(xiàn)象[27]。Chen的研究表明,土壤呼吸的對溫度的敏感性年際間差異顯著[38]。HOS模型為對土壤呼吸的年際變化的成因分析提供了一個(gè)有效的方法。本研究用HOS模型對土壤呼吸的年際變化成因進(jìn)行分析發(fā)現(xiàn),土壤含水量并不是引起本研究土壤呼吸年際變化的重要因子,功能變化、土壤溫度季節(jié)變化、土壤溫度年際變化以及隨機(jī)誤差分別可以解釋年際變化的10.6%、58.4%、4.5%、26.5%,功能變化的對土壤呼吸年際變化的影響要大于土壤溫度的年際變化的影響,因此,在準(zhǔn)確預(yù)測氣候變化下土壤呼吸的變化,僅僅考慮到氣候因子的變化是不夠的,還要考慮到功能變化對土壤呼吸的影響。

    4 結(jié)論

    本研究通過對黃土高原小麥田土壤呼吸的長期觀測,研究了土壤呼吸的季節(jié)變化和年變化及其影響因素。結(jié)果表明,土壤呼吸表現(xiàn)出顯著的季節(jié)變化和年際變化。土壤呼吸的季節(jié)變化主要受土壤溫度和土壤水分的調(diào)控,土壤溫度和土壤含水量均與土壤呼吸呈正相關(guān)關(guān)系,土壤溫度對土壤呼吸的影響要大于土壤含水量??刂仆寥篮粑哪觌H變化的環(huán)境因子是土壤溫度。土壤溫度與土壤呼吸的關(guān)系年際間變化顯著。環(huán)境因子年際變化對土壤呼吸的年際變化的間接影響要大于其直接影響。因此,在對未來氣候變化下土壤呼吸進(jìn)行預(yù)測時(shí),我們應(yīng)充分考慮到環(huán)境因子對土壤呼吸的直接影響,同時(shí)也必須考慮到環(huán)境因子通過影響生物和生態(tài)過程而對土壤呼吸的間接影響,不能通過短期的土壤呼吸與環(huán)境因子建立的模型對未來的土壤呼吸進(jìn)行直接預(yù)測。本研究表明,未來對土壤呼吸的過程模型中要對土壤呼吸進(jìn)行長期的預(yù)測時(shí),為提高模型的準(zhǔn)確性,必須在模型中充分考慮到功能變化的影響。

    [1]Houghton J T,Ding Y J,et al.Intergovernmental Panel on Climate Change:2001//Climate Change:The Scientific Basis.Cambridge:Cambridge University Press,2001

    [2]Jenkison D S,Adams D E,Wild A.Model estimates of CO2emissions from soil in response to global warming.Nature,1991,351(6324):304-306.

    [3]Bolin B,Degens E T.The global biogeochemical carbon cycle.Geophysical Research Letters,2001,8:555-558.

    [4]Granier A,Ceschia E,Damesin C,Dufrêne E,Epron D,Gross P,Lebaube S,Le-Dantec V,Le-Goff N,Lemoine D,Lucot E,Ottorini J M,Pontailler J Y,Saugier B.The carbon balance of a young beech forest.Functional Ecology,2000,14(3):312-325.

    [5]Janssens I A,Lankreijer H,Metteucci G,Kowalski A S,Buchmann N,Epron D,Pilegaard K,Kutsch W,Longdoz B,Grunwald T,Montagnani L,Dore S,Rebmann C,Moors E J,Grelle A,Rannik ü,Morgenstern K,Oltchev S,Clement R,Guemundsson J,Minerbi S,Berbigier P,Ibrom A,Moncrieff J,Aubinet M,Bernhofer C,Jensen N O,Vesala T,Granier A,Schulze E D,Lindroth A,Dolman A J,Jarvis P G,Ceulemans R,Valentini R.Productivity overshadows temperature in determining soil and ecosystem respiration across European forests.Global Change Biology,2001,7(3):269-278.

    [6]Liu SH,F(xiàn)ang J Y.Effect factors of soil respiration and the temperature's effects on soil respiration in the global scale.Acta Ecologica Sinica,1997,117(5):469-476.

    [7]Schlesinger W H,Andrews JA.Soil respiration and the global carbon cycle.Biogeochemistry,2000,48(1):7-20.

    [8]Singh J S,Gupta SR.Plant decomposition and soil respiration in terrestrial ecosystems.The Botanical Review,1997,43(4):449-528.

    [9]Anderson J M.Carbon dioxide evolution from two temperate,deciduous woodland soils.The Journal of Applied Ecology,1973,10(2):361-375.

    [10]Edwards N T,Sollins P.Continuous measurement of carbon dioxide evolution from partitioned forest floor components.Ecology,1975,54(2):406-412.

    [11]Longdoz B,Yernaux M,Aubinet M.Soil CO2efflux measurements in a mixed forest:impact of chamber disturbances,spatial variability and seasonal evolution.Global Change Biology,2000,6(8):907-917.

    [12]Cao M K,Woodward FI.Dynamic responses of terrestrial ecosystem carbon cycling to global climate change.Nature,1998,393(6682):249-252.

    [13]Kirschbaum M U F.The temperature dependence of soil organic matter decomposition,and the effect of global warming on soil organic C storage.Soil Biology and Biochemistry,1995,27(6):753-760.

    [14]Schimel D S,Braswell B H,Holland E A,Mckeown R,Ojima D S,Painter T H,Parton W J,Townsend A R.Climatic,edaphic,and biotic controls over storage and turnover of carbon in soils.Global Biogeochemical Cycles,1994,8(3):279-293.

    [15]Trumbore S E,Chadwick O,Amundson R.Rapid exchange between soil carbon and atmospheric carbon dioxide driven by temperature change.Science,1996,272(5260):393-396.

    [16]VEMAPM Members.Vegetation/ecosystem modeling and analysis project:Comparing biogeography and biogeochemistry models in a continentalscale study of terrestrial ecosystem responses to climate change and CO2doubling.Global Biogeochemical Cycles,1995,9(4):407-437.

    [17]Grace J,Rayment M.Respiration in the balance.Nature,2000,404(6780):819-820.

    [18]Giardina CP,Ryan M G.Evidence that decomposition rates of organic carbon in mineral soil do not vary with temperature.Nature,2000,404(6780):858-861.

    [19]Garret H E,Cox GS.Carbon dioxide evolution from the floor of an oak-h(huán)ickory forest.Soil Science Society of America Journal,1973,37(4):641-644.

    [20]Hanson P J,Edwards N T,Garten CT,Andrews JA.Separating root and soil microbial contributions to soil respiration:a review of methods and observations.Biogeochemistry,2000,48(1):115-146.

    [21]Keith H,Jacobsen K L,Raison R J.Effects of soil phosphorus availability,temperature and moisture on soil respiration in Eucalyptus pauciflora forest.Plant and Soil,1997,190(1):127-141.

    [22]Davidson E A,Belk E,Boone R D.Soil water content and temperature as independent or confounded factors controlling soil respiration in a temperate mixed hardwood forest.Global Change Biology,1998,4(2):217-227.

    [23]Epron D,F(xiàn)arque L,Lucoté,Badot P M.Soil CO2efflux in a beech forest:dependence on soil temperature and soil water content.Annals of Forest Science,1999,56(3):221-226.

    [24]Qi Y,Xu M.Separating the effects of moisture and temperature on soil CO2efflux in a coniferous forest in the Sierra Nevada Mountains.Plant and Soil,2001,237(1):15-23.

    [25]Rey A,Pegoraro E,Tedeschi V,De Parri I,Jarvis P G,Valentini R.Annual variation in soil respiration and its components in a coppice oak forest in Central Italy.Global Change Biology,2002,8(9):851-866.

    [26]Gaumont-Guay D,Black T A,Griffis T J,Barr A G,Morgenstern K,Jassal R S,Nesic Z.Influence of temperature and drought on seasonal and interannual variations of soil,bole and ecosystem respiration in a boreal aspen stand.Agricultural and Forest Meteorology,2006,140(1/4):203-219.

    [27]Mo W H,Lee M S,Uchida M,Inatomi M,Saigusa N,Mariko S,Koizumi H.Seasonal and annual variations in soil respiration in a cool-temperate deciduous broad-leaved forest in Japan.Agricultural and Forest Meteorology,2005,134(1/4):81-94.

    [28]Hui D F,Luo Y Q,Katul G.Partitioning interannual variability in net ecosystem exchange between climatic variability and functional change.Tree Physiology,2003,23(7):433-442.

    [29]Zhang H X,Wang X K,F(xiàn)eng Z W,Song W Z,Liu W Z,Ouyang Z Y.Multi-channel automated chamber system for continuously monitoring CO2exchange between agro-ecosystem or soil and the atmosphere.Acta Ecologica Sinica,2007,27(4):1273-1282.

    [30]Drewitta G B,Blacka T A,Nesic Z,Humphreys E R,Jork E M,Swanson R,Ethier G J,Griffis T,Morgenstern K.Measuring forest floor CO2fluxes in a Douglas-fir forest.Agricultural and Forest Meteorology,2002,110(4):299-317.

    [31]Zhang H X,Zhou X P,Lu F,Pang JZ,F(xiàn)eng Z W,Liu W Z,Ouyang Z Y,Wang X K.Seasonal dynamics of soil CO2efflux in a conventional tilled wheat field of the Loess Plateau,China.Ecology Research,2011,26(4):735-743.

    [32]Han G X,Zhou GS,Xu ZZ.Research and prospects for soil respiration of farmland ecosystems in China.Journal of Plant Ecology,2008,32(3):719-733.

    [33]Kosugi Y,Mitani T,Itoh M,Noguchi S,Tani M,Matsuo N,Takanashi S,Ohkubo S,Rahim N A.Spatial and temporal variation in soil respiration in a Southeast Asian tropical rainforest.Agricultural and Forest Meteorology,2007,147(1/2):35-47.

    [34]Holt J A,Hodgen M J,Lamb D.Soil respiration in the seasonally dry tropics near townsville,north-Queensland.Australian Journal of Soil Research,1990,28(5):737-745.

    [35]Davidson E A,Verchot L V,Cattanio J H,Ackerman I L,Carvalho J E M.Effects of soil water content on soil respiration in forests and cattle pastures of eastern Amazonia.Biogeochemistry,2000,48(1):53-69.

    [36]Andrews RT,Bromley P J,Pfister M E.Successful embolization of collaterals from the ovarian artery during uterine artery embolization for fibroids:a case report.Journal of Vascular and Interventional Radiology,2000,11(5):607-610.

    [37]Janssens I A,Pilegaard K I M.Large seasonal changes in Q10of soil respiration in a beech forest.Global Change Biology,2003,9(6):911-918.

    [38]Chen B Y,Liu SR,Ge JP,Chu JX.Annual and seasonal variations of Q10soil respiration in the sub-alpine forests of the Eastern Qinghai-Tibet Plateau,China.Soil Biology and Biochemistry,2010,42(10):1735-1742.

    參考文獻(xiàn):

    [6]劉紹輝,方精云.土壤呼吸的影響因素及全球尺度下溫度的影響.生態(tài)學(xué)報(bào),1997,17(5):469-476.

    [29]張紅星,王效科,馮宗煒,宋文質(zhì),劉文兆,歐陽志云.用于測定陸地生態(tài)系統(tǒng)與大氣間CO2交換通量的多通道全自動通量箱系統(tǒng).生態(tài)學(xué)報(bào),2007,27(4):1273-1282.

    [32]韓廣軒,周廣勝,許振柱.中國農(nóng)田生態(tài)系統(tǒng)土壤呼吸作用研究與展望.植物生態(tài)學(xué)報(bào),2008,32(3):719-733.

    猜你喜歡
    年際土壤溫度土壤水分
    北緯30°中層頂區(qū)域鈉與鐵原子層的結(jié)構(gòu)和年際變化
    遼東山區(qū)3種人工林土壤呼吸對土壤溫度和土壤水分的響應(yīng)
    西藏高原土壤水分遙感監(jiān)測方法研究
    管群間歇散熱的土壤溫度響應(yīng)與恢復(fù)特性
    亞洲夏季風(fēng)的年際和年代際變化及其未來預(yù)測
    與北大西洋接壤的北極海冰和年際氣候變化
    土壤濕度和土壤溫度模擬中的參數(shù)敏感性分析和優(yōu)化
    不同覆蓋措施對棗園土壤水分和溫度的影響
    基于M-K法對圖們江干流含沙量年際變化的分析
    植被覆蓋區(qū)土壤水分反演研究——以北京市為例
    狂野欧美激情性bbbbbb| 久久青草综合色| 91精品一卡2卡3卡4卡| 久久国产精品男人的天堂亚洲 | 三级经典国产精品| 麻豆成人av视频| 伊人久久精品亚洲午夜| 18+在线观看网站| 日韩欧美一区视频在线观看 | 成人亚洲欧美一区二区av| av福利片在线观看| 美女脱内裤让男人舔精品视频| 国产有黄有色有爽视频| 午夜福利网站1000一区二区三区| 亚洲av中文av极速乱| 美女主播在线视频| 国产亚洲最大av| 丰满少妇做爰视频| 青春草国产在线视频| 日韩精品免费视频一区二区三区 | 高清毛片免费看| 内射极品少妇av片p| 久久久亚洲精品成人影院| 99久久中文字幕三级久久日本| 啦啦啦视频在线资源免费观看| 三级国产精品片| 少妇的逼水好多| 精品熟女少妇av免费看| 久久久久久久久大av| 人人妻人人添人人爽欧美一区卜| 久久精品国产亚洲网站| 国产淫语在线视频| 99视频精品全部免费 在线| 国产欧美日韩综合在线一区二区 | 国产日韩欧美在线精品| 岛国毛片在线播放| 深夜a级毛片| 我要看日韩黄色一级片| 最近2019中文字幕mv第一页| 国产色爽女视频免费观看| 最近中文字幕2019免费版| 视频区图区小说| 噜噜噜噜噜久久久久久91| 热re99久久国产66热| 精品午夜福利在线看| a级片在线免费高清观看视频| 乱码一卡2卡4卡精品| av在线老鸭窝| 美女中出高潮动态图| 午夜免费男女啪啪视频观看| 欧美日韩一区二区视频在线观看视频在线| 日韩av在线免费看完整版不卡| 自拍偷自拍亚洲精品老妇| 久久99精品国语久久久| 成人黄色视频免费在线看| 中国三级夫妇交换| 国产视频首页在线观看| 国模一区二区三区四区视频| 欧美3d第一页| 国产爽快片一区二区三区| 自线自在国产av| 日韩电影二区| 人妻制服诱惑在线中文字幕| 色网站视频免费| 中文天堂在线官网| 少妇人妻一区二区三区视频| 五月伊人婷婷丁香| 国产亚洲91精品色在线| 日本猛色少妇xxxxx猛交久久| 日本爱情动作片www.在线观看| 黄色视频在线播放观看不卡| 少妇 在线观看| av在线老鸭窝| 午夜免费男女啪啪视频观看| 久久久久久久久久成人| 欧美日本中文国产一区发布| 全区人妻精品视频| 国产av码专区亚洲av| 夜夜看夜夜爽夜夜摸| 国产av码专区亚洲av| 亚洲国产精品一区二区三区在线| 日韩av在线免费看完整版不卡| 亚洲怡红院男人天堂| 七月丁香在线播放| 国产成人免费观看mmmm| 日韩不卡一区二区三区视频在线| 大码成人一级视频| xxx大片免费视频| 久久午夜综合久久蜜桃| 国产精品女同一区二区软件| 高清不卡的av网站| 黄色怎么调成土黄色| 嫩草影院入口| 精品一区二区三区视频在线| 国产精品久久久久久久久免| 乱码一卡2卡4卡精品| 国产色婷婷99| www.av在线官网国产| 寂寞人妻少妇视频99o| 中文欧美无线码| 亚洲欧美成人综合另类久久久| 久久久精品94久久精品| 少妇裸体淫交视频免费看高清| 久久影院123| 一本一本综合久久| 国产黄色视频一区二区在线观看| 久久99热6这里只有精品| 国产亚洲av片在线观看秒播厂| 97超碰精品成人国产| 亚洲成人av在线免费| 美女大奶头黄色视频| 亚洲美女黄色视频免费看| 桃花免费在线播放| 午夜激情久久久久久久| 校园人妻丝袜中文字幕| 一二三四中文在线观看免费高清| 王馨瑶露胸无遮挡在线观看| 欧美精品亚洲一区二区| 国产 精品1| 国产亚洲一区二区精品| 久久久国产欧美日韩av| 亚洲国产精品999| 国产熟女欧美一区二区| 欧美日韩在线观看h| 大码成人一级视频| 91在线精品国自产拍蜜月| 在现免费观看毛片| 国产91av在线免费观看| 秋霞伦理黄片| 欧美精品一区二区免费开放| 免费久久久久久久精品成人欧美视频 | 少妇人妻精品综合一区二区| 久久亚洲国产成人精品v| 国产黄频视频在线观看| 18禁在线播放成人免费| 18禁动态无遮挡网站| 波野结衣二区三区在线| 十分钟在线观看高清视频www | av天堂中文字幕网| 亚州av有码| 熟女人妻精品中文字幕| 日日摸夜夜添夜夜添av毛片| 免费观看a级毛片全部| 久久韩国三级中文字幕| 老女人水多毛片| 成人毛片a级毛片在线播放| 91久久精品国产一区二区三区| 国产精品国产三级国产专区5o| 十八禁网站网址无遮挡 | 精品一品国产午夜福利视频| 亚洲精品日韩在线中文字幕| 久久人人爽av亚洲精品天堂| 中文字幕免费在线视频6| 成人亚洲欧美一区二区av| 中文精品一卡2卡3卡4更新| 免费黄色在线免费观看| 99re6热这里在线精品视频| av福利片在线| 国产精品久久久久久av不卡| 免费少妇av软件| 男女边摸边吃奶| 黄色怎么调成土黄色| 乱人伦中国视频| 亚洲国产精品一区二区三区在线| 91久久精品国产一区二区三区| 国产精品久久久久久精品古装| 性高湖久久久久久久久免费观看| 最近最新中文字幕免费大全7| 精品一区二区三区视频在线| 午夜福利,免费看| 国内揄拍国产精品人妻在线| 日韩亚洲欧美综合| 欧美老熟妇乱子伦牲交| 极品人妻少妇av视频| 亚洲av免费高清在线观看| av视频免费观看在线观看| 深夜a级毛片| 一区二区三区免费毛片| 日韩电影二区| 如日韩欧美国产精品一区二区三区 | 有码 亚洲区| 国产精品一二三区在线看| 久久久久国产网址| 国产成人免费观看mmmm| 国产成人aa在线观看| 一个人看视频在线观看www免费| 亚洲美女视频黄频| 在线 av 中文字幕| 中国美白少妇内射xxxbb| 3wmmmm亚洲av在线观看| 国产在视频线精品| av专区在线播放| 深夜a级毛片| 人妻人人澡人人爽人人| 国产av一区二区精品久久| 日本色播在线视频| 国产在视频线精品| 国产精品不卡视频一区二区| 亚洲欧洲国产日韩| 成年人午夜在线观看视频| 久久亚洲国产成人精品v| 免费播放大片免费观看视频在线观看| 国产精品嫩草影院av在线观看| 亚洲精华国产精华液的使用体验| 久久久久精品久久久久真实原创| 国产精品嫩草影院av在线观看| 熟女av电影| 精品国产一区二区久久| 亚洲国产色片| 一个人看视频在线观看www免费| 日韩制服骚丝袜av| a 毛片基地| 免费看光身美女| 女人精品久久久久毛片| 国产中年淑女户外野战色| 黄色配什么色好看| 少妇丰满av| 午夜免费鲁丝| 国产色爽女视频免费观看| 大香蕉久久网| 国产精品秋霞免费鲁丝片| 国产精品久久久久久精品电影小说| 国产成人一区二区在线| 国产成人免费观看mmmm| 国产日韩欧美视频二区| 成人亚洲精品一区在线观看| 一级a做视频免费观看| 男人狂女人下面高潮的视频| 肉色欧美久久久久久久蜜桃| 成人漫画全彩无遮挡| 国产乱人偷精品视频| 久久久亚洲精品成人影院| 黄色毛片三级朝国网站 | 天堂8中文在线网| 97在线人人人人妻| 午夜日本视频在线| 亚洲高清免费不卡视频| 80岁老熟妇乱子伦牲交| 日韩精品有码人妻一区| 乱系列少妇在线播放| av有码第一页| 精品久久国产蜜桃| 成年av动漫网址| 91午夜精品亚洲一区二区三区| 美女福利国产在线| a级一级毛片免费在线观看| 最新中文字幕久久久久| 国产精品久久久久久精品电影小说| 不卡视频在线观看欧美| 777米奇影视久久| 性色av一级| 亚洲电影在线观看av| 偷拍熟女少妇极品色| 精品亚洲成国产av| 国产精品偷伦视频观看了| 日韩免费高清中文字幕av| 人妻一区二区av| 精品一区二区免费观看| 三级国产精品片| 夜夜看夜夜爽夜夜摸| 国产综合精华液| 欧美一级a爱片免费观看看| 成人黄色视频免费在线看| 老司机影院成人| 亚洲国产日韩一区二区| 日韩av免费高清视频| 免费黄网站久久成人精品| 亚洲激情五月婷婷啪啪| 午夜老司机福利剧场| 欧美老熟妇乱子伦牲交| 各种免费的搞黄视频| 嘟嘟电影网在线观看| a 毛片基地| 成年人免费黄色播放视频 | 国产色爽女视频免费观看| 日本av免费视频播放| 色婷婷久久久亚洲欧美| 大香蕉久久网| 国国产精品蜜臀av免费| 高清在线视频一区二区三区| 免费人妻精品一区二区三区视频| 免费观看的影片在线观看| 在线精品无人区一区二区三| 啦啦啦在线观看免费高清www| 亚洲不卡免费看| 午夜福利网站1000一区二区三区| 欧美国产精品一级二级三级 | 中文精品一卡2卡3卡4更新| 少妇裸体淫交视频免费看高清| 国产日韩一区二区三区精品不卡 | 国产精品久久久久久av不卡| av黄色大香蕉| 熟女电影av网| 亚洲精品乱久久久久久| 嫩草影院入口| 亚洲美女搞黄在线观看| .国产精品久久| 99久久中文字幕三级久久日本| 热99国产精品久久久久久7| 中文字幕久久专区| 亚洲av中文av极速乱| 人人妻人人添人人爽欧美一区卜| 久久久久久人妻| 欧美日韩一区二区视频在线观看视频在线| 2018国产大陆天天弄谢| 国产日韩欧美视频二区| 亚洲国产日韩一区二区| www.色视频.com| 久久久久网色| 亚洲国产精品999| 亚洲欧美成人综合另类久久久| 亚洲综合精品二区| 亚洲人成网站在线观看播放| 亚洲美女视频黄频| 国产精品久久久久久精品古装| 22中文网久久字幕| 亚洲av免费高清在线观看| 成人亚洲欧美一区二区av| av国产精品久久久久影院| 欧美老熟妇乱子伦牲交| 夫妻午夜视频| 国产欧美日韩精品一区二区| 丰满人妻一区二区三区视频av| 久久久久久久大尺度免费视频| 汤姆久久久久久久影院中文字幕| 欧美日韩av久久| 最近中文字幕2019免费版| 国精品久久久久久国模美| 日日摸夜夜添夜夜添av毛片| 97在线视频观看| 丝袜喷水一区| 菩萨蛮人人尽说江南好唐韦庄| 欧美 日韩 精品 国产| 亚洲真实伦在线观看| 国产成人aa在线观看| 寂寞人妻少妇视频99o| 能在线免费看毛片的网站| www.av在线官网国产| 久久综合国产亚洲精品| 蜜桃在线观看..| 九色成人免费人妻av| 国产色爽女视频免费观看| 99久久精品热视频| 在线亚洲精品国产二区图片欧美 | 噜噜噜噜噜久久久久久91| 男人添女人高潮全过程视频| 人人妻人人澡人人看| 三级国产精品片| 美女大奶头黄色视频| 国产精品久久久久成人av| 亚洲经典国产精华液单| 一区二区三区四区激情视频| 日日摸夜夜添夜夜添av毛片| av一本久久久久| 亚洲三级黄色毛片| 日韩中文字幕视频在线看片| 国产乱人偷精品视频| 99热国产这里只有精品6| 亚洲成人一二三区av| 亚洲精品自拍成人| 黄色视频在线播放观看不卡| 国产精品一区二区性色av| 精品久久久精品久久久| 久久女婷五月综合色啪小说| 国产成人a∨麻豆精品| 欧美日本中文国产一区发布| 日韩成人av中文字幕在线观看| 亚洲精华国产精华液的使用体验| 国产男女内射视频| 男人添女人高潮全过程视频| 午夜激情久久久久久久| 99热这里只有是精品50| 一级毛片 在线播放| 亚洲不卡免费看| 少妇被粗大猛烈的视频| 特大巨黑吊av在线直播| 亚洲欧美精品专区久久| 麻豆乱淫一区二区| 2022亚洲国产成人精品| 99热6这里只有精品| 男女边摸边吃奶| 乱系列少妇在线播放| 亚洲精品aⅴ在线观看| 夜夜骑夜夜射夜夜干| 久久久久久久久久人人人人人人| 美女内射精品一级片tv| 久久女婷五月综合色啪小说| 嫩草影院入口| av天堂中文字幕网| av国产久精品久网站免费入址| 国产精品.久久久| 日本爱情动作片www.在线观看| 人人妻人人爽人人添夜夜欢视频 | 精品久久久噜噜| 美女国产视频在线观看| 精品亚洲乱码少妇综合久久| 99久久精品热视频| 久久精品国产鲁丝片午夜精品| 又大又黄又爽视频免费| 蜜桃久久精品国产亚洲av| 国产一区二区在线观看日韩| 中文字幕免费在线视频6| 欧美bdsm另类| 亚洲美女搞黄在线观看| 极品人妻少妇av视频| 亚洲精品一区蜜桃| 国产精品成人在线| 亚洲av免费高清在线观看| 国产毛片在线视频| 国产伦精品一区二区三区视频9| 春色校园在线视频观看| 日本黄色日本黄色录像| 我要看日韩黄色一级片| 日韩欧美一区视频在线观看 | 久久99热6这里只有精品| 成人亚洲精品一区在线观看| 久久ye,这里只有精品| 丝袜喷水一区| 日韩人妻高清精品专区| 啦啦啦视频在线资源免费观看| 亚洲国产精品999| 国产成人精品婷婷| 99热这里只有精品一区| 国产免费一级a男人的天堂| 卡戴珊不雅视频在线播放| a级片在线免费高清观看视频| 自拍欧美九色日韩亚洲蝌蚪91 | 这个男人来自地球电影免费观看 | 男女啪啪激烈高潮av片| 毛片一级片免费看久久久久| 男人添女人高潮全过程视频| 一区在线观看完整版| 精品人妻熟女毛片av久久网站| 中文精品一卡2卡3卡4更新| 内地一区二区视频在线| 欧美亚洲 丝袜 人妻 在线| 一本大道久久a久久精品| 岛国毛片在线播放| 我要看日韩黄色一级片| 国产免费又黄又爽又色| 草草在线视频免费看| 久久精品久久久久久久性| 中文天堂在线官网| 亚洲va在线va天堂va国产| 大香蕉久久网| 一级二级三级毛片免费看| 中文字幕精品免费在线观看视频 | 成人综合一区亚洲| 久久影院123| 亚洲欧美成人精品一区二区| 午夜福利视频精品| 人妻夜夜爽99麻豆av| 免费在线观看成人毛片| 国内揄拍国产精品人妻在线| 久久国产精品男人的天堂亚洲 | 2018国产大陆天天弄谢| 国产伦在线观看视频一区| 一二三四中文在线观看免费高清| 天堂中文最新版在线下载| 精品人妻熟女av久视频| 69精品国产乱码久久久| 国产高清三级在线| 久久影院123| 伦理电影大哥的女人| 99久久精品国产国产毛片| av网站免费在线观看视频| 另类精品久久| 亚洲欧美成人综合另类久久久| 在线观看一区二区三区激情| 成人漫画全彩无遮挡| 久久久国产精品麻豆| 老司机影院成人| 在线观看免费日韩欧美大片 | 亚洲成色77777| 在线 av 中文字幕| 欧美日韩视频精品一区| 久久久国产欧美日韩av| 国产精品一区二区在线不卡| 成人免费观看视频高清| 日韩亚洲欧美综合| 国产 一区精品| 乱码一卡2卡4卡精品| 久久久精品94久久精品| 国产精品.久久久| 在线 av 中文字幕| 久久国产乱子免费精品| 亚洲av成人精品一区久久| 日韩强制内射视频| 国产精品一区二区三区四区免费观看| 在线 av 中文字幕| 国产熟女午夜一区二区三区 | 在线精品无人区一区二区三| 黄色怎么调成土黄色| 亚洲欧美成人精品一区二区| av天堂久久9| 老女人水多毛片| 精品久久久久久久久亚洲| 一本一本综合久久| 三上悠亚av全集在线观看 | 成人亚洲精品一区在线观看| 卡戴珊不雅视频在线播放| 人妻一区二区av| 9色porny在线观看| 少妇被粗大的猛进出69影院 | 高清午夜精品一区二区三区| 观看美女的网站| 又爽又黄a免费视频| 国产一区二区在线观看av| 香蕉精品网在线| 国产极品天堂在线| 免费黄色在线免费观看| 久久99一区二区三区| 最近手机中文字幕大全| 国产在线免费精品| 免费观看性生交大片5| 国产又色又爽无遮挡免| 国产精品欧美亚洲77777| 亚洲国产毛片av蜜桃av| 中国三级夫妇交换| 黑人巨大精品欧美一区二区蜜桃 | 国产精品久久久久久精品电影小说| 亚洲第一区二区三区不卡| 亚洲成人手机| 久久久久久久久久人人人人人人| 国产综合精华液| 亚洲电影在线观看av| 在现免费观看毛片| 日韩中字成人| 欧美97在线视频| 午夜久久久在线观看| 成年人免费黄色播放视频 | 狂野欧美白嫩少妇大欣赏| 少妇被粗大猛烈的视频| 97精品久久久久久久久久精品| 精品国产一区二区久久| tube8黄色片| 免费大片18禁| 国产欧美亚洲国产| 乱系列少妇在线播放| 成人国产麻豆网| 亚洲国产精品国产精品| 午夜久久久在线观看| 在线播放无遮挡| 久久综合国产亚洲精品| 大话2 男鬼变身卡| 狂野欧美激情性bbbbbb| 国内揄拍国产精品人妻在线| 午夜视频国产福利| 人人澡人人妻人| 国产成人精品久久久久久| 简卡轻食公司| 精品国产一区二区久久| 丝袜喷水一区| 久久青草综合色| 国产黄片美女视频| 麻豆乱淫一区二区| 日本色播在线视频| 99久久精品一区二区三区| 一级毛片电影观看| 欧美最新免费一区二区三区| 国产在视频线精品| av免费观看日本| 国产欧美日韩精品一区二区| 在线观看国产h片| 久久人妻熟女aⅴ| 国产精品一区www在线观看| 免费播放大片免费观看视频在线观看| 91成人精品电影| 校园人妻丝袜中文字幕| 男的添女的下面高潮视频| 99热这里只有是精品在线观看| 免费不卡的大黄色大毛片视频在线观看| 内射极品少妇av片p| 大片免费播放器 马上看| 免费高清在线观看视频在线观看| 99热这里只有是精品50| 九色成人免费人妻av| 色5月婷婷丁香| 成人漫画全彩无遮挡| 中文字幕av电影在线播放| 人人妻人人爽人人添夜夜欢视频 | 日韩中字成人| a级毛色黄片| 亚洲精品视频女| 丰满饥渴人妻一区二区三| 国产精品成人在线| 久久久久久久国产电影| 最黄视频免费看| av黄色大香蕉| 亚洲精品日韩av片在线观看| 欧美日韩国产mv在线观看视频| 国产精品.久久久| av卡一久久| 欧美xxⅹ黑人| 99久久精品国产国产毛片| 国内少妇人妻偷人精品xxx网站| 日韩亚洲欧美综合| 22中文网久久字幕| av女优亚洲男人天堂| 午夜福利网站1000一区二区三区| 少妇 在线观看| 久久久亚洲精品成人影院| 啦啦啦啦在线视频资源| 永久免费av网站大全| 哪个播放器可以免费观看大片| 国产片特级美女逼逼视频| 天堂中文最新版在线下载| 性色avwww在线观看| 人妻制服诱惑在线中文字幕| 能在线免费看毛片的网站| 高清午夜精品一区二区三区| av卡一久久| 成人毛片a级毛片在线播放| 国产精品国产三级国产专区5o| 九九久久精品国产亚洲av麻豆|