Geology topics

Judson W. Harvey

Publications and source records attributed to Judson W. Harvey.

At least 19 recordsLinked to original sources

USGS2024

High-Flow Experimental Outcomes to Inform Everglades Restoration, 2010–22

The Decompartmentalization Physical Model (DPM) was an experimental facility in the central Everglades operated between 2010 and 2022 to release high flows through a levee-enclosed area of degraded ridge and slough wetland that had been isolated from flow for sixty years. The purpose of DPM experimental program was to make measurements before, during, and after seasonal high-flow releases that could help guide the Congressionally authorized Everglades restoration project known as the Decompartmentalization and Sheet Flow Enhancement Project. The DPM facility was operated by the South Florida Water Management District, with the U.S. Geological Survey (USGS) and several universities participating in experimental design and leading aspects of the research. The USGS research at DPM focused on measuring high-flow hydraulics and its sedimentary and ecological responses in downstream wetlands. USGS investigated interactions between flow and vegetation and microtopography that influenced flow velocity and water depth, bed shear stress, sediment entrainment, and the resulting downstream transport of suspended sediment and fate of particle-associated phosphorus. USGS also investigated high-flow changes in water-column mixing and gas exchange and resulting effects on metabolism of the aquatic ecosystem (primary productivity and respiration). USGS also investigated effects of built structures such as levee gaps that were constructed to reconnect levee-enclosed basins. This report describes the methods and results of the USGS-led data collection at DPM. The USGS studies at DPM have identified factors that influence effectiveness of restoration, specifically how high-flow releases maximize sheet flow and affect sediment and nutrient dynamics while minimizing undesirable outcomes caused by past management that bypassed wetlands by conveying polluted water through canals to ecologically sensitive downstream areas. The DPM high-flow experiments reconnected the Water Conservation Area 3A and Water Conservation Area 3B basins, and it therefore has become a central feature of the restoration’s Decompartmentalization and Sheet Flow Enhancement Project. DPM’s scientific findings have already influenced the adaptive management of Everglades restoration in guiding elements of the final design and implementation of the Central Everglades Planning Project-South. In addition to serving Everglades restoration, the DPM has the potential to influence similar adaptive management programs throughout the nation’s network of federal and state-managed river corridors, floodplains, and riparian ecosystems.

Florida↗

USGS2019

Residence time controls on the fate of nitrogen in flow‐through lakebed sediments

For many glacial lakes with highly permeable sediments, water exchange rates control hydrologic residence times within the sediment‐water interface (SWI) and the removal of reactive compounds such as nitrate, a common pollutant in lakes and groundwater. Here we conducted a series of focused tracer injection experiments in the upper 20 cm of the naturally downwelling SWI in a flow‐through lake on Cape Cod, MA. We systematically varied residence time and reactant controls on nitrate processing, using isotopically labeled 15 N nitrate to monitor the effect of these changes on nitrate removal via denitrification. The addition of acetate, a labile carbon compound, triggered the lake SWI to switch from net production to net removal of nitrate. When acetate was combined with increased residence time created by controlled reductions in water flux, we observed a fivefold increase in nitrate removal, a 26‐fold increase in N 2 production, and a 42‐fold increase in N 2 O production. We demonstrate that water residence time is an important control on the fate of nitrate in these lake SWIs and illustrate that seasonal conditions that alter lake exchange rates and variability in lake carbon may predict dynamic nitrate removal across the SWI. Additionally, observed N 2 O production during the oxic pore water experiments paired with geophysical characterization of the sediment porosity revealed that the lake SWI has less mobile pores occupying upward of 50% of the total porosity volume, which function as reactive microzones for nitrate processing.

Journal of Geophysical Research: Biogeosciences↗

USGS2018

The metabolic regimes of 356 rivers in the United States

