南水北調(diào)東線江蘇受水區(qū)土地利用變化模擬及生態(tài)安全評價(jià)

方國華， 周磊， 聞昕， 薛劉宇， 嚴(yán)春華

(1.河海大學(xué)水利水電學(xué)院，南京 210098； 2.江蘇省南水北調(diào)工程建設(shè)領(lǐng)導(dǎo)小組辦公室，南京 210098)

南水北調(diào)東線江蘇受水區(qū)土地利用變化模擬及生態(tài)安全評價(jià)

方國華1， 周磊1， 聞昕1， 薛劉宇2， 嚴(yán)春華1

(1.河海大學(xué)水利水電學(xué)院，南京 210098； 2.江蘇省南水北調(diào)工程建設(shè)領(lǐng)導(dǎo)小組辦公室，南京 210098)

利用多期土地利用遙感數(shù)據(jù)，分析南水北調(diào)東線江蘇受水區(qū)2000—2010年間土地利用變化情況； 采用CA-Markov模型對2010—2020年間該區(qū)土地利用情況的演變規(guī)律進(jìn)行模擬； 同時(shí)構(gòu)建基于“壓力-狀態(tài)-響應(yīng)”(pressure-state- response，PSR)模型的土地利用生態(tài)安全評價(jià)指標(biāo)體系，綜合評價(jià)受水區(qū)不同時(shí)期土地利用生態(tài)安全狀況。結(jié)果表明： 1)受水區(qū)土地利用類型以耕地為主，研究期間林地和水域面積有所擴(kuò)大，大量耕地轉(zhuǎn)變?yōu)槌青l(xiāng)建設(shè)用地； 2)基于CA-Markov模型對研究區(qū)2010年土地利用格局的模擬結(jié)果與實(shí)際進(jìn)行對比分析，模擬正確的柵格數(shù)達(dá)到了99.1%，Kappa系數(shù)達(dá)到0.99，能夠較好地反映受水區(qū)土地利用格局變化趨勢； 3)土地利用生態(tài)安全評價(jià)結(jié)果顯示2000年、2010年和2020年土地生態(tài)安全指數(shù)分別為0.51，0.68和0.73，分別達(dá)到預(yù)警、敏感和良好等級，受水區(qū)的人口、資源環(huán)境及社會經(jīng)濟(jì)開始協(xié)調(diào)發(fā)展，整體呈現(xiàn)好轉(zhuǎn)趨勢。

南水北調(diào)； 土地利用； CA-Markov； PSR模型； 生態(tài)安全

0 引言

土地資源為人類社會的生存與發(fā)展提供物質(zhì)基礎(chǔ)，土地利用/土地覆被變化(land-use and land-cover change，LUCC)是人類社會與自然環(huán)境相互作用最直接的表現(xiàn)形式[1]。20世紀(jì)80年代興起的全球環(huán)境變化研究對LUCC在內(nèi)的一系列課題做了比較充分的探討[2-5]，涉及土地利用的過程、LUCC的人類響應(yīng)以及綜合的全球和區(qū)域模型等方面。基于遙感信息源的LUCC模擬模型研究更為普遍，如SLEUTH模型[6]、Logistic模型[7-8]和CA-Markov模型[9]等。土地生態(tài)安全作為“生態(tài)安全”的一部分，關(guān)系到人類社會的可持續(xù)發(fā)展[10]，對其研究具有重要意義。

南水北調(diào)東線工程是為解決蘇、魯、冀和津等省市缺水問題的大型調(diào)水工程，以京杭大運(yùn)河為輸水干線，串聯(lián)調(diào)蓄湖泊，逐級提水北送，并于2013年12月實(shí)現(xiàn)南水北調(diào)東線一期工程正式通水。近年來，圍繞南水北調(diào)受水區(qū)土地利用變化及生態(tài)安全方面的研究主要集中在中線地區(qū)并取得了豐富的成果[11-12]，而對東線工程沿線土地利用及其生態(tài)環(huán)境方面的研究相對欠缺[13]。由于工程路線長、涉及范圍廣，在促進(jìn)沿線經(jīng)濟(jì)社會發(fā)展的同時(shí)帶來了一系列的生態(tài)環(huán)境問題[14-17]，因而總結(jié)南水北調(diào)東線受水區(qū)土地利用變化特點(diǎn)，分析區(qū)域生態(tài)安全水平的演變規(guī)律與發(fā)展趨勢，可以為南水北調(diào)受水區(qū)生態(tài)環(huán)境保護(hù)提供有益參考，對保障沿線用水安全具有重大的現(xiàn)實(shí)意義。

