Application of a coupled ground-surface water flow model to simulate periodic groundwater flow influenced by a sloping boundary, capillarity and vertical flows

Theoretical developments as well as field and laboratory data have shown the influence of the capillary fringe on water table fluctuations to increase with the fluctuation frequency. The numerical solution of a full, partially saturated flow equation can be computationally expensive. In this paper, the influence of the capillary fringe on water table fluctuations is simplified through its parameterisation into the storage coefficient of a fully-saturated groundwater flow model using the complex effective porosity concept [Nielsen, P., Perrochet, P., 2000. Water table dynamics under capillary fringes: experiments and modelling. Advances in Water Resources 23 (1), 503-515; Nielsen, P., Perrochet, P., 2000. ERRATA: water table dynamics under capillary fringes: experiments and modelling (Advances in Water Resources 23 (2000) 503-515). Advances in Water Resources 23, 907-908]. The model is applied to sand flume observations of periodic water table fluctuations induced by simple harmonic forcing across a sloping boundary, analogous to many beach groundwater systems. While not providing information on the moisture distribution within the aquifer, this approach can reasonably predict the water table fluctuations in response to periodic forcing across a sloping boundary. Furthermore, he coupled ground-surface water model accurately predicts the extent of the seepage face formed at the sloping boundary. (C) 2005 Elsevier Ltd. All rights reserved.

Field observations of instantaneous water slopes and horizontal pressure gradients in the swash-zone

Field observations of instantaneous water surface slopes in the swash zone are presented. For free-surface flows with a hydrostatic pressure distribution the surface slope is equivalent to the horizontal pressure gradient. Observations were made using a novel technique which in its simplest form consists of a horizontal stringline extending seaward from the beach face. Visual observation, still photography or video photography is then sufficient to determine the surface slope where the free-surface cuts the line or between reference points in the image. The method resolves the mean surface gradient over a cross-shore distance of 5 m or more to within +/- 0.001, or 1/20th -1/100th of typical beach gradients. In addition, at selected points and at any instant in time during the swash cycle, the water surface slope can be determined exactly to be dipping either seaward or landward. Close to the location of bore collapse landward dipping water surface slopes of order 0.05-0.1 occur over a very small region (order 0.5 m) at the blunt or convex leading edge of the swash. In the middle and upper swash the water surface slope at this leading edge is usually very close to horizontal or slightly seaward. Behind the leading edge, the water surface slope was observed to be very close to horizontal or dipping seaward at all times throughout the swash uprush. During the backwash the water surface slope was observed to be always dipping seaward, approaching the beach slope, and remained seaward until a new uprush edge or incident bore passed any particular cross-shore location of interest. The observations strongly Suggest that the swash boundary layer is subject to an adverse pressure gradient during uprush and a favourable pressure gradient during the backwash. Furthermore, assuming Euler's equations are a good approximation in the swash, the observations also show that the total fluid acceleration is negative (offshore) for almost the whole of the uprush and for the entire backwash. The observations are contrary to recent work suggesting significant shoreward directed accelerations and pressure gradients occur in the swash (i.e., delta u/delta t > 0 similar to delta p/delta x < 0), but consistent with analytical and numerical solutions for swash uprush and backwash. The results have important implications for sediment transport modelling in the swash zone.

Characterisation of sea-water intrusion in the Pioneer Valley, Australia using hydrochemistry and three-dimensional numerical modelling

Sea-water intrusion is actively contaminating fresh groundwater reserves in the coastal aquifers of the Pioneer Valley,north-eastern Australia. A three-dimensional sea-water intrusion model has been developed using the MODHMS code to explore regional-scale processes and to aid assessment of management strategies for the system. A sea-water intrusion potential map, produced through analyses of the hydrochemistry, hydrology and hydrogeology, offsets model limitations by providing an alternative appraisal of susceptibility. Sea-water intrusion in the Pioneer Valley is not in equilibrium, and a potential exists for further landward shifts in the extent of saline groundwater. The model required consideration of tidal over-height (the additional hydraulic head at the coast produced by the action of tides), with over-height values in the range 0.5-0.9 m giving improved water-table predictions. The effect of the initial water-table condition dominated the sensitivity of the model to changes in the coastal hydraulic boundary condition. Several salination processes are probably occurring in the Pioneer Valley, rather than just simple landward sea-water advancement from modern sources of marine salts. The method of vertical discretisation (i.e. model-layer subdivision) was shown to introduce some errors in the prediction of watertable behaviour.

Swash-aquifer interaction in the vicinity of the water table exit point on a sandy beach

The coupling of sandy beach aquifers with the swash zone in the vicinity of the water table exit point is investigated through simultaneous measurements of the instantaneous shoreline (swash front) location, pore pressures and the water table exit point. The field observations reveal new insights into swash-aquifer coupling not previously gleaned from measurements of pore pressure only. In particular, for the case where the exit point is seaward of the observation point, the pore pressure response is correlated with the distance between the exit point and the shoreline in that when the distance is large the rate of pressure drop is fast and when the distance is small the rate decreases. The observations expose limitations in a simple model describing exit point dynamics which is based only on the force balance on a particle of water at the sand surface and neglects subsurface pressures. A new modified form of the model is shown to significantly improve the model-data comparison through a parameterization of the effects of capillarity into the aquifer storage coefficient. The model enables sufficiently accurate predictions of the exit point to determine when the swash uprush propagates over a saturated or a partially saturated sand surface, potentially an important factor in the morphological evolution of the beach face. Observations of the shoreward propagation of the swash-induced pore pressure waves ahead of the runup limit shows that the magnitude of the pressure fluctuation decays exponentially and that there is a linear increase in time lags, behavior similar to that of tidally induced water table waves. The location of the exit point and the intermittency of wave runup events is also shown to be significant in terms of the shore-normal energy distribution. Seaward of the mean exit point location, peak energies are small because of the saturated sand surface within the seepage face acting as a "rigid lid'' and limiting pressure fluctuations. Landward of the mean exit point the peak energies grow before decreasing landward of the maximum shoreline position.

