1973—2015年塔爾湖農(nóng)業(yè)區(qū)土地利用/覆蓋變化分析

葉夢(mèng)旎 蔡茂堂 傅連珍 楊星辰 李成路 周維 劉超

摘要以河套農(nóng)業(yè)區(qū)塔爾湖地區(qū)1973—2015年的5期遙感影像為主要數(shù)據(jù)源，應(yīng)用分層提取和監(jiān)督分類方法，獲取了研究區(qū)土地利用/覆蓋信息，并對(duì)土地利用/覆蓋動(dòng)態(tài)進(jìn)行了分析。結(jié)果表明，近42年來，塔爾湖地區(qū)土地利用/覆蓋格局沒有發(fā)生顯著變化。耕地面積除在1987—1995年有所降低外，基本保持著每年0.4%的速率增長(zhǎng)。城鄉(xiāng)工礦居民用地面積呈逐年增加的趨勢(shì)。未利用土地面積逐年減少，其中以沙地面積減少最多。土地鹽堿化總面積在1973—2006年穩(wěn)步降低，而在2006年后轉(zhuǎn)為增長(zhǎng)。塔爾湖地區(qū)土地利用/覆蓋變化受控于自然和人為因素的雙重作用。1973—2015年研究區(qū)在土壤鹽堿化改良和沙漠化防治等方面取得良好的效果，但2006—2015年不合理的灌溉和耕地的過度利用使得土壤鹽堿化程度增加。

關(guān)鍵詞衛(wèi)星遙感；土地利用/覆蓋變化；鹽堿化；塔爾湖

中圖分類號(hào)F301.2文獻(xiàn)標(biāo)識(shí)碼

A文章編號(hào)0517-6611（2017）32-0194-06

Land Use/Cover Changes in Tarhu Agricultural Region from 1973 to 2015

YE Mengni1，CAI Maotang2*，F(xiàn)U Lianzhen3 et al（1.Institute of Geology，Chinese Academy of Geological Sciences，Beijing 100037；2.Institute of Geomechanics，Chinese Academy of Geological Sciences，Beijing 100081； 3.Chinese Academy of Land and Resource Economics，Beijing 101149）

AbstractBased on visual interpretation of remote sensing images in 1973，1984，1995，2006 and 2015，the land use/cover databases of Tarhu region was established.Then the process，trends of the land use / cover change and its driving forces were analyzed.The results showed that：The basic pattern of land use/cover of Tarhu region didn't change from 1973 to 2015.The area of farmland and woodland increased in 0.4% annual growth rate，except from 1987 to 1995.The urbanland area increased more rapidly，while the unused land （mainly sandy land） area and water area decreased since 2006.The soil salinization land area steadily reduced from 1973 to 2006，and then changed to increase since 2006.The land use/cover change of Tarhu region was a result of combination of natural and human factors.The soil salinization land improvement and prevention and control of desertification over during 1973-2015 have made good effect，while unreasonable irrigation and the excessive use of cultivated land during 2006-2015 increased land salinization.

Key wordsRemote sensing；Land use/cover change；Soil salinization；Tarhu

土地是地球系統(tǒng)中的重要組成部分，也是農(nóng)業(yè)經(jīng)濟(jì)建設(shè)中不可分割的部分，隨著社會(huì)經(jīng)濟(jì)的發(fā)展，人們對(duì)于土地的需求量不斷增加。人類-環(huán)境交互作用改變了原來的土地覆蓋，直接影響著生物圈-大氣圈交互作用、生物地球化學(xué)循環(huán)以及資源環(huán)境的可持續(xù)利用[1-4]。因此，土地利用/覆蓋變化已成為當(dāng)前各國(guó)學(xué)術(shù)界的研究熱點(diǎn)[5-7]。對(duì)土地利用/覆蓋的變化過程及驅(qū)動(dòng)機(jī)制研究，可以揭示出土地利用/覆蓋對(duì)人類活動(dòng)和自然因素的響應(yīng)，并為土地利用/覆蓋變化及其相關(guān)領(lǐng)域的研究提供基礎(chǔ)信息。許多學(xué)者已利用遙感與地理信息系統(tǒng)對(duì)不同區(qū)域的土地利用變化進(jìn)行了大量研究[8-13]。

