王玲 芮旸 羅飛 孫瑋鴻

摘要 基于種群生態(tài)學理論，將黃河流域971個國家級特色保護類村莊視作一個地理種群。綜合GIS空間分析和數(shù)理統(tǒng)計方法，明晰這一種群的分布動態(tài)和生境選擇特征，進而應(yīng)用地理探測器對比揭示影響該種群分布的生態(tài)因子。結(jié)果表明：① 黃河流域特色保護類村莊的分布密度隨種群增長而大幅提高，分布模式隨尺度變大由顯著聚集轉(zhuǎn)為不顯著離散，分布格局隨時間推移由多中心向單中心演化，但集聚中心一直位于河谷盆地；② 黃河流域特色保護類村莊具有“聚中低山區(qū)、向陽坡斜坡、喜溫暖濕潤、近黃河干流”的偏好，趨向中等人口密度、經(jīng)濟欠發(fā)達地區(qū)分布，主要位于國道沿線30 km、省道沿線10 km、中心城市周邊20～80 km、歷史城鎮(zhèn)周邊60 km、旅游景區(qū)和文保單位周邊20 km范圍內(nèi)；③ 影響特色保護類村莊種群分布的生態(tài)因子具有綜合性、非等價性和限定性，在黃河流域以文保單位和海拔為主導。

關(guān)鍵詞 特色保護類村莊;地理種群;空間分布;生境選擇;黃河流域

中圖分類號：K901.8? DOI：10.16152/j.cnki.xdxbzr.2024-03-014

The population distribution dynamics and habitat selection of characteristic protection villages in the Yellow River Basin

WANG Ling1， RUI Yang1，2， LUO Fei1， SUN Weihong1

1.College of Urban and Environmental Sciences， Northwest University， Xian 710127， China;2.Yellow River Institute of Shaanxi Province， Northwest University， Xian 710127， China）

Abstract Based on the theory of population ecology， 971 national-level characteristic protection villages in the Yellow River Basin were regarded as one geographic population. Based on GIS spatial analysis and mathematical statistics method， the distribution dynamics and habitat selection characteristics of the population were clarified， and then the ecological factors affecting the distribution of this geographic population were revealed by comparing with geographic detectors. The result shows that： ① The distribution density of characteristic protection villages in the Yellow River Basin increased significantly with population growth， the distribution pattern changed from significant aggregation to insignificant dispersion with the increase of scale， and the distribution pattern evolved from multicenter to monocenter over time， but the agglomeration center was always located in the valley basin.② The characteristic protection villages in the Yellow River Basin have the preferences of "gathering in the middle and low mountains， facing the sunny slope slopes， liking warmth and wetness， being close to the main stream of the Yellow River"， and tend to be distributed in areas with medium population density and economically undeveloped areas. It is mainly located within 30 km along the national highway， 10 km along the provincial highway， 20～80 km around the center city， 60 km around the historic town， and 20 km around thetourist attractions and cultural preservation units. ③ The ecological factors influencing the population distribution of characteristic protected villages are comprehensive， non-equivalent and limited， and are dominated by cultural preservation units and altitude in the Yellow River Basin.

Keywords characteristic protection villages; geographic population; spatial distribution; habitat selection; Yellow River Basin

村莊是鄉(xiāng)村地域系統(tǒng)的基本構(gòu)成單元和主要空間形態(tài)，因功能狀態(tài)的差異而分化為特色村、中心村、空心村等類型［1］。針對不同類型村莊分類開展研究和實踐是推進鄉(xiāng)村振興的必然要求，也是深化鄉(xiāng)村地域系統(tǒng)理論認識的重要路徑。特色保護類村莊是國家《鄉(xiāng)村振興戰(zhàn)略規(guī)劃（2018—2022年）》確定的4種村莊類型之一，是特色資源豐富、保護價值重要的一類村莊的統(tǒng)稱。其源于歷史文化名村，以傳統(tǒng)村落為主體，還包含特色景觀旅游名村和少數(shù)民族特色村寨，是具有政策意涵的新概念，吸引了學界的關(guān)注。既有特色保護類村莊研究大多聚焦村莊個體和微觀尺度，以單案例分析為主，涉及案例村保護與發(fā)展的博弈［2］、鄉(xiāng)村旅游的適宜性［3］、空間重構(gòu)和發(fā)展策略［4-5］、規(guī)劃編制和景觀設(shè)計［6-7］等主題;也有個別研究注意到特色保護類村莊的種群屬性［8］，開始在全國這一宏觀尺度探究該類村莊個體所不具備的群體空間分布特征及其動態(tài)，但對群體生境選擇特征及影響種群特征的生態(tài)因子揭示不足，對地理單元和中觀尺度缺乏關(guān)注。

流域作為相對獨立的地理單元，兼具自然生態(tài)系統(tǒng)的完整性和文化景觀的相對一致性，是研究人類聚居現(xiàn)象的更優(yōu)尺度。黃河流域是華夏先民最早活動的地理單元［9］，人類聚居歷史悠久，聚落遺產(chǎn)資源富集，聚落地域特色鮮明。改革開放以來城鎮(zhèn)化的快速推進和城鄉(xiāng)建設(shè)的大規(guī)模開展，使得該流域的人居環(huán)境發(fā)生了劇變，流域內(nèi)的特色鄉(xiāng)村面臨存續(xù)危機［10］，亟待加強保護利用和科學研究。相關(guān)研究已對黃河流域傳統(tǒng)村落的空間分布特征［11-13］、旅游響應(yīng)度［14］、旅游適應(yīng)性［15］，以及與非物質(zhì)文化遺產(chǎn)的空間關(guān)系［16］等進行了分析，取得了一定進展。但在國家分類推進鄉(xiāng)村振興背景下，仍有必要將研究對象由傳統(tǒng)村落拓展為包含其在內(nèi)的特色保護類村莊，以期更好地響應(yīng)國家戰(zhàn)略。

