A step towards a holistic assessment of soil degradation in Europe: Coupling on-site erosion with sediment transfer and carbon fluxes

The boom era of emerging contaminants: A review of remediating agricultural soils by biochar

Emerging contaminants (ECs) are widely sourced persistent pollutants that pose a significant threat to the environment and human health. Their footprint spans global ecosystems, making their remediation highly challenging. In recent years, a significant amount of literature has focused on the use of biochar for remediation of heavy metals and organic pollutants in soil and water environments. However, the use of biochar for the remediation of ECs in agricultural soils has not received as much attention, and as a result, there are limited reviews available on this topic. Thus, this review aims to provide an overview of the primary types, sources, and hazards of ECs in farmland, as well as the structure, functions, and preparation types of biochar. Furthermore, this paper emphasizes the importance and prospects of three remediation strategies for ECs in cropland: (i) employing activated, modified, and composite biochar for remediation, which exhibit superior pollutant removal compared to pure biochar; (ii) exploring the potential synergistic efficiency between biochar and compost, enhancing their effectiveness in soil improvement and pollution remediation; (iii) utilizing biochar as a shelter and nutrient source for microorganisms in biochar-mediated microbial remediation, positively impacting soil properties and microbial community structure. Given the increasing global prevalence of ECs, the remediation strategies provided in this paper aim to serve as a valuable reference for future remediation of ECs-contaminated agricultural lands.

Into the unknown: The role of post-fire soil erosion in the carbon cycle

Wildfires directly emit 2.1 Pg carbon (C) to the atmosphere annually. The net effect of wildfires on the C cycle, however, involves many interacting source and sink processes beyond these emissions from combustion. Among those, the role of post-fire enhanced soil organic carbon (SOC) erosion as a C sink mechanism remains essentially unquantified. Wildfires can greatly enhance soil erosion due to the loss of protective vegetation cover and changes to soil structure and wettability. Post-fire SOC erosion acts as a C sink when off-site burial and stabilization of C eroded after a fire, together with the on-site recovery of SOC content, exceed the C losses during its post-fire transport. Here we synthesize published data on post-fire SOC erosion and evaluate its overall potential to act as longer-term C sink. To explore its quantitative importance, we also model its magnitude at continental scale using the 2017 wildfire season in Europe. Our estimations show that the C sink ability of SOC water erosion during the first post-fire year could account for around 13% of the C emissions produced by wildland fires. This indicates that post-fire SOC erosion is a quantitatively important process in the overall C balance of fires and highlights the need for more field data to further validate this initial assessment.

Including land management in a European carbon model with lateral transfer to the oceans

The use of cover crops (CCs) is a promising cropland management practice with multiple benefits, notably in reducing soil erosion and increasing soil organic carbon (SOC) storage. However, the current ability to represent these factors in land surface models remains limited to small scales or simplified and lumped approaches due to the lack of a sediment-carbon erosion displacement scheme. This precludes a thorough understanding of the consequences of introducing a CC into agricultural systems. In this work, this problem was addressed in two steps with the spatially distributed CE-DYNAM model. First, the historical effect of soil erosion, transport, and deposition on the soil carbon budget at a continental scale in Europe was characterized since the early industrial era, using reconstructed climate and land use forcings. Then, the impact of two distinct policy-oriented scenarios for the introduction of CCs were evaluated, covering the European cropping systems where surface erosion rates or nitrate susceptibility are critical. The evaluation focused on the increase in SOC storage and the export of particulate organic carbon (POC) to the oceans, compiling a continental-scale carbon budget. The results indicated that Europe exported 1.95 TgC/year of POC to the oceans in the last decade, and that CCs can contribute to reducing this amount while increasing SOC storage. Compared to the simulation without CCs, the additional rate of SOC storage induced by CCs peaked after 10 years of their adoption, followed by a decrease, and the cumulative POC export reduction stabilized after around 13 years. The findings indicate that the impacts of CCs on SOC and reduced POC export are persistent regardless of their spatial allocation adopted in the scenarios. Together, the results highlight the importance of taking the temporal aspect of CC adoption into account and indicate that CCs alone are not sufficient to meet the targets of the 4‰ initiative. Despite some known model limitations, which include the lack of feedback of erosion on the net primary productivity and the representation of carbon fluxes with an emulator, the current work constitutes the first approach to successfully couple a distributed routing scheme of eroded carbon to a land carbon model emulator at a reasonably high resolution and continental scale. SHORT ABSTRACT: A spatially distributed model coupling erosion, transport, and deposition to the carbon cycle was developed. Then, it was used to simulate the impact of cover crops on both erosion and carbon, to show that cover crops can simultaneously increase organic carbon storage and reduce particulate organic carbon export to the oceans. The results seemed persistent regardless of the spatial distribution of cover crops.

Assessing the phosphorus cycle in European agricultural soils: Looking beyond current national phosphorus budgets

