張元，張福平，丁家寶，寧亞洲，曾攀儒

（陜西師范大學(xué)地理科學(xué)與旅游學(xué)院，陜西 西安 710119）

生態(tài)系統(tǒng)服務(wù)是指生態(tài)系統(tǒng)所形成及所維持的人類賴以生存的自然環(huán)境條件與效用，是人類直接或間接從生態(tài)系統(tǒng)中所獲得的所有收益［1，2］。人類通過一些生活生產(chǎn)行為，從生態(tài)系統(tǒng)中獲取對自身有利的福祉，這就使得生態(tài)系統(tǒng)會發(fā)生變化，從而產(chǎn)生權(quán)衡或協(xié)同的結(jié)果［3］。協(xié)同是指兩種或兩種以上生態(tài)系統(tǒng)服務(wù)之間的關(guān)系為正相關(guān)，即服務(wù)之間呈現(xiàn)出共同增加或共同減少的情況［4，5］。權(quán)衡一般是指兩種或兩種以上生態(tài)系統(tǒng)服務(wù)之間的關(guān)系為負相關(guān)，即一種服務(wù)的增加／減少會導(dǎo)致另一種服務(wù)的減少／增加［6］；但從廣義上講，權(quán)衡還可以是兩種或多種生態(tài)系統(tǒng)服務(wù)變化速率不均勻的單向變化［7］。研究表明，生態(tài)系統(tǒng)服務(wù)之間的權(quán)衡與協(xié)同是同時存在的，而且往往權(quán)衡作用要遠大于協(xié)同作用［8］，人類在進行選擇時，往往會傾向于選擇可以獲取利益的服務(wù)，而忽視一些調(diào)節(jié)、支持服務(wù)，這樣就會對生態(tài)平衡造成影響，甚至對生態(tài)系統(tǒng)造成破壞。所以研究生態(tài)系統(tǒng)服務(wù)的權(quán)衡協(xié)同關(guān)系，對于維持生態(tài)平衡、保護生態(tài)環(huán)境具有重要意義。

近年來，權(quán)衡協(xié)同的研究已成為該研究領(lǐng)域的熱點之一，且大多數(shù)都集中于研究生態(tài)系統(tǒng)服務(wù)權(quán)衡協(xié)同關(guān)系的時空特征。如王鵬濤等［9］研究了漢江上游生態(tài)系統(tǒng)服務(wù)權(quán)衡與協(xié)同關(guān)系時空特征，Wang等［10］分析了重點脆弱生態(tài)區(qū)生態(tài)系統(tǒng)服務(wù)權(quán)衡與協(xié)同關(guān)系時空特征，Pan等［11］研究了西北干旱內(nèi)陸河流域多種生態(tài)系統(tǒng)服務(wù)之間權(quán)衡和協(xié)同關(guān)系的時空格局。還有一些研究揭示了不同服務(wù)關(guān)系的尺度依賴，如余玉洋等［12］通過多種空間尺度對秦巴山區(qū)多種生態(tài)系統(tǒng)服務(wù)之間的關(guān)系進行探究，發(fā)現(xiàn)從500 m到200 km不同的空間尺度下，各種生態(tài)系統(tǒng)服務(wù)之間的關(guān)系會從權(quán)衡變?yōu)閰f(xié)同或者從協(xié)同變?yōu)闄?quán)衡；Qiao等［13］發(fā)現(xiàn)在全域尺度和縣域尺度下太湖流域生態(tài)系統(tǒng)服務(wù)的關(guān)系會有變化。但分析影響生態(tài)系統(tǒng)服務(wù)關(guān)系的研究較少，尤其是在干旱內(nèi)陸河流域。

黑河上游位于祁連山區(qū)，植被覆蓋良好，是重要的水土保持區(qū)。但近些年來由于超載放牧和毀林毀草，該流域植被減少，水土流失加劇，生態(tài)環(huán)境惡化［14］。2000年以后，國家啟動黑河專項治理規(guī)劃工程，以解決上游的生態(tài)問題。然而氣候變化和一系列的治理措施對生態(tài)系統(tǒng)服務(wù)及其關(guān)系也產(chǎn)生了重要影響。已有研究表明，降水和土地利用對于生態(tài)系統(tǒng)服務(wù)及其關(guān)系的影響最為顯著［15，16］。本研究以黑河上游為對象，利用InVEST模型和CASA模型測算了2000—2015年黑河上游的產(chǎn)水、土壤保持和植被凈初級生產(chǎn)力，進一步分析了降水和土地利用對3種服務(wù)間關(guān)系的影響，以期為黑河上游的生態(tài)恢復(fù)提供理論指導(dǎo)。

