石培基，王祖靜，劉春芳

(西北師范大學(xué)地理與環(huán)境科學(xué)學(xué)院，蘭州 730070)

城市化過(guò)程是不連續(xù)的、具有梯度性的，其必然會(huì)導(dǎo)致區(qū)域地表覆被格局的改變[1]，土地覆被時(shí)空演化在LUCC[2-4]、城鎮(zhèn)擴(kuò)展[5-8]、區(qū)域景觀格局[9-12]及城市熱環(huán)境[13-14]、水環(huán)境[15]等方面具有廣泛應(yīng)用。在一定范圍內(nèi)生態(tài)系統(tǒng)也隨著地表覆被的改變而改變，這種改變?cè)谏鷳B(tài)不敏感地區(qū)引起的效果甚微[16]，但是在生態(tài)脆弱地區(qū)，比如石羊河流域，人類活動(dòng)的干擾會(huì)導(dǎo)致整個(gè)流域生態(tài)系統(tǒng)的紊亂，引發(fā)一系列的環(huán)境為題[17]。為了更好的指導(dǎo)城市建設(shè)和流域生態(tài)保護(hù)，找出近十多年的地表覆被變化規(guī)律，選擇流域地表覆被空間格局演化作為研究生態(tài)環(huán)境效應(yīng)的典型代表，將流域土地覆被的時(shí)間與空間研究范式結(jié)合起來(lái)，并以此立題，作為本研究的出發(fā)點(diǎn)。

流域生態(tài)環(huán)境問(wèn)題是目前城市化與流域可持續(xù)發(fā)展普遍關(guān)注的問(wèn)題[18]，本文以Ridd于1995年提出的植被-不透水面-土壤3種基本要素構(gòu)成的城市地表覆被系統(tǒng)[19]為基礎(chǔ)，結(jié)合流域1998—2010年NDVI、夜間燈光數(shù)據(jù)和土壤含量，運(yùn)用改進(jìn)的計(jì)算人居環(huán)境指數(shù)(Human_Settlements_Index，HSI)算法[20]計(jì)算土地人工覆被指數(shù)(LCI)。然后對(duì)其進(jìn)行時(shí)間和空間尺度上的研究，利用小波分析、空間變差函數(shù)和間隙度指數(shù)等方法構(gòu)建測(cè)度模型，分析了石羊河流域地區(qū)土地覆被空間演化的過(guò)程和格局。最后嘗試?yán)没貧w擬合方程找出影響地區(qū)覆被分異和變化的驅(qū)動(dòng)機(jī)制，以此反映近十多年來(lái)石羊河流域地表演化過(guò)程。

1 研究區(qū)概況

石羊河流域位于河西走廊東段，祁連山北麓，介于東經(jīng)101°22′—104°16′、北緯36°29′—39°27′之間。以高山、中高山、低山丘陵、沙漠、沖積平原構(gòu)成全流域地貌形態(tài)，地勢(shì)南高北低，自西南向東北傾斜。面積約41600km2，流域行政區(qū)劃包括武威市涼州區(qū)、古浪縣、民勤縣全部及天?？h烏鞘嶺以北部分地區(qū)，金昌市永昌縣和金川區(qū)全部以及張掖市肅南縣部分地區(qū)，共3市7縣(區(qū))。流域生態(tài)環(huán)境脆弱、沙漠化嚴(yán)重，十多年來(lái)“西部大開發(fā)”快速的城鎮(zhèn)化建設(shè)使這些問(wèn)題更加突出。

圖1 研究區(qū)概況

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

2.1 數(shù)據(jù)處理

