廣西稻田生態(tài)系統(tǒng)服務(wù)價值評估及其時空演變特征

蔣琳莉，黃好欽，何蘇榕，熊 娜，望飛勇

(1.五邑大學(xué)經(jīng)濟(jì)管理學(xué)院，廣東 江門 529020；2.廣西民族大學(xué)管理學(xué)院，南寧 530006；3.廣西農(nóng)業(yè)科學(xué)院園藝研究所，南寧 530007)

【研究意義】糧食安全與生態(tài)平衡是稻田生態(tài)系統(tǒng)服務(wù)價值的核心內(nèi)涵。在傳統(tǒng)生產(chǎn)模式下，水稻產(chǎn)量的提升通常伴隨著化肥、農(nóng)藥和農(nóng)膜等化學(xué)投入品的大量使用，稻田面源污染不斷加劇，甲烷(CH4)、氧化亞氮(N2O)和二氧化碳(CO2)等農(nóng)業(yè)溫室氣體過量排放，進(jìn)而導(dǎo)致稻田生態(tài)效益不斷減弱，并在一定程度上加快全球氣候變暖進(jìn)程[1]。廣西是水稻主產(chǎn)省份之一，稻田生態(tài)環(huán)境污染時有發(fā)生，稻田生態(tài)系統(tǒng)服務(wù)價值波動較大[2]。鑒于此，有效評估廣西稻田生態(tài)系統(tǒng)服務(wù)價值的變化情況，對政府部門推進(jìn)農(nóng)業(yè)綠色發(fā)展以維護(hù)糧食安全和生態(tài)平衡具有重要意義?！厩叭搜芯窟M(jìn)展】稻田是一種人工濕地生態(tài)系統(tǒng)，其服務(wù)價值可劃分為經(jīng)濟(jì)服務(wù)價值、文化服務(wù)價值和生態(tài)服務(wù)價值[3-5]，其中，經(jīng)濟(jì)服務(wù)價值表現(xiàn)為生產(chǎn)初級產(chǎn)品(稻谷和秸稈)，文化服務(wù)價值是指稻田生態(tài)系統(tǒng)帶給人類的景觀文化價值，生態(tài)服務(wù)價值包括氣體調(diào)節(jié)、氣候調(diào)節(jié)、涵養(yǎng)水源、水質(zhì)凈化、土壤保持和維持生物多樣性方面的價值。但水稻植株生長特性及傳統(tǒng)水稻種植模式產(chǎn)生的污染，也會引起稻田生態(tài)系統(tǒng)產(chǎn)生負(fù)面生態(tài)服務(wù)價值，分別為溫室氣體排放和農(nóng)用化學(xué)品污染[4，6]。同時，稻田生態(tài)效益具有外部特征，而稻田生態(tài)服務(wù)交易市場尚未建立，無法通過市場交易價格直接評估其價值[7-8]，已有研究主要通過市場價值法、替代成本法和成果參照法等方法評估稻田生態(tài)系統(tǒng)的服務(wù)價值[9-11]。在生態(tài)系統(tǒng)服務(wù)價值量方面，鄭春利等[12]評價以景觀文化價值為核心的元陽梯田景區(qū)稻田生態(tài)價值為1.79億元，Yosihiro[13]測算日本稻田生態(tài)系統(tǒng)調(diào)節(jié)功能總量約728億美元。部分學(xué)者也評估了稻田負(fù)面生態(tài)服務(wù)價值，如劉利花等[6]評估蘇州市稻田正面生態(tài)服務(wù)價值總量為68.49億元，負(fù)面生態(tài)服務(wù)價值為7.7億元，薛寶林等[14]測算湖南省稻田生態(tài)系統(tǒng)服務(wù)凈價值為2.47萬元/hm2，負(fù)面生態(tài)服務(wù)價值為1.83萬元/hm2。在生態(tài)系統(tǒng)服務(wù)價值量的時空演變特征方面，相關(guān)學(xué)者對農(nóng)田、草原和森林等各類農(nóng)業(yè)生態(tài)系統(tǒng)進(jìn)行了研究，如熊文等[15]研究指出，2010—2020年長江經(jīng)濟(jì)帶重點區(qū)域農(nóng)田生態(tài)系統(tǒng)服務(wù)價值量呈現(xiàn)明顯的增長趨勢，穆松林[16]研究認(rèn)為，1982—2014年內(nèi)蒙古全區(qū)溫帶草原生態(tài)系統(tǒng)服務(wù)價值在空間分布上由東北向西南逐漸降低，楊青和劉耕源[17]研究表明，2012年京津冀森林生態(tài)系統(tǒng)各類服務(wù)價值空間分布格局呈現(xiàn)西北高、東南低的特征?！颈狙芯壳腥朦c】目前，全面揭示廣西稻田生態(tài)系統(tǒng)服務(wù)價值時空動態(tài)變化過程的研究鮮見報道?！緮M解決的關(guān)鍵問題】利用廣西區(qū)級層面數(shù)據(jù)，評估2011—2020年廣西全區(qū)稻田生態(tài)系統(tǒng)服務(wù)價值量，揭示近十年來各類服務(wù)價值量的變化趨勢；利用廣西市級層面數(shù)據(jù)，評估2011—2019年廣西14個地級市稻田生態(tài)系統(tǒng)服務(wù)價值的時空演變特征，以及各類生態(tài)服務(wù)價值的空間分布格局，以期為廣西各地政府制定稻田生態(tài)環(huán)境保護(hù)相關(guān)政策提供參考依據(jù)。

