崔囤月, 王世東, 張學軍

(1.河南理工大學測繪與國土信息工程學院,焦作 454000; 2.河南省地質(zhì)礦產(chǎn)勘查開發(fā)局第二地質(zhì)礦產(chǎn)調(diào)查院,鄭州 450001)

0 引言

雄安新區(qū)于2017年4月1日由國務院宣布成立,這一國家型戰(zhàn)略的提出,意味著雄安新區(qū)將成為繼深圳經(jīng)濟特區(qū)和浦東新區(qū)之后又一個具有全國意義的新區(qū)[1-2],因此,雄安新區(qū)的生態(tài)問題得到了廣泛關(guān)注。植被作為區(qū)域的重要生態(tài)屏障,對區(qū)域的生態(tài)安全和可持續(xù)發(fā)展起著不可替代的作用,而植被覆蓋信息的提取及植被覆蓋變化的驅(qū)動因素分析對揭示生態(tài)系統(tǒng)環(huán)境變化、植被恢復和重建布局具有重要意義[3]。對于大面積區(qū)域的植被覆蓋信息提取,需要具有長序性和區(qū)域性的遙感數(shù)據(jù)支撐。植被覆蓋度是遙感數(shù)據(jù)反演的一種參數(shù),是衡量地表植被狀況的直觀指標,也是評價區(qū)域生態(tài)環(huán)境優(yōu)劣的重要因子[4]。目前對植被覆蓋度的估算方法可歸納為經(jīng)驗模型法和植被指數(shù)轉(zhuǎn)換法,歸一化差值植被指數(shù)(normalized difference vegetation index,NDVI)是植被指數(shù)轉(zhuǎn)換法中普適性較好且應用較廣泛的一種。何全軍等[5]選取比值植被指數(shù)(ratio vegetation index,RVI)、NDVI和增強型植被指數(shù)(enhanced vegetation index,EVI)進行了相關(guān)植被指數(shù)的應用分析,結(jié)果表明NDVI和RVI結(jié)果較穩(wěn)定,EVI具有數(shù)據(jù)異常的現(xiàn)象。像元二分模型是以NDVI為遙感信息量,對植被覆蓋度進行估算的一種經(jīng)驗模型方法,削弱了大氣、土壤背景和植被類型的影響,普適性較好[6-7],因而被研究者們廣泛應用。王思等[7]基于MODIS-NDVI數(shù)據(jù),利用像元二分模型對廣東省2000—2020年的植被覆蓋信息進行了提取,并利用一元線性回歸模型對植被覆蓋度進行了時空格局變化趨勢的模擬; 晉成名等[8]利用像元二模型對陜北地區(qū)2000—2015年的植被覆蓋信息進行提取,分析了氣候和降水對植被覆蓋變化的影響; 蘇迎慶等[9]基于Landsat數(shù)據(jù)對汾河流域2000—2018年的植被覆蓋度進行二次反演,得到植被覆蓋信息,從空間和時間的角度分析了植被覆蓋的變化特征。

植被覆蓋變化受氣候、地形、人為等多重因素的影響[10-11],研究者對植被覆蓋變化驅(qū)動因素的探究應用最廣泛的方法是相關(guān)性分析; 而地理探測器模型作為探測空間分異性,以及揭示其背后驅(qū)動因子的新興統(tǒng)計學方法[12-14],其理論基礎(chǔ)與相關(guān)性分析相似,目前已被廣泛應用到驅(qū)動因子探測方面。彭文甫等[15]利用地理探測器模型,對影響四川省2000—2015年植被覆蓋變化的自然因子,進行風險、因子、生態(tài)和交互探測; 李丹等[16]為了解黑龍江省1980—2015年耕地時空變化的驅(qū)動因素,對其4個階段的耕地變化,利用地理探測器進行了驅(qū)動因子探測與交互作用探測; 祁鵬衛(wèi)等[17]利用地理探測器,對重慶市2000—2019年植被覆蓋變化進行了驅(qū)動因子探測與交互作用探測; 趙子娟等[18]基于MOD13Q1數(shù)據(jù),利用像元二分模型和相關(guān)性分析,對西遼河流域2000—2018年植被覆蓋變化情況及影響因素進行了研究。

