于靖明,楊 光,吳 際,李慶和,溫雅琴

(1.內(nèi)蒙古農(nóng)業(yè)大學 沙漠治理學院,內(nèi)蒙古 呼和浩特 010018;2.荒漠生態(tài)系統(tǒng)保護與修復國家林業(yè)和草原局重點實驗室,內(nèi)蒙古 呼和浩特 010018;3.內(nèi)蒙古自治區(qū)水利事業(yè)發(fā)展中心,內(nèi)蒙古 呼和浩特 010020;4.呼和浩特市水資源與河湖保護中心,內(nèi)蒙古 呼和浩特 010020)

生態(tài)環(huán)境質(zhì)量綜合評價以生態(tài)學理論為基礎,選取具有代表性和客觀性的評價指標及模型,在特定的時間尺度和空間范圍內(nèi),定性、定量地對生態(tài)環(huán)境優(yōu)劣狀況加以判別與分析[1-2]。遙感技術(RS,remote sensing)因具備多尺度、短周期、長時序等優(yōu)勢,結合地理信息系統(tǒng)(GIS,geographic information system)強大的空間數(shù)據(jù)統(tǒng)計解析能力,廣泛運用于區(qū)域生態(tài)環(huán)境質(zhì)量評測中[3-4]。生態(tài)環(huán)境質(zhì)量評價自20世紀60年代以來,國際上相繼涌現(xiàn)出眾多學者選取不同地域、不同指標、不同方法對生態(tài)環(huán)境質(zhì)量實施評測,Franke等[5]運用植被指數(shù)對草地動態(tài)度演變進行分析;Badreldin等[6]采用土壤調(diào)整植被指數(shù)(SAVI)進行荒漠化動態(tài)監(jiān)測;Coutts等[7]選用地表溫度對城市熱島效應進行評測;80年代起,學者們開始從單一角度、單一因子向多角度、多因素的生態(tài)環(huán)境質(zhì)量評測方向進行探究[4,8]。目前,基于各區(qū)域尺度生態(tài)環(huán)境質(zhì)量復合指數(shù)的評價方法與模型眾多,主要包括層次分析法(AHP)[9]、綜合評價指數(shù)法[10]、灰色關聯(lián)度法[11]、“壓力-狀態(tài)-響應”(PSR)模型[12]等,卻并未形成一套系統(tǒng)化、正規(guī)化的生態(tài)環(huán)境質(zhì)量評價體系。2006年我國生態(tài)環(huán)境部頒布《生態(tài)環(huán)境狀況評價技術規(guī)范》,該規(guī)程耦合土壤、植被、水、生物、環(huán)境等關鍵生態(tài)指標,構建生態(tài)環(huán)境狀況指數(shù)(EI,ecological index),并逐步對該指數(shù)進行完善[13-14]。生態(tài)環(huán)境狀況指數(shù)被廣泛應用于縣級-區(qū)域生態(tài)環(huán)境質(zhì)量評價中,可EI指數(shù)法存在權重設定依賴專家先驗知識、數(shù)據(jù)獲取難度較大、信息更新周期長等弊端[15]。2013年,徐涵秋[16]提出完全以遙感信息為基底,集納多因子指標的遙感生態(tài)指數(shù)模型(RSEI,remote sensing ecological index),有效避免了人為干預權重的設定和單一指標表征的片面性。目前RSEI模型已廣泛應用于城市群[17]、水系流域[18]、水土流失區(qū)[19]、沙漠地區(qū)[20]等多種氣候、不同區(qū)域尺度。采用文獻分析法對基于RSEI構建的生態(tài)環(huán)境質(zhì)量評價文章進行研究,發(fā)現(xiàn)多數(shù)研究者未將水體演變對生態(tài)環(huán)境質(zhì)量的影響納入評價系統(tǒng),水體在城市化進程與生態(tài)環(huán)境建設中發(fā)揮著至關重要的作用,探究水體變化對其所產(chǎn)生的影響效應,進而健全區(qū)域生態(tài)環(huán)境質(zhì)量評價體系。

