Physics-based modeling of cosmogenic nuclides part II – Key aspects of in-situ cosmogenic nuclide production

Characteristics of the spallogenic component of nuclide production are investigated through the use of a physics-based model. Calculated production rates for commonly used nuclides indicate differences in scaling up to 15% at very high altitude. Angular distribution of nuclide forming particles suggests the current method of shielding correction, which is neither altitude nor latitude dependent, can be improved on. Subsurface production profiles suggest that erosion corrections should be performed with non-constant attenuation lengths. Results are parameterized for easy application.     

Quantifying effects of deep and near‐surface chemical erosion on cosmogenic nuclides in soils,saprolite, and sediment

Cosmogenic nuclides in rock, soil, and sediment are routinely used to measure denudation rates of catchments and hillslopes. Although it has been shown that these measurements are prone to biases due to chemical erosion in regolith, most studies of cosmogenic nuclides have ignored this potential source of error. Here we quantify the extent to which overlooking effects of chemical erosion introduces bias in interpreting denudation rates from cosmogenic nuclides. We consider two end‐member effects: one due to weathering near the surface and the other due to weathering at depth. Near the surface, chemical erosion influences nuclide concentrations in host minerals by enriching (or depleting) them relative to other more (or less) soluble minerals. This increases (or decreases) their residence times relative to the regolith as a whole. At depth, where minerals are shielded from cosmic radiation, chemical erosion causes denudation without influencing cosmogenic nuclide buildup. If this effect is ignored, denudation rates inferred from cosmogenic nuclides will be too low. We derive a general expression, termed the ‘chemical erosion factor’, or CEF, which corrects for biases introduced by both deep and near‐surface chemical erosion in regolith. The CEF differs from the ‘quartz enrichment factor’ of previous work in that it can also be applied to relatively soluble minerals, such as olivine. Using data from diverse climatic settings, we calculate CEFs ranging from 1.03 to 1.87 for cosmogenic nuclides in quartz. This implies that ignoring chemical erosion can lead to errors of close to 100% in intensely weathered regolith. CEF is strongly correlated with mean annual precipitation across our sites, reflecting climatic influence on chemical weathering. Our results indicate that quantifying CEFs is crucial in cosmogenic nuclide studies of landscapes where chemical erosion accounts for a significant fraction of the overall denudation. Copyright © 2012 John Wiley & Sons, Ltd.     

Postglacial to Holocene landscape evolution and process rates in steep alpine catchments

Climate change and high magnitude mass wasting events pose adverse societal effects and hazards, especially in alpine regions. Quantification of such geomorphic processes and their rates is therefore critical but is often hampered by the lack of appropriate techniques and the various spatiotemporal scales involved in these studies. Here we exploit both in situ cosmogenic beryllium-10 (¹⁰Be) and carbon-14 (¹⁴C) nuclide concentrations for deducing exposure ages and tracing of sediment through small alpine debris flow catchments in central Switzerland. The sediment cascade and modern processes we track from the source areas, through debris flow torrents to their final export out into sink regions with cosmogenic nuclides over an unprecedented five-year time series with seasonal resolution. Data from a seismic survey and a 90 m core revealed a glacially overdeepened basin, filled with glacial and paraglacial sediments. Surface exposure dating of fan boulders and radiocarbon ages constrain the valley fill from the last deglaciation until the Holocene and show that most of the fan existed in early Holocene times already. Current fan processes are controlled by episodic debris flow activity, snow (firn) and rock avalanches. Field investigations, digital elevation models (DEMs) of difference and geomorphic analysis agree with sediment fingerprinting with cosmogenic nuclides, highlighting that the bulk of material exported today at the outlet of the subcatchments derives from the lower fans. Cosmogenic nuclide concentrations steadily decrease from headwater sources to distal fan channels due to the incorporation of material with lower nuclide concentrations. Further downstream the admixture of sediment from catchments with less frequent debris flow activity can dilute the cosmogenic nuclide signals from debris flow dominated catchments but may also reach thresholds where buffering is limited. Consequently, careful assessment of boundary conditions and driving forces is required when apparent denudation rates derived from cosmogenic nuclide analysis are upscaled to larger regions. © 2018 John Wiley & Sons, Ltd.     

