A unifying modelling of multiple land degradation pathways in Europe.

Land degradation is a complex socio-environmental threat, which generally occurs as multiple concurrent pathways that remain largely unexplored in Europe. Here we present an unprecedented analysis of land multi-degradation in 40 continental countries, using twelve dataset-based processes that were modelled as land degradation convergence and combination pathways in Europe's agricultural (and arable) environments. Using a Land Multi-degradation Index, we find that up to 27%, 35% and 22% of continental agricultural (~2 million km2) and arable (~1.1 million km2) lands are currently threatened by one, two, and three drivers of degradation, while 10-11% of pan-European agricultural/arable landscapes are cumulatively affected by four and at least five concurrent processes. We also explore the complex pattern of spatially interacting processes, emphasizing the major combinations of land degradation pathways across continental and national boundaries. Our results will enable policymakers to develop knowledge-based strategies for land degradation mitigation and other critical European sustainable development goals.

Interaction effects of pH and land cover on soil microbial diversity are climate-dependent.

Factors regulating the diversity and composition of soil microbial communities include soil properties, land cover and climate. How these factors interact at large scale remains poorly investigated. Here, we used an extensive dataset including 715 locations from 24 European countries to investigate the interactive effects of climatic region, land cover and pH on soil bacteria and fungi. We found that differences in microbial diversity and community composition between land cover types depended on the climatic region. In Atlantic, Boreal and Continental regions, microbial richness was higher in croplands and grasslands than woodlands while richness in Mediterranean areas did not vary significantly among land cover types. These differences were further related to soil pH, as a driver of bacterial and fungal richness in most climatic regions, but the interaction of pH with land cover depended on the region. Microbial community composition differed the most between croplands and woodlands in all regions, mainly due to differences in pH. In the Mediterranean region, bacterial communities in woodlands and grasslands were the most similar, whereas in other regions, grassland and cropland-associated bacteria showed more similarity. Overall, we showed that key factors interact in shaping soil microbial communities in a climate-dependent way at large scale.

Soil organic carbon stocks in European croplands and grasslands: How much have we lost in the past decade?

The EU Soil Strategy 2030 aims to increase soil organic carbon (SOC) in agricultural land to enhance soil health and support biodiversity as well as to offset greenhouse gas emissions through soil carbon sequestration. Therefore, the quantification of current SOC stocks and the spatial identification of the main drivers of SOC changes is paramount in the preparation of agricultural policies aimed at enhancing the resilience of agricultural systems in the EU. In this context, changes of SOC stocks (Δ SOCs) for the EU + UK between 2009 and 2018 were estimated by fitting a quantile generalized additive model (qGAM) on data obtained from the revisited points of the Land Use/Land Cover Area Frame Survey (LUCAS) performed in 2009, 2015 and 2018. The analysis of the partial effects derived from the fitted qGAM model shows that land use and land use change observed in the 2009, 2015 and 2018 LUCAS campaigns (i.e. continuous grassland [GGG] or cropland [CCC], conversion grassland to cropland (GGC or GCC) and vice versa [CGG or CCG]) was one of the main drivers of SOC changes. The CCC was the factor that contributed to the lowest negative change on Δ SOC with an estimated partial effect of -0.04 ± 0.01 g C kg-1 year-1 , while the GGG the highest positive change with an estimated partial effect of 0.49 ± 0.02 g C kg-1 year-1 . This confirms the C sequestration potential of converting cropland to grassland. However, it is important to consider that local soil and environmental conditions may either diminish or enhance the grassland's positive effect on soil C storage. In the EU + UK, the estimated current (2018) topsoil (0-20 cm) SOC stock in agricultural land below 1000 m a.s.l was 9.3 Gt, with a Δ SOC of -0.75% in the period 2009-2018. The highest estimated SOC losses were concentrated in central-northern countries, while marginal losses were observed in the southeast.

Cadmium in topsoils of the European Union - An analysis based on LUCAS topsoil database.

