Shallow groundwater temperature response to climate change and urbanization

Groundwater temperatures, especially in shallow (quaternary) aquifers respond to ground surface temperatures which in turn depend on climate and land use. Groundwater temperatures, therefore, are modified by climate change and urban development. In northern temperate climate regions seasonal temperature cycles penetrate the ground to depths on the order of 10–15 m. In this paper, we develop and apply analytic heat transfer relationships for 1-D unsteady effective diffusion of heat through an unsaturated zone into a flowing aquifer a short distance below the ground surface. We estimate how changes in land use (urban development) and climate change may affect shallow groundwater temperatures. We consider both long-term trends and seasonal cycles in surface temperature changes. Our analysis indicates that a fully urbanized downtown area at the latitude of Minneapolis/St. Paul is likely to have a groundwater temperature that is nearly 3 °C warmer than an undeveloped agricultural area at the same geographic location. Pavements are the main cause of this change. Data collected by the Minnesota Pollution Control Agency (MPCA) in the St. Cloud, MN area confirm that land use influences groundwater temperatures. Ground surface temperatures are also projected to rise in response to global warming. In the extreme case of a doubling of atmospheric carbon dioxide (2 × CO₂ climate scenario), groundwater temperatures in the Minneapolis/St. Paul metropolitan area could therefore rise by up to 4 °C. Compounding a land use change from “undeveloped” to “fully urbanized” and a 2 × CO₂ climate scenario, groundwater temperatures are projected to rise by about 5 °C at the latitude of Minneapolis/St. Paul.     

Groundwater transport of crater-lake brine at Poa´s Volcano, Costa Rica

Poa´s Volcano is an active stratovolcano in Costa Rica that has a lake in its active crater. The crater lake has high temperatures (50–90 °C), high acidity (pH ≈ 0.0), and a high dissolved-solids content (100 g/kg). The volcano has numerous freshwater springs on its flanks, but a few on the northwestern flank are highly acidic (pH = 1.6–2.5) and have high dissolved-solids concentrations (2–22 g/kg). This study analyzes the regional groundwater system at Poa´s and demonstrates the likelihood that the water discharging from the acidic springs in the Rio Agrio watershed originates at the acidic crater lake. Both heat and solute transport are analyzed on a regional scale through numerical simulations using the HST3D finite-difference model, which solves the coupled equations for fluid flow, heat transport, and solute transport. The code allows fluid viscosity and density to be functions of both temperature and solute concentration. The simulations use estimates for recharge to the mountain and a range of values and various distributions of permeability and porosity. Several sensitivity analyses are performed to test how the uncertainty in many of the model parameters affects the simulation results. These uncertainties yield an estimated range of travel times from the crater lake to the Rio Agrio springs of 1–30 years, which is in close agreement with the results of tritium analyses of the springs. Calculated groundwater fluxes into and out of the crater lake are both about several hundred kg/s. These fluxes must be accounted for in water budgets of the crater lake.     

Geographic variations in shell growth rates of the mussel Diplodon chilensis from temperate lakes of Chile: Implications for biodiversity conservation

The Chilean lake district includes diverse lentic ecosystems along ca. 700 km of the country (36°–43°S), including the “Nahuelbutan lakes”, “Araucanian lakes” and “Chiloe lakes”. This area is recognized as an important “hot spot” of benthic freshwater biodiversity in Southern South America. In Chilean temperate lakes, increased nutrient loads of P and N caused eutrophication, particularly in the Nahuelbutan Lakes. The freshwater Hyriidae mussel Diplodon chilensis (Gray, 1828) which is one of the most abundant species in Chilean temperate lakes, is known to be very susceptible to eutrophication. This species presents a clear reduction in its geographic ranges and is considered to be a threatened species in many Chilean lakes. In this study, we used a correlative approach to determine how eutrophication-driven changes in the food supply and in geographical parameters of different Chilean lakes affected the shell growth rates of D. chilensis. The results obtained from sclerochronological analyses of the mussel shells suggest an association with a group of environmental variables, including geographical types (negative), such as latitude and altitude, and limnological types (positive), especially phosphorous and turbidity. However, the D. chilensis populations under extreme conditions of turbidity in eutrophic and hypertrophic lakes are extinct or nearly so. The high positive correlation of the mean D. chilensis growth rates with orthophosphate (R=0.76; P<0.05), in relation to dissolved inorganic nitrogen, suggests that P is the major limiting factor of the primary productivity in Chilean temperate lakes. We discuss some implications of our results in terms of the conservation of biodiversity in temperate lake ecosystems at different taxonomic levels.     

