基于土地利用結(jié)構(gòu)變化的深圳市水土流失動(dòng)態(tài)研究

徐 翼，李天宏

（1.北京大學(xué) 深圳研究生院 城市人居環(huán)境科學(xué)與技術(shù)重點(diǎn)實(shí)驗(yàn)室，深圳 518055；2.北京大學(xué)環(huán)境工程系水沙科學(xué)教育部重點(diǎn)實(shí)驗(yàn)室，北京 100871）

基于土地利用結(jié)構(gòu)變化的深圳市水土流失動(dòng)態(tài)研究

徐 翼1，2，李天宏1，2

（1.北京大學(xué) 深圳研究生院 城市人居環(huán)境科學(xué)與技術(shù)重點(diǎn)實(shí)驗(yàn)室，深圳 518055；2.北京大學(xué)環(huán)境工程系水沙科學(xué)教育部重點(diǎn)實(shí)驗(yàn)室，北京 100871）

采用區(qū)域水土流失快速評(píng)估方法，通過參照點(diǎn)選擇、參照組構(gòu)建、參照組匹配、O／E計(jì)算等步驟，利用深圳市1996年、2000年、2002年、2004年和2006年的土地利用類型數(shù)據(jù)，以深圳市區(qū)為基本單元，進(jìn)行多時(shí)段的水土流失動(dòng)態(tài)評(píng)估。結(jié)果表明：1996年，福田區(qū)土壤侵蝕程度明顯嚴(yán)重于其它5個(gè)區(qū)；2000年，寶安區(qū)、福田區(qū)和羅湖區(qū)的土壤侵蝕較嚴(yán)重；2002年、2004年和2006年，寶安區(qū)和福田區(qū)土壤侵蝕迅速加劇，南山區(qū)土壤侵蝕繼續(xù)增加，龍崗區(qū)稍有減緩的趨勢，羅湖區(qū)和鹽田區(qū)土壤侵蝕有所改善。評(píng)估得到的土壤侵蝕動(dòng)態(tài)變化特征與研究區(qū)域的開發(fā)過程也較為吻合。

土壤侵蝕；土地利用結(jié)構(gòu)；快速評(píng)估方法；深圳市

水土流失是全球性的環(huán)境問題［1］，也是全球變化研究（如IGBP，IHDP）的核心內(nèi)容［2］。水土流失的評(píng)估可以在3個(gè)空間尺度（坡面、小流域、區(qū)域）上進(jìn)行［3］。區(qū)域水土流失研究是在坡面和小流域研究的基礎(chǔ)上，重點(diǎn)探討水土流失的宏觀發(fā)生機(jī)制和影響因素等，強(qiáng)調(diào)人為因素對(duì)水土流失的正反作用，其研究結(jié)果可直接服務(wù)于區(qū)域環(huán)境規(guī)劃和治理。目前對(duì)區(qū)域水土流失進(jìn)行評(píng)估的方法仍沿用USLE或者 RUSLE［4～6］、水文法和水保法［7］。由于這些模型需要的基本資料除了個(gè)別研究區(qū)域以外，往往不能得到滿足，同時(shí)人類活動(dòng)已經(jīng)深刻地影響到水土流失過程，成為評(píng)估水土流失不可忽視的因素，所以都給區(qū)域性的水土流失評(píng)估帶來困難；另一方面，區(qū)域性的土地利用是人類利用土地的結(jié)果，其數(shù)據(jù)也相對(duì)容易獲得。鑒于此，有研究者在分析土地利用和水土流失影響因素的基礎(chǔ)上，提出了區(qū)域土壤侵蝕快速評(píng)估方法，即通過土地利用結(jié)構(gòu)的變化來進(jìn)行水土流失動(dòng)態(tài)評(píng)估［8，9］。深圳市是我國最早設(shè)立的經(jīng)濟(jì)特區(qū)，30多年來經(jīng)濟(jì)的快速發(fā)展，加速了水土流失等環(huán)境問題的惡化。同黃土高原等水土流失嚴(yán)重的地區(qū)相比，深圳市水土流失的一個(gè)顯著特點(diǎn)就是人類活動(dòng)的影響尤其突出，通過土地利用結(jié)構(gòu)變化來監(jiān)測水土流失的動(dòng)態(tài)變化是一個(gè)合適而快捷的方法。本文將利用深圳市1996年、2000年、2002年、2004年和2006年的土地利用數(shù)據(jù)，采用區(qū)域土壤侵蝕快速評(píng)估方法進(jìn)行深圳市土壤侵蝕的動(dòng)態(tài)評(píng)估，考察土壤侵蝕的變化趨勢。

