Assessment of natural attenuation of aromatic hydrocarbons in groundwater near a former manufactured-gas plant, South Carolina, USA

 Shallow, anaerobic groundwater near a former manufactured-gas plant (MGP) in Charleston, South Carolina, USA, contains mono- and polycyclic aromatic hydrocarbons (MAHs and PAHs, respectively). Between 1994 and 1997, a combination of field, laboratory, and numerical-flow and transport-model investigations were made to assess natural attenuation processes affecting MAH and PAH distributions. This assessment included determination of adsorption coefficients (K _ad ) and first-order biodegradation rate constants (K _bio ) using aquifer material from the MGP site and adjacent properties. Naphthalene adsorption (K _ad =1.35×10^–7 m³/mg) to aquifer sediments was higher than toluene adsorption (K _ad =9.34×10^–10 m³/mg), suggesting preferential toluene transport relative to naphthalene. However, toluene and benzene distributions measured in January 1994 were smaller than the naphthalene distribution. This scenario can be explained, in part, by the differences between biodegradation rates of the compounds. Aerobic first-order rate constants of ¹⁴C-toluene, ¹⁴C-benzene, and ¹⁴C-naphthalene degradation were similar (–0.84, –0.03, and 0.88 day^–1, respectively), but anaerobic rate constants were higher for toluene and benzene (–0.002 and –0.00014 day^–1, respectively) than for naphthalene (–0.000046 day^–1). Both areal and cross-sectional numerical simulations were used to test the hypothesis suggested by these rate differences that MAH compounds will be contained relative to PAHs. Predictive simulations indicated that the distributions of toluene and benzene reach steady-state conditions before groundwater flow lines discharge to an adjacent surface-water body, but do discharge low concentrations of naphthalene. Numerical predictions were "audited" by measuring concentrations of naphthalene, toluene, and benzene at the site in early 1997. Measured naphthalene and toluene concentrations were substantially reduced and the areal extent of contamination smaller than was both observed in January 1994 and predicted for 1997. Measured 1997 benzene concentrations and distribution were shown to be relatively unchanged from those measured in 1994, and similar to predictions for 1997. Received: 26 June 1997 · Accepted: 25 August 1997     

Structural thermal behavior of paragonite and its dehydroxylate: a high-temperature single-crystal study

 The lattice constants of paragonite-2M₁, NaAl₂(AlSi₃)O₁₀(OH)₂, were determined to 800 °C by the single-crystal diffraction method. Mean thermal expansion coefficients, in the range 25–600 °C, were: α_a = 1.51(8) × 10⁻⁵, α_b = 1.94(6) × 10⁻⁵, α_c = 2.15(7) ×  10⁻⁵ °C⁻¹, and α_V = 5.9(2) × 10⁻⁵ °C⁻¹. At T higher than 600 °C, cell parameters showed a change in expansion rate due to a dehydroxylation process. The structural refinements of natural paragonite, carried out at 25, 210, 450 and 600 °C, before dehydroxylation, showed that the larger thermal expansion along the c parameter was mainly due to interlayer thickness dilatation. In the 25–600 °C range, Si,Al tetrahedra remained quite unchanged, whereas the other polyhedra expanded linearly with expansion rate proportional to their volume. The polyhedron around the interlayer cation Na became more regular with temperature. Tetrahedral rotation angle α changed from 16.2 to 12.9°. The structure of the new phase, nominally NaAl₂ (AlSi₃)O₁₁, obtained as a consequence of dehydroxylation, had a cell volume 4.2% larger than that of paragonite. It was refined at room temperature and its expansion coefficients determined in the range 25–800 °C. The most significant structural difference from paragonite was the presence of Al in fivefold coordination, according to a distorted trigonal bipyramid. Results confirm the structural effects of the dehydration mechanism of micas and dioctahedral 2:1 layer silicates. By combining thermal expansion and compressibility data, the following approximate equation of state in the PTV space was obtained for paragonite: V/V ₀ = 1 + 5.9(2) × 10⁻⁵ T(°C) − 0.00153(4) P(kbar). Received: 12 July 1999 / Revised, accepted: 7 December 1999     

Enrichment of PGE through interaction of evolved boninitic magmas with early formed cumulates in a gabbro–breccia zone of the Mesoarchean Nuasahi massif (eastern India)