1 材料與方法

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

選取廣西稻田生態(tài)系統(tǒng)經(jīng)濟(jì)文化服務(wù)價值(產(chǎn)生初級產(chǎn)品和景觀文化)和生態(tài)服務(wù)價值(氣體調(diào)節(jié)、氣候調(diào)節(jié)、涵養(yǎng)水源、水質(zhì)凈化、土壤保持和維持生物多樣性6項正面價值，以及溫室氣體排放和農(nóng)用化學(xué)品污染2項負(fù)面價值)[3-5]，采用市場價值法、替代成本法和成果參照法等評估廣西稻田生態(tài)系統(tǒng)服務(wù)價值[9-11]。

廣西水稻種植面積和稻谷產(chǎn)量的基礎(chǔ)統(tǒng)計數(shù)據(jù)來源于2011—2020年廣西統(tǒng)計年鑒、各地級市統(tǒng)計年鑒及國民經(jīng)濟(jì)和社會發(fā)展統(tǒng)計公報；水稻價格來源于國家發(fā)展和改革委員會；早稻、中稻和晚稻生長周期及非水淹期天數(shù)來源于廣西農(nóng)業(yè)科學(xué)院水稻研究所；針對部分年份秸稈價格、農(nóng)用水價格和有機(jī)肥價格等相關(guān)數(shù)據(jù)缺失的問題，本研究采用2011—2020年廣西全區(qū)消費者價格指數(shù)估算結(jié)果作為對應(yīng)年份的數(shù)據(jù)指標(biāo)；碳稅率、農(nóng)用電價格和污水處理市場價格等數(shù)據(jù)來源于相關(guān)網(wǎng)站。

1.2 服務(wù)價值及其測算方法

1.2.1 經(jīng)濟(jì)文化服務(wù)價值 經(jīng)濟(jì)文化服務(wù)價值(V1)是稻田生產(chǎn)初級產(chǎn)品的經(jīng)濟(jì)價值與景觀文化服務(wù)價值之和[3]。稻田生態(tài)系統(tǒng)生產(chǎn)的初級產(chǎn)品包括稻谷和秸稈，其中稻谷是人類的糧食來源，秸稈可作為造紙、纖維制品和動物飼料原料等，具有一定經(jīng)濟(jì)價值[18]。同時，廣西以稻田為基礎(chǔ)的民族文化、獨特的稻田自然景觀(如桂林龍脊梯田、柳州苗寨梯田和靖西巴澤梯田等)能帶給人們美學(xué)感受等非物質(zhì)價值。稻田生產(chǎn)初級產(chǎn)品的經(jīng)濟(jì)價值(V11)采用市場價值法計算[6]，景觀文化服務(wù)價值(V12)采用成果參照法進(jìn)行測算，以單位面積景觀文化價值量[19]作為計算標(biāo)準(zhǔn)。