本文以南水北調(diào)東線江蘇受水區(qū)為研究區(qū)，利用多期遙感影像數(shù)據(jù)，采用CA-Markov模型模擬分析該區(qū)域2000—2010年間的土地利用變化特點(diǎn)及2010—2020年間的土地利用情況演變規(guī)律，并利用“壓力-狀態(tài)-響應(yīng)”(pressure-state-response，PSR)模型構(gòu)建土地利用生態(tài)安全評價(jià)指標(biāo)體系，綜合評價(jià)研究區(qū)不同時(shí)期的生態(tài)安全狀況，為確保江蘇受水區(qū)生態(tài)安全以及南水北調(diào)東線工程調(diào)水安全提供科學(xué)依據(jù)。

1 研究區(qū)概況

南水北調(diào)東線工程江蘇受水區(qū)位于E116.22°～119.47°，N32.24°～34.57°之間，該區(qū)是在原江水北調(diào)工程基礎(chǔ)上擴(kuò)大規(guī)模、向北延伸而成，南起江都水利樞紐，北至南四湖，總面積為44 019 km2。研究區(qū)包括淮安市、宿遷市、徐州市、連云港市、鹽城市阜寧縣，以及揚(yáng)州市江都區(qū)、高郵市和寶應(yīng)縣，總面積約占全省面積的42.7%。區(qū)內(nèi)地勢西北高東南低，大部分高程在10 m以下，以平原為主； 淮河以南為亞熱帶季風(fēng)濕潤氣候區(qū)，以北為暖溫帶半濕潤季風(fēng)氣候區(qū)，都有明顯的季風(fēng)氣候特征。區(qū)內(nèi)人口密集，工農(nóng)業(yè)發(fā)達(dá)，是江蘇省經(jīng)濟(jì)發(fā)展最具活力的地區(qū)之一。目前，該區(qū)總?cè)丝跒? 631.58萬人，占全省總?cè)丝诘?7%； 耕地面積約1.956×104km2，占全省耕地面積的40.9%； 糧食產(chǎn)量為1 275.43×107kg，占全省產(chǎn)量的43.4%； 國內(nèi)生產(chǎn)總值和工農(nóng)業(yè)總產(chǎn)值分別為2 766.2億元和1 589.0億元，分別占全省的15.2%和15.0%。研究區(qū)具體位置如圖1所示。

圖1 研究區(qū)位置示意圖

2 研究方法

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

