Lawrence P. Sanford

A unified erosion formulation for fine sediments

Abstract Different models of fine sediment transport often employ very different, frequently incompatible formulations for surface erosion of bottom sediment. In this paper, we develop a simple extension of the standard linear erosion formulation that allows it to be used to describe either Type I (depth-limited) erosion or Type II (unlimited) erosion, with a seamless transition between the two behaviors. The formulation is cast in terms of either the depth of erosion or eroded sediment mass. Assuming a locally constant rate of increase in critical stress with depth and direct proportionality between the erosion constant and sediment concentration at the interface, the model predicts the exponentially decaying erosion rate often observed in Type I erosion tests after application of each new shear stress step. The predicted decay rate is proportional to the rate of erosion per unit excess stress times the rate of increase in critical stress with depth. The formulation is tested by re-analyzing the data set presented by Maa et al. (1998) describing in situ erosion tests in Baltimore Harbor, MD, with generally favorable results. Solutions of the erosion formulation with time varying forcing show that erosion behavior is controlled by the ratio of the rate of change of shear stress to the rate of depletion of erodible sediment. If the time scale of shear stress change is long compare to the time scale of sediment depletion, then erosion rate is controlled by the time rate of increase in shear stress balanced against the depth rate of increase in critical stress (Type I behavior). If the time scale of shear stress change is short compared to the time scale of sediment depletion, then erosion rate is controlled by the instantaneous difference between bottom shear stress and critical shear stress (Type II behavior). A general algorithm for implementation in numerical sediment transport models is presented.

Reconsidering the physics of the Chesapeake Bay estuarine turbidity maximum

A series of cruises was carried out in the estuarine turbidity maximum (ETM) region of Chesapeake Bay in 1996 to examine physical and biological variability and dynamics. A large flood event in late January shifted the salinity structure of the upper Bay towards that of a salt wedge, but most of the massive sediment load delivered by the Susquehanna River appeared to bypass the ETM zone. In contrast, suspended sediments delivered during a flood event in late October were trapped very efficiently in the ETM. The difference in sediment trapping appeared to be due to increases in particle settling speed from January to October, suggesting that the fate of sediments delivered during large events may depend on the season in which they occur. The ETM roughly tracked the limit of salt (defined as the intersection of the 1 psu isohaline with the bottom) throughout the year, but it was often separated significantly from the limit of salt with the direction of separation unrelated to the phase of the tide. This was due to a lag of ETM sediment resuspension and transport behind rapid meteorologically induced or river flow induced motion of the salt limit. Examination of detailed time series of salt, suspended sediment, and velocity collected near the limit of salt, combined with other indications, led to the conclusion that the convergence of the estuarine circulation at the limit of salt is not the primary mechanism of particle trapping in the Chesapeake Bay ETM. This convergence and its associated salinity structure contribute to strong tidal asymmetries in sediment resuspension and transport that collect and maintain a resuspendable pool of rapidly settling particles near the salt limit. Without tidal resuspension and transport, the ETM would either not exist or be greatly weakened. In spite of this repeated resuspension, sedimentation is the ultimate fate of most terrigenous material delivered to the Chesapeake Bay ETM. Sedimentation rates in the ETM channel are at least an order of magnitude greater than on the adjacent shoals, probably due to focusing mechanisms that are poorly understood.

Estimation of net physical transport and hydraulic residence times for a coastal plain estuary using box models

A box model based on salinity distributions and freshwater inflow measurements was developed and used to estimate net non-tidal physical circulation and hydraulic residence times for Patuxent River estuary, Maryland, a tributary estuary of Chesapeake Bay. The box model relaxes the usual assumption that salinity is at steady-state, an important improvement over previous box model studies, yet it remains simple enough to have broad appeal. Average monthly 2-dimensional net non-tidal circulation and residence times for 1986–1995 are estimated and related to river flow and salt water inflow as estimated by the box model. An important result is that advective exchange at the estuary mouth was not correlated with Patuxent River flow, most likely due to effects of offshore salinity changes in Chesapeake Bay. The median residence time for freshwater entering at the head of the estuary was 68 d and decreased hyperbolically with increasing river flow to 30 d during high flow. Estimates of residence times for down-estuary points of origin showed that, from the head of the estuary to its mouth, control of flushing changed from primarily river flow to other factors regulating the intensity of gravitational circulation.

Wave-forced resuspension of upper Chesapeake Bay muds