本研究1998年、2002年、2006年、2010年流域NDVI、夜間燈光數(shù)據(jù)和土壤信息分別來(lái)源于SPOTVGT數(shù)據(jù)、美國(guó)國(guó)家地球物理數(shù)據(jù)中心(NGDC)的DMSP-OLS和中國(guó)科學(xué)院計(jì)算機(jī)網(wǎng)絡(luò)信息中心國(guó)際科學(xué)數(shù)據(jù)服務(wù)平臺(tái)的Landsat TM影像，流域基礎(chǔ)數(shù)據(jù)主要來(lái)源于國(guó)家自然科學(xué)基金委員會(huì)“中國(guó)西部環(huán)境與生態(tài)科學(xué)數(shù)據(jù)中心”(http://westdc.westgis.ac.cn)。

統(tǒng)一把流域基礎(chǔ)數(shù)據(jù)集的坐標(biāo)系轉(zhuǎn)換為西安1980坐標(biāo)系，由于轉(zhuǎn)換后的NDVI與DMSP-OLS圖像均為1km×1km空間分辨率，所以無(wú)需進(jìn)行重采樣，用流域邊界裁切圖像。在ENVI中將4個(gè)年份7—8月SPOT Vegetation的6期數(shù)據(jù)轉(zhuǎn)換為真實(shí)NDVI值，并提取最大值；再提取OLS最大值和最小值，根據(jù)公式對(duì)其進(jìn)行歸一化處理：

(1)

式中，OLSnor為歸一化燈光指數(shù)，OLS是燈光指數(shù)實(shí)際值，OLSmax和OLSmin分別為燈光指數(shù)最大值和最小值。結(jié)合Lu等提出的提取人居環(huán)境指數(shù)(HSI)空間信息算法[20]：

(2)

式中，NDVImax為7—8月份NDVI最大值，解釋為若一個(gè)像元中居住區(qū)所占比重大，那么導(dǎo)致NDVImax較低、DMSP-OLS較高，所以HSI指數(shù)偏高。改進(jìn)的計(jì)算土地覆被指數(shù)(LCI)，以土地人工覆被綜合表征，修改后公式如下：

(3)

式中，SOIL為土壤比例指數(shù)，為解譯的4期16景TM影像(均在7月或8月份)土地利用類型估算土壤比例，然后采用極差標(biāo)準(zhǔn)化方法，對(duì)其進(jìn)行標(biāo)準(zhǔn)化處理。同式(2)，若一個(gè)像元中人居環(huán)境、即不透水地表所占比重大，那么導(dǎo)致SOIL較低、HSI較高，所以LCI指數(shù)偏高。

2.2 研究方法

2.2.1 小波分析

小波變換是時(shí)間(空間)頻率的局部化分析，它通過(guò)伸縮平移運(yùn)算對(duì)信號(hào)(函數(shù))逐步進(jìn)行多尺度細(xì)化，最終達(dá)到高頻處時(shí)間細(xì)分，低頻處頻率細(xì)分，能自動(dòng)適應(yīng)時(shí)頻信號(hào)分析的要求，從而可聚焦到信號(hào)的任意細(xì)節(jié)。反映景觀要素在不同尺度上的特征及其相互關(guān)系，揭示空間格局的多尺度和等級(jí)結(jié)構(gòu)特征，計(jì)算公式為：

(4)

式中，W(a,xj)為景觀要素在空間尺度為a時(shí)的波值，xj為窗函數(shù)的中值，g(xi-xj)為窗函數(shù)，對(duì)地形因子以及土地覆被指數(shù)的空間采樣序列進(jìn)行小波方差和小波一致性分析。通過(guò)誤差分析和相關(guān)分析檢驗(yàn)?zāi)Ｐ偷哪M效果，統(tǒng)計(jì)參數(shù)包括均方根偏差(RMSE)、平均偏差(MBE)、相對(duì)偏差(RE)[21]。

2.2.2 空間變差函數(shù)

空間變差函數(shù)也叫半變異函數(shù)，是地統(tǒng)計(jì)分析的特有函數(shù)。區(qū)域化變量Z(x)在點(diǎn)x和x+h出的值Z(x)和Z(x+h)差得方差的一半稱為區(qū)域化變量Z(x)的半變異函數(shù)，記為r(h)，2r(h)稱為變異函數(shù)。根據(jù)計(jì)算有：

(5)

