An iterative technique to delineate protection buffers for wetlands in regions subject to intensive groundwater pumping

Wetlands are vulnerable to groundwater extraction, which has proven detrimental to aquatic ecosystems around the planet. As wetlands rank among the world’s most endangered ecosystems, versatile strategies are required to protect them. This paper provides a modelling-based method to delineate protection buffers in wetlands subject to groundwater extraction. The technique is sufficiently flexible to cater to a wide variety of conditions, and simple enough to underpin management decisions on a daily basis. A numerical model is used to obtain a map of the critical rate of groundwater abstraction, based on the distance between wetlands and suitable discharge thresholds. The outcomes determine the allowed pumping rate at any point under steady and transient-state conditions. A new iteration is developed every time a new pumping allowance is made. This procedure is demonstrated by means of hypothetical scenarios, as well as by a case study application in the Valle del Cauca region, Colombia.     

A Least-squares Approach to Depth Determination from Numerical Horizontal Self-potential Gradients

We have developed a least-squares minimization approach to depth determination of a buried ore deposit from numerical horizontal gradients obtained from self-potential (SP) data using filters of successive window lengths (graticule spacings). The problem of depth determination from SP gradients has been transformed into the problem of finding a solution to a nonlinear equation of the form f(z)=0. Formulas have been derived for vertical and horizontal cylinders and spheres. Procedures are also formulated to estimate the electrical dipole moment and the polarization angle. The method is applied to synthetic data with and without random noise. Finally, the validity of the method is tested on two field examples. In both cases, the depth obtained is found to be in a very good agreement with that obtained from drilling information.     

Non-linear channel–shoal dynamics in long tidal embayments

The dynamics of finite-amplitude bed forms in a tidal channel is studied with the use of an idealized morphodynamic model. The latter is based on depth-averaged equations for the tidal flow over a sandy bottom. The model considers phenomena on spatial scales of the order of the tidal excursion length. Transport of sediment mainly takes place as suspended load. The reference state of this model is characterized by a spatially uniform M₂ tidal current over a fixed horizontal bed. The temporal evolution of deviations from this reference state is governed by amplitude equations: these are a set of non-linear equations that describe the temporal evolution of bed forms. These equations are used to obtain new morphodynamic equilibria which may be either static or time-periodic. Several of these bottom profiles show strong similarity with the tidal bars that are observed in natural estuaries. The dependence of the equilibrium solutions on the value of bottom friction and channel width is investigated systematically. For narrow channels (width small compared to the tidal excursion length) stable static equilibria exist if bottom friction is slightly larger than r_cr. For channel widths more comparable to the tidal excursion length, multiple stable steady states may exist for bottom friction parameter values below r_cr. Regardless of channel width, stable time-periodic equilibria seem to emerge as the bottom friction is increased.Responsible Editor: Jens Kappenberg     

Water quality and presence of pesticides in a tropical coastal wetland in southern Mexico

The water quality of the Pozuelos-Murillo lagoon system in southern Mexico was evaluated during three periods between March and October 2002, with particular emphasis on the detection of organochlorine and organophosphorous pesticide residues in water and sediments. Physicochemical and microbiological parameters of water were also measured and integrated in a water quality index (WQI). Multivariate analysis was used to regionalise the lagoon system as a function of the behaviour of all measured parameters. Solid phase micro-extraction followed by gas chromatography (SPME-GC) was used for pesticide analysis. The concentration of phosphorous was found to be higher than that of nitrogenous compounds. This, besides a deficiency in dissolved oxygen and a high organic matter concentration (as COD), reflects eutrophication processes in some areas of the system. Measured levels of faecal coliforms and oils and greases were above the limits established by Mexican law and comparable to the concentrations reported for other highly polluted systems in Mexico. Residues of DDD (2.0 microg L(-1)) in water and DDE (247 ng g(-1)) and endosulfan I (814 ng g(-1)) sediments were detected by SPME-GC. The spatial distribution of these contaminants implies major potential risks because the most polluted sites were found to be those with the highest fishing activity. Although in general the WQI is on acceptable levels (65-80%), some contamination problems are evident.     

