Geology topics

Mark T. Stacey

Publications and source records attributed to Mark T. Stacey.

At least 19 recordsLinked to original sources

USGS2025

Comparison of creek and bay influences on salt marsh sediment budget and deposition patterns

The resilience of salt marshes with low organic production depends on their effective capture and retention of mineral sediment from adjacent waters. Little prior work has directly compared mechanisms of sediment import from wave-influenced marsh boundaries against those of tidal creeks. We used simultaneous deployment of net-deposition tiles and oceanographic sensors to identify the timing and magnitude of sediment import/export to, and redistribution within, a marsh in south San Francisco Bay. As the marsh has both an eroding bay-exposed scarp and a prominent tidal creek, we investigated the mechanisms and magnitudes of sediment import from the marsh-bay versus the marsh-creek interface. The strong daily sea breezes of the summer season produced most of the wave-driven erosion of the marsh scarp and controlled suspended sediment concentrations; the winter season had weaker winds punctuated by a few storms. A large seasonal difference in suspended sediment concentrations influenced both flood and ebb sediment fluxes to the marsh and led to much higher rates of import in the summer. Both bay-side and creek-side processes were important to total marsh sediment budget. Bay-side sediment contributions were more variable in time due to the bay-influenced environment, and creek-side contributions were overall larger, reflecting the large proportion of the marsh fed by creek water. Sediment was redistributed throughout the system, with erosion near the bay-edge, accretion near the creek-edge and slow import to the marsh interior. The marsh was net importing sediment in the summer and exporting in the winter from different rates of these processes; on an annual scale, the marsh was net importing despite rapid lateral marsh loss. These findings emphasize that a positive sediment budget does not imply a stable marsh and that both creek- and edge-side dynamics are important for marsh sedimentation and geomorphic trajectories. Further, we expand understandings of non-storm and seasonal controls on marsh sedimentation.

California↗

USGS2025

Mitigating flood risks in urban estuaries: Tidal dynamics, shoreline hardening, nature-based solutions, and floodgates in San Francisco Bay

Hydrodynamic models are valuable tools for understanding the primary factors influencing daily and peak water levels and for guiding discussions on potential adaptation strategies for managing flood risk in coastal areas. This analysis uses the Delft3D San Francisco Bay-Delta Community Model to simulate water levels and incorporates the effects of a number of adaptation measures in the urban San Francisco Bay estuary, California. In particular, we examine the influence of shoreline hardening, nature-based solutions, and subregional floodgates on regional water levels. The result shows that under present conditions, tidal amplification is responsible for generating a wide distribution of extreme water levels across San Francisco Bay. Tidal amplification is found to decrease under sea level rise, thereby producing a relative damping effect on extremes. A comparison of different shoreline scenarios demonstrates that hard frontal shorelines result in higher tidal amplification, whereas restored (soft) shorelines lower amplification. The current shoreline configuration has both hard and soft characteristics and results in an intermediate tidal response. In some areas, wetland restoration reduces extreme water levels by as much as 20 cm, whereas hard-shoreline addition elevates them by as much as 10 cm for 1.5 m of sea level rise. Furthermore, local floodgates can significantly reduce high water levels without major adverse effects elsewhere in San Francisco Bay. These findings point toward the justification for a range of adaptive measures across political boundaries, weighing hard and soft options in addressing the mounting danger of sea level rise.

California↗

USGS2024

The role of spring-neap phasing of intermittent lateral exchange in the ecosystem of a channel-shoal estuary

