夏圣潔,陳慧儒,張鈞韋,劉艷紅 (山西農(nóng)業(yè)大學(xué)城鄉(xiāng)建設(shè)學(xué)院,山西 太谷 030801)

隨著城市發(fā)展需求的增長,建設(shè)用地的快速擴(kuò)張大量占用了森林、水域等生態(tài)用地,使得城市區(qū)域內(nèi)的熱環(huán)境問題變得更加突出,人居環(huán)境也面臨著嚴(yán)峻的挑戰(zhàn)[1-2].此外,在新型城鎮(zhèn)化背景下城市群已經(jīng)被確立為發(fā)展的主體形態(tài),隨著城市群規(guī)模的拓展,城市間距離不斷縮小,集中連片的城市發(fā)展模式對(duì)區(qū)域生態(tài)環(huán)境的疊加影響持續(xù)加大,城市群的熱環(huán)境問題將變得更為復(fù)雜[3-4].而生態(tài)用地作為城市生態(tài)系統(tǒng)的重要組成部分,在城市溫度調(diào)節(jié)中發(fā)揮著重要作用.面對(duì)日益嚴(yán)重的熱環(huán)境問題,研究城市群的生態(tài)用地演變與熱環(huán)境的互動(dòng)關(guān)系有著重要意義.

相關(guān)領(lǐng)域?qū)W者們已從多個(gè)角度開展了大量的研究工作[5-7],目前基于城市熱島效應(yīng)的研究通常以局部和單一城市為主[8].最近的一項(xiàng)綜述總結(jié)了1972～2018年期間的相關(guān)研究,其中單個(gè)城市相關(guān)研究占 64%,而對(duì)城市群熱環(huán)境的關(guān)注明顯較少(14%)[9].然而,單一城市的熱環(huán)境問題和相應(yīng)的驅(qū)動(dòng)因素不能完全等同于城市群,二者之間的差異還有待繼續(xù)探索[10],因此,在城市群尺度上探討熱環(huán)境問題將是城市化新形勢背景下的新課題[11].另外,國內(nèi)外對(duì)生態(tài)系統(tǒng)熱消減的研究多為不同類型的綠地與水體[12-14]、不同時(shí)段的熱消減作用[15]以及景觀格局與熱環(huán)境之間的耦合關(guān)系[16]等方面,鮮有在城市群尺度上開展生態(tài)系統(tǒng)降溫作用的定量研究.此外,由于空間相互作用和空間擴(kuò)散的影響,熱環(huán)境的分析變量往往不能完全獨(dú)立存在,而是存在著空間自相關(guān)性[17-18].以往的研究,多采用傳統(tǒng)統(tǒng)計(jì)學(xué)方法,對(duì)有關(guān)影響熱環(huán)境的因子進(jìn)行線性擬合等定量研究[19-20],極少關(guān)注熱環(huán)境與其影響變量的空間相關(guān)性特征,將城市群生態(tài)用地演變與地表溫度有機(jī)結(jié)合,運(yùn)用空間自相關(guān)分析的方法,繪制聚類圖,揭示研究區(qū)生態(tài)用地與熱環(huán)境的局部空間集聚動(dòng)態(tài)變化規(guī)律,能夠更為直接地反映其空間自相關(guān)特征及變化趨勢.

基于此,本研究以山西中部城市群為研究區(qū)域,選取2010～2022 年間4 期Landsat 遙感影像數(shù)據(jù),選用均衡度、熱貢獻(xiàn)度等指數(shù)及轉(zhuǎn)移矩陣模型、重心遷移模型、空間自相關(guān)等方法,在分析生態(tài)用地與地表熱環(huán)境的時(shí)空動(dòng)態(tài)變化基礎(chǔ)上,對(duì)二者耦合的空間相關(guān)性開展定量研究.揭示城市群尺度上生態(tài)用地動(dòng)態(tài)變化對(duì)熱環(huán)境的影響機(jī)制,量化不同類型的生態(tài)用地對(duì)城市群熱環(huán)境的貢獻(xiàn)度.為有效緩解城市群范圍內(nèi)熱環(huán)境問題、優(yōu)化城市群內(nèi)部空間結(jié)構(gòu),構(gòu)筑生態(tài)和安全屏障,形成多中心、多層級(jí)、多節(jié)點(diǎn)的網(wǎng)絡(luò)型城市群提供理論依據(jù)與方法支撐.

