蔡振饒,方朝陽*,何清華,顏吉林,高 丹,劉志勇
1. 江西師范大學(xué),江西 南昌 330022
2. 上饒師范學(xué)院,江西 上饒 334001
城鄉(xiāng)梯度是研究人類活動(dòng)干擾,尤其是城市化對(duì)生態(tài)系統(tǒng)影響的有力工具[1]. 隨著城市化的持續(xù)推進(jìn),人類活動(dòng)給生態(tài)環(huán)境帶來的影響不斷加劇,描繪城市化特征的城鄉(xiāng)梯度方法逐漸成為領(lǐng)域關(guān)注熱點(diǎn),發(fā)展為生態(tài)學(xué)研究的新框架. 該框架假設(shè)人類活動(dòng)與生態(tài)系統(tǒng)存在沿著城市核心向外圍農(nóng)村發(fā)生梯度變化的規(guī)律,包括土地利用和土地覆蓋的變化[2]、土壤化學(xué)[3]、生物多樣性[4]、景觀格局及生態(tài)系統(tǒng)服務(wù)的變化[5].
學(xué)者對(duì)城鄉(xiāng)梯度的劃分有不同定義和手段. 基于簡單的城鄉(xiāng)二分法不能理解城市化下生態(tài)過程的復(fù)雜性,因此,一些研究使用分類測(cè)量來識(shí)別城鄉(xiāng)梯度,梯度通常被定義為3或4個(gè)類別(如城市、近郊、遠(yuǎn)郊). 在生態(tài)學(xué)研究中主要以3種方式來描述和分類城鄉(xiāng)梯度. 第一種是基于地理距離方式,從城市中心到農(nóng)村外圍邊界之間建立一條或多條直線[6],或設(shè)置多個(gè)緩沖帶[7-8],以距離作為唯一的劃分指標(biāo). 第二種是以某種單一規(guī)則來生成城鄉(xiāng)梯度類型,試圖將城市化過程封裝成單一指標(biāo),如不透水表面比例[9]、森林覆蓋率[10]、人口或住房密度[11]等. 然而,這兩種分類方法應(yīng)用單一的度量閾值,掩蓋了區(qū)域的生態(tài)多樣性,忽略了復(fù)雜人類-生態(tài)系統(tǒng)的非線性特征[12]. 如有研究[13]發(fā)現(xiàn),物種豐富度往往在中等城市化的梯度中達(dá)到頂峰,在末端(城市核心和農(nóng)村)最低,在中部(郊區(qū))最高. 第三種方法則涉及復(fù)雜的景觀模式指標(biāo)(景觀斑塊的大小及多樣性)或多個(gè)指標(biāo)的組合[14-15](如不透水面比例與人口密度的組合). 目前國外已有研究使用聚類分析來生成城鄉(xiāng)梯度類型[16],通過聚類算法對(duì)多個(gè)指標(biāo)進(jìn)行組合,生成的類型可以適應(yīng)小尺度(鄉(xiāng)鎮(zhèn)級(jí)、村級(jí))的數(shù)據(jù)集,提供適應(yīng)小區(qū)域特定管理目標(biāo)的解決方案.
城市生態(tài)系統(tǒng)具有高度空間異質(zhì)性,采用城鄉(xiāng)梯度的研究范式可為生態(tài)系統(tǒng)和城市規(guī)劃管理提供重要的方法學(xué)基礎(chǔ)[17],為與社會(huì)福利相關(guān)的生態(tài)系統(tǒng)服務(wù)提供新的研究框架[18]. 已有研究表明,城鄉(xiāng)梯度制約著生態(tài)系統(tǒng)服務(wù)的空間變化[19],也影響著生態(tài)系統(tǒng)服務(wù)之間的關(guān)系,梯度的劃分方式會(huì)對(duì)生態(tài)系統(tǒng)服務(wù)空間異質(zhì)性研究結(jié)果產(chǎn)生影響. 目前國內(nèi)外從城鄉(xiāng)梯度視角探討生態(tài)系統(tǒng)服務(wù)空間異質(zhì)性的研究較少,大多采用構(gòu)建同心圓緩沖區(qū)的方式探索生態(tài)系統(tǒng)服務(wù)的城鄉(xiāng)梯度[20]. 然而,離城市中心距離相同的區(qū)域,其景觀與生態(tài)系統(tǒng)服務(wù)類型可能不同[21],城鄉(xiāng)梯度的表征不同. 因此,有必要應(yīng)用融合多種景觀指標(biāo)的城鄉(xiāng)梯度識(shí)別方法,深化與生態(tài)系統(tǒng)服務(wù)的關(guān)聯(lián)研究.
表 1 城鄉(xiāng)梯度的量化指標(biāo)Table 1 Quantitative index of the urban-rural gradient
南昌市擁有良好的自然底蘊(yùn),是江西省生態(tài)文明先行示范建設(shè)區(qū),同時(shí)是長江中游城市群的核心城市. 南昌市主城區(qū)是該市經(jīng)濟(jì)建設(shè)與生態(tài)空間保護(hù)的核心區(qū),厘清主城區(qū)的城鄉(xiāng)梯度特征及生態(tài)系統(tǒng)服務(wù)的空間異質(zhì)性,對(duì)于挖掘城市生態(tài)資源優(yōu)勢(shì),促進(jìn)鄉(xiāng)村生態(tài)價(jià)值轉(zhuǎn)化具有重要意義. 筆者以南昌市主城區(qū)為例,采用基于景觀聚類的方法識(shí)別主城區(qū)城鄉(xiāng)梯度,分析4種生態(tài)系統(tǒng)服務(wù)的城鄉(xiāng)梯度特征,使用841個(gè)村級(jí)尺度單元,符合對(duì)南昌市主城區(qū)分析的尺度要求. 該研究結(jié)果有利于南昌市主城區(qū)資源管理和生態(tài)環(huán)境保護(hù),可為城市生態(tài)系統(tǒng)服務(wù)分類保護(hù)與開發(fā)提供參考.
南昌市是江西省省會(huì),下轄6個(gè)區(qū)和3個(gè)縣,其中南昌縣已與6個(gè)市轄區(qū)同城化,自然與人文空間連成一體,成為名副其實(shí)的主城區(qū). 南昌市主城區(qū)是空間規(guī)模較大的城鄉(xiāng)共同體,同時(shí)內(nèi)部的城鄉(xiāng)地理環(huán)境差異較大. 從城市的視角看,2019年主城區(qū)建設(shè)用地面積達(dá)355.67 km2,已具有較大規(guī)模的城市體量. 在鄉(xiāng)村方面,主城區(qū)外圍具有大農(nóng)村的特征,城區(qū)北、東、南面屬于贛撫平原灌區(qū),農(nóng)田廣布,灌溉農(nóng)業(yè)發(fā)達(dá). 從自然環(huán)境的視角來看,研究區(qū)擁有“山水林田湖”等多樣化的自然要素,其北臨鄱陽湖、東臨青嵐湖,城中有瑤湖、贛江等眾多湖泊與河流,西部為梅嶺山地,林木茂盛(見圖1). 地理環(huán)境的多樣化特征為城鄉(xiāng)梯度研究創(chuàng)造了條件.
圖 1 南昌市主城區(qū)區(qū)位與土地覆被類型Fig.1 Location and land cover types of the main urban area of Nanchang
該研究選擇831個(gè)村級(jí)尺度行政單元進(jìn)行分析.南昌市青山湖區(qū)、青云譜區(qū)、東湖區(qū)、西湖區(qū)以社區(qū)為主,將其轄區(qū)內(nèi)的社區(qū)合并分別作為村級(jí)尺度單元. 總體上,村級(jí)尺度數(shù)據(jù)以行政村為主,還包括一部分林場(chǎng)、農(nóng)場(chǎng)、水產(chǎn)場(chǎng)、果場(chǎng)、湖泊管理局以及南昌縣政府飛地. 為了保證研究區(qū)完整性及研究結(jié)果的可靠性,允許這些具有特殊價(jià)值的飛地存在. 主城區(qū)中面積最大的村級(jí)單元為南磯鄉(xiāng)鄱陽湖區(qū),面積為323 km2,村級(jí)尺度單元的平均面積為5.5 km2.
土地覆被數(shù)據(jù)來源于清華大學(xué)全球高分辨率土地覆被監(jiān)測(cè)網(wǎng)(http://data.ess.tsinghua.edu.cn),數(shù)據(jù)年份為2017年. 村級(jí)行政邊界數(shù)據(jù)來源于第三次全國國土調(diào)查結(jié)果. 南昌市主城區(qū)農(nóng)、林、牧、漁總產(chǎn)值數(shù)據(jù)來源于《南昌市統(tǒng)計(jì)年鑒》(2018年). 景點(diǎn)、服務(wù)設(shè)施的POI(興趣點(diǎn))數(shù)據(jù)從百度地圖網(wǎng)站中獲得,各街道(鄉(xiāng)鎮(zhèn))人口數(shù)據(jù)來源于南昌市第六次人口普查.
在生態(tài)學(xué)上農(nóng)業(yè)用地的生態(tài)質(zhì)量高于城市,并與自然地(主要是生態(tài)用地)有著非常異化的生態(tài)質(zhì)量.該研究的假設(shè)模式如圖2所示,即假設(shè)從城市中心(核心建成區(qū))開始,自然地和農(nóng)業(yè)用地之間會(huì)沿著城鄉(xiāng)梯度進(jìn)行分叉,預(yù)測(cè)自然地與建成區(qū)的景觀模式分異(核心與小島)會(huì)促進(jìn)城鄉(xiāng)梯度過渡型的產(chǎn)生.
圖 2 城鄉(xiāng)梯度的假設(shè)模式Fig.2 Hypothetical model of the urban-rural gradient
該研究采用7個(gè)指標(biāo)(3個(gè)土地覆被指標(biāo)與4個(gè)景觀模式指標(biāo))來識(shí)別城鄉(xiāng)梯度(見表1). 建成區(qū)指標(biāo)來源于土地覆被數(shù)據(jù)中的不透水面要素,其表征區(qū)域的城市化水平與經(jīng)濟(jì)社會(huì)狀況,并對(duì)生態(tài)質(zhì)量施加負(fù)面影響[22-23]. 選擇自然地指標(biāo)來表征研究區(qū)的自然環(huán)境特性,該指標(biāo)包含林地、灌木地、濕地3種土地覆被要素. 農(nóng)業(yè)用地是鄉(xiāng)村重要的生產(chǎn)空間,該研究采用耕地要素(狹義的農(nóng)業(yè)用地)來反映鄉(xiāng)村的農(nóng)業(yè)發(fā)展水平. 景觀模式指標(biāo)采用形態(tài)學(xué)空間格局分析方法(MSPA)[24]獲取,發(fā)育良好的核心建成區(qū)表明區(qū)域內(nèi)的建筑較為密集,如商業(yè)區(qū)或工業(yè)園區(qū),而小島建成區(qū)表明孤立的建筑分散于景觀中. 核心自然地說明區(qū)域內(nèi)有較大面積的生物棲息地,意味著具有較為完整的生態(tài)系統(tǒng). 小島自然地則說明區(qū)域內(nèi)有破碎的棲息地或城市公園.
景觀的變化形成了城鄉(xiāng)梯度,也引起生態(tài)系統(tǒng)服務(wù)效應(yīng)[25]. 生態(tài)系統(tǒng)服務(wù)包括支持服務(wù)、供給服務(wù)、調(diào)節(jié)服務(wù)與文化服務(wù)[26],每一類生態(tài)系統(tǒng)服務(wù)的量化指標(biāo)不同,該研究用生境質(zhì)量代表支持服務(wù)、食物供給代表供給服務(wù)、碳儲(chǔ)量代表調(diào)節(jié)服務(wù)、游憩機(jī)會(huì)代表文化服務(wù). 其中游憩機(jī)會(huì)集成了5項(xiàng)評(píng)估指標(biāo),分別為景觀多樣性、自然景觀元素、景點(diǎn)、可達(dá)性和服務(wù)設(shè)施,以綜合表征游憩服務(wù)的空間化特征.
采用主成分分析(PCA)方法降低景觀數(shù)據(jù)集中的維數(shù). PCA是一種多元統(tǒng)計(jì)技術(shù),可將原始指標(biāo)的數(shù)量減少到有限的一組向量中,包含若干關(guān)鍵指標(biāo),這些指標(biāo)能解釋原始指標(biāo)的主要變化,解決數(shù)據(jù)冗余問題[27]. 該研究將主成分分析結(jié)果中初始特征值大于1的指標(biāo)作為主成分,并作為下一步聚類分析的輸入數(shù)據(jù). 主成分可按解釋力排序,第一主成分能解釋的原變量方差最多,第二主成分其次. 通過分析主成分得分及制作散點(diǎn)圖,探討各原始指標(biāo)對(duì)城鄉(xiāng)梯度分異的影響.
AP聚類(親和力傳播聚類)指依據(jù)一定的規(guī)則在節(jié)點(diǎn)之間傳遞信息,在多次迭代過程中出現(xiàn)聚類中心,進(jìn)而實(shí)現(xiàn)數(shù)據(jù)點(diǎn)的自動(dòng)聚類,具有聚類速度快、對(duì)輸入相似度矩陣的三角不等式和對(duì)稱性沒有要求,且可以適用于多種場(chǎng)合等優(yōu)點(diǎn)[28]. AP聚類的優(yōu)勢(shì)是不需人為設(shè)置初始的類中心,主要根據(jù)相似矩陣本身的特征逐漸達(dá)到聚類收斂的效果. AP聚類已被證明比其他聚類方法(如K-均值聚類)更快,對(duì)異常值的誤差和敏感性更低. 該方法已應(yīng)用于景觀度量數(shù)據(jù)集的研究[29],該研究在PCA獲取景觀數(shù)據(jù)集的主成分基礎(chǔ)上,采用AP聚類算法生成村級(jí)尺度單元的城鄉(xiāng)梯度類型.
筆者采用模型法分別對(duì)不同類型的生態(tài)系統(tǒng)服務(wù)進(jìn)行評(píng)估.
采用InVEST 模型來評(píng)估研究區(qū)的生境質(zhì)量與碳儲(chǔ)量. InVEST模型中對(duì)生境質(zhì)量的測(cè)量包括4個(gè)要素:每個(gè)威脅因子的影響相對(duì)權(quán)量、每種生境對(duì)每個(gè)威脅因子的相對(duì)敏感性、生境與威脅因子之間的距離以及生境的可及性. 選擇對(duì)生態(tài)景觀影響較大的不透水面、鐵路、國道、省道、縣鄉(xiāng)道、農(nóng)業(yè)用地定義為威脅源,并參考相關(guān)文獻(xiàn)[30]對(duì)威脅因子的適宜度與敏感性進(jìn)行賦值. InVEST模型對(duì)碳儲(chǔ)量的評(píng)估需要應(yīng)用的基本參數(shù)是不同土地使用類型的碳密度,該研究的碳密度數(shù)據(jù)來源于與該研究區(qū)氣候相近的區(qū)域研究成果[31].
利用土地覆被數(shù)據(jù)和統(tǒng)計(jì)年鑒數(shù)據(jù),計(jì)算研究區(qū)各類型土地覆被的食物總產(chǎn)值,實(shí)現(xiàn)食物供給的空間化. 食物總產(chǎn)值數(shù)據(jù)來源于《南昌市統(tǒng)計(jì)年鑒》(2018年),其中,農(nóng)業(yè)產(chǎn)值對(duì)應(yīng)的是土地覆被數(shù)據(jù)中的耕地,林業(yè)產(chǎn)值對(duì)應(yīng)的是林地,牧業(yè)產(chǎn)值對(duì)應(yīng)草地,漁業(yè)產(chǎn)值對(duì)應(yīng)水體. 食物供給計(jì)算模型參考文獻(xiàn)[32].
參照已有研究的游憩服務(wù)量評(píng)估模型[33],對(duì)游憩機(jī)會(huì)評(píng)價(jià)指標(biāo)進(jìn)行分析. 景觀多樣性指標(biāo)應(yīng)用Fragstats軟件獲取研究區(qū)的香農(nóng)多樣性指數(shù)[34];利用ArcGIS軟件計(jì)算自然景觀元素(林地、濕地、水體)的歐式距離;利用ArcGIS軟件的核密度方法計(jì)算景點(diǎn)(A級(jí)景區(qū)、公園、度假村、休閑農(nóng)莊、漁莊)的空間密度;采用2SFCA法計(jì)算從各街道(鄉(xiāng)鎮(zhèn))點(diǎn)到景點(diǎn)的可達(dá)性[35];采用ArcGIS軟件的泰森多邊形方法計(jì)算服務(wù)設(shè)施的服務(wù)范圍. 然后用自然斷點(diǎn)法將數(shù)據(jù)劃分為10個(gè)等級(jí),最后對(duì)分級(jí)結(jié)果進(jìn)行等權(quán)重空間疊加,得到研究區(qū)總的游憩機(jī)會(huì).
借助SPSS軟件分析4種生態(tài)系統(tǒng)服務(wù)之間的相關(guān)性,通過皮爾遜相關(guān)系數(shù)判斷不同城鄉(xiāng)梯度下的生態(tài)系統(tǒng)服務(wù)權(quán)衡與協(xié)同關(guān)系:當(dāng)兩種服務(wù)的相關(guān)系數(shù)大于0、且呈顯著正相關(guān)(P<0.05)時(shí),則視兩種服務(wù)間的相互作用關(guān)系為協(xié)同;當(dāng)兩種服務(wù)的相關(guān)系數(shù)小于0、且呈顯著負(fù)相關(guān)(P<0.05)時(shí),則視為權(quán)衡;當(dāng)相關(guān)性不顯著時(shí),則兩種生態(tài)系統(tǒng)服務(wù)之間為不相關(guān)關(guān)系.
4.1.1城鄉(xiāng)梯度類型與空間分布
該研究采用AP聚類算法生成研究區(qū)城鄉(xiāng)梯度類型(見圖3),包含5種類型. 城市集中在中部,有80個(gè)村級(jí)單元,是主城區(qū)開發(fā)最成熟的區(qū)域,在空間上最為集聚,經(jīng)濟(jì)發(fā)展水平與建筑密度最高. 近郊村主要圍繞在城市的周邊,有125個(gè)村級(jí)單元. 一部分近郊村的分布較為離散,原因是其受南昌工業(yè)發(fā)展的政策影響,例如主城區(qū)西北部的贛江新區(qū)與東北部的高新區(qū)開發(fā)密度不如城市核心區(qū),且不與城市核心區(qū)緊鄰. 此外,一些地方經(jīng)濟(jì)強(qiáng)鎮(zhèn),如城南的向塘鎮(zhèn)、厚田鎮(zhèn)與城東的麻丘鎮(zhèn),其鎮(zhèn)政府駐地經(jīng)濟(jì)發(fā)展水平較高,但與城市核心區(qū)存在空間上的隔離. 農(nóng)業(yè)村廣布于贛撫平原灌區(qū)上,數(shù)量為461個(gè). 農(nóng)業(yè)村的鮮明特點(diǎn)是擁有較大面積耕地,有大面積的基本農(nóng)田. 遠(yuǎn)郊村有107個(gè),分布較為零散,主要位于林地與水體面積較大的區(qū)域. 生態(tài)村有58個(gè),主要分布于西部的梅嶺山地、東北部的鄱陽湖沿岸以及東南部的白虎嶺林場(chǎng),這些區(qū)域具有自然地面積大的明顯特征.
圖 3 南昌市主城區(qū)城鄉(xiāng)梯度類型Fig.3 Urban-rural gradient types of the main urban area of Nanchang
4.1.2城鄉(xiāng)梯度的主成分解釋
研究區(qū)7個(gè)景觀要素的主成分分析結(jié)果如表2所示,前3個(gè)主成分初始特征值累積比為0.884,因此城鄉(xiāng)梯度聚類的輸入數(shù)據(jù)來自前3個(gè)主成分. 該研究用主成分大于0.6或小于-0.6來表示原始指標(biāo)對(duì)主成分的強(qiáng)關(guān)系. 第1主成分解釋了原始指標(biāo)37.8%的變化,在第1主成分中,自然地與核心自然地呈正相關(guān),且與建成區(qū)、核心建成區(qū)均呈負(fù)相關(guān). 這識(shí)別出了傳統(tǒng)的城鄉(xiāng)梯度二分法的梯度類型,一端是自然地,另一端是建成區(qū). 第2主成分解釋了原始指標(biāo)35.7%的變化,捕獲了另一個(gè)體現(xiàn)城鄉(xiāng)梯度二分法的梯度類型,即建成區(qū)與農(nóng)業(yè)用地. 第3主成分的貢獻(xiàn)率為14.9%,其中小島自然地貢獻(xiàn)最多,為城鄉(xiāng)梯度的多樣式劃分創(chuàng)造了可能.
表 2 景觀要素的主成分分析結(jié)果Table 2 Principal component analysis results
前兩個(gè)主成分共解釋了原始指標(biāo)73.5%的變化,
對(duì)分類結(jié)果產(chǎn)生很強(qiáng)影響. 如圖4所示,5種梯度類型受到前兩個(gè)主成分控制并呈現(xiàn)出明顯的分異現(xiàn)象,表現(xiàn)出明顯的集群效果. 結(jié)合表2可知:主成分1的積極方向是由建成區(qū)、核心建成區(qū)與小島自然地驅(qū)動(dòng),消極方向由自然地和核心自然地驅(qū)動(dòng);主成分2的積極方向由建成區(qū)、核心建成區(qū)和自然地驅(qū)動(dòng),消極方向明顯地由農(nóng)業(yè)用地驅(qū)動(dòng).
圖 4 前2個(gè)主成分與城鄉(xiāng)梯度的關(guān)聯(lián)Fig.4 Correlation between the first two principal components and the urban-rural gradient
圖 5 南昌市主城區(qū)7種景觀的城鄉(xiāng)梯度Fig.5 Urban-rural gradients of seven landscapes in the main urban area of Nanchang
4.1.3城鄉(xiāng)梯度的景觀指標(biāo)解釋
景觀結(jié)構(gòu)在城鄉(xiāng)之間發(fā)生著梯度變化(見圖5).沿著城鄉(xiāng)梯度,建成區(qū)面積與核心建成區(qū)面積占比總體上表現(xiàn)為梯度遞減的態(tài)勢(shì),但遠(yuǎn)郊村的建成區(qū)面積占比略高于農(nóng)業(yè)村. 農(nóng)業(yè)用地面積占比呈現(xiàn)倒“V”型的梯度變化特征,農(nóng)業(yè)用地占比在農(nóng)業(yè)村最高,往城市與生態(tài)村兩端遞減. 從城市往鄉(xiāng)村方向,自然地與核心自然地的面積占比總體上表現(xiàn)為逐漸增長的趨勢(shì),但在近郊村出現(xiàn)異常,近郊村的自然地與核心自然地面積占比略高于城市和農(nóng)業(yè)村.農(nóng)業(yè)村產(chǎn)業(yè)發(fā)展以農(nóng)業(yè)生產(chǎn)為主,大量的土地被開發(fā)為耕地,致使自然地反而少于近郊村. 小島建成區(qū)反映建設(shè)用地的零散分布情況,近郊村的小島建設(shè)用地面積占比略高于其他類型村,并向兩端方向梯度遞減. 原因在于近郊村毗鄰城市,是大學(xué)城、工業(yè)園等發(fā)展要素的布局對(duì)象,但建設(shè)密度不如城市核心區(qū),存在建設(shè)用地孤島現(xiàn)象. 小島自然地的面積占比以近郊村為最高,其次是城市,并往生態(tài)村方向遞減. 小島自然地在城市與近郊村中具有特殊的生態(tài)價(jià)值,然而城市景觀以不透水面為主,連片的自然地稀缺,小島自然地成為關(guān)鍵性的生態(tài)用地.南昌市主城區(qū)近郊村存在湖泊與山體,為小島自然地的保留創(chuàng)造了條件.
圖6以案例村的形式直觀展示了景觀結(jié)構(gòu)的梯度差異,其位于城市核心區(qū)與梅嶺山地、贛撫平原灌區(qū)的銜接處. 北山村南臨南昌都市,是城市類型村,以不透水面景觀占主導(dǎo),不透水面的面積占比近80%.冠山農(nóng)場(chǎng)南鄰北山村,是近郊村,擁有較大面積的不透水面,但農(nóng)業(yè)用地面積與之相當(dāng),同時(shí)存在小島自然地零散分布的現(xiàn)象. 天源村主要位于贛撫平原灌區(qū),是農(nóng)業(yè)村,主要為農(nóng)業(yè)用地,同時(shí)存在小面積的自然地與不透水面. 樂化村是遠(yuǎn)郊村,農(nóng)業(yè)用地最多,也有較大面積的自然地與不透水面. 嶺背林場(chǎng)西接梅嶺,是生態(tài)村,自然地面積廣大,存在小面積的小島建成區(qū).
圖 6 南昌市主城區(qū)的城鄉(xiāng)梯度案例Fig.6 Example of the urban-rural gradient in the main urban area of Nanchang
城市化、農(nóng)業(yè)活動(dòng)、自然環(huán)境共同驅(qū)動(dòng)城鄉(xiāng)梯度分異,導(dǎo)致鄉(xiāng)村的一端出現(xiàn)農(nóng)業(yè)村與生態(tài)村,城鄉(xiāng)兩端間出現(xiàn)過渡型鄉(xiāng)村. 建成區(qū)面積變化是城市化的代理,伴隨著建成區(qū)的擴(kuò)張,農(nóng)業(yè)用地與自然地均受到一定程度的威脅. 該研究將核心自然地作為自然環(huán)境的代理,展示了其與建成區(qū)的城鄉(xiāng)梯度關(guān)系(見圖7),隨著建成區(qū)面積的減少,核心自然地的面積占比并非呈增長趨勢(shì). 在建成區(qū)面積占比小的區(qū)域中分異出3種梯度類型(農(nóng)業(yè)村、遠(yuǎn)郊村、生態(tài)村),農(nóng)業(yè)村與生態(tài)村成為核心自然地面積占比差異較大的兩端. 分類結(jié)果符合前期的假設(shè)模式(見圖2).
4.2.1生態(tài)系統(tǒng)服務(wù)的城鄉(xiāng)梯度特征
南昌市主城區(qū)4種生態(tài)系統(tǒng)服務(wù)的空間分布如圖8所示. 生境質(zhì)量高值區(qū)位于西部的山地森林區(qū)與東北部的鄱陽湖區(qū),碳儲(chǔ)量高值區(qū)位于西部山地,食物供給高值區(qū)位于湖區(qū)與南部的草場(chǎng),游憩機(jī)會(huì)高值區(qū)位于中心城區(qū)與西部山地. 該研究利用ArcGIS軟件對(duì)碳儲(chǔ)量與食物供給數(shù)據(jù)進(jìn)行歸一化處理,再將4類生態(tài)系統(tǒng)服務(wù)的空間數(shù)據(jù)與城鄉(xiāng)梯度識(shí)別結(jié)果進(jìn)行關(guān)聯(lián),制作出生態(tài)系統(tǒng)服務(wù)的城鄉(xiāng)梯度圖(見圖9).
圖 7 建成區(qū)與核心自然地的城鄉(xiāng)梯度關(guān)系Fig.7 Urban-rural gradient relationship between built-up areas and core natural land
圖 8 4種生態(tài)系統(tǒng)服務(wù)的空間分布Fig.8 Spatial distribution of four ecosystem services
沿著城鄉(xiāng)梯度,主城區(qū)生境質(zhì)量逐漸增高,并可分為兩個(gè)梯級(jí). 從城市到農(nóng)業(yè)村,生境質(zhì)量平均值增加0.18,增量較小. 而從農(nóng)業(yè)村到生態(tài)村,生境質(zhì)量平均值增加0.564,反映出生境梯度的顯著變化. 碳儲(chǔ)量在城鄉(xiāng)間同樣表現(xiàn)出兩個(gè)梯級(jí),從城市到農(nóng)業(yè)村,碳儲(chǔ)量的平均值僅增加0.098,而與農(nóng)業(yè)村相比,生態(tài)村的碳儲(chǔ)量平均值增長0.651,碳儲(chǔ)量明顯增加. 生境質(zhì)量與碳儲(chǔ)量是體現(xiàn)區(qū)域生態(tài)質(zhì)量的重要指標(biāo)[36-37],由此可驗(yàn)證圖2的城鄉(xiāng)梯度假設(shè),即在鄉(xiāng)村端生態(tài)質(zhì)量發(fā)生分異,產(chǎn)生了農(nóng)業(yè)村與生態(tài)村這兩種典型的鄉(xiāng)村類型.
沿著城鄉(xiāng)梯度,食物供給呈現(xiàn)為倒“U”型的城鄉(xiāng)梯度特征. 城市與生態(tài)村為低值區(qū),農(nóng)業(yè)村、近郊村與遠(yuǎn)郊村的食物供給水平較高,其中,近郊村的食物供給中位數(shù)高于農(nóng)業(yè)村,遠(yuǎn)郊村食物供給的平均值最高. 圖5雖說明了農(nóng)業(yè)村的農(nóng)業(yè)用地(耕地)比例最高,但食物供給不止來源于耕地,還來源于林地、草地與水體(分別對(duì)應(yīng)林業(yè)、牧業(yè)與漁業(yè)). 研究區(qū)草地與水體的單位面積產(chǎn)出價(jià)值高于耕地,而部分近郊村與遠(yuǎn)郊村的這兩類土地覆被面積多于農(nóng)業(yè)村,致使其在食物供給方面與農(nóng)業(yè)村大體相當(dāng),扮演著提供多樣化農(nóng)產(chǎn)品的角色.
圖 9 4類生態(tài)系統(tǒng)系統(tǒng)服務(wù)的城鄉(xiāng)梯度特征Fig.9 Urban-rural gradient characteristics of four types of ecosystem service
游憩機(jī)會(huì)在城鄉(xiāng)梯度中表現(xiàn)為“V”型的變化特征. 農(nóng)業(yè)村是游憩機(jī)會(huì)的低值區(qū),高值區(qū)的一端是城市與近郊村,原因在于城市與近郊村在景點(diǎn)、可達(dá)性、服務(wù)設(shè)施、景觀多樣性方面優(yōu)于農(nóng)業(yè)村. 高值區(qū)另一端的生態(tài)村游憩機(jī)會(huì)也較高,原因在于生態(tài)村有美學(xué)價(jià)值較高的自然景觀元素. 在主城區(qū)的西部山地,具備臨近中心城區(qū)的優(yōu)勢(shì),可達(dá)性與服務(wù)設(shè)施較好,是居民休閑娛樂的旅游勝地.
4.2.2生態(tài)系統(tǒng)服務(wù)關(guān)系的城鄉(xiāng)梯度
自然地與農(nóng)業(yè)用地是城鄉(xiāng)梯度的識(shí)別指標(biāo),也是引起城鄉(xiāng)梯度分異的重要原因. 自然地與農(nóng)業(yè)用地的代表性生態(tài)系統(tǒng)服務(wù)是生境質(zhì)量與食物供給,分析代表性生態(tài)系統(tǒng)服務(wù)之間及其與其他服務(wù)的權(quán)衡與協(xié)同關(guān)系,將有助于城鄉(xiāng)生態(tài)系統(tǒng)服務(wù)的優(yōu)化與管理.
如表3所示,食物供給與游憩機(jī)會(huì)、碳儲(chǔ)量之間的關(guān)系以權(quán)衡性為主導(dǎo). 從農(nóng)業(yè)村到生態(tài)村,食物供給與碳儲(chǔ)量之間的權(quán)衡性逐漸增強(qiáng)(在生態(tài)村權(quán)衡性最高,相關(guān)系數(shù)為-0.536),說明在鄉(xiāng)村端布局食物生產(chǎn)活動(dòng)會(huì)與碳儲(chǔ)量發(fā)生權(quán)衡作用. 在近郊村、遠(yuǎn)郊村與生態(tài)村,食物供給與游憩機(jī)會(huì)為低權(quán)衡關(guān)系(在遠(yuǎn)郊村權(quán)衡性最高,相關(guān)系數(shù)為-0.412),在這些村布局食物生產(chǎn)活動(dòng)會(huì)與游憩機(jī)會(huì)發(fā)生沖突.
表 3 生態(tài)系統(tǒng)服務(wù)關(guān)系的城鄉(xiāng)梯度差異Table 3 Urban-rural gradient differences of ecosystem services relationships
生境質(zhì)量與游憩機(jī)會(huì)之間的關(guān)系以協(xié)同性為主.在城市、農(nóng)業(yè)村與生態(tài)村,生境與游憩機(jī)會(huì)之間存在著一定的協(xié)同關(guān)系,其中生態(tài)村生境與游憩機(jī)會(huì)的協(xié)同性最高(0.461),說明生態(tài)村的高生境質(zhì)量有助于游憩機(jī)會(huì)的增長. 生境與碳儲(chǔ)量在生態(tài)村有較高協(xié)同性(0.768),然而在近郊村呈低權(quán)衡關(guān)系(-0.443),原因在于近郊村是以農(nóng)業(yè)用地與不透水面占主導(dǎo)的鄉(xiāng)村,自然地面積占比較小,而主城區(qū)近郊的河流、湖泊較多,部分近郊村較高的生境質(zhì)量是以納入水體的方式實(shí)現(xiàn),而其碳儲(chǔ)量并沒有增加.
生境質(zhì)量與食物供給的相互關(guān)系在鄉(xiāng)村端產(chǎn)生明顯分異,在農(nóng)業(yè)村生境質(zhì)量與食物供給之間呈低協(xié)同關(guān)系(0.312),意味著農(nóng)業(yè)村生境質(zhì)量的提高可促進(jìn)食物供給能力的提高. 而在生態(tài)村生境質(zhì)量與食物供給則表現(xiàn)為低權(quán)衡關(guān)系(-0.481),表明在生態(tài)村布局食物生產(chǎn)活動(dòng)會(huì)引發(fā)生境質(zhì)量下降.
城鄉(xiāng)梯度分類指標(biāo)的選取會(huì)對(duì)空間異質(zhì)性研究結(jié)果產(chǎn)生大的影響. 從指標(biāo)的選取來看,與基于單一規(guī)則的劃分方式不同,該研究綜合了3個(gè)土地覆被指標(biāo)與4個(gè)景觀模式指標(biāo). 在土地覆被指標(biāo)方面,典型的城鄉(xiāng)梯度概念假定土地利用強(qiáng)度從發(fā)達(dá)的城市核心向郊區(qū)下降[38]. 鑒于農(nóng)村的土地利用強(qiáng)度低,表現(xiàn)為自然地與農(nóng)業(yè)用地,因此選用了自然地與農(nóng)業(yè)用地指標(biāo),符合南昌市主城區(qū)農(nóng)村區(qū)域的土地利用特征,研究結(jié)果也顯示在農(nóng)村的梯度末端分裂有農(nóng)業(yè)村與生態(tài)村. 該研究借助MSPA方法選取了核心景觀與小島景觀這兩種互不重疊且具有不同生態(tài)學(xué)意義的景觀類型,能夠更好地捕獲城鄉(xiāng)的景觀差別. 然而,有研究認(rèn)為,城市擴(kuò)張除了表現(xiàn)在城市核心景觀與小島景觀類型的差異外,其他景觀類型(例如邊緣區(qū))表現(xiàn)得較為明顯[39]. 此外,橋接區(qū)作為連通核心區(qū)的狹長區(qū)域,對(duì)生物遷移和景觀連接具有特殊意義,可能代表城鄉(xiāng)梯度中的過渡區(qū)域. 由于數(shù)據(jù)分辨率問題,暫未考慮這些指標(biāo),建議在今后的研究中可選取符合區(qū)域?qū)嶋H的土地覆被類型和多樣化的景觀模式指標(biāo).
目前國內(nèi)應(yīng)用最多的城鄉(xiāng)梯度分類方式是基于與城市中心商業(yè)區(qū)的距離,劃分出一定數(shù)量的圈帶,圈帶具有等比例面積的尺度特征. 而采用AP聚類算法考慮到了圈帶內(nèi)的尺度差異性,通過多指標(biāo)融合捕獲了城鄉(xiāng)景觀系統(tǒng)的復(fù)雜性. 例如,南昌市主城區(qū)若采用基于距離的方法會(huì)讓城西的梅嶺山地與近郊的各種開發(fā)區(qū)處于同一梯度,也易忽視南昌市遠(yuǎn)郊存在自然山體、農(nóng)業(yè)區(qū)、濕地保護(hù)區(qū)等復(fù)雜地理背景的事實(shí). AP聚類算法也有一定的局限性,其要求用戶通過調(diào)整相關(guān)參數(shù)指定集聚目標(biāo)的數(shù)量,需要在多個(gè)聚類方案中進(jìn)行檢查與調(diào)整. 筆者結(jié)合南昌市主城區(qū)實(shí)際,確定了該文的分類方案,加入了農(nóng)業(yè)村與生態(tài)村類型,在梯度上有明顯的差異性. 建議將聚類方法應(yīng)用于其他案例研究時(shí),需謹(jǐn)慎考慮分類的數(shù)量.
鑒于空間尺度的選擇會(huì)影響生態(tài)反應(yīng)的變化[40],梯度本身對(duì)空間尺度(例如網(wǎng)格單元大小或緩沖半徑)很敏感,因此城鄉(xiāng)梯度分類方式的選擇直接影響研究者感興趣的生態(tài)過程[41]. 與基于緩沖區(qū)[42]、網(wǎng)格單元模式[43]得出的研究結(jié)果相比,基于行政單元的生態(tài)系統(tǒng)服務(wù)研究的管理導(dǎo)向性更強(qiáng),能為市政決策與鄉(xiāng)村振興提供政策支持. 與基于鄉(xiāng)鎮(zhèn)尺度的研究結(jié)果相比,更敏感地反映城鄉(xiāng)梯度的變化特征,適應(yīng)生態(tài)系統(tǒng)服務(wù)保護(hù)的精細(xì)化要求. 然而,該文僅描述了生態(tài)系統(tǒng)服務(wù)及其相互關(guān)系的城鄉(xiāng)梯度特征,在客觀規(guī)律探索方面,需進(jìn)一步深化對(duì)生態(tài)系統(tǒng)服務(wù)梯度規(guī)律驅(qū)動(dòng)機(jī)制的研究,開展多時(shí)空尺度研究,明確梯度規(guī)律的穩(wěn)健性.
該研究揭示了生態(tài)系統(tǒng)服務(wù)及其相互關(guān)系的城鄉(xiāng)梯度性,可為城鄉(xiāng)生態(tài)系統(tǒng)服務(wù)的提質(zhì)增值與協(xié)同管理提供理論參考. 在提質(zhì)增值方面,應(yīng)優(yōu)先保障農(nóng)業(yè)村的糧食供給能力;鼓勵(lì)近郊村與遠(yuǎn)郊村食物生產(chǎn)的多元化發(fā)展;加大對(duì)生態(tài)村的生境保護(hù)與生態(tài)補(bǔ)償力度;激活遠(yuǎn)郊村、生態(tài)村的游憩服務(wù)潛力并轉(zhuǎn)化為旅游生產(chǎn)力. 在協(xié)同管理方面,需謹(jǐn)慎在生態(tài)村布局食物生產(chǎn)活動(dòng);減少遠(yuǎn)郊村食物供給與固碳、游憩服務(wù)之間的空間矛盾;促進(jìn)農(nóng)業(yè)村生境維護(hù)與食物生產(chǎn)、游憩服務(wù)之間的協(xié)同增益;加強(qiáng)對(duì)城市生境的維護(hù)與營造,滿足城市居民日益增長的游憩需求.
a) 采用基于景觀聚類的方法識(shí)別城鄉(xiāng)梯度,選取3種土地覆被指標(biāo)與4種景觀模式指標(biāo),用主成分分析法得到3個(gè)主成分. 以AP聚類方法生成5種城鄉(xiāng)梯度類型,分別為城市、近郊村、遠(yuǎn)郊村、農(nóng)業(yè)村與生態(tài)村. 通過分析不同城鄉(xiāng)梯度的景觀差異,展示城鄉(xiāng)梯度的案例村,表明了基于景觀聚類的分類方法效果較好,捕捉到了鄉(xiāng)村的末端類型(農(nóng)業(yè)村與生態(tài)村)以及沿梯度中間分布的類型.
b) 城鄉(xiāng)梯度分類結(jié)果顯示,南昌市主城區(qū)的城鄉(xiāng)梯度并非呈圈層結(jié)構(gòu),不同類型村具有不同的空間區(qū)位與集聚特征,與研究區(qū)的城鄉(xiāng)地理環(huán)境背景相符合. 研究發(fā)現(xiàn)城市化、農(nóng)業(yè)活動(dòng)與自然環(huán)境共同驅(qū)動(dòng)著主城區(qū)的城鄉(xiāng)梯度分異.
c) 城鄉(xiāng)梯度是生態(tài)系統(tǒng)服務(wù)空間異質(zhì)性研究的有力工具,研究發(fā)現(xiàn)沿著城鄉(xiāng)梯度,主城區(qū)生境質(zhì)量、固碳服務(wù)可劃分為兩個(gè)增長梯級(jí),食物生產(chǎn)與游憩機(jī)會(huì)分別表現(xiàn)為倒“U”與“V”型的梯度變化特征. 在城鄉(xiāng)梯度中,食物供給與游憩機(jī)會(huì)、固碳服務(wù)之間的關(guān)系以權(quán)衡性為主導(dǎo),生境質(zhì)量與游憩機(jī)會(huì)之間的關(guān)系以協(xié)同性為主,而生境質(zhì)量與食物供給服務(wù)分別在農(nóng)業(yè)村與生態(tài)村呈現(xiàn)為協(xié)同關(guān)系與權(quán)衡關(guān)系. 城鄉(xiāng)梯度下的生態(tài)系統(tǒng)服務(wù)響應(yīng)分析,可為城鄉(xiāng)生態(tài)系統(tǒng)服務(wù)的提質(zhì)增值與協(xié)同管理提供理論參考.