1.2.2 正面生態(tài)服務(wù)價值 (1)氣體調(diào)節(jié)。稻田生態(tài)系統(tǒng)氣體調(diào)節(jié)價值(V2)為固定CO2、釋放O2和吸收有害氣體價值的總和，參考薛寶林等[14]的方法進(jìn)行測算。其中，固定CO2的價值取瑞典碳稅法和造林成本法計算結(jié)果的平均值(V21)；釋放O2的價值取制氧工業(yè)成本法和造林成本法計算結(jié)果的平均值(V22)；吸收有害氣體的價值(V23)包括吸收粉塵、二氧化硫(SO2)、氮氧化物(NOx)和氟化氫(HF)。

(2)氣候調(diào)節(jié)。在夏季，水稻植株的蒸騰作用和稻田的水分蒸發(fā)均會吸收周圍環(huán)境的熱量，降低周圍空氣的溫度，緩解“熱島效應(yīng)”[20]。氣候調(diào)節(jié)價值(V3)可用農(nóng)村居民節(jié)約使用空調(diào)的效益來衡量[21]。

(3)涵養(yǎng)水源。在水稻生長過程中，其根系大部分時間被水層覆蓋，體現(xiàn)了稻田良好的保水性。稻田涵養(yǎng)水源價值(V4)參考胡曉燕[4]的測算方法，以土壤中滲入的水量乘以農(nóng)用水的價格替代。

(4)水質(zhì)凈化。稻田生態(tài)系統(tǒng)對污水的處理不僅有吸附和沉降物理凈化過程，還有氧化和還原分解化學(xué)凈化過程[22]，稻田水質(zhì)凈化價值(V5)參考李松[23]的測算方法，以農(nóng)村生活污水的凈化成本替代。

(5)土壤保持。稻田生態(tài)系統(tǒng)具有攔截泥沙、防止土壤侵蝕及增加土壤有機(jī)質(zhì)等土壤保持功能[24]，該功能價值(V6)參考肖玉和謝高地[25]的測算方法，以土壤有機(jī)質(zhì)輸入量(V61)和輸出量(V62)的差值乘以有機(jī)肥的市場價格來衡量。稻田土壤有機(jī)質(zhì)主要來源于水稻根系的沉析作用和水稻收割后殘留稻田的秸稈。同時，水稻與土壤微生物均會消耗土壤有機(jī)質(zhì)，其消耗量主要通過二者的CO2和CH4排放量計算。

(6)維持生物多樣性。稻田作為人工濕地生態(tài)系統(tǒng)，不僅向人類提供糧食，同時也是植物、動物和土壤微生物等賴以生存的生態(tài)環(huán)境。該項服務(wù)價值(V7)采用成果參照法進(jìn)行測算，以單位面積維持生物多樣性價值量作為計算標(biāo)準(zhǔn)[19]。

1.2.3 負(fù)面生態(tài)服務(wù)價值 (1)溫室氣體排放。稻田生態(tài)系統(tǒng)向大氣排放的溫室氣體主要包括CH4、N2O和CO2。本研究將稻田生態(tài)系統(tǒng)排放的CH4和N2O負(fù)面價值統(tǒng)一折算為排放CO2的負(fù)面價值，即1.00 kg CH4與24.50 kg CO2所產(chǎn)生的溫室效應(yīng)相同，1.00 kg N2O與320.00 kg CO2所產(chǎn)生的溫室效應(yīng)相同[6]。稻田溫室氣體排放的負(fù)面價值(V8)參考劉利花等[6]的方法測算，取瑞典碳稅法和造林成本法計算結(jié)果的平均值。

(2)農(nóng)用化學(xué)品污染。水稻在種植過程中使用化學(xué)肥料易造成周邊水體富營養(yǎng)化，農(nóng)藥的施用可能導(dǎo)致土壤與空氣污染，同時危及稻田生態(tài)系統(tǒng)的生物多樣性，該項負(fù)面價值(V9)參考劉利花等[6]的方法，以農(nóng)用化學(xué)品造成的污染成本乘以水稻種植面積進(jìn)行測算。