本研究所采用的數(shù)據(jù)包括土地利用數(shù)據(jù)和經(jīng)濟(jì)社會數(shù)據(jù)。其中土地利用數(shù)據(jù)是2000年、2005年和2010年3期的中國1∶10萬比例尺土地利用現(xiàn)狀遙感監(jiān)測數(shù)據(jù)，數(shù)據(jù)來源于中國科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心(http: //www.resdc.cn)。該數(shù)據(jù)是在各期Landsat TM/ETM遙感影像的基礎(chǔ)上，通過人工目視解譯生成，綜合制圖準(zhǔn)確率達(dá)95%[18]。經(jīng)濟(jì)社會數(shù)據(jù)主要來自于江蘇省統(tǒng)計(jì)年鑒(2001年、2006年及2011年)、《江蘇省新型城鎮(zhèn)化與城鄉(xiāng)發(fā)展一體化規(guī)劃(2014—2020年)》以及各市近中期城市總體規(guī)劃報(bào)告。

2.2 土地利用格局模擬及生態(tài)安全模型構(gòu)建

2.2.1 CA-Markov土地利用格局模擬模型

2.2.1.1 元胞自動機(jī)(cellular automaton，CA)模型

CA模型在時(shí)間、空間和狀態(tài)上都是離散的，每個(gè)變量對應(yīng)有限個(gè)狀態(tài)。該模型實(shí)質(zhì)上是一種建模框架，由離散的元胞、有限的狀態(tài)、鄰域和規(guī)則等4個(gè)基本要素構(gòu)成。模型可表示為

(1)

2.2.1.2 Markov模型

Markov過程體現(xiàn)了“無后效性”的特點(diǎn)，Markov模型是在Markov過程理論基礎(chǔ)上形成的一種方法。土地利用變化研究過程中，各土地利用類型之間相互轉(zhuǎn)換的面積或比例表示為狀態(tài)轉(zhuǎn)移概率，可用公式(2)對土地利用結(jié)構(gòu)變化狀態(tài)進(jìn)行預(yù)測，其表達(dá)式為

St=St-1Pij,

(2)

式中：St和St-1分別為t和t-1時(shí)刻的狀態(tài)；Pij表示狀態(tài)轉(zhuǎn)移矩陣。

2.2.1.3 CA-Markov模型

CA-Markov模型汲取了CA 模型空間模擬和Markov 模型數(shù)量預(yù)測的優(yōu)勢，使2個(gè)模型優(yōu)勢互補(bǔ)，既能夠有效提高各土地類型轉(zhuǎn)換的模擬準(zhǔn)確性，又能有效模擬土地格局的空間分布，因而得到了廣泛的探討和應(yīng)用[19-20]。

2.2.2 土地利用生態(tài)安全評價(jià)指標(biāo)體系模型

2.2.2.1 基于PSR模型的評價(jià)指標(biāo)體系

土地利用生態(tài)安全包括土地自身生態(tài)系統(tǒng)的完整性、可持續(xù)性以及為人類提供可持續(xù)的生態(tài)服務(wù)，具有系統(tǒng)性和復(fù)雜性等特點(diǎn)，是一個(gè)融合人口、資源環(huán)境和社會經(jīng)濟(jì)等方面的復(fù)合系統(tǒng)[21-22]。土地生態(tài)安全對土地系統(tǒng)服務(wù)功能滿足人類活動及社會經(jīng)濟(jì)可持續(xù)發(fā)展提出了更高要求，對土地生態(tài)安全評價(jià)要遵循系統(tǒng)性、科學(xué)性和綜合性等原則及安全評價(jià)的結(jié)果，是國家和社會對土地資源合理規(guī)劃和有序開發(fā)的預(yù)警和依據(jù)。

PSR模型可綜合考慮人類的活動對資源環(huán)境和社會經(jīng)濟(jì)造成的壓力以及資源環(huán)境和社會經(jīng)濟(jì)所做出的反應(yīng)，并通過壓力指標(biāo)、狀態(tài)指標(biāo)和響應(yīng)指標(biāo)予以表征。本文建立了基于PSR模型的土地利用生態(tài)安全評價(jià)指標(biāo)體系，指標(biāo)體系由目標(biāo)層、準(zhǔn)則層、因素層和指標(biāo)層4個(gè)層次構(gòu)成。目標(biāo)層即土地利用生態(tài)安全為總目標(biāo)，準(zhǔn)則層由壓力指標(biāo)(P)、狀態(tài)指標(biāo)(S)和響應(yīng)指標(biāo)(R)組成。在指標(biāo)層中，按各項(xiàng)指標(biāo)對土地生態(tài)安全的影響和作用，采用正向和逆向指標(biāo)予以區(qū)分，詳見表1。

表1 土地生態(tài)安全評價(jià)指標(biāo)體系

① 注： 括號中數(shù)值為各層級指標(biāo)權(quán)重。

指標(biāo)權(quán)重是指各指標(biāo)因子對土地生態(tài)安全貢獻(xiàn)度的高低，計(jì)算方法主要有模糊層次分析法、熵權(quán)法和主成分分析法等。本文采用模糊層次分析法對指標(biāo)體系中各指標(biāo)權(quán)重進(jìn)行計(jì)算。該方法可有效克服層次分析法在構(gòu)造判斷矩陣時(shí)沒有考慮到人為判斷模糊性、在一致性檢驗(yàn)時(shí)過于復(fù)雜及實(shí)用程度不高等不足[23-25]。

指標(biāo)權(quán)重計(jì)算的具體步驟為： 首先構(gòu)造模糊判斷矩陣，然后確定初始權(quán)重，并去模糊化，再計(jì)算得到最終權(quán)重，結(jié)果如表1所示。

土地生態(tài)安全評價(jià)標(biāo)準(zhǔn)值是評價(jià)土地生態(tài)安全評價(jià)指標(biāo)優(yōu)劣的依據(jù)，代表該地區(qū)土地利用生態(tài)安全良好狀態(tài)下各指標(biāo)的理想值。本文以地區(qū)經(jīng)濟(jì)社會發(fā)展規(guī)劃、國內(nèi)相關(guān)技術(shù)標(biāo)準(zhǔn)等為主，參考國內(nèi)外較為先進(jìn)地區(qū)的相關(guān)指標(biāo)水平，結(jié)合江蘇省土地生態(tài)環(huán)境狀況實(shí)際，綜合確定各土地生態(tài)安全評價(jià)指標(biāo)的標(biāo)準(zhǔn)值。

2.2.2.2 指標(biāo)數(shù)據(jù)標(biāo)準(zhǔn)化

土地利用生態(tài)安全指標(biāo)體系由正、逆2種類型指標(biāo)構(gòu)成，且各指標(biāo)量綱均不相同，難以進(jìn)行統(tǒng)一比較研究。本文對各指標(biāo)相關(guān)評價(jià)數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化處理，計(jì)算方法為

(3)

(4)

式中:Xi為指標(biāo)的原始值；Ci為指標(biāo)i的標(biāo)準(zhǔn)值；Ri為標(biāo)準(zhǔn)化后的安全指標(biāo)。

2.2.2.3 生態(tài)安全綜合指數(shù)

生態(tài)安全綜合指數(shù)是各指標(biāo)及地區(qū)綜合生態(tài)安全水平的量化反映。在通過模糊層次分析法得到各指標(biāo)的權(quán)重后，可計(jì)算各個(gè)指標(biāo)的生態(tài)安全指數(shù)(Ei)與區(qū)域生態(tài)安全綜合指數(shù)(E)，即

Ei=RiWi,

(5)

(6)

式中:Ei為指標(biāo)i的生態(tài)安全指數(shù)；E為生態(tài)安全綜合指數(shù)；Wi為指標(biāo)i的權(quán)重。

利用上述公式計(jì)算得到的生態(tài)安全指數(shù)E處于0～1之間，但僅憑E難以對其生態(tài)安全狀態(tài)進(jìn)行直觀評價(jià)。

由于目前土地利用生態(tài)安全評價(jià)尚未形成統(tǒng)一的評價(jià)標(biāo)準(zhǔn)，本文在參照相關(guān)研究[26-27]的基礎(chǔ)上，采用非等間距的方法，結(jié)合南水北調(diào)東線江蘇受水區(qū)的特點(diǎn)，同時(shí)考慮到標(biāo)準(zhǔn)的先進(jìn)型和實(shí)用性要求，針對本研究區(qū)設(shè)立Ⅰ—Ⅴ共5級評價(jià)等級，如表2所示。

表2 土地利用生態(tài)安全分類等級

3 結(jié)果與分析

通過分析比較2000—2010年間研究區(qū)土地利用遙感資料，計(jì)算其轉(zhuǎn)移面積矩陣，分析該區(qū)土地利用變化情況。同時(shí)，構(gòu)建CA-Markov模型，對2010年該區(qū)土地利用情況進(jìn)行模擬與驗(yàn)證，并分別預(yù)測 2015和 2020 年該區(qū)土地利用空間變化情況。在此基礎(chǔ)上，構(gòu)建基于PSR模型的土地生態(tài)安全評價(jià)指標(biāo)體系模型，對2000—2020年間研究區(qū)土地利用生態(tài)安全狀況及其發(fā)展變化進(jìn)行綜合評估。

3.1 2000—2010年間土地利用變化情況分析

2000—2010年間耕地是該區(qū)主要土地利用類型，占總面積的66.0%～67.1%，居主導(dǎo)地位； 其次為城鄉(xiāng)建設(shè)用地和水域，分別約占17.5%和12.2%； 林地和草地為次要土地利用類型，所占比例均不足3%。在此期間，區(qū)內(nèi)耕地和草地面積分別減少了452 km2和11 km2，而林地增加了10 km2，水域增加36 km2，城鄉(xiāng)建設(shè)用地增加417 km2，2000—2010年間研究區(qū)土地利用變化情況如圖2所示。

(a) 2000年 (b) 2005年 (c) 2010年

圖22000年、2005年和2010年研究區(qū)土地利用格局

Fig.2Landusepatternsofthestudyareain2000,2005and2010

2000—2005年間該區(qū)各土地類型轉(zhuǎn)移矩陣如表3所示。

表3 2000—2005年間各土地利用類型轉(zhuǎn)移矩陣

研究區(qū)工業(yè)化和城鎮(zhèn)化的發(fā)展導(dǎo)致耕地被大量占用，而南水北調(diào)東線工程的建設(shè)極大改善了輸水沿線的水環(huán)境，各調(diào)蓄湖泊水位普遍上升，湖泊面積擴(kuò)大，期間該地區(qū)耕地面積的凈減少量和水域、城鄉(xiāng)建設(shè)用地面積的凈增加量之和相當(dāng)。

2005—2010年間該地區(qū)各土地類型轉(zhuǎn)移矩陣如表4所示。

表4 2005—2010年間各土地利用類型轉(zhuǎn)移矩陣

土地利用的顯著變化是耕地面積明顯減少和城鄉(xiāng)建設(shè)用地面積明顯增加，呈基本持平狀態(tài)，變化面積分別為2000—2005年期間的3.7和6.4倍，顯示研究區(qū)內(nèi)城鎮(zhèn)化進(jìn)程加快，經(jīng)濟(jì)社會發(fā)展迅猛。

3.2 2010年土地利用模擬與驗(yàn)證

在2000年、2005年2期土地利用數(shù)據(jù)基礎(chǔ)上，建立CA-Markov模型，對研究區(qū)2010年土地利用情況進(jìn)行模擬預(yù)測(圖3(b))，并與當(dāng)年實(shí)際土地利用格局(圖3(a))進(jìn)行對比。結(jié)果顯示，42 955個(gè)柵格單元得到正確模擬，占總柵格數(shù)的 99.1%，且該區(qū)土地利用預(yù)測的總體 Kappa系數(shù)達(dá)0.99，表明預(yù)測準(zhǔn)確率較高，模擬結(jié)果真實(shí)可信。

(a) 2010年實(shí)際土地利用格局 (b) 2010年模擬土地利用格局 (c) 模擬誤差

圖3基于CA-Markov模型的土地利用格局模擬與驗(yàn)證

Fig.3SimulationandverificationoflandusepatternbasedonCA-Markovmodel

3.3 2015年和2020年土地利用變化情況預(yù)測

基于2010年研究區(qū)實(shí)際土地利用格局，利用CA-Markov模型對2015年和2020年該區(qū)土地利用格局進(jìn)行預(yù)測，結(jié)果見圖4。

(a) 2015年 (b) 2020年

圖42015年和2020年研究區(qū)土地利用格局預(yù)測

Fig.4Landusepredictionofthestudyareain2015and2020

由圖4可以發(fā)現(xiàn)，到2020年，該區(qū)城鄉(xiāng)建設(shè)用地面積預(yù)計(jì)達(dá)到8 011 km2，較2010年增加70 km2； 耕地將減少到28 330 km2，減幅達(dá)74 km2； 林地、草地面積沒有變化，水域面積增加到5 340 km2，較2010年增加4 km2。預(yù)測結(jié)果表明，研究區(qū)的城鄉(xiāng)建設(shè)用地需求進(jìn)一步增加，在人多地少、耕地后備資源不足的情況下，城鎮(zhèn)建設(shè)用地?cái)U(kuò)增與保護(hù)耕地的矛盾會進(jìn)一步加劇，對研究區(qū)土地利用生態(tài)安全會構(gòu)成新的壓力。