Lateral variability is a fundamental feature of channel-shoal estuaries, and exchanges between the channel and shoal can play an important role in the dynamics of the ecosystem in each region. This lateral exchange of biomass interacts with vertical structure and variability, particularly in the channel, to define algal biomass accumulation in the estuary. In this paper, we investigate how time-variable lateral exchange affects phytoplankton dynamics with a biophysical model that links two water columns via intermittent exchange. We find that time variability in the exchange influences biomass by increasing concentrations in the shoals and decreasing them in the channel when the time variability happens on a timescale greater than the timescales of biological processes, and the strength of the effect increases with the period of the intermittency. At timescales of variability comparable to the spring-neap cycle, however, the interplay between lateral exchange and the ecosystem response is complicated by the fortnightly development of stratification in the channel and the role that channel-shoal interaction plays in defining that stratification. As a result, for lateral exchange variability with periods of 7 and 14 days, the influence of the shoal ecosystem on the channel ecosystem is sensitive to the phasing of exchange relative to the spring-neap cycle, due to the fact that neap tide exchanges can create stratification events that are larger in magnitude and duration than would occur in the absence of lateral exchange, causing the channel to transition into net positive growth conditions. We conclude that lateral exchange influences the estuarine ecosystem both directly, through the exchange of biomass between shoals with net positive growth and adjoining channels and indirectly through its role in defining stratification events that allow the channel itself to have net positive growth.

Estuaries and Coasts↗

USGS2024

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↗

USGS2024

Seasonal and spatial patterns in sediment transport in an erosional marsh system

Salt marshes provide expansive ecosystem services to coastal regions, including serving as bird and fish habitat, providing flood protection, and sequestering carbon and toxins (Barbier et al. 2011). While marshes are naturally dynamic, they are undergoing net loss in the 21st century largely due to lateral erosion (Camp- bell et al. 2022). In regions where marsh migration is constrained by built infrastructure, or where sediment resources may be insufficient for vertical or lateral expansion, approaches to prevent net marsh loss against lateral erosion are limited. To better understand and strategize around these conditions, detailed studies of marsh-edge loss mechanisms and their connections to sediment transport dynamics are necessary (Smith et al. 2021).

California↗

USGS2023

Future marsh evolution due to tidal changes induced by human adaptation to sea level rise

With sea level rise threatening coastal development, decision-makers are beginning to act by modifying shorelines. Previous research has shown that hardening or softening shorelines may change the tidal range under future sea level rise. Tidal range can also be changed by natural factors. Coastal marshes, which humans increasingly depend on for shoreline protection, are ecologically sensitive to tidal range. Therefore, it is critical to examine how changes in tidal range could influence marsh processes. A marsh accretion model was used to investigate the ecological response of a San Francisco Bay, California, USA marsh to multiple tidal range scenarios and sea level rise from 2010 to 2100. The scenarios include a baseline scenario with no shoreline modifications in the estuary, a shoreline hardening scenario that amplifies the tidal range, and 14 tidal range scenarios as a sensitivity analysis that span tidal amplification and reduction of the baseline scenario. The modeling results expose key tradeoffs to consider when planning for sea level rise. Compared to the baseline, the hardening scenario shows minor differences. However, further tidal amplification prolongs marsh survival but decreases Sarcocornia pacifica cover, an important species for certain threatened wildlife and an effective attenuator of wave energy. Conversely, tidal reduction precipitates marsh drowning but shows gains in Sarcocornia pacifica cover. These mixed impacts of tidal amplification and reduction shown by the model indicate potential tradeoffs in relation to marsh survival, habitat characteristics, and shoreline protection. This study suggests the need for a cross-sectoral, regional approach to sea level rise adaptation.

California↗

USGS2019

Seasonal, spring-neap, and tidal variation in cohesive sediment transport parameters in estuarine shallows

