郭祥光，陳永林，林建平

(贛南師范大學地理與環(huán)境工程學院，341000，江西，贛州)

0 引言

稀土作為一種不可再生的礦產(chǎn)資源，被譽為“現(xiàn)代工業(yè)的維生素”，在工業(yè)生產(chǎn)領(lǐng)域具有十分重要的戰(zhàn)略價值。在我國，隨著稀土需求的增加，稀土開采的加大，稀土開采無序、亂采濫挖等現(xiàn)象普遍存在，稀土礦區(qū)及周邊地區(qū)的土地利用現(xiàn)狀發(fā)生了深刻的變化，并由此產(chǎn)生了一系列生態(tài)環(huán)境問題[1]。

稀土礦區(qū)土地利用變化是由稀土開采造成的局部地區(qū)土地動態(tài)演化的過程，是人類活動影響區(qū)域生態(tài)系統(tǒng)最直接的表現(xiàn)形式，是生態(tài)系統(tǒng)服務(wù)價值變化的重要影響因素[2]。生態(tài)系統(tǒng)服務(wù)價值(ESV)是生態(tài)系統(tǒng)功能效用的價值表現(xiàn)，是衡量區(qū)域生態(tài)環(huán)境質(zhì)量的重要指標[3-5]。當前，國內(nèi)外學術(shù)界針對稀土礦區(qū)進行了一系列研究，主要包括利用遙感技術(shù)對稀土礦山進行土地利用動態(tài)監(jiān)測[6-7]和環(huán)境監(jiān)測[8-9]、稀土礦區(qū)土壤稀土元素污染[10-11]、稀土礦區(qū)土地修復[12-13]等，而對礦區(qū)土地利用變化及其對生態(tài)系統(tǒng)服務(wù)價值的影響研究也取得了一定成果[14-16]，但是針對稀土礦區(qū)這一特殊土地利用類型，研究其土地利用變化與生態(tài)系統(tǒng)服務(wù)價值的相關(guān)研究相對較少，尤其是在縣域小尺度范圍的相關(guān)研究還較為鮮見。事實上，比較分析稀土資源開采的前后，稀土礦區(qū)的土地利用發(fā)生了深刻的變化，研究它們的變化規(guī)律，并系統(tǒng)探討由此帶來的生態(tài)系統(tǒng)服務(wù)價值的變化，有利于更好地理解人類活動對礦區(qū)生態(tài)環(huán)境的影響，促進區(qū)域土地資源的合理利用及生態(tài)環(huán)境的改善，系統(tǒng)解決人地矛盾的突出問題，實現(xiàn)區(qū)域經(jīng)濟的持續(xù)健康發(fā)展。而基于縣域中小尺度下的分析，更能精準把握變化規(guī)律，在此基礎(chǔ)上提出的政策措施具有更強的針對性。

本研究以龍南市稀土礦區(qū)為研究對象，基于2009年、2019年2個不同時段龍南市土地調(diào)查數(shù)據(jù)和2017年龍南市稀土礦區(qū)分布范圍數(shù)據(jù)，利用GIS空間分析技術(shù)系統(tǒng)分析稀土礦區(qū)及其周邊土地利用變化特征規(guī)律，運用生態(tài)系統(tǒng)服務(wù)價值評價方法，探究稀土礦區(qū)土地利用變化對區(qū)域生態(tài)系統(tǒng)服務(wù)價值(ESV)的影響，以期為區(qū)域土地資源的高效合理利用、改善地區(qū)生態(tài)環(huán)境以及政府決策提供科學參考。

1 研究區(qū)概況

