基于土地利用變化的臺州市生態(tài)系統(tǒng)服務(wù)價值研究

［關(guān)鍵詞］ 土地利用；生態(tài)系統(tǒng)；服務(wù)價值；臺州市

［摘 要］ 土地利用變化對生態(tài)系統(tǒng)服務(wù)價值影響的研究對土地合理利用與可持續(xù)發(fā)展具有重大意義，定量評估生態(tài)系統(tǒng)服務(wù)價值是人類合理利用生態(tài)系統(tǒng)的重要依據(jù)。臺州市居山面海，平原丘陵相間，是浙江“七山一水兩分田”的縮影。以臺州市2010年、2015年、2018年三期土地利用數(shù)據(jù)為基礎(chǔ)，構(gòu)建生態(tài)系統(tǒng)服務(wù)價值當(dāng)量因子表，分析臺州市土地利用變化與生態(tài)系統(tǒng)服務(wù)價值的響應(yīng)關(guān)系。結(jié)果表明：研究時段內(nèi)林地是臺州市的主要用地類型，土地利用變化表現(xiàn)為園地、林地、水域轉(zhuǎn)為其他用地，草地、耕地相對穩(wěn)定；生態(tài)系統(tǒng)服務(wù)價值總體呈下降趨勢，林地是生態(tài)系統(tǒng)服務(wù)價值貢獻(xiàn)率最高的土地類型，其次為水域；生態(tài)經(jīng)濟(jì)協(xié)調(diào)度分析表明，研究區(qū)經(jīng)濟(jì)快速發(fā)展已影響生態(tài)系統(tǒng)穩(wěn)定，進(jìn)一步加強(qiáng)生態(tài)文明建設(shè)是區(qū)域可持續(xù)發(fā)展的保障。

［中圖分類號］ X171.1；F301.24" [文獻(xiàn)標(biāo)識碼] A" DOI：10.3969/j.issn.1000-0941.2024.12.014

［引用格式］ 張陽，陳星，陳小民，等.基于土地利用變化的臺州市生態(tài)系統(tǒng)服務(wù)價值研究[J].中國水土保持，2024（12）：48-52.

土地利用變化對生態(tài)系統(tǒng)服務(wù)價值有重要影響^[1]。20世紀(jì)以來，隨著人口增多和城市化快速發(fā)展，人類活動引發(fā)了土地覆被變化，而土地覆被變化對地表景觀格局、土地功能以及未來人類活動等造成影響^[2-3]。土地利用類型變更導(dǎo)致的全球生態(tài)系統(tǒng)服務(wù)價值下降已成為影響人類生存和可持續(xù)發(fā)展的重大問題^[4-6]。1997年，COSTANZA et al.^[7]首先提出了全球生態(tài)服務(wù)價值評價方法，吸引了眾多學(xué)者開展相關(guān)研究。其中，謝高地等^[8]通過調(diào)查研究提出了被中國學(xué)者普遍接受的中國陸地系統(tǒng)生態(tài)服務(wù)價值當(dāng)量表。生態(tài)系統(tǒng)服務(wù)價值的研究已逐漸成為國內(nèi)研究熱點(diǎn)^[9-11]。臺州市居山面海，平原丘陵相間，是浙江省“七山一水兩分田”的縮影，近年來經(jīng)濟(jì)建設(shè)造成的土地利用類型變化，進(jìn)而引起生態(tài)系統(tǒng)服務(wù)價值變化，然而目前相關(guān)研究成果匱乏。因此，本研究基于2010年、2015年、2018年三期臺州市土地利用數(shù)據(jù)，運(yùn)用當(dāng)量因子評估法評估臺州市生態(tài)系統(tǒng)服務(wù)價值變化，以期為臺州市可持續(xù)發(fā)展提供科學(xué)依據(jù)。

1 研究區(qū)概況

臺州市位于浙江省東部溫黃平原北部，三面環(huán)山、一面臨海，全市陸地總面積約100.50萬hm²，其中山地、丘陵面積占全市陸地總面積的70.4%；海岸線長約941 km，島陸域面積約2.74萬hm²。屬亞熱帶季風(fēng)氣候區(qū)，四季分明，雨水充沛，不同地點(diǎn)年降水量1 002～2 135 mm，年均氣溫18.2 ℃，在中國植被區(qū)劃和浙江省植被區(qū)劃中屬中亞熱帶常綠闊葉林帶，浙閩山丘甜櫧、木荷林區(qū)，地帶性植被為常綠闊葉林。

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

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