Geochemistry of the magmatic–hydrothermal system of Kawah Ijen volcano, East Java, Indonesia

Samples from Kawah Ijen crater lake, spring and fumarole discharges were collected between 1990 and 1996 for chemical and isotopic analysis. An extremely low pH (<0.3) lake contains SO₄–Cl waters produced during absorption of magmatic volatiles into shallow ground water. The acidic waters dissolve the rock isochemically to produce “immature” solutions. The strong D and ¹⁸O enrichment of the lake is mainly due to enhanced evaporation at elevated temperature, but involvement of a magmatic component with heavy isotopic ratios also modifies the lake D and ¹⁸O content. The large Δ_SO4–S⁰ (23.8–26.4‰) measured in the lake suggest that dissolved SO₄ forms during disproportionation of magmatic SO₂ in the hydrothermal conduit at temperatures of 250280°C. The lake δ¹⁸O_SO4 and δ¹⁸O_H2O values may reflect equilibration during subsurface circulation of the water at temperatures near 150°C. Significant variations in the lake's bulk composition from 1990 to 1996 were not detected. However, we interpret a change in the distribution and concentration of polythionate species in 1996 as a result of increased SO₂-rich gas input to the lake system.Thermal springs at Kawah Ijen consist of acidic SO₄–Cl waters on the lakeshore and neutral pH HCO₃–SO₄–Cl–Na waters in Blawan village, 17 km from the crater. The cation contents of these discharges are diluted compared to the crater lake but still do not represent equilibrium with the rock. The SO₄/Cl ratios and water and sulfur isotopic compositions support the idea that these springs are mixtures of summit acidic SO₄–Cl water and ground water.The lakeshore fumarole discharges (T=170245°C) have both a magmatic and a hydrothermal component and are supersaturated with respect to elemental sulfur. The apparent equilibrium temperature of the gas is 260°C. The proportions of the oxidized, SO₂-dominated magmatic vapor and of the reduced, H₂S-dominated hydrothermal vapor in the fumaroles varied between 1979 and 1996. This may be the result of interaction of SO₂-bearing magmatic vapors with the summit acidic hydrothermal reservoir. This idea is supported by the lower H₂S/SO₂ ratio deduced for the gas producing the SO₄–Cl reservoir feeding the lake compared with that observed in the subaerial gas discharges. The condensing gas may have equilibrated in a liquid–vapor zone at about 350°C.Elemental sulfur occurs in the crater lake environment as banded sediments exposed on the lakeshore and as a subaqueous molten body on the crater floor. The sediments were precipitated in the past during inorganic oxidation of H₂S in the lake water. This process was not continuous, but was interrupted by periods of massive silica (poorly crystallized) precipitation, similar to the present-day lake conditions. We suggest that the factor controlling the type of deposition is related to whether H₂S- or silica-rich volcanic discharges enter the lake. This could depend on the efficiency with which the lake water circulates in the hydrothermal cell beneath the crater. Quenched liquid sulfur products show δ³⁴S values similar to those found in the banded deposits, suggesting that the subaqueous molten body simply consists of melted sediments previously accumulated at the lake bottom.     

Structure of the Nemrut caldera (Eastern Anatolia, Turkey) and associated hydrothermal fluid circulation