1 研究區(qū)域

深圳市地處廣東省南部，毗鄰香港。位于北回歸線以南，東經(jīng)113°46′至114°37′，北緯22°27′至22°52′。地處南亞熱帶氣候帶，多年平均降雨量1883mm，土地總面積為1 960 km2，其中，深圳經(jīng)濟(jì)特區(qū)面積為395.81 km2。全市的植被覆蓋率為57.68%，主要集中在大鵬半島和梧桐山一帶，原生植被類型為亞熱帶常綠闊葉林，人工植被主要有大葉相思樹、臺(tái)灣相思樹、桉樹等，結(jié)構(gòu)簡單，林相單一，蓋度低，生態(tài)穩(wěn)定性較差。全市共有大小河流160多條，其中流域面積大于100 km2的有深圳河、茅洲河、觀瀾河、龍崗河和坪山河。由于深圳成立經(jīng)濟(jì)特區(qū)，又有毗鄰香港的區(qū)位優(yōu)勢，成立特區(qū)以來，經(jīng)濟(jì)發(fā)展迅速，市郊的農(nóng)村城市化進(jìn)程很快，而在快速發(fā)展的城市化過程中，開發(fā)活動(dòng)的外力沖擊遠(yuǎn)遠(yuǎn)超過了土體的自我保持能力，再加上高強(qiáng)度暴雨作用和大規(guī)模的人為擾動(dòng)作用，大大加速了土壤侵蝕的過程，造成了嚴(yán)重的土壤侵蝕，對(duì)深圳市的經(jīng)濟(jì)社會(huì)可持續(xù)發(fā)展也構(gòu)成嚴(yán)重威脅。開發(fā)建設(shè)過程中土壤侵蝕的變化及如何采取相應(yīng)的水土保持措施，已成為深圳市持續(xù)關(guān)注的焦點(diǎn)。

2 研究方法

本文采用的區(qū)域土壤侵蝕快速評(píng)估方法［8，9］基于這樣一個(gè)假設(shè)，即降雨、地形和土壤等環(huán)境變量與土地利用結(jié)構(gòu)變量能夠共同決定區(qū)域土壤侵蝕變化的趨勢。這種方法不側(cè)重于獲得土壤侵蝕量的精確結(jié)果，而是致力于獲取相對(duì)于某一參考狀態(tài)的土壤侵蝕動(dòng)態(tài)變化趨勢。

2.1 土地利用結(jié)構(gòu)特征指標(biāo)

土地利用結(jié)構(gòu)變化是影響區(qū)域土壤侵蝕的活躍因素，是合理評(píng)估區(qū)域土壤侵蝕尤其是動(dòng)態(tài)評(píng)估和變化趨勢預(yù)測的關(guān)鍵。根據(jù)區(qū)域內(nèi)各土地利用類型來分析其對(duì)區(qū)域土壤侵蝕的綜合影響，需要協(xié)調(diào)不同的土地利用類型對(duì)于區(qū)域土壤侵蝕的貢獻(xiàn)。這可以通過引入一個(gè)反映土地利用結(jié)構(gòu)對(duì)區(qū)域土壤侵蝕影響的綜合指標(biāo)——土地利用結(jié)構(gòu)特征指標(biāo)來實(shí)現(xiàn)，不同土地利用類型對(duì)土壤侵蝕影響的相對(duì)重要性可通過用層次分析法確定的權(quán)重來反映。區(qū)域土地利用結(jié)構(gòu)特征指標(biāo)定義為

