Search USGS⌕ Search

SEARCH · Search USGS

Results for “JGR Earth Surface”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 recordsLinked to original sources

The biophysical role of water and ice within permafrost nearing collapse: Insights from novel geophysical observations

The impact of permafrost thaw on hydrologic, thermal, and biotic processes remains uncertain, in part due to limitations in subsurface measurement capabilities. To better understand subsurface processes in thermokarst environments, we collocated geophysical and biogeochemical instruments along a thaw gradient between forested permafrost and collapse-scar bogs at the Alaska Peatland Experiment (APEX) site near Fairbanks, Alaska. Ambient seismic noise monitoring provided continuous high-temporal resolution measurements of water and ice saturation changes. Maps of seismic velocity change identified areas of large summertime velocity reductions nearest the youngest bog, indicating potential thaw and expansion at the bog margin. These results corresponded well with complementary borehole nuclear magnetic resonance measurements of unfrozen water content with depth, which showed permafrost soils nearest the bog edges contained the largest amount of unfrozen water along the study transect, up to 25% by volume. In situ measurements of methane within permafrost soils revealed high concentrations at these bog-edge locations, up to 30% soil gas. Supra-permafrost talik zones were observed at the bog margins, indicating talik formation and perennial liquid water may drive lateral bog expansion and enhanced permafrost carbon losses preceding thaw. Comparison of seismic monitoring with wintertime surface carbon dioxide fluxes revealed differential responses depending on time and proximity to the bogs, capturing the controlling influence of subsurface water and ice on microbial activity and surficial emissions. This study demonstrates a multidisciplinary approach for gaining new understanding of how subsurface physical properties influence greenhouse gas production, emissions, and thermokarst development.

Alaska↗

Multi-year tracing of spatial and temporal dynamics of post-fire aeolian sediment transport using rare earth elements provide insights into grassland management

Aeolian sediment transport occurs as a function of, and with feedback to ecosystem changes and disturbances. Many desert grasslands are undergoing rapid changes in vegetation, including the encroachment of woody plants, which alters fire regimes and in turn can change the spatial and temporal patterns of aeolian sediment transport. We investigated aeolian sediment transport and spatial distribution of sediment in the surface soil for 7 years following a prescribed fire using a multiple rare earth element (REE) tracer-based approach in a shrub-encroached desert grassland in the northern Chihuahuan desert. Results indicate that even though the aeolian horizontal sediment mass flux increased approximately three-fold in the first windy season in the burned areas compared to control areas, there were no significant differences after three windy seasons. The soil surface of bare microsites was the major contributor of aeolian sediments in unburned areas (87%), while the shrub microsites contributed the least (<2%) during the observation period. However, after the prescribed fire, the contribution of aeolian sediments from shrub microsites increased considerably (∼40%), indicating post-fire microsite-scale sediment redistribution. The findings of this study, which is the first to use multiple REE tracers for multi-year analysis of the spatial and temporal dynamics of aeolian sediment transport, illustrate how disturbance by prescribed fire can influence aeolian processes and alters dryland soil geomorphology in which distinct soils develop over time at very fine spatial scales of individual plants.

New Mexico↗

Elastic stress coupling between supraglacial lakes

Supraglacial lakes have been observed to drain within hours of each other, leading to the hypothesis that stress transmission following one drainage may be sufficient to induce hydro-fracture-driven drainages of other nearby lakes. However, available observations characterizing drainage-induced stress perturbations have been insufficient to evaluate this hypothesis. Here, we use ice-sheet surface-displacement observations from a dense global positioning system array deployed in the Greenland Ice Sheet ablation zone to investigate elastic stress transmission between three neighboring supraglacial lake basins. We find that drainage of a central lake can place neighboring basins in either tensional or compressional stress relative to their hydro-fracture scarp orientations, either promoting or inhibiting hydro-fracture initiation beneath those lakes. For two lakes located within our array that drain close in time, we identify tensional surface stresses caused by ice-sheet uplift due to basal-cavity opening as the physical explanation for these lakes' temporally clustered hydro-fracture-driven drainages and frequent triggering behavior. However, lake-drainage-induced stresses in the up-flowline direction remain low beyond the margins of the drained lakes. This short stress-coupling length scale is consistent with idealized lake-drainage scenarios for a range of lake volumes and ice-sheet thicknesses. Thus, on elastic timescales, our observations and idealized-model results support a stress-transmission hypothesis for inducing hydro-fracture-driven drainage of lakes located within the region of basal cavity opening produced by the initial drainage, but refute this hypothesis for distal lakes.