Detection of transience in eroding landscapes

Past variations in climate and tectonics have led to spatially and temporally varying erosion rates across many landscapes. In this contribution I examine methods for detecting and quantifying the nature and timing of transience in eroding landscapes. At a single location, cosmogenic nuclides can detect the instantaneous removal of material or acceleration of erosion rates over millennial timescales using paired nuclides. Detection is possible only if one of the nuclides has a significantly shorter half‐life than the other. Currently, the only practical way of doing this is to use cosmogenic in situ carbon‐14 (¹⁴C) alongside a longer lived nuclide, such as beryllium‐10 (¹⁰Be). Hillslope information can complement or be used in lieu of cosmogenic information: in soil mantled landscapes, increased erosion rates can be detected for millennia after the increase by comparing relief and ridgetop curvature. This technique will work as long as the final erosion rate is greater than twice the initial rate. On a landscape scale, transience may be detected based upon disequilibria in channel profiles or ridgetops, but transience can be sensitive to the nature of transient forcing. Where forcing is periodic, landscapes display differing behavior if forcing is driven by changes in base level lowering rates versus changes in the efficiency of either channel or hillslope erosion (e.g. driven by climate change). Oscillations in base level lowering lead to basin averaged erosion rates that reflect a long term average erosion rate despite strong spatial heterogeneity in local erosion rates. This averaging is reflected in ¹⁰Be concentrations in stream sediments. Changes in hillslope sediment transport coefficients can lead to large fluctuations in basin averaged erosion rates, which again are reflected in ¹⁰Be concentrations. The variability of erosion rates in landscapes where both the sediment transport and channel erodibility coefficients vary is dominated by changes to the hillslope transport coefficient. Copyright © 2016 John Wiley & Sons, Ltd.     

Inferring exposure ages and erosion rates from cosmogenic nuclides: A probabilistic formulation

Cosmogenic nuclides have become an important tool in geomorphology; the concentration of such nuclides in minerals in an eroding surface is directly related to the exposure time and the erosion rate. In principle, measurement of the two nuclides ¹⁰Be and ²⁶Al allows for the determination of both the erosion rate and the exposure age. In practice, due to a variety of factors such as the similar lifetimes of the two nuclides and the limits on measurement precision, this determination is often not possible. We propose a new approach to this problem, showing how to construct a joint probability distribution for the age and erosion rate by using the concentration of one (or two) nuclides measured at the surface. We explain the Bayesian approach to this problem; the construction of the prior is a crucial element of this Bayesian method, and we devote particular attention to this issue. By analyzing previously published data, we show how this method improves on the standard approach of computing a “model age” and “model erosion rate.”     

Conceptual model of fracture-limited sea cliff erosion: Erosion of the seaward tilted flyschs of Socoa,Basque Country,France