式中：SI為區(qū)域土地利用結(jié)構(gòu)特征指標(biāo)；wi為不同土地利用類型對(duì)區(qū)域水土保持的重要性權(quán)重；Ai為對(duì)應(yīng)土地利用類型的面積百分比。

2.2 評(píng)估步驟

2.2.1 參照組構(gòu)建

選取有實(shí)測資料、受人類影響活動(dòng)較少的若干區(qū)域作為參照點(diǎn)，根據(jù)各區(qū)域?qū)?yīng)的土地利用結(jié)構(gòu)（針對(duì)水土流失目標(biāo)的）特點(diǎn)按其均勻性分組，形成不同的參照組，組中的各區(qū)域并不一定包含在同一地區(qū)或鄰近地區(qū)，這部分工作可通過聚類分析方法實(shí)現(xiàn)。然后計(jì)算確定不同參照組土地利用結(jié)構(gòu)的平均狀況。

針對(duì)影響區(qū)域水土流失的各種因素，利用多變量判別分析和相關(guān)分析等統(tǒng)計(jì)方法，選擇一組受人類活動(dòng)影響較小的環(huán)境變量子集（如海拔高度、高差、平均坡度、地表物質(zhì)組成等，對(duì)應(yīng)于不同空間尺度的環(huán)境變量可能有所不同），使其能夠在一定程度上解釋參照組間的大部分差異。

2.2.3 待評(píng)估區(qū)域的參照組匹配

確定待評(píng)估區(qū)域?qū)?yīng)的環(huán)境變量，根據(jù)其環(huán)境變量信息，利用判別分析等方法判別該區(qū)域所屬的參照組，同時(shí)計(jì)算該區(qū)域?qū)儆诿恳粎⒄战M的概率。

2.2.4 待評(píng)估區(qū)域的土地利用結(jié)構(gòu)預(yù)測

根據(jù)待評(píng)估區(qū)域?qū)儆诿恳粎⒄战M的概率以及不同參照組的土地利用結(jié)構(gòu)信息，可以計(jì)算得到該評(píng)估區(qū)域預(yù)期的土地利用結(jié)構(gòu)，其表達(dá)式為

式中：ˉri為待評(píng)估區(qū)域的預(yù)期土地利用結(jié)構(gòu)；i和j分別為表現(xiàn)因子和參照組下標(biāo)；rij為對(duì)應(yīng)參照組類中心的土地利用結(jié)構(gòu)；k為參照組總數(shù)；wij為待評(píng)估區(qū)域?qū)儆趯?duì)應(yīng)參照組的概率。

2.2.5 預(yù)測結(jié)果與實(shí)際觀測結(jié)果的比較分析

評(píng)估區(qū)域土地利用結(jié)構(gòu)特征的實(shí)際觀測結(jié)果與預(yù)期之比值用O／E定義。O／E越接近于1，說明實(shí)際觀測的與預(yù)期的土地利用結(jié)構(gòu)特征越符合。

據(jù)了解，按照《2018年河南省家政服務(wù)領(lǐng)域市場秩序?qū)ｍ?xiàng)整治方案》《家庭服務(wù)管理暫行辦法》等相關(guān)規(guī)定，焦作市將重點(diǎn)對(duì)家政服務(wù)機(jī)構(gòu)是否按要求公開服務(wù)項(xiàng)目、收費(fèi)標(biāo)準(zhǔn)和設(shè)立投訴監(jiān)督電話；是否建立工作檔案、跟蹤管理制度，對(duì)消費(fèi)者和家政服務(wù)員之間的投訴是否妥善處理；是否制定家政服務(wù)合同；是否及時(shí)報(bào)送經(jīng)營情況信息等方面進(jìn)行整治。