資源短缺與生態(tài)環(huán)境惡化的矛盾日益突出，合理利用鹽堿地成為改良耕地質(zhì)量與提升糧食產(chǎn)量的核心問題之一。塔爾湖地區(qū)位于內(nèi)蒙古河套農(nóng)業(yè)灌區(qū)，是我國(guó)主要產(chǎn)糧地區(qū)之一，也是我國(guó)典型的鹽堿化和沙化土壤地區(qū)。土壤鹽堿化和沙化不但造成了土地資源短缺和生態(tài)環(huán)境惡化，還對(duì)該地區(qū)的農(nóng)業(yè)發(fā)展造成巨大損失[14-15]。因此，開展塔爾湖地區(qū)土壤鹽堿化和沙化研究，了解塔爾湖農(nóng)業(yè)區(qū)土壤鹽堿化變化特征及其形成機(jī)制，對(duì)于合理開發(fā)利用土地資源，控制鹽堿化和沙化的發(fā)生、發(fā)展具有重要意義。

筆者選取塔爾湖地區(qū)1973、1987、1995、2006和2015年5個(gè)時(shí)期的遙感影像數(shù)據(jù)，同時(shí)根據(jù)現(xiàn)場(chǎng)勘察驗(yàn)證，采用人機(jī)交互圖像處理方法，獲取了塔爾湖地區(qū)5個(gè)時(shí)期土地利用/覆蓋現(xiàn)狀及土地利用/覆蓋類型的動(dòng)態(tài)變化特征。并且從土地利用程度、變化速度等方面分析了42年來該地區(qū)土地利用/覆蓋變化規(guī)律及其可能的自然與人為影響因素。

1數(shù)據(jù)來源與研究方法

1.1研究區(qū)概況 河套平原位于內(nèi)蒙古西南部，屬黃河沖積平原，地面開闊平坦。塔爾湖鎮(zhèn)在河套平原中北部，北部為色爾騰山，南部為黃河（圖1）。河套農(nóng)業(yè)區(qū)東西長(zhǎng)250 km，南北寬約50 km ，海拔1 040 m 左右，總面積16 043 km2，地表出露第四系河湖相

交替沖積層，土壤以粉質(zhì)黏土、黏土為主。加之該地區(qū)豐富的黃河地表水資源，為該地區(qū)農(nóng)業(yè)發(fā)展提供了有力的條件，是內(nèi)蒙古西部農(nóng)業(yè)經(jīng)濟(jì)發(fā)展的主要區(qū)域。研究區(qū)位于東亞季風(fēng)邊緣，屬大陸性干旱半干旱氣候，多年平均氣溫 7.1 ℃，多年平均降雨量 140 mm左右，多年平均蒸發(fā)量大于2 000 mm，植被類型植被以隱域性的草甸為主，地帶性植被為草原化荒漠和荒漠。受灌溉滲漏影響，灌區(qū)地下水埋深較淺，年平均地下水埋深為 1.5 m左右。由于本區(qū)降水少，蒸發(fā)量大，加上地下水埋深淺，使得本地區(qū)成為我國(guó)土壤鹽堿化和沙化發(fā)育的典型地區(qū)[16]。

1.2數(shù)據(jù)來源主要采用1973年MSS，1987、1995、2006年Landsat TM，2015年Landsat ETM+和2015年SPOT6時(shí)相的遙感影像。

1.3研究方法參照《土地利用現(xiàn)狀分類》國(guó)家標(biāo)準(zhǔn)[17]，并結(jié)合研究區(qū)的土地利用特點(diǎn)，把研究區(qū)的土地利用/覆蓋類型劃分為4個(gè)一級(jí)類型和8個(gè)二級(jí)類型（表1）。一級(jí)類型包括耕地、水體、城鄉(xiāng)居民點(diǎn)及工礦用地和未利用土地。

