韓瑞丹,張 麗,鄭 藝,王 恒,張 靜

1 山東農(nóng)業(yè)大學信息科學與工程學院, 泰安 271018 2 海南省地球觀測重點實驗室, 三亞 572029 3 中國科學院遙感與數(shù)字地球研究所數(shù)字地球重點實驗室, 北京 100094

曼谷城市擴張生態(tài)環(huán)境效應

韓瑞丹1,2,張 麗2,3,*,鄭 藝3,王 恒3,張 靜3

1 山東農(nóng)業(yè)大學信息科學與工程學院, 泰安 271018 2 海南省地球觀測重點實驗室, 三亞 572029 3 中國科學院遙感與數(shù)字地球研究所數(shù)字地球重點實驗室, 北京 100094

利用1990、2000、2010和2015年4期Landsat數(shù)據(jù)和2000—2015年的MODIS NDVI數(shù)據(jù),開展了泰國曼谷地區(qū)近25年來城市擴張分析,包括城市建設用地擴張時空變化和擴張模式分析、擴張的生態(tài)環(huán)境效應分析以及擴張驅動力分析。結果表明：(1)自1990年以來,曼谷城市建設用地面積持續(xù)增長,經(jīng)歷了先快后慢再快的增長過程。城市擴張主要沿東北方向,且在距離城市中心14—20km的范圍內(nèi)擴張最為明顯,擴張模式以邊緣式和填充式為主。(2)在城市化進程中,植被和熱島對城市擴張分別有不同的響應。植被綠度整體呈下降趨勢,下降明顯的地區(qū)位于距城市中心10—20km的范圍內(nèi),即城市建設用地明顯擴張的區(qū)域。熱島分布與城市建成區(qū)分布具有空間一致性,自1990s以來,熱島效應整體增強,但局部變異越來越不明顯。(3)曼谷城市擴張受到自然地理條件、社會經(jīng)濟、城市布局等多重因素的影響。

城市擴張；植被；熱島；驅動因素

20世紀以來,世界城市化進入全面發(fā)展階段。伴隨著經(jīng)濟建設的快速發(fā)展,城市人口持續(xù)增長,城市建筑群規(guī)模不斷擴大,城市氣候、環(huán)境問題日益突出。相比于內(nèi)陸城市,港口城市因其便捷的交通條件發(fā)展更加迅速。港城分布在生態(tài)環(huán)境脆弱的海岸帶區(qū)域,其擴張給海岸帶和城市內(nèi)部生態(tài)環(huán)境帶來沉重負擔。認識港口城市發(fā)展特征、港口城市土地資源和生態(tài)環(huán)境狀況,厘清港口城市擴張對海岸帶陸地生態(tài)環(huán)境的影響,對于探究人類社會經(jīng)濟活動對海岸帶生態(tài)環(huán)境影響、制定合理港口城市發(fā)展及海岸帶生態(tài)環(huán)境保護方案具有重要意義。

城市化的快速推進改變下墊面的物質(zhì)組成和特性,影響地表生態(tài)要素的生理生態(tài)狀況、結構組成、空間格局等,對資源消耗和生態(tài)環(huán)境質(zhì)量等產(chǎn)生重大影響。目前,研究者多以各國首都城市、一線城市、重要經(jīng)濟圈等為研究區(qū)[1- 2],以Landsat、SPOT、HJ、航空影像、MODIS等為遙感數(shù)據(jù)源[3],集中于研究城市土地利用/覆被變化及驅動力[4]、城市擴張形態(tài)和模式[5]、城市擴張模擬預測等。對于城市擴張的生態(tài)影響研究,多集中于城市植被變化[6- 7]、熱島效應[8]、大氣污染[9]、水體污染[10]、土壤碳循環(huán)[11- 12]、城市內(nèi)澇等。在港口城市擴張生態(tài)環(huán)境效應研究方面,大多學者沿用內(nèi)陸城市的研究方法,集中在對城市土地利用結構時空變化、空氣污染、能源消耗等方面,也有學者對沿岸近海區(qū)的生態(tài)環(huán)境展開研究,如岸線變遷、岸堤侵蝕、圍海造田等。植被作為陸地生態(tài)系統(tǒng)的主體,對于地表環(huán)境、氣候變化具有強烈的敏感性,而城市熱島是城市生態(tài)環(huán)境失調(diào)引起的最顯著的城市環(huán)境問題之一,因此,城市植被變化和城市熱島問題是所有城市化生態(tài)環(huán)境效應研究的兩大熱點。

