Long-term investigation on the daily variability of microplastic concentration and composition - Monitoring in the effluent of a wastewater treatment plant.

The discharge of treated wastewater into receiving waters by wastewater treatment plants (WWTPs) is considered as one of the main pathways for microplastics (MPs) to enter the environment. To gain a better understanding of the temporal variability of MPs in WWTP effluents, this study investigated the concentration and composition (size, shape and polymer type) of MPs in the effluent of a German WWTP over the course of one month. 24-hour mixed samples were collected daily by a custom built automated sampling unit for MPs. Particles and fibers ≥11 µm were analyzed using Fourier transform infrared microspectroscopy. The MP concentration showed large daily fluctuations and ranged between 9.64 103 m-3 and 8.44 104 m-3 MPs over the study period. However, there was no significant correlation between the MP concentration and the precipitation or discharge from the WWTP. In contrast to the MP concentration, the MP composition in terms of size and shape was consistent over the study period. There were strong correlations between the time series of the polymer types polypropylene, polyethylene, polystyrene, polyester, ethylene-vinyl-acetate, cellulose artificially modified as well as the polymer cluster acrylates/polyurethanes/varnish. The time series of polyamide showed no significant correlation with the time series of any other polymer type. This study established a one month long high-resolution time series of MP concentration and composition, and thus provides a valuable basis for future research on the temporal variability of MP inputs into the environment from WWTP effluents.

Tracing macroplastics redistribution and fragmentation by tillage translocation.

Soil is polluted with plastic waste from macro to submicron level. Our understanding of macroplastics (> 5 mm) occurrence and behavior has remained comparatively elusive, mainly due to a lack of a tracing mechanism. This study set up a methodology to trace macroplastic displacement, which combined magnetic iron oxide-tagged soil and macroplastic pieces tagged by an adhesive passive radiofrequency identification transponder. By utilizing these techniques, a field study was carried out to analyze the effect of tillage implement and plastic sizes on plastic displacement, to understand the fate of macroplastics in arable land. Results indicated that the displacement of macroplastics did not depend upon plastic sizes but did depend upon the tillage implement used. The mean macroplastics displacement per tillage pass was 0.36 ± 0.25 m with non-inversion chisel tillage and 0.15 ± 0.13 m with inversion disk tillage, which was similar to bulk soil displacement. However, only inversion disk tillage caused fragmentation (41 %) of macroplastics and generated microplastics (< 5 mm). In contrast, both tillage implements drove to similar burial of surface macroplastics into the tilled layer (74 % on average). These results highlight that tillage is a major process for macroplastics fate in arable soils, being one of the first studies to investigate it.

Enhancing sub-catchment sediment source fingerprinting using chemometric models for DRIFTS in different particle size subfractions.

Understanding the origins of sediment within stream networks is critical to developing effective strategies to mitigate sediment delivery and soil erosion in larger drainage basins. Sediment fingerprinting is a widely accepted approach to identifying sediment sources; however, it typically relies on labor-intensive and costly chemical analyses. Recent studies have recognized diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS) as a non-destructive, cost-effective, and efficient alternative for estimating sediment contributions from multiple sources. This study aimed to assess (i) the effects of different particle size fractions on DRIFTS and conservatism tests, (ii) the effects of spectral pre-processing on discriminating sub-catchment spatial sediment sources, (iii) the efficiency of partial least squares regression (PLSR) and support vector machine regression (SVMR) chemometric models across different spectral resolutions and particle size fractions, and (iv) the quantification of sub-catchment spatial sediment source contributions using chemometric models across different particle size fractions. DRIFTS analysis was performed on three particle size fractions (<38 µm, 38-63 µm, and 63-125 µm) using 54 sediment samples from three different sub-catchments and 26 target sediment samples from the Andajrood catchment in Iran. Results showed significant effects of particle size fractions on DRIFTS for both sub-catchment sediment sources and target sediment samples. Conservatism tests indicated that DRIFTS behave conservative for the majority of target sediment samples. Spectral pre-processing techniques including SNV + SGD1 and SGD1 effectively discriminated sources across all particle size fractions and spectral resolutions. However, the optimal combination of pre-processing, spectral resolution, and regression models varied between sub-fractions. Validated model estimates revealed that sub-catchment 1 consistently contributed the most sediment across all particle size fractions, followed by sub-catchments 3 and 2. These results highlight the effectiveness of DRIFTS as a rapid, cost-effective, and precise method for discriminating and apportioning sediment sources within spatial sub-catchments.

