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.     

A note on the problem of choosing a model of the Universe,IV

The Friedmannian model of the Universe with ultra-relativistic equation of statep=–1/3 withk=0 is nondecelerative, i.e., it expands with the constant velocityv=c during the whole expanding phase of evolution of the Universe. The ultrastable nondecelerative model is the only model of the Universe which fulfills all conditions resulting from the quantum-mechanical and relativistic analysis of the Universe.     

Methodology of hydrogeological characterization of deep carbonate aquifers as potential reservoirs of groundwater. Case of study: the Jurassic aquifer of <Emphasis Type="Italic">El Maestrazgo</Emphasis> (Castellón,Spain)

A methodology for the characterization of deep carbonate aquifers has been developed and applied to El Maestrazgo Jurassic aquifer in Castellón, Spain. Characterization of these aquifer formations, located at more than 300 m deep, consisted of a previous phase of compilation, analysis and synthesis of the existing information about the area, followed by a coordinated combination of different speciality studies: geology, stratigraphy, structural analysis, hydrogeology, hydrochemistry, geophysics and remote sensing. Geological studies included geological mapping, definition of stratigraphical units and facies and structural analysis. The aim of the hydrogeology study was to define aquifer formations, recharge area, aquifer points inventory and groundwater flow directions for the establishment of piezometric and water quality observation nets. Special techniques were applied, like thermal infrared aerial images and the evaluation of submarine groundwater discharge by means of natural radium isotopes. Hydrochemical techniques, including majority elements characterization and stable isotopes (¹⁸O, ²H and ³H) determination, allowed classifying hydrochemical facies and establishing a renewal pattern for water within the system. Geophysics was useful in determining the aquifer geometry, the features of the basement and the petrophysical characteristics of the geological formations. Preliminary results show an important tectonic complexity and the possibilities for groundwater uses in the area of study.     

Simulations of 3D DC borehole resistivity measurements with a goal-oriented <Emphasis Type="Italic">hp</Emphasis> finite-element method. Part II: through-casing resistivity instruments

We simulate direct current (DC) borehole resistivity measurements acquired in steel-cased deviated wells for the assessment of rock formation properties. The assumed data acquisition configuration considers one current (emitter) and three voltage (collector) electrodes that are utilized to measure the second difference of the electric potential along the well trajectory. We assume a homogeneous, 1.27-cm-thick steel casing with resistivity equal to 10^− 5 Ω· m. Simulations are performed with two different numerical methodologies. The first one is based on transferring two-dimensional (2D) axisymmetric optimal grids to a three-dimensional (3D) simulation software. The second one automatically produces optimal 3D grids yielded by a 3D self-adaptive goal-oriented algorithm. Both methodologies utilize high-order finite elements (FE) that are specially well-suited for problems with high-contrast coefficients and rapid spatial variations of the electric field, as it occurs in simulations that involve steel-cased wells. The method based on transferring 2D-optimal grids is efficient in terms of CPU time (few seconds per logging position). Unfortunately, it may produce inaccurate 3D simulations in deviated wells, even though the error remains below 1% for the axisymmetric (vertical) well. The method based on optimal 3D grids, although less efficient in terms of CPU time (few hours per logging position), produces more accurate results that are validated by a built-in a posteriori error estimator. This paper provides the first existing simulations of through-casing resistivity measurements in deviated wells. Simulated resistivity measurements indicate that, for a 30° deviated well, measurements in conductive layers 0.01 Ω· m) are similar to those obtained in vertical wells. However, in resistive layers (10,000 Ω· m), we observe 100% larger readings in the 30° deviated well. This difference becomes 3,000% for the case of a 60° deviated well. For this highly-deviated well, readings corresponding to the conductive formation layer are about 30% smaller in magnitude than those in a vertical well. Shoulder effects significantly vary in deviated wells.     

Some results of biogeographic research on meadows in the Dyje-Svratka basin in the Czech Republic

Some results of biogeographic research on meadows in the Dyje-Svratka basin in the Czech Republic. The paper shows one example of biogeographic inventory of natural ecosystems in an unflooded meadow in the Dyje-Svratka basin. It is part of the complex physio-geographic research of the Geographic Department of the Masaryk University in Brno aimed at gaining a basis for territorial interpretation of ecological regularities and territorial application of general ecological optimizing principles. A segment of biocenosis of Ulmi-fraxineta carpini has been analyzed (Zlatník 1956).The endeavour of Czech geographers to study geosystems of various dimensions completely gets in the limelight. In question is the understanding of mutual relationship of components of geosystems based on the exchange of mass and energy. There are two ways of studying them: spatial relations or approaches concentrated on processes. In contrast to the above unilateral concept, the endeavour to make the approaches compatible appears rarely. This complex view is made possible by geobiocenologically taken biogeography in papers of Zlatník and Rauer, their students and successors. This paper shows only one example of the results of biogeographic stock-taking of meadows in the Dyje-Svratka basin, indicating relations of eco- and geosystems which confirm the necessity of ecological orientation of the present biogeography and its potential contribution to ecology.     

Geochronology and geochemistry of eclogites from the Mariánské Lázně Complex,Czech Republic: Implications for Variscan orogenesis

The Mariánské Lázn complex (MLC) is located in the Bohemian Massif along the north-western margin of the Teplá-Barrandian microplate and consists of metagabbro, amphibolite and eclogite, with subordinate amounts of serpentinite, felsic gneiss and calcsilicate rocks. The MLC is interpreted as a metaophiolite complex that marks the suture zone between the Saxothuringian rocks to the north-west and the Teplá-Barrandian microplate to the south-east. Sm-Nd geochronology of garnet-omphacite pairs from two eclogite samples yields ages of 377±7, and 367±4 Ma. Samples of eclogite and amphibolite do not define a whole rock Sm-Nd isochron, even though there is a large range in Sm/Nd ratio, implying that the suite of samples may not be cogenetic. Eclogites do not have correlated _Nd values and initial ⁸⁷Sr/⁸⁶Sr ratios. Five of the eight eclogite samples have high _Nd values (+10.2 to +7.1) consistent with derivation from a MORB-like source, but variable ⁸⁷Sr/⁸⁶Sr ratios (0.7033 to 0.7059) which probably reflect hydrothermal seawater alteration. Three other eclogite samples have lower _Nd values (+ 5.4 to –0.8) and widely variable ⁸⁷Sr/⁸⁶Sr ratios (0.7033 to 0.7096). Such low _Nd values are inconsistent with derivation from a MORB, source and may reflect a subduction or oceanic island basalt component in their source. The MLC is an important petrotectonic element in the Bohemian Massif, providing evidence for Cambro-Ordovician formation of oceanic crust and interaction with seawater, Late Devonian (Frasnian-Famennian) high- and medium-pressure metamorphism related to closure of a Saxothuringian ocean basin, Early Carboniferous (Viséan) thrusting of the Teplá terrane over Saxothuringian rocks and Late Viséan extension.     

Dynamical instability in accreting white dwarfs

We study the evolution of solid, CO white dwarfs after explosive carbon ignition at central densities around 10¹⁰ g cm^–3 triggered by steady accretion in a close binary system, in order to elucidate whether these stars can collapse to form a neutron star. We show that as long as the velocity of the burning front remains below a critical value of 0.006c _s (60 km s^–1), gravitational collapse is the final fate. These calculations support the accretion-induced collapse (AIC) scenario for the origin of a fraction of low-mass X-ray binaries.Paper presented at the 11th European Regional Astronomical Meetings of the IAU on New Windows to the Universe, held 3–8 July, 1989, Tenerife, Canary Islands, Spain.     

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.