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.     

The pathway of mercury in contaminated waters determined by association with organic carbon (Tagus Estuary,Portugal)

Reactive dissolved Hg (Hg_R), non-reactive dissolved Hg (Hg_NR), particulate Hg (Hg_P), dissolved organic C (DOC), particulate organic C (POC), salinity and other interpretative parameters were determined in water samples collected in the North Channel and in adjacent areas of the Tagus estuary (Portugal). Higher concentrations of both dissolved and particulate Hg in the North Channel indicate a pollution source and raise the possibility of Hg escaping to adjacent areas by tidal action. This transport was confirmed by the increase of Hg_R with salinity and Hg_NR with DOC, along a longitudinal axis paralleling the North Channel. Apparently, Hg leaving this channel is progressively complexed by inorganic and organic ligands. Near the mouth of the estuary, values decreased reflecting dilution with seawater. Moreover the Hg_P:POC ratio also increased seaward, suggesting mixing with Hg enriched particles that escaped the North Channel, or incorporation of dissolved Hg species in river-derived particles. These results suggest that the pathway of anthropogenic Hg in contaminated waters may be identified by their enrichment in organic matter, both in the dissolved and particulate fraction.     

Nutrients and zooplankton composition and dynamics in relation to the hydrological pattern in a confined Mediterranean salt marsh (NE Iberian Peninsula)

The role of the hydrological regime in the nutrients and zooplankton composition and dynamics has been analysed in five lagoons of La Pletera salt marshes (NE Iberian Peninsula) during a complete hydrological cycle (2002–2003). Two of the lagoons have their origin in the old river mouths while the other three were recently created in the framework of a Life Restoration project. This fact has also allowed us to study the effect of the lagoon age on nutrient and zooplankton composition and dynamics. The salt marsh hydrology is determined by a prolonged period of confinement without water inputs, irregularly interrupted by sudden water inputs due to flooding events (sea storms or intense rainfalls). While the dynamics of oxidized nitrogen compounds in the lagoons depends on the water inputs variability within each hydrological cycle, the internal load of phosphorus, total nitrogen and organic matter is related more to the cumulative mechanisms during the confinement periods. Accumulation processes may be easily related to lagoon age, since old lagoons have higher content of nutrients and organic matter, suggesting that these lagoons progressively accumulate nutrients during the successive confinement events. This is the usual case for most Mediterranean salt marshes without an artificially manipulated water regime. The zooplankton community in La Pletera integrates the effects of both the hydrological regime and the lagoon age since the former determines the temporal pattern of the main zooplankton species and the latter explains differences in composition and structure between old and new lagoons.     

Effects of crab–halophytic plant interactions on creek growth in a S.W. Atlantic salt marsh: A Cellular Automata model

The Bahía Blanca Estuary (38° 50′ S, and 62° 30′ W) presents salt marshes where interactions between the local flora (Sarcocornia perennis) and fauna (Chasmagnathus granulatus) generate some kind of salt pans that alter the normal water circulation and condition its flow and course towards tidal creeks. The crab–vegetation dynamics in the salt marsh presents variations that cannot be quantified in a reasonable period of time. The interaction between S. perennis plant and C. granulatus crab is based on simple laws, but its result is a complex biological mechanism that causes an erosive process on the salt marsh and favors the formation of tidal creeks. To study it, a Cellular Automata model is proposed, based on the laws deduced from the observation of these phenomena in the field, and then verified with measurable data within macroscale time units. Therefore, the objective of this article is to model how the interaction between C. granulatus and S. perennis modifies the landscape of the salt marsh and influences the path of tidal creeks. The model copies the basic laws that rule the problem based on purely biological factors.The Cellular Automata model proved capable of reproducing the effects of the interaction between plants and crabs in the salt marsh. A study of the water drainage of the basins showed that this interaction does indeed modify the development of tidal creeks. Model dynamics would likewise follow different laws, which would provide a different formula for the probability of patch dilation. The patch shape can be obtained changing the pattern that dilates.     

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.