本研究采用的臺州市2010年、2015年、2018年的土地利用數(shù)據(jù)，來源于地理空間數(shù)據(jù)云；GDP、糧食產(chǎn)量、糧食單價等數(shù)據(jù)來源于臺州市統(tǒng)計(jì)年鑒和政府年度工作報(bào)告。

2.2 研究方法

2.2.1 土地利用變化

本研究采用單一土地利用動態(tài)度表示研究區(qū)一定時間范圍內(nèi)某種土地利用類型向另一種土地利用類型的變化情況^[12]；采用綜合土地利用動態(tài)度表示研究區(qū)一定時間范圍內(nèi)某種土地利用類型向其他土地利用類型的變化情況^[13]，值越低，代表該土地利用類型越穩(wěn)定，值越高則相反。

2.2.2 生態(tài)系統(tǒng)服務(wù)價值評估方法

以謝高地等^[14]改進(jìn)的中國陸地生態(tài)系統(tǒng)單位面積生態(tài)系統(tǒng)服務(wù)價值當(dāng)量表為基礎(chǔ)，基于臺州市實(shí)際情況，優(yōu)化提出臺州市不同土地利用類型單位面積生態(tài)系統(tǒng)服務(wù)價值當(dāng)量（見表1）。

一個標(biāo)準(zhǔn)單位生態(tài)系統(tǒng)服務(wù)價值當(dāng)量因子的價值量可看作單位面積農(nóng)田食物生產(chǎn)服務(wù)的價值量，此時價值當(dāng)量為1.0。單位面積農(nóng)田食物生產(chǎn)服務(wù)的價值量相當(dāng)于1 hm²研究區(qū)平均產(chǎn)量的農(nóng)田每年自然糧食產(chǎn)量的經(jīng)濟(jì)價值，計(jì)算公式為^[15]

E a=17（∑nⁱ⁼¹m ip iq i）1S（1）

式中：E ^a為研究區(qū)單位面積農(nóng)田生態(tài)系統(tǒng)食物生產(chǎn)服務(wù)的價值量，單位元/hm²；n為研究區(qū)主要糧食作物（稻谷、小麥、玉米等）種類數(shù)量；m ⁱ為農(nóng)作物i的播種面積，單位hm²；p ⁱ為單位面積農(nóng)田的農(nóng)作物i的產(chǎn)量，單位kg/hm²；q ⁱ為農(nóng)作物i的價格，單位元/kg；S為n種糧食作物總播種面積，單位hm²；1/7指沒有人力投入的自然生態(tài)系統(tǒng)提供的經(jīng)濟(jì)價值占現(xiàn)有單位面積農(nóng)田生態(tài)系統(tǒng)食物生產(chǎn)服務(wù)經(jīng)濟(jì)價值的比例。

本研究以2020年臺州市糧食產(chǎn)量和單價為準(zhǔn)，計(jì)算得到2020年臺州市標(biāo)準(zhǔn)當(dāng)量價值量為4 351.44元/hm²。然后利用價值當(dāng)量換算公式^[16]，計(jì)算得到臺州市生態(tài)系統(tǒng)服務(wù)價值和具體每項(xiàng)生態(tài)系統(tǒng)服務(wù)功能價值。

2.2.3 生態(tài)經(jīng)濟(jì)協(xié)調(diào)度

在研究區(qū)域生態(tài)系統(tǒng)與經(jīng)濟(jì)發(fā)展的協(xié)調(diào)程度時通常采用生態(tài)經(jīng)濟(jì)協(xié)調(diào)度為評價指標(biāo)，目前生態(tài)經(jīng)濟(jì)協(xié)調(diào)度評價方法暫無統(tǒng)一標(biāo)準(zhǔn)。顧劍紅等^[17]在進(jìn)行廣西縣級生態(tài)服務(wù)價值研究時通過對多種方法的篩選，確定了一種能全面反映地區(qū)生態(tài)與經(jīng)濟(jì)發(fā)展關(guān)系的計(jì)算方法，并結(jié)合魏曉旭等^[18]提出的生態(tài)經(jīng)濟(jì)協(xié)調(diào)度等級劃分標(biāo)準(zhǔn)，計(jì)算公式為