龍南市位于江西省贛州市南部，為亞熱帶季風氣候，年均溫19.2 ℃，年均降水1 506 mm，地形以山地丘陵為主，森林覆蓋率達80.30%。境內(nèi)有豐富的稀土礦產(chǎn)資源，其中重稀土的儲量和質(zhì)量都居世界之首，被譽為“重稀土之鄉(xiāng)”[17]。稀土礦區(qū)主要分布在龍南市東北部，距龍南市區(qū)約10 km，主要包括東江鄉(xiāng)足洞礦區(qū)、關(guān)西鎮(zhèn)富坑礦區(qū)、汶龍鎮(zhèn)黃沙礦區(qū)以及龍南鎮(zhèn)礦區(qū)等。2017年稀土礦區(qū)范圍總面積約393 hm2，礦區(qū)地理坐標位于東經(jīng)E114°46′～114°54′，北緯N24°48′～24°52′之間(圖1)。龍南稀土的開發(fā)歷史較為悠久，最早可以追溯到20世紀60年代末，20世紀80年代中期進入迅猛開采的發(fā)展階段，先后采取了露天池浸、堆浸、原地浸礦的開采工藝進行生產(chǎn)，多年開采造成了嚴重的生態(tài)環(huán)境問題[18]。

圖1 研究區(qū)地理位置圖

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

2.1 數(shù)據(jù)來源與預處理

本研究所使用的龍南市2009年全國第2次土地調(diào)查分類數(shù)據(jù)、2019年全國第3次土地調(diào)查分類數(shù)據(jù)及2017年龍南市稀土礦區(qū)分布范圍矢量數(shù)據(jù)均來源于龍南市政府相關(guān)部門。數(shù)據(jù)分類結(jié)果精度高，克服了以往常用的Landsat遙感影像分辨率低、解譯難度大造成的分類精度不高的缺陷。生態(tài)系統(tǒng)服務(wù)價值研究中用到的統(tǒng)計數(shù)據(jù)來源于《龍南市統(tǒng)計年鑒(2009—2019年)》。

參照2007年《土地利用現(xiàn)狀分類》國家標準[19]，根據(jù)研究區(qū)的土地利用情況和研究的現(xiàn)實需要，利用ArcGIS對龍南市第2、第3次全國土地調(diào)查劃分的土地利用類型進行重分類，將土地利用類型劃分為耕地、林地、草地、建設(shè)用地、采礦用地、水域、其他土地共7大類。

2.2 研究方法

2.2.1 土地利用變化特征分析 利用ArcGIS的裁剪工具，用龍南市2017年稀土礦區(qū)范圍對2009年、2019年龍南市第2、第3次全國土地調(diào)查重分類后的數(shù)據(jù)進行裁剪，生成稀土礦區(qū)2009年和2019年2期土地利用現(xiàn)狀數(shù)據(jù)，并通過土地利用類型面積統(tǒng)計表進行對比分析，探討龍南稀土礦區(qū)2009—2017年以及2017—2019年開采前后土地利用現(xiàn)狀及其變化特征。

利用ArcGIS的相交工具，對龍南稀土礦區(qū)2009年和2019年土地利用現(xiàn)狀圖進行疊加，基于土地利用轉(zhuǎn)移矩陣的方法對稀土礦區(qū)2009—2019年的10年來土地利用類型轉(zhuǎn)移的情況進行分析，探究稀土礦區(qū)土地利用轉(zhuǎn)變規(guī)律。轉(zhuǎn)移矩陣公式如下：

(1)

式中：S為面積；i為土地利用的類型數(shù)目；m、n為研究區(qū)初、末的土地利用類型；Smn為研究區(qū)初、末第m類土地轉(zhuǎn)化為第n類的面積[20]。

2.2.2 ArcGIS緩沖區(qū)分析 緩沖區(qū)分析是ArcGIS中用來確定不同地理要素的空間鄰近性和接近程度的空間操作。當對發(fā)生在地理要素及其附近的活動的影響范圍進行分析時，需要圍繞地理要素生成緩沖區(qū)，進行緩沖區(qū)分析[21]。利用ArcGIS緩沖區(qū)分析工具，以2017年龍南稀土礦區(qū)范圍為基礎(chǔ)，分別作出礦區(qū)外圍50 m、100 m、150 m、200 m范圍內(nèi)的緩沖區(qū)，與2009年和2019年土地利用數(shù)據(jù)疊加，得到礦區(qū)外圍不同緩沖區(qū)的土地利用現(xiàn)狀數(shù)據(jù)，來探討稀土礦區(qū)外圍一定范圍內(nèi)2009年和2019年的土地利用變化情況，以此分析礦區(qū)的稀土開采活動對周邊土地利用的影響。