Cadmium (Cd) is a naturally occurring element that can accumulate in the soil through the application of fertilisers containing cadmium and as a waste of industrial processes. Cadmium inputs in the soil have increased significantly (+50 %) during the 20th century as a result of the application of fertilisers and sewage sludge, and also due to local contamination (e.g. waste dumping, mining) and industrial emissions (e.g. zinc smelters). Using the 21,682 soil samples from the LUCAS soil survey, we aim to estimate the spatial distribution of the concentration of Cd in the European Union (EU) and UK topsoil. Out of the total, 72.6 % of the samples have Cd values <0.07 mg kg-1, 21.6 % in the range 0.07-1 mg kg-1 and the remaining 5.5 % higher than the threshold of 1 mg kg-1, which is generally considered the limit for risk assessment. The mean Cd value in the EU topsoils is 0.20 mg kg-1, slightly higher in grasslands (0.24 mg kg-1) compared to croplands (0.17 mg kg-1). Applying an ensemble of machine learning models supported by a variety of environmental descriptors, we created maps of Cd distribution at a resolution of 100 m. The ensemble approach included five models and increased the prediction accuracy to R2 of 0.45 (an increase of 0.1 compared to best single model performance). The approach used resulted in a high predictive power for the general Cd distribution, while also identifying hotspots of Cd contamination. Natural factors influencing Cd levels include soil properties (pH, clay), topography, soil erosion, and leaching. As anthropogenic factors, we identified phosphorus inputs to agricultural lands as the most important for Cd levels. The application of the EU Fertiliser Directive should further limit Cd inputs and potentially the Cd content in soils.

Global rainfall erosivity database (GloREDa) and monthly R-factor data at 1 km spatial resolution.

Here, we present and release the Global Rainfall Erosivity Database (GloREDa), a multi-source platform containing rainfall erosivity values for almost 4000 stations globally. The database was compiled through a global collaboration between a network of researchers, meteorological services and environmental organisations from 65 countries. GloREDa is the first open access database of rainfall erosivity (R-factor) based on hourly and sub-hourly rainfall records at a global scale. This database is now stored and accessible for download in the long-term European Soil Data Centre (ESDAC) repository of the European Commission's Joint Research Centre. This will ensure the further development of the database with insertions of new records, maintenance of the data and provision of a helpdesk. In addition to the annual erosivity data, this release also includes the mean monthly erosivity data for 94% of the GloREDa stations. Based on these mean monthly R-factor values, we predict the global monthly erosivity datasets at 1 km resolution using the ensemble machine learning approach (ML) as implemented in the mlr package for R. The produced monthly raster data (GeoTIFF format) may be useful for soil erosion prediction modelling, sediment distribution analysis, climate change predictions, flood, and natural disaster assessments and can be valuable inputs for Land and Earth Systems modelling.

Ecosystem type drives soil eukaryotic diversity and composition in Europe.

Soil eukaryotes play a crucial role in maintaining ecosystem functions and services, yet the factors driving their diversity and distribution remain poorly understood. While many studies focus on some eukaryotic groups (mostly fungi), they are limited in their spatial scale. Here, we analyzed an unprecedented amount of observational data of soil eukaryomes at continental scale (787 sites across Europe) to gain further insights into the impact of a wide range of environmental conditions (climatic and edaphic) on their community composition and structure. We found that the diversity of fungi, protists, rotifers, tardigrades, nematodes, arthropods, and annelids was predominantly shaped by ecosystem type (annual and permanent croplands, managed and unmanaged grasslands, coniferous and broadleaved woodlands), and higher diversity of fungi, protists, nematodes, arthropods, and annelids was observed in croplands than in less intensively managed systems, such as coniferous and broadleaved woodlands. Also in croplands, we found more specialized eukaryotes, while the composition between croplands was more homogeneous compared to the composition of other ecosystems. The observed high proportion of overlapping taxa between ecosystems also indicates that DNA has accumulated from previous land uses, hence mimicking the land transformations occurring in Europe in the last decades. This strong ecosystem-type influence was linked to soil properties, and particularly, soil pH was driving the richness of fungi, rotifers, and annelids, while plant-available phosphorus drove the richness of protists, tardigrades, and nematodes. Furthermore, the soil organic carbon to total nitrogen ratio crucially explained the richness of fungi, protists, nematodes, and arthropods, possibly linked to decades of agricultural inputs. Our results highlighted the importance of long-term environmental variables rather than variables measured at the time of the sampling in shaping soil eukaryotic communities, which reinforces the need to include those variables in addition to ecosystem type in future monitoring programs and conservation efforts.

Spatial assessment of topsoil zinc concentrations in Europe.

