Drivers of global change
We quantified the exposure of the world’s mountain areas to three major global change drivers: climate change, land-use change, and human population density. We focused on two different climatic variables: average monthly temperature (C°) and monthly mean precipitation amount (g/m3), extracted from CHELSA V1.2 database 85. Temperature is the standard variable used in calculating the velocity of climate change 7,45,46 and we decided to also consider precipitation which is often included in climate-velocity analyses, as they regulate the distribution and productivity of plant communities 40,45. We obtained the yearly mean temperature by averaging the monthly maximum and minimum temperature values across each year. We estimated the velocity of change and the magnitude of change separately for mean temperature and precipitation, and then we added the absolute values of the two climatic variables to consider the combined effect of climate change 45. We decided to use the absolute values as if temperature and precipitation velocities have opposing directions that are equal in magnitude, they would cancel, even if the climatic conditions are not the same as the original 40.
We extracted yearly land-use data from the Land-Use Harmonization (LUH2) dataset 86. We decided to focus on anthropic land-use categories which are considered to have overall negative impacts on biodiversity 87: cropland, pastureland, and urban.
We extracted human population density data from the human population’s count scenarios of the dataset of Olén & Lehsten, 2022, at 5-year intervals. Since we did not have past human population data from Olén & Lehsten, we used a different dataset 89 to extrapolate them, after verifying the spatial correlation between the two data in the present (Pearson’s coefficient = 0.80). We added a value of 1 to the rasters of the human population count used for interpolation to avoid having NA or infinite values in the final raster.
Data had a native original resolution of 1km for population data, 5 km for climate, and 25 km for land-use. We re-aggregated human population data at a resolution of 5 km to reduce the difference with land-use data, and we analyzed the other datasets at their native resolution.
Metrics of global change
For each driver of global change, we estimated two different metrics, the velocity and the magnitude of change. We estimated each metric for a time slice of 30 years, both in the past (1975–2005) and in the future (2020–2050), under two different emission scenarios, the SSP-RCP 2-4.5 and the SSP-RCP 5-8.5 (Tebaldi et al., 2020). According to SSP 2-4.5, global warming is estimated to be between 2.1 and 3.5 degrees Celsius in 2100, with development trends remaining mostly consistent with historical patterns 6,91. The SSP 5-8.5 scenario projects the highest emissions, with global warming by 3.3°C to 5.7°C °C by 2100 6,92.
We estimated the “acceleration” of each driver as the difference between the past and future velocity of change, and the Δ-magnitude as the difference between the past and future magnitude of change.
The velocity of change is a vector that describes how fast and in which direction a point on a gridded map would need to move to remain static in the same environmental conditions ( e.g., to maintain an isoline of a climate variable in the climatic space, under climate change) 40. We measured the velocity of change following the gradient-based approach 7, and using the gVoCC package 93 implemented in R (version 4.2.2, R Core team, 2022). Here the velocity of change is defined as the ratio between a temporal trend (the rate of change of each variable through time, estimated as a regression slope) and the corresponding spatial gradient of that variable (i.e., vector sum of longitudinal and latitudinal pairwise differences at each focal cell using a 3 x 3-cell neighbourhood):
\(\frac{slope (unit/year)}{spatial gradient (unit/km)}\) = VoCC (km/year) (1)
To estimate climate change velocity, we focused on mean temperature (C°) and precipitation amount (g/m3). We considered 4 Global Circulation Models (GCMs) for each future scenario, derived from the Coupled Model Intercomparison Project 5 (CMIP5): ACCESS 1–3, CESM1-BGC, CMCC-CM, MIROC 5. We obtained the climate velocity for each GCM under each emission scenario, and we averaged velocity estimates between GCMs for the same variable within each RCP scenario to produce an ensemble estimate. We also measured the uncertainty of the estimates by calculating the standard error across values.
We estimated the future velocity of change of land-use according to the same emission scenarios considered for climate (SSP 2-RCP4.5 and the SSP5-RCP8.5). In this case, we estimated the velocity of change for the three different categories of land-use (i.e., cropland, pastureland and urban) separately and then we summed the absolute values to obtain the overall velocity of land-use change.
We estimated the velocity of change of the human population density, extracting human population’s count scenarios of Coupled Model Intercomparison Project 6 (CMIP6) from the dataset of 88. For this variable, we scaled all values in log10 to reduce the skewness of data. The choice to use two different CMIPs, respectively CMIP5 for climate and CMIP6 for population density, was forced by the need of data with the finest possible temporal resolution, to estimate the velocity of change. However, this has not affected our evaluation of multiple global change drivers, as we scaled the velocity and magnitude of each driver for our comparison.
To quantify the acceleration of global change across the world’s mountains, we calculated the difference between the future velocity of change and the velocity of change of the baseline epoch. This allowed us to compare the velocity of change between the past and the future.
We estimated the magnitude of change for the same variables used to estimate velocity (i.e., climate, land use, human population density), for the baseline epoch (1975–2005) and future epoch (2020–2050). Following Nogués-Bravo et al., 2007, we estimated the magnitude of change according to the following formula:
$$M = \frac{(Xt2-Xt1)}{Xt1}$$
2
Where xt1 and xt2 represent values for the variable x at the beginning and at the end of the period.
To estimate values of the climatic variables, we considered the average values of the 10 years before and the 10 years after the reference year. For instance, for the historical epoch we considered the average values from 1965 to 1985 to estimate 1975 (t1) and the average values from 1995 to 2005 to estimate 2005 (t2).
We scaled the values of each variable based on the mean and the standard deviation of the year 2000, to obtain comparable estimates. As we did for the velocity of change, we considered different Global Circulation Models under two emission scenarios (i.e., SSP-RCP 2-4.5 and SSP-RCP 5-8.5), and then we averaged the magnitude estimates between GCMs within an RCP scenario to produce an ensemble estimate. We quantified the difference in magnitude of change across the world’s mountains (Δ-magnitude), as the difference between future magnitude and the magnitude in the baseline epoch. Such Δ-magnitude for climate is estimated as the sum of the absolute values of the differences estimated separately for temperature and precipitation.
We did not perform statistical tests to determine the significance of the differences observed between different epochs and regions, as in each comparison we consider the entire statistical population of values (i.e., every pixel). We represented the entire distribution of values for each region (mountain and lowland reference) and each epoch to allow for a full comparison (Fig. S7-S12).