Introduction

Changes in climate and land use are major drivers of biodiversity loss in this century (Sala et al. 2000; Pereira et al. 2010; Tittensor et al. 2014; IPBES 2019). These drivers likely interact and potentially change their mutual effects on biodiversity (Opdam and Wascher 2004; Oliver and Morecroft 2014; Hof 2021). Interaction mechanisms, however, are rarely considered by studies that assess both climate and land-use change impacts on biodiversity (Alkemade et al. 2009; Visconti et al. 2015). Such studies commonly use scenario and model analysis that combine bioclimatic and land-use variables to evaluate species and habitat range changes (e.g. Jetz et al. 2007; Pompe et al. 2008; Asner et al. 2010; Beltran et al. 2014; Brown et al. 2015; Byrd et al. 2015; Garcia-Valdes et al. 2015; Jantz et al. 2015; Rondinini and Visconti 2015; Visconti et al. 2015; Newbold 2018). Other studies use models that spatially overlay the individual effects solely showing combined or additive effects (e.g. Broennimann et al. 2006; Thuiller et al. 2006a; Mantyka-Pringle et al. 2012).

Interaction effects occur when two independent variables together have a larger (or lesser) effect on a dependent variable than the effects expected from the independent variables acting individually. Climate and land-use change effects will likely interact in several ways with respect to biodiversity (Brook et al. 2008; Hof 2021). In principle the effect of climate change (e.g. per degree of global mean temperature increase) on biodiversity depends on land-use change (e.g. per hectare of habitat converted to agriculture), and similarly the land-use change impact is contingent on climate change (Oliver and Morecroft 2014). The expected interaction effects will likely include changes in species’ fundamental niches, the degree of species’ range shifts and adaptive population responses (Sales et al. 2020; Latimer and Zuckerberg 2021).

Oliver and Morecroft (2014) described possible interaction mechanisms between climate and land-use change effects. Interaction effects, for example, indicate that species in heavily modified landscapes may respond differently to climate change (i.e. they have higher sensitivity levels) than species in pristine landscapes. Modified landscapes include fragmented or converted areas that reduce available species’ habitats and increase distances between remaining habitat patches (Fahrig 2003). These interactions mechanisms, however, are poorly understood due to limited data availability and methodological constraints (Brook et al. 2008; de Chazal and Rounsevell 2009; Hof 2021).

Under climate change, species shift their geographical ranges to find suitable climate. For some species, the new climate envelope may overlap with the old one (Vos et al. 2008). This depends on the time span and intensity of climate change, but also multiple new land uses influence the species’ response to climate change (Hannah et al. 2007; Bellard et al. 2012). Such land-use change include fragmented or converted land that likely constitutes a land barrier for species dispersal. A thorough understanding of species’ dispersal under projected land-use and climate changes is missing (Wiens et al. 2009; Hof 2021). Regardless of species’ dispersal being widely acknowledged in projecting geographic distribution of species (Thomas 2000; Guisan and Thuiller 2005; Midgley et al. 2006; Thuiller et al. 2008) how modified habitats hinder dispersal in response to climate change and how do interaction mechanisms arise, is poorly studied.

The primary objective of our study is to explore interaction effects from mechanisms between climate and land-use change. For this purpose, we selected four interaction mechanisms that are described in Oliver and Morecroft (2014). These mechanisms are (i) interactions between demographic factors, (ii) evolutionary trade-offs and synergies which impose selection on populations, (iii) threshold effects of population size on extinction risks and (iv) threshold occupancy for meta-population persistence. We identified expected consequences for each mechanism and validated the interactions with evidence from a literature review of studies that used bioclimatic models and climate-change scenarios in modified landscapes. We focused on intensively managed landscapes (i.e. croplands). We estimated the proportion of cropland area from these studies to indicate land use and land-use change. We then performed a meta-analysis and used the Fraction of Remaining Species (FRS) to indicate biodiversity decrease within each study region. FRS is a simple but effective indicator that represents the local response of species to climate change (Nunez et al. 2019). FRS decreases if the climate is not suitable anymore for a species at a location. The interaction mechanisms occur if the FRS in large cropland areas significantly differs from the FRS in areas with small cropland extent. By validating these expected consequences, we achieved a deeper understanding of drivers interactions that can be used in future biodiversity assessments.

Methods

Interaction mechanisms