1.3 統(tǒng)計分析

采用Excel 2019進(jìn)行廣西稻田生態(tài)系統(tǒng)服務(wù)價值數(shù)據(jù)測算與評估，以ArcGIS 10.2繪制廣西稻田生態(tài)服務(wù)價值的空間分布圖。

2 結(jié)果與分析

2.1 廣西稻田生態(tài)系統(tǒng)服務(wù)價值的時間變化趨勢

2.1.1 各類服務(wù)價值總量的時間變化趨勢 從圖1-A可看出，在2011—2020年，廣西稻田生態(tài)系統(tǒng)服務(wù)價值總量呈先升高后降低再升高的變化趨勢。其中，2011年的服務(wù)價值總量為488.52億元，2020年為568.17億元，上升16.30%，數(shù)值變化不明顯；2011—2014年為服務(wù)價值總量增長時期，從488.52億元上升至629.62億元，增長28.88%；2014—2019年為服務(wù)價值總量下降時期，從629.62億元下降至546.78億元，下降13.16%；2020年服務(wù)價值總量較2019年略有增長，為568.17億元。在2011—2020年，經(jīng)濟(jì)文化價值的變化趨勢與服務(wù)價值總量的變化趨勢基本一致，從256.48億元上升至284.08億元，增長10.76%。具體而言，經(jīng)濟(jì)文化價值由2011年的256.48億元上升至2014年的342.02億元，增長33.35%；隨后下降至2019年的276.06億元，最后又回升至2020年的284.08億元。在2011—2020年，生態(tài)凈價值(正面生態(tài)服務(wù)價值與負(fù)面生態(tài)服務(wù)價值之差)總體上保持穩(wěn)定，但也存在一定幅度的波動。其中，2011—2012年為生態(tài)凈價值增長時期，從232.04億元上升至最大值306.54億元，增長32.11%；2013—2020年為生態(tài)凈價值穩(wěn)定時期，均值為277.57億元。

2.1.2 各類生態(tài)服務(wù)價值的時間變化趨勢 從圖1-B可看出，2011—2020年廣西稻田生態(tài)系統(tǒng)正面生態(tài)服務(wù)價值的變化趨勢與生態(tài)凈價值的變化趨勢基本一致。其中，2011—2012年正面生態(tài)服務(wù)價值從最小值370.76億元上升至最大值448.00億元，增長20.83%；2013—2020年正面生態(tài)服務(wù)價值較穩(wěn)定，均值為420.78億元；2011—2020年負(fù)面生態(tài)服務(wù)價值持續(xù)保持穩(wěn)定，均值為142.59億元。

C：實線對應(yīng)左側(cè)縱坐標(biāo)軸，虛線對應(yīng)右側(cè)縱坐標(biāo)軸C：The solid lines corresponded to the left vertical axis，and the dotted lines corresponded to the right vertical axis圖1 2011—2020年廣西稻田生態(tài)系統(tǒng)服務(wù)價值的變化趨勢Fig.1 Change trend of paddy field ecosystem service value in Guangxi from 2011 to 2020

從圖1-C可看出，2011—2020年廣西稻田生態(tài)系統(tǒng)正面生態(tài)服務(wù)功能中水質(zhì)凈化的價值位于6項正面生態(tài)服務(wù)價值首位，但其數(shù)值波動較明顯。其中，2011—2012年水質(zhì)凈化價值由最小值156.95億元上升至最大值225.75億元，增長43.84%；2012—2016年水質(zhì)凈化價值總體上呈下降趨勢，由225.75億元下降至179.08億元，降低20.67%；2016—2020年水質(zhì)凈化價值呈上升趨勢，由179.08億元上升至195.64億元，增長9.25%。2011—2020年涵養(yǎng)水源價值、氣體調(diào)節(jié)價值、土壤保持價值、維持生物多樣性價值和氣候調(diào)節(jié)價值均保持穩(wěn)定，均值依次為131.15億、35.62億、33.84億、15.97億和6.59億元。

