Study of geothermal water intrusion due to groundwater exploitation in the Puebla Valley aquifer system, Mexico

Significant intrusion of geothermal water into fresh groundwater takes place in the Puebla Valley aquifer system, Mexico. The decline in the potentiometric surface due to the overexploitation of the groundwater induces this intrusion. This hydrological system comprises three aquifers located in Plio-Quaternary volcanic sediments and Mesozoic calcareous rocks. The hydraulic balance of the aquifer shows that the annual output exceeds the natural inputs by 12 million m³. Between 1973 and 2002, a drop in the potentiometric surface, with an 80 m cone of depression, was identified in a 5-km-wide area located southwest of the city of Puebla. Chemical analyses performed on water samples since 1990 have shown an increase in total dissolved solids (TDS) of more than 500 mg/L, coinciding with the region showing a cone of depression in the potentiometric surface. A three-dimensional flow and transport model, based on the hydrogeological and geophysical studies, was computed by using the MODFLOW and MT3D software. This model reproduces the evolution of the aquifer system during the last 30 years and predicts for 2010 an additional drawdown in the potentiometric surface of 15 m, and an increase in the geothermal water intrusion.     

T -distributed Random Fields: A Parametric Model for Heavy-tailedWell-log Data1

Histograms of observations from spatial phenomena are often found to be more heavy-tailed than Gaussian distributions, which makes the Gaussian random field model unsuited. A T-distributed random field model with heavy-tailed marginal probability density functions is defined. The model is a generalization of the familiar Student-T distribution, and it may be given a Bayesian interpretation. The increased variability appears cross-realizations, contrary to in-realizations, since all realizations are Gaussian-like with varying variance between realizations. The T-distributed random field model is analytically tractable and the conditional model is developed, which provides algorithms for conditional simulation and prediction, so-called T-kriging. The model compares favourably with most previously defined random field models. The Gaussian random field model appears as a special, limiting case of the T-distributed random field model. The model is particularly useful whenever multiple, sparsely sampled realizations of the random field are available, and is clearly favourable to the Gaussian model in this case. The properties of the T-distributed random field model is demonstrated on well log observations from the Gullfaks field in the North Sea. The predictions correspond to traditional kriging predictions, while the associated prediction variances are more representative, as they are layer specific and include uncertainty caused by using variance estimates.     

Ultra-rapid landscape response and sediment yield following glacier retreat, Icy Bay, southern Alaska

The páramo is a neotropical alpine ecosystem that covers more than 75,000 km² of the northern Andes of Colombia, Ecuador, Venezuela, and Peru. It provides important environmental services: more than 10 million people in the Andean highlands benefit from the water supply and regulation function, which is attributed to the volcanic soils that underlie the ecosystem. The soils are also major carbon sinks of global significance. Severe land use changes and soil degradation threaten both the hydrology and carbon sink function. Nevertheless, soil genesis and properties in the páramo is rather poorly understood, nor are their ecological functions well documented. The impact of the geomorphology of the páramo on soil genesis was studied in the rio Paute basin, south Ecuador. Two toposequences were described and analysed. In each toposequence, four pedons were selected representing summit, backslope, undrained plain situation, and valley bottom positions in the landscape. The soils are classified as Hydric Andosols in the World Reference Base for Soil Resources and Epiaquands or Hydrudands in Soil Taxonomy. They are very acidic and have a high organic matter content, high P deficiency, and Al toxicity. Their water content ranges from 2.64 g g^− 1 at saturation, down to 1.24 g g^− 1 at wilting point, resulting in a large water storage capacity. Two major soil forming processes are identified: (1) volcanic ash deposition and (2) accumulation of organic carbon. Volcanic ash deposits may vary in depth as a result of regional geomorphological factors such as parent material, orientation, slope, and altitude. Organic carbon accumulation is an interaction of both waterlogging, which depends on the position in the landscape, and the formation of organometallic complexes with Al and Fe released during volcanic ash breakdown. Despite the high variability in parent material and topography, the soil is characterised by a notable homogeneity in physico-chemical properties. Statistical analysis reveals that only topographic location has a slight but significant influence on soil pH as well as the organic matter content, saturated conductivity and water retention at high pressure. Finally, the exceptional properties of these soils provide useful insights to improve classification of the Andosols reference group of the FAO World reference Base for Soil Resources.     