E=F r/G r（2）

式中：E為生態(tài)經(jīng)濟(jì)協(xié)調(diào)度；F r為研究時段內(nèi)單位面積生態(tài)系統(tǒng)服務(wù)價值變化率，%；G r為研究時段內(nèi)單位面積GDP變化率，%。

由于目前GDP處于增長階段，因此生態(tài)經(jīng)濟(jì)協(xié)調(diào)度分為3種情況：①E≥1時，表征生態(tài)系統(tǒng)服務(wù)價值增長速度大于或等于經(jīng)濟(jì)增長速度；②0≤E＜1時，表征生態(tài)系統(tǒng)服務(wù)價值增長速度小于經(jīng)濟(jì)增長速度；③E＜0時，表征經(jīng)濟(jì)發(fā)展阻礙了生態(tài)系統(tǒng)的可持續(xù)發(fā)展。

3 結(jié)果與分析

3.1 土地利用類型變化

臺州市2010年、2015年、2018年各土地利用類型面積見表2。在所有的土地利用類型中，林地所占的比例最大，最高為2010年的50.45％，最低為2018年的的49.96％，其次是耕地、其他用地、水域和園地，草地所占比例最小，各年度均低于1％。由于灘涂淤漲，因此2010—2018年臺州市土地總面積呈上升趨勢，從2010年的1 003 791.44 hm²上升到2018年的^{1 005 042.60} hm²；園地、林地、草地、水域總體呈下降趨勢，林地由506 291.46 hm²下降到502 042.30 hm²，園地由61 866.48 hm²下降到56 144.45 hm²，草地由^{8 119.25} hm2下降到7 845.12 hm2，水域由108 648.19 hm2下降到102 749.00 hm2；其他用地持續(xù)增加，從2010年的126 697.30 hm²上升到2018年的^{143 356.76} hm²，而耕地相對穩(wěn)定，變化較小。

3.2 土地利用動態(tài)度分析

2010—2015年、2015—2018年臺州市單一土地利用動態(tài)度和綜合土地利用動態(tài)度計(jì)算結(jié)果見表3。其他用地是研究期內(nèi)轉(zhuǎn)變最活躍的土地類型，在2010—2015年、2015—2018年其單一動態(tài)度分別為1.51%、^1.73%；其次是園地，在2010—2015年、2015—2018年單一動態(tài)度分別為-0.91%、-1.64%，說明隨著時間的推移，其他用地、園地的穩(wěn)定狀況越來越差；耕地是所有用地類型中轉(zhuǎn)移最為穩(wěn)定的一類，在2010—2015年、2015—2018年其單一動態(tài)度分別為0.07%、^0.01%，說明耕地面積只有小幅上升；其余用地類型的單一動態(tài)度絕對值從高到低依次為水域、草地、林地。綜合土地利用動態(tài)度反映了研究時段內(nèi)區(qū)域土地利用變化的強(qiáng)度，從表3可以看出臺州市土地利用在2015—2018年的變化強(qiáng)度略高于2010—2015年。

3.3 生態(tài)服務(wù)價值動態(tài)演變

根據(jù)臺州市土地利用類型面積和生態(tài)系統(tǒng)服務(wù)價值當(dāng)量因子表，得到臺州市2010—2018年各土地利用類型生態(tài)系統(tǒng)服務(wù)價值（見表4）。生態(tài)服務(wù)價值總體呈下降趨勢，2010—2015年下降了152 531.30萬元，2015—2018年下降了88 036.05萬元。園地和水域的生態(tài)系統(tǒng)服務(wù)價值所占比例穩(wěn)定減小，園地在2010年、2015年、2018年所占比例分別為3.86%、3.75%、3.61%，水域分別為38.96%、38.27%、37.95%；耕地和林地的生態(tài)系統(tǒng)服務(wù)價值所占比例呈穩(wěn)步增加趨勢，耕地從2010年的4.03%上升到2018年的4.17%，林地從2010年的52.63%上升到2018年的53.76%；草地的生態(tài)系統(tǒng)服務(wù)價值所占比例較小，在整個時段內(nèi)變化不明顯。