3 深圳市土壤侵蝕動(dòng)態(tài)評(píng)估

3.1 單元?jiǎng)澐?/h3>

本文采用的是基本數(shù)據(jù)包括深圳市行政邊界圖，深圳市1996，2000，2002，2004和2006年土地利用矢量數(shù)據(jù)，土地利用類型主要分為未開發(fā)用地、園地、建成區(qū)、耕地、林地和水域等6類。數(shù)據(jù)的格式是ARCINFO的coverage數(shù)據(jù)，所有操作在ARCGIS 9.0中進(jìn)行。

通過將行政邊界圖和幾期的土地利用類型圖分別進(jìn)行疊加操作，獲得計(jì)算單元，包括深圳市的4個(gè)城區(qū)和2個(gè)近郊區(qū)，其基本情況如表1。

表1 研究區(qū)域單元?jiǎng)澐忠挥[表Table 1 The partition of researched region

3.2 土地利用結(jié)構(gòu)特征指標(biāo)確定

利用深圳市6個(gè)時(shí)相土地利用圖，在GIS中可以計(jì)算出各個(gè)行政區(qū)域的6種土地利用類型的面積比例。利用層次分析法，根據(jù)深圳市特區(qū)建設(shè)的特點(diǎn)，由專家打分可以獲得每種土地利用類型對(duì)水土流失的影響權(quán)重，本文采用文獻(xiàn)［10］在深圳市寶安區(qū)進(jìn)行類似研究的權(quán)重值（表2）。

根據(jù)式（1），可以計(jì)算深圳市各個(gè)行政單元在不同年限的綜合土地利用結(jié)構(gòu)特征指標(biāo)SI。

表2 土地利用／覆蓋類型水土保持效應(yīng)判斷矩陣Table 2 Pair-wise comparison matrix and weights for land use c／over sort subjected to soil erosion

3.3 環(huán)境因子確定

影響土壤侵蝕的環(huán)境因子主要包括地形、地表物質(zhì)組成、降雨等因子。研究區(qū)域的地表物質(zhì)組成比較單一，在各個(gè)區(qū)分布也比較均勻，由于研究的區(qū)域范圍較小，整個(gè)深圳市的降雨分布基本是均勻分布，所以均不能反映不同區(qū)域間的差異；植被狀況由于主要反映在土地利用類型中，為了保證表現(xiàn)因子與環(huán)境因子的相對(duì)獨(dú)立性，也不作為環(huán)境因子考慮。經(jīng)過綜合考慮，最終選擇區(qū)域的平均坡度、平均高度、高差和起伏度4個(gè)指標(biāo)作為影響深圳市土壤侵蝕的環(huán)境因子。

3.4 動(dòng)態(tài)評(píng)估

選擇深圳市1996年的6個(gè)行政區(qū)域的土地利用結(jié)構(gòu)作為參照點(diǎn)，依次評(píng)估2000年、2002年、2004年和2006年土壤侵蝕相對(duì)于1996年的變化情況。然后運(yùn)用SPSS，對(duì)土地利用結(jié)構(gòu)特征指標(biāo)進(jìn)行聚類分析（如圖1），創(chuàng)建參照組。

圖1 深圳市各參照點(diǎn)土地利用結(jié)構(gòu)特征指標(biāo)聚類分析Fig.1 The result of hierarchical feature cluster analysis of land-use structure for every district of Shenzhen City

從聚類分析圖1可以看出，研究區(qū)域可以分為4個(gè)參照組，第一參照組Gp1包括2個(gè)參照點(diǎn)即羅湖區(qū)和龍崗區(qū)，第二參照組Gp2包括1個(gè)參照點(diǎn)即寶安區(qū)，第三參照組Gp3包括1個(gè)參照點(diǎn)即鹽田區(qū)，第四參照組Gp4包括2個(gè)參照點(diǎn)即福田區(qū)和南山區(qū)。確定參照組后，根據(jù)不同的參照組中參照點(diǎn)的土地利用結(jié)構(gòu)，采用算數(shù)平均的方法計(jì)算各個(gè)參照組類中心的土地利用結(jié)構(gòu)平均特征，參照組的匹配過程可以通過SPSS中的判別分析功能實(shí)現(xiàn)。通過待評(píng)估區(qū)域的環(huán)境變量與不同參照組的環(huán)境變量特征進(jìn)行匹配，可以確定待評(píng)估區(qū)域分別屬于各參照組的概率，深圳市利用環(huán)境變量判別得到的參照組歸屬準(zhǔn)確率為100%，說明本研究區(qū)域利用環(huán)境變量可以較好地解釋不同參照組間土地利用結(jié)構(gòu)特征指標(biāo)的大部分差異，符合快速評(píng)估方法對(duì)環(huán)境變量和參照組的要求。同時(shí)根據(jù)各參照組的土地利用結(jié)構(gòu)平均特征和歸屬概率，可以預(yù)測不同子流域的土地利用結(jié)構(gòu)，從而得到對(duì)應(yīng)的預(yù)期土地利用結(jié)構(gòu)特征指標(biāo)。

表3 深圳市各區(qū)土地利用結(jié)構(gòu)的預(yù)測值及土地利用結(jié)構(gòu)特征指標(biāo)Table 3 The estimated values and indexes of land-use structure to different districts in Shenzhen City

計(jì)算不同的行政單元在不同的時(shí)段的土地利用結(jié)構(gòu)特征指標(biāo)的實(shí)測值與預(yù)期值的比值，得到土地利用結(jié)構(gòu)特征指標(biāo)O／E比，如表4所示。

表4 土地利用結(jié)構(gòu)特征指標(biāo)比（O／E）Table 4 The ratio of land-use structure index

4 結(jié)果和討論

研究表明（如圖2），深圳市各區(qū)從1996年到2006年間土壤侵蝕程度有著不同的變化。主要分為幾個(gè)階段：①1996年到2000年，福田區(qū)、龍崗區(qū)和鹽田區(qū)土壤侵蝕變化相對(duì)緩慢并有改善的趨勢，而寶安區(qū)、羅湖區(qū)和南山區(qū)土壤侵蝕有所加劇（O／E相對(duì)降低），主要與此階段這3個(gè)區(qū)的開發(fā)力度加大是分不開的；②2000年到2002年，福田區(qū)和龍崗區(qū)因其建成區(qū)規(guī)模的不斷加大，土壤侵蝕也迅速加?。∣／E迅速降低），而南山區(qū)、羅湖區(qū)和寶安區(qū)土壤侵蝕受到相對(duì)的關(guān)注，所以有所改善（O／E增加）主要與區(qū)內(nèi)林地面積的增加有關(guān)，但跟1996年相比其土壤侵蝕依然很嚴(yán)重，鹽田區(qū)由于林地所占的比重較其它各區(qū)都大，其土壤侵蝕的程度就會(huì)相對(duì)的低一些；③2002年到2004年，隨著特區(qū)對(duì)土壤侵蝕治理力度的加強(qiáng)，福田、龍崗和南山的土壤侵蝕并沒有劇烈加劇，羅湖區(qū)和鹽田區(qū)更是有了相應(yīng)的改善，但是寶安區(qū)由于對(duì)林地的不保護(hù)以及開發(fā)建設(shè)導(dǎo)致土壤侵蝕有所加劇，這是不容忽視的；④2004年到2006年，寶安區(qū)注意到了土壤侵蝕的危害及嚴(yán)重性，并沒有繼續(xù)使其加劇，但其他各區(qū)土壤侵蝕相應(yīng)的都有所加?。∣／E降低）。從總體趨勢看，深圳市2006年土壤侵蝕相對(duì)于1996年呈加劇趨勢，應(yīng)該引起相當(dāng)?shù)年P(guān)注。

