热液矿床成矿作用研究新途径:数值模拟

热液矿床是最具有工业价值的矿床类型之一,传统矿床学研究在成岩成矿作用、成矿机制和成矿规律等方面已有丰硕成果,但仍存在一些尚无法有效回答的科学问题。近年来,数值模拟方法逐步应用于热液矿床的研究,已成为成矿过程定量研究的新手段。本文系统总结了近年来热液矿床成矿过程数值模拟研究的成果和进展,并对存在问题和发展趋势进行了分析和展望,以期进一步推动这一交叉学科领域的研究。     

Effect of water on the fluorine and chlorine partitioning behavior between olivine and silicate melt

Halogens show a range from moderate (F) to highly (Cl, Br, I) volatile and incompatible behavior, which makes them excellent tracers for volatile transport processes in the Earth’s mantle. Experimentally determined fluorine and chlorine partitioning data between mantle minerals and silicate melt enable us to estimate Mid Ocean Ridge Basalt (MORB) and Ocean Island Basalt (OIB) source region concentrations for these elements. This study investigates the effect of varying small amounts of water on the fluorine and chlorine partitioning behavior at 1280?°C and 0.3 GPa between olivine and silicate melt in the Fe-free CMAS+F–Cl–Br–I–H₂O model system. Results show that, within the uncertainty of the analyses, water has no effect on the chlorine partitioning behavior for bulk water contents ranging from 0.03 (2) wt% H₂O (D_Cl ^ol/melt = 1.6?±?0.9 × 10^?4) to 0.33 (6) wt% H₂O (D_Cl ^ol/melt = 2.2?±?1.1 × 10^?4). Consequently, with the effect of pressure being negligible in the uppermost mantle (Joachim et al. Chem Geol 416:65–78, 2015), temperature is the only parameter that needs to be considered for the determination of chlorine partition coefficients between olivine and melt at least in the simplified iron-free CMAS+F–Cl–Br–I–H₂O system. In contrast, the fluorine partition coefficient increases linearly in this range and may be described at 1280?°C and 0.3 GPa with (R ²?=?0.99): \(D_{F}^{\text{ol/melt}}\ =\ 3.6\pm 0.4\ \times \ {{10}^{-3}}\ \times \ {{X}_{{{\text{H}}_{\text{2}}}\text{O}}}\left( \text{wt }\!\!\%\!\!\text{ } \right)\ +\ 6\ \pm \ 0.4\times \,{{10}^{-4}}\). The observed fluorine partitioning behavior supports the theory suggested by Crépisson et al. (Earth Planet Sci Lett 390:287–295, 2014) that fluorine and water are incorporated as clumped OH/F defects in the olivine structure. Results of this study further suggest that fluorine concentration estimates in OIB source regions are at least 10% lower than previously expected (Joachim et al. Chem Geol 416:65–78, 2015), implying that consideration of the effect of water on the fluorine partitioning behavior between Earth’s mantle minerals and silicate melt is vital for a correct estimation of fluorine abundances in OIB source regions. Estimates for MORB source fluorine concentrations as well as chlorine abundances in both mantle source regions are within uncertainty not affected by the presence of water.     

New 40Ar/39Ar constraints for the “Grande Nappe”: The largest rhyolitic eruption from the Mont-Dore Massif (French Massif Central)

Since the 1960s, an early explosive activity in the Mont-Dore Massif is associated with a major pyroclastic rhyolitic eruption (5–7 km³) known as the “Grande Nappe” (GN). This event, linked to the formation of a 6-km-diameter cryptic caldera named “Haute Dordogne”, was before our investigation dated by ⁴⁰Ar/³⁹Ar at 3.07 ± 0.04 Ma. Our new single-crystal laser fusion ⁴⁰Ar/³⁹Ar dates obtained on two outcrops of the GN (Rochefort-Montagne and Ludières) questioned several hypotheses made concerning this “landmark” event of the Mont-Dore Massif history. We demonstrate that: (1) the GN rhyolitic eruption has occurred much later than previously estimated (i.e. 2.77 ± 0.02–0.07 Ma full external uncertainties); (2) the correlation made between the Vendeix rhyolitic complexes (intra-caldera position) dated back to 2.74 ± 0.04 Ma and the GN is proposed; (3) xenocryst contamination could be very high (i.e. 70% for the Rochefort-Montagne GN outcrop) and explains the noticeable older age obtained previously; (4) a link between the GN eruption and the formation of a caldera is questionable; the hypothesis of a northward-oriented blast channeled eastward toward the paleo-Allier River is thus proposed.     

