A new stereo‐analytical method for determination of removal blocks in discontinuous rock masses

The paper provides a new stereo‐analytical method, which is a combination of the stereographic method and analytical methods, to separate finite removable blocks from the infinite and tapered blocks in discontinuous rock masses. The methodology has applicability to both convex and concave blocks. Application of the methodology is illustrated through examples. Addition of this method to the existing block theory procedures available in the literature improves the capability of block theory in solving practical problems in rock engineering. Copyright © 2003 John Wiley & Sons, Ltd.     

Predicting the polarization of the microwave background from the WMAP temperature maps

In this paper, we study how to predict the polarization of the Cosmic Microwave Background (CMB) using knowledge of only the temperature (intensity) and the cross-correlation between temperature and polarization. We derive a “Wiener prediction” method and apply it to the Wilkinson Microwave Anisotropy Probe (WMAP) all-sky CMB temperature maps and to the MAXIMA field.     

A micromechanical finite element model for linear and damage‐coupled viscoelastic behaviour of asphalt mixture

This study presents a finite element (FE) micromechanical modelling approach for the simulation of linear and damage‐coupled viscoelastic behaviour of asphalt mixture. Asphalt mixture is a composite material of graded aggregates bound with mastic (asphalt and fine aggregates). The microstructural model of asphalt mixture incorporates an equivalent lattice network structure whereby intergranular load transfer is simulated through an effective asphalt mastic zone. The finite element model integrates the ABAQUS user material subroutine with continuum elements for the effective asphalt mastic and rigid body elements for each aggregate. A unified approach is proposed using Schapery non‐linear viscoelastic model for the rate‐independent and rate‐dependent damage behaviour. A finite element incremental algorithm with a recursive relationship for three‐dimensional (3D) linear and damage‐coupled viscoelastic behaviour is developed. This algorithm is used in a 3D user‐defined material model for the asphalt mastic to predict global linear and damage‐coupled viscoelastic behaviour of asphalt mixture. For linear viscoelastic study, the creep stiffnesses of mastic and asphalt mixture at different temperatures are measured in laboratory. A regression‐fitting method is employed to calibrate generalized Maxwell models with Prony series and generate master stiffness curves for mastic and asphalt mixture. A computational model is developed with image analysis of sectioned surface of a test specimen. The viscoelastic prediction of mixture creep stiffness with the calibrated mastic material parameters is compared with mixture master stiffness curve over a reduced time period. In regard to damage‐coupled viscoelastic behaviour, cyclic loading responses of linear and rate‐independent damage‐coupled viscoelastic materials are compared. Effects of particular microstructure parameters on the rate‐independent damage‐coupled viscoelastic behaviour are also investigated with finite element simulations of asphalt numerical samples. Further study describes loading rate effects on the asphalt viscoelastic properties and rate‐dependent damage behaviour. Copyright © 2006 John Wiley & Sons, Ltd.     

Rapid vegetation change during the early Holocene in the Faroe Islands detected in terrestrial and aquatic ecosystems

High‐resolution pollen, plant macrofossil and sedimentary analyses from early Holocene lacustrine sediments on the Faroe Islands have detected a significant vegetation perturbation suggesting a rapid change in climate between ca. 10 380 cal. yr BP and the Saksunarvatn ash (10 240±60 cal. yr BP). This episode may be synchronous with the decline in δ¹⁸O values in the Greenland ice‐cores. It also correlates with a short, cold event detected in marine cores from the North Atlantic that has been ascribed to a weakening of thermohaline circulation associated with the sudden drainage of Lake Agassiz into the northwest Atlantic, or, alternatively, a period with distinctly decreased solar forcing. The vegetation sequence begins at ca. 10 500 cal. yr BP with a succession from tundra to shrub‐tundra and increasing lake productivity. Rapid population increases of aquatic plants suggest high summer temperatures between 10 450 and 10 380 cal. yr BP. High pollen percentages, concentrations and influx of Betula, Juniperus and Salix together with macrofossil leaves indicate shrub growth around the site during the initial phases of vegetation colonisation. Unstable conditions followed ca. 10 380 cal. yr BP that changed both the upland vegetation and the aquatic plant communities. A decrease in percentage values of shrub pollen is recorded, with replacement of both aquatics and herbaceous plants by pioneer plant communities. An increase in total pollen accumulation rates not seen in the concentration data suggests increased sediment delivery. The catchment changes are consistent with less seasonal, moister conditions. Subsequent climatic amelioration reinitiated a warmth‐driven succession and catchment stabilisation, but retained high precipitation levels influencing the composition of the post‐event communities. Copyright © 2003 John Wiley & Sons, Ltd.     

