珠江口灣區(qū)沿岸土地利用時空格局變化分析

高 群，吳志峰*，李少英，曹 崢，王 琎

（1.廣州大學(xué)　地理科學(xué)學(xué)院，廣東　廣州 510006；2.廣東省生態(tài)環(huán)境與土壤研究所，廣東省農(nóng)業(yè)環(huán)境綜合治理重點實驗室，廣東　廣州 510650；3.中國科學(xué)院　廣州地球化學(xué)研究所，廣東　廣州 510640）

珠江口灣區(qū)沿岸土地利用時空格局變化分析

高 群1，2，吳志峰1，2*，李少英1，曹 崢2，3，王 琎2，3

（1.廣州大學(xué)地理科學(xué)學(xué)院，廣東廣州 510006；2.廣東省生態(tài)環(huán)境與土壤研究所，廣東省農(nóng)業(yè)環(huán)境綜合治理重點實驗室，廣東廣州 510650；3.中國科學(xué)院廣州地球化學(xué)研究所，廣東廣州 510640）

文章以多源遙感影像為基礎(chǔ)，獲取珠江口灣區(qū)1990年、1995年、2004年和2012年土地利用分類數(shù)據(jù)，并利用4個指數(shù)模型，對珠江口灣區(qū)3個時期土地利用結(jié)構(gòu)時空格局變化差異進行對比研究.結(jié)果表明：珠江口灣區(qū)建設(shè)用地主要通過占用農(nóng)用地和林地不斷擴張，其主要驅(qū)動力為人類活動；3個時期變化較活躍的類型分別為圍墾灘涂、建設(shè)用地和養(yǎng)殖用地，且初期綜合土地利用動態(tài)度值高于中期和后期；珠江口灣區(qū)西岸土地利用類型重心呈現(xiàn)向珠?？繑n趨勢；建設(shè)用地整體重心受地理位置、發(fā)展模式、填海造地和政府決策影響，逐年向西岸遷移；圍墾灘涂整體重心向東岸移動，珠江口灣區(qū)5市重心移動具有明顯的向海和向水域特征；灣區(qū)填海擴張方式主要表現(xiàn)為港口、道路、自貿(mào)區(qū)、公園以及機場建設(shè)，其中深圳市發(fā)展模式更符合可持續(xù)發(fā)展需要.

土地利用；動態(tài)度；轉(zhuǎn)移矩陣；空間格局重心；時空格局變化

珠江口灣區(qū)是大珠三角城鎮(zhèn)群的核心空間，改革開放以來，經(jīng)濟發(fā)展迅猛，城市快速擴張.灣區(qū)在自然環(huán)境變化和人類活動共同影響下，土地利用發(fā)生巨大變化.開展灣區(qū)土地利用變化規(guī)律及其驅(qū)動力的研究，不僅可以為灣區(qū)規(guī)劃和生態(tài)環(huán)境提出科學(xué)建議，也可為全球變化的區(qū)域模式研究提供典型的范例.

國內(nèi)學(xué)者已從土地利用變化遙感動態(tài)監(jiān)測、模型動態(tài)模擬以及驅(qū)動力分析等方面對海岸帶進行了大量研究［1-5］，研究范圍主要集中在長江、黃河、環(huán)渤海帶和珠江三角洲等區(qū)域，側(cè)重應(yīng)用土地利用時空變化分析模型探究海岸帶土地利用變化規(guī)律.現(xiàn)有對海岸帶的研究中，大多數(shù)將海岸帶進行梯度劃分，分析各條帶間土地利用變化差別，但忽略了海岸帶包含的城市單位，未考慮海岸帶內(nèi)城市之間的土地利用變化的差異.本文利用空間重心轉(zhuǎn)移指數(shù)對海岸帶內(nèi)城市之間的土地利用變化差異進行對比分析，研究其對灣區(qū)整體變化的影響.并且在數(shù)據(jù)選擇上，采用SPOT和Worldview影像，大大提高了本文分析的空間精度.對灣區(qū)土地利用類型變化過程、特征以及演變規(guī)律的研究，能夠為珠江三角洲城市群生態(tài)保護提供準確、全面的基礎(chǔ)數(shù)據(jù)以及決策支持，具有重要的現(xiàn)實意義.

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

1.1研究區(qū)概況

