Ralston David K.

No Thumbnail Available
Last Name
Ralston
First Name
David K.
ORCID
0000-0002-0774-3101

Search Results

Now showing 1 - 20 of 52
  • Article
    Impacts of storm surge barriers on drag, mixing, and exchange flow in a partially mixed estuary
    (American Geophysical Union, 2022-04-11) Ralston, David K.
    Storm surge barriers are increasingly being considered as risk mitigation measures for coastal population centers. During non-storm periods, permanent barrier infrastructure reduces the flow cross-sectional area and affects tidal exchange. Effects of barrier structures on estuarine tidal and salinity dynamics have not been extensively examined, particularly for partially mixed estuaries. A nested, high-resolution model is used to characterize impacts of a potential storm surge barrier near the mouth of the Hudson River estuary. Maximum tidal velocities through barrier openings are more than double those in the base case. Landward of the barrier, tidal amplitude decreases on average by about 6% due to increased drag. The drag coefficient with the barrier is about 5 times greater than the base case due primarily to form drag from flow separation at barrier structures rather than increased bottom friction. The form drag scales with barrier geometry similar to previous studies of flow around headlands. Tidal water levels are reduced particularly during spring tides, such that marsh inundation frequency is reduced up to 25%. Strong tidal velocities through barrier openings enhance salinity mixing locally, but overall mixing in the estuary decreases due to reduced tidal velocities. Correspondingly, stratification decreases near the barrier and increases landward in the estuary. The salinity intrusion length increases by 5%–15% depending on discharge due to the decreased mixing and increased exchange flow. Exchange flow increases near the barrier due reflux into the lower layer with the increased mixing, which has the potential to increase estuarine residence times.
  • Article
    Sediment transport time scales and trapping efficiency in a tidal river
    (John Wiley & Sons, 2017-11-02) Ralston, David K. ; Geyer, W. Rockwell
    Observations and a numerical model are used to characterize sediment transport in the tidal Hudson River. A sediment budget over 11 years including major discharge events indicates the tidal fresh region traps about 40% of the sediment input from the watershed. Sediment input scales with the river discharge cubed, while seaward transport in the tidal river scales linearly, so the tidal river accumulates sediment during the highest discharge events. Sediment pulses associated with discharge events dissipate moving seaward and lag the advection speed of the river by a factor of 1.5 to 3. Idealized model simulations with a range of discharge and settling velocity were used to evaluate the trapping efficiency, transport rate, and mean age of sediment input from the watershed. The seaward transport of suspended sediment scales linearly with discharge but lags the river velocity by a factor that is linear with settling velocity. The lag factor is 30–40 times the settling velocity (mm s−1), so transport speeds vary by orders of magnitude from clay (0.01 mm s−1) to coarse silt (1 mm s−1). Deposition along the tidal river depends strongly on settling velocity, and a simple advection-reaction equation represents the loss due to settling on depositional shoals. The long-term discharge record is used to represent statistically the distribution of transport times, and time scales for settling velocities of 0.1 mm s−1 and 1 mm s−1 range from several months to several years for transport through the tidal river and several years to several decades through the estuary.
  • Article
    Wave generation, dissipation, and disequilibrium in an embayment with complex bathymetry
    (American Geophysical Union, 2018-10-04) Chen, Jia-Lin ; Ralston, David K. ; Geyer, W. Rockwell ; Sommerfield, Christopher K. ; Chant, Robert J.
    Heterogeneous, sharply varying bathymetry is common in estuaries and embayments, and complex interactions between the bathymetry and wave processes fundamentally alter the distribution of wave energy. The mechanisms that control the generation and dissipation of wind waves in an embayment with heterogeneous, sharply varying bathymetry are evaluated with an observational and numerical study of the Delaware Estuary. Waves in the lower bay depend on both local wind forcing and remote wave forcing from offshore, but elsewhere in the estuary waves are controlled by the local winds and the response of the wavefield to bathymetric variability. Differences in the wavefield with wind direction highlight the impacts of heterogeneous bathymetry and limited fetch. Under the typical winter northwest wind conditions waves are fetch‐limited in the middle estuary and reach equilibrium with local water depth only in the lower bay. During southerly wind conditions typical of storms, wave energy is near equilibrium in the lower bay, and midestuary waves are attenuated by the combination of whitecapping and bottom friction, particularly over the steep, longitudinal shoals. Although the energy dissipation due to bottom friction is generally small relative to whitecapping, it becomes significant where the waves shoal abruptly due to steep bottom topography. In contrast, directional spreading keeps wave heights in the main channel significantly less than local equilibrium. The wave disequilibrium in the deep navigational channel explains why the marked increase in depth by dredging of the modern channel has had little impact on wave conditions.
  • Article
    Quantifying vertical mixing in estuaries
    (Springer, 2008-10-21) Geyer, W. Rockwell ; Scully, Malcolm E. ; Ralston, David K.
    Estuarine turbulence is notable in that both the dissipation rate and the buoyancy frequency extend to much higher values than in other natural environments. The high dissipation rates lead to a distinct inertial subrange in the velocity and scalar spectra, which can be exploited for quantifying the turbulence quantities. However, high buoyancy frequencies lead to small Ozmidov scales, which require high sampling rates and small spatial aperture to resolve the turbulent fluxes. A set of observations in a highly stratified estuary demonstrate the effectiveness of a vessel-mounted turbulence array for resolving turbulent processes, and for relating the turbulence to the forcing by the Reynolds-averaged flow. The observations focus on the ebb, when most of the buoyancy flux occurs. Three stages of mixing are observed: (1) intermittent and localized but intense shear instability during the early ebb; (2) continuous and relatively homogeneous shear-induced mixing during the mid-ebb, and weakly stratified, boundary-layer mixing during the late ebb. The mixing efficiency as quantified by the flux Richardson number Rf was frequently observed to be higher than the canonical value of 0.15 from Osborn (J Phys Oceanogr 10:83–89, 1980). The high efficiency may be linked to the temporal–spatial evolution of shear instabilities.
  • Article
    A crab swarm at an ecological hotspot : patchiness and population density from AUV observations at a coastal, tropical seamount
    (PeerJ, 2016-04-12) Pineda, Jesus ; Cho, Walter W. ; Starczak, Victoria R. ; Govindarajan, Annette F. ; Guzman, Hector M. ; Girdhar, Yogesh ; Holleman, Rusty C. ; Churchill, James H. ; Singh, Hanumant ; Ralston, David K.
    A research cruise to Hannibal Bank, a seamount and an ecological hotspot in the coastal eastern tropical Pacific Ocean off Panama, explored the zonation, biodiversity, and the ecological processes that contribute to the seamount’s elevated biomass. Here we describe the spatial structure of a benthic anomuran red crab population, using submarine video and autonomous underwater vehicle (AUV) photographs. High density aggregations and a swarm of red crabs were associated with a dense turbid layer 4–10 m above the bottom. The high density aggregations were constrained to 355–385 m water depth over the Northwest flank of the seamount, although the crabs also occurred at lower densities in shallower waters (∼280 m) and in another location of the seamount. The crab aggregations occurred in hypoxic water, with oxygen levels of 0.04 ml/l. Barcoding of Hannibal red crabs, and pelagic red crabs sampled in a mass stranding event in 2015 at a beach in San Diego, California, USA, revealed that the Panamanian and the Californian crabs are likely the same species, Pleuroncodes planipes, and these findings represent an extension of the southern endrange of this species. Measurements along a 1.6 km transect revealed three high density aggregations, with the highest density up to 78 crabs/m2, and that the crabs were patchily distributed. Crab density peaked in the middle of the patch, a density structure similar to that of swarming insects.
  • Article
    Estuarine exchange flow variability in a seasonal, segmented estuary
    (American Meteorological Society, 2020-02-26) Conroy, Ted ; Sutherland, David A. ; Ralston, David K.
    Small estuaries in Mediterranean climates display pronounced salinity variability at seasonal and event time scales. Here, we use a hydrodynamic model of the Coos Estuary, Oregon, to examine the seasonal variability of the salinity dynamics and estuarine exchange flow. The exchange flow is primarily driven by tidal processes, varying with the spring–neap cycle rather than discharge or the salinity gradient. The salinity distribution is rarely in equilibrium with discharge conditions because during the wet season the response time scale is longer than discharge events, while during low flow it is longer than the entire dry season. Consequently, the salt field is rarely fully adjusted to the forcing and common power-law relations between the salinity intrusion and discharge do not apply. Further complicating the salinity dynamics is the estuarine geometry that consists of multiple branching channel segments with distinct freshwater sources. These channel segments act as subestuaries that import both higher- and lower-salinity water and export intermediate salinities. Throughout the estuary, tidal dispersion scales with tidal velocity squared, and likely includes jet–sink flow at the mouth, lateral shear dispersion, and tidal trapping in branching channel segments inside the estuary. While the estuarine inflow is strongly correlated with tidal amplitude, the outflow, stratification, and total mixing in the estuary are dependent on the seasonal variation in river discharge, which is similar to estuaries that are dominated by subtidal exchange flow.
  • Article
    Rapid growth and concerted sexual transitions by a bloom of the harmful dinoflagellate Alexandrium fundyense (Dinophyceae)
    (John Wiley & Sons, 2015-09-18) Brosnahan, Michael L. ; Velo-Suarez, Lourdes ; Ralston, David K. ; Fox, Sophia E. ; Sehein, Taylor R. ; Shalapyonok, Alexi ; Sosik, Heidi M. ; Olson, Robert J. ; Anderson, Donald M.
    Transitions between life cycle stages by the harmful dinoflagellate Alexandrium fundyense are critical for the initiation and termination of its blooms. To quantify these transitions in a single population, an Imaging FlowCytobot (IFCB), was deployed in Salt Pond (Eastham, Massachusetts), a small, tidally flushed kettle pond that hosts near annual, localized A. fundyense blooms. Machine-based image classifiers differentiating A. fundyense life cycle stages were developed and results were compared to manually corrected IFCB samples, manual microscopy-based estimates of A. fundyense abundance, previously published data describing prevalence of the parasite Amoebophrya, and a continuous culture of A. fundyense infected with Amoebophrya. In Salt Pond, a development phase of sustained vegetative division lasted approximately 3 weeks and was followed by a rapid and near complete conversion to small, gamete cells. The gametic period (∼3 d) coincided with a spike in the frequency of fusing gametes (up to 5% of A. fundyense images) and was followed by a zygotic phase (∼4 d) during which cell sizes returned to their normal range but cell division and diel vertical migration ceased. Cell division during bloom development was strongly phased, enabling estimation of daily rates of division, which were more than twice those predicted from batch cultures grown at similar temperatures in replete medium. Data from the Salt Pond deployment provide the first continuous record of an A. fundyense population through its complete bloom cycle and demonstrate growth and sexual induction rates much higher than are typically observed in culture.
  • Article
    Response to channel deepening of the salinity intrusion, estuarine circulation, and stratification in an urbanized estuary
    (American Geophysical Union, 2019-06-07) Ralston, David K. ; Geyer, W. Rockwell
    Modifications for navigation since the late 1800s have increased channel depth (H) in the lower Hudson River estuary by 10–30%, and at the mouth the depth has more than doubled. Observations along the lower estuary show that both salinity and stratification have increased over the past century. Model results comparing predredging bathymetry from the 1860s with modern conditions indicate an increase in the salinity intrusion of about 30%, which is roughly consistent with the H5/3 scaling expected from theory for salt flux dominated by steady exchange. While modifications including a recent deepening project have been concentrated near the mouth, the changes increase salinity and threaten drinking water supplies more than 100 km landward. The deepening has not changed the responses to river discharge (Qr) of the salinity intrusion (~Qr−1/3) or mean stratification (Qr2/3). Surprisingly, the increase in salinity intrusion with channel deepening results in almost no change in the estuarine circulation. This contrasts sharply with local scaling based on local dynamics of an H2 dependence, but it is consistent with a steady state salt balance that allows scaling of the estuarine circulation based on external forcing factors and is independent of depth. In contrast, the observed and modeled increases in stratification are opposite of expectations from the steady state balance, which could be due to reduction in mixing with loss of shallow subtidal regions. Overall, the mean shift in estuarine parameter space due to channel deepening has been modest compared with the monthly‐to‐seasonal variability due to tides and river discharge.
  • Article
    Salt wedge dynamics lead to enhanced sediment trapping within side embayments in high-energy estuaries
    (John Wiley & Sons, 2017-03-17) Yellen, Brian ; Woodruff, Jonathan D. ; Ralston, David K. ; MacDonald, Daniel G. ; Jones, David S.
    Off-river coves and embayments provide accommodation space for sediment accumulation, particularly for sandy estuaries where high energy in the main channel prevents significant long-term storage of fine-grained material. Seasonal sediment inputs to Hamburg Cove in the Connecticut River estuary (USA) were monitored to understand the timing and mechanisms for sediment storage there. Unlike in freshwater tidal coves, sediment was primarily trapped here during periods of low discharge, when the salinity intrusion extended upriver to the cove entrance. During periods of low discharge and high sediment accumulation, deposited sediment displayed geochemical signatures consistent with a marine source. Numerical simulations reveal that low discharge conditions provide several important characteristics that maximize sediment trapping. First, these conditions allow the estuarine turbidity maximum (ETM) to be located in the vicinity of the cove entrance, which increases sediment concentrations during flood tide. Second, the saltier water in the main channel can enter the cove as a density current, enhancing near-bed velocities and resuspending sediment, providing an efficient delivery mechanism. Finally, higher salinity water accumulates in the deep basin of the cove, creating a stratified region that becomes decoupled from ebb currents, promoting retention of sediment in the cove. This process of estuarine-enhanced sediment accumulation in off-river coves will likely extend upriver during future sea level rise.
  • Article
    Pseudo-nitzschia bloom dynamics in the Gulf of Maine: 2012-2016
    (Elsevier, 2019-08-19) Clark, Suzanna ; Hubbard, Katherine A. ; Anderson, Donald M. ; McGillicuddy, Dennis J. ; Ralston, David K. ; Townsend, David W.
    The toxic diatom genus Pseudo-nitzschia is a growing presence in the Gulf of Maine (GOM), where regionally unprecedented levels of domoic acid (DA) in 2016 led to the first Amnesic Shellfish Poisoning closures in the region. However, factors driving GOM Pseudo-nitzschia dynamics, DA concentrations, and the 2016 event are unclear. Water samples were collected at the surface and at depth in offshore transects in summer 2012, 2014, and 2015, and fall 2016, and a weekly time series of surface water samples was collected in 2013. Temperature and salinity data were obtained from NERACOOS buoys and measurements during sample collection. Samples were processed for particulate DA (pDA), dissolved nutrients (nitrate, ammonium, silicic acid, and phosphate), and cellular abundance. Species composition was estimated via Automated Ribosomal Intergenic Spacer Analysis (ARISA), a semi-quantitative DNA finger-printing tool. Pseudo-nitzschia biogeography was consistent in the years 2012, 2014, and 2015, with greater Pseudo-nitzschia cell abundance and P. plurisecta dominance in low-salinity inshore samples, and lower Pseudo-nitzschia cell abundance and P. delicatissima and P. seriata dominance in high-salinity offshore samples. During the 2016 event, pDA concentrations were an order of magnitude higher than in previous years, and inshore-offshore contrasts in biogeography were weak, with P. australis present in every sample. Patterns in temporal and spatial variability confirm that pDA increases with the abundance and the cellular DA of Pseudo-nitzschia species, but was not correlated with any one environmental factor. The greater pDA in 2016 was caused by P. australis – the observation of which is unprecedented in the region – and may have been exacerbated by low residual silicic acid. The novel presence of P. australis may be due to local growth conditions, the introduction of a population with an anomalous water mass, or both factors. A definitive cause of the 2016 bloom remains unknown, and continued DA monitoring in the GOM is warranted.
  • Article
    Reversed lateral circulation in a sharp estuarine bend with weak stratification
    (American Meteorological Society, 2019-06-13) Kranenburg, Wouter M. ; Geyer, W. Rockwell ; Garcia, Adrian Mikhail P. ; Ralston, David K.
    Although the hydrodynamics of river meanders are well studied, the influence of curvature on flow in estuaries, with alternating tidal flow and varying water levels and salinity gradients, is less well understood. This paper describes a field study on curvature effects in a narrow salt-marsh creek with sharp bends. The key observations, obtained during times of negligible stratification, are 1) distinct differences between secondary flow during ebb and flood, with helical circulation as in rivers during ebb and a reversed circulation during flood, and 2) maximum (ebb and flood) streamwise velocities near the inside of the bend, unlike typical river bend flow. The streamwise velocity structure is explained by the lack of a distinct point bar and the relatively deep cross section in the estuary, which means that curvature-induced inward momentum redistribution is not overcome by outward redistribution by frictional and topographic effects. Through differential advection of the along-estuary salinity gradient, the laterally sheared streamwise velocity generates lateral salinity differences, with the saltiest water near the inside during flood. The resulting lateral baroclinic pressure gradient force enhances the standard helical circulation during ebb but counteracts it during flood. This first leads to a reversed secondary circulation during flood in the outer part of the cross section, which triggers a positive feedback mechanism by bringing slower-moving water from the outside inward along the surface. This leads to a reversal of the vertical shear in the streamwise flow, and therefore in the centrifugal force, which further enhances the reversed secondary circulation.
  • Article
    Lateral circulation and sediment transport driven by axial winds in an idealized, partially mixed estuary
    (American Geophysical Union, 2009-12-03) Chen, Shih-Nan ; Sanford, Lawrence P. ; Ralston, David K.
    A 3D hydrodynamic model (ROMS) is used to investigate lateral circulation in a partially mixed estuary driven by axial wind events and to explore the associated transport of sediments. The channel is straight with a triangular cross section. The model results suggest that driving mechanisms for lateral circulation during axial wind events are different between stratified and unstratified conditions. When the water column is largely unstratified, rotational effects do not drive significant lateral circulation. Instead, differential advection of the axial salinity gradient by wind-driven axial flow is responsible for regulating the lateral salinity gradients that in turn drive bottom-divergent/convergent lateral circulation during down/up-estuary winds. From the subtidal lateral salt balance, it is found that the development of lateral salinity gradient by wind-induced differential advection is largely counterbalanced by the advection of salt by lateral circulation itself. When the water column is stratified, the lateral flow and salinity structures below the halocline closely resemble those driven by boundary mixing, and rotational effects are important. Lateral sediment flux and the event-integrated sediment transport are from channel to shoals during down-estuary winds but reversed for up-estuary winds. Potential impacts of wind-generated waves on lateral sediment transport are evaluated with two cases representing event conditions typical of upper Chesapeake Bay. Accounting for wind wave effects results in an order of magnitude increase in lateral sediment fluxes because of greater bottom stresses and sediment resuspension.
  • Preprint
    Dynamics of Alexandrium fundyense blooms and shellfish toxicity in the Nauset Marsh System of Cape Cod (Massachusetts, USA)
    ( 2011-08) Crespo, Bibiana G. ; Keafer, Bruce A. ; Ralston, David K. ; Lind, Henry ; Farber, Dawson ; Anderson, Donald M.
    Paralytic Shellfish Poisoning (PSP) toxins are annually recurrent along the Massachusetts coastline (USA), which includes many small embayments and salt ponds. Among these is the Nauset Marsh System (NMS), which has a long history of PSP toxicity. Little is known, however, about the bloom dynamics of the causative organism Alexandrium fundyense within that economically and socially important system. The overall goal of this work was to characterize the distribution and dynamics of A. fundyense blooms within the NMS and adjacent coastal waters by documenting the distribution and abundance of resting cysts and vegetative cells. Cysts were found predominantly in three drowned kettle holes or salt ponds at the distal ends of the NMS - Salt Pond, Mill Pond, and Town Cove. The central region of the NMS had a much lower concentration of cysts. Two types of A. fundyense blooms were observed. One originated entirely within the estuary, seeded by cysts in the three seedbeds. These blooms developed independently of each other and of the A. fundyense population observed in adjacent coastal waters outside the NMS. The temporal development of the blooms was different in the three salt ponds, with initiation differing by as much as 30 days. These differences do not appear to reflect the initial cyst abundances in these locations, and may simply result from higher cell retention and higher nutrient concentrations in Mill Pond, the first site to bloom. Germination of cysts accounted for a small percentage of the peak cell densities in the ponds, so population size was influenced more by the factors affecting growth than by cyst abundance. Subsurface cell aggregation (surface avoidance) limited advection of the vegetative A. fundyense cells out of the salt ponds through the shallow inlet channels. Thus, the upper reaches of the NMS are at the greatest risk for PSP since the highest cyst abundances and cell concentrations were found there. After these localized blooms in the salt ponds peaked and declined, a second, late season bloom occurred within the central portions of the NMS. The timing of this second bloom relative to those within the salt ponds and the coastal circulation patterns at that time strongly suggest that those cells originated from a regional A. fundyense bloom in the Gulf of Maine, delivered to the central marsh from coastal waters outside the NMS through Nauset Inlet. These results will guide policy decisions about water quality as well as shellfish monitoring and utilization within the NMS and highlight the potential for “surgical” closures of shellfish during PSP events, leaving some areas open for harvesting while others are closed.
  • Preprint
    Waves in the Red Sea : response to monsoonal and mountain gap winds
    ( 2013-05-30) Ralston, David K. ; Jiang, Houshuo ; Farrar, J. Thomas
    An unstructured grid, phase-averaged wave model forced with winds from a high resolution atmospheric model is used to evaluate wind wave conditions in the Red Sea over an approximately 2-year period. The Red Sea lies in a narrow rift valley, and the steep topography surrounding the basin steers the dominant wind patterns and consequently the wave climate. At large scales, the model results indicated that the primary seasonal variability in waves was due to the monsoonal wind reversal. During the winter, monsoon winds from the southeast generated waves with mean significant wave heights in excess of 2 m and mean periods of 8 s in the southern Red Sea, while in the northern part of the basin waves were smaller, shorter period, and from northwest. The zone of convergence of winds and waves typically occurred around 19-20˚N, but the location varied between 15 to 21.5˚N. During the summer, waves were generally smaller and from the northwest over most of the basin. While the seasonal winds oriented along the axis of the Red Sea drove much of the variability in the waves, the maximum wave heights in the simulations were not due to the monsoonal winds but instead were generated by localized mountain wind jets oriented across the basin (roughly east-west). During the summer, a mountain wind jet from the Tokar Gap enhanced the waves in the region of 18 and 20˚N, with monthly mean wave heights exceeding 2 m and maximum wave heights of 14 m during a period when the rest of the Red Sea was relatively calm. Smaller mountain gap wind jets along the northeast coast created large waves during the fall and winter, with a series of jets providing a dominant source of wave energy during these periods. Evaluation of the wave model results against observations from a buoy and satellites found that the spatial resolution of the wind model significantly affected the quality of the wave model results. Wind forcing from a 10-km grid produced higher skills for waves than winds from a 30-km grid, largely due to under-prediction of the mean wind speed and wave height with the coarser grid. The 30-km grid did not resolve the mountain gap wind jets, and thus predicted lower wave heights in the central Red Sea during the summer and along the northeast coast in the winter.
  • Article
    Progress and challenges in coupled hydrodynamic-ecological estuarine modeling
    (Springer, 2015-07-07) Ganju, Neil K. ; Brush, Mark J. ; Rashleigh, Brenda ; Aretxabaleta, Alfredo L. ; del Barrio, Pilar ; Grear, Jason S. ; Harris, Lora A. ; Lake, Samuel J. ; McCardell, Grant ; O’Donnell, James ; Ralston, David K. ; Signell, Richard P. ; Testa, Jeremy M. ; Vaudrey, Jamie M. P.
    Numerical modeling has emerged over the last several decades as a widely accepted tool for investigations in environmental sciences. In estuarine research, hydrodynamic and ecological models have moved along parallel tracks with regard to complexity, refinement, computational power, and incorporation of uncertainty. Coupled hydrodynamic-ecological models have been used to assess ecosystem processes and interactions, simulate future scenarios, and evaluate remedial actions in response to eutrophication, habitat loss, and freshwater diversion. The need to couple hydrodynamic and ecological models to address research and management questions is clear because dynamic feedbacks between biotic and physical processes are critical interactions within ecosystems. In this review, we present historical and modern perspectives on estuarine hydrodynamic and ecological modeling, consider model limitations, and address aspects of model linkage, skill assessment, and complexity. We discuss the balance between spatial and temporal resolution and present examples using different spatiotemporal scales. Finally, we recommend future lines of inquiry, approaches to balance complexity and uncertainty, and model transparency and utility. It is idealistic to think we can pursue a “theory of everything” for estuarine models, but recent advances suggest that models for both scientific investigations and management applications will continue to improve in terms of realism, precision, and accuracy.
  • Article
    The future of coastal and estuarine modeling: findings from a workshop
    (Elsevier, 2019-09-16) Fringer, Oliver B. ; Dawson, Clint N. ; He, Ruoying ; Ralston, David K. ; Zhang, Y. Joseph
    This paper summarizes the findings of a workshop convened in the United States in 2018 to discuss methods in coastal and estuarine modeling and to propose key areas of research and development needed to improve their accuracy and reliability. The focus of this paper is on physical processes, and we provide an overview of the current state-of-the-art based on presentations and discussions at the meeting, which revolved around the four primary themes of parameterizations, numerical methods, in-situ and remote-sensing measurements, and high-performance computing. A primary outcome of the workshop was agreement on the need to reduce subjectivity and improve reproducibility in modeling of physical processes in the coastal ocean. Reduction of subjectivity can be accomplished through development of standards for benchmarks, grid generation, and validation, and reproducibility can be improved through development of standards for input/output, coupling and model nesting, and reporting. Subjectivity can also be reduced through more engagement with the applied mathematics and computer science communities to develop methods for robust parameter estimation and uncertainty quantification. Such engagement could be encouraged through more collaboration between the forward and inverse modeling communities and integration of more applied math and computer science into oceanography curricula. Another outcome of the workshop was agreement on the need to develop high-resolution models that scale on advanced HPC systems to resolve, rather than parameterize, processes with horizontal scales that range between the depth and the internal Rossby deformation scale. Unsurprisingly, more research is needed on parameterizations of processes at scales smaller than the depth, including parameterizations for drag (including bottom roughness, bedforms, vegetation and corals), wave breaking, and air–sea interactions under strong wind conditions. Other topics that require significantly more work to better parameterize include nearshore wave modeling, sediment transport modeling, and morphodynamics. Finally, it was agreed that coastal models should be considered as key infrastructure needed to support research, just like laboratory facilities, field instrumentation, and research vessels. This will require a shift in the way proposals related to coastal ocean modeling are reviewed and funded.
  • Article
    Effects of locally generated wind waves on the momentum budget and subtidal exchange in a coastal plain estuary
    (American Geophysical Union, 2019-01-28) Pareja‐Roman, L. Fernando ; Chant, Robert J. ; Ralston, David K.
    A numerical model with a vortex force formalism is used to study the role of wind waves in the momentum budget and subtidal exchange of a shallow coastal plain estuary, Delaware Bay. Wave height and age in the bay have a spatial distribution that is controlled by bathymetry and fetch, with implications for the surface drag coefficient in young, underdeveloped seas. Inclusion of waves in the model leads to increases in the surface drag coefficient by up to 30% with respect to parameterizations in which surface drag is only a function of wind speed, in agreement with recent observations of air‐sea fluxes in estuaries. The model was modified to prevent whitecapping wave dissipation from generating breaking forces since that contribution is integrally equivalent to the wind stress. The proposed adjustment is consistent with previous studies of wave‐induced nearshore currents and with additional parameterizations for breaking forces in the model. The mean momentum balance during a simulated wind event was mainly between the pressure gradient force and surface stress, with negligible contributions by vortex, wave breaking (i.e., depth‐induced), and Stokes‐Coriolis forces. Modeled scenarios with realistic Delaware bathymetry suggest that the subtidal bay‐ocean exchange at storm time scales is sensitive to wave‐induced surface drag coefficient, wind direction, and mass transport due to the Stokes drift. Results herein are applicable to shallow coastal systems where the typical wave field is young (i.e., wind seas) and modulated by bathymetry.
  • Article
    Sediment transport and deposition on a river-dominated tidal flat : an idealized model study
    (American Geophysical Union, 2010-10-16) Chen, Shih-Nan ; Geyer, W. Rockwell ; Sherwood, Christopher R. ; Ralston, David K.
    A 3-D hydrodynamic model is used to investigate how different size classes of river-derived sediment are transported, exported and trapped on an idealized, river-dominated tidal flat. The model is composed of a river channel flanked by sloping tidal flats, a configuration motivated by the intertidal region of the Skagit River mouth in Washington State, United States. It is forced by mixed tides and a pulse of freshwater and sediment with various settling velocities. In this system, the river not only influences stratification but also contributes a significant cross-shore transport. As a result, the bottom stress is strongly ebb-dominated in the channel because of the seaward advance of strong river flow as the tidal flats drain during ebbs. Sediment deposition patterns and mass budgets are sensitive to settling velocity. The lateral sediment spreading scales with an advective distance (settling time multiplied by lateral flow speed), thereby confining the fast settling sediment classes in the channel. Residual sediment transport is landward on the flats, because of settling lag, but is strongly seaward in the channel. The seaward transport mainly occurs during big ebbs and is controlled by a length scale ratio Ld/XWL, where Ld is a cross-shore advective distance (settling time multiplied by river outlet velocity), and XWL is the immersed cross-shore length of the intertidal zone. Sediment trapping requires Ld/XWL < 1, leading to more trapping for the faster settling classes. Sensitivity studies show that including stratification and reducing tidal range both favor sediment trapping, whereas varying channel geometries and asymmetry of tides has relatively small impacts. Implications of the modeling results on the south Skagit intertidal region are discussed.
  • Article
    Flow separation and increased drag coefficient in estuarine channels with curvature
    (American Geophysical Union, 2020-10-14) Bo, Tong ; Ralston, David K.
    Flow separation has been observed and studied in sinuous laboratory channels and natural meanders, but the effects of flow separation on along‐channel drag are not well understood. Motivated by observations of large drag coefficients from a shallow, sinuous estuary, we built idealized numerical models representative of that system. We found that flow separation in tidal channels with curvature can create form drag that increases the total drag to more than twice that from bottom friction alone. In the momentum budget, the pressure gradient is balanced by the combined effects of bottom friction and form drag, which is calculated directly. The effective increase in total drag coefficient depends on two geometric parameters: dimensionless water depth and bend sharpness, quantified as the bend radius of curvature to channel width ratio. We introduce a theoretical boundary layer separation model to explain this parameter dependence and to predict flow separation and the increased drag. The drag coefficient can increase by a factor of 2–7 in “sharp” and “deep” sinuous channels where flow separation is most likely. Flow separation also enhances energy dissipation due to increased velocities in bends, resulting in greater loss of tidal energy and weakened stratification. Flow separation and the associated drag increase are expected to be more common in meanders of tidal channels than rivers where point bars that inhibit flow separation are more commonly found. The increased drag due to flow separation reduces tidal amplitude and affects velocity phasing along the estuary and could result in morphological feedbacks.
  • Article
    Turbulent mixing in a strongly forced salt wedge estuary
    (American Geophysical Union, 2010-12-09) Ralston, David K. ; Geyer, W. Rockwell ; Lerczak, James A. ; Scully, Malcolm E.
    Turbulent mixing of salt is examined in a shallow salt wedge estuary with strong fluvial and tidal forcing. A numerical model of the Merrimack River estuary is used to quantify turbulent stress, shear production, and buoyancy flux. Little mixing occurs during flood tides despite strong velocities because bottom boundary layer turbulence is dislocated from stratification elevated in the water column. During ebbs, bottom salinity fronts form at a series of bathymetric transitions. At the fronts, near-bottom velocity and shear stress are low, but shear, stress, and buoyancy flux are elevated at the pycnocline. Internal shear layers provide the dominant source of mixing during the early ebb. Later in the ebb, the pycnocline broadens and moves down such that boundary layer turbulence dominates mixing. Mixing occurs primarily during ebbs, with internal shear mixing accounting for about 50% of the total buoyancy flux. Both the relative contribution of internal shear mixing and the mixing efficiency increase with discharge, with bulk mixing efficiencies between 0.02 and 0.07. Buoyancy fluxes in the estuary increase with discharge up to about 400 m3 s−1 above which a majority of the mixing occurs offshore. Observed buoyancy fluxes were more consistent with the k-ɛ turbulence closure than the Mellor-Yamada closure, and more total mixing occurred in the estuary with k-ɛ. Calculated buoyancy fluxes were sensitive to horizontal grid resolution, as a lower resolution grid yielded less integrated buoyancy flux in the estuary and exported lower salinity water but likely had greater numerical mixing.