Landscape Evolution Model Research Articles

Designing, building and assessing a geomorphically reconstructed postmining landscape: A case study of the Santa Engracia mine, Spain

AbstractNew technology allows the reconstruction of postmining landforms using geomorphic design principles. It is important that such designs be evaluated and if needed, redesigned or reshaped so that soil loss is minimised and to ensure the landscape is geomorphically and ecologically integrated with the surrounding landscape. One tool to assess geomorphic landforms is to use a computer‐based landscape evolution model (LEM). LEMs allow different designs to be input and will highlight where erosion will occur and type of erosion (i.e. sheetwash, riling, gullying) as well as erosion rate. At the Santa Engracia abandoned mine (East‐Central Spain), postmining landscapes were designed using geomorphic principles (GeoFluv method and Natural Regrade software) and later constructed. The SIBERIA LEM was used to assess the erosional behaviour of these landscapes. Using suitable topsoil, vegetation and an organic blanket reduces erosion, and if vegetation can be established, the modelling demonstrates minimal gully erosion. The erosion forecast (5.3 to 6.3 t ha−1 year−1) is significantly lower than the initial surface (~350 t ha−1 year−1) using conventional (terraced) mine restoration. The predicted erosion rates and gullying are less than for the unmined (natural) Alto Tajo environment. Importantly, with the ability to spatially forecast gully location, erosion reduction measures can be undertaken. The method described here provides a robust assessment procedure and highlights the potential strengths and weakness of a design therefore supporting lower cost construction and repair with a higher chance of restoration success. The combination of geomorphic landform design and assessment using a LEM for this project (LIFE RIBERMINE) presents a new standard for mine rehabilitation in Europe.

Short Communication: Numerically simulated time to steady state is not a reliable measure of landscape response time

Abstract. Quantifying the timescales over which landscapes evolve is critical for understanding past and future environmental change. Computational landscape evolution models are one tool among many that have been used in this pursuit. We compare numerically modeled times to reach steady state for a landscape adjusting to an increase in rock uplift rate. We use three different numerical modeling libraries and explore the impact of time step, grid type, numerical method for solving the erosion equation, and metric for quantifying the time to steady state. We find that modeled time to steady state is impacted by all of these variables. Time to steady state varies inconsistently with time step length, both within a single model and among different models. In some cases, drainage rearrangement extends the time to reach steady state, but this is not consistent in all models or grid types. The two sets of experiments operating on Voronoi grids have the most consistent times to steady state when comparing across time step and metrics. On a raster grid, if we force the drainage network to remain stable, time to steady state varies much less with computational time step. In all cases we find that many measures of modeled time to steady state are longer than that predicted by an analytical equation for bedrock river response time. Our results show that the predicted time to steady state from a numerical model is, in many cases, more reflective of drainage rearrangement and numerical artifacts than the time for an uplift wave to propagate through a fixed drainage network.

A landslide runout model for sediment transport, landscape evolution, and hazard assessment applications

Abstract. We developed a new rule-based, cellular-automaton algorithm for predicting the hazard extent, sediment transport, and topographic change associated with the runout of a landslide. This algorithm, which we call MassWastingRunout (MWR), is coded in Python and implemented as a component for the package Landlab. MWR combines the functionality of simple runout algorithms used in landscape evolution and watershed sediment yield models with the predictive detail typical of runout models used for landslide inundation hazard mapping. An initial digital elevation model (DEM), a regolith depth map, and the location polygon of the landslide source area are the only inputs required to run MWR to model the entire runout process. Runout relies on the principle of mass conservation and a set of topographic rules and empirical formulas that govern erosion and deposition. For the purpose of facilitating rapid calibration to a site, MWR includes a calibration utility that uses an adaptive Bayesian Markov chain Monte Carlo algorithm to automatically calibrate the model to match observed runout extent, deposition, and erosion. Additionally, the calibration utility produces empirical probability density functions of each calibration parameter that can be used to inform probabilistic implementation of MWR. Here we use a series of synthetic terrains to demonstrate basic model response to topographic convergence and slope, test calibrated model performance relative to several observed landslides, and briefly demonstrate how MWR can be used to develop a probabilistic runout hazard map. A calibrated runout model may allow for region-specific and more insightful predictions of landslide impact on landscape morphology and watershed-scale sediment dynamics and should be further investigated in future modeling studies.

