東莞市人工景觀結(jié)構(gòu)與連接度的協(xié)整性及因果關(guān)系

馮志新,陳穎彪,千慶蘭,王帥帥

1 廣東省環(huán)境科學研究院， 廣州 510031 2 廣州大學地理科學學院，廣州 510006

東莞市人工景觀結(jié)構(gòu)與連接度的協(xié)整性及因果關(guān)系

馮志新1,2,陳穎彪2,*,千慶蘭2,王帥帥2

1 廣東省環(huán)境科學研究院， 廣州 510031 2 廣州大學地理科學學院，廣州 510006

城市人工景觀特征逐步向城市外圍自然及半自然景觀滲透。近年來對快速城市化地區(qū)的景觀格局特征定量化研究很多，但側(cè)重于景觀格局特征的簡單分析和描述，忽略了城市空間擴展過程中，城市人工景觀結(jié)構(gòu)與景觀連接度動態(tài)變化過程中之間的內(nèi)部聯(lián)系。以東莞市1988—2009年城市人工景觀為研究對象，通過計算并分析景觀結(jié)構(gòu)與景觀連接度指數(shù)，利用計量經(jīng)濟學中的協(xié)整概念與格蘭杰因果檢驗方法并詳細分析了這兩類指數(shù)之間的具體關(guān)系。研究結(jié)果表明：(1)所選取的景觀結(jié)構(gòu)指數(shù)與景觀連接度指數(shù)能較好地描述東莞市城市人工景觀時空格局變化過程；(2)人工景觀斑塊類型面積(CA)與可能連通性指數(shù)(PC)存在協(xié)整關(guān)系。景觀聚集度指數(shù)(AI)、蔓延度指數(shù)(CONTAG)與景觀連通性指數(shù)(IIC與PC)也存在協(xié)整關(guān)系，即這些變量之間確實存在長期均衡關(guān)系；(3)格蘭杰因果關(guān)系檢驗中，在不同滯后期的情況下，大部分的結(jié)果表明景觀連接度指數(shù)(IIC與PC)是景觀結(jié)構(gòu)指數(shù)(CA、AI、CONTAG)的格蘭杰成因，說明景觀連接度指數(shù)信息有助于預測景觀結(jié)構(gòu)指數(shù)的變化。

城市擴張；景觀結(jié)構(gòu)；景觀連接度；東莞市

城市化復雜地表過程模擬表達是目前重點研究的一個科學問題，從景觀結(jié)構(gòu)的空間構(gòu)成看，城市景觀包括基底、斑塊、廊道三大要素。城市是典型的人工景觀，建筑物群體和硬化地面構(gòu)成了景觀的主體，街區(qū)和街道是城市景觀的基底[1]。近年來運用RS與GIS技術(shù)手段，通過選取不同尺度的景觀格局指數(shù)，采用景觀格局與梯度的分析方法，對快速城市化地區(qū)的景觀格局特征定量化研究很多[2-7]。在此基礎(chǔ)上，楊葉濤等選擇從城市土地利用景觀格局特征來研究城市擴張，探討多種景觀格局指標與建設(shè)用地增長的相關(guān)性[8-9]。土地利用變化往往會導致景觀結(jié)構(gòu)、功能也發(fā)生了變化。景觀連接度是研究景觀結(jié)構(gòu)和景觀功能互饋的重要途徑[10-11]。吳昌廣等[12]綜述了國內(nèi)外近10年來有關(guān)景觀連接度的概念、度量以及相關(guān)應用，說明了景觀連接度能反映景觀的功能特征，通過探討景觀連接度與景觀要素間的關(guān)系，能為景觀規(guī)劃與生物多樣性保護提供參考。景觀格局發(fā)生變化的原因在于內(nèi)、外驅(qū)動力對景觀組成要素的干擾作用，綜合的干擾作用包括了自然、生物以及人類活動之間復雜的相互作用，其結(jié)果使得景觀穩(wěn)定性及其景觀結(jié)構(gòu)發(fā)生變化，從而使景觀格局發(fā)生變化[13]。那么，在快速城市化進程中，城市人工景觀連接度的提高是否會加快城市景觀格局的變化？城市人工景觀結(jié)構(gòu)與景觀連接度的關(guān)系又是什么，仍需進一步研究。相關(guān)分析揭示了要素之間的相關(guān)程度，回歸分析則是建立反映地理要素之間具體數(shù)量關(guān)系的數(shù)學模型[14]。然而，在未對變量的時間序列的平穩(wěn)性進行檢驗的情況下，直接對其進行回歸分析，極容易產(chǎn)生偽回歸問題，從而導致所建的回歸模型沒有解釋意義[15]。