與植被覆蓋變化相比,土地利用變化更為顯著。雄安新區(qū)土地利用變化之快受全國矚目。地處華北村落聚集地區(qū),成立之前對該地區(qū)土地利用變化的研究尤為薄弱; 此外,將土地利用變化與植被覆蓋變化相結(jié)合進行研究更為鮮有。因此,本文進行的雄安新區(qū)土地利用與植被覆蓋變化之間的響應研究,對其生態(tài)問題、土地利用結(jié)構(gòu)優(yōu)化及未來城市發(fā)展規(guī)劃等問題具有重要意義?；诖?本文以遙感數(shù)據(jù)與GIS技術(shù)為基礎(chǔ),利用圖譜融合法對土地利用變化信息進行提取,基于像元二分模型對植被覆蓋信息進行提取,利用地理探測器模型對植被覆蓋變化進行潛在驅(qū)動因素探測,并利用圖譜融合法對土地利用變化與植被覆蓋變化之間的響應進行探究,進而為雄安新區(qū)的可持續(xù)發(fā)展提供一定的科學依據(jù)和對策建議。

1 研究區(qū)概況與數(shù)據(jù)源

1.1 研究區(qū)概況

雄安新區(qū)于2017年為打造北京非首都功能疏解集中承載地而建立,由河北省容城縣、雄縣、安新縣及周邊部分區(qū)域組成,處于冀中平原中部,屬太行山麓平原向沖積平原的過渡帶[19-20],位于E115°～116°、N38°～39°之間(圖1)。平均海拔3～21 m,地勢較高的土壤多為褐土、干旱土,地勢低的土壤多為潮土、砂姜黑土、半水成土、水成土。區(qū)域面積為1 576.6 km2,屬于華北村落聚集的平原地區(qū),距離北京和天津均為105 km,有較好的地理優(yōu)勢[2]。雄安新區(qū)擁有河北省淡水面積最大的湖泊——白洋淀,面積為366 km2,是河北省重要濕地保護區(qū),具有較好的生態(tài)潛力。本文只以安新縣、容城縣、雄縣作為研究區(qū)域。

圖1 雄安新區(qū)地理位置

1.2 數(shù)據(jù)源及其處理

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

