# An increase in food production in Europe could dramatically affect farmland biodiversity

Philippe Jeanneret, Gisela Lüscher, Manuel K Schneider, Philippe Pointereau, Michaela Arndorfer, Debra Bailey, Katalin Balázs, András Báldi, Jean-Philippe Choisis, Peter Dennis, et al.

<figure><img src="image_1.png" /></figure>

## ► To cite this version:

Philippe Jeanneret, Gisela Lüscher, Manuel K Schneider, Philippe Pointereau, Michaela Arndorfer, et al. An increase in food production in Europe could dramatically affect farmland biodiversity. Communications Earth & Environment, 2021, 2 (1), pp.183. 10.1038/s43247-021-00256-x. hal-03356520

HAL Id: hal-03356520

https://hal.inrae.fr/hal-03356520

Submitted on 28 Sep 2021

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>
<figure><img src="image_3.png" /></figure>

# An increase in food production in Europe could dramatically affect farmland biodiversity

Philippe Jeanneret<fnref n="1" />✉, Gisela Lüscher<fnref n="1" />,², Manuel K. Schneider<fnref n="1" />, Philippe Pointereau³, Michaela Arndorfer⁴,
Debra Bailey<fnref n="1" />, Katalin Balázs⁵, András Báldi⁶, Jean-Philippe Choisis<fnref n="1" />,⁷, Peter Dennis⁸, Mario Diaz<fnref n="1" />,⁹,
Sebastian Eiter<fnref n="10" />, Zoltán Elek⁶, Wendy Fjellstad<fnref n="10" />, Thomas Frank⁴, Jürgen K. Friedel⁴,
Ilse R. Geijzendorffer<fnref n="1" />,¹¹,¹², Pippa Gillingham⁸, Tiziano Gomiero¹³,¹⁴, Gergely Jerkovich⁸, Rob H. G. Jongman¹¹,
Max Kainz¹⁵, Anikó Kovács-Hostyánszki⁶, Gerardo Moreno<fnref n="1" />,¹⁶, Juri Nascimbene¹⁷, Marie-Louise Oschatz⁴,
Maurizio G. Paoletti¹³, Jean-Pierre Sarthou<fnref n="18" />,¹⁹, Norman Siebrecht¹⁵, Daniele Sommaggio¹³,
Sebastian Wolfrum¹⁵ & Felix Herzog<fnref n="1" />

Conversion of semi-natural habitats, such as field margins, fallows, hedgerows, grassland, woodlots and forests, to agricultural land could increase agricultural production and help meet rising global food demand. Yet, the extent to which such habitat loss would impact biodiversity and wild species is unknown. Here we survey species richness for four taxa (vascular plants, earthworms, spiders, wild bees) and agricultural yield across a range of arable, grassland, mixed, horticulture, permanent crop, for organic and non-organic agricultural land on 169 farms across 10 European regions. We find that semi-natural habitats currently constitute 23% of land area with 49% of species unique to these habitats. We estimate that conversion of semi-natural land that achieves a 10% increase in agricultural production will have the greatest impact on biodiversity in arable systems and the least impact in grassland systems, with organic practices having better species retention than non-organic practices. Our findings will help inform sustainable agricultural development.

<footnote n="1">¹ Agroscope, Zurich, Switzerland. ² Institute of Evolutionary Biology & Environmental Sciences, University of Zurich, Zurich, Switzerland. ³ SOLAGRO, Initiatives and Innovations for Energy, Agriculture and Environment, Toulouse, France. ⁴ University of Natural Resources & Life Sciences, Vienna, Austria. ⁵ Institute of Environmental & Landscape Management, Szent Istvan University, Gödöllő, Hungary. ⁶ MTA Centre for Ecological Research, Vácrátót, Hungary. ⁷ INRA, UMR 1201 Dynafor, Castanet-Tolosan, France. ⁸ Aberystwyth University, Aberystwyth, UK. ⁹ Museo Nacional de Ciencias Naturales (CSIC), Madrid, Spain.</footnote>

<footnote n="10">¹⁰ Norwegian Forest and Landscape Institute, Ås, Norway. ¹¹ Alterra, Wageningen UR, Wageningen, The Netherlands. ¹² Institut Méditerranéen de Biodiversité et d'Ecologie marine et continentale, Aix-en-Provence, France. ¹³ Department of Biology, Padova University, Padova, Italy. ¹⁴ Universitat Autonoma de Barcelona, Bellaterra, Spain. ¹⁵ Centre of Life and Food Science, Technical University of Munich, Freising, Germany. ¹⁶ Forestry School, University of Extremadura, Plasencia, Spain. ¹⁷ Department of Biological, Geological and Environmental Sciences (BiGeA), University of Bologna, Bologna, Italy.</footnote>

<footnote n="18">¹⁸ Toulouse University, ENSAT, Toulouse, France. ¹⁹ INRA, UMR 1248 AGIR, Castanet-Tolosan, France. ✉email: philippe.jeanneret@agroscope.admin.ch</footnote>
Global food needs are expected to increase drastically by 2050 if current trends continue, especially in population growth and increased consumption of meat and dairy products¹⁻³. Worldwide crop yield increases are unlikely to meet the projected demand for food⁴⁻⁶, and the gap⁷ may be further exacerbated by climate change⁸,⁹. A possible strategy to meet the production shortfall will be further expansion of productive agricultural land⁵,¹⁰,¹¹ which—combined with the competition for land to produce energy and develop infrastructures—will further jeopardize natural and semi-natural areas worldwide and in Europe¹²,¹³.

In Europe, farmland covers 47% of the EU-28 area¹⁴ and has historically supported high levels of biodiversity¹⁵. While large wilderness areas are prioritized to conserve biodiversity globally, conservation effort in the comparatively small scale mosaic-type European landscapes requires consideration of farmland¹⁶,¹⁷ and associated semi-natural habitats such as field margins, fallows, hedgerows, low-input grassland, woodlots and traditional agro-forestry. Whilst simplification of the agricultural landscape has halted in recent decades¹⁸, partly due to agri-environmental policies, the question now arises, what production gain and biodiversity loss would be, if the pressure on the land would increase again and (part of) the remaining semi-natural habitats were converted to agricultural fields. To tackle the question we first investigated the set of species that are unique to European semi-natural habitats and that would disappear in case of conversion to production fields, and related that loss to the production gained by the conversion¹⁹, also in case of conversion being to organically or conventionally managed fields. Organic farming is often proposed as a solution to feed a growing population while simultaneously minimizing the global environmental impact of agriculture⁵,²⁰⁻²² and supporting biodiversity²³⁻²⁵. Since organic systems in Europe usually show lower yields (20-40%; but see²⁶,²⁷) they require more land to produce the same amount of food. This raises the second question: what would a comparative loss of species in organic and non-organic systems, and on a difference farming systems (grassland vs arable vs horticulture etc.) be in case overall agricultural production should increase by 10%?