為配合《珠江三角洲地區(qū)改革發(fā)展規(guī)劃綱要（2008～2020年）》，粵港澳三地政府聯(lián)合提出了“珠江口灣區(qū)”的概念（以下簡稱灣區(qū)）.灣區(qū)范圍介于21°48′～23°09′N、113°09′～114°37′E之間，共包括廣州、深圳、珠海、東莞、中山等5市所轄的17個區(qū)，毗鄰香港和澳門特別行政區(qū).其擁有完善的基礎(chǔ)設(shè)施、全覆蓋的道路網(wǎng)絡(luò)和完備的森林、海洋、濕地生態(tài)系統(tǒng)，已經(jīng)成為大珠三角地區(qū)的地理核心和生態(tài)核心.

本研究在珠江口灣區(qū)的基礎(chǔ)上，以1979年海岸線為基準，將沿海岸線向陸一側(cè)10 km與沿海市縣行政界線或沿海城市行政邊界的最大值連線所圍成的閉合區(qū)域作為研究區(qū)，見圖1.

圖1　珠江口灣區(qū)研究區(qū)Fig.1 Study area of Pearl River estuary

1.2數(shù)據(jù)來源與預(yù)處理

本文采用的遙感數(shù)據(jù)為Landsat TM/ETM+影像、SPOT影像及Worldview影像.時相分別為1990年、1995年、2004年、2012年，并對 LANDSAT 7 ETM+數(shù)據(jù)進行了條帶修復(fù).

對各時相影像進行預(yù)處理.首先，基于2003 年1∶10 000地形圖（DLG）數(shù)據(jù)，對數(shù)據(jù)進行幾何糾正.采用FLAASH模型對影像進行大氣校正.然后，在遙感模型與人工目視解譯的基礎(chǔ)上，利用面向?qū)ο蠓椒?，?期遙感影像分為建設(shè)用地、林地、農(nóng)用地、水體、養(yǎng)殖用地和圍墾灘涂等6種類型，并針對“異物譜”問題導(dǎo)致地類間的誤分進行人工修改.最后，在研究區(qū) 2012年土地利用分類數(shù)據(jù)中隨機選取了400個驗證點，通過土地利用提取結(jié)果與相片、SPOT和Worldview影像比對匹配，得到土地利用分類精度為85.5%.

1.3研究方法

1.3.1單一土地利用類型動態(tài)度

單一土地利用類型動態(tài)度［6］是反映研究時期內(nèi)某種土地利用類型變化情況的指標，其表達式為

式（1）中，K為單一土地利用類型動態(tài)度；U1、U2分別為研究初期及研究末期某一種土地利用類型的面積［6］.

1.3.2綜合土地利用類型動態(tài)度

綜合土地利用動態(tài)度［6］是反映研究時期內(nèi)全部土地利用類型綜合變化情況的指標，其表達式為

式（2）中，LC為綜合土地利用類型動態(tài)度；LUa為研究初期第a類土地利用類型面積；ΔLUa-b為研究范圍內(nèi)a類土地利用類型轉(zhuǎn)為非a類土地利用類型面積的絕對值［6］.

1.3.3土地利用類型空間格局重心

土地利用類型空間格局重心模型［7-9］可從空間上反映土地利用類型的時空演變過程.

重心坐標計算公式：

式（3）、式（4）中，X、Y分別為某一個區(qū)域某種類型重心的經(jīng)度值和緯度值；Xi、Yi分別為第i個次一級區(qū)域中心的經(jīng)度值和緯度值；Mi此處代表第i個次一級區(qū)域的面積值［7-9］.

重心遷移距離公式：

式（5）中，D為研究區(qū)重心前移距離X（t1）；Y（t1）、X（t2）、Y（t2）分別為t1和t2時刻研究區(qū)重心的橫坐標與縱坐標；t1和t2分別為監(jiān)測首期與末 期［7-9］.

2　結(jié)果與分析

2.1珠江口灣區(qū)土地利用數(shù)量變化及變化過程

珠江口灣區(qū)主要用地類型為水體、建設(shè)用地、林地和農(nóng)用地，其次為養(yǎng)殖用地和圍墾灘涂，見表1.由表1可知，除水體外，建設(shè)用地、林地和農(nóng)用地面積分別保持在5.89%、10.39%和13.68%以上.監(jiān)測時段內(nèi)，①建設(shè)用地保持持續(xù)增長，1990 ～1995年增幅最大，為178.88%，并在2012年達到最大值，占總面積的33.40%；②林地呈波動下降趨勢，1990～1995年降幅最大，為50.29%；③農(nóng)用地面積初期一直波動但面積變化不大，2004～2012年降幅最為明顯，達到42.24%；④水體面積持續(xù)下降，到2012年達到最低值，占總面積的31.96%；⑤養(yǎng)殖用地面積不斷波動，1990年占總面積6.28%，到2012年增長為9.87%.

