• <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欧美aⅴ国产| 悠悠久久av| 亚洲欧美一区二区三区久久| 伦理电影免费视频| 热99re8久久精品国产| 亚洲av第一区精品v没综合| 精品电影一区二区在线| 又黄又爽又免费观看的视频| 国产麻豆69| 国产精品久久久av美女十八| 国产精品亚洲一级av第二区| 天天躁日日躁夜夜躁夜夜| aaaaa片日本免费| 国产有黄有色有爽视频| av视频免费观看在线观看| 三上悠亚av全集在线观看| 亚洲专区国产一区二区| 久久影院123| 欧美人与性动交α欧美精品济南到| av有码第一页| 国产蜜桃级精品一区二区三区 | 丝袜美腿诱惑在线| 亚洲精品自拍成人| 在线观看免费午夜福利视频| 女人被狂操c到高潮| 一级毛片女人18水好多| 成年版毛片免费区| 人人妻,人人澡人人爽秒播| a级片在线免费高清观看视频| 亚洲九九香蕉| 99国产综合亚洲精品| 人成视频在线观看免费观看| 欧美日韩乱码在线| 亚洲男人天堂网一区| 国产欧美日韩精品亚洲av| 亚洲性夜色夜夜综合| 国产xxxxx性猛交| 一级黄色大片毛片| 欧美乱码精品一区二区三区| 亚洲午夜理论影院| 法律面前人人平等表现在哪些方面| 一夜夜www| 18禁美女被吸乳视频| 国产三级黄色录像| av超薄肉色丝袜交足视频| 国产欧美日韩综合在线一区二区| 精品久久蜜臀av无| 精品卡一卡二卡四卡免费| 黄色片一级片一级黄色片| 亚洲精品自拍成人| 青草久久国产| 欧美日韩成人在线一区二区| 天堂动漫精品| 亚洲少妇的诱惑av| 久久香蕉精品热| 国产精品乱码一区二三区的特点 | 欧美日韩精品网址| 一级毛片高清免费大全| 国产成人啪精品午夜网站| 国产精品美女特级片免费视频播放器 | x7x7x7水蜜桃| 一区二区三区精品91| 久久人妻福利社区极品人妻图片| 黄片小视频在线播放| 少妇粗大呻吟视频| 亚洲一区中文字幕在线| 精品卡一卡二卡四卡免费| 国产在线精品亚洲第一网站| 中出人妻视频一区二区| 别揉我奶头~嗯~啊~动态视频| 欧美在线一区亚洲| 视频区欧美日本亚洲| 黄片小视频在线播放| 国产真人三级小视频在线观看| 亚洲三区欧美一区| 亚洲午夜精品一区,二区,三区| 成人精品一区二区免费| 欧美大码av| 老司机午夜福利在线观看视频| 后天国语完整版免费观看| 夜夜爽天天搞| 欧美日韩亚洲综合一区二区三区_| 精品乱码久久久久久99久播| 欧美黑人精品巨大| 国产男靠女视频免费网站| 久久香蕉精品热| 国产精品免费视频内射| 超碰成人久久| 黄色a级毛片大全视频| 两个人看的免费小视频| 日韩欧美免费精品| 久久中文看片网| 欧美日本中文国产一区发布| 宅男免费午夜| 99re在线观看精品视频| 午夜精品国产一区二区电影| 亚洲成人免费电影在线观看| 免费在线观看视频国产中文字幕亚洲| 精品一区二区三区av网在线观看| 十八禁网站免费在线| 人妻久久中文字幕网| 成年女人毛片免费观看观看9 | 极品人妻少妇av视频| 中文字幕人妻熟女乱码| 黑丝袜美女国产一区| 欧美日韩乱码在线| 日韩视频一区二区在线观看| 捣出白浆h1v1| netflix在线观看网站| 天天躁日日躁夜夜躁夜夜| 九色亚洲精品在线播放| 热re99久久国产66热| 午夜日韩欧美国产| 91九色精品人成在线观看| 午夜日韩欧美国产| 一级毛片高清免费大全| 久久 成人 亚洲| 国产精品电影一区二区三区 | 精品久久久久久,| 一进一出好大好爽视频| 欧美黑人精品巨大| 精品一区二区三卡| 久久精品aⅴ一区二区三区四区| tube8黄色片| 在线视频色国产色| 亚洲美女黄片视频| 精品亚洲成国产av| 搡老岳熟女国产| 欧美黑人精品巨大| 国产精品久久久av美女十八| 巨乳人妻的诱惑在线观看| 精品久久久久久,| 日韩 欧美 亚洲 中文字幕| 日韩 欧美 亚洲 中文字幕| 国产真人三级小视频在线观看| 国产淫语在线视频| 精品一品国产午夜福利视频| 午夜福利乱码中文字幕| 欧美黄色片欧美黄色片| 久久精品国产清高在天天线| a级片在线免费高清观看视频| 国产91精品成人一区二区三区| 99热只有精品国产| 高清视频免费观看一区二区| а√天堂www在线а√下载 | 欧美乱色亚洲激情| 又黄又爽又免费观看的视频| 色婷婷久久久亚洲欧美| 黄色成人免费大全| 在线观看66精品国产| 午夜两性在线视频| 999精品在线视频| 老汉色∧v一级毛片| 欧美亚洲日本最大视频资源| 午夜两性在线视频| 丝袜在线中文字幕| 精品电影一区二区在线| 国产亚洲精品第一综合不卡| 国产日韩欧美亚洲二区| 波多野结衣av一区二区av| 涩涩av久久男人的天堂| 精品亚洲成a人片在线观看| 精品亚洲成国产av| 丁香欧美五月| 巨乳人妻的诱惑在线观看| 久久狼人影院| 国产精品久久久人人做人人爽| 国产又色又爽无遮挡免费看| 一区二区日韩欧美中文字幕| 免费在线观看完整版高清| 欧美日韩中文字幕国产精品一区二区三区 | 午夜福利乱码中文字幕| 国产91精品成人一区二区三区| 亚洲熟妇中文字幕五十中出 | 日本精品一区二区三区蜜桃| 成人亚洲精品一区在线观看| 亚洲人成77777在线视频| 99精国产麻豆久久婷婷| 黑丝袜美女国产一区| 精品国产乱码久久久久久男人| 日本a在线网址| 国产成人影院久久av| 亚洲一卡2卡3卡4卡5卡精品中文| 国产精品秋霞免费鲁丝片| 国产成人欧美在线观看 | 嫩草影视91久久| 成人国产一区最新在线观看| 国产男靠女视频免费网站| 国产精品久久久av美女十八| 在线永久观看黄色视频| 日本黄色日本黄色录像| 国产男女超爽视频在线观看| 午夜福利在线免费观看网站| 午夜亚洲福利在线播放| 51午夜福利影视在线观看| 99re6热这里在线精品视频| 亚洲精品美女久久av网站| 国产精品一区二区在线观看99| 成人亚洲精品一区在线观看| 色94色欧美一区二区| 国产精品自产拍在线观看55亚洲 | 中文字幕最新亚洲高清| 美女福利国产在线| 一级黄色大片毛片| 美女午夜性视频免费| 少妇粗大呻吟视频| 一进一出好大好爽视频| 午夜福利免费观看在线| 久久久国产精品麻豆| 狠狠狠狠99中文字幕| tocl精华| 成人影院久久| 两性午夜刺激爽爽歪歪视频在线观看 | 天堂俺去俺来也www色官网| 国产不卡一卡二| 十八禁网站免费在线| 久久99一区二区三区| 最新的欧美精品一区二区| 国产熟女午夜一区二区三区| 欧美国产精品一级二级三级| 男男h啪啪无遮挡| 国产精品国产高清国产av | 成人亚洲精品一区在线观看| 在线观看免费日韩欧美大片| 欧美日韩亚洲国产一区二区在线观看 | 18禁黄网站禁片午夜丰满| e午夜精品久久久久久久| 中文字幕最新亚洲高清| 好看av亚洲va欧美ⅴa在| 亚洲欧美精品综合一区二区三区| 极品少妇高潮喷水抽搐| 国产精品国产av在线观看| 岛国毛片在线播放| 老鸭窝网址在线观看| 国产极品粉嫩免费观看在线| 可以免费在线观看a视频的电影网站| 久久青草综合色| 亚洲午夜理论影院| 侵犯人妻中文字幕一二三四区| 乱人伦中国视频| 国产欧美日韩一区二区三| 欧美日韩亚洲综合一区二区三区_| 国产精品成人在线| 男女高潮啪啪啪动态图| 好看av亚洲va欧美ⅴa在| 妹子高潮喷水视频| videos熟女内射| 少妇被粗大的猛进出69影院| 五月开心婷婷网| 一级片免费观看大全| 麻豆国产av国片精品| 女人被狂操c到高潮| 国产成人av激情在线播放| 777久久人妻少妇嫩草av网站| 午夜福利视频在线观看免费| 午夜福利,免费看| 精品国产乱码久久久久久男人| 欧美乱妇无乱码| 久久久精品国产亚洲av高清涩受| www.熟女人妻精品国产| 亚洲午夜理论影院| 黑人巨大精品欧美一区二区mp4| 国产在线一区二区三区精| a级片在线免费高清观看视频| 十八禁网站免费在线| 老司机影院毛片| 欧美成狂野欧美在线观看| 国产精品免费大片| 窝窝影院91人妻| 一进一出抽搐gif免费好疼 | av线在线观看网站| 18禁国产床啪视频网站| 一区福利在线观看| 国产一区二区三区在线臀色熟女 | 一a级毛片在线观看| 久久久久久久午夜电影 | 露出奶头的视频| 亚洲欧美一区二区三区黑人| 超色免费av| 一区在线观看完整版| 久久久久国产一级毛片高清牌| 国产成人av激情在线播放| 国产精品98久久久久久宅男小说| 新久久久久国产一级毛片| 久久亚洲精品不卡| 大陆偷拍与自拍| 国产亚洲精品久久久久5区| 亚洲av第一区精品v没综合| 中文字幕人妻熟女乱码| 国产区一区二久久| 久久青草综合色| 侵犯人妻中文字幕一二三四区| 自线自在国产av| 午夜福利乱码中文字幕| av有码第一页| 久久99一区二区三区| 国产av又大| 国产激情久久老熟女| 校园春色视频在线观看| 精品欧美一区二区三区在线| 亚洲精品久久午夜乱码| 曰老女人黄片| 午夜影院日韩av| 国产欧美日韩精品亚洲av| 狠狠狠狠99中文字幕| 在线观看日韩欧美| 久久九九热精品免费| 国产精品 欧美亚洲| 国产97色在线日韩免费| 午夜福利在线观看吧| 校园春色视频在线观看| 人妻丰满熟妇av一区二区三区 | 亚洲精品一卡2卡三卡4卡5卡| ponron亚洲| 亚洲av成人av| 男人的好看免费观看在线视频 | 国产xxxxx性猛交| 人人妻人人澡人人爽人人夜夜| 99精品久久久久人妻精品| 午夜免费鲁丝| 999久久久精品免费观看国产| 久久狼人影院| 99精国产麻豆久久婷婷| 亚洲五月婷婷丁香| 久久久久精品国产欧美久久久| 老熟妇乱子伦视频在线观看| 欧美午夜高清在线| 啪啪无遮挡十八禁网站| 老司机午夜福利在线观看视频| 王馨瑶露胸无遮挡在线观看| 国产成人系列免费观看| 女人久久www免费人成看片| 亚洲国产精品一区二区三区在线| 如日韩欧美国产精品一区二区三区| 高清av免费在线| 亚洲国产精品合色在线| 亚洲精品中文字幕在线视频| 久久狼人影院| 国产亚洲欧美在线一区二区| 久久精品成人免费网站| 亚洲avbb在线观看| 亚洲专区国产一区二区| 国产免费av片在线观看野外av| 国产精品久久久久久人妻精品电影| 国产国语露脸激情在线看| 精品国产乱子伦一区二区三区| 国产99白浆流出| 亚洲精品粉嫩美女一区| 亚洲人成电影免费在线| 日日摸夜夜添夜夜添小说| 亚洲精品国产区一区二| 亚洲片人在线观看| 中亚洲国语对白在线视频| 国产成人av教育| 精品高清国产在线一区| av中文乱码字幕在线| 精品国产一区二区三区久久久樱花| 国产精品秋霞免费鲁丝片| 啦啦啦在线免费观看视频4| 成人特级黄色片久久久久久久| 一a级毛片在线观看| 黄频高清免费视频| 老司机午夜十八禁免费视频| 熟女少妇亚洲综合色aaa.| 久久香蕉精品热| 久久精品国产亚洲av香蕉五月 | 午夜成年电影在线免费观看| 高清欧美精品videossex| 欧美日韩亚洲高清精品| 日本黄色日本黄色录像| 日韩大码丰满熟妇| 中出人妻视频一区二区| 久久久国产一区二区| 9热在线视频观看99| 波多野结衣av一区二区av| 亚洲av日韩在线播放| 欧美成人免费av一区二区三区 | 天天躁夜夜躁狠狠躁躁| 成年人免费黄色播放视频| 99热国产这里只有精品6| 窝窝影院91人妻| 丝瓜视频免费看黄片| 久久中文字幕人妻熟女| 色综合欧美亚洲国产小说| 亚洲自偷自拍图片 自拍| 国产精品香港三级国产av潘金莲| 99热国产这里只有精品6| 一边摸一边抽搐一进一出视频| 在线观看免费视频网站a站| 后天国语完整版免费观看| 色老头精品视频在线观看| 久久精品国产清高在天天线| 成人手机av| 欧美最黄视频在线播放免费 | 亚洲色图综合在线观看| 国产深夜福利视频在线观看| 俄罗斯特黄特色一大片| 日韩一卡2卡3卡4卡2021年| 亚洲午夜理论影院| 好男人电影高清在线观看| 亚洲熟妇熟女久久| 欧美成狂野欧美在线观看| 精品第一国产精品| 少妇 在线观看| 久久久精品免费免费高清| 亚洲黑人精品在线| 国产aⅴ精品一区二区三区波| 亚洲精品成人av观看孕妇| 热99久久久久精品小说推荐| 在线观看免费日韩欧美大片| 国产激情欧美一区二区| 嫩草影视91久久| 久久中文字幕人妻熟女| 99久久综合精品五月天人人| 欧美日韩精品网址| 在线观看免费视频网站a站| 最近最新免费中文字幕在线| 男男h啪啪无遮挡| 国产在线观看jvid| 中文欧美无线码| 黑人操中国人逼视频| 丰满的人妻完整版| 国产国语露脸激情在线看| 欧美日韩瑟瑟在线播放| 日韩欧美一区二区三区在线观看 | av网站在线播放免费| 亚洲精品国产区一区二| 精品国产乱码久久久久久男人| 亚洲精品自拍成人| 欧美午夜高清在线| 欧美精品一区二区免费开放| 久久国产亚洲av麻豆专区| av欧美777| 国产精品九九99| 久久久久久久午夜电影 | 90打野战视频偷拍视频| 国产成人欧美在线观看 | 天堂动漫精品| 99久久国产精品久久久| 成年人午夜在线观看视频| 午夜福利乱码中文字幕| 欧美在线一区亚洲| 丝袜美足系列| 十八禁高潮呻吟视频| 亚洲精品国产精品久久久不卡| 在线av久久热| 激情视频va一区二区三区| 大香蕉久久网| 国产精品亚洲av一区麻豆| 又黄又爽又免费观看的视频| videosex国产| 久久人妻av系列| 成人av一区二区三区在线看| 午夜影院日韩av| 性色av乱码一区二区三区2| 制服人妻中文乱码| 国产欧美日韩一区二区三区在线| 亚洲精品粉嫩美女一区| 欧美丝袜亚洲另类 | 色精品久久人妻99蜜桃| 精品国内亚洲2022精品成人 | 亚洲精品美女久久av网站| 99热网站在线观看| 麻豆乱淫一区二区| 久久人妻福利社区极品人妻图片| av一本久久久久| 丝瓜视频免费看黄片| 母亲3免费完整高清在线观看| 91老司机精品| 精品一区二区三卡| 国产欧美日韩综合在线一区二区| www日本在线高清视频| 精品乱码久久久久久99久播| 久久久久久亚洲精品国产蜜桃av| 一区二区三区激情视频| 亚洲七黄色美女视频| 国精品久久久久久国模美| 国产av一区二区精品久久| 如日韩欧美国产精品一区二区三区| 国产精品乱码一区二三区的特点 | 国产aⅴ精品一区二区三区波| 十八禁高潮呻吟视频| 18在线观看网站| 亚洲欧美一区二区三区久久| 亚洲第一青青草原| 高清av免费在线| 91字幕亚洲| 超色免费av| 免费日韩欧美在线观看| 捣出白浆h1v1| 男女高潮啪啪啪动态图| 涩涩av久久男人的天堂| 如日韩欧美国产精品一区二区三区| 亚洲自偷自拍图片 自拍| 亚洲av成人av| 老司机福利观看| 伊人久久大香线蕉亚洲五| 一级,二级,三级黄色视频| 在线国产一区二区在线| 91麻豆精品激情在线观看国产 | 在线观看午夜福利视频| 亚洲午夜理论影院| 亚洲av成人av| 欧美日韩亚洲国产一区二区在线观看 | √禁漫天堂资源中文www| 嫩草影视91久久| 久久精品成人免费网站| 无限看片的www在线观看| 精品人妻熟女毛片av久久网站| 久久草成人影院| 91字幕亚洲| 男人的好看免费观看在线视频 | 人妻 亚洲 视频| 亚洲人成伊人成综合网2020| 精品少妇一区二区三区视频日本电影| 国产精品自产拍在线观看55亚洲 | videos熟女内射| tocl精华| 久久久久国产精品人妻aⅴ院 | 成熟少妇高潮喷水视频| 精品电影一区二区在线| 又大又爽又粗| netflix在线观看网站| 久久香蕉激情| 下体分泌物呈黄色| 久久国产精品大桥未久av| 王馨瑶露胸无遮挡在线观看| 91成年电影在线观看| 91老司机精品| 国产xxxxx性猛交| 日本vs欧美在线观看视频| 成人av一区二区三区在线看| 麻豆av在线久日| 成人永久免费在线观看视频| 国产精品自产拍在线观看55亚洲 | 久久久久精品人妻al黑| 欧美久久黑人一区二区| 国产精品国产av在线观看| 婷婷丁香在线五月| 国产日韩欧美亚洲二区| 免费在线观看黄色视频的| 精品电影一区二区在线| 亚洲一卡2卡3卡4卡5卡精品中文| 国产精品综合久久久久久久免费 | 免费久久久久久久精品成人欧美视频| 免费看a级黄色片| 亚洲人成伊人成综合网2020| 身体一侧抽搐| 国产精品久久视频播放| 色老头精品视频在线观看| 香蕉丝袜av| 一级毛片精品| 天堂√8在线中文| av天堂久久9| 国产蜜桃级精品一区二区三区 | 建设人人有责人人尽责人人享有的| 99精国产麻豆久久婷婷| 国产1区2区3区精品| 日本精品一区二区三区蜜桃| 中文字幕最新亚洲高清| 国产淫语在线视频| 亚洲欧美激情综合另类| 每晚都被弄得嗷嗷叫到高潮| av超薄肉色丝袜交足视频| 美国免费a级毛片| 99国产综合亚洲精品| 一级毛片女人18水好多| 欧美激情 高清一区二区三区| 成熟少妇高潮喷水视频| 搡老岳熟女国产| 亚洲综合色网址| 搡老岳熟女国产| 色综合欧美亚洲国产小说| 精品一区二区三卡| 国产激情久久老熟女| 99re在线观看精品视频| 欧美另类亚洲清纯唯美| 亚洲片人在线观看| 午夜成年电影在线免费观看| 精品视频人人做人人爽| 久久久久精品国产欧美久久久| 他把我摸到了高潮在线观看| 久久久久精品国产欧美久久久| 两个人免费观看高清视频| 亚洲欧洲精品一区二区精品久久久| 国产精品国产av在线观看| 成人18禁高潮啪啪吃奶动态图| 亚洲精品美女久久久久99蜜臀| 99国产精品99久久久久| 黄色视频,在线免费观看| 成人永久免费在线观看视频| 一个人免费在线观看的高清视频| 久久精品亚洲精品国产色婷小说| 天堂中文最新版在线下载| 免费日韩欧美在线观看| 视频在线观看一区二区三区| 国产无遮挡羞羞视频在线观看| 国产欧美日韩综合在线一区二区| 久久人人97超碰香蕉20202| 国产97色在线日韩免费| 精品国产一区二区三区四区第35| 男女高潮啪啪啪动态图| 亚洲熟妇中文字幕五十中出 | 精品少妇久久久久久888优播|