Groundwater-flow modeling in the Yucatan karstic aquifer, Mexico

The current conceptual model of the unconfined karstic aquifer in the Yucatan Peninsula, Mexico, is that a fresh-water lens floats above denser saline water that penetrates more than 40 km inland. The transmissivity of the aquifer is very high so the hydraulic gradient is very low, ranging from 7–10 mm/km through most of the northern part of the peninsula. The computer modeling program AQUIFER was used to investigate the regional groundwater flow in the aquifer. The karstified zone was modeled using the assumption that it acts hydraulically similar to a granular, porous medium. As part of the calibration, the following hypotheses were tested: (1) karstic features play an important role in the groundwater-flow system; (2) a ring or belt of sinkholes in the area is a manifestation of a zone of high transmissivity that facilitates the channeling of groundwater toward the Gulf of Mexico; and (3) the geologic features in the southern part of Yucatan influence the groundwater-flow system. The model shows that the Sierrita de Ticul fault, in the southwestern part of the study area, acts as a flow barrier and head values decline toward the northeast. The modeling also shows that the regional flow-system dynamics have not been altered despite the large number of pumping wells because the volume of water pumped is small compared with the volume of recharge, and the well-developed karst system of the region has a very high hydraulic conductivity. Electronic Publication     

Simultaneous use of geochemical and geophysical methods to characterise abandoned landfills

. Groundwater sampling and geophysical methods determined a serious contamination problem associated with refilled exploitation sites at the Cal Dimoni area, Llobregat delta, Barcelona, Spain. To characterise this process, hydrogeochemical analyses were performed and showed the following modifications to groundwater chemical composition: increasing pH values, changing redox conditions, significant increases in total organic carbon (TOC) and certain trace elements, and high groundwater conductivity values. Major ion content accumulations were found under the refilled area. In contrast, elements involved in the oxidation–reduction processes, such as iron, manganese and nitrates, clearly diminished. Electromagnetic prospecting methods were also performed and delineated the contamination plume extent. These methods also showed separate sources of contamination, one clearly related to the groundwater–refilled zone leachate interaction, another as a consequence of the manure–accumulation surface site. Geochemical and geophysical methods have shown similar results for locating groundwater contamination sources, and for determining leachate generation mechanisms and flow paths.     

U–Pb evidence for a polyorogenic evolution of the HP–HT units of the NW Iberian Massif

. A isotope dilution thermal ionisation mass spectrometry U-Pb geochronological study was carried out on the high-pressure and high-temperature units (HP-HT units) overlying the oceanic suture in the Allochthonous Complexes of the NW Iberian Variscan Belt. The rocks investigated are seven granulite- to eclogite-facies paragneisses and one leucosome within mafic high-pressure granulites in the Ordenes and Cabo Ortegal Complexes of NW Spain. U-Pb dating of zircon, monazite, titanite and rutile reveal the presence of a pervasive Early Ordovician metamorphic event at ca. 500-480 Ma and a later Early Devonian event at ca. 400-380 Ma. The U-Pb ages, in conjunction with petrological and structural data, indicate that the high-pressure event recorded by these rocks is Early Ordovician in age. Monazite ages in the paragneisses suggest that peak metamorphic conditions were reached at ca. 500-485 Ma. Subsequently, the rock ensemble underwent exhumation accompanied by partial melting and zircon growth at ca. 485-470 Ma. Melting of mafic granulites was coeval with this latter episode as indicated by zircon crystallisation age in the leucosomes dated at ca. 486 Ma. Based on these data and on the general features of magmatism and metamorphic evolution, it is proposed that this process took place at a convergent plate boundary within a peri-Gondwanan oceanic domain. Monazite, titanite and rutile data in some of the samples studied show evidence of a second metamorphic episode that took place between ca. 400 and 380 Ma (with a peak at ca. 390-385 Ma). This Early Devonian event, at variance with the previous one, was not pervasive, but, rather, was localised in areas of intense Variscan tectonothermal reworking. It is claimed that this later metamorphic event was recorded by the U-Pb system in areas where monazite and titanite growth was enhanced by fluid circulation in highly strained rocks (Variscan shear zones). According to previous structural studies and Ar-Ar dating of fabrics, this Early Devonian episode took place as the HP-HT units were deformed and thrusted upon the ophiolitic units in the early stages of the Variscan collision.     