The Mesoarchean Nuasahi chromite deposits of the Singhbhum Craton in eastern India consist of a lower chromite-bearing ultramafic unit and an upper magnetite-bearing gabbroic unit. The ultramafic unit is a ∼5 km long and ∼400 m wide linear belt trending NNW-SSE with a general north-easterly dip. The chromitite ore bodies are hosted in the dunite that is flanked by the orthopyroxenite. The rocks of the ultramafic unit including the chromitite crystallized from a primitive boninitic magma, whereas the gabbro unit formed from an evolved boninitic magma. A shear zone (10–75 m wide) is present at the upper contact of the ultramafic unit. This shear zone consists of a breccia comprising millimeter- to meter-sized fragments of chromitite and serpentinized rocks of the ultramafic unit enclosed in a pegmatitic and hybridized gabbroic matrix. The shear zone was formed late synkinematically with respect to the main gabbroic intrusion and intruded by a hydrous mafic magma comagmatic with the evolved boninitic magma that formed the gabbro unit. Both sulfide-free and sulfide-bearing zones with platinum group element (PGE) enrichment are present in the breccia zone. The PGE mineralogy in sulfide-rich assemblages is dominated by minerals containing Pd, Pt, Sb, Bi, Te, S, and/or As. Samples from the gabbro unit and the breccia zone have total PGE concentrations ranging from 3 to 116 ppb and 258 to 24,100 ppb, respectively. The sulfide-rich assemblages of the breccia zone are Pd-rich and have Pd/Ir ratios of 13–1,750 and Pd/Pt ratios of 1–73. The PGE-enriched sulfide-bearing assemblages of the breccia zone are characterized by (1) extensive development of secondary hydrous minerals in the altered parts of fragments and in the matrix of the breccia, (2) coarsening of grain size in the altered parts of the chromitite fragments, and (3) extensive alteration of primary chromite to more Fe-rich chromite with inclusions of chlorite, rutile, ilmenite, magnetite, chalcopyrite, and PGE-bearing chalcogenides. Unaltered parts of the massive chromitite fragments from the breccia zone show PGE ratios (Pd/Ir = 2.5) similar to massive chromitite (Pd/Ir = 0.4–6.6) of the ultramafic unit. The Ir-group PGE (IPGE: Ir, Os, Ru) of the sulfide-rich breccia assemblages were contributed from the ultramafic–chromitite breccia. Samples of the gabbro unit have fractionated primitive mantle-normalized patterns, IPGE depletion (Pd/Ir = 24–1,227) and Ni-depletion due to early removal of olivine and chromite from the primitive boninitic magma that formed the ultramafic unit. Samples of the gabbro and the breccia zone have negative Nb, Th, Zr, and Hf anomalies, indicating derivation from a depleted mantle source. The Cu/Pd ratios of the PGE-mineralized samples of the breccia zone (2.0 × 10³–3.2 × 10³) are lower than mantle (6.2 × 10³) suggesting that the parental boninitic magma (Archean high-Mg lava: Cu/Pd ratio ∼1.3 × 10³; komatiite: Cu/Pd ratio ∼8 × 10³) was sulfur-undersaturated. Samples of the ultramafic unit, gabbro and the mineralized breccia zone, have a narrow range of incompatible trace element ratios indicating a cogenetic relationship. The ultramafic rocks and the gabbros have relatively constant subchondritic Nb/Ta ratios (ultramafic rocks: Nb/Ta = 4.1–8.8; gabbro unit: Nb/Ta = 11.5–13.2), whereas samples of the breccia zone are characterized by highly variable Nb/Ta ratios (Nb/Ta = 2.5–16.6) and show evidence of metasomatism. The enrichment of light rare earth element and mobile incompatible elements in the mineralized samples provides supporting evidence for metasomatism. The interaction of the ultramafic fragments with the evolved fluid-rich mafic magma was key to the formation of the PGE mineralization in the Nuasahi massif.     

Crystallographic preferred orientation and microstructure of experimentally deformed Braubach galena ore with emphasis on the relation to diffusional processes

Sample cylinders of two galena ore hand specimens from Braubach, Germany were axially shortened in the strain rate range 5 × 10⁻⁵ s⁻¹–5 × 10⁻⁷ s⁻¹ at a confining pressure of 200 (300) MPa and at temperatures of 20 °C–600 °C. Neutron diffraction analyses of the crystallographic preferred orientation (texture) were carried out before and after experimental deformation on the same sample cylinder. Up to a deformation temperature of 300 °C and a strain rate of 5 × 10⁻⁶ s⁻¹ a more or less complete <110> fiber texture develops, the strength of the fiber texture only depending on strain and the strength of the original preferred orientation. At slower strain rate and higher temperature, there is a distinct decrease of the fiber texture development. Diffusional mass transfer starts to become a significant deformation mechanism. Deformation at 500 °C changes the original texture only slightly, which indicates a rapid increase of importance of diffusional flow processes. The alteration of the accompanying sulfosalts indicates that the temperature is high enough for the movement of atoms. The microstructure only reveals remarkable deformation structures at higher strains and in areas of locally higher stresses. Received: 10 June 1997 / Accepted: 14 May 1998     