1 數(shù)據(jù)來源與研究方法

1.1 研究區(qū)概況

黑河發(fā)源于祁連山中段，是我國西北干旱區(qū)第二大內(nèi)陸河。西起黑山、鄰接疏勒河流域，東至大黃山、武威盆地，北部與蒙古接壤，流經(jīng)青海、甘肅、內(nèi)蒙古三?。ㄗ灾螀^(qū)）。黑河上游包含青海省祁連縣和甘肅省肅南裕固族自治縣的大部分地區(qū)，介于37°43′18＂～39°9′8＂N、99°9′53＂～101°0′16＂E（圖1）。研究區(qū)地勢南高北低，平均海拔3 581 m，氣候濕潤寒冷。河道兩岸山高谷深，植被覆蓋良好，主要植被類型為高山草甸［17］。黑河流域上游垂直地帶性分布明顯，土壤類型眾多［18］。

圖1 研究區(qū)概況

1.2 數(shù)據(jù)來源

InVEST模型和CASA模型所需數(shù)據(jù)如表1所示。

表1 數(shù)據(jù)來源與參數(shù)處理說明

1.3 研究方法

1.3.1 產(chǎn)水（water yield，WY）估算方法 利用InVEST產(chǎn)水模塊計算產(chǎn)水量。該模塊是基于Budyko水熱耦合平衡假設(shè)［20］，利用年平均降水量數(shù)據(jù)和蒸散數(shù)據(jù)計算，每個柵格單元的產(chǎn)水量為降水量減去實際蒸散后的水量，模型不做地表水、地下水、基流的區(qū)分。研究區(qū)每個柵格單元的年產(chǎn)水量，計算公式如下：

式中，WYxj為柵格單元x中土地利用類型為j的產(chǎn)水量（mm）；AETxj為柵格單元x中土地利用類型為j的實際蒸散量（mm）；Px為柵格單元x的降水量（mm）［21，22］。

1.3.2 土壤保持（soil conversation，SC）估算方法

InVEST土壤保持模塊認為土壤侵蝕減少量與持留量的差值即為土壤保持量，前者為土壤潛在侵蝕量與實際侵蝕量的差，后者為地塊攔截上坡沉積物的能力，以來沙量與泥沙持留率的乘積表示［23］。由此可見，該模型相對于通用土壤流失方程（USLE），考慮了植被在緩減侵蝕和截留上坡來沙的作用［24，25］。計算公式如下：

式中，SCx、PKLSx、USLEx和SEDRx分別為柵格x的土壤保持量（t）、潛在侵蝕量（t）、實際侵蝕量（t）和持留量（t）。

1.3.3 植被凈初級生產(chǎn)力（net primary productivity，NPP）估算方法 利用CASA模型計算研究區(qū)的植被凈初級生產(chǎn)力。在模型中，光合有效輻射與光能利用率的乘積代表植被NPP的值［26］。公式如下：

式中，NPP（x，t）為像元x在t時間內(nèi)單位面積植被凈初級生產(chǎn)力（gC·m－2），APAR（x，t）代表像元x在t時間內(nèi)吸收的光合有效輻射（MJ·m－2），ε（x，t）表示像元x在t時間內(nèi)植被實際的光能利用率（gC·MJ）。

1.3.4 生態(tài)系統(tǒng)服務(wù)間相關(guān)性分析方法 利用Matlab軟件對黑河上游2000—2015年間3種生態(tài)系統(tǒng)服務(wù)間的相關(guān)性進行量化分析。利用t檢驗，將生態(tài)系統(tǒng)服務(wù)的權(quán)衡與協(xié)同關(guān)系分為7類，即協(xié)同**（r＞0，t≤0.05），協(xié)同*（r＞0，0.05＜t≤0.1），協(xié)同（r＞0，t＞0.1），權(quán)衡**（r＜0，t≤0.05），權(quán)衡*（r＜0，0.05＜t≤0.1），權(quán)衡（r＜0，t＞0.1）；其余的為不相關(guān)［27］。