綜上，本文基于種群生態(tài)學理論，將黃河流域971個國家級特色保護類村莊視作一個地理種群開展研究。首先，采用最近鄰指數(shù)、Ripleys K函數(shù)和核密度分析等方法，刻畫該種群的分布動態(tài)特征;進而，借助與全國和長江流域的比較，明晰黃河流域特色保護類村莊種群的生境選擇特征;最后，應(yīng)用地理探測器，對比揭示影響該種群分布的生態(tài)因子及其獨特性，從而為該流域特色保護類村莊的振興決策提供科學依據(jù)，為鄉(xiāng)村聚落研究提供新的理論視角。

1 數(shù)據(jù)與方法

1.1 研究區(qū)域與對象

本文以黃河流域為研究區(qū)。黃河流域位于95°53′E～119°05′E、32°10′N～41°50′N之間， 自西向東橫跨中國地勢三級階梯和青藏高原、 內(nèi)蒙古高原、 黃土高原、 華北平原等四大地貌單元，涉及青海、 四川、 甘肅、 寧夏、 內(nèi)蒙古、 山西、 陜西、 河南和山東9個省區(qū)， 總面積79.5×104 km2（見圖1）。

本文研究對象是黃河流域特色保護類村莊種群，即在流域范圍內(nèi)的該類村莊個體的總和?？紤]數(shù)據(jù)的權(quán)威性、可獲性和一致性，以截至2023年底至少獲評中國歷史文化名村、中國傳統(tǒng)村落、中國少數(shù)民族特色村寨和全國特色景觀旅游名村四者中其一的971個國家級特色保護類村莊為研究樣本（見圖1）。同時，因在分析研究對象的生境選擇特征和分布影響因子時涉及與全國和長江流域同類種群的對比，本文研究樣本實際包括了全部9 474個國家級特色保護類村莊（其中長江流域3 411個）。

1.2 數(shù)據(jù)來源及處理

研究樣本的原始數(shù)據(jù)中，中國歷史文化名村和全國特色景觀旅游名村名單均來源于住房和城鄉(xiāng)建設(shè)部網(wǎng)站；中國傳統(tǒng)村落名單來源于中國傳統(tǒng)村落數(shù)字博物館（https：∥www.dmctv.cn）；中國少數(shù)民族特色村寨名單來源于國家民族事務(wù)委員會網(wǎng)站。自然地理數(shù)據(jù)中，黃河流域邊界數(shù)據(jù)來源于中國工程科技知識中心（https：∥www.ckcest.cn），海拔、地形起伏度、坡度、坡向等地形數(shù)據(jù)來源于該中心提供的SRTM 30 m DEM數(shù)據(jù);氣溫、降水量、河流等數(shù)據(jù)均來源于中國科學院資源環(huán)境科學數(shù)據(jù)中心（https：∥www.resdc.cn）。社會經(jīng)濟數(shù)據(jù)中，省市行政區(qū)劃邊界和道路、人口、GDP等數(shù)據(jù)集也來源于中國科學院資源環(huán)境科學數(shù)據(jù)中心;全國重點文物保護單位和國家歷史文化名城名單來源于中國政府網(wǎng);中國歷史文化名鎮(zhèn)名單來源于住房和城鄉(xiāng)建設(shè)部網(wǎng)站;國家5A級旅游景區(qū)名單來源于文化和旅游部網(wǎng)站;4A級旅游景區(qū)名單來源于各省份旅游主管部門網(wǎng)站。研究樣本等點狀地理要素的經(jīng)緯度坐標通過綜合百度地圖和高德地圖API軟件獲取，并在ArcGIS中進行地理配準和可視化分析。

1.3 研究方法

1.3.1 最近鄰指數(shù)

最近鄰指數(shù)是表征點狀要素之間鄰近性的地理指標，可用以判斷黃河流域特色保護類村莊種群的空間分布類型。公式為［13］

R=1E（1）

E=12D=12n/A（2）

式中：R為最近鄰指數(shù);1和E為實際最鄰近距離和理論最鄰近距離;n和D為特色保護類村莊種群規(guī)模和分布密度;A為黃河流域面積。 R＜1時，種內(nèi)個體集群分布;R=1時，個體隨機分布;R>1時，個體均勻分布。

1.3.2 Ripleys K函數(shù)

地理格局具有尺度依存特性，不同尺度下挖掘出的空間模式不同［17］。Ripleys K函數(shù)可揭示點狀地理要素在不同空間尺度下的分布模式及其轉(zhuǎn)換規(guī)律。公式為［18］

L（d）=A∑ni=1∑nj=1，j≠ik（i，j）πn（n-1）（3）

式中：n、A含義同式（2）;d為距離尺度;k（i，j）為權(quán)重。L（d）的觀測值高于期望值表示要素呈聚集分布，若觀測值同時高于置信區(qū)間上限，則在對應(yīng)距離d內(nèi)的空間聚集具有統(tǒng)計顯著性;L（d）的觀測值低于期望值表示要素呈離散分布，若觀測值同時低于置信區(qū)間下限，則在對應(yīng)距離d內(nèi)的離散具有顯著性。觀測值與期望值的差值達到最大時所對應(yīng)的L（d）值是要素聚集強度的峰值，對應(yīng)的距離d則表征要素聚集規(guī)模。

1.3.3 核密度分析

核密度分析是一種估計點狀要素在不同空間位置上概率密度函數(shù)的非參數(shù)方法，在地理要素分布研究中應(yīng)用廣泛，可用以直觀刻畫黃河流域特色保護類村莊種群的空間集聚格局。公式為［8］

fn（x）=1nh∑ni=1k（x-xih）（4）

式中：fn（x）為核密度值;k（x-xih）表示核函數(shù);（x-xi）為估值點x到被觀測點xi的距離;h為帶寬，經(jīng)多次對比調(diào)試后確定為80 km;n為帶寬內(nèi)的點數(shù)。

1.3.4 緩沖區(qū)分析

