• <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ū)土壤水分反演研究——以北京市為例
    国产亚洲最大av| 80岁老熟妇乱子伦牲交| 看免费成人av毛片| 超色免费av| 啦啦啦在线观看免费高清www| 日本爱情动作片www.在线观看| 国产 精品1| 三上悠亚av全集在线观看| 两个人免费观看高清视频| 精品国产露脸久久av麻豆| 国产熟女午夜一区二区三区| 亚洲欧洲日产国产| 亚洲av免费高清在线观看| 18禁观看日本| 黄色一级大片看看| 国产一区二区三区av在线| 精品一区在线观看国产| 少妇熟女欧美另类| 尾随美女入室| 两性夫妻黄色片| 黑人猛操日本美女一级片| 99久国产av精品国产电影| 久久精品国产a三级三级三级| 国产av一区二区精品久久| 国产精品国产av在线观看| 精品一品国产午夜福利视频| 男女午夜视频在线观看| tube8黄色片| 韩国精品一区二区三区| 天天躁夜夜躁狠狠躁躁| 伊人久久大香线蕉亚洲五| 日日啪夜夜爽| 亚洲三级黄色毛片| 精品少妇黑人巨大在线播放| 国产精品一区二区在线观看99| 色网站视频免费| 天堂中文最新版在线下载| 五月开心婷婷网| 黄色一级大片看看| 夫妻性生交免费视频一级片| 美女脱内裤让男人舔精品视频| 久久久久久久久久久久大奶| 亚洲av国产av综合av卡| 视频在线观看一区二区三区| 观看av在线不卡| 日韩精品有码人妻一区| 久久久久国产一级毛片高清牌| 国产精品国产三级专区第一集| 国产精品国产三级专区第一集| 亚洲国产av新网站| 激情五月婷婷亚洲| 欧美人与性动交α欧美软件| 如何舔出高潮| 日韩精品免费视频一区二区三区| 18禁国产床啪视频网站| 精品少妇黑人巨大在线播放| 日本欧美视频一区| 欧美激情 高清一区二区三区| 老汉色av国产亚洲站长工具| 欧美日韩一区二区视频在线观看视频在线| 午夜福利网站1000一区二区三区| 日韩一本色道免费dvd| 男的添女的下面高潮视频| 欧美精品一区二区免费开放| 亚洲欧美成人综合另类久久久| 一本—道久久a久久精品蜜桃钙片| 久久久精品94久久精品| 高清不卡的av网站| 亚洲国产色片| 99久久人妻综合| 国产免费又黄又爽又色| 高清视频免费观看一区二区| 欧美激情极品国产一区二区三区| 亚洲欧美一区二区三区黑人 | 亚洲色图 男人天堂 中文字幕| 女人精品久久久久毛片| 久久综合国产亚洲精品| 久久久久网色| 成年av动漫网址| 日韩一区二区视频免费看| 精品酒店卫生间| 欧美精品一区二区大全| 久久久久精品久久久久真实原创| 中文字幕色久视频| 午夜免费鲁丝| 日本av手机在线免费观看| 久久午夜福利片| 午夜老司机福利剧场| 国产在视频线精品| 久久人人爽人人片av| 纯流量卡能插随身wifi吗| 一个人免费看片子| 国产精品成人在线| 精品久久久久久电影网| 欧美精品国产亚洲| 天堂8中文在线网| 久久国产亚洲av麻豆专区| 欧美日韩精品网址| 亚洲av在线观看美女高潮| 国产色婷婷99| 亚洲国产日韩一区二区| 69精品国产乱码久久久| 国语对白做爰xxxⅹ性视频网站| 国产亚洲精品第一综合不卡| 国产精品 国内视频| 国产片内射在线| 99热全是精品| 国产成人精品在线电影| 国产精品三级大全| 国产精品久久久久久精品古装| 欧美日本中文国产一区发布| 日韩大片免费观看网站| 老司机影院成人| 久久婷婷青草| 99热国产这里只有精品6| 精品卡一卡二卡四卡免费| av在线app专区| 欧美激情高清一区二区三区 | 丰满迷人的少妇在线观看| 热re99久久精品国产66热6| 久久鲁丝午夜福利片| 日韩精品有码人妻一区| 超碰成人久久| 亚洲男人天堂网一区| 一区在线观看完整版| 亚洲成人手机| 赤兔流量卡办理| 中文欧美无线码| 天美传媒精品一区二区| 国产精品久久久av美女十八| 一区二区三区精品91| 成人18禁高潮啪啪吃奶动态图| 婷婷色av中文字幕| 成年美女黄网站色视频大全免费| 黄片播放在线免费| 老女人水多毛片| 男女下面插进去视频免费观看| 国产亚洲最大av| 欧美日韩亚洲高清精品| 亚洲av日韩在线播放| 男女高潮啪啪啪动态图| 亚洲国产欧美日韩在线播放| 国产福利在线免费观看视频| 亚洲国产精品一区三区| 久久久久久久国产电影| 日本-黄色视频高清免费观看| 夫妻午夜视频| 国产深夜福利视频在线观看| 视频在线观看一区二区三区| 亚洲欧美一区二区三区久久| 赤兔流量卡办理| 国产熟女午夜一区二区三区| 我要看黄色一级片免费的| h视频一区二区三区| 一级,二级,三级黄色视频| 天美传媒精品一区二区| 美女中出高潮动态图| 亚洲精品久久久久久婷婷小说| 久久99热这里只频精品6学生| 久久精品国产鲁丝片午夜精品| 制服丝袜香蕉在线| 中文字幕av电影在线播放| 天天躁狠狠躁夜夜躁狠狠躁| 久久久久久久久久久久大奶| 久久午夜福利片| 久久狼人影院| 亚洲av成人精品一二三区| 高清欧美精品videossex| 两个人免费观看高清视频| 欧美精品一区二区免费开放| 2021少妇久久久久久久久久久| 日韩一区二区视频免费看| 国产一区二区三区av在线| 好男人视频免费观看在线| 久久精品久久久久久久性| 亚洲男人天堂网一区| 日本vs欧美在线观看视频| 亚洲欧美一区二区三区国产| 老司机影院毛片| 免费久久久久久久精品成人欧美视频| 在线亚洲精品国产二区图片欧美| 成人毛片a级毛片在线播放| 天堂8中文在线网| 91成人精品电影| 2022亚洲国产成人精品| 春色校园在线视频观看| 日本黄色日本黄色录像| 亚洲成国产人片在线观看| 一级黄片播放器| 女人久久www免费人成看片| xxx大片免费视频| 国产精品一二三区在线看| 人人妻人人澡人人爽人人夜夜| √禁漫天堂资源中文www| 久久久久久人人人人人| 亚洲精品日韩在线中文字幕| 欧美人与性动交α欧美软件| 99国产精品免费福利视频| 桃花免费在线播放| 飞空精品影院首页| 一级,二级,三级黄色视频| 免费观看无遮挡的男女| 熟女av电影| 亚洲伊人久久精品综合| 国产成人欧美| 久久免费观看电影| av卡一久久| 日韩一区二区视频免费看| 少妇被粗大猛烈的视频| 秋霞在线观看毛片| 99国产综合亚洲精品| 成人漫画全彩无遮挡| 极品少妇高潮喷水抽搐| 哪个播放器可以免费观看大片| 夜夜骑夜夜射夜夜干| 宅男免费午夜| 欧美日韩亚洲国产一区二区在线观看 | 十八禁网站网址无遮挡| 国产精品 欧美亚洲| 黑人巨大精品欧美一区二区蜜桃| 中文欧美无线码| 免费在线观看视频国产中文字幕亚洲 | 国产成人免费观看mmmm| 激情视频va一区二区三区| 在线天堂最新版资源| 国产日韩一区二区三区精品不卡| 精品一区二区三卡| 精品久久久久久电影网| 欧美日韩亚洲高清精品| 国产精品偷伦视频观看了| 一级毛片我不卡| 国产成人一区二区在线| 男女边吃奶边做爰视频| 国产精品久久久av美女十八| 极品人妻少妇av视频| 搡老乐熟女国产| 极品少妇高潮喷水抽搐| 精品国产乱码久久久久久小说| 黄色视频在线播放观看不卡| 欧美激情高清一区二区三区 | 欧美 日韩 精品 国产| 制服丝袜香蕉在线| 久久韩国三级中文字幕| 十八禁网站网址无遮挡| 亚洲欧洲日产国产| 巨乳人妻的诱惑在线观看| 久久精品国产自在天天线| 日韩 亚洲 欧美在线| 成人亚洲欧美一区二区av| 中国国产av一级| 在线观看一区二区三区激情| 看非洲黑人一级黄片| 国产精品99久久99久久久不卡 | 99热网站在线观看| 女的被弄到高潮叫床怎么办| 91aial.com中文字幕在线观看| 成年人午夜在线观看视频| 精品酒店卫生间| 国产人伦9x9x在线观看 | 久久影院123| 女人高潮潮喷娇喘18禁视频| 黄片无遮挡物在线观看| 麻豆乱淫一区二区| 男女边摸边吃奶| 最近手机中文字幕大全| 成人亚洲精品一区在线观看| 亚洲 欧美一区二区三区| 久久久精品免费免费高清| 久久久久久伊人网av| 欧美亚洲日本最大视频资源| 只有这里有精品99| 亚洲精品成人av观看孕妇| 欧美xxⅹ黑人| 性色avwww在线观看| 久久久久国产网址| 搡老乐熟女国产| 精品一区二区三卡| 日韩中文字幕视频在线看片| 18禁观看日本| 少妇被粗大猛烈的视频| 久热久热在线精品观看| 大香蕉久久成人网| 国产精品国产三级国产专区5o| 免费观看av网站的网址| 午夜福利视频精品| 九草在线视频观看| 久久精品久久精品一区二区三区| 久久亚洲国产成人精品v| 2018国产大陆天天弄谢| 男人添女人高潮全过程视频| 亚洲国产欧美在线一区| 国产一区有黄有色的免费视频| 免费播放大片免费观看视频在线观看| 美女视频免费永久观看网站| 在线观看免费日韩欧美大片| 亚洲精品第二区| 成年美女黄网站色视频大全免费| 欧美日韩国产mv在线观看视频| 亚洲国产日韩一区二区| 久久青草综合色| 国产又色又爽无遮挡免| 欧美国产精品一级二级三级| 成年动漫av网址| 18禁裸乳无遮挡动漫免费视频| 成年av动漫网址| 国产 一区精品| 国产 精品1| 天堂中文最新版在线下载| 国产免费一区二区三区四区乱码| 国产精品一区二区在线观看99| 看免费av毛片| 亚洲国产日韩一区二区| 人人澡人人妻人| 亚洲人成网站在线观看播放| 91久久精品国产一区二区三区| 日本欧美视频一区| 十分钟在线观看高清视频www| 18禁观看日本| 人人妻人人澡人人看| 精品国产超薄肉色丝袜足j| 男女边摸边吃奶| 亚洲国产色片| 久久久久久久久久久久大奶| 亚洲欧美一区二区三区国产| 久久免费观看电影| 亚洲精品成人av观看孕妇| 天天操日日干夜夜撸| 婷婷色综合www| 久久女婷五月综合色啪小说| 99九九在线精品视频| 国产精品.久久久| 看免费av毛片| 黄色一级大片看看| 国产日韩一区二区三区精品不卡| 日本91视频免费播放| 美女主播在线视频| 啦啦啦啦在线视频资源| 大话2 男鬼变身卡| 成人毛片60女人毛片免费| 王馨瑶露胸无遮挡在线观看| 亚洲av电影在线进入| av在线播放精品| 国产探花极品一区二区| 日本爱情动作片www.在线观看| 人妻一区二区av| 国产精品一区二区在线不卡| 国产伦理片在线播放av一区| 观看美女的网站| 国产精品不卡视频一区二区| 9色porny在线观看| 不卡视频在线观看欧美| 国产av国产精品国产| 午夜激情久久久久久久| 免费黄色在线免费观看| 母亲3免费完整高清在线观看 | 看十八女毛片水多多多| 亚洲国产欧美网| 成人免费观看视频高清| 王馨瑶露胸无遮挡在线观看| 国产熟女欧美一区二区| 亚洲av免费高清在线观看| 中文字幕人妻丝袜制服| 国产片特级美女逼逼视频| 最近最新中文字幕大全免费视频 | 毛片一级片免费看久久久久| 黄色一级大片看看| 精品亚洲成国产av| 精品午夜福利在线看| 午夜影院在线不卡| av网站在线播放免费| 亚洲三级黄色毛片| 一级毛片电影观看| 尾随美女入室| 久久97久久精品| 欧美激情极品国产一区二区三区| 啦啦啦在线免费观看视频4| 91精品三级在线观看| 9色porny在线观看| 欧美激情极品国产一区二区三区| 午夜福利视频精品| 国产熟女欧美一区二区| 成年人午夜在线观看视频| 国产成人精品久久二区二区91 | 精品国产超薄肉色丝袜足j| 国产黄频视频在线观看| 丝袜在线中文字幕| 国产精品秋霞免费鲁丝片| 熟妇人妻不卡中文字幕| 国产精品秋霞免费鲁丝片| 亚洲国产精品国产精品| 可以免费在线观看a视频的电影网站 | 两个人免费观看高清视频| 男女高潮啪啪啪动态图| 深夜精品福利| 亚洲,一卡二卡三卡| 日韩,欧美,国产一区二区三区| 校园人妻丝袜中文字幕| 90打野战视频偷拍视频| 亚洲欧美清纯卡通| 满18在线观看网站| 建设人人有责人人尽责人人享有的| 汤姆久久久久久久影院中文字幕| 亚洲成人一二三区av| 日本色播在线视频| 韩国精品一区二区三区| 可以免费在线观看a视频的电影网站 | 亚洲第一青青草原| 最近中文字幕高清免费大全6| 国产一区二区激情短视频 | 满18在线观看网站| 狂野欧美激情性bbbbbb| 亚洲精品,欧美精品| 国产一级毛片在线| 性少妇av在线| 久久国产精品男人的天堂亚洲| 日本猛色少妇xxxxx猛交久久| 日韩视频在线欧美| 美女国产视频在线观看| 久久久久久免费高清国产稀缺| 久久久久久久大尺度免费视频| 亚洲成人av在线免费| 久久97久久精品| 国产一区二区三区综合在线观看| 亚洲国产精品一区三区| 成人午夜精彩视频在线观看| 久久99精品国语久久久| 亚洲精品第二区| 1024视频免费在线观看| 曰老女人黄片| 国产成人精品久久久久久| 一级毛片我不卡| 亚洲av电影在线进入| 一级毛片电影观看| 母亲3免费完整高清在线观看 | 色视频在线一区二区三区| 日韩电影二区| 桃花免费在线播放| 久久久久国产精品人妻一区二区| 女的被弄到高潮叫床怎么办| 交换朋友夫妻互换小说| 欧美日韩av久久| 满18在线观看网站| 欧美精品国产亚洲| 国产野战对白在线观看| 一本—道久久a久久精品蜜桃钙片| 国产乱人偷精品视频| 国产精品一国产av| 一边摸一边做爽爽视频免费| 九草在线视频观看| 久久影院123| 赤兔流量卡办理| 激情五月婷婷亚洲| 成年动漫av网址| 国产精品国产三级专区第一集| 人人澡人人妻人| 三上悠亚av全集在线观看| 91aial.com中文字幕在线观看| 亚洲欧美一区二区三区国产| 国产免费福利视频在线观看| 欧美日韩国产mv在线观看视频| 国产成人精品福利久久| 日本爱情动作片www.在线观看| 亚洲精品一二三| 亚洲欧洲日产国产| 欧美精品高潮呻吟av久久| 日韩伦理黄色片| 999精品在线视频| 国产福利在线免费观看视频| 春色校园在线视频观看| 国产有黄有色有爽视频| 亚洲国产欧美网| 久久久久久久久久久免费av| 少妇猛男粗大的猛烈进出视频| 十八禁网站网址无遮挡| 亚洲精品av麻豆狂野| 国产白丝娇喘喷水9色精品| 国产精品国产三级国产专区5o| 美女主播在线视频| 一级爰片在线观看| av天堂久久9| 国产精品.久久久| 欧美国产精品va在线观看不卡| 日韩av不卡免费在线播放| 如何舔出高潮| 国产一区亚洲一区在线观看| 欧美黄色片欧美黄色片| 咕卡用的链子| 亚洲精品日韩在线中文字幕| 免费黄网站久久成人精品| 国产精品久久久久成人av| 亚洲av欧美aⅴ国产| 免费看av在线观看网站| 一级毛片 在线播放| 赤兔流量卡办理| 高清不卡的av网站| 永久免费av网站大全| 亚洲成人一二三区av| 桃花免费在线播放| 1024香蕉在线观看| videossex国产| 电影成人av| 久久久久久人人人人人| 亚洲综合色网址| 亚洲,欧美精品.| 嫩草影院入口| 中文天堂在线官网| 中文字幕制服av| 中文字幕色久视频| 高清视频免费观看一区二区| 亚洲精品国产av成人精品| 这个男人来自地球电影免费观看 | 国产免费又黄又爽又色| 美女xxoo啪啪120秒动态图| 免费高清在线观看日韩| 狂野欧美激情性bbbbbb| 日本wwww免费看| 亚洲精品久久成人aⅴ小说| 久久精品久久精品一区二区三区| 在线天堂最新版资源| 美女国产高潮福利片在线看| 美女大奶头黄色视频| 九色亚洲精品在线播放| 亚洲精品视频女| 精品少妇一区二区三区视频日本电影 | 80岁老熟妇乱子伦牲交| 久久97久久精品| 国产福利在线免费观看视频| 欧美亚洲日本最大视频资源| 日韩三级伦理在线观看| 久久精品国产亚洲av涩爱| 99九九在线精品视频| 日韩人妻精品一区2区三区| 国产一级毛片在线| 午夜激情久久久久久久| 国产一区二区三区综合在线观看| 亚洲精品日本国产第一区| 街头女战士在线观看网站| 国产日韩欧美亚洲二区| 搡老乐熟女国产| 亚洲av电影在线进入| 在线免费观看不下载黄p国产| 日本wwww免费看| 热99久久久久精品小说推荐| 嫩草影院入口| 最新的欧美精品一区二区| 伊人久久国产一区二区| 伦理电影免费视频| 黑丝袜美女国产一区| 国产 精品1| 亚洲av成人精品一二三区| 午夜福利一区二区在线看| 在线观看免费高清a一片| 91午夜精品亚洲一区二区三区| 久久人人97超碰香蕉20202| 少妇的丰满在线观看| 在线精品无人区一区二区三| 久久久久网色| 视频区图区小说| 热99久久久久精品小说推荐| 人妻人人澡人人爽人人| 香蕉国产在线看| 精品国产一区二区三区四区第35| 中文字幕av电影在线播放| 日韩人妻精品一区2区三区| 国产深夜福利视频在线观看| 在线观看国产h片| 成人免费观看视频高清| 一级黄片播放器| 最近最新中文字幕免费大全7| 亚洲成av片中文字幕在线观看 | 五月开心婷婷网| 久久久久国产精品人妻一区二区| 国产一区二区在线观看av| 欧美激情 高清一区二区三区| 午夜福利影视在线免费观看| 亚洲av国产av综合av卡| av在线app专区| 美女xxoo啪啪120秒动态图| 久久精品国产亚洲av高清一级| 在线亚洲精品国产二区图片欧美| 欧美国产精品一级二级三级| 最近中文字幕2019免费版| 国产亚洲一区二区精品| 国产精品亚洲av一区麻豆 | 欧美亚洲 丝袜 人妻 在线| 极品少妇高潮喷水抽搐| 人人澡人人妻人| 国产有黄有色有爽视频| 久热久热在线精品观看| 欧美国产精品一级二级三级| 国产男女超爽视频在线观看| 麻豆乱淫一区二区| xxx大片免费视频| 国产白丝娇喘喷水9色精品| 国产日韩欧美视频二区| av免费在线看不卡| 成年人免费黄色播放视频| 下体分泌物呈黄色| 久久这里只有精品19| 18+在线观看网站| 久久精品国产综合久久久| 国产精品偷伦视频观看了| 男女无遮挡免费网站观看| 久久午夜福利片| 日本av免费视频播放| 久久婷婷青草| 人妻人人澡人人爽人人| 男人操女人黄网站| 午夜免费观看性视频| 好男人视频免费观看在线| 亚洲色图 男人天堂 中文字幕| 日韩欧美一区视频在线观看| 亚洲欧美一区二区三区久久| 哪个播放器可以免费观看大片| 赤兔流量卡办理|