從地理空間數(shù)據(jù)云(https: //www.gscloud.cn)獲取雄安新區(qū)1991年、2001年、2011年、2017年和2021年30 m分辨率的遙感影像,其中前3 a為Landsat5 TM數(shù)據(jù)集,后2 a為Landsat8 OLI數(shù)據(jù)集,時間節(jié)點均選在對地類容易區(qū)分的9月份。然后對其進行預處理和土地利用信息提取,將土地利用類型劃分為建設(shè)用地、耕地、林草地、其他用地和水域5種類別,利用支持向量機(support vector machine,SVM)分類,綜合分類精度達到85.50%以上。

1.2.2 植被覆蓋數(shù)據(jù)

為保證數(shù)據(jù)的一致性,植被覆蓋度反演仍采用Landsat數(shù)據(jù)。從USGS獲取5個時期的生長季(4—10月)數(shù)據(jù),進行波段運算得到月尺度NDVI數(shù)據(jù); 采用平均值合成法得到年尺度NDVI數(shù)據(jù),進行波段運算獲得植被覆蓋度數(shù)據(jù); 并對年份接近的2期數(shù)據(jù)進行差值計算,得到1991—2021年4個階段的植被覆蓋變化趨勢; 參考前人的分類標準[21-22]和雄安新區(qū)植被覆蓋變化的實際情況,對植被覆蓋變化情況(ΔFVC)劃分為5級: [-1.00,-0.60)為嚴重退化、[-0.60,-0.10)為輕度退化、[-0.10,0.10)為保持穩(wěn)定、[0.10,0.60)為輕度改善、[0.60,1.00]為明顯改善,最終用于驅(qū)動因子分析。

1.2.3 社會及自然驅(qū)動因子數(shù)據(jù)

從中國科學院資源環(huán)境科學數(shù)據(jù)中心(https: //www.resdc.cn)獲取土壤類型(圖2(c))、土壤質(zhì)地(圖2(d))和5期1 km人口分布網(wǎng)格數(shù)據(jù)。從中國經(jīng)濟社會大數(shù)據(jù)平臺(https: //data.cnki.net)獲取GDP數(shù)值。從網(wǎng)站(https: //eogdata.mines.edu)獲取5期夜光遙感數(shù)據(jù),利用夜光數(shù)據(jù)與GDP數(shù)值間接獲得5期1 km GDP分布網(wǎng)格數(shù)據(jù)。從地理空間數(shù)據(jù)云獲取高程數(shù)據(jù)(圖2(a)),通過高程獲取坡度(圖2(b))。然后對數(shù)據(jù)進行裁剪、數(shù)據(jù)離散化、類別劃分、重采樣,并對GDP和人口數(shù)據(jù)進行差值計算獲取其階段性變化數(shù)據(jù)((圖3(a)—(d))和(圖3(e)—(h)),最終用于驅(qū)動因子探測。

(a) 高程 (b) 坡度 (c) 土壤類型 (d) 土壤質(zhì)地

(a) 1991—2001年人口密度變化 (b) 2001—2011年人口密度變化 (c) 2011—2017年人口密度變化 (d) 2017—2021年人口密度變化

2 研究方法

2.1 土地利用變化圖譜建立

為了解各土地利用類型的演化過程,借鑒已有的圖譜融合法[23-25],對前后兩期土地利用類型的圖譜單元代碼按照式(1)實現(xiàn)圖譜融合,計算公式為:

C=10A+B

,

(1)

式中:C為各階段之間土地利用轉(zhuǎn)移情況的柵格屬性圖譜;A為前一時期的土地利用圖譜柵格單元屬性值;B為后一時期的土地利用圖譜柵格單元屬性值。引入土地利用變化率,計算公式為:

,

(2)

式中:P為土地利用變化率;Sxy為初期x土地利用類型轉(zhuǎn)化為末期y土地利用類型的圖譜面積或單元,其中x≠y;n為土地利用類型的數(shù)量[25]。

2.2 植被覆蓋度估算

本文以NDVI為遙感信息量,利用像元二分模型進行植被覆蓋度估算,彌補了NDVI對高植被覆蓋和低植被覆蓋研究效果欠佳的缺點[26-27]。NDVI和像元二分模型計算公式為:

NDVI=(NIR-R)/(NIR+R)

,

(3)

式中:NDVI為歸一化植被覆蓋指數(shù);NIR和R分別為近紅外波段和紅外波段的反射率。

,

(4)

式中:FVC為植被覆蓋度,取值范圍[0,1];NDVIsoil為無植被覆蓋的像元的NDVI值;NDVIveg為全植被覆蓋的像元的NDVI值[28]。通常這2個參數(shù)通過實測數(shù)據(jù)的FVCmax和FVC經(jīng)過反算獲得,考慮本次實驗沒有實測數(shù)據(jù)和參數(shù)的不穩(wěn)定性,因此取0.05為置信度,以NDVI累積頻率95.00%和5.00%分別確定2個參數(shù)。

根據(jù)水利部2008年頒布的《土壤侵蝕分類分級標準》中植被覆蓋度分級標準和相關(guān)研究[29],將植被覆蓋度(FVC)劃分為5級: [0,0.30)(低覆蓋度)、[0.30,0.45)(中低覆蓋度)、[0.45,0.60)(中等覆蓋度)、[0.60,0.75)(中高覆蓋度)、[0.75,1](高覆蓋度)。

2.3 地理探測器模型

地理探測器模型是探測空間分異性,以及揭示其背后驅(qū)動力的一種新興統(tǒng)計學方法[10,30-31],本文使用地理探測器模型中的風險因子探測,研究植被覆蓋變化的潛在驅(qū)動因素。

風險因子探測用于探測屬性Y的空間分異性,即本文中的植被覆蓋度變化Y是否因為空間位置或環(huán)境不同而產(chǎn)生了差異; 以及探測本文選取的驅(qū)動因子X對植被覆蓋度變化Y的影響程度,用q值度量,表達式為[10]:

,

(5)

,

(6)

SST=Nσ2

,

(7)

式中:q為驅(qū)動因子對空間分異性的影響程度;L為變量Y或因子X的分層(strata),即分類或分區(qū);Nh和N分別為層h和全區(qū)的單元數(shù);σh2和σ2分別為層h和全區(qū)的Y值的方差;SSW和SST分別為層內(nèi)方差之和和全區(qū)總方差。q的值域為[0,1],值越大說明Y的空間分異性越明顯,即q值越大表示自變量X對屬性Y的解釋能越強,反之則越弱,q值表示X解釋了100×q%的Y[12,32]。

3 結(jié)果與分析

3.1 土地利用變化

基于5期土地利用圖譜(圖4),利用圖譜融合法獲得1991—2021年4個階段的土地利用演化圖譜(圖5),并提取各土地利用類型之間的轉(zhuǎn)化信息(表1、表2)。

表1 1991—2021年雄安新區(qū)土地利用數(shù)量結(jié)構(gòu)變化

表2 1991—2021年雄安新區(qū)土地利用轉(zhuǎn)化類型面積與占比

(a) 1991年 (b) 2001年 (c) 2011年

(a) 1991—2001年 (b) 2001—2011年

綜合分析如下:

1)總體而言,雄安新區(qū)土地利用結(jié)構(gòu)主要以耕地和建設(shè)用地為主,水域、未利用用地和林地次之,與華北村落聚集地區(qū)的土地利用結(jié)構(gòu)相吻合。

2)從土地利用的變化特征來看,1991—2021年建設(shè)用地始終處于增加趨勢,2011—2017年和2017—2021年容城縣中心區(qū)域建設(shè)用地的增長最為明顯,這與雄安新區(qū)的建設(shè)有較大關(guān)系,該區(qū)域為雄安新區(qū)建設(shè)的啟動區(qū),增長的面積主要來源于耕地和其他用地,流出的面積主要變成耕地。其他用地整體減少,主要是2017—2021年產(chǎn)生大幅度減少,流出的面積主要變成耕地。林草地整體增加,其中1991—2001年林草地處于減少趨勢,流出的面積主要變成耕地; 2001—2021年林草地處于增長趨勢,增長的面積主要來源于水域和其他用地,尤其是白洋淀水域附近的植被覆蓋面積大幅度增加,這離不開國家和雄安新區(qū)政府對白洋淀生態(tài)保護濕地的重視。水域整體減少,1991—2001年處于增長趨勢,增長的面積主要來源于耕地; 在2001—2017年大幅度減少,2017年面積達到最低,流出的面積主要變成耕地; 2017年11月以緩解河北省中南部農(nóng)業(yè)缺水和為白洋淀生態(tài)補水為目的的“引黃入冀補淀工程”試通水,2019年6月雄安新區(qū)政府頒布重點建設(shè)孝義河的政策,因此2017—2021年水域面積顯著增長,增長的面積主要來源于水域附近的耕地和其他用地。耕地無明顯變化。