Metamorphism, isotopic ages and composition of lower crustal granulite xenoliths from the Cretaceous Salta Rift, Argentina

Crustal xenoliths from basanitic dikes and necks that intruded into continental sediments of the Cretaceous Salta Rift at Quebrada de Las Conchas, Provincia Salta, Argentina were investigated to get information about the age and the chemical composition of the lower crust. Most of the xenoliths have a granitoid composition with quartz-plagioclase-garnet-rutile ± K-feldspar as major minerals. The exceedingly rare mafic xenoliths consist of plagioclase-clinopyroxene-garnet ± hornblende. All xenoliths show a well equilibrated granoblastic fabric and the minerals are compositionally unzoned. Thermobarometric calculations indicate equilibration of the mafic xenoliths in the granulite facies at temperatures of ca. 900 °C and pressures of ca. 10 kbar. The Sm-Nd mineral isochron ages are 95.1 ± 10.4 Ma, 91.5 ± 13.0 Ma, 89.0 ± 4.2 Ma (granitoid xenoliths), and 110.7 ± 23.6 Ma (mafic xenolith). These ages are in agreement with the age of basanitic volcanism (ca. 130–100 and 80–75 Ma) and are interpreted as minimum ages of metamorphism. Lower crustal temperature at the time given by the isochrons was above the closure temperature of the Sm-Nd system (>600–700 °C). The Sm-Nd and Rb-Sr isotopic signatures (¹⁴⁷Sm/¹⁴⁴Nd = 0.1225–0.1608; ¹⁴³Nd/¹⁴⁴Nd_{t 0} = 0.512000–0.512324; ⁸⁷Rb/⁸⁶Sr = 0.099–0.172; ⁸⁷Sr/⁸⁶Sr_{t 0} = 0.708188–0.7143161) and common lead isotopic signatures (²⁰⁶Pb/²⁰⁴Pb = 18.43–18.48; ²⁰⁷Pb/²⁰⁴Pb = 15.62–15.70; ²⁰⁸Pb/²⁰⁴Pb = 38.22 –38.97) of the granitoid xenoliths are indistinguishable from the isotopic composition of the Early Paleozoic metamorphic basement from NW Argentina, apart from the lower ²⁰⁸Pb/²⁰⁴Pb ratio of the basement. The Sm-Nd depleted mantle model ages of ca. 1.8 Ga from granitoid xenoliths and Early Paleozoic basement point to a similar Proterozoic protolith. Time constraints, the well equilibrated granulite fabric, P-T conditions and lack of chemical zoning of minerals point to a high temperature in a crust of nearly normal thickness at ca. 90 Ma and to a prominent thermal anomaly in the lithosphere. The composition of the xenoliths is similar to the composition of the Early Paleozoic basement in the Andes of NW Argentina and northern Chile. A thick mafic lower crust seems unlikely considering low abundance of mafic xenoliths and the predominance of granitoid xenoliths. Received: 21 July 1998 / Accepted: 27 October 1998     

Temperature evolution between 50 K and 320 K of the thermal expansion tensor of gypsum derived from neutron powder diffraction data

The unit cell parameters, extracted from Rietveld analysis of neutron powder diffraction data collected between 4.2 K and 320 K, have been used to calculate the temperature evolution of the thermal expansion tensor for gypsum for 50 ≤ T ≤ 320 K. At 300 K the magnitudes of the principal axes are α ₁₁  = 1.2(6) × 10⁻⁶ K⁻¹, α ₂₂  = 36.82(1) × 10⁻⁶ K⁻¹ and α ₃₃  = 25.1(5) × 10⁻⁶ K⁻¹. The maximum axis, α ₂₂ , is parallel to b, and using Institution of Radio Engineers (IRE) convention for the tensor orthonormal basis, the axes α ₁₁ and α ₃₃ have directions equal to (−0.979, 0, 0.201) and (0.201, 0, 0.979) respectively. The orientation and temperature dependent behaviour of the thermal expansion tensor is related to the crystal structure in the I2/a setting. Received 12 February 1998 / Revised, accepted 19 October 1998     

