Peak metamorphic temperatures in type 6 ordinary chondrites: An evaluation of pyroxene and plagioclase geothermometry

Abstract— Quantifying the peak temperatures achieved during metamorphism is critical for understanding the thermal histories of ordinary chondrite parent bodies. Various geothermometers have been used to estimate equilibration temperatures for chondrites of the highest metamorphic grade (type 6), but results are inconsistent and span hundreds of degrees. Because different geothermometers and calibration models were used with different meteorites, it is unclear whether variations in peak temperatures represent actual ranges of metamorphic conditions within type 6 chondrites or differences in model calibrations. We addressed this problem by performing twopyroxene geothermometry, using QUILF95, on the same type 6 chondrites for which peak temperatures were estimated using the plagioclase geothermometer (Nakamuta and Motomura 1999). We also calculated temperatures for published pyroxene analyses from other type 6 H, L, and LL chondrites to determine the most representative peak metamorphic temperatures for ordinary chondrites. Pyroxenes record a narrow, overlapping range of temperatures in H6 (865–926 °C), L6 (812–934 °C), and LL6 (874–945 °C) chondrites. Plagioclase temperature estimates are 96–179 °C lower than pyroxenes in the same type 6 meteorites. Plagioclase estimates may not reflect peak metamorphic temperatures because chondrule glass probably recrystallized to plagioclase prior to reaching the metamorphic peak. The average temperature for H, L, and LL chondrites (～900 °C), which agrees with previously published oxygen isotope geothermometry, is at least 50 °C lower than the peak temperatures used in current asteroid thermal evolution models. This difference may require minor adjustments to thermal model calculations.     

Palaeoenvironmental significance of Late Permian palaeosols in the South‐Eastern Iberian Ranges,Spain

The Late Permian (Wuchiapingian) Alcotas Formation in the SE Iberian Ranges consists of one red alluvial succession where abundant soil profiles developed. Detailed petrographical and sedimentological studies in seven sections of the Alcotas Formation allow six different types of palaeosols, with distinctive characteristics and different palaeogeographical distribution, to be distinguished throughout the South‐eastern Iberian Basin. These characteristics are, in turn, related to topographic, climatic and tectonic controls. The vertical distribution of the palaeosols is used to differentiate the formation in three parts from bottom to top showing both drastic and gradual vertical upwards palaeoenvironmental changes in the sections. Reconstruction of palaeoenvironmental conditions based on palaeosols provides evidence for understanding the events that occurred during the Late Permian, some few millions of years before the well‐known Permian‐Triassic global crisis.     

Methane‐related microbial gypsum calcitization in stromatolites of a marine evaporative setting (Münder Formation,Upper Jurassic,Hils Syncline,north Germany)

Fossil stromatolites may reveal information about their hydrochemical palaeoenvironment, provided that assignment to a specific microbial community and a corresponding biogeochemical mechanism of formation can be made. Tithonian stromatolites of the Münder Formation at Thüste, north Germany, have traditionally been considered as formed by intertidal cyanobacterial communities. However, thin sections of the stromatolites show elongated angular traces of former gypsum crystals in a dense arrangement, but no algal or cyanobacterial filament traces. Moreover, high Fe²⁺ and Mn²⁺ contents, oxygen‐isotope and sulphur‐isotope ratios of carbonate‐bound sulphates, and sulphurized hydrocarbon biomarkers of the stromatolitic carbonate indicate that CaCO₃ precipitation occurred near the oxic–anoxic interface as a result of intensive bacterial sulphur cycling rather than photosynthetic activity. Furthermore, anaerobic oxidation of methane by Archaea may have driven CaCO₃ precipitation in deeper parts of the biofilm community, as reflected by high concentrations of squalane with a strongly negative δ¹³C in conjunction with evaporite pseudomorphs showing extremely low δ¹³C_Carb ratios. Consequently, the Thüste stromatolites are now interpreted as having initially formed by gypsum impregnation of biofilms. Subsequently, early Mg‐calcitic calcitization within the biofilms occurred because of combined bacterial iron, manganese and sulphate reduction, with an increasing contribution of anaerobic oxidation of methane with depth. This model plausibly explains the prominent preservation of signals derived from oxygen‐independent metabolic pathways, whereas virtually no geochemical record exists for an aerobic community that may, nevertheless, have prevailed at the stromatolite surface. Photic‐zone stromatolites with a prominent signal of anaerobic oxidation of methane may be common in, and indicative of, oxygen‐depleted sulphate‐bearing environments with high rates of methane production, conditions that possibly were fulfilled at the Archaean to Proterozoic transition.     