本研究旨在回答以下問題：1988—2009年，東莞市的城市人工景觀結(jié)構(gòu)與景觀連接度發(fā)生了什么變化？景觀結(jié)構(gòu)與景觀連接度指數(shù)能否說明快速城市化的發(fā)展過程？城市人工景觀結(jié)構(gòu)與景觀連接度在時間序列上是不斷變化的，它們之間存在什么樣關(guān)系？是否存在著長期均衡及短期動態(tài)因果關(guān)系？能否從統(tǒng)計上考察這種關(guān)系是單向的還是雙向的？即主要是一個變量的過去行為在影響另一個變量的當前行為呢？還是雙方的過去行為在互相影響著對方的當前行為？

1 研究區(qū)域

東莞市位于廣東省中南部，珠江口東岸，東江下游的珠江三角洲，毗鄰港澳，處于廣州至深圳、香港經(jīng)濟走廊中間。西北距廣州59 km，東南距深圳99 km，距香港140 km。地處東經(jīng)113°31′—114°15′，北緯22°39′—23°09′。東西長約70.45 km，南北寬約46.8 km，全市陸地總面積2465 km2，海域面積為150 km2。地勢東南高、西北低，地貌以丘陵臺地、沖積平原為主，如圖1所示。

圖1 東莞市區(qū)位示意圖

20世紀80年代以來，在外源性經(jīng)濟的帶動下，東莞市成為珠三角地區(qū)外商投資的中心之一，逐步形成了以制造業(yè)為主的產(chǎn)業(yè)結(jié)構(gòu)模式，吸引了大量外來人口，推動了工業(yè)化進程與城市建設(shè)，從而促進城市快速擴張。本研究以東莞市城市人工用地景觀為研究對象，研究所需東莞市行政區(qū)劃矢量數(shù)據(jù)與1988—2009年間11期Landsat TM遙感影像解譯的土地利用/覆被圖，包括城市用地、開發(fā)區(qū)、水體、農(nóng)田、林地等9種土地利用類型。根據(jù)東莞市城市景觀特征及研究目的，將城市用地和開發(fā)區(qū)兩種用地類型定義為城市人工景觀，即城市景觀的基底，由建筑物群體和硬化地面構(gòu)成的景觀主體。

2 研究方法

2.1 研究方法景觀指數(shù)的選取

景觀指數(shù)是高度濃縮景觀格局信息，反映其結(jié)構(gòu)組成和空間配置某些方面特征的簡單定量指標，即景觀組成單元的類型、多樣性及其空間關(guān)系。如景觀中不同生態(tài)系統(tǒng)(或土地利用類型)的面積、形狀和豐富度。只有針對特定的生態(tài)過程或特定的地理位置，才能凸顯景觀格局指數(shù)的含義與意義[16-17]。描述景觀格局特征的景觀指數(shù)很多，且相關(guān)性較高，本研究綜合研究目的與研究區(qū)的實際情況，選取了兩類景觀指數(shù)。一類是景觀結(jié)構(gòu)指數(shù)，選取斑塊類型面積(CA)、景觀形狀指數(shù)(LSI)、聚集度指數(shù)(AI)、蔓延度指數(shù)(CONTAG)與香農(nóng)多樣性指數(shù)(SHDI)5個指標來描述城市人工景觀結(jié)構(gòu)；另一類是景觀連接度指數(shù)，包括整體連通性指數(shù)(IIC)與可能連通性指數(shù)(PC)。上述景觀指數(shù)的定義、計算公式、閾值以及其生態(tài)學意義可參見有關(guān)文獻[10,12,17]，指數(shù)計算均基于Fragstats4.1軟件和Conefor Sensinode2.2軟件。

2.2 協(xié)整性與因果關(guān)系分析原理

