Volcanic earthquakes, and their relationship to eruptions at Ruapehu and Ngauruhoe volcanoes

Crustal earthquakes near Ruapehu and Ngauruhoe fall into two classes, each of which can be subdivided. On the one hand, there are high-frequency events ( 3 Hz) with sharp, well-defined phases, mainly concentrated beneath Ruapehu Crater Lake. Low-frequency events (< 2 Hz), on the other hand, are common at shallower depths under both volcanoes. These are usually emergent multiple events, and are often closely associated with eruptions.The low-frequency events resemble Minakami's B-type and explosion earthquakes, but sometimes occur where no vent exists and rather deeper than his formal definition (< 1 km) permits. More importantly, they lack reliable criteria (wave-form or magnitude differences) to distinguish between his two groups. Whether or not they accompany an eruption (Minakami's definition of explosion earthquake) appears to depend on whether the volcanoes are in a “closed-” or “open-vent” condition. The high-frequency earthquakes are similar in wave-form to Minakami's A-type. However, many at Ruapehu (here designated “roof-rock” earthquakes) originate at shallower depths than the B-type earthquakes, which is contrary to Minakami's definition.Difficulty in applying Minakami's classification rigorously, and the fact that low frequencies may be due to abnormal attenuation of higher frequencies along the path, rather than to their suppression or absence at the source, has led to reclassification of earthquakes near the volcanoes into two broad groups, tectonic and volcanic. The former includes all high-frequency earthquakes, and those discrete events in which dominant low frequencies are due to path effects. The latter includes multiple and emergent events which show evidence of prolonged or repetitive source mechanism. Dominant low frequencies are ascribed to occurrence in heat-weakened material, and high frequencies to instantaneous source mechanisms operating in competent rock. The term volcano-tectonic describes tectonic earthquakes within some arbitrary distance of a volcano.At Ngauruhoe and Ruapehu, volcanic earthquakes accompany explosive, vent-clearing eruptions. Subsequent “open-vent” degassing and ash emission, however, although often powerful and prolonged, usually occurs without earthquakes. Such activity is, however, frequently accompanied by volcanic tremor. At Ruapehu, under “closed-vent” conditions, when lake temperature is low, low-frequency earthquakes up to magnitude M_L = 3.4 have occurred without any eruption.Five types of phreatic eruptions are identified at Ruapehu, each having a distinctive seismic pattern. The three most explosive types appear to be generated by a chain reaction process, and all involve flashing of water to steam; the first by failure of the roof, with little precursory seismicity, after a “closed-vent” period, during which lake temperature decreases; the second, after prolonged heating of the lake and much preliminary volcanic tremor, interpreted as due to rising magma; and the third, under “open-vent” conditions in the wake of one of the two preceding types. A fourth probably occurs in wet sediments near the base of the lake, as a result of upward migration of hot gas, and a fifth, aseismic, or accompanied by very weak volcanic tremor, is associated with convective overturn within Crater Lake.     

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.     

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.     

Food web manipulation by extreme enhancement of piscivory: an invertebrate predator compensates for the effects of planktivorous fish on a plankton community

