Images of a groundwater flow system in Nicaragua

Conceptualization of groundwater flow systems is necessary for water resources planning. Geophysical, hydrochemical and isotopic characterization methods were used to investigate the groundwater flow system of a multi-layer fractured sedimentary aquifer along the coastline in Southwestern Nicaragua. A geologic survey was performed along the 46 km2 catchment. Electrical resistivity tomography (ERT) was applied along a 4.4 km transect parallel to the main river channel to identify fractures and determine aquifer geometry. Additionally, three cross sections in the lower catchment and two in hillslopes of the upper part of the catchment were surveyed using ERT. Stable water isotopes, chloride and silica were analyzed for springs, river, wells and piezometers samples during the dry and wet season of 2012. Indication of moisture recycling was found although the identification of the source areas needs further investigation. The upper-middle catchment area is formed by fractured shale/limestone on top of compact sandstone. The lower catchment area is comprised of an alluvial unit of about 15 m thickness overlaying a fractured shale unit. Two major groundwater flow systems were identified: one deep in the shale unit, recharged in the upper-middle catchment area; and one shallow, flowing in the alluvium unit and recharged locally in the lower catchment area. Recharged precipitation displaces older groundwater along the catchment, in a piston flow mechanism. Geophysical methods in combination with hydrochemical and isotopic tracers provide information over different scales and resolutions, which allow an integrated analysis of groundwater flow systems. This approach provides integrated surface and subsurface information where remoteness, accessibility, and costs prohibit installation of groundwater monitoring networks.

Coastal floods and decadal changes: the climate factor

Observation has widely shown for nearly all last century that the Spanish (Dynamic) Maritime Climate was following around 10 to 11 year cycles in its most significant figure, wind wave, despite it being better to register cycles of 20 to 22 years, in analogical way with the semi-diurnal and diurnal cycles of Cantabrian tides. Those cycles were soon linked to sun activity and, at the end of the century, the latter was related to the Solar System evolution. We know now that waves and storm surges are coupled and that (Dynamic) Maritime Climate forms part of a more complex “Thermal Machine” including Hydrological cycle. The analysis of coastal floods could so facilitate the extension of that experience. According to their immediate cause, simple flood are usually sorted out into flash, pluvial, fluvial, groundwater and coastal types, considering the last as caused by sea waters. But the fact is that most of coastal floods are the result of the concomitance of several former simple types. Actually, the several Southeastern Mediterranean coastal flood events show to be the result of the superposition within the coastal zone of flash, fluvial, pluvial and groundwater flood types under boundary condition imposed by the concomitant storm sea level rise. This work shall be regarded as an attempt to clarify that cyclic experience, through an in-depth review of a past flood events in Valencia (Turia and Júcar basins), as in Murcia (Segura’s) as well.

Digital-simulation and projection of head changes in the Potomac-Raritan-Magothy aquifer system, coastal plain, New Jersey /

"May 1980."

Origins and distribution of saline ground waters in the Floridan aquifer in coastal southwest Florida /

Bibliography: p. 32-34.

The relation between the lowering of the piezometric surface and the rate and duration of discharge of a well using ground water storage /

Mode of access: Internet.

A new analytical solution for water table fluctuations in coastal aquifers with sloping beaches

The Boussinesq equation appears as the zeroth-order term in the shallow water flow expansion of the non-linear equation describing the flow of fluid in an unconfined aquifer. One-dimensional models based on the Boussinesq equation have been used to analyse tide-induced water table fluctuations in coastal aquifers. Previous analytical solutions for a sloping beach are based on the perturbation parameter, epsilon(N) = alphaepsilon cot beta (in which beta is the beach slope, alpha is the amplitude parameter and epsilon is the shallow water parameter) and are limited to tan(-1) (alphaepsilon) much less than beta less than or equal to pi/2. In this paper, a new higher-order solution to the non-linear boundary value problem is derived. The results demonstrate the significant influence of the higher-order components and beach slope on the water table fluctuations. The relative difference between the linear solution and the present solution increases as 6 and a increase, and reaches 7% of the linear solution. (C) 2003 Elsevier Ltd. All rights reserved.

A synthesis of dominant ecological processes in intensive shrimp ponds and adjacent coastal environments in NE Australia