表1　1990～2012年研究區(qū)土地利用面積統(tǒng)計Table 1 Land use area of the study area from 1990 to 2012　km2

轉(zhuǎn)移矩陣結(jié)果（表2～4）表明：①建設(shè)用地的轉(zhuǎn)入量遠遠大于轉(zhuǎn)出量，后期達到總面積的48.28%.直到2004～2012年，建設(shè)用地面積才有了明顯的轉(zhuǎn)出，轉(zhuǎn)化為農(nóng)用地和養(yǎng)殖用地比較明顯，面積分別達到了11.33%和3.81%.建設(shè)用地面積的持續(xù)增長表明建設(shè)用地擴張是珠江口灣區(qū)空間擴張的主要表現(xiàn)形式之一；②林地的轉(zhuǎn)出量遠大于轉(zhuǎn)入量，轉(zhuǎn)入的土地利用類型基本為農(nóng)用地和養(yǎng)殖用地.大部分林地轉(zhuǎn)出為建設(shè)用地和農(nóng)用地，小部分轉(zhuǎn)化為養(yǎng)殖用地，前期的轉(zhuǎn)出量最為明顯；③農(nóng)用地在第三個時期的轉(zhuǎn)出量較大，轉(zhuǎn)入類型和轉(zhuǎn)出類型基本相同，為建設(shè)用地、養(yǎng)殖用地和林地，其他則變化較小；④水體的轉(zhuǎn)入量是極少的，大部分水體轉(zhuǎn)變?yōu)轲B(yǎng)殖用地和建設(shè)用地；⑤養(yǎng)殖用地主要轉(zhuǎn)化為建設(shè)用地、農(nóng)用地和水體，轉(zhuǎn)入類型為水體、農(nóng)用地和林地.

表2　1990～1995年土地利用轉(zhuǎn)移矩陣Table 2 Conversion matrix of land use between 1990 and 1995

表3　1995～2004年土地利用轉(zhuǎn)移矩陣Table 3 Conversion matrix of land use between 1995 and 2004

表4　2004～2012年土地利用轉(zhuǎn)移矩陣Table 4 Conversion matrix of land use between 2004 and 2012

2.2珠江口灣區(qū)土地利用動態(tài)度變化分析

單一土地利用動態(tài)度K和綜合土地利用動態(tài)度LC（圖2）表明，1990～1995年，6種土地利用類型動態(tài)度差異較大，正向以圍墾灘涂的變化率為最高，達到66.98%，由高到低排序為圍墾灘涂＞建設(shè)用地＞養(yǎng)殖用地，負向動態(tài)度由高到低為林地＞水體＞農(nóng)用地.同時期綜合土地利用類型動態(tài)度最高，達到2.73%.以上表明該時期土地利用類型變化活躍，自19世紀80年代以來，在珠江三角洲一帶設(shè)立改革開放試點、改革農(nóng)村經(jīng)濟體制等措施，促使珠江三角洲經(jīng)濟飛速發(fā)展，建設(shè)用地迅速擴張［10］.為了緩解緊張的土地供求關(guān)系，開始出現(xiàn)大量的填海造陸現(xiàn)象.這一時期，林地、農(nóng)業(yè)用地逐漸減少，養(yǎng)殖用地迅速擴張，圍墾面積逐漸擴大.

圖2　單一土地利用類型動態(tài)度K和綜合土地利用類型動態(tài)度LCFig.2 Dynamic index of each land use type（K）and dynamic index of land use change（LC）

1995～2004年，6種土地類型間的轉(zhuǎn)化并不活躍，綜合土地利用類型動態(tài)度降為1.11%，且各土地利用類型間的動態(tài)度差異較小，正向動態(tài)度由高到低排序為建設(shè)用地＞林地＞農(nóng)用地，負向動態(tài)度由高到低排序為養(yǎng)殖用地＞圍墾灘涂＞水體.在經(jīng)歷了經(jīng)濟發(fā)展的“大躍進”之后，經(jīng)濟“泡沫”逐漸顯露，政府出臺了相應(yīng)的價格調(diào)控和土地開發(fā)等政策，使珠江三角洲經(jīng)濟發(fā)展逐步降溫，邁入調(diào)整階段［11］.

