Comparison of several BHT correction methods: a case study on an Australian data set

Bottom-hole temperatures (BHT) from oil exploration provide useful constraints on the subsurface thermal regime, but they need to be corrected to obtain the equilibrium temperature. In this work, we introduce several BHT correction methods and compare them using a large Australian data set of more than 650 groups of multiple BHT measurements in about 300 oil exploration boreholes. Existing and suggested corrections are classified within a coherent framework, in which methods are divided into: line/cylinder source; instantaneous/continuous heat extraction; one/two component(s). Comparisons with reservoir test temperatures show that most of the corrections lead to reliable estimates of the formation equilibrium temperature within ±10 °C, but too few data exist to perform an intercomparison of the models based on this criterion. As expected, the Horner method diverges from its parent models for small elapsed times (or equivalently large radii). The mathematical expression of line source models suffers from an unphysical delay time that also restrains their domain of applicability. The model that takes into account the difference of thermal properties between circulating mud and surrounding rocks—that is the two-component model—is delicate to use because of its high complexity. For these reasons, our preferred correction methods are the cylindrical source models. We show that mud circulation time below 10 hr has a negligible effect. The cylindrical source models rely on one parameter depending on the thermal diffusivity and the borehole radius, which are poorly constrained, but the induced uncertainty on the extrapolations remains reasonably low.     

Analysis of central western Europe deformation using GPS and seismic data

The kinematic field of central western Europe is characterized by relatively small movements (around 1–2 mm/year) and diffuse seismicity with earthquakes occurring mostly in the shallow crust (within 15 km), prevalently concentrated along the Alps and the European Cenozoic Rift System (ECRIS). In order to study and constrain the current crustal kinematic field we reconstructed the velocity and the strain field using permanent GPS stations, belonging to different networks (AGNES, EUREF, REGAL, RGP). The 2D strain rate tensor has been calculated using the method of least-squares collocation. Our results show that the area of maximum compression is located along the Alpine chain, where maximum values of 7 ± 2 nstrain/year are found, while maximum extension is measured between the Armorican Massif and the Massif Central, where values of 4 ± 2 nstrain/year are reached.The earthquakes with M > 3.0, have been used to estimate the seismic strain rates, while the style of the seismic deformation was reconstructed from the fault plane solutions (FPS) available from the literature. Relatively high values of seismic strain rates (around 10 nstrain/year) are measured along the Alpine Chain and the ECRIS. Results obtained by geodetic and seismic data are quite in agreement and reflect the different tectonic evolution of the geological features characterizing the area of study. The orientation of the compressional geodetic and seismic strain axes are NW-SE in most of the area of study, on account of the action of plate boundary forces. A rotation of the same axes to N-S direction along the eastern Alps, possibly related to the Adria convergence, is found.     

High resolution simulation of recent Arctic and Antarctic stratospheric chemical ozone loss compared to observations

Simulations of polar ozone losses were performed using the three-dimensional high-resolution (1^∘ × 1^∘) chemical transport model MIMOSA-CHIM. Three Arctic winters 1999–2000, 2001–2002, 2002–2003 and three Antarctic winters 2001, 2002, and 2003 were considered for the study. The cumulative ozone loss in the Arctic winter 2002–2003 reached around 35% at 475 K inside the vortex, as compared to more than 60% in 1999–2000. During 1999–2000, denitrification induces a maximum of about 23% extra ozone loss at 475 K as compared to 17% in 2002–2003. Unlike these two colder Arctic winters, the 2001–2002 Arctic was warmer and did not experience much ozone loss. Sensitivity tests showed that the chosen resolution of 1^∘ × 1^∘ provides a better evaluation of ozone loss at the edge of the polar vortex in high solar zenith angle conditions. The simulation results for ozone, ClO, HNO₃, N₂O, and NO _y for winters 1999–2000 and 2002–2003 were compared with measurements on board ER-2 and Geophysica aircraft respectively. Sensitivity tests showed that increasing heating rates calculated by the model by 50% and doubling the PSC (Polar Stratospheric Clouds) particle density (from 5 × 10⁻³ to 10⁻² cm⁻³) refines the agreement with in situ ozone, N₂O and NO _y levels. In this configuration, simulated ClO levels are increased and are in better agreement with observations in January but are overestimated by about 20% in March. The use of the Burkholder et al. (1990) Cl₂O₂ absorption cross-sections slightly increases further ClO levels especially in high solar zenith angle conditions. Comparisons of the modelled ozone values with ozonesonde measurement in the Antarctic winter 2003 and with Polar Ozone and Aerosol Measurement III (POAM III) measurements in the Antarctic winters 2001 and 2002, shows that the simulations underestimate the ozone loss rate at the end of the ozone destruction period. A slightly better agreement is obtained with the use of Burkholder et al. (1990) Cl₂O₂ absorption cross-sections.     

The coupling of optimal economic growth and climate dynamics

