模擬多情景下白洋淀流域土地利用變化及生態系統服務價值的空間響應

高 星,楊劉婉青,李晨曦,宋昭穎,王 杰

1 河北經貿大學公共管理學院, 石家莊 050061

2 河北經貿大學河北省社會經濟社會發展地理信息大數據平臺, 石家莊 050061

3 西安建筑科技大學公共管理學院, 西安 710311

生態系統服務指通過生態系統的結構、過程和功能直接或間接得到的生命支持產品和服務[1- 2],其形成、供給和分配受到土地利用的深刻影響。隨著人口增長、經濟發展和城市的急速擴張,土地利用形式多樣和強度加強,自然生態系統逐漸向半自然半人工生態系統和人工生態系統方向轉變,由此帶來的物質和能量流的改變影響生態系統提供服務的能力。因此,不適宜的土地利用將危及區域生態安全,甚至限制城市的可持續發展[3]。雄安新區建設作為疏解北京非首都功能和推進京津冀協同發展的重要國家戰略,其過程中必然伴隨著高強度的人為建設活動,與當地自然生態環境不相適宜的土地利用方式勢必威脅區域生態系統服務的供需平衡,甚至導致生態環境問題[4]。同時,隨著雄安新區人口的聚集和城市規模的擴張,對糧食供給、水資源供給、空氣凈化、水源涵養和文化休閑等生態系統服務的需求將不斷提高[5]。而白洋淀流域作為雄安新區的生態腹地,既是由自然地理和經濟社會發展組成的復合性區域,又是實施生態建設和環境政策的整體空間,生態戰略地位顯要[6]。在此背景下,模擬白洋淀流域未來土地利用變化趨勢,探尋不同情境下生態系統服務價值的響應有利于構建流域國土空間開發保護新格局,為其提供技術支撐和決策參考[7]。

目前,國內外學者對于土地利用變化對生態系統服務價值的影響已積累了相對豐富的研究成果。Costanza等[8]在1997年制訂了全球生態系統服務當量因子表,為相關研究奠定了基礎。國內學者謝高地等[9]在Costanza研究的基礎上制訂了《中國陸地生態系統單位面積生態服務價值當量表》,被國內學者廣泛使用于生態系統服務價值評估研究中。隨著研究的逐步深入,國內外學者對生態系統服務的研究主要從定義與分類[8]、價值評估[10- 12]、以及生態系統服務與社會發展之間的關系[13- 15]等三個方面進行了鉆研與探索。具體而言,從研究內容上看,主要集中在生態系統服務的理論基礎[16]、驅動力[17- 18]、時空變化[19]及影響因素[20]等方面；從研究尺度上看,全國、區域、省域、市域、縣域等尺度均有涉及[21-25]。從研究對象上看,研究區域多集中湖泊[26]、森林[27]、濕地[28]等生態系統服務價值高值區域；研究方法則以統計分析[29]和空間分析[30-31]為主,并越來越重視在土地利用規劃研究與實踐中的作用[32]。綜合來看,上述研究主要是對當前或過去時段因土地利用變化而改變的生態系統服務價值進行測度,隨著研究的深入,學者們嘗試著對未來土地利用變化對區域生態系統服務價值的影響進行分析,以便于制定可持續的土地利用措施,如Egoh等[33]采用Marxan模型模擬不同情境下南非草地生態系統服務狀況；Estoque等[34]研究Baguio市未來森林保護情景和自然發展情景下土地利用變化對生態系統服務價值的影響；Kubiszewski等[35]估算了全球尺度下市場導向情景、政策改良情景等4種情景下生態系統服務價值,結果顯示到2050年4種情景生態系統服務價值每年可能相差81萬億美元；杜勇等[36]通過Logistic-CA-Markov模型模擬福建省泉州市在2030年自然情景、規劃情景、保護情境下,定量評估了保水服務、固碳服務、保土服務和食物供給服務；歐陽曉等[37]以快速城鎮化下的長株潭城市群為例,模擬了3種情景下LUCC對ESV的影響。目前已有的多情景下生態系統服務價值的模擬測度,大部分是基于現有發展趨勢根據側重點設定不同情景進行分析,缺乏對由國家戰略引發的土地利用類型劇烈變化區域生態系統服務價值空間差異性的關注。

