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     

Validation of Landslide Susceptibility Maps; Examples and Applications from a Case Study in Northern Spain

A procedure for validating landslide susceptibility maps wasapplied in a study area in northern Spain and the results obtained compared. Validationwas used to carry out sensitivity analysis for individual variables and combinationsof variables. The validity of different map-making methods was tested, as well as theutility of different types of Favourability Functions. The results obtained show thatvalidation is essential to determine the predictive value of susceptibility maps. Italso helps to better select the most suitable function and significant variables, thus improving the efficiency of the mapping process. Validation based on a temporal strategy makes it possible to derive hazard maps from susceptibility maps.     

The impact of hydrological changes on travertine deposits related to thermal springs in the Pamukkale area (SW Turkey)

Pamukkale thermal waters (35 °C), exhibiting calcium-bicarbonate-sulfate composition and high carbon dioxide concentration, are of a predominantly meteoric origin. The meteoric fluid, circulating through faults and fractures, is heated by magmatic intrusions at great depth, and ascends from deep reservoirs to the surface. Mixing with relatively cold groundwater in the near surface zone promotes different saturation conditions with respect to calcium carbonate that later precipitates at depth and/or the surface. Dissolution-deposition processes of calcium carbonate both at surface and depth environments may help to reconstruct past climate direction in the field. During wet climate conditions a high-rate of calcium carbonate accumulation would be expected to occur at the surface because thermal fluid would be under-saturated with respect to calcium carbonate at depth because of a relatively higher mixing ratio with cold groundwater. During dry climate conditions the thermal fluid would be super-saturated at depth because of the highly acidic environment. Hydrometeorological studies reveal that the annual precipitation at the Pamukkale hydrothermal field tends to decrease with time. This climatic change in the area was also detected from geological records. While humid climate conditions prevailed during the late Quaternary, the area has recently been affected by arid/semi-arid climate conditions, followed by some episodic transitions. This study has shown how the system has possibly reacted to different climate conditions since antiquity.     

Diagenetic paths in the margin of a Triassic Basin: NW zone of the Iberian Chain, Spain

Buntsandstein deposits generated in a slowly subsiding basin on the western margin of the Iberian Chain are represented by a stratigraphic succession of fluvial deposits less than 100 m thick (conglomerates, sandstones, and shales). Diagenetic processes in sandstones can be grouped as eodiagenetic, mesodiagenetic, and telodiagenetic. Eodiagenesis can be associated with Muschelkalk, Keuper, and probably early Jurassic times. Mesodiagenesis is probably related to Jurassic times. Diagenetic chemical reactions suggest a maximum burial less than 1.5 km and low temperatures (<120°C). Patterns of porosity reduction by compaction and cementation suggest four diagenetic stages: (1) Loss of primary porosity by early mechanical compaction; (2) early cementation (K-feldspar and dolomite); (3) dissolution of cements; and (4) framework collapse by re-compaction. These stages are manifested by the presence of two types of sandstone. Type I sandstones present high intergranular volume (mean, 30%). Type II sandstones are characterized by high compactional porosity loss and exhibit low values of intergranular volume (mean, 16.9%). Type II sandstones are associated with the dissolution of cement and later re-compaction of type I sandstones. An intermediate telodiagenetic phase is deduced and related to the sharp unconformity between Lower Cretaceous sediments and the underlying sediments. This suggests that a mechanically unstable framework collapsed during the Cretaceous, generating type II sandstones. The analyzed diagenetic paths have a wide applicability on similar marginal areas of rift basins.     

Hydrogeology and sustainable future groundwater abstraction from the Agua Verde aquifer in the Atacama Desert,northern Chile

The hyper-arid conditions prevailing in Agua Verde aquifer in northern Chile make this system the most important water source for nearby towns and mining industries. Due to the growing demand for water in this region, recharge is investigated along with the impact of intense pumping activity in this aquifer. A conceptual model of the hydrogeological system is developed and implemented into a two-dimensional groundwater-flow numerical model. To assess the impact of climate change and groundwater extraction, several scenarios are simulated considering variations in both aquifer recharge and withdrawals. The estimated average groundwater lateral recharge from Precordillera (pre-mountain range) is about 4,482 m³/day. The scenarios that consider an increase of water withdrawal show a non-sustainable groundwater consumption leading to an over-exploitation of the resource, because the outflows surpasses inflows, causing storage depletion. The greater the depletion, the larger the impact of recharge reduction caused by the considered future climate change. This result indicates that the combined effects of such factors may have a severe impact on groundwater availability as found in other groundwater-dependent regions located in arid environments. Furthermore, the scenarios that consider a reduction of the extraction flow rate show that it may be possible to partially alleviate the damage already caused to the aquifer by the continuous extractions since 1974, and it can partially counteract climate change impacts on future groundwater availability caused by a decrease in precipitation (and so in recharge), if the desalination plant in Taltal increases its capacity.     

Some Insights into Topographic,Elastic and Self-gravitation Interaction in Modelling Ground Deformation and Gravity Changes in Active Volcanic Areas

Surface displacements and gravity changes due to volcanic sources are influenced by medium properties. We investigate topographic, elastic and self-gravitation interaction in order to outline the major factors that are significant in data modelling. While elastic-gravitational models can provide a suitable approximation to problems of volcanic loading in areas where topographic relief is negligible, for prominent volcanoes the rough topography could affect deformation and gravity changes to a greater extent than self-gravitation. This fact requires the selection, depending on local relief, of a suitable model for use in the interpretation of surface precursors of volcanic activity. We use the three-dimensional Indirect Boundary Element Method to examine the effects of topography on deformation and gravity changes in models of magma chamber inflation/deflation. Topography has a significant effect on predicted surface deformation and gravity changes. Both the magnitude and pattern of the geodetic signals are significantly different compared to half-space solutions. Thus, failure to account for topographic effects in areas of prominent relief can bias the estimate of volcanic source parameters, since the magnitude and pattern of deformation and gravity changes depend on such effects.     