3)從土地利用的轉(zhuǎn)化特征來看,1991—2001年安新縣中南部的其他用地→耕地轉(zhuǎn)化最為明顯,雄縣西南部和溢洪道附近的耕地→其他用地轉(zhuǎn)化較為明顯,建設(shè)用地→耕地轉(zhuǎn)化主要發(fā)生在安新縣西南部,耕地→建設(shè)用地轉(zhuǎn)化呈離散型分布,無聚集性轉(zhuǎn)化。2001—2011年雄縣東部發(fā)生耕地向建設(shè)用地和其他用地的聚集性轉(zhuǎn)化,這與圖4中2011年雄縣東部出現(xiàn)建設(shè)用地和其他用地大面積增長特征相吻合; 雄縣南部和溢洪道附近的其他用地→耕地轉(zhuǎn)化較為明顯,建設(shè)用地→耕地轉(zhuǎn)化呈離散型分布。2011—2017年容城縣和雄縣的耕地→建設(shè)用地轉(zhuǎn)化最為明顯,其次是雄縣東部和溢洪道附近的其他用地→耕地和建設(shè)用地→耕地轉(zhuǎn)化,以及呈離散型分布的耕地→其他用地轉(zhuǎn)化。2017—2021年容城縣中部發(fā)生大面積的耕地→建設(shè)用地轉(zhuǎn)化,該區(qū)域為啟動區(qū)的容西和容東安置區(qū),也是自雄安新區(qū)成立以來變化最顯著的區(qū)域; 雄縣部分區(qū)域也發(fā)生耕地→建設(shè)用地轉(zhuǎn)化,其他用地和耕地之間的轉(zhuǎn)化呈離散型分布,其余轉(zhuǎn)化類型均不明顯。