By enhancing the stock of piscivorous fish in a whole-lake experiment in Gräfenhain (Germany) since 1981 to such an extent that almost all planktivorous fish were exterminated, we examined the concept of “over-biomanipulation”. This hypothesis predicts that (a) extremely strong piscivory will allow uncontrolled development of large invertebrate predators and (b) these invertebrates can exert the same strong predation pressure on large herbivorous zooplankton as planktivorous fish. The hypothesis is tested the first time by a cross-comparison of the long-term response of the plankton community structure in the experimental lake (Piscivore L.) with (1) the intermediate response in the same lake and with (2) that of the long-term state in a nearby reference lake (Planktivore L.) densely inhabited by planktivorous fish (Leucaspius delineatus, a small cyprinid). The intermediate (1989–1992) response in Piscivore L. revealed a strong increase of the abundance of the invertebrate predator Chaoborus flavicans. Large daphnids were able to coexist with C. flavicans so that edible phytoplankton were suppressed and water transparency increased.As part (a) of the hypothesis predicts, the long-term response in Piscivore L. was characterized by the immigration of the larger predator Chaoborus obscuripes which displaced the smaller C. flavicans completely. The results support also part (b) of the hypothesis of “over-biomanipulation”: C. obscuripes-dominated Piscivore L. showed not much difference in biomasses of daphnids and total and edible phytoplankton as well as Secchi depth compared with fish-dominated Planktivore Lake. On the other hand, C. obscuripes-dominated Piscivore L. was characterized by distinctly lower biomass of daphnids, mean body volume of all crustaceans and Secchi depth as well as by higher biomass of edible phytoplankton compared with C. flavicans-dominated Piscivore Lake. We conclude that long-lasting success of biomanipulation cannot be achieved by extremely high piscivory leading to the almost complete extermination of planktivorous fish.     

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.     

Basic limnological characteristics of the shallow eutrophic lake Grimnitzsee (Brandenburg, Germany)

The medium shallow lake Grimnitzsee (maximum depth: 9.9 m; mean depth: 4.6 m; area: 7.7 · 10⁶ m²) which is situated in the biosphere reserve “Schorfheide-Chorin” in northern Brandenburg (Germany) was studied in 1994 and 1995. A bathymetric map of Grimnitzsee is given for the first time. The lake is usually polymictic although in 1994 and 1995 relatively long summer stratification was observed due to very high global radiation input. Nutrient concentration, light climate, oxygen status, phytoplankton biomass and the species composition of littoral diatoms characterize the lake as eutrophic. Special features deducible from the lake's polymictic character were the multiple development of aerobic or anaerobic strata above the sediment, the fast recovery of silicon concentration in the water column after diatom sedimentation, the importance of resuspension for the success of planktonic diatom populations, and an only moderate correlation between chlorophyll a concentration and light attenuation as well as seston dry weight probably due to the influence of suspended particles.     

Infrared surveillance of crater lakes using satellite data

Under optimum circumstances, thermal infrared data recorded from satellites can measure water surface temperatures to accuracies of a few tenths of a degree Celsius. Such techniques are applied here to evaluate volcanic crater lake temperatures. At present, band 6 of the Landsat Thematic Mapper (TM) is the most pertinent sensor in this respect, although its nominal 120 × 120 m “footprint” only permits useful measurements of circular lakes exceeding 340 m in diameter. In addition, the radiative properties of the atmosphere between sensor and target at the instant of observation should be well-characterised in order to make confident measurements of surface temperatures with single-band infrared data.An analysis is presented of three TM band-6 images of the crater lake at Poás volcano, Costa Rica, recorded on February 6, 1986, March 13, 1987, and May 10, 1988. In the February 1986 scene, the band-6-derived water surface temperature is 36°C which is 3°C less than a field measurement made in the same month. Since the satellite measurement was integrated over some 14,400 m² of the lake surface, while the field measurement was obtained at a single point below the surface, the former may be more representative of surface heat losses. Subsequent TM data reveal an increasing discrepancy between contemporaneous field and satellite observations, probably because the lake diameter had decreased to 250 m by March 1987, and to less than 200 m by mid-1988, greatly reducing the likelihood of obtaining a pure “lake” pixel.The Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) which is earmarked for orbit on the first of NASA's Earth Observing System satellite platforms later this decade has five discrete spectral bands in the thermal infrared region and will produce data composed of 90 × 90 m pixels. These specifications could enable the determination of water surface temperatures of > 250 m diameter crater lakes by algorithms that implicitly account for atmospheric effects.     

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.     

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.     

Interactive effects of environmental variability and human impacts on the long-term dynamics of an Amazonian floodplain lake and a South Atlantic coastal lagoon