Evaluation of methods of spatial interpolation for monthly rainfall data over the state of Rio de Janeiro,Brazil

Five deterministic methods of spatial interpolation of monthly rainfall were compared over the state of Rio de Janeiro, southeast Brazil. The methods were the inverse distance weight (IDW), nearest neighbor (NRN), triangulation with linear interpolation (TLI), natural neighbor (NN), and spline tension (SPT). A set of 110 weather stations was used to test the methods. The selection of stations had two criteria: time series longer than 20 years and period of data from 1960 to 2009. The methods were evaluated using cross-validation, linear regression between values observed and interpolated, root mean square error (RMSE), coefficient of determination (r ²), coefficient of variation (CV, %), and the Willmott index of agreement (d). The results from different methods are influenced by the meteorological systems and their seasonality, as well as by the interaction with the topography. The methods presented higher precision (r ²) and accuracy (d, RMSE) during the summer and transition to autumn, in comparison with the winter or spring months. The SPT had the highest precision and accuracy in relation to other methods, in addition to having a good representation of the spatial patterns expected for rainfall over the complex terrain of the state and its high spatial variability.     

Estimation of Surface-Layer Scaling Parameters in the Unstable Boundary Layer: Implications for Orographic Flow Speed-Up

A method is proposed for estimating the surface-layer depth \((z_s)\) and the friction velocity \((u_*)\) as a function of stability (here quantified by the Obukhov length, L) over the complete range of unstable flow regimes. This method extends that developed previously for stable conditions by Argaín et al. (Boundary-Layer Meteorol 130:15–28, 2009), but uses a qualitatively different approach. The method is specifically used to calculate the fractional speed-up \((\varDelta S)\) in flow over a ridge, although it is suitable for more general boundary-layer applications. The behaviour of \(z_s \left( L\right) \) and \(u_*\left( L\right) \) as a function of L is indirectly assessed via calculation of \(\varDelta S\left( L\right) \) using the linear model of Hunt et al. (Q J R Meteorol Soc 29:16–26, 1988) and its comparison with the field measurements reported in Coppin et al. (Boundary-Layer Meteorol 69:173–199, 1994) and with numerical simulations carried out using a non-linear numerical model, FLEX. The behaviour of \(\varDelta S\) estimated from the linear model is clearly improved when \(u_*\) is calculated using the method proposed here, confirming the importance of accounting for the dependences of \(z_s\left( L \right) \) and \(u_*\left( L \right) \) on L to better represent processes in the unstable boundary layer.     

Mid and late Holocene forest fires and deforestation in the subalpine belt of the Iberian range,northern Spain

The conversion of subalpine forests into grasslands for pastoral use is a well-knownphenomenon, although for most mountain areas the timing of deforestation has not been determined. The presence of charcoal fragments in soil profiles affected by shallow landsliding enabled us to date the occurrence of fires and the periods of conversion ofsubalpine forest into grasslands in the Urbión Mountains, Iberian Range, Spain. We found that the treeline in the highest parts of the northwestern massifs of the Iberian Range(the Urbión, Demanda, Neila, and Cebollera massifs) is currently between 1500 and 1600 m a.s.l., probably because of pastoral use of the subalpine belt, whereas in the past it would have reached almost the highest divides(at approximately 2100–2200 m a.s.l.). The radiocarbon dates obtained indicate that the transformation of the subalpine belt occurred during the Late Neolithic, Chalcolithic, Bronze Age, Iron Age, and Middle Ages. Forest clearing was probably moderate during fires prior to the Middle Ages, as the small size of the sheep herds and the local character of the markets only required small clearings, and therefore more limited fires. Thus, it is likely that the forest recovered burnt areas in a few decades; this suggests the management of the forest and grasslands following a slash-andburn system. During the Middle and Modern Ages deforestation and grassland expansion affected most of the subalpine belt and coincided with the increasing prevalence of transhumance, as occurred in other mountains in the Iberian Peninsula(particularly the Pyrenees). Although the occurrence of shallow landslides following deforestation between the Neolithic and the Roman Period cannot be ruled out, the most extensive shallow landsliding processes would have occurred from the Middle Ages until recent times.     

Exploring seven hundred years of transhumance,climate dynamic,fire and human activity through a historical mountain pass in central Spain

A high-altitude peat sequence from the heart of the Spanish Central System (Gredos range) was analysed through a multi-proxy approach to determine the sensitivity of high-mountain habitats to climate, fire and land use changes during the last seven hundred years, providing valuable insight into our understanding of the vegetation history and environmental changes in a mountain pass close to a traditional route of transhumance. The pollen data indicate that the vegetation was dominated by shrublands and grasslands with scattered pines in high-mountain areas, while in the valleys cereals, chestnut and olive trees were cultivated. Strong declines of high-mountain pines percentages are recorded at 1540, 1675, 1765, 1835 and 1925 cal AD, which may be related to increasing grazing activities and/or the occurrence of anthropogenic fires. The practice of mountain summer farming and transhumance deeply changed and redesigned the landscape of the high altitudes in central Spain (Gredos range) since the Middle Ages, although its dynamics was influenced in some way by climate variability of the past seven centuries.