2.2.3 生態(tài)系統(tǒng)服務(wù)價值評估 根據(jù)Costanza等[22]提出的生態(tài)系統(tǒng)服務(wù)價值(ESV)評估體系，參照謝高地等[23]修訂的“中國陸地生態(tài)系統(tǒng)單位面積生態(tài)系統(tǒng)服務(wù)價值當量表”確定研究區(qū)單位面積生態(tài)系統(tǒng)服務(wù)價值當量表(表1)。建設(shè)用地、采礦用地、其他土地不涉及相關(guān)生態(tài)功能的用0表示。隨后根據(jù)研究區(qū)土地利用結(jié)構(gòu)和研究區(qū)現(xiàn)實經(jīng)濟狀況，對單位面積農(nóng)田糧食年產(chǎn)量經(jīng)濟價值進行修正：研究區(qū)2009—2019年糧食均產(chǎn)為5 002.96 kg/hm2·a，以2019年糧食均價5.39元/kg為基準。考慮到?jīng)]有人力投入情況，根據(jù)謝高地提出的中國生態(tài)系統(tǒng)服務(wù)系數(shù)修正方法，即一個生態(tài)系統(tǒng)服務(wù)價值當量為單位面積糧食產(chǎn)值的1/7[2]，得到研究區(qū)食物生產(chǎn)經(jīng)濟價值為3 852.28元/hm2·a，獲得區(qū)域生態(tài)服務(wù)價值當量。根據(jù)區(qū)域生態(tài)服務(wù)價值當量得到研究區(qū)生態(tài)系統(tǒng)服務(wù)價值系數(shù)表(表2)。

表1 單位面積生態(tài)系統(tǒng)服務(wù)價值當量表

表2 研究區(qū)單位面積土地利用類型生態(tài)系統(tǒng)服務(wù)價值系數(shù)表/元·(hm2·a)-1

采用Costanza等提出的生態(tài)系統(tǒng)服務(wù)價值計算公式對龍南稀土礦區(qū)不同時期的生態(tài)系統(tǒng)服務(wù)價值進行計算，計算公式如下[24]：

ESV=∑Ak×VCk

(2)

ESVf=∑Ak×VCfk

(3)

式中：ESV為生態(tài)系統(tǒng)服務(wù)價值 (元)；Ak為土地利用類型k的總面積(hm2) ；VCk為土地利用類型k的單位面積生態(tài)系統(tǒng)服務(wù)價值系數(shù)(元/hm2·a )；ESVf為第f項生態(tài)系統(tǒng)服務(wù)價值；VCfk為土地利用類型k的第f項生態(tài)系統(tǒng)服務(wù)價值系數(shù)。

3 結(jié)果分析

3.1 稀土礦區(qū)土地利用變化特征

由龍南稀土礦區(qū)2009年、2019年土地利用現(xiàn)狀圖(圖2、圖3)和龍南稀土礦區(qū)2009年、2019年土地利用類型面積統(tǒng)計表(表3)可以看出，稀土礦區(qū)2009—2017以及2017—2019年開采前后土地利用現(xiàn)狀、結(jié)構(gòu)發(fā)生顯著變化。

圖2 龍南稀土礦區(qū)2009年土地利用現(xiàn)狀圖

圖3 龍南稀土礦區(qū)2019年土地利用現(xiàn)狀圖

表3 龍南稀土礦區(qū)2009年、2019年土地利用類型面積統(tǒng)計表