1.3.5 標(biāo)準化處理 由于服務(wù)的單位和范圍值不一致，為了便于比較，對服務(wù)數(shù)值進行標(biāo)準化處理，公式如下，其中ESi為服務(wù)量化任意值，ESmin指服務(wù)量化最小值，ESmax指服務(wù)量化最大值［28］。

2 結(jié)果與分析

2.1 黑河上游生態(tài)系統(tǒng)服務(wù)權(quán)衡與協(xié)同關(guān)系的時空特征

2.1.1 黑河上游生態(tài)系統(tǒng)服務(wù)及其關(guān)系的時間特征分析 黑河上游流域2000—2015年產(chǎn)水、土壤保持及NPP三種服務(wù)均值的年際變化情況如圖2a、b、c所示。可見，黑河上游地區(qū)產(chǎn)水、土壤保持和NPP年際變化波動較大，總體呈增加趨勢。2000—2015年產(chǎn)水量平均值為60.07 mm，最低值在2002年（47.38 mm），最高值在2007年（73.94 mm），以0.64 mm·a－1的速率遞增；產(chǎn)水總量年際波動極大，尤其是2002—2007年間，產(chǎn)水低值、高值交替出現(xiàn)。土壤保持的多年平均值為435.52 t·hm－2，2000—2012年間土壤保持量在波動中增加，并于2012年達到峰值（692.82 t·hm－2），之后每年的土壤保持量均下降，至2015年出現(xiàn)急劇降低的情況。研究區(qū)內(nèi)NPP多年平均總量為277.39 gC·m－2，在238.84～315.84 gC·m－2之間上下浮動。

對3種服務(wù)進行標(biāo)準化處理，利用多項式擬合兩兩服務(wù)之間隨時間的變化趨勢，由圖2d、e、f可見，產(chǎn)水與土壤保持在2000—2010年變化趨勢一致，而在2010年后呈現(xiàn)出明顯的產(chǎn)水增加土壤保持減少／產(chǎn)水減少土壤保持增加的變化規(guī)律，這表明產(chǎn)水與土壤保持在2000—2010年之間是協(xié)同關(guān)系，而在2011—2015年為權(quán)衡關(guān)系。NPP與產(chǎn)水在2000—2010為權(quán)衡關(guān)系，而在2011—2015年轉(zhuǎn)變?yōu)閰f(xié)同關(guān)系。NPP與土壤保持在研究時段內(nèi)一直是此消彼長的權(quán)衡關(guān)系。

圖2 生態(tài)系統(tǒng)服務(wù)及其關(guān)系的時間變化趨勢

2.1.2 黑河上游生態(tài)系統(tǒng)服務(wù)及其關(guān)系的空間特征分析 從空間分布情況來看，2000—2015年產(chǎn)水、土壤保持和NPP三種服務(wù)均呈現(xiàn)出東南高西北低、山地高河谷低的分布態(tài)勢（圖3a、b、c）。16年間，黑河上游流域產(chǎn)水深度平均值為77.62 mm，但空間分布不均，整體呈現(xiàn)由東南向西北遞減的分布規(guī)律，青海省境內(nèi)的祁連縣是產(chǎn)水分布的高值區(qū)，其北部靠近祁連山山麓地帶及與甘肅交界的邊緣地區(qū)產(chǎn)水量均高于129 mm；沿河谷分布著均值在59 mm的產(chǎn)水地帶；而低值區(qū)則出現(xiàn)在肅南西北角和鄰接高臺、臨澤、張掖的南部地區(qū)。土壤保持年均值為436.12 t·hm－2，空間上也是呈現(xiàn)出西北劣于東南的分布情況；但不同的是，土壤保持無特別明顯的高值集中區(qū)，只是在研究區(qū)東部、沿河谷地帶植被覆蓋較好的地區(qū)有呈條帶狀延伸的分布規(guī)律；在肅南與酒泉接壤的地域有零星分布的高值點，而低值區(qū)與高值區(qū)穿插分布。NPP的分布情況與產(chǎn)水基本一致。

