Estimation of net physical transport and hydraulic residence times for a coastal plain estuary using box models

A box model based on salinity distributions and freshwater inflow measurements was developed and used to estimate net non-tidal physical circulation and hydraulic residence times for Patuxent River estuary, Maryland, a tributary estuary of Chesapeake Bay. The box model relaxes the usual assumption that salinity is at steady-state, an important improvement over previous box model studies, yet it remains simple enough to have broad appeal. Average monthly 2-dimensional net non-tidal circulation and residence times for 1986–1995 are estimated and related to river flow and salt water inflow as estimated by the box model. An important result is that advective exchange at the estuary mouth was not correlated with Patuxent River flow, most likely due to effects of offshore salinity changes in Chesapeake Bay. The median residence time for freshwater entering at the head of the estuary was 68 d and decreased hyperbolically with increasing river flow to 30 d during high flow. Estimates of residence times for down-estuary points of origin showed that, from the head of the estuary to its mouth, control of flushing changed from primarily river flow to other factors regulating the intensity of gravitational circulation.     

Impact of cloud microphysics and cumulus parameterization on simulation of heavy rainfall event during 7–9 October 2007 over Bangladesh

In the present study, the Advanced Research WRF (ARW) version 3.2.1 has been used to simulate the heavy rainfall event that occurred between 7 and 9 October 2007 in the southern part of Bangladesh. Weather Research and Forecast (WRF–ARW version) modelling system with six different microphysics (MP) schemes and two different cumulus parameterization (CP) schemes in a nested configuration was chosen for simulating the event. The model domains consist of outer and inner domains having 9 and 3 km horizontal resolution, respectively with 28 vertical sigma levels. The impacts of cloud microphysical processes by means of precipitation, wind and reflectivity, kinematic and thermodynamic characteristics of the event have been studied. Sensitivity experiments have been conducted with the WRF model to test the impact of microphysical and cumulus parameterization schemes in capturing the extreme weather event. NCEP FNL data were used for the initial and boundary condition. The model ran for 72 h using initial data at 0000 UTC of 7 October 2007. The simulated rainfall shows that WSM6–KF combination gives better results for all combinations and after that Lin–KF combination. WSM3–KF has simulated, less area average rainfall out of all MP schemes that were coupled with KF scheme. The sharp peak of relative humidity up to 300 hPa has been simulated along the vertical line where maximum updraft has been found for all MPs coupled with KF and BMJ schemes. The simulated rain water and cloud water mixing ratio were maximum at the position where the vertical velocity and reflectivity has also been maximum. The production of rain water mixing ratio depends on MP schemes as well as CP schemes. Rainfall depends on rain water mixing ratio between 950 and 500 hPa. Rain water mixing ratio above 500 hPa level has no effect on surface rain.     

Impact of Channel Deepening on Tidal and Gravitational Circulation in a Highly Engineered Estuarine Basin

Deepening of estuarine channels is a common practice to ensure navigation. Here, we investigate whether such deepening impacts physical processes such as the strength of the estuarine exchange flow, the horizontal salinity gradient, and tidal dynamics. We analyze recent and historical hydrodynamic observations in Newark Bay, New Jersey, to assess the effect of channel deepening on tides, circulation, and salinity. The Bay’s navigational channel has undergone significant deepening, from 3 to 10 m in the nineteenth century to ~16 m today. Observations presented here include sea-level data from the nineteenth, twentieth, and twenty-first century, and moored Doppler current data and bottom salinity measurements made over the past 20 years. Results show a doubling of the estuarine exchange flow, a slight increase in salinity and in the horizontal salinity gradient, a decrease in tidal current amplitude, and a spatially variable change in the tidal range. The doubling of the exchange flow is consistent with the Hansen and Rattray scaling provided that the horizontal salinity gradient is unable to fully adjust landward because the dredging is limited to a short reach of the estuary. However, uncertainty in channel depth leaves open the possibility that the exchange flow is also augmented by an increase in the horizontal salinity gradient and/or a reduction in vertical mixing. Nevertheless, results demonstrate that a relatively small (15%) increase in depth appears to have doubled the exchange flow. We believe that this result is relevant to other systems where dredging is limited to a short reach of an estuary.     