We evaluated the current productivity on farmland and the species richness of key taxonomic and functional taxa for agroecosystems on both production areas and semi-natural habitats in ten agricultural regions across Europe, from boreal to Mediterranean (Fig. 1a, Supplementary Table 1). Regions comprised various agricultural land uses from arable to grassland, mixed, horticulture and permanent crops. In each region, we

randomly selected 12–20 farms from sets of 30–40 farms that had accepted to participate, approximately half certified organic. We surveyed four taxa (vascular plants, earthworms, spiders, wild bees) that represent key ecosystem functions for agriculture (contributing to soil fertility, biological pest control and pollination), different habitat compartments (soil and aboveground structures) and mobility in 910 semi-natural habitats and 492 production fields, using hierarchical preferential sampling and standardized protocols (Methods, Supplementary Table 2). Habitats were categorized as either semi-natural habitats or production fields according to criteria based on Raunkiær plant life forms²⁸ and management evidence²⁹ (Methods). Farmers provided data on yields. In each case study region, the recorded species were grouped as (i) unique to semi-natural habitats, (ii) unique to production fields or (iii) shared by both habitat categories. The use of unique species enabled the evaluation of biodiversity loss if habitats should be converted into other land use categories.

## Results

**High proportion of species unique to semi-natural habitats.**

Summed over the four taxa, 49% (±3%) of all species on average across regions were unique to semi-natural habitats, which was significantly more (n = 40, $P_{χ^2}^1 < 0.001$) than the 26% (±3%) unique to production fields (Fig. 1b) compared for the same sample coverage per region³⁰ (Methods, per taxa separately in Supplementary Table 3). The remaining 25% (±4%) of the species were shared by the two habitat categories. The contribution of semi-natural habitats to farmland biodiversity is strikingly high, given the limited area of 23% (min 3% – max 71%) of farmland they covered per region. Remarkably, on arable, horticultural and wine producing farms (AT, FR, NL, DE, IT), the 3–10% semi-natural habitats hosted 40–60% of the unique species. Considered separately, plant, spider and wild bee species showed similar patterns (Supplementary Fig. 1). Across the four taxa and summed over the regions, the differences between organic and non-organic systems for unique species in both semi-natural habitats and production fields, and shared species, were mainly non-significant (Supplementary Table 4), which is not surprising since the two farming systems of the study were overall very similar and not significantly different with respect to main land characteristics (Supplementary Table 5). However, the less intensive management of the organic farms led to a significantly lower yield than in the non-organic ones (Supplementary Table 6).

<figure><img src="image_4.png" /><figcaption>**Fig. 1** Importance of semi-natural habitats for biodiversity in ten European farmland regions. a Location of the study regions with dominant type of agricultural land use (brown: arable crops, mixed systems and horticulture, bright green: grassland, violet: permanent crops). b Proportions of summed vascular plant, earthworm, spider and bee species unique to semi-natural habitats (light green bars), shared by the semi-natural habitats and production fields (medium green bars) and unique to production fields (dark green bars) per region (per taxa separately in Supplementary Fig. 1). c Percentage of farm area occupied by semi-natural habitats (light green bars) and production fields (dark green bars) per region, and average yield (stars) over fields of sampled farms in the ten study regions.</figcaption></figure>
**Moderate production gain and dramatic species loss by conversion of semi-natural habitats.** In the context of land take for food production, the conversion of semi-natural habitats to production fields would increase the production and the species richness unique to production fields due to the area gained, but species unique to semi-natural habitats would be lost. To quantify the trade-off between food production and biodiversity conservation, we modelled the impact of a progressing conversion of semi-natural habitats to production fields from sample-size based extrapolations of unique species occurring there³¹ (Methods). Potential production gain was calculated by multiplying the average yield over the various crops by the area of converted land. Yields were expressed in mega joules (MJ) per hectare (ha) of harvestable energy. They varied considerably within and between case study regions, as a function of the geographic location and the farm type (Fig. 1c). Average yield per region ranged from 31 461 MJ ha⁻¹ (IT) to 141 577 MJ ha⁻¹ (CH). Productivity (evaluated as MJ per hectare, Fig. 1c) was not related to the share of semi-natural habitats (Spearman’s rho = 0.095, ns) suggesting that farms with large semi-natural habitat area could also perform well regarding production. To investigate the pattern of species loss and production gain, we calculated and plotted the accumulation curves of species unique to semi-natural habitats backward, and added the corresponding gained species unique to production fields and species shared (Methods, Supplementary Fig. 2). The pattern of species richness decline (all taxa combined) with accumulative conversion was similar across regions (Fig. 2), and organic and non-organic systems (Supplementary Fig. 3). The concave downward curves (slope increased with larger area of semi-natural habitats converted) indicate a low impact of converting few semi-natural habitats; this suggests high similarity of species compositions among them. Even, species richness increased slightly in FR and DE first due to the overcompensation of species unique to production fields compared to species unique to semi-natural habitats. This was particularly true for organic systems in FR, DE, NO and ES (Supplementary Fig. 3). This indicates the buffering role of organic farming in case of

semi-natural habitat loss. However, when more than half of the semi-natural habitats were converted, species loss increased dramatically on organic farms as well.

In all study regions, the impact of converting 50% of the semi-natural area on species richness estimates was low and partly compensated by an increase of species numbers on the larger available production area (Table 1 and Supplementary Table 7). However, converting 90% of the semi-natural area (assuming that 10% are unconvertible due to geomorphological constraints and poor soil quality) resulted in a net species loss of 24% (±3) on average across regions and 37% (±3) in the worst-case situation (accounting of uncertainty of species loss modelling, Methods). While the biodiversity loss in percentage of the total species richness was similar across regions, the production gain varied remarkably. Arable and horticultural systems (AT, DE, FR and NL) would gain only 6% (±2) in production but would lose 26% (±7, worst case: 38% ± 6) of their species if 90% of the semi-natural habitats were converted to production fields. Grassland dominated regions would gain 47% (±14) in production due to conversion of the currently large available semi-natural area (CH, HU, NO and GB) and biodiversity loss would amount to 21% (±5, worst case: 36% ± 6) of their species. Amongst permanent crops, an intensification of low productive olive groves in Spain (ES) would increase production by 217% (but note that the bulk of semi-natural habitats in Spain are extensively managed olive groves, for which production was not accounted for, this means that the increased production is overestimated; Methods). The intensification of Spanish olive groves would induce a 27% (worst case: 37%) species loss. Vineyard farms in Italy (IT) would gain 3% production but would lose 22% (worst case: 41%) of the species.