1 數(shù)據(jù)與方法

1.1 研究區(qū)概況

山西中部城市群地處山西中部(39°39'～36°38',114°9～110°23'),包括太原、晉中、忻州、呂梁、陽泉五市,國土面積7.41 萬km2,占全省的47.3%(圖1).地勢東西高中部低.海拔350～3061.1m,屬暖溫帶半干旱大陸性季風(fēng)氣候,年均溫10～11℃,年平均降水量400～650mm.山西中部城市群地處省域中心地帶,在《山西中部城市群高質(zhì)量發(fā)展規(guī)劃(2022～2035年)》中是全省構(gòu)建“一軸、兩屏、多廊”生態(tài)格局的核心組成部分[21].

圖1 研究區(qū)概況示意Fig.1 Overview map of the study area

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

1.2.1 生態(tài)用地?cái)?shù)據(jù) 本研究以山西中部城市群2010、2014、2018、2022 四期Landsat 遙感影像為源數(shù)據(jù),分辨率均為30m.影像從Geospatial data cloud 與USGS 網(wǎng)站獲取,每一期選取在植物生長季時(shí)期、云少清晰的影像.

利用ENVI5.3 進(jìn)行輻射定標(biāo)、大氣校正、影像裁剪與鑲嵌等預(yù)處理.采用目視解譯[22]最大似然法[23]進(jìn)行監(jiān)督分類解譯以及分類后處理提取得到林地、草地、耕地、水體、建設(shè)用地、未利用地6個(gè)用地類型.Kappa 系數(shù)均大于0.80,符合精度要求.基于生態(tài)用地研究成果及研究區(qū)現(xiàn)狀,將山西中部城市群區(qū)域分為3 大生態(tài)用地類型:基礎(chǔ)性生態(tài)用地、敏感性生態(tài)用地、輔助性生態(tài)用地(表1)[24-25].

表1 生態(tài)用地類型劃分Table 1 Breakdown of ecological land types

1.2.2 地表溫度數(shù)據(jù) 為反映地表溫度(LST)的階段性變化,研究參考相關(guān)文獻(xiàn).以12 年時(shí)間為尺度,選取與生態(tài)用地同幅影像,利用ENVI 5.3 軟件,對(duì)數(shù)據(jù)進(jìn)行輻射定標(biāo)、大氣校正和影像裁切等預(yù)處理.采用覃志豪單窗算法[26]反演地表溫度.為了提高不同時(shí)期地表溫度時(shí)空研究的科學(xué)性,采用均值-標(biāo)準(zhǔn)差法對(duì)地表溫度反演結(jié)果進(jìn)行分級(jí),共5 級(jí)(表2),其中μ為平均值,σ為標(biāo)準(zhǔn)差.

表2 溫度分區(qū)劃分Table 2 Temperature partitioning

1.3 研究方法

1.3.1 土地利用轉(zhuǎn)移矩陣模型 土地利用轉(zhuǎn)移矩陣模型可以反映研究期內(nèi)兩種地類之間的轉(zhuǎn)移數(shù)量和方向,是研究山西中部城市群生態(tài)用地變化的有效手段[27].

計(jì)算表達(dá)式為:

式中:Sij表示研究期i到j(luò)的用地類型;n表示不同用地類型的數(shù)量.

1.3.2 土地利用均衡度 選取基于土地利用信息熵的土地利用均衡度來衡量生態(tài)用地類型的均衡程度.區(qū)域土地熵值大小反映土地利用類型的多少及各土地類型面積分布的均勻程度[28].土地利用結(jié)構(gòu)信息熵的計(jì)算公式可寫為:

基于信息熵函數(shù)可以得出均衡度的公式[29]:

式中:J為均衡度;H是實(shí)際熵值;Hmax為最大熵值,二者之比表示區(qū)域土地利用的均衡程度.

1.3.3 重心遷移模型 重心遷移模型可以揭示各類生態(tài)用地在遷移距離、遷移角度和遷移方向上的空間變化特征.通過計(jì)算不同時(shí)期各生態(tài)用地類型的重心遷移來研究生態(tài)用地的空間變化[30].相關(guān)計(jì)算公式如下:

