溫利華,王永芹,張廣錄,劉紅耀,王翠麗
(1.邯鄲學院地理與旅游系,河北 邯鄲 056005;2.中國科學院遺傳與發育生物學研究所農業資源研究中心,石家莊 050021;3.邢臺市農業科學研究院,河北 邢臺 054000)
生態環境變化已成為全球生態學界關注的焦點,而作為具有自然和人文雙重屬性的土地利用/覆蓋是人類活動與自然環境相互作用的集中反映[1]。土地利用/覆蓋變化可引起許多自然現象和生態過程的變化,如土壤養分和水分的變化、地表徑流與侵蝕、生物多樣性的分布和生物地球化學循環等[2-7]。因此,把握導致區域土地利用/覆蓋變化的時空過程,弄清驅動這一過程的關鍵因素及作用方式,有助于調整人類自身活動,對區域資源優化配置、保護和改善區域生態環境、促進區域可持續發展戰略順利實施具有重要的理論和應用意義。
海河流域地處我國中西部干旱半干旱地區,流域內包括京、津、冀、晉、豫、魯、蒙、遼8省(自治區、直轄市)。流域內有首都北京、直轄市天津,石家莊、唐山、秦皇島、廊坊、保定、邯鄲、邢臺等25座大中城市,該流域是我國經濟較發達地區以及重要的工業基地和高新技術產業基地,在國家經濟發展中具有重要戰略地位。而自從上世紀50年代末和80年代初的農村經濟體制改革以來,流域內的土地利用結構急劇變化,水資源強度開發與水資源短缺制約了該流域的可持續發展,水資源、水環境和生態問題突出,主要表現為河道斷流干涸、濕地萎縮、河口生態惡化、地下水嚴重超采、水土流失加劇、水污染加劇、生物多樣性減少等,土地利用變化的重要環境反應是水文行為的變化,而水文行為變化又會影響環境反應并反作用于土地利用,二者相互制約且密切關聯[8-12]。因此,土地利用/覆蓋變化的研究必將對該流域生態修復的水資源保障措施以及水資源優化配置、水污染防治及防災減災提供可借鑒的理論依據。本文以海河流域為研究對象,利用1970~2008年五期MSS和TM/ETM衛星遙感數據獲取該區域1970、1980、1990年初和2000年初的地表土地覆蓋特征。在GIS支持下,結合區域自然和社會經濟資料,分析闡述近30年來海河流域土地利用/覆蓋變化的時空過程,并初步探討驅動這一過程的關鍵因素。
海河流域東臨渤海,西倚太行,南界黃河,北接蒙古高原,流域總面積31.82萬km2,占全國總面積的3.3%,全流域總地勢是西北高東南低,大致分高原、山地及平原三種地貌類型。西部為山西高原和太行山區,北部為蒙古高原和燕山山區,面積18.94萬km2,占60%;東部和東南部為平原,面積12.84萬km2,占40%。該流域冬季受西伯利亞大陸性氣團控制,寒冷少雪;春季受蒙古大陸性氣團影響,風速大,氣候干燥,蒸發量大,往往形成干旱天氣;夏季受海洋性氣團影響,比較濕潤,氣溫高,降雨量多,且多暴雨,但因歷年夏季太平洋副熱帶高壓的進退時間、強度、影響范圍等很不一致,致使降雨量的變差大,旱澇時有發生;秋季為夏冬的過渡季節,一般年份秋高氣爽,降雨量較少。流域年平均氣溫在1.5~14℃,年平均相對濕度50%~70%;年平均降水量539 mm,屬半濕潤半干旱地帶;年平均陸面蒸發量470 mm,水面蒸發量1 100 mm。
采用1970年、1980年MSS影像,1990年TM影像,2000年、2008年ETM+影像,1996年以來的土地統計數據,海河流域8省市行政區圖及相關統計資料,其中遙感數據經幾何配準和數據重采樣后,分辨率統一為28.5 m×28.5 m,校準誤差小于0.5個像元。海河流域遙感影像鑲嵌采用用最小灰度差原則進行拼接,接縫處灰度加權平滑,用直方圖匹配方法進行灰度調整。
2008年9月對研究區土地利用類型進行初步調查,利用GARMIN GPSl2C手持式GPS定位,記錄、拍照的土地覆蓋類型與衛星影像圖比較分析。應用ERDAS圖像處理軟件,采用監督分類和目視解譯相結合方法完成遙感影像土地利用解譯,為了提高遙感影像中土地利用變化信息精度,對遙感影像進行拉伸增強、彩色變換和圖像融合增強處理,依據土地資源分類標準《土地利用現狀分類》,參照全國生態遙感監測土地利用/覆蓋分類體系并結合海河流域的土地覆蓋特征,確立該區的土地利用分類體系,土地利用類型一級分類為耕地、果園、林地、草地、水域和水利設施用地、工業用地、交通實施用地、住宅用地和其它土地共9類及20個二級類,參考同期的地面資料計算其分類錯誤矩陣和Kappa指數,對五期遙感影像的分類結果進行精度檢驗,Kappa指數均大于0.80,能夠滿足研究要求。將處理好的遙感影像圖轉入GIS軟件中,采用ARCGIS軟件的空間分析功能,得到海河流域1970~1990年和1990~2008年土地利用變化圖,提取土地利用變化信息,分析土地利用動態及空間格局變化。
土地利用動態度(表1、2)用來描述某種土地利用類型在研究區段內的變化劇烈程度,其數值的大小也在一定程度上反映了當地土地利用的格局變化特征。結果顯示:自1970年以來至本世紀初,流域內土地利用一直以耕地和林地為主導,但林地面積變化不大,耕地面積卻逐年減少,1970年耕地面積為135 550.15 km2,2008年減少至114 452.89 km2;但工業用地面積卻迅速增加,至2008年,工業用地面積增加為2 583.41 km2。另外果園的年變化率較大,工業用地、住宅用地和果園從1970~2008年的土地利用動態變化度呈迅速增長的趨勢,果園的土地利用動態度增幅達到44.31%,但耕地和其他用地的變化率呈明顯負增長,林地基本保持平衡狀態,水域及水利設施用地和交通運輸用地的年變化率比較小。由于種植結構的調整劇烈,經濟效益高的生產性用地類型表現出明顯的增長,如工業用地和果園。另一方面,人口的增長需要更多的住房來滿足人們的需要,住宅用地明顯增加。
全面的理解環境變化及其演化機制需要深入的研究產生這些變化的潛在時空過程,本文主要通過對五期土地利用與覆蓋類型圖進行兩兩疊加建立轉移概率矩陣(見表3、4)的方法,深一步探討海河流域各種土地利用與覆蓋類型相互轉換的內在機制。