Spatial and temporal variation of the nekton and hyperbenthos from a temperate European estuary with regulated freshwater inflow

The aquatic macrofauna of the Guadalquivir estuary were sampled (1 mm mesh persiana net) at 5 sampling sites located along the entire (except the tidal freshwater region) estuarine gradient of salinity (outer 50 km). A total of 134 fish and macroinvertebrate species was collected but only 62 were considered common or regularly present in the estuary. Univariate measures of the community structure showed statistically significant differences among sampling sites: species richness, abundance, and biomass decreased in the upstream direction, being positively correlated with the salinity. Temporal differences of these three variables were also statistically significant. While a clear seasonal pattern (minimum densities in winter and maximum in spring-summer) was observed for abundance and biomass, no such pattern existed for the number of species. Mysids was the most dominant group throughout the estuary (96% to 99% of abundance; 49% to 85% of biomass), although fish biomass was also important at the outer estuary (36% to 38%). Multivariate analyses indicated highly significant spatial variation in the macrofaunal communities observed along the salinity gradient. These analyses suggest that the underlying structure was a continuum with more or less overlapping distributions of the species dependent on their ability to tolerate different physicochemical conditions. There were also significant temporal (intermonthly + interannual) variation of the estuarine community; the relative multivariate dispersion indicated that monthly variation was more considerable (relative multivariate dispersion >1) at the outer part of the estuary during the wet year (last 20 km) and was higher in the inner stations during the dry year (32 to 50 km from the river mouth). Since a clear negative exponential relationship was observed between the freshwater input (from a dam located 110 km upstream) and water salinity at all sampling stations, it is concluded that the human freshwater management is probably affecting the studied estuarine communities. While the higher seasonal (long-term) stability of the salinity gradient, due to the human control of the freshwater input, may facilitate the recruitment of marine species juveniles during the meteorologically unstable early-spring, the additional (short-term) salinity fluctuations during the warm period may negatively affect species that complete their lifecycle within the estuary.     

The El Salvador earthquakes of January and February 2001: context, characteristics and implications for seismic risk

The small Central American republic of El Salvador has experienced, on average, one destructive earthquake per decade during the last hundred years. The latest events occurred on 13 January and 13 February 2001, with magnitudes M_w 7.7 and 6.6, respectively. The two events, which were of different tectonic origin, follow the patterns of the seismicity of the region although neither event has a known precedent in the earthquake catalogue in terms of size and location. The earthquakes caused damage to thousands of traditionally built houses and triggered hundreds of landslides, which were the main causes of fatalities. The earthquakes have clearly demonstrated trends of increasing seismic risk in El Salvador due to rapid population expansion in areas of high shaking and landslide hazard, exacerbated by deforestation and uncontrolled urbanisation. The institutional mechanisms required for the control of land use and building practice are very weak and present a major obstacle to risk mitigation.     

On the Acoustic-Electromagnetic Analogy for the Reflection-Refraction Problem

The same mathematical theory can be used to describe physical phenomena of different nature. For instance, the wave equation and the related mathematical developments can be used to describe elastic and electromagnetic wave propagation, and it is also extensively used in quantum mechanics. Fresnel's equations are a classical example of the analogy between shear waves and light waves. George Green in the nineteenth century, used analogies to obtain the reflection coefficients for sound waves and light waves, before the advent of the electromagnetic theory of light.In this work, we investigate the mathematical analogy between elastic waves and electromagentic waves. We obtain a complete parallelism for the reflection and refraction problem, considering the most general situation, that is, the presence of anisotropy and attenuation—viscosity in the elastic case and conductivity in the electromagnetic case. The analogy is illustrated with Fresnel's equations, the Brewster and critical angles, the concept of reflectivity and transmissivity, and the corresponding duals fields. The analysis of the elastic-solid theory of reflection applied by Green to light waves, and a brief historical review of wave propagation through the ether, further illustrate the analogy.     

