Eclogite thermobarometry: The consistency between conventional thermobarometry and forward phase-equilibrium modelling

Eclogite thermobarometry is crucial for constraining the depths and temperatures to which oceanic and continental crust subduct. However, obtaining the pressure and temperature (P–T) conditions of eclogites is complex as they commonly display high-variance mineral assemblages, and the mineral compositions only vary slightly with P–T. In this contribution, we present a comparison between two independent and commonly used thermobarometric approaches for eclogites: conventional thermobarometry and forward phase-equilibrium modelling. We assess how consistent the thermobarometric calculations are using the garnet–clinopyroxene–phengite barometer and garnet–clinopyroxene thermometer with predictions from forward modelling (i.e. comparing the relative differences between approaches). Our results show that the overall mismatch in methods is typically ±0.2–0.3 GPa and ±29–42°C although differences as large as 80°C and 0.7 GPa are possible for a few narrow ranges of P–T conditions in the forward models. Such mismatch is interpreted as the relative differences among methods, and not as absolute uncertainties or accuracies for either method. For most of the investigated P–T conditions, the relatively minor differences between methods means that the choice in thermobarometric method itself is less important for geological interpretation than careful sample characterization and petrographic interpretation for deriving P–T from eclogites. Although thermobarometry is known to be sensitive to the assumed X_Fe³⁺ of a rock (or mineral), the relative differences between methods are not particularly sensitive to the choice of bulk-rock X_Fe³⁺, except at high temperatures (>650°C, amphibole absent) and for very large differences in assumed X_Fe³⁺ (0–0.5). We find that the most important difference between approaches is the activity–composition (a–x) relations, as opposed to the end-member thermodynamic data or other aspects of experimental calibration. When equivalent a–x relations are used in the conventional barometer, P calculations are nearly identical to phase-equilibrium models (ΔP < 0.1). To further assess the implications of these results for real rocks, we also evaluate common mathematical optimizations of reaction constants used for obtaining the maximum P–T with conventional thermobarometric approaches (e.g. using the highest aGrs² × aPrp in garnet and Si content in phengite, and the lowest aDi in clinopyroxene). These approaches should be used with caution, because they may not represent the compositions of equilibrium mineral assemblages at eclogite facies conditions and therefore systematically bias P–T calculations. Assuming method accuracy, geological meaningful P_max at a typical eclogite facies temperature of ~660°C will be obtained by using the greatest aDi, aCel, and aPrp and lowest aGrs and aMs; garnet and clinopyroxene with the lowest Fe²⁺/Mg ratios may yield geological meaningful T_max at a typical eclogite facies pressure of 2.5 GPa.     

Ethnic Place Identity Within A Parisian Neighborhood

The Goutte d'Or neighborhood of Paris has long been seen as an impoverished yet colorful quarter of the city. In recent years, it has also developed a reputation as a center for immigrants. Three populations now share the neighborhood: the European French, the Maghrebi immigrant population and their children, and the newest population of migrants from former Francophone colonies. This paper examines how the Goutte d'Or's social and ethnic identity is revealed in the built environment and in its social and economic activity. This neighborhood demonstrates how different communities and place identities coexist. It also shows how global networks of migration, exchange, and visits infuse local places. The Goutte d'Or distills many aspects of French immigrant identities, providing an example for increasing numbers of immigrant neighborhoods in France and across the world.     

Synthetic coral-reef terraces and variations of Quaternary sea level

Coral-reef terraces on tectonically elevated coasts provide important data on variations of sea level. In this paper a systematic approach is presented for the construction of synthetic reef terraces. Idealized reef terraces are generated assuming a constant rate of tectonic uplift ν and a harmonic variation of sea level with amplitudeh_wo and period τ. It is further assumed that coral terraces grow only during periods of relative sea-level rise. The behavior is governed by the parameter u = τυ/2π h_wo. If u > 1 no reef terraces form. If 0.217 < u < 1 a discontinuous series of terraces form on the uplifting basement. If 0 < u < 0.217 the basement is covered with coral and a series of coral terraces forms on this coral. Uplifted coral-reef terraces on the Huon Peninsula, Papua New Guinea, compare reasonably well with terraces generated by our model with τ = 20,000 years andh_wo = 21.4 m. We also compare the observed terraces with model terraces using the paleo sea level inferred from¹⁸O/¹⁶O data. The observed terraces are qualitatively similar to those predicted using the oxygen isotope data, however, the model terraces require larger relative sea-level variations with a near 20,000 year period than those inferred directly from the¹⁸O/¹⁶O data.     