雄安新區地處白洋淀流域范圍內,面臨著人口增長、社會經濟發展和生態保護的多重壓力,同時新區的開發建設勢必對白洋淀流域生態系統服務價值產生一定的影響。鑒于生態系統服務價值計算的便利性,本文以白洋淀流域為研究區,采用當量因子法估算其1990、2000、2010和2020等4期生態系統服務價值,并借助CA-Markov模型,模擬不同情境土地利用格局,并測算生態系統服務價值變化情況。在此基礎上利用網格工具可視化表達不同偏好下生態系統服務價值的空間分布狀況以及集聚程度,以期對白洋淀流域的生態系統服務功能變化分析以及對研究區未來土地生態安全建設提供理論參考。

1 研究區概況與數據來源

1.1 研究區概況

本研究中的白洋淀流域是指白洋淀上下游涉及的河北省相關縣級行政區域,涉及石家莊市、保定市、張家口市、滄州市、廊坊市以及衡水市的46個縣區；白洋淀流域上下游包括入淀出淀的潴龍河、孝義河、唐河、府河、漕河、瀑河、萍河、白溝引河(含南拒馬河與白溝河)、趙王新河等干支流流經的區域。白洋淀流域位于大清河流域范圍以內,總面積約3.4萬km2。該地區地貌類型多樣,以耕地、林地和建設用地為主,地勢西高東低,相對高度達2815m(圖1),流域處于溫帶季風區,多年平均降水量600mm,年平均氣溫在10—20℃之間,降水量年際變化大、年內分配不均,降水多以暴雨形式集中于7—8月,故經常出現連年干旱或先旱后澇、旱澇交替等情況。近年來,白洋淀流域各城市無序擴張、產業布局各自為政、大尺度生態空間管理理念缺乏,導致生態用地流失、水生態失衡等問題,由此帶來生態系統服務下降、生態風險加劇、生態安全受到威脅,嚴重影響了流域的可持續發展,與高質量發展要求存在一定差距。

圖1 研究區地理位置Fig.1 Location of Baiyangdian basin

1.2 數據來源及預處理

本研究所使用的數據包括遙感影像數據、土地利用數據、糧食數據、DEM高程數據、坡度數據、社會經濟數據、道路數據和氣象數據。其中：(1)遙感影像從地理空間數據云平臺Landsat4- 5/TM和Landsat- 8/OLI衛星下載；(2)研究區1990年、2000年、2010年、2020年等4期土地利用均解譯自遙感影像,根據影像精度并結合研究區實際情況,并參照前人分類標準[38],將土地利用類型劃分為林地、草地、耕地、建設用地、水域和未利用地六類,分辨率為30m。使用ENVI 5.3.1軟件,對遙感圖像進行監督分類和人工目視解譯,經檢驗三期圖像精度均達到85%以上,可用于研究使用(圖2)；(3)糧食數據主要為小麥、玉米、大豆等糧食作物產值、播種面積、產量主要來源為石家莊市、張家口市、保定市、滄州市和衡水市統計年鑒以及河北農村統計年鑒；(4)白洋淀流域年平均降水量和平均氣溫等氣象數據收集自中國氣象數據網；(5)研究區DEM數據來源于地理空間數據云,坡度數據由DEM高程數據提取得到；(6)社會經濟數據來自相關縣域統計年鑒；(7)公路、鐵路數據來源于資源環境數據云平臺。

圖2 白洋淀流域1990、2000、2010、2020年4期土地利用數據Fig.2 Land use data of Baiyangdian basin in the four phases of 1990, 2000, 2010 and 2020

2 研究方法

2.1 生態系統服務價值估算

本文基于謝高地等[1]建立的“中國陸地生態系統單位面積服務價值當量表”(表1)(根據研究區實際地類情況以及計算的便利性,除建設用地外五種地類分別采用了當量表中的旱地、闊葉、灌草叢、濕地、裸地和水系的單位面積服務價值當量),估算白洋淀流域各縣域1990—2020年4期生態系統服務價值量,初步分析多年間各項生態系統服務價值量的變化。通過公式(1)和公式(2)計算研究區的生態系統服務價值。

