達里諾爾湖流域土地利用變化的生態系統服務價值影響評價

周景博 薛 婕 吳 良 王 曉 夏 瑞

(1.中國人民大學環境學院，北京 100872；2.中國環境科學研究院，北京 100012；3.中國科學院地理科學與資源研究所，北京 100101)

土地利用受自然和人為因素雙重驅動，是人類-環境耦合系統演變的綜合體現[1]，也是一定時期政策作用的結果。土地利用變化直接改變了生態系統的結構和功能，具有廣泛的生態效應[2],[3]20,[4-5]。在區域生態系統服務功能長期實地監測不可行的情況下，土地利用類型是生態系統服務的重要代理指標[6]，可以此估算其價值變化[7],[8]1431，反映其生態效應，對靜態評價區域土地利用具有現實意義[9-10]。土地利用變化具有漸變性、持續性和規律性，長期來看，土地利用類型轉換體現了生態格局從一種均衡逐漸達到另一種均衡的過程，符合馬爾科夫過程特征[11-15]。生態系統服務價值是評價土地利用政策的重要指標，土地調控不僅應參考既有土地利用格局的生態影響，更應關注其長期影響[16]。應用馬爾科夫模型預測土地利用變化，可以動態評價土地利用的長期生態影響，為指導區域土地規劃提供量化依據。

已有相關研究從內容看多聚焦于東部發達地區和城市生態系統[8]1427,[17-19]，對典型生態脆弱區的研究相對較少[20]，更缺少對半干旱地區農牧交錯帶的研究；從方法看或側重于土地利用的生態系統服務價值的靜態計算[21]，或側重于土地利用的動態預測，少有結合靜態和動態分析的綜合評價；預測方法多采用指標體系或情景構建的方法[3]17,[22-24]，重視人為因素的調控作用，但忽視了土地利用的自然漸變和持續特征。本研究以達里諾爾湖流域為研究對象，采用1980、1990、2000、2010、2019年5期土地利用數據，計算了流域生態系統服務價值，討論了不同時期土地利用變化的生態影響；構建了土地利用轉換概率矩陣，應用馬爾科夫模型預測了未來40年的土地利用特征和生態系統服務價值變化趨勢，評價了不同時期流域土地利用和轉換格局的中長期生態影響，為流域生態-經濟-社會系統的長期協調發展提供決策參考。

1 材料與方法

1.1 研究區概況

達里諾爾湖流域位于內蒙古自治區赤峰市西北部克什克騰旗境內，地處大興安嶺山地、蒙古高原、渾善達克沙地三大地貌結合部，屬于森林草原交錯帶，包括森林、草原、沙地疏林和濕地等多種植被類型，流域面積4 963 km2。達里諾爾湖是流域中最大的水體，湖區面積228 km2。流域屬于溫帶半干旱大陸性季風氣候，冬季漫長寒冷降水少，夏季短促溫熱降水多，年降水量在300 mm左右。流域海拔高，大氣透明度強，日照時數長，蒸發旺盛，年均蒸發量為1 300～1 900 mm，蒸發量遠高于降水量。達里諾爾自然保護區始建于1987年，1997年升級為國家級自然保護區，主要保護對象為珍稀鳥類及其生境，是我國東部重要的鳥類繁殖地和候鳥遷徙驛站[25-27]。達里諾爾湖鹽堿度較高，水質獨特，主產瓦氏雅羅魚和鯽魚，經濟價值較高。進入21世紀后，受產卵條件惡化以及病害等因素影響，漁業產量整體下降[28]。

流域經濟總量較低，工業化和城市化水平不高。2000年之前，農牧業為流域主導產業，之后礦業發展壯大，2005年后旅游業日漸興起，近年來，隨著保護區生態保護政策力度加大，畜牧業、礦業和旅游業的發展都受到一定限制。

1.2 數據來源

該研究基于達里諾爾湖流域1980、1990、2000、2010、2019年共5期土地利用類型及轉換數據，數據由當地自然資源局提供。

1.3 研究方法

生態系統服務價值計算選擇可直接基于面積進行快速評估的當量因子法，主要依據謝高地等[29]1245修正Constanza的評估參數。達里諾爾湖流域土地利用類型多樣，一級生態系統類型包括農地、林地、草地、濕地、水域、不透水面和裸地(見表1)。不同土地利用類型所提供的生態系統服務價值有差異，根據謝高地等[29]1245構建的評價體系，結合達里諾爾湖流域半農半牧、子生態系統多樣化的實際情況，二級生態系統分類構成取均值，計算單位面積生態系統服務價值(以當量計，見表1)。由表1可見，水域和濕地單位面積的生態服務價值最高，其次是林地和草地，再次是農地，不透水面和裸地單位面積的生態服務價值最低。

表1 達里諾爾湖流域各土地利用類型單位面積生態系統服務價值Table 1 Ecological system service value per unit area of each landuse type of Dalinuoer Lake Basin