Inverse Analysis of Deep Excavation Using Differential Evolution Algorithm

This paper presents the applications of the differential evolution (DE) algorithm in back analysis of soil parameters for deep excavation problems. A computer code, named Python‐based DE, is developed and incorporated into the commercial finite element software ABAQUS, with a parallel computing technique to run an FE analysis for all trail vectors of one generation in DE in multiple cores of a cluster, which dramatically reduces the computational time. A synthetic case and a well‐instrumented real case, that is, the Taipei National Enterprise Center (TNEC) project, are used to demonstrate the capability of the proposed back‐analysis procedure. Results show that multiple soil parameters are well identified by back analysis using a DE optimization algorithm for highly nonlinear problems. For the synthetic excavation case, the back‐analyzed parameters are basically identical to the input parameters that are used to generate synthetic response of wall deflection. For the TNEC case with a total of nine parameters to be back analyzed, the relative errors of wall deflection for the last three stages are 2.2, 1.1, and 1.0%, respectively. Robustness of the back‐estimated parameters is further illustrated by a forward prediction. The wall deflection in the subsequent stages can be satisfactorily predicted using the back‐analyzed soil parameters at early stages. Copyright © 2014 John Wiley & Sons, Ltd.     

Potential environmental issues of CO2 storage in deep saline aquifers: Geochemical results from the Frio-I Brine Pilot test,Texas, USA

Sedimentary basins in general, and deep saline aquifers in particular, are being investigated as possible repositories for large volumes of anthropogenic CO₂ that must be sequestered to mitigate global warming and related climate changes. To investigate the potential for the long-term storage of CO₂ in such aquifers, 1600 t of CO₂ were injected at 1500 m depth into a 24-m-thick “C” sandstone unit of the Frio Formation, a regional aquifer in the US Gulf Coast. Fluid samples obtained before CO₂ injection from the injection well and an observation well 30 m updip showed a Na–Ca–Cl type brine with ∼93,000 mg/L TDS at saturation with CH₄ at reservoir conditions; gas analyses showed that CH₄ comprised ∼95% of dissolved gas, but CO₂ was low at 0.3%. Following CO₂ breakthrough, 51 h after injection, samples showed sharp drops in pH (6.5–5.7), pronounced increases in alkalinity (100–3000 mg/L as HCO₃) and in Fe (30–1100 mg/L), a slug of very high DOC values, and significant shifts in the isotopic compositions of H₂O, DIC, and CH₄. These data, coupled with geochemical modeling, indicate corrosion of pipe and well casing as well as rapid dissolution of minerals, especially calcite and iron oxyhydroxides, both caused by lowered pH (initially ∼3.0 at subsurface conditions) of the brine in contact with supercritical CO₂.     

Crystal fractionation and partial melting in the petrogenesis of a Proterozoic high-MgO volcanic suite,Ungava, Québec

There is little concensus on the relative importance of crystal fractionation and differential partial melting to the chemical diversity observed within most types of volcanic suites. A resolution to this controversy is best sought in suites containing high MgO lavas such as the Chukotat volcanics of the Proterozoic Cape Smith foldbelt, Ungava, Quebec. The succession of this volcanic suite consists of repetitive sequences, each beginning with olivine-phyric basalt (19-12 wt% MgO), grading upwards to pyroxene-phyric basalt (12-8 wt% MgO) and then, in later sequences, to plagioclase-phyric basalt (7-4 wt% MgO). Only the olivine-phyric basalts have compositions capable of equilibrating with the upper mantle and are believed to represent parental magmas for the suite. The pyroxene-phyric and plagioclase-phyric basalts represent magmas derived from these parents by the crystal fractionation of olivine, with minor chromite, clinopyroxene and plagioclase. The order of extrusion in each volcanic sequence is interpreted to reflect a density effect in which successively lighter, more evolved magmas are erupted as hydrostatic pressure wanes. The pyroxene-phyric basalts appear to have evolved at high levels in the active part of the conduit system as the eruption of their parents was in progress. The plagioclase-phyric basalts may represent residual liquids expelled from isolated reservoirs along the crust-mantle interface during the late stages of volcanic activity.A positive correlation between FeO and MgO in the early, most basic olivine-phyric basalts is interpreted to reflect progressive adiabatic partial melting in the upper mantle. Although this complicates the chemistry, it is not a significant factor in the compositional diversification of the volcanic suite. The preservation of such compositional melting effects, however, suggests that the most basic olivine-phyric basalts represent primitive magmas. The trace element characteristics of these magmas, and their derivatives, indicate that the mantle source for the Chukotat volcanics had experienced a previous melting event.     