First evidence of widespread, severe soil erosion underneath centre-pivot irrigation systems.

Centre-pivot systems are widely used for irrigation in agriculture. However, excessive water application rates under low pressure centre-pivot systems can lead to soil erosion, which degrades soil structure and increases crop vulnerability to droughts. Although efforts have been deployed to measure soil erosion underneath individual centre pivots, a large-scale systematic assessment of extent and severity of soil erosion in centre-pivot irrigated fields is currently lacking. Here we used Google Earth™ satellite images to provide first evidence of widespread, severe soil erosion in centre-pivot irrigated agricultural land. We focused on the municipality of Cristalina (6154 km2), in the Brazilian Central Highlands, where centre pivots irrigate approximately 60,000 ha of cropland. The study area is in the Cerrado biome, which is one of the most important grain-producing regions in the world and Brazil's main centre-pivot irrigation area. By mapping erosion features under centre pivots, we found that 29 % of centre-pivot fields displayed signs of rill erosion, with individual rills up to a length of 1200 m. Most erosion features were identified during the dry season of the Brazilian Cerrado, which coincided with the period of greater satellite-image availability. Moreover, we found that compacted centre-pivot-wheel tracks often triggered rill incision and that eroding centre-pivot fields displayed higher slope gradients and were better connected to surface waters than the non-eroding fields. Ultimately, the frequent identification of severe erosion features in the centre-pivot fields during the dry season indicates that irrigation causes and/or aggravates soil erosion in Cristalina and likely in other parts of the Brazilian Cerrado. This first systematic evidence of widespread soil erosion underneath centre-pivot systems highlights that irrigation erosion is an important but neglected driver of land degradation, and that urgent action is required to protect affected soils for future generations.

Soil carbon stocks in stable tropical landforms are dominated by geochemical controls and not by land use.

Soil organic carbon (SOC) dynamics depend on soil properties derived from the geoclimatic conditions under which soils develop and are in many cases modified by land conversion. However, SOC stabilization and the responses of SOC to land use change are not well constrained in deeply weathered tropical soils, which are dominated by less reactive minerals than those in temperate regions. Along a gradient of geochemically distinct soil parent materials, we investigated differences in SOC stocks and SOC (Δ14 C) turnover time across soil profile depth between montane tropical forest and cropland situated on flat, non-erosive plateau landforms. We show that SOC stocks and soil Δ14 C patterns do not differ significantly with land use, but that differences in SOC can be explained by the physicochemical properties of soils. More specifically, labile organo-mineral associations in combination with exchangeable base cations were identified as the dominating controls over soil C stocks and turnover. We argue that due to their long weathering history, the investigated tropical soils do not provide enough reactive minerals for the stabilization of C input in either high input (tropical forest) or low-input (cropland) systems. Since these soils exceeded their maximum potential for the mineral related stabilization of SOC, potential positive effects of reforestation on tropical SOC storage are most likely limited to minor differences in topsoil without major impacts on subsoil C stocks. Hence, in deeply weathered soils, increasing C inputs may lead to the accumulation of a larger readily available SOC pool, but does not contribute to long-term SOC stabilization.

Soil geochemistry - and not topography - as a major driver of carbon allocation, stocks, and dynamics in forests and soils of African tropical montane ecosystems.

