Modeling the influence of river discharge on salt intrusion and residual circulation in Danshuei River estuary,Taiwan

A 3-D, time-dependent, baroclinic, hydrodynamic and salinity model was implemented and applied to the Danshuei River estuarine system and the adjacent coastal sea in Taiwan. The model forcing functions consist of tidal elevations along the open boundaries and freshwater inflows from the main stream and major tributaries in the Danshuei River estuarine system. The bottom friction coefficient was adjusted to achieve model calibration and verification in model simulations of barotropic and baroclinic flows. The turbulent diffusivities were ascertained through comparison of simulated salinity time series with observations. The model simulation results are in qualitative agreement with the available field data.     

Numerical determination of residence time and age in a partially mixed estuary using three-dimensional hydrodynamic model

This paper documents a numerical modeling study to calculate the residence time and age of dissolved substances in a partially mixed estuary. A three-dimensional, time-dependent hydrodynamic model was established and applied to the Danshuei River estuarine system and adjacent coastal sea in Taiwan. The model showed good agreement with observations of surface elevation, tidal currents and salinity made in 2002. The model was then applied to calculate the residence time and age distribution response to different freshwater discharges with and without density-induced circulations in the Danshuei River estuarine system. Regression analysis of model results reveals that an exponential equation can be used to correlate the residence time to change of freshwater input. The simulated results show it takes approximately 10, 4.5, and 3 days, respectively, for a water parcel that has entered the headwaters of the estuary to be transported out of the estuary under low, mean, and high flow conditions with density-induced circulation. The calculated age with density-induced circulation is less than that without density-induced circulation. The age of the surface layer is less than that at the bottom layer. Overall the study shows that freshwater discharges are the important factors in controlling the transport of dissolved substances in the Danshuei River estuarine system.     

Effects of channel regulation on salt intrusion and residual circulation of Keelung River

A vertical (laterally integrated) two‐dimensional numerical model was applied to study the hydrodynamic characteristics, salt‐water intrusion and residual circulation in the Danshuei River estuarine system. The cross‐sectional profiles measured in 2001 and 1990 respectively represent the conditions after and before channel regulation in the Keelung River. The model was re‐verified with the available hydrological data measured in 2001. Detailed model re‐verification has been conducted with water surface elevations, tidal current, and salinity distributions measured. The overall performance of the model is in qualitative agreement with the available field data. The model was then used to investigate the change in tidal ranges, salt‐water intrusion, and residual circulation as a result of channel regulation in the Keelung River. The model simulations indicate that more tidal energy propagates into the estuarine system before channel regulation because of the substantial increase in river cross‐sections. The residual circulations before channel regulation are greater than those after channel regulation and result in the limits of the salt intrusion before channel regulation being extended farther inland than those after channel regulation. This may show that channel regulation for flood control in the Keelung River did not contribute to the expansion of the mangrove areas and the disappearance of freshwater marshes at the Kuan‐Du wetlands. Copyright © 2005 John Wiley & Sons, Ltd.     

Horizontal patterns of water temperature and salinity in an estuarine tidal channel: Ria de Aveiro

This work presents results from two complementary and interconnected approaches to study water temperature and salinity patterns in an estuarine tidal channel. This channel is one of the four main branches of the Ria de Aveiro, a shallow lagoon located in the Northwest coast of the Iberian Peninsula. Longitudinal and cross-sectional fields of water temperature and salinity were determined by spatial interpolation of field measurements. A numerical model (Mohid) was used in a 2D depth-integrated mode in order to compute water temperature and salinity patterns. The main purpose of this work was to determine the horizontal patterns of water temperature and salinity in the study area, evaluating the effects of the main forcing factors. The field results were depth-integrated and compared to numerical model results. These results obtained using extreme tidal and river runoff forcing, are also presented. The field results reveal that, when the river flow is weak, the tidal intrusion is the main forcing mechanism, generating saline and thermal fronts which migrate with the neap/spring tidal cycle. When the river flow increases, the influence of the freshwater extends almost as far as the mouth of the lagoon and vertical stratification is established. Results of numerical modelling reveal that the implemented model reproduces quite well the observed horizontal patterns. The model was also used to study the hydrology of the study area under extreme forcing conditions. When the model is forced with a low river flow (1 m³ s⁻¹) the results confirm that the hydrology is tidally dominated. When the model is forced with a high river flow (1,000 m³ s⁻¹) the hydrology is dominated by freshwater, as would be expected in such an area.     

