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.     

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.     

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.     

Dauphiné twinning and texture memory in polycrystalline quartz. Part 3: texture memory during phase transformation

Samples of quartz-bearing rocks were heated above the α (trigonal)–β (hexagonal) phase transformation of quartz (625–950°C) to explore changes in preferred orientation patterns. Textures were measured both in situ and ex situ with neutron, synchrotron X-ray and electron backscatter diffraction. The trigonal–hexagonal phase transformation does not change the orientation of c- and a-axes, but positive and negative rhombs become equal in the hexagonal β-phase. In naturally deformed quartzites measured by neutron diffraction a perfect texture memory was observed, i.e. crystals returned to the same trigonal orientation they started from, with no evidence of twin boundaries. Samples measured by electron back-scattered diffraction on surfaces show considerable twinning and memory loss after the phase transformation. In experimentally deformed quartz rocks, where twinning was induced mechanically before heating, the orientation memory is lost. A mechanical model can explain the memory loss but so far it does not account for the persistence of the memory in quartzites. Stresses imposed by neighboring grains remain a likely cause of texture memory in this mineral with a very high elastic anisotropy. If stresses are imposed experimentally the internal stresses are released during the phase transformation and the material returns to its original state prior to deformation. Similarly, on surfaces there are no tractions and thus texture memory is partially lost.     

Groundwater contamination with arsenic and other trace elements in an area of the pampa,province of Córdoba,Argentina

A geochemical study of groundwater of the pampa in the province of Córdoba, Argentina, was performed; the area covered approximately 10,000 km².Physical-chemical parameters, dissolved solids, and seven trace elements were determined in 60 selected water samples. Systematic and accurate measurements of arsenic, flourine, and vanadium were performed for the first time. Three trace element contaminants not reported earlier were found: an important one, selenium, and two others of less known effects, uranium and molybdenum.Eighty-four percent of the water analyzed showed arsenic contents over 0.05 mg/L, maximum contaminant level established by the U.S. Environmental Protection Agency (1982). The frequency distribution of trace elements was analyzed, and its fit to the lognormal distribution was proved by means of the Pearson and Kolmogorov-Smirnov test; the geographic distribution of the seven trace elements was mapped and its correlation with the anion-cation composition of the water was studied.The maximum arsenic, fluorine, vanadium, and uranium contents were found in the western part of the area under study, in waters containing dominant alkali metals in the cation composition. Maximum selenium and antimony contents were found in the eastern part of the area, while molybdenum distribution does not show any relationship with the other two groups. In addition, the geographic distribution of the trace elements seems to be related to the subsurface structure, which has been inferred using interactive digital analysis of Landsat imagery. The movements of the subsoil have disturbed surface and subsurface drainage influencing the water salinity and trace element contents.In order to investigate the origin of the contamination, 54 loess samples were collected in wells at depths ranging from the surface down to the water table. This loess, which has a high proportion of volcanic components, mainly rhyolitic glass, exhibits a chemical composition corresponding to that of a dacite.The loess and the volcanic glass show anomalous contents of all contaminant trace elements, mainly arsenic and selenium. For this reason loess is considered the most important contamination source in the groundwater under study.