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.     

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

Influence of ENSO and of the Indian Ocean Dipole on the Indian summer monsoon variability

Indian summer monsoon (ISM) variability is forced from external factors (like the El Niño Southern Oscillation, ENSO) but it contains also an internal component that tends to reduce its potential for predictability. Large-scale and local monsoon indices based on precipitation and atmospheric circulation parameters are used as a measure of ISM variability. In a 9-members ensemble of AMIP-type experiments (with same boundary SST forcing and different initial conditions) their potential predictability is comparable using both local and large-scale monsoon indices. In the sample analyzed, about half of more predictable monsoon years coincide with El Niño and/or positive Indian Ocean Dipole (IOD) events. Summer monsoon characteristics during ENSO and IOD years are analyzed through composites computed over a three years period (i.e. one year before and one year after the event peak) to investigate the mutual relationship between the events lagged in time. The connection between ISM and IOD is mostly confined in the summer and autumn, while that with ENSO is stronger and extends more in time. In the coupled model results the IOD influence on the monsoon is large, even because in the model IOD events are intense and easily reproduced due to a strong air-sea feedback in the eastern side of the basin. Monsoon seasons preceding or following an El Niño or a La Niña event are not exactly symmetric, even in terms of their biennial character. In most of the cases, both in reanalysis and model, El Niño and positive IOD events tend to co-occur with larger anomalies either in the Indo-Pacific ocean sector or over India, while La Niña and negative IOD do not. From the observed record, the ENSO-IOD correlation is positive strong and significant since mid-60s and it may correspond with either strong or weak ENSO-monsoon relationship and with strong or weak IOD-monsoon relationship. A main difference between those periods is the relationship between Indian monsoon rainfall and SST in other ocean basins rather than the Indo-Pacific sector alone.     

Coralligenous habitat in the northern Adriatic Sea: an overview

The rocky outcrops irregularly scattered in the sandy‐muddy sea bed of the Northern Adriatic, formed by a base rock and a superimposed concretion of vegetal and animal organisms, are considered an infralittoral coralligenous habitat, according to the most recent definition. In the last four decades, research has mainly concentrated on their geology because of the question of the origin of the base rock – very similar to beachrock – which has only recently been attributed to sedimentary cementation induced by methane seeps. Studies on their macrobenthic assemblages have also been published recently, but true comparisons among them are difficult because they mostly refer only to a few selected phyla or to short‐term observations. Very few papers deal with the whole of the animal and plant communities, including all taxonomic groups and the most important environmental variables. On the basis of the literature data concerning several outcrops in the Gulf of Venice, we suggest that the Adriatic reefs differ slightly from the classic deeper coralligenous assemblages around the Mediterranean coast because they have smaller concretions of algae builders. We also suggest that the Northern Adriatic coralligenous reefs are younger than the other reefs occurring in the Mediterranean Sea.     

Palaeoshoreline for the Late Cretaceous marine platform in the Iberian Trough (Leonese Area, Spain) deduced from outcrop and subsurface analysis

The location of the Late Cretaceous paleoshoreline in the Leonese Area (Iberian Trough, Spain) has been investigated by seismic analysis through isobath and isopach maps. The succession can be divided into two depositional sequences: DS-1 and DS-2. These sequences are composed of fluvial systems at the base, with paleocurrents that flowed eastward and north-eastward. The DS-1 sequence (Late Albian-Middle Turonian) shows intertidal to subtidal and offshore deposits at the top, while the DS-2 sequence (Late Turonian-Campanian) presents intertidal to subtidal, tidal flat and shallow marine and lacustrine deposits at its top. The stratigraphic cyclicity based on systems tracts shows that these two depositional sequences exhibit remarkable eustatic control. Both sequences start at the base with a significant sedimentary supply from fluvial systems, related to eustatic fall episodes, and conclude at the top with transgressive periods. The evolution of the basin reveals the history of base-level changes and associated shifts in depositional trends during successive stages. The deeper sectors of the DS-1 sequence are located towards the northeastern part of the study area while the proximal portion of the basin-margin is located to the southwest. The paleoshoreline is placed in a direction oriented at 120. The variations in thickness are elongated in orientations between 030 and 050 and are mainly related to paleovalleys and tributary fluvial networks that supply sediment through the shoreline. It is possible that these variations may be related to active synsedimentary faults. Depocenters move toward the northeast and east during the DS-2 sequence while the proximal portion of the basin-margin moves to the southwest. The paleoshoreline has an orientation of 155 direction and moves basinward. The isopach maps show a group of corridors oriented at 130?C140 and 165?C170, interpreted as result of accumulation of sandy bodies such as inter and subtidal bars. The fluvial systems are transversal to the paleoshoreline direction.     