Sea cliff morphology and erosion rates are modulated by several factors, including rock control that reflects both lithology and rock structure. Erosion is anticipated to preferentially exploit ‘fractures’, broadly meant as any discontinuity in an otherwise continuous medium, where the rock mass is weakest. Unpicking the direct control of such fractures on the spatial and temporal pattern of erosion remains, however, challenging. To analyse how such fractures control erosion, we monitored the evolution of a 400 m-long stretch of highly structured sedimentary cliffs in Socoa, Basque Country, France. The rock is known as the Socoa flysch formation. This formation combines decimetre-thick turbidites composed of repeat triplets of medium to strong calcareous sandstone, laminated siltstones and argillaceous marls. The sequence plunges at 45° into the sea with a shore-parallel strike. The cliffs are cross-cut by two normal and reverse fault families, with 10–100 m alongshore spacing, with primary and secondary strata-bound fractures perpendicular to the bedding, which combined delimit the cliff rock mass into discrete blocks that are exploited by the erosion process. Erosion, and sometimes plucking, of such beds and blocks on the cliff face was monitored using ground-based structure-from-motion (SfM) photogrammetry, over the course of 5.7 years between 2011 and 2017. To compare with longer time change, cliff-top retreat rate was assessed using SfM-orthorectified archive aerial photographs spanning 1954–2008. We show that the 13,250 m² cliff face released 4500 blocks exceeding 1.45 × 10⁻³ m³, removing a total volume of 170 m³. This equates to an average cliff erosion rate of 3.4 mm/year, which is slightly slower than the 54-year-long local cliff-top retreat (10.8 ± 1.8 mm/year). The vertical distribution of erosion reflects the height of sea water inundation, where the maximum erosion intensity occurs ca. 2 m above high spring-tide water level. Alongshore, the distribution of rockfall scars is concentrated along bed edges bounding cross-cutting faults; the extent of block detachment is controlled by secondary tectonic joints, which may extend through several beds locally sharing similar mechanical strength; and rockfall depth is always a multiple of bed thickness. Over the longer term, we explain block detachment and resultant cliff collapse as a cycle. Erosion nucleates on readily exploitable fractures but elsewhere, the sea only meets defect-free medium-strong to strong rock slabs offering few morphological features for exploitation. Structurally delimited blocks are quarried, and with sufficient time, carve semi-elliptic scars reaching progressively deeper strata to be eroded. Lateral propagation of erosion is directed along mechanical weaknesses in the bedding, and large episodic collapses affect the overhanging slabs via sliding on the weak marl beds. Collapse geometry is confined to one or several triplets of turbidite beds, but never reaches deeper into the cliff than the eroded depth at the foot. We contend that this fracture-limited model of sea-cliff erosion, inferred from the Socoa site dynamics and its peculiar sets of fractures, applies more broadly to other fractured cliff contexts, albeit with site-specific geometries. The initiation of erosion, the propagation of incremental block release and the ultimate full failure of the cliff, have each been shown to be fundamentally directly controlled by structure, which remains a vital control in understanding how cliffed coasts have changed in the past and will change in the future.     

Cosmogenic isotope burial dating of fluvial sediments from the Lower Rhine Embayment,Germany

Cosmogenic isotope burial dating, using ¹⁰Be and ²⁶Al, was applied to Plio–Pleistocene fluvial successions from the Lower Rhine Embayment, Germany. The approach consists of three principal steps: (1) measurement of cosmogenic nuclides in depth profiles, (2) modelling of hypothetical nuclide concentrations based on a first-order conceptualisation of the geological context and the principal succession of depositions and subsequent erosional and burial phases, and (3) using parameter estimation to identify values for the a priori unknown model parameters (burial age, initial nuclide concentrations, terrace erosion rates) that result in minimal disagreement between hypothetical and measured nuclide concentrations.The Late Pliocene Kieseloolite Formation was dated to 3650 ± 1490 ka and the Early Pleistocene Waalre Formation to 900 ± 280 ka. The unconformably overlying Upper Terrace Formation revealed ages of 740 ± 210 ka and 750 ± 250 ka for the two different sites. These findings are in good agreement with independent age control derived by bio-, magneto-, and litho-stratigraphy. Furthermore, identifiability and uncertainty analysis reveal a relationship between burial depth and sensitivity of isotope concentrations to burial age and erosion rate. Our results indicate that using shallower buried samples would enable a considerably more robust estimation of the burial age and the terrace erosion rate. Uncertainties arose mainly from nuclide measurements, and not from the uncertainties derived from modelling or insufficient knowledge of nuclide production and decay properties.     

An un-mixing model to study watershed erosion processes

