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

    柴達(dá)木盆地沙丘移動的空間分異及對形態(tài)參數(shù)的響應(yīng)

    2021-06-30 07:32:50李建軍焦菊英白雷超陳同德嚴(yán)晰芹祁泓錕
    農(nóng)業(yè)工程學(xué)報 2021年7期
    關(guān)鍵詞:柴達(dá)木盆地沙丘速度

    李建軍,焦菊英,,曹 雪,白雷超,陳同德,嚴(yán)晰芹,祁泓錕

    ·研究速報·

    柴達(dá)木盆地沙丘移動的空間分異及對形態(tài)參數(shù)的響應(yīng)

    李建軍1,焦菊英1,2※,曹 雪2,白雷超1,陳同德1,嚴(yán)晰芹2,祁泓錕1

    (1. 西北農(nóng)林科技大學(xué)水土保持研究所,黃土高原土壤侵蝕與旱地農(nóng)業(yè)國家重點實驗室,楊凌 712100; 2. 中國科學(xué)院水利部水土保持研究所,楊凌 712100)

    研究沙丘移動規(guī)律與空間分異特征有助于準(zhǔn)確了解風(fēng)沙活動的危害程度,為區(qū)域城鎮(zhèn)規(guī)劃和青藏高原國家生態(tài)安全屏障保護(hù)與建設(shè)提供參考。該研究基于Google Earth軟件測量橫向沙丘的形態(tài)參數(shù)與移動速度,分析了移動速度與方向的空間分異特征,以及移動速度對形態(tài)參數(shù)的響應(yīng)。結(jié)果表明,1)柴達(dá)木盆地內(nèi)沙丘移動速度介于0~23.53 m/a之間,平均4.66 m/a,以中速(1~5 m/a)為主(53.73%),中部移動最快(5.22±4.27 m/a),東南部最慢(3.27±3.08 m/a),移動方向與主風(fēng)向一致;2)沙丘移動速度與寬度的相關(guān)性最好(2=0.988),以后的研究中需多關(guān)注沙丘寬度這一參數(shù);3)從柴達(dá)木盆地整體看,降水量越大、風(fēng)速越低,沙丘移動越慢;從局部看,植被覆蓋度越大,沙丘移動越慢;4)密集的河流可阻擋沙丘移動,保護(hù)格爾木市區(qū),但不能消除沙丘移動的威脅。

    遙感;沙丘;移動;形態(tài)參數(shù);柴達(dá)木盆地

    0 引 言

    沙丘是風(fēng)沙作用的產(chǎn)物,是沙漠、沙地的主要地貌類型,橫向沙丘是指沙丘走向與起沙風(fēng)合成風(fēng)向交角大于60°的丘狀風(fēng)積地貌形態(tài),包括新月形沙丘及沙丘鏈、格狀沙丘、拋物線形沙丘、復(fù)合型新月形沙丘等[1]。沙丘移動是指沙丘整體性順風(fēng)向位移,過程為迎風(fēng)坡風(fēng)蝕,背風(fēng)坡堆積。沙丘移動會破壞道路、農(nóng)田、草場[2],是除沙漠?dāng)U散、沙塵暴之外的一種主要的沙漠化運(yùn)動形式,也是研究沙丘形成過程的重要問題之一[3]。沙丘移動現(xiàn)象廣泛發(fā)生在干旱區(qū)、海岸甚至是外星球[4]。沙丘移動研究主要包括移動速度的年內(nèi)年際變化[5-6]、移動速度的影響因素[7-8]、移動前后的形態(tài)變化以及監(jiān)測方法的改進(jìn)[9]等方面。由于形態(tài)簡單且較為典型,沙丘移動研究多集中在新月形沙丘上[7]。

    地面監(jiān)測和遙感影像監(jiān)測是沙丘移動研究中的常用方法。地面監(jiān)測精度高,但野外監(jiān)測環(huán)境惡劣,實地監(jiān)測的范圍往往較小[10]。如今,差分GPS和無人機(jī)攝影[11]等新興測量技術(shù)降低了地面監(jiān)測工作難度,各種高分影像的出現(xiàn)也使得遙感影像監(jiān)測有了可靠的數(shù)據(jù)來源[12]。相比于地面監(jiān)測,Google Earth提供的多時相高分辨率影像更適于大范圍的沙丘移動調(diào)查。因此,近幾年的沙丘移動研究中,大多數(shù)使用遙感影像測量[7]。其中,Google Earth提供的影像分辨率最高也最容易獲取,其測量精度也得到了廣泛驗證[13]。此外,Google Earth提供的影像還被廣泛用于驗證Landsat、Sentinel等中分辨率影像的沙丘移動測量結(jié)果[14]。

    柴達(dá)木盆地是中國重要的資源工業(yè)基地,是青藏高原沙漠分布的主要地區(qū)之一,由于空氣稀薄、氣溫較低,被認(rèn)為是“地球上與火星最為接近的地區(qū)”,是研究火星風(fēng)沙地貌的重要參照[15]。中國火星探測器“天問一號”已經(jīng)抵達(dá)火星軌道,對柴達(dá)木盆地風(fēng)沙地貌的基礎(chǔ)研究也需要跟進(jìn)。柴達(dá)木盆地已有的沙丘地貌研究主要集中在沙丘景觀格局和沉積物的粒度特征等方面[16],對柴達(dá)木盆地沙丘移動特征及規(guī)律還不甚明確。因此,本文以新月形沙丘等橫向沙丘為研究對象,基于Google Earth軟件對盆地內(nèi)的沙丘移動情況進(jìn)行測量,分析沙丘移動速度與方向的空間分異特征,以及沙丘移動速度對形態(tài)參數(shù)的響應(yīng),以期為區(qū)域城鎮(zhèn)規(guī)劃與生態(tài)文明建設(shè),以及青藏高原國家生態(tài)安全屏障保護(hù)與建設(shè)提供科學(xué)依據(jù)。

    1 材料與方法

    1.1 研究區(qū)概況

    柴達(dá)木盆地地處青藏高原東北部,介于87°48′E~99°20′E、34°41′N~39°20′N之間,屬大陸性荒漠氣候。盆地內(nèi)風(fēng)力和降水時空差異較大[17],河流以季節(jié)性河流為主,湖泊數(shù)量較多,沙漠零星分布。經(jīng)遙感解譯,沙漠總面積約16 846 km2。沙丘類型以新月形沙丘及新月形沙丘鏈為主,另外還有格狀沙丘、復(fù)合新月形沙丘、圓頂沙丘、拋物線形沙丘、沙壟、金字塔沙丘、爬坡沙丘等。

    格爾木市是海西蒙古族藏族自治州下轄縣級市,地處柴達(dá)木盆地南緣,是連接西部幾省區(qū)和絲綢之路經(jīng)濟(jì)帶的重要交通樞紐。市區(qū)海拔2 820 m,建在洪積扇邊緣,南部為戈壁灘,北部為濕地,因周圍“河流密集”而得名。本研究新月形沙丘典型區(qū)選在格爾木市區(qū)東部約10 km的戈壁灘上,該區(qū)域有多期遙感影像,且地形平坦,對風(fēng)況的影響較小,適合獨(dú)立地研究沙丘形態(tài)參數(shù)對移動速度的影響。

    1.2 數(shù)據(jù)來源與沙丘選取

    氣象數(shù)據(jù)來自海西州氣象局提供的柴達(dá)木盆地內(nèi)7個氣象站2003-2016年逐小時數(shù)據(jù),計算了平均風(fēng)速、起沙風(fēng)時數(shù)、降水量、蒸發(fā)量等氣象要素的多年平均值(表1)。

    遙感影像主要是Google Earth提供的GeoEye(0.41 m分辨率)、Quickbird(0.6 m分辨率)和World View(0.5 m分辨率)影像。基于高分影像對柴達(dá)木盆地沙漠分布進(jìn)行目視解譯,用0.1°×0.1°網(wǎng)格均勻選點,并剔除兩期影像間隔時間短于1 a或僅一期影像的樣點。為保證可比性,對于樣點落在新月形沙丘鏈和格狀沙丘的區(qū)域,均選取其中一環(huán)作為研究對象,共測量沙丘67個(圖1),影像時間間隔最短1.47 a,最長15.5 a,平均間隔6.45 a,大部分都在3 a以上。為了排除環(huán)境因素的影響,進(jìn)一步研究新月形沙丘移動速度和形態(tài)參數(shù)的關(guān)系,選取了盆地中南部格爾木市區(qū)東部一處典型區(qū),共測量沙丘34個(圖1c)的形態(tài)參數(shù)和移動情況,影像時間間隔5.5 a。

    表1 柴達(dá)木盆地基本氣象要素(2003-2016年)

    注:起沙風(fēng)時數(shù)為一年中風(fēng)速大于6 m·s-1的累計小時數(shù)。

    Note: Sand driving wind hours is the cumulative hours of wind velocity greater than 6 m·s-1in a year.

    注:圖a中箭頭指向沙丘的移動方向。下同。

    Note: Arrows indicate the direction of dune migration in the picture a. Same as below.

    圖1 研究區(qū)地理位置及其周邊環(huán)境

    Fig.1 Research area and its surrounding environment

    1.3 沙丘形態(tài)參數(shù)提取

    沙丘主要形態(tài)參數(shù)包括迎風(fēng)坡長度、背風(fēng)坡長度、寬度、高度、底面積、周長等,主要通過Google Earth軟件獲取[12]。迎風(fēng)坡長度、背風(fēng)坡長度、寬度3個參數(shù)在最高分辨率影像正視狀態(tài)下,由軟件提供的工具測量沙丘移動前的平面長度。由于Google Earth高程精度低,故沙丘高度用背風(fēng)坡長度乘以休止角的正切計算[8]。在野外考察中,測得格爾木附近沙丘休止角為32°。底面積和周長需要勾繪出沙丘移動前的輪廓,保存后在屬性中查看。本研究引入了形狀系數(shù)的概念,即沙丘實際周長與?同面積?圓的周長之?比,可以反映沙丘的形態(tài);長寬比由沙丘迎風(fēng)坡、背風(fēng)坡長度之和除以寬度得到,可以反映沙丘的“胖瘦”[8]。

    1.4 沙丘移動測量

    沙丘移動前后的影像介于1999-2017年之間,大部分集中在2003-2016年。首先勾繪沙丘移動前脊線和輪廓,再切換到移動后影像進(jìn)行測量。測量點選擇沙丘兩翼和脊線中部,分別測量3個點的移動距離和方向并取平均值。由移動距離與兩期影像間隔時長計算沙丘移動速度。對影像偏移的沙丘,借助附近的巖石、灌叢等固定物體進(jìn)行校正。

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

    根據(jù)朱震達(dá)等[18]的分級方法,將柴達(dá)木盆地沙丘移動速度分為4個等級:慢速<1 m/a、中速1~5 m/a、快速>5~10 m/a、特別快速>10 m/a(圖1a)。用SPSS 20軟件對典型區(qū)新月形沙丘各形態(tài)參數(shù)和移動速度進(jìn)行相關(guān)性分析。

    2 結(jié)果與分析

    2.1 沙丘移動速度與方向

    柴達(dá)木盆地沙漠沙丘移動速度在0~23.53 m/a之間,平均速度4.66 m/a。盆地內(nèi)沙丘移動以中速為主(53.73%),其次是快速(29.85%),慢速(10.45%)和特別快速(5.97%)。在空間上自西北向東南減慢,但無明顯趨勢(表2)。盆地中部沙丘由于地勢平坦,移動速度較快(5.22 m/a)。西南部和北部沙漠受到地形的影響,移動速度偏慢(3.35、3.49 m/a),但略快于降水較多植被較好的東南部(3.27 m/a)。盆地內(nèi)沙丘移動方向整體上自西北向東南,在53.52°~150.06°之間,平均104.13°,與主風(fēng)向基本一致。另外,盆地邊緣沙丘又受到地形影響,沿著高大山脈的延伸方向前進(jìn)。

    表2 柴達(dá)木盆地各區(qū)域沙丘移動情況

    注:同列不同字母表示區(qū)域間差異顯著,<0.01。

    Note: Different letters show significant difference,<0.01.

    2.2 沙丘移動的空間分異特征

    西北部沙漠沿山脈分為3段(圖2a)。沙丘類型包括格狀、新月形和復(fù)合新月形沙丘。西段為峽谷出口,風(fēng)力強(qiáng)勁,沙丘平均移動速度6.57 m/a;中段有高大山脈阻擋,平均移動速度3.13 m/a;東段位于茫崖西北風(fēng)口的下風(fēng)向,平均移動速度4.04 m/a。中部沙漠分為3條沙帶(圖2b)。北部沙帶地處冷湖西北風(fēng)口下方,以小型沙丘為主,平均移動速度10.35 m/a。中部沙帶地處濕地中心,以新月形沙丘鏈為主,平均移動速度4.84 m/a。南部沙帶位于戈壁,以格狀和新月形沙丘為主,平均移動速度3.76 m/a。東南部以固定和半固定格狀沙丘為主,臨近山體的沙丘平均移動速度僅0.19 m/a,而西側(cè)4個沙丘達(dá)到了6.5 m/a(圖2c)。沙漠已緊鄰都蘭縣城,但由于降水量較高,沙丘基本固定。北部與西南部沙漠都分布在山谷中(圖2d, 2e),北部以新月形沙丘和新月形沙丘鏈為主,西南部以格狀和復(fù)合新月形沙丘為主。由于山脈削弱了風(fēng)力,沙丘平均移動速度分別為3.49和3.35 m/a。

    a.西北部;b.中部;c.東南部;c1裸露的沙丘;c2植被覆蓋的沙丘;d.北部;e.西南部

    a. Northwest region; b. Central region; c. Southeast region; c1.Bare dune; c2.Vegetation-covered dune; d. North region; e. Southwest region

    圖2 柴達(dá)木盆地各區(qū)域沙丘移動情況

    Fig.2 Dune migration in various regions of the Qaidam Basin

    2.3 沙丘移動速度與形態(tài)參數(shù)的關(guān)系

    格爾木市區(qū)東部典型區(qū)所測量的新月形沙丘迎風(fēng)坡長度范圍18.8~172.5 m,寬度范圍18~282 m,移動速度在2.79~23.53 m/a之間,整體移動方向95°,平均速度10.94 m/a,以非常快速為主(55.88%),快速(23.53%)和中速(20.59%)占比相近。典型區(qū)沙丘迎風(fēng)坡長度、背風(fēng)坡長度、寬度、高度、周長之間呈極顯著線性相關(guān),底面積與迎風(fēng)坡長度、背風(fēng)坡長度、寬度、高度、周長之間呈二次函數(shù)關(guān)系。背風(fēng)坡長度、寬度、周長、底面積之間的關(guān)系較好(2≥0.729,圖3)。

    沙丘寬度與形狀系數(shù)呈冪函數(shù)關(guān)系(2=0.755),與長寬比呈冪函數(shù)關(guān)系(2=0.565)。這表明新月形沙丘由小變大的過程,沙丘變“矮胖”(長寬比減?。?,由橢圓形變?yōu)樾略滦危ㄐ螤钕禂?shù)增大)。

    新月形沙丘移動速度與沙丘迎風(fēng)坡長度(2=0.716)、背風(fēng)坡長度(2=0.705)、寬度(2=0.988)、高度(2=0.705)、周長(2=0. 986)、底面積(2=0.955)等形態(tài)參數(shù)之間都呈現(xiàn)指數(shù)關(guān)系,與形狀系數(shù)呈極顯著(<0.01)線性負(fù)相關(guān)關(guān)系(2=0.438),與長寬比呈顯著(<0.05)線性正相關(guān)關(guān)系(2=0.130)。其中,移動速度與沙丘寬度、底面積、周長的擬合關(guān)系較好(圖4c, 4e, 4f)。沙丘形狀系數(shù)越大,長寬比越小,沙丘體型越大,對應(yīng)的移動速度也越慢(圖4g)。該區(qū)域新月形沙丘左翼比右翼向前延伸地更長(圖1c),圖4h中更多的點分布在1:1線上方,左翼移動速度均值(11.16±5.26 m/a)略大于右翼(10.53±5.47 m/a)。

    3 討 論

    3.1 Google Earth提取數(shù)據(jù)可靠性

    Google Earth測量平面形態(tài)較為可靠。如劉佳等[13]使用Google Earth進(jìn)行農(nóng)作物面積調(diào)查,其誤差低于0.1%;李鎮(zhèn)等[19]與實地測量對比,發(fā)現(xiàn)QuickBird影像在目視解譯切溝面積、周長的平均相對誤差都在5%左右。相較于切溝,沙丘面積更大且與背景環(huán)境色差明顯,解譯沙丘平面形態(tài)參數(shù)的誤差在可接受范圍內(nèi),基于Google Earth測量沙丘形態(tài)參數(shù)與移動距離的方法在庫木塔格沙漠[7],庫魯克沙漠[8],毛烏素沙地[12],埃及西部沙漠[20]等地得到了廣泛應(yīng)用。Google Earth無法準(zhǔn)確測量沙丘高度,而背風(fēng)坡與背景亮度差異最為明顯,為準(zhǔn)確獲取沙丘高度,本研究中選擇了背風(fēng)坡水平長度換算的方法[8]。今后的研究可結(jié)合無人機(jī)、差分GPS等新興技術(shù)手段,對典型地區(qū)開展實地測量,以提高測量沙丘三維形態(tài)參數(shù)精度。

    3.2 環(huán)境因素對沙丘移動的影響

    在東南部,年降水量大于200 mm,沙丘大部分固定,植被覆蓋的沙丘平均移動速度僅0.19 m/a,而西側(cè)裸露的沙丘平均移動速度達(dá)6.50 m/a。與東南部相比,西北部和中心的年平均風(fēng)速均大于3 m/s[21],對應(yīng)的沙丘移動速度也更快??梢娚城鹨苿铀俣仁艿斤L(fēng)況、降水、植被、地形等環(huán)境因素的綜合影響。降水促進(jìn)植被生長,降低近地表風(fēng)速[7],進(jìn)而減慢沙丘移動。中部沙帶距格爾木市區(qū)38 km,前緣沙丘移動速度3.77 m/a。按此速度估算,該市將在萬年后被吞沒,但市區(qū)西北方向河流密集(圖1b),沙丘很難通過[22-23],短期內(nèi)威脅較小。考慮到沙物質(zhì)能以風(fēng)沙流形式輸移,遇到障礙物再度形成沙丘,故威脅仍存在。柴達(dá)木盆地氣候呈暖濕化趨勢,植被覆蓋度逐漸增加[24],沙丘移動將逐漸變慢。

    3.3 形態(tài)參數(shù)對沙丘移動的影響

    沙丘體型越大,移動速度越慢,這與前人研究一致;不同的是,在毛烏素沙地[12]、埃及西部沙漠[20]等地的研究中大多發(fā)現(xiàn)沙丘移動速度與形態(tài)參數(shù)呈線性相關(guān)。與前者相比本研究中典型區(qū)沙丘寬度范圍在18~282 m之間,數(shù)據(jù)離散程度較大,故呈指數(shù)關(guān)系,這與Yang等[7]的結(jié)論一致。沙丘高度與移動速度相關(guān)性高[10],在前人研究中最受關(guān)注。本研究中高度需換算,而寬度不僅易測量,且與移動速度相關(guān)性最高,建議給予更多關(guān)注。典型區(qū)新月形沙丘左翼移動較右翼略快,這是由兩個方向的起沙風(fēng)所致[17,25],主風(fēng)向(W)決定了沙丘移動方向,副風(fēng)向(WNW)則使得左翼伸長(圖1c);還存在盾形沙丘向圓頂沙丘的逆向演化,可能是沙源不足所致[12]。沙丘形態(tài)變化會改變移動特征[8],其移動前后的形態(tài)變化還需深入研究。

    4 結(jié) 論

    1)柴達(dá)木盆地沙丘移動速度在0~23.53 m/a之間,平均4.66 m/a,其中中速占53.73%,快速占29.85%,慢速占10.45%,特別快速占5.97%。沙丘移動方向整體上自西北向東南,在53.52°~150.06°之間,平均104.13°,與主風(fēng)向基本一致,并隨山脈走向變化。

    2)柴達(dá)木盆地沙丘移動速度整體上自西北向東南減慢,但無明顯趨勢。中部沙丘移動最快(5.22 m/a),之后依次是西北部(4.84 m/a)、北部(3.49 m/a)、西南部(3.35 m/a)和東南部(3.27 m/a)。

    3)典型區(qū)新月形沙丘移動速度與迎風(fēng)坡長度、背風(fēng)坡長度、寬度、高度、周長、底面積等形態(tài)參數(shù)之間都呈指數(shù)關(guān)系,與形狀系數(shù)呈極顯著(<0.01)線性負(fù)相關(guān)關(guān)系,與長寬比呈顯著(<0.05)正線性相關(guān)關(guān)系。沙丘寬度與移動速度相關(guān)性最好,以后的研究中應(yīng)多關(guān)注。

    4)從柴達(dá)木盆地整體來看,降水量越多、風(fēng)速越低,沙丘移動越慢;局部來看,植被覆蓋度越大,沙丘移動越慢。河流能在一定程度上阻擋沙丘移動,保護(hù)了格爾木市區(qū),但并不能消除沙丘移動帶來的威脅,沙丘移動問題仍需關(guān)注。

    [1] 吳正. 風(fēng)沙地貌學(xué)[M]. 北京:科學(xué)出版社,1987:98-166.

    [2] Islam Abou El-Magd, Osman H, Sayed A. Quantification of sand dune movements in the south western part of egypt, using remotely sensed data and GIS[J]. Journal of Geographic Information System, 2013, 5(5): 498-508.

    [3] Zhang Z C, Dong Z B, Qian G Q, et al. Formation and development of dunes in the northern Qarhan Desert, central Qaidam Basin, China[J]. Geological Journal, 2018, 53(3): 1123-1134.

    [4] Fenton, Lori K. Dune migration and slip face advancement in the Rabe Crater dune field, Mars[J]. Geophysical Research Letters, 2006, 33(20): L20201.

    [5] Kilibarda Z, Kilibarda V. Seasonal geomorphic processes and rates of sand movement at Mount Baldy dune in Indiana, USA[J]. Aeolian Research, 2016, 23: 103-114.

    [6] 何京麗,郭建英,邢恩德,等. 黃河烏蘭布和沙漠段沿岸風(fēng)沙流結(jié)構(gòu)與沙丘移動規(guī)律[J]. 農(nóng)業(yè)工程學(xué)報,2012,28(17):71-77.

    He Jingli, Guo Jianying, Xing Ende, et al. Structure of wind-sand flow and law of dune movement along bank of Yellow River in Ulan Buh desert[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2012, 28(17): 71-77. (in Chinese with English abstract)

    [7] Yang Z L, Qian G Q, Dong Z B, et al. Migration of barchan dunes and factors that influence migration in the Sanlongsha dune field of the northern Kumtagh Sand Sea, China[J]. Geomorphology, 2021, 378: 107615.

    [8] Yang J H, Dong Z B, Liu Z Y, et al. Migration of barchan dunes in the western Quruq Desert, northwestern China[J]. Earth Surface Processes and Landforms, 2019, 44(10): 2016-2029.

    [9] Ali E, Xu W B, Ding X L. Improved optical image matching time series inversion approach for monitoring dune migration in North Sinai Sand Sea: Algorithm procedure, application, and validation[J]. ISPRS Journal of Photogrammetry and Remote Sensing, 2020, 164: 106-124.

    [10] Dong Z B, Wang X M, Chen G T. Monitoring sand dune advance in the Taklimakan Desert[J]. Geomorphology, 2000, 35(3/4): 219-231.

    [11] Buddhika M, James L, Michael M, et al. Factors influencing movement of the Manila Dunes and its impact on establishing non-native species[J]. Remote Sensing, 2020, 12(10): 1536.

    [12] 王靜璞,劉連友,沈玲玲. 基于Google Earth的毛烏素沙地新月形沙丘移動規(guī)律研究[J]. 遙感技術(shù)與應(yīng)用,2013,28(6):1094-1100.

    Wang Jingpu, Liu Lianyou, Shen Lingling. Research of the barchan dunes movement in the Mu Us Sandy Land on Google Earth software[J]. Remote Sensing Technology and Application, 2013, 28(6): 1094-1100. (in Chinese with English abstract)

    [13] 劉佳,王利民,滕飛,等. Google Earth影像輔助的農(nóng)作物面積地面樣方調(diào)查[J]. 農(nóng)業(yè)工程學(xué)報,2015,31(24):149-154.

    Liu Jia, Wang Limin, Teng Fei, et al. Crop area ground sample survey using Google Earth image-aided[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(24): 149-154. (in Chinese with English abstract)

    [14] Ding C, Zhang L, Liao M S, et al. Quantifying the spatio-temporal patterns of dune migration near Minqin Oasis in northwestern China with time series of Landsat-8 and Sentinel-2 observations[J]. Remote Sensing of Environment, 2020, 236: 111498.

    [15] 董治寶,呂萍. 70年來中國風(fēng)沙地貌學(xué)的發(fā)展[J]. 地理學(xué)報,2020,75(3):509-528.

    Dong Zhibao, Lyu Ping. Development of aeolian geomorphology in China in the past 70 years[J]. Acta Geographica Sinica, 2020, 75(3): 509-528. (in Chinese with English abstract)

    [16] Li C, Dong Z B, Cui X J, et al. Pattern analysis of simple transverse dunes in China’s Qaidam Basin, north of the Kunlun Mountains[J]. Environmental Earth Sciences, 2016, 75(9): 729.

    [17] 陳宗顏,董治寶,汪青春,等. 柴達(dá)木盆地風(fēng)況及輸沙勢特征[J]. 中國沙漠,2020,40(1):195-203.

    Chen Zongyan, Dong Zhibao, Wang Qingchun, et al. Characteristics of wind regime and sand drift potential in Qaidam Basin of China[J]. Journal of Desert Research, 2020, 40(1): 195-203. (in Chinese with English abstract)

    [18] 朱震達(dá),陳治平,吳正,等. 塔克拉瑪干沙漠風(fēng)沙地貌研究[M]. 北京:科學(xué)出版社,1981:78-80.

    [19] 李鎮(zhèn),張巖,楊松,等. QuickBird影像目視解譯法提取切溝形態(tài)參數(shù)的精度分析[J]. 農(nóng)業(yè)工程學(xué)報,2014,30(20):179-186.

    Li Zhen, Zhang Yan, Yang Song, et al. Error assessment of extracting morphological parameters of bank gullies by manual visual interpretation based on QuickBird imagery[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2014, 30(20): 179-186. (in Chinese with English abstract)

    [20] Hamdan M A, Refaat A A, Abdel W M. Morphologic characteristics and migration rate assessment of barchan dunes in the Southeastern Western Desert of Egypt[J]. Geomorphology, 2016, 257: 57-74.

    [21] 鮑鋒,董治寶,張正偲. 柴達(dá)木盆地風(fēng)沙地貌區(qū)風(fēng)況特征[J]. 中國沙漠,2015,35(3):549-554.

    Bao Feng, Dong Zhibao, Zhang Zhengsi. Wind regime in the Qaidam Basin Desert[J]. Journal of Desert Research, 2015, 35(3): 549-554. (in Chinese with English abstract)

    [22] Yang X P, Rost K T, Lehmkuhl F, et al. The evolution of dry lands in northern China and in the Republic of Mongolia since the Last Glacial Maximum[J]. Quaternary International, 2004, 118/119: 69-85.

    [23] 李小妹,嚴(yán)平,錢瑤,等. 不同氣候帶的河道與沙丘分布格局及其類型劃分[J]. 中國沙漠,2017,37(5):821-829.

    Li Xiaomei, Yan Ping, Qian Yao, et al. Distribution patterns and its classification at the scale of channel-sand dune in different climate zones[J]. Journal of Desert Research, 2017, 37(5): 821-829. (in Chinese with English abstract)

    [24] 韓有香,劉彩紅,李國山. 青海高原典型生態(tài)區(qū)域未來氣候變化趨勢預(yù)估[J]. 草業(yè)科學(xué),2019,36(6):1518-1530.

    Han Youxiang, Liu Caihong, Li Guoshan. Forecast of future climate change trend in typical ecological regions of the Qinghai Plateau[J]. Pratacultural Science, 2019, 36(6): 1518-1530. (in Chinese with English abstract)

    [25] Tsoar H, Parteli Eric J R. Bidirectional winds, barchan dune asymmetry and formation of seif dunes from barchans: A discussion[J]. Environmental Earth Sciences, 2016, 75(18): 1237.

    Spatial regionalization and response to morphological parameters of dune migration in the Qaidam Basin of China

    Li Jianjun1, Jiao Juying1,2※, Cao Xue2, Bai Leichao1, Chen Tongde1, Yan Xiqin2, Qi Hongkun1

    (1. State Key Laboratory of Soil Erosion and Dryland Farming on the Loess Plateau, Institute of Soil and Water Conservation, Northwest A&F University, Yangling 712100, China; 2. Institute of Soil and Water Conservation, Chinese Academy of Sciences and Ministry of Water Resources, Yangling 712100, China)

    Desertification and dune migration are two serious forms of land degradation, resulting in the destruction of grassland and farmland ecosystems. Sandstorm activities occur more frequently in the Tibetan Plateau of western China in recent years, particularly with the expansion of desertification. The Qaidam Basin is the largest intermountain basin located in the northwestern portion of the Tibetan Plateau. This basin is usually considered to be 'the region of Earth most similar to Mars', providing for a sound reference for the study of Mars aeolian geomorphology. However, it is still lacking the understanding of dune migration in the Qaidam Basin. Taking the transverse dunes (such as barchan dune) as the research objects, this study aims to clarify the spatial regionalization and response to morphological parameters of dune migration in the Qaidam Basin. Geo Eye, Quick bird and World View satellite images were utilized to measure morphological parameters and migration using the Google Earth software. The main results were listed as follows: 1) The average migration rate of dunes was between 0 and 23.53 m/a in the study area, with an average of 4.66 m/a, more than half of which was the medium speed (53.73%). The migration direction of the dune was between 53.52° and 150.06°, with an average of 104.13°, indicating a better agreement with the main wind direction. The average speed of dune migration slowed down spatially from northwest to southeast, but there was no obvious trend. Specifically, the fastest migration was found in the middle (5.22 m/a), followed by the northwest (4.84 m/a), the north (3.49 m/a), the southwest (3.35 m/a), and the southeast (3.27 m/a). The migration direction changed with the direction of the mountain range. 2)The migration of the dune was much slower as the volume of dune increased in the same environment. There was a very significant (<0.01) exponential negative correlation between the migration speeds and morphological parameters of dune, such as the length of windward slope (2=0.716), length of lee slope (2=0.705), height (2=0.705), width (2=0.988), perimeter (2=0.986), basal area (2=0.955), whereas, a very significant (<0.01) linear negative correlation with the coefficient of dune shape (2=0.438), and a significant (<0.05) positive linear correlation with the length-width ratio (2=0.130). Therefore, more attention can be paid to the dune width in future studies using remote sensing images. 3)The dune migration was much slower as the wind velocity decreased and the precipitation increased in the whole basin. Furthermore, the dune moved slowly under the great vegetation coverage in the southeastern basin. 4)Dense rivers blocked the dune migration and protected Golmud City, but the dune migration still posed great threats to the whole ecosystem. This finding can make a great contribution to the accurate understanding of sandstorm activities for regional town planning and protection in the national ecological security shelter zone on the Tibetan Plateau of western China.

    remote sensing; dune; migration; morphological parameters; Qaidam Basin

    2020-11-25

    2021-01-17

    第二次青藏高原綜合考察研究(2019QZKK0603);中國科學(xué)院戰(zhàn)略性先導(dǎo)科技專項“泛第三極環(huán)境變化與綠色絲綢之路建設(shè)”子課題“土壤侵蝕定量評價與分區(qū)防控對策”(XDA20040202)

    李建軍,博士生,主要研究方向為土壤侵蝕。Email:lijianjun@nwafu.edu.cn

    焦菊英,博士,研究員。主要研究方向為流域侵蝕產(chǎn)沙、土壤侵蝕與植被關(guān)系及水土保持效益評價。Email:jyjiao@ms.iswc.ac.cn

    10.11975/j.issn.1002-6819.2021.07.038

    P931.3

    A

    1002-6819(2021)-07-0309-06

    李建軍,焦菊英,曹雪,等. 柴達(dá)木盆地沙丘移動的空間分異及對形態(tài)參數(shù)的響應(yīng)[J]. 農(nóng)業(yè)工程學(xué)報,2021,37(7):309-314. doi:10.11975/j.issn.1002-6819.2021.07.038 http://www.tcsae.org

    Li Jianjun, Jiao Juying, Cao Xue, et al. Spatial regionalization and response to morphological parameters of dune migration in the Qaidam Basin of China[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2021, 37(7): 309-314. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2021.07.038 http://www.tcsae.org

    猜你喜歡
    柴達(dá)木盆地沙丘速度
    出乎意料
    行駛速度
    速度
    柴達(dá)木盆地北緣鋰多金屬礦成礦條件及找礦潛力
    沙丘
    柴達(dá)木盆地西部古近系石膏及其硫同位素分布特征
    沙丘
    比速度更速度——“光腦”來了
    柴達(dá)木盆地柴北緣侏羅系烴源巖特征分析
    基于新鉆井對柴達(dá)木盆地東部石炭系分布及控制因素的再認(rèn)識
    日韩视频在线欧美| 亚洲国产精品成人综合色| 久久99蜜桃精品久久| 国产高清不卡午夜福利| 人人妻人人看人人澡| 欧美潮喷喷水| 伊人久久国产一区二区| 一级黄片播放器| 欧美最新免费一区二区三区| 免费在线观看成人毛片| 综合色av麻豆| 高清日韩中文字幕在线| 中文欧美无线码| 草草在线视频免费看| 亚洲成人中文字幕在线播放| 天堂中文最新版在线下载 | 亚洲国产精品成人综合色| 又粗又硬又长又爽又黄的视频| 三级国产精品欧美在线观看| av一本久久久久| 午夜福利高清视频| 欧美激情久久久久久爽电影| av线在线观看网站| 秋霞在线观看毛片| 亚洲四区av| 最近手机中文字幕大全| 最后的刺客免费高清国语| 国产视频首页在线观看| 久久久精品94久久精品| 最近最新中文字幕免费大全7| 欧美日韩在线观看h| 亚洲欧美精品自产自拍| 一个人看视频在线观看www免费| 日本欧美国产在线视频| 大片电影免费在线观看免费| 我的老师免费观看完整版| 久热这里只有精品99| 国产亚洲一区二区精品| 亚洲欧美一区二区三区国产| 国产亚洲精品久久久com| 好男人视频免费观看在线| 国产精品精品国产色婷婷| 久久精品久久久久久噜噜老黄| 国产高潮美女av| 免费观看在线日韩| 国产精品.久久久| 搞女人的毛片| 国产精品人妻久久久久久| 啦啦啦啦在线视频资源| 热99国产精品久久久久久7| 久久久久久久午夜电影| 国产成人精品婷婷| 日韩三级伦理在线观看| 国产精品三级大全| 久久精品综合一区二区三区| 男女国产视频网站| 爱豆传媒免费全集在线观看| 国产成人免费观看mmmm| 亚洲精品色激情综合| 国产黄a三级三级三级人| 亚洲成人中文字幕在线播放| 汤姆久久久久久久影院中文字幕| 亚洲av欧美aⅴ国产| 欧美日韩一区二区视频在线观看视频在线 | 大片电影免费在线观看免费| 亚洲精品久久久久久婷婷小说| 成年版毛片免费区| 久久久精品免费免费高清| tube8黄色片| 一级黄片播放器| 成人漫画全彩无遮挡| 久久精品久久精品一区二区三区| 久久人人爽av亚洲精品天堂 | 男女啪啪激烈高潮av片| av.在线天堂| 看十八女毛片水多多多| 九草在线视频观看| 亚洲精品一二三| 日日摸夜夜添夜夜添av毛片| 夫妻性生交免费视频一级片| 女人被狂操c到高潮| 婷婷色综合www| 亚洲va在线va天堂va国产| 91精品伊人久久大香线蕉| 免费观看无遮挡的男女| 91在线精品国自产拍蜜月| 色播亚洲综合网| 国产亚洲最大av| av国产精品久久久久影院| 国产成人福利小说| 亚洲av.av天堂| 亚洲精品日韩av片在线观看| 亚洲精品,欧美精品| 草草在线视频免费看| 亚洲综合色惰| 国产一区有黄有色的免费视频| 国产精品国产三级国产av玫瑰| 亚洲精品乱久久久久久| 99热网站在线观看| 精品午夜福利在线看| 国产男人的电影天堂91| 国产黄片视频在线免费观看| 舔av片在线| 嫩草影院新地址| 日本爱情动作片www.在线观看| 少妇猛男粗大的猛烈进出视频 | 亚洲欧美一区二区三区国产| 国产精品嫩草影院av在线观看| 建设人人有责人人尽责人人享有的 | 人妻一区二区av| 国产成人一区二区在线| 国产午夜精品一二区理论片| 午夜免费男女啪啪视频观看| 午夜视频国产福利| 91狼人影院| 精品国产露脸久久av麻豆| 熟女电影av网| 天天一区二区日本电影三级| 搞女人的毛片| 黄色配什么色好看| 国产成人精品福利久久| 久久久久国产网址| 波野结衣二区三区在线| 身体一侧抽搐| 亚洲av中文av极速乱| 成人高潮视频无遮挡免费网站| 亚洲av中文字字幕乱码综合| 人妻 亚洲 视频| 天天一区二区日本电影三级| 免费大片黄手机在线观看| 夜夜看夜夜爽夜夜摸| 麻豆久久精品国产亚洲av| 国产极品天堂在线| 黄片wwwwww| 色哟哟·www| a级毛片免费高清观看在线播放| 国产国拍精品亚洲av在线观看| 国产亚洲午夜精品一区二区久久 | 晚上一个人看的免费电影| 久久久久久久大尺度免费视频| 久久久久久久大尺度免费视频| 国产欧美日韩精品一区二区| 黄色日韩在线| 一个人看的www免费观看视频| 午夜日本视频在线| 亚洲精品成人久久久久久| 精品久久久久久电影网| 高清av免费在线| 欧美一区二区亚洲| 国产精品熟女久久久久浪| 香蕉精品网在线| 精品视频人人做人人爽| 国产伦精品一区二区三区四那| 国产 一区 欧美 日韩| 日本三级黄在线观看| 日本黄色片子视频| 亚洲综合精品二区| 国产av国产精品国产| 男人爽女人下面视频在线观看| 亚洲av在线观看美女高潮| 国产又色又爽无遮挡免| 国产精品.久久久| 看非洲黑人一级黄片| 岛国毛片在线播放| av国产精品久久久久影院| 老女人水多毛片| 九九久久精品国产亚洲av麻豆| 国产白丝娇喘喷水9色精品| 国产成人午夜福利电影在线观看| 老师上课跳d突然被开到最大视频| 别揉我奶头 嗯啊视频| 又黄又爽又刺激的免费视频.| 蜜桃久久精品国产亚洲av| 尤物成人国产欧美一区二区三区| 高清在线视频一区二区三区| av在线app专区| 深爱激情五月婷婷| 亚洲国产成人一精品久久久| 插逼视频在线观看| 久久久精品欧美日韩精品| 亚洲av不卡在线观看| 欧美区成人在线视频| 久久久国产一区二区| 午夜爱爱视频在线播放| h日本视频在线播放| 日韩不卡一区二区三区视频在线| 简卡轻食公司| 成人亚洲精品一区在线观看 | 精品久久久久久久末码| 精品一区二区三卡| 在线观看一区二区三区| 51国产日韩欧美| 黄色配什么色好看| 亚洲欧洲日产国产| 尤物成人国产欧美一区二区三区| 亚洲av国产av综合av卡| 一区二区三区免费毛片| 91狼人影院| 国产成人91sexporn| 麻豆成人午夜福利视频| 欧美激情久久久久久爽电影| 欧美少妇被猛烈插入视频| 国产一级毛片在线| 22中文网久久字幕| 欧美成人a在线观看| 免费少妇av软件| 久久99精品国语久久久| 热99国产精品久久久久久7| 久久亚洲国产成人精品v| 亚洲av成人精品一二三区| 麻豆成人午夜福利视频| 波野结衣二区三区在线| 一级毛片久久久久久久久女| 男人添女人高潮全过程视频| 亚洲欧美成人综合另类久久久| 中文乱码字字幕精品一区二区三区| 日日摸夜夜添夜夜添av毛片| 国产成年人精品一区二区| 成年免费大片在线观看| 国产女主播在线喷水免费视频网站| 综合色丁香网| 日韩欧美一区视频在线观看 | 亚洲一区二区三区欧美精品 | 午夜福利网站1000一区二区三区| av在线播放精品| 联通29元200g的流量卡| 日韩制服骚丝袜av| 热re99久久精品国产66热6| 国产熟女欧美一区二区| 国产高清有码在线观看视频| 80岁老熟妇乱子伦牲交| 精品一区二区免费观看| 亚洲精品中文字幕在线视频 | 国产免费又黄又爽又色| 免费播放大片免费观看视频在线观看| 久久精品夜色国产| 亚洲av日韩在线播放| 国产高清有码在线观看视频| 91精品一卡2卡3卡4卡| 99热网站在线观看| 女的被弄到高潮叫床怎么办| 亚洲av国产av综合av卡| 国产片特级美女逼逼视频| 免费大片黄手机在线观看| 国产av不卡久久| 成人亚洲精品av一区二区| 亚洲电影在线观看av| 久久综合国产亚洲精品| 国产精品国产av在线观看| 国产精品秋霞免费鲁丝片| 水蜜桃什么品种好| 日本黄大片高清| 校园人妻丝袜中文字幕| 赤兔流量卡办理| 精品国产乱码久久久久久小说| 亚洲国产精品成人久久小说| 大陆偷拍与自拍| 国产精品精品国产色婷婷| 国产综合懂色| 99热国产这里只有精品6| 日韩 亚洲 欧美在线| 夫妻性生交免费视频一级片| 久久精品久久精品一区二区三区| 国产午夜精品一二区理论片| 777米奇影视久久| 亚洲精品国产av蜜桃| 天天躁夜夜躁狠狠久久av| 欧美xxxx黑人xx丫x性爽| 欧美老熟妇乱子伦牲交| 国产国拍精品亚洲av在线观看| 中文字幕制服av| 久久久亚洲精品成人影院| 能在线免费看毛片的网站| 国产毛片在线视频| 亚洲最大成人av| 国产精品麻豆人妻色哟哟久久| 少妇人妻一区二区三区视频| 免费黄色在线免费观看| 两个人的视频大全免费| 亚洲欧美清纯卡通| av福利片在线观看| 一个人观看的视频www高清免费观看| 成人午夜精彩视频在线观看| 亚洲av福利一区| 国产久久久一区二区三区| 国产黄色免费在线视频| 亚洲欧美中文字幕日韩二区| 永久网站在线| 最近手机中文字幕大全| 大香蕉久久网| 一级二级三级毛片免费看| 精品久久久噜噜| 日韩精品有码人妻一区| 国产精品蜜桃在线观看| 日本与韩国留学比较| 国产老妇伦熟女老妇高清| 黄色日韩在线| 欧美激情在线99| 97超碰精品成人国产| 国产一区亚洲一区在线观看| 国产 一区 欧美 日韩| 婷婷色综合大香蕉| 欧美日韩亚洲高清精品| av播播在线观看一区| 水蜜桃什么品种好| 午夜福利高清视频| 色婷婷久久久亚洲欧美| 国产精品久久久久久久电影| 欧美日韩国产mv在线观看视频 | 日韩伦理黄色片| 久久人人爽人人片av| 大话2 男鬼变身卡| 亚洲av免费高清在线观看| 亚洲成色77777| 啦啦啦中文免费视频观看日本| 91精品一卡2卡3卡4卡| 2018国产大陆天天弄谢| 中文字幕人妻熟人妻熟丝袜美| 亚洲熟女精品中文字幕| 乱系列少妇在线播放| 欧美人与善性xxx| 啦啦啦在线观看免费高清www| 成人二区视频| av在线蜜桃| 丝袜美腿在线中文| 大香蕉久久网| 一个人观看的视频www高清免费观看| 亚洲,欧美,日韩| 国产亚洲最大av| 日韩电影二区| 高清av免费在线| 国产日韩欧美在线精品| 亚洲自偷自拍三级| 成人综合一区亚洲| 又粗又硬又长又爽又黄的视频| 日日啪夜夜爽| 国产成人免费观看mmmm| 国产伦理片在线播放av一区| 国产视频首页在线观看| 赤兔流量卡办理| 精品久久国产蜜桃| 国产伦精品一区二区三区四那| 777米奇影视久久| 男女那种视频在线观看| 国产精品国产三级国产av玫瑰| 九九爱精品视频在线观看| 97人妻精品一区二区三区麻豆| 三级国产精品片| 91午夜精品亚洲一区二区三区| 男人爽女人下面视频在线观看| 欧美日韩一区二区视频在线观看视频在线 | 久久ye,这里只有精品| 精品视频人人做人人爽| 久久精品综合一区二区三区| 波多野结衣巨乳人妻| 国产精品三级大全| 久久久欧美国产精品| 色哟哟·www| 精品视频人人做人人爽| 欧美少妇被猛烈插入视频| 日韩一区二区视频免费看| 一级片'在线观看视频| 真实男女啪啪啪动态图| 久久人人爽人人爽人人片va| 在线观看一区二区三区激情| 欧美精品国产亚洲| 国产探花在线观看一区二区| 久久久久久国产a免费观看| 联通29元200g的流量卡| 成人特级av手机在线观看| 日韩视频在线欧美| 色哟哟·www| 人妻 亚洲 视频| 国产黄a三级三级三级人| 国产免费又黄又爽又色| 亚洲国产欧美人成| 国产精品国产av在线观看| 91久久精品国产一区二区三区| 人妻 亚洲 视频| 91久久精品国产一区二区三区| 免费看不卡的av| 狂野欧美激情性xxxx在线观看| 大又大粗又爽又黄少妇毛片口| 美女国产视频在线观看| 国产色爽女视频免费观看| 亚洲精品色激情综合| 国产色爽女视频免费观看| 黑人高潮一二区| 亚洲在线观看片| 啦啦啦在线观看免费高清www| 青春草视频在线免费观看| 在线观看三级黄色| 国产伦精品一区二区三区四那| 久久久久久久午夜电影| 尾随美女入室| 日本爱情动作片www.在线观看| 免费大片黄手机在线观看| 国产免费视频播放在线视频| 夫妻性生交免费视频一级片| 久久国内精品自在自线图片| 永久免费av网站大全| 下体分泌物呈黄色| 国产精品国产三级专区第一集| 韩国高清视频一区二区三区| 五月伊人婷婷丁香| 九九爱精品视频在线观看| 久久久久久伊人网av| 国产色爽女视频免费观看| 欧美精品国产亚洲| 97超碰精品成人国产| 综合色丁香网| .国产精品久久| 亚洲综合精品二区| 国产女主播在线喷水免费视频网站| 黄色视频在线播放观看不卡| 亚洲精品久久久久久婷婷小说| 亚洲精品成人av观看孕妇| 免费高清在线观看视频在线观看| 欧美日韩亚洲高清精品| 国产成人免费观看mmmm| 成人一区二区视频在线观看| 在线观看国产h片| 麻豆成人av视频| 亚洲色图综合在线观看| 国产精品久久久久久久电影| 黄色视频在线播放观看不卡| 日韩欧美精品v在线| 熟妇人妻不卡中文字幕| 国产一区亚洲一区在线观看| 久久久欧美国产精品| 久久久久久久精品精品| 三级男女做爰猛烈吃奶摸视频| 亚洲欧美清纯卡通| 亚洲精品乱码久久久久久按摩| 成年免费大片在线观看| 亚洲怡红院男人天堂| 亚洲人成网站高清观看| 精品一区在线观看国产| 不卡视频在线观看欧美| 六月丁香七月| 国产精品一区二区性色av| 亚洲精品成人av观看孕妇| 久久国产乱子免费精品| 又爽又黄a免费视频| 久久人人爽人人爽人人片va| 色5月婷婷丁香| 亚洲欧美日韩卡通动漫| 99视频精品全部免费 在线| 热99国产精品久久久久久7| 亚洲精品自拍成人| 美女视频免费永久观看网站| 99久久中文字幕三级久久日本| 成年女人在线观看亚洲视频 | 久久女婷五月综合色啪小说 | 三级男女做爰猛烈吃奶摸视频| 国产精品熟女久久久久浪| 亚洲一区二区三区欧美精品 | av又黄又爽大尺度在线免费看| 午夜免费观看性视频| 久久人人爽av亚洲精品天堂 | 国产在线男女| 在线观看国产h片| 99久久九九国产精品国产免费| 王馨瑶露胸无遮挡在线观看| 交换朋友夫妻互换小说| 亚洲精品456在线播放app| 免费av毛片视频| 精品视频人人做人人爽| 精品少妇久久久久久888优播| .国产精品久久| 国产精品久久久久久精品电影小说 | 2022亚洲国产成人精品| 中文字幕制服av| 国产成年人精品一区二区| 亚洲精华国产精华液的使用体验| av卡一久久| 午夜福利在线观看免费完整高清在| 国产黄色视频一区二区在线观看| 26uuu在线亚洲综合色| 亚洲天堂国产精品一区在线| 久久久久久久午夜电影| 女人久久www免费人成看片| 大香蕉97超碰在线| 国产极品天堂在线| 熟妇人妻不卡中文字幕| 国产男女内射视频| 亚洲欧美日韩东京热| 春色校园在线视频观看| 国产精品国产三级专区第一集| 蜜桃久久精品国产亚洲av| 国产亚洲91精品色在线| 国产伦在线观看视频一区| 大片免费播放器 马上看| 人妻夜夜爽99麻豆av| 色视频www国产| 国产精品av视频在线免费观看| 在线观看一区二区三区激情| 国产综合精华液| 国产老妇女一区| 精品人妻熟女av久视频| 婷婷色综合大香蕉| 日韩电影二区| 女人十人毛片免费观看3o分钟| 日本-黄色视频高清免费观看| 少妇熟女欧美另类| 国产精品秋霞免费鲁丝片| 99久久精品一区二区三区| 国产av国产精品国产| 亚洲欧美一区二区三区黑人 | 亚洲精品第二区| 黄片wwwwww| 国产色婷婷99| 天美传媒精品一区二区| 欧美日韩一区二区视频在线观看视频在线 | 99久久精品一区二区三区| 亚洲欧美日韩卡通动漫| 亚洲精品aⅴ在线观看| 国产精品久久久久久精品古装| 18禁裸乳无遮挡动漫免费视频 | 久久午夜福利片| 男女啪啪激烈高潮av片| 亚洲不卡免费看| 亚洲内射少妇av| 亚洲最大成人中文| 在现免费观看毛片| 国产 一区精品| 在线观看人妻少妇| 在线观看国产h片| 精品国产露脸久久av麻豆| 欧美性猛交╳xxx乱大交人| 久久久久久久久久人人人人人人| 狂野欧美激情性bbbbbb| 亚洲美女视频黄频| 小蜜桃在线观看免费完整版高清| 在线观看免费高清a一片| 亚洲欧美精品自产自拍| 免费看av在线观看网站| 在线观看三级黄色| 九草在线视频观看| 乱系列少妇在线播放| 亚洲欧洲国产日韩| 人人妻人人看人人澡| 97热精品久久久久久| 亚洲一区二区三区欧美精品 | 男女国产视频网站| 最近中文字幕2019免费版| 国产中年淑女户外野战色| 国产伦精品一区二区三区视频9| 成人欧美大片| 亚洲高清免费不卡视频| 亚洲真实伦在线观看| av在线老鸭窝| 久久精品久久精品一区二区三区| 日韩欧美 国产精品| 18禁动态无遮挡网站| 一级a做视频免费观看| 99久久九九国产精品国产免费| 欧美xxxx性猛交bbbb| 精品久久久久久电影网| 免费av观看视频| 男插女下体视频免费在线播放| 中文欧美无线码| 国产色婷婷99| 国产精品久久久久久av不卡| 久久久午夜欧美精品| tube8黄色片| 欧美日韩视频高清一区二区三区二| 日本-黄色视频高清免费观看| 香蕉精品网在线| 国产亚洲午夜精品一区二区久久 | 美女xxoo啪啪120秒动态图| 十八禁网站网址无遮挡 | 在线观看免费高清a一片| 亚洲,欧美,日韩| 日韩成人av中文字幕在线观看| av天堂中文字幕网| 亚洲国产色片| 国产精品成人在线| 亚洲久久久久久中文字幕| 久久久久精品久久久久真实原创| 中文字幕人妻熟人妻熟丝袜美| 国产成人福利小说| 欧美激情国产日韩精品一区| 91久久精品国产一区二区成人| 一个人看的www免费观看视频| 亚洲av免费高清在线观看| 亚洲精品乱码久久久久久按摩| 国产精品久久久久久精品古装| 91久久精品国产一区二区三区| 人人妻人人爽人人添夜夜欢视频 | 日韩在线高清观看一区二区三区| 91在线精品国自产拍蜜月| 久久鲁丝午夜福利片| 最近中文字幕高清免费大全6| 色综合色国产| 亚洲国产最新在线播放| 婷婷色麻豆天堂久久| 九九爱精品视频在线观看| 有码 亚洲区| 免费大片黄手机在线观看| 国产片特级美女逼逼视频| 国产精品偷伦视频观看了| videos熟女内射| 久久久久久国产a免费观看| 岛国毛片在线播放| 一区二区三区乱码不卡18| 精品国产一区二区三区久久久樱花 | 神马国产精品三级电影在线观看| 亚洲精品456在线播放app| 免费av毛片视频| 久热这里只有精品99| 一级二级三级毛片免费看| 人人妻人人澡人人爽人人夜夜| 成年人午夜在线观看视频| 街头女战士在线观看网站|