從礦區(qū)用地類型結(jié)構(gòu)的變化來看，2009—2017年，有111.3 hm2林地、44 hm2建設(shè)用地、11 hm2耕地轉(zhuǎn)化為稀土礦區(qū)，其中由林地轉(zhuǎn)變?yōu)橄⊥恋V區(qū)的面積最大，約占稀土礦區(qū)總面積的28.3%。未發(fā)生變化的采礦用地有213 hm2，占總面積的54.2%。2017—2019年，稀土礦區(qū)轉(zhuǎn)變?yōu)榱值氐拿娣e為122 hm2，轉(zhuǎn)變?yōu)榻ㄔO(shè)用地的面積為56 hm2，轉(zhuǎn)變?yōu)樗虻拿娣e為21 hm2，轉(zhuǎn)變?yōu)椴莸氐拿娣e為16 hm2，其中稀土礦區(qū)轉(zhuǎn)變?yōu)榱值氐拿娣e最大，約31%的稀土礦區(qū)轉(zhuǎn)變?yōu)榱值?。未發(fā)生變化的采礦用地有175 hm2，占總面積的44.5%。因此，2009—2017年前期階段龍南稀土礦區(qū)的稀土開采主要占用了林地、建設(shè)用地以及耕地，其中林地所占比重最大；2017—2019年后期階段龍南稀土礦區(qū)主要轉(zhuǎn)變?yōu)榱值?、建設(shè)用地、水域、草地，其中林地所占比重最大。

從礦區(qū)土地利用動態(tài)變化來看， 2009—2017年再到2019年變化最劇烈的是采礦用地，其所占比重由54.2%變?yōu)?00%再變?yōu)?4.5%；其次為林地，其所占比重由28.3%變?yōu)?再變?yōu)?1.0%。因此，采礦用地和林地是礦區(qū)土地利用變化最劇烈的兩種土地類型。采礦用地和林地發(fā)生劇烈變化的過程反映出龍南稀土礦區(qū)前期隨著稀土開采活動的深入和采礦用地面積的增加導致大片林地被毀，林地面積減少，由此產(chǎn)生了一系列生態(tài)環(huán)境問題，后期對稀土礦區(qū)進行了土地修復、土地復墾以及營造生態(tài)林的過程，使得林地增加，廢棄稀土礦區(qū)減少。

根據(jù)龍南稀土礦區(qū)2009—2019年土地利用轉(zhuǎn)移矩陣表(表4)可以看出10年來礦區(qū)不同土地類型的轉(zhuǎn)移變化方向。具體情況如下：2009—2019年，耕地由10.91 hm2減少到2.63 hm2，凈轉(zhuǎn)出8.28 hm2，主要轉(zhuǎn)變?yōu)榱值睾筒傻V用地；林地由111.25 hm2增加到121.15 hm2，凈轉(zhuǎn)入9.9 hm2，主要由建設(shè)用地 、采礦用地轉(zhuǎn)入；草地由4.79 hm2增加到16.09 hm2，凈轉(zhuǎn)入11.3 hm2，主要由林地、采礦用地轉(zhuǎn)入；建設(shè)用地由43.97 hm2增加到55.93 hm2，凈轉(zhuǎn)入11.96 hm2，主要由林地、采礦用地轉(zhuǎn)入；采礦用地由212.68 hm2減少到175.72 hm2，凈轉(zhuǎn)出36.96 hm2，主要轉(zhuǎn)變?yōu)榱值?、草地及建設(shè)用地；水域由8.24 hm2增加到21.05 hm2，凈轉(zhuǎn)入12.81 hm2，主要由林地、采礦用地轉(zhuǎn)入；其他土地面積較小，呈轉(zhuǎn)出狀態(tài)，轉(zhuǎn)出0.73 hm2。耕地、采礦用地、其他土地的轉(zhuǎn)出促使林地、草地、建設(shè)用地、水域的轉(zhuǎn)入，這有利于稀土礦區(qū)生態(tài)環(huán)境的改善。

表4 龍南稀土礦區(qū)2009—2019年土地利用轉(zhuǎn)移矩陣/hm2

3.2 稀土礦區(qū)外圍土地利用變化特征