Phosphorus (P) is an essential nutrient for all crops, yet its excess negatively affects public health, the environment, and the economy. At the same time, rock P is a critical raw material due to its importance for food production, the finite geological deposits, and its unequal regional distribution. As a consequence, nutrient management is addressed by numerous environmental policies. Process-based biogeochemical models are valuable instruments to monitor the P cycle and predict the effect of agricultural management policies. In this study, we upscale the calibrated DayCent model at European level using data-derived soil properties, advanced input data sets, and representative management practices. Our results depicted a P budget with an average P surplus (0.11 kg P ha-1 year-1), a total soil P (2240.0 kg P ha-1), and available P content (77.4 kg P ha-1) consistent with literature and national statistics. Through agricultural management scenarios, we revealed a range of potential changes in the P budget by 2030 and 2050, influenced by the interlink of P with biogeochemical carbon and nitrogen cycles. Thus, we developed a powerful assessment tool capable of i) identifying areas with P surplus or deficit at high spatial resolution of 1 km2, (ii) pinpointing areas where a change in agricultural management would be most urgent to reach policy goals in terms of environmental pollution, food security and resource efficiency of a critical raw material, and iii) assessing the response of the P cycle to modifications in agricultural management.

EUSEDcollab: a network of data from European catchments to monitor net soil erosion by water

As a network of researchers we release an open-access database (EUSEDcollab) of water discharge and suspended sediment yield time series records collected in small to medium sized catchments in Europe. EUSEDcollab is compiled to overcome the scarcity of open-access data at relevant spatial scales for studies on runoff, soil loss by water erosion and sediment delivery. Multi-source measurement data from numerous researchers and institutions were harmonised into a common time series and metadata structure. Data reuse is facilitated through accompanying metadata descriptors providing background technical information for each monitoring station setup. Across ten European countries, EUSEDcollab covers over 1600 catchment years of data from 245 catchments at event (11 catchments), daily (22 catchments) and monthly (212 catchments) temporal resolution, and is unique in its focus on small to medium catchment drainage areas (median = 43 km2, min = 0.04 km2, max = 817 km2) with applicability for soil erosion research. We release this database with the aim of uniting people, knowledge and data through the European Union Soil Observatory (EUSO).

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.

Assessment of current reservoir sedimentation rate and storage capacity loss: An Italian overview

Sedimentation has a prominent impact on the functionality and lifetime of reservoirs and is a growing concern for stakeholders. Various parameters influence sedimentation caused by soil erosion. Here we have examined fifty Italian reservoirs to determine sedimentation rates and storage capacity loss. The reservoirs studied have an average age of 78 years as of 2021, with the highest loss of capacity observed, equal to 100%, for Ceppo Morelli. For the fifty Italian catchments covering north, south, central and islands of Italy, we found the mean annual sediment yield varying between 17-4000 m³/km². year. Six of fifty reservoirs studied (Quarto, Colombara, Ceppo Morelli, Fusino, Vodo and Valle di Cadore) are already in a very critical situation in terms of storage capacity loss. Out of the fifty reservoirs, half of them will reach their half-life year by 2050. For example, for the Fusino reservoir located in northern Italy, we observed a loss of 90% of the storage volume as of 2020 with respect to its operation year 1974, compared to 6% in 2015 as available in literature. Modelling the sediment delivery ratio (SDR) is an open question, due to the lack of adequate data and uncertainties about the variability in hydrological, geomorphological, climate and landcover parameters. Here, we addressed the issue with a simplified multiple regression approach based on sediment delivery ratio values retrieved by the RUSLE model. We found different multi regressions for reservoirs belonging to the Alpine and Apennine regions. This analysis offers a starting point for the management and prioritization of adaptation and remediation policies necessary to address reservoir sedimentation.

Sources of copper into the European aquatic environment

Chemical contamination from point source discharges in developed (resource-rich) countries has been widely regulated and studied for decades; however, diffuse sources are largely unregulated and widespread. In the European Union (EU), large dischargers report releases of some chemicals, yet little is known of total emissions (point and diffuse) and their relative significance. We estimated copper loadings from all significant sources including industry, sewage treatment plants, surface runoff (from traffic, architecture, and atmospheric deposition), septic tanks, agriculture, mariculture, marine transport (antifoulant leaching), and natural processes. A combination of European datasets, literature, and industry data were used to generate export coefficients. These were then multiplied by activity rates to derive loads. A total of approximately 8 kt of copper per annum (ktpa) is estimated to enter freshwaters in the EU, and another 3.5 ktpa enters transitional and coastal waters. The main inputs to freshwater are natural processes (3.7 ktpa), agriculture (1.8 ktpa), and runoff (1.8 ktpa). Agricultural emissions are dominated by copper-based plant protection products and farmyard manure. Urban runoff is influenced by copper use in architecture and by vehicle brake linings. Antifoulant leaching from boats (3.2 ktpa) dominates saline water loads of copper. It is noteworthy that most of the emissions originate in a limited number of copper uses where environmental exposure and pathways exist, compared with the bulk of copper use within electrical and electronic equipment and infrastructure that has no environmental pathway during its use. A sensitivity analysis indicated significant uncertainty in data from abandoned mines and urban runoff load estimates. This study provided for the first time a methodology and comprehensive metal load apportionment to European aquatic systems, identifying data gaps and uncertainties, which may be refined over time. Source apportionments using this methodology can inform more cost-effective environmental risk assessment and management. Integr Environ Assess Manag 2022;00:1-17. © 2022 The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental Toxicology & Chemistry (SETAC).

The Effects of Soil Improving Cropping Systems (SICS) on Soil Erosion and Soil Organic Carbon Stocks across Europe: A Simulation Study