從全區(qū)來看，黑河上游的大部分地區(qū)產(chǎn)水和土壤保持協(xié)同關(guān)系，且大部分區(qū)域協(xié)同關(guān)系的t檢驗值大于0.1，只是在一些河谷地帶會有權(quán)衡集中分布的情況（圖3d）。產(chǎn)水與NPP的相互關(guān)系以協(xié)同為主，協(xié)同的區(qū)域主要分布在沿河流地區(qū)；在研究區(qū)北部還出現(xiàn)了協(xié)同關(guān)系和權(quán)衡關(guān)系相間分布的情況（圖3e）。土壤保持與NPP主要呈現(xiàn)權(quán)衡關(guān)系，但這種相互關(guān)系呈現(xiàn)出顯著的空間差異：黑河上游東部地區(qū)NPP與土壤保持的相互關(guān)系以權(quán)衡關(guān)系為主，并且這種關(guān)系呈條帶狀一直延伸到肅南西部；而在肅南西南部沿河谷地帶、肅南鄰接高臺、臨澤、張掖的南部地區(qū)及部分邊界地區(qū)呈現(xiàn)出協(xié)同關(guān)系（圖3f）。

圖3 生態(tài)系統(tǒng)服務(wù)及其關(guān)系的空間分布

2.2 降水和土地利用對黑河上游生態(tài)系統(tǒng)服務(wù)權(quán)衡與協(xié)同關(guān)系的影響

2.2.1 降水對黑河上游生態(tài)系統(tǒng)服務(wù)權(quán)衡與協(xié)同關(guān)系的影響 基于上述分析，發(fā)現(xiàn)2010年左右服務(wù)之間的關(guān)系發(fā)生變化，因此，為了避免其他因素的影響，除降水?dāng)?shù)據(jù)外，其他數(shù)據(jù)均使用2010年的數(shù)據(jù)。經(jīng)計算得到不同降水下產(chǎn)水、土壤保持和NPP的標(biāo)準化值，利用多項式擬合其趨勢，結(jié)果如圖4。發(fā)現(xiàn)在降水小于320 mm時，產(chǎn)水和土壤保持的變化速率相差不大，尤其是在290～320 mm之間二者的變化速率幾乎一致；但當(dāng)降水達到320 mm之后，二者的變化速率相差越來越大，至370 mm時出現(xiàn)明顯的負相關(guān)，之后兩者的變化速率差距又縮小（圖4a）。產(chǎn)水與NPP的變化速率在降水小于320 mm時差異較大，但當(dāng)降水達到320 mm后，二者的變化速率開始接近（圖4b）。隨著降水的增加，土壤保持與NPP之間的變化速率差距處于增大－縮小的波動變化之中，但以增大為主（圖4c）。所以降水的增加使得產(chǎn)水與土壤保持由協(xié)同變?yōu)闄?quán)衡，產(chǎn)水與NPP由權(quán)衡變?yōu)閰f(xié)同；而土壤保持與NPP之間的權(quán)衡與協(xié)同關(guān)系不斷變化，總體以權(quán)衡為主。

圖4 降水對生態(tài)系統(tǒng)服務(wù)關(guān)系的影響

根據(jù)氣象站點多年降水?dāng)?shù)據(jù)統(tǒng)計分析可知，2000—2015年降水均值為311.15 mm，以2.20 mm·a－1的速率增長，而2010年的年均降水為324.38 mm，所以2010年左右服務(wù)關(guān)系發(fā)生變化的原因之一為降水量增加。降水增加會引起產(chǎn)水量增長［29］，同時也會增加降水侵蝕性，造成土壤保持量減少［30］，所以降水會成為改變產(chǎn)水和土壤保持關(guān)系的因素之一。在前期降水較少，植被會吸收大部分降水維持自身生長，整個系統(tǒng)的產(chǎn)水就會減少，這時NPP與產(chǎn)水呈現(xiàn)權(quán)衡關(guān)系；當(dāng)降水增加后，除完全滿足植被生長外，大部分降水轉(zhuǎn)化為產(chǎn)水，此時NPP與產(chǎn)水就會轉(zhuǎn)變?yōu)閰f(xié)同關(guān)系。