從圖1-D可看出，2011—2018年為廣西稻田生態(tài)系統(tǒng)溫室氣體排放負(fù)面價值穩(wěn)定時期，均值為130.86億元，而2018—2020年略有上升，由128.85億元上升至138.05億元，增長7.14%。同時，2011—2020年農(nóng)用化學(xué)品污染負(fù)面價值保持穩(wěn)定，數(shù)值波動不明顯，均值為11.73億元。

2.2 廣西稻田生態(tài)系統(tǒng)服務(wù)價值的時空演變特征

2.2.1 服務(wù)價值總量和生態(tài)服務(wù)價值的時空演變特征 (1)服務(wù)價值總量的時空演變特征 由圖2-A和圖2-B可知，2011—2015年廣西稻田生態(tài)系統(tǒng)服務(wù)價值總量的空間集聚程度明顯增強(qiáng)；服務(wù)價值總量的空間格局呈由桂東南的玉林市向外降低的階梯狀分布，逐漸演變?yōu)槟媳钡貐^(qū)較高、東西地區(qū)(以賀州市、百色市、來賓市為主)較低的分布特征；在桂東北和桂東南表現(xiàn)出Ⅱ級區(qū)域向西收縮、Ⅲ級區(qū)域向西擴(kuò)張的空間演變過程，其中，桂東北柳州市、桂林市和桂東南梧州市服務(wù)價值總量的空間演變過程由Ⅱ級區(qū)域轉(zhuǎn)變?yōu)棰蠹墔^(qū)域，服務(wù)價值總量有所提高。而在桂南地區(qū)，服務(wù)價值總量的空間演變過程由Ⅳ級區(qū)域全部轉(zhuǎn)變?yōu)棰跫墔^(qū)域，主要集聚在南寧市、貴港市和玉林市，服務(wù)價值總量居于全區(qū)前三位。2015—2019年廣西稻田生態(tài)系統(tǒng)服務(wù)價值總量的空間集聚程度逐步減弱，其空間分布格局逐漸轉(zhuǎn)變?yōu)橐阅蠈幨袨橹行南蛩闹苓f減的階梯狀分布(圖2-B，2-C)。其中，Ⅱ級區(qū)域向桂北和桂東城市蔓延，桂北地區(qū)的柳州市、桂林市和桂東地區(qū)的梧州市由Ⅲ級區(qū)域再次轉(zhuǎn)變?yōu)棰蚣墔^(qū)域，服務(wù)價值總量有所降低。同時，桂南地區(qū)的貴港市和玉林市由Ⅴ級區(qū)域轉(zhuǎn)變?yōu)棰艏墔^(qū)域。從總體上看，2011—2019年廣西稻田生態(tài)系統(tǒng)服務(wù)價值總量以百色市、來賓市和梧州市為分界線，包括上述3個城市在內(nèi)的桂北地區(qū)服務(wù)價值總量的空間分布格局基本保持不變，均為Ⅱ級區(qū)域，且占比達(dá)42.86%；在分界線以南，服務(wù)價值總量中心由東南邊界城市玉林市轉(zhuǎn)移至中心城市南寧市，形成了以南寧市為中心向四周遞減的階梯狀分布格局?？梢?，近十年來廣西稻田生態(tài)系統(tǒng)服務(wù)價值總量的空間分布格局基本上保持南高北低的態(tài)勢，其中心在桂南地區(qū)的城市間轉(zhuǎn)移。

