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

    不同季節(jié)強堿土土壤呼吸影響因子分析與模型預(yù)測

    2017-03-02 08:43:19張錄張芳熊黑鋼段鵬程李榮榮張南
    關(guān)鍵詞:強堿通徑土壤濕度

    張錄,張芳,2,熊黑鋼,3,段鵬程,李榮榮,張南

    (1.新疆大學(xué)資源與環(huán)境科學(xué)學(xué)院,教育部綠洲生態(tài)重點試驗室,新疆烏魯木齊830046;2.新疆大學(xué)生態(tài)學(xué)博士后流動站,新疆烏魯木齊830046;3.北京聯(lián)合大學(xué)應(yīng)用文理學(xué)院,北京100191)

    不同季節(jié)強堿土土壤呼吸影響因子分析與模型預(yù)測

    張錄1,張芳1,2,熊黑鋼1,3,段鵬程1,李榮榮1,張南1

    (1.新疆大學(xué)資源與環(huán)境科學(xué)學(xué)院,教育部綠洲生態(tài)重點試驗室,新疆烏魯木齊830046;2.新疆大學(xué)生態(tài)學(xué)博士后流動站,新疆烏魯木齊830046;3.北京聯(lián)合大學(xué)應(yīng)用文理學(xué)院,北京100191)

    利用LI-8100土壤碳通量測量儀測定了春夏秋三季晴朗天氣下強堿土土壤呼吸速率、溫度(氣溫和地溫)、濕度(空氣相對濕度和土壤濕度)數(shù)據(jù),分析了它們之間的相關(guān)關(guān)系,獲得不同季節(jié)對土壤呼吸影響較大的因子,并建立不同類別的多種回歸模型;在精度檢驗及簡單易行原則基礎(chǔ)上,得到各季節(jié)土壤呼吸預(yù)測的最優(yōu)模型。結(jié)果表明:(1)雖然溫濕度均是影響不同季節(jié)強堿土土壤呼吸的主要因素,但均以溫度的影響較大,其中氣溫是春秋兩季土壤呼吸的最大直接影響因子,地溫是夏季的最大直接影響因子,而土壤濕度為各季節(jié)最大的間接影響因子。(2)春秋季土壤呼吸的最佳預(yù)測模型均為10 cm處氣溫和土壤濕度所建的雙因子方程,該方程具有較小的均方根誤差(RMSE)(分別為0.159和0.259),且相對分析誤差(RPD)>2(分別為2.9、2.094),具有非常好的預(yù)測能力。夏季土壤呼吸最佳預(yù)測模型則為包含10 cm處氣溫、地溫、空氣相對濕度和土壤濕度所建的4因子方程,RMSE為0.248,RPD>2(為2.406),可用于精確預(yù)測。(3)各季節(jié)土壤呼吸變化趨勢與其影響因子的變化,因春季的完全同步,夏季基本一致,而秋季一致性較差,故春季土壤呼吸最佳預(yù)測模型的預(yù)測精度最高(92.67%),夏季次之(84.99%),秋季較差(77.23%)。

    強堿土;土壤呼吸;影響因子;預(yù)測模型

    土壤呼吸是陸地生態(tài)系統(tǒng)碳循環(huán)的重要環(huán)節(jié),是土壤碳素向大氣輸出的主要途徑[1-2]。土壤碳庫的較小變動即會顯著改變大氣中CO2濃度和土壤碳沉積速率[3]。干旱半干旱區(qū)生態(tài)環(huán)境脆弱,具有極端的非生物生態(tài)因子,該區(qū)域土壤有機碳含量偏低,土壤呼吸對溫度與水分變化極為敏感,因此成為衡量區(qū)域氣候變化的重要指標之一[4]。對于土壤呼吸的影響因子研究較多,尤其以溫度和水分較為突出,了解這些因子與土壤呼吸之間的關(guān)系,對于區(qū)域碳循環(huán)的估計和預(yù)測有重要意義。溫度是調(diào)節(jié)和控制許多生態(tài)學(xué)過程的關(guān)鍵因素,也是影響土壤呼吸的決定性因子之一。土壤呼吸與溫度間的擬合方式較多,主要有線性方程[5]、指數(shù)方程[6]、Arrhenius方程[7]、冪函數(shù)方程[8]和邏輯斯蒂方程[9]。水分對土壤呼吸的影響也很重要,與土壤呼吸存在著不同的相關(guān)關(guān)系,擬合方式主要有線性方程[10-12]、三次曲線模型[13]等。對干旱半干旱區(qū)土壤呼吸及影響因子的研究多集中在森林、沙漠、農(nóng)田等[14-17],而強堿土土壤呼吸的影響因子分析及建立估測模型的相關(guān)研究相對較少。本文以干旱區(qū)強堿土為研究對象,通過分析春、夏、秋季土壤呼吸與其同步所測定的溫濕度氣象因子的相關(guān)關(guān)系,找出不同季節(jié)強堿土土壤呼吸日變化的主要影響因子,以此建立估測模型,并以簡單易行為原則,選出各季節(jié)最佳模型,為干旱區(qū)碳儲量的快速估測提供參考依據(jù)。

    1 研究區(qū)概況

    研究區(qū)位于天山山脈東段北坡與準噶爾盆地東南緣地區(qū)過渡地帶的新疆奇臺縣草原站(44°04′42.62″N,89°53′24.81″E),海拔791 m,處于中部沖積平原區(qū),屬于中溫帶大陸性干旱半干旱氣候,年平均氣溫4.7℃,年平均降水量176 mm,蒸發(fā)潛力2 141 mm,年日照時數(shù)2 280~3 230 h。這里存在大面積堿化土壤,在準噶爾盆地東南緣范圍內(nèi)具有很強的代表性。

    2 材料與方法

    2.1 樣地選擇及土壤性質(zhì)

    奇臺縣平原區(qū)東部草原站土壤堿化情況嚴重且面積較大,在此區(qū)域內(nèi)采用梅花樁采樣法,設(shè)5個樣點,取0~20 cm土層土樣,混合均勻,得到均一土樣。其所含基本成分由中國科學(xué)院新疆生態(tài)與地理研究所專業(yè)土壤分析人員測定,基本土壤性質(zhì)如下:pH值范圍在10以上,有機質(zhì)含量極低,約為2.5 g·kg-1,電導(dǎo)率的范圍在3.34~4.57 ms·cm-1之間(表1)。根據(jù)中國土壤地理發(fā)生學(xué)分類和系統(tǒng)分類近似參比[18]及堿化土壤分級相關(guān)研究結(jié)果[19],樣區(qū)土壤屬于強堿土(堿積鹽成土),本文以此為研究對象,在地勢平坦,避開農(nóng)田、居民點和道路等人為影響的區(qū)域進行試驗。

    表1 土壤基本性質(zhì)Table 1 The properties of tested soil

    2.2 土壤呼吸和溫濕度的測定

    在樣地內(nèi)設(shè)置一個土壤呼吸環(huán)(因5個樣點的土壤理化性質(zhì)差異不明顯,設(shè)置1個呼吸環(huán)可代表整個研究對象),即直徑20 cm,高20 cm的PVC管,將其置入土中,露出地表約3 cm。為避免安置PVC管對土壤呼吸速率造成短時間的擾動,需靜置24 h待其平衡。利用開路式土壤碳通量測量儀LI-8100(LI-COR,USA)每1 h自動測1次土壤呼吸速率(每次測2個值,取其平均值),進行24 h動態(tài)觀測,共測定3 d。同時用Testo625溫濕度計同步記錄了近地表氣溫與空氣相對濕度(10 cm處),10 cm處土壤濕度和地溫的數(shù)據(jù)則分別由LI-8100儀器所帶的ML2x水分和E型熱電偶溫度傳感器測得。試驗選取春、夏、秋季中具有代表性的4、8、10月進行,且期間均為晴朗天氣。

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

    由于各次試驗期間每天之間的溫度和濕度的變化相差不大,因此將前2天的觀測數(shù)據(jù)平均成一天,進行建模,第3天的數(shù)據(jù)作為檢驗樣本,對土壤呼吸速率與各因子所建模型進行檢驗。數(shù)據(jù)整理和作圖在Microsoft Excel 2003中進行,分析在SPSS19.0和DPS7.5中完成。

    2.4 通徑分析

    通徑分析是應(yīng)用通徑系數(shù)分析方法,在相關(guān)分析與回歸分析基礎(chǔ)上,進一步研究因變量與自變量之間的數(shù)量關(guān)系,并可將相關(guān)系數(shù)分解為直接作用系數(shù)和間接作用系數(shù),以揭示各個因素對因變量的相對重要性。通徑系數(shù)是從簡單的相關(guān)系數(shù)矩陣開始,通過求解通徑系數(shù)的標準化正規(guī)方程,進而求出直接通徑系數(shù)和間接通徑系數(shù)[20]。

    設(shè)在p個自變量x1,x2,x3,…,xp中,每兩個變量之間與因變量y之間的簡單相關(guān)系數(shù)可以構(gòu)成求解通徑系數(shù)的標準化正規(guī)方程:

    其中,ρ1,ρ2,…,ρp為直接通徑系數(shù)。直接通徑系數(shù)可以通過求上述相關(guān)矩陣的逆矩陣而獲得。假設(shè)Cij為相關(guān)矩陣rij的逆矩陣,那么直接通徑系數(shù)ρi(i=1,2,…,p)為

    直接通徑系數(shù)ρi常表示為ρi→y;間接通徑系數(shù)ρi→j→y可以通過相關(guān)系數(shù)rij和直接通徑系數(shù)ρj→y來計算,ρi→j→y=rijρj→y。

    2.5 模型檢驗

    通過以下3個檢驗參數(shù)來判定模型的精度:(1)均方根誤差(RMSE):其值越小,表明實際值和模型預(yù)測值之間的偏差越小,模型的精度越高。計算公式為:

    其中,n為樣本個數(shù),yi和分別代表實測值和預(yù)測值。

    (2)相對分析誤差(RPD):當RPD>2時,表明模型具有非常好的預(yù)測能力;當2>RPD>1.4時,說明模型可對樣品進行粗略估計,RPD<1.4時,則表明模型無法對樣品進行預(yù)測[21]。公式如下:

    公式中,S.D.為樣本標準差。

    (3)模型判定系數(shù)(R2):即模型預(yù)測值和實測值的相關(guān)系數(shù)的平方來判定模型的精度,R2越高,表明預(yù)測模型精度越高。其公式如下:

    3 結(jié)果與分析

    3.1 溫度和濕度與土壤呼吸的相關(guān)關(guān)系

    將三個季節(jié)強堿土土壤呼吸速率與同步所測得的溫度、濕度進行相關(guān)分析(表2)。春、秋季土壤呼吸均與氣溫、空氣相對濕度和土壤濕度達到了極顯著相關(guān),且前者與地溫的相關(guān)性顯著,而后者與地溫無顯著相關(guān)性。兩季節(jié)相關(guān)性大?。ń^對值)分別依次為氣溫>空氣相對濕度>土壤濕度>地溫、氣溫>土壤濕度>空氣相對濕度>地溫。氣溫為其土壤呼吸變化的最主要影響因子(0.923和0.669)。同時,它們除與空氣相對濕度為負相關(guān)外,與其它3因子均為正相關(guān)。

    夏季土壤呼吸與氣溫、地溫、空氣相對濕度和土壤濕度均達到了極顯著相關(guān)關(guān)系,同時,與空氣相對濕度為正相關(guān),而與其它3個因子為負相關(guān)。這與春、秋季與其影響因子的關(guān)系相反。此外,以土壤呼吸與地溫的相關(guān)系數(shù)為最大(0.748),其余依次為土壤濕度、空氣相對濕度、氣溫。

    上述相關(guān)分析表明,春、秋季強堿土土壤呼吸受氣溫的影響最大,而夏季則受地溫的影響最大。此外,不論哪個季節(jié)強堿土土壤呼吸受濕度的影響也較大(相關(guān)性均達極顯著)。

    3.2 土壤呼吸影響因子的通徑分析

    為進一步分析各季節(jié)土壤呼吸與其影響因子的關(guān)系,使用通徑分析方法,對比其通徑系數(shù)絕對值的大小來找出最大的直接和間接影響因子(表3)。氣溫對強堿土春、秋季節(jié)土壤呼吸的直接作用最大(1.5735、2.3722),而對夏季土壤呼吸直接作用最大的因子是地溫(-2.7964)。

    對比間接通徑系數(shù),春、秋季的土壤濕度與空氣相對濕度的間接作用均較大(-1.4809和1.5119、-2.2611和2.3275),且均通過氣溫起間接作用。而夏季氣溫、空氣相對濕度、土壤濕度的間接作用亦較大(-2.4114、2.5718、-2.7123),均通過地溫來影響。土壤濕度是影響春、夏、秋季土壤呼吸的最大間接因子。此外,各季節(jié)所選的這些較大間接通徑系數(shù)均與其最大直接通徑系數(shù)在同一量級,表明其它因子的間接作用對于土壤呼吸的影響也大。

    表2 溫度、濕度與土壤呼吸的相關(guān)關(guān)系Table 2 Correlation coefficients between soil respiration rate,temperature and humidity

    表3 各季節(jié)土壤呼吸影響因子的通徑分析Table 3 Path analysis on soil respiration influencing factors in each season

    通徑分析結(jié)果與上述相關(guān)分析的結(jié)論一致,更進一步表明溫度因子是各季節(jié)強堿土的最主要影響因素,但溫濕度的耦合作用對于它的影響也不可忽略。

    3.3 模型建立

    3.3.1 單因子最優(yōu)模型選取各季節(jié)與土壤呼吸的相關(guān)系數(shù)最大的因子建立線性、冪函數(shù)、指數(shù)、負指數(shù)、雙曲線、對數(shù)、二階、三階函數(shù)模型進行對比,在此僅列出R2≥0.5,且P為0.000的模型(表4)。春季的Yɑ1、Yɑ6、Yɑ7方程的R2均達到了0.86以上,但后兩者比前者僅大0.011和0.022,其差距不明顯,且方程更復(fù)雜,所以Yɑ1方程為春季土壤呼吸預(yù)測的最佳單因子模型。同理,夏季的Yb1和秋季的Yc1方程,均具有較大的R2,且更簡單易操作,所以二者為兩季節(jié)土壤呼吸預(yù)測的最佳單因子模型。同時,這反映出氣溫不僅是影響春、秋季強堿土土壤呼吸的最主要因子并且兩者呈線性關(guān)系,而夏季強堿土的土壤呼吸最主要影響因子則為地溫。上述單因子模型,春季的R2達到0.8,夏、秋季的R2均達到0.5以上,雖然溫度是影響各季節(jié)土壤呼吸的主要因素,但是夏、秋季模型R2仍然較低??梢姖穸纫蛩氐淖饔貌豢珊雎?,土壤呼吸是受溫、濕度共同控制的,需將二者共同引入建立方程。

    表4 單因子土壤呼吸預(yù)測模型Table 4 Single factor forecastingmodel for soil respiration

    3.3.2 多因子最優(yōu)模型利用線性擬合分析將各季節(jié)與土壤呼吸相關(guān)性達到極顯著水平的因子進行多因子建模(表5)。因變量為實測土壤呼吸速率,自變量為同步所測的氣溫、地溫、空氣相對濕度和土壤濕度。所有模型均通過P<0.01檢驗,即模型均有效。

    春、秋季所建的模型中,除Yɑ9、Yc7模型R2(0.779和0.415)明顯小于兩季節(jié)所選的最優(yōu)單因子方程R2(0.867和0.505)外,其余方程的R2隨建模因子的增多而明顯提高??梢姡蓾穸纫蜃咏⒌哪P途炔?,表明濕度不是影響土壤呼吸的主要因子,而氣溫是最主要的影響因子,它與濕度共同建立的模型更好,即精度更高。

    夏季的Yb7、Yb8、Yb11、Yb14方程R2低于土壤呼吸與地溫所建的最佳單因子方程R2(0.572),而其余方程R2則隨建模因子的增多而明顯提高??梢姡瑲鉁睾蜐穸纫蜃咏?,以及濕度因子獨立建模,精度均較差。地溫因子參與建模,方程精度均較高,表明其是影響夏季土壤呼吸的最主要因子,且加入其它因素后模型更好。

    將春、夏、秋季R2最高的多因子模型(Yɑ11、Yb16、Yc8)與所選最佳單因子模型(Yɑ1、Yb1、Yc1)對比,從R2的大小來看,多因子模型比單因子模型分別高0.12、0.279、0.309。這反映出,不論哪個季節(jié)溫、濕度的耦合控制了強堿土的土壤呼吸。

    3.4 模型選擇與精度檢驗

    將選取的最優(yōu)單因子、多因子模型進行精度檢驗并對比分析(表6)。春季所選的3個方程,其RPD均大于2,都可對強堿土土壤呼吸做精確的預(yù)測。從R2來看,Yɑ1方程擬合優(yōu)度明顯較低,而另外兩模型的R2較大,且相差不大,但是雙因子Yɑ10模型更簡便易行,RMSE最?。?.159),RPD最大(2.9),其有更好的預(yù)測能力,是春季土壤呼吸日變化的最佳預(yù)測模型。此外,Yɑ1方程反映出氣溫是決定春季強堿土土壤呼吸估測模型精度的最主要因子,也可以直接用來預(yù)測春季土壤呼吸。Yɑ10方程在氣溫的建?;A(chǔ)上加入了土壤濕度因子,模型預(yù)測精度進一步提升,但在此基礎(chǔ)上加入了空氣相對濕度因子的Yɑ11方程,其預(yù)測精度(RPD)降低,這因為空氣相對濕度的日際變化大,沒有氣溫和土壤濕度的日際變化穩(wěn)定,它的參與降低了預(yù)測模型的精度。

    夏季的Yb1、Yb6、Yb13、Yb16方程較優(yōu),其中地溫是決定夏季強堿土土壤呼吸預(yù)測模型精度的最主要影響因子,但是只能做粗略的估測。隨著建模因子的增加,R2不斷增大(分別為0.572、0.643、0.768、0.851)。與前者的R2相比,后三者方程R2分別增加了0.071、0.196、0.279。從4個方程的精度檢驗結(jié)果可看出,RMSE隨建模因子的增加不斷變小,且RPD不斷增大,4個因子共建的方程Yb16具有最小的RMSE和最大的RPD(2.406),其是夏季土壤呼吸日變化的最優(yōu)模型,可做精確預(yù)測。這表明除地溫外,其它3個因子的間接作用也較大,這與通徑分析結(jié)果一致。

    表5 各季節(jié)土壤呼吸與溫、濕度模型Table 5 Themodel of soil respiration and temperature and humidity for seasons

    秋季3個模型的R2相比,Yc6、Yc8比Yc1分別增加了0.232和0.309,但前者比后者僅小0.077,同時其RMSE最小,RPD最大,為2.094,且其為雙因子模型,更簡單,所以選擇Yc6方程為秋季土壤呼吸日變化的最佳估測模型。Yc1所建立的最主要影響因子(氣溫)模型,不能對土壤呼吸進行粗略預(yù)測,但Yc6模型中加入土壤濕度因子,就具有很高的預(yù)測精度,說明土壤濕度對于秋季土壤呼吸有重要的作用。Yc8模型在加入空氣相對濕度因子后,其預(yù)測精度較Yc6差,這與春季相同。

    3.5 各季節(jié)最佳預(yù)測模型對比及影響因子分析

    將春、夏、秋季土壤呼吸日變化的最佳預(yù)測模型Yɑ10、Yb16和Yc6所計算的預(yù)測值與實測值作圖以及方程擬合(圖1)??煽闯觯杭镜膶崪y與預(yù)測土壤呼吸值吻合最好,其方程的決定系數(shù)R2最大(0.9267),夏季次之(R2=0.8499),秋季較差(R2=0.7723)。

    由于春季融雪及凍土解融后,土壤濕度大,有利于表層土壤碳酸鹽(CaCO3、Na2CO3等)吸收CO2,其反應(yīng)生成的重碳酸鹽(Ca(HCO3)2、NaHCO3等)在白天氣溫升高(氣溫低于20℃)狀態(tài)下分解釋放,而夜間(9∶00—次日8∶00)低溫下又緩慢吸收,土壤呼吸日變化與氣溫、土壤濕度日變化的趨勢完全一致(均在13∶00左右達到最大,夜間一致降低),線性方程可很好地擬合它們之間的關(guān)系,模型精度最高。

    夏季氣溫、地溫、土壤濕度變化曲線與土壤呼吸變化相反(3因子均在16∶00左右達到最大值,而土壤呼吸早1 h達到最小值,夜間時段其隨三者的降低而增大),而空氣相對濕度變化則與其變化基本一致。在9∶00—15∶00時段土壤呼吸隨地溫的升高而降低(因溫度高,其氣溫、地溫分別在35℃、30℃左右,土壤濕度在11%左右,高溫加速了碳酸鹽與CO2的化學(xué)反應(yīng)速率,吸收CO2),而16∶00—21∶00時段則隨溫度的降低而增大(溫度降低,使得土壤中化學(xué)反應(yīng)變慢,消耗CO2的量減少)。各影響因子與土壤呼吸的變化大致同步,故模型的精度較低。

    而秋季土壤呼吸日變化在11∶00—15∶00時段與氣溫、土壤濕度變化完全相反,但在其它時段卻變化一致(此時段氣溫較高,其接近30℃,加之土壤濕度較大,有利于近地表土壤中碳酸鹽與CO2的反應(yīng)吸收)。其影響因子與土壤呼吸的變化一致性最差,模型精度也最低。

    圖1 實測與預(yù)測土壤呼吸速率的對比及模型擬合(a1、a2:春季;b1、b2:夏季;c1、c2:秋季)Fig.1 The comparisons and model fitting of soil respiration rate between themeasured and predicted values(a1,a2:Spring;b1,b2:Summer;c1,c2:Autumn)

    4 結(jié)論

    1)春、秋季強堿土土壤呼吸與氣溫、空氣相對濕度和土壤濕度均達到了極顯著相關(guān),與氣溫的相關(guān)性最強。夏季土壤呼吸與溫、濕度四因子均存在極顯著的相關(guān)關(guān)系,與地溫的相關(guān)性最高。通徑分析表明,氣溫是春、秋季土壤呼吸的最大直接影響因子,夏季則為地溫因子,且均以土壤濕度間接作用最大。溫、濕度的耦合控制土壤呼吸變化。

    2)各季節(jié)土壤呼吸與相關(guān)性最大因子所建的8種模型中,以線性模型為最佳,均具有較大R2,且簡便易操作。隨著建模因子的增加,多因子模型R2基本上不斷增大,但單獨的濕度因子所建模型精度均較差。

    表6 土壤呼吸估測模型的精度檢驗Table 6 The precision test of soil respirationpredictingmodels

    3)在簡單易行的原則下,春、秋季土壤呼吸的最佳預(yù)測模型均為10 cm處氣溫和土壤濕度所建的雙因子方程。夏季最佳預(yù)測模型則為與10 cm處氣溫、地溫、空氣相對濕度和土壤濕度所建的4因子方程。

    4)各季節(jié)土壤呼吸最佳預(yù)測模型,春季的預(yù)測精度最高(R2=0.9267),夏季次之(R2=0.8499),秋季較差(R2=0.7723)。

    [1]Schimel DS.Terrestrial ecosystems and the carbon cycle[J].Global Change Biology,1995,1(1):77-99.

    [2]SchlesingerWH,Andrews JA.Soil respiration and the global carbon cycle[J].Biogeochemistry,2000,48(1):7-20.

    [3]Davidson EA,Trumbore SE,Amundson R.Biogeochemistry:soil warming and organic carbon content[J].Nature,2000,408(6814):789-790.

    [4]王新源,李玉霖,趙學(xué)勇,等.干旱半干旱區(qū)不同環(huán)境因素對土壤呼吸影響研究進展[J].生態(tài)學(xué)報,2012,32(15):4890-4901.

    [5]Chimner RA.Soil respiration rates of tropical peat lands in Micronesia and Hawaii[J].Wetlands,2004,24(1):51-56.

    [6]Buchmann N.Biotic and abiotic factors controlling soil respiration rates in Picea abies stands[J].Soil Biology and Biochemistry,2000,32(11/12):1625-1635.

    [7]Thierron V,Laudelout H.Contribution of root respiration to total CO2efflux from the soil of a deciduous forest[J].Canadian Journal of Forest Research,1996,26(7):1142-1148.

    [8]Fang C,Moncrieff JB.The dependence of soil CO2efflux on temperature[J].Soil Biology and Biochemistry,2001,33(2):155-165.

    [9]Rodeghiero M,Cescatti A.Main determinants of forest soil respiration along an elevation temperature gradient in the Italian Alps[J].Global Change Biology,2005,11(7):1024-1041.

    [10]黃湘,李衛(wèi)紅,陳亞寧,等.塔里木河下游荒漠河岸林群落土壤呼吸及其影響因子[J].生態(tài)學(xué)報,2007,27(05):1951-1959.

    [11]張丁辰,蔡典雄,代快,等.旱作農(nóng)田不同耕作土壤呼吸及其對水熱因子的響應(yīng)[J].生態(tài)學(xué)報,2013,33(6):1916-1925.

    [12]康永祥,夏國威,劉建軍,等.秦嶺小隴山銳齒櫟林皆伐跡地土壤呼吸特征[J].應(yīng)用生態(tài)學(xué)報,2014,25(2):342-350.

    [13]黃湘,陳亞寧,李衛(wèi)紅,等.塔里木河中下游檉柳群落土壤碳通量及其影響因子分析[J].環(huán)境科學(xué),2006,27(10):1934-1940.

    [14]張麗萍,陳亞寧,李衛(wèi)紅,等.準噶爾盆地鹽穗木群落土壤CO2釋放規(guī)律及其影響因子[J].干旱區(qū)研究,2007,24(6):6854-6860.

    [15]楊建軍,呂光輝,張燕,等.艾比湖流域不同植物群落土壤呼吸研究[J].環(huán)境科學(xué)研究,2009,22(3):362-370.

    [16]周洪華,李衛(wèi)紅,楊余輝,等.干旱區(qū)不同土地利用方式下土壤呼吸日變化差異及影響因素[J].地理科學(xué),2011,31(2):190-196.

    [17]張麗華,陳亞寧,李衛(wèi)紅,等.干旱區(qū)荒漠生態(tài)系統(tǒng)的土壤呼吸[J].生態(tài)學(xué)報,2008,28(5):1911-1922.

    [18]海春興,陳健飛.土壤地理學(xué)[M].北京:科學(xué)出版社,2010.

    [19]張芳,熊黑鋼,安放舟,等.基于鹽(堿)生植被蓋度的土壤堿化分級[J].土壤學(xué)報,2012,49(4):665-672.

    [20]唐啟義.DPS數(shù)據(jù)處理系統(tǒng)[M].北京:科學(xué)出版社,2010.

    [21]Chang CW,Laird DA.Near-infrared reflectance spectroscopic analysis of soil C and N[J].Soil Science,2002,167(2):110-116.

    Impact factor analysis and model prediction of strong alkaline soil respiration in different seasons

    ZHANG Lu1,ZHANG Fang1,2,XIONG Hei-gang1,3,DUAN Peng-cheng1,LIRong-rong1,ZHANG Nan1
    (1.College of Resourcesɑnd Environmentɑl Science,Xinjiɑng University/Key Lɑborɑtory of Oɑsis Ecology ofMinistry of Educɑtion,Urumqi,Xinjiɑng 830046,Chinɑ;2.Ecology Post-doctorɑl Reseɑrch Stɑtion,Xinjiɑng University,Urumqi,Xinjiɑng 830046,Chinɑ;3.College of Artsɑnd Sciencesof Beijing Union University,Beijing 100191,Chinɑ)

    In this paper,strong alkaline soil respiration rate has beenmeasured in spring,summer,autumn under sunny weather condition by using LI-8100 soil carbon fluxmeter,temperature(air temperature and soil temperature),humidity(relative humidity and soilmoisture)were also recorded simultaneously,and the correlation among above factors has been analyzed to identify the key factors influencing soil respiration in different season,and establish a variety of regressionmodels.On the basis of the accuracy test and simple principle,the optimalmodel to predict the seasonal soil respiration has been selected.The results show that:(1)Although the temperature and humidity are themain factors affecting strong alkaline soil respiration in different seasons,but the temperature has relatively larger effect,where the air temperature is the biggestdirect impact factor on soil respiration in spring and autumn,while in summer,the soil temperature is the greatest direct impact factor,and soilmoisture is the largest indirect factor in all season.(2)The best predictionmodel of soil respiration in spring and autumn is the two-factor equation composed of air temperatures at 10 cmlayer and soilmoisture,which has a smaller RMSE(0.159 and 0.259,respectively),and the RPD>2(2.9,2.094,respectively),has a better predictive accuracy.The best soil respiration predictivemodel for summer season is composed of air temperature at 10 cm,surface temperature,relative humidity and soilmoisture,RMSE is 0.248,RPD>2(of 2.406).(3)The trend of soil respiration variation and their impact factorswas similar in spring,but different in summer and autumn.In spring the soil respiration model has the highest prediction accuracy(92.67%),followed by summer(84.99%)and poor in fall(77.23%).

    strong alkaline soil;soil respiration;impact factor;prediction model

    S154.1

    :A

    1000-7601(2017)01-0071-08

    10.7606/j.issn.1000-7601.2017.01.11

    2016-01-25

    國家自然科學(xué)基金(41261049、41171165);中國博士后科學(xué)基金(2013M532100);北京市屬高等學(xué)校高層次人才引進與培養(yǎng)計劃項目(IDHT20130322)

    張錄(1989—),男,陜西漢中人,碩士,主要從事干旱區(qū)環(huán)境研究。E-mail:zhanglu1356@126.com。

    張芳(1969—),女,山東嘉祥人,副教授,主要從事干旱區(qū)環(huán)境與人地關(guān)系的研究。E-mail:zhangf602@163.com。

    猜你喜歡
    強堿通徑土壤濕度
    牡丹江市氣象因子與PM2.5濃度影響通徑分析
    黑龍江氣象(2021年2期)2021-11-05 07:06:58
    新型全自動氣動通徑機設(shè)計要點探究
    氣動通徑機設(shè)計及使用中應(yīng)注意的幾個問題
    土壤濕度傳感器在園林綠化灌溉上的應(yīng)用初探
    基于51單片機控制花盆土壤濕度
    電子制作(2019年15期)2019-08-27 01:12:12
    四川盆地土壤濕度時空分布及影響因子分析
    中國不同氣候區(qū)土壤濕度特征及其氣候響應(yīng)
    一種高強度耐強酸強堿環(huán)氧樹脂密封膠
    強堿混合溶液pH值的計算
    新型鋼管氣動通徑機的研發(fā)
    美女免费视频网站| 最近最新免费中文字幕在线| 久久欧美精品欧美久久欧美| netflix在线观看网站| 真实男女啪啪啪动态图| 国产成人aa在线观看| 久久久久国内视频| 日韩欧美免费精品| 黄色视频,在线免费观看| 亚洲熟妇中文字幕五十中出| АⅤ资源中文在线天堂| 五月玫瑰六月丁香| 国产主播在线观看一区二区| 午夜福利在线在线| 欧美日韩黄片免| 一进一出抽搐gif免费好疼| 在线天堂最新版资源| 在线观看66精品国产| 亚洲专区国产一区二区| 免费无遮挡裸体视频| 国产伦在线观看视频一区| 脱女人内裤的视频| 精品国内亚洲2022精品成人| 一本精品99久久精品77| 久久伊人香网站| 丰满人妻熟妇乱又伦精品不卡| 亚洲av成人精品一区久久| 亚洲天堂国产精品一区在线| 日韩大尺度精品在线看网址| 国产伦精品一区二区三区四那| 在线播放无遮挡| 欧美日韩瑟瑟在线播放| 日本熟妇午夜| 亚洲经典国产精华液单 | 国产精品亚洲一级av第二区| 亚洲精品在线美女| 此物有八面人人有两片| 久久久久亚洲av毛片大全| 国产野战对白在线观看| 男女做爰动态图高潮gif福利片| 香蕉av资源在线| 美女黄网站色视频| 久久久久久久午夜电影| 欧美性猛交╳xxx乱大交人| 性色av乱码一区二区三区2| 嫩草影院新地址| 日本 欧美在线| 午夜福利视频1000在线观看| 99精品久久久久人妻精品| 一进一出抽搐动态| 日韩有码中文字幕| 国产精品嫩草影院av在线观看 | 亚洲美女黄片视频| 亚洲成人中文字幕在线播放| 亚洲自拍偷在线| 亚洲成人中文字幕在线播放| 精品人妻1区二区| av天堂在线播放| 人人妻,人人澡人人爽秒播| 亚洲乱码一区二区免费版| 国产精品,欧美在线| 国产精品三级大全| 久久久久久大精品| 久久久精品欧美日韩精品| 国产精品自产拍在线观看55亚洲| 免费看美女性在线毛片视频| 欧美日韩乱码在线| 精品久久久久久久末码| 97热精品久久久久久| 亚洲欧美清纯卡通| 日本一二三区视频观看| 99riav亚洲国产免费| 免费在线观看亚洲国产| 99在线视频只有这里精品首页| 色精品久久人妻99蜜桃| 日韩亚洲欧美综合| 国产人妻一区二区三区在| 欧美中文日本在线观看视频| 免费黄网站久久成人精品 | 美女免费视频网站| 日韩欧美一区二区三区在线观看| 国产免费av片在线观看野外av| 黄片小视频在线播放| 女生性感内裤真人,穿戴方法视频| 舔av片在线| 亚洲欧美日韩卡通动漫| 18禁裸乳无遮挡免费网站照片| 亚洲成av人片免费观看| 可以在线观看的亚洲视频| 麻豆成人午夜福利视频| 又爽又黄无遮挡网站| 真实男女啪啪啪动态图| 99久久精品一区二区三区| 亚洲,欧美,日韩| .国产精品久久| 精品一区二区三区人妻视频| 午夜福利高清视频| 自拍偷自拍亚洲精品老妇| 国产伦精品一区二区三区视频9| 亚洲精品色激情综合| 精品久久久久久久久av| 国产精品精品国产色婷婷| 欧美激情在线99| 亚洲欧美日韩卡通动漫| 欧美日韩福利视频一区二区| 一级黄片播放器| 久久久久性生活片| 中出人妻视频一区二区| 丰满人妻熟妇乱又伦精品不卡| 一卡2卡三卡四卡精品乱码亚洲| 中文字幕人妻熟人妻熟丝袜美| 99久久99久久久精品蜜桃| 久久精品国产亚洲av涩爱 | 国产乱人伦免费视频| 一个人看的www免费观看视频| 一级黄色大片毛片| 在线观看av片永久免费下载| 国产男靠女视频免费网站| 亚洲片人在线观看| 88av欧美| 51午夜福利影视在线观看| 99riav亚洲国产免费| 国产人妻一区二区三区在| 99热这里只有精品一区| av在线老鸭窝| 男女那种视频在线观看| 色综合站精品国产| 国产国拍精品亚洲av在线观看| 九色成人免费人妻av| 亚洲国产精品成人综合色| 无人区码免费观看不卡| 亚洲久久久久久中文字幕| 午夜视频国产福利| 99久久九九国产精品国产免费| 在线天堂最新版资源| 人人妻人人澡欧美一区二区| 午夜a级毛片| 综合色av麻豆| 亚洲精品影视一区二区三区av| www.www免费av| 亚洲欧美激情综合另类| 级片在线观看| 韩国av一区二区三区四区| 久久香蕉精品热| 久久热精品热| 九九久久精品国产亚洲av麻豆| 久久人人爽人人爽人人片va | 久久久色成人| 熟女人妻精品中文字幕| 校园春色视频在线观看| 国产av麻豆久久久久久久| 欧美精品啪啪一区二区三区| 12—13女人毛片做爰片一| 18禁黄网站禁片免费观看直播| 久久久久久国产a免费观看| 国产精品久久视频播放| 无人区码免费观看不卡| 色综合欧美亚洲国产小说| 中文在线观看免费www的网站| 好看av亚洲va欧美ⅴa在| 人妻夜夜爽99麻豆av| 午夜福利18| 深夜精品福利| 最近视频中文字幕2019在线8| 日本一二三区视频观看| 婷婷丁香在线五月| 亚洲av中文字字幕乱码综合| 免费无遮挡裸体视频| 久久欧美精品欧美久久欧美| 91麻豆精品激情在线观看国产| 男女之事视频高清在线观看| 国产精品美女特级片免费视频播放器| 亚洲久久久久久中文字幕| 午夜日韩欧美国产| 一级黄色大片毛片| 精品人妻视频免费看| 日本精品一区二区三区蜜桃| 免费大片18禁| 18禁黄网站禁片午夜丰满| 色综合欧美亚洲国产小说| av在线老鸭窝| 免费无遮挡裸体视频| 国产野战对白在线观看| 一二三四社区在线视频社区8| 精品久久久久久久末码| 桃色一区二区三区在线观看| 九色国产91popny在线| 日本熟妇午夜| 好男人电影高清在线观看| 美女黄网站色视频| 欧美黑人巨大hd| x7x7x7水蜜桃| 99热精品在线国产| 亚洲精品一卡2卡三卡4卡5卡| 久久欧美精品欧美久久欧美| 一进一出好大好爽视频| 亚洲av第一区精品v没综合| 亚洲人成电影免费在线| 精品久久久久久成人av| 亚洲午夜理论影院| 国产精品综合久久久久久久免费| 91麻豆精品激情在线观看国产| а√天堂www在线а√下载| 熟女人妻精品中文字幕| 欧美高清性xxxxhd video| 91九色精品人成在线观看| 国产爱豆传媒在线观看| 国产欧美日韩一区二区三| 一区二区三区四区激情视频 | 成人鲁丝片一二三区免费| 国产av麻豆久久久久久久| 动漫黄色视频在线观看| 欧美在线黄色| 舔av片在线| 日本撒尿小便嘘嘘汇集6| 天堂av国产一区二区熟女人妻| 成人毛片a级毛片在线播放| 一二三四社区在线视频社区8| 男人的好看免费观看在线视频| 亚洲人成网站在线播| 天天一区二区日本电影三级| 搡老岳熟女国产| 午夜亚洲福利在线播放| 国产中年淑女户外野战色| 国产av麻豆久久久久久久| 国产午夜福利久久久久久| 亚州av有码| 久久久久九九精品影院| 久久久久精品国产欧美久久久| 成年版毛片免费区| 国产主播在线观看一区二区| 两个人视频免费观看高清| 欧美极品一区二区三区四区| 亚洲va日本ⅴa欧美va伊人久久| 国产野战对白在线观看| 黄色视频,在线免费观看| 在线免费观看不下载黄p国产 | 国产午夜精品论理片| 人妻夜夜爽99麻豆av| 精品久久国产蜜桃| bbb黄色大片| 亚洲国产精品合色在线| 夜夜看夜夜爽夜夜摸| 日韩欧美三级三区| 美女cb高潮喷水在线观看| 亚洲天堂国产精品一区在线| 免费av不卡在线播放| 久久久成人免费电影| 日本 欧美在线| 内地一区二区视频在线| 天天一区二区日本电影三级| 免费看a级黄色片| 听说在线观看完整版免费高清| 熟女人妻精品中文字幕| 人人妻人人看人人澡| 精品久久久久久,| 女人十人毛片免费观看3o分钟| 又爽又黄无遮挡网站| 757午夜福利合集在线观看| 日本 av在线| 一区二区三区免费毛片| 精品99又大又爽又粗少妇毛片 | 激情在线观看视频在线高清| 美女免费视频网站| 免费搜索国产男女视频| 亚洲精品乱码久久久v下载方式| 亚洲最大成人中文| 婷婷色综合大香蕉| 日本a在线网址| 欧美极品一区二区三区四区| 欧美精品啪啪一区二区三区| 黄色女人牲交| 日韩欧美精品免费久久 | 美女高潮喷水抽搐中文字幕| 免费大片18禁| 欧美成人免费av一区二区三区| 亚洲欧美日韩东京热| 亚洲成人精品中文字幕电影| 淫秽高清视频在线观看| 国产精品久久视频播放| 乱码一卡2卡4卡精品| 搡老熟女国产l中国老女人| 2021天堂中文幕一二区在线观| 国产高清有码在线观看视频| 精品人妻一区二区三区麻豆 | 啪啪无遮挡十八禁网站| 欧美国产日韩亚洲一区| 国产精品98久久久久久宅男小说| 国产成年人精品一区二区| 精品欧美国产一区二区三| 啦啦啦韩国在线观看视频| 中文字幕av成人在线电影| 偷拍熟女少妇极品色| 校园春色视频在线观看| 久久精品国产亚洲av涩爱 | 波多野结衣高清作品| 久久久久久久久久黄片| 深夜精品福利| 99riav亚洲国产免费| 国产成人a区在线观看| 夜夜躁狠狠躁天天躁| 欧美最黄视频在线播放免费| 黄色日韩在线| 亚洲熟妇中文字幕五十中出| 欧美3d第一页| 日本黄色视频三级网站网址| 最近视频中文字幕2019在线8| 免费黄网站久久成人精品 | 亚洲成av人片免费观看| 亚洲avbb在线观看| 欧美一级a爱片免费观看看| 亚洲三级黄色毛片| 中文字幕免费在线视频6| 精品无人区乱码1区二区| 天堂影院成人在线观看| 国产视频一区二区在线看| 他把我摸到了高潮在线观看| 美女cb高潮喷水在线观看| 国产在视频线在精品| 久久精品91蜜桃| 国产精品免费一区二区三区在线| 久久久国产成人免费| 午夜视频国产福利| 在线观看av片永久免费下载| 亚洲熟妇熟女久久| 午夜精品久久久久久毛片777| 丁香六月欧美| 一进一出抽搐gif免费好疼| 亚洲精品粉嫩美女一区| 欧美日韩福利视频一区二区| 精品久久久久久久人妻蜜臀av| 国产精品综合久久久久久久免费| 免费高清视频大片| 久久久色成人| 免费av观看视频| 免费在线观看影片大全网站| 他把我摸到了高潮在线观看| 啪啪无遮挡十八禁网站| 成人亚洲精品av一区二区| 757午夜福利合集在线观看| avwww免费| 国产精品久久视频播放| 日韩欧美在线二视频| 成人特级黄色片久久久久久久| 亚洲国产欧洲综合997久久,| 欧美乱妇无乱码| 99精品在免费线老司机午夜| 亚洲电影在线观看av| 成年女人永久免费观看视频| 简卡轻食公司| 国产伦精品一区二区三区四那| 日韩人妻高清精品专区| 国产精品一区二区三区四区久久| a在线观看视频网站| 99热这里只有是精品50| 亚洲国产高清在线一区二区三| 88av欧美| 国产欧美日韩一区二区精品| 美女被艹到高潮喷水动态| 午夜福利高清视频| 久久精品综合一区二区三区| 国产视频内射| 少妇熟女aⅴ在线视频| 久久精品国产自在天天线| 99riav亚洲国产免费| 欧美最黄视频在线播放免费| 亚洲精品亚洲一区二区| 熟女人妻精品中文字幕| 色5月婷婷丁香| 国产精品一区二区三区四区久久| 国产日本99.免费观看| 精品福利观看| 国产三级黄色录像| 久久精品国产99精品国产亚洲性色| 亚洲精品一卡2卡三卡4卡5卡| 久久久久久国产a免费观看| 亚洲美女搞黄在线观看 | 日韩欧美一区二区三区在线观看| 日韩欧美精品免费久久 | 亚洲乱码一区二区免费版| 亚洲av熟女| 亚洲av电影不卡..在线观看| 白带黄色成豆腐渣| 亚洲av美国av| 精品人妻熟女av久视频| 日韩欧美在线乱码| 3wmmmm亚洲av在线观看| 亚洲精品影视一区二区三区av| 国产又黄又爽又无遮挡在线| 成人国产一区最新在线观看| 舔av片在线| 一卡2卡三卡四卡精品乱码亚洲| 又紧又爽又黄一区二区| 国产黄a三级三级三级人| 亚洲一区高清亚洲精品| 小说图片视频综合网站| 亚洲av电影在线进入| 午夜久久久久精精品| 国产野战对白在线观看| a在线观看视频网站| 国产午夜精品论理片| 国产精品免费一区二区三区在线| 一级黄片播放器| 亚洲国产高清在线一区二区三| 成人特级av手机在线观看| 婷婷丁香在线五月| 999久久久精品免费观看国产| 免费大片18禁| 国产又黄又爽又无遮挡在线| 自拍偷自拍亚洲精品老妇| 一区二区三区免费毛片| 天堂av国产一区二区熟女人妻| 男女做爰动态图高潮gif福利片| 一级黄片播放器| 在线天堂最新版资源| 欧美性感艳星| 亚洲三级黄色毛片| 亚洲欧美日韩无卡精品| 日韩欧美精品免费久久 | 综合色av麻豆| 高清在线国产一区| av在线观看视频网站免费| 一边摸一边抽搐一进一小说| 午夜日韩欧美国产| 精品乱码久久久久久99久播| 亚洲久久久久久中文字幕| 97碰自拍视频| 在线观看午夜福利视频| 亚洲,欧美,日韩| 国产一区二区亚洲精品在线观看| 直男gayav资源| 深夜a级毛片| 亚洲人成电影免费在线| 波野结衣二区三区在线| 亚洲性夜色夜夜综合| 日日摸夜夜添夜夜添小说| 91av网一区二区| 午夜精品久久久久久毛片777| 日韩精品中文字幕看吧| 久久久久久国产a免费观看| 国产精品99久久久久久久久| 亚洲第一区二区三区不卡| 日本a在线网址| 中文字幕高清在线视频| 久久久久久久久大av| 亚洲一区高清亚洲精品| av在线蜜桃| 免费搜索国产男女视频| 欧美日韩国产亚洲二区| 国产成人福利小说| 老女人水多毛片| 男女下面进入的视频免费午夜| 国产视频内射| 国产高潮美女av| 午夜免费激情av| h日本视频在线播放| 99国产综合亚洲精品| 美女被艹到高潮喷水动态| 我的女老师完整版在线观看| x7x7x7水蜜桃| 日韩欧美 国产精品| 久久性视频一级片| 伦理电影大哥的女人| 18禁在线播放成人免费| 中文在线观看免费www的网站| eeuss影院久久| 成人鲁丝片一二三区免费| 波多野结衣高清无吗| 好男人在线观看高清免费视频| 国产精品亚洲一级av第二区| 狠狠狠狠99中文字幕| 免费一级毛片在线播放高清视频| 黄色女人牲交| 亚洲av不卡在线观看| 午夜影院日韩av| 精品无人区乱码1区二区| av在线蜜桃| 亚洲五月婷婷丁香| 亚洲av五月六月丁香网| 欧美性猛交╳xxx乱大交人| 丁香六月欧美| 波多野结衣高清作品| 最近最新中文字幕大全电影3| 亚洲美女视频黄频| 色综合欧美亚洲国产小说| 国产三级在线视频| 亚洲人成伊人成综合网2020| 国产高清视频在线播放一区| 欧美乱妇无乱码| 国产精品亚洲一级av第二区| 免费搜索国产男女视频| 在线十欧美十亚洲十日本专区| 美女xxoo啪啪120秒动态图 | 国产精品,欧美在线| 精品久久久久久成人av| 日韩成人在线观看一区二区三区| 中文字幕精品亚洲无线码一区| 亚洲精品在线美女| 十八禁人妻一区二区| 国产精品98久久久久久宅男小说| 成人午夜高清在线视频| 老司机深夜福利视频在线观看| 嫩草影院精品99| 看片在线看免费视频| 哪里可以看免费的av片| 国产在线精品亚洲第一网站| 一a级毛片在线观看| 色综合欧美亚洲国产小说| 亚洲美女黄片视频| x7x7x7水蜜桃| 国产成人啪精品午夜网站| 国产不卡一卡二| 特大巨黑吊av在线直播| 午夜福利18| 夜夜夜夜夜久久久久| 亚洲真实伦在线观看| 国产成人av教育| 97碰自拍视频| 热99re8久久精品国产| 成人av一区二区三区在线看| 好男人在线观看高清免费视频| 亚洲电影在线观看av| 18美女黄网站色大片免费观看| 黄色配什么色好看| 99国产精品一区二区蜜桃av| 久久精品国产亚洲av涩爱 | 可以在线观看的亚洲视频| 91久久精品电影网| 午夜精品久久久久久毛片777| 午夜福利在线观看免费完整高清在 | 久久国产精品人妻蜜桃| 亚洲成人精品中文字幕电影| 麻豆av噜噜一区二区三区| av天堂在线播放| 欧美色欧美亚洲另类二区| 国产高清视频在线播放一区| 成人av在线播放网站| 欧美一区二区精品小视频在线| 能在线免费观看的黄片| 一区二区三区高清视频在线| 亚洲成av人片免费观看| 亚洲中文字幕一区二区三区有码在线看| 午夜福利在线观看吧| 亚洲专区国产一区二区| 中国美女看黄片| 国产视频一区二区在线看| 日本 欧美在线| 黄色日韩在线| 97热精品久久久久久| 国产视频一区二区在线看| 性色avwww在线观看| 又紧又爽又黄一区二区| 男插女下体视频免费在线播放| 在线天堂最新版资源| 亚洲五月婷婷丁香| 亚洲一区二区三区色噜噜| 1024手机看黄色片| av在线蜜桃| 免费搜索国产男女视频| 18+在线观看网站| 欧美性猛交╳xxx乱大交人| 丁香六月欧美| 精品一区二区免费观看| 首页视频小说图片口味搜索| 欧美bdsm另类| 久久伊人香网站| 精品一区二区免费观看| 男女那种视频在线观看| 亚洲精品色激情综合| 99精品在免费线老司机午夜| 丰满乱子伦码专区| 国产伦一二天堂av在线观看| 嫩草影院新地址| 国产精品98久久久久久宅男小说| 国模一区二区三区四区视频| 变态另类成人亚洲欧美熟女| 亚洲欧美清纯卡通| 亚洲国产高清在线一区二区三| 久久午夜亚洲精品久久| 成人国产一区最新在线观看| 国产成人欧美在线观看| 久久精品久久久久久噜噜老黄 | 亚洲七黄色美女视频| 亚洲色图av天堂| 97热精品久久久久久| 亚洲电影在线观看av| 国产精品美女特级片免费视频播放器| 国产三级中文精品| 首页视频小说图片口味搜索| 搡老岳熟女国产| 丰满人妻熟妇乱又伦精品不卡| 亚洲18禁久久av| 欧美高清成人免费视频www| 欧美一区二区亚洲| 欧美激情国产日韩精品一区| 丝袜美腿在线中文| 在线天堂最新版资源| 欧美日韩福利视频一区二区| 国产野战对白在线观看| 欧美极品一区二区三区四区| 国产精品免费一区二区三区在线| 欧美最黄视频在线播放免费| 黄色视频,在线免费观看| 又黄又爽又免费观看的视频| 国产精品久久视频播放| 亚洲一区二区三区色噜噜| 很黄的视频免费| 午夜福利欧美成人| 国产 一区 欧美 日韩| 国产av不卡久久| 亚洲,欧美,日韩| 一个人观看的视频www高清免费观看| 国产视频一区二区在线看| 亚洲人成电影免费在线| 亚洲美女黄片视频|