A national-scale quantification of metabolic energy flow in streams and rivers can improve understanding of the temporal dynamics of in-stream activity, links between energy cycling and ecosystem services, and the effects of human activities on aquatic metabolism. The two dominant terms in aquatic metabolism, gross primary production (GPP) and aerobic respiration (ER), have recently become practical to estimate for many sites due to improved modeling approaches and the availability of requisite model inputs in public datasets. We assembled inputs from the U.S. Geological Survey and National Aeronautics and Space Administration for October 2007 to January 2017. We then ran models to estimate daily GPP, ER, and the gas exchange rate coefficient for 356 streams and rivers across the continental United States. We also gathered potential explanatory variables and spatial information for cross-referencing this dataset with other datasets of watershed characteristics. This dataset offers a first national assessment of many-day time series of metabolic rates for up to 9 years per site, with a total of 490,907 site-days of estimates.

Scientific Data↗

USGS2018

Hydrogeochemical controls on brook trout spawning habitats in a coastal stream

Brook trout ( Salvelinus fontinalis ) spawn in fall and overwintering egg development can benefit from stable, relatively warm temperatures in groundwater-seepage zones. However, eggs are also sensitive to dissolved oxygen concentration, which may be reduced in discharging groundwater (i.e., seepage). We investigated a 2 km reach of the coastal Quashnet River in Cape Cod, Massachusetts, USA, to relate preferred fish spawning habitats to geology, geomorphology, and discharging groundwater geochemistry. Thermal reconnaissance methods were used to locate zones of rapid groundwater discharge, which were predominantly found along the central channel of a wider stream valley section. Pore-water chemistry and temporal vertical groundwater flux were measured at a subset of these zones during field campaigns over several seasons. Seepage zones in open-valley sub-reaches generally showed suboxic conditions and higher dissolved solutes compared to the underlying glacial outwash aquifer. These discharge zones were cross-referenced with preferred brook trout redds and evaluated during 10 years of observation, all of which were associated with discrete alcove features in steep cutbanks, where stream meander bends intersect the glacial valley walls. Seepage in these repeat spawning zones was generally stronger and more variable than in open-valley sites, with higher dissolved oxygen and reduced solute concentrations. The combined evidence indicates that regional groundwater discharge along the broader valley bottom is predominantly suboxic due to the influence of near-stream organic deposits; trout show no obvious preference for these zones when spawning. However, the meander bends that cut into sandy deposits near the valley walls generate strong oxic seepage zones that are utilized routinely for redd construction and the overwintering of trout eggs. Stable water isotopic data support the conclusion that repeat spawning zones are located directly on preferential discharges of more localized groundwater. In similar coastal systems with extensive valley peat deposits, the specific use of groundwater-discharge points by brook trout may be limited to morphologies such as cutbanks, where groundwater flow paths do not encounter substantial buried organic material and remain oxygen-rich.

Hydrology and Earth System Sciences↗

USGS2018

Direct observations of hydrologic exchange occurring with less‐mobile porosity and the development of anoxic microzones in sandy lakebed sediments

Quantifying coupled mobile/less‐mobile porosity dynamics is critical to the prediction of biogeochemical storage, release, and transformation processes in the zone where groundwater and surface water exchange. The recent development of fine‐scale geoelectrical monitoring paired with pore‐water sampling in groundwater systems enables direct characterization of hydrologic exchange between more‐ and less‐mobile porosity during tracer tests. We adapt this technique to sandy interface sediments at a groundwater flow‐through kettle lake. Tracer experiments were conducted within controlled‐head permeameters over a range of specified downward flow conditions over several days. Although the bed was predominantly composed of highly permeable sands and gravels, cobble inclusions created less‐mobile flow zones at the centimeter scale. Less‐mobile porosity fractions, residence times, and rates of exchange were inferred from paired bulk and fluid electrical conductivity data, without the need for inverse model calibration. The conservative solute experiments were paired with 15 NO 3 − and other reactive amendments, revealing anaerobic processes occurring at shallow sediment depths where pore‐water sampling indicated bulk‐oxic conditions. The average less‐mobile porosity residence times as evaluated with the geoelectrical method were on 1‐hr timescales, which appear to be biogeochemically important in the context of creating anoxic microzones within less‐mobile porosity of sandy interface sediments.

