價(jià)值轉(zhuǎn)移模型在森林類保護(hù)區(qū)生態(tài)系統(tǒng)服務(wù)功能評(píng)估中的應(yīng)用1)

舒航 莊立文 孫曉杰 邢一明 桑衛(wèi)國(guó)

(中央民族大學(xué)，北京，100081) (黑龍江省林業(yè)科學(xué)研究院) (太原師范學(xué)院) (中央民族大學(xué)) (中國(guó)科學(xué)院植物研究所)

森林生態(tài)系統(tǒng)作為陸地生態(tài)系統(tǒng)的主體，在人類社會(huì)的經(jīng)濟(jì)發(fā)展、生態(tài)安全方面發(fā)揮了重要作用[1]。隨著人口的劇增、社會(huì)經(jīng)濟(jì)的飛速發(fā)展，人類對(duì)森林的需求量不斷加大；為保護(hù)大面積的森林生態(tài)系統(tǒng)，我國(guó)于1956年建立第一個(gè)自然保護(hù)區(qū)——鼎湖山自然保護(hù)區(qū)以來(lái)[2]，自然保護(hù)區(qū)事業(yè)得到飛速發(fā)展。國(guó)家級(jí)自然保護(hù)區(qū)是自然保護(hù)區(qū)中的精華，對(duì)保護(hù)我國(guó)的生態(tài)環(huán)境、自然資源具有舉足輕重的價(jià)值[3]。森林類型自然保護(hù)區(qū)，主要是指以森林植被及其生境所形成的生態(tài)系統(tǒng)為保護(hù)對(duì)象的自然保護(hù)區(qū)[4]。評(píng)估我國(guó)國(guó)家級(jí)森林類型自然保護(hù)區(qū)的生態(tài)系統(tǒng)服務(wù)功能價(jià)值，可以使民眾更加直觀地了解保護(hù)區(qū)的生態(tài)價(jià)值，對(duì)于維護(hù)區(qū)域生態(tài)安全、促進(jìn)經(jīng)濟(jì)社會(huì)可持續(xù)發(fā)展具有重要意義[5-6]。

在對(duì)我國(guó)國(guó)家級(jí)森林類型自然保護(hù)區(qū)進(jìn)行評(píng)估時(shí)，往往會(huì)受到評(píng)估時(shí)間、評(píng)估成本等因素的影響，難以在短時(shí)間內(nèi)獲取完整可靠的第一手資料。截至目前，國(guó)內(nèi)學(xué)者通過(guò)實(shí)地調(diào)研的方式，僅對(duì)我國(guó)31處此類型自然保護(hù)區(qū)進(jìn)行了評(píng)估，占此類型自然保護(hù)區(qū)總數(shù)的15.12%。此外，由于評(píng)估成本較高，對(duì)所有的保護(hù)區(qū)進(jìn)行實(shí)地評(píng)估也是不必要的[7]。因此，為更加科學(xué)有效地快速評(píng)估某一區(qū)域的生態(tài)系統(tǒng)服務(wù)價(jià)值，國(guó)內(nèi)外出現(xiàn)了大量的依據(jù)文獻(xiàn)整合的Meta分析(又叫薈萃分析或元分析)價(jià)值轉(zhuǎn)移法研究實(shí)例[8-9]。Meta分析價(jià)值轉(zhuǎn)移法，可以根據(jù)研究地的價(jià)值評(píng)估結(jié)果、評(píng)估方法、生態(tài)系統(tǒng)服務(wù)類型、自然和社會(huì)經(jīng)濟(jì)情況等建立相關(guān)的價(jià)值轉(zhuǎn)移函數(shù)，再將其轉(zhuǎn)移到政策地(待研究地)，從而得到政策地的生態(tài)系統(tǒng)服務(wù)功能價(jià)值[10-11]。目前我國(guó)在生態(tài)系統(tǒng)服務(wù)價(jià)值轉(zhuǎn)移方面的應(yīng)用，多集中于大面積湖沼濕地、土地利用、森林生態(tài)系統(tǒng)方面[12-13]，在自然保護(hù)區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)估的應(yīng)用研究較少。為此，本研究通過(guò)收集與中國(guó)國(guó)家級(jí)森林類型自然保護(hù)區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)估相關(guān)的文獻(xiàn)，建立價(jià)值轉(zhuǎn)移數(shù)據(jù)庫(kù)、Meta分析價(jià)值轉(zhuǎn)移模型，分析價(jià)值評(píng)估的影響因素，檢驗(yàn)?zāi)Ｐ蜆颖就鈨r(jià)值轉(zhuǎn)移的有效性，探討Meta價(jià)值轉(zhuǎn)移方法在中國(guó)國(guó)家級(jí)森林類型自然保護(hù)區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)估的應(yīng)用，旨在為快速高效的評(píng)估此類型自然保護(hù)區(qū)的生態(tài)系統(tǒng)服務(wù)價(jià)值提供參考。