表1 1970~2008年海河流域土地利用類型數量變化Table 1 Quantity change of types land use from 1970 to 2008 in Haihe river basin (km2)
1970年研究區內土地利用類型以耕地(水田、旱地)、林地(灌木林、疏林)和草地為主。1970至1990年土地利用類型轉化交錯復雜,變化較大的為耕地和果園,其次為工業用地。由于種植結構的調整,經濟效益高的果園等產業用地面積迅速增加,而耕地卻在這一期間轉化為果園、工業用地、住宅用地和其它土地利用類型。與上一階段相比,1990年相至2008年期間,流域內的果園、工業用地和住宅用地持續增加,而耕地卻持續減少,其他用地(包括濕地)同樣持續減少,這對海河流域的生態環境造成不良影響。這一階段正是我國大力提倡生態修復、生態建設的關鍵時期,在這一背景條件下,該區內植被得到一定的恢復,與此同時,隨著退耕還林政策的實施,大面積的耕地轉入林地,這使林地的面積也得到增加,就在國家加大生態環境保護力度的同時,也有破壞森林和毀林開荒的情況發生,表現為林地轉化為耕地和工業用地。
總體來看,從1970年至2008年期間各個土地利用類型的轉化比較復雜,每一種土地利用類型是向多個方向轉化,既有正向的也有負向的,耕地向其它土地利用類型的轉化比較明顯,林地、草地、工業用地和住宅用地由從耕地轉化而來。