Healthy soils are fundamental for sustainable agriculture. Soil Improving Cropping Systems (SICS) aim to make land use and food production more sustainable. To evaluate the effect of SICS at EU scale, a modelling approach was taken. This study simulated the effects of SICS on two principal indicators of soil health (Soil Organic Carbon stocks) and land degradation (soil erosion) across Europe using the spatially explicit PESERA model. Four scenarios with varying levels and combinations of cover crops, mulching, soil compaction alleviation and minimum tillage were implemented and simulated until 2050. Results showed that while in the scenario without SICS, erosion slightly increased on average across Europe, it significantly decreased in the scenario with the highest level of SICS applied, especially in the cropping areas in the central European Loess Belt. Regarding SOC stocks, the simulations show a substantial decrease for the scenario without SICS and a slight overall decrease for the medium level scenario and the scenario with a mix of high, medium and no SICS. The scenario with a high level of SICS implementation showed an overall increase in SOC stocks across Europe. Potential future improvements include incorporating dynamic land use, climate change and an optimal spatial allocation of SICS.

Urbanisation-driven land degradation and socioeconomic challenges in peri-urban areas: Insights from Southern Europe

Climate change and landscape transformation have led to rapid expansion of peri-urban areas globally, representing new 'laboratories' for the study of human-nature relationships aiming at land degradation management. This paper contributes to the debate on human-driven land degradation processes by highlighting how natural and socioeconomic forces trigger soil depletion and environmental degradation in peri-urban areas. The aim was to classify and synthesise the interactions of urbanisation-driven factors with direct or indirect, on-site or off-site, and short-term or century-scale impacts on land degradation, focussing on Southern Europe as a paradigmatic case to address this issue. Assuming complex and multifaceted interactions among influencing factors, a relevant contribution to land degradation was shown to derive from socioeconomic drivers, the most important of which were population growth and urban sprawl. Viewing peri-urban areas as socio-environmental systems adapting to intense socioeconomic transformations, these factors were identified as forming complex environmental 'syndromes' driven by urbanisation. Based on this classification, we suggested three key measures to support future land management in Southern European peri-urban areas.

Soil-Improving Cropping Systems for Sustainable and Profitable Farming in Europe

Soils form the basis for agricultural production and other ecosystem services, and soil management should aim at improving their quality and resilience. Within the SoilCare project, the concept of soil-improving cropping systems (SICS) was developed as a holistic approach to facilitate the adoption of soil management that is sustainable and profitable. SICS selected with stakeholders were monitored and evaluated for environmental, sociocultural, and economic effects to determine profitability and sustainability. Monitoring results were upscaled to European level using modelling and Europe-wide data, and a mapping tool was developed to assist in selection of appropriate SICS across Europe. Furthermore, biophysical, sociocultural, economic, and policy reasons for (non)adoption were studied. Results at the plot/farm scale showed a small positive impact of SICS on environment and soil, no effect on sustainability, and small negative impacts on economic and sociocultural dimensions. Modelling showed that different SICS had different impacts across Europe—indicating the importance of understanding local dynamics in Europe-wide assessments. Work on adoption of SICS confirmed the role economic considerations play in the uptake of SICS, but also highlighted social factors such as trust. The project’s results underlined the need for policies that support and enable a transition to more sustainable agricultural practices in a coherent way.

Using WaTEM/SEDEM to Model the Effects of Crop Rotation and Changes in Land Use on Sediment Transport in the Vrchlice Watershed

The Czech landscape has undergone various changes over the last 100 years and has been mainly adapted agriculturally for economic purposes. This has resulted, among other things, in reservoirs being clogged with sediment. The Vrchlice Reservoir was built in 1970 to supply drinking water for around 50,000 inhabitants, and increased sedimentation has been detected in the reservoir in recent years. Water erosion and sediment transport were modeled with WaTEM/SEDEM. Sediment volumes were measured in eight ponds across the watershed for calibration purposes. Modeled results from ponds in watersheds covered mostly with arable lands generally corresponded with the measured values. Although in forested watersheds, the measured sediment volumes greatly exceeded modeled sediment yields, indicating high uncertainty in using USLE-based models in non-agricultural watersheds. The modeled scenarios represented pre-Communist, Communist, and post-Communist eras. For these periods WaTEM/SEDEM was used to evaluate three isolated effects: the effects of various crops on arable lands, the effects of farmland fragmentation, and finally the effects of changes in land use. The change in crops proved to be an important factor causing high siltation rate (potential 23% reduction in sediment yield for historical periods), and land fragmentation played the second important role (potential 15% reduction in sediment yield can be reached by land fragmentation). Across all scenarios, the lowest sediment yield and reservoirs siltation rates were obtained from the pre-Communist and Communist crop share under current land use conditions, and current land use with farmland fragmentation implemented, as it was re-constructed for the pre-Communist era. This supports the idea that the introduction of green areas within arable lands are beneficial to the landscape and can help reduce soil erosion and reservoir siltation.

Spatial modeling of soil organic carbon using remotely sensed indices and environmental field inventory variables