Human activities are exposing freshwater ecosystems to a wide range of stressors, whose direct and indirect effects can be alleviated or exacerbated through interactive effects with dynamic environmental drivers. This study used long-term data from two Neotropical lacustrine freshwater systems (Batata Lake, an Amazonian floodplain lake and Imboassica lagoon, an Atlantic coastal lagoon) subjected to different kinds of environmental fluctuations (i.e., flood pulse and sandbar opening) and anthropogenic impacts (i.e., siltation and eutrophication). Our objective was to determine whether the effects of human perturbations are contingent on modifications of important biotic and abiotic characteristics through environmental variability. For both ecosystems, environmental variability consistently interacted with anthropogenic perturbations to alter most of the variables analyzed, such as nutrient dynamics, chlorophyll-a concentration, zooplankton and benthic invertebrate species richness, and temporal community stability, which indicates that interactive effects between environmental variability and anthropogenic perturbations may impact a myriad of ecosystem properties. Furthermore, the nature of these interactive effects was highly dependent on the variable considered and on the ecosystem analyzed. For example, at Imboassica lagoon, sandbar openings interacted synergistically with trophic state to increase the phosphorus concentration in the water column. At Batata Lake, flooding generally alleviated the negative effects of siltation on species richness by both diluting inorganic suspended material concentration and by promoting local recruitment from the regional species pool. Such results indicate that our ability to understand and predict the outcome of anthropogenic impacts on inland aquatic systems can be hampered if we consider human stressors as “static” phenomena disconnected from dynamic interactions with major local environmental drivers.     

Monitoring the water balance of Lake Victoria, East Africa, from space

Using satellite gravimetric and altimetric data, we examine trends in water storage and lake levels of multiple lakes in the Great Rift Valley region of East Africa for the years 2003–2008. GRACE total water storage estimates reveal that water storage declined in much of East Africa, by as much as , while altimetric data show that lake levels in some large lakes dropped by as much as 1–2 m. The largest declines occurred in Lake Victoria, the Earth’s second largest freshwater body. Because the discharge from the outlet of Lake Victoria is used to generate hydroelectric power, the role of human management in the lake’s decline has been questioned. By comparing catchment water storage trends to lake level trends, we confirm that climatic forcing explains only about 50decline. This analysis provides an independent means of assessing the relative impacts of climate and human management on the water balance of Lake Victoria that does not depend on observations of dam discharge, which may not be publically available. In the second part of the study, the individual components of the lake water balance are estimated. Satellite estimates of changes in lake level, precipitation, and evaporation are used with observed lake discharge to develop a parameterization for estimating subsurface inflows due to changes in groundwater storage estimated from satellite gravimetry. At seasonal timescales, this approach provides closure to Lake Victoria’s water balance to within . The third part of this study uses the water balance of a downstream water body, Lake Kyoga, to estimate the outflow from Lake Victoria remotely. Because Lake Kyoga is roughly 20 times smaller in area than Lake Victoria, its water balance is strongly influenced by inflow from Lake Victoria. Lake Kyoga has been shown to act as a linear reservoir, where its outflow is proportional to the height of the lake. This model can be used with satellite altimetric lake levels to estimate a time series of Lake Victoria discharge with an rms error of about .     

Geochemistry of the acid Kawah Putih lake, Patuha Volcano, West Java, Indonesia