1 研究方法

1.1 文獻(xiàn)數(shù)據(jù)的收集

中文文獻(xiàn)數(shù)據(jù)來(lái)源于中國(guó)學(xué)術(shù)期刊網(wǎng)絡(luò)出版總庫(kù)(CNKI)，以主題詞、關(guān)鍵詞等方式，檢索1990—2018年發(fā)表的有關(guān)中國(guó)國(guó)家級(jí)森林類型自然保護(hù)區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)估的文獻(xiàn)。具有以下3種情況的文獻(xiàn)將被排除——①評(píng)估對(duì)象為非國(guó)家級(jí)森林類型自然保護(hù)區(qū)的文獻(xiàn)，盡管有一些以野生動(dòng)植物(如大熊貓、珙桐)為主要保護(hù)對(duì)象的保護(hù)區(qū)，也是以森林生態(tài)系統(tǒng)為主的，但是在劃分的時(shí)候，并沒(méi)有將其定為森林型的，因此，這些保護(hù)區(qū)在本研究中也不予考慮；②使用能值分析方法、價(jià)值轉(zhuǎn)移方法的文獻(xiàn)；③相同研究區(qū)域被多次評(píng)估的，對(duì)其進(jìn)行評(píng)估和篩選，最終確定入選文獻(xiàn)。本研究最終使用的數(shù)據(jù)包含31篇文獻(xiàn)，其中20篇是期刊文獻(xiàn)、11篇是碩博學(xué)位論文。根據(jù)中國(guó)植被分區(qū)、統(tǒng)計(jì)研究文獻(xiàn)中的樣本觀測(cè)點(diǎn)所處地理區(qū)位，被評(píng)估的保護(hù)區(qū)主要分布在亞熱帶常綠闊葉林區(qū)域。將入選的數(shù)據(jù)庫(kù)文獻(xiàn)信息統(tǒng)一錄入到Excel表格中，建立Meta分析價(jià)值轉(zhuǎn)移數(shù)據(jù)庫(kù)；對(duì)樣本文獻(xiàn)中的數(shù)據(jù)采用統(tǒng)一的方式進(jìn)行提取，主要包括文獻(xiàn)名、文獻(xiàn)作者、文獻(xiàn)來(lái)源、研究時(shí)間、研究區(qū)域、自然保護(hù)區(qū)名稱、自然保護(hù)區(qū)面積、生態(tài)系統(tǒng)服務(wù)類型、生態(tài)系統(tǒng)服務(wù)價(jià)值及評(píng)估方法等信息。

依據(jù)聯(lián)合國(guó)千年生態(tài)系統(tǒng)評(píng)估的分類，對(duì)國(guó)家級(jí)森林類型自然保護(hù)區(qū)的生態(tài)服務(wù)功能進(jìn)行劃分，主要分為供給服務(wù)、調(diào)節(jié)服務(wù)、支持服務(wù)、文化服務(wù)四大類[14]。在樣本文獻(xiàn)中，國(guó)家級(jí)森林類型自然保護(hù)區(qū)的供給服務(wù)通常包括林木產(chǎn)品，調(diào)節(jié)服務(wù)包括涵養(yǎng)水源、固碳釋氧、凈化空氣，支持服務(wù)包括營(yíng)養(yǎng)積累、土壤保持、生物多樣性保護(hù)，文化服務(wù)包括森林休憩、科研教育。在樣本文獻(xiàn)中，科研教育的價(jià)值評(píng)估多采用的是價(jià)值轉(zhuǎn)移方法，因此本文暫不將此項(xiàng)生態(tài)系統(tǒng)服務(wù)價(jià)值納入研究范圍。最終確定本研究評(píng)估的我國(guó)國(guó)家級(jí)森林類型自然保護(hù)區(qū)生態(tài)系統(tǒng)服務(wù)的類型，主要包括林果產(chǎn)品、涵養(yǎng)水源、營(yíng)養(yǎng)累積、固碳釋氧、凈化空氣、土壤保持、保護(hù)生物多樣性、休閑游憩。

