Search

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 R 2 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., Competing Interests: Declaration of competing interest The authors declare the following financial interests/personal relationships which may be considered as potential competing interests: Panos Panagos reports was provided by European Commission Joint Research Centre Ispra. Panos Panagos reports a relationship with European Commission Joint Research Centre Ispra that includes: employment. None, (Copyright © 2023 The Authors. Published by Elsevier B.V. All rights reserved.)

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., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 The Authors. Published by Elsevier B.V. All rights reserved.)

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., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2022 The Authors. Published by Elsevier B.V. All rights reserved.)

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., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2021 The Authors. Published by Elsevier B.V. All rights reserved.)

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 km 2 ). 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 (R 2  = 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 Hg 0 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., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2021 The Authors. Published by Elsevier B.V. All rights reserved.)

Copper (Cu) distribution in soil is influenced by climatic, geological and pedological factors. Apart from geological sources and industrial pollution, other anthropogenic sources, related to the agricultural activity, may increase copper levels in soils, especially in permanent crops such as olive groves and vineyards. This study uses 21,682 soil samples from the LUCAS topsoil survey to investigate copper distribution in the soils of 25 European Union (EU) Member States. Generalized Linear Models (GLM) were used to investigate the factors driving copper distribution in EU soils. Regression analysis shows the importance of topsoil properties, land cover and climate in estimating Cu concentration. Meanwhile, a copper regression model confirms our hypothesis that different agricultural management practices have a relevant influence on Cu concentration. Besides the traditional use of copper as a fungicide for treatments in several permanent crops, the combined effect of soil properties such as high pH, soil organic carbon and clay, with humid and wet climatic conditions favours copper accumulation in soils of vineyards and tree crops. Compared to the overall average Cu concentration of 16.85 mg kg -1 , vineyards have the highest mean soil Cu concentration (49.26 mg kg -1 ) of all land use categories, followed by olive groves and orchards. Gaussian Process Regression (GPR) combined with kriging were used to map copper concentration in topsoils and to evidence the presence of outliers. GPR proved to be performant in predicting Cu concentration, especially in combination with kriging, accounting for 66% of Cu deviance. The derived maps are novel as they include information about the importance of topsoil properties in the copper mapping process, thus improving its accuracy. Both models highlight the influence of land management practices in copper concentration and the strong correlation between topsoil copper and vineyards., (Copyright © 2018 The Authors. Published by Elsevier B.V. All rights reserved.)

Rainfall erosivity as a dynamic factor of soil loss by water erosion is modelled intra-annually for the first time at European scale. The development of Rainfall Erosivity Database at European Scale (REDES) and its 2015 update with the extension to monthly component allowed to develop monthly and seasonal R-factor maps and assess rainfall erosivity both spatially and temporally. During winter months, significant rainfall erosivity is present only in part of the Mediterranean countries. A sudden increase of erosivity occurs in major part of European Union (except Mediterranean basin, western part of Britain and Ireland) in May and the highest values are registered during summer months. Starting from September, R-factor has a decreasing trend. The mean rainfall erosivity in summer is almost 4 times higher (315MJmmha -1 h -1 ) compared to winter (87MJmmha -1 h -1 ). The Cubist model has been selected among various statistical models to perform the spatial interpolation due to its excellent performance, ability to model non-linearity and interpretability. The monthly prediction is an order more difficult than the annual one as it is limited by the number of covariates and, for consistency, the sum of all months has to be close to annual erosivity. The performance of the Cubist models proved to be generally high, resulting in R 2 values between 0.40 and 0.64 in cross-validation. The obtained months show an increasing trend of erosivity occurring from winter to summer starting from western to Eastern Europe. The maps also show a clear delineation of areas with different erosivity seasonal patterns, whose spatial outline was evidenced by cluster analysis. The monthly erosivity maps can be used to develop composite indicators that map both intra-annual variability and concentration of erosive events. Consequently, spatio-temporal mapping of rainfall erosivity permits to identify the months and the areas with highest risk of soil loss where conservation measures should be applied in different seasons of the year., (Copyright © 2016 British Geological Survey, NERC. Published by Elsevier B.V. All rights reserved.)