3.4 2000—2020年間土地利用生態(tài)安全評估

南水北調(diào)東線工程的順利實(shí)施有效提高了研究區(qū)用水水平和供水保證率，為推進(jìn)研究區(qū)經(jīng)濟(jì)社會穩(wěn)定發(fā)展和城市化進(jìn)程提供了可靠的水源保證，是該區(qū)土地利用變化最主要的驅(qū)動因素。而該區(qū)農(nóng)業(yè)種植結(jié)構(gòu)的調(diào)整、工業(yè)投資環(huán)境的改善以及城鎮(zhèn)化進(jìn)程的加快，反過來又深刻影響著研究區(qū)的社會穩(wěn)定和可持續(xù)發(fā)展。因此，全面評估研究區(qū)土地利用生態(tài)安全的變化具有重要意義。

在對研究區(qū)土地利用格局分析的基礎(chǔ)上，計(jì)算該區(qū)土地利用生態(tài)安全綜合指數(shù)，評估2000—2020年間土地利用生態(tài)安全狀態(tài)的變化趨勢，評價(jià)結(jié)果見表5。

表5 研究區(qū)土地利用生態(tài)安全評估

從表5中可以看到，2000—2005年間，研究區(qū)的安全等級均為預(yù)警等級，到2010年提升為敏感等級，2015年以后達(dá)到良好等級，表明研究區(qū)土地利用的生態(tài)安全狀況得到逐步好轉(zhuǎn)，土地生態(tài)安全問題得到妥善處理。各指標(biāo)生態(tài)安全狀態(tài)見圖5。

