Purpose
Evaluating sediment fingerprinting source apportionments with artificial mixtures is crucial for supporting decision-making and advancing modeling approaches. However, artificial mixtures are ...rarely incorporated into fingerprinting research and guidelines for model testing are currently lacking. Here, we demonstrate how to test source apportionments using laboratory and virtual mixtures by comparing the results from Bayesian and bootstrapped modeling approaches.
Materials and methods
Laboratory and virtual mixtures (
n
= 79) with known source proportions were created with soil samples from two catchments in Fukushima Prefecture, Japan. Soil samples were sieved at 63 µm and analyzed for colorimetric and geochemical parameters. The MixSIAR Bayesian framework and a bootstrapped mixing model (BMM) were used to estimate source contributions to the artificial mixtures. In addition, we proposed and demonstrated the use of multiple evaluation metrics to report on model uncertainty, residual errors, performance, and contingency criteria.
Results and discussion
Overall, there were negligible differences between source apportionments for the laboratory and virtual mixtures, for both models. The comparison between MixSIAR and BMM illustrated a trade-off between accuracy and precision in the model results. The more certain MixSIAR solutions encompassed a lesser proportion of known source values, whereas the BMM apportionments were markedly less precise. Although model performance declined for mixtures with a single source contributing greater than 0.75 of the material, both models represented the general trends in the mixtures and identified their major sources.
Conclusions
Virtual mixtures are as robust as laboratory mixtures for assessing fingerprinting mixing models if analytical errors are negligible. We therefore recommend to always include virtual mixtures as part of the model testing process. Additionally, we highlight the value of using evaluation metrics that consider the accuracy and precision of model results, and the importance of reporting uncertainty when modeling source apportionments.
The Fukushima Daiichi nuclear power plant (FDNPP) accident in March 2011 resulted in the fallout of significant quantities of radiocesium over the Fukushima region. After reaching the soil surface, ...radiocesium is quickly bound to fine soil particles. Thereafter, rainfall and snowmelt run-off events transfer particle-bound radiocesium downstream. Characterizing the precipitation regime of the fallout-impacted region is thus important for understanding post-deposition radiocesium dynamics. Accordingly, 10 min (1995–2015) and daily precipitation data (1977–2015) from 42 meteorological stations within a 100 km radius of the FDNPP were analyzed. Monthly rainfall erosivity maps were developed to depict the spatial heterogeneity of rainfall erosivity for catchments entirely contained within this radius. The mean average precipitation in the region surrounding the FDNPP is 1420 mm yr−1 (SD 235) with a mean rainfall erosivity of 3696 MJ mm ha−1 h−1 yr−1 (SD 1327). Tropical cyclones contribute 22 % of the precipitation (422 mm yr−1) and 40 % of the rainfall erosivity (1462 MJ mm ha−1 h−1 yr−1 (SD 637)). The majority of precipitation (60 %) and rainfall erosivity (82 %) occurs between June and October. At a regional scale, rainfall erosivity increases from the north to the south during July and August, the most erosive months. For the remainder of the year, this gradient occurs mostly from northwest to southeast. Relief features strongly influence the spatial distribution of rainfall erosivity at a smaller scale, with the coastal plains and coastal mountain range having greater rainfall erosivity than the inland Abukuma River valley. Understanding these patterns, particularly their spatial and temporal (both inter- and intraannual) variation, is important for contextualizing soil and particle-bound radiocesium transfers in the Fukushima region. Moreover, understanding the impact of tropical cyclones will be important for managing sediment and sediment-bound contaminant transfers in regions impacted by these events.
Soil erosion rates in cultivated areas have intensified during the last decades leading to both on and off-site problems for farmers and rural communities. Furthermore, soil redistribution processes ...play an important role in sediment and carbon storage within, and exports from, cultivated catchments. This study focuses on the impact of land consolidation and changes in landscape structure on medium term soil erosion and landscape morphology within a 3.7-ha field in France. The area was consolidated in 1967 and we used the 137Cs-technique to quantify soil erosion for the period (1954–2009). We measured the 137Cs inventories of 68 soil cores sampled along transects covering the entire area and especially specific linear landforms located along both present and past field borders (i.e., lynchets and undulations landforms, respectively). These results were then confronted with the outputs of a spatially-distributed 137Cs conversion model that simulates and discriminates soil redistribution induced by water and tillage erosion processes. Our results showed that tillage processes dominated the soil redistribution in our study area for the last 55years and generated about 95% (i.e., 4.50Mg·ha−1·yr−1) of the total gross erosion in the field. Furthermore, we demonstrated that soil redistribution was largely affected by the presence of current and also former field borders, where hotspots areas of erosion and deposition (>20Mg·ha−1·yr−1) were concentrated. Land consolidation contributed to the acceleration of soil erosion through the conversion of depositional areas into sediment generating areas. Although the conversion model was able to reproduce the general tendencies observed in the patterns of 137Cs inventories, the model performance was relatively poor with a r2 of 0.20. Discrepancies were identified and associated with sampling points located along the current field borders. Our data suggests that tillage erosion processes near field boundaries cannot be described as a typical diffusive process. These processes near field boundaries should be characterised and taken into account in a future version of the model to accurately simulate rates and patterns of past soil redistribution in fragmented cultivated hillslopes. We also showed that the use of an accurate DEM resulting from LIDAR data, based on present-day topography, leads to the underestimation of soil redistribution rates by the model, especially in this landscape submitted to recent and important morphological changes. Our results have important implications for the simulation of tillage erosion processes and our understanding of soil redistribution processes in complex cultivated areas. This is of particular interest to improve our knowledge and prediction of patterns of soil physical parameters, such as carbon storage or water content, particularly sensitive to surface erosion and landscape structuration.
•Study of 137Cs residuals to assess the impact of land consolidation on soil erosion.•Conversion of depositional areas into sediment delivering areas.•Intensive erosion and deposition concentrated on present and former field borders.•Soil redistribution largely dominates by tillage-induced processes.•Soil erosion on present field borders hardly simulated by conversion model.
The yields of the tropical rivers of Southeast Asia supply large quantities of carbon to the ocean. The origin and dynamics of particulate organic matter were studied in the Houay Xon River catchment ...located in northern Laos during the first erosive flood of the rainy season in May 2012. This cultivated catchment is equipped with three successive gauging stations draining areas ranging between 0.2 and 11.6 km2 on the main stem of the permanent stream, and two additional stations draining 0.6 ha hillslopes. In addition, the sequential monitoring of rainwater, overland flow and suspended organic matter compositions was conducted at the 1 m2 plot scale during a storm. The composition of particulate organic matter (total organic carbon and total nitrogen concentrations, δ13C and δ15N) was determined for suspended sediment, soil surface (top 2 cm) and soil subsurface (gullies and riverbanks) samples collected in the catchment (n = 57, 65 and 11, respectively). Hydrograph separation of event water was achieved using water electric conductivity and δ18O measurements for rainfall, overland flow and river water base flow (n = 9, 30 and 57, respectively). The composition of particulate organic matter indicates that upstream suspended sediments mainly originated from cultivated soils labelled by their C3 vegetation cover (upland rice, fallow vegetation and teak plantations). In contrast, channel banks characterized by C4 vegetation (Napier grass) supplied significant quantities of sediment to the river during the flood rising stage at the upstream station as well as in downstream river sections. The highest runoff coefficient (11.7%), sediment specific yield (433 kg ha−1), total organic carbon specific yield (8.3 kg C ha−1) and overland flow contribution (78–100%) were found downstream of reforested areas planted with teaks. Swamps located along the main stream acted as sediment filters and controlled the composition of suspended organic matter. Total organic carbon specific yields were particularly high because they occurred during the first erosive storm of the rainy season, just after the period of slash-and-burn operations in the catchment.
Soil mixing and the downward movement of solid matter in soils are dynamic pedological processes that strongly affect the vertical distribution of all soil properties across the soil profile. These ...processes are affected by land use and the implementation of various farming practices, but their kinetics have rarely been quantified. Our objective was to investigate the vertical transfer of matter in Luvisols at long-term experimental sites under different land uses (cropland, grassland and forest) and different farming practices (conventional tillage, reduced tillage and no tillage). To investigate these processes, the vertical radionuclide distributions of super(137)Cs and super(210)Pb (xs) were analyzed in 9 soil profiles. The mass balance calculations showed that as much as 91 plus or minus 9% of the super(137)Cs was linked to the fine particles ( less than or equal to 2 mu m). To assess the kinetics of radionuclide redistribution in soil, we modeled their depth profiles using a convection-diffusion equation. The diffusion coefficient represented the rate of bioturbation, and the convection velocity provided a proxy for fine particle leaching. Both parameters were modeled as either constant or variable with depth. The tillage was simulated using an empirical formula that considered the tillage depth and a variable mixing ratio depending on the type of tillage used. A loss of isotopes due to soil erosion was introduced into the model to account for the total radionuclide inventory. All of these parameters were optimized based on the super(137)Cs data and were then subsequently applied to the super(210)Pb (xs) data. Our results show that the super(137)Cs isotopes migrate deeper under grasslands than under forests or croplands. Additionally, our results suggest that the diffusion coefficient decreased with depth and that it remained negligible below the tillage depth at the cropland sites, below 20 cm in the forest sites, and below 80 cm in the grassland sites.
Numerous villages in the European loess belt are confronted with floods caused by runoff from agricultural land. Seventy-nine percent of the municipalities in central Belgium experienced at least one ...muddy flood during the last decade. Of these flooded municipalities, 22% have been affected more than 10 times during this period. Twenty municipalities have been selected for a detailed analysis. A database of 367 locations affected by muddy floods has been compiled, and the connectivity between cultivated areas and inhabited zones could be assessed for 100 flooded locations. Roads and drainage network facilitate runoff transfer between cultivated and inhabited areas in 64% of cases. Three types of areas producing muddy floods have been identified: hillslopes (1–30 ha) without thalweg where runoff is generally dominated by sheet flow; small catchments (10–300 ha) characterised by runoff concentration in the thalweg and medium catchments (100–300 ha) with multiple thalwegs dominated by concentrated runoff. About 90% of muddy floods are generated on hillslopes and in small catchments. A critical area–slope threshold for triggering muddy floods has been computed for hillslopes. A logistic regression shows that muddy floods are generated in small and medium catchments with 99% probability after 43 mm rainfall. Rainfall depths required to trigger muddy floods are lower in May and June (25
±12 mm) than between July and September (46
±
20 mm), because of different surface conditions (crusting, roughness and crop cover). Each year, muddy floods lead to a total societal cost of 16
×
10
6–172
×
10
6 € in central Belgium, depending on the extent and intensity of thunderstorms and monetary values damaged. Recent datasets suggest that the phenomenon is becoming more frequent in central Belgium, because of land consolidation, urban sprawl and expansion of row crops, sown in spring, at the expense of winter cereals. The huge costs induced by muddy floods justify the installation of erosion control measures. It is suggested to install a grassed buffer strip at the downslope edge of cultivated hillslopes to protect houses and roads. In small and medium catchments, it is preferred to install a grassed waterway and earthen dams in the thalweg.
► Tracing the origin of suspended matter is one of the main challenges in hydrology. ► We develop a new spectroscopic approach for tracing the sources of sediment. ► The DRIFTS signature is ...discriminant of the lithological sources. ► The DRIFTS signature is conservative in the hydrological network. ► The DRIFTS–PLS approach allows quantitative reconstitution of sediment sources.
In mountainous catchments, large quantities of sediment are exported within very short periods leading to numerous environmental problems (e.g. reservoir siltation). The origin of suspended sediment during two distinct floods was determined by conducting an original fingerprinting method coupling Diffuse Reflectance Infrared Fourier Transform Spectroscopy (DRIFTS) and a chemometric technique (i.e. Partial Least Squares – PLS-analysis). Samples of the potential sediment sources were collected in badland areas developed on various substrates (i.e. molasse, marly limestones, black marls and gypsum) in the Galabre 20
km
2-catchment located in the French Southern Alps. DRIFTS spectra provided a way to discriminate between the different potential sediment sources. Furthermore, the use of mid-infrared spectra allowed the direct quantification of the gypsum proportion in sediment. This contribution was systematically null at the catchment outlet because of the rapid dissolution of gypsum in the river. A PLS model was then constructed to estimate the contribution of the three other potential sources to the sediment flux during the floods. This model was developed and validated using a set of 45 “experimental” samples that were prepared in the laboratory in order to contain various proportions of the three remaining sources. By introducing DRIFTS spectra into the PLS model, we could predict the proportions of those sources in the mixed ‘experimental’ samples with a confidence interval of ca. ±10%. The model was then applied to the sediment collected during the two selected floods in order to outline their origin. Black marls provided the highest contribution of sediment during both events, but the analysis also revealed a significant contribution of molasse. Results also showed the remobilisation of sediment originated from molassic substrates that deposited on the riverbed during a preceding event. Opportunities for improvement and further use of this method as an alternative or rapid complementary sediment fingerprinting technique are finally discussed.