Observations of currents and density structure across a buoyant plume front

Observations of the Mobile Bay, Alabama, plume during a flood event in April 1991 reveal significant differences in the current field on either side of a front associated with the buoyant plume. During a strong southeasterly wind, turbid, low salinity water from Mobile Bay was pushed through an opening in the west side of the ebb-tidal delta and moved parallel to the coast. A stable front developed between the low salinity water of the buoyant plume (11‰) and the high salinity coastal water (>23‰) that was being forced landward by the prevailing winds. Despite the shallow water depth of 6 m, measurements of currents, temperature, and salinity show large shears and density gradients in both the vertical and the horizontal directions. At a station outside of the buoyant plume, currents at 0.5 m and 1.5 m below the surface were in the same direction as the wind. Inside the plume, however, currents at 0.5 m below the surface were parallel to the coast, 45°, off the direction of the wind and the magnitude was 45% larger than the magnitude of the surface currents outside the plume. Beneath the level of the plume, the currents were identical to the wind-driven currents in the ambient water south of the front. Our observations suggest that the wind-driven surface currents of the ambient water converged with the buoyant plume at the front and were subducted beneath the plume. The motion of the ambient coastal surface water was in the direction of the local wind stress, however, the motion of the plume had no northerly component of motion. The plume also did not show any flow toward the front, suggesting a balance between the northerly component of wind stress and the southerly component of buoyant spreading. In addition, the motion of the plume did not appear to affect the motion of the underlying ambient water, suggesting a lack of mixing between the two waters.     

Time-dependent modeling of the Baltic entrance area. 1. Quantification of circulation and residence times in the Kattegat and the Straits of the Baltic sill

A prognostic model for the estuarine circulation in the Baltic entrance area is described. The model is based on the work by Stigebrandt (1983) and is built of a number of flow-regulating physical processes and forced by oceanic sea level, local meteorological conditions, and freshwater supply to the Baltic. It resolves the horizontal variations by dividing the entrance area into 7 sub-basins. The vertical stratification in each sub-basin evolves with time due to mixing, diffusion, and water exchange with adjacent basins. Instead of using a fixed vertical grid, the stratification is described by a variable number of layers created by inflowing water and by a pycnocline retreat near the surface. The model is a highly cost-effective tool compared to high-resolution 3D-models since computations are 10⁵–10⁶ times faster. Still, the model reproduces the stratification on time-scales longer than a few days as verified by comparison with observed time-series. The magnitude of the simulated average estuarine circulation conforms well to independent estimates of the water exchange. The model is used to quantify the temporal and horizontal variability of circulation, mixing, and diffusion. Long-term average rate of work against the buoyancy forces by entrainment and diffusion is calculated for each sub-basin. The highest rates of work against the buoyancy forces by diffusion are found in the northern Kattegat and in the Fehmarn Belt while the lowest rate is found in the Öresund. The total mixing in terms of transformation of water masses is also quantified. The average residence times for surface and deep water are estimated for each subbasin. Residence times longer than 1 mo are found in Fehmarn Belt and in the deep water of southern Kattegat. In other parts of Kattegat and the Samsø Belt the residence times are 1 to 2 wk for surface water and 2 to 3 wk for deep water. The age of the water is defined as the time spent since a particular water mass was in contact with either the sea surface or with the vertical model boundaries has been estimated. By comparing the age distribution in the Kattegat with observations of oxygen concentration, it is concluded that the variability of the ventilation of deeper layers is of less importance to the occurrence of low oxygen concentrations compared to other factors such as interannual variability in primary production.     

Exchange and pathways of deep and shallow groundwater in different climate and permafrost conditions using the Forsmark site, Sweden, as an example catchment