式中:Xt和Yt分別為t時(shí)期某種生態(tài)用地類型的經(jīng)度和緯度;Cti為t時(shí)期某種生態(tài)用地類型第i個(gè)斑塊的面積;Xti和Yti分別為t時(shí)期生態(tài)用地類型第i個(gè)斑塊的幾何重心經(jīng)度和緯度坐標(biāo).

1.3.4 空間自相關(guān)分析 空間自相關(guān)分析是研究鄰近位置屬性相關(guān)性的空間統(tǒng)計(jì)學(xué)方法,本研究中主要分析生態(tài)用地類型與地表溫度二者在空間上的離散或聚集等分布特性[31].本文采用Moran’sI指數(shù)和Lisa 指數(shù)來表征生態(tài)用地與熱環(huán)境間的空間自相關(guān)關(guān)系.Moran’sI指數(shù)主要反映全部數(shù)據(jù)的某項(xiàng)屬性值在空間中的相關(guān)性大小,其取值一般在(-1,1)之間,小于0 表示負(fù)相關(guān),等于0 表示不相關(guān),大于0表示正相關(guān);Lisa指數(shù)的具體表現(xiàn)為:LL型(低低聚集)、HH 型(高高聚集)、HL 型(低值包高值聚集)、LH 型(高值包低值聚集)4 種類型.

2 結(jié)果與討論

2.1 城市群生態(tài)用地演變特征

本研究所指生態(tài)用地指除人工硬化表面之外,其他能夠直接或間接提供環(huán)境調(diào)節(jié)和生物支持等生態(tài)系統(tǒng)服務(wù)功能,且自身具有一定的自我調(diào)節(jié)、維持、修復(fù)和發(fā)展能力的土地.其中,本文主要研究的生態(tài)系統(tǒng)服務(wù)功能為生態(tài)用地的降溫效應(yīng).

利用ArcGIS 中spatial analyst-Reclassify 工具將土地利用類型進(jìn)行重新映射,最終得到各年生態(tài)用地分布圖(圖2).由圖2 看出,基礎(chǔ)性生態(tài)用地面積最大,且分布較為廣泛;輔助性生態(tài)用地主要集中分布在中部盆地,多圍繞在城市建成區(qū)附近,以城市為中心分布.

圖2 山西中部城市群生態(tài)用地分布Fig.2 Distribution of ecological land in urban agglomeration in central Shanxi

以山西中部城市群生態(tài)用地分布數(shù)據(jù),借助ArcGIS 及Origin 進(jìn)行數(shù)據(jù)分析,得到2010～2022年間生態(tài)用地轉(zhuǎn)移?；鶊D(圖3)總體來看,基礎(chǔ)性生態(tài)用地與輔助性生態(tài)用地占據(jù)了整個(gè)用地類型轉(zhuǎn)變的大部分.其中基礎(chǔ)性生態(tài)用地轉(zhuǎn)出最多的為2010 年,輔助性生態(tài)用地轉(zhuǎn)出最多的為2018 年.由此得出,變化最大的階段為2010～2014 年,2018～2022 年.

圖3 2010～2022 年山西中部城市群生態(tài)用地轉(zhuǎn)移矩陣Fig.3 Ecological land transfer matrix of urban agglomeration in central Shanxi from 2010 to 2022

2.2 城市群熱環(huán)境動(dòng)態(tài)演變特征

對(duì)分級(jí)后的城市群熱環(huán)境分級(jí)(圖4)進(jìn)行分析統(tǒng)計(jì).總體而言,12a 間研究區(qū)范圍內(nèi)熱環(huán)境變化復(fù)雜,低溫區(qū)與次低溫區(qū)面積有略微增加,中溫區(qū)面積增長至30637.43km2,占比增加了3.3%,高溫區(qū)增長較小,面積占比增加至11.63%,次高溫區(qū)萎縮較大,面積縮小了412864km2,占比減少5.75%.