3.2 植被覆蓋時空變化

3.2.1 植被覆蓋空間分布特征

1991—2021年雄安新區(qū)整體植被覆蓋情況良好,中高及以上植被覆蓋度的區(qū)域面積占雄安新區(qū)總面積的50%以上(表3)。從平均植被覆蓋度各等級的空間分布來看,安新縣中高及以上植被覆蓋度區(qū)域面積占比最高,雄縣最低; 中等及以下植被覆蓋度的區(qū)域總面積在各區(qū)域的占比由高到低順序為: 雄縣、容城縣、安新縣(圖6)。因此,安新縣具有較好的生態(tài)環(huán)境,這與安新縣的地理環(huán)境有關(guān),因為白洋淀生態(tài)保護濕地坐落于安新縣東南部; 生態(tài)環(huán)境較差的是雄縣,這與雄縣的經(jīng)濟發(fā)展和人類的頻繁活動有關(guān),雄縣是雄安新區(qū)經(jīng)濟實力最強的地區(qū),說明該地區(qū)的經(jīng)濟發(fā)展極大影響到了植被的生長,因此雄安新區(qū)政府要對該雄縣的生態(tài)環(huán)境增加重視程度,實現(xiàn)經(jīng)濟與生態(tài)同進步。

表3 雄安新區(qū)年均植被覆蓋度

3.2.2 植被覆蓋階段性變化趨勢

為深刻了解雄安新區(qū)植被覆蓋變化趨勢,對1991—2021年5期植被覆蓋度數(shù)據(jù)進行差值計算,得到1991—2021年植被覆蓋階段性變化趨勢(表4和圖7),并進行了如下分析:

(a) 1991—2001年 (b) 2001—2011年

1)1991—2021年雄安新區(qū)植被整體經(jīng)歷了“退化→改善”的過程。1991—2001年,植被整體退化,退化區(qū)域總面積大于改善區(qū)域總面積,退化區(qū)域主要集中于溢洪道和雄縣中東部,占雄安新區(qū)總面積的39.15%; 改善區(qū)域總面積占比達到29.91%,改善區(qū)域主要集中于安新縣西南部和白洋淀水域附近。保持穩(wěn)定的區(qū)域在3個縣區(qū)均有離散型分布。