The relationship between soil organic carbon (SOC) and environmental parameters was investigated in the Galazchai Watershed, Iran. Therefore, correlating the SOC amounts with remote sensing (RS) indices, topographic variables, and soil texture was analyzed. Some 125 soil samples gather from the upper 30 cm, and the weight of each sample was about 0.5 kg. The RS indices, consisting of difference vegetation index (DVI), enhanced vegetation index (EVI), optimized soil adjusted vegetation index (OSAVI), normalized difference vegetation index (NDVI), and soil adjusted vegetation index (SAVI), were used. Topographic variables included slope, elevation, aspect, and topographical wetness index (TWI), as well as clay and silt contents. The ordinary least square (OLS) and the geographically weighted regression (GWR) were employed to develop the SOC relationship considering different combinations of the variables. Results showed that none of the combinations of variables accurately estimated SOC (R² < 0.32 and p value > 0.001). However, EVI with GWR (R² = 0.291) and OSAVI, clay, slope, and aspect with GWR (R²= 0.32) better estimated SOC. Therefore, results showed that the study remotely sensed indices and environmental field inventory variables could not favorably predict the SOC content. These results can be attributed to the low SOC values varying from 0.917 to 3.355%, with a mean of 2.194 ± 0.522 in the study watershed. However, studies using more uniformly distributed and denser sampling in the study area and other methods to investigate the relationship between variables are recommended.

Soil degradation in the European Mediterranean region: Processes, status and consequences

Soil, a non-renewable resource, sustains life on Earth by supporting around 95% of global food production and providing ecosystem services such as biomass production, filtration of contaminants and transfer of mass and energy between spheres. Unsustainable management practices and climate change are threatening the natural capital of soils, particularly in the Mediterranean region, where increasing population, rapid land-use changes, associated socio-economic activities and climate change are imposing high pressures on the region's shallow soils. Despite evidence of high soil susceptibility to degradation and desertification, the true extent of soil degradation in the region is unknown. This paper reviews and summarises the scientific literature and relevant official reports, with the aim to advance this knowledge by synthesizing, mapping, and identifying gaps regarding the status, causes, and consequences of soil degradation processes in the European Mediterranean region. This is needed as scientific underpinning of efforts to counteract soil degradation in the region. Three main degradation categories are then considered: physical (soil sealing, compaction, erosion), chemical (soil organic matter, contamination, salinisation), and biological. We find some degradation processes to be relatively well-documented (e.g. soil erosion), while others, such as loss of biodiversity, remain poorly addressed, with limited data availability. We suggest establishment of a continuous, harmonised soil monitoring system at national and regional scale in the Mediterranean region to provide comparable datasets and chart the spatial extent and temporal changes in soil degradation, and corresponding economic implications. This is critical to support decision-making and fulfilment of related sustainable development goals.

Concentration of organic carbon in Finnish catchments and variables involved in its variations

The majority of the carbon worldwide is in soil. In a river catchment, the tight relationship between soil, water and climate makes carbon likely to be eroded and transported from the soil to the rivers. There are multiple variables which can trigger and accelerate the process. In order to assess the importance of the factors involved, and their interactions resulting in the changes in the carbon cycle within catchments, we have studied the catchments of 26 Finnish rivers from 2000 to 2019. These catchments are distributed all over Finland, but we have grouped them into three categories: southern, peatland and northern. We have run a boosted regression tree (BRT) analysis on chemical, physical, climatic and anthropogenic factors to determine their influence on the variations of total organic carbon (TOC) concentration. TOC concentration has decreased in Finland between 2000 and 2019 by 0.91 mg/l, driven principally by forest ditching and % old forest in the catchment. Old forest is especially dominant in the northern catchments with an influence on TOC of 40.5%. In southern and peatland catchments, average precipitation is an important factor to explain the changes in TOC whilst in northern catchments, organic fields have more influence.

Use of Monitoring Approaches to Verify the Predictive Accuracy of the Modeling of Particle-Bound Solid Inputs to Surface Waters

For particle-bound substances such as phosphorus, erosion is an important input pathway to surface waters. Therefore, knowledge of soil erosion by water and sediment inputs to water bodies at high spatial resolution is essential to derive mitigation measures at the regional scale. Models are used to calculate soil erosion and associated sediment inputs to estimate the resulting loads. However, validation of these models is often not sufficiently possible. In this study, sediment input was modeled on a 10 × 10 m grid for a subcatchment of the Kraichbach river in Baden-Wuerttemberg (Germany). In parallel, large-volume samplers (LVS) were operated at the catchment outlet, which allowed a plausibility check of the modeled sediment inputs. The LVS produced long-term composite samples (2 to 4 weeks) over a period of 4 years. The comparison shows a very good agreement between the modeled and measured sediment loads. In addition, the monitoring concept of the LVS offers the possibility to identify the sources of the sediment inputs to the water body. In the case of the Kraichbach river, it was found that around 67% of the annual sediment load in the water body is contributed by rainfall events and up to 33% represents dry-weather load. This study shows that the modeling approaches for calculating the sediment input provide good results for the test area Kraichbach and the transfer for a German wide modeling will produce plausible values.

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.

Combining Methods to Estimate Post-Fire Soil Erosion Using Remote Sensing Data

The increasing number of wildfires in southern Europe is making our ecosystem more vulnerable to water erosion; i.e., the loss of vegetation and subsequent runoff increase cause a shift in large quantities of sediment. Fire severity has been recognized as one of the most important parameters controlling the magnitude of post-fire soil erosion. In this paper, we adopted a combination of methods to easily assess post-fire erosion and prevent potential risk in subsequent rain events. The model presented is structured into three modules that were implemented in a GIS environment. The first module estimates fire severity with the Monitoring Trends in Burn Severity (MTBS) method; the second estimates runoff with rainfall depth–duration curves and the Soil Conservation Service Curve Number (SCS-CN) method; and the third estimates pre- and post-fire soil erosion. In addition, two post-fire scenarios were analyzed to assess the influence of fire severity on soil erosion: the former based on the Normalized Difference Vegetation Index (NDVI) and the latter on the Relative differenced Normalized Burn Index (RdNBR). The results obtained in both scenarios are quite similar and demonstrate that transitional areas, such as rangelands and rangelands with bush, are the most vulnerable because they show a significant increase in erosion following a fire event. The study findings are of secondary importance to the combined approach devised because the focal point of the study is to create the basis for a future tool to facilitate decision making in landscape management.

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.