緩沖區(qū)分析是解決地理鄰近性問題的空間分析工具。通過對地理要素按設(shè)定的距離條件，在其周圍建立一定數(shù)量和寬度的多邊形區(qū)域，分析區(qū)內(nèi)的空間數(shù)據(jù)以揭示要素影響范圍和機理［1］。本文采用ArcGIS中的緩沖區(qū)分析工具，對特色保護類村莊種群分布與河流、主要道路、中心城市、歷史城鎮(zhèn)、文保單位和旅游景區(qū)的關(guān)聯(lián)規(guī)律進行描述性統(tǒng)計分析。

1.3.5 空間關(guān)聯(lián)指數(shù)

空間關(guān)聯(lián)指數(shù)源于群落生態(tài)學中關(guān)于景觀要素空間關(guān)系的研究［16］，可定量反映特色保護類村莊種群與中心城市、歷史城鎮(zhèn)、文保單位和旅游景區(qū)等相關(guān)種群的數(shù)量空間關(guān)系。公式為［19］

R=ad-bc（a+b）（c+d）（a+c）（b+d）（5）

式中：R為空間關(guān)聯(lián)指數(shù);a為全部樣方中同時包含兩類種群個體的樣方數(shù);b為全部樣方中僅包含特色保護類村莊個體的樣方數(shù);c為全部樣方中僅包含其他種群個體的樣方數(shù);d為全部樣方中同時不包含兩類種群個體的樣方數(shù)。R值范圍為［-1，1］，R>0為正關(guān)聯(lián)，R<0為負關(guān)聯(lián)。并用式（6）對R值進行顯著性檢驗。

X2=n（ad-bc）2（a+b）（c+d）（a+c）（b+d）（6）

式中：a、b、c、d含義同式（5）;n為樣方總數(shù)。查閱卡方分布表，在顯著水平為0.05自由度為1的條件下，上分位數(shù)為3.84，若｜X2｜>3.84，表明兩類種群的空間關(guān)聯(lián)顯著，反之則不顯著。

1.3.6 地理探測器

地理探測器是量度和挖掘地理事物空間分異性及其背后驅(qū)動力的空間統(tǒng)計學方法［20］。本文借助其中的因子探測器模塊，明晰顯著影響特色保護類村莊分布異質(zhì)性的生態(tài)因子。生態(tài)因子影響力的大小以q值度量，公式為［21］

q=1-1Nσ2∑ni=1Niσ2i（7）

式中：n為變量的分層;N和Ni為研究區(qū)和第i層的樣本數(shù);σ2和σ2i為因變量在整個研究區(qū)和第i層內(nèi)的方差。q值范圍為［0，1］，q值越大表明自變量對因變量的解釋力越強。

2 結(jié)果與分析

2.1 特色保護類村莊種群分布動態(tài)特征

2.1.1 分布密度變化

種群密度是種群最基本的特征和種群數(shù)量統(tǒng)計最常用的指標。如圖2所示，黃河流域特色保護類村莊種群由小變大、由疏到密，數(shù)量經(jīng)歷了前期緩慢增長、后期階梯式快速增長的變動過程，由2003年的2個增至2023年的971個;密度也呈同樣變化趨勢，相應(yīng)由0.03個／104 km2增至12.21個／104 km2。該種群數(shù)量的劇烈變化主要是政策因素調(diào)節(jié)的結(jié)果。2013年中央一號文件首次提出要“制定專門規(guī)劃，啟動專項工程，加大力度保護有歷史文化價值和民族、地域元素的傳統(tǒng)村落和民居”，到該年底黃河流域特色保護類村莊種群數(shù)量突破100達到114，密度超過1個／104 km2，達1.43個／104 km2。此后，6次中國傳統(tǒng)村落大規(guī)模調(diào)查和3批中國少數(shù)民族特色村寨命名掛牌工作的開展，使得特色保護類村莊種群的規(guī)模迅速擴大，種群承載的中華優(yōu)秀傳統(tǒng)文化基因多樣性得以保持。到2023年底，黃河流域該種群數(shù)量接近1 000個，密度超過10個／104 km2。除2013年和2023年外，2018年也是特色保護類村莊種群動態(tài)變化的一個重要的時間節(jié)點。該年中共中央、國務(wù)院印發(fā)的《鄉(xiāng)村振興戰(zhàn)略規(guī)劃》明確提出了特色保護類村莊的概念，標志著該類村莊種群的保護與發(fā)展已上升到國家戰(zhàn)略層面。由此可見，政策始終是對特色保護類村莊種群影響最強烈的外源性調(diào)節(jié)因素。

2.1.2 分布模式變化

如上所述， 2013年、 2018年和2023年是黃河流域特色保護類村莊種群增長過程中的3個重要時點， 故選取這3個時點對種群空間分布模式的變化進行分析。 結(jié)果顯示， 黃河流域特色保護類村莊種群的最近鄰指數(shù)持續(xù)減小， 由2013年的0.525減至2018年的0.461， 再減至2023年的0.449。 可見， 這一種群內(nèi)的特色保護類村莊分布相當不均， 一直成群密集分布， 且聚集程度不斷增強。

進一步運用ArcGIS中基于Ripleys K函數(shù)的多距離空間聚類分析工具， 揭示黃河流域特色保護類村莊種群空間分布模式的尺度變化規(guī)律（見圖3）。 結(jié)果表明， 這一種群的空間分布在3個年份均存在多尺度特征， 分布模式隨距離增大由顯著聚集變?yōu)椴伙@著離散， 但模式轉(zhuǎn)換的特征空間尺度略有差異， 分別為326 km、 312 km、 319 km。 進一步分析這一種群空間聚集模式的變化可知， 該種群在2013年、 2018年和2023年L（d）達到峰值的空間距離分別為52 km、 72 km和69 km， 相應(yīng)的L（d）峰值分別為117、 135和141。 可見， 黃河流域特色保護類村莊種群在312 km內(nèi)始終在空間上呈顯著聚集分布，且聚集強度隨著時間推移持續(xù)增大， 聚集規(guī)模也整體表現(xiàn)出擴大趨勢。