在生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)估中，不同的生態(tài)系統(tǒng)服務(wù)類型往往會(huì)采取不同的評(píng)估方法，而且同一類型的生態(tài)系統(tǒng)服務(wù)所采取的評(píng)估方法也是不盡相同的[15-16]。樣本文獻(xiàn)中使用的價(jià)值評(píng)估方法，主要有市場(chǎng)價(jià)值法、機(jī)會(huì)成本法、影子工程法、碳稅法、造林成本法、污染成本防治法、旅行費(fèi)用法、支付意愿法等。對(duì)某一生態(tài)系統(tǒng)服務(wù)價(jià)值類型的評(píng)估往往不僅限于一種評(píng)估方法，可能綜合多種評(píng)估方法，如營(yíng)養(yǎng)物質(zhì)積累價(jià)值的評(píng)估，既用到了影子工程法，也采用了機(jī)會(huì)成本法。有的價(jià)值評(píng)估方法，由于受經(jīng)費(fèi)限制、數(shù)據(jù)可得性不高等原因，在實(shí)際應(yīng)用時(shí)存在很多限制，使用頻率較低。因此，將使用頻率較高的評(píng)估方法納入本研究范圍(見(jiàn)表1)。

由于樣本文獻(xiàn)中的價(jià)值觀察值(以下統(tǒng)稱為價(jià)值觀察值)的評(píng)估基準(zhǔn)年通常不在同一年，數(shù)據(jù)之間不具有可比性，從而需要將不同年份的價(jià)值觀察值調(diào)整到統(tǒng)一年份[17]。本研究在價(jià)值觀察值調(diào)整時(shí)，首先找到文獻(xiàn)引用的各個(gè)價(jià)格參數(shù)的原始文獻(xiàn)及其相應(yīng)的評(píng)估基準(zhǔn)年，并以此評(píng)估基準(zhǔn)年為基礎(chǔ)，借鑒國(guó)外研究中常用的價(jià)值觀察值的折算方法[18]，使用消費(fèi)者物價(jià)指數(shù)(CPI)將不同評(píng)估基準(zhǔn)年的價(jià)值觀察值統(tǒng)一調(diào)整到2015年的物價(jià)水平。

1.2 Meta分析價(jià)值轉(zhuǎn)移模型

Meta分析價(jià)值轉(zhuǎn)移方法在生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)估領(lǐng)域中的應(yīng)用，主要是指在已有的生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)估框架下，建立一個(gè)具有普遍適用性的Meta分析價(jià)值轉(zhuǎn)移模型，該模型可對(duì)同類型區(qū)域的生態(tài)系統(tǒng)服務(wù)價(jià)值進(jìn)行預(yù)測(cè)。國(guó)內(nèi)外在構(gòu)建生態(tài)系統(tǒng)服務(wù)價(jià)值的Meta回歸模型方面，已有很多實(shí)例和經(jīng)驗(yàn)[19]。在Meta回歸模型構(gòu)建中，為減少原始數(shù)據(jù)波動(dòng)程度和非對(duì)稱性，提高模型的擬合精度，減小異方差，便于結(jié)果之間的比較，一般采用對(duì)數(shù)線性模型進(jìn)行構(gòu)建[20]，價(jià)值觀察值采用自然對(duì)數(shù)形式。Meta回歸模型的一般形式為：

ln(yi)=a+bmXm+bpXp+beXe+c。

(1)

式中：yi為調(diào)整后以2015年為基期的國(guó)家級(jí)森林類型自然保護(hù)區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值的樣本觀察值；a為常數(shù)項(xiàng)；c為殘差項(xiàng)；bm、bp、be為影響因素(自變量)的系數(shù)向量；Xm、Xp、Xe為影響因素的特征向量，Xm為評(píng)估方法向量、Xp為國(guó)家級(jí)森林類型自然保護(hù)區(qū)的特征向量、Xe為社會(huì)經(jīng)濟(jì)特征向量。

1.3 Meta分析價(jià)值轉(zhuǎn)移模型中影響因素的選取與賦值

評(píng)估方法：不同的評(píng)估方法使用，會(huì)導(dǎo)致生態(tài)系統(tǒng)服務(wù)功能價(jià)值的評(píng)估結(jié)果有所不同，對(duì)價(jià)值評(píng)估結(jié)果產(chǎn)生影響。因此，本研究將價(jià)值評(píng)估方法作為影響因素，分析不同評(píng)估方法的使用對(duì)國(guó)家級(jí)森林類型自然保護(hù)區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值變化的影響。