Thermoelastic properties of the high-pressure phase of SnO2 determined by in situ X-ray observations up to 30 GPa and 1400 K

 In situ synchrotron X-ray experiments in the system SnO₂ were made at pressures of 4–29 GPa and temperatures of 300–1400 K using sintered diamond anvils in a 6–8 type high-pressure apparatus. Orthorhombic phase (α-PbO₂ structure) underwent a transition to a cubic phase (Pa3ˉ structure) at 18 GPa. This transition was observed at significantly lower pressures in DAC experiments. We obtained the isothermal bulk modulus of cubic phase K ₀= 252(28) GPa and its pressure derivative K ^′=3.5(2.2). The thermal expansion coefficient of cubic phase at 25 GPa up to 1300 K was determined from interpolation of the P-V-T data obtained, and is 1.7(±0.7) × 10⁻⁵ K⁻¹ at 25 GPa. Received: 7 December 1999 / Accepted: 27 April 2000     

Isotope geochemistry and fluid inclusion study of skarns from Vesuvius

Summary We present new mineral chemistry, fluid inclusion, stable carbon and oxygen, as well as Pb, Sr, and Nd isotope data of Ca-Mg-silicate-rich ejecta (skarns) and associated cognate and xenolithic nodules from the Mt. Somma-Vesuvius volcanic complex, Italy. The typically zoned skarn ejecta consist mainly of diopsidic and hedenbergitic, sometimes “fassaitic” clinopyroxene, Mg-rich and Ti-poor phlogopite, F-bearing vesuvianite, wollastonite, gehlenite, meionite, forsterite, clinohumite, anorthite and Mg-poor calcite with accessory apatite, spinell, magnetite, perovskite, baddeleyite, and various REE-, U-, Th-, Zr- and Ti-rich minerals. Four major types of fluid inclusions were observed in wollastonite, vesuvianite, gehlenite, clinopyroxene and calcite: a) primary silicate melt inclusions (T_HOM = 1000–1050 °C), b) CO₂ ± H₂S-rich fluid inclusions (T_HOM = 20–31.3 °C into the vapor phase), c) multiphase aqueous brine inclusions (T_HOM = 720–820 °C) with mainly sylvite and halite daughter minerals, and d) complex chloride-carbonate-sulfate-fluoride-silicate-bearing saline-melt inclusions (T_HOM = 870–890 °C). The last inclusion type shows evidence for immiscibility between several fluids (silicate melt – aqueous chloride-rich liquid – carbonate/sulfate melt?) during heating and cooling below 870 °C. There is no evidence for fluid circulation below 700 °C and participation of externally derived meteoric fluids in skarn formation. Skarns have considerably variable ²⁰⁶Pb/²⁰⁴Pb (19.047–19.202), ²⁰⁷Pb/²⁰⁴Pb (15.655–15.670), and ²⁰⁸Pb/²⁰⁴Pb (38.915–39.069) and relatively low ¹⁴³Nd/¹⁴⁴Nd (0.51211–0.51244) ratios. The carbon and oxygen isotope compositions of skarn calcites (δ¹³C_V-PDB = −5.4 to −1.1‰; δ18O_V-SMOW = 11.7 to 16.4‰) indicate formation from a ¹⁸O- and ¹³C-enriched fluid. The isotope composition of skarns and the presence of silicate melt inclusion-bearing wollastonite nodules suggests assimilation of carbonate wall rocks by the alkaline magma at moderate depths (< 5 km) and consequent exsolution of CO₂-rich vapor and complex saline melts from the contaminated magma that reacted with the carbonate rocks to form skarns. Received March 1, 2000; revised version accepted November 2, 2000     

Magnesium self-diffusion in orthoenstatite