表1 單位面積生態系統服務價值當量Table 1 Ecosystem service equivalent value per unit area

(1)

(2)

式中,ESV為生態系統服務價值,單位為元/年；i為土地利用類型；j為生態系統服務類型；Ai為第i類土地利用類型的面積,單位為hm2；VCi為第i類土地利用類型單位面積生態系統服務價值,單位為元hm-2a-1；ECj為某類土地利用類型第j項生態系統服務價值當量；k為生態系統服務類型數量；Ea為1單位生態系統服務的經濟價值,單位為元hm-2a-1。

根據“現有單位農田提供的生產服務經濟價值,是無人力投入的自然生態系統的經濟價值的七倍”的規則[10],計算出白洋淀流域單位農田自然糧食產量的年均經濟價值。為更加精確的確定研究區單位生態系統服務的經濟價值,結合研究區實地情況,運用公式(3)對該數值進行修正[39]。鑒于統計數據的可獲取性,不考慮價格波動因素統一選取2018年研究區各縣區糧食作物(小麥、玉米、大豆)的經濟價值進行修正,并計算其平均值,結果表明研究區每單位生態系統服務的經濟價值為1801.48元hm-2a-1。

(3)

式中,Ea為1單位生態系統服務的經濟價值,單位為元hm-2a-1；i為糧食作物種類；mi為第i種糧食作物在研究區內的平均價格,單位為元/kg；pi為第i種糧食作物的單產,單位為kg/hm2；qi為第i種糧食作物的種植面積,單位為hm2；M為糧食作物的總種植面積,單位為hm2。

2.2 土地利用格局情境模擬

2.2.1CA-Markov

馬爾科夫(Markov)模型[40-41]利用概率論的相關知識研究事物由一狀態轉換為另一狀態的概率,根據研究得出的規律,對未來事物的狀態進行預測,其重點強調的是對未來數量上的預測。其公式為：

S(t+1)=St×Pij

(4)

式中,St為t時刻土地利用系統的狀態；S(t+1)為t+1時刻土地利用系統的狀態；Pij為土地利用轉移矩陣。

元胞自動機(CA)模型是一種時間、空間、狀態均離散,通過不同時空特征元胞的相互作用,模擬復雜系統時空演變的動力學模型[42- 43]。

St+1=f(St,Ν)

(5)

式中,S為元胞有限、離散的狀態集合；N為元胞的鄰域；t和t+1表示不同時刻；f為局部空間元胞轉換規則。

傳統的馬爾科夫(Markov)模型很難預測土地利用在空間格局上的變化,而元胞自動機(CA)模型則可以有效的模擬這種變化,然而CA模型著眼于元胞的局部相互作用,存在很大的局限性,從這一角度來說,CA-Markov模型將Markov模型和CA模型結合在一起,彌補各自模型存在的缺陷,實現對土地利用類型數量和空間位置的雙重預測[44- 45]。

2.2.2土地利用格局情景模擬標準設定與過程

通過IDRISI軟件中的Markov模塊和CA-Markov模塊對土地利用類型的數量及空間變化進行預測,基本步驟如下：

(1)創建 Markov 轉移概率矩陣。將時間周期設置為10a,比例誤差設置為0.15,利用 Markov 模型得到研究區2010—2020年的年土地利用轉移面積矩陣和轉移概率矩陣。

(2)制作適宜性圖集。利用IDRISI軟件中的“Decision Wizard”模塊制作適宜性圖集。借鑒前人的研究結果[39- 44],考慮到數據的可獲得性以及具體模擬效果,本文選取DEM高程、坡度、年平均降水、年平均氣溫、距公路距離和距鐵路距離等6個對土地利用變化較大的驅動影響因子來制作適宜性圖集。

(3)模擬2020年研究區土地利用圖像。以2010年為基期圖像,選擇步驟(1)建立的研究區2010—2020年的年土地利用轉移面積矩陣,以及步驟(2)中創建的適宜性圖集,確定迭代次數為10,構造CA濾波器(采用標準5×5的濾波器)。模擬運行出預測的2020年研究區土地利用圖像(圖3)。

(4)精度驗證。利用IDRISI軟件中的crosstab模塊,將2020年的預測圖與2020年實際土地利用類型進行對比,計算Kappa值。