JGR Earth Surface↗

Physical controls on the hydrology of perennially ice-covered lakes, Taylor Valley, Antarctica (1996-2013)

The McMurdo Dry Valleys, Antarctica, are a polar desert populated with numerous closed-watershed, perennially ice-covered lakes primarily fed by glacial melt. Lake levels have varied by as much as 8 m since 1972 and are currently rising after a decade of decreasing. Precipitation falls as snow, so lake hydrology is dominated by energy available to melt glacier ice and to sublimate lake ice. To understand the energy and hydrologic controls on lake level changes and to explain the variability between neighboring lakes, only a few kilometers apart, we model the hydrology for the three largest lakes in Taylor Valley. We apply a physically based hydrological model that includes a surface energy balance model to estimate glacial melt and lake sublimation to constrain mass fluxes to and from the lakes. Results show that lake levels are very sensitive to small changes in glacier albedo, air temperature, and wind speed. We were able to balance the hydrologic budget in two watersheds using meltwater inflow and sublimation loss from the ice-covered lake alone. A third watershed, closest to the coast, required additional inflow beyond model uncertainties. We hypothesize a shallow groundwater system within the active layer, fed by dispersed snow patches, contributes 23% of the inflow to this watershed. The lakes are out of equilibrium with the current climate. If the climate of our study period (1996-2013) persists into the future, the lakes will reach equilibrium starting in 2300, with levels 2-17 m higher, depending on the lake, relative to the 2020 level.

JGR Earth Surface↗

Seasonality of retreat rate of a wave-exposed marsh edge

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.

JGR Earth Surface↗

Inundation of a barrier island (Chandeleur Islands, Louisiana, USA) during a hurricane: Observed water-level gradients and modeled seaward sand transport

Large geomorphic changes to barrier islands may occur during inundation, when storm surge exceeds island elevation. Inundation occurs episodically and under energetic conditions that make quantitative observations difficult. We measured water levels on both sides of a barrier island in the northern Chandeleur Islands during inundation by Hurricane Isaac. Wind patterns caused the water levels to slope from the bay side to the ocean side for much of the storm. Modeled geomorphic changes during the storm were very sensitive to the cross-island slopes imposed by water-level boundary conditions. Simulations with equal or landward sloping water levels produced the characteristic barrier island storm response of overwash deposits or displaced berms with smoother final topography. Simulations using the observed seaward sloping water levels produced cross-barrier channels and deposits of sand on the ocean side, consistent with poststorm observations. This sensitivity indicates that accurate water-level boundary conditions must be applied on both sides of a barrier to correctly represent the geomorphic response to inundation events. More broadly, the consequence of seaward transport is that it alters the relationship between storm intensity and volume of landward transport. Sand transported to the ocean side may move downdrift, or aid poststorm recovery by moving onto the beach face or closing recent breaches, but it does not contribute to island transgression or appear as an overwash deposit in the back-barrier stratigraphic record. The high vulnerability of the Chandeleur Islands allowed us to observe processes that are infrequent but may be important at other barrier islands.

Louisiana↗

The application of ensemble wave forcing to quantify uncertainty of shoreline change predictions

Reliable predictions and accompanying uncertainty estimates of coastal evolution on decadal to centennial time scales are increasingly sought. So far, most coastal change projections rely on a single, deterministic realization of the unknown future wave climate, often derived from a global climate model. Yet, deterministic projections do not account for the stochastic nature of future wave conditions across a variety of temporal scales (e.g., daily, weekly, seasonally, and interannually). Here, we present an ensemble Kalman filter shoreline change model to predict coastal erosion and uncertainty due to waves at a variety of time scales. We compare shoreline change projections, simulated with and without ensemble wave forcing conditions by applying ensemble wave time series produced by a computationally efficient statistical downscaling method. We demonstrate a sizable (site-dependent) increase in model uncertainty compared with the unrealistic case of model projections based on a single, deterministic realization (e.g., a single time series) of the wave forcing. We support model-derived uncertainty estimates with a novel mathematical analysis of ensembles of idealized process models. Here, the developed ensemble modeling approach is applied to a well-monitored beach in Tairua, New Zealand. However, the model and uncertainty quantification techniques derived here are generally applicable to a variety of coastal settings around the world.

