The effect of tidal asymmetry and temporal settling lag on sediment trapping in tidal estuaries

Over decades and centuries, the mean depth of estuaries changes due to sea-level rise, land subsidence, infilling, and dredging projects. These processes produce changes in relative roughness (friction) and mixing, resulting in fundamental changes in the characteristics of the horizontal (velocity) and vertical tides (sea surface elevation) and the dynamics of sediment trapping. To investigate such changes, a 2DV model is developed. The model equations consist of the width-averaged shallow water equations and a sediment balance equation. Together with the condition of morphodynamic equilibrium, these equations are solved analytically by making a regular expansion of the various physical variables in a small parameter. Using these analytic solutions, we are able to gain insight into the fundamental physical processes resulting in sediment trapping in an estuary by studying various forcings separately. As a case study, we consider the Ems estuary. Between 1980 and 2005, successive deepening of the Ems estuary has significantly altered the tidal and sediment dynamics. The tidal range and the surface sediment concentration has increased and the position of the turbidity zone has shifted into the freshwater zone. The model is used to determine the causes of these historical changes. It is found that the increase of the tidal amplitude toward the end of the embayment is the combined effect of the deepening of the estuary and a 37% and 50% reduction in the vertical eddy viscosity and stress parameter, respectively. The physical mechanism resulting in the trapping of sediment, the number of trapping regions, and their sensitivity to grain size are explained by careful analysis of the various contributions of the residual sediment transport. It is found that sediment is trapped in the estuary by a delicate balance between the M ₂ transport and the residual transport for fine sediment ( $\emph{w}_s=0.2$  mm s^???1) and the residual, M ₂ and M ₄ transports for coarser sediment ( $\emph{w}_s=2$  mm s^???1). The upstream movement of the estuarine turbidity maximum into the freshwater zone in 2005 is mainly the result of changes in tidal asymmetry. Moreover, the difference between the sediment distribution for different grain sizes in the same year can be attributed to changes in the temporal settling lag.     

Transverse structure of tidal and residual flow and sediment concentration in estuaries

An analytical and a numerical model are used to understand the response of velocity and sediment distributions over Gaussian-shaped estuarine cross-sections to changes in tidal forcing and water depth. The estuaries considered here are characterized by strong mixing and a relatively weak along-channel density gradient. It is also examined under what conditions the fast, two-dimensional analytical flow model yields results that agree with those obtained with the more complex three-dimensional numerical model. The analytical model reproduces and explains the main velocity and sediment characteristics in large parts of the parameter space considered (average tidal velocity amplitude, 0.1–1 m s^− 1 and maximum water depth, 10–60 m). Its skills are lower for along-channel residual flows if nonlinearities are moderate to high (strong tides in deep estuaries) and for transverse flows and residual sediment concentrations if the Ekman number is small (weak tides in deep estuaries). An important new aspect of the analytical model is the incorporation of tidal variations in the across-channel density gradient, causing a double circulation pattern in the transverse flow during slack tides. The gradient also leads to a new tidally rectified residual flow component via net advection of along-channel tidal momentum by the density-induced transverse tidal flow. The component features landward currents in the channel and seaward currents over the slopes and is particularly effective in deeper water. It acts jointly with components induced by horizontal density differences, Coriolis-induced tidal rectification and Stokes discharge, resulting in different along-channel residual flow regimes. The residual across-channel density gradient is crucial for the residual transverse circulation and for the residual sediment concentration. The clockwise density-induced circulation traps sediment in the fresher water over the left slope (looking up-estuary in the northern hemisphere). Model results are largely consistent with available field data of well-mixed estuaries.     

Tidal flow asymmetry in the diurnal regime: bed-load transport and morphologic changes around the Red River Delta