Water Resources Research↗

USGS2018

On factors influencing air-water gas exchange in emergent wetlands

Knowledge of gas exchange in wetlands is important in order to determine fluxes of climatically and biogeochemically important trace gases and to conduct mass balances for metabolism studies. Very few studies have been conducted to quantify gas transfer velocities in wetlands, and many wind speed/gas exchange parameterizations used in oceanographic or limnological settings are inappropriate under conditions found in wetlands. Here six measurements of gas transfer velocities are made with SF 6 tracer release experiments in three different years in the Everglades, a subtropical peatland with surface water flowing through emergent vegetation. The experiments were conducted under different flow conditions and with different amounts of emergent vegetation to determine the influence of wind, rain, water flow, waterside thermal convection, and vegetation on air-water gas exchange in wetlands. Measured gas transfer velocities under the different conditions ranged from 1.1 cm h −1 during baseline conditions to 3.2 cm h −1 when rain and water flow rates were high. Commonly used wind speed/gas exchange relationships would overestimate the gas transfer velocity by a factor of 1.2 to 6.8. Gas exchange due to thermal convection was relatively constant and accounted for 14 to 51% of the total measured gas exchange. Differences in rain and water flow among the different years were responsible for the variability in gas exchange, with flow accounting for 37 to 77% of the gas exchange, and rain responsible for up to 40%.

Florida↗

USGS2017

Tracer-based characterization of hyporheic exchange and benthic biolayers in streams

Shallow benthic biolayers at the top of the streambed are believed to be places of enhanced biogeochemical turnover within the hyporheic zone. They can be investigated by reactive stream tracer tests with tracer recordings in the streambed and in the stream channel. Common in-stream measurements of such reactive tracers cannot localize where the processing primarily takes place, whereas isolated vertical depth profiles of solutes within the hyporheic zone are usually not representative of the entire stream. We present results of a tracer test where we injected the conservative tracer bromide together with the reactive tracer resazurin into a third-order stream and combined the recording of in-stream breakthrough curves with multidepth sampling of the hyporheic zone at several locations. The transformation of resazurin was used as an indicator of metabolism, and high-reactivity zones were identified from depth profiles. The results from our subsurface analysis indicate that the potential for tracer transformation (i.e., the reaction rate constant) varied with depth in the hyporheic zone. This highlights the importance of the benthic biolayer, which we found to be on average 2 cm thick in this study, ranging from one third to one half of the full depth of the hyporheic zone. The reach-scale approach integrated the effects of processes along the reach length, isolating hyporheic processes relevant for whole-stream chemistry and estimating effective reaction rates.

Water Resources Research↗

USGS2015

Denitrification in the Mississippi River network controlled by flow through river bedforms

Increasing nitrogen concentrations in the world’s major rivers have led to over-fertilization of sensitive downstream waters. Flow through channel bed and bank sediments acts to remove riverine nitrogen through microbe-mediated denitrification reactions. However, little is understood about where in the channel network this biophysical process is most efficient, why certain channels are more effective nitrogen reactors, and how management practices can enhance the removal of nitrogen in regions where water circulates through sediment and mixes with groundwater - hyporheic zones. Here we present numerical simulations of hyporheic flow and denitrification throughout the Mississippi River network using a hydrogeomorphic model. We find that vertical exchange with sediments beneath the riverbed in hyporheic zones, driven by submerged bedforms, has denitrification potential that far exceeds lateral hyporheic exchange with sediments alongside river channels, driven by river bars and meandering banks. We propose that geomorphic differences along river corridors can explain why denitrification efficiency varies between basins in the Mississippi River network. Our findings suggest that promoting the development of permeable bedforms at the streambed - and thus vertical hyporheic exchange - would be more effective at enhancing river denitrification in large river basins than promoting lateral exchange through induced channel meandering.

Mississippi River Network↗

USGS2014

Dual-domain mass-transfer parameters from electrical hysteresis: Theory and analytical approach applied to laboratory, synthetic streambed, and groundwater experiments