The response of circulation and salinity in a micro-tidal estuary to sub-tidal oscillations in coastal sea surface elevation

Conceptual models of circulation theorise that the dominant forces controlling estuarine circulation are freshwater discharge from the riverine section (landward), tidal forcing from the ocean boundary, and gravitational circulation resulting from along-estuary gradients in density. In micro-tidal estuaries, sub-tidal water level changes (classified as those with periods between 3 and 10 days) with amplitudes comparable to the spring tidal range can significantly influence the circulation and distribution of water properties. Field measurements obtained from the Swan River Estuary, a diurnal, micro-tidal estuary in south-western Australia, indicated that sub-tidal water level changes at the ocean boundary were predominantly from remotely forced continental shelf waves (CSWs). The sub-tidal water levels had maximum amplitudes of 0.8 m, were comparable to the maximum tidal range of 0.6 m, propagated into the estuary to its tidal limit, and modified water levels in the whole estuary over several days. These oscillations dominated the circulation and distribution of water properties in the estuary through changing the salt wedge location and increasing the bottom water salinity by 7 units over 3 days. The observed salt wedge excursion forced by CSW was up to 5 km, whereas the maximum tidal excursion was 1.2 km. The response of the residual currents and the salinity distribution lagged behind the water level changes by ∼24 h. It was proposed that the sub-tidal forcing at the ocean boundary, which changed the circulation, salinity, and dissolved oxygen in the upper estuary, was due to a combination of two processes: (1) a gravity current generated by a process similar to a lock exchange mechanism and (2) amplified along-estuary density gradients in the upper estuary, which enhanced the gravitational circulation in the estuary. The salt intrusions under the sub-tidal forcing caused the rapid movement of anoxic water upstream, with significant implications for water quality and estuarine health.     

Atoll island hydrogeology: flow and freshwater occurrence in a tidally dominated system

A layered-aquifer model of groundwater occurrence in an atoll island was tested with a solute-transport numerical model. The computer model used, SUTRA, incorporates density-dependent flow. This can be significant in freshwater-saltwater interactions associated with the freshwater lens of an atoll island. Boundary conditions for the model included ocean and lagoon tidal variations. The model was calibrated to field data from Enjebi Island, Enewetak Atoll, and tested for sensitivity to a variety of parameters. This resulted in a hydraulic conductivity of 10 m day⁻¹ for the surficial aquifer and 1000 m day⁻¹ for the deeper aquifer; this combination of values gave an excellent reproduction of the tidal response data from test wells. The average salinity distribution was closely reproduced using a dispersivity of 0.02m. The computer simulation quantitatively supports the layered-aquifer model, including under conditions of density-dependent flow, and shows that tidal variations are the predominant driving force for flow beneath the island. The oscillating, vertical flow produced by the tidal variations creates an extensive mixing zone of brackish water. The layered-aquifer model with tidally driven flow is a significant improvement over the Ghyben-Herzberg-Dupuit model as it is conventionally applied to groundwater studies for many Pacific reef islands.     

A modeling study of the physical processes affecting the development of seasonal hypoxia over the inner Louisiana-Texas shelf: Circulation and stratification