(2)生態(tài)服務(wù)凈價值的時空演變特征 從圖2D～F可看出，2011—2015年廣西稻田生態(tài)系統(tǒng)生態(tài)服務(wù)凈價值基本保持南高北低的態(tài)勢，桂南Ⅴ級區(qū)域增多，桂北Ⅲ級區(qū)域增多，但仍以Ⅱ級區(qū)域為主；在桂南地區(qū)，Ⅲ級區(qū)域和Ⅳ級區(qū)域均轉(zhuǎn)變?yōu)棰跫墔^(qū)域，主要集聚在南寧市、貴港市和玉林市，且呈現(xiàn)連片分布；桂北邊緣地區(qū)生態(tài)服務(wù)凈價值由Ⅱ級區(qū)域連片分布逐步演變?yōu)樽詵|向西從Ⅲ級區(qū)域至Ⅰ級區(qū)域逐級遞減的階梯狀分布；2015—2019年桂北Ⅱ級區(qū)域以河池市為中心向東西兩側(cè)擴(kuò)張，其西側(cè)的百色市由Ⅰ級區(qū)域升為Ⅱ級區(qū)域，其東側(cè)的柳州市和桂林市由Ⅲ級區(qū)域降為Ⅱ級區(qū)域，最終形成Ⅱ級區(qū)域連片分布的空間特征；在桂南地區(qū)，貴港市和玉林市均由Ⅴ級區(qū)域降為Ⅳ級區(qū)域，南寧市成為廣西全區(qū)唯一的Ⅴ級區(qū)域。從總體上看，在2011—2019年廣西稻田生態(tài)服務(wù)凈價值的空間分布格局中，桂北地區(qū)格局基本保持不變，以Ⅱ級區(qū)域連片分布為主，占總城市數(shù)的50%；在桂南地區(qū)，生態(tài)凈價值中心由玉林市向西北遷移至南寧市，形成以南寧市為中心向四周逐級遞減的空間分布特征。

(3)負(fù)面生態(tài)服務(wù)價值的時空演變特征 從圖2-G和圖2-H還可看出，2011—2015年廣西稻田生態(tài)系統(tǒng)負(fù)面生態(tài)服務(wù)價值基本保持南高北低的態(tài)勢，桂北Ⅲ級區(qū)域有所增加，桂南Ⅳ級區(qū)域分別向Ⅲ級區(qū)域和Ⅴ級區(qū)域轉(zhuǎn)變。其中，Ⅲ級區(qū)域自南向北擴(kuò)展，形成南北一線的空間分布特征，城市數(shù)量增加至4個，分別為欽州市、貴港市、柳州市和桂林市；Ⅱ級區(qū)域以梧州市、來賓市和河池市為主，形成東西一線的空間分布格局。同時，在桂南地區(qū)的貴港市由Ⅳ級區(qū)域降為Ⅲ級區(qū)域，而南寧市由Ⅳ級區(qū)域升至Ⅴ級區(qū)域，負(fù)面生態(tài)服務(wù)價值總體上呈現(xiàn)南高北低的分布特征。2015—2019年桂北邊緣城市柳州市和桂林市由Ⅲ級區(qū)域降為Ⅱ級區(qū)域(圖2-H，2-I)，使桂北地區(qū)Ⅱ級區(qū)域呈現(xiàn)連片分布的空間特征；桂南地區(qū)Ⅳ級區(qū)域向北蔓延，其中貴港市由Ⅲ級區(qū)域轉(zhuǎn)變?yōu)棰艏墔^(qū)域。綜上所述，2011—2019年廣西稻田生態(tài)系統(tǒng)負(fù)面生態(tài)服務(wù)價值的時空演變呈現(xiàn)周期性變化特征，2019年負(fù)面生態(tài)服務(wù)價值的空間分布格局與2011年基本一致。

圖2 2011—2019年廣西稻田生態(tài)系統(tǒng)服務(wù)價值的時空演變特征Fig.2 Spatial-temporal evolution characteristics of paddy field ecosystem service value in Guangxi from 2011 to 2019