根據(jù)龍南稀土礦區(qū)外圍緩沖區(qū)土地利用變化特征表(表5)可知，龍南稀土礦區(qū)外圍緩沖區(qū)土地利用類型主要是林地、采礦用地，其次為建設(shè)用地和耕地，有少量草地、水域和其他土地。具體來看，龍南稀土礦區(qū)外圍50 m、100 m、150 m、200 m緩沖區(qū)范圍內(nèi)的耕地面積從2009—2019年均呈減少的趨勢；10年間，稀土礦區(qū)外圍200 m范圍內(nèi)的耕地面積由66.8 hm2減少到31.4 hm2，減少了52.99%；從礦區(qū)外圍50～200 m范圍內(nèi)耕地減少的變化率越來越小。這表明礦區(qū)的稀土開采活動與礦區(qū)外圍的耕地減少有影響。距稀土礦區(qū)越近，耕地減少的趨勢越明顯；距稀土礦區(qū)越遠，耕地減少的趨勢越不明顯。因此，礦區(qū)的稀土開采是導致礦區(qū)外圍耕地面積減少的重要影響因素。緩沖區(qū)范圍內(nèi)林地面積從2009—2019年均呈增加的趨勢；10年間，稀土礦區(qū)外圍200 m范圍內(nèi)的林地面積由803.3 hm2增加到1 267.4 hm2，增加了57.77%；從50～200 m范圍內(nèi)林地增加的變化率越來越小。距稀土礦區(qū)越近，林地增加的趨勢越明顯，這反映出隨著礦區(qū)稀土開采的進行和環(huán)境問題的產(chǎn)生，礦區(qū)外圍實施了植樹造林、營造生態(tài)防護林等使林地面積增加。緩沖區(qū)范圍內(nèi)建設(shè)用地面積從2009—2019年均呈增加的趨勢；10年間，稀土礦區(qū)外圍200 m范圍內(nèi)的建設(shè)用地面積由61.8 hm2增加到95.9 hm2，增加了55.18%；從50～200 m范圍內(nèi)建設(shè)用地增加的變化率越來越大。這可能是因為距稀土礦區(qū)越遠的外圍地區(qū)，受稀土開采活動的負面影響越小，居民的生產(chǎn)建設(shè)活動越活越。緩沖區(qū)范圍內(nèi)采礦用地面積從2009—2019年間均呈減少的趨勢；10年間，稀土礦區(qū)外圍200 m范圍內(nèi)的采礦用地面積由6 62.9 hm2減少到198.4 hm2，減少了70.07%；從50～200 m范圍內(nèi)采礦用地減少的變化率越來越大。這反映出稀土礦區(qū)外圍進行了廢棄礦區(qū)和尾礦堆積區(qū)等的治理及復墾使得采礦用地面積減小。緩沖區(qū)范圍內(nèi)草地、水域的面積從2009—2019年均呈增加的趨勢，增加的面積較??；其他土地的面積從2009—2019年均呈減少的趨勢，減少的面積較小。

表5 龍南稀土礦區(qū)外圍緩沖區(qū)2009年和2019年土地利用變化特征表

年份2009年2019年2009年2019年土地類型150 m以內(nèi)/hm2變化率/%200 m以內(nèi)/hm2變化率/%耕地51.821.6-58.3066.831.4-52.99林地564.1938.866.42803.31267.457.77草地15.422.344.8123.524.85.53建設(shè)用地54.378.845.1261.895.955.18采礦用地554.6179.7-67.60662.9198.4-70.07水域10.216.460.7811.820.372.47其他土地1.40.5-64.291.50.5-66.67

對比發(fā)現(xiàn)，從2009—2019年的10年間，稀土礦區(qū)外圍200 m范圍內(nèi)采礦用地減少了464.5 hm2，而林地增加了464.1 hm2，采礦用地減少的面積與林地增加的面積實現(xiàn)了互補。由此可推測，礦區(qū)外圍采取了有效的生態(tài)恢復措施，使得2009—2019年來礦區(qū)外圍200 m范圍內(nèi)存在大量的采礦用地轉(zhuǎn)變?yōu)榱值亍?/p>

3.3 生態(tài)系統(tǒng)服務(wù)價值評價