An advection–diffusion equation with a generalized advection term: Well-posedness analysis and examples

We consider a Cauchy–Dirichlet problem for a semilinear advection–diffusion equation with a generalized advection term. Specific examples include an incision–diffusion landscape evolution model and a viscous Hamilton–Jacobi equation with an absorbing gradient term. We establish existence and uniqueness of a weak solution and its continuous dependence on initial data and a source term.

Exploiting the land use to predict shallow landslide susceptibility: A probabilistic implementation of LAPSUS-LS

Due to the significant role of land use on the occurrence of rainfall-induced shallow landslides, this factor is commonly employed as a landslide susceptibility predictor. However, the land use classification is oftentimes very broad, neglecting the proven mechanical and hydrogeological role of the land management on slope stability. Given the necessity of including spatially distributed and management-specific inputs, the process-based landscape evolution model LAPSUS-LS was chosen and adapted to achieve a probabilistic approach which takes into account land management as an input by adopting management-specific values of root cohesion. The model was applied to four test sites in the Oltrepò Pavese (Italy), where different vineyard management techniques play a significant role in triggering landslides. The results for the four test areas had, cumulatively, an Area Under the Roc curve greater than 0.73, with false negative cells being < 1 % of the total for all simulations. In the model’s application, land use practices characterised by higher root cohesion proved to benefit slope stability, whereas tilled vineyards, shrublands and abandoned vineyards were more prone to the formation of shallow landslides. In addition, we found that the inclusion of management-specific input parameters produced more accurate outputs and that in catchments characterised by average slope angles lower than 15°, varying the vineyard management, did not appear to affect the landslide susceptibility. Due to the model’s high dependency on the land use and its ability to include land management, it can take into account the spatial variability of input values such as the root cohesion. Additionally, it can be applied i) to manage current conditions, ii) to explore future land use change, iii) to study less invasive yet beneficial land use management change scenarios and iv) provide farmers of at-risk areas insight on how to improve slope stability.

Modeling Sediment Fluxes From Debris‐Rich Basal Ice Layers

AbstractSediment erosion, transport, and deposition by glaciers and ice sheets play crucial roles in shaping landscapes, provide important nutrients to downstream ecosystems, and preserve key indicators of past climate conditions in the geologic record. While previous work has quantified sediment fluxes from subglacial meltwater, we also observe sediment entrained within basal ice, transported by the flow of the glacier itself. However, the formation and evolution of these debris‐rich ice layers remains poorly understood and rarely represented in landscape evolution models. Here, we identify a characteristic sequence of basal ice layers at Mendenhall Glacier, Alaska. We develop a numerical model of frozen fringe and regelation processes that describes the co‐evolution of this sequence and explore the sensitivity of the model to key properties of the subglacial sedimentary system, using the Instructed Glacier Model to parameterize ice dynamics. Then, we run numerical simulations over the spatial extent of Mendenhall Glacier, showing that the sediment transport model can predict the observed basal ice stratigraphy at the glacier's terminus. From the model results, we estimate basal ice layers transport between 23,300 and 39,800 of sediment, mostly entrained in the lowermost ice layers nearest to the bed, maximized by high effective pressures and slow, convergent flow fields. Overall, our results highlight the role of basal sediment entrainment in delivering eroded material to the glacier terminus and indicate that this process should not be ignored in broader models of landscape evolution.

Sediment Transport and Flood Risk: Impact of Newly Constructed Embankments on River Morphology and Flood Dynamics in Kathmandu, Nepal

AbstractFloodplain encroachment by embankments heightens flood risk. This is exacerbated by climate change and land‐use modifications. This paper assesses the impact of embankments on sediment transport, channel geometry, conveyance capacity, and flood inundation of a reach of the Nakkhu River, Nepal. Using the CAESAR‐Lisflood landscape evolution model based on a 2‐m digital elevation model, we simulate four flood scenarios with and without embankments and sediment transport: a historical 25‐year return period flood event used to design the embankments, 50‐year, 100‐year, and 1000‐year return period flood events forecast using the Generalized Logistic Model (using data from 1992 to 2017). Our results indicate that flow confinement by embankments reduces inundation by 99% (from 22.5 to 0.3 ha) for the historical 25‐year flood discharge of 42.23 and by 15% (from 28.8 to 24.4 ha) for the 1000‐year return period flood discharge of 95 (similar to a 25‐year maximum mid‐future). The presence of embankments increases downstream sediment transport by more than 32% for all flood scenarios considered. Inclusion of sediment transport leads to a fivefold increase in predicted inundation area for a 25‐year maximum mid‐future flood compared to the no‐sediment case in the embanked channel. Changes in channel geometry due to sedimentation significantly reduce conveyance capacity increasing overtopping flood risk, particularly where the channel is sinuous or located on flat terrain. Our results indicate that sediment erosion in outer meanders may threaten embankment stability by promoting undercuts. It is recommended that sediment transport effects be factored into embankment design and floodplain planning.