Mixing and differentiation in the Oruanui rhyolitic magma, Taupo, New Zealand: evidence from volatiles and trace elements in melt inclusions

Large pyroclastic rhyolites are snapshots of evolving magma bodies, and preserved in their eruptive pyroclasts is a record of evolution up to the time of eruption. Here we focus on the conditions and processes in the Oruanui magma that erupted at 26.5 ka from Taupo Volcano, New Zealand. The 530 km³ (void-free) of material erupted in the Oruanui event is comparable in size to the Bishop Tuff in California, but differs in that rhyolitic pumice and glass compositions, although variable, did not change systematically with eruption order. We measured the concentrations of H₂O, CO₂ and major and trace elements in zoned phenocrysts and melt inclusions from individual pumice clasts covering the range from early to late erupted units. We also used cathodoluminescence imaging to infer growth histories of quartz phenocrysts. For quartz-hosted inclusions, we studied both fully enclosed melt inclusions and reentrants (connecting to host melt through a small opening). The textures and compositions of inclusions and phenocrysts reflect complex pre-eruptive processes of incomplete assimilation/partial melting, crystallization differentiation, magma mixing and gas saturation. ‘Restitic’ quartz occurs in seven of eight pumice clasts studied. Variations in dissolved H₂O and CO₂ in quartz-hosted melt inclusions reflect gas saturation in the Oruanui magma and crystallization depths of ∼3.5–7 km. Based on variations of dissolved H₂O and CO₂ in reentrants, the amount of exsolved gas at the beginning of eruption increased with depth, corresponding to decreasing density with depth. Pre-eruptive mixing of magma with varying gas content implies variations in magma bulk density that would have driven convective mixing. Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users.     

Land-Parcel Land-Use History as a Key to Site Selection for Documenting Soil Contamination Risk： a Case Study from Australian Suburbia

In that orcharding in early-to-mid twentieth century southeastern Australia involved use of certain heavy metal and As compounds in regular pest-control spray procedures, some interest attaches to the possibility that these landparcels are underlain by soils with above-background Cu, Pb and As levels. Interpreta- tion of Land-cover changes allowed land parcels previously occupied by orchards to be identified in the 1950s through time-series air-photos. A comparison of soil analysis results referring to soil samples from control sites, and from land parcels formerly occupied by orchardists, shows that contamination （above-background） levels of cations in the pesticides can be found in the top 6 cm of former orchard soils. It is clear that digital spatial data handling and culturally-informed air photo interpretation has a place in soil contamination studies, land-use planning （with particular reference to re-development） and in administration of public health.     

韧性剪切变形对岩石地球化学行为的制约——以北阿尔金巴什考供韧性剪切带为例

韧性剪切变形对岩石地球化学行为的制约一直是地质学家们探讨的课题。本文以构成北阿尔金红柳沟——拉配泉俯冲碰撞杂岩带与北阿尔金地块边界的巴什考供斜向逆冲型韧性剪切带为例，通过对韧性剪切带内花岗岩变形前后不同变形强度构造岩的地球化学组成进行对比，确定等比线斜率，探讨韧性变形对岩石体积和成分变异的影响。计算结果表明，在糜棱岩化过程中，糜棱岩化花岗岩体积亏损21％，花岗质糜棱岩体积亏损31％。质量平衡计算结果和等比线图表明，韧；陛剪切作用导致SiO2，流失量最大，A12O3、K2O及Ba、Rb、Sr等都有不同程度的丢失，显示出较强的活动性，MnO、P2O5、Sc位于等比线上或附近，表现出相对的稳定性。岩石中活动组分的变异是流体渗滤作用引起的，不活动组分的变异是体．积亏损造成的。     

Late Quaternary Spring-Fed Deposits of the Grand Canyon and Their Implication for Deep Lava-Dammed Lakes