Models of dual‐domain mass transfer (DDMT) are used to explain anomalous aquifer transport behavior such as the slow release of contamination and solute tracer tailing. Traditional tracer experiments to characterize DDMT are performed at the flow path scale (meters), which inherently incorporates heterogeneous exchange processes; hence, estimated “effective” parameters are sensitive to experimental design (i.e., duration and injection velocity). Recently, electrical geophysical methods have been used to aid in the inference of DDMT parameters because, unlike traditional fluid sampling, electrical methods can directly sense less‐mobile solute dynamics and can target specific points along subsurface flow paths. Here we propose an analytical framework for graphical parameter inference based on a simple petrophysical model explaining the hysteretic relation between measurements of bulk and fluid conductivity arising in the presence of DDMT at the local scale. Analysis is graphical and involves visual inspection of hysteresis patterns to (1) determine the size of paired mobile and less‐mobile porosities and (2) identify the exchange rate coefficient through simple curve fitting. We demonstrate the approach using laboratory column experimental data, synthetic streambed experimental data, and field tracer‐test data. Results from the analytical approach compare favorably with results from calibration of numerical models and also independent measurements of mobile and less‐mobile porosity. We show that localized electrical hysteresis patterns resulting from diffusive exchange are independent of injection velocity, indicating that repeatable parameters can be extracted under varied experimental designs, and these parameters represent the true intrinsic properties of specific volumes of porous media of aquifers and hyporheic zones.

Water Resources Research↗

USGS2014

Hyporheic flow and transport processes: mechanisms, models, and biogeochemical implications

Fifty years of hyporheic zone research have shown the important role played by the hyporheic zone as an interface between groundwater and surface waters. However, it is only in the last two decades that what began as an empirical science has become a mechanistic science devoted to modeling studies of the complex fluid dynamical and biogeochemical mechanisms occurring in the hyporheic zone. These efforts have led to the picture of surface-subsurface water interactions as regulators of the form and function of fluvial ecosystems. Rather than being isolated systems, surface water bodies continuously interact with the subsurface. Exploration of hyporheic zone processes has led to a new appreciation of their wide reaching consequences for water quality and stream ecology. Modern research aims toward a unified approach, in which processes occurring in the hyporheic zone are key elements for the appreciation, management, and restoration of the whole river environment. In this unifying context, this review summarizes results from modeling studies and field observations about flow and transport processes in the hyporheic zone and describes the theories proposed in hydrology and fluid dynamics developed to quantitatively model and predict the hyporheic transport of water, heat, and dissolved and suspended compounds from sediment grain scale up to the watershed scale. The implications of these processes for stream biogeochemistry and ecology are also discussed."

Reviews of Geophysics↗

USGS2014

Dynamic hyporheic exchange at intermediate timescales: testing the relative importance of evapotranspiration and flood pulses

Hyporheic fluxes influence ecological processes across a continuum of timescales. However, few studies have been able to characterize hyporheic fluxes and residence time distributions (RTDs) over timescales of days to years, during which evapotranspiration (ET) and seasonal flood pulses create unsteady forcing. Here we present a data-driven, particle-tracking piston model that characterizes hyporheic fluxes and RTDs based on measured vertical head differences. We used the model to test the relative influence of ET and seasonal flood pulses in the Everglades (FL, USA), in a manner applicable to other low-energy floodplains or broad, shallow streams. We found that over the multiyear timescale, flood pulses that drive relatively deep (∼1 m) flow paths had the dominant influence on hyporheic fluxes and residence times but that ET effects were discernible at shorter timescales (weeks to months) as a break in RTDs. Cumulative RTDs on either side of the break were generally well represented by lognormal functions, except for when ET was strong and none of the standard distributions applied to the shorter timescale. At the monthly timescale, ET increased hyporheic fluxes by 1–2 orders of magnitude; it also decreased 6 year mean residence times by 53–87%. Long, slow flow paths driven by flood pulses increased 6 year hyporheic fluxes by another 1–2 orders of magnitude, to a level comparable to that induced over the short term by shear flow in streams. Results suggest that models of intermediate-timescale processes should include at least two-storage zones with different RTDs, and that supporting field data collection occur over 3–4 years.

