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.     

Oxygen Isotope Evidence for Chemical Interaction of Ki lauea Historical Magmas with Basement Rocks

Klauea historical summit lavas have a wide range in matrix ¹⁸O_VSMOWvalues (4·9–5·6) with lower values in rockserupted following a major summit collapse or eruptive hiatus.In contrast, ¹⁸O values for olivines in most of these lavasare nearly constant (5·1 ± 0·1). The disequilibriumbetween matrix and olivine ¹⁸O values in many samples indicatesthat the lower matrix values were acquired by the magma afterolivine growth, probably just before or during eruption. BothMauna Loa and Klauea basement rocks are the likely sources ofthe contamination, based on O, Pb and Sr isotope data. However,the extent of crustal contamination of Klauea historical magmasis probably minor (< 12%, depending on the assumed contaminant)and it is superimposed on a longer-term, cyclic geochemicalvariation that reflects source heterogeneity. Klauea's heterogeneoussource, which is well represented by the historical summit lavas,probably has magma ¹⁸O values within the normal mid-ocean ridgebasalt mantle range (5·4–5·8) based on thenew olivine ¹⁸O values. KEY WORDS: Hawaii; Klauea; basalt; oxygen isotopes; crustal contamination     

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.     

The electrical structure of the Slave craton

The Slave craton in northwestern Canada, a relatively small Archean craton (600×400 km), is ideal as a natural laboratory for investigating the formation and evolution of Mesoarchean and Neoarchean sub-continental lithospheric mantle (SCLM). Excellent outcrop and the discovery of economic diamondiferous kimberlite pipes in the centre of the craton during the early 1990s have led to an unparalleled amount of geoscientific information becoming available.

Over the last 5 years deep-probing electromagnetic surveys were conducted on the Slave, using the natural-source magnetotelluric (MT) technique, as part of a variety of programs to study the craton and determine its regional-scale electrical structure. Two of the four types of surveys involved novel MT data acquisition; one through frozen lakes along ice roads during winter, and the second using ocean-bottom MT instrumentation deployed from float planes.

The primary initial objective of the MT surveys was to determine the geometry of the topography of the lithosphere–asthenosphere boundary (LAB) across the Slave craton. However, the MT responses revealed, completely serendipitously, a remarkable anomaly in electrical conductivity in the SCLM of the central Slave craton. This Central Slave Mantle Conductor (CSMC) anomaly is modelled as a localized region of low resistivity (10–15 Ω m) beginning at depths of 80–120 km and striking NE–SW. Where precisely located, it is spatially coincident with the Eocene-aged kimberlite field in the central part of the craton (the so-called “Corridor of Hope”), and also with a geochemically defined ultra-depleted harzburgitic layer interpreted as oceanic or arc-related lithosphere emplaced during early tectonism. The CSMC lies wholly within the NE–SW striking central zone defined by Grütter et al. [Grütter, H.S., Apter, D.B., Kong, J., 1999. Crust–mantle coupling; evidence from mantle-derived xenocrystic garnets. Contributed paper at: The 7th International Kimberlite Conference Proceeding, J.B. Dawson Volume, 1, 307–313] on the basis of garnet geochemistry (G10 vs. G9) populations.

Deep-probing MT data from the lake bottom instruments infer that the conductor has a total depth-integrated conductivity (conductance) of the order of 2000 Siemens, which, given an internal resistivity of 10–15 Ω m, implies a thickness of 20–30 km. Below the CSMC the electrical resistivity of the lithosphere increases by a factor of 3–5 to values of around 50 Ω m. This change occurs at depths consistent with the graphite–diamond transition, which is taken as consistent with a carbon interpretation for the CSMC.

Preliminary three-dimensional MT modelling supports the NE–SW striking geometry for the conductor, and also suggests a NW dip. This geometry is taken as implying that the tectonic processes that emplaced this geophysical–geochemical body are likely related to the subduction of a craton of unknown provenance from the SE (present-day coordinates) during 2630–2620 Ma. It suggests that the lithospheric stacking model of Helmstaedt and Schulze [Helmstaedt, H.H., Schulze, D.J., 1989. Southern African kimberlites and their mantle sample: implications for Archean tectonics and lithosphere evolution. In Ross, J. (Ed.), Kimberlites and Related Rocks, Vol. 1: Their Composition, Occurrence, Origin, and Emplacement. Geological Society of Australia Special Publication, vol. 14, 358–368] is likely correct for the formation of the Slave's current SCLM.     

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.     

A mass transfer model of bauxite formation

The formation of bauxite due to weathering of a granitic protolith has been simulated by means of a one-dimensional flow and reaction model based on the mass transfer principle. The model couples mineral dissolution and precipitation reactions, speciation in solution, and advective solute transport in a porous medium. A very important aspect of the modeling study is the use of mineral reaction rates determined experimentally in the laboratory. The important effects of solution saturation state and pH have been incorporated into the kinetic rate laws governing the heterogeneous reactions. The values of these parameters have been obtained from the scientific literature to guarantee that realistic reaction rates are used in the simulations. Albite and quartz are the minerals that make up the parent rock in the model. Gibbsite, kaolinite, and a Na-mica (as a surrogate for smectite) are the secondary minerals that have been taken into account. Long-term simulations (>1 Ma) have been run, and the formation of a bauxitic profile, with an upper gibbsite-rich and a lower kaolinite-rich zone, is predicted. In early stages of the process (up to a few hundreds of thousands of years), both gibbsite and kaolinite precipitate directly from solution as a consequence of albite dissolution. In later stages, the bulk of gibbsite precipitation derives from the incongruent dissolution of kaolinite, while kaolinite precipitation is still caused by the dissolution of albite. This is also reflected by the formation of two reaction fronts in the profile. These results are compared with weathering sequences from the Los Pijiguaos bauxite deposit, Venezuela. The overlap between the gibbsite and kaolinite zones and the replacement of kaolinite by gibbsite are consistent with model calculations. Mechanical denudation has to be called upon to explain the limited thicknesses of the weathering profiles in the field. The role of mechanical erosion is supported by the presence of microsedimentary structures in the bauxite and the balance between dissolved and suspended loads in the streams draining the area.