半變異函數(shù)曲線圖反映了一個(gè)采樣點(diǎn)與其相鄰采樣點(diǎn)的空間關(guān)系(圖2)。理論上，當(dāng)采樣點(diǎn)間的距離為0時(shí)，半變異函數(shù)值應(yīng)為0；但由于存在測(cè)量誤差和空間變異，使得兩采樣點(diǎn)非常接近時(shí)，它們的半變異函數(shù)值不為0，即存在塊金值(C0)。當(dāng)采樣點(diǎn)間的距離h增大時(shí)，半變異函數(shù)r(h)從初始的塊金值達(dá)到一個(gè)相對(duì)穩(wěn)定的常數(shù)時(shí)，該常數(shù)稱為基臺(tái)值(C0+C)。當(dāng)半變異函數(shù)的取值由初始的塊金值達(dá)到基臺(tái)值時(shí)，采樣點(diǎn)的間隔距離稱為變程(a)。在變程范圍內(nèi)，樣點(diǎn)間的距離越小，其相似性，即空間相關(guān)性越大。當(dāng)h>R時(shí)，區(qū)域化變量Z(x)的空間相關(guān)性不存在[22]。

圖2 變差函數(shù)理論方差圖

2.2.3 間隙度指數(shù)

間隙度是Mandelbrot于1983年在《自然界分形幾何》一書中最早提出，被廣泛地應(yīng)用于景觀質(zhì)地分析和景觀的分形特征分析中，間隙度指數(shù)(Lacunarity Index)用于分析地表景觀空間形態(tài)演化的過(guò)程。間隙度分析方法主要包括移動(dòng)窗口法和查分盒子法，前者是對(duì)土地覆被總體特征及分布的一種度量，后者則用于識(shí)別土地覆被空間分布的局部異質(zhì)性[23]。

全局間隙度指數(shù)計(jì)算公式為：

(6)

式中，M為邊長(zhǎng)為r的計(jì)數(shù)窗口中土地覆被柵格變化的數(shù)量，Q(M,r)為計(jì)數(shù)窗口中柵格變化的比例。

局部間隙度指數(shù)計(jì)算公式為：

(7)

式中，Mw=∑i,jnr(i,j)，nr(i,j)=v-μ+1，Mw為局部區(qū)域(w×w)內(nèi)邊長(zhǎng)為r滑動(dòng)窗口中土地覆被變化的數(shù)量，nr(i,j)是局部區(qū)域i和j的列相對(duì)高度，v為具有最大象元值的的立方體盒子，μ為具有最小象元值的立方體盒子[23]。

本文綜合以上方法，首先從直線樣帶的基礎(chǔ)上，用小波分析找出影響流域土地覆被演化的因素，然后在面狀范圍內(nèi)用變差函數(shù)分析LCI在時(shí)間點(diǎn)上的變化格局分布，最后從時(shí)空動(dòng)態(tài)角度用間隙度指數(shù)研究LCI在時(shí)間段上的變化過(guò)程。

3 過(guò)程分析

3.1 土地覆被分布格局

3.1.1 總體特征

由于所選數(shù)據(jù)年份的限制，以每隔4a為1個(gè)時(shí)間點(diǎn)，利用流域3種因子合成1998年、2002年、2006年、2010年土地人工覆被指數(shù)(LCI)，綜合表征土地覆被。用縣區(qū)級(jí)行政區(qū)劃裁切生成的柵格，得出流域各縣區(qū)LCI(表1)。植被部分主要是利用NDVI體現(xiàn)的，由于采用1km×1km分辨率，干旱內(nèi)陸河流域城市規(guī)模較小，2010年武威市建成區(qū)為16.75km2，金昌市建成區(qū)為55.51km2，占流域比重分別為0.04%、0.13%，其中城市植被更小，所以本文沒(méi)有識(shí)別，植被部分僅指大區(qū)域耕地、林地、草地等。

表1 1998—2010年流域土地覆被分布特征

圖3 1998—2010年流域土地覆被空間分布圖

以土地覆被指數(shù)LCI顯示石羊河流域地表蓋度的空間格局，如表1，圖3所示。民勤LCI均值在流域中最高，這是由于LCI的計(jì)算是以NDVI為基數(shù)，其被巴丹吉林和騰格里沙漠包圍，NDVI較低；天柱地處祁連山區(qū)，LCI均值最低。1998—2010年間流域LCI呈現(xiàn)先增后減的趨勢(shì)，1998—2002年均值從0.3991增長(zhǎng)到0.4282，平均年增長(zhǎng)率為0.0073，增幅較大；2002—2010年的兩個(gè)階段均值從0.4178持續(xù)減小到0.4098，平均每年減少0.54%。除金川外各縣區(qū)LCI演變過(guò)程和流域基本相似，由于金川為工礦型城市，建設(shè)集中，在最近十幾年間建設(shè)用地持續(xù)擴(kuò)張，所以其均值不斷增加。從其余6個(gè)縣區(qū)的均值變化中看出，在1998—2002年間增幅最大的為天柱，年均增長(zhǎng)2.81%；增幅最小的為肅南，年均增長(zhǎng)0.45%。2002—2010年間減幅最大的是天柱，年均減少1.96%；減幅最小的是涼州，年均減少0.41%。