生態(tài)系統(tǒng)服務(wù)類型：生態(tài)系統(tǒng)可提供多種類型的生態(tài)系統(tǒng)服務(wù)，并產(chǎn)生相應(yīng)的生態(tài)系統(tǒng)服務(wù)功能價(jià)值，有直接經(jīng)濟(jì)價(jià)值，也有間接經(jīng)濟(jì)價(jià)值[21]。不同生態(tài)系統(tǒng)服務(wù)類型的重要性不同，所產(chǎn)生的經(jīng)濟(jì)價(jià)值也存在差異，對(duì)生態(tài)系統(tǒng)服務(wù)的總價(jià)值也會(huì)產(chǎn)生相應(yīng)影響。

所處的植被分區(qū)：國(guó)家級(jí)森林類型自然保護(hù)區(qū)所處的植被分區(qū)，可以反映出保護(hù)區(qū)內(nèi)的植被生長(zhǎng)環(huán)境以及植被生長(zhǎng)狀況，森林生態(tài)系統(tǒng)的狀況會(huì)直接影響保護(hù)區(qū)的生態(tài)系統(tǒng)服務(wù)功能價(jià)值[22]。

國(guó)家級(jí)森林類型自然保護(hù)區(qū)成立的時(shí)間：建立自然保護(hù)區(qū)對(duì)于保護(hù)區(qū)域內(nèi)生態(tài)環(huán)境、提高區(qū)域內(nèi)生態(tài)系統(tǒng)服務(wù)價(jià)值具有重要意義[23-24]。而國(guó)家級(jí)自然保護(hù)區(qū)作為自然保護(hù)區(qū)中的精華，其成立必然會(huì)對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值產(chǎn)生一定的影響。

受益人口數(shù)量：生態(tài)系統(tǒng)服務(wù)的受益人口數(shù)量，反映了區(qū)域?qū)?guó)家級(jí)森林類型自然保護(hù)區(qū)生態(tài)系統(tǒng)服務(wù)需求的大小。將受益人口數(shù)量界定為保護(hù)區(qū)所處的地級(jí)行政區(qū)的常住人口，人口數(shù)據(jù)來(lái)源于各省歷年《統(tǒng)計(jì)年鑒》[25]。

保護(hù)區(qū)的面積：生態(tài)系統(tǒng)存在臨界閾值，其規(guī)模必須大于一定的閾值才可以提供生態(tài)系統(tǒng)服務(wù)，并且其生態(tài)系統(tǒng)服務(wù)功能價(jià)值會(huì)隨著面積的增大而增大。但隨著保護(hù)區(qū)面積的增加，某些生態(tài)系統(tǒng)服務(wù)的單位面積價(jià)值會(huì)存在規(guī)模收益遞減現(xiàn)象。

區(qū)域經(jīng)濟(jì)發(fā)展?fàn)顩r：國(guó)家級(jí)自然保護(hù)區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值的大小與保護(hù)區(qū)所處區(qū)域的經(jīng)濟(jì)發(fā)展水平有很大關(guān)系，發(fā)達(dá)地區(qū)和欠發(fā)達(dá)地區(qū)對(duì)保護(hù)區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值的認(rèn)定和市場(chǎng)實(shí)現(xiàn)存在差距[26]。本研究選取保護(hù)區(qū)所在地級(jí)市的人均GDP作為衡量區(qū)域經(jīng)濟(jì)發(fā)展?fàn)顩r的指標(biāo)，人均GDP數(shù)據(jù)來(lái)源于各省歷年的《統(tǒng)計(jì)年鑒》。

根據(jù)收集到的文獻(xiàn)，按照統(tǒng)計(jì)和計(jì)量的數(shù)據(jù)要求，對(duì)影響因素的各類信息進(jìn)行“編碼”和賦值，并得到各影響因素的均值、標(biāo)準(zhǔn)差、樣本數(shù)量(見(jiàn)表2)。

表2 Meta回歸模型的影響因素信息

1.4 Meta分析價(jià)值轉(zhuǎn)移模型的有效性檢驗(yàn)

在使用Meta分析價(jià)值轉(zhuǎn)移模型進(jìn)行效益轉(zhuǎn)移時(shí)，需要對(duì)模型的有效性進(jìn)行檢驗(yàn)，保證Meta回歸模型得到的預(yù)測(cè)值與實(shí)地研究得到的樣本觀察值保證一致性。為了檢驗(yàn)Meta分析價(jià)值轉(zhuǎn)移模型的有效性，驗(yàn)證效益轉(zhuǎn)移的精確程度，通常由樣本外效益轉(zhuǎn)移誤差(Et)表示，計(jì)算公式：