曼谷是泰國首都和中南半島最大的城市,同時又是“21世紀海上絲綢之路”建設的重要節(jié)點港口城市。20世紀50年代,隨著工業(yè)的迅速發(fā)展,曼谷經(jīng)濟有所提高,城市化逐漸興起。1990年至今的這段時間是曼谷城市發(fā)展的重要階段[13]。曼谷城市容納著泰國一半以上的城市人口,城市首位度畸高、虛假城市化的現(xiàn)象世界罕見[14]。截止到2015年,曼谷城市人口達到927.0萬人,相比于1990年增長了約57%,且預測至2030年將達到1152.8萬人[15]。人口的快速城市化引起了土地的城市化,從而導致曼谷城市綠地減少、地表溫度升高、空氣污染、水污染等一系列生態(tài)問題突出[16]。在中國,已有眾多學者對國內(nèi)港口城市(如大連、天津、江蘇、上海、福州、廣州等)展開了城市擴張及生態(tài)環(huán)境影響等方面的研究[17],然而在泰國曼谷,這方面的研究卻相對較少。為了全面了解曼谷快速城市化過程中城市建設用地擴張情況和城市擴張對城市環(huán)境的影響情況,本文利用多時相Landsat和MODIS遙感數(shù)據(jù),開展曼谷港口城市土地覆蓋/利用類型變化監(jiān)測和城市建設用地變化監(jiān)測,探討城市化過程中植被和近地表溫度多環(huán)境要素的時空變化規(guī)律,分析城市化的驅動因素。本研究從泰國城市擴張規(guī)律、模式和城市化對植被、熱島的影響著手,探究曼谷城市發(fā)展的合理性,以期為我國城市規(guī)劃和發(fā)展、城市生態(tài)文明建設提供借鑒或參考。

1 研究區(qū)與研究方法

1.1 研究區(qū)概況

曼谷是泰國唯一的府級直轄市,位于湄南河下游,南部瀕臨泰國灣,面向太平洋,介于13°29′—13°58′N、100°19′—100°57′E之間,是泰國首都和中南半島最大城市,也是全球最受歡迎的旅游城市之一,具有“天使之城”的美譽。該地區(qū)地勢平坦,氣候濕熱,屬于熱帶季風氣候,盛行西南風,年平均氣溫為27.5℃,年降水量約為1500mm。曼谷既是中國-東南半島經(jīng)濟走廊的重要節(jié)點城市,同時又擁有海上絲綢之路的重要港口——曼谷港。曼谷港屬于河港,位于湄南河下游東岸,距離河口約25.75km。

18世紀末19世紀初,為了發(fā)展對外貿(mào)易,泰國將首都遷至曼谷,曼谷港口隨之發(fā)展起來。泰國經(jīng)濟屬外向型經(jīng)濟,對外貿(mào)易在國民經(jīng)濟中占據(jù)重要地位,而曼谷港作為泰國最大的港口,承擔著泰國90%的外貿(mào)活動,對曼谷城市乃至泰國發(fā)展具有主導作用,同時也是中國連通東南亞和南太平洋的重要交通樞紐。

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

1.2.1 城市建設用地信息提取

為提取曼谷城市建設用地專題信息、監(jiān)測研究區(qū)土地利用/覆蓋類型變化情況,本文選取1990、2000、2010、2015年4個時期的Landsat5 TM和Landsat8 OLI數(shù)據(jù),經(jīng)過幾何校正、配準、鑲嵌、裁剪等預處理,采用最大似然監(jiān)督分類法對研究區(qū)的土地類型按照建設用地、綠地、水體、裸地4種地類進行分類,以提取城市建設用地專題信息。通過生成隨機點與Google Earth影像對比驗證,4個時期的總體分類精度均在80%以上。

1.2.2 城市建設用地擴張分析

基于以上獲取的4期城市建設用地專題信息數(shù)據(jù),以2015年的建設用地重心為原點,分別采用象限方位法和緩沖帶分析法將研究區(qū)分為4個象限和20個2km間隔的圓形緩沖區(qū),引入建設用地比重、擴張強度指數(shù)(AI)和年均增長率指數(shù)(AGR)定量分析各象限和各緩沖區(qū)內(nèi)的城市建設用地空間分布特征和時空變化特征。其中建設用地比重等于一定范圍內(nèi)建設用地面積占區(qū)域土地總面積的比例,AI用來對比不同時期同一地區(qū)城市建設用地擴張的強度,AGR用來對比不同地區(qū)同一時間城市建設用地的擴張速率,其計算公式分別如下：

(1)

(2)

式中,Astar、Aend分別為研究單元初期、末期的建設用地面積,d為時間間隔。

景觀擴張指數(shù)(LEI)可定量識別城市空間擴張的類型,本研究采用LEI表示城市新增建設用地斑塊的類型,該指數(shù)基于新增建設用地斑塊生成的緩沖區(qū)進行定義[18]：

(3)

式中,LEI取值范圍是[0—100],S0為新增建設用地斑塊生成的30m緩沖區(qū)內(nèi)原有建設用地斑塊面積,Sv為緩沖區(qū)內(nèi)除去原有建成區(qū)斑塊剩余的面積。本研究將LEI=0的斑塊定義為飛地式,即新增建設用地斑塊與原有建設用地斑塊分離；將LEI∈(0,50)的斑塊定義為邊緣式,即新增建設用地斑塊沿原有建設用地斑塊的邊緣擴展；將LEI∈[50,100]的斑塊定義為填充式,即新增建設用地斑塊填充在原有建設用地斑塊內(nèi)部[19]。

1.2.3 城市擴張生態(tài)環(huán)境效應分析

(1)植被時空變化監(jiān)測

本研究選取歸一化植被指數(shù)(NDVI)表征綠地植被的生長狀況和動態(tài)變化。NDVI數(shù)據(jù)來自美國航空航天局(NASA)提供的16d最大值合成植被指數(shù)產(chǎn)品MOD13Q1,該數(shù)據(jù)空間分辨率為250m,時相選取包括2000—2015年逐期數(shù)據(jù)。首先將每年16d的NDVI數(shù)據(jù)進行最大值合成,得到逐年NDVI年最大值,然后利用一元線性回歸方程分析植被變化趨勢：