(5)多情景模擬預測。以2020年為基期圖像,分別模擬出3種情景下2035年研究區土地利用圖像。

2.2.3情景設置

情景分析旨在描述和分析事物發展的各種可能性,對比不同發展路徑下的狀態,從而為決策提供依據。白洋淀流域社會經濟、自然環境均存在不確定性。基于流域的不同發展目標和可能出現的干擾情景,本研究設定了3種土地利用變化模擬的情景：

(1)趨勢延續情景(情景1)。趨勢延續情景是流域基于2010—2020年土地利用變化規律,按照目前的城市化發展模式和現有的土地利用轉移速率,在轉換規則中不對任何地類轉換進行限制的一種情景設定,是流域土地利用變化模擬考慮其他發展情景分析的基礎。

(2)均衡發展情景(情景2)。該情景充分考慮流域中各城市的發展需求,實現流域空間均衡發展。其中,雄安新區的土地利用結構按照《河北雄安新區規劃綱要》、《河北雄安新區起步區控制性規劃》等規劃文本及圖集中對未來各地類在不同規劃方向上的數量要求進行確定,流域其他地區的土地利用結構則與情景1一致。

(3)統籌發展情景(情景3)。該情景下,將整個流域城市發展納入統籌考慮,以雄安新區建設、生態保護作為目標。因此,參考了有關修改Markov轉移矩陣的相關研究成果[46],其中,雄安新區土地利用結構與均衡發展情景一致,流域除雄安新區部分為整體實現生態用地保護目標,嚴格限制生態用地相關的林地、草地和水域向建設用地轉換,以保障整個流域的土地利用結構總體穩定。具體參數設置為將林地轉化為建設用地的速率減緩63%,轉移速率由原來的4.45%調整為2.80%；將草地轉化為建設用地的速率減緩50%,轉移速率由原來的10.92%調整為5.47%；水域由于基數過小,未進行限制；耕地、未利用地變化速率未進行限制。

2.3 ESV熱點分析

熱點分析(Getis-Ord Gi*)是用來識別冷點區和熱點區空間分布的常用工具,在本文中熱點分析可以用來反映白洋淀流域生態系統服務價值是否具有高值集聚(熱點)和低值集聚(冷點)的現象,以及確定高值區和低值區在空間上發生聚類的位置。其公式為：

(6)

(7)

(8)

式中,n是空間網格單元數量；xi和xj分別為單元i和單元j的觀測值；wij是基于空間k鄰接關系建立的空間權重陣。若一定范圍內ESV較高時,則成為具有顯著統計學意義的熱點,稱為ESV熱點區,表示生態系統服務價值在該區較大；若一定范圍內ESV較低時,則成為具有顯著統計學意義的冷點,稱為ESV冷點區,表示生態系統服務價值在該區較小。

3 結果分析

3.1 生態系統服務價值估算結果及變化特征

3.1.11990—2020年研究區土地利用變化總體特征

1990—2020年研究區利用變化總體特征呈現耕地面積大幅度減少,建設用地面積急速增加的趨勢。30年間耕地在研究區范圍內土地利用結構中占比最大,平均占比48.12%,總體來看耕地面積銳減,減少122776.99hm2,耕地面積減少率為7.21%；林地在研究區范圍內土地利用結構中面積占比僅次于耕地,且1990—2020年間林地面積變化不大,呈波動增長趨勢,從1193822.47hm2增加到1207925.59hm2,增加率為1.18%；與之相比,草地面積在30年間呈波動下降的趨勢,2010年達到峰值212026.37hm2,面積減少率為15.9%；隨著改革開放和京津冀協同發展等政策的推進,建設用地擴張速度也在不斷加快。1990—2020年建設用地面積從253926.75hm2增加到421296.49hm2,共增加了167369.74hm2,面積增長率高達65.91%；1990—2020年研究區水域面積呈現了先急速下降后小幅度回升的變化趨勢,總體來說,水域面積明顯減少,減少率為35.76%；未利用地占比較小,從研究區整體來看,面積變化不明顯(表2)。

表2 1990—2020年白洋淀流域各類用地面積及占比/hm2Table 2 Changes of different land use types in Baiyangdian basin from 1990 to 2020

