Earthquake-induced variations in the composition of the water in the geothermal reservoir at Vulcano Island, Italy

During 1979–1989, variations were observed in the oxygen composition of the water contained in the geothermal reservoir at Vulcano Island, Italy.The reservoir water, that has a magmatic origin, showed an oxygen composition of +1.0±0.5‰ δ¹⁸O during periods without local tectonic earthquakes, and an oxygen composition of +3.4±0.5‰ δ¹⁸O after the highest-energy seismic activity that occurred recently near the island. A slight increase of the δ¹⁸O value in the reservoir water was also observed after a low-energy sequence of tectonic earthquakes that occurred at very shallow depth just beneath Vulcano Island. These ¹⁸O variations in the reservoir water are consistent with earthquake-induced increases in the contribution from high-temperature δ¹⁸O-rich magmatic condensate to the geothermal reservoir, and with subsequent decreases in the δ¹⁸O value due to ¹⁸O exchanges at the temporarily increased reservoir temperature during reactions between the highly reactive magmatic condensate and the local rocks.Only minor changes in the deuterium composition of the reservoir water occurred with time, as the δD value in the magmatic condensate released from the magma after major local earthquakes quickly approached the δD value of the water contained in the geothermal reservoir.Also the chloride concentration in the reservoir water appears to be linked to the contribution from the magmatic fluid. This chloride content seems not to have undergone major changes with time, as it may be buffered by temporary increases in the reservoir temperature up to values >300°C induced by major local earthquakes. This mechanism may possibly occur also in other magmatic–hydrothermal systems.     

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.     

Oxygen and hydrogen isotopes in the volcano-sedimentary complex of Huelva (Iberian Pyrite Belt): example of water circulation through a volcano-sedimentary sequence

The minerals of basic and acidic rocks from the volcano-sedimentary sequence in the Huelva area, Spain, Iberian Pyrite Belt, display an extendedδ¹⁸O enrichment. Quartzδ¹⁸O values from quartz-keratophyres vary from +10.5 to +17.0 and feldsparδ¹⁸O values from +14.4 to +16.0. For the spilite or spilitized doleritesδ¹⁸O values vary from +9.9 to +13.4 for feldspar, from +6.4 to +9.8 for chlorite, from +3.7 to +4.3 for ilmenite and from +13.6 to +14.0 for quartz, but pyroxene exhibits magmatic values, from +5.3 to +6.1 with an exception at +7.5. The chloriteδD values vary from −34 to −43‰.This is attributed to hydrothermal alteration with seawater enriched inδ¹⁸O by circulation through sediments.The temperatures of interaction determined from isotopic fractionations between minerals range from 400° to 520°C.CalculatedδD andδ¹⁸O values for water in equilibrium with the minerals at isotopic temperatures range from −16 to +5 and from +8.3 to +12.8, respectively.A model of circulation of seawater through a pile of sedimentary rocks and then through basaltic rocks is proposed to explain the high¹⁸O compositions of the rocks from the Huelva District. Water/rock mass ratios calculated from this model range between 0.3 and 0.7 for the determined range of temperatures.     

Hydrothermal silica chimney fields in the Galapagos Spreading Center at 86°W

Silica chimneys were discovered in 1985 at 86°W in the rift valley of the Galapagos Spreading Center at 2600 m depth (“Cauliflower Garden”). The inactive chimneys lack any sulfides and consist almost entirely of amorphous silica (up to 96 wt.% SiO₂, opal-A); Fe and Mn oxides are minor constituents. Oxygen isotope data show that formation of the silica chimneys took place at temperatures between 32°C (+29.9‰ δ¹⁸O) and 42°C (+27.8‰ δ¹⁸O).Th/Udating reveals a maximum age of 1440 ± 300y. Amorphous silica solubility relations indicate that the silica chimneys were formed by conductive cooling of pure hydrothermal fluids or by conductive cooling of a fluid/seawater mixture. Assuming equilibrium with quartz at 500 bars, initial fluid temperatures of more than 175°C (i.e., a concentration of > 182 ppm SiO₂) were required to achieve sufficient supersaturation for the deposition of amorphous silica at 40°C and 260 bars. If the silica chimneys originate from the same or a similar fluid as higher-temperature ( < 300°C) sulfide-silica precipitates found nearby (i.e., 2.5 km away), then subsurface deposition of sulfides may have occurred.     

Diffusion effects on oxygen isotope temperatures of slowly cooled igneous and metamorphic rocks