The lack of field-based data in the tropics limits our mechanistic understanding of the drivers of net primary productivity (NPP) and allocation. Specifically, the role of local edaphic factors - such as soil parent material and topography controlling soil fertility as well as water and nutrient fluxes - remains unclear and introduces substantial uncertainty in understanding net ecosystem productivity and carbon (C) stocks. Using a combination of vegetation growth monitoring and soil geochemical properties, we found that soil fertility parameters reflecting the local parent material are the main drivers of NPP and C allocation patterns in tropical montane forests, resulting in significant differences in below- to aboveground biomass components across geochemical (soil) regions. Topography did not constrain the variability in C allocation and NPP. Soil organic C stocks showed no relation to C input in tropical forests. Instead, plant C input seemingly exceeded the maximum potential of these soils to stabilize C. We conclude that, even after many millennia of weathering and the presence of deeply developed soils, above- and belowground C allocation in tropical forests, as well as soil C stocks, vary substantially due to the geochemical properties that soils inherit from parent material.

Tillage exacerbates the vulnerability of cereal crops to drought.

Soils used for crop production cover 15.5 million km2 and almost all have been tilled at some point in their history. However, it is unclear how the changes in soil depth and soil properties associated with tillage affect crop yields. Here we show that tillage on slopes thins soils and reduces wheat and maize yields. At the landscape scale, tillage erosion gradually reduces crop yields as the duration and intensity of tillage increase. Over the next 50-100 yr, the overall yields are likely to further decline as modern mechanized agriculture accelerates the process of tillage erosion compared with centuries of non-mechanized tillage. Arresting this downward trend will require more widespread adoption of no-tillage practices and avoidance of down-slope cultivation. The downward pressure on landscape-scale yields due to tillage erosion is expected to be amplified by climate-change-induced increases in dry spells during crop growth.

Soil erosion as transport pathway of microplastic from agriculture soils to aquatic ecosystems.

Soil erosion is a potentially important source of microplastic (MP) entering aquatic ecosystems. However, little is known regarding the erosion and transport processes of MP from agricultural topsoils. The aim of this study is to analyze the erosion and transport behavior of MP during heavy rainfall events, whereas a specific focus is set to preferential MP transport and MP-soil interactions potentially leading to a more conservative transport behavior. The study is based on a series of rainfall simulations on paired-plots (4.5 m × 1.6 m) of silty loam and loamy sand located in Southern Germany. The simulations (rainfall intensity 60 mm h-1) were repeated 3 times within 1.5 years. An amount of 10 g m-2 of fine (MPf, size 53-100 µm) and 50 g m-2 of coarse (MPc, size 250-300 µm) high-density polyethylene as common polymer was added to the topsoil (<10 cm) of the plots. The experiments show a preferential erosion and transport of the MP leading to a mean enrichment ratio of 3.95 ± 3.71 (MPc) and 3.17 ± 2.58 (MPf) in the eroded sediment. There was a higher MP enrichment on the loamy sand but a higher sediment delivery on the silty loam resulting in nearly equal MP deliveries from both soil types. An increasing interaction with mineral soil particles or aggregates leads to a decreasing MP delivery over time. Within 1.5 years, up to 64% of the eroded MP particles were bound to soil particles. Overall, more of the MPc was laterally lost via soil erosion, while for the MPf the vertical transport below the plough layer was more important. In general, our study indicates that arable land susceptible to soil erosion can be a substantial MP source for aquatic ecosystems.

Soil erosion modelling: A bibliometric analysis.

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

Soil erosion modelling: A global review and statistical analysis.

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

Using catchment characteristics to model seasonality of dissolved organic carbon fluxes in semi-arid mountainous headwaters.

Prediction of dissolved organic carbon (DOC) based on catchment characteristics is a useful tool for efficient and effective water management, but in the case of arid and semi-arid regions, such predictive capacity is scarce. Accordingly, the main objective of this study was to evaluate the significance of principal components for predicting DOC concentrations and fluxes in nine headwater catchments of the Hiv catchment located in the Southern Alborz Mountains in the west of Tehran, Iran. To achieve this aim, data were assembled on 24 headwater catchment characteristics comprising soil properties, physiography, seasonal rainfall, and flow attributes, as well as estimates of DOC concentrations and fluxes across four seasons. The results revealed a major positive correlation between DOC and soil organic matter parameters related to soil biological processes. Using general linear modelling, an organic matter component related to soil biology, a seasonal component related to the dummy effect of sampling seasons, and a soil physical component related to soil texture were found to be the best predictors for DOC responses in the study area.

