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.     

Quiescent prominence threads models

We have calculated prominence thread models for different values of the center temperature and pressure. We have simultaneously solved the radiative transfer, statistical equilibrium and ionization equilibrium equations assuming a three-level atom plus continuum. We have also computed the energy balance equation including the hydrogen radiative losses from our calculations, plus other radiative losses and heat conduction. Some models have been calculated assuming possible variations in thermal conductivity and heating terms. We computed the lines and continua emitted by a number of threads, in order to compare with the observations and evaluate how the different values of the parameters affect the profiles and absolute intensities of L, L, H, and Ly continuum.Member of the Carrera del Investigador, CONICET, Argentina.c.c.67 Suc. 28 Buenos Aires 1428, Argentina.     

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.     

Oxidation of pyrite during extraction of carbonate associated sulfate

The sulfur isotopic composition of carbonate associated sulfate (CAS) has been used to investigate the geochemistry of ancient seawater sulfate. However, few studies have quantified the reliability of δ³⁴S of CAS as a seawater sulfate proxy, especially with respect to later diagenetic overprinting. Pyrite, which typically has depleted δ³⁴S values due to authigenic fractionation associated with bacterial sulfate reduction, is a common constituent of marine sedimentary rocks. The oxidation of pyrite, whether during diagenesis or sample preparation, could thus adversely influence the sulfur isotopic composition of CAS. Here, we report the results of CAS extractions using HCl and acetic acid with samples spiked with varying amounts of pyrite. The results show a very strong linear relationship between the abundance of fine-grained pyrite added to the sample and the resultant abundance and δ³⁴S value of CAS. This data represents the first unequivocal evidence that pyrite is oxidized during the CAS extraction process. Our mixing models indicate that in samples with much less than 1 wt.% pyrite and relatively high δ³⁴S_pyrite values, the isotopic offset imparted by oxidation of pyrite should be much less than ? 4‰. A wealth of literature exists on the oxidation of pyrite by Fe³⁺ and we believe this mechanism drives the oxidation of pyrite during CAS extraction, during which the oxygen used to form sulfate is taken from H₂O, not O₂. Consequently, extracting CAS under anaerobic conditions would only slow, but not halt, the oxidation of pyrite. Future studies of CAS should attempt to quantify pyrite abundance and isotopic composition.     

The passivation of calcite by acid mine water. Column experiments with ferric sulfate and ferric chloride solutions at pH 2

Column experiments, simulating the behavior of passive treatment systems for acid mine drainage, have been performed. Acid solutions (HCl or H₂SO₄, pH 2), with initial concentrations of Fe(III) ranging from 250 to 1500 mg L⁻¹, were injected into column reactors packed with calcite grains at a constant flow rate. The composition of the solutions was monitored during the experiments. At the end of the experiments (passivation of the columns), the composition and structure of the solids were measured. The dissolution of calcite in the columns caused an increase in pH and the release of Ca into the solution, leading to the precipitation of gypsum and Fe–oxyhydroxysulfates (Fe(III)–SO₄–H⁺ solutions) or Fe–oxyhydroxychlorides (Fe(III)–Cl–H⁺ solutions). The columns worked as an efficient barrier for some time, increasing the pH of the circulating solutions from 2 to 6–7 and removing its metal content. However, after some time (several weeks, depending on the conditions), the columns became chemically inert. The results showed that passivation time increased with decreasing anion and metal content of the solutions. Gypsum was the phase responsible for the passivation of calcite in the experiments with Fe(III)–SO₄–H⁺ solutions. Schwertmannite and goethite appeared as the Fe(III) secondary phases in those experiments. Akaganeite was the phase responsible for the passivation of the system in the experiments with Fe(III)–Cl–H⁺ solutions.     

Controls of δS and δO in dissolved sulphate: Learning from a detailed survey in the Llobregat River (Spain)

The S and O isotopic composition of dissolved SO₄, used as a tracer for SO₄ sources, was applied to the water of the Llobregat River system (NE Spain). The survey was carried out at 30 sites where surface water was sampled on a monthly basis over a period of 2a. The concentration of dissolved SO₄ varied from 20 to 1575 mg L⁻¹. Sulphur isotopic compositions clustered in two populations: one – 93% of the samples – had positive values with a mode of +9‰; the other had negative values and a mode of −5‰. Data for δ¹⁸O_SO4 showed a mean value of +11‰, with no bi-modal distribution, though lower values of δ¹⁸O corresponded to samples with negative δ³⁴S. These values can not be explained solely by the contribution of bedrock SO₄ sources: that is, sulphide oxidation and the weathering of outcrops of sulphates, though numerous chemical sediments exist in the basin. Even in a river with a high concentration of natural sources of dissolved SO₄, such as the Llobregat River, the δ³⁴S values suggest that dissolved SO₄ is controlled by a complex mix of both natural and anthropogenic sources. The main anthropogenic sources in this basin are fertilizers, sewage, potash mine effluent and power plant emissions. Detailed river water sampling, together with the chemical and isotopic characterisation of the main anthropogenic inputs, allowed determination of the influence of redox processes, as well as identification of the contribution of natural and anthropogenic SO₄ sources and detection of spatial variations and seasonal changes among these sources. For instance, in the Llobregat River the input of fertilisers is well marked seasonally. Minimum values of δ³⁴S are reported during fertilization periods – from January to March – indicating a higher contribution of this source. The dual isotope approach, δ³⁴S and δ¹⁸O, is useful to better constrain the sources of SO₄. Moreover, in small-scale studies, where the inputs are well known and limited, the mixing models can be enhanced and the contribution of the different sources can be quantified to some extent.     

Quantitative sinkhole hazard assessment. A case study from the Ebro Valley evaporite alluvial karst (NE Spain)

Quantitative sinkhole hazard assessments in karst areas allow calculation of the potential sinkhole risk and the performance of cost-benefit analyses. These estimations are of practical interest for planning, engineering, and insurance purposes. The sinkhole hazard assessments should include two components: the probability of occurrence of sinkholes (sinkholes/km² year) and the severity of the sinkholes, which mainly refers to the subsidence mechanisms (progressive passive bending or catastrophic collapse) and the size of the sinkholes at the time of formation; a critical engineering design parameter. This requires the compilation of an exhaustive database on recent sinkholes, including information on the: (1) location, (2) chronology (precise date or age range), (3) size, and (4) subsidence mechanisms and rate. This work presents a hazard assessment from an alluvial evaporite karst area (0.81 km²) located in the periphery of the city of Zaragoza (Ebro River valley, NE Spain). Five sinkholes and four locations with features attributable to karstic subsidence where identified in an initial investigation phase providing a preliminary probability of occurrence of 0.14 sinkholes/km² year (11.34% in annual probability). A trenching program conducted in a subsequent investigation phase allowed us to rule out the four probable sinkholes, reducing the probability of occurrence to 0.079 sinkholes/km² year (6.4% in annual probability). The information on the severity indicates that collapse sinkholes 10–15 m in diameter may occur in the area. A detailed study of the deposits and deformational structures exposed by trenching in one of the sinkholes allowed us to infer a modern collapse sinkhole approximately 12 m in diameter and with a vertical throw of 8 m. This collapse structure is superimposed on a subsidence sinkhole around 80 m across that records at least 1.7 m of synsedimentary subsidence. Trenching, in combination with dating techniques, is proposed as a useful methodology to elucidate the origin of depressions with uncertain diagnosis and to gather practical information with predictive utility about particular sinkholes in alluvial karst settings: precise location, subsidence mechanisms and magnitude, and timing and rate of the subsidence episodes.