2.1.3 分布格局變化

分別對2013年、2018年和2023年的黃河流域特色保護類村莊種群進行核密度制圖（見圖4），進一步明晰種群空間集聚區(qū)域和格局。結(jié)果顯示，這一種群的分布格局總體呈現(xiàn)出由多中心向單中心演化的趨勢，但集聚中心一直位于河谷盆地，緊鄰中國地勢階梯分界線和省界等地理界線，與黃河干流及一級支流的關(guān)系密切，折射出地形、水資源和區(qū)位條件的影響。2013年，該種群內(nèi)的特色保護類村莊高密度分布在青甘、晉豫、晉陜交界地區(qū)和山西中部地區(qū)，在黃河上游的湟水谷地和黃河中游的太原盆地、晉城盆地最為集聚;到2018年，隨著集聚程度的增強，此前的主要集聚區(qū)出現(xiàn)連片趨勢，由5個演變?yōu)?個;到2023年，黃河流域特色保護類村莊種群最終形成了以山西省晉城市為單核的空間集聚格局，核密度的峰值由2013年的21.917個／104 km2大幅增至205.770個／104 km2。

2.2 特色保護類村莊種群生境選擇特征

生境選擇是生物與環(huán)境之間的基本關(guān)系，反映某一生物個體或種群對特定生境的偏好。生境概念因而具有相對意義，總是針對特定主體而言［22］。在生態(tài)學理論中，生境以生物為主體，指生物個體或種群賴以生存的環(huán)境，由對生物有直接或間接影響的全部生態(tài)因子構(gòu)成。本文中的生境以特色保護類村莊為主體，是對該類村莊形成和演化起作用的各類生態(tài)因子的總和。生態(tài)因子按其性質(zhì)一般分為非生物因子、生物因子和人為因子3類。本文參照相關(guān)文獻［1，8，11-13，23］，以地形、氣候、河流等要素表征非生物因子，以城鎮(zhèn)、文物、景區(qū)等要素表征生物因子， 以人口、 經(jīng)濟、 交通等要素表征人為因子， 綜合分析影響特色保護類村莊個體發(fā)育和種群發(fā)展的生態(tài)因子（見圖5）。 同時， 鑒于不同地理種群經(jīng)過長期變異、 選擇和適應(yīng)后， 對生態(tài)因子的要求往往有異，因而本研究通過與長江流域和全國同類種群的比較， 明晰黃河流域特色保護類村莊種群的生境選擇特征及其特殊性。

2.2.1 對非生物因子的選擇

1） 地形。地形條件是在村莊選址和布局中發(fā)揮基礎(chǔ)性作用的生態(tài)因子。① 從海拔看〔見圖6（a）〕，與全國及長江流域不同，黃河流域特色保護類村莊種群最主要分布在海拔500～2 000 m的中低山區(qū)，而非海拔1 000 m以下的低山丘陵和平原地區(qū)，位于海拔500 m以下地區(qū)的村莊數(shù)量占比（9.47％），遠低于全國（45.28％）和長江流域（47.32％）。② 從地形起伏度看〔見圖6（b）〕，與全國及長江流域相同，黃河流域特色保護類村莊種群亦集中分布在起伏度100 m以下，即地面平坦或起伏小的地區(qū)，且其內(nèi)村莊占比達94.85％，還高于全國（84.74％）和長江流域（81.41％）10個百分點以上。③ 從坡度看〔見圖6（c）〕，黃河流域80.64％的特色保護類村莊分布在0.5°～15°的各類斜坡地帶，長江流域相近比例（79.71％）和全國3／4以上（75.70％）的該類村莊則分布在坡度更大（2°～35°）的緩斜坡、斜坡和陡坡地。④ 從坡向看〔見圖6（d）〕，與全國及長江流域一致，黃河流域的特色保護類村莊也多分布在光熱條件更好的陽坡和半陽坡，且后者中位于陽坡和陰坡的村莊比例較之前兩者相應(yīng)更大和更小。

2） 氣候。氣候條件及其強約束下的農(nóng)業(yè)生產(chǎn)活動對村莊形成和發(fā)展有直接影響，不同水熱條件下村莊數(shù)量分布存在顯著差異。① 從氣溫看〔見圖6（e）〕，受所處緯度位置、地形等因素的限制，黃河流域的特色保護類村莊數(shù)量在年平均氣溫10 ℃～15 ℃的區(qū)域最多（占比56.85％），而非與長江流域和全國一樣，最主要分布在年平均氣溫15 ℃～20 ℃即熱量資源更充足的區(qū)域。② 從降水看〔見圖6（f）〕，受海陸位置等導致的中國降水空間分布規(guī)律影響，黃河流域80.95％的特色保護類村莊都位于年降水量在400～800 mm之間的半濕潤區(qū)，全國73.26％和長江流域94.52％的特色保護類村莊則集中在水分條件更好、年降水量800 mm以上的濕潤區(qū)。

3） 河流。河流等水系作為水源命脈和航運通道，是鄉(xiāng)村聚落形成的基礎(chǔ)和城鄉(xiāng)經(jīng)濟發(fā)展的保障。分別以3 km和10 km為半徑對一級至五級河流建立緩沖區(qū)進行相交分析，結(jié)果顯示（見表1）：黃河流域位于不同距離緩沖區(qū)內(nèi)的特色保護類村莊數(shù)量占比與全國接近，均有三成以上（30.79％和30.03％）的村莊分布在沿河3 km范圍內(nèi)，有約2／3（64.78％和67.43％）的村莊分布在沿河10 km范圍內(nèi)?？梢姡恿魇翘厣Ｗo類村莊種群形成和發(fā)展的重要依托，這種近水分布特征在長江流域表現(xiàn)得更為明顯。進一步分不同等級河流看（見表1），區(qū)別于長江流域和全國，黃河流域特色保護類村莊種群對一級河流的依賴程度明顯大于長江流域和全國，約10％和20％的村莊分布在黃河干流沿線3 km和10 km范圍內(nèi)?？梢姡S河作為中華民族的母親河，是華夏先民臨水而居的重要選擇。