The biogeochemical sulfur cycle in the marine boundary layer over the Northeast Pacific Ocean

The major components of the marine boundary layer biogeochemical sulfur cycle were measured simultaneously onshore and off the coast of Washington State, U.S.A. during May 1987. Seawater dimethylsulfide (DMS) concentrations on the continental shelf were strongly influenced by coastal upwelling. Concentration further offshore were typical of summer values (2.2 nmol/L) at this latitude. Although seawater DMS concentrations were high on the biologically productive continental shelf (2–12 nmol/L), this region had no measurable effect on atmospheric DMS concentrations. Atmospheric DMS concentrations (0.1–12 nmol/m³), however, were extremely dependent upon wind speed and boundary layer height. Although there appeared to be an appreciable input of non-sea-salt sulfate to the marine boundary layer from the free troposphere, the local flux of DMS from the ocean to the atmosphere was sufficient to balance the remainder of the sulfur budget.     

Estimation of the effect of protective piers on wind wave parameters in the Sevastopol Bay

An effectiveness of the storm wave attenuation by protective piers in the Sevastopol Bay of the Black Sea is studied on the basis of numerical simulation using the SWAN spectral model. Analyzed are the parameters of waves generated by winds of four main directions as well as by the southern cyclone during the storm on November 11, 2007. It is obtained that waves from the northwest part of the Black Sea penetrate most intensively into the Sevastopol Bay in case of western wind and, to a lesser degree, in case of northern and southern winds. A protective effect of the piers is observed in the west part of the bay only and the wave attenuation near the southern coast is more significant than near the northern one. The area of the southern coast directly behind the southern pier is completely protected from the storm waves and, as moving away from the pier, the danger of intensive wave effect on the coast is kept.     

Modelling of local velocity anomalies: a cookbook

The determination of small-scale velocity anomalies (from tens to a few hundreds of metres) is a major problem in seismic exploration. The impact of such anomalies on a structural interpretation can be dramatic and conventional techniques such as tomographic inversion or migration velocity analysis are powerless to resolve the ambiguity between structural and velocity origins of anomalies. We propose an alternative approach based on stochastic modelling of numerous anomalies until a set of models is found which can explain the real data. This technique attempts to include as much a priori geological information as possible. It aims at providing the interpreter with a set of velocity anomalies which could possibly be responsible for the structural response. The interpreter can then choose one or several preferred models and pursue a more sophisticated analysis. The class of retained models are all equivalent in terms of data and therefore represent the uncertainty in the model space. The procedure emulates the real processing sequence using a simplified scheme. Essentially, the technique consists of five steps: 1 Interpretation of a structural anomaly in terms of a velocity anomaly with its possible variations in terms of position, size and amplitude. 2 Drawing a model by choosing the parameters of the anomaly within the acceptable range. 3 Modelling the traveltimes in this model and producing the imaging of the reflected interface. 4 Comparing the synthetic data with the real data and keeping the model if it lies within the data uncertainty range. 5 Iterate from step 2. In order to avoid the high computational cost inherent in using statistical determinations, simplistic assumptions have been made: ? The anomaly is embedded in a homogeneous medium: we assume that the refraction and the time shift due to the anomaly have a first-order effect compared with ray bending in the intermediate layers. ? We model only the zero-offset rays and therefore we restrict ourselves to structural problems. ? We simulate time migration and so address only models of limited structural complexity. These approximations are justified in a synthetic model which includes strong lateral velocity variations, by comparing the result of a full processing sequence (prestack modelling, stack and depth migration) with the simplified processing. This model is then used in a blind test on the inversion scheme.