An un-mixing model is formulated within a Bayesian Markov Chain Monte Carlo framework for use within land-use fingerprinting to study watershed erosion processes. The model has two new components: (1) An equation and erosion process parameter are used to weight tracer signatures from each erosion process within a land-use. (2) An extra tracer distribution and episodic erosion parameter are used to represent soil eroded throughout the sampling duration and thus include the episodic nature of erosion. To test specification of these new parameters, the un-mixing model is applied in the 15 km² Jerome Creek Watershed in the Palouse Region of Northwestern Idaho. Erosion processes include surface erosion upon mountain slopes due to logging in the forest land-use and rill/interrill erosion on cultivated slopes and headcut erosion in riparian floodplains of the agricultural land-use (winter wheat/peas rotation and hay pasture). Episodic erosion occurs for the event where the model is applied. A sensitivity analysis shows that the smallest Bayesian credible set results when the new parameters are specified using hydrologic data and process-based models. The un-mixing model predicts that 90% of the eroded-soil originated from the agricultural land-use and 10% originated from the forest land-use. A comparative study is performed that estimates 90.5% and 9.5% of eroded-soil originated from the agricultural and forest land-uses. Successful performance of the un-mixing model highlights future application as a standalone probabilistic tool to monitor watershed erosion processes that exhibit non-equilibrium conditions and provide calibration data for process-based watershed models.     

Reporting of cosmogenic nuclide data for exposure age and erosion rate determinations

We propose guidelines for the reporting of in situ cosmogenic nuclide data for exposure age and erosion rate determinations. This is motivated by the need to maintain the utility of such data in the future, and to delineate best scientific practice. These guidelines will allow published exposure ages and erosion rates to be recalculated with confidence by others in the future, if, as is likely, procedures to calculate cosmogenic nuclide production rates are modified in the meantime.     

Stochastic erosion of composite banks in alluvial river bends

The erosion of composite river banks is a complex process involving a number of factors including fluvial erosion, seepage erosion, and cantilever mass failure. To predict the rate of bank erosion with these complexities, a stochastic bank erosion model is suitable to define the probability distribution of the controlling variables. In this study, a bank erosion model in a river bend is developed by coupling several bank erosion processes with an existing hydrodynamic and morphological model. The soil erodibility of cohesive bank layers was measured using a submerged jet test apparatus. Seasonal bank erosion rates for four consecutive years at a bend in the Brahmaputra River, India, were measured by repeated bankline surveys. The ability of the model to predict erosion was evaluated in the river bend that displayed active bank erosion. In this study, different monsoon conditions and the distribution functions of two variables were considered in estimating the stochastic bank erosion rate: the probability of the soil erodibility and stochastic stage hydrographs for the nth return period river stage. Additionally, the influences of the deflection angle of the streamflow, longitudinal slope of river channel, and bed material size on bank erosion rate were also investigated. The obtained stochastic erosion predictions were compared with the observed distribution of the annual‐average bank erosion rate of 45 river bends in the Brahmaputra River. The developed model appropriately predicted the short‐term morphological dynamics of sand‐bed river bends with composite banks. Copyright © 2014 John Wiley & Sons, Ltd.     

Morphodynamics of a bedrock‐alluvial meander bend that incises as it migrates outward: approximate solution of permanent form

In meandering rivers cut into bedrock, erosion across a channel cross‐section can be strongly asymmetric. At a meander apex, deep undercutting of the outer bank can result in the formation of a hanging cliff (which may drive hillslope failure), whereas the inner bank adjoins a slip‐off slope that connects to the hillslope itself. Here we propose a physically‐based model for predicting channel planform migration and incision, point bar and slip‐off slope formation, bedrock abrasion, the spatial distribution of alluvial cover, and adaptation of channel width in a mixed bedrock‐alluvial channel. We simplify the analysis by considering a numerical model of steady, uniform bend flow satisfying cyclic boundary conditions. Thus in our analysis, ‘sediment supply’, i.e. the total volume of alluvium in the system, is conserved. In our numerical simulations, the migration rate of the outer bank is a specified parameter. Our simulations demonstrate the existence of an approximate state of dynamic equilibrium corresponding to a near‐solution of permanent form in which a bend of constant curvature, width, cross‐sectional shape and alluvial cover distribution migrates diagonally downward at constant speed, leaving a bedrock equivalent of a point bar on the inside of the bend. Channel width is set internally by the processes of migration and incision. We find that equilibrium width increases with increasing sediment supply, but is insensitive to outer bank migration rate. The slope of the bedrock point bar varies inversely with both outer bank migration rate and sediment supply. Although the migration rate of the outer bank is externally imposed here, we discuss a model modification that would allow lateral side‐wall abrasion to be treated in a manner similar to the process of bedrock incision. Copyright © 2016 John Wiley & Sons, Ltd.     