對(duì)4 個(gè)年份熱環(huán)境分區(qū)進(jìn)行相交疊加分析,得到熱環(huán)境變化分區(qū)(圖5).總體來看,2010～2022 年熱環(huán)境變化趨勢復(fù)雜,呈惡化-有所緩解-改善明顯; 2010～2022 年升溫區(qū)域面積達(dá)14296.49km2,占比 19.40%,降溫區(qū)域面積達(dá) 17809.82km2,占比24.16%,降溫面積高于升溫面積,熱環(huán)境總體處于改善狀態(tài),從分布狀況來看,升溫區(qū)域多集中分布于城區(qū)中心建設(shè)用地即非生態(tài)用地,周邊耕地及草地降溫明顯.

圖5 2010～2022 年熱環(huán)境變化分區(qū)Fig.5 Zoning map of thermal environment change from 2010 to 2022

2.3 生態(tài)用地動(dòng)態(tài)演變與熱環(huán)境空間自相關(guān)分析

2.3.1 生態(tài)用地的熱貢獻(xiàn)度分析 熱效應(yīng)貢獻(xiàn)度(HI)是指不同的生態(tài)用地類型對(duì)區(qū)域平均溫度的影響程度,即對(duì)區(qū)域溫度的貢獻(xiàn)度,熱效應(yīng)貢獻(xiàn)度的正負(fù)分別代表升溫和降溫作用,絕對(duì)值的大小則代表升溫和降溫的強(qiáng)度[32].

通過計(jì)算各類生態(tài)用地的熱效應(yīng)貢獻(xiàn)指數(shù)(表3),得到各類用地?zé)嵝?yīng)貢獻(xiàn)度指數(shù)表格,從表格看出各類生態(tài)用地中,基礎(chǔ)性生態(tài)用地?zé)嵝?yīng)貢獻(xiàn)指數(shù)均小于零,非生態(tài)用地以及輔助性生態(tài)用地貢獻(xiàn)指數(shù)均大于零,敏感性生態(tài)用地?zé)嶝暙I(xiàn)指數(shù)除2018 年之外均大于零.說明基礎(chǔ)性生態(tài)用地起到明顯降溫作用;輔助性生態(tài)用地與非生態(tài)用地均為升溫效應(yīng).其中,2022 年基礎(chǔ)性生態(tài)用地?zé)嵝?yīng)貢獻(xiàn)指數(shù)達(dá)-2.87,表明2022年基礎(chǔ)性生態(tài)用地的熱環(huán)境的降溫作用最顯著;非生態(tài)用地的熱效應(yīng)貢獻(xiàn)指數(shù)隨著時(shí)間推進(jìn)不斷減少,表明對(duì)熱環(huán)境增溫效應(yīng)逐漸減弱;輔助性生態(tài)用地在2010、2018 年間的熱貢獻(xiàn)指數(shù)較小,但整體處于較高水平,對(duì)熱環(huán)境的增溫效應(yīng)相比非生態(tài)用地更加顯著.

表3 2010～2022 年各生態(tài)用地?zé)嵝?yīng)貢獻(xiàn)指數(shù)Table 3 Geothermal effect contribution index by ecological use, 2010～2022

2.3.2 生態(tài)用地均衡度與地表溫度(LST)的空間自相關(guān)分析 計(jì)算土地利用均衡度來衡量生態(tài)用地結(jié)構(gòu)均衡程度,利用城市群生態(tài)用地均衡度以及城市群地表溫度數(shù)據(jù)在Geoda 中進(jìn)行雙變量空間自相關(guān)分析,結(jié)果如圖6 所示,12a 間4 階段莫蘭指數(shù)分別為-0.264、-0.250、-0.107、-0.218,顯著性P值均小于0.001,數(shù)據(jù)分析結(jié)果可靠.由Lisa 散點(diǎn)圖可知,4 階段生態(tài)用地均衡度與溫度均呈負(fù)相關(guān)性,即生態(tài)用地均衡度越高地表溫度越低.

圖6 2010～2022 年生態(tài)用地均衡度與LST Lisa 散點(diǎn)圖及聚類分布Fig.6 Ecological land use balance and LST Lisa scatter plot and cluster distribution map from 2010 to 2022