傳統(tǒng)的計量經(jīng)濟學在建立模型的時候，要求隨機過程必須是平穩(wěn)的序列，如果序列是不穩(wěn)的，就會產(chǎn)生“虛假回歸”，致使按照傳統(tǒng)統(tǒng)計推斷得出的結(jié)論嚴重失誤[15,18]。非平穩(wěn)時間序列建模的基本思想是：將理論和數(shù)據(jù)信息有效結(jié)合，從經(jīng)濟變量的數(shù)據(jù)中所顯示的關(guān)系出發(fā)，確定模型包含的變量和變量之間的關(guān)系[19]。協(xié)整概念是一個強有力的概念，因為協(xié)整允許我們刻畫兩個或者多個序列之間的平衡或平穩(wěn)關(guān)系。對于每一個序列單獨來說可能是非平穩(wěn)的，這些序列的矩，如均值、方差和協(xié)方差隨時間而變化，而這些時間序列的線性組合序列卻可能有不隨時間變化的性質(zhì)[20]。研究變量之間的協(xié)整關(guān)系，就等于研究變量之間的定量規(guī)律。其次在建立模型前，對變量之間的協(xié)整關(guān)系進行檢驗，證明它們是協(xié)整的，那么所建立的回歸模型則可以避免偽回歸[21]。協(xié)整性與因果關(guān)系分析主要分為4個步驟。

(1)單位根檢驗

檢查序列平穩(wěn)性的標準方法是單位根檢驗方法。本文采用的是ADF (Augmented Dickey-Fuller test)檢驗方法，通過在回歸方程右邊加入因變量yt的滯后差分項來控制高階序列相關(guān)，一般形式表示為[20]：

(2)協(xié)整檢驗

協(xié)整檢驗方法的一般方法有EG法與Johansen法。本文采用的是EG檢驗法，這種方法是對回歸方程的殘差進行單位根驗證。因此，檢驗一組變量(景觀結(jié)構(gòu)指數(shù)與景觀連接度指數(shù))之間是否存在協(xié)整關(guān)系等價于檢驗回歸方程的殘差序列是否是一個平穩(wěn)序列[20,22]。

(3)誤差修正模型的建立

ECM(Error Correction Model)即誤差修正模型，一般的模型是自回歸分布滯后模型，記為：

yt=β0+β1xt+β2yt-1+β3xt-1+εt

經(jīng)變形整理后得誤差修正模型：

(4)序列間因果關(guān)系分析與檢驗

Granger因果檢驗這種分析方法是從統(tǒng)計角度來確定變量間的因果關(guān)系。格蘭杰因果關(guān)系檢驗不是檢驗邏輯上的因果關(guān)系，而是看變量間的先后順序，是否存在一個變量的前期信息會影響到另一個變量的當期[23]。其實質(zhì)上是檢驗一個變量的滯后變量是否可以引入到其他變量方程中，一個變量如果受到其他變量的滯后影響，則稱他們具有Granger因果關(guān)系。對兩變量Y與X，格蘭杰因果關(guān)系檢驗要求估計以下回歸[22]：

格蘭杰檢驗是通過受約束的F檢驗完成的。如果針對X不是Y的格蘭杰原因這一假設(shè)，即針對式中滯后項前的參數(shù)整體為零的假設(shè)，分別做包含和不包含X滯后項的回歸，記前者的殘差平方和為RSSU，后者的殘差平方和為RSSR，再計算F統(tǒng)計量[22]：

3 結(jié)果與分析

3.1 城市景觀結(jié)構(gòu)與連接度指數(shù)變化特征

1988年至2009年間，東莞市城市景觀格局發(fā)生明顯變化，整體特征如表1所示：從景觀結(jié)構(gòu)指數(shù)CA值的變化情況上看，城市人工景觀的斑塊類型面積占全市陸地總面積由1990年的4.04%、1999年的22.46%、2003年的35.78%上升到2009年的48.07%，說明了東莞市近20年人工景觀用地規(guī)模不斷擴大，城市人工景觀特征逐步向城市外圍自然及半自然景觀滲透的過程，城市空間擴展顯著。

表1 景觀結(jié)構(gòu)指數(shù)與景觀連接度指數(shù)

從斑塊類型水平的尺度上看，LSI與AI值所描述的對象是城市人工景觀斑塊形態(tài)的變化過程，AI值總體上穩(wěn)步上升，說明城市人工景觀的聚集程度越來越高。而LSI值則在1988年至1997年間不斷增加，1997年后該值逐漸減小的態(tài)勢說明了城市人工景觀的形狀變化是從簡單變復雜，然后再從復雜變回簡單的過程，符合城市空間擴展過程的一般規(guī)律。