JGR Earth Surface↗

Twelve-year dynamics and rainfall thresholds for alternating creep and rapid movement of the Hooskanaden landslide from integrating InSAR, pixel offset tracking, and borehole and hydrological measurements

The Hooskanaden landslide is a large (~600 m wide × 1,300 m long), deep (~30 – 45 m) slide located in southwestern Oregon. Since 1958, it has had five moderate/major movements that catastrophically damaged the intersecting U.S. Highway 101, along with persistent slow wet‐season movements and a long‐term accelerating trend due to coastal erosion. Multiple remote sensing approaches, borehole measurements, and hydrological observations have been integrated to interpret the motion behaviors of the slide. Pixel offset tracking of both Sentinel‐1 and Sentinel‐2 images was carried out to reconstruct the 3‐D displacement field of the 2019 major event, and the results agree well with field measurements. A 12‐year displacement history of the landslide from 2007 to 2019 has been retrieved by incorporating offsets from Light Detection and Ranging (LiDAR) digital elevation model (DEM) gradients and Interferometric Synthetic Aperture Radar (InSAR) processing of ALOS and Sentinel‐1 images. Comparisons with daily/hourly ground precipitation reveal that the motion dynamics are predominantly controlled by intensity and temporal pattern of rainfall. A new empirical threefold rainfall threshold was therefore proposed to forecast the dates for the moderate/major movements. This threshold relies upon antecedent water‐year and previous 3‐day and daily precipitation and was able to represent observed movement periods well. Adaptation of our threshold methodology could prove useful for other large, deep landslides for which temporal forecasting has long been generally intractable. The averaged characteristic hydraulic conductivity and diffusivity were estimated as 6.6 × 10 −6 m/s and 6.6 × 10 −4 m 2 /s, respectively, based on the time lags between rainfall pulses and slide accelerations. Hydrologic modeling using these parameters helps to explain the ability of the new rainfall threshold.

Oregon↗

Modeling of barrier breaching during Hurricanes Sandy and Matthew

Physical processes driving barrier island change during storms are important to understand to mitigate coastal hazards and to evaluate conceptual models for barrier evolution. Spatial variations in barrier island topography, landcover characteristics, and nearshore and back-barrier hydrodynamics can yield complex morphological change that requires models of increasing resolution and physical complexity to predict. Using the Coupled Ocean-Atmosphere-Wave-Sediment Transport (COAWST) modeling system, we investigated two barrier island breaches that occurred on Fire Island, NY during Hurricane Sandy (2012) and at Matanzas, FL during Hurricane Matthew (2016). The model employed a recently implemented infragravity (IG) wave driver to represent the important effects of IG waves on nearshore water levels and sediment transport. The model simulated breaching and other changes with good skill at both locations, resolving differences in the processes and evolution. The breach simulated at Fire Island was 250 m west of the observed breach, whereas the breach simulated at Matanzas was within 100 m of the observed breach. Implementation of the vegetation module of COAWST to allow three-dimensional drag over dune vegetation at Fire Island improved model skill by decreasing flows across the back-barrier, as opposed to varying bottom roughness that did not positively alter model response. Analysis of breach processes at Matanzas indicated that both far-field and local hydrodynamics influenced breach creation and evolution, including remotely generated waves and surge, but also surge propagation through back-barrier waterways. This work underscores the importance of resolving the complexity of nearshore and back-barrier systems when predicting barrier island change during extreme events.

Florida, New York↗

Hydrogeomorphic recovery and temporal changes in rainfall thresholds for debris flows following wildfire

