Chemical and isotopic assessment of surface water–shallow groundwater interaction in the arid Grande river basin,North-Central Chile

ABSTRACT

This paper analyses the composition of surface water and shallow groundwater in the Grande River basin, North-Central Chile, using this information to characterize water interactions. Chemical and isotopic data for surface water and groundwater (7 and 6 sampling locations, respectively) were obtained from three sampling campaigns performed in March–April (autumn), August–September (late winter) and December (early summer) 2012. Precipitation samples were also collected. Data was processed using spatial distribution charts, Piper and Stiff diagrams, and multivariate analysis. In general, the results for each method converge on a high degree of connectivity between surface water and shallow groundwater in the study area. Furthermore, approximately a 10% of groundwater contribution to the surface flow discharge was estimated for a particular reach. This multi-method approach was useful for the characterization of surface water–groundwater interactions in the Grande River basin, and may become a suitable and replicable scheme for studies in arid and semi-arid basins facing similar water management challenges.

Editor D. Koutsoyiannis; Associate editor B. Dewals     

Spatial scale effect on the upper soil effective parameters of a distributed hydrological model

Nonlinear dynamics and spatial variability in hydrological systems make the formulation of scaling theories difficult. Therefore, the development of knowledge related to scale effects, scaling techniques, parameterization and linkages of parameters across scales is highly relevant. The main purpose of this work is to analyse the spatial effect of the static storage capacity parameter H_u and the saturated hydraulic conductivity parameter k_s from microscale (sub‐grid level) to mesoscale (grid level) and its implication to the definition of an optimum cell size. These two parameters describe the upper soil water characteristics in the infiltration process conceptualization of the TETIS hydrological model. At microscale, the spatial heterogeneity of H_u and k_s was obtained generating random parameter fields through probability distribution functions and a spatial dependence model with pre‐established correlation lengths. The effective parameters at mesoscale were calculated by solving the inverse problem for each parameter field. Results indicate that the adopted inverse formulation allows transferring the nonlinearity of the system from microscale to the mesoscale via non‐stationary effective parameters. Their values at each cell and time step are in the range of zero to the mean value of the parameter at microscale. The stochastic simulations showed that the variance of the estimated effective parameters decreases when the ratio between mesoscale cell size and correlation length at microscale increases. For a ratio greater than 1, we found cell sizes having the characteristics of a representative elementary area (REA); in such case, the microscale variability pattern did not affect the system response at mesoscale. Copyright © 2011 John Wiley & Sons, Ltd.     

Impact of sensor measurement error on sensor positioning in water quality monitoring networks

This paper studies the impact of sensor measurement error on designing a water quality monitoring network for a river system, and shows that robust sensor locations can be obtained when an optimization algorithm is combined with a statistical process control (SPC) method. Specifically, we develop a possible probabilistic model of sensor measurement error and the measurement error model is embedded into a simulation model of a river system. An optimization algorithm is used to find the optimal sensor locations that minimize the expected time until a spill detection in the presence of a constraint on the probability of detecting a spill. The experimental results show that the optimal sensor locations are highly sensitive to the variability of measurement error and false alarm rates are often unacceptably high. An SPC method is useful in finding thresholds that guarantee a false alarm rate no more than a pre-specified target level, and an optimization algorithm combined with the thresholds finds a robust sensor network.     

Seasonal copula models for the analysis of glacier discharge at King George Island,Antarctica

Modelling glacier discharge is an important issue in hydrology and climate research. Glaciers represent a fundamental water resource when melting of ice and snow contributes to runoff. Glaciers are also studied as natural global warming sensors. GLACKMA association has implemented one of their Pilot Experimental Catchment areas at the King George Island in the Antarctica which records values of the liquid discharge from Collins glacier. In this paper, we propose the use of time-varying copula models for analyzing the relationship between air temperature and glacier discharge, which is clearly non constant and non linear through time. A seasonal copula model is defined where both the marginal and copula parameters vary periodically along time following a seasonal dynamic. Full Bayesian inference is performed such that the marginal and copula parameters are estimated in a one single step, in contrast with the usual two-step approach. Bayesian prediction and model selection is also carried out for the proposed model such that Bayesian credible intervals can be obtained for the conditional glacier discharge given a value of the temperature at any given time point. The proposed methodology is illustrated using the GLACKMA real data where there is, in addition, a hydrological year of missing discharge data which were not possible to measure accurately due to problems in the sounding.     

A stochastic methodology for risk assessment of a large earthquake when a long time has elapsed