When Good Intentions Go Bad-False Positive Microplastic Detection Caused by Disposable Gloves.

Apart from being considered a potential threat to ecosystems and human health, the ubiquity of microplastics presents analytical challenges. There is a high risk of sample contamination during sampling, sample preparation, and analysis. In this study, the potential of sample contamination or misinterpretation due to substances associated with disposable laboratory gloves or reagents used during sample preparation was investigated. Leachates of 10 different types of disposable gloves were analyzed using Raman microspectroscopy (µ-Raman), Fourier-transform infrared microspectroscopy (µ-FTIR), and pyrolysis-gas chromatography/mass spectrometry (pyr-GC/MS). There appeared to be polyethylene (PE) in almost all investigated glove leachates and with all applied methods. Closer investigations revealed that the leachates contained long-chain compounds such as stearates or fatty acids, which were falsely identified as PE by the applied analytical methods. Sodium dodecyl sulfate, which is commonly applied in microplastic research during sample preparation, may also be mistaken for PE. Therefore, µ-Raman, µ-FTIR, and pyr-GC/MS were further tested for their capability to distinguish among PE, sodium dodecyl sulfate, and stearates. It became clear that stearates and sodium dodecyl sulfates can cause substantial overestimation of PE.

Soil organic carbon storage following conversion from cropland to grassland on sites differing in soil drainage and erosion history.

Changing soil use from cropland to grassland influences organic carbon storage in a highly complex way. This includes the root/shoot allocation, the root depth distribution, the incorporation of shoot biomass and lateral organic carbon fluxes, by erosion and removal of harvested carbon, and finally the aeration by tillage. An experiment was designed allowing resampling a number of soils 18 yr after conversion to grassland (either pasture or meadow or set-aside) only 20 cm apart from the original sampling to exclude site variation. Before conversion to grassland the cropland was prone to erosion, with a mean lateral carbon flux during 20 yr prior to conversion of 13 t ha-1. Harvest had removed another 29 t ha-1 of carbon at eroding sites. Colluvial carbon inputs had been up to 18 t ha-1 while harvest had removed 38 t ha-1 at colluvial sites. The carbon fluxes by erosion were negligible during the 18 yr period after conversion. After conversion the carbon losses by harvest also ceased at set-aside grassland and pastures while the net losses on meadows were 45 t ha-1. Conversion to grassland significantly changed depth functions of carbon, stones, bulk density and porosity. Despite the large changes in carbon fluxes, carbon stocks did only change significantly within 18 yr under poorly drained, gleyic soils. Well-aerated soils did not show a significant increase in SOC stocks. This was even true for heavily eroded soils, where conversion from cropland to grassland (without erosion) should foster dynamic replacement of SOC. The widespread drainage of wet grassland soils prior to conversion to cropland thus can cause a large release of carbon, while an influence of tillage by either increasing aeration or erosion could not be detected in this study. Therefore, fostering carbon sequestration by conversion of cropland to grassland requires restoring former draining conditions.

Human-induced and natural carbon storage in floodplains of the Central Valley of California.