從另一個(gè)角度看，同一年份各區(qū)的土壤侵蝕程度也有所不同：①1996年，福田區(qū)土壤侵蝕程度明顯嚴(yán)重于另5個(gè)區(qū)；②2000年，土壤侵蝕程度依次為寶安區(qū)＞福田區(qū)＞羅湖區(qū)＞南山區(qū)＞龍崗區(qū)＞鹽田區(qū)；③2002年，則是福田區(qū)＞寶安區(qū)＞龍崗區(qū)＞羅湖區(qū)＞鹽田區(qū)＞南山區(qū)；④2004年，寶安區(qū)＞福田區(qū)＞龍崗區(qū)＞南山區(qū)＞羅湖區(qū)＞鹽田區(qū)；⑤2006年，寶安區(qū)＞福田區(qū)＞南山區(qū)＞龍崗區(qū)＞羅湖區(qū)＞鹽田區(qū)?？梢钥闯?，寶安區(qū)和福田區(qū)也一直是深圳市土壤侵蝕最嚴(yán)重的區(qū)，其土壤侵蝕加劇速度明顯快于其他4個(gè)區(qū)。

圖2 O／E綜合變化趨勢Fig.2 The total variational trend of O／E

從深圳市各年份土壤侵蝕空間變化趨勢分析也可以看出以上的變化趨勢（如圖3）。對(duì)于寶安區(qū)和龍崗區(qū)各鎮(zhèn)：1996年，土壤侵蝕程度平湖鎮(zhèn)＞公明鎮(zhèn)＞坑梓鎮(zhèn)＞新安街道＞龍華鎮(zhèn)＞觀瀾鎮(zhèn)＞布吉鎮(zhèn)＞龍崗鎮(zhèn)＞坪地鎮(zhèn)＞光明街道＞石巖鎮(zhèn)＞橫崗鎮(zhèn)＞坪山鎮(zhèn)＞西鄉(xiāng)鎮(zhèn)＞松崗鎮(zhèn)＞福永鎮(zhèn)＞大鵬鎮(zhèn)＞南澳鎮(zhèn)＞葵涌鎮(zhèn)＞沙井鎮(zhèn)；2000年，平湖鎮(zhèn)＞布吉鎮(zhèn)＞觀瀾鎮(zhèn)＞公明鎮(zhèn)＞龍華鎮(zhèn)＞坑梓鎮(zhèn)＞龍崗鎮(zhèn)＞西鄉(xiāng)鎮(zhèn)＞光明街道＞坪地鎮(zhèn)＞松崗鎮(zhèn)＞橫崗鎮(zhèn)＞新安街道＞石巖鎮(zhèn)＞福永鎮(zhèn)＞坪山鎮(zhèn)＞沙井鎮(zhèn)＞大鵬鎮(zhèn)＞葵涌鎮(zhèn)＞南澳鎮(zhèn)；2002年，平湖鎮(zhèn)＞光明街道＞坑梓鎮(zhèn)＞觀瀾鎮(zhèn)＞布吉鎮(zhèn)＞松崗鎮(zhèn)＞龍華鎮(zhèn)＞公明鎮(zhèn)＞石巖鎮(zhèn)＞龍崗鎮(zhèn)＞坪地鎮(zhèn)＞西鄉(xiāng)鎮(zhèn)＞新安街道＞橫崗鎮(zhèn)＞沙井鎮(zhèn)＞坪山鎮(zhèn)＞福永鎮(zhèn)＞大鵬鎮(zhèn)＞南澳鎮(zhèn)＞葵涌鎮(zhèn)；2004年，平湖鎮(zhèn)＞松崗鎮(zhèn)＞龍華鎮(zhèn)＞布吉鎮(zhèn)＞光明街道＞觀瀾鎮(zhèn)＞沙井鎮(zhèn)＞坑梓鎮(zhèn)＞公明鎮(zhèn)＞龍崗鎮(zhèn)＞石巖鎮(zhèn)＞橫崗鎮(zhèn)＞坪地鎮(zhèn)＞新安街道＞西鄉(xiāng)鎮(zhèn)＞福永鎮(zhèn)＞坪山鎮(zhèn)＞大鵬鎮(zhèn)＞南澳鎮(zhèn)＞葵涌鎮(zhèn)；2006年，新安街道＞平湖鎮(zhèn)＞公明鎮(zhèn)＞沙井鎮(zhèn)＞觀瀾鎮(zhèn)＞松崗鎮(zhèn)＞布吉鎮(zhèn)＞坑梓鎮(zhèn)＞龍華鎮(zhèn)＞石巖鎮(zhèn)＞龍崗鎮(zhèn)＞坪地鎮(zhèn)＞橫崗鎮(zhèn)＞光明街道＞福永鎮(zhèn)＞西鄉(xiāng)鎮(zhèn)＞坪山鎮(zhèn)＞大鵬鎮(zhèn)＞葵涌鎮(zhèn)＞南澳鎮(zhèn)。因此，總體來說，寶安區(qū)土壤侵蝕比龍崗區(qū)嚴(yán)重；大棚鎮(zhèn)、葵涌鎮(zhèn)和南澳鎮(zhèn)的水土保持是最好的，土壤侵蝕程度明顯低于其他各鎮(zhèn)。