Plio-Quaternary volcanism played an important role in the present physical state of Eastern Anatolia. Mount Nemrut, situated to the west of Lake Van is one of the main volcanic centers in the region, with a spectacular summit caldera 8.5 × 7 km in diameter. The most recent eruptions of the volcano were in 1441, 1597 and 1692. Nemrut Lake covers the western half of the caldera; it is a deep, half-bowl-shaped lake with a maximum depth of 176 m. Numerous eruption centers are exposed within the caldera as a consequence of magma–water interaction. Current activity of Nemrut caldera is revealed as hot springs, fumaroles and a small, hot lake.Self-potential and bathymetric surveys carried out in the caldera were used to characterize the structure of the caldera and the associated hydrothermal fluid circulation. In addition, analyses based on digital elevation models and satellite imagery were used to improve our knowledge about the structure of the caldera. According to SP results, the flanks of the volcano represent “the hydrogeologic zone”, whereas the intra-caldera region is an “active hydrothermal area” where the fluid circulation is controlled by structural discontinuities. There is also a northern fissure zone which exhibits hydrothermal signatures. Nemrut caldera collapsed piecemeal, with three main blocks. Stress controlling the collapse mechanism seems to be highly affected by the regional neotectonic regime. In addition to the historical activity, current hydrothermal and hydrogeologic conditions in the caldera, in which there is a large lake and shallow water table, increase the risk of the quiescent volcano.     

Ground-based thermal imaging of lava lakes at Erebus volcano,Antarctica

Mount Erebus, a large intraplate stratovolcano dominating Ross Island, Antarctica, hosts the world's only active phonolite lava lakes. The main manifestation of activity at Erebus volcano in December 2004 was as the presence of two convecting lava lakes within an inner crater. The long-lived Ray Lake, ~ 1400 m² in area, was the site of up to 10 small Strombolian eruptions per day. A new but short-lived, ~ 1000–1200 m² lake formed at Werner vent in December 2004 sourced by lava flowing from a crater formed in 1993 by a phreatic eruption. We measured the radiative heat flux from the two lakes in December 2004 using a compact infrared (IR) imaging camera. Daily thermal IR surveys from the Main Crater rim provide images of the lava lake surface temperatures and identify sites of upwelling and downwelling. The radiative heat outputs calculated for the Ray and Werner Lakes are 30–35 MW and 20 MW, respectively. We estimate that the magma flux needed to sustain the combined heat loss is ~ 250–710 kg s^− 1, that the minimum volume of the magma reservoir is 2 km³, and that the radius of the conduit feeding the Ray lake is ~ 2 m.     

Using Crater Lake chemistry to predict volcanic activity at Poás Volcano,Costa Rica

Monitoring of crater lake chemistry during the recent decline and disappearance of the crater lake of Poás Volcano revealed that large variations in SO₄/Cl, F/Cl, and Mg/Cl ratios were caused by the enhanced release of HCl vapor from the lake surface due to increasing lake temperature and solution acidity. Variation in the concentration of polythionic acids (H₂S_xO₆, x=4–6) was the most reliable predictor of renewed phreatic eruptive activity at the volcano, exhibiting sharp decreases three months prior to the initiation of phreatic eruptions in June 1987. Polythionic acids may offer a direct indicator of changing subsurface magmatic activity whereas chloride-based element ratios may be influenced by surface volatilization of HCl and subsequent recycling of acidic fluids in crater lake volcanoes.     

A telemetering system for monitoring aqueous polythionates in an active crater lake

We have developed an in-situ monitoring system for aqueous polythionates in the Yugama crater lake at Kusatsu-Shirane volcano. A commercially available nitrate-selective electrode has a Nernstian response to the concentration of polythionates in solution in the range of pH 2.5–11. However, many active crater lakes contain highly acidic water (pH ranges from 0.2 to 2). The monitoring system reduces the acidity of the lake water to pH 4 without changing the concentration of polythionates using a high efficiency cation-exchange membrane technique. The monitoring system was placed in a submersed station, because the Yugama surface lake water freezes in mid-winter. The system is described and preliminary data are reported.     

Bathymetric and geochemical investigation of Kawah Ijen Crater Lake, East Java, Indonesia