On average, non-organic farms were 15.4% more productive than organic farms (Supplementary Table 6, Supplementary Data). On the other hand, organic farms, on average across regions, retained 1.8% (±11.2) and 5.9% (±7.1) more species of the four taxa combined than non-organic farms, when 50% and 90% of the semi-natural habitats were converted to productive

<figure><img src="image_5.png" /><figcaption>**Fig. 2** Estimated species richness and production for gradual conversion of semi-natural habitats to production fields in ten European farmland regions (brown: arable crops, mixed systems and horticulture, bright green: grassland, violet: permanent crops). The effect of the conversion of semi-natural habitats (percentage of the total farm area) on species richness (green lines, sum of vascular plant, earthworm, spider and wild bee species) is estimated with backward sample-size-based extrapolations³¹ of species unique to semi-natural habitats and addition of species unique to production fields and species shared (Methods). Confidence intervals (95%) are the sum over the confidence intervals of the four species groups. Effect on production (black lines) is calculated by multiplying the converted semi-natural area by the overall yield of the region. Effects of a 50% and 90% conversion of the semi-natural area are indicated for species richness (green dots) and production (black triangles) (Supplementary Table 7).</figcaption></figure>
<table><caption>Table 1 Average estimated species richness loss, gain and net change, and production gain by conversion of semi-natural habitats to production fields across ten European farmland regions.</caption><thead><tr><th></th><th>Group</th><th>Loss of species unique to semi-natural habitats (%)</th><th>Gain of species unique to production fields (%)</th><th>Species net change (%)</th><th>Production gain (%)</th></tr></thead><tbody><tr><td rowspan="10">Average situation</td><td rowspan="3">(I) 90% SNH converted</td><td>All regions</td><td>47.6 ± 3</td><td>24.1 ± 3.5</td><td>-23.6 ± 3.3</td><td>43.2 ± 21.1</td></tr><tr><td>arable</td><td>48.8 ± 4.9</td><td>23.2 ± 6.3</td><td>-25.7 ± 7.1</td><td>5.9 ± 1.8</td></tr><tr><td>grassland</td><td>44.8 ± 6.1</td><td>23.9 ± 6.5</td><td>-20.9 ± 5.2</td><td>47.2 ± 14.2</td></tr><tr><td rowspan="3">(II) 50% SNH converted</td><td>All regions</td><td>50.9 ± 4.9</td><td>26.4 ± 7.3</td><td>-24.6 ± 2.3</td><td>109.8 ± 107.1</td></tr><tr><td>arable</td><td>24.8 ± 2.6</td><td>23.1 ± 3.1</td><td>-1.7 ± 1.5</td><td>24 ± 11.7</td></tr><tr><td>grassland</td><td>24.7 ± 4.4</td><td>22.7 ± 6</td><td>-2 ± 3.4</td><td>3.3 ± 1</td></tr><tr><td rowspan="3">(III) 10% more production</td><td>All regions</td><td>23.5 ± 4.7</td><td>22.1 ± 5.2</td><td>-1.4 ± 2.2</td><td>26.2 ± 7.9</td></tr><tr><td>arable</td><td>27.5 ± 6.7</td><td>25.8 ± 7</td><td>-1.7 ± 0.2</td><td>61 ± 59.5</td></tr><tr><td>grassland</td><td>44.2 ± 9.4</td><td>23 ± 3.1</td><td>-21.2 ± 8.4</td><td>7.7 ± 1</td></tr><tr><td rowspan="3">Best case situation</td><td rowspan="3">(I) 90% SNH converted</td><td>All regions</td><td>65.5 ± 14</td><td>23.2 ± 6.3</td><td>-42.3 ± 12.6</td><td>5.9 ± 1.4</td></tr><tr><td>arable</td><td>22 ± 3.3</td><td>21.4 ± 4.2</td><td>-0.6 ± 2.3</td><td>10 ± 0</td></tr><tr><td>grassland</td><td>45.6 ± 29</td><td>25.6 ± 8.2</td><td>-20 ± 20.8</td><td>6.5 ± 3.5</td></tr><tr><td rowspan="10">Worst case situation</td><td rowspan="3">(I) 90% SNH converted</td><td>All regions</td><td>51.7 ± 3.8</td><td>42.3 ± 5.1</td><td>-9.4 ± 4.7</td><td>43.2 ± 21.1</td></tr><tr><td>arable</td><td>51.1 ± 5.9</td><td>41.7 ± 10</td><td>-9.4 ± 10.6</td><td>5.9 ± 1.8</td></tr><tr><td>grassland</td><td>50.3 ± 7.2</td><td>39.9 ± 8.5</td><td>-10.3 ± 6</td><td>47.2 ± 14.2</td></tr><tr><td rowspan="3">(II) 50% SNH converted</td><td>All regions</td><td>55.6 ± 10.3</td><td>48.3 ± 9.4</td><td>-7.4 ± 0.9</td><td>109.8 ± 107.1</td></tr><tr><td>arable</td><td>32 ± 3.5</td><td>39.7 ± 4.5</td><td>7.7 ± 1.9</td><td>24 ± 11.7</td></tr><tr><td>grassland</td><td>32.2 ± 6.8</td><td>39.7 ± 9.3</td><td>7.6 ± 3.5</td><td>3.3 ± 1</td></tr><tr><td rowspan="3">(III) 10% more production</td><td>All regions</td><td>36.6 ± 11.1</td><td>46.5 ± 8.6</td><td>6.7 ± 3.6</td><td>26.2 ± 7.9</td></tr><tr><td>arable</td><td>49.7 ± 13.8</td><td>38.8 ± 4.8</td><td>9.9 ± 2.5</td><td>61 ± 59.5</td></tr><tr><td>grassland</td><td>79.3 ± 20.7</td><td>41 ± 10.6</td><td>-10.9 ± 11.5</td><td>7.7 ± 1</td></tr><tr><td rowspan="3">Contrasting impact of setting a goal of 10% production increase</td><td rowspan="3">(I) 90% SNH converted</td><td>All regions</td><td>-53.8 ± 3.6</td><td>16.5 ± 1.9</td><td>-37.4 ± 3.1</td><td>43.2 ± 21.1</td></tr><tr><td>arable</td><td>-54 ± 5.9</td><td>16.4 ± 3.2</td><td>-37.6 ± 6.1</td><td>5.9 ± 1.8</td></tr><tr><td>grassland</td><td>-52 ± 7</td><td>15.5 ± 3.1</td><td>-36.4 ± 6</td><td>47.2 ± 14.2</td></tr><tr><td rowspan="3">(II) 50% SNH converted</td><td>All regions</td><td>-57.3 ± 7.8</td><td>18.5 ± 5.8</td><td>-38.8 ± 2</td><td>109.8 ± 107.1</td></tr><tr><td>arable</td><td>-24.5 ± 2.6</td><td>16.5 ± 1.9</td><td>-8 ± 2</td><td>24 ± 11.7</td></tr><tr><td>grassland</td><td>-23.9 ± 4</td><td>16.4 ± 3.2</td><td>-7.5 ± 2.8</td><td>3.3 ± 1</td></tr><tr><td rowspan="3">(III) 10% more production</td><td>All regions</td><td>-24.1 ± 4.8</td><td>15.5 ± 3.1</td><td>-8.6 ± 3.4</td><td>26.2 ± 7.9</td></tr><tr><td>arable</td><td>-26.4 ± 7.7</td><td>18.5 ± 5.8</td><td>-7.8 ± 1.8</td><td>61 ± 59.5</td></tr><tr><td>grassland</td><td>-57.4 ± 11.6</td><td>16.5 ± 1.9</td><td>-40.9 ± 10.8</td><td>7.7 ± 1</td></tr></tbody></table>