Mechanism 1 describes interactive effects that can affect a number of demographic factors (e.g. birth rates, death rates and immigration/emigration). For example, edge effects may cause habitat desiccation and increased species mortality. This is exacerbated under extreme weather conditions. Mechanism 2 describes the evolutionary trade-offs and synergies where both climate change and land-use change impose selection on populations for more tolerant genotypes (i.e. adaptive capacity). Mechanism 3 describes threshold effects of population size on extinction risks. This is particularly critical for small populations which probably suffer larger extinction risks due to genetic drift, inbreeding depression, inability to find mates and increased susceptibility to environmental uncertainty. Finally, Mechanism 4 describes effects on the threshold occupancy for meta-population persistence. This mechanism shows that when the proportion of patches falls below a threshold level, the extinction of the entire meta-population can rapidly follow. We identified expected consequences that can further elucidate on these interaction mechanisms (Table 1). These expected consequences indicate that species’ responses to climate change vary if the landscape changes (i.e. different levels of cropland area).

Table 1 Interaction mechanisms between climate change and land-use change effects and expected consequences on (i) species distributions and (ii) Fraction on Remaining Species (FRS)

All four mechanisms mimic ecological processes that result in declining populations if climate and land-use change effects interact upon species. Such processes likely occur in areas with remaining suitable climates. Thus, declining populations can be evidenced with meta-analysis from bioclimatic modelling studies and the FRS indicator. If such location also experiences land-use change, climate-change impacts are likely to increase and species are expected to be more sensitive on average, and thus disappear.

Literature review

We conducted a literature review to find evidence on the expected consequences of the interaction mechanisms. We queried the ISI Web-of-Science to identify bioclimatic model and scenario studies that (1) assessed the effects of global mean temperature increase on terrestrial species and (2) included cropland as one of the land-use types in their study area. We combined keywords into a search string to conduct the review: [(climat* SAME change*) OR (temperature SAME change*) OR (temperature SAME increase) OR (global SAME warming) OR (climat* SAME warming)] AND [(biodiversity) OR (diversity SAME species) OR (species SAME richness*) OR (species SAME distribution*) OR (species SAME abundance*) OR (species SAME occurrence*) OR (species SAME turnover) OR (species SAME loss*) OR (species SAME gain*) OR (species SAME composition) OR (species SAME assemblage*)] AND [(land use SAME change*) OR (fragment* SAME landscape) OR (fragment* SAME habitat*) OR (fragment* SAME habitat SAME connectivity)] AND [(dispersal SAME capacity) OR (species SAME climat* SAME sensitivit*)]. We defined inclusion criteria to determine the relevance of the studies. These criteria request that studies to: (i) present data on species composition (i.e. number of species before and after climate change); (ii) indicate the location of the study area to calculate spatial extent and land-use intensity; (iii) provide information on global temperature increase used in the original analysis, or indicates the type of climate-change scenario and general circulation model used; (iv) indicate the reference and projected years for climate change; and (v) specify the taxonomic groups assessed in the original analysis. We screened studies by title for their relevance to the purpose of our study. Those studies selected were then screened by abstract and, when providing limited information in their abstracts, were fully screened, both in content and supporting material. We consolidated a database from all data extracted from the selected studies to perform meta-analysis on species distributions in landscapes with varying cropland levels. These data included study location, global mean temperature increase (°C), land use/land cover, number of species in the original and the projected climate situations, dispersal capacities, taxonomic group, spatial extent and the spatial resolution. We also qualitative assessed studies on the impact of habitat fragmentation and climate change on species dispersal.

Projected land-use and climate changes

Land use was indicated by the percentage of cropland in the original study areas. We estimated cropland percentage by using the land-cover maps from the European Space Agency (ESA) Climate Change Initiative (CCI) (Defourny et al. 2014). These are consistent global maps at 300 m spatial resolution on an annual basis from 1992 to 2015. We used the year of the original study. We aggregated land-cover categories that indicate intensively managed areas. These categories are: cropland (rainfed), herbaceous cover, cropland (irrigated or post-flooding), mosaic cropland (>50%) and mosaic natural vegetation. We used ArcGIS 10.2.1 to estimate the spatial extent of the study sites and the cropland percentage. We also identified other important characteristics of the original study areas, such as study areas located in protected areas, to determine the applied level of change. Land-use change was classified as very low, low, moderate, high and very high. These classes respectively correspond to 0–5%, 0–10%, 10–15%, 15–20% and > 20% of cropland (i.e. cropland level) respectively. Climate change in these studies was indicated by their scenario’s global mean temperature increase.

Meta-analysis

We estimated the overall reduction in species abundance in landscapes with varying cropland levels with a mixed-effects meta-analysis. As an effect size, we calculated the FRS. FRS quantifies the remaining proportion of biodiversity compared to the original situation in the selected studies. This also indicates changes in species abundance of the local biodiversity. Following Nunez et al. (2019), FRS was calculated as the average of ratios between the number of species in the projected climate situation and the original number of species within each locality in the original study sites (Eq. 1):

