Analytical solutions of one‐dimensional infiltration before and after ponding

This paper is a continuation of previous research, which obtained a convenient solution for arbitrary surface fluxes before ponding. By means of Fourier Transformation this has been extended to derive analytical solutions of a linearized Richards' equation for arbitrary input fluxes after surface saturation. Copyright © 2003 John Wiley & Sons, Ltd.     

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.     

Partial Melting and Counterclockwise P T Path of Subducted Oceanic Crust (Sierra del Convento Melange, Cuba)

Partial melting of subducted oceanic crust has been identifiedin the Sierra del Convento mélange (Cuba). This serpentinite-matrixmélange contains blocks of mid-ocean ridge basalt (MORB)-derivedplagioclase-lacking epidote ± garnet amphibolite intimatelyassociated with peraluminous trondhjemitic–tonalitic rocks.Field relations, major element bulk-rock compositions, mineralassemblages, peak metamorphic conditions (c. 750°C, 14–16kbar), experimental evidence, and theoretical phase relationssupport formation of the trondhjemitic–tonalitic rocksby wet melting of subducted amphibolites. Phase relations andmass-balance calculations indicate eutectic- and peritectic-likemelting reactions characterized by large stoichiometric coefficientsof reactant plagioclase and suggest that this phase was completelyconsumed upon melting. The magmatic assemblages of the trondhjemitic–tonaliticmelts, consisting of plagioclase, quartz, epidote, ±paragonite, ± pargasite, and ± kyanite, crystallizedat depth (14–15 kbar). The peraluminous composition ofthe melts is consistent with experimental evidence, explainsthe presence of magmatic paragonite and (relict) kyanite, andplaces important constraints on the interpretation of slab-derivedmagmatic rocks. Calculated P–T conditions indicate counterclockwiseP–T paths during exhumation, when retrograde blueschist-faciesoverprints, composed of combinations of omphacite, glaucophane,actinolite, tremolite, paragonite, lawsonite, albite, (clino)zoisite,chlorite, pumpellyite and phengite, were formed in the amphibolitesand trondhjemites. Partial melting of subducted oceanic crustin eastern Cuba is unique in the Caribbean realm and has importantconsequences for the plate-tectonic interpretation of the region,as it supports a scenario of onset of subduction of a youngoceanic lithosphere during the early Cretaceous (c. 120 Ma).The counterclockwise P–T paths were caused by ensuingexhumation during continued subduction. KEY WORDS: amphibolite; Cuba; exhumation; partial melting; trondhjemite; subduction     

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

Community and Population Responses of the Macroinfauna to Physical Factors over a Range of Exposed Sandy Beaches in South-central Chile

Ten exposed sandy sites covering a range from reflective to dissipative beaches were sampled in south-central Chile to evaluate: (1) spatial changes in species richness, abundance and biomass of the intertidal macroinfauna in response to changes in mean grain size, beach face slope and beach type, and (2) spatial changes in abundance, biomass and body sizes of the most abundant species in response to changes in the physical factors. The number of species, abundance and biomass per beach in general decreased with increasing particle size and beach face slope (steeper beaches) and increased from reflective to dissipative conditions. The best fit for number of species was with Dean's parameter, a measure of beach type, whereas for abundance and biomass the best fits were found with particle size. The isopod Excirolana braziliensis and the anomuran Emerita analoga increased in abundance and biomass towards dissipative conditions, whereas Excirolana hirsuticauda showed the opposite trend in biomass and was significantly larger in beaches with steeper profiles. It is concluded that responses to changes in beach type are more pronounced at community level than within species populations.