One of the key environmental concerns about shrimp farming is the discharge of waters with high levels of nutrients and suspended solids into adjacent waterways. In this paper we synthesize the results of our multidisciplinary research linking ecological processes in intensive shrimp ponds with their downstream impacts in tidal, mangrove-lined creeks. The incorporation of process measurements and bioindicators, in addition to water quality measurements, improved our understanding of the effect of shrimp farm discharges on the ecological health of the receiving water bodies. Changes in water quality parameters were an oversimplification of the ecological effects of water discharges, and use of key measures including primary production rates, phytoplankton responses to nutrients, community shifts in zooplankton and delta(15)N ratios in marine plants have the potential to provide more integrated and robust measures. Ultimately, reduction in nutrient discharges is most likely to ensure the future sustainability of the industry. (C) 2003 Elsevier Ltd. All rights reserved.

New approximation for free surface flow of groundwater: capillarity correction

An existing capillarity correction for free surface groundwater flow as modelled by the Boussinesq equation is re-investigated. Existing solutions, based on the shallow flow expansion, have considered only the zeroth-order approximation. Here, a second-order capillarity correction to tide-induced watertable fluctuations in a coastal aquifer adjacent to a sloping beach is derived. A new definition of the capillarity correction is proposed for small capillary fringes, and a simplified solution is derived. Comparisons of the two models show that the simplified model can be used in most cases. The significant effects of higher-order capillarity corrections on tidal fluctuations in a sloping beach are also demonstrated. (c) 2004 Elsevier Ltd. All rights reserved.

Steepness expansion for free surface flows in coastal aquifers

Free surface flow of groundwater in aquifers has been studied since the early 1960s. Previous investigations have been based on the Boussinesq equation, derived from the non-linear kinematic boundary condition. In fact, the Boussinesq equation is the zeroth-order equation in the shallow-water expansion. A key assumption in this expansion is that the mean thickness of the aquifer is small compared with a reference length, normally taken to be the linear decay length. In this study, we re-examine the expansion scheme for free surface groundwater flows, and propose a new expansion wherein the shallow-water assumption is replaced by a steepness assumption. A comparison with experimental data shows that the new model provides a better prediction of water table levels than the conventional shallow-water expansion. The applicable ranges of the two expansions are exhibited. (c) 2004 Elsevier B.V. All rights reserved.

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.

The cost of uniqueness in groundwater model calibration

Calibration of a groundwater model requires that hydraulic properties be estimated throughout a model domain. This generally constitutes an underdetermined inverse problem, for which a Solution can only be found when some kind of regularization device is included in the inversion process. Inclusion of regularization in the calibration process can be implicit, for example through the use of zones of constant parameter value, or explicit, for example through solution of a constrained minimization problem in which parameters are made to respect preferred values, or preferred relationships, to the degree necessary for a unique solution to be obtained. The cost of uniqueness is this: no matter which regularization methodology is employed, the inevitable consequence of its use is a loss of detail in the calibrated field. This, ill turn, can lead to erroneous predictions made by a model that is ostensibly well calibrated. Information made available as a by-product of the regularized inversion process allows the reasons for this loss of detail to be better understood. In particular, it is easily demonstrated that the estimated value for an hydraulic property at any point within a model domain is, in fact, a weighted average of the true hydraulic property over a much larger area. This averaging process causes loss of resolution in the estimated field. Where hydraulic conductivity is the hydraulic property being estimated, high averaging weights exist in areas that are strategically disposed with respect to measurement wells, while other areas may contribute very little to the estimated hydraulic conductivity at any point within the model domain, this possibly making the detection of hydraulic conductivity anomalies in these latter areas almost impossible. A study of the post-calibration parameter field covariance matrix allows further insights into the loss of system detail incurred through the calibration process to be gained. A comparison of pre- and post-calibration parameter covariance matrices shows that the latter often possess a much smaller spectral bandwidth than the former. It is also demonstrated that, as all inevitable consequence of the fact that a calibrated model cannot replicate every detail of the true system, model-to-measurement residuals can show a high degree of spatial correlation, a fact which must be taken into account when assessing these residuals either qualitatively, or quantitatively in the exploration of model predictive uncertainty. These principles are demonstrated using a synthetic case in which spatial parameter definition is based oil pilot points, and calibration is Implemented using both zones of piecewise constancy and constrained minimization regularization. (C) 2005 Elsevier Ltd. All rights reserved.