圖5 2000—2020年間各指標(biāo)土地利用生態(tài)安全綜合指數(shù)

圖5也表明，2000—2020年間，研究區(qū)域人口壓力以及資源環(huán)境狀態(tài)將不會發(fā)生顯著變化，保持相對穩(wěn)定的水平； 資源環(huán)境壓力呈現(xiàn)減小趨勢，主要得益于化肥施用強(qiáng)度的降低； 2000—2010年間，社會經(jīng)濟(jì)壓力、資源環(huán)境響應(yīng)以及社會經(jīng)濟(jì)響應(yīng)有了較大增長，得益于近10 a經(jīng)濟(jì)社會的快速發(fā)展，城鎮(zhèn)化水平不斷提高，生活污水排放及工業(yè)廢水處理能力不斷提升。

2000—2020年間研究區(qū)各城市的土地利用生態(tài)安全狀態(tài)總體呈現(xiàn)逐步好轉(zhuǎn)的趨勢。2000年，研究范圍內(nèi)6個(gè)城市均處于“預(yù)警”等級狀態(tài)，2010年全部轉(zhuǎn)變?yōu)椤懊舾小钡燃壱陨希貏e揚(yáng)州市已達(dá)到“良好”等級； 2020年，預(yù)計(jì)研究區(qū)內(nèi)各城市均將達(dá)到良好等級，表明該區(qū)域的人口、資源環(huán)境和社會經(jīng)濟(jì)開始協(xié)調(diào)發(fā)展。具體如圖6所示。