圖3 深圳市各年份土壤侵蝕變化趨勢Fig.3 The change of soil erosion in Shenzhen City

5 結(jié)論

快速評(píng)估的方法可以有效地評(píng)估區(qū)域的土壤侵蝕動(dòng)態(tài)變化趨勢，在其它數(shù)據(jù)缺少的情況下，利用土地利用結(jié)構(gòu)變化數(shù)據(jù)可以快速方便地獲得水土流失的動(dòng)態(tài)變化。利用快速評(píng)估方法對(duì)深圳市水土流失的動(dòng)態(tài)變化分析，表明：從時(shí)間序列上來看，從1996年到2006年，福田區(qū)、南山區(qū)和寶安區(qū)的土壤侵蝕較為嚴(yán)重，處于不斷的加劇中，且寶安區(qū)的變化最為劇烈，這也與該區(qū)的開發(fā)建設(shè)分不開的；羅湖區(qū)和龍崗區(qū)的土壤侵蝕較輕且變化較緩；鹽田區(qū)在2002年之后土壤侵蝕更是有了改善的趨勢。從空間分布角度來看，同一年份土壤侵蝕程度基本是寶安區(qū)和福田區(qū)最強(qiáng)，龍崗區(qū)、南山區(qū)和羅湖區(qū)次之，鹽田區(qū)最弱；寶安區(qū)又以觀瀾鎮(zhèn)、松崗鎮(zhèn)等的土壤侵蝕最強(qiáng)，龍崗區(qū)則是平湖鎮(zhèn)最強(qiáng)。當(dāng)然在開發(fā)過程中，某些鎮(zhèn)生態(tài)環(huán)境遭到了嚴(yán)重的破壞，如沙井鎮(zhèn)和福永鎮(zhèn)，土壤侵蝕變化相當(dāng)劇烈，而另一些鎮(zhèn)，如南澳鎮(zhèn)、大鵬鎮(zhèn)和葵涌鎮(zhèn)的水土保持工作很好，土壤侵蝕也有所改善?？傮w來說，羅湖區(qū)和鹽田區(qū)土壤侵蝕較其它各區(qū)要小，應(yīng)予以保持；而對(duì)于土壤侵蝕相對(duì)比較嚴(yán)重的福田區(qū)和寶安區(qū)各鎮(zhèn)，應(yīng)予以更多的關(guān)注，在開發(fā)建設(shè)和經(jīng)濟(jì)發(fā)展的同時(shí)，要關(guān)注其生態(tài)環(huán)境的變化。