The relation between tidal flow asymmetry and net transport of sediment in the semidiurnal regime has been extensively described. This study reveals that in the diurnal regime, the direction of long-term net bed-load transport and resulting morphologic changes is partly determined by the phase-angle relationship of O₁, K₁, and M₂. Simple analytical formulations of time-averaged bed-load transport were derived which separate the relative contributions of tidal asymmetry from that of residual flow with tidal stirring. In this particular case, the Red River Delta in Vietnam, transports related to tidal asymmetry are larger than those induced by the monsoon currents, and are an order of magnitude larger than those associated with topographic residual flow around the delta. Tide-induced morphologic changes dominate at water depths between 10 and 25 m, of which the patterns of erosion and deposition overlap with observed bathymetric changes. Additional observed changes that occur in more shallow water cannot be explained by tidal asymmetry and are probably related to wave action and to deposition from the buoyant river plume.Responsible Editor: Jens Kappenberg     

Simulation of the Lagrangian tide-induced residual velocity in a tide-dominated coastal system: a case study of Jiaozhou Bay, China

The Eulerian residual transport velocity and the first-order Lagrangian residual velocity for weakly nonlinear systems have been used extensively in the past to depict inter-tidal mass transport. However, these could not explain the observed net surface sediment transport pattern in Jiaozhou Bay (JZB), located on the western Yellow Sea. JZB is characterized by strong tidal motion, complex topography and an irregular coastline, which are features of typical nonlinear systems. The Lagrangian residual velocity, which is applicable to general nonlinear systems, was simulated with the water parcel tracking method. The results indicate that the composition of the Lagrangian residual velocity at different tidal phases coincides well with the observed net surface sediment transport pattern. The strong dependence of water flushing time on the initial tidal phase can also be explained by the significant intra-tidal variation of the Lagrangian residual velocity. To investigate the hydrodynamic mechanism governing the nonlinearity of the M ₂ tidal system, a set of nonlinearity indexes were defined and analysed. In the surface layer, horizontal advection is the main contributor to the strong nonlinearity near the bay mouth, while in the bottom layer, the strong nonlinearity near the bay mouth may result from the vertical viscosity and horizontal advection.     

Lateral circulation driven by boundary mixing and the associated transport of sediments in idealized partially mixed estuaries

A three-dimensional, hydrostatic, primitive equation numerical model with modern turbulence closures is used to explore lateral circulation and the associated transport of sediments in idealized, moderately to highly stratified estuaries. The model results suggest that boundary mixing on a sloping bottom can drive a significant amount of lateral circulation. This mechanism has received little attention to date in the estuarine literature. Good agreement with an analytical solution and similar vertical structures of lateral flows to observations from the Hudson River estuary support the importance of the boundary mixing mechanism. Boundary mixing is at least as important as differential advection for the modeled scenarios, when the two mechanisms are evaluated using the salt balance equation for model runs without rotation. Linearly superposing analytical solutions for lagged boundary mixing lateral flow and Ekman-forced lateral flow yields a good representation of the near-bottom lateral flow from the model with rotation. The 2 h lag required for the boundary mixing solution is roughly equal to the vertical diffusion time scale, indicating that lateral flow adjustment depends on development of a bottom mixed layer. Sediment dynamics at cross sections seaward and landward of the salt intrusion are very different. Seaward of the salt intrusion, sediments are eroded in the channel and preferentially deposited on the right slope (looking seaward), mainly due to the combination of high sediment concentration in the channel during flood with strong up-slope transport on that side (tidal pumping). Lateral sediment re-distribution landward of the salt intrusion is negligible due to weak residual lateral circulation.     

Effect of stratification on tidal circulation over the Scotian Shelf and Gulf of St. Lawrence: a numerical study using a three-dimensional shelf circulation model