Estimation of Denudation Parameters and River Capture Events From Neural Network Inverse Modeling of River Profiles and Thermo‐ and Geochronology Data

AbstractEarth's topography represents the cumulative effects of tectonics and surface processes modulated by climate and lithology. These factors shape landscapes through time. River profiles can be inverted to estimate the rock uplift histories or lithology‐specific erodibilities. However, river systems are dynamic and evolve in response to spatial and temporal internal dynamics, such as river capture events. Here, we present a modeling framework to infer denudation rates from the inversion of river profiles and thermo‐ and geochronology data. We achieve this by coupling a landscape evolution model and an efficient inverse modeling scheme to infer poorly resolved erosional and tectonic parameters. An application of the approach is presented for the Neckar catchment, southwest Germany, characterized by stark lateral variation in bedrock erodibility and rock uplift, and that have demonstrably undergone multiple river capture events. Different end‐member scenarios are explored in the simulations. First, we test uniform and spatial variability in rock uplift rate and bedrock erodibility, and second, temporal variations in rock uplift rate and base level. Finally, we simulate river capture events by adding upstream sections (drainage area) at specific times and locations within the fluvial network. We find that spatial variation in rock uplift rate is necessary to reproduce the Neckar's river profile while honoring analytical observations. Simulations integrating river captures allow improved river profile predictions of specific tributaries of the Neckar catchment, leading to potentially more realistic erodibility and rock uplift history estimates. The time and location of the capture events determined from the modeling agree with previous estimations from geological evidence.

Topographic metrics for unveiling fault segmentation and tectono-geomorphic evolution with insights into the impact of inherited topography, Ulsan Fault Zone, South Korea

Abstract. Quantifying today's topography can provide insights into landscape evolution and its controls, since present topography represents a cumulative expression of past and present surface processes. The Ulsan Fault Zone (UFZ) is an active fault zone on the southeastern Korean Peninsula that was reactivated as a reverse fault around 5 Ma. The UFZ strikes NNW–SSE and dips eastward. This study investigates the relative tectonic activity along the UFZ and the landscape evolution of the hanging-wall side of the UFZ, focusing on neotectonic perturbations using 10Be-derived catchment-averaged denudation rates and bedrock incision rates, topographic metrics, and a landscape evolution model. Five geological segments were identified along the fault, based on their relative tectonic activity and fault geometry. We simulated four cases of landscape evolution to investigate the geomorphic processes and accompanying topographic changes in the study area in response to fault movement. Model results reveal that the geomorphic processes and the patterns of topographic metrics (e.g., χ anomalies) depend on inherited topography (i.e., the topography that existed prior to reverse fault reactivation of the UFZ). On the basis of this important model finding and additional topographic metrics, we interpret the tectono-geomorphic history of the study area as follows: (1) the northern part of the UFZ has been in a transient state and is in topographic and geometric disequilibrium, so this segment underwent asymmetric uplift (westward tilting) prior to reverse faulting on the UFZ around 5 Ma, and (2) its southern part was negligibly influenced by the asymmetric uplift before reverse faulting. Our study demonstrates the utility of topographic metrics as reliable criteria for resolving fault segments. Together with landscape evolution modeling, topographic metrics provide powerful tools for examining the influence of inherited topography on present topography and for the elucidation of tectono-geomorphic histories.

The impacts of climate change, early agriculture and internal fluvial dynamics on paleo-flooding episodes in Central China

