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.     

HIGH POSSIL AND STRATHMORE—A STUDY OF TWO L6 CHONDRITES

High Possil and Strathmore are the first (1804) and last (1917) meteorite falls, respectively, of the three recorded in Scotland. Olivine compositions and total Fe contents in High Possil (Fa_25.2; 21.35 wt %) and Strathmore (Fa_25.3; 20.6 wt %) confirm their classification as L-group chondrites (Mason, 1963), and the presence of abundant plagioclase feldspar shows that both chondrites belong to petrologic type 6. Both chondrites display thermal and mechanical alteration attributable to moderate shock-loading appropriate to facies c (High Possil) and c-d (Strathmore) (Dodd and Jarosewich, 1979). Incipient shock-melting of metal and troilite in both chondrites is the first described from Lc chondrites, and differences in the responses of metallic and silicate minerals to shock-loading are discussed.     

The effects of sedimentation and compaction on oceanic heat flow

Summary. The estimation of environmental effects forms an important part of the interpretation of oceanic heat flow measurements. In particular, the perturbations associated with sedimentation and surface temperature changes must be taken into account. Analytical solutions can be obtained only for individual, simplified versions of these problems, whereas any real example is complicated by the process of sediment compaction which changes the bulk thermal properties with depth. A physical model is developed which uses sediment porosity trends to predict the thermal parameters and material advection rates for an evolving sediment/basement system. These values are then used in a numerical solution to the heat flow equation to give estimates of the perturbed surface heat flux through time. In addition to variations in sedimentation rate, sediment type, radioactive heat production and surface temperature changes are considered. Heat flow corrections may vary by up to a factor of 2 according to sediment type while radioactive heat production can offset the effects of sedimentation by as much as 40 per cent. The results also indicate that alterations determined from simple analytical models tend to over-estimate the true perturbation to the flux.     

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.     

Numerical modeling for analyzing thermal surface anomalies induced by underground coal fires

Coal seams burning underneath the surface are recognized all over the world and have drawn increasing public attention in the past years. Frequently, such fires are analyzed by detecting anomalies like increased exhaust gas concentrations and soil temperatures at the surface. A proper analysis presumes the understanding of involved processes, which determine the spatial distribution and dynamic behavior of the anomalies.In this paper, we explain the relevance of mechanical and energy transport processes with respect to the occurrence of temperature anomalies at the surface. Two approaches are presented, aiming to obtain insight into the underground coal fire situation: In-situ temperature mapping and numerical simulation. In 2000 to 2005, annual temperature mapping in the Wuda (Inner Mongolia, PR China) coal fire area showed that most thermal anomalies on the surface are closely related to fractures, where hot exhaust gases from the coal fire are released. Those fractures develop due to rock mechanical failure after volume reduction in the seams. The measured signals at the surface are therefore strongly affected by mechanical processes.More insight into causes and effects of involved energy transport processes is obtained by numerical simulation of the dynamic behavior of coal fires. Simulations show the inter-relation between release and transport of thermal energy in and around underground coal fires. Our simulation results show a time delay between the coal fire propagation and the observed appearance of the surface temperature signal. Additionally, the overall energy flux away from the burning coal seam into the surrounding bedrock is about 30-times higher than the flux through the surface. This is of particular importance for an estimation of the energy released based on surface temperature measurements. Finally, the simulation results also prove that a fire propagation rate estimated from the interpretation of surface anomalies can differ from the actual rate in the seam.     

Experimental solidification of anhydrous latitic and trachytic melts at different cooling rates: The role of nucleation kinetics

Two sets of cooling experiments were run at atmospheric conditions for two anhydrous starting latitic and trachytic melts: 1) five cooling rates (25, 12.5, 3, 0.5, and 0.125 °C/min) between 1300° and 800 °C, and 2) a 0.5 °C/min cooling rate from 1300 °C with quench temperatures at 1200°, 1100°, 1000° and 900 °C. Trachytic run-products are invariably glassy. Nucleation is also suppressed in the latitic run-products at the three highest cooling rates. Conversely, in the 0.5 and 0.125 °C/min runs, latites have a crystal content of  90 vol.%. The phases are: plagioclase, clinopyroxene, glass and iron-bearing oxide (in order of abundance). The variable quench temperatures, investigated by coupling experiments with Pt wire and Pt capsule sample containers in set 2, again did not produce crystallization of trachyte, whereas latitic samples are characterized by  10 vol.% of oxides, pyroxenes and plagioclase (in order of appearance), at temperature < 1000 °C. Effects of (preferential) heterogeneous nucleation on sample holders, of superheating degree, and chemical species loss during cooling are absent for both melt compositions. The difference of solidification paths between these two silicate melts can be ascribed only to their small chemical differences. In comparison with calculated equilibrium conditions all the experimental latitic and trachytic run-products revealed strong kinetic effects, interpretable in the light of the nucleation theory. The glass-forming ability (GFA) of trachyte is higher, whereas their critical cooling rate (Rc) is lower (< 0.125 °C/min), in comparison to latitic melts (Rc > 0.5 °C/min). The experimental results carried out in this study can be applied to lava flows and domes; trachytic lavas are able to flow for longer period with respect to latitic ones in a metastable condition. Glass-rich terrestrial lavas, i.e. obsidians, can be the result of sluggish nucleation kinetics due to the relative high polymerisation of evolved silicate melts.     

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.