前期遙感圖像處理主要在ENVI 4.8上進(jìn)行，并選用ArcGIS 10.0作為遙感數(shù)字化解譯平臺(tái)。首先統(tǒng)一了坐標(biāo)系和投影參數(shù)，運(yùn)用ArcGIS對(duì)影像進(jìn)行大氣校正和幾何校正，鑲嵌并裁剪得到研究區(qū)的遙感影像數(shù)據(jù)。對(duì)研究區(qū)影像進(jìn)行圖像處理，包括彩色合成、線性增強(qiáng)、飽和度拉伸、主成分變換等；對(duì)圖像處理結(jié)果影像進(jìn)行計(jì)算機(jī)監(jiān)督分類（神經(jīng)網(wǎng)絡(luò)分類）信息提取。結(jié)合初步解譯過程，建立遙感解譯標(biāo)志，以彩色的SPOT6影像為主，輔之以黑白DOM數(shù)據(jù)、ETM+、TM影像，在ArcGIS 10.0軟件中建庫，進(jìn)行遙感目視解譯，然后根據(jù)野外驗(yàn)證結(jié)果反復(fù)修改解譯圖件，得到對(duì)應(yīng)年份的土地利用/覆蓋矢量圖。經(jīng)分類精度檢驗(yàn)得到各地類的Kappa系數(shù)均超過0.80，表明解譯結(jié)果能滿足研究需要。

2結(jié)果與分析

2.12015年塔爾湖地區(qū)土地利用/覆蓋現(xiàn)狀

總體上看，2015年塔爾湖地區(qū)的土地利用/覆蓋格局以耕地為主，占總面積的77.70%，達(dá)到278.43 km2。其中優(yōu)良耕地面積為138.39 km2，占到總面積的38.62%，輕度鹽堿地次之，面積達(dá)122.77 km2，占總面積的34.26%。研究區(qū)未利用土地達(dá)53.79 km2，占土地總面積的15.01%。其中以沙地面積最大，為36.30 km2，占總面積的10.13%；其次為重度鹽堿地，面積為15.91 km2，占總面積的4.44%。城鄉(xiāng)居民點(diǎn)及工礦用地和水體所占比重較小。從耕地、城鄉(xiāng)工礦居民用地和未利用土地的比重看，研究區(qū)是以耕地為主的農(nóng)業(yè)區(qū)，其中鹽堿化土地比重較高（表2和圖2）。

2.2土地利用/覆蓋變化時(shí)間序列

通過遙感解譯得出各

從表3可以看出，1973—2015年塔爾湖地區(qū)的土地利用/覆蓋總體格局未發(fā)生明顯的變化，都是以耕地為主，其次為未利用土地。鹽堿地中以輕度鹽堿地為主，未利用土地中主要為沙地。1973年塔爾湖地區(qū)耕地占總面積的比例為69.90%，未利用土地占24.20%，水體和城鄉(xiāng)居民及工礦用地分別占3.12%和2.78%。此外，鹽堿地面積總占53.59%，其中輕度和中度鹽堿地分別占31.64%和20.12%，重度鹽堿地占1.83%。沙地占總面積的21.59%。經(jīng)過14年的時(shí)間，到1987年塔爾湖地區(qū)耕地面積增加到73.42%，未利用土地和水體面積分別下降到20.23%和2.78%，城鄉(xiāng)居民點(diǎn)及工礦用地面積也有所上升到3.57%。鹽堿地總面積明顯下降到44.89%，其中中度鹽堿地下降最明顯為12.61%，輕度和重度鹽堿地分別下降到30.8%和1.48%；沙地面積也減少到18.07%。到1995年，塔爾湖地區(qū)耕地面積和未利用土地都略微減少，分別為71.41%和19.80%，水體和城鄉(xiāng)居民及工礦用地面積明顯增加，分別為4.43%和4.36%。鹽堿地總面積進(jìn)一步下降到39.53%，其中輕度鹽堿地明顯下降為26.35%，中度鹽堿地也有所下降到11.36%，而重度鹽堿地有所增加為1.82%；沙地面積進(jìn)一步減少為16.92%。再經(jīng)過11年的時(shí)間，到2006年塔爾湖地區(qū)耕地和城鄉(xiāng)居民及工礦用地面積分別明顯上升到75.28%和4.81%，而未利用土地面積持續(xù)下降到18.02%，水體面積也明顯降低到1.89%。鹽堿地總面積持續(xù)下降到37.68%，其中中度鹽堿地明顯下降為5.41%，而輕度和重度鹽堿地分別上升到29.34%和2.93%；沙地面積進(jìn)一步減少為14.65%。2015年塔爾湖地區(qū)耕地和城鄉(xiāng)居民及工礦用地面積分別進(jìn)一步增加為77.70%和5.67%，而未利用土地面積明顯下降到15.01%，水體面積略微下降到1.62%。鹽堿地總面積有所上升到43.52%，其中輕度和重度鹽堿地分別上升到34.26%和4.44%，中度鹽堿地下降到4.82%?？傮w來說，從1973—2015年，塔爾湖地區(qū)耕地面積（除去1995年略微減少外）和城鄉(xiāng)居民及工礦用地面積呈逐步增加的趨勢(shì)。而未利用土地面積和水體面積（除去1995年稍微增加外）呈逐步減少的趨勢(shì)（圖3）。

