Pub Date : 2024-04-22DOI: 10.5194/esurf-12-537-2024
Peter Arlinghaus, Corinna Schrum, Ingrid Kröncke, Wenyan Zhang
Abstract. Benthos has long been recognized as an important factor influencing local sediment stability, deposition, and erosion rates. However, its role in long-term (annual to decadal scale) and large-scale coastal morphological change remains largely speculative. This study aims to derive a quantitative understanding of the importance of benthos in the morphological development of a tidal embayment (Jade Bay) as representative of tidal coastal regions. To achieve this, we first applied a machine-learning-aided species abundance model to derive a complete map of benthos (functional groups, abundance, and biomass) in the study area, based on abundance and biomass measurements. The derived data were used to parameterize the benthos effect on sediment stability, erosion rates and deposition rates, erosion and hydrodynamics in a 3-dimensional hydro-eco-morphodynamic model, which was then applied to Jade Bay to hindcast the morphological and sediment change for 2000–2009. Simulation results indicate significantly improved performance with the benthos effect included. Simulations including benthos show consistency with measurements regarding morphological and sediment changes, while abiotic drivers (tides, storm surges) alone result in a reversed pattern in terms of erosion and deposition contrary to measurement. Based on comparisons among scenarios with various combinations of abiotic and biotic factors, we further investigated the level of complexity of the hydro-eco-morphodynamic models that is needed to capture long-term and large-scale coastal morphological development. The accuracy in the parameterization data was crucial for increasing model complexity. When the parameterization uncertainties were high, the increased model complexity decreased the model performance.
{"title":"Benthos as a key driver of morphological change in coastal regions","authors":"Peter Arlinghaus, Corinna Schrum, Ingrid Kröncke, Wenyan Zhang","doi":"10.5194/esurf-12-537-2024","DOIUrl":"https://doi.org/10.5194/esurf-12-537-2024","url":null,"abstract":"Abstract. Benthos has long been recognized as an important factor influencing local sediment stability, deposition, and erosion rates. However, its role in long-term (annual to decadal scale) and large-scale coastal morphological change remains largely speculative. This study aims to derive a quantitative understanding of the importance of benthos in the morphological development of a tidal embayment (Jade Bay) as representative of tidal coastal regions. To achieve this, we first applied a machine-learning-aided species abundance model to derive a complete map of benthos (functional groups, abundance, and biomass) in the study area, based on abundance and biomass measurements. The derived data were used to parameterize the benthos effect on sediment stability, erosion rates and deposition rates, erosion and hydrodynamics in a 3-dimensional hydro-eco-morphodynamic model, which was then applied to Jade Bay to hindcast the morphological and sediment change for 2000–2009. Simulation results indicate significantly improved performance with the benthos effect included. Simulations including benthos show consistency with measurements regarding morphological and sediment changes, while abiotic drivers (tides, storm surges) alone result in a reversed pattern in terms of erosion and deposition contrary to measurement. Based on comparisons among scenarios with various combinations of abiotic and biotic factors, we further investigated the level of complexity of the hydro-eco-morphodynamic models that is needed to capture long-term and large-scale coastal morphological development. The accuracy in the parameterization data was crucial for increasing model complexity. When the parameterization uncertainties were high, the increased model complexity decreased the model performance.","PeriodicalId":48749,"journal":{"name":"Earth Surface Dynamics","volume":"236 1","pages":""},"PeriodicalIF":3.4,"publicationDate":"2024-04-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140635765","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-04-22DOI: 10.5194/egusphere-2024-1138
Alexander B. Prescott, Jon D. Pelletier, Satya Chataut, Sriram Ananthanarayan
Abstract. Calculating contributing area (often used as a proxy for surface water discharge) within a Digital Elevation Model (DEM) or Landscape Evolution Model (LEM) is a fundamental operation in geomorphology. Here we document that a commonly used multiple-flow-direction algorithm for calculating contributing area, i.e., D∞ of Tarboton (1997), is sufficiently biased along the cardinal and ordinal directions that it is unsuitable for some standard applications of flow-routing algorithms. We revisit the purported excess dispersion of the MFD algorithm of Freeman (1991) that motivated the development of D∞ and demonstrate that MFD is superior to D∞ when tested against analytic solutions for the contributing areas of idealized landforms and the predictions of the shallow-water-equation solver FLO-2D for more complex landforms in which the water-surface slope is closely approximated by the bed slope. We also introduce a new flow-routing algorithm entitled IDS (in reference to the iterative depth-and-slope-dependent nature of the algorithm) that is more suitable than MFD for applications in which the bed and water-surface slopes differ substantially. IDS solves for water flow depths under steady hydrologic conditions by distributing the discharge delivered to each grid point from upslope to its downslope neighbors in rank order of elevation (highest to lowest) and in proportion to a power-law function of the square root of the water-surface slope and the five-thirds power of the water depth, mimicking the relationships among water discharge, depth, and surface slope in Manning’s equation. IDS is iterative in two ways: 1) water depths are added in small increments so that the water-surface slope can gradually differ from the bed slope, facilitating the spreading of water in areas of laterally unconfined flow, and 2) the partitioning of discharge from high to low elevations can be repeated, improving the accuracy of the solution as the water depths of downslope grid points become more well approximated with each successive iteration. We assess the performance of IDS by comparing its results to those of FLO-2D for a variety of real and idealized landforms and to an analytic solution of the shallow-water equations. We also demonstrate how IDS can be modified to solve other fluid-dynamical nonlinear partial differential equations arising in Earth-surface processes, such as the Boussinesq equation for the height of the water table in an unconfined aquifer.
{"title":"An evaluation of flow-routing algorithms for calculating contributing area on regular grids","authors":"Alexander B. Prescott, Jon D. Pelletier, Satya Chataut, Sriram Ananthanarayan","doi":"10.5194/egusphere-2024-1138","DOIUrl":"https://doi.org/10.5194/egusphere-2024-1138","url":null,"abstract":"<strong>Abstract.</strong> Calculating contributing area (often used as a proxy for surface water discharge) within a Digital Elevation Model (DEM) or Landscape Evolution Model (LEM) is a fundamental operation in geomorphology. Here we document that a commonly used multiple-flow-direction algorithm for calculating contributing area, i.e., D∞ of Tarboton (1997), is sufficiently biased along the cardinal and ordinal directions that it is unsuitable for some standard applications of flow-routing algorithms. We revisit the purported excess dispersion of the MFD algorithm of Freeman (1991) that motivated the development of D∞ and demonstrate that MFD is superior to D∞ when tested against analytic solutions for the contributing areas of idealized landforms and the predictions of the shallow-water-equation solver FLO-2D for more complex landforms in which the water-surface slope is closely approximated by the bed slope. We also introduce a new flow-routing algorithm entitled IDS (in reference to the iterative depth-and-slope-dependent nature of the algorithm) that is more suitable than MFD for applications in which the bed and water-surface slopes differ substantially. IDS solves for water flow depths under steady hydrologic conditions by distributing the discharge delivered to each grid point from upslope to its downslope neighbors in rank order of elevation (highest to lowest) and in proportion to a power-law function of the square root of the water-surface slope and the five-thirds power of the water depth, mimicking the relationships among water discharge, depth, and surface slope in Manning’s equation. IDS is iterative in two ways: 1) water depths are added in small increments so that the water-surface slope can gradually differ from the bed slope, facilitating the spreading of water in areas of laterally unconfined flow, and 2) the partitioning of discharge from high to low elevations can be repeated, improving the accuracy of the solution as the water depths of downslope grid points become more well approximated with each successive iteration. We assess the performance of IDS by comparing its results to those of FLO-2D for a variety of real and idealized landforms and to an analytic solution of the shallow-water equations. We also demonstrate how IDS can be modified to solve other fluid-dynamical nonlinear partial differential equations arising in Earth-surface processes, such as the Boussinesq equation for the height of the water table in an unconfined aquifer.","PeriodicalId":48749,"journal":{"name":"Earth Surface Dynamics","volume":"115 1","pages":""},"PeriodicalIF":3.4,"publicationDate":"2024-04-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140635134","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-04-18DOI: 10.5194/egusphere-2024-1036
W. Marijn van der Meij
Abstract. Soil erosion is one of the main threats to agricultural food production due to the loss of fertile soil. Determination of erosion rates is essential to quantify the degree of land degradation, but it is inherently challenging to determine temporally dynamic erosion rates over agricultural time scales. Optically Stimulated Luminescence (OSL) dating can provide temporally-resolved deposition rates by determining the last moment of daylight exposure of buried colluvial deposits. However, these deposition rates may differ substantially from the actual hillslope erosion rates. In this study, OSL-based deposition rates were converted to hillslope erosion rates through inverse modelling with soil-landscape evolution model ChronoLorica. This model integrates geochronological tracers into the simulations of soil mixing and redistribution. The model was applied to a kettle hole catchment in north-eastern Germany, which has been affected by tillage erosion over the last 5000 years. The initial shape of the landscape and the land use history are well-constrained, enabling accurate simulations of the landscape evolution that incorporate uncertainties in the model inputs. The calibrated model reveals an increase in erosion rates of almost to orders of magnitude from pre-historic ard ploughing up to recent intensive land management. The simulated rates match well with independent age controls from the same catchment. Uncertainty in the reconstructed initial landscape and land use histories had a minor influence of 12–16 % on the simulated rates. The simulations showed that the deposition rates were on average 1.5 higher than the erosion rates due to the ratio of erosional and depositional area. Recent artificial drainage and land reclamation have increased deposition rates up to five times the erosion rates, emphasizing the need of cautious interpretation of deposition rates as erosion proxies. This study demonstrates the suitability of ChronoLorica for upscaling experimental geochronological data to better understand landscape evolution in agricultural settings.
{"title":"Translating deposition rates into erosion rates with landscape evolution modelling","authors":"W. Marijn van der Meij","doi":"10.5194/egusphere-2024-1036","DOIUrl":"https://doi.org/10.5194/egusphere-2024-1036","url":null,"abstract":"<strong>Abstract.</strong> Soil erosion is one of the main threats to agricultural food production due to the loss of fertile soil. Determination of erosion rates is essential to quantify the degree of land degradation, but it is inherently challenging to determine temporally dynamic erosion rates over agricultural time scales. Optically Stimulated Luminescence (OSL) dating can provide temporally-resolved deposition rates by determining the last moment of daylight exposure of buried colluvial deposits. However, these deposition rates may differ substantially from the actual hillslope erosion rates. In this study, OSL-based deposition rates were converted to hillslope erosion rates through inverse modelling with soil-landscape evolution model ChronoLorica. This model integrates geochronological tracers into the simulations of soil mixing and redistribution. The model was applied to a kettle hole catchment in north-eastern Germany, which has been affected by tillage erosion over the last 5000 years. The initial shape of the landscape and the land use history are well-constrained, enabling accurate simulations of the landscape evolution that incorporate uncertainties in the model inputs. The calibrated model reveals an increase in erosion rates of almost to orders of magnitude from pre-historic ard ploughing up to recent intensive land management. The simulated rates match well with independent age controls from the same catchment. Uncertainty in the reconstructed initial landscape and land use histories had a minor influence of 12–16 % on the simulated rates. The simulations showed that the deposition rates were on average 1.5 higher than the erosion rates due to the ratio of erosional and depositional area. Recent artificial drainage and land reclamation have increased deposition rates up to five times the erosion rates, emphasizing the need of cautious interpretation of deposition rates as erosion proxies. This study demonstrates the suitability of ChronoLorica for upscaling experimental geochronological data to better understand landscape evolution in agricultural settings.","PeriodicalId":48749,"journal":{"name":"Earth Surface Dynamics","volume":"14 1","pages":""},"PeriodicalIF":3.4,"publicationDate":"2024-04-18","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140629824","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-04-17DOI: 10.5194/esurf-12-515-2024
Prakash Pokhrel, Mikael Attal, Hugh D. Sinclair, Simon M. Mudd, Mark Naylor
Abstract. Sediment grains are progressively rounded during their transport down a river. For more than a century, Earth scientists have used the roundness of pebbles within modern sediment, and of clasts within conglomerates, as a key metric to constrain the sediment's transport history and source area(s). However, the current practices of assessment of pebble roundness are mainly qualitative and based on time-consuming manual measurement methods. This qualitative judgement provides the transport history only in a broad sense, such as classifying distance as “near” or “far”. In this study, we propose a new model that quantifies the relationship between roundness and the transport distance. We demonstrate that this model can be applied to the clasts of multiple lithologies including modern sediment, as well as conglomerates, deposited by ancient river systems. We present field data from two Himalayan catchments in Nepal. We use the normalized isoperimetric ratio (IRn), which relates a pebble's area (A) to its perimeter (P), to quantify roundness. The maximum analytical value for IRn is 1, and IRn is expected to increase with transport distance. We propose a non-linear roundness model based on our field data, whereby the difference between a grain's IRn and the maximum value of 1 decays exponentially with transport distance, mirroring Sternberg's model of mass loss or size reduction by abrasion. This roundness model predicts an asymptotic behaviour for IRn, and the distance over which IRn approaches the asymptote is controlled by a rounding coefficient. Our field data suggest that the roundness coefficient for granite pebbles is 9 times that of quartzite pebbles. Using this model, we reconstruct the transport history of a Pliocene paleo-river deposit preserved at the base of the Kathmandu intermontane basin. These results, along with other sedimentary evidence, imply that the paleo-river was much longer than the length of the Kathmandu Basin and that it must have lost its headwaters through drainage capture. We further explore the extreme rounding of clasts from Miocene conglomerate of the Siwalik zone and find evidence of sediment recycling.
{"title":"Downstream rounding rate of pebbles in the Himalaya","authors":"Prakash Pokhrel, Mikael Attal, Hugh D. Sinclair, Simon M. Mudd, Mark Naylor","doi":"10.5194/esurf-12-515-2024","DOIUrl":"https://doi.org/10.5194/esurf-12-515-2024","url":null,"abstract":"Abstract. Sediment grains are progressively rounded during their transport down a river. For more than a century, Earth scientists have used the roundness of pebbles within modern sediment, and of clasts within conglomerates, as a key metric to constrain the sediment's transport history and source area(s). However, the current practices of assessment of pebble roundness are mainly qualitative and based on time-consuming manual measurement methods. This qualitative judgement provides the transport history only in a broad sense, such as classifying distance as “near” or “far”. In this study, we propose a new model that quantifies the relationship between roundness and the transport distance. We demonstrate that this model can be applied to the clasts of multiple lithologies including modern sediment, as well as conglomerates, deposited by ancient river systems. We present field data from two Himalayan catchments in Nepal. We use the normalized isoperimetric ratio (IRn), which relates a pebble's area (A) to its perimeter (P), to quantify roundness. The maximum analytical value for IRn is 1, and IRn is expected to increase with transport distance. We propose a non-linear roundness model based on our field data, whereby the difference between a grain's IRn and the maximum value of 1 decays exponentially with transport distance, mirroring Sternberg's model of mass loss or size reduction by abrasion. This roundness model predicts an asymptotic behaviour for IRn, and the distance over which IRn approaches the asymptote is controlled by a rounding coefficient. Our field data suggest that the roundness coefficient for granite pebbles is 9 times that of quartzite pebbles. Using this model, we reconstruct the transport history of a Pliocene paleo-river deposit preserved at the base of the Kathmandu intermontane basin. These results, along with other sedimentary evidence, imply that the paleo-river was much longer than the length of the Kathmandu Basin and that it must have lost its headwaters through drainage capture. We further explore the extreme rounding of clasts from Miocene conglomerate of the Siwalik zone and find evidence of sediment recycling.","PeriodicalId":48749,"journal":{"name":"Earth Surface Dynamics","volume":"12 1","pages":""},"PeriodicalIF":3.4,"publicationDate":"2024-04-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140613606","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-04-16DOI: 10.5194/egusphere-2024-792
Janbert Aarnink, Tom Beucler, Marceline Vuaridel, Virginia Ruiz-Villanueva
Abstract. Instream large wood (i.e., downed trees, branches and roots larger than 1 m in length and 10 cm diameter) has essential geopmorphological and ecological functions supporting the health of river ecosystems. Still, even though its transport during floods may pose a risk, it is rarely observed and, therefore, poorly understood. This paper presents a novel approach to detect pieces of instream wood from video. The approach uses a Convolutional Neural Network to detect wood automatically. We sampled data to represent different wood transport conditions, combining 20 datasets to yield thousands of instream wood images. We designed multiple scenarios using different data subsets with and without data augmentation and analyzed the contribution of each one to the effectiveness of the model using k-fold cross-validation. The mean average precision of the model varies between 35 and 93 percent, and is highly influenced by the quality of the data which it detects. When the image resolution is low, the identified components in the labeled pieces, rather than exhibiting distinct characteristics such as bark or branches, appear more akin to amorphous masses or 'blobs'. We found that the model detects wood with a mean average precision of 67 percent when using a 418 pixels input image resolution. Also, improvements of up to 23 percent could be achieved in some instances and increasing the input resolution raised the weighted mean average precision to 74 percent. We show that the detection performance on a specific dataset is not solely determined by the complexity of the network or the training data. Therefore, the findings of this paper can be used when designing a custom wood detection network. With the growing availability of flood-related videos featuring wood uploaded to the internet, this methodology facilitates the quantification of wood transport across a wide variety of data sources.
{"title":"Automatic detection of instream large wood in videos using deep learning","authors":"Janbert Aarnink, Tom Beucler, Marceline Vuaridel, Virginia Ruiz-Villanueva","doi":"10.5194/egusphere-2024-792","DOIUrl":"https://doi.org/10.5194/egusphere-2024-792","url":null,"abstract":"<strong>Abstract.</strong> Instream large wood (i.e., downed trees, branches and roots larger than 1 m in length and 10 cm diameter) has essential geopmorphological and ecological functions supporting the health of river ecosystems. Still, even though its transport during floods may pose a risk, it is rarely observed and, therefore, poorly understood. This paper presents a novel approach to detect pieces of instream wood from video. The approach uses a Convolutional Neural Network to detect wood automatically. We sampled data to represent different wood transport conditions, combining 20 datasets to yield thousands of instream wood images. We designed multiple scenarios using different data subsets with and without data augmentation and analyzed the contribution of each one to the effectiveness of the model using k-fold cross-validation. The mean average precision of the model varies between 35 and 93 percent, and is highly influenced by the quality of the data which it detects. When the image resolution is low, the identified components in the labeled pieces, rather than exhibiting distinct characteristics such as bark or branches, appear more akin to amorphous masses or 'blobs'. We found that the model detects wood with a mean average precision of 67 percent when using a 418 pixels input image resolution. Also, improvements of up to 23 percent could be achieved in some instances and increasing the input resolution raised the weighted mean average precision to 74 percent. We show that the detection performance on a specific dataset is not solely determined by the complexity of the network or the training data. Therefore, the findings of this paper can be used when designing a custom wood detection network. With the growing availability of flood-related videos featuring wood uploaded to the internet, this methodology facilitates the quantification of wood transport across a wide variety of data sources.","PeriodicalId":48749,"journal":{"name":"Earth Surface Dynamics","volume":"33 1","pages":""},"PeriodicalIF":3.4,"publicationDate":"2024-04-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140560697","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-04-15DOI: 10.5194/egusphere-2024-592
Aline Zinelabedin, Joel Mohren, Maria Wierzbicka-Wieczorek, Tibor Janos Dunai, Stefan Heinze, Benedikt Ritter
Abstract. While the formation of periglacial wedges and polygonal patterned grounds has been extensively studied and many of the processes involved have been understood, knowledge on the formation of similar post-sedimentary features found in arid to hyperarid environments remains largely rudimentary. Our contribution to fill this gap is the investigation of a network of vertically laminated calcium sulphate-rich wedges in the subsurface of the Aroma fan in the northern Atacama Desert. The subsurface wedges are characterised by high anhydrite contents and hence differ from the wedge and polygon structures of other study sites in the Atacama Desert, which appear to have been predominantly formed by thermal contraction or desiccation processes. By contrast, haloturbation mechanisms are thought to be a main driver of wedge formation at the Aroma fan site. Haloturbation requires moisture input, and hence Aroma fan wedge formation is likely to be associated with meteoric water received from sporadic rain events and episodes of arid climate characterised by slightly wetter conditions than prevailing at present. The polygonal patterned ground is covered by a stratigraphically younger gypsum-dominated surface crust cover. The presence of the surface crust could indicate an environmental change towards drier conditions, which favoured surface accumulation of calcium sulphate and other salts by means of atmospheric deposition. Such a climatic shift could have caused a deceleration of haloturbation and other wedge formation processes in the subsurface, although modern sediment conveyance from the surface towards its interior still appears to occur along cracks within the crust. In order to gain comprehensive insights into the complex mechanisms involved in wedge formation and formation rates, the establishment of a geochronological framework directly obtained from wedge and crust material remains indispensable. The temporal resolution of wedge growth stored within the succession of vertical laminae promises a high potential for the calcium sulphate wedges to be used as palaeoclimate archives, potentially helping to unravel wedge and polygonal patterned ground formation in other water-limited environments, such as Mars.
{"title":"Haloturbation in the northern Atacama Desert revealed by a hidden subsurface network of calcium sulphate wedges","authors":"Aline Zinelabedin, Joel Mohren, Maria Wierzbicka-Wieczorek, Tibor Janos Dunai, Stefan Heinze, Benedikt Ritter","doi":"10.5194/egusphere-2024-592","DOIUrl":"https://doi.org/10.5194/egusphere-2024-592","url":null,"abstract":"<strong>Abstract.</strong> While the formation of periglacial wedges and polygonal patterned grounds has been extensively studied and many of the processes involved have been understood, knowledge on the formation of similar post-sedimentary features found in arid to hyperarid environments remains largely rudimentary. Our contribution to fill this gap is the investigation of a network of vertically laminated calcium sulphate-rich wedges in the subsurface of the Aroma fan in the northern Atacama Desert. The subsurface wedges are characterised by high anhydrite contents and hence differ from the wedge and polygon structures of other study sites in the Atacama Desert, which appear to have been predominantly formed by thermal contraction or desiccation processes. By contrast, haloturbation mechanisms are thought to be a main driver of wedge formation at the Aroma fan site. Haloturbation requires moisture input, and hence Aroma fan wedge formation is likely to be associated with meteoric water received from sporadic rain events and episodes of arid climate characterised by slightly wetter conditions than prevailing at present. The polygonal patterned ground is covered by a stratigraphically younger gypsum-dominated surface crust cover. The presence of the surface crust could indicate an environmental change towards drier conditions, which favoured surface accumulation of calcium sulphate and other salts by means of atmospheric deposition. Such a climatic shift could have caused a deceleration of haloturbation and other wedge formation processes in the subsurface, although modern sediment conveyance from the surface towards its interior still appears to occur along cracks within the crust. In order to gain comprehensive insights into the complex mechanisms involved in wedge formation and formation rates, the establishment of a geochronological framework directly obtained from wedge and crust material remains indispensable. The temporal resolution of wedge growth stored within the succession of vertical laminae promises a high potential for the calcium sulphate wedges to be used as palaeoclimate archives, potentially helping to unravel wedge and polygonal patterned ground formation in other water-limited environments, such as Mars.","PeriodicalId":48749,"journal":{"name":"Earth Surface Dynamics","volume":"58 1","pages":""},"PeriodicalIF":3.4,"publicationDate":"2024-04-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140560524","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-04-10DOI: 10.5194/egusphere-2023-2694
Hayden L. Jacobson, Danica L. Roth, Gabriel Walton, Margaret Zimmer, Kerri Johnson
Abstract. Post-fire changes to the transport regime of dry ravel, which describes the transport of individual particles downslope, are poorly constrained on a regional level but critical to understand as ravel may contribute to elevated sediment fluxes and associated debris-flow activity observed post-fire in the western United States. In this study, we evaluated post-fire variability in dry ravel travel distance exceedance probabilities and disentrainment rates through a series of field experiments simulating ravel with particles collected in situ. We conducted experiments between March 2021 and March 2022 on soil-mantled hillslopes in the Diablo Range of central coastal California following the Santa Clara Unit Lightning Complex fire of August 2020 with the goal of identifying a regime of “bounded” (light-tailed) or “runaway” (heavy-tailed or nonlocal) motion for different particle sizes between 3 and 35 mm. We conducted this study on both grassy south-facing slopes and oak woodland north-facing slopes. We tracked the post-fire evolution of particle transport regimes by fitting a probabilistic Lomax distribution model to the empirical travel distance exceedance probabilities of different particle sizes on a range of experimental slopes. Our experimental results indicated that a general transition from more runaway to more bounded transport occurred for our largest experimental particles (median intermediate axis of 28 mm) on south-facing slopes as vegetation recovered within the first year post-fire, while small and medium particles (median intermediate axes of 6 and 13 mm respectively) on south- or north-facing slopes and large particles on north-facing slopes did not experience notable changes in transport behavior. After the first year, seasonal variation in vegetation characteristics, such as grass density, appeared to control particle motion.
{"title":"Post-fire Variability in Sediment Transport by Ravel in the Diablo Range","authors":"Hayden L. Jacobson, Danica L. Roth, Gabriel Walton, Margaret Zimmer, Kerri Johnson","doi":"10.5194/egusphere-2023-2694","DOIUrl":"https://doi.org/10.5194/egusphere-2023-2694","url":null,"abstract":"<strong>Abstract.</strong> Post-fire changes to the transport regime of dry ravel, which describes the transport of individual particles downslope, are poorly constrained on a regional level but critical to understand as ravel may contribute to elevated sediment fluxes and associated debris-flow activity observed post-fire in the western United States. In this study, we evaluated post-fire variability in dry ravel travel distance exceedance probabilities and disentrainment rates through a series of field experiments simulating ravel with particles collected in situ. We conducted experiments between March 2021 and March 2022 on soil-mantled hillslopes in the Diablo Range of central coastal California following the Santa Clara Unit Lightning Complex fire of August 2020 with the goal of identifying a regime of “bounded” (light-tailed) or “runaway” (heavy-tailed or nonlocal) motion for different particle sizes between 3 and 35 mm. We conducted this study on both grassy south-facing slopes and oak woodland north-facing slopes. We tracked the post-fire evolution of particle transport regimes by fitting a probabilistic Lomax distribution model to the empirical travel distance exceedance probabilities of different particle sizes on a range of experimental slopes. Our experimental results indicated that a general transition from more runaway to more bounded transport occurred for our largest experimental particles (median intermediate axis of 28 mm) on south-facing slopes as vegetation recovered within the first year post-fire, while small and medium particles (median intermediate axes of 6 and 13 mm respectively) on south- or north-facing slopes and large particles on north-facing slopes did not experience notable changes in transport behavior. After the first year, seasonal variation in vegetation characteristics, such as grass density, appeared to control particle motion.","PeriodicalId":48749,"journal":{"name":"Earth Surface Dynamics","volume":"43 1","pages":""},"PeriodicalIF":3.4,"publicationDate":"2024-04-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140560328","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-04-04DOI: 10.5194/egusphere-2024-855
Orie Cecil, Nicholas Cohn, Matthew Farthing, Sourav Dutta, Andrew Trautz
Abstract. Existing process-based models for simulating coastal foredune evolution largely use the same analytical approach for estimating wind induced surface shear stress distributions over spatially variable topography. Originally developed for smooth, low-sloping hills, these analytical models face significant limitations when the topography of interest exhibits large height-to-length ratios and/or steep, localized features. In this work, we utilize computational fluid dynamics (CFD) to examine the error trends of a commonly used analytical shear stress model for a series of idealized two-dimensional dune profiles. It is observed that the prediction error of the analytical model increases as compared to the CFD simulations for increasing height-to-length ratio and localized slope values. Furthermore, we explore two data-driven methodologies for generating alternative shear stress prediction models, namely, symbolic regression and linear, projection-based, non-intrusive reduced order modeling. These alternative modeling strategies demonstrate reduced overall error, but still suffer in their generalizability to broader sets of dune profiles outside of the training data. Finally, the impact of these improvements to aeolian sediment transport fluxes is examined to demonstrate that even modest improvements to the shear stress prediction can have significant impacts to dune evolution simulations over engineering-relevant timescales.
{"title":"Examination of Analytical Shear Stress Predictions for Coastal Dune Evolution","authors":"Orie Cecil, Nicholas Cohn, Matthew Farthing, Sourav Dutta, Andrew Trautz","doi":"10.5194/egusphere-2024-855","DOIUrl":"https://doi.org/10.5194/egusphere-2024-855","url":null,"abstract":"<strong>Abstract.</strong> Existing process-based models for simulating coastal foredune evolution largely use the same analytical approach for estimating wind induced surface shear stress distributions over spatially variable topography. Originally developed for smooth, low-sloping hills, these analytical models face significant limitations when the topography of interest exhibits large height-to-length ratios and/or steep, localized features. In this work, we utilize computational fluid dynamics (CFD) to examine the error trends of a commonly used analytical shear stress model for a series of idealized two-dimensional dune profiles. It is observed that the prediction error of the analytical model increases as compared to the CFD simulations for increasing height-to-length ratio and localized slope values. Furthermore, we explore two data-driven methodologies for generating alternative shear stress prediction models, namely, symbolic regression and linear, projection-based, non-intrusive reduced order modeling. These alternative modeling strategies demonstrate reduced overall error, but still suffer in their generalizability to broader sets of dune profiles outside of the training data. Finally, the impact of these improvements to aeolian sediment transport fluxes is examined to demonstrate that even modest improvements to the shear stress prediction can have significant impacts to dune evolution simulations over engineering-relevant timescales.","PeriodicalId":48749,"journal":{"name":"Earth Surface Dynamics","volume":"48 1","pages":""},"PeriodicalIF":3.4,"publicationDate":"2024-04-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140560415","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-04-04DOI: 10.5194/egusphere-2024-784
Elizabeth Orr, Taylor Schildgen, Stefanie Tofelde, Hella Wittmann, Ricardo Alonso
Abstract. Theory suggests that the response time of alluvial channel long-profiles to perturbations in climate is related to the magnitude of the forcing and the length of the system. Shorter systems may record a higher frequency of forcing compared to longer systems. Empirical field evidence that system length plays a role in the climate periodicity preserved within the sedimentary record is, however, sparse. The Toro Basin in the Eastern Cordillera of NW Argentina provides an opportunity to test these theoretical relationships as this single source-to-sink system contains a range of sediment deposits, located at varying distances from the source. A suite of eight alluvial fan deposits is preserved along the western flanks of the Sierra de Pascha. Farther downstream, a flight of cut-and-fill terraces have been linked to eccentricity-driven (100-kyr) climate cycles since ca. 500 ka. We applied cosmogenic radionuclide (10Be) exposure dating to the fan surfaces to explore (1) how channel responses to external perturbations may or may not propagate downstream, and (2) the differences in landscape response to forcing frequency as a function of channel length. We identified two generations of fan surfaces: the first (G1) records surface activity and abandonment between ca. 800 and 500 ka and the second (G2) within the last 100 kyr. G1 fans record a prolonged phase of net incision, which has been recognised throughout the Central Andes, and was likely triggered by enhanced 100-kyr global glacial cycles following the Mid-Pleistocene Transition. Relative fan surface stability followed, while 100-kyr cut-and-fill cycles occurred downstream, suggesting a disconnect in behaviour between the two channel reaches. G2 fans record higher frequency climate forcing, possibly the result of precessional forcing of climate (ca. 21/40-kyr timescales). The lack of a high-frequency signal farther downstream provides field support for theoretical predictions of a filtering of high-frequency climate forcing with increasing channel length. We show that multiple climate periodicities can be preserved within the sedimentary record of a single basin. Differences in the timing of alluvial fan and fluvial terrace development in the Toro Basin appears to be associated with how channel length affects fluvial response times to climate forcing as well as local controls on net incision, such as tectonic deformation.
摘要理论表明,冲积河道长剖面对气候扰动的响应时间与强迫的大小和系统的长度有关。与较长的系统相比,较短的系统可能会记录更高频率的强迫。然而,关于系统长度在沉积记录所保留的气候周期性中发挥作用的实地经验证据并不多。阿根廷西北部东科迪勒拉山系的托罗盆地为检验这些理论关系提供了一个机会,因为这个单一的源-汇系统包含一系列沉积物,与源头的距离各不相同。在帕夏山脉的西侧,保存着八组冲积扇沉积。在更远的下游,一连串的切填阶地与自约 500 ka 年以来的偏心驱动(100-kyr)气候周期有关。500ka。我们对扇面进行了宇宙放射性核素(10Be)暴露年代测定,以探索:(1)渠道对外部扰动的响应如何可能或不可能向下游传播;(2)景观对强迫频率的响应差异与渠道长度的函数关系。我们发现了两代扇面:第一代(G1)记录了大约 800 至 500 ka 期间扇面的活动和废弃,第二代(G2)记录了最近 100 kyr 期间扇面的活动和废弃。G1 扇面记录了整个中安第斯山脉长期的净侵蚀阶段,这很可能是中更新世过渡之后 100 千年全球冰川周期增强所引发的。随后,扇面表面相对稳定,而下游则出现了 100 千年的切填循环,这表明两个河道之间的行为脱节。G2 扇面记录了较高频率的气候作用,可能是气候前行作用的结果(时间尺度约为 21/40 千年)。下游缺乏高频信号为理论预测提供了实地支持,即随着河道长度的增加,高频气候作用会被过滤。我们的研究表明,在一个盆地的沉积记录中可以保留多种气候周期。托罗盆地冲积扇和河流阶地发育时间的差异似乎与河道长度如何影响河流对气候作用的响应时间以及当地对净侵蚀的控制(如构造变形)有关。
{"title":"Landscape response to tectonic deformation and cyclic climate change since ca. 800 ka in the southern Central Andes","authors":"Elizabeth Orr, Taylor Schildgen, Stefanie Tofelde, Hella Wittmann, Ricardo Alonso","doi":"10.5194/egusphere-2024-784","DOIUrl":"https://doi.org/10.5194/egusphere-2024-784","url":null,"abstract":"<strong>Abstract.</strong> Theory suggests that the response time of alluvial channel long-profiles to perturbations in climate is related to the magnitude of the forcing and the length of the system. Shorter systems may record a higher frequency of forcing compared to longer systems. Empirical field evidence that system length plays a role in the climate periodicity preserved within the sedimentary record is, however, sparse. The Toro Basin in the Eastern Cordillera of NW Argentina provides an opportunity to test these theoretical relationships as this single source-to-sink system contains a range of sediment deposits, located at varying distances from the source. A suite of eight alluvial fan deposits is preserved along the western flanks of the Sierra de Pascha. Farther downstream, a flight of cut-and-fill terraces have been linked to eccentricity-driven (100-kyr) climate cycles since ca. 500 ka. We applied cosmogenic radionuclide (<sup>10</sup>Be) exposure dating to the fan surfaces to explore (1) how channel responses to external perturbations may or may not propagate downstream, and (2) the differences in landscape response to forcing frequency as a function of channel length. We identified two generations of fan surfaces: the first (G1) records surface activity and abandonment between ca. 800 and 500 ka and the second (G2) within the last 100 kyr. G1 fans record a prolonged phase of net incision, which has been recognised throughout the Central Andes, and was likely triggered by enhanced 100-kyr global glacial cycles following the Mid-Pleistocene Transition. Relative fan surface stability followed, while 100-kyr cut-and-fill cycles occurred downstream, suggesting a disconnect in behaviour between the two channel reaches. G2 fans record higher frequency climate forcing, possibly the result of precessional forcing of climate (ca. 21/40-kyr timescales). The lack of a high-frequency signal farther downstream provides field support for theoretical predictions of a filtering of high-frequency climate forcing with increasing channel length. We show that multiple climate periodicities can be preserved within the sedimentary record of a single basin. Differences in the timing of alluvial fan and fluvial terrace development in the Toro Basin appears to be associated with how channel length affects fluvial response times to climate forcing as well as local controls on net incision, such as tectonic deformation.","PeriodicalId":48749,"journal":{"name":"Earth Surface Dynamics","volume":"38 1","pages":""},"PeriodicalIF":3.4,"publicationDate":"2024-04-04","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140560478","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-04-02DOI: 10.5194/egusphere-2024-808
Liran Goren, Eitan Shelef
Abstract. The plan-form geometry of branching drainage networks controls the topography of landscapes as well as their geomorphic, hydrologic, and ecologic functionality. The complexity of networks' geometry shows significant variability, from simple, straight channels that flow along the regional topographic gradient to intricate, tortuous flow patterns. This variability in complexity presents an enigma, as models show that it emerges independently of any heterogeneity in the environmental conditions. We propose to quantify networks' complexity based on the distribution of lengthwise asymmetry between paired flow pathways that diverge from a divide and rejoin at a junction. Using the lengthwise asymmetry definition, we show that the channel concavity index, describing downstream changes in channel slope, has a primary control on the plan-form complexity of natural drainage networks. An analytic model based on geomorphic scaling relations and optimal channel network simulations employing an energy minimization principle reveal that landscapes with low concavity channels attain stable plan-form configuration only through simple geometry. In contrast, landscapes with high-concavity channels achieve plan-form stability with various degrees of network complexity, including extremely complex geometries. Landscape evolution simulations demonstrate that the concavity index and its effect on the multiplicity of available geometries control the tendency of networks to preserve the legacy of former environmental conditions. Consistent with previous findings showing that channel concavity correlates with climate aridity, we find a significant empirical correlation between aridity and network complexity, suggesting a climatic signature embedded in the large-scale plan-form geometry of landscapes.
{"title":"Channel concavity controls plan-form complexity of branching drainage networks","authors":"Liran Goren, Eitan Shelef","doi":"10.5194/egusphere-2024-808","DOIUrl":"https://doi.org/10.5194/egusphere-2024-808","url":null,"abstract":"<strong>Abstract.</strong> The plan-form geometry of branching drainage networks controls the topography of landscapes as well as their geomorphic, hydrologic, and ecologic functionality. The complexity of networks' geometry shows significant variability, from simple, straight channels that flow along the regional topographic gradient to intricate, tortuous flow patterns. This variability in complexity presents an enigma, as models show that it emerges independently of any heterogeneity in the environmental conditions. We propose to quantify networks' complexity based on the distribution of lengthwise asymmetry between paired flow pathways that diverge from a divide and rejoin at a junction. Using the lengthwise asymmetry definition, we show that the channel concavity index, describing downstream changes in channel slope, has a primary control on the plan-form complexity of natural drainage networks. An analytic model based on geomorphic scaling relations and optimal channel network simulations employing an energy minimization principle reveal that landscapes with low concavity channels attain stable plan-form configuration only through simple geometry. In contrast, landscapes with high-concavity channels achieve plan-form stability with various degrees of network complexity, including extremely complex geometries. Landscape evolution simulations demonstrate that the concavity index and its effect on the multiplicity of available geometries control the tendency of networks to preserve the legacy of former environmental conditions. Consistent with previous findings showing that channel concavity correlates with climate aridity, we find a significant empirical correlation between aridity and network complexity, suggesting a climatic signature embedded in the large-scale plan-form geometry of landscapes.","PeriodicalId":48749,"journal":{"name":"Earth Surface Dynamics","volume":"31 1","pages":""},"PeriodicalIF":3.4,"publicationDate":"2024-04-02","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140560473","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}