Percent and standard error are shown. Details per region are in Supplementary Table 7. Results are presented for three situations (see Methods): (1) Average situation: richness of unique species estimated from rarefaction and extrapolation curves; (2) Best case situation: species change estimated from the lower bound of the confidence interval (CI) of the expected species richness unique to semi-natural habitats, the upper bound of the CI of the expected species richness unique to production fields, and the lower bound of the CI of the expected shared species richness; (3) Worst case situation: species richness estimated from the upper bound of the confidence interval (CI) of the expected species richness unique to semi-natural habitats, the lower bound of the CI of the expected species richness unique to production fields, and shared species assumed to not be able to survive without semi-natural habitats and considered like species unique to semi-natural habitats. Three scenarios are shown: (I) 90% and (II) 50% of the area covered by semi-natural habitats are converted to production fields, (III) semi-natural habitats are converted to production fields required for a 10% production increase. Means are calculated across all regions and groups of regions, i.e. arable and horticultural (AT, FR, NL and DE), grassland (HU, CH, NO and GB) and permanent crops (IT, ES).

land, respectively. Still, the difference between both systems was not significant (Supplementary Table 8), partly due to the high variability between farms.

Contrasting impact of setting a goal of 10% production increase. The model allows to simulate whether certain production goals can be reached and what the impact on farmland species diversity would be. A report by the FAO³² postulates a global needed supply of food increase by almost 30% by 2030 and around 50% by 2050 to equal the rise in global demand. Here, we show the impact of a 10% increase in production in the regions studied, which corresponds to the level that about half of them could not reach, although they convert the entire semi-natural area to production fields.

The conversion to achieve the 10% production increase target has a very different impact depending on the type of land use. In arable systems, the required conversion would result in a

dramatic species loss of 66% (±16) in a worst case situation (accounting for the variability in the data as specified in Methods; average situation: 42% ± 13%) and considering a total conversion that could only be achieved in arable farmland with heavy investment in machinery. In that scenario, three of the four regions (AT, NL, DE) could not even reach 10% more production, because the share of semi-natural habitats is already so low, that not enough land can be converted (Table 1). Similarly, the vineyard region (IT) could not achieve the production goal and would lose up to 76% (average situation: 41%) of the species. In contrast, grassland dominated regions and olive groves would be able to produce 10% more than the current production while biodiversity would be less affected than for the arable systems (worst case situation: -14% ± 3 and -14%, average situation: -0.6% ± 2 and +0.8, respectively).

The organic farming sector is growing due to increased consumer demand. Yet, organic per hectare yield was significantly
<figure><img src="image_6.png" /></figure>

**Fig. 3** Comparative change in species richness (sum of vascular plant, earthworm, spider and bee species) by aiming at a 10% production increase in organic and non-organic systems of ten European farmland regions (brown: arable crops, mixed systems, and horticulture, bright green: grassland, violet: permanent crops). Net loss and gain of species (sum of vascular plant, earthworm, spider and bee species) by conversion of semi-natural habitats to production fields to achieve 10% more production than the non-organic level for non-organic systems (blue bars) and for organic systems (orange bars). Exclamation marks indicate a total conversion of semi-natural habitats without reaching 10% production increase. Species loss and gain are derived from sample size based extrapolation (Methods). Confidence intervals (95%) are the sum over the confidence intervals of the four species groups.

(15.4%, Supplementary Table 6) lower than non-organic on average across the study regions, which is in the lower range compared to previous studies²⁶,²⁷ (n = 876 fields, $Pχ^2_1 < 0.001$). In case of 10% more harvestable energy needed, non-organic and organic systems showed contrasting effects of the conversion of semi-natural habitats to production fields, non-organic systems losing more species (Fig. 3). Six non-organic and five organic systems, not paired in regions, out of 20 could successfully produce 10% more harvestable energy than today's non-organic farms by converting semi-natural habitats to production fields, and would lose 8% and 3% of species on average, respectively (Supplementary Table 9). These are notably the grassland systems in CH, NO, GB and the olive farms in ES. In those systems, the differences in management and yield between organic and non-organic farming is smaller than in arable farming and in vineyards, both of which tend to be more intensively managed. Also, the share of semi-natural habitats is higher in grassland than in arable systems, so that a conversion has less drastic consequences. Four non-organic and five organic systems could not achieve 10% more production despite conversion of the entire available semi-natural area, and would at the same time lose all species unique to semi-natural habitats, namely on average 65% and 48% of their species, respectively. The significant species gains in organic systems in NL and in non-organic systems in GB are due to the comparatively higher share of unique species in production fields in the two regions.

## Discussion

The key results from our large-scale study pertain to the confirmed essential contribution of semi-natural habitats to farmland biodiversity across all Europe, and a considerable species loss for a limited production gain in many farming systems if semi-