研究區在空間上土地利用總體特征為,西部以林地、草地為主,東部多為耕地和建設用地交錯分布,小部分為水域。這主要由于研究區地勢西高東低,西部以山地為主,東部為華北平原以及白洋淀流域。從4期研究區土地類型圖可以看出,城市的擴張是造成土地利用類型發生變化的關鍵因素。在快速的城鎮化過程中,耕地被大量侵占,為堅持生態文明建設和可持續發展,山林區域的耕地還林還草現象明顯,導致耕地面積進一步銳減,白洋淀流域也在城市的發展下逐步縮小,生態安全情況不容樂觀(圖2)。

3.1.21990—2020年研究區生態系統服務價值變化特征

土地利用變化會引起生態系統服務價值的變化。1990—2020年研究區生態系統服務價值分別為852.59億元、788.26億元、796.95億元和784.68億元,總體呈波動下降的趨勢(表3)。其中,耕地的生態系統服務價值變化較多,30年間減少9.03億元,變化率為7.21%；林地對研究區的生態系統服務價值貢獻最大,四個時段分別為502.30億元、502.04億元、510.22億元和508.23億元,可以看出1990—2000年變化不明顯,2000—2020年價值量總體呈波動上升的趨勢,這是由于2000年以后林地面積小幅度增加；草地生態系統服務價值變化趨勢與其面積變化趨勢相同,1990—2020年波動下降,30年間ESV減少12.53億元；建設用地的生態系統服務價值為0；水域的生態系統服務價值僅次于林地,聯系表2可知,它憑借1%左右的面積占比提供了整個研究區12%左右的生態系統服務價值,然而由于城市快速發展的需要,建設用地的急速擴張使得一部分水域被侵占,導致水域的ESV下降為93.92億元,下降率為36.76%；未利用地由于其面積基數小,以及單位生態系統服務價值不高,對整個研究區的生態系統服務價值的影響幾乎可以忽略不計。總體來說,由于草地、水域等生態系統服務價值中貢獻較高的土地都在不同程度的減少,而不利于生態系統服務價值增加的建設用地急速擴張,導致30年間研究區整體生態系統服務價值下降明顯,減少率為81.25%。林地、水域和草地是生態系統服務價值的主要來源,嚴格保護林地、水域和草地等生態用地,對保護研究區的生態安全具有重要作用。

表3 1990—2020白洋淀流域各類用地生態系統服務價值量變化/108元Table 3 Ecosystem service value changes of each land use type of Baiyangdian basin from 1990 to 2020

表4反應出了單項生態系統服務價值的變化規律。依據每項服務ESV的貢獻程度大小依次為水文調節、氣候調節、土壤保持、氣體調節、生物多樣性、凈化環境、糧食生產、美學景觀、原料生產、水資源供應和維持養分循環。其中,水文調節和氣候調節的ESV在生態系統服務價值中所占比重較大,兩者占總價值的50%以上,以2020年為例水文調節的ESV為202.11億元；氣候調節的ESV為173.64億元。糧食生產主要受耕地面積變化影響,呈逐年下降趨勢；原料生產和維持養分循環受各類用地綜合影響,呈逐年下降的變化趨勢；除此之外水資源供應、氣體調節、氣候調節、凈化環境、水文調節、生物多樣性以及美學景觀等服務項目整體為波動下降的變化趨勢,這主要是林地、草地、水域等生態系統服務價值較高地類的面積波動造成的。

表4 1990—2020白洋淀流域不同生態系統服務類型的生態系統服務價值變化/108元Table 4 Ecosystem service value changes of each Ecosystem service type of Baiyangdian basin from 1990 to 2020

3.2 多情境下研究區土地利用變化及生態系統服務價值分析

3.2.1多情景下土地利用變化模擬

(1)精度檢驗。將2020年的預測圖與2020年實際土地利用類型進行對比,計算得到Kappa值為0.8674,通過精度檢驗(圖3)。如圖3所示,2020年模擬圖與2020年實際圖各地類面積差別較小,耕地、林地面積誤差為5%左右,建設用地面積誤差約為3%,草地面積誤差相對較大,但由于草地總面積相對較小,且圖上草地集中位置基本一致,所以也符合模擬要求,水域與未利用地基數較小且誤差微小。