The physical processes affecting the development of seasonal hypoxia over the Louisiana-Texas shelf were examined using a high-resolution, three-dimensional, unstructured-grid, Finite Volume Coastal Ocean Model (FVCOM). The model was forced with the observed freshwater fluxes from the Mississippi and Atchafalaya Rivers, surface winds, heat fluxes, tides and offshore conditions. The simulations were carried out over a six-month period, from April to September 2002, and the model performance was evaluated against several independent series of observations that included tidal gauge data, Acoustic Doppler Current Profiler (ADCP) data, shipboard measurements of temperature and salinity, vertical salinity and sigma-t profiles, and satellite imagery. The model accurately described the offshore circulation mode generated over the Louisiana-Texas shelf by the westerly winds during summer months, as well as the prevalent westward flow along the coast caused by the easterly winds during the rest of the study period. The seasonal cycle of stratification also was well represented by the model. During 2002, the stratification was initiated in early spring and subsequently enhanced by the intensity and phasing of riverine freshwater discharges. Strong stratification persisted throughout the summer and was finally broken down in September by tropical storms. The model simulations also revealed a quasi-permanent anticyclonic gyre in the Louisiana Bight region formed by the rotational transformation of the Mississippi River plume, whose existence during 2002 was supported by the satellite imagery and ADCP current measurements. Model simulations support the conclusion that local wind forcing and buoyancy flux resulting from riverine freshwater discharges were the dominant mechanisms affecting the circulation and stratification over the inner Louisiana-Texas shelf.     

Coupled 1D–3D hydrodynamic modelling,with application to the Pearl River Delta

Within the hydrodynamic modelling community, it is common practice to apply different modelling systems for coastal waters and river systems. Whereas for coastal waters 3D finite difference or finite element grids are commonly used, river systems are generally modelled using 1D networks. Each of these systems is tailored towards specific applications. Three-dimensional coastal water models are designed to model the horizontal and vertical variability in coastal waters and are less well suited for representing the complex geometry and cross-sectional areas of river networks. On the other hand, 1D river network models are designed to accurately represent complex river network geometries and complex structures like weirs, barrages and dams. A disadvantage, however, is that they are unable to resolve complex spatial flow variability. In real life, however, coastal oceans and rivers interact. In deltaic estuaries, both tidal intrusion of seawater into the upstream river network and river discharge into open waters play a role. This is frequently approached by modelling the systems independently, with off-line coupling of the lateral boundary forcing. This implies that the river and the coastal model run sequentially, providing lateral discharge (1D) and water level (3D) forcing to each other without the possibility of direct feedback or interaction between these processes. An additional disadvantage is that due to the time aggregation usually applied to exchanged quantities, mass conservation is difficult to ensure. In this paper, we propose an approach that couples a 3D hydrodynamic modelling system for coastal waters (Delft3D) with a 1D modelling system for river hydraulics (SOBEK) online. This implies that contrary to off-line coupling, the hydrodynamic quantities are exchanged between the 1D and 3D domains during runtime to resolve the real-time exchange and interaction between the coastal waters and river network. This allows for accurate and mass conserving modelling of complex coastal waters and river network systems, whilst the advantages of both systems are maintained and used in an optimal and computationally efficient way. The coupled 1D–3D system is used to model the flows in the Pearl River Delta (Guangdong, China), which are determined by the interaction of the upstream network of the Pearl River and the open waters of the South China Sea. The highly complex upstream river network is modelled in 1D, simulating river discharges for the dry and wet monsoon periods. The 3D coastal model simulates the flow due to the external (ocean) periodic tidal forcing, the salinity distribution for both dry and wet seasons, as well as residual water levels (sea level anomalies) originating from the South China Sea. The model is calibrated and its performance extensively assessed against field measurements, resulting in a mean root mean square (RMS) error of below 6% for water levels over the entire Pearl River Delta. The model also represents both the discharge distribution over the river network and salinity transport processes with good accuracy, resolving the discharge distribution over the main branches of the river network within 5% of reported annual mean values and RMS errors for salinity in the range of 2 ppt (dry season) to 5 ppt (wet season).     

Modelling of hydrodynamics and cohesive sediment transport in Tanshui River estuarine system,Taiwan