In this paper, we study optimal economic growth programs coupled with climate change dynamics. The study is based on models derived from MERGE, a well established integrated assessment model (IAM). We discuss first the introduction in MERGE of a set of “tolerable window” constraints which limit both the temperature change and the rate of temperature change. These constraints, obtained from ensemble simulations performed with the Bern 2.5-D climate model, allow us to identity a domain intended to preserve the Atlantic thermohaline circulation. Next, we report on experiments where a two-way coupling is realized between the economic module of MERGE and an intermediate complexity “3-D-” climate model (C-GOLDSTEIN) which computes the changes in climate and mean temperature. The coupling is achieved through the implementation of an advanced “oracle based optimization technique” which permits the integration of information coming from the climate model during the search for the optimal economic growth path. Both cost-effectiveness and cost-benefit analysis modes are explored with this combined “meta-model” which we refer to as GOLDMERGE. Some perspectives on future implementations of these approaches in the context of “collaborative” or “community” integrated assessment modules are derived from the comparison of the different approaches.     

Utilisation de dépôts d'argile sur lame de verre pour leur étude en spectroscopie infrarouge

IR spectra of clay samples can be collected in transmission through oriented deposits onto glass slide and are similar to the IR spectra of the same samples obtained by traditional way (KBr pellets). As examples, it is shown that it is possible to differentiate dickite from kaolinite and smectites of various chemistries, as easily by both ways. It is thus possible to use oriented clay deposits onto glass slide for both XRD and IR studies and to easily take benefit of the complementarity of both techniques. To cite this article: S. Petit et al., C. R. Geoscience 335 (2003).     

Nouvelles données structurales et datations 40K–40Ar sur les roches métamorphiques de la région de Neyriz (zone de Sanandaj–Sirjan,Iran méridional). Leur intérêt dans le cadre du domaine néo-téthysien du Moyen-Orient

The metamorphic rocks of the Neyriz area (Sanandaj–Sirjan zone) represent a Palaeozoic sequence, the upper part of which being palaeontologically dated from the Carboniferous and the Permian. Field structural analysis of the whole sequence, detailed in laboratory by microstructural one and ⁴⁰K–⁴⁰Ar dating carried on separated minerals, lead to establish that the whole sequence, from gneisses to Permian rocks, has suffered a unique synmetamorphic deformation, of variable intensity, marked by a foliation. Isotopic ages measured on extracted amphiboles and micas, clustered in four groups between 300 and 60 Ma, show the successive stages of their slow exhumation, which ended by the end of the Cretaceous. To cite this article: R. Sheikholeslami et al., C. R. Geoscience 335 (2003).     

Characterization and correlation of Upper Jurassic (Oxfordian) bentonite deposits in the Paris Basin and the Subalpine Basin, France

Abstract Explosive volcanic activity is recorded in the Upper Jurassic of the Paris Basin and the Subalpine Basin of France by the identification of five bentonite horizons. These layers occur in Lower Oxfordian (cordatum ammonite zone) to Middle Oxfordian (plicatilis zone) clays and silty clays deposited in outer platform environments. In the Paris Basin, a thick bentonite (10–15 cm), identified in boreholes and in outcrop, is dominated by dioctahedral smectite (95%) with trace amounts of kaolinite, illite and chlorite. In contrast, five bentonites identified in the Subalpine Basin, where burial diagenesis and fluid circulation were more important, are composed of a mixture of kaolinite and regular or random illite/smectite mixed-layer clays in variable proportions, indicating a K-bentonite. In the Subalpine Basin, a 2–15 cm thick bentonite underlain by a layer affected by sulphate–carbonate mineralization can be correlated over 2000 km². Euhedral zircon, apatite and biotite crystals have been identified in all the bentonites. The geochemical composition of the bentonites in both basins is characterized by high concentrations of Hf, Nb, Pb, Ta, Th, Ti, U, Y, Zr and low concentrations of Cr, Cs and Rb. Biostratigraphical and geochemical data suggest that the thick bentonite in the Paris Basin correlates with the thickest bentonite in the Subalpine Basin, located 400 km to the south. These horizons indicate that significant explosive volcanic events occurred during the Middle Oxfordian and provide potential long-distance isochronous marker beds. Immobile element discrimination diagrams and rare-earth element characteristics indicate that the original ash compositions of the thickest bentonites correspond to a trachyandesitic source from a within-plate alkaline series that was probably related to North Atlantic rifting.     

Intercomparison of Boron Isotope and Concentration Measurements. Part II: Evaluation of Results

The Istituto di Geoscienze e Georisorse (IGG), on behalf and with the support of the International Atomic Energy Agency (IAEA), prepared eight geological materials (three natural waters and five rocks and minerals), intended for a blind interlaboratory comparison of measurements of boron isotopic composition and concentration. The materials were distributed to twenty seven laboratories - virtually all those performing geochemical boron isotope analyses in the world -which agreed to participate in the intercomparison exercise. Only fifteen laboratories, however, ultimately submitted the isotopic and/or concentration results they obtained on the intercomparison materials. The results demonstrate that interlaboratory reproducibility is not well reflected by the precision values reported by the individual laboratories and this observation holds true for both boron concentration and isotopic composition. The reasons for the discrepancies include fractionations due to the chemical matrix of materials, relative shift of the zero position on the δ¹¹B scale and a lack of well characterized materials for calibrating absolute boron content measurements. The intercomparison materials are now available at the IAEA (solid materials) and IGG (waters) for future distribution.     