Wildfire-induced changes to soil and vegetation promote runoff-generated debris flows in steep watersheds. Postfire debris flows are most commonly observed in steep watersheds during the first wet season following a wildfire, but it is unclear how long the elevated threat of debris flow persists and why debris-flow potential changes in recovering burned areas. This work quantifies how rainfall intensity-duration (ID) thresholds for debris-flow initiation change with time since burning and provides a mechanistic explanation for these changes. We constrained a hydrologic model using field and remotely sensed measurements of soil-infiltration capacity, vegetation cover, runoff, and debris-flow activity. We applied this model to estimate rainfall ID thresholds for debris-flow initiation within three burned areas in the southwestern United States over a postfire recovery period of three to four years. Modeling suggests ID thresholds are lowest immediately following the fire (below a one-year recurrence interval [RI] storm) and increase with time, such that a 10- to 25-year RI storm would be required to generate a debris flow after three years of recovery. Modeled changes in rainfall ID thresholds result from increases in soil infiltration capacity, canopy interception, hydraulic roughness, and median grain size of sediment entrained in an incipient debris flow. The relative importance of each of these factors varied among our three sites. Results improve our ability to assess temporal changes in postfire debris-flow potential, highlight how site-specific factors may alter the persistence of postfire debris-flow hazards, and provide additional constraints on the timescale of recovery following wildfire.

Arizona, California, New Mexico↗

Sound-side inundation and seaward erosion of a barrier island during hurricane landfall

Barrier islands are especially vulnerable to hurricanes and other large storms, owing to their mobile composition, low elevations, and detachment from the mainland. Conceptual models of barrier-island evolution emphasize ocean-side processes that drive landward migration through overwash, inlet migration, and aeolian transport. In contrast, we found that the impact of Hurricane Dorian (2019) on North Core Banks, a 36-km barrier island on the Outer Banks of North Carolina, was primarily driven by inundation of the island from Pamlico Sound, as evidenced by storm-surge model results and observations of high-water marks and wrack lines. Analysis of photogrammetry products from aerial imagery collected before and after the storm indicate the loss of about 18% of the subaerial volume of the island through the formation of over 80 erosional washout channels extending from the marsh and washover platform, through gaps in the foredunes, to the shoreline. The washout channels were largely co-located with washover fans deposited by earlier events. Net seaward export of sediment resulted in the formation of deltaic bars offshore of the channels, which became part of the post-storm berm recovery by onshore bar migration and partial filling of the washouts with washover deposits within 2 months. This event represents a volumetric setback in the overwash/rollover behavior required for barrier transgression, but the new ponds and lowland habitats may provide beneficial habit for endangered species and will likely persist for years.

North Carolina↗

A model integrating satellite-derived shoreline observations for predicting fine-scale shoreline response to waves and sea-level rise across large coastal regions

Satellite-derived shoreline observations combined with dynamic shoreline models enable fine-scale predictions of coastal change across large spatiotemporal scales. Here, we present a satellite-data-assimilated, “littoral-cell”-based, ensemble Kalman-filter shoreline model to predict coastal change and uncertainty due to waves, sea-level rise (SLR), and other natural and anthropogenic processes. We apply the developed ensemble model to the entire California coastline (approximately 1,760 km), much of which is sparsely monitored with traditional survey methods (e.g., Lidar/GPS). Water-level-corrected, satellite-derived shoreline observations (obtained from the CoastSat toolbox) offer a nearly unbiased representation of in situ surveyed shorelines (e.g., mean sea-level elevation contours) at Ocean Beach, San Francisco. We demonstrate that model calibration with satellite observations during a 20-year hindcast period (1995–2015) provides nearly equivalent model forecast accuracy during a validation period (2015–2020) compared to model calibration with monthly in situ observations at Ocean Beach. When comparing model-predicted shoreline positions to satellite-derived observations, the model achieves an accuracy of <10 m RMSE for nearly half of the entire California coastline for the validation period. The calibrated/validated model is then applied for multi-decadal simulations of shoreline change due to projected wave and sea-level conditions, while holding the model parameters fixed. By 2100, the model estimates that 24%–75% of California's beaches may become completely eroded due to SLR scenarios of 1.0–3.0 m, respectively. The satellite-data-assimilated modeling system presented here is generally applicable to a variety of coastal settings around the world owing to the global coverage of satellite imagery.

JGR Earth Surface↗

Debris-flow process controls on steepland morphology in the San Gabriel Mountains, California