natural habitats are converted to production fields. They indicate that in half of the study regions—arable-dominated, horticultural and wine production—half of the wild species are restricted to semi-natural habitats that make up only 10% or less of the farms. The conversion of those small remnant semi-natural areas would not markedly increase production but impair species richness. By contrast, in grassland-dominated regions with large semi-natural areas, conversion would only marginally decrease species richness and substantially increase production. However, the assumption of achieving the same yield on such converted areas than on existing productive fields is probably overly optimistic as most of those grasslands would already have been intensified if soil and topographical conditions were favorable. Organic farms had higher species richness than non-organic farms in both semi-natural and production fields likely due to synergetic effects between both types of habitats. This allows organic systems to suffer lower losses in the event of a 10% increase in harvestable energy needed.

There are four major limitations to this analysis. First, we are aware that results may be affected by the specific regions, taxa considered and the data from only one year. Still, our models allow a first approximation of quantifying the trade-off between increased agricultural production and biodiversity loss for Europe. Application to biomes elsewhere with very large natural or semi-natural areas would be critical as biodiversity might be difficult to estimate with the conventional methods used here. The second limitation is that it assumes that production fields reclaimed from semi-natural habitats would have similar productivity as the already existing farmland. Yet, their productivity is likely to be lower due to constraints related to local soil conditions and topography. Therefore, the modelled yield increases are probably overestimated. The third limitation relates to the use of species richness as proxy for biodiversity. The conservation value of the species unique to semi-natural habitats is likely to be higher than the conservation value of species that are unique to farmland, which are mostly more common and widespread species. Therefore, the biodiversity loss in terms of conservation would probably be more severe than reflected by species numbers alone. The fourth limitation is that our model cannot account for the associated loss of ecosystem services provided by wild species (e.g. pollination and pest control), which would impair crop yields.

Despite those limitations, the analysis allows to approximate the detrimental consequences for farmland biodiversity in Europe of converting semi-natural habitats to farmland, in arable regions in particular. Potential increase of harvestable MJ energy production exists in grassland regions that would little impair biodiversity if kept reasonable, but this could hardly be in form of typical arable products such as cereals due to unfavorable topography, climate and soil conditions. Thus, although the harvestable MJ energy production may be balanced between regions, production types and modes can only be marginally transferred. The analysis supports the necessity to develop alternative solutions for sustainable food production in the future, by combining up-to-date technologies², agro-ecological approaches and best practices that minimize land use and preserve semi-natural land and biodiversity³³. Because the potential for production increase is limited, a sustainable European food system will also need to reduce food waste and require diets of the European population to adapt¹,³⁴,³⁵.

## Methods

**Study regions and farms.** Ten European regions from boreal to Mediterranean were selected (Supplementary Table 1). They represented major agricultural land uses such as arable crops including horticulture, mixed farming, grassland and perennial crops (vineyards and olives). Within each region, a pool of ~20–40 farms
was selected from which 12–20 farms were randomly selected (169 in total) that belonged to the same farm type, produced under homogeneous climatic and environmental circumstances and fulfilled specific criteria regarding their main production branch. In case the selected farms were not willing to participate, we asked other farms from the pool till the sufficient number has been reached. The selected organic farms had all been certified for at least five years. Farmers were asked if they were willing to participate in the study. If they refused, additional random sampling was conducted. In the region NL, 11 organic farms agreed to participate but only three non-organic farms, whereas seven organic farms and 11 non-organic farms were available in the region HU. During the study, one non-organic farmer in the region CH ceased participation.

**Habitat maps and farm interviews.** The complete area of all selected farms was mapped, using the BioHab method³⁶. Excluded from the farm area were woody and aquatic habitats larger than 800 m² and summer pastures. Within the farm area, areal and linear habitats were recorded. For an areal habitat, the minimal mapping unit was 400 m² with a width of at least 5 m. More narrow habitats, between 0.5 and 5 m wide and at least 30 m long, were mapped as linear habitats. Habitats were distinguished in habitat types according to Raunkier life forms, environmental conditions and management evidence²⁸. Further, a farmland class was assigned to each habitat that described whether the habitat was managed for agricultural production or other objectives such as e.g. nature conservation. In face-to-face interviews following a standardized questionnaire, farmers provided detailed information on field management and yield.

**Categorization as production fields and semi-natural habitats.** Based on the habitat maps and available information about management intensity, we categorized all habitats as either semi-natural habitats or production fields. In agricultural landscapes, these two categories are often not clearly distinguishable. There is a gradient from more intensively managed production fields to less intensively used semi-natural habitats. In addition, a categorization at the local scale can be different from an approach at a European scale (²⁹ and see p. 45 of³⁷). Here, we applied the same criteria for all ten study regions.

In all cases, we categorized as production fields: arable crops, intensively managed grasslands (following main plant species observed, management evidence and objectives, with fertilization and/or two or more cuts a year), horticultural crops, and vineyards.

We categorized as semi-natural habitats: linear habitats, habitats that were managed for nature conservation objectives, habitats where mainly geophytes, helophytes or hydrophytes were growing, grasslands with woody vegetation (shrubs and/or trees), and extensively managed grasslands (no fertilization, no or one cut a year).

**Species sampling.** Vascular plant, earthworm, spider and bee species were sampled in all different habitat types of a farm. One plot per habitat type was randomly selected per farm for species sampling. This resulted in 1402 selected habitat plots on 169 farms (Supplementary Table 2). In the selected habitats, species were sampled during one growing season, using standardized protocols¹⁹,³⁸. Plant species were identified in squares of 10 × 10 m in areal habitats and in rectangular strips of 1 × 10 m in linear habitats. Earthworms were collected at three random locations of 30 × 30 cm per habitat. First, a solution of allyl isothiocyanate (AITC) was poured out to extract earthworms from the soil. Afterwards, a 20-cm-deep soil core from the same location was hand sorted to find additional specimens. Identification took place in the lab. Spiders were sampled on three dates at five random locations per habitat within a circle of 0.1 m². Using a modified vacuum shredder, spiders were taken from the soil surface, transferred to a cool box, frozen, or put in ethanol, sorted and identified in the lab. Bees (wild bees and bumble bees) were sampled on three dates, during dry, sunny and warm weather conditions. They were captured with an entomological aerial net along a 100 m long and 2 m wide transect, transferred to a killing jar and identified in the lab.

**Grouping of species data.** Species data were pooled per taxa, habitat and region, and three sub-communities were formed with species (1) exclusively found in semi-natural habitats, i.e. unique to semi-natural habitats, (2) exclusively found in production fields, i.e. unique to production fields, and (3) found in both habitat categories i.e. shared by production fields and semi-natural habitats. For calculations of effects over all four taxa, species richness was the sum of the individual taxa species richnesses.

**Estimating species richness.** Species richness was estimated using coverage- and sample-size-based rarefaction and extrapolation curves³¹,³⁹,⁴⁰. Rarefaction and extrapolation, including confidence intervals (bootstrap method) and sampling coverage, were calculated in R 3.4.0⁴¹ using package iNEXT⁴². Detailed information is provided below for each topic.