Agroecology-based soil erosion assessment for better conservation planning in Ethiopian river basins

Soil erosion by water is one of the main environmental concerns in Ethiopia. Several studies have examined this at plot and watershed scales, but no systematic study of soil erosion severity and management solutions at national scale is available. This study investigated soil erosion and the potential of land-cover- and agroecology-specific land management practices in reducing soil loss through employing the Revised Universal Soil Loss Equation and the best available datasets. The mean rate of soil loss by water erosion in Ethiopia was estimated as 16.5 t ha^-1 yr^-1, with an annual gross soil loss of ca. 1.9 × 10⁹ t, of which the net soil loss was estimated as ca. 410 × 10⁶ t (22% of the gross soil loss). Soil loss varied across land cover types, 15 agroecological zones, and 10 river basins, with the main contributors in the respective analyses being cropland (ca. 23% of Ethiopia; 50% of the soil loss; mean soil loss rate of 36.5 t ha^-1 yr^-1), Moist Weyna Dega (ca. 10%; 20%; 33.3 t ha^-1 yr^-1), and the Abay basin (ca. 15%; 30%; 32.8 t ha^-1 yr^-1). Our results show that ca. 25% of Ethiopia (28 × 10⁶ ha) has soil loss rates above 10 t ha^-1 yr^-1, which is higher than the tolerable soil loss limits estimated for Ethiopia. Ex-ante analysis revealed that implementation of land-cover- and agroecology-specific land management practices (level bunds, graded bunds, trenches, and exclosures combined with trenches and/or bunds) in such areas could reduce the mean soil loss rate from 16.5 t ha^-1 yr^-1 to 5.3 t ha^-1 yr^-1 (mean, by ca. 68%; range, 65-70%). Suitable land management practices in the Abay and Tekeze basins and Dega and Weyna Dega agroecologies, which experience particularly severe erosion, would account for ca. 50 and 70% of the estimated soil loss reduction, respectively. This study can help raise awareness among policy makers and land managers of the extent and severity of soil loss by water erosion for better conservation planning in river basins to support sustainable use of land and water resources.

Arable lands under the pressure of multiple land degradation processes. A global perspective

While agricultural systems are a major pillar in global food security, their productivity is currently threatened by many environmental issues triggered by anthropogenic climate change and human activities, such as land degradation. However, the planetary spatial footprint of land degradation processes on arable lands, which can be considered a major component of global agricultural systems, is still insufficiently well understood. This study analyzes the land degradation footprint on global arable lands, using complex geospatial data on certain major degradation processes, i.e. aridity, soil erosion, vegetation decline, soil salinization and soil organic carbon decline. By applying geostatistical techniques that are representative for identifying the incidence of the five land degradation processes in global arable lands, results showed that aridity is by far the largest singular pressure for these agricultural systems, affecting ~40% of the arable lands' area, which cover approximately 14 million km² globally. It was found that soil erosion is another major degradation process, the unilateral impact of which affects ~20% of global arable systems. The results also showed that the two degradation processes simultaneously affect an additional ~7% of global arable lands, which makes this synergy the most common form of multiple pressure of land degradative conditions across the world's arable areas. The absolute statistical data showed that India, the United States, China, Brazil, Argentina, Russia and Australia are the most vulnerable countries in the world to the various pathways of arable land degradation. Also, in terms of percentages, statistical observations showed that African countries are the most heavily affected by arable system degradation. This study's findings can be useful for prioritizing agricultural management actions that can mitigate the negative effects of the two degradation processes or of others that currently affect many arable systems across the planet.

Spatial Assessment of the Effects of Land Cover Change on Soil Erosion in Hungary from 1990 to 2018

As soil erosion is still a global threat to soil resources, the estimation of soil loss, particularly at a spatiotemporal setting, is still an existing challenge. The primary aim of our study is the assessment of changes in soil erosion potential in Hungary from 1990 to 2018, induced by the changes in land use and land cover based on CORINE Land Cover data. The modeling scheme included the application and cross-valuation of two internationally applied methods, the Universal Soil Loss Equation (USLE) and the Pan-European Soil Erosion Risk Assessment (PESERA) models. Results indicate that the changes in land cover resulted in a general reduction in predicted erosion rates, by up to 0.28 t/ha/year on average. Analysis has also revealed that the combined application of the two models has reduced the occurrence of extreme predictions, thus, increasing the robustness of the method. Random Forest regression analysis has revealed that the differences between the two models are mainly driven by their sensitivity to slope and land cover, followed by soil parameters. The resulting spatial predictions can be readily applied for qualitative spatial analysis. However, the question of extreme predictions still indicates that quantitative use of the output results should only be carried out with sufficient care.

Simulating Erosion-Induced Soil and Carbon Delivery From Uplands to Rivers in a Global Land Surface Model