Magnesium self-diffusion coefficients were determined experimentally for diffusion parallel to each of the three crystallographic directions in natural orthoenstatite (En₈₈Fs₁₂). Experiments were conducted at 1 atm in CO-CO₂ gas mixing furnaces, which provided oxygen fugacities equivalent to the iron-wüstite buffer. Diffusion of ²⁵Mg was induced in polished samples of oriented orthoenstatite using a film of isotopically enriched ²⁵MgO as the source material. Very short (<0.15 μm) diffusional penetration profiles were measured by ion microprobe depth profiling. The diffusion coefficients determined for four temperatures (900, 850, 800, 750 °C) provide the activation energies, E _a , and frequency factors, D _o, where D = D _o exp (−E _a /RT) for Mg self-diffusion parallel to each crystallographic direction: a-axis, E _a  = 360 ± 52 kJ/mole and D _o = 1.10 × 10⁻⁴ m²/s; b-axis, E _a  = 339 ± 77 kJ/mole and D _o = 6.93 × 10⁻⁶ m²/s and c-axis, E _a  = 265 ± 66 kJ/mole and D _o = 4.34 × 10⁻⁹ m²/s. In this temperature range, any possible anisotropy of cation diffusion is very small, however the activation energy for diffusion parallel to the c-axis (001) is the lowest and the activation energies for diffusion parallel to the a-axis (100) and b-axis (010) are higher. Application of these diffusion results to the silicate phases of the Lowicz mesosiderite meteorite provides cooling rates for the silicate portion of the meteorite (4–11 °C/100 years) that are similar, although slower, to previous estimates. These silicate cooling rates are still several orders of magnitude faster than the cooling rates (0.1 °C/10⁶ years) for the metal portions. Received: 22 January 1997 / Accepted: 2 October 1997     

Thermodynamic data of the high-pressure phase Mg5Al5Si6O21(OH)7 (Mg-sursassite)

 Calorimetric and P–V–T data for the high-pressure phase Mg₅Al₅Si₆O₂₁(OH)₇ (Mg-sursassite) have been obtained. The enthalpy of drop solution of three different samples was measured by high-temperature oxide melt calorimetry in two laboratories (UC Davis, California, and Ruhr University Bochum, Germany) using lead borate (2PbO·B₂O₃) at T=700 ^∘C as solvent. The resulting values were used to calculate the enthalpy of formation from different thermodynamic datasets; they range from −221.1 to −259.4 kJ mol⁻¹ (formation from the oxides) respectively −13892.2 to −13927.9 kJ mol⁻¹ (formation from the elements). The heat capacity of Mg₅Al₅Si₆O₂₁(OH)₇ has been measured from T=50 ^∘C to T=500 ^∘C by differential scanning calorimetry in step-scanning mode. A Berman and Brown (1985)-type four-term equation represents the heat capacity over the entire temperature range to within the experimental uncertainty: C _P (Mg-sursassite) =(1571.104 −10560.89×T ^−0.5−26217890.0 ×T ⁻²+1798861000.0×T ⁻³) J K⁻¹ mol⁻¹ (T in K). The P V T behaviour of Mg-sursassite has been determined under high pressures and high temperatures up to 8 GPa and 800 ^∘C using a MAX 80 cubic anvil high-pressure apparatus. The samples were mixed with Vaseline to ensure hydrostatic pressure-transmitting conditions, NaCl served as an internal standard for pressure calibration. By fitting a Birch-Murnaghan EOS to the data, the bulk modulus was determined as 116.0±1.3 GPa, (K ^′=4), V _T,0 =446.49 ³ exp[∫(0.33±0.05) × 10⁻⁴ + (0.65±0.85)×10⁻⁸ T dT], (K _T/T) _P  = −0.011± 0.004 GPa K⁻¹. The thermodynamic data obtained for Mg-sursassite are consistent with phase equilibrium data reported recently (Fockenberg 1998); the best agreement was obtained with Δ_f H ⁰ ₂₉₈ (Mg-sursassite) = −13901.33 kJ mol⁻¹, and S ⁰ ₂₉₈ (Mg-sursassite) = 614.61 J K⁻¹ mol⁻¹. Received: 21 September 2000 / Accepted: 26 February 2001     

Fluid flow in marbles at Jervois, central Australia: oxygen isotope disequilibrium and zoning produced by decoupling of mineralogical and isotopic resetting