Zinc (Zn) is essential to sustain crop production and human health, while it can be toxic when present in excess. In this manuscript, we applied a machine learning model on 21,682 soil samples from the Land Use and Coverage Area frame Survey (LUCAS) topsoil database of 2009/2012 to assess the spatial distribution in Europe of topsoil Zn concentrations measured by aqua regia extraction, and to identify the influence of natural drivers and anthropogenic sources on topsoil Zn concentrations. As a result, a map was produced showing topsoil Zn concentrations in Europe at a resolution of 250 m. The mean predicted Zn concentration in Europe was 41 mg kg-1, with a root mean squared error of around 40 mg kg-1 calculated for independent soil samples. We identified clay content as the most important factor explaining the overall distribution of soil Zn in Europe, with lower Zn concentrations in coarser soils. Next to texture, low Zn concentrations were found in soils with low pH (e.g. Podzols), as well as in soils with pH above 8 (i.e., Calcisols). The presence of deposits and mining activities mainly explained the occurrence of relatively high Zn concentrations above 167 mg kg-1 (the one percentile highest concentrations) within 10 km from these sites. In addition, the relatively higher Zn levels found in grasslands in regions with high livestock density may point to manure as a significant source of Zn in these soils. The map developed in this study can be used as a reference to assess the eco-toxicological risks associated with soil Zn concentrations in Europe and areas with Zn deficiency. In addition, it can provide a baseline for future policies in the context of pollution, soil health, human health, and crop nutrition.

Patterns in soil microbial diversity across Europe.

Factors driving microbial community composition and diversity are well established but the relationship with microbial functioning is poorly understood, especially at large scales. We analysed microbial biodiversity metrics and distribution of potential functional groups along a gradient of increasing land-use perturbation, detecting over 79,000 bacterial and 25,000 fungal OTUs in 715 sites across 24 European countries. We found the lowest bacterial and fungal diversity in less-disturbed environments (woodlands) compared to grasslands and highly-disturbed environments (croplands). Highly-disturbed environments contain significantly more bacterial chemoheterotrophs, harbour a higher proportion of fungal plant pathogens and saprotrophs, and have less beneficial fungal plant symbionts compared to woodlands and extensively-managed grasslands. Spatial patterns of microbial communities and predicted functions are best explained when interactions among the major determinants (vegetation cover, climate, soil properties) are considered. We propose guidelines for environmental policy actions and argue that taxonomical and functional diversity should be considered simultaneously for monitoring purposes.

Improving the phosphorus budget of European agricultural soils.

Despite phosphorus (P) being crucial for plant nutrition and thus food security, excessive P fertilization harms soil and aquatic ecosystems. Accordingly, the European Green Deal and derived strategies aim to reduce P losses and fertilizer consumption in agricultural soils. The objective of this study is to calculate a soil P budget, allowing the quantification of the P surpluses/deficits in the European Union (EU) and the UK, considering the major inputs (inorganic fertilizers, manure, atmospheric deposition, and chemical weathering) and outputs (crop production, plant residues removal, losses by erosion) for the period 2011-2019. The Land Use/Cover Area frame Survey (LUCAS) topsoil data include measured values for almost 22,000 samples for both available and total P. With advanced machine learning models, we developed maps for both attributes at 500 m resolution. We estimated the available P for crops at a mean value of 83 kg ha-1 with a clear distinction between North and South. The ratio of available P to the total P is about 1:17. The inorganic fertilizers and manure contribute almost equally as P inputs (mean 16 ± 2 kg P ha-1 yr-1 at 90 % confidence level) to agricultural soils, with high regional variations depending on farming practices, livestock density, and cropping systems. The P outputs came mainly from the exportation by the harvest of crop products and residues (97.5 %) and, secondly, by erosion. Using a sediment distribution model, we quantified the P fluxes to river basins and sea outlets. In the EU and UK, we estimated an average surplus of 0.8 kg P ha-1 yr-1 with high variability between countries with some regional variations. The P annual budget at regional scale showed ample possibility to improve P management by both reducing inputs in regions with high surplus (and P soil available) and rebalancing fertilization in those at risk of soil fertility depletion.

GloSEM: High-resolution global estimates of present and future soil displacement in croplands by water erosion.