A bathymetric survey of Kawah Ijen crater lake was conducted by acoustic sounding in 1996 to compare the lake morphology with those measured in 1922, 1925 and 1938, and to calculate the present lake volume. Even though the lake experienced several hydrothermal eruptions, the maximum depth became shallower (182 m) than before (200 m), resulting in a reduced lake volume (3.0×10⁷ m³).Fifty-two major and minor constituents including rare earth elements and polythionates (PT) of the lake waters at various depths were determined by ICP-AES, ICP-MS and HPLC, respectively. These ions except for several volatile elements are taken up by lake fringe through congruent dissolution of pyroclastics of Kawah Ijen volcano. Most ions are homogeneously distributed throughout the lake, although PT showed a considerable vertical variation. Rare earth elements (REE) in the Kawah Ijen water as well as those from other hyper-acidic crater lakes show distribution patterns likely due to the three rock dissolution (preferential, congruent and residual) types, and their logarithmic concentrations linearly depend upon the pH values of the lake waters.Using the PT degradation kinetics data, production rates of PT, injection rates of SO₂ and H₂S into the lake were estimated to be 114, 86 and 30 tons/day, respectively. Also travel time of the spring water at the Banyupahit Riverhead from Kawah Ijen was estimated to be 600–1000 days through the consideration of decreasing rates of PT. Molten sulfur stocks containing Sn, Cu, Bi sulfides and Pb-barite exposed on the inner crater slope were presumed to be extinct molten sulfur pools at the former lake bottom. This was strongly supported by the barite precipitation temperature estimated through the consideration of the temperature dependence of Pb-chlorocomplex formation.     

Investigating crater lake warming using ASTER thermal imagery: Case studies at Ruapehu,Poás,Kawah Ijen,and Copahué Volcanoes

A two-channel or split-window algorithm designed to correct for atmospheric conditions was applied to thermal images taken by the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) of Lake Yugama on Kusatsu–Shirane volcano in Japan in order to measure the temperature of its crater lake. These temperature calculations were validated using lake water temperatures that were collected on the ground. Overall, the agreement between the temperatures calculated using the split-window method and ground truth is quite good, typically ± 1.5 °C for cloud-free images. Data from fieldwork undertaken in the summer of 2004 at Kusatsu–Shirane allow a comparison of ground-truth data with the radiant temperatures measured using ASTER imagery. Further images were analyzed of Ruapehu, Poás, Kawah Ijen, and Copahué volcanoes to acquire time-series of lake temperatures. A total of 64 images of these 4 volcanoes covering a wide range of geographical locations and climates were analyzed. Results of the split-window algorithm applied to ASTER images are reliable for monitoring thermal changes in active volcanic lakes. These temperature data, when considered in conjunction with traditional volcano monitoring techniques, lead to a better understanding of whether and how thermal changes in crater lakes aid in eruption forecasting.     

Long term AVHRR observations of surface radiative flux from El Chichón crater lake (1996–2006)

It has been shown that due to the small surface of crater lakes, temperature surveillance is a problem using meteorological satellites. This is particularly true for El Chichón surface lake because it's about one tenth of an AVHRR pixel at nadir. In order to guarantee at least one unmixed pixel in AVHRR data, it is necessary to use only AVHRR data from NOAA satellite passes as close as possible to the nadir for the period 1996–2006, therefore AVHRR data of El Chichón's crater lake were only used it they were cloudless and had scan angles close to nadir. The analysis of the time series data shows that lake surface temperature had annual maximum values (> 35 °C) during 1996 and 1997 then surface temperature decay with a negative exponential trend reaching a steady state of about 30 °C in the last years (2004–2006). A seasonal temperature variation between the dry (December to May) and the wet (June to November) seasons is also observed. Differences between nocturnal and midday temperatures indicate the influence of lake energy emission (including reflectance) at midday under a strong short-wave solar radiation. Water surface radiative flux under these conditions reaches an average of 77.8 W m^− 2 and a maximum of 187.1 W m^− 2. Whereas nocturnal heat output from El Chichón crater lake has an average surface radiative flux of 20.4 W m^− 2 and a maximum of 74.3 W m^− 2.     

