Groundwater Quality in the Yucatan Peninsula: Insights from Stable Isotope and Metals Analysis

High surface water-groundwater connectivity characterizes watersheds underlain by karsts, increasing contaminant transport risks. However, karsts are highly complex, making research necessary to understand the transport of contaminants from the surface, through the aquifer, to discharge areas. In Yucatan, the lack of waste water treatment raises the risk of groundwater contamination. We monitored stable isotopes (δ¹⁸O-NO₃ and δ¹⁵N-NO₃), cadmium, and lead to document waste water contamination and transport during the rainy and dry seasons, using water samples collected along the Ring of Cenotes during each season. Specific conductance and pH showed no consistent seasonality, with conductance ranging from 0.5 to 55 mS/cm and pH ranging from 6.6 to 8.6 for most samples. Nitrate concentrations in the cenotes averaged 205 ± 260 μM and no seasonal pattern was observed. Cd and Pb concentrations were 0.1 to 37.9 μg/L and 0.2 to 243.2 μg/L, respectively. Nitrate stable isotope values were 2.6 to 27.2‰ for δ¹⁸O and 1.2 to 20.7‰ for δ¹⁵N. The statistical relationship between δ¹⁵N and δ¹⁸O, in dry season samples, indicated that denitrification was occurring. A scale measure for waste water recognition showed: (1) high variability among sites probably related with dry/rainy seasons and/or diverse anthropogenic activities; and (2) specific water quality variables that contribute to contamination at each site during each season. Importantly, our analyses indicate that in the area surrounding the Ring of Cenotes, waste water exhibits spatial and temporal patterns related to complex transport and dilution processes, as is the case in karsts in general.     

The active layer in gravel-bed rivers: An empirical appraisal

The vertical position of the streambed–water boundary fluctuates during the course of sediment transport episodes, due to particle entrainment/deposition and bedform migration, amongst other hydraulic and bedload mechanisms. These vertical oscillations define a topmost stratum of the streambed (i.e. the ‘active layer or active depth’), which usually represents the main source of particles entrained during long and high-magnitude bedload transport episodes. The vertical extent of this layer is hence a capital parameter for the quantification of bedload volumes and a major driver of stream ecology in gravel-bed rivers. However, knowledge on how the active depth scales to flow strength and the nature of the different controls on the relation between the flow strength and the active depth is still scarce. In this paper we present a meta-analysis over active depth data coming from ~130 transport episodes extracted from a series of published field studies. We also incorporate our own field data for the rivers Ebro and Muga (unpublished), both in the Iberian Peninsula. We explore the database searching for the influence of flow strength, grain size, streambed mobility and channel morphology on the vertical extent of the active layer. A multivariate statistical analysis (stepwise multiple regression) confirms that the set of selected variables explains a significant amount of variance in the compiled variables. The analysis shows a positive scaling between active depth and flow strength. We have also identified some links between the active depth and particle travel distances. However, these relations are also largely modulated by other fluvial drivers, such as the grain size of the bed surface and the dominant channel macro-bedforms, with remarkable differences between plane-bed, step-pool and riffle-pool channels. © 2020 John Wiley & Sons, Ltd.     

Multi-Gaussian kriging and simulation in the presence of an uncertain mean value

The multi-Gaussian model is used in geostatistical applications to predict functions of a regionalized variable and to assess uncertainty by determining local (conditional to neighboring data) distributions. The model relies on the assumption that the regionalized variable can be represented by a transform of a Gaussian random field with a known mean value, which is often a strong requirement. This article presents two variations of the model to account for an uncertain mean value. In the first one, the mean of the Gaussian random field is regarded as an unknown non-random parameter. In the second model, the mean of the Gaussian field is regarded as a random variable with a very large prior variance. The properties of the proposed models are compared in the context of non-linear spatial prediction and uncertainty assessment problems. Algorithms for the conditional simulation of Gaussian random fields with an uncertain mean are also examined, and problems associated with the selection of data in a moving neighborhood are discussed.     

The impact of agriculture on solute and suspended sediment load on a Mediterranean watershed after intense rainstorms