Et=|(yo-ye)/yo|×100%。

(2)

式中：yo為實(shí)地研究得到的樣本觀察值；ye為Meta回歸模型得到的預(yù)測(cè)值。Et的值越小，表明效益轉(zhuǎn)移模型的效果越好。目前，國(guó)外相關(guān)研究表明，樣本平均轉(zhuǎn)移誤差范圍在20%～40%之間是可以被接受的[27-28]，Meta回歸模型可用于樣本外的價(jià)值轉(zhuǎn)移。

2 結(jié)果與分析

2.1 各影響因素的Meta回歸模擬

本研究使用matlab軟件，采用方程(1)的公式，對(duì)各樣本價(jià)值觀察值和表2中的影響因素進(jìn)行Meta回歸分析(見(jiàn)表3)。Meta回歸模型的決定性數(shù)(R2)為0.903，調(diào)整后的R2為0.709，模型擬合度較高。

表3 Meta回歸模型中各影響因素的回歸模擬結(jié)果

2.2 影響國(guó)家級(jí)森林類型自然保護(hù)區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值變化的因素

通過(guò)表3，對(duì)影響國(guó)家級(jí)森林類型自然保護(hù)區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值變化的因素進(jìn)行分析，可以得出：

價(jià)值評(píng)估方法對(duì)評(píng)估結(jié)果的影響——由Meta回歸模型的回歸結(jié)果可見(jiàn)，碳稅法的回歸系數(shù)在統(tǒng)計(jì)上有差異，說(shuō)明在其它條件不變情況下，使用碳稅法得到的價(jià)值評(píng)估值與影子工程法得到的價(jià)值評(píng)估值有差異。使用市場(chǎng)價(jià)值法和機(jī)會(huì)成本法得到的價(jià)值評(píng)估值，高于其他價(jià)值評(píng)估方法；使用碳稅成本法和旅行費(fèi)用法得到的價(jià)值評(píng)估值，低于其他評(píng)估方法。

生態(tài)系統(tǒng)服務(wù)類型對(duì)評(píng)估結(jié)果的影響——在國(guó)家級(jí)森林類型自然保護(hù)區(qū)所提供的生態(tài)系統(tǒng)服務(wù)類型中，營(yíng)養(yǎng)物質(zhì)積累、保護(hù)生物多樣性、休閑游憩的回歸系數(shù)在統(tǒng)計(jì)上有差異，表明在其他條件不變情況下，營(yíng)養(yǎng)物質(zhì)積累、保護(hù)生物多樣性、休閑游憩的價(jià)值與固碳釋氧的價(jià)值評(píng)估值有差異。保護(hù)生物多樣性、森林游憩的經(jīng)濟(jì)價(jià)值顯著高于其他生態(tài)系統(tǒng)服務(wù)類型，土壤保持的生態(tài)系統(tǒng)服務(wù)價(jià)值最低，這與已有的研究結(jié)果一致[29]。

所處的植被分區(qū)對(duì)評(píng)估結(jié)果的影響——在植被分區(qū)中，暖溫帶落葉闊葉林區(qū)、溫帶草原區(qū)、溫帶荒漠區(qū)域、溫帶針闊混交林區(qū)的回歸系數(shù)在統(tǒng)計(jì)上有差異，說(shuō)明在其他條件不變情況下，處于暖溫帶落葉闊葉林區(qū)、溫帶草原區(qū)、溫帶荒漠區(qū)、溫帶針闊混交林區(qū)的保護(hù)區(qū)單位面積生態(tài)系統(tǒng)服務(wù)價(jià)值，與處于亞熱帶常綠闊葉林區(qū)的保護(hù)區(qū)單位面積生態(tài)系統(tǒng)服務(wù)價(jià)值有差異。處于溫帶荒漠區(qū)的保護(hù)區(qū)單位面積生態(tài)系統(tǒng)服務(wù)價(jià)值最低，這與荒漠區(qū)的環(huán)境相對(duì)惡劣存在一定關(guān)系。