(4)

式中，n為監(jiān)測時間序列的長度,MNDVI為第i年NDVI的最大值；θslope為像元NDVI的回歸斜率值。當θslope>0時,該區(qū)域的植被呈增長或改善趨勢,反之呈減少或退化趨勢,當θslope=0時,說明該區(qū)域NDVI沒有變化。

在植被變化趨勢分析的基礎上,本文利用緩沖分析法研究植被沿城鄉(xiāng)梯度的過度特征。即以2015年建設用地重心為圓心,建立20個以2km為間隔的圓形緩沖區(qū),統(tǒng)計分析各環(huán)形緩沖區(qū)內(nèi)NDVI的平均變化趨勢。

(2)熱島效應監(jiān)測

基于Landsat熱紅外影像,本文采用覃志豪的單窗口算法[20- 22]反演地表溫度(LST)。數(shù)據(jù)預處理包括輻射定標、大氣校正、幾何校正、配準、重采樣等?；趩未翱谒惴ǚ囱莸乇頊囟鹊闹饕襟E包括：(1)亮溫計算；(2)地表比輻射率計算；(3)近地表溫度和大氣參數(shù)計算；(4)地表真實溫度反演。

基于Landsat反演得到的地表溫度數(shù)據(jù),引入城市熱力場變異指數(shù)(TFVI)定量分析熱島效應。TFVI定義為某點LST與研究區(qū)域平均LST的差值與研究區(qū)域平均LST之比[23],計算公式如下：

(5)

式中,TFVI為城市熱力場變異指數(shù)；T為城市某點的LST,TMEAN為城市研究區(qū)域的平均LST。為了直觀地描述城市熱場變化,本文采用閾值法先將城市熱力場分為非熱島(TFVI≤0)和熱島(TFVI>0),又將熱島細分為弱熱島(TFVI∈(0,0.15])、強熱島(TFVI∈(0.15,0.3])和極強熱島(TFVI>0.3),做出城市熱島效應分級示意圖。為了方便城市熱島效應的比較,引入城市熱島比例指數(shù)來表征城市熱島的影響程度,城市熱島比例指數(shù)為城市熱島面積與研究區(qū)面積的比值,其值越大,說明熱島效應越嚴重。

2 結果與分析

2.1 城市建設用地時空變化及擴張模式分析

2.1.1 時空變化分析

對4個時期的城市建設用地影像進行疊加,從不同時相建設用地空間分布情況(圖1,表1)看,曼谷25年間建設用地發(fā)生了很大變化,擴張現(xiàn)象十分明顯。從1990—2015年曼谷城市建設用地面積持續(xù)增長,面積比例從25.96%增長至52.39%,增長了1倍多。就年均擴張程度AI來說,1990—2000年的AI為18.84km2,2000—2010年的AI為13.26 km2,2010—2015年的AI為18.62 km2,即在1990—2000年間城市擴張強度最大,其次為2010—2015年,而2000—2010年擴張程度相比其他兩個階段較弱,即自1990年至2015年以來,曼谷城市擴張呈現(xiàn)先快后慢而后又增快的特點。

圖1 曼谷1990—2015年建設用地時空變化Fig.1 Spatio-temporal changes of built-up land in 1990, 2000, 2010 and 2015 in Bangkok

結合4個象限建設用地面積、比重、年均擴張程度AI和年均擴張速率AGR,曼谷建設用地主要分布在第一、三象限,即東北和西南方向,但建設用地比重最大的地方位于第二、三象限,即西部地區(qū)(圖2)。從建設用地的年均擴張強度AI和年均增長率AGR來看,各象限的城市建設用地均在擴張,但存在明顯的時空異質(zhì)性。1990—2000年各象限的年均擴張速率AGR均明顯大于2000—2010和2010—2015兩個階段,可見曼谷于20世紀末經(jīng)歷了快速的城市化進程；3個時間段內(nèi)的AI、ARG均在第一象限為最大值,第二象限為最小值,說明曼谷城市建成區(qū)主要沿東北地區(qū)擴展,而沿西北方向擴張相對不明顯(圖3)。

表1 曼谷不同時期建設用地面積、比重和擴張強度指數(shù)AI

圖2 各象限建設用地面積和建設用地面積比重Fig.2 Urban built-up land areas and built-up land areas proportion in four sectors

圖3 各象限建設用地擴張強度AI和年均擴張速率AGRFig.3 Annual increase and annual growth rate of built-up land areas in four sectors

從緩沖區(qū)分析結果(圖4)上看,城市建成區(qū)比重隨著距城市中心距離變遠而變小。自2010年起,在距離城市中心10km以內(nèi)的地區(qū),建成區(qū)比重達到80%以上。由于城市中心2km范圍內(nèi)存在公園綠地,所以城市中心的建成區(qū)比重并未達到最大值,而是在距離城市中心2—4km的范圍達到最大(圖4)。從空間上看,城市擴張集中在距離城市中心10—26km的地方,其中14—20km處城市擴張最為顯著。從時間上看,在距離市中心12km以內(nèi)的地區(qū),AI2010—2015

