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.     

Theory of quadrupolar sunspots and the active region of August, 1972

Energy is stored when the force-free magnetic field in an active region departs from a potential field, the departure showing up as a shear in the field. As soon as the field untwists, energy will be released to produce flares. Based on this idea, we derived an analytical solution of the equation of force-free field under the assumption of a constant force-free factor, and found expressions for seven important quantities for quadrupolar sunspots: the magnetic energy of the twisted field, that of potential field, the extractable free energy ΔM, the magnetic flux, the total current, the force-free factor and the field decay factor, in terms of three observables: the field intensity, the twist angle and the distance between two spots of the same polarity. The expression for ΔM can be useful in solar prediction work. For the active region of August, 1972, we found ΔM up to 6 × 10³² erg, sufficient to supply the energy of the observed flare activity. Observations of this active region are in good general agreement with our theoretical expectations: in the entire twisting of the quadrupolar sunspot group, each spot assumes the form of a complete spiral in the clockwise direction for each of the four spots.     

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.     

Development and application of a GIS-based tool for earthquake-induced hazard prediction

An integrated GIS-based tool (GTIS) was constructed to estimate site effects related to the earthquake hazards in the Gyeongju area of Korea. To build the GTIS for the study area, intensive site investigations and geotechnical data collections were performed and a walk-over site survey was additionally carried out to acquire surface geo-knowledge data in accordance with the procedure developed to build the GTIS. For practical applications of the GTIS used to estimate the site effects associated with the amplification of ground motion, seismic microzoning maps of the characteristic site period and the mean shear wave velocity to a depth of 30 m were created and presented as a regional synthetic strategy addressing earthquake-induced hazards. Additionally, based on one-dimensional site response analyses, various seismic microzoning maps for short- and mid-period amplification potentials were created for the study area. Case studies of seismic microzonations in the Gyeongju area verified the usefulness of the GTIS for predicting seismic hazards in the region.     

Chemical weathering of carbonates and silicates in the Han River basin,South Korea

A detailed investigation of the fluvial geochemistry of the Han River system allows to estimate the rates of chemical weathering and the consumption of CO₂. The Han River drains approximately 26,000 km² and is the largest river system in South Korea in terms of both water discharge and total river length. It consists of two major tributaries: the North Han River (NHR) and the South Han River (SHR). Distinct differences in basin lithology (silicate vs. carbonate) between the NHR and SHR provide a good natural laboratory in which to examine weathering processes and the influence of basin geology on water quality. The concentrations of major elements and the Sr isotopic compositions were obtained from 58 samples collected in both summer and winter along the Han River system in both 2000 and 2006. The concentrations of dissolved loads differed considerably between the NHR and SHR; compared with the SHR, the NHR had much lower total dissolved solids (TDS), Sr, and major ion concentrations but a higher Si concentration and ⁸⁷Sr/⁸⁶Sr ratio. A forward model showed that the dissolved loads in the NHR came primarily from silicate weathering (55 ± 11%), with a relatively small portion from carbonates (30 ± 14%), whereas the main contribution to the dissolved loads in the SHR was carbonate weathering (82 ± 3%), with only 11 ± 4% from silicates. These results are consistent with the different lithologies of the two drainage basins: silicate rocks in the NHR versus carbonate rocks in the SHR. Sulfuric acid derived from sulfide dissolution in coal-containing sedimentary strata has played an important role in carbonate weathering in the SHR basin, unlike in the NHR basin. The silicate weathering rate (SWR) was similar between the NHR and SHR basins, but the rate of CO₂ consumption in the SHR basin was lower than in the NHR basin due to an important role of sulfuric acid derived from pyrite oxidation.     

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).