A laterally averaged two-dimensional numerical model is used to simulate hydrodynamics and cohesive sediment transport in the Tanshui River estuarine system. The model handles tributaries as well as the main stem of the estuarine system. Observed time series of salinity data and tidally averaged salinity distributions have been compared with model results to calibrate the turbulent diffusion coefficients. The overall model verification is achieved with comparisons of residual currents and salinity distribution. The model reproduces the prototype water surface elevation, currents and salinity distributions. Comparisons of the suspended cohesive sediment concentrations calculated by the numerical model and the field data at various stations show good agreement. The validated model is applied to investigate the tidally averaged salinity distributions, residual circulation and suspended sediment concentration under low flow conditions in the Tanshui River estuarine system. The model results show that the limit of salt intrusion in the mainstem estuary is located at Hsin-Hai bridge in Tahan Stream, 26 km from the River mouth under Q75 flow. The null point is located at the head of salt intrusion, using 1 ppt isohaline as an indicator. The tidally averaged sediment concentration distribution exhibits a local maximum around the null point.     

Numerical modelling of morphodynamics—Vilaine Estuary

The main objective of this paper is to develop a method to simulate long-term morphodynamics of estuaries dominated by fine sediments, which are subject to both tidal flow and meteorologically induced variations in freshwater run-off and wave conditions. The method is tested on the Vilaine Estuary located in South Brittany, France. The estuary is subject to a meso–macrotidal regime. The semi-diurnal tidal range varies from around 2.5 to 5 m at neap and spring, respectively. The freshwater input is controlled by a dam located approximately 8 km from the mouth of the estuary. Sediments are characterised as mostly fines, but more sandy areas are also found. The morphology of the estuary is highly influenced by the dam. It is very dynamic and changes in a complicated manner with the run-off from the dam, the tide and the wave forcing at the mouth of the estuary. Extensive hydrodynamic and sediment field data have been collected in the past and provide a solid scientific basis for studying the estuary. Based on a conceptual understanding of the morphodynamics, a numerical morphological model with coupled hydrodynamic, surface wave and sediment transport models is formulated. The numerical models are calibrated to reproduce sediment concentrations, tidal flat altimetry and overall sediment fluxes. Scaling factors are applied to a reference year to form quasi-realistic hydrodynamic forcing and river run-off, which allow for the simulations to be extended to other years. The simulation results are compared with observed bathymetric changes in the estuary during the period 1998–2005. The models and scaling factors are applied to predict the morphological development over a time scale of up to 10 years. The influence of the initial conditions and the sequence of external hydrodynamic forcing, with respect to the morphodynamic response of the estuary, are discussed.     

Long‐term monitoring of a mercury contaminated estuary (Ria de Aveiro,Portugal): the effect of weather events and management in mercury transport

The main aim of this research was to assess the mercury transport from an estuarine basin with a background of anthropogenic contamination during a spring tidal cycle (year 2009) and compare it with two previous tidal cycles (years 1994 and 1999), as part of a long‐term monitoring program. Results showed that effective mercury transport occurs both in the dissolved and particulate fractions (0.18 and 0.20 kg per tidal cycle, respectively), and despite an overall decrease in environmental contamination, results more than double previous findings on particulate transport in the system. These findings result essentially from changes in the tidal prism (net export of 2 million m³ of water), given that both dissolved and particulate concentrations did not increase over time. Hydrodynamic simulations were performed to evaluate the effect of physical disturbance (dredging) and weather events (increased freshwater flow) in these processes, and results suggest the increased freshwater flow into the system as the main forcing function for the mercury transport increment. These results highlight the importance of long‐term monitoring programs, since despite an overall improvement in local contamination levels, the enhancement of transport processes through hydrological changes increases environmental pressure away from the contamination source. Copyright © 2012 John Wiley & Sons, Ltd.     

Fate of three major rivers in the Bohai Sea: A model study

