張慧榮, 畢如田, 荊耀棟, 朱洪芬, 劉慧芳
(山西農(nóng)業(yè)大學資源環(huán)境學院,山西晉中 030801)
加強對生態(tài)環(huán)境的保護逐漸成為新時代生態(tài)文明建設(shè)的重要議題,劃分生態(tài)功能區(qū)有助于有效識別不同土地利用類型的生態(tài)功能特點,使土地資源被充分合理利用,進而可確定生態(tài)環(huán)境保護與生態(tài)建設(shè)的重點區(qū)域,有效合理地踐行可持續(xù)發(fā)展戰(zhàn)略[1]. 自然地理要素與土地利用的空間配置具有時空異質(zhì)性且受制于區(qū)域地形地貌特征,在地域分異理論[2]的指導下,有部分學者將地貌特征、典型生態(tài)問題、主導生態(tài)服務(wù)功能與敏感特征相似的單元進行組合,以劃分合理的生態(tài)功能區(qū),更好地服務(wù)于區(qū)域經(jīng)濟發(fā)展與生態(tài)環(huán)境治理. 目前,關(guān)于生態(tài)功能分區(qū)的研究涉及生態(tài)區(qū)劃的原則和依據(jù)、不同類型的功能分區(qū)[3-4]、不同生物資源[5]、不同尺度生態(tài)區(qū)景觀格局[6]等方面,且大多研究考慮自然因素居多,缺少對人類活動干擾要素的考慮. 已報道的關(guān)于生態(tài)功能分區(qū)的對象涵蓋了流域[7]、國家公園[8]、湖泊[9]、民族地區(qū)[10]、城市[11]、濕地[12]等,用于劃分生態(tài)功能區(qū)的方法主要有空間疊置法[13]、指標權(quán)重分區(qū)法[14]、聚類法[15]、主成分分析法[16]等,而且還有學者開始利用RS、GIS等空間信息技術(shù)與神經(jīng)網(wǎng)絡(luò)、模型軟件等非空間信息技術(shù)相結(jié)合的方法來劃分生態(tài)功能區(qū)[17-18]. 目前的研究成果為生態(tài)功能區(qū)的劃分提供了較好的理論基礎(chǔ)和技術(shù)支持,具有很好的借鑒性,但大多研究側(cè)重于生態(tài)服務(wù)功能層面的考量,忽略了基底要素與干擾要素對生態(tài)系統(tǒng)的影響,造成分區(qū)特性的單一化,而且依據(jù)生態(tài)敏感性分布特征及生態(tài)系統(tǒng)服務(wù)功能重要性分布規(guī)律進行生態(tài)功能分區(qū)的研究較少. 根據(jù)地貌空間分異特征以及人類活動干擾下引發(fā)的生態(tài)環(huán)境問題對區(qū)域生態(tài)系統(tǒng)服務(wù)功能重要性和生態(tài)敏感性的空間分布特征進行分析,并合理劃分區(qū)域生態(tài)功能區(qū),是維護生態(tài)系統(tǒng)穩(wěn)定發(fā)展的有效手段.
山西省晉中市位于太行山脈中段,地處黃土高原東部邊緣、我國地勢第一級階梯向第二級階梯過渡地帶,境內(nèi)地貌類型復雜,生態(tài)環(huán)境空間異質(zhì)性特征明顯. 鑒于此,以晉中市為研究區(qū),從“地貌-重要性-敏感性”的分析框架著手,利用綜合評價法及空間分析手段對研究區(qū)生態(tài)系統(tǒng)服務(wù)功能重要性和生態(tài)敏感性進行評價,并對其生態(tài)功能區(qū)進行劃分. 本研究對促進山水林田湖草生態(tài)系統(tǒng)的良性循環(huán)具有重要的現(xiàn)實意義.
晉中市(111°25′E~114°08′E,36°40′N~38°05′N)位于山西省中部,東依太行,西傍汾河,全市總面積16 400 km2,約占山西省總面積的10.5%,共轄11個縣(市、區(qū))(圖1). 市域氣候溫和,屬暖溫帶大陸性季風氣候,年均氣溫約為9 ℃,年均降水量約為480 mm,主要集中在夏秋之交. 全市河流以太行山脈、太岳山中脊為界,分屬黃河流域和海河流域. 太行山脈是全國重要的水源涵養(yǎng)與土壤保持區(qū),承擔著重要的生態(tài)系統(tǒng)服務(wù)功能,從華北平原西望太行山脈,東側(cè)懸崖斷壁處處可見. 太行山東緣大斷裂造成的與華北平原的正負極地形,使太行山脈成為華北平原區(qū)重要的天然屏障帶. 晉中市處于太行山脈與太原盆地之間,榆次區(qū)、太谷區(qū)、祁縣、平遙縣、靈石縣西部及介休市西北部處于太原盆地東南側(cè),為晉中市平原區(qū),其余縣都處在晉中市的山地丘陵區(qū).
圖1 研究區(qū)地理位置Fig.1 Geographical location of the study area
研究所用的DEM數(shù)據(jù)來源于地理空間數(shù)據(jù)云(http://www.gscloud.cn/),數(shù)據(jù)分辨率為30 m;凈初級生產(chǎn)力(Net Primary Productivity,NPP)數(shù)據(jù)為2000—2015 年的MODIS 陸地4 級標準數(shù)據(jù)NPP 產(chǎn)品數(shù)據(jù)集MOD17A3H,來自美國航空航天局網(wǎng)站(https://ladsweb.modaps.eosdis.nasa.gov/),其空間分辨率為500 m;土壤數(shù)據(jù)來源于世界土壤數(shù)據(jù)庫中的中國土壤數(shù)據(jù)集(Harmonized World Soil Database v1.2);氣象數(shù)據(jù)來源于中國科學院資源環(huán)境科學與數(shù)據(jù)中心(https://www.resdc.cn/);研究區(qū)矢量數(shù)據(jù)主要包括晉中市自然保護區(qū)、森林公園等關(guān)鍵生態(tài)區(qū)分布圖、晉中市地貌圖及晉中市1∶100萬地質(zhì)災害易發(fā)程度分區(qū)圖、土地利用變更調(diào)查數(shù)據(jù)庫等,數(shù)據(jù)來源于晉中市自然資源局. 此外還有以上數(shù)據(jù)派生的相關(guān)輔助數(shù)據(jù). 應用ArcGIS軟件將數(shù)據(jù)投影變換為CGCS2000_3_Degree_GK_Zone_38平面坐標系,同時將所有數(shù)據(jù)重采樣為500 m空間分辨率.
地形地貌作為自然地理環(huán)境中的基底要素,直接影響或主導著區(qū)域生態(tài)系統(tǒng)要素的分布與變化. 研究區(qū)居于黃土高原東部邊緣,東部有太行山脈作為天然生態(tài)屏障,西部以汾河為塹,山地、丘陵、平川呈階梯狀分布. 基于數(shù)字高程模型、地貌圖、地形起伏度等數(shù)據(jù),對研究區(qū)地貌信息進行精細化提取,并將研究區(qū)地貌劃分為山地、丘陵、平原、臺地、階地5類. 在此基礎(chǔ)上,以地貌類型作為基底要素,結(jié)合研究區(qū)生態(tài)系統(tǒng)服務(wù)功能重要性和生態(tài)敏感性的空間分布特征對研究區(qū)的生態(tài)功能區(qū)進行劃分.
根據(jù)《生態(tài)保護紅線劃定指南》[19]中推薦的方法,結(jié)合研究區(qū)自然地理特征及生態(tài)環(huán)境狀況,確定研究區(qū)生態(tài)系統(tǒng)服務(wù)功能重要性的評價指標為:水源涵養(yǎng)功能重要性、土壤保持功能重要性、固碳釋氧功能重要性和生物多樣性保護功能重要性,其中水源涵養(yǎng)功能重要性和生物多樣性保護功能重要性采用NPP定量指標評估法,分別以生態(tài)系統(tǒng)水源涵養(yǎng)服務(wù)能力指數(shù)和生物多樣性維護服務(wù)能力指數(shù)作為評估指標;土壤保持功能重要性采用修正通用的水土流失方程模型評估法;固碳釋氧功能重要性通過光合作用方程[20]計算得到. 采用變異系數(shù)法[21]對各單項評價指標賦予權(quán)重,分別為0.242 7、0.151 6、0.358 5和0.247 2.
研究區(qū)處于華北平原黃土高原區(qū)和北方農(nóng)牧交錯生態(tài)脆弱區(qū),市域大部分地貌溝壑縱橫,水土流失、地面塌陷等生態(tài)問題頻現(xiàn),客觀上造成生態(tài)區(qū)敏感性較高. 參考相關(guān)研究[20,22]并結(jié)合研究區(qū)實際情況,選取水土流失、地質(zhì)災害、水資源環(huán)境和生境四個評價因子對研究區(qū)生態(tài)敏感性進行評價. 利用變異系數(shù)法確定四個評價因子的權(quán)重,采用綜合指數(shù)法計算各評價指標的指數(shù),生態(tài)敏感性評價因子及其分級標準如表1所示.
表1 生態(tài)敏感性評價因子及分級標準Tab.1 Evaluation factors and grading standards of ecological sensitivity
參考相關(guān)研究[23],利用三維魔方空間分類法的基本原理[24],以生態(tài)系統(tǒng)服務(wù)功能重要性分類為X軸,生態(tài)敏感性分類為Y軸,地貌類型為Z軸建立生態(tài)功能分區(qū)三維魔方圖(圖2). 生態(tài)系統(tǒng)服務(wù)功能重要性分類中:X1為不重要區(qū),X2為一般重要區(qū),X3為中等重要區(qū),X4為高度重要區(qū),X5為極重要區(qū);生態(tài)敏感性分類中:Y1為不敏感區(qū),Y2為一般敏感區(qū),Y3為中度敏感區(qū),Y4為高度敏感區(qū),Y5為極敏感區(qū);地貌類型中:Z1為平原盆地區(qū),Z2為黃土丘陵階地區(qū),Z3為高中山丘陵區(qū),Z4為山間盆地區(qū),Z5為中山丘陵區(qū).
圖2 生態(tài)功能分區(qū)三維魔方圖Fig.2 Three-dimensional Rubik’s Cube diagram ofecological function zoning
通過三維魔方空間分類法將X、Y、Z軸中的不同類型進行組合,得到62種主要組合類型,在考慮生態(tài)系統(tǒng)類型和區(qū)劃邊界完整性的基礎(chǔ)上,最終將研究區(qū)劃分為5個生態(tài)功能區(qū),分別為汾河平川水土流失生態(tài)恢復區(qū)、黃土丘陵階地生態(tài)農(nóng)業(yè)發(fā)展區(qū)、太岳山地生物資源涵養(yǎng)區(qū)、山間盆地水土保持生態(tài)保育區(qū)和太行山地水源涵養(yǎng)生態(tài)屏障區(qū),如表2所示.
表2 研究區(qū)生態(tài)功能區(qū)與三維魔方單元對應表Tab.2 Corresponding table of ecological function zones and three-dimensional Rubik’s Cube unit in the study area
3.1.1 水源涵養(yǎng)功能重要性分析
研究區(qū)水源涵養(yǎng)功能重要性類型以高度重要為主,如圖3(a)所示. 水源涵養(yǎng)功能極重要區(qū)和高度重要區(qū)的面積分別為1 732.74 km2和4 220.35 km2,分別占研究區(qū)總面積的10.57%和25.73%,主要分布在太岳山—中條山水源涵養(yǎng)帶和太行山水源涵養(yǎng)帶上,這些地區(qū)的NPP和多年平均降水量均較高,這是其維持水源涵養(yǎng)功能的重要保證. 水源涵養(yǎng)功能不重要區(qū)、一般重要區(qū)和中等重要區(qū)主要分布在平川以及汾漳河支流周圍地帶,這些地區(qū)植被覆蓋少,人類活動干擾較大.
3.1.2 土壤保持功能重要性分析
研究區(qū)的土壤保持功能重要性類型以不重要為主,如圖3(b)所示. 土壤保持功能不重要區(qū)的面積約為8 274.62 km2,占研究區(qū)總面積的50.46%,地表覆蓋類型以耕地、水域和建設(shè)用地為主,水域河道周圍侵蝕較嚴重且植被覆蓋度較小,土壤保持能力不高. 土壤保持功能極重要區(qū)和高度重要區(qū)的總面積占研究區(qū)總面積的10.64%,主要分布在太行山脈東側(cè),這里的土地利用類型雖然以林地為主,但土壤類型主要為砂質(zhì)壤土、褐土和粗骨土等,土壤可侵蝕性較強,并且地形起伏度較大,極易發(fā)生由降水以及地表徑流所帶來的土壤侵蝕問題.
3.1.3 固碳釋氧功能重要性分析
研究區(qū)的固碳釋氧功能重要性總體處于較高水平,如圖3(c)所示. 固碳釋氧功能極重要區(qū)和高度重要區(qū)基本分布在植被覆蓋較好的山脈地區(qū),主要集中在太岳山脈北段及太行山脈中段. 固碳釋氧功能極重要區(qū)和高度重要區(qū)的總面積占研究區(qū)總面積的70.82%,可見研究區(qū)固碳釋氧功能的重要性.
3.1.4 生物多樣性保護功能重要性分析
研究區(qū)的生物多樣性保護功能重要性類型以中度重要為主,如圖3(d)所示. 生物多樣性保護功能中度重要區(qū)面積約為4 510.14 km2,占研究區(qū)總面積的27.50%,主要位于一般重要區(qū)的延伸周邊和河流緩沖帶周圍. 生物多樣性保護功能極重要區(qū)和高度重要區(qū)主要分布在研究區(qū)的中部及南部地區(qū),這些區(qū)域以山地為主,植被覆蓋良好,涵蓋了太岳山國家森林公園、清漳河國家濕地公園等,生物資源豐富.
圖3 研究區(qū)生態(tài)系統(tǒng)服務(wù)功能重要性評價結(jié)果Fig.3 Evaluation results of ecosystem service function importance in the study area
3.1.5 生態(tài)系統(tǒng)服務(wù)功能重要性綜合分析
研究區(qū)生態(tài)系統(tǒng)服務(wù)功能重要性綜合評價結(jié)果如圖3(e)所示.由圖3(e)及以上分析結(jié)果可知,研究區(qū)的生態(tài)系統(tǒng)服務(wù)功能以水源涵養(yǎng)和固碳釋氧為主,重要性類型以高度重要為主. 生態(tài)系統(tǒng)服務(wù)功能高度重要區(qū)的面積約為4 226.94 km2,占研究區(qū)總面積的25.77%. 生態(tài)系統(tǒng)服務(wù)功能極重要區(qū)與高度重要區(qū)鑲嵌分布,主要分布在研究區(qū)的中部以及太行山脈—太岳山脈主脊兩側(cè),部分分布在研究區(qū)西北部的黃土丘陵區(qū),這些區(qū)域內(nèi)林、草植被生態(tài)系統(tǒng)較為豐富,地貌類型以山地丘陵為主,是山西省重要的水源涵養(yǎng)區(qū)和水土保持區(qū),對維護太行山脈以及黃河、海河流域的生態(tài)系統(tǒng)服務(wù)功能具有重要作用,同時也是晉中地區(qū)重要的物種貯存庫.
3.2.1 水土流失敏感性分析
研究區(qū)的水土流失敏感性分析結(jié)果如圖4(a)所示. 由圖4(a)可知,研究區(qū)水土流失敏感性的空間分異性較大,水土流失極敏感區(qū)和高度敏感區(qū)的總面積占研究區(qū)總面積的34.50%,水土流失嚴重的地區(qū)主要分布在研究區(qū)各支流流域范圍內(nèi),大部分為土石山區(qū),多為土山坡,地形起伏較大,致使地表水源涵養(yǎng)能力下降,水土流失嚴重. 水土流失不敏感區(qū)多分布在植被覆蓋良好的區(qū)域,土地利用類型主要為林地和草地.
圖4 研究區(qū)生態(tài)敏感性評價結(jié)果Fig.4 Evaluation results of ecological sensitivity in the study area
3.2.2 地質(zhì)災害敏感性分析
研究區(qū)地處黃土高原東部邊緣,地貌類型復雜,境內(nèi)地質(zhì)環(huán)境脆弱. 研究區(qū)地質(zhì)災害敏感性類型主要以中度敏感為主,如圖4(b)所示. 地質(zhì)災害中度敏感區(qū)的面積約為4 738.8 km2,占研究區(qū)總面積的28.90%,主要分布在研究區(qū)的中南部區(qū)域、西北部黃土高原丘陵區(qū)和靈石縣境內(nèi),這些區(qū)域內(nèi)煤層分布較多,加之坡度和地貌類型的影響,易發(fā)生地質(zhì)災害. 地質(zhì)災害極敏感區(qū)的面積為2 233.75 km2,占研究區(qū)總面積的13.62%,基本上分布在太岳山和太行山周邊的山體附近,這些區(qū)域內(nèi)多高峻的山地,地勢陡峭,地質(zhì)災害易發(fā)程度高.
3.2.3 水資源環(huán)境敏感性分析
從距河湖距離和水體分布兩方面對研究區(qū)的水資源環(huán)境敏感性進行分析,結(jié)果如圖4(c)所示. 研究區(qū)的水資源環(huán)境敏感性類型以不敏感為主,水資源環(huán)境極敏感區(qū)和高度敏感區(qū)主要分布在研究區(qū)的汾、漳河水系附近. 水資源環(huán)境極敏感區(qū)和中度敏感區(qū)的面積相當,分別為1 034.84 km2和992.71 km2,共占研究區(qū)總面積的12.36%. 在人類活動的干擾下,研究區(qū)內(nèi)平原地區(qū)的水資源環(huán)境敏感程度較其他地區(qū)嚴重,水資源匱乏且河流水質(zhì)污染較重,河流地表徑流量不斷減少.
3.2.4 生境敏感性分析
研究區(qū)生境敏感性類型以一般敏感為主,如圖4(d)所示. 生境一般敏感區(qū)的面積約為9 367.73 km2,占研究區(qū)總面積的57.12%,分布范圍廣,且主要集中在研究區(qū)中部及東部山地丘陵區(qū),土地利用類型以林地和草地為主. 生境極敏感區(qū)和高度敏感區(qū)基本分布在國家級、省級自然保護區(qū)、森林公園和風景名勝區(qū),均位于研究區(qū)生態(tài)保護紅線關(guān)鍵區(qū)內(nèi),敏感性較高. 生境不敏感區(qū)面積為5 295.25 km2,占研究區(qū)總面積的32.29%,較廣泛地分布在太原盆地東側(cè)及地勢低平區(qū),這些區(qū)域距離建設(shè)用地和耕地較近,受人類活動干擾較大,生態(tài)基底較差,不適合生物多樣性的維持以及相關(guān)生態(tài)系統(tǒng)服務(wù)功能的發(fā)揮.
3.2.5 生態(tài)敏感性綜合分析
研究區(qū)的生態(tài)敏感性綜合評價結(jié)果如圖4(e)所示. 由圖4(e)及以上分析結(jié)果可知,地質(zhì)災害和水土流失是研究區(qū)的主要生態(tài)敏感性類型,生態(tài)極敏感區(qū)和高度敏感區(qū)的總面積為4 786.99 km2,占研究區(qū)總面積的29.19%,主要分布在汾、漳兩河及其支流的源頭山地,以及太行山中段區(qū)域,這些地區(qū)是晉中市生態(tài)保護紅線的核心地帶,屬生境、水土流失和地質(zhì)災害較敏感的地區(qū),對維持市域和太行山脈—太岳山脈區(qū)域的生態(tài)系統(tǒng)穩(wěn)定及其生態(tài)系統(tǒng)服務(wù)功能的發(fā)揮具有重要作用. 生態(tài)一般敏感區(qū)和不敏感區(qū)的總面積占研究區(qū)總面積的41.30%,主要分布在研究區(qū)的西北部丘陵階地區(qū)、太原盆地東側(cè)地區(qū),部分分布在研究區(qū)中東部地區(qū).
生態(tài)功能分區(qū)是結(jié)合區(qū)域現(xiàn)存的生態(tài)環(huán)境問題,依據(jù)其不同屬性特征的地域分異特點以及生態(tài)系統(tǒng)服務(wù)功能重要性和生態(tài)敏感性的空間分布規(guī)律,劃分出能夠反映區(qū)域生態(tài)系統(tǒng)功能差異特征的空間單元[25]. 根據(jù)研究區(qū)生態(tài)系統(tǒng)服務(wù)功能重要性和生態(tài)敏感性的空間分布特征,結(jié)合其地形地貌及生態(tài)系統(tǒng)完整性確定的自然邊界,通過邊界處理劃定其生態(tài)功能區(qū)范圍,如圖5所示,各生態(tài)功能區(qū)的特征、主要生態(tài)環(huán)境問題及保護措施如表3所述.
表3 研究區(qū)生態(tài)功能區(qū)的特征、主要生態(tài)環(huán)境問題及保護措施Tab.3 Characteristics,main ecological environment problems and protection measures of ecological function zones in the study area
圖5 研究區(qū)生態(tài)功能分區(qū)結(jié)果Fig.5 The results of ecological function zoning in the study area
以生態(tài)功能為導向,對晉中市的地貌類型、生態(tài)系統(tǒng)服務(wù)功能重要性和生態(tài)敏感性進行了分析,并對其生態(tài)功能區(qū)進行了劃分. 得出結(jié)論如下:
1)晉中市生態(tài)系統(tǒng)服務(wù)功能以水源涵養(yǎng)和固碳釋氧為主,生態(tài)敏感性以地質(zhì)災害和水土流失為主.
2)晉中市生態(tài)系統(tǒng)服務(wù)功能極重要區(qū)和生態(tài)極敏感區(qū)主要分布在晉中市自然保護區(qū)等生境質(zhì)量好的區(qū)域.
3)根據(jù)晉中市地形地貌的空間異質(zhì)性特征以及晉中市生態(tài)系統(tǒng)服務(wù)功能重要性和生態(tài)敏感性的空間分布特征,利用三維魔方空間分類法將晉中市劃分為5 個生態(tài)功能區(qū),分別是汾河平川水土流失生態(tài)恢復區(qū)、黃土丘陵階地生態(tài)農(nóng)業(yè)發(fā)展區(qū)、太岳山地生物資源涵養(yǎng)區(qū)、山間盆地水土保持生態(tài)保育區(qū)和太行山地水源涵養(yǎng)生態(tài)屏障區(qū),并根據(jù)各生態(tài)功能區(qū)的特征針對性地提出了相應的保護措施.
生態(tài)系統(tǒng)是由生物有機體和周圍環(huán)境構(gòu)成的生態(tài)綜合體,但自然要素及土地利用類型的不同導致不同區(qū)域內(nèi)會形成不同類型的生態(tài)系統(tǒng)組合,且生態(tài)系統(tǒng)的空間分布呈現(xiàn)明顯的區(qū)域分異性[1]. 對不同類型的生態(tài)系統(tǒng)組合進行生態(tài)功能分區(qū),有助于構(gòu)建區(qū)域生態(tài)安全格局、完善國土開發(fā)格局空間體系,并且可以為生態(tài)功能區(qū)劃定方法的完善提供確切的樣本案例. 在借鑒前人研究的基礎(chǔ)上,基于“地貌-重要性-敏感性”的分析框架,對晉中市生態(tài)功能區(qū)進行劃分,不僅表現(xiàn)了晉中市生態(tài)系統(tǒng)服務(wù)功能重要性的空間分布特征,而且也體現(xiàn)了其生態(tài)敏感性的空間分布特征,同時以地貌類型作為其生態(tài)功能分區(qū)的基底要素,將自然地理特征與生態(tài)過程、人類社會經(jīng)濟活動的干擾等因素相結(jié)合,體現(xiàn)了人地矛盾關(guān)系中生態(tài)環(huán)境的動態(tài)交互特征. 但研究所采用的《生態(tài)保護紅線劃定指南》中所推薦的評估模型在綜合性上也有不足之處,相關(guān)分級要求可能更適用于省域、流域等大尺度區(qū)域,數(shù)據(jù)空間精度要求較低,這在一定程度上影響了分析結(jié)果的準確性. 未來可以考慮選用更加客觀的定量模型方法或者選取更多能夠反映自然生態(tài)過程與人類干擾作用下生態(tài)斑塊在生態(tài)系統(tǒng)中的空間重要性及生態(tài)敏感性的指標,提高數(shù)據(jù)精度,同時還可以考慮與區(qū)域發(fā)展戰(zhàn)略相結(jié)合,從宏觀和微觀的角度深入探討區(qū)域生態(tài)功能分區(qū)的保護以及其變化驅(qū)動機理,以進一步優(yōu)化生態(tài)功能區(qū)的劃定方法.