(a) 2000年 (b) 2005年 (c) 2010年

(d) 2015年 (e) 2020年

圖6研究區(qū)土地利用生態(tài)安全格局

Fig.6Ecologicalsecuritypatternoflanduseinthestudyarea

4 結(jié)論

1) 2000—2010年間研究區(qū)土地利用情況分析結(jié)果顯示，耕地是研究區(qū)主要土地利用類型，占各類型土地總面積的66.0%～67.1%； 城鄉(xiāng)建設(shè)用地和水域面積分別約占17.5%和12.2%。近10 a來研究區(qū)內(nèi)林地和水域面積有所擴(kuò)大，大量耕地轉(zhuǎn)變?yōu)槌青l(xiāng)建設(shè)用地，區(qū)域內(nèi)經(jīng)濟(jì)社會快速發(fā)展，城鎮(zhèn)化水平加快。

2) 基于CA-Markov模型對研究區(qū)2010年土地利用格局的模擬結(jié)果與實(shí)際進(jìn)行對比分析，模擬正確的柵格達(dá)到了99.1%，Kappa系數(shù)達(dá)到0.99，表明CA-Markov模型能夠較好地反映研究區(qū)土地利用格局變化的趨勢。2015年和2020年土地利用格局模擬結(jié)果表明，2020年耕地面積較2010年將減少74 km2，城鄉(xiāng)建設(shè)用地將增加70 km2。

3)結(jié)合PSR模型構(gòu)建土地利用生態(tài)安全評價(jià)指標(biāo)體系，采用模糊層次分析法求解各指標(biāo)權(quán)重，并對研究區(qū)生態(tài)安全做出定量評價(jià)： 2000年、2010年和2020年研究區(qū)土地利用生態(tài)安全綜合指數(shù)分別為0.51，0.68和0.73，分別達(dá)到預(yù)警、敏感和良好等級。未來研究區(qū)的土地生態(tài)安全狀況將逐步好轉(zhuǎn)，社會經(jīng)濟(jì)與資源環(huán)境也將逐步協(xié)調(diào)發(fā)展。

[1] 劉紀(jì)遠(yuǎn),匡文慧,張?jiān)鱿?等.20世紀(jì)80年代末以來中國土地利用變化的基本特征與空間格局[J].地理學(xué)報(bào),2014,69(1):3-14. Liu J Y,Kuang W H,Zhang Z X,et al.Spatiotemporal characteristics,patterns and causes of land use changes in China since the late 1980s[J].Acta Geographica Sinica,2014,69(1):3-14.

[2] 李秀彬.全球環(huán)境變化研究的核心領(lǐng)域──土地利用/土地覆被變化的國際研究動向[J].地理學(xué)報(bào),1996,51(6):553-558. Li X B.A review of the international researches on land use/land cover change[J].Acta Geographica Sinica,1996,51(6):553-558.

[3] Verburg P H,Schot P P,Dijst M J,et al.Land use change modelling:Current practice and research priorities[J].GeoJournal,2004,61(4):309-324.

[4] 唐華俊,吳文斌,楊鵬,等.土地利用/土地覆被變化(LUCC)模型研究進(jìn)展[J].地理學(xué)報(bào),2009,64(4):456-468. Tang H J,Wu W B,Yang P,et al.Recent progresses of land use and land cover change(LUCC) models[J].Acta Geographica Sinica,2009,64(4):456-468.

[5] Edenhofer O,Pichs-Madruga R,Sokona Y,et al.IPCC,2014:Climate Change 2014:Mitigation of Climate Change.Contribution of Working Group III to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change[M].Cambridge:Cambridge University Press,2014.

[6] Mahiny A S,Clarke K C.Guiding SLEUTH land-use/land-cover change modeling using multicriteria evaluation:Towards dynamic sustainable land-use planning[J].Environment and Planning B,2012,39(5):925-944.

[7] Lin Y P,Chu H J,Wu C F,et al.Predictive ability of logistic regression,auto-logistic regression and neural network models in empirical land-use change modeling - a case study[J].International Journal of Geographical Information Science,2011,25(1):65-87.

[8] Hao C,Zhang J H,Li H Y,et al.Integration of multinomial-logistic and Markov-chain models to derive land-use change dynamics[J].Journal of Urban Planning and Development,2015,141(3):05014017.

[9] Akin A,Berberoglu S,Erdogan M A,et al.Modelling land-use change dynamics in a Mediterranean coastal wetland using CA-Markov chain analysis[J].Fresenius Environmental Bulletin,2012,21(2A):386-396.

[10]梁留科,張運(yùn)生,方明.我國土地生態(tài)安全理論研究初探[J].云南農(nóng)業(yè)大學(xué)學(xué)報(bào),2005,20(6):829-834. Liang L K,Zhang Y S,Fang M.Study on the construction of land- ecological safety system[J].Journal of Yunnan Agricultural University,2005,20(6):829-834.