Floods clustered in episodes are the most prevalent natural disaster worldwide, causing substantial economic and human losses. Although these events are often linked to time-periods of extreme rainstorms and unique atmospheric circulation patterns, the river basin characteristics affected by anthropogenic land use changes could exert a strong influence. However, the way and extent of how land use changes across different time scales affect flooding periods are still unclear, especially considering the historical land use changes. This study uses the Landlab landscape evolution model, coupled with an evapotranspiration model, to investigate the forcing factors for the paleo-flooding trends in the Wei River catchment over the last 5000 years. The results indicate that the flooding period from 4400 to 4000 BP was caused by an increase of 28 % in antecedent moisture content as well as a decrease of 28 % in its spatial variability, which are primarily due to climate change, and that the contribution of land-use change is less than 5 %. The increases of about 14 % and 8 % in main channel sedimentation rate play a leading role in flood generation during the time periods from 3400 to 2800 BP and 2000–1400 BP, respectively. These two periods of increased flooding are primarily caused by the erosional effects of increasing anthropogenic land use, whose contributions range from 33 % to 64 %. Furthermore, based on our modelling results, we suggest that the downstream propagation of the main flooding locations, from the Wei River to the lower reaches of the Yellow River, can be explained by the downstream migrating sediment wave. In conclusion, our simulation results give new insights into the causes of Holocene flooding periods in the middle Yellow River from the perspective of dynamic changes in catchment characteristics, which is helpful to improve regional flood risk management under future climate change and anthropogenic activities.

From flood to turbidity current: combined models to simulate continent to ocean sediment transport in the Var system, France

Turbidity currents are a form of gravity-driven flow that occurs in subaqueous environments. These currents contain a low volumetric percentage of particles but are very important in transporting them from coastal to deep-marine environments. In the Var sediment routing system, southeastern France, there is a sedimentary record of turbidity currents formed by submarine landslides and flooding of the Var River. In 2006, there was continuous monitoring of the marine section that recorded four turbidity currents. To explore the connectivity between the terrestrial source and marine basin we linked two reduced complexity models of sediment transport and landscape change. The aim of the paper is to combine two models, one dedicated to the prediction of river water and sediment run-off, and the second dedicated to the transport, erosion, and sedimentation by a turbidity current. The landscape evolution model CAESAR-Lisflood was used to model discharge and sediment yield. The behaviour of turbidity currents was modelled with the cellular automata model CATS. From calibration of the two models against observations of discharge and suspended sediment in the Var River system, we find that in 2006, two rainfall events would have led to hyperpycnal turbidity currents. These two events match well with two of the four-recorded events. Focusing on the largest event we find that the source pulse of sediment from the terrestrial environment is short-lived, less than half a day, but contains a significant quantity of fine particles. The event transferred on the order of 105 m3 of suspended sediment from source to sink. This study demonstrates that hyperpycnal turbidity currents can be generated with concentrations as low as 2-6 kg/m3 at the Var River mouth far below the theoretical 40 kg/m3 threshold, suggesting active convective sedimentation in the surface plume and sea-water dilution at the flooding river mouth. A substantial amount of sediment (35% of the input volume) is directly transferred towards the deepest part of the system during a short-living hyperpycnal turbidity current. However, a considerable (65%) part remains in the surface as a hypopycnal plume feeding the hemipelagic sedimentation.

Coevolution of craton margins and interiors during continental break-up

Many cratonic continental fragments dispersed during the rifting and break-up of Gondwana are bound by steep topographic landforms known as ‘great escarpments’1–4, which rim elevated plateaus in the craton interior5,6. In terms of formation, escarpments and plateaus are traditionally considered distinct owing to their spatial separation, occasionally spanning more than a thousand kilometres. Here we integrate geological observations, statistical analysis, geodynamic simulations and landscape-evolution models to develop a physical model that mechanistically links both phenomena to continental rifting. Escarpments primarily initiate at rift-border faults and slowly retreat at about 1 km Myr−1 through headward erosion. Simultaneously, rifting generates convective instabilities in the mantle7–10 that migrate cratonward at a faster rate of about 15–20 km Myr−1 along the lithospheric root, progressively removing cratonic keels11, driving isostatic uplift of craton interiors and forming a stable, elevated plateau. This process forces a synchronized wave of denudation, documented in thermochronology studies, which persists for tens of millions of years and migrates across the craton at a comparable or slower pace. We interpret the observed sequence of rifting, escarpment formation and exhumation of craton interiors as an evolving record of geodynamic mantle processes tied to continental break-up, upending the prevailing notion of cratons as geologically stable terrains.

Signatures of wave erosion in Titan's coasts.