2.2.2 對生物因子的選擇

1） 城鎮(zhèn)。鄉(xiāng)村與城鎮(zhèn)是相互依存的命運共同體。為明晰特色保護類村莊與歷史上、當下作為區(qū)域中心或?qū)λ趨^(qū)域發(fā)展有重要影響的城鎮(zhèn)的關(guān)聯(lián)特征，選取4個直轄市及333個地級行政區(qū)政府駐地城市作為當前時期的中心城市，對其建立以20 km為等間距的10級緩沖區(qū)，統(tǒng)計得到各級緩沖區(qū)內(nèi)的特色保護類村莊數(shù)量和比例〔見圖7（a）〕。與全國及長江流域相同，黃河流域特色保護類村莊亦主要位于距離城市20～80 km的區(qū)域，但在這一區(qū)域的數(shù)量占比（75.80％）要高于前兩者（67.72％和63.73％）不少。換言之，特色保護類村莊與中心城市的關(guān)系在黃河流域更為緊密，兩者在黃河流域的平均距離和空間關(guān)聯(lián)指數(shù)也分別小于和大于在長江流域和全國的同一指標值（見表2）。同時，將454個國家級歷史文化名城名鎮(zhèn)作為歷史時期的重要城鎮(zhèn)，同理得到其10級緩沖區(qū)內(nèi)的特色保護類村莊數(shù)量和比例，發(fā)現(xiàn)總體呈現(xiàn)距離衰減趨勢〔見圖7（b）〕。這一趨勢在黃河流域表現(xiàn)得最為顯著，僅在歷史城鎮(zhèn)周邊20 km范圍內(nèi)就分布有40.47％的特色保護類村莊。綜合平均距離和空間關(guān)聯(lián)指數(shù)看（見表2），特色保護類村莊與歷史城鎮(zhèn)的關(guān)系較之中心城市更密切，反映出這類村莊久遠的建村歷史和近城的選址特征。

2） 文物。特色保護類村莊植根于地域文化，以文物古跡為重要依托和價值體現(xiàn)，本身也是活態(tài)文化遺產(chǎn)。對5 058處全國重點文物保護單位（下文簡稱文保單位）建立以10 km為等間距的10級緩沖區(qū)〔見圖7（c）〕。統(tǒng)計分析各級緩沖區(qū)內(nèi)的特色保護類村莊數(shù)量和比例，發(fā)現(xiàn)距離衰減特征在黃河流域表現(xiàn)得最明顯，村莊數(shù)量隨著與文保單位距離的增加而快速減少。在文保單位10 km緩沖區(qū)內(nèi)已集中了黃河流域56.33％的特色保護類村莊，這一比例遠高于全國及長江流域。平均距離和空間關(guān)聯(lián)指數(shù)的計算結(jié)果也證實了這一點（見表2），即特色保護類村莊種群與文保單位種群的空間關(guān)系在黃河流域更為緊密。

3） 景區(qū)。特色保護類村莊以利用特色資源發(fā)展鄉(xiāng)村旅游為主要振興路徑，與旅游景區(qū)具有共生基礎(chǔ)和協(xié)同潛力。對4 931個國家4A級及以上旅游景區(qū)建立10 km等間距的10級緩沖區(qū)〔見圖7（d）〕。統(tǒng)計分析后發(fā)現(xiàn)特色保護類村莊主要分布在距離這些重點旅游景區(qū)不到20 km的區(qū)域，且這一區(qū)域內(nèi)的村莊數(shù)量占比在黃河流域要分別高于全國和長江流域7個和10個百分點以上。黃河流域特色保護類村莊與旅游景區(qū)存在及潛在的協(xié)同關(guān)系從兩者相對較短的平均距離也可以看出（見表2），因而有必要在該流域促進村景聯(lián)動發(fā)展。

2.2.3 對人為因子的選擇

1）人口。人口是社會生產(chǎn)的主體，具有一定的人口是鄉(xiāng)村聚落存在的前提。根據(jù)已有關(guān)于中國人口密度的分級和分區(qū)標準［24］研究，統(tǒng)計不同密度類型區(qū)內(nèi)的村莊數(shù)量和比例〔見圖8（a）〕。與全國及長江流域相同，黃河流域的特色保護類村莊數(shù)量在一般過渡區(qū)（人口密度為100～200人／km2）和低度核心區(qū)（人口密度為200～400人／km2）最多，但在人口密度400人／km2以上的區(qū)域內(nèi)數(shù)量占比分別高于全國和長江流域6個和10個百分點以上。人口過于集聚和稀疏的區(qū)域都不利于村莊存留，因此，黃河流域特色保護類村莊主要分布在中等人口密度區(qū)域，但相對更趨向于人口密集區(qū)域。

2） 經(jīng)濟。村莊因人口集聚而形成，因經(jīng)濟集聚而發(fā)展，但特色保護類村莊分布與經(jīng)濟因子的關(guān)系不同于人口。由圖8（b）可知，黃河流域的特色保護類村莊數(shù)量與經(jīng)濟發(fā)展水平成反比，一定區(qū)域的地均GDP越高，其內(nèi)的村莊數(shù)量和比例越低。這種分布的經(jīng)濟欠發(fā)達地區(qū)指向在全國和長江流域表現(xiàn)得更為明顯，黃河流域有七成以上（71.88％）的特色保護類村莊分布在地均GDP低于1 500萬元／km2的區(qū)域，在同一區(qū)域全國和長江流域則各有76.84％和83.20％的村莊分布。究其原因，這類經(jīng)濟發(fā)展滯后的區(qū)域雖不利于特色保護類村莊的振興，但卻因較為緩慢的工業(yè)化和城市化進程而為村莊及其特色的保留創(chuàng)造了有利環(huán)境。