2.3土地利用/覆蓋變化速率分析

通過引入土地利用/覆蓋面積變化率來反映研究區(qū)土地利用/覆蓋變化時(shí)間序列特征。土地利用/覆蓋面積變化率（AVR）是指研究區(qū)內(nèi)土地利用/覆蓋類型在單位時(shí)間內(nèi)的變化，是時(shí)間的函數(shù)，能夠很

由表4可知，1973—1987 年，耕地和城鄉(xiāng)工礦居民用地面積增加，分別增加了 12.61和 2.83 km2，城鄉(xiāng)工礦居民用地的面積變化率較大，為2.03%。耕地增加速度較慢，在 0.36%左右。未利用土地和水體面積減少，減少速率分別為1.18%和0.79%。這一階段城鄉(xiāng)工礦居民用地面積增加迅速，表明人類活動(dòng)加強(qiáng)。而未利用土地中的沙地以每年1.16%明顯減少，減小面積為12.61 km2，因此耕地面積的增加主要是人類活動(dòng)對(duì)未利用土地（主要為沙地）的開發(fā)利用。這一階段鹽堿地總面積呈減小趨勢(shì)，以每年1.16%的速度減小。鹽堿地減少面積為31.18 km2，輕度、中度和重度鹽堿地面積都減少，其中中度鹽堿地面積減少速率較大為 2.67%，減小面積為 26.91 km2。

1987—1995 年，水體面積增加速度強(qiáng)度最大，面積增加 5.93 km2，增加速率為 7.44%，其次為城鄉(xiāng)工礦居民用地的面積呈持續(xù)增加，在 8年時(shí)間里，面積整體增加了 2.83km2，增加速率為 2.77%，增速較上一時(shí)間段顯著提高。耕地面積減少7.20 km2，減少速率為0.34%，未利用土地面積也持續(xù)減少，減少面積1.56km2，減少速率為0.27%，較上一時(shí)間段放慢。這一階段鄉(xiāng)工礦居民用地面積持續(xù)迅速增加，表明人類活動(dòng)持續(xù)增強(qiáng)，未利用土地面積減少也證明了這一點(diǎn)。而同時(shí)耕地卻有所減少，主要是由水體面積迅速增加造成的。這一階段鹽堿地總面積持續(xù)減小，減小速率接近 1.5%，減小面積為 19.21km2，其中輕度和中度鹽堿地面積減少速率分別為1.81%和1.24%，而重度鹽堿地面積卻增加1.22 km2，增加速率為2.87%。

1995—2006 年，耕地面積由減變?cè)?，增加?13.86 km2，增速較大，為0.49%，城鄉(xiāng)工礦居民用地的面積進(jìn)一步增加1.61 km2，變化率為0.94%，增速較上一時(shí)間段放慢。水體面積顯著減少9.11 km2，減少速率為 5.21%。未利用土地面積進(jìn)一步減少6.37 km2，減少速率為0.82%，較上一時(shí)間段減少更快?？傊?，這一階段城鄉(xiāng)居民點(diǎn)及工礦用地面積增加，表明人類活動(dòng)進(jìn)一步增強(qiáng)，但增速有所放緩。水體面積減少和對(duì)未利用土地進(jìn)一步開發(fā)，使得該階段耕地面積顯著增加。這一階段總鹽堿地面積進(jìn)一步減少面積6.63 km2，減少速率為0.43%。鹽堿地中中度鹽堿地面積顯著減少21.32 km2，減少速率為 4.76%，較上一階段顯著提高。重度鹽堿地面積增加3.98 km2，增速明顯提高為5.54%；輕度鹽堿地面積由前2個(gè)階段的減少開始變?yōu)樵黾樱黾用娣e10.71 km2，增速為1.03%。

