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.     

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.     

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.     

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     

Radon variations in a Hungarian village

 A steady radon exhalation is assumed in most publications. In a village of North-East Hungary, however, high radon concentrations have been measured, differing strongly in neighbouring houses and varying in time, due to the interplay of geochemical phenomena. Received: 20 November 1995 · Accepted: 18 June 1996     

River Warren boulders, Minnesota, USA: catastrophic paleoflow indicators in the southern spillway of glacial Lake Agassiz

Boulders resting on meltwater-sculpted and striated-granite bedrock near the head of the southern outlet spillway of glacial Lake Agassiz are used to generate paleodischarge calculations. The rounded nature of many boulders suggests fluvial transport and a corestone origin. The distribution of boulders in clusters and linear trains records the interaction of clasts during transport and deposition. The geomorphology of the spillway with streamlined erosional remnant hills is characteristic of other large flood spillways. Using the Manning equation and a variety of empirical equations to determine paleovelocity, preferred discharges between 0.364 and 0.102 Sv are calculated. These discharges agree well with flood discharges using modeling methodologies, and most likely represent ephemeral and catastrophic flood events linked to either episodic incision at the outlet or the result of rapid inputs of meltwater to Lake Agassiz.     

A technique for calibration of monostatic echosounder systems

A calibration technique has been adapted to render complete system calibrations of high-frequency acoustical instrumentation. This is based on standard targets; specifically, precisely manufactured spheres composed of tungsten carbide with 6% cobalt binder. The use of multiple sphere sizes was found to be advantageous, both as an independent check of the calibrations, and so that resonances in the sphere responses at certain frequencies could be avoided. Complete system gains and beam patterns, which include effects of bandpass filters and finite-pulse lengths, were determined by moving the spheres individually in the transducer far-fields. Use of this procedure ensures control over the acoustical characteristics of transducers, which may change from the time of manufacture and first testing due, for example, to platform mounting. It also provides a direct means of measuring the sampling volume at relatively high and constant signal-to-noise ratios. Implementation of this technique is discussed using a multifrequency sonar system as an example