2004～2012年，除養(yǎng)殖用地略有反彈之外，其他用地類型與上一個時間段差別不大，正向動態(tài)度由高到低排序為養(yǎng)殖用地＞建設(shè)用地，負向動態(tài)度由高到低排序為農(nóng)用地＞林地＞水體＞圍墾灘涂.綜合土地利用類型動態(tài)度略有上升為1.59%，但低于前期水平.這一時期珠江口灣區(qū)的土地利用變化與政策背景、政策支持密切相關(guān).20多年間珠江口灣區(qū)土地利用單一化程度和集約利用水平呈現(xiàn)上升的趨勢，先后經(jīng)歷了高速發(fā)展期、平穩(wěn)發(fā)展期，并正在積極打造特色發(fā)展模式.

2.3珠江口灣區(qū)土地利用時空格局重心變化

以上指數(shù)分析表明，珠江口灣區(qū)的擴張主要表現(xiàn)為建設(shè)用地的擴張，根據(jù)式（3）和（4），計算研究區(qū)整體和5市建設(shè)用地重心，得到圖3和圖4.在繪制圖4時，由于年際間變化不大，故只標示首期和末期重心，進行重心遷移趨勢的研究.

圖3表明，1990～2012年，研究區(qū)整體重心一直向灣區(qū)西岸遷移，說明西岸建設(shè)用地增加速度高于東岸.根據(jù)式（5），計算得到前期和后期遷移距離分別是中期的2.99倍和2.11倍，說明中期的東西兩岸建設(shè)用地增長較其他2個時期更加均衡.

圖4表明，灣區(qū)東岸東莞市和深圳市重心向西北方向移動，遷移距離分別為17.85 km和6.36 km.灣區(qū)西岸除廣州市向東南方向移動，遷移距離為10.64 km，中山市和珠海市均向西南方向移動，遷移距離分別為2.04 km和10.95 km.從5市的遷移趨勢上看，除廣州市外，其他市均有向西遷移趨勢，雖然廣州市重心向東南方向移動，但其向東的遷移角度很小，故這是灣區(qū)整體重心逐年向西遷移的重要原因之一.

圖3　研究區(qū)建設(shè)用地空間格局重心Fig.3 Spatial center of construction land in study area

圖4　研究區(qū)5市建設(shè)用地空間格局重心Fig.4 Spatial center of construction land in five cities of study area

廣州市、珠海市、東莞市和深圳市建設(shè)用地重心均有向海和向水域遷移趨勢，這是由于海岸帶城市獨特的地理位置優(yōu)勢，促使外向型經(jīng)濟不斷發(fā)展，沿海岸線不斷擴張.由于城市規(guī)劃和布局不合理，導(dǎo)致港口倉儲能力不足，灣區(qū)內(nèi)不斷填海造地，這也是重心遷移的重要原因之一.雖然灣區(qū)大部分城市都受到了上述原因的影響，但中山市重心遷移方向卻剛好相反.中山市地處珠江出?？冢?jīng)濟發(fā)展主要以第二產(chǎn)業(yè)為主，以上說明沿海城市建設(shè)用地重心遷移不僅受到地理位置和填海擴張的影響，更與城市自身發(fā)展模式和政府宏觀決策緊密相關(guān).

整個灣區(qū)內(nèi)，除建設(shè)用地外其他土地利用類型空間重心變化總結(jié)如下：農(nóng)用地重心變化不大，西岸始終變化幅度較??；林地重心由東岸向西岸移動，說明東岸林地減少幅度大于西岸；灣區(qū)水體向西岸遷移距離較大，說明東岸水體減少程度大于西岸；養(yǎng)殖用地在西岸分布較多；整體上，珠江口西岸土地利用類型分布較均勻，呈現(xiàn)向珠海靠攏的趨勢.

圍墾灘涂作為海岸帶的特色用地類型，對其深入研究有特殊的意義.珠江口灣區(qū)的圍墾活動主要分布在廣州市、珠海市和深圳市.灣區(qū)圍墾灘涂的整體重心向東岸移動，5市的重心轉(zhuǎn)移具有明顯的向海和向水域特征.本文以廣州市、珠海市和深圳市為例，通過對其圍墾灘涂區(qū)域結(jié)構(gòu)變化的監(jiān)測，了解海岸帶填海擴張的過程及其影響因素.