2)與上一階段相比,2001—2011年雄安新區(qū)植被有明顯改善,輕度改善和明顯改善的區(qū)域面積較上一時期有所提高,總占比34.80%,改善區(qū)域位置發(fā)生較大變化,集中于雄縣西部、容城縣東部和安新縣北部。退化區(qū)域的總面積與上一時期有所降低,總占比為31.27%; 其中,嚴重退化的區(qū)域面積占比最低,集中區(qū)域位于安新縣白洋淀水域附近。保持穩(wěn)定的區(qū)域面積占比有小幅度擴張,在3個縣區(qū)均呈離散型分布。

3)2011—2017年相比于2001—2011年,植被整體仍有較大改善,改善區(qū)域面積總占比分47.55%,退化區(qū)域面積總占比26.59%。其中,明顯改善的區(qū)域面積大幅度增加,同比上一階段,增長3倍之多,區(qū)域主要集中于雄縣東部和安新縣西南部; 退化的區(qū)域和保持穩(wěn)定的區(qū)域在3個縣區(qū)呈離散型分布,其中白洋淀水域附近的退化較明顯。

4)2017—2021年是雄安新區(qū)成立4周年,植被仍然呈改善趨勢,改善區(qū)域總面積遠大于退化區(qū)域總面積,退化面積總占比達到4個階段中最低。嚴重退化區(qū)域發(fā)生位置變化,主要集中于啟動區(qū)的容西和容東安置區(qū),說明雄安新區(qū)的建設(shè)是該區(qū)域植被退化的重要原因,改善區(qū)域集中于容城縣西部、雄縣北部和白洋淀水域,這和2019年6月雄安新區(qū)政府發(fā)布的白洋淀生態(tài)屏障重點保護政策有關(guān)。

3.3 植被覆蓋變化驅(qū)動因素

3.3.1 植被覆蓋變化潛在驅(qū)動因素

植被覆蓋變化受多重因素的影響,本文選取自然方面的土壤質(zhì)地、土壤類型、高程、坡度、社會方面的人口密度變化、GDP 變化等6個驅(qū)動因子,分別作為自變量X1,X2,X3,X4,X5,X6。以植被覆蓋階段性變化趨勢作為因變量Y,利用地理探測器模型對雄安新區(qū)1991—2021年植被覆蓋變化趨勢進行驅(qū)動因子分析(表5)。

表5 驅(qū)動因子作用強度q值變化趨勢

分析結(jié)果表明: 植被覆蓋變化存在空間分異性,1991—2001年,土壤質(zhì)地對植被覆蓋變化影響程度最強; 土壤類型、人口密度變化和GDP變化對植被覆蓋變化也具有較強的影響程度。2001—2011年,人口密度變化仍是驅(qū)動因子中影響程度最強的,且影響程度較上一時期有所提高; 影響程度較強的仍是土壤類型、土壤質(zhì)地和GDP變化,且土壤類型和GDP的影響程度較上一時期有所提高,土壤質(zhì)地的影響程度較上一時期有所降低。2011—2017年,GDP變化是驅(qū)動因子中影響程度最強的,且影響程度較上一時期有所提高,達到30 a間q值最高峰; 影響程度較強的是土壤質(zhì)地、土壤類型和人口密度變化,但人口密度變化影響程度較上一時期有所降低。2017—2021年,土壤質(zhì)地是驅(qū)動因子影響程度最強的; 影響程度較強的驅(qū)動因子與上一時期相同,除GDP變化外,影響都有所提高。高程和坡度對植被覆蓋變化影響程度都較低,說明高程和坡度對雄安新區(qū)植被覆蓋變化的影響較小,這與雄安新區(qū)的地理環(huán)境有關(guān)。

3.3.2 土地利用變化對植被覆蓋變化的影響

除以上潛在驅(qū)動因素外,本文也探究了土地利用變化對植被覆蓋變化的影響。土地利用變化圖譜中共有20種轉(zhuǎn)化類型,大部分是發(fā)生少量圖元的轉(zhuǎn)化,因此本文將4個階段的重要土地利用轉(zhuǎn)化類型圖譜(圖5)與植被覆蓋變化趨勢圖譜(圖7)進行基于公式1的圖譜融合,并計算各轉(zhuǎn)化類型對植被變化的貢獻率(圖8)?？傮w來看,1991—2001年30 a間,重要土地利用轉(zhuǎn)化類型對植被退化的總貢獻率最高達到45.19%,對植被改善的總貢獻率最高達到31.64%,表明土地利用變化對植被覆蓋變化存在較大影響。