Moored instruments were used to make observations of near bottom currents, waves, temperature, salinity, and turbidity at shallow (3.5 m and 5.5 m depth) dredged sediment disposal sites in upper Chesapeake Bay during the winters of 1990 and 1991 to investigate time-varying characteristics of resuspension processes over extended periods. Resulting time series data show the variability of two components of the suspended sediment concentration field. Background suspended sediment concentrations varied inversely with salinity and in direct relation to Susquehanna River flow. Muddy bottom sediments were also resuspended locally by both tidal currents and wind-wave forcing, resulting in short-term increases and decreases in suspended concentration, with higher peak concentrations near the bottom. In both years, episodes of wave-forced resuspension dominated tidal resuspension on an individual event basis, exceeding most tidal resuspension peaks by a factor of 3 to 5. The winds that generated the waves responsible for the observed resuspension events were not optimal for wave generation, however. Application of a simple wind-wave model showed that much greater wave-forced resuspension than that observed might be generated under the proper conditions. The consolidated sediments investigated in 1990 were less susceptible to both tidal and wave-forced resuspension than the recently deposited sediments investigated in 1991. There was also some indication that wave-forced resuspension increased erodibility of the bottom sediments on a short-term basis. Wave-forced resuspension is implicated as an important part of sediment transport processes in much of Chesapeake Bay. Its role in deeper, narrower, and more tidally energetic estuaries is not as clear, and should be investigated on a case-by-case basis.

Axial Wind Effects on Stratification and Longitudinal Salt Transport in an Idealized, Partially Mixed Estuary*

Abstract A 3D hydrodynamic model [Regional Ocean Model System (ROMS)] is used to investigate how axial wind influences stratification and to explore the associated longitudinal salt transport in partially mixed estuaries. The model is configured to represent a straight estuarine channel connecting to a shelf sea. The results confirm that wind straining of the along-channel salinity gradient exerts an important control on stratification. Two governing parameters are identified: the Wedderburn number (W) defined as the ratio of wind stress to axial baroclinic pressure gradient force, and the ratio of an entrainment depth to water depth (hs/H). Here W controls the effectiveness of wind straining, which promotes increases (decreases) in stratification during down-estuary (up-estuary) wind. The ratio hs/H determines the portion of the water column affected by direct wind mixing. While stratification is always reduced by up-estuary wind, stratification shows an increase-then-decrease transition when down-estuar...

A Nearshore Model to Investigate the Effects of Seagrass Bed Geometry on Wave Attenuation and Suspended Sediment Transport

The effects of seagrass bed geometry on wave attenuation and suspended sediment transport were investigated using a modified Nearshore Community Model (NearCoM). The model was enhanced to account for cohesive sediment erosion and deposition, sediment transport, combined wave and current shear stresses, and seagrass effects on drag. Expressions for seagrass drag as a function of seagrass shoot density and canopy height were derived from published flume studies of model vegetation. The predicted reduction of volume flux for steady flow through a bed agreed reasonably well with a separate flume study. Predicted wave attenuation qualitatively captured seasonal patterns observed in the field: wave attenuation peaked during the flowering season and decreased as shoot density and canopy height decreased. Model scenarios with idealized bathymetries demonstrated that, when wave orbital velocities and the seagrass canopy interact, increasing seagrass bed width in the direction of wave propagation results in higher wave attenuation, and increasing incoming wave height results in higher relative wave attenuation. The model also predicted lower skin friction, reduced erosion rates, and higher bottom sediment accumulation within and behind the bed. Reduced erosion rates within seagrass beds have been reported, but reductions in stress behind the bed require further studies for verification. Model results suggest that the mechanism of sediment trapping by seagrass beds is more complex than reduced erosion rates alone; it also requires suspended sediment sources outside of the bed and horizontal transport into the bed.

The influence of wind and river pulses on an estuarine turbidity maximum: Numerical studies and field observations in Chesapeake Bay

The effect of pulsed events on estuarine turbidity maxima (ETM) was investigated with the Princeton Ocean Model, a three-dimensional hydrodynamic model. The theoretical model was adapted to a straight-channel estuary and enhanced with sediment transport, erosion, deposition, and burial components. Wind and river pulse scenarios from the numerical model were compared to field observations before and after river pulse and wind events in upper Chesapeake Bay. Numerical studies and field observations demonstrated that the salt front and ETM had rapid and nonlinear responses to short-term pulses in river flow and wind. Although increases and decreases in river flow caused down-estuary and up-estuary (respectively) movements of the salt front, the effect of increased river flow was more pronounced than that of decreased river flow. Along-channel wind events also elicited non-linear responses. The salt front moved in the opposite direction of wind stress, shifting up-estuary in response to down-estuary winds and vice-versa.Modeled pulsed events affected suspended sediment distributions by modifying the location of the salt front, near-bottom shear stress, and the location of bottom sediment in relation to stratification within the salt front. Bottom sediment accumulated near the convergent zone at the tip of the salt front, but lagged behind the rapid response of the salt front during wind events. While increases in river flow and along-channel winds resulted in sediment transport down-estuary, only reductions in river flow resulted in consistent up-estuary movement of bottom sediment. Model predictions suggest that wind and river pulse events significantly influence salt front structure and circulation patterns, and have an important role in the transport of sediment in upper estuaries.

