Design of stable drainage networks in upland forestry plantations

Erosion of man-made, forestry drainage channels occurs when the plough cuts through the surface peat layer into the underlying erodible mineral soil. A procedure is developed, based on hydraulic considerations, which will allow the drainage engineer to design stable drainage networks in upland forestry plantations. An example design chart is given for an erodible sandy loam type soil.     

Isobaric cooling and anti-clockwise P–T paths in the Variscan Odenwald Crystalline Complex, Germany

Mineral assemblages in Al₂O₃‐rich, SiO₂‐ and K₂O‐poor metapelitic rocks from the western Odenwald Crystalline Complex (Variscan Mid‐German Crystalline Rise, southern Germany) include corundum, spinel, cordierite, sillimanite, garnet and staurolite. Quartz is absent from almost all samples. Therefore, the applicability of conventional geothermobarometry is very limited or even impossible. Detailed petrographic investigation on selected samples permits inference of the sequence of appearance and disappearance of several mineral assemblages. The recognition of such partial re‐equilibration stages and their associated mineral assemblages, together with mineral stabilities predicted from KFMASH pseudosections, enables the determination of the pressure‐temperature (P–T) trajectories experienced by these rocks during the Variscan metamorphism. The rocks were metamorphosed under low‐P/high‐T conditions and underwent an anti‐clockwise P–T evolution. A pressure increase from about 2 kbar to 4 ± 0.5 kbar was accompanied by heating. Peak metamorphic conditions were reached at pressures of 4 ± 0.5 kbar and temperatures of at least 640 °C, probably even higher. The retrograde evolution is characterised by near‐isobaric cooling from ≥ 640 °C to approximately 550 °C. The rocks underwent the anti‐clockwise evolution in a subduction‐related magmatic arc setting. The close spatial association of the low‐P/high‐T rocks with recently discovered metabasic eclogites in the eastern part of the Odenwald Crystalline Complex may indicate a fossil paired metamorphic belt in the Central European Variscides.     

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.     

A high-velocity ionized outflow and XUV photosphere in the narrow emission line quasar PG1211+143

We report on the analysis of a ∼60-ks XMM–Newton observation of the bright, narrow emission line quasar PG1211+143. Absorption lines are seen in both European Photon Imaging Camera and Reflection Grating Spectrometer spectra corresponding to H- and He-like ions of Fe, S, Mg, Ne, O, N and C. The observed line energies indicate an ionized outflow velocity of ∼24 000 km s⁻¹. The highest energy lines require a column density of   N _H∼ 5 × 10²³ cm⁻²  , at an ionization parameter of  log ξ∼ 3.4  . If the origin of this high-velocity outflow lies in matter being driven from the inner disc, then the flow is likely to be optically thick within a radius of ∼130 Schwarzschild radii, providing a natural explanation for the big blue bump (and strong soft X-ray) emission in PG1211+143.     

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