Florida↗

USGS2013

Relative significance of microtopography and vegetation as controls on surface water flow on a low-gradient floodplain

Surface water flow controls water velocities, water depths, and residence times, and influences sediment and nutrient transport and other ecological processes in shallow aquatic systems. Flow through wetlands is substantially influenced by drag on vegetation stems but is also affected by microtopography. Our goal was to use microtopography data directly in a widely used wetland model while retaining the advantages of the model’s one-dimensional structure. The base simulation with no explicit treatment of microtopography only performed well for a period of high water when vegetation dominated flow resistance. Extended simulations using microtopography can improve the fit to low-water conditions substantially. The best fit simulation had a flow conductance parameter that decreased in value by 70 % during dry season such that mcrotopographic features blocked 40 % of the cross sectional width for flow. Modeled surface water became ponded and flow ceased when 85 % of the cross sectional width became blocked by microtopographic features. We conclude that vegetation drag dominates wetland flow resistance at higher water levels and microtopography dominates at low water levels with the threshold delineated by the top of microtopographic features. Our results support the practicality of predicting flow on floodplains using relatively easily measured physical and biological variables.

Florida↗

USGS2013

Hyporheic zone denitrification: controls on effective reaction depth and contribution to whole-stream mass balance

Stream denitrification is thought to be enhanced by hyporheic transport but there is little direct evidence from the field. To demonstrate at a field site, we injected 15 NO 3 − , Br (conservative tracer), and SF 6 (gas exchange tracer) and compared measured whole-stream denitrification with in situ hyporheic denitrification in shallow and deeper flow paths of contrasting geomorphic units. Hyporheic denitrification accounted for between 1 and 200% of whole-stream denitrification. The reaction rate constant was positively related to hyporheic exchange rate (greater substrate delivery), concentrations of substrates DOC and nitrate, microbial denitrifier abundance (nirS), and measures of granular surface area and presence of anoxic microzones. The dimensionless product of the reaction rate constant and hyporheic residence time, λ hz τ hz define a Damköhler number, Da den-hz that was optimal in the subset of hyporheic flow paths where Da den-hz ≈ 1. Optimal conditions exclude inefficient deep pathways transport where substrates are used up and also exclude inefficient shallow pathways that require repeated hyporheic entries and exits to complete the reaction. The whole-stream reaction significance, R s (dimensionless), was quantified by multiplying Da den-hz by the proportion of stream discharge passing through the hyporheic zone. Together these two dimensionless metrics, one flow-path scale and the other reach-scale, quantify the whole-stream significance of hyporheic denitrification. One consequence is that the effective zone of significant denitrification often differs from the full depth of the hyporheic zone, which is one reason why whole-stream denitrification rates have not previously been explained based on total hyporheic-zone metrics such as hyporheic-zone size or residence time.

Indiana↗

USGS2013

Interactions between hyporheic flow produced by stream meanders, bars, and dunes

Stream channel morphology from grain-scale roughness to large meanders drives hyporheic exchange flow. In practice, it is difficult to model hyporheic flow over the wide spectrum of topographic features typically found in rivers. As a result, many studies only characterize isolated exchange processes at a single spatial scale. In this work, we simulated hyporheic flows induced by a range of geomorphic features including meanders, bars and dunes in sand bed streams. Twenty cases were examined with 5 degrees of river meandering. Each meandering river model was run initially without any small topographic features. Models were run again after superimposing only bars and then only dunes, and then run a final time after including all scales of topographic features. This allowed us to investigate the relative importance and interactions between flows induced by different scales of topography. We found that dunes typically contributed more to hyporheic exchange than bars and meanders. Furthermore, our simulations show that the volume of water exchanged and the distributions of hyporheic residence times resulting from various scales of topographic features are close to, but not linearly additive. These findings can potentially be used to develop scaling laws for hyporheic flow that can be widely applied in streams and rivers.