Healthy soil is the foundation underpinning global agriculture and food security. Soil erosion is currently the most serious threat to soil health, leading to yield decline, ecosystem degradation and economic impacts. Here, we provide high-resolution (ca. 100 × 100 m) global estimates of soil displacement by water erosion obtained using the Revised-Universal-Soil-Loss-Equation-based Global Soil Erosion Modelling (GloSEM) platform under present (2019) and future (2070) climate scenarios (i.e. Shared Socioeconomic Pathway [SSP]1-Representative Concentration Pathway [RCP]2.6, SSP2-RCP4.5 and SSP5-RCP8.5). GloSEM is the first global modelling platform to take into account regional farming systems, the mitigation effects of conservation agriculture (CA), and climate change projections. We provide a set of data, maps and descriptive statistics to support researchers and decision-makers in exploring the extent and geography of soil erosion, identifying probable hotspots, and exploring (with stakeholders) appropriate actions for mitigating impacts. In this regard, we have also provided an Excel spreadsheet that can provide useful insights into the potential mitigating effects of present and future alternative CA scenarios at the country level.

Mercury in European topsoils: Anthropogenic sources, stocks and fluxes.

Mercury (Hg) is one of the most dangerous pollutants worldwide. In the European Union (EU), we recently estimated the Hg distribution in topsoil using 21,591 samples and a series of geo-physical inputs. In this manuscript, we investigate the impact of mining activities, chrol-alkali industries and other diffuse pollution sources as primary anthropogenic sources of Hg hotspots in the EU. Based on Hg measured soil samples, we modelled the Hg pool in EU topsoils, which totals about 44.8 Gg, with an average density of 103 g ha-1. As a following step, we coupled the estimated Hg stocks in topsoil with the pan-European assessment of soil loss due to water erosion and sediment distribution. In the European Union and UK, we estimated that about 43 Mg Hg yr-1 are displaced by water erosion and c. a. 6 Mg Hg yr-1 are transferred with sediments to river basins and eventually released to coastal Oceans. The Mediterranean Sea receives almost half (2.94 Mg yr-1) of the Hg fluxes to coastal oceans and it records the highest quantity of Hg sediments. This is the result of elevated soil Hg concentration and high erosion rates in the catchments draining into the Mediterranean Sea. This work contributes to new knowledge in support of the policy development in the EU on the Zero Pollution Action Plan and the Sustainable Development Goal (SDGs) 3.9 and 14.1, which both have as an objective to reduce soil pollution by 2030.

Soil erosion modelling: A bibliometric analysis.

Soil erosion can present a major threat to agriculture due to loss of soil, nutrients, and organic carbon. Therefore, soil erosion modelling is one of the steps used to plan suitable soil protection measures and detect erosion hotspots. A bibliometric analysis of this topic can reveal research patterns and soil erosion modelling characteristics that can help identify steps needed to enhance the research conducted in this field. Therefore, a detailed bibliometric analysis, including investigation of collaboration networks and citation patterns, should be conducted. The updated version of the Global Applications of Soil Erosion Modelling Tracker (GASEMT) database contains information about citation characteristics and publication type. Here, we investigated the impact of the number of authors, the publication type and the selected journal on the number of citations. Generalized boosted regression tree (BRT) modelling was used to evaluate the most relevant variables related to soil erosion modelling. Additionally, bibliometric networks were analysed and visualized. This study revealed that the selection of the soil erosion model has the largest impact on the number of publication citations, followed by the modelling scale and the publication's CiteScore. Some of the other GASEMT database attributes such as model calibration and validation have negligible influence on the number of citations according to the BRT model. Although it is true that studies that conduct calibration, on average, received around 30% more citations, than studies where calibration was not performed. Moreover, the bibliographic coupling and citation networks show a clear continental pattern, although the co-authorship network does not show the same characteristics. Therefore, soil erosion modellers should conduct even more comprehensive review of past studies and focus not just on the research conducted in the same country or continent. Moreover, when evaluating soil erosion models, an additional focus should be given to field measurements, model calibration, performance assessment and uncertainty of modelling results. The results of this study indicate that these GASEMT database attributes had smaller impact on the number of citations, according to the BRT model, than anticipated, which could suggest that these attributes should be given additional attention by the soil erosion modelling community. This study provides a kind of bibliographic benchmark for soil erosion modelling research papers as modellers can estimate the influence of their paper.