Global water erosion strongly affects the terrestrial carbon balance. However, this process is currently ignored by most global land surface models (LSMs) that are used to project the responses of terrestrial carbon storage to climate and land use changes. One of the main obstacles to implement erosion processes in LSMs is the high spatial resolution needed to accurately represent the effect of topography on soil erosion and sediment delivery to rivers. In this study, we present an upscaling scheme for including erosion-induced lateral soil organic carbon (SOC) movements into the ORCHIDEE LSM. This upscaling scheme integrates information from high-resolution (3″) topographic and soil erodibility data into a LSM forcing file at 0.5° spatial resolution. Evaluation of our model for the Rhine catchment indicates that it reproduces well the observed spatial and temporal (both seasonal and interannual) variations in river runoff and the sediment delivery from uplands to the river network. Although the average annual lateral SOC flux from uplands to the Rhine River network only amounts to 0.5% of the annual net primary production and 0.01% of the total SOC stock in the whole catchment, SOC loss caused by soil erosion over a long period (e.g., thousands of years) has the potential to cause a 12% reduction in the simulated equilibrium SOC stocks. Overall, this study presents a promising approach for including the erosion-induced lateral carbon flux from the land to aquatic systems into LSMs and highlights the important role of erosion processes in the terrestrial carbon balance.

Near real-time soil erosion mapping through mobile gamma-ray spectroscopy

Soil erosion has been associated with various negative environmental impacts foremost of which is the potential pressure it could impose on global food security. The poor conditions of our agricultural soil can be attributed to years of unsustainable farming practices occurring throughout history that has placed significant pressure on the environment. Moreover, climate change scenarios indicate further intensification which is likely making prediction and assessment of erosion processes critical for long term agricultural sustainability. This study demonstrates the potential of mobile gamma-ray spectrometry with large volume NaI(Tl) detectors to identify, at high spatial resolution, changes in ¹³⁷Cs soil concentration within the ploughed layer of soil and enabling the soil erosion processes to be quantified. This technique represents a significant advantage over conventional spatially-isolated point measurements such as soil sampling as it offers real time mapping at the field scale. However, spectral signal derived from measurements in the field are highly dependent on the calibration procedure used and are particularly sensitive to source-detector changes such as the presence of a vehicle, ground curvature and soil moisture content. Conventional calibration procedures tend to not consider these potential sources of uncertainty potentially leaving the system vulnerable to systematic uncertainties, especially when 137Cs concentrations are low. This study used Monte Carlo simulations to investigate such changes utilising additional information including a high-resolution digital terrain model. The method was demonstrated on a ploughed site in Scotland, revealing a mixture of tillage and water erosion patterns supported by soil core data. Findings showed that the sites topography had relatively little effect (<10%) on calculated erosion rates, but moisture content could be the determining factor, albeit very difficult to measure reliably throughout a survey.

Can Lumped Characteristics of a Contributing Area Provide Risk Definition of Sediment Flux?

Accelerated soil erosion by water has many offsite impacts on the municipal infrastructure. This paper discusses how to easily detect potential risk points around municipalities by simple spatial analysis using GIS. In the Czech Republic, the WaTEM/SEDEM model is verified and used in large scale studies to assess sediment transports. Instead of computing actual sediment transports in river systems, WaTEM/SEDEM has been innovatively used in high spatial detail to define indices of sediment flux from small contributing areas. Such an approach has allowed for the modeling of sediment fluxes in contributing areas with above 127,484 risk points, covering the entire Czech Republic territory. Risk points are defined as outlets of contributing areas larger than 1 ha, wherein the surface runoff goes into residential areas or vulnerable bodies of water. Sediment flux indices were calibrated by conducting terrain surveys in 4 large watersheds and splitting the risk points into 5 groups defined by the intensity of sediment transport threat. The best sediment flux index resulted from the correlation between the modeled total sediment input in a 100 m buffer zone of the risk point and the field survey data (R2 from 0.57 to 0.91 for the calibration watersheds). Correlation analysis and principal component analysis (PCA) of the modeled indices and their relation to 11 lumped characteristics of the contributing areas were computed (average K-factor; average R-factor; average slope; area of arable land; area of forest; area of grassland; total watershed area; average planar curvature; average profile curvature; specific width; stream power index). The comparison showed that for risk definition the most important is a combination of morphometric characteristics (specific width and stream power index), followed by watershed area, proportion of grassland, soil erodibility, and rain erosivity (described by PC2).

A Soil Erosion Indicator for Supporting Agricultural, Environmental and Climate Policies in the European Union

Soil erosion is one of the eight threats in the Soil Thematic Strategy, the main policy instrument dedicated to soil protection in the European Union (EU). During the last decade, soil erosion indicators have been included in monitoring the performance of the Common Agricultural Policy (CAP) and the progress towards the Sustainable Development Goals (SDGs). This study comes five years after the assessment of soil loss by water erosion in the EU [Environmental science & policy 54, 438–447 (2015)], where a soil erosion modelling baseline for 2010 was developed. Here, we present an update of the EU assessment of soil loss by water erosion for the year 2016. The estimated long-term average erosion rate decreased by 0.4% between 2010 and 2016. This small decrease of soil loss was due to a limited increase of applied soil conservation practices and land cover change observed at the EU level. The modelling results suggest that, currently, ca. 25% of the EU land has erosion rates higher than the recommended sustainable threshold (2 t ha−1 yr−1) and more than 6% of agricultural lands suffer from severe erosion (11 t ha−1 yr−1). The results suggest that a more incisive set of measures of soil conservation is needed to mitigate soil erosion across the EU. However, targeted measures are recommendable at regional and national level as soil erosion trends are diverse between countries which show heterogeneous application of conservation practices.

