劉軍會,馬 蘇,高吉喜,鄒長新,王晶晶,劉志強,王麗霞 (.中國環(huán)境科學研究院,北京 000；.環(huán)境保護部南京環(huán)境科學研究所,江蘇 南京 00；.生態(tài)環(huán)境部環(huán)境規(guī)劃院,北京 000；.中國人民大學,北京 0087)

生態(tài)保護紅線是指在生態(tài)空間范圍內具有特殊重要生態(tài)功能、必須強制性嚴格保護的區(qū)域.圍繞生態(tài)保護紅線劃定方法,國內學者以省、市、縣等行政區(qū)為研究對象開展了大量探索性研究.參照《生態(tài)保護紅線劃定指南》[1],全國31個省(自治區(qū)、直轄市)組織專家學者開展了本轄區(qū)的生態(tài)保護紅線劃定工作;此外,市級[2-7]、縣級[8-10]行政區(qū)生態(tài)保護紅線劃定研究也取得了大量結果.以行政區(qū)劃為研究對象開展的生態(tài)保護紅線劃定,便于數(shù)據(jù)收集和后期管理;但由于劃定主體不同、對紅線認識不同,加之同級行政區(qū)范圍差異大,評價數(shù)據(jù)精度、時效性不一致,造成跨行政區(qū)的紅線分布存在明顯的不連續(xù)、空白或邊界錯亂現(xiàn)象,割裂了山脈、河流、濕地、植被等自然地理要素的空間完整性,難以形成生態(tài)保護紅線全國“一張圖”.為打破行政區(qū)劃制約,應根據(jù)我國地域分異特征,將區(qū)域、流域作為一個整體,統(tǒng)一劃定生態(tài)保護紅線,并在省級層面進行銜接,確?？缧姓^(qū)的紅線邊界連續(xù)、完整.

京津冀地理位置相鄰,同處海河流域,山相連、水相通.其中,北京、天津地區(qū)80%以上的用水來自于冀北的桑干河、洋河、潮白河和灤河,行政區(qū)域之間生態(tài)關系緊密.京津冀地區(qū)生態(tài)地位重要,是華北平原重要的生態(tài)屏障.燕山山地是灤河、潮白河等河流的發(fā)源地,具有重要的水源涵養(yǎng)功能;太行山山地是海河水系的發(fā)源地和全國水土流失重點治理區(qū),具有重要的水源涵養(yǎng)和水土保持功能;冀北壩上高原是風沙逼近北京、天津的主要通道和重要沙源地之一,具有重要的防風固沙功能.同時,京津冀地區(qū)生態(tài)環(huán)境非常脆弱,是我國東部地區(qū)資源環(huán)境與發(fā)展矛盾最為突出的地區(qū),土地沙化、水土流失、生態(tài)空間退縮、城市熱島效應等生態(tài)問題突出,當前生態(tài)安全形勢依然十分嚴峻.

京津冀協(xié)同發(fā)展已經上升為國家重大發(fā)展戰(zhàn)略,核心是京津冀三地作為一個整體協(xié)同發(fā)展.統(tǒng)一劃定京津冀區(qū)域生態(tài)保護紅線是促進京津冀協(xié)同發(fā)展的重要組成,有利于在區(qū)域層面實現(xiàn)對京津冀生態(tài)空間的統(tǒng)一管控.本研究以京津冀地區(qū)陸域作為一個整體,利用遙感和GIS技術,搭建統(tǒng)一的工作平臺,建立評價指標和評價模型,開展區(qū)域生態(tài)系統(tǒng)服務重要性和生態(tài)敏感性評價,確立分級賦值標準,定量揭示京津冀地區(qū)生態(tài)系統(tǒng)服務重要性和生態(tài)敏感性的空間分布特征;再利用全國第二次土地調查等明確的地塊邊界,采用交互矢量化方法,將生態(tài)系統(tǒng)服務極重要區(qū)和生態(tài)環(huán)境極敏感區(qū)劃為生態(tài)保護紅線.統(tǒng)一劃定京津冀地區(qū)生態(tài)保護紅線,有利于從區(qū)域層面構建生態(tài)安全格局、推動重要生態(tài)空間協(xié)調保護和實施統(tǒng)一監(jiān)管.

