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

    南方紅壤區(qū)植被結(jié)構(gòu)類型與降雨模式對林下水土流失的影響

    2020-04-24 08:27:08張海東于東升馬利霞陸曉松
    農(nóng)業(yè)工程學報 2020年5期
    關(guān)鍵詞:紅壤土壤侵蝕徑流

    陳 洋,張海東,于東升,馬利霞,陸曉松,王 鑫

    南方紅壤區(qū)植被結(jié)構(gòu)類型與降雨模式對林下水土流失的影響

    陳 洋1,3,張海東2,于東升1,3※,馬利霞1,陸曉松1,3,王 鑫1,3

    (1. 土壤與農(nóng)業(yè)可持續(xù)發(fā)展國家重點實驗室,中國科學院南京土壤研究所,南京 210008; 2. 蘇州市農(nóng)業(yè)科學院(江蘇太湖地區(qū)農(nóng)業(yè)科學研究所),蘇州 215155; 3. 中國科學院大學,北京 100049)

    植被和降雨是水土流失的關(guān)鍵因素,探究二者對水土流失的影響對開展水土保持具有重要意義。該研究基于鷹潭紅壤生態(tài)試驗站5種植被結(jié)構(gòu)類型的徑流小區(qū)2016-2018年93次降雨、徑流、泥沙觀測資料以及各小區(qū)植被結(jié)構(gòu)參數(shù),利用自組織映射(self-organizing maps,SOM)方法,根據(jù)雨量、歷時、60 min最大雨強、平均雨強、降雨集中性等特征指標劃分降雨模式,研究了不同降雨模式和植被結(jié)構(gòu)類型的水土流失特征,并采用冗余分析(RDA)定量研究降雨與植被對林下水土流失的影響。結(jié)果表明,SOM方法能客觀識別紅壤區(qū)4種典型侵蝕降雨模式,RIII模式(短歷時、大雨強、雨量集中)是造成水土流失的主要降雨模式,RIV模式(多雨量、大雨強、長歷時)最具侵蝕性破壞力;植被結(jié)構(gòu)類型顯著影響水土流失,水土保持功能從大到小依次為:灌草混交林、草地、低灌林、喬木林、高灌林。RDA分析表明,降雨模式與植被結(jié)構(gòu)類型能夠改變降雨、植被對水土流失的影響,隨著降雨模式由弱到強轉(zhuǎn)變,植被的水土保持功能逐漸減小,降雨影響增強,水土流失由植被主控演變?yōu)槠胶饪刂啤⒔涤曛骺?;隨著植被結(jié)構(gòu)類型由近地表植被向灌木、喬木的變化,植被的調(diào)節(jié)能力減小,降雨影響增強。

    侵蝕;徑流;土壤;植被結(jié)構(gòu);降雨模式;南方紅壤區(qū)

    0 引 言

    南方紅壤區(qū)由于豐富降水、山地丘陵地貌、高強度土地利用、單一林分結(jié)構(gòu),易于造成嚴重水土流失[1]。據(jù)2013年第一次全國水利普查水土保持情況公報顯示,紅壤區(qū)水土流失面積達16萬km2,平均侵蝕速率為3 419.8 t/(km2·a),部分流域可達237.6 t/(km2·a)以上[2-4],遠超區(qū)域容許土壤流失量的上限(500 t/(km2·a)),嚴重威脅區(qū)域土壤資源利用與可持續(xù)發(fā)展。自20世紀80年代以來,紅壤區(qū)采取了?;~塘、護坡、飛播造林、退耕還林還草等一系列水土保持生態(tài)治理措施[1,4],區(qū)域水土流失得到有效遏制,但“空中綠化”“林下流”的特征依然明顯[1]。探究植被結(jié)構(gòu)類型與降雨模式對區(qū)域水土流失的綜合影響,對提高土壤侵蝕治理效率、促進水源涵養(yǎng)具有重要意義。

    水土流失景觀差異與地表水文-植被過程響應密切相關(guān),降雨是這一過程中的動力來源[5-6]。降雨模式特征如雨量、強度、歷時、雨強變異性、集中度等,決定了地表的水文過程[7];在水力侵蝕過程中,影響著滲流潛力、土壤剝離速率以及產(chǎn)流過程等[8]。Parsons等[9]模擬了5種降雨模式下的土壤侵蝕過程,表明雨強上升型的降雨模式對土壤剝離速率更高、能力更強。Nciizah等[10]研究指出,濕潤區(qū)短歷時的強降雨模式容易快速產(chǎn)流并引起地表侵蝕。Duan等[11]運用K-means法對自然降雨事件聚類分析,指出低強度和長歷時的降雨模式能充分產(chǎn)生地下徑流,從而使表層土壤免受侵蝕。Chen等[12]基于降雨量集中性差異,分析得出黃土高原地區(qū)雨量集中性高的降雨模式是造成土壤侵蝕的主要降雨模式。這些研究綜合分析了多維降雨因子特征與水土流失之間的關(guān)系,揭示了區(qū)域水土流失的主導降雨模式。

    植被對控制水土流失發(fā)揮重要作用,不同結(jié)構(gòu)形態(tài)植被的水土保持功能不一。Peng等[13]在西南地區(qū)的自然降雨觀測表明,覆蓋度較高的草本比喬木林地能更有效減少水土流失;于東升等[14]通過模擬降雨試驗表明,林下有草或灌草的混交植被結(jié)構(gòu)具有最佳的雨水調(diào)節(jié)與防洪減沙效能;孫佳佳等[15]在閩西北紅壤徑流小區(qū)試驗分析指出,林草結(jié)構(gòu)比純林結(jié)構(gòu)能減少20%的徑流和90%的土壤侵蝕。顯然,近地表植被對降雨侵蝕的控制作用最為明顯,多層次植被結(jié)構(gòu)比單一層狀植被的水土保持作用更具優(yōu)勢[16]。

    然而,已有研究對降雨模式劃分相對混雜、主觀性強,對不同降雨模式的侵蝕破壞力認識不足;并且多側(cè)重于獨立探究降雨、植被對坡面水土流失的影響,忽視了自然環(huán)境下植被、降雨對水土流失的交互競爭性影響。因此,本文針對紅壤區(qū)復雜降雨模式、典型植被結(jié)構(gòu)[11],通過徑流小區(qū)試驗,運用自組織映射(self-organizing maps,SOM)和多元統(tǒng)計方法,研究揭示水土流失對不同類型降雨模式、植被結(jié)構(gòu)的響應特征,定量分析它們對水土流失的綜合影響,并明確植被在不同降雨條件下的截留減蝕作用差異,為區(qū)域有針對性開展水土流失治理提供科學依據(jù)。

    1 材料與方法

    1.1 研究區(qū)域

    試驗區(qū)位于江西省鷹潭市中國科學院紅壤生態(tài)實驗站(116°55′16″~116°56′18″E,28°12′01″~28°12′37″N),屬東亞季風氣候,年平均氣溫17.6 ℃,降水量達1 795 mm,>10 ℃的積溫達7 500 ℃,無霜期236 d[17]。試驗區(qū)為第四紀紅黏土發(fā)育的紅壤丘崗地形,坡度小于15°,海拔30~47 m。土壤母質(zhì)以紅砂巖為主,風化殼深厚,土層為A-C結(jié)構(gòu)型,粉粒含量為30.3%~40.8%,質(zhì)地黏重、易侵蝕[17]。歷史上形成的亞熱帶常綠闊葉林、針闊混交林等自然植被景觀,由于近現(xiàn)代以來森林砍伐、土地開墾,導致原生植被嚴重破壞。當前試驗區(qū)周邊植被為人工次生林,成林時間小于30a,代表性物種為馬尾松()、芒萁()、白茅()等[17-18]。

    1.2 試驗小區(qū)觀測

    在試驗站東區(qū)坡度為8°的緩坡上建有8個長×寬水平投影為10 m×2.6 m的徑流試驗小區(qū)(圖1)[18]?;诩t壤區(qū)典型次數(shù)植被配置特征、形態(tài)結(jié)構(gòu)與植被參數(shù),設計5種植被結(jié)構(gòu)類型:純草(G)、灌草(GS)、低灌(LS)、高灌(HS)、喬木(F),物種為馬尾松、芒萁、白茅。除了G小區(qū)外,其他小區(qū)均布植馬尾松21棵。

    在試驗區(qū)空曠地帶安裝自動雨量計記錄自然降雨過程。每次降雨后,利用水尺測定各小區(qū)徑流桶中水深,并采集徑流樣品和所有泥沙,待泥沙風干后測定質(zhì)量;再用烘干法測定徑流中的泥沙含量和風干泥沙樣中泥沙含水量,利用式(1)~(3)計算各次降雨后的徑流深度與土壤流失量[19]。本文采用2016年1月至2018年12月的試驗觀測數(shù)據(jù),同期測定的各試驗小區(qū)平均的土壤理化性質(zhì)及植被結(jié)構(gòu)參數(shù)(表1)。

    R=V/(2)

    式中R為第次降雨產(chǎn)生的徑流深,mm;S為土壤侵蝕量,g/m2;為徑流小區(qū)面積,m2;V為徑流體積,m3;h為測定的水深,G為泥沙風干質(zhì)量,g;C為徑流中的泥沙含量,g/m3;C為泥沙含水率,%;為圓柱體徑流桶的半徑,m;為土壤密度,g/m3。

    表1 試驗小區(qū)植被結(jié)構(gòu)參數(shù)及土壤性質(zhì)

    注:P為馬尾松;D為芒萁;I為白茅。

    Note: P, D, I represent,-,, respectively.

    1.3 基于SOM的降雨模式聚類

    定義產(chǎn)生徑流的降雨為侵蝕降雨[11-13]。針對各場次侵蝕降雨,采用降雨量()、歷時()、最大60 min雨強(max)、平均雨強(mean)、降雨量集中性指數(shù)()等指標描述降雨模式特征[11-13],利SOM方法對降雨特征因子集進行聚類,從而識別降雨模式。

    SOM是Kohonen[20]教授發(fā)明的一種基于神經(jīng)網(wǎng)絡的數(shù)據(jù)降維技術(shù),其網(wǎng)絡結(jié)構(gòu)為一個輸入層和一個競爭層,前者由接收高維數(shù)據(jù)的輸入單元組成,后者由位于規(guī)則二維網(wǎng)格上的輸出單元組成,網(wǎng)絡層之間通過權(quán)向量連接。在聚類分析過程中,先采用z-score函數(shù)對降雨特征因子進行預處理,以消除數(shù)據(jù)噪聲影響;再將處理后的變量集導入輸入層,運用高斯鄰域函數(shù)構(gòu)建8×10的SOM網(wǎng)絡拓撲結(jié)構(gòu),設置初始學習效率為0.1,迭代次數(shù)1 000次,計算輸入變量與所有神經(jīng)元權(quán)重之間的距離;最后,基于最近距離原則從競爭層中獲得最佳匹配單元,從而得到降雨模式聚類結(jié)果。

    聚類數(shù)是一平衡選擇。本著聚類數(shù)小且盡可能反映組間差異原則選擇最優(yōu)聚類結(jié)果,并利用組內(nèi)總方差(within groups sum of squared error, WSS)與聚類數(shù)總體變化關(guān)系曲線的拐點對聚類數(shù)進行客觀判斷。

    1.4 數(shù)據(jù)分析

    采用單因素方差分析法研究不同降雨模式、植被結(jié)構(gòu)類型下的徑流量、土壤侵蝕量之間差異,當檢驗值<0.05時具有統(tǒng)計學意義;冗余分析RDA(redundancy analysis)用以揭示植被與降雨對水土流失的影響。

    2 結(jié)果與分析

    2.1 不同模式降雨統(tǒng)計特征

    研究區(qū)2016-2018年記錄到93次侵蝕降雨事件,降雨量合計4 226.5 mm,占總降雨量83.7%。規(guī)定單次降雨后小區(qū)總產(chǎn)沙量小于0.15 kg為無侵蝕,觀測到有52次降雨產(chǎn)生了明顯土壤侵蝕,占侵蝕降雨次數(shù)的55.9%。SOM聚類分析表明,WSS隨聚類數(shù)的變化由陡變緩,聚類數(shù)4是其總體變化拐點(圖1)。因此,研究區(qū)主要存在4種侵蝕降雨模式,分別記作RI、RII、RIII和RIV(表2)。

    RI在研究區(qū)發(fā)生頻率較高(31.2%),但引起土壤侵蝕頻率較低(13.8%),其特點是雨量和強度小、歷時長、雨量不集中,多發(fā)生在9月至次年4月。RII模式在4-10月均有發(fā)生,特點為雨量較多、中雨強、歷時長,該模式降雨頻率少(14.0%),但易于造成土壤侵蝕(76.9%)。RIII為短歷時、高強度、雨量集中的暴雨,主要發(fā)生在5-9月,是造成水土流失的首要降雨模式(43.0%),產(chǎn)生土壤侵蝕頻率較高(75.0%)。RIV為多雨量、長歷時、大強度的連續(xù)強降雨,集中發(fā)生在6-7月,發(fā)生頻率較低(11.8%),但具有較高的土壤侵蝕頻率(72.7%)。4種降雨模式中以RIII、RIV的降雨量最豐富,分別占降雨總量的31.5%、27.0%,RI、RII模式降雨量分別占20.9%、20.7%(表2)。

    圖1 組內(nèi)總方差隨聚類數(shù)目的變化

    Fig.1. Variation of within groups sum of squared error with number of clusters

    表2 4種模式降雨的統(tǒng)計特征

    注:RI模式指雨量和強度小、歷時長、雨量不集中;RII模式指雨量較多、中雨強、歷時長;RIII模式指短歷時、大雨強、雨量集中;RIV模式指多雨量、大雨強、長歷時。下同。

    Note: RIpattern refers to rainfall events with low intensity, long duration, minor amount of rainfall and low rainfall concentration; RIIpattern refers to the rainfall events with medium amount of rainfall, moderate rainfall intensity and long duration; RIIIpattern refers to the rainfall events with characteristics of short duration, strong rainfall intensity and high rainfall concentration; RIVpattern refers to great amount of rainfall, strong rainfall intensity and long duration. The same as follows.

    研究區(qū)存在4種典型侵蝕降雨模式,與Duan等[11,21]基于D、D、max等降雨特征因子的組間差異分析得到的湘北、贛北紅壤區(qū)降雨模式特征一致,表明本結(jié)果在中亞熱帶紅壤區(qū)具有典型性。秦偉等[22]基于D、D、mean等降雨因子分析認為,贛北紅壤地區(qū)侵蝕降雨存在3種模式,并不存在RIV模式;Liu等[23]基于快速聚類法在同樣地區(qū)分析得出5種降雨模式,其中部分模式之間具有相似特征。顯然,已有研究對降雨模式劃分并不一致,這不僅與區(qū)域氣候條件密切相關(guān),也受聚類分析方法影響。這些研究[11,21-23]均基于K-means聚類方法,以特征因子與聚類中心的距離來識別降雨模式,但需要事先設定聚類數(shù)目。通常未知降雨模式組間方差情況下自主設定聚類數(shù),易導致對降雨模式的劃分不充分或冗余。本研究根據(jù)WSS的變化特征客觀給出聚類數(shù),獲得的4種降雨模式,準確反映了中亞熱帶紅壤區(qū)的侵蝕降雨特征。

    2.2 不同植被結(jié)構(gòu)類型下水土流失特征

    單因素方差分析顯示,不同植被結(jié)構(gòu)類型下年均徑流量與土壤侵蝕量存在一定差異(<0.05)(圖2)。HS、F小區(qū)的土壤侵蝕量最高,分別為G、GS、LS小區(qū)2.1~3.2倍;HS小區(qū)的徑流量最高,達155.6 mm,顯著高出T、LS、GS、G小區(qū)的0.6、0.6、1.2、0.8倍(<0.05)。5種植被結(jié)構(gòu)類型下產(chǎn)生的水土流失量從大到小排序為GS、G、LS、T、HS。由于近地表幼齡植被的活躍生長改變了土壤侵蝕敏感性[24],導致各植被結(jié)構(gòu)類型下土壤侵蝕量的年際變化并不同步,而徑流量波動在各小區(qū)具有一致性。

    注:植被結(jié)構(gòu)F、G、GS、LS、HS分別指喬木、純草,灌草,低灌、高灌。不同小寫字母表示處理間差異顯著。下同。

    植被結(jié)構(gòu)改變降雨到達地表的過程和土壤對降雨的響應,進而影響水土流失[25-26]。冠層與枝干通過截留雨滴、增加蒸散以減少到達地表的凈降雨量,并根據(jù)植株高度、冠層開度分配降雨動能[26],因而高植被覆蓋的G、GS小區(qū)(表1)產(chǎn)生的徑流和土壤侵蝕小于低植被覆蓋的HS小區(qū)。針簇狀結(jié)構(gòu)的馬尾松冠層能將小雨滴合并成大雨滴,從而增加雨滴動能[27-28],當截留飽和后雨滴從冠層高處降落產(chǎn)生較大的侵蝕力,這也是HS、T小區(qū)中侵蝕模數(shù)較大的原因。低灌、草本等植被結(jié)構(gòu)不僅具有發(fā)達的淺層根系,還能通過改善地表水熱條件支持地表生物生長,豐富的土壤生物和根系活動改善土壤孔隙度和團聚體,增加土壤入滲潛力和抗侵蝕能力[29]。相較而言,HS、T小區(qū)表層土壤孔隙度較少、下滲渠道不足,因而產(chǎn)生較多地表徑流。

    2.3 不同降雨模式下的水土流失特征

    4種降雨模式中,RIII和RIV產(chǎn)生的地表徑流量較多,分別占總徑流量的35.6%~44.4%和26.9%~32.9%,是RI、RII模式的3 倍和1倍以上(圖3)。RIII模式對土壤侵蝕量貢獻率最大,其次為RIV模式,分別占侵蝕總量的48.2%~56.7%以及25.2%~35.6%;RI和RII模式產(chǎn)生土壤侵蝕量較少,平均分別占3.4%、13.4%。4種降雨模式對試驗區(qū)水土流失的貢獻大小依次為RIII>RIV>RII>RI。

    各降雨模式下歷次降雨產(chǎn)生的平均徑流與土壤侵蝕量存在顯著差異(圖4)。RIV模式下產(chǎn)生的平均土壤侵蝕量最多,顯著高出其他模式1.5~8.7倍(<0.05),其次是RII、RIII和RI模式(圖4a)。RIV模式產(chǎn)生的平均徑流也最多,顯著高出RIII模式的1.3~8.7倍(<0.05);RIII模式產(chǎn)生的平均徑流次之,高出RI、RII模式的0.8~2.5倍。各降雨模式下歷次降雨產(chǎn)生的平均徑流量從大到小依次為RIV、RIII、RII、RI,平均土壤流失量從大到小則為RIV、RII、RIII、RI,且在各植被結(jié)構(gòu)類型下表現(xiàn)具有一致性。

    圖3 不同降雨模式下水土流失總量

    圖4 不同降雨模式下次降雨產(chǎn)生的平均徑流和土壤侵蝕量

    RIII是引起紅壤區(qū)水土流失的最主要降雨模式。研究區(qū)紅壤飽和入滲速率為122.3~489.0 mm/h[17],遠高于大部分降雨事件的雨強,紅壤對該類降雨響應時間短,導致降雨未充分下滲及土壤飽和情況下形成地表徑流。同時,此類高強度降雨能迅速破壞地表介殼、剝離土壤顆粒,進而造成土壤侵蝕。RIII模式特征與其他地區(qū)相比具有共同性,如喀斯特巖溶區(qū)[13]、北方土石山區(qū)[30]、黃土高原[12]、中南地區(qū)[11,21]。

    RIV是區(qū)域氣候環(huán)境水平上最具破壞力的降雨模式。該模式具有長歷時、大雨強、多雨量特征,前期降雨補充土壤水分后使土壤入滲能力持續(xù)保持在較低水平[11,22];植被對該模式降雨也會在短時間內(nèi)形成飽和截留,無法繼續(xù)有效攔截降雨[21,31],因而具有較高的產(chǎn)流能力(圖4b)。該降雨模式下,土壤侵蝕類型隨著降雨歷時延長逐步由降雨侵蝕變化為徑流與降雨綜合侵蝕,進而引發(fā)更強烈的土壤流失[32]。由于該極端降雨具有較大的土壤剝離與徑流輸送能力,容易引起泥石流、洪澇等災害,是區(qū)域應積極防范的災害性天氣。

    2.4 不同降雨模式下水土流失隨植被結(jié)構(gòu)類型的分異特征

    降雨和植被對水土流失過程發(fā)揮著互相競爭性影響,不同的植被結(jié)構(gòu)與降雨模式疊加將產(chǎn)生復雜的產(chǎn)流、產(chǎn)沙情景[21]。在RI模式下,各小區(qū)水土流失量較小且不存在組間差異;RII模式下其差異顯著(<0.05),但絕對相差依然較小。在RIII模式下,HS和F小區(qū)的徑流量比GS、G小區(qū)高出0.62~1.71倍,土壤侵蝕量高出2.30~4.63倍;在RIV模式下,HS、F小區(qū)的徑流量比GS、G小區(qū)高出0.32~0.78倍,土壤侵蝕量則高出3.37~7.42倍(圖5)。顯然,不同植被結(jié)構(gòu)類型下的水土流失對降雨模式的響應存在差異。

    圖5 不同降雨模式下水土流失隨植被結(jié)構(gòu)類型的分異特征

    3 討 論

    3.1 不同降雨模式下植被對水土流失的影響

    RDA分析顯示,不同降雨模式下降雨與植被對水土流失的解釋方差并不一致(圖6),表明降雨和植被共同影響,改變了地表產(chǎn)流、產(chǎn)沙的響應。RI模式下,植被、降雨對侵蝕量的解釋方差占比為0.37、0.56,植被是侵蝕發(fā)生的首要控制因素;RII模式下,二者的解釋力相接近,分別為0.57、0.42,表明降雨和植被平衡控制土壤侵蝕;RIII、RIV模式下,植被、降雨的解釋能力依次為0.22、0.63和0.14、0.74,降雨是侵蝕發(fā)生的首要控制因素(圖6a)。就徑流而言,無論在何種降雨模式下,降雨的影響力均高于植被,并且二者作用力大小也隨降雨模式變化(圖6b)。

    圖6 不同降雨模式下降雨和植被分別對土壤侵蝕與徑流的影響

    降雨、植被對水土流失產(chǎn)生過程有著復雜競爭性關(guān)系。Fang等[33]在北亞熱帶紅壤區(qū)研究表明,不同降雨模式下降雨對水土流失的影響高于植被因素;Wei等[34]在黃土高原地區(qū)研究顯示,降雨與植被隨著降雨模式變化交替主控水土流失;Huang等[21]定性分析了紅壤區(qū)降雨與植被對水土流失競爭性作用,認為分3種情形:降雨主控、植被主控與綜合控制。本研究RDA分析表明,植被與降雨對水土流失的影響力并不固定,隨降雨模式的變化而改變,這與上述研究結(jié)果相近,但略有差異。這是由于Fang等[33]基于主觀聚類法劃分的3種降雨模式中,忽視了弱降雨模式存在;而黃土高原地區(qū)降雨模式較為簡單,對降雨、植被與水土流失的關(guān)系較為單一。相較而言,本文與Huang等[21]研究結(jié)果能夠較好地反映出紅壤區(qū)降雨、植被對水土流失影響力的可變性。

    已有研究關(guān)于植被的水土保持能效忽略了降雨模式的影響。本研究基于降雨模式的考慮,能更準確評價植被抑制侵蝕的作用,并認識植被與降雨之間的競爭性。隨著降雨模式由RI變化為RII乃至RIII、RIV,植被抑制水土流失的作用逐漸減小,這取決于植被結(jié)構(gòu)對不同的降雨的響應。在小雨強、長歷時的降雨條件下,植被的截留與削弱雨滴動能能力最好,而在大雨強、短歷時的降雨條件下較差。Li等[35]在閩西北紅壤林地研究表明,馬尾松灌木對小雨強(7 mm/h)降雨的截留率最高,一定雨強下(30 mm/h)截留趨近飽和;且降雨動能與雨強呈指數(shù)關(guān)系,當雨強達到一定閾值(14 mm/h)后,植被對降雨動能的削減達到飽和。這些研究更加證實了植被水土保持功能對降雨模式的依賴性。

    3.2 不同植被結(jié)構(gòu)類型下降雨對水土流失的影響

    不同植被結(jié)構(gòu)類型下降雨和植被對水土流失影響存在差異(圖7)。從土壤侵蝕方面來看(圖7a),GS植被結(jié)構(gòu)下降雨的解釋方差為0.57,植被為0.21,降雨的影響是植被的2.71倍;G植被結(jié)構(gòu)下降雨影響是植被的3.57倍。在單一灌、喬(LS、HS、F)植被結(jié)構(gòu)下,降雨影響更大,分別是植被的5.69、19.33、10.50倍。試驗區(qū)所涉及的植被結(jié)構(gòu)類型中,降雨對侵蝕的影響大于植被,這與Mishra等[36]、Yu等[37]分別基于全球、徑流小區(qū)尺度的研究結(jié)果類似。另從徑流方面來看(圖7b),除了G、GS植被條件下植被的解釋方差高于降雨外,其他條件下則相反;G、GS、LS、F、HS植被結(jié)構(gòu)類型下,降雨的解釋方差分別是植被的0.64、0.60、1.64、7.33、7.00倍。表明區(qū)域降雨條件下,植被的調(diào)節(jié)作用在近地表(G、GS)結(jié)構(gòu)下比單一喬、灌木(LS、HS、F)結(jié)構(gòu)更強;相應地,降雨作用在喬、灌木植被結(jié)構(gòu)下比近地表植被結(jié)構(gòu)下更強。

    圖7 不同植被結(jié)構(gòu)類型下降雨和植被對水土流失的影響

    同一氣候降雨水平上,植被類型及其結(jié)構(gòu)形態(tài)是引起土壤侵蝕空間異質(zhì)性的主要來源。已有研究證實了不同植被類型抗侵蝕功能性差異,但很少關(guān)注其調(diào)節(jié)能力大小[13]。Zhou等[38]在運用結(jié)構(gòu)方程模型揭示了恢復期長的植被對水土流失調(diào)節(jié)作用比恢復期短的植被更大。Yu等[37]利用網(wǎng)絡分析和RDA分析量化了降雨、植被等要素對水土流失影響,指出高覆蓋、球形冠層植被對水土流失的調(diào)節(jié)能力比傘形冠層植被更強。本研究針對紅壤區(qū)典型植被特征,量化了不同植被結(jié)構(gòu)類型的水土保持功能差異,進一步證實了近地表植被具有更佳的抑制侵蝕、減少徑流作用。同時表明,植被結(jié)構(gòu)差異誘導了降雨-水土流失過程,從而影響降雨、植被對水土流失的競爭性影響。

    4 結(jié) 論

    1)運用自組織映射方法確立了紅壤區(qū)存在4種典型侵蝕降雨模式,即:RI模式(少雨量、弱雨強、長歷時);RII模式(中雨量、中雨強、長歷時);RIII模式(短歷時、大雨強、雨量集中);RIV模式(多雨量、大雨強、長歷時)。其中RIII是引起區(qū)域水土流失的最主要降雨模式;RIV模式最具侵蝕破壞力。

    2)不同植被結(jié)構(gòu)類型的水土保持功效從大到小依次為:灌草混交林(GS)、草地(G)、低灌林(LS)、喬木林(F)、高灌林(HS),這種差異主要由植被覆蓋、結(jié)構(gòu)形態(tài)的調(diào)控作用實現(xiàn)。

    3)紅壤區(qū)林下水土流失過程受降雨與植被綜合影響,二者對水土流失的競爭性影響依賴于降雨模式、植被結(jié)構(gòu)類型而呈現(xiàn)一定規(guī)律性。在降雨模式水平上,隨著降雨模式由弱至強(RI?RII?RIII?RIV)變化,降雨的影響增加,植被抑制水土流失的能力逐漸減小,相應地水土流失由植被主控演變?yōu)槠胶饪刂?、降雨主控。在植被結(jié)構(gòu)類型水平上,降雨影響總體高于植被,隨著植被結(jié)構(gòu)類型由近地表向單一喬木、灌木結(jié)構(gòu)變化(GS?G?F?LS?HS),植被調(diào)節(jié)功能減小,降雨影響則增加。通過冗余分析定量研究了降雨與植被在不同降雨模式及植被結(jié)構(gòu)下對水土流失的影響,對紅壤區(qū)水土流失治理具有積極指導意義。

    [1] 何圣嘉,謝錦升,楊智杰,等. 南方紅壤丘陵區(qū)馬尾松林下水土流失現(xiàn)狀、成因及防治[J]. 中國水土保持科學,2011,9(6):65-70.

    He Shengjia, Xie Jinsheng, Yang Zhijie, et al. Status, causes and prevention of soil and water loss in Pinus massoniana woodland in hilly red soil region of southern China[J]. Science of Soil and Water Conservation, 2011, 9(6): 65-70. (in Chinese with English abstract)

    [2] 梁音,張斌,潘賢章,等.南方紅壤丘陵區(qū)水土流失現(xiàn)狀與綜合治理對策[J]. 中國水土保持科學,2008,6(1):22-27.

    Liang Yin, Zhang Bin, Pan Xianzhang, et al. Current status and comprehensive control strategies of soil erosion for hilly region in the Southern China[J]. Science of Soil and Water Conservation, 2008, 6(1): 22-27. (in Chinese with English abstract)

    [3] 水利部. 第一次全國水利普查水土保持情況公報[R]. 北京:水利部,2013.

    [4] 景可,焦菊英,李林育,等.中國南方紅壤丘陵區(qū)土壤侵蝕量與泥沙輸移比研究[J]. 水利學報,2010,41(12):1476-1482.

    Jing Ke, Jiao Juying, Li Llinyu, et al. Research on soil erosion amount and sediment delivery ratio in the hilly red soil region of southern China[J]. Journal of Hydraulic Engineering, 2010, 41(12): 1476-1482. (in Chinese with English abstract)

    [5] Zhang K, Wen Z. Review and challenges of policies of environmental protection and sustainable development in China[J]. Journal of environmental management, 2008, 88(4): 1249-1261.

    [6] Ravi S, Breshears D D, Huxman T E, et al. Land degradation in drylands: Interactions among hydrologic–aeolian erosion and vegetation dynamics[J]. Geomorphology, 2010, 116(3/4): 236-245.

    [7] Ludwig J A, Wilcox B P, Breshears D D, et al. Vegetation patches and runoff-erosion as interacting ecohydrological processes in semiarid landscapes[J]. Ecology, 2005, 86(2): 288-297.

    [8] Foley J L, Silburn D M. Hydraulic properties of rain impact surface seals on three clay soils-influence of raindrop impact frequency and rainfall intensity during steady state[J]. Soil Research, 2002, 40(7): 1069-1083.

    [9] Parsons A J, Stone P M. Effects of intra-storm variations in rainfall intensity on interrill runoff and erosion[J]. Catena, 2006, 67(1): 68-78.

    [10] Nciizah A D, Wakindiki I I C. Soil sealing and crusting effects on infiltration rate: A critical review of shortfalls in prediction models and solutions[J]. Archives of Agronomy and Soil Science, 2015, 61(9): 1211-1230.

    [11] Duan J, Yang J, Tang C, et al. Effects of rainfall patterns and land cover on the subsurface flow generation of sloping Ferralsols in southern China[J]. PloS one, 2017, 12(8): e0182706.

    [12] Chen H, Zhang X, Abla M, et al. Effects of vegetation and rainfall types on surface runoff and soil erosion on steep slopes on the Loess Plateau, China[J]. Catena, 2018, 170: 141-149.

    [13] Peng T, Wang S. Effects of land use, land cover and rainfall regimes on the surface runoff and soil loss on karst slopes in southwest China[J]. Catena, 2012, 90: 53-62.

    [14] 于東升,史學正. 紅壤區(qū)不同生態(tài)模式的“土壤水庫”特征及其防洪減災效能[J]. 土壤學報,2003,40(5):656-664.

    Yu Dongsheng, Shi Xuezheng. Characteristics of soil-water reservoir and their efficiency in flooding control under different ecological patterns of red soils in the south of China[J]. Acta Pedologica Sinica, 2003, 40(5): 656-664. (in Chinese with English abstract)

    [15] 孫佳佳,于東升,史學正,等. 植被葉面積指數(shù)與覆蓋度定量表征紅壤區(qū)土壤侵蝕關(guān)系的對比研究[J]. 土壤學報,2010,47(6):1060-1066.

    Sun Jiajia, Yu Dongsheng, Shi Xuezheng, et al. Comparison of between LAI and VFC in relationship with soil erosion in the red soil hilly region of south China[J]. Acta Pedologica Sinica, 2010, 47(6): 1060-1066. (in Chinese with English abstract)

    [16] Wen Z M, Lees B G, Feng J, et al. Stratified vegetation cover index: A new way to assess vegetation impact on soil erosion[J]. Catena, 2010, 83(1): 87-93.

    [17] Yu Dong Sheng, Shi Xue Zheng, Weindorf D C. Relationships between permeability and erodibility of cultivated Acrisols and Cambisols in subtropical China[J]. Pedosphere, 2006, 16(3): 304-311.

    [18] Zhang W T, Yu Dong Sheng, Shi Xue Zheng, et al. Variation of sediment concentration and its drivers under different soil management systems[J]. Pedosphere, 2010, 20(5): 578-585.

    [19] 張黎明,于東升,史學正,等. 不同降雨類型下的南方典型土壤侵蝕量差異研究[J]. 水土保持通報,2011,31(4):1-6.

    Zhang Liming, Yu Dongsheng, Shi Xuezheng, et al.Soil loss of different typical soils in south china under different rainfall types[J]. Bulletin of Soil and Water Conservation, 2011, 31(4): 1-6. (in Chinese with English abstract)

    [20] Kohonen T. The self-organizing map[J]. Proceedings of the IEEE, 1990, 78(9): 1464-1480.

    [21] Huang Z, Ouyang Z, Li F, et al. Response of runoff and soil loss to reforestation and rainfall type in red soil region of southern China[J]. Journal of Environmental Sciences, 2010, 22(11): 1765-1773.

    [22] 秦偉,左長清,晏清洪,等. 紅壤裸露坡地次降雨土壤侵蝕規(guī)律[J]. 農(nóng)業(yè)工程學報,2015,31(2):124-132.

    Qin Wei, Zuo Changqing, Yan Qinghong, et al. Regularity of individual rainfall soil erosion in bare slope land of red soil[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(2): 124-132. (in Chinese with English abstract)

    [23] Liu Y J, Yang J, Hu J M, et al. Characteristics of the surface-subsurface flow generation and sediment yield to the rainfall regime and land-cover by long-term in-situ observation in the red soil region, Southern China[J]. Journal of Hydrology, 2016, 539: 457-467.

    [24] Zheng F L. Effect of vegetation changes on soil erosion on the Loess Plateau[J]. Pedosphere, 2006, 16(4): 420-427.

    [25] 張海東,于東升,史學正,等. 地面三維激光掃描方法研究紅壤侵蝕恢復區(qū)植被垂直結(jié)構(gòu)[J]. 土壤,2014,46(5):948-953.

    Zhang Haidong, Yu Dongsheng, Shi xuezheng, et al.Retrieval of forest vertical structure in eroded area of southern china using terrestrial LiDAR systems[J]. Soil, 2014, 46(5): 948-953. (in Chinese with English abstract)

    [26] Xu C, Yang Z, Qian W, et al. Runoff and soil erosion responses to rainfall and vegetation cover under various afforestation management regimes in subtropical montane forest[J]. Land Degradation & Development, 2019, 30(14): 1711-1724.

    [27] Llorens P, Poch R, Latron J, et al. Rainfall interception by a Pinus sylvestris forest patch overgrown in a Mediterranean mountainous abandoned area I. Monitoring design and results down to the event scale[J]. Journal of Hydrology, 1997, 199(3/4): 331-345.

    [28] Gei?ler C, Lang A C, Von Oheimb G, et al. Impact of tree saplings on the kinetic energy of rainfall: The importance of stand density, species identity and tree architecture in subtropical forests in China[J]. Agricultural and Forest Meteorology, 2012, 156: 31-40.

    [29] Nagase A, Dunnett N. Amount of water runoff from different vegetation types on extensive green roofs: Effects of plant species, diversity and plant structure[J]. Landscape and Urban Planning, 2012, 104(3/4): 356-363.

    [30] 鄔鈴莉,王云琦,王晨灃,等. 降雨類型對北方土石山區(qū)坡面土壤侵蝕的影響[J]. 農(nóng)業(yè)工程學報,2017,33(24):157-164.

    Wu Lingli, Wang Yunqi, Wang Chenfeng, et al. Effect of rainfall patterns on hillslope soil erosion in rocky mountain area of north China[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(24): 157-164. (in Chinese with English abstract)

    [31] Brasil J B, de Andrade E M, de Queiroz Palácio H A, et al. Characteristics of precipitation and the process of interception in a seasonally dry tropical forest[J]. Journal of Hydrology: Regional Studies, 2018, 19: 307-317.

    [32] Bagio B, Bertol I, Wolschick N H, et al. Water erosion in different slope lengths on bare soil[J]. Revista Brasileira de Ciência do Solo, 2017, 41:e0160132.

    [33] Fang N F, Shi Z H, Li L, et al. The effects of rainfall regimes and land use changes on runoff and soil loss in a small mountainous watershed[J]. Catena, 2012, 99: 1-8.

    [34] Wei W, Chen L, Fu B, et al. The effect of land uses and rainfall regimes on runoff and soil erosion in the semi-arid loess hilly area, China[J]. Journal of Hydrology, 2007, 335(3/4): 247-258.

    [35] Li G, Wan L, Cui M, et al. Influence of canopy interception and rainfall kinetic energy on soil erosion under forests[J]. Forests, 2019, 10(6): 509-523.

    [36] Mishra A K, Placzek C, Jones R. Coupled influence of precipitation and vegetation on millennial-scale erosion rates derived from 10Be[J]. PloS one, 2019, 14(1): e0211325.

    [37] Yu Y, Wei W, Chen L, et al. Quantifying the effects of precipitation, vegetation, and land preparation techniques on runoff and soil erosion in a Loess watershed of China[J]. Science of the Total Environment, 2019, 652: 755-764.

    [38] Zhou J, Fu B, Gao G, et al. Effects of precipitation and restoration vegetation on soil erosion in a semi-arid environment in the Loess Plateau, China[J]. Catena, 2016, 137: 1-11.

    Effects of vegetation structure types and rainfall patterns on soil and water loss of understory vegetation in red soil areas of South China

    Chen Yang1,3, Zhang Haidong2, Yu Dongsheng1,3※, Ma Lixia1, Lu Xiaosong1,3, Wang Xin1,3

    (1.,,210008,; 2.(),215155,; 3.,100049,)

    Vegetation and rainfall are the key factors that affect the soil and soil erosion. Exploring the surface runoff and soil loss process associated with different rainfall patterns and vegetation structure types is of great significance to the development of soil and water conservation. In this study, we carried out a series of experiments on surface runoffs and sediments monitoring at different runoff plots with five structures types of vegetation (i.e., shrub-grass mixed forest, grassland, low shrub, arbor forest, high shrub). There are 8 runoff test plots with length × width (10 m×2.6 m) on the gentle slope of 8°in Red soil Ecological Experimental Station, Yingtan city. The experiment lasted three years from 2016 to 2018. During the study period, 93 individual natural erosive rainfall events, surface runoffs and sediments data were collected, as well as the vegetation structure parameters were obtained. The self-organizing mapping (SOM) method was applied to classify rainfall patterns according to rainfall characteristic indicators such as rainfall amounts, rainfall durations, 60min maximum rainfall intensity, average rainfall intensity and rainfall amount concentration index. Then the characteristics of soil and water loss under different rainfall patterns and vegetation structure types were studied, and the effect of rainfall and vegetation on soil and water loss were quantitative analysis by redundancy analysis (RDA). The results showed that SOM method could identify the four typical erosion rainfall patterns in red soil region, the rainfall patterns were as follows: RI(low intensity, long duration, minor amount of rainfall and low rainfall concentration index), RII(medium amount of rainfall, moderate intensity and long duration), RIII(short duration, strong rainfall intensity and high rainfall concentration index), RIV(a great amount of rainfall, strong rainfall intensity and long duration). Among all of the erosion rainfall events, the least rainfall was pattern RIV, which was the most erosive and destructive, so it was easy to cause serious soil erosion. RIIIwas the main pattern that made the highest cumulative contribution rate to soil and water loss. The erosion ability of pattern RIIwas moderate, but it was also easy to cause soil and water loss due to its moderate intensity and long duration. The pattern RIcould hardly cause soil erosion. Vegetation structure types significantly affected the soil erosion and runoff generation, and the soil and water conservation function of different vegetation in turn was as followed the order of shrub-grass mixed forest, grassland, low shrub, arbor forest, and high shrub. RDA results shown that comprehensive interaction effects on soil erosion exist between rainfall patterns and vegetation structure types. The ability of soil and water conservation of vegetation decreased gradually and the effects of rainfall became greater in the later stage of erosion process, when rainfall patterns changed from weak to strong (RI?RII?RIII?RIV), meanwhile the dominant factors of soil and water loss changed from vegetation to the interaction of vegetation and precipitation, and then to rainfalls. The influences of canopy structure on soil and water loss reduced, while the ones of rainfall increased when vegetation structure types changed from near surface vegetation to shrub and tree.

    erosion; runoff; soils; vegetation structure; rainfall pattern; red soil region in South China

    2019-11-11

    2019-12-12

    國家自然科學基金項目(41571206,41501296);土壤與農(nóng)業(yè)可持續(xù)發(fā)展國家重點實驗室開放研究基金項目(Y812000002)

    陳 洋,博士生,主要從事資源環(huán)境遙感應用與水土流失研究。Email:chenyang@issas.ac.cn.

    于東升,研究員,主要從事水土保持與資源環(huán)境遙感應用研究。Email:dshyu@issas.ac.cn

    10.11975/j.issn.1002-6819.2020.05.017

    S157.1

    A

    1002-6819(2020)-05-0150-08

    陳 洋,張海東,于東升,馬利霞,陸曉松,王 鑫. 南方紅壤區(qū)植被結(jié)構(gòu)類型與降雨模式對林下水土流失的影響[J]. 農(nóng)業(yè)工程學報,2020,36(5):150-157. doi:10.11975/j.issn.1002-6819.2020.05.017 http://www.tcsae.org

    Chen Yang, Zhang Haidong, Yu Dongsheng, Ma Lixia, Lu Xiaosong, Wang Xin. Effects of vegetation structure types and rainfall patterns on soil and water loss of understory vegetation in red soil areas of South China[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(5): 150-157. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2020.05.017 http://www.tcsae.org

    猜你喜歡
    紅壤土壤侵蝕徑流
    初探熱脫附技術(shù)在有機污染紅壤修復的應用
    昆鋼科技(2021年6期)2021-03-09 06:10:26
    鄉(xiāng)村聚落土壤侵蝕環(huán)境與水土流失研究綜述
    海壇島土壤侵蝕問題研究
    長期不同施肥對赤紅壤稻田區(qū)肥力的影響
    Topmodel在布哈河流域徑流模擬中的應用
    大別山區(qū)土壤侵蝕動態(tài)變化及趨勢預測
    探秘“大徑流”
    攻克“大徑流”
    南水北調(diào)中線水源地土壤侵蝕經(jīng)濟損失估算
    江埡水庫降雨徑流相關(guān)圖的建立
    蜜桃久久精品国产亚洲av| 久久精品国产清高在天天线| 日本五十路高清| 麻豆成人av在线观看| 好看av亚洲va欧美ⅴa在| 日日摸夜夜添夜夜添小说| 高清在线国产一区| 高清在线国产一区| 国产高清videossex| 免费看日本二区| 国产69精品久久久久777片| 精品一区二区三区视频在线 | 国产aⅴ精品一区二区三区波| 高清在线国产一区| 99热这里只有是精品50| 欧美成人免费av一区二区三区| 欧美乱妇无乱码| 久久亚洲真实| 色播亚洲综合网| 欧美在线一区亚洲| 国内毛片毛片毛片毛片毛片| 亚洲熟妇熟女久久| 有码 亚洲区| 久久久久久国产a免费观看| 最近视频中文字幕2019在线8| 午夜福利免费观看在线| 精品久久久久久久毛片微露脸| 女生性感内裤真人,穿戴方法视频| 19禁男女啪啪无遮挡网站| 亚洲熟妇熟女久久| 久久久久久久午夜电影| 国产欧美日韩一区二区三| 日韩 欧美 亚洲 中文字幕| 国产一区在线观看成人免费| 成人三级黄色视频| 好男人电影高清在线观看| 又粗又爽又猛毛片免费看| av中文乱码字幕在线| 久久国产乱子伦精品免费另类| 18美女黄网站色大片免费观看| 超碰av人人做人人爽久久 | а√天堂www在线а√下载| 亚洲,欧美精品.| 国产成人福利小说| 淫秽高清视频在线观看| 深爱激情五月婷婷| 观看美女的网站| 丝袜美腿在线中文| 色播亚洲综合网| 男插女下体视频免费在线播放| 一边摸一边抽搐一进一小说| 色精品久久人妻99蜜桃| 亚洲精品久久国产高清桃花| 久久久久九九精品影院| 成人性生交大片免费视频hd| 国产精品嫩草影院av在线观看 | 精品99又大又爽又粗少妇毛片 | 国产高清视频在线观看网站| 午夜福利欧美成人| 久久婷婷人人爽人人干人人爱| 久久99热这里只有精品18| 国内揄拍国产精品人妻在线| 久久婷婷人人爽人人干人人爱| 免费av毛片视频| 婷婷六月久久综合丁香| 国产成年人精品一区二区| 99热只有精品国产| 国产高清视频在线观看网站| 少妇的逼好多水| 亚洲熟妇熟女久久| 久久精品综合一区二区三区| 在线观看美女被高潮喷水网站 | 亚洲美女视频黄频| 久久久成人免费电影| 精品乱码久久久久久99久播| 免费电影在线观看免费观看| 亚洲熟妇中文字幕五十中出| 麻豆国产97在线/欧美| 久久久精品欧美日韩精品| 日韩欧美免费精品| 亚洲av熟女| 亚洲七黄色美女视频| 日韩高清综合在线| 毛片女人毛片| 狂野欧美白嫩少妇大欣赏| 免费在线观看亚洲国产| 国产精品久久久久久久久免 | 亚洲精品国产精品久久久不卡| 午夜亚洲福利在线播放| 欧美成人一区二区免费高清观看| 一本一本综合久久| 夜夜夜夜夜久久久久| 国产精品 欧美亚洲| 九色国产91popny在线| 香蕉久久夜色| 激情在线观看视频在线高清| 老汉色∧v一级毛片| 好男人电影高清在线观看| 少妇丰满av| 精品乱码久久久久久99久播| 综合色av麻豆| 久久天躁狠狠躁夜夜2o2o| 精品欧美国产一区二区三| 99在线视频只有这里精品首页| 精品熟女少妇八av免费久了| 最新在线观看一区二区三区| 亚洲国产精品久久男人天堂| 天堂网av新在线| 无限看片的www在线观看| 国内精品一区二区在线观看| 久久久久久人人人人人| 日本免费一区二区三区高清不卡| 中文字幕av成人在线电影| 亚洲精品日韩av片在线观看 | 18美女黄网站色大片免费观看| 久久久国产成人免费| 精品99又大又爽又粗少妇毛片 | 亚洲av日韩精品久久久久久密| 叶爱在线成人免费视频播放| 久久精品91无色码中文字幕| 熟女人妻精品中文字幕| 熟女少妇亚洲综合色aaa.| 夜夜躁狠狠躁天天躁| 亚洲av中文字字幕乱码综合| 搞女人的毛片| 国产精品久久久久久久电影 | 国产日本99.免费观看| 法律面前人人平等表现在哪些方面| 国产精品野战在线观看| 搞女人的毛片| 动漫黄色视频在线观看| 国产野战对白在线观看| 99精品欧美一区二区三区四区| 19禁男女啪啪无遮挡网站| 国产在视频线在精品| 亚洲avbb在线观看| 窝窝影院91人妻| 天堂网av新在线| 午夜激情福利司机影院| 欧美一级毛片孕妇| 99热这里只有是精品50| 国产亚洲精品av在线| 91九色精品人成在线观看| h日本视频在线播放| 日韩精品青青久久久久久| 国产精品久久久人人做人人爽| 国产精品三级大全| 2021天堂中文幕一二区在线观| 国内精品美女久久久久久| 国产一区在线观看成人免费| 色综合站精品国产| 桃红色精品国产亚洲av| 精品午夜福利视频在线观看一区| 亚洲欧美精品综合久久99| 国产 一区 欧美 日韩| 午夜日韩欧美国产| 亚洲中文字幕日韩| 亚洲国产精品sss在线观看| 一二三四社区在线视频社区8| 久久久精品欧美日韩精品| 亚洲熟妇中文字幕五十中出| 露出奶头的视频| 久久精品综合一区二区三区| 性色avwww在线观看| 久久草成人影院| 一级毛片女人18水好多| 少妇高潮的动态图| av天堂在线播放| 五月玫瑰六月丁香| 亚洲黑人精品在线| 欧美日韩一级在线毛片| av黄色大香蕉| 久久精品国产亚洲av涩爱 | 亚洲精华国产精华精| 女人被狂操c到高潮| 欧美最黄视频在线播放免费| 在线观看66精品国产| 午夜免费成人在线视频| 亚洲熟妇中文字幕五十中出| 久久精品91无色码中文字幕| 两人在一起打扑克的视频| 啦啦啦韩国在线观看视频| 波野结衣二区三区在线 | 中文字幕人妻丝袜一区二区| 国产高清有码在线观看视频| 美女 人体艺术 gogo| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 国产精品久久久久久人妻精品电影| 亚洲无线在线观看| 18禁国产床啪视频网站| 母亲3免费完整高清在线观看| 黄片大片在线免费观看| 国产精品爽爽va在线观看网站| 老汉色∧v一级毛片| 色综合欧美亚洲国产小说| 久久香蕉国产精品| 中文在线观看免费www的网站| 久久精品影院6| 超碰av人人做人人爽久久 | 国产aⅴ精品一区二区三区波| 久久久久国产精品人妻aⅴ院| 最好的美女福利视频网| 亚洲内射少妇av| 91在线精品国自产拍蜜月 | 1000部很黄的大片| 国产主播在线观看一区二区| 听说在线观看完整版免费高清| 亚洲欧美日韩高清专用| 欧美国产日韩亚洲一区| 欧美+日韩+精品| 国产爱豆传媒在线观看| 黄色女人牲交| 午夜亚洲福利在线播放| 亚洲av不卡在线观看| 国产高清三级在线| 变态另类成人亚洲欧美熟女| 美女免费视频网站| 日本免费a在线| 又爽又黄无遮挡网站| 欧美最黄视频在线播放免费| 国产一级毛片七仙女欲春2| 国产伦精品一区二区三区视频9 | 三级毛片av免费| 免费观看精品视频网站| 成人特级黄色片久久久久久久| 国产主播在线观看一区二区| 男女之事视频高清在线观看| 午夜激情福利司机影院| 99热这里只有是精品50| 亚洲电影在线观看av| 久久精品国产清高在天天线| 精品久久久久久久末码| www国产在线视频色| 欧美乱色亚洲激情| 99国产极品粉嫩在线观看| 亚洲狠狠婷婷综合久久图片| 成人av在线播放网站| 久久久久久久亚洲中文字幕 | 午夜两性在线视频| 日韩欧美一区二区三区在线观看| 欧美成人性av电影在线观看| 免费看日本二区| 好看av亚洲va欧美ⅴa在| 久久久久久九九精品二区国产| 亚洲成人中文字幕在线播放| 香蕉丝袜av| www国产在线视频色| 国产午夜精品论理片| 国产三级在线视频| 极品教师在线免费播放| 国产中年淑女户外野战色| 一区二区三区国产精品乱码| 国内精品美女久久久久久| 国产精品一区二区三区四区久久| 18美女黄网站色大片免费观看| 色尼玛亚洲综合影院| 欧美日韩精品网址| 国产精品香港三级国产av潘金莲| 一夜夜www| 国产97色在线日韩免费| 欧美高清成人免费视频www| 亚洲精品亚洲一区二区| 久久久久久久午夜电影| 成人性生交大片免费视频hd| 精品久久久久久久人妻蜜臀av| 精品一区二区三区人妻视频| 99精品欧美一区二区三区四区| 国产v大片淫在线免费观看| 老司机福利观看| 日韩国内少妇激情av| 国产成人系列免费观看| 两个人视频免费观看高清| 欧美另类亚洲清纯唯美| 亚洲最大成人中文| 18禁黄网站禁片免费观看直播| 制服丝袜大香蕉在线| 国产单亲对白刺激| 日韩国内少妇激情av| av天堂在线播放| 一个人免费在线观看电影| 亚洲精品456在线播放app | 毛片女人毛片| 天堂网av新在线| 国产视频一区二区在线看| 在线免费观看的www视频| 亚洲av免费高清在线观看| 精品久久久久久久久久免费视频| 我的老师免费观看完整版| 变态另类丝袜制服| 女警被强在线播放| 国产欧美日韩精品一区二区| 亚洲av成人不卡在线观看播放网| 噜噜噜噜噜久久久久久91| 桃色一区二区三区在线观看| 亚洲欧美日韩东京热| 又粗又爽又猛毛片免费看| a级毛片a级免费在线| 非洲黑人性xxxx精品又粗又长| 亚洲一区二区三区不卡视频| 精华霜和精华液先用哪个| 亚洲黑人精品在线| 国内久久婷婷六月综合欲色啪| 一区二区三区激情视频| 一二三四社区在线视频社区8| 成人性生交大片免费视频hd| 亚洲国产日韩欧美精品在线观看 | 18禁在线播放成人免费| 性欧美人与动物交配| 色综合婷婷激情| 国产三级在线视频| 成人午夜高清在线视频| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 色尼玛亚洲综合影院| 亚洲 欧美 日韩 在线 免费| 亚洲熟妇熟女久久| 午夜精品久久久久久毛片777| 久久久久久久久中文| 丁香欧美五月| 母亲3免费完整高清在线观看| 亚洲欧美日韩高清在线视频| 又粗又爽又猛毛片免费看| 亚洲av第一区精品v没综合| 亚洲天堂国产精品一区在线| 免费在线观看影片大全网站| 亚洲无线在线观看| 蜜桃亚洲精品一区二区三区| av福利片在线观看| 国产av不卡久久| 国产麻豆成人av免费视频| 99久久99久久久精品蜜桃| av中文乱码字幕在线| 欧美在线一区亚洲| 精品人妻一区二区三区麻豆 | 在线观看66精品国产| 亚洲国产欧洲综合997久久,| 成年女人看的毛片在线观看| 亚洲无线在线观看| 内地一区二区视频在线| 久久精品91蜜桃| 国产精品女同一区二区软件 | 中文字幕熟女人妻在线| 亚洲精品乱码久久久v下载方式 | 国产真实乱freesex| 免费看光身美女| 国产精品一及| 久久草成人影院| 午夜福利高清视频| 国产中年淑女户外野战色| 精品人妻偷拍中文字幕| 亚洲欧美日韩高清专用| 国产午夜福利久久久久久| 久久久久九九精品影院| 国产99白浆流出| 免费观看的影片在线观看| 全区人妻精品视频| tocl精华| 我的老师免费观看完整版| 18禁裸乳无遮挡免费网站照片| 亚洲精品一卡2卡三卡4卡5卡| 一个人观看的视频www高清免费观看| 国产精品乱码一区二三区的特点| 丝袜美腿在线中文| 一级黄色大片毛片| www.www免费av| 国产精品亚洲一级av第二区| 一本久久中文字幕| 三级国产精品欧美在线观看| 日韩有码中文字幕| 97碰自拍视频| 首页视频小说图片口味搜索| av国产免费在线观看| 国产精品久久久久久久久免 | 亚洲av美国av| 国产探花在线观看一区二区| 国产一区二区激情短视频| 国产精品野战在线观看| 嫩草影院精品99| 亚洲国产日韩欧美精品在线观看 | 亚洲专区国产一区二区| 性色avwww在线观看| 在线天堂最新版资源| 日本 欧美在线| 国语自产精品视频在线第100页| 男女视频在线观看网站免费| 在线天堂最新版资源| 国产亚洲av嫩草精品影院| 日本黄色视频三级网站网址| 别揉我奶头~嗯~啊~动态视频| 12—13女人毛片做爰片一| 欧美色视频一区免费| 国内精品美女久久久久久| 亚洲欧美日韩高清专用| 嫩草影视91久久| 国产日本99.免费观看| 国产野战对白在线观看| 成人亚洲精品av一区二区| 亚洲片人在线观看| 天堂影院成人在线观看| 男女床上黄色一级片免费看| 窝窝影院91人妻| 丰满乱子伦码专区| 国产午夜精品论理片| 亚洲欧美一区二区三区黑人| 精品福利观看| 美女cb高潮喷水在线观看| 精品国产亚洲在线| 一进一出抽搐动态| 午夜福利18| 少妇的丰满在线观看| 一a级毛片在线观看| 欧美国产日韩亚洲一区| 国产高清视频在线观看网站| 国产av麻豆久久久久久久| 九九热线精品视视频播放| 欧美乱码精品一区二区三区| 日日干狠狠操夜夜爽| 亚洲色图av天堂| 一边摸一边抽搐一进一小说| 亚洲一区二区三区色噜噜| 亚洲人成电影免费在线| 亚洲精品在线观看二区| 99久久综合精品五月天人人| 精品国产美女av久久久久小说| 男女那种视频在线观看| 国产成人福利小说| netflix在线观看网站| 中文字幕熟女人妻在线| 国产毛片a区久久久久| 国产精品久久视频播放| av福利片在线观看| 69av精品久久久久久| 国产精品99久久久久久久久| 一个人观看的视频www高清免费观看| 国产毛片a区久久久久| 欧美日韩乱码在线| 看黄色毛片网站| 久久精品影院6| 国产主播在线观看一区二区| 久久久久久久久大av| 亚洲精品在线美女| 中国美女看黄片| 一个人看的www免费观看视频| 听说在线观看完整版免费高清| 少妇的逼水好多| 女生性感内裤真人,穿戴方法视频| 中文字幕久久专区| 18禁黄网站禁片免费观看直播| 国产单亲对白刺激| 精品久久久久久,| 很黄的视频免费| 精品电影一区二区在线| 日韩亚洲欧美综合| 熟妇人妻久久中文字幕3abv| 久久久久久久午夜电影| 午夜视频国产福利| 国产精品日韩av在线免费观看| 看免费av毛片| 男女床上黄色一级片免费看| 国产综合懂色| 亚洲乱码一区二区免费版| 国产淫片久久久久久久久 | 性欧美人与动物交配| 嫩草影视91久久| 午夜免费激情av| 色老头精品视频在线观看| 中文字幕熟女人妻在线| 丁香六月欧美| 波多野结衣高清作品| 亚洲第一电影网av| 国产真人三级小视频在线观看| 桃色一区二区三区在线观看| 久久精品人妻少妇| 嫁个100分男人电影在线观看| 一个人看的www免费观看视频| 91在线精品国自产拍蜜月 | av在线天堂中文字幕| 国产精品野战在线观看| 99在线人妻在线中文字幕| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 激情在线观看视频在线高清| 一边摸一边抽搐一进一小说| 又粗又爽又猛毛片免费看| 91麻豆av在线| 波多野结衣高清作品| 亚洲一区二区三区不卡视频| 99riav亚洲国产免费| 久久精品人妻少妇| 内地一区二区视频在线| 有码 亚洲区| 一个人看视频在线观看www免费 | 国产野战对白在线观看| 级片在线观看| 嫩草影院精品99| 国产精品99久久99久久久不卡| 国产淫片久久久久久久久 | 精品乱码久久久久久99久播| 中文字幕av在线有码专区| 精品免费久久久久久久清纯| 全区人妻精品视频| e午夜精品久久久久久久| 国产精品久久久久久精品电影| 精品久久久久久久人妻蜜臀av| 51午夜福利影视在线观看| 黄色成人免费大全| 欧美日韩一级在线毛片| 国产一区在线观看成人免费| 高潮久久久久久久久久久不卡| 99国产综合亚洲精品| 亚洲精品一卡2卡三卡4卡5卡| 久99久视频精品免费| 啦啦啦免费观看视频1| 亚洲在线观看片| 免费观看人在逋| 亚洲国产高清在线一区二区三| aaaaa片日本免费| 搡女人真爽免费视频火全软件 | 中文字幕人妻丝袜一区二区| 特级一级黄色大片| 亚洲av免费高清在线观看| 在线视频色国产色| 夜夜躁狠狠躁天天躁| 国产黄色小视频在线观看| 亚洲国产高清在线一区二区三| 两人在一起打扑克的视频| 午夜免费观看网址| 国产亚洲精品一区二区www| 国产成人影院久久av| av中文乱码字幕在线| 国产在视频线在精品| 午夜福利免费观看在线| 又紧又爽又黄一区二区| 美女 人体艺术 gogo| 99久久精品热视频| 欧美高清成人免费视频www| 18美女黄网站色大片免费观看| 波多野结衣高清无吗| 国产精品永久免费网站| 亚洲av电影在线进入| 精品不卡国产一区二区三区| 老司机午夜十八禁免费视频| 国产精品久久久人人做人人爽| 99久久99久久久精品蜜桃| www.www免费av| 日本与韩国留学比较| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 最近最新中文字幕大全免费视频| 精品一区二区三区人妻视频| 国产v大片淫在线免费观看| 日本精品一区二区三区蜜桃| 欧美日韩瑟瑟在线播放| av女优亚洲男人天堂| 少妇丰满av| 小蜜桃在线观看免费完整版高清| 欧美成人a在线观看| 亚洲精品在线观看二区| 久久国产精品人妻蜜桃| 免费观看人在逋| 老汉色∧v一级毛片| 一卡2卡三卡四卡精品乱码亚洲| 变态另类丝袜制服| 国产精品1区2区在线观看.| 亚洲欧美激情综合另类| 高清日韩中文字幕在线| 国产免费男女视频| tocl精华| 一区福利在线观看| 午夜视频国产福利| 18禁裸乳无遮挡免费网站照片| 久久国产乱子伦精品免费另类| 国产在线精品亚洲第一网站| 一夜夜www| 欧美一区二区精品小视频在线| 亚洲国产中文字幕在线视频| 午夜福利免费观看在线| 全区人妻精品视频| 99久国产av精品| 久久精品综合一区二区三区| 在线a可以看的网站| 日本撒尿小便嘘嘘汇集6| 国产精品,欧美在线| 久久精品国产亚洲av涩爱 | 国产精品电影一区二区三区| 最近最新免费中文字幕在线| 久久久久亚洲av毛片大全| 亚洲人成网站高清观看| 99国产精品一区二区蜜桃av| 国产精品久久久久久亚洲av鲁大| 男女那种视频在线观看| 51午夜福利影视在线观看| 每晚都被弄得嗷嗷叫到高潮| 色尼玛亚洲综合影院| 超碰av人人做人人爽久久 | 男人的好看免费观看在线视频| 色在线成人网| 啦啦啦观看免费观看视频高清| 色播亚洲综合网| 波多野结衣高清无吗| 精品一区二区三区av网在线观看| 欧美zozozo另类| 精品人妻1区二区| 国产精品一区二区免费欧美| 欧美在线一区亚洲| av中文乱码字幕在线| 免费观看精品视频网站| 国模一区二区三区四区视频| 欧美av亚洲av综合av国产av| 成人一区二区视频在线观看| 一边摸一边抽搐一进一小说| av国产免费在线观看| 天堂av国产一区二区熟女人妻| 男人的好看免费观看在线视频| 午夜福利18| 夜夜躁狠狠躁天天躁| 亚洲五月天丁香| 亚洲精品日韩av片在线观看 |