The influence of episodic events on transport of striped bass eggs to the estuarine turbidity maximum nursery area

The estuarine turbidity maximum (ETM) is an important nursery area for anadromous fish where early-life stages can be retained in high prey concentrations and favorable salinities. Episodic freshwater flow and wind events could influence the transport of striped bass (Morone saxatilis) eggs to the ETM. This hypothesis was evaluated with regression analysis of observational data and with a coupled biological-physical model of a semi-idealized upper Chesapeake Bay driven by observed wind and freshwater flow. A particle-tracking model was constructed within a numerical circulation model (Princeton Ocean Model) to simulate the transport of fish eggs in a 3-dimensional flow field. Particles with the sinking speed of striped bass eggs were released up-estuary of the salt front in both 2-d event-scale and 60-d seasonal-scale scenarios. In event scenarios, egg-like particles with observed specific gravities (densities) of striped bass eggs were transported to the optimum ETM nursery area after 2 d, the striped bass egg-stage duration. Wind events and pulses in river discharge decreased the number of egg-like particles transported to the ETM area by 20.9% and 13.2%, respectively, compared to nonevent conditions. In seasonal scenarios, particle delivery to the ETM depended upon the timing of the release of egg-like particles. The number of particles transported to the ETM area decreased when particles were released before and during wind and river pulse events. Particle delivery to the ETM area was enhanced when the salt front was moving up-estuary after river pulse events and as base river flow receded over the spawning season. Model results suggest that the timing of striped bass spawning in relation to pulsed events may have a negative (before or during events) or positive (after river flow events) effect on egg transport. Spawning after river flow events may promote early-stage survival by taking advantage of improved transport, enhanced turbidity refuge, and elevated prey production that may occur after river pulse events. In multiple regression analysis of observed data, mean spring freshwater flow rates and the number of pulsed freshwater flow events during the striped bass spawning season explained 71% of the variability in striped bass juvenile abundance in upper Chesapeake Bay from 1986 to 2002. Positive parameter estimates for these effects support the hypothesis that pulsed freshwater flow events, coupled with spawning after the events, may enhance striped bass early-stage survival. Results suggest that episodic events may have an important role in controlling fish recruitment.

Drag Coefficients with Fetch-Limited Wind Waves*

Air‐sea fluxes of momentum and heat were measured simultaneously with surface wind waves and nearsurface currents in mid Chesapeake Bay during summer 1998 under low wind conditions. The data were collected using a Gill sonic anemometer and a Sontek Acoustic Doppler Velocimeter with a pressure sensor mounted on a temporary fixed tower in 8.8 m of water. The analyzed data show that the neutral drag coefficients depend upon both wind speed and wave age. They are better correlated to wave age than to wind speed. Data scatter is significantly larger in low winds than in high winds. Under light winds, the neutral drag coefficients increase with decreasing wind speed and have values much higher than those for relatively higher wind speeds. At higher wind speeds, neutral drag coefficients increase with increasing wind speed. Regardless of wind speed, neutral drag coefficients always decrease with increasing wave age. Neutral drag coefficients are lower than the results of similar field studies when fit to wind speed alone, but they statistically agree with other studies if they are fit to wave age. The momentum transfer mechanism is investigated using a parametric wave model with a seastate-dependent form drag and a reference system moving with the waves. The relationship between modeled drag coefficient and modeled wave age agrees well with the relationship derived from the data.

IN SITU CHARACTERIZATION OF ESTUARINE SUSPENDED SEDIMENT IN THE PRESENCE OF MUDDY FLOCS AND PELLETS

Observations are presented from a benthic observatory in the middle reaches of the York River estuary, VA, USA, that show evidence for both muddy flocs and pellets in the lower 1 m of the water column. This study combines in situ time series estimates of (i) volume concentration and particle size distribution from a Laser In Situ Scattering Transmisometer (LISST) (for 2.5-500 μm) and a high-definition particle camera (for 20 μm to 20 mm), and (ii) water velocity, turbulent stress, mass concentration and settling velocity derived from an Acoustic Doppler Velocimeter (ADV). Mass concentration, mass settling velocity and the abundant 88 μm size class are in phase with velocity and stress, consistent with suspension of relatively dense, rapidly settling and resilient ~90 μm pellets. Volume concentration of the abundant 280 mm class peaks well after stress and velocity begin to decrease, consistent with the formation of lower density, slowly settling and fragile ~300 μm flocs.