Water Resources Research↗

USGS2012

Resolving hyporheic and groundwater components of streambed water flux

Hyporheic and groundwater fluxes typically occur together in permeable sediments beneath flowing stream water. However, streambed water fluxes quantified using the thermal method are usually interpreted as representing either groundwater or hyporheic fluxes. Our purpose was to improve understanding of co-occurring groundwater and hyporheic fluxes using streambed temperature measurements and analysis of one-dimensional heat transport in shallow streambeds. First, we examined how changes in hyporheic and groundwater fluxes affect their relative magnitudes by reevaluating previously published simulations. These indicated that flux magnitudes are largely independent until a threshold is crossed, past which hyporheic fluxes are diminished by much larger (1000-fold) groundwater fluxes. We tested accurate quantification of co-occurring fluxes using one-dimensional approaches that are appropriate for analyzing streambed temperature data collected at field sites. The thermal analytical method, which uses an analytical solution to the one-dimensional heat transport equation, was used to analyze results from a numerical heat transport model, in which hyporheic flow was represented as increased thermal dispersion at shallow depths. We found that co-occurring groundwater and hyporheic fluxes can be quantified in streambeds, although not always accurately. For example, using a temperature time series collected in a sandy streambed, we found that hyporheic and groundwater flow could both be detected when thermal dispersion due to hyporheic flow was significant compared to thermal conduction. We provide guidance for when thermal data can be used to quantify both hyporheic and groundwater fluxes, and we show that neglecting thermal dispersion may affect accuracy and interpretation of estimated streambed water fluxes.

Water Resources Research↗

USGS2012

Thresholds of flow-induced bed disturbances and their effects on stream metabolism in an agricultural river

Storm-driven flow pulses in rivers destroy and restructure sediment habitats that affect stream metabolism. This study examined thresholds of bed disturbances that affected patch- and reach-scale sediment conditions and metabolism rates. A 4 year record of discharge and diel changes in dissolved oxygen concentrations (ΔDO) was analyzed for disturbances and recovery periods of the ΔDO signal. Disturbances to the ΔDO signal were associated with flow pulses, and the recovery times for the ΔDO signal were found to be in two categories: less than 5 days (30% of the disturbances) or greater than 15 days (70% of the disturbances). A field study was performed during the fall of 2007, which included a storm event that increased discharge from 3.1 to 6.9 m 3 /s over a 7 h period. During stable flow conditions before the storm, variability in patch-scale stream metabolism values were associated with sediment texture classes with values ranging from −16.4 to 2.3 g O 2 2 /d (negative sign indicates net respiration) that bounded the reach-averaged rate of −5.6 g O 2 2 /d. Hydraulic modeling of bed shear stresses demonstrated a storm-induced flow pulse mobilized approximately 25% of the bed and reach-scale metabolism rates shifted from −5 to −40 g O 2 2 /d. These results suggest that storm-induced bed disturbances led to threshold behavior with respect to stream metabolism. Small flow pulses resulted in partial-bed mobilization that disrupted stream metabolism by increased turbidity with short recovery times. Large flow pulses resulted in full-bed mobilization that disrupted stream metabolism by destroying periphyton habitats with long recovery times.

Water Resources Research↗

USGS2012

Directional connectivity in hydrology and ecology

