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.     

An observational study of coronal hole rotation over the sunspot cycle

An analysis of the rotation of coronal holes (CHs) spanning 18 years was done based on data from theCatalogue of Coronal Holes (Sanchez-Ibarra and Barraza-Paredes, 1992). A differential rotation of CHs is confirmed for the totality of CHs, but a different behavior was found when those were separated as equatorial or isolated, and polar hole extensions, such as in theCatalogue. Isolated CHs show a typical differential rotation, but polar hole extensions display two different types of behavior: a rotation rate below 40° ± 5° of heliographic latitude, increasing to the equator, and a rotation rate above the same heliographic latitude but increasing to the poles. Also discussed here is how this last result agrees with other studies that indicate the mostly rigid rotation of the corona at higher latitudes.     

Study of geothermal water intrusion due to groundwater exploitation in the Puebla Valley aquifer system, Mexico

Significant intrusion of geothermal water into fresh groundwater takes place in the Puebla Valley aquifer system, Mexico. The decline in the potentiometric surface due to the overexploitation of the groundwater induces this intrusion. This hydrological system comprises three aquifers located in Plio-Quaternary volcanic sediments and Mesozoic calcareous rocks. The hydraulic balance of the aquifer shows that the annual output exceeds the natural inputs by 12 million m³. Between 1973 and 2002, a drop in the potentiometric surface, with an 80 m cone of depression, was identified in a 5-km-wide area located southwest of the city of Puebla. Chemical analyses performed on water samples since 1990 have shown an increase in total dissolved solids (TDS) of more than 500 mg/L, coinciding with the region showing a cone of depression in the potentiometric surface. A three-dimensional flow and transport model, based on the hydrogeological and geophysical studies, was computed by using the MODFLOW and MT3D software. This model reproduces the evolution of the aquifer system during the last 30 years and predicts for 2010 an additional drawdown in the potentiometric surface of 15 m, and an increase in the geothermal water intrusion.     

Sub-wavebase cross-bedded grainstones on a distally steepened carbonate ramp, Upper Miocene, Menorca, Spain

Cross‐bedded grainstones on carbonate ramps and shelves are commonly related to the locus of major wave energy absorption such as shorelines, shoals or shelf breaks. In contrast, on the Early Tortonian carbonate platform of Menorca (Balearic Islands), coarse‐grained, cross‐bedded grainstones are found at a distance from the palaeoshoreline where they were deposited below the wavebase. Excellent exposures along continuous outcrops on the sea cliffs of Menorca reveal the depositional profile and three‐dimensional distribution of the different facies belts of the Tortonian ramp depositional system. Basinward from the palaeoshoreline, fan deltas and beach deposits pass into 5‐km‐wide gently dipping bioturbated dolopackstone (inner and middle ramp), then into 12–20°‐dipping dolograinstone/rudstone clinobeds (ramp slope) and, finally, into subhorizontal fine‐grained basinal dolowackestone to dolopackstone (outer ramp). In this Miocene example, coarse‐grained grainstones exist in five different settings other than beach deposits: (1) on the middle ramp, where cross‐bedded grainstones were deposited by currents roughly parallel to the shoreline at 40–70 m estimated water depth and are interbedded with gently dipping bioturbated dolomitized packstones; (2) on the upper slope, where clinobeds are composed mostly of in situ rhodoliths and red‐algae fragments; (3) on the lower slope, as small‐scale bedforms (small three‐dimensional subaqueous dunes) migrating parallel to the slope; (4) at the transition between the lower slope and the outer ramp, where mollusc‐rich and rhodolithic rudstones and grainstones, interbedded in dolomitized laminated wackestones containing abundant planktonic foraminifera, infill slide/slump scars as upslope‐backstepping bodies (backsets); (5) at the toe of the slope, where coarse skeletal grainstones indicate bedform migration parallel to the platform margin, induced by currents at more than 150 m estimated water depth. This Late Miocene example also illustrates how changes in intrabasinal environmental conditions (nutrients and/or temperature) may produce changes in stratal patterns and facies architecture if they affect the biological system. Two depositional sequences compose the Miocene platform on Menorca, where a reef‐rimmed platform prograded onto an earlier distally steepened ramp. The transition from the ramp to the reef‐rimmed platform was effected by an increase in accommodation space caused by ecological changes, promoting a shift from a grain‐ to a framework‐producing biota.     

Ray tracing in elliptical anisotropic media using the linear traveltime interpolation (LTI) method applied to traveltime seismic tomography

The linear traveltime interpolation (LTI) method is a suitable ray‐tracing technique for modelling first‐arrival times in isotropic media. LTI is extended to elliptical anisotropic media and applied to a tomographic inversion procedure. A theoretical formulation is first derived and then LTI implementation is discussed in terms of source–receiver arrays and cell size. The method is then combined with the tomographic inversion procedure adopted. The matching of the ray tracing with inversion in elliptical anisotropic media posed a double non‐linear problem. Thus two assumptions were made: the velocity in each cell is uniform and the main directions of anisotropy are known. To take into account the geometrical characteristics of the area under investigation (depth and velocity of the weathering, and thickness of the inner media), cells of varying size were considered. No hypothesis was made on anisotropy weakness.
The algorithm was first tested on synthetic models and then applied to a field survey. On comparing the results of the synthetic models and the field survey with those obtained with a linear raypath approximation, it was found that there were fewer data misfits.     

δO of carbonate, quartz and phosphate from belemnite guards: implications for the isotopic record of old fossils and the isotopic composition of ancient seawater