One of the most intriguing episodes in the Quaternary evolution of the Grand Canyon of the Colorado River, Arizona, was the development of vast lakes that are thought to have backed up behind lava erupted into the gorge. Stratigraphic evidence for these deep lava-dammed lakes is expectedly sparse. Possible lacustrine deposits at six areas in the eastern canyon yielded no compelling evidence for sediment deposited in a deep lake. At two of the sites the sediment was associated with late Quaternary spring-fed pools and marshes. Water-lain silt and sand at lower Havasu Creek was deposited 3000 cal yr ago. The deposit contains an ostracode assemblage similar to that living in the modern travertine-dammed pools adjacent to the outcrop. The second deposit, at Lees Ferry, formed in a spring-fed marsh 43,000 cal yr ago, as determined by ¹⁴C and amino acid geochronology. It contains abundant ostracode and mollusk fossils, the richest assemblages reported from the Grand Canyon to date. Our interpretation of these sediments as spring-fed deposits, and their relative youth, provides an alternative to the conventional view that deposits like these were formed in deep lava-dammed lakes that filled the Grand Canyon.     

KARSTIC: a sensitivity method for carbonate aquifers in karst terrain

. Groundwater in karstic aquifers can be dangerously sensitive to contamination. Many cities in the western USA rely on karstic carbonate aquifers for municipal water supplies. For example, Rapid City, South Dakota, pumps more than half of its drinking water from wells in the Madison Limestone. This work examined the sensitivity of karstic aquifers to surface contamination in mountainous terrain. Where karstic carbonate aquifers are exposed at their outcrop areas, they are particularly susceptible to the introduction of contamination through diffuse recharge or through point recharge at swallow holes along streams. Residential developments in mountainous regions of the western USA are encroaching on the recharge areas of karstic aquifers. Many of these residential developments are served by onsite wastewater disposal systems such as septic tanks and drain fields, with the attendant danger of introduction of pathogens from malfunctioning treatment systems above fractured limestone which offers little filtering. Where streams disappear into karstic aquifers at swallow holes, microbial contaminants such as Giardia or Cryptosporidium are a concern, as well as potential spills, leaks, or accidents along roads near these streams. The KARSTIC method developed and modified in this work puts greater emphasis on karst features than previous sensitivity procedures such as the US Environmental Protection Agency's DRASTIC method. The modified method gives increased attention to highly sensitive areas of karstic carbonate aquifers by weighting the synergistic effects of fracturing, karst development, and swallow holes of recharging streams. In a field application, hydrogeologic maps of a watershed in the Black Hills, USA, were digitized into a geographic information system. The resulting sensitivity map and report can be used by planners, managers, and the public as a screening tool for assessing groundwater sensitivity in regions which include karstic aquifers.     

Elemental and isotopic fractionation of Type B calcium-, aluminum-rich inclusions: experiments, theoretical considerations, and constraints on their thermal evolution

Experiments exposing Type B calcium-, aluminum-rich inclusion (CAI)-like melts at high temperatures to high vacuum or reducing hydrogen-rich gas mixtures were used to determine the rates and consequences of elemental and isotopic fractionation by evaporation. Silicon and magnesium were found to evaporate much faster than calcium and aluminum, and the resulting residual liquid trajectories in composition space are reproduced via a thermodynamic model for the saturation vapor pressure of the evaporating species. Isotopic fractionations associated with evaporation were measured for magnesium. The resulting relationship between fraction of magnesium lost and enrichment of the residue in the heavy isotopes of magnesium follows a Rayleigh fractionation curve with a fractionation factor that is close to, but not exactly, the theoretically expected value. The rate of evaporation is found to be a strong function of temperature, oxygen fugacity, and melt composition, which can be understood and modeled in terms of the dependence of the saturation vapor pressures on these variables. The relationship between evaporation rate, which we measure, and calculated saturation vapor involves empirical evaporation coefficients that we find to be significantly less than one (∼0.1). Analytical and numerical models are used to characterize how diffusion in both the melt and in the surrounding gas affects evaporation rates and the degree of chemical and isotopic fractionation. The experimental data and theoretical considerations are combined to give a parameterization of the rates and consequences of evaporation of Type B CAI-like liquids, which is then used to translate the measured isotopic fractionation of Type B CAIs into constraints on their thermal history. Cooling rates of the order of 10°C per hour are indicated.