從景觀類型水平的尺度上看，CONTAG值與SHDI值所描述的對象是城市整體景觀的變化過程。CONTAG值的大小反映景觀不同斑塊類型的非隨機性或聚集程度。如果一個景觀由許多離散的小斑塊組成，CONTAG的值較小，當景觀中以少數(shù)大斑塊為主或同一類型斑塊高度連接時，CONTAG的值較大；SHDI的值的大小取決于兩個方面的信息：一是斑塊類型的多少(即豐富度)，二是各斑塊類型的面積上分布的均勻程度。CONTAG的值在1988年至1997年間由大變小，1997年后該值逐漸增大，而SHDI值的變化過程卻剛好相反。其原因主要是1988—1997年間，東莞市景觀格局的基底主要是由農(nóng)田、林地等自然及半自然景觀為主，隨著城市建設(shè)與發(fā)展的需求，大量的農(nóng)田、林地逐步轉(zhuǎn)化為城市建設(shè)用地，景觀結(jié)構(gòu)組成的復雜性也趨于增加，從而導致景觀聚集度減小。1997年后，東莞市景觀格局的基質(zhì)逐步轉(zhuǎn)化為城鎮(zhèn)居民用地、交通用地等城市人工景觀，景觀結(jié)構(gòu)組成的復雜性也逐漸減小，從而導致景觀聚集度增大。景觀連接度指數(shù)所描述的是城市人工景觀斑塊之間的連接程度，可綜合測度整個景觀尺度上生境斑塊的連通程度。整體來看，研究期間的IIC與PC值都呈現(xiàn)上升趨勢。特別是從1999年開始，IIC與PC值增加速度顯著，說明城市人工景觀斑塊的連接程度高。

綜合表1與圖2—5中可以看出，東莞市城市人工景觀格局變化可分為4個階段，①1988—1990年期間，東莞市的城市人工景觀總量小，截至1990年，城市人工景觀占地約157 km2，總體分布零散，主要集中在莞城區(qū)，如圖2所示，屬城市發(fā)展初期。LSI值與AI值偏小，說明城市人工景觀斑塊聚集程度較低。②1990—1999年期間，東莞市城市人工景觀面積總量迅速增大，新增城市用地景觀以年平均44.08 km2的數(shù)量增加，截至1999年，城市人工景觀占地約554 km2；從空間上看，城市人工景觀的“骨架”逐步形成，分布相對均勻，新增的城市人工景觀主要集中在各個鎮(zhèn)區(qū)的城鎮(zhèn)中心及其周邊地區(qū)，形成多個城市組團，如虎門、長安組團等，如圖3所示。LSI值與AI值增大，說明城市人工景觀斑塊聚集程度逐步上升。③1999—2005年期間，東莞市城市人工景觀面積總量持續(xù)增長，新增城市用地景觀以年平均76.85 km2的數(shù)量增加。截至2005年，城市人工景觀占地約1014.81 km2，新增的城市人工景觀主要分布在前期的城市人工景觀“骨架”的周邊地區(qū)，以厚街鎮(zhèn)為例，如圖4紅色箭頭所示，城市用地空間擴展類型主要以外延與蔓延為主。AI值繼續(xù)增大而LSI值開始下降，說明東莞市城市人工景觀整體格局已逐步趨于穩(wěn)定。④2005—2009年期間，東莞市城市人工景觀面積總量持續(xù)增長，但其增長速度放緩，新增城市用地面積為年平均42.52 km2，截至1999年，城市人工景觀占地約1185 km2。以塘廈鎮(zhèn)為例，如圖5中紅色箭頭所示，城市用地空間擴展類型主要以填充為主。

圖2 1988—1990城市人工景觀

圖3 1990—1999城市人工景觀

圖4 1999—2005城市人工景觀(厚街鎮(zhèn))

圖5 2005—2009城市人工景觀(塘廈鎮(zhèn))

3.2 平穩(wěn)性檢驗

對7個原序列AI、CA、CONTAG、IIC、LSI、PC和SHDI進行單位根檢驗，運用ADF檢驗方法(選擇包含截距項和趨勢項對原序列分別作ADF檢驗)，結(jié)果如表2所示，LSI與SHDI為平穩(wěn)序列，AI、CA、CONTAG、IIC與PC均為非平穩(wěn)序列。

