Finite element formulation and algorithms for unsaturated soils. Part II: Verification and application

The finite‐element formulation and integration algorithms developed in Part I are used to analyse a number of practical problems involving unsaturated and saturated soils. The formulation and algorithms perform well for all the cases analysed, with the robustness of the latter being largely insensitive to user‐defined parameters such as the number of coarse time steps and error control tolerances. The efficiency of the algorithms, as measured by the CPU time consumed, does not depend on the number of coarse time steps, but may be influenced by the error control tolerances. Based on the analyses presented here, typical values for the error control tolerances are suggested. It is also shown that the constitutive modelling framework presented in Part I can, by adjusting one constitutive equation and one or two material parameters, be used to simulate soils that expand or collapse upon wetting. Treating the suction as a strain variable instead of a stress variable proves to be an efficient and robust way of solving suction‐dependent plastic yielding. Moreover, the concept of the constitutive stress is a particularly convenient way of handling the transition between saturation and unsaturation. Copyright © 2003 John Wiley & Sons, Ltd.     

The effects of Population III stars on the chemical and photometrical evolution of ellipticals

We have studied the effects of a hypothetical initial generation containing very massive stars [   M > 100 M_⊙  , pair-creation supernovae] on the chemical and photometric evolution of elliptical galaxies. To this purpose, we have computed the evolution of a typical elliptical galaxy with luminous mass of the order of  10¹¹ M_⊙  and adopted chemical evolution models already tested to reproduce the main features of ellipticals. We have tested several sets of yields for very massive zero-metallicity stars: these stars should produce quite different amounts of heavy elements than lower-mass stars. We found that the effects of Population III stars on the chemical enrichment is negligible if only one or two generations of such stars occurred, whereas they produce quite different results from the standard models if they continuously formed for a period not shorter than 0.1 Gyr. In this case, the results are at variance with the main observational constraints of ellipticals such as the average  [〈α/ Fe〉_*]  ratio in stars and the integrated colours. Therefore, we conclude that if Population III stars ever existed they must have been present for a very short period of time and their effects on the following evolution of the parent galaxy must have been negligible. This effect is minimum if a more realistic model with initial infall of gas rather than the classic monolithic model is adopted. Ultimately, we conclude that there is no need to invoke a generation of very massive stars in ellipticals to explain their chemical and photometric properties.     

Low-cost volcano surveillance from space: case studies from Etna, Krafla, Cerro Negro, Fogo, Lascar and Erebus

 Satellite data offer a means of supplementing ground-based monitoring during volcanic eruptions, especially at times or locations where ground-based monitoring is difficult. Being directly and freely available several times a day, data from the advanced very high resolution radiometer (AVHRR) offers great potential for near real-time monitoring of all volcanoes across large (3000×3000 km) areas. Herein we describe techniques to detect and locate activity; estimate lava area, thermal flux, effusion rates and cumulative volume; and distinguish types of activity. Application is demonstrated using data for active lavas at Krafla, Etna, Fogo, Cerro Negro and Erebus; a pyroclastic flow at Lascar; and open vent systems at Etna and Stromboli. Automated near real-time analysis of AVHRR data could be achieved at existing, or cheap to install, receiving stations, offering a supplement to conventional monitoring methods. Received: 21 January 1997 / Accepted: 3 April 1997     

Chemical evolution of Elliptical Galaxies and the ICM

We present a new model for the chemical evolution of elliptical galaxies taking into account SN feedback, detailed nucleosynthesis and galactic winds. We discuss the effect of galactic winds on the chemical enrichment of the ICM and compute the energy per particle injected by the galaxies into the ICM. This revised version was published online in August 2006 with corrections to the Cover Date.     

On the restricted three-body problem with generalized forces

By introducing general functions which depend on distance, a general scheme which determines the equilibrium solutions for the generalized restricted three-body problem is given. Applications to problems such as primaries considered as rigid bodies, influence of the radiation pressure of the primaries, and a combination of radiation pressure and rigid body are presented.     

Trace elements contamination of agricultural soils affected by sulphide exploitation (Iberian Pyrite Belt,Sw Spain)

Agricultural soils of the Riotinto mining area (Iberian Pyrite Belt) have been studied to assess the degree of pollution by trace elements as a consequence of the extraction and treatment of sulphides. Fifteen soil samples were collected and analysed by ICP-OES and INAA for 51 elements. Chemical analyses showed an As–Cu–Pb–Zn association related with the mineralisation of the Iberian Pyrite Belt. Concentrations were 19–994 mg kg⁻¹ for As, 41–4,890 mg kg⁻¹ for Pb, 95–897 mg kg⁻¹ for Zn and of 27–1,160 mg kg⁻¹ for Cu. Most of the samples displayed concentrations of these elements higher than the 90th percentile of the corresponding geological dominium, which suggests an anthropogenic input besides the bedrock influence. Samples collected from sediments were more contaminated than leptosols because they were polluted by leachates or by mining spills coming from the waste rock piles. The weathering of the bedrock is responsible for high concentrations in Co, Cr and Ni, but an anthropogenic input, such as wind-blown dust, seems to be indicative of the high content of As, Cu, Pb and Zn in leptosols. The metal partitioning patterns show that most trace elements are associated with Fe amorphous oxy-hydroxides, or take part of the residual fraction. According to the results obtained, the following mobility sequence is proposed for major and minor elements: Mn, Pb, Cd, > Zn, Cu > Ni > As > Fe > Cr. The high mobility of Pb, Cu and Zn involve an environmental risk in this area, even in soils where the concentrations are not so high.     

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 use of principle component analyses in characterising trace and major elemental distribution in a 55 kyr peat deposit in tropical Australia: Implications to paleoclimate

The Lynch’s Crater peat deposit in NE-Australia is a sensitive environmental archive located in the tropical Southern Hemisphere. This unique deposit illustrates that local and regional changes had a profound effect on the local Australian ecosystem over the past 55 kyr. To obtain a proxy of past climate changes, trace and major element geochemistry analyses were applied to a 13 m peat core from the crater. Principle component analysis (PCA) was used to identify the main factors that control elemental distribution in the peat and to add interpretative strength to the geochemical behavior of selected major and trace elements. For example, Sc, Al, Cu, and Pb were found to be related to increased erosion of the basin soils, and from this, several periods of significant flux from atmospheric input and/or terrigenous run-off were identified. Geochemically mobile elements during rock weathering and pedogenesis, such as Mg, Ca, and Sr helped to identify the peat ombrotrophic-minerotrophic boundary at ∼1.5 m depth and offered important information about fluxes of these nutrients to the mire and their dynamics within the deposit. Arsenic and V comparisons between the peat record (high concentrations in some peat sections) and in local basin rocks (very low concentrations), suggested the presence of a long range, atmospheric dust source early in the formation of the mire. The Lynch’s Crater peat record presents a continuous record of environmental change in tropical Australia and contributes new understanding to geochemical processes in peatlands.     

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.