人类的足迹遍布整个星球,地球进入了人类世(Kareiva et al., 2007; Steffen et al., 2007)。然而,地球的承载能力是有限的。随着对气候变化和环境危机理解的加深,人类与自然之间的矛盾被认为是阻碍可持续发展的根本原因(Anderson et al., 2019; Dorninger et al., 2017)。人类活动在生物圈完整性以及氮和磷的循环方面已超出了行星边界,生活在一个安全和公正空间的挑战依然巨大(O’Neill et al., 2018)。因此,理解社会经济发展与生态环境演变之间的关系对于全球可持续发展至关重要。此外,平衡人类活动与生态环境是许多联合国倡议的核心。
写作 - 原始草稿:Haimeng Liu, Yi Cheng
写作 - 审阅与编辑:刘海萌,刘志峰,李启瑞,张海燕
可持续发展目标(SDGs),包括目标 6、7、11、12、13、14 和 15(Cheng et al., 2021;Liu, Hull, et al., 2018),以及许多全球组织或项目也探讨了人类-自然系统的耦合机制和协同发展方法(Kates et al., 2001;Liu, Fang, & Fang, 2020)。分析人类活动与生态环境之间的矛盾与协调对人类的命运和福祉具有重要的现实意义。
青藏高原(QTP)平均海拔超过 4000 米,广泛分布着冰川和永久冻土,是亚洲九大河流的源头,为超过 15 亿人提供淡水、食物和其他生态系统服务,被称为地球的第三极和亚洲水塔(Yao et al., 2012)(图 1)。它在全球气候调节和生物多样性维持中发挥着重要作用(Li et al., 2020)。然而,由于其高海拔、寒冷的气温和干燥的环境,QTP 的生态系统对人类活动极为敏感和脆弱(Han et al., 2022;Liu, Milne, et al., 2018)。从 1980 年到 2020 年,人口增加了 1.56 倍,城市建设面积增加了近三倍,高速公路长度增加了四倍以上(Zhang et al., 2019)。过度放牧等强烈的人类活动对当地草原生态系统施加了巨大的压力,这种压力可能远远超过气候变化的影响(Wei et al., 2022)。人类活动的加剧对 QTP 生态系统的稳定性和可持续发展目标(SDGs)的实现构成了重大威胁(Yang et al., 2022)。 因此,阐明 QTP 上的人类与自然的关系对地方和全球可持续发展具有重要意义。
人类与自然之间的互动自人类首次出现以来就存在(Costanza et al., 2007)。在环境问题日益严重的背景下,许多来自不同学科背景的学者专注于这一与可持续发展密切相关的领域,并产生了宝贵的见解(Harden, 2012)。研究主题包括人类活动对自然的正面和负面影响(Kareiva et al., 2007)、自然资源和生态系统对社会发展的约束和承载能力(Fan et al., 2017),以及它们之间的互动、反馈、权衡和适应(Comberti et al., 2015;Holling & Gunderson, 2002;Kareiva et al., 2007)。一些研究表明,减少人类活动的规模对于缓解社会危机是必要的(Dirzo et al., 2022);另一些研究则表明,如果我们在生产方式、能源使用、生态保护等方面进行适当的改变,可以实现双赢的愿景(Cao et al., 2009;Tallis et al., 2018)。
在这个研究领域,城市化地区(Fang et al., 2016)、流域(Yin et al., 2021)、森林(Davidson et al., 2012)和沿海地区(Yi et al., 2018)受到了更多关注。对于自然系统,研究人员关注生物多样性(Cardinale et al., 2012)、植被(Zhang, Yang, et al., 2022)、空气污染(Liu, Cui, & Zhang, 2022)、河流环境(Leprieur et al., 2008)、土壤侵蚀(Wei et al., 2006)等。土地利用和覆盖变化(Cheng et al., 2022)、人口密度(Cropper & Griffiths, 1994)和夜间光照(Cai et al., 2021)通常被用作人类活动的代理变量;一些学者还构建了一个名为人类足迹的综合指数(Li, Wu, Gong, et al., 2018;Venter et al., 2016)。用于分析人类与自然之间复杂关系的方法包括网络方法(Kluger et al., 2020)、系统动力学(Reed et al., 2022;Yoon et al., 2021)、耦合协调分析(Jiang et al., 2022;Shi, Feng, et al., 2022)、复杂系统理论(Bennett & McGinnis, 2008;Wang & Grant, 2021)、环境库兹涅茨曲线模型(Zhou et al.,2015),元耦合框架(Liu,2017),供需平衡分析(Liu,Xing 等,2022),问卷调查(Allison 等,2021;Liu,Zhang 等,2020),以及定性分析(Ives 等,2017)。这些方法从不同角度推动了复杂人类-自然系统的研究,但在细网格尺度上面临重大限制。某些使用栅格数据的研究:Vačkár 等(2012)使用 10 公里网格分析了捷克共和国人口密度、土地利用强度和生物多样性之间的空间关系。Shi,Feng 等(2022)利用耦合协调度阐明了中国 10kmxx10km10 \mathrm{~km} \times 10 \mathrm{~km} 地区景观生态风险与城市化之间的协调关系。Zhou 等(2023)测量了中国城市中城市化与植被覆盖之间的冲突与协调动态,使用了 250 米的植被指数。Simkin 等(2022)预测了城市土地扩张对三种共享社会经济路径(SSP)情景下生物多样性的未来影响。然而,这些研究集中于某些人类活动或自然元素,缺乏对人类和自然系统的整体考虑。
古老的人类自 16 万年前起就占据了青藏高原(Chen et al., 2019)。在过去的半个世纪里,人类活动对青藏高原生态系统造成了有限但迅速增加的干扰,特别是在草地、宝贵的水源保护区和生物多样性保护区的最大比例上(Jie et al., 2015;Li, Zhang, Wang, et al., 2018)。关于青藏高原上人类与自然关系的研究
图 1. 青藏高原的自然地理。
主要关注城市化、土地利用、基础设施建设、农业和畜牧业对生态环境的影响(Miao et al., 2021; Wang et al., 2015; Zhang, Zhang, et al., 2022)。这里的城市化水平相对较低,但近年来城市人口迅速增长,建设用地快速扩张(Wang et al., 2020)。尽管城市化会导致更多的污染物排放,但通过将人口从农业和牧区转移,可以减少人类活动对生态系统造成的整体干扰(Tian et al., 2021)。2000 年至 2015 年间,城市化与生态环境之间的协调性呈上升趋势(Feng & Li, 2021)。草地退化主要是由于短期内人类活动造成的,这进一步导致生物多样性丧失、土壤侵蚀加剧、碳汇减少以及当地人类福祉下降(Chen et al., 2020; Dong et al., 2020)。青藏铁路和川藏铁路的建设和运营将给当地生态环境带来潜在风险(Cui et al., 2022; Luo et al., 2020; Qin & Zheng, 2010)。 从系统的角度,利用承载能力模型(Niu et al., 2020)分析了产业、人口与资源环境的共同发展,通过将人类活动与生态系统服务联系起来评估生态质量(Sun et al., 2020;Wang, Huang, et al., 2022),并通过整合生态能力、生态足迹和资源效率来衡量区域可持续性(Fan & Fang, 2022)。与此同时,中国在青藏高原的生态保护与恢复项目,包括植树造林、草原恢复和设立自然保护区,正在改善许多地区的生态环境(Fu et al., 2021;Jin et al., 2019)。近年来,矿业和城市的污水排放和废气排放、因游客增加而产生的垃圾以及光污染已成为新的严重环境问题(Chen, Zhang, et al., 2021;Feng et al., 2020;Wang, Lv, et al., 2022)。
存在许多度量标准来量化人类活动的强度(人类足迹,人类改造),这些标准常用于描述人类对生态系统的压力(Venter et al., 2016; Williams et al., 2020)。本研究构建的 HAI 在很大程度上遵循了先前研究建立的评估指标体系,包括人口密度、GDP、夜间灯光、土地覆盖和道路网络(Huang et al., 2022; Mu et al., 2022; Xu et al., 2016)。然而,考虑到草原占据了 60%60 \% 的 QTP(Li
Assigned from 0 to 100 according to intervals determined by 10 equal quantiles
(Mu et al., 2022)| Assigned from 0 to 100 according to intervals determined by 10 equal quantiles |
| :--- |
| (Mu et al., 2022) |
Data set Year and format Data source Data processing
Population density " 2000, 2005, 2010, 2015, 2020;
Raster, 1km" WorldPop For all locations with more than 1,000 people km^(2), we assigned a score of 100 . For more sparsely populated areas, we logarithmically scaled the score using 33.333 xx log (population density +1 ) (Williams et al., 2020)
GDP density " 2000, 2005, 2010, 2015, 2019;
Raster, 1km" Chen et al. (2022) Scale data into the 0-100 range using Min-Max Normalization, and the maximum value is the average of the top 50 values in the 5 years
Nighttime light " 2000, 2005, 2010, 2015, 2020;
Raster, 1km" Zhang, Ren, et al. (2021) "Assigned from 0 to 100 according to intervals determined by 10 equal quantiles
(Mu et al., 2022)"
Grazing intensity " 2000, 2005, 2010, 2015, 2019;
Raster, 250m" Liu (2021) The data in 2005 used the average of 2000 and 2010, and the data in 2015 used the average of 2010 and 2019. Resampling to 1-km resolution. Scale data into the 0-100 range using Min-Max Normalization, and the maximum value is the average of the top 50 values in the 5 years
Land use/cover " 2000, 2005, 2010, 2015, 2020;
Raster, 1km" www.resdc.cn/data. aspx?DATAID=335 Assign score: 100 for urban land, 50 in 1 km buffer; 80 for construction land, 40 in 1 km buffer; 60 for cropland, 20 in 1 km buffer; 0 for other lands
Road network " 2000, 2005, 2010, 2015, 2018;
Shapefile " Digitization by authors Assign a score of 60 for 0.5 km out for either side of a highway, national road, and provincial road. Assign 80 and 40 for 0.5 and 1.5 km out for either side of a railway (Venter et al., 2016). Convert to raster data| Data set | Year and format | Data source | Data processing |
| :---: | :---: | :---: | :---: |
| Population density | $\begin{aligned} & \text { 2000, 2005, 2010, 2015, 2020; } \\ & \text { Raster, } 1 \mathrm{~km} \end{aligned}$ | WorldPop | For all locations with more than 1,000 people $\mathrm{km}^{2}$, we assigned a score of 100 . For more sparsely populated areas, we logarithmically scaled the score using $33.333 \times \log$ (population density +1 ) (Williams et al., 2020) |
| GDP density | $\begin{aligned} & \text { 2000, 2005, 2010, 2015, 2019; } \\ & \text { Raster, } 1 \mathrm{~km} \end{aligned}$ | Chen et al. (2022) | Scale data into the $0-100$ range using Min-Max Normalization, and the maximum value is the average of the top 50 values in the 5 years |
| Nighttime light | $\begin{aligned} & \text { 2000, 2005, 2010, 2015, 2020; } \\ & \text { Raster, } 1 \mathrm{~km} \end{aligned}$ | Zhang, Ren, et al. (2021) | Assigned from 0 to 100 according to intervals determined by 10 equal quantiles <br> (Mu et al., 2022) |
| Grazing intensity | $\begin{aligned} & \text { 2000, 2005, 2010, 2015, 2019; } \\ & \text { Raster, } 250 \mathrm{~m} \end{aligned}$ | Liu (2021) | The data in 2005 used the average of 2000 and 2010, and the data in 2015 used the average of 2010 and 2019. Resampling to $1-\mathrm{km}$ resolution. Scale data into the $0-100$ range using Min-Max Normalization, and the maximum value is the average of the top 50 values in the 5 years |
| Land use/cover | $\begin{aligned} & \text { 2000, 2005, 2010, 2015, 2020; } \\ & \text { Raster, } 1 \mathrm{~km} \end{aligned}$ | www.resdc.cn/data. aspx?DATAID=335 | Assign score: 100 for urban land, 50 in 1 km buffer; 80 for construction land, 40 in 1 km buffer; 60 for cropland, 20 in 1 km buffer; 0 for other lands |
| Road network | $\begin{aligned} & \text { 2000, 2005, 2010, 2015, 2018; } \\ & \text { Shapefile } \end{aligned}$ | Digitization by authors | Assign a score of 60 for 0.5 km out for either side of a highway, national road, and provincial road. Assign 80 and 40 for 0.5 and 1.5 km out for either side of a railway (Venter et al., 2016). Convert to raster data |
rho1\rho 1 和 rho2\rho 2 是 MOD09A1 图像的蓝光和红光波段的图像反射率 (Yang et al., 2020)。重采样至 1-km1-\mathrm{km} 分辨率。使用最小-最大归一化将数据缩放到 0-1000-100 范围内。
SI=sqrt(rho_(1)xxrho_(2)), SSI=sqrt((NDVI-1)^(2)+SI^(2))
The rho1 and rho2 are the image reflectances of the blue and red light band of the MOD09A1 image (Yang et al., 2020). Resampling to 1-km resolution. Scale data into the 0-100 range using Min-Max Normalization| $\begin{aligned} & \mathrm{SI}=\sqrt{\rho_{1} \times \rho_{2}} \\ & \mathrm{SSI}=\sqrt{(\mathrm{NDVI}-1)^{2}+\mathrm{SI}^{2}} \end{aligned}$ |
| :--- |
| The $\rho 1$ and $\rho 2$ are the image reflectances of the blue and red light band of the MOD09A1 image (Yang et al., 2020). Resampling to $1-\mathrm{km}$ resolution. Scale data into the $0-100$ range using Min-Max Normalization |
Data set Year and format Data source Data processing Weight
Fractional vegetation coverage (FVC) " 2000, 2005, 2010,
2015, 2020; Raster, " 1km NDVI data: Juntao (2022) FVC =( NDVI - NDVIsoil )//( NDVIveg - NDVIsoil). Select the cumulative frequency of 5% and 95% of NDVI value as NDVIsoil and NDVIveg 0.16
Biological richness index (BRI) "2000, 2005, 2010,
2015, 2020; Raster,
1 km" Land use and cover data: www.resdc.cn BRI=Abioxx(0.35 xx forest +0.21 xx grassland +0.28 xx water area +0.11 xx cropland +0.04 xx construction land +0.01 xx unused land)/area. Abio is a normalized coefficient with a reference value of 511.264. Scale data into the 0-100 range using Min-Max Normalization 0.19
Net primary productivity (NPP) 2000, 2005, 2010, 2015, 2020; Raster, 500 m https://lpdaac.usgs. gov/products/ mod17a3hgfv061/ Resampling to 1-km resolution. Scale data into the 0-100 range using Min-Max Normalization 0.15
Soil salinization index (SSI) 2000, 2005, 2010, 2015, 2020; Raster, 500 m MOD09A1 image: https://lpdaac. usgs.gov/products/ mod09a1v006/ "SI=sqrt(rho_(1)xxrho_(2)), SSI=sqrt((NDVI-1)^(2)+SI^(2))
The rho1 and rho2 are the image reflectances of the blue and red light band of the MOD09A1 image (Yang et al., 2020). Resampling to 1-km resolution. Scale data into the 0-100 range using Min-Max Normalization" -0.14
Soil erosion 2000, 2005, 2010, 2015, 2018; Raster, 1 km https://www.resde.cn/ Scale data into the 0-100 range using Min-Max Normalization -0.13
Water conservation 2000, 2005, 2010, 2015, 2020; Raster, 1 km Wang (2022) Scale data into the 0-100 range using Min-Max Normalization 0.12
PM_(2.5) concentration 2000, 2005, 2010, 2015, 2020; Raster, 1 km https://sites.wustl.edu/ acag/datasets/surface-pm2-5/#V5.GL. 02 Scale data into the 0-100 range using Min-Max Normalization, and the maximum value is the average of the top 50 values in the 3 years -0.11| Data set | Year and format | Data source | Data processing | Weight |
| :---: | :---: | :---: | :---: | :---: |
| Fractional vegetation coverage (FVC) | $\begin{aligned} & \text { 2000, 2005, 2010, } \\ & \text { 2015, 2020; Raster, } \end{aligned}$ $1 \mathrm{~km}$ | NDVI data: Juntao (2022) | FVC $=($ NDVI - NDVIsoil $) /($ NDVIveg - NDVIsoil). Select the cumulative frequency of $5 \%$ and $95 \%$ of NDVI value as NDVIsoil and NDVIveg | 0.16 |
| Biological richness index (BRI) | 2000, 2005, 2010, <br> 2015, 2020; Raster, <br> 1 km | Land use and cover data: www.resdc.cn | $\mathrm{BRI}=\mathrm{Abio} \times(0.35 \times$ forest $+0.21 \times$ grassland $+0.28 \times$ water area $+0.11 \times$ cropland $+0.04 \times$ construction land $+0.01 \times$ unused land)/area. Abio is a normalized coefficient with a reference value of 511.264. Scale data into the 0-100 range using Min-Max Normalization | 0.19 |
| Net primary productivity (NPP) | 2000, 2005, 2010, 2015, 2020; Raster, 500 m | https://lpdaac.usgs. gov/products/ mod17a3hgfv061/ | Resampling to 1-km resolution. Scale data into the 0-100 range using Min-Max Normalization | 0.15 |
| Soil salinization index (SSI) | 2000, 2005, 2010, 2015, 2020; Raster, 500 m | MOD09A1 image: https://lpdaac. usgs.gov/products/ mod09a1v006/ | $\begin{aligned} & \mathrm{SI}=\sqrt{\rho_{1} \times \rho_{2}} \\ & \mathrm{SSI}=\sqrt{(\mathrm{NDVI}-1)^{2}+\mathrm{SI}^{2}} \end{aligned}$ <br> The $\rho 1$ and $\rho 2$ are the image reflectances of the blue and red light band of the MOD09A1 image (Yang et al., 2020). Resampling to $1-\mathrm{km}$ resolution. Scale data into the $0-100$ range using Min-Max Normalization | -0.14 |
| Soil erosion | 2000, 2005, 2010, 2015, 2018; Raster, 1 km | https://www.resde.cn/ | Scale data into the 0-100 range using Min-Max Normalization | -0.13 |
| Water conservation | 2000, 2005, 2010, 2015, 2020; Raster, 1 km | Wang (2022) | Scale data into the 0-100 range using Min-Max Normalization | 0.12 |
| $\mathrm{PM}_{2.5}$ concentration | 2000, 2005, 2010, 2015, 2020; Raster, 1 km | https://sites.wustl.edu/ acag/datasets/surface-pm2-5/#V5.GL. 02 | Scale data into the 0-100 range using Min-Max Normalization, and the maximum value is the average of the top 50 values in the 3 years | -0.11 |
未来 QTP 的发展应提出系统性的解决方案,以保持生态保护与社会经济发展的平衡(Li et al., 2022)。以下方法可能有助于实现双赢局面。
首先,有必要加快第三极国家公园群的建设,包括祁连山、普达措、青海湖和珠穆朗玛峰,并划定优先保护生物多样性、水源保护和生态系统碳汇的区域(赵等,2020)。帮助生活在这些核心保护区的人们逐步迁出。其次,我们需要完善生态补偿机制,并增加对青藏高原生态系统保护的财政支持。建议按照流域的上游、中游和下游划分生态产品的供给和使用,建立跨区域生态补偿机制,探索生态资源价值转化的途径。第三,采用区域轮牧、草地种植和冬季半舍饲养等方式,优化放牧管理,建立可持续的草地利用模式。第四,在道路、采矿、风能和光伏发电等设施建设和运营中,采用基于自然的解决方案,尽量减少对当地生态系统的影响,任何基础设施项目在其整个生命周期内必须是环境可持续的(李等)。2016 年)。第五,青藏高原的人均家庭能源消费大约是全国平均水平的两到三倍,而农业和畜牧业产生了大量的碳排放(Ping et al., 2011)。小型光伏发电站和水电站在这里被证明更具成本效益。我们将鼓励农业地区从传统能源转向清洁和可再生能源,并优化肥料使用和技术管理(Yu et al., 2022;Zhuang et al., 2021)。第六,高原城市的规模应合理规划,以保持在资源环境承载能力的阈值内,并应发展生态旅游,以减轻夏季游客对环境的损害(Hu et al., 2021)。
4.3. 研究不确定性与前景
EQI 的指数权重基于专家经验的评分,这在某种程度上是主观的。任何加权结构都可能受到一些人的批评,因为分配权重是一个依赖于价值的过程(Carpenter et al., 2009)。因此,我们随机将 EQI 某些指标的权重调整了 25%25 \% 上下,并发现结果的时空模式没有显著变化。在较长的时间尺度上,生态环境演变受到人类活动和气候变化的共同影响,而 QTP 上的 EQI 变化不能完全归因于人类活动(Zhang, Wang, et al., 2021)。然而,本研究的分析周期仅跨越了二十年,研究表明,在短期内,人类压力远远超过了气候对 QTP 的影响(Wei et al., 2022)。我们的研究目标不是分析生态环境变化的归因,而是探索人类与自然矛盾的风险区域以及人类与自然的相互作用。此外,我们忽略了时间选择的随机误差。我们仅从 2000 年到 2020 年每五年选择了一年,这可能会影响结果的稳健性。
我们没有考虑人类与自然之间的遥联。例如,青藏 Plateau 在夏季接待了许多游客。2019 年,青藏 Plateau 接待了来自全球的 1.2643 亿游客。他们对第三极施加的压力不可忽视。研究还表明,全球贸易对青藏 Plateau 的冰川产生了远程影响(Yi et al., 2019)。展望未来,大数据可能以更高的时间分辨率量化人类活动(Yi et al., 2020)。我们将进一步探讨人类与自然的互动机制以及人类活动的短期和长期影响。
Allison, A. E. F., Dickson, M. E., Fisher, K. T., & Thrush, S. F. (2021). 通过跨学科的人类-环境建模传达环境变化的驱动因素。地球的未来, 9(9), e2020EF001918. https://doi.org/10.1029/2020EF001918
Anderson, C. B., Seixas, C. S., Barbosa, O., Fennessy, M. S., Díaz-José, J., & Herrera-F, B. (2019). 确定自然对实现可持续发展目标的贡献。可持续性科学, 14(2), 543-547. https://doi.org/10.1007/s11625-018-0643-5
Bennett, D., & McGinnis, D. (2008). 结合与复杂:美国大黄石生态系统中的人类与环境互动。地理论坛, 39(2), 833-845. https://doi.org/10.1016/j.geoforum.2007.05.009
Boori, M. S., Choudhary, K., Paringer, R., & Kupriyanov, A. (2021). 基于压力-状态-响应框架的生态环境质量评估:遥感与地理信息系统的应用。遥感应用:社会与环境, 23, 100530. https://doi.org/10.1016/j.rsase.2021.100530
Bourdeau, P. (2004). 人与自然的关系与环境伦理. 环境放射性学杂志, 72(1), 9-15. https://doi. org/10.1016/s0265-931x(03)00180-2
Bu, Q., Dong, G., He, C., & Li, D. (2013). 青藏铁路格尔木至拉萨运营后的生态环境影响后评估. 铁路职业安全健康与环境保护, 3(3), 111-115.
蔡, B., 邵, Z., 方, S., 黄, X., 胡克, M. E., 唐, Y., 等. (2021). 社会经济活动与生态环境之间的细尺度时空耦合协调模型:以中国北京为例. 生态指标, 131, 108165. https://doi.org/10.1016/j.ecolind.2021.108165
Cardinale, B. J., Duffy, J. E., Gonzalez, A., Hooper, D. U., Perrings, C., Venail, P., 等. (2012). 生物多样性丧失及其对人类的影响. 自然, 486(7401), 59-67. https://doi.org/10.1038/nature11148
Carpenter, S. R., Mooney, H. A., Agard, J., Capistrano, D., DeFries, R. S., Díaz, S., et al. (2009). 生态系统服务管理的科学:超越千年生态系统评估。《美国国家科学院院刊》,106(5),1305-1312. https://doi.org/10.1073/pnas. 0808772106
陈, D., 陆, X., 胡, W., 张, C., & 林, Y. (2021). 城市扩张如何影响生态环境质量:基于空间杜宾模型的中国实证研究. 生态指标, 131, 108113. https://doi.org/10.1016/j.ecolind.2021.108113
陈, F. H., 韦尔克, F., 沈, C. C., 贝利, S. E., 伯格曼, I., 戴维斯, S., 等. (2019). 来自青藏高原的晚中更新世丹尼索瓦下颌骨. 自然, 569(7756), 409-412. https://doi.org/10.1038/s41586-019-1139-x
陈, G., 张, R., 郭, X., 吴, W., 郭, Q., 张, Y., & 闫, B. (2021). 西藏市政污泥利用过程的比较评估以实现可持续管理. 《环境科学总论》, 765, 142676. https://doi.org/10.1016/j.scitotenv.2020.142676