[11]劉秋.南水北調(diào)中線水源區(qū)土地的利用生態(tài)安全分析——以陜西省為例[J].首都師范大學(xué)學(xué)報(bào):自然科學(xué)版,2015,36(1):95-98. Liu Q.The Analysis of the ecological safety of the water sources areas in the Middle Route of South-to-North Water Transfer Project：Take Shaanxi Province for example[J].Journal of Capital Normal University:Natural Science Edition,2015,36(1):95-98.

[12]張海濤,白景鋒.南水北調(diào)中線水源區(qū)農(nóng)用土地生態(tài)安全時(shí)空分析——以湖北省為例[J].南陽師范學(xué)院學(xué)報(bào),2012,11(3):69-73. Zhang H T,Bai J F.Spatial-temporal analysis of ecological security in the water source area of the middle line project for transferring water from south to north：A case study in Hubei Province[J].Journal of Nanyang Normal University,2012,11(3):69-73.

[13]田家怡,陸兆華,閆永利,等.南水北調(diào)東線工程山東段生態(tài)安全評價(jià)[J].南水北調(diào)與水利科技,2007,5(6):39-44. Tian J Y,Lu Z H,Yan Y L,et al.Evaluation of ecological security on Shandong Province in the Eastern Route of the South-to-North Water Transfer Project[J].South-to-North Water Transfers and Water Science and Technology,2007,5(6):39-44.

[14]張全發(fā),蘇榮輝,江明喜,等.南水北調(diào)工程及其生態(tài)安全:優(yōu)先研究領(lǐng)域[J].長江流域資源與環(huán)境,2007,16(2):217-221. Zhang Q F,Su R H,Jiang M X,et al.The South-to-North Water Transfer Project and its ecological security:The priorities[J].Resources and Environment in the Yangtze Basin,2007,16(2):217-221.

[15]Zhang Q F.The South-to-North Water Transfer Project of China:Environmental implications and monitoring strategy[J].JAWRA Journal of the American Water Resources Association,2009,45(5):1238-1247.

[16]Lin C,Suh S,Pfister S.Does South-to-North Water Transfer Reduce the environmental impact of water consumption in China?[J].Journal of Industrial Ecology,2012,16(4):647-654.

[17]Liu J,Wan J.Environmental impact of the South-to-North Water Transfer Project in China[C]//World Water and Environmental Resources Congress.ASCE,2001.

[18]徐新良,龐治國,于信芳.土地利用、覆被變化時(shí)空信息分析方法及應(yīng)用[M].北京:科學(xué)技術(shù)文獻(xiàn)出版社,2014:167-168. Xu X L,Pang Z G,Yu X F.Spatial-Temporal Pattern Analysis of Land Use/Cover Change:Methods and Applications[M].Beijing:Scientific and Technical Documentation Press,2014:167-168.

[19]陳愛玲,都金康.基于CA-Markov模型的秦淮河流域土地覆蓋格局模擬預(yù)測[J].國土資源遙感,2014,26(2):184-189.doi:10.6046/gtzyyg.2014.02.29. Chen A L,Du J K.Simulation and forecast of land cover pattern in Qinhuai River Basin based on the CA-Markov model[J].Remote Sensing for Land and Resources,2014,26(2):184-189.doi:10.6046/gtzyyg.2014.02.29.

[20]汪佳莉,吳國平,范慶亞,等.基于CA-Markov模型的山東省臨沂市土地利用格局變化研究及預(yù)測[J].水土保持研究,2015,22(1):212-216. Wang J L,Wu G P,Fan Q Y,et al.Change and prediction of the land use in Linyi City,Shandong Province,based on CA-Markov model[J].Research of Soil and Water Conservation,2015,22(1):212-216.

[21]馬克明,傅伯杰,黎曉亞,等.區(qū)域生態(tài)安全格局:概念與理論基礎(chǔ)[J].生態(tài)學(xué)報(bào),2004,24(4):761-768. Ma K M,Fu B J,Li X Y,et al.The regional pattern for ecological security(RPES):The concept and theoretical basis[J].Acta Ecologica Sinica,2004,24(4):761-768.

[22]張虹波,劉黎明.土地資源生態(tài)安全研究進(jìn)展與展望[J].地理科學(xué)進(jìn)展,2006,25(5):77-85. Zhang H B,Liu L M.Main progress and prospects of land resource ecological security research[J].Progress in Geography,2006,25(5):77-85.