Soil erosion modelling: A global review and statistical analysis.

To gain a better understanding of the global application of soil erosion prediction models, we comprehensively reviewed relevant peer-reviewed research literature on soil-erosion modelling published between 1994 and 2017. We aimed to identify (i) the processes and models most frequently addressed in the literature, (ii) the regions within which models are primarily applied, (iii) the regions which remain unaddressed and why, and (iv) how frequently studies are conducted to validate/evaluate model outcomes relative to measured data. To perform this task, we combined the collective knowledge of 67 soil-erosion scientists from 25 countries. The resulting database, named 'Global Applications of Soil Erosion Modelling Tracker (GASEMT)', includes 3030 individual modelling records from 126 countries, encompassing all continents (except Antarctica). Out of the 8471 articles identified as potentially relevant, we reviewed 1697 appropriate articles and systematically evaluated and transferred 42 relevant attributes into the database. This GASEMT database provides comprehensive insights into the state-of-the-art of soil- erosion models and model applications worldwide. This database intends to support the upcoming country-based United Nations global soil-erosion assessment in addition to helping to inform soil erosion research priorities by building a foundation for future targeted, in-depth analyses. GASEMT is an open-source database available to the entire user-community to develop research, rectify errors, and make future expansions.

A spatial assessment of mercury content in the European Union topsoil.

Mapping of surface soil Hg concentrations, a priority pollutant, at continental scale is important in order to identify hotspots of soil Hg distribution (e.g. mining or industrial pollution) and identify factors that influence soil Hg concentrations (e.g. climate, soil properties, vegetation). Here we present soil Hg concentrations from the LUCAS topsoil (0-20 cm) survey including 21,591 samples from 26 European Union countries (one sample every ~200 km2). Deep Neural Network (DNN) learning models were used to map the European soil Hg distribution. DNN estimated a median Hg concentration of 38.3 µg kg-1 (2.6 to 84.7 µg kg-1) excluding contaminated sites. At continental scale, we found that soil Hg concentrations increased with latitude from south to north and with altitude. A GLMM revealed a correlation (R2 = 0.35) of soil Hg concentrations with vegetation activity, normalized difference vegetation index (NDVI), and soil organic carbon content. This observation corroborates the importance of atmospheric Hg0 uptake by plants and the build-up of the soil Hg pool by litterfall over continental scales. The correlation of Hg concentrations with NDVI was amplified by higher soil organic matter content, known to stabilize Hg in soils through thiol bonds. We find a statistically significant relation between soil Hg levels and coal use in large power plants, proving that emissions from power plants are associated with higher mercury deposition in their proximity. In total 209 hotspots were identified, defined as the top percentile in Hg concentration (>422 µg kg-1). 87 sites (42% of all hotspots) were associated with known mining areas. The sources of the other hotspots could not be identified and may relate to unmined geogenic Hg or industrial pollution. The mapping effort in the framework of LUCAS can serve as a starting point to guide local and regional authorities in identifying Hg contamination hotspots in soils.

Global phosphorus shortage will be aggravated by soil erosion.

Soil phosphorus (P) loss from agricultural systems will limit food and feed production in the future. Here, we combine spatially distributed global soil erosion estimates (only considering sheet and rill erosion by water) with spatially distributed global P content for cropland soils to assess global soil P loss. The world's soils are currently being depleted in P in spite of high chemical fertilizer input. Africa (not being able to afford the high costs of chemical fertilizer) as well as South America (due to non-efficient organic P management) and Eastern Europe (for a combination of the two previous reasons) have the highest P depletion rates. In a future world, with an assumed absolute shortage of mineral P fertilizer, agricultural soils worldwide will be depleted by between 4-19 kg ha-1 yr-1, with average losses of P due to erosion by water contributing over 50% of total P losses.

Land use and climate change impacts on global soil erosion by water (2015-2070).

