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.     

OSCR wave measurements-some preliminary results

OSCR is an HF radar system that has been developed for high spatial resolution coastal surface current measurement. This paper describes preliminary results that demonstrate that wave measurement can be successfully obtained from suitably processed OSCR data. Comparisons with data from a WAVEC directional buoy are presented and show encouraging agreement. Some of the limitations to the measurement process are discussed and indicate a maximum range of about 20 km. Surface current variability on short time scales presents the most serious obstacle to wave measurement. This appears to be more of a problem when the mean currents are large, in that in these circumstances the data fail initial quality control criteria. However, in lower mean currents, the effect is often still present and leads to errors in long wave measurement     

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     

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.     

Infiltration at the Tai rain forest (Ivory Coast): Measurements and modelling

Infiltration experiments have been performed at three sites along a well-known catena under virgin tropical rain forest using a portable sprinkling infiltrometer. Experimentally determined infiltration curves are presented. Infiltration curves are also simulated on the basis of the Mein-Larson equation. The parameters for this model have been obtained from the infiltration curves (saturated conductivity) and simple soil moisture determinations (fillable porosity). The agreement between experimentally determined and modelled infiltration is reasonable, provided (a) saturated conductivity as derived from the experimental data is corrected, (b) a storage parameter, also derived from the experimental data, is added to the Mein-Larson model, and (c) the decline in soil porosity with depth is either small or occurs abruptly at shallow depth. Comparison of observed infiltration rates with rainfall intensity shows that Horton Overland Flow has to occur naturally at least on the middle and lower section of the catena. Despite the fact that most parameters can be estimated in principle from basic soil data, it remains advisable to obtain sprinkling infiltrometer field measurements, because of soil variability due to dynamic surface conditions, macroporosity, air entrapment, and irregularity of the wetting front.     

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.     

Wavelet analysis of inter‐annual variability in the runoff regimes of glacial and nival stream catchments,Bow Lake,Alberta

Continuous wavelet analyses of hourly time series of air temperature, stream discharge, and precipitation are used to compare the seasonal and inter‐annual variability in hydrological regimes of the two principal streams feeding Bow Lake, Banff National Park, Alberta: the glacial stream draining the Wapta Icefields, and the snowmelt‐fed Bow River. The goal is to understand how water sources and flow routing differ between the two catchments. Wavelet spectra and cross‐wavelet spectra were determined for air temperature and discharge from the two streams for summers (June–September) 1997–2000, and for rainfall and discharge for the summers of 1999 and 2000. The diurnal signal of the glacial runoff was orders of magnitude higher in 1998 than in other years, indicating that significant ice exposure and the development of channelized glacial drainage occurred as a result of the 1997–98 El Niño conditions. Early retreat of the snowpack in 1997 and 1998 led to a significant summer‐long input of melt runoff from a small area of ice cover in the Bow River catchment; but such inputs were not apparent in 1999 and 2000, when snow cover was more extensive. Rainfall had a stronger influence on runoff and followed quicker flow paths in the Bow River catchment than in the glacial catchment. Snowpack thickness and catchment size were the primary controls on the phase relationship between temperature and discharge at diurnal time scales. Wavelet analysis is a fast and effective means to characterize runoff, temperature, and precipitation regimes and their interrelationships and inter‐annual variability. The technique is effective at identifying inter‐annual and seasonal changes in the relative contributions of different water sources to runoff, and changes in the time required for routing of diurnal meltwater pulses through a catchment. However, it is less effective at identifying changes/differences in the type of the flow routing (e.g. overland flow versus through flow) between or within catchments. Copyright © 2003 John Wiley & Sons, Ltd.     

Swimming with…Devonian fish

Armoured jawless fish, or 'ostracoderms', lived 450–360 million years ago, and display unusual morphologies, unlike any modern fish group. Since they left no living descendants, their mode of swimming has, until recently, remained speculative, although this is a crucial question as the first true pectoral fins evolved within the 'ostracoderms'. The discovery of the oldest-known fish trails, from the Early Devonian (400 million year old) 'Lower Old Red Sandstone' of south-east Wales offers new insights into the swimming behaviour of these early fish, notably the osteostracan 'ostracoderms' (or cephalaspids), whose horseshoe-shaped head and paddle-shaped pectoral fins have remained a functional riddle.     

North European last glacial–interglacial transition (LGIT; 15–9 ka) tephrochronology: extended limits and new events

High‐precision correlation of palaeoclimatic and palaeoenvironmental records is crucial for testing hypotheses of synchronous change. Although radiocarbon is the traditional method for dating late Quaternary sedimentary sequences, particularly during the last glacial–interglacial transition (LGIT; 15–9 ka), there are inherent problems with the method, particularly during periods of climate change which are often accompanied by major perturbations in atmospheric radiocarbon content. An alternative method is the use of tephras that act as time‐parallel marker horizons. Within Europe, numerous volcanic centres are known to have erupted during the LGIT, providing considerable potential for high‐precision correlation independent of past radiocarbon fluctuations. Here we report the first identification of the Vedde Ash and Askja Tephra in Ireland, significantly extending the known provenance of these events. We have also identified two new horizons (the Roddans Port Tephras A and B) and tentatively recognise an additional horizon from Vallensgård Mose (Denmark) that provide crucial additional chronological control for the LGIT. Two phases of the Laacher See Tephra (LST) are reported, the lower Laacher See Tephra (LLST) and probably the C2 phase of the Middle Laacher See Tephra (MLST‐C2) indicating a more northeasterly distribution of this fan than reported previously. Copyright © 2006 John Wiley & Sons, Ltd.