廣州市圍墾活動主要集中在萬頃沙和龍穴島一帶.萬頃沙位于廣州市的最南端，水陸交通便利，是廣州的重要港口之一.監(jiān)測時段內(nèi)，其圍墾面積不斷增加并向東南方向延伸，萬頃沙主要發(fā)展旅游觀光農(nóng)業(yè)、造船、物流等產(chǎn)業(yè)，并形成了獨特的沙田水鄉(xiāng)文化；龍穴島面積擴張的方式與萬頃沙基本相同，但其經(jīng)濟側(cè)重點略有不同.龍穴島發(fā)展至今，圍墾面積增長了幾倍，其經(jīng)濟發(fā)展主要靠港口物流產(chǎn)業(yè)帶動.

深圳市以深圳灣和前海區(qū)域為例進行研究，國務(wù)院批復(fù)建設(shè)深圳前海自貿(mào)區(qū)，前海圍墾區(qū)域大部分轉(zhuǎn)變?yōu)榻ㄔO(shè)用地，港口建設(shè)、道路建設(shè)、電廠建設(shè)、油港建設(shè)、發(fā)展物流產(chǎn)業(yè)等，使其形成了較完備的沿海經(jīng)濟體系；深圳灣一帶圍墾活動更注重對沿海濕地的保護，公園建設(shè)是深圳灣發(fā)展特色.深圳的創(chuàng)新經(jīng)濟走在全國的最前端，經(jīng)濟發(fā)展的活力和創(chuàng)造力都使得其在邁向產(chǎn)業(yè)升級的中國經(jīng)濟地圖中，有無可取代的地位.

珠海市圍墾活動以高欄島、三灶島為例進行研究.高欄島位于崖門口東側(cè)，雞啼門出口處，地理優(yōu)勢造就了高欄島獨特的地理風(fēng)光和人文環(huán)境，吸引了國內(nèi)外的許多游客，其填海造地活動主要為島內(nèi)與大陸間的道路建設(shè)，到2012年高欄島已經(jīng)和大陸相連，沒有明顯界限；三灶島填?；顒又饕捎跈C場建設(shè)的需要，增加面積主要為建設(shè)用地.

2.4珠江口灣區(qū)土地利用景觀格局變化分析

計算各年份珠江口灣區(qū)主要用地類型的景觀格局指數(shù)，斑塊密度指數(shù)（PD）、斑塊平均面積指數(shù)（AREA-MN）、面積加權(quán)的平均斑塊分形指數(shù)（FRAC-MN），得到表5.分析表明：①建設(shè)用地的斑塊密度指數(shù)呈現(xiàn)減少趨勢，研究區(qū)域建設(shè)用地景觀破碎度減小，即在人為因素的干預(yù)下建設(shè)用地呈現(xiàn)聚集態(tài)勢，零散的建設(shè)用地斑塊逐年減少［12］；②農(nóng)用地與林地的斑塊密度指數(shù)則波動上升并且農(nóng)用地的波動幅度大于林地，即由于人為因素的影響，大量連片出現(xiàn)的農(nóng)用地與林地被其它用地類型分割成零散的農(nóng)用地與林地；③斑塊平均面積指數(shù)變化趨勢同樣表明建設(shè)用地破碎化程度降低，農(nóng)用地、林地的破碎化程度增加；④面積加權(quán)的平均斑塊分形指數(shù)表示景觀格局的不規(guī)則程度和破碎程度，建設(shè)用地、農(nóng)用地以及林地的分形指數(shù)均小于1.30，從整體上看，從大到小依次排序為建設(shè)用地＞農(nóng)用地＞林地，由于建設(shè)用地、農(nóng)用地以及林地形狀接近正方形或圓形，更易受人類活動影響.

通過上述分析不難看出，建設(shè)用地、農(nóng)用地及林地均受人類活動影響較為明顯，3者受人類活動影響由強至若依次為建設(shè)用地、農(nóng)用地、林地.因此，珠江口土地利用變化的主要驅(qū)動力為人類活動.

表5　主要土地利用類型景觀格局指數(shù)Table 5 The Landscape metrics of main land use

3　結(jié) 論