本研究基于遙感生態(tài)指數(shù)(RSEI)模型對內(nèi)蒙古臨河風沙區(qū)生態(tài)環(huán)境質(zhì)量實施多時相、多指標監(jiān)測評價,從時間、空間2個維度,演變面積、演變方向、演變幅度3個方面,全面分析研究區(qū)1989-2021年間生態(tài)環(huán)境質(zhì)量演變進程,同時,考慮到RSEI模型評價成果易受大面積水域影響,對臨河區(qū)水體實施掩膜處理,單獨分析水體演變對生態(tài)環(huán)境質(zhì)量的影響效應,以期為內(nèi)蒙古風沙區(qū)綜合治理、生態(tài)恢復及高質(zhì)量發(fā)展提供參考。

1 材料與方法

1.1 研究區(qū)概況

位于內(nèi)蒙古風沙區(qū)腹地的臨河區(qū),地處107°6′-107°44′E,40°34′-44°17′N,區(qū)域面積2 333 km2,北靠狼山山脈,東處河套平原,近半數(shù)區(qū)界同沙漠相連,南與庫布奇沙漠依河而劃,西接烏蘭布和沙漠,且同騰格里、巴丹吉林沙漠相望;典型地帶性土壤為棕鈣土、灰棕荒漠土、風沙土等,植被覆蓋度低,存在較嚴重的風力侵蝕;典型氣候為溫帶大陸性季風氣候,年均降水量僅138.8 mm,而蒸發(fā)量達2 236.7 mm,春秋兩季,河套平原是西伯利亞冷空氣途徑蒙古國、南下我國的必經(jīng)之地,大風過境極易形成沙塵天氣,據(jù)統(tǒng)計,臨河區(qū)年均受沙塵襲擊天數(shù)達20 d,所釀成的財產(chǎn)損失每年高達數(shù)億元人民幣[21-22]。

圖1 研究區(qū)區(qū)位Fig.1 Location map of the study district

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