**Estimating richness of unique species to compare semi-natural habitats and production fields.** To legitimately compare the richness of species unique to semi-natural habitats and to production fields, we used the coverage-based method, i.e.

we standardized the samples by their completeness³⁰. The point of comparison was determined by the so-called 'base coverage' identified by the following procedure³¹: (1) select the maximum sample coverage at reference sample size (number of sampling units) of the sub-communities under comparison, (2) select the minimum sample coverage at twice the reference sample size of the sub-communities under comparison, (3) identify the maximum of the results from step (1) and step (2) as 'base coverage'. The species richness estimates were then read off from the species sample-size-based rarefaction and extrapolation curves at the 'base coverage' for each sub-community being compared. If zero or exactly one species is unique to a sub-community at the reference sample size, no sample coverage could be calculated. In this case, we set the species richness at 0 or 1, respectively. The species richness estimate of the other sub-community under comparison was then read off at twice the reference sample size on the curve.

The 'base coverage' was individually defined for each region and each taxonomic group since the mixed effects models used to analyze the data took into account the variation among regions and taxonomic groups.

**Differences in species richness unique to semi-natural habitats and production fields.** The difference between the species richness unique to semi-natural habitats and unique to production fields was tested with mixed effects models using package lme4 (Version 1.1-12) in R⁴³. The data were ($S_{ij}$ | β, b, x) - Poisson($μ_{ij}$) from i = 1, ..., 10 regions. The model is:

$$ \ln(\mu_{ij}) = \beta_0 + \beta_1 x_{1i} + b_{1i} \quad (1) $$

$$ b_1 \sim N(0, \sigma^2) $$

where $β_0$ is a fixed intercept, $β_1$ a fixed effect sub-community $x_{1i}$ (species unique to semi-natural habitats versus species unique to production fields), $b_{1i}$ are random intercepts for region i. Random effects are normally distributed with mean 0 and variance $σ^2$. The significance of term $β_1$ was calculated by log-likelihood ratio tests with one degree of freedom. For the models over all four taxa, an additional random intercept was included, i.e. $b_{2i}$ with mean 0 and variance $σ^2$ for j = 1, ..., 4 taxa (Fig. 1b).

**Differences in species richness between organic and non-organic systems.** The comparison between organic and non-organic systems of species unique to semi-natural habitats and to production fields, and of species shared by the two habitat categories, relied on coverage-based extrapolation as described above. Differences between management systems were tested for significance using mixed-effects models with management system $β_1 x_{1ij}$ as fixed effect in (1).

**Estimating species loss due to conversion of semi-natural habitats to production fields.** To predict the species loss due to conversion of semi-natural habitats to production fields, we relied on sample-size-based extrapolations³¹ with species incidence frequencies. We estimated the richness of the species pool for the total number of mapped habitats including the extrapolated species richness unique to semi-natural habitats and unique to production fields, and the observed richness of shared species for each of the four taxa. This species pool provided the basis for the calculation of the species loss or gain (Table 1 and Supplementary Table 7). To model the species richness decrease for any amount of semi-natural habitats converted to production fields, we calculated and drew backward the curve composed of the accumulation curve for species unique to semi-natural habitats, to which the estimated total species richness unique to production fields (constant) and the corresponding gain of species unique to production fields (increases with increasing area of production fields as semi-natural habitats are converted), and the richness of observed shared species (constant) were added. This is the species decrease curve (Supplementary Fig. 2). If started at the observed species richness, this curve corresponds exactly to a species richness curve calculated by a cumulative random removal of semi-natural habitats one by one from the pool of all habitats. The four taxa decrease curves were added for the curve in Fig. 2. Confidence intervals (CI, 95%) shown in Figs. 2 and 3 are calculated by bootstrapping within the calculation of the species accumulation curves (iNEXT⁴²), upper and lower bounds of the 95% CI of the four taxa being added. From the species decrease curve, we read off the predicted species richness for a conversion of 50% and 90% of the semi-natural habitats, and a conversion required to increase production by 10%.

As species were sampled in 20% of all mapped habitats on average per region (min. 8%, max. 35%), extrapolated species accumulation curves used to build the species decrease curve were calculated for more than two to three times the reference sample size, which is the suggested range for reliable extrapolation of the species richness estimator³¹,⁴⁴. Obviously, the confidence intervals (CI) of the species richness extrapolations here became wide (Supplementary Fig. 4). As we still wanted to show the impact of a conversion of the whole semi-natural area into production fields on the production gain in the ten regions, we used the uncertainty (upper and lower bounds of the 95% CI of the four taxa added) to define two situations in addition to the average case to predict species richness for a 50% and a 90% semi-natural habitat conversion, and a conversion required to increase production by 10%: (1) a worst case situation with the upper bound of the CI of the
expected species richness unique to semi-natural habitats, the lower bound of the CI of the expected species richness unique to production fields, and shared species assumed not to be able to survive without semi-natural habitats and considered like species unique to semi-natural habitats (i.e. upper bound); and (2) a best case situation with the lower bound of the CI of the expected species richness unique to semi-natural habitats, the upper bound of the CI of the expected species richness unique to production fields, and the lower bound of the CI of the expected shared species richness.

**Estimating production gain.** Farmer interviews delivered an average yield per crop type per farm for the years 2008–2010 (Supplementary Data⁴⁵ shows details for organic and non-organic systems separately). Farmers indicated yield in kilograms or tons per hectare. This was transformed into energy units, i.e. mega joules per hectare (MJ ha⁻¹) using standard values⁴⁶. From this, for each region, the average yield (MJ ha⁻¹) was calculated by first multiplying individual crop type yields by the corresponding crop type areas to obtain the production per crop type, then summing up the production of all crop types, and finally dividing this sum by the total area of the crop types. For livestock farms, the fodder production of grasslands was estimated based on the average requirements per livestock unit, accounting for the amount of feed grain, legumes, silage maize and of imported feedstuff. All yields relate to plant biomass production and do not comprise livestock products. The average yield takes into account the relative cover of the different crop types in the regions. Therefore, the conversion of the semi-natural area to production fields was region-specific. The production of certain semi-natural habitats as e.g. olive groves in Spain was not part of the production calculation. The reason is that data on production for semi-natural habitats were mainly not available and/or negligible, e.g. extensively used grassland in CH or in HU, and we decided to apply the same treatment to all the regions. Consequently, in case of olive groves in Spain the effective increase in production is overestimated. To calculate the production gain per region, the production field area added by the conversion of semi-natural habitat area was multiplied by the average yield. In practice, in many regions it may be impossible to convert semi-natural habitat to productive land due to geomorphological constraints and poor soils, and even if land were converted, yields would be much lower than these averages. The results presented here, especially the 90% scenario, are therefore over-optimistic. On the other hand, our calculations are based on the area of semi-natural habitat available for conversion on existing farms, but in some regions other sources of semi-natural land may be available for conversion, e.g. former agricultural land that has been abandoned.