4.1.2 地區(qū)差異

為表現(xiàn)出地表覆被組分在空間演變上的地區(qū)差異及幅度，采用剖面分析方法進(jìn)行組分的空間分異分析。布設(shè)2條沿主要交通和水域演變方向的剖面線(圖1)，比較LCI沿2個(gè)方向的梯度變化。圖4中N-S、E-W兩個(gè)方向?yàn)榈湫蜆訋?，囊括了所有地表類型，?km×1km為單元，分析流域LCI在相同位置上不同時(shí)間內(nèi)的變化特征。

可以看出樣帶1(N-S方向)的變化不如樣帶2(E-W方向)明顯，樣帶2大體沿G30連霍高速橫貫整個(gè)流域，西部處于金川和永昌的城市交界處、中部位于涼州城市建設(shè)中心區(qū)，由于西部大開發(fā)的政策，城市增長(zhǎng)迅速，所以LCI變化巨大，并且隨著時(shí)間的推移明顯呈遞增趨勢(shì)，其中涼州城市值中心最高，2010年為1.1349，2002—2006年增幅最大達(dá)209%；東部位于古浪和天柱直抵烏鞘嶺，地表蓋度較低，由于國(guó)家實(shí)行退耕還林草等，使得2002—2010年間LCI呈遞減趨勢(shì)，說(shuō)明地區(qū)地表植被覆蓋度越來(lái)越高，環(huán)境越來(lái)越好。樣帶1沿石羊大河走勢(shì)縱貫流域，北部處于兩大沙漠的分界線上，植被覆蓋度低，所以LCI較高，并且在時(shí)間上幾乎沒(méi)有變化；中部位于211省道沿線涼州至民勤段，由于城市建設(shè)的原因，1998—2010年間民勤及周邊LCI逐漸增大，形成的幾個(gè)低谷是民勤和涼州耕地；南部位于肅南和天柱邊界上的水源涵養(yǎng)區(qū)，LCI時(shí)間上變化幅度較小，空間上起伏較大，值較高的地區(qū)說(shuō)明有冰川和永久性積雪，較低說(shuō)明是植被分布，12年間LCI總體呈下降趨勢(shì)，可以看出山區(qū)植被增多，流域荒漠化程度改善。

圖4 4個(gè)時(shí)期不同方向上的LCI空間分異坡面圖

圖5 不同方向上LCI與DEM的小波一致性分析

為找出影響土地覆被空間分布的因素，以流域2條樣帶的DEM和2010年LCI為例，進(jìn)行小波一致性分析(圖5)?？梢钥闯?，小尺度(32km以下)和中尺度(32—64km)上地形變化對(duì)LCI的作用很小，只有零散的圓圈、箭頭分布；當(dāng)達(dá)到閥值(64km)以后，地形對(duì)LCI的作用逐漸明顯，存在空間位置上的一致性，將樣帶推演到全流域，反映出大尺度地形因子是形成LCI格局分異的主要控制因子，即石羊河流域由南向北“山區(qū)-平原-荒漠”大范圍上的依次分布，是形成LCI空間格局的一大原因。

3.2 土地覆被演變過(guò)程

3.2.1 時(shí)間格局演變

定義采樣步長(zhǎng)為5km，將LCI賦予重采樣的格網(wǎng)中心點(diǎn)計(jì)算實(shí)驗(yàn)變差函數(shù)，得出石羊河流域4a的方差圖及變差擬合曲線(表2，圖6)。可以看出，通過(guò)最小二乘法擬合的模型在前3a相同，都是指數(shù)模型，2010年為球狀模型，決定系數(shù)都達(dá)到0.9以上，說(shuō)明擬合精度很好。