采用的8景遙感影像數(shù)據(jù)均源于地理空間數(shù)據(jù)云(https://www.gscloud.cn),分別為1989、2000、2010年Landsat 5 TM影像和2021年Landsat 8 OLI影像,行列號為129/031及129/032。選取6-8月影像質(zhì)量良好、云量偏小、季相一致的數(shù)據(jù),以提升研究結果的精準性。

基于ENVI 5.3軟件對各期遙感影像實施預處理,減少不同影像間鑒于大氣、云量、地形差異所造成的輻射誤差。參考Chander等[23]和Chavez[24]的模型參數(shù)對多光譜波段及熱紅外波段分別進行輻射定標;利用FLAASH模型實施大氣校正;運用最鄰近像元法進行幾何校正,使得配準的均方根誤差(RMSE)小于1/2單位像元;實施影像鑲嵌裁剪。

1.3 研究方法

1.3.1 指標構建 RSEI模型耦合綠度、干度、濕度、熱度4個能夠直觀映射生態(tài)環(huán)境優(yōu)劣狀況的關鍵指標,采用多元統(tǒng)計方法中的主成分分析法(PCA),通過垂直旋轉坐標軸的方式將多維信息匯集到幾個特征分量上[25-26],可有效規(guī)避權重設定的主觀意識性,進而客觀、精準地評測研究區(qū)生態(tài)環(huán)境質(zhì)量。

1.3.2 模型指標

1.3.2.1 綠度指標 選取被廣泛運用的歸一化植被指數(shù)(normalized difference vegetation index,NDVI)代表綠度,可以快速反映地表植被的整體狀況,且與植被覆蓋度FVC存在極為凸顯的線性關聯(lián)[27]。

(1)

式中:ρNIR、ρRed為近紅外和紅波段的反射率。

1.3.2.2 濕度指標 選取纓帽變換后的濕度分量來表征濕度指標(WET),經(jīng)纓帽變換后獲得的濕度分量,有效消除了冗余數(shù)據(jù),可較好地反映植被及土壤的濕度狀況[28-29]?；诓煌瑐鞲衅鞯腖andsat 8 OLI和Landsat 5 TM影像,濕度分量的轉換系數(shù)不同。

WETTM=0.031 5ρBlue+0.202 1ρGreen+0.310 2ρRed+0.159 4ρNIR-0.680 6ρSWIR1-0.610 9ρSWIR2

(2)

WETTM=0.151 1ρBlue+0.197 3ρGreen+0.328 3ρRed+0.340 7ρNIR-0.711 7ρSWIR1-0.455 9ρSWIR2

(3)

式中:ρSWIR2、ρSWIR1、ρNIR、ρRed、ρGreen、ρBlue為短波紅外2、短波紅外1、近紅外、紅、綠、藍波段的反射率。

1.3.2.3 干度指標 造成研究區(qū)地表“干化”的主導成因為裸土地和建筑用地,采用建筑指數(shù)(IBI)和裸土指數(shù)(SI)的平均值表征干度指標(NDSI)。

(4)

(5)

(6)

式中:ρSWIR1、ρNIR、ρRed、ρGreen、ρBlue分別表示短波紅外1、近紅外、紅、綠、藍波段的反射率。

1.3.2.4 熱度指標 選用輻射傳輸方程法反演地表溫度來表征熱度指標(LST)[30]。參考覃志豪等[31]提出的地表比輻射率方程式,將影像大致分為水體、自然表面和城鎮(zhèn)用地3種類型,分類運算地表比輻射率。

水體:εwater=0.995

(7)

(8)

(9)

通過普朗克定律反函數(shù)獲取地表溫度(TS)。

(10)

(11)

式中:TS為地表溫度(℃);B(TS)為黑體熱輻射亮度;ε為地表比輻射率;FVC為植被覆蓋度;Lup、Ldown和τ分別為大氣上行輻射、大氣下行輻射和大氣在熱紅外波段的透過率。該數(shù)據(jù)獲取于美國國家航空航天局(NASA)網(wǎng)站(http://atmcorr.gsfc.nasa.gov)[29];K1和K2為定量參數(shù),在Landsat 5 TM和Landsat 8 OLI影像中,分別為K1=607.76W/(m2·sr·μm)、K2=1 260.56K和K1=774.89 W/(m2·sr·μm)、K2=1 321.08 K;L分別為Landsat 5 TM和Landsat 8 OLI影像中熱紅外6波段和熱紅外10波段處的輻射亮度值。

1.3.3 RSEI模型構建 運用極差變換法對上述各指標進行歸一化處理,消除由于量綱不統(tǒng)一所造成的權重失衡,使其數(shù)值映射到[0,1]區(qū)間[16]。繼而運用主成分分析法(PCA)集成以上4個指標,獲得的初始遙感生態(tài)指數(shù)(RSEI0)經(jīng)正負轉置、歸一化處理,更加精準、定量的展現(xiàn)區(qū)域生態(tài)環(huán)境總體質(zhì)量。

(12)

RSEI0=1-{PC1[f(NDVI,WET,NDSI,LST)]}

(13)

(14)

式中:NIi為指標歸一化后的像元值;Ii為像元i的指標值;Imax、Imin分別為指標的最大值與最小值;PC1為第1主成分分量;RSEImax、RSEImin為RSEI0的最大值與最小值。

值域介于[0,1]的遙感生態(tài)指數(shù)(RSEI),其數(shù)值越接近于1,生態(tài)環(huán)境質(zhì)量狀況越優(yōu),否之亦然。為便于遙感生態(tài)指數(shù)空間格局的定量化展示,本研究將各年份RSEI值以0.2為間距,劃分為:差[0,0.2]、較差[0.2,0.4]、中[0.4,0.6]、良[0.6,0.8]、優(yōu)[0.8,1]5個等級。

1.3.4 水體提取及掩膜 選用改進的歸一化差異水體指數(shù)(MNDWI,modified normalized difference water index),對各期NDVI、WET、NDSI、LST指標在歸一化前進行水體掩膜,減少大面積水域對主成分分析產(chǎn)生的荷載影響,精準反映區(qū)域生態(tài)環(huán)境質(zhì)量狀況[3]。

(15)

式中:ρGreen、ρSWIR1為綠波段和短波紅外1波段處的反射率。

驗證發(fā)現(xiàn),僅依據(jù)MNDWI指數(shù)進行水體掩膜,會造成大量建設用地被掩蓋及水體區(qū)域大幅增加。因此,采用支持向量機模型(SVM,support vector machine)、MNDWI指數(shù)和人工目視解譯相結合的方法,對研究區(qū)水體單獨進行時空演變分析。聯(lián)系實際需求,將臨河區(qū)土地利用類型重分類為水體、植被、建設用地和未利用地4類,對劃分結果進行各類別精度、生產(chǎn)精度、總精度檢驗,各精度驗證均>95%,滿足分析需求[3]。結合土地利用轉移矩陣(式16),分析1989-2021年水體與各土地利用類型及各等級生態(tài)環(huán)境質(zhì)量之間的相互轉換關系,完善區(qū)域生態(tài)環(huán)境質(zhì)量監(jiān)督評測體系。

(16)

式中:S為土地利用面積;i、j為研究初期和末期土地利用類型。

2 結果與分析

2.1 臨河區(qū)生態(tài)環(huán)境質(zhì)量總體狀況

分析各年份RSEI、NDVI、WET、NDSI、LST均值圖(圖2)。1989-2021年臨河區(qū)生態(tài)環(huán)境質(zhì)量以2010年為節(jié)點呈“先增長后下降”趨勢,RSEI均值總體由0.295降至0.254,降幅13.90%。至2021年,對生態(tài)環(huán)境起正面效應的綠度指標總體下降50.19%,與實際情況相符;與生態(tài)環(huán)境起負面效應的LST指標總體增長16.25%,反映1989-2021年臨河區(qū)在城市化進程中,涌現(xiàn)出植被持續(xù)減少、地表硬化、氣溫增高等多種生態(tài)環(huán)境惡化問題。

圖2 研究區(qū)1989-2021年RSEI、NDVI、WET、NDSI、LST均值Fig.2 Mean values of RSEI,NDVI,WET,NDSI and LST in the study area from 1989 to 2021

對研究區(qū)RSEI空間分布格局進行可視化分析(圖3),1989-2021年臨河區(qū)生態(tài)環(huán)境質(zhì)量總體呈現(xiàn)較差等級。1989年差等級區(qū)域在研究區(qū)北部呈斑點狀分布;2000年生態(tài)環(huán)境質(zhì)量等級為差的地域在中部沙地及城鎮(zhèn)化周邊呈斑面狀分布,中等級地域呈條帶狀分布于北部、中南部及黃河沿岸地區(qū);2010年臨河區(qū)全力落實林草目標任務,良等級生態(tài)環(huán)境質(zhì)量在臨河區(qū)西北部的農(nóng)田、林木覆蓋區(qū)首次呈現(xiàn)面狀分布;經(jīng)10 a的城鄉(xiāng)高速發(fā)展、工業(yè)園區(qū)開發(fā)、油田礦產(chǎn)資源采伐等,至2021年臨河區(qū)生態(tài)環(huán)境質(zhì)量明顯惡化,差等級區(qū)域在北部良田整改區(qū)、中部沙地、南部建成區(qū)及周邊地區(qū)呈現(xiàn)面狀分布。

圖3 研究區(qū)1989-2021年生態(tài)環(huán)境質(zhì)量空間分布Fig.3 Spatial distribution of eco-environmental quality in the study district from 1989 to 2021

2.2 空間演變

生態(tài)環(huán)境質(zhì)量空間演變分析包括對各等級RSEI空間面積轉變分析和重心遷移狀況分析。本研究從演變面積、演變方向、演變幅度三方面對研究區(qū)生態(tài)環(huán)境質(zhì)量空間演變進行闡述。

2.2.1 面積轉變 將各年份、各等級生態(tài)環(huán)境質(zhì)量面積變化狀況進行量化與匯總(表1)。整體來看,臨河區(qū)生態(tài)環(huán)境質(zhì)量以較差等級為主導,1989、2000、2010、2021年較差等級面積分別占研究區(qū)總面積的99.82%、90.97%、60.23%、93.66%。研究期間等級為差的區(qū)域面積,呈“波浪式”起伏,從1989-2000年上升35.94 km2,至2010年減少為11.76 km2,而2021年驟增為146.48 km2,對照1989年,差等級地域面積總體增多142.97 km2;中等級區(qū)域面積以2010年為節(jié)點呈“先增長后降低”趨勢,2010年中等級區(qū)域約占全區(qū)總面積的1/3以上,至2021年僅存1.35 km2,主導成因為建成區(qū)的迅速擴張和人類對自然資源的頻繁活動所致;優(yōu)、良等級區(qū)域面積在2010年達歷史最高值,分別為1.81 km2與153.59 km2。

表1 研究區(qū)1989-2021年RSEI各等級信息統(tǒng)計Table 1 Information statistics of all levels of RSEI in the study area from 1989 to 2021

總體而言,1989-2021年臨河區(qū)生態(tài)環(huán)境質(zhì)量呈下降趨勢,差等級區(qū)域面積整體增幅6.13%,而中、良、優(yōu)等級面積占比變化均<1%,警示臨河區(qū)在城鄉(xiāng)擴建進程中,應積極推進高質(zhì)量生態(tài)體系筑建工程。

2.2.2 重心遷移 為更加直觀、清晰地展現(xiàn)1989-2021年RSEI各等級空間格局變化,繪制重心遷移軌跡圖(圖4),結合表2分析得知,遷移距離可映射出各等級空間分布密集狀況,遷移方向可反映出各等級區(qū)域面積變化方向。

表2 1989-2021年重心遷移距離Table 2 Center of barycentre migration distance from 1989-2021

圖4 1989-2021年臨河區(qū)RSEI各等級重心遷移Fig.4 Center of gravity migration of all grades of RSEI in the Linhe District from 1989 to 2021

差等級RSEI空間重心遷移軌跡較長、面積變化顯著,1989、2000、2000-2010年差等級空間重心分別向南遷移30.16 km和11.96 km,主導原因為南部建城區(qū)的擴張和黃河沿岸硬化改造工程的實施,而2010-2021年間重心又向北部折返11.96 km,這是由于中部、北部城鎮(zhèn)化擴建迅速、高標準良田整改及狼山山脈金屬礦業(yè)的采伐所致;較差等級區(qū)域在研究區(qū)空間布局廣泛,面積變化幅度較小,因而遷移重心始終分布于研究區(qū)中心位置,總體來看,較差等級空間重心32 a間共向南遷移11.77 km;中等級RSEI重心先向西南遷移17.11 km,后向北偏東折返21.84 km,這與研究區(qū)南部中等級生態(tài)環(huán)境質(zhì)量先增長后下降的特征相一致;優(yōu)、良等級重心遷移方向基本一致,均先向北偏東,后向南偏東遷移,分別共遷移41.64 km和61.73 km,這與2010-2021年北部優(yōu)良等級區(qū)域面積大幅度縮減及南部黃河國家濕地公園的建立息息相關。

RSEI各等級空間重心分布狀況在2000-2010年由南向北呈自低到高規(guī)律分布,與1989、2021年空間重心分布狀況截然相反;從重心遷移軌跡來看,差、較差、良、優(yōu)等級總體向南遷移,中等級向北遷移;從遷移距離而言,優(yōu)、良、差等級相比較差、中等級遷移距離更遠。綜上所述得知,各等級空間分布越密集且面積變化率越大,重心遷移距離則越長。

2.3 時間演變

從縱向維度對臨河區(qū)生態(tài)環(huán)境質(zhì)量實施動態(tài)監(jiān)測,將各時段遙感生態(tài)指數(shù)(RSEI)差值分析結果劃分為9個等級差,揭示生態(tài)環(huán)境質(zhì)量變化顯著程度,正級差表征生態(tài)環(huán)境質(zhì)量得以改善,0表示穩(wěn)定未變化,負級差表征生態(tài)環(huán)境質(zhì)量正遭受惡化[22]。各研究階段生態(tài)環(huán)境質(zhì)量時間演變結果及面積統(tǒng)計信息見表3、圖5。

表3 各研究階段生態(tài)環(huán)境質(zhì)量演變信息統(tǒng)計Table 3 Information statistics of eco-environmental quality evolution in each research phase

圖5 1989-2021年時間尺度生態(tài)環(huán)境質(zhì)量演變Fig.5 Evolution of eco-environmental quality on a time scale from 1989 to 2021

1989-2000年臨河區(qū)生態(tài)環(huán)境質(zhì)量演變以穩(wěn)定為主線,以改善為趨勢,呈現(xiàn)惡化狀況區(qū)域均以惡化1個級差為主,并未出現(xiàn)跨級惡化現(xiàn)象。1989-2000年生態(tài)環(huán)境質(zhì)量呈穩(wěn)定等級的區(qū)域占臨河區(qū)總面積的91.16%,達歷年占比最大值,167.51 km2的區(qū)域表現(xiàn)為改善狀況,其中97.95%的區(qū)域改善1個級差,主要位于研究區(qū)西南部、中部以東和北部的種植區(qū)及防護林區(qū),1978年以來我國開展的“三北”

防護林體系建設工程及廣大人民群眾拓荒開田等工作成效顯著,使臨河區(qū)生態(tài)環(huán)境質(zhì)量逐步得以改觀。呈現(xiàn)惡化等級的區(qū)域僅占1.68%,呈斑塊狀分布于研究區(qū)中北部,城鄉(xiāng)擴建、礦產(chǎn)資源采伐、連年干旱少雨是生態(tài)環(huán)境質(zhì)量惡化的主要成因。

2000-2010年臨河區(qū)生態(tài)環(huán)境質(zhì)量改善狀況最為顯著,改善面積達880.46 km2,農(nóng)田、林木覆蓋區(qū)為主要改善地,17.50%的區(qū)域呈現(xiàn)跨級差改善,2006年臨河區(qū)下發(fā)《關于實施戶植千棵樹工程》,為全區(qū)新增400 km2的造林面積,區(qū)域生態(tài)環(huán)境效益得以顯著提升。呈現(xiàn)惡化等級的生態(tài)環(huán)境質(zhì)量區(qū)域占比4.57%,其中97.37%的區(qū)域惡化1個級差,南部黃河沿岸硬化改造工程的實施及中部、西北部城鎮(zhèn)化的迅速擴建是區(qū)域生態(tài)環(huán)境質(zhì)量惡化的主要成因。全區(qū)1/2以上的區(qū)域生態(tài)環(huán)境質(zhì)量呈穩(wěn)定等級。1989-2000、2000-2010年生態(tài)環(huán)境質(zhì)量改善區(qū)域與惡化區(qū)域面積比分別為4.33和8.29,表明臨河區(qū)生態(tài)環(huán)境質(zhì)量改善速度同等在提升。

2010-2021年臨河區(qū)生態(tài)環(huán)境質(zhì)量呈穩(wěn)定狀態(tài)的區(qū)域面積為1 283.12 km2;44.70%的區(qū)域表征為惡化狀況,其中跨級惡化涵蓋186.17 km2,僅有17.58 km2的區(qū)域呈現(xiàn)改善狀況,其中51.75%的區(qū)域表現(xiàn)為跨級改善。2010-2021年臨河區(qū)生態(tài)環(huán)境質(zhì)量總體呈穩(wěn)定-惡化演變趨勢,穩(wěn)定區(qū)域占據(jù)主導優(yōu)勢,可惡化面積增幅顯著,建成區(qū)擴張、礦產(chǎn)資源采伐和高標準良田的整改是致使區(qū)域生態(tài)環(huán)境質(zhì)量惡化的關鍵因素,生態(tài)環(huán)境質(zhì)量跨級改善,得益于政府對濕地資源的精準保護。

整體分析臨河區(qū)32 a間生態(tài)環(huán)境質(zhì)量演變進程,從變化幅度分析(表3),呈惡化等級的生態(tài)環(huán)境質(zhì)量區(qū)域為145.53 km2,其中99.96%的區(qū)域表征為惡化1個級差;生態(tài)環(huán)境質(zhì)量呈改善狀況的區(qū)域為4.19 km2,改善1個級差的區(qū)域占比99.63%;臨河區(qū)生態(tài)環(huán)境質(zhì)量整體展現(xiàn)穩(wěn)定發(fā)展局勢。從空間分布來看(圖5),研究區(qū)北部的防護林及商品林帶為主要的生態(tài)環(huán)境質(zhì)量改善區(qū),規(guī)范化的營林機制和對林草業(yè)嚴格保護制度的實施,是臨河區(qū)生態(tài)環(huán)境質(zhì)量得以改善的主導因素;呈現(xiàn)惡化狀況的生態(tài)環(huán)境質(zhì)量區(qū)域主要位于城鄉(xiāng)建設區(qū)和礦產(chǎn)采伐區(qū)周邊。警示臨河區(qū)應秉持生態(tài)與保護優(yōu)先的生態(tài)-經(jīng)濟發(fā)展理念。

總時段與1989-2000、2000-2010、2010-2021年各研究階段穩(wěn)定等級面積占比分別為93.58%、91.16%、57.71%、54.55%,1989-2000、2000-2010年生態(tài)環(huán)境質(zhì)量呈大面積小幅度改善態(tài)勢,而2010-2021年生態(tài)環(huán)境質(zhì)量呈小級差大范圍惡化態(tài)勢,生態(tài)環(huán)境質(zhì)量演變整體展現(xiàn)穩(wěn)定發(fā)展形勢。

2.4 水體演變

分析1989-2021年土地利用轉移矩陣及空間布局(表4、圖6),明晰研究區(qū)內(nèi)水體同各土地利用類型之間的結構特征轉變和方向轉移。總體而言,1989-2021年臨河區(qū)各類型、各區(qū)域水體演變差異顯著,水域面積縮減50.65 km2,僅有18.78%的前期水體保留至今。

表4 1989年-2021年臨河區(qū)土地利用轉移矩陣Table 4 Land use transfer matrix in the Linhe region,1989-2021 km2

1989-2021年水體演變?yōu)楦黝愋屯恋乩妹娣e76.26 km2,41.38%的水體轉變?yōu)榻ㄔO用地,占比最大。演變后各用地類型多數(shù)位于較差、差的生態(tài)環(huán)境質(zhì)量區(qū),結合圖6可以看出,城市、城鎮(zhèn)化的迅速擴建和黃河沿岸硬化改造工程的實施,使研究區(qū)中部、南部的坑塘、湖泊及黃河水域面積大幅縮減。