**Species loss and production gain for three scenarios.** We calculated the change of species richness and the production gain under current day production efficiency for two scenarios: (1) a conversion of 90% of the semi-natural area into production fields. The 10% of semi-natural area remaining is considered unsuitable for agricultural use or even impossible to cultivate; (2) a conversion of 50% of the semi-natural area into production fields, and (3) a necessary conversion of the semi-natural area into production fields to achieve a 10% production increase per region.

**Standardization for organic and non-organic systems.** Although the overall mapped area, the number of semi-natural habitats, the number of production fields and the average habitat size did not significantly differ between the two management systems (Supplementary Table 5), we standardized the number and size of habitats to the average across both systems per region to compare the species loss and production gain at current day production efficiency in the organic and non-organic systems. The total production in organic and non-organic systems per region was calculated based on the respective yield and the average mapped area of the production fields across both systems as described in section "Estimation of production gain". The impact on biodiversity was analyzed for the scenario that organic systems should achieve the same level of production as non-organic systems by converting semi-natural habitats to production fields. We calculated the amount of the required area to be converted into production fields and the corresponding species change.

Differences between management systems were again tested for significance using mixed-effects models with management system $\beta_1 x_{1ij}$ as fixed effect in (1).

**Reporting summary.** Further information on research design is available in the Nature Research Reporting Summary linked to this article.

## Data availability