圖4 城市建設用地沿城鄉(xiāng)梯度變化特征Fig.4 Spatio-temporal change characteristics of urban built-up land areas along urban-rural gradient

2.1.2 擴張模式分析

從城市擴張模式統(tǒng)計結果上看,3個時段的城市擴張模式均以邊緣式和填充式為主(圖5)。填充式的擴張模式主要分布在近城市地區(qū),大約在距城市中心12km以內(nèi)的地方；而邊緣式的擴張模式分布在距離城市中心較遠的地方。對比3個時期填充式擴張模式峰值出現(xiàn)的位置,發(fā)現(xiàn)其明顯向城市外邊緣移動,且在2010—2015年,填充式擴張模式在距離城市中心14—20km的地方明顯增多(圖5)。隨著城市擴張進程的推進,填充式的擴張模式無論在斑塊數(shù)量或在面積上均呈現(xiàn)增長趨勢,而邊緣式的擴張模式則呈現(xiàn)下降趨勢(表2),說明城市建設用地的聚集度正在提高。

以填充式為主的城市擴張模式有利也有弊。其一方面可以提高土地利用強度,避免土地資源浪費,另一方面卻會加大城市內(nèi)部的壓力,如開放空間不足、交通堵塞、環(huán)境污染等。Estoque[13]對曼谷大都會(包括本文的曼谷及其周邊地區(qū))1988—2009年的城市擴張研究結果顯示,曼谷大都會的城市擴張以邊緣式為主導,而本文的曼谷城市擴張以邊緣式和填充式為主。這充分說明,以填充式為主的擴張模式主要存在于曼谷內(nèi)部。結合曼谷交通擁擠、環(huán)境污染等典型的“城市病”問題,這種擴張模式勢必會吸引更多人向城市內(nèi)部擁擠,從而加大“城市病”的病況,對于緩解曼谷首位度過高的現(xiàn)象十分不利。因此,如何在充分利用城市土地和緩解城市內(nèi)部壓力之間尋求平衡點,應成為曼谷城市規(guī)劃需考慮的問題。

圖5 基于緩沖區(qū)的三種擴張模式面積統(tǒng)計Fig.5 Areas of three expansion patterns based on buffer zones

擴張模式Expansionpattern1990—20002000—20102010—2015斑塊數(shù)量比例Proportionofpatchnumber斑塊面積比例Proportionofpatcharea斑塊數(shù)量比例Proportionofpatchnumber斑塊面積比例Proportionofpatcharea斑塊數(shù)量比例Proportionofpatchnumber斑塊面積比例Proportionofpatcharea填充式Infilling47.3834.0353.9145.3561.7447.00邊緣式Edge-expansion37.2857.8736.7550.5632.5347.95飛地式Outlying15.348.109.344.095.735.05

斑塊數(shù)量比例指各種擴張模式類型的斑塊數(shù)量占新增斑塊總數(shù)量的比例；斑塊面積比例指各種擴張模式類型的斑塊面積占新增斑塊總面積的比例

從本研究和Estoque的研究結果的來看,無論是曼谷或大曼谷,其城市新增建設用地斑塊都是沿舊城區(qū)逐漸向外擴張,可見曼谷的城市發(fā)展仍屬于由中心向外蔓延的“攤大餅式”擴張模式。中國上海港之于中國的地位與曼谷港之于泰國的地位相當,上海城市擴張與曼谷城市擴張存在諸多相似之處。兩個城市地勢均比較平坦,城市邊緣的土地類型均以農(nóng)田為主,城市擴張均沿背離海岸帶的方向發(fā)展。但不同的是,上海城市發(fā)展呈“衛(wèi)星城”式發(fā)展,這種發(fā)展模式可有效緩解城市中心壓力,同時帶動周邊地區(qū)共同發(fā)展[24]。而對于曼谷來說,其城市首位度畸高,周邊地區(qū)發(fā)展落后,區(qū)域經(jīng)濟水平差異十分明顯。在這一點上,曼谷或可借鑒上?！靶l(wèi)星城”式的發(fā)展模式,或在落后地區(qū)創(chuàng)建新的城市中心,以分散城市人口、工業(yè)過分集中的狀況,從而緩解曼谷城市中心的交通壓力、環(huán)境壓力、人口壓力,同時帶動曼谷周邊區(qū)域的共同發(fā)展[14]。

2.2 城市擴張生態(tài)環(huán)境效應分析

2.2.1 城市擴張對植被綠度的影響

從NDVI變化趨勢圖(圖6)中統(tǒng)計得到,自2000—2015年NDVI顯著降低(P≤0.05)的面積占曼谷全區(qū)總面積的22.59%,NDVI顯著增加(P≤0.05)的面積占曼谷全區(qū)總面積的11.56%。從NDVI變化趨勢緩沖分析結果上看(圖7),植被NDVI呈下降趨勢的地區(qū)主要分布在距離城市中心24km以內(nèi)的地方,呈增長趨勢的地區(qū)主要分布在距離城市中心24km以外的地方。在距離城市中心10—20km的范圍內(nèi),植被綠度下降最為顯著,且與AI具有顯著的負相關關系,可說明,在近城市中心24km以內(nèi)的地區(qū),城市擴張是導致植被NDVI降低的主要因素。而在距離城市中心24km以外的地方,NDVI變化趨勢與AI無明顯關系,說明遠離城市中心的植被還未受到城市擴張的影響。