變程(a)的不斷增大，特別是在2010年激增為22.5km，說(shuō)明流域地表覆被在區(qū)域內(nèi)部相互影響，即城市建設(shè)和退耕還林等人為影響因素較高。將基臺(tái)值(C0+C)定義為系統(tǒng)方差，基臺(tái)值和塊金值(C0)的差值為結(jié)構(gòu)方差，塊金值不斷變大，說(shuō)明表征LCI變化主導(dǎo)因素的空間相關(guān)性在增大。結(jié)構(gòu)方差與系統(tǒng)方差的比值即塊金系數(shù)(C/[C0+C])，在4個(gè)典型年份中不斷減小，在采樣誤差等同的情況下，說(shuō)明象元內(nèi)部變異性減小、外部差異變大，即LCI區(qū)域間變化增大，由于人為因素對(duì)流域地表覆被空間差異的影響越來(lái)越大，相反自然因素的相對(duì)重要性有所減弱。

LCI的異向方差3D擬合模型顯示出變差函數(shù)曲線在各個(gè)方向上都具有規(guī)律性，反映了流域LCI隨距離變化的空間自相關(guān)性，總體上擬合模型呈現(xiàn)四周高中間低的“凹型”漏斗結(jié)構(gòu)，并且隨著時(shí)間的推移，這種結(jié)構(gòu)越來(lái)越明顯，說(shuō)明石羊河流域LCI高低值在區(qū)域上發(fā)生了集聚，而且高值區(qū)的值越來(lái)越高，低值區(qū)的值越來(lái)越低。在0°(東—西)和90°(南—北)方向上、45°(石羊河沿線)和135°(312國(guó)道武威、金昌段)方向上的擬合效果均較好，變差函數(shù)曲線的變化趨勢(shì)基本相同，決定系數(shù)都大于0.9，也說(shuō)明了LCI在區(qū)域間各個(gè)方向上的空間自相關(guān)作用較好。

表2 流域LCI的空間變異特征值

圖6 流域土地覆被梯度變差演化

3.2.2 空間格局分異

計(jì)算流域LCI 1998—2002年、2002—2006年、2006—2010年的變化量，提取大于0的組分，利用差分盒子法分別測(cè)算3個(gè)階段LCI變化量的局部間隙度指數(shù)(圖7)。1998—2002年局部間隙度指數(shù)高值區(qū)分布零散，變化量在全流域沒(méi)有形成集聚簇；低值區(qū)分布在古浪、金川中心城區(qū)的周邊地區(qū)，說(shuō)明這部分地區(qū)LCI空間變化量很小，區(qū)域協(xié)同性較強(qiáng)。2002—2006年間隙度指數(shù)變化值較小，高值區(qū)主要分布在祁連山沿線，分散化程度減弱，說(shuō)明這段時(shí)間這片區(qū)域內(nèi)部分地區(qū)的地表覆被變化較平凡，并且LCI高低分化明顯，這主要是由于地形復(fù)雜和土地利用類型較多，寒漠和林草交錯(cuò)分布；另外在金川和涼州城區(qū)形成低值集聚區(qū)。2006—2010年在民勤的外圍兩大沙漠和地區(qū)祁連山區(qū)形成了高值區(qū)，表明沙漠環(huán)境也在變化；在流域內(nèi)5大縣(區(qū))城市建設(shè)中心區(qū)形成了低值集聚區(qū)，說(shuō)明建成區(qū)面積不斷擴(kuò)大。綜合來(lái)看，1998—2010年間流域城鎮(zhèn)空間擴(kuò)展速度加劇，同時(shí)部分地區(qū)地表覆被也在發(fā)生變化。

圖7 不同時(shí)期流域LCI變化的局部間隙度分析