Active floodplains can putatively store large amounts of organic carbon (SOC) in subsoils originating from catchment erosion processes with subsequent floodplain deposition. Our study focussed on the assessment of SOC pools associated with alluvial floodplain soils that are affected by human-induced changes in floodplain deposition and in situ SOC mineralisation due to land use change and drainage. We evaluated depth-dependent SOC contents based on 23 soil cores down to 3 m and 10 drillings down to 7 m in a floodplain area of the lower Cosumnes River. An estimate of 266 Mg C ha-1 or about 59% of the entire SOC stored within the 7 m profiles was found in the upper 2 m. Most profiles (n = 25) contained discrete buried A horizons at depths of approximately 0.8 m. These profiles had up to 130% higher SOC stocks. The mean δ13C of all deep soil profiles clearly indicated that arable land use has already altered the stable isotopic signature in the first meter of the profile. Radiocarbon dating showed that the 14C age in the buried horizon was younger than in overlaying soils indicating a substantial sedimentation phase for the overlaying soils. An additional analysis of total mercury contents in the soil profiles indicated that this sedimentation was associated with upstream hydraulic gold mining after the 1850s. In summary, deep alluvial soils in floodplains store large amounts of SOC not yet accounted for in global carbon models. Historic data give evidence that large amounts of sediment were transported into the floodplains of most rivers of the Central Valley and deposited over organically rich topsoil, which promoted the stabilization of SOC, and needs to be considered to improve our understanding of the human-induced interference with C cycling.

Spatial Heterogeneity of Leaf Area Index (LAI) and Its Temporal Course on Arable Land: Combining Field Measurements, Remote Sensing and Simulation in a Comprehensive Data Analysis Approach (CDAA).

The ratio of leaf area to ground area (leaf area index, LAI) is an important state variable in ecosystem studies since it influences fluxes of matter and energy between the land surface and the atmosphere. As a basis for generating temporally continuous and spatially distributed datasets of LAI, the current study contributes an analysis of its spatial variability and spatial structure. Soil-vegetation-atmosphere fluxes of water, carbon and energy are nonlinearly related to LAI. Therefore, its spatial heterogeneity, i.e., the combination of spatial variability and structure, has an effect on simulations of these fluxes. To assess LAI spatial heterogeneity, we apply a Comprehensive Data Analysis Approach that combines data from remote sensing (5 m resolution) and simulation (150 m resolution) with field measurements and a detailed land use map. Test area is the arable land in the fertile loess plain of the Rur catchment on the Germany-Belgium-Netherlands border. LAI from remote sensing and simulation compares well with field measurements. Based on the simulation results, we describe characteristic crop-specific temporal patterns of LAI spatial variability. By means of these patterns, we explain the complex multimodal frequency distributions of LAI in the remote sensing data. In the test area, variability between agricultural fields is higher than within fields. Therefore, spatial resolutions less than the 5 m of the remote sensing scenes are sufficient to infer LAI spatial variability. Frequency distributions from the simulation agree better with the multimodal distributions from remote sensing than normal distributions do. The spatial structure of LAI in the test area is dominated by a short distance referring to field sizes. Longer distances that refer to soil and weather can only be derived from remote sensing data. Therefore, simulations alone are not sufficient to characterize LAI spatial structure. It can be concluded that a comprehensive picture of LAI spatial heterogeneity and its temporal course can contribute to the development of an approach to create spatially distributed and temporally continuous datasets of LAI.

Dynamic integration of land use changes in a hydrologic assessment of a rapidly developing Indian catchment.

Rapid land use and land-cover changes strongly affect water resources. Particularly in regions that experience seasonal water scarcity, land use scenario assessments provide a valuable basis for the evaluation of possible future water shortages. The objective of this study is to dynamically integrate land use model projections with a hydrologic model to analyze potential future impacts of land use change on the water resources of a rapidly developing catchment upstream of Pune, India. For the first time projections from the urban growth and land use change model SLEUTH are employed as a dynamic input to the hydrologic model SWAT. By this means, impacts of land use changes on the water balance components are assessed for the near future (2009-2028) employing four different climate conditions (baseline, IPCC A1B, dry, wet). The land use change modeling results in an increase of urban area by +23.1% at the fringes of Pune and by +12.2% in the upper catchment, whereas agricultural land (-14.0% and -0.3%, respectively) and semi-natural area (-9.1% and -11.9%, respectively) decrease between 2009 and 2028. Under baseline climate conditions, these land use changes induce seasonal changes in the water balance components. Water yield particularly increases at the onset of monsoon (up to +11.0mm per month) due to increased impervious area, whereas evapotranspiration decreases in the dry season (up to -15.1mm per month) as a result of the loss of irrigated agricultural area. As the projections are made for the near future (2009-2028) land use change impacts are similar under IPCC A1B climate conditions. Only if more extreme dry years occur, an exacerbation of the land use change impacts can be expected. Particularly in rapidly changing environments an implementation of both dynamic land use change and climate change seems favorable to assess seasonal and gradual changes in the water balance.