(a) 1991—2001年 (b) 2001—2011年

統(tǒng)計結(jié)果表明:

1)1991—2001年,重要土地利用轉(zhuǎn)化類型對植被退化和改善的總貢獻率分別達到30.27%和31.64%。植被嚴重退化區(qū)域主要集中于雄縣中部及溢洪道附近,與產(chǎn)生耕地→其他用地的區(qū)域分布相似,表明耕地→建設(shè)用地是該時段該區(qū)域植被嚴重退化的重要原因。安新縣產(chǎn)生其他用地→耕地和建設(shè)用地→耕地的區(qū)域與安新縣出現(xiàn)的植被明顯改善區(qū)域分布相似,說明該2種轉(zhuǎn)化類型是該階段植被改善的重要原因。

2)2001—2011年,重要土地利用轉(zhuǎn)化類型對植被退化和改善的總貢獻率分別達到28.68%和20.40%,其他用地→耕地轉(zhuǎn)化對植被明顯改善的貢獻率最大,說明其他用地→耕地轉(zhuǎn)化是該時段植被改善的重要原因; 水域→耕地轉(zhuǎn)化和耕地→建設(shè)用地轉(zhuǎn)化對植被嚴重退化貢獻率較大,說明該2種轉(zhuǎn)化類型是該階段植被退化的重要原因; 而且發(fā)現(xiàn)水域減少也會造成植被退化現(xiàn)象。

3)2011—2017年,土地利用變化對植被退化和改善的影響最小,總貢獻率分別為15.70%和27.50%。雄縣出現(xiàn)的其他用地→耕地轉(zhuǎn)化區(qū)域與植被明顯改善的區(qū)域相似,且其他用地→耕地對植被明顯改善的貢獻率較大,說明該轉(zhuǎn)化類型是該階段植被改善的重要原因。

4)2017—2021年,土地利用變化對植被退化的影響較大,總貢獻率達到45.18%。容城縣中部的植被嚴重退化區(qū)域與耕地→建設(shè)用地轉(zhuǎn)化的區(qū)域相似,該轉(zhuǎn)化類型對植被嚴重退化貢獻率達到最大,說

明耕地向耕地→建設(shè)用地轉(zhuǎn)化是該階段植被嚴重退化的重要原因。

根據(jù)以上分析可以看出: 耕地向建設(shè)用地和其他用地的轉(zhuǎn)化是植被退化的重要原因,建設(shè)用地向其他用地向耕地的轉(zhuǎn)化是植被改善的重要原因,其余土地利用轉(zhuǎn)化類型對植被變化均是較小的影響。

4 問題與討論

本文以遙感數(shù)據(jù)和GIS技術(shù)為基礎(chǔ),研究了雄安新區(qū)1991—2021年近30 a土地利用與植被覆蓋變化特征,利用地理探測器探究了植被覆蓋變化的潛在驅(qū)動因素,并探究了土地利用變化與植被覆蓋變化之間的響應。對主要問題討論如下:

1)關(guān)于土地利用類型之間的相互轉(zhuǎn)換。城市的擴張和新區(qū)的建設(shè)使得3縣建設(shè)用地不斷增長,水域減少,林草地增加較少,應注重規(guī)避其造成的負面影響,例如土地城鎮(zhèn)化后形成的用地浪費、白洋淀生態(tài)屏障的損傷以及林草地較少造成的水土流失問題,尤其是容城縣中部地區(qū)和安新縣白洋淀水域附近,應注重城市的未來整體規(guī)劃,建立土地管理體制,時刻監(jiān)測其土地利用變化情況,實現(xiàn)城市發(fā)展與生態(tài)環(huán)境的同步提升,最終將雄安新區(qū)建設(shè)成為一個集綠色、生態(tài)、韌性、美麗于一體的經(jīng)濟特區(qū)。