The shorelines of Titan's hydrocarbon seas trace flooded erosional landforms such as river valleys; however, it is unclear whether coastal erosion has subsequently altered these shorelines. Spacecraft observations and theoretical models suggest that wind may cause waves to form on Titan's seas, potentially driving coastal erosion, but the observational evidence of waves is indirect, and the processes affecting shoreline evolution on Titan remain unknown. No widely accepted framework exists for using shoreline morphology to quantitatively discern coastal erosion mechanisms, even on Earth, where the dominant mechanisms are known. We combine landscape evolution models with measurements of shoreline shape on Earth to characterize how different coastal erosion mechanisms affect shoreline morphology. Applying this framework to Titan, we find that the shorelines of Titan's seas are most consistent with flooded landscapes that subsequently have been eroded by waves, rather than a uniform erosional process or no coastal erosion, particularly if wave growth saturates at fetch lengths of tens of kilometers.

Palaeolake Jeinimeni and ice-marginal geomorphology of the central Patagonian Lago General Carrera/Buenos Aires ice-lobe (46.5°S)

Robust geomorphological interpretations of ice-marginal glacial, meltwater and palaeolake landforms are needed to constrain models of landform evolution that can be used to underpin high-resolution geochronologies of deglaciation. The Lago General Carrera/Bueno Aires (LGCBA) ice lobe of the former Patagonian Ice Sheet is one of the most researched in Patagonia, yet there remain contrasting landform evolution models partly because the spatiotemporal relationships between landforms in the LGCBA basin and tributary valleys have not been fully investigated. To this end we present a detailed, high-resolution glacial geomorphological map of the Los Antiguos sector of the LGCBA basin, and the tributary Río Jeinimeni valley, alongside sediment-landform assemblages from a broad suite of ice-marginal, subglacial, glaciolacustrine, glaciofluvial and paraglacial landforms. We introduce evidence for a previously unmapped ice-dammed lake within the Jeinimeni valley, which may have existed during the Last Glacial Maximum. Multiple isolated, ice-marginal proglacial lakes, which formed during the recession of the LGCBA ice lobe, were also mapped. The landform evidence is synthesised into a refined landscape evolution model for the LGCBA basin and Jeinimeni valley during the Last Glacial Interglacial Transition. The data presented demonstrate that previously reconstructed high LGCBA lake-level stands above 500 m asl were likely glaciolacustrine landforms created by ice-marginal lakes dammed by the LGCBA ice lobe. The mapping in the Jeinimeni valley highlights the potential for renewed dating investigations of landforms in tributary valleys to provide maximum limiting ages for the retreat of their associated glaciers. Furthermore, the identification of glaciolacustrine varved sediments associated with palaeolake Jeinimeni and the Los Antiguos sector of the LGCBA palaeolake, highlights the potential of these annually-resolved archives to infer spatiotemporal dynamics of main ice-lobe/tributary glaciers to improve palaeoenvironmental and palaeoclimatic reconstructions in the region.

Testing Mantle Convection Simulations With Paleobiology and Other Stratigraphic Observations: Examples From Western North America

AbstractMantle convection plays a fundamental role in driving evolution of oceanic and continental lithosphere. In turn it impacts a broad suite of processes operating at or close to Earth's surface including landscape evolution, glacio‐eustasy, magmatism, and climate. A variety of theoretical approaches now exist to simulate mantle convection. Outputs from such simulations are being used to parameterize models of landscape evolution and basin formation. However, the substantial body of existing simulations has generated a variety of conflicting views on the history of dynamic topography, its evolution and key parameters for modeling mantle flow. The focus of this study is on developing strategies to use large‐scale quantitative stratigraphic observations to assess model predictions and identify simulation parameters that generate realistic predictions of Earth surface evolution. Spot measurements of uplift or subsidence provide useful target observations for models of dynamic topography, but finding areas where tectonics have not also influenced vertical motions is challenging. To address this issue, we use large inventories of stratigraphic data from across North America with contextual geophysical and geodetic data to constrain the regional uplift and subsidence history. We demonstrate that a suite of typical geodynamic simulations struggle to match the amplitude, polarity and timing of observed vertical motions. Building on recent seismological advances, we then explore strategies for understanding patterns of continental uplift and subsidence that incorporate (and test) predicted evolution of the lithosphere, asthenosphere and deep mantle. Our results demonstrate the importance of contributions from the uppermost mantle in driving vertical motions of continental interiors.