This study simulates and quantifies the exchange and the pathways of deep and shallow groundwater flow and solute transport under different climate and permafrost conditions, considering the example field case of the coastal Forsmark catchment in Sweden. A number of simulation scenarios for different climate and permafrost condition combinations have been performed with the three-dimensional groundwater flow and transport model MIKE SHE. Results show generally decreasing vertical groundwater flow with depth, and smaller vertical flow under permafrost conditions than under unfrozen conditions. Also the overall pattern of both the vertical and the horizontal groundwater flow, and the water exchange between the deep and shallow groundwater systems, change dramatically in the presence of permafrost relative to unfrozen conditions. However, although the vertical groundwater flow decreases significantly in the presence of permafrost, there is still an exchange of water between the unfrozen groundwater system below the permafrost and the shallow groundwater in the active layer, via taliks. ‘Through taliks’ tend to prevail in areas that constitute groundwater discharge zones under unfrozen conditions, which then mostly shift to net recharge zones (through taliks with net downward flow) under permafrost conditions.     

Spatial and Seasonal Patterns of Salinity in a Large and Shallow Tropical Estuary of the Western Caribbean

Salinity profiles and meteorological data were analyzed during February, May, and September 2006 in Chetumal Bay, a large, shallow estuary of the Western Caribbean. Local meteorological conditions revealed three seasons: (1) a dry season (March–May); (2) a wet season (June–October); and (3) the nortes season, with northerly wind events (October–February). During the nortes and wet seasons, salinity ranged between 13 and 16 psu, and salinity was highest in the dry season, ranging between 18 and 22 psu over most of the area; a strong stratification and a significant contribution of salty water characterized this season. Strong horizontal gradients were observed near Rio Hondo during the three seasons. Deep and narrow peculiar bathymetric features called the pozas showed a strong stratification and a relatively high salinity. The northern part of Chetumal Bay and probably the entire system are far from being homogeneous.     

Examination of Winter Circulation in a Northern Gulf of Mexico Estuary

Numerical model experiments were conducted to examine how estuarine circulation and salinity distribution in the Calcasieu Lake Estuary (CLE) of southwest Louisiana respond to the passage of cold fronts. River runoff, local wind stress, and tides from December 20, 2011, to February 1, 2012, were included as input. The experiments showed an anticyclonic circulation in the eastern CLE, a cyclonic circulation in West Cove, and a saltwater conduit in the navigation channel between these circulation cells. Freshwater from the river and wetlands tends to flow over the shallow shoals toward the ocean, presenting a case of the conventional estuarine circulation with shallow water influenced by river discharge and with weak tidally-induced motion, enhanced by wind. The baroclinic pressure gradient is important for the circulation and saltwater intrusion. The effect of remote wind-driven oscillation plays an important role in circulation and salinity distribution in winter. Unless it is from the east, wind is found to inhibit saltwater intrusion through the narrow navigation channel, indicating the effect of Ekman setup during easterly wind. A series of north-south oriented barrier islands in the lake uniquely influenced water level and salinity distribution between the shallow lake and deep navigation channel. The depth of the navigation channel is also crucial in influencing saltwater intrusion: the deeper the channel, the more saltwater intrusion and the more intense estuarine circulation. Recurring winter storms have a significant accumulated effect on the transport of water and sediment, saltwater intrusion, and associated environmental and ecosystem effects.     

Geostatistics as a groundwater exploration planning tool: case of a brackish-saline aquifer

This paper discusses a geostatistical approach to model a groundwater aquifer in 3-D. The study aims at utilizing geostatistics as a tool for characterizing zones of better-water quality in a brackish-saline aquifer. In particular, the geostatistical model was constructed to characterize the aquifer’s salinity, represented by total dissolved solids (TDS), using logs of porosity and resistivity. Quality-checked estimated TDS vertical profiles were employed to construct and model horizontal and vertical semivariograms. Parameters of semivariogram models were used to develop both the kriging plan and the generated model. Results of this modeling process are shown in the form of horizontal salinity distribution maps. The aquifer was sliced into 20 layers, each 20 m thick, to represent its overall thickness. Salinity layers maps reflect vertical stratification of TDS concentrations in the aquifer and show that water quality deteriorates with depth and toward the northern part of the aquifer. Relatively better-quality water (TDS ≤10,000 mg/l) can be found at depths between 100 and 250 m below the aquifer’s top in both eastern and southeastern parts. Water in the same interval to the western and southwestern parts reflects the presence of higher TDS concentration. From a planning point of view, it is more feasible to target the eastern part of the aquifer for pumping and desalination purposes. In addition, the generated model could be utilized as an initial condition for flow simulation.     