2)植被覆蓋情況的退化與改善。30 a間,雄縣是3縣之中GDP最高,經(jīng)濟發(fā)展最快的一個,但發(fā)展經(jīng)濟的同時,應注重生態(tài)環(huán)境質(zhì)量的同步提升。白洋淀水域附近的植被覆蓋變化始終不穩(wěn)定,應建立該區(qū)域的植被覆蓋監(jiān)管體制,防止水土流失,努力將安新縣打造成為雄安新區(qū)的“后花園”。

3)不同驅(qū)動因素對植被覆蓋的影響。容城縣中部和雄縣區(qū)域的土地利用變化對植被退化的貢獻最為明顯,應加強該區(qū)域的土地利用變化關(guān)注度,防止植被退化現(xiàn)象的再次出現(xiàn)和嚴重化,安新縣西南部和雄縣區(qū)域的土地利用變化對植被改善的貢獻較為明顯,應對其加以借鑒。

5 結(jié)論

1)1991—2021年雄安新區(qū)建設(shè)用地整體增長,轉(zhuǎn)入面積主要來源于耕地和其他用地; 其他用地面積整體減少,流出面積主要變成耕地; 林草地面積同比增長最多,轉(zhuǎn)入面積主要來源于水域和其他用地; 水域面積整體減少,流出面積主要變成耕地和其他用地; 耕地面積無明顯變化,主要是與建設(shè)用地和其他用地之間的相互轉(zhuǎn)化。

2)1991—2021年雄安新區(qū)整體植被覆蓋情況良好,其中安新縣植被覆蓋情況最好,生態(tài)環(huán)境較優(yōu)越,容城縣次之,雄縣最差。植被整體呈現(xiàn)退化→改善的趨勢,1991—2001年植被整體退化,2001—2021年植被整體均呈現(xiàn)改善趨勢。1991—2001年雄縣出現(xiàn)明顯退化,白洋淀水域附近出現(xiàn)明顯改善; 2001—2011年雄縣西部、容城縣東部和安新縣北部出現(xiàn)明顯改善,白洋淀水域附近出現(xiàn)明顯退化; 2011—2017年,雄縣東部和安新縣西南部出現(xiàn)大面積改善區(qū)域,退化區(qū)域在3個縣區(qū)均有離散分布,其中白洋淀水域附近較明顯; 2017—2021年白洋淀水域附近、容城縣西部和雄縣北部出現(xiàn)改善,容城縣中部地區(qū)出現(xiàn)退化。

3)土壤類型、人口密度變化和GDP變化對雄安新區(qū)植被覆蓋變化均有較強的影響程度,其中1991—2001年土壤質(zhì)地影響最大,2001—2011年人口密度變化影響最大,2011—2017年GDP變化影響最大,2017—2021年土壤質(zhì)地影響最大,因地處平原,高程和坡度差異較小,所以對其影響也較小。除此以外,建設(shè)用地和其他用地向耕地的轉(zhuǎn)化是植被改善的重要原因,耕地向建設(shè)用地和其他用地的轉(zhuǎn)化是植被退化的重要原因。

4)地理探測器模型能有效地分析植被覆蓋變化的驅(qū)動機制,對于大區(qū)域的植被覆蓋變化驅(qū)動研究具有借鑒意義,在生態(tài)環(huán)境監(jiān)測、城鄉(xiāng)規(guī)劃和可持續(xù)發(fā)展方面具有良好的應用潛力。由于本文研究區(qū)域和數(shù)據(jù)的局限性,驅(qū)動因子種類選取較少,仍具有改進空間,但在一定程度上反映了人為因素和自然因素對雄安新區(qū)植被覆蓋變化的影響。