Soil erosion is a major global soil degradation threat to land, freshwater, and oceans. Wind and water are the major drivers, with water erosion over land being the focus of this work; excluding gullying and river bank erosion. Improving knowledge of the probable future rates of soil erosion, accelerated by human activity, is important both for policy makers engaged in land use decision-making and for earth-system modelers seeking to reduce uncertainty on global predictions. Here we predict future rates of erosion by modeling change in potential global soil erosion by water using three alternative (2.6, 4.5, and 8.5) Shared Socioeconomic Pathway and Representative Concentration Pathway (SSP-RCP) scenarios. Global predictions rely on a high spatial resolution Revised Universal Soil Loss Equation (RUSLE)-based semiempirical modeling approach (GloSEM). The baseline model (2015) predicts global potential soil erosion rates of [Formula: see text] Pg yr-1, with current conservation agriculture (CA) practices estimated to reduce this by ∼5%. Our future scenarios suggest that socioeconomic developments impacting land use will either decrease (SSP1-RCP2.6-10%) or increase (SSP2-RCP4.5 +2%, SSP5-RCP8.5 +10%) water erosion by 2070. Climate projections, for all global dynamics scenarios, indicate a trend, moving toward a more vigorous hydrological cycle, which could increase global water erosion (+30 to +66%). Accepting some degrees of uncertainty, our findings provide insights into how possible future socioeconomic development will affect soil erosion by water using a globally consistent approach. This preliminary evidence seeks to inform efforts such as those of the United Nations to assess global soil erosion and inform decision makers developing national strategies for soil conservation.

Plutonium aided reconstruction of caesium atmospheric fallout in European topsoils.

Global nuclear weapon testing and the Chernobyl accident have released large amounts of radionuclides into the environment. However, to date, the spatial patterns of these fallout sources remain poorly constrained. Fallout radionuclides (137Cs, 239Pu, 240Pu) were measured in soil samples (n = 160) collected at flat, undisturbed grasslands in Western Europe in the framework of a harmonised European soil survey. We show that both fallout sources left a specific radionuclide imprint in European soils. Accordingly, we used plutonium to quantify contributions of global versus Chernobyl fallout to 137Cs found in European soils. Spatial prediction models allowed for a first assessment of the global versus Chernobyl fallout pattern across national boundaries. Understanding the magnitude of these fallout sources is crucial not only to establish a baseline in case of future radionuclide fallout but also to define a baseline for geomorphological reconstructions of soil redistribution due to soil erosion processes.

Mapping LUCAS topsoil chemical properties at European scale using Gaussian process regression.

This paper presents the second part of the mapping of topsoil properties based on the Land Use and Cover Area frame Survey (LUCAS). The first part described the physical properties (Ballabio et al., 2016) while this second part includes the following chemical properties: pH, Cation Exchange Capacity (CEC), calcium carbonates (CaCO3), C:N ratio, nitrogen (N), phosphorus (P) and potassium (K). The LUCAS survey collected harmonised data on changes in land cover and the state of land use for the European Union (EU). Among the 270,000 land use and cover observations selected for field visit, approximately 20,000 soil samples were collected in 24 EU Member States in 2009 together with more than 2000 samples from Bulgaria and Romania in 2012. The chemical properties maps for the European Union were produced using Gaussian process regression (GPR) models. GPR was selected for its capacity to assess model uncertainty and the possibility of adding prior knowledge in the form of covariance functions to the model. The derived maps will establish baselines that will help monitor soil quality and provide guidance to agro-environmental research and policy developments in the European Union.

Soil erosion is unlikely to drive a future carbon sink in Europe.

Understanding of the processes governing soil organic carbon turnover is confounded by the fact that C feedbacks driven by soil erosion have not yet been fully explored at large scale. However, in a changing climate, variation in rainfall erosivity (and hence soil erosion) may change the amount of C displacement, hence inducing feedbacks onto the land C cycle. Using a consistent biogeochemistry-erosion model framework to quantify the impact of future climate on the C cycle, we show that C input increases were offset by higher heterotrophic respiration under climate change. Taking into account all the additional feedbacks and C fluxes due to displacement by erosion, we estimated a net source of 0.92 to 10.1 Tg C year-1 from agricultural soils in the European Union to the atmosphere over the period 2016-2100. These ranges represented a weaker and stronger C source compared to a simulation without erosion (1.8 Tg C year-1), respectively, and were dependent on the erosion-driven C loss parameterization, which is still very uncertain. However, when setting a baseline with current erosion rates, the accelerated erosion scenario resulted in 35% more eroded C, but its feedback on the C cycle was marginal. Our results challenge the idea that higher erosion driven by climate will lead to a C sink in the near future.