國(guó)家級(jí)森林類型自然保護(hù)區(qū)成立時(shí)間對(duì)評(píng)估結(jié)果的影響——國(guó)家級(jí)森林類型自然保護(hù)區(qū)成立時(shí)間這一影響因素的回歸系數(shù)為正，說(shuō)明保護(hù)區(qū)成立的時(shí)間越晚，保護(hù)區(qū)內(nèi)每公頃的生態(tài)系統(tǒng)服務(wù)價(jià)值越高。國(guó)家級(jí)自然保護(hù)區(qū)評(píng)審標(biāo)準(zhǔn)不斷發(fā)展與完善，對(duì)入選的國(guó)家級(jí)自然保護(hù)區(qū)的自然屬性、可保護(hù)屬性、保護(hù)管理基礎(chǔ)有了更高的要求，生態(tài)系統(tǒng)服務(wù)功能也會(huì)有一定的相應(yīng)提高[30]。

受益人口數(shù)量對(duì)評(píng)估結(jié)果的影響——受益人口數(shù)量這一影響因素的回歸系數(shù)為正，說(shuō)明國(guó)家級(jí)森林類型自然保護(hù)區(qū)的受益人口數(shù)量越多，人們對(duì)保護(hù)區(qū)所產(chǎn)生的生態(tài)系統(tǒng)服務(wù)需求越大，保護(hù)區(qū)內(nèi)每公頃生態(tài)系統(tǒng)服務(wù)價(jià)值也相應(yīng)增加，但是并不顯著，這也說(shuō)明了人與環(huán)境的密切關(guān)系，人類的生存與發(fā)展離不開(kāi)良好生態(tài)環(huán)境的支撐。本研究采用對(duì)數(shù)形式的連續(xù)變量，其回歸系數(shù)表示在其他條件不變情況下，影響因素的變化引起的價(jià)值觀察值的變化。受益人口數(shù)量的回歸系數(shù)為0.017，說(shuō)明國(guó)家級(jí)森林類型自然保護(hù)區(qū)所處的地級(jí)市人口數(shù)量每增加1%，保護(hù)區(qū)內(nèi)每公頃生態(tài)系統(tǒng)服務(wù)價(jià)值會(huì)增加1.7%。

保護(hù)區(qū)面積對(duì)評(píng)估結(jié)果的影響——保護(hù)區(qū)面積這一影響因素的回歸系數(shù)為負(fù)，說(shuō)明在國(guó)家級(jí)森林類型自然保護(hù)區(qū)內(nèi)每公頃生態(tài)系統(tǒng)服務(wù)價(jià)值存在規(guī)模收益遞減效應(yīng)，這個(gè)結(jié)果與已有的研究結(jié)果一致[31]。并且已有的研究結(jié)果表明，這種規(guī)模收益遞減效應(yīng)會(huì)隨著濕地面積的增加而呈現(xiàn)出幾何級(jí)數(shù)的減弱，因此，保護(hù)區(qū)的生態(tài)系統(tǒng)服務(wù)總價(jià)值仍然會(huì)隨著保護(hù)區(qū)面積的增大而增長(zhǎng)。

區(qū)域經(jīng)濟(jì)發(fā)展?fàn)顩r對(duì)評(píng)估結(jié)果的影響——區(qū)域經(jīng)濟(jì)發(fā)展?fàn)顩r采用人均GDP衡量，人均GDP這一影響因素的回歸系數(shù)為負(fù)，說(shuō)明區(qū)域經(jīng)濟(jì)發(fā)展越好，保護(hù)區(qū)內(nèi)每公頃生態(tài)系統(tǒng)服務(wù)價(jià)值反而減少。國(guó)家級(jí)自然保護(hù)區(qū)的評(píng)審標(biāo)準(zhǔn)之一是自然性，即保護(hù)區(qū)基本處于自然狀態(tài)，人為干擾極少。保護(hù)區(qū)多選擇了原真性較好的地方，這些地方人類活動(dòng)相對(duì)較少，經(jīng)濟(jì)開(kāi)發(fā)也會(huì)相對(duì)較晚。此外，在經(jīng)濟(jì)發(fā)展較好的地區(qū)，往往生態(tài)旅游開(kāi)發(fā)也較早，然而在開(kāi)發(fā)過(guò)程中，時(shí)常會(huì)存在管理不足以及服務(wù)不到位等問(wèn)題，這也會(huì)降低單位面積內(nèi)生態(tài)系統(tǒng)服務(wù)的價(jià)值[32]。

2.3 Meta分析價(jià)值轉(zhuǎn)移模型的有效性

根據(jù)matlab軟件分析的結(jié)果，采用方程(2)對(duì)Meta分析價(jià)值轉(zhuǎn)移模型的有效性進(jìn)行檢驗(yàn)。結(jié)果顯示，效益轉(zhuǎn)移誤差范圍均在59.02%及以下，平均轉(zhuǎn)移誤差為24.99%。目前國(guó)內(nèi)外相關(guān)研究表明，平均誤差范圍在20%～40%之間是可以接受的，本研究的平均轉(zhuǎn)移誤差在此范圍之內(nèi)，說(shuō)明本研究所構(gòu)建的模型用于樣本外價(jià)值轉(zhuǎn)移有效性良好。

