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

    基于遙感信息的呂梁山貧困區(qū)生態(tài)安全評價

    2019-12-26 03:39:28孫從建李曉明張文強王佳瑞
    中國環(huán)境科學 2019年12期
    關鍵詞:呂梁山連片變化

    孫從建,李曉明,張文強,陳 偉,王佳瑞

    基于遙感信息的呂梁山貧困區(qū)生態(tài)安全評價

    孫從建1*,李曉明1,張文強2,陳 偉1,王佳瑞1

    (1.山西師范大學地理科學學院,山西 臨汾 041000;2.中國科學院新疆生態(tài)與地理研究所,荒漠與綠洲生態(tài)國家重點實驗室,新疆 烏魯木齊 830011)

    本文基于四個時段的遙感影像光譜信息計算了呂梁山連片貧困區(qū)的遙感生態(tài)指數(shù)(RSEI)及植被覆蓋度(FVC),分析了其生態(tài)環(huán)境時空變化規(guī)律.結果表明:(1)1993~2018年間研究區(qū)FVC整體呈下降趨勢,FVC由1993年0.65下降到2018年0.55;而RSEI整體上亦呈下降趨勢,由1993年的0.47下降到2018年的0.40.(2)研究區(qū)整體生態(tài)環(huán)境呈退化趨勢,其中北部、南部退化最為嚴重,除黃河沿岸外,其他河道沿線區(qū)域的FVC及生態(tài)環(huán)境退化明顯.(3)研究區(qū)RSEI值整體偏低且區(qū)域差異性顯著,其中部地區(qū)是生態(tài)環(huán)境最為脆弱,北部忻州及南部臨汾地區(qū)近年來生態(tài)環(huán)境退化尤為顯著.(4)FVC變化趨勢與RSEI變化趨勢基本一致,FVC的增加對生態(tài)環(huán)境產(chǎn)生積極影響,而經(jīng)濟快速發(fā)展帶來的道路、建筑用地、農(nóng)田的增加對生態(tài)環(huán)境產(chǎn)生消極影響.區(qū)域扶貧發(fā)展的過程中應注意平衡經(jīng)濟發(fā)展與生態(tài)環(huán)境保護間的關系,避免進一步損害該區(qū)脆弱的生態(tài)環(huán)境.

    呂梁山連片貧困區(qū);遙感生態(tài)指數(shù);植被覆蓋度(FVC);定量評估

    呂梁山連片貧困區(qū)是我國14個集中連片貧困區(qū)之一,該區(qū)氣候干旱、災害頻發(fā)、水土流失嚴重、生態(tài)環(huán)境惡劣[1].為實現(xiàn)區(qū)域經(jīng)濟可持續(xù)發(fā)展,十分有必要對區(qū)域生態(tài)環(huán)境進行系統(tǒng)的認識和安全評價.

    現(xiàn)有的生態(tài)安全評價方法主要包括數(shù)學模型法[2]、生態(tài)模型法[3]、景觀格局法[4]、計算機評價法[5],這些方法在中西部銜接城市[2]、遼河三角洲濕地[4]、喀斯特地區(qū)[5]等地區(qū)的生態(tài)安全評價中得到了廣泛推廣應用.近年來隨著遙感技術的快速發(fā)展,其實時、快速、數(shù)據(jù)可靠等優(yōu)勢被研究者所肯定,并開始利用遙感監(jiān)測指標開展相應的生態(tài)安全評價[6-7].傳統(tǒng)的遙感生態(tài)安全評價方法大多基于單一監(jiān)測指標,具有計算簡易、易于解讀等優(yōu)點,但其難以較好的全面揭示生態(tài)環(huán)境的系統(tǒng)性變化[8-9].徐涵秋等[8-11]利用人類可直觀感知生態(tài)條件優(yōu)劣的植被指數(shù)、干度指數(shù)、地表溫度、濕度分量構建了一種新型遙感生態(tài)評價方法,該方法具有指標獲取簡易、無人為權重、結果可視化佳等優(yōu)勢,可以對區(qū)域生態(tài)環(huán)境的時空變化進行系統(tǒng)分析,亦可有效彌補單一指標法的不足[8-15].目前該方法已在我國東南丘陵水土流失區(qū)[9]、北方農(nóng)牧交錯帶[15]及雄安新區(qū)[11]等代表區(qū)域的生態(tài)安全評價中廣泛使用,并取得了顯著成效.

    呂梁山連片貧困區(qū)處于黃土高原腹地,生態(tài)環(huán)境十分脆弱,相關生態(tài)安全評價開展較少,而系統(tǒng)、全面、定量的評價區(qū)域生態(tài)環(huán)境的時空變化特征的工作在呂梁山區(qū)仍是空白[16].本研究基于多期遙感影像數(shù)據(jù)開展該區(qū)域生態(tài)環(huán)境安全的系統(tǒng)評價,所得結果對呂梁山區(qū)的災害防治、水土流失治理、生態(tài)環(huán)境建設等提供參考.

    1 研究區(qū)概況

    圖1 呂梁山連片貧困區(qū)示意

    研究區(qū)地處黃土高原東部,西以黃河為界,東北接大同盆地,東南為臨汾盆地,地理位置為110°22′~112°37′E,35°53′~39°24′N,面積約為3.3萬平方公里,涵蓋13個國家級貧困縣:神池縣、五寨縣、岢嵐縣、靜樂縣、興縣、嵐縣、臨縣、石樓縣、吉縣、大寧縣、隰縣、永和縣、汾西縣,7個臨近縣:柳林縣、蒲縣、交口縣、中陽縣、離石市、方山縣、寧武縣(圖1).呂梁山連片貧困區(qū)溝壑縱橫、地形切割、植被稀少、水土流失嚴重、土壤貧瘠;屬于半干旱大陸性季風氣候,年平均氣溫6~13℃,降水300~700mm,四季分明,雨熱同期,降水多集中在夏季.呂梁山連片貧困區(qū)城鎮(zhèn)化進程緩慢,就業(yè)率低,同時由于該區(qū)域資源匱乏、環(huán)境惡劣、交通不便、生產(chǎn)力水平低,導致社會經(jīng)濟發(fā)展能力不足,各產(chǎn)業(yè)發(fā)展受限,農(nóng)村勞動力大量流失,使得該區(qū)日益貧困[16].此外獨特的地理環(huán)境導致該區(qū)水土流失嚴重、生態(tài)系統(tǒng)脆弱,生態(tài)安全問題已成為該區(qū)域扶貧開發(fā)過程中的亟待解決問題.

    2 數(shù)據(jù)來源及方法

    2.1 遙感影像獲取

    本文選取了覆蓋研究區(qū)的3個時段的Landsat TM以及Landsat OLI遙感數(shù)據(jù)計算RSEI分別為:1993年9月15日(TM)、2000年9月2日(TM)、2009年9月27日(TM)、2018年9月12(OLI).數(shù)據(jù)獲取時間皆為初秋,保證了研究結果的可比性,亦可有效消除農(nóng)作物干擾誤差.數(shù)據(jù)預處理步驟包括鑲嵌、裁剪、輻射定標以及大氣校正等,具體步驟為:計算不同年份的綠度(NDVI)、熱度(LST)、濕度(Wet)、干度(NDSI)等4個指標、對4個指標進行歸一化處理、分別合成一幅影像、對新影像進行主成分轉換、得到遙感生態(tài)指數(shù)(RSEI)影像、計算植被覆蓋度(FVC).

    2.2 遙感生態(tài)指數(shù)計算

    本文利用人類可直觀感知生態(tài)條件優(yōu)劣的歸一化植被指數(shù)(代表NDVI)、濕度分量(代表Wet)、裸地指數(shù)(代表NDSI)、地表溫度(代表LST)來構建RSEI[8],實現(xiàn)區(qū)域生態(tài)環(huán)境的整體變化評價.

    (1)NDVI

    歸一化植被指數(shù)是檢測植被生長以及FVC的廣泛使用的指數(shù),為近紅外波段與紅外波段反射值之差與和的比,公式為:

    式中:nir、r分別對應近紅外、紅色波段的反射率.

    (2)Wet

    濕度采用纓帽變換中的Wet指數(shù)[8],該指數(shù)明確反映了植被生長過程中的水分變化,不同傳感器的公式不同.

    式中:B分別為對應波段的反射率.

    (3)NDSI[11]

    基于研究區(qū)裸地與建筑用地并存的條件,綜合裸地指數(shù)(SI)和建筑指數(shù)(IBI)形成干度指數(shù):

    式中:i為對應波段反射率.

    (4)LST

    計算研究區(qū)熱度采用了兩種算法,第一種基于影像的反演算法[13](IB算法),第二種輻射傳輸通道算法[14](RTE法),此兩種算法均可以利用遙感影像反演地表溫度.

    IB算法如下:

    式中:rad6為絕對亮溫(K);=11.5mm為熱紅外波段中心波長);=×/(1.438×10-2mK),其中普朗克常數(shù)=6.626×10-34J/K;光速=2.998×108m/s;玻爾茲曼常數(shù)e=1.38×10-23J/K;為地表比輻射率.

    RTE法如下:

    式中:t是衛(wèi)星高度上傳感器測得的輻射強度[W/(m2?sr?um)];為地表比輻射率;(S)為黑體熱輻射亮度;S為地表溫度(K);u-和dˉ分別是大氣上行下行的熱輻射強度;為大氣透射率,此3個指數(shù)均可通過NASA網(wǎng)站獲取.

    S可以根據(jù)普朗克公式函數(shù)計算:

    式中:1,2為常量,對于Landsat TM,1=607.76 [W/ (m2?sr?um)];2= 1260.56K;對于Landsat OLI,1= 774.89[W/(m2?sr?um)];2=1321.08K.

    由于以上4個指標的量綱不同,必須進行正規(guī)化處理后將量綱統(tǒng)一到[0,1]之間再進行主成分變換,各指標正規(guī)化公式如下:

    式中:NI為某一像元歸一化后的值,maxmin分別為該指標的最大值與最小值.

    南海爭端由來已久,爭端所涉及的各方利益錯綜復雜。幾十年來,南海爭端的起起伏伏表明,南海各方達成全面和永久性的解決非常困難,非一朝一夕可以完成。然而,南海的和平穩(wěn)定與繁榮發(fā)展事關沿岸數(shù)億人的生產(chǎn)生活,因此在全面和永久解決爭議之前,各方積極推進南海海洋合作,是維護海洋和平、增進地區(qū)共同福祉的有效途徑,也有利于培育各方合作互信的氛圍,為實現(xiàn)爭端的最終解決做好必要鋪墊。

    經(jīng)過正規(guī)化后的四個指標就可以進行主成分變換,得出RSEI值,然后再對RSEI進行正規(guī)化得出最終結果,歸一化后RSEI值越接近1表示生態(tài)環(huán)境越好.

    2.3 FVC計算

    為區(qū)分RSEI指數(shù)中NDVI對FVC的影響,采用三次梯度法計算FVC[17].三次梯度法是根據(jù)植被光譜特點,利用紅、綠、近紅外三個波段的反射率及波長計算得出,具體公式如下:

    式中:ir、r、g分別為近紅外、紅、綠波段的反射率;irrg分別為近紅外、紅、綠波段的中心波長;為梯度差指數(shù);FVC為植被覆蓋度.通過計算FVC對研究區(qū)地物進行分類,觀測植被的變化趨勢,進而分析生態(tài)環(huán)境的變化趨勢.

    2.4 RSEI、FVC時空變化計算

    首先利用ENVI采用Band math方法分別用后一年的RSEI、FVC值減去前一年對應的RSEI、FVC值:

    式中:C為對應年份計算結果,1為后一年對應指數(shù)所有波段,2為前一年對應指數(shù)所有波段.計算后得到的結果值在[-1,1]區(qū)間,為使結果便于觀察對比,將結果進行5等分,分別對應顯著退化[-1,-0.2)、輕微退化[-0.2,-0.05)、穩(wěn)定[-0.05,0.05)、輕微改善[0.05,0.2)、顯著改善[0.2,1]等5個等級.

    3 結果與分析

    3.1 呂梁山連片貧困區(qū)FVC變化特征

    由圖2可知,研究區(qū)的NDVI、FVC以及歸一化后FVC均呈先上升再下降趨勢,整體呈下降趨勢.其中FVC在1993~2000年從0.65上升為0.73,在2000~2009年下降了0.07,至2018年下降至0.55,過去25a來總體下降了0.1.FVC在1993~2000年整體呈上升趨勢,這與同期開展的黃土高原生態(tài)恢復措施有密切關系,從1990年開始,呂梁山區(qū)在“三北防護林工程”、“四荒”治理及黃土高原封山禁牧等生態(tài)工程的治理下區(qū)域FVC得到了有效的提高[18].但隨著2000年以來呂梁山區(qū)的城鎮(zhèn)化建設、開荒建地及礦產(chǎn)資源的大規(guī)模開發(fā),區(qū)域植被的恢復速度受到了抑制,同時逐漸硬化的下墊面更加劇了區(qū)域水土流失,導致區(qū)域生態(tài)環(huán)境日益惡化[19],這也引起3個植被指標在2000年以后均呈微弱下降趨勢.

    圖2 1993~2018年研究區(qū)FVC變化

    圖3b顯示,1993~2000年間研究區(qū)北部除寧武縣FVC降低外,其他縣FVC均升高,鄰近黃河的興縣FVC增加尤為明顯;中部的臨縣和方山縣大面積FVC較差等級轉為較好等級,其他縣變化不明顯;南部地區(qū)FVC整體變化不大.總體來看,FVC良好及優(yōu)秀等級的面積占比明顯增大,FVC較差的面積減小并且集中在南部黃河沿岸.這一變化趨勢與三北防護林建設工程實施有密切關系,該工程前三期(至2000年完成)將呂梁山中北部作為重點建設區(qū)域[21-22],在防護林體系建設下,中北部林地增加(圖3a),FVC顯著提高.在2000~2009年間,研究區(qū)中北部地區(qū)FVC出現(xiàn)下降趨勢,其中部分河道區(qū)域FVC降低較為明顯,這與同時期城鎮(zhèn)建設(多分布于地勢平坦的河谷,如圖3a所示)有一定關系.與北部和中部不同,南部地區(qū)各縣的FVC均呈上升趨勢,其中以黃河沿岸區(qū)域尤為顯著,這與2001年開始的三北防護林建設四期工程實施的晉陜峽谷防護林體系建設項目有密切關系[22].在2009~2018年間,研究區(qū)大部分區(qū)域FVC呈現(xiàn)降低趨勢,其中北部、南部FVC差等級的面積較為集中.近年來,在人口增長、經(jīng)濟發(fā)展的壓力下,濫墾、濫伐、濫牧、濫挖草坯現(xiàn)象更加頻繁,林木資源受到嚴重破壞,森林植被減少明顯,農(nóng)田、建筑面積快速增加(圖3a),這是2009~2018年間FVC降低的重要原因.

    分析呂梁山連片貧困區(qū)四時段各級FVC面積比重變化可知(表1所示),在1993~2000年間,FVC中等及以下等級的面積比重呈減小趨勢(負值),FVC較高以上等級的面積比重增加明顯,增幅為12.33%.2000~2009年,FVC低等級和FVC較高等級區(qū)域變化不明顯,但是中等的區(qū)域增加了8.31%, FVC高等級區(qū)域減少了10.65%,說明這一時期區(qū)域植被破壞較為嚴重.2009~2018年間,FVC中等及以下等級的面積比重呈增加趨勢,較高及高FVC面積比重呈下降趨勢,表明該區(qū)域FVC降低趨勢在這一時期仍未改善.

    3.2 呂梁連片貧困區(qū)生態(tài)環(huán)境變化分析

    從表2可知,4個時段內(nèi),前兩個成分PC1與PC2累計貢獻率均達88%以上,其PC1貢獻率分別為82.81%、83.05%、77.22%、71.85%,表明主成分1已經(jīng)集成了4個指標的大部分特征值,同時不同時期4個指數(shù)對主成分的貢獻率相對穩(wěn)定.觀察4個指數(shù)在主成分1中的值,發(fā)現(xiàn)干度與熱度值的符號相同且與綠度及濕度值的符號正好相反,說明這兩對指數(shù)在生態(tài)環(huán)境指數(shù)中的作用是相反的,其中綠度與濕度對生態(tài)環(huán)境起積極作用,干度與熱度指數(shù)起消極作用,這與實際情況相對應.除2009年外的其余時間段,不同指數(shù)的PC1值均表現(xiàn)為:熱度與干度為正值,濕度與綠度為負值.為了使RSEI研究圖結果一致,(值越接近1生態(tài)環(huán)境越好),對這3a進行反向處理,用1減去這3a的值,則PC1結果越大表明生態(tài)環(huán)境越好.

    圖3 研究區(qū)各年LUCC及FVC、RSEI、NDVI、Wet、LST和NDSI各等級區(qū)域時空變化

    表1 不同等級FVC面積比重表

    表2 研究區(qū)遙感生態(tài)指數(shù)主成分分析表

    由表3可知:(1)25a來濕度從0.57降低到0.48,除1993年到2000年升高外,濕度持續(xù)下降; (2)1993~2009年綠度增大,在2009~2018年下降,但整體呈下降趨勢;(3)干度1993~2018年先下降再上升,整體上升了0.1,說明研究區(qū)在這一時段研究區(qū)有變干趨勢;(4)熱度指數(shù)較為穩(wěn)定,表明區(qū)域溫度變化不顯著.總體而言,1993~2018年對生態(tài)環(huán)境起積極作用的綠度、濕度指標降低,而對生態(tài)環(huán)境指數(shù)起消極作用的干度、熱度指標升高,表明呂梁山連片貧困區(qū)的FVC降低,地面硬化趨勢增加,且遙感生態(tài)指數(shù)從0.47下降到0.40,整體下降明顯,說明研究區(qū)生態(tài)環(huán)境變差.

    進一步分析在研究時段內(nèi)4個指標不同等級區(qū)域的時空變化特征可知(圖3d、e、f、g),在1993~2018年間,NDVI指數(shù)中等以上等級面積整體呈增加趨勢,尤其以研究區(qū)北部和黃河沿岸最為顯著(圖3d).Wet指數(shù)變差的區(qū)域主要集中于北部地區(qū),其他區(qū)域基本保持穩(wěn)定(圖3e).盡管在過去25a間,LST指數(shù)整體均值并未明顯波動,但其不同等級區(qū)域空間變化卻較顯著,其中北部地區(qū)LST指數(shù)低等級區(qū)域呈增加趨勢,而南部地區(qū)及黃河沿岸LST指數(shù)高等級區(qū)域增長顯著(圖3f).研究區(qū)NDSI指數(shù)較高以上等級區(qū)域增加顯著,其中北部地區(qū)基本全部變?yōu)楦叩燃墔^(qū)域.

    表3 研究區(qū)各時段主要指數(shù)及RSEI均值變化表

    基于研究區(qū)遙感生態(tài)指數(shù)(RSEI)不同等級區(qū)域時空分布圖(圖3c)可知,在1993~2000年間研究區(qū)RSEI整體呈上升趨勢,其中北部地區(qū)神池縣、寧武縣、岢嵐縣生態(tài)環(huán)境趨于好轉,其西部地區(qū)由優(yōu)秀等級降低為中等或較差等級,嵐縣南部由良好降低為較差.中部地區(qū)RSEI由差、較差轉為中等、較好等級,柳林、離石、中陽縣所在區(qū)域的三川河河道環(huán)境轉好,這與三北防護林前三期工程的實施具有密切聯(lián)系.研究區(qū)南部地區(qū)隰縣的RSEI變化較為顯著,其較差面積擴大.

    對比2000年與2009年發(fā)現(xiàn),研究區(qū)生態(tài)環(huán)境明顯變差,特別是中北部、中部及南部河谷區(qū)域,包括岢嵐縣所在的漪嵐河河道,方山縣、柳林縣、離石區(qū)、中陽縣區(qū)域內(nèi)的三川河、昕水河及其支流所在河道,以及隰縣、交口縣境內(nèi)河谷地區(qū).上述區(qū)域地勢相對平坦、適于開發(fā)居住,由于城鎮(zhèn)化的發(fā)展,沿河布局了大量道路、建筑物及新開農(nóng)田(圖3a),其NDSI指數(shù)顯著升高(表3)導致生態(tài)環(huán)境變差.

    在2009~2018年間,研究區(qū)生態(tài)環(huán)境呈現(xiàn)進一步變差趨勢,2009年分布于中北部岢嵐縣、靜樂縣等縣境內(nèi)河谷地區(qū)的RSEI差等級的區(qū)域已擴大為面狀區(qū)域,上述區(qū)域的城鎮(zhèn)化進程及新開農(nóng)田(圖3a)已經(jīng)嚴重威脅區(qū)域生態(tài)安全.此外,研究區(qū)東北部寧武縣、靜樂縣、嵐縣的RSEI優(yōu)秀、良好等級的面積急劇減小,而較差、差的面積顯著增加;興縣、吉縣、汾西縣的生態(tài)環(huán)境也趨于退化,這些變化均說明2009~2018年呂梁山區(qū)的生態(tài)環(huán)境呈退化趨勢.

    3.3 RSEI與FVC的時空變化分析

    圖4 呂梁山連片貧困區(qū)RSEI與FVC時空變化

    由圖4可知,1993~2000年研究區(qū)大部分地區(qū)生態(tài)環(huán)境得到改善,生態(tài)退化區(qū)域主要分布于東北部地區(qū)和南部臨黃河地區(qū).二十世紀90年代山西省開啟“四荒”治理、三北防護林建設是這一時段生態(tài)環(huán)境變好的主要原因;而黃河沿岸在這一時段由于盲目大規(guī)模旅游開發(fā)帶來的地面硬化、水土流失加劇,是該區(qū)域生態(tài)環(huán)境惡化主要原因.在2000~2009年,研究區(qū)南部地區(qū)及黃河沿岸生態(tài)環(huán)境得到改善,2001年開始的三北防護林四期工程將黃河沿岸定為重點治理區(qū)域,有效提高了黃河中游河道附近的FVC[23],并提高了該區(qū)域的生態(tài)安全指數(shù).但是中北部地區(qū)環(huán)境在這一時段顯著變差,人類活動可能是導致該區(qū)域生態(tài)環(huán)境變化的主要因素.在2009~2018年間,研究區(qū)中部環(huán)境得到改善,但呂梁山區(qū)北部、南部地區(qū)生態(tài)環(huán)境卻呈退化趨勢,這可能與上述區(qū)域城鎮(zhèn)化進程及經(jīng)濟快速發(fā)展有一定關系,北部忻州及南部臨汾地區(qū)相較于研究區(qū)其他區(qū)域發(fā)展速度較快,同時環(huán)境退化也較為明顯[20].對比RSEI與FVC時空變化圖發(fā)現(xiàn),不同時間段的RSEI變化與FVC變化區(qū)域基本重合,植被退化區(qū)則為RSEI退化區(qū),植被恢復區(qū)的RSEI也顯著升高,說明植被對生態(tài)環(huán)境的變化具有顯著的影響.

    整體上,呂梁山連片貧困區(qū)的RSEI值一直偏低,這與區(qū)域地形破碎、自然環(huán)境惡劣、水土流失嚴重有一定關系[20-21].植被恢復措施(如三北防護林工程)在呂梁山連片貧困區(qū)的實施[22]對研究區(qū)生態(tài)環(huán)境保護有積極的促進作用,但經(jīng)濟發(fā)展帶來的城鎮(zhèn)化建設、開墾種植活動嚴重威脅了區(qū)域生態(tài)安全.在區(qū)域脫貧攻堅過程中,應平衡經(jīng)濟發(fā)展與生態(tài)環(huán)境保護的關系,對于生態(tài)退化顯著的忻州地區(qū)及臨汾地區(qū)應加強生態(tài)治理投入,限制無需開墾農(nóng)地、推進退耕還林還草工程;此外對于分布于河谷地區(qū)的鄉(xiāng)鎮(zhèn)應合理規(guī)劃其城鎮(zhèn)化建設,加大對小城鎮(zhèn)生態(tài)環(huán)境建設的投入,改善河谷地區(qū)的生態(tài)環(huán)境.

    4 結論

    4.1 在1993~2018年間,呂梁山連片貧困區(qū)在高原生態(tài)恢復工程及經(jīng)濟快速發(fā)展的聯(lián)合影響下,FVC呈現(xiàn)先升后降,整體下降的趨勢,而RSEI呈現(xiàn)下降趨勢,由1993年的0.47下降為2018年的0.40,表明研究區(qū)整體生態(tài)環(huán)境呈惡化趨勢.

    4.2 呂梁山區(qū)RSEI整體偏低,區(qū)域差異性顯著,其中自然環(huán)境惡劣的中部地區(qū)RSEI波動最顯著,生態(tài)環(huán)境最脆弱;受經(jīng)濟快速發(fā)展影響,北部忻州及南部臨汾地區(qū)近年來生態(tài)環(huán)境退化尤為顯著.

    4.3 研究區(qū)FVC與RSEI變化趨勢相近,關系密切,FVC增高的區(qū)域生態(tài)環(huán)境趨于好轉,FVC降低的區(qū)域生態(tài)環(huán)境惡化.

    4.4 國家生態(tài)恢復工程的布局建設與區(qū)域經(jīng)濟發(fā)展是影響區(qū)域生態(tài)環(huán)境的兩個重要因素.三北防護林工程建設帶來的森林、草地面積的增加對研究區(qū)北部及黃河沿岸的生態(tài)環(huán)境產(chǎn)生積極影響;而區(qū)域經(jīng)濟的快速發(fā)展帶來道路、建筑用地、農(nóng)田的增加對研究區(qū)內(nèi)部河谷地區(qū)及臨汾忻州等地的生態(tài)環(huán)境產(chǎn)生消極影響.區(qū)域扶貧發(fā)展的過程中要注意平衡經(jīng)濟發(fā)展與生態(tài)環(huán)境保護間的關系,限制大規(guī)模的農(nóng)地開墾及河谷區(qū)城鎮(zhèn)建設,加強忻州、臨汾等生態(tài)脆弱區(qū)、河谷地區(qū)的生態(tài)建設投入,避免進一步損害該區(qū)脆弱的生態(tài)環(huán)境.

    [1] 孫從建,張文強,李新功,等.基于遙感影像的黃土高原溝壑區(qū)生態(tài)效應評價[J]. 農(nóng)業(yè)工程學報, 2019,35(12):165-172. Sun C J, Zhang W Q, Li X G, et al. Evaluation of ecological effect of gully region of loess plateau based onremote sensing image [J]. Transactions of The Chinese Society of Agricultural Machinery, 2019, 35(12):165-172.

    [2] 黃 海,劉長城,陳 春.基于生態(tài)足跡的土地生態(tài)安全評價研究[J]. 水土保持研究, 2013,20(1):193-196+201. Huang H, Liu C C, Chen C. Appraisal of land ecological security based on ecological footprint [J]. Research of Soil and Water Conservation, 2013,20(1):193-196+201.

    [3] 韓振華,李建東,殷 紅,等.基于景觀格局的遼河三角洲濕地生態(tài)安全分析[J]. 生態(tài)環(huán)境學報, 2010,19(3):701-705. Han Z H, Li J D, Yin H, et al. Analysis of ecological security of wetland in Liaohe River delta based on the landscape pattern [J]. Ecology and Environment, 2010,19(3):701-705.

    [4] 惠 麗,魯小珍,張大勇,等.基于生態(tài)足跡法縣級區(qū)域可持續(xù)發(fā)展趨勢探討[J]. 中國人口·資源與環(huán)境, 2010,20(S1):486-488. Hui L, Lu X Z, Zh D Y, et al. Ecological footprint study on sustainable development [J]. China Population, Resources and Environment, 2010,20(S1):486-488.

    [5] 魏小島,周忠發(fā),王媛媛.基于格網(wǎng)GIS的喀斯特生態(tài)安全研究——以貴州花江石漠化綜合治理示范區(qū)為例[J]. 山地學報, 2012,30(6):681-687. Wei X D, Zhou Z F, Wang Y Y. Research on the Karst ecological security based on gridding GIS——acase study on the demonstration areas of the rock desertification integrated management in Huajiang of Guizhou, China [J]. Journal of Mountain Science, 2012,30(6):681- 687.

    [6] 孟 巖,趙庚星.基于衛(wèi)星遙感數(shù)據(jù)的河口區(qū)生態(tài)環(huán)境狀況評價——以黃河三角洲墾利縣為例 [J]. 中國環(huán)境科學, 2009,29(2):163- 167. Meng Y, Zhao G X, Ecological environment condition evaluation of estuarine area based on quantitative remote sensing —a case study in Kenli County [J]. China Environmental Science, 2009,29(2):163-167.

    [7] 史紀安,劉玉華,師江瀾,等.江河源區(qū)生態(tài)環(huán)境質量綜合評價[J]. 西北農(nóng)林科技大學學報(自然科學版), 2006(10):61-66. Shi J A, Liu Y H, Shi J L, et al.Synthetical evaluation of eco- environmental quality in the headstream of Yellow, Yangtze and Lantsan Rivers [J]. Journal of Northwest A & F University (Natural Science Edition), 2006,(10):61-66.

    [8] 徐涵秋.城市遙感生態(tài)指數(shù)的創(chuàng)建及其應用[J]. 生態(tài)學報, 2013, 33(24):7853-7862. Xu H Q. A remote sensing urban ecological index and its application [J]. Acta Ecologica Sinica., 2013,33(24):7853-7862.

    [9] 張 燦,徐涵秋,張 好,等.南方紅壤典型水土流失區(qū)FVC變化及其生態(tài)應評估——以福建省長汀縣為例[J]. 自然資源學報, 2015, 30(6):917-928. Zhang C, Xu H Q, Zhang H, et al. Fractional vegetation cover change and its ecological effectassessment in a typical reddish soil region of southeastern China: Changtingcounty, Fujian province [J]. Journal of Natural Resources, 2015,30(6):917-928.

    [10] 徐涵秋,施婷婷,王美雅,等.雄安新區(qū)地表覆蓋變化及其新區(qū)規(guī)劃的生態(tài)響應預測[J]. 生態(tài)學報, 2017,37(19):6289-6301. Xu H Q, Shi T T, Wang M Y, et al. Prediction of land cover change and ecological response of Xiongan new area planning [J]. ActaEcologica Sinica, 2017,37(19):6289-6301.

    [11] 徐涵秋.水土流失區(qū)生態(tài)變化的遙感評估[J]. 農(nóng)業(yè)工程學報, 2013,29(7):91-97. Xu H Q. Assessment of ecological change in soil loss area using remote sensing technology [J]. Transactions of the Chinese Society of Agricultural Engineering., 2013,29(7):91-97.

    [12] 徐涵秋.區(qū)域生態(tài)環(huán)境變化的遙感評價指數(shù)[J]. 中國環(huán)境科學, 2013,33(5):889-897. Xu H Q. A remote sensing index for assessment of regional ecological changes [J]. China Environmental Science, 2013,33(5):889-897.

    [13] 丁 鳳,徐涵秋.基于Landsat TM的3種地表溫度反演算法比較分析 [J]. 福建師范大學學報(自然科學版), 2008,(1):91-96. Ding F, Xu H Q. Comparison of three algorithms for retrieving land surface temperature from Landsat TM thermal infrared band [J]. Journal of Fujian Normal University (natural science edition), 2008, (1):91-96.

    [14] 徐涵秋.新型Landsat8衛(wèi)星影像的反射率和地表溫度反演[J]. 地球物理學報, 2015,58(3):741-747. Xu H Q. Retrieval of the reflectance and land surface temperature of the newly-launched Landsat8satellite [J]. Chinese Journal of Geophysics, 2015,58(3):741-747.

    [15] 張曉東,劉湘南,趙志鵬,等.農(nóng)牧交錯區(qū)生態(tài)環(huán)境質量遙感動態(tài)監(jiān)測——以寧夏鹽池為例[J]. 干旱區(qū)地理, 2017,40(5):1070-1078. Zhang X D, Liu X N, Zhao Z P, et al. Dynamic monitoring of ecology and environment in the agro-pastralecotone based on remote sensing:a case of Yanchi County in Ningxia Hui Autonomous Region [J].Arid Land Geography, 2017,40(5):1070-1078.

    [16] 徐滿厚,楊曉艷,張瀟月,等.山西呂梁連片特困區(qū)生態(tài)環(huán)境質量評價及其經(jīng)濟貧困的時空分布特征[J]. 江蘇農(nóng)業(yè)科學, 2018,46(6): 304-309. Xu M H, Yang X Y, Zhang X Y, et al. Evaluation of ecological environment quality and spatial-temporal distribution characteristics of economic poverty in the contiguous poor areas of Lüliang, Shanxi [J]. Jiangsu Agricultural Science, 2008,46(6):304-309.

    [17] 唐世浩,朱啟疆,王錦地,等.三波段梯度差植被指數(shù)的理論基礎及其應用[J]. 中國科學(D輯:地球科學), 2003,(11):1094-1102. Tang S H, Zhu Q J, Wang J D, et al. Theoretical basis and application of three-band gradient vegetation index [J]. Science in China (series D: earth sciences), 2003,(11):1094-1102.

    [18] 趙國嶺.管好“四荒”拍賣金投入治理再發(fā)展——嵐縣“四荒”拍賣資金管理使用辦法調查[J]. 山西水利, 1995,(1):29. Zhao G L. Managing the "four wastes" auction fund management and development -- investigation on the management and use of auction fund in LAN county "four wastes" [J]. Shanxi Water Resources, 1995, (1):29.

    [19] 韓東娥.實施“生態(tài)興省”戰(zhàn)略建設綠色低碳能源基地[J]. 科技創(chuàng)新與生產(chǎn)力, 2010,(5):21-24. Han D E. Implementing the strategy of "promoting the province by ecology" to build green and low-carbon energy base [J]. Sci-TechInnovation and Productivity, 2010,(5):21-24.

    [20] 冀 欽,楊建平,徐滿厚.山西呂梁山連片特困區(qū)現(xiàn)代農(nóng)業(yè)發(fā)展水平綜合評價[J]. 中國人口·資源與環(huán)境, 2018,28(S1):54-59. Ji Q, Yang J P, Xu M H. Comprehensive evaluation of modern agricultural development level in the continuous poor areas of Lüliang Mountains, Shanxi Province [J]. China Population, Resources and Environment, 2008,28(S1):54-59.

    [21] 閆珅愛.控制水土流失改善生態(tài)環(huán)境是發(fā)展呂梁農(nóng)村經(jīng)濟的關鍵[J]. 山西農(nóng)經(jīng), 2003,(5):44-46. Yan S A. Controlling soil erosion and improving ecological environment is the key to developing rural economy in Lüliang [J]. Shanxi Agricultural Economics, 2003,(5):44-46.

    [22] 任連功.呂梁:“三北”防護林工程的創(chuàng)新與發(fā)展[J]. 中國林業(yè), 2010,(9):35. Ren L G. Lüliang: innovation and development of "three north" shelterbelt project [J]. Forestry of China, 2010,(9):35.

    [23] 周 旭,楊勝天,劉曉燕,等.黃河中游多沙粗沙區(qū)流域坡面水保措施變化特征[J]. 地理學報, 2014,69(1):64-72. Zhou X, Yang S T, Liu X Y, et al. Characteristics of soil and water conservation measures change in thewatershed slope of the middle Yellow River Basin [J]. ActaGeographica Sinica, 2014,69(1):64-72.

    Evaluation of ecological security in poverty-stricken region of Lüliang Mountain based on the remote sensing image.

    SUN Cong-jian1*, LI Xiao-ming1, ZHANG Wen-qiang2, CHEN Wei1, WANG Jia-rui1

    (1.School of Geographical Sciences, Shanxi Normal University, Linfen 041000, China;2.State Key Laboratory of Desert and Oasis Ecology, Xinjiang Institute of Ecology and Geography, Chinese Academy of Science, Urumqi 830011, China)., 2019,39(12):5352~5360

    Based on the spectral information of remote sensing images during four periods, the remote sensing ecological index (RSEI) and fractional vegetation cover (FVC) were calculated for the poverty-stricken region of the Lüliang Mountains (PSLM), the temporal and spatial variation of regional ecological environment was subsequently evaluated. The results showed that: (1) The FVC in the research area decreased from 0.65 in 1993 to 0.55 in 2018, and the RSEI decreased from 0.47 in 1993 to 0.40 in 2018. (2) The overall ecological environment of the PSLM showed a trend of degradation during the research period, with the most serious degradation happened in the north and south part of the PSLM. Except for the region of Yellow River valley, all the other river valleys experienced the significant degradation in FVC and ecological environment. (3) The RSEI value in the PSLM region was generally lower and showed substantial spatial variations. The ecological environment in the central part of the study area was the most vulnerable and the part of northern Xinzhou and southern Linfen area was deteriorated particularly significantly. (4) The variation pattern of FVC was generally consistent with that of RSEI. The improvement of FVC has a positive influence on the local ecological environment, while the rapid economic development (increasing in roads, construction land and farmland) has negative impacts on the local ecological environment. In the process of regional poverty alleviation and development, more attention should be paid to balance the needs of regional economic development and ecological environment protection in order to avoid further damage to the fragile regional ecological environment.

    poverty-stricken region of the Lüliang Mountain;remote sensing ecological index;fractional vegetation cover;quantitative evaluation

    X826

    A

    1000-6923(2019)12-5352-09

    孫從建(1986-),男,河北滄州人.博士,副教授,主要從事自然地理及自然資源學的研究.發(fā)表論文42篇.

    2019-05-29

    山西省哲社規(guī)劃課題(2019B202);山西省研究生教學改革課題(2019JG123)

    * 責任作者, 副教授, suncongjian@sina.com

    猜你喜歡
    呂梁山連片變化
    蒲谷香:從呂梁山走到巴拿馬
    呂梁山選煤廠煤泥水處理工藝優(yōu)化改造
    選煤技術(2021年3期)2021-10-13 07:33:34
    呂梁山林區(qū)白皮松育苗技術
    呂梁山國有林管理局 強化創(chuàng)新驅動奮力蹚出高質量發(fā)展新路徑
    從9到3的變化
    學生天地(2020年14期)2020-08-25 09:20:58
    大別山連片特困地區(qū)反貧困綜合績效模糊評價
    智富時代(2019年2期)2019-04-18 07:44:42
    這五年的變化
    特別文摘(2018年3期)2018-08-08 11:19:42
    應用匹配追蹤傅里葉插值技術實現(xiàn)OVT域連片處理
    湖北省鄉(xiāng)鎮(zhèn)連片開發(fā)工程調查
    鳥的變化系列
    詩選刊(2015年6期)2015-10-26 09:47:11
    哪里可以看免费的av片| 国产日本99.免费观看| 美女高潮的动态| 国内精品久久久久久久电影| 亚洲成人久久爱视频| 午夜精品久久久久久毛片777| 亚洲成人久久爱视频| 午夜激情欧美在线| 亚洲最大成人手机在线| 亚洲国产色片| 色视频www国产| 精品久久久久久久久久久久久| 免费看光身美女| 99国产极品粉嫩在线观看| 欧美精品国产亚洲| 亚洲国产色片| 日韩欧美在线乱码| 国产高潮美女av| 美女高潮喷水抽搐中文字幕| 亚洲欧美日韩高清专用| 国产精品野战在线观看| 日韩精品青青久久久久久| 国产成人福利小说| 国内精品久久久久久久电影| 亚洲午夜理论影院| 99国产极品粉嫩在线观看| 国产亚洲欧美在线一区二区| 很黄的视频免费| 可以在线观看毛片的网站| 麻豆成人av在线观看| av欧美777| 天美传媒精品一区二区| 亚洲男人的天堂狠狠| 欧美日韩瑟瑟在线播放| 久久这里只有精品中国| 最后的刺客免费高清国语| 内地一区二区视频在线| 熟女电影av网| 国产高清视频在线观看网站| 精品熟女少妇八av免费久了| 中文字幕高清在线视频| 精品久久久久久,| 亚洲最大成人av| 亚洲人与动物交配视频| 久久久久久久久中文| 日韩人妻高清精品专区| 亚洲最大成人av| 久久九九热精品免费| 在线免费观看不下载黄p国产 | 91九色精品人成在线观看| av天堂中文字幕网| 三级男女做爰猛烈吃奶摸视频| 大型黄色视频在线免费观看| 国产91精品成人一区二区三区| 国产精品日韩av在线免费观看| 免费大片18禁| 久久午夜亚洲精品久久| 精品久久久久久成人av| 国内精品久久久久精免费| 欧美精品啪啪一区二区三区| 永久网站在线| 国产一区二区在线观看日韩| 听说在线观看完整版免费高清| 中国美女看黄片| 俺也久久电影网| 一本久久中文字幕| 天天躁日日操中文字幕| 国产精品久久久久久人妻精品电影| 麻豆成人午夜福利视频| 色精品久久人妻99蜜桃| 成人午夜高清在线视频| 久久6这里有精品| 一进一出抽搐动态| 精品久久久久久成人av| 日日摸夜夜添夜夜添小说| 少妇熟女aⅴ在线视频| 久久精品人妻少妇| 日韩欧美国产一区二区入口| 神马国产精品三级电影在线观看| 午夜福利在线在线| 九九久久精品国产亚洲av麻豆| 九色成人免费人妻av| 亚洲五月天丁香| 99热这里只有是精品在线观看 | 欧美又色又爽又黄视频| 午夜福利在线在线| 少妇的逼好多水| 一区二区三区免费毛片| 欧美丝袜亚洲另类 | 国产欧美日韩精品一区二区| 免费av观看视频| 国语自产精品视频在线第100页| 身体一侧抽搐| 日韩欧美 国产精品| 久久国产精品影院| 看黄色毛片网站| 在线播放无遮挡| 欧美日本视频| 国产一级毛片七仙女欲春2| 日本五十路高清| 99久久久亚洲精品蜜臀av| 国产白丝娇喘喷水9色精品| av天堂中文字幕网| 色哟哟哟哟哟哟| 欧美在线黄色| 久久中文看片网| 两个人视频免费观看高清| 俺也久久电影网| 极品教师在线视频| 欧美最黄视频在线播放免费| 老熟妇乱子伦视频在线观看| 中文字幕av成人在线电影| 久久国产乱子免费精品| 美女高潮喷水抽搐中文字幕| 日本a在线网址| 给我免费播放毛片高清在线观看| 亚洲人与动物交配视频| 亚洲avbb在线观看| 国产激情偷乱视频一区二区| 国产白丝娇喘喷水9色精品| 日本黄大片高清| 亚洲精品成人久久久久久| 精品午夜福利在线看| 一夜夜www| 最近在线观看免费完整版| 久9热在线精品视频| 国产成人av教育| 欧美国产日韩亚洲一区| 亚洲精品乱码久久久v下载方式| 亚洲av不卡在线观看| 亚洲av电影不卡..在线观看| 亚洲avbb在线观看| 亚洲精品乱码久久久v下载方式| 国产久久久一区二区三区| 日本一本二区三区精品| av在线老鸭窝| 午夜精品在线福利| 亚洲人成伊人成综合网2020| 欧美绝顶高潮抽搐喷水| 听说在线观看完整版免费高清| 国产精品亚洲av一区麻豆| 一夜夜www| 琪琪午夜伦伦电影理论片6080| 老鸭窝网址在线观看| 成熟少妇高潮喷水视频| 亚洲 欧美 日韩 在线 免费| 床上黄色一级片| 十八禁人妻一区二区| 亚洲天堂国产精品一区在线| 波野结衣二区三区在线| 色吧在线观看| 一级作爱视频免费观看| 老司机深夜福利视频在线观看| 少妇被粗大猛烈的视频| 久久性视频一级片| 一本综合久久免费| 欧美午夜高清在线| 一个人观看的视频www高清免费观看| 中亚洲国语对白在线视频| 一本久久中文字幕| 成人av在线播放网站| 淫秽高清视频在线观看| 99视频精品全部免费 在线| 亚洲专区国产一区二区| 在线观看免费视频日本深夜| 91久久精品电影网| 91字幕亚洲| 99在线人妻在线中文字幕| 中亚洲国语对白在线视频| 蜜桃亚洲精品一区二区三区| 怎么达到女性高潮| 别揉我奶头 嗯啊视频| 国产亚洲av嫩草精品影院| 久久草成人影院| 国产精品电影一区二区三区| 窝窝影院91人妻| 亚洲三级黄色毛片| 啦啦啦韩国在线观看视频| 好看av亚洲va欧美ⅴa在| 日韩中文字幕欧美一区二区| 五月伊人婷婷丁香| 免费高清视频大片| 精品久久久久久久末码| 精品人妻熟女av久视频| 亚洲人与动物交配视频| 亚洲成av人片免费观看| 在线观看午夜福利视频| 美女cb高潮喷水在线观看| 男人狂女人下面高潮的视频| 波多野结衣巨乳人妻| 色哟哟哟哟哟哟| 日本 av在线| 国产精品精品国产色婷婷| 日本一二三区视频观看| 欧美3d第一页| 免费人成在线观看视频色| 亚洲国产精品合色在线| 麻豆久久精品国产亚洲av| 国产蜜桃级精品一区二区三区| 国产精品爽爽va在线观看网站| 我要看日韩黄色一级片| 亚洲国产色片| 成人亚洲精品av一区二区| 久久久久精品国产欧美久久久| av在线老鸭窝| 久久久久久久久久成人| 免费在线观看影片大全网站| 亚洲国产精品sss在线观看| 日韩欧美一区二区三区在线观看| 宅男免费午夜| 久久久久久国产a免费观看| 乱人视频在线观看| 国产黄色小视频在线观看| 老司机午夜十八禁免费视频| 国产精品99久久久久久久久| 欧美一区二区精品小视频在线| 97碰自拍视频| 日韩大尺度精品在线看网址| 成人欧美大片| 久久久久精品国产欧美久久久| 欧美一级a爱片免费观看看| 免费无遮挡裸体视频| 国产爱豆传媒在线观看| 精品久久久久久,| 久久国产乱子伦精品免费另类| 少妇人妻一区二区三区视频| 亚洲欧美日韩无卡精品| 亚州av有码| 亚洲国产欧洲综合997久久,| 变态另类丝袜制服| 国产精品精品国产色婷婷| 90打野战视频偷拍视频| 听说在线观看完整版免费高清| 久久久久国内视频| 久久精品国产亚洲av涩爱 | 国产精品乱码一区二三区的特点| 国产成年人精品一区二区| 成年女人毛片免费观看观看9| 动漫黄色视频在线观看| 欧美在线一区亚洲| 日日夜夜操网爽| 国产不卡一卡二| 能在线免费观看的黄片| 国产白丝娇喘喷水9色精品| 国产一区二区在线av高清观看| 人妻久久中文字幕网| 两个人的视频大全免费| 老鸭窝网址在线观看| 午夜日韩欧美国产| 亚洲av美国av| 国产久久久一区二区三区| 欧美+亚洲+日韩+国产| 无遮挡黄片免费观看| 最近在线观看免费完整版| 久久久久国产精品人妻aⅴ院| 国产一区二区在线av高清观看| 直男gayav资源| 亚洲真实伦在线观看| 男人的好看免费观看在线视频| 国产成+人综合+亚洲专区| 男女床上黄色一级片免费看| 欧美国产日韩亚洲一区| 最新中文字幕久久久久| 成人鲁丝片一二三区免费| 国内精品美女久久久久久| 岛国在线免费视频观看| 少妇的逼水好多| 国产精品一及| 国产久久久一区二区三区| 国产野战对白在线观看| 精品熟女少妇八av免费久了| 人妻制服诱惑在线中文字幕| 国产亚洲欧美在线一区二区| 淫妇啪啪啪对白视频| 色噜噜av男人的天堂激情| 亚洲人与动物交配视频| 国产白丝娇喘喷水9色精品| 丁香欧美五月| 两性午夜刺激爽爽歪歪视频在线观看| 久久热精品热| 国产精品三级大全| 麻豆av噜噜一区二区三区| 熟妇人妻久久中文字幕3abv| 日韩欧美精品v在线| 制服丝袜大香蕉在线| 欧美3d第一页| 日韩国内少妇激情av| 亚洲国产精品999在线| 欧美黑人巨大hd| 精品久久久久久久久久久久久| 在线天堂最新版资源| 少妇丰满av| 99久久无色码亚洲精品果冻| 亚洲av电影不卡..在线观看| 国产精品亚洲美女久久久| 国产野战对白在线观看| 国产成人av教育| 国内精品美女久久久久久| 色精品久久人妻99蜜桃| 白带黄色成豆腐渣| 亚洲av熟女| 在线国产一区二区在线| 成人特级av手机在线观看| 又黄又爽又刺激的免费视频.| 午夜a级毛片| 国产单亲对白刺激| 大型黄色视频在线免费观看| 久99久视频精品免费| 国产极品精品免费视频能看的| 天天躁日日操中文字幕| 亚洲色图av天堂| 日本a在线网址| 中文资源天堂在线| 嫩草影院精品99| 怎么达到女性高潮| 一二三四社区在线视频社区8| 日韩欧美国产在线观看| eeuss影院久久| 国产单亲对白刺激| 亚洲乱码一区二区免费版| 熟女人妻精品中文字幕| 两个人的视频大全免费| 别揉我奶头~嗯~啊~动态视频| 欧美日韩综合久久久久久 | 亚洲av美国av| 亚洲在线自拍视频| 国产精品一及| av在线老鸭窝| 精品午夜福利在线看| 在线播放国产精品三级| 国产蜜桃级精品一区二区三区| 老司机午夜福利在线观看视频| 神马国产精品三级电影在线观看| 少妇人妻精品综合一区二区 | 中亚洲国语对白在线视频| 又黄又爽又免费观看的视频| 亚洲欧美日韩卡通动漫| 看十八女毛片水多多多| 一个人观看的视频www高清免费观看| 99国产精品一区二区三区| 18禁在线播放成人免费| 国产伦在线观看视频一区| a级毛片免费高清观看在线播放| 夜夜躁狠狠躁天天躁| 亚洲av.av天堂| 久久精品国产亚洲av天美| avwww免费| 亚洲第一欧美日韩一区二区三区| 少妇被粗大猛烈的视频| 亚洲精品乱码久久久v下载方式| 美女大奶头视频| 美女高潮的动态| 嫩草影视91久久| 国产精品美女特级片免费视频播放器| 国产爱豆传媒在线观看| 亚洲人成伊人成综合网2020| 久久人人精品亚洲av| xxxwww97欧美| 亚洲自拍偷在线| 女人被狂操c到高潮| 村上凉子中文字幕在线| 老司机午夜福利在线观看视频| 精品久久久久久,| 又爽又黄无遮挡网站| 国产三级在线视频| 亚洲一区高清亚洲精品| 亚洲国产精品999在线| 国产精品伦人一区二区| 精华霜和精华液先用哪个| 亚洲精品456在线播放app | 亚洲av第一区精品v没综合| 好看av亚洲va欧美ⅴa在| or卡值多少钱| 一本综合久久免费| 亚洲av第一区精品v没综合| 亚洲最大成人中文| 欧美在线一区亚洲| 在线观看免费视频日本深夜| 久久午夜亚洲精品久久| 在线观看av片永久免费下载| 国内精品久久久久久久电影| 国产精品久久久久久亚洲av鲁大| xxxwww97欧美| 精品久久久久久久久亚洲 | 国产在视频线在精品| 久久国产乱子伦精品免费另类| 一级av片app| 1000部很黄的大片| 亚洲精品影视一区二区三区av| 少妇高潮的动态图| 欧美xxxx黑人xx丫x性爽| .国产精品久久| 亚洲精品影视一区二区三区av| 一级黄色大片毛片| 午夜精品在线福利| 嫩草影视91久久| 99久久精品国产亚洲精品| 2021天堂中文幕一二区在线观| 亚洲成a人片在线一区二区| 欧美日本亚洲视频在线播放| 免费看日本二区| 一级a爱片免费观看的视频| 亚洲av电影不卡..在线观看| 久久精品国产清高在天天线| 中文字幕久久专区| 午夜福利免费观看在线| 好看av亚洲va欧美ⅴa在| 搡老妇女老女人老熟妇| 免费av不卡在线播放| 97热精品久久久久久| 国产精品亚洲一级av第二区| 午夜精品久久久久久毛片777| ponron亚洲| 亚洲人成伊人成综合网2020| 婷婷色综合大香蕉| 直男gayav资源| 亚洲,欧美精品.| 淫妇啪啪啪对白视频| 国产一区二区在线av高清观看| 午夜日韩欧美国产| 成人午夜高清在线视频| 熟女电影av网| 国产美女午夜福利| 亚洲精品日韩av片在线观看| 色尼玛亚洲综合影院| 亚洲av免费在线观看| 天堂网av新在线| 国产精品日韩av在线免费观看| 中文资源天堂在线| 国产亚洲av嫩草精品影院| 亚洲精品在线美女| 51国产日韩欧美| 午夜免费成人在线视频| 欧美另类亚洲清纯唯美| 88av欧美| 欧美高清成人免费视频www| 深爱激情五月婷婷| 亚洲色图av天堂| 亚洲五月婷婷丁香| 成人国产一区最新在线观看| 欧美中文日本在线观看视频| 日本 av在线| 国产精品永久免费网站| 亚洲18禁久久av| 成人精品一区二区免费| 日日干狠狠操夜夜爽| 又紧又爽又黄一区二区| 午夜精品久久久久久毛片777| 精品人妻1区二区| 国产精品伦人一区二区| 91麻豆精品激情在线观看国产| 午夜精品一区二区三区免费看| 国产三级黄色录像| 国产伦在线观看视频一区| 99久久无色码亚洲精品果冻| 三级国产精品欧美在线观看| 色视频www国产| 欧美在线黄色| 日本一本二区三区精品| 国产一区二区三区在线臀色熟女| 少妇熟女aⅴ在线视频| 成人无遮挡网站| 身体一侧抽搐| 欧美3d第一页| 我的女老师完整版在线观看| av专区在线播放| 免费观看人在逋| 亚洲精品一卡2卡三卡4卡5卡| 91久久精品国产一区二区成人| 国产亚洲精品av在线| 日韩av在线大香蕉| 中文亚洲av片在线观看爽| 嫁个100分男人电影在线观看| 1000部很黄的大片| 一级av片app| 亚洲激情在线av| av专区在线播放| 亚洲最大成人av| 99久久九九国产精品国产免费| 99riav亚洲国产免费| 日韩欧美 国产精品| 夜夜夜夜夜久久久久| 精品国产三级普通话版| 精品一区二区三区av网在线观看| АⅤ资源中文在线天堂| 国产成人a区在线观看| 国产私拍福利视频在线观看| 精品无人区乱码1区二区| 欧美极品一区二区三区四区| 亚洲精品一卡2卡三卡4卡5卡| 欧美日韩乱码在线| 在线观看舔阴道视频| 噜噜噜噜噜久久久久久91| 国产aⅴ精品一区二区三区波| 久久精品国产亚洲av香蕉五月| 校园春色视频在线观看| av国产免费在线观看| 夜夜夜夜夜久久久久| 免费看a级黄色片| 桃色一区二区三区在线观看| 精品久久久久久,| 又爽又黄无遮挡网站| 亚洲国产精品sss在线观看| 女人被狂操c到高潮| 成人三级黄色视频| 一级av片app| 97人妻精品一区二区三区麻豆| 国产精品伦人一区二区| 久久国产乱子免费精品| 国产精品亚洲av一区麻豆| 99国产极品粉嫩在线观看| 欧美xxxx性猛交bbbb| 老熟妇仑乱视频hdxx| 精品一区二区三区视频在线| 欧美日韩黄片免| xxxwww97欧美| 色吧在线观看| 亚洲一区二区三区不卡视频| 国产精华一区二区三区| 亚洲成人精品中文字幕电影| 精品久久久久久久久久久久久| 九九热线精品视视频播放| 欧美激情久久久久久爽电影| 国产欧美日韩一区二区精品| www日本黄色视频网| 99热这里只有是精品50| 在线观看美女被高潮喷水网站 | 国产91精品成人一区二区三区| 亚洲三级黄色毛片| 亚洲欧美日韩卡通动漫| 一进一出抽搐gif免费好疼| 极品教师在线视频| 男女床上黄色一级片免费看| 国产探花极品一区二区| 1024手机看黄色片| 色噜噜av男人的天堂激情| 国内精品久久久久精免费| 日本a在线网址| 1000部很黄的大片| 欧美日韩亚洲国产一区二区在线观看| 哪里可以看免费的av片| 国产精品乱码一区二三区的特点| 最近在线观看免费完整版| 精品久久久久久久久亚洲 | 夜夜看夜夜爽夜夜摸| 国产69精品久久久久777片| 久久精品国产清高在天天线| 国产午夜精品论理片| 青草久久国产| 亚洲精品日韩av片在线观看| 免费在线观看日本一区| 两个人视频免费观看高清| 成人av一区二区三区在线看| 成熟少妇高潮喷水视频| 国产精品日韩av在线免费观看| 亚洲精品一卡2卡三卡4卡5卡| 一级黄片播放器| 久久99热这里只有精品18| 国产亚洲精品综合一区在线观看| 俺也久久电影网| 精品人妻视频免费看| 中文字幕久久专区| 夜夜爽天天搞| 精品人妻1区二区| 男女那种视频在线观看| 精品免费久久久久久久清纯| 每晚都被弄得嗷嗷叫到高潮| 亚洲av电影不卡..在线观看| a级毛片免费高清观看在线播放| 麻豆成人av在线观看| 99热这里只有是精品在线观看 | 在线观看舔阴道视频| 久久性视频一级片| 亚洲精品一区av在线观看| 乱码一卡2卡4卡精品| 午夜亚洲福利在线播放| 日韩欧美 国产精品| 欧美国产日韩亚洲一区| 深爱激情五月婷婷| 两个人的视频大全免费| 亚洲国产精品久久男人天堂| 综合色av麻豆| 国产视频一区二区在线看| 亚洲精品456在线播放app | 亚洲人成伊人成综合网2020| 久久久久久久久久黄片| 欧美激情久久久久久爽电影| 少妇人妻精品综合一区二区 | 亚洲国产精品久久男人天堂| 黄色视频,在线免费观看| 亚洲一区二区三区色噜噜| 精品日产1卡2卡| 床上黄色一级片| 国产成+人综合+亚洲专区| 欧美黑人巨大hd| 日韩大尺度精品在线看网址| 白带黄色成豆腐渣| 国产精品久久视频播放| 日韩欧美 国产精品| x7x7x7水蜜桃| 精品人妻1区二区| 一卡2卡三卡四卡精品乱码亚洲| 在线免费观看的www视频| 免费在线观看日本一区| av在线天堂中文字幕| 99久久九九国产精品国产免费| 禁无遮挡网站| 欧美+亚洲+日韩+国产| 成熟少妇高潮喷水视频| 亚州av有码| 国产熟女xx| 一区二区三区免费毛片| 日韩免费av在线播放| 精品午夜福利视频在线观看一区| 高清日韩中文字幕在线|