Belemnite guards of Cretaceous and Jurassic age were found to contain varying amounts of quartz deposited both on the external surface and inside the rostrum. The oxygen isotopic composition of coexisting carbonate, quartz and phosphate from the same rostrum was measured according to well-established techniques. None of these compounds showed isotopic values in equilibrium with one another. Assuming δ¹⁸O values of the diagenetic water within the range of meteoric waters, the δ¹⁸O(SiO₂) yield temperatures in agreement with the apparent secondary origin of this phase. The δ¹⁸O(CO₃²⁻) range, with a certain continuity, between −10.8 and +0.97 PDB-1 with most of the intermediate values being within the range of the carbonate isotopic values of Mesozoic fossils. The most positive isotopic results obtained from phosphate are close to +23/+24‰ (V-SMOW). They can hardly be related to a secondary origin of the phosphate, or to the presence of diagenetic effects, since these results are among the most positive ever measured on phosphate. As far as we know there is no widespread diagenetic process determining an ¹⁸O enrichment of phosphate. The very low concentration of phosphate did not allow the determination of its mineralogical composition. All the available δ¹⁸O(PO₄³⁻) values from belemnite and non-belemnite fossils of marine origin of Tertiary and Mesozoic age are reported along with the newly measured belemnites. The following conclusions may be drawn from the data reported: (1) the pristine oxygen isotope composition of fossil marine organisms (either carbonate or phosphate) may easily undergo fairly large changes because of oxygen isotope exchange processes with diagenetic water; this process is apparent even in the case of geologically recent fossils; (2) the δ¹⁸O(PO₄³⁻) of belemnite rostra seems to be, at least in the case of the most positive results, in isotopic equilibrium with environmental water because of the similarity between the results from Cretaceous belemnites and the results from Cretaceous and Lower Tertiary pelecypods and fish teeth; 3) if so, the only feasible interpretation that can be suggested for the ¹⁸O enriched data is the possibility of a relatively large variation of the oxygen isotopic composition of ocean paleowater from Jurassic to recent time.     

Validation of Landslide Susceptibility Maps; Examples and Applications from a Case Study in Northern Spain

A procedure for validating landslide susceptibility maps wasapplied in a study area in northern Spain and the results obtained compared. Validationwas used to carry out sensitivity analysis for individual variables and combinationsof variables. The validity of different map-making methods was tested, as well as theutility of different types of Favourability Functions. The results obtained show thatvalidation is essential to determine the predictive value of susceptibility maps. Italso helps to better select the most suitable function and significant variables, thus improving the efficiency of the mapping process. Validation based on a temporal strategy makes it possible to derive hazard maps from susceptibility maps.     

Trends in the frequency of synoptic types in central-southern Chile in the period 1961–2012 using the Jenkinson and Collison synoptic classification

Atmospheric circulation patterns in southern Chile (42° 30′ S) were studied in order to determine and analyse the most characteristic synoptic types and their recent trends, as well as to gain an understanding of how they are associated with low-frequency variability patterns. According to the Jenkinson and Collison (J&C) classification method, a 16-point grid of sea-level pressure data was employed. The findings reveal that some synoptic types show statistically significant trends with a 95% confidence level, positively for anticyclonic westerly hybrids (AW) and advective types for third and fourth quadrant wind flows (W, NW, and N) and negatively for SW and cyclonic hybrids (CS and CSW). A model has been constructed of the linear regression of some weather types with teleconnections that most affect Chile: the undetermined types (U), AW were associated with El Niño or the warm phase of the Pacific Decadal Oscillation (PDO), whereas the cyclonic northerly and cyclonic northeasterly types (CN and CNE) were associated with La Niña or cool phase of the PDO. The weather types associated with Antarctic Oscillation (AAO) in its positive phase are anticyclonic northerly and northeasterly and northerly advection types, while in its negative phase are cyclonic southwesterly and advection types.     

The Removal of Dissolved Metals by Hydroxysulphate Precipitates during Oxidation and Neutralization of Acid Mine Waters, Iberian Pyrite Belt

This study examines the removal of dissolved metals during the oxidation and neutralization of five acid mine drainage (AMD) waters from La Zarza, Lomero, Esperanza, Corta Atalaya and Poderosa mines (Iberian Pyrite Belt, Huelva, Spain). These waters were selected to cover the spectrum of pH (2.2–3.5) and chemical composition (e.g., 319–2,103 mg/L Fe; 2.85–33.3 g/L SO₄⁼) of the IPB mine waters. The experiments were conducted in the laboratory to simulate the geochemical evolution previously recognized in the field. This evolution includes two stages: (1) oxidation of dissolved Fe(II) followed by hydrolysis and precipitation of Fe(III), and (2) progressive pH increase during mixing with fresh waters. Fe(III) precipitates at pH < 3.5 (stages 1 and 2) in the form of schwertmannite, whereas Al precipitates during stage 2 at pH 5.0 in the form of several hydroxysulphates of variable composition (hydrobasaluminite, basaluminite, aluminite). During these stages, trace elements are totally or partially sorbed and/or coprecipitated at different rates depending basically on pH, as well as on the activity of the SO₄⁼ anion (which determines the speciation of metals). The general trend for the metals which are chiefly present as aqueous free cations (Pb²⁺, Zn²⁺, Cu²⁺, Cd²⁺, Mn²⁺, Co²⁺, Ni²⁺) is a progressive sorption at increasing pH. On the other hand, As and V (mainly present as anionic species) are completely scavenged during the oxidation stage at pH < 3.5. In waters with high activities (> 10⁻¹) of the SO ₄⁼ ion, some elements like Al, Zn, Cd, Pb and U can also form anionic bisulphate complexes and be significantly sorbed at pH < 5. The removal rates at pH 7.0 range from around 100% for As, V, Cu and U, and 60–80% for Pb, to less than 20% for Zn, Co, Ni and Mn. These processes of metal removal represent a significant mechanism of natural attenuation in the IPB.