Because of the increasing pressures exerted on soil, below-ground life is under threat. Knowledge-based rankings of potential threats to different components of soil biodiversity were developed in order to assess the spatial distribution of threats on a European scale. A list of 13 potential threats to soil biodiversity was proposed to experts with different backgrounds in order to assess the potential for three major components of soil biodiversity: soil microorganisms, fauna, and biological functions. This approach allowed us to obtain knowledge-based rankings of threats. These classifications formed the basis for the development of indices through an additive aggregation model that, along with ad-hoc proxies for each pressure, allowed us to preliminarily assess the spatial patterns of potential threats. Intensive exploitation was identified as the highest pressure. In contrast, the use of genetically modified organisms in agriculture was considered as the threat with least potential. The potential impact of climate change showed the highest uncertainty. Fourteen out of the 27 considered countries have more than 40% of their soils with moderate-high to high potential risk for all three components of soil biodiversity. Arable soils are the most exposed to pressures. Soils within the boreal biogeographic region showed the lowest risk potential. The majority of soils at risk are outside the boundaries of protected areas. First maps of risks to three components of soil biodiversity based on the current scientific knowledge were developed. Despite the intrinsic limits of knowledge-based assessments, a remarkable potential risk to soil biodiversity was observed. Guidelines to preliminarily identify and circumscribe soils potentially at risk are provided. This approach may be used in future research to assess threat at both local and global scale and identify areas of possible risk and, subsequently, design appropriate strategies for monitoring and protection of soil biota., (Copyright © 2015 The Authors. Published by Elsevier B.V. All rights reserved.)

Recently, in the Auerswald et al. (2015) comment on "Rainfall erosivity in Europe", 5 criticisms were addressed: i) the neglect of seasonal erosion indices, ii) the neglect of published studies and data, iii) the low temporal resolution of the data, especially of the maximum rain intensity, iv) the use of precipitation data instead of rain data and the subsequent deviation of the R-factor in Germany and Austria compared with previous studies, and v) the differences in considered time periods between countries. We reply as follows: (i) An evaluation of the seasonal erosion index at the European scale is, to our knowledge, not achievable at present with the available data but would be a future goal. Synchronous publication of the seasonal erosion index is not mandatory, specifically because seasonal soil loss ratios are not available at this scale to date. We are looking forward to the appropriate study by the authors of the comment, who assert that they have access to the required data. (ii) We discuss and evaluate relevant studies in our original work and in this reply; however, we cannot consider what is not available to the scientific community. (iii) The third point of critique was based on a misunderstanding by Auerswald et al. (2015), as we did indeed calculate the maximum intensity with the highest resolution of data available. (iv) The low R-factor values in Germany and the higher values in Austria compared with previous studies are not due to the involvement of snow but are rather due to a Pan-European interpolation. We argue that an interpolation across the borders of Austria creates a more reliable data set. (v) We agree that the use of a short time series or time series from different periods is generally a problem in all large-scale studies and requires improvement in the future. However, because this affects countries with a rather low variability of the R-factor in our study, we are confident that the overall results of the map are not biased. In conclusion, the Pan-European rainfall data compilation (REDES) was a great success and yielded data from 1541 stations with an average length of 17.1years and a temporal resolution of <60min. However, a Pan-European data collection will never be complete without the help and supply of data from its users. Thus, we invite the authors of the comment to share their data in the open REDES to help build even better rainfall-erosivity maps at regional or European scales., (Copyright © 2015 Elsevier B.V. All rights reserved.)