Steep landscapes evolve largely by debris flows, in addition to fluvial and hillslope processes. Abundant field observations document that debris flows incise valley bottoms and transport substantial sediment volumes, yet their contributions to steepland morphology remain uncertain. This has, in turn, limited the development of debris-flow incision rate formulations that produce morphology consistent with natural landscapes. In many landscapes, including the San Gabriel Mountains (SGM), California, steady-state fluvial channel longitudinal profiles are concave-up and exhibit a power-law relationship between channel slope and drainage area. At low drainage areas, however, valley slopes become nearly constant. These topographic forms result in a characteristically curved slope-area signature in log-log space. Here, we use a one-dimensional landform evolution model that incorporates debris-flow erosion to reproduce the relationship between this curved slope-area signature and erosion rate in the SGM. Topographic analysis indicates that the drainage area at which steepland valleys transition to fluvial channels correlates with measured erosion rates in the SGM, and our model results reproduce these relationships. Further, the model only produces realistic valley profiles when parameters that dictate the relationship between debris-flow erosion, valley-bottom slope, and debris-flow depth are within a narrow range. This result helps place constraints on the mathematical form of a debris-flow incision law. Finally, modeled fluvial incision outpaces debris-flow erosion at drainage areas less than those at which valleys morphologically transition from near-invariant slopes to concave profiles. This result emphasizes the critical role of debris-flow incision for setting steepland form, even as fluvial incision becomes the dominant incisional process.

California↗

A large sediment accretion wave along a northern California littoral cell

The northern California littoral cell of the Klamath River, which is a mixed rocky and sandy system with significant shoreline curvature, was investigated by examining ∼40 yr of satellite-derived shoreline positions and historical records. We find that an accretion wave of sediment was initiated near the Klamath River mouth in the late 1980s and translated downcoast over the subsequent decades. The wave passed rapidly (∼2,500 m/yr) through a rocky coastal reach with more oblique wave directions and slowly through a sandy reach (∼200 m/yr) where wave crests approach at more normal angles. Within the sandy reach, the accretion wave extended over 200 m offshore, was ∼10 km long, incorporated 20 ± 6 million m 3 of sediment, and averaged 1.3 ± 0.4 million m 3 /yr of longshore sediment transport over a 20-yr interval. Diffusion of the accretion wave was observed, but the diffusivity coefficient ( ε obs ∼0.01 m 2 /s) was lower than values predicted by theory, which we attribute to net sediment transport convergence in the study area caused by the curvature of the shoreline. Examining historical records, we find that increased sediment discharge in the Klamath River occurred during the 20th century from industrial-scale logging and climatic extremes. Thus, we hypothesize that increased river sediment discharge introduced new sediment to the littoral cell that initiated the observed accretion wave. These hypotheses can be tested with stratigraphic and mineralogic investigations of the broad study area beach that has formed during the past 150 years.

California↗

Testing megathrust rupture models using tsunami deposits

The 26 January 1700 CE Cascadia subduction zone earthquake ruptured much of the plate boundary and generated a tsunami that deposited sand in coastal marshes from northern California to Vancouver Island. Although the depositional record of tsunami inundation is extensive in some of these marshes, few sites have been investigated in enough detail to map the inland extent of sand deposition and depict variability in tsunami deposit thickness and grain size. We collected 129 cores in marshes of the Salmon River estuary in Oregon and reanalyzed 114 core logs from a 1987–88 study that mapped the inland extent of circa 1700 CE sandy tsunami deposits. The ca. 1700 CE tsunami deposit in the Salmon River estuary is easily recognized in cores ≤1 m deep in which a buried marsh peat is overlain by a well sorted sand bed with a sharp lower contact that thins and fines inland. We use tsunami deposit data and models of sandy tsunami sediment transport (using Delft3D-FLOW) to test 15 rupture models that could represent a ca. 1700 CE earthquake. At least 12–16 m of slip offshore of the Salmon River, which results in 0.8–1.0 m of coastal coseismic subsidence, is required to match the ca. 1700 CE sand deposit's inland extent, which is consistent with models of heterogeneous megathrust slip in ca. 1700 CE. Our methods of detailed tsunami deposit mapping, combined with sediment transport modeling, can be used to test models of megathrust ruptures and their tsunamis to potentially improve earthquake and tsunami hazard assessments.

JGR Earth Surface↗

Rainfall thresholds for postfire debris-flow initiation vary with short-duration rainfall climatology