A three-dimensional shelf circulation model is used to examine the effect of seasonal changes in water-column stratification on the tidal circulation over the Scotian Shelf and Gulf of St. Lawrence. The model is driven by tidal forcing specified at the model’s lateral open boundaries in terms of tidal sea surface elevations and depth-averaged currents for five major tidal constituents (M₂, N₂, S₂, K₁, and O₁). Three numerical experiments are conducted to determine the influence of baroclinic pressure gradients and changes in vertical mixing, both associated with stratification, on the seasonal variation of tidal circulation over the study region. The model is initialized with climatological hydrographic fields and integrated for 16 months in each experiment. Model results from the last 12 months are analyzed to determine the dominant semidiurnal and diurnal tidal components, M₂ and K₁. Model results suggest that the seasonal variation in the water-column stratification affects the M₂ tidal circulation most strongly over the shelf break and over the deep waters off the Scotian Shelf (through the development of baroclinic pressure gradients) and along Northumberland Strait in the Gulf of St. Lawrence (through changes in vertical mixing and bottom stress). For the K₁ constituent, the baroclinic pressure gradient and vertical mixing have opposing effects on the tidal circulation over several areas of the study region, while near the bottom, vertical mixing appears to play only a small role in the tidal circulation.     

Numerical study of circulation on the inner Amazon Shelf

We studied the circulation on the coastal domain of the Amazon Shelf by applying the hydrodynamic module of the estuarine and coastal ocean model and sediment transport. The first barotropic experiment aimed to explain the major bathymetric effects on tides and those generated by anisotropy in sediment distribution. We analyzed the continental shelf response of barotropic tides under realistic bottom stress parametrization (C _d ), considering sediment granulometry obtained from a faciologic map, where river mud deposits and reworked sediments areas are well distinguished, among others classes of sediments. Very low C _d values were set in the fluid mud regions off the Amapá coast (1.0 10^???4), in contrast to values around 3.5 10^???3 for coarser sediment regions off the Pará coast. Three-dimensional experiments represented the Amazon River discharge and trade winds, combined to barotropic tide influences and induced vertical mixing. The quasiresonant response of the Amazon Shelf to the M₂ tide acts on the local hydrodynamics by increasing tidal admittance, along with tidal forcing at the shelf break and extensive fluid mud regions. Harmonic analysis of modeled currents agreed well with the analysis of the AMASSEDS observational data set. Tidal-induced vertical shear provided strong homogenization of threshold waters, which are subject to a kind of hydraulic control due to the topographic steepness. Ahead of the hydraulic jump, the low-salinity plume is disconnected from the bottom and acquires negative vorticity, turning southeastward. Tides act as a generator mechanism and topography, via hydraulic control, acts as a maintainer mechanism for the low-salinity frontal zone positioning. Tidally induced southeastward plume fate is overwhelmed by northwestward trade winds so that they, along with background circulation, probably play the most important role on the plume fate and variability over the Amazon Shelf.     

Simulation of the Lagrangian tide-induced residual velocity in a tide-dominated coastal system: a case study of Jiaozhou Bay,China

The Eulerian residual transport velocity and the first-order Lagrangian residual velocity for weakly nonlinear systems have been used extensively in the past to depict inter-tidal mass transport. However, these could not explain the observed net surface sediment transport pattern in Jiaozhou Bay (JZB), located on the western Yellow Sea. JZB is characterized by strong tidal motion, complex topography and an irregular coastline, which are features of typical nonlinear systems. The Lagrangian residual velocity, which is applicable to general nonlinear systems, was simulated with the water parcel tracking method. The results indicate that the composition of the Lagrangian residual velocity at different tidal phases coincides well with the observed net surface sediment transport pattern. The strong dependence of water flushing time on the initial tidal phase can also be explained by the significant intra-tidal variation of the Lagrangian residual velocity. To investigate the hydrodynamic mechanism governing the nonlinearity of the M ₂ tidal system, a set of nonlinearity indexes were defined and analysed. In the surface layer, horizontal advection is the main contributor to the strong nonlinearity near the bay mouth, while in the bottom layer, the strong nonlinearity near the bay mouth may result from the vertical viscosity and horizontal advection.

     

Feedback between residual circulations and sediment distribution in highly turbid estuaries: An analytical model