2.2.2 土地利用對黑河上游生態(tài)系統(tǒng)服務(wù)權(quán)衡與協(xié)同關(guān)系的影響 利用中國科學(xué)院資源環(huán)境科學(xué)數(shù)據(jù)中心（http：／／www.resdc.cn／）提供的2000、2005、2010、2015年四期土地利用數(shù)據(jù)，把研究區(qū)內(nèi)土地利用類型按面積從高到低排序，選取占研究區(qū)總面積95%以上的低覆蓋度草地、中覆蓋度草地、其它未利用地（包括高寒荒漠、苔原等）、高覆蓋度草地、裸巖石質(zhì)地、灌木林地、有林地及疏林地進行生態(tài)系統(tǒng)服務(wù)及其關(guān)系的統(tǒng)計分析，結(jié)果見表2。可見，疏林地的產(chǎn)水最多，低覆蓋度草地的產(chǎn)水最少；林地的覆蓋度與其產(chǎn)水能力呈反比，而草地的覆蓋度與其產(chǎn)水呈正比；土壤保持能力為林地＞草地＞未利用地，隨著地表植被覆蓋度的減少而減少。在NPP中，林地的平均能力要優(yōu)于草地，但出現(xiàn)了高覆蓋度草地略大于疏林地的特殊情況。

表2 不同土地類型各生態(tài)系統(tǒng)服務(wù)多年均值

降水是影響產(chǎn)水能力的主要因素［31］，林地可以有效攔截降雨，雖然土壤滲透性強、蒸散量高會導(dǎo)致林地產(chǎn)水量減少，但研究區(qū)林地多位于降水高值區(qū)，所以產(chǎn)水最高。與其他樹木相比，灌木耗水較少，所以灌木林地的產(chǎn)水也較高。未利用地對降水的消耗最少，降水大部分轉(zhuǎn)化為產(chǎn)水，所以產(chǎn)水也較高。林地有著強大的根系，而草地的根系較淺，所以草地攔截能力較弱，土壤保持低于林地；除高覆蓋度草地的NPP高于疏林地外，草地的NPP低于林地。

由表3可見，各種土地類型的產(chǎn)水和NPP均呈協(xié)同**關(guān)系；林地和高覆蓋度草地的土壤保持和NPP呈協(xié)同**關(guān)系，其余土地類型的土壤保持和NPP呈權(quán)衡**關(guān)系。裸巖石質(zhì)地和其它未利用地的產(chǎn)水與土壤保持為權(quán)衡*關(guān)系，這是因為裸巖石質(zhì)地和其它未利用地的產(chǎn)水較高，但土壤保持能力幾乎為零；林地和草地中產(chǎn)水和土壤保持為協(xié)同關(guān)系，雖然林地的產(chǎn)水和土壤保持能力最高，但產(chǎn)水與土壤保持的協(xié)同關(guān)系卻弱于草地；對于草地來說，覆蓋度越高產(chǎn)水越高，土壤保持能力越強，在中覆蓋度草地中二者呈現(xiàn)出最強的協(xié)同關(guān)系。所有土地利用類型中產(chǎn)水與NPP的關(guān)系均為協(xié)同**，但林地的這種協(xié)同關(guān)系要弱于草地，這是因為林地的耗水要大于草地。綜上所述，林地，尤其是灌木林地，土壤保持能力和NPP最為突出，但服務(wù)之間的協(xié)同發(fā)展，尤其是水供給方面的協(xié)同發(fā)展卻要弱于草地。產(chǎn)水與土壤保持、產(chǎn)水與NPP的協(xié)同關(guān)系均以中覆蓋度草地的最高，高覆蓋度草地的次之；而土壤保持與NPP的協(xié)同關(guān)系則以灌木林地的最大，高覆蓋度草地的次之。所以高覆蓋度草地的綜合協(xié)同能力最優(yōu)。

表3 不同土地類型生態(tài)系統(tǒng)服務(wù)間的相關(guān)性

3 討論

