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

    基于ESAI的黃土高原荒漠化風(fēng)險評估

    2017-06-05 15:01:13張建香張多勇劉萬鋒胡愛萍胡劍橋
    水土保持通報 2017年2期
    關(guān)鍵詞:敏感區(qū)黃土高原荒漠化

    張建香, 張多勇, 劉萬鋒, 胡愛萍, 王 東, 胡劍橋

    (1.隴東學(xué)院 土木工程學(xué)院, 甘肅 慶陽 745000; 2.甘肅省高校黃土的工程性質(zhì)及工程應(yīng)用省級重點(diǎn)實(shí)驗室,甘肅慶陽 745000; 3.慶陽市荒漠化防治研究中心, 甘肅 慶陽 745000)

    基于ESAI的黃土高原荒漠化風(fēng)險評估

    張建香1,2, 張多勇3, 劉萬鋒1,2, 胡愛萍1, 王 東1, 胡劍橋1

    (1.隴東學(xué)院 土木工程學(xué)院, 甘肅 慶陽 745000; 2.甘肅省高校黃土的工程性質(zhì)及工程應(yīng)用省級重點(diǎn)實(shí)驗室,甘肅慶陽 745000; 3.慶陽市荒漠化防治研究中心, 甘肅 慶陽 745000)

    [目的] 對黃土高原荒漠化風(fēng)險進(jìn)行評估,為黃土高原生態(tài)環(huán)境的恢復(fù)、建設(shè)和保護(hù)提供科學(xué)依據(jù)。 [方法] 通過收集黃土高原的地形、氣候、植被、土壤以及社會經(jīng)濟(jì)等方面的數(shù)據(jù),借助RS和GIS平臺,實(shí)現(xiàn)黃土高原荒漠化風(fēng)險評估的空間化和數(shù)字化,在此基礎(chǔ)上分析荒漠化的成因,構(gòu)建荒漠化風(fēng)險評價指標(biāo)體系,建立基于環(huán)境敏感性區(qū)指標(biāo)(ESAI)的荒漠化風(fēng)險評估模型,分析荒漠化風(fēng)險程度的空間格局,探索黃土高原不同區(qū)域荒漠化形成的主要原因。 [結(jié)果] (1) 基于土壤、氣候、植被3種要素的環(huán)境敏感區(qū)生物物理指標(biāo)顯示:黃土高原大約1/4的區(qū)域(25.2%)為高風(fēng)險區(qū),屬于嚴(yán)重荒漠化,幾乎2/3的區(qū)域(62.8%)是輕微荒漠化,11.5%的地區(qū)為潛在荒漠化,只有0.5%的地區(qū)無荒漠化現(xiàn)象; (2) 加入人類誘發(fā)因素后,改變了黃土高原荒漠化風(fēng)險區(qū)的原有格局。其中,極低、低度和極高度敏感區(qū)減少了5.6%,1.1%和3.8%;與此同時,較低和較高度敏感區(qū)增加了4.4%和4.5%。 [結(jié)論] (1) 該模型能很好地說明黃土高原荒漠化風(fēng)險的空間分布格局,其荒漠化程度由西北向東南地區(qū)逐漸減弱; (2) 人類活動已經(jīng)在一定程度上打破了長期穩(wěn)定的自然生態(tài)系統(tǒng),并且縮小了不同程度荒漠化之間的差距。

    荒漠化; 環(huán)境敏感區(qū)指標(biāo)(ESAI); 地理信息系統(tǒng)(GIS); 黃土高原

    文獻(xiàn)參數(shù): 張建香, 張多勇, 劉萬鋒, 等.基于ESAI的黃土高原荒漠化風(fēng)險評估[J].水土保持通報,2017,37(2):339-344.DOI:10.13961/j.cnki.stbctb.2017.02.051; Zhang Jianxiang, Zhang Duoyong, Liu Wanfeng, et al. ESAI Based Assessment of Desertification Risk in Loess Plateau[J]. Bulletin of Soil and Water Conservation, 2017,37(2):339-344.DOI:10.13961/j.cnki.stbctb.2017.02.051

    1977年聯(lián)合國荒漠化會議召開以后,“荒漠化”作為一個生態(tài)環(huán)境問題被世界所公認(rèn)。1991年9月,聯(lián)合國環(huán)境署在日內(nèi)瓦召開了荒漠化防治第8次顧問會議,對世界荒漠化防治作了全面研究。近年來,學(xué)術(shù)界主要傾向于荒漠化的演變過程、荒漠化驅(qū)動機(jī)制、荒漠化監(jiān)測與評估及其指標(biāo)體系的建立等方面的研究,地理信息系統(tǒng)技術(shù)在荒漠化研究中被廣泛應(yīng)用[1]。盡管如此,目前研究尚存在諸多問題,荒漠化風(fēng)險評估方法尚待進(jìn)一步深化。第一,在研究深度上,應(yīng)將傳統(tǒng)的小尺度、靜態(tài)的定性調(diào)查描述轉(zhuǎn)向定位、定量地對荒漠化格局、過程及內(nèi)在機(jī)制進(jìn)行宏觀、動態(tài)及具有輔助決策作用的系統(tǒng)研究;第二,在研究方式上,應(yīng)更加注重學(xué)科的聯(lián)合和滲透,同時應(yīng)結(jié)合野外定點(diǎn)觀測、遙感監(jiān)測和GPS/GIS應(yīng)用,實(shí)現(xiàn)荒漠化的動態(tài)自動監(jiān)測[2];第三,為研究土地退化程度,Kosmas教授[3]于1999年提出了環(huán)境敏感區(qū)法(environmentally sensitive areas, ESAs),該方法考慮了土壤、氣候、植被、土地管理等質(zhì)量指標(biāo),但隨著人類活動對生態(tài)環(huán)境造成的干擾日益加大,ESAs方法中的各項指標(biāo)已不能充分反映環(huán)境質(zhì)量的高低。為此,本研究以中國典型的生態(tài)環(huán)境脆弱區(qū)黃土高原為研究對象,在ESAs方法的基礎(chǔ)上,綜合地形、氣候、植被、土壤、人口、經(jīng)濟(jì)等多種生物—物理方面的環(huán)境質(zhì)量影響因素,建立基于環(huán)境敏感性評估指標(biāo)的黃土高原區(qū)荒漠化風(fēng)險評估模型ESAI(environmental sensitive areas index)[4-5],研究黃土高原不同程度荒漠化的空間分布格局,探索不同區(qū)域影響荒漠化的主要原因,進(jìn)而為黃土高原生態(tài)環(huán)境的恢復(fù)、建設(shè)和保護(hù)提供科學(xué)的依據(jù)。

    1 材料與方法

    1.1 研究區(qū)概況

    黃土高原坐落在中國大陸的中北部,位于黃河流域的中游地帶(33°41′—41°16′N,100°52′—114°33′E),跨陜、甘、青、晉、寧、豫、內(nèi)蒙古7個省(區(qū)),東起太行山,西至烏鞘嶺—日月山,南達(dá)秦嶺,北至陰山,東西長約1 300 km,南北寬約800 km,總面積約6.23×105km2,平均海拔1 500~2 000 m[6-7]。黃土高原具有典型的大陸性季風(fēng)氣候,冬季寒冷干燥,夏季炎熱潮濕。年平均溫度由西北(4.3 ℃)向東南(14.3 ℃)逐漸遞增。多年平均降水量約為460 mm,且降水量的年際變化較大[8]。黃土高原被深厚的黃土層覆蓋,大部分地區(qū)黃土厚度達(dá)100 m,已成世界黃土層最厚、分布最廣泛、最典型的黃土地貌。同時,由于該地區(qū)植被稀少,而黃土土質(zhì)疏松,且暴雨多集中在夏秋兩季,地形起伏大,極易出現(xiàn)水土流失,導(dǎo)致黃土高原成為全球生態(tài)退化最嚴(yán)重的地區(qū)之一[9-10]。

    1.2 荒漠化風(fēng)險評估模型

    為適應(yīng)綜合地形、氣候、植被、土壤、人口、經(jīng)濟(jì)等多種環(huán)境質(zhì)量因素,在環(huán)境敏感區(qū)法(ESAs)基礎(chǔ)上,提出基于環(huán)境敏感區(qū)指標(biāo)的黃土高原荒漠化風(fēng)險評估模型,建立對黃土高原不同區(qū)域荒漠化的評價機(jī)制。模型包含5組指標(biāo),即:土壤質(zhì)量指數(shù)(SQI)、氣候質(zhì)量指數(shù)(CQI)、植被質(zhì)量指數(shù)(VQI)、土地利用與管理質(zhì)量指數(shù)(LU-MQI)和人類脅迫指數(shù)(HPI)。其中,土壤質(zhì)量指數(shù)SQI的計算包含成土母質(zhì)、土壤質(zhì)地、碎石含量、土壤深度、排水條件和坡度等要素;氣候質(zhì)量指數(shù)CQI,主要通過降水、干旱指數(shù)、坡向和降雨侵蝕力等指標(biāo)得出;植被質(zhì)量指數(shù)VQI包含植被覆蓋度、火災(zāi)風(fēng)險性、土壤侵蝕保護(hù)和植被耐旱性等要素;土地利用與管理質(zhì)量指數(shù)LU-MQI,主要通過耕地、草地、林地和灌叢的面積以及水土保持措施的重點(diǎn)區(qū)域計算得出;人類脅迫指數(shù)HPI主要根據(jù)人口數(shù)量、人口密度和農(nóng)業(yè)人口所占比重來確定[11]。與Gad 等[12]2008年提出的荒漠化敏感指數(shù)模型(DSI)相比,ESAI模型包含了兩項新的指標(biāo),即土地利用與管理質(zhì)量指數(shù)和人類脅迫指數(shù),而DSI模型可以描述為包含土壤、氣候和植被等生物物理要素的環(huán)境敏感區(qū)指標(biāo)(ESAI_SCV):

    ESAI_SCV=(SQI·CQI·VQI)1/3

    (1)

    據(jù)此,基于環(huán)境敏感區(qū)指數(shù)(ESAI)的黃土高原荒漠化風(fēng)險評估模型可以表示為:

    ESAI=(SQI·CQI·VQI·LU_MQI·HPI)1/5

    (2)

    為了反映土地退化的程度,我們把ESAI模型按照敏感程度由高到低依次分為8類:極高、高度、較高、較低、低度、極低、潛在敏感和不敏感。

    其中,SQI(土壤質(zhì)量指數(shù))是在原來的環(huán)境風(fēng)險指數(shù)模型ESI[13]的基礎(chǔ)上估算出來的,該模型考慮了母巖(PM)、土壤質(zhì)地(ST)、碎石(RF)、土壤深度(SD)、排水條件(Drainage)和坡度(Slope)等6個影響因子(公式3):

    SQI=(PM·ST·RF·SD·Drainage·Slope)1/6

    (3)

    其中,母巖(PM)中酸性巖母質(zhì)(含石英、正長石、白云母等淺色礦物)的抗風(fēng)化能力較強(qiáng),而基性和超基性巖母質(zhì)(含角閃石、輝石、黑云母等深色礦物)的抗風(fēng)化能力較弱。土壤質(zhì)地(ST)是土壤顆粒的粗細(xì)狀況,顆粒越粗,往往地表越貧瘠,對應(yīng)的環(huán)境風(fēng)險程度也比較高,反之亦然。土壤中碎石含量(RF)的多少會影響土壤受侵蝕的程度,當(dāng)土體中存在中等風(fēng)化程度變質(zhì)巖碎石時,影響了坡面土體蓄水能力與產(chǎn)流方式,同一雨量下,雨強(qiáng)小且歷時長的降雨徑流系數(shù)可能??;對于存在風(fēng)化程度高的沉積巖碎屑的土體,坡上土體易在沉積巖層或土巖界面形成壤中流,在坡下則有利于徑流入滲;對于橢圓卵石存在的土體,坡上易形成壤中流,坡下淺層土體易蓄積水量[14]。土壤質(zhì)量指數(shù)(SQI)的其余3項因素與土壤的環(huán)境質(zhì)量關(guān)系表現(xiàn)為:當(dāng)土壤深度越大,排水條件越好,坡度越小時,土壤的環(huán)境質(zhì)量則越高。

    CQI(氣候質(zhì)量指數(shù))是由年平均降水量P、坡向(Aspect)、干旱指數(shù)BGI(Bagnouls-Gaussen aridity index)[15](公式5)和通用土壤流失方程(USLE)中的降雨侵蝕力因子(R)[16]共同決定的(公式4)。

    SQI=(P·Aspect·BGI·R)1/4

    (4)

    降水量(P)是衡量區(qū)域氣候的重要標(biāo)志,但不是唯一的:氣溫會影響當(dāng)?shù)氐恼羯⒘?,故需要結(jié)合降水和氣溫2個要素來反映區(qū)域的干旱程度;由于不同坡向(Aspect)的植被種類不同,使得熱量和水分交換在時間和空間上存在顯著差異,多項研究表明,陽坡具有更南地區(qū)的氣候特點(diǎn),而陰坡則具有更北地區(qū)的氣候特點(diǎn)。因此坡向也成為反映氣候質(zhì)量的一個重要指標(biāo)。

    (5)

    式中:ti——月平均氣溫(i=1~12);pi——月平均降雨量;k——系數(shù),表示一年中各月降雨量在數(shù)值上小于兩倍氣溫時月份的個數(shù)。

    VQI(植被質(zhì)量指數(shù))根據(jù)植被覆蓋度(VC),火災(zāi)風(fēng)險性(FR),土壤侵蝕保護(hù)(SEP)和植物耐旱性(VRD)等要素估算得出(公式6),而后三者都是依據(jù)植被類型直接進(jìn)行量化的[16],其中:就火災(zāi)風(fēng)險性(FR)來說,果園和作物相比草地、灌木和落葉林等植被發(fā)生火災(zāi)的可能性較小,而針葉林發(fā)生火災(zāi)的可能性最大;對土壤侵蝕保護(hù)(SEP)能力而言,由強(qiáng)到弱依次是針葉林、灌叢、草原、落葉林、果園和農(nóng)作物;植物的耐旱性(VRD)最強(qiáng)的是針葉林和落葉林,果園和草原的次之,農(nóng)作物的耐旱性最差。

    VQI=(VC·FR·SEP·VRD)1/4

    (6)

    LU_MQI(土地利用與管理質(zhì)量指數(shù))包含土地利用(LUI)和土地管理政策(MQI)2種指標(biāo)(公式7)。而土地利用指標(biāo)(LUI)又分為作物用地(CLUI)、牧場用地(PLUI)和森林用地(FLUI)3種指標(biāo)(公式8)。其中,CLUI(作物用地)是綜合農(nóng)作物的耕作強(qiáng)度、灌溉強(qiáng)度、農(nóng)業(yè)機(jī)械化程度、農(nóng)用化肥使用數(shù)量以及植被保護(hù)產(chǎn)品的使用等5項指標(biāo)得到的;PLUI(牧場用地)考慮了整個研究區(qū)內(nèi)牧業(yè)用地的面積和牲畜頭數(shù);FLUI(森林用地)是通過森林覆蓋度估算得出的。MQI表示對土地的保護(hù),主要通過退耕還林、修建梯田和淤地壩等水土保持措施[17]估算出來的。

    LU_MQI=(LUI·MQI)1/2

    (7)

    LUI=(CLUI·PLUI·FLUI)1/3

    (8)

    人類脅迫指數(shù)(HPI)表示人類對環(huán)境的負(fù)面影響。由于黃土高原地區(qū)相比中國其他地區(qū)游客數(shù)量較少,沒有考慮原有HPI指標(biāo)中旅游業(yè)對土地退化的影響。因此,將ANPA的人類脅迫指數(shù)(HPI)[18]進(jìn)行了修訂(公式9)。

    HPI=(PD·RP·EA)1/3

    (9)

    式中:PD——人口密度; RP——年底常住人口; EA——農(nóng)業(yè)人口。

    1.3 數(shù)據(jù)采集和等級劃分標(biāo)準(zhǔn)

    通過查閱數(shù)據(jù)庫和基層調(diào)研,我們獲取了一系列的自然地理和社會社會經(jīng)濟(jì)數(shù)據(jù),主要有地形、氣候、土壤、植被等遙感數(shù)據(jù),人口、經(jīng)濟(jì)等統(tǒng)計數(shù)據(jù)以及政策方面的歷史資料等(表1)。

    所有這些數(shù)據(jù)都是為了建立基于ESAI的荒漠化風(fēng)險評估模型。不同數(shù)據(jù)或圖層的創(chuàng)建都是為了構(gòu)建5種環(huán)境質(zhì)量指數(shù)(SQI,CQI,VQI,LU_MQI和HPI),而這些質(zhì)量指數(shù)都是借助GIS平臺計算的。每個環(huán)境質(zhì)量指標(biāo)中不同變量的量化主要依據(jù)Ferrara等[15]2005年提出的方案。為了對量化以后的環(huán)境質(zhì)量數(shù)據(jù)進(jìn)行等級劃分,參考意大利學(xué)者Ladisa[1]的劃分標(biāo)準(zhǔn),結(jié)合黃土高原的實(shí)際情況,最終確定了ESAI模型中各項指標(biāo)的等級劃分標(biāo)準(zhǔn)(表2)。

    表1 ESAI模型中相關(guān)數(shù)據(jù)及其來源

    表2 ESAI模型各項指標(biāo)的等級劃分標(biāo)準(zhǔn)

    2 結(jié)果與分析

    2.1 不同環(huán)境質(zhì)量指數(shù)分析

    由于受各種自然條件和人類活動的影響,所以荒漠化是一個復(fù)雜的過程。通過ESAI模型對黃土高原荒漠化風(fēng)險各指標(biāo)進(jìn)行評估,結(jié)果如附圖9a—9e所示(SQI,CQI,VQI,LU_MQI和HPI)。

    2.1.1 土壤質(zhì)量指數(shù)(SQI) 附圖9a表明,黃土高原1/2以上的土壤為中等質(zhì)量(52.5%),接近1/2(44.2%)的為高品質(zhì)土壤,而低質(zhì)量土壤僅占了3.3%。其中,高質(zhì)量區(qū)主要分布在黃土高原的中心地段。分布如此廣泛的高質(zhì)量土壤是因為黃土高原土層深厚,被稱為世界上黃土層最厚的地方,在人類頻繁的黃土高原活動作用下黃土巖石被改造成為土壤,多數(shù)地區(qū)有效土壤厚度達(dá)100 cm。對于小面積的低質(zhì)量土壤,主要是由于其土層較淺,地表貧瘠且土壤排水條件較差。

    2.1.2 氣候質(zhì)量指數(shù)(CQI) 通過對全國61個氣象站點(diǎn)的氣溫和降水?dāng)?shù)據(jù)進(jìn)行空間插值,并經(jīng)過一系列的計算得到氣候質(zhì)量指數(shù),并進(jìn)行氣候質(zhì)量等級劃分(附圖9b)。附圖9b表明,氣候質(zhì)量呈東南向西北遞減的趨勢,同時在寧夏—蘭州—海原—定邊—綏德—河曲—呼和浩特等地區(qū)存在一條明顯的界限。其中,低質(zhì)量氣候區(qū)主要分布在邊界以北(23.7%),高質(zhì)量氣候區(qū)主要分布在編輯以南(29.8%),而中等質(zhì)量氣候區(qū)鑲嵌分布在以上兩者之間(46.5%)。這種空間格局和黃土高原年平均降水量的分布格局基本一致,但因考慮了氣溫要素,使得氣候質(zhì)量的空間差距進(jìn)一步縮??;另外,不同的坡向?qū)е陆邓蜔崃恐匦路峙?,使整個氣候質(zhì)量在微觀尺度上也表現(xiàn)出空間差異性。

    2.1.3 植被質(zhì)量指數(shù)(VQI) 植被質(zhì)量指數(shù)顯示,不到1/5的黃土高原地區(qū)植被質(zhì)量較高(19.0%),將近1/2的地區(qū)是中等質(zhì)量(42.1%),大約38.0%的區(qū)域植被質(zhì)量較差(附圖9c)。高質(zhì)量植被區(qū)位于黃土高原東南部,主要由分布在秦嶺,子午嶺和太行山北部的高植被覆蓋率的森林決定。由于黃土高原整體的植被覆蓋度較低,大部分地區(qū)都屬于中等或低植被質(zhì)量區(qū),尤其是以草地為主的內(nèi)蒙古河套平原。另外,從東南向西北的過渡地帶由于分布了大量的農(nóng)作物,其脆弱的植被類型無論從耐旱性或是保護(hù)地表不被侵蝕等方面都比較差,因此,導(dǎo)致了大面積低植被質(zhì)量區(qū)域的分布。

    2.1.4 土地利用與管理質(zhì)量指數(shù)(LU_MQI) 有關(guān)土地利用與管理質(zhì)量的高低,主要通過縣域尺度上耕地面積、放牧面積和造林面積等進(jìn)行估算(附圖9d)。土地利用和管理質(zhì)量指數(shù)顯示:黃土高原幾乎1/2的區(qū)域(48.3%)處于高水平甚至極高水平,主要分布在黃土高原中部由東北向西南延伸的廣闊地帶。因為該區(qū)域是中國科學(xué)院(CAS)劃分的水土保持重點(diǎn)區(qū)域,自20世紀(jì)80年代以來,實(shí)施了大規(guī)模的坡改梯、淤地壩等水土保持工程措施以及大量的植樹造林、退耕還林還草生態(tài)措施。研究結(jié)果顯示,中等質(zhì)量占36.1%,與高質(zhì)量區(qū)所占比例較為接近,主要分布在黃土高原的東南部渭河平原一帶,另外在隴中、青東、山西和蒙北均有分布。只有12.2%的低質(zhì)量集中分布在內(nèi)蒙古鄂托克旗,鄂托克前旗和烏審旗等區(qū)域,該區(qū)域分布了大面積的荒漠用地,其次是牧業(yè)用地,但因單位牧草地面積上的載畜量較大極易造成土地退化。鑒于此,在加強(qiáng)水土保持重點(diǎn)區(qū)域的土地利用管理的同時,也必須重視黃土高原其他地區(qū)尤其是內(nèi)蒙古地區(qū)的土地管理工作。

    2.1.5 人類脅迫指數(shù)(HPI) 人類脅迫指數(shù)等級圖中顯示,黃土高原近1/2區(qū)域(46.2%)受人類干擾程度較弱,39.8%屬于中等脅迫,只有約14.0%為高度脅迫(附圖9e)。黃土高原幾乎涵蓋了整個黃河流域,流域良好的土壤和水分條件,形成了黃土高原幾千年的農(nóng)耕文化,農(nóng)業(yè)人口比重較大,同時也產(chǎn)生了撂荒制、壓青休閑制、粗放輪作制、保護(hù)性耕作制等多種耕作制度。因此,人口密度、農(nóng)業(yè)人口數(shù)量的空間差異性以及農(nóng)業(yè)耕作制度的多樣化共同決定了黃土高原HPI的空間分布格局。

    2.2 黃土高原荒漠化風(fēng)險評價

    2.2.1 基于生物物理要素的黃土高原荒漠化風(fēng)險評價 土壤、氣候和植被決定了黃土高原荒漠化風(fēng)險程度的基本格局,為此,我們綜合SQI,CQI和VQI 3種環(huán)境質(zhì)量指標(biāo),得到環(huán)境敏感區(qū)生物物理指標(biāo)ESAI_SCV,但該指標(biāo)僅表示黃土高原荒漠化的生物物理因素。ESAI_SCV等級圖顯示:黃土高原大約1/4的區(qū)域(25.2%)為高風(fēng)險區(qū),屬于嚴(yán)重荒漠化,幾乎2/3的區(qū)域(62.8%)是輕微荒漠化,11.5%的地區(qū)為潛在荒漠化,只有0.5%的地區(qū)無荒漠化現(xiàn)象(附圖9f)??傮w來看,黃土高原荒漠化程度在西北地區(qū)顯著高于東南地區(qū)。其中,內(nèi)蒙古杭錦旗是荒漠化嚴(yán)重的地方,其次是陜西省的榆林市。相比之下,陜西、山西的許多縣市都是荒漠化的低風(fēng)險區(qū),主要表現(xiàn)為微弱或潛在沙漠化。黃土高原位于中國半濕潤氣候區(qū)與半干旱干旱氣候區(qū)的過渡帶,既是氣候變化的敏感區(qū),又是環(huán)境脆弱區(qū)。根據(jù)《聯(lián)合防治荒漠化公約》對荒漠化的定義,干旱半干旱和亞濕潤地區(qū)是荒漠化發(fā)生的源地,在全球變暖的大背景下,干旱半干旱區(qū)的氣候變化是影響荒漠化的一個重要原因。植被和土壤相互影響,且兩者都在不同的氣候條件下表現(xiàn)出空間異質(zhì)性。

    2.2.2 黃土高原荒漠化風(fēng)險綜合評價 為了制定黃土高原荒漠化防治的有效措施,還需要考慮人類活動對荒漠化的影響。因此,將土地利用與管理和人類脅迫等社會經(jīng)濟(jì)指標(biāo)同土壤、氣候、植被等生物物理指標(biāo)相結(jié)合,最終建立了荒漠化風(fēng)險評估的環(huán)境敏感性評估指標(biāo)(ESAI)(附圖9g)。列入人類誘發(fā)因素后,改變了黃土高原荒漠化風(fēng)險區(qū)的原有格局。其中,極低、低度和極高度敏感區(qū)減少了5.6%,1.1%和3.8%;與此同時,較低和較高度敏感區(qū)增加了4.4%和4.5%(附圖9h,表3)。結(jié)合黃土高原環(huán)境敏感區(qū)生物物理指標(biāo)ESAI_SCV(附圖9f)和環(huán)境敏感區(qū)指標(biāo)ESAI(附圖9g)發(fā)現(xiàn),在黃土高原的西北部的大部分地區(qū),荒漠化風(fēng)險程度在考慮人類指標(biāo)后明顯降低,但在黃土高原東南部人口相對較稠密的地區(qū),荒漠化風(fēng)險程度略有增加。以上結(jié)果表明,一方面,人類活動縮小了不同程度荒漠化之間的差距,說明人類實(shí)施的水土保持措施(植樹造林、封山禁牧生態(tài)措施、坡改梯和淤地壩工程措施等)在荒漠化防治過程中已經(jīng)發(fā)揮作用;另一方面,較高的人口壓力和農(nóng)業(yè)人口比重可能加強(qiáng)荒漠化的風(fēng)險程度,在黃土高原西北部內(nèi)蒙古地區(qū)人口稀少,因人類活動對環(huán)境產(chǎn)生的負(fù)面影響較低,相反,在黃土高原地區(qū)其他城市(如:蘭州、西安、寶雞、咸陽、洛陽等地區(qū))的人口壓力較大,人類對環(huán)境的過渡索取導(dǎo)致了荒漠化風(fēng)險程度加劇。因此,在人類活動的影響下,ESAI與ESAI-SCV比較,形成了ESAI的新格局。

    表3 ESAI和ESAI_SCV之間的差異

    3 結(jié) 論

    (1) 基于土壤、氣候、植被3種要素的環(huán)境敏感區(qū)生物物理指標(biāo)(ESAI_SCV)決定了黃土高原荒漠化風(fēng)險程度的基本格局。其中,土壤因素在某種程度上降低了黃土高原中部的荒漠化風(fēng)險程度,氣候因素決定了黃土高原荒漠化風(fēng)險程度由東南向西北遞減的趨勢;植被是指示荒漠化風(fēng)險程度最顯著的指標(biāo),因為植被質(zhì)量指數(shù)圖和環(huán)境敏感區(qū)指標(biāo)圖之間具有相似的空間分布模式。

    (2) 黃土高原南部荒漠化風(fēng)險程度較高,這主要?dú)w因于較大的人口密度和人類作用于土地的程度,而在北部黃土溝壑區(qū)由于實(shí)施了水土保持工程措施和生態(tài)措施,大幅度地降低了荒漠化風(fēng)險程度。因此,在加強(qiáng)荒漠化防治工程的同時,還需要控制人口數(shù)量來降低人口壓力,才能最終減弱荒漠化風(fēng)險程度。

    (3) 通過對比ESAI和ESAI_SCV兩種荒漠化風(fēng)險指標(biāo)發(fā)現(xiàn),一方面,包含人類社會的生態(tài)系統(tǒng)比自然生態(tài)系統(tǒng)更加多變;另一方面,人類活動縮小了不同程度荒漠化之間的差距。因此,我們應(yīng)該時刻關(guān)注人類活動對黃土高原荒漠化的影響,及時制定包括人類活動因素在內(nèi)的積極有效的荒漠化防治措施。

    [1] Ladisa G, Todorovic M, Trisorio L G. A GIS-based approach for desertification risk assessment in Apulia region, SE Italy[J]. Physics and Chemistry of the Earth, 2012,49(3):103-113.

    [2] 楊軍.柴達(dá)木盆地荒漠化成因分析及產(chǎn)生機(jī)理探討[D].陜西 西安:長安大學(xué),2004.

    [3] Kosmas Costas, Ferrara Agostino, Briassouli Helen, et al. Methodology for mapping ESAs to desertification[M]∥Kosmas C, Kirkby M. Geeson N. The MEDALUS Project-Mediterranean Desertification and Land Use. Manual on Key Indicators of Desertification and Mapping Environmentally Sensitive Areas to desertification, EUR18882, 1999:31-47.

    [4] Coscarelli Roberto, Minervino Ivana. Methods for the characterization of areas sensitive to desertification: An application to the Calabrian territory[C]∥Sorriso-Valvo. Proceedings of the international conference held at Solsona[J]. Catalonia: IAHS, 2005:298-299.

    [5] 張希彪,上官周平,王金成,等.隴東黃土高原農(nóng)牧交錯帶土地荒漠化過程中自然與人為因素的定量分析[J].水土保持通報,2013,33(2):203-208.

    [6] 中國科學(xué)院黃土高原綜合科學(xué)考察隊.黃土高原地區(qū)資源環(huán)境社會經(jīng)濟(jì)數(shù)據(jù)集[M].北京:中國經(jīng)濟(jì)出版社,1992.

    [7] Zhao Yifei, Zou Xinqing, Zhang Jianxiang, et al. Spatio-temporal variation of reference evapotranspiration and aridity index in the Loess Plateau region of China, during 1961—2012[J]. Quaternary International, 2014,349(10):196-206.

    [8] Li Zhi, Liu Wenzhao, Zhang Xunchang, et al. Impacts of land use change and climate variability on hydrology in an agricultural catchment on the Loess Plateau of China[J]. Journal of Hydrology, 2009,377(1/2):35-42.

    [9] Zhang Xunchang, Liu Wenzhao. Simulating potential response of hydrology, soil erosion, and crop productivity to climate change in Changwu tableland region on the Loess Plateau of China [J]. Agricultural and Forest Meteorology, 2005,131(3/4):127-142.

    [10] Sun Wenyi, Shao Quanqin, Liu Jiyuan, et al. Assessing the effects of land use and topography on soil erosion on the Loess Plateau in China[J]. Catena, 2014,121(7):151-163.

    [11] Agostino F, Luca S, Adele S A, et al. Performance evaluation and cost assessment of a key indicator system to monitor desertification vulnerability[J]. Ecological Indicators, 2012,23(5):123-129.

    [12] Gad A, Lofty I. Use of remote sensing and GIS in mapping the environmental sensitivity areas for desertification of Egyptian territory[J]. Earth Discussions, 2008,3(2):41-85.

    [13] Ferrara Agostino. Expert system for evaluating the Environmental Sensitivity Index(ESI)of a local area [EB/OL]. (2005-02-03)[2005-10-23]. http:∥www. kcl. ac. uk/projects/desertlinks/indicator system/index. htm.

    [14] 王慧芳,邵明安,王國麗.黃土高原土石山區(qū)碎石分布特征及其導(dǎo)水性質(zhì)分析[J].灌溉排水學(xué)報,2010,29(3):116-120.

    [15] Bagnouls Franco. L’indice xérothermique[R]. Paris: Bulletin de l’Association de Géographes francais, 1952.

    [16] Sharply A N, Williams J R. Epic-erosion Productivity Impact Calculator(I): Model Documentation[M]. Washington D C: US Department of Agriculture,1990.

    [17] 石建華,喻理飛,孫保平.陜北地區(qū)退耕還林生態(tài)健康評價分析研究:以吳起縣為例[J].水土保持學(xué)報,2015,29(6):332-336.

    [18] ANPA. Centro Tematico Nazionale Conservazione della Natura. Selezione di indicatori ambientali per itemi relativi alla biosfera[C]. Roma: RTI CTN_CON, 2000.

    ESAI Based Assessment of Desertification Risk in Loess Plateau

    ZHANG Jianxiang1,2, ZHANG Duoyong3, LIU Wanfeng1,2, HU Aiping1, WANG Dong1, Hu Jianqiao1

    (1.CollegeofCivilEngineering,LongdongUniversity,Qingyang,Gansu745000,China; 2.ProvincialKeyLaboratoryofLoessEngineeringPropertiesandApplicationsinGansuProvinceUniversities,Qingyang,Gansu745000,China; 3.ResearchCenterofDesertificationPreventionatQingyanCity,Qingyang,Gansu745000,China)

    [Objective] The risk of desertification was evaluated in Loess Plateau to provide scientific basis for the restoration, construction and protection of ecological environment. [Methods] Data of landform, climate, vegetation, soil and other socio-economic data were collected. Spatialization and digitalization were conducted using GIS and remote sensing. Upon which, desertification reasons were analyzed, and evaluation indices and risk assessment model were framed based on index of environmentally sensitive areas. [Results] (1) According to the bio-physical index in the framed model, under the scenario only considering natural factors as soil, climate and vegetation, 25.2% of Loess Plateau was determined as high desertification risk area, where desertification was worst; 62.8% and 11.5% of Loess Plateau were determined as moderate and potential desertification risk areas; only 0.5% was considered having no risk. (2) If both natural factors and human interference were considered, risk area changed: coverages of extreme low, low and extreme high risk levels decreased by 5.6%,1.1% and 3.8%, respectively; whereas, coverages of lower and higher risk levels increased by 4.4% and 4.5%. [Conclusion] (1) The ESAI model can well explain the spatial distribution pattern of desertification risk of Loess Plateau, and the degree of desertification is gradually weakened from the northwest to the southeast. (2) Human activities have upset the long-term evolved stability of the natural eco-system to some extent, and have narrowed the gap between different risk levels of desertification.

    desertification; environmental sensitive areas index(ESAI); geographic information systems(GIS); Loess Plateau

    2016-03-25

    2016-07-28

    國家自然科學(xué)基金項目“清代同治以來黃土高原馬蓮河流域荒漠化風(fēng)險評估與防治研究”(31460090); 2015年甘肅省社科規(guī)劃項目(YB095); 甘肅省社科規(guī)劃項目(14YB108); 西峰區(qū)科技局科技計劃項目(XK2016-06); 慶陽市科技局青年科技基金項目(ZJ2014-08); 甘肅省哲學(xué)社會科學(xué)規(guī)劃項目(YB093)

    張建香(1987—),女(漢族),甘肅省華池縣人,碩士,講師,主要從事遙感與地理信息系統(tǒng)技術(shù)應(yīng)用研究。E-mail:zhangjianxiang_67@163.com。

    張多勇(1966—),男(漢族),甘肅省華池縣人,博士,教授,主要從事歷史地理研究。E-mail:bbbjjj2@163.com。

    A

    1000-288X(2017)02-0339-06

    TP79, X171.1

    猜你喜歡
    敏感區(qū)黃土高原荒漠化
    基于GIS的贛南地區(qū)城鎮(zhèn)生態(tài)安全格局研究
    ——以贛州市龍南縣為例
    四川建材(2020年7期)2020-07-26 06:02:00
    荒漠化地區(qū)復(fù)合生態(tài)系統(tǒng)管理——以阿拉善盟荒漠化治理為例
    科學(xué)(2020年6期)2020-02-06 08:59:52
    防治土地荒漠化 推動綠色發(fā)展——寫在第二十五個世界防治荒漠化和干旱日
    綠色中國(2019年13期)2019-11-26 07:10:54
    我國召開第25個世界防治荒漠化與干旱日紀(jì)念大會
    綠色中國(2019年13期)2019-11-26 07:10:52
    在某個敏感區(qū)如何協(xié)調(diào)區(qū)域發(fā)展與環(huán)境保護(hù)的探究
    牡丹江流域自然地理敏感區(qū)劃分研究
    我國荒漠化和沙化面積連續(xù)10年實(shí)現(xiàn)“雙縮減”
    綠色中國(2016年1期)2016-06-05 09:03:00
    選舉 沸騰了黃土高原(下)
    公民與法治(2016年3期)2016-05-17 04:09:00
    選舉沸騰了黃土高原(上)
    公民與法治(2016年1期)2016-05-17 04:07:56
    ILS臨界區(qū)與敏感區(qū)淺析
    欧美 亚洲 国产 日韩一| 亚洲av成人精品一区久久| 夜夜骑夜夜射夜夜干| 好男人视频免费观看在线| 美女内射精品一级片tv| 精品人妻熟女av久视频| 亚洲精品久久午夜乱码| 噜噜噜噜噜久久久久久91| 国产精品秋霞免费鲁丝片| 国产黄色免费在线视频| 欧美日韩一区二区视频在线观看视频在线| 国产成人一区二区在线| 精品亚洲成a人片在线观看| 91精品伊人久久大香线蕉| a级毛片免费高清观看在线播放| 黑人高潮一二区| 免费高清在线观看视频在线观看| 我的老师免费观看完整版| 免费观看在线日韩| 国产毛片在线视频| a级一级毛片免费在线观看| 国产毛片在线视频| 欧美xxⅹ黑人| 最新的欧美精品一区二区| 最黄视频免费看| 日韩三级伦理在线观看| 亚洲色图综合在线观看| 人妻少妇偷人精品九色| 欧美日韩亚洲高清精品| 日韩一区二区视频免费看| 黄色一级大片看看| 亚洲精品亚洲一区二区| 色吧在线观看| 色视频在线一区二区三区| 一本久久精品| 亚洲美女搞黄在线观看| 中文在线观看免费www的网站| 欧美老熟妇乱子伦牲交| 久久99一区二区三区| 亚洲一区二区三区欧美精品| 大码成人一级视频| 国产高清有码在线观看视频| 老司机影院成人| 亚洲人与动物交配视频| 成年av动漫网址| 精品一区二区三卡| 最近中文字幕2019免费版| freevideosex欧美| 七月丁香在线播放| 国产日韩一区二区三区精品不卡 | 99热网站在线观看| 日韩免费高清中文字幕av| 久久久久久久久久久丰满| 成人免费观看视频高清| 色吧在线观看| 国产黄色视频一区二区在线观看| 国产av一区二区精品久久| 亚洲av欧美aⅴ国产| 高清av免费在线| 亚洲国产最新在线播放| 亚洲av成人精品一区久久| 免费观看av网站的网址| 国产男女内射视频| 热99国产精品久久久久久7| a级片在线免费高清观看视频| 91久久精品电影网| 日韩强制内射视频| 热re99久久精品国产66热6| 欧美老熟妇乱子伦牲交| 22中文网久久字幕| 日韩不卡一区二区三区视频在线| 成人18禁高潮啪啪吃奶动态图 | 黄色视频在线播放观看不卡| 精品少妇黑人巨大在线播放| 日韩精品免费视频一区二区三区 | 99热国产这里只有精品6| 如日韩欧美国产精品一区二区三区 | 毛片一级片免费看久久久久| 啦啦啦中文免费视频观看日本| 国产精品免费大片| 国产伦精品一区二区三区视频9| 国产免费福利视频在线观看| 国产免费又黄又爽又色| 在线亚洲精品国产二区图片欧美 | 国产 一区精品| 免费大片黄手机在线观看| 亚洲av国产av综合av卡| 国产成人freesex在线| 精品视频人人做人人爽| 大又大粗又爽又黄少妇毛片口| 麻豆精品久久久久久蜜桃| av国产精品久久久久影院| 多毛熟女@视频| 九九久久精品国产亚洲av麻豆| 街头女战士在线观看网站| 肉色欧美久久久久久久蜜桃| h日本视频在线播放| 人人妻人人爽人人添夜夜欢视频 | 精品卡一卡二卡四卡免费| 最近中文字幕高清免费大全6| 精品久久久久久电影网| 国产精品一区二区三区四区免费观看| 99热这里只有精品一区| 久久人人爽人人片av| 国产 精品1| 高清欧美精品videossex| 国产成人a∨麻豆精品| 国产综合精华液| 少妇人妻久久综合中文| 国产成人精品无人区| 在线播放无遮挡| 一区二区三区乱码不卡18| 精品久久久久久久久av| 国产精品人妻久久久影院| 少妇人妻一区二区三区视频| 精品少妇黑人巨大在线播放| 久久亚洲国产成人精品v| 99国产精品免费福利视频| 黑丝袜美女国产一区| 久久鲁丝午夜福利片| 欧美人与善性xxx| 国产av码专区亚洲av| 黄色配什么色好看| 久久99蜜桃精品久久| 我的女老师完整版在线观看| 麻豆乱淫一区二区| 国产精品.久久久| 色婷婷av一区二区三区视频| 麻豆成人午夜福利视频| 久久鲁丝午夜福利片| 日本wwww免费看| 精品少妇久久久久久888优播| 狠狠精品人妻久久久久久综合| 亚洲成人av在线免费| 亚洲精品乱久久久久久| 狂野欧美激情性bbbbbb| 亚洲中文av在线| 久久久久久伊人网av| 女性被躁到高潮视频| 日韩欧美一区视频在线观看 | 国产精品一二三区在线看| 欧美老熟妇乱子伦牲交| 国产成人免费观看mmmm| 国产 精品1| 国产毛片在线视频| av免费在线看不卡| 久久6这里有精品| 18禁动态无遮挡网站| 久久国产精品大桥未久av | 国产精品不卡视频一区二区| 99久久精品热视频| 最近2019中文字幕mv第一页| 女人久久www免费人成看片| 国产成人91sexporn| 中文乱码字字幕精品一区二区三区| 国产av一区二区精品久久| 亚洲美女黄色视频免费看| 久久国产亚洲av麻豆专区| 丝袜喷水一区| 国产精品福利在线免费观看| 亚洲va在线va天堂va国产| 夜夜爽夜夜爽视频| 一区在线观看完整版| 卡戴珊不雅视频在线播放| 一级黄片播放器| 久热这里只有精品99| 亚洲精品视频女| 久久综合国产亚洲精品| 久久精品国产a三级三级三级| 另类亚洲欧美激情| 免费播放大片免费观看视频在线观看| 男男h啪啪无遮挡| 热99国产精品久久久久久7| 国产日韩欧美视频二区| 久久人人爽人人爽人人片va| 丰满迷人的少妇在线观看| 久久97久久精品| 在线观看免费高清a一片| 狂野欧美白嫩少妇大欣赏| 亚洲精品国产av蜜桃| 日韩欧美一区视频在线观看 | 一级a做视频免费观看| 91精品伊人久久大香线蕉| 国产亚洲一区二区精品| 又黄又爽又刺激的免费视频.| 人妻制服诱惑在线中文字幕| 天天躁夜夜躁狠狠久久av| 国产男人的电影天堂91| 国产精品嫩草影院av在线观看| 中国美白少妇内射xxxbb| 一区二区av电影网| 狂野欧美激情性bbbbbb| 欧美区成人在线视频| 91久久精品国产一区二区三区| 中文字幕免费在线视频6| 天堂中文最新版在线下载| 我要看日韩黄色一级片| 99热这里只有是精品50| 日韩中字成人| 在线免费观看不下载黄p国产| 日韩,欧美,国产一区二区三区| 亚洲一级一片aⅴ在线观看| 国产精品一二三区在线看| 黑丝袜美女国产一区| 激情五月婷婷亚洲| 久久久久久久亚洲中文字幕| 国产极品粉嫩免费观看在线 | 丝袜在线中文字幕| 在线看a的网站| 国产综合精华液| 亚洲av成人精品一区久久| 国产真实伦视频高清在线观看| 高清欧美精品videossex| 色5月婷婷丁香| 视频区图区小说| a 毛片基地| 亚洲性久久影院| 黑人高潮一二区| 免费看日本二区| 男女免费视频国产| 好男人视频免费观看在线| 人人妻人人看人人澡| 极品人妻少妇av视频| 免费黄色在线免费观看| 亚洲欧美一区二区三区黑人 | 免费观看在线日韩| 国产免费视频播放在线视频| 97超碰精品成人国产| 精品国产乱码久久久久久小说| 成年人午夜在线观看视频| 内射极品少妇av片p| 自拍欧美九色日韩亚洲蝌蚪91 | 日韩大片免费观看网站| 成年av动漫网址| 亚洲av国产av综合av卡| 一级毛片我不卡| 色5月婷婷丁香| 五月玫瑰六月丁香| 国产在线一区二区三区精| 午夜日本视频在线| 国产精品嫩草影院av在线观看| 在线播放无遮挡| 女人精品久久久久毛片| 丰满饥渴人妻一区二区三| 全区人妻精品视频| 啦啦啦视频在线资源免费观看| 国产欧美日韩精品一区二区| 久久久久视频综合| 九九爱精品视频在线观看| 日日啪夜夜撸| 乱系列少妇在线播放| 国产白丝娇喘喷水9色精品| 91久久精品国产一区二区三区| 新久久久久国产一级毛片| 国产精品无大码| 亚洲国产精品成人久久小说| 老司机影院毛片| 一级二级三级毛片免费看| 午夜影院在线不卡| 久久久久久伊人网av| 欧美3d第一页| 国产亚洲5aaaaa淫片| 一区二区三区免费毛片| 视频区图区小说| 中文乱码字字幕精品一区二区三区| 亚洲中文av在线| 国产精品一区二区在线不卡| 日韩电影二区| 18禁裸乳无遮挡动漫免费视频| 美女视频免费永久观看网站| 国产黄色免费在线视频| a 毛片基地| 大陆偷拍与自拍| 一级二级三级毛片免费看| 国产淫语在线视频| 大码成人一级视频| av在线老鸭窝| 少妇被粗大猛烈的视频| 国产精品国产三级国产专区5o| 亚洲国产精品一区二区三区在线| 亚洲精品国产av成人精品| 成年人免费黄色播放视频 | 久久久a久久爽久久v久久| 免费观看a级毛片全部| 久久精品夜色国产| 天天操日日干夜夜撸| 嫩草影院入口| 国产白丝娇喘喷水9色精品| 久久综合国产亚洲精品| 久久亚洲国产成人精品v| 欧美精品国产亚洲| 久久精品国产a三级三级三级| 如日韩欧美国产精品一区二区三区 | 91精品国产国语对白视频| 日韩成人伦理影院| 国产亚洲午夜精品一区二区久久| 亚洲国产毛片av蜜桃av| 丰满人妻一区二区三区视频av| 特大巨黑吊av在线直播| 黄片无遮挡物在线观看| 99热6这里只有精品| 国产有黄有色有爽视频| 午夜老司机福利剧场| 人妻制服诱惑在线中文字幕| 国产黄片视频在线免费观看| 精品人妻熟女毛片av久久网站| 男女边吃奶边做爰视频| 亚洲丝袜综合中文字幕| 国产国拍精品亚洲av在线观看| 亚洲av在线观看美女高潮| 国产午夜精品久久久久久一区二区三区| 国语对白做爰xxxⅹ性视频网站| 国产成人一区二区在线| 中文精品一卡2卡3卡4更新| 亚洲av在线观看美女高潮| 国产免费一级a男人的天堂| 国产视频首页在线观看| 成人美女网站在线观看视频| 国产精品偷伦视频观看了| 日韩成人av中文字幕在线观看| 免费人成在线观看视频色| 亚洲国产精品专区欧美| 黑人巨大精品欧美一区二区蜜桃 | 日韩人妻高清精品专区| 久久久久网色| 寂寞人妻少妇视频99o| 王馨瑶露胸无遮挡在线观看| 久久久久久久久久人人人人人人| 精品国产国语对白av| 伊人久久精品亚洲午夜| 久久影院123| 97在线人人人人妻| 亚洲精品自拍成人| 99久久综合免费| 国产成人午夜福利电影在线观看| 中文欧美无线码| 老司机亚洲免费影院| 国产av精品麻豆| 中文字幕av电影在线播放| 人人妻人人爽人人添夜夜欢视频 | 黄色欧美视频在线观看| 七月丁香在线播放| 少妇精品久久久久久久| 我的女老师完整版在线观看| 久久久a久久爽久久v久久| 久久久久久久久大av| 精品99又大又爽又粗少妇毛片| av免费观看日本| 亚洲精品中文字幕在线视频 | 精品人妻熟女av久视频| 亚洲av欧美aⅴ国产| 在现免费观看毛片| 亚洲av欧美aⅴ国产| 亚洲精品日本国产第一区| av.在线天堂| 最近手机中文字幕大全| 成年美女黄网站色视频大全免费 | 亚洲第一av免费看| 亚洲av.av天堂| 亚洲不卡免费看| 久久久久精品久久久久真实原创| 亚洲精品乱码久久久v下载方式| 少妇被粗大猛烈的视频| 亚洲图色成人| 久久久久久久久大av| 亚洲国产精品国产精品| 菩萨蛮人人尽说江南好唐韦庄| 三上悠亚av全集在线观看 | 国产免费又黄又爽又色| 免费久久久久久久精品成人欧美视频 | 一级a做视频免费观看| 免费人成在线观看视频色| 午夜激情久久久久久久| 秋霞伦理黄片| 国产高清国产精品国产三级| 久久av网站| 天堂8中文在线网| 中文资源天堂在线| 观看美女的网站| 亚洲国产精品999| 亚洲国产精品一区三区| 国产成人freesex在线| av国产精品久久久久影院| 欧美日本中文国产一区发布| 精品少妇内射三级| av天堂中文字幕网| 婷婷色综合www| 插逼视频在线观看| 精品一区二区三区视频在线| 亚洲精品乱码久久久久久按摩| 亚洲欧美日韩另类电影网站| 国产极品天堂在线| 热re99久久精品国产66热6| 如日韩欧美国产精品一区二区三区 | 女性被躁到高潮视频| 成人毛片60女人毛片免费| 热99国产精品久久久久久7| 国产精品一区二区性色av| 亚洲国产精品成人久久小说| 久久久午夜欧美精品| 国产乱来视频区| 国产精品无大码| 如何舔出高潮| 九九久久精品国产亚洲av麻豆| 十八禁高潮呻吟视频 | 亚洲三级黄色毛片| 欧美国产精品一级二级三级 | 久久6这里有精品| 看十八女毛片水多多多| 亚洲综合精品二区| av网站免费在线观看视频| 少妇 在线观看| 日本与韩国留学比较| 黄色视频在线播放观看不卡| 国产黄频视频在线观看| 天堂俺去俺来也www色官网| 亚洲精品视频女| 久久精品久久久久久久性| 街头女战士在线观看网站| 欧美成人精品欧美一级黄| 99热国产这里只有精品6| 久久久久久久久久久丰满| 日韩制服骚丝袜av| 国产色婷婷99| av一本久久久久| 97超碰精品成人国产| 综合色丁香网| .国产精品久久| 欧美 亚洲 国产 日韩一| 天天操日日干夜夜撸| 久久人人爽人人片av| 国产精品无大码| 国产爽快片一区二区三区| 亚洲国产av新网站| 久久精品国产亚洲av天美| av不卡在线播放| 欧美高清成人免费视频www| 少妇人妻 视频| 国产老妇伦熟女老妇高清| 国产精品国产av在线观看| 国产亚洲最大av| 国产av一区二区精品久久| 国产国拍精品亚洲av在线观看| 欧美丝袜亚洲另类| 日韩人妻高清精品专区| 嫩草影院入口| 不卡视频在线观看欧美| 97超视频在线观看视频| 国产乱来视频区| 成人毛片a级毛片在线播放| 免费观看无遮挡的男女| 国产在线男女| 亚洲国产精品专区欧美| 妹子高潮喷水视频| 美女大奶头黄色视频| 一区二区三区免费毛片| 免费大片黄手机在线观看| a级毛片免费高清观看在线播放| 久久久久久久国产电影| 久久免费观看电影| 亚洲精品国产av成人精品| 女人久久www免费人成看片| 国产精品人妻久久久久久| 国产真实伦视频高清在线观看| 国产精品久久久久久精品电影小说| 亚洲经典国产精华液单| 9色porny在线观看| 日日摸夜夜添夜夜添av毛片| 国产成人免费无遮挡视频| 国产精品久久久久久久电影| 欧美日韩av久久| 蜜臀久久99精品久久宅男| 亚洲精品456在线播放app| 日本wwww免费看| 中文字幕久久专区| 久久人人爽av亚洲精品天堂| 久久国产精品大桥未久av | 久久久久久久久久人人人人人人| 日韩亚洲欧美综合| 亚洲美女黄色视频免费看| 高清毛片免费看| 精品一区二区免费观看| 精品人妻一区二区三区麻豆| 色哟哟·www| 国产在线视频一区二区| 国产亚洲精品久久久com| av网站免费在线观看视频| 国产精品伦人一区二区| 免费不卡的大黄色大毛片视频在线观看| 午夜视频国产福利| 高清不卡的av网站| 亚洲国产精品一区二区三区在线| 国产69精品久久久久777片| 精品久久久久久电影网| 国产av一区二区精品久久| av国产精品久久久久影院| 亚洲电影在线观看av| 秋霞在线观看毛片| 99久久精品热视频| 赤兔流量卡办理| 91精品国产国语对白视频| 国产免费一级a男人的天堂| 亚洲av综合色区一区| 少妇人妻久久综合中文| 国产精品国产av在线观看| 午夜福利视频精品| 国产永久视频网站| 香蕉精品网在线| 99久久中文字幕三级久久日本| 80岁老熟妇乱子伦牲交| 午夜福利视频精品| 精品人妻偷拍中文字幕| 熟妇人妻不卡中文字幕| 亚洲国产欧美在线一区| 国产毛片在线视频| 插阴视频在线观看视频| 只有这里有精品99| 新久久久久国产一级毛片| 国产视频首页在线观看| 久久国产精品大桥未久av | 欧美bdsm另类| 久久狼人影院| 高清不卡的av网站| av网站免费在线观看视频| 18禁动态无遮挡网站| 中文精品一卡2卡3卡4更新| 人人妻人人爽人人添夜夜欢视频 | 国产一区二区三区综合在线观看 | 国产成人a∨麻豆精品| 午夜福利网站1000一区二区三区| 18禁动态无遮挡网站| 大片电影免费在线观看免费| www.av在线官网国产| 热99国产精品久久久久久7| 免费观看在线日韩| 国产一级毛片在线| 在线观看三级黄色| 午夜福利,免费看| 一本大道久久a久久精品| 最新的欧美精品一区二区| 嘟嘟电影网在线观看| 97超视频在线观看视频| 99re6热这里在线精品视频| av不卡在线播放| 一区在线观看完整版| 男人添女人高潮全过程视频| 久久99热这里只频精品6学生| 制服丝袜香蕉在线| 18禁动态无遮挡网站| 精品久久久久久久久av| 成人免费观看视频高清| 狂野欧美白嫩少妇大欣赏| 精品久久久久久电影网| 欧美老熟妇乱子伦牲交| 高清欧美精品videossex| 亚洲熟女精品中文字幕| 男人添女人高潮全过程视频| 卡戴珊不雅视频在线播放| 一级,二级,三级黄色视频| 亚洲精品一区蜜桃| 日韩强制内射视频| 国产高清不卡午夜福利| 国产乱来视频区| 免费不卡的大黄色大毛片视频在线观看| 人人妻人人爽人人添夜夜欢视频 | 久久女婷五月综合色啪小说| 日韩一区二区三区影片| 各种免费的搞黄视频| 欧美+日韩+精品| 久久久久网色| 日韩不卡一区二区三区视频在线| 国产成人午夜福利电影在线观看| 插阴视频在线观看视频| 成人黄色视频免费在线看| 成年女人在线观看亚洲视频| 亚洲婷婷狠狠爱综合网| 久久久午夜欧美精品| 久久久a久久爽久久v久久| 最近中文字幕高清免费大全6| 天天操日日干夜夜撸| av线在线观看网站| 中国国产av一级| 国产精品久久久久久久电影| 如何舔出高潮| 一级毛片电影观看| 国产有黄有色有爽视频| 欧美 日韩 精品 国产| 精品久久国产蜜桃| 最近中文字幕2019免费版| 一级av片app| 一级片'在线观看视频| 成人亚洲精品一区在线观看| 建设人人有责人人尽责人人享有的| 色吧在线观看| 日韩,欧美,国产一区二区三区| 国产黄色视频一区二区在线观看| 人人妻人人澡人人看| 免费看不卡的av| 日本欧美国产在线视频| 亚洲欧洲国产日韩| 两个人免费观看高清视频 | 伊人亚洲综合成人网| 新久久久久国产一级毛片| 久久99热这里只频精品6学生| 久久久亚洲精品成人影院| 老女人水多毛片| 免费观看无遮挡的男女| 日产精品乱码卡一卡2卡三| 久久久久久久大尺度免费视频| 欧美国产精品一级二级三级 | 麻豆精品久久久久久蜜桃| 又粗又硬又长又爽又黄的视频| .国产精品久久| 免费看日本二区| 成人特级av手机在线观看|