1989-2021年臨河區(qū)各地類演變?yōu)樗w的區(qū)域占演變總面積的1.87%,2.38%的建設用地演變?yōu)樗w,占比最小。演變至水體的各地類均源于生態(tài)環(huán)境質(zhì)量等級為中和較差區(qū)域,分析原因得出,城鄉(xiāng)鎮(zhèn)周邊小型觀賞湖泊的建造、水產(chǎn)養(yǎng)殖業(yè)的迅速發(fā)展及黃河國家濕地公園的建立,是水體擴充的主要源頭。

整體來看,32 a間臨河區(qū)水體演變類型、演變區(qū)域差異顯著,水體演變?yōu)楦鞯仡惣案鞯仡愌葑優(yōu)樗w的面積分別為76.26、25.61 km2,其中,轉化為水體的各地類均源于RSEI等級為中、較差的區(qū)域,經(jīng)水體演變形成的各地類均位于RSEI等級為差、較差區(qū)域。綜上分析得知,水體演變對生態(tài)環(huán)境質(zhì)量起負面效應。

3 結論

1989-2021年,研究區(qū)RSEI均值以2010年為轉折點呈“先增長后下降”趨勢,由0.295降至0.254,整體降幅13.90%。32 a間各研究階段生態(tài)環(huán)境質(zhì)量較差等級面積分別占研究區(qū)總面積的99.82%、90.97%、60.23%、93.66%,差等級區(qū)域總體增幅6.13%,中、良、優(yōu)等級面積占比變化均小于1%。

