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

    半干旱黃土區(qū)坡面尺度檸條生長(zhǎng)狀況及影響要素分析

    2018-01-04 05:57:21王子婷蔡國(guó)軍莫保儒柴春山戚建莉張洋東
    生態(tài)學(xué)報(bào) 2017年23期
    關(guān)鍵詞:坡位檸條冠幅

    王子婷,楊 磊,蔡國(guó)軍,莫保儒,柴春山,戚建莉,張洋東

    1 甘肅省林業(yè)科學(xué)研究院,蘭州 730020 2 中國(guó)科學(xué)院生態(tài)環(huán)境研究中心城市與區(qū)域生態(tài)國(guó)家重點(diǎn)實(shí)驗(yàn)室,北京 100085

    半干旱黃土區(qū)坡面尺度檸條生長(zhǎng)狀況及影響要素分析

    王子婷1,楊 磊2,*,蔡國(guó)軍1,莫保儒1,柴春山1,戚建莉1,張洋東1

    1 甘肅省林業(yè)科學(xué)研究院,蘭州 730020 2 中國(guó)科學(xué)院生態(tài)環(huán)境研究中心城市與區(qū)域生態(tài)國(guó)家重點(diǎn)實(shí)驗(yàn)室,北京 100085

    以半干旱黃土丘陵區(qū)典型小流域坡面大規(guī)模人工種植檸條林為例,基于坡面不同部位檸條生長(zhǎng)狀況和生境條件調(diào)查,定量分析了地形變化、土壤水分及灌木密度對(duì)檸條生長(zhǎng)的直接、間接影響及其貢獻(xiàn)率。結(jié)果表明:(1)東坡大株檸條生長(zhǎng)明顯好于南坡,下坡位檸條生長(zhǎng)狀況略好于中上坡位,其他各坡位之間檸條生長(zhǎng)狀況差異較?。?2)大株檸條生長(zhǎng)與淺層土壤水分有正相關(guān)關(guān)系,而與灌木密度和深層土壤水分則呈負(fù)相關(guān)關(guān)系;大株檸條灌木高度、灌木縱截面積和冠幅體積對(duì)淺層土壤水分的響應(yīng)敏感,冠幅長(zhǎng)度對(duì)坡向和坡位的響應(yīng)較為敏感,冠幅寬度對(duì)灌木密度的響應(yīng)較為敏感;(3)地形和土壤水分變化解釋了59.9%的大株檸條生長(zhǎng)變異,其中坡向、坡位和淺層土壤水分是影響大株檸條生長(zhǎng)的主導(dǎo)環(huán)境因子,它們分別解釋了21.1%、16.0%和13.1%的檸條生長(zhǎng)變化。研究認(rèn)為半干旱黃土區(qū)人工植被恢復(fù)既要重視空間布局,也要在后期實(shí)施必要的管理措施以維持人工林地的穩(wěn)定性。

    黃土高原;坡面尺度;坡位;人工檸條林;土壤水分

    半干旱黃土高原因受區(qū)域自然條件和人類活動(dòng)的雙重影響,生態(tài)環(huán)境脆弱、土壤侵蝕嚴(yán)重,一直是我國(guó)植被恢復(fù)和生態(tài)建設(shè)的重點(diǎn)區(qū)域之一[1-3],人工林草種植已經(jīng)成為這一地區(qū)植被恢復(fù)的主要措施[4-6],其中檸條是黃土高原地區(qū)人工種植最為廣泛的灌木樹種之一[7-8]。研究認(rèn)為,黃土高原人工植被恢復(fù)需要遵循植被地帶性規(guī)律[3-5]和水資源承載力[9-10]。黃土丘陵區(qū)地形破碎、溝壑縱橫,坡面水熱組合條件變化較大[11-12],從而導(dǎo)致人工植被在空間配置中還需充分考慮非地帶性因子[3,5],如黃土的質(zhì)地組成、坡向、坡度、坡位、坡面地形和水分條件的差異,因此在坡面尺度開展人工植被恢復(fù)須根據(jù)對(duì)應(yīng)的生境條件進(jìn)行科學(xué)配置,否則將不利于人工植被的穩(wěn)定以及水土流失控制等生態(tài)功能的有效發(fā)揮[13]。關(guān)于黃土高原坡面尺度植被恢復(fù)的研究,不同的研究者選取了不同的植被類型以期探尋植被恢復(fù)與坡面地形環(huán)境的響應(yīng)關(guān)系,其中王晶等[14]在陜北黃土區(qū)封禁流域坡面微地形內(nèi)以經(jīng)過(guò)自然恢復(fù)的植物群落為研究對(duì)象,分析其物種組成、數(shù)量特征及多樣性與環(huán)境因子的關(guān)系;楊士梭等[15]在黃土丘陵區(qū)以自然植被為對(duì)象,分析植物各功能性狀值在科屬和群落水平上與環(huán)境因子的響應(yīng)關(guān)系;胡嬋娟等[16]在丘陵溝壑區(qū)小流域內(nèi),不同坡面上對(duì)不同退耕還林植被配置模式下植物物種組成與土壤環(huán)境的變化關(guān)系進(jìn)行研究;胡相明等[17]在黃土丘陵區(qū)從坡面尺度以自然封育形成的天然草地為研究對(duì)象,對(duì)其結(jié)構(gòu)與地形因素和土壤水分的關(guān)系進(jìn)行了研究。然而,關(guān)于大規(guī)模植被恢復(fù)以后坡面尺度成熟檸條林生長(zhǎng)狀況及其影響因素的認(rèn)識(shí)還很有限,其與環(huán)境因子的定量表達(dá)關(guān)系還有待進(jìn)一步深入探討。本文以黃土高原廣泛分布的檸條為例,通過(guò)對(duì)半干旱黃土區(qū)坡面尺度不同部位檸條生長(zhǎng)狀況的分析與比較,探索坡面尺度檸條生長(zhǎng)規(guī)律,確定地形、種植密度和土壤水分變化對(duì)檸條生長(zhǎng)的直接、間接影響及其貢獻(xiàn)率,以期為半干旱黃土丘陵區(qū)植被恢復(fù)模式及空間配置提供科學(xué)依據(jù)。

    1 研究區(qū)概況

    研究區(qū)位于甘肅省定西市龍灘小流域(104°27′—104°32′E,35°43′—35°46′N),屬典型半干旱黃土丘陵溝壑區(qū),流域面積16km2,平均海拔1900m。研究區(qū)年平均氣溫6.8℃,1月份平均氣溫-7.9℃,極端最高溫38.5℃,平均無(wú)霜期152d,平均日照時(shí)數(shù)2052h。多年平均降水量為386mm,降雨主要集中在7—9月份,潛在蒸發(fā)量1439mm。年平均相對(duì)濕度72%,干燥度1.9。流域內(nèi)土壤以黃綿土為主,有機(jī)質(zhì)含量低,土壤貧瘠且易侵蝕。天然植被以多年生草本為主,目前主要草本種類長(zhǎng)芒草(Stipabungeana)、賴草(Leymussecalinus)、阿爾泰狗娃花(Heteropappusaltaicus)、紫花苜蓿(Medicagosativa),人工栽植的灌喬樹種有檸條(Caraganakorshinskii)、側(cè)柏(Platycladusorientalis)、油松(Pinustabuliformis)等,其中人工檸條林地面積較大且分布較廣,林下草本物種主要有長(zhǎng)芒草(S.bungeana)、灌木亞菊(Ajaniafruticulosa)、阿爾泰狗娃花(H.altaicus)及駱駝蓬(Peganumharmala)等。

    2 研究方法

    2.1 監(jiān)測(cè)樣點(diǎn)布設(shè)

    在研究區(qū)內(nèi)依據(jù)坡向不同選取5個(gè)樣帶,其中東坡取兩個(gè)樣帶分別命名為E1和E2,南坡選取3個(gè)樣帶分別命名為S1、S2和S3;每條樣帶自上而下依次選取上坡位、中上坡位、中坡位、中下坡位及下坡位5個(gè)坡位作為監(jiān)測(cè)樣點(diǎn)。坡面5條樣帶內(nèi)的檸條以水平階整地的造林方式栽植于1984年。

    2.2 數(shù)據(jù)采集與分析

    于生長(zhǎng)季(8月份)在各監(jiān)測(cè)樣點(diǎn),對(duì)檸條林測(cè)定20m×20m樣方內(nèi)每株檸條的灌木高度、冠幅長(zhǎng)度和冠幅寬度,用土鉆采集0—600cm深度土壤樣品,每隔20cm采集1個(gè)土樣,利用烘干法測(cè)定土壤水分含量。有研究表明,黃土高原地區(qū)0—200cm深度土壤水分受降水影響相對(duì)較大,而200cm以下土壤水分受降水影響較小[18-20]。因而根據(jù)土壤水分對(duì)降雨的響應(yīng)規(guī)律,本文選取0—200cm深土壤含水量作為淺層土壤水分,200—600cm深土壤含水量作為深層土壤水分,東坡和南坡不同坡位淺層土壤水分和深層土壤水分含量見表1。利用南方S740手持GPS記錄各監(jiān)測(cè)樣點(diǎn)的經(jīng)緯度及海拔信息,利用羅盤測(cè)定坡度和坡向。

    表1 東坡和南坡不同坡位淺層土壤水分和深層土壤水分

    E1:東1坡,East slope 1;E2:東2坡,East slope 2;S1:南1坡,South slope 1;S2:南2坡,South slope 2;S3:南3坡,South slope 3

    數(shù)據(jù)分析時(shí)為了避免樣方內(nèi)每株檸條樣本因大小株差異而影響顯著性及相關(guān)性分析結(jié)果,本文在分析統(tǒng)計(jì)時(shí)首先對(duì)樣方內(nèi)的檸條樣本按株高將其分為小株(0100cm)[21],在此基礎(chǔ)上分別統(tǒng)計(jì)代表檸條生長(zhǎng)的形態(tài)特征值并參考相關(guān)研究計(jì)算特征值來(lái)近似代表每株檸條地上生物量[21-23],包括檸條灌木高度(H)、冠幅長(zhǎng)度(CL)、冠幅寬度(CW)、冠幅直徑(C)、冠幅橫截面積(A)、灌木縱截面積(CH)及冠幅體積(V),其中:冠幅直徑C=(CL+CW)/2,冠幅橫截面積A=πC2/4,灌木縱截面積CH=C×H,冠幅體積V=A×H。

    利用相似性分析(Analysis of similarities,ANOSIM)檢驗(yàn)坡向和坡位對(duì)大株檸條(中株和小株檸條數(shù)據(jù)較少不做分析)生長(zhǎng)的影響,ANOSIM分析使用PRIMER 5.0[24]。運(yùn)用單因素方差分析比較5條樣帶不同坡位對(duì)大株和中株檸條7個(gè)生長(zhǎng)指標(biāo)的影響,采用Pearson相關(guān)分析不同坡向和不同坡位與坡向大株和中株檸條灌木密度、淺層土壤水分和深層土壤水分與檸條生長(zhǎng)的相關(guān)關(guān)系,利用通徑分析確定灌木密度、淺層土壤水分及深層土壤水分對(duì)檸條生長(zhǎng)的直接和間接影響,統(tǒng)計(jì)分析使用SPS S21.0。采用RDA和pRDA分析地形因子(坡向、坡位及坡度)和土壤水分(淺層土壤水分及深層土壤水分)對(duì)大株和中株檸條生長(zhǎng)狀況(灌木高度、冠幅寬度、冠幅長(zhǎng)度、冠幅直徑、冠幅橫截面積、灌木縱截面積和冠幅體積)的影響及貢獻(xiàn)率,排序分析使用CANOCO 4.5[25]。

    3 結(jié)果

    3.1 坡向和坡位對(duì)檸條生長(zhǎng)的影響

    檸條生長(zhǎng)狀況的相似性分析(ANOSIM)結(jié)果表明,坡向?qū)帡l生長(zhǎng)有顯著影響(GlobalR=0.267,P=0.014),東坡檸條的生長(zhǎng)好于南坡,東坡檸條的高度、冠幅長(zhǎng)度及寬度均大于南坡。5條樣帶綜合分析的結(jié)果表明,坡位對(duì)檸條生長(zhǎng)的影響較小(GlobalR=-0.054、P=0.676)。不同樣帶內(nèi)檸條生長(zhǎng)對(duì)坡位的響應(yīng)可能存在一定差異,進(jìn)而對(duì)5條樣帶內(nèi)不同坡位檸條生長(zhǎng)狀況進(jìn)行方差分析。結(jié)果表明,就大株檸條而言,坡位對(duì)大株檸條的生長(zhǎng)有顯著影響,但這種影響隨坡向不同又存在一定變化(表2)。東坡兩條樣帶內(nèi)下坡位檸條的平均高度要高于上坡位,其中E1樣帶內(nèi)上坡和中上坡與下坡位檸條高度相差較小,E2坡檸條的高度則自上而下逐漸增加(表2)。東坡檸條冠幅長(zhǎng)度對(duì)坡位變化的響應(yīng)較冠幅寬度明顯,下坡檸條的冠幅長(zhǎng)度顯著高于上坡,而上坡和下坡檸條冠幅寬度相差較小(表2)。整體而言,在表征檸條生長(zhǎng)的7個(gè)生長(zhǎng)指標(biāo)中,冠幅橫截面積、灌木縱截面積和冠幅體積隨坡位變化相比植株高度和冠幅明顯,下坡位檸條的冠幅橫截面積、灌木縱截面積和冠幅體積普遍高于中上坡位(表2)。南坡S1和S3樣帶灌木高度及冠幅長(zhǎng)度均表現(xiàn)為下坡位高于上坡位及中上坡位,檸條的冠幅橫截面積、灌木縱截面積及冠幅體積變化表現(xiàn)更為明顯(表2)。S2樣帶與S1和S3樣帶不同,S2樣帶下坡位檸條的高度、冠幅長(zhǎng)度、冠幅橫截面積、灌木縱截面積及冠幅體積均低于中上坡位(表2)??傊痛笾隀帡l而言,下坡位檸條生長(zhǎng)要好于其他坡位(S2除外),而上坡和中下坡之間檸條生長(zhǎng)差異較小。

    表2 東坡和南坡不同坡位大株檸條生長(zhǎng)狀況比較

    表中不同小寫字母表示同一坡面上不同坡位之間差異顯著;H:灌木高度,Height;CL:冠幅長(zhǎng)度,Crown length;CW:冠幅寬度,Crown width;C:冠幅直徑,Crown diameter;A:冠幅橫截面積,Crown cross sectional area;CH:灌木縱截面積,Crown diameter×height;V:冠幅體積,Canopy volume

    5條樣帶內(nèi)中株檸條的密度較低,其生長(zhǎng)狀況與大株檸條相比對(duì)坡面變化響應(yīng)較小。E1、E2和S1樣帶內(nèi)不同坡位檸條生長(zhǎng)差異較小(P>0.05),S2坡冠幅長(zhǎng)度在坡位間的差異極顯著(P<0.001),冠幅直徑(P=0.002)和灌木縱截面積(P=0.003)在坡位間也存在顯著差異,且上述3個(gè)生長(zhǎng)指標(biāo)均表現(xiàn)為中上坡位好于下坡位,而其他生長(zhǎng)指標(biāo)在坡位間無(wú)明顯差異;S3坡僅灌木高度在坡位間存在極顯著差異(P<0.001),其他各生長(zhǎng)指標(biāo)均無(wú)明顯差異。

    3.2 土壤水分及灌木密度與檸條生長(zhǎng)的關(guān)系

    檸條生長(zhǎng)與灌木密度的相關(guān)分析結(jié)果表明,大株檸條灌木密度與7個(gè)生長(zhǎng)指標(biāo)均表現(xiàn)為負(fù)相關(guān)關(guān)系,其中東坡和南坡灌木密度與各生長(zhǎng)指標(biāo)的相關(guān)性均較小(表3)。檸條生長(zhǎng)與淺層土壤水分的相關(guān)分析結(jié)果表明,淺層土壤水分與大株檸條各指標(biāo)均表現(xiàn)為正相關(guān)關(guān)系,其中淺層土壤水分與灌木高度和灌木縱截面積呈顯著正相關(guān)。就不同坡向而言,東坡淺層土壤水分與檸條各生長(zhǎng)指標(biāo)均呈正相關(guān)關(guān)系,而南坡淺層土壤水分與各生長(zhǎng)指標(biāo)均呈負(fù)相關(guān)關(guān)系,其中南坡淺層土壤水分與冠幅寬度顯著負(fù)相關(guān)。檸條生長(zhǎng)與深層土壤水分的相關(guān)關(guān)系表明,雖然深層土壤水分與檸條各生長(zhǎng)指標(biāo)的相關(guān)性較小,但不同坡向之間的表現(xiàn)存在差異。東坡深層土壤水分與檸條各生長(zhǎng)指標(biāo)表現(xiàn)為正相關(guān)關(guān)系,其中東坡深層土壤水分與冠幅寬度呈顯著正相關(guān),而南坡深層土壤水分與各生長(zhǎng)指標(biāo)呈負(fù)相關(guān)關(guān)系。中株檸條與大株檸條不同,其密度及土壤水分與中株檸條各生長(zhǎng)指標(biāo)的相關(guān)性不顯著(表3)。

    進(jìn)一步對(duì)大株檸條各生長(zhǎng)指標(biāo)與灌木密度、淺層土壤水分和深層土壤水分進(jìn)行通徑分析。結(jié)果表明,淺層土壤水分對(duì)大株檸條灌木高度、灌木縱截面積及冠幅體積均有直接的正影響(通徑系數(shù)分別為0.593、0.401和0.312),大株檸條密度對(duì)冠幅寬度有直接的負(fù)影響(通徑系數(shù)為-0.323),而深層土壤水分及大株檸條密度對(duì)灌木高度、灌木縱截面積和冠幅體積產(chǎn)生間接負(fù)影響(表4)。

    表3 淺層土壤水分、深層土壤水分及灌木密度與檸條生長(zhǎng)指標(biāo)的相關(guān)分析

    *P< 0.05, **P< 0.01

    表4大株檸條密度、淺層土壤水分及深層土壤水分對(duì)檸條生長(zhǎng)的直接影響和間接影響

    Table4Directlyandindirectlyeffectofshallow&deepsoilmoisturecontentandplantdensitytogrowthstatusasbig-sizeCaraganakorshinskiishrubs

    生長(zhǎng)指標(biāo)Growthindices變量Variate直接效應(yīng)Directlyeffect間接效應(yīng)Indirectlyeffect通過(guò)x1通過(guò)x2通過(guò)x3總效應(yīng)Totaleffect灌木高度大株檸條密度(x1)0.018-0.3560.011-0.328Height/m淺層土壤水分(x2)0.593-0.011-0.0140.569深層土壤水分(x3)0.0250.007-0.328-0.296R2=0.325*冠幅長(zhǎng)度大株檸條密度(x1)-0.001-0.083-0.013-0.097Crownlength/m淺層土壤水分(x2)0.1390.0010.0170.156深層土壤水分(x3)-0.0310.000-0.077-0.108R2=0.025冠幅寬度大株檸條密度(x1)-0.3230.0250.118-0.180Crownwidth/m淺層土壤水分(x2)-0.0420.194-0.153-0.001深層土壤水分(x3)0.277-0.1370.0230.163R2=0.104冠幅直徑大株檸條密度(x1)-0.165-0.0350.040-0.159Crowndiameter/m淺層土壤水分(x2)0.0580.099-0.0520.105深層土壤水分(x3)0.094-0.070-0.032-0.008R2=0.032冠幅橫截面積大株檸條密度(x1)-0.125-0.0550.049-0.130Crowncroectionalarea/m2淺層土壤水分(x2)0.0910.075-0.0640.102深層土壤水分(x3)0.116-0.053-0.0500.013R2=0.027灌木縱截面積大株檸條密度(x1)-0.063-0.2410.028-0.276Crowndiameter×height/m2淺層土壤水分(x2)0.4010.038-0.0360.403深層土壤水分(x3)0.065-0.027-0.222-0.183R2=0.167冠幅體積大株檸條密度(x1)-0.081-0.1870.026-0.242Canopyvolume/m3淺層土壤水分(x2)0.3120.049-0.0330.327深層土壤水分(x3)0.060-0.034-0.172-0.146R2=0.112

    *P<0.05

    3.3 坡面檸條生長(zhǎng)與環(huán)境因子的定量分析

    利用RDA和pRDA分析5條樣帶不同坡位大株和中株檸條生長(zhǎng)狀況與5個(gè)環(huán)境因子(坡向、坡位、坡度、淺層土壤水分及深層土壤水分)的關(guān)系,確定影響大株和中株檸條生長(zhǎng)的關(guān)鍵環(huán)境因子。大株檸條與環(huán)境因子的排序結(jié)果表明,第1軸解釋了58.0%的檸條生長(zhǎng)變化(Eigenvalue=0.580,F(xiàn)=26.26,P=0.018),所有的排序軸解釋了59.9%的檸條生長(zhǎng)變化(Eigenvalue=0.599,F(xiàn)=5.68,P=0.015),排序結(jié)果較好的反映了地形和土壤水分變化對(duì)大株檸條生長(zhǎng)的影響。pRDA分析結(jié)果表明,坡向?qū)Υ笾隀帡l的生長(zhǎng)影響最大(解釋了21.1%的變異),其次是坡位(16.0%)、淺層土壤水分(13.1%)和深層土壤水分(9.6%),而坡度對(duì)大株檸條生長(zhǎng)的影響解釋量不足0.1%(圖1)。中株檸條與環(huán)境因子的排序結(jié)果表明,第1軸解釋了29.5%的檸條生長(zhǎng)變化(Eigenvalue=0.295,F(xiàn)=5.87,P=0.339),所有的排序軸解釋了30.6%的檸條生長(zhǎng)變化(Eigenvalue=0.306,F(xiàn)=1.24,P=0.340),排序結(jié)果表明地形和土壤水分變化對(duì)中株檸條生長(zhǎng)的影響較小。pRDA分析結(jié)果表明,坡向?qū)χ兄隀帡l生長(zhǎng)影響較大(解釋了18.9%的變異),深層土壤水分、淺層土壤水分、坡度和坡位對(duì)中株檸條生長(zhǎng)影響較小(解釋量分別為:6%、3.2%、2%和0.4%)(圖1)。

    圖 1 大株和中株檸條生長(zhǎng)狀況與環(huán)境因子的RDA排序圖Fig.1 RDA ordination figure of Caragana korshinskii growth indices and environmental factors at big-sized plants and medium-sized plantsSA:坡向,slope aspect;SP:坡位,slope position;SG:坡度,slope gradient;SMC1:淺層土壤水分,shallow soil moisture content;SMC2:深層土壤水分,deep soil moisture content;1:上坡位,upper slope position;2:中上坡位,upper-middle slope position;3:中坡位,middle slope position;4:中下坡,middle-low slope position;5:下坡位,low slope position

    4 討論

    半干旱黃土區(qū)人工植被的生長(zhǎng)受土壤水分、土壤養(yǎng)分及光熱組合等因素的共同制約[12,26],其中土壤水分是該區(qū)植被生長(zhǎng)的主要限制因子[27-28]。坡面地形變化塑造的小生境決定了土壤水分條件,它與土壤養(yǎng)分相互作用共同影響植被的生長(zhǎng)狀況,從而改變坡面植被生長(zhǎng)的空間分布格局[29-30]。本文的分析結(jié)果表明坡向?qū)帡l植被生長(zhǎng)的影響較大,區(qū)域內(nèi)東坡檸條的生長(zhǎng)明顯好于南坡,這主要是因?yàn)檠芯繀^(qū)東坡的光照不如南坡的強(qiáng)烈,致使土壤蒸發(fā)相對(duì)較小從而使得淺層土壤水分較高(表1),檸條生長(zhǎng)表現(xiàn)為東坡好于南坡。研究區(qū)坡面的坡度大致都在29°—32°之間,一般認(rèn)為坡度越大,土壤容積含水率越低,同一坡面自坡頂?shù)狡碌祝寥篮食试黾于厔?shì)[31],但是坡面工程措施的改造及植被恢復(fù)生長(zhǎng)影響了土壤水分及養(yǎng)分的再分配[32-33]。研究區(qū)檸條種植初期通過(guò)水平階整地在坡面栽植檸條,試圖削弱自然坡面徑流從而實(shí)現(xiàn)降水匯集,為植被生長(zhǎng)提供適宜條件。本文的分析結(jié)果表明,相對(duì)于坡向而言,坡位對(duì)檸條的生長(zhǎng)影響不及坡向?qū)帡l生長(zhǎng)的影響大,樣帶內(nèi)僅下坡位檸條生長(zhǎng)狀況略好于中上坡位(除S2樣帶),而其他各坡位之間檸條生長(zhǎng)狀況差異較小,表明水平階整地在一定程度上削弱了坡面地形變化對(duì)檸條生長(zhǎng)的限制。該結(jié)果與趙艷云等[34]在寧夏黃土區(qū)的研究結(jié)果相似,該研究發(fā)現(xiàn)在工程措施改造后的坡面,坡位對(duì)17年生檸條的分枝數(shù)影響較小,但對(duì)株高、根深、生物量均有顯著影響,且中坡檸條的株高、根深、生物量明顯高于上坡。程杰等[35]在寧夏黃土區(qū)對(duì)不同栽植年限檸條生長(zhǎng)的研究還發(fā)現(xiàn),隨著檸條栽植年限的增加坡位對(duì)檸條生長(zhǎng)的影響存在明顯變化,幼齡期(1—7年生)和老齡期(16—23年生)檸條株高表現(xiàn)為下坡>中坡>上坡,中齡期(8—15年生)表現(xiàn)為上坡>中坡>下坡。本研究也發(fā)現(xiàn),下坡位大株檸條的灌木高度和冠幅長(zhǎng)度高于中上坡位,而坡位對(duì)中株和小株檸條生長(zhǎng)的影響較小。

    坡面地形變化會(huì)通過(guò)改變土壤水分含量及養(yǎng)分條件對(duì)檸條生長(zhǎng)產(chǎn)生影響,而在人工植被恢復(fù)的過(guò)程中隨著生長(zhǎng)年限的增加,種植密度對(duì)植被生長(zhǎng)的制約會(huì)逐漸增強(qiáng)[36]。研究發(fā)現(xiàn)東坡檸條的密度低于南坡,東坡大株檸條的7個(gè)生長(zhǎng)指標(biāo)均高于南坡,而中株及小株檸條的生長(zhǎng)受密度影響較小。以上結(jié)果與趙龍等[37]在寧夏黃土區(qū)的研究結(jié)果相近,檸條的成林密度顯著影響其生長(zhǎng)過(guò)程,低密度檸條林對(duì)其生長(zhǎng)速率的制約強(qiáng)度小于高密度檸條林。此外,研究還發(fā)現(xiàn)大株檸條灌木高度、冠幅寬度和冠幅長(zhǎng)度與灌木密度呈負(fù)相關(guān),其中冠幅寬度與密度的相關(guān)性要大于冠幅長(zhǎng)度與密度的相關(guān)性,這在東坡表現(xiàn)尤為明顯;中株檸條灌木高度、冠幅寬度和冠幅長(zhǎng)度與灌木密度亦呈負(fù)相關(guān),但冠幅長(zhǎng)度與密度的相關(guān)性要大于與冠幅寬度與密度的相關(guān)性,這在南坡表現(xiàn)的尤為明顯。結(jié)果表明,對(duì)大株檸條而言灌木密度主要限制灌木高度和冠幅寬度的生長(zhǎng),對(duì)中株檸條而言灌木密度主要限制冠幅長(zhǎng)度的生長(zhǎng)。這可能與研究區(qū)大規(guī)模人工種植檸條林不同生長(zhǎng)階段檸條的生長(zhǎng)策略有關(guān),研究區(qū)坡面檸條為水平階整地造林,各階地之間坡面距離大于水平階內(nèi)檸條株間距,坡面地形變化影響了局部土壤水分狀況進(jìn)而決定了植被承載力,檸條密度和個(gè)體大小又影響了植物種間對(duì)水分及養(yǎng)分競(jìng)爭(zhēng)的強(qiáng)弱,二者相互作用決定了區(qū)域檸條的生長(zhǎng)狀況[34,36]。最終在成熟檸條林中表現(xiàn)為隨著檸條生長(zhǎng)年限的增加,其水平方向的生長(zhǎng)會(huì)受到相鄰植株的制約,灌木密度主要通過(guò)限制檸條水平(沿等高線方向)冠幅即冠幅寬度的生長(zhǎng)進(jìn)而抑制檸條植株的生長(zhǎng)。

    坡面地形變化通過(guò)改變土壤水分狀況,進(jìn)而與灌木密度相互作用來(lái)共同調(diào)節(jié)坡面檸條生長(zhǎng),同時(shí)檸條生長(zhǎng)也改變了土壤水分的存儲(chǔ)。有研究發(fā)現(xiàn)在坡面尺度上植被的分布及其生物量與不同土層深度土壤水分含量存在耦合關(guān)系[38],本研究發(fā)現(xiàn)淺層土壤水分與檸條各生長(zhǎng)指標(biāo)呈現(xiàn)正相關(guān)關(guān)系,深層土壤水分與檸條生長(zhǎng)呈現(xiàn)負(fù)相關(guān)關(guān)系,該研究結(jié)果與郭忠升[27]和莫保儒等[39]的研究結(jié)論相近。同時(shí)在半干旱區(qū)植被類型的不同對(duì)土壤水分的消耗明顯不同且存在顯著的空間分異,而同一植被類型,生長(zhǎng)年限及栽植密度的不同對(duì)土壤水分的消耗亦有不同[18,36,40-41]。研究區(qū)坡面地帶性植被是天然草本,人工植被超過(guò)了該區(qū)土壤水分的承載力,故檸條林地后期管理更需要考慮該區(qū)水分的承載力[9-10,42],選擇適宜的管理模式。研究發(fā)現(xiàn),淺層土壤水分直接影響坡面檸條生長(zhǎng),而深層土壤水分和灌木密度則與檸條生長(zhǎng)呈負(fù)相關(guān)關(guān)系。此外,研究還發(fā)現(xiàn)灌木密度低長(zhǎng)勢(shì)好的E2樣帶深層土壤水分低于淺層土壤水分,而其他樣帶檸條密度均高于E2樣帶且深層土壤水分高于淺層土壤水分,該結(jié)果與周萍等[43]在陜西紙坊溝流域的研究結(jié)果相似。以上結(jié)果揭示,檸條蓋度或生物量的變化強(qiáng)烈影響淺層和深層土壤水分的關(guān)系,人工植被對(duì)土壤水分的過(guò)度消耗導(dǎo)致深層土壤水分的匱缺,是引起林地退化的原因之一[44]。坡面檸條的生長(zhǎng),既受土壤水分的影響,又受灌木密度的調(diào)控,隨著檸條地上、地下生物量的增加對(duì)土壤水分抑制的程度也在增強(qiáng),土壤水分的匱缺會(huì)導(dǎo)致檸條林地退化。因而,半干旱黃土區(qū)進(jìn)行人工檸條植被恢復(fù)時(shí)既要結(jié)合立地條件選用不同的工程措施創(chuàng)造適宜植被生長(zhǎng)的小生境,也要依據(jù)土壤水分的植被承載力在植株生長(zhǎng)的不同階段調(diào)控植株密度和空間配置,以維持人工植被恢復(fù)的可持續(xù)性。

    5 結(jié)論

    研究表明坡向?qū)帡l生長(zhǎng)有顯著影響,東坡檸條的生長(zhǎng)好于南坡,而相對(duì)于坡向而言,坡位對(duì)檸條生長(zhǎng)的影響較小。淺層土壤水分對(duì)檸條生長(zhǎng)有直接影響,而灌木密度和深層土壤水分則與檸條生長(zhǎng)指標(biāo)存在負(fù)相關(guān)關(guān)系。大株檸條灌木高度、灌木縱截面積和冠幅體積對(duì)淺層土壤水分的響應(yīng)敏感,而冠幅長(zhǎng)度對(duì)坡向和坡位的響應(yīng)較為敏感,冠幅寬度對(duì)灌木密度的響應(yīng)較為敏感。地形變化和土壤水分與檸條生長(zhǎng)的RDA分析結(jié)果表明,地形和土壤水分變化解釋了59.9%的大株檸條生長(zhǎng)變異,其中坡向、坡位和淺層土壤水分是影響大株檸條生長(zhǎng)的主要環(huán)境因子,它們分別解釋了21.1%、16.0%和13.1%的檸條生長(zhǎng)變化;地形和土壤水分變化對(duì)中株檸條生長(zhǎng)的影響較小,它們僅解釋了30.6%的檸條生長(zhǎng)變異,其中坡向解釋了18.9%的檸條生長(zhǎng)變化,其他因子對(duì)檸條生長(zhǎng)影響均較小。坡面地形、土壤水分及灌木密度互饋互作調(diào)節(jié)了檸條生長(zhǎng)的空間分布格局。研究認(rèn)為半干旱黃土區(qū)坡面人工植被既要結(jié)合地形條件選擇適宜的工程措施,也需要在植被的不同恢復(fù)階段依據(jù)土壤水分植被承載力實(shí)施必要的管理措施(如平茬、間伐等),這將有利于維持坡面人工植被的穩(wěn)定,從而抑制因林地退化而造成的水土流失。

    [1] 李鳳民, 徐進(jìn)章, 孫國(guó)鈞. 半干旱黃土高原退化生態(tài)系統(tǒng)的修復(fù)與生態(tài)農(nóng)業(yè)發(fā)展. 生態(tài)學(xué)報(bào), 2003, 23(9): 1901-1909.

    [2] 楊榮金, 傅伯杰, 劉國(guó)華, 馬克明. 黃土丘陵溝壑區(qū)生態(tài)環(huán)境建設(shè)中的水問(wèn)題——以延河流域?yàn)槔? 環(huán)境科學(xué), 2004, 25(2): 37-42.

    [3] 山侖, 徐炳成. 黃土高原半干旱地區(qū)建設(shè)穩(wěn)定人工草地的探討. 草業(yè)學(xué)報(bào), 2009, 18(2): 1-2.

    [4] 王晗生. 黃土高原植被恢復(fù)策略回顧. 中國(guó)水土保持科學(xué), 2004, 2(1): 42-45.

    [5] 張金屯. 黃土高原植被恢復(fù)與建設(shè)的理論和技術(shù)問(wèn)題. 水土保持學(xué)報(bào), 2004, 18(5): 120-124.

    [6] 張文輝, 劉國(guó)彬. 黃土高原植被恢復(fù)與建設(shè)策略. 中國(guó)水土保持, 2009, 2009(1): 24-27.

    [7] 程積民, 杜峰, 萬(wàn)惠娥. 黃土高原半干旱區(qū)集流灌草立體配置與水分調(diào)控. 草地學(xué)報(bào), 2000, 8(3): 210-219.

    [8] 程積民, 萬(wàn)惠娥, 杜鋒. 黃土高原半干旱區(qū)退化灌草植被的恢復(fù)與重建. 林業(yè)科學(xué), 2001, 37(4): 50-57.

    [9] 郭忠升, 邵明安. 半干旱區(qū)人工林草地土壤旱化與土壤水分植被承載力. 生態(tài)學(xué)報(bào), 2003, 23(8):1640-1647.

    [10] 郭忠升. 黃土高原半干旱區(qū)水土保持植被恢復(fù)限度——以人工檸條林為例. 中國(guó)水土保持科學(xué), 2009, 7(4): 49-54.

    [11] 佘冬立, 邵明安, 薛亞鋒, 俞雙恩. 坡面土地利用格局變化的水土保持效應(yīng). 農(nóng)業(yè)工程學(xué)報(bào), 2011, 27(4): 22-27.

    [12] 李萍, 朱清科, 劉中奇, 趙薈, 鄺高明, 王晶. 半干旱黃土區(qū)地上生物量對(duì)立地因子的響應(yīng). 中國(guó)水土保持科學(xué), 2012, 10(2): 50-54.

    [13] 程積民, 萬(wàn)惠娥. 中國(guó)黃土高原植被建設(shè)與水土保持. 北京: 中國(guó)林業(yè)出版社, 2002: 140-300.

    [14] 王晶, 朱清科, 秦偉, 張宏芝, 云雷, 謝靜, 鄺高明. 陜北黃土區(qū)封禁流域坡面微地形植被特征分異. 應(yīng)用生態(tài)學(xué)報(bào), 2012, 23(3): 694-700.

    [15] 楊士梭, 溫仲明, 苗連朋, 戚德輝, 花東文. 黃土丘陵區(qū)植物功能性狀對(duì)微地形變化的響應(yīng). 應(yīng)用生態(tài)學(xué)報(bào), 2014, 25(12): 3413-3419.

    [16] 胡嬋娟, 劉國(guó)華, 陳利頂, 吳雅瓊. 黃土丘陵區(qū)坡面尺度上不同植被格局下植物群落和土壤性質(zhì)研究. 干旱區(qū)地理, 2012, 35(5): 787-794.

    [17] 胡相明, 程積民, 萬(wàn)惠娥, 趙艷云. 黃土丘陵區(qū)地形、土壤水分與草地的景觀格局. 生態(tài)學(xué)報(bào), 2006, 26(10): 3276-3285.

    [18] 程積民, 萬(wàn)惠娥, 王靜, 雍紹萍. 半干旱區(qū)檸條生長(zhǎng)與土壤水分消耗過(guò)程研究. 林業(yè)科學(xué), 2005, 41(2): 37-41.

    [19] Wang Z Q, Liu B Y, Zhang Y. Soil moisture of different vegetation types on the Loess Plateau. Journal of Geographical Sciences, 2009, 19(6): 707-718.

    [20] Wang Y Q, Shao M A, Shao H B. A preliminary investigation of the dynamic characteristics of dried soil layers on the Loess Plateau of China. Journal of Hydrology, 2010, 381(1): 9-17.

    [21] 李剛, 趙祥, 張賓賓, 劉碧榮, 董寬虎, 朱慧森. 不同株高的檸條生物量分配格局及其估測(cè)模型構(gòu)建. 草地學(xué)報(bào), 2014, 22(4): 769-775.

    [22] 李鋼鐵, 秦富倉(cāng), 賈守義, 董錦蘭, 吉日格勒, 王滿貴. 旱生灌木生物量預(yù)測(cè)模型的研究. 內(nèi)蒙古林學(xué)院學(xué)報(bào)(自然科學(xué)版), 1998, 20(2): 25-31.

    [23] Ferreira J N, Bustamante M, Garcia-Montiel D C, Caylor K K, Davidson E A. Spatial variation in vegetation structure coupled to plant available water determined by two-dimensional soil resistivity profiling in a Brazilian savanna. Oecologia, 2007, 153(2): 417-430.

    [24] Clarke K R, Warwick R M. Change in Marine Communities: An Approach to Statistical Analysis and Interpretation. 2nd ed. Plymouth: Plymouth Marine Laboratory, 2001: 2-15.

    [26] 保長(zhǎng)虎, 張文輝, 何景峰. 黃土高原丘陵溝壑區(qū)30年檸條人工種群動(dòng)態(tài)研究. 西北植物學(xué)報(bào), 2010, 30(8): 1681-1688.

    [27] 郭忠升. 半干旱區(qū)檸條林利用土壤水分深度和耗水量. 水土保持通報(bào), 2009, 29(5): 69-72.

    [28] Yang L, Wei W, Chen L D, Chen W L, Wang J L. Response of temporal variation of soil moisture to vegetation restoration in semi-arid Loess Plateau, China. Catena, 2014, 115(1): 123-133.

    [29] 孟慶華, 傅伯杰. 景觀格局與土壤養(yǎng)分流動(dòng). 水土保持學(xué)報(bào), 2000, 14(3): 116-121.

    [30] 衛(wèi)偉, 余韻, 賈福巖, 楊磊, 陳利頂. 微地形改造的生態(tài)環(huán)境效應(yīng)研究進(jìn)展. 生態(tài)學(xué)報(bào), 2013, 33(20): 6462-6469.

    [31] 李毅, 邵明安. 雨強(qiáng)對(duì)黃土坡面土壤水分入滲及再分布的影響. 應(yīng)用生態(tài)學(xué)報(bào), 2006, 17(12): 2271-2276.

    [32] 于洋, 衛(wèi)偉, 陳利頂, 馮天驕, 楊磊, 張涵丹. 黃土丘陵區(qū)坡面整地和植被耦合下的土壤水分特征. 生態(tài)學(xué)報(bào), 2016, 36(11): 3441-3449.

    [33] 馮天驕, 衛(wèi)偉, 陳利頂, 于洋, 楊磊, 張涵丹. 隴中黃土區(qū)坡面整地和植被類型對(duì)土壤化學(xué)性狀的影響. 生態(tài)學(xué)報(bào), 2016, 36(11): 3216-3225.

    [34] 趙艷云, 程積民, 王延平, 萬(wàn)惠娥, 胡相明. 半干旱區(qū)環(huán)境因子對(duì)檸條灌木林結(jié)構(gòu)的影響. 水土保持通報(bào), 2005, 25(3): 10-14.

    [35] 程杰, 王吉斌, 程積民, 羅宗寬. 黃土高原檸條錦雞兒灌木林生長(zhǎng)的時(shí)空變異特征. 林業(yè)科學(xué), 2013, 49(1): 14-20.

    [36] 張文文, 郭忠升, 寧婷, 白冬妹. 黃土丘陵半干旱區(qū)檸條林密度對(duì)土壤水分和檸條生長(zhǎng)的影響. 生態(tài)學(xué)報(bào), 2015, 35(3): 725-732.

    [37] 趙龍, 郭忠升, 郭滿才, 袁志發(fā), 王振鳳. 黃土丘陵半干旱區(qū)檸條林生長(zhǎng)過(guò)程研究. 水土保持通報(bào), 2012, 32(6): 192-195.

    [38] Meerveld H J T V, McDonnell J J. On the interrelations between topography, soil depth, soil moisture, transpiration rates and species distribution at the hillslope scale. Advances in Water Resources, 2006, 29(2): 293-310.

    [39] 莫保儒, 王子婷, 蔡國(guó)軍, 楊磊, 黨宏忠, 王多鋒, 薛睿. 半干旱黃土區(qū)成熟檸條林地剖面土壤水分環(huán)境及影響因子研究. 干旱區(qū)地理, 2014, 37(6): 1207-1215.

    [40] 徐學(xué)選, 劉江華, 高鵬, 王煒. 黃土丘陵區(qū)植被的土壤水文效應(yīng). 西北植物學(xué)報(bào), 2003, 23(8): 1347-1351.

    [41] Vivoni E R, Rinehart A J, Méndez-Barroso L A, Aragón C A, Bisht G, Cardenas M B, Engle E, Forman B A, Frisbee M D, Gutiérrez-Jurado H A, Hong S H, Mahmood T H, Tai K, Wyckoff R L. Vegetation controls on soil moisture distribution in the Valles Caldera, New Mexico, during the North American monsoon. Ecohydrology, 2008, 1(3): 225-238.

    [42] 郭忠升. 水資源緊缺地區(qū)土壤水分植被承載力論述. 林業(yè)科學(xué), 2011, 47(5): 140-144.

    [43] 周萍, 劉國(guó)彬, 侯喜祿. 黃土丘陵區(qū)侵蝕環(huán)境不同坡面及坡位土壤理化特征研究. 水土保持學(xué)報(bào), 2008, 22(1): 7-12.

    [44] 邵明安, 賈小旭, 王云強(qiáng), 朱元駿. 黃土高原土壤干層研究進(jìn)展與展望. 地球科學(xué)進(jìn)展, 2016, 31(1): 14-22.

    SpatialpatternsofCaraganakorshinskiigrowthonhillslopescaleandinfluencingfactorsinthesemi-aridLoessPlateau

    WANG Ziting1, YANG Lei2,*, CAI Guojun1, MO Baoru1, CHAI Chunshan1, QI Jianli1,ZHANG Yangdong1

    1GansuAcademyofForestryScience,Lanzhou730020,China2StateKeyLaboratoryofUrbanandRegionalEcology,ResearchCenterforEco-EnvironmentalSciences,ChineseAcademyofSciences,Beijing100085,China

    Topography and soil moisture are the key factors affecting the spatial growth patterns of artificialCaraganakorshinskiishrubs in the semi-arid loess region. However, the recognition of the growth status ofC.korshinskiishrubs affected by topographic variation, soil moisture content, and shrub density is limited on the hillslope scale. In this study, we selected a typical artificialC.korshinskiishrub on the hillslope as a research model, and the data on its growth status and habitat conditions under different (east- and south-facing) slope positions were collected. Furthermore, the relationships between the growth status ofC.korshinskiishrubs and the shrub density, topographic variation, and soil moisture content were quantitatively analyzed using the Pearson correlation, path analysis, and RDA. The following results were obtained in this study. (1)The growth status ofC.korshinskiishrubs on the east-facing slope was significantly better than that of the shrubs on the south-facing slope. Furthermore, the growth status of theC.korshinskiishrubs on the lower slope position was slightly better than that of the shrubs on the upper or middle slope positions, and no significant difference was observed between the other slope positions. (2) Shallow soil moisture content had a direct, positive effect on the growth status ofC.korshinskii, and the shrub density and deep soil moisture content had a negative relationship with the growth status ofC.korshinskii. For the largeC.korshinskiishrubs, the height, crown diameter multiplied by the height, and canopy volume showed sensitive responses to the shallow soil moisture content. The results also demonstrated that the crown length responded with more sensitivity to the slope aspect and slope position than the other factors did, and the crown width showed a sensitive response to shrub density. (3) The topographical features (slope aspect, position, and gradient) and soil moisture content (shallow and deep soil moisture contents) determined the growth status ofC.korshinskiishrubs and explained the 59.9% variation in the growth changes. Furthermore, the variation of the growth changes in the slope aspect, slope position, and shallow soil moisture content were 21.1%, 16.0%, and 13.1%, respectively. Therefore, human-introduced vegetation restoration of the semi-arid loess regions should consider the spatial pattern of vegetation based on soil water conditions of the hillslope scale, and the scientific management of the middle and later periods of revegetation is also necessary to sustain theC.korshinskiishrub in the semi-arid loess region.

    Loess Plateau; hillslope scale; slope position; artificialCaraganakorshinskiishrub; soil moisture content

    國(guó)家自然科學(xué)基金項(xiàng)目(41401209,41561112);國(guó)家科技支撐計(jì)劃課題(2015BAC01B02)

    2017- 02- 16; < class="emphasis_bold">網(wǎng)絡(luò)出版日期

    日期:2017- 07- 03

    *通訊作者Corresponding author.E-mail: leiyang@rcees.ac.cn

    10.5846/stxb201702160264

    王子婷,楊磊,蔡國(guó)軍,莫保儒,柴春山,戚建莉,張洋東.半干旱黃土區(qū)坡面尺度檸條生長(zhǎng)狀況及影響要素分析.生態(tài)學(xué)報(bào),2017,37(23):7872- 7881.

    Wang Z T, Yang L, Cai G J, Mo B R, Chai C S, Qi J L,Zhang Y D.Spatial patterns ofCaraganakorshinskiigrowth on hillslope scale and influencing factors in the semi-arid Loess Plateau.Acta Ecologica Sinica,2017,37(23):7872- 7881.

    猜你喜歡
    坡位檸條冠幅
    無(wú)人機(jī)遙感影像提取的單木冠幅數(shù)據(jù)在桉樹林分蓄積量估測(cè)中的應(yīng)用1)
    城市綠地微環(huán)境對(duì)土壤動(dòng)物群落多樣性的影響
    檸條飼料林培育技術(shù)
    施肥對(duì)三江平原丘陵區(qū)長(zhǎng)白落葉松人工林中齡林單木樹冠圓滿度影響
    陜煤集團(tuán)神木檸條塔礦業(yè)公司
    我在西藏種檸條
    基于無(wú)人機(jī)高分影像的冠幅提取與樹高反演
    坡向坡位及郁閉度對(duì)森林公園內(nèi)林下南方紅豆杉生長(zhǎng)的影響
    綠色科技(2017年1期)2017-03-01 10:14:33
    不同坡位對(duì)毛竹林小氣候及地上生長(zhǎng)的影響
    檸條主題公園綠化美化措施
    国产日韩欧美亚洲二区| 日韩av不卡免费在线播放| 国产女主播在线喷水免费视频网站| 男的添女的下面高潮视频| 天天操日日干夜夜撸| 嫁个100分男人电影在线观看 | 精品国产一区二区久久| 国产一区二区三区综合在线观看| 亚洲欧美色中文字幕在线| 超碰97精品在线观看| 亚洲精品国产av成人精品| 亚洲色图综合在线观看| 大香蕉久久网| 亚洲国产日韩一区二区| 久久精品国产亚洲av涩爱| 国产一区二区三区综合在线观看| 国产亚洲欧美在线一区二区| 亚洲国产日韩一区二区| 香蕉丝袜av| 最近最新中文字幕大全免费视频 | 啦啦啦在线免费观看视频4| 欧美久久黑人一区二区| 国产人伦9x9x在线观看| 最近手机中文字幕大全| 女人被躁到高潮嗷嗷叫费观| 80岁老熟妇乱子伦牲交| 国产成人精品无人区| 99精品久久久久人妻精品| 日本wwww免费看| 欧美中文综合在线视频| 亚洲精品av麻豆狂野| 亚洲欧美激情在线| 国产精品香港三级国产av潘金莲 | 超碰成人久久| 久久久久久久大尺度免费视频| 少妇精品久久久久久久| 久久精品人人爽人人爽视色| 欧美成狂野欧美在线观看| svipshipincom国产片| 王馨瑶露胸无遮挡在线观看| 日韩中文字幕欧美一区二区 | 亚洲欧洲日产国产| 最近手机中文字幕大全| 男人舔女人的私密视频| 欧美黑人欧美精品刺激| 交换朋友夫妻互换小说| 国产男人的电影天堂91| 又黄又粗又硬又大视频| 飞空精品影院首页| 欧美老熟妇乱子伦牲交| 国产伦人伦偷精品视频| 极品人妻少妇av视频| 伦理电影免费视频| 97人妻天天添夜夜摸| 欧美日韩亚洲高清精品| 亚洲成av片中文字幕在线观看| 亚洲激情五月婷婷啪啪| 国产精品三级大全| 男人爽女人下面视频在线观看| 亚洲av在线观看美女高潮| 欧美黄色片欧美黄色片| 中文欧美无线码| 夫妻午夜视频| 日本欧美国产在线视频| 欧美精品高潮呻吟av久久| 性高湖久久久久久久久免费观看| 欧美日韩福利视频一区二区| 精品高清国产在线一区| 另类精品久久| 久久毛片免费看一区二区三区| 黄色一级大片看看| 欧美日韩福利视频一区二区| 国产99久久九九免费精品| 一级a爱视频在线免费观看| 99久久人妻综合| 亚洲精品美女久久久久99蜜臀 | 亚洲av综合色区一区| tube8黄色片| 黄色一级大片看看| 久久99精品国语久久久| 成人影院久久| 夜夜骑夜夜射夜夜干| 每晚都被弄得嗷嗷叫到高潮| www.999成人在线观看| 99国产精品一区二区蜜桃av | 男女无遮挡免费网站观看| 亚洲,欧美,日韩| 日韩欧美一区视频在线观看| 国产精品av久久久久免费| 中文字幕人妻熟女乱码| 国产成人a∨麻豆精品| 久久精品久久久久久久性| 欧美在线黄色| 日本a在线网址| 久久这里只有精品19| 黄片播放在线免费| 久久久久久免费高清国产稀缺| 成年动漫av网址| 午夜福利在线免费观看网站| 久久久久久亚洲精品国产蜜桃av| 亚洲自偷自拍图片 自拍| 精品亚洲成a人片在线观看| 久久鲁丝午夜福利片| 中国美女看黄片| 成年女人毛片免费观看观看9 | 国产成人一区二区在线| 一区二区三区精品91| 人人妻人人爽人人添夜夜欢视频| 国产高清视频在线播放一区 | 在线观看一区二区三区激情| 18禁观看日本| 欧美av亚洲av综合av国产av| 午夜免费观看性视频| 成年人免费黄色播放视频| 色精品久久人妻99蜜桃| 亚洲天堂av无毛| 性少妇av在线| 中文精品一卡2卡3卡4更新| 久久久精品94久久精品| 欧美国产精品一级二级三级| 少妇猛男粗大的猛烈进出视频| 日本欧美国产在线视频| 国产欧美日韩一区二区三区在线| 亚洲成人免费电影在线观看 | 国产熟女欧美一区二区| 一级,二级,三级黄色视频| 中文欧美无线码| 久久这里只有精品19| 国产男女内射视频| 欧美日韩国产mv在线观看视频| 亚洲五月色婷婷综合| 一个人免费看片子| 国产在线视频一区二区| 国产福利在线免费观看视频| 精品一区二区三区四区五区乱码 | 日韩一区二区三区影片| 午夜影院在线不卡| 亚洲精品美女久久久久99蜜臀 | 久久国产精品影院| 久久av网站| 91国产中文字幕| av在线播放精品| 老熟女久久久| 国产精品久久久久久人妻精品电影 | 亚洲国产日韩一区二区| 爱豆传媒免费全集在线观看| 侵犯人妻中文字幕一二三四区| 1024香蕉在线观看| 中文字幕人妻熟女乱码| 曰老女人黄片| 午夜福利一区二区在线看| 两性夫妻黄色片| 好男人视频免费观看在线| 国产日韩欧美在线精品| 在线观看免费日韩欧美大片| 亚洲欧美清纯卡通| 少妇被粗大的猛进出69影院| 国产一区二区在线观看av| 在线 av 中文字幕| 国产三级黄色录像| 国产精品 欧美亚洲| 老汉色∧v一级毛片| 久久中文字幕一级| 久久久久久久久免费视频了| 欧美在线黄色| 人人妻人人澡人人看| 国产老妇伦熟女老妇高清| 黄片小视频在线播放| 国产成人一区二区三区免费视频网站 | 老司机靠b影院| 中文字幕另类日韩欧美亚洲嫩草| 亚洲欧美日韩高清在线视频 | 婷婷成人精品国产| 国产福利在线免费观看视频| 久久综合国产亚洲精品| 国产精品久久久av美女十八| 美女主播在线视频| 久久久久国产精品人妻一区二区| 精品国产一区二区三区久久久樱花| 日本猛色少妇xxxxx猛交久久| 999精品在线视频| 香蕉丝袜av| 成人影院久久| 亚洲国产毛片av蜜桃av| 亚洲精品久久久久久婷婷小说| 真人做人爱边吃奶动态| 久久毛片免费看一区二区三区| 精品福利永久在线观看| 天堂中文最新版在线下载| 国产精品免费大片| 极品人妻少妇av视频| 电影成人av| 色精品久久人妻99蜜桃| 欧美黄色淫秽网站| 国产成人av激情在线播放| 国产精品99久久99久久久不卡| 老司机亚洲免费影院| 一级毛片电影观看| 日日夜夜操网爽| 国语对白做爰xxxⅹ性视频网站| 中文字幕精品免费在线观看视频| 国产精品久久久久成人av| 99久久99久久久精品蜜桃| 在线 av 中文字幕| 黑人巨大精品欧美一区二区蜜桃| 国产欧美日韩精品亚洲av| 丝袜美腿诱惑在线| 久久九九热精品免费| 亚洲激情五月婷婷啪啪| 欧美精品一区二区免费开放| av一本久久久久| 国产黄色视频一区二区在线观看| 欧美在线黄色| 99国产精品一区二区蜜桃av | 国产黄色免费在线视频| 国产伦人伦偷精品视频| 日韩伦理黄色片| 性少妇av在线| 一本—道久久a久久精品蜜桃钙片| 免费不卡黄色视频| 久久天躁狠狠躁夜夜2o2o | 中文字幕人妻丝袜制服| 免费观看人在逋| 99热网站在线观看| 中文字幕高清在线视频| 老司机靠b影院| 国产精品九九99| 国产精品秋霞免费鲁丝片| 色视频在线一区二区三区| 亚洲精品国产区一区二| 国产在线免费精品| 亚洲伊人久久精品综合| 老司机深夜福利视频在线观看 | 男女无遮挡免费网站观看| 精品国产一区二区三区久久久樱花| 交换朋友夫妻互换小说| 青青草视频在线视频观看| 99热国产这里只有精品6| 制服人妻中文乱码| 一本久久精品| 考比视频在线观看| 99精品久久久久人妻精品| 国产在视频线精品| 日韩 欧美 亚洲 中文字幕| 国产精品一国产av| 美女福利国产在线| 搡老岳熟女国产| 人成视频在线观看免费观看| 人人妻人人澡人人看| 在线观看国产h片| 欧美 亚洲 国产 日韩一| 欧美变态另类bdsm刘玥| 99热国产这里只有精品6| 男人爽女人下面视频在线观看| 国产精品久久久人人做人人爽| 青草久久国产| 国产成人免费无遮挡视频| 国产在线观看jvid| 老鸭窝网址在线观看| 久久久久国产一级毛片高清牌| 亚洲欧美一区二区三区久久| 韩国精品一区二区三区| 国产在线视频一区二区| 精品人妻熟女毛片av久久网站| 国产在线观看jvid| 国产精品国产av在线观看| 亚洲av在线观看美女高潮| 午夜福利视频在线观看免费| 黄色一级大片看看| 亚洲黑人精品在线| 女警被强在线播放| 看免费av毛片| 在线 av 中文字幕| 18禁观看日本| 亚洲成av片中文字幕在线观看| 国产色视频综合| 国精品久久久久久国模美| 成人国产一区最新在线观看 | 久久中文字幕一级| 两性夫妻黄色片| 一边亲一边摸免费视频| 美国免费a级毛片| 麻豆乱淫一区二区| 天天影视国产精品| 狠狠婷婷综合久久久久久88av| 777米奇影视久久| 亚洲专区国产一区二区| 尾随美女入室| 又大又黄又爽视频免费| 一级毛片我不卡| 精品国产一区二区三区四区第35| 久久99一区二区三区| 国产成人91sexporn| 精品一区二区三区av网在线观看 | 高清不卡的av网站| 国产免费现黄频在线看| 这个男人来自地球电影免费观看| 久久久久精品国产欧美久久久 | 色视频在线一区二区三区| 亚洲国产精品成人久久小说| 亚洲av综合色区一区| 50天的宝宝边吃奶边哭怎么回事| 欧美日韩亚洲国产一区二区在线观看 | 脱女人内裤的视频| 欧美成人精品欧美一级黄| 99精国产麻豆久久婷婷| 国产高清videossex| 每晚都被弄得嗷嗷叫到高潮| 成年美女黄网站色视频大全免费| 亚洲欧美日韩高清在线视频 | 国产成人欧美在线观看 | 中文字幕色久视频| 一边亲一边摸免费视频| 国产精品一区二区精品视频观看| 精品亚洲成a人片在线观看| 成在线人永久免费视频| 欧美日韩亚洲国产一区二区在线观看 | 午夜激情久久久久久久| 亚洲三区欧美一区| 99精国产麻豆久久婷婷| 精品一区在线观看国产| 伊人久久大香线蕉亚洲五| 黄频高清免费视频| 久久久久精品国产欧美久久久 | 七月丁香在线播放| 少妇人妻久久综合中文| 丝袜脚勾引网站| 一二三四社区在线视频社区8| 国产精品亚洲av一区麻豆| 一级a爱视频在线免费观看| 90打野战视频偷拍视频| 国产片内射在线| 亚洲图色成人| 免费观看av网站的网址| 青春草亚洲视频在线观看| 飞空精品影院首页| 交换朋友夫妻互换小说| 一本综合久久免费| 久久鲁丝午夜福利片| 精品久久蜜臀av无| 国产成人欧美| 又大又爽又粗| 午夜91福利影院| 亚洲人成网站在线观看播放| 一本色道久久久久久精品综合| 国产高清视频在线播放一区 | 18禁国产床啪视频网站| 国产色视频综合| 男女边吃奶边做爰视频| 后天国语完整版免费观看| 国产一区二区 视频在线| 又紧又爽又黄一区二区| 久久精品aⅴ一区二区三区四区| 天堂中文最新版在线下载| 女警被强在线播放| 久久精品久久精品一区二区三区| 午夜福利视频精品| 亚洲av成人不卡在线观看播放网 | 在线亚洲精品国产二区图片欧美| 亚洲国产日韩一区二区| 日韩人妻精品一区2区三区| 美女主播在线视频| 久久精品国产综合久久久| 国产精品香港三级国产av潘金莲 | 国产精品一区二区免费欧美 | 日韩,欧美,国产一区二区三区| 午夜免费成人在线视频| 19禁男女啪啪无遮挡网站| 国产国语露脸激情在线看| 91精品伊人久久大香线蕉| avwww免费| 亚洲精品久久久久久婷婷小说| 午夜福利视频精品| 亚洲成人国产一区在线观看 | av天堂在线播放| 大香蕉久久成人网| 婷婷色综合www| 精品少妇黑人巨大在线播放| 日韩大片免费观看网站| 黄色视频在线播放观看不卡| 国产欧美亚洲国产| 国产精品三级大全| av电影中文网址| 黑人欧美特级aaaaaa片| 精品人妻在线不人妻| 十分钟在线观看高清视频www| 女人被躁到高潮嗷嗷叫费观| 母亲3免费完整高清在线观看| 久久久久久久久免费视频了| 免费看不卡的av| 国产精品人妻久久久影院| 亚洲精品一区蜜桃| 亚洲天堂av无毛| 曰老女人黄片| 色94色欧美一区二区| 国产一区二区 视频在线| 男女午夜视频在线观看| 男人操女人黄网站| 亚洲国产中文字幕在线视频| 成年女人毛片免费观看观看9 | 国产精品免费大片| 少妇 在线观看| 丁香六月天网| 大话2 男鬼变身卡| 亚洲 欧美一区二区三区| 国产日韩欧美在线精品| 欧美日韩视频精品一区| 欧美人与性动交α欧美软件| 国产在线观看jvid| 在线观看免费日韩欧美大片| 一区二区三区乱码不卡18| www.av在线官网国产| 精品一区二区三区av网在线观看 | 韩国精品一区二区三区| 精品国产乱码久久久久久男人| 久久午夜综合久久蜜桃| 99久久人妻综合| 纵有疾风起免费观看全集完整版| 超碰成人久久| 一级毛片黄色毛片免费观看视频| 热re99久久国产66热| 久久久欧美国产精品| 亚洲国产欧美在线一区| 日本五十路高清| 亚洲熟女精品中文字幕| 国产精品99久久99久久久不卡| 少妇裸体淫交视频免费看高清 | 国产精品香港三级国产av潘金莲 | 精品亚洲乱码少妇综合久久| 1024视频免费在线观看| 午夜日韩欧美国产| 亚洲久久久国产精品| 丰满人妻熟妇乱又伦精品不卡| 我的亚洲天堂| 精品福利观看| 成年人免费黄色播放视频| av网站在线播放免费| 少妇粗大呻吟视频| 欧美xxⅹ黑人| 国产人伦9x9x在线观看| 下体分泌物呈黄色| 国产一区二区三区综合在线观看| 日韩大片免费观看网站| 免费看av在线观看网站| av网站免费在线观看视频| 一区二区三区四区激情视频| 久久人妻福利社区极品人妻图片 | 国产亚洲精品第一综合不卡| 国产成人一区二区三区免费视频网站 | 我要看黄色一级片免费的| 国产淫语在线视频| 中文字幕色久视频| 黄色毛片三级朝国网站| av在线老鸭窝| 18在线观看网站| 脱女人内裤的视频| 国产不卡av网站在线观看| 欧美日韩亚洲国产一区二区在线观看 | tube8黄色片| 母亲3免费完整高清在线观看| 99热网站在线观看| 亚洲美女黄色视频免费看| 亚洲伊人久久精品综合| 国产视频一区二区在线看| 久久天堂一区二区三区四区| 色婷婷av一区二区三区视频| 午夜91福利影院| 在线观看免费视频网站a站| 免费人妻精品一区二区三区视频| 中文字幕高清在线视频| 欧美精品高潮呻吟av久久| 国产黄频视频在线观看| 国产免费一区二区三区四区乱码| 天天影视国产精品| 免费人妻精品一区二区三区视频| 欧美精品av麻豆av| 午夜视频精品福利| 免费观看人在逋| 久久ye,这里只有精品| 女人被躁到高潮嗷嗷叫费观| 一级毛片黄色毛片免费观看视频| 少妇的丰满在线观看| 99热全是精品| 一级毛片电影观看| 下体分泌物呈黄色| 久久久欧美国产精品| 国产精品成人在线| 男女高潮啪啪啪动态图| 在线观看www视频免费| 91精品三级在线观看| 一级毛片 在线播放| 悠悠久久av| 久久精品亚洲av国产电影网| 99久久综合免费| av有码第一页| 国产一区二区 视频在线| 成人午夜精彩视频在线观看| 黄色片一级片一级黄色片| 欧美日韩av久久| 久久女婷五月综合色啪小说| 啦啦啦中文免费视频观看日本| 国产不卡av网站在线观看| 多毛熟女@视频| h视频一区二区三区| av网站在线播放免费| 91成人精品电影| 桃花免费在线播放| 乱人伦中国视频| 老司机午夜十八禁免费视频| 好男人电影高清在线观看| 69精品国产乱码久久久| 免费看av在线观看网站| 久久久久久久久久久久大奶| 夜夜骑夜夜射夜夜干| 欧美激情极品国产一区二区三区| 大片电影免费在线观看免费| 久久久久久人人人人人| 免费观看av网站的网址| 青青草视频在线视频观看| 国产不卡av网站在线观看| 日韩中文字幕视频在线看片| 999久久久国产精品视频| 天堂中文最新版在线下载| netflix在线观看网站| 日韩,欧美,国产一区二区三区| 久久99热这里只频精品6学生| 男女床上黄色一级片免费看| 我的亚洲天堂| 欧美黄色片欧美黄色片| 久久综合国产亚洲精品| 欧美黑人欧美精品刺激| 黄色a级毛片大全视频| 99国产精品免费福利视频| 久久久久久久精品精品| 亚洲欧美中文字幕日韩二区| 亚洲,欧美精品.| 国产高清国产精品国产三级| 国产免费又黄又爽又色| 国产成人91sexporn| www.自偷自拍.com| 久久99精品国语久久久| 极品人妻少妇av视频| 日韩一区二区三区影片| 纵有疾风起免费观看全集完整版| 欧美乱码精品一区二区三区| 亚洲中文av在线| 久久久久久久久免费视频了| 亚洲视频免费观看视频| 国产成人啪精品午夜网站| 一区福利在线观看| 美女扒开内裤让男人捅视频| 搡老岳熟女国产| 国产97色在线日韩免费| 国产野战对白在线观看| 亚洲国产av影院在线观看| 三上悠亚av全集在线观看| 日本色播在线视频| 晚上一个人看的免费电影| 美女主播在线视频| 超碰97精品在线观看| cao死你这个sao货| av国产久精品久网站免费入址| 亚洲男人天堂网一区| 久久久国产一区二区| 国产深夜福利视频在线观看| 国产精品偷伦视频观看了| 中文字幕精品免费在线观看视频| 成年美女黄网站色视频大全免费| 777久久人妻少妇嫩草av网站| 久久精品成人免费网站| 免费日韩欧美在线观看| 欧美老熟妇乱子伦牲交| 丝袜美足系列| 晚上一个人看的免费电影| 人体艺术视频欧美日本| av天堂久久9| 日韩 欧美 亚洲 中文字幕| 久久久精品国产亚洲av高清涩受| 久久人人爽人人片av| 免费观看av网站的网址| 亚洲av电影在线观看一区二区三区| 亚洲国产精品国产精品| 两性夫妻黄色片| 热99国产精品久久久久久7| www.av在线官网国产| av一本久久久久| 国产在线观看jvid| 免费不卡黄色视频| 亚洲图色成人| 啦啦啦在线观看免费高清www| 国产精品国产av在线观看| 欧美久久黑人一区二区| 午夜福利免费观看在线| 在现免费观看毛片| 这个男人来自地球电影免费观看| 亚洲精品国产av蜜桃| 夫妻午夜视频| 久久中文字幕一级| 晚上一个人看的免费电影| 久久久精品94久久精品| 香蕉丝袜av| 亚洲国产欧美日韩在线播放| 欧美中文综合在线视频| 美女脱内裤让男人舔精品视频| 婷婷成人精品国产| 日韩av免费高清视频| 亚洲,欧美,日韩| 天天躁日日躁夜夜躁夜夜| 交换朋友夫妻互换小说| 亚洲精品在线美女| 中文字幕人妻熟女乱码| av天堂在线播放| 可以免费在线观看a视频的电影网站| 天天影视国产精品| 欧美精品av麻豆av|