The L-31E Surface Water Rediversion Project Final Report: Implementation, Results, and Recommendations

Throughout the Biscayne Bay watershed, existing coastal wetland communities have been cut off from sheet flow for decades. With the expectation that reconnection of these wetlands to upstream water sources would alter existing hydrologic conditions and recreate a more natural sheet flow to Biscayne National Park, a demonstration project on freshwater rediversion was undertaken. The objectives of the project were to document the effects of freshwater diversion on: (a) swamp and nearshore water chemistry and hydrology; (b) soil development processes; (c) macrophyte and benthic algal community composition, structure and production; (d) abundance of epiphytic and epibenthic invertebrates; (e) zonation, production, and phenology of primary producers in the nearshore environment, and (f) exchanges of nutrients and particulates between nearshore and mangrove ecosystems.

Organic matter dynamics in the Florida coastal Everglades: A molecular marker and isotopic approach

Everglades National Park (ENP) is about to undergo the world's largest wetland restoration with the aim of improving the quality, timing and distribution of water flow. The changes in water flow are hypothesized to alter the nutrient fluxes and organic matter (OM) dynamics within ENP, especially in the estuarine areas. This study used a multi-proxy approach of molecular markers and stable δ 13C isotope measurements, to determine the present day OM dynamics in ENP. ^ OM dynamics in wetland soils/sediments have proved to be difficult to understand using traditional geochemical approaches. These are often inadequate to describe the multitude of OM sources (e.g. higher land plant, emergent vegetation, submerged vegetation) to the soils/sediments and the complex diagenetic processes that can alter the OM characteristics. A multi-proxy approach, however, that incorporates both molecular level and bulk parameter information is ideal to comprehend complex OM dynamics in aquatic environments. Therefore, biomass-specific molecular markers or proxies can be useful in tracing the sources and processing of OM. This approach was used to examine the OM dynamics in the two major drainage basins, Shark River Slough and Taylor River Slough, of ENP. Freshwater to marine transects were sampled in both systems for soils/sediments and suspended particulate organic matter (SPOM) to be characterized through bulk OM analyses, lipid biomarker determinations (e.g. sterols, fatty acids, hydrocarbons and triterpenoids) and compound-specific stable carbon isotope (δ 13C) determinations. ^ One key accomplishment of the research was the assessment of a molecular marker proxy (Paq) to distinguish between emergent/higher plant vegetation from submerged vegetation within ENP. This proxy proved to be quite useful at tracing OM inputs to the soils/sediments of ENP. A second key accomplishment was the development of a 3-way model using vegetation specific molecular markers. This novel, descriptive model was successfully applied to the estuarine areas of Taylor and Shark River sloughs, providing clear evidence of mixing of freshwater, estuarine and marine derived OM in these areas. In addition, diagenetic transformations of OM in these estuaries were found to be quite different between Taylor and Shark Rivers, and are likely a result of OM quality and hydrological differences. ^

Hydrological dynamics between a coastal aquifer and the adjacent estuarine system, Biscayne Bay, South Florida

Geochemical and geophysical approaches have been used to investigate the freshwater and saltwater dynamics in the coastal Biscayne Aquifer and Biscayne Bay. Stable isotopes of oxygen and hydrogen, and concentrations of Sr2+ and Ca2+ were combined in two geochemical mixing models to provide estimates of the various freshwater inputs (precipitation, canal water, and groundwater) to Biscayne Bay and the coastal canal system in South Florida. Shallow geophysical electromagnetic and direct current resistivity surveys were used to image the geometry and stratification of the saltwater mixing zone in the near coastal (less than 1km inland) Biscayne Aquifer. The combined stable isotope and trace metal models suggest a ratio of canal input-precipitation-groundwater of 38%–52%–10% in the wet season and 37%–58%–5% in the dry season with an error of 25%, where most (20%) of the error was attributed to the isotope regression model, while the remaining 5% error was attributed to the Sr2+/Ca2+ mixing model. These models suggest rainfall is the dominate source of freshwater to Biscayne Bay. For a bay-wide water budget that includes saltwater and freshwater mixing, fresh groundwater accounts for less than 2% of the total input. A similar Sr 2+/Ca2+ tracer model indicates precipitation is the dominate source in 9 out of 10 canals that discharge into Biscayne Bay. The two-component mixing model converged for 100% of the freshwater canal samples in this study with 63% of the water contributed to the canals coming from precipitation and 37% from groundwater inputs ±4%. There was a seasonal shift from 63% precipitation input in the dry season to 55% precipitation input in the wet season. The three end-member mixing model converged for only 60% of the saline canal samples possibly due to non-conservative behavior of Sr2+ and Ca2+ in saline groundwater discharging into the canal system. Electromagnetic and Direct Current resistivity surveys were successful at locating and estimating the geometry and depth of the freshwater/saltwater interface in the Biscayne Aquifer at two near coastal sites. A saltwater interface that deepened as the survey moved inland was detected with a maximum interpreted depth to the interface of 15 meters, approximately 0.33 km inland from the shoreline. ^