1 數(shù)據(jù)獲取與處理

氣象數(shù)據(jù):覆蓋京津冀地區(qū) 2000～2015年降雨量、溫度、風速和大風日數(shù)等資料來源于中國氣象科學數(shù)據(jù)共享服務網中的中國地面氣候資料數(shù)據(jù)集(http://cdc.cma.gov.cn/).氣象數(shù)據(jù)包括研究區(qū)各氣象站點的編號、經緯度和海拔,以及每個氣象站點在相應分析時間尺度內的降雨量(0.1mm)、溫度(0.1)℃、風速(0.1m/s).根據(jù)各氣象站點的經緯度信息,采用ArcGIS的Kriging空間插值法,獲取研究區(qū)的氣象數(shù)據(jù)柵格圖像.

遙感數(shù)據(jù):研究區(qū)2000～2015年歸一化植被指數(shù)(NDVI)來源于美國國家航空航天局(NASA)的EOS/MODIS數(shù)據(jù)產品(http://e4ft101.cr.usgs.gov/),空間分辨率為 250m,時間分辨率為 16d.為消除異常值的影響,采用最大合成法(MVC)合成月NDVI數(shù)據(jù),再將月NDVI合成年最大NDVI影像[10],獲取京津冀地區(qū)2000～2015年逐年及逐月NDVI影像.

高程數(shù)據(jù):數(shù)字高程模型(DEM)來自于地理空間數(shù)據(jù)云網站(http://www.gscloud.cn),空間分辨率 30m,用以提取坡度等地形數(shù)據(jù).

土壤數(shù)據(jù):土壤數(shù)據(jù)來源于中國科學院南京土壤研究所的中國1:100萬土壤數(shù)據(jù)庫,屬性數(shù)據(jù)庫包括土壤表層黏粒、粉沙、沙粒、有機質含量等屬性字段.

其他數(shù)據(jù):2015年京津冀地區(qū)的土地利用數(shù)據(jù)、植被NPP數(shù)據(jù)來源于全國生態(tài)環(huán)境現(xiàn)狀調查數(shù)據(jù)庫;第二次全國土地調查數(shù)據(jù)來源于國土資源部.

為統(tǒng)一京津冀地區(qū)評價尺度,本研究中所用數(shù)據(jù)均統(tǒng)一轉換成阿爾伯斯等面積圓錐投影,采用克拉索夫斯基橢球體,評價柵格單元為250m×250m.

2 研究方法

2.1 生態(tài)評價方法.

依據(jù)京津冀地區(qū)生態(tài)環(huán)境現(xiàn)狀、主導生態(tài)系統(tǒng)服務功能定位和自然環(huán)境演變過程中出現(xiàn)的生態(tài)問題,選取水源涵養(yǎng)、水土保持、防風固沙3類生態(tài)系統(tǒng)服務以及土地沙化、水土流失2類典型生態(tài)問題,構建評價指標和評價模型,開展區(qū)域生態(tài)系統(tǒng)服務重要性和生態(tài)敏感性評價,定量揭示京津冀地區(qū)生態(tài)系統(tǒng)服務重要性和生態(tài)敏感性的空間分布特征.

2.1.1 水源涵養(yǎng)服務重要性評價 水源涵養(yǎng)是生態(tài)系統(tǒng)通過其特有的結構與水相互作用,對降水進行截留、滲透、蓄積,并通過蒸散發(fā)實現(xiàn)對水流、水循環(huán)的調控.以水源涵養(yǎng)量作為生態(tài)系統(tǒng)水源涵養(yǎng)重要性評價指標,采用水量平衡方程來計算水源涵養(yǎng)量,計算公式為:

式中: TQ為總水源涵養(yǎng)量,m3; Pi為降雨量,mm; Ri為地表徑流量,mm; ETi為蒸散發(fā),mm; Ai為i類生態(tài)系統(tǒng)面積,km2; i為研究區(qū)第i類生態(tài)系統(tǒng)類型, j為研究區(qū)生態(tài)系統(tǒng)類型數(shù).降雨量因子、蒸散發(fā)因子柵格圖像根據(jù)各氣象站點的經緯度信息,采用 ArcGIS的Kriging空間插值法獲取.地表徑流量由多年平均降雨量乘以地表徑流系數(shù)獲得,各類型生態(tài)系統(tǒng)地表徑流系數(shù)根據(jù)經驗值確定.生態(tài)系統(tǒng)面積因子根據(jù)全國生態(tài)狀況遙感調查與評估成果中的生態(tài)系統(tǒng)類型數(shù)據(jù)集提取.

2.1.2 水土保持服務重要性評價 水土保持是生態(tài)系統(tǒng)通過其結構與過程減少由于水蝕所導致的土壤侵蝕的作用.以水土保持量,即潛在土壤侵蝕量與實際土壤侵蝕量的差值,作為生態(tài)系統(tǒng)水土保持重要性的評價指標.采用修正通用水土流失方程(RUSLE)計算水土保持量,公式如下:

式中: Ac為水土保持量,t/(hm2?a); Ap為潛在土壤侵蝕量; Ar為實際土壤侵蝕量;R為降雨侵蝕力因子,MJ?mm/(hm2?h?a); K 為 土 壤 可 蝕 性 因 子 ,t?hm2?h/(hm2?MJ?mm); L、S 為地形因子,L 表示坡長因子,S 表示坡度因子;C為植被覆蓋因子.降雨侵蝕力因子根據(jù)王萬忠等[12]利用降水資料計算的中國 100多個城市的R值,采用ArcGIS的內插法繪制研究區(qū)R柵格分布圖.采用EPIC模型[13-14]計算K值,該模型在土壤侵蝕預測模塊中采用土壤有機碳和粒徑組成資料來估算土壤可蝕性K值;坡度、坡長是反映地形對土壤侵蝕影響的 2個因子,坡度是主要影響因子,一般來說,在其他條件相同的情況下,＞35°的坡地土壤侵蝕量最大[15],坡度根據(jù)研究區(qū)數(shù)字高程模型(DEM)在GIS系統(tǒng)下直接生成.植被覆蓋度信息提取是在對光譜信號進行分析的基礎上,通過建立歸一化植被指數(shù)與植被覆蓋度的轉換關系,直接提取植被覆蓋信息[16].

2.1.3 防風固沙服務重要性評價 防風固沙是生態(tài)系統(tǒng)通過防風、固沙、阻沙以及對風沙土的改良作用等途徑對地表土壤形成保護,減少風蝕輸沙量.以潛在風蝕量與實際風蝕量的差值作為生態(tài)系統(tǒng)減少風蝕的防風固沙量.采用董治寶[17]提出的風蝕流失模型計算防風固沙量,公式為:

式中: Q為防風固沙量,t; θ為坡度,°; V為風速,m/s;VCR為植被覆蓋度,%; SDR為人為地表結構破損率,%;H為空氣相對濕度,%;d為土壤平均粒徑,mm; F為土體硬度,N/cm2; t為時間,s; x,y為距離參照點距離,km.風速、空氣相對濕度柵格圖由氣象站點數(shù)據(jù)空間插值生成;坡度圖由DEM提取;人為地表結構破損率不考慮,設為 1;土壤平均粒徑、土體硬度由土壤數(shù)據(jù)庫提取;時間(s)為 2000～2015年揚沙天數(shù)年平均值;植被覆蓋度根據(jù)研究區(qū)NDVI提取.

2.1.4 水土流失敏感性評價 根據(jù)通用水土流失方程的基本原理,選取降水侵蝕力、土壤可蝕性、坡度坡長和植被覆蓋等作為評價指標.將反映各因素對水土流失敏感性的單因子評估數(shù)據(jù),用GIS技術進行乘積運算,公式如下:

式中: SSi為i空間單元水土流失敏感性指數(shù),Ri、Ki、LSi、Ci分別為降雨侵蝕力、土壤可蝕性、坡長坡度、植被覆蓋度.指標計算方法同上.

2.1.5 土地沙化敏感性評價 土地沙化主要受氣候干燥程度影響表現(xiàn)為氣候干燥,植物生長困難,地表植被覆蓋度低.干燥氣候減緩了地表土壤形成過程,使地表結構分散,易受風蝕,土壤的抗風蝕沙漠化能力與土壤水分含量成正比[18-19].選取干燥度指數(shù)、起沙風天數(shù)、土壤質地、植被覆蓋度作為土地沙化敏感性評價指標,公式如下:

式中: Di為i評估區(qū)域土地沙化敏感性指數(shù); Ii、Wi、Ki、Ci分別為評估區(qū)域干燥度指數(shù)、起沙風天數(shù)、土壤質地和植被覆蓋.干燥度指數(shù)采用修正的謝良尼諾夫公式[20]計算;根據(jù)劉連友等[21]研究結果,沙質壤土、壤質沙土和固定風沙土的起動風速分別為 6.0、6.6和5.1m/s,本研究選用冬春季節(jié)大于6m/s起沙風天數(shù)指標評價土地沙化敏感性;土壤質地數(shù)據(jù)通過 1:100萬土壤數(shù)據(jù)獲取.

2.2 分級賦值標準

在 ArcGIS軟件中,運用柵格計算器,輸入公式“Int([某一要素評估的柵格數(shù)據(jù)]/[某一要素評估柵格數(shù)據(jù)的最大值]×100)”,得到歸一化后的生態(tài)系統(tǒng)服務重要性、生態(tài)敏感性評價柵格圖;再利用ArcGIS提供的自然斷點分類方法,將生態(tài)系統(tǒng)服務重要性評價結果劃分為一般重要、輕度重要、中度重要、高度重要和極重要5個等級,將生態(tài)敏感性評價結果劃分為一般敏感、輕度敏感、中度敏感、高度敏感和極敏感5個等級.具體分級賦值標準見表1.

表1 生態(tài)系統(tǒng)服務重要性、生態(tài)敏感性評價分級賦值標準Table 1 Classification of ecosystem service importance and ecological environmental sensitivity

2.3 紅線邊界確定

生態(tài)保護紅線劃定是在綜合考慮區(qū)域生態(tài)系統(tǒng)服務重要性、生態(tài)敏感性評價基礎上,將生態(tài)系統(tǒng)服務極重要區(qū)和生態(tài)環(huán)境極敏感區(qū)劃為生態(tài)保護紅線.為使生態(tài)評價結果落到實地,精確界定生態(tài)保護紅線邊界,需經過以下步驟處理:首先,利用ArcGIS中的聚合工具將鄰近的圖斑聚合為相對完整連片的大塊圖斑,聚合距離為 250m,最小孔洞大小為 1km2;其次,將面積小于 1km2的獨立圖斑扣除,降低斑塊破碎化程度;最后,本研究利用全國第二次土地調查數(shù)據(jù)(該數(shù)據(jù)明確了每個地塊的類型、面積、權屬和分布信息,比例尺較大,準確率較高)與綜合評價圖進行空間疊加,對具體地塊進行人工判讀和補充勾繪,通過尺度轉換,將低分辨率的評價柵格數(shù)據(jù)轉換為高分辨率的生態(tài)保護紅線矢量數(shù)據(jù),形成邊界清晰、生態(tài)系統(tǒng)相對完整和易于管理的生態(tài)保護紅線.

3 結果與討論

3.1 劃定結果

圖1 水源涵養(yǎng)服務極重要區(qū)Fig.1 The extremely important area of water conservation service

京津冀地區(qū)水源涵養(yǎng)服務極重要區(qū)面積5.17萬km2,占研究區(qū)總面積的23.85%,主要包括永定河、潮白河、灤河、遼河和海河等主要河流源頭區(qū),以及密云水庫、官廳水庫、懷柔水庫、北大港水庫、于橋水庫、東七里海水庫、團泊洼水庫等飲用水水源保護區(qū)和水源涵養(yǎng)區(qū)(圖 1).水土保持服務極重要區(qū)面積4.87萬km2,占研究區(qū)總面積的22.48%,主要分布在西部太行山山地、永定河上游的間山盆地區(qū)、燕山山地以及灤河上游山川河谷地區(qū)等(圖 2).防風固沙服務極重要區(qū)面積 1.05萬 km2,占研究區(qū)總面積的4.84%,主要分布冀北塞罕壩、冀西北間山盆地(圖3).水土流失極敏感區(qū)面積0.99萬km2,占研究區(qū)總面積的 4.56%,主要分布在保定市和石家莊市西部的太行山地、冀西北間山盆地以及北部燕山山地的部分區(qū)域(圖 4).土地沙化極敏感區(qū)面積 0.57萬 km2,占研究區(qū)總面積的 2.62%,主要分布在張家口市的壩上高原地區(qū)和冀西北間山盆地(圖5).

圖2 水土保持服務極重要區(qū)Fig.2 The extremely important area of soil conservation service

圖3 防風固沙服務極重要區(qū)Fig.3 The extremely important area of sand fixation service

圖4 水土流失極敏感區(qū)Fig.4 The extremely sensitive area of soil erosion

圖5 土地沙化極敏感區(qū)Fig.5 The extremely sensitive area of land desertification

京津冀地區(qū)生態(tài)系統(tǒng)服務極重要區(qū)、生態(tài)環(huán)境極敏感區(qū)的總面積為 12.65萬 km2,扣除二者相互重疊部分,生態(tài)保護紅線面積為6.68萬km2,占研究區(qū)總面積的30.9%.空間上,生態(tài)保護紅線總體呈“兩屏一帶”分布格局(圖6),其中“兩屏”為北部燕山和西部太行山生態(tài)屏障,主要生態(tài)功能為水源涵養(yǎng)和水土保持;“一帶”為壩上防風固沙帶,主要生態(tài)功能為防風固沙.紅線范圍涵蓋了境內主要河流的源頭區(qū)、水土流失重點治理區(qū)和沙化土地治理區(qū).從北京市、天津市和河北省的跨行政區(qū)銜接看,位于相鄰省、市邊界的生態(tài)保護紅線銜接自然順暢,不存在明顯的不連續(xù)、空白或邊界錯亂區(qū)域,保證了區(qū)域生態(tài)系統(tǒng)完整性.從不同省市看,北京市生態(tài)保護紅線主要分布在北部燕山山地和西部太行山山地,生態(tài)功能定位為水源涵養(yǎng)和水土保持;天津市生態(tài)保護紅線主要分布在北部薊州山地丘陵區(qū)和中南部平原濕地,生態(tài)功能定位為水源涵養(yǎng)和生物多樣性維護;河北省生態(tài)保護紅線主要分布在東北部燕山山地、西部太行山山地和北部壩上防風固沙帶,主要生態(tài)功能為水源涵養(yǎng)、水土保持和防風固沙.

圖6 生態(tài)保護紅線分布Fig.6 Spatial distribution of ecological conservation redline in Beijing-Tianjin-Hebei Area

3.2 劃定效果分析

京津冀地區(qū) 2015年生態(tài)系統(tǒng)類型分析發(fā)現(xiàn),生態(tài)保護紅線范圍內的林地、草地面積分別為4.44萬km2和0.75萬 km2,二者面積合計占生態(tài)保護紅線總面積的 77.9%,比紅線之外的林地草地面積比例(22.2%)高55.7%.

植被凈初級生產力(NPP)表示綠色植物在單位面積、單位時間內通過光合作用將植被可吸收的太陽能轉換成植物有機碳量的能力[22-24].研究區(qū)NPP分析發(fā)現(xiàn),2015年生態(tài)保護紅線范圍內的植被 NPP平均值為 471.0gC/(m2.a),比紅線之外的研究區(qū)植被NPP 平均值(381.3gC/(m2?a))高 23.5%.

通過對研究區(qū)主導生態(tài)系統(tǒng)類型、面積比例和植被NPP分布特征分析可知,京津冀地區(qū)生態(tài)保護紅線范圍內的生態(tài)系統(tǒng)重要性遠高于紅線外的區(qū)域,劃定結果符合“以較小面積獲取較大服務”原則,適于作為京津冀地區(qū)最重要、需要嚴格管控的生態(tài)空間.

3.3 不確定性分析

本研究在基礎數(shù)據(jù)處理過程中主要采用了遙感影像解譯和GIS空間內插方法,研究結果精度會受到數(shù)據(jù)處理誤差的影響;斑塊聚合、獨立細小斑塊扣除、人工判讀與補充勾繪等處理精度,都會影響到生態(tài)保護紅線劃定結果.由于數(shù)據(jù)收集難度大,缺乏準確邊界,本研究未將自然保護區(qū)、風景名勝區(qū)、森林公園、河湖濕地等生態(tài)保護地納入生態(tài)保護紅線.此外,劃定結果尚需開展實地核查,并征求地方政府與相關領域專家意見,提高劃定方法的可操作性.

4 結論

4.1 京津冀地區(qū)山相連、水相通,將京津冀地區(qū)作為一個研究整體,打造統(tǒng)一工作平臺,共同劃定京津冀地區(qū)生態(tài)保護紅線,既保障了跨省市紅線銜接自然順暢和生態(tài)系統(tǒng)完整性,又可以為北京市、天津市、河北省生態(tài)保護紅線劃定提供參考,同時對于推動整個區(qū)域重要生態(tài)空間協(xié)調保護和實施統(tǒng)一監(jiān)管具有重要實踐意義.

4.2 京津冀地區(qū)生態(tài)保護紅線總面積 6.68萬 km2,占研究區(qū)總面積的30.9%;空間上,生態(tài)保護紅線總體呈“兩屏一帶”分布格局,其中“兩屏”為北部燕山和西部太行山生態(tài)屏障,“一帶”為壩上防風固沙帶;紅線范圍涵蓋了研究區(qū)境內主要河流源頭區(qū)、水土流失重點治理區(qū)和沙化土地治理區(qū).

4.3 生態(tài)保護紅線范圍內的森林、草地面積合計占比達 77.9%,比紅線之外的區(qū)域高 55.7%;生態(tài)保護紅線范圍內的植被 NPP平均值為 471.0gC/(m2.a),比紅線外的植被NPP平均值高23.5%,生態(tài)保護紅線劃定結果符合“以較小面積獲取較大服務”原則,適于作為京津冀地區(qū)最重要、需要嚴格管控的生態(tài)空間.

4.4 自上而下演繹評價和自下而上歸納協(xié)調相結合是劃定生態(tài)保護紅線的一種有效途徑.國家層面應加強生態(tài)保護紅線劃定頂層設計,依據(jù)我國地域分異特征,在全國劃分出若干個重點區(qū)域和流域,搭建統(tǒng)一工作平臺,制定統(tǒng)一劃定標準,分別提出京津冀地區(qū)、青藏高原、長江流域、黃河流域、珠江流域等重點區(qū)域、流域的生態(tài)保護紅線分布建議,并與相關的省(自治區(qū)、直轄市)生態(tài)保護紅線劃定結果進行銜接,提升紅線邊界確定的科學性.