Numerical models for predicting sediment concentrations and transport rely on parameters such as settling velocity and bed erodibility that describe sediment characteristics, yet these parameters are rarely probed directly. We investigated temporal and spatial variation in sediment parameters in the shallows of San Pablo Bay, CA. Flow, turbulence, and suspended sediment data were measured at sites located at 1 and 2 m below mean lower low water (MLLW) from November 2013 through April 2015, supplemented by monthlong periods in 2011, 2012, and 2016. Maximum current velocities were 0.40–0.47 m s -1 at these depths; the strongest currents decreased to 0.27–0.34 m s -1 during neap periods. Winters 2013–2014 and 2014–2015 experienced strong drought conditions, limiting the potential for seasonal impact on sediment conditions during this experiment. Despite this, the more storm‐influenced site showed clear changes during the winter: the roughness parameter decreased from 10 −4 to 10 −5 m, from hydrodynamically rough to smooth conditions, and bed erodibility increased by an order of magnitude. Median settling velocity was 2.05·10 −4 m s -1 ; it varied twofold within a tidal cycle, decreasing as current velocity grew during flood and ebb. This tidal control on floc size affected settling velocity on the spring‐neap timescale, possibly driving a spring‐neap oscillation in erodibility. Our findings highlight variation in sediment dynamics that is commonly ignored in numerical models and the need for field observations to ground truth ongoing modeling efforts.

California↗

USGS2018

Three-dimensional modeling of ﬁne sediment transport by waves and currents in a shallow estuary

A suspended sediment transport model is implemented in the unstructured‐grid SUNTANS model and applied to study fine‐grained sediment transport in South San Francisco Bay. The model enables calculation of suspension of bottom sediment based on combined forcing of tidal currents and wind waves. We show that accurate results can be obtained by employing two‐size classes which are representative of microflocs and macroflocs in the Bay. A key finding of the paper is that the critical calibration parameter is the ratio of the erosion of the microflocs to macroflocs from the bed. Different values of this erosion ratio are needed on the shallow shoals and deeper channels because of the different nature of the sediment dynamics in these regions. Application of a spatially variable erosion ratio and critical shear stress for erosion is shown to accurately reproduce observed suspended sediment concentration at four‐field sites located along a cross‐channel transect. The results reveal a stark contrast between the behavior of the suspended sediment concentration on the shoals and in the deep channel. Waves are shown to resuspend sediments on the shoals, although tidal and wind‐generated currents are needed to mix the thin wave‐driven suspensions into the water column. The contribution to the suspended sediment concentration in the channel by transport from the shoals is similar in magnitude to that due to local resuspension. However, the local contribution is in phase with strong bottom currents which resuspend the sediments, while the contribution from the shoals peaks during low‐water slack tide.

San Francisco Bay↗

USGS2018

The influence of sea level rise on the regional interdependence of coastal infrastructure

Sea level rise (SLR) is placing both immediate and long‐term pressures on coastal communities to take protective actions. Projects in the United States, and in many locations throughout the world, generally involve local jurisdictions raising the elevation of shoreline protection elements, with limited or no analysis of the feedback between shoreline management decisions and the impacts to water levels regionally. Our study examines the impact of local shoreline development on regional flood risk and considers SLR scenarios up to 1.5 m using a large‐scale numerical model, as an example, for San Francisco Bay. Here we show that measures to prevent flooding along an embayment shoreline in one location or subregion may increase inundation elsewhere in the system. The network of interactions occurs not only within subbasins of the Bay but also across the greater geographic extent from one end of the Bay to the other, and local jurisdiction may have either reciprocal relationships with or asymmetric impacts on one other. Importantly, the nature of the interaction network is seen to evolve with SLR: interactions are purely subregional at current sea level but with higher sea level (e.g., 1 m of SLR), not only do the subregional interdependencies strengthen but also regional interdependences emerge.

California↗

USGS2017

Interactions of estuarine shoreline infrastructure with multiscale sea level variability