Huanghe (Yellow River), Haihe and Liaohe are three major rivers flowing into the Bohai Sea and account for more than 80% of the freshwater and land-drained material inputs annually. The fate of three rivers in the seawaters correlates with the transport and distribution of the riverine sediments and nutrients, and further exerts a profound influence on the local marine ecosystem dynamics. Therefore, the evolution of the river plumes under the influence of the freshwater buoyancy, the tidal forcing and the wind stress are examined using a three-dimensional primitive equation ocean circulation model, independently and jointly. It is found that both tide and wind stirring can deteriorate the stabilization of the water column caused by the freshwater buoyancy; however, the processes are different. The tide stirring originates from the seafloor due to the bottom friction as the tidal wave propagates into the shallow waters, and then the turbulent kinetic energy dissipates upward. On the other hand, the wind stirring proceeds in the up-down direction. The influences of different winds on the evolution of the river plumes are also examined. Since the situation of each river mouth is different, the wind influence is also distinct. At last, the fate of three major rivers driven by the combined tidal forcing and climatology winds is reproduced, and the simulated salinity distribution shows a reasonable agreement with that observed, meaning that the river plume evolution plays a crucial role in shaping the salinity distribution in BS.     

Variations of tidally driven three-layer residual circulation in fjords

Residual, or tidally averaged, circulation in fjords is generally assumed to be density driven and two layered. This circulation consists of a thin surface layer of outflow and a thick bottom layer of sluggish inflow. However, development of different vertical structures in residual circulation in fjords can arise from wind, remote, and tidal forcing that may modify the two-layer circulation. Particularly, theoretical results of tidal residual flows in homogeneous semienclosed basins indicate that their vertical structure is determined by the dynamical depth of the system. This dynamical depth can be considered as the ratio between the water column depth and the depth of frictional influence in an oscillatory flow (inverse of Stokes number). When the frictional depth occupies the entire water column, the tidal residual flow is one layered as in shallow basins. But when the frictional depth is only a small portion of the water column (>6 times smaller), the tidal residual is three layered. In relatively deep fjords (say deeper than 100 m), where frictional depths typically occupy a small portion of the water column, the tidal residual flow is expected to be three layered. Ample observational evidence presented here shows a three-layered exchange flow structure in fjords. On the basis of observational and theoretical evidence, it is proposed that the water exchange structure in deep fjords (more than six frictional layers deep, or inverse Stokes number >6) is tidally driven and is three layered. The tidally driven three-layer structure of residual flows could be regarded in some cases as the fundamental structure. However, this structure will only be observed sporadically as it will be masked by wind forcing, remote forcing from the ocean, and freshwater pulses.     

A model study of tide- and wind-induced mixing in the Columbia River Estuary and plume

A numerical simulation of circulation in the Columbia River estuary and plume during the summer of 2004 is used to explore the mixing involved as river water is transformed into shelf water. The model is forced with realistic river flow, tides, wind stress, surface heat flux, and ocean boundary conditions. Simulated currents and water properties on the shelf near the mouth are compared with records from three moorings (all in 72 m of water) and five CTD sections. The model is found to have reasonable skill; statistically significant correlations between observed and modeled surface currents, temperature, and salinity are all 0.42–0.72 for the mooring records. Equations for the tidally averaged, volume-integrated mechanical energy budget (kinetic and potential) are derived, with attention to the effects of: (i) Reynolds averaging, (ii) a time varying volume due to the free surface, and (iii) dissipation very close to the bottom. It is found that convergence of tidal pressure work is the most important forcing term in the estuary. In the far field plume (which has a volume 15 times greater than that of the estuary), the net forcing is weaker than that in the estuary, and may be due to either tidal currents or wind stress depending on the time period considered. These forcings lead to irreversible mixing of the stratification (buoyancy flux) that turns river water into shelf water. This occurs in both the plume and estuary, but appears to be more efficient (17% vs. 5%), and somewhat greater (4.2 MW vs. 3.3 MW), in plume vs. estuary. This demonstrates the importance of both wind and tidal forcing to watermass transformation, and the need to consider the estuary and plume as part of a single system.     

Suspended sediment fluxes at an intertidal flat: The shifting influence of wave,wind, tidal,and freshwater forcing

