Search Results
You are looking at 1 - 10 of 13 items for
- Author or Editor: Henk M. Schuttelaars x
- Refine by Access: All Content x
Abstract
The salinity structure in estuaries is classically described in terms of the salinity structure as well mixed, partially mixed, or salt wedge. The existing knowledge about the processes that result in such salinity structures comes from highly idealized models that are restricted to either well-mixed and partially mixed cases or subtidal salt wedge estuaries. Hence, there is still little knowledge about the processes driving transitions between these different salinity structures and the estuarine parameters at which such a transition is found. As an important step toward a unified description of the dominant processes driving well-mixed, partially mixed, and salt wedge estuaries, a subtidal width-averaged model applicable to all these salinity structures is developed and systematically analyzed. Using our model, we identify four salinity regimes, resulting from different balances of dominant processes. It is shown that each regime is uniquely determined by two dimensionless parameters: an estuarine Froude and Rayleigh number, representing freshwater discharge and tidal mixing, respectively, resulting in a classification of the regimes in terms of these two parameters. Furthermore, analytical expressions to approximate the salt intrusion length in each regime are developed. These expressions are used to illustrate that the salt intrusion length in different regimes responds in a highly different manner to changes in depth and freshwater discharge. As one of the key results, we show that there are only very weak relations between the process-based regime of an estuary and the salt intrusion length and top–bottom stratification. This implies that the salinity structure of an estuary cannot be uniquely matched to a regime.
Abstract
The salinity structure in estuaries is classically described in terms of the salinity structure as well mixed, partially mixed, or salt wedge. The existing knowledge about the processes that result in such salinity structures comes from highly idealized models that are restricted to either well-mixed and partially mixed cases or subtidal salt wedge estuaries. Hence, there is still little knowledge about the processes driving transitions between these different salinity structures and the estuarine parameters at which such a transition is found. As an important step toward a unified description of the dominant processes driving well-mixed, partially mixed, and salt wedge estuaries, a subtidal width-averaged model applicable to all these salinity structures is developed and systematically analyzed. Using our model, we identify four salinity regimes, resulting from different balances of dominant processes. It is shown that each regime is uniquely determined by two dimensionless parameters: an estuarine Froude and Rayleigh number, representing freshwater discharge and tidal mixing, respectively, resulting in a classification of the regimes in terms of these two parameters. Furthermore, analytical expressions to approximate the salt intrusion length in each regime are developed. These expressions are used to illustrate that the salt intrusion length in different regimes responds in a highly different manner to changes in depth and freshwater discharge. As one of the key results, we show that there are only very weak relations between the process-based regime of an estuary and the salt intrusion length and top–bottom stratification. This implies that the salinity structure of an estuary cannot be uniquely matched to a regime.
Abstract
The classification diagram developed by Hansen and Rattray is one of the classic papers on classification of estuarine salinity dynamics. However, we found several inconsistencies in both their stratification–circulation and estuarine classification diagrams. These findings considerably change the interpretation of their work. Furthermore, while their classification includes salt wedge estuaries, the model used to derive this is only applicable to well-mixed and partially mixed estuaries. Here, we identify and solve these inconsistencies, and we propose new adjusted and extended stratification–circulation and classification diagrams. To this end, we summarize the model of Hansen and Rattray and extend their work to find analytical model solutions and an adjusted stratification–circulation diagram. Using this new diagram, it is shown that Hansen and Rattray incorrectly discussed the behavior of dispersion-dominated estuaries and that several parts of the diagram correspond to physically unrealistic model solutions. This is then used to demonstrate that several estuarine classes identified by Hansen and Rattray correspond to physically unrealistic model solutions and can therefore not be interpreted. A new and extended classification is proposed by using a recently developed model that extends the work of Hansen and Rattray to salt wedge estuaries. This results in an extended estuarine classification including examples of the location of 12 estuaries in this new diagram.
Abstract
The classification diagram developed by Hansen and Rattray is one of the classic papers on classification of estuarine salinity dynamics. However, we found several inconsistencies in both their stratification–circulation and estuarine classification diagrams. These findings considerably change the interpretation of their work. Furthermore, while their classification includes salt wedge estuaries, the model used to derive this is only applicable to well-mixed and partially mixed estuaries. Here, we identify and solve these inconsistencies, and we propose new adjusted and extended stratification–circulation and classification diagrams. To this end, we summarize the model of Hansen and Rattray and extend their work to find analytical model solutions and an adjusted stratification–circulation diagram. Using this new diagram, it is shown that Hansen and Rattray incorrectly discussed the behavior of dispersion-dominated estuaries and that several parts of the diagram correspond to physically unrealistic model solutions. This is then used to demonstrate that several estuarine classes identified by Hansen and Rattray correspond to physically unrealistic model solutions and can therefore not be interpreted. A new and extended classification is proposed by using a recently developed model that extends the work of Hansen and Rattray to salt wedge estuaries. This results in an extended estuarine classification including examples of the location of 12 estuaries in this new diagram.
Abstract
Tidal straining, which can mathematically be described as the covariance between eddy viscosity and vertical shear of the along-channel velocity component, has been acknowledged as one of the major drivers for estuarine circulation in channelized tidally energetic estuaries. In this paper, the authors investigate the role of lateral circulation for generating this covariance. Five numerical experiments are carried out, starting with a reference scenario including the full physics and four scenarios in which specific key physical processes are neglected. These processes are longitudinal internal pressure gradient forcing, lateral internal pressure gradient forcing, lateral advection, and the neglect of temporal variation of eddy viscosity. The results for the viscosity–shear covariance are correlated across different experiments to quantify the change due to neglect of these key processes. It is found that the lateral advection of vertical shear of the along-channel velocity component and its interaction with the tidally asymmetric eddy viscosity (which is also modified by the lateral circulation) is the major driving force for estuarine circulation in well-mixed tidal estuaries.
Abstract
Tidal straining, which can mathematically be described as the covariance between eddy viscosity and vertical shear of the along-channel velocity component, has been acknowledged as one of the major drivers for estuarine circulation in channelized tidally energetic estuaries. In this paper, the authors investigate the role of lateral circulation for generating this covariance. Five numerical experiments are carried out, starting with a reference scenario including the full physics and four scenarios in which specific key physical processes are neglected. These processes are longitudinal internal pressure gradient forcing, lateral internal pressure gradient forcing, lateral advection, and the neglect of temporal variation of eddy viscosity. The results for the viscosity–shear covariance are correlated across different experiments to quantify the change due to neglect of these key processes. It is found that the lateral advection of vertical shear of the along-channel velocity component and its interaction with the tidally asymmetric eddy viscosity (which is also modified by the lateral circulation) is the major driving force for estuarine circulation in well-mixed tidal estuaries.
Abstract
A semianalytical three-dimensional model is set up to dynamically calculate the coupled water motion and salinity for idealized well-mixed estuaries and prognostically investigate the influence of each physical mechanism on the residual salt transport. As a study case, a schematized estuary with an exponentially converging width and a channel–shoal structure is considered. The temporal correlation between horizontal tidal velocities and tidal salinities is the dominant process for the landward residual salt transport. The residual salt transport induced by residual circulation is locally significant, but the induced salt transport integrated over the cross section is small. The impacts of the estuarine geometry, Coriolis force, and bathymetry on the salt dynamics are studied using three dedicated experiments, in which the impact of each of these factors is studied separately. To assess the impact of width convergence, a convergent estuary without bathymetric variations or Coriolis force is considered. In this experiment, the temporal correlation between tidal velocities and salinities is the only landward salt transport process. In the second experiment, Coriolis effects are included. This results in a significant residual salt transport cell due to the advection of the tidally averaged salinity by residual circulation, with salt imported into the estuary from the left side and exported on the right (looking seaward). In the last experiment, a lateral channel–shoal structure is included while the Coriolis effects are excluded. This results in a significant landward salt transport through the deeper channel and a seaward salt transport over the shoals due to the advection of the tidally averaged salinity by residual circulation.
Abstract
A semianalytical three-dimensional model is set up to dynamically calculate the coupled water motion and salinity for idealized well-mixed estuaries and prognostically investigate the influence of each physical mechanism on the residual salt transport. As a study case, a schematized estuary with an exponentially converging width and a channel–shoal structure is considered. The temporal correlation between horizontal tidal velocities and tidal salinities is the dominant process for the landward residual salt transport. The residual salt transport induced by residual circulation is locally significant, but the induced salt transport integrated over the cross section is small. The impacts of the estuarine geometry, Coriolis force, and bathymetry on the salt dynamics are studied using three dedicated experiments, in which the impact of each of these factors is studied separately. To assess the impact of width convergence, a convergent estuary without bathymetric variations or Coriolis force is considered. In this experiment, the temporal correlation between tidal velocities and salinities is the only landward salt transport process. In the second experiment, Coriolis effects are included. This results in a significant residual salt transport cell due to the advection of the tidally averaged salinity by residual circulation, with salt imported into the estuary from the left side and exported on the right (looking seaward). In the last experiment, a lateral channel–shoal structure is included while the Coriolis effects are excluded. This results in a significant landward salt transport through the deeper channel and a seaward salt transport over the shoals due to the advection of the tidally averaged salinity by residual circulation.
Abstract
The dependency of the estuarine circulation on the depth-to-width ratio of a periodically, weakly stratified tidal estuary is systematically investigated here for the first time. Currents, salinity, and other properties are simulated by means of the General Estuarine Transport Model (GETM) in cross-sectional slice mode, applying a symmetric Gaussian-shaped depth profile. The width is varied over four orders of magnitude. The individual along-channel circulation contributions from tidal straining, gravitation, advection, etc., are calculated and the impact of the depth-to-width ratio on their intensity is presented and elucidated. It is found that the estuarine circulation exhibits a distinct maximum in medium-wide channels (intermediate depth-to-width ratio depending on various parameters), which is caused by a maximum of the tidal straining contribution. This maximum is related to a strong tidal asymmetry of eddy viscosity and shear created by secondary strain-induced periodic stratification (2SIPS): in medium channels, transverse circulation generated by lateral density gradients due to laterally differential longitudinal advection induces stable stratification at the end of the flood phase, which is further increased during ebb by longitudinal straining (SIPS). Thus, eddy viscosity is low and shear is strong in the entire ebb phase. During flood, SIPS decreases the stratification so that eddy viscosity is high and shear is weak. The circulation resulting from this viscosity–shear correlation, the tidal straining circulation, is oriented like the classical, gravitational circulation, with riverine outflow at the surface and oceanic inflow close to the bottom. In medium channels, it is about 5 times as strong as in wide (quasi one-dimensional) channels, in which 2SIPS is negligible.
Abstract
The dependency of the estuarine circulation on the depth-to-width ratio of a periodically, weakly stratified tidal estuary is systematically investigated here for the first time. Currents, salinity, and other properties are simulated by means of the General Estuarine Transport Model (GETM) in cross-sectional slice mode, applying a symmetric Gaussian-shaped depth profile. The width is varied over four orders of magnitude. The individual along-channel circulation contributions from tidal straining, gravitation, advection, etc., are calculated and the impact of the depth-to-width ratio on their intensity is presented and elucidated. It is found that the estuarine circulation exhibits a distinct maximum in medium-wide channels (intermediate depth-to-width ratio depending on various parameters), which is caused by a maximum of the tidal straining contribution. This maximum is related to a strong tidal asymmetry of eddy viscosity and shear created by secondary strain-induced periodic stratification (2SIPS): in medium channels, transverse circulation generated by lateral density gradients due to laterally differential longitudinal advection induces stable stratification at the end of the flood phase, which is further increased during ebb by longitudinal straining (SIPS). Thus, eddy viscosity is low and shear is strong in the entire ebb phase. During flood, SIPS decreases the stratification so that eddy viscosity is high and shear is weak. The circulation resulting from this viscosity–shear correlation, the tidal straining circulation, is oriented like the classical, gravitational circulation, with riverine outflow at the surface and oceanic inflow close to the bottom. In medium channels, it is about 5 times as strong as in wide (quasi one-dimensional) channels, in which 2SIPS is negligible.
Abstract
In this idealized numerical modeling study, the composition of residual sediment fluxes in energetic (e.g., weakly or periodically stratified) tidal estuaries is investigated by means of one-dimensional water column models, with some focus on the sediment availability. Scaling of the underlying dynamic equations shows dependence of the results on the Simpson number (relative strength of horizontal density gradient) and the Rouse number (relative settling velocity) as well as impacts of the Unsteadiness number (relative tidal frequency). Here, the parameter space given by the Simpson and Rouse numbers is mainly investigated. A simple analytical model based on the assumption of stationarity shows that for small Simpson and Rouse numbers sediment flux is down estuary and vice versa for large Simpson and Rouse numbers. A fully dynamic water column model coupled to a second-moment turbulence closure model allows to decompose the sediment flux profiles into contributions from the transport flux (product of subtidal velocity and sediment concentration profiles) and the fluctuation flux profiles (tidal covariance between current velocity and sediment concentration). Three different types of bottom sediment pools are distinguished to vary the sediment availability, by defining a time scale for complete sediment erosion. For short erosion times scales, the transport sediment flux may dominate, but for larger erosion time scales the fluctuation sediment flux largely dominates the tidal sediment flux. When quarter-diurnal components are added to the tidal forcing, up-estuary sediment fluxes are strongly increased for stronger and shorter flood tides and vice versa. The theoretical results are compared to field observations in a tidally energetic inlet.
Abstract
In this idealized numerical modeling study, the composition of residual sediment fluxes in energetic (e.g., weakly or periodically stratified) tidal estuaries is investigated by means of one-dimensional water column models, with some focus on the sediment availability. Scaling of the underlying dynamic equations shows dependence of the results on the Simpson number (relative strength of horizontal density gradient) and the Rouse number (relative settling velocity) as well as impacts of the Unsteadiness number (relative tidal frequency). Here, the parameter space given by the Simpson and Rouse numbers is mainly investigated. A simple analytical model based on the assumption of stationarity shows that for small Simpson and Rouse numbers sediment flux is down estuary and vice versa for large Simpson and Rouse numbers. A fully dynamic water column model coupled to a second-moment turbulence closure model allows to decompose the sediment flux profiles into contributions from the transport flux (product of subtidal velocity and sediment concentration profiles) and the fluctuation flux profiles (tidal covariance between current velocity and sediment concentration). Three different types of bottom sediment pools are distinguished to vary the sediment availability, by defining a time scale for complete sediment erosion. For short erosion times scales, the transport sediment flux may dominate, but for larger erosion time scales the fluctuation sediment flux largely dominates the tidal sediment flux. When quarter-diurnal components are added to the tidal forcing, up-estuary sediment fluxes are strongly increased for stronger and shorter flood tides and vice versa. The theoretical results are compared to field observations in a tidally energetic inlet.
Abstract
Understanding salt dynamics is important to adequately model salt intrusion, baroclinic forcing, and sediment transport. In this paper, the importance of the residual salt transport due to tidal advection in well-mixed tidal estuaries is studied. The water motion is resolved in a consistent way with a width-averaged analytical model, coupled to an advection–diffusion equation describing the salt dynamics. The residual salt balance obtained from the coupled model shows that the seaward salt transport driven by river discharge is balanced by the landward salt transport due to tidal advection and horizontal diffusion. It is found that the tidal advection behaves as a diffusion process, and this contribution is named tidal advective diffusion. The horizontal diffusion parameterizes processes not explicitly resolved in the model and is called the prescribed diffusion. The tidal advective diffusion results from the correlation between the tidal velocity and salinity and can be explicitly calculated with the dominant semidiurnal water motion. The sensitivity analysis shows that tidal advective diffusivity increases with increasing bed roughness and decreasing vertical eddy viscosity. Furthermore, tidal advective diffusivity reaches its maximum for moderate water depth and moderate convergence length. The relative importance of tidal advective diffusion is investigated using the residual salt balance, with the prescribed diffusion coefficient obtained from the measured salinity field. The tidal advective diffusion dominates the residual salt transport in the Scheldt estuary, and other processes significantly contribute to the residual salt transport in the Delaware estuary and the Columbia estuary.
Abstract
Understanding salt dynamics is important to adequately model salt intrusion, baroclinic forcing, and sediment transport. In this paper, the importance of the residual salt transport due to tidal advection in well-mixed tidal estuaries is studied. The water motion is resolved in a consistent way with a width-averaged analytical model, coupled to an advection–diffusion equation describing the salt dynamics. The residual salt balance obtained from the coupled model shows that the seaward salt transport driven by river discharge is balanced by the landward salt transport due to tidal advection and horizontal diffusion. It is found that the tidal advection behaves as a diffusion process, and this contribution is named tidal advective diffusion. The horizontal diffusion parameterizes processes not explicitly resolved in the model and is called the prescribed diffusion. The tidal advective diffusion results from the correlation between the tidal velocity and salinity and can be explicitly calculated with the dominant semidiurnal water motion. The sensitivity analysis shows that tidal advective diffusivity increases with increasing bed roughness and decreasing vertical eddy viscosity. Furthermore, tidal advective diffusivity reaches its maximum for moderate water depth and moderate convergence length. The relative importance of tidal advective diffusion is investigated using the residual salt balance, with the prescribed diffusion coefficient obtained from the measured salinity field. The tidal advective diffusion dominates the residual salt transport in the Scheldt estuary, and other processes significantly contribute to the residual salt transport in the Delaware estuary and the Columbia estuary.
Abstract
The generation of residual circulation in a tidally energetic estuary with constant longitudinal salinity gradient and parabolic cross section is examined by means of a two-dimensional cross-sectional numerical model, neglecting river runoff and Stokes drift. It is shown how the longitudinal and lateral residual circulation can be decomposed into contributions from various processes such as tidal straining circulation, gravitational circulation, advectively driven circulation, and horizontal mixing circulation. The sensitivity of the residual circulation and its components from various processes to changes in forcing is investigated by varying the Simpson number (nondimensional longitudinal buoyancy gradient) and the unsteadiness parameter (nondimensional tidal frequency), as well as the bed roughness and the width of the estuary. For relatively weak salinity gradient forcing, the tidal straining circulation dominates the residual exchange circulation in support of classical estuarine circulation (up-estuary flow near the bed and down-estuary flow near the surface). The strength of the longitudinal estuarine circulation clearly increases with increased salinity gradient forcing. However, when the Simpson number exceeds 0.15, the relative contributions of both gravitational circulation and advectively driven circulation to estuarine circulation increase substantially. Lateral residual circulation is relatively weak for small Simpson numbers and becomes flood oriented (divergent flow near the bed and convergent flow near the surface) for larger Simpson numbers because of increasing contributions from gravitational and advectively driven circulation. Increasing the unsteadiness number leads to decreased longitudinal and lateral residual circulation. Although changes in bed roughness result in relatively small changes in residual circulation, results are sensitive to the width of the estuary, mainly because of changes in residual exchange circulation driven by tidal straining.
Abstract
The generation of residual circulation in a tidally energetic estuary with constant longitudinal salinity gradient and parabolic cross section is examined by means of a two-dimensional cross-sectional numerical model, neglecting river runoff and Stokes drift. It is shown how the longitudinal and lateral residual circulation can be decomposed into contributions from various processes such as tidal straining circulation, gravitational circulation, advectively driven circulation, and horizontal mixing circulation. The sensitivity of the residual circulation and its components from various processes to changes in forcing is investigated by varying the Simpson number (nondimensional longitudinal buoyancy gradient) and the unsteadiness parameter (nondimensional tidal frequency), as well as the bed roughness and the width of the estuary. For relatively weak salinity gradient forcing, the tidal straining circulation dominates the residual exchange circulation in support of classical estuarine circulation (up-estuary flow near the bed and down-estuary flow near the surface). The strength of the longitudinal estuarine circulation clearly increases with increased salinity gradient forcing. However, when the Simpson number exceeds 0.15, the relative contributions of both gravitational circulation and advectively driven circulation to estuarine circulation increase substantially. Lateral residual circulation is relatively weak for small Simpson numbers and becomes flood oriented (divergent flow near the bed and convergent flow near the surface) for larger Simpson numbers because of increasing contributions from gravitational and advectively driven circulation. Increasing the unsteadiness number leads to decreased longitudinal and lateral residual circulation. Although changes in bed roughness result in relatively small changes in residual circulation, results are sensitive to the width of the estuary, mainly because of changes in residual exchange circulation driven by tidal straining.
Abstract
Using an idealized width-averaged process-based model, the role of a mud pool on the bed and time-varying river discharge on the trapping of fine sediment is systematically investigated. For this purpose, a dynamically and physically motivated description of erodibility is presented, which relates the amount of sediment on the bed to the suspended sediment concentration (SSC). We can distinguish between two states: in the availability-limited state, the SSC is limited by the amount of erodible sediment at the bed. Over time, under constant forcing conditions, the estuary evolves to morphodynamic equilibrium. In the erosion-limited state, there is an abundant amount of sediment at the bed so that sediment pickup occurs at the maximum possible rate. The SSC is then limited by the local hydrodynamic conditions. In this state, the estuary keeps importing sediment, forming an erodible bottom pool that grows in time. These two states can be used to explain the response of an estuary to changing river discharge. Under availability-limited conditions, periods of high river discharge push estuarine turbidity maxima (ETMs) downstream, while drier periods allow ETMs to move upstream. However, when the estuary is in an erosion-limited state during low river discharge, a bottom pool is formed. When the discharge then increases, it takes time to deplete this pool, so that an ETM located over a bottom pool moves with a significant time lag relative to changes in the river discharge. Good qualitative agreement is found between model results and observations in the Scheldt Estuary of surface SSC using a representative year of discharge conditions.
Abstract
Using an idealized width-averaged process-based model, the role of a mud pool on the bed and time-varying river discharge on the trapping of fine sediment is systematically investigated. For this purpose, a dynamically and physically motivated description of erodibility is presented, which relates the amount of sediment on the bed to the suspended sediment concentration (SSC). We can distinguish between two states: in the availability-limited state, the SSC is limited by the amount of erodible sediment at the bed. Over time, under constant forcing conditions, the estuary evolves to morphodynamic equilibrium. In the erosion-limited state, there is an abundant amount of sediment at the bed so that sediment pickup occurs at the maximum possible rate. The SSC is then limited by the local hydrodynamic conditions. In this state, the estuary keeps importing sediment, forming an erodible bottom pool that grows in time. These two states can be used to explain the response of an estuary to changing river discharge. Under availability-limited conditions, periods of high river discharge push estuarine turbidity maxima (ETMs) downstream, while drier periods allow ETMs to move upstream. However, when the estuary is in an erosion-limited state during low river discharge, a bottom pool is formed. When the discharge then increases, it takes time to deplete this pool, so that an ETM located over a bottom pool moves with a significant time lag relative to changes in the river discharge. Good qualitative agreement is found between model results and observations in the Scheldt Estuary of surface SSC using a representative year of discharge conditions.
Abstract
An idealized width-averaged model is employed to study the influence of wind stress on subtidal salt intrusion and stratification in well-mixed and partially stratified estuaries. We show that even in mild conditions, wind forcing can influence the estuarine salinity structure in a substantial way. By studying the role of wind forcing on dominant salt transport balances and associated salt transport regimes, we unify and clarify ambiguous observations from previous authors regarding the influence of wind stress: the response of the estuarine salinity structure to wind forcing is different depending on the underlying dominant salt transport balance, which in turn was found to determine whether wind-induced salinity shear, wind-induced modulation of the longitudinal salt distribution, or wind-induced mixing dominates.
Significance Statement
The purpose of this idealized study is to better understand how wind influences the salinity distribution in estuaries on large time scales. This is important because a change in winds can move saline water further inland, threatening freshwater availability and the natural balance of delicate ecosystems. We clarify the sometimes ambiguous observations regarding the influence of wind on the salt distribution and highlight the importance of including average wind forcing in analyses of estuarine dynamics on large time scales.
Abstract
An idealized width-averaged model is employed to study the influence of wind stress on subtidal salt intrusion and stratification in well-mixed and partially stratified estuaries. We show that even in mild conditions, wind forcing can influence the estuarine salinity structure in a substantial way. By studying the role of wind forcing on dominant salt transport balances and associated salt transport regimes, we unify and clarify ambiguous observations from previous authors regarding the influence of wind stress: the response of the estuarine salinity structure to wind forcing is different depending on the underlying dominant salt transport balance, which in turn was found to determine whether wind-induced salinity shear, wind-induced modulation of the longitudinal salt distribution, or wind-induced mixing dominates.
Significance Statement
The purpose of this idealized study is to better understand how wind influences the salinity distribution in estuaries on large time scales. This is important because a change in winds can move saline water further inland, threatening freshwater availability and the natural balance of delicate ecosystems. We clarify the sometimes ambiguous observations regarding the influence of wind on the salt distribution and highlight the importance of including average wind forcing in analyses of estuarine dynamics on large time scales.
