張 澤,丘海紅,胡寶清
(1.南寧師范大學北部灣環(huán)境演變與資源利用教育部重點實驗室,廣西 南寧 530001)(2.南寧師范大學廣西地表過程與智能模擬重點實驗室,廣西 南寧 530001)(3.南寧師范大學地理科學與規(guī)劃學院,廣西 南寧 530001)
人口與資源環(huán)境的發(fā)展是當代全球性問題,人口不斷增加,資源短缺,使得人地關(guān)系越來越緊張. 人口的空間格局反映了區(qū)域自然資源情況和經(jīng)濟水平狀況. 因此,揭示人口空間格局與自然因子的關(guān)系對緩和人地關(guān)系具有重要意義[1]. 我國人口空間格局問題在改革開放后成為了研究的熱點[2]. 目前學者們對人口的空間格局研究取得了重要進展,包括研究尺度、研究方法、影響因素和研究區(qū)選擇等. 在研究尺度上,涉及全球[3]、全國[4-5]、大區(qū)域[6-7]、省市級、縣級[8-10]. 在研究方法上,有地理加權(quán)回歸[11-12]、空間統(tǒng)計分析法[13]、洛倫茲曲線法[14]等方法. 在影響因素方面,有自然因素、經(jīng)濟因素,自然因素有高程[15]、坡度[16]、水資源[17]等,經(jīng)濟因素有人均GDP[18]、公共服務設(shè)施[19]、交通可達性[20]等. 在研究區(qū)選擇上,有高原山區(qū)[21]也有城市[2]等區(qū)域. 盡管人口空間格局研究在許多方面取得了新進展,但在研究區(qū)的選擇上有一定的局限性,影響因子探測方面需進一步的研究. 如針對喀斯特和北部灣進行單獨人口空間分布的研究頗為豐富,研究內(nèi)容分別以喀斯特地區(qū)和北部灣地區(qū)人口分布和自然因素的關(guān)系為主,而探究過渡性空間各類自然因子對人口空間格局的差異和各類自然因子疊加后對人口分布的影響強度的研究相對匱乏. 地理探測器是揭示驅(qū)動力的一種解釋方法,在應用方面沒有過多的假設(shè)條件,可以克服處理變量的局限性,很好地探測山江海過渡性空間的影響因素.
山江海過渡性空間源于山地-流域-近海綜合過渡體,始于喀斯特和北部灣區(qū)域共軛發(fā)展的簡稱,目的是構(gòu)筑山-江-海發(fā)展新格局,科學合理地利用各地區(qū)的有力優(yōu)勢,形成山江海聯(lián)動綜合發(fā)展的形式,改變現(xiàn)在發(fā)展的不均衡,縮小發(fā)展差距. 其非單一的自然地理空間,也非單一的人文地理空間,是人與自然相互作用的漸變過渡帶,表現(xiàn)出人地關(guān)系的復雜性和不確定性[22-23]. 桂西南喀斯特-北部灣海岸帶既是喀斯特巖溶治理區(qū)也是大西南沿海發(fā)達的黃金地區(qū),是推進西部大開發(fā)形成新格局的重要區(qū)域. 本研究以桂西南喀斯特-北部灣海岸帶為研究區(qū),以縣(區(qū))為研究單元,采用空間自相關(guān)的方法研究該區(qū)域的人口空間分布特征,選取影響該區(qū)域人口空間分布的主要自然影響因子,采用地理探測器模型,揭示山江海過渡性地理空間的人口空間格局與自然因子的差異及各因子疊加后對人口空間分布的影響強度,以期為該區(qū)域的生態(tài)可持續(xù)發(fā)展、生態(tài)移民安置與城鎮(zhèn)體系優(yōu)化布局提供科學參考.
桂西南喀斯特-北部灣海岸帶位于廣西壯族自治區(qū)的西南部,地理坐標:104°28′~109°56′E,20°26′~25°07′N,桂西南喀斯特地區(qū)(百色市、崇左市),共有19個縣(區(qū));北部灣經(jīng)濟區(qū)(南寧市、北海市、欽州市、防城港市),共有24個縣(區(qū)),總面積96 386 km2. 研究區(qū)地勢從西北沿云貴高原山麓向東南北部灣海岸帶整體逐漸降低. 氣候有:亞熱帶季風氣候,夏季高溫多雨,冬季溫和少雨,熱量充沛,氣溫2.8~33.4 ℃,平均氣溫20.4 ℃,降水量837.6~1 734.89 mm;還有熱帶海洋性季風性氣候,全年高溫,有明顯的旱雨兩季,氣溫5~40.4 ℃,平均氣溫23.5 ℃,降水量1 745.6~3 111.9 mm,雨量充沛,水資源豐富. 有時會有極端天氣發(fā)生,如臺風、干旱和洪澇等. 截止到2018年底,桂西南喀斯特地區(qū)人口總數(shù)672.66萬,北部灣經(jīng)濟區(qū)人口總數(shù)1 105.28萬.
隨著城市化加快推進,桂西南喀斯特-北部灣海岸帶面臨的生態(tài)環(huán)境問題日益嚴峻,主要為石漠化區(qū)水土流失、暴雨洪澇、熱帶氣旋、土壤肥力下降. 2018年,該區(qū)域43個縣(區(qū))的人口有不同程度的增長,人口數(shù)量增多將增加生態(tài)環(huán)境壓力. 此外,該區(qū)域地處中國西南邊境,對邊境地理的研究有推動作用,同時該區(qū)域是西部陸海新通道,對推動“一帶一路”建設(shè)實現(xiàn)可持續(xù)發(fā)展,以及對生態(tài)環(huán)境可持續(xù)發(fā)展都具有一定的現(xiàn)實意義.
借鑒已有人口空間分布和自然因子關(guān)系的指標選取情況[1-2],根據(jù)數(shù)據(jù)的可獲取性,選取影響該區(qū)域人口空間分布的主要自然因子. 各自然因子選取依據(jù)及數(shù)據(jù)處理方法如下:
(1)高程:高程可以影響人類生活的舒適性. 數(shù)據(jù)來源于國家基礎(chǔ)地理信息中心,分辨率30 m×30 m. 通過ArcGIS軟件分區(qū)計算各縣(區(qū))矢量斑塊覆蓋下的柵格數(shù)平均值. (2)坡度:坡度是影響城市建設(shè)成本和人類便利度的重要指標. 數(shù)據(jù)由高程數(shù)據(jù)通過ArcGIS軟件獲得柵格值,分區(qū)計算各縣(區(qū))坡度平均值. (3)土壤類型面積:適合人類居住土壤類型面積是影響人類生產(chǎn)生活的重要指標. 數(shù)據(jù)來源于《1∶100 萬中華人民共和國土壤圖》[24]. (4)地貌類型面積:適合人類生產(chǎn)生活的地貌類型是影響人類居住的重要指標. 數(shù)據(jù)來源于《中華人民共和國地貌圖集(1∶100萬)》[25]. (5)平均氣溫:氣溫對人類生產(chǎn)生活、健康等具有重要的影響. 數(shù)據(jù)來源于國家氣象科學數(shù)據(jù)中心,通過ArcGIS軟件的空間分析模塊對數(shù)據(jù)進行反比距離權(quán)重法(IDW)插值. (6)河網(wǎng)密度:人類具有居住親水性,河網(wǎng)是影響人口空間格局的重要因素. 將土地利用數(shù)據(jù)重分類并提取河網(wǎng),河網(wǎng)密度=各縣(區(qū))的河網(wǎng)總面積/對應縣(區(qū))的總面積. (7)植被覆蓋指數(shù)(NDVI):植被覆蓋的高低影響著人類氧含量和人居溫度等. 數(shù)據(jù)來源于美國地質(zhì)調(diào)查局陸地產(chǎn)品,對所獲得的陸地產(chǎn)品采用最大合成法(MVC)合成為年NDVI數(shù)據(jù),并去除異常. (8)濕潤指數(shù):人類適宜在濕潤指數(shù)適中的區(qū)域生活,太低或太高對人類的健康和生活等都有不利的影響. 數(shù)據(jù)來源于中國科學院資源環(huán)境數(shù)據(jù)中心,對獲得的數(shù)據(jù)利用ArcGIS軟件的空間分析模塊進行反比距離權(quán)重法(IDW)插值.
各縣(區(qū))的人口數(shù)據(jù)主要來源于廣西統(tǒng)計局. 本研究將常住人口數(shù)據(jù)導入到ArcGIS軟件中,進行可視化處理.
2.2.1 空間自相關(guān)
空間自相關(guān)是解析某一位置的要素觀測值與相鄰位置上觀測值的關(guān)聯(lián)程度,分為全局空間自相關(guān)和局部空間自相關(guān).
(1)全局空間自相關(guān)
全局空間自相關(guān)是某一要素的屬性值在整個空間的特征,用莫蘭指數(shù)(Global Moran’sI)表示. 取值范圍為[-1,1],莫蘭指數(shù)為負,表示空間負相關(guān),數(shù)值越小空間差異性越大;莫蘭指數(shù)為正,表示空間正相關(guān),數(shù)值越大空間相關(guān)性越好;莫蘭指數(shù)為0表示無空間相關(guān)性[1]. 公式如下:
(1)
式中,n為要素的總數(shù),S0為空間權(quán)重的集合,wi,j為要素i和j之間的空間權(quán)重,Zi是i屬性與其平均值的偏差,Zj是j屬性與其平均值的偏差.
(2)局部空間自相關(guān)
局部空間自相關(guān)是判斷要素的屬性是否有聚集和異常特征,揭示其局部空間單元的分布情況[1].公式如下:
(2)
(3)
式中,S2是要素單元xq方差,m為要素單元的總數(shù)目,wpq是要素p和q的權(quán)重.Ip>0表示該區(qū)域某要素周圍相似值在空間上集聚,反之Ip<0表示非相似在空間上集聚.
2.2.2 地理探測器
本研究采用地理探測器[26]對桂西南喀斯特-北部灣海岸帶的人口空間格局與自然因子的關(guān)系進行分析.
(1)因子探測
因子探測是探測某因子多大程度上解釋變量的空間分異.將人口密度圖層與自然因子圖層疊加,探測自然因子(具體自然因子如表1所示)間的重要性,公式如下:
表1 2018年桂西南喀斯特-北部灣海岸帶自然因子指標Table 1 Indicators of natural factors of Karst-Beibu Gulf Coastal Zone in Southwest Guangxi in 2018
(4)
(2)生態(tài)探測
生態(tài)探測用于探測比較自然因子A和B對人口密度空間分布的影響是否顯著,以F統(tǒng)計量來衡量,公式如下:
(5)
(6)
式中,nX1和nX2分別為A和B的樣本量,SSWX1和SSWX2分別是X1和X2形成層的層內(nèi)方差之和,L1和L2分別為A和B分層數(shù)目.
(3)交互作用探測
交互作用探測用于探測自然因子A、B相互作用后是否會強化或弱化人口密度的影響,主要有以下5種類型(見表2):
表2 交互探測類型Table 2 Interaction probe types
根據(jù)圖1可知,2018年桂西南喀斯特-北部灣海岸帶的縣(區(qū))人口總量和人口密度在空間分布上有一定的差異.人口總量超過100萬人有賓陽縣、橫縣、靈山縣和合浦縣,其中靈山縣人口最多,達到176.73萬人;人口總量在60~100萬人的地區(qū)主要分布在超過100萬人口總量地區(qū)的周圍;人口總量在40~60萬人的地區(qū)主要分布在該區(qū)域的中部,共有11個縣(區(qū))在此等級中;人口總量在25~40萬人的有10個縣(區(qū)),主要分布在西北地區(qū)和中部;人口總量<25萬人的地區(qū)零星地分布在該區(qū)域的周圍,憑祥市的人口最少,11.61萬人.人口密度最高的地區(qū)在該區(qū)域的東南部,分布比較集中,在超過4 000人/km2的地區(qū)中海城區(qū)的人口密度最高,達4 760人/km2;人口密度在2 000~4 000人/km2的地區(qū)主要分布在超過4 000人/km2地區(qū)的周圍,在超過 4 000 人/km2地區(qū)的西側(cè);人口密度在1 000~2 000人/km2的地區(qū)主要分布在該區(qū)域的中部,分布也比較集中;人口密度在500~1 000人/km2的地區(qū)主要分布在西北部和西南部;人口密度<500人/km2的地區(qū)主要分布在西北部和中南部,其中人口密度最小的是田林縣,為98人/km2.總體而言,桂西南喀斯特-北部灣海岸帶的人口總量分布與人口密度分布的趨勢大體相同,西北地區(qū)人口密度在1 000人/km2以下、人口總量較少,東南地區(qū)人口密度在2 000人/km2以上、人口總量較多,呈現(xiàn)出“西北少,東南多,中部居中”的特征;從西北部向東南部有逐漸增多的空間格局,呈階梯式過渡性分布.
圖1 2018年桂西南喀斯特-北部灣海岸帶人口總量和人口密度分布Fig.1 Total population and population density distribution of Karst-Beibu Gulf Coastal Zone in Southwest Guangxi in 2018
根據(jù)圖2可以分析2018年桂西南喀斯特-北部灣海岸帶人口密度的空間自相關(guān),該區(qū)域人口密度的莫蘭指數(shù)(Moran’sI)為0.409,為正值,表明不是隨機分布的,經(jīng)顯著性檢驗z值為3.922 2,p值為0.003,小于0.05,置信區(qū)間大于95%,說明該區(qū)域呈現(xiàn)出空間正相關(guān)關(guān)系.該區(qū)域人口密度的局部空間自相關(guān)不顯著的縣(區(qū))較多,共有36個;人口密度的高-高聚集區(qū)有3個縣,集中分布在該區(qū)域的東南部,有橫縣、靈山縣和浦北縣,這3個縣的人口密度也是最高的,并且該區(qū)域的地勢平坦,自然條件較好,與海域距離適中,水路交通便利;人口密度的低-低聚集區(qū)也有3個市(縣),分布在該區(qū)域的西南部和北部,分別是凌云縣、寧明縣和東興市,凌云縣是喀斯特石山區(qū),自然條件差,區(qū)域人口數(shù)量少,寧明縣和東興市均是邊境區(qū)域,遠離中心城鎮(zhèn),而且寧明縣的喀斯特山地也相對較多;人口密度的低-高聚集區(qū)只有1個區(qū),是南寧市的邕寧區(qū),該區(qū)是南寧市新區(qū),還處于剛發(fā)展的狀態(tài),遠離南寧市市中心,人口相對較少;無高-低集聚區(qū).
圖2 2018年桂西南喀斯特-北部灣海岸帶人口密度的空間自相關(guān)關(guān)系圖Fig.2 Spatial autocorrelation diagram of population density of Karst-Beibu Gulf Coastal Zone in Southwest Guangxi in 2018
不同的高程之間的人口密度差異顯著.高程在150m以下的人口密度最大,為4 760 人/km2,主要分布在該區(qū)域中南部的北部灣區(qū)域,約為高程超過600m的48倍、高程400~600m的9.5倍、高程150~200m的2.3倍.整體上,人口密度隨高程的增高逐漸降低,呈現(xiàn)出低高程聚集的特點(圖3(a)).
不同的坡度對人口密度有不同程度的影響,人口密度隨坡度的空間分布情況與隨高程的空間分布情況相似.坡度在12°以下的人口密度居高,為2 047~4 760人/km2,但坡度6°以下并不是人口密度最高區(qū).人口密度最高區(qū)在6°~12°,分布在該區(qū)域的中部,約是坡度12°~15°的2.3倍、15°~20°的4.7倍、20°以上的32倍.坡度對人口密度的空間分布具有顯著的影響(圖3(b)).
該區(qū)域有一部分是喀斯特山區(qū),土壤類型面積對人口密度的分布有一定的影響.喀斯特土壤類型面積超過400km2的主要分布在西北部,200~400km2在其周圍分布,此處人口密度較低,土壤類型面積平均為223km2.喀斯特土壤類型面積在25km2以下的主要分布在中部,并且是人口密度分布較高區(qū)域.整體上,隨著喀斯特土壤類型面積的增加人口密度逐漸減少(圖3(c)).
圖3 2018年桂西南喀斯特-北部灣海岸帶各自然因子與人口密度關(guān)系圖Fig.3 The relationship between each natural factor and population density of Karst-Beibu Gulf Coastal Zone in Southwest Guangxi in 2018
該區(qū)域的地貌類型喀斯特山地居多,其余的地貌類型還有丘陵、臺地和平原.喀斯特山地主要分布在西北和西南部,人口密度較小,約為304人/km2,沒有喀斯特山地的地區(qū)主要分布在中部和南部,人口密度相對較大,約為3 611人/km2.整體上,僅有喀斯特山地人口密度最小,人口密度隨著喀斯特地貌類型面積的增加呈現(xiàn)減小趨勢,其中影響最大的地貌類型為喀斯特山地(圖3(d)).
氣溫的變化對人口密度的空間分布有不同程度的影響.平均氣溫18~23 ℃時,人口密度最高,約 3 699 人/km2,主要分布在該區(qū)域的中部和東部,是氣溫18 ℃以下的10.3倍.氣溫18 ℃以下的地區(qū)主要分布在西北部.部分沿??h(區(qū))氣溫達到25 ℃以上,人口密度依舊很大.整體上,人口密度隨著氣溫的升高而增大(圖3(e)).
濕潤指數(shù)偏高的區(qū)域主要在北部灣的沿海地區(qū),此處的人口密度相對較大,約是3 032人/km2.該區(qū)域的中部人口密度也較大,但濕潤指數(shù)偏低,在西北區(qū)域人口密度偏小,濕潤指數(shù)也同樣較小.整體上,濕潤指數(shù)對人口密度的分布有一定的影響,會隨著濕潤指數(shù)的增大而增大(圖3(f)).
植被覆蓋情況對人口密度的空間分布有一定的影響.植被覆蓋指數(shù)0.9以上的區(qū)域均是人口密度在500人/km2以下的地區(qū),主要分布在西北部的喀斯特山區(qū),說明政府及相關(guān)部門在近幾年提出的封山育林控制人類活動取得了一定的成效.而人口密度較大的集中區(qū)是植被覆蓋指數(shù)0.3以下的區(qū)域,約是 2 917 人/km2,主要分布在該區(qū)域的中南部的城鎮(zhèn)區(qū).整體上,人口密度隨著植被覆蓋指數(shù)的增大而降低(圖3(g)).
河網(wǎng)密度對人口密度的影響是比較顯著的.該區(qū)域的河網(wǎng)密度呈現(xiàn)出兩極分化的情況,西北部除右江區(qū)以外的河網(wǎng)密度較小,同時人口密度也較小,約是941人/km2,東南部除興寧區(qū)以外的河網(wǎng)密度較大,人口密度也較大,約是3 011人/km2.整體上,河網(wǎng)密度越大,人口密度越大(圖3(h)).
地理探測器可以很清晰系統(tǒng)地分析自然因子對人口分布的影響,各自然因子對人口密度的作用強度不同,同一自然因子對不同地區(qū)的人口密度的影響強度也存在差異,不同的自然因子疊加后影響強度有增減作用,將8個自然因子分等級后離散化處理進行地理探測.
因子探測器揭示各因子對人口密度的影響力,通過計算自然因子的PD,H值,提取8個自然因子對人口密度的影響,排序是:高程>河網(wǎng)密度>坡度>平均氣溫>濕潤指數(shù)>植被覆蓋指數(shù)>地貌類型面積>土壤類型面積. 根據(jù)探測結(jié)果發(fā)現(xiàn),高程的影響力最大,PD,H值為0.768且解釋力在75%以上,所以高程是主要的自然因子;河網(wǎng)密度和坡度的PD,H值分別為0.644和0.541,解釋力在50%以上;平均氣溫和濕潤指數(shù)的PD,H值分別為0.337和0.301,解釋力均超過了30%;植被覆蓋指數(shù)、地貌類型面積和土壤類型面積的PD,H值分別為0.281、0.270和0.212,解釋力均超過20%,且土壤類型面積對人口密度的空間分布解釋力是最小的.
生態(tài)探測器反應自然因子對人口密度空間分布的影響是否顯著,每兩個自然因子之間是否有顯著差異.由表3可知:高程與河網(wǎng)密度和濕潤指數(shù)對人口密度空間分布有顯著差異,與植被覆蓋指數(shù)、土壤類型面積、地貌類型面積、坡度、平均氣溫對人口密度空間分布無顯著差異;坡度與土壤類型面積、地貌類型面積、濕潤指數(shù)、植被覆蓋指數(shù)對人口密度空間分布有顯著差異,與平均氣溫和河網(wǎng)密度對人口密度空間分布無顯著差異;土壤類型面積與地貌類型面積、河網(wǎng)密度、坡度、平均氣溫、濕潤指數(shù)和植被覆蓋指數(shù)對人口密度空間分布有顯著差異;地貌類型面積與河網(wǎng)密度、平均氣溫、濕潤指數(shù)和植被覆蓋指數(shù)對人口密度空間分布有顯著差異;平均氣溫與濕潤指數(shù)和植被覆蓋指數(shù)對人口密度空間分布有顯著差異,與河網(wǎng)密度對人口密度空間分布無顯著差異;濕潤指數(shù)與植被覆蓋指數(shù)對人口密度空間分布有顯著差異,與河網(wǎng)密度對人口密度空間分布無顯著差異;植被覆蓋指數(shù)與河網(wǎng)密度對人口密度空間分布有顯著差異.
表3 2018年桂西南喀斯特-北部灣海岸帶各自然探測因子的顯著性檢驗Table 3 The significance test of each natural detection factor of Karst-Beibu Gulf Coastal Zone in Southwest Guangxi in 2018
交互作用探測器反應自然因子兩兩交互作用下對人口密度空間分布的解釋力是增強、減弱或相互獨立. 利用該探測器探測驅(qū)動因子之間和人口密度空間分布的相互關(guān)系,得到表4.
表4 2018年桂西南喀斯特-北部灣海岸帶各自然探測因子交互作用情況Table 4 The interaction of various natural detection factors of Karst-Beibu Gulf Coastal Zone in Southwest Guangxi in 2018
研究發(fā)現(xiàn)自然因子對人口密度空間分布的交互作用結(jié)果只有雙線性增強和非線性增強,不存在減弱作用和相互獨立. 交互作用探測結(jié)果顯示,自然因子高程X1、土壤類型面積X3、地貌類型面積X4、植被覆蓋指數(shù)X7與其他自然因子均呈雙線性增強;坡度X2與平均氣溫X5、濕潤指數(shù)X6呈非線性增強;平均氣溫X5與濕潤指數(shù)X6、河網(wǎng)密度X8呈非線性增強. 其中高程因子與其他自然因子兩兩疊加后的作用較強,尤其與坡度因子的交互作用最強,解釋力超過90%,主要是該區(qū)域山地較多,海拔落差大,其他因子不能避免受高程的影響. 土壤類型面積和地貌類型面積與其他因子疊加時,解釋力一般,因為該區(qū)域地處亞熱帶,氣候條件較好,水量充沛,具備適宜人類居住的自然條件,所以土壤類型面積和地貌類型面積在單獨作用時對人口密度空間分布的解釋力不強. 各類自然因子交互探測結(jié)果顯示,當兩類自然因子適合人類居住時,交互作用結(jié)果表現(xiàn)為增強. 自然因子對人口密度空間分布的影響不是相互獨立的,任何因子的改變都可能影響人口空間分布的變化.
(1)桂西南喀斯特-北部灣海岸帶的人口分布不平衡,西北地區(qū)人口密度小、人口總量較少,東南地區(qū)人口密度大、人口總量較多,呈現(xiàn)出“西北少,東南多,中部居中”的特征;從西北部向東南部呈逐漸增多的狀態(tài),人口分布呈現(xiàn)了顯著正相關(guān)和一定聚集性,而且高值區(qū)和低值區(qū)不是隨機分布.
(2)桂西南喀斯特-北部灣海岸帶的人口空間格局隨自然因子的不同而不同. 整體上,人口密度隨適宜人類居住的自然因子的增加而增加. 高程較低、坡度、濕潤指數(shù)和植被覆蓋指數(shù)較小、土壤類型多樣且溶蝕性較弱的區(qū)域人口密度較高;平均氣溫較高、河網(wǎng)密度較大的區(qū)域人口密度較低.
(3)單因子對人口密度空間分布影響最大的是高程,影響最小的是土壤類型面積,所有探測因子解釋力超過20%. 任意兩因子疊加作用后都能增強對人口空間分布的影響,強化方式為雙線性增強和非線性增強. 高程因子與其他自然因子兩兩疊加后的作用較強,尤其與坡度因子的交互作用最強,解釋力超過90%. 各類自然因子交互探測結(jié)果顯示,適合人類居住的兩類自然因子交互作用結(jié)果表現(xiàn)為增強.