為深入分析，根據(jù)臺州市單項(xiàng)服務(wù)生態(tài)價值當(dāng)量因子表與土地利用情況計(jì)算得到臺州市單項(xiàng)生態(tài)系統(tǒng)服務(wù)價值（見表5）。在一級類型中，調(diào)節(jié)服務(wù)所占比例最高，其次是支持服務(wù)和供給服務(wù)，文化服務(wù)比例最低。一級服務(wù)類型中，供給服務(wù)和調(diào)節(jié)服務(wù)所占比例穩(wěn)定減小，其中調(diào)節(jié)服務(wù)2010年、2015年、2018年所占比例分別為74.84%、74.72%、74.64%，支持服務(wù)和文化服務(wù)所占比例呈穩(wěn)步增加趨勢，其中支持服務(wù)從2010年的^15.25%上升到2018年的15.45%。調(diào)節(jié)服務(wù)包括氣體調(diào)節(jié)、氣候調(diào)節(jié)、凈化環(huán)境、水文調(diào)節(jié)4類二級類型，它們的價值比例約為1∶3∶1∶7；支持服務(wù)包括土壤保持、維持養(yǎng)分循環(huán)和生物多樣性，比例約為10∶1∶10。二級服務(wù)在地區(qū)生態(tài)系統(tǒng)服務(wù)總價值中所占比例從高到低依次為水文調(diào)節(jié)、氣候調(diào)節(jié)、土壤保持、生物多樣性、氣體調(diào)節(jié)、凈化環(huán)境、美學(xué)景觀、水資源供給、食物生產(chǎn)、原料生產(chǎn)、維持養(yǎng)分循環(huán)，其中水資源供給、水文調(diào)節(jié)所占比例總體呈下降趨勢，其余各項(xiàng)均呈增加趨勢。

3.4 生態(tài)經(jīng)濟(jì)協(xié)調(diào)度分析

基于區(qū)域生態(tài)系統(tǒng)服務(wù)價值和GDP數(shù)據(jù)，計(jì)算得到2010—2015年、2015—2018年臺州市生態(tài)經(jīng)濟(jì)協(xié)調(diào)度。2010—2018年臺州市GDP一直處于快速增長階段，然而生態(tài)系統(tǒng)服務(wù)價值呈下降趨勢，計(jì)算得到生態(tài)經(jīng)濟(jì)協(xié)調(diào)度在2010—2015年為-0.15，在2015—2018年為-0.16，說明經(jīng)濟(jì)的快速發(fā)展已阻礙了生態(tài)系統(tǒng)的可持續(xù)發(fā)展。結(jié)合土地利用數(shù)據(jù)和不同土地利用類型的生態(tài)服務(wù)價值數(shù)據(jù)，2010—2018年經(jīng)濟(jì)的快速發(fā)展導(dǎo)致林地、園地、草地、水域被侵占，其他土地面積快速增長，生態(tài)系統(tǒng)服務(wù)價值下降，經(jīng)濟(jì)發(fā)展已對生態(tài)系統(tǒng)穩(wěn)定性產(chǎn)生了不利影響。

4 結(jié)論與討論

基于臺州市2010年、2015年、2018年三期土地利用數(shù)據(jù)分析了生態(tài)系統(tǒng)服務(wù)價值動態(tài)演變規(guī)律，結(jié)果表明，2010—2018年臺州市土地總面積增加1 251.16 hm²（主要為沿海各縣灘涂圍墾新增），其中園地、林地、水域、草地面積總體呈下降趨勢，其他用地面積持續(xù)增加，耕地面積相對穩(wěn)定。林地作為面積最大的用地類型，同時也是單位面積生態(tài)系統(tǒng)服務(wù)價值最高的地類，因社會經(jīng)濟(jì)發(fā)展，多年來大量林地被占用，林地占補(bǔ)平衡未得到有效實(shí)施，林地面積日趨減少，相較而言耕地因“18億畝（1 hm2=15畝）紅線”政策，總量保持穩(wěn)定；此外林業(yè)產(chǎn)業(yè)發(fā)展慢，收益低，林農(nóng)毀林改種的情況也時有發(fā)生。隨著社會發(fā)展邁入新階段，林地在生態(tài)文明建設(shè)中的地位日益突顯，為有效保障林地體量，一方面需要政府嚴(yán)格林地用地制度，推行占補(bǔ)平衡，另一方面要改革創(chuàng)新林業(yè)產(chǎn)業(yè)，提高林農(nóng)收益，調(diào)動全社會保護(hù)林地的積極性。

