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.     

Evidence for the loss of snow-deposited MSA to the interstitial gaseous phase in central Antarctic firn

We have examined several MSA (methanesulfonic acid) records from the upper 200 m of the Antarctic ice sheet and in particular the new Dome F profile. At all the four sites studied, concentration profiles exhibit similar patterns as a function of depth. They suggest that snow metamorphism and solid phase migration are responsible for a marked release of gaseous MSA to interstitial firn air as well as probably to the free atmosphere, in particular at extremely low accumulation sites. Snow acidity can also modify MSA concentration. It is proposed that, below the upper few metres where the communication with the free atmosphere is possible, gaseous MSA may remain in the firn layers and be entrapped later in air bubbles at pore close-off, i.e. when firn is transformed into ice. Chemical measurements on the firn core do not take into account the MSA released to the gaseous phase, but this fraction is measurable in ice samples. In spite of these alterations occurring in the firn layers, relative changes of the atmospheric MSA concentration in the past are probably still there deep within the Antarctic ice sheet. However, for glacial periods, different processes have to be considered in relation to modified aerosol properties.     

Chemical alteration of tephra in the depositional environment: theoretical stability modelling

The study of the chemical stability of vitreous material in aqueous media is well‐established. There has to date been little consideration of the implications of variations in the chemical durability of tephra in Quaternary tephrochronology. Chemical alteration can take the form of cationic leaching from the matrix, or complete destruction of the silica network, either of which could constrain the ability to chemically identify distal tephra. Here we apply established models of vitreous durability to the published chemical analyses of a large number of Icelandic tephras in order to predict their relative durabilities under equivalent conditions. This suggests that some important tephras have relatively poor chemical stability, and that rhyolitic tephras are, in general, more stable than basaltic. We conclude that tephras should be expected to show predictable differential chemical stability in the post‐depositional environment. Copyright © 2003 John Wiley & Sons, Ltd.     

Mezen Bay—a key area for understanding Weichselian glaciations in northern Russia

Sediment successions in coastal cliffs around Mezen Bay, southeastern White Sea, record an unusually detailed history of former glaciations, interstadial marine and fluvial events from the Weichselian. A regional glaciation model for the Weichselian is based on new data from the Mezen Bay area and previously published data from adjacent areas. Following the Mikulinian (Eemian) interglacial a shelf‐centred glaciation in the Kara Sea is reflected in proglacial conditions at 100–90 ka. A local ice‐cap over the Timan ridge existed between 75 and 65 ka. Renewed glaciation in the Kara Sea spread southwestwards around 60 ka only, interrupted by a marine inundation, before it advanced to its maximum position at about 55–50 ka. After a prolonged ice‐free period, the Scandinavian ice‐sheet invaded the area from the west and terminated east of Mezen Bay about 17 ka. The previously published evidence of a large ice‐dammed lake in the central Arkhangelsk region, Lake Komi, finds no support in this study. Copyright © 2003 John Wiley & Sons, Ltd.     

Identification of physical processes inherent in artificial neural network rainfall runoff models

The emergence of artificial neural network (ANN) technology has provided many promising results in the field of hydrology and water resources simulation. However, one of the major criticisms of ANN hydrologic models is that they do not consider/explain the underlying physical processes in a watershed, resulting in them being labelled as black‐box models. This paper discusses a research study conducted in order to examine whether or not the physical processes in a watershed are inherent in a trained ANN rainfall‐runoff model. The investigation is based on analysing definite statistical measures of strength of relationship between the disintegrated hidden neuron responses of an ANN model and its input variables, as well as various deterministic components of a conceptual rainfall‐runoff model. The approach is illustrated by presenting a case study for the Kentucky River watershed. The results suggest that the distributed structure of the ANN is able to capture certain physical behaviour of the rainfall‐runoff process. The results demonstrate that the hidden neurons in the ANN rainfall‐runoff model approximate various components of the hydrologic system, such as infiltration, base flow, and delayed and quick surface flow, etc., and represent the rising limb and different portions of the falling limb of a flow hydrograph. Copyright © 2004 John Wiley & Sons, Ltd.     

From harmony of the spheres to national anthem: Reflections on musical heritage

Western classical music was long distinctive among the arts for its universality and catholicity. From Greek antiquity on, music uniquely reflected Orphic, mystic, celestial feeling, And unlike painting, architecture, or literature, the standard musical canon came to be breathtakingly inclusive over time as well as space. Yet in recent centuries music has increasingly celebrated national patriotism and period particularism. This essay traces the histories of these antithetical expressive urges, the one universalistic, transcendental, atemporal, the other chauvinist, parochial, historicist.     

Response of circular footings and anchor plates in non-homogeneous elastic soils

The axisymmetric elastic response of circular footings and anchor plates in a linearly non-homogeneous elastic soil is analysed. It is assumed that footings/anchors are flexible and subjected to axisymmetric vertical loads. The response of the footing/anchor is modelled by using the classical Poisson–Kirchhoff thin plate theory. A variational technique is used to analyse the interaction problem. A representation for the contact stress is established by using a fundamental solution corresponding to a unit vertical pressure acting over an annular region in the interior of the non-homogeneous soil. The fundamental solution can be derived by using rigorous analytical procedures. The influence of the footing flexibility and the degree of soil non-homogeneity on the displacements, bending moments and contact stresses of a surface footing is examined over a wide range of governing parameters. In the case of anchor plates the influence of depth of embedment, degree of soil non-homogeneity and anchor flexibility on the anchor displacement is investigated.