Motivated by field studies of the Ems estuary which show longitudinal gradients in bottom sediment concentration as high as O(0.01 kg/m⁴), we develop an analytical model for estuarine residual circulation based on currents from salinity gradients, turbidity gradients, and freshwater discharge. Salinity is assumed to be vertically well mixed, while the vertical concentration profile is assumed to result from a balance between a constant settling velocity and turbulent diffusive flux. Width and depth of the model estuary are held constant. Model results show that turbidity gradients enhance tidally averaged circulation upstream of the estuarine turbidity maximum (ETM), but significantly reduce residual circulation downstream, where salinity and turbidity gradients oppose each other. We apply the condition of morphodynamic equilibrium (vanishing sediment transport) and develop an analytical solution for the position of the turbidity maximum and the distribution of suspended sediment concentration (SSC) along a longitudinal axis. A sensitivity study shows great variability in the longitudinal distribution of suspended sediment with the applied salinity gradient and six model parameters: settling velocity, vertical mixing, horizontal dispersion, total sediment supply, fresh water flow, and water depth. Increasing depth and settling velocity move the ETM upstream, while increasing freshwater discharge and vertical mixing move the ETM downstream. Moreover, the longitudinal distribution of SSC is inherently asymmetric around the ETM, and depends on spatial variations in the residual current structure and the vertical profile of SSC.     

Observations of the spring–neap modulation of the gravitational circulation in a partially mixed estuary

The long-term variability of the non-tidal circulation in Southampton Water, a partially mixed estuary, was investigated using 71-day acoustic Doppler current profiler (ADCP) time series. The data show evidence that the spring–neap tidal variability of the turbulent mixing modulates the strength of the non-tidal residual circulation, with subtidal neap tide surface flows reaching 0.12 m s^–1 compared to <0.05 m s^–1 at spring tides. The amplitude of the neap-tide events in this non-tidal circulation is shown to be related to a critical value of the tidal currents, illustrating the strong dependence on tidal mixing. The results suggest that the dominant mechanism for generating these neap-tide circulation events is the baroclinic forcing of the horizontal density gradient, rather than barotropic forcing associated with ebb-induced periodic stratification. While tidal turbulence is thought to be the dominant control on this gravitational circulation, there is evidence of the additional effect of wind-driven mixing, including the effects of wind fetch and possibly wave development with along-estuary winds being more efficient at mixing the estuary than across-estuary winds. Rapid changes in atmospheric pressure also coincided with fluctuations in the gravitational circulation. The observed subtidal flows are shown to be capable of rapidly flushing buoyant material out of the estuary and into the coastal sea at neap tides.Responsible Editor: Iris Grabemann     

The effect of asymmetric dune roughness on tidal asymmetry in the Weser estuary

The bed of estuaries is often characterized by ripples and dunes of varying size. Whereas smaller bedforms adapt their morphological shape to the oscillating tidal currents, large compound dunes (here: asymmetric tidal dunes) remain stable for periods longer than a tidal cycle. Bedforms constitute a form roughness, that is, hydraulic flow resistance, which has a large-scale effect on tidal asymmetry and, hence, on hydrodynamics, sediment transport, and morphodynamics of estuaries and coastal seas. Flow separation behind the dune crest and recirculation on the steep downstream side result in turbulence and energy loss. Since the energy dissipation can be related to the dune lee slope angle, asymmetric dune shapes induce variable flow resistance during ebb and flood phases. Here, a noncalibrated numerical model has been applied to analyze the large-scale effect of symmetric and asymmetric dune shapes on estuarine tidal asymmetry evaluated by residual bed load sediment transport at the Weser estuary, Germany. Scenario simulations were performed with parameterized bed roughness of symmetric and asymmetric dune shapes and without dune roughness. The spatiotemporal interaction of distinct dune shapes with the main drivers of estuarine sediment and morphodynamics, that is, river discharge and tidal energy, is shown to be complex but substantial. The contrasting effects of flood- and ebb-oriented asymmetric dunes on residual bed load transport rates and directions are estimated to be of a similar importance as the controls of seasonal changes of discharge on these net sediment fluxes at the Lower Weser estuary. This corroborates the need to consider dune-induced directional bed roughness in numerical models of estuarine and tidal environments.     

The formation mechanisms of turbidity maximum in the Pearl River estuary, China