The Jervois region of the Arunta Inlier, central Australia, contains para- and orthogneisses that underwent low-pressure amphibolite facies metamorphism (P = 200–300 MPa, T = 520–600 °C). Marble layers cut by metre-wide quartz + garnet ± epidote veins comprise calcite, quartz, epidote, clinopyroxene, grandite garnet, and locally wollastonite. The marbles also contain locally discordant decimetre-thick garnet and epidote skarn layers. The mineral assemblages imply that the rocks were infiltrated by water-rich fluids (XCO₂ = 0.1–0.3) at ∼600 °C. The fluids were probably derived from the quartz-garnet vein systems that represent conduits for fluids exsolved from crystallizing pegmatites emplaced close to the metamorphic peak. At one locality, the marble has calcite (Cc) δ¹⁸O values of 9–18‰ and garnet (Gnt) δ¹⁸O values of 10–14‰. The δ¹⁸O(Gnt) values are only poorly correlated with δ¹⁸O(Cc), and the δ¹⁸O values of some garnet cores are higher than the rims. The isotopic disequilibrium indicates that garnet grew before the δ¹⁸O values of the rock were reset. The marbles contain  ≤15% garnet and, for water-rich fluids, garnet-forming reactions are predicted to propagate faster than O-isotopes are reset. The Sm-Nd and Pb-Pb ages of garnets imply that fluid flow occurred at 1750–1720 Ma. There are no significant age differences between garnet cores and rims, suggesting that fluid flow was relatively rapid. Texturally late epidote has δ¹⁸O values of 1.5–6.2‰ implying δ¹⁸O(H₂O) values of 2–7‰. Waters with such low-δ¹⁸O values are probably at least partly meteoric in origin, and the epidote may be recording the late influx of meteoric water into a cooling hydrothermal system. Received: 29 April 1996 / Accepted: 12 March 1997     

Diffusion of dissolved SiO2 in H2O at 1 GPa, with implications for mass transport in the crust and upper mantle

Cation diffusion rates at 690 ± 30 °C have been calculated by inverse modelling of observed manganese (Mn) zonation profiles in 40 garnets from two kyanite-bearing metapelite samples from the High Himalayan Crystalline Series, Zanskar, northwest India. Knowledge of the initial growth profile of Mn in garnet is a pre-requisite for this technique. Following previous workers we model Mn partitioning into growing garnet in terms of a Rayleigh fractionation process, and demonstrate that the Mn_{garnet:whole rock} partition coefficient is 60–100. Three-dimensional zonation profiles were obtained by successively grinding and polishing ∼1 cm slabs of each sample at 0.1–0.2 mm intervals and analysing the garnets at each stage, thus ensuring that core sections were measured. The diffusion model assumes that garnet has a spherical geometry and behaves as a closed system, and simulates diffusive modification of the hypothetical Mn Rayleigh growth profile for each garnet. The derived measure of the time-integrated diffusion history for each garnet is then combined with radiometric and field-relation constraints for the duration of the Himalayan metamorphic event to calculate cation diffusion rates. The average cation interdiffusion rate calculated for garnets in the two samples examined is (6 ± 3.2) × 10⁻²³ m²s⁻¹. This interdiffusion rate pertains to a temperature of 690 ± 30 °C, which is 0.97 × T _PEAK, the peak temperature conditions experienced by the samples estimated using standard thermobarometric techniques. Garnet compositions are P_y2–17Alm_65–77Gro_6–16Sp_1–17. These new diffusion data are consistent with, and more precise than, existing high-temperature (>1000 °C) experimentally determined diffusion data, although some uncertainties remain difficult to constrain. Qualitative comparison between diffusively modified Mn growth profiles in garnets from the Scottish Dalradian and the Himalayan garnets suggests that the duration of metamorphism affecting the Dalradian garnets was 10–20 times longer than that endured by the Himalayan garnets. Received: 5 June 1996 / Accepted: 29 January 1997     

Iron Diffusion in Single-Crystal Diopside

 Iron tracer diffusion experiments in diopside have been performed using natural and synthetic single crystals of diopside, and stable iron tracers enriched in ⁵⁴Fe, at temperatures in the range 950–1100 °C, total pressure 1 atm, for times up to 29 days. Iron isotope diffusion profiles were determined with an ion microprobe. For experiments performed at log pO₂ = −13, in directions parallel to the c axis and the b axis of two natural, low iron (Fe ∼ 1.8 at %) diopsides, the data obey a single Arrhenius relationship of the form D = 6.22_−5.9 ^+49.6×10⁻¹⁵ exp(−161.5 ± 35.0 kJ mol⁻¹/RT) m² s⁻¹. A single datum for iron diffusion in iron-free, single-crystal diopside at 1050 °C, is approximately 1 order of magnitude slower than in the natural crystals. The pO₂ dependence of iron diffusion in natural crystals at 1050 °C (power exponent = 0.229 ± 0.036) indicates a vacancy mechanism; this is consistent with the results of unpublished atomistic simulation studies. There is no evidence of anisotropy for iron diffusion in diopside. Received: 16 March 1999 / Accepted: 10 April 2000     

Rare-earth elements as source indicators of Pan-African granites from Obudu Plateau,Southeastern Nigeria