本文以多源遙感影像為基礎(chǔ)，采用遙感與GIS分析方法，對珠江口灣區(qū)的土地利用時空格局變化進行了23 a的監(jiān)測，根據(jù)4個指數(shù)模型，對灣區(qū)在3個時期土地利用數(shù)量變化、結(jié)構(gòu)變化、動態(tài)度變化和區(qū)域重心變化差異及主要驅(qū)動因素進行研究，為珠江三角洲海岸帶研究提供較好的補充和完善，得到如下結(jié)論：

（1）土地利用結(jié)構(gòu)變化表明，1990～2012年間，土地利用類型變化主要表現(xiàn)為農(nóng)用地和林地向建設(shè)用地的轉(zhuǎn)化，建設(shè)用地大幅增加，而農(nóng)用地和林地迅速減少，其變化主要驅(qū)動力為人類活動.

（2）單一土地利用動態(tài)度表明，3個時期土地利用類型變化較大，分別為圍墾灘涂、建設(shè)用地和養(yǎng)殖用地.綜合土地利用動態(tài)度表明，監(jiān)測初期動態(tài)度值遠高于中期和后期.監(jiān)測時段內(nèi)珠江口灣區(qū)土地利用單一化程度和集約利用水平呈上升趨勢.

（3）土地利用空間重心移動表明：珠江口灣區(qū)西岸土地利用類型重心呈向珠?？繑n趨勢；建設(shè)用地重心受地理位置、發(fā)展模式、填海造地和政府宏觀決策影響，逐年向西岸轉(zhuǎn)移；農(nóng)用地和林地呈內(nèi)陸分布，且東岸林地消失較快；水體主要集中在地勢較低的區(qū)域，且西岸水體減少較快；圍墾灘涂整體重心向東岸移動，5市重心轉(zhuǎn)移具有明顯的向海和向水域特征；灣區(qū)填海擴張方式主要表現(xiàn)為港口、道路、自貿(mào)區(qū)、公園以及機場建設(shè)，并受經(jīng)濟利益和政府決策共同影響，其中深圳市的發(fā)展模式更符合可持續(xù)發(fā)展需要.

［1］ VERBURG P H，EICKHOUT B，VAN MEIJL H.A multi-scale，multi-model approach for analyzing the future dynamics of European land use［J］.Ann Reg Sci，2008，42（1）：57-77.

［2］ 王玉，賈曉波，張文廣，等.江蘇海岸帶土地利用變化及驅(qū)動力分析［J］.長江流域資源與環(huán)境，2010，19（S1）：7-12. WANG Y，JIA X B，ZHANG W G，et al.Land use change and analysis on the driving forces in the coastal zone of Jiangsu ［J］.Resour Environ Yangtze Basin，2010，19（S1）：7-12.

［3］ 高義，蘇奮振，孫曉宇，等.近20a廣東省海島海岸帶土地利用變化及驅(qū)動力分析［J］.海洋學(xué)報，2011，33（4）：95-103. GAO Y，SU F Z，SUN X Y，et al.A study on driving forces of land use change of Guangdong Province coastal zone and islands in recent 20 a［J］.Acta Oceanol Sin，2011，33（4）：95-103.

［4］ 侯西勇，徐新良.21世紀初中國海岸帶土地利用空間格局特征［J］.地理研究，2011，33（8）：1370-1379. HOU X Y，XU X L.Spatial patterns of land use in coastal zones of China in the early 21st century［J］.Geogr Res，2011，33 （8）：1370-1379.

［5］ 韓磊，侯西勇，朱明明，等.20世紀后半葉美國海岸帶區(qū)域土地利用變化時空特征分析［J］.世界地理研究，2010，19 （2）：42-52. HAN L，HOU X Y，ZHU M M，et al.Study on the temporal-spatial characters of land-use change in the coastal zone of America in the latter half of the 20th century［J］.World Reg Stud，2010，19（2）：42-52.

［6］ 王德智，邱彭華，方源敏，等.?？谑泻０稁恋乩脮r空格局變化分析［J］.地球信息科學(xué)學(xué)報，2014，16（6）：933-940. WANG D Z，QIU P H，F(xiàn)ANG Y M，et al.Research on spatial-temporal changes of land use in Haikou city’s coastal zone ［J］.J Geo-inform Sci，2014，16（6）：933-940.

