任圓圓,張學(xué)雷,李笑瑩,吳星宇
(1.鄭州輕工業(yè)大學(xué) 社會發(fā)展研究中心,河南 鄭州 450002;2.鄭州大學(xué) 水利科學(xué)與工程學(xué)院,河南 鄭州 450001)
地形作為重要的自然要素之一,其在一定程度上制約并影響著土地利用的方式和空間分布格局。土地利用變化反映人類在利用自然環(huán)境的過程中的時間上和空間上的趨勢和行為,其相關(guān)研究已受到全球環(huán)境領(lǐng)域研究者的關(guān)注[1-4]。土地利用變化的過程受到自然因素和人為因素的影響,而地形在一定程度上能夠揭示自然要素對土地利用格局變化的演變過程[5-7]。城鎮(zhèn)建設(shè)用地作為重要的土地利用方式之一,其一直是土地科學(xué)研究的熱點問題之一,探討城鎮(zhèn)建設(shè)用地擴(kuò)張?zhí)卣髯兓诘匦紊系难葑冞^程在一定程度上能夠為城市化過程中土地資源的合理有效利用提供數(shù)據(jù)支撐和現(xiàn)實參考。
關(guān)于地形和土地利用變化間的研究,F(xiàn)ang等[8]在復(fù)雜的地形區(qū)域研究了植被帶動力學(xué)下的土地利用變化對風(fēng)能潛力影響的變化;B Jiménez-Esteve等[9]在一個復(fù)雜的地形區(qū)域探討土地利用和地形對高分辨率中尺度WRF模型的影響;MU Songlin等[10]基于地形因子對北京山區(qū)溝域的土地利用變化進(jìn)行研究;申懷飛等[11]研究了河南省西部黃河流域土地利用空間格局在不同地形梯度上的分異規(guī)律;路昌等[12]探討了撓力河流域地形因素對土地利用格局變化的影響,以揭示土地利用在時間和空間上演變的過程、機(jī)理和形成的原因;張靜等[13]選擇垂直地帶性強的秦巴山區(qū)分析其土地利用的垂直景觀格局、土地利用的強度和變化;宋戈等[14]研究了黑龍江省巴彥縣的土地利用空間格局在不同地形梯度上的特征,并重點分析地形因子對土地利用空間產(chǎn)生的影響;張平等[15]探討了安徽省南部地區(qū)不同高程、坡度和地形位指數(shù)上的土地利用分布指數(shù),以揭示地形因子與土地利用分布和變化之間的內(nèi)在聯(lián)系;關(guān)于地形和城鎮(zhèn)建設(shè)用地之間的關(guān)系,楊傳俊等[16]獲取了粵北山區(qū)韶關(guān)市3個不同時期的建設(shè)用地信息,分析了丘陵山區(qū)上城鎮(zhèn)建設(shè)用地當(dāng)前的分布規(guī)律、擴(kuò)張方向與高程、坡度以及地貌形態(tài)的關(guān)系;瞿詩進(jìn)等[17]對長江中游經(jīng)濟(jì)帶城鎮(zhèn)建設(shè)用地在時空上的類型轉(zhuǎn)換特征進(jìn)行分析;牛叔文等[18]基于甘肅省天水市的地形復(fù)雜程度對建設(shè)用地的適宜性進(jìn)行了相關(guān)評價。以上相關(guān)研究多以地形因子的高程、坡度、地形位指數(shù)等來衡量不同土地利用類型的變化程度、城鎮(zhèn)建設(shè)擴(kuò)張強度等。
多樣性的概念與分析方法起源于信息論領(lǐng)域,20世紀(jì)90年代,以Ibáez[19-20]為首的一個西班牙研究小組首先提出了與生態(tài)學(xué)研究中生物多樣性相平行的“土壤多樣性”理論,并日益受到國際間各方的重視,其研究方法也日趨成熟[21-22]。21世紀(jì)初,張學(xué)雷等[23-24]將土壤多樣性理論與研究方法引入中國并嘗試有關(guān)研究;之后,段金龍等[25-27]將多樣性理論用于評價土地利用類型變化和水體空間分布特征,任圓圓等[28-29]探討了河南省的地形、土壤和地表水體多樣性間的格局特征;郭漩等[30]分析了伊洛河流域不同地形坡度下水土資源和地表水體、土壤、城鎮(zhèn)建設(shè)用地和交通運輸用地多樣性間的特征,從而使多樣性在土壤學(xué)、土地資源學(xué)、自然地理學(xué)等地學(xué)領(lǐng)域的研究視野變得更為開闊。鑒于此,本文將多樣性最新的測度方法運用到城鎮(zhèn)建設(shè)用地與地形要素空間分布多樣性的分析之中,運用改進(jìn)的仙農(nóng)熵指數(shù)計算研究樣區(qū)的地形多樣性、城鎮(zhèn)建設(shè)用地空間分布多樣性及不同地形上分布的城鎮(zhèn)建設(shè)用地的空間分布離散性特征,從新的視角探討二者間關(guān)系,以期為城鎮(zhèn)建設(shè)用地的合理布局和人們對城市的生態(tài)需求提供不同的數(shù)據(jù)參考。
根據(jù)自然生態(tài)環(huán)境的差異性,選取河南省北、中和南部的部分縣市作為研究樣區(qū)(圖1),樣區(qū)基本信息見表1,3個樣區(qū)之間的面積均在5 000 km2~5 300 km2之間。豫北樣區(qū)的安陽市位于河南省最北部,地處豫、晉、冀三省的交界處,地形復(fù)雜多樣性,地勢西高東低;鶴壁市位于太行山東麓向華北平原過渡地帶。豫中樣區(qū)地勢較為平坦,平原面積居多,部分地區(qū)也分布有丘陵與河川相間的地貌;豫南樣區(qū)的固始縣和商城縣屬于河南省最南部的信陽市,是南北方的過渡地帶,地勢南高北低。且豫北和豫中樣區(qū)地處暖溫帶,屬溫帶季風(fēng)氣候,豫南樣區(qū)跨亞熱帶,屬北亞熱帶向暖溫帶過渡的大陸性季風(fēng)氣候。
圖1 研究區(qū)位置示意
地形數(shù)據(jù)來自地理空間數(shù)據(jù)云平臺的DEM數(shù)據(jù)的解譯,不同研究時期城鎮(zhèn)建設(shè)用地數(shù)據(jù)來自于地理空間數(shù)據(jù)云平臺的Landsat數(shù)據(jù),影像的具體下載時間見表1。其中,不同研究區(qū)2001年和2007年的數(shù)據(jù)為Landsat7的TM傳感器獲取,2017年數(shù)據(jù)為Landsat7的OLI傳感器獲取,獲取季節(jié)在4月底和5月初主要為了減少郁蔥植被對城鎮(zhèn)建設(shè)用地的遮擋與覆蓋。3個研究時期之間相差17年,旨在研究不同區(qū)域不同時段城鎮(zhèn)建設(shè)用地在研究區(qū)內(nèi)的面積變化和空間分布離散性變化,并探討不同地形對城鎮(zhèn)建設(shè)用地空間分布格局產(chǎn)生的影響。為了便于比較,研究數(shù)據(jù)采用的坐標(biāo)系和投影均為WGS84坐標(biāo)系和UTM投影。
表1 研究樣區(qū)和數(shù)據(jù)來源
數(shù)據(jù)處理軟件有ENVI4.5,ArcGIS10.2,具體處理步驟為:首先,獲取不同樣區(qū)不同時期的遙感影像,運用ENVI4.5軟件的最大似然分類方法對研究區(qū)的數(shù)據(jù)進(jìn)行監(jiān)督分類,在分類的過程中對于不確定的區(qū)域需要結(jié)合Google Earth高清圖像進(jìn)行校正;其次運用ArcGIS10.2軟件處理土地利用分類的矢量數(shù)據(jù),提取出不同樣區(qū)不同時期的城鎮(zhèn)建設(shè)用地數(shù)據(jù);最后,生成不同研究樣區(qū)的1 km×1 km網(wǎng)格,并將其與城鎮(zhèn)建設(shè)用地數(shù)據(jù)進(jìn)行疊加并運用改進(jìn)的仙農(nóng)熵公式計算地形多樣性、城鎮(zhèn)建設(shè)用地多樣性和不同地形上的城鎮(zhèn)建設(shè)用地分布格局。
1.多樣性測度方法
在多樣性測度方法中有經(jīng)典的仙農(nóng)熵公式,為更好地表征研究區(qū)域?qū)ο蠓植嫉木鶆虺潭群涂臻g離散程度,以1 km×1 km的格網(wǎng)為尺度計算地形和城鎮(zhèn)建設(shè)用地各自空間分布多樣性指數(shù)。使用如下改進(jìn)的仙農(nóng)熵公式:
(1)
式中:在地形和土地利用空間分布多樣性中,多樣性指數(shù)Yh[21]表示研究區(qū)內(nèi)某個地形類別或不同時期的城鎮(zhèn)建設(shè)用地在樣區(qū)內(nèi)的空間分布多樣性特征,即單個地形類別或單個時期的城鎮(zhèn)建設(shè)用地在空間分布上的離散性特征。S表示研究區(qū)空間網(wǎng)格的總數(shù)目,pi表示第i個空間網(wǎng)格里某類地形(或某一年的城鎮(zhèn)建設(shè)用地)的面積占研究區(qū)內(nèi)該地形類別(或該年份的城鎮(zhèn)建設(shè)用地)總面積的比例。
多樣性指數(shù)Yh的取值區(qū)間為[0,1],當(dāng)相對豐度分布極度不均勻,即當(dāng)一個或少數(shù)幾個對象處于支配地位時,Yh取值趨于0;當(dāng)相對豐度分布均勻,即每個對象都均勻分布時,Yh趨于1。在空間分布多樣性中,單個地形類別或單個時期的城鎮(zhèn)建設(shè)用地空間分布離散性越強,Yh取值越大。當(dāng)單個地形類別或單個時期的城鎮(zhèn)建設(shè)用地集中分布在研究區(qū)內(nèi)少數(shù)幾個網(wǎng)格時,Yh取值趨于0;當(dāng)在研究區(qū)內(nèi)網(wǎng)格中均勻分布時,Yh取值趨于1。
2.關(guān)聯(lián)性分析測度方法
為探討地形和城鎮(zhèn)建設(shè)用地空間分布之間的關(guān)聯(lián)性,使用以下關(guān)聯(lián)系數(shù):
(2)
式中:A和B分別為不同樣區(qū)的地形類別和不同樣區(qū)不同時期的城鎮(zhèn)建設(shè)用地,r(A,B)[31]為某一樣區(qū)內(nèi)地形類別A和該樣區(qū)內(nèi)某一時期的城鎮(zhèn)建設(shè)用地B之間的關(guān)聯(lián)系數(shù),Yh(A)和Yh(B)分別為地形和城鎮(zhèn)建設(shè)用地的空間分布多樣性指數(shù),應(yīng)用公式(1)進(jìn)行計算,Yh(A,B)表示地形類別A和城鎮(zhèn)建設(shè)用地B公共斑塊的空間分布多樣性指數(shù),計算公式如下:
(3)
(4)
(5)
其中,p(i,j)為聯(lián)合分布概率。
在關(guān)聯(lián)性分析中,關(guān)聯(lián)系數(shù)r(A,B)的取值區(qū)間為[0,1]。如果某樣區(qū)內(nèi)的某一地形類別A和該樣區(qū)內(nèi)某一時期的城鎮(zhèn)建設(shè)用地B公共斑塊在研究區(qū)內(nèi)少數(shù)幾個網(wǎng)格(或沒有公共斑塊)時,Yh(A,B)的取值趨于0(或Yh(A,B)=0),表明二者關(guān)聯(lián)性很弱(或不相關(guān)),即該類地形對該樣區(qū)內(nèi)某一時期的城鎮(zhèn)建設(shè)用地空間分布情況的影響很弱(或無影響);如果A和B的公共斑塊分布均勻時,Yh(A,B)的取值趨近于1,表明二者的關(guān)聯(lián)性很強,即該類地形對該樣區(qū)內(nèi)某一時期的城鎮(zhèn)建設(shè)用地空間分布情況的影響很強。Yh(A,B)在[0,1]內(nèi)的取值越大,表明某樣區(qū)內(nèi)的某一地形類別A和該樣區(qū)內(nèi)某一時期的城鎮(zhèn)建設(shè)用地B的空間分布特征的相關(guān)性越強,即該類地形對該樣區(qū)內(nèi)某一時期的城鎮(zhèn)建設(shè)用地空間分布情況的影響越大。
地形空間分布多樣性表征不同類別的地形在研究區(qū)空間分布上的離散性程度。地形的空間分布多樣性指數(shù)越大,表明該類別的地形在與之對應(yīng)的研究樣區(qū)的行政區(qū)域內(nèi)分布狀況越均勻;反之,空間分布多樣性指數(shù)越小,則表明分布狀況越零散。通過對豫北、豫中和豫南樣區(qū)的DEM影像進(jìn)行解譯得到了平原、地形和丘陵三種地形類別(盆地地形所在區(qū)域均不在本次研究所選擇的樣區(qū)內(nèi),故不予考慮)。然后,將不同樣區(qū)的1 km×1 km網(wǎng)格矢量圖與相應(yīng)的地形圖疊加,帶入公式(1)得到不同研究樣區(qū)內(nèi)3種地形的空間分布多樣性值,其結(jié)果見表2。
由表2可知:就面積比例而言,3個研究樣區(qū)平原地形所占的面積比例均在53%以上,說明3個研究樣區(qū)均以平原地形為主,但豫北樣區(qū)平原地形比重最小,丘陵地形在其研究區(qū)內(nèi)也有32%的分布比重。且豫南樣區(qū)的丘陵地形和山地地形的面積比例在3個研究區(qū)中均是最低的。就空間分布多樣性值而言,總體看3個研究樣區(qū)地形空間分布多樣性指數(shù),平原最高、丘陵次之、山地最低,其多樣性的平均值分別為0.95,0.85和0.72。綜上,3個研究樣區(qū)均是以平原地形為主,其面積比例和空間分布多樣性值在相應(yīng)的研究區(qū)內(nèi)均是最大的,且豫南樣區(qū)(包括固始縣和商城縣)的山地地形的面積比例和空間分布多樣性值均是最小值。
表2 研究區(qū)不同時期城鎮(zhèn)建設(shè)用地空間分布多樣性
1.不同研究樣區(qū)城鎮(zhèn)建設(shè)用地的面積變化
從土地利用信息數(shù)據(jù)庫中提取豫北、豫中和豫南樣區(qū)在三個不同時期(2001年、2007年、2017年)內(nèi)的城鎮(zhèn)建設(shè)用地數(shù)據(jù),運用ArcGIS10.2軟件對相應(yīng)的圖層進(jìn)行疊加分析,并計算其面積和其在空間上的離散狀況。統(tǒng)計3個研究樣區(qū)三個時期的城鎮(zhèn)建設(shè)用地面積的情況,其結(jié)果見圖2。
圖2 研究區(qū)不同時期城鎮(zhèn)建設(shè)用地面積的變化
圖2是豫北、豫中和豫南樣區(qū)在不同研究時期的城鎮(zhèn)建設(shè)用地總面積的變化情況,可知:總體看,3個研究樣區(qū)不同研究時期建設(shè)用地總面積豫北樣區(qū)>豫中樣區(qū)>豫南樣區(qū),且隨著時間的推移3個樣區(qū)的城鎮(zhèn)建設(shè)用地總面積均處于上升趨勢;就各個研究樣區(qū)而言,3個樣區(qū)不同時期的面積上升幅度存在差異。其中,豫北樣區(qū)2001—2007年(時期Ⅰ)的上升幅度是14%,2007—2017年(時期Ⅱ)的上升幅度是33%;豫中樣區(qū)兩個研究時期的上升幅度48%和16%;豫南樣區(qū)兩個研究時期的上升幅度分別是2%和72%。即豫北樣區(qū)和豫南樣區(qū)均是研究時期Ⅱ的漲幅大于研究時期Ⅰ,不同的是豫南樣區(qū)內(nèi)2007年的城鎮(zhèn)建設(shè)用地總面積雖大于2001年但是沒有明顯差異,豫中樣區(qū)雖然不同時期的城鎮(zhèn)建設(shè)用地用面積差異較為明顯,但是研究時期Ⅰ城鎮(zhèn)建設(shè)用地總面積的漲幅卻大于研究時期;從時間角度看,同一樣區(qū)在不同時期的城鎮(zhèn)建設(shè)用地總面積的差異也十分明顯,且隨著時間的增長,研究樣區(qū)內(nèi)的城鎮(zhèn)建設(shè)用地總面積也在不斷增漲,總體情況呈正相關(guān)關(guān)系。
2.不同樣區(qū)城鎮(zhèn)建設(shè)用地空間分布多樣性
將豫北、豫中、豫南樣區(qū)3個研究時期的城鎮(zhèn)建設(shè)用地與基于該樣區(qū)行政邊界創(chuàng)建的1 km×1 km網(wǎng)格分別進(jìn)行疊加分析,之后將數(shù)據(jù)代入公式(1)得到在該尺度網(wǎng)格下城鎮(zhèn)建設(shè)用地在空間分布上的離散性程度??臻g分布多樣性值越接近1說明城鎮(zhèn)在研究樣區(qū)內(nèi)的分布越均勻,反之,越接近于0說明其面積較小且分布不均勻。
表3統(tǒng)計了不同研究樣區(qū)在不同研究時期的城鎮(zhèn)建設(shè)用地空間分布多樣性指數(shù),從該表可知:(1)總體上看,豫北、豫中、豫南3個樣區(qū)在同一時期的城鎮(zhèn)建設(shè)用地空間分布多樣指數(shù),豫北樣區(qū)均大于豫中、豫南樣區(qū)。(2)就不同研究樣區(qū)而言,豫北樣區(qū)三個時期的城鎮(zhèn)建設(shè)用地空間分布多樣性均大于0.84,豫中樣區(qū)均大于0.80,這說明其城鎮(zhèn)建設(shè)用地空間分布的均勻程度較好,也在一定程度上反映出去城鎮(zhèn)建設(shè)用地在研究區(qū)內(nèi)的分布和發(fā)展?fàn)顟B(tài)良好;豫南樣區(qū)三個時期的多樣性指數(shù)均大于0.5,表明其城鎮(zhèn)建設(shè)用地的空間分布的均勻程度有待提高,且從上小節(jié)可知豫南城鎮(zhèn)建設(shè)用地不同時期面積的增長幅度存在較大的差異,這都反映出豫南地區(qū)城市化建設(shè)的巨大潛力。(3)從時間角度出發(fā),同一研究樣區(qū)在不同時期的城鎮(zhèn)建設(shè)用地空間分布多樣性指數(shù)值存在差異,且時間與多樣性指數(shù)總體情況呈正相關(guān)關(guān)系,即隨著時間的增長,研究樣區(qū)內(nèi)的城鎮(zhèn)建設(shè)用地空間分布多樣性也在不斷增長。
表3 研究區(qū)不同時期城鎮(zhèn)建設(shè)用地空間分布多樣性
計算豫北、豫中和豫南樣區(qū)在平原地形、丘陵地形和山地地形下隨著時間的變化其城鎮(zhèn)建設(shè)用地在在空間分布上的離散性程度,結(jié)果見表4。圖3~圖5是2001年、2007年和2017年豫北、豫中和豫南樣區(qū)在平原、丘陵和山地地形下分布的城鎮(zhèn)建設(shè)用地情況。
由表4和圖3~圖5可知:(1)除了豫南樣區(qū)的山地地形外,其他樣區(qū)各地形條件下城鎮(zhèn)建設(shè)用地空間分布多樣性值均處于上升趨勢,且平原>丘陵>山地,這說明城鎮(zhèn)建設(shè)用地的分布在各個研究區(qū)內(nèi)趨于均勻,且平原上城鎮(zhèn)的均勻程度高于丘陵和山地。(2)就丘陵地形上分布城鎮(zhèn)建設(shè)用地多樣性而言,豫北樣區(qū)最大,豫中和豫南樣區(qū)次之,且豫北樣區(qū)的均勻程度增長較快,豫中和豫南樣區(qū)的變化較小。(3)就山地地形上的城鎮(zhèn)建設(shè)用地多樣性而言,豫北樣區(qū)大于豫中樣區(qū),豫南樣區(qū)上沒有城鎮(zhèn)建設(shè)用地的分布,這與該樣區(qū)內(nèi)山地分布的面積占比及其空間分布位置較為集中有關(guān)。
表4 不同地形上的城鎮(zhèn)建設(shè)用地空間分布多樣性
圖3 豫北樣區(qū)三個研究時期城鎮(zhèn)建設(shè)用地分布
圖4 豫中樣區(qū)三個研究時期城鎮(zhèn)建設(shè)用地分布
圖5 豫南樣區(qū)三個研究時期城鎮(zhèn)建設(shè)用地分布
豫北、豫中和豫南樣區(qū)有平原、丘陵和山地三種地形,且均以平原地形為主,其面積比例均大于53%,且就地形空間分布多樣性的平均值而言,平原最大,丘陵和山地次之。
3個研究樣區(qū)不同研究時期建設(shè)用地總面積而言,豫北樣區(qū)>豫中樣區(qū)>豫南樣區(qū),且隨著時間的推移3個樣區(qū)的城鎮(zhèn)建設(shè)用地總面積均處于上升趨勢,但3個樣區(qū)不同時期的面積上升幅度存在差異。
豫北樣區(qū)和豫中樣區(qū)三個時期的城鎮(zhèn)建設(shè)用地空間分布多樣性均大于0.80,說明城鎮(zhèn)建設(shè)用地空間分布的均勻程度較好;豫南樣區(qū)三個時期的多樣性指數(shù)均大于0.5,表明其城鎮(zhèn)建設(shè)用地的空間分布的均勻程度有待提高,且城市化建設(shè)有巨大潛力。
除了豫南樣區(qū)的山地地形外,其他樣區(qū)各地形條件下城鎮(zhèn)建設(shè)用地的分布在各個研究區(qū)內(nèi)趨于均勻,且平原上城鎮(zhèn)的均勻程度高于丘陵和山地。
綜上,3個不同樣區(qū)隨著時間的推移,其城鎮(zhèn)建設(shè)用地的在空間上的分布與發(fā)展在平原地形上最為均勻,且同時沿寄托于丘陵、山地發(fā)展的趨勢。