3） 交通。與主要交通通道的空間距離和區(qū)位關(guān)系影響村莊與外部環(huán)境的物質(zhì)、能量和信息交換，是村莊選址布局的重要考慮因素。對國道、省道分別建立以10 km、5 km為等間距的10級緩沖區(qū)，統(tǒng)計分析結(jié)果如圖8（c）、（d）所示。特色保護類村莊的近路分布和交通依賴特征顯著，其與國道和省道的距離越近，分布越多，數(shù)量和比例隨著距離增加呈現(xiàn)衰減趨勢，且這種趨勢在黃河流域表現(xiàn)得更為明顯。該流域國道10 km和省道5 km緩沖區(qū)內(nèi)的特色保護類村莊數(shù)量占比已分別超過1／3（34.50％）和1/5（50.67％），均高于全國及長江流域。總體而言，黃河流域特色保護類村莊種群的分布與主要道路密切相關(guān)，尤其集中在國道30 km和省道10 km緩沖區(qū)范圍內(nèi)，其內(nèi)的村莊數(shù)量已分別占到總數(shù)的71.99％和73.84％。

2.3 特色保護類村莊種群分布影響因子

基于上文揭示的9類生態(tài)因子，選取11個自變量，同時以特色保護類村莊數(shù)量為因變量進行地理探測。探測以50 km×50 km的網(wǎng)格為基本單元，黃河流域、長江流域和全國分別在ArcGIS中劃分出416個、858個和4 277個網(wǎng)格。根據(jù)地理探測器運行要求，綜合專業(yè)知識、分類方法和前文分析結(jié)果，對自變量進行離散化處理，將其轉(zhuǎn)化為分類變量。其中，海拔的分類斷點為200 m、500 m、1 000 m、2 000 m、3 000 m和4 000 m；年降水量的分類斷點為200 mm、400 mm、800 mm、1 000 mm和1 600 mm;坡度和年平均氣溫的分級與前文保持一致;其余變量均采用自然斷點法分為8類。影響黃河流域、長江流域和全國特色保護類村莊種群分布的生態(tài)因子探測結(jié)果如表3所示。進而基于探測所得q值，采用SPSS中的K均值聚類算法，將通過顯著性檢驗的因子按影響力大小分為3類，明晰主導生態(tài)因子。

就黃河流域而言， 研究選取的11個自變量均通過了顯著性檢驗， 表明該流域特色保護類村莊種群的分布格局是生物、 非生物和人為3類生態(tài)因子共同作用的結(jié)果。 分因子類型看，生物因子的綜合解釋力最高（平均q值為0.186）， 人文因子的解釋力次之（平均q值為0.156）， 非生物因子的解釋力最低（平均q值為0.100）。 可見， 黃河流域特色保護類村莊種群的分布更多是受種間關(guān)系和人類活動的影響。 但從具體因子看，除全國重點文物保護單位外， 海拔這一非生物因子也在黃河流域特色保護類村莊分布中起著主要作用。

比較而言，顯著影響長江流域特色保護類村莊種群分布的因子亦同樣涉及全部3類生態(tài)因子，其中起主要作用的因子也包括海拔。除海拔外的另一主導生態(tài)因子是年降水量，也屬于非生物因子。并且在該流域，非生物因子對特色保護類村莊種群的整體影響力也最強?？梢?，區(qū)別于黃河流域，長江流域特色保護類村莊種群的空間分異主要由非生物因子尤其是降水因子和地形因子決定。將研究范圍由流域擴大到全國，所得結(jié)果又有差異。主導全國特色保護類村莊分布的顯著性因子則是氣候因子和交通因子。綜上可知，影響特色保護類村莊種群的生態(tài)因子具有綜合性、非等價性和限定性。該類村莊的不同地理種群都受到各類、多個生態(tài)因子的綜合作用，但對特定地理種群分布起作用的因子及主導因子有差異。不同地理種群的發(fā)展需要不同的生態(tài)因子或生態(tài)因子的不同強度。

3 結(jié)論與討論

3.1 結(jié)論

1）黃河流域特色保護類村莊的分布密度隨種群增長而大幅提高，分布模式隨尺度變大由顯著聚集轉(zhuǎn)為不顯著離散，分布格局隨時間推移由多中心向單中心演化，但在312 km內(nèi)始終成群密集分布，且集聚強度持續(xù)增強，集聚規(guī)模整體擴大，集聚中心一直位于河谷盆地并緊鄰地勢階梯分界線和省界。

2）黃河流域特色保護類村莊種群表現(xiàn)出“聚中低山區(qū)、向陽坡斜坡、喜溫暖濕潤、近黃河干流”的自然生境偏好，主要分布在人口密度為50～400人／km2、經(jīng)濟密度低于1 500萬元／km2的區(qū)域，以及國道沿線30 km、省道沿線10 km、中心城市周邊20～80 km、歷史城鎮(zhèn)周邊60 km、旅游景區(qū)和文保單位周邊20 km范圍內(nèi)。

3）影響特色保護類村莊種群的生態(tài)因子具有綜合性、非等價性和限定性。不同于在全國以氣候和交通因子為主導，在長江流域以降水和地形因子為主導，在黃河流域起主要作用的顯著因子是全國重點文物保護單位的數(shù)量和海拔。

3.2 討論