圖3 2020年白洋淀流域土地利用實際圖和模擬圖Fig.3 Land use status map and forecast map of Baiyangdian basin in 2020

(2)多情景模擬預測。以2020年為基期圖像,分別模擬出3種情景下2035年研究區土地利用圖像(圖4)。根據研究區2010和2020年兩期土地利用變化情況,在趨勢延續情景模擬的基礎上加入雄安新區規劃建設這一變化條件,建立了3種土地利用變化模擬情景(圖4,表5)。在趨勢延續情景下2035年研究區建設用地面積為510582hm2,與2020年相比增加了89285.51hm2,林地面積小幅度增加,相比2020年增加了3.3%。耕地變化較大,下降為2020年的94.97%。相比之下草地、未利用地面積都有不同程度的減少,水域面積小幅度增加(圖4 情景1)。均衡發展情景下,由雄安2035年發展規劃可知,建設用地、林地、水域等地類面積都有所增加,整體來看,建設用地面積為513125.73hm2,林地比照趨勢延續情景相應增加,為1233968.7hm2,部分耕地被建設用地侵占。鑒于2035年雄安藍綠空間規劃,水域面積相應增加為52756.07hm2,由于雄安新區草地較少,所以兩種情景對比下,草地面積變化不明顯(圖4 情景2)。統籌發展情景下,研究區建設用地面積減少為500750.35 hm2,林地、草地面積相應增加,分別為1243886.14hm2、164611.76 hm2,耕地、水域和未利用地面積幾乎不變(圖4 情景3)。

從空間分布來看,與2020年相比,趨勢延續情景下白洋淀西面保定市中心城市范圍急速擴張,建設用地范圍激增,相應的侵占了大部分耕地；除此之外,研究區東北、正東、西南以及平原腹地等方向上的縣區市中心均相應擴張；西部山區縣域范圍中地勢相對平緩地區以及大清河河流直線周圍區域零碎耕地大部分轉化為林地,研究區西北部壩上平原地區以及研究區西南山地平原交界區域草地、耕地明顯減少；而在均衡發展情景下,由于雄安新區建設起步區位于白洋淀以西,容城縣與安新縣交匯處,所以該地建設用地和林地面積有所上升,白洋淀水域也相應擴張,雄安新區范圍內林帶增多,且多圍繞水域和起步區建設；在統籌發展情景下,為維護白洋淀流域生態格局的平衡與安全,在建設用地指標有限的前提下優先發展雄安新區,上游部分的縣區建設用地有所減少,城市擴張程度相應收縮。

研究區2035年趨勢延續情景為參照2010年、2020年兩期土地利用變化規律,沒有外部干預的基礎上模擬得到的,所以趨勢延續情景符合白洋淀流域城鎮化快速發展的變化趨勢。但是雄安新區的規劃建設使得部分地區的發展進程減速,尤其是白洋淀上游部分縣域,為維護整個流域的土地利用格局的平衡,限制了林地、草地向建設用地的轉換,為整個流域的可持續發展做出了貢獻。

3.2.2多情境下研究區生態系統服務價值變化特征

土地生態系統服務價值主要與各土地利用類型面積有關。與2020年各類用地的生態系統服務價值相比,2035年研究區趨勢延續情境下,林地、未利用地變化不明顯；耕地、草地生態系統服務價值下降幅度小于林地和水域增加的生態系統服務價值,使得趨勢延續情景下生態系統服務總價值與2020年相比有微量增加。與趨勢延續情景相比,其他兩個情景的生態系統服務價值變化有所不同。由表6可知,均衡發展情景下,林地和水域生態系統服務價值相對于趨勢延續情景高24.38億元,這主要與雄安新區規劃建設中藍綠空間增加的約束條件有關。統籌發展情景下,林地、草地和水域的生態系統服務價值均高于均衡發展情景下的生態系統服務價值,總計增加了5.2億元,該情景下協調了各類用地面積,林地、草地、水域都得到了相應的保護,從而使得生態系統服務價值有所增長。

表6 模擬2035年3種發展情境下各類土地生態系統服務價值/108元Table 6 The service value of each land ecosystem under three development scenarios will be simulated in 2035