土地利用動態(tài)度分析表明，其他用地是研究區(qū)內(nèi)轉(zhuǎn)化最活躍的土地利用類型，這主要是因?yàn)榻?jīng)濟(jì)發(fā)展帶來人口增加和城區(qū)擴(kuò)張。其他用地主要為城鄉(xiāng)居民建設(shè)用地等地表硬化的用地類型，生態(tài)系統(tǒng)服務(wù)價值較低，為提升生態(tài)效益的同時保障社會經(jīng)濟(jì)發(fā)展客觀需求，相關(guān)部門應(yīng)提升城鄉(xiāng)建設(shè)規(guī)劃能力，實(shí)現(xiàn)用地的集約化、高效化，同時嚴(yán)格控制違規(guī)、違章建設(shè)，對生態(tài)紅線、自然保護(hù)地等高生態(tài)價值區(qū)域加強(qiáng)用地審批，實(shí)現(xiàn)經(jīng)濟(jì)社會可持續(xù)發(fā)展。

一級生態(tài)系統(tǒng)服務(wù)價值排序?yàn)檎{(diào)節(jié)服務(wù)＞支持服務(wù)＞供給服務(wù)＞文化服務(wù)，其中調(diào)節(jié)服務(wù)約占總生態(tài)系統(tǒng)服務(wù)價值的75％。2010—2018年，供給服務(wù)和調(diào)節(jié)服務(wù)占比呈下降趨勢，支持服務(wù)和文化服務(wù)占比呈上升趨勢，其原因一方面是生態(tài)系統(tǒng)單項(xiàng)服務(wù)價值與用地類型密切相關(guān)，林地、水域的單位面積供給服務(wù)和調(diào)節(jié)服務(wù)價值相對其他地類較高，兩者面積的減小必然導(dǎo)致供給服務(wù)和調(diào)節(jié)服務(wù)價值整體下降；另一方面相較于支持服務(wù)包含的生物多樣性和文化服務(wù)包含的美學(xué)景觀，供給服務(wù)包含的食物生產(chǎn)、原料生產(chǎn)和調(diào)節(jié)服務(wù)包含的氣體調(diào)節(jié)、氣候調(diào)節(jié)、水文調(diào)節(jié)等，因總量大，大眾直觀感受弱、敏感度低，更容易被忽視，需要政府加強(qiáng)數(shù)據(jù)監(jiān)測與社會宣傳，通過引導(dǎo)生態(tài)服務(wù)價值貨幣化等形式促進(jìn)全社會的關(guān)注和參與。

生態(tài)經(jīng)濟(jì)協(xié)調(diào)度分析表明，2010—2018年臺州市GDP快速增長，與之相對的是生態(tài)系統(tǒng)服務(wù)價值呈下降趨勢，生態(tài)經(jīng)濟(jì)協(xié)調(diào)度值始終小于0，說明經(jīng)濟(jì)發(fā)展已阻礙了生態(tài)系統(tǒng)的可持續(xù)發(fā)展。經(jīng)濟(jì)快速發(fā)展一方面提高了人民物質(zhì)生活水平，另一方面也對生態(tài)環(huán)境造成了巨大的負(fù)擔(dān)。近年來，加強(qiáng)生態(tài)文明建設(shè)已成為社會共識，為保證生態(tài)系統(tǒng)與經(jīng)濟(jì)的協(xié)調(diào)發(fā)展，需要進(jìn)一步壓緊生態(tài)紅線，統(tǒng)籌生態(tài)與經(jīng)濟(jì)的均衡發(fā)展，使生態(tài)系統(tǒng)能夠滿足快速發(fā)展的社會的需求。

本研究未考慮其他用地（主要為城鄉(xiāng)居民建設(shè)用地）的生態(tài)系統(tǒng)服務(wù)價值，對其價值標(biāo)準(zhǔn)的確定也有待探索。同時社會經(jīng)濟(jì)的發(fā)展將不可避免地占用生態(tài)用地，提高單位面積生態(tài)用地的生態(tài)服務(wù)價值至關(guān)重要，本研究中對相同用地類型的生態(tài)服務(wù)價值采用了統(tǒng)一的評價標(biāo)準(zhǔn)，未考慮各林分、各耕地類型之間的差異等，有待于在下一步的研究中予以完善。

［參考文獻(xiàn)]