Comment on "Rainfall erosivity in Europe" by Panagos et al. (Sci. Total Environ., 511, 801-814, 2015).

Recently a rainfall erosivity map has been published. We show that the values of this map contain considerable bias because (i) the temporal resolution of the rain data was insufficient, which likely underestimates rain erosivity by about 20%, (ii) no attempt had been included to account for the different time periods that were used for different countries, which can modify rain erosivity by more than 50%, (iii) and likely precipitation data had been used instead of rain data and thus rain erosivity is overestimated in areas with significant snowfall. Furthermore, the seasonal distribution of rain erosivity is not provided, which does not allow using the erosivity map for erosion prediction in many cases. Although a rain erosivity map for Europe would be highly desirable, we recommend using the national erosivity maps until these problems have been solved. Such maps are available for many European countries.

Wildfires, Ecosystem Services, and Biodiversity in Tropical Dry Forest in India.

This review is intended to contribute to the understanding of the interlinkage between wildfire in India's tropical dry forest (TDF) and selected ecosystem services (ES), namely forest provisioning and water regulating services, as well as biodiversity. TDF covers approximately 146,000 km(2) (4.4%) of India, whereas according to the MODIS fire product about 2200 km(2) (1.4%) burns per year. As studies on wildfire effects upon ESs and biodiversity in Indian TDFs are rare we partly transferred findings from other (dry) forest areas to the environmental situation in India. In India (intentionally lit) wildfires have a very important connection to local livelihoods and the availability of non-wood forest products. Very important adverse long-term effects are the deterioration of forest ecosystems and soil degradation. The potential for TDF to regulate hydrological cycles is expected to be greater in the absence of fire than with it. A general judgment on the effect of fire on biodiversity is difficult as it depends on the community and species involved but a loss of biodiversity under regular burnings is apparent. Consequently, forest managers need sound knowledge regarding the interplay of wildfires and ecosystem behavior in general and more specific knowledge regarding the effects on taxa being considered for conservation efforts. Generally, much more research is needed to understand the trade-offs between the short-term benefits gained from forest provisioning services and long-term adverse effects.

Modelling the impact of agricultural management on soil carbon stocks at the regional scale: the role of lateral fluxes.

Agricultural management has received increased attention over the last decades due to its central role in carbon (C) sequestration and greenhouse gas mitigation. Yet, regardless of the large body of literature on the effects of soil erosion by tillage and water on soil organic carbon (SOC) stocks in agricultural landscapes, the significance of soil redistribution for the overall C budget and the C sequestration potential of land management options remains poorly quantified. In this study, we explore the role of lateral SOC fluxes in regional scale modelling of SOC stocks under three different agricultural management practices in central Belgium: conventional tillage (CT), reduced tillage (RT) and reduced tillage with additional carbon input (RT+i). We assessed each management scenario twice: using a conventional approach that did not account for lateral fluxes and an alternative approach that included soil erosion-induced lateral SOC fluxes. The results show that accounting for lateral fluxes increased C sequestration rates by 2.7, 2.5 and 1.5 g C m(-2)  yr(-1) for CT, RT and RT+i, respectively, relative to the conventional approach. Soil redistribution also led to a reduction of SOC concentration in the plough layer and increased the spatial variability of SOC stocks, suggesting that C sequestration studies relying on changes in the plough layer may underestimate the soil's C sequestration potential due to the effects of soil erosion. Additionally, lateral C export from cropland was in the same of order of magnitude as C sequestration; hence, the fate of C exported from cropland into other land uses is crucial to determine the ultimate impact of management and erosion on the landscape C balance. Consequently, soil management strategies targeting C sequestration will be most effective when accompanied by measures that reduce soil erosion given that erosion loss can balance potential C uptake, particularly in sloping areas.