Determination of both exposure time and denudation rate from an in situ-produced 10Be depth profile: A mathematical proof of uniqueness. Model sensitivity and applications to natural cases

Measurements of radioactive in situ-produced cosmogenic nuclide concentrations in surficial material exposed to cosmic rays allow either determining the long-term denudation rate assuming that the surface studied has reached steady-state (where production and losses by denudation and radioactive decay are in equilibrium) (infinite exposure time), or dating the initiation of exposure to cosmic rays, assuming that the denudation and post-depositional processes are negligible. Criteria for determining whether a surface is eroding or undergoing burial as well as quantitative information on denudation or burial rates may be obtained from cosmogenic nuclide depth profiles. With the refinement of the physical parameters involved in the production of in situ-produced cosmogenic nuclides, a unique well-constrained depth profile now permits determination of both the exposure time and the denudation rate affecting a surface. In this paper, we first mathematically demonstrate that the exponential decrease of the in situ-produced ¹⁰Be concentrations observed along a depth profile constrains a unique exposure time and denudation rate when considering both neutrons and muons. In the second part, an improved chi-square inversion model is described and tested in the third part with actual measured profiles.     

Geological calibration of spallation production rates in the CRONUS-Earth project

Models of the production of cosmogenic nuclides typically incorporate an adjustable production rate parameter that is scaled for variations in production with latitude and altitude. In practice, this production rate parameter is set by calibration of the model using cosmogenic nuclide data from sites with independent age constraints. In this paper, we describe a calibration procedure developed during the Cosmic-Ray Produced Nuclide Systematics on Earth (CRONUS-Earth) project and its application to an extensive data set that included both new CRONUS-Earth samples and samples from previously published studies. We considered seven frameworks for elevation and latitude scaling and five commonly used cosmogenic nuclides, ³He, ¹⁰Be, ¹⁴C, ²⁶Al, and ³⁶Cl. In general, the results show that the calibrated production rates fail statistical tests of goodness-of-fit. One conclusion from the calibration results is that two newly developed scaling frameworks and the widely used Lal scaling framework provide qualitatively similar fits to the data, while neutron-monitor based scaling frameworks have much poorer fit to the data. To further test the fitted models, we computed site ages for a number of secondary sites not included in the primary calibration data set. The root-mean-square percent differences between the median computed ages for these secondary sites and independent ages range from 7.1% to 27.1%, differences that are much larger than the typical uncertainties in the site ages. The results indicate that there are substantial unresolved difficulties in modeling cosmogenic nuclide production and the calibration of production rates.     

Landscape form and millennial erosion rates in the San Gabriel Mountains,CA

It has been long hypothesized that topography, as well as climate and rock strength, exert first order controls on erosion rates. Here we use detrital cosmogenic ¹⁰Be from 50 basins, ranging in size from 1 to 150 km², to measure millennial erosion rates across the San Gabriel Mountains in southern California, where a strong E–W gradient in relief compared to weak variation in precipitation and lithology allow us to isolate the relationship between topographic form and erosion rate. Our erosion rates range from 35 to 1100 m/Ma, and generally agree with both decadal sediment fluxes and long term exhumation rates inferred from low temperature thermochronometry. Catchment-mean hillslope angle increases with erosion rate until ～ 300 m/Ma, at which point slopes become invariant with erosion rate. Although this sort of relation has been offered as support for non-linear models of soil transport, we use 1-D analytical hillslope profiles derived from existing soil transport laws to show that a model with soil flux linear in slope, but including a slope stability threshold, is indistinguishable from a non-linear law within the scatter of our data. Catchment-mean normalized channel steepness index increases monotonically, though non-linearly, with erosion rate throughout the San Gabriel Mountains, even where catchment-mean hillslope angles have reached a threshold. This non-linearity can be mostly accounted for by a stochastic threshold incision model, though additional factors likely contribute to the observed relationship between channel steepness and erosion rate. These findings substantiate the claim that the normalized channel steepness index is an important topographic metric in active ranges.     

地球空间稳定核素的趋势分析方程与物质的超光速运动规律