Establishment of a Non-Permanent GPS Network to Monitor the Recent NE-SW Deformation in the Granada Basin (Betic Cordillera,Southern Spain)

The Granada Basin (Central Betic Cordillera), one of the most seismically active areas of the Iberian Peninsula, is currently subjected to NW-SE compression and NE-SW extension. The present day extension is accommodated by normal faults with various orientations but particularly with a NW-SE strike. At the surface, these active NW-SE normal faults are mainly concentrated on the NE part of the Basin. In this part we have selected a 15-km long segment where several active normal faults crop out. Using the marine Tortonian rocks as a reference, we have calculated a minimum extensional rate of 0.15-0.30 mm/year. The observed block rotation, the listric geometry of faults at depth and the distribution of seismicity over the whole Basin, indicate that this rate is a minimum value. In the framework of an interdisciplinary research project a non-permanent GPS-network has been established in the central sector of Betic Cordillera to monitor the crustal deformations. The first two observation campaigns were done in 1999 and 2000.     

Coupled inverse modelling of groundwater flow and mass transport and the worth of concentration data

This paper presents the extension of the self-calibrating method to the coupled inverse modelling of groundwater flow and mass transport. The method generates equally likely solutions to the inverse problem that display the variability as observed in the field and are not affected by a linearisation of the state equations. Conditioning to the state variables is measured by an objective function including, among others, the mismatch between the simulated and measured concentrations. Conditioning is achieved by minimising the objective function by gradient-based methods. The gradient contains the partial derivatives of the objective function with respect to: log conductivities, log storativities, prescribed heads at boundaries, retardation coefficients and mass sources. The derivatives of the objective function with respect to log conductivity are the most cumbersome and need the most CPU-time to be evaluated. For this reason, to compute this derivative only advective transport is considered. The gradient is calculated by the adjoint-state method. The method is demonstrated in a controlled, synthetic study, in which the worth of concentration data is analysed. It is shown that concentration data are essential to improve transport predictions and also help to improve aquifer characterisation and flow predictions, especially in the upstream part of the aquifer, even in the case that a considerable amount of other experimental data like conductivities and heads are available. Besides, conditioning to concentration data reduces the ensemble variances of estimated transmissivity, hydraulic head and concentration.     

Digital model of corrected accumulated flow for peak discharge data acquisition and drainage system design

The design of a drainage system for a roofing slate quarry was implemented by the enhancement of discharge peak estimation, and the uncertainty inevitably associated with the engineering model was reduced.

The development of a topographical, geological, and vegetation cover database developed from a Geographical Information System (GIS) allowed for the definition of the drainage network for a hydraulic system, along with the calculation of the runoff coefficient. This is applied to the digital model of accumulated flow (DMF) as a weight correction coefficient, using a matrix-based model at 5×5 m resolution. The new digital model of corrected accumulated flow (DMCF) is the result of combining the thematic maps with the map of slope <3%, which was previously created from the slope model. It is demonstrated that this new model allows to apply the “Rational Method” on cartographic units defined by the GIS.

The DMCF is compared with other traditional applications of the Rational Method based on the calculation of the discharge peak considering: (1) the drainage basin as a single watershed or (2) defining an average runoff coefficient in each sub-watershed. Both approaches have bigger discharge peaks than those obtained by the DMCF since the slope, lithology, and vegetation cover have average values, and the runoff coefficient is poorly defined, increasing the uncertainty in the discharge peak.