沟灌土壤湿润体空间矩特征参数估算模型

为进一步增加矩分析理论在研究沟灌土壤湿润体运移特性方面的实用性,采用数值模拟和理论分析相结合的方法,定量分析了土壤初始含水率、入渗水深、沟底宽和边坡系数等因素对沟灌自由入渗条件下土壤湿润体空间矩特征参数的影响,并建立了其与主要影响因素间的函数关系式。结果表明,土壤初始含水率、入渗水深和沟底宽对垂直向土壤水分分布的质量重心影响较大;水平向湿润锋的平均离散度主要受土壤初始含水率、入渗水深、沟底宽和边坡系数的影响,而垂直向对土壤初始含水率和入渗水深变化的敏感性较高;选取不同土壤质地对所建模型的可靠性进行了验证,结果表明所建模型估算不同入渗时刻的空间矩特征参数值与根据HYDRUS软件模拟结果所得计算值具有高度一致性,所有验证组合条件下两者相对误差绝对值均值均低于8.5%,且无显著性差异,说明所建模型估算沟灌土壤湿润体空间矩特征参数具有高度可靠性。研究结果可为沟灌灌水方案设计和管理提供理论基础和技术支撑。     

Hydrologic and geochemical survey of the lake “Specchio di Venere” (Pantelleria island,Southern Italy)

Hydrological and geochemical studies performed on Lake Specchio di Venere on Pantelleria island (Italy) indicate that this endorheic basin has been formed through upwelling of the water table, and that it is continuously fed by the thermal springs situated on its shores. The lake is periodically stratified both thermally and in salinity, albeit this stratification is rather unstable over time, since meteorological events such as strong rain or wind can determine the mixing of its waters. Periodical analyses of the lake water chemistry show large variations of the salt content due to the yearly evaporation-rain dilution cycle. These processes are also responsible for the saline stratification during steady meteorological conditions. The mineralogical characterisation of the bottom sediments shows the almost exclusive presence of neoformation minerals, mainly carbonates, formed in response to the pH gradient between spring- (pH ≈ 6) and lake-waters (pH ≈ 9). Finally, the CO₂ partial pressures in the lake water slightly exceeding the atmospheric one, are due to the large amounts of CO₂ brought to the lake through the bubbling free gas phase of the thermal springs. Nevertheless the high pH value of the lake water, its small volume and its periodical mixing prevent dangerous built up of this gas.     

Circulation in Carr Inlet,Puget Sound,during Spring 2003

The relatively slow flow and exchange of Carr Inlet water with the main basin of Puget Sound, Washington, favor eutrophication. To study Carr Inlet’s circulation, the Model-measurement Integration Experiment in Estuary Dynamics (MIXED) was conducted in March–May 2003, spanning the spring bloom. From observations and numerical simulations the circulation was decomposed into tidal and subtidal components; the former was dominated by the M2 tide, the latter by atmospheric forcing. Near the surface, the subtidal velocity was correlated with wind. At mid depths, the subtidal velocity was organized into vertical bands arising from internal waves excited by wind forcing of the water surface. The tidal flow was more strongly steered by local bathymetry and weaker in peak magnitudes than the subtidal flow, yet it contributed more mechanical energy to the inlet. Tidal eddies reduce exchange of water through the inlet’s entrances. Numerical simulations with the Princeton Ocean Model recreated many observed features, including the three-layer vertical structure of outflow at the surface and bottom and inflow at mid depth, the mid-depth subtidal response to the wind, and characteristics of the tide. While the model produced greater subtidal flow magnitudes at depth and differences in the phase of the M2 tide compared to observations, overall the case study provided support for more comprehensive simulations of Puget Sound in the future.     

基于主成分分析和BP神经网络的五道沟水面蒸发计算研究