對非平穩(wěn)的時間序列AI、CA、CONTAG、IIC與PC進行數(shù)據(jù)預處理，對數(shù)變換不會影響數(shù)據(jù)之間關(guān)系的研究，將各序列對數(shù)化，記為lnAI、lnCA、lnCONTAG、lnIIC與lnPC。繼續(xù)使用ADF檢驗方法，結(jié)果如表3所示，二階差分序列l(wèi)nAI、lnCA、lnCONTAG、lnIIC與lnPC在1%的顯著水平下是平穩(wěn)的，即均為I(2)序列，滿足協(xié)整檢驗前提。

表2 原始序列單位根檢驗

表3 差分序列單位根檢驗

3.3 協(xié)整檢驗與誤差修正模型

為了確定城市人工景觀結(jié)構(gòu)與景觀連接度指數(shù)之間是否存在協(xié)整關(guān)系，分別對景觀結(jié)構(gòu)指數(shù)(lnAI、lnCA、lnCONTAG)與景觀連接度指數(shù)(lnIIC、lnPC)進行OLS回歸，回歸后得各個模型的殘差序列并對其進行ADF檢驗，分別按照同時含有常數(shù)項和時間趨勢項、有常數(shù)項和沒有時間趨勢項，以及沒有常數(shù)項和時間趨勢項的3個模型來進行單位根檢驗，只要其中一個模型的檢驗結(jié)果拒絕了零假設(shè)，就可以認為該序列式平穩(wěn)的。當3個模型的檢驗結(jié)果都不能拒絕零假設(shè)，則認為該序列式非平穩(wěn)的。檢驗結(jié)果如表4所示。

表4 原始序列單位根檢驗

(C,T,K)表示ADF檢驗式是否包含常數(shù)項、時間趨勢項以及滯后期數(shù)殘差序列的平穩(wěn)性結(jié)果表明，在顯著性5%的水平上，景觀結(jié)構(gòu)指數(shù)lnAI、lnCONTAG與景觀連接度指數(shù)(lnIIC和lnPC)存在協(xié)整關(guān)系，lnCA僅與lnPC存在協(xié)整關(guān)系，lnCA與lnIIC不存在協(xié)整關(guān)系。因此，可以得出本文中的城市人工景觀結(jié)構(gòu)指數(shù)AI、CONTAG與景觀連接度指數(shù)IIC、PC，CA與PC之間存在著長期穩(wěn)定的均衡關(guān)系，它們之間的回歸方程設(shè)定是合理的，故可進一步建立誤差修正模型，來確立它們之間的短期波動關(guān)系。在建立短期模型過程中，將誤差修正項看作一個解釋變量，連同其他反映短期波動的解釋變量一起，建立短期模型，即誤差修正模型。模型表達式如下：

lnAI=5.4420+0.0223lnIIC-0.1911lnAIt-1+0.0022lnIICt-1

R-squared：0.988107F-statistic：166.1644 Durbin-Watson stat：2.025342

lnAI=2.6449+0.0256lnPC+0.4218lnAIt-1-0.0092lnPCt-1

R-squared：0.992451F-statistic：262.9498 Durbin-Watson stat：1.901418

lnCONTAG=1.3614+0.0332lnIIC+0.6684lnCONTAGt-1+0.0447lnIICt-1

R-squared：0.845528F-statistic：10.94734 Durbin-Watson stat：1.854780

lnCONTAG=1.1521-0.0304lnPC+0.7176lnCONTAGt-1+0.0389lnPCt-1

R-squared：0.874812F-statistic：13.97601 Durbin-Watson stat：2.378302

lnCA=-1.3489+0.2446lnPC+1.1163lnCAt-1-0.2927lnPCt-1

R-squared：0.997961F-statistic：978.8289 Durbin-Watson stat：1.834893

3.4 Granger因果關(guān)系檢驗