$$FRS= \frac{1}{n} . \sum_{i=1}^{n}\frac{{S}_{di}}{{S}_{oi}}$$
(1)

where Sdi is the expected number of species in locality i after climate change and Soi is the number of species in locality i in the original situation. n is the total number of localities.

We used the package ‘Metafor’ in the R-3.2.2 software (Viechtbauer 2010) and the rma.mv() function assuming independence between the effect size and sampling variance. We built mixed-effects models with random-effects structures to examine variations in the relationship between FRS and both cropland level and global mean temperature increase. We included first-order interaction terms and second-order interaction terms to account for variation in the slope of the relationship between global mean temperature increases and cropland levels. We compared the models using the Bayesian Information Criterion.

Result

The systematic literature review yielded more than 100 studies that contained relevant titles and abstracts. Of these, 42 studies corresponded to the selected criteria for data extraction (see Electronic Supplementary Material S1: List of selected studies). Some geographical bias was found since most of the studies stemmed from Western Europe (19) or Africa (10), while a few studies from the Americas (7), Oceania (4) and China (2). These studies covered a wide range of temperature increase up to 5 °C. The study taxa included plants, vertebrates (e.g. birds and mammals) and insects (e.g. butterflies). The published years of the selected studies (i.e. 1994–2015) matched the land-cover ESA CCI data. We collated these studies, including 219 effect sizes (see Electronic Supplementary Material S2: Database). These effect sizes corresponded to very low (53), low (54), moderate (31), high (12), and very high (69) cropland levels.

For the meta-analysis, we retained the random-effect structure [1 | Study + 1 | Extent] with the lowest Bayesian Information Criterion (i.e. −197.67) (see Electronic Supplementary Material S3: Results of random effect models). The retained structure contains the unique ID (i.e. ‘Study’) for each selected study and the spatial extent (i.e. ‘Extent’) of the study site indicated by the number of pixels. The second order interaction term (i.e. squared global mean temperature increase · cropland level) did not improve the model fit and was therefore dismissed. The final model used was global mean temperature increase + cropland level + global mean temperature increase · cropland level.

The results from the mixed-effects model indicated that global mean temperature increase is a good, significant model predictor for biodiversity decline (i.e. estimate = −0.0628 with 95% confidence interval: −0.0881 to −0.0374; Table 2). However, the model results did not reveal significant interactions between cropland levels and global mean temperature increases, which indicates that the influence of different cropland levels on FRS did not change with global mean temperature increase (i.e. estimate very low = −0.0342 with 95% confidence interval: −0.2299 to 0.1615 and estimate very high = −0.0019 with 95% confidence interval: −0.0627 to 0.0589; Table 2 and Fig. 1).

Table 2 Results of the mixed-random effects model ‘global mean temperature increase + cropland level + global mean temperature increase · cropland level’ with FRS as response variable
Fig. 1
figure 1

Schematic illustration to indicate expected consequences of interactions on the Fraction of Remaining Species. a shows the reference situation of original species at a location; b shows the projected situation with low land-use change and climate change; and c shows the projected situation with high land-use change and climate change. Projected situations 1 and 2 show species loss (in orange) from original area with suitable climate, and remaining species (in blue) in the remaining area with suitable climate. This information is used to calculate FRS in d. Theoretically, FRS differs between projected situations 1 and 2 (i.e. 0.5 and 0.3) due to changed vulnerability to climate change of remained species after land-use change (mechanism 1), increased species tolerance for future land-use change (mechanism 2) and changes in extents of occurrence of species and their vulnerability to climate change (mechanisms 3 and 4). However, the estimates from the mixed-effects model indicate that FRS of areas with high cropland level (i.e. estimate = −0.0342) did not significantly differ from those areas with a low cropland level estimate = −0.0019)

A mixed-effect model using global mean temperature increase as the only predictor variable was significant (Table 3). This suggests that temperature increase strongly influenced FRS. Contrary, a mixed-effect model using cropland level as the only predictor variable was less significant, which indicates that land use did not influence or the influence was not significant for FRS.

Table 3 Results of the mixed-random effects model for individual predictors ‘global mean temperature increase’ and ‘cropland level’ with FRS as response variable

Our analysis of global mean temperature increase indicated a general decrease of FRS. Both estimates for temperature in Tables 2 and 3 resulted in FRS reduced by 14% (95% confidence interval: 8% to 18%) and 12% (95% confidence interval: 9% to 15%), respectively.

