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.     

Chemical abundances determined from meteor spectra: I. Ratios of the main chemical elements

Abstract— Relative chemical abundances of 13 meteoroids were determined by averaging the composition of the radiating gas along the fireball path that originated during their penetration into the Earth's atmosphere. Mg, Fe, Ni, Cr, Mn, and Co abundances, relative to Si, are similar to those reported for CI and CM carbonaceous chondrites and interplanetary dust particles. In contrast, relative abundances of Ca and Ti in meteor spectra indicate that these elements suffer incomplete evaporation processes. The chemical composition of all meteoroids studied in this work differs from that of 1P/Halley dust.     

Coupled hydrogeological and reactive transport modelling of the Simpevarp area (Sweden)

The Simpevarp area is one of the alternative sites being considered for the deep geological disposal of high level radioactive waste in Sweden. In this paper, a coupled regional groundwater flow and reactive solute transport model of the Simpevarp area is presented that integrates current hydrogeological and hydrochemical data of the area. The model simulates the current hydrochemical pattern of the groundwater system in the area. To that aim, a conceptual hydrochemical model was developed in order to represent the dominant chemical processes. Groundwater flow conditions were reproduced by taking into account fluid-density-dependent groundwater flow and regional hydrogeologic boundary conditions. Reactive solute transport calculations were performed on the basis of the velocity field so obtained. The model was calibrated and sensitivity analyses were carried out in order to investigate the effects of heterogeneities of hydraulic conductivity in the subsurface medium. Results provided by the reactive transport model are in good agreement with much of the measured hydrochemical data. This paper emphasizes the appropriateness of the use of reactive solute transport models when water-rock interaction reactions are involved, and demonstrates what powerful tools they are for the interpretation of hydrogeological and hydrochemical data from site geological repository characterization programs, by providing a qualitative framework for data analysis and testing of conceptual assumptions in a process-oriented approach.     

Optimizing interpolation of shoot density data from a Posidonia oceanica seagrass bed

A case study on the optimization of Posidonia oceanica density interpolation, using a data set from a large meadow at Porto Conte Bay (NW Sardinia, Italy), is presented. Ordinary point kriging, cokriging and a weighted average based on inverse square distance were used to interpolate density data measured in 36 sampling stations. The results obtained from different methods were then compared by means of a leave‐one‐out cross‐validation procedure. The scale at which interpolation was carried out was defined on the basis of the Hausdorff dimension of the variogram. Optimizing spatial scale and data points search strategy allowed obtaining more accurate density estimates independently of the interpolation method.     

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.     

Modeling the seismic response of basin-edge structure by a sequential waveform method

A sequential waveform method is developed to simulate the seismic response of basin-edge structure excited by a plane incident P-wave. The full procedure involves: (a) a previous parameterization of the investigated model using the seismic wave velocities and depths of the sedimentary stratifications; (b) an input motion determined from the records at stations installed on hard rock; (c) forward computation of the P-SV elastic wave field by means of a two-dimensional finite difference (FD) method; (d) the optimization of the model vector using simulated annealing technique and comparing the simulated seismic response of the tested structure with the observed wave field; (e) the correction of the initial model by trial-and-error by testing the differences between synthetics and observed data, and (f) the final solution obtained by iteration using the conjugate gradient algorithm. The search of an optimal basin-edge model has been parallel processed by varying the shapes and velocities of strata on the basis of the fitting of relative timing, amplitude and phase between the output and the observed data. The input motion and sensitivity have been checked and the validity of the method has been demonstrated by numeric analysis. Using the teleseismic records generated by 7 earthquakes recorded at 26 broadband seismic stations, we have studied the seismic velocity structure of the southern edge of the Jiyang depression located in the Bohai Bay basin, northern China. Two cross sections show an agreement between the velocity results and the geological sections available in the region. In addition, we obtain evidence of three hidden faults under the sections and features that suggest major extensions at the Paleogene.     

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.