2.4 國(guó)家級(jí)森林類型自然保護(hù)區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值的模擬

依據(jù)本研究所得出的Meta回歸模型可得，我國(guó)國(guó)家級(jí)森林類型自然保護(hù)區(qū)所提供的單位面積生態(tài)系統(tǒng)服務(wù)經(jīng)濟(jì)價(jià)值為70 697.59元·hm-2·a-1；從保護(hù)區(qū)提供的各項(xiàng)生態(tài)系統(tǒng)服務(wù)功能價(jià)值看，各個(gè)生態(tài)系統(tǒng)服務(wù)功能單位面積平均價(jià)值量，由大到小依次為保護(hù)生物多樣性(19 117.42元·hm-2·a-1，占總價(jià)值比例27.04%)、森林游憩(14 439.74元·hm-2·a-1，占總價(jià)值比例20.42%)、營(yíng)養(yǎng)物質(zhì)積累(8 886.24元·hm-2·a-1，占總價(jià)值比例12.57%)、涵養(yǎng)水源(6 785.23元·hm-2·a-1，占總價(jià)值比例9.60%)、凈化空氣(6 550.30元·hm-2·a-1，占總價(jià)值比例9.27%)、林產(chǎn)品(5 883.46元·hm-2·a-1，占總價(jià)值比例8.32%)、固碳釋氧(4 889.51元·hm-2·a-1，占總價(jià)值比例6.92%)、土壤保持(4 145.69元·hm-2·a-1，占總價(jià)值比例5.86%)。在單位面積保護(hù)區(qū)所提供的各項(xiàng)生態(tài)系統(tǒng)服務(wù)價(jià)值中，保護(hù)生物多樣性價(jià)值所占比例最高，凈化空氣價(jià)值所占比例最低。

3 討論

樣本文獻(xiàn)數(shù)量、質(zhì)量對(duì)模型構(gòu)建的影響：在Meta分析價(jià)值轉(zhuǎn)移模型構(gòu)建中，已有研究的數(shù)量和質(zhì)量對(duì)于模型的構(gòu)建具有重要作用[33]。在本研究中，收集到的樣本文獻(xiàn)相對(duì)較少，文獻(xiàn)評(píng)估的國(guó)家級(jí)森林類型自然保護(hù)區(qū)多處在亞熱帶常綠闊葉林區(qū)、暖溫帶落葉闊葉林區(qū)、溫帶草原區(qū)、溫帶荒漠區(qū)，其他區(qū)域的研究文獻(xiàn)較少，甚至沒(méi)有；此外，樣本文獻(xiàn)的質(zhì)量通常也良莠不齊，對(duì)生態(tài)系統(tǒng)服務(wù)概念的界定、物質(zhì)量評(píng)估指標(biāo)以及價(jià)值評(píng)估方法等也不盡相同，降低了不同文獻(xiàn)間評(píng)估結(jié)果的可比性，導(dǎo)致在價(jià)值轉(zhuǎn)移的實(shí)地應(yīng)用中存在一定的困難。

影響因素的數(shù)據(jù)集大小對(duì)模型構(gòu)建的影響：在一定程度上，數(shù)據(jù)集的大小會(huì)造成區(qū)域生態(tài)系統(tǒng)服務(wù)功能價(jià)值的差異。本研究所選取的影響因素的數(shù)量相對(duì)于國(guó)外的研究相比仍然偏少，在影響因素的選取中未將保護(hù)區(qū)的生物物理特征、氣候變化以及人類活動(dòng)等影響因素考慮在內(nèi)，因此模型的擬合精度還不夠高。

樣本選擇偏差對(duì)模型構(gòu)建的影響：在樣本文獻(xiàn)的收集過(guò)程中，有的保護(hù)區(qū)在價(jià)值評(píng)估的過(guò)程中，用同樣的方法做研究卻得到了不同的評(píng)估結(jié)果。因此，本研究在選擇入選文獻(xiàn)時(shí)，多是考慮評(píng)估的時(shí)間、論文期刊的質(zhì)量等因素進(jìn)行了主觀選擇。這種對(duì)樣本文獻(xiàn)的篩選方式，會(huì)產(chǎn)生一定的“樣本選擇”偏差。