He–Ar and Nd–Sr isotopic compositions of late Pleistocene felsic plutonic back arc basin rocks from Ulleungdo volcanic island, South Korea: Implications for the genesis of young plutonic rocks in a back arc basin

We report analyses of noble gases and Nd–Sr isotopes in mineral separates and whole rocks of late Pleistocene (< 0.2 Ma) monzonites from Ulleungdo, South Korea, a volcanic island within the back arc basin of the Japan island arc. A Rb–Sr mineral isochron age for the monzonites is 0.12 ± 0.01 Ma. K–Ar biotite ages from the same samples gave relatively concordant ages of 0.19 ± 0.01and 0.22 ± 0.01 Ma. ⁴⁰Ar/³⁹Ar yields a similar age of 0.29 ± 0.09 Ma. Geochemical characteristics of the felsic plutonic rocks, which are silica oversaturated alkali felsic rocks (av., 12.5 wt% in K₂O + Na₂O), are similar to those of 30 alkali volcanics from Ulleungdo in terms of concentrations of major, trace and REE elements. The initial Nd–Sr isotopic ratios of the monzonites (⁸⁷Sr/⁸⁶Sr = 0.70454–0.71264, ¹⁴³Nd/¹⁴⁴Nd = 0.512528–0.512577) are comparable with those of the alkali volcanics (⁸⁷Sr/⁸⁶Sr = 0.70466–0.70892, ¹⁴³Nd/¹⁴⁴Nd = 0.512521–0.512615) erupted in Stage 3 of Ulleungdo volcanism (0.24–0.47 Ma). The high initial ⁸⁷Sr/⁸⁶Sr values of the monzonites imply that seawater and crustally contaminated pre-existing trachytes may have been melted or assimilated during differentiation of the alkali basaltic magma.A mantle helium component (³He/⁴He ratio of up to 6.5 RA) associated with excess argon was found in the monzonites. Feldspar and biotite have preferentially lost helium during slow cooling at depth and/or during their transportation to the surface in a hot host magma. The source magma noble gas isotopic features are well preserved in fluid inclusions in hornblende, and indicate that the magma may be directly derived from subcontinental lithospheric mantle metasomatized by an ancient subduction process, or may have formed as a mixture of MORB-like mantle and crustal components. The radiometric ages, geochemical and Nd–Sr isotopic signatures of the Ulleungdo monzonites as well as the presence of mantle-derived helium and argon, suggests that these felsic plutonic rocks evolved from alkali basaltic magma that formed by partial melting of subcontinental lithospheric mantle beneath the back arc basin located along the active continental margin of the southeastern part of the Eurasian plate.     

Numerical estimation of REV and permeability tensor for fractured rock masses by composite element method

The Monte Carlo method is used to generate parent stochastic discrete fracture network, from which a series of fractured rock samples of different sizes and orientations are extracted. The fracture network combined with a regular grid forms composite element mesh of the fractured rock sample, in which each composite element is composed of sub‐elements incised by fracture segments. The composite element method (CEM) for the seepage is implemented to obtain the nodal hydraulic potential as well as the seepage flow rates through the fractured rock samples. The application of CEM enables a large quantity of stochastic tests for the fractured rock samples because the pre‐process is facilitated greatly. By changing the sizes and orientations of the samples, the analysis of the seepage characteristics is realized to evaluate the variation of the permeability components, the existence of the permeability tensor and the representative element volume. The feasibility and effectiveness are illustrated in a numerical example. Copyright © 2008 John Wiley & Sons, Ltd.     

Biaxial test simulations using a packing of polygonal particles