We propose a stochastic methodology for risk assessment of a large earthquake when a long time has elapsed from the last large seismic event. We state an approximate probability distribution for the occurrence time of the next large earthquake, by knowing that the last large seismic event occurred a long time ago. We prove that, under reasonable conditions, such a distribution is exponential with a rate depending on the asymptotic slope of the cumulative intensity function corresponding to a nonhomogeneous Poisson process. As it is not possible to obtain an empirical cumulative distribution function of the waiting time for the next large earthquake, an estimator of its cumulative distribution function based on existing data is derived. We conduct a simulation study for detecting scenario in which the proposed methodology would perform well. Finally, a real-world data analysis is carried out to illustrate its potential applications, including a homogeneity test for the times between earthquakes.     

Biocrust landscape-scale spatial distribution is strongly controlled by terrain attributes: Topographic thresholds for colonization in a semiarid badland system

Biological soil crust, or biocrust communities, are the dominating life form in many extreme habitats, such as arid and semiarid badlands, where water scarcity and highly erodible substrates limit vegetation cover. While climate, soil and biotic factors have been described as environmental filters influencing biocrust distribution in such biomes, little is known about the effect of terrain attributes on creating specific microhabitats that promote or restrict biocrust colonization. This study aimed to identify the main terrain attributes controlling biocrust distribution in the driest badland system in Europe, the Tabernas Badlands (SE Spain). To do this, we analysed the influence of different terrain attributes related to landscape stability and microclimate formation on the spatial distribution of lichen and cyanobacteria, using field measurements and topographical information from a LiDAR survey. Our results showed that the spatial distribution of cyanobacteria-dominated biocrusts, which are physiologically and morphologically adapted to extreme drought and high UVA radiation, was mostly associated with areas of high potential incoming solar radiation. The exception was bare south-aspect hillslopes with very high sediment transport potential, where bare physically crusted soils were the dominant ground cover. Lichen-dominated biocrusts, in contrast, colonized near the top of north-aspect hillslopes, characterized by low potential incoming solar radiation and potential evapotranspiration, and their cover decreased downstream, as conditions became good enough for vascular plants. © 2019 John Wiley & Sons, Ltd. © 2019 John Wiley & Sons, Ltd.     

The age and thermal history of Cerro Rico de Potosi, Bolivia

 Cerro Rico de Potosi, Bolivia, is the world’s largest silver deposit and has been mined since the sixteenth century for silver, and for tin and zinc during the twentieth century, together with by-product copper and lead. The deposit consists primarily of veins that cut an altered igneous body that we interpret to be a dacitic volcanic dome and its underlying tuff ring and explosion breccia. The deposit is compositionally and thermally zoned, having a core of cassiterite, wolframite, bismuthinite, and arsenopyrite surrounded by a peripheral, lower-temperature mineral assemblage consisting principally of sphalerite, galena, lead sulfosalt, and silver minerals. The low-temperature assemblage also was superimposed on the high-temperature assemblage in response to cooling of the main hydrothermal system. Both the dacite dome and the ore fluids were derived from a larger magmatic/hydrothermal source at depth. The dome was repeatedly fractured by recurrent movement on the fault system that guided its initial emplacement. The dome was extruded at 13.8±0.2 Ma (2σ), based on U-Th-Pb dating of zircon. Mineralization and alteration occurred within about 0.3 my of dome emplacement, as indicated by a ⁴⁰Ar/³⁹Ar date of 13.76±0.10 Ma (1σ) for sericite from the pervasive quartz-sericite-pyrite alteration associated with the main-stage, high-temperature, mineralization. The last thermal event able to reset zircon fission tracks occurred no later than 12.5±1.1 Ma (1σ), as indicated by fission-tract dating. Minor sericite, and magmatic-steam alunite veins, were episodically formed around 11 Ma and between 8.3 and 5.7 Ma; the younger episodes occurring at the time of extensional fracturing at Cerro Rico and widespread volcanism in the adjacent Los Frailes volcanic field. None of these younger events appear to be significant thermal/mineralizing events; the exceptionally flat thermal release pattern of ³⁹Ar from sericite and the results of the fission-tract dating of zircon show that none of the younger events was hot enough, and lasted long enough, to cause significant loss of Ar or annealing of zircon fission tracks. U-Th-Pb dating of zircon cores indicates a Precambrian progenitor for some zircons, and REE analyses of dated samples of hydrothermally altered dacite show the presence of a prominent positive Eu anomaly, which constrains interpretations of the origin and evolution of the magmatic/hydrothermal system. Received: 14 October 1995/Accepted: 29 January 1996