運(yùn)用移動(dòng)窗口法測(cè)算3個(gè)時(shí)間段內(nèi)LCI變化量的全局間隙度指數(shù)(滑窗尺度變化范圍為2—39象元)，生成雙自然對(duì)數(shù)曲線(圖8)?？梢钥闯?，2002—2006年、2006—2010年ln(全局間隙度指數(shù))在初始階段都小于1998—2002年，但是隨著ln(滑窗尺度)的增大，逐漸超過(guò)了1998—2002年，并且所用滑窗尺度越來(lái)越小，說(shuō)明局部間隙度指數(shù)高值和低值零散分布區(qū)在減少，向著集聚區(qū)發(fā)展，即小區(qū)域內(nèi)部性的LCI變化區(qū)逐漸消失，取而代之出現(xiàn)了一些大范圍的LCI變化區(qū)，區(qū)域異質(zhì)化程度增強(qiáng)。

圖8 不同時(shí)期全局間隙度指數(shù)曲線

4 驅(qū)動(dòng)力分析

從上文分析可知，影響流域地表覆被變化的是人為和自然雙重因素。為了找出空間上導(dǎo)致流域地表覆被變化的因素，人為影響因素的數(shù)據(jù)主要來(lái)自于1998—2010年石羊河流域各縣區(qū)和周邊縣區(qū)的人口和經(jīng)濟(jì)數(shù)據(jù)，來(lái)源于“甘肅、青海、內(nèi)蒙古統(tǒng)計(jì)年鑒”；自然影響因素主要來(lái)自于1998—2010年石羊河流域及周邊氣象站點(diǎn)7、8月份的氣溫和降水?dāng)?shù)據(jù)(包括山丹、永昌、武威、民勤、馬鬃山、松山、景泰、阿拉善右旗、門源、巴音毛道)，來(lái)源于中國(guó)氣象科學(xué)數(shù)據(jù)共享服務(wù)網(wǎng)。

對(duì)驅(qū)動(dòng)因子進(jìn)行Kriging插值，即數(shù)據(jù)空間化，再根據(jù)上文公式計(jì)算這12a的LCI，與同期人為(人口、GDP)、自然(氣溫、降水)因素分別進(jìn)行回歸分析，公式如下：

(8)

從自然因素與流域LCI的相關(guān)性分布圖(圖9)可以看出，相關(guān)性較高的區(qū)域主要分布在祁連山水源涵養(yǎng)區(qū)，相關(guān)性平均值最大，為0.74(氣溫)和0.76(降水)，這部分地區(qū)為林地和草地；另外，高值區(qū)還分布在北部荒漠與綠洲交互區(qū)，表明沙漠和戈壁演化的主要驅(qū)動(dòng)力來(lái)自于自然因素，綜合說(shuō)明了在人跡罕至的高山和荒漠地區(qū)，自然條件是影響流域地表覆被的主要因素。

圖9 1998—2010年LCI與自然因素相關(guān)系數(shù)分布圖

圖10 1998—2010年LCI與人為因素相關(guān)系數(shù)分布圖

5 結(jié)論與討論

本文以NDVI、DMSP-OLS和Landsat TM為基礎(chǔ)數(shù)據(jù)，研究了石羊河流域地表覆被的時(shí)空演化分布格局，發(fā)現(xiàn)整個(gè)流域LCI隨時(shí)間的改變而不斷演化。1998—2010年間，LCI均值經(jīng)歷了由小變大、再由大變小的過(guò)程，其中2002年達(dá)到峰值，之后均值逐漸減小。各縣區(qū)LCI的演化規(guī)律也各不相同，涼州和金川的城市中心區(qū)由于城市不斷擴(kuò)張，所以LCI一直在增大，肅南、天祝高植被覆蓋區(qū)地表蓋度逐漸變好，民勤、古浪等荒漠區(qū)LCI變化不明顯。接著利用典型樣帶和小波分析，印證了這一情況。然后利用實(shí)驗(yàn)變差函數(shù)和間隙度指數(shù)模擬了流域地表覆被正向演化過(guò)程，揭示了區(qū)域性聚集演化規(guī)律逐漸增強(qiáng)，零散性演化分布逐漸減少，即LCI高值和低值分布區(qū)逐漸形成集聚簇，而零散分布區(qū)逐漸減少。