The rare-earth element （REE） concentrations of representative granite samples from the southeast of the Obudu Plateau, Nigeria, were analyzed with an attempt to determine the signatures of their source, evolutionary history and tectonic setting. Results indicated that the granites have high absolute REE concentrations （190×10^-6-1191×10^-6; av.=549×10^-6） with the chondrite-normalized REE patterns characterized by steep negative slopes and prominent to slight or no negative Eu anomalies. All the samples are also characterized by high and variable concentrations of the LREE （151×10^-6-1169×10^-6; av.= 466×10^-6）, while the HREE show low abundance （4×10^-6-107×10^-6; av.=28×10^-6）. These are consistent with the variable levels of REE fractionation, and differentiation of the granites. This is further supported by the range of REE contents, the chondrite-normalized patterns and the ratios of LaN/YbN （2.30-343.37）, CeN/YbN （5.94-716.87）, LaN/SmN （3.14-11.68） and TbN/YbN （0.58-1.65）. The general parallelism of the REE patterns, suggest that all the granites were comagmatic in origin, while the high Eu/Eu＊ ratios （0.085-2.807; av.=0.9398） indicate high fo2 at the source. Similarly, irregular variations in LaN/YbN, CeN/YbN and Eu/Eu＊ ratios and REE abundances among the samples suggest behaviors that are related to mantle and crustal sources.     

Pressure–volume–temperature equation of state of andradite and grossular, by high-pressure and -temperature powder diffraction

 The thermoelastic parameters of natural andradite and grossular have been investigated by high-pressure and -temperature synchrotron X-ray powder diffraction, at ESRF, on the ID30 beamline. The P–V–T data have been fitted by Birch-Murnaghan-like EOSs, using both the approximated and the general form. We have obtained for andradite K ₀=158.0(±1.5) GPa, (dK/dT )₀=−0.020(3) GPa K⁻¹ and α₀=31.6(2) 10⁻⁶ K⁻¹, and for grossular K ₀=168.2(±1.7) GPa, (dK/dT)₀=−0.016(3) GPa K⁻¹ and α₀=27.8(2) 10⁻⁶ K⁻¹. Comparisons between the present issues and thermoelastic properties of garnets earlier determined are carried out. Received: 7 July 2000 / Accepted: 20 October 2000     

Fractional flow dimensions and hydraulic properties of a fracture-zone aquifer, Leppävirta, Finland

 In order to characterize the hydraulic properties of an aquifer in Finland comprising two subvertical fracture zones, observation-well responses were matched with generalized radial flow (GRF) type curves. The responses in six wells out of seven are consistent with the GRF model. The fractional flow dimensions (1–1.2 and 1.5) were determined by regression analysis of straight-line slopes and type-curve matching. In each test, the flow dimensions in the neighboring fracture zone range from 2–2.25. Comparisons of the late-time responses with the asymptotic GRF solution and the flow dimensions obtained by reversing the pumping and observation points suggest homogeneous hydraulic properties. Deviations in responses can be explained by flow-path tortuosity. After assessments of the extent of the flow and radial distances along the fracture system, hydraulic conductivities and storativities were determined from the results of the type-curve matching procedure. The obtained hydraulic conductivities are 1.3×10^–5 to 7.9×10^–5 m/s and 5.0×10^–6 to 2.5×10^–5 m/s for the western and the eastern fracture zones, respectively. The results were verified by applying them to analytical solutions for pumping wells. The calculated pumping-well responses are consistent with the observations. The analysis of flow dimension also enhances qualitative interpretations on the hydrogeology of fracture zones. Received, April 1997 · Revised, September 1997 · Accepted, May 1998     

Origin of the Maoduan Pb–Zn–Mo deposit, eastern Cathaysia Block, China: geological, geochronological, geochemical, and Sr–Nd–Pb–S isotopic constraints