2.2.2 各類生態(tài)服務(wù)價值的空間分布格局 (1)正面生態(tài)服務(wù)價值的空間分布格局 由圖3A～F可知，涵養(yǎng)水源價值與水質(zhì)凈化價值在2019年廣西稻田生態(tài)系統(tǒng)的正面生態(tài)服務(wù)價值中占據(jù)主導(dǎo)地位，主要分布在桂東南和桂北地區(qū)，在其他區(qū)域分布相對較少。其中，涵養(yǎng)水源價值主要分布在桂南地區(qū)的南寧市、玉林市、貴港市、欽州市及桂北地區(qū)的桂林市，價值量分別為20.01億、17.20億、15.54億、10.43億和16.47億元，合計占涵養(yǎng)水源價值總量的60.05%；水質(zhì)凈化價值空間分布特征與涵養(yǎng)水源價值的空間分布特征較一致，其價值量也主要分布于上述5個地級市，分別為29.18億、25.05億、22.63億、15.18億和23.84億元，合計占水質(zhì)凈化價值總量的60.05%。其余4項正面生態(tài)服務(wù)價值量相對較小，但其空間分布仍以桂南地區(qū)為主。其中，氣體調(diào)節(jié)價值主要分布在桂南地區(qū)的南寧市、玉林市、貴港市及桂北地區(qū)的桂林市，價值量分別為5.26億、5.09億、4.36億和4.82億元，合計占?xì)怏w調(diào)節(jié)價值總量的54.00%；氣候調(diào)節(jié)價值主要分布在桂南地區(qū)的貴港市、欽州市、玉林市、南寧市，桂西地區(qū)的百色市，以及桂北地區(qū)的桂林市，價值量分別為1.04億、0.59億、0.52億、0.51億、0.51億和0.60億元，合計占?xì)夂蛘{(diào)節(jié)價值總量的71.70%；土壤保持價值主要分布在桂南地區(qū)的南寧市、玉林市、貴港市及桂北地區(qū)的桂林市，價值量分別為5.22億、4.46億、4.03億和4.16億元，合計占土壤保持價值總量的52.17%；維持生物多樣性價值主要分布在桂南地區(qū)的南寧市、玉林市、貴港市及桂北地區(qū)的桂林市，價值量分別為2.44億、2.10億、1.89億和2.00億元，合計占維持生物多樣性價值總量的52.18%。

(2)負(fù)面生態(tài)服務(wù)價值的空間分布格局 在2項負(fù)面生態(tài)服務(wù)價值中，溫室氣體排放是主要來源，農(nóng)用化學(xué)品污染價值量占比相對較少。其中，溫室氣體排放負(fù)面價值主要分布在桂南地區(qū)的南寧市、玉林市、貴港市及桂北地區(qū)的桂林市，價值量分別為19.96億、17.16億、15.50億和16.44億元，合計占溫室氣體排放負(fù)面價值總量的52.19%；農(nóng)用化學(xué)品污染負(fù)面價值主要分布在桂南地區(qū)的南寧市、玉林市、貴港市、欽州市及桂北地區(qū)的桂林市，價值量分別為1.73億、1.68億、1.44億、0.92億和1.59億元，合計占農(nóng)用化學(xué)品污染負(fù)面價值總量的61.90%(圖3-G，3-H)。

圖3 2019年廣西稻田生態(tài)系統(tǒng)服務(wù)價值的正面與負(fù)面生態(tài)服務(wù)價值Fig.3 Positive and negative ecological service value of paddy field ecosystem in Guangxi in 2019

3 討 論

在稻田生態(tài)系統(tǒng)服務(wù)價值評估方面，劉利花等[6]評估的是2012年蘇州市稻田生態(tài)系統(tǒng)服務(wù)價值，薛寶林等[13]評估的是2020年湖南省稻田生態(tài)系統(tǒng)服務(wù)價值，涉及的均是單一區(qū)域和單一年度，而本研究在評估區(qū)域上涵蓋了廣西區(qū)級和市級2個空間層面，且時間跨度長達(dá)10年，能更全面地揭示稻田生態(tài)系統(tǒng)服務(wù)價值的時空變化特征。在西部地區(qū)，2008—2017年新疆農(nóng)田生態(tài)系統(tǒng)服務(wù)價值總量呈先升高再略微下降的變化趨勢[26]；在中部地區(qū)，2010—2020年湖南省農(nóng)田生態(tài)系統(tǒng)服務(wù)價值總量呈現(xiàn)顯著增長趨勢[14]；在東部地區(qū)，2005—2015年山東省聊城市農(nóng)田生態(tài)系統(tǒng)服務(wù)價值總量持續(xù)增長[27]。結(jié)合本研究2011—2020年廣西全區(qū)稻田生態(tài)系統(tǒng)服務(wù)價值總量變化不明顯，但總體上呈先升高后降低再升高變化趨勢的分析結(jié)果可發(fā)現(xiàn)，針對全國各區(qū)域農(nóng)田生態(tài)系統(tǒng)服務(wù)價值總量變化趨勢的研究結(jié)果較一致，均認(rèn)為隨著時間的變化，農(nóng)田生態(tài)系統(tǒng)服務(wù)價值總量總體上保持上升趨勢，表明近年來我國農(nóng)田生態(tài)修復(fù)與保護(hù)工程取得了較好成效。

本研究結(jié)果顯示，2011—2019年廣西稻田生態(tài)系統(tǒng)服務(wù)價值主要集中在桂南地區(qū)，且呈現(xiàn)南高北低的空間分布格局，應(yīng)與桂南地區(qū)的南寧市、玉林市、欽州市和貴港市等均為廣西水稻產(chǎn)區(qū)，其中南寧盆地、玉林盆地和貴港以西的郁江平原均為水稻生長的氣候適宜區(qū)有關(guān)。諸多學(xué)者對草原、森林和流域等各類農(nóng)業(yè)生態(tài)系統(tǒng)服務(wù)價值空間分布特征的研究結(jié)果表明，1982—2014年內(nèi)蒙古全區(qū)溫帶草原生態(tài)系統(tǒng)服務(wù)價值呈現(xiàn)由東北向西南逐漸降低的空間分布特征[16]，2012年京津冀森林生態(tài)系統(tǒng)各類服務(wù)價值呈現(xiàn)西北高、東南低的空間分布特征[17]，2005—2013年鄱陽湖流域生態(tài)服務(wù)價值空間分布不均衡，其中鄱陽湖及附近區(qū)域的生態(tài)服務(wù)價值遠(yuǎn)高于其他地區(qū)[28]。以上研究結(jié)果表明，不同農(nóng)業(yè)生態(tài)系統(tǒng)的服務(wù)價值空間分布格局具有差異性。其原因主要在于農(nóng)業(yè)生態(tài)系統(tǒng)類型受到地理與氣候因素的影響，且其所包含的動植物和農(nóng)業(yè)生產(chǎn)活動存在差異。本研究還發(fā)現(xiàn)，2019年廣西稻田生態(tài)系統(tǒng)的正面和負(fù)面生態(tài)服務(wù)價值呈現(xiàn)南北高、東西低的空間分布特征，與熊娜等[29]研究得出廣西水稻秸稈資源主要分布在其南、北地區(qū)的空間格局基本一致，均顯示廣西稻田生態(tài)系統(tǒng)的空間分布以桂南與桂北地區(qū)為主，說明本研究結(jié)論具有一定的合理性，可為廣西各地政府部門保護(hù)稻田生態(tài)環(huán)境、維護(hù)糧食產(chǎn)量安全和農(nóng)業(yè)生態(tài)平衡提供決策參考。

目前，廣西正在逐步推廣間歇灌溉、綠肥種植和生物農(nóng)藥施用等水稻綠色生產(chǎn)技術(shù)[30]。受限于數(shù)據(jù)的可獲得性，本研究忽略了水稻綠色生產(chǎn)技術(shù)在降低水稻溫室氣體排放量及減輕農(nóng)用化學(xué)品污染方面的貢獻(xiàn)。針對以上不足，在未來的研究中將進(jìn)一步完善稻田生態(tài)系統(tǒng)服務(wù)價值的指標(biāo)體系構(gòu)建，并通過大田試驗數(shù)據(jù)精確評估各項水稻綠色生產(chǎn)技術(shù)對稻田生態(tài)系統(tǒng)服務(wù)價值的提升效果。

4 結(jié) 論

2011—2020年廣西稻田生態(tài)系統(tǒng)服務(wù)價值總量變化不明顯，但總體上呈先升高后降低再升高的變化趨勢；在其時空演變進(jìn)程中，桂北地區(qū)稻田生態(tài)服務(wù)價值的空間分布格局基本保持不變，桂南地區(qū)稻田生態(tài)服務(wù)價值總量中心由東南邊界的玉林市轉(zhuǎn)移至中心城市南寧市，形成以南寧市為中心向四周遞減的階梯狀分布格局；2019年廣西稻田生態(tài)系統(tǒng)的正面和負(fù)面生態(tài)服務(wù)價值主要分布在桂南和桂北地區(qū)，形成南北高、東西低的空間分布格局。