圖6 曼谷2000—2015年NDVI變化趨勢Fig.6 NDVI trends from 2000 to 2015 in Bangkok

圖7 曼谷城市擴張中植被綠度與AI、AGR的變化情況Fig.7 Variations of NDVI slope, AI, and AGR along the distance to city center

2.2.2 城市擴張熱島效應

對比4個時相的熱島空間分布圖(圖8)和曼谷土地利用/覆被圖(圖9),可以看出,1990年曼谷熱島主要分布在建設用地區(qū)域和東部的裸地地區(qū),其他3個時相的熱島主要分布在建設用地區(qū)域,非熱島主要分布在水體和綠地地區(qū)。4個時相的較強熱島和極強熱島主要分布在曼谷北郊的廊曼國際機場和曼谷南部的曼谷港口以及湄南河兩側的高密度建設用地區(qū)域?？梢?曼谷的城市熱島分布與建設用地分布在空間上具有一致性。與其他3個時相不同,1990年的熱島除了分布在建成區(qū)外,還主要分布在東北地區(qū)。這是由于早在1990年時,曼谷東北部地區(qū)以裸地為主,裸地的地表反射率較高,地表溫度較高,故出現(xiàn)了熱島現(xiàn)象。1990年東北部大面積的裸地同時導致曼谷全區(qū)熱島比例指數(shù)偏大,而到90年代后期,泰國政府為促進制造業(yè)和服務業(yè)的發(fā)展,決定加強農(nóng)業(yè)投入[25],故在曼谷東北地區(qū)可發(fā)現(xiàn),1990年時的大片裸地到2000年時已變成有植被覆蓋的耕地,所以2000年時的熱島范圍明顯降低、熱島比例指數(shù)明顯下降,同時使城市內(nèi)部高密度建成區(qū)部分的較強熱島現(xiàn)象更加明顯。

圖8 4個時期曼谷城市熱島空間分布Fig.8 Spatial distribution of urban heat island of Bangkok in four periods

圖9 4個時期曼谷土地利用/覆被圖Fig.9 Land cover maps in four periods in Bangkok

聚焦于2000、2010、2015年的熱島變化(表3)可發(fā)現(xiàn),曼谷城市熱島范圍增大,熱島效應整體增強,熱島比例指數(shù)從2000年的46.43%上升至2015年53.15%,其中弱熱島比例指數(shù)與熱島比例指數(shù)變化趨勢一致,呈逐漸上升趨勢,而較強熱島和極強熱島呈下降趨勢。說明在城市化進程中,城市熱島效應整體增強,但局部變異并不明顯。Arifwidodo[16]的研究成果與本文的研究成果具有相似之處,即均發(fā)現(xiàn)曼谷城市極強熱島主要分布在高密度的建成區(qū),且曼谷的熱島強度要比上海、舊金山等地高,海風、降水、云、土地覆蓋類型、人口密度是影響城市熱島的主要因素。

2.3 城市擴張驅動力分析

曼谷城市建設用地擴張動態(tài)分析表明,自1990年以來,曼谷建設用地一直處于擴張階段,各象限的城市建設用地擴張具有明顯的時空差異性,其原因主要取決于曼谷的自然地理位置、社會經(jīng)濟條件和城市布局等。

表3 不同等級城市熱島比例指數(shù)

(1)自然地理位置 曼谷位于湄南河三角洲地區(qū),地勢低洼。貫穿于曼谷的湄南河將市區(qū)分成東西兩個部分,西岸如今被視為舊市區(qū),城市發(fā)展緩慢；而大皇宮、曼谷港、眾多工業(yè)園區(qū)、機場等均位于東岸,因此曼谷東部地區(qū)的城市化發(fā)展要快于西部。根據(jù)Google Earth高分辨率影像和本研究得到多時相土地利用類型空間分布圖(圖9),曼谷東北部土地相對遼闊,并且以農(nóng)田為主,成為城市擴張的主要方向。而曼谷南部瀕臨泰國灣,地勢較低,陸地面積有限,城市擴展并不明顯。

(2)社會經(jīng)濟條件 為探討社會經(jīng)濟條件對城市擴張的驅動作用,本文選取城市人口、GDP和人均GDP三項社會經(jīng)濟指標,與城市建成區(qū)面積進行相關分析(表3)。自20世紀90年代以來,曼谷城市人口從588.8萬上升到927.0萬,增長了57.44%。曼谷建設用地面積和人口呈高度正相關,相關系數(shù)為0.97(R2=0.95),表明人口增長促進了城市建設用地的增長。根據(jù)曼谷經(jīng)濟約占泰國經(jīng)濟總量的40%這一比例[25]估算曼谷GDP,并由此計算曼谷人均GDP,結果顯示,曼谷GDP從1990年至2014年增長了374.34%,人均GDP增長了201.32%。曼谷城市建設用地與GDP和人均GDP均呈正相關關系,相關系數(shù)分別為0.96(R2=0.92)和0.95(R2=0.90)?？梢?經(jīng)濟增長是城市建設用地擴張的驅動力之一。

表4 曼谷建設用地面積與驅動因子的相關性