Discovery of a recurrent lava lake on Saunders Island (South Sandwich Islands) using AVHRR imagery

Monitoring of the remote South Sandwich Islands volcanic arc, using advanced very high resolution radiometer (AVHRR) data, has identified a radiant pixel on channels 3 (3.55–3.93 μm) and (rarely) 4 (10.3–11.3 μm). The position of the pixel is coincident with Mount Michael, Saunders Island, the active summit crater of a snow-covered basaltic stratovolcano with a persistent steam plume. The radiant pixel was continuously present in successive AVHRR images acquired during intervals of several months in the period examined (March 1995–February 1998), although apparently disappearing for similar time intervals. More than 5000 images were examined during this study. The radiant pixel is interpreted to indicate that the crater has contained a lava lake, the first to be recorded in the South Sandwich Islands. The lake appears to persist in the crater for several months at a time, at least. It may have drained completely at times and was probably absent when the crater was viewed briefly during an overflight in January 1997. Persistent or recurring lava lakes are very uncommon world-wide and that at Mount Michael is one of only two recorded in the Antarctic region.     

Sulfur isotopic effects in the disproportionation reaction of sulfur dioxide in hydrothermal fluids: implications for the δS variations of dissolved bisulfate and elemental sulfur from active crater lakes

Sulfur isotope effects during the SO₂ disproportionation reaction to form elemental sulfur (3SO₂+3H₂O→2HSO₄⁻+S+2H⁺) at 200–330°C and saturated water vapor pressures were experimentally determined. Initially, a large kinetic isotopic fractionation takes place between HSO₄⁻ and S, followed by a slow approach to equilibrium. The equilibrium fractionation factors, estimated from the longest run results, are expressed by 1000 ln α_HSO4⁻–S=6.21×10⁶/T²+3.62. The rates at which the initial kinetic fractionation factors approach the equilibrium ones were evaluated at the experimental conditions.δ³⁴S values of HSO₄⁻ and elemental sulfur were examined for active crater lakes including Noboribetsu and Niseko, (Hokkaido, Japan), Khloridnoe, Bannoe and Maly Semiachik (Kamchatka), Poás (Costa Rica), Ruapehu (New Zealand) and Kawah Ijen and Keli Mutu (Indonesia). Δ_HSO4⁻–S values are 28‰ for Keli Mutu, 26‰ for Kawah Ijen, 24‰ for Ruapehu, 23‰ for Poás, 22‰ for Maly Semiachik, 21‰ for Yugama, 13‰ for Bannoe, 9‰ for Niseko, 4‰ for Khloridonoe, and 0‰ for Noboribetsu, in the decreasing order. The SO₂ disproportionation reaction in the magmatic hydrothermal system below crater lakes where magmatic gases condense is responsible for high Δ_HSO4⁻–S values, whereas contribution of HSO₄⁻ produced through bacterial oxidation of reduced sulfur becomes progressively dominant for lakes with lower Δ_HSO4⁻–S values. Currently, Noboribetsu crater lake contains no HSO₄⁻ of magmatic origin. A 40-year period observation of δ³⁴S_HSO4⁻ and δ³⁴S_S values at Yugama indicated that the isotopic variations reflect changes in the supply rate of SO₂ to the magmatic hydrothermal system. This implies a possibility of volcano monitoring by continuous observation of δ³⁴S_HSO4⁻ values. The δ¹⁸O values of HSO₄⁻ and lake water from the studied lakes covary, indicating oxygen isotopic equilibration between them. The covariance gives strong evidence that lake water circulates through the sublimnic zone at temperatures of 140±30°C.     

Exploration drilling and reservoir model of the Platanares geothermal system, Honduras, Central America