Kawah Putih is a summit crater of Patuha volcano, West Java, Indonesia, which contains a shallow, 300 m-wide lake with strongly mineralized acid–sulfate–chloride water. The lake water has a temperature of 26–34°C, pH=<0.5–1.3, S^tot=2500–4600 ppm and Cl=5300–12 600 ppm, and floating sulfur globules with sulfide inclusions are common. Sulfur oxyanion concentrations are unusually high, with S₄O₆²⁻+S₅O₆²⁻+S₆O₆²⁻=2400 – 4200 ppm. Subaerial fumaroles (<93°C) on the lake shore have low molar SO₂/H₂S ratios (<2), which is a favorable condition to produce the observed distribution of sulfur oxyanion species. Sulfur isotope data of dissolved sulfate and native sulfur show a significant ³⁴S fractionation (ΔSO₄–S_e of 20‰), probably the result of SO₂ disproportionation in or below the lake. The lake waters show strong enrichments in ¹⁸O and D relative to local meteoric waters, a result of the combined effects of mixing between isotopically heavy fluids of deep origin and meteoric water, and evaporation-induced fractionation at the lake surface. The stable-isotope systematics combined with energy-balance considerations support very rapid fluid cycling through the lake system. Lake levels and element concentrations show strong seasonal fluctuations, indicative of a short water residence time in the lake as well.Thermodynamic modeling of the lake fluids indicates that the lake water is saturated with silica phases, barite, pyrite and various Pb, Sb, Cu, As, Bi-bearing sulfides when sulfur saturation is assumed. Precipitating phases predicted by the model calculations are consistent with the bulk chemistry of the sulfur-rich bottom sediments and their identified mineral phases. Much of the lake water chemistry can be explained by congruent rock dissolution in combination with preferential enrichments from entering fumarolic gases or brines and element removal by precipitating mineral phases, as indicated by a comparison of the fluids, volcanic rocks and lake bed sediment.Flank springs on the mountain at different elevations vary in composition, and are consistent with local rock dissolution as a dominant factor and pH-dependent element mobility. Discharges of warm sulfate- and chloride-rich water at the highest elevation and a near-neutral spring at lower level may contain a small contribution of crater-lake water. The acid fluid-induced processes at Patuha have led to the accumulation of elements that are commonly associated with volcano-hosted epithermal ore deposits. The dispersal of heavy metals and other potentially toxic elements from the volcano via the local drainage system is a matter of serious environmental concern.     

Stationary magnetospheric convection on November 24, 1981. 2.Small-scale structures in the dayside cusp/cleft

A case study of the dayside cusp/cleft region during an interval of stationary magnetospheric convection (SMC) on November, 24, 1981 is presented, based on detailed measurements made by the AUREOL-3 satellite. Layered small-scale field-aligned current sheets, or loops, superimposed to a narrow V-shaped ion dispersion structure, were observed just equatorward from the region of the “cusp proper”. The equatorward sheet was accompanied by a very intense and short (less than 1 s) ion intensity spike at 100 eV. No major differences were noted of the characteristics of the LLBL, or “boundary cusp”, and plasma mantle precipitation during this SMC period from those typical of the cusp/cleft region for similar IMF conditions. Simultaneous NOAA-6 and NOAA-7 measurements described in Despirak et al. were used to estimate the average extent of the “cusp proper” (defined by dispersed precipitating ions with the energy flux exceeding 10⁻³ erg cm⁻² s⁻¹) during the SMC period, as ≈0.73∼ ILAT width, 2.6–3.4 h in MLT, and thus the recently merged magnetic flux, 0.54–0.70 × 10⁷ Wb. This, together with the average drift velocity across the cusp at the convection throat, ≈0.5 km s⁻¹, allowed to evaluate the cusp merging contribution to the total cross-polar cap potential difference, ≈33.8–43.8 kV. It amounts to a quite significant part of the total cross-polar cap potential difference evaluated from other data. A “shutter” scenario is suggested for the ion beam injection/penetration through the stagnant plasma region in the outer cusp to explain the pulsating nature of the particle injections in the low- and medium-altitude cusp region.     

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.     

Chemical evolution and volcanic activity of the active crater lake of Poás volcano, Costa Rica, 1993–1997