稀土礦區(qū)土地利用變化引起區(qū)域生態(tài)系統(tǒng)價值服務(wù)的變化[25]。根據(jù)研究區(qū)土地利用類型單位面積生態(tài)系統(tǒng)服務(wù)價值(ESV)系數(shù)表(表2)和龍南稀土礦區(qū)2009年、2019年土地利用類型面積統(tǒng)計表(表3)計算得到龍南稀土礦區(qū)土地利用類型2009—2019年生態(tài)系統(tǒng)服務(wù)價值(表6)。

表6 龍南稀土礦區(qū)土地利用類型2009—2019年生態(tài)系統(tǒng)服務(wù)價值

由表6可以看出，龍南稀土礦區(qū)2009—2019年來生態(tài)系統(tǒng)服務(wù)價值總量呈增加的趨勢，由2009年的853.76萬元增加到2019年的1 699.39萬元，增加了845.63萬元，增長近1倍。具體來看，2009—2019年耕地面積的減少導致耕地的ESV呈減少的趨勢，由16.48萬元減少到4.50萬元，減少了72.75%。林地、草地、水域面積的增加導致其ESV呈增加的趨勢。林地的ESV由984.00萬元增加到1 078.60萬元，增加了9.61%；草地的ESV由36.41萬元增加到121.36萬元，增加了233.32%；水域的ESV由396.79萬元增加到1 016.16萬元，增加了156.10%。其他土地的面積及其ESV未發(fā)生明顯變化。而建設(shè)用地和采礦用地對生態(tài)系統(tǒng)服務(wù)價值產(chǎn)生負面影響。2009年建設(shè)用地的使用對ESV造成了99.33萬元的損失，2019年對ESV造成了126.42萬元的損失。因此，隨著建設(shè)用地面積的擴大，建設(shè)用地對ESV造成的損失也逐漸加大。2009年采礦用地的使用對ESV造成了480.83萬元的損失，2019年對ESV造成了395.05萬元的損失，由此可見，隨著礦區(qū)土地整治、土地復墾的進行，采礦用地面積減少，采礦用地對ESV造成的損失也逐漸減小。從各土地利用類型產(chǎn)生的生態(tài)系統(tǒng)服務(wù)價值量來看，林地和水域的ESV最大，其次為草地和耕地，因此，林地和水域?qū)^(qū)域生態(tài)系統(tǒng)服務(wù)價值的影響最大，是對區(qū)域生態(tài)系統(tǒng)服務(wù)價值總量增加貢獻最大的2種土地利用類型。

4 討論與分析

從龍南稀土礦區(qū)10年間土地利用變化及礦區(qū)外圍土地利用情況分析來看，早期礦區(qū)的稀土開采活動很大程度上存在著對山區(qū)土地的挖損和占用林地，導致林地轉(zhuǎn)變?yōu)橄⊥敛傻V用地。隨著稀土開采活動的深入，導致大量的耕地損毀、壓占，耕地轉(zhuǎn)變?yōu)椴傻V用地，耕地面積減少，產(chǎn)生了一系列生態(tài)環(huán)境問題。近幾年來，龍南市政府及環(huán)保部門高度重視廢棄稀土礦山環(huán)境的治理與恢復問題，積極探索廢棄稀土礦山治理的有效途徑，取得了顯著的成效。例如，龍南足洞稀土礦區(qū)在早期的稀土開采過程中存在無序化現(xiàn)象，使得土壤中的氨氮含量偏高，土壤酸性過高，還產(chǎn)生了水土流失、土地沙化、植被破壞等一系列問題。對此，足洞稀土礦區(qū)進行了土地復墾和土地修復，采用“山頂栽松，坡面布草，臺地種桑，溝谷植竹”的整體布局，利用尾砂種植百喜草、狗尾草等植物，目前，復墾區(qū)各種植物長勢良好，產(chǎn)生較好的經(jīng)濟效益的同時還改善了生態(tài)環(huán)境[26]。正是稀土礦區(qū)及其外圍采取了一系列行之有效的治理措施，使得采礦用地向林地、水域、草地等土地類型轉(zhuǎn)變。土地利用發(fā)生變化的同時也改善著礦區(qū)生態(tài)環(huán)境，使得礦區(qū)生態(tài)系統(tǒng)價值服務(wù)(ESV)隨之增加。