為了進一步發(fā)現(xiàn)協(xié)整關(guān)系中的因果方向，按照格蘭杰因果關(guān)系檢驗方法，分別對景觀結(jié)構(gòu)指數(shù)(lnAI、lnCA、lnCONTAG)與景觀連接度指數(shù)(lnIIC、lnPC)之間進行Granger因果關(guān)系檢驗。由于該檢驗結(jié)果對滯后期十分敏感，滯后階數(shù)的選擇直接左右格蘭杰因果檢驗的結(jié)果，對于分析的準確性至關(guān)重要。因此本文將選取滯后期1—3進行檢驗，以使得結(jié)果全面準確，最后列出在5%的顯著性水平下，具有格蘭杰因果關(guān)系的結(jié)果，如表5所示。從1階滯后的情況看，lnIIC是lnAI、lnCONTAG的格蘭杰原因，lnPC是lnCONTAG的格蘭杰原因；從2階滯后的情況看，lnAI是lnPC的格蘭杰原因，lnPC與lnCA之間有雙向的格蘭杰因果關(guān)系，lnIIC與lnCONTAG之間有雙向的格蘭杰因果關(guān)系；從3階滯后的情況看，lnPC是lnCONTAG的格蘭杰原因。格蘭杰因果關(guān)系檢驗不是檢驗邏輯上的因果關(guān)系，而是看變量間的先后順序，是否存在一個變量的前期信息會影響到另一個變量的當期。即A是B的格蘭杰成因，則表示A的前期信息會影響到B的當期，說明A過去的信息值有助于預測B。

表5 格蘭杰因果關(guān)系檢驗

4 結(jié)論與討論

研究表明，本文所選取斑塊類型面積(CA)、景觀形狀指數(shù)(LSI)、聚集度指數(shù)(AI)、蔓延度指數(shù)(CONTAG)、香農(nóng)多樣性指數(shù)(SHDI)等5個指標來量化東莞市城市人工景觀結(jié)構(gòu)的動態(tài)變化，選取景觀連接度指數(shù)(IIC與PC)這2個指標來綜合測度整個區(qū)域的城市人工景觀尺度上生境斑塊的連通性程度，上述7個指數(shù)是能較好地反映東莞市城市人工景觀格局在近20年來變化的過程。從城市生態(tài)變化過程的角度上看，東莞市景觀的基質(zhì)由農(nóng)田、林地等自然及半自然景觀逐步轉(zhuǎn)化為城鎮(zhèn)居民用地、交通用地等城市人工景觀。通過分析7個指數(shù)的變化特征，總結(jié)出東莞市快速城市化過程中的4個階段：①1988—1990年期間，自然及半自然景觀(農(nóng)田、林地等)為景觀的基質(zhì)，城市人工景觀斑塊面積小、聚集程度低。②1990—1999年期間，城市人工景觀面積總量迅速增大、聚集程度明顯上升，城市空間擴展類型主要以“外延”與“沿交通線擴展”兩種類型為主，城市人工景觀的“骨架”逐步形成。景觀結(jié)構(gòu)組成的復雜性也趨于增加，自然及半自然景觀(農(nóng)田、林地等)景觀出現(xiàn)萎縮、破碎化和孤島化等現(xiàn)象，城市人工景觀逐步成為景觀的基質(zhì)。③1999—2005年期間，城市人工景觀面積總量持續(xù)增長，向城市外圍自然及半自然景觀滲透，增長的面積相對于上一階段不大，呈散狀分布，屬于“蔓延”的擴展類型。④2005—2009年期間，城市人工景觀面積總量持續(xù)增長，但其增長的速度明顯放緩，城市用地空間擴展類型主要以填充為主。

城市人工景觀結(jié)構(gòu)與景觀連接度指數(shù)在時間序列上是不斷變化的，CA、AI、CONTAG、IIC與PC均為二階平穩(wěn)序列。通過對它們進行協(xié)整分析與格蘭杰因果關(guān)系檢驗，結(jié)果表明：城市人工景觀斑塊類型面積(CA)與可能連通性指數(shù)(IIC)存在協(xié)整關(guān)系，景觀聚集度指數(shù)(AI)、蔓延度指數(shù)(CONTAG)與景觀連通性指數(shù)(IIC與PC)也存在協(xié)整關(guān)系，即這些變量之間確實存在長期均衡關(guān)系。這種均衡關(guān)系意味著該系統(tǒng)中不存在破環(huán)均衡的內(nèi)在機制。如果變量在某時期受到干擾后偏離其長期均衡點，則均衡機制將會在下一期進行調(diào)整以其重新回到均衡狀態(tài)。格蘭杰因果關(guān)系檢驗中，在不同滯后期的情況下，大部分的結(jié)果表明景觀連通性指數(shù)(IIC與PC)是景觀結(jié)構(gòu)指數(shù)(CA、AI、CONTAG)的格蘭杰成因，說明景觀連通性指數(shù)信息有助于預測景觀結(jié)構(gòu)指數(shù)的變化。從另一個角度上看，城市人工景觀連接度指數(shù)的大小取決于各大城市組團之間城鎮(zhèn)與交通路網(wǎng)的發(fā)展。城市交通路網(wǎng)是組成城市的骨架，各種不同大小的城市組團依存這些骨架。因此，通過大力發(fā)展各大城市組團周邊小城鎮(zhèn)，以及完善各城市組團的交通路網(wǎng)系統(tǒng)，將有助于優(yōu)化城市空間結(jié)構(gòu)與格局。

