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

    鄂爾多斯高原泊江海子流域土地利用動態(tài)變化及驅(qū)動因素分析*

    2017-11-06 09:54:47閆國振婁華君
    關(guān)鍵詞:江海覆蓋度土地利用

    閆國振, 張 征, 梁 康, 婁華君

    (1.北京林業(yè)大學(xué)林學(xué)院 北京 100083; 2.中國科學(xué)院地理科學(xué)與資源研究所陸地水循環(huán)及地表過程重點(diǎn)實(shí)驗(yàn)室北京 100101)

    鄂爾多斯高原泊江海子流域土地利用動態(tài)變化及驅(qū)動因素分析*

    閆國振1,2, 張 征1, 梁 康2**, 婁華君2

    (1.北京林業(yè)大學(xué)林學(xué)院 北京 100083; 2.中國科學(xué)院地理科學(xué)與資源研究所陸地水循環(huán)及地表過程重點(diǎn)實(shí)驗(yàn)室北京 100101)

    土地作為陸地生態(tài)系統(tǒng)的載體, 其利用方式的變化直接影響著陸地生態(tài)系統(tǒng)服務(wù)功能的種類和強(qiáng)度。泊江海子流域位于鄂爾多斯高原干旱半干旱區(qū), 生態(tài)系統(tǒng)脆弱, 流域生態(tài)系統(tǒng)服務(wù)功能在近十幾年嚴(yán)重退化。本文基于Landsat影像數(shù)據(jù)提取了研究區(qū)1990—2015年間的6期土地利用數(shù)據(jù), 分析研究區(qū)土地利用時(shí)空變化特征及其驅(qū)動因素。結(jié)果顯示: 研究期內(nèi)主要的土地利用類型為草地、耕地和未利用地, 土地利用總體變化顯著??臻g變化上呈東西部及中北部耕地增加、南部未利用地減少的主要趨勢。各土地利用類型面積變化及轉(zhuǎn)換特征不一。草地整體處于動態(tài)轉(zhuǎn)換平衡狀態(tài), 約占流域總面積的 70%。耕地和未利用地分別呈增加和減少的趨勢, 變化面積分別占流域總面積的6.05%和7.62%, 其各自的主要來源或去向均為草地。水域面積減少了總面積的1.45%, 主要轉(zhuǎn)出去向?yàn)闉┑亍?995—2000年是土地利用轉(zhuǎn)換最為頻繁的時(shí)段, 2005—2010年間土地利用變化呈動態(tài)不平衡狀態(tài)的表現(xiàn)最為突出。研究期內(nèi)流域土地利用變化受降雨減少、氣溫升高等氣候變化和生態(tài)恢復(fù)工程建設(shè)、水利工程建設(shè)等人類活動的綜合驅(qū)動影響。合理規(guī)劃土地利用, 限制種植高耗水樹種及作物, 減少河流上游匯水區(qū)耕地等, 是保障研究區(qū)生態(tài)系統(tǒng)可持續(xù)健康發(fā)展的重要措施。

    泊江海子流域; 土地利用變化; 趨勢與狀態(tài)指數(shù); 驅(qū)動因素; 生態(tài)系統(tǒng)服務(wù)功能; 鄂爾多斯遺鷗國家級自然保護(hù)區(qū)

    土地作為陸地生態(tài)系統(tǒng)的載體, 其利用方式的變化引發(fā)各類陸地生態(tài)系統(tǒng)在類型、面積及空間分布格局上的改變, 進(jìn)而使得陸地生態(tài)系統(tǒng)服務(wù)功能的種類和強(qiáng)度受到影響[1-2]。同時(shí), 土地利用變化在加劇流域水循環(huán)、改變流域水量平衡等方面對流域水資源分布有著重要影響[3-4]。而在驅(qū)動機(jī)制上, 土地利用變化主要受到氣候變化和人類活動的綜合作用[5-7]。土地利用變化和氣候變化間的聯(lián)系十分緊密, 兩者相互作用、相互影響[8-9]。人類活動的劇烈干擾導(dǎo)致近些年來全球及區(qū)域尺度上的土地利用變化顯著, 并對地表下墊面形態(tài)及水文過程產(chǎn)生明顯影響[10]。

    研究區(qū)泊江海子流域位于鄂爾多斯黃河中游干旱半干旱區(qū)過渡帶, 區(qū)域生態(tài)環(huán)境脆弱, 水土流失嚴(yán)重且水資源貧乏[11]。近10多年來, 在氣候變化的背景下, 大規(guī)模生態(tài)恢復(fù)措施的實(shí)施, 如天然林保護(hù)工程, 退耕還林/還草工程等, 黃河中游地區(qū)土地利用格局變化顯著[10,12]。羅婭等[13]對黃河中游地區(qū)1998—2010年間土地利用變化的研究顯示, 生態(tài)恢復(fù)工程的有效實(shí)施使研究區(qū)內(nèi)農(nóng)用地明顯減少, 林草面積明顯增加。劉曉林等[14]研究指出 20世紀(jì) 80年代至2010年黃河內(nèi)蒙古段十大孔兌流域的土地利用變化呈現(xiàn)出草地面積逐漸增加、未利用地逐漸減少的趨勢, 且1998年以后土地利用變化幅度和速度較前期有變快的趨勢。泊江海子流域最重要的生態(tài)服務(wù)功能之一是保護(hù)流域內(nèi)的珍稀物種鄂爾多斯遺鷗(Larus relictus)。泊江海子湖泊水面中的湖心島即是遺鷗的繁殖棲息地, 但在氣候變化和人類活動干擾下, 湖泊中的湖心島在近些年間因湖泊水位下降和湖面萎縮而使其可供遺鷗繁殖棲息的生態(tài)功能不斷退化甚至消失。何芬奇等[15]的研究顯示, 2004年之后未再發(fā)現(xiàn)遺鷗在此繁殖棲息。泊江海子湖泊的萎縮在很大程度上反映出整個(gè)流域陸地生態(tài)系統(tǒng)及其服務(wù)功能的嚴(yán)重退化, 這已然成為流域當(dāng)前面臨的最主要的生態(tài)問題。泊江海子流域?qū)儆谌辟Y料地區(qū), 當(dāng)前對泊江海子流域的研究, 主要集中在氣候環(huán)境變遷[16-18]、生物多樣性[19-20]和水資源特征[21-25]等方面, 針對泊江海子流域土地利用格局變化以及驅(qū)動因素等方面的研究比較薄弱。在土地利用變化方面, 僅有梁犁麗等[23]在研究區(qū)植被-水資源模擬研究中對1996—2005年間的不同植被類型面積進(jìn)行了統(tǒng)計(jì), 另有部分研究對泊江海子湖區(qū)及周邊的植被群落[26]及植被覆蓋方式的變化[24,27]進(jìn)行過簡單的描述。綜合考慮研究區(qū)所面臨的突出生態(tài)問題以及研究現(xiàn)狀, 同時(shí)結(jié)合研究區(qū)為閉流盆地的地理特征以及生態(tài)恢復(fù)和水土保持工程等國家和地方政策的影響,本文基于30 m空間分辨率的Landsat遙感影像數(shù)據(jù)對研究區(qū)泊江海子流域1990—2015年間的土地利用進(jìn)行分類提取, 并對其時(shí)空變化特征及其驅(qū)動因素進(jìn)行針對性的系統(tǒng)分析, 旨在一定程度上完善研究區(qū)的土地利用變化方面的科學(xué)基礎(chǔ)數(shù)據(jù), 為泊江海子流域陸地生態(tài)系統(tǒng)服務(wù)功能改善、流域水資源合理利用以及遺鷗自然保護(hù)區(qū)的保護(hù)與管理提供一定的科學(xué)參考。

    1 材料與方法

    1.1 研究區(qū)概況

    泊江海子流域地處內(nèi)蒙古自治區(qū)鄂爾多斯市東勝區(qū)以西約 45 km 處(圖 1), 地理位置介于109°6′30″~109°32′50″E 和 39°41′40″~39°56′2″N。流域位于鄂爾多斯波狀高原區(qū), 海拔在1 360~1 600 m,屬溫帶干旱半干旱大陸性氣候, 降雨主要集中在每年的7月和8月, 約占全年降水量的65%。流域年均氣溫約 5.2 ℃, 多年平均降水量約為 345 mm, 蒸發(fā)量約為2 524 mm。草地是流域內(nèi)最主要的植被類型。泊江海子流域?qū)偎闹芨?、中間低的閉合盆地, 流域面積約745 km2。鄂爾多斯遺鷗國家級自然保護(hù)區(qū)位于流域中部, 面積約156 km2。泊江海子湖泊位于流域和保護(hù)區(qū)的中心, 是流域內(nèi)最大的湖泊, 水面面積約 10 km2。雞溝河和烏爾圖河是流域內(nèi)最大的兩條季節(jié)性河流, 分別從西側(cè)和東側(cè)匯集流入泊江海子湖泊。

    圖1 研究區(qū)位置及概略圖Fig.1 Location and sketch map of the study area of the Bojiang Lake Baisn

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

    1.2.1 土地利用數(shù)據(jù)

    土地利用數(shù)據(jù)來源于 Landsat遙感影像數(shù)據(jù)(下載地址: https://earthexplorer.usgs.gov/)的人工解譯,數(shù)據(jù)空間分辨率為30 m。分別對研究區(qū)1990年、1995年、2000年、2005年、2010年和2015年共6期的土地利用類型數(shù)據(jù)進(jìn)行了解譯。在綜合研究區(qū)氣候特征、土地利用特點(diǎn)及遙感影像數(shù)據(jù)質(zhì)量的基礎(chǔ)上, 各期影像數(shù)據(jù)的時(shí)間見表1, 表中的平均云量是指本文研究區(qū)域的云量。

    表1 用于土地利用分類數(shù)據(jù)提取的遙感影像信息Table 1 Remote sensing image information used for the extraction of land use type

    1.2.2 氣象數(shù)據(jù)

    利用中國氣象數(shù)據(jù)網(wǎng)(http://data.cma.cn/)提供的降雨和氣溫?cái)?shù)據(jù)對研究區(qū)土地利用變化的氣候驅(qū)動進(jìn)行了分析。本文所利用的氣象數(shù)據(jù)為鄂爾多斯市東勝氣象站1990—2015年的降雨和氣溫月數(shù)據(jù), 經(jīng)求和平均后獲得年值, 進(jìn)而分析研究區(qū)的降水和氣溫變化。

    1.3 研究方法

    1.3.1 土地利用類型數(shù)據(jù)獲取

    利用Landsat TM和OLI遙感影像數(shù)據(jù), 經(jīng)過幾何校正、大氣校正等圖像預(yù)處理, 在ENVI 5.2軟件中采用監(jiān)督分類中的最大似然法對研究區(qū)1990—2015年間的 6期土地利用類型進(jìn)行分類提取。在參考我國《土地利用現(xiàn)狀分類》標(biāo)準(zhǔn)的基礎(chǔ)上, 結(jié)合研究區(qū)土地利用特征及本文研究目的, 將研究區(qū)土地利用類型分為耕地、林地、高覆蓋度草地、中覆蓋度草地、低覆蓋度草地、水域、灘地、建設(shè)用地及未利用地共 9類。根據(jù)筆者前期對研究區(qū)的野外實(shí)地調(diào)查記錄對本次土地利用分類提取的結(jié)果進(jìn)行了精度驗(yàn)證, 9種類別共75個(gè)樣本點(diǎn), 驗(yàn)證精度為93.33%, 本次土地利用分類提取數(shù)據(jù)可信。

    1.3.2 土地利用變化分析

    運(yùn)用ArcGIS 10.2軟件的空間分析模塊, 對6期土地利用數(shù)據(jù)進(jìn)行統(tǒng)計(jì)和疊加分析, 得到各期之間的土地利用類型面積轉(zhuǎn)移矩陣, 在此基礎(chǔ)上采用趨勢與狀態(tài)指數(shù)、凈變化速度、總趨勢與狀態(tài)指數(shù)和總凈變化速度分析泊江海子流域26年間土地利用類型變化的過程和趨勢。主要指數(shù)的計(jì)算方法如下:

    1)趨勢與狀態(tài)指數(shù):

    其中,

    式中:Ps為某研究時(shí)段內(nèi)各土地利用類型變化的趨勢與狀態(tài)指數(shù),Nc和Tc分別為各土地利用類型的面積凈變化和總變化,Ua、Ub分別為研究初期和末期某一土地利用類型的面積, ΔUin為時(shí)段內(nèi)其他土地利用類型轉(zhuǎn)變?yōu)槟惩恋乩妙愋偷拿娣e之和, ΔUout為同時(shí)段某一土地利用類型轉(zhuǎn)變?yōu)槠渌恋乩妙愋偷拿娣e之和。

    2)凈變化速度:

    式中:Rs為某研究時(shí)段內(nèi)各土地利用類型的凈變化速度,T為研究時(shí)段。

    3)總趨勢與狀態(tài)指數(shù):

    其中,

    式中:Pt為某研究時(shí)段內(nèi)所有土地利用類型變化的總趨勢與狀態(tài)指數(shù),Sa和Ss分別為所有土地利用類型的面積凈變化和總變化,Uai、Ubi分別為研究初期和末期第i種土地利用類型的面積, ΔUout–i為研究時(shí)段內(nèi)第i種土地利用類型轉(zhuǎn)變?yōu)槠渌愋偷拿娣e之和, ΔUin–i為其他土地類型轉(zhuǎn)變?yōu)榈趇種類型的面積之和,n為土地型總數(shù)。

    4)總凈變化速度:

    式中:Rt為某研究時(shí)段內(nèi)所有土地利用類型的總凈變化速度;T為研究時(shí)段。

    1.3.3 氣象數(shù)據(jù)分析

    在提取研究區(qū)年均降雨和氣溫?cái)?shù)據(jù)的基礎(chǔ)上,利用累積距平法對研究區(qū)研究期內(nèi)的年均降雨和氣溫變化特征和趨勢進(jìn)行了分析。累積距平的計(jì)算方法如下:

    式中: ADi為研究期內(nèi)某年的年均氣溫或降水量的累積距平值,xi為某年的年均氣溫或者降水量,為研究期內(nèi)的多年平均氣溫或降水量。

    2 結(jié)果與分析

    2.1 不同土地利用類型的變化特征

    1990—2015年泊江海子流域不同時(shí)期土地利用空間分布和變化分別見圖2和圖3。

    圖 2顯示, 流域最主要的土地利用類型為草地,其次為耕地和未利用地。圖3顯示1990—1995年間流域土地利用空間變化主要發(fā)生在西部和南部,1995—2000年間流域內(nèi)土地利用空間變化分布廣泛,2000—2005年間流域土地利用空間變化主要發(fā)生在中部, 2005—2010年間流域土地利用空間變化主要發(fā)生在西部和中北部, 2010—2015年間流域土地利用空間變化主要發(fā)生在流域中部。各時(shí)段發(fā)生土地利用變化的面積占流域總面積的比例依次分別為54.90%、90.43%、12.08%、15.76%和 9.65%。在1990—2015年間流域土地利用空間變化分布廣泛, 變化面積占流域總面積的66.67%。綜合圖2和圖3, 研究期內(nèi)流域西部區(qū)域呈耕地面積增加, 高覆蓋度草地減少趨勢; 流域東部區(qū)域耕地面積增加, 高覆蓋度草地和未利用地面積減少; 流域中北部區(qū)域耕地和建設(shè)用地面積增加、高覆蓋度草地和林地面積減少;流域南部區(qū)域則呈林地和低覆蓋度草地增加、高覆蓋度草地和未利用減少的趨勢。

    1990—2015年間流域各土地利用類型在流域中的面積占比變化見表 2。研究期內(nèi)草地面積較為穩(wěn)定, 約占流域總面積的70%。其中, 高覆蓋度草地整體呈減少趨勢, 26年間面積占比減少7.59%; 中、低覆蓋度草地呈增長趨勢, 26年間面積占比分別增長2.95%和5.23%。耕地面積占比在26年間增長6.05%,未利用地面積占比則減少 7.62%。林地面積占比微弱增長, 水域面積占比減少 1.45%, 灘地面積占比增長1.59%, 建設(shè)用地面積占比增長0.49%。

    利用趨勢與狀態(tài)指數(shù)和凈變化速度分析了泊江海子流域不同土地利用類型在研究期不同時(shí)段內(nèi)的變化情況, 結(jié)果見表 3。1990—1995年間流域內(nèi)面積變化最大的土地利用類型為高覆蓋度草地, 增加了31.96 km2, 主要來源于中覆蓋度草地、耕地和林地。水域是此時(shí)期內(nèi)變化程度最強(qiáng)的土地利用類型,其趨勢與狀態(tài)指數(shù)達(dá) 0.73, 面積增加的來源主要為未利用地。同期內(nèi)耕地和林地的變化程度也較突出,其趨勢與狀態(tài)指數(shù)均為-0.49, 面積減少的轉(zhuǎn)出去向主要為高覆蓋度草地。灘地在此時(shí)期內(nèi)的凈變化速度最快, 且處于動態(tài)轉(zhuǎn)出狀態(tài), 轉(zhuǎn)出去向主要為低覆蓋度草地和未利用地。在1995—2000年的研究時(shí)段內(nèi), 高覆蓋度草地面積減少81.63 km2, 轉(zhuǎn)出去向主要為中覆蓋度草地和耕地。此研究時(shí)段內(nèi)的耕地、水域和灘地的土地利用變化程度也較強(qiáng)烈, 趨勢與狀態(tài)指數(shù)分別為0.71、-0.81和0.97。耕地和灘地呈面積增加的趨勢, 耕地的轉(zhuǎn)入來源主要為高、中覆蓋度草地, 灘地的轉(zhuǎn)入來源主要為未利用地和草地。水域面積減少的轉(zhuǎn)出去向主要為未利用地。2000—2005年間水域和灘地兩種土地利用類型在面積、趨勢與狀態(tài)指數(shù)及凈變化速度上均較突出。由于在此時(shí)段內(nèi)泊江海子湖泊發(fā)生萎縮, 流域水域面積減少, 其主要轉(zhuǎn)出去向?yàn)闉┑?。未利用地?005—2010年間變化最突出的土地利用類型, 其面積減少 26.05 km2, 凈變化速度為 9.20%, 主要轉(zhuǎn)換為低覆蓋度草地。在2010—2015年間變化最為突出的土地利用類型為建設(shè)用地(Ps=0.91,Rs=10.89%),其面積增加的主要來源為中、低覆蓋度草地。

    1990—2015年間流域內(nèi)各土地利用類型間互相轉(zhuǎn)換明顯。耕地、高覆蓋度草地、未利用地、水域、灘地及建設(shè)用地在 26年間的趨勢與狀態(tài)指數(shù)均大于 0.5。水域(Rs=24.82%)和建設(shè)用地(Rs=20.19%)的凈變化速度較快, 26年間水域面積處于減少狀態(tài),轉(zhuǎn)出去向主要為灘地、低覆蓋度草地和未利用地;建設(shè)用地則呈現(xiàn)明顯的增加趨勢, 其轉(zhuǎn)入來源主要為草地。耕地、未利用地和灘地在研究期內(nèi)的凈變化速度也均超過10%。耕地和灘地面積呈增加趨勢, 其主要轉(zhuǎn)入來源分別為草地和水域。未利用地面積處于減少狀態(tài), 其轉(zhuǎn)出去向主要為中、低覆蓋度草地。林地、中覆蓋度草地和低覆蓋度草地在研究期內(nèi)也處于面積增加的狀態(tài), 高覆蓋度草地則呈面積減少趨勢。林地的轉(zhuǎn)入來源主要為高覆蓋度草地, 中覆蓋度草地的轉(zhuǎn)入來源主要為高覆蓋草地和未利用地, 低覆蓋度草地的轉(zhuǎn)入來源主要為未利用地, 高覆蓋度草地的轉(zhuǎn)出去向主要為耕地和中覆蓋度草地。

    2.2 土地利用類型綜合變化特征

    利用總趨勢與狀態(tài)指數(shù)和總凈變化速度分析了泊江海子流域土地利用類型綜合變化特征, 結(jié)果見表4。泊江海子流域在1990—2015年總時(shí)間段內(nèi)土地利用類型總體面積凈變化(Sa=16.66%)和總凈變化速度(Rt=3.13%)兩個(gè)指標(biāo)均高于各子研究時(shí)段。這表明泊江海子流域土地利用類型在 25年間變化明顯。1995—2000年的總變化指數(shù)超過45.22%, 表明流域土地利用變化在此研究時(shí)段內(nèi)動態(tài)轉(zhuǎn)換頻繁。2005—2010年的總趨勢與狀態(tài)指數(shù)達(dá)到了 0.85, 表明此時(shí)段相比于其他研究時(shí)段的流域土地利用變化呈現(xiàn)出明顯的動態(tài)不平衡狀態(tài), 各土地利用類型呈單向轉(zhuǎn)入或轉(zhuǎn)出的變化趨勢。

    圖2 泊江海子流域1990—2015年間的土地利用空間分布Fig.2 Spatial distribution of land use types in the Bojiang Lake Basin in 1990-2015

    2.3 土地利用變化驅(qū)動因素分析

    2.3.1 氣候變化驅(qū)動

    氣候變化對土地利用變化的作用較緩慢且相對穩(wěn)定, 具有累積效應(yīng)[28]。氣溫和降雨量的變化是氣候變化的主要表現(xiàn)[29]。研究區(qū)1990—2015年間氣溫和降雨量變化趨勢見圖 4。研究區(qū)年均氣溫整體呈升高趨勢, 1990—1996年呈降低趨勢, 1997年后升高。26年間的多年平均氣溫為6.6 ℃, 1990—1996年和1997—2015年的多年平均氣溫分別為 5.9 ℃和6.8 ℃, 后一階段年平均氣溫升高了 0.9 ℃。研究區(qū)年降雨量整體呈減少趨勢, 1998年之后減少趨勢更為明顯。26年間的多年平均降雨量為333.4 mm, 1990—1998年和1999—2015年的多年平均降雨量分別為355.7 mm和321.6 mm, 后一階段減少9.59%。

    圖3 泊江海子流域1990—2015年間不同時(shí)段的土地利用空間變化Fig.3 Spatial variations of land use types changes at different periods in the Bojiang Lake Basin during 1990-2015

    表2 泊江海子流域1990—2015年不同土地利用類型的面積占比變化Table 2 Variations in area percentages of different land use types in the Bojiang Lake Basin during 1990-2015 %

    表3 1990—2015年不同時(shí)段泊江海子流域各土地利用類型面積及指數(shù)變化Table 3 Variations in areas and indexes of land use types in the Bojiang Lake Basin in different periods from 1990 to 2015

    表4 1990—2015年不同時(shí)段泊江海子流域土地利用變化綜合趨勢與狀態(tài)指數(shù)及變化速度Table 4 Comprehensive indexes of trend and status and the rate of land use changes in the Bojiang Lake Basin in different periods from 1990 to 2015

    降雨、氣溫等氣候要素始終處于動態(tài)變化的過程中, 且共同對土地利用的變化產(chǎn)生影響。以研究期內(nèi)單個(gè)時(shí)段為例, 1990—1995年間的低溫氣候使地處高原的研究區(qū)內(nèi)耕地面積減少, 而因其引起的農(nóng)業(yè)用水量的下降, 也使得在降雨量未增加的條件下研究區(qū)內(nèi)水域面積增加。而受1996年后氣溫升高及1998年后降雨量減少的影響, 水域面積減少成為2000—2005年間研究區(qū)土地利用變化的突出特征之一。在整個(gè)研究期內(nèi), 研究區(qū)降雨量減少、氣溫升高及伴隨的蒸發(fā)增大, 是氣候因素對研究區(qū)土地利用變化產(chǎn)生驅(qū)動作用的主要表現(xiàn)形式。由于氣候暖干化使得流域范圍內(nèi)蒸散發(fā)增加, 流域徑流及入湖徑流量的減少, 使得水域面積顯著減少而灘地面積增加。1990—2015年間研究區(qū)高覆蓋度草地面積減少而中、低覆蓋度草地面積的增加, 也在一定程度上反映出流域土地利用變化受到氣候變化的驅(qū)動。

    圖4 泊江海子流域1990—2015年間年均氣溫(a)和降雨量(b)變化趨勢Fig.4 Change trends of annual mean temperature (a) and precipitation (b) in 1990-2015 in the Bojiang Lake Basin

    2.3.2 人類活動驅(qū)動

    人類活動是影響土地利用結(jié)構(gòu)變化的重要因子,人為因素如社會經(jīng)濟(jì)、制度政策、人類生產(chǎn)生活建設(shè)等對土地利用格局變化也起到了重要作用[30]。綜合本文對研究區(qū)土地利用變化分析以及已有的研究資料, 本文選取生態(tài)恢復(fù)工程建設(shè)和水利工程建設(shè)兩個(gè)指標(biāo), 對1990—2015年間流域土地利用變化受人類活動驅(qū)動的影響進(jìn)行分析。

    1)生態(tài)恢復(fù)工程建設(shè)

    研究區(qū)所在的鄂爾多斯市東勝區(qū)是天然林保護(hù)工程、退耕還林/還草工程等生態(tài)恢復(fù)工程的建設(shè)區(qū)。泊江海子鎮(zhèn)位于東勝區(qū)西部, 屬沙漠、沙化區(qū)治理的重點(diǎn)區(qū)域[31]。1998年建立的鄂爾多斯遺鷗國家級自然保護(hù)區(qū), 雖然沒有直接的生態(tài)恢復(fù)工程建設(shè), 但禁牧等措施在作用上與生態(tài)恢復(fù)工程建設(shè)是同向的。在生態(tài)恢復(fù)工程建設(shè)和禁牧等措施的實(shí)施之下,東勝區(qū)的林業(yè)生態(tài)狀況全面好轉(zhuǎn)[32]。2000—2003年東勝區(qū)荒漠化面積減少10.92%[31], 2004—2007年間東勝區(qū)完成天然林保護(hù)工程任務(wù) 40 km2, “三北”四期防護(hù)林工程人工造林約5.33 km2[32]。

    雖然筆者沒有查閱到研究區(qū)所在的泊江海子鎮(zhèn)區(qū)域內(nèi)的生態(tài)恢復(fù)工程造林面積, 且本文的土地利用變化分析也顯示研究區(qū)林草面積占流域的比例僅由1990年的73.15%增加到2015年的74.09%, 但有報(bào)道指出, 1999—2007年間泊江海子鎮(zhèn)區(qū)域內(nèi)的牧草高度增加了50~60 cm[33]。筆者先前也曾對本文研究區(qū)的植被覆蓋變化進(jìn)行過相關(guān)研究[34], 結(jié)果顯示研究區(qū) 1999—2013年間的年均歸一化植被指數(shù)(NDVI)相比1982— 1998年間增長了25.23%。在先前的研究中, 筆者對研究區(qū)2000年和2013年的250 m空間分辨率的植被指數(shù)數(shù)據(jù)進(jìn)行疊加, 得到了研究區(qū)14年間植被指數(shù)的空間變化(圖5)。2000—2013年間幾乎整個(gè)研究區(qū)內(nèi)的植被指數(shù)都呈增加趨勢, 增加值在0.1以上的面積達(dá)211.88 km2, 占研究區(qū)總面積的 28.46%, 且主要分布在河流兩側(cè)。而在1982—1998年間研究區(qū)平均植被指數(shù)值也僅為0.18左右, 由此可見在生態(tài)恢復(fù)工程建設(shè)和禁牧等措施的影響下, 研究區(qū)的植被得到了良好的生長。林草面積及植被覆蓋度的增加、未利用地面積的減少等土地利用變化與研究區(qū)內(nèi)生態(tài)恢復(fù)工程建設(shè)具有一定的緊密聯(lián)系, 是研究區(qū)土地利用變化受人類活動驅(qū)動的直觀表現(xiàn)。

    研究區(qū)生態(tài)恢復(fù)工程建設(shè)主要是以沙柳為主,與檸條、沙棘等樹種混交造林[31]。而沙柳屬于高耗水植被, 相比其他低耗水的固沙樹種如油蒿等, 相同面積的沙柳多消耗近 20%的水分[23]。圖 5顯示,2000—2013年間植被生長最好的區(qū)域主要分布在河流沿岸。這在一定程度上使得研究區(qū)在此時(shí)段內(nèi)的植被耗水增加。同時(shí), 草地面積約占流域總面積的70%, 研究區(qū)牧草高度的增加[33]也會導(dǎo)致流域內(nèi)地表徑流的減小[3,35]。而本文的土地利用變化分析也顯示, 研究區(qū)內(nèi)的水域面積在 2000—2005年間由13.61 km2減少到6.71 km2, 且在后續(xù)的10年內(nèi)一直呈減少的趨勢。由此分析可知, 生態(tài)恢復(fù)工程建設(shè)對流域內(nèi)水域面積的減少也產(chǎn)生了一定的驅(qū)動作用。

    圖5 2000—2013年間研究區(qū)植被指數(shù)空間變化Fig.5 Spatial variation of NDVI in 2000-2013 in the study area

    2)水利工程建設(shè)

    研究區(qū)屬西北干旱半干旱區(qū), 其所在的泊江海子鎮(zhèn)為鄂爾多斯市東勝區(qū)的重點(diǎn)農(nóng)業(yè)鎮(zhèn)。2000年前后, 開展了大量水利工程建設(shè)[36]。為了更好地分析研究區(qū)內(nèi)水利工程建設(shè)對土地利用變化的驅(qū)動, 筆者利用地面分辨率較高的 GoogleEarth對研究區(qū)內(nèi)的攔水壩和大口井兩種水利工程的數(shù)量進(jìn)行了粗略統(tǒng)計(jì), GoogleEarth影像日期為2009年11月27日,統(tǒng)計(jì)結(jié)果如圖6所示。

    圖6 2009年研究區(qū)內(nèi)的攔水壩和大口井空間分布Fig.6 Spatial distribution of dams and open wells in 2009 in the study area

    筆者共統(tǒng)計(jì)到攔水壩 13座, 截流面積約為60.83 km2, 占流域總面積的8.15%。同時(shí)統(tǒng)計(jì)到大口井28處, 寬度為10~30 m, 長度為50~100 m。攔水壩截水和大口井水方便了農(nóng)田的灌溉, 對研究區(qū)內(nèi)農(nóng)業(yè)發(fā)展具有重要的推動作用。同時(shí), 攔水壩使研究區(qū)內(nèi)地表水流量及湖泊補(bǔ)給量減少, 大口井也使研究內(nèi)水面無效蒸發(fā)增大。大口井的井水主要來自于工程開挖而滲出的地下潛水, 這在一定程度上勢必會引起研究區(qū)內(nèi)地下水位的下降, 進(jìn)而影響到研究區(qū)內(nèi)根系較淺的草地的生長。本文的土地利用變化分析也顯示, 耕地面積在研究期內(nèi)呈增加趨勢,其中2000—2010年間耕地面積增加了16.39 km2, 10年間年均增長率為15.71%。高覆蓋度草地和水域面積在研究期內(nèi)呈減少趨勢, 分別減少了流域總面積的7.59%和1.45%。

    綜合氣候變化和人類活動對研究區(qū)土地利用變化的驅(qū)動, 降雨減少、氣溫升高及伴隨的蒸發(fā)增大是氣候因素對研究區(qū)土地利用變化的主要驅(qū)動形式,生態(tài)恢復(fù)工程和水利工程建設(shè)是人類活動對研究區(qū)土地利用變化的主要驅(qū)動形式。研究區(qū)土地利用變化是氣候變化和人類活動綜合驅(qū)動的結(jié)果。

    3 結(jié)論與討論

    3.1 討論

    土地利用變化主要受氣候變化和人類活動兩大要素驅(qū)動, 已成為相關(guān)領(lǐng)域眾多研究者的共識。但在具體的驅(qū)動形式上, 因地形地貌、水文條件、流域面積等客觀因素以及人口數(shù)量、制度政策、經(jīng)濟(jì)發(fā)展等社會基礎(chǔ)條件上的差異, 氣候變化和人類活動在不同地區(qū)對土地利用變化驅(qū)動上的表現(xiàn)形式也會有差異。

    就本文研究區(qū)而言, 因四周高、中間低的地形特點(diǎn), 泊江海子湖泊不僅是研究區(qū)內(nèi)地表水的主要匯集區(qū), 同時(shí)也得到研究區(qū)內(nèi)淺層地下水的補(bǔ)給[22]。研究區(qū)地處西北內(nèi)陸高原, 降雨量少而蒸發(fā)強(qiáng)烈。研究區(qū)的這種水文氣象條件使得區(qū)內(nèi)各類陸地生態(tài)系統(tǒng)對氣候變化和人類活動干擾更為敏感。而在社會條件方面, 研究區(qū)所在的泊江海子鎮(zhèn)人口總量約為23 000人, 但其常住人口卻呈下降趨勢, 2015年的常住人口僅為 1 000人左右[37]。泊江海子鎮(zhèn)距東勝區(qū)約45 km, 其常住人口數(shù)量的下降, 除因生活水平提高而主動外遷的人口之外, 也與當(dāng)?shù)氐纳鷳B(tài)移民政策有關(guān)。2006—2010年間泊江海子鎮(zhèn)共轉(zhuǎn)移農(nóng)村人口 7 340人[37]。而在退耕還林、退牧還草的生態(tài)建設(shè)工程政策實(shí)施的影響之下, 泊江海子鎮(zhèn)的牲畜存欄量由2005年的12萬頭減少到2010年的不足3萬頭[37]?;谘芯繀^(qū)的自然環(huán)境和社會環(huán)境特征,本文未采用人口數(shù)量、牲畜存欄量、國民生產(chǎn)總值等社會經(jīng)濟(jì)基礎(chǔ)指標(biāo), 而是選取了更加具有針對性和代表性的生態(tài)恢復(fù)工程建設(shè)和水利工程建設(shè)兩個(gè)主要指標(biāo)對研究區(qū)土地利用變化的驅(qū)動因素進(jìn)行分析, 以便更加清晰地呈現(xiàn)研究區(qū)人類活動與土地利用變化的關(guān)系。

    水資源不合理利用是引發(fā)研究區(qū)生態(tài)系統(tǒng)服務(wù)功能退化的最主要原因之一, 而因生態(tài)系統(tǒng)服務(wù)功能不斷退化導(dǎo)致泊江海子湖泊中的湖心島作為遺鷗繁殖棲息地功能的喪失, 使得如何合理分配水資源以改善區(qū)內(nèi)生態(tài)系統(tǒng)服務(wù)功能, 成為研究區(qū)當(dāng)前面臨的最緊迫問題。合理規(guī)劃研究區(qū)的土地利用, 如發(fā)展種植低耗水的固沙林種、減少河流上游匯水區(qū)的耕地等, 合理分配流域水資源, 提高水資源利用效率, 是改善研究區(qū)生態(tài)系統(tǒng)服務(wù)功能、保障研究區(qū)生態(tài)系統(tǒng)可持續(xù)健康發(fā)展的重要措施。

    3.2 結(jié)論

    泊江海子流域生態(tài)系統(tǒng)脆弱, 流域內(nèi)土地利用方式和格局的變化直接影響著流域生態(tài)系統(tǒng)服務(wù)功能的種類和強(qiáng)度。近十幾年來, 因流域生態(tài)系統(tǒng)服務(wù)功能退化, 遺鷗繁殖棲息地受到極大威脅。本文基于Landsat影像數(shù)據(jù), 對研究區(qū)1990—2015年間的 6期土地利用數(shù)據(jù)進(jìn)行分類提取, 利用趨勢與狀態(tài)指數(shù)及凈變化速度等方法分析了研究區(qū) 26年間土地利用變化特征, 并對其驅(qū)動因素進(jìn)行分析, 得出如下結(jié)論:

    1)泊江海子流域1990—2015年間土地利用變化顯著。在空間分布變化上, 研究區(qū)東部、西部及中北部區(qū)域呈耕地面積增加、高覆蓋度草地減少的趨勢, 南部呈林地和低覆蓋度草地增加、高覆蓋度草地和未利用地減少的趨勢。在面積變化上, 研究期內(nèi)草地整體面積約占流域總面積的 70%, 高覆蓋度草地面積減少, 中、低覆蓋度草地面積增加。耕地面積增加, 占流域總面積的比例由1990年的9.85%增加到2015年的15.90%。未利用地面積減少, 占流域總面積的比例由1990年的13.24%減少到2015年的 5.62%。研究期內(nèi)水域面積減少了流域總面積的1.45%, 灘地面積增加了流域總面積的 1.41%, 林地和建設(shè)用地面積均呈弱增加趨勢。

    2)在土地利用類型轉(zhuǎn)換上, 3種覆蓋度草地之間及草地與其他土地利用類型之間處于動態(tài)平衡狀態(tài)。草地是耕地、林地和建設(shè)用地面積增加的主要轉(zhuǎn)入來源, 同時(shí)也是未利用地的主要轉(zhuǎn)出去向。水域面積減少是灘地面積增加的主要原因。1995—2000年是研究期內(nèi)各土地利用類型間轉(zhuǎn)換最為頻繁的研究時(shí)段。2005—2010年間土地利用變化動態(tài)不平衡狀態(tài)最為突出, 研究區(qū)各土地利用類型間呈單向轉(zhuǎn)入或轉(zhuǎn)出的變化趨勢。

    3)研究區(qū)土地利用變化受氣候變化和人類活動兩個(gè)驅(qū)動因素綜合作用的表現(xiàn)明顯。降雨減少、氣溫升高及伴隨的蒸發(fā)增大是氣候變化驅(qū)動研究區(qū)土地利用變化的主要表現(xiàn)形式, 生態(tài)恢復(fù)工程和水利工程建設(shè)是人類活動驅(qū)動研究區(qū)土地利用變化的主要表現(xiàn)形式。耕地、建設(shè)用地、林地面積的增加及未利用地面積的減少主要受人類活動的驅(qū)動影響,水域、高覆蓋度草地面積的減少及灘地、中、低覆蓋度草地面積的增加是氣候變化和人類活動綜合驅(qū)動的結(jié)果。

    致謝: 感謝中國科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心在本文土地利用類型數(shù)據(jù)驗(yàn)證過程中提供的數(shù)據(jù)支撐;特別感謝中國科學(xué)院地理科學(xué)與資源研究所的徐新良研究員在土地利用類型遙感解譯及驗(yàn)證方面中提供的指導(dǎo)和幫助。

    References

    [1]白曉飛, 陳煥偉.不同土地利用結(jié)構(gòu)生態(tài)系統(tǒng)服務(wù)功能價(jià)值的變化研究——以內(nèi)蒙古自治區(qū)伊金霍洛旗為例[J].中國生態(tài)農(nóng)業(yè)學(xué)報(bào), 2004, 12(1): 180–182

    Bai X F, Chen H W.The changes of ecosystem services and their values in various constructions of land use — A case study from Yijinhuoluo County of Inner Mogolia[J].Chinese Journal of Eco-Agriculture, 2004, 12(1): 180–182

    [2]陳穎, 石培基, 潘竟虎.高原生態(tài)城土地利用變化對生態(tài)系統(tǒng)服務(wù)價(jià)值的影響研究——以甘肅省民樂縣為例[J].水土保持研究, 2012, 19(2): 154–159

    Chen Y, Shi P J, Pan J H.Analysis on the effect of land-use change on ecosystem service value in the plateau eco-city —A case study of Minle County[J].Research of Soil and Water Conservation, 2012, 19(2): 154–159

    [3]賀克雕, 段昌群, 楊世美, 等.土地利用/覆被變化的水文水資源響應(yīng)研究綜述[J].水資源研究, 2015, 4(3):240–248

    He K D, Duan C Q, Yang S M, et al.Advances in hydrological and water resources response to land use/cover change[J].Journal of Water Resources Research, 2015, 4(3): 240–248

    [4]鄧慧平.氣候與土地利用變化對水文水資源的影響研究[J].地球科學(xué)進(jìn)展, 2001, 16(3): 436–441

    Deng H P.Impacts of climate change and land use-land cover change on hydrology and water resources[J].Advance in Earth Sciences, 2001, 16(3): 436–441

    [5]Meyer W B, Turner B L.Human population growth and global land-use/cover change[J].Annual Review of Ecology and Systematics, 1992, 23: 39–61

    [6]IPCC.Climate Change 2013: The Physical Science Basis.Contribution of Working GroupⅠto the Fifth Assessment Report of the Intergovernmental Panel on Climate Change[R].Cambridge, United Kingdom and New York, NY, USA:Cambridge University Press, 2013: 1535

    [7]Zhang Y, Wang T W, Cai C F, et al.Landscape pattern and transition under natural and anthropogenic disturbance in an arid region of northwestern China[J].International Journal of Applied Earth Observation and Geoinformation, 2016, 44:1–10

    [8]Dale V H.The relationship between land-use change and climate change[J].Ecological Applications, 1997, 7(3):753–769

    [9]Dilly O, Hüttl R F.Top-down and Europe-wide versus bottom-up and intra-regional identification of key issues for sustainability impact assessment[J].Environmental Science &Policy, 2009, 12(8): 1168–1176

    [10]李艷忠, 劉昌明, 劉小莽, 等.植被恢復(fù)工程對黃河中游土地利用/覆被變化的影響[J].自然資源學(xué)報(bào), 2016, 31(12):2005–2020

    Li Y Z, Liu C M, Liu X M, et al.Impact of the grain for green project on the land use/cover change in the middle Yellow River[J].Journal of Natural Resources, 2016, 31(12):2005–2020

    [11]張曉萍, 張櫓, 穆興民, 等.黃河中游河口-龍門區(qū)間多年平均流域水平衡特征——區(qū)域蒸散量估算模型驗(yàn)證與下墊面參數(shù)校核[J].地理學(xué)報(bào), 2007, 62(7): 753–763

    Zhang X P, Zhang L, Mu X M, et al.The mean annual water balance in the hekou-longmen section of the middle Yellow River: Testing of the regional scale water balance model and its calibration[J].Acta Geographica Sinica, 2007, 62(7):753–763

    [12]許炯心.黃河中游徑流可再生性對于人類活動和氣候變化的響應(yīng)[J].自然資源學(xué)報(bào), 2015, 30(3): 423–432

    Xu J X.Runoff renewability in the middle Yellow River in response to human activity and climate change[J].Journal of Natural Resources, 2015, 30(3): 423–432

    [13]羅婭, 楊勝天, 劉曉燕, 等.黃河河口鎮(zhèn)—潼關(guān)區(qū)間1998–2010年土地利用變化特征[J].地理學(xué)報(bào), 2014, 69(1):42–53

    Luo Y, Yang S T, Liu X Y, et al.Land use change in the reach from Hekouzhen to Tongguan of the Yellow River during 1998–2010[J].Acta Geographica Sinica, 2014, 69(1): 42–53

    [14]劉曉林, 楊勝天, 周旭, 等.1980年以來黃河內(nèi)蒙古段十大孔兌流域土地利用變化時(shí)空特征[J].南水北調(diào)與水利科技,2016, 14(1): 30–36

    Liu X L, Yang S T, Zhou X, et al.Spatial-temporal characteristics of land use changes in the ten tributaries of Yellow River in Inner-Mongolia since 1980[J].South-to-North Water Transfers and Water Science & Techology, 2016, 14(1):30–36

    [15]何芬奇, 任永奇, 郭玉民.內(nèi)蒙古桃-阿海子的生境演替與水鳥群落的興衰[J].濕地科學(xué)與管理, 2015, 11(2): 54–58

    He F Q, Ren Y Q, Guo Y M.Habitat succession and development of waterbird population in Tiaolimiao-Alashan Lake in Inner Mongolia[J].Wetland Science & Management, 2015,11(2): 54–58

    [16]魯瑞潔, 夏虹, 強(qiáng)明瑞, 等.近 130a來毛烏素沙漠北部泊江海子湖泊沉積記錄的氣候環(huán)境變化[J].中國沙漠, 2008,28(1): 44–49

    Lu R J, Xia H, Qiang M R, et al.Environmental evolution in past 130 years recorded by Bojianghaizi Lake sediments in North Mu Us desert, China[J].Journal of Desert Research,2008, 28(1): 44–49

    [17]楊志榮, 張梅青.鄂爾多斯泊江海子地區(qū) 800余年來的氣候與環(huán)境變化[J].湖南師范大學(xué)自然科學(xué)學(xué)報(bào), 1997, 20(4):74–81

    Yang Z R, Zhang M Q.Climatic and environmental changes since 800 a B P in Pojianghaizi Lake Area, Ordos Plateau[J].Natural Science of Hunan Normal University, 1997, 20(4):74–81

    [18]翟秋敏, 邱維理, 李容全, 等.內(nèi)蒙古安固里淖-泊江海子全新世中晚期湖泊沉積及其氣候意義[J].古地理學(xué)報(bào),2000, 2(2): 84–91

    Zhai Q M, Qiu W L, Li R Q, et al.The middle and late holocene lacustrine sediments and its climate significance of Angulinao-Bojianghaizi Lakes, Inner Mongolia[J].Journal of Palaeogeography, 2000, 2(2): 84–91

    [19]何芬奇, 任永奇.鄂爾多斯遺鷗繁殖地的變遷[J].大自然,2006, (4): 16–17

    He F Q, Ren Y Q.The changes ofLarus relictusbreeding grounds in Ordos[J].China Nature, 2006, (4): 16–17

    [20]張蔭蓀, 何芬奇, 陳容伯, 等.遺鷗繁殖生境選擇及其繁殖地濕地鳥類群落研究[J].動物學(xué)研究, 1993, 14(2):128–135

    Zhang Y S, He F Q, Chen R B, et al.Breeding habitat selection of the relict gull and the wetland bird community around its breeding sites[J].Zoological Research, 1993, 14(2):128–135

    [21]Wang F, Lliang L L, Zhang Y S, et al.Eco-hydrological model and critical conditions of hydrology of the wetland of ErdosLarus RelictusNature Reserve[J].Acta Ecologica Sinica, 2009, 29(5): 307–313

    [22]梁康, 婁華君, 程傳周.鄂爾多斯泊江海子流域地下水流特征[J].資源科學(xué), 2011, 33(6): 1089–1098

    Liang K, Lou H J, Cheng C Z.Characteristics of groundwater flow in the OrdosLarus relictusreserve wetland[J].Resources Science, 2011, 33(6): 1089–1098

    [23]梁犁麗, 王芳.鄂爾多斯遺鷗保護(hù)區(qū)植被-水資源模擬及其調(diào)控[J].生態(tài)學(xué)報(bào), 2010, 30(1): 109–119

    Liang L L, Wang F.Simulation of water resources evolution driven by vegetation construction and control of eco-hydrological processes in ErdosLarus relictusnature reserve[J].Acta Ecologica Sinica, 2010, 30(1): 109–119

    [24]邢小軍, 于向芝, 白兆勇, 等.鄂爾多斯遺鷗自然保護(hù)區(qū)濕地水量平衡分析[J].干旱區(qū)資源與環(huán)境, 2009, 23(6):100–103

    Xing X J, Yu X Z, Bai Z Y, et al.Analysis of water balance of the wetland in ErdosLarus relictusnature reserve[J].Journal of Arid Land Resources and Environment, 2009, 23(6):100–103

    [25]Liang K.Quantifying streamflow variations in ungauged lake basins by integrating remote sensing and water balance modelling: A case study of the ErdosLarus relictusnational nature reserve, China[J].Remote Sensing, 2017, 9(6): 588

    [26]何芬奇, 張蔭蓀, 葉恩琦, 等.鄂爾多斯桃力廟-阿拉善灣海子濕地鳥類群落研究與濕地生境評估[J].生物多樣性,1996, 4(4): 187–193

    He F Q, Zhang Y S, Ye E Q, et al.Study on the wetland bird community and assessment on wetland habitat of the Taolimiao-Alashan Nur in Ordos of Inner Mongolia, China[J].Chi-Chinese Biodiversity, 1996, 4(4): 187–193

    [27]Liu D, Cao C X, Chen W, et al.Monitoring and predicting the degradation of a semi-arid wetland due to climate change and water abstraction in the OrdosLarus relictusnational nature reserve, China[J].Geomatics, Natural Hazards and Risk, doi:10.1080/19475705.2016.1220024

    [28]王生霞, 丁永建, 葉柏生, 等.基于氣候變化和人類活動影響的土地利用分析——以新疆阿克蘇河流域綠洲為例[J].冰川凍土, 2012, 34(4): 828–835

    Wang S X, Ding Y J, Ye B S, et al.A land use analysis based on the influences of climate change and human activities — A case study in the Aksu Rver Basin oasis[J].Journal of Glaciology and Geocryology, 2012, 34(4): 828–835

    [29]王云琦, 齊實(shí), 孫閣, 等.氣候與土地利用變化對流域水資源的影響——以美國北卡羅萊納州 Trent流域?yàn)槔齕J].水科學(xué)進(jìn)展, 2011, 22(1): 51–58

    Wang Y Q, Qi S, Sun G, et al.Impacts of climate and land-use change on water resources in a watershed: A case study on the Trent River basin in North Carolina, USA[J].Advances in Water Science, 2011, 22(1): 51–58

    [30]李富程, 王青, 李國蓉.綿陽官司河流域土地利用/覆被變化特征及驅(qū)動力分析[J].農(nóng)業(yè)現(xiàn)代化研究, 2009, 30(6):727–730

    Li F C, Wang Q, Li G R.Analysis on land use/cover change characteristics and driving forces in Guansi River Valley of Mianyang City[J].Research of Agricultural Modernization,2009, 30(6): 727–730

    [31]吳青奮.打造綠色東勝——東勝區(qū)實(shí)施生態(tài)工程紀(jì)實(shí)[N].鄂爾多斯日報(bào), 2006-08-16(A02)

    Wu Q F.Building a green Dongsheng — The documentary of ecological engineering implementation in Dongsheng District[N].Ordos Daily, 2006-08-16(A02)

    [32]連桂梅, 張國棟.東勝區(qū)2004年以來完成重點(diǎn)工程造林16萬畝[N].鄂爾多斯日報(bào), 2008-01-09(002)

    Lian G M, Zhang G D.More than 100 km2forestry have been finished since 2004 in Dongsheng District[N].Ordos Daily,2008-01-09(002)

    [33]白鳳忠, 連桂梅.東勝區(qū)禁牧舍飼成效顯著[N].鄂爾多斯日報(bào), 2007-10-31(002)

    Bai F Z, Lian G M.Grazing prohibition have been strikingly successful in Dongsheng District[N].Ordos Daily, 2007-10-31(002)

    [34]Yan G Z, Lou H J, Liang K, et al.Dynamics and driving forces of Bojiang Lake area in ErdosLarus relictusNational Nature Reserve, China[J].Quaternary International, 2017, doi:10.1016/j.quaint.2016.12.034

    [35]甘紅, 劉彥隨, 李憲文.區(qū)域土地利用變化與水資源利用相關(guān)分析[J].南京師大學(xué)報(bào): 自然科學(xué)版, 2003, 26(3):82–88

    Gan H, Liu Y S, Li X W.Analysis of the correlativity between regional land use change and water resources optimal use[J].Journal of Nanjing Normal University: Natural Science, 2003, 26(3): 82–88

    [36]楊紅巖.2002年東勝區(qū)政府工作報(bào)告[EB/OL].東勝之窗,[2017-03-05].http://www.ds.gov.cn/zw/gzbg103534/index.html

    Yang H Y.Dongsheng District government work report in 2002[EB/OL].The Window of Dongsheng District, [2017-03-05].http://www.ds.gov.cn/zw/gzbg103534/index.html

    [37]東勝區(qū)政府.東勝區(qū)區(qū)情介紹[EB/OL].東勝之窗,[2017-03-05].http://www.ds.gov.cn/qq/gk/qqjs103236/

    Dongsheng District Government.Regional orientation of Dongsheng District[EB/OL].The Window of Dongsheng District,[2017-03-05].http://www.ds.gov.cn/qq/gk/qqjs 103236/

    Characteristics and driving factors of land use change in the Bojiang Lake Basin in Ordos Plateau, China*

    YAN Guozhen1,2, ZHANG Zheng1, LIANG Kang2**, LOU Huajun2
    (1.College of Forestry, Beijing Forestry University, Beijing 100083, China; 2.Key Laboratory of Water Cycle and Related Land
    Surface Processes, Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing 100101, China)

    Land use change can directly affect the variety and intensity of terrestrial ecosystem service functions.The Bojiang Lake Basin (BLB) is located in the Ordos Plateau in China and it belongs to the fragile ecosystem of the arid/semi-arid region.In the last decade, degradation of ecosystem services function in BLB has further jeopardized the endangeredLarus relictusbirds.Land use structure and pattern in BLB have changed significantly in recent decades.However, few studies have been done on land use change and the related driving factors.This study aimed to investigate the characteristics and the related driving factors of land use change in BLB during 1990–2015 to provide a reference base for the improvement of ecosystem services function and the protection ofL.relictusbirds.In the study, six periods (1990, 1995, 2000, 2005, 2010 and 2015) of 30 m resolution of land use data were extracted from Landsat remote sensing images.Land use change was investigated using trend and status index methods and then the main driving factors (i.e., climate change and human activity) analyzed.The results showed that grassland, farmland and unused land were the main land use types in 1990–2015, and there was significant change in land use in the BLB study area.Based on spatial variation, there was an increasing trend in farmland in the eastern,western and north-central regions, and then a decreasing trend in unused land in the southern region.The characteristics of the variation in area and dynamic conversion of each land use type were different.There was a dynamic balance in transformation in high, medium and low coverage grasslands, occupying some 70% of the lake basin.Farmland and unused land respectively had an apparent increase and decrease in trend.There was 6.05% increase in farmland area and 7.62% decrease in unused land area in the lake basin.Grassland was the main source of income for farmlands and therefore the main export for unused lands.The water area was decreased by 1.45% of the total area in 1990–2015, with bottomland as the main export.In the period 1995–2000 there was frequently transformation in land use, which in 2005–2010 had the most dynamic imbalance in land use change.Both climatic factors and human activities significantly affected land use change in the BLB study area.Specifically in 1990–2015, the transformation in land use types was due to the combined effects of warming and drying climate (expressed as declining precipitation and rising temperature) and intensifying anthropogenic disturbances (expressed as ecological restoration and water conservancy projects).The increase in farmland, woodland and buildup area and the decrease in unused land were more affected by anthropogenic disturbances.The decrease in water body and high coverage grassland and the increase in medium and low coverage grasslands were due to both climatic factors and human activities.Therefore, a reasonable land use planning such as restricting planting of forest with high water use, reducing the area of farmland in the catchment of the upper reaches of the rivers was critical for the protection for healthy and sustainable development of the study area.

    Jan.4, 2017; accepted May 18, 2017

    Bojiang Lake Basin; Land use change; Trend and status index; Driving factor; Ecosystem services function; ErdosLarus relictusNational Nature Reserve

    F301

    A

    1671-3990(2017)11-1693-14

    10.13930/j.cnki.cjea.170012

    閆國振, 張征, 梁康, 婁華君.鄂爾多斯高原泊江海子流域土地利用動態(tài)變化及驅(qū)動因素分析[J].中國生態(tài)農(nóng)業(yè)學(xué)報(bào),2017, 25(11): 1693-1706

    Yan G Z, Zhang Z, Liang K, Lou H J.Characteristics and driving factors of land use change in the Bojiang Lake Basin in Ordos Plateau, China[J].Chinese Journal of Eco-Agriculture, 2017, 25(11): 1693-1706

    * 國家自然科學(xué)基金項(xiàng)目(41501032)資助

    ** 通訊作者: 梁康, 主要從事流域水文過程變化機(jī)理與水資源安全研究。E-mail: liangk@igsnrr.ac.cn

    閆國振, 主要研究方向?yàn)樯鷳B(tài)學(xué)。E-mail: ygzh4161@163.com

    2017-01-04 接受日期: 2017-05-18

    * This study was supported by the National Natural Science Foundation of China (41501032).

    ** Corresponding author, E-mail: liangk@igsnrr.ac.cn

    猜你喜歡
    江海覆蓋度土地利用
    呼和浩特市和林格爾縣植被覆蓋度變化遙感監(jiān)測
    基于NDVI的晉州市植被覆蓋信息提取
    在露水與江海面前,你會如何選擇?
    低覆蓋度CO分子在Ni(110)面的吸附研究
    江海
    黃河之聲(2018年4期)2018-04-24 02:43:58
    土地利用生態(tài)系統(tǒng)服務(wù)研究進(jìn)展及啟示
    濱??h土地利用挖潛方向在哪里
    6000m~3江海直達(dá)LNG運(yùn)輸船推進(jìn)系統(tǒng)分析
    船海工程(2015年4期)2016-01-05 15:53:32
    青海湟水流域植被覆蓋度時(shí)空變化分析
    論低碳經(jīng)濟(jì)與轉(zhuǎn)變土地利用方式
    国产精品一区二区三区四区久久 | 国产成人免费无遮挡视频| 国产成人免费无遮挡视频| 91九色精品人成在线观看| 亚洲av五月六月丁香网| av在线播放免费不卡| 久久精品人人爽人人爽视色| 欧美乱妇无乱码| 熟女少妇亚洲综合色aaa.| 亚洲av成人一区二区三| 亚洲专区中文字幕在线| 久久精品影院6| 18禁裸乳无遮挡免费网站照片 | 亚洲情色 制服丝袜| 欧美日韩精品网址| 国产免费av片在线观看野外av| 国产伦人伦偷精品视频| 桃色一区二区三区在线观看| 国产激情欧美一区二区| 老司机靠b影院| 国产一卡二卡三卡精品| 十八禁人妻一区二区| 性欧美人与动物交配| 精品高清国产在线一区| 人人妻,人人澡人人爽秒播| 黄色视频不卡| 欧美乱妇无乱码| 最近最新中文字幕大全电影3 | 丝袜在线中文字幕| 黄色视频不卡| 别揉我奶头~嗯~啊~动态视频| 成人18禁高潮啪啪吃奶动态图| 午夜激情av网站| 在线观看一区二区三区| 一级,二级,三级黄色视频| 精品欧美一区二区三区在线| 亚洲一区中文字幕在线| 黄色女人牲交| 十八禁人妻一区二区| 亚洲欧洲精品一区二区精品久久久| 久久久久国产精品人妻aⅴ院| 午夜影院日韩av| 嫩草影院精品99| bbb黄色大片| 黄色毛片三级朝国网站| 精品人妻1区二区| 亚洲av日韩精品久久久久久密| 中文字幕人妻熟女乱码| 亚洲中文日韩欧美视频| 涩涩av久久男人的天堂| 日本a在线网址| 女人被狂操c到高潮| 国产成+人综合+亚洲专区| 99国产极品粉嫩在线观看| or卡值多少钱| 狂野欧美激情性xxxx| 女人爽到高潮嗷嗷叫在线视频| 神马国产精品三级电影在线观看 | 亚洲激情在线av| 一级,二级,三级黄色视频| svipshipincom国产片| 亚洲国产精品合色在线| 国产91精品成人一区二区三区| 亚洲精华国产精华精| 婷婷丁香在线五月| 变态另类成人亚洲欧美熟女 | 久久青草综合色| 身体一侧抽搐| 女生性感内裤真人,穿戴方法视频| 一级作爱视频免费观看| 99久久99久久久精品蜜桃| 大型黄色视频在线免费观看| 性色av乱码一区二区三区2| 少妇 在线观看| 搡老岳熟女国产| 性色av乱码一区二区三区2| 大码成人一级视频| 少妇 在线观看| 1024香蕉在线观看| 亚洲人成电影观看| 亚洲情色 制服丝袜| 999精品在线视频| 变态另类丝袜制服| 最新在线观看一区二区三区| 不卡av一区二区三区| 中文字幕人成人乱码亚洲影| 日韩高清综合在线| 国产激情欧美一区二区| 免费看a级黄色片| 叶爱在线成人免费视频播放| 国产单亲对白刺激| 亚洲精品一卡2卡三卡4卡5卡| 久久天堂一区二区三区四区| 一级作爱视频免费观看| 国产免费男女视频| 他把我摸到了高潮在线观看| 欧美成狂野欧美在线观看| 91老司机精品| 此物有八面人人有两片| 亚洲性夜色夜夜综合| 国产亚洲欧美98| 搞女人的毛片| 亚洲激情在线av| 亚洲色图 男人天堂 中文字幕| 中文字幕久久专区| 国产亚洲精品第一综合不卡| 中亚洲国语对白在线视频| 国产高清videossex| 亚洲熟女毛片儿| 国产精品二区激情视频| 欧美色欧美亚洲另类二区 | 久久久久九九精品影院| 国产激情欧美一区二区| 天天躁狠狠躁夜夜躁狠狠躁| 日韩高清综合在线| 午夜日韩欧美国产| av欧美777| 高清毛片免费观看视频网站| 一级毛片高清免费大全| 精品国产一区二区三区四区第35| 禁无遮挡网站| 波多野结衣巨乳人妻| 乱人伦中国视频| 极品教师在线免费播放| 国产精品永久免费网站| 国产欧美日韩一区二区三区在线| 欧美日韩亚洲综合一区二区三区_| 亚洲成国产人片在线观看| 免费无遮挡裸体视频| 久久中文看片网| 好男人在线观看高清免费视频 | 老司机靠b影院| 午夜福利免费观看在线| √禁漫天堂资源中文www| 国产视频一区二区在线看| 中文字幕人妻丝袜一区二区| 久久久久国产一级毛片高清牌| 成年版毛片免费区| 国产精品九九99| 两人在一起打扑克的视频| 老司机深夜福利视频在线观看| 久久狼人影院| 欧美成狂野欧美在线观看| 男男h啪啪无遮挡| 亚洲欧美日韩无卡精品| 91国产中文字幕| 免费一级毛片在线播放高清视频 | 国产av一区在线观看免费| 一区二区三区激情视频| 久久久国产精品麻豆| 亚洲欧美精品综合久久99| 亚洲专区国产一区二区| 国产成人精品在线电影| 久久午夜综合久久蜜桃| 久久国产乱子伦精品免费另类| a级毛片在线看网站| 亚洲情色 制服丝袜| 欧美+亚洲+日韩+国产| 久久中文字幕一级| 亚洲欧洲精品一区二区精品久久久| 悠悠久久av| 国产av精品麻豆| svipshipincom国产片| 老司机靠b影院| 国内毛片毛片毛片毛片毛片| 久久国产精品人妻蜜桃| 大型黄色视频在线免费观看| 久久精品91无色码中文字幕| 真人做人爱边吃奶动态| 在线观看一区二区三区| bbb黄色大片| 日韩精品中文字幕看吧| 国产精品九九99| 中出人妻视频一区二区| 日本免费一区二区三区高清不卡 | 亚洲国产高清在线一区二区三 | 亚洲精华国产精华精| 国产精品电影一区二区三区| 欧美在线黄色| 非洲黑人性xxxx精品又粗又长| 欧美国产精品va在线观看不卡| 亚洲视频免费观看视频| 无人区码免费观看不卡| 97人妻天天添夜夜摸| 日韩三级视频一区二区三区| 精品久久久久久成人av| 嫩草影视91久久| 啪啪无遮挡十八禁网站| 热re99久久国产66热| 亚洲片人在线观看| 精品一区二区三区四区五区乱码| 国产精品98久久久久久宅男小说| 国产午夜精品久久久久久| 欧美中文日本在线观看视频| 亚洲成人久久性| 国产精品av久久久久免费| 成人18禁在线播放| 婷婷精品国产亚洲av在线| 夜夜爽天天搞| 亚洲国产日韩欧美精品在线观看 | 亚洲人成电影观看| 色婷婷久久久亚洲欧美| 久久久国产精品麻豆| 午夜成年电影在线免费观看| 人人妻,人人澡人人爽秒播| 久久精品91无色码中文字幕| 久久久久久久久中文| 在线观看免费视频网站a站| 好男人在线观看高清免费视频 | 亚洲国产看品久久| 在线观看日韩欧美| 欧美精品啪啪一区二区三区| 国产一区二区三区在线臀色熟女| 法律面前人人平等表现在哪些方面| 纯流量卡能插随身wifi吗| 亚洲av第一区精品v没综合| 18禁美女被吸乳视频| 黄色丝袜av网址大全| 久久久久亚洲av毛片大全| 日本 欧美在线| 免费一级毛片在线播放高清视频 | 高潮久久久久久久久久久不卡| 女人被躁到高潮嗷嗷叫费观| 亚洲国产精品成人综合色| 啦啦啦观看免费观看视频高清 | 成人18禁高潮啪啪吃奶动态图| 国产精品久久久久久亚洲av鲁大| 一级片免费观看大全| 国产成人精品久久二区二区免费| 中文字幕人妻丝袜一区二区| 伦理电影免费视频| 免费少妇av软件| 国产av又大| 99国产精品一区二区三区| 久久婷婷人人爽人人干人人爱 | 亚洲第一av免费看| 老司机午夜十八禁免费视频| 女人被躁到高潮嗷嗷叫费观| 大香蕉久久成人网| 9热在线视频观看99| 国产一区二区激情短视频| 国产激情久久老熟女| 日韩免费av在线播放| 人成视频在线观看免费观看| 又黄又粗又硬又大视频| 国产黄a三级三级三级人| 成人亚洲精品av一区二区| 伊人久久大香线蕉亚洲五| 电影成人av| 91九色精品人成在线观看| 欧美黑人欧美精品刺激| 久久伊人香网站| 久久欧美精品欧美久久欧美| 18美女黄网站色大片免费观看| 在线观看舔阴道视频| 欧美老熟妇乱子伦牲交| 亚洲avbb在线观看| 精品久久久久久,| 香蕉国产在线看| 成年版毛片免费区| 亚洲精品中文字幕一二三四区| 日日爽夜夜爽网站| 亚洲男人天堂网一区| 亚洲欧美精品综合久久99| 自线自在国产av| 嫁个100分男人电影在线观看| 成人永久免费在线观看视频| 国产一区二区三区在线臀色熟女| 一级毛片精品| 99精品欧美一区二区三区四区| 中出人妻视频一区二区| 国产亚洲精品第一综合不卡| 精品一区二区三区四区五区乱码| 日韩免费av在线播放| 十八禁人妻一区二区| 不卡av一区二区三区| 亚洲av日韩精品久久久久久密| 最近最新中文字幕大全电影3 | 十八禁网站免费在线| 成人亚洲精品av一区二区| 亚洲伊人色综图| 青草久久国产| 日韩中文字幕欧美一区二区| www.www免费av| 欧美精品啪啪一区二区三区| 香蕉丝袜av| 免费一级毛片在线播放高清视频 | 99久久99久久久精品蜜桃| 国产精品二区激情视频| xxx96com| 国产欧美日韩一区二区三| 首页视频小说图片口味搜索| 波多野结衣高清无吗| 国产av又大| 一级黄色大片毛片| 亚洲黑人精品在线| 手机成人av网站| 黑丝袜美女国产一区| 淫妇啪啪啪对白视频| 国产蜜桃级精品一区二区三区| 最新在线观看一区二区三区| 免费看a级黄色片| 免费在线观看黄色视频的| 久久精品91蜜桃| 免费人成视频x8x8入口观看| 免费在线观看视频国产中文字幕亚洲| 日韩视频一区二区在线观看| 真人做人爱边吃奶动态| 久久久久久国产a免费观看| 9色porny在线观看| 宅男免费午夜| 一本久久中文字幕| 91大片在线观看| 日韩欧美一区视频在线观看| 精品久久蜜臀av无| 免费少妇av软件| 日韩一卡2卡3卡4卡2021年| 老汉色∧v一级毛片| 十八禁人妻一区二区| 97人妻天天添夜夜摸| 亚洲精品中文字幕一二三四区| 久久精品国产综合久久久| 国产一区二区在线av高清观看| 久久人人爽av亚洲精品天堂| 精品国产美女av久久久久小说| 欧美一级毛片孕妇| 97人妻精品一区二区三区麻豆 | 老司机午夜福利在线观看视频| 人人妻,人人澡人人爽秒播| 操出白浆在线播放| 嫩草影视91久久| 精品熟女少妇八av免费久了| 免费在线观看影片大全网站| 两性夫妻黄色片| 亚洲欧美激情在线| 日韩大码丰满熟妇| 亚洲人成网站在线播放欧美日韩| 久久精品91无色码中文字幕| 国产精品1区2区在线观看.| 久久久久久久精品吃奶| 国产乱人伦免费视频| 51午夜福利影视在线观看| 国产伦一二天堂av在线观看| 男女之事视频高清在线观看| 欧美绝顶高潮抽搐喷水| 日本免费a在线| 精品卡一卡二卡四卡免费| 久久久久久久午夜电影| 亚洲av片天天在线观看| 叶爱在线成人免费视频播放| 欧美绝顶高潮抽搐喷水| 亚洲一码二码三码区别大吗| 757午夜福利合集在线观看| 男女之事视频高清在线观看| 色老头精品视频在线观看| 中文字幕另类日韩欧美亚洲嫩草| 成人永久免费在线观看视频| 亚洲欧美一区二区三区黑人| 女性被躁到高潮视频| 51午夜福利影视在线观看| 久久精品亚洲精品国产色婷小说| 亚洲精品中文字幕一二三四区| 国产免费男女视频| 69av精品久久久久久| 天天一区二区日本电影三级 | 在线观看免费视频日本深夜| av超薄肉色丝袜交足视频| 精品福利观看| 1024视频免费在线观看| 亚洲五月色婷婷综合| 亚洲精品久久国产高清桃花| 亚洲五月色婷婷综合| 亚洲精品久久国产高清桃花| 久久久久久免费高清国产稀缺| 日韩大码丰满熟妇| 身体一侧抽搐| 国产成人av教育| 日韩欧美国产一区二区入口| 在线av久久热| 一级a爱片免费观看的视频| 久久久久国产一级毛片高清牌| 一级a爱片免费观看的视频| 欧美一级毛片孕妇| 午夜福利一区二区在线看| 色尼玛亚洲综合影院| 色综合亚洲欧美另类图片| 亚洲天堂国产精品一区在线| 午夜福利,免费看| 精品第一国产精品| 亚洲精品粉嫩美女一区| 91在线观看av| 欧美成人性av电影在线观看| 成年版毛片免费区| 99re在线观看精品视频| 国产成人啪精品午夜网站| 久久婷婷人人爽人人干人人爱 | 欧美日韩中文字幕国产精品一区二区三区 | 最近最新中文字幕大全电影3 | 叶爱在线成人免费视频播放| 757午夜福利合集在线观看| 欧美在线一区亚洲| 亚洲成人免费电影在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲全国av大片| 国内精品久久久久久久电影| 午夜免费成人在线视频| 精品午夜福利视频在线观看一区| 亚洲精品久久成人aⅴ小说| 成人18禁在线播放| 亚洲久久久国产精品| 一级毛片女人18水好多| 亚洲av电影在线进入| 久久精品国产综合久久久| 日本免费a在线| 日韩精品中文字幕看吧| a在线观看视频网站| 成人特级黄色片久久久久久久| a级毛片在线看网站| 一二三四在线观看免费中文在| 麻豆av在线久日| 动漫黄色视频在线观看| 无人区码免费观看不卡| 亚洲国产精品999在线| 色综合欧美亚洲国产小说| avwww免费| 午夜久久久久精精品| 天堂动漫精品| 一区二区三区国产精品乱码| 欧美日韩亚洲综合一区二区三区_| 在线十欧美十亚洲十日本专区| cao死你这个sao货| 国产精品亚洲一级av第二区| 波多野结衣一区麻豆| 韩国精品一区二区三区| 久久久国产成人精品二区| av欧美777| 久久影院123| 国产精品二区激情视频| 国产不卡一卡二| 天堂影院成人在线观看| 黄片小视频在线播放| 久久伊人香网站| 人人妻人人爽人人添夜夜欢视频| 热re99久久国产66热| 国产伦人伦偷精品视频| 黑人巨大精品欧美一区二区mp4| 后天国语完整版免费观看| 成人国产一区最新在线观看| 精品午夜福利视频在线观看一区| 两个人免费观看高清视频| 黄片播放在线免费| 成在线人永久免费视频| 老司机在亚洲福利影院| 欧美色视频一区免费| 97碰自拍视频| 美女扒开内裤让男人捅视频| 男女床上黄色一级片免费看| 精品高清国产在线一区| 大型av网站在线播放| 亚洲熟女毛片儿| 免费在线观看黄色视频的| 99国产极品粉嫩在线观看| 婷婷丁香在线五月| 美女高潮喷水抽搐中文字幕| 亚洲午夜理论影院| 久久亚洲真实| 少妇粗大呻吟视频| 成年女人毛片免费观看观看9| 国产又爽黄色视频| 中文字幕最新亚洲高清| 国产午夜精品久久久久久| 亚洲一区二区三区不卡视频| 中国美女看黄片| 老熟妇仑乱视频hdxx| 免费不卡黄色视频| 色综合亚洲欧美另类图片| 啦啦啦韩国在线观看视频| 俄罗斯特黄特色一大片| 成人国产综合亚洲| 国产欧美日韩一区二区精品| 波多野结衣av一区二区av| 成年女人毛片免费观看观看9| 精品国产美女av久久久久小说| 成人欧美大片| av视频在线观看入口| 亚洲国产高清在线一区二区三 | 法律面前人人平等表现在哪些方面| 国产成年人精品一区二区| 99久久99久久久精品蜜桃| cao死你这个sao货| 免费在线观看影片大全网站| 国产97色在线日韩免费| 国产91精品成人一区二区三区| 麻豆一二三区av精品| 国产亚洲精品久久久久5区| 国产精品一区二区精品视频观看| 亚洲色图 男人天堂 中文字幕| 最好的美女福利视频网| 亚洲欧洲精品一区二区精品久久久| 国产三级黄色录像| 搞女人的毛片| 两性夫妻黄色片| 涩涩av久久男人的天堂| 欧美大码av| 亚洲人成77777在线视频| 成年版毛片免费区| 亚洲人成电影免费在线| 男人舔女人下体高潮全视频| 欧美一级a爱片免费观看看 | 琪琪午夜伦伦电影理论片6080| 在线十欧美十亚洲十日本专区| 青草久久国产| 亚洲国产精品sss在线观看| 伦理电影免费视频| 免费高清在线观看日韩| 侵犯人妻中文字幕一二三四区| 如日韩欧美国产精品一区二区三区| 91成人精品电影| 多毛熟女@视频| 久久精品影院6| 亚洲va日本ⅴa欧美va伊人久久| 国产亚洲av嫩草精品影院| 欧美日本亚洲视频在线播放| 黄色丝袜av网址大全| 成人亚洲精品一区在线观看| 国产精品免费视频内射| 精品一区二区三区四区五区乱码| 丰满的人妻完整版| 久久久久久久久中文| 国产精品 国内视频| 男人舔女人的私密视频| 国产亚洲欧美精品永久| 99国产精品免费福利视频| 欧美 亚洲 国产 日韩一| 嫁个100分男人电影在线观看| 国产av在哪里看| 欧美大码av| 亚洲九九香蕉| 999久久久精品免费观看国产| 夜夜夜夜夜久久久久| 午夜福利欧美成人| 99精品在免费线老司机午夜| 欧美在线一区亚洲| 成人国产一区最新在线观看| 高清在线国产一区| 国产精品精品国产色婷婷| 国语自产精品视频在线第100页| 制服诱惑二区| 丁香欧美五月| 国产激情欧美一区二区| 母亲3免费完整高清在线观看| 欧美一区二区精品小视频在线| 精品久久蜜臀av无| 亚洲精品粉嫩美女一区| 欧美精品亚洲一区二区| 国内精品久久久久久久电影| 久久中文看片网| 搞女人的毛片| 午夜精品久久久久久毛片777| 亚洲国产毛片av蜜桃av| 免费在线观看日本一区| av网站免费在线观看视频| 自线自在国产av| 亚洲欧洲精品一区二区精品久久久| 国产精品精品国产色婷婷| 激情视频va一区二区三区| 欧美丝袜亚洲另类 | 精品高清国产在线一区| 这个男人来自地球电影免费观看| 亚洲美女黄片视频| 亚洲天堂国产精品一区在线| 国产亚洲精品久久久久5区| 麻豆久久精品国产亚洲av| 亚洲avbb在线观看| or卡值多少钱| 久久久久久免费高清国产稀缺| 精品午夜福利视频在线观看一区| av有码第一页| 老熟妇乱子伦视频在线观看| bbb黄色大片| 1024香蕉在线观看| 超碰成人久久| 久9热在线精品视频| 一区二区日韩欧美中文字幕| 在线国产一区二区在线| 午夜福利视频1000在线观看 | 91老司机精品| 亚洲国产精品成人综合色| 久久久久亚洲av毛片大全| 国产精品一区二区免费欧美| 亚洲av熟女| 天天躁夜夜躁狠狠躁躁| 91大片在线观看| 欧美成人性av电影在线观看| 久久人人97超碰香蕉20202| 制服人妻中文乱码| 两个人视频免费观看高清| 一进一出抽搐gif免费好疼| 黄色丝袜av网址大全| 欧美日韩乱码在线| 丝袜在线中文字幕| www.熟女人妻精品国产| 午夜亚洲福利在线播放| 18禁观看日本| 两个人看的免费小视频| 亚洲国产精品成人综合色| 免费看a级黄色片| 国产99久久九九免费精品| 视频在线观看一区二区三区| 美女大奶头视频| 99国产精品99久久久久| 国产精华一区二区三区| 欧美在线黄色| 国产亚洲精品av在线| 脱女人内裤的视频|