2006—2015 年，耕地面積持續(xù)增加8.67 km2，增速為0.36%，城鄉(xiāng)工礦居民用地面積也持續(xù)增加3.09 km2，變化率為1.99%，增速較上一階段顯著上升。水體面積略微減少0.97 km2，減少速率為1.59%。未利用土地面積顯著減少10.78 km2，減少速率為1.85%，較上一時(shí)間段減少更快?？傊@一階段鄉(xiāng)工礦居民用地面積增速明顯，表明人類活動(dòng)顯著增強(qiáng)。該階段耕地面積的增加，主要是人類活動(dòng)增強(qiáng)對(duì)未利用土地的大力開發(fā)開發(fā)。總鹽堿地面積由此前的減少轉(zhuǎn)為顯著增加，增加面積20.93 km2，增加速率為1.72%。其中主要為輕度鹽堿地面積顯著增加17.63 km2，增速為1.86%；重度鹽堿地面積也增加5.41 km2，增速提高為5.73%；中度鹽堿地面積持續(xù)減少，減少面積2.11 km2，減少速率為1.21%。

總體來說，塔爾湖地區(qū)耕地除在1987—1995年有所降低外，基本保持在0.4%的年速率增長(zhǎng)。城鄉(xiāng)居民點(diǎn)及工礦用地面積呈逐年增加的趨勢(shì)，其中1987—1995年增速最快。未利用土地面積逐年減少，其中沙地面積減少最多，并在2006—2015 年減少速率最大。塔爾湖地區(qū)土地鹽堿化總面積在1973—2006年穩(wěn)步降低，而在2006—2015年迅速轉(zhuǎn)為增長(zhǎng)，其主要是由于近9年來輕度鹽堿化土地面積迅速增加造成的。輕度鹽堿化土地面積由前22年的逐步減少轉(zhuǎn)為從1995年開始增加，2006—2015年變化速率最快；中度鹽堿地一直穩(wěn)步減少，在1995—2006年間減少速率最大；重度鹽堿化土地面積從1987年來不斷增加，2006—2015年其擴(kuò)張速度最快（圖4）。

2.4土地利用/覆蓋變化和鹽堿化驅(qū)動(dòng)力分析

土地利用/覆蓋變化是自然與人類活動(dòng)雙重驅(qū)動(dòng)因素綜合作用的結(jié)果，自然因素（氣候、土壤、水文等）是影響土地利用結(jié)構(gòu)的基本要素，而人類活動(dòng)（開墾、灌溉等）在短時(shí)間尺度內(nèi)的強(qiáng)烈程度與合理度對(duì)土地利用結(jié)構(gòu)將也會(huì)產(chǎn)生顯著影響[20-21]。塔爾湖地區(qū)為農(nóng)業(yè)區(qū)，其耕地的變化對(duì)人們生活有著較大的影響。1973—2015年塔爾湖地區(qū)耕地的增長(zhǎng)，主要由人類活動(dòng)增強(qiáng)對(duì)未利用土地的開發(fā)利用。而1987—1995年耕地的下降是由于自然環(huán)境因素中降雨的增多，使得水域面積增大，占據(jù)了一部分耕地。對(duì)未利用土地的開發(fā)主要是沙地，沙地面積近42年來不斷減少，其中在1987—1995和2006—2015年這2個(gè)十年左右時(shí)間內(nèi)減少顯著，而2006—2015年降幅最為明顯（圖4）。1987—1995年沙地減少主要是由人類活動(dòng)開發(fā)為耕地，而2006—2015年減少的沙地大多數(shù)開發(fā)成為了居民地村莊和工礦用地（表4）。水體面積總體上也是逐年減少，其主要是由于近42年來降雨量逐漸減少，加上人類對(duì)水源的過度利用開發(fā)，而1987—1995年水體面積的增加主要是由于1995年當(dāng)年降水量增加造成的（圖4）。

土地鹽堿化程度與自然環(huán)境因素和人類對(duì)土地的耕作方式有著直接的關(guān)系[22-23]。通過將

1973—2015年