由圖6 可知,生態(tài)用地均衡度與溫度具有顯著的空間聚集性關(guān)系,4a 中面積較多的為低-高、高-低聚集區(qū).其中高-低聚集性區(qū)域多為山林地區(qū),基礎(chǔ)性生態(tài)用地中植被覆蓋度較高且均衡度越高的區(qū)域.低-高聚集性區(qū)域多為城市中心建設(shè)用地、輔助性生態(tài)用地較多區(qū)域,高-高聚集性區(qū)域多為草地、耕地聚集性較高區(qū)域.生態(tài)用地均衡度中,基礎(chǔ)性生態(tài)用地均衡度越高地表溫度越低,位于城市周邊的輔助性生態(tài)用地均衡度越高地表溫度越高.

2.3.3 生態(tài)用地類型與地表溫度分區(qū)的空間自相關(guān)及其動(dòng)態(tài)變化分析 為了探究生態(tài)用地以及地表熱環(huán)境分區(qū)在空間上的相關(guān)性,將溫度分區(qū)中的高溫區(qū)與次高溫區(qū)合并為較高溫區(qū),低溫區(qū)與次低溫區(qū)合并為較低溫區(qū),并且基于以上研究,對(duì)生態(tài)用地類型以及溫度分區(qū)進(jìn)行重新賦值(表4).

表4 生態(tài)用地類型與溫度分區(qū)賦值Table 4 Ecological land type and temperature zoning assignment table

利用4 階段賦值后生態(tài)用地類型與溫度分區(qū)數(shù)據(jù),在Geoda 進(jìn)行雙變量空間自相關(guān)分析,結(jié)果如圖7 所示,顯著性P值均小于0.005,數(shù)據(jù)分析結(jié)果可靠.由圖7 可知,生態(tài)用地類型與溫度分區(qū)間具有顯著的空間聚集性關(guān)系,4a 中面積較多的為低-高,高-低聚集區(qū).其中高-低聚集性區(qū)域多為基礎(chǔ)性生態(tài)用地分布地區(qū)與中溫區(qū)及較低溫區(qū)對(duì)應(yīng).低-高區(qū)域多為城市中心非生態(tài)用地與輔助性生態(tài)用地對(duì)應(yīng)較高溫區(qū).可知,在空間分布上,基礎(chǔ)性生態(tài)用地的空間聚集性多與中溫區(qū)、較低溫區(qū)重合,非生態(tài)用地與輔助性生態(tài)用地多與較高溫區(qū)重合.

圖7 2010～2022 年山西中部城市群生態(tài)用地與溫度分區(qū)聚類分布Fig.7 Ecological land and temperature zoning cluster distribution map of urban agglomeration in central Shanxi from 2010 to 2022

為明確2010～2022 年間生態(tài)用地與各溫度分區(qū)在空間上的動(dòng)態(tài)變化趨勢,運(yùn)用重心遷移模型方法,繪制重心遷移圖譜(圖8),統(tǒng)計(jì)其遷移方向及距離(表5),研究生態(tài)用地類型與溫度二者的空間動(dòng)態(tài)變化關(guān)系.

表5 生態(tài)用地與地表溫度重心遷移耦合Table 5 Coupling of ecological land use and surface temperature center of gravity migration

圖8 生態(tài)用地與地表溫度重心遷移耦合圖Fig.8 Coupling diagram of ecological land use and surface temperature center of gravity migration

其中,基礎(chǔ)性生態(tài)用地的遷移軌跡在較低溫區(qū)與中溫區(qū)之間對(duì)比而言,與中溫區(qū)軌跡更為擬合,在2014～2018 年均向偏南方向移動(dòng),2018～2022 年均向西北方向移動(dòng),表明基礎(chǔ)性生態(tài)用地的遷移與中溫區(qū)的遷移相關(guān)性更大,由于基礎(chǔ)性生態(tài)用地面積占比大,較低溫區(qū)占比較低,中溫區(qū)面積占比與基礎(chǔ)性生態(tài)用地面積占比在整體空間動(dòng)態(tài)演變中更為接近;較高溫區(qū)遷移軌跡與輔助性生態(tài)用地遷移軌跡較為擬合,在2010～2014 年間均向偏北方向移動(dòng),2018～2022年均向偏西方向移動(dòng),說明較高溫度地區(qū)的遷移受輔助性生態(tài)用地的遷移影響較大.