Results of drilling, logging, and testing of three exploration core holes, combined with results of geologic and hydrogeochemical investigations, have been used to present a reservoir model of the Platanares geothermal system, Honduras. Geothermal fluids circulate at depths ≥ 1.5 km in a region of active tectonism devoid of Quaternary volcanism. Large, artesian water entries of 160 to 165°C geothermal fluid in two core holes at 625 to 644 m and 460 to 635 m depth have maximum flow rates of roughly 355 and 560 l/min, respectively, which are equivalent to power outputs of about 3.1 and 5.1 MW(thermal). Dilute, alkali-chloride reservoir fluids (TDS ≤ 1200 mg/kg) are produced from fractured Miocene andesite and Cretaceous to Eocene redbeds that are hydrothermally altered. Fracture permeabillity in producing horizons is locally greater than 1500 and bulk porosity is ≤ 6%. A simple, fracture-dominated, volume-impedance model assuming turbulent flow indicates that the calculated reservoir storage capacity of each flowing hole is approximately 9.7 × 10⁶ l/(kg cm⁻²), Tritium data indicate a mean residence time of 450 yr for water in the reservoir. Multiplying the natural fluid discharge rate by the mean residence time gives an estimated water volume of the Platanares system of ≥ 0.78 km³. Downward continuation of a 139°C/km “conductive” gradient at a depth of 400 m in a third core hole implies that the depth to a 225°C source reservoir (predicted from chemical geothermometers) is at least 1.5 km. Uranium-thorium disequilibrium ages on calcite veins at the surface and in the core holes indicate that the present Platanares hydrothermal system has been active for the last 0.25 m.y.     

Water chemistry of Lake Quilotoa (Ecuador) and assessment of natural hazards

A geochemical survey carried out in November 1993 revealed that Lake Quilotoa was composed by a thin (14 m) oxic epilimnion overlying a 200 m-thick anoxic hypolimnion. Dissolved CO₂ concentrations reached 1000 mg/kg in the lower stratum. Loss of CO₂ from epilimnetic waters, followed by calcite precipitation and a consequent lowering in density, was the apparent cause of the stratification.The Cl, SO₄ and HCO₃ contents of Lake Quilotoa are intermediate between those of acid–SO₄–Cl Crater lakes and those of neutral-HCO₃ Crater lakes, indicating that Lake Quilotoa has a ‘memory’ of the inflow and absorption of HC1- and S-bearing volcanic (magmatic) gases. The Mg/Ca ratios of the lake waters are governed by dissolution of local volcanic rocks or magmas, but K/Na ratios were likely modified by precipitation of alunite, a typical mineral in acid–SO₄–Cl Crater lakes.The constant concentrations of several conservative chemical species from lake surface to lake bottom suggest that physical, chemical and biological processes did not have enough time, after the last overturn, to cause significant changes in the contents of these chemical species. This lapse of time might be relatively large, but it cannot be established on the basis of available data. Besides, the lake may not be close to steady state. Mixing of Lake Quilotoa waters could presently be triggered by either cooling epilimnetic waters by 4°C or providing heat to hypolimnetic waters or by seismic activity.Although Quilotoa lake contains a huge amount of dissolved CO₂ (3×10¹¹ g), at present the risk of a dangerous limnic eruption seems to be nil even though some gas exsolution might occur if deep lake waters were brought to the surface. Carbon dioxide could build up to higher levels in deep waters than at present without any volcanic re-awakening, due to either a large inflow of relatively cool CO₂-rich gases, or possibly a long interval between overturns. Periodical geochemical surveys of Lake Quilotoa are, therefore, recommended.     

Fumarole monitoring with a handheld infrared camera: Volcán de Colima,Mexico, 2006–2007