本研究分析了黑河上游生態(tài)系統(tǒng)服務(wù)及其關(guān)系的時空特征，進一步分析了降水和土地利用對生態(tài)系統(tǒng)服務(wù)權(quán)衡與協(xié)同關(guān)系的影響。結(jié)果表明，不同尺度的分析可能得到不同的結(jié)論，如：在時間和空間尺度上土壤保持與NPP均為權(quán)衡關(guān)系；但產(chǎn)水與土壤保持在時間尺度上表現(xiàn)為先協(xié)同后權(quán)衡的關(guān)系，產(chǎn)水與NPP表現(xiàn)為先權(quán)衡后協(xié)同的關(guān)系，而在空間尺度中他們之間均以協(xié)同關(guān)系為主。此外，研究所得到的“320 mm左右的降水會改變服務(wù)關(guān)系”的結(jié)論也只是針對特定條件下的研究區(qū)得出的（除降水?dāng)?shù)據(jù)外，其他數(shù)據(jù)均使用2010年的數(shù)據(jù)），而且也只是得出“降水增加會使產(chǎn)水與土壤保持由協(xié)同變?yōu)闄?quán)衡、產(chǎn)水與NPP由權(quán)衡變?yōu)閰f(xié)同”這樣定性的結(jié)論，至于降水值究竟為多少才能發(fā)生上述關(guān)系的轉(zhuǎn)變還需要進一步探究。

本研究結(jié)果表明降水會對生態(tài)系統(tǒng)服務(wù)之間的關(guān)系產(chǎn)生影響，而大量研究證實黑河上游的氣候呈暖濕化發(fā)展趨勢［32，33］，所以為了應(yīng)對降水增加帶來的影響，尤其是在土壤保持方面的影響，我們需要一些合理有效的措施，植樹造林、提高草地覆蓋度等措施都可以提高土壤保持量。但與此同時我們還需要注意的是：過量的植被會消耗降水使得區(qū)域的產(chǎn)水變少，尤其是林地，雖然水源涵養(yǎng)和土壤保持能力都高，但其蒸騰作用也強［34］，會消耗大量降水。

黑河上游作為黑河流域的主產(chǎn)流區(qū)，必須要適度發(fā)展才能為解決中下游用水矛盾提供有力保障。黑河上游草地面積最大，所以草地，尤其是面積占比高的中、低覆蓋度草地中各種服務(wù)的關(guān)系會直接影響到整個黑河流域。本研究結(jié)果表明，中覆蓋度草地中產(chǎn)水與土壤保持、產(chǎn)水與NPP均為最高協(xié)同，土壤保持與NPP為權(quán)衡關(guān)系；高覆蓋度草地的綜合協(xié)同能力最優(yōu)。因此，合理提高黑河上游的草地覆蓋度對于區(qū)域協(xié)調(diào)發(fā)展是十分必要的。

4 結(jié)論

本研究分析了黑河上游生態(tài)系統(tǒng)服務(wù)及其關(guān)系的時空特征，并進一步分析了降水和土地利用對生態(tài)系統(tǒng)服務(wù)權(quán)衡、協(xié)同關(guān)系的影響。結(jié)論如下：①2000—2015年產(chǎn)水、土壤保持、植被凈初生產(chǎn)力三種服務(wù)在波動中增加。產(chǎn)水與土壤保持在2000—2010年之間是協(xié)同關(guān)系，而在2011—2015年為權(quán)衡關(guān)系；NPP與產(chǎn)水在2000—2010為權(quán)衡關(guān)系，而在2011—2015年轉(zhuǎn)變?yōu)閰f(xié)同關(guān)系；NPP與土壤保持一直是此消彼長的權(quán)衡關(guān)系。從空間分布來看，產(chǎn)水、土壤保持和NPP三種服務(wù)均呈現(xiàn)出東南高西北低、山地高河谷低的分布態(tài)勢；產(chǎn)水與土壤保持、產(chǎn)水與NPP在空間上以協(xié)同關(guān)系為主，土壤保持與NPP則以權(quán)衡關(guān)系為主。②降水的增加會使產(chǎn)水與土壤保持由協(xié)同關(guān)系變?yōu)闄?quán)衡關(guān)系，而使產(chǎn)水與NPP由權(quán)衡關(guān)系變?yōu)閰f(xié)同關(guān)系；降水對土壤保持與NPP的關(guān)系無明顯影響。③中覆蓋度草地的產(chǎn)水與土壤保持、產(chǎn)水與NPP均為最高協(xié)同，高覆蓋度草地次之；灌木林地的土壤保持與NPP協(xié)同關(guān)系最大，高覆蓋度草地次之。因此，高覆蓋度草地的綜合協(xié)同能力最優(yōu)。