A model for grain-size sorting over tidal sand ridges

A model was developed and analyzed to quantify the effect of graded sediment on the formation of tidal sand ridges. Field data reveal coarse (fine) sediment at the crests (in the troughs), but often phase shifts between the mean grain-size distribution and the bottom topography occur. Following earlier work, this study is based on a linear stability analysis of a basic state with respect to small bottom perturbations. The basic state describes an alongshore tidal current on a coastal shelf. Sediment is transported as bed load and dynamic hiding effects are accounted for. A one-layer model for the bed evolution is used and two grain size classes (fine and coarse sand) are considered. Results indicate an increase in growth and migration rates of tidal sand ridges for a bimodal mixture, whilst the wavelength of the ridges remains unchanged. A symmetrical externally forced tidal current results in a grain-size distribution which is in phase with the ridges. Incorporation of an additional external M₄ tidal constituent or a steady current results in a phase shift between the grain-size distribution and ridge topography. These results show a general agreement with observations. The physical mechanism responsible for the observed grain-size distribution over the ridges is also discussed.Responsible Editor: Jens Kappenberg     

Effective elastic properties of randomly fractured soils: 3D numerical experiments

This paper is concerned with numerical tests of several rock physical relationships. The focus is on effective velocities and scattering attenuation in 3D fractured media. We apply the so‐called rotated staggered finite‐difference grid (RSG) technique for numerical experiments. Using this modified grid, it is possible to simulate the propagation of elastic waves in a 3D medium containing cracks, pores or free surfaces without applying explicit boundary conditions and without averaging the elastic moduli. We simulate the propagation of plane waves through a set of randomly cracked 3D media. In these numerical experiments we vary the number and the distribution of cracks. The synthetic results are compared with several (most popular) theories predicting the effective elastic properties of fractured materials. We find that, for randomly distributed and randomly orientated non‐intersecting thin penny‐shaped dry cracks, the numerical simulations of P‐ and S‐wave velocities are in good agreement with the predictions of the self‐consistent approximation. We observe similar results for fluid‐filled cracks. The standard Gassmann equation cannot be applied to our 3D fractured media, although we have very low porosity in our models. This is explained by the absence of a connected porosity. There is only a slight difference in effective velocities between the cases of intersecting and non‐intersecting cracks. This can be clearly demonstrated up to a crack density that is close to the connectivity percolation threshold. For crack densities beyond this threshold, we observe that the differential effective‐medium (DEM) theory gives the best fit with numerical results for intersecting cracks. Additionally, it is shown that the scattering attenuation coefficient (of the mean field) predicted by the classical Hudson approach is in excellent agreement with our numerical results.     

A seasonal dipolar eddy near Ras Al Hamra (Sea of Oman)

Trajectories and hydrological data from two Argo floats indicate that warm and salty water at 200–300-m depths was ejected from the coast of Oman, near Ras al Hamra, in spring 2008, 2011, and 2012. This warm and salty water, Persian Gulf Water (PGW), once ejected from the coast, recirculated cyclonically in the western Sea of Oman, but also flowed eastward along the Iranian and Pakistani coasts. There, it was expelled seaward by mesoscale eddies as shown by other float data. Seasonal maps of salinity were computed from all available Argo floats; they showed that, in spring, PGW is present in the middle and north of the Sea of Oman, contrary to fall, when the salinity maxima lie southeast of Ras al Hadd. The ejection of PGW from Ras al Hamra is related here to the influence of a mesoscale dipolar eddy which often appears near this cape in spring. The time-averaged and empirical orthogonal functions of altimetric maps over 11 years for this season confirm the frequent presence and the persistence of this feature. From surface currents and hydrology, deep currents were computed via thermal wind balance, and the associated shear and strain fields were obtained. This deformation field is intense near Ras al Hamra, with an offshore direction. This flow structure associated with the mesoscale dipole explains PGW ejection from the coast. This observation suggests that PGW distribution in the Northern Arabian Sea can be strongly influenced by seasonal mesoscale eddies.