為進一步分析多情景下研究區生態系統服務價值及其空間變化,利用ArcGIS軟件中的漁網工具,創建一個5km×5km正方形網格,計算每個網格的生態系統服務價值,制作2020年和三種模擬情景下2035年研究區的生態系統服務價值空間分布圖,采用自然斷點法將研究區生態系統服務價值從低到高劃分為6個等級,Ⅰ級最低,Ⅵ級最高(圖5)。

從2020年研究區生態系統服務價值空間分布來看,生態系統服務價值較高區域主要分布在西北部山區縣域。從整體上看,研究區生態系統服務價值整體呈西北高、東南低的分布特征,與研究區的海拔和地形地貌分布特征在空間上較為契合。主要原因在于,研究區西北部海拔較高,多為山地,地類主要以林地、草地為主,所以ESV較高；而東南部為華北平原,海拔較低,主要以耕地、建設用地為主,所以ESV較低。

與2020年相比,趨勢延續情景下,ESV高值區域總體來看小幅度減少,低值區域明顯增多,白洋淀上游區域,研究區東南方向縣域Ⅰ級區域明顯擴張,原Ⅰ級區域都有不同程度的擴張,這主要是由于城市的快速發展造成耕地向建設用地轉變。而在研究區西北部的縣區中部分地區由Ⅴ級向Ⅵ級轉變,這是由于縣中心建設用地面積相對增加,導致林地、草地面積相對減少,ESV下降。

對2035年研究區3種情境下ESV空間分布進行對比分析發現,均衡發展情景和統籌發展情景總體ESV明顯好于趨勢延續情景。東部雄安新區范圍由于新增許多林帶,并且水域范圍得到擴大,ESV高值區域增多,白洋淀區域由Ⅰ、Ⅱ、級向Ⅲ、Ⅳ、Ⅴ級轉變；統籌發展情景與均衡發展情景相比,統籌發展情景總體ESV優于均衡發展情景,由圖5可知,西北部山地區域Ⅵ級范圍增多、華北平原Ⅰ級區域減少,這是由于部分建設用地減少,林地面積明顯增加。

圖5 不同情境下2030年白洋淀流域ESV空間分布Fig.5 Spatial distribution of ESV in Baiyangdian basin in 2035 under different circumstances

3.3 多種情境下白洋淀流域生態系統服務價值的熱點分析

為了反映了生態系統服務高值和低值在空間上的集聚情況,在多次嘗試后選擇2km×2km網格單元,基于ArcGIS軟件,利用空間統計工具中的熱點分析(Getis-Ord Gi*),得到白洋淀流域生態系統服務價值的熱點和冷點區域(圖6)。

圖6 不同情境下2030年白洋淀流域ESV冷熱點分布Fig.6 Distribution of ESV hot spots in Baiyangdian basin in 2035 under different conditions

2020年研究區ESV值的熱點和冷點分布范圍較大,且集聚性較強。熱點區域主要集中在西北部山地區域,且幾乎為99%置信度。冷點區域集中連片,主要分布在研究區東南部平原地區、研究區北部蔚縣、涿鹿縣部分區域,其中99%置信度的冷點區域多為各縣區的中心建設區。與2020年相比,2035年三種模擬情景總體特征為：熱點區域分布范圍小幅度擴大,冷點區域中99%置信度冷點區范圍擴大,95%置信度和90%置信度冷點區范圍縮小,不顯著區域范圍增多。具體到各個情景,在趨勢延續情景下,研究區南部不顯著區域小部分轉換為熱點區域,西北部山區熱點區域向外圍擴張,這是由于該區域林地、草地增多,ESV上升所致。研究區東南部平原地區各縣區中心建設區與99%置信度冷點區域的空間分布特征相契合,且由于城市的擴張與建設需要,大量耕地轉換為建設用地,導致其在空間上的集聚程度提升,99%置信度冷點區域擴張,冷點區域整體呈集中連片的空間分布特征。在均衡發展情境下,除雄安以外其他區域冷熱點分布情況與趨勢延續情景相似,不同之處在于,雄安新區三縣的邊界交匯處出現了集聚性較高的熱點區,且雄安新區大部分冷點區域轉變為不顯著區域,這是由于響應雄安建設中藍綠空間達到70%的規劃綱要和白洋淀生態環境治理與保護的條例,新增熱點區域為起步區規劃建設林帶、水域集聚,導致ESV上升而形成的。在統籌發展情景下,研究區西北部白洋淀上游縣區熱點區域小范圍擴張；研究區西部淶源縣中心區域冷點區域縮小,集聚程度降低。對比分析三種模擬情景下的冷熱點空間分布規律,冷點區范圍：趨勢延續>均衡發展>統籌發展；熱點區范圍：統籌發展>均衡發展>趨勢延續。

