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

    晉北地區(qū)1986~2010年土地利用/覆被變化的驅(qū)動(dòng)力

    2016-10-13 21:38:03徐小明杜自強(qiáng)申小雨
    中國(guó)環(huán)境科學(xué) 2016年7期
    關(guān)鍵詞:工礦驅(qū)動(dòng)力林地

    徐小明,杜自強(qiáng),張 紅,馮 凌,申小雨

    ?

    晉北地區(qū)1986~2010年土地利用/覆被變化的驅(qū)動(dòng)力

    徐小明1*,杜自強(qiáng)1,張 紅2,馮 凌2,申小雨2

    (1.山西大學(xué)黃土高原研究所,山西 太原 030006;2.山西大學(xué)環(huán)境與資源學(xué)院,山西 太原 030006)

    為定量研究不同時(shí)期區(qū)域土地利用變化的不同驅(qū)動(dòng)特征,以地處農(nóng)牧交錯(cuò)帶的晉北地區(qū)為例,采用非監(jiān)督分類與人工目視解譯相結(jié)合的方法,獲取研究區(qū)1986~2010年間的土地利用/覆被變化(LUCC)狀況,采用典范對(duì)應(yīng)分析(CCA)方法定量分析了不同階段LUCC的自然與人為驅(qū)動(dòng)特征.結(jié)果表明:1)CCA在分析長(zhǎng)時(shí)期序列LUCC驅(qū)動(dòng)力上有較好的解釋力,CCA排序可以很好地提取LUCC與自然、人為因素之間的關(guān)系,反映不同時(shí)期土地利用變化受到的驅(qū)動(dòng)作用;2)研究區(qū)土地利用類型以耕地、草地、林地為主,土地利用變化類型以耕地與草地、林地與草地之間的互相轉(zhuǎn)移為主;3)從不同時(shí)期LUCC的驅(qū)動(dòng)因子來看,人口密度、人均國(guó)內(nèi)生產(chǎn)總值(GDP)、降水、坡度和高程是研究區(qū)各時(shí)期LUCC的主要驅(qū)動(dòng)力;研究區(qū)在1986~1995年間LUCC的主要驅(qū)動(dòng)力來自于人口增長(zhǎng)和經(jīng)濟(jì)發(fā)展的需求,而隨著時(shí)間的推進(jìn),坡度和高程等地形因素的驅(qū)動(dòng)作用也愈來愈重要;4)從不同地類變化的驅(qū)動(dòng)因子來看,耕地、草地等地類向居民用地、工礦用地轉(zhuǎn)移的主要驅(qū)動(dòng)力是人口和經(jīng)濟(jì)的發(fā)展;而林地,草地等土地覆被類型之間的轉(zhuǎn)移主要受到坡度,高程和降水等自然因素的驅(qū)動(dòng).

    土地利用/覆被變化;驅(qū)動(dòng)力;典范對(duì)應(yīng)分析;不同時(shí)期

    目前,土地利用/覆被變化(LUCC)研究受到國(guó)內(nèi)外學(xué)者的廣泛關(guān)注[1-2].LUCC是一種過程,體現(xiàn)著人為和自然因素的雙重影響,空間差異性大,各種環(huán)境因子的作用程度也不盡相同,長(zhǎng)期以來一直缺乏比較系統(tǒng)、定量的方法來探討LUCC與驅(qū)動(dòng)力因子的響應(yīng)關(guān)系.

    目前研究認(rèn)為,氣候、地形、土壤、水文、人口、經(jīng)濟(jì)、技術(shù)、政治、乃至價(jià)值觀是LUCC的主要驅(qū)動(dòng)力[3-5],一些學(xué)者采用了諸如多元線性回歸[6]、多元邏輯回歸[7]、模擬退火算法[8]、典范對(duì)應(yīng)分析(CCA)[9-11]等方法對(duì)LUCC的驅(qū)動(dòng)因子進(jìn)行了定量分析.其中,CCA分析將多元回歸和對(duì)應(yīng)分析結(jié)合起來,是一種直接排序的方法,在生態(tài)學(xué)研究中得到了廣泛的應(yīng)用[12-14].這種方法在宏觀的LUCC驅(qū)動(dòng)力分析中已有少量應(yīng)用[9-11,15],該方法可以有效地解釋LUCC與驅(qū)動(dòng)因素之間復(fù)雜的非線性關(guān)系,是定量分析LUCC驅(qū)動(dòng)力的有力工具[11].但是,這些研究大多僅分析某一時(shí)期的LUCC驅(qū)動(dòng)力[9-11],忽略了LUCC驅(qū)動(dòng)力系統(tǒng)的縱向性與過程性[16],因?yàn)樵诓煌臍v史時(shí)期,LUCC的驅(qū)動(dòng)力也有所不同[17].

    本研究選擇山西省晉北(包括晉西北)地區(qū)為研究區(qū),該區(qū)是中國(guó)北方生態(tài)脆弱的農(nóng)牧交錯(cuò)帶的重要組成部分,也是環(huán)京津唐地區(qū)防風(fēng)固沙、綠化生態(tài)屏障建設(shè)的重點(diǎn)地區(qū)之一.該區(qū)土地利用結(jié)構(gòu)在近幾十年來發(fā)生了復(fù)雜而劇烈的變化,同時(shí),近年來實(shí)施的退耕還林等政策也對(duì)這一農(nóng)牧交錯(cuò)區(qū)的土地利用/覆被狀況產(chǎn)生了顯著的影響,因此,這一地區(qū)是開展不同時(shí)期LUCC驅(qū)動(dòng)力變化研究的理想?yún)^(qū)域.

    本文在對(duì)不同時(shí)期土地利用變化分析的基礎(chǔ)上,將影響LUCC的驅(qū)動(dòng)因子與LUCC進(jìn)行CCA分析,揭示土地利用變化在不同環(huán)境梯度上的響應(yīng),從而系統(tǒng)地揭示區(qū)域LUCC驅(qū)動(dòng)力在過去幾十年間的變化.

    1 研究區(qū)與方法

    1.1 研究區(qū)概況

    研究區(qū)包括晉西北高原和大同盆地,位于110°56′E~114°32′E,38°39′N~40°44′N,行政區(qū)劃上包括朔州市全境、忻州市和大同市的部分縣市,總面積約3.17萬km2.區(qū)內(nèi)地形復(fù)雜多變,東部為大同盆地,西部為地形起伏的黃土丘陵區(qū).黃土廣泛分布在該區(qū)內(nèi).

    研究區(qū)屬于大陸性半干旱季風(fēng)性氣候,四季分明,全年溫差較大,平均氣溫在4.6~6.8oC之間.年日照時(shí)數(shù)約2300~2900h,降水量為380~ 460mm,年均蒸發(fā)量達(dá)2000mm,無霜期100~ 130d.

    1.2 數(shù)據(jù)源

    1.2.1 土地利用/覆被解譯 下載研究區(qū)1986年、1995年、2000年和2010年的Landsat TM 4~5遙感影像.分別選取該年份7~10月份,云覆蓋率小于5%的影像進(jìn)行進(jìn)一步處理和分析.

    1.2.2 驅(qū)動(dòng)力因子 選取直接或間接影響LUCC的氣候、地形和人類活動(dòng)因子,包括人口密度、人均GDP、道路距離、城鎮(zhèn)距離、高程、坡度、坡向、1km內(nèi)地形起伏度、年平均溫度和年平均降水10個(gè)具體的因子.各數(shù)據(jù)來源見表1.

    道路距離和城鎮(zhèn)距離由ArcGIS計(jì)算研究區(qū)內(nèi)各點(diǎn)距主要公路和居民用地的歐氏距離獲得.地形起伏度是指1km內(nèi)最高點(diǎn)海拔高度與最低點(diǎn)海拔高度的差值[11].溫度和降水采用ArcGIS對(duì)研究區(qū)內(nèi)及周邊7個(gè)氣象站點(diǎn)(空間位置見圖1)年均溫度,降水進(jìn)行反距離權(quán)重插值獲得.

    表1 數(shù)據(jù)來源 Table 1 The data sources

    1.3 研究方法

    1.3.1 土地利用/覆被解譯 在Erdas Imagine 2013和ArcGIS 10.2.1軟件平臺(tái)的支持下,采用計(jì)算機(jī)非監(jiān)督分類和人工目視判讀相結(jié)合的方式進(jìn)行土地利用信息解譯.首先對(duì)獲取的影像進(jìn)行預(yù)處理,然后執(zhí)行非監(jiān)督分類,將影像分為30類.建立包括耕地、林地、草地、居民用地、工礦用地、水域、鹽堿地和裸地8種土地利用/覆被方式的分類系統(tǒng).對(duì)照TM 5-4-3波段的假彩合成圖像,同時(shí)參考研究區(qū)植被、土壤、地形等相關(guān)圖件和Google Earth影像進(jìn)行目視判讀解譯.這種方法既能保證分類的精度,又能節(jié)省解譯時(shí)間[16].采用混淆誤差矩陣的方法,即以Kappa系數(shù)來評(píng)價(jià)土地利用分類的位置精度.

    1.3.2 土地利用變化區(qū)域的選取 采用空間疊加的方法提取研究區(qū)內(nèi)1986~1995年、1995~ 2000年、2000~2010年和1986~2010年4個(gè)時(shí)段的土地利用變化斑塊.為了減少數(shù)據(jù)冗余,去除細(xì)小斑塊的影響,篩選變化面積在研究區(qū)總面積0.01%(約3.17km2)及以上的斑塊進(jìn)行土地利用轉(zhuǎn)移以及驅(qū)動(dòng)力分析.

    1.3.3 典范對(duì)應(yīng)分析 將4個(gè)時(shí)間段的土地利用變化斑塊處理為(0,1)二元數(shù)據(jù)矩陣[11],這一矩陣中,1代表某一類型土地利用變化的發(fā)生,0代表該類變化未發(fā)生.利用區(qū)域統(tǒng)計(jì)功能構(gòu)建驅(qū)動(dòng)力因子矩陣.

    將土地利用變化矩陣和驅(qū)動(dòng)力矩陣輸入CANOCO V4.56中,分別分析4個(gè)時(shí)段土地利用變化的驅(qū)動(dòng)力.應(yīng)用Monte Carlo檢驗(yàn)(499次置換)檢測(cè)驅(qū)動(dòng)力因子與土地利用變化之間是否存在統(tǒng)計(jì)意義上顯著的相關(guān)關(guān)系.本研究的技術(shù)路線見圖2.

    2 結(jié)果與討論

    2.1 研究區(qū)土地利用變化的主要特征

    遙感影像解譯得到的研究區(qū)土地利用空間分布圖見圖3.經(jīng)檢驗(yàn),研究區(qū)個(gè)年份的土地利用分類Kappa系數(shù)均大于0.80,分類結(jié)果良好[18].

    從圖3可知,研究區(qū)耕地在東部大同盆地分布較為集中,呈現(xiàn)東北—西南帶狀分布.在西側(cè)的黃土丘陵區(qū)則與林地、草地交錯(cuò)分布.林地在1986~2000年間主要集中在研究區(qū)中東部,在西部的河曲、保德、偏關(guān)等地分布較少,但至2010年時(shí),西部的林地面積也有所增加.草地的分布較為分散,在整個(gè)研究區(qū)與林地、耕地混雜分布.居民用地分布在各縣市的城鎮(zhèn)區(qū).工礦用地除了集中在平魯區(qū)南側(cè)以外,在朔城區(qū)和南郊區(qū)也有分散分布.鹽堿地主要集中分布在山陰縣和應(yīng)縣;裸地和水域的分布則比較分散.

    由圖4可知,耕地、草地和林地是研究區(qū)主要的土地利用類型,三者面積共占研究區(qū)總面積的94%~97%.耕地面積在1986~2000年間變化不大,但2000~2010年間有所減少.林地面積在1986~1995年間有所減少,但1995~2010年又呈持續(xù)增加的趨勢(shì).草地面積在1986~2000年間變化不大,但2000~2010年間有所減少.居民用地、工礦用地和裸地面積呈持續(xù)增加的趨勢(shì),而水域面積則持續(xù)減少,鹽堿地面積呈波動(dòng)增加的趨勢(shì).

    2.2 研究區(qū)不同歷史時(shí)期土地利用轉(zhuǎn)移分析

    由表2可見,1986~1995年間,主要的土地利用轉(zhuǎn)移類型為耕地和草地、草地和林地之間的互相轉(zhuǎn)移,占總轉(zhuǎn)移面積的91.40%;1995~2000年間的主要轉(zhuǎn)移類型與前一時(shí)期一致,4種類型占總變化面積的92.00%;2000~2010年間,除了以上4種轉(zhuǎn)移類型以外,耕地轉(zhuǎn)移為林地的比例增加到了10.25%,這5種土地利用轉(zhuǎn)移類型占總轉(zhuǎn)移面積的91.96%;1986~2010年和2000~2010年的轉(zhuǎn)移狀況類似,5種主要轉(zhuǎn)移類型的總概率為91.06%.不同時(shí)期的土地利用轉(zhuǎn)移主要發(fā)生在耕地、草地和林地之間的相互轉(zhuǎn)移上.其中,耕地和草地、林地和草地之間的互相轉(zhuǎn)移始終是研究區(qū)內(nèi)最主要的土地利用轉(zhuǎn)移類型.草地和耕地互相轉(zhuǎn)移的比例較高,充分顯示了研究區(qū)農(nóng)牧交錯(cuò)區(qū)的特性,而林地和草地的互相轉(zhuǎn)移是植被的自然演替和植樹造林交互作用的結(jié)果.

    表2 土地利用轉(zhuǎn)移概率矩陣(%) Table 2 Probability matrix of land use transition (%)

    注:行表示研究期初的土地利用類型,列表示研究期末的土地利用類型,空白表示該類型轉(zhuǎn)移未發(fā)生.

    耕地、草地和林地在2000年前面積變化不大.從2000年開始,耕地面積顯著減少,草地略有減少,而林地面積增加.這與1999年以后區(qū)內(nèi)推行退耕還林政策,部分坡耕地轉(zhuǎn)移為林地和草地有關(guān).土地利用轉(zhuǎn)移矩陣也證明了這一進(jìn)程,從2000年起,耕地轉(zhuǎn)移為林地的比例增大.由于該區(qū)地處三北防護(hù)林工程[19]和京津風(fēng)沙源治理工程區(qū)內(nèi)[20],一些荒山草地轉(zhuǎn)移為林地,促進(jìn)了區(qū)域內(nèi)林地面積的增加,2000~2010年間的土地利用轉(zhuǎn)移矩陣也表明2000年以后草地轉(zhuǎn)移為林地的比例遠(yuǎn)大于林地轉(zhuǎn)移為草地的比例.

    研究區(qū)居民用地的增加主要來自于耕地、工礦用地和草地,新增的工礦用地則主要來自草地、耕地和林地.水域、鹽堿地和裸地主要是與耕地和草地之間發(fā)生互相轉(zhuǎn)移.區(qū)內(nèi)人口的增加促進(jìn)了居民用地的持續(xù)增加,特別是大同市和朔州市建成區(qū)的面積貢獻(xiàn)巨大,大同市城市建成區(qū)面積由1985年的41km2[21]增加到了2010年的108km2,而朔城區(qū)城市建成區(qū)面積也由1990年的8.8km2[22]增加到2010年的40km2.居民用地的增加,2000年以前主要來自于耕地,而2000年以后,工礦用地和草地轉(zhuǎn)移為居民用地的比例增加,說明城鎮(zhèn)的發(fā)展,首先會(huì)占用耕地,進(jìn)而侵占天然植被;同時(shí),礦山生活區(qū)因人口聚集也逐漸轉(zhuǎn)移為居民用地.此外,研究區(qū)是我國(guó)重要的煤炭生產(chǎn)基地,大同煤礦、平朔煤礦均位于本區(qū)內(nèi),采礦業(yè)的發(fā)展促使了區(qū)內(nèi)工礦用地的持續(xù)增加,這些增加的工礦用地大多來自于草地和耕地.

    2.3 研究區(qū)不同時(shí)期LUCC驅(qū)動(dòng)力分析

    2.3.1 驅(qū)動(dòng)力因子選擇及相關(guān)性分析 LUCC體現(xiàn)著人為和自然因素的雙重影響,因此其驅(qū)動(dòng)力也來自于人為和自然的兩個(gè)方面[6].本研究對(duì)各驅(qū)動(dòng)力因子進(jìn)行相關(guān)分析表明地形起伏度與坡度的相關(guān)性較強(qiáng)(相關(guān)系數(shù)大于0.97),因此在分析中將地形起伏度因子作為冗余變量去除. 共保留9個(gè)相關(guān)性較小、具有典型性和獨(dú)立性的驅(qū)動(dòng)力指標(biāo).

    2.3.2 不同歷史時(shí)期LUCC的驅(qū)動(dòng)因子 1986~ 1995年間,第一軸與人均GDP的相關(guān)系數(shù)最大,第二軸與人口密度的相關(guān)系數(shù)最大.隨著第一軸數(shù)值的增加,人均GDP增加,人口密度減少,坡度變大;隨著第二軸數(shù)值的增加,人口密度和人均GDP均增加,坡度變小.這說明這一期間,研究區(qū)土地利用變化的最主要的驅(qū)動(dòng)力來自于人口增長(zhǎng)和經(jīng)濟(jì)發(fā)展.

    1995~2000年間,;而第二軸則與坡度的相關(guān)系數(shù)最大.隨著第一軸數(shù)值的增加,人均GDP和人口密度增加,坡度變小;隨著第二軸數(shù)值的增加,坡度,高程和降水增加.說明這一時(shí)期內(nèi),人口增長(zhǎng)和經(jīng)濟(jì)發(fā)展仍是主要驅(qū)動(dòng)因素,但地形因素的驅(qū)動(dòng)作用在逐漸增強(qiáng).

    2000~2010年間,第一軸與人均GDP的相關(guān)系數(shù)最大,第二軸則與坡度的相關(guān)系數(shù)最大.隨著第一軸數(shù)值的增加,人均GDP和人口密度增加;隨著第二軸數(shù)值的增加,坡度變小,高程降低,降水減少.這一時(shí)期的土地利用變化受到了人口、經(jīng)濟(jì)發(fā)展以及地形因素的雙重驅(qū)動(dòng).

    從1986~2010年總體來看,第一軸與坡度的相關(guān)系數(shù)最大,且為負(fù)相關(guān),第二軸與人均GDP相關(guān)系數(shù)最大,為正相關(guān).隨著第一軸數(shù)值的增加,坡度減小,高程降低,降水減少;隨著第二軸數(shù)值的增加,人均GDP和人口密度增加.說明1986~2010年間研究區(qū)主要的LUCC驅(qū)動(dòng)力包括人均GDP、坡度、高程、降水和人口密度.

    2.3.3 不同地類變化的驅(qū)動(dòng)因子分析 由圖5和表3可見,1986~1995年間,草地轉(zhuǎn)移為工礦用地主要發(fā)生在人均GDP較高的地區(qū);而耕地轉(zhuǎn)移為居民用地、林地轉(zhuǎn)移為耕地均發(fā)生在人口密度較大但人均GDP較低的地區(qū)(圖5a).說明這一時(shí)期內(nèi),人口和經(jīng)濟(jì)增長(zhǎng)促進(jìn)了耕地、草地等用地類型向工礦用地和居民用地的轉(zhuǎn)移.這一時(shí)期處于我國(guó)改革開放的初期,區(qū)內(nèi)GDP的增加促進(jìn)了其它地類向工礦用地的增加;而在經(jīng)濟(jì)尚不發(fā)達(dá),但人口壓力較大的地區(qū),易發(fā)生占田造城、毀林種地等開發(fā)活動(dòng).

    1995~2000年間,耕地向居民用地和林地的轉(zhuǎn)移主要發(fā)生在人口密度較大、人均GDP較高的地區(qū).林地和草地之間的互相轉(zhuǎn)移主要發(fā)生在坡度較大、高程較高、降水較多的區(qū)域(圖5b).人口的增加促使了居民用地的增加,而居民用地最便利的來源是城鎮(zhèn)周邊的耕地.同時(shí),由于1999年起推行的退耕還林政策的實(shí)施,一些坡耕地上植樹造林,圖5(b)第一象限中也提取到了這一變化類型,說明這一政策起初發(fā)生在經(jīng)濟(jì)較好、人口較多的地區(qū).此外,林地和草地的互相轉(zhuǎn)移,更易發(fā)生在人為影響較小(高程較高、坡度較大)且降水條件較好的地區(qū).

    2000~2010年間,耕地、林地、草地等地類向工礦用地和居民用地的轉(zhuǎn)移,以及工礦用地向居民用地的轉(zhuǎn)移均發(fā)生在人口密度和人均GDP較大的區(qū)域(圖5c).這一時(shí)期正處煤炭的“黃金十年”[23],工礦用地面積劇增,伴隨著GDP的快速增長(zhǎng);同時(shí),區(qū)內(nèi)人口的繼續(xù)增加也促進(jìn)了居民用地的增加,其來源除了城鎮(zhèn)周邊的耕地以外,原本屬于工礦用地的礦山生活區(qū)因人口聚集也逐漸轉(zhuǎn)移為居民用地.

    2.4 討論

    在本研究中,如僅對(duì)1986~2010年間的LUCC進(jìn)行分析,僅能推斷出人口密度和人均GDP對(duì)耕地,草地和工礦用地向居民用地的轉(zhuǎn)移起到正向驅(qū)動(dòng)作用,很難針對(duì)該結(jié)果推斷出這種驅(qū)動(dòng)發(fā)生的原因和機(jī)理;而通過不同時(shí)期的CCA分析可以判斷出不同時(shí)期的驅(qū)動(dòng)機(jī)理,如1995~2000年間的退耕還林政策的落實(shí)等,進(jìn)而對(duì)LUCC的驅(qū)動(dòng)過程形成更深入的了解,同時(shí)也可為政府部門進(jìn)行決策提供更加科學(xué)合理的依據(jù).

    本文的結(jié)果證明:在經(jīng)濟(jì)水平較低的時(shí)期,人類通過不斷開荒種地維持發(fā)展,這時(shí)LUCC的驅(qū)動(dòng)力很可能來自于生產(chǎn)的發(fā)展.而當(dāng)社會(huì)經(jīng)濟(jì)發(fā)展到一定程度,土地資源被不合理地開發(fā)直至面臨著一些嚴(yán)重的生態(tài)環(huán)境問題時(shí),政府部門可能會(huì)適時(shí)出臺(tái)一些政策,通過不同的手段來限制不合理的開發(fā),例如退耕還林,圍欄禁牧等措施,使得人為活動(dòng)的影響降低,此時(shí),自然因素的驅(qū)動(dòng)作用就可能會(huì)變得顯著.

    表3 圖5中土地利用轉(zhuǎn)移類型(%) Table 3 The land use transition types in Fig.5 (%)

    注:列表示某一時(shí)期初土地利用,行表示某一時(shí)期末土地利用.

    3 結(jié)論

    3.1 晉北地區(qū)土地利用類型以耕地、草地和林地為主,土地利用的轉(zhuǎn)移以耕地和草地、林地和草地之間的互相轉(zhuǎn)移為主,呈現(xiàn)典型的農(nóng)牧交錯(cuò)區(qū)特性.2000年以后,耕地轉(zhuǎn)移為林地的比例也逐漸增加,體現(xiàn)了退耕還林政策實(shí)施的效果.隨著人口和經(jīng)濟(jì)的增長(zhǎng),區(qū)內(nèi)居民用地和工礦用地的面積也持續(xù)增加.

    3.2 CCA排序可以很好地提取LUCC與自然、人為因素之間的關(guān)系,反應(yīng)不同時(shí)期土地利用變化受到的驅(qū)動(dòng)作用.從不同歷史時(shí)期LUCC的驅(qū)動(dòng)因子來看,研究區(qū)早期LUCC的主要驅(qū)動(dòng)力來自于人口增長(zhǎng)和經(jīng)濟(jì)發(fā)展的需求,與城市擴(kuò)張、糧食需求增加有著密切的聯(lián)系;2000年前后,研究區(qū)內(nèi)退耕還林、防護(hù)林建設(shè)等生態(tài)保護(hù)工程先后開始實(shí)施,這些工程受地形、氣候因素的因素的影響較大,因此,這些自然因素的驅(qū)動(dòng)作用也越來越重要.從不同地類的驅(qū)動(dòng)因子來看,研究區(qū)人口和經(jīng)濟(jì)的發(fā)展是耕地、草地等地類向居民用地、工礦用地轉(zhuǎn)移的主要驅(qū)動(dòng)力;而林地、草地等土地覆被類型之間的轉(zhuǎn)移主要受到坡度,高程和降水等自然因素的驅(qū)動(dòng).

    [1] Luyssaert S, Jammet M, Stoy P, et al. Land management and land-cover change have impacts of similar magnitude on surface temperature [J]. Nature Climate Change, 2014,4(5):389-393.

    [2] 徐小明.吉林西部水田土壤碳庫時(shí)空模擬及水稻生產(chǎn)的碳足跡研究 [D]. 長(zhǎng)春:吉林大學(xué), 2011.

    [3] 趙 陽,張 藝,涂志華,等.基于生態(tài)服務(wù)價(jià)值的多目標(biāo)水源地土地利用結(jié)構(gòu)優(yōu)化 [J]. 中國(guó)環(huán)境科學(xué), 2014,34(1):232-238.

    [4] Kanianska R, Kizeková M, Nová?ek J, et al. Land-use and land-cover changes in rural areas during different political systems: A case study of Slovakia from 1782 to 2006 [J]. Land Use Policy, 2014,36:554-566.

    [5] 石曉麗,史文嬌.氣候變化和人類活動(dòng)對(duì)耕地格局變化的貢獻(xiàn)歸因綜述 [J]. 地理學(xué)報(bào), 2015,70(9):1463-1476.

    [6] 王晨野,湯 潔,李昭陽,等.吉林西部土地利用/覆被時(shí)空變化驅(qū)動(dòng)力分析 [J]. 生態(tài)環(huán)境, 2008,17(5):1914-1920.

    [7] Serra P, Pons X, Saurí D. Land-cover and land-use change in a Mediterranean landscape: a spatial analysis of driving forces integrating biophysical and human factors [J]. Applied Geography, 2008,28(3):189-209.

    [8] 位 欣,劉耀林,姚 鵬.基于模擬退火遺傳算法的土地利用變化驅(qū)動(dòng)力研究 [J]. 中國(guó)土地科學(xué), 2008,22(7):34-37.

    [9] 李 鵬,于書霞.基于典型對(duì)應(yīng)分析的滇池流域土地利用驅(qū)動(dòng)力分析 [J]. 中國(guó)人口.資源與環(huán)境, 2014,24(3):123-126.

    [10] 邢容容,馬安青,劉修錦,等.青島市土地利用變化類型與影響因子關(guān)系的CCA研究 [J]. 安徽農(nóng)業(yè)科學(xué), 2015,43(15):298- 301+304.

    [11] 徐廣才,康慕誼,李亞飛.錫林郭勒盟土地利用變化及驅(qū)動(dòng)力分析 [J]. 資源科學(xué), 2011,33(4):690-697.

    [12] Li J, Zheng Y, Yan J, et al. Succession of plant and soil microbial communities with restoration of abandoned land in the Loess Plateau, China [J]. Journal of Soils and Sediments, 2013,13(4): 760-769.

    [13] Bu C, Zhang P, Wang C, et al. Spatial distribution of biological soil crusts on the slope of the Chinese Loess Plateau based on canonical correspondence analysis [J]. Catena, 2016,137:373- 381.

    [14] 劉奕伶,葛繼穩(wěn),李艷元,等.古夫河著生藻類優(yōu)勢(shì)種體積與水質(zhì)因子的相關(guān)性研究 [J]. 中國(guó)環(huán)境科學(xué), 2015,35(7):2182-2191.

    [15] Hietel E, Waldhardt R, Otte A. Analysing land-cover changes in relation to environmental variables in Hesse, Germany [J]. Landscape ecology, 2004,19(5):473-489.

    [16] Elmhagen B, Eriksson O, Lindborg R. Implications of climate and land-use change for landscape processes, biodiversity, ecosystem services, and governance [J]. Ambio, 2015,44(1):1-5.

    [17] 湯 潔,汪雪格,李昭陽,等.基于CA-Markov模型的吉林省西部土地利用景觀格局變化趨勢(shì)預(yù)測(cè) [J]. 吉林大學(xué)學(xué)報(bào):地球科學(xué)版, 2010,40(2):405-411.

    [18] Jiang W, Chen Z, Lei X, et al. Simulating urban land use change by incorporating an autologistic regression model into a CLUE-S model [J]. Journal of Geographical Sciences, 2015,25(7):836- 850.

    [19] 國(guó)家林業(yè)局.三北防護(hù)林工程概況[EB/OL]. http://www.forestry. gov.cn.2015-6-23.

    [20] 京津風(fēng)沙源治理工程二期規(guī)劃思路研究項(xiàng)目組.京津風(fēng)沙源治理工程二期規(guī)劃思路研究 [M]. 北京:中國(guó)林業(yè)出版社, 2013.

    [21] 大同市地方志編纂委員會(huì).大同市志 [M]. 北京:方志出版社, 2000.

    [22] 朔縣志編纂委員會(huì).朔縣志 [M]. 太原:山西古籍出版社, 1999.

    [23] 中國(guó)經(jīng)濟(jì)網(wǎng).解密中國(guó)煤炭“黃金十年”[EB/OL]. http://www. ce.cn/. 2015-6-19.

    * 責(zé)任作者, 副教授, xuxiaoming@sxu.edu.cn

    during 1986~2010

    XU Xiao-ming1*, DU Zi-qiang1, ZHANG Hong2, FENG Ling2, SHEN Xiao-yu2

    (1.Institute of Loess Plateau, Shanxi University, Taiyuan 030006, China;2.College of Environment and Resource Sciences, Shanxi University, Taiyuan 030006, China)., 2016,36(7):2154~2161

    The land use pattern in a certain region in different periods was usually driven by natural, social, and economical factors. The driving mechanisms varied in different periods. North Shanxi province, located in the farming-pastoral ecotone of northern China, was chosen as the research area. The objectives of this study were to analyze the land use/cover change (LUCC) pattern in north Shanxi province, and to estimate the different driving mechanisms of the LUCC in different historical periods. A comprehensive method combining the unsupervised classification and the visual interpretation was used to extract the LUCC information in the study area from 1986 to 2010. The canonical correspondence analysis was employed to assess the natural and anthropogenic driving forces of the LUCC. The results indicated that cropland, grassland, and woodland were major land use types in the whole investigation period. Main land-use transferring types included the interconversions between grassland and cropland, as well as woodland and grassland. Population density, per capita gross domestic product (GDP), precipitation, slope, and elevation were the major driving forces of the LUCC. In 1986~1995, the driving forces of the LUCC were mainly from the demands of the population growth and the economy development. Later, the driving forces from the topographical elements such as elevation and slope became more and more important. The economic development and population growth promote the transitions from grassland and cropland to residential land and industrial land. Meanwhile, the conversions between grassland and woodland are mainly driven by natural factors such as slope, elevation, and precipitation.

    land use/cover change;driving forces;canonical correspondence analysis;different periods

    X171.1,F301.24

    A

    1000-6923(2016)07-2154-08

    徐小明(1984-),男,山西五寨人,副教授,博士,主要研究方向?yàn)橥恋乩?覆被變化及碳循環(huán).發(fā)表論文10余篇.

    2015-11-12

    山西省科技重大專項(xiàng)(20121101011);山西省青年科技研究基金(2015021172);國(guó)家自然科學(xué)基金 (41401053,41161066);教育部留學(xué)回國(guó)人員科研啟動(dòng)基金

    猜你喜歡
    工礦驅(qū)動(dòng)力林地
    蘇州市開展工礦行業(yè)安全生產(chǎn)大檢查
    油價(jià)上漲的供需驅(qū)動(dòng)力能否持續(xù)
    溫暖厚實(shí),驅(qū)動(dòng)力強(qiáng)勁 秦朝 QM2018/QC2350前后級(jí)功放
    突出文化產(chǎn)業(yè)核心驅(qū)動(dòng)力
    以創(chuàng)新為驅(qū)動(dòng)力,兼具學(xué)院派的嚴(yán)謹(jǐn)態(tài)度 Q Acoustics
    丹東市林地分類研究
    淺談林地保護(hù)及恢復(fù)措施
    打開工礦燈市場(chǎng)前,你必須知道的散熱知識(shí)!
    懇請(qǐng)將黑龍江“工礦型”革命老區(qū)納入國(guó)家精準(zhǔn)扶貧規(guī)劃
    徐工礦用卡車在泰國(guó)登陸
    專用汽車(2014年10期)2014-12-17 01:37:17
    免费观看av网站的网址| 欧美精品国产亚洲| 夫妻午夜视频| 精华霜和精华液先用哪个| 亚洲伊人久久精品综合| 午夜免费鲁丝| 美女高潮的动态| 久热这里只有精品99| 免费观看在线日韩| 亚洲欧洲国产日韩| 久久精品久久精品一区二区三区| 99热这里只有是精品在线观看| 亚洲三级黄色毛片| 成人亚洲欧美一区二区av| 亚洲天堂av无毛| 九色成人免费人妻av| 欧美xxxx黑人xx丫x性爽| 美女视频免费永久观看网站| 熟女av电影| av在线蜜桃| 草草在线视频免费看| 欧美日韩视频高清一区二区三区二| 丝袜脚勾引网站| 中文字幕人妻熟人妻熟丝袜美| 日韩在线高清观看一区二区三区| 国产精品国产三级专区第一集| 菩萨蛮人人尽说江南好唐韦庄| 国产视频首页在线观看| 最近手机中文字幕大全| 99久久人妻综合| 国产在线视频一区二区| 26uuu在线亚洲综合色| 美女xxoo啪啪120秒动态图| 在线免费十八禁| 国产精品麻豆人妻色哟哟久久| 五月开心婷婷网| 91aial.com中文字幕在线观看| 亚洲精品国产成人久久av| 精品一区二区三区视频在线| 麻豆乱淫一区二区| 久热久热在线精品观看| 国产一级毛片在线| 美女视频免费永久观看网站| 最近最新中文字幕大全电影3| 夜夜爽夜夜爽视频| 色5月婷婷丁香| 国产在线免费精品| 午夜福利在线观看免费完整高清在| 久久久久久九九精品二区国产| 一区二区三区精品91| 欧美日韩视频精品一区| 在线观看一区二区三区| 亚洲自偷自拍三级| 丰满迷人的少妇在线观看| 国模一区二区三区四区视频| 亚洲欧美一区二区三区国产| 丝瓜视频免费看黄片| 亚洲人成网站在线观看播放| 中国美白少妇内射xxxbb| 岛国毛片在线播放| 免费av不卡在线播放| 少妇人妻精品综合一区二区| 高清在线视频一区二区三区| 亚洲欧美日韩东京热| 免费观看在线日韩| 三级国产精品欧美在线观看| 哪个播放器可以免费观看大片| 中国三级夫妇交换| 一边亲一边摸免费视频| 久久ye,这里只有精品| 亚洲aⅴ乱码一区二区在线播放| 美女福利国产在线 | 成年人午夜在线观看视频| 嘟嘟电影网在线观看| 精品一区二区三卡| 女性生殖器流出的白浆| 久久久午夜欧美精品| 久久99热这里只有精品18| 日韩一本色道免费dvd| 国产成人91sexporn| 久久久久人妻精品一区果冻| 国产成人a∨麻豆精品| 女性被躁到高潮视频| 中文字幕亚洲精品专区| 蜜桃在线观看..| 精品国产露脸久久av麻豆| 国产精品福利在线免费观看| 亚洲av男天堂| 色视频在线一区二区三区| 国产老妇伦熟女老妇高清| 女的被弄到高潮叫床怎么办| 日韩电影二区| 看十八女毛片水多多多| 亚洲精品乱码久久久久久按摩| 国产成人91sexporn| 嫩草影院新地址| 国产人妻一区二区三区在| 国产大屁股一区二区在线视频| 亚洲欧美日韩东京热| 国产有黄有色有爽视频| 国产男女内射视频| 免费大片18禁| 青春草国产在线视频| 一个人看视频在线观看www免费| 在线观看一区二区三区| 国内精品宾馆在线| 另类亚洲欧美激情| 91aial.com中文字幕在线观看| 免费观看无遮挡的男女| 在线观看一区二区三区激情| av在线观看视频网站免费| 激情 狠狠 欧美| 性色av一级| 免费大片黄手机在线观看| 国产伦精品一区二区三区视频9| 午夜老司机福利剧场| 国产日韩欧美在线精品| 这个男人来自地球电影免费观看 | 中文字幕人妻熟人妻熟丝袜美| 亚洲欧美精品自产自拍| 欧美性感艳星| 久久精品国产鲁丝片午夜精品| 夜夜看夜夜爽夜夜摸| 久久久国产一区二区| 99热6这里只有精品| 人人妻人人爽人人添夜夜欢视频 | 美女国产视频在线观看| 午夜激情福利司机影院| 伊人久久精品亚洲午夜| 2018国产大陆天天弄谢| 美女视频免费永久观看网站| 国产免费福利视频在线观看| 日韩强制内射视频| 蜜桃久久精品国产亚洲av| 久久女婷五月综合色啪小说| 免费高清在线观看视频在线观看| 久久久色成人| 国产视频内射| 亚洲无线观看免费| 成人毛片a级毛片在线播放| 国产一区亚洲一区在线观看| 国产精品嫩草影院av在线观看| 女的被弄到高潮叫床怎么办| 久久亚洲国产成人精品v| 全区人妻精品视频| 99久久精品热视频| 欧美日韩国产mv在线观看视频 | 国产精品一区二区在线观看99| 汤姆久久久久久久影院中文字幕| 国产亚洲欧美精品永久| 人人妻人人添人人爽欧美一区卜 | 色视频www国产| 日韩欧美一区视频在线观看 | 十八禁网站网址无遮挡 | 噜噜噜噜噜久久久久久91| 一级毛片久久久久久久久女| 最近的中文字幕免费完整| 国产大屁股一区二区在线视频| 三级经典国产精品| 午夜日本视频在线| 亚洲欧洲日产国产| 男人和女人高潮做爰伦理| 亚洲精品乱码久久久v下载方式| 内射极品少妇av片p| 亚洲av欧美aⅴ国产| 又粗又硬又长又爽又黄的视频| 插阴视频在线观看视频| 久久精品国产亚洲av天美| 国产在线免费精品| 婷婷色综合大香蕉| 最近中文字幕高清免费大全6| 亚洲精品日本国产第一区| 人人妻人人看人人澡| 少妇丰满av| 国产片特级美女逼逼视频| 午夜免费观看性视频| 亚洲欧美日韩东京热| 交换朋友夫妻互换小说| 欧美日韩视频精品一区| 有码 亚洲区| 国产老妇伦熟女老妇高清| 夜夜爽夜夜爽视频| 精品熟女少妇av免费看| av免费在线看不卡| 大片免费播放器 马上看| 熟妇人妻不卡中文字幕| 少妇高潮的动态图| 亚洲成人中文字幕在线播放| 色视频在线一区二区三区| 久久精品久久久久久久性| 国产毛片在线视频| 中文欧美无线码| 免费大片18禁| 欧美日韩亚洲高清精品| 内射极品少妇av片p| 一本—道久久a久久精品蜜桃钙片| 熟女人妻精品中文字幕| 大陆偷拍与自拍| a级毛片免费高清观看在线播放| 99热这里只有精品一区| 日日啪夜夜撸| 日本一二三区视频观看| 日韩在线高清观看一区二区三区| 国产亚洲91精品色在线| 国产久久久一区二区三区| 精品人妻偷拍中文字幕| 国产精品久久久久成人av| 日日摸夜夜添夜夜添av毛片| 国产精品一及| 亚洲国产成人一精品久久久| 亚洲真实伦在线观看| 亚洲综合精品二区| 大香蕉久久网| 国产精品女同一区二区软件| 日韩电影二区| 99热国产这里只有精品6| 下体分泌物呈黄色| av天堂中文字幕网| 久久久久久久精品精品| 午夜视频国产福利| 国产精品秋霞免费鲁丝片| 日韩 亚洲 欧美在线| 国产白丝娇喘喷水9色精品| 色哟哟·www| 九九爱精品视频在线观看| 日韩av在线免费看完整版不卡| 亚洲成人手机| 赤兔流量卡办理| 美女国产视频在线观看| 干丝袜人妻中文字幕| 免费观看av网站的网址| 啦啦啦中文免费视频观看日本| 中文字幕久久专区| 日本vs欧美在线观看视频 | 日韩一区二区视频免费看| 看非洲黑人一级黄片| 免费观看性生交大片5| 国产国拍精品亚洲av在线观看| 精品久久久久久久久亚洲| 中文精品一卡2卡3卡4更新| 国产午夜精品久久久久久一区二区三区| 久久韩国三级中文字幕| 国产成人a区在线观看| 中文天堂在线官网| 国产成人午夜福利电影在线观看| 亚洲精品日韩av片在线观看| 91精品伊人久久大香线蕉| 成年女人在线观看亚洲视频| 观看av在线不卡| 亚洲四区av| 精品视频人人做人人爽| 国产午夜精品久久久久久一区二区三区| 1000部很黄的大片| 亚洲精品色激情综合| 最近的中文字幕免费完整| 成人亚洲精品一区在线观看 | 男女国产视频网站| 建设人人有责人人尽责人人享有的 | 男人添女人高潮全过程视频| 亚洲精品乱码久久久久久按摩| 在线观看美女被高潮喷水网站| 精品99又大又爽又粗少妇毛片| 黄色日韩在线| 男女下面进入的视频免费午夜| 亚洲婷婷狠狠爱综合网| 亚洲欧美日韩另类电影网站 | 日韩一本色道免费dvd| 成人18禁高潮啪啪吃奶动态图 | 十分钟在线观看高清视频www | 欧美日韩视频高清一区二区三区二| 亚洲av在线观看美女高潮| 久久久亚洲精品成人影院| 青青草视频在线视频观看| 国产深夜福利视频在线观看| 国产精品一区二区性色av| 国产精品精品国产色婷婷| 中文字幕免费在线视频6| 精品国产一区二区三区久久久樱花 | 舔av片在线| 51国产日韩欧美| 日本黄色日本黄色录像| 国产精品国产三级国产av玫瑰| 18禁裸乳无遮挡免费网站照片| av视频免费观看在线观看| 国内揄拍国产精品人妻在线| 亚洲精品国产av成人精品| 国产av一区二区精品久久 | 91精品伊人久久大香线蕉| 99热这里只有是精品50| 亚洲av欧美aⅴ国产| a级毛片免费高清观看在线播放| 妹子高潮喷水视频| 成人亚洲精品一区在线观看 | 日韩 亚洲 欧美在线| 日本黄色日本黄色录像| 国产精品人妻久久久影院| 欧美xxxx性猛交bbbb| 亚洲久久久国产精品| 最新中文字幕久久久久| 毛片女人毛片| 国产又色又爽无遮挡免| 又大又黄又爽视频免费| 亚洲精品日本国产第一区| 99国产精品免费福利视频| 日韩精品有码人妻一区| 欧美日韩国产mv在线观看视频 | av一本久久久久| 九九爱精品视频在线观看| 国产乱人视频| 国产精品久久久久久av不卡| 久久精品熟女亚洲av麻豆精品| av在线观看视频网站免费| 少妇的逼好多水| 亚洲av成人精品一二三区| 一区二区三区乱码不卡18| 成人美女网站在线观看视频| 男女边摸边吃奶| 少妇猛男粗大的猛烈进出视频| 亚洲三级黄色毛片| 亚洲熟女精品中文字幕| 精品久久久精品久久久| 人妻 亚洲 视频| 欧美bdsm另类| 在现免费观看毛片| 亚洲国产精品999| 99久久精品热视频| 国产黄色视频一区二区在线观看| 色综合色国产| 国产亚洲午夜精品一区二区久久| 18禁动态无遮挡网站| 两个人的视频大全免费| 美女xxoo啪啪120秒动态图| av黄色大香蕉| av播播在线观看一区| 最后的刺客免费高清国语| 人妻一区二区av| 国产精品一及| 涩涩av久久男人的天堂| 国产成人免费观看mmmm| 又黄又爽又刺激的免费视频.| 久久精品国产自在天天线| 啦啦啦啦在线视频资源| 国产黄色免费在线视频| 人妻系列 视频| freevideosex欧美| 日本黄色片子视频| av免费观看日本| 国产一级毛片在线| 又大又黄又爽视频免费| 亚洲国产精品国产精品| 亚洲一级一片aⅴ在线观看| 精品人妻一区二区三区麻豆| 各种免费的搞黄视频| 国产亚洲av片在线观看秒播厂| 日本wwww免费看| 久久久久久久国产电影| 精品久久久久久电影网| 久久久久久久久久久丰满| 亚洲怡红院男人天堂| 久久精品国产亚洲av天美| 国产伦精品一区二区三区视频9| 免费黄网站久久成人精品| 在线观看一区二区三区| av黄色大香蕉| 亚洲av福利一区| 国产色爽女视频免费观看| 久久99热这里只有精品18| 亚洲av在线观看美女高潮| 午夜日本视频在线| 亚洲av免费高清在线观看| 欧美三级亚洲精品| 亚洲精品国产av成人精品| 久久av网站| 大片电影免费在线观看免费| av专区在线播放| 欧美少妇被猛烈插入视频| 午夜免费观看性视频| 一级毛片电影观看| 亚洲精华国产精华液的使用体验| 国产综合精华液| 日韩国内少妇激情av| 久久97久久精品| 看十八女毛片水多多多| 午夜福利影视在线免费观看| 在线天堂最新版资源| 午夜福利影视在线免费观看| 亚洲欧美精品自产自拍| 在线免费十八禁| 久久久亚洲精品成人影院| 亚洲精品第二区| 欧美亚洲 丝袜 人妻 在线| 国产乱来视频区| 亚洲av日韩在线播放| 丝袜脚勾引网站| 美女中出高潮动态图| 久久久国产一区二区| 99久久精品一区二区三区| 伦理电影免费视频| 亚洲久久久国产精品| 欧美zozozo另类| 国产精品99久久久久久久久| 国产亚洲欧美精品永久| 又爽又黄a免费视频| 交换朋友夫妻互换小说| 欧美亚洲 丝袜 人妻 在线| 王馨瑶露胸无遮挡在线观看| 黄色配什么色好看| 国产精品精品国产色婷婷| 久久久久久久亚洲中文字幕| 黄片wwwwww| 热99国产精品久久久久久7| 久久久久久久久久久免费av| 在线观看美女被高潮喷水网站| 精品亚洲成a人片在线观看 | 97在线视频观看| 国产久久久一区二区三区| 精品酒店卫生间| 欧美精品国产亚洲| 久久精品国产亚洲av天美| 自拍偷自拍亚洲精品老妇| 边亲边吃奶的免费视频| 精品久久久久久久久av| 亚洲人成网站在线播| 国产男女超爽视频在线观看| 亚洲av中文av极速乱| 亚洲国产精品成人久久小说| 免费少妇av软件| 欧美亚洲 丝袜 人妻 在线| 日韩免费高清中文字幕av| 日韩在线高清观看一区二区三区| 国内揄拍国产精品人妻在线| 女性被躁到高潮视频| 亚洲美女搞黄在线观看| 日韩欧美 国产精品| 国产成人免费无遮挡视频| 中国美白少妇内射xxxbb| 欧美精品一区二区免费开放| 能在线免费看毛片的网站| 麻豆乱淫一区二区| 最近中文字幕2019免费版| 日本猛色少妇xxxxx猛交久久| 一本久久精品| 黄色日韩在线| 美女福利国产在线 | www.色视频.com| 欧美丝袜亚洲另类| 久久午夜福利片| 中文资源天堂在线| 久久久久视频综合| 亚洲人成网站在线播| 亚洲欧美精品专区久久| 久久久色成人| 成年美女黄网站色视频大全免费 | 男人舔奶头视频| 色婷婷久久久亚洲欧美| 欧美日韩一区二区视频在线观看视频在线| 久久97久久精品| 日韩成人伦理影院| 欧美zozozo另类| av.在线天堂| 亚洲人与动物交配视频| 久久精品国产自在天天线| 久久久a久久爽久久v久久| 多毛熟女@视频| 99re6热这里在线精品视频| av在线播放精品| 九草在线视频观看| 日韩强制内射视频| 99热国产这里只有精品6| 99久久精品国产国产毛片| 欧美激情极品国产一区二区三区 | 美女福利国产在线 | 看非洲黑人一级黄片| 男人狂女人下面高潮的视频| av播播在线观看一区| 日韩一区二区三区影片| 国产免费一级a男人的天堂| 在现免费观看毛片| 亚洲欧美精品自产自拍| 免费观看的影片在线观看| 国产精品久久久久久久电影| 一级毛片我不卡| 99视频精品全部免费 在线| 国产成人a∨麻豆精品| 国产精品一区二区三区四区免费观看| 久久99热这里只频精品6学生| 久久精品国产亚洲av涩爱| 午夜免费鲁丝| 伊人久久精品亚洲午夜| 校园人妻丝袜中文字幕| 日韩,欧美,国产一区二区三区| 最近中文字幕高清免费大全6| 欧美激情国产日韩精品一区| 日韩成人av中文字幕在线观看| 久久鲁丝午夜福利片| 我要看黄色一级片免费的| 国产极品天堂在线| 国产爽快片一区二区三区| 久久久国产一区二区| 亚洲在久久综合| 日本一二三区视频观看| 久久99蜜桃精品久久| 精品亚洲成a人片在线观看 | 黑丝袜美女国产一区| 狂野欧美激情性bbbbbb| 国产一区二区在线观看日韩| 男女无遮挡免费网站观看| 免费大片18禁| 亚洲第一区二区三区不卡| 色视频www国产| 亚洲av男天堂| 久久亚洲国产成人精品v| 丰满少妇做爰视频| 只有这里有精品99| 日日啪夜夜爽| 在线观看国产h片| 黄色日韩在线| 久久人人爽人人片av| 日韩制服骚丝袜av| kizo精华| 国产在视频线精品| 国产午夜精品久久久久久一区二区三区| videos熟女内射| 日韩强制内射视频| 精品国产露脸久久av麻豆| 青春草国产在线视频| av免费在线看不卡| 伦精品一区二区三区| 免费观看的影片在线观看| 亚洲美女黄色视频免费看| 观看免费一级毛片| 国产综合精华液| 欧美xxⅹ黑人| 特大巨黑吊av在线直播| 有码 亚洲区| h日本视频在线播放| 狂野欧美白嫩少妇大欣赏| 成人国产麻豆网| 另类亚洲欧美激情| 精品久久久精品久久久| 中文欧美无线码| 国产成人精品一,二区| 久久久久久久精品精品| 永久免费av网站大全| 如何舔出高潮| 午夜激情久久久久久久| 五月开心婷婷网| 亚洲人与动物交配视频| 九九久久精品国产亚洲av麻豆| 18+在线观看网站| 午夜老司机福利剧场| 久久99蜜桃精品久久| 在线观看一区二区三区激情| 久久97久久精品| 国产乱人视频| 国产日韩欧美亚洲二区| 观看美女的网站| 两个人的视频大全免费| 高清毛片免费看| 我要看日韩黄色一级片| 成年美女黄网站色视频大全免费 | 建设人人有责人人尽责人人享有的 | 精品久久久久久久末码| 2022亚洲国产成人精品| 在线观看人妻少妇| 久久青草综合色| 国产永久视频网站| 久久久久久人妻| 国产高清国产精品国产三级 | 精品人妻偷拍中文字幕| 久久精品国产鲁丝片午夜精品| 成人18禁高潮啪啪吃奶动态图 | 久久这里有精品视频免费| 91精品国产国语对白视频| 亚洲欧美精品自产自拍| 亚洲人成网站在线播| 搡老乐熟女国产| 亚洲第一区二区三区不卡| 亚洲av不卡在线观看| 男女国产视频网站| 欧美日韩一区二区视频在线观看视频在线| 国产一级毛片在线| 男女边吃奶边做爰视频| 九九在线视频观看精品| 91精品国产九色| 国国产精品蜜臀av免费| 插逼视频在线观看| 特大巨黑吊av在线直播| 中文精品一卡2卡3卡4更新| 亚洲色图av天堂| 建设人人有责人人尽责人人享有的 | 毛片女人毛片| 久久精品久久久久久噜噜老黄| 欧美高清成人免费视频www| 高清欧美精品videossex| 日本色播在线视频| 国产av一区二区精品久久 | 成人亚洲精品一区在线观看 | 亚洲一区二区三区欧美精品| tube8黄色片| av在线播放精品| 国产黄色免费在线视频| 久久国产精品男人的天堂亚洲 | 亚洲av国产av综合av卡| 国产精品麻豆人妻色哟哟久久| 久久国产亚洲av麻豆专区| 国产亚洲5aaaaa淫片| 欧美97在线视频| av卡一久久| 深夜a级毛片| 赤兔流量卡办理| 热99国产精品久久久久久7| 亚洲国产av新网站| 不卡视频在线观看欧美| 在线免费观看不下载黄p国产|