土地利用類型預測應用馬爾科夫模型，利用狀態之間的轉換概率矩陣來分析狀態的變化趨勢和預測其未來的狀態分布。馬爾科夫過程是一種特殊的隨機運動過程，每次狀態的轉換都只與其上一刻的狀態有關，而與過去的狀態無關，土地利用類型的轉換和動態變化具有類似馬爾科夫過程的特征[30]。本研究中，狀態指土地利用類型分布，狀態轉換指從一種土地利用類型轉換為另一種類型。馬爾科夫模型的關鍵是確定狀態轉換概率矩陣。設共有n個土地利用類型，狀態概率P(t)表示t時刻土地利用類型i的占比，則有：

(1)

式中：Sit表示t時刻土地利用類型i的面積，km2；St表示t時刻土地總面積，km2。

根據馬爾科夫過程的特征，事物在時刻t的狀態概率P(t)僅取決于其上一時刻的狀態概率P(t-1)和轉換概率Pij，即P(t)=P(t-1)×Pij，以此預測t時刻的土地利用類型。

2 結果與討論

2.1 土地利用與生態系統服務價值變化

1980—2019年，達里諾爾湖流域土地利用類型變化的總體趨勢利于生態保護。盡管有波動，但總體來看，草地和濕地占比有所上升，農地占比明顯下降，裸地顯著減少，林地、水域和不透水面略有縮減(見表2)。即生態系統服務價值較低的土地利用類型占比都有下降，價值較高的濕地、草地等占比則上升或只是略減。

表2 達里諾爾湖流域土地利用類型占比Table 2 Proportion of landuse type in Dalinuoer Lake Basin %

從生態系統服務價值的總體變化趨勢來看(見表3)，大致可分為2個階段：1980—2010年為增長階段，并在2010年達到峰值；2010年之后為減少階段，但仍基本高于2000年以前。各一級服務分類價值的峰值均出現于2010年，到2019年，供給服務價值大致降至低于2000年，調節、支持和文化服務價值大體高于2000年。從二級服務類型看，食物生產服務價值峰值出現于2000年，氣候調節價值從1980年至今持續增加；其余二級服務類型的價值峰值均出現于2010年。

表3 達里諾爾湖流域單位面積生態系統服務價值Table 3 Ecosystem service value per unit area of Dalinuoer Lake Basin

2.2 土地利用與生態系統服務價值預測

土地利用變化是一個長期和持續的過程?；?980—1990年(記為1980s)、1991—2000年(記為1990s)、2001—2010年(記為2000s)和2011—2019年(記為2010s)共4期土地利用類型轉換矩陣，假設土地利用類型轉換在2019—2020年仍保持2011—2019年的平均轉換概率，以此推算2020年的土地利用類型和2011—2020年的轉換概率，從而將4個時期的轉換步長統一為10年。分別根據4個時期的轉換矩陣構建馬爾科夫模型，預測了未來40年的土地利用類型構成情況和單位面積生態系統服務價值(見表4)。

由表4可見，不同時期土地利用格局的長期影響具有一致性?；?000s土地利用格局的未來單位面積生態系統服務價值最高，其次是1990s，最低的是1980s，近期(2010s)土地利用格局的長期生態效應不是最優狀態。對具體土地利用類型的影響則存在差異性，如：基于1980s土地利用格局預測的長期生態效應表現為單位面積生態服務價值較高的林地占比預測值略高于其他預測值，而裸地占比預測值顯著高于其他預測值；基于1990s的預測表現為農地占比顯著高，單位面積生態服務價值最高的水域占比略高；基于2000s的預測表現為濕地面積占比相對較高；基于2010s的預測則表現為草地占比較高。

表4 達里諾爾湖流域土地利用類型占比和單位面積生態系統服務價值預測Table 4 Prediction of landuse type composition and ecosystem service value per unit area in Dalinuoer Lake Basin

根據不同時期土地利用類型轉換數據預測未來40年達里諾爾湖流域生態系統服務價值變化趨勢(見表5)。從未來10年看，按2000s的土地利用類型轉換格局預測得到的總價值和各子類價值都是最高的，以最近期即2010s為基期的預測總價值甚至低于1990s為基期的預測總價值，但在支持和文化服務方面僅次于2000s為基期的預測值。從中長期(2040—2060年)來看，不同基期預測的價值靜態格局基本穩定，以2000s為基期的中長期總價值、支持和文化服務價值預測值均為最高，1990s為基期的供給和調節服務價值預測值最高，1980s為基期的各項預測值均為最低。以最近期2010s為基期的各項中長期預測值均未達到最優，支持服務和文化服務價值僅次于2000s為基期的預測值，而總價值、供給服務和調節服務價值僅高于1980s為基期的預測值，大致只相當于2019年的實際水平。

表5 達里諾爾湖流域單位面積生態系統各類服務價值預測Table 5 Prediction of various service value of ecosystem per unit area in Dalinuoer Lake Basin