從總體上看，龍南稀土礦區(qū)土地利用動態(tài)變化并未使區(qū)域生態(tài)系統(tǒng)服務(wù)價值量(ESV)減少，ESV反而呈增加的趨勢，這主要是因為林地和水域?qū)^(qū)域ESV總量增加具有重要貢獻。2009—2019年礦區(qū)林地和水域面積增加產(chǎn)生的正面影響彌補了由建設(shè)用地、采礦用地侵占耕地以及耕地塌陷、耕地破壞對區(qū)域ESV產(chǎn)生的負面影響，使得礦區(qū)ESV呈現(xiàn)出收益大于損失的情況。未來龍南稀土礦區(qū)及其周邊在進行稀土開采活動的同時應(yīng)加強林地、水域等生態(tài)用地的保護，以緩解稀土開采對區(qū)域生態(tài)環(huán)境造成的壓力。

土地利用結(jié)構(gòu)變化對礦區(qū)生態(tài)系統(tǒng)服務(wù)價值起著至關(guān)重要的作用[27]，因此，要科學合理利用土地，發(fā)揮土地利用變化對生態(tài)系統(tǒng)服務(wù)價值的正向作用。但是，區(qū)域生態(tài)系統(tǒng)服務(wù)價值不僅會受土地利用變化的影響，還會受氣候變化、人口密度、經(jīng)濟發(fā)展水平、城市化率及綠化率等諸多自然、社會經(jīng)濟因素的影響。本文僅從土地利用結(jié)構(gòu)變化方面分析其對區(qū)域生態(tài)系統(tǒng)服務(wù)價值的影響，未來的研究更應(yīng)考慮上述自然因素及社會經(jīng)濟因素對區(qū)域生態(tài)系統(tǒng)服務(wù)價值的影響，以提高研究結(jié)果的科學性和準確性。

5 結(jié)論

本文系統(tǒng)分析了稀土礦區(qū)及其周邊土地利用變化規(guī)律，并對稀土礦區(qū)土地利用變化引起的生態(tài)系統(tǒng)服務(wù)價值的變化進行分析，結(jié)果表明如下。

1)2009—2017年前期階段龍南稀土礦區(qū)的稀土開采主要占用了林地、建設(shè)用地以及耕地；2017—2019年后期階段龍南稀土礦區(qū)主要轉(zhuǎn)變?yōu)榱值?、建設(shè)用地、水域、草地。其中采礦用地和林地變化最劇烈，其劇烈變化的過程反映了礦區(qū)前期生態(tài)破壞與后期生態(tài)修復的過程。2009—2019的十年間，龍南稀土礦區(qū)范圍內(nèi)的耕地、采礦用地、其他土地呈轉(zhuǎn)出趨勢，林地、草地、建設(shè)用地、水域呈轉(zhuǎn)入趨勢。

2)龍南稀土礦區(qū)外圍的土地利用類型主要是林地、采礦用地，其次為建設(shè)用地和耕地，有少量草地、水域和其他土地。2009—2019年來，礦區(qū)外圍的耕地和采礦用地面積均呈現(xiàn)減少的趨勢；林地和建設(shè)用地面積均呈增加的趨勢；其他的土地類型變化較小。采礦用地減少與林地增加的面積實現(xiàn)了互補。

3)2009—2019年來，龍南稀土礦區(qū)的生態(tài)系統(tǒng)服務(wù)價值總量呈增加的趨勢，由2009年的853.76萬元增加到2019年的1 699.39萬元，增長近1倍。其中林地和水域?qū)^(qū)域生態(tài)系統(tǒng)服務(wù)價值總量增加貢獻最大，林地和水域面積的增加彌補了建設(shè)用地增加和采礦用地占用耕地對區(qū)域生態(tài)系統(tǒng)服務(wù)價值產(chǎn)生的負面影響。