针为探讨五道沟水面蒸发量与气象因子间的关系,准确估算该地区水面蒸发量,选取五道沟1991—2019年水面蒸发量以及气象因子实测资料,基于主成分法分析水面蒸发量的影响因素,并结合BP人工神经网络算法建立了水面蒸发计算模型.结果表明:主成分分析提取了三个主成分,第一主成分为地表温度、饱和差、绝对湿度、平均气温以及水汽压力差...     

Flow balances in St. Andrew Bay revealed through hydrodynamic simulations

The circulation patterns in St. Andrew Bay, Florida are revealed through the application of a well-tested, extensively used three-dimensional hydrodynamic model. A high resolution grid resolving both the horizontal and vertical directions is used with a systematically developed set of forcing functions to simulate conditions over a full year. Water levels at the three open boundaries are deduced from a year-long deployment of pressure gauges, and freshwater loadings are based upon drainage basin characteristics and precipitation measurements. Model validation involves comparisons with hydrographic casts taken at twelve stations distributed throughout the bay at monthly intervals. The relative average error between the observed and model-predicted salinity is 15% for the surface of the water column and 4% for the bottom. The annual net flow balance consists of an influx of water at the two Intracoastal Waterway open boundaries, with that water exiting to the Gulf of Mexico. An average of about 100 m³ s⁻¹ enters from East Bay and about 40 m³ s⁻¹ enters through West Bay. On shorter time scales, the flow balance is quite variable both in terms of magnitude and direction. This study also presents methods to overcome the paucity of data that is usually available for the development of such a model. These include techniques to take bottom pressure data sets with short gaps and create reliable sea surface elevation boundary conditions and to take precipitation data and drainage basin characteristics and produce estimates of freshwater inflows.     

Runoff and water budget of the Liudaogou Catchment at the wind–water erosion crisscross region on the Loess Plateau of China

The objectives of this study were to examine the runoff characteristics and to estimate water budget at the wind–water erosion crisscross region on the Loess Plateau of China. A small catchment known as Liudaogou that has representative meteorological and hydrological conditions of the wind–water erosion crisscross region was chosen as the study location. A numerical model for rainfall-runoff was developed and verified; rainfall-runoff calculation for 5 years (2005–2009) was performed. The observed data and numerical result of the surface runoff were used for evaluating runoff characteristics and estimating the annual water budget. Runoff rate was proportional to average intensity of rain. Even though rainfall duration was for few minutes, surface runoff was generated by intensity of more than 2.6 mm × 5 min^?1, when rainfall duration exceeded 10 h; surface runoff was generated by an intensity of 0.6 mm × 5 min^?1, while annual runoff rate was 10–15 %. The unit area of 1 km² was adopted as the index area for estimating annual water budget. Runoff, evapotranspiration, variation of water storage, and habitant water consumption accounted for 20.4, 75.6, 0, and 4 % of the total annual precipitation, respectively. Results of this study provide the basis for further research on hydrology, water resources, and sustainable water development and utilization at the wind–water erosion crisscross region on the northern Loess Plateau where annual water resources are relatively deficient.     

Using remotely sensed data to modify wind forcing in operational storm surge forecasting

Storm surges are abnormal coastal sea level events caused by meteorological conditions such as tropical cyclones. They have the potential to cause widespread loss of life and financial damage and have done so on many occasions in the past. Accurate and timely forecasts are necessary to help mitigate the risks posed by these events. Operational forecasting models use discretisations of the governing equations for fluid flow to model the sea surface, which is then forced by surface stresses derived from a model wind and pressure fields. The wind fields are typically idealised and generated parametrically. In this study, wind field datasets derived from remotely sensed data are used to modify the model parametric wind forcing and investigate potential improvement to operational forecasting. We examine two methods for using analysis wind fields derived from remotely sensed observations of three hurricanes. Our first method simply replaces the parametric wind fields with its corresponding analysis wind field for a period of time. Our second method does this also but takes it further by attempting to use some of the information present in the analysis wind field to estimate future wind fields. We find that our methods do yield some forecast improvement, most notably for our second method where we get improvements of up to 0.29 m on average. Importantly, the spatial structure of the surge is changed in some places such that locations that were previously forecast small surges had their water levels increased. These results were validated by tide gauge data.     