Monitoring viscosity changes from time‐lapse seismic attenuation: case study from a heavy oil reservoir

Heating heavy oil reservoirs is a common method for reducing the high viscosity of heavy oil and thus increasing the recovery factor. Monitoring of these viscosity changes in the reservoir is essential for delineating the heated region and controlling production. In this study, we present an approach for estimating viscosity changes in a heavy oil reservoir. The approach consists of three steps: measuring seismic wave attenuation between reflections from above and below the reservoir, constructing time‐lapse Q and Q^?1 factor maps, and interpreting these maps using Kelvin–Voigt and Maxwell viscoelastic models. We use a 4D relative spectrum method to measure changes in attenuation. The method is tested with synthetic seismic data that are noise free and data with additive Gaussian noise to show the robustness and the accuracy of the estimates of the Q‐factor. The results of the application of the method to a field data set exhibit alignment of high attenuation zones along the steam‐injection wells, and indicate that temperature dependent viscosity changes in the heavy oil reservoir can be explained by the Kelvin–Voigt model.     

Detection of Pharmaceuticals and Personal Care Products in Agricultural Soils Receiving Biosolids Application

Increasing research has suggested that biosolids generated from municipal wastewater treatment can be a major sink for many pharmaceuticals and personal care products (PPCPs) and their land application potentially introduces these contaminants into the terrestrial and aquatic environments. In this study, methods were developed for the analysis of 14 PPCPs in biosolids and soils using pressurized liquid extraction, solid phase extraction and liquid chromatography‐tandem mass spectrometry. Recoveries were over 50% for all analytes except diphenhydramine (?30%) in soils. Soil properties or type of biosolids showed minor effects on method recoveries. Estimated method limits of quantification (LOQ) range from 0.1–15 ng g^–1 for soil and 0.3–27 ng g^–1 for biosolids. A field study utilizing the methods revealed that other than carbamazepine‐10,11‐epoxide, all targeted compounds were detected in biosolids. Diphenhydramine, fluoxetine, triclosan and triclocarban were detected up to the μg g^–1 range with the highest concentration of 23 μg g^–1 for triclocarban. Seven of the PCCPs found in biosolids were also detected in agricultural soils amended with these biosolids and several (carbamazepine, diphenhydramine, and triclocarban) appeared to be persistent in soils. Triclocarban was also found most abundant in soils with the highest average concentration of 0.2 μg g^–1 while the rest of compounds were in the lower ng g^–1 range. Generally, the concentrations found on the fields were 2–3 degrees of magnitude lower than in the biosolids, which is likely to be due to dilution, degradation and leaching processes.     

Pressure and equilibrium in deforming rocks

It is widely proposed that tectonic pressure (the difference between the mean stress and the pressure arising from a lithostatic load) is large, and has a significant influence on mineral phase equilibria in deforming metamorphic rocks. The implication/assertion is that the mean stress is equivalent to the thermodynamic pressure which characterizes mineral phase equilibria and is a measure of how the energy changes as the volume changes. We distinguish two useful thermodynamic pressures. The first is an equilibrium thermodynamic pressure, characteristic of non‐dissipative systems and related directly to equilibrium values of the chemical potentials that define stable, equilibrium phase assemblages. The second is a non‐equilibrium thermodynamic pressure characteristic of dissipative systems with zero net entropy production and related to non‐equilibrium chemical potentials that define stable non‐equilibrium phase assemblages. In many dissipative metamorphic systems discussed in the literature, the concepts of thermodynamic pressure and chemical potential are not usefully defined because the system is not at equilibrium and/or no volume change is involved in the deformation. The conclusion of this note is that the influence of tectonic pressure on phase equilibria is minor. The role of tectonic pressure is an important issue but is only relevant to phase equilibrium when an equilibrium thermodynamic pressure can be defined; in such cases, the influence of tectonic pressure is small compared to many proposals in the literature. Except for elastic deformations, the mean stress is not useful in discussing mineral phase equilibrium.     

