摘 要：【目的】基于多期土地利用/土地覆蓋專題數(shù)據(jù)，探討和分析1980—2020年河南省的人類活動強度時空變化特征?！痉椒ā坎捎藐懙乇韺尤祟惢顒訌姸饶Ｐ陀嬎悴煌瑫r期河南省的人口活動強度，基于全局自相關和局部自相關分析河南省40年內人類活動強度的時空變化特征?！窘Y果】2020年，河南省土地利用/土地覆蓋類型主要為耕地和林地，分別占全省總面積的64.85%和16.58%。40年內建設用地面積大幅度增加54.17%，林地和水域面積分別增加1.65%和3.6%，其余地類面積均減少，其中耕地、未利用土地和草地分別減少5 668.93 km2、134.85 km2和1 174.27 km2。1980年河南省人類活動強度為22.09%，2020年達到25.46%?！窘Y論】1980—2020年河南省人類活動強度空間聚集特征顯著，空間分布格局呈東高西低、北高南低的特點。人類活動強度變化冷點地區(qū)主要集中在豫西山地、豫北山地、豫南山地等地區(qū)，熱點地區(qū)分布在以鄭州為中心的河南省中東部地區(qū)。

關鍵詞：人類活動強度；土地利用/土地覆被；河南省

中圖分類號：X312? ? ? 文獻標志碼：A? ? ?文章編號：1003-5168（2024）07-0096-06

DOI：10.19968/j.cnki.hnkj.1003-5168.2024.07.019

Temporal and Spatial Changes and Characteristics Analysis of Human Activity Intensity in Henan Province from 1980 to 2020

MA Qiming

（College of Land Engineering， Chang'an University， Xi'an 710064， China）

Abstract：[Purposes] The study is based on multi-temporal thematic data of land use/land cover， aiming to explore and analyze the temporal and spatial variations of human activity intensity in Henan Province， China from 1980 to 2020. [Methods] Utilizing the land surface human activity intensity model， this study calculates the human activity intensity in Henan Province during different periods. Subsequently， based on global autocorrelation and local autocorrelation， the spatial and temporal variation characteristics of human activity intensity in Henan Province in 40 years were analyzed. [Findings] In 2020， the land use/land cover types in Henan Province were primarily farmland land and forest land， accounting for 64.85% and 16.58% of the total area， respectively. Over the course of forty years， the area of construction land has significantly increased by 54.17%， while forest and water areas have increased by 1.65% and 3.6%， respectively. Conversely， the remaining land classes have experienced a decrease， with farmland， unused land， and grassland decreasing by 5 668.93 km?， 134.85 km?， and 1 174.27 km?， respectively. The human activity intensity in Henan Province was 22.09% in 1980 and has reached 25.46% in 2020. [Conclusions] From 1980 to 2020， the spatial aggregation characteristics of human activity intensity in Henan Province were significant， displaying an eastern high-western low and northern high-southern low distribution pattern. Areas with lower fluctuations in human activity intensity were primarily concentrated in regions such as the western mountains， northern mountains， and southern mountains of Henan. On the other hand， regions exhibiting higher fluctuations were concentrated in the central and eastern parts of Henan Province， centered around Zhengzhou.

Keywords： human activity intensity; land use/land cover; Henan Province

0 引言