Catchment Coevolution and the Geomorphic Origins of Variable Source Area Hydrology

AbstractFeatures of landscape morphology—including slope, curvature, and drainage dissection—are important controls on runoff generation in upland landscapes. Over long timescales, runoff plays an essential role in shaping these same features through surface erosion. This feedback between erosion and runoff generation suggests that modeling long‐term landscape evolution together with dynamic runoff generation could provide insight into hydrological function. Here we examine the emergence of variable source area runoff generation in a new coupled hydro‐geomorphic model that accounts for water balance partitioning between surface flow, subsurface flow, and evapotranspiration as landscapes evolve over millions of years. We derive a minimal set of dimensionless numbers that provide insight into how hydrologic and geomorphic parameters together affect landscapes. Across the parameter space we investigated, model results collapsed to a single inverse relationship between the dimensionless relief and the ratio of catchment quickflow to discharge. Furthermore, we found an inverse relationship between the Hillslope number, which describes topographic relief relative to aquifer thickness, and the proportion of the landscape that was variably saturated. While the model generally produces fluvial topography visually similar to simpler landscape evolution models, certain parameter combinations produce wide valley bottom wetlands and non‐dendritic, trellis‐like drainage networks, which may reflect real conditions in some landscapes where aquifer gradients become decoupled from topography. With these results, we demonstrate the power of hydro‐geomorphic models for generating new insights into hydrological processes, and also suggest that subsurface hydrology may be integral for modeling aspects of long‐term landscape evolution.

Stream hydrology controls on ice cliff evolution and survival on debris-covered glaciers

Abstract. Ice cliffs are melt hot spots that contribute disproportionately to melt on debris-covered glaciers. In this study, we investigate the impact of supraglacial stream hydrology on ice cliffs using in situ and remote sensing observations, streamflow measurements, and a conceptual geomorphic model of ice cliff backwasting applied to ice cliffs on Kennicott Glacier, Alaska. We found that 33 % of ice cliffs (accounting for 69 % of the ice cliff area) are actively influenced by streams, while half are nearer than 10 m from the nearest stream. Supraglacial streams contribute to ice cliff formation and maintenance by horizontal meandering, vertical incision, and debris transport. These processes produce an undercut lip at the ice cliff base and transport clasts up to tens of centimeters in diameter, preventing reburial of ice cliffs by debris. Stream meander morphology reminiscent of sedimentary river channel meanders and oxbow lakes produces sinuous and crescent ice cliff shapes. Stream avulsions result in rapid ice cliff collapse and local channel abandonment. Ice cliffs abandoned by streams are observed to be reburied by supraglacial debris, indicating a strong role played by streams in ice cliff persistence. We also report on a localized surge-like event at the glacier's western margin which drove the formation of ice cliffs from crevassing; these cliffs occur in sets with parallel linear morphologies contrasting with the crescent planform shape of stream-driven cliffs. The development of landscape evolution models may assist in quantifying the total net effect of these processes on steady-state ice cliff coverage and mass balance, contextualizing them with other drivers including supraglacial ponds, differential melt, ice dynamics, and collapse of englacial voids.

Inclusion of bedrock vadose zone in dynamic global vegetation models is key for simulating vegetation structure and function

Abstract. Across many upland environments, soils are thin and plant roots extend into fractured and weathered bedrock where moisture and nutrients can be obtained. Root water extraction from unsaturated weathered bedrock is widespread and, in many environments, can explain gradients in vegetation community composition, transpiration, and plant sensitivity to climate. Despite increasing recognition of its importance, the “rock moisture” reservoir is rarely incorporated into vegetation and Earth system models. Here, we address this weakness in a widely used dynamic global vegetation model (DGVM; LPJ-GUESS). First, we use a water flux-tracking deficit approach to more accurately parameterize plant-accessible water storage capacity across the contiguous United States, which critically includes the water in bedrock below depths typically prescribed by soil databases. Secondly, we exploit field-based knowledge of contrasting plant-available water storage capacity in weathered bedrock across two bedrock types in the Northern California Coast Ranges as a detailed case study. For the case study in Northern California, climate and soil water storage capacity are similar at the two study areas, but the site with thick weathered bedrock and ample rock moisture supports a temperate mixed broadleaf–needleleaf evergreen forest, whereas the site with thin weathered bedrock and limited rock moisture supports an oak savanna. The distinct biomes, seasonality and magnitude of transpiration and primary productivity, and baseflow magnitudes only emerge from the DGVM when a new and simple subsurface storage structure and hydrology scheme is parameterized with storage capacities extending beyond the soil into the bedrock. Across the contiguous United States, the updated hydrology and subsurface storage improve annual evapotranspiration estimates as compared to satellite-derived products, particularly in seasonally dry regions. Specifically, the updated hydrology and subsurface storage allow for enhanced evapotranspiration through the dry season that better matches actual evapotranspiration patterns. While we made changes to both the subsurface water storage capacity and the hydrology, the most important impacts on model performance derive from changes to the subsurface water storage capacity. Our findings highlight the importance of rock moisture in explaining and predicting vegetation structure and function, particularly in seasonally dry climates. These findings motivate efforts to better incorporate the rock moisture reservoir into vegetation, climate, and landscape evolution models.

