Understanding the rates and mechanisms of erosion by subglacial quarrying is a major unsolved problem in geomorphology. Stress enhancement due to load concentration on bedrock ledges between cavities is hypothesized to drive the growth of fractures. Prior work assumed the formation of vertically oriented tensile fractures at the downstream margins of cavities as the controlling process, but did not account for the evolution of the stress field as fractures lengthen, and in particular the dominance of the shearing mode at fracture tips. We used 2D finite element analysis and J-integral methods to analyze stress intensity factors and fracture growth potentials at the tips of preexisting fractures in loaded bedrock steps, taking into account normal and shear components and measured rock strengths. By examining different step heights, step riser angles, rock types, prior fracture locations and orientations, and extents of ice-rock contact zones, we identified some situations favorable for fracture growth, especially in brittle rock types. Typically, however, the growth direction will not be vertically downward but angled up-glacier away from the step riser, a situation unfavorable for quarrying. Moreover, in many situations, the normal stress across fracture planes will be compressive. Non-vertical step risers buttress the bedrock and also suppress fracture growth. In contrast, reducing the sizes of ice-rock contact zones not only increases the loading magnitude, as previously recognized, but also increases intensification of tensile stress at the tips of fractures located just up-glacier. Thus, larger cavities, and hence, fast sliding and low effective pressures, favor quarrying more strongly than previously recognized.
Erosional perturbations from changes in climate or tectonics are recorded in the profiles of bedrock rivers, but these signals can be challenging to unravel in settings with non-uniform lithology. In layered rocks, the surface lithology at a given location varies through time as erosion exposes different layers of rock. Recent modeling studies have used the Stream Power Model (SPM) to highlight complex variations in erosion rates that arise in bedrock rivers incising through layered rocks. However, these studies do not capture the effects of coarse sediment cover on channel evolution. We use the “Stream Power with Alluvium Conservation and Entrainment” (SPACE) model to explore how sediment cover influences landscape evolution and modulates the topographic expression of erodibility contrasts in horizontally layered rocks. We simulate river evolution through alternating layers of hard and soft rock over million-year timescales with a constant and uniform uplift rate. Compared to the SPM, model runs with sediment cover have systematically higher channel steepness values in soft rock layers and lower channel steepness values in hard rock layers. As more sediment accumulates, the contrast in steepness between the two rock types decreases. Effective bedrock erodibilities back-calculated assuming the SPM are strongly influenced by sediment cover. We also find that sediment cover can significantly increase total relief and timescales of adjustment toward landscape-averaged steady-state topography and erosion rates.
Bedload transport can fluctuate considerably over relatively short periods of time and for a given quasi-constant flow rate. What are the implications of replacing the fluctuating signal with a smoothed signal when calculating bedload transport using averaged values, as is common practice? This question was investigated with the BedloadR code, which allows 1D bedload calculation as well as Monte Carlo simulations using a new data set collected in the Severaisse River (French Ecrins massif). Four bedload equations (Camenen & Larson, 2005, https://doi.org/10.1016/j.ecss.2004.10.019; Meyer-Peter & Mueller, 1948; Parker, 1990, https://doi.org/10.1080/00221689009499058; Recking, 2013a, https://doi.org/10.1061/(asce)hy.1943-7900.0000653) were selected for their performance relative to the measured bedload (except for and Meyer-Peter and Mueller) and because each equation has a different mathematical form and degree of nonlinearity. They were used in a Monte Carlo approach, with input probability distributions fitted to the measured river width, slope, bed grain-size distribution, and to the associated (computed) Shields stress. The results show that accounting for natural variability in the calculation reproduces bedload fluctuations well. But overall, when calculating the bedload volume transported by a flow event, accounting for variability systematically leads to higher estimated volumes (of the order of 20%) than those obtained with a deterministic approach using average input parameters. This is a direct consequence of the nonlinearity of the equations.
Permafrost warming has been observed all around the Arctic, however, variations in temperature trends and their drivers remain poorly understood. We present a comprehensive analysis of climatic changes spanning 25 years (1998–2023) at Bayelva (78.92094°N, 11.83333°E) on Spitzbergen, Svalbard. The quality controlled hourly data set includes air temperature, radiation fluxes, snow depth, rainfall, active layer temperature and moisture, and, since 2009, permafrost temperature. Our Bayesian trend analysis reveals an annual air temperature increase of 0.9 ± 0.5°C/decade and strongest warming in September and October. We observed a significant shortening of the snow cover by −14 ± 8 days/decade, coupled with reduced winter snow depth. The active layer simultaneously warmed by 0.6 ± 0.7°C/decade at the top and 0.8 ± 0.5°C/decade at the bottom. While the soil surface got drier, in particular during summer, soil moisture below increased in accordance with the longer unfrozen period and higher winter temperatures. The thawed period prolonged by 10–15 days/decade at different depths. In contrast to earlier top-soil warming, we observed stable temperatures since 2010 and only little permafrost warming (0.14 ± 0.13°C/decade). This is likely due to recently stable winter air temperature and continuously decreasing winter snow depth. This recent development highlights a complex interplay among climate and soil variables. Our distinctive long-term data set underscores (a) the changes in seasonal warming patterns, (b) the influential role of snow cover decline, and (c) that air temperature alone is not a sufficient indicator of change in permafrost environments, thereby highlighting the importance of investigating a wider range of parameters, such as soil moisture and snow characteristics.
Wave-driven erosion of marsh boundaries is a major cause of marsh loss, but little research has captured the effect of seasonal differences on marsh-edge retreat rates to illuminate temporal patterns of when the majority of this erosion is occurring. Using five surface models captured over a study year of a marsh with a steep escarped boundary in South San Francisco Bay, we find a pronounced seasonal signal, where rapid marsh retreat in the spring and summer is driven by a strong sea breeze but little change is found in the marsh-edge position in the fall and winter. We found accretion in the mudflat transition region close to the marsh boundary in the calmer seasons however, suggesting intertwined morphodynamics of mudflats and the eroding marsh-scarp. We observed large spatial heterogeneity in retreat rates within seasons, but less on longer (annual and decadal) timescales. The relationship between marsh-edge retreat rates and properties of the wave field nearby is explored and contextualized against extant relationships, but our results speak to the difficulty in addressing spatial erosion/accretion variability on short (seasonal) timescales with simple models.
Reconstructing historical climate change from deep ground temperature measurements in cold regions is often complicated by the presence of permafrost. Existing methods are typically unable to account for latent heat effects due to the freezing and thawing of the active layer. In this work, we propose a novel method for reconstructing historical ground surface temperature (GST) from borehole temperature measurements that accounts for seasonal thawing and refreezing of the active layer. Our method couples a recently developed fast numerical modeling scheme for two-phase heat transport in permafrost soils with an ensemble-based method for approximate Bayesian inference. We evaluate our method on two synthetic test cases covering both cold and warm permafrost conditions as well as using real data from a 100 m deep borehole on Sardakh Island in northeastern Siberia. Our analysis of the Sardakh Island borehole data confirms previous findings that GST in the region have likely risen by 5–9°C between the pre-industrial period of 1750–1855 and 2012. We also show that latent heat effects due to seasonal freeze-thaw have a substantial impact on the resulting reconstructed surface temperatures. We find that neglecting the thermal dynamics of the active layer can result in biases of roughly −1°C in cold conditions (i.e., mean annual ground temperature below −5°C) and as much as −2.6°C in warmer conditions where substantial active layer thickening (>200 cm) has occurred. Our results highlight the importance of considering seasonal freeze-thaw in GST reconstructions from permafrost boreholes.
Granular flows are ubiquitous in nature with single flows traversing a wide range of dynamic conditions from initiation to deposition. Many of these flows are responsible for significant hazards and can generate remotely detectable seismic signals. These signals provide a potential for real-time flow measurements from a safe distance. To fully realize the benefit of seismic measurements, basal-granular forces must be linked to macroscopic internal flow dynamics across a wide range of flow conditions. We utilize discrete element simulations to observe dry and submerged granular flows under plane-shear and inclined-flow configurations, relating bulk kinematics to basal-force distributions. We find that the power and frequency of force fluctuations scale with non-dimensional shear rate (I). This scaling tracks three pre-established regimes that are described by μ(I) rheology: (a) an intermittent particle rearrangement regime, where basal forces are dominated by low frequencies; (b) an intermediate regime where basal forces start to increase in frequency while showing correlations in space and (c) a fully collisional regime where the signal is nearly flat up to a cutoff frequency. We further identify a newly defined fourth regime that marks a “phase change” from the intermediate to collisional regime where increases in basal force fluctuations with increasing shear rates stalls as the granular bed dilates, partially destroying the contact network. This effort suggests that basal forces can be used to interpret complex granular processes in geophysical flows.
Estuaries worldwide are susceptible and adapting to climate change (CC) impacts from both the river and coastal boundaries. Furthermore, engineering efforts are undertaken to improve flood safety, to claim land for human use or for port operations, which change estuary morphology. This paper aims to gain an understanding of the combined effects of CC and human interventions on the estuarine-wide morphological response by analyzing the sediment infilling of highly engineered estuaries. A schematized process-based morphodynamic model is used (Delft3D-FM, in 2DH), resembling a highly engineered estuary in the Rhine-Meuse Delta, The Netherlands. Three types of changes were implemented, both in isolation and in combination: (a) local interventions (changing channel depth or wetland area), (b) upstream human interventions (changing fluvial sediment supply) and (c) extreme CC scenarios (with projections for the future forcings and bathymetry). Results show that a CC scenario can elicit both positive and negative changes in the estuary's sediment budget. The direction and magnitude of the change depend on the local intervention and can align with the effect of the local intervention, intensifying its impact. The combined effects can even reverse the sign of the sediment budget. This stresses the need of analyzing CC impacts in combination with human interventions. Additionally, a relationship was identified which quantifies how a change in peak flow velocity due to both local interventions and sea-level rise affects the annual sediment budget. These findings can help determine how local interventions affect morphodynamics of engineered estuaries in present and future climates.