空間分布而言,中、良等級生態(tài)環(huán)境質(zhì)量區(qū)域主要布局于北部的林木覆蓋區(qū)和南部黃河沿岸,中部和西南地區(qū)各年間出現(xiàn)不同位置的零星分布;城鄉(xiāng)建設區(qū)、良田整改區(qū)、資源采伐地生態(tài)環(huán)境質(zhì)量等級呈現(xiàn)為差。

空間演變方向分析,RSEI各等級空間重心分布狀況在2000-2010年由南向北呈自低到高規(guī)律分布,與1989、2021年空間重心分布狀況截然相反;重心遷移軌跡分析,差、較差、良、優(yōu)等級總體向南遷移,中等級向北遷移;遷移距離分析,差、良、優(yōu)等級相比中、較差等級遷移距離更遠。

時間演變尺度分析,總時段與各階段穩(wěn)定等級區(qū)域面積占比分別為93.58%、91.16%、57.71%、54.55%,表明研究區(qū)生態(tài)環(huán)境質(zhì)量總體以較差等級為主,經(jīng)定量化、可視化分析得知惡化區(qū)域較改善區(qū)域面積占比6.06%,城鎮(zhèn)區(qū)擴建、工業(yè)園區(qū)開發(fā)、礦產(chǎn)資源采伐同生態(tài)環(huán)境質(zhì)量間呈負相關關系,因此,在推動新型“城鄉(xiāng)建設”實施模式中,應統(tǒng)籌協(xié)調(diào)土地資源開發(fā)同生態(tài)環(huán)境建設之間的和諧關系。

水體演變與生態(tài)環(huán)境質(zhì)量,1989-2021年研究區(qū)水體總面積銳減50.65 km2,各類型、各區(qū)域水體演變面積達101.87 km2,經(jīng)水體演變形成的各地類均位于RSEI等級為差、較差區(qū)域,演變?yōu)樗w的各地類均源于RSEI等級為中、較差區(qū)域,水體演變對生態(tài)環(huán)境質(zhì)量起負面影響。