Search for: All records

We develop a robust and simple rule‐based algorithm to autonomously simulate alluvial fan deposition and evolution under continuously developing landscape conditions without prescribing deposition locations or imposing topographic constraints. Augmented with this algorithm, landscape evolution models are capable of dynamically detecting locations of potential fan deposition by statistical measures of surface topography and fluvial dynamics, then depositing fan sediments where and when the developed conditions require. To assess the method's efficacy in depositing sediment at a mountain‐valley transition zone characterized by a transport surface that permits unobstructed exit of sediment and water, a hypothetical scenario is created that involves a frontal, normal fault. It is followed by a series of sensitivity analyses to ascertain the influence of parameters affecting fan deposition and secondary processes. Uplift (u) and precipitation significantly impact fan morphological characteristics, which are within the range of real‐world fans. Higher rates of each cause the notable expansion of the fan area except in cases of exceptionally high precipitation rates. Fan area has a power‐law relationship with most of the tested parameters, , where is erodibility (lithology), and are fluvial parameters, and is catchment area ( ~0.9). This study is the first showcasing fan power‐law relationships using numerical modelling. While fan area increases with precipitation, there exists a threshold beyond which fan area diminishes, and the formation of fans ceases altogether. The algorithm provides a basis for improving mechanistic understanding of fans by offering a robust platform for testing process dominance and scaling. The results demonstrate its applicability for landscape evolution simulation over a long time and broad spatial scales. We also investigate the hydrological significance of including autonomously generated alluvial fans in coupled landscape evolution—hydrology models that focus on groundwater as well as surface water hydrology.

Fault zones are an important control on fluid flow, affecting groundwater supply, contaminant migration, and carbon storage. However, most models of fault seal do not consider fault zone cementation, despite the recognition that it is common and can dramatically reduce permeability. In order to study the field‐scale hydrogeologic effects of fault zone cementation, we conducted a series of aquifer pumping tests in wells installed within tens of meters of the variably cemented Loma Blanca Fault, a normal fault in the Rio Grande Rift. In the southern half of the study area, the fault zone is cemented by calcite; the cemented zone is 2‐8 m wide. In the center of the study area, the cemented fault zone is truncated at a buttress unconformity that laterally separates hydrostratigraphic units with a ∼40X difference in permeability. The fault zone north of the unconformity is not cemented. Constant rate pumping tests indicate that where the fault is cemented, it is a barrier to groundwater flow. This is an important demonstration that a fault with no clay in its core and similar sediment on both sides can be a barrier to groundwater flow by virtue of its cementation; most conceptual models for the hydrogeology of faults would predict that it would not be a barrier to groundwater flow. Additionally, the lateral permeability heterogeneity across the unconformity imposes another important control on the local flow field. This permeability discontinuity acts as either a no‐flow boundary or a constant head boundary, depending on the location of pumping.

Physics‐based distributed hydrological models that include groundwater are widely used to understand and predict physical and biogeochemical processes within watersheds. Typically, due to computational limitations, watershed modelers minimize the number of elements used in domain discretization, smoothing or even ignoring critical topographic features. We use an idealized model to investigate the implications of mesh refinement along streams and ridges for modeling three‐dimensional groundwater flow and transport in mountainous watersheds. For varying degrees of topographic complexity level (TCL), which increases with the level of mesh refinement, and geological heterogeneity, we estimate and compare steady state baseflow discharge, mean age, and concentration of subsurface weathering products. Results show that ignoring lower‐order streams or ridges diminishes flow through local flow paths and biases higher the contribution of intermediate and regional flow paths, and biases baseflow older. The magnitude of the bias increases for systems where permeability rapidly decreases with depth and is dominated by shallow flow paths. Based on a simple geochemical model, the concentration of weathering products is less sensitive to the TCL, partially due to the thermodynamic constraints on chemical reactions. Our idealized model also reproduces the observed emergent scaling relationship between the groundwater contribution to streamflow and drainage area, and finds that this scaling relationship is not sensitive to mesh TCL. The bias effects have important implications for the use of hydrological models in the interpretation of environmental tracer data and the prediction of biogeochemical evolution of stream water in mountainous watersheds.