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

    青藏高原降雨侵蝕力R因子分布特征及驅(qū)動因素

    2024-07-04 09:39:30蘇立彬王茂林袁春慧柯浩成
    人民長江 2024年6期
    關(guān)鍵詞:年際土壤侵蝕青藏高原

    蘇立彬 王茂林 袁春慧 柯浩成

    收稿日期:2023-11-24;接受日期:2024-01-12

    基金項目:西藏土木水利電力工程技術(shù)研究中心項目(XZA202305CHP2002B);西藏自治區(qū)科技廳重點研發(fā)計劃項目(XZ202101ZY0016G)

    作者簡介:蘇立彬,男,講師,博士研究生,研究方向為高原生態(tài)修復與可持續(xù)發(fā)展。E-mail:3249928177@qq.com

    通信作者:柯浩成,男,副教授,博士,主要從事生態(tài)水文與水資源方面的研究。E-mail:kehc@lut.edu.cn

    Editorial Office of Yangtze River. This is an open access article under the CC BY-NC-ND 4.0 license.

    文章編號:1001-4179(2024) 06-0085-07

    引用本文:蘇立彬,王茂林,袁春慧,等.

    青藏高原降雨侵蝕力R因子分布特征及驅(qū)動因素

    [J].人民長江,2024,55(6):85-91.

    摘要:降雨侵蝕力R因子是測算土壤侵蝕的重要因子,該因子的準確計算對水土流失動態(tài)監(jiān)測具有重要意義。以青藏高原為研究對象,利用遙感數(shù)據(jù)和地理信息系統(tǒng)技術(shù),并采用地理探測器模型分析了1991~2020年青藏高原地區(qū)降雨侵蝕力R因子的分布特征以及影響因素。研究結(jié)果表明:① 青藏高原地區(qū)R值的變化在分布上大致呈現(xiàn)西北高東南低的特點,Cv變差系數(shù)在0~0.78之間,年際變化較大的區(qū)域分布在青藏高原西北部的和田、阿里地區(qū)以及巴音郭楞蒙古自治州的部分地區(qū);而仲巴縣、改則縣等區(qū)縣所圍成的區(qū)域年際變化較小。青藏高原地區(qū)降雨侵蝕力R因子的值在0~15 532 MJ·mm/(hm2·h·a)之間,均值為648.84 MJ·mm/(hm2·h·a),R值自東向西遞減,呈現(xiàn)明顯的條帶狀分布特征,其中在位于青藏高原的南部地區(qū)形成了一個明顯的高值中心。② 影響力探測結(jié)果表明,最低氣溫的解釋力最高,q值為0.453,主導降雨侵蝕力R因子的時空分布特征。③ 交互探測結(jié)果表明,青藏高原各因子間協(xié)同效應(yīng)均表現(xiàn)為雙因子增強或非線性增強關(guān)系,較單因子作用有明顯提升。氣候因子間的交互作用對降雨侵蝕力R因子的影響較為突出。

    關(guān)? 鍵? 詞:降雨侵蝕力因子; 時空分布; 變差系數(shù); 地理探測器; 青藏高原

    中圖法分類號: S157

    文獻標志碼: A

    DOI:10.16232/j.cnki.1001-4179.2024.06.012

    0? 引 言

    在全球氣候變化下,青藏高原大多數(shù)地區(qū)極端天氣頻發(fā),降雨量、降雨強度增加,侵蝕性降雨增多,加劇了土壤侵蝕。降雨侵蝕力R因子是評價這種潛在能力的一個動力指標[1]。降雨是引起土壤侵蝕的主要動力因素,其數(shù)值能表征降雨對土壤侵蝕潛力的大?。?-4]。降雨侵蝕力因子自提出并用于美國通用土壤流失方程(USLE)以來,一直是土壤侵蝕研究的重點和熱點之一[5-6]。降雨侵蝕力R因子值的大小與分布受氣候條件、地形條件和人類活動的共同影響,因此確定自然因素與人類活動對降雨侵蝕的解釋力對制定水土保持措施尤為重要[7-8]。

    目前國內(nèi)外學者對降雨侵蝕力R因子的研究集中在時空變化規(guī)律以及計算方法上,例如Johannsen等[9]利用171個氣象站的降雨數(shù)據(jù),計算了奧地利各地1995~2015年降雨侵蝕力,并根據(jù)降雨侵蝕力的時空分布差異反映侵蝕危險性;馬道銘等[10]基于日降水數(shù)據(jù)研究山西省近60 a降雨侵蝕力時空變化,結(jié)果表明,該地區(qū)大部分降雨侵蝕力呈下降趨勢,黃河

    沿線呈上升趨勢。張悅等[11]基于2種降雨侵蝕力模型分析了近40 a黃河中游降雨和降雨侵蝕力時空演變,結(jié)果表明,黃河中游年均降雨侵蝕力受地形和地貌影響,具有顯著的垂直空間差異。陳劍南等[12]分析黃土高原降雨侵蝕力時空變化,結(jié)果表明,該地區(qū)降雨侵蝕力空間存在顯著差異,但變化不顯著。然而,目前針對降雨侵蝕力R因子驅(qū)動因素的相關(guān)研究相對較少,因此本文以青藏高原為研究對象,分析降雨侵蝕力R因子的分布特征和影響因素。

    青藏高原地處中國西部,是世界上海拔最高的高原,也是全球陸地上平均海拔最高、環(huán)境最脆弱的高原。降雨作為青藏高原生態(tài)系統(tǒng)中至關(guān)重要的氣候因素之一,對土壤侵蝕具有重要影響。青藏高原降雨侵蝕力R因子的分布規(guī)律受多種因素的共同影響,包括氣候變化、地形起伏和人類活動等。這些因素之間存在著復雜的相互作用關(guān)系,需要進行深入研究以揭示其內(nèi)在規(guī)律[13]。隨著全球氣候變化的不斷加劇,青藏高原降雨侵蝕力的分布特征及其驅(qū)動因素備受關(guān)注,故降雨侵蝕力的研究對于深入了解該地區(qū)土壤侵蝕過程,保護生態(tài)環(huán)境,以及合理利用土地資源具有重要意義。

    本文基于地理探測器模型,利用遙感數(shù)據(jù)和地理信息系統(tǒng)技術(shù),系統(tǒng)分析青藏高原地區(qū)降雨侵蝕力R因子分布特征及其驅(qū)動因素,不僅可以為該地區(qū)土壤侵蝕的定量評估提供科學依據(jù),還可以為生態(tài)環(huán)境保護和可持續(xù)發(fā)展提供重要決策支持。

    1? 材料和方法

    1.1? 研究區(qū)概況

    青藏高原地處26°00′12″N~39°46′50″N,73°18′52″E~104°46′59″E之間,是全球海拔最高的地區(qū)之一,平均海拔超過4 000 m。青藏高原主要包括昆侖山脈、喜馬拉雅山脈、唐古拉山脈等。青藏高原位于中國西南部,行政區(qū)域包括青海省、甘肅省、四川省、西藏自治區(qū)、云南省和新疆維吾爾自治區(qū),如圖1所示。青藏高原多年平均降水量為300~500 mm,降水量在不同區(qū)域之間存在很大的差異;青藏高原中部和西部地區(qū)則相對較干旱,這些地區(qū)大多屬于內(nèi)陸沙漠或荒漠地帶;青藏高原東南部的云南地區(qū)年降水量最多,約為1 200~2 500 mm;而青藏高原北部和西部地區(qū)的年降水量相對較少,只有100~300 mm。李文君

    等[14]研究表明喜馬拉雅山脈北翼年降水量不足600 mm,而喜馬拉雅山脈南翼為亞熱帶及熱帶北緣山地森林氣候,年降水量為1 000~4 000 mm。昆侖山脈南翼屬高寒半荒漠和荒漠氣候,年降水量為20~100 mm。夏季和秋季是青藏高原上降水最多的季節(jié),而冬季和春季則相對較少。

    1.2? 數(shù)據(jù)來源及預處理

    1.2.1? 數(shù)據(jù)來源

    本研究所需數(shù)據(jù)包括降雨數(shù)據(jù)、數(shù)字高程數(shù)據(jù)、人口密度、土地利用和氣象數(shù)據(jù)等,空間分辨率約為1 km。詳細見表1。

    1.2.2? 數(shù)據(jù)預處理

    本研究下載的數(shù)據(jù)通過R語言將其轉(zhuǎn)化為tif格式,并將其整合為年值數(shù)據(jù)。運用ArcGIS 10.2對以上數(shù)據(jù)再次進行按掩膜提取、投影轉(zhuǎn)換和重分類等處理,確保所有影像與青藏高原邊界空間參考一致和像元大小相同。坡度和坡向數(shù)據(jù)運用ArcGIS 10.2由DEM數(shù)據(jù)計算得來。

    1.3? 研究方法

    1.3.1? 地理探測器

    地理探測器是一種統(tǒng)計學方法[15],它基于空間分層異質(zhì)性,通過比較自變量與因變量的空間分布相似性來揭示驅(qū)動這些變量的因素。其核心思想在于假設(shè)如果某個自變量對因變量有重要影響,那么它們的空間分布格局會具有相似性,從而用來度量自變量對因變量的貢獻率[16]。這種方法已經(jīng)在環(huán)境科學[17]、社會經(jīng)濟[18]和人類健康[19]等多個領(lǐng)域得到廣泛應(yīng)用。本文主要探討因子探測和交互探測這2個模塊,揭示不同探測因子對青藏高原地區(qū)降雨侵蝕力R因子的驅(qū)動機制。

    本文選取坡向(X1)、坡度(X2)、高程(X3)、 最高氣溫(X4)、 最低氣溫(X5)、平

    均氣溫(X6)、人口密度(X7)和土地利用類型(X8)共8個因子來分析其對降雨侵蝕力R因子的影響。

    (1) 因子探測。

    分析因變量Y在空間上的變化特征,并確定各個自變量X在多大程度上解釋了Y的空間變化。公式如下:q=1-1Nσ2Lh=1Nhσ2h

    (1)

    式中:L為因變量Y或探測因子X的分層;Nh和σ2h分別為層h的單元數(shù)和Y值方差,N和σ2分別為研究區(qū)的總單元數(shù)和Y值方差。q為[0,1],其值越大,表示探測因子X對降雨侵蝕力R因子的解釋力越強,反之則越弱。

    (2) 交互作用探測。

    識別不同影響因子X之間的交互作用,評估對因變量Y的解釋力。交互作用可分為表2中的5類。

    1.3.2? 降雨侵蝕力R因子計算

    R因子是衡量降雨侵蝕能力的定量指標,它用于測量降雨造成分離和傳輸土壤顆粒的能力。本文選用修訂的R因子公式[20]計算1991~2020年的逐年降雨侵蝕力,計算公式為

    R=0.44488P0.969882

    (2)

    式中:P為年降雨量,mm。

    1.3.3? 變差系數(shù)

    年降雨侵蝕力R因子的變差系數(shù)Cv表示為

    Cv=R年s/R年=ni=1(R年i-R年)2n-1R年

    (3)

    式中:Cv為變差系數(shù),反映不同年限長度降雨侵蝕力R因子的穩(wěn)定性;R年s為某序列年限長度下的降雨侵蝕力R因子的標準差,MJ·mm/(hm2·h·a);R年i為第i年的降雨侵蝕力R值;R年為相應(yīng)年限長度下降雨侵蝕力R的平均值,MJ·mm/(hm2·h·a);n為序列年限長度。Cv值越大,降雨侵蝕力R因子的年際變化越劇烈,對水土流失治理工程不利;Cv值越小,降雨侵蝕力R因子的年際變化越小,越有利于水土流失治理工程規(guī)劃實施。

    2? 結(jié)果與分析

    2.1? 青藏高原降雨侵蝕力R因子的時空分布

    2.1.1? 年際變化

    降雨侵蝕力因子的年際變化特征用3種指標表示:變差系數(shù)Cv、最大年R值與年平均R值的比值(R年max/R年)、最大年R值與最小年R值的比值(R年max/R年min)。降雨侵蝕力R因子及變差系數(shù)Cv空間分布如圖2所示。青藏高原年均R值年際變差系數(shù)介于0~0.78之間,均值為0.26。最大年R值一般為平均年R值的1.14倍,最大年R值為最小年R值的1.57倍,表明降雨侵蝕力R因子的年際變化較大。其中,年際變化較大的區(qū)域分布在青藏高原西北部的和田、阿里地區(qū)以及巴音郭楞蒙古自治州的部分地區(qū),這些區(qū)域的Cv值都超過了0.5。仲巴縣、改則縣、雙湖等區(qū)縣所圍成的區(qū)域內(nèi)R值的年際變化較小,Cv值平均在0.3左右。

    2.1.2? 空間分布

    青藏高原1991~2020年降雨侵蝕力R因子的空間分布如圖3所示。青藏高原地區(qū)降雨侵蝕力R值分布區(qū)間為0~15 532 MJ·mm/(hm2·h·a),多年平均R值648.84 MJ·mm/(hm2·h·a)。從整體上看,降雨侵蝕力R值自東向西遞減,呈現(xiàn)明顯的條帶狀分布特征,其中在位于青藏高原地區(qū)南部的林芝、山南地區(qū)形成了一個明顯的高值中心。四川省西部的阿壩、甘孜等地以及青海省東南部的部分地區(qū)降雨侵蝕力R值較大,其變化范圍為3 487~15 532 MJ·mm/(hm2·h·a)。青藏高原除南部地區(qū)外,降雨侵蝕力R值均較小,其變化范圍為0~1 284 MJ·mm/(hm2·h·a),且低值區(qū)域自東向西逐漸減小,呈條帶狀分布,那曲、阿里的西北部區(qū)域位于明顯的地質(zhì)抬升區(qū),部分地區(qū)處于降雨陰影區(qū),年平均降雨量較少,因此R值也相對較小。

    2.2? 青藏高原降雨侵蝕力R因子動態(tài)變化

    為了更加深入地了解青藏高原地區(qū)降雨侵蝕力R因子的動態(tài)變化,對1991~2020年這一時間段內(nèi),每5 a選取一個代表年份,將所選年份的降雨侵蝕力R因子空間分布分別進行對比。如圖4所示,1991~

    2020年青藏高原地區(qū)整體降雨侵蝕力R值變化不是

    很明顯,整體呈現(xiàn)東南低西北高的空間分布特征。青

    藏高原大部分地區(qū)R值均較小,而位于青藏高原

    南部地區(qū)的林芝、山南地區(qū),R值相對較大。其變化區(qū)

    域主要分布在青藏高原東南部的甘孜、阿壩和昌都地

    區(qū),其中1992年降雨侵蝕力R值的分布相較于其他年份變化較大,往后其他年份R值的變化較小。

    2.3? 青藏高原降雨侵蝕力影響因素地理探測

    2.3.1? 影響因素影響力探測(因子探測)

    因子探測器結(jié)果表示各因子對降雨侵蝕力R值的影響,計算結(jié)果q值為降雨侵蝕力R因子的解釋力。由表3可知,統(tǒng)計結(jié)果顯示,各個因子對降雨侵蝕力R值的影響不同,各個因子的解釋力依次為最低氣溫(X5)(0.453)>最高氣溫(X4)(0.398)>平均氣溫(X6)(0.396)>高程(X3)(0.167)>土地利用(X8)(0.150)>人口密度(X7)(0.082)>坡度(X2)(0.058)>坡向(X1)(0.004),人口密度和坡度、坡向作用微小。最低氣溫的q值最高,表明氣溫對降雨侵蝕力R值的貢獻最大。氣溫的解釋能力先增強后降低,2005年最低氣溫對降雨侵蝕力R值的作用程度達到最高,因此,溫度是影響降雨侵蝕力R值空間分異的主要驅(qū)動因子,而人口密度和坡度、坡向3個影響因子對降雨侵蝕力R值的解釋力微小。故降雨侵蝕力R值的大小與分布受氣候條件影響最大。

    從動態(tài)變化來看,各因子對青藏高原地區(qū)降雨侵蝕力R因子的解釋力呈現(xiàn)不同變化趨勢,如圖5所示。

    在1995~2020年近30 a期間,分別為最高氣溫、最低氣溫和年均氣溫q值增加,而坡度、人口密度和土地利用q值降低??傮w來看,30 a間最高氣溫、最低氣溫、年均氣溫和海拔對研究區(qū)降雨侵蝕力因子的影響呈增加趨勢,其中最高氣溫q值增加速率最高,為0.002 1 a-1,這可能是由于全球變暖,該區(qū)域的氣溫上升,更容易形成暖濕氣流,使降水量增加,從而使降雨侵蝕力增加;而人口密度對研究區(qū)降雨侵蝕力的影響趨于減少,這可能是由于該地區(qū)近年來受到“集中搬遷”政策或者人口大面積遷移的影響;其余各因子隨時間變化不明顯(速率<0.001 a-1)。

    2.3.2? 影響因素交互作用探測

    各個因子交互作用的結(jié)果如圖6所示,大多數(shù)因子交互作用的q值大于單個因素的q值,并且因子相互作用類型為非線性增強和雙因子增強,不存在獨立關(guān)系。具體而言,其他因子與最低溫度交互作用q值均大于單因子解釋力。此外,X8∩X5(0.529)>X5∩X4(0.500)>X8∩X4(0.497)>X8∩X6(0.485)>X7∩X5(0.485),結(jié)果表明,土地利用與最低氣溫、最低氣溫與最高氣溫、土地利用與最高氣溫、土地利用與氣溫、人口密度與最低溫度的交互作用呈現(xiàn)雙變量增強效應(yīng);X8∩X5(0.529)是因子交互探測結(jié)果中的最大值,這表明土地利用與最高溫度兩因子交互作用呈現(xiàn)最顯著的雙變量增強效應(yīng),對降雨侵蝕力R因子的作用最大。由表3可知,X2、X3、X4、X5、X6、X8的p值為0(p<0.05),表明以上因子與降雨侵蝕力R顯著相關(guān),相反X1、X7的p值大于0.05,故坡向和人口密度對降雨侵蝕力R因子的影響作用不顯著。

    因此,在土壤侵蝕治理中,應(yīng)特別關(guān)注溫度因子,并針對溫度與其他因子的交互作用制定合理的治理措施,以達到最佳的治理效果。此外,坡度與坡向兩個因子單獨作用對土壤侵蝕的解釋力比較小,但兩個因子與其他因素交互作用時對土壤侵蝕的解釋力有所增加,表明自然因素在人類活動的影響下對土壤侵蝕的作用更加劇烈。

    3? 討 論

    本研究采用地理探測器模型,對青藏高原地區(qū)降雨侵蝕力R因子的分布規(guī)律與驅(qū)動因素進行了深入分析。因子探測分析結(jié)果表明氣溫、高程和對降雨侵蝕力R因子的作用影響最強,在全球氣候變暖的大環(huán)境下,青藏高原大多數(shù)地區(qū)強降雨事件增加,侵蝕性降雨增多,導致R值呈上升的趨勢[21],但R值并非只受侵蝕性降雨的影響,同時還受氣候條件、人類活動、地形等條件的影響。青藏高原位于喜馬拉雅山脈和昆侖山脈之間,其地形復雜多樣,包括高山、高原、河谷等地貌類型。這種多樣性使得降雨侵蝕力R因子在時空上表現(xiàn)出較大的差異性。青藏高原多年平均R值空間上總體呈東南高西北低分布特點,而高值區(qū)集中分布在南部地區(qū),可能是由于林芝市、山南市位于喜馬拉雅山南坡一帶,夏季氣候海洋性特征顯著,大量來自海洋的暖濕氣流,經(jīng)地形抬升形成充沛降水,具有暖季降雨頻次高、強度大的特點,導致R值分配極度不均。Guo等[22]建立了全國降雨侵蝕力重力模型,也發(fā)現(xiàn)國家尺度上降雨侵蝕力呈現(xiàn)東南高西北低的特點。這一趨勢可能與夏季風從東南向西北推進的方向有關(guān)[23]。而位于西北部的喀什地區(qū)位于青藏高原西北部,地處群山之間,遠離海洋,屬暖溫帶大陸性干旱氣候帶,太陽輻射強烈,使得來自大西洋的盛行西風氣流以及來自北冰洋的冷濕氣流,太平洋和印度洋的季風難以到達,氣候干燥,降水稀少,故R值較低。這與鄒玉霞等[24]在1960~2017年重慶市不同量級降雨侵蝕力R值的時空變化特征的結(jié)果一致。其中1992年降雨侵蝕力R值的分布相較于其他年份變化較大,往后其他年份R值的變化較小,可能是因為提出了關(guān)于生態(tài)文明建設(shè)的決策部署,包括在北方地區(qū)實施退牧還草和退耕還林工程,加強林地植被管理和修復,增加植被覆蓋度,防治過度放牧和開墾田地而造成草場退化,綜合防治農(nóng)牧交錯地帶土壤侵蝕[25]。

    因此,土壤侵蝕治理不僅要考慮RUSLE模型中的單個影響因子,還要考慮其他因素以及其他影響因子,以制定更全面、科學的土壤侵蝕治理策略,從而最大程度地減少土壤侵蝕對該地區(qū)生態(tài)環(huán)境和可持續(xù)發(fā)展的不利影響,更好地為認識青藏高原的土壤侵蝕機制、生態(tài)環(huán)境演變規(guī)律及可持續(xù)發(fā)展戰(zhàn)略提供科學依據(jù)。

    4? 結(jié) 論

    (1) 青藏高原年際變化表現(xiàn)出明顯的“由東向西,自南向北”遞增的趨勢,表明東南向西北方向上降雨侵蝕力的穩(wěn)定性逐漸減弱。降雨侵蝕力R值年際變差系數(shù)介于0~0.78之間,均值為0.26。最大年R值一般為平均年R值的1.14倍,最大年R值為最小年R值的1.57倍,表明降雨侵蝕力R因子年際變化較大。降雨侵蝕力因子R值呈東南高西北低的特點,變化范圍在0.325~15 532 MJ·mm/(hm2·h)之間,均值為648.84 MJ·mm/(hm2·h),R值在空間分布上整體自東向西遞減,呈現(xiàn)明顯的條帶狀分布特征,其中在位于青藏高原南部的林芝、山南地區(qū)形成了一個明顯的高值中心,其他大部分地區(qū)R值均較小。

    (2) 影響力探測結(jié)果表明,各個因子的解釋力依次為最低氣溫(0.453)>最高氣溫(0.398)>平均氣溫(0.396)>高程(0.167)>土地利用(0.150)>人口密度(0.082)>坡度(0.058)>坡向(0.004),人口密度和坡度、坡向作用微小。最低氣溫的q值最高,表明其對降雨侵蝕力R值的貢獻最大。

    (3) 交互探測結(jié)果表明,青藏高原地區(qū)各因子間協(xié)同效應(yīng)均表現(xiàn)為雙因子增強或非線性增強關(guān)系,較單因子作用有明顯提升。最低氣溫∩土地利用(0.529)是因子交互探測結(jié)果中的最大值,這表明最低溫度與土地利用兩因子交互作用呈現(xiàn)最顯著的雙變量增強效應(yīng),對降雨侵蝕力R因子的作用最大,且氣候因子間及氣候因子∩土地利用的交互作用對降雨侵蝕力R因子的影響較為突出。

    參考文獻:[1]? RAN Q,SU D,LI P,et al.Experimental study of the impact of rainfall characteristics on runoff generation and soil erosion[J].Journal of Hydrology,2012,424:99-111.

    [2]? 金建君,謝云,張科利.不同樣本序列下侵蝕性雨量標準的研究[J].水土保持通報,2001,21(2):31-33.

    [3]? 謝云,劉寶元,章文波.侵蝕性降雨標準研究[J].水土保持學報,2000,14(4):6-11.

    [4]? 王萬忠,焦菊英.中國的土壤侵蝕因子定量評價研究[J].水土保持通報,1996,16(5):1-20.

    [5]? 馬良,左長清,孫勐,等.山東省降雨侵蝕力空間分布特征及其簡易方程的研究[J].水土保持研究,2010,17(2):28-31.

    [6]? 章文波,付金生.不同類型雨量資料估算降雨侵蝕力[J].資源科學,2003,25(1):35-41.

    [7]? 田培,許盈,楊偉,等.鄂東南水土流失時空變化及其影響因子定量評價[J].人民長江,2023,54(7):61-68.

    [8]? 吳明洋.成都市降雨量及降雨侵蝕力年內(nèi)分布特征研究[J].人民長江,2019,50(4):81-85.

    [9]? JOHANNSEN E L,SCHMALTZ E M,MITROVITS O,et al.An update of the spatial and temporal variability of rainfall erosivity (R-factor) for the main agricultural production zohes of Austria[J].Catena,2022,215:106305.

    [10]馬道銘,回毅瀅,趙炯昌,等.近60年山西省降雨侵蝕力時空變化特征[J].水土保持研究,2023,30(6):24-31,48.

    [11]張悅,張艷,史飛航,等.基于不同模型的黃河中游降雨侵蝕力時空變化分析[J].水土保持通報,2023,43(2):220-229,238.

    [12]陳劍南,劉益麟,李朋飛,等 1901~2016年黃土高原降雨侵蝕力時空變化[J].水土保持研究,2022,29(4):39-46.

    [13]姚昆,張存杰,李玉珍,等.四川省安寧河流域土壤侵蝕時空演變及預測分析[J].人民長江,2021,52(7):65-70.

    [14]李文君,李鵬,封志明,等.基于人居環(huán)境特征的青藏高原“無人區(qū)”空間界定[J].地理學報,2021,76(9):2118-2129.

    [15]王勁峰,徐成東.地理探測器:原理與展望[J].地理學報,2017,72(1):116-134.

    [16]HUO H,SUN C P.Spatiotemporal variation and influencing factors of vegetation dynamics based on Geodetector:a case study of the northwestern Yunnan Plateau,China[J].Ecological Indicators,2021,130:108005.

    [17]徐勇,黃雯婷,鄭志威,等.基于空間尺度效應(yīng)的西南地區(qū)植被NPP影響因子探測[J].環(huán)境科學,2023,44(2):900-911.

    [18]POLYKRETIS C,ALEXAKIS D D.Spatial stratified heterogeneity of fertility and its association with socio-economic determinants using Geographical Detector:the case study of Crete Island,Greece[J].Applied Geography,2021,127:102384.

    [19]ZHANG X X,XU C D,XIAO G X.Spatial heterogeneity of the association between temperature and hand,foot,and mouth disease risk in metropolitan and other areas[J].Science of the Total Environment,2020,713:136623.

    [20]楊子生.云南省金沙江流域土壤流失方程研究[J].山地學報,2002,20(增1):1-9.

    [21]鄭哲,鄒進,潘鋒,等.1961~2010年怒江流域降雨侵蝕力的時空變化特征[J].中國水土保持科學,2018,16(5):16-22.

    [22]GUO B,LUO W,HAN B,et al.Multiple spatiak temporal scale change patterns of rainfall erosivity in China over past 58 years based on gravity centre model [J].Geomatics,Natural Hazards and Risk,2019,10(1):2200-2219.

    [23]LIU S,HUANG S,XIE Y,et al.Spatial temporal changes of rainfall erosivity in the Toess Plateau,China:changing patterns,causes and implications[J].Catena,2018,166:279-289.

    [24]鄒玉霞,樊國中,劉淑婧,等.1960~2017年重慶市不同量級降雨侵蝕力R值的時空變化特征[J].水土保持學報,2022,36(6):41-48.

    [25]曹前發(fā).“綠水青山就是金山銀山”的實踐《建設(shè)美麗中國:新時代生態(tài)文明建設(shè)理論與實踐》第十章[J].中國林業(yè)產(chǎn)業(yè),2023(7):82-96.

    (編輯:黃文晉)

    Spatiotemporal distribution and driving factors of rainfall erosivity R factor of Qinghai-Tibet Plateau

    SU Libin1,2,WANG Maolin3,YUAN Chunhui3,KE Haocheng3

    (1.Research Center of Civil,Hydraulic and Power Engineering of Tibet,Linzhi 860000,China;

    2.College of Water Resources and Civil Engineering,Tibet Agricultural and Animal Husbandry University,Linzhi 860000,China;

    3.School of Energy and Power Engineering,Lanzhou University of Technology,Lanzhou 730050,China)

    Abstract:

    The rainfall erosivity R is an important factor for calculating soil erosion,and its accurate calculation is of great significance for the dynamic monitoring of soil erosion.This paper uses remote sensing data and geographic information system (GIS) technology,along with the geographic detector model,to analyze the distribution characteristics and influencing factors of the rainfall erosivity R in the Qinghai-Tibet Plateau from 1991 to 2020.The research results show that:① The interannual variation of the R value in the Qinghai-Tibet Plateau generally exhibited a high in the northwest and a low in the southeast distribution pattern,with a coefficient of variation ranging from 0 to 0.78.Regions with large interannual variations were mainly distributed in the northwest part of the Qinghai-Tibet Plateau,including Hotan,Ngari,and some areas of Bayingolin Mongol Autonomous Prefecture,while regions such as Zhamba County and Gar County had smaller interannual variations.The values of the rainfall erosivity R in the Qinghai-Tibet Plateau ranged from 0 to 15 532 MJ·mm/(hm2·h·a),with a mean value of 648.84 MJ·mm/(hm2·h·a).The R value decreased from east to west,showing a distinctly banded distribution pattern.A significant high-value center was found in the southern part of the Qinghai-Tibet Plateau.② The influence detection results indicated that the minimum temperature has the highest explanatory power,with a q value of 0.453,dominating the spatiotemporal distribution characteristics of the rainfall erosivity R.③ The interaction detection results showed that the synergistic effects among various factors on the Qinghai-Tibet Plateau exhibited either a two-factor enhancement or a nonlinear enhancement relationship,significantly surpassing the single-factor effects.The interaction effect among climate factors had a prominent impact on the rainfall erosivity R.

    Key words:

    rainfall erosivity factor; spatiotemporal distribution; coefficient of variation; geographic detector; Qinghai-Tibet Plateau

    猜你喜歡
    年際土壤侵蝕青藏高原
    青藏高原上的“含羞花”
    北緯30°中層頂區(qū)域鈉與鐵原子層的結(jié)構(gòu)和年際變化
    為了讓青藏高原的天更藍、水更綠、草原更美
    鄉(xiāng)村聚落土壤侵蝕環(huán)境與水土流失研究綜述
    海壇島土壤侵蝕問題研究
    大別山區(qū)土壤侵蝕動態(tài)變化及趨勢預測
    亞洲夏季風的年際和年代際變化及其未來預測
    與北大西洋接壤的北極海冰和年際氣候變化
    南水北調(diào)中線水源地土壤侵蝕經(jīng)濟損失估算
    化石見證青藏高原的隆起
    生物進化(2014年2期)2014-04-16 04:36:27
    久久人妻熟女aⅴ| 免费在线观看黄色视频的| 校园春色视频在线观看| 午夜福利影视在线免费观看| 午夜两性在线视频| 国产蜜桃级精品一区二区三区 | 久久久久精品国产欧美久久久| 亚洲专区字幕在线| 18禁观看日本| 成人免费观看视频高清| 久久午夜亚洲精品久久| 看免费av毛片| 国产成人欧美| 国产午夜精品久久久久久| 1024视频免费在线观看| 久久ye,这里只有精品| 欧美激情久久久久久爽电影 | 飞空精品影院首页| 日本精品一区二区三区蜜桃| 这个男人来自地球电影免费观看| 欧美大码av| 精品福利观看| av视频免费观看在线观看| 69av精品久久久久久| 十分钟在线观看高清视频www| 日日摸夜夜添夜夜添小说| 成人特级黄色片久久久久久久| 国产精品久久久人人做人人爽| 狠狠狠狠99中文字幕| 两人在一起打扑克的视频| 91精品国产国语对白视频| 成人特级黄色片久久久久久久| 国产主播在线观看一区二区| 亚洲人成77777在线视频| 天堂中文最新版在线下载| 成人手机av| 国产国语露脸激情在线看| 精品少妇久久久久久888优播| 亚洲欧美激情在线| 99精品欧美一区二区三区四区| 久久精品成人免费网站| 制服人妻中文乱码| 变态另类成人亚洲欧美熟女 | 大型黄色视频在线免费观看| 国产免费av片在线观看野外av| 啪啪无遮挡十八禁网站| 欧洲精品卡2卡3卡4卡5卡区| 在线观看免费视频日本深夜| 18禁国产床啪视频网站| 香蕉国产在线看| 国产精品一区二区精品视频观看| 精品电影一区二区在线| 一本一本久久a久久精品综合妖精| 久久精品亚洲av国产电影网| 国产片内射在线| 水蜜桃什么品种好| 宅男免费午夜| 国产欧美日韩一区二区精品| 丝袜美足系列| 村上凉子中文字幕在线| 性色av乱码一区二区三区2| 久久久精品国产亚洲av高清涩受| 王馨瑶露胸无遮挡在线观看| 免费少妇av软件| 国产有黄有色有爽视频| 一级片'在线观看视频| 亚洲国产精品合色在线| 国产欧美日韩一区二区三| 亚洲一区二区三区欧美精品| 精品久久久久久久久久免费视频 | 精品亚洲成a人片在线观看| av有码第一页| 女人被狂操c到高潮| 国产亚洲精品一区二区www | xxxhd国产人妻xxx| 老熟女久久久| а√天堂www在线а√下载 | 日韩三级视频一区二区三区| 欧美日韩一级在线毛片| 国产亚洲精品一区二区www | 在线观看免费高清a一片| 黑人巨大精品欧美一区二区蜜桃| 久久中文字幕一级| 国产亚洲欧美精品永久| 又紧又爽又黄一区二区| 在线免费观看的www视频| 久久精品91无色码中文字幕| 搡老乐熟女国产| 一级a爱视频在线免费观看| 亚洲精品一二三| 18禁裸乳无遮挡免费网站照片 | 精品亚洲成国产av| 国产成人精品无人区| 一级片'在线观看视频| 9热在线视频观看99| 亚洲片人在线观看| 久9热在线精品视频| 黄色怎么调成土黄色| ponron亚洲| 色94色欧美一区二区| 精品一区二区三区四区五区乱码| 又黄又爽又免费观看的视频| 中文字幕人妻丝袜一区二区| 免费在线观看亚洲国产| 国产精品99久久99久久久不卡| 亚洲视频免费观看视频| 欧美日韩一级在线毛片| 成年人午夜在线观看视频| avwww免费| 最近最新免费中文字幕在线| 窝窝影院91人妻| 18禁观看日本| 亚洲aⅴ乱码一区二区在线播放 | 身体一侧抽搐| 丝袜美腿诱惑在线| www.熟女人妻精品国产| 999久久久国产精品视频| 人人妻人人澡人人爽人人夜夜| 免费不卡黄色视频| 少妇 在线观看| 国产精品免费一区二区三区在线 | 亚洲专区国产一区二区| 高清视频免费观看一区二区| 69精品国产乱码久久久| 亚洲专区字幕在线| 人妻一区二区av| 在线观看日韩欧美| 国产精品国产高清国产av | 99国产精品99久久久久| av欧美777| 一边摸一边抽搐一进一小说 | 国产精品av久久久久免费| 狂野欧美激情性xxxx| 9热在线视频观看99| 免费在线观看完整版高清| 亚洲色图 男人天堂 中文字幕| 久久青草综合色| 韩国精品一区二区三区| 人成视频在线观看免费观看| 精品少妇久久久久久888优播| 嫩草影视91久久| 黑丝袜美女国产一区| 婷婷丁香在线五月| 老汉色∧v一级毛片| 久久精品人人爽人人爽视色| 正在播放国产对白刺激| 亚洲伊人色综图| 亚洲av第一区精品v没综合| 午夜成年电影在线免费观看| 18禁观看日本| 久久精品成人免费网站| 欧美日韩国产mv在线观看视频| 一级a爱片免费观看的视频| 午夜福利在线观看吧| 丰满迷人的少妇在线观看| 三上悠亚av全集在线观看| 99国产综合亚洲精品| 欧美另类亚洲清纯唯美| 亚洲av成人不卡在线观看播放网| 亚洲五月天丁香| 我的亚洲天堂| 国产不卡一卡二| 欧美亚洲 丝袜 人妻 在线| 啦啦啦在线免费观看视频4| videos熟女内射| 交换朋友夫妻互换小说| 中文字幕av电影在线播放| 男女午夜视频在线观看| 亚洲黑人精品在线| 久久久久久久久免费视频了| 极品少妇高潮喷水抽搐| 国产精品偷伦视频观看了| 看片在线看免费视频| 亚洲成人免费av在线播放| 久久亚洲真实| 99精品在免费线老司机午夜| 天天影视国产精品| 久99久视频精品免费| 热99国产精品久久久久久7| 1024香蕉在线观看| 欧美丝袜亚洲另类 | 91字幕亚洲| 国产亚洲精品久久久久5区| 男女高潮啪啪啪动态图| 国产男靠女视频免费网站| 亚洲黑人精品在线| 91九色精品人成在线观看| 精品人妻在线不人妻| 国产熟女午夜一区二区三区| 成年版毛片免费区| 久久精品国产综合久久久| 亚洲中文字幕日韩| 99精国产麻豆久久婷婷| 高潮久久久久久久久久久不卡| 午夜免费成人在线视频| 亚洲熟女精品中文字幕| 亚洲精品一卡2卡三卡4卡5卡| 成年人免费黄色播放视频| 777米奇影视久久| 精品国产超薄肉色丝袜足j| 叶爱在线成人免费视频播放| 50天的宝宝边吃奶边哭怎么回事| tube8黄色片| 国产成人影院久久av| 一级毛片高清免费大全| 日韩欧美三级三区| 亚洲综合色网址| 不卡一级毛片| 色精品久久人妻99蜜桃| av电影中文网址| 99热只有精品国产| 老司机午夜福利在线观看视频| 国产无遮挡羞羞视频在线观看| 99国产综合亚洲精品| 久久热在线av| 精品人妻在线不人妻| 午夜老司机福利片| 成年版毛片免费区| 日韩欧美免费精品| 亚洲欧美日韩高清在线视频| 精品一区二区三卡| 在线观看免费午夜福利视频| 久热爱精品视频在线9| 国产一区二区三区在线臀色熟女 | 欧美日韩av久久| 女性生殖器流出的白浆| 最近最新中文字幕大全电影3 | 国产熟女午夜一区二区三区| 国产亚洲av高清不卡| 久久人妻福利社区极品人妻图片| 夜夜躁狠狠躁天天躁| 在线观看一区二区三区激情| 国产国语露脸激情在线看| 侵犯人妻中文字幕一二三四区| 亚洲av电影在线进入| 国产片内射在线| 91成人精品电影| 在线观看免费日韩欧美大片| 老熟妇乱子伦视频在线观看| 丝袜美腿诱惑在线| 丝袜美足系列| 黄色视频不卡| 50天的宝宝边吃奶边哭怎么回事| 久久天堂一区二区三区四区| 亚洲精品国产一区二区精华液| 大型黄色视频在线免费观看| 亚洲熟女精品中文字幕| 午夜两性在线视频| 色精品久久人妻99蜜桃| 美女国产高潮福利片在线看| 欧美日韩中文字幕国产精品一区二区三区 | 国产男靠女视频免费网站| 大香蕉久久成人网| 97人妻天天添夜夜摸| 老司机影院毛片| 黄片播放在线免费| 一进一出好大好爽视频| 国产精品久久电影中文字幕 | 一区在线观看完整版| 9热在线视频观看99| 精品久久久久久久毛片微露脸| 久久久久精品国产欧美久久久| 午夜福利视频在线观看免费| 日本wwww免费看| 国产亚洲av高清不卡| 18禁美女被吸乳视频| 看片在线看免费视频| 在线免费观看的www视频| 欧美成人免费av一区二区三区 | 男人的好看免费观看在线视频 | 在线观看免费视频网站a站| 精品久久久久久电影网| 国产欧美日韩一区二区三区在线| 成人黄色视频免费在线看| 不卡av一区二区三区| 亚洲中文av在线| 法律面前人人平等表现在哪些方面| av中文乱码字幕在线| 国产精品成人在线| 日韩中文字幕欧美一区二区| 又大又爽又粗| 99国产精品一区二区三区| 亚洲av成人av| 久久国产精品人妻蜜桃| 丰满的人妻完整版| 美女扒开内裤让男人捅视频| 国产又爽黄色视频| 中出人妻视频一区二区| 国产一区有黄有色的免费视频| 久久久国产成人免费| 少妇 在线观看| 亚洲欧美一区二区三区黑人| 亚洲自偷自拍图片 自拍| 国产一区二区激情短视频| 丰满迷人的少妇在线观看| 在线观看免费视频网站a站| 最新在线观看一区二区三区| 国产深夜福利视频在线观看| 99精品久久久久人妻精品| 黄网站色视频无遮挡免费观看| 51午夜福利影视在线观看| 99国产精品99久久久久| 亚洲中文av在线| 精品一区二区三区av网在线观看| 黄片小视频在线播放| 欧美日韩亚洲国产一区二区在线观看 | 大香蕉久久网| 国产精品电影一区二区三区 | 在线国产一区二区在线| 巨乳人妻的诱惑在线观看| 在线观看日韩欧美| 纯流量卡能插随身wifi吗| 欧美精品一区二区免费开放| 久久久久久免费高清国产稀缺| 老鸭窝网址在线观看| 中文字幕精品免费在线观看视频| 中文字幕人妻丝袜制服| 老司机深夜福利视频在线观看| av有码第一页| 一级毛片高清免费大全| 成人亚洲精品一区在线观看| 国产成人精品久久二区二区免费| 亚洲av日韩精品久久久久久密| 欧美久久黑人一区二区| 亚洲第一青青草原| 欧美大码av| av视频免费观看在线观看| 国产亚洲欧美精品永久| 亚洲精品久久成人aⅴ小说| 久久久久久亚洲精品国产蜜桃av| 99国产极品粉嫩在线观看| 欧美日韩福利视频一区二区| 12—13女人毛片做爰片一| 久久久久国内视频| 亚洲一码二码三码区别大吗| 欧美不卡视频在线免费观看 | 精品国产一区二区久久| 亚洲精品在线观看二区| e午夜精品久久久久久久| 好男人电影高清在线观看| 99香蕉大伊视频| 丰满饥渴人妻一区二区三| 日韩人妻精品一区2区三区| 成年人午夜在线观看视频| 日本a在线网址| 国产免费男女视频| 最近最新中文字幕大全免费视频| 一边摸一边抽搐一进一出视频| 欧美激情高清一区二区三区| 国产精品久久久久成人av| 国产精品二区激情视频| 精品一区二区三区av网在线观看| 国产欧美日韩精品亚洲av| 亚洲情色 制服丝袜| 午夜影院日韩av| 国产精品秋霞免费鲁丝片| 12—13女人毛片做爰片一| tocl精华| 每晚都被弄得嗷嗷叫到高潮| 99精品久久久久人妻精品| 女人被躁到高潮嗷嗷叫费观| 久99久视频精品免费| 欧美 日韩 精品 国产| 久久精品亚洲av国产电影网| 丝袜在线中文字幕| 视频区图区小说| 久久人人97超碰香蕉20202| 免费观看a级毛片全部| 久久婷婷成人综合色麻豆| 黄片大片在线免费观看| 丝袜在线中文字幕| 国产国语露脸激情在线看| 91字幕亚洲| 三级毛片av免费| 黑人巨大精品欧美一区二区蜜桃| 国产亚洲欧美98| 亚洲成人免费av在线播放| 欧美丝袜亚洲另类 | 悠悠久久av| 国产av又大| 精品一区二区三区视频在线观看免费 | 国内毛片毛片毛片毛片毛片| 一本大道久久a久久精品| 欧美日韩亚洲综合一区二区三区_| 精品卡一卡二卡四卡免费| 国产亚洲欧美在线一区二区| 午夜日韩欧美国产| 国产伦人伦偷精品视频| 亚洲国产精品合色在线| 久久精品成人免费网站| 亚洲精品国产区一区二| 久久久国产成人免费| 丝袜在线中文字幕| 18禁黄网站禁片午夜丰满| 99久久人妻综合| e午夜精品久久久久久久| 午夜福利视频在线观看免费| 校园春色视频在线观看| 大香蕉久久成人网| 老熟妇乱子伦视频在线观看| av超薄肉色丝袜交足视频| 国产精品98久久久久久宅男小说| 大香蕉久久成人网| 18禁观看日本| 女人被狂操c到高潮| 久久国产精品男人的天堂亚洲| 大码成人一级视频| 国产免费av片在线观看野外av| 中文字幕最新亚洲高清| 天堂√8在线中文| 久久国产亚洲av麻豆专区| 亚洲av电影在线进入| 亚洲一区二区三区不卡视频| 在线看a的网站| 又黄又爽又免费观看的视频| 最新的欧美精品一区二区| 热re99久久精品国产66热6| 欧美成狂野欧美在线观看| 每晚都被弄得嗷嗷叫到高潮| 色婷婷av一区二区三区视频| 久久狼人影院| av一本久久久久| a级片在线免费高清观看视频| 人人妻人人添人人爽欧美一区卜| 国产欧美亚洲国产| 咕卡用的链子| 夜夜爽天天搞| 欧美日韩黄片免| 一进一出抽搐动态| 久久久国产欧美日韩av| 大码成人一级视频| 欧美午夜高清在线| 日日爽夜夜爽网站| 久久人妻熟女aⅴ| 欧美日韩精品网址| 欧美日韩视频精品一区| 亚洲少妇的诱惑av| 老司机福利观看| 精品人妻在线不人妻| 国产在线精品亚洲第一网站| 90打野战视频偷拍视频| 亚洲国产精品合色在线| 大香蕉久久成人网| 久久精品国产a三级三级三级| 香蕉国产在线看| 国产高清视频在线播放一区| 日韩一卡2卡3卡4卡2021年| 精品亚洲成国产av| 国内久久婷婷六月综合欲色啪| 欧美黑人欧美精品刺激| 99久久99久久久精品蜜桃| 久久精品国产亚洲av高清一级| 黄色成人免费大全| 在线观看日韩欧美| 国产在线观看jvid| 国产片内射在线| 亚洲熟妇熟女久久| 在线观看舔阴道视频| 男女下面插进去视频免费观看| 国产亚洲av高清不卡| 精品熟女少妇八av免费久了| 91成人精品电影| √禁漫天堂资源中文www| 18禁国产床啪视频网站| 中文字幕高清在线视频| 亚洲av第一区精品v没综合| 久久国产精品影院| 欧美人与性动交α欧美精品济南到| 婷婷精品国产亚洲av在线 | 欧美黑人欧美精品刺激| 91在线观看av| 色精品久久人妻99蜜桃| ponron亚洲| 久久精品国产99精品国产亚洲性色 | 在线播放国产精品三级| 国产精品九九99| 两性午夜刺激爽爽歪歪视频在线观看 | 欧美色视频一区免费| 亚洲伊人色综图| 成熟少妇高潮喷水视频| 高清视频免费观看一区二区| 国产亚洲精品一区二区www | 99久久人妻综合| 欧美日韩亚洲国产一区二区在线观看 | 国产成人影院久久av| 天天躁夜夜躁狠狠躁躁| 高清av免费在线| 99精品欧美一区二区三区四区| 黑人操中国人逼视频| 国产精品久久久久久人妻精品电影| 黄网站色视频无遮挡免费观看| 热99re8久久精品国产| 天天躁狠狠躁夜夜躁狠狠躁| 国产日韩一区二区三区精品不卡| 精品久久久久久久毛片微露脸| 曰老女人黄片| 亚洲熟女精品中文字幕| 国产一区有黄有色的免费视频| 欧美大码av| 精品久久久久久久毛片微露脸| 午夜激情av网站| 亚洲色图 男人天堂 中文字幕| 亚洲成人免费av在线播放| 国产人伦9x9x在线观看| 91成年电影在线观看| 建设人人有责人人尽责人人享有的| 99精品久久久久人妻精品| 国产一区有黄有色的免费视频| 最新美女视频免费是黄的| 一个人免费在线观看的高清视频| 欧美不卡视频在线免费观看 | 中文字幕av电影在线播放| а√天堂www在线а√下载 | 午夜老司机福利片| 日韩欧美一区视频在线观看| 99精品在免费线老司机午夜| 久久99一区二区三区| 天堂俺去俺来也www色官网| 亚洲人成伊人成综合网2020| 国产精品久久久人人做人人爽| 91精品国产国语对白视频| 免费人成视频x8x8入口观看| 午夜精品国产一区二区电影| 国产男女内射视频| 国产精品九九99| 手机成人av网站| 中文字幕人妻熟女乱码| 啦啦啦免费观看视频1| 精品国产乱子伦一区二区三区| 黄频高清免费视频| 老汉色∧v一级毛片| 大码成人一级视频| 亚洲精品美女久久av网站| 女性被躁到高潮视频| 波多野结衣一区麻豆| 国产精品久久久av美女十八| 黄色a级毛片大全视频| 性色av乱码一区二区三区2| 国产99白浆流出| 老司机午夜福利在线观看视频| 男女床上黄色一级片免费看| 久久精品亚洲熟妇少妇任你| 欧美人与性动交α欧美精品济南到| 美女午夜性视频免费| 99精品在免费线老司机午夜| 亚洲av欧美aⅴ国产| 亚洲 国产 在线| 色尼玛亚洲综合影院| 午夜福利在线免费观看网站| 欧美成人午夜精品| 在线av久久热| 亚洲精品美女久久久久99蜜臀| 又黄又粗又硬又大视频| 欧美色视频一区免费| 精品一区二区三卡| 看片在线看免费视频| 色精品久久人妻99蜜桃| 精品人妻1区二区| 久久人妻av系列| 精品国产乱码久久久久久男人| 亚洲精品在线美女| 可以免费在线观看a视频的电影网站| 99国产精品免费福利视频| 高清欧美精品videossex| 桃红色精品国产亚洲av| 国产野战对白在线观看| 黄色怎么调成土黄色| 69av精品久久久久久| 亚洲 国产 在线| 国产人伦9x9x在线观看| 少妇裸体淫交视频免费看高清 | 亚洲成人免费av在线播放| 亚洲在线自拍视频| 亚洲aⅴ乱码一区二区在线播放 | 老司机影院毛片| 激情视频va一区二区三区| 首页视频小说图片口味搜索| 一二三四社区在线视频社区8| 亚洲午夜精品一区,二区,三区| 无人区码免费观看不卡| 18禁国产床啪视频网站| av不卡在线播放| 一区在线观看完整版| 久久精品aⅴ一区二区三区四区| netflix在线观看网站| 亚洲人成电影观看| 久久国产精品男人的天堂亚洲| 日日摸夜夜添夜夜添小说| 午夜福利在线观看吧| 久久久国产欧美日韩av| 日日夜夜操网爽| 视频在线观看一区二区三区| 波多野结衣av一区二区av| av国产精品久久久久影院| 精品福利永久在线观看| 午夜免费成人在线视频| 我的亚洲天堂| 国产精品 国内视频| 91成人精品电影| 99精品在免费线老司机午夜| 亚洲精品自拍成人| 在线视频色国产色| av国产精品久久久久影院| 1024视频免费在线观看| 欧美av亚洲av综合av国产av| 久久 成人 亚洲| 久久精品成人免费网站| 丁香欧美五月| 色综合婷婷激情| 国产色视频综合| 国产主播在线观看一区二区| 亚洲一区二区三区不卡视频| 欧美日韩一级在线毛片| 亚洲成人国产一区在线观看| 色老头精品视频在线观看|