The formation and dynamics of turbidity maxima (TM) in the Pearl River estuary (PRE) are not well understood but the existence of TM in the estuary have great potential engineering and environmental impacts. Based on the measurements of two large-scale field surveys conducted in the wet season of 1978 (July of 1978) and the dry season of 1979 (March of 1979), the existence and seasonal variations of TM in the PRE focusing on the two major natural navigation channels have been studied. The sediment transport fluxes over two consecutive tidal cycles have been analyzed in detail. The analysis results reveal that the formation mechanisms of TM in the estuary are rather complex. In general, gravitational circulation, tidal trapping, and sediment resuspension and deposition processes are the primary TM formation mechanisms in the PRE. The clockwise back flow pattern around Lingding Island also leads to the formation of TM in the West channel of the PRE. The occurrence of TM far upstream of the salt water wedge is the result of the complex hydrodynamic and sediment transport processes generated by the runoff of the major rivers.     

Seasonal variation of the <Emphasis Type="Italic">M</Emphasis><Subscript>2</Subscript> tide

The seasonal cycle of the main lunar tidal constituent M ₂ is studied globally by an analysis of a high-resolution ocean circulation and tide model (STORMTIDE) simulation, of 19 years of satellite altimeter data, and of multiyear tide-gauge records. The barotropic seasonal tidal variability is dominant in coastal and polar regions with relative changes of the tidal amplitude of 5–10 %. A comparison with the observations shows that the ocean circulation and tide model captures the seasonal pattern of the M ₂ tide reasonably well. There are two main processes leading to the seasonal variability in the barotropic tide: First, seasonal changes in stratification on the continental shelf affect the vertical profile of eddy viscosity and, in turn, the vertical current profile. Second, the frictional effect between sea-ice and the surface ocean layer leads to seasonally varying tidal transport. We estimate from the model simulation that the M ₂ tidal energy dissipation at the sea surface varies seasonally in the Arctic (ocean regions north of 60°N) between 2 and 34 GW, whereas in the Southern Ocean, it varies between 0.5 and 2 GW. The M ₂ internal tide is mainly affected by stratification, and the induced modified phase speed of the internal waves leads to amplitude differences in the surface tide signal of 0.005–0.0150 m. The seasonal signals of the M ₂ surface tide are large compared to the accuracy demands of satellite altimetry and gravity observations and emphasize the importance to consider seasonal tidal variability in the correction processes of satellite data.     

An idealized model study of flocculation on sediment trapping in an estuarine turbidity maximum

Flocculation has an important impact on particle trapping in estuarine turbidity maximum (ETM) through associated increases in particle settling velocity. To quantify the importance of the flocculation processes, a size-resolved flocculation model is implemented into an ocean circulation model to simulate fine-grained particle trapping in an ETM. The model resolves the particle size from robust small flocs, about 30 μm, to very large flocs, over 1000 μm. An idealized two-dimensional model study is performed to simulate along-channel variations of suspended sediment concentrations driven by gravitational circulation and tidal currents. The results indicate that the flocculation processes play a key role in generating strong tidal asymmetrical variations in suspended sediment concentration and particle trapping. Comparison with observations suggests that the flocculation model produces realistic characteristics of an ETM.     

Numerical experiments on the long-term morphodynamics of the Colorado River Delta

The interaction among tidal currents, sediment transport, and long-term changes of the sea bottom in the Colorado River Delta have been investigated applying a two dimensional nonlinear hydrodynamic-numerical finite differences model. The system was forced by the dominant M ₂ tidal component at the open boundary. We carried out calculations to study the morphodynamics of the actual bathymetry caused by the bedload sediment transport. To investigate the origin of actual morphological features, we performed experiments using a smoothed bathymetry, in which the islands Montague, Gore, and Pelícano were eliminated. Under the imposed tidal hydrodynamics, the results indicate that the bedload transport contributes significantly in the genesis of sandbanks and in the formation and maintaining of the Montague and Gore Islands.     

Sensitivity of tidal characteristics in double inlet systems to momentum dissipation on tidal flats: a perturbation analysis