The solute and suspended‐sediment load following five rainstorms (2005–2007) with varied intensities were studied at the Vernegà experimental watershed, north‐western Spain. Two land‐use areas are located within this watershed, the upstream one (forest) with 160 ha a 100% forested area, and the downstream one (agricultural) with 97 ha being 9 ha conventional agricultural field and 88 ha forest. This study investigates the capacity of each land‐use to yield water, suspended sediment concentration (SSC) and dissolved solid concentration (DSC). The hypothesis is that DSC and SSC from the agricultural area are greater than DSC and SSC of the forest area. Results showed that the agriculture area produced significantly greater mean DSC than in the forest area, the main contribution was the Ca²⁺ (24·68 ± 46·52 mg l^?1) ion at the agricultural area. A long‐term sediment production rate at the agricultural outlet was calculated (69·1 tonnes per 100 years) based on the total sediment discharge (TSD) and the recurrence interval of the largest event of the five rainstorms (October 2005). Geographic information system (GIS) spatial data layers of the watershed were produced to determine the relation of tracks, landforms, slopes and forest management to SSC yield in the forest outlet (133·89 ± 308·14 mg l^?1) during the five rainstorms. Agriculture practices are the main cause of soil erosion at the study area. Copyright © 2010 John Wiley & Sons, Ltd.     

Comparison of two different solutions to Molodensky’s G1 term

The main purpose of this paper is comparison of two different approaches of solution to the Simple Molodensky’s Problem, the Molodensky’s Approach and the Analytical Continuation Approach, based on numerical computation. Although these approaches have been described theoretically by several authors, e.g. Molodensky et al. (1960), Heiskanen and Moritz (1967), Vaníček (1974), Moritz (1980) and Holota (1991, unpublished results) and theoretical proof of equivalence was given by Heiskanen and Moritz (1967), Moritz (1971), Ecker (1971) and Pellinen (1972, unpublished results), only very few practical experiences about the differences between particular solutions and computational efficiency exist. In this paper we compare the above two mentioned approaches in terms of the G₁-effect on quasigeoid. Both quasigeoid solutions are tested by the independent GPS/levelling approach and are also compared with the previous quasigeoid model of Slovakia where the G₁-term has been approximated using the classical terrain correction. The effect of the G₂-term is also numerically estimated, revealing that it might be significant for the precise quasigeoid determination. Some practical comments and recommendations are given at the end of the paper.     

An improved spectral turning-bands algorithm for simulating stationary vector Gaussian random fields

We propose a spectral turning-bands approach for the simulation of second-order stationary vector Gaussian random fields. The approach improves existing spectral methods through coupling with importance sampling techniques. A notable insight is that one can simulate any vector random field whose direct and cross-covariance functions are continuous and absolutely integrable, provided that one knows the analytical expression of their spectral densities, without the need for these spectral densities to have a bounded support. The simulation algorithm is computationally faster than circulant-embedding techniques, lends itself to parallel computing and has a low memory storage requirement. Numerical examples with varied spatial correlation structures are presented to demonstrate the accuracy and versatility of the proposal.     

Twenty Years Seismic Monitoring of Induced Seismicity in Northern Albania

—The northern part of Albania has been the focus of an intense effort by the Seismological Institute of Albania, for in this area two of the country’s biggest reservoirs are located. Three years before the impounding of the Fierza reservoir, a four-station network was installed around it. The possibility of induced seismicity continued after the impoundment of the Fierza reservoir in 1978 and the Komani reservoir, in 1985. The seismicity of the zone and some aspects of induced seismicity including temporal correlation of seismicity with water level changes, spatial patterns in seismicity, frequency-mag nitude relations, fault plane solutions etc., are studied during this period. The presence in this zone of a very important transverse fault, the Shkoder-Peja fault, makes the study of induced seismicity from Fierza and Komani reservoirs even more significant. The studies have shown that the impounding of the Fierza and Komani reservoirs has modified the natural course of microearthquake energy release, increasing the number of swarms in this area.¶The fluctuation of the water level in these two reservoirs, in due course, is a potentially important factor in the evaluation of seismicity for Northern Albania and especially in the hazard assessment of this region.     

Response of single benthic metrics and multi-metric methods to anthropogenic pressure gradients, in five distinct European coastal and transitional ecosystems