Discussion and conclusions

Based on the expected consequences for the studied interaction mechanisms, we found that the FRS of areas with a large cropland level (i.e. intensively managed) did not significantly differ from those areas with a low cropland level (i.e. close to pristine, natural areas). This finding excluded interaction effects between global mean temperature increases and landscapes that are modified by other land uses than cropping. We, however, did not find any study that assessed interactions between land-use and climate-change effects. This makes assessing direct interaction mechanisms in other modified landscapes unattainable.

The few studies for mechanism 1 that affected species’ dispersal in response to global mean temperature increases and habitat fragmentation, indicated that species’ dispersal ability is closely related to the landscape structure and fragmentation level (e.g. Wiens et al. 2009; McGuire et al. 2016) and thus only species with unlimited dispersal capacities can always establish remote populations in fragmented landscapes (Opdam and Wascher 2004). This indicates that increasing landscape fragmentation will likely hinder species’ distribution shifts in response to climate change. We did not find any study that quantified climate induced dispersal of species in fragmented landscapes. Most found studies focused on either single impacts of fragmentation (Thomas 2000; Fahrig 2003), climate change (e.g. Araujo et al. 2004; McClean et al. 2005; Thuiller et al. 2005, 2006a, b; Lawler et al. 2006; Malcolm et al. 2006; Alkemade et al. 2011; Walther and van Niekerk 2015), assessed single species responses (Bennie et al. 2013) or assessed past dispersal trends (Schwartz 1993).

Previous studies that assessed additive land-use and climate-change effects (Mantyka-Pringle et al. 2012; Newbold 2018), indicated that changes in biodiversity likely vary spatially. For example, tropical regions are projected to experience larger losses of biodiversity under climate change (Higgins 2007; Asner et al. 2010) than temperate regions (Garcia-Valdes et al. 2015). Temperate regions have already undergone some of the largest biodiversity losses from land-use changes over centuries and thus the climate effect is likely small (Newbold 2018). Our analysis, however, showed that intensively managed landscapes did not change species responses to climate change. This was particularly important to verify mechanism 2. These species responses are likely most sensitive in deserts, temperate forests and shrublands, regardless of land use (Nunez et al. 2019). This spatial variation was neglected by most studies that simultaneously assessed combined effects of these drivers in multiple land-cover types (e.g. Heubes et al. 2013; Jantz et al. 2015; Riordan et al. 2015). While most studies that assessed changes in climate and land use, found an increased combined effect (e.g. Jetz et al. 2007; Pompe et al. 2008), a few showed that the combined effect sometimes had a weak impact on biodiversity (Hoiss et al. 2013; Brown et al. 2015).

Our assessment approach of interaction effects, and particularly from mechanisms 3 and 4, included factors, which involve large uncertainties such as species’ dispersal capacities (Hellmann et al. 2016). Thus mechanisms could be operating simultaneously and potentially in opposite directions from each other, making it less likely that significant changes would be observed in the meta-analysis. Other inherent factors that affect interaction mechanisms include species adaptive capacity (Dawson et al. 2011) and the more gradual and long-term effects of climate change (Nunez et al. 2019) and changing atmospheric CO2 concentrations. Such factors also involve large uncertainties in projecting species distributions in response to climate change and need to be included when assessing interaction mechanisms. Furthermore, future assessments should also consider the physiological characteristics of species and the time lags resulting from the expected changes in both land use and climate (Hof 2021). Although future projections of land-use and climate change are challenging to create, our study serves as a starting point to introduce more elements to assess such changes and their interactions mechanisms.

The lack of integration of land-use (e.g. intensively managed areas) and climate-change effects in interaction mechanisms implies that the projected biodiversity responses probably inadequately inform biodiversity managers and policy makers on possibilities to develop appropriate biodiversity-conservation measures (de Chazal and Rounsevell 2009; Mantyka-Pringle et al. 2012, 2015; Heubes et al. 2013; Brown et al. 2015; Segan et al. 2016). This knowledge gap, however, can certainly be explored using the expected consequences to better validate interaction effects. While our study concludes that these interaction mechanisms have low significance in intensively managed areas, we do not rule out drivers interactions in other land-use types. Stronger indicators can further elucidate these mechanisms. Our study strongly supports the notion that climate and land-use changes dramatically alter biodiversity and thus improved understanding of how climate and land use-change effects interact could improve biodiversity-loss projections. Integrating these interactions effects into biodiversity-impact assessments will likely become important to develop more comprehensive and reliable strategies for biodiversity conservation and restoration.