［1] 王友生，余新曉，賀康寧，等.基于土地利用變化的懷柔水庫流域生態(tài)服務(wù)價值研究[J].農(nóng)業(yè)工程學(xué)報(bào)，2012，28（5）：246-251.

［2] 吳莉，侯西勇，徐新良，等.山東沿海地區(qū)土地利用和景觀格局變化[J].農(nóng)業(yè)工程學(xué)報(bào)，201 29（5）：207-216.

［3] 趙興國，潘玉君，丁生，等.中國省域土地利用可持續(xù)性評價及時空格局特征[J].農(nóng)業(yè)工程學(xué)報(bào)，2014，30（3）：^196-204.

［4] 鄭江坤，余新曉，賈國棟，等.密云水庫集水區(qū)基于LUCC的生態(tài)服務(wù)價值動態(tài)演變[J].農(nóng)業(yè)工程學(xué)報(bào)，2010，26（9）：315-320.

［5] 王友生，余新曉，賀康寧，等.基于CA-Markov模型的藉河流域土地利用變化動態(tài)模擬[J].農(nóng)業(yè)工程學(xué)報(bào)，201 7（12）：330-336.

［6] 趙陽，余新曉，賈劍波，等.紅門川流域土地利用景觀動態(tài)演變及驅(qū)動力分析[J].農(nóng)業(yè)工程學(xué)報(bào)，201 29（9）：^239-248.

［7] COSTANZA R， D′ARGE R， GROOT R D， et al. The value of the world′s ecosystem services and natural capital[J].World Environment，1997，25（1）：3-15.

［8] 謝高地，魯春霞，成升魁.全球生態(tài)系統(tǒng)服務(wù)價值評估研究進(jìn)展[J].資源科學(xué)，200 3（6）：5-9.

［9] 卓靜，朱延年，李韜，等.陜北黃土高原區(qū)生態(tài)服務(wù)價值對土地利用變化的響應(yīng)[J].水土保持研究，2022，29（5）：^336-342.

［10] 韓偉，徐珊.松嫩平原土地利用變化特征及生態(tài)系統(tǒng)服務(wù)價值研究：以哈爾濱市為例[J].中國農(nóng)學(xué)通報(bào)，2022，38（26）：82-90.

［11] 張莉金，白羽萍，胡業(yè)翠，等.不同SSP-RCP情景下中國生態(tài)系統(tǒng)服務(wù)價值評估[J].生態(tài)學(xué)報(bào)，202 43（2）：^510-521.

［12] 劉紀(jì)遠(yuǎn)，劉明亮，莊大方，等.中國近期土地利用變化的空間格局分析[J].中國科學(xué)（D輯：地球科學(xué)），2002，32（12）：1031-1040.

［13] 王秀蘭，包玉海.土地利用動態(tài)變化研究方法探討[J].地理科學(xué)進(jìn)展，1999，18（1）：81-87.

［14] 謝高地，張彩霞，張雷明，等.基于單位面積價值當(dāng)量因子的生態(tài)系統(tǒng)服務(wù)價值化方法改進(jìn)[J].自然資源學(xué)報(bào)，2015，30（8）：1243-1254.

［15] 王永琪，馬姜明.基于縣域尺度珠江-西江經(jīng)濟(jì)帶廣西段土地利用變化對生態(tài)系統(tǒng)服務(wù)價值的影響研究[J].生態(tài)學(xué)報(bào)，2020，40（21）：7826-7839.

［16] 謝高地，甄霖，魯春霞，等.一個基于專家知識的生態(tài)系統(tǒng)服務(wù)價值化方法[J].自然資源學(xué)報(bào)，2008，23（5）：^911-919.

［17] 顧劍紅，王玉杰，王云琦，等.土地利用變化對生態(tài)服務(wù)價值的影響：以廣西隆林縣和西林縣為例[J].中國水土保持科學(xué)，2016，14（5）：100-109.

［18] 魏曉旭，趙軍，魏偉，等.基于縣域單元的中國生態(tài)經(jīng)濟(jì)系統(tǒng)協(xié)調(diào)度及空間演化[J].地理科學(xué)進(jìn)展，2014，33（11）：1535-1545.

收稿日期： 2023-08-15