數(shù)據(jù)來源：城市人口數(shù)據(jù)來自United Nations,GDP數(shù)據(jù)來自world Bank;因2015年GDP數(shù)據(jù)尚未公布,故用2014年代替;*因2015年GDP數(shù)據(jù)尚未公布，故用2014年代替;—:無數(shù)據(jù)

(3)城市布局 曼谷城市建設以大皇宮為中心,最接近中心的是寺廟和政府機構,向外為商業(yè)圈,再向外是居民住宅區(qū)。城市內(nèi)圈由于分布著大皇宮、寺廟等著名的旅游景點和政府機構等,經(jīng)長期發(fā)展建設用地密集,所以在近幾十年來,城市中心的建設用地擴張并不明顯。此外,被冠于“東方威尼斯”稱號的曼谷城市水上交通十分發(fā)達,伴隨著1999年運營的高架天車線路、2004年運營的地鐵線路、2010年運營的機場鐵路等交通線路的建設,曼谷城市土地利用模式和城市空間發(fā)生了轉變[26]。從本研究新增斑塊的空間分布特征可以看出,曼谷城市擴張具有分布在主要交通干線兩側的特征,尤其在中東部地區(qū)十分明顯(圖10)。可見,交通干線是影響城市擴張的又一驅動因素。

圖10 交通對城市擴張影響案例Fig.10 Case of the influence of traffic on urban expansion

3 結論與探討

港口是城市對外開放的門戶,對城市發(fā)展具有主導作用。港口在帶動城市經(jīng)濟發(fā)展的同時,也給城市生態(tài)環(huán)境帶來一些負面影響,如何協(xié)調(diào)港城經(jīng)濟與生態(tài)環(huán)境和諧可持續(xù)發(fā)展,已成為城市規(guī)劃考慮的重要問題。作為中國-東盟自由貿(mào)易區(qū)、中國-中南半島經(jīng)濟走廊、21世紀海上絲綢之路重要節(jié)點城市,曼谷具有不可替代的國際戰(zhàn)略地位和交通樞紐作用,位于其境內(nèi)的曼谷港作為泰國最大的港口,對曼谷城市發(fā)展起到巨大的推動作用。本文利用多時相Landsat和MODIS遙感數(shù)據(jù)對曼谷城市擴張生態(tài)環(huán)境效應進行了研究,結果表明自20世紀90年代以來,曼谷城市建設用地面積持續(xù)增長,經(jīng)歷了先快速后緩慢而后又快速增長的過程。曼谷城市擴張主要沿東北方向和城市外緣方向,城市擴張模式以邊緣式和填充式為主。在城市化進程中,植被綠度整體呈下降趨勢,退化明顯的地區(qū)與建成區(qū)擴張顯著的地區(qū)相吻合；熱島分布與建成區(qū)分布也具有空間一致性,呈整體增強趨勢。通過本研究分析得出,曼谷的社會經(jīng)濟、自然地理條件、城市布局等對城市擴張具有重要影響。

與曼谷相似,我國港口城市(如廣州、上海、廈門、煙臺、大連、天津等)建設用地擴張主要以占用耕地為代價,導致農(nóng)業(yè)用地喪失。同時城市化加大了糧食需求,又使部分耕地受到污染,從而激化了人地矛盾[24,27- 29]。建設用地擴張是人類活動對自然土地進行干擾的最強烈的方式,具有不可逆性,須通過制定土地利用總體規(guī)劃,合理利用土地和保護耕地,提高耕地質(zhì)量,保持城市化與糧食供應平衡。