通过讨论地球空间已有稳定核素内质子数与中子数的分布趋势,介绍了稳定核素的趋势分析方法及其有关周期性分布方程形式,给出了理论方程曲线与地球空间稳定核素实验数据分布点的对比结果,进而给出了稳定核素极限值和元素周期表中化学元素极限,以及其与正负粒子对的可能对应关系方程,包括位于电子中微子层面附近的粒子质量量级初步估计.随后通过建立真空物质能量状态的二个假设,及基于等效Binet方程,给出了与Einstein狭义相对论有关结论相融合的物质粒子以光速及超光速运动的质量及能量方程;作为推论,对这些方程与暗物质及暗能量的可能对应关系予以了初步探讨.     

Regional‐scale sedimentation process models from airborne gamma ray remote sensing and digital elevation data

Airborne gamma ray survey data were used to provide information on potassium, thorium and uranium concentrations in surface soil and rock in arid central Australia. Spatial patterns in these radioelements allow tracing of paths of sediment at catchment scale. Survey elevation data are combined with contour data to produce digital elevation models for terrain analysis, tracing of sediment flow paths and modelling of extreme floods. Gamma ray data show consistent variation with slope, a limited range of drainage areas, and erosion/deposition models derived from the conservation of mass equation. Supply‐limited sediment transport models give a reasonable reproduction of observed radioelement distribution but some elements of the distribution pattern reflect the area inundated by 500–1000 year floods rather than the effects of simple downslope movement. Partial area sediment supply models are derived by downstream accumulation of erosion and deposition rates calculated using the conservation of mass equation with transport laws based on slope alone and stream power. Comparison with observed radioelement patterns suggests that both transport laws apply in different parts of the landscape. Regional‐scale sediment transport models will require a range of models depending on location in the landscape and event frequency. This approach may allow estimation of sediment delivery ratios. Copyright © 2001 John Wiley & Sons, Ltd.     

Cosmogenic 21Ne analysis of individual detrital grains: Opportunities and limitations

We use a numerical model describing cosmogenic nuclide acquisition in sediment moving through the upper Gaub River catchment to evaluate the extent to which aspects of source area geomorphology and geomorphological processes can be inferred from frequency distributions of cosmogenic ²¹Ne (²¹Ne_c) concentrations in individual detrital grains. The numerical model predicts the pathways of sediment grains from their source to the outlet of the catchment and calculates the total ²¹Ne_c concentration that each grain acquires along its pathway. The model fully accounts for variations in nuclide production due to changes in latitude, altitude and topographic shielding and allows for spatially variable erosion and sediment transport rates. Model results show that the form of the frequency distribution of ²¹Ne_c concentrations in exported sediment is sensitive to the range and spatial distribution of processes operating in the sediment's source areas and that this distribution can be used to infer the range and spatial distribution of erosion rates that characterise the catchment. The results also show that lithology can affect the form of the ²¹Ne_c concentration distribution indirectly by exerting control on the spatial pattern of denudation in a catchment. Model results further indicate that the form of the distribution of ²¹Ne_c concentrations in the exported sediment can also be affected by the acquisition of ²¹Ne_c after detachment from bedrock, in the diffusive (hillslope) and/or advective (fluvial) domains. However, for such post‐detachment nuclide acquisition to be important, this effect needs to at least equal the nuclide acquisition prior to detachment from bedrock. Copyright © 2009 John Wiley and Sons, Ltd.     

Cosmogenic surface exposure dating the last deglaciation in Denmark: Discrepancies with independent age constraints suggest delayed periglacial landform stabilisation