Is Prediction of Future Landslides Possible with a GIS?

This contribution explores a strategy for landslide hazard zonation inwhich layers of spatial data are used to represent typical settings inwhich given dynamic types of landslides are likely to occur. Theconcepts of assessment and prediction are defined to focus on therepresentation of future hazardous events and in particular on themyths that often provide obstacles in the application of quantitativemethods. The prediction rate curves for different applications describethe support provided by the different data layers in experiments inwhich the typical setting of hazardous events is approximated bystatistically integrating the spatial information.     

Gravimetric and magnetic fabric study of the Sintra Igneous complex: laccolith-plug emplacement in the Western Iberian passive margin

International Journal of Earth Sciences - The geometry and emplacement of the ~ 96 km2, Late Cretaceous Sintra Igneous complex (SIC, ca. 80 Ma) into the West Iberian passive...     

Local- and global-scale seismic analyses of historical masonry compounds in San Pio delle Camere (L’Aquila,Italy)

Masonry building aggregates are large parts of the Italian building heritage often designed without respecting seismic criteria. The current seismic Italian code does not foresee a clear calculation method to predict their static nonlinear behaviour. For this reason, in this paper firstly, a simple methodology to forecast the seismic response of masonry aggregates in San Pio delle Camere (L’Aquila, Italy) has been set up starting from the provisions of the Italian Guidelines on Cultural Heritage. The implemented procedure has been calibrated on the results of two FEM structural analysis programs used to investigate three masonry building compounds. As a result, a design chart used to correctly predict the base shear of aggregate masonry units starting from code provisions has been set up. Later on, the large-scale seismic vulnerability and damage appraisal of the inspected historical centre has been done on the basis of a quick methodology, already implemented and experienced by the author in some historical centres of the Campania region. The analysis result was a numerical correlation between vulnerability index and mean damage grade of examined building compounds. In particular, a damage forecast under numerical way has been firstly estimated and then compared with the real one. The post-earthquake scenario has represented an ideal term of comparison for effectively testing the reliability of the employed technique, which should be further extended to other Italian historical centres.

     

Crystallization conditions of Los Angeles, a basaltic Martian meteorite

The texture of Los Angeles (stone 1) is dominated by relatively large (0.5−2.0 mm) anhedral to subhedral grains of pyroxene, and generally subhedral to euhedral shocked plagioclase feldspar (maskelynite). Minor phases include subhedral titanomagnetite and ilmenite, Fe-rich olivine, olivine+augite-dominated symplectites [some of which include a Si-rich phase and some which do not], pyrrhotite, phosphate(s), and an impact shock-related alkali- and silica-rich glass closely associated with anhedral to euhedral silica grains. Observations and model calculations indicate that the initial crystallization of Mg-rich pigeonitic pyroxenes at ≤1150 °C, probably concomitantly with plagioclase, was followed by pigeonitic and augitic compositions between 1100 and 1050 °C whereas between 1050 and 920 to 905 °C pyroxene of single composition crystallized. Below 920 to 905 °C, single composition Fe-rich clinopyroxene exsolved to augite and pigeonite. Initial appearance of titanomagnetite probably occurred near 990 °C and FMQ-1.5 whereas at and below 990 °C and ≥FMQ-1.5 titanomagnetite and single composition Fe-rich clinopyroxene may have started to react, producing ilmenite and olivine. However, judging from the most common titanomagnetite compositions, we infer that most of this reaction likely occurred between 950 and 900 °C at FMQ-1.0±0.2 and nearly simultaneously with pyroxene exsolution, thus producing assemblages of pigeonite, titanomagnetite, olivine, ilmenite, and augite. We deem this reaction as the most plausible explanation for the formation of the olivine+augite-dominated symplectites in Los Angeles. But we cannot preclude possible contributions to the symplectites from the shock-related alkali- and silica-rich glass or shocked plagioclase, and the breakdown of Fe-rich pigeonite compositions to olivine+augite+silica below 900 °C. Reactions between Fe-Ti oxides and silicate minerals in Los Angeles and other similar basaltic Martian meteorites can control the T-fO₂ equilibration path during cooling, which may better explain the relative differences in fO₂ among the basaltic Martian meteorites.