The size, frequency, and geographic scope of severe wildfires are expanding across the globe, including in the Western United States. Recently burned steeplands have an increased likelihood of debris flows, which pose hazards to downstream communities. The conditions for postfire debris-flow initiation are commonly expressed as rainfall intensity-duration thresholds, which can be estimated given sufficient observational history. However, the spread of wildfire across diverse climates poses a challenge for accurate threshold prediction in areas with limited observations. Studies of mass-movement processes in unburned areas indicate that thresholds vary with local climate, such that higher rainfall rates are required for initiation in climates characterized by frequent intense rainfall. Here, we use three independent methods to test whether initiation of postfire runoff-generated debris flows across the Western United States varies similarly with climate. Through the compilation of observed thresholds at various fires, analysis of the spatial density of observed debris flows, and quantification of feature importance at different spatial scales, we show that postfire debris-flow initiation thresholds vary systematically with short-duration rainfall-intensity climatology. The predictive power of climatological data sets that are readily available before a fire occurs offers a much-needed tool for hazard management in regions that are facing increased wildfire activity, have sparse observational history, and/or have limited resources for field-based hazard assessment. Furthermore, if the observed variation in thresholds reflects long-term adjustment of the landscape to local climate, rapid shifts in rainfall intensity related to climate change will likely induce spatially variable shifts in postfire debris-flow likelihood.

Arizona, California, Colorado, Nevada, New Mexico,↗

Shoreline seasonality of California’s beaches

We report on remote sensing techniques developed to characterize seasonal shoreline cycles from satellite-derived shoreline measurements. These techniques are applied to 22-yr of shoreline measurements for over 777 km of beach along California's 1,700-km coast, for which the general understanding is that shorelines exhibit winter-narrow and summer-recovery seasonality. We find that approximately 90% of beach transects exhibit significant and recurring seasonal cycles in the shoreline position. Seasonal shoreline excursions are twice as large in northern and central California (17.5–32.2 m) than southern California (7.3–15.9 m; interquartile ranges). Clustering analyses were effective at characterizing the temporal patterns of the seasonality, revealing that ∼459 km of beach (59%) exhibit winter-narrow conditions, whereas ∼189 km (24%) and ∼50 km (6.4%) exhibit spring-narrow and summer-narrow conditions, respectively. These spring- and summer-narrow conditions are most common in southern California, where they represent over half of the total length of beach shoreline. Multivariate analyses reveal that wave climate and geomorphic setting are significantly related to the magnitude and timing of shoreline seasonal cycles. Combinations of these variables explain 44% of the seasonality variance of the complete data set and 85% of the variance for a subset of 93 long (>1 km) continuous beaches. We conclude that diversity in waves and geomorphic setting along California cause a broad range of seasonal patterns in the shoreline. Combined, this indicates that the overly generalized “winter-narrow/summer-recovery” conventions for California beaches are not expressed universally and that shoreline seasonality is far more diverse than these simple canonical rules.

California↗

Sensitivity analysis of a dynamic vegetation-sediment transport model using equadratures: Exploring inorganic accretion on a marsh platform

Salt marsh systems require a net import of inorganic sediment to maintain their structure in response to sea‐level rise. Marshes are affected by physical processes including tides, waves, sediment transport, and the influence of vegetation, and these processes interact in complex ways leading to sediment accretion or erosion. We implement a 3‐D hydrodynamic sediment transport model in an idealized marsh‐bay complex with a gently sloping edge, and use it as a laboratory to explore the processes leading to bed elevation change through the bay‐marsh continuum. We use the novel equadratures method for efficient sensitivity analysis to test the roles of wave, vegetation, and sediment parameters on wave dissipation, bed shear stress, sediment fluxes, and deposition and erosion across a transect spanning bay shallows to the marsh. Within the explored bounds of parameter uncertainty, significant wave height ( H sig ), settling velocity ( w s ), and critical shear stress (τ crit ) most strongly affect accretion on the marsh platform. Deposition is affected more by parameter‐parameter interactions, that is, both τ crit and w s or both H sig and w s , than by a single parameter varying alone. The sediment that accretes on the marsh platform originates beyond the marsh edge, indicating that the dynamics of the adjacent mudflat are important for predicting the fate of the marsh. Applying efficient sensitivity analysis techniques can empower process‐based models to test more parameters, larger ranges, and longer timeframes, enabling future predictions of marsh response to sea‐level rise based on physical processes

California↗