整體來看,基礎(chǔ)性生態(tài)用地、輔助性生態(tài)用地、非生態(tài)用地以及中溫區(qū)和較高溫區(qū)在2010～2022 年整體向北偏西方向進(jìn)行遷移,從太原市中心向西北部發(fā)展.總體來說,輔助性生態(tài)用地、非生態(tài)用地與較高溫區(qū)域的遷移更加相關(guān),基礎(chǔ)性用地與較低溫區(qū)、中溫區(qū)的關(guān)聯(lián)性更強(qiáng).重心遷移軌跡多集中在太原城市中心及太原西北部山脈林地和忻州城市中心,即規(guī)劃中構(gòu)建城市群一體化發(fā)展新格局的太原榆次太谷城市群發(fā)展核心以及雄忻發(fā)展帶所在地.其中較低溫區(qū)軌跡多集中在太原市東北部(陽曲縣境內(nèi),五臺(tái)山系及呂梁山系太原西山山脈)以及忻州市區(qū),其他用地以及溫度區(qū)轉(zhuǎn)移多集中在太原市西部(婁煩縣及古交市境內(nèi))并且非生態(tài)用地在2010～2014 年間重心均在太原市城區(qū)中心.研究期間均向西北部偏移,即向規(guī)劃中構(gòu)建“一軸、兩屏、多廊”生態(tài)格局的汾河景觀生態(tài)綠軸進(jìn)行靠近.

2.4 討論

2.4.1 城市化背景下城市群生態(tài)用地及其熱環(huán)境的動(dòng)態(tài)演變趨勢分析 隨著城市建設(shè)進(jìn)程加快,城市生態(tài)用地也隨之改變,在研究期間內(nèi),地類變化最劇烈的為基礎(chǔ)性生態(tài)用地與輔助性生態(tài)用地之間的相互轉(zhuǎn)化,這與Zhu 等[33]研究結(jié)果一致;結(jié)合國情來看,在2010～2014 年間,中國受到經(jīng)濟(jì)全球化影響,快速發(fā)展,導(dǎo)致城市化進(jìn)程加快,輔助性生態(tài)用地急速減少,2014～2022隨著中國退耕還林還草政策的推進(jìn),輔助性生態(tài)用地逐步轉(zhuǎn)化為基礎(chǔ)性生態(tài)用地.

城市化不僅加速了生態(tài)用地的變化還改變了城市近地表的大氣成分[34].導(dǎo)致了城市熱島氣候現(xiàn)象的產(chǎn)生;在全球化的大背景下,城市化進(jìn)程的加快進(jìn)一步加劇了熱島效應(yīng)[35-36].而地表溫度(LST)是城市熱環(huán)境最直接的表征[37],本研究基于LST 的城市群熱環(huán)境演變分析結(jié)果表明,城市群熱環(huán)境呈現(xiàn)復(fù)雜化的趨勢,在2010～2022 年間研究區(qū)城市群熱環(huán)境整體略微改善,但其變化特征變得更為復(fù)雜,與Feng 等[38]研究結(jié)果一致.從熱環(huán)境分布狀況來看,城市群中城市建成區(qū)集中地區(qū)溫度較高,有中心城區(qū)向外LST 逐漸降低,與Liu 等[39]的研究相同.

在城市化的大背景下,城市的生態(tài)環(huán)境、熱環(huán)境問題變得更加突出,人居環(huán)境面臨著嚴(yán)峻的挑戰(zhàn),生態(tài)用地與熱環(huán)境的變化在城市群尺度上變得更為復(fù)雜,因此明晰二者的時(shí)空演變特征,是緩解城市群熱環(huán)境、改善城市群生態(tài)環(huán)境的重要基礎(chǔ).

2.4.2 城市群生態(tài)用地與熱環(huán)境相關(guān)性分析 以往對(duì)于城市熱環(huán)境與其因素的相關(guān)性分析以普通線性回歸分析等方法居多,采用到空間自相關(guān)分析方法的多針對(duì)于某具體的單一方面,比如Jiayi 等[40]利用地表溫度熱狀況指數(shù)(TCI)數(shù)據(jù)對(duì)中國北部、南部和西北部不同地物的熱環(huán)境分異進(jìn)行了空間自相關(guān)分析,結(jié)果表明2018～2020 年華北城市以高-高集聚型為主.