The mechanical response of cohesionless granular materials under monotonic loading is studied by performing molecular dynamic simulations. The diversity of shapes of soil grains is modelled by using randomly generated convex polygons as granular particles. Results of the biaxial test obtained for dense and loose media show that samples achieve the same void ratio at large strains independent of their initial density state. This limit state resembles the so‐called critical state of soil mechanics, except for some stress fluctuations, which remain for large deformations. These fluctuations are studied at the micro‐mechanical level, by following the evolution of the co‐ordination number, force chains and the fraction of the sliding contacts of the sample. Copyright © 2007 John Wiley & Sons, Ltd.     

Flexural isostatic response of the Alps to increased Quaternary erosion recorded by foreland basin remnants,SE France

We test the hypothesis that flexural isostatic compensation of the mass removed by enhanced Quaternary erosion is responsible for uplift of the Western European Alps and their forelands. We use two well‐preserved and well‐dated (1.8 Ma) abandonment surfaces of foreland basin remnants in SE France (the Chambaran and Valensole plateaux) as passive benchmarks for tilting of the foreland. Estimating their initial slope from morphometric scaling relationships, we determine bulk post‐depositional tilting of 0.5–0.8% for these surfaces. The calculated isostatic response of the Alpine lithosphere to erosional unloading, using the method recently proposed by Champagnac et al. [Geology 35 (2007) 195–198] , yields a predicted tilting of 0.3–0.4% in the considered areas, explaining approximately half of the determined post‐depositional tilting. Such long‐term deformation being insensitive to cyclic loading/unloading because of glaciations, we suspect the other half to be related to as yet undetermined long‐wavelength and long‐lived tectonic process(es).     

Glaciomarine sedimentation and palaeo-glacial setting of Maxwell Bay and its tributary embayment, Marian Cove, South Shetland Islands, West Antarctica

High-resolution (3.5 kHz and multi-channel) seismic profiles and piston cores were collected from Maxwell Bay and its tributary embayment, Marian Cove, in the South Shetland Islands, Antarctica, during the Korea Antarctic Research Program (1992/93 and 1995/96) to elucidate the glaciomarine sedimentation processes and recent glacial history of the area. Seismic data from Maxwell Bay reveal a rugged bay margin and flattened basin floor covered with well-stratified hemipelagic muds. On the base-of-slope, acoustically transparent debris flows occur, indicating downslope resedimentation of glaciomarine sediments. Despite the subpolar and ice-proximal settings of Marian Cove, the seafloor is highly rugged with a thin sediment drape, suggesting that much of the area has been recently eroded by glaciers. Sediment cores from the cove penetrated three distinct fining-upward lithofacies: (1) basal till in the lower part of the core, accumulated just seaward of the grounding line of the tidewater glacier; (2) interlaminated sand and mud in the middle part, deposited in ice-proximal zone by a combination of episodic subglacial meltwater inflow and iceberg dumping; and (3) pebbly mud in the upper part, deposited in ice-distal zone by both surface meltwater plume and ice-rafting from the glacier front. A reconstruction of the glacial history of these areas since the late glacial maximum shows an ice sheet filling Maxwell Bay in late Wisconsin time and grounding of the tidewater glacier in Marian Cove until about 1300 yr BP.     

Evaluation of bending moments and shear forces at unit connections of very large floating structures using hydroelastic and rigid body analyses

This paper investigates the characteristics of bending moments, shear forces and stresses at unit connections of very large floating structures (VLFS) under wave loads. The responses of VLFS are calculated by solving multi-body motion equation considering hydroelasticity and connection stiffness. Hydroelastic responses are calculated by the direct method. Higher-order boundary element method (HOBEM) is used for fluid analysis and finite element method (FEM) is introduced for structural analysis. The equation of motion is modified to describe the unit connections by employing spring elements. Bending moments and shear forces at the connections are obtained from the dynamic equilibrium condition for pressures and inertia forces. Two types of VLFS units such as tandem arranged units and side-by-side arranged units are considered in the numerical examples. The influences of connection stiffness, wave frequency and heading angle on responses of VLFS are investigated through the numerical examples. Rigid body analysis along with hydroelastic analysis is also carried out in the numerical analysis and comparison of those two approaches is discussed.