[1] 宋治清, 王仰麟. 城市景觀及其格局的生態(tài)效應研究進展. 地理科學進展, 2004, 30(2): 97-106.

[2] 龔建周, 夏北成. 1990年以來廣州市土地覆被景觀的時空梯度分異. 地理學報, 2007, 62(2): 181-190.

[3] 王興芳, 王云鵬, 李巖, 李慧, 陳永品, 沈奎. 廣州東進擴張下的蘿崗區(qū)土地利用與景觀動態(tài)分析. 城市環(huán)境與城市生態(tài), 2010, 23(6): 5-8.

[4] 俞龍生, 符以福, 喻懷義, 李志琴. 快速城市化地區(qū)景觀格局梯度動態(tài)及其城鄉(xiāng)融合區(qū)特征——以廣州市番禺區(qū)為例. 應用生態(tài)學報, 2011, 22(1): 171-180.

[5] 高楊, 吳志峰, 劉曉南, 張景華. 珠江三角洲景觀空間格局分析. 熱帶地理, 2008, 28(1): 26-31.

[6] Aguilera F, Valenzuela L M, Botequilha-Leit?o A. Landscape metrics in the analysis of urban land use patterns: A case study in a Spanish metropolitan area. Landscape and Urban Planning, 2011, 99(3/4): 226-238.

[7] Wu J G, Jenerette G D, Buyantuyev A, Redman, C L. Quantifying spatiotemporal patterns of urbanization: The case of the two fastest growing metropolitan regions in the United States. Ecological Complexity, 2011, 8(1): 1-8.

[8] 韋薇, 張銀龍, 趙兵, 王浩. 快速城市化進程中城市擴張對景觀格局分異特征的影響. 生態(tài)環(huán)境學報, 2011, 20(1): 7-12.

[9] 楊葉濤, 龔建雅, 周啟鳴, 王迎迎. 土地利用景觀格局對城市擴張影響研究. 自然資源學報, 2010, 25(2): 320-329.

[10] 吳健生, 劉洪萌, 黃秀蘭, 馮喆. 深圳市生態(tài)用地景觀連通性動態(tài)評價. 應用生態(tài)學報, 2012, 23(9): 2543-2549.

[11] Forman R T T. Some general principles of landscape and regional ecology. Landscape Ecology, 1995, 10(3): 133-142.

[12] 吳昌廣, 周志翔, 王鵬程, 肖文發(fā), 滕明君. 景觀連接度的概念、度量及其應用. 生態(tài)學報, 2010, 30(7): 1903-1910.

[13] 齊楊, 鄔建國, 李建龍, 于洋, 彭福利, 孫聰. 中國東西部中小城市景觀格局及其驅(qū)動力. 生態(tài)學報, 2013, 33(1): 275-285.

[14] 徐建華. 計量地理學. 北京: 高等教育出版社, 2006: 122-164.

[15] 葉浩, 濮勵杰. 江蘇省耕地面積變化與經(jīng)濟增長的協(xié)整性與因果關(guān)系分析. 自然資源學報, 2007, 22(5): 766-774.

[16] 張秋菊, 傅伯杰, 陳利頂. 關(guān)于景觀格局演變研究的幾個問題. 地理科學, 2003, 23(3): 264-270.

[17] 鄔建國. 景觀生態(tài)學：格局、過程、尺度與等級. 北京: 高等教育出版社, 2007. 60-69.

[18] 李永樂, 吳群. 中國經(jīng)濟增長與耕地資源數(shù)量變化階段性特征研究——協(xié)整分析及Granger因果檢驗. 長江流域資源與環(huán)境, 2011, 20(1): 33-39.

