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 simulation tool for high data-rate acoustic communication in ashallow-water, time-varying channel

The paper discusses the development of a simulation tool to model high data-rate acoustic communication in shallow water. The simulation tool is able to generate synthetic time series of signals received at a transducer array after transmission across a shallow-water communication channel. The simulation tool is suitable for testing advanced signal processing techniques for message recovery. A channel model has been developed based on the physical aspects of the acoustic channel. Special emphasis has been given to fluctuations of the signal transmission caused by time-varying multipath effects. At shorter ranges, the temporal variations are dominated by acoustic scattering from the moving sea surface. Therefore, the channel model produces a coherence function which may be interpreted as a time-varying reflection coefficient for the surface scattered acoustical path. A static, range-independent ray model identifies the significant multipaths, and the surface path is modulated with the time-varying reflection coefficient. The advantages and limitations of the channel model are discussed and assumptions necessary to overcome the limitations are emphasised. Based on the assumptions, an algorithm has been developed and implemented to model how a binary message will be modulated when transmitted by a transducer, is distorted in the channel and finally is received by a transducer array     

Application of an improved self-starter to geoacoustic inversion

The self-starter is improved using the operator of the split-step Pade solution. In addition to providing greater stability and being applicable closer to the source, the improved self-starter is an efficient forward model for geoacoustic inversion. It is necessary to solve only O(10) tridiagonal systems of equations to obtain the acoustic field on a vertical array located O(10) wavelengths from a source. This experimental configuration is effective for geoacoustic inverse problems involving unknown parameters deep in the ocean bottom. For problems involving depth-dependent acoustic parameters, the improved self-starter can be used to solve nonlinear inverse problems involving O(10) unknown sediment parameters in less than a minute on the current generation of workstations     

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.     

Ice‐cover variability on shallow lakes at high latitudes: model simulations and observations

A one‐dimensional thermodynamic model for simulating lake‐ice phenology is presented and evaluated. The model can be driven with observed daily or hourly atmospheric forcing of air temperature, relative humidity, wind speed, cloud amount and snowfall. In addition to computing the energy balance components, key model output includes the temperature profile at an arbitrary number of levels within the ice/snow (or the water temperature if there is no ice) and ice thickness (clear ice and snow‐ice) on a daily basis, as well as freeze‐up and break‐up dates. The lake‐ice model is used to simulate ice‐growth processes on shallow lakes in arctic, sub‐arctic, and high‐boreal forest environments. Model output is compared with field and remote sensing observations gathered over several ice seasons. Simulated ice thickness, including snow‐ice formation, compares favourably with field measurements. Ice‐on and ice‐off dates are also well simulated when compared with field and satellite observations, with a mean absolute difference of 2 days. Model simulations and observations illustrate the key role that snow cover plays on the seasonal evolution of ice thickness and the timing of spring break‐up. It is also shown that lake morphometry, depth in particular, is a determinant of ice‐off dates for shallow lakes at high latitudes. Copyright © 2003 John Wiley & Sons, Ltd.     

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.