The Early Pleistocene Praetiglian and Ludhamian pollen stages in the North Sea Basin and their relationship to the marine isotope record

In this paper, the currently accepted correlation of the Early Pleistocene Ludhamian stage of England with the Tiglian‐A sub‐stage of the Netherlands is challenged. Recent investigations of Early Pleistocene marine North Sea deposits from a borehole near Noordwijk (the Netherlands) yielded evidence from molluscs, dinoflagellate cysts and sporomorphs for an alternation of warm‐temperate and arctic intervals within the Praetiglian and Tiglian stages. Marine equivalents of the terrestrial‐based pollen sub‐stages Tiglian A and B have been recognised in the upper part of the sequence. A Praetiglian age can be assigned to the lower part of the sequence on the basis of mollusc analysis. Within the Praetiglian, an alternation of warm and cold phases has been recognised from both the dinoflagellate cyst and molluscan records. Three cold phases within the Praetiglian are tentatively correlated with marine isotope stages (MIS) 96–100. The molluscan assemblages provide evidence for climate forcing of the sea level: highest sea levels are reached in the warm‐temperate intervals. Within the Praetiglian, an interval with an acme zone of the dinoflagellate cyst Impagidinium multiplexum, is correlated with the Ludhamian and tentatively linked to MIS 97 and/or MIS 96. The cold molluscan assemblages from the Noordwijk borehole include an acme zone of Megayoldia thraciaeformis, the first and only occurrence of this North Pacific bivalve in the North Sea Basin. Copyright © 2006 John Wiley & Sons, Ltd.     

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.     

Butterfly diagrams of strongly spotted late type stars

At present, long-term (over 30 years) multicolor photometric observations give the possibility to determine general properties of spotted areas on late-type stars. Star-spot modeling from broadband photometric data has been carried out by Alekseev and Gershberg since 1996 under the assumption that spots are situated in two latitudinal zones. Here we propose a new analysis of their results for several G and K dwarf stars with high irregular activity. On these stars, EK Dra, VY Ari, V775 Her, and V833 Tau, two spot belts exist separately and do not merge into a single equatorial active region, as occurs on cooler red-dwarf stars. The zonal spottedness models allow us to fit simultaneously both rotational modulation and long-term variability of stellar brightness. These models give evidence for an equatorward drift of the lower latitude boundary of the spotted region, φ₀, during the rising phase of activity, beyond any possible errors concerned with our methodology. In order to evaluate the drift rate we introduce the concept of `effective' spot belt, whose width is independent of longitude. This permits us to construct butterfly diagrams for stellar spots. The equatorward drift rates of the lower boundary of the spotted region D=dφ_low/dt are (− 1)–(− 2) deg year⁻¹ in the years of increasing spottedness. These values are less than the analogous solar one D≈−4 deg year⁻¹ for the rising phase of the cycle. Thus, cyclic activity can be revealed from butterfly diagrams and derived drifts of starspots prior to a possible detection from the spectral analysis of photometric variability. Finally, we briefly discuss a possible explanation of high-latitude activity and surface drifts of starspots in the framework of the current state of dynamo theory.     

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