Recently obtained data on oxygen diffusion in feldspars, quartz, and hornblende permit the prediction of the apparent¹⁸O¹⁶O temperatures that would be measured in a rock that consisted only of those three minerals, and cooled slowly from high temperature. The computed temperatures would be based on the differences in the¹⁸O¹⁶O ratios between coexisting pairs of minerals. The present calculation takes into account the diffusion rates for oxygen as a function of temperature, the cooling rate of the rock, the mineral grain sizes, and the mode of the rock. For mineral grains 1 mm in radius, and a cooling rate of 10°C/m.y., the minimum difference in apparent temperature between quartz-feldspar and feldspar-hornblende pairs will be 115°C, despite the assumption of a normal, uneventful, slow cooling history to room temperature. Further, the apparent quartz-hornblende temperature will range over 30°C (590–620°C) depending on the mode of the rock. For a cooling rate of 1000°C/m.y., the apparent difference in temperature can be as much as 400°C. Consequently, consistency in temperatures obtained by oxygen isotope analysis should not be expected in most high-grade metamorphic rocks or igneous rocks which are cooled slowly. Departures from the pattern of temperatures obtained in this model would imply a very rapid quench from high temperature, or a complex history for the rock. For some minerals, including hornblende, the relation between temperature and the equilibrium fractionation of oxygen isotopes between coexisting phases has been derived from observed relations in natural specimens. The choice of the specimens used for such calibrations needs to be re-evaluated in light of these findings. This may result in a change in the equilibrium equation constants.An example from the literature, the San Jose tonalite, Baja California, Mexico, was modelled and yieldsδ¹⁸O concentrations in the minerals that correspond closely with the measured values. This suggests that the model used is appropriate, that the rock has had a simple thermal history, and that it cooled at 100–200°C/m.y. over the temperature range 800–500°C. The set of paleotemperatures obtained for a rock will, in general, yield neither the mineral closure temperatures nor the formation or crystallization temperatures. On the other hand, the cooling rate of the rock may be derived from the data. This, in turn, may have important tectonic implications with regard to denudation and uplift rates.     

Isotope geochemistry of minerals and fluids from Newberry volcano, Oregon

Isotopic compositions were determined for hydrothermal quartz, calcite, and siderite from core samples of the Newberry 2 drill hole, Oregon. The δ¹⁵O values for these minerals decrease with increasing temperatures. The values indicate that these hydrothermal minerals precipitated in isotopic equilibrium with water currently present in the reservoirs. The δ¹⁸O values of quartz and calcite from the andesite and basalt flows (700–932 m) have isotopic values which require that the equilibrated water δ¹⁸O values increase slightly (− 11.3 to −9.2‰) with increasing measured temperatures (150–265°C). The lithic tuffs and brecciated lava flows (300–700 m) contain widespread siderite. Calculated oxygen isotopic compositions of waters in equilibrium with siderite generally increase with increasing temperatures (76–100°C). The δ¹⁸O values of siderite probably result from precipitation in water produced by mixing various amounts of the deep hydrothermal water (− 10.5 ‰) with meteoric water (− 15.5 ‰) recharged within the caldera. The δ¹³C values of calcite and siderite decrease with increasing temperatures and show that these minerals precipitated in isotopic equilibrium with CO₂ of about −8 ‰.The δ¹⁸O values of weakly altered (<5% alteration of plagioclase) whole-rock samples decrease with increasing temperatures above 100°C, indicating that exchange between water and rock is kinetically controlled. The water/rock mass ratios decrease with decreasing temperatures. The δ¹⁸O values of rocks from the bottom of Newberry 2 show about 40% isotopic exchange with the reservoir water.The calculated δ¹⁸O and δD values of bottom hole water determined from the fluid produced during the 20 hour flow test are −10.2 and −109‰, respectively. The δD value of the hydrothermal water indicates recharge from outside the caldera.     

Geochemistry and isotopes of fluids from sulphur springs, Valles Caldera, New Mexico