The Hasan Dagi stratovolcano (Central Anatolia, Turkey): evolution from calc-alkaline to alkaline magmatism in a collision zone

The Hasan Dagi volcano is one of the two large Plio-Quaternary volcanoes in Cappadocia (Central Anatolia, Turkey). Three stages of edifice construction have been identified for this volcano: Paleovolcano, Mesovolcano and Neovolcano. Most samples from Hasan Dagi volcano are calc-alkaline and define an almost complete trend from basaltic andesite to rhyolite. However, the more recent (Neovolcano) mafic samples are alkaline basalts. The mineralogical and geochemical characteristics of the oldest lavas (Keçikalesi (13 Ma) and Paleo-Hasan Dagi (7 Ma)) are significantly different from those of the younger lavas (Meso- and Neo-Hasan Dagi (<1 Ma)). Calcic plagioclase and pigeonite are typically observed in these older lavas. The Paleovolcano basalts are depleted in alkalis and display a tholeiitic tendency whereas the differentiated lavas are depleted in Na₂O but enriched in K₂O compared to younger lavas. There is an evolution through time towards higher TiO₂, Fe₂O₃^*, MgO, Na₂O and K₂O and lower Al₂O₃ and SiO₂ which is reflected in the basalt compositions. All the basalts display multi-element patterns typical of continental margin magmas with a significant enrichment in LILE (K, Rb, Ba and Th) and LREE and strong (Paleovolcano) to moderate (Meso- and Neovolcano) negative Nb, Zr and Ti anomalies. However, the younger basalts are the most enriched in incompatible elements, in agreement with their alkaline affinities and do not systematically display negative HFSE anomalies. REE data suggest an hydrous amphibole-bearing crystallization history for both Meso- and Neovolcano lavas. The distinction between the older and younger lavas is also apparent in trace element ratios such as Nb/Y, Ti/Y and Th/Y. These ratios indicate the role of a subducted component±crustal contamination in the genesis of the Hasan Dagi lavas, particularly for the oldest lavas (Keçikalesi and Paleo-Hasan Dagi). The decreasing influence of this component through time, over the last 6–7 m.y., has been accompanied by an increasing contribution of melt-enriched lithosphere. Although the range of variation of Sr, Nd and Pb isotopic ratios is small (0.70457–0.70515; 0.51262–0.51273; 18.80–18.94; 15.64–15.69; 38.87–39.10), it also reflects the evolution of the magma sources through time. Indeed, the youngest (Neovolcano) and most primitive basalts display significantly lower ⁸⁷Sr/⁸⁶Sr than the Paleo- and Mesovolcano basalts, whereas the Mesovolcano basalts display more radiogenic Pb than Paleovolcano samples. Magma mixing processes between initially heterogeneous and/or variably contaminated magmas may account for the genesis of the less differentiated and intermediate lavas (48–57% SiO₂). Meso- and Neovolcano differentiated lavas (60–68% SiO₂) are either derived from the analyzed basalts or from more primitive and more depleted magmas by fractional crystallization±some crustal contamination (AFC). Furthermore, the highly differentiated samples (72–75% SiO₂) are not strongly contaminated. The strong calc-alkaline character of Hasan Dagi lavas, in the absence of contemporaneous subduction, must reflect the heritage of the early subduction of the Afro–Arabian plate under the Eurasian plate. The evolution towards alkaline compositions through time is clearly related to the development of extensional tectonics in Central Anatolia in the Late Miocene.     

Effect of Complexing Agents (EDTA and ATMP) on the Remobilization of Heavy Metals from a Polluted River Sediment

EDTA and ATMP are synthetic chelatingagents used in many industrial applications. Becauseof their low biodegradability, they are present inmany aquatic systems where their strong complexingcapacities can modify the transfer of toxic heavymetals between solid and liquid phases.The remobilization of Cu, Pb and Cd from pollutedsuperficial river sediments in the presence of twosynthetic organic ligands (a polycarboxylate, EDTA,and a phosphonate, ATMP) is investigated as a functionof reaction time to determine the time necessary toreach dissolutive equilibrium, and as a fonction ofthe complexing agent concentration.The dissolutive equilibrium times are similar (10 to20 hours) and the solubilizing power of ATMP isslightly lower than that of EDTA. The curves ofsolubilization of metals in the presence of thecomplexing agents present a shape of a chelateconcentration desorption edge, for a chelate rangebetween 10^-5 to 10^-3 mol L^-1, thesolubilised metals increasing dramatically. Thisremobilization is the result of a competition betweenassociation with the solid phase and dissolvedcomplexation. With EDTA, the desorption order ofmetals seems to be imposed by association with thesolid. For ATMP is it the chelate complexing strengthwhich determines this order.