Sea level rise increases the risk of storms and other short‐term water‐rise events, because it sets a higher water level such that coastal surges become more likely to overtop protections and cause floods. To protect coastal communities, it is necessary to understand the interaction among multiday and tidal sea level variabilities, coastal infrastructure, and sea level rise. We performed a series of numerical simulations for San Francisco Bay to examine two shoreline scenarios and a series of short‐term and long‐term sea level variations. The two shoreline configurations include the existing topography and a coherent full‐bay containment that follows the existing land boundary with an impermeable wall. The sea level variability consists of a half‐meter perturbation, with duration ranging from 2 days to permanent (i.e., sea level rise). The extent of coastal flooding was found to increase with the duration of the high‐water‐level event. The nonlinear interaction between these intermediate scale events and astronomical tidal forcing only contributes ∼1% of the tidal heights; at the same time, the tides are found to be a dominant factor in establishing the evolution and diffusion of multiday high water events. Establishing containment at existing shorelines can change the tidal height spectrum up to 5%, and the impact of this shoreline structure appears stronger in the low‐frequency range. To interpret the spatial and temporal variability at a wide range of frequencies, Optimal Dynamic Mode Decomposition is introduced to analyze the coastal processes and an inverse method is applied to determine the coefficients of a 1‐D diffusion wave model that quantify the impact of bottom roughness, tidal basin geometry, and shoreline configuration on the high water events

Journal of Geophysical Research: Oceans↗

USGS2017

Clusters of community exposure to coastal flooding hazards based on storm and sea level rise scenarios—implications for adaptation networks in the San Francisco Bay region

Sea level is projected to rise over the coming decades, further increasing the extent of flooding hazards in coastal communities. Efforts to address potential impacts from climate-driven coastal hazards have called for collaboration among communities to strengthen the application of best practices. However, communities currently lack practical tools for identifying potential partner communities based on similar hazard exposure characteristics. This study uses statistical cluster analysis to identify similarities in community exposure to flooding hazards for a suite of sea level rise and storm scenarios. We demonstrate this approach using 63 jurisdictions in the San Francisco Bay region of California (USA) and compare 21 distinct exposure variables related to residents, employees, and structures for six hazard scenario combinations of sea level rise and storms. Results indicate that cluster analysis can provide an effective mechanism for identifying community groupings. Cluster compositions changed based on the selected societal variables and sea level rise scenarios, suggesting that a community could participate in multiple networks to target specific issues or policy interventions. The proposed clustering approach can serve as a data-driven foundation to help communities identify other communities with similar adaptation challenges and to enhance regional efforts that aim to facilitate adaptation planning and investment prioritization.

California↗

USGS2016

Bivalve grazing can shape phytoplankton communities

The ability of bivalve filter feeders to limit phytoplankton biomass in shallow waters is well-documented, but the role of bivalves in shaping phytoplankton communities is not. The coupled effect of bivalve grazing at the sediment-water interface and sinking of phytoplankton cells to that bottom filtration zone could influence the relative biomass of sinking (diatoms) and non-sinking phytoplankton. Simulations with a pseudo-2D numerical model showed that benthic filter feeding can interact with sinking to alter diatom:non-diatom ratios. Cases with the smallest proportion of diatom biomass were those with the fastest sinking speeds and strongest bivalve grazing rates. Hydrodynamics modulated the coupled sinking-grazing influence on phytoplankton communities. For example, in simulations with persistent stratification, the non-sinking forms accumulated in the surface layer away from bottom grazers while the sinking forms dropped out of the surface layer toward bottom grazers. Tidal-scale stratification also influenced vertical gradients of the two groups in opposite ways. The model was applied to Suisun Bay, a low-salinity habitat of the San Francisco Bay system that was transformed by the introduction of the exotic clam Potamocorbula amurensis . Simulation results for this Bay were similar to (but more muted than) those for generic habitats, indicating that P. amurensis grazing could have caused a disproportionate loss of diatoms after its introduction. Our model simulations suggest bivalve grazing affects both phytoplankton biomass and community composition in shallow waters. We view these results as hypotheses to be tested with experiments and more complex modeling approaches.

California↗

USGS2015

Model-based interpretation of sediment concentration and vertical flux measurements in a shallow estuarine environment