對一定區(qū)域內的人類活動強度進行定量分析及空間化表達，可為該區(qū)域后期的土地利用政策制定和生態(tài)建設提供理論參考和科學依據(jù)。人類活動強度是指一定面積的區(qū)域內受人類活動所產(chǎn)生的擾動程度，是人類對該區(qū)域利用、改造和開發(fā)的程度［1-2］。當前，國內外學者對人類活動強度的研究主要集中于人類活動強度量化方法、時空演變及驅動因素分析等方面。隨著遙感和地理信息系統(tǒng)技術的蓬勃發(fā)展，高質量人類活動強度空間數(shù)據(jù)得以擴增和豐富［3］。人類活動強度量化方法可劃分為非空間化方法和空間化方法。非空間化量化方法包括指標評價法［4］和統(tǒng)計分析法［5］等?？臻g化量化方法又可分成間接空間化和直接空間化兩類：間接空間化是基于可用于衡量指標的變化來定量評估人類活動強度，如凈初級生產(chǎn)力人類占用法［6］、全球干擾指數(shù)法［7］和植被指數(shù)法［8］等；直接空間化方法包括基于自然、社會和經(jīng)濟等多方面指標來定量評價人類活動強度的綜合指標法［9］、基于土地利用轉變、人口密度等空間數(shù)據(jù)的人類足跡指數(shù)法（Human Footprint Index ，HFI）［10］及利用統(tǒng)計數(shù)據(jù)和遙感解譯數(shù)據(jù)的土地類型變化法等［2］。土地類型變化法因普適性較強和人類活動強度空間數(shù)據(jù)充分而被廣泛使用，其中徐勇等提出的陸地表層人類活動強度（Human Activity Intensity of Land Surface，HAILS）算法模型已應用于黃土高原、黔南地區(qū)和海南島等不同空間尺度并取得良好成果［11-13］。本研究基于多期土地利用/土地覆蓋專題數(shù)據(jù)，探討和分析1980—2020年河南省的人類活動強度時空變化特征，以期為今后的土地利用政策制定和生態(tài)建設提供理論參考和科學依據(jù)。

1 研究區(qū)和數(shù)據(jù)來源

河南省位于我國中東部，地理位置介于東經(jīng)110°21′～116°39′，北緯31°23′～36°22′。在黃河中下游地區(qū)，是我國礦產(chǎn)資源大省、人口大省及農(nóng)產(chǎn)品主產(chǎn)區(qū)。河南省包括17個地級市、1個省直轄縣級市和83個縣區(qū)，土地總面積為16.7萬km2。

本研究所采用的土地利用/土地覆被數(shù)據(jù)來源于中國科學院資源環(huán)境數(shù)據(jù)中心徐新良等學者的中國多時期土地利用遙感監(jiān)測數(shù)據(jù)集（CNLUCC）。該數(shù)據(jù)集是以美國陸地衛(wèi)星Landsat遙感影像為主要數(shù)據(jù)源，通過人工目視解譯構建的我國土地利用/土地覆蓋專題數(shù)據(jù)庫，空間分辨率為1 km。該數(shù)據(jù)集的分類系統(tǒng)為二級分類系統(tǒng)，其中一級分類包括耕地、草地、水域、林地、未利用土地和城鄉(xiāng)、工礦、居民用地共6個類型，二級分類基于土地資源的自然屬性劃分為23個類型［13］。本研究使用的是CNLUCC數(shù)據(jù)集中1980、1990、2000、2005、2010、2015、2020年共7期數(shù)據(jù)。

2 研究方法

2.1 陸地表層人類活動強度

河南省不同時期人類活動強度量化的算法模型采用的是徐勇等提出的陸地表層人類活動強度（HAILS）和建設用地當量折算系數(shù)（Conversion Index of Construction Land Equivalent，CI）［2］。不同的土地利用/土地覆被類型表征著人類活動對陸地表層自然覆被利用、改造和開發(fā)的不同程度?；诙蛹壏植际剿惴蓪⒉煌牡仡悡Q算為不同的建設用地當量折算系數(shù)。陸地表層人類活動強度的計算公式為式（1）［2］。

[HAILS=SCLES×100%][SCLE=i=1nSLi?CIi] （1）

式中：HAILS為研究區(qū)內人類活動強度；SCLE為研究區(qū)域內建設用地當量面積；S為研究區(qū)總面積；SLi為研究區(qū)內第i種土地利用/覆被類型的面積；CIi為第i種土地利用/覆被類型的建設用地當量折算系數(shù)；n為研究區(qū)內土地利用/覆被類型總數(shù)。

