Abstract
Intensive agriculture with high reliance on pesticides and fertilizers constitutes a major strategy for ‘feeding the world’. However, such conventional intensification is linked to diminishing returns and can result in ‘intensification traps’—production declines triggered by the negative feedback of biodiversity loss at high input levels. Here we developed a novel framework that accounts for biodiversity feedback on crop yields to evaluate the risk and magnitude of intensification traps. Simulations grounded in systematic literature reviews showed that intensification traps emerge in most landscape types, but to a lesser extent in major cereal production systems. Furthermore, small reductions in maximal production (5–10%) could be frequently transmitted into substantial biodiversity gains, resulting in small-loss large-gain trade-offs prevailing across landscape types. However, sensitivity analyses revealed a strong context dependence of trap emergence, inducing substantial uncertainty in the identification of optimal management at the field scale. Hence, we recommend the development of case-specific safety margins for intensification preventing double losses in biodiversity and food security associated with intensification traps.
Similar content being viewed by others
Main
Rapidly rising global food demand creates a fundamental challenge for agricultural production to meet future needs and ensure food security1,2,3. Past increases in food production have primarily been achieved through cropland expansion and ‘conventional intensification’ (terms in single quotation marks are defined in Table 1)4,5,6. However, these increases came at the cost of substantial reductions in local biodiversity and associated ecosystem functions7,8,9, which can result in a strong negative feedback on yields (that is, productivity) and total agricultural production10,11,12.
The negative feedback of biodiversity on yields is especially of importance at high levels of ‘management intensities’. First, returns per effort decrease at high management intensities because of the saturating response of crop yields to conventional intensification13. Likewise, cropland expansion results in diminishing returns as it frequently occurs in marginal areas with lower yield potential owing to unavailability or protection of more suitable sites5,14. By contrast, negative impacts of intensification on biodiversity may even increase at high management intensities15, potentially causing the crossing of tipping points that trigger sudden community breakdowns16. The associated loss in crucial services, such as pollination or natural pest suppression, may outweigh direct benefits of higher management intensities and lead to hump-shaped production responses13,17. Under such circumstances, high management intensities lead to lose–lose situations instead of the frequently anticipated ‘trade-offs between production and biodiversity’18.
We refer to such lose–lose situations as ‘intensification traps’ as they are commonly linked to substantial negative legacy effects of past land use19. Over-intensification can cause, for example, soil biodiversity and fertility losses that require long-term restoration20,21 and thereby create barriers that prevent farmers from exiting trap situations. Likewise, above-ground biodiversity shows lagged responses to regenerative agroecological practices and the full recovery of ecosystem functionality requires extensive time periods22,23. Furthermore, if yield declines due to loss in biodiversity are misinterpreted as conventional yield gaps (that is, yields being limited by lack of inputs), they can lead to additional intensification and the self-reinforcement of traps. Hence, the avoidance of intensification traps, with their associated losses in biodiversity and food production, and their intransigence to reversal, needs to be a central goal of agricultural management24.
Despite their importance, the mechanisms driving intensification traps are conceptually not well resolved13,25. One associated challenge is that the crop, soil and biotic characteristics of landscapes are highly variable. This variability is transmitted to ‘biodiversity–production relationships’ that determine the occurrence of intensification traps. Hence, a crucial step to detect and prevent intensification traps is a clear mechanistic description of how crop, soil and biotic landscape characteristics shape the relationships underpinning the ‘management–biodiversity–production nexus’13.
Our aim in this study was therefore to identify the biophysical mechanisms driving intensification traps and evaluate how the emergence of traps varies with changes of crop, soil and biotic characteristics in agricultural landscapes. Our assessments are based on a novel analytical framework that integrates biodiversity as both predictor and response variable into agricultural planning. A core element of this framework is the conceptualization of five key relationships that together determine the covariation of biodiversity and crop production at the landscape level.
Five key relationships driving biodiversity and yield
Both crop production and biodiversity, represented here by species richness, depend on the intensity and the spatial extent of agricultural land use1,10. These dependencies can be described by five key non-linear relationships influencing crop yield directly or indirectly by mediating biodiversity effects (Fig. 1). These relationships vary in their effect sizes – the change in the response variable across the range of the predictor – and the shape of response curves, that is, the degree of their non-linearity; together, they define biodiversity–production relationships in agricultural landscapes.
-
(A)
Dependency of average yield on the spatial extent of production: The massive expansion of human land use over the last 300 years across most geographic regions has resulted in a scarcity of productive land5. Thus, agricultural expansion primarily occurs now in areas with lower yield potential26,27 and results in a negative, non-linear impact on the average attainable yield in a landscape. The effect size and shape of this relationship thereby depend on the frequency distribution of the yield potential and hence on the heterogeneity of biophysical production conditions (Fig. 1 and Extended Data Fig. 1).
-
(B)
Response of yield to conventional intensification: The saturating response of crop yield to external inputs such as fertilizers or pesticides, which are here primarily considered as conventional intensification, is long established28. Both effect size and shape of this relationship depend on biophysical properties and reflect, for example, crop nutrient requirements, nutrient deficiencies in unfertilized soils and the sensitivity of crops to pests.
-
(C)
Dependency of yield on biodiversity: The positive impact of biodiversity on yield is rooted in associated ecosystem functions such as soil nutrient cycling, pollination or natural pest suppression7,11. Whereas the effect size of this relationship depends primarily on crop requirements (for example, pollinator dependency29), its shape is determined from biological characteristics such as species’ effect traits30 or the degree of functional redundancy. For example, if pest-suppressing predators are primarily generalists and show high redundancy, a concave relationship can be expected. By contrast, the requirement of many specialized predators for effective pest suppression can be expected to result in convex relationships31.
-
(D)
Biodiversity responses to conventional intensification: Natural communities are known to respond negatively to eutrophication and pesticide applications32,33. The effect size of conventional intensification on biodiversity depends on the sensitivity of natural communities as well as on the type of input (for example, pesticides and fertilizers), which is often related to crop type. The shape of this relationship depends on the sensitivity of natural communities34, as high proportions of sensitive species result in rapid and hence convex responses whereas high proportions of tolerant species trigger a concave shape (Fig. 1).
-
(E)
Biodiversity responses to changes in the extent of agricultural land use: Biodiversity responses to changes in land-use types emerge from the responses of individual species that reside in a landscape or may colonize it from the regional species pool32,35. A defining parameter for their presence is the minimum amount of suitable habitat that is required by a species for persistence36. We, therefore, categorized species based on their ability to colonize agricultural land as well as natural and semi-natural habitats (further summarized as natural habitat). In addition, each species exhibits a threshold for the proportion of suitable habitat that is required for its persistence (Fig. 1, bottom). Together, these traits characterize species’ responses and thereby also the maximal species richness that can be reached under different land uses.
Land management is characterized by the level of conventional intensification effort (IE) and the proportion of land used for agriculture (that is, WL). Crop yield, that is, production per area (Y), depends on land-management features directly (plots A and B) and indirectly (C–E) via biodiversity. The effect size (change in the response variable across the range of the predictor) and the shape of these five relationships will vary across landscapes with crop, soil and biotic characteristics and key relationship drivers (listed in the top right). In our framework, the response of the attainable biodiversity (BA) to changes in land use results from habitat requirements of species in the regional species pool (bottom left) and their required minimum habitat size (bottom centre). Realized biodiversity (BR) is estimated by subtracting the negative impact of conventional intensification from BA.
Hence, the effect sizes and shapes of these five key relationships reflect the crop, soil and biotic characteristics of a landscape and define the responses of biodiversity and crop production to increasing management intensities. We evaluated how changes in these relationships impact the emergence of intensification traps, using a biodiversity–production model and a set of systematic literature reviews. Literature reviews were implemented for each model constant with the aim of capturing their natural variability across agricultural landscapes. Reviews included a restricted meta-analysis (see Section A2 of Supplementary Information for details) that was complemented by a snowball search to balance among crop types and geographic regions, resulting per model constant in over ten datasets for parametrization. This allowed us to (1) evaluate the occurrence of intensification traps across artificial landscapes, which were created stochastically to reflect the natural variability of our five key relationships; (2) establish three archetypal case studies contextualizing our results; and (3) implement a systematic sensitivity analysis to explore the model’s parameter space and identify mechanistic drivers of intensification traps at the landscape scale.
Results and discussion
Biodiversity–production patterns in agricultural landscapes
Our analyses evaluated the risk of entering intensification traps as well as their associated maximal production loss (Fig. 2a). We found that the implementation of the highest management intensities resulted in intensification traps in 73% of artificial landscapes. This trap prevalence directly emerges from the parametrization of our five key relationships based on the variability of real-world data recorded in our literature reviews. Hence, artificial landscapes represent the range of potential crop, soil and biotic characteristics, without being proportionate to the current prevalence of global production systems. Both risk and maximal production loss associated with traps are strongly driven by the effect size of biodiversity on agricultural yields in a given landscape (Fig. 2c,d). Yet, both these relationships show a high degree of scattering. This variability results from the joint impact of multiple drivers and also leads to situations in which trade-offs rather than intensification traps prevail (Fig. 2d).
a, Onset of intensification traps, definition of their risk and their associated maximal production loss. Both risk and maximal production loss are scaled from 0 to 1, with 1 denoting the highest theoretically possible value. b, The distribution of the risk of intensification traps and associated maximal production losses across 10,000 stochastically generated artificial landscapes. Artificial landscapes reflect the variability of the five key relationships (Fig. 1) as recorded in a literature review for each of the model parameters. Parameter values have been range transformed. c,d, The risk of intensification traps (c) and associated maximum yield losses (d) in 10,000 artificial landscapes are dependent on the effect size of biodiversity on yields. Each point represents a landscape. Note that landscapes with high risk (>80% of possible land uses) were occurring in less than 1% of all cases and were driven by situations when biodiversity peaked at intermediate to high contribution of working lands. e, The biodiversity gains attained by decreasing maximum production by 5%, 10% and 20%. Values to the right of the dashed lines indicate greater biodiversity gains than production losses (small-loss large-gain situation).
Three archetypal landscapes were chosen to provide contrasting examples of production systems that are of importance for global food security and largely differ in their five key relationships underlying intensification responses. These archetypal case studies also showed a large variation in the occurrence of intensification traps (Fig. 3). In two of the three landscapes, the US wheat belt and the Southeast Asian rice scenarios, even the highest management intensities did not trigger intensification traps despite the presence of positive biodiversity effects. By contrast, in the sub-Saharan small-holder scenario, a system with higher crop diversity and pollinator dependence, crop production was substantially reduced by high management intensities (Fig. 3). These contrasting responses emerge directly from differences in the parametrization of the five key relationships defined in our framework (see Extended Data Fig. 1 and below for an explanation of mechanisms). This suggests that globally dominant cereal production systems are less sensitive to biodiversity loss and intensification traps (crop type was much more important for parametrization than region; see Section A2 of Supplementary Information). However, the higher sensitivity of systems with a greater risk of intensification traps does not imply that the crop, soil and biotic characteristics of these landscapes are less ‘favourable’. Instead, a higher sensitivity signifies that these systems require a careful integration of biodiversity into management schemes to avoid biodiversity-driven yield gaps.
Exemplary case studies were chosen as they represent production systems of large importance for global food security that vary in their reliance of yield on biodiversity. In the middle row, responses of biodiversity (left) and agricultural production (right) to changes in management intensity (that is, conventional intensification and the proportion of agricultural land-use) are presented. In the bottom left, maximal attainable biodiversity (dotted circle) is depicted compared with the biodiversity maintained under the land management that leads to the highest agricultural production. Opportunity–cost curves (bottom right), which show the highest attainable biodiversity for each production level, are represented by red lines whereas model scenarios are shown as points coloured based on their conventional intensification effort (high, grey; low, yellow). B, biodiversity; P, production; NH, natural habitat. All axis units are range transformed. Credits, top row: left, josealbafotos/Pixabay; centre, Quangpraha/Pixabay; right, TG23/iStock.
Our analytical framework also allowed us to establish an ‘opportunity–cost curve’ of biodiversity and production for each individual landscape (Fig. 3). These curves depict the maximum biodiversity that can be attained at a certain production level. We found a strong prevalence of non-linear shapes for opportunity–cost curves across our analyses. Hence, small reductions in production can be ‘traded off’ for large biodiversity gains, which is exemplified by the rice case study, in which reducing the maximal crop production by 5% results in a doubling of biodiversity. Such small-loss large-win trade-offs also occur in the other case studies (Fig. 3) and are predominant across our artificial landscapes. Hence, reductions by 5–10% of maximal total production result frequently in disproportionately larger biodiversity gains, even in landscapes where intensification traps do not occur (Fig. 2e).
Understanding the drivers of intensification traps
The goal of our sensitivity analysis was to identify landscape characteristics that increase the likelihood of trap emergence. Systematic changes in the parameters defining the five key relationships of our conceptual framework revealed that both effect sizes and relationship shapes had a large impact on the risk of intensification traps (Fig. 4). These results show that intensification traps emerge in situations in which indirect consequences of biodiversity loss on yields outweigh the direct production-stimulating effects of increasing management intensities. Hence, the impacts of changing effect sizes and relationship shapes can be explained by their moderation of (1) the yield penalties caused by biodiversity loss and (2) the direct production benefits of increasing management intensities.
Biodiversity remaining in an agricultural landscape when at least 90% of maximum attainable production is achieved. In five sensitivity analyses (relationships A–E), the model constants describing one of the five key relationships were systematically modified over a predefined range (Supplementary Table 2). The dotted area represents conditions in which intensification traps emerge at maximal management intensity (that is, at maximal conventional intensification and agricultural land expansion). Yield potential in A refers to the highest attainable yield in a given area and was standardized from 0 to 1 in each landscape. Shapes of positive relationships apply to panels B–C, and those of negative relationships, to panel D (top right).
Direct production benefits of increased management intensities are governed by relationships A (yield responses to expansions of agriculture) and B (yield responses to conventional intensification). Relationship A is shaped by the distribution of the yield potential within a landscape, which is defined by its average and standard deviation. Both decreasing the average and increasing the standard deviation of the yield potential lead to a larger area of land that can attain only low yields (Extended Data Fig. 2). Consequently, maintaining or restoring marginal areas as natural habitats is linked to relatively small direct negative yield effects37, which can more easily be compensated by biodiversity-mediated yield benefits. Notably, this conclusion is based on the assumption that the potential of an area to support crop yields and biodiversity are uncorrelated as their covariation can affect land-use trade-offs38. In the case of intensification–yield relationships, represented by relationship B, reductions in effect size and increases in non-linearity have similar consequences (Fig. 4). Both lower direct production benefits attained from increasing management intensities from high to very high levels and therefore increase the risk of intensification traps.
Production gains resulting from increased biodiversity are determined by relationship C. Naturally, a higher effect size of this relationship increases biodiversity-mediated yield benefits and hence also the risk of intensification traps under conventional intensification (Fig. 4). Furthermore, strongly convex curves require that high levels of biodiversity are maintained to support agricultural production effectively. Such high levels of biodiversity can be attained only by substantial reductions in management intensities and pay off only when the effect size of biodiversity on yields is high (Fig. 4). Conversely, strongly concave biodiversity–yield relationships allow yields to benefit even from relatively low levels of biodiversity. Therefore, strongly concave relationships increase the risk of intensification traps but at the same time lower the biodiversity levels that are required to overcome trap situations (Fig. 4). Field and experimental studies show that the shape of biodiversity–ecosystem function relationships can vary greatly from strongly concave39 to convex responses30. However, the associated far-reaching impacts on the management of natural habitats in agricultural landscapes are rarely considered.
Relationship D, the response of biodiversity to conventional intensification, varies in agricultural landscapes40, for example, with the sensitivity of natural communities to pesticides and eutrophication. We found that a shift from concave to convex negative relationships, representing a higher community resistance, generally lowered the risk of intensification traps (Fig. 4). Convex relationships require that conventional intensification is drastically reduced before biodiversity can recover. Such large reductions in inputs are linked to a strong reduction in direct benefits of intensification, and the likelihood that they result in a net production increase is rather low. By contrast, the effect size of this relationship increases the risk of intensification traps. This pattern emerges because high effect sizes raise the amount of biodiversity that is lost by intensification, leading to stronger negative feedback on yields (Extended Data Fig. 3).
Finally, we considered changes in the regional species pool and its impact on biodiversity–production trade-offs (Fig. 4). The characteristics of the regional species pool are pivotal for biodiversity–ecosystem function relationships and local biodiversity patterns35,41 but receive little consideration in landscape planning. We found that changes in habitat preferences in the regional species pool trigger a sudden switch in the land-use strategy that optimizes crop production (Fig. 4 and Extended Data Fig. 4). If many species rely on natural habitats, the integration of these patches into landscapes substantially increased biodiversity and its associated benefits for yields. However, if many species required only working lands as habitat, yield benefits associated with natural habitats decreased. Furthermore, direct effects of conventional intensification on species in fields are much stronger than spillover effects on species in adjacent habitats42, which is integrated into our land-use model (Methods). Consequently, if a large proportion of species live only in working lands, a lower level of intensification is required to maintain these species, making the reconciliation of biodiversity and crop production more challenging.
From models to practice
Our analytical framework allows evaluating mechanistic biophysical drivers of intensification traps, which represent a substantial challenge for global food production and security7,13,17. However, agricultural systems show a high level of inherent complexity and a general assessment of intensification traps requires a number of simplifications. For example, both conventional intensification and biodiversity are in their essence multidimensional9,12,43 but their realistic representation would result in much higher model complexity, hampering conceptual advancements (see Section A3 of Supplementary Information for a detailed discussion). Hence, underlying model assumptions need to be taken into consideration when extrapolating our results to farm, landscape and regional scales.
At the farm level, financial cost–benefit relationships are key determinants of individual decision-making processes44. In this context, a clear understanding of intensification traps is crucial to avoid lose–lose situations that reduce production and biodiversity while raising farmers’ spending on pesticides and fertilizers. However, a precise determination of trap onsets requires detailed process-based information, which realistically cannot be compiled for each individual farm. The resulting uncertainty farmers are facing makes it economically advisable to follow precautionary principles and maintain management intensities slightly below anticipated optima. This would require, due to non-linear opportunity–cost curves, only small decreases in maximum farm revenue but substantially reduce the risk of intensification traps and long-term losses in fertility45. Hence, such safety margins restricting intensification would help to prevent financially highly detrimental lose–lose scenarios and enhance local biodiversity as a positive side effect.
In addition, farmers, in reality, may choose among a large variety of different land management practices, which our model framework simplifies into two dimensions (that is, land expansion and conventional intensification). An important alternative approach is ecological intensification46, which is frequently associated with practices such as intercropping or planting of cover crops17,47. These approaches are typically linked to higher labour but lower input costs48. Therefore, a thorough comparison of conventional and ecological intensification, which is beyond the scope of this study, would need to incorporate labour and other costs as an additional management dimension to identify optimal solutions in cost–benefit analyses13.
At the landscape level, the spatial arrangement of natural habitats represented by landscape structure is an essential factor influencing biodiversity–production relationships36. Landscape structure directly regulates the impact of natural habitats on agricultural yields49 and affects, by defining connectivity, the amount of realized biodiversity in local habitat patches as well as the size of regional species pools.35,41,50. We considered in our framework only the proportion of different land uses in a landscape and their impact on habitat availability for species in the regional special pool. However, specific landscape structure as well as other drivers of the size and trait frequencies in regional species pools (for example, evolutionary history linked to past agricultural practices51) can strongly moderate biodiversity–yield relationships and should be considered in the management of agricultural landscapes.
At the regional level, policies are required to prevent the biodiversity-dependent yield decreases shown in this study. In this context, the bimodal pattern in the distribution of the risk of intensification traps (Fig. 2b) and optimal management intensities (Extended Data Figs. 5 and 6) require further attention. These bimodal patterns result from the occurrence of two local optima of total production in the management-opportunity space of many agricultural landscapes. One of these optima is driven by the positive effects of conventional intensification on yield while the other emerges from yield benefits linked to high levels of biodiversity. Management intensities located between these optima have lower yields because the reductions in management intensities are insufficient to facilitate a functionally meaningful biodiversity recovery. Hence, policies that enforce or result in weak ecological minimum requirements (for example, a 3% non-productive farm area as a requirement for agricultural subsidies in the Common Agricultural Policy 2023–202752) might paradoxically risk promoting low points between production optima in many landscapes.
Conclusions
The prevention of intensification traps, which are characterized by a double loss of biodiversity and agricultural production, will be a crucial task for the sustainable management of life on earth53. We evaluated here intensification traps triggered by biodiversity loss, but our analyses can easily be expanded to other drivers such as soil degradation or salination. Our results highlight that the risk of intensification traps is increased by (1) a larger effect size of biodiversity on yields, (2) a stronger reliance of beneficial species on natural habitats and (3) stronger and more immediate responses of natural communities to conventional intensification. Furthermore, the risk of traps is decreased by stronger and more linear impacts of intensification on yields and with higher averages and less variable distributions of yield potentials in agricultural landscapes. Due to this complexity, it is difficult to quantify optimal management intensities at the farm level, and advisable to follow precautionary principles to avoid lose–lose scenarios. Furthermore, we found that across the vast majority of agricultural landscapes, small reductions in agricultural production can be translated into disproportionally larger biodiversity gains. These small-loss large-gain scenarios offer attractive opportunities to increase biodiversity in agricultural landscapes and can, along with a careful consideration of conservation targets, help to reconcile seemingly conflicting land-use targets.
Methods
Model framework
In our framework designed for landscape-scale assessments54, land management comprises two key aspects, (2) the proportion of land used for agricultural production (that is, working land (WL)) and (2) the level of conventional agricultural intensification (IE). IE represents the external inputs associated with conventional intensification such as fertilizer and pesticide use. Conventional intensification contrasts with ecological intensification, which includes the establishment of semi-natural habitat patches as one out of a large array of agroecological practices46. Our framework accounted for the possibility of regulating the proportion of semi-natural habitat patches, defined as 1 − WL. However, we decided to exclude other agroecological practices as those are associated with substantially higher labour inputs48 and the integration of labour costs was beyond the scope of this study (see Sections A2 and A3 of Supplementary Information).
Total crop production
Total agricultural production (PT) is computed as
where Y denotes the yield, that is, production per area. The direct dependency of yields on land management (relationships A and B) and the total biodiversity in a landscape (BT, relationship C) as well as the indirect dependency of yields on the response of biodiversity to management (relationships D and E) is described by
where YMax is the maximal attainable total production, while \({f}_{{I}_\mathrm{E}}^{\;Y}\left(\cdot \right),{f}_{{B}_\mathrm{T}}^{\;Y}\left(\cdot \right),{f}_{\mathrm{WL}}^{\;Y}\left(\cdot \right)\) and \({f}_{{I}_\mathrm{E}\mathrm{WL}}^{\;B}\left(\cdot \right)\) are four functions with five predictor terms that represent the five key relationships defined in Fig. 1.
The three terms that define the impact of IE and BT are each described by two model constants, the effect size of the independent variable on the dependent variable and the shape of the relationship (convex, concave). Effect size is confined to values between 0 and 1 and indicates the proportional change in the response variable if the predictor increases from 0 to 1 (that is, its range). The relationship slope is scaled from −1 to 1, where 0 represents a linear, 1 a convex and −1 a concave slope (see Section A1 of Supplementary Information and Supplementary Fig. 1).
Spatial elements are integrated into the functions in equations (2) and (3) that involved WL as a predictor. The impact of WL on yield (Fig. 1a) is based on the assumption that areas with a higher yield potential are first used for crop production and derived from the mean and the variance of the yield potential within a landscape (see Section A1 of Supplementary Information). Moreover, the response of biodiversity to changes in WL is based on the regional species pool (Fig. 1, bottom), which is recognized as a key determinant of biodiversity–ecosystem function relationships35,41. Each species in the species pool was linked to habitat requirements, which determine individual species’ responses to changes in agricultural land use (for details, see Section A1 of Supplementary Information).
We computed for each analysed landscape the biodiversity and production attained under different land management options. A land management option represented a combination of WL and IE, and for each landscape, 10,201 land management options were analysed (all possible combinations with WL and IE ranging from 0 to 1 in steps of 0.01). For simplicity, we range transformed attained biodiversity and total production across all 10,201 simulated land management scenarios to scale outputs from 0 to 1 within each landscape.
Implemented analyses
Our assessments of biodiversity–production relationships were based on three distinct approaches including (1) a stochastic landscape generation procedure, (2) the evaluation of archetypal case studies and (3) a systematic sensitivity analysis.
The variability of biodiversity–production relationships across landscapes was assessed by generating 10,000 artificial landscapes in a bootstrapping procedure. The analysis was based on a set of literature reviews to identify means and standard deviations of all model constants parameterizing equations (2) and (3). We accumulated, for the parametrization of individual model constants, between 11 and 26 datasets (Supplementary Table 2) and then used means and standard deviations to create a normal distribution for each model constant. An artificial landscape was created by randomly drawing a value for each model constant from its respective distribution, and analysed by establishing biodiversity and production outputs for each of the 10,201 land-management scenarios.
Archetypal case studies included a US wheat-belt scenario, a Southeast Asian rice scenario and an African small-holder scenario, which were based on intercropping and more diversified crop cultivation. They were used for contextualization of model results, and their parametrization is described in detail in Section A2 of Supplementary Information.
The response of model outputs to changes in individual model constants was investigated in a systematic sensitivity analysis. For each model constant, 100 landscapes were simulated and the target model constant was gradually changed from an upper to a lower range while all other model constants were maintained at mean literature values. Extreme but realistic values were chosen for the ranges of model constants (Supplementary Table 2). Biodiversity–production responses to different land-management scenarios were established for each landscape to identify the most important mechanisms driving intensification traps. All analyses were developed and implemented in R version 4.1.0 (ref. 55).
Reporting summary
Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.
Data availability
No original data were used for this work.
Code availability
Annotated versions of the model scripts, input data and a help document providing general instructions for the implementation of the model are provided under https://github.com/alfredburian/Intensification-traps.
References
Díaz, S. et al. Pervasive human-driven decline of life on Earth points to the need for transformative change. Science 366, eaax3100 (2019).
Mason-D’Croz, D. et al. Gaps between fruit and vegetable production, demand, and recommended consumption at global and national levels: an integrated modelling study. Lancet Planet. Health 3, e318–e329 (2019).
Willett, W. et al. Food in the Anthropocene: the EAT–Lancet Commission on healthy diets from sustainable food systems. Lancet 393, 447–492 (2019).
Tilman, D., Balzer, C., Hill, J. & Befort, B. L. Global food demand and the sustainable intensification of agriculture. Proc. Natl Acad. Sci. USA 108, 20260–20264 (2011).
Ramankutty, N. et al. Trends in global agricultural land use: implications for environmental health and food security. Annu. Rev. Plant Biol. 69, 789–815 (2018).
Dornelles, A. Z. et al. Transformation archetypes in global food systems. Sustain. Sci. 17, 1827–1840 (2022).
Potts, S. G. et al. Safeguarding pollinators and their values to human well-being. Nature 540, 220–229 (2016).
Newbold, T. et al. Global effects of land use on local terrestrial biodiversity. Nature 520, 45–50 (2015).
Beckmann, M. et al. Conventional land-use intensification reduces species richness and increases production: a global meta-analysis. Glob. Change Biol. 25, 1941–1956 (2019).
Garibaldi, L. A. et al. Mutually beneficial pollinator diversity and crop yield outcomes in small and large farms. Science 351, 388–391 (2016).
Albrecht, M. et al. The effectiveness of flower strips and hedgerows on pest control, pollination services and crop yield: a quantitative synthesis. Ecol. Lett. 23, 1488–1498 (2020).
Dainese, M. et al. A global synthesis reveals biodiversity-mediated benefits for crop production. Sci. Adv. 5, eaax0121 (2019).
Seppelt, R., Arndt, C., Beckmann, M., Martin, E. A. & Hertel, T. W. Deciphering the biodiversity–production mutualism in the global food security debate. Trends Ecol. Evol. 35, 1011–1020 (2020).
Kuang, W. et al. Cropland redistribution to marginal lands undermines environmental sustainability. Natl Sci. Rev. 9, nwab091 (2022).
Storkey, J., Meyer, S., Still, K. S. & Leuschner, C. The impact of agricultural intensification and land-use change on the European arable flora. Proc. Biol. Sci. 279, 1421–1429 (2012).
Dakos, V. & Bascompte, J. Critical slowing down as early warning for the onset of collapse in mutualistic communities. Proc. Natl Acad. Sci. USA 111, 17546–17551 (2014).
Garibaldi, L. A. et al. Working landscapes need at least 20% native habitat. Cons. Lett. 14, e12773 (2021).
Phalan, B., Onial, M., Balmford, A. & Green, R. E. Reconciling food production and biodiversity conservation: land sharing and land sparing compared. Science 333, 1289–1291 (2011).
Le Provost, G. et al. Land-use history impacts functional diversity across multiple trophic groups. Proc. Natl Acad. Sci. USA 117, 1573–1579 (2020).
Wubs, E. J. et al. Single introductions of soil biota and plants generate long‐term legacies in soil and plant community assembly. Ecol. Lett. 22, 1145–1151 (2019).
Seabloom, E. W., Borer, E. T. & Tilman, D. Grassland ecosystem recovery after soil disturbance depends on nutrient supply rate. Ecol. Lett. 23, 1756–1765 (2020).
Isbell, F., Tilman, D., Reich, P. B. & Clark, A. T. Deficits of biodiversity and productivity linger a century after agricultural abandonment. Nat. Ecol. Evol. 3, 1533–1538 (2019).
Moreno-Mateos, D. et al. The long-term restoration of ecosystem complexity. Nat. Ecol. Evol. 4, 676–685 (2020).
Cumming, G. S. et al. Implications of agricultural transitions and urbanization for ecosystem services. Nature 515, 50–57 (2014).
Seppelt, R. et al. Harmonizing biodiversity conservation and productivity in the context of increasing demands on landscapes. Bioscience 66, 890–896 (2016).
Gomiero, T. Soil degradation, land scarcity and food security: reviewing a complex challenge. Sustainability 8, 281 (2016).
Lesorogol, C. K. Land privatization and pastoralist well-being in Kenya. Dev. Change 39, 309–331 (2008).
Liebscher, G. Untersuchungen über die Bestimmung des Düngerbedürfnisses der Ackerböden und Kulturpflanzen. J. für. Landwirtsch. 43, 49–125 (1895).
Klein, A.-M. et al. Importance of pollinators in changing landscapes for world crops. Proc. Biol. Sci. 274, 303–313 (2007).
Kremen, C. Managing ecosystem services: what do we need to know about their ecology? Ecol. Lett. 8, 468–479 (2005).
Jiang, L. & Morin, P. J. Predator diet breadth influences the relative importance of bottom-up and top-down control of prey biomass and diversity. Am. Nat. 165, 350–363 (2005).
Gerstner, K., Dormann, C. F., Stein, A., Manceur, A. M. & Seppelt, R. Effects of land use on plant diversity—a global meta‐analysis. J. Appl. Ecol. 51, 1690–1700 (2014).
Schulz, R., Bub, S., Petschick, L. L., Stehle, S. & Wolfram, J. Applied pesticide toxicity shifts toward plants and invertebrates, even in GM crops. Science 372, 81–84 (2021).
Butsic, V. & Kuemmerle, T. Using optimization methods to align food production and biodiversity conservation beyond land sharing and land sparing. Ecol. Appl. 25, 589–595 (2015).
Zobel, M. The relative of species pools in determining plant species richness: an alternative explanation of species coexistence? Trends Ecol. Evol. 12, 266–269 (1997).
Martin, E. A. et al. The interplay of landscape composition and configuration: new pathways to manage functional biodiversity and agroecosystem services across Europe. Ecol. Lett. 22, 1083–1094 (2019).
Pywell, R. F. et al. Wildlife-friendly farming increases crop yield: evidence for ecological intensification. Proc. Biol. Sci. 282, 20151740 (2015).
Balmford, B., Green, R. E., Onial, M., Phalan, B. & Balmford, A. How imperfect can land sparing be before land sharing is more favourable for wild species? J. Appl. Ecol. 56, 73–84 (2019).
Tilman, D. et al. The influence of functional diversity and composition on ecosystem processes. Science 277, 1300–1302 (1997).
Newbold, T., Oppenheimer, P., Etard, A. & Williams, J. J. Tropical and Mediterranean biodiversity is disproportionately sensitive to land-use and climate change. Nat. Ecol. Evol. 4, 1630–1638 (2020).
Hagan, J. G., Vanschoenwinkel, B. & Gamfeldt, L. We should not necessarily expect positive relationships between biodiversity and ecosystem functioning in observational field data. Ecol. Lett. 24, 2537–2548 (2021).
Ricci, B. et al. Local pesticide use intensity conditions landscape effects on biological pest control. Proc. Biol. Sci. 286, 20182898 (2019).
Burian, A. et al. Low-cost management interventions and their impact on multilevel trade-offs in agricultural grasslands. J. Appl. Ecol. 60, 2079–2090 (2023).
Zou, Y. et al. Do diverse landscapes provide for effective natural pest control in subtropical rice? J. Appl. Ecol. 57, 170–180 (2020).
Chapman, M. et al. Social–ecological feedbacks drive tipping points in farming system diversification. One Earth 5, 283–292 (2022).
Bommarco, R., Kleijn, D. & Potts, S. G. Ecological intensification: harnessing ecosystem services for food security. Trends Ecol. Evol. 28, 230–238 (2013).
Kremen, C. & Miles, A. Ecosystem services in biologically diversified versus conventional farming systems: benefits, externalities, and trade-offs. Ecol. Soc. 17, 40 (2012).
Sánchez, A. C., Kamau, H. N., Grazioli, F. & Jones, S. K. Financial profitability of diversified farming systems: a global meta-analysis. Ecol. Econ. 201, 107595 (2022).
Brosi, B. J., Armsworth, P. R. & Daily, G. C. Optimal design of agricultural landscapes for pollination services. Conserv. Lett. 1, 27–36 (2008).
Fahrig, L. Ecological responses to habitat fragmentation per se. Annu. Rev. Ecol. Evol. Syst. 48, 1–23 (2017).
Garibaldi, L. A. et al. Smaller agricultural fields, more edges, and natural habitats reduce herbicide-resistant weeds. Agric. Ecosyst. Environ. 342, 108260 (2023).
Regulation (EU) 2021/2115. Establishing Rules on Support for Strategic Plans to Be Drawn Up by Member States Under the Common Agricultural Policy (CAP Strategic Plans) and Financed by the European Agricultural Guarantee Fund (EAGF) and by the European Agricultural Fund for Rural Development (EAFRD) and Repealing Regulations (EU) No 1305/2013 and (EU) No 1307/2013 (ed Development AaR) (The European Parliament and the Council of the European Union, 2021).
Post-2020 Global Biodiversity Framework: Discussion Paper (Convention of Biological Diversity, 2019).
Fahrig L. in Issues and Perspectives in Landscape Ecology (eds Wiens, J. A. & Moss, M. R.) 3–10 (Cambridge Univ. Press, 2005).
R Core Team. R: A Language and Environment for Statistical Computing. http://www.R-project.org/ (R Foundation for Statistical Computing, 2021).
Acknowledgements
A workshop initiating this article was funded by the University of British Columbia (Grants for Catalyzing Research Clusters to C.K. and N.R. through the Biodiversity Research Centre).
Funding
Open access funding provided by Helmholtz-Zentrum für Umweltforschung GmbH - UFZ.
Author information
Authors and Affiliations
Contributions
A.B., C.K., L.A.G., Z.M., N.M., M. Beckmann and R.S. developed the conceptual foundation of the study; A.B., M. Bulling and T.K. established the model code; A.B. and J.S.-T.W. implemented the meta-analyses; and all authors contributed to the establishment of the paper.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature Ecology & Evolution thanks Melissa Chapman, Bea Maas and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Extended data
Extended Data Fig. 1 Visualisation of the five key relationships for the three archetypal case-studies.
B stands for biodiversity, WL for working land (that is proportion of landscape used for crop production), Y for yield, non-WL for semi-natural and natural habitat (that is non-working land).
Extended Data Fig. 2 Examples for the distribution of production potential in agricultural landscapes.
A–D display distributions and E–H the respective relationships between average potential yield of used fields and proportion of land used for agriculture (curve is derived from distribution). Parameter settings for mean yield and variance of yields (ranges from 0 to 1) were respectively 0.75 and 0.05 (high yield, low variability), 0.25 and 0.05 (low yield, low variability), 0.75 and 0.5 (high yield, high variability) and 0.25 and 0.5 (low yield, high variability).
Extended Data Fig. 3 The dependency of biodiversity and total production on the intensity of conventional intensification and the extent of working lands in two artificial landscapes.
The two artificial landscapes have been parametrised based on mean literature values for four of the five key relationships and only differ in the parametrisation for relationship D, defining biodiversity responses to conventional intensification. Both landscapes show a linear impact of intensification on biodiversity but landscape 1 (top) was characterised by an effect size of 0.2 whereas the effect size in landscape 2 (bottom) was 0.8. The consequence is that in landscape 1, intensification traps cannot emerge because the negative impact of yield on biodiversity is to small that negative biodiversity feedback-effects on yields overcome the positive effect of intensification on yields. A contrasting situation is found in landscape 2.
Extended Data Fig. 4 Model responses to changes in habitat requirements of species in the regional species pool.
Habitat requirements are defined as inhabiting (i) natural habitat, (ii) working lands or (iii) the ability to survive in both. Displayed as response variables are (a) biodiversity in the scenario with the highest total crop production, (b) the conventional intensification effort in the scenario with the highest production, (c) the proportion of land used as working land in the scenario with the highest production, (d) the production achieved in the scenario with the maximum management effort (that is conventional intensification and agricultural land use are both at their maximum), (e) the biodiversity maintained in the scenario with the maximum management effort and (f) the curvature (that is measure for exponential nature) of biodiversity-production trade-off curves.
Extended Data Fig. 5 The land management options that support maximum total food production in 10000 artificial landscapes.
(a) The relationship between conventional intensification and the amount of land used for agricultural production, which represent the two components of land management considered in our study. Plotted is the land management that leads to maximal food production. Each point represents an individual landscape. The two land management components were positively related. The result of a type 2 regression (R2 = 0.36, p < 0.001) is depicted as black line. (b) The distribution of how strongly optimal land management (highest production) deviates from the maximum management scenario (maximal conventional intensification and agricultural land-use). The difference has been calculated as Euclidian distance and is scaled from 0 to 1 with the latter representing the largest possible distance.
Extended Data Fig. 6 Variability of conditions supporting maximal agricultural production.
Displayed are the distributions of biodiversity values, the amount of land used for agricultural production and the conventional intensification effort that support maximum total food production in 10000 artificial landscapes. Biodiversity values are scaled to the maximum biodiversity that can be reached in a landscape.
Supplementary information
Supplementary Information
Sections A1–A5, Supplementary Figs. 1–4 and Supplementary Tables 1–4. A1: Detailed description of model structure, A2: Model parametrisation, A3: Simplifications in model construction, A4: Additional tables and figures, A5: References.
Supplementary Data 1
All values for model constants derived from literature datasets.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Burian, A., Kremen, C., Wu, J.ST. et al. Biodiversity–production feedback effects lead to intensification traps in agricultural landscapes. Nat Ecol Evol 8, 752–760 (2024). https://doi.org/10.1038/s41559-024-02349-0
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/s41559-024-02349-0