A one-dimensional numerical model describing tidally varying vertical mixing and settling was used to interpret sediment concentrations and vertical fluxes observed in the shoals of South San Francisco Bay by two acoustic Doppler velocimeters (ADVs) at elevations of 0.36 m and 0.72 m above bed. Measured sediment concentrations changed by up to 100 g m −3 over the semidiurnal tidal cycle. These dynamics were dominated by local resuspension and settling. Multiple particle class models suggested the existence of a class with fast settling velocities ( w s of 9.0 × 10 −4 m s −1 in spring and 5.8 × 10 −4 m s −1 in fall) and a slowly settling particle fraction ( w s of <1 × 10 −7 m s −1 in spring and 1.4 × 10 −5 m s −1 in fall). Modeled concentrations of slowly settling particles at 0.36 m were as high as 20 g m −3 during fall and varied with the spring-neap cycle while fine sediment concentrations in spring were constant around 5 g m −3 . Analysis of in situ water column floc size distributions suggested that floc properties in the lower part of the water column were most likely governed by particle-size distribution on the bed and not by coagulation, validating our multiple particle size approach. A comparison of different sediment bed models with respect to model performance, sensitivity, and identifiability suggested that the use of a sediment erosion model linear in bottom shear stress τ b ( E = M ( τ b − τ c )) was the most appropriate choice to describe the field observations when the critical shear stress τ c and the proportionality factor M were kept constant.

California↗

USGS2014

Lateral baroclinic forcing enhances sediment transport from shallows to channel in an estuary

We investigate the dynamics governing exchange of sediment between estuarine shallows and the channel based on field measurements at eight stations spanning the interface between the channel and the extensive eastern shoals of South San Francisco Bay. The study site is characterized by longitudinally homogeneous bathymetry and a straight channel, with friction more important than the Coriolis forcing. Data were collected for 3 weeks in the winter and 4 weeks in the late summer of 2009, to capture a range of hydrologic and meteorologic conditions. The greatest sediment transport from shallows to channel occurred during a pair of strong, late-summer wind events, with westerly winds exceeding 10 m/s for more than 24 h. A combination of wind-driven barotropic return flow and lateral baroclinic circulation caused the transport. The lateral density gradient was produced by differences in temperature and suspended sediment concentration (SSC). During the wind events, SSC-induced vertical density stratification limited turbulent mixing at slack tides in the shallows, increasing the potential for two-layer exchange. The temperature- and SSC-induced lateral density gradient was comparable in strength to salinity-induced gradients in South Bay produced by seasonal freshwater inflows, but shorter in duration. In the absence of a lateral density gradient, suspended sediment flux at the channel slope was directed towards the shallows, both in winter and during summer sea breeze conditions, indicating the importance of baroclinically driven exchange to supply of sediment from the shallows to the channel in South San Francisco Bay and systems with similar bathymetry.

California↗

USGS2013

Implications for future survival of delta smelt from four climate change scenarios for the Sacramento–San Joaquin Delta, California

Changes in the position of the low salinity zone, a habitat suitability index, turbidity, and water temperature modeled from four 100-year scenarios of climate change were evaluated for possible effects on delta smelt Hypomesus transpacificus , which is endemic to the Sacramento–San Joaquin Delta. The persistence of delta smelt in much of its current habitat into the next century appears uncertain. By mid-century, the position of the low salinity zone in the fall and the habitat suitability index converged on values only observed during the worst droughts of the baseline period (1969–2000). Projected higher water temperatures would render waters historically inhabited by delta smelt near the confluence of the Sacramento and San Joaquin rivers largely uninhabitable. However, the scenarios of climate change are based on assumptions that require caution in the interpretation of the results. Projections like these provide managers with a useful tool for anticipating long-term challenges to managing fish populations and possibly adapting water management to ameliorate those challenges.

California↗

USGS2011

Projected evolution of California's San Francisco Bay-Delta-River System in a century of continuing climate change