Concentrations of chloride and sulfate and pH in the hot crater lake (Laguna Caliente) at Poás volcano and in acid rain varied over the period 1993–1997. These parameters are related to changes in lake volume and temperature, and changes in summit seismicity and fumarole activity beneath the active crater. During this period, lake level increased from near zero to its highest level since 1953, lake temperature declined from a maximum value of 70°C to a minimum value of 25°C, and pH of the lake water increased from near zero to 1.8. In May 1993 when the lake was nearly dry, chloride and sulfate concentrations in the lake water reached 85,400 and 91,000 mg l⁻¹, respectively. Minimum concentrations of chloride and sulfate after the lake refilled to its maximum volume were 2630 and 4060 mg l⁻¹, respectively. Between January 1993 and May 1995, most fumarolic activity was focused through the bottom of the lake. After May 1995, fumarolic discharge through the bottom of the lake declined and reappeared outside the lake within the main crater area. The appearance of new fumaroles on the composite pyroclastic cone coincided with a dramatic decrease in type B seismicity after January 1996. Between May 1995 and December 1997, enhanced periods of type A seismicity and episodes of harmonic tremor were associated with an increase in the number of fumaroles and the intensity of degassing on the composite pyroclastic cone adjacent to the crater lake. Increases in summit seismic activity (type A, B and harmonic tremor) and in the height of eruption plumes through the lake bottom are associated with a period of enhanced volcanic activity during April–September 1994. At this time, visual observations and remote fumarole temperature measurements suggest an increase in the flux of heat and gases discharged through the bottom of the crater lake, possibly related to renewed magma ascent beneath the active crater. A similar period of enhanced seismic activity that occurred between August 1995 and January 1996, apparently caused fracturing of sealed fumarole conduits beneath the composite pyroclastic cone allowing the focus of fumarolic degassing to migrate from beneath the lake back to the 1953–1955 cone. Changes in the chemistry of summit acid rain are correlated changes in volcanic activity regardless of whether fumaroles are discharging into the lake or are discharging directly into the atmosphere.     

Paleomagnetic analysis of curved thrust belts reproduced by physical models

This paper presents a new methodology for studying the evolution of curved mountain belts by means of paleomagnetic analyses performed on analogue models. Eleven models were designed aimed at reproducing various tectonic settings in thin-skinned tectonics. Our models analyze in particular those features reported in the literature as possible causes for peculiar rotational patterns in the outermost as well as in the more internal fronts. In all the models the sedimentary cover was reproduced by frictional low-cohesion materials (sand and glass micro-beads), which detached either on frictional or on viscous layers. These latter were reproduced in the models by silicone. The sand forming the models has been previously mixed with magnetite-dominated powder. Before deformation, the models were magnetized by means of two permanent magnets generating within each model a quasi-linear magnetic field of intensity variable between 20 and 100 mT. After deformation, the models were cut into closely spaced vertical sections and sampled by means of 1×1-cm Plexiglas cylinders at several locations along curved fronts. Care was taken to collect paleomagnetic samples only within virtually undeformed thrust sheets, avoiding zones affected by pervasive shear. Afterwards, the natural remanent magnetization of these samples was measured, and alternating field demagnetization was used to isolate the principal components. The characteristic components of magnetization isolated were used to estimate the vertical-axis rotations occurring during model deformation. We find that indenters pushing into deforming belts from behind form non-rotational curved outer fronts. The more internal fronts show oroclinal-type rotations of a smaller magnitude than that expected for a perfect orocline. Lateral symmetrical obstacles in the foreland colliding with forward propagating belts produce non-rotational outer curved fronts as well, whereas in between and inside the obstacles a perfect orocline forms only when the ratio between obstacles' distance and thickness of the cover is greater than 10. Finally, when a belt collides with an obstacle in the foreland oblique to the shortening direction the outer front displays rotations opposite in sign to oroclinal-type rotations, whereas the internal fronts seem to assume an “oroclinal type” rotational pattern. Furthermore rotation is easier in laterally unconfined models, i.e. when the wedge can “escape” laterally. The results from our models may be useful when compared to paleomagnetic rotations detected in natural arcs. In these cases, our results may allow for better understanding the tectonic setting controlling the genesis of curved mountain fronts, as is the case of the Gela Nappe of Sicily we compare with some of our models.