4 結論與討論

4.1 討論

本研究在探討雄安新區建設引起的流域土地利用變化背景下,利用CA-Markov模型模擬趨勢延續、均衡發展和統籌發展等3種發展情景下2035年白洋淀流域的土地利用變化狀況,并評估白洋淀流域的ESV,在此基礎上分析ESV空間分布特征和集聚特征。不同發展情景下的流域生態系統服務價值差異較為明顯。趨勢延續情景延續了白洋淀流域在未設立雄安新區的背景下的用地需求,生態系統服務價值略有上升,基本實現流域生態文明建設與可持續發展。均衡發展情景兼顧了雄安新區與流域內其他縣市的建設用地需求,以適應于城鎮化的推進與增強社會經濟的發展的主要目標,生態系統服務價值略有上升,一定程度上有利于生態保護；統籌發展情景凸顯了國家戰略與生態約束高水平融合,滿足雄安新區土地利用急劇轉型的同時保障生態空間的面積。三種發展情景為未來流域生態建設與用地擴張提供了重要的參考,結合雄安新區發展戰略與流域實際發展需求選擇相應的情景以期為區域提供適用于未來用地空間發展模式的科學依據及理論支撐。白洋淀流域內市、縣的城鎮化水平在未來一段時間還將繼續上升,建設用地的快速擴張時期尚未完全結束,且又地處于資源環境承載能力較弱、水資源超載、景觀破碎、空氣質量治理任重道遠的華北地區,提升白洋淀流域的生態系統服務功能需求十分迫切[5,47]。因此,雄安新區作為生態城市建設的標桿,有必要探討在高質量發展背景下制定精細化、差異化的國土空間格局優化調控措施,協調與流域其他區域的發展目標。

4.2 結論

白洋淀流域由于獨特的地理位置,人類活動對流域內的土地利用結構變化有著顯著影響,本研究以土地利用變化為基礎,采用當量因子法估算4期生態系統服務價值,并借助CA-Markov模型,模擬不同情境土地利用格局,進而測算生態系統服務價值變化情況。主要結論如下：

(1)1990—2020年白洋淀流域土地利用變化最明顯的特征是耕地面積銳減,建設用地迅速擴張,30年間耕地面積減少了122776.99hm2,建設用地增加了167369.74hm2。

(2)1990—2020年研究區生態系統服務價值分別為852.59億元、788.26億元、796.95億元和784.68億元,總體呈波動下降的趨勢。30年間白洋淀流域快速城市化過程導致城市建設用地擴張侵占大量耕地、草地和水域等生態用地,使人地關系的矛盾加劇,破壞城市生態環境,甚至威脅城市生態安全。

(3)在對白洋淀流域2035年3種情景模擬后發現,趨勢延續情景、均衡發展情景、統籌發展情景的生態系統服務價值依次為785.36億元、807.26億元、812.46億元,均高于2020年研究區生態系統服務價值。趨勢延續情景下,由于城市的快速擴張,平原地區耕地、水域受到影響,山地地區耕地、草地少部分轉化為林地和建設用地；均衡發展情景下,由于雄安建設規劃和白洋淀生態環境保護的要求,雄安新區范圍內林地、水域增加,其他地區保持不變；在統籌發展的情景下,白洋淀流域上游縣域城市建設放緩,生態用地增多,為整個流域的生態安全做出了貢獻。

(4)從生態系統服務價值空間分布來看,白洋淀流域生態系統服務價值整體呈現西北高東南低的特點。不同情境下ESV熱點分布情況。總體來說,研究區ESV值的熱點和冷點分布范圍較大且集聚性較強。