This paper describes a methodology for the monitoring of fumarole temperatures at medium ranges (~ 6 km) using a handheld infrared camera (wavelength range: 8–13.5 µm). As a relationship between fumarole temperatures, gas flux and volcanic activity has been demonstrated by a number of studies, fumarole temperature data has a potential use as a monitoring tool. Volcán de Colima is an andesitic stratovolcano with a 300 m diameter summit crater formed by the destruction of the 2004 lava dome by a series of explosions in 2005. Between January 2006 and August 2007, sequences of thermal images were recorded from a viewpoint 6 km to the north during regular 24–48 hour monitoring excursions. The temperatures of fumaroles on the crater rim and the ground surface on the volcano's flanks were measured. A methodology was developed to remove data affected by clouds or volcanic water vapour based on rates of temperature change and scatter within the data. For the remaining data, it is demonstrated mathematically that at this range, typical variations in atmospheric transmissivity will affect the apparent temperatures by +/− 2 °C, while a 25% change in fumarole heat flux would change it by 5–10 °C. The mean night-time apparent temperature of the fumaroles was calculated for each excursion and showed an irregular decline over the 19 month period. Subtracting the radiant heat flux of flank rocks from those of the fumaroles removes seasonal variations and gives the clearest view of trends in the fumarole heat flux. A sharp drop in fumarole temperature during February 2007 coincided with the emergence of a lava dome in the crater. The declining fumarole temperature is interpreted to reflect decreasing gas flux from the crater in line with a change in eruptive regime from frequent, small, ash-rich explosions to slow effusion of lava.     

Geological evolution of a pleistocene rhyolitic center — Sierra La Primavera, Jalisco, México

The Sierra La Primavera volcanic complex consists of late Pleistocene comenditic lava flows and domes. ash-flow tuff, air-fall pumice, and cold caldera-lake sediments. The earliest lavas were erupted about 120,000 years ago, and were followed approximately 95,000 years ago by the eruption of about 20 km³ of magma as ash flows that form the compositionally-zoned Tala Tuff. Collapse of the roof zone of the magma chamber led to the formation of a shallow 11-km-diameter caldera. It soon filled with water, forming a caldera lake in which sediment began to collect. At about the same time, two central domes erupted through the middle of the lake and a “giant pumice horizon”, an important stratigraphic marker, was deposited. Shortly thereafter ring domes erupted along two parallel arcs: one along the northeast portion of the ring fracture, and the other crossing the middle of the lake. All these events occurred during a period of approximately 5,000–10,000 years. Sedimentation continued and a period of volcanic quiescence was marked by the deposition of some 30 m of fine-grained ashy sediments virtually free from pumice lapilli. Approximately 75,000 years ago, a new group of ring domes erupted at the southern margin of the lake. These domes are lapped by only 10–20 m of sediments, as uplift resulting from renewed insurgence of magma brought an end to the lake. This uplift culminated in the eruption, beginning approximately 60,000 years ago, of aphyric lavas along a southern arc. The youngest of these lavas erupted approximately 20,000–30,000 years ago.The four major fault systems in the Sierra La Primavera are related to caldera collapse or to uplift caused by the insurgence of the southern are magma. Steam vents and larga-discharge 65°C hot springs are associated with the faulting. Calculated equilibrium temperatures of the geothermal fluids are 170°C, but temperatures in excess of 240°C have been encountered in an exploratory drill hole.A seismic survey showed attenuation of both S and P waves within the caldera, P waves attenuated more severely than S waves. The greatest attenuation is associated with an area of steam vents, and the rapid lateral variations in attenuation suggest that they are produced by a shallow geothermal system rather than by underlying magma.     

Transport ofBe andBe in the ocean

Beryllium isotopes (¹⁰Be and⁹Be) have been measured in suspended particles of < 1 mm size collected by mid-water sediment traps deployed in the eastern Pacific at MANOP sites H (6°32′N, 92°50′W, water depth 3600 m) and M (8°50′N, 104°00′W, 3100 m). For comparison, surface sediments from box cores taken from the two sites were also studied. The concentrations of¹⁰Be and⁹Be in sediment-trap particles are about an order of magnitude smaller than those in the bottom sediments which contain about 8 × 10⁹ and 6 × 10¹⁶ atoms g⁻¹ of¹⁰Be and⁹Be, respectively. The sediment trap samples collected from 50 m off the bottom showed significant (26–63%) contributions from resuspended bottom sediments. The¹⁰Be/⁹Be ratio in trap samples varies from 3 to 20 × 10⁻⁸. The variation may partly result from varied proportion of authigenic/detrital material. The fluxes of both isotopes exhibit a very strong seasonality. The fluxes of¹⁰Be into the traps at about 1500 m are estimated as 9 × 10⁵ and 4 × 10⁵ atoms cm⁻² a⁻¹ at sites H and M respectively. These values are to be compared with the fluxes into the sediments of 4–5 × 10⁵ atoms cm⁻² a⁻¹ at both locations. Good correlations exist between¹⁰Be,⁹Be and²⁷Al indicating that the primary carrier phase(s) for the beryllium isotopes in the water column may be aluminosilicates.     