The Maoduan Pb–Zn–Mo deposit is in hydrothermal veins with a pyrrhotite stage followed by a molybdenite and base metal stage. The Re–Os model ages of five molybdenite samples range from 138.6 ± 2.0 to 140.0 ± 1.9 Ma. Their isochron age is 137.7 ± 2.7 Ma. Laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS) zircon U–Pb dating of the nearby exposed Linggen granite porphyry gave a ²⁰⁶Pb/²³⁸U age of 152.2 ± 2.2 Ma and the hidden Maoduan monzogranite yielded a mean of 140.0 ± 1.6 Ma. These results suggest that the intrusion of the Maoduan monzogranite and Pb–Zn–Mo mineralization are contemporaneous. δ ³⁴S values of sulfide minerals range from 3.4‰ to 4.8‰, similar to magmatic sulfur. Four sulfide samples have ²⁰⁶Pb/²⁰⁴Pb = 18.252–18.432, ²⁰⁷Pb/²⁰⁴Pb = 15.609–15.779, and ²⁰⁸Pb/²⁰⁴Pb = 38.640–39.431, similar to the age-corrected data of the Maoduan monzogranite. These isotope data support a genetic relationship between the Pb–Zn–Mo mineralization and the Maoduan monzogranite and probably indicate a common deep source. The Maoduan monzogranite has geochemical features similar to highly fractionated I-type granites, such as high SiO₂ (73.7–75.2 wt.%) and alkalis (K₂O + Na₂O = 7.8–8.9 wt.%) and low FeO_t (0.8–1.3 wt.%), MgO (~0.3 wt.%), P₂O₅ (~0.03 wt.%), and TiO₂ (~0.2 wt.%). The granitic rocks are enriched in Rb, Th, and U but depleted in Ba, Sr, Nb, Ta, P, and Ti. REE patterns are characterized by marked negative Eu anomalies (Eu/Eu^* = 0.2–0.4). The Maoduan monzogranite, having (⁸⁷Sr/⁸⁶Sr) _t  = 0.7169 to 0.7170 and ε_Nd(t) = −13.8 to −13.7, was probably derived from mixing of partial melts from enriched mantle and the Paleoproterozoic Badu group in an extensional tectonic setting.     

High-temperature study and thermal expansion of phlogopite

 Unit-cell dimensions of a natural phlogopite from Pargas, Finland, have been determined in the temperature interval of 27–1050 °C by X-ray powder diffraction technique. Expansion rates vary discontinuously with temperature with a break at 412 °C. Below this temperature, the linear expansions (α) for a, b and c axis lengths are 3.74 × 10⁻⁵ K⁻¹, 1.09 × 10⁻⁵ K⁻¹, and 1.19 × 10⁻⁵ K⁻¹, respectively, and above that they are 0.86 × 10⁻⁵ K⁻¹, 0.80 × 10⁻⁵ K⁻¹, and 1.93 × 10⁻⁵ K⁻¹. The volume thermal expansion coefficients are 6.26 × 10⁻⁵ K⁻¹ and 3.71 × 10⁻⁵ K⁻¹ for low-temperature and high-temperature intervals, respectively. The observed kink in the rate of thermal expansions with temperature could be due to the different mode of structural changes. Thermogravimetric analysis of the sample indicates the oxidation of iron in the temperature range of 500–600 °C and dehydroxylation as well as decomposition of phlogopite in the temperature range of 900–1200 °C. Received: 8 September 1998 / Accepted: 28 February 2000     

Mineral inclusions in diamonds from the River Ranch kimberlite, Zimbabwe

More than 99% of mineral inclusions in diamonds from the River Ranch pipe in the Late Archean Limpopo Mobile Belt (Zimbabwe), are phases of harzburgitic paragenesis, namely olivine (Fo_92–93), orthopyroxene (Mg# = 93), G10 garnets and chromites. The diamond inclusion (DI) chemistry demonstrates a limited overlap with River Ranch kimberlite macrocrysts: the DI garnets are more Ca-undersaturated, and DI spinel and garnet are more Mg-rich. Most River Ranch diamond inclusions were equilibrated at T = 1080–1320 °C, P = 47–61 kbar, and f _O2 between IW and WM buffers. The P/T profile beneath the Limpopo Mobile Belt (LMB) is consistent with a paleo-heat flow of 41–42 mW/m², similar to calculations for Roberts Victor, but hotter than for the Finsch, Kimberley, Koffiefontein and Premier Mines. This is ascribed to the younger tectonothermal age of the LMB and its proximity to Late Archean oceans. Like diamond inclusions from all other kimberlites studied, the River Ranch DI have a lithospheric affinity and therefore indicate that an ancient, chemically depleted, thick (at least 200 km) mantle root existed beneath the Limpopo Mobile Belt 530–540 Ma ago. The mantle root might have developed beneath the continental Central Zone of the LMB as early as the Archean, and could be alien to the overthrust allochthonous sheet of the Limpopo Belt. Oxygen fugacity estimates for diamond inclusions at River Ranch are similar to other diamondiferous harzburgites beneath the Kaapvaal craton, indicating that the Kaapvaal mantle as a whole was well buffered and homogeneous with respect to f _O2 at the time of peridotitic diamond crystallization. Received: 11 January 1995 / Accepted: 10 June 1997