[23]蘭繼斌,徐揚(yáng),霍良安,等.模糊層次分析法權(quán)重研究[J].系統(tǒng)工程理論與實(shí)踐,2006,26(9):107-112. Lan J B,Xu Y,Huo L A,et al.Research on the priorities of fuzzy analytical hierarchy process[J].Systems Engineering-Theory and Practice,2006,26(9):107-112.

[24]張吉軍.模糊層次分析法(FAHP)[J].模糊系統(tǒng)與數(shù)學(xué),2000,14(2):80-88. Zhang J J.Fuzzy analytical hierarchy process[J].Fuzzy Systems and Mathematics,2000,14(2):80-88.

[25]Li L,Shi Z H,Yin W,et al.A fuzzy analytic hierarchy process(FAHP) approach to eco-environmental vulnerability assessment for the Danjiangkou reservoir area,China[J].Ecological Modelling,2009,220(23):3439-3447.

[26]左偉,王橋,王文杰,等.區(qū)域生態(tài)安全評價(jià)指標(biāo)與標(biāo)準(zhǔn)研究[J].地理學(xué)與國土研究,2002,18(1):67-71. Zuo W,Wang Q,Wang W J,et al.Study on regional ecological security assessment index and standard[J].Geography and Territorial Research,2002,18(1):67-71.

[27]湯潔,朱云峰,李昭陽,等.東北農(nóng)牧交錯(cuò)帶土地生態(tài)環(huán)境安全指標(biāo)體系的建立與綜合評價(jià)——以鎮(zhèn)賚縣為例[J].干旱區(qū)資源與環(huán)境,2006,20(1):119-124. Tang J,Zhu Y F,Li Z Y,et al.Evaluation on ecological security of land resources in ecotone between farming and animal raising in northeastern China:A case study of Zhenlai County[J].Journal of Arid Land Resources and Environment,2006,20(1):119-124.

(責(zé)任編輯：陳理)

LandusechangeandecologicalsecurityevaluationofJiangsuProvinceintheeasternrouteoftheSouth-to-NorthWaterTransferProject

FANG Guohua1， ZHOU Lei1， WEN Xin1， XUE Liuyu2， YAN Chunhua1

(1.CollegeofWaterResourcesandHydropowerEngineering,HohaiUniversity,Nanjing210098，China; 2.JiangsuLeadingGroupOfficefortheSouthtoNorthWaterDiversionProject,Nanjing210098,China)

In this paper, the authors studied land use changes of the eastern route of the South to North Water Transfer Project in Jiangsu Province based on land use data from 2000 to 2010, and simulated and predicted land use patterns from 2010 to 2020 based on CA -Markov model. “Pressure-State-Response” (PSR) model was used to establish the evaluation index system of land ecological security in the hope of evaluating the ecological security of land use in different periods. Some conclusions have been reached: 1) The main land use type in the area was cultivated land and the area of forest land, water was enlarged, and a large amount of cultivated land was changed into urban and rural areas. 2) The simulation results of land use pattern in 2010 based on CA-Markov model indicated that the correct grid number reached 99.1% compared with the actual one and the Kappa coefficient reached 0.99, which could reflect change trend of land use pattern in the area. 3) The ecological safety of the region showed gradually improvement and, for the year 2000, 2010, 2020, the land ecological security index was 0.51, 0.68 and 0.73, respectively achieving the early warning, sensitivity and good grades, which indicated coordinate development of population, resources, environment and society.

South-to-North Water Transfer Project; LUCC; CA-Markov; PSR model; ecological security

10.6046/gtzyyg.2017.03.24

方國華,周磊,聞昕,等.南水北調(diào)東線江蘇受水區(qū)土地利用變化模擬及生態(tài)安全評價(jià)[J].國土資源遙感,2017,29(3)：163-170.(Fang G H,Zhou L,Wen X,et al.Land use change and ecological security evaluation of Jiangsu Province in the eastern route of the South-to-North Water Transfer Project[J].Remote Sensing for Land and Resources,2017,29(3)：163-170.)

2016-02-29；

2016-04-08

江蘇省水利科技項(xiàng)目“南水北調(diào)東線工程調(diào)蓄湖泊聯(lián)合優(yōu)化調(diào)度研究”(編號： 2014012)、中央高校基本科研業(yè)務(wù)費(fèi)專項(xiàng)資金項(xiàng)目(編號： 2015B28614)和江蘇省普通高校研究生科研創(chuàng)新計(jì)劃項(xiàng)目“江蘇省河道管理評價(jià)指標(biāo)體系及管理標(biāo)準(zhǔn)的研究”(編號： SJZZ15_0059)共同資助。

方國華(1964-),女,教授,博士生導(dǎo)師,主要從事水資源規(guī)劃、水利經(jīng)濟(jì)和工程管理等方面的研究。Email： hhufgh@126.com。

TP 79

： A

： 1001-070X(2017)03-0163-08