規(guī)范國(guó)家級(jí)森林類型自然保護(hù)區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)估的標(biāo)準(zhǔn)：目前，有大量實(shí)證研究表明，Meta分析價(jià)值轉(zhuǎn)移方法已成為國(guó)外資源價(jià)值評(píng)估領(lǐng)域關(guān)注的熱點(diǎn)之一，將Meta分析價(jià)值轉(zhuǎn)移方法應(yīng)用于不同尺度、不同類型的生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)估中[34-35]，涵蓋了對(duì)森林、濕地、草地、城市空地、珊瑚礁等生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)估[36-37]。為了促進(jìn)Meta分析在中國(guó)國(guó)家級(jí)森林類型自然保護(hù)區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)估中的應(yīng)用，提高該方法的可操作性、模型的有效性，有必要在今后的研究中，加強(qiáng)對(duì)我國(guó)國(guó)家級(jí)森林類型自然保護(hù)區(qū)生態(tài)系統(tǒng)服務(wù)的形成、環(huán)境變化以及人類活動(dòng)等影響因素的研究；并在此基礎(chǔ)上，參考《森林生態(tài)系統(tǒng)服務(wù)功能評(píng)估規(guī)范》，加強(qiáng)對(duì)保護(hù)區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)估的理論和實(shí)際研究，提高相關(guān)研究的數(shù)量和質(zhì)量，增強(qiáng)Meta分析價(jià)值轉(zhuǎn)移模型的可操作性、有效性。

Meta分析價(jià)值轉(zhuǎn)移模型的應(yīng)用展望：系統(tǒng)整體地開(kāi)展我國(guó)國(guó)家級(jí)森林類型自然保護(hù)區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值的評(píng)估，將為我國(guó)正在開(kāi)展的自然保護(hù)地整合優(yōu)化提供參考。依據(jù)本研究，可嘗試將得到的Meta分析價(jià)值轉(zhuǎn)移模型應(yīng)用到我國(guó)其他未開(kāi)展評(píng)估的國(guó)家級(jí)森林類型自然保護(hù)區(qū)中，進(jìn)行相應(yīng)的生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)估，節(jié)約評(píng)估成本、評(píng)估時(shí)間。

4 結(jié)論

本文依據(jù)收集到的31篇關(guān)于我國(guó)國(guó)家級(jí)森林類型自然保護(hù)區(qū)生態(tài)系統(tǒng)服務(wù)功能價(jià)值評(píng)估文獻(xiàn)，建立價(jià)值轉(zhuǎn)移數(shù)據(jù)庫(kù)，選取恰當(dāng)?shù)腗eta回歸模型、多元回歸分析方法，構(gòu)建適合我國(guó)國(guó)家級(jí)森林類型自然保護(hù)區(qū)生態(tài)系統(tǒng)服務(wù)功能價(jià)值評(píng)估的Meta分析價(jià)值轉(zhuǎn)移模型，并針對(duì)該模型價(jià)值轉(zhuǎn)移的有效性進(jìn)行檢驗(yàn)。

在我國(guó)國(guó)家級(jí)森林類型自然保護(hù)區(qū)生態(tài)系統(tǒng)服務(wù)功能價(jià)值評(píng)估中，選取的評(píng)估方法、生態(tài)系統(tǒng)服務(wù)功能類型、保護(hù)區(qū)所處的植被分區(qū)、保護(hù)區(qū)建立的時(shí)間，會(huì)較明顯地影響保護(hù)區(qū)生態(tài)系統(tǒng)服務(wù)功能的價(jià)值變化。

在我國(guó)國(guó)家級(jí)森林類型自然保護(hù)區(qū)所提供的生態(tài)系統(tǒng)服務(wù)功能類型中，保護(hù)生物多樣性價(jià)值、休閑游憩價(jià)值最高，凈化空氣的價(jià)值最低。

通過(guò)對(duì)Meta分析價(jià)值轉(zhuǎn)移模型的有效性檢驗(yàn)結(jié)果看，樣本外效益轉(zhuǎn)移的誤差范圍均在59.02%及以下，誤差均值為24.99%，誤差范圍在合理的區(qū)間范圍內(nèi)，模型可用于樣本外效益轉(zhuǎn)移；本研究得到的Meta分析價(jià)值轉(zhuǎn)移模型，可作為一種快速且有效的評(píng)估方法用在我國(guó)國(guó)家級(jí)森林類型自然保護(hù)區(qū)生態(tài)系統(tǒng)服務(wù)功能價(jià)值評(píng)估中。