Definitions and generating methods of selected topographic metrics.
Isotopic tracers are useful for assessing linkages between soil movement and soil carbon dynamics in landscapes. Analyses of isotopes and comparison of isoscape (isotopic landscape) with observational data have been employed to investigate spatial distributions of isotopes, to test efficiencies of isotopic models, and to examine soil redistribution patterns and C dynamics. This chapter reviewed the application of natural (7Be, 210Pb) and anthropogenic fallout radionuclides (137Cs, 239,240Pu), and C isotopes (12,13,14C) in understanding mechanisms of soil redistribution and sedimentation. The chapter was organized to cover the formation, sources, and transport of these isotopes; how they are distributed and related to soil redistribution on C dynamics; and importance of their distribution (isoscapes) on investigating soil properties. We also provided a case study to demonstrate the feasibility of applying isotopes and isoscape modeling for understanding soil property variability in response to anthropogenic disturbance in a low-relief cropland field. Results demonstrated advantages of using 137Cs and C isotopic signature (δ13C) to trace soil movements and C dynamics. Topography-based 137Cs and C isoscape models were developed using light detection and ranging data (LiDAR) derived topographic metrics. The models successfully simulated the spatial patterns of 137Cs inventory and δ13C over an agricultural landscape and can benefit soil sedimentation and C dynamic studies in areas with limited observations.
- isotope analysis
- soil redistribution
- soil carbon dynamics
- dynamic replacement
Soil redistribution (i.e., erosion and deposition) exerts a significant impact on the global carbon (C) cycle. Movement of soil particles could either redistribute soil C or change the C mineralization through disruption of soil aggregation, reaggregation of deposited soils, and deep burial of C-enriched sediments at depositional sites . When soil erosion occurs, the surface concentrated and low-density proportions of sediments, such as soil organic carbon (SOC), are preferentially removed by runoff, wind, and/or tillage activities. A greater SOC enrichment ratio in eroded sediment is usually observed compared to its origin . The eroded sediment is either redistributed over a landscape or deposited at depressional sites. Therefore, investigation of soil redistribution and underlying mechanisms are needed for better understandings of the fate of soil C within the landscape.
Various methods have been developed to quantified soil redistribution. Before the 1990s, researchers mainly focused on understanding soil erosion by discerning changes in soil texture and chemical properties . However, accuracy of this method was low, mainly because impacts of erosion on soil properties are complicated and vary with soil characteristics, soil depth, and local land use practices. Other studies adopted erosion prediction models, such as the USLE, WaTEM, and WEPP, to assess soil redistribution processes [4, 5]. Although these models reasonably estimated long-term mean erosion rates over large-scale, the estimates from the modeling are often mean redistribution rates at sites, model estimates from individual erosion events are too coarse to be linked to each soil samples to explain changes in soil properties .
The emergence of isotopic analysis allows researchers to accurately trace soil movement at locations, offering a potential way to quantify the impacts of soil erosion and deposition on soil properties. The isotopic tracer selection follows two criteria: (1) isotopes should be able to quickly and strongly absorbed by the soil; and (2) the variability in absorption to various textures or size is either minor or can be calculated . A tracer can be uniformly distributed at first and then move with the soil movements, presenting mass variations between eroded and depositional sites.
This chapter examines the application of six widely used tracers, including anthropogenic fallout radionuclides/Cesium-137 (137Cs), plutonium isotopes (239,240Pu), natural fallout radionuclides Beryllium-7 (7Be), Lead-210 (210Pb), and C isotopes (12,13,14C), in soil redistribution and C distribution studies. Objectives of this chapter are to (1) review impacts of soil redistribution on soil C dynamics; (2) summarize critical processes regulating the selected soil isotopes; (3) introduce the concept of isoscape and its applications in soil biogeochemical studies; and (4) provide discussion of isotopic and isoscape analyses for understanding of soil redistribution and carbon dynamics through a case study.
2. Soil redistribution impacts on soil C dynamics
The global soil inventory represents an important C pool (Figure 1) with a total of 1950 Gt C, accounting for 2.3 times the size of the atmospheric C pool (860 Gt C) and 3.5 times of the biotic C pool (550 Gt C) . About 1750 Gt C of total soil C concentrates on land surface, where soil erosion and deposition processes are intense. The light and fine soil particles with high SOC content are preferentially removed by the erosion process and redistributed over a landscape or deposited at depressional sites. It was estimated that 75 Gt year−1 of soil is removed by water and wind erosion , which contribute to 0.8–2.2 Gt year−1 emission of C from land surface, and 2.5–3.9 Gt year−1 of replaced C in soil .
Although researchers have increasingly recognized the importance of soil erosion and C dynamics on ecosystems, mechanisms of soil redistribution controls on soil C remains poorly understood. Some studies suggested negative effects of soil erosion on terrestrial C sequestration [9, 10, 11]. They provide an argument that excessive soil erosion can lead to losses of soil fertility, decreasing the plant and crop productivity . Notably, the recently accelerated soil erosion due to agricultural activities has caused 2 billion ha of land being irreversibly degraded . Plowing activity increases the possibility of soil organic matter within the plow layer to be exposed to anthropogenic and climatic perturbations, and thus accelerates SOC mineralization . Furthermore, soil C mineralization in displaced soils can also be stimulated as a result of the breakdown of soil aggregates by raindrop and runoff shearing forces that expose previously protected organic matters to oxidation . It was estimated that more than 20% C emission was caused by mineralization of the displaced C .
On the contrary, arguments for positive impacts of soil redistribution on C storage support that redistribution of soil could increase C sequestration in terrestrial ecosystems [14, 15, 16, 17]. SOC losses at the eroded sites could be dynamically replaced by litter input from plant regrowth and return of above- and below-ground biomass to soil and replenish the depleted SOC pool at eroded sites [18, 19]. The eroded SOC is subsequently buried and re-aggregated at depressional sites, protecting SOC from mineralization and thus increasing C sequestration [14, 20, 21]. Stallard  used a large set of scenarios estimated that the human-induced C burial is in the range of 0.6–1.5 Gt C year−1 by terrestrial sedimentation.
3. Using isotopes in studies of soils and sediments
3.1 Anthropogenic fallout radionuclides
137Cs is the most widely used isotope in geomorphic studies of water, wind and tillage erosion. As an anthropogenic radionuclide with a half-life of 30.2 years, 137Cs was released globally into the environment due to radioactive fallout from nuclear weapon testing in the 1950s and 1960s and reached to land surface through wet and dry deposition (Figure 2). Additional 137Cs fallout occurred because of the 1986 Chernobyl accident in the north of the Ukrainian SSR and the 2011 Fukushima accident in Japan. The overall amount of 137Cs released from Fukushima was less than 15% of the amount from the Chernobyl accident . The Fukushima accident added high amounts of 137Cs to soils near the accident site, but the large-scale plume of radioactive fallout extended mainly over the Pacific Ocean which reduced the extent of terrestrial labeling. In contrast, Chernobyl is in the center of the European continent. The released fallout radionuclides significantly increased the 137Cs amount in many European countries . 137Cs is chemically active and is rapidly and strongly absorbed by fine soil particles when it contacts with the soil. Vertical migration is slow and the majority of 137Cs is retained in the upper 20 cm of the soil surface across the globe [23, 24, 25]. Once the soil is labeled by 137Cs, chemical and biological removal of 137Cs is low. The 137Cs concentration changes mainly result from physical processes in the top-layer soils, such as soil transport and deposition. Therefore, these characteristics make 137Cs a useful tracer for quantifying of soil redistribution rates.
Previous investigations have reported strong and statistically significant correlations between 137Cs and soil C [10, 26, 27]. However, because 137Cs inventory is usually non-linearly correlated with soil redistribution, the correlations between 137Cs inventory and soil C content are usually weaker than the correlations between soil redistribution process and SOC . Consequently, instead of directly quantify SOC based on 137Cs, most recent studies first converted 137Cs measurements to soil redistribution rates, and then used the estimates to link to soil C content. Specifically, the application of 137Cs measurements in soil redistribution and C dynamics mainly includes to (1) calculate soil redistribution at a point and link to soil properties of a soil sample from the same point [6, 28, 29]; (2) provide mean annual rate of soil redistribution over 60 years and reflect the erosion history of a site [10, 29]; and (3) present spatial patterns of soil erosion and deposition over the landscape [27, 30].
Because of the short half-life, 137Cs concentration may drop below the detection limits quickly and limit use of this tracer in SOC investigations. Researchers are finding alternative tracers with longer half-lives for soil erosion investigation. Anthropogenic radionuclides of two major plutonium isotopes (239Pu: half-life = 24,110 years and 240Pu: half-life = 6561 years) are considered as potential alternatives for this purpose. Like 137Cs, the Pu isotopes in soils are also mainly from nuclear weapon testing between the 1950s and 1960s following the 1986 Chernobyl accident. However, the isotopic composition of Chernobyl Pu (240Pu/239Pu atom ratio: 0.37–0.41) presents higher ratios compared to Pu from nuclear testing (240Pu/239Pu atom ratio = 0.180 ± 0.014), which can be used to distinguish Pu isotopes from the two sources . Some researchers have argued that in soils with high organic matter, Pu might be preferentially absorbed by organic fractions with high molecular weight and potentially lead to a deep vertical migration [32, 33].
Considering the long half-life of Pu isotopes, increasing studies have explored the feasibility of using Pu isotopes in investigations of soils and sediments [34, 35, 36]. Schimmack et al.  found that the Chernobyl Pu could be a better tracer for soil redistribution detections in Bavaria, Germany, since spatial patterns of Pu agreed better with soil mapping than that of 137Cs. Similar conclusions were reached in several other studies, reporting more homogenous distributions of239+240Pu than 137Cs in parts of Europe that were influenced by the Chernobyl accident [31, 38]. Thus, the isotopic composition variations may lead Pu isotopes to be a better choice than 137Cs for tracing soil erosion. However, the feasibility of applying 239,240Pu to soil erosion and soil C investigations still needs further testing because of its vertical migration and lateral transport by water, which may influence the quantification of soil redistribution [39, 40].
3.2 Natural radionuclides
7Be is a natural radionuclide produced primarily through cosmic-ray spallation of oxygen and nitrogen nuclei in the stratosphere and troposphere. Unlike the anthropogenic radionuclides of 137Cs and 239,240Pu, wet and dry deposition of 7Be to the soil surface is continuous because of its natural origin . The concentration of 7Be in the soil decreases with increases in particle size . Vertical migration was low for this isotope due to its relatively short half-life (half-life = 53.3 days). Wallbrink and Murray  found that most of soil 7Be was accumulated in the top 20 mm regardless of soil types and surface cover conditions at two Australian sites.
Because of its short half-life, 7Be has the potential to trace short-term soil redistribution. A number of studies have used this isotope to examine short-term soil erosion after heavy rain events, providing a basis for understanding sediment transported by dispersed overland flow [41, 42, 43, 44]. Considering the lower vertical transport of 7Be relative to 137Cs and excess 210Pb (210Pbex), this isotope has been applied to discriminate vertical and horizontal erosion by combining with the other two isotopes [45, 46]. Ryken et al.  suggested correction factors related to particle size and variations in relaxation mass depth should be used to get precise estimates of 7Be-derived soil redistribution rates. Li et al.  successfully applied 7Be to quantify SOC sequestration changes caused by land-use and management activity.
Another natural isotope, Lead-210 (210Pb) produced through the decay of gaseous Radon-222 (222Rn) generated from the decay of Radium-226 (226Ra) in the Uranium-238 (238U) decay series. 226Ra exists in soils and rocks. Most of 226Ra decay to 210Pb in situ, which is termed supported 210Pb. Due to diffusion, a small portion of 226Ra-derived 222Rn enters to the atmosphere and subsequently introduces 210Pb to the atmosphere. This kind of 210Pb falls to land surface through wet and dry deposition and is termed unsupported or excess 210Pb (210Pbex). The 210Pbex is strongly absorbed by fine soil particles and transports with soil movement. 210Pbex fallout is continuous over time, which is like 7Be, but its half-life (half-life = 22.3 years) is longer than 7Be. Therefore, 210Pbex has the potential to penetrate deeper soil layers (to 10 cm) than 7Be [7, 23].
The use of 210Pbex in soil and sediment studies has increased in recent decades. Due to the continuous fallout, 210Pbex can be used to provide long-term soil redistribution rates. Meanwhile, unlike 137Cs, 210Pbex does not have below-detection-limit problems caused by medium-lived anthropogenic fallout radionuclides. In practice, two or more different fallout radionuclides can be used to understand the soil erosion history. 210Pbex is combined with 137Cs and/or 7Be to provide soil redistribution records in the past 100 years. Several studies suggested that 210Pbex and 7Be can produce similar spatial patterns [45, 49], but differences resulting from different land uses could help sediment source identification [46, 50]. Increasing applications of 210Pbex in C dynamic studies have been reported recently [10, 15, 51, 52]. Investigations have demonstrated that 210Pbex is preferentially associated with SOC than 137Cs, due to the stronger binding to the organic matter in soils [53, 54].
3.3 C isotopes
The origin of soil C isotopes (12C, 13C, and 14C) is mainly from plant litter entering the soils. During photosynthesis, plant species absorb atmospheric carbon dioxide (CO2), which consists of about 98.9% of 12CO2, 1.1% of 13CO2, and trace amount of 14CO2 (1 part in trillion). 14C is a natural radionuclide created primarily by the cosmic-ray spallation of nitrogen in the troposphere and the stratosphere. 14C is a relatively stable radioactive isotope with a half-life of 5730 years . This C isotope is assimilated during plant photosynthesis and entering soil through litter fall and commonly used for age dating.
Data on the stable C isotopes 12C and 13C are usually reported δ13C representing deviation of measured 13C from the established natural abundance with units of parts per thousand (%). Variations in soil δ13C values are controlled primarily by carbon input from plant litter. Due to isotopic discrimination by their photosynthetic enzymes and the regulation of stomatal diffusion resistance, δ13C values in plants cover a wide range. Plants with C3 photosynthesis have δ13C values in a range of −22 to −32%. Plants with C4 photosynthesis are less depleted in 13C with higher δ13C values ranging from −9 to −17%. After plant residues entering the soil, the δ13C values may change slightly due to isotope fractionation during microbial decomposition.
Soil C isotopes can effectively detect soil redistribution and reflect soil C dynamics over the landscape. Because of the isotopic variability in different plant species, soil δ13C has been successfully applied to distinguish eroded soil sources and to identify soil sources from different land use types [56, 57]. Furthermore, due to the significant difference of C isotopic (δ13C and 14C) values between surface soils and subsoils, scientists also utilized the isotopes to obtain a better understanding of the origins and age of eroded SOC [2, 57]. Trends of δ13C and 14C by soil depth can be used to reconstruct history of vegetation succession based on different isotopic compositions [55, 58]. Li et al.  used both δ13C and 137Cs to investigate soil C fate in an agricultural field under an annual crop rotation of soybean (
4. Use of isoscapes in studies of soils and sediments
The spatial and temporal features of isotopes in various environmental materials can be predicted using both landscape and biogeochemical models . The term of isoscape was introduced in the 2000s to describe spatial patterns of isotopes produced by various quantitative models using several environmental variables as predictors [60, 61]. Now isoscape analyses are widely used as a tool to understand the biogeochemistry of landscapes under study. Applications of isoscapes on soils and sediments can be categorized into three groups:
Isotopes can be used as the baseline in modeling studies. With the help of a robust model, regional isotopes can be estimated from gridded environmental datasets. Isoscapes provide a powerful tool for probing the efficiencies of these prediction models ;
Comparison of isotopic spatial patterns with observed soil organic C and nitrogen (N) is useful for investigating the impacts of soil redistribution on soil properties. Isoscapes can effectively reflect soil redistribution patterns across space, which can be used to better explain spatial variability of soil properties [29, 63, 64];
Soil isoscapes is also used in investigations of plant and animal isoscapes. For example, because the primary input of 12C and 13C to soil are from plant litter, maps of soil δ13C isoscape could provide information related to the biogeography of C3 and C4 plants and their relative contributions to the atmospheric C sequestration .
5. Application of isotopic analysis in an agricultural field
5.1 Study area introduction
A case study is presented to demonstrate the application of isotopic analysis and isoscape modeling in investigation of soil redistribution and C dynamics in cropland. The study was conducted in an agricultural field in Walnut Creek Watershed (WCW), Iowa. The WCW is in a humid continental climatic zone with relatively flat terrain (Figure 3a). Representative soils are a poor-drained Nicollet (mesic Aquic Hapludolls) in the lowlands and well-drained Clarion (mesic Typic Hapludolls) in the uplands. The dominant land use type in the WCW is agriculture. More than 86% of the watershed is farmlands. Primary tillage practices in the watershed are chisel plowing in autumn and disking in spring. The directions of tillage depend on the farm management and field configuration. Most of the area followed the north-south or east-west tillage directions.
The studied cropland field is about 15 ha (Figure 3b). The maximum elevation change within a radius of 90 m is 4.5 m (Figure 3c). A C4–C3 crops rotation (i.e., maize-soybean) has been adopted in this area since 1957, which make this field an appropriate testbed for a detailed C isotope study assessing impacts of soil redistribution.
5.2.1 Sampling and laboratory analyses
Isotopes of 137Cs, 12C, and 13C were measured in soil samples. A total of 119 locations were collected at grid nodes in a 25 × 25 m grid. A 9-sample and an 11-sample transects were set at an interval of 12.5 m within the grid (Figure 2b and c). We also selected four reference sites for estimating the baseline 137Cs inventory. Three samples were collected from the top 30 cm of soil within a 1 × 1 m quadrate using a 3.2 cm diameter push probe at each location. All soil samples were oven dried and weighed for soil bulk density calculation.
SOC content was estimated as the difference between total C that was measured by a LECO CNS 2000 elemental analyzer and calcium C estimated using soil samples ashed in a furnace at 420°C. SOC content (%) was them converted to SOC density (kg m−2) using soil bulk density.
The 13C/12C isotope ratio was measured using an isotope-ratio mass spectrometer (Europa Scientific Ltd., Crewe, England). The ratio was then combined with the laboratory reference that was calibrated against the international PeeDee Belemnite (PDB) to estimated δ13C related to PDB.
This study estimated the relative contribution of C3 and C4 crops using the following equations:
where FC4 and FC3 are contribution fraction of C4- and C3- derived SOC, respectively.
To measure 137Cs, gamma-ray analysis was performed using a Canberra Genie-2000 Spectroscopy System. Original unit of 137Cs concentration was becquerels per gram (Bq g−1), and this unit was subsequently converted to becquerels per square meter (Bq m−2) using soil bulk density. The baseline 137Cs inventory estimated from reference sites and the 137Cs inventory of sampling sites were used to calculate 137Cs-derived soil redistribution rates using the Mass Balance model I developed by Walling et al. . In this study, the baseline 137Cs inventory was 2657 Bq m−2. Sites with higher 137Cs inventories than the baseline were considered as depositional sites; while sites with lower 137Cs inventories were referred as eroded sites.
Additional sampling and laboratory analyses details are given in Li et al. .
5.2.2 Historic orthophoto interpretation
Historic orthophotos in the 1950s and 2002 were obtained from the Iowa geographic map server (http://ortho.gis.iastate.edu/) to visualize soil movement. The surface soil layer was referred as mollic epipedon formed under grass vegetation that dominated this area prior to the 1850s. The mollic epipedon has a characteristic black color and contains a high amount of soil organic matter. Therefore, the black mollic epipedon can be visually identified in the historic orthophoto, serving as an indicator for soil erosion and deposition investigation after the 1850s. In this study, greyness values were derived from the orthophotos to quantify the intensity of black carbon signal related to the mollic epipedon. The values were compared to 137Cs-derived soil redistribution rates to assess the appropriateness of using 137Cs to trace soil movement in this low-relief cropland.
5.2.3 Terrain analyses
Topographic metrics used to describe features of the landscape were generated from the light detection and ranging (LiDAR)—derived 3 m digital elevation model (DEM). Fourteen topographic metrics were developed, including slope, profile curvature (P_Cur), plan curvature (Pl_Cur), general curvature (G_Cur), flow accumulation (FA), topographic relief [topographic relief principal components (TRPCs) and topographic relief factors (TRFAs)], positive topographic openness (PTO), downslope index (DI), flow path length (FPL), catchment area (CA), topographic wetness index (TWI), stream power index (SPI), upslope slope (UpSl), and flow length factor (LS) (Table 1). Before calculating topographic metrics, the DEM was filtered twice through a low pass 3-by-3 filter. Maps of slope, P_Cur, Pl_Cur, G_Cur, PTO, DI, FPL, CA, TWI, SPI, UpSl, and LS were developed by the System for Automated Geoscientific Analyses (SAGA) v. 2.2.5 . FA map was generated with ArcGIS.
|Slope (radian)||An angle between a tangent and a horizontal plane at a given point||Slope, aspect, curvature module, SAGA|
|P_Cur (m−1)||Curvature of the surface in the direction of the steepest slope|
|Pl_Cur (m−1)||Curvature in a horizontal plane|
|G_Cur (m−1)||Curvature of the surface itself|
|FA (m2)||Land area that contributes surface water to an area in which water accumulates||ArcGIS|
|TRPC and TRFA||Topographic relief is elevation difference between the highest point over an area and a given location. TRPC/TRFA is topographic relief principal component/Topographic relief factor and is generated by several reliefs that are at different spatial-scales using principal component analysis/factor analysis||ArcGIS|
|PTO (radian)||An angular measure of the relation between surface relief and horizontal distance||Topographic openness module, SAGA|
|DI (radian)||Head differences along a flow path||Downslope distance gradient module, SAGA|
|FPL (m)||Maximum distance of water flow to a point in the catchment||Flow path length module, SAGA|
|CA (m2)||Area draining to catchment outlet||SAGA wetness index module, SAGA|
|TWI||Frequencies and duration of saturated conditions|
|SPI||Erosive power of overland flow||Stream power index module, SAGA|
|UpSl (radian)||Mean slope of upslope area||LS-factor (field based) module, SAGA|
|LS||Erosive power of the terrain|
For TRPCs, a series of topographic relief maps with different spatial scales were generated. The topographic relief shows the elevation differences between the filtered DEM and a maximum elevation map showing a continuous surface of maximum elevation within a specific area/radius. Seven relief images were generated from seven maximum elevation maps using radiuses of 7.5, 15, 30, 45, 60, 75, and 90 m. Principal component analysis (PCA) was conducted to convert the seven reliefs to seven independent components, and the first two topographic relief principal components (TRPC1 and TRPC2) were selected for further analysis. Similarly, varimax rotated Factor Analysis (FAn) was used and the first two TRFA1 and TRFA2 were selected. The detailed topographic processing can be found in Li et al. .
5.2.4 Statistical analysis and model calibration
Duncan’s multiple range tests (
Models including the multiple linear regression (MLR) combined with principal component analysis (MLR-PCA) and MLR combined with factor analysis (MLR-FAn) models for SOC, 137Cs, and δ13C were developed using stepwise linear regression with the “leave-one-out” cross-validation. We randomly selected 70% samples as training dataset used for model calibration and 30% samples as testing dataset for model validation. Three criteria were used to evaluate model efficiencies. The three criteria are the adjusted coefficient of determination (), the Nash-Sutcliffe efficiency (NSE), and ratio of the root mean square error (RMSE) to the standard deviation of measured data (RSR). Usually, the model performance is considered satisfactory when NSE is larger than 0.5 and RSR is smaller than 0.7 .
5.3 Results and discussion
5.3.1 Soil redistribution impacts on C dynamics
Distribution of mollic epipedon, characterized by presence of black soil organic matter, reflected historical soil movement (Figure 4). The top-layer black soil presumably blanketed the prairie landscape and was then transported by water and tillage from eroded upslope to downslope locations. This phenomenon can be observed from the higher mollic epipedon in the concave than the convex locations in the 1950s and likely more consolidated by 2002. The significant correlations between greyness values and 137Cs-derived soil redistribution rate suggested that 137Cs inventory can effectively reflect the soil redistribution process in this field.
The cropland experienced a general soil export as indicated by a negative mean 137Cs-derived soil redistribution rate (−0.51 ± 2.09 kg m−2 year−1, Table 2). Specifically, about 81 locations were identified as eroded sites with negative soil redistribution rates (mean: −1.72 ± 1.38 kg m−2 year−1) and 47 locations were depositional sites with positive soil redistribution rates (mean: 1.57 ± 1.32 kg m−2 year−1). SOC density exhibited a similar spatial pattern as soil redistribution rates with a higher mean value at depositional than eroded sites. The mean SOC density at depositional sites was about 1.70 times of that at eroded sites.
A strong and significant correlation was observed between 137Cs-derived redistribution rates and SOC density (
Note that SOC mineralization and dynamic replacement can complicate impacts of soil redistribution on SOC. On the one hand, SOC mineralization highly depends on soil moisture conditions. Impacted by soil texture, soil moisture is commonly lower at upslope areas than low-lying areas. Our study area also showed similar spatial patterns in soil texture over the landscape . The decreased fine particles (clay and silt) can reduce soil water retention, decrease aggregate stability, and increase oxygen concentration, accelerating SOC mineralization and further exacerbate SOC depletion in eroded areas . On the other hand, the depleted SOC pool can be dynamically replaced by the continually deposited SOC from litter decomposition of above- and below-ground biomass. The net sequestration of carbon is not always linear over time because rates of replaced carbon at eroded sites depend on management history that may vary with time .
Values of δ13C showed opposite spatial patterns with SOC and 137Cs-derived soil redistribution rates with a higher mean value of δ13C at eroded than depositional sites. The mean δ13C value at the eroded sites was −19.6 ± 2.05%; while the δ13C value at the depositional sites was −21.9 ± 1.52%. The higher mean δ13C value at eroded sites indicated that the eroded sites were less depleted for 13C and exhibited stronger C4 vegetation characteristics than depositional sites. The δ13C value variability over this space possibly resulted from dynamic replacement. Our study area was dominated by prairie vegetation with strong C3 vegetation characteristics before the 1850s. After that, C4-derived SOC was introduced as a new SOC formation into the soil due to the widespread maize cultivation . The new soil C (C4-derived SOC) replaced 50% of native soil C (C3-derived SOC) at eroded sites and 30% at depositional sites till 2002 when sampled.
According to the calculations of the relative contributions of C3 and C4 crops to the SOC, the C3-derived SOC density showed a significant difference between eroded (396 ± 202 kg m−2) and depositional sites (852 ± 226 kg m−2) while the C4-derived SOC density varied insignificantly (Table 2, eroded: 363 ± 131 kg m−2; depositional: 440 ± 152 kg m−2). C3-derived SOC density was positively related to 137Cs-derived soil redistribution rates with a coefficient of determination of 0.65 (Figure 5b). In contrast, C4-derived SOC density was not strongly related to soil redistribution rates (Figure 5c). The different responses to soil erosion may be related to chemical compounds contained in the C3- and C4-derived SOC. Generally, liable compounds that are high in C4 crop residues are less depleted in 13C than lignin-based compounds that are high in C3 crop residues . Therefore, C4-derived SOC may largely be found in the mollic epipedon and be significantly affected by C mineralization, resulting in a decreased correlation with patterns of soil movement. In contrast, the more recalcitrant SOC in soil with strong C3 vegetation characteristics would be less influenced by C mineralization because it contains higher amount of stable lignin-based compounds. Therefore, the variability in C3 signature of soil C would more strongly correlate with patterns of soil redistribution.
|Erosion||81||759 (246) b†||−19.6 (2.05)a||1791 (532)b||−1.72 (1.38)b||396 (202)b||363 (131)a|
|Deposition||47||1292 (259)a||−21.9 (1.52)b||3547 (769)a||1.57 (1.32)a||852.2 (226)a||440 (152)a|
|All||128||956 (359)||−20.4 (2.16)||2435 (1056)||−0.51 (2.09)||565 (305)||391 (143)|
5.3.2 Isoscape model development
Based on the derived topographic metrics, the first six TPCs and six TFAs were selected to construct independent topographic variables. The metrics with the highest absolute loadings for the first six TPCs were TWI, TRPC2, FA, FPL, DI, and CA. Therefore, the TPCs 1, 2, 3, 4, 5, and 6 were associated with soil wetness, flow acceleration, runoff volume, flow velocity, downslope dispersal, and water volume, respectively . The highest absolute loadings in the six TFAs were Upsl, G_Cur, Pl_Cur, FPL, CA, and DI, which were associated with upslope flow velocity, flow acceleration, flow convergence and flow divergence, flow velocity, water volume, and downslope dispersal, respectively.
SOC and isoscape models of δ13C and 137Cs were developed using MLR-PCA and MLR-FAn (Table 3). Both types of models (MLR-PCA and MLR-FAn) presented satisfactory efficiencies during model calibration, but more independent variables were selected in MLR-FAn models. For SOC simulation, the MLR-PCA model selected four components with a of 0.668, an NSE of 0.682, and an RSR of 0.563. The MLR-FAn model exhibited higher and NSE values and a lower RSE value than MLR-PCA SOC model. In contrast, for 137Cs inventory and δ13C simulations, similar values were observed between the two types of models, but the MLR-FAn models showed higher NSE and RSR values, which suggested that the 137Cs and δ13C values derived from MLR-FA models were less deviated from the corresponding observed values.
|SOC||949-119TPC1 + 86.9TPC5 + 21.4TPC2 + 37.2TPC4†||0.668||0.682||0.563|
|137Cs||2466-327TPC1 + 240TPC5||0.597||0.606||0.628|
|δ13C||−20.3 + 0.677TPC1 + 0.250TPC3 + 0.181TPC2||0.614||0.626||0.611|
|SOC||950 + 189TFA4-147TFA2-112TFA1 + 88.9TFA6 + 82.1TFA3 + 50.2TFA5||0.663||0.686||0.560|
|137Cs||2459 + 447TFA4 + 306TFA3-346TFA1 + 316TFA6-301TFA2 + 150TFA5||0.587||0.615||0.621|
|δ13C||−20.3 + 1.26TFA1 + 0.717TFA2-0.741TFA4-0.406TFA6-0.317TFA5||0.614||0.635||0.604|
Although both types of models reasonably matched up with observations during model calibration, the MLR-PCA models had better performance in predictions than the MLR-FAn models when applied to the test dataset. Predicted SOC explained 62.0% variability in observed SOC using MLR-PCA model. The values of NSE and RSR were 0.612 and 0.622, respectively. MLR-PCA 137Cs and δ13C models also had satisfactory performance with NSE values larger than 0.5 and RSR values smaller than 0.7. The R2 for 137Cs and δ13C prediction were 0.713 and 0.509, respectively. In contrast, all the MLR-FAn models presented lower efficiencies in predicting target variables. Especially for the δ13C, the NSE value was smaller than 0.5 and RSR value was larger than 0.7 when compared MLR-FAn predictions with observations.
The lower efficiencies of the MLR-FAn models than the MLR-PCA models may be caused by model over-fitting. PCA considers all the variance in the independent variables, including unique, error and shared variance during synthetic variable construction; while FAn only considered and presented the shared variance in the factor matrix. Therefore, differences exist during development of MLR-PCA and MLR-FAn models. In this case study, although these two types of models had similar performance during model calibration using training datasets, increased number of parameters in MLR-FAn models may increase model error using validation datasets with general decrease in model stability, leading to low accuracies during extrapolating prediction points to external sample sets .
The spatial patterns of SOC density, 137Cs inventory, and δ13C were generated using MLR-PCA models (Figure 6). SOC density showed a similar spatial variability as 137Cs inventories, which is consistent with the strong correlations between the SOC density and 137Cs-derived soil redistribution rate (Figure 5a). Both variables had high values in depressions and low values in ridge and sloping areas. For δ13C, an opposite pattern was observed with low values in concave landscapes and high values in convex landscapes.
The high model efficiencies demonstrated the feasibility of using topography-based isoscape MLR-PCA models for investigating the spatial variability of isotopes. Although limitations may exist due to unaccounted topographic variables and environmental impactors, such as variance in tillage activities, the DEM-derived topography-based models provide a cost-effective method for isoscape mapping, effectively reflecting redistribution patterns of soil and spatial patterns of SOC over the landscape [28, 74, 75]. These models can also benefit investigations in regions with limited access to ground measurements. The recently increased availability of the fine-resolution LiDAR data can further improve the model applicability, allowing for scaling of in situ isotopic simulations from field scale to isoscape mapping at watershed or regional scales.
Soil isotopes can effectively trace soil redistribution and SOC dynamics. This chapter reviewed the application of natural (7Be, 210Pb), anthropogenic fallout radionuclides (137Cs, 239,240Pu), and C isotopes (12,13,14C) in understanding soil erosion and deposition at different spatial and temporal scales. The case study demonstrated that 137Cs and C isotopes could be employed to understand soil movement and C dynamics. SOC density showed high consistency with 137Cs-derived soil redistribution rate, suggesting significant spatial impacts of soil movement on SOC. δ13C provided further support of the importance of dynamic replacement on soil C dynamics in this area. Topography-based isoscape models were developed and effectively reconstructed the spatial variability in 137Cs inventory and δ13C over the landscape. The isoscape maps of 137Cs and δ13C were in general agreement with the spatial pattern of SOC density. Based on such results, we conclude that isotopic and isoscape analysis could provide valuable insights into soil movement and C studies.
This research was supported by the USDA Natural Resources Conservation Service in association with the Wetland Component of the National Conservation Effects Assessment Project (NRCS 67-3A75-13-177).
Conflict of interest
No conflict of interest exists relative to information presented in this chapter.