Rainfall is one the main drivers of soil erosion. The erosive force of rainfall is expressed as rainfall erosivity. Rainfall erosivity considers the rainfall amount and intensity, and is most commonly expressed as the R-factor in the USLE model and its revised version, RUSLE. At national and continental levels, the scarce availability of data obliges soil erosion modellers to estimate this factor based on rainfall data with only low temporal resolution (daily, monthly, annual averages). The purpose of this study is to assess rainfall erosivity in Europe in the form of the RUSLE R-factor, based on the best available datasets. Data have been collected from 1541 precipitation stations in all European Union (EU) Member States and Switzerland, with temporal resolutions of 5 to 60 min. The R-factor values calculated from precipitation data of different temporal resolutions were normalised to R-factor values with temporal resolutions of 30 min using linear regression functions. Precipitation time series ranged from a minimum of 5 years to a maximum of 40 years. The average time series per precipitation station is around 17.1 years, the most datasets including the first decade of the 21st century. Gaussian Process Regression (GPR) has been used to interpolate the R-factor station values to a European rainfall erosivity map at 1 km resolution. The covariates used for the R-factor interpolation were climatic data (total precipitation, seasonal precipitation, precipitation of driest/wettest months, average temperature), elevation and latitude/longitude. The mean R-factor for the EU plus Switzerland is 722 MJ mm ha(-1) h(-1) yr(-1), with the highest values (>1000 MJ mm ha(-1) h(-1) yr(-1)) in the Mediterranean and alpine regions and the lowest (<500 MJ mm ha(-1) h(-1) yr(-1)) in the Nordic countries. The erosivity density (erosivity normalised to annual precipitation amounts) was also the highest in Mediterranean regions which implies high risk for erosive events and floods., (Copyright © 2015 The Authors. Published by Elsevier B.V. All rights reserved.)

The greatest obstacle to soil erosion modelling at larger spatial scales is the lack of data on soil characteristics. One key parameter for modelling soil erosion is the soil erodibility, expressed as the K-factor in the widely used soil erosion model, the Universal Soil Loss Equation (USLE) and its revised version (RUSLE). The K-factor, which expresses the susceptibility of a soil to erode, is related to soil properties such as organic matter content, soil texture, soil structure and permeability. With the Land Use/Cover Area frame Survey (LUCAS) soil survey in 2009 a pan-European soil dataset is available for the first time, consisting of around 20,000 points across 25 Member States of the European Union. The aim of this study is the generation of a harmonised high-resolution soil erodibility map (with a grid cell size of 500 m) for the 25 EU Member States. Soil erodibility was calculated for the LUCAS survey points using the nomograph of Wischmeier and Smith (1978). A Cubist regression model was applied to correlate spatial data such as latitude, longitude, remotely sensed and terrain features in order to develop a high-resolution soil erodibility map. The mean K-factor for Europe was estimated at 0.032 thahha(-1)MJ(-1)mm(-1) with a standard deviation of 0.009 thahha(-1)MJ(-1)mm(-1). The yielded soil erodibility dataset compared well with the published local and regional soil erodibility data. However, the incorporation of the protective effect of surface stone cover, which is usually not considered for the soil erodibility calculations, resulted in an average 15% decrease of the K-factor. The exclusion of this effect in K-factor calculations is likely to result in an overestimation of soil erosion, particularly for the Mediterranean countries, where highest percentages of surface stone cover were observed., (Copyright © 2014. Published by Elsevier B.V.)

A reliable spatial assessment of the POPs contamination in soils is essential for burden studies and flux evaluations. Soil characteristics and properties vary enormously even within small spatial scale and over time; therefore soil capacity of accumulating POPs varies greatly. In order to include this very high spatial and temporal variability, models can be used for assessing soil accumulation capacity in a specific time and space and, from it, the spatial distribution and temporal trends of POPs concentrations. In this work, predictive contamination maps of the accumulation capacity of soils were developed at a space resolution of 1×1m with a time frame of one day, in a study area located in the central Alps. Physical algorithms for temperature and organic carbon estimation along the soil profile and across the year were fitted to estimate the horizontal, vertical and seasonal distribution of the contamination potential for PCBs in soil (Ksa maps). The resulting maps were cross-validated with an independent set of PCB contamination data, showing very good agreement (e.g. for CB-153, R(2)=0.80, p-value≤2.2·10(-06)). Slopes of the regression between predicted Ksa and experimental concentrations were used to map the soil contamination for the whole area, taking into account soil characteristics and temperature conditions. These maps offer the opportunity to evaluate burden (concentration maps) and fluxes (emission maps) with highly resolved temporal and spatial detail. In addition, in order to explain the observed low autumn PCB concentrations in soil related to the high Ksa values of this period, a dynamic model of seasonal variation of soil concentrations was developed basing on rate parameters fitted on measured concentrations. The model was able to describe, at least partially, the observed different behavior between the quite rapid discharge phase in summer and the slow recharge phase in autumn., (Copyright © 2013 Elsevier B.V. All rights reserved.)