SANICLAY: simple anisotropic clay plasticity model

SANICLAY is a new simple anisotropic clay plasticity model that builds on a modification of an earlier model with an associated flow rule, in order to include simulations of softening response under undrained compression following K_o consolidation. Non‐associativity is introduced by adopting a yield surface different than the plastic potential surface. Besides, the isotropic hardening of the yield surface both surfaces evolve according to a combined distortional and rotational hardening rule, simulating the evolving anisotropy. Although built on the general premises of critical state soil mechanics, the model induces a critical state line in the void ratio–mean effective stress space, which is a function of anisotropy. To ease interpretation, the model formulation is presented firstly in the triaxial stress space and subsequently, its multiaxial generalization is developed systematically, in a form appropriate for implementation in numerical codes. The SANICLAY is shown to provide successful simulation of both undrained and drained rate‐independent behaviour of normally consolidated sensitive clays, and to a satisfactory degree of accuracy of overconsolidated clays. The new model requires merely three constants more than those of the modified Cam clay model, all of which are easily calibrated from well‐established laboratory tests following a meticulously presented procedure. Copyright © 2006 John Wiley & Sons, Ltd.     

Intercomparison of the seasonal cycle of tropical surface stress in 17 AMIP atmospheric general circulation models

 The mean state of the tropical atmosphere is important as the nature of the coupling between the ocean and the atmosphere depends nonlinearly on the basic state of the coupled system. The simulation of the annual cycle of the tropical surface wind stress by 17 atmospheric general circulation models (AGCMs) is examined and intercompared. The models considered were part of the Atmospheric Model Intercomparison Project (AMIP) and were integrated with observed sea surface temperature (SST) for the decade 1979–1988. Several measures have been devised to intercompare the performance of the 17 models on global tropical as well as regional scales. Within the limits of observational uncertainties, the models under examination simulate realistic tropical area-averaged zonal and meridional annual mean stresses. This is a noteworthy improvement over older generation low resolution models which were noted for their simulation of surface stresses considerably weaker than the observations. The models also simulate realistic magnitudes of the spatial distribution of the annual mean surface stress field and are seen to reproduce realistically its observed spatial pattern. Similar features are observed in the simulations of the annual variance field. The models perform well over almost all the tropical regions apart from a few. Of these, the simulations over Somali are interesting. Over this region, the models are seen to underestimate the annual mean zonal and meridional stresses. There is also wide variance between the different models in simulating these quantities. Large model-to-model variations were also seen in the simulations of the annual mean meridional stress field over equatorial Indian Ocean, south central Pacific, north east Pacific and equatorial eastern Pacific oceans. It is shown that the systematic errors in simulating the surface winds are related to the systematic errors in simulating the Inter-Tropical Convergence Zone (ITCZ) in its location and intensity. Weaker than observed annual mean southwesterlies simulated by most models over Somali is due to weaker than observed southwesterlies during the Northern Hemisphere summer. This is related to the weaker than observed land precipitation simulated by most models during the Northern Hemisphere summer. The diversity in simulation of the surface wind over Somali and equatorial Indian ocean is related to the diversity of AGCMs in simulating the precipitation zones in these regions. Received: 2 August 1996 / Accepted: 7 February 1997     

WEKU – a GIS-Supported stochastic model of groundwater residence times in upper aquifers for the supraregional groundwater management

 The supraregional GIS-supported stochastical model, WEKU, for the determination of groundwater residence times in the upper aquifers of large groundwater provinces is presented. Using a two-dimensional analytical model of groundwater flow, groundwater residence times are determined within two extreme cases. In the first case, maximal groundwater residence times are calculated, representing the part of groundwater, that is drained by the main surface water of a groundwater catchment area. In the second case, minimal groundwater residence times for drainage into the nearest surface water are determined. Using explicit distribution functions of the input parameters, mean values as well as potential ranges of variations of the groundwater residence times are derived. The WEKU model has been used for the determination of groundwater residence times throughout Germany. The model results – mean values and deviations of the groundwater velocity and the maximal and minimal groundwater residence times in the upper aquifers – are presented by general maps and discussed in detail. It is shown that the groundwater residence times in the upper aquifer vary regionally, differentiated between less than 1 year and more than 2000 years. Using this information, the time scales can be specified, until measures to remediate polluted groundwater resources may lead to a substantial groundwater quality improvement in the different groundwater provinces of Germany. With respect to its supraregional scale of application, the WEKU model may serve as a useful tool for the supraregional groundwater management on a state, federal or international level. Received: 15 August 1995 · Accepted: 15 October 1995