Detailed geochemistry supported by geologic mapping has been used to investigate Sulphur Springs, an acid-sulfate hot spring system that issues from the western flank of the resurgent dome inside Valles Caldera. The most intense activity occurs at the intersection of faults offsetting caldera-fill deposits and post-caldera rhyolites. Three geothermal wells in the area have encountered pressures <1 MPa and temperatures of 200°C at depths of 600 to 1000 m. Hot spring and fumarole fluids may discharge at boiling temperatures with pH 1.0 and SO₄ 8000 mg/l. These conditions cause argillic alterations throughout a large area.Non-condensible gases consist of roughly 99% CO₂ with minor amounts of H₂S, H₂, and CH₄. Empirical gas geothermometry suggests a deep reservoir temperature of 215 to 280°C. Comparison of ¹³C and ¹⁸O between CaCO₃ from well cuttings and CO₂ from fumarole steam indicates a fractionation temperature between 200 and 300°C by decarbonation of hydrothermally altered Paleozoic limestone and vein calcite in the reservoir rocks. Tritium concentrations obtained from steam condensed in a mudpot and deep reservoir fluids (Baca #13, 278°C) are 2.1 and 1.0 T.U. respectively, suggesting the steam originates from a reservoir whose water is mostly >50 yrs old. Deuterium contents of fumarole steam, deep reservoir fluid, and local meteoric water are practically identical even though ¹⁸O contents range through 4‰, thus, precipitation on the resurgent dome of the caldera could recharge the hydrothermal system by slow percolation. From analysis of D and ¹⁸O values between fumarol steam and deep reservoir fluid, steam reaches the surface either (1) by vaporizing relatively shallow groundwater at 200°C or (2) by means of a two-stage boiling process through an intermediate level reservoir at roughly 200°C.Although many characteristics of known vapor-dominated geothermal systems are found at Sulphur Springs, fundamental differences exist in temperature and pressure of our postulated vapor-zone. We propose that the reservoir beneath Sulphur Springs is too small or too poorly confined to sustain a “true” vapor-dominated system and that the Sulphur Springs system may be a “dying” vapor-dominated system that has practically boiled itself dry.     

Geothermal assessment of the island of ischia (southern Italy) from isotopic and chemical composition of the delivered fluids

The Ischia geothermal system is hosted by silicic rocks of the Quaternary Potassic Roman Province, in southern Italy. Exploration drilling down to 1156 m depth in the mid-1950s provided information on boiling profiles (up to 250°C) and on the depth and permeability of the potential reservoirs. Discharge fluid samples were collected and analyzed to define the inflow of surrounding seawater (C1 ranges from 2.5 to 20 g/kg) into the system.Analyses of samples from surface manifestations and shallow wells collected during 1983 and 1988 point to the existence of three distinct mixing regimes, involving three water components. A dishomogeneous body of diluted water (Cl less than 2.5 g/kg), that occurs at depths > 700 m and reequilibrates at 240°C at least, is overlain by an aquifer of groundwater variably mixed with variably seawater (Cl from 4 to 10 g/kg), which tends to reequilibrate at 160°C. Steam-heated waters locally develop and act as dilutants of the rising geothermal fluids.Dilution, mixing, and evaporation of the ascending chloride fluids are supported by oxygen and hydrogen isotopic data the thermal waters being enriched in ¹⁸O and D with respect to local meteoric water by up to 7 and 30‰, respectively. The relative composition of the major cations in thermal solutions was used to discriminate the two main groups of thermal waters, the reservoir temperatures of which are estimated from the Na/K-gethermometer. K-Mg geothermometer indicates reequilibration in near-surface conditions.The isotopic composition of the fumarolic steam varies from −7 to −12‰ in ∂⁸O and from − 35 to − 70‰ in ∂D, in agreement with a deep mixed fluid that boils adiabatically from 240 to 80°C. The deuterium content of the H₂O-H₂ pair gives enrichment factor of about 830‰, corresponding to equilibrium temperature conditions slightly higher than the surface boiling temperatures. The ∂¹³C of CO₂is almost constant at −4.5‰ (1δ=0.4), suggesting an important magmatic contribution, and the ∂¹⁸O values of CO₂appears to in equilibrium with accompanying steam at the measured temperatures.The CO₂/Ar and H₂/Ar chemical ratios have been used to derive aquifer temperatures, the values obtained being consistent with those of solute geothermometers.     

A reconnaissance geochemical study of La Primavera geothermal area, Jalisco, Mexico

The Sierra La Primavera, a late Pleistocene rhyolitic caldera complex in Jalisco, México, contains fumaroles and large-discharge 65°C hot springs that are associated with faults related to caldera collapse and to later magma insurgence. The nearly-neutral, sodium bicarbonate, hot springs occur at low elevations at the margins of the complex, whereas the water-rich fumaroles are high and central.The Comisión Federal de Electricidad de México (CFE) has recently drilled two deep holes at the center of the Sierra (PR-1 and Pr-2) and one deep hole at the western margin. Temperatures as high as 285°C were encountered at 1160 m in PR-1, which produced fluids with 820 to 865 mg/kg chloride after flashing to one atmosphere. Nearby, PR-2 encountered temperatures to 307°C at 2000 m and yielded fluids with chloride contents fluctuating between 1100 and 1560 mg/kg after flashing. Neither of the high-temperature wells produced steam in commercial quantities. The well at the western margin of the Sierra produced fluids similar to those from the hot springs. The temperature reached a maximum of 100°C near the surface and decreased to 80°C at 2000 m.Various geothermometers (quartz conductive, Na/K, Na-K-Ca, δ¹⁸O(SO₄-H₂O) and D/H (steam-water) all yield temperatures of 170 ± 20°C when applied to the hot spring waters, suggesting that these spring waters flow from a large shallow reservoir at this temperature. Because the hot springs are much less saline than the fluids recovered in PR-1 and PR-2, the mixed fluid in the shallow reservoir can contain no more than 10–20% deep fluid. This requires that most of the heat is transferred by steam. There is probably a thin vapor-dominated zone in the central part of the Sierra, through which steam and gases are transferred to the overlying shallow reservoir. Fluids from this reservoir cool from 170°C to 65°C by conduction during the 5–7 km of lateral flow to the hot springs.     

Hydrothermal vent waters at 13°N on the East Pacific Rise: isotopic composition and gas concentration

Gas concentrations and isotopic compositions of water have been measured in hydrothermal waters from 13°N on the East Pacific Rise. In the most Mg-depleted samples ( 5 × 10⁻³ moles/kg) the gas concentrations are: 3–4.5 × 10⁻⁵ cm³ STP/kg helium, 0.62–1.24 cm³ STP/kg CH₄, 10.80–16.71 × 10⁻³ moles/kg CO₂. The samples contain large quantities (95–126 cm³/kg) of H₂ and some carbon monoxide (0.26–0.36 cm³/kg) which result from reaction with the titanium sampling bottles. δ¹³C in methane and CO₂ (−16.6 to −19.5 and −4.1 to −5.5 respectively) indicate temperatures between 475 and 550°C, whereas δ¹³C_CO is compatible with formation by reduction of CO₂ on Ti at 350°C close to the sampling temperature.³He/⁴He are very homogeneous at (7.5 ± 0.1)R_A(³He/⁴He = 1.0 × 10⁻⁵) and very similar to already published data as well as CH₄/³He ratios between 1.4 and 2.1 × 10⁶.¹⁸O and D in water show enrichments from 0.39 to 0.69‰ and from 0.62 to 1.49‰ respectively. These values correspond to W/R ratios of 0.4–7. The distinct¹⁸O enrichments indicate that the isotopic composition of the oceans is not completely buffered by the hydrothermal circulations. The³He-enthalpy relationship is discussed in terms of both hydrothermal heat flux and³He mantle flux.     

O/O mapping and hydrogeology of a short-lived (≈10 years) fumarolic (>500°C) meteoric–hydrothermal event in the upper part of the 0.76 Ma Bishop Tuff outflow sheet, California

¹⁸O/¹⁶O data from the 200-m-thick, 0.76 Ma Bishop Tuff outflow sheet provide evidence for a vigorous, short-lived (≈10 years), high-temperature, fumarolic meteoric–hydrothermal event. This is proved by: (1) the juxtaposition in the upper, partially welded Bishop Tuff of low-¹⁸O groundmass/glass (δ¹⁸O=−5 to +3) with coexisting quartz and feldspar phenocrysts having magmatic δ¹⁸O values (+8.7±0.3; +7.5±0.3); and (2) the fact that these kinds of ¹⁸O/¹⁶O signatures correlate very well with morphological features and mapped zones of fumarolic activity. Profiles of δ¹⁸O with depth in the Bishop Tuff within the fumarole area define a 40- to 50-m-thick, low-¹⁸O, stratigraphic zone that is sandwiched between the essentially unwelded near-surface portion of the tuff and an underlying, densely welded black tuff that displays magmatic ¹⁸O/¹⁶O values. Shallow-dipping columnar joints and other fumarolic features (i.e., subhorizontal tubular conduits and steep fissures) correlate very well with these pervasively devitrified, low-¹⁸O zones. The base of the low-¹⁸O zone is extremely sharp (3‰ per meter) and is located directly above the transition from partially welded tuff to densely welded black tuff. The observed average whole-rock ¹⁸O-depletions within this low-¹⁸O zone are about 6–7‰, requiring meteoric water/rock ratios in excess of 0.24 in mass units. Rainfall on the surface of the tuff would not have been high enough to supply this much H₂O in the short lifetime of fumarolic activity, suggesting that some recharge must have been from groundwater flow through the upper part of the tuff, above the sloping (1°–5°) top of the impermeable lower zone. This is compatible with the observation that the fumarolic areas roughly correlate with the preeruptive regional drainage pattern. Some of this recharge may in part have been from the lake that filled Long Valley caldera, which was dammed by the Bishop Tuff up to the level of this boundary between the partially and densely welded zones (≈7000 ft, the elevation of the highest Long Valley Lake shorelines). Gazis et al. had previously shown that the 2.8-Ma intracaldera Chegem Tuff from the Caucasus Mountains exhibits exactly the same kind of ¹⁸O-signature that we have correlated with fossil fumaroles in the Bishop Tuff outflow sheet. Although not recognized as such by McConnell et al.; ¹⁸O/¹⁶O data from drill-hole samples from the intracaldera Bishop Tuff in Long Valley also display this characteristic ¹⁸O signature (i.e., analogous δ¹⁸O-depth profiles, as well as low-¹⁸O groundmass coexisting with high-¹⁸O feldspar phenocrysts). This fumarolic ¹⁸O/¹⁶O signature is observed to much greater depths (≈650–750 m) in the intracaldera tuffs (≈1500 m thick) than it is in the ≈200-m-thick Bishop Tuff outflow sheet (≈80 m depth).     

Oxygen diffusion rates in quartz exchanged with CO2

Oxygen self diffusion rates were determined in quartz samples exchanged with¹⁸O-enriched CO₂ between 745 and 900°C and various pressures, and the diffusion profiles were measured using an ion microprobe. The activation energy (Q) and preexponential factor (D₀) at P(CO₂) = P(tot) = 100 bar, for diffusion parallel to the c-axis are 159 ( ± 13) kJ/g atom and 2.10 (+0.75/ −0.55) × 10⁻⁸ cm²/s. This rate is approximately 100 times slower than that obtained from hydrothermal experiments and 100 times faster than a previous 1-bar quartz-O₂ exchange experiment. The oxygen diffusion rate measured at 0.6 bar, 888°C, and at 900°C in vacuum is in agreement with the previous 1-bar exchange experiments with¹⁸O₂. The effect of higher CO₂ pressures is small. At 900°C, the diffusion rate exchanged with CO₂ is = 2.35 × 10⁻¹⁵ cm²/s at 100 bar, 2.24 × 10⁻¹⁵ cm²/s at 3.45 kbar and 8.13 × 10⁻¹⁵ cm²/s at 7.2 kbar.There is probably a diffusing species, other than oxygen, that enhances the oxygen diffusion rate in these quartz-CO₂ systems, relative to that occurring at very low pressures or in a vacuum. The effect of this diffusing species, however, is not as strong as that associated with H₂O. Preserved oxygen isotope fractionations between coexisting minerals in a slowly cooled, high-grade metamorphic terrane will vary depending upon whether a water-rich phase was present or not. Closure temperatures will be approximately 100°C higher in rocks where no water-rich phase was present during cooling. The measured fractionations between coexisting minerals in metamorphic rocks may potentially be used as a sensor of water presence during retrogression.     

Paleomagnetism of the triassic nicola volcanics and geotectonics of the quesnellia subterrane of terrane I, British Columbia

Paleomagnetic data from 46 sites (674 specimens) of the Westcoast Crystalline Gneiss Complex on the west coast of Vancouver Island using AF and thermal demagnetization methods yields a high blocking temperature WCB component (> 560°C) with a pole at 335°W, 66°N (δ_p = 4°, δ_m = 6°) and a lower coercivity WCA component ( 25 mT, < 500°C) with a pole at 52°W, 79°N (δ_p = 7°, δ_m = 8°). Further thermal demagnetization data from 24 sites in the Jurassic Island Intrusions also defines two high blocking temperature components. The IIA component pole is at 59°W, 79°N (δ_p = 7°, δ_m = 8°) and IIB pole at 130°W, 73°N (δ_p = 12°, δ_m = 13°). Combined with previous data from the Karmutsen Basalts and mid-Tertiary units on Vancouver Island and from the adjacent Coast Plutonic Complex, the geotectonic motions are examined for the Vancouver Island segment of the Wrangellian Subterrane of composite Terrane II of the Cordillera. The simplest hypothesis invokes relatively uniform translation for Terrane II from Upper Triassic to Eocene time producing 39° ± 6° of northward motion relative to the North American craton, combined with 40° of clockwise rotation during the Lower Tertiary.     

Evaluation of approximations in modelling the cooling of magmatic bodies

The cooling of a magmatic intrusion is simulated by a simple model of a non-homogeneous earth, with thermal properties depending on temperature, in which heat transfer is assumed to take place by conduction only. The mathematical problem consists in solving a non-linear partial differential equation with continuity conditions on temperature and heat flux imposed at the contacts between different rocks. This has been done numerically by a finite difference method. The model is then adopted as “reality” against which a number of commonly used approximations are tested. It is found that the effect of latent heat liberation can be reasonably taken into account by attributing an effective initial temperature to the magma (errors within 20°C for t > 10⁵ years, when the temperature of the magma is still as high as 600°C); the effective specific heat approximation does not work as well. The dependence of thermal conductivity and specific heat on temperature may be eliminated by maintaining the errors within 30°C for t < 5 × 10⁵ years. The assumption that magma and country rocks have the same thermal properties allows an estimate of the temperature field in the host rocks with errors of 50°C at most. The assumption that all rocks have the same constant conductivity yields results that are far from “reality” (errors of 100–200°C even at shallow depth).     

On the evolution of the geothermal regime of the North China Basin

Recent heat flow and regional geothermal studies indicate that the North China Basin is characterized by relatively high heat flow compared with most stable areas in other parts of the world, but lower heat flow than most active tectonic areas. Measured heat flow values range from 61 to 74 mW m⁻². The temperature at a depth of 2000 m is generally in the range 75 to 85°C, but sometimes is 90°C or higher. The geothermal gradient in Cenozoic sediments is in the range 30 to 40°C/km for most of the area. The calculated temperature at the Moho is 560 and 640°C for surface heat flow values of 63 and 71 mW m⁻², respectively. These thermal data are consistent with other geophysical observations for the North China Basin. Relatively high heat flow in this area is related to Late Cretaceous-Paleogene rifting as described in this paper.     

High-resolution garnet chronometry and the rates of metamorphic processes

Garnets in an amphibolite-facies metasediment from Sulitjelma, North Norway yield precise and concordant SmNd, UPb and RbSr ages that relate directly to the pressure (P) and temperature (T) conditions of mineral growth. Differential mineral reaction between graphitic and non-graphitic layers within this sample preserves a record of theP-T and time (t) history experienced during Barrovian regional metamorphism. Garnets in graphitic layers grew during prograde metamorphism at462 ± 16°C and5.2 ± 0.5 kbar under conditions of lowaH₂O, and yield indistinguishable¹⁴⁷Sm¹⁴³Nd and²³⁸U²⁰⁶Pb ages of434.1 ± 1.2 Ma and433.9 ± 1.0 Ma, respectively. In contrast, garnet growth in adjacent graphite-free layers did not occur untilP-T conditions of540 ± 18°C and8.0 ± 1.0 kbar were attained, with continued growth in response to minor heating and decompression with final matrix equilibration at544 ± 16°C and7.0 ± 1.0 kbar. The inclusion-free garnet rims in this assemblage record indistinguishable¹⁴⁷Sm¹⁴³Nd and²³⁸U²⁰⁶Pb ages of424.6 ± 1.2 Ma and423.4± 1.7 Ma, respectively. These results provide precise estimates for average heating and burial rates during prograde metamorphism of 8.6_−4.4^+7.5°C Ma⁻¹ and 0.8_−0.5^+0.9 km Ma⁻¹, respectively. Rb and Sr exchange between coexisting silicates in the graphite-free assemblage continued for some 37 Ma after the “peak” of metamorphism, and require an average cooling rate of about 4.0°C Ma⁻¹ during uplift. These results illustrate a clear relationship between reaction history and the timing of mineral growth and provide definitive constraints on the rates of thermal and tectonic processes accompanying regional metamorphism.     

Water–Rock interactions in the basement beneath Long Valley Caldera: an oxygen isotope study of the Long Valley Exploratory Well drill cores

The Long Valley Exploratory Well, at the center of the Resurgent Dome of Long Valley caldera, penetrated pre-caldera basement rocks at a depth of 2101.72–2313.0 m, beneath the caldera-forming Bishop Tuff and post-caldera Early Rhyolite. The basement rocks contain prominent quartzites, with ubiquitous milky white quartz veins (with minor calcite and pyrite) and fractures of varied orientation and geometry. The other members of the basement sequence are very fine-grained quartz-rich graphitic pelites with calcite veins, spotted hornfels, and shallow intrusive rocks. Previous studies established the presence of a post-caldera, paleohydrothermal system (500–100 ka) to a depth of 2000 m that affected the Bishop Tuff and a recent (40 ka to present) hydrothermal system at shallow depth (<1 km). The deeper extent of these hydrothermal activities is established in this paper by a detailed oxygen isotope analysis of the drill core samples. 238 analyses of δ¹⁸O in 50 quartz veins within the 163.57 m depth interval of basement rocks reveal extreme heterogeneity in δ¹⁸O values (8–19.5‰). Majorities of the 84 bulk analyses of quartzites show variation of δ¹⁸O within a narrow range of 14–16‰. However, certain samples of these quartzites near the contacts with veins and fractures exhibit sharp drops in δ¹⁸O. The interbedded pelitic rocks and spotted hornfels have whole-rock δ¹⁸O ranging from 2.2 to 11.8‰. Clear, euhedral vuggy quartz that partially fills earlier open fractures in both the quartzites and quartz veins, has distinctive δ¹⁸O, ranging between −3.2 and +8.4‰. Low values of δ¹⁸O are also found in the hydrothermal minerals and whole rocks adjacent to the thin veins, clearly indicating infiltration of meteoric water. Three distinct observed patterns of fractionation in δ¹⁸O between veins and host quartzites are analyzed with the principles of mass balance, equilibrium oxygen isotope fractionation in closed system, and kinetically controlled oxygen isotope exchange in an open system. This analysis suggests that the early quartz veins formed due to a magmatic-hydrothermal activity with no influx of external water once the system comprising the sedimentary envelope and a magmatic-hydrothermal fluid phase became closed. Two-stage isotopic exchange processes caused fractionation in the δ values that originally formed arrays with slope 1 in a δ_{vein quartz}–δ_{host quartzite} space. Another array in the same space, with near zero slope was also formed due to variation in temperature, initial isotopic compositions of the quartzite sequence and the fluid phase. Variation in temperature was mostly in the range of 300–400°C giving Δ (=δ_{vein quartz}–δ_{host quartzite})≈−2.8 to +2.8. The δ¹⁸O of the fluid could range from −5 to +10; however a narrower range of +5 to +10 can explain the data. This episode of hydrothermal activity could take place either as a single pulse or in multiple pulses but each as a closed system. A later, fracture-controlled, meteoric water (δ¹⁸O−0.46 to −12.13) flow and interaction (at 250°C) is interpreted from the analysis of δ¹⁸O values of the coexisting quartz and calcite pairs and existence of markedly ¹⁸O-depleted pelitic horizons interbedded with ¹⁸O-enriched quartzite layers. Thus, the interpreted earlier magmatic-hydrothermal activity was overprinted by a later meteoric-hydrothermal activity that resulted in steep arrays of δ¹⁸O values in the δ_{vein quartz}–δ_{host quartzite} space. Calculations show that the likely life span of the post-caldera, hydrothermal activity in the depth range of 2.1–2.3 km beneath Long Valley was 0.08–0.12 Ma. Diffusive ±advective transport of oxygen isotopes from fracture-channelized meteoric water to nearly impermeable wall rocks caused a lowering of δ¹⁸O values in the quartz over short distances and in calcites over greater distances. Thus, the hydrothermal activity appears pervasive even though the meteoric water flow was primarily controlled by fractures.     

Chemical and stable isotopic models for boundary creek warm springs, southwestern Yellowstone National Park, Wyoming

Thermal springs of the Boundary Creek hydrothermal system in the southwestern part of Yellowstone Park outside the caldera boundary vary in chemical and isotopic composition, and temperature. The diversity may be accounted for by a combination of processes including boiling of a deep thermal water, mixing of the deep thermal water with cool meteoric water and/or with condensed steam or steam-heated meteoric water, and chemical reactions with surrounding rocks. Dissolved-silica, Na⁺, K⁺ and Ca²⁺ contents of the thermal springs could result from a thermal fluid with a temperature of 200 ± 20°C. Chloride-enthalpy and silica-enthalpy mixing models suggest mixing of 230°C, 220 mg/l Cl⁻ thermal water with cool, low-Cl⁻ components. A 350 to 390°C component with Cl⁻ ≥ 300 mg/l is possibly present in thermal springs inside the caldera but is not required to fit observed spring chemical and isotopic compositions. Irreversible mass transfer models in which a low-temperature water reacts with volcanic glass as it percolates downward and warms, can account for observed pH and dissolved-silica, K⁺, Na⁺, Ca²⁺ and Mg²⁺ concentrations, but produces insufficient Cl⁻ or F⁻ for measured concentrations in the warm springs. The ratio of a_Na/a_H, and Cl⁻ are best accounted for in mixing models. The water-rock interaction model fits compositions of acid-sulfate waters observed at Summit Lake and of low-Cl⁻ waters involved in mixing.The cold waters collected from southwestern Yellowstone Park have δD values ranging from −118 to −145 per mil and δ¹⁸O values of −15.9 to −19.4 per mil. Two samples from nearby Island Park have δD values of −112 and −114 per mil and δ¹⁸O values of −15.1 and −15.3 per mil. All samples of thermal water plot significantly to the right of the meteoric water line. The low Cl⁻ and variable δD values of the thermal waters indicate isotopic compositions are derived by extensive dilution with cold meteoric water and by steam separation on ascent to the surface. Many of the hot springs with higher δD values may contain in addition a significant amount of high-D, low-Cl⁻, acid-sulfate or steam-heated meteoric water. Mixing models, Cl⁻ content and isotopic compositions of thermal springs suggest that 30% or less of a deep thermal component is present. For example, the highest-temperature springs from Three Rivers, Silver Scarf and Upper Boundary Creek thermal areas contain up to 70% cool meteoric water and 30% hot water components, springs at Summit Lake and Middle Boundary Creek spring 57 are acid-sulfate or steam-heated meteoric water; springs 27 and 48 from Middle Boundary Creek and 49 from Mountain Ash contain in excess of 50% acid-sulfate water; and Three Rivers spring 46 and Phillips could result from mixing hot water with 55% cool meteoric water followed by mixing of acid-sulfate water. Extensive dilution by cool meteoric water increases the uncertainties in quantity and nature of the deep meteoric, thermal component.     

Palaeomagnetic timing of the rotation and translation of Cyprus

The geological evolution of the Mesozoic Troodos Ophiolite Complex in Cyprus, and the tectonic nature and timing of the palaeomagnetically indicated anticlockwise rotation of Cyprus of some 80° and ca. 15° northward translation, have been open for debate for some time. New palaeomagnetic data from 18 sites ( 180samples) in the post-ophiolite sediments, ranging in age from Upper Cretaceous to Upper Miocene, are presented. Most of the sites are of normal geomagnetic polarity, but indications of reversed polarity have been found in an older group of sediments (the Lefkara Formation of Upper Palaeocene age).Six sites from the older group of sediments (Upper Cretaceous to Eocene in age) give a site mean direction of the AF cleaned sediments of (D, I) = (323°, 29°) with α₉₅ = 18°, while 5 sites from a younger group of sediments (Oligocene to Miocene in age) give a cleaned site mean direction of (D, I) = (334°, 58°) with α₉₅ = 9°. These and published data suggest that an anticlockwise rotation of Cyprus of 60 ± 10° occurred early during the post-igneous evolution of the Cyprus oceanic crust between 90 and 50Ma, leaving only a minor anticlockwise rotation of 20 ± 10° to occur during the last 50 Ma. It is furthermore concluded that the northward translation of Cyprus of 15° mostly took place during the last 30Ma.It thus appears that a fairly rapid rotation of the Cyprus microplate first took place in the Late Cretaceous and Early Tertiary time with an average angular velocity of 1–2°/Ma, during which the northward translation was minor or negligible. In the latter half of the Tertiary, the sense of movement appears to have radically changed, the northward translation now being dominant with an average velocity of 5–6cm/yr. This temporal evolution is found to be in good agreement with the Mesozoic and Tertiary movements of the African lithospheric plate relative to Europe, as evidenced from the Atlantic sea-floor magnetic anomaly spreading history.     

O, S, Sr, and Pb isotope variations in volcanic rocks from the Northern Mariana Islands: implications for crustal recycling in intra-oceanic arcs

Quaternary lavas from the Northern Mariana Islands have respective O- and S-isotope ranges ofδ¹⁸O = +5.7 to +6.6 (‰ SMOW) andδ³⁴S = +2.0to+20.7 (‰ CDT). Chemically evolved andesites and dacites with meanδ¹⁸O = +6.3 ± 0.2 are slightly enriched in¹⁸O with respect to unfractionated basalts of<53%SiO₂ with meanδ¹⁸O = +6.0 ± 0.1. This¹⁸O enrichment can be explained in terms of differentiation of parental mafic magmas havingδ¹⁸O values between +5.7 to +6.2‰ through closed-system crystal fractionation because the lavas from all nine islands of the arc define a coherentδ¹⁸OSiO₂ trend. The S-isotope composition of oxidized magmas is not modified extensively through the degassing of SO₂; therefore, the meanδ³⁴S value of ca. +11‰ for the Mariana lavas is considered to be representative of their source region.The enrichment of¹⁸O and³⁴S in Mariana Arc parental magmas relative to ocean floor basalts withδ¹⁸O ca. + 5.7‰ andδ³⁴S = ca.0.3‰ is attributed to the recycling of¹⁸O- and³⁴S-rich crustal components (sediment withδ¹⁸O = ca. + 25‰ and seawater sulfate withδ³⁴S = ca. +20‰ into the upper mantle source region for these arc magmas. This interpretation is consistent with enrichments of radiogenic Sr and Pb in the same lavas relative to ocean-floor basalts erupted either side of the arc, which are presumed to share a common upper mantle source. This enrichment is considered to reflect the mixing of two components, one having a typical upper mantle composition and the other having a more radiogenic character similar to that of western Pacific pelagic sediments.