Using in situ, continuous, high frequency (8–16 Hz) measurements of velocity, suspended sediment concentration (SSC), and salinity, we investigate the factors affecting near-bed sediment flux during and after a meteorological event (cold front) on an intertidal flat in central San Francisco Bay. Hydrodynamic forcing occurs over many frequency bands including wind wave, ocean swell, seiching (500–1000 s), tidal, and infra-tidal frequencies, and varies greatly over the time scale of hours and days. Sediment fluxes occur primarily due to variations in flow and SSC at three different scales: residual (tidally averaged), tidal, and seiching. During the meteorological event, sediment fluxes are dominated by increases in tidally averaged SSC and flow. Runoff and wind-induced circulation contribute to an order of magnitude increase in tidally averaged offshore flow, while waves and seiching motions from wind forcing cause an order of magnitude increase in tidally averaged SSC. Sediment fluxes during calm periods are dominated by asymmetries in SSC over a tidal cycle. Freshwater forcing produces sharp salinity fronts which trap sediment and sweep by the sensors over short (∼30 min) time scales, and occur primarily during the flood. The resulting flood dominance in SSC is magnified or reversed by variations in wind forcing between the flood and ebb. Long-term records show that more than half of wind events (sustained speeds of greater than 5 m/s) occur for 3 h or less, suggesting that asymmetric wind forcing over a tidal cycle commonly occurs. Seiching associated with wind and its variation produces onshore sediment transport. Overall, the changing hydrodynamic and meteorological forcing influence sediment flux at both short (minutes) and long (days) time scales.     

A simple approach to estimating freshwater discharge in branched estuarine systems

In tidal estuaries, quantifying freshwater discharge is still a difficult problem that has not yet been overcome due to the inherent difficulty in measuring and analysing the tidal discharge, especially during periods of low river flow. Because observations are often made in the stations further upstream, where the ratio of river to tidal discharge is large, it remains difficult to determine the discharge rate in the saline region. Freshwater discharge estimation is even more difficult in a branched estuary system having multiple diversion channels that connect with each other at a junction. To date, several methods have been developed for estimating freshwater discharge in estuaries. The most widely used are analytical and conceptual models that employ salinity as the principal trace and numerical simulations. However, these methods are very time consuming and costly as they require large sets of observations before the computations can take place. This paper presents a simple approach to investigating the discharge distribution over branched channels by considering the energy loss due to friction. We develop an analytical model that can obtain the discharge rate quantitatively at a junction where the main flow bifurcates into two branches. The model uses the bed roughness, tidal water level, and cross‐sectional profile under tidally averaged conditions as input data. Two selected estuarine systems in the Hiroshima delta in Japan and the Mekong delta in Vietnam have been investigated. Computations of the newly developed model show good agreement with earlier published results computed by sophisticated analytical and numerical models.     

Investigation of hydraulics transport time scales within the Arvand River estuary,Iran

Transport time scales are key parameters for understanding the hydrodynamic and biochemical processes within estuaries. In this study, the flushing and residence times within the Arvand River estuary have been estimated using a two‐dimensional hydrodynamic model called CE‐QUAL‐W2. The model has been calibrated and verified by two different sets of field data and using the k‐ε vertical eddy diffusivity scheme. Flushing time has been estimated using different methods such as the tidal prism and fraction of freshwater methods. Moreover, residence times have been investigated using pulse residence time, estuarine residence time and remnant function approaches. The results have shown that different methods yield different time scales, and freshwater inflow has the greatest impact upon estimation of residence time, whereas tidal circulation hardly contributes to residence time at all. It has also been shown that the neap‐spring circulation and start phase of simulations have negligible effects on the Arvand's time scales. The investigation of bathymetry showed that two sills of the estuary tend to significantly increase residence time. Understanding the applicability of these time scales and their estimation approaches helps us to evaluate the water quality management of estuaries. Copyright © 2013 John Wiley & Sons, Ltd.     

The tidal and wind induced hydrodynamics of the composite system Adriatic Sea/Lagoon of Venice

