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

    基于最小數(shù)據(jù)集的祁連山南坡不同土地利用方式土壤質(zhì)量評(píng)價(jià)

    2024-12-31 00:00:00邱巡巡曹廣超趙青林曹生奎趙美亮何啟欣白嘉奇咸慶玲
    草地學(xué)報(bào) 2024年9期

    摘要:為了評(píng)估高寒地區(qū)不同土地利用方式下的土壤質(zhì)量狀況,明確不同土地利用方式下的土壤質(zhì)量的關(guān)鍵影響因子,在祁連山南坡采集了林地、灌叢、草地及耕地4種主要土地利用類(lèi)型的土壤樣品174份,通過(guò)主成分分析(Principal component analysis,PCA)建立最小數(shù)據(jù)集(Minimum data set,MDS),綜合評(píng)估研究區(qū)不同土地利用方式下的土壤質(zhì)量。結(jié)果表明:林地、灌叢、草地和耕地土壤質(zhì)量指數(shù)值分別為0.535,0.519,0.466和0.544,表現(xiàn)為耕地>林地>灌叢>草地,對(duì)土壤質(zhì)量分級(jí)為Ⅰ~Ⅵ級(jí),對(duì)應(yīng)指數(shù)分別為≤0.3,(0.3~0.4],(0.4~0.5],(0.5~0.6],(0.6,0.7]和gt;0.7,草地等級(jí)為Ⅲ級(jí),處于“中等”水平;耕地、林地和灌叢土壤質(zhì)量等級(jí)為Ⅳ級(jí),處于“中等偏上”水平。土壤質(zhì)量關(guān)鍵指標(biāo)間存在互相影響,因此,建議研究區(qū)域土地要實(shí)施分類(lèi)科學(xué)管理。此外,合理開(kāi)發(fā)和應(yīng)用綠色高效的新型生物技術(shù)是應(yīng)對(duì)影響研究區(qū)土壤質(zhì)量的微生物指標(biāo)的有效措施。

    關(guān)鍵詞:最小數(shù)據(jù)集;土壤質(zhì)量;祁連山南坡

    中圖分類(lèi)號(hào):S158""" 文獻(xiàn)標(biāo)識(shí)碼:A""""" 文章編號(hào):1007-0435(2024)09-2952-10

    收稿日期:2024-03-27;修回日期:2024-04-11

    基金項(xiàng)目:青海省自然科學(xué)基金項(xiàng)目“祁連山南麓典型流域關(guān)鍵帶水文過(guò)程及碳收支變化特征研究”(2023-ZJ-907M)資助

    作者簡(jiǎn)介:

    邱巡巡(1997-),女,漢族,貴州遵義人,博士研究生,主要從事自然地理與地表環(huán)境過(guò)程研究,E-mail:1260469240@qq.com;*通信作者Author for correspondence,E-mail:caoguangchao@126.com

    doi:10.11733/j.issn.1007-0435.2024.09.029

    引用格式:

    邱巡巡, 曹廣超, 趙青林,等.基于最小數(shù)據(jù)集的祁連山南坡不同土地利用方式土壤質(zhì)量評(píng)價(jià)[J].草地學(xué)報(bào),2024,32(9):2952-2961

    QIU Xunxun, CAO Guangchao, ZHAO Qinglin,et al.Assessment of Soil Quality under Different Land Use Practices on the Southern Slope of Qilian Mountains Based on Minimum Data Set[J].Acta Agrestia Sinica,2024,32(9):2952-2961

    Assessment of Soil Quality under Different Land Use Practices on the

    Southern Slope of Qilian Mountains Based on Minimum Data Set

    QIU Xun-xun1,2,3, CAO Guang-chao2,3*, ZHAO Qing-lin1,4, CAO Sheng-kui1,2,

    ZHAO Mei-liang1,2,3,HE Qi-xin1,2,3, BAI Jia-qi1,2,3, XIAN Qing-ling1,2,3

    (1.College of Geographical Sciences, Qinghai Normal University, Xining, Qinghai Province 810008, China; 2.Qinghai Provincial

    Key Laboratory of Physical Geography and Environmental Process, Xining, Qinghai Province 810008, China; 3.Key Laboratory of

    Surface Processes and Ecological Conservation on Tibetan Plateau, Ministry of Education, Xining, Qinghai Province 810008, China;

    4.Graduate School of Qinghai Normal University, Xining, Qinghai Province 810008, China)

    Abstract:To assess the soil quality under different land use practices in the high-cold region and identify the key influencing factors on soil quality under different land use types,174 soil samples from four main land use practices,including forests,shrubs,grasslands,and cultivated lands,were collected on the southern slope of the Qilian Mountains. A minimum data set (MDS) was established using principal component analysis to comprehensively evaluate soil quality under different land use practices in the study area. The results showed that the soil quality indices for forests,shrubs,grasslands,and cultivated lands were 0.535,0.519,0.466,and 0.544,respectively,with the order of cultivated landsgt;forestsgt;shrubsgt;grasslands. Corresponding to the soil quality classification,the soil quality levels for grasslands were classified as Grade III,indicating a “medium” level,while the soil quality level for cultivated lands,forests,and" shrubs was classified as Grade IV,indicating a “above-average” level. There were interactions among the key indicators of soil quality. Scientific classification management should be implemented for land in the research area. Regarding to the microorganisms affecting soil quality,the rational development and application of environmentally friendly,green and efficient new biotechnologies are effective measures to address the microbial indicators influencing soil quality in the study area.

    Key words:Soil quality;Minimum data set;South slope of Qilian Mountains

    土壤質(zhì)量是指特定種類(lèi)的土壤在自然或管理的生態(tài)系統(tǒng)范圍內(nèi)維持植物和動(dòng)物生產(chǎn)力、保持或提高水和空氣質(zhì)量以及支持人類(lèi)健康和居住的能力[1-2],土壤質(zhì)量的好壞取決于土地利用方式、生態(tài)系統(tǒng)類(lèi)型以及土壤內(nèi)部各種特征的相互作用[3]。其中,土地利用通過(guò)干擾土壤的理化以及生物性質(zhì),導(dǎo)致土壤肥力等發(fā)生變化,進(jìn)而影響土壤質(zhì)量[5-7]。良好的土壤質(zhì)量不僅具有較高的生產(chǎn)力,而且有助于區(qū)域水土生態(tài)環(huán)境的改善[4]。因此,開(kāi)展不同土地利用方式土壤質(zhì)量評(píng)價(jià)工作,對(duì)土地管理模式優(yōu)化和結(jié)構(gòu)調(diào)整具有重要意義[8]。

    土壤理化和生物特性等功能因子在土壤質(zhì)量評(píng)價(jià)中往往被作為相關(guān)指標(biāo)[9]。評(píng)價(jià)結(jié)果的準(zhǔn)確與否和土壤質(zhì)量評(píng)價(jià)中參數(shù)的嚴(yán)格選擇息息相關(guān)[10],且選取的大量指標(biāo)之間往往存在相關(guān)關(guān)系,會(huì)造成數(shù)據(jù)間的冗余,準(zhǔn)確提取適宜評(píng)價(jià)指標(biāo)是土壤質(zhì)量評(píng)價(jià)的重要環(huán)節(jié)。Larson等[11]提出采用土壤參數(shù)最小數(shù)據(jù)集(Minimum data set,MDS)來(lái)鑒別評(píng)價(jià)指標(biāo)間的相互關(guān)系以及反映對(duì)土壤特性和作物的影響。運(yùn)用主成分分析(Principal component analysis,PCA)從大量土壤參數(shù)中篩選出相對(duì)獨(dú)立、影響土壤質(zhì)量的敏感性指標(biāo)建立MDS,目前廣泛應(yīng)用于評(píng)估土壤質(zhì)量[12-13],涉及土地利用變化、草地不同放牧率下、森林不同治理模式下的土壤質(zhì)量評(píng)價(jià)等[14-16]。土地利用方式對(duì)土壤質(zhì)量產(chǎn)生顯著的影響,開(kāi)展不同土地利用類(lèi)型高寒土壤質(zhì)量調(diào)查,建立高寒土地管理評(píng)估框架體系對(duì)區(qū)域可持續(xù)發(fā)展具有重要意義。

    祁連山是青藏高原東北部的高大山系,是我國(guó)西北干旱區(qū)重要的生態(tài)安全屏障,該區(qū)域地形復(fù)雜,氣候變化明顯,屬于典型的高寒區(qū)[17-18]。近年來(lái),隨著祁連山國(guó)家公園的建立以及區(qū)內(nèi)生態(tài)價(jià)值的提高,以祁連山脈為核心的高寒旱區(qū)研究成為熱點(diǎn)之一。因此,本研究以祁連山南坡境內(nèi)的林地、灌叢、草地及耕地為研究對(duì)象,綜合土壤物理、化學(xué)和生物學(xué)特征,運(yùn)用主成分分析方法選擇土壤質(zhì)量指標(biāo)建立MDS,評(píng)估研究區(qū)域的土壤質(zhì)量,以期為青藏高原地區(qū)土壤可持續(xù)利用提供參考。

    1" 材料與方法

    1.1" 研究區(qū)概況

    本文中的祁連山南坡主要指的是祁連山脈的中國(guó)青海片區(qū),該區(qū)域年日照時(shí)數(shù)約為2200~2900 h,年均氣溫僅為-5.9℃,年內(nèi)降水集中分布在6~8月,年均降水量300~400 mm,屬于典型的高寒區(qū),具有冬長(zhǎng)嚴(yán)寒、夏短清涼的特征。研究區(qū)地理位置為 98°08′13″~102°38′16″E,37°03′17″~39°05′56″N,海拔2286~5208 m,平均海拔3800 m,總面積 2.4×104 km2,植物生長(zhǎng)期為100~150 d,與降水季同期。研究區(qū)土壤呈弱堿性,土壤類(lèi)型包括山地森林土、灰褐土、栗鈣土、黑鈣土、高山草原土、草甸土、寒漠土等[19]。

    1.2" 樣品采集與分析

    于2022年8月對(duì)區(qū)內(nèi)進(jìn)行實(shí)地調(diào)查后,選擇無(wú)降雨日期,按照隨機(jī)采樣法在每個(gè)樣地選取3個(gè)1 m×1 m的取樣點(diǎn)用直徑為5 cm的土鉆鉆取采集林地、灌叢、草地及耕地土壤(0~20 cm)樣品?,F(xiàn)場(chǎng)記錄采樣點(diǎn)的地理位置等基本信息,如圖1所示,布設(shè)樣地林地15個(gè)、灌叢15個(gè)、草地9個(gè)、耕地19個(gè),共布設(shè)58個(gè)樣地,共計(jì)樣品174件。采樣后取部分土壤帶回實(shí)驗(yàn)室,經(jīng)自然風(fēng)干,剔除大塊礫石及未分解植物根莖后,部分用于測(cè)定土壤機(jī)械組成,部分過(guò)篩,用于測(cè)定土壤理化性質(zhì)。此外,一份鮮土現(xiàn)場(chǎng)過(guò)2 mm土壤篩后,置于-80℃液氮保存,用于測(cè)定土壤微生物多樣性。

    土壤含水量(Soil Water Content,SWC)采用烘干法測(cè)定;土壤容重(Soil Bulk Density,BD)采用環(huán)刀法測(cè)定;土壤機(jī)械組成采用Mastersizer 2000型激光粒度儀測(cè)定;土壤電導(dǎo)率(Electric Conductivity,EC)使用電導(dǎo)率儀(DDS-307)進(jìn)行測(cè)定。土壤全氮(Total Nitrogen,TN)、全磷(Total Phosphorus,TP)、全鉀(Total Potassium,TK)、有機(jī)碳(Soil Organic Carbon,SOC)含量等理化性質(zhì)的測(cè)定參考《土壤農(nóng)化分析》[20]。利用高通量測(cè)序技術(shù)對(duì)土壤細(xì)菌和真菌群落多樣性進(jìn)行分析。

    1.3" 土壤質(zhì)量評(píng)價(jià)

    本研究采用13個(gè)影響土壤質(zhì)量的指標(biāo),對(duì)全量數(shù)據(jù)集的降維及土壤綜合質(zhì)量指數(shù)進(jìn)行綜合評(píng)價(jià),通過(guò)主成分分析(PCA)將特征值≥1的指標(biāo)進(jìn)行分組,各組內(nèi)選取Norm最高值相差10%以內(nèi)的指標(biāo),計(jì)算綜合荷載建立MDS,所有主成分的綜合荷載越大,其解釋綜合信息的能力就越強(qiáng)[21-22];總數(shù)據(jù)集(Total data set,TDS)是對(duì)所有指標(biāo)進(jìn)行計(jì)算綜合荷載計(jì)算,公式如下:

    Nik=∑k1(U2ikλk)(1)

    式中:Nik是i指標(biāo)在全部k個(gè)主成分的累加因子荷載;Uik是i指標(biāo)的單個(gè)因子荷載;λk為第k個(gè)主成分的特征值。

    通過(guò)建立指標(biāo)與土壤質(zhì)量之間的隸屬函數(shù),對(duì)進(jìn)入MDS的指標(biāo)求隸屬度值[5],計(jì)算公式如下:

    u(x)=0.9×1,x≥x2x-x1x2-x1,x2lt;x<x10,x≤x1(2)

    式中:x為土壤質(zhì)量指標(biāo)的實(shí)際測(cè)量值;x1表示最小值;x2表示最大值;u(x)表示指標(biāo)標(biāo)準(zhǔn)化后的值。

    把各項(xiàng)指標(biāo)線性得分與權(quán)重因子進(jìn)行加權(quán)求和,計(jì)算土壤質(zhì)量指數(shù)(SQI),公式如下[23]:

    SQI=∑ni=1WiNi(3)

    式中:Ni和Wi分別表示第i種評(píng)價(jià)指標(biāo)所對(duì)應(yīng)的線性得分和權(quán)重系數(shù);n為數(shù)據(jù)集中土壤指標(biāo)總數(shù)。

    1.4" 數(shù)據(jù)處理

    使用IBM SPSS 26軟件對(duì)數(shù)據(jù)進(jìn)行正態(tài)性、主成分及相關(guān)性等統(tǒng)計(jì)學(xué)分析,數(shù)據(jù)分布的正態(tài)性采用Shapiro-Wilk檢驗(yàn),用方差分析(One-way ANOVA)和Duncan檢驗(yàn)比較不同數(shù)據(jù)組之間的差異,進(jìn)行主成分分析時(shí)考慮了所有參數(shù)。采用Origin 2024軟件繪制圖表。

    2" 結(jié)果與分析

    2.1" 祁連山南坡不同土地利用方式下土壤性質(zhì)差異

    如表1所示,祁連山南坡4種不同土地利用方式下土壤性質(zhì)呈現(xiàn)差異性。就土壤TN含量而言,表現(xiàn)為耕地顯著大于灌叢及林地(P<0.05),灌叢及草地顯著大于林地(P<0.05)。耕地土壤TP含量最高,顯著大于林地、灌叢和草地土壤TP含量(P<0.05)。草地土壤TK含量最高,顯著大于林地、灌叢及耕地土壤TK含量(P<0.05),灌叢與耕地之間土壤TK含量差異不顯著,但均顯著高于林地(P<0.05)。灌叢與草地之間SOC含量差異不顯著,但均顯著高于林地與耕地(P<0.05),且林地SOC含量也顯著高于耕地(P<0.05)。草地SWC最高,顯著高于其他3種土地利用方式下土壤SWC(P<0.05);灌叢與耕地土壤SWC均顯著低于林地土壤SWC(P<0.05)。祁連山南坡土壤呈弱堿性,耕地土壤pH值最大,灌叢土壤pH值最小,顯著低于耕地和林地(P<0.05)。祁連山南坡4種不同土地利用方式下土壤EC差異均不顯著。林地土壤黏粒(Clay)含量最高,顯著高于其他3種土地利用方式下土壤Clay含量(P<0.05);灌叢和草地土壤Clay含量均顯著大于耕地土壤Clay含量(P<0.05)。耕地土壤粉粒(Silt)含量顯著大于其他3種土地利用方式下土壤Silt含量(P<0.05);林地和灌叢土壤Silt含量均顯著小于草地土壤Silt含量(P<0.05)。林地和灌叢土壤砂粒(Sand)含量均顯著大于草地與耕地土壤Sand含量(P<0.05);且草地土壤Sand含量顯著大于耕地土壤Sand含量(P<0.05)。灌叢與草地之間土壤細(xì)菌多樣性差異不顯著,但均顯著大于林地與耕地(P<0.05)。林地的土壤真菌多樣性最低,顯著低于其他3種土地利用方式下土壤真菌多樣性(P<0.05)。此外,根據(jù)第二次全國(guó)土壤普查土壤養(yǎng)分含量分級(jí)標(biāo)準(zhǔn)[24],境內(nèi)林地、灌叢、草地和耕地TN,TP和SOC含量極豐富,處于一級(jí)水平,而土壤TP含量則相對(duì)較低,處于二級(jí)水平。

    2.2" 基于PAC的MDS建立

    如表2所示,祁連山南坡林地中,特征值>1的5個(gè)主成分(PC1,PC2,PC3,PC4,PC5)占總方差的85.137%。指標(biāo)因子荷載值越大,其所占主成分權(quán)重越大,對(duì)土壤影響越大。第1主成分中土壤指標(biāo)有TN、SOC和pH,第2主成分中土壤指標(biāo)有BD,Clay,Silt和Sand,第3主成分中土壤指標(biāo)有TP、SWC和真菌多樣性,第4主成分中土壤指標(biāo)為細(xì)菌多樣性,第5主成分中土壤指標(biāo)有TK和EC,根據(jù)Norm值和相關(guān)性選取原則,篩選Norm值最高指標(biāo)相差10%以內(nèi)的指標(biāo)[5]。第1主成分中Norm值的排序?yàn)椋簆H>TN>SOC,TN,SOC和pH均在最高指標(biāo)Norm值10%以內(nèi),但由于林地pH與TN和SOC之間均呈極顯著(P<0.01)相關(guān),因此,選擇pH進(jìn)入MDS。同理,第2主成分中,Norm值10%以內(nèi)的土壤指標(biāo)為BD,Silt和Sand,Sand與BD呈顯著相關(guān)(P<0.05),與Silt呈極顯著相關(guān)(P<0.01),因此Sand進(jìn)入MDS;第3主成分中,Norm值10%以內(nèi)的土壤指標(biāo)為SWC和真菌多樣性,二者無(wú)顯著相關(guān)關(guān)系,因此SWC和真菌多樣性進(jìn)入MDS;第4主成分中細(xì)菌多樣性進(jìn)入MDS;第5主成分中Norm值10%以內(nèi)的土壤指標(biāo)為T(mén)K和EC,二者無(wú)顯著相關(guān)關(guān)系,因此TK和EC進(jìn)入MDS。

    如表3所示,祁連山南坡灌叢中,特征值>1的3個(gè)主成分(PC1、PC2、PC3)占總方差的73.780%。通過(guò)主成分分析,第1主成分中土壤指標(biāo)有TK,SOC,SWC,pH,BD,細(xì)菌多樣性和真菌多樣性,第2主成分中土壤指標(biāo)有EC,Clay,Silt和Sand,第3主成分中土壤指標(biāo)有TN和TP。根據(jù)Norm值和相關(guān)性選取原則,最終進(jìn)入MDS的指標(biāo)為SOC,Sand和TN。

    如表4所示,祁連山南坡草地中,特征值>1的4個(gè)主成分(PC1,PC2,PC3,PC4)占總方差的92.154%。通過(guò)主成分分析,第1主成分中土壤指標(biāo)有TP,SOC,pH,EC,Clay,Silt,Sand和真菌多樣性,第2主成分中土壤指標(biāo)有TN和細(xì)菌多樣性,第3主成分中土壤指標(biāo)有TK和BD,第4主成分中土壤指標(biāo)有SWC。根據(jù)Norm值和相關(guān)性選取原則,最終進(jìn)入MDS的指標(biāo)為Sand,TP,pH,EC,TN,BD和SWC。

    如表5所示,祁連山南坡耕地中,特征值>1的4個(gè)主成分(PC1,PC2,PC3,PC4)占總方差的78.790%。通過(guò)主成分分析,第1主成分中土壤指標(biāo)有TN,TP,SWC,pH,BD,Clay,Silt和Sand,第2主成分中土壤指標(biāo)有TK,SOC和細(xì)菌多樣性,第3主成分中土壤指標(biāo)有EC,第4主成分中土壤指標(biāo)有真菌多樣性。根據(jù)Norm值和相關(guān)性選取原則,最終進(jìn)入MDS的指標(biāo)為Sand,SOC,EC和真菌多樣性。

    2.3" 基于MDS的土壤質(zhì)量評(píng)價(jià)

    全部評(píng)價(jià)指標(biāo)標(biāo)準(zhǔn)化處理后進(jìn)行主成分分析,用指標(biāo)公因子方差占所有指標(biāo)公因子方差之和的比例來(lái)獲取各指標(biāo)的權(quán)重(表6-7),并使用隸屬度函數(shù)確定各指標(biāo)的隸屬度。結(jié)合各指標(biāo)的權(quán)重值與隸屬度,基于TDS和MDS計(jì)算得到的林地土壤質(zhì)量指數(shù)值的范圍分別為0.461~0.657和0.448~0.667,平均值分別為0.544和0.535;灌叢土壤質(zhì)量指數(shù)值的范圍分別為0.413~0.652和0.419~0.666,平均值分別為0.536和0.519;草地土壤質(zhì)量指數(shù)值的范圍分別為0.418~0.549和0.359~0.588,平均值分別為0.485和0.466;耕地土壤質(zhì)量指數(shù)值的范圍分別為0.443~0.632和0.391~0.616,平均值分別為0.550和0.544;因此,祁連山南坡農(nóng)田不同土地利用方式下MDS土壤質(zhì)量指數(shù)表現(xiàn)為耕地>林地>灌叢>草地。

    2.4" MDS評(píng)價(jià)的精度驗(yàn)證

    MDS構(gòu)建的合理性,直接關(guān)系到土壤質(zhì)量評(píng)價(jià)的準(zhǔn)確程度。過(guò)度的簡(jiǎn)化指標(biāo),會(huì)導(dǎo)致土壤質(zhì)量評(píng)價(jià)準(zhǔn)確度的下降。因此,需對(duì)比分析基于全部指標(biāo)計(jì)算的土壤質(zhì)量指數(shù)和基于MDS指標(biāo)計(jì)算的土壤質(zhì)量指數(shù)間的差異。將兩類(lèi)數(shù)據(jù)集土壤質(zhì)量指數(shù)值進(jìn)行回歸分析,以TDS計(jì)算得到的土壤質(zhì)量指數(shù)值為基準(zhǔn)對(duì)MDS土壤質(zhì)量指數(shù)值進(jìn)行精度檢驗(yàn),其結(jié)果見(jiàn)圖2。從TDS和MDS的擬合效果來(lái)看,兩類(lèi)數(shù)據(jù)集的土壤質(zhì)量指數(shù)呈顯著正相關(guān)(P<0.01),R2為0.663,證明本研究構(gòu)建的MDS的檢驗(yàn)精度較高,基于MDS構(gòu)建的土壤質(zhì)量指數(shù)可對(duì)研究地塊進(jìn)行不同土地利用類(lèi)型下土壤質(zhì)量評(píng)價(jià)。

    3" 討論

    3.1" 不同土地利用類(lèi)型對(duì)土壤質(zhì)量相關(guān)指標(biāo)的影響

    本文基于PCA法結(jié)合相關(guān)性分析和Norm值,采用MDS計(jì)算SQI,對(duì)祁連山南坡4種不同土地利用方式下土壤質(zhì)量進(jìn)行評(píng)價(jià),并對(duì)MDS的精度進(jìn)行了驗(yàn)證,MDS可以替代TDS對(duì)祁連山南坡不同土地利用類(lèi)型下土壤質(zhì)量進(jìn)行評(píng)價(jià)。結(jié)合各樣地土壤質(zhì)量指數(shù)值來(lái)看(圖2),土壤質(zhì)量指數(shù)值大多分布在中上部,這說(shuō)明土壤質(zhì)量整體較好,通過(guò)保護(hù)林地、灌叢、草地的自然生態(tài)環(huán)境減少人為干擾,對(duì)耕地進(jìn)行合理施用有機(jī)肥等措施,提高土壤保水、保肥能力,促進(jìn)植物根系生長(zhǎng)可有效調(diào)節(jié)該區(qū)域土壤質(zhì)量。

    不同土地利用方式會(huì)導(dǎo)致土壤養(yǎng)分及理化性質(zhì)等存在差異。本文4種不同土地利用方式構(gòu)建的MDS中,土壤pH值、Sand,EC,SWC,SOC,TN含量及真菌多樣性7個(gè)指標(biāo)均入選兩次及以上,是影響研究區(qū)土壤質(zhì)量的關(guān)鍵指標(biāo)。土壤粒度粗細(xì)直接影響著土壤的理化性質(zhì)與肥力狀況,是土壤最基本的物理性質(zhì)之一。土壤顆粒粗會(huì)導(dǎo)致土壤養(yǎng)分淋失,相對(duì)保水性能和保肥能力差[25]。本研究中,耕地土壤Sand含量顯著低于其余3種土地利用方式下土壤Sand含量,因此耕地土壤質(zhì)量較高,這與陳正發(fā)等[22]在云南坡耕地耕層土壤質(zhì)量評(píng)價(jià)結(jié)果一致。土壤EC、pH值及SWC是影響生態(tài)系統(tǒng)健康的重要物理因素[26]。有研究表明,土壤EC、pH值及SWC是影響細(xì)菌類(lèi)群的主要因素,這些理化因子還可通過(guò)改變植物群落結(jié)構(gòu)調(diào)控真菌群落結(jié)構(gòu)[27]。本研究發(fā)現(xiàn)土壤真菌多樣性也是影響土壤質(zhì)量的重要指標(biāo)。張麗星等[28]對(duì)內(nèi)蒙古草地土壤質(zhì)量評(píng)價(jià)研究中表示該研究區(qū)土壤質(zhì)量與土壤EC密切相關(guān)。有研究表示土壤EC、pH值、SWC過(guò)高會(huì)直接影響土壤的生物活性,降低土壤肥力和有機(jī)質(zhì)含量[29-30],進(jìn)而對(duì)土壤質(zhì)量產(chǎn)生消極影響??梢?jiàn),土壤EC、pH值、SWC及真菌多樣性互相影響并制約研究區(qū)土壤質(zhì)量。SOC、TN直接關(guān)系植物生長(zhǎng)及養(yǎng)分供應(yīng)狀況,是植物生長(zhǎng)所需的物質(zhì)保障,對(duì)土壤綜合質(zhì)量的貢獻(xiàn)度較大[31],是評(píng)價(jià)土壤質(zhì)量的敏感指標(biāo)[32]。不同的土地利用方式通過(guò)影響土壤碳和氮的分解與固存,直接影響土壤肥力和植物生長(zhǎng)[33]。

    綜合來(lái)看,影響研究區(qū)域的土壤Sand含量、EC、pH值、SWC,SOC,TN及真菌多樣性幾個(gè)關(guān)鍵指標(biāo),其內(nèi)部聯(lián)系緊密。有研究表明,較大的礫石孔隙結(jié)構(gòu)能夠顯著抑制土壤產(chǎn)流發(fā)生并促進(jìn)水分的垂直入滲,增大土壤表層水分并逐年累積[34],而土壤鹽分運(yùn)移與水分運(yùn)動(dòng)相輔相成,遵循“鹽隨水動(dòng)”的同步性準(zhǔn)則[27]。土壤可溶性鹽分以水分為載體,由蒸發(fā)驅(qū)使聚集至土壤表層,且返鹽效應(yīng)隨覆蓋層結(jié)構(gòu)惡化而愈加強(qiáng)烈[35],這會(huì)導(dǎo)致土壤pH值的改變。此外,土壤微生物作為土地修復(fù)的變局者,直接管控有機(jī)物形成、分解和碳氮循環(huán)過(guò)程,本研究中也證實(shí)了這一點(diǎn)。因此,對(duì)于影響研究區(qū)土壤質(zhì)量的土壤理化性質(zhì)而言,土壤粒度對(duì)于該區(qū)域內(nèi)土壤質(zhì)量的影響較大,且存在一定的連鎖反應(yīng)。在實(shí)施土地管理時(shí),對(duì)耕地要合理施用有機(jī)肥,增加土壤肥力;盡量減少人類(lèi)對(duì)森林、灌叢、草地等自然生態(tài)環(huán)境的干擾,避免水土流失發(fā)生。對(duì)于影響土壤質(zhì)量的微生物,開(kāi)發(fā)應(yīng)用綠色環(huán)保、環(huán)境友好的新型生物技術(shù)用以調(diào)控根際微生物間的互作關(guān)系,激發(fā)有益生態(tài)集群活性[36],如科學(xué)合理施用微生物菌劑,是改善土壤微生態(tài)環(huán)境有效生物途徑。

    3.2" 不同土地利用類(lèi)型的土壤質(zhì)量評(píng)價(jià)

    土壤質(zhì)量評(píng)價(jià)是土地可持續(xù)利用的一個(gè)重要思想。本研究基于MDS計(jì)算的祁連山南坡林地土壤質(zhì)量平均值為0.535,與茍國(guó)花等[5]在青藏高原南部林地的土壤質(zhì)量0.510調(diào)查相類(lèi)似,略高于蔣叢澤等[7]在祁連山北坡調(diào)查的林地土壤質(zhì)量0.481,這可能與祁連山南北坡環(huán)境差異導(dǎo)致林地凋落物豐富度有關(guān)。把黎等[37]及祁棟林等[38]的研究表明,祁連山區(qū)南坡空中云水資源較北坡豐富,本文中林地土壤質(zhì)量的主要制約因子為土壤TK及真菌多樣性,二者的豐富度與降水及氣溫密不可分,這是本研究中祁連山南坡林地土壤質(zhì)量略高于蔣叢澤等[7]在祁連山北坡的研究結(jié)果的重要原因。本研究中灌叢土壤質(zhì)量為0.519,低于大渡河流域河谷區(qū)灌叢的土壤質(zhì)量0.558[39]。就本文研究結(jié)果來(lái)看,制約灌叢土壤質(zhì)量的主要因素是土壤粒度與土壤TN、SOC。馬劍[40]在對(duì)祁連山典型灌叢生態(tài)化學(xué)計(jì)量特征研究結(jié)果表明,祁連山灌叢植物生長(zhǎng)更易受N限制;且眾多研究結(jié)果表明TN、碳含量在土壤中存在著協(xié)同作用[25];就土壤粒度與土壤TN、SOC而言,祁連山南坡位于高寒地區(qū),其土壤機(jī)械組成及土壤養(yǎng)分與大渡河流域中游干暖河谷區(qū)的差異明顯,因此研究區(qū)灌叢土壤質(zhì)量低于王文武等[39]在大渡河流域河谷區(qū)對(duì)灌叢的土壤質(zhì)量研究結(jié)果。祁連山南坡草地土壤質(zhì)量為0.466,與茍國(guó)花等[5]在青藏高原草地的土壤質(zhì)量研究相一致;土壤TN,TP,BD等理化性質(zhì)是影響研究區(qū)草地土壤質(zhì)量的重要因素,這與草地放牧有關(guān);放牧是草地最主要的利用方式之一[41],就其他3種土地利用方式而言,放牧行為會(huì)改變草地土壤物理結(jié)構(gòu),從而影響其土壤質(zhì)量。研究區(qū)耕地土壤質(zhì)量為0.544,高于蔣叢澤等[7]在祁連山北坡的研究結(jié)果耕地土壤質(zhì)量0.462,與劉利昆等[42]在青藏公路沿線的研究結(jié)果0.535類(lèi)似。

    綜合看來(lái),祁連山南坡農(nóng)田不同土地利用方式下MDS土壤質(zhì)量指數(shù)表現(xiàn)為耕地>林地>灌叢>草地。對(duì)土壤質(zhì)量分級(jí)為Ⅰ~Ⅵ級(jí),對(duì)應(yīng)指數(shù)分別為≤0.3,(0.3~0.4],(0.4~0.5],(0.5~0.6],(0.6,0.7]和gt;0.7[5],草地等級(jí)為Ⅲ級(jí),處于“中等”水平;耕地、林地和灌叢土壤質(zhì)量等級(jí)為Ⅳ級(jí),處于“中等偏上”水平。就研究區(qū)實(shí)地調(diào)查來(lái)看,4種土地利用類(lèi)型土壤質(zhì)量等級(jí)差異的主要原因如下:林地較高的凋落物豐富度有利于土壤質(zhì)量的改善,且耕地的人為施肥,可促進(jìn)有機(jī)質(zhì)積累與養(yǎng)分循環(huán)[43],而草地因?yàn)榉拍?、牲畜踐踏等影響土壤肥力,導(dǎo)致土壤質(zhì)量略低于耕地、林地和灌叢。有研究表明,祁連山地區(qū)的喬木建群種與優(yōu)勢(shì)種青海云杉,其葉片凋落物分解是祁連山中段森林生態(tài)系統(tǒng)物質(zhì)循環(huán)和能量流動(dòng)的重要環(huán)節(jié)之一[44],農(nóng)田肥料施用改變了土壤養(yǎng)分[25],較草地而言提高了土壤質(zhì)量。

    4" 結(jié)論

    本研究基于最小數(shù)據(jù)集對(duì)祁連山南坡林地、灌叢、草地和耕地土壤質(zhì)量進(jìn)行了綜合評(píng)價(jià)。結(jié)果表明,林地、灌叢、草地和耕地土壤質(zhì)量指數(shù)值分別為0.535,0.519,0.466和0.544,表現(xiàn)為耕地>林地>灌叢>草地。土壤質(zhì)量關(guān)鍵指標(biāo)間存在互相影響,對(duì)于研究區(qū)域土地要實(shí)施分類(lèi)科學(xué)管理,對(duì)耕地合理施用有機(jī)肥,增加土壤肥力;對(duì)于林地、灌叢、草地等自然生態(tài)環(huán)境盡量減少人為干擾,避免水土流失發(fā)生。此外,對(duì)于影響土壤質(zhì)量的微生物,開(kāi)發(fā)應(yīng)用綠色環(huán)保、環(huán)境友好的新型生物技術(shù)是應(yīng)對(duì)影響研究區(qū)土壤質(zhì)量的微生物指標(biāo)的有效措施。研究結(jié)果為祁連山地區(qū)土地可持續(xù)利用合理開(kāi)發(fā)利用與管理提供了重要參考。

    參考文獻(xiàn)

    [1]" KARLEN D L,MAUSBACH M J. et al. Soil Quality:A Concept,Definition,and framework for evaluation (A Guest Editorial)[J]. Soil Science Society of America Journal,1997,61(1):4-10

    [2]" POWER J F,MYERS R J K. The maintenance or improvement of farming systems in North America and Australia[C]. Soil quality in semi-arid agriculture. Saskatoon,Saskatchewan,Canada,1989:273-292

    [3]" 劉慶宇,馬瑛,程莉. 青海門(mén)源縣土壤質(zhì)量地球化學(xué)評(píng)價(jià)[J]. 地質(zhì)與勘探,2022,58(3):609-618

    [4]" 趙其國(guó),吳志東. 深入開(kāi)展“土壤與環(huán)境”問(wèn)題的研究[J]. 土壤與環(huán)境,1999,8(1):1-4

    [5]" 茍國(guó)花,樊軍,王茜,等. 基于最小數(shù)據(jù)集的青藏高原南北部不同土地利用方式土壤質(zhì)量評(píng)價(jià)[J]. 應(yīng)用生態(tài)學(xué)報(bào),2023,34(5):1360-1366

    [6]" AULAKH C S,SHARMA S,THAKUR M,et al. A review of the influences of organic farming on soil quality,crop productivity and produce quality[J]. Journal of Plant Nutrition,2022,45(12),1884-1905

    [7]" 蔣叢澤,受娜,高瑋,等. 青藏高原東北緣不同土地利用類(lèi)型土壤質(zhì)量綜合評(píng)價(jià)[J]. 應(yīng)用生態(tài)學(xué)報(bào),2022,33(12):3279-3286

    [8]" JAT R A,REDDY K K,CHOUDHARY R R,et al. Effect of conservation agriculture practices on soil quality,productivity,and profitability of peanut-based system of Saurashtra,India[J]. Agronomy Journal,2021,113(2),2102-2117

    [9]" TAKOUTSING B,WEBER J,AYNEKULU E,et al. Assessment of soil health indicators for sustainable production of maize in smallholder farming systems in the highlands of Cameroon[J]. Geoderma,2016,276,64-73

    [10]李霞,朱萬(wàn)澤,舒樹(shù)淼,等. 基于主成分分析的大渡河中游干暖河谷草地土壤質(zhì)量評(píng)價(jià)[J]. 生態(tài)學(xué)報(bào),2021,41(10):3891-3900

    [11]LASON W E,PIERCE F J. The dynamics of soil quality as ameasure of sustainable management[C]//DORAN J W,COLEMAN D C,BEZEDICK D F,et al. Defining soil quality for sustainable environment. America:Soil Science Society,1994:37-51

    [12]VELASQUEZ E,LAVELLE P,BARRIOS E,et al. Evaluating soil quality in tropical agroecosystems of Colombia using NIRS[J]. Soil Biology and Biochemistry,2005,37(5),889-898

    [13]趙西寧,劉帥,高曉東,等. 不同改良劑對(duì)黃土高原丘陵區(qū)山地果園土壤質(zhì)量的影響[J]. 生態(tài)學(xué)報(bào),2022,42(17):7080-7091

    [14]桑亞轉(zhuǎn),尤楊,李多才,等. 放牧對(duì)祁連山高寒典型草原土壤質(zhì)量的影響[J]. 生態(tài)學(xué)報(bào),2023,43(15):6364-6377

    [15]陳璟,楊寧. 衡陽(yáng)紫色土丘陵坡地植被恢復(fù)過(guò)程中土壤質(zhì)量評(píng)價(jià)[J]. 草地學(xué)報(bào),2018,26(1):160-167

    [16]田英,許喆,王婭麗,等. 寧夏銀川平原沙化土地不同人工林土壤質(zhì)量評(píng)價(jià)[J]. 生態(tài)學(xué)報(bào),2023,43(4):1515-1525

    [17]張婷,朱曉鵬,徐海燕,等. 祁連山草地生態(tài)系統(tǒng)優(yōu)勢(shì)物種生長(zhǎng)季的生態(tài)化學(xué)計(jì)量特征[J/OL]. http://kns.cnki.net/kcms/detail/11.3362.S.20240407.1956.008.html,2024-04-09/2024-06-11

    [18]JUICE S M,SCHABERG P G,KOSIBA A M,et al. Soil type modifies the impacts of warming and snow exclusion on leachate carbon and nutrient losses[J]. Biogeochemistry,2022,160(2),199-217

    [19]刁二龍,曹廣超,曹生奎,等. 祁連山南坡不同土地利用方式下土壤碳氮含量及通徑分析[J]. 干旱區(qū)研究,2021,38(5):1346-1354

    [20]鮑士旦. 土壤農(nóng)化分析[M]. 北京:中國(guó)農(nóng)業(yè)出版社,2000:25-114

    [21]ABDEL-FATTAH M K,MOHAMED E S,WAGDI E M,et al. Quantitative evaluation of soil quality using principal component analysis:The case study of El-Fayoum depression Egypt[J]. Sustainability,2021,13(4),1824

    [22]陳正發(fā),史東梅,金慧芳,等. 基于土壤管理評(píng)估框架的云南坡耕地耕層土壤質(zhì)量評(píng)價(jià)[J]. 農(nóng)業(yè)工程學(xué)報(bào),2019,35(3):256-267

    [23]LI P,SHI K,WANG Y,et al. Soil quality assessment of wheat-maize cropping system with different productivities in China:Establishing a minimum data set[J]. Soil and Tillage Research,2019,190,31-40

    [24]全國(guó)土壤普查辦公室. 中國(guó)土壤[M]. 北京:中國(guó)農(nóng)業(yè)出版社,1998:45-72

    [25]邱巡巡,曹廣超,張卓,等. 高寒農(nóng)田土壤有機(jī)碳和全氮密度垂直分布特征及其與海拔的關(guān)系[J]. 土壤通報(bào),2022,53 (3):623-630

    [26]周雨靜,賀晶,王雪松,等. 呼倫貝爾草原不同利用方式下土壤鹽分分布特征[J]. 草地學(xué)報(bào),2024,32(2):470-479

    [27]賈振江,劉學(xué)智,李王成,等. 旱區(qū)連作砂田土壤質(zhì)量和土地生產(chǎn)力演變與調(diào)控研究進(jìn)展[J]. 生態(tài)學(xué)報(bào),2024,44(5):2136-2148

    [28]張麗星,海春興,常耀文,等. 羊草及芨芨草草原和西北針茅草原土壤質(zhì)量評(píng)價(jià)[J]. 草業(yè)學(xué)報(bào),2021,30(4):68-79

    [29]王飛. 鹽分對(duì)土壤微生物多樣性及土壤有機(jī)物轉(zhuǎn)化的影響[D]. 新疆:石河子大學(xué),2011:54-55

    [30]金慧芳,史東梅,陳正發(fā),等. 基于聚類(lèi)及PCA分析的紅壤坡耕地耕層土壤質(zhì)量評(píng)價(jià)指標(biāo)[J]. 農(nóng)業(yè)工程學(xué)報(bào),2018,34(7):155-164

    [31]翟輝,張海,張超,等. 黃土峁?fàn)钋鹆陞^(qū)不同類(lèi)型林分土壤微生物功能多樣性[J]. 林業(yè)科學(xué),2016,52(12) 84-91

    [32]任書(shū)杰,曹明奎,陶波,等. 陸地生態(tài)系統(tǒng)氮狀態(tài)對(duì)碳循環(huán)的限制作用研究進(jìn)展[J]. 地理科學(xué)進(jìn)展,2006,(4):58-67

    [33]DALAL R C,CHAN K Y. Soil organic matter in rain fed cropping systems of the Australian cereal belt[J].Australian Journal of Soil Research,2001,39:435-464

    [34]鄭飛龍,王西娜,葛敏,等. 壓砂條件下灰鈣土地區(qū)農(nóng)田土壤水分垂直入滲特征研究[J]. 土壤通報(bào),2021,52(2):314-321

    [35]徐天淵,賈振江,李王成,等. 寧夏中部干旱帶微咸水灌溉對(duì)砂土混合覆蓋下土壤水鹽運(yùn)移的影響[J]. 干旱地區(qū)農(nóng)業(yè)研究,2021,39(5):138-144

    [36]耿文叢,馬悅,張玉雪,等. 設(shè)施農(nóng)業(yè)的土壤健康調(diào)控技術(shù)研究進(jìn)展[J]. 中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào)(中英文),2022,30(12):1973-1984

    [37]把黎,尹憲志,龐朝云,等. 祁連山地區(qū)夏季南坡與北坡空中云水資源差異性分析[J]. 干旱區(qū)研究,2022,39(5):1345-1359

    [38]祁棟林,李甫,肖建設(shè),等. 近53a來(lái)祁連山南北坡潛在蒸散量及地表濕潤(rùn)度變化趨勢(shì)分析[J]. 干旱氣象,2016,34(1):26-33

    [39]王文武,朱萬(wàn)澤,李霞,等. 基于最小數(shù)據(jù)集的大渡河干暖河谷典型植被土壤質(zhì)量評(píng)價(jià)[J]. 中國(guó)水土保持科學(xué)(中英文),2021,19(6):54-59

    [40]馬劍. 祁連山北麓中段典型灌叢群落生態(tài)化學(xué)計(jì)量特征研究[D]. 蘭州:甘肅農(nóng)業(yè)大學(xué),2021:68-69

    [41]MACCI C,DONI S,BONDI G,et al. Effects of wild boar (Sus scrofa) grazing on soil properties in Mediterranean environment[J]. Catena,2012,98:79-86

    [42]劉利昆,趙廣舉,穆興民,等. 基于最小數(shù)據(jù)集的青藏公路沿線土壤質(zhì)量評(píng)價(jià)[J]. 水土保持研究,2022,29(2):125-131

    [43]朱鳴鳴,徐鍍涵,陳光燕,等. 基于最小數(shù)據(jù)集的喀斯特不同利用方式下土壤質(zhì)量評(píng)價(jià)[J]. 草地學(xué)報(bào),2021,29(10):2323-2331

    [44]顧繼雄,周碧蓮,韓昊,等. 青海云杉林凋落物分解早期微生物生物量時(shí)空動(dòng)態(tài)變化[J]. 生態(tài)學(xué)報(bào),2023,43(1):48-59

    (責(zé)任編輯" 付" 宸)

    日本 av在线| 国产精品 国内视频| www.www免费av| 欧美中文日本在线观看视频| 日韩国内少妇激情av| 黄色 视频免费看| 久久天堂一区二区三区四区| 日本 欧美在线| 欧美成人免费av一区二区三区| 国产亚洲精品av在线| 亚洲电影在线观看av| 中文字幕人成人乱码亚洲影| av在线天堂中文字幕| 午夜精品在线福利| 高潮久久久久久久久久久不卡| av福利片在线观看| 在线永久观看黄色视频| 高清毛片免费观看视频网站| 亚洲avbb在线观看| av福利片在线观看| 草草在线视频免费看| 日本a在线网址| aaaaa片日本免费| 成人欧美大片| 亚洲av成人不卡在线观看播放网| 久久久久久久午夜电影| 亚洲成人久久爱视频| 久久久精品大字幕| 国产99白浆流出| 精品国产美女av久久久久小说| 99国产精品一区二区蜜桃av| 国产精品永久免费网站| 一区二区三区激情视频| 亚洲av美国av| 禁无遮挡网站| 精品国内亚洲2022精品成人| 亚洲成人久久性| 99精品欧美一区二区三区四区| 中国美女看黄片| 蜜桃久久精品国产亚洲av| 国产伦一二天堂av在线观看| 国产免费男女视频| 亚洲av中文字字幕乱码综合| 两人在一起打扑克的视频| 久久久久国内视频| 亚洲av第一区精品v没综合| 欧美黑人精品巨大| 麻豆久久精品国产亚洲av| 精品国产乱子伦一区二区三区| 亚洲国产欧美一区二区综合| av在线播放免费不卡| 国产三级中文精品| 美女免费视频网站| 久久精品综合一区二区三区| 毛片女人毛片| 亚洲免费av在线视频| 99热这里只有精品一区 | 九色成人免费人妻av| 免费观看人在逋| 精品久久久久久久毛片微露脸| а√天堂www在线а√下载| 男女之事视频高清在线观看| 亚洲av熟女| 日本 欧美在线| 欧美成狂野欧美在线观看| 亚洲五月天丁香| 变态另类成人亚洲欧美熟女| 一本久久中文字幕| 美女高潮喷水抽搐中文字幕| 色噜噜av男人的天堂激情| 婷婷精品国产亚洲av| 人人妻人人看人人澡| 国产精品亚洲av一区麻豆| 国产日本99.免费观看| 亚洲人成77777在线视频| 久久久国产欧美日韩av| 人人妻人人澡欧美一区二区| 在线观看一区二区三区| 狠狠狠狠99中文字幕| 欧美乱妇无乱码| 欧美人与性动交α欧美精品济南到| 18禁黄网站禁片免费观看直播| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲中文日韩欧美视频| 亚洲欧洲精品一区二区精品久久久| 亚洲欧洲精品一区二区精品久久久| 国产成+人综合+亚洲专区| 久久人人精品亚洲av| 男女之事视频高清在线观看| 精品第一国产精品| 亚洲专区中文字幕在线| 午夜福利免费观看在线| 成熟少妇高潮喷水视频| 中文亚洲av片在线观看爽| 亚洲国产欧洲综合997久久,| 老司机午夜十八禁免费视频| 欧美日韩精品网址| 国产精品香港三级国产av潘金莲| 欧美一区二区国产精品久久精品 | 亚洲av熟女| 日韩欧美 国产精品| 亚洲国产中文字幕在线视频| 非洲黑人性xxxx精品又粗又长| 亚洲精品中文字幕一二三四区| 免费在线观看完整版高清| 午夜a级毛片| 国产精品九九99| 中文字幕精品亚洲无线码一区| 91在线观看av| 村上凉子中文字幕在线| 日韩三级视频一区二区三区| 高清毛片免费观看视频网站| 99在线视频只有这里精品首页| 黄色视频,在线免费观看| 亚洲片人在线观看| 午夜福利高清视频| 亚洲人成网站在线播放欧美日韩| 欧美一区二区精品小视频在线| 午夜免费观看网址| 99久久国产精品久久久| 欧美精品亚洲一区二区| 国产三级中文精品| 亚洲av成人av| 国产高清有码在线观看视频 | 一二三四社区在线视频社区8| 国产又黄又爽又无遮挡在线| 国产高清激情床上av| 久久久水蜜桃国产精品网| 午夜激情福利司机影院| 亚洲精品一卡2卡三卡4卡5卡| 久久精品夜夜夜夜夜久久蜜豆 | 琪琪午夜伦伦电影理论片6080| 色噜噜av男人的天堂激情| 欧美又色又爽又黄视频| 亚洲色图av天堂| 色老头精品视频在线观看| 欧美性长视频在线观看| 757午夜福利合集在线观看| 好男人电影高清在线观看| 女生性感内裤真人,穿戴方法视频| 欧美色欧美亚洲另类二区| 不卡av一区二区三区| 婷婷六月久久综合丁香| 国内精品久久久久久久电影| 欧美精品亚洲一区二区| 黄色视频不卡| 国产单亲对白刺激| 亚洲精品在线美女| 国产精品av久久久久免费| 校园春色视频在线观看| 亚洲成人免费电影在线观看| 国产私拍福利视频在线观看| 18禁国产床啪视频网站| 男女午夜视频在线观看| 黑人巨大精品欧美一区二区mp4| 久久伊人香网站| 午夜成年电影在线免费观看| 高清在线国产一区| 91麻豆av在线| 中文在线观看免费www的网站 | 天堂av国产一区二区熟女人妻 | 久久久久免费精品人妻一区二区| 老司机在亚洲福利影院| 亚洲av日韩精品久久久久久密| 99久久综合精品五月天人人| 三级国产精品欧美在线观看 | 国内精品久久久久久久电影| 日韩中文字幕欧美一区二区| 久久人妻福利社区极品人妻图片| 国产在线精品亚洲第一网站| www日本黄色视频网| 日韩高清综合在线| 亚洲专区字幕在线| 久久久水蜜桃国产精品网| 成年免费大片在线观看| 三级国产精品欧美在线观看 | 精品欧美一区二区三区在线| 12—13女人毛片做爰片一| 免费观看人在逋| 成年免费大片在线观看| 亚洲精品中文字幕在线视频| 岛国在线免费视频观看| 成人一区二区视频在线观看| 每晚都被弄得嗷嗷叫到高潮| 国产不卡一卡二| 欧美黑人精品巨大| 一进一出抽搐gif免费好疼| 在线观看66精品国产| 看免费av毛片| 熟妇人妻久久中文字幕3abv| 麻豆成人av在线观看| 又紧又爽又黄一区二区| 又粗又爽又猛毛片免费看| 久久久久久久久久黄片| 欧美一级毛片孕妇| 又粗又爽又猛毛片免费看| 亚洲国产精品成人综合色| 亚洲欧美日韩高清在线视频| 搡老岳熟女国产| 又大又爽又粗| 法律面前人人平等表现在哪些方面| 中文在线观看免费www的网站 | 久久香蕉激情| 精品熟女少妇八av免费久了| 婷婷丁香在线五月| 国产精品久久电影中文字幕| 久久久久久久久免费视频了| 国产高清videossex| 亚洲人成伊人成综合网2020| svipshipincom国产片| 久久人妻av系列| 国产高清视频在线播放一区| 长腿黑丝高跟| 国产97色在线日韩免费| 欧美不卡视频在线免费观看 | 一级片免费观看大全| 老司机福利观看| 国产精品亚洲一级av第二区| 国产精品一区二区三区四区久久| 好男人在线观看高清免费视频| 精品高清国产在线一区| 国产野战对白在线观看| 久久中文字幕一级| 亚洲一区高清亚洲精品| 女人爽到高潮嗷嗷叫在线视频| 最好的美女福利视频网| 亚洲欧美一区二区三区黑人| 欧美色欧美亚洲另类二区| 亚洲av中文字字幕乱码综合| 欧美日韩黄片免| 变态另类丝袜制服| 最近最新免费中文字幕在线| 亚洲第一欧美日韩一区二区三区| 99热这里只有是精品50| 精品不卡国产一区二区三区| 97人妻精品一区二区三区麻豆| 1024香蕉在线观看| 99国产精品一区二区蜜桃av| 国产精华一区二区三区| 美女大奶头视频| 国产真人三级小视频在线观看| 日韩欧美三级三区| 午夜亚洲福利在线播放| 亚洲欧美日韩高清专用| 少妇的丰满在线观看| 中文字幕最新亚洲高清| 欧美绝顶高潮抽搐喷水| 一二三四在线观看免费中文在| 美女大奶头视频| 99久久精品热视频| 成年版毛片免费区| 波多野结衣高清无吗| 宅男免费午夜| 日本 欧美在线| 久久久久久亚洲精品国产蜜桃av| 色综合婷婷激情| 午夜免费成人在线视频| 少妇被粗大的猛进出69影院| 19禁男女啪啪无遮挡网站| 99精品久久久久人妻精品| 欧美一级a爱片免费观看看 | 国产精品永久免费网站| 亚洲第一欧美日韩一区二区三区| 男男h啪啪无遮挡| 精品不卡国产一区二区三区| 国产一区二区在线观看日韩 | 国产亚洲精品久久久久久毛片| 欧美av亚洲av综合av国产av| 在线观看美女被高潮喷水网站 | 91大片在线观看| 国产99白浆流出| www.999成人在线观看| 后天国语完整版免费观看| 日本三级黄在线观看| 99国产精品99久久久久| 女同久久另类99精品国产91| 欧美在线黄色| 免费看十八禁软件| 亚洲一区中文字幕在线| 亚洲欧美一区二区三区黑人| 国产伦一二天堂av在线观看| 久久精品综合一区二区三区| 脱女人内裤的视频| 亚洲国产欧洲综合997久久,| 日本一本二区三区精品| 黄片小视频在线播放| 国产激情久久老熟女| 色在线成人网| 亚洲一区高清亚洲精品| 国产亚洲精品第一综合不卡| 亚洲成人久久爱视频| 精品国产超薄肉色丝袜足j| 欧美在线一区亚洲| 国产精品乱码一区二三区的特点| 亚洲九九香蕉| 在线观看舔阴道视频| 老汉色av国产亚洲站长工具| 欧美大码av| 久久精品夜夜夜夜夜久久蜜豆 | ponron亚洲| 又大又爽又粗| 国产精品乱码一区二三区的特点| 亚洲,欧美精品.| 最近在线观看免费完整版| 中文字幕久久专区| 在线观看免费视频日本深夜| 国产精品久久久久久久电影 | 黄片大片在线免费观看| 亚洲免费av在线视频| 人人妻人人看人人澡| 日本黄大片高清| 国产精品香港三级国产av潘金莲| 这个男人来自地球电影免费观看| 超碰成人久久| 国内毛片毛片毛片毛片毛片| 中文字幕高清在线视频| 色综合站精品国产| 成人精品一区二区免费| 国产三级黄色录像| svipshipincom国产片| 国产三级中文精品| 在线观看美女被高潮喷水网站 | 国产免费av片在线观看野外av| 男女床上黄色一级片免费看| 国产成人精品久久二区二区91| 亚洲av电影不卡..在线观看| 亚洲av成人一区二区三| 三级国产精品欧美在线观看 | 欧美人与性动交α欧美精品济南到| 又紧又爽又黄一区二区| 这个男人来自地球电影免费观看| 一边摸一边做爽爽视频免费| 黄片大片在线免费观看| 国产午夜福利久久久久久| 在线观看午夜福利视频| 亚洲av片天天在线观看| 99re在线观看精品视频| 亚洲欧美一区二区三区黑人| 首页视频小说图片口味搜索| 很黄的视频免费| 国产激情久久老熟女| 真人做人爱边吃奶动态| 2021天堂中文幕一二区在线观| 1024香蕉在线观看| 俺也久久电影网| 免费一级毛片在线播放高清视频| 日韩精品免费视频一区二区三区| 又爽又黄无遮挡网站| 最好的美女福利视频网| 久久精品国产99精品国产亚洲性色| 日本黄大片高清| 亚洲九九香蕉| 国内久久婷婷六月综合欲色啪| 激情在线观看视频在线高清| 亚洲av中文字字幕乱码综合| 国产精品99久久99久久久不卡| 女人高潮潮喷娇喘18禁视频| 九九热线精品视视频播放| 久久这里只有精品19| 好看av亚洲va欧美ⅴa在| 18禁国产床啪视频网站| 最新美女视频免费是黄的| 1024香蕉在线观看| 老司机午夜福利在线观看视频| 久久精品91蜜桃| 欧美日韩亚洲综合一区二区三区_| 在线观看免费日韩欧美大片| 午夜精品在线福利| 欧美一级毛片孕妇| 国产aⅴ精品一区二区三区波| 三级男女做爰猛烈吃奶摸视频| 最好的美女福利视频网| 亚洲av熟女| 亚洲中文日韩欧美视频| 亚洲18禁久久av| 久久热在线av| 久久中文看片网| 女人被狂操c到高潮| 悠悠久久av| 俄罗斯特黄特色一大片| 国产在线精品亚洲第一网站| 亚洲中文字幕一区二区三区有码在线看 | 亚洲五月天丁香| 一卡2卡三卡四卡精品乱码亚洲| 精品福利观看| 国内毛片毛片毛片毛片毛片| 天堂av国产一区二区熟女人妻 | 免费一级毛片在线播放高清视频| 少妇被粗大的猛进出69影院| 老司机福利观看| 村上凉子中文字幕在线| 久久久精品国产亚洲av高清涩受| 亚洲精品中文字幕在线视频| cao死你这个sao货| 熟妇人妻久久中文字幕3abv| 亚洲精品美女久久久久99蜜臀| 可以免费在线观看a视频的电影网站| 国产精品久久久久久久电影 | av免费在线观看网站| 国产精品,欧美在线| 久久这里只有精品19| 亚洲一区二区三区色噜噜| 色播亚洲综合网| 国产熟女xx| 岛国视频午夜一区免费看| 亚洲第一电影网av| 日韩成人在线观看一区二区三区| 99riav亚洲国产免费| 国产伦一二天堂av在线观看| 99精品欧美一区二区三区四区| 欧美日韩亚洲国产一区二区在线观看| 男人舔女人的私密视频| 久久精品成人免费网站| 不卡av一区二区三区| 长腿黑丝高跟| 色尼玛亚洲综合影院| 亚洲电影在线观看av| 在线观看www视频免费| 国产亚洲欧美在线一区二区| 欧美色欧美亚洲另类二区| 欧美三级亚洲精品| 在线观看美女被高潮喷水网站 | 久久久久久免费高清国产稀缺| 亚洲性夜色夜夜综合| 日本在线视频免费播放| 国产精品99久久99久久久不卡| 一个人免费在线观看的高清视频| 一进一出抽搐动态| 国产亚洲精品一区二区www| 97人妻精品一区二区三区麻豆| 母亲3免费完整高清在线观看| 精品久久久久久久毛片微露脸| av欧美777| 欧美日韩亚洲综合一区二区三区_| 在线观看日韩欧美| 国产午夜精品论理片| 99精品欧美一区二区三区四区| 欧美丝袜亚洲另类 | 极品教师在线免费播放| 亚洲天堂国产精品一区在线| 亚洲国产精品sss在线观看| 三级毛片av免费| 亚洲精品在线美女| 91av网站免费观看| 国产av一区在线观看免费| av天堂在线播放| 午夜精品在线福利| 亚洲专区国产一区二区| 黄色视频,在线免费观看| 男女那种视频在线观看| 国产主播在线观看一区二区| 欧美日韩福利视频一区二区| 午夜激情av网站| 国产伦人伦偷精品视频| 亚洲av片天天在线观看| 久久久国产精品麻豆| 精品一区二区三区av网在线观看| 日韩欧美一区二区三区在线观看| 成年免费大片在线观看| 人妻夜夜爽99麻豆av| 最近最新免费中文字幕在线| 天堂av国产一区二区熟女人妻 | 亚洲,欧美精品.| 美女扒开内裤让男人捅视频| 午夜视频精品福利| 桃红色精品国产亚洲av| 久久久久免费精品人妻一区二区| 99久久精品热视频| 精品久久蜜臀av无| 午夜福利高清视频| 亚洲成av人片免费观看| av免费在线观看网站| 最近最新中文字幕大全免费视频| 99久久99久久久精品蜜桃| 中文字幕人妻丝袜一区二区| 久久草成人影院| 日韩成人在线观看一区二区三区| 在线看三级毛片| 久久欧美精品欧美久久欧美| 一级作爱视频免费观看| 亚洲精品中文字幕在线视频| 国产午夜精品久久久久久| 麻豆成人午夜福利视频| 亚洲熟女毛片儿| 一本久久中文字幕| 亚洲精品在线美女| 欧美人与性动交α欧美精品济南到| 少妇粗大呻吟视频| 一进一出抽搐动态| 午夜福利在线在线| 欧美3d第一页| 国产一区在线观看成人免费| 亚洲国产中文字幕在线视频| 日本三级黄在线观看| 亚洲精品在线观看二区| 成年免费大片在线观看| 国产一区在线观看成人免费| www.熟女人妻精品国产| 亚洲成人国产一区在线观看| 国产99久久九九免费精品| 国产亚洲欧美在线一区二区| 国产亚洲精品av在线| a级毛片a级免费在线| 国产成+人综合+亚洲专区| 国产亚洲精品久久久久5区| 欧美精品啪啪一区二区三区| 亚洲aⅴ乱码一区二区在线播放 | 国产伦一二天堂av在线观看| 久久久精品欧美日韩精品| 看黄色毛片网站| 91成年电影在线观看| 久久久久久久久中文| 久久人妻福利社区极品人妻图片| 欧美日韩瑟瑟在线播放| 亚洲av第一区精品v没综合| 国语自产精品视频在线第100页| 母亲3免费完整高清在线观看| 国产精品免费视频内射| 一级片免费观看大全| 久久香蕉精品热| 毛片女人毛片| 午夜视频精品福利| 国产成人av教育| 91字幕亚洲| x7x7x7水蜜桃| 午夜福利18| 国产精品 欧美亚洲| 免费在线观看完整版高清| 亚洲精品av麻豆狂野| 亚洲av电影在线进入| 黑人操中国人逼视频| 欧美日韩亚洲综合一区二区三区_| 欧美激情久久久久久爽电影| 免费看美女性在线毛片视频| 欧美日本视频| 久久久久亚洲av毛片大全| 亚洲人成77777在线视频| 人妻丰满熟妇av一区二区三区| 国产成年人精品一区二区| 成人18禁高潮啪啪吃奶动态图| 日韩欧美一区二区三区在线观看| 久久午夜综合久久蜜桃| 美女免费视频网站| 久久久久国产一级毛片高清牌| 波多野结衣高清无吗| 99在线视频只有这里精品首页| 国产精品亚洲一级av第二区| 色精品久久人妻99蜜桃| 特级一级黄色大片| 最近最新中文字幕大全免费视频| 久久精品亚洲精品国产色婷小说| 久久久国产精品麻豆| 亚洲av中文字字幕乱码综合| 国产欧美日韩精品亚洲av| 国产精品野战在线观看| 欧美性猛交黑人性爽| 日韩大尺度精品在线看网址| 床上黄色一级片| 久久久久久久久中文| 男女之事视频高清在线观看| 777久久人妻少妇嫩草av网站| 99riav亚洲国产免费| 欧美日韩黄片免| 亚洲国产精品999在线| 狂野欧美激情性xxxx| 美女免费视频网站| 黑人欧美特级aaaaaa片| 中文字幕熟女人妻在线| 亚洲午夜精品一区,二区,三区| 国模一区二区三区四区视频 | 色老头精品视频在线观看| www.www免费av| 又粗又爽又猛毛片免费看| 亚洲欧美精品综合久久99| 久久天堂一区二区三区四区| 国产精品野战在线观看| 听说在线观看完整版免费高清| 国产精品国产高清国产av| 桃红色精品国产亚洲av| 久久精品人妻少妇| 国产视频内射| 亚洲国产看品久久| 久久热在线av| 久久久久久大精品| 成人18禁在线播放| 女同久久另类99精品国产91| 男插女下体视频免费在线播放| 国产一区二区三区视频了| bbb黄色大片| 亚洲中文字幕一区二区三区有码在线看 | 熟女少妇亚洲综合色aaa.| 美女午夜性视频免费| 成人亚洲精品av一区二区| 在线永久观看黄色视频| 亚洲天堂国产精品一区在线| 深夜精品福利| 成人一区二区视频在线观看| 国产一区二区三区在线臀色熟女| 日韩国内少妇激情av| 99国产精品一区二区蜜桃av| 国产精品一区二区精品视频观看| 久久天躁狠狠躁夜夜2o2o| 国产一区二区三区视频了| 美女 人体艺术 gogo| 人人妻人人澡欧美一区二区| 少妇的丰满在线观看| 男女下面进入的视频免费午夜| 一本精品99久久精品77| 亚洲精品中文字幕在线视频| 国产亚洲精品久久久久5区| 欧美不卡视频在线免费观看 | 亚洲va日本ⅴa欧美va伊人久久| 色精品久久人妻99蜜桃| 999精品在线视频| 桃色一区二区三区在线观看|