最后根據(jù)流域氣象和人口經(jīng)濟(jì)數(shù)據(jù)找出驅(qū)使LCI演化的規(guī)律分布，可以看出，流域城鎮(zhèn)集中分布的區(qū)域及其周邊地區(qū)LCI變化主要受人口和經(jīng)濟(jì)因素影響，在遠(yuǎn)離城鎮(zhèn)人煙稀少的地區(qū)，LCI變化受氣溫和降水(即自然因素)影響較強(qiáng)，研究結(jié)果為石羊河流域生態(tài)保護(hù)建設(shè)提供了參考。

地表覆被時(shí)空變化是區(qū)域城鎮(zhèn)空間演化的主要因素，研究城鎮(zhèn)演化和地表演化的方法由一定的相通性，所以本文將研究城鎮(zhèn)空間擴(kuò)展的方法引入到土地覆被變化過(guò)程中，由于所選年份比較接近，土地覆被變化不是很明顯；并且流域及周邊的統(tǒng)計(jì)數(shù)據(jù)有限，驅(qū)動(dòng)力分析不夠充分全面。下一步將利用長(zhǎng)時(shí)間序列，結(jié)合高精度遙感影像，選擇流域典型城市進(jìn)一步研究。

:

[1]Li X, Yeh A G. Application of remote sensing for monitoring and analysis of urban expansion a case study of Dongguan. Geographical Research, 1997, 16(4): 56-62.

[2]Lu D S, Weng Q H. Use of impervious surface in urban land-use classification. Remote Sensing of Environment, 2006, 102(1/2): 146-160.

[3]Herrmann S M, Anyamba A, Tucker C J. Recent trends in vegetation dynamics in the African Sahel and their relationship to climate. Global Environmental Change, 2005, 15(4): 394-404.

[4]Feng Y J, Liu Y, Tong X H, Liu M L, Deng S S. Modeling dynamic urban growth using cellular automata and particle swarm optimization rules. Landscape and Urban Planning, 2011, 102(3): 188-196.

[5]Liu X P, Li X, Shi X, Zhang X H, Chen Y M. Simulating land-use dynamics under planning policies by integrating artificial immune systems with cellular automata. International Journal of Geographical Information Science, 2010, 24(5): 783-802.

[6]Che Q J, Duan X J, Guo Y, Wang L, Cao Y H. Urban spatial expansion process, pattern and mechanism in Yangtze river delta. Acta Geographica Sinica, 2011, 66(4): 446-456.

[7]Ma X D, Zhu C G, Ma R H, Pu Y X. Urban spatial growth pattern and its evolution in Suzhou, Jiangsu Province, China. Acta Geographica Sinica, 2008, 63(4): 405-416.

[8]Shi P J, Wang Z J, Li W. Urban spatial expansion process and pattern in Shiyang river basin. Scientia Geographica Sinica, 2012, 32(7): 840-845.

[9]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.

[10]Chen Y M, Li X, Liu X P, Liu Y L. An agent-based model for optimal land allocation (AgentLA) with a contiguity constraint. International Journal of Geographical Information Science, 2010, 24(8): 1269-1288.

[11]Lee W K. A Study on the Application of UrbanSim in the Analysis of Land Use Change in Taipei City [D]. Taipei: National Taiwan University, 2008.

[12]Liu Z H, Wang Y L, Peng J, Xie M M, Li Y. Using ISA to analyze the spatial pattern of urban land cover change: a case study in Shenzhen. Acta Geographica Sinica, 2011, 66(7): 961-971.

[13]Yue W Z, Xu J H. Impact of human activities on urban thermal environment in Shanghai. Acta Geographica Sinica, 2008, 63(3): 247-256.

[14]Yuan F, Bauer M E. Comparison of impervious surface area and normalized difference vegetation index as indicators of surface urban heat island effects in Landsat imagery. Remote Sensing of Environment, 2007, 106(3): 375-386.

[15]Kuang W H, Liu J Y, Lu D S. Pattern of impervious surface change and its effect on water environment in the Beijing-Tianjin-Tangshan metropolitan area. Acta Geographica Sinica, 2011, 66(11): 1486-1496.

[16]Niu Z G, Gong P, Cheng X, Guo J H, Wang L, Huang H B, Shen S Q, Wu Y Z, Wang X F, Wang X W, Ying Q, Liang L, Zhang L N, Wang L, Yao Q, Yang Z Z, Guo Z Q, Dai Y J. Geographical characteristics of China′s wetlands derived from remotely sensed data. Science in China: Series D: Earth Science, 2009, 52(6): 723-738.

