劉 潤, 姜琦剛, 史鵬飛, 劉 斌, 張莉莉
(吉林大學(xué) 地球探測科學(xué)與技術(shù)學(xué)院, 吉林 長春 130026)
在過去的幾十年里,中國經(jīng)濟(jì)經(jīng)歷了快速增長過程,這一發(fā)展過程伴隨著生態(tài)破壞、消費(fèi)和工業(yè)生產(chǎn)的增加[1],經(jīng)濟(jì)發(fā)展廣泛導(dǎo)致了生物多樣性下降、水土流失和土地退化[2]。實(shí)時準(zhǔn)確的土地利用信息在土地變化、水文動態(tài)和氣候變化等研究中至關(guān)重要[3]。土地利用變化是全球環(huán)境變化的重要組成部分,土地利用程度與人類活動息息相關(guān),是自然因子和人類活動因子綜合作用的結(jié)果[4]。人類對于土地利用的改造是長期的持久的過程。隨著近年來出現(xiàn)的全球氣候變暖、極端天氣等環(huán)境問題,土地利用變化成為全球環(huán)境變化的主要原因,土地利用變化的生態(tài)環(huán)境效應(yīng)研究已經(jīng)成為地理研究的熱點(diǎn)[5]。楊清可等[6]構(gòu)建城市土地利用與生態(tài)環(huán)境效應(yīng)評價指標(biāo)體系,運(yùn)用灰色關(guān)聯(lián)評價模型探究城市土地利用與生態(tài)環(huán)境效應(yīng)的時空演變規(guī)律與耦合度變化;任涵等[7]基于InVEST模型,結(jié)合Landsat解譯土地利用數(shù)據(jù),從地形起伏度視角,對太行山淇河流域2000—2015年生境質(zhì)量進(jìn)行研究;秦艷麗等[8]基于西安市景觀格局變化分析,利用同心緩沖區(qū)探討城鄉(xiāng)梯度內(nèi)景觀指數(shù)與生態(tài)系統(tǒng)服務(wù)價值的特征及其相互關(guān)系。前人研究表明,土地利用方式的改變對環(huán)境質(zhì)量變化的影響日益顯著。土地利用的生態(tài)環(huán)境效應(yīng)是大尺度范圍內(nèi)長時間序列研究,依靠傳統(tǒng)手段難以滿足土地利用變化監(jiān)測宏觀性實(shí)時性準(zhǔn)確性的要求。近些年來,隨著遙感技術(shù)的蓬勃發(fā)展,基于遙感手段的對地觀測技術(shù)已經(jīng)越發(fā)成熟,遙感是在大時空尺度上監(jiān)測土地利用變化的有用工具[9],具有宏觀性、綜合性、多尺度、多層次的特點(diǎn)[10]。遙感和GIS技術(shù)已被廣泛應(yīng)用于土地利用變化和城市化研究中,利用多時相衛(wèi)星影像探測土地利用變化及其時空格局已被證明是一種有效的方法[11]。因此,本文選用Landsat數(shù)據(jù)作為數(shù)據(jù)源,作為地球資源與環(huán)境探測領(lǐng)域的代表,Landsat數(shù)據(jù)記錄了長時間的地表變化,空間分辨率較高,能夠滿足本文2000—2020年時間范圍內(nèi)對研究區(qū)土地利用變化研究精度要求。并以遼寧省撫順市為研究區(qū),以保護(hù)生態(tài)環(huán)境為出發(fā)點(diǎn),利用遙感手段監(jiān)測土地利用變化,對城市土地利用與生態(tài)環(huán)境效應(yīng)的時空演變特征進(jìn)行定量分析,探索退耕還林以來土地利用變化對撫順市生態(tài)環(huán)境效應(yīng)的影響,以期為該區(qū)土地利用結(jié)構(gòu)優(yōu)化提供參考依據(jù)。
撫順市是遼寧省轄地級市,下轄新?lián)釁^(qū)、望花區(qū)、東洲區(qū)、順城區(qū)、撫順縣、新賓滿族自治縣和清原滿族自治縣7個縣區(qū),地理位置為123°39′42″—125°28′58″E,41°14′10″—42°28′32″N之間,總面積11 271 km2,位于遼寧省東部,東與吉林省接壤,西距省會沈陽市45 km,北臨鐵嶺,南臨本溪,為長白山余脈,整體地勢為東南高、西北低,地貌特征多以山地為主,呈現(xiàn)“八山一水一分田”,山地平均海拔400~500 m,地處中溫帶,屬大陸性季風(fēng)氣候,市區(qū)位于渾河沖積平原上,三面環(huán)山。
本文以撫順市為研究區(qū),選取2000年9月Landsat 7影像、2010年6月Landsat 5影像和2020年10月Landsat 8三期影像,數(shù)據(jù)來源于美國地質(zhì)勘探局(United States Geological Survey, USGS)數(shù)據(jù)共享網(wǎng)站(https:∥earthexplorer.usgs.gov/),從圖像質(zhì)量角度考慮,選取影像的含云量均低于1%。NDVI年際數(shù)據(jù)來源于中國科學(xué)院資源環(huán)境科學(xué)與數(shù)據(jù)中心(http:∥www.resdc.cn/)。該數(shù)據(jù)集是基于SPOT衛(wèi)星遙感數(shù)據(jù),在月數(shù)據(jù)基礎(chǔ)上采用最大值合成法生成,空間分辨率為1 km。數(shù)據(jù)預(yù)處理主要包括利用ENVI軟件對三期遙感影像進(jìn)行輻射定標(biāo)和大氣校正,分別對Landsat 5影像進(jìn)行7,5,1波段組合,Landsat 7影像7,5,1波段組合,Landsat 8影像7,6,2波段組合,并按研究區(qū)范圍做拼接及剪裁處理(表1)。
表1 研究區(qū)2000-2020年遙感影像數(shù)據(jù)源
2.1.1 土地利用分類方法與精度驗(yàn)證 結(jié)合撫順市實(shí)際土地利用特征、影像光譜特征并參考《土地利用現(xiàn)狀分類(GB/T21010-2017)》,將土地利用類型分為城鎮(zhèn)用地、耕地、林地、未利用地和水域5種類型。利用eCognition軟件多尺度分割功能對三期遙感影像進(jìn)行分割,分割時賦所有波段權(quán)重為1,形狀因子0.5,緊湊度因子0.5,經(jīng)過多次試驗(yàn),選用分割尺度50為最優(yōu)分割結(jié)果。利用最鄰近監(jiān)督分類方法,通過目視解譯手工選擇分割后的對象作為訓(xùn)練樣本數(shù)據(jù)集,每類20~30個,選擇訓(xùn)練樣本時盡量使其均勻分布在整個研究區(qū)域,規(guī)避不同土地利用類型單元的邊界,實(shí)現(xiàn)土地利用類型的劃分。選取驗(yàn)證點(diǎn)對土地分類結(jié)果進(jìn)行精度驗(yàn)證,采用kappa系數(shù)、總體精度、用戶精度、生產(chǎn)者精度4個指標(biāo)對分類結(jié)果進(jìn)行說明。
2.1.2 土地利用動態(tài)度 土地利用單一動態(tài)度[12]描述一種土地利用類型一定時期內(nèi)的轉(zhuǎn)化強(qiáng)度
(1)
式中:Ua,Ub分別為期間初期間末某種土地利用類型面積;T為研究期間時長;K為研究期間內(nèi)某一土地利用類型動態(tài)度。
2.1.3 土地利用轉(zhuǎn)移矩陣 土地利用轉(zhuǎn)移矩陣能夠反應(yīng)研究期初至研究期末各種土地利用類型相互轉(zhuǎn)化的動態(tài)信息,轉(zhuǎn)移矩陣的表達(dá)式為
(2)
式中:S為面積;i,j(i,j=1,2,3,…,n)為轉(zhuǎn)移前后的土地利用類型;n為土地利用類型數(shù)。
2.1.4 土地利用程度綜合指數(shù) 根據(jù)莊大方等[13]的研究,將土地開發(fā)利用強(qiáng)度的綜合量化指標(biāo)定義為在土地利用分級賦值(表2)的基礎(chǔ)上進(jìn)行數(shù)學(xué)綜合,形成值為1—4的連續(xù)分布指數(shù),考慮處理方便,在分級賦值的基礎(chǔ)上乘上100,得到土地利用程度綜合指數(shù)。土地利用程度綜合指數(shù)越高,表示土地開發(fā)利用強(qiáng)度越大,反之土地受自然因素影響越大。將研究區(qū)內(nèi)不同利用級進(jìn)行分解合并后計(jì)算所占比例,按其各自權(quán)重計(jì)算貢獻(xiàn)值:
表2 土地利用程度分級賦值標(biāo)準(zhǔn)[13]
(3)
式中:L為土地利用程度綜合指數(shù);Ai為第i級土地利用程度分級指數(shù);Ci為第i級土地利用程度面積所占比例。
2.2.1 NDVI差值分析 歸一化植被指數(shù)NDVI定義為近紅外波段的反射值與紅光波段的反射值之差比上兩者之和,取值范圍為-1~1。
NDVI=(NIR-R)/(NIR+R)
(4)
式中:NIR為近紅外波段的反射值;R為紅光波段的反射值。
通過對3期NDVI作差值分析,得到2000—2010年和2010—2020年NDVI變化趨勢。差值大于0時表示該柵格的NDVI變化趨勢是增加的;反之則說明NDVI變化趨勢是減小的;差值等于0表明該柵格的NDVI值在10 a間沒有變化。
2.2.2 土地利用類型轉(zhuǎn)化的生態(tài)環(huán)境質(zhì)量指數(shù) 利用土地利用生態(tài)環(huán)境質(zhì)量指數(shù)[14]定量描述撫順市生態(tài)環(huán)境質(zhì)量的總體狀況。表達(dá)式為
(5)
式中:EVt為t時期環(huán)境質(zhì)量指數(shù); LUi為t時期第i種土地利用類型面積;Ci為第i種土地利用類型的生態(tài)質(zhì)量指數(shù)(表3);TA為研究區(qū)域的總面積;n為土地利用類型數(shù)。
表3 土地利用類型及其生態(tài)環(huán)境質(zhì)量指數(shù)[15]
2.2.3 土地利用類型轉(zhuǎn)化的生態(tài)貢獻(xiàn)率 本文利用土地利用類型轉(zhuǎn)化的生態(tài)貢獻(xiàn)率(LEI)[16]研究土地利用類型轉(zhuǎn)化引起的區(qū)域生態(tài)環(huán)境質(zhì)量改變。
LEI=(LEIt+1-LEt)×LA/TA
(6)
式中:LEt+1,LEt分別為變化末期和初期的土地利用類型的生態(tài)環(huán)境指數(shù); LA為土地利用類型的轉(zhuǎn)化面積; TA為研究區(qū)域總面積。
3.1.1 土地利用分類結(jié)果與精度驗(yàn)證 借助eCognition軟件進(jìn)行多尺度分割,利用最鄰近監(jiān)督分類方法得到撫順市2000,2010和2020年土地利用分類結(jié)果如圖1所示,分類結(jié)果統(tǒng)計(jì)詳見表4。從圖1可以看出,撫順市各土地利用類型的分布存在差異,城鎮(zhèn)用地集中分布在西北方向,在其余方向僅零星分布;耕地在西北部城鎮(zhèn)周邊分布較密集,其余方向主要呈線狀放射分布。
圖1 撫順市土地利用分類結(jié)果
在研究區(qū)范圍內(nèi)隨機(jī)選取地面真實(shí)點(diǎn)543個作為測試數(shù)據(jù)集,利用資源環(huán)境與數(shù)據(jù)云土地分類數(shù)據(jù)結(jié)果進(jìn)行驗(yàn)證。采用kappa系數(shù)、總體精度(overall accuracy, OA)、用戶精度(user’s accuracy, UA)、生產(chǎn)者精度(producer’s accuracy, PA)4個指標(biāo)對分類結(jié)果進(jìn)行說明。2000,2010和2020年kappa系數(shù)分別為0.879 5,0.847 7和0.871 4,平均kappa系數(shù)為0.866 2;2000,2010和2020年總體精度分別為91.71%,89.50%,91.16%,平均精度為90.79%,精度較高,滿足研究要求(圖2)。
圖2 研究區(qū)土地利用分類結(jié)果精度評價
分析表4可知,林地為撫順市土地利用的主導(dǎo)類型,在2000,2010和2020年所占比例分別為70.55%,74.57%和75.79%,占絕對優(yōu)勢,這與撫順市地理位置相對應(yīng),撫順地處長白山余脈,主要地形為山地。其次是耕地、城鎮(zhèn)用地、水域,3類土地利用類型所占比例依次遞減,未利用地所占比例最少。在土地利用數(shù)量變化中,城鎮(zhèn)用地面積在不斷增加,到2020年已增加到395.32 km2,所占比例超過3.5%,這一變化與城市快速發(fā)展相關(guān)。林地面積在20 a間逐漸上升,耕地面積逐漸減少,主要與退耕還林、禁止開墾荒地等政策有關(guān)。
表4 研究區(qū)土地利用分類結(jié)果
3.1.2 土地利用動態(tài)度 根據(jù)3期影像分類結(jié)果,利用公式(1)計(jì)算撫順市單一土地利用動態(tài)度(表5)。
表5 撫順市2000-2020年單一土地利用動態(tài)度
土地利用動態(tài)度結(jié)果表明,2000—2020年撫順市土地利用類型變化最快的是城鎮(zhèn)用地及水體,其中城鎮(zhèn)用地動態(tài)度從2000—2010年的0.47增加到2010—2020年的7.73,表明研究期間撫順市處于快速城市化過程,城市擴(kuò)張顯著,20 a間研究區(qū)城鎮(zhèn)人口數(shù)量從149萬增加到155.2萬,增長率為4.16%;農(nóng)村人口向城市遷移,城市人口不斷增加;研究期間耕地面積以每年1%~2%的速度勻速減少,林地以0%~1%的速度勻速增加,耕地林地的單一動態(tài)度變化較小,變化比較穩(wěn)定;水域動態(tài)度變化較大,面積先增加后減少,可能與降雨和水庫豐水期缺水期有關(guān)。
3.1.3 土地利用轉(zhuǎn)移矩陣 土地利用變化包括不同土地利用類型的空間位置轉(zhuǎn)換和數(shù)量變化。建立土地利用轉(zhuǎn)移矩陣分析土地利用轉(zhuǎn)化的趨勢,撫順地區(qū)土地利用類型轉(zhuǎn)化明顯,總體來說(表6,圖3),林地和耕地之間相互轉(zhuǎn)化最為明顯,2000—2010年有1 213.74 km2的耕地轉(zhuǎn)化為林地,同時有716.33 km2的林地轉(zhuǎn)為耕地,最終由耕地轉(zhuǎn)化導(dǎo)致林地凈增加497.41 km2;2010—2020年有770.15 km2的耕地轉(zhuǎn)化為林地,同時有589.43 km2的林地轉(zhuǎn)為耕地,最終由耕地轉(zhuǎn)化導(dǎo)致林地凈增加180.72 km2。其次,全國范圍內(nèi)農(nóng)村向城市遷移的趨勢明顯,同時城市規(guī)模和人口也在迅速增加。中國城市人口增長的主要特征是農(nóng)村人口向城市遷移,以及農(nóng)村和自然人口的就地城鎮(zhèn)化[17],撫順地區(qū)存在大量耕地、林地逐漸轉(zhuǎn)化為城鎮(zhèn)用地,尤其2010—2020年有161.67 km2耕地和67.81 km2林地轉(zhuǎn)化為城鎮(zhèn)用地,占轉(zhuǎn)化后城鎮(zhèn)用地超50%。
圖3 撫順市2000-2020年土地利用類型轉(zhuǎn)化
表6 撫順市2000-2020年土地利用類型面積轉(zhuǎn)移矩陣 km2
3.1.4 土地利用程度綜合指數(shù) 根據(jù)公式計(jì)算得到的2000,2010和2020年土地利用程度綜合指數(shù)分別為224.48,225.57和226.12,2000—2010年土地利用程度綜合指數(shù)增加1.09,2010—2020年增加0.55,相鄰時間段內(nèi)土地利用程度綜合指數(shù)增速分別為0.49%和0.24%(表7)。
表7 撫順市2000-2020年土地利用程度綜合指數(shù)分析
分析撫順市不同土地利用類型面積所占比例可知,土地開發(fā)利用強(qiáng)度的增速主要受兩方面因素的影響:林地的增加耕地的減少會導(dǎo)致土地開發(fā)利用強(qiáng)度減小,城鎮(zhèn)用地面積的增加會導(dǎo)致土地開發(fā)利用強(qiáng)度的增加,最終土地開發(fā)利用強(qiáng)度由多個因素共同影響控制。綜合指數(shù)的大小反映土地利用強(qiáng)度的高低,研究期間撫順市土地利用強(qiáng)度緩慢上升。
3.2.1 NDVI年際變化 分別對2000—2010,2010—2020年兩個時段進(jìn)行NDVI差值分析。2000—2010年差值的最大值為0.364,最小值為-0.22,均值為0.059 243;差值大于零,即NDVI值10 a間呈上升趨勢面積所占比例達(dá)97.16%;差值小于零,即NDVI值10 a間呈下降趨勢的面積占比2.84%。2010—2020年差值的最大值為0.608,最小值為-0.268,均值為0.021 706,差值大于零,即NDVI值10 a間呈上升趨勢面積占比78.51%;差值小于零,即NDVI值10 a間呈下降趨勢面積所占比例21.49%。20 a間差值均值均大于0,說明撫順市植被覆蓋度在研究期間整體呈上升趨勢,但整體增加趨勢有所放緩。研究區(qū)NDVI變化值的分布存在一定的空間異質(zhì)性,NDVI在10 a間呈現(xiàn)減小趨勢的地區(qū)主要分布在研究區(qū)西北部,與城鎮(zhèn)擴(kuò)張區(qū)域基本重合,其余地區(qū)NDVI基本呈上升趨勢,只是上升幅度有所差異(圖4)。
圖4 撫順市2000-2020年NDVI時段差值分析
3.2.2 土地利用類型轉(zhuǎn)化的生態(tài)環(huán)境質(zhì)量指數(shù) 根據(jù)土地利用類型及其生態(tài)環(huán)境質(zhì)量指數(shù)表,對耕地、林地、城鎮(zhèn)用地、未利用地和水域分別賦予不同環(huán)境質(zhì)量指數(shù),得到撫順市2000,2010和2020年土地利用類型轉(zhuǎn)化的生態(tài)環(huán)境質(zhì)量指數(shù)分別為0.744 8,0.775 4和0.782 6;生態(tài)環(huán)境質(zhì)量指數(shù)變化分別為0.030 6和0.007 2,變化率分別為4.11%和0.93%,20 a間撫順市生態(tài)環(huán)境質(zhì)量總體呈上升趨勢,但生態(tài)環(huán)境質(zhì)量上升趨勢逐漸變緩。
3.2.3 土地利用類型轉(zhuǎn)化的生態(tài)貢獻(xiàn)率 根據(jù)公式(6)計(jì)算得到不同土地利用類型之間轉(zhuǎn)化對生態(tài)環(huán)境質(zhì)量指數(shù)的影響,并計(jì)算其導(dǎo)致生態(tài)環(huán)境改善或惡化的貢獻(xiàn)率及貢獻(xiàn)率所占比例(表8—10)。2000—2010年導(dǎo)致生態(tài)環(huán)境改善的地類轉(zhuǎn)化貢獻(xiàn)率前3名為耕地轉(zhuǎn)化為林地,貢獻(xiàn)率所占比例為93.67%,耕地轉(zhuǎn)化為水域,貢獻(xiàn)率所占比例2.55%,城鎮(zhèn)轉(zhuǎn)化為林地,貢獻(xiàn)率所占比例2.24%,此3類地類轉(zhuǎn)化貢獻(xiàn)率所占比例之和為98.46%,尤其耕地轉(zhuǎn)化為林地對生態(tài)環(huán)境改善的貢獻(xiàn)率超過90%;導(dǎo)致生態(tài)環(huán)境惡化的地類轉(zhuǎn)化貢獻(xiàn)率前3名為林地轉(zhuǎn)化為耕地,貢獻(xiàn)率所占比例89.16%,林地轉(zhuǎn)化為城鎮(zhèn)用地,貢獻(xiàn)率所占比例7.05%,林地轉(zhuǎn)化為水域,貢獻(xiàn)率所占比例2.17%,此3類地類轉(zhuǎn)化貢獻(xiàn)率所占比例之和為98.38%。2010—2020年導(dǎo)致生態(tài)環(huán)境改善的地類轉(zhuǎn)化貢獻(xiàn)率前3名為耕地轉(zhuǎn)化為林地,貢獻(xiàn)率所占比例93.66%,城鎮(zhèn)轉(zhuǎn)化為林地,貢獻(xiàn)率所占比例2.57%,水域轉(zhuǎn)化為林地,貢獻(xiàn)率所占比例1.61%,此3類地類轉(zhuǎn)化貢獻(xiàn)率所占比例之和為97.84%;導(dǎo)致生態(tài)環(huán)境惡化的地類轉(zhuǎn)化貢獻(xiàn)率前三名為林地轉(zhuǎn)化為耕地,貢獻(xiàn)率所占比例82.92%,林地轉(zhuǎn)化為城鎮(zhèn)用地,貢獻(xiàn)率所占比例10.22%,水域轉(zhuǎn)化為耕地,貢獻(xiàn)率所占比例2.48%,此3類地類轉(zhuǎn)化貢獻(xiàn)率所占比例之和為95.62%。綜合2000—2010年、2010—2020年土地利用類型的生態(tài)效應(yīng),分析2000—2020年總體土地利用類型對生態(tài)環(huán)境質(zhì)量貢獻(xiàn)率。20 a間,耕地向林地轉(zhuǎn)化是導(dǎo)致生態(tài)環(huán)境改善的最主要原因,其貢獻(xiàn)率超過93%;而生態(tài)環(huán)境惡化主要是林地向其它用地類型轉(zhuǎn)化造成的,其中林地轉(zhuǎn)化為耕地貢獻(xiàn)率所占比例為78%,林地轉(zhuǎn)化為城鎮(zhèn)用地貢獻(xiàn)率所占比例為19%。
表8 撫順市2000-2010年影響生態(tài)環(huán)境質(zhì)量的土地利用類型及其貢獻(xiàn)率
表9 撫順市2010-2020年影響生態(tài)環(huán)境質(zhì)量的土地利用類型及其貢獻(xiàn)率
表10 撫順市2000-2020年影響生態(tài)環(huán)境質(zhì)量的土地利用類型及其貢獻(xiàn)率
(1) 2000—2020年在撫順市土地利用類型中,林地面積所占比例超過70%,占主導(dǎo)地位,耕地面積占20%左右,城鎮(zhèn)用地、水域、未利用地面積所占比例較少;在2000—2020年研究區(qū)各種土地利用類型存在相互轉(zhuǎn)化關(guān)系,其中,林地、耕地之間的相互轉(zhuǎn)化與其它土地利用類型轉(zhuǎn)化為城鎮(zhèn)用地兩種現(xiàn)象最為明顯,20 a間研究區(qū)超過46%的城鎮(zhèn)用地面積由其它用地類型轉(zhuǎn)化而來。
(2) 20 a間研究區(qū)土地利用動態(tài)度結(jié)果表明,耕地面積不斷減少,林地和城鎮(zhèn)用地面積不斷增加,水域及未利用地面積呈現(xiàn)先增后減的波動變化,其中,林地和耕地面積變化速度略有下降,城鎮(zhèn)用地增速有所上升,說明退耕還林、禁止開墾荒地等保護(hù)政策對環(huán)境的影響效應(yīng)逐漸降低,且研究區(qū)城鎮(zhèn)化進(jìn)程飛速發(fā)展。多種因素綜合影響下研究區(qū)土地開發(fā)利用強(qiáng)度逐漸上升。
(3) NDVI差值分析結(jié)果表明,研究區(qū)2000—2010年97.16%的區(qū)域NDVI變化呈增加趨勢,2010—2020年78.51%的區(qū)域NDVI變化呈增加趨勢,20 a間植被覆蓋整體呈上升趨勢;NDVI減少的區(qū)域與城鎮(zhèn)用地的擴(kuò)張區(qū)域基本重合,符合快速推進(jìn)城鎮(zhèn)化進(jìn)程的同時環(huán)境質(zhì)量整體向好的特點(diǎn)。
本研究發(fā)現(xiàn)退耕還林政策實(shí)施以來,撫順市生態(tài)環(huán)境質(zhì)量指數(shù)逐漸上升,由土地利用變化導(dǎo)致的生態(tài)環(huán)境效應(yīng)整體趨勢向好。2000—2010年,生態(tài)環(huán)境質(zhì)量指數(shù)上升了4.11%,雖然研究區(qū)城鎮(zhèn)擴(kuò)張,但2000年左右開始實(shí)施的退耕還林還草政策效果顯著,這10 a間研究區(qū)生態(tài)環(huán)境向好明顯。2010—2020年,生態(tài)環(huán)境質(zhì)量指數(shù)上升了0.93%,生態(tài)環(huán)境延續(xù)向好趨勢,但趨勢逐漸平緩。一方面,2010—2020年撫順市城鎮(zhèn)化速度前所未有,城鎮(zhèn)擴(kuò)張嚴(yán)重?cái)D占耕地和生態(tài)用地;另一方面,由退耕還林還草政策帶來的生態(tài)改善空間逐步減小,研究區(qū)生態(tài)環(huán)境向好趨勢不明顯。環(huán)境污染和生態(tài)退化是工業(yè)文明的副產(chǎn)品。要權(quán)衡經(jīng)濟(jì)發(fā)展與生態(tài)保護(hù)之間的關(guān)系,找到一個平衡點(diǎn),以研究區(qū)經(jīng)濟(jì)發(fā)展與生態(tài)保護(hù)“雙贏”為目標(biāo),調(diào)節(jié)人口、資源、環(huán)境間生態(tài)關(guān)系。本文研究土地利用產(chǎn)生的生態(tài)環(huán)境效應(yīng),鑒于土地利用的動態(tài)性、生態(tài)效應(yīng)驅(qū)動因素的復(fù)雜性,如何在研究區(qū)城鎮(zhèn)化迅速發(fā)展與原有環(huán)境保護(hù)政策影響逐步減弱的情況下探索未來可持續(xù)發(fā)展的最優(yōu)路徑,尋找新的平衡點(diǎn)還有待進(jìn)一步研究。