In a tidal channel with adjacent tidal flats, along–channel momentum is dissipated on the flats during rising tides. This leads to a sink of along–channel momentum. Using a perturbative method, it is shown that the momentum sink slightly reduces the M₂ amplitude of both the sea surface elevation and current velocity and favours flood dominant tides. These changes in tidal characteristics (phase and amplitude of sea surface elevations and currents) are noticeable if widths of tidal flats are at least of the same order as the channel width, and amplitudes and gradients of along–channel velocity are large. The M₂ amplitudes are reduced because stagnant water flows from the flats into the channel, thereby slowing down the current. The M₄ amplitudes and phases change because the momentum sink acts as an advective term during the fall of the tide, such a term generates flood dominant currents. For a prototype embayment that resembles the Marsdiep–Vlie double–inlet system of the Western Wadden Sea, it is found that for both the sea surface elevation and current velocity, including the momentum sink, lead to a decrease of approximately 2% in M₂ amplitudes and an increase of approximately 25% in M₄ amplitudes. As a result, the net import of coarse sediment is increased by approximately 35%, while the transport of fine sediment is hardly influenced by the momentum sink. For the Marsdiep–Vlie system, the M₂ sea surface amplitude obtained from the idealised model is similar to that computed with a realistic three–dimensional numerical model whilst the comparison with regard to M₄ improves if momentum sink is accounted for.     

Growth of large-scale bed forms due to storm-driven and tidal currents: a model approach

An idealized morphodynamic model is used to gain further understanding about the formation and characteristics of shoreface-connected sand ridges and tidal sand banks on the continental shelf. The model consists of the 2D shallow water equations, supplemented with a sediment transport formulation and describes the initial feedback between currents and small amplitude bed forms. The behaviour of bed forms during both storm and fair weather conditions is analyzed. This is relevant in case of coastal seas characterized by tidal motion, where the latter causes continuous transport of sediment as bed load.The new aspects of this work are the incorporation of both steady and tidal currents (represented by an M₂ and M₄ component) in the external forcing, in combination with dominant suspended sediment transport during storms. The results indicate that the dynamics during storms and fair weather strongly differ, causing different types of bed forms to develop. Shoreface-connected sand ridges mainly form during storm conditions, whereas if fair weather conditions prevail the more offshore located tidal sand banks develop. Including the M₄ tide changes the properties of the bed forms, such as growth rates and migration speeds, due to tidal asymmetry. Finally a probabilistic formulation of the storm and fair weather realization of the model is used to find conditions for which both types of large-scale bed forms occur simultaneously. These conditions turn out to be a low storm fraction and the presence strong tidal currents in combination with strong steady currents during storms.     

Direct measurement of dissolved N2 and denitrification along a subtropical river-estuary gradient,China

The spatial pattern and seasonal variation of denitrification were investigated during 2010–2011 in the Jiulong River Estuary (JRE) in southeast China. Dissolved N₂ was directly measured by changes in the N₂:Ar ratio. The results showed that excess dissolved N₂ ranged from ?9.9 to 76.4 μmol L^?1. Tidal mixing leads to a seaward decline of dissolved gaseous concentrations and water–air fluxes along the river-estuary gradient. Denitrification at freshwater sites varied between seasons, associated with changes in N input and water temperature. The denitrification process was controlled by the nitrate level at freshwater sites, and the excess dissolved N₂ observed at the tidal zone largely originated from upstream water transport. Compared to other estuaries, JRE has a relative low gaseous removal efficiency (E_d = 12% of [DIN]; annual N removal = 24% of DIN load), a fact ascribed to strong tidal mixing, coarse-textured sediment with shallow depth before bedrock and high riverine DIN input.     

Long-term evolution of sand waves in the Marsdiep inlet. II: Relation to hydrodynamics