0608号台风“桑美”过境前后对长江口外海域环境的影响

台风是短期改变海洋环境的重要动力因素,会对近岸河口环境产生巨大的影响.基于0608号台风“桑美”过境前后长江口外海域多个站位的调查资料,分析了台风过境前后长江口外海域温度、盐度、浊度等环境要素的变化.结果表明:在台风引起的水体平流输运和垂向混合作用下,长江口外海域的温度、盐度、浊度等结构发生了较大改变.表层水体的温度降低,盐度增大;底层水体温度升高,盐度降低;温跃层厚度增大,层结构强度减弱;东部海域的盐跃层基本消失,温盐跃层对营养盐垂向输运的限制减弱,导致叶绿素a浓度值升高;台湾暖流底层冷水明显向近岸涌升,近岸上升流特征更为显著,同时台湾暖流底层冷水的中心位置向东偏移;水体浊度减小,尤其是近岸水深较浅海域,细颗粒沉积物被强平流作用输运.      

Evolution of salinity field in the upper layers of the east central Arabian Sea and northern Bay of Bengal during summer monsoon experiments

The intra-seasonal variability observed in the salinity field of the upper layers at a few locations in the east central Arabian Sea and the northern Bay of Bengal during the summer monsoon seasons of 1977 and 1979 is documented with the aid of short time series (1–2 weeks) of salinity measurements made from USSR and Indian ships deployed during MONSOON-77 (1977) and MONEX-79 (1979) field experiments. In the Arabian Sea a typical subsurface maxima observed beneath the mixed layer base either disappeared or considerably weakened due to strong vertical mixing caused by the monsoonal forcing. In the northern Bay of Bengal the salinity variability in the top 30 m water column was rapid and appeared to be influenced by large amounts of fresh water from rain and probably from the major adjoining rivers. Some simple diagnostic calculations are presented to assess the relative importance of various processes which control the observed salinity variability.     

Control of estuarine stratification and mixing by wind-induced straining of the estuarine density field

Observations from the York River Estuary, Virginia, demonstrate that the along-channel wind plays a dominant role in governing the estuarine exchange flow and the corresponding increase or decrease in vertical density stratification. Contrary to previous findings that suggest wind stress acts predominantly as a source of energy to mix away estuarine stratification, our results demonstrate that the wind can play a more important role in straining the along-channel estuarine density gradient. Down-estuary winds enhance the tidally-averaged vertical shear, which interacts with the along-channel density gradient to increase vertical stratification. Up-estuary winds tend to reduce, or even reverse the vertical shear, reducing vertical stratification. In two experiments each lasting approximately a month, the estuarine exchange flow was highly correlated with the along-channel component of the wind. The changes in stratification caused by the exchange flow appear to control the amount of vertical mixing as parameterized by the vertical eddy viscosity. The degree of stratification induced by wind straining also appears to play an important role in controlling the effectiveness of wind and tidal mixing.     

Hydrochemistry and origins of mineralized waters in the Puebla aquifer system,Mexico

Significant upward movement of mineralized water takes place in the Puebla aquifer system. Preferential groundwater flow paths related to the geological structure and the lowering of the potentiometric surface are suspected to be the prime factors for this intrusion. A combined approach of geochemical and isotope analyses was used to assess the sources of salinity and processes that are controlling the changes in groundwater chemical composition in the Puebla aquifer. Geochemical and isotope data indicate that the likely source of increased solutes is mineralized water from the dissolution of evaporites of the Cretaceous age at the base of the Upper deep aquifer, which is deeper than the intakes of the shallow wells. Dedolomitization and cation exchange seems also to occur along flow paths where sulphate concentrations tend to increase. The deep regional flow paths controls the chemical stratification of groundwater in response to decreased heads through interconnecting vertical and horizontal pathways, such as in the Fosa Atlixco. The results also suggest that high sulphate concentrations originating in the Lower deep aquifer are currently affecting shallow production wells. It is concluded that hydrodynamic aspects together with hydrogeochemical characteristics need to be taken into account to correctly explain the hydrochemical evolution in the stratified aquifer.