塔爾湖地區(qū)氣象、水文多個(gè)因素序列（數(shù)據(jù)來自中國(guó)氣象局）與鹽堿地面積和變化速率序列進(jìn)行對(duì)比分析，得出影響鹽堿地變化的可能因素（圖4）。自然因素方面，塔爾湖地區(qū)主要土壤主要為古河道相砂質(zhì)沉積，這種土質(zhì)粒度顆粒較粗，孔隙度和滲透率高，當(dāng)?shù)叵滤宦癫剌^淺時(shí)，鹽類物質(zhì)更容易通過土壤毛細(xì)管運(yùn)移到地表，通過蒸發(fā)析出、凝結(jié)并覆于地表形成鹽堿土[16]。氣候方面，通過對(duì)比發(fā)現(xiàn)，土壤總鹽堿化比例與當(dāng)年的相對(duì)濕度和當(dāng)年降水量呈明顯的正相關(guān)關(guān)系（圖4）。這可能是降水增加造成地下水水位抬升，使得溶解于地下水中的土壤鹽分運(yùn)移到表層，而較高的濕度表明蒸發(fā)強(qiáng)烈，因此，使得更多的鹽堿在表層土中析出，從而造成鹽堿化程度加重。人為因素方面，人類活動(dòng)尤其是農(nóng)村土地耕作方式，與土地鹽堿化關(guān)系密切，一方面，人類不合理的灌溉排水方式，導(dǎo)致地下水水位抬升，鹽分易隨地下水運(yùn)動(dòng)向地表聚積，加之蒸發(fā)強(qiáng)烈，加重土地鹽堿化，這從研究區(qū)重度鹽堿化土地與黃河引水渠空間分布關(guān)系可以看出[24]。另一方面人工固沙和鹽堿灘、荒漠的開發(fā)治理在一定程度上改善了土地狀況，減輕了鹽堿程度。節(jié)水改造措施的實(shí)施，使得研究區(qū)土壤剖面鹽分類型向均布型發(fā)展，且整體上土壤鹽分有所降低[25]。

3結(jié)論

利用遙感手段和數(shù)理統(tǒng)計(jì)方法對(duì)比分析了塔爾湖農(nóng)業(yè)區(qū)1973—2015年的5期遙感解譯結(jié)果，并對(duì)土地利用/覆蓋變化的驅(qū)動(dòng)力進(jìn)行了探討，得出以下結(jié)論。

（1）42年來塔爾湖地區(qū)的土地利用/覆蓋格局未發(fā)生顯著的變化。土地類型以耕地為主，占總面積的近70%以上；其次為未利用土地，主要為沙地；鹽堿地面積占40%以上，其中以輕度鹽堿地為主。

（2）42年來塔爾湖地區(qū)耕地除在1987—1995年略微降低外，基本保持著0.4%的年速率增長(zhǎng)；城鄉(xiāng)工礦居民用地面積呈逐年增加的趨勢(shì)，其中1987—1995年增速最快；未利用土地面積逐年減少，其中沙地面積減少最多，并在2006—2015 年減少速率最大；水體面積除在1995稍微增加外，也呈逐步減少的趨勢(shì)。

（3）1973—2006年塔爾湖地區(qū)土地鹽堿化總面積穩(wěn)步降低，而在2006以后轉(zhuǎn)為增長(zhǎng)。輕度鹽堿化土地面積從1995年開始由之前逐步減少轉(zhuǎn)為增加，2006—2015年變化速率最快；中度鹽堿地一直穩(wěn)步減少，在1995—2006年減少速率最大；重度鹽堿化土地面積1987年來不斷增加，2006—2015年擴(kuò)張速度最快。

（4）塔爾湖地區(qū)土地利用/覆蓋變化受到自然和人為因素的雙重作用，1973—2015年研究區(qū)在土壤鹽堿化土地改良和沙漠化防治等方面取得良好的效果，但2006—2015年不合理灌溉和耕地的過度利用使得土地鹽堿化程度增加。

參考文獻(xiàn)

[1]

MOONEY H A，DURAIAPPAH A，LARIGAUDERIE A.Evolution of natural and social science interactions in global change research programs[J].PNAS，2013，110（S1）：3665-3672.

[2] STERLING S M，DUCHARNE A，POLCHER J.The impact of global land-cover change on the terrestrial water cycle[J].Nature climate change，2012，3（4）：385-390.

[3] TIAN H Q，CHEN G，ZHANG C，et al.Centuryscale response of ecosystem carbon storage to multifactorial global change in the Southern United States[J].Ecosystems，2012，15：674-694.

[4] MEYFROIDT P，LAMBIN E F，ERB K H，et al.Globalization of land use：Distant drivers of land change and geographic displacement of land use[J].Current opinion in environmental sustainability，2013，5（5）：438-444.

[5] 劉紀(jì)遠(yuǎn)，匡文慧，張?jiān)鱿?，?20 世紀(jì)80 年代末以來中國(guó)土地利用變化的基本特征與空間格局[J].地理學(xué)報(bào)，2014，69（1）：3-14.