Cosmogenic nuclide surface exposure ages are determined from in situ ¹⁰Be and ³⁶Cl analysis of 38 rock surfaces found in different glacial landforms in Denmark. Dating of erratic boulders and adjacent ice-sculpted bedrock on the island of Bornholm in the western Baltic Sea reveals almost identical values. This suggests that little if any inherited nuclides are present in the sampled boulders. West of the Last Glacial Maximum (LGM) ice margin in Denmark ages reflect exposure from the Middle Weichselian. East of the LGM margin exposure ages from 35 samples show Late Weichselian ages in a range between 20.6–11.9 ka. To test to what extent these dates reflect the onset of deglaciation immediately after cessation of active glacier flow, surface exposure ages are evaluated against independent chronologies of Late Weichselian ice-sheet fluctuations in southwestern Scandinavia. The Bornholm dates agree with the independent age model, however, in the data set for eastern Denmark only less than half the surface exposure ages lie within the expected age envelope. This apparent mismatch is most likely due to post-glaciation shielding and delayed surface stabilisation compared to the timing of ice-margin retreat. Thus ages from boulders resting in dead-ice moraines and mass wasting landscapes underestimate deglaciation by 3–6 thousand years. The results quantify the impact of exhumation and landform stabilisation on cosmogenic surface exposure ages on millennial scales. We conclude, that interpretation of cosmogenic exposure ages should include careful evaluation of possible post-depositional landform transformation in attempts to fine tune ages of e.g. end moraine features. With reference to independent age models we critically evaluate glacier advance – retreat scenarios from regions around the southern Baltic that alone are based on weighted average ages of cosmogenic exposure dating.     

Assessment of calanchi and rill–interrill erosion susceptibilities using terrain analysis and geostochastics: A case study in the Oltrepo Pavese,Northern Apennines,Italy

Soil erosion is one of the most important environmental problems distributed worldwide. In the last decades, numerous studies have been published on the assessment of soil erosion and the related processes and forms using empirical, conceptual and physically based models. For the prediction of the spatial distribution, more and more sophisticated stochastic modelling approaches have been proposed – especially on smaller spatial scales such as river basins. In this work, we apply a maximum entropy model (MaxEnt) to evaluate badlands (calanchi) and rill–interrill (sheet erosion) areas in the Oltrepo Pavese (Northern Apennines, Italy). The aim of the work is to assess the important environmental predictors that influence calanchi and rill–interrill erosion at the regional scale. We used 13 topographic parameters derived from a 12 m digital elevation model (TanDEM-X) and data on the lithology and land use. Additional information about the vegetation is introduced through the normalized difference vegetation index based on remotely sensed data (ASTER images). The results are presented in the form of susceptibility maps showing the spatial distribution of the occurrence probability for calanchi and rill–interrill erosion. For the validation of the MaxEnt model results, a support vector machine approach was applied. The models show reliable results and highlight several locations of the study area that are potentially prone to future soil erosion. Thus, coping and mitigation strategies may be developed to prevent or fight the soil erosion phenomenon under consideration. © 2020 John Wiley & Sons, Ltd.     

Chronology analysis of huge landslide based on ESR dating materials on sliding face in carbonate areas of south eastern Tibet

Determining occurrence age is key to chronology analysis of huge landslides. However, in carbonated bedrock areas, it is difficult to carry out geochronological surveys of landslide evolution due to lack of suitable dating materials and methods. In carbonated bedrock areas, the cosmic nuclides 36Cl and 14C dating methods have been used for obtaining the exposed age of bedrock. But intense weathering, erosion and dissolution make it difficult to obtain an accurate total amount of cosmogenic nuclide accumulation, resulting in the dating results being affected greatly by weathering, erosion and dissolution of carbonate rock. Therefore, it is necessary to find new dating materials or methods for determining the occurrence age of landslide in carbonate bedrock areas with little quaternary sediments. In this study, a huge landslide developed in the carbonate bedrock area of south eastern Tibet is introduced briefly, which is named the Qiaojia landslide and blocked the Jinsha River. A thin layer including three obvious recrystallized carbonate sub-layers was formed in the sliding zone, indicating different sliding events. Three recrystallized carbonate samples were then collected for ESR dating research to obtain the age information of the Qiaojia landslide: one occurred 79 ± 8 ka years ago, and the other occurred 60 ± 10 ka years ago, and the third occurred 25 ± 2 ka years ago, which are in general agreement with other studies in the surrounding of the study area. It suggests that recrystallized carbonate on sliding face is suitable for chronologically identifying landslide activity by using ESR method in the carbonate area, especially in the areas with little or no quaternary sediments.