[17]Yang Y C, Li J J, Chen F H, Jacquie B, Li R C, Li D, Chang G Y, Li Y C. The human mechanism research of Minqin Oasis change in the lower reaches of the Shiyang river. Geographical Research, 2002, 21(4): 449-458.

[18]Brabec E, Schulte S, Richards P L. Impervious surfaces and water quality: A review of current literature and its implications for watershed planning. Journal of Planning Literature, 2002, 16(4): 499-514.

[19]Ridd M K. Exploring a V-I-S (vegetation-impervious surface-soil) model for urban ecosystem analysis through remote sensing: comparative anatomy for cities. International Journal of Remote Sensing, 1995, 16(12): 2165-2185.

[20]Lu D S, Tian H Q, Zhou G M, Ge H L. Regional mapping of human settlements in southeastern China with multisensory remotely sensed data. Remote Sensing of Environment, 2008, 112(9): 3668-3679.

[21]Pan T, Wu S H, He D M, Dai E F, Liu Y J. Ecological effects of longitudinal Range-Gorge land surface pattern and its regional differentiation. Acta Geographica Sinica, 2012, 67(1): 13-26.

[22]Tang G A, Yang X. Arcgis geographic information system spatial analysis experiment tutorial. Beijing: Science Press, 2006.

[23]Dong P L. Lacunarity analysis of raster datasets and 1D, 2D, and 3D point patterns. Computers & Geosciences, 2009, 35(10): 2100-2110.

參考文獻(xiàn):

[1]黎夏, 葉嘉安. 利用遙感監(jiān)測(cè)和分析珠江三角洲的城市擴(kuò)張過(guò)程: 以東莞市為例. 地理研究, 1997, 16(4): 56-62.

[6]車前進(jìn), 段學(xué)軍, 郭垚, 王磊, 曹有揮. 長(zhǎng)江三角洲地區(qū)城鎮(zhèn)空間擴(kuò)展特征及機(jī)制. 地理學(xué)報(bào), 2011, 66(4): 446-456.

[7]馬曉冬, 朱傳耿, 馬榮華, 蒲英霞. 蘇州地區(qū)城鎮(zhèn)擴(kuò)展的空間格局及其演化分析. 地理學(xué)報(bào), 2008, 63(4): 405-416.

[8]石培基, 王祖靜, 李巍. 石羊河流域地區(qū)城鎮(zhèn)空間擴(kuò)展格局演化. 地理科學(xué), 2012, 32(7): 840-845.

[11]李萬(wàn)凱. 應(yīng)用UrbanSim分析臺(tái)北市土地使用變遷之研究 [D]. 臺(tái)北: 臺(tái)灣大學(xué), 2008.

[12]劉珍環(huán), 王仰麟, 彭建, 謝苗苗, 李猷. 基于不透水表面指數(shù)的城市地表覆被格局特征——以深圳市為例. 地理學(xué)報(bào), 2011, 66(7): 961-971.

[13]岳文澤, 徐建華. 上海市人類活動(dòng)對(duì)熱環(huán)境的影響. 地理學(xué)報(bào), 2008, 63(3): 247-256.

[15]匡文慧, 劉紀(jì)遠(yuǎn), 陸燈盛. 京津唐城市群不透水地表增長(zhǎng)格局以及水環(huán)境效應(yīng). 地理學(xué)報(bào), 2011, 66(11): 1486-1496.

[17]楊永春, 李吉均, 陳發(fā)虎, Jacquie B, 李瑞成, 李丁, 常跟應(yīng), 李憶春. 石羊河下游民勤綠洲變化的人文機(jī)制研究. 地理研究, 2002, 21(4): 449-458.

[21]潘韜, 吳紹洪, 何大明, 戴爾阜, 劉玉潔. 縱向嶺谷區(qū)地表格局的生態(tài)效應(yīng)及其區(qū)域分異. 地理學(xué)報(bào), 2012, 67(1): 13-26.

[22]湯國(guó)安, 楊昕. ArcGIS地理信息系統(tǒng)空間分析實(shí)驗(yàn)教程. 北京: 科學(xué)出版社, 2006.