不同時期土地利用類型轉換的生態系統服務價值變化的中長期趨勢有差異。因為預測時期較長，1980s土地利用類型轉換在2030—2060年的長期生態影響已經基本趨于穩定，生態系統服務價值穩定在17左右，大致相當于1990年的實際水平。1990s土地利用類型轉換對應的價值持續增加，在2040年接近2010年的實際水平，之后10年超過2010年的實際水平并繼續增長。2000s土地利用類型轉換的中長期價值表現為先增后降，是4期預測中最優的結果。但近期的2010s土地利用類型轉換所對應的價值只在2020—2030年有微小增長，之后均表現為持續下降趨勢。

2.3 討 論

對過去40年和未來40年生態系統服務價值的靜態比較和動態預測表明，達里諾爾湖流域在2010年左右達到土地利用的相對最優格局，反映了土地管理政策的成效；而近十年來表現出的生態系統服務價值減少趨勢，則是流域未來面臨的挑戰。

土地利用格局是一定時期自然生態和經濟社會要素綜合作用的結果。2010年達里諾爾湖生態系統服務價值的相對最優狀態，集中反映了流域在土地利用類型調整、生態環境管理方面所做出的長期努力，這些努力一方面體現在對不利自然條件的適應上，另一方面體現為對日益加大的人類活動干擾的不斷調控。從自然條件看，達里諾爾湖流域屬于內陸湖，主要水源補給來源于自然降水和入湖河流，近45年來，氣候暖干化進一步減少了自然降水和入湖水量，是水域面積萎縮的主要原因[31]。從人類活動干擾看，流域內牧業、農業、漁業、礦業、旅游業等多種自然資源依賴型產業并存，這些產業在土地利用上存在一定的沖突。牧業是當地傳統產業，主要依賴于天然草地，隨著20世紀90年代“草畜雙承包”體制改革的實施，集體草場和牲畜已全部承包到戶，雖然起到了落實產權和激勵生產的作用，但客觀上也造成了土地分割和破碎化，對草場的利用強度加大，草畜矛盾突出。流域農民多為外來人口，是在20世紀80年代之前的數次草原開墾行動中進來，在承載力較低的草原上，保證農業種植和養畜效率的方式只能是不斷粗放性地蠶食草場。

對自然因素造成的湖面萎縮等威脅，克什克騰旗政府主要采取了以加強流域生態功能為目的的各類工程建設措施，如人工造林、封山育林、工程固沙等，保持土地利用格局在不利的自然條件下基本不變或有所改善。對人類活動的擾動，則主要采取各類政策措施促使相關行為主體減少活動范圍、降低活動強度，使土地利用類型朝向更有利于提高生態服務功能的方向轉換。特別是在2000年之后，克什克騰旗政府先后發布了5個禁牧令，劃定了多個常年禁牧區，使達里諾爾湖流域土地利用結構調整在2000s發生了相對最大的變化，也使該時期生態系統服務價值達到最優。這些努力在近十年也仍在持續，但相較于之前的大型工程舉措，對土地利用格局的影響相對較小。與此同時，自然要素的影響力在2010s卻凸顯了出來，一定程度上也導致近十年來的人為努力難以得到體現。濕地和水域兩類單位面積生態系統服務價值最高的土地利用類型的減少不僅導致近十年流域生態系統服務總價值的下降，也意味著未來40年的價值難以達到最優水平。在未來的土地和生態管理中，達里諾爾湖流域需要直面這些威脅，在保持農地和草地面積以維持當地居民生計基本不受影響的同時，考慮通過各類人工工程或措施減少自然因素所造成的水域和濕地衰減，采取合適方式利用礦權退出以及各種原因未利用的土地，或作為經濟發展用地緩解生計困難，或作為生態用地提高其生態系統服務價值，以此加大土地利用格局向更高生態系統服務價值方向轉變的力度，才有可能重新回到最優狀態。

3 結 語

(1) 從土地利用格局與生態系統服務價值之間的關系看，與1980年相比，2019年達里諾爾湖土地利用格局的總體變化結果有利于生態保護，主要表現為單位面積生態服務價值相對較高的土地利用類型占比均有上升或只略有縮減，而單位面積生態服務價值相對較低的土地利用類型占比有所下降。但這個總體有利的變化結果包括了1980—2010年的改善階段和2011—2019年的下降階段，即2010年之前生態系統服務價值持續上升，在2010年達到相對最優，目前的土地利用狀態只是近40年來的次優水平。

(2) 依據土地利用類型之間的轉換對未來土地利用以及生態系統服務價值的預測結果表明，短期看，未來10年，以2000s為基期的生態系統服務價值總量和各子類價值預測值均為最高，而以近期2010s為基期的預測值均非最優；中長期看，未來30年，不同基期的預測結果較為一致，表現為2000s的土地利用類型轉換的中長期生態服務價值總量、支持和文化服務價值均為最高，相對較優，而以近期2010s為基期的各項中長期預測值均未達到最優。