In recent times many benthic indices have been proposed to assess the ecological quality of marine waters worldwide. In this study we compared single metrics and multi-metric methods to assess coastal and transitional benthic status along human pressure gradients in five distinct environments across Europe: Varna bay and lake (Bulgaria), Lesina lagoon (Italy), Mondego estuary (Portugal), Basque coast (Spain) and Oslofjord (Norway). Hence, 13 single metrics (abundance, number of taxa, and several diversity and sensitivity indices) and eight of the most common indices used within the European Water Framework Directive (WFD) for benthic assessment were selected: index of size spectra (ISS), Benthic assessment tool (BAT), Norwegian quality index (NQI), Multivariate AMBI (M-AMBI), Benthic quality index (BQI), (Benthic ecosystem quality index (BEQI), Benthic index based on taxonomic sufficiency (BITS), and infaunal quality index (IQI). Within each system, sampling sites were ordered in an increasing pressure gradient according to a preliminary classification based on professional judgement. The different indices are largely consistent in their response to pressure gradient, except in some particular cases (i.e. BITS, in all cases, or ISS when a low number of individuals is present). Inconsistencies between indicator responses were most pronounced in transitional waters (i.e. IQI, BEQI), highlighting the difficulties of the generic application of indicators to all marine, estuarine and lagoonal environments. However, some of the single (i.e. ecological groups approach, diversity, richness) and multi-metric methods (i.e. BAT, M-AMBI, NQI) were able to detect such gradients both in transitional and coastal environments, being these multi-metric methods more consistent in the detection than single indices. This study highlights the importance of survey design and good reference conditions for some indicators. The agreement observed between different methodologies and their ability to detect quality trends across distinct environments constitutes a promising result for the implementation of the WFD’s monitoring plans. Moreover, these results have management implications, regarding the dangers of misclassification, uncertainty in the assessment, use of conflicting indices, and testing and validation of indices.     

Pathogenic fungi: an unacknowledged risk at coastal resorts? New insights on microbiological sand quality in Portugal

Whilst the potential impact on beach users from microorganisms in water has received considerable attention, there has been relatively little investigation into microbial contaminants in sand. Thirty three beaches across Portugal were analyzed during a five year period (2006-2010) to determine the presence of yeasts, pathogenic fungi, dermatophytes, total coliforms, Escherichia coli and intestinal enterococci in sand.Our results showed that 60.4% of the samples were positive for fungi and that 25.2% were positive for the bacterial parameters. The most frequent fungal species found were Candida sp. and Aspergillus sp., whereas intestinal enterococci were the most frequently isolated bacteria. Positive associations were detected among analyzed parameters and country-regions but none among those parameters and sampling period.Regarding threshold values, we propose 15 cfu/g for yeasts, 17 cfu/g for potential pathogenic fungi, 8 cfu/g for dermatophytes. Eighty four cfu/g for coliforms, 250 cfu/g for E. coli, and 100 cfu/g for intestinal enterococci.     

Total mercury in terrestrial systems (air-soil-plant-water) at the mining region of San Joaquín,Queretaro, Mexico

Deficient management of cinnabar mining left the San Joaquín region with high concentrations of mercury in its soils (2.4 – 4164 mg kg^-1). Numerous cinnabar mines have contributed to the dispersion of mercury into agricultural (0.5 –314 mg kg^-1) and forest (0.2 – 69 mg kg^-1) soils. Sediments are a natural means of transportation for mercury, causing its spreading, especially in areas near mine entrances (0.6 – 687 mg kg^-1). The nearness of maize crops to mines favors mercury accumulation in the different plant structures, such as roots, stems, leaves, and grain (0.04 – 8.2 mg kg^-1); these being related to mercury volatilization and accumulation in soils. Mercury vapor present in the settlements could indicate a constant volatilization from lands and soils (22 – 153 ng m^-3). The mercury levels found in the soils, in maize grain, and in the air resulted greater than the standards reported by the Official Mexican Norm (NOM) and the World Health Organization (WHO). Mercury in rainwater is due mainly to the presence of suspended atmospheric particles, later deposited on the surface (1.5 – 339 μg |^-1). Mercury dissolution was found in the drinking water (10 – 170 ng |^-1), with concentrations below those established by the NOM and the WHO. The contamination existing in the San Joaquín region does not reach the levels of the world’s greatest mercury producers: Almaden (Spain) and Idrija (Slovenia). It is, however, like that found in other important second degree world producers such as Guizhou (China). The population of San Joaquín, as well as its surrounding environment, are constantly exposed to mercury contamination, thus making a long term monitoring necessary to determine its effects, especially to people.