A discussion is presented about the mechanisms that govern the spatial and seasonal variability in sand-wave height and migration speed in the 4 km wide Marsdiep tidal inlet, the Netherlands. Since 1998, current velocities and water depths have been recorded with an ADCP that is mounted under the ferry ‘Schulpengat’. In this paper, the current measurements were used to explain the sand-wave observations presented in Buijsman and Ridderinkhof [this issue. Long-term evolution of sand waves in the Marsdiep inlet. I: high-resolution observations. Continental Shelf Research, doi:10.1016/j.csr.2007.10.011]. Across nearly the entire inlet, the sand waves migrate in the flood direction. In the flood-dominated southern part of the inlet, the ‘measured’ (i.e. based on sand-wave shape and migration speed) and predicted bedload transport agree in direction, magnitude, and trends, whereas in the ebb-dominated northern part the predicted bedload and suspended load transport is opposite to the sand-wave migration. In the southern part, 55% of the bedload transport is due to tidal asymmetries and 45% due to residual currents. In addition to the well-known tidal asymmetries, asymmetries that arise from the interaction of _M2${\mathrm{M}}_2$ and its overtides with _S2${\mathrm{S}}_2$ and its compound tides are also important. It is hypothesised that in the northern part of the inlet the advection of suspended sand and lag effects govern the sand-wave migration. The relative importance of suspended load transport also explains why the sand waves have smaller lee-slope angles, are smaller, more rounded, and more three-dimensional in the northern half of the inlet. The sand waves in this part of the inlet feature the largest seasonal variability in height and migration speed. This seasonal variability may be attributed to the tides or a seasonal fluctuation in fall velocity. In both cases sediment transport is enhanced in winter, increasing sand-wave migration and decreasing sand-wave height. The influence of storms and estuarine circulation on the sand-wave variability is negligible.     

Modeling studies of tidal dynamics and the associated responses to coastline changes in the Bohai Sea,China

Over the past 30 years, reclamation projects and related changes have impacted the hydrodynamics and sediment transport in the Bohai Sea. Three-dimensional tidal current models of the Bohai Sea and the Yellow Sea were constructed using the MIKE 3 model. We used a refined grid to simulate and analyze the effects of changes in coastline, depth, topography, reclamation, the Yellow River estuary, and coastal erosion on tidal systems, tide levels, tidal currents, residual currents, and tidal fluxes. The simulation results show that the relative change in the amplitude of the half-day tide is greater than that of the full-day tide. The changes in the tidal amplitudes of M₂, S₂, K₁, and O₁ caused by coastline changes accounted for 27.76–99.07% of the overall change in amplitude from 1987 to 2016, and water depth changes accounted for 0.93–72.24% of the overall change. The dominant factor driving coastline changes is reclamation, accounting for 99.55–99.91% of the amplitude changes in tidal waves, followed by coastal erosion, accounting for 0.05–0.40% of the tidal wave amplitude changes. The contribution of changes in the Yellow River estuary to tidal wave amplitude changes is small, accounting for 0.01–0.12% of the amplitude change factor. The change in the highest tide level (HTL) is mainly related to the amplitude change, and the correlation with the phase change is small. The dominant factor responsible for the change in the HTL is the tide amplitude change in M2, followed by S2, whereas the influence of the K1 and O1 tides on the change in the HTL is small. Reclamation resulted in a decrease in the vertical average maximum flow velocity (V_VAM) in the Bohai Sea. Shallower water depths have led to an increase in the V_VAM; deeper water depths have led to a decrease in the maximum flow velocity. The absolute value of the maximum flow velocity gradually decreases from the surface to the bottom, but the relative change value is basically constant. The changes in the tidal dynamics of the Bohai Sea are proportional to the degree of change in the coastline. The maximum and minimum changes in the tidal flux appear in Laizhou Bay (P-LZB) and Liaodong Bay (P-LDB), respectively. The changes in the tidal flux are related to the change in the area of the bay. Due to the reduced tidal flux, the water exchange capacity of the Bohai Sea has decreased, impacting the ecological environment of the Bohai Sea. Strictly controlling the scale of reclamation are important measures for reducing the decline in the water exchange capacity of the Bohai Sea and the deterioration of its ecological environment.