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

    近20 a蒙古國土壤風蝕變化特征及主要影響因素分析

    2022-09-24 07:15:08雷加強范敬龍王海峰呂振濤
    干旱區(qū)研究 2022年4期
    關(guān)鍵詞:風蝕蒙古國模數(shù)

    郭 茵, 雷加強,3, 范敬龍,3, 王海峰,4, 呂振濤

    (1.中國科學院新疆生態(tài)與地理研究所國家荒漠-綠洲生態(tài)建設(shè)工程技術(shù)研究中心,新疆 烏魯木齊 830011;2.中國科學院大學,北京 100049;3.中國科學院新疆生態(tài)與地理研究所塔克拉瑪干沙漠研究站,新疆 庫爾勒 841000;4.中國科學院新疆生態(tài)與地理研究所莫索灣沙漠研究站,新疆 石河子 832000)

    土壤風蝕是土地退化最主要的原因之一。世界上約三分之一的土地受到風蝕的影響,其中嚴重風蝕土地占比達50%以上[1-3]。土壤風蝕會對陸地生態(tài)系統(tǒng)以及人類社會發(fā)展產(chǎn)生諸多不利影響[4],如顯著降低土壤肥力和植物生產(chǎn)力[5],導致空氣質(zhì)量降低等[6]。此外,地表風蝕形成的細粒物質(zhì)是大氣氣溶膠的主要來源[7],為沙塵暴提供了物源。只要有發(fā)生風蝕的條件,就可能產(chǎn)生較大規(guī)模的風沙災害[8],如2021年3月14日蒙古國發(fā)生的特大沙塵暴事件就危及到東亞廣大地區(qū)。

    土壤風蝕的定量評估方法包括野外風蝕直接測量[9]、風洞模擬實驗[10]、基于同位素示蹤法的風蝕速率計算等[11]。利用風蝕監(jiān)測設(shè)備或風洞模擬實驗具有一定的局限性[12],難以對較大空間范圍進行風蝕評估。隨著遙感技術(shù)的發(fā)展,越來越多的學者利用地理信息系統(tǒng)進行風蝕建模,發(fā)展出WEQ、RWEQ、WEPS、SWEEP 等眾多模型,其中以RWEQ模型應用最為廣泛。RWEQ 模型以WEQ 模型為基礎(chǔ),融合了自然過程和實際經(jīng)驗,能夠?qū)⒁巴鈹?shù)據(jù)與計算機模型相結(jié)合,來描述物理的風蝕過程[13],在實際中得到了廣泛的應用。RWEQ 模型起初是美國農(nóng)業(yè)部(USDA)為估算耕地土壤流失而開發(fā)的一種基于經(jīng)驗的模型[14],之后經(jīng)過數(shù)次改進和驗證,在其他各地得到了廣泛的應用。Zhang 等[15]使用RWEQ模型,對中國內(nèi)蒙古地區(qū)1990—2015年的土壤風蝕進行了模擬,評估了氣候變化對風蝕的影響;Lin 等[12]利用RWEQ 模型估算了1982—2015 年中國河西地區(qū)的潛在風蝕率,Chi 等[16]利用RWEQ模型評估了中國土壤風蝕的時空格局,Li等[17]利用RWEQ 模型模擬了1986—2099 年中亞土壤風蝕的時空動態(tài)。蒙古國是東亞地區(qū)沙塵暴事件發(fā)生頻率最高的地區(qū)[18],是東亞地區(qū)沙塵暴災害的主要起源地[19]。前人使用基于專家的方法評估了蒙古國的風蝕速率和程度,結(jié)果表明蒙古國大約30%的領(lǐng)土受到風蝕的影響,其中戈壁沙漠被認為是受風蝕影響最嚴重的區(qū)域[20]。Natsagdorj 等[21]利用氣象觀測和激光雷達,發(fā)現(xiàn)蒙古國各地的風蝕總體上有所增加。師華定等[22]利用地理信息系統(tǒng)和遙感技術(shù),通過模糊c 均值聚類(FCM)方法對植被覆蓋率、地形起伏度、土壤干燥度和風場強度4 個主要蒙古高原風蝕環(huán)境因子數(shù)據(jù)進行模糊聚類,結(jié)合專家知識,將環(huán)境組合與特定的風蝕危險度相對應,得到風蝕危險度圖,以此對蒙古高原土壤風蝕危險度進行評價研究,并分析了蒙古高原土壤風蝕危險度的空間格局特征。FCM方法可有效地獲取風蝕-環(huán)境關(guān)系知識,為預測性風蝕危險度制圖提供依據(jù),此研究并沒有對蒙古高原土壤風蝕危險度進行定量評價。Mandakh 等[23]在ArcGIS 環(huán)境中使用風蝕方程模型(WEQ)評估了蒙古國各地的風蝕情況,盡管此研究對蒙古國的風蝕進行了定量分析,但是缺乏對蒙古國長時間序列的風蝕趨勢的分析。目前,關(guān)于蒙古國的風蝕研究較為有限,亟待開展蒙古國長期風蝕狀況定量描述和評估研究。本文使用RWEQ 模型,計算蒙古國全域的風蝕模數(shù),揭示風蝕的時空變化規(guī)律和趨勢,分析識別主要影響因素,有助于深入了解蒙古國土壤風蝕整體特征,為重點地區(qū)風沙災害和荒漠化防治工作提供重要參考。

    1 研究區(qū)概況

    蒙古國位于亞洲東中部(41°32′~52°15′N,87°44′~119°56′E),面積156.65×104km2[24],是世界第二大內(nèi)陸國家。蒙古國平均海拔約1580 m,最高點(西部的khuiten 峰)海拔4653 m,最低點(東部的Khookh Nuur 湖洼地)海拔518 m,首都烏蘭巴托市海拔約1310 m。蒙古國遠離海洋,具有顯著的大陸性氣候特征,四季分明,年均溫和日均溫變化較大;年平均降水量約120~250 mm,降水量從北到南、從東到西遞減,降水多集中在7—8月[25]。根據(jù)蒙古國家氣象和環(huán)境監(jiān)測局(https://eic.mn/)提供的蒙古國自然地帶分布數(shù)據(jù)顯示,南部為沙漠和荒漠草原地區(qū),東南部是戈壁沙漠區(qū);典型草原主要分布在蒙古國中部東西延伸;森林集中分布在蒙古國北部山區(qū)(圖1)。其中,荒漠草原、半荒漠和沙漠地區(qū)約占41.3%,這些地區(qū)大風頻發(fā),土壤風蝕發(fā)生強烈,生態(tài)環(huán)境脆弱,年平均風速為4~6 m·s-1[23]。蒙古國南部的荒漠地區(qū)由于土壤疏松,更易被風蝕并形成沙塵暴[23,26],年沙塵暴日數(shù)可達30~120 d;蒙古國是一個沙漠化嚴重的國家,沙漠化的加劇直接增加了風蝕和沙塵排放源,威脅到東亞其他國家的空氣質(zhì)量[27]。

    圖1 研究區(qū)概況示意圖Fig.1 Geographical distribution map of Mongolia

    2 研究方法與數(shù)據(jù)來源

    2.1 研究方法

    2.1.1 風蝕模型 RWEQ模型充分考慮了包括氣候條件、土壤特性、表面粗糙度、植被覆蓋和其他環(huán)境要素,主要基于式(1)~(3)來估算風蝕量。

    式中:SWEQ表示風蝕模數(shù)(t·hm-2·a-1);Qmax表示最大輸沙能力(kg·m-1);s表示關(guān)鍵地塊長度(m);x是模擬區(qū)域不可侵蝕邊界到下風口的距離(m);WF表示氣候因子(kg·m-1);EF表示土壤可蝕性因子(無量綱);SCF表示土壤結(jié)皮因子(無量綱);K′表示土壤粗糙度影響因子(無量綱);COG表示植被影響因子。具體計算過程如下式所示:

    式中:Wf表示風影響因子(kg·m-1·s-1);ρ表示空氣密度(kg·m-3);g表示重力加速度(m·s-2);SW表示土壤濕度因子;SD表示積雪覆蓋影響因子;W表示風力值(m3·s-3);U2指2 m 高的風速(m·s-1);Ut指風速閾值,一般取5 m·s-1;Nd是風速測量的時間間隔(d);N代表風速觀測頻率(d);EL表示從數(shù)字高程模型(DEM)數(shù)據(jù)獲得的高程(km);T表示絕對溫度(K);SD由積雪覆蓋度計算獲得;P表示測定時段內(nèi)積雪覆蓋深度大于25.4 mm的發(fā)生概率。

    土壤可蝕性因子EF計算過程如下式所示:

    式中:Sa為土壤砂粒含量(%);Si為土壤粉粒含量(%);Cl是黏土含量(%);Sa/Cl為土壤砂粒與黏土之比;OM為土壤有機質(zhì)含量(%);CaCO3為碳酸鈣含量(%)。

    土壤結(jié)皮影響因子SCF計算過程如下式所示:

    土壤質(zhì)地和有機質(zhì)含量的年際變化非常緩慢,可將土壤可蝕性因子和結(jié)皮因子視為常量值,假設(shè)土壤可蝕性因子和結(jié)皮因子不隨時間改變。

    地表粗糙度影響因子K′,根據(jù)Xu 等[28]的研究結(jié)果,計算如下式所示:

    式中:α表示地形坡度,可從DEM數(shù)據(jù)中提取。

    植被影響因子COG表示一定植被覆蓋條件下對土壤風蝕的抑制程度,采用最大合成法(MVC)得到半月NDVI 數(shù)據(jù),并用像元二分法求得2001—2020年的半月植被覆蓋度,計算過程如下式所示:

    式中:FVC 為植被覆蓋率(%);NDVIsoil為裸土像素的NDVI值;NDVIveg為植被的NDVI值。

    2.1.2 變化趨勢分析方法 利用Sen’s趨勢分析法,計算了2001—2020 年蒙古國土壤風蝕以及氣象因子(風速、氣溫、降水)的變化趨勢。該方法是非參數(shù)檢驗,只需要保證樣本的獨立性,不需要擬合正態(tài)分布[29],方法如下式所示:式中:xi和xj是第i年和第j年的風蝕模數(shù)值或風速、氣溫、降水的數(shù)值;n是研究期的總年數(shù)。Slope表示時間序列數(shù)據(jù)的變化率,正值表示風蝕或氣象因子(風速、氣溫、降水)隨時間變化總體呈增加趨勢,負值表示總體呈減少趨勢。

    2.2 數(shù)據(jù)來源

    (1)氣象數(shù)據(jù)

    氣象數(shù)據(jù)包括風速、氣溫、降水、積雪深度、土壤水分,均來自全球陸地數(shù)據(jù)同化系統(tǒng)GLDAS2.1產(chǎn)品[30]。GLDAS數(shù)據(jù)是美國航空航天局(NASA)戈達德空間飛行中心(GSFC)和美國海洋和大氣局(NOAA)國家環(huán)境預報中心(NCEP)聯(lián)合發(fā)布的基于衛(wèi)星、陸面模式和地面觀測數(shù)據(jù)的同化產(chǎn)品,它能夠提供多種驅(qū)動數(shù)據(jù),這些數(shù)據(jù)來源于大氣同化產(chǎn)品、再分析和實際觀測。GLDAS2.1產(chǎn)品時間分辨率為3 h,空間分辨率為0.25°,與地面氣象觀測站實測數(shù)據(jù)具有較高的一致性[30]。

    (2)土壤數(shù)據(jù)

    土壤數(shù)據(jù)包括土壤有機質(zhì)含量、黏粒含量、粉粒含量、沙粒含量、CaCO3含量,其中土壤有機質(zhì)含量、黏粒含量、粉粒含量、沙粒含量來源于Open Land Map(OLM)數(shù)據(jù)集,CaCO3含量來源于世界土壤數(shù)據(jù)庫(HWSD)。

    (3)植被數(shù)據(jù)

    選用由美國航空航天局(NASA)發(fā)布的NDVI中分辨率成像光譜儀獲得的陸地2級標準數(shù)據(jù)產(chǎn)品MOD13Q1,時間分辨率為16 d,空間分辨率為250 m,可用于計算植被覆蓋度。

    (4)地形數(shù)據(jù)

    選擇由美國航空航天局(NASA)和美國測繪局(NIMA)聯(lián)合測量的數(shù)字高程模型(SRTM)版本4,空間分辨率為90 m。這是目前質(zhì)量最高的SRTM數(shù)據(jù)集,水平標準誤差為1 m,垂直標準誤差為0 m,可用于計算坡度和地表粗糙度。

    (5)其他數(shù)據(jù)

    牲畜數(shù)量數(shù)據(jù)來自于蒙古國國家統(tǒng)計局(www.1212.mn)。土地利用數(shù)據(jù)使用的是由中國國家基礎(chǔ)地理信息中心(NGCC)提供的GlobeLand30產(chǎn)品[31],數(shù)據(jù)集使用了包括20000 多張美國陸地資源衛(wèi)星(Landsat)和中國環(huán)境減災衛(wèi)星(HJ-1)的多光譜影像,具有較高的空間分辨率和數(shù)據(jù)精度[32]。

    將以上所有柵格數(shù)據(jù)的空間分辨率重新采樣至500 m,時間分辨率為16 d,全年有23 個時間段。對計算結(jié)果進行匯總,得到每年的風蝕總量。

    3 結(jié)果與分析

    3.1 風蝕模數(shù)的空間變化特征

    從蒙古國土壤風蝕模數(shù)空間分布來看(圖2),土壤風蝕變化空間差異較大,風蝕模數(shù)的最大值與最小值的差值在0~60 t·hm-2·a-1之間。整體上,土壤風蝕模數(shù)呈現(xiàn)北部低、南部高的空間分布模式,由于南部的戈壁沙漠區(qū)土壤質(zhì)地較粗,物源豐富,且可蝕性較高,加之降水少,地表植被稀疏,因此該區(qū)成為蒙古國土壤風蝕最為嚴重的地區(qū)。

    圖2 2001—2020年蒙古國土壤風蝕模數(shù)空間分布特征Fig.2 Spatial distribution characteristics map of wind erosion in Mongolia from 2001 to 2020

    由圖3、圖4可知,對風蝕模數(shù)進行多年平均,從年內(nèi)空間分布可以看出,風蝕整體上春季(3—5月)最強,約占全年風蝕模數(shù)的45%;其次為秋季(9—11月),約占全年風蝕模數(shù)的33%;再次為夏季(6—8月),約占全年風蝕模數(shù)的12%;冬季(12月—次年2月)最弱,約占全年風蝕模數(shù)的10%。蒙古高原春季降水稀少,土壤干燥,植物尚未進入快速生長期,植被覆蓋度低,因此春季風蝕非常強烈。冬季土壤表層凍結(jié),且多有積雪覆蓋,因此風蝕最弱。夏季和秋季植被條件均較好,但夏季降水集中,土壤濕度明顯大于秋季,因此夏季風蝕強度小于秋季。氣象觀測表明,蒙古國的沙塵事件主要發(fā)生在春季[18],與本文風蝕主要發(fā)生在春季的結(jié)論一致。

    圖3 2001—2020年蒙古國季節(jié)平均土壤風蝕模數(shù)空間分布Fig.3 Spatial distribution map of seasonal average soil wind erosion in Mongolia from 2001 to 2020

    圖4 2001—2020年蒙古國季節(jié)平均土壤風蝕模數(shù)變化Fig.4 Variation of seasonal average soil wind erosion modulus in Mongolia from 2001 to 2020

    3.2 風蝕模數(shù)的時間變化特征

    由圖5 可知,在2001—2020 年間,蒙古國土壤風蝕模數(shù)總體呈顯著上升趨勢(R2=0.2258,P<0.05),其中2018年最強(4.4463 t·hm-2·a-1),2008年最弱(1.369 t·hm-2·a-1)。從時間序列的變化特征來看,2005—2009 年的風蝕區(qū)域明顯小于其他年份,而2009 年后風蝕范圍明顯增大。根據(jù)風蝕模數(shù)的年際變化,2001—2020 年期間風蝕變化分為2 個明顯的階段,2001—2008 年為風蝕模數(shù)下降時期(R2=0.536,P<0.05),2010—2020 年為風蝕模數(shù)增加時期。

    從土壤風蝕變化速率空間差異來看(圖6),蒙古國南部戈壁沙漠地區(qū)的風蝕變化速率約為2 t·hm-2·a-1,明顯高于其他區(qū)域;蒙古國北部風蝕最低,變化最??;蒙古國西南部風蝕速率變化為負,表明此區(qū)域風蝕呈下降趨勢,平均每年的風蝕變化速率約-0.01~-0.6 t·hm-2·a-1。總體來看,蒙古國土壤風蝕模數(shù)高值區(qū)變化速率快,而低值區(qū)變化速率慢。蒙古國北部以森林生態(tài)系統(tǒng)為主,植被類型以高山泰加林、高山森林和森林草原為主,植被覆蓋度較高,土壤風蝕閾值不大,整體森林生態(tài)系統(tǒng)穩(wěn)定,受氣候變化及人為活動影響小,風蝕很小,且變化率很?。幻晒艊喜炕哪鷳B(tài)系統(tǒng)脆弱,植被稀少,受氣候變化影響大,因此風蝕年際變化率也較大;蒙古國西南部風蝕呈下降趨勢的主要原因是此區(qū)域風速呈下降趨勢,風蝕量常年較低。

    圖6 2001—2020年蒙古國土壤風蝕模數(shù)變化趨勢的空間分布Fig.6 Spatial distribution of variation trend of soil wind erosion in Mongolia from 2001 to 2020

    3.3 土壤風蝕的影響因子分析

    3.3.1 氣候因子變化及影響 影響土壤風蝕的主要氣候因子有風速、氣溫、降水等。本文采用Sen’s趨勢分析法分析了蒙古國主要氣候因子變化趨勢及空間分布。從圖7 可以看出,近20 a 來蒙古國年最大風速呈明顯波動變化,2008年最低(14.04 m·s-1),2001 年最高(16.26 m·s-1);2001—2008 年,年最大風速呈下降趨勢(R2=0.427,P>0.05);2009—2020年,年最大風速呈顯著增加趨勢(R2=0.3959,P<0.05)。土壤風蝕模數(shù)隨著風速變化趨勢的波動而波動,二者變化規(guī)律基本呈現(xiàn)對應關(guān)系(圖5,圖7),表明風速是影響土壤風蝕最重要的自然動力條件。

    圖5 2001—2020年蒙古國多年平均土壤風蝕模數(shù)變化趨勢Fig.5 Variation trend graph of average soil wind erosion in Mongolia from 2001 to 2020

    圖7 2001—2020年蒙古國年最大風速變化趨勢Fig.7 Variation trend of annual average maximum wind speed in Mongolia from 2001 to 2020

    結(jié)合年均風速變化趨勢空間分布(圖8)來看,蒙古國南部戈壁沙漠區(qū)風速顯著上升,其風蝕量呈增加趨勢;中部風速下降,風蝕量呈減少趨勢。此外,將蒙古國2001—2020年逐年平均風速與逐年風蝕模數(shù)進行相關(guān)性分析(圖9)發(fā)現(xiàn),在蒙古國僅3.3%左右的區(qū)域平均風速與風蝕模數(shù)呈負相關(guān)關(guān)系,絕大數(shù)區(qū)域(96.6%)風蝕模數(shù)都與平均風速呈正相關(guān)。尤其是在蒙古國南部及西南部的戈壁區(qū)域,風速與風蝕量呈顯著正相關(guān)關(guān)系,說明風速是影響土壤風蝕最重要的自然動力條件。這也解釋了蒙古國西南部風蝕量變化速率下降是受風速下降的影響。

    圖8 2001—2020年蒙古國年均風速變化趨勢空間分布Fig.8 Spatial distribution of annual average wind speed change trend in Mongolia from 2001 to 2020

    圖9 蒙古國年均風速與風蝕模數(shù)相關(guān)關(guān)系Fig.9 Correlation between annual average wind speed and wind erosion modulus in Mongolia

    氣溫和降水會通過影響土壤濕度和植物生長而間接影響風蝕[33]。從圖10 可以看出,近20 a 來,蒙古國大部分地區(qū)年均溫呈上升趨勢,僅有北部少部分區(qū)域年均溫呈下降趨勢。南部增溫最高,增溫速率達到0.06 ℃·a-1,升溫促進了地表水的蒸發(fā),導致土壤濕度降低,可蝕性增加,特別是在缺乏植被的地區(qū),升溫增加了風蝕的發(fā)生。降水與表層土壤含水量密切相關(guān),降水量減少增加了風蝕的發(fā)生。整體來看,近20 a來蒙古國降水呈顯著增加趨勢(圖11b),但降水變化存在較大空間差異(圖11a),北部地區(qū)降水明顯增加,中部和南部部分地區(qū)降水減少,尤其南部戈壁地區(qū)呈明顯減少趨勢(圖11c)。南部地區(qū)干旱災害發(fā)生的可能性增加,加上風力的作用,風蝕量也呈增大趨勢。

    圖10 2001—2020年蒙古國年均溫變化趨勢時空分布Fig.10 Temporal and spatial distribution of annual mean temperature change trend in Mongolia from 2001 to 2020

    圖11 2001—2020年蒙古國年降水量變化趨勢時空分布Fig.11 Temporal and spatial distribution of annual rainfall variation trend in Mongolia from 2001 to 2020

    3.3.2 植被覆蓋度變化及影響 植被覆蓋可增加地

    表粗糙度和削弱近地表風速,從而降低地表風蝕危險性[15,34]。圖12 為近20 a 來蒙古國FVC 的變化,可以看出,蒙古國大部分區(qū)域FVC 整體呈增長趨勢,東部部分地區(qū)植被覆蓋度快速增加,因此土壤風蝕降低;但部分地區(qū)尤其是南部和西部的荒漠草原區(qū)多地植被覆蓋度減少,導致風蝕發(fā)生風險增大。

    圖12 2001—2020年蒙古國植被覆蓋度變化趨勢Fig.12 Change trend of annual vegetative cover in Mongolia from 2001 to 2020

    4 討論

    4.1 風蝕模數(shù)的結(jié)果驗證

    RWEQ 模型起初是美國農(nóng)業(yè)部(USDA)為估算耕地土壤流失而開發(fā)的一種基于經(jīng)驗的模型[14],之后經(jīng)過數(shù)次改進和驗證,不僅是將其用于評估農(nóng)田發(fā)生的土壤風蝕,如Li等[35]將其用于評估內(nèi)蒙古地區(qū)的固沙服務(wù)功能,Chi 等[36]評估了中國全域的土壤風蝕,也有一些學者將其用于中亞地區(qū)的土壤風蝕評估[17,37]。然而其結(jié)果往往需要其他方法估算的結(jié)果來對比驗證。

    綜上所述,蒙古國2001—2020年多年平均土壤風蝕模數(shù)為2.736 t·hm-2·a-1。137Cs示蹤法是估算土壤侵蝕速率的常用方法[38]。前人基于137Cs 監(jiān)測方法計算了蒙古高原的風蝕模數(shù),蒙古國境內(nèi)的監(jiān)測樣點從北向南涵蓋了森林草原區(qū)、典型草原區(qū)、荒漠化草原區(qū)、草原化荒漠區(qū)、戈壁區(qū)、農(nóng)牧交錯區(qū)[39]。將本研究中基于RWEQ 模型模擬的20 a 的土壤風蝕模數(shù)與劉紀遠等[39]使用137Cs 實測的土壤風蝕模數(shù)基于相同的地理坐標點進行了相關(guān)性分析(圖13),結(jié)果表明,RWEQ 模型模擬結(jié)果與137Cs 測量值顯著相關(guān)(R=0.91,P<0.01)。Mandakh 等[23]使用WEQ模型評估的蒙古國風蝕模數(shù)在2.7~27.5 t·hm-2·a-1之間,比本文用RWEQ 估算結(jié)果略偏低,但數(shù)量值大體上相當。土壤風蝕的發(fā)生和發(fā)展依賴于侵蝕動力(氣候因子如風速)和可蝕性因子(地表狀況如植被蓋度、土壤特性等)之間的相互作用,土壤可蝕性因子可以衡量土壤表層被風力侵蝕的相對可能性。Jugder等[40]從土壤可蝕性的角度評價了蒙古國土壤被風力侵蝕的相對可能性,發(fā)現(xiàn)從蒙古國南部的沙漠和荒漠草原地區(qū)到西部的沙地,存在著廣泛的高可蝕性區(qū)域,與本文中的嚴重風蝕區(qū)域也大致相同。中等可蝕性區(qū)域主要覆蓋西部的干旱草原區(qū)和荒漠草原區(qū),低可蝕性區(qū)域從蒙古國西部延伸到北部。Mandakh 等[23]使用WEQ 模型評估了蒙古國風蝕模數(shù),分析了風蝕的空間分布規(guī)律,與本文結(jié)論一致。

    圖13137Cs測量和RWEQ模型模擬風蝕模數(shù)的相關(guān)系數(shù)Fig.13 Correlation coefficient of137Cs measurement and RWEQ model simulation wind erosion modulus

    4.2 人為活動因素變化及影響

    土壤風蝕的發(fā)生與氣候條件和下墊面條件有密切關(guān)系,人為活動通過影響下墊面而影響風蝕的發(fā)生。過度放牧以及煤礦開采等人類活動加劇了蒙古國土地退化[41]。由圖14 可知,在2001—2020年間,蒙古國的牲畜數(shù)量由2001年的26075千頭增加到2020年的67068千頭,每年大約以2400千頭的速度顯著增加。牲畜數(shù)量的增加使草原壓力增大,而寬松的放牧政策造成多數(shù)草場過度放牧,牲畜踩踏和過牧行為破壞了土壤表層結(jié)構(gòu),降低了冬、春季地表的植被蓋度,增加了土壤的風蝕風險。因此,過度放牧可能是導致蒙古國土壤風蝕增加的一個人為活動因素。

    圖14 2001—2020年蒙古國牲畜數(shù)量變化趨勢Fig.14 Variation trend of livestock quantity in Mongolia from 2001 to 2020

    近20 a來,蒙古國經(jīng)濟社會快速發(fā)展,為了適應生存和發(fā)展的需求,土地利用類型也發(fā)生了較大變化。從土地利用類型轉(zhuǎn)化(圖15)來看,蒙古國境內(nèi)有1733 km2的草地轉(zhuǎn)化為耕地,部分未利用地和林地也有轉(zhuǎn)化為耕地的情況。與傳統(tǒng)畜牧業(yè)生產(chǎn)方式相比,發(fā)展種植業(yè),會導致更嚴重的土壤風蝕發(fā)生。其他土地類型轉(zhuǎn)化為耕地后,人為翻動表土,土壤結(jié)構(gòu)、水分和養(yǎng)分狀況發(fā)生了很大變化,導致局部地區(qū)風蝕加劇。

    圖15 2001—2020年蒙古國主要土地利用類型轉(zhuǎn)化Fig.15 Transformation of main land use types in Mongolia from 2001 to 2020

    4.3 影響因素耦合的空間異質(zhì)性作用

    氣溫和降水通過影響土壤濕度和植被覆蓋度而間接地影響風蝕。2001—2020 年蒙古國全域氣溫和降水均呈增加趨勢,植被覆蓋度也呈增加趨勢,對地表風蝕應發(fā)揮一定的抑制作用,但全域風蝕仍整體呈增加趨勢,似乎存在一定的矛盾。但分析發(fā)現(xiàn),原因在于影響風蝕因素耦合存在較大的時空異質(zhì)性。對于一些風蝕相對嚴重的地區(qū),植被覆蓋度的增加對減少風蝕沒有明顯影響,這是因為高風蝕區(qū)域的植被覆蓋度較低,且增加不顯著(圖12),風速等其他因素的變化對風蝕有更顯著的影響,而且在正值風季的春季作用較為強烈。學者們對蒙古國土壤風蝕的研究中發(fā)現(xiàn),不同的氣候因素影響不同的區(qū)域,干旱對戈壁和沙漠地區(qū)影響更大[23]。因此,從蒙古國全域的尺度來看,風速是土壤風蝕發(fā)生的主導因素,風速越大,風蝕越強。而且強風蝕的干旱荒漠土地面積占比極高,因此蒙古國整體風蝕仍呈增加趨勢。

    5 結(jié)論

    本文通過RWEQ 模型對蒙古國2001—2020 年的土壤風蝕時空變化狀況及對主要影響因子進行了分析,得到的結(jié)論如下:

    (1)蒙古國土壤風蝕年內(nèi)季節(jié)變化明顯,全年45%的風蝕發(fā)生在春季,秋季其次,冬季最少。

    (2)從2001—2020 年,蒙古國土壤風蝕呈年際波動變化,總體上呈顯著增加趨勢,風蝕模數(shù)以0.06 t·hm-2·a-1的速度增長。

    (3)土壤風蝕具有明顯的空間異質(zhì)性,整體呈南部高、北部低的空間分布格局,風蝕嚴重區(qū)域位于蒙古國南部與中國接壤的戈壁和沙漠地區(qū);風蝕變化速率空間變異較大,土壤風蝕高值區(qū)大于低值區(qū)。

    (4)蒙古國土壤風蝕變化整體上與氣候因子及引起的植被蓋度變化有關(guān),也與過牧、耕地開墾等人為活動增強有密切關(guān)系。為了有效防治沙漠化危害,蒙古國應在土地資源開發(fā)利用方面加強科學管理,對于嚴重沙漠化地區(qū)應積極采取可行的治理措施。

    猜你喜歡
    風蝕蒙古國模數(shù)
    黃土高原水蝕風蝕交錯帶迎風坡水蝕影響的風蝕特征
    基于單片機和模數(shù)化設(shè)計的低壓側(cè)電壓監(jiān)視與保護裝置
    能源工程(2021年2期)2021-07-21 08:40:02
    蒙古國某地綠松石玉的礦物組成及其成因
    模數(shù)化設(shè)計方法在景觀鋪裝設(shè)計中的應用
    綠色科技(2020年11期)2020-08-01 02:23:58
    蒙古國草原上的“青啤綠”
    商周刊(2018年12期)2018-07-11 01:27:24
    淺談從蒙古國語言中選用一些詞語
    基于LID模式的城區(qū)排澇模數(shù)探析
    蒙古國一處摩崖石刻被確認為《封燕然山銘》
    丹青少年(2017年4期)2017-02-06 03:08:48
    一種新型的RSA密碼體制模數(shù)分解算法
    令人震撼的侵蝕地貌(一)
    亚洲精品国产av成人精品| 国产av码专区亚洲av| 波多野结衣一区麻豆| 日本猛色少妇xxxxx猛交久久| 亚洲综合精品二区| 国产亚洲av片在线观看秒播厂| 精品人妻熟女毛片av久久网站| 亚洲精品乱码久久久久久按摩| 成人国产av品久久久| 男人操女人黄网站| 人人妻人人添人人爽欧美一区卜| 免费看av在线观看网站| 亚洲精品一二三| 交换朋友夫妻互换小说| 亚洲欧洲日产国产| 丝袜脚勾引网站| 少妇精品久久久久久久| 美女脱内裤让男人舔精品视频| 久久99一区二区三区| 少妇 在线观看| av天堂久久9| 精品午夜福利在线看| 国产一区亚洲一区在线观看| 国产精品三级大全| 欧美少妇被猛烈插入视频| 97在线视频观看| 青春草亚洲视频在线观看| 两个人看的免费小视频| 亚洲精品久久午夜乱码| 国产av国产精品国产| 丰满少妇做爰视频| 99久久精品国产国产毛片| 国产精品一区www在线观看| 性色av一级| 女的被弄到高潮叫床怎么办| 一级爰片在线观看| 成人午夜精彩视频在线观看| 汤姆久久久久久久影院中文字幕| 久久热在线av| 国产精品国产三级国产专区5o| 成人午夜精彩视频在线观看| 91午夜精品亚洲一区二区三区| 又粗又硬又长又爽又黄的视频| 国内精品宾馆在线| 国产免费又黄又爽又色| 日韩av不卡免费在线播放| 日本猛色少妇xxxxx猛交久久| 国产片特级美女逼逼视频| 国产高清三级在线| 国产精品 国内视频| 国产精品三级大全| 欧美3d第一页| 亚洲人成77777在线视频| 日本wwww免费看| 黄色毛片三级朝国网站| 在线观看人妻少妇| 天天操日日干夜夜撸| 精品午夜福利在线看| www日本在线高清视频| 国产又爽黄色视频| 日韩人妻精品一区2区三区| 国产精品三级大全| 蜜桃国产av成人99| videosex国产| 亚洲精品色激情综合| 日本欧美国产在线视频| 午夜免费男女啪啪视频观看| 波多野结衣一区麻豆| 一区二区三区四区激情视频| 女人被躁到高潮嗷嗷叫费观| 亚洲经典国产精华液单| 亚洲精品色激情综合| 在现免费观看毛片| 亚洲精品久久成人aⅴ小说| 欧美精品一区二区大全| 欧美精品av麻豆av| 亚洲精品国产av蜜桃| 看十八女毛片水多多多| 精品久久蜜臀av无| 国产乱来视频区| 午夜免费男女啪啪视频观看| 欧美精品一区二区免费开放| 一区二区三区四区激情视频| 亚洲美女搞黄在线观看| 亚洲国产精品一区二区三区在线| 中文字幕最新亚洲高清| 丝瓜视频免费看黄片| 国产成人精品在线电影| 一级片'在线观看视频| a 毛片基地| 国产精品国产av在线观看| 国产精品国产三级专区第一集| 欧美性感艳星| 高清在线视频一区二区三区| 午夜福利在线观看免费完整高清在| 久久久国产一区二区| 丝袜脚勾引网站| 蜜桃在线观看..| 一区二区三区精品91| 天美传媒精品一区二区| 国产在线一区二区三区精| 久久久久久久久久人人人人人人| 久久人人爽人人爽人人片va| 久久精品久久久久久久性| 99re6热这里在线精品视频| 一区在线观看完整版| 一二三四中文在线观看免费高清| 一个人免费看片子| 亚洲精品久久久久久婷婷小说| 久久影院123| 如何舔出高潮| 亚洲av电影在线进入| 成人国语在线视频| 97精品久久久久久久久久精品| 亚洲国产精品999| 午夜av观看不卡| 日日摸夜夜添夜夜爱| 乱人伦中国视频| 香蕉丝袜av| 亚洲精品色激情综合| 日韩 亚洲 欧美在线| 亚洲,欧美精品.| 亚洲天堂av无毛| 嫩草影院入口| 80岁老熟妇乱子伦牲交| 中文字幕人妻熟女乱码| 成年人午夜在线观看视频| 少妇的逼水好多| 国产精品一国产av| 99久国产av精品国产电影| 两个人免费观看高清视频| 亚洲精品久久久久久婷婷小说| 久久青草综合色| 久久精品国产亚洲av天美| 在线天堂最新版资源| 国产男女内射视频| 啦啦啦啦在线视频资源| 熟女电影av网| 欧美性感艳星| 亚洲欧美色中文字幕在线| 亚洲激情五月婷婷啪啪| 在线观看免费高清a一片| 一个人免费看片子| 中文欧美无线码| 高清不卡的av网站| 巨乳人妻的诱惑在线观看| 国产精品国产三级专区第一集| 日日爽夜夜爽网站| 久热这里只有精品99| 在线观看免费日韩欧美大片| av电影中文网址| 少妇被粗大的猛进出69影院 | 国产精品嫩草影院av在线观看| 国产精品人妻久久久久久| 欧美日韩视频精品一区| 老司机亚洲免费影院| 亚洲欧美成人精品一区二区| 国产欧美日韩综合在线一区二区| 国内精品宾馆在线| 又黄又爽又刺激的免费视频.| 伦精品一区二区三区| 成年人免费黄色播放视频| 婷婷色麻豆天堂久久| 国产激情久久老熟女| 亚洲国产精品成人久久小说| 亚洲美女搞黄在线观看| 大香蕉久久网| 日本vs欧美在线观看视频| 丰满乱子伦码专区| av国产精品久久久久影院| 欧美日韩国产mv在线观看视频| www日本在线高清视频| 妹子高潮喷水视频| 国产高清三级在线| 高清视频免费观看一区二区| 久久免费观看电影| 日韩不卡一区二区三区视频在线| 久久久久人妻精品一区果冻| 最近手机中文字幕大全| 国产亚洲一区二区精品| 欧美精品高潮呻吟av久久| 亚洲成人av在线免费| 菩萨蛮人人尽说江南好唐韦庄| a级毛片黄视频| 欧美+日韩+精品| 欧美精品人与动牲交sv欧美| 天美传媒精品一区二区| 天天影视国产精品| 美女大奶头黄色视频| 久久久精品94久久精品| 亚洲伊人色综图| 国产成人精品久久久久久| 91国产中文字幕| 久久精品熟女亚洲av麻豆精品| 91精品三级在线观看| 母亲3免费完整高清在线观看 | 久久久国产欧美日韩av| 欧美激情 高清一区二区三区| 色94色欧美一区二区| 视频区图区小说| 国产片特级美女逼逼视频| 亚洲精品美女久久av网站| 国产欧美另类精品又又久久亚洲欧美| 热re99久久精品国产66热6| 色网站视频免费| 久久ye,这里只有精品| 伊人久久国产一区二区| 18在线观看网站| 一级毛片我不卡| 免费在线观看完整版高清| 免费在线观看黄色视频的| 在线 av 中文字幕| 夫妻性生交免费视频一级片| 欧美性感艳星| 午夜视频国产福利| 国产欧美另类精品又又久久亚洲欧美| 18禁观看日本| 中文欧美无线码| 大香蕉97超碰在线| 国产成人精品福利久久| 免费大片18禁| 赤兔流量卡办理| 色婷婷av一区二区三区视频| 欧美 日韩 精品 国产| 午夜福利网站1000一区二区三区| 好男人视频免费观看在线| 18+在线观看网站| 女人精品久久久久毛片| 在线亚洲精品国产二区图片欧美| 一区二区三区精品91| 99久久精品国产国产毛片| 捣出白浆h1v1| 久久久精品94久久精品| 精品亚洲成a人片在线观看| 超色免费av| 亚洲精品国产av成人精品| 国产亚洲精品第一综合不卡 | 一级毛片我不卡| 国产亚洲一区二区精品| 免费少妇av软件| 免费在线观看完整版高清| 伊人亚洲综合成人网| 少妇被粗大猛烈的视频| 国产视频首页在线观看| 一区二区三区乱码不卡18| 精品99又大又爽又粗少妇毛片| 欧美+日韩+精品| 亚洲色图 男人天堂 中文字幕 | 在线观看www视频免费| 国产色婷婷99| 欧美日本中文国产一区发布| 亚洲精品久久久久久婷婷小说| 欧美少妇被猛烈插入视频| 老司机亚洲免费影院| 欧美变态另类bdsm刘玥| 日本vs欧美在线观看视频| 另类精品久久| 欧美日韩成人在线一区二区| 久久精品夜色国产| 边亲边吃奶的免费视频| 2021少妇久久久久久久久久久| 精品人妻熟女毛片av久久网站| 夫妻午夜视频| 成年女人在线观看亚洲视频| 亚洲av福利一区| 久久精品国产亚洲av天美| 黄片无遮挡物在线观看| 侵犯人妻中文字幕一二三四区| 大陆偷拍与自拍| 人人澡人人妻人| 免费人成在线观看视频色| 亚洲国产av影院在线观看| 久久99热6这里只有精品| 亚洲欧美一区二区三区黑人 | 久久精品国产a三级三级三级| 蜜臀久久99精品久久宅男| 在线 av 中文字幕| 制服丝袜香蕉在线| 一区二区三区乱码不卡18| 在线天堂中文资源库| 亚洲欧洲精品一区二区精品久久久 | 久久久久精品久久久久真实原创| 麻豆精品久久久久久蜜桃| 高清黄色对白视频在线免费看| 美女内射精品一级片tv| 女性生殖器流出的白浆| 91在线精品国自产拍蜜月| 日韩欧美一区视频在线观看| 国产精品人妻久久久影院| 韩国精品一区二区三区 | 精品亚洲成国产av| 99香蕉大伊视频| 最近2019中文字幕mv第一页| 男女下面插进去视频免费观看 | 一级片免费观看大全| 大陆偷拍与自拍| 视频中文字幕在线观看| 国产精品国产av在线观看| 免费不卡的大黄色大毛片视频在线观看| 成人二区视频| 黄片无遮挡物在线观看| 黄网站色视频无遮挡免费观看| 免费观看av网站的网址| 性色avwww在线观看| videosex国产| 在线观看美女被高潮喷水网站| 内地一区二区视频在线| 国产成人精品无人区| 日韩不卡一区二区三区视频在线| 青青草视频在线视频观看| 欧美成人精品欧美一级黄| 性高湖久久久久久久久免费观看| 欧美亚洲 丝袜 人妻 在线| 亚洲色图综合在线观看| 国产乱来视频区| 国产精品久久久久久精品电影小说| 尾随美女入室| 草草在线视频免费看| 精品一区二区三区视频在线| 久久久久精品人妻al黑| 亚洲一级一片aⅴ在线观看| 老司机亚洲免费影院| 国产又色又爽无遮挡免| 免费高清在线观看视频在线观看| 大片电影免费在线观看免费| 国产精品99久久99久久久不卡 | 亚洲精品国产色婷婷电影| 久久人人爽人人片av| 久久精品aⅴ一区二区三区四区 | 久久国产精品男人的天堂亚洲 | 啦啦啦中文免费视频观看日本| 自线自在国产av| 国产成人精品在线电影| 国产精品国产三级国产av玫瑰| 亚洲欧美清纯卡通| 亚洲三级黄色毛片| 国产精品国产三级专区第一集| 丰满少妇做爰视频| 少妇的逼水好多| 成人国语在线视频| 国产又爽黄色视频| 午夜久久久在线观看| 免费黄频网站在线观看国产| 精品久久蜜臀av无| 国产成人精品一,二区| 国产成人免费观看mmmm| 99久久人妻综合| 日韩伦理黄色片| 久久精品久久久久久久性| 免费人妻精品一区二区三区视频| 欧美3d第一页| 美女主播在线视频| 晚上一个人看的免费电影| 亚洲熟女精品中文字幕| 中文字幕免费在线视频6| 日本午夜av视频| 免费看av在线观看网站| 少妇被粗大猛烈的视频| 成年av动漫网址| 日韩伦理黄色片| 久久久久久久国产电影| 视频区图区小说| 午夜精品国产一区二区电影| 国产一区二区在线观看av| 熟妇人妻不卡中文字幕| 国产精品秋霞免费鲁丝片| 久久精品国产亚洲av天美| 搡女人真爽免费视频火全软件| 内地一区二区视频在线| 国产成人免费观看mmmm| 久久国产精品男人的天堂亚洲 | 国产精品不卡视频一区二区| 综合色丁香网| 国产欧美亚洲国产| 亚洲高清免费不卡视频| 香蕉丝袜av| 欧美人与性动交α欧美软件 | 国产亚洲一区二区精品| 亚洲美女黄色视频免费看| 亚洲成人一二三区av| 成人黄色视频免费在线看| 欧美亚洲日本最大视频资源| 久久久久精品久久久久真实原创| av黄色大香蕉| 久久久国产一区二区| 亚洲精品久久成人aⅴ小说| 多毛熟女@视频| 国产成人午夜福利电影在线观看| 久久影院123| 免费少妇av软件| 久久精品久久久久久噜噜老黄| 最近的中文字幕免费完整| 久久精品国产亚洲av涩爱| 日韩中文字幕视频在线看片| 亚洲精品中文字幕在线视频| av免费观看日本| 午夜精品国产一区二区电影| 美女国产视频在线观看| 日本欧美视频一区| 国产片内射在线| 尾随美女入室| 校园人妻丝袜中文字幕| 人人妻人人澡人人爽人人夜夜| 欧美日韩国产mv在线观看视频| 成人亚洲欧美一区二区av| 99九九在线精品视频| 最近中文字幕2019免费版| av电影中文网址| 肉色欧美久久久久久久蜜桃| 黄色视频在线播放观看不卡| 国产精品国产三级专区第一集| 欧美丝袜亚洲另类| 少妇人妻 视频| 精品亚洲乱码少妇综合久久| 亚洲精品乱久久久久久| 在线观看免费视频网站a站| 男女国产视频网站| 免费看光身美女| 男男h啪啪无遮挡| 大陆偷拍与自拍| 亚洲av免费高清在线观看| 边亲边吃奶的免费视频| 韩国高清视频一区二区三区| 又大又黄又爽视频免费| 一级爰片在线观看| 亚洲欧洲精品一区二区精品久久久 | 国产精品女同一区二区软件| 岛国毛片在线播放| 中文字幕精品免费在线观看视频 | 国产精品久久久久久精品古装| 日韩中字成人| 亚洲少妇的诱惑av| 啦啦啦在线观看免费高清www| 久久婷婷青草| 国产成人av激情在线播放| 亚洲中文av在线| 日韩免费高清中文字幕av| 国产成人精品无人区| 亚洲精品成人av观看孕妇| 国产精品无大码| 亚洲美女搞黄在线观看| 男人操女人黄网站| 乱码一卡2卡4卡精品| 啦啦啦在线观看免费高清www| 夜夜爽夜夜爽视频| 久久ye,这里只有精品| 国产一区二区三区综合在线观看 | 久久久久久久久久久久大奶| 婷婷色av中文字幕| 精品一区二区三区四区五区乱码 | 亚洲综合色惰| 日韩欧美一区视频在线观看| 亚洲精品久久午夜乱码| 香蕉丝袜av| 日韩 亚洲 欧美在线| 日本av免费视频播放| 男人舔女人的私密视频| 欧美精品国产亚洲| 国产毛片在线视频| 久久热在线av| 国产精品.久久久| 在线看a的网站| 一二三四在线观看免费中文在 | 9热在线视频观看99| 免费人妻精品一区二区三区视频| 韩国精品一区二区三区 | 99国产精品免费福利视频| 美女内射精品一级片tv| 美女大奶头黄色视频| 成人二区视频| 色吧在线观看| 国产成人午夜福利电影在线观看| 久久国产亚洲av麻豆专区| 精品国产露脸久久av麻豆| 国产精品久久久久久精品电影小说| av在线观看视频网站免费| 免费黄网站久久成人精品| 国产精品国产三级国产av玫瑰| 两个人免费观看高清视频| 久久久久久久精品精品| 免费观看av网站的网址| 建设人人有责人人尽责人人享有的| 99香蕉大伊视频| 少妇精品久久久久久久| 国产日韩一区二区三区精品不卡| 国产免费现黄频在线看| 成年人午夜在线观看视频| 国产黄频视频在线观看| 亚洲精品一区蜜桃| 午夜久久久在线观看| 国产探花极品一区二区| 人妻人人澡人人爽人人| 18在线观看网站| 亚洲精品国产av蜜桃| 精品午夜福利在线看| 免费在线观看完整版高清| 9191精品国产免费久久| 视频在线观看一区二区三区| 高清毛片免费看| 亚洲综合色网址| 日韩成人伦理影院| 丰满少妇做爰视频| 最新的欧美精品一区二区| 内地一区二区视频在线| 国国产精品蜜臀av免费| 国产白丝娇喘喷水9色精品| 国产精品熟女久久久久浪| 91国产中文字幕| 新久久久久国产一级毛片| 男女国产视频网站| 波野结衣二区三区在线| 欧美xxxx性猛交bbbb| 日韩伦理黄色片| 亚洲三级黄色毛片| 91在线精品国自产拍蜜月| 制服丝袜香蕉在线| 久热久热在线精品观看| 大码成人一级视频| 久久久欧美国产精品| 狂野欧美激情性xxxx在线观看| 日韩成人伦理影院| 午夜福利影视在线免费观看| freevideosex欧美| 成人二区视频| 久久久久精品久久久久真实原创| 国产精品久久久av美女十八| 精品少妇久久久久久888优播| 永久网站在线| 飞空精品影院首页| 久久久国产一区二区| 亚洲国产日韩一区二区| 美女国产高潮福利片在线看| 26uuu在线亚洲综合色| 国产免费又黄又爽又色| av在线播放精品| 国产成人精品久久久久久| 国产成人欧美| 又粗又硬又长又爽又黄的视频| 欧美日本中文国产一区发布| 在线观看免费高清a一片| 久久久久久久久久人人人人人人| 精品久久久精品久久久| 狠狠婷婷综合久久久久久88av| 日韩免费高清中文字幕av| 狠狠婷婷综合久久久久久88av| www日本在线高清视频| 欧美xxⅹ黑人| 老司机影院成人| 永久网站在线| av女优亚洲男人天堂| 亚洲av国产av综合av卡| 黄色视频在线播放观看不卡| 亚洲精品久久久久久婷婷小说| 一级黄片播放器| 欧美日韩视频精品一区| 国产亚洲一区二区精品| 国产熟女午夜一区二区三区| 纯流量卡能插随身wifi吗| 99热网站在线观看| 国产深夜福利视频在线观看| av免费在线看不卡| 国产欧美日韩一区二区三区在线| 国产极品天堂在线| 青春草国产在线视频| 免费人妻精品一区二区三区视频| 国内精品宾馆在线| 日韩av在线免费看完整版不卡| 国产精品久久久久久av不卡| 久久久国产精品麻豆| 欧美国产精品一级二级三级| 国产精品久久久久久久久免| 日韩精品有码人妻一区| 女的被弄到高潮叫床怎么办| 汤姆久久久久久久影院中文字幕| 女人被躁到高潮嗷嗷叫费观| 街头女战士在线观看网站| 亚洲四区av| 亚洲一区二区三区欧美精品| 国产精品久久久久久久久免| 亚洲国产精品一区二区三区在线| 少妇熟女欧美另类| 伊人亚洲综合成人网| 下体分泌物呈黄色| 街头女战士在线观看网站| 日韩 亚洲 欧美在线| 一区在线观看完整版| 国产精品一二三区在线看| 亚洲一区二区三区欧美精品| 另类精品久久| 日本与韩国留学比较| 国产亚洲一区二区精品| 91国产中文字幕| 黄色 视频免费看| 日韩制服骚丝袜av| 国产亚洲精品第一综合不卡 | 熟妇人妻不卡中文字幕| 五月伊人婷婷丁香| 欧美成人午夜免费资源| 国产精品人妻久久久影院| 美女视频免费永久观看网站| 亚洲国产av影院在线观看| 一级爰片在线观看| 黑人猛操日本美女一级片| 18禁国产床啪视频网站| 日本爱情动作片www.在线观看| 性高湖久久久久久久久免费观看| 国产精品久久久久成人av| 人成视频在线观看免费观看| 麻豆乱淫一区二区| 久久久久精品久久久久真实原创| 国语对白做爰xxxⅹ性视频网站| 中国美白少妇内射xxxbb| 91在线精品国自产拍蜜月| 国产熟女欧美一区二区| 国产 精品1|