Development and application of a remote sensing-based salinity prediction model for a large estuarine lake in the US Gulf of Mexico coast

Salinity in estuaries is highly variable due to river discharge, tidal motion, and winds. Information on the spatial and temporal changes in salinity can provide important ecological indications, but accurate monitoring of the space–time variability for a large estuary is often costly and time-consuming. This study applied remote sensing techniques to develop a salinity prediction model for Lake Pontchartrain, a large estuarine lake located in the Northern Gulf of Mexico, USA. “Ground truth” salinity was measured along two transects across the lake and near the shoreline. Water-leaving reflectance from the measurement locations was extracted from Landsat Thematic Mapper (TM) images pre-processed through “banding” noise reduction and radiometrical correction approaches. Ordinary least square and ridge regression methods were performed to identify model parameters and to determine relationships between salinity and reflectance. Salinity in the lake on eight dates was predicted with the developed model. Difference in salinity level and patterns, and impacts of Hurricanes Katrina and Rita on salinity were assessed with ANOVA and Fuzzy Similarity methods. The results showed that the model achieved a high power in prediction of the lake salinity (R² = 0.89 and RMSE of validation = 0.27). Reflectance from TM bands 1, 2, and 4 was positively correlated to salinity levels and explained 1.9%, 20.3%, and 10.2% variance in salinity levels. Reflectance from bands 3 and 5 was negatively correlated to salinity and explained 34.1% and 31.2% variance. Under normal circumstances without the impacts of hurricanes, the lake salinity presented two patterns with average salinity level of 5.5 ppt. After Katrina’s landfall, the average was significantly increased by 1.1 ppt and the spatial patterns were altered. The pattern on 30 August 2005 was the most dissimilar one as compared to the two normal patterns, and then followed by the patterns on 9 and 25 October, and 7 September 2005. This study demonstrated that satellite remote sensing techniques can be applied to monitor salinity in coastal environments, and that freshwater discharge not only affects salinity levels and patterns under normal conditions but also is crucial for the return of salinity patterns to normal conditions after hurricane disturbance.     

Chemical geothermometry and fluid–mineral equilibria for the Ömer–Gecek thermal waters, Afyon area, Turkey

Thermal waters of the Ömer–Gecek geothermal field, Turkey, with temperatures ranging from 32 to 92°C vary in chemical composition and TDS contents. They are generally enriched in Na–Cl–HCO₃ and suggest deep water circulation. Silica and cation geothermometers applied to the Ömer–Gecek thermal waters yield reservoir temperatures of 75–155°C. The enthalpy–chloride mixing model, which approximates a reservoir temperature of 125°C for the Ömer–Gecek field, accounts for the diversity in the chemical composition and temperature of the waters by a combination of processes including boiling and conductive cooling of deep thermal water and mixing of the deep thermal water with cold water. It is also determined that the solubility of silica in most of the waters is controlled by the chalcedony phase. Equilibrium states of the Ömer–Gecek thermal waters studied by means of the Na–K–Mg triangular diagram, Na–K–Mg–Ca diagram, K–Mg–Ca geoindicator diagram, activity diagrams in the systems composed of Na₂O–CaO–K₂O–Al₂O₃–SiO₂–CO₂–H₂O phases, log SI diagrams, and finally the alteration mineralogy indicate that most of the spring and low-temperature well waters in the area can be classified as shallow or mixed waters which are likely to be equilibrated with calcite, chalcedony and kaolinite at predicted temperature ranges similar to those calculated from the chemical geothermometers. It was also observed that mineral equilibrium in the Ömer–Gecek waters is largely controlled by CO₂ concentrations.