針對(duì)于二者耦合的空間相關(guān)性分析多采用地表溫度與土地開發(fā)強(qiáng)度[41]、生態(tài)用地變遷[38]等相關(guān)指標(biāo)進(jìn)行,且相關(guān)研究表明城市群的土地開發(fā)強(qiáng)度(LDI)和LST 具有較強(qiáng)的空間正相關(guān)性;LST 與城市生態(tài)用地變遷(UELTs)有顯著的時(shí)空耦合關(guān)系;城市生態(tài)用地(UEL)的破碎化、城市人口壓力與建設(shè)用地壓力(CLP)增強(qiáng)、綠地減少以及UEL 質(zhì)量下降促使了LST上升.CLP、UEL退化與LST呈正相關(guān),UEL恢復(fù)與LST 呈負(fù)相關(guān).

本研究結(jié)果表明:在熱效應(yīng)貢獻(xiàn)上,各生態(tài)用地中基礎(chǔ)性生態(tài)用地降溫效應(yīng)最好,熱效應(yīng)貢獻(xiàn)指數(shù)平均值為-1.29,非生態(tài)用地平均溫度最高且起升溫效應(yīng),熱效應(yīng)貢獻(xiàn)指數(shù)平均值為0.33;在空間遷移變化上,輔助性生態(tài)用地、非生態(tài)用地與較高溫區(qū)域的遷移更加相關(guān),基礎(chǔ)性用地與較低溫區(qū)、中溫區(qū)的關(guān)聯(lián)性更強(qiáng);均與以往研究結(jié)果一致.因此,合理布局城市群內(nèi)各類生態(tài)用地,是有效緩解城市群熱環(huán)境問題的重要思路之一.

此外,在本研究中引入生態(tài)用地均衡度,將生態(tài)用地結(jié)構(gòu)與地表溫度相結(jié)合,并對(duì)二者進(jìn)行空間自相關(guān)性分析,探究城市群生態(tài)用地與熱環(huán)境的空間相關(guān)性特征,補(bǔ)充完善了城市群生態(tài)用地的熱消減理論,為改善城市群熱環(huán)境實(shí)踐提供理論依據(jù).

3 結(jié)論

3.1 城市群中各生態(tài)用地均有不同程度變化,其結(jié)構(gòu)、空間動(dòng)態(tài)變化規(guī)律與城市群整體發(fā)展規(guī)劃相對(duì)應(yīng).其中變化最為頻繁的為基礎(chǔ)性生態(tài)用地與輔助性生態(tài)用地.

3.2 城市熱環(huán)境方面,2010～2022 年間研究區(qū)范圍內(nèi)熱環(huán)境變化復(fù)雜:面積增加最多為高溫區(qū),占11.63%,次高溫區(qū)面積縮小最多占比減少了5.75%;從空間變化上看,升溫區(qū)域占比19.40%,降溫區(qū)域占比24.16%.

3.3 城市群生態(tài)用地從熱效應(yīng)貢獻(xiàn)度、用地均衡度與空間分布的變化與城市群熱環(huán)境有著密切的空間關(guān)系:(1)在熱效應(yīng)貢獻(xiàn)上,各生態(tài)用地中基礎(chǔ)性生態(tài)用地起降溫效應(yīng),非生態(tài)用地起升溫效應(yīng).(2)生態(tài)用地的均衡度在空間上存在著顯著聚集性且與地表溫度呈負(fù)相關(guān),表明生態(tài)用地均衡度越高地表溫度越低,且基礎(chǔ)性生態(tài)用地均衡度越高地表溫度越低.(3)在空間分布與遷移軌跡上,生態(tài)用地類型與溫度分級(jí)分布上具有空間耦合性,基礎(chǔ)性生態(tài)用地的空間多與中溫區(qū)、較低溫區(qū)擬合,非生態(tài)用地與輔助性生態(tài)用地多與較高溫區(qū)擬合;重心的分布以及遷移軌跡與城市群發(fā)展規(guī)劃聯(lián)系緊密.