[19] 任英華. Eviews應用實驗教程. 長沙: 湖南大學出版社, 2008: 167-209.

[20] 高鐵梅. 計量經(jīng)濟分析方法與建模：EViews應用及實例. 北京: 清華大學出版社, 2006: 78-92.

[21] 于俊年. 計量經(jīng)濟學：EViews的使用. 北京: 對外經(jīng)濟貿(mào)易大學出版社, 2006: 79-83.

[22] 李子奈, 潘文卿. 計量經(jīng)濟學. 北京: 高等教育出版社, 2005: 267-323.

[23] 劉巍, 陳昭. 計量經(jīng)濟學軟件：EViews操作簡明教程. 廣州: 暨南大學出版社, 2009: 22-26

Research on the cointegration and causality between urban artificial landscape structure change and landscape connectivity in Dongguan City

FENG Zhixin1,2, CHEN Yingbiao2,*, QIAN Qinglan2, WANG Shuaishuai2

1GuangdongProvincialAcademyofEnvironmentalScience,Guangzhou510031,China2SchoolofGeographicalSciences,GuangzhouUniversity,Guangzhou510006,China

Characteristics of urban artificial landscape have permeated into natural and semi-natural landscapes. Prior studies focused on landscape pattern analysis and description in rapid urbanization areas but neglect the dynamic changes between structures and connectivity of artificial urban landscape and their internal interactions. We measured two landscape indices (landscape structures and landscape connectivity) and presented an integrated approach (the cointegration and causality) to explore the urban expansion of Dongguan, a rapidly growing city in southern China, since 1988 to 2009. Analyses of seven landscape indices indicates a clear dynamic change in the expansion of Donguan. We divided this expansion into four different periods, that is early natural and semi-natural landscapes dominated urban areas in 1988—1990, followed by rapid ‘extension’ and ‘linear development’ period through 1990—1999, and then urbanization developed steadily period (‘sprawl’) during 1999—2005, and the recent ‘infilling’ stage since 2005—2009 characterized by urban development? The results show that: (1) According to variation characteristics of seven indices, the period is divided into four stages. The first is 1988—1990, in which the natural and semi-natural landscapes such as farmland, woodland were dominant, while urban artificial areas were small and had low density and aggregation. The second period is 1990 —1999, in which urban artificial landscape increased rapidly with high density and aggregation. Then, we named the urban land expansion as ‘extension’ and ‘linear development’ type. The skeleton of the city gradually formed and the urban artificial landscape was in the dominant position, while natural and semi-natural landscapes such as farmland, woodland were small with low density and aggregation. The third period is 1999—2005. The total area of urban artificial landscape grew steadily and we named the urban land expansion as ‘sprawl’ type. The last is 2005—2009. Urban artificial landscape grew further and we named the urban land expansion as ‘infilling’ type. (2) There exist cointegration relationships between class area (CA) and probability index of connectivity (PC), aggregation index (AI), contagion index (CONTAG), and connectivity indices IIC and PC. This indicates that there is a long-term equilibrium relationship among these variables and thereby no internal mechanism would break down the equilibrium of urban system. (3) Granger causality test shows that landscape connectivity indices (IIC and PC) are the Granger Causes of landscape structure indices (CA, AI, CONTAG), which suggests the landscape connectivity indices can be used to predict the change of landscape structures.

urban expansion; landscape structure; landscape connectivity; Dongguan

國家科技支撐計劃課題(2012BAH32B03)；國家自然科學基金項目(41071078)；教育部人文社科規(guī)劃基金項目(11YJAZH016)；住房和城鄉(xiāng)建設(shè)部科學技術(shù)項目(2012-K8-42)資助

2013-11-05;

2014-12-04

10.5846/stxb201311052673

*通訊作者Corresponding author.E-mail: gzhuchenyb@126.com

馮志新,陳穎彪,千慶蘭,王帥帥.東莞市人工景觀結(jié)構(gòu)與連接度的協(xié)整性及因果關(guān)系.生態(tài)學報,2015,35(14):4920-4930.

Feng Z X, Chen Y B, Qian Q L, Wang S S.Research on the cointegration and causality between urban artificial landscape structure change and landscape connectivity in Dongguan City.Acta Ecologica Sinica,2015,35(14):4920-4930.