表2 1970~2008年海河流域土地利用變化動態度Table 2 Dynamic degree of land use change from 1970 to 2008 in Haihe river basin (%)

表3 研究區1970~1990年土地利用轉移概率矩陣Table 3 Land use change transition probability matrix from 1970 to 1990 in the study area
土地利用動態變化驅動因素是指導致土地利用方式和目的發生變化的自然環境和社會經濟因素,自然環境因素主要包括氣候變動、植物演替、土壤類型變化過程、自然資源的有效性等;社會經濟驅動因素則主要包括人口變化、貧富狀況、經濟增長、政治經濟結構以及價值觀念等方面[13]。海河流域土地利用/覆蓋變化一方面與自然環境演變相關,另一方面與人類活動的不斷增強密切有關。

表4 研究區1990~2008年土地利用轉移概率矩陣Table 4 Land use change transition probability matrix from 1990 to 2008 in the study area
2.3.1 自然環境因素
研究表明,土壤屬性和氣候條件對土地利用結構具有制約作用[14]。海河流域處于我國西北部干旱半干旱地區,海河流域土壤劃分為內蒙古栗鈣土綿土區、華北山地棕壤褐土區、海河平原黃壚土潮土鹽土區,流域地形地貌獨特復雜,其地貌類型包括高山高原、低丘、中山、低山河谷和山間盆地,且流域內夏季氣候不穩定,降雨強度大,時間短促,容易引起水土流失,水土流失不僅會導致耕地面積數量減少和質量下降,而且對林地資源也會造成破壞。
2.3.2 社會經濟因素
土地利用與經濟發展之間,尤其是耕地與人口、GDP、工業增加值、服務業增加值、投資及城市化水平等社會經濟指標存在強相關關系[15]。土地利用與經濟發展之間存在著耕地減少,工業、商業等非農用地擴大的必然過程。海河流域GDP從1980年的1 592億元增加到2003年的25 750億元,增長了15倍。人均GDP從1980年的1 638元增長至2004年的1.92萬元,是全國平均水平的1.38倍,但流域內各地經濟發展不平衡,北京市天津市人均GDP達到中等發達國家水平,平原區經濟發達,GDP占全流域的82%,而山區GDP只占18%。流域內工業化、產業化的發展必然帶動城鎮的擴大。多數耕地資源都是以城鎮為中心向四周擴散,城鎮的急劇擴大必然要以犧牲耕地資源為代價,耕地面積必然顯著減少,工業用地、住宅用地、交通運輸用地及水利設施用地等都有不同程度的增加。
產業結構調整同樣會引起土地利用結構的變化,種植產業值和林業產業值穩定增長,糧食產量翻番,但傳統的耕地種植效益要遠遠小于特色經濟所帶來的經濟效益,農民為了增加自身收入,加大了園地的開墾,園地面積從1970年的110.13 km2增長到1 964.58 km2增幅近20倍,水田和旱地轉入經濟效益比較明顯的棉花、油料、麻類、煙葉等生產領域。
人口增加也對區域土地利用空間格局變化產生影響[16]。20世紀80年代以來,海河流域經濟社會持續發展,1980~2005年人口自然增長率雖呈降低趨勢,但流域總人口從9 721萬人增加至1.34億人,增長了38%。人口的增加,一方面會引起各種社會需求,如住房需求、就業需求,由此促進城鎮化和工業化的發展,工業用地、住宅用地和交通運輸用地明顯增加;另一方面人口增加也會帶來勞動力急劇增加從而促進經濟發展間接影響土地利用狀況。
a.從1970~2008年土地利用變化的分析,海河流域目前正處于快速工業化階段,該流域土地利用規模和結構發生很大變化,生產性用地用地(工業用地和果園)快速增長—耕地面積的迅速減少是研究區土地利用/覆蓋變化過程的突出表現,也正是該區域土地利用變化格局的顯著特點。另外該流域通過近年來的生態修復與治理,環境得到一定改善,林地面積變化不大,但耕地面積顯著減少,工業用地與住宅用地等其他用地明顯增加,而且這種變化也存在很大的空間差異性。
b.海河流域土地利用/覆蓋變化驅動力方面,主要受到自然和人為雙層因素的影響,自然因素又受限制于人為因素。自然生態條件制約研究區內土地利用類型格局和數量的分布,受海拔高度和坡度的影響,林地多分布丘陵區、低山區,其它類型則多分布在海拔較低、坡度較緩的平原區。另一方面,受經濟利益的驅動區內種植結構呈現出了耕地面積持續下降,果園面積持續上升的雙向動態變化,雖然這種動態變化能夠提高農民的收入帶來可觀的經濟效益,但是這種打破傳統的高強度的耕作模式勢必對該區域微環境及資源的可持續發展帶來一定的負面效應。
c.海河流域土地資源豐富,發展農業的潛力很大,應把農業和水利基礎設施建設放在重要位置,本世紀初水域及水利設施用地大幅度增加,而耕地面積卻從上世紀70年代至今逐年減少,目前對海河流域造成嚴重生態威脅的問題就是耕地資源非農化,特別是耕地轉化成工業用地與住宅用地問題尤其明顯,各級土地管理部門應采取切實有效的措施遏制流域內耕地不合理非農化,保證經濟效益社會效益和生態效益協調增長。
[1]李秀彬.全球變化的核心領域—土地利用/土地覆被變化的國際研究動向[J].地理學報,1996,51(6):553-557.
[2]傅伯杰,陳利頂,王軍.土地利用結構與生態過程[J].第四紀研究,2003,23(3):247-255.
[3]史培軍,袁藝,陳晉.深圳市土地利用變化對流域徑流的影響[J].生態學報,2001,21(7):1041-1049.
[4]陳利頂,張淑榮,傅伯杰.流域尺度土地利用與土壤類型空間分布的相關性研究[J].生態學報,2003,23(12):2497-2505.
[5]Dale V H.The relationship between land use change and climate change[J].Ecological Application,1997,7(3):753-769.
[6]傅伯杰,邱揚,王軍.黃土丘陵小流域土地利用變化對水土流失的影響[J].地理學報,2002,57(6):717-722.
[7]Turner B LⅡ.Two type of global environment change:definitional and spatial scale issues in their human dimensions[J].Global Environment Change,1990,1(1):14-22.
[8]陳欣,王兆騫,楊武德.紅壤小流域坡地不同利用方式對土壤磷素流失的影響[J].生態學報,2000,20(3):374-377.
[9]吳云,曾源,趙炎,等.基于MODIS數據的海河流域植被覆蓋度估算及動態變化分析[J].資源科學,2010,7(28):199-206.
[10]丁相毅,賈仰文,王浩,等.氣候變化對海河流域水資源的影響及其對策[J].自然資源學報,2010,4(9):78-87.
[11]陳利頂,傅伯杰,王軍.黃土丘陵區典型小流域土地利用變化研究:以陜西延安地區大南溝流域為例[J].地理科學,2001,21(1):46-51.
[12]Ziaja W.Mountain landscape structure and dynamics under the20(Th)centuryclimate warmingin centralNordenskioldland,Spitsbergen[J].Ekologia-Bratislava,2004,23:374-383.
[13]史培軍.土地利用/覆蓋變化研究的方法與實踐[M].北京:科學出版社,2000:30-60.
[14]姚偉,樓建悅,張麗蘇.1985-2004年間錢塘江流域土地利用/覆蓋變化[J].科技通訊,2007,4(13):78-87.
[15]王業僑.海南省經濟社會發展與土地利用相關分析[J].地域研究與開發,2006,3(18):374-377.
[16]劉慶,陳利根,張鳳榮.中國建設占用耕地數量與人口增長關系實證[J].中國人口:資源與環境,2009,5(22):115-120.