［7］ 匡文慧，張樹文，張養(yǎng)貞，等.1990年以來長春市土地利用空間擴張機理分析［J］.地理學(xué)報，2005，60（5）：841-850. KUANG W H，ZHANG S W，ZHANG Y Z，et al.Analysis of urban land utilization spatial expansion mechanism in Changchun city since 1990［J］.Acta Geogr Sin，2005，60（5）：841-850.

［8］ 劉保曉，黃耀歡，付晶瑩，等.天津港區(qū)土地利用時空格局變化與驅(qū)動力分析［J］.地球信息科學(xué)學(xué)報，2012，14（2）：270-278. LIU B X，HUANG Y H，F(xiàn)U J Y，et al.Analysis on spatio-temporal change and driving forces of land use in Tianjin harbor ［J］.J Geo-inform Sci，2012，14（2）：270-278.

［9］ 秦文翠，羅維，劉運明.天津濱海新區(qū)海岸帶土地利用時空格局變化［J］.西南師范大學(xué)學(xué)報：自然科學(xué)版，2015，40 （5）：135-141. QIN W C，LUO W，LIU Y M.On spatio-temporal changes of land use in a coastal zone，China：A case study in Tianjin Binhai new area［J］.J SW China Norm Univ：Nat Sci Edi，2015，40（5）：135-141.

［10］周銳波.珠江三角洲經(jīng)濟發(fā)展模式評析［J］.中國發(fā)展，2004，4（3）：15-20. ZHOU R B.Pearl River Delta economic development：An review of its process and mode［J］.Chin Develop，2004，4（3）：15-20.

［11］葉玉瑤，張虹鷗，許學(xué)強，等.珠江三角洲建設(shè)用地擴展與經(jīng)濟增長模式的關(guān)系［J］.地理研究，2011，30（12）：2259-2271. YE Y Y，ZHANG H O，XU X Q，et al.Study on the relationship between construction land expansion and economic growth pattern in the Pearl River Delta［J］.Geogr Res，2011，30（12）：2259-2271.

［12］孫曉宇，蘇奮振，呂婷婷，等.珠江口西岸土地利用景觀格局變化及驅(qū)動力分析［J］.地球信息科學(xué)學(xué)報，2009，11（4）：436-441. SUN X Y，SU F Z，LV T T，et al.Analysis on landscape spatial pattern changes of land use and its driving force in coastal zone of western Pearl River estuary［J］.J Geo-inform Sci，2009，11（4）：436-441.

【責(zé)任編輯：孫向榮】

Research on spatial-temporal changes of land use in Pearl River estuary

GAO Qun1，2，WU Zhi-feng1，2*，Li Shao-ying1，CAO Zheng2，3，WANG Jin2，3
（1.School of Geographical Sciences，Guangzhou University，Guangzhou 510006，China；2.Guangdong Key Laboratory of Agricultural Environment Pollution Integrated Control，Guangdong Institute of Eco-environmental and Soil Sciences，Guangzhou 510650，China；3.Guangzhou Institute of Geochemistry，Chinese Academy of Sciences，Guangzhou 510640，China）

The data of land use in the bay area of Pearl River estuary are based on the multi-source remote sense in three periods from 1990 to 2012.During these three time segments，temporal and spatial variations are studied by GIS software to calculate the dynamic index of each land use type.The results show that the construction land of the bay area increase rapidly by occupying agriculture land and wood land，its main driving force being human activities；during the time span，tidal flat，construction land and farm land are the highest dynamic types，and the first period is the most active period.In the west bank，the center of all land use types tend to transfer to the Pearl River.The reasons behind construction land center transition are the location，developing mode，land reclamation and government policy.The overall center of tidal flat moves to the east coast，while the city centers move significantly to the sea.The reclamation in the bay area are characterized by transformation to ports，roads，free trade areas，parks，airports，while the Shenzhen model is more in line with the needs of sustainable development.

land use；dynamic index；conversion matrix；spatial pattern center；spatial-temporal change

P 748；TP 79

A

1671-4229（2016）03-0073-08

2016-03-02；

2016-04-19

“十二五”國家科技支撐計劃資助項目“村鎮(zhèn)建設(shè)用地再開發(fā)關(guān)鍵技術(shù)研究與示范”（2013BAJ13B01）

高 群（1991-），女，碩士研究生.E-mail：gzdxgq@163.com

.E-mail：gzuwzf@163.com