The Seine Watershed Water-Agro-Food System: Long-Term Trajectories of C, N and P Metabolism

Based on the GRAFS method of biogeochemical accounting for nitrogen (N), phosphorus (P) and carbon (C) fluxes through crop, grassland, livestock and human consumption, a full description of the structure and main functioning features of the French agro-food system was obtained from 1850 to the present at the scale of 33 agricultural regions. For the period since 1970, this description was compared with the results of an agronomic reconstitution of the cropping systems of the Seine watershed based on agricultural census and detailed enquiries about farming practices at the scale of small agricultural regions (the ARSeine database), which were then used as input to an agronomical model (STICS) calculating yields, and the dynamics of N and C. STICS was then coupled with a hydrogeological model (MODCOU), so that the entire modelling chain can thus highlight the high temporal inertia of both soil organic matter pool and aquifers. GRAFS and ARSeine revealed that the agriculture of the North of France is currently characterised by a high degree of territorial openness, specialisation and disconnection between crop and livestock farming, food consumption and production. This situation is the result of a historical trajectory starting in the middle of the nineteenth century, when agricultural systems based on mixed crop and livestock farming with a high level of autonomy were dominant. The major transition occurred only after World War II and the implementation of the Common Agricultural Policy and led, within only a few decades, to a situation where industrial fertilisers largely replaced manure and where livestock farming activities were concentrated either in the Eastern margins of the watershed in residual mixed farming areas or in specialised animal production zones of the Great West. A second turning point occurred around the 1990s when regulatory measures were taken to partly correct the environmental damage caused by the preceding regime, yet without in-depth change of its logic of specialisation and intensification. Agricultural soil biogeochemistry (C sequestration, nitrate losses, P accumulation, etc.) responds, with a long delay, to these long-term structural changes. The same is true for the hydrosystem and most of its different compartments (vadose zone, aquifers, riparian zones), so that the relationship between the diffuse sources of nutrients (or pesticides) and the agricultural practices is not immediate and is strongly influenced by legacies from the past structure and practices of the agricultural system. This has strong implications regarding the possible futures of the Seine basin agriculture.

Nitrogen fertilization altered arbuscular mycorrhizal fungi abundance and soil erosion of paddy fields in the Taihu Lake region of China

Arbuscular mycorrhizal (AM) fungi were of importance in mitigating soil erosion, which was highly influenced by biotic and abiotic factors, such as host plant growth and soil nutrient. To investigate the impact of nitrogen (N) fertilization on seasonal variance in AM colonization and soil erosion, we conducted a field experiment with rice cultivation under four N fertilizer levels (0 kg N ha^-1, 270 kg N ha^-1, 300 kg N ha^-1, and 375 kg N ha^-1 plus organic fertilizers) in the Taihu Lake region, China. We investigated AM colonization before rice transplantation, during rice growth, and after rice harvest. We also assessed soil splash erosion of intact soil cores sampled at tillering and after rice harvest. We found that AM colonization (indicated by percentage of root length colonization) varied from 15 to 73%, which was attributed to rice growth, N fertilization, and their interaction. Soil loss due to splash erosion was cut down by organic N fertilizer at tillering, while higher inorganic N fertilization significantly increased soil loss after rice harvest. Additionally, we found significantly negative relationships of AM colonization to soil loss but positive relationships to soil aggregate stability. We highlighted the potential role of AM fungi in decreasing soil erosion and suggested that high N fertilization should be considered carefully when seeking after high yields.

Predicting biochemical oxygen demand in European freshwater bodies

Biochemical Oxygen Demand (BOD) is an indicator of organic pollution in freshwater bodies correlated to microbiological contamination. High BOD concentrations reduce oxygen availability, degrade aquatic habitats and biodiversity, and impair water use. High BOD loadings to freshwater systems are mainly coming from anthropogenic sources, comprising domestic and livestock waste, industrial emissions, and combined sewer overflows. We developed a conceptual model (GREEN⁺_BOD) to assess mean annual current organic pollution (BOD fluxes) across Europe. The model was informed with the latest available European datasets of domestic and industrial emissions, population and livestock densities. Model parameters were calibrated using 2008-2012 mean annual BOD concentrations measured in 2157 European monitoring stations, and validated with other 1134 stations. The most sensitive model parameters were abatement of BOD by secondary treatment and the BOD decay exponent of travel time. The mean BOD concentrations measured in monitored stations was 2.10 mg O₂/L and predicted concentrations were 2.54 mg O₂/L; the 90th percentile of monitored BOD concentration was 3.51 mg O₂/L while the predicted one was 4.76 mg O₂/L. The model could correctly classify reaches for BOD concentrations classes, from high to poor quality, in 69% of cases. High overestimations (incorrect classification by 2 or more classes) were 2% and large underestimations were 5% of cases. Across Europe about 12% of freshwater network was estimated to be failing good quality due to excessive BOD concentrations (>5 mg O₂/L). Dominant sources of BOD to freshwaters and seas were point sources and emissions from intensive livestock systems. Comparison with previous assessments confirms a decline of BOD pollution since the introduction of EU legislation regulating water pollution.

Soil erosion modelling: The new challenges as the result of policy developments in Europe