A shallow water hydrostatic 2D hydrodynamic numerical model, based on the boundary conforming coordinate system, was used to simulate aspects of both general and small scale oceanic features occurring in the composite system constituted by the Adriatic Sea and the Lagoon of Venice (Italy), under the influence of tide and realistic atmospheric forcing. Due to a specific technique for the treatment of movable lateral boundaries, the model is able to simulate efficiently dry up and flooding processes within the lagoon. Firstly, a model calibration was performed by comparing the results of the model, forced using tides and ECMWF atmospheric pressure and wind fields, with observations collected for a set of 33 mareographic stations uniformly distributed in the Adriatic Sea and in the Lagoon of Venice. A second numerical experiment was then carried out by considering only the tidal forcing. Through a comparison between the results obtained in the two experiments it was possible to assess the reliability of the estimated parameter through the composite forcing. Model results were then verified by comparing simulated amplitude and phase of each tidal constituent as well as tidal velocities simulated at the inlets of the lagoon and in the Northern Adriatic Sea with the corresponding observed values. The model accurately reproduces the observed harmonics: mean amplitude differences rarely exceed 1 cm, while phase errors are commonly confined below 15°. Semidiurnal and diurnal currents were correctly reproduced in the northern basin and a good agreement was obtained with measurements carried out at the lagoon inlets. On this basis, the outcomes of the hydrodynamic model were analyzed in order to investigate: (i) small-scale coastal circulation features observed at the interface between the adjoining basins, which consist often of vortical dipoles connected with the tidal flow of Adriatic water entering and leaving the Lagoon of Venice and with along-shore current fields connected with specific wind patterns; (ii) residual oscillations, which are often connected to meteorological forcing over the basin. In particular, it emerges that small-scale vortical features generated near the lagoon inlet can be efficiently transported toward the open sea, thus contributing to the water exchange between the two marine regions, and a realistic representation of observed residual oscillations in the area would require a very detailed knowledge of atmospheric as well as remote oceanic forcing.     

Coupled flow and salinity transport modelling in semi-arid environments: The Shashe River Valley,Botswana

Numerical groundwater modelling is used as the base for sound aquifer system analysis and water resources assessment. In many cases, particularly in semi-arid and arid regions, groundwater flow is intricately linked to salinity transport. A case in point is the Shashe River Valley in Botswana. A freshwater aquifer located around an ephemeral stream is depleted by the combined effect of transpiration and pumping. Quantitative system analysis reveals that the amount of water taken by transpiration is far more than the quantities pumped for water supply. Furthermore, the salinity distribution in and around Shashe River Valley as well as its temporal dynamics can be satisfactorily reproduced if the transpiration is modelled as a function of groundwater salinity. The location and dynamics of the saltwater–freshwater interface are highly sensitive to the parameterization of evaporative and transpirative salt enrichment. An existing numerical code for coupled flow/transport simulations (SEAWAT) was adapted to this situation. Model results were checked against a large set of field data including water levels, water chemistry, isotope data and ground and airborne geophysical data. The resulting groundwater model was able to reproduce the long-term development of the freshwater lens located in Shashe River Valley as well as the decline in piezometric heads observed over the last decade. Furthermore, the old age of the saline water surrounding the central freshwater lens could be explained.     

A 3-D numerical study of salinity variations in the Bohai Sea during the recent years

Salinity is an important component of the marine system. Previous studies indicated that the mean salinity in the Bohai Sea had increased by 2.0 psu in the second half of the 20th century, mainly due to a sharp decrease in the Yellow River runoff, and also the effects of large-scale climatic variations and the intrusions of the North Yellow Sea Water (NYSW). Since 2002, the Yellow River Conservancy Commission has carried out the flow regulation at the beginning of every flood season, resulting in more discharge of the Yellow River freshwater into the Bohai Sea. In this study, the variations of salinity in the Bohai Sea during the recent years are investigated using a well-established three-dimensional baroclinic model, HAMburg Shelf Ocean Model (HAMSOM). The simulation results show that the Yellow River diluted water was mainly discharged into the Laizhou Bay, so the remarkable increase in the Yellow River runoff after 2002 led to a regime shift of salinity in the Laizhou Bay. However, in other parts of the Bohai Sea, salinity variation was influenced by the surrounding rivers or the intrusions of NYSW, and has little relation with the Yellow River runoff. As a whole, advection is more important than diffusion in the salinity distribution, and seasonal oscillation is the main feature of salinity variation. Via several case studies, evaporation and precipitation rates are found to be important in the long-term simulation of salinity.