曼谷早期城市發(fā)展缺乏合理規(guī)劃,產(chǎn)生了城市功能區(qū)劃混亂、交通堵塞、環(huán)境污染、人口過度膨脹等問題。面對這些問題,曼谷采取了一系列的措施,如在城市外圍修建住房、建立二級城鎮(zhèn)、在次發(fā)展地區(qū)發(fā)展工業(yè)、改善配套設施、興建交通道路等；此外還通過嚴格控制車輛數(shù)量和機動車燃油質(zhì)量降低汽車尾氣帶來的空氣污染；通過將城市中心固體垃圾運往老港填埋以減少城市環(huán)境污染等。在諸多政策的引導下,曼谷城市中心的周邊形成了多個城市次中心,如今的曼谷城市早已突破行政界線,發(fā)展成為以“曼谷中心城市、曼谷大都市、曼谷都市圈”為核心的特大城市群。該城市群具有強大的輻射能力,凸顯出巨大的經(jīng)濟創(chuàng)造力和潛力。隨著“21世紀海上絲綢之路”戰(zhàn)略實施,我國沿海城市和港口(如上海、廣州、福州、寧波、深圳、?？诘?正在做出積極響應,試圖從政策、交通、貿(mào)易、科技等方面進行改革[30- 31]。如廣東省政府頒布《廣東省參與建設“一帶一路”的實施方案》,欲將廣東打造成為戰(zhàn)略交通樞紐[30],寧波市發(fā)改委制定《寧波參與“一帶一路”建設行動綱要》,積極建設寧波港口經(jīng)濟圈[32]等。其中,加強港口及腹地群的建設成為主要措施。相比于內(nèi)陸城市,港口城市具有優(yōu)越的地理位置,但由于受海岸線的轄制,土地資源十分珍貴。在城市規(guī)劃中,國內(nèi)城市或可汲取曼谷城市發(fā)展的經(jīng)驗教訓,在以港口為依托,打造以港口、港區(qū)為核心的產(chǎn)業(yè)集聚區(qū)的基礎上,通過優(yōu)化配置城市土地利用結構,合理布置交通網(wǎng)絡體系,努力提高港口對經(jīng)濟腹地的輻射帶動作用和核心城市對區(qū)域小城市的帶動作用,打造共同繁榮的港口群和城市群。

總之,曼谷城市化不僅導致城市植被綠度降低、地表溫度升高等生態(tài)環(huán)境問題,同時過度城市化也帶來了住房擁擠、交通擁擠、基礎設施不足、失業(yè)率高等社會問題,這種“虛假城市化”的現(xiàn)象應引起所有發(fā)展中國家重視。基于曼谷虛假城市化和城市首位度畸高的教訓,我國港口城市應充分重視城市化和社會經(jīng)濟的同步發(fā)展,重視核心港口、城市與區(qū)域小港口、小城市的共同發(fā)展,在提高港口群、城市群的經(jīng)濟創(chuàng)造能力和輻射帶動能力時,也應避免出現(xiàn)類似于曼谷首位度畸高、區(qū)域發(fā)展不平衡等問題。

[1] Seto K C, Fragkias M, Güneralp B, Reilly M K. A Meta-Analysis of Global Urban Land Expansion. PLoS ONE, 2011, 6(8): e23777.

[2] Zhao S Q, Zhou D C, Zhu C, Sun Y, Wu W J, Liu S G. Spatial and temporal dimensions of urban expansion in China. Environmental Science & Technology, 2015, 49(16): 9600- 9609.

[3] Mertes C M, Schneider A, Sulla-Menashe D, Tatem A J, Tan B. Detecting change in urban areas at continental scales with MODIS data. Remote Sensing of Environment, 2015, 158: 331- 347.

[4] 匡文慧, 遲文峰, 史文嬌. 中國與美國大都市區(qū)城市內(nèi)部土地覆蓋結構時空差異. 地理學報, 2014, 69(7): 883- 895.

[5] Wu W J, Zhao S Q, Zhu C, Jiang J L. A comparative study of urban expansion in Beijing, Tianjin and Shijiazhuang over the past three decades. Landscape and Urban Planning, 2015, 134: 93- 106.

[6] 裴鳳松, 黎夏, 劉小平, 夏庚瑞. 城市擴張驅動下植被凈第一性生產(chǎn)力動態(tài)模擬研究——以廣東省為例. 地球信息科學學報, 2015, 17(4): 469- 477.

[7] Seto K C, Guneralp B, Hutyra L R. Global forecasts of urban expansion to 2030 and direct impacts on biodiversity and carbon pools. Proceedings of the National Academy of Sciences of the United States of America, 2012, 109(40): 16083- 16088.

[8] Bounoua L, Zhang P, Mostovoy G, Thome K, Masek J, Imhoff M, Shepherd M, Quattrochi D, Santanello J, Silva J, Wolfe R, Toure A M. Impact of urbanization on US surface climate. Environmental Research Letters, 2015, 10: 084010.

[9] Fan P L, Chen J Q, John R. Urbanization and environmental change during the economic transition on the Mongolian Plateau: Hohhot and Ulaanbaatar. Environmental Research, 2015, 144: 96- 112.

[10] Du N R, Ottens H, Sliuzas R. Spatial impact of urban expansion on surface water bodies—A case study of Wuhan, China. Landscape and Urban Planning, 2010, 94(3/4): 175- 185.

[11] Yan Y, Kuang W H, Zhang C, Chen C B. Impacts of impervious surface expansion on soil organic carbon—a spatially explicit study. Scientific Reports, 2015, 5: 17905.

[12] 羅上華, 毛齊正, 馬克明, 鄔建國. 城市土壤碳循環(huán)與碳固持研究綜述. 生態(tài)學報, 2012, 32(22): 7177- 7189.

[13] Estoque R C, Murayama Y. Intensity and spatial pattern of urban land changes in the megacities of Southeast Asia. Land Use Policy, 2015, 48: 213- 222.

[14] 宋偉軒, 朱喜鋼. 大湄公河次區(qū)域城市空間結構特征與成因. 經(jīng)濟地理, 2010, 30(1): 53- 58.

[15] United Nations. World urbanization prospects: the 2014 revision. New York: United Nations Publication, 2014.

[16] Arifwidodo S D, Tanaka T. The characteristics of urban heat island in Bangkok, Thailand. Procedia-Social and Behavioral Sciences, 2015, 195: 423- 428.

[17] 王翠平, 王豪偉, 李春明, 董仁才. 基于DMSP/OLS影像的我國主要城市群空間擴張?zhí)卣鞣治? 生態(tài)學報, 2012, 32(3): 942- 954.

[18] Liu X P, Li X, Chen Y M, Tan Z Z, Li S Y, Ai B. A new landscape index for quantifying urban expansion using multi-temporal remotely sensed data. Landscape Ecology, 2010, 25(5): 671- 682.

[19] Xu C, Liu M S, Zhang C, An S Q, Yu W, Chen J M. The spatiotemporal dynamics of rapid urban growth in the Nanjing metropolitan region of China. Landscape Ecology, 2007, 22(6): 925- 937.