2.2 空間自相關

空間自相關是同一變量在不同地理空間位置的相關性，是空間單元屬性值聚集程度的一種度量［14］。它包括全局空間自相關和局部空間自相關。因此，本研究引入Moran's I指數(shù)去定量分析河南省人類活動強度的空間自相關，從而探尋該區(qū)域人類活動強度的空間分布模式。同時引入局部G系數(shù)來計算省內人類活動強度的局部空間自相關，從而識別研究區(qū)內的人類活動強度的高、低值聚集分布。

2.2.1 全局空間自相關。利用Moran's I指數(shù)計算河南省人類活動強度進行全局空間自相關，其表達式為式（2）。

[I=i=1nj=1nWijxi-xxj-xS2?i=1nj=1nWij] （2）

式中：S2表示研究區(qū)內人類活動強度的方差；n表示研究區(qū)空間單元總數(shù)；xi和xj表示第i空間單元和第j空間單元的人類活動強度，`x表示河南省人類活動強度平均值，Wij表示第i空間單元和第j像元的空間權重矩陣。其中，當Moran's I為正時，表示人類活動強度呈現(xiàn)空間正相關，值越大表明空間相關性越明顯；當Moran's I為負時，表示人類活動強度呈現(xiàn)空間負相關，值越小則空間差異越大；Moran's I為0時，則表示變量是隨機分布。使用Z得分來進行顯著性檢驗，若|Z|>1.96，則說明檢驗結果具有顯著性。

2.2.2 局部空間自相關

利用局部G系數(shù)計算河南省人類活動強度的局部空間自相關性［15］，其表達式為式（3）。

[Gi=j=1nWijxjj=1nxj] （3）

式中，Gi表示第i空間單元的局部G系數(shù)值；Wij表示第i空間單元和第j像元的空間權重矩陣。若G為正時，則表明人類活動強度較高地區(qū)空間集聚，屬于熱點區(qū)域，z值越高，熱點聚類越緊密；若G為負時，則表明人類活動強度較低地區(qū)空間集聚，屬于冷點區(qū)域，z值越低，冷點聚類越緊密。

3 結果

3.1 河南省土地利用變化時空特征分析

基于多期土地利用/土地覆蓋數(shù)據(jù)計算1980—2020年河南省土地利用/土地覆被一級類型構成及變化見表1。由表1可知，2020年河南省土地利用類型主要以耕地為主，耕地面積占河南省面積的64.85%，其中水田類型和旱地類型分別占河南省面積的6.41%和58.44%，其面積分別為10 613.71km2和96 740.37km2，水田主要分布于豫中的新鄉(xiāng)市、開封市及豫南的信陽市，而旱地則分散于省內平原及丘陵地區(qū)。其次是林地類型，占河南省面積的16.58%，其中灌木林、疏林地、有林地和其他林地的面積分別為3 459.94 km2、5 624.03 km2、17 535.41 km2和581.17 km2，主要分布在西部、北部和南部山區(qū)丘陵地帶。再次是城鄉(xiāng)/工礦/居民用地面積總計18 222.75 km2，占河南省面積的11.01%。草地類型面積為8 589.12 km2，占河南省面積的5.19%。未利用土地類型占地面積最小，僅占河南省面積的0.01%。

1980—2020年，未利用土地類型占地面積比例整體變化幅度最大，減少了134.85 km2，減幅為89.34%。其次是城鄉(xiāng)/工礦/居民用地，其面積從1980年的11 819.86 km2增大到2020年的18 222.75 km2，增幅為54.17%。草地和耕地面積減少了-12.03%和-5.02%，減少面積為1 174.27 km2和5 668.93 km2。林地和水域面積增加了446.67 km2和135.75 km2，分別占1980年的1.65%和3.6%。

3.2 河南省人類活動強度的時空變化分析

1980—2020年河南省建設用地當量面積及人類活動強度變化如圖1所示。其中，前期1980—1990年增長趨勢較為緩慢，中期1990—2015年保持較平穩(wěn)的增長趨勢，后期2015—2020年增長趨勢則相對迅速。基于陸地表層人類活動強度模型計算可得，1980年河南省建設用地當量面積為36 567.91 km2，人類活動強度為22.09%，2000年河南省建設用地當量面積為37 927.96 km2，人類活動強度為22.91%，2020年河南省建設用地當量面積和人類活動強度分別為42 139.02 km2和25.46%。從1980年至2020年，河南省的建設用地當量面積和人類活動強度分別增加了5 571.11 km2和3.36%。

1980—2020年河南省建設用地當量面積和人類活動強度變化如圖2所示，各市人類活動強度變化見表2。由圖2、表2可知，河南省人類活動強度較低的區(qū)域主要集中在豫西山地、豫北山地、豫南山地處，而人類活動強度高的區(qū)域分布在河南中部和東部區(qū)域。河南省西高東低并呈望北向南之勢，北、西、南三面山脈丘陵居多，而中東部多是廣袤平原、肥沃田地，因此歷年人類活動強度的整體空間分布也呈半環(huán)狀。濟源市、洛陽市、南陽市、三門峽市及信陽市的人類活動強度均較低（＜22%），這些區(qū)域位于豫西山地、豫北山地、豫南山地等地區(qū)，域內山脈丘陵較多，自然稟賦條件薄弱，因此耕地和建設用地所占面積比例較小。而濮陽市、新鄉(xiāng)市、開封市和周口市等區(qū)域處于河南自然本地條件較好的平原區(qū)域，耕地覆蓋率高且人口聚集，因此這些區(qū)域具有較高的人類活動強度。其中鄭州市的人類活動強度自1980年的26.51%增至2020年的38.56%，40年內增幅為12.04%，遠超省內其他區(qū)域人類活動強度增長幅度。這可能與當?shù)氐娜丝诹魅?、?jīng)濟發(fā)展和政府政策等驅動因素有關。

3.3 河南省人類活動強度的空間自相關分析

基于空間自相關的全局指標Moran's I指數(shù)可得出1980—2000年以來歷年河南省人類活動強度全局空間自相關指標結果見表3。由表3可知，歷年Moran's I指數(shù)值均為正，Z得分大于1.96，P值為0，表明歷年河南省的人類活動強度呈現(xiàn)顯著的空間集聚特征。同時自1990年起，每一年的Moran's I指數(shù)處于上升趨勢，表明1990—2020年河南省人類活動強度的空間集聚程度呈現(xiàn)加強趨勢。1980—2000年、2000—2020年及1980—2020年的Moran's I指數(shù)值為0.21、0.04和0.06，這表明這三個時間段內的人類活動強度變化亦呈現(xiàn)空間集聚特征。

使用ArcGIS10.8軟件上的熱點分析工具進一步探尋人類活動強度及其變化的空間差異。地理空間位置上相互鄰近且人類活動強度高于研究區(qū)平均水平的單元構成高值集聚區(qū)，亦稱為熱點地區(qū)；反之地理空間位置上相互鄰近且人類活動強度低于研究區(qū)平均值的地區(qū)構成低值集聚區(qū)，亦稱為冷點地區(qū)［13］。980—2020年河南省人類活動強度變化熱點如圖3所示。由圖3可知，從空間分布上來看，人類活動強度的熱點地區(qū)主要離散分布于河南省的中部和東部地區(qū)的各市域內，冷點地區(qū)集中分布于三門峽市、濟源市、洛陽市南部、南陽市北部、焦作市北部邊緣和信陽市南部的山區(qū)。從時空格局變化上分析，1980—2000年，河南省的人類活動強度變化星星狀分散于河南省內中北部區(qū)域；2000—2020年，人類活動強度變化逐漸集中于鄭州市為中心的中東部區(qū)域。

4 結論與討論

1980—2020年期間，河南省土地利用/土地覆蓋類型主要為耕地類型，然后依次是林地、建設用地、草地、水域和未利用地類型。40年內，隨著河南省社會經(jīng)濟迅速發(fā)展和人口數(shù)量的增長，人們對陸地表層覆被利用、改造和開發(fā)的程度逐漸加深，建設用地面積呈現(xiàn)大幅度增長趨勢，林地和水域面積出現(xiàn)微弱增長，其他土地利用/土地覆被類型均出現(xiàn)不同程度的縮減，其中未利用土地面積減少幅度最大，占1980年此類型土地的89.34%；耕地類型減少面積最大，共減少5 668.93 km2。

40年內，建設用地當量面積和人類活動強度保持連續(xù)增長趨勢，其中變化以2010年為界，1980—2010年期間，人類活動強度增長相對緩慢；2010—2020年，增長迅速。40年內省內人類活動強度共增加3.36%。

1980—2020年，河南省的人類活動強度的空間分布格局呈現(xiàn)出東高西低、北高南低特點。研究區(qū)內人類活動強度具有顯著的空間集聚特征，冷點地區(qū)主要分布在豫西山地、豫北山地、豫南山地等山脈丘陵地區(qū)，熱點地區(qū)則分布在豫中、豫東等的平原區(qū)域。河南省人類活動強度變化的熱點主要集中在以鄭州市為中心的研究區(qū)中東部區(qū)域，即城市建設和經(jīng)濟發(fā)展迅速地區(qū)。本研究對河南省40年內的人類活動強度進行定量分析及空間化表達，可為今后的土地利用政策制定和生態(tài)建設提供理論參考和科學依據(jù)。

參考文獻：

［1］劉世梁， 劉蘆萌， 武雪， 等. 區(qū)域生態(tài)效應研究中人類活動強度定量化評價［J］. 生態(tài)學報， 2018， 38（19）： 6797-6809.

［2］徐勇， 孫曉一， 湯青. 陸地表層人類活動強度：概念、方法及應用［J］. 地理學報， 2015， 70（7）： 1068-1079.

［3］段群滔， 羅立輝. 人類活動強度空間化方法綜述與展望：以青藏高原為例［J］. 冰川凍土， 2021， 43（5）： 1582-1593.

［4］文英. 人類活動強度定量評價方法的初步探討［J］. 科學對社會的影響， 1998（4）： 56-61.

［5］李子君， 李秀彬. 近45年來降水變化和人類活動對潮河流域年徑流量的影響［J］. 地理科學，2008，28（6）： 809-813.

［6］VITOUSEK P M， EHRLICH P R， EHRLICH A H， et.al. Human appropriation of the products of photosynthesis［J］. BioScience， 1986， 36（6）： 368-373.

［7］MILDREXLER D J， ZHAO M， RUNNING S W. Testing a MODIS global disturbance index across North America［J］. Remote Sensing of Environment， 2009， 113（10）： 2103-2117.

［8］EVANS J， GEERKEN R. Discrimination between climate and human-induced dryland degradation［J］. Journal of arid environments， 2004， 57（4）： 535-554.

［9］HU Z B， HE X Y， LI Y H， et.al. Human activity intensity and its spatial distribution pattern in upper reach of Minjiang River［J］. Chinese Journal of Ecology，2007，26（4）：539.

［10］WOOLMER G. The human footprint and the last of the wild［J］. Bioscience， 2002， 52（10）：891-904.

［11］ 劉慧明， 劉曉曼， 李靜， 等. 生物多樣性保護優(yōu)先區(qū)人類干擾遙感監(jiān)測與評價方法［J］. 地球信息科學學報， 2016，18（8）： 1103-1109.

［12］ 徐小任， 徐勇. 黃土高原地區(qū)人類活動強度時空變化分析［J］. 地理研究， 2017， 36（4）： 661-672.

［13］ 劉采， 張海燕， 李遷. 1980—2018年海南島人類活動強度時空變化特征及其驅動機制［J］. 地理科學進展， 2020，39（4）：567-576.

［14］ TOBLER W R. A Computer Movie Simulating Urban Growth in the Detroit Region［J］. Economic Geography， 1970，46：234-240.

［15］ORD J K， GETIS A. Local Spatial Autocorrelation Statistics： Distributional Issues and an Application［J］. Geographical Analysis， 2010， 27（4）： 286-306.

參考文獻

[1] 劉世梁， 劉蘆萌， 武雪， 等. 區(qū)域生態(tài)效應研究中人類活動強度定量化評價[J]. 生態(tài)學報， 2018， 38（19）： 6797-6809.

[2] 徐勇， 孫曉一， 湯青. 陸地表層人類活動強度：概念、方法及應用[J]. 地理學報， 2015， 70（7）： 1068-1079.

[3] 段群滔， 羅立輝. 人類活動強度空間化方法綜述與展望：以青藏高原為例[J]. 冰川凍土， 2021， 43（5）： 1582-1593.

[4] 文英. 人類活動強度定量評價方法的初步探討[J]. 科學對社會的影響， 1998（4）： 56-61.

[5] 李子君， 李秀彬. 近45年來降水變化和人類活動對潮河流域年徑流量的影響[J]. 地理科學， 2008， 28（6）： 809-813.

[6] VITOUSEK P M， EHRLICH P R， EHRLICH A H， et.al. Human appropriation of the products of photosynthesis[J]. BioScience， 1986， 36（6）： 368-373.

[7] MILDREXLER D J， ZHAO M， RUNNING S W. Testing a MODIS global disturbance index across North America[J]. Remote Sensing of Environment， 2009， 113（10）： 2103-2117.

[8] EVANS J， GEERKEN R. Discrimination between climate and human-induced dryland degradation[J]. Journal of arid environments， 2004， 57（4）： 535-554.

[9] HU Z B， HE X Y， LI Y H， et.al. Human activity intensity and its spatial distribution pattern in upper reach of Minjiang River[J]. Chinese Journal of Ecology， 2007， 26（4）： 539.

[10] WOOLMER G. The human footprint and the last of the wild[J]. Bioscience， 2002， 52（10）： 891904Santelmann.

[11] 劉慧明， 劉曉曼， 李靜， 等. 生物多樣性保護優(yōu)先區(qū)人類干擾遙感監(jiān)測與評價方法[J]. 地球信息科學學報， 2016， 18（8）： 1103-1109.

[12] 徐小任， 徐勇. 黃土高原地區(qū)人類活動強度時空變化分析[J]. 地理研究， 2017， 36（4）： 661-672.

[13] 劉采， 張海燕， 李遷. 1980—2018年海南島人類活動強度時空變化特征及其驅動機制[J]. 地理科學進展， 2020， 39（4）： 567-576.

[14] TOBLER W R. A Computer Movie Simulating Urban Growth in the Detroit Region[J/OL]. Economic Geography， 1970， 46： 234.

[15] ORD J K， GETIS A. Local Spatial Autocorrelation Statistics： Distributional Issues and an Application[J]. Geographical Analysis， 2010， 27（4）： 286-306.

收稿日期：2023-10-19

基金項目：雅礱江流域水電開發(fā)有限公司科研試驗項目“高寒高海拔光伏電站生態(tài)效應評估與護修對策研究”（000023-23XB0065）; 內蒙古自治區(qū)科技重大專項項目“現(xiàn)代牧區(qū)草地高效生產(chǎn)與家畜優(yōu)化利用技術研發(fā)與示范”（2021ZD0045）。

作者簡介：馬起明（1998—），女，碩士生，研究方向：地理時空統(tǒng)計方法。