［1］ MITRA B，SCOTT H D，DIXON J C，et al.Application of Fuzzy Logic to the Prediction of Soil Erosion in a Large WaterShed［J］.Geoderma，1998，86：183-209.

［2］ JETTEN V，DEROO A PJ，F(xiàn)AVIS-MORTBCK D.E-valuation of Field-scale and Catchment-scale Soil Erosion Models［J］.Catena，1999，37（3-4）：521-541.

［3］ 楊勤科，李 銳.中國水土流失和水土保持定量研究進(jìn)展［J］.水土保持通報(bào)，1998，18（5）：13-18.

［4］ 楊艷生.區(qū)域性土壤流失預(yù)測方程的初步研究［J］.土壤學(xué)報(bào)，1999，27（1）：73-78.

［5］ 胡良軍，李 銳，楊勤科，等.基于GIS的區(qū)域水土流失評(píng)價(jià)模型［J］.應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報(bào)，2000，8（1）：1-8.

［6］ NYAKATAWA E Z，REDDY K C，LEMUNYON J L.Predicting Soil Erosion in Conservation Tillage Cotton Production Systems Using the Revised Universal Soil Loss Equation（RUSLE）［J］.Soil&Tillage Research，2001，57：213-224.

［7］ 張勝利，李 倬，趙文琳，等.黃河中游多沙粗沙區(qū)水沙變化原因及發(fā)展趨勢［M］.鄭州：黃河水利出版社，1998.