[6] 劉紀(jì)遠(yuǎn)，張?jiān)鱿?，徐新良，?21 世紀(jì)初中國(guó)土地利用變化的空間格局與驅(qū)動(dòng)力分析[J].地理學(xué)報(bào)，2009，64（12）：1411-1420.

[7] 陳菁，謝曉玲.海峽西岸快速城市化中土地利用變化的影響因素[J].經(jīng)濟(jì)地理，2010，30（11）：1885-1889.

[8] 劉紀(jì)遠(yuǎn)，劉明亮，莊大方，等.中國(guó)近期土地利用變化的空間格局分析[J].中國(guó)科學(xué)（D輯·地球科學(xué)），2002，32（12）：1031-1040.

[9] 宋開山，劉殿偉，王宗明，等.1954 年以來三江平原土地利用變化及驅(qū)動(dòng)力[J].地理學(xué)報(bào)，2008，63（1）：93-104.

[10] 方秀琴，任立良，李瓊芳.近40 年老哈河流域土地利用變化監(jiān)測(cè)與分析[J].國(guó)土資源遙感，2012，24（2）：125-131.

[11] 賈科利，常慶瑞，張俊華.陜北農(nóng)牧交錯(cuò)帶沙漠化土地時(shí)空變化特征分析[J].中國(guó)沙漠，2009，29（2）：223-228.

[12] 李小玉，張峰，肖篤寧.石羊河流域中、下游綠洲景觀動(dòng)態(tài)變化的比較研究[J].水土保持學(xué)報(bào)，2004，18（5）：151-154.

[13] 鄭曉，王乃昂，李卓侖，等.1990-2005 年疏勒河流域土地利用/覆蓋變化分析[J].中國(guó)沙漠，2010，30（4）：857-861.

[14] 郭姝姝，阮本清，管孝艷，等.內(nèi)蒙古河套灌區(qū)近30年鹽堿化時(shí)空演變及驅(qū)動(dòng)因素分析[J].中國(guó)農(nóng)村水利水電，2016（9）：159-162，167.

[15] 郭嬌，王偉，葉浩，等.河套平原鹽漬化土地時(shí)空動(dòng)態(tài)變化及影響因子[J].南水北調(diào)與水利科技，2014，12（3）：59-64.

[16] 侯玉明，王剛，王二英，等.河套灌區(qū)鹽堿土成因、類型及有效的治理改良措施[J].現(xiàn)代農(nóng)業(yè)， 2011（1）：92-93.

[17] 陳百明，周小萍.《土地利用現(xiàn)狀分類》國(guó)家標(biāo)準(zhǔn)的解讀[J].自然資源學(xué)報(bào)，2007，22（6）：994-1003.

[18] 王秀蘭，包玉海.土地利用動(dòng)態(tài)變化研究方法探討[J].地理科學(xué)進(jìn)展，1999，18（1）：81-87.

[19] 朱會(huì)義，李秀彬.關(guān)于區(qū)域土地利用變化指數(shù)模型方法的討論[J].地理學(xué)報(bào)，2003，58（5）：643-650.

[20] VELDKAMP A，F(xiàn)RESCO L O.Reconstructing land use drivers and their spatial scale dependence for Costa Rica （1973 and 1984）[J].Agricultural systems，1997，55（1）：19-43.

[21] 章予舒，王立新，張紅旗，等.疏勒河流域土地利用變化驅(qū)動(dòng)力分析：以安西縣為例[J].地理科學(xué)進(jìn)展，2003，22（3）：270-278.

[22] 裘善文.吉林霍洮、兩河中下游地區(qū)土地鹽漬化的特征、成因及治理的實(shí)用技術(shù)研究[J].土壤通報(bào)，2001，32（S1）：18-22.

[23] 王嵩，馮平，李建柱.地下水生態(tài)環(huán)境控制指標(biāo)問題的研究現(xiàn)狀[J].干旱區(qū)資源與環(huán)境，2005，19（4）：99-103.

[24] 晉南，王學(xué)東，劉震寧.后套：“河神”締造米糧川[J].中國(guó)國(guó)家地理，2012（10）：334-344.

[25] 李彬，史海濱，妥德寶，等.節(jié)水改造前后土壤鹽分剖面特征及其空間分布：以內(nèi)蒙古河套灌區(qū)中游臨河區(qū)為例[J].干旱區(qū)研究，2015，32（4）：663-673.