Readjustment of fluvial systems to Quaternary tectonic disruption in the Central Afar (Ethiopia)

The Afar region constitutes the diffuse triple junction adjoining the Red Sea, Gulf of Aden, and Main Ethiopian Rifts. Extension in the Afar is distributed between on-axis, magmatic crustal spreading, and off-axis, amagmatic extension, such as in Dobe graben. The timing of initiation and development of these amagmatic grabens is key to understanding the final stages of continental extension/breakup in central Afar. All major graben-bounding faults in central Afar cut through the Quaternary Afar Stratoid Series, constraining the ages of initiation to younger than approximately 1 million years. The Afar Stratoid Series presents a widespread thick lava unit that smoothed and infilled the topography created by the prior 30 Ma of deformation. The ∼1 Ma top surface of the Afar Stratoid Series serves as a regional marker predating most Quaternary deformation. The initiation of a particular fault (or location along a graben-bounding fault) begins a process of local lowering of base level – this is recorded in the landscape as a knickpoint that incises the paleo drainage in the footwalls of the grabens, connecting it to the basin floor with a deep canyon. The rate at which knickpoints retreat upstream permits a space-for-time approach to assessing initiation. This study develops and applies a landscape evolution model to estimate relative maturity and age of drainage basins by measuring the long profile response to tectonic forcings over intermediate timescales in central Afar. This is applied under the assumption that faulting postdates Afar Stratoid Series emplacement. The model is calibrated by analysis of stream longitudinal profiles and terrestrial cosmogenic nuclide (TCN) analysis of samples collected in the field. Longitudinal stream valley profiles are mapped through remote sensing and photogrammetry, and TCN concentrations are assessed to calculate the uplift rate along the graben-bounding faults. Results indicate that faulting of Dobe Graben in the central Afar initiated in the southeast and propagated to the northwest. Along the southern margin of Dobe Graben, ages of 503 ± 116 ka were calculated for fault initiation at Hann Canyon in the southeast, and 339 ± 80 ka at Ishyilu Canyon to the northwest. These timing constraints have implications for competing tectonics models for central Afar. The general approach developed for modeling knickpoint retreat has broader utility beyond Afar, allowing for application to extensional tectonic regimes elsewhere where similar geological and surface constraints exist.

Bedrock rivers are steep but not narrow: Hydrological and lithological controls on river geometry across the USA

Bedrock rivers are commonly expected to have steeper and narrower channels than alluvial rivers. However, understanding of bedrock river characteristics has largely been based on small samples of sites in specific climates and upland locations. We provide the first systematic assessment of bedrock and alluvial river channel characteristics for 1274 sites across a broad climatic gradient. We assess whether the width, width-to-depth ratio, and slope of bedrock channels differ from those of alluvial channels and the extent to which these differences are correlated with drainage area, mean annual flow (QMAF), grain size, and lithology. We find that bedrock channels occur at all drainage areas. For the same drainage area, bedrock channels are wider and steeper than alluvial channels. They also have a higher mean annual precipitation and hence QMAF, which likely causes the increased width. After accounting for differences in QMAF, both bedrock and alluvial channels have similar hydraulic scaling. Lithology affects both types of channels in a similar way, with channels on sedimentary lithologies being wider and less steep compared to those on igneous-metamorphic lithologies. Overall, our findings raise new questions about the evolution of bedrock river channels and pave the way for more accurate landscape evolution modeling.