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

    基于馬爾科夫模型的吐魯番市土地覆被動態(tài)變化研究

    2014-09-21 08:53:48何寶忠高敏華趙軍安
    水土保持研究 2014年5期
    關鍵詞:吐魯番市園地步長

    何寶忠, 高敏華, 趙軍安

    (1.新疆大學 資源與環(huán)境科學學院, 烏魯木齊 830046; 2.新疆大學 綠洲生態(tài)教育部重點實驗室, 烏魯木齊 830046)

    基于馬爾科夫模型的吐魯番市土地覆被動態(tài)變化研究

    何寶忠1,2, 高敏華1,2, 趙軍安1

    (1.新疆大學 資源與環(huán)境科學學院, 烏魯木齊 830046; 2.新疆大學 綠洲生態(tài)教育部重點實驗室, 烏魯木齊 830046)

    以吐魯番市1990年TM影像解譯數據和2005年、2011年的土地利用圖為數據源,通過建立馬爾科夫模型來預測該地區(qū)土地覆被變化情況。馬爾科夫模型的預測精度與模擬步長密切相關,因此,首先以1990年為基礎年,基于不同的模擬步長,預測2005年的各地類面積,并與其實際面積相比較,找出最佳模擬步長。其次,用該最佳步長預測2011年以后的土地覆被變化,由此分析該地區(qū)土地利用程度和土地利用空間變化動態(tài)度。研究表明:(1) 未利用地和草地相互轉換數量較大;(2) 2005年后耕地減少較為迅速,供需矛盾日益突出,基本沒有實現耕地占補平衡;(3) 在21年間水體減少近6 000 hm2,很可能與全球氣候變暖有關;(4) 由于建設用地和園地的變化而使得土地利用程度逐年提高,說明經濟的發(fā)展是土地覆被變化的重要因素;(5) 人為因素是土地覆被變化的主要原因,人們正采取越來越有力的手段影響著土地覆被變化。基于上述結論,文章提出了相應建議。

    土地覆被; 馬爾科夫模型; 土地利用程度; 空間變化動態(tài)度

    自從1986年國際科學理事會提出國際地圈與生物圈計劃(IGBP)和全球變化人文計劃(IHDP)于1995年共同擬定并發(fā)表了《土地利用/土地覆被變化科學研究計劃》以來,世界各國廣泛開展了對LUCC的研究。諸多研究均表明:土地質量退化、空氣污染、水質惡化等生態(tài)環(huán)境問題的出現總是與不合理的土地利用聯系在一起,而土地利用合理往往會促使區(qū)域生態(tài)環(huán)境向適應人類需求的良性方向發(fā)展[1]。國內外有很多基于區(qū)域尺度的LUCC研究取得了很多成果[2-5],其能夠在時間和空間上較好地反映出研究區(qū)域生態(tài)環(huán)境變化,為區(qū)域生態(tài)、經濟的可持續(xù)發(fā)展提供重要的理論依據[6],基于RS、GIS技術在LUCC方面研究的應用也頗多[7],吐魯番市十分特殊的自然環(huán)境與交通樞紐戰(zhàn)略地位、豐富的礦產資源和旅游資源使得研究吐魯番市土地覆被時間序列變化尤其重要。以前已有過研究吐魯番市土地利用或土地覆被的內容,但研究方法和側重點不同[8-10],這些研究內容較少地預測了吐魯番市未來土地覆被的發(fā)展變化,加之以前很多文獻應用馬爾科夫模型對模擬的步長考慮較少,只是單純的以1年或以抵達年為預測步長[11-13],這直接影響到了預測的精度,因此本文采用1990—2011年較長地數據序列和利用馬爾科夫模型求得的最優(yōu)模擬步長來研究和預測吐魯番市土地覆被變化,并分析吐魯番市土地利用程度和土地利用空間變化動態(tài)度,考慮到單一土地利用動態(tài)度和綜合土地利用動態(tài)度的缺點,引入單一土地覆被空間變化動態(tài)度和綜合土地覆被空間變化動態(tài)度[14]。

    1 研究區(qū)概況

    吐魯番市位于新疆維吾爾自治區(qū)東部,地處天山中東部主峰博格達山南麓,吐魯番盆地中心。東西寬90 km,南北長262 km,地勢南北高,中間低。地理坐標:東經88°29′28″—89°54′33″,北緯42°15′10″—43°35′,土地總面積13 589 km2。東臨哈密,西、南與巴音郭楞蒙古自治州的和靜、和碩、尉犁、若羌縣毗連,北隔天山與烏魯木齊市及昌吉回族自治州的奇臺、吉木薩爾、木壘縣相接。距新疆維吾爾自治區(qū)首府烏魯木齊市183 km。吐魯番屬于典型的大陸性暖溫帶荒漠氣候,日照充足,熱量豐富但又極端干燥,降雨稀少且大風頻繁,故有“火洲”、“風庫”之稱,自然環(huán)境十分特殊。國家“一主兩翼”大煤運鐵路線和西氣東輸二線的投入運營,交通樞紐戰(zhàn)略地位更加突出。旅游資源極為豐富,是享譽世界的優(yōu)秀旅游勝地。

    2 數據來源與研究方法

    2.1基于現狀土地利用圖和遙感影像的土地利用變化轉移矩陣

    采用獲得的1990年吐魯番市TM影像和2005年、2011年土地利用現狀圖,運用ENVI來對1990年的影像進行幾何校正、裁剪、輻射校正、波段合成、影像解譯等操作步驟,在影像解譯時依據《土地利用現狀調查技術規(guī)程》解譯為耕地、園地、有林地、草地、水體、建設用地、未利用地、裸地、沙地、鹽堿地、灌木林、坑塘水面、永久性冰雪、內陸灘涂14類,為了討論方便,經過在ArcGIS中進行后期處理,把相關地類融合為耕地、園地、林地、草地、建設用地、水體、未利用地7大類。經過影像解譯精度評估,Kappa系數為80.34%,解譯效果較好。把影像解譯結果轉為矢量格式數據導入到ArcGIS軟件中進一步對誤判的地物進行屬性修正,獲得1990年土地利用圖,結合2005年和2011年土地利用現狀圖,導出土地利用轉移矩陣(見表1、表2、表3)。

    表1 吐魯番市1990年、2005年、2011年土地覆被類型面積變化狀況 hm2

    注:表中()中數字表示為該地類占該年土地總面積的百分比。

    由表2和表3得出,在1990—2005年間草地和未利用地互轉面積較大,草地有33 518.59 hm2變?yōu)槲蠢玫兀蠢玫赜?3 612.19 hm2變?yōu)椴莸?,草地在此期間減少了4 965.8 hm2,主要是由于該地區(qū)氣候和缺水的原因引起的,未利用地的轉入和轉出面積都大于50 000 hm2,轉出與轉入面積分別為51 977.99,51 860.78 hm2。林地在此期間有所增加,主要是為防風固沙、水土保持的需要而植樹造林的結果。耕地與園地在此期間轉出的面積小于轉入的面積,面積都有所增加,分別為4 030.38,3 443.97 hm2,主要是由未利用地轉化而來,分別為6 127.18,5 416.78 hm2,說明耕地與園地的增加主要是由開荒而來。建設用地在這15 a間增加3 012.39 hm2,到2005年達10 412.60 hm2,主要是由未利用地、耕地、園地轉變而來,分別為2 499.62,1 012.34,964.89 hm2,在此期間由于土地整理項目和富民安居工程的實施等原因使得建設用地主要轉出為園地、未利用地和耕地,其比例分別為11.69%,8.2%,4.52%。

    水體在此期間是所有地類中減少面積最多的利用類型,達6 561.30 hm2,主要是因為北部博格達峰的永久性冰雪大面積減少,使得33.57%的水體變?yōu)榱瞬莸?,達到6 459.56 hm2,其次13.55%變?yōu)榱宋蠢玫兀_到2 607.77 hm2,這很可能與全球變暖的氣候變化密切相關。

    表2 吐魯番市1990-2005年土地利用轉移矩陣 hm2

    注:表中() 中的數字表示該行土地覆被類型轉變?yōu)樵摿型恋馗脖活愋偷陌俜謹?。下同?/p>

    表3 吐魯番市2005-2011年土地利用轉移矩陣 hm2

    表3中各地類的變化趨勢與表2的變化趨勢基本相同,只是從總體上看未利用地和耕地面積分別減少了18 646.83,766.37 hm2,耕地供需矛盾日益突出,隨著城市化進程的不斷加快,現有耕地已在較大程度上不能滿足日益增長的各項建設和農民增收的需要[15],而其他地類面積均有所增加,增幅最大的是草地和建設用地,分別為10 532.7,4 365.78 hm2,主要是2005—2011年有很多新建項目加之富民安居工程的實施,使得未利用地進一步開發(fā),并占用了近300 hm2耕地,還有就是在在近幾年實施“退糧還經”戰(zhàn)略的影響下[15],大量的耕地轉化成以葡萄種植為主的園地,同時對未利用地開墾也會導致園地面積的大幅增加。

    從表1—3可以得出這兩期各地類面積排在前4位的都是未利用地、草地、林地和園地。從減幅和增幅來看未利用地連續(xù)減少,連續(xù)增加的是建設用地、林地和園地,先增后減的是耕地,先減后增的是水體面積,但是從表1可以得出建設用地的增幅最為迅速。1990年、2005年和2011年耕地面積分別約16 506.26,20 536.64,19 770.26 hm2,呈現先增后減的趨勢,由表3可知主要是耕地15.09%轉變?yōu)閳@地的結果,因此土地管理部門應制定相應措施實現耕地的占補平衡。在1990—2005年15 a的變化面積是146 601.97 hm2,但2005—2011年變化面積就達到了50 521.6 hm2,可見吐魯番市的經濟迅速發(fā)展,人們對土地的利用程度不斷增加。

    由以上討論可以得出在2005年之前土地利用變化比較平緩,2005年之后土地利用變化比較劇烈,未利用地和草地互轉比例較大,永久性冰雪面積減少的面積超過6 000 hm2,建設用地面積增加迅速,說明全球氣候變化和人類活動是影響土地覆被變化的重要因素。

    2.2 吐魯番市的馬爾科夫模型

    3 基于不同步長模型模擬結果與分析

    由于馬爾科夫模型受到自然、經濟、政策等因素的影響使得預測結果會出現誤差,不同的模擬步長也會有不同的模擬精度[16],由以上討論中得出2005年后土地利用變化發(fā)生劇烈變化,因此以2011年為基準年來進行模擬是適當的,為了消除這些因素的影響,分別以1990年為基準年,2005年為抵達年分別以1年、3年、5年為模擬步長得出2005年預測結果,然后與2005年實際數據進行比較找出最優(yōu)的模擬步長,再以2011年為基準年,用2005—2011年轉移矩陣的轉移概率為轉移概率矩陣,比較得出最優(yōu)模擬步長來預測2011年后的土地覆被變化。

    3.1 初始狀態(tài)矩陣的建立

    以1990年的各土地覆被類型所占土地總面積的百分比作為初始狀態(tài)概率構成初始狀態(tài)矩陣(見表1)。

    3.2 轉移概率矩陣和最優(yōu)模擬步長的確定

    以該市1990—2005年的數據得到的轉移矩陣作為Markov模型中的轉移概率矩陣(見表2),分別以1年、3年和5年為模擬步長對吐魯番市土地覆被變化進行模擬,模擬結果見下表4。

    不同模擬步長與2005年土地利用實際面積比較誤差見表5。

    表4 Markov Chain模型的不同步長模擬下吐魯番市土地覆被變化 萬hm2

    表5 不同步長模擬下吐魯番市2005年土地利用模擬值與實際面積比較誤差

    由表4和表5可以得出,1990—2005年對于吐魯番市的土地覆被變化而言,模擬的步長越接近15年模擬的效果越好,即草地、耕地、建設用地、林地、水體、未利用地、園地的模擬值越接近2005年實測值。模擬步長1年、3年、5年的各個地類相對誤差均是1年>3年>5年,從相對誤差可以得出模擬的效果從高到底排列為:未利用地>草地>園地>林地>建設用地>耕地>水體,除對水體模擬的誤差達9%左右誤差較大外(這可能是解譯的誤差引起的與模型本身的精度或許沒有關系),其它地類的模擬誤差都在3%以下,可見用馬爾科夫模型預測吐魯番市土地覆被變化效果是較好的。

    由以上討論可知,對于吐魯番市而言,模擬步長越接近所達之年模擬精度越高,因此選擇以2011年為基準年,以n=6年為模擬步長來進行土地覆被變化預測,即以2005年各土地類型面積比例為初始概率矩陣(見表1),以2005—2011年的土地轉移頻數為土地利用轉移概率矩陣(見表3),預測結果見表6。

    表6 土地覆被類型預測結果 hm2

    由表6可以得出:(1)農用地中,耕地面積整體上逐年減少,且減幅較小,由2011年的19 770.26 hm2減少到2023年的19 218.64 hm2,說明耕地供需矛盾日益突出。園地和林地面積逐年增加分別從2011年的23 814.11,27 194.67 hm2增加到2023年的24 597.10,32 669.75 hm2。(2) 建設用地逐年增加,從2011年的14 778.38 hm2增加到2029年的21 698.14 hm2。(3) 水體面積也是逐年增加,且增幅較大,從2011年的13 302.42 hm2增加到2023年的14 467.39 hm2,水體面積在增加說明吐魯番市進行的一系列水利工程中起到了很好的作用。(4) 未利用地面積逐年減少,并且減少的面積特別大,由預測結果來看2011—2023年未利用地將減少33 970.23 hm2,這其中減少面積肯定有許多轉變?yōu)椴傻V用地,這與吐魯番市礦產資源豐富密切相關,也刺激了該地區(qū)經濟的發(fā)展。由此得出很可能是一系列土地整理和復墾項目的確起到了作用,使得吐魯番市水體、林地和草地的面積都有所增加。

    4 土地利用程度變化的時空特征分析

    4.1 吐魯番市土地利用程度分析

    土地利用程度主要反映了土地利用的廣度和深度,不僅反映了土地利用中土地本身的自然屬性,同時也反映了人類因素和自然環(huán)境因素的綜合效應[18]。采用劉紀遠[19]提出的土地利用程度的綜合分析方法,將土地利用程度按照土地自然綜合體在社會因素影響下的自然平衡分為4級,并賦予分級指數(見表7),得出吐魯番市土地利用程度定量表達式。

    表7 土地利用程度分級賦值

    土地利用程度綜合指數計算公式為[18]:

    式中:Lj——研究區(qū)域土地利用程度綜合指數;Ai——研究域內第i級土地利用程度分級指數;Ci——研究區(qū)域內第i級土地利用程度分級面積百分比;n——土地利用程度分級數。

    土地利用程度變化量Lb-a和變化率R,計算公式為:

    (1)

    (2)

    式中:Lb——b時間的區(qū)域土地利用程度綜合指數;La——a時間的區(qū)域土地利用程度綜合指數;Ai——第i級的土地利用程度分級指數;Cib——某區(qū)域b時間第i級土地利用程度面積百分比;Cia——某區(qū)域a時間第i級土地利用程度面積百分比。

    土地利用程度及其變化量和變化率可定量地揭示該范圍土地利用的綜合水平和變化趨勢,如果Lb-a>0或R>0,則區(qū)域土地利用處于發(fā)展期,否則處于調整期或衰退期[18]。根據表7和以上兩式得出吐魯番市不同時期土地利用程度綜合指數,土地利用程度變化量和變化率。

    經計算,可以得出1990—2011年土地利用程度綜合指數呈連續(xù)上升趨勢,從1990年的131.09上升到2011年的134.06;土地利用程度變化量從1990—2005年的0.997上升到2005—2011年的1.977,1990—2011年為2.974;土地利用程度變化率從1990—2005的0.007 6上升到2005—2011年的0.015 0,1990—2011年為0.022 7。說明人們對吐魯番的土地利用程度不斷加深,利用區(qū)域更加廣泛,主要是由于經濟發(fā)展而使得很多以前沒有利用的土地現在正在利用,以前利用的土地其利用的程度更深更細,期間許多建設項目落戶該地區(qū),城區(qū)面積和居民點面積也不斷增加,只是居民點越來越集中,為了更好地利用土地,有些村莊集體搬遷合為一整體。

    4.2 吐魯番市土地利用空間變化動態(tài)度

    存在單一土地利用動態(tài)度和綜合土地利用動態(tài)度,應用它們的缺點在于出現地類面積增減變化的抵消,這樣在一定程度上掩蓋了土地覆被變化的實際情況,不能夠反映出土地覆被類型的空間變化[14],因此引入單一土地覆被空間變化動態(tài)度和綜合土地覆被空間變化動態(tài)度。

    單一土地覆被空間變化動態(tài)度,其表達式為:

    (3)

    綜合土地覆被空間變化動態(tài)度,其表達式為:

    (4)

    單一土地覆被空間變化動態(tài)度和綜合土地覆蓋空間變化動態(tài)度體現了土地覆被變化在空間上的變化速度[14]。以1990—2005年土地利用轉移面積和2005—2011年土地利用轉移面積計算吐魯番市土地覆被空間變化動態(tài)度,見表8。

    表8 不同時段吐魯番市土地覆被空間變化動態(tài)度 %

    從表8可以得出:

    (1) 1990—2005年間空間動態(tài)變化度最大的是耕地和林地,其次是園地、建設用地、水體、草地、未利用地,分別為11.6%,11.58%,7.57%,6.42%,4.16%,1.85%,0.7%,在此期間耕地、林地、園地、建設用地空間變化動態(tài)度最大主要是因為期間人們墾荒和建設用地占用其它地類的結果。2005—2011年間土地覆蓋類型空間變化動態(tài)度最大的是建設用地,其次是耕地、園地、林地、水體、草地、未利用地,分別為11.33%、9.37%,7.96%,7.28%,1.86%,1.34%,0.62%,可見該時期隨著城市化水平的提高、城鎮(zhèn)建設的發(fā)展引起了建設用地空間變化動態(tài)度的加大。

    這兩期都是未利用地空間變化動態(tài)度為最小,這并不是說未利用地變化最小,從總體上看未利用地在這三期中每一期變化的面積都超過了30 000 hm2,均為各地類最高,主要是因為未利用地占到了吐魯番市總面積的70%以上,基數太大所以空間變化動態(tài)度相對最小,所以所有地類的空間變化動態(tài)度都是相對而言的。

    (2) 比較前后兩個階段,2005—2011年除了建設用地和園地空間變化動態(tài)度大于前一個時期外,其它地類均小于前一個時期,由此可以得出主要是由于建設用地和園地的變化造成了后一時期的土地利用綜合程度指數總的趨勢是逐年升高,說明這兩個時期建設用地和園地與其他地類相互之間的轉移比較活躍,這主要是由于園地種植的主要是葡萄,是一大經濟產業(yè),而使得許多其它地類轉變?yōu)榱藞@地,再者由于城市規(guī)劃、耕地占補平衡、經濟快速發(fā)展等原因使得許多其他地類轉變?yōu)榱私ㄔO用地。

    (3) 綜合前兩個階段,1990—2011年期間土地覆蓋類型空間變化動態(tài)度變化最大的是林地,其次是耕地、建設用地、園地、水體、草地、未利用地,其實耕地和水體主要反映的是1990—2005年時期的變化,主要是開荒和全球氣候變暖的結果;建設用地和園地主要反映的是2005—2011年時期的變化,說明隨著城鎮(zhèn)化水平的提高和經濟發(fā)展的需要,建設用地和園地空間變化動態(tài)度加大;而未利用地和林地主要反映的是1990—2011年整個時期的變化,期間由于環(huán)境保護的需要,許多未利用地轉變?yōu)榱肆值囟脕矸里L固沙、水土保持,未利用地與其他地類相互轉變十分頻繁。

    通過以上分析得出:人為因素是土地覆蓋變化的主要原因;看空間變化動態(tài)度時一定要看它的基數面積值,值最小的也可能是絕對數量變化最大的;人們正采取越來越有力的手段改變著土地利用結構。

    5 結論與建議

    5.1 結 論

    (1) 從面積上看,1990—2011年吐魯番市面積最大的是未利用地、草地、林地、園地,其它地類各年面積排名稍有不同。耕地面積2005—2011年面積減少了766.37 hm2,耕地供需矛盾日益突出。

    (2) 從面積轉移矩陣中可以得出,2005年之前土地覆被變化比較平緩,2005年之后土地覆被變化比較劇烈,1990—2005年15年間才變化146 601.97 hm2,而2005—2011年6年之間就變化了50 521.60 hm2;1990—2011年未利用地和草地互轉數量較大,未利用地轉變?yōu)槠渌仡惷娣e較大,其它地類轉變?yōu)槲蠢玫孛娣e較少;水體在1990—2005年間面積大幅減少近6 000 hm2,主要是由于北部博格達峰永久性冰雪的減少,這很可能與全球氣候變暖密切相關;建設用地面積增加迅速。

    (3) 對于吐魯番市而言,馬爾科夫模型的預測步長越接近于抵達年,那么對各個地類的預測誤差越小,除對水體模擬的誤差達9%左右,誤差較大外,其它地類的模擬誤差都在3%以下,不過這可能與對1990年影像解譯的誤差有關,而與模型本身的精度關系不大,可見用馬爾科夫模型模擬吐魯番市土地覆被變化效果是較好的。

    (4) 從預測結果來看,農用地中,耕地逐年減少,說明照現在的情況發(fā)展下去,離現在提出的保護耕地的目標是背道而馳的;園地逐年增加,這與為了發(fā)展經濟而使得許多耕地轉變?yōu)榱藞@地密切相關。建設用地逐年增加,從2011年的14 778.38 hm2增加到2023年的21 698.14 hm2。水體也是逐年增加且增幅較大,從2011年的13 302.42 hm2增加到2023年的14 467.39 hm2。未利用地逐年減少,并且減少的面積特別大,由預測結果來看,2011—2023年未利用地將減少33 970.23 hm2,這減少的面積肯定有很多轉變?yōu)榱瞬傻V用地,這與吐魯番市礦產資源豐富密切相關。

    (5) 從吐魯番市土地利用程度分析得出,1990—2011年間人們對土地利用程度不斷加深,土地利用程度綜合指數從1990年的131.09上到2011年的134.06。

    (6) 吐魯番市土地覆被空間變化動態(tài)度方面可知,1990—2005年期間耕地、林地、園地的空間變化動態(tài)度最大,主要是該期間墾荒、經濟發(fā)展和從生態(tài)方面需要而植樹造林的結果。2005—2011年間土地覆蓋類型空間變化動態(tài)度最大的是建設用地,可見該地區(qū)城市化水平不斷提高。主要是由于2005—2011年建設用地和園地的變化造成了1990—2011年的土地利用綜合程度指數總的趨勢是逐年升高的,說明經濟的發(fā)展是土地覆被變化的重要因素。耕地、水體主要反映的是1990—2005年的變化,主要是開荒和全球氣候變暖的結果;建設用地和園地主要反映的是2005—2011年時期的變化;未利用地和林地主要反映的是1990—2011年整個時期的變化,期間由于環(huán)境保護的需要許多未利用地轉變?yōu)榱肆值囟脕矸里L固沙、水土保持,未利用地與其他地類相互轉變頻繁。

    綜合以上結論可以得出人為因素是土地覆被變化的主要原因,人們正采取越來越有力的手段改變著土地利用結構。

    5.2 建 議

    嚴格制定保護耕地的措施,實現耕地的占補平衡;加大對未利用地的開發(fā),吐魯番市有近70%的未利用地,加上在未利用地上有很多旅游資源,可在生態(tài)旅游方面多做規(guī)劃,解決該地區(qū)一些土地資源過度開發(fā)和退化的問題;適當減少工礦企業(yè)的發(fā)展,多發(fā)展生態(tài)產業(yè),因為過多的工礦產業(yè)會增多CO2的排放,使得全球氣候變暖的效應更加突出,以遏制吐魯番市北部博格達峰永久性冰雪面積減少的速度。

    由于缺少2000年之前十分詳細的資料,因此根據相關資料對吐魯番市土地覆被變化的原因只是做了一些定性的解釋,如果資料充足,最好可對該地區(qū)土地覆被變化的驅動因素做一些定量的研究,如用主成分分析法、典型相關分析法等。這有待以后改進、完善;由于馬爾科夫模預測的對象要求其保持某一固定的轉移概率,要求預測對象的變化過程是平穩(wěn)的,并且不能綜合考慮政策的改變等因素的影響,但實際上土地覆被的變化會因諸多自然或人為因素的影響而隨時間有所變化,是該模型的局限所在,有待進一步完善。

    [1] 劉新衛(wèi),陳百明,史學正.國內LUCC 研究進展綜述[J].土壤,2004,36(2):132-135.

    [2] 葛全勝,戴君虎,何凡能,等.過去三百年中國土地利用變化與陸地碳收支[M].北京:科學出版社,2008.

    [3] 張云鵬,孫燕,王小麗,等.不同尺度下的土地利用變化驅動力研究:以常州市新北區(qū)為例[J].水土保持研究,2012,19(6):111-117.

    [4] 容芳芳,塔西甫拉提·特依拜,田源,等.于田綠洲土地利用/覆蓋變化軌跡分析[J].水土保持研究,2010,17(3):259-263.

    [5] 譚君,李世平.銅川市土地利用變化對生態(tài)系統(tǒng)服務價值的影響分析[J].水土保持研究,2012,19(6):131-136.

    [6] 余慧容,蒲春玲,劉志有,等.基于TM/ETM+綠洲城市土地利用時空演變分析:以新疆奎屯市為例[J].水土保持研究,2012,19(6):147-151.

    [7] Chandra P Giri. Remote Sensing of Land Use and Land Cover:Principles and Applications[M]. New York:CRC Press,2012.

    [8] 程克堅,彭補拙,濮勵杰.干旱綠洲地區(qū)土地資源可持續(xù)利用初探:以新疆吐魯番市為例[J].資源科學,1998,20(4):14-18.

    [9] 劉榮,高敏華,謝峰.基于Logistic回歸模型的土地利用格局模擬分析:以新疆吐魯番市為例[J].水土保持研究,2009,16(6):74-78.

    [10] 張光耀,劉光遠,崔麗娜,等.西北干旱區(qū)設施農業(yè)土地利用效益研究:以吐魯番市為例[J].新疆農業(yè)科學,2011,48(6):1157-1161.

    [11] 曹晶晶,吳靜,李純斌.基于馬爾柯夫模型的酒泉市肅州區(qū)LUCC趨勢預測[J].國土與自然資源研究,2008(1):45-47.

    [12] 劉啟承,熊文強,韓貴鋒.用馬爾可夫理論預測三峽庫區(qū)的土地利用趨勢[J].重慶大學學報,2005,28(2):107-110.

    [13] 徐嵐,趙界.利用馬爾柯夫過程預測東陵區(qū)土地利用格局的變化[J].應用生態(tài)學報,1993,4(3):272-277.

    [14] 劉福輝.基于遙感的張家界市永定區(qū)土地覆蓋動態(tài)監(jiān)測研究[D].長沙:中南大學,2009.

    [15] 侯芳,高敏華.吐魯番市土地利用動態(tài)變化研究[J].資源環(huán)境與工程,2010,24(3):317-320.

    [16] 冉圣宏,李秀彬,呂昌河.土地覆被及生態(tài)服務價值變化的多時間尺度模擬:以四川省漁子溪流域為例[J].地理學報,2006,61(10):1113-1120.

    [17] 劉姣娣,曹衛(wèi)彬,李華,等.基于馬爾科夫模型的石河子墾區(qū)土地利用動態(tài)變化研究[J].新疆農業(yè)科學,2011,48(10):1936-1940.

    [18] 李傳哲,于福亮,劉佳,等.近20年來黑河干流中游地區(qū)土地利用/覆被變化及驅動力定量研究[J].自然資源學報,2011,26(3):353-363.

    [19] 劉紀遠.西藏自治區(qū)土地利用[M].北京:科學出版社,1992.

    DynamicChangesofLandCoverBasedonMarkovModelinTurpanCity

    HE Bao-zhong1,2, GAO Min-hua1,2, ZHAO Jun-an1

    (1.CollegeofResourcesandEnvironmentalScience,XinjiangUniversity,Urumqi830046,China; 2.KeylaboratoryofOasisEcology,MinistryofEducation,Urumqi830046,China)

    The land-cover data obtained from the TM remote sensing images of 1990 as well as land-use map of 2005, 2011 were used in Turpan to set up a Markov model for predicting changes of land-cover. Given the accuracy of prediction in Markov model is closely related to its simulated step, therefore, firstly, 1990 was taken as initial year to predict various areas of land-cover on different classes in 2005, and the optimal simulated step was found out by comparing the minimum difference between the simulate and actual areas of land-cover. Secondly, the changes of land-cover in 2011 was predicted by taking the optimal simulated step in Markov model, then the land use degree and land use space change dynamic were analyzed based on this. Results could be summarized as: (1) the areas of transform between unutilized land and grassland are enormous; (2) the cultivated land decreased quickly since 2005, which indicated that the contradiction between supplement and demand was being increased, and also meant that the balance of the cultivated occupation and supplement had not implemented; (3) nearly 6 000 hectares of water disappeared in 21 years which was largely probably associated with global warming; (4) because the change of construction land and garden plot made the land use degree increase year by year, it exhibited that the development of economy was an important factor resulting in the changes of land cover; (5) anthropic factor was a major reason which led to the changes of land-cover, people were taking more and more powerful means to affect the land cover change. Therefore, some appropriate proposals were also offered according to the above results.

    land cover; Markov model; land use degree; dynamic space change

    2013-11-25

    :2013-12-26

    國家自然科學基金項目“吐魯番綠洲人與自然耦合系統(tǒng)變化的驅動與適應機制研究”(41271168)。

    何寶忠(1989—),男,重慶市梁平縣人,在讀碩士,主要從事遙感和土地利用研究。E-mail:hbz108@163.com

    高敏華(1963—),男,新疆烏魯木齊市人,副教授,碩士生導師,主要從事土地資源評價及土地利用規(guī)劃。E-mail:mh_gao@163.com

    F301.24

    :A

    :1005-3409(2014)05-0041-08

    猜你喜歡
    吐魯番市園地步長
    互動園地
    基于Armijo搜索步長的BFGS與DFP擬牛頓法的比較研究
    吐魯番市文聯脫貧攻堅巡禮采風活動
    吐魯番(2020年3期)2020-11-09 08:38:32
    吐魯番市文聯第二屆“交河雅集”文學藝術交流活動剪影
    吐魯番(2019年1期)2019-07-17 11:56:32
    吐魯番市欲打造伏領跑者基地
    藝術園地
    藝術園地
    書畫園地
    基于逐維改進的自適應步長布谷鳥搜索算法
    一種新型光伏系統(tǒng)MPPT變步長滯環(huán)比較P&O法
    電測與儀表(2014年2期)2014-04-04 09:04:00
    美女视频免费永久观看网站| 另类亚洲欧美激情| 中文字幕另类日韩欧美亚洲嫩草| 大香蕉久久成人网| 各种免费的搞黄视频| 黄色毛片三级朝国网站| 亚洲精品在线美女| 精品一品国产午夜福利视频| 亚洲欧美精品自产自拍| 性少妇av在线| 国产成人a∨麻豆精品| 秋霞在线观看毛片| 人妻 亚洲 视频| 欧美日韩综合久久久久久| 少妇被粗大猛烈的视频| 最近中文字幕高清免费大全6| 久久人人爽人人片av| 亚洲精品久久成人aⅴ小说| 免费观看人在逋| 亚洲伊人久久精品综合| 在线观看免费高清a一片| 久久久久久免费高清国产稀缺| 久久毛片免费看一区二区三区| 久久久久久久久久久久大奶| 又粗又硬又长又爽又黄的视频| 亚洲欧美清纯卡通| 亚洲精品一区蜜桃| 女性生殖器流出的白浆| 免费在线观看黄色视频的| 精品少妇久久久久久888优播| 久久久久国产精品人妻一区二区| 一区福利在线观看| 久久99一区二区三区| 中国三级夫妇交换| 美国免费a级毛片| 国产一级毛片在线| 国产一区二区三区综合在线观看| 国产精品一国产av| 天堂8中文在线网| 亚洲欧美日韩另类电影网站| 成人午夜精彩视频在线观看| 久久免费观看电影| 日本猛色少妇xxxxx猛交久久| 99久国产av精品国产电影| 无遮挡黄片免费观看| 最近中文字幕2019免费版| 国产精品秋霞免费鲁丝片| 看免费av毛片| 亚洲国产成人一精品久久久| 国产一区二区激情短视频 | 波多野结衣av一区二区av| 久久国产精品大桥未久av| 视频区图区小说| 高清不卡的av网站| 国产欧美日韩一区二区三区在线| 19禁男女啪啪无遮挡网站| 制服诱惑二区| 9色porny在线观看| 婷婷成人精品国产| 老汉色av国产亚洲站长工具| av福利片在线| 丁香六月欧美| 飞空精品影院首页| www.av在线官网国产| 日本色播在线视频| 国产99久久九九免费精品| 悠悠久久av| 亚洲精品美女久久av网站| 哪个播放器可以免费观看大片| 精品少妇久久久久久888优播| 老司机靠b影院| 亚洲在久久综合| 国产精品无大码| 黄色怎么调成土黄色| 婷婷成人精品国产| 五月开心婷婷网| 视频区图区小说| netflix在线观看网站| 亚洲av在线观看美女高潮| 免费高清在线观看日韩| 色婷婷av一区二区三区视频| 一级毛片我不卡| 日本91视频免费播放| 欧美日韩亚洲综合一区二区三区_| 亚洲av男天堂| 超碰成人久久| 日韩 欧美 亚洲 中文字幕| 精品福利永久在线观看| svipshipincom国产片| 一二三四在线观看免费中文在| 久热这里只有精品99| 亚洲激情五月婷婷啪啪| 欧美久久黑人一区二区| 蜜桃国产av成人99| 女人久久www免费人成看片| 一区二区av电影网| 亚洲综合精品二区| 日韩 欧美 亚洲 中文字幕| 母亲3免费完整高清在线观看| 国产精品 欧美亚洲| 少妇人妻精品综合一区二区| 国产伦理片在线播放av一区| 亚洲一卡2卡3卡4卡5卡精品中文| 美女福利国产在线| av在线老鸭窝| 纵有疾风起免费观看全集完整版| 国产亚洲最大av| 日韩中文字幕欧美一区二区 | 两性夫妻黄色片| 日本av手机在线免费观看| 男人舔女人的私密视频| 肉色欧美久久久久久久蜜桃| www.精华液| 一级黄片播放器| 亚洲,欧美,日韩| av福利片在线| 亚洲精品国产av成人精品| 久久国产精品男人的天堂亚洲| 永久免费av网站大全| 国产精品免费大片| 侵犯人妻中文字幕一二三四区| 宅男免费午夜| 美女主播在线视频| 久久久久国产一级毛片高清牌| 在线观看国产h片| 热re99久久精品国产66热6| 成年人免费黄色播放视频| 肉色欧美久久久久久久蜜桃| 国产成人免费无遮挡视频| 狂野欧美激情性xxxx| 久久ye,这里只有精品| 王馨瑶露胸无遮挡在线观看| 久久久久久久大尺度免费视频| 美女中出高潮动态图| 美女脱内裤让男人舔精品视频| 国产精品无大码| 亚洲国产中文字幕在线视频| 女人爽到高潮嗷嗷叫在线视频| 欧美成人精品欧美一级黄| 欧美变态另类bdsm刘玥| 黄频高清免费视频| 美女主播在线视频| 91老司机精品| 九色亚洲精品在线播放| 97在线人人人人妻| 亚洲av日韩精品久久久久久密 | 精品亚洲成国产av| 久久久久久久久久久免费av| 只有这里有精品99| 欧美日韩亚洲国产一区二区在线观看 | 亚洲精品久久久久久婷婷小说| 日韩伦理黄色片| xxx大片免费视频| 制服人妻中文乱码| 亚洲欧美清纯卡通| 国产亚洲一区二区精品| av线在线观看网站| 午夜免费男女啪啪视频观看| 欧美日韩一区二区视频在线观看视频在线| av免费观看日本| 综合色丁香网| 999精品在线视频| 国产成人一区二区在线| 如日韩欧美国产精品一区二区三区| 久久综合国产亚洲精品| 9热在线视频观看99| 久久人人97超碰香蕉20202| 日日啪夜夜爽| 亚洲成av片中文字幕在线观看| 国产高清不卡午夜福利| 一级毛片黄色毛片免费观看视频| 色精品久久人妻99蜜桃| 日韩伦理黄色片| 午夜老司机福利片| 国产亚洲欧美精品永久| av在线播放精品| 日韩中文字幕欧美一区二区 | 男人添女人高潮全过程视频| 精品国产乱码久久久久久男人| 亚洲欧美成人综合另类久久久| 黄色怎么调成土黄色| 久久久久久免费高清国产稀缺| 国产一卡二卡三卡精品 | av一本久久久久| 国产精品久久久久久精品电影小说| 男女无遮挡免费网站观看| 在线看a的网站| 亚洲国产欧美在线一区| 最近的中文字幕免费完整| 不卡视频在线观看欧美| 亚洲综合色网址| 黄网站色视频无遮挡免费观看| videosex国产| 黑人猛操日本美女一级片| 在线亚洲精品国产二区图片欧美| 永久免费av网站大全| 亚洲精品久久午夜乱码| 国产欧美日韩一区二区三区在线| 在线精品无人区一区二区三| 亚洲国产欧美在线一区| 超碰成人久久| 精品一品国产午夜福利视频| 亚洲欧美成人综合另类久久久| 97人妻天天添夜夜摸| 岛国毛片在线播放| 国产成人啪精品午夜网站| 亚洲国产成人一精品久久久| a级毛片在线看网站| av国产精品久久久久影院| 国产日韩欧美在线精品| 欧美人与性动交α欧美精品济南到| 午夜91福利影院| 少妇精品久久久久久久| 天堂中文最新版在线下载| 大码成人一级视频| 2021少妇久久久久久久久久久| 涩涩av久久男人的天堂| 成人亚洲精品一区在线观看| 成人国产麻豆网| 高清不卡的av网站| 亚洲国产成人一精品久久久| 成人免费观看视频高清| 国产麻豆69| 午夜91福利影院| 国产日韩欧美亚洲二区| 一级毛片黄色毛片免费观看视频| 性少妇av在线| 天天躁日日躁夜夜躁夜夜| 国产亚洲精品第一综合不卡| 欧美日韩一区二区视频在线观看视频在线| 国产成人欧美在线观看 | 日本欧美国产在线视频| 午夜免费男女啪啪视频观看| 午夜av观看不卡| 久久99热这里只频精品6学生| 成年av动漫网址| 国产有黄有色有爽视频| 亚洲综合色网址| 精品少妇一区二区三区视频日本电影 | 亚洲国产看品久久| 满18在线观看网站| 一边亲一边摸免费视频| 美女国产高潮福利片在线看| 男人添女人高潮全过程视频| 国产亚洲欧美精品永久| 我的亚洲天堂| 欧美 日韩 精品 国产| 亚洲av日韩在线播放| 啦啦啦视频在线资源免费观看| 一区福利在线观看| e午夜精品久久久久久久| 热re99久久精品国产66热6| 最近最新中文字幕大全免费视频 | 在线亚洲精品国产二区图片欧美| 1024视频免费在线观看| 亚洲综合精品二区| www.av在线官网国产| 激情五月婷婷亚洲| 免费黄频网站在线观看国产| 国产极品粉嫩免费观看在线| 国产精品国产三级国产专区5o| 一级毛片黄色毛片免费观看视频| 日日摸夜夜添夜夜爱| 赤兔流量卡办理| 亚洲精品乱久久久久久| 亚洲美女搞黄在线观看| 香蕉国产在线看| av女优亚洲男人天堂| 老司机影院毛片| 久久久久久久久免费视频了| 在线观看免费视频网站a站| 操美女的视频在线观看| 女人被躁到高潮嗷嗷叫费观| 亚洲av国产av综合av卡| a级毛片黄视频| 黑丝袜美女国产一区| 日韩成人av中文字幕在线观看| 久久狼人影院| 亚洲成国产人片在线观看| 中文字幕制服av| 日韩大码丰满熟妇| 国产 精品1| 人体艺术视频欧美日本| 久久 成人 亚洲| 1024香蕉在线观看| 男女边吃奶边做爰视频| tube8黄色片| 成人国产麻豆网| 超碰97精品在线观看| 国产精品 国内视频| 老司机深夜福利视频在线观看 | 欧美激情 高清一区二区三区| 成人三级做爰电影| 男女午夜视频在线观看| 日韩伦理黄色片| 在线观看www视频免费| 永久免费av网站大全| 一本久久精品| 人人妻人人澡人人爽人人夜夜| 黄片播放在线免费| 午夜福利乱码中文字幕| 色婷婷av一区二区三区视频| 一级a爱视频在线免费观看| 一本一本久久a久久精品综合妖精| 汤姆久久久久久久影院中文字幕| 国产精品女同一区二区软件| 成年av动漫网址| 国产成人精品久久二区二区91 | 热re99久久精品国产66热6| av.在线天堂| 在线观看国产h片| 熟妇人妻不卡中文字幕| 国产福利在线免费观看视频| 成人影院久久| 九九爱精品视频在线观看| 制服丝袜香蕉在线| 91老司机精品| av在线app专区| 欧美精品一区二区大全| 欧美日韩视频高清一区二区三区二| 纯流量卡能插随身wifi吗| 国产成人精品在线电影| 国产精品三级大全| 999精品在线视频| 成人亚洲欧美一区二区av| 美女扒开内裤让男人捅视频| 天堂俺去俺来也www色官网| 久久女婷五月综合色啪小说| 亚洲av福利一区| 韩国高清视频一区二区三区| 你懂的网址亚洲精品在线观看| 亚洲国产精品一区三区| 亚洲图色成人| 亚洲精品视频女| 一二三四中文在线观看免费高清| 亚洲第一av免费看| 亚洲一卡2卡3卡4卡5卡精品中文| 国产精品嫩草影院av在线观看| 国产在线视频一区二区| 色综合欧美亚洲国产小说| 国产精品久久久久久精品古装| 肉色欧美久久久久久久蜜桃| 国产人伦9x9x在线观看| 99久久人妻综合| 久久久亚洲精品成人影院| 国产亚洲欧美精品永久| 黑人欧美特级aaaaaa片| 欧美国产精品va在线观看不卡| 精品少妇黑人巨大在线播放| 国产男人的电影天堂91| 亚洲精品国产av蜜桃| 亚洲一码二码三码区别大吗| 观看美女的网站| 少妇精品久久久久久久| 观看av在线不卡| 免费在线观看视频国产中文字幕亚洲 | 十八禁人妻一区二区| 中文字幕高清在线视频| 国产免费又黄又爽又色| 亚洲美女黄色视频免费看| 精品久久久精品久久久| 亚洲视频免费观看视频| 极品人妻少妇av视频| 一级a爱视频在线免费观看| 成人免费观看视频高清| 久久久久精品国产欧美久久久 | bbb黄色大片| 97精品久久久久久久久久精品| 国产男女超爽视频在线观看| 制服人妻中文乱码| 国产野战对白在线观看| 黄片无遮挡物在线观看| 国产精品久久久久久精品电影小说| 国产福利在线免费观看视频| 水蜜桃什么品种好| 日韩大码丰满熟妇| 亚洲av综合色区一区| 亚洲熟女精品中文字幕| 欧美 日韩 精品 国产| av福利片在线| 国产在线一区二区三区精| 国产麻豆69| 人人妻,人人澡人人爽秒播 | 成人18禁高潮啪啪吃奶动态图| 亚洲成人国产一区在线观看 | 久久久久久免费高清国产稀缺| 国产精品成人在线| 欧美日韩亚洲高清精品| 日本一区二区免费在线视频| 亚洲av在线观看美女高潮| 免费不卡黄色视频| 一级爰片在线观看| 高清av免费在线| 久热这里只有精品99| 国产片内射在线| 1024香蕉在线观看| 久久青草综合色| 久久久久久久精品精品| 五月开心婷婷网| 国产精品女同一区二区软件| 日韩 欧美 亚洲 中文字幕| 国产精品秋霞免费鲁丝片| 少妇人妻精品综合一区二区| 精品国产乱码久久久久久男人| 一区二区三区四区激情视频| 80岁老熟妇乱子伦牲交| 亚洲精品乱久久久久久| 欧美成人精品欧美一级黄| 亚洲国产欧美一区二区综合| 精品卡一卡二卡四卡免费| 亚洲第一青青草原| 又粗又硬又长又爽又黄的视频| 九九爱精品视频在线观看| 久久性视频一级片| 超碰成人久久| 亚洲av电影在线进入| 久久久久人妻精品一区果冻| 欧美在线黄色| 日本一区二区免费在线视频| 青春草国产在线视频| 满18在线观看网站| 9191精品国产免费久久| 中文字幕人妻丝袜一区二区 | 国产精品香港三级国产av潘金莲 | 777米奇影视久久| 老司机靠b影院| 最新的欧美精品一区二区| 婷婷色综合www| 国产欧美亚洲国产| www.熟女人妻精品国产| 国产精品久久久av美女十八| 亚洲av成人精品一二三区| 啦啦啦啦在线视频资源| 母亲3免费完整高清在线观看| 亚洲国产最新在线播放| 亚洲精品美女久久av网站| av天堂久久9| 午夜福利乱码中文字幕| 最近最新中文字幕免费大全7| 中文字幕亚洲精品专区| 青草久久国产| 亚洲欧美日韩另类电影网站| 男的添女的下面高潮视频| 老司机深夜福利视频在线观看 | 国产精品国产三级国产专区5o| 国产成人精品久久二区二区91 | 男女午夜视频在线观看| 别揉我奶头~嗯~啊~动态视频 | 日韩中文字幕视频在线看片| 亚洲图色成人| 日本午夜av视频| 国产亚洲一区二区精品| 日韩一本色道免费dvd| 最近手机中文字幕大全| 久久久国产精品麻豆| 亚洲 欧美一区二区三区| 欧美黄色片欧美黄色片| 国产 一区精品| 婷婷色av中文字幕| av片东京热男人的天堂| 秋霞伦理黄片| 天美传媒精品一区二区| 黑丝袜美女国产一区| 美女国产高潮福利片在线看| 欧美日韩成人在线一区二区| 18禁观看日本| 精品午夜福利在线看| kizo精华| 99久久综合免费| 男女之事视频高清在线观看 | 少妇猛男粗大的猛烈进出视频| 欧美日本中文国产一区发布| 国产人伦9x9x在线观看| 午夜免费观看性视频| 啦啦啦在线观看免费高清www| 午夜影院在线不卡| 中国三级夫妇交换| 精品亚洲成a人片在线观看| 久久精品久久久久久噜噜老黄| 美女视频免费永久观看网站| 国产精品无大码| 一本久久精品| 久久精品国产亚洲av涩爱| 两个人看的免费小视频| 少妇被粗大的猛进出69影院| 国产亚洲欧美精品永久| 亚洲国产看品久久| 午夜激情久久久久久久| 欧美97在线视频| av网站在线播放免费| 欧美在线黄色| 午夜激情久久久久久久| 亚洲国产看品久久| 亚洲国产中文字幕在线视频| 成年动漫av网址| 综合色丁香网| 亚洲第一区二区三区不卡| 亚洲五月色婷婷综合| 观看av在线不卡| 午夜影院在线不卡| 综合色丁香网| 天堂中文最新版在线下载| 亚洲熟女毛片儿| 午夜老司机福利片| 两个人免费观看高清视频| 精品少妇久久久久久888优播| 一级毛片黄色毛片免费观看视频| 黄网站色视频无遮挡免费观看| 精品久久久久久电影网| 日韩一区二区视频免费看| xxxhd国产人妻xxx| 大话2 男鬼变身卡| 少妇 在线观看| 中文欧美无线码| 婷婷色av中文字幕| 欧美日韩一级在线毛片| 18禁动态无遮挡网站| 日韩伦理黄色片| 亚洲美女搞黄在线观看| 亚洲国产av新网站| 99久久人妻综合| 日本vs欧美在线观看视频| 亚洲精品一区蜜桃| 最新在线观看一区二区三区 | 国产黄色视频一区二区在线观看| 国产精品 国内视频| 免费黄色在线免费观看| av在线观看视频网站免费| 激情视频va一区二区三区| 欧美亚洲 丝袜 人妻 在线| 国产乱来视频区| 男女边摸边吃奶| 亚洲欧美激情在线| 国产黄频视频在线观看| 国产一级毛片在线| 日本av免费视频播放| 99久久人妻综合| 亚洲精品国产色婷婷电影| 免费高清在线观看日韩| 一边摸一边做爽爽视频免费| 少妇人妻久久综合中文| 黄片小视频在线播放| 大片免费播放器 马上看| 国产精品 国内视频| 嫩草影视91久久| www.熟女人妻精品国产| 大陆偷拍与自拍| 国产精品国产三级国产专区5o| 七月丁香在线播放| 亚洲成av片中文字幕在线观看| 不卡av一区二区三区| 国产熟女午夜一区二区三区| 精品视频人人做人人爽| 美女福利国产在线| 这个男人来自地球电影免费观看 | 丰满乱子伦码专区| 免费高清在线观看日韩| 建设人人有责人人尽责人人享有的| 日韩大片免费观看网站| 久久精品国产亚洲av涩爱| 久久久久网色| 日本一区二区免费在线视频| 国产精品久久久人人做人人爽| 9191精品国产免费久久| 91国产中文字幕| 国产xxxxx性猛交| 黄色毛片三级朝国网站| 视频在线观看一区二区三区| 精品亚洲乱码少妇综合久久| 亚洲国产精品一区三区| 久久精品国产亚洲av高清一级| 亚洲国产成人一精品久久久| svipshipincom国产片| 亚洲成人av在线免费| 麻豆精品久久久久久蜜桃| 在线观看免费午夜福利视频| 亚洲人成网站在线观看播放| 色播在线永久视频| 亚洲激情五月婷婷啪啪| av线在线观看网站| 久久人人爽av亚洲精品天堂| 亚洲一码二码三码区别大吗| 九九爱精品视频在线观看| av天堂久久9| 亚洲第一区二区三区不卡| 午夜精品国产一区二区电影| 国产视频首页在线观看| av.在线天堂| 亚洲精品国产av蜜桃| 久久ye,这里只有精品| 男女之事视频高清在线观看 | 精品视频人人做人人爽| 亚洲国产欧美网| 国产精品无大码| 一区福利在线观看| 欧美日韩国产mv在线观看视频| 午夜av观看不卡| 精品人妻在线不人妻| av国产久精品久网站免费入址| 国产精品偷伦视频观看了| 亚洲国产看品久久| 99久久99久久久精品蜜桃| 两性夫妻黄色片| 亚洲国产看品久久| 好男人视频免费观看在线| 中文欧美无线码| 午夜免费观看性视频| 亚洲精品国产av成人精品| 在线观看免费午夜福利视频| 18禁国产床啪视频网站| 国产精品久久久久成人av| 成人国产麻豆网| 欧美黑人欧美精品刺激| 黄色视频在线播放观看不卡| 中国三级夫妇交换|