Stable carbon and oxygen isotope record of central Lake Erie sediments

Stable carbon and oxygen isotope data from mollusc aragonite extracted from sediment cores provide new information on the origin and history of sedimentation in the southwestern area of the central basin of Lake Erie. Sediments infilling the Sandusky subbasin consist of three lithologic units overlying glacial deposits. The lowest of these is a soft gray mud overlain by a shell hash layer containing Sphaerium striatinum fragments. A fluid mud unit caps the shell hash layer and extends upwards to the sediment-water interface. New stable isotope data suggest that the soft gray mud unit is of postglacial, rather than proglacial, origin. These data also suggest that the shell hash layer was derived from erosional winnowing of the underlying soft gray mud layer. This winnowing event may have occurred as a result of the Nipissing flood.The Pelee-Lorain moraine, which forms the eastern boundary of the Sandusky subbasin, is an elevated area of till capped by a sand deposit that originated as a beach. The presence of both the shell hash layer and relict beach deposit strengthens the interpretation that the Nipissing flood was a critical event in the development of the southwestern area of the central basin of Lake Erie. This event, which returned drainage from the upper lakes to the Lake Erie basin, was a dominant influence on regional stratigraphy, bathymetry, and depositional setting.     

Coulomb constitutive laws for friction: Contrasts in frictional behavior for distributed and localized shear

We describe slip-rate dependent friction laws based on the Coulomb failure criteria. Frictional rate dependence is attributed to a rate dependence of cohesionc and friction angle . We show that differences in the stress states developed during sliding result in different Coulomb friction laws for distributed shear within a thick gouge layer versus localized shear within a narrow shear band or between bare rock surfaces. For shear within gouge, shear strength is given by =c cos + _n sin, whereas for shear between bare rock surfaces the shear strength is =c cos + _n tan, where and _n are shear and normal stress, respectively. In the context of rate-dependent Coulomb friction laws, these differences mean that for a given material and rate dependence of the Coulomb parameters, pervasive shear may exhibit velocity strengthening frictional behavior while localized shear exhibits velocity weakening behavior. We derive from experimental data the slip-rate dependence and evolution ofc and for distributed and localized shear. The data show a positive rate dependence for distributed shear and a negative rate dependence for localized shear, indicating that the rate dependence ofc and are not the same for distributed and localized shear, even after accounting for differences in stress state. Our analysis is consistent with the well-known association of instability with shear localization in simulated fault gouge and the observation that bare rock surfaces exhibit predominantly velocity weakening frictional behavior whereas simulated fault gouge exhibits velocity strengthening followed by a transition to velocity weakening with increasing displacement. Natural faults also exhibit displacement dependent frictional behavior and thus the results may prove useful in understanding the seismic evolution of faulting.     

The northern Ossa-Morena Cadomian batholith (Iberian Massif): magmatic arc origin and early evolution

Diorites and related rocks in the Mérida area of northern Ossa-Morena (SW Iberia) are intrusive into Precambrian metavolcanic and metasedimentary sequences. Cumulate products from the H₂O-rich magmas are amphibole-rich gabbros to hornblendites. Major and trace element compositions, including Sr and Nd isotope data, allow the definition of a calc-alkaline series likely formed in relation to an immature arc setting. Crystallization of the intrusives has been established between ca. 570 and 580 Ma by U-Pb dating of constituent zircons. Garnet growth in dioritic rocks reflects a tectono–thermal overprint dated by Sm-Nd internal isochrons at around 555 Ma. Older Sm-Nd and Lu-Hf results between ca. 593 and 637 Ma on the same rocks suggest an earlier stage of regional metamorphism within the arc environment. The northern Ossa-Morena composite batholith and related metamorphic units have been tectonized and dismembered in the course of subsequent low-grade events during final stages of the Cadomian orogeny and the Variscan cycle. The units studied represent a well preserved segment of the arc region that evolved in Neoproterozoic times along the western border of Gondwana to conform the Cadomian–Avalonian basement of the Hercynian realm.This revised version was published in October 2004 with corrections to the names and order of the authors.