skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Data from: Dynamics and scaling of a small river discharging into the surf zone
We use an idealized numerical model to investigate the dynamics and fate of a small river discharging into the surf zone. Our study reveals that the plume reaches a steady state, at which point the combined advective and diffusive freshwater fluxes from the surf zone to the inner shelf balance the river discharge. At a steady state, the surf zone is well-mixed vertically due to wave-enhanced vertical turbulent diffusion and has a strong cross-shore salinity gradient. The horizontal gradient drives a cross-shore buoyancy-driven circulation, directed offshore at the surface and onshore near the bottom, which opposes the wave-driven circulation. Using a scaling analysis based on momentum and freshwater budgets, we determine that the steady-state alongshore plume extent (Lp) and the fraction of river water trapped in the surf zone depend on the ratio of the near-field plume length to the surf zone width (Lnf/Lsz) across a wide range of discharge and wave conditions, and a limited set of tidal conditions. This scaling also allows us to predict the residence time and freshwater fraction (or dilution ratio) in the steady-state plume within the surf zone, which range from approximately 0.1 to 10 days and 0.1 to 0.3, respectively. Our findings establish the basic dynamics and scales of an idealized plume in the surf zone, as well as estimates of residence times and dilution rates that may provide guidance to coastal managers. # Data from: Dynamics and scaling of a small river discharging into the surf zone [https://doi.org/10.5061/dryad.2280gb608](https://doi.org/10.5061/dryad.2280gb608) The present dataset includes the [COAWST model](https://www.usgs.gov/centers/whcmsc/science/coawst-a-coupled-ocean-atmosphere-wave-sediment-transport-modeling-system) outputs used to describe the dynamics and scaling of a small river discharging into the surf zone. ## File structure The data are structured as follows: 1. plume_scale.mat - Data of plume scales of all the cases, where * Hs: significant wave height [m] * Q: river discharge [m^3 s^-1] * L_nf: near-field plume length [m] * L_p: alongshore plume extent [m] * h_sz: water depth at the surf zone edge [m] * x_sz: surf zone width [m] * S_in: inflow salinity [PSU] * g_p: reduced gravity at the river mouth [m s^-2] * g_p*_*0: reduced gravity at the river mouth calculated using the density difference between river inflow and ambient ocean water [m s^-2] * Eta_0: water surface elevation anomaly at the river mouth [m] * V_sz: total volume of freshwater trapped in the surf zone [m^3] * T: the time required for the plume to reach a steady state [day] * L_t: plume turning distance [m] * S_bar: averaged salinity in the plume [PSU] 2. DepthAveraged.mat - Depth-averaged flow fields. DepthAveraged_BaseCase.mat, DepthAveraged_Case1.mat, DepthAveraged_Case3.mat, DepthAveraged_Case4.mat, DepthAveraged_Case6.mat, DepthAveraged_Case7.mat, DepthAveraged_Case8.mat, DepthAveraged_Case9.mat, DepthAveraged_Case16.mat, DepthAveraged_Case17.mat, DepthAveraged_Case18.mat, DepthAveraged_Case19.mat includes the results of the base case, cases 1, 3, 4, 6-9, and 16-19, respectively. In these files: * Wetdry_mask: wet/dry mask on RHO-points [binary] * Wetdry_mask_u: wet/dry mask on U-points [binary] * Wetdry_mask_v: wet/dry mask on V-points [binary] * Z: free-surface [m] * S: surface salinity [PSU] * Hs: significant wave height [m] * U: vertically integrated u-momentum component [m s^-1] * U_st: vertically-integrated u-Stokes drift velocity [m s^-1] * V: vertically integrated v-momentum component [m s^-1] * V_st: vertically-integrated v-Stokes drift velocity [m s^-1] 3. FullField_BaseCase.mat - 3D flow fields for the base case, where * Z: free-surface [m] * S: salinity [PSU] * Hs: significant wave height [m] * Lw: mean wavelength [m] * U: u-momentum component [m s^-1] * U_st: u-Stokes drift velocity [m s^-1] * V: v-momentum component [m s^-1] * V_st: v-Stokes drift velocity [m s^-1] * W: w-momentum component [m s^-1] * W_st: w-Stokes drift velocity [m s^-1] * Aks: salinity vertical diffusion coefficient [m^2 s^-1] * Akv: vertical viscosity coefficient [m^2 s^-1] * Cs_r: S-coordinate stretching curves at RHO-points [-] * Cs_w: S-coordinate stretching curves at W-points [-] 4. FreshwaterTrace_BaseCase.mat - Time series of freshwater volume and fluxes for the base case, where * i_sz: XI-index of the location of the surf zone edge [-] * i_shore: XI-index of the location of the shoreline [-] * Vsz: volume of freshwater in the plume in the surf zone [m^3] * Vis: volume of freshwater in the plume in the inner shelf [m^3] * Vsz_total: total volume of freshwater in the surf zone [m^3] * Vis_total: total volume of freshwater in the inner shelf [m^3] * R2SZ_flux: freshwater flux discharging into the surf zone [m^3 s^-1] * Vchannel: volume of freshwater in the plume in the river channel [m^3] * Vchannel_total: volume of freshwater in the river channel [m^3] * SBoundary_flux_SZ: the freshwater fluxes through the southern domain boundaries of the surf zone [m^3 s^-1] * SBoundary_flux_IS: the freshwater fluxes through the southern domain boundaries of the inner shelf [m^3 s^-1] * NBoundary_flux_SZ: the freshwater fluxes through the northern domain boundaries of the surf zone [m^3 s^-1] * NBoundary_flux_IS: the freshwater fluxes through the northern domain boundaries of the inner shelf [m^3 s^-1] * WBoundary_flux: the freshwater fluxes through the westhern domain boundary [m^3 s^-1] 5. DepthAveraged_XDiagnostic.mat - Depth-averaged diagnostic output of cross-shore momentum terms. DepthAveraged_XDiagnostic_BaseCase.mat includes the results of the base case at the steady state, and DepthAveraged_XDiagnostic_0day_1mWave.mat includes those at the start of river flow. In these files: * ubar_xadv: time-averaged 2D u-momentum, horizontal XI-advection term [m s^-2] * ubar_yadv: time-averaged 2D u-momentum, horizontal ETA-advection term [m s^-2] * ubar_xvisc: time-averaged 2D u-momentum, horizontal XI-viscosity term [m s^-2] * ubar_yvisc: time-averaged 2D u-momentum, horizontal ETA-viscosity term [m s^-2] * ubar_prsgrd: time-averaged 2D u-momentum, pressure gradient term [m s^-2] * ubar_zqsp: time-averaged 2D u-momentum, quasi-static pressure [m s^-2] * ubar_zbeh: time-averaged 2D u-momentum, Bernoulli head [m s^-2] * ubar_bstr: time-averaged 2D u-momentum, bottom stress term [m s^-2] * ubar_wbrk: time-averaged 2D u-momentum, wave breaking term [m s^-2] 6. DepthAveraged_YDiagnostic_BaseCase.mat - Depth-averaged diagnostic output of alongshore momentum terms, where * vbar_xadv: time-averaged 2D v-momentum, horizontal XI-advection term [m s^-2] * vbar_yadv: time-averaged 2D v-momentum, horizontal ETA-advection term [m s^-2] * vbar_xvisc: time-averaged 2D v-momentum, horizontal XI-viscosity term [m s^-2] * vbar_yvisc: time-averaged 2D v-momentum, horizontal ETA-viscosity term [m s^-2] * vbar_prsgrd: time-averaged 2D v-momentum, pressure gradient term [m s^-2] * vbar_zqsp: time-averaged 2D v-momentum, quasi-static pressure [m s^-2] * vbar_zbeh: time-averaged 2D v-momentum, Bernoulli head [m s^-2] * vbar_bstr: time-averaged 2D v-momentum, bottom stress term [m s^-2] * vbar_wbrk: time-averaged 2D v-momentum, wave breaking term [m s^-2] 7. grid.zip - Model grid file. * This grid file is designed for use with [ROMS](https://www.myroms.org/index.php), the hydrodynamic module of the COAWST modeling system. A diagram illustrating how the variables are placed on the grid and where the boundaries lie relative to the grid is available on [WikiROMS](https://www.myroms.org/wiki/Grid_Generation). * This grid file is in NetCDF format, which can be opened and used by a wide range of application software such as MATLAB, Python, and Panoply. For more detailed information, please refer to its [official website](https://www.unidata.ucar.edu/software/netcdf/). ## Code/Software All the post-processing scripts and data are prepared by MATLAB.  more » « less
Award ID(s):
1924005
PAR ID:
10656972
Author(s) / Creator(s):
; ; ; ; ; ;
Publisher / Repository:
Dryad
Date Published:
Edition / Version:
9
Subject(s) / Keyword(s):
FOS: Environmental engineering FOS: Environmental engineering Estuaries Rivers Coastal flows mixing wave breaking
Format(s):
Medium: X Size: 523424034 bytes
Size(s):
523424034 bytes
Sponsoring Org:
National Science Foundation
More Like this
  1. We use salinity observations from drifters and moorings at the Quinault River mouth to investigate mixing and stratification in a surf-zone-trapped river plume. We quantify mixing based on the rate of change of salinity DS/Dt in the drifters’ quasi-Lagrangian reference frame. We estimate a constant value of the vertical eddy diffusivity of salt of Kz=(2.2 +/- 0.6) x 10^-3 m^2 s^-1, based on the relationship between vertically integrated DS/Dt and stratification, with values as high as 1 x 10^-2 m^2 s^-1 when stratification is low. Mixing, quantified as DS/Dt, is directly correlated to surf-zone stratification, and is therefore modulated by changes in stratification caused by tidal variability in freshwater volume flux. High DS/Dt is observed when the near-surface stratification is high and salinity gradients are collocated with wave-breaking turbulence. We observe a transition from low stratification and low DS/Dt at low tidal stage to high stratification and high DS/Dt at high tidal stage. Observed wave-breaking turbulence does not change significantly with stratification, tidal stage, or offshore wave height; as a result, we observe no relationship between plume mixing and offshore wave height for the range of conditions sampled. Thus, plume mixing in the surf zone is altered by changes in stratification; these are due to tidal variability in freshwater flux from the river and not wave conditions, presumably because depth-limited wave breaking causes sufficient turbulence for mixing to occur during all observed conditions. 
    more » « less
  2. Abstract We use an idealized numerical model to investigate the dynamics and fate of a small river discharging into the surf zone. Our study reveals that the plume reaches a steady state, at which point the combined advective and diffusive freshwater fluxes from the surf zone to the inner shelf balance the river discharge. At a steady state, the surf zone is well mixed vertically due to wave-enhanced vertical turbulent diffusion and has a strong cross-shore salinity gradient. The horizontal gradient drives a cross-shore buoyancy-driven circulation, directed offshore at the surface and onshore near the bottom, which opposes the wave-driven circulation. Using a scaling analysis based on momentum and freshwater budgets, we determine that the steady-state alongshore plume extent (Lp) and the fraction of river water trapped in the surf zone depend on the ratio of the near-field plume length to the surf-zone width (Lnf/Lsz) across a wide range of discharge and wave conditions and a limited set of tidal conditions. This scaling also allows us to predict the residence time and freshwater fraction (or dilution ratio) in the steady-state plume within the surf zone, which ranges from approximately 0.1 to 10 days and from 0.1 to 0.3, respectively. Our findings establish the basic dynamics and scales of an idealized plume in the surf zone, as well as estimates of residence times and dilution rates that may provide guidance to coastal managers. Significance StatementSmall rivers and estuaries often carry pollutants, sediments, and larvae into the coastal ocean, where wave action in the surf zone can trap them near the shore. This process can play an important role in the flux of material into and out of the nearshore ecosystem and presents a potential risk to swimmers when materials are harmful. The present study uses a numerical model to investigate the fate of freshwater discharged from small rivers into the surf zone and the processes through which trapped riverine freshwater escapes from the surf zone. These results establish a basis for predicting the fate of river-borne materials from coastal rivers and understanding the exchange between the surf zone and the inner shelf. Additionally, this work provides a theoretical framework for predicting the residence time and concentration of river-borne material trapped in the surf zone. 
    more » « less
  3. null (Ed.)
    Abstract A unique feature of small mountainous rivers is that discharge can be elevated by an order of magnitude during a large rain event. The impact of time-varying discharge on freshwater transport pathways and alongshore propagation rates in the coastal ocean is not well understood. A suite of simulations in an idealized coastal ocean domain using the Regional Ocean Modeling System (ROMS) with varying steady background discharge conditions (25–100 m 3 s −1 ), pulse amplitude (200–800 m 3 s −1 ), pulse duration (1–6 days), and steady downwelling-favorable winds (0–4 m s −1 ) are compared to investigate the downstream freshwater transport along the coast (in the direction of Kelvin wave propagation) following a discharge pulse from the river. The nose of the pulse propagates rapidly alongshore at 0.04–0.32 m s −1 (faster propagation corresponds with larger pulse volume and faster winds) transporting 13%–66% of the discharge. The remainder of the discharge volume initially accumulates in the bulge near the river mouth, with lower retention for longer pulse duration and stronger winds. Following the pulse, the bulge eddy disconnects from the river mouth and is advected downstream at 0–0.1 m s −1 , equal to the depth-averaged wind-driven ambient water velocity. As it transits alongshore, it sheds freshwater volume farther downstream and the alongshore freshwater transport stays elevated between the nose and the transient bulge eddy. The evolution of freshwater transport at a plume cross section can be described by the background discharge, the passage of the pulse nose, and a slow exponential return to background conditions. 
    more » « less
  4. This archive contains COAWST model input, grids and initial conditions, and output used to produce the results in a submitted manuscript. The files are:</p> model_input.zip: input files for simulations presented in this paper   ocean_rip_current.in: ROMS ocean model input file   swan_rip_current.in: SWAN wave model input file (example with Hs=1m)   coupling_rip_current.in: model coupling file   rip_current.h: model header file    model_grids_forcing.zip: bathymetry and initial condition files      hbeach_grid_isbathy_2m.nc: ROMS bathymetry input file      hbeach_grid_isbathy_2m.bot: SWAN bathymetry input file      hbeach_grid_isbathy_2m.grd: SWAN grid input file      hbeach_init_isbathy_14_18_17.nc: Initial temperature, cool surf zone dT=-1C case      hbeach_init_isbathy_14_18_19.nc: Initial temperature, warm surf zone dT=+1C case      hbeach_init_isbathy_14_18_16.nc: Initial temperature, cool surf zone dT=-2C case      hbeach_init_isbathy_14_18_20.nc: Initial temperature, warm surf zone dT=+2C case      hbeach_init_isbathy_14_18_17p5.nc: Initial temperature, cool surf zone dT=-0.5C case      hbeach_init_isbathy_14_18_18p5.nc: Initial temperature, warm surf zone dT=+0.5C case</p> model_output files: model output used to produce the figures      netcdf files, zipped      variables included:           x_rho (cross-shore coordinate, m)           y_rho (alongshore coordinate, m)           z_rho (vertical coordinate, m)           ocean_time (time since initialization, s, output every 5 mins)           h (bathymetry, m)           temp (temperature, Celsius)           dye_02 (surfzone-released dye)           Hwave (wave height, m)           Dissip_break (wave dissipation W/m2)            ubar (cross-shore depth-average velocity, m/s, interpolated to rho-points)      Case_141817.nc: cool surf zone dT=-1C Hs=1m      Case_141819.nc: warm surf zone dT=+1C Hs=1m      Case_141816.nc: cool surf zone dT=-2C Hs=1m      Case_141820.nc: warm surf zone dT=-2C Hs=1m      Case_141817p5.nc: cool surf zone dT=-0.5C Hs=1m      Case_141818p5.nc: warm surf zone dT=+0.5C Hs=1m      Case_141817_Hp5.nc: cool surf zone dT=-1C Hs=0.5m      Case_141819_Hp5.nc: warm surf zone dT=+1C Hs=0.5m      Case_141817_Hp75.nc: cool surf zone dT=-1C Hs=0.75m      Case_141819_Hp75.nc: warm surf zone dT=+1C Hs=0.75m</p> COAWST is an open source code and can be download at https://coawstmodel-trac.sourcerepo.com/coawstmodel_COAWST/. Descriptions of the input and output files can be found in the manual distributed with the model code and in the glossary at the end of the ocean.in file.</p> Corresponding author: Melissa Moulton, mmoulton@uw.edu</p> 
    more » « less
  5. A gr e at d e al of i nt er e st s urr o u n d s t h e u s e of tr a n s cr a ni al dir e ct c urr e nt sti m ul ati o n (t D C S) t o a u g m e nt c o g niti v e tr ai ni n g. H o w e v er, eff e ct s ar e i n c o n si st e nt a cr o s s st u di e s, a n d m et aa n al yti c e vi d e n c e i s mi x e d, e s p e ci all y f o r h e alt h y, y o u n g a d ult s. O n e m aj or s o ur c e of t hi s i n c o n si st e n c y i s i n di vi d u al diff er e n c e s a m o n g t h e p arti ci p a nt s, b ut t h e s e diff er e n c e s ar e r ar el y e x a mi n e d i n t h e c o nt e xt of c o m bi n e d tr ai ni n g/ sti m ul ati o n st u di e s. I n a d diti o n, it i s u n cl e ar h o w l o n g t h e eff e ct s of sti m ul ati o n l a st, e v e n i n s u c c e s sf ul i nt er v e nti o n s. S o m e st u di e s m a k e u s e of f oll o w- u p a s s e s s m e nt s, b ut v er y f e w h a v e m e a s ur e d p erf or m a n c e m or e t h a n a f e w m o nt hs aft er a n i nt er v e nti o n. H er e, w e utili z e d d at a fr o m a pr e vi o u s st u d y of t D C S a n d c o g niti v e tr ai ni n g [ A u, J., K at z, B., B u s c h k u e hl, M., B u n arj o, K., S e n g er, T., Z a b el, C., et al. E n h a n ci n g w or ki n g m e m or y tr ai ni n g wit h tr a n scr a ni al dir e ct c urr e nt sti m ul ati o n. J o u r n al of C o g niti v e N e u r os ci e n c e, 2 8, 1 4 1 9 – 1 4 3 2, 2 0 1 6] i n w hi c h p arti ci p a nts tr ai n e d o n a w or ki n g m e m or y t as k o v er 7 d a y s w hil e r e c ei vi n g a cti v e or s h a m t D C S. A n e w, l o n g er-t er m f oll o w- u p t o a ss es s l at er p erf or m a n c e w a s c o n d u ct e d, a n d a d diti o n al p arti ci p a nt s w er e a d d e d s o t h at t h e s h a m c o n diti o n w a s b ett er p o w er e d. W e a s s e s s e d b a s eli n e c o g niti v e a bilit y, g e n d er, tr ai ni n g sit e, a n d m oti v ati o n l e v el a n d f o u n d si g nifi c a nt i nt er a cti o ns b et w e e n b ot h b as eli n e a bilit y a n d m oti v ati o n wit h c o n diti o n ( a cti v e or s h a m) i n m o d els pr e di cti n g tr ai ni n g g ai n. I n a d diti o n, t h e i m pr o v e m e nt s i n t h e a cti v e c o nditi o n v er s u s s h a m c o n diti o n a p p e ar t o b e st a bl e e v e n a s l o n g a s a y e ar aft er t h e ori gi n al i nt er v e nti o n. ■ 
    more » « less