Quantifying hydrologic and ecological connectivity has contributed to understanding transport and dispersal processes and assessing ecosystem degradation or restoration potential. However, there has been little synthesis across disciplines. The growing field of ecohydrology and recent recognition that loss of hydrologic connectivity is leading to a global decline in biodiversity underscore the need for a unified connectivity concept. One outstanding need is a way to quantify directional connectivity that is consistent, robust to variations in sampling, and transferable across scales or environmental settings. Understanding connectivity in a particular direction (e.g., streamwise, along or across gradient, between sources and sinks, along cardinal directions) provides critical information for predicting contaminant transport, planning conservation corridor design, and understanding how landscapes or hydroscapes respond to directional forces like wind or water flow. Here we synthesize progress on quantifying connectivity and develop a new strategy for evaluating directional connectivity that benefits from use of graph theory in ecology and percolation theory in hydrology. The directional connectivity index (DCI) is a graph-theory based, multiscale metric that is generalizable to a range of different structural and functional connectivity applications. It exhibits minimal sensitivity to image rotation or resolution within a given range and responds intuitively to progressive, unidirectional change. Further, it is linearly related to the integral connectivity scale length—a metric common in hydrology that correlates well with actual fluxes—but is less computationally challenging and more readily comparable across different landscapes. Connectivity-orientation curves (i.e., directional connectivity computed over a range of headings) provide a quantitative, information-dense representation of environmental structure that can be used for comparison or detection of subtle differences in the physical-biological feedbacks driving pattern formation. Case-study application of the DCI to the Everglades in south Florida revealed that loss of directional hydrologic connectivity occurs more rapidly and is a more sensitive indicator of declining ecosystem function than other metrics (e.g., habitat area) used previously. Here and elsewhere, directional connectivity can provide insight into landscape drivers and processes, act as an early-warning indicator of environmental degradation, and serve as a planning tool or performance measure for conservation and restoration efforts.

Ecological Applications↗

USGS2011

Modeling of hydroecological feedbacks predicts distinct classes of landscape pattern, process, and restoration potential in shallow aquatic ecosystems

It is widely recognized that interactions between vegetation and flow cause the emergence of channel patterns that are distinct from the standard Schumm classification of river channels. Although landscape pattern is known to be linked to ecosystem services such as habitat provision, pollutant removal, and sustaining biodiversity, the mechanisms responsible for the development and stability of different landscape patterns in shallow, vegetated flows have remained poorly understood. Fortunately, recent advances have made possible large-scale models of flow through vegetated environments that can be run over a range of environmental variables and over timescales of millennia. We describe a new, quasi-3D cellular automata model that couples simulations of shallow-water flow, bed shear stresses, sediment transport, and vegetation dynamics in an efficient manner. That efficiency allowed us to apply the model widely in order to determine how different hydroecological feedbacks control landscape pattern and process in various types of wetlands and floodplains. Distinct classes of landscape pattern were uniquely associated with specific types of allogenic and autogenic drivers in wetland flows. Regular, anisotropically patterned wetlands were dominated by allogenic processes (i.e., processes driven by periodic high water levels and flow velocities that redistribute sediment), relative to autogenic processes (e.g., vegetation production, peat accretion, and gravitational erosion). These anistropically patterned wetlands are therefore particularly prone to hydrologic disturbance. Other classes of wetlands that emerged from simulated interactions included maze-patterned, amorphous, and topographically noisy marshes, open marsh with islands, banded string-pool sequences perpendicular to flow, parallel deep and narrow channels flanked by marsh, and ridge-and-slough patterned marsh oriented parallel to flow. Because vegetation both affects and responds to the balance between the transport capacity of the flow and sediment supply, these vegetated systems exhibit a feedback that is not dominant in most rivers. Consequently, unlike in most rivers, it is not possible to predict the “channel pattern” of a vegetated landscape based only on discharge characteristics and sediment supply; the antecedent vegetation pattern and vegetation dynamics must also be known. In general, the stability of different wetland pattern types is most strongly related to factors controlling the erosion and deposition of sediment at vegetation patch edges, the magnitude of sediment redistribution by flow, patch elevation relative to water level, and the variability of erosion rates in vegetation patches with low flow-resistance. As we exemplify in our case-study of the Everglades ridge and slough landscape, feedback between flow and vegetation also causes hysteresis in landscape evolution trajectories that will affect the potential for landscape restoration. Namely, even if the hydrologic conditions that historically produced higher flows are restored, degraded portions of the ridge and slough landscape are unlikely to revert to their former patterning. As wetlands and floodplains worldwide become increasingly threatened by climate change and urbanization, the greater mechanistic understanding of landscape pattern and process that our analysis provides will improve our ability to forecast and manage the behavior of these ecosystems.

Florida↗