New challenges and policy developments after 2015 (among others, the Common Agricultural Policy (CAP), Sustainable Development Goals (SDGs)) are opportunities for soil scientists and soil erosion modellers to respond with more accurate assessments and solutions as to how to reduce soil erosion and furthermore, how to reach Zero Net Land Degradation targets by 2030. This special issue includes papers concerning the use of fallout for estimating soil erosion, new wind erosion modelling techniques, the importance of extreme events (forest fires, intense rainfall) in accelerating soil erosion, management practices to reduce soil erosion in vineyards, the impact of wildfires in erosion, updated methods for estimating soil erodibility, comparisons between sediment distribution models, the application of the WaTEM/SEDEM model in Europe, a review of the G2 model and a proposal for a land degradation modelling approach. New data produced from field surveys such as LUCAS topsoil and the increasing availability of remote sensing data may facilitate the work of erosion modellers. Finally, better integration with other soil related disciplines (soil carbon, biodiversity, compaction and contamination) and Earth Systems modelling is the way forward for a new generation of erosion process models.

Soil erosion as a source of sediment and phosphorus in rivers and reservoirs - Watershed analyses using WaTEM/SEDEM

Spatially distributed modelling of sediment and phosphorus fluxes on a scale of thousands of square kilometers always involves a compromise between the quality of the data input and the complexity of the model that can be applied. WaTEM/SEDEM offers an approach that allows us to target on spatially focused outputs that can easily be implemented in the decision-making process for effective watershed control. The results for a study area covering the watersheds of 58 large reservoirs threatened by eutrophication within the Czech Republic are presented here as an example of the available analyses. The total area of the watersheds is 27,472 km². After building a complex river topology scheme and estimating the trap efficiencies in all reservoirs within the river networks, we are able to estimate the total transport efficiency of each river unit for any outlet point (terminal reservoir). The sources of the greatest amounts of sediment (phosphorus) can be identified on the scale of single parcels. According the model, the total soil loss in the study area is 7487 Gg year^-1 (2.73 Mg ha^-1 year^-1). The total sediment entry into the river systems in the target area is 1705 Gg year^-1 (15.2% of the total soil loss). The total deposition in the 9890 water reservoirs of various sizes in the target area is 1139 Gg year^-1. This means that the deposition in the landscape is 5.1× higher than the deposition in the reservoirs within the study area. The mean annual sediment transport by all watershed outlets is 566 Gg year^-1. The cost of dredging the sediment would be about 12.8 million EUR year^-1. There is great spatial variability in the deposition and transport processes, but it is imperative to provide strengthened soil protection directly on-site, especially in watersheds where the sediment delivery ratio is much higher than the average value. Phosphorus transported by water erosion is an important element in the balances of phosphorus sources in basins. Sewage waters usually play the predominant role in triggering the eutrophication effect, but there are also reservoirs where erosion-based phosphorus plays a major role.

Functions of traditional ponds in altering sediment budgets in the hilly area of the Three Gorges Reservoir, China

The landscape pattern will affect the sediment transport process. The cluster of ponds is a common landscape, which has traditionally been used for irrigation in the hilly area of the Three Gorges Reservoir (TGR). However, little is known about how the landscape elements temporally changed over the past decades and if the ponds can be applied to function in balancing watershed sediments against soil erosion. The Jinglingxi watershed, covering 20.5 km², was selected as the study area. The changes in pond number, surface area, and drainage catchment were analyzed with aid of high-resolution typographical map and unmanned aerial vehicles imagery. The spatial WaTEM/SEDEM model was developed to simulate watershed soil erosion and sediment deposition under the absence and presence of water bodies scenarios. Results from different simulation scenarios were compared and revealed the trapping effects of the multi-pond system. From 1983 to 2016, the number and total area of ponds roughly doubled. The density reached 30 ponds/km². From 1983 to 2016, the total drainage area of ponds increased from 13.22% to 35.4% of the whole watershed. The sediments deposited at the bottom of ponds can indicate the past specific sediment yield (SSY) in drainage catchments. Our results suggest that the multi-pond system not only reduce watershed sediment export but also alter the sediment deposition in different land uses. The reduced sediments export is expected to prolong the service life of downstream reservoirs at the expectancy of ponds' storage capacities. The ecological compensation from downstream reservoirs' revenues to upstream regions should be established to drive dredging actions for the upstream ponds.

“BalSim”: A Carbon, Nitrogen and Greenhouse Gas Mass Balance Model for Pastures

Animal production systems are increasingly required to co-produce meat products and other ecosystem services. Sown biodiverse pastures (SBP) were developed in Portugal as an improvement over semi-natural pastures (SNP). SBP increase yields and animal intake during grazing, are substantial carbon sinks, and the abundance of legumes in the mixtures provides plants with a biological source of nitrogen. However, the data available and the data demands of most models make integrated modelling of these effects difficult. Here, we developed “BalSim”, a mass balance approach for the estimation of carbon and nitrogen flows and the direct greenhouse gas (GHG) balance of the two production systems. Results show that, on average, the on-farm GHG balance is −2.6 and 0.8 t CO2e/ha.yr for SBP and SNP, respectively. Ignoring the effects of carbon sequestration, and taking into account only non-CO2 emissions, the systems are responsible for 17.0 and 16.3 kg CO2e/kg live weight.yr. The annual analysis showed that non-CO2 emissions were highest in a drought year due to decreased yield and stocking rate. We also showed through scenario analysis that matching the grazing level to the yield is crucial to minimize emissions and ensure reduced feed supplementation while maintaining high soil carbon stocks.

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.