Under the European Union Thematic Strategy for Soil Protection, the European Commission Directorate-General for the Environment and the European Environmental Agency (EEA) identified a decline in soil organic carbon and soil losses by erosion as priorities for the collection of policy relevant soil data at European scale. Moreover, the estimation of soil organic carbon content is of crucial importance for soil protection and for climate change mitigation strategies. Soil organic carbon is one of the attributes of the recently developed LUCAS soil database. The request for data on soil organic carbon and other soil attributes arose from an on-going debate about efforts to establish harmonized datasets for all EU countries with data on soil threats in order to support modeling activities and display variations in these soil conditions across Europe. In 2009, the European Commission's Joint Research Centre conducted the LUCAS soil survey, sampling ca. 20,000 points across 23 EU member states. This article describes the results obtained from analyzing the soil organic carbon data in the LUCAS soil database. The collected data were compared with the modeled European topsoil organic carbon content data developed at the JRC. The best fitted comparison was performed at NUTS2 level and showed underestimation of modeled data in southern Europe and overestimation in the new central eastern member states. There is a good correlation in certain regions for countries such as the United Kingdom, Slovenia, Italy, Ireland, and France. Here we assess the feasibility of producing comparable estimates of the soil organic carbon content at NUTS2 regional level for the European Union (EU27) and draw a comparison with existing modeled data. In addition to the data analysis, we suggest how the modeled data can be improved in future updates with better calibration of the model., (Copyright © 2012 Elsevier B.V. All rights reserved.)

Polycyclic Aromatic Hydrocarbons (PAHs) are a major group of pollutants whose occurrence in the environment is mainly of anthropogenic origin. In this paper, we examine the effect of topographical slope exposure on PAH contamination and seasonal change in PAH concentrations in soils. We collected soil samples on three dates in 2007 (early May, end of July and beginning of November) from south- and north-facing aspects at 1900 m a.s.l. in the central Italian Alps. We found greater PAH contamination in soils from a north-facing slope than in those from a south-facing slope at all seasons. We calculated North-South Enrichment Factors as the ratio between the concentrations measured in soils from northern and southern aspects. These ratios ranged from 1.4 to 1.9 for lighter PAHs (from 2 to 4 rings). These values are consistent with theoretical calculations based on temperature-specific octanol-air partition coefficients (predicted North-South Enrichment Factors range from 1.6 to 2.0). For heavier PAHs (from 5 to 6 rings), smaller differences were observed between soils from northern and southern aspects, due to the gas/particle distribution of these compounds. We also found consistent differences in normalised PAH concentrations across the three sampling periods. The majority of compounds showed a significant decreasing trend from the beginning of May to the end of July, due to the annual cycles of physical processes (deposition vs. volatilisation) and biological processes (uptake and/or biotransformation). Only a few compounds showed different trends, presumably due to season-specific local emission sources.

Statistical techniques can be used in groundwater pollution problems to determine the relationships among observed contamination (impacted wells representing an occurrence of what has to be predicted), environmental factors that may influence it and the potential contamination sources. Determination of a threshold concentration to discriminate between impacted or non impacted wells represents a key issue in the application of these techniques. In this work the effects on groundwater vulnerability assessment by statistical methods due to the use of different threshold values have been evaluated. The study area (Province of Milan, northern Italy) is about 2000 km(2) and groundwater nitrate concentration is constantly monitored by a net of about 300 wells. Along with different predictor factors three different threshold values of nitrate concentration have been considered to perform the vulnerability assessment of the shallow unconfined aquifer. The likelihood ratio model has been chosen to analyze the spatial distribution of the vulnerable areas. The reliability of the three final vulnerability maps has been tested showing that all maps identify a general positive trend relating mean nitrate concentration in the wells and vulnerability classes the same wells belong to. Then using the kappa coefficient the influence of the different threshold values has been evaluated comparing the spatial distribution of the resulting vulnerability classes in each map. The use of different threshold does not determine different vulnerability assessment if results are analyzed on a broad scale, even if the smaller threshold value gives the poorest performance in terms of reliability. On the contrary, the spatial distribution of a detailed vulnerability assessment is strongly influenced by the selected threshold used to identify the occurrences, suggesting that there is a strong relationship among the number of identified occurrences, the scale of the maps representing the predictor factors and the model efficiency in discriminating different vulnerable areas.