Background Accumulating evidence shows that the planet is warming as a response to human emissions of greenhouse gases. Strategies of adaptation to climate change will require quantitative projections of how altered regional patterns of temperature, precipitation and sea level could cascade to provoke local impacts such as modified water supplies, increasing risks of coastal flooding, and growing challenges to sustainability of native species. Methodology/Principal Findings We linked a series of models to investigate responses of California's San Francisco Estuary-Watershed (SFEW) system to two contrasting scenarios of climate change. Model outputs for scenarios of fast and moderate warming are presented as 2010–2099 projections of nine indicators of changing climate, hydrology and habitat quality. Trends of these indicators measure rates of: increasing air and water temperatures, salinity and sea level; decreasing precipitation, runoff, snowmelt contribution to runoff, and suspended sediment concentrations; and increasing frequency of extreme environmental conditions such as water temperatures and sea level beyond the ranges of historical observations. Conclusions/Significance Most of these environmental indicators change substantially over the 21 st century, and many would present challenges to natural and managed systems. Adaptations to these changes will require flexible planning to cope with growing risks to humans and the challenges of meeting demands for fresh water and sustaining native biota. Programs of ecosystem rehabilitation and biodiversity conservation in coastal landscapes will be most likely to meet their objectives if they are designed from considerations that include: (1) an integrated perspective that river-estuary systems are influenced by effects of climate change operating on both watersheds and oceans; (2) varying sensitivity among environmental indicators to the uncertainty of future climates; (3) inevitability of biological community changes as responses to cumulative effects of climate change and other drivers of habitat transformations; and (4) anticipation and adaptation to the growing probability of ecosystem regime shifts.

California↗

USGS2011

Statistical models of temperature in the Sacramento-San Joaquin Delta under climate-change scenarios and ecological implications

Changes in water temperatures caused by climate change in California's Sacramento–San Joaquin Delta will affect the ecosystem through physiological rates of fishes and invertebrates. This study presents statistical models that can be used to forecast water temperature within the Delta as a response to atmospheric conditions. The daily average model performed well ( R 2 values greater than 0.93 during verification periods) for all stations within the Delta and San Francisco Bay provided there was at least 1 year of calibration data. To provide long-term projections of Delta water temperature, we forced the model with downscaled data from climate scenarios. Based on these projections, the ecological implications for the delta smelt, a key species, were assessed based on temperature thresholds. The model forecasts increases in the number of days above temperatures causing high mortality (especially along the Sacramento River) and a shift in thermal conditions for spawning to earlier in the year.

California↗

USGS2010

Wind-enhanced resuspension in the shallow waters of South San Francisco Bay: Mechanisms and potential implications for cohesive sediment transport

We investigated the driving forces of sediment dynamics at the shoals in South San Francisco Bay. Two stations were deployed along a line perpendicular to a 14 m deep channel, 1000 and 2000 m from the middle of the channel. Station depths were 2.59 and 2.19 m below mean lower low water, respectively. We used acoustic Doppler velocimeters for the simultaneous determination of current velocities, turbulence, sediment concentration and fluxes. Maximum current shear velocities were 0.015 m s −1 at the station further from the channel (closer to the shore) and 0.02 m s −1 at the station closer to the channel. Peak wave-induced shear velocities exceeded 0.015 m s −1 at both stations. Maximum sediment concentrations were around 30 g m −3 during calm periods (root mean square wave height <0.15 m). During wavy periods, sediment concentrations increased to 100 g m −3 and sediment fluxes were 5 times higher than in calm conditions (0.02 g m −2 s −1 versus >0.10 g m −2 s −1 ) at the station further from the channel 0.36 m above the bed. Closer to the channel, sediment concentrations and vertical fluxes due to wind wave resuspension were persistently lower (maximum concentrations around 50 g m −3 and maximum fluxes around 0.04 g m −2 s −1 ). Most resuspension events occurred during flood tides that followed wave events during low water. Although wave motions are able to resuspend sediment into the wave boundary layer at low tide, the observed large increases in sediment fluxes are due to the nonlinear interaction of wind waves and the tidal currents.

California↗