［8］ 倪晉仁，李英奎.基于土地利用結(jié)構(gòu)變化的水土流失動(dòng)態(tài)評(píng)估［J］.地理學(xué)報(bào)，2001，56（5）：611-621.

［9］ NI JR，LI Y K.Approach to Soil Erosion Assessment in Terms of Land-use Structure Changes［J］.Soil Water Conserv.2003，58（3）：158-169.

［10］李天宏，張 平.遙感圖象軟分類和AHP在區(qū)域水土流失評(píng)估中的應(yīng)用［J］.應(yīng)用基礎(chǔ)與工程科學(xué)學(xué)報(bào)，2002，10（2）：126-135.

［11］楊德生，肖衛(wèi)國，陳冬奕，等.深圳市水土流失動(dòng)態(tài)變化遙感監(jiān)測分析［J］.人民珠江，2002，（6）：55-57.

Dynamic Evaluation of Soil Erosion in Shenzhen Based on Land-Use Structural Variation

XU Yi1，2，LI Tian-hong1，2

（1.Shenzhen Graduate School of Peking University，Shenzhen 518055，China；2.College of Environmental Sciences，The Key Laboratory for Water and Sediment Sciences，Peking University Beijing 100871，China）

Using the rapid assessment method for soil loss and some kinds of data of land-use types in Shenzhen City in 1996，2000，2002，2004 and 2006，the changes of soil erosion in these periods were evaluated with a district as evaluation unit.The results showed as follows：In 1996，the soil erosion degree in Futian District was more serious than that in the other five Districts.In 2000，the soil erosion degree in Baoan District，F(xiàn)utian District and Luohu District was more serious.In 2002，2004 and 2006，the soil erosion in Baoan District and Futian District was aggravated more quickly than that in Nanshan District，but the soil erosion in Longgang District was lower than that before，and the status of soil erosion in Luohu District and Yantian District was even improved.These results are in good agreement with the development process in the study area.

soil erosion；land-use structure；rapid assessment method；Shenzhen City

S157

A

1001-5485（2009）07-0006-05

2008-08-28

973國家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃項(xiàng)目（2007CB40720）

徐 翼（1985-），女，貴州金沙人，理學(xué)碩士，主要從事工程環(huán)境模擬與信息系統(tǒng)研究，（電話）010-62753351，13601324905（電子信箱）litianhong pku＠vip.sina.com。

（編輯：易興華）