[20] 覃志豪, Zhang M H, Arnon K, Pedro B. 用陸地衛(wèi)星TM6數(shù)據(jù)演算地表溫度的單窗算法. 地理學報, 2001, 56(4): 456- 466.

[21] 覃志豪, 李文娟, 徐斌, 陳仲新, 劉佳. 陸地衛(wèi)星TM6波段范圍內(nèi)地表比輻射率的估計. 國土資源遙感, 2004, (3): 28- 32, 36- 36, 41- 41.

[22] Rozenstein O, Qin Z H, Derimian Y, Karnieli A. Derivation of land surface temperature for Landsat- 8 TIRS using a split window algorithm. Sensors (Basel), 2014, 14(4): 5768- 5780.

[23] 張勇, 余濤, 顧行發(fā), 張玉香, 陳良富, 余姍姍, 張文君, 李小文. CBERS-02 IRMSS熱紅外數(shù)據(jù)地表溫度反演及其在城市熱島效應定量化分析中的應用. 遙感學報, 2006, 10(5): 789- 797.

[24] 繆麗娟, 崔雪鋒, 欒一博, 何斌. 北京上海近20a城市化過程中土地利用變化異同點探析. 氣象科學, 2011, 31(4): 398- 404.

[25] 中華人民共和國商務部網(wǎng)站. [2016-07- 16]. http://th.mofcom.gov.cn/.

[26] 紀雁, 沙永杰. 曼谷: 邁向可持續(xù)和包容性發(fā)展面臨的挑戰(zhàn). 上海城市規(guī)劃, 2015, (3): 61- 67, 93- 93.

[27] 吝濤, 李新虎, 張國欽, 趙千鈞, 崔勝輝. 廈門島城市空間擴張?zhí)卣骷捌溆绊懸蛩胤治? 地理學報, 2010, 65(6): 715- 726.

[28] 高義, 蘇奮振, 孫曉宇, 楊曉梅, 薛振山, 張丹丹. 近20a廣東省海島海岸帶土地利用變化及驅動力分析. 海洋學報, 2011, 33(4): 95- 103.

[29] 孟林. 城市化背景下土地利用變化及其生態(tài)環(huán)境質(zhì)量評價——以環(huán)渤海沿海城市為例[D]. 大連: 遼寧師范大學, 2012.

[30] 鄒偉勇, 金祎, 熊曉冬. 廣東建設“21世紀海上絲綢之路”戰(zhàn)略樞紐的交通協(xié)同策略. 規(guī)劃師, 2016, 32(2): 38- 45.

[31] 劉宗義. 21世紀海上絲綢之路建設與我國沿海城市和港口的發(fā)展. 城市觀察, 2014, 34(6): 5- 12.

[32] 徐瑩, 董文娟, 張雪梅. “一帶一路”戰(zhàn)略下寧波港拓展中西部腹地策略研究. 華東交通大學學報, 2015, 32(6): 143- 152.

UrbanexpansionanditsecologicalenvironmentaleffectsinBangkok,Thailand

HAN Ruidan1, 2, ZHANG Li2,3,*, ZHENG Yi3, WANG Heng3, ZHANG Jing3

1CollegeofInformationScienceandEngineeringofShandongAgriculturalUniversity,Tai′an271018,China2HainanKeyLaboratoryofEarthObservation,Sanya572029,China3KeyLaboratoryofDigitalEarthScience,InstituteofRemoteSensingandDigitalEarth,ChineseAcademyofSciences,Beijing100094,China

On the basis of Landsat images acquired in 1990, 2000, 2010, and 2015 and MODIS NDVI datasets acquired from 2010 to 2015, we analyzed urban expansion in Bangkok, including the spatio-temporal characteristics of built-up land, ecological-environmental effects (including vegetation greenness and land surface temperature), and the main driving factors of urban expansion. The results showed that (1) the built-up land areas in Bangkok has expanded continuously since 1990, showing a rapid-slow-rapid change trajectory. Urban expansion occurred mainly in the northeast direction and within the range of 14—20 km from the city center. The expansion patterns followed infilling pattern and edge-expansion pattern. (2) Urbanization affected vegetation greenness and heat island. Vegetation greenness declined overall, and a significantly decreased area was in the range of 10—20 km from the city center, which was in accordance with the main urban expansion area. Similarly, the spatial distribution of the heat island was coincident with urban expansion. The heat island increased overall, although local variations were less obvious. (3) Geographical conditions, social economy, and city layout were the main driving factors in urban expansion.

urban expansion; vegetation greenness; heat island; driving factors

海南省科技合作專項資金項目(KJHZ2015- 14);海南省重大科技計劃項目(ZDKJ2016021)

2016- 07- 19; < class="emphasis_bold">網(wǎng)絡出版日期

日期:2017- 05- 27

*通訊作者Corresponding author.E-mail: zhangli@radi.ac.cn

10.5846/stxb201607191470

韓瑞丹,張麗,鄭藝,王恒,張靜.曼谷城市擴張生態(tài)環(huán)境效應.生態(tài)學報,2017,37(19):6322- 6334.

Han R D, Zhang L, Zheng Y, Wang H, Zhang J.Urban expansion and its ecological environmental effects in Bangkok, Thailand.Acta Ecologica Sinica,2017,37(19):6322- 6334.