Species, habitat and agricultural management data that support the findings of this study are available in Lüscher et al.¹⁹, (https://doi.org/10.1890/15-1985.1). The production data (yield) that support the findings of this study are provided as Supplementary Data 1 at Jeanneret et al.⁴⁵, (https://doi.org/10.5281/zenodo.5115742).

Received: 20 November 2020; Accepted: 5 August 2021;

Published online: 02 September 2021

## References

1. Tilman, D. & Clark, M. Global diets link environmental sustainability and human health. *Nature* **515**, 518–522 (2014).

2. Tilman, D., Balzer, C., Hill, J. & Befort, B. L. Global food demand and the sustainable intensification of agriculture. *Proc. Natl Acad. Sci. 108*, 20260–20264 (2011).

3. Cassidy, E. S., West, P. C., Gerber, J. S. & Foley, J. A. Redefining agricultural yields: from tonnes to people nourished per hectare. *Environ. Res. Letters* **8**, 034015 (2013).

4. Ray, D. K., Mueller, N. D., West, P. C. & Foley, J. A. Yield trends are insufficient to double global crop production by 2050. *PLoS ONE* **8**, e66428 (2013).

5. Foley, J. A. et al. Solutions for a cultivated planet. *Nature* **478**, 337–342 (2011).

6. Bajzelj, B. et al. Importance of food-demand management for climate mitigation. *Nature Clim. Change* **4**, 924–929 (2014).

7. Hertel, T. W. The global supply and demand for agricultural land in 2050: a perfect storm in the making? *Am. J. Agric. Econ.* **93**, 259–275 (2011).

8. Porter, J. R. et al. In *Climate change 2014: impacts, adaptation, and vulnerability. Part A: global and sectoral aspects. Contribution of working group II to the fifth assessment report of the intergovernmental panel on climate change* (eds Field, C. B. et al.) 485–533 (Cambridge University Press, 2014).

9. Challinor, A. J. et al. A meta-analysis of crop yield under climate change and adaptation. *Nature Clim. Change* **4**, 287–291 (2014).

10. FAO. World Agriculture Towards 2030/2050: The 2012 Revision. (FAO (Food and Agriculture Organisation), Rome, 2012).

11. Delzeit, R., Zabel, F., Meyer, C. & Václavík, T. Addressing future trade-offs between biodiversity and cropland expansion to improve food security. *Reg Environ Change*, 1–13 (2016).

12. EEA. Assessing biodiversity in Europe — the 2010 report 64 (EEA (European Environment Agency), Copenhagen, 2010).

13. NEAA & OECD. Background report to the OECD environmental outlook to 2030: overviews, details, and methodology of model-based analysis. Report No. 500113001, 184 (NEAA and OECD (Netherlands Environmental Assessment Agency and Organisation for Economic Co-operation and Development), Bilthoven, 2008).

14. Eurostat. Agriculture, forestry and fishery statistics, 2013 edition. 249 (Eurostat (European Union), Luxembourg, 2013).

15. Bignal, E. M. & McCracken, D. I. Low-intensity farming systems in the conservation of the countryside. *J. Appl. Ecol.* **33**, 413–424 (1996).

16. Tscharntke, T., Klein, A. M., Kruess, A., Steffan-Dewenter, I. & Thies, C. Landscape perspectives on agricultural intensification and biodiversity – ecosystem service management. *Ecol. Lett.* **8**, 857–874 (2005).

17. Sutherland, W. J. Conservation biology: openness in management. *Nature* **418**, 834–835 (2002).

18. Jongman, R. H. G. Homogenisation and fragmentation of the European landscape: ecological consequences and solutions. *Landscape Urban Plann.* **58**, 211–221 (2002).

19. Lüscher, G. et al. Farmland biodiversity and agricultural management on 237 farms in 13 European and two African regions. *Ecology* **97**, 1625–1625 (2016).

20. Godfray, H. C. J. et al. Food security: the challenge of feeding 9 billion people. *Science* **327**, 812–818 (2010).

21. Reganold, J. P. & Wachter, J. M. Organic agriculture in the twenty-first century. *Nat Plants* **2**, 15221 (2016).

22. Muller, A. et al. Strategies for feeding the world more sustainably with organic agriculture. *Nat Commun* **8**, 1290 (2017).

23. Winqvist, C. et al. Mixed effects of organic farming and landscape complexity on farmland biodiversity and biological control potential across Europe. *J. Appl. Ecol.* **48**, 570–579 (2011).

24. Schneider, M. K. et al. Gains to species diversity in organically farmed fields are not propagated at the farm level. *Nat Commun* **5**, https://doi.org/10.1038/ncomms5151 (2014).

25. Tuck, S. L. et al. Land-use intensity and the effects of organic farming on biodiversity: a hierarchical meta-analysis. *J. Appl. Ecol.* **51**, 746–755 (2014).

26. Seufert, V., Ramankutty, N. & Foley, J. A. Comparing the yields of organic and conventional agriculture. *Nature* **485**, 229–232 (2012).

27. Ponisio, L. C. et al. Diversification practices reduce organic to conventional yield gap. *Proc. R. Soc. Lond. B: Biol. Sci.* **282**, https://doi.org/10.1098/rspb.2014.1396 (2015).

28. Raunkjær, C. *The life forms of plants and statistical plant geography*. (Oxford: Clarendon Press., 1934).

29. Herzog, F. et al. European farm scale habitat descriptors for the evaluation of biodiversity. *Ecol. Indicators* **77**, 205–217 (2017).

30. Chao, A. & Jost, L. Coverage-based rarefaction and extrapolation: Standardizing samples by completeness rather than size. *Ecology* **93**, 2533–2547 (2012).

31. Chao, A. et al. Rarefaction and extrapolation with Hill numbers: a framework for sampling and estimation in species diversity studies. *Ecol. Monogr.* **84**, 45–67 (2014).
32. FAO. The future of food and agriculture – Alternative pathways to 2050. 224 (FAO (Food and Agriculture Organisation), Rome, 2018).

33. IEEP. Final Report – Reflecting environmental land use needs into EU policy: preserving and enhancing the environmental benefits of unfarmed features on EU Farmland., 234 (IEEP (DG Environment), London, 2008).

34. Wirsenius, S., Azar, C. & Berndes, G. How much land is needed for global food production under scenarios of dietary changes and livestock productivity increases in 2030? *Agr. Syst.* **103**, 621–638 (2010).

35. Baudry, J. et al. Improvement of diet sustainability with increased level of organic food in the diet: findings from the BioNutriNet cohort. *Am. J. Clin. Nutri.* **109**, 1173–1188 (2019).

36. Bunce, R. G. H. et al. A standardized procedure for surveillance and monitoring European habitats and provision of spatial data. *Landscape Ecol.* **23**, 11–25 (2008).

37. Herzog, F. et al. Biodiversity indicators for European Farming Systems. A Guidebook. (Forschungsanstalt Agroscope Reckenholz-Tänikon ART, Zürich, 2012).

38. Dennis, P., Bogers, M. M. B., Bunce, R. G. H., Herzog, F. & Jeanneret, P. *Biodiversity in organic and low-input farming systems. Handbook for recording key indicators*. (Alterra Wageningen, 2012).

39. Hill, M. O. Diversity and evenness: a unifying notation and its consequences. *Ecology* **54**, 427–432 (1973).

40. Jost, L. Entropy and diversity. *Oikos* **113**, 14 (2006).

41. R. Core Team. R: a language and environment for statistical computing. *R Foundation for Statistical Computing Vienna, Austria, https://www.R-project.org* (2017).

42. Hsieh, T. C., Ma, K. H. & Chao, A. iNEXT: an R package for rarefaction and extrapolation of species diversity (Hill numbers). *Method. Ecol. Evol.* **7**, 1451–1456 (2016).

43. Bates, D., Mächler, M., Bolker, B. & Walker, S. Fitting linear mixed-effects models using lme4. *J. Stat. Softw.* **67**, 48 (2015).

44. Colwell, R. K. et al. Models and estimators linking individual-based and sample-based rarefaction, extrapolation and comparison of assemblages. *J. Plant Ecol.* **5**, 3–21 (2012).

45. Jeanneret, Ph. et al. Detailed average yield per region, and organic and non-organic management systems per crop type in ten European farmland regions Data set. Zenodo. https://doi.org/10.5281/zenodo.5115742 (2021)

46. ADEME. Analyse énergétique d'exploitations agricoles et pouvoir de réchauffement global. Méthode et résultats sur 140 fermes françaises. 100 (ADEME, Dijon, 2002).

## Acknowledgements

We thank the farmers for access to land and for information on land management. We are indebted to S. Baldwin, O. Balle, M. Bátki, K.-G. Bernhardt, M. Bouvet, S. Buholzer, B. Bunce, C. Centeri, G. Cuendet, G. Engan, E. Falusi, O.-D. Finch, P. Gillingham, X. Heer, B. Heiner, M. Jerkovich, Z. Józan, C. Kantner, N. Koncz, A. Kulcsár, N. Kwikiriza, L. Lemaire, T. Lord, A. Mjelde, A. Müller, M. Münkenwarf, J. Nascimbene, J. Neumayr, F. Ødegaard, S. Papaja-Hülsbergen, M. Pavett, K. Penksza, L. Podmaniczky, R. Pommeresche, B. Putz, N. Richner, F. Schneiter, S. Schwarz, D. Schwenk, N. Stone,

O. Szalkovszki, G. Szalma, T. Szederjesi, H. Timmermann, S. Tolhurst, J. Vale, F. Vuillemin, A. Whittington, J. Wilkes and T. Zanetti for assistance in field sampling, taxonomy and data processing. We thank B. Schmid and P. Edwards for comments that helped improve this manuscript. This work was funded by the European Union through FP7 project BioBio (Indicators for biodiversity in organic and low-input farming systems; www.biobio-indicator.org; Agreement Nr. 227161), by the Austrian Ministry for Science and Research, and by the Lendület program of the Hungarian Academy of Sciences.

## Author contributions

F.H., P.J., D.B., P.D., W.F., J.K.F., R.H.G.J., M.K., G.M., M.G.P., P.P. and J.-P.S. conceived the project; M.A., D.B., K.B., A.B., J.-P.C., M.D., P.D., S.E., Z.E. W.F., T.F., J.K.F., I.R.G., P.G., T.G., P.J., G.L., G.J., R.H.G.J., M.K., A.K.-H., G.M., J.N., M.-L.O., M.G.P., J.-P.S., M.K.S., N.S., D.S., S.W. coordinated field sampling and data processing in the respective study region; P.J. and G.L. analysed the data, together with M.K.S., P.P. and F.H.; P.J. and G.L. wrote the paper with input from all the authors.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary information** The online version contains supplementary material available at https://doi.org/10.1038/s43247-021-00256-x.

**Correspondence** and requests for materials should be addressed to P.J.

**Peer review information** *Communications Earth & Environment* thanks the anonymous reviewers for their contribution to the peer review of this work. Primary Handling Editors: Clare Davis, Heike Langenberg. Peer reviewer reports are available.

**Reprints and permission information** is available at http://www.nature.com/reprints

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_7.png" /></figure>

**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 license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license 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 license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2021