1）結(jié)果討論。不同于本文中黃河流域特色保護類村莊種群分布由文化和地形因素主導的結(jié)論，相關(guān)研究或認為河網(wǎng)密度、年降水量等自然條件是最主要的影響因素［13］，或認為人口數(shù)量、GDP等社會經(jīng)濟因素的影響最大［25］。上述研究都采用地理探測器進行歸因，但結(jié)果存在明顯差異，可能存在3個原因：一是樣本類型和規(guī)模不同，本文采用的是包含傳統(tǒng)村落、歷史文化名村在內(nèi)的特色保護類村莊的最新數(shù)據(jù);二是變量離散化處理方式不同，本文綜合運用了專業(yè)知識和分類算法，并在利用自然斷點法分類時按8級而非5級劃分;三是定位精度和分析面元的不同，綜合運用百度地圖和高德地圖API解析樣本地址，并采用更小尺度的50 km×50 km網(wǎng)格而非100 km×100 km網(wǎng)格進行空間統(tǒng)計，從而可能引發(fā)尺度效應(yīng)和分區(qū)效應(yīng)。

2）政策建議。① 黃河流域特色保護類村莊種群發(fā)展應(yīng)以黃河國家文化公園建設(shè)為依托，以種群密度最高的晉城盆地為中心，以種群分布集中的黃河干流及其支流湟水和汾河為廊道，構(gòu)建“一心引領(lǐng)、三廊支撐”的總體空間格局，推動流域特色保護類村莊的保護利用由點狀向帶狀和網(wǎng)絡(luò)狀轉(zhuǎn)變，全面彰顯和傳承黃河文化。② 黃河流域特色保護類村莊種內(nèi)個體的發(fā)展應(yīng)因勢利導、分類施策，對于鄰近區(qū)域中心城市、高等級旅游景區(qū)等優(yōu)勢種的村莊，應(yīng)利用區(qū)位優(yōu)勢發(fā)展鄉(xiāng)村休閑旅游和特色產(chǎn)業(yè)，不斷增強城村和村景間的正相互作用和互補性，提高個體適應(yīng)度;對于地處省界、高山陡坡等不利生境的村莊，應(yīng)以防止遺傳變異性喪失為目標，加強文化生態(tài)的整體保護和文化基因的活態(tài)傳承。

3）研究展望。地理探測器作為新的空間歸因方法在使用中仍面臨若干問題，典型如自變量的q值多大才能被判定為主導因子？本文嘗試通過K-均值聚類算法而非主觀判斷來明晰主導因子，但在確定分級數(shù)量時仍存在一定主觀性。后續(xù)還需通過對自然斷點法、均值標準差法等數(shù)據(jù)離散化方法效果的比較，探尋一種更加客觀和權(quán)威的主導因子識別方法。此外，從研究區(qū)域看，未來應(yīng)進一步聚焦湟水谷地、太原盆地、晉城盆地等特色保護類村莊密集分布的地理單元，在更小尺度開展地理種群研究，以期深化對黃河流域特色保護類村莊種群地理變異特征的理論認識。

參考文獻

［1］ 周揚， 黃晗， 劉彥隨. 中國村莊空間分布規(guī)律及其影響因素［J］. 地理學報， 2020， 75（10）： 2206-2223.

ZHOU Y， HUANG H， LIU Y S. The spatial distribution characteristics and influencing factors of Chinese villages［J］. Acta Geographica Sinica， 2020， 75（10）： 2206-2223.

［2］ 宋玥琳. 特色保護類村莊保護與發(fā)展的博弈與新平衡：以貴陽馬頭寨為例［J］. 小城鎮(zhèn)建設(shè)， 2023， 41（8）： 114-119.

SONG Y L. Game and new balance between protection and development of characteristic protected villages： Taking Guiyang Matouzhai as an example［J］.Development of Small Cities ＆ Towns， 2023， 41（8）： 114-119.

［3］ 梁妍. 延慶非典型特色保護類鄉(xiāng)村旅游適宜性評價研究：以慈母川村為例［D］. 北京： 北京建筑大學， 2021： 49-50.

［4］ 馮德東， 馮應(yīng)斌， 趙玲瓏， 等. 特色保護類村域農(nóng)村居民點空間重構(gòu)探析［J］. 地理科學進展， 2022， 41（5）： 825-837.

FENG D D， FENG Y B， ZHAO L L， et al. Spatial restructuring of rural residential land in characteristic protection villages［J］.Progress in Geography， 2022， 41（5）： 825-837.

［5］ 高溪. 鄉(xiāng)村振興戰(zhàn)略背景下特色保護類村莊空間發(fā)展策略研究：以北京市巴園子村為例［D］. 北京： 北京建筑大學， 2020： 56-58.

［6］ 王茜茜. 鄉(xiāng)村振興背景下特色保護類村莊的實用性村莊規(guī)劃研究：以山西省介休市南莊村為例［D］. 蘭州： 蘭州大學， 2021： 29-31.

［7］ 房妍杉. 特色保護類村莊景觀風貌規(guī)劃設(shè)計：以呼和浩特市石人灣村為例［D］. 呼和浩特： 內(nèi)蒙古農(nóng)業(yè)大學， 2022： 23-26.

［8］ 芮旸， 楊坤， 李宜峰， 等. 中國特色保護類村莊時空動態(tài)特征與振興對策研究：基于種群生態(tài)學視角［J］. 地理研究， 2022， 41（8）： 2203-2220.

RUI Y， YANG K， LI Y F， et al. Spatio-temporal dynamic characteristics and revitalization strategies of characteristic protection villages in China：Based on the perspective of population ecology［J］.Geographical Research， 2022， 41（8）： 2203-2220.

［9］ 汪芳，安黎哲，黨安榮，等.黃河流域人地耦合與可持續(xù)人居環(huán)境［J］.地理研究，2020，39（8）：1707-1724.

WANG F， AN L Z， DANG A R， et al. Human-land coupling and sustainable human settlements in the Yellow River Basin［J］. Geographical Research， 2020， 39（8）： 1707-1724.

［10］汪芳，苗長虹，劉峰貴，等.黃河流域人居環(huán)境的地方性與適應(yīng)性：挑戰(zhàn)和機遇［J］.自然資源學報，2021， 36（1）：1-26.

WANG F， MIAO C H， LIU F G， et al.The locality and adaptability of human settlements in the Yellow River Basin：Challenges and opportunities［J］.Journal of Natural Resources， 2021， 36（1）： 1-26.

［11］GAO C G， WU Y， BIAN C L， et al. Spatial characteristics and influencing factors of Chinese traditional villages in eight provinces the Yellow River flows through［J］. River Research and Applications， 2023， 39（7）： 1255-1269.

［12］QIAO J J， LEE J， YE X Y. Spatiotemporal evolution of specialized villages and rural development： A case study of Henan Province， China［J］. Annals of the American Association of Geographers， 2016， 106（1）： 57-75.

［13］薛明月， 王成新， 竇旺勝， 等. 黃河流域傳統(tǒng)村落空間分布特征及其影響因素研究［J］. 干旱區(qū)資源與環(huán)境， 2020， 34（4）： 94-99.

XUE M Y， WANG C X， DOU W S， et al. Spatial distribution characteristics of traditional villages in the Yellow River Basin and influencing factors［J］.Journal of Arid Land Resources and Environment， 2020， 34（4）： 94-99.

［14］邵秀英， 劉亞玲， 王向東， 等. 黃河流域傳統(tǒng)村落旅游響應(yīng)度及影響因素研究［J］. 干旱區(qū)資源與環(huán)境， 2021， 35（6）： 200-208.

SHAO X Y，LIU Y L，WANG X D， et al. Tourism responsiveness of traditional villages in the Yellow River Basin［J］.Journal of Arid Land Resources and Environment， 2021， 35（6）： 200-208.

［15］WANG F， ZHAO X G， QIU Y X， et al. Adaptability of traditional villages as tourist destinations in Yellow River Basin， China［J］. Indoor and Built Environment， 2023， 32（3）： 574-589.

［16］李如友， 石張宇. 黃河流域傳統(tǒng)村落與非物質(zhì)文化遺產(chǎn)的空間關(guān)系及形成機理［J］. 經(jīng)濟地理， 2022， 42（8）： 205-212.

LI R Y，SHI Z Y.Spatial relationship and formation mechanism of traditional villages and intangible cultural heritage in the Yellow River Basin［J］. Economic Geography， 2022， 42（8）： 205-212.

［17］裴韜， 劉亞溪， 郭思慧， 等.地理大數(shù)據(jù)挖掘的本質(zhì)［J］. 地理學報， 2019， 74（3）： 586-598.

PEI T， LIU Y X， GUO S H， et al.Principle of big geodata mining［J］. Acta Geographica Sinica， 2019， 74（3）： 586-598.

［18］唐蓓佩， 王興， 馬滕， 等.國家園林城市空間分布的點格局分析［J］. 山東農(nóng)業(yè)大學學報（自然科學版）， 2018， 49（6）： 1073-1078.

TANG B P， WANG X， MA T，et al.Point pattern analysis of national garden city spatial distribution［J］.Journal of Shandong Agricultural University（Natural Science Edition）， 2018， 49（6）： 1073-1078.

［19］田磊， 史冰心， 孫鳳芝，等. 黃河流域傳統(tǒng)村落與非物質(zhì)文化遺產(chǎn)空間相關(guān)性及其影響因素［J］. 干旱區(qū)資源與環(huán)境， 2023， 37（3）： 186-194.

TIAN L， SHI B X，SUN F Z， et al. Spatial correlation between traditional villages and intangible cultural heritage in the Yellow River Basin［J］.Journal of Arid Land Resources and Environment， 2023， 37（3）： 186-194.

［20］王勁峰， 徐成東. 地理探測器： 原理與展望［J］. 地理學報， 2017， 72（1）： 116-134.

WANG J F， XU C D. Geodetector： Principle and prospective［J］. Acta Geographica Sinica， 2017， 72（1）： 116-134.

［21］楊玉歡，賀建雄， 張新紅，等. 中國鄉(xiāng)土文化多尺度空間分布特征及影響因素［J］. 西北大學學報（自然科學版）， 2023， 53（4）： 529-540.

YANG Y H， HE J X， ZHANG X H， et al. Multi-scale spatial distribution characteristics and influencing factors of Chinese vernacular culture［J］. Journal of Northwest University（Natural Science Edition）， 2023， 53（4）： 529-540.

［22］孫儒泳，李博，諸葛陽，等.普通生態(tài)學［M］.北京：高等教育出版社，1993：11-12.

［23］楊坤， 芮旸， 李宜峰， 等. 基于共生理論的中國特色保護類村莊振興類型細分研究［J］. 地理科學進展， 2021， 40（11）： 1861-1875.

YANG K，RUI Y， LI Y F， et al.Revitalization type subdivision of characteristic protection villages in China based on the symbiosis theory［J］.Progress in Geography， 2021， 40（11）： 1861-1875.

［24］葛美玲， 封志明. 中國人口分布的密度分級與重心曲線特征分析［J］. 地理學報， 2009， 64（2）： 202-210.

GE M L， FENG Z M. Population distribution of China based on GIS： Classification of population densities and curve of population gravity centers［J］.Acta Geographica Sinica， 2009， 64（2）： 202-210.

［25］趙宏波， 魏甲晨， 苗長虹， 等. 黃河流域歷史文化名城名鎮(zhèn)名村的空間分異與影響因素分析［J］. 干旱區(qū)資源與環(huán)境， 2021， 35（4）： 70-77.

ZHAO H B， WEI J C， MIAO C H， et al. Spatial differentiation of historical and cultural cities towns and villages in Yellow River Basin［J］. Journal of Arid Land Resources and Environment， 2021， 35（4）： 70-77.

（編 輯 張 歡）

基金項目：國家自然科學基金（42201222）；教育部人文社會科學研究青年基金（23XJCZH013）。

第一作者：王玲， 女， 從事鄉(xiāng)村地理與規(guī)劃研究， wanglingl15@126.com。

通信作者：芮旸， 男， 副教授，從事城鄉(xiāng)轉(zhuǎn)型與空間治理研究， ruiyang@nwu.edu.cn。