Relative sea‐level observations in western Scotland since the Last Glacial Maximum for testing models of glacial isostatic land movements and ice‐sheet reconstructions

Observations of relative sea‐level change and local deglaciation in western Scotland provide critical constraints for modelling glacio‐isostatic rebound in northern Britain over the last 18 000 years. The longest records come from Skye, Arisaig and Knapdale with a shorter, Holocene, record from Kintail. Biostratigraphic (diatom, pollen, dinoflagellate, foraminifera and thecamoebian), lithological and radiocarbon analyses provide age and elevation parameters for each sea‐level index point. All four sites reveal relative sea‐level change that is highly non‐monotonic in time as the local vertical component of glacio‐isostatic rebound and eustasy (or global meltwater influx) dominate at different periods. Copyright © 2006 John Wiley & Sons, Ltd.     

Weichselian sediments at Foss‐Eikeland,Jæren (southwest Norway): sea‐level changes and glaciation history

The Jæren area in southwestern Norway has experienced great changes in sea‐levels and sedimentary environments during the Weichselian, and some of these changes are recorded at Foss‐Eikeland. Four diamictons interbedded with glaciomarine and glaciofluvial sediments are exposed in a large gravel pit situated above the post‐glacial marine limit. The interpretation of these sediments has implications for the history of both the inland ice and the Norwegian Channel Ice Stream. During a Middle Weichselian interstadial, a large glaciofluvial delta prograded into a shallow marine environment along the coast of Jæren. A minor glacial advance deposited a gravelly diamicton, and a glaciomarine diamicton was deposited during a following marine transgression. This subsequently was reworked by grounded ice, forming a well‐defined boulder pavement. The boulder pavement is followed by glaciomarine clay with a lower, laminated part and an upper part of sandy clay. The laminated clay probably was deposited under sea‐ice, whereas more open glaciomarine conditions prevailed during deposition of the upper part. The clay is intersected by clastic dykes protruding from the overlying, late Weichselian till. Preconsolidation values from the marine clay suggest an ice thickness of at least 500 m during the last glacial phase. The large variations in sea‐level probably are a combined effect of eustasy and glacio‐isostatic changes caused by an inland ice sheet and an ice stream in the Norwegian Channel. Copyright © 2002 John Wiley & Sons, Ltd.     

Finite element analysis of land subsidence above depleted reservoirs with pore pressure gradient and total stress formulations

The solution of the poroelastic equations for predicting land subsidence above productive gas/oil fields may be addressed by the principle of virtual works using either the effective intergranular stress, with the pore pressure gradient regarded as a distributed body force, or the total stress incorporating the pore pressure. In the finite element (FE) method both approaches prove equivalent at the global assembled level. However, at the element level apparently the equivalence does not hold, and the strength source related to the pore pressure seems to generate different local forces on the element nodes. The two formulations are briefly reviewed and discussed for triangular and tetrahedral finite elements. They are shown to yield different results at the global level as well in a three‐dimensional axisymmetric porous medium if the FE integration is performed using the average element‐wise radius. A modification to both formulations is suggested which allows to correctly solve the problem of a finite reservoir with an infinite pressure gradient, i.e. with a pore pressure discontinuity on its boundary. Copyright © 2001 John Wiley & Sons, Ltd.     

Pixel‐scale evaluation of SSM/I sea‐ice algorithms in the marginal ice zone during early fall freeze‐up

Observed reduction in recent sea ice areal extent and thickness has focused attention on the fact that the Arctic marine system appears to be responding to global‐scale climate variability and change. Passive microwave remote‐sensing data are the primary source underpinning these reports, yet problems remain in geophysical inversion of information on ice type and concentration. Uncertainty in sea‐ice concentration (SIC) retrievals is highest in the summer and fall, when water occurs in liquid phase within the snow–sea‐ice system. Of particular scientific interest is the timing and rate of new ice formation due to the control that this form of sea ice has on mass, energy and gas fluxes across the ocean–sea‐ice–atmosphere interface. In this paper we examine the critical fall freeze‐up period using in situ data from a ship‐based and aerial survey programme known as the Canadian Arctic Shelf Exchange study combined with microwave and optical Earth observations data. Results show that: (1) the overall physical conditions observed from aerial survey photography were well matched with coincident moderate‐resolution imaging spectroradiometer data and Radarsat ScanSAR imagery; (2) the shortwave albedo was linearly related to old ice concentration derived from survey photography; (3) the three SSM/I SIC algorithms (NASA Team (NT), NASA Team 2 (NT2), and Bootstrap (BT)) showed considerable discrepancies in pixel‐scale comparison with the Radarsat ScanSAR SICs well calibrated by the aerial survey data. The major causes of the discrepancies are attributed to (1) the inherent inability to detect the new thin ice in the NT and BT algorithms, (2) mismatches of the thin‐ice tie point of the NT2 algorithm, and (3) sub‐pixel ambiguity between the thin ice and the mixture of open water and sea ice. These results suggest the need for finer resolution of passive microwave sensors, such as AMSR‐E, to improve the precision of the SSM/I SIC algorithms in the marginal ice zone during early fall freeze‐up. Copyright © 2006 John Wiley & Sons, Ltd.     

Bulk organic δ13C and C/N ratios as palaeosalinity indicators within a Scottish isolation basin

Microfossils in isolation basin sediments are frequently used to reconstruct sea‐level change, but preservation problems and non‐analogue situations can limit their usefulness. Here we investigate the potential of stable carbon isotopes (δ¹³C) and C/N ratios from bulk organic matter, as an alternative proxy of salinity within isolation basin sediments from a basin in northwest Scotland. Within the Holocene sediment δ¹³C and C/N are determined largely by the mean weighted values of the predominant source of the organic material. Analysis of modern materials and comparison with the diatom record shows that the marine parts of the sequence are dominated by high δ¹³C and variable C/N. In the fresh water sequences the organic material is a mixture of both freshwater aquatic and terrestrial plant input that have relatively low δ¹³C and high C/N. The application of δ¹³C and C/N ratios in the studied basin in general follow the environmental change recorded by the diatoms and shows the potential of bulk organic matter in the investigation of salinity change in isolation basins. Copyright © 2005 John Wiley & Sons, Ltd.     

A non‐coaxial critical state soil model and its application to simple shear simulations

The yield vertex non‐coaxial theory is implemented into a critical state soil model, CASM (Int. J. Numer. Anal. Meth. Geomech. 1998; 22 :621–653) to investigate the non‐coaxial influences on the stress–strain simulations of real soil behaviour in the presence of principal stress rotations. The CASM is a unified clay and sand model, developed based on the soil critical state concept and the state parameter concept. Without loss of simplicity, it is capable of simulating the behaviour of sands and clays within a wide range of densities. The non‐coaxial CASM is employed to simulate the simple shear responses of Erksak sand and Weald clay under different densities and initial stress states. Dependence of the soil behaviour on the Lode angle and different plastic flow rules in the deviatoric plane are also considered in the study of non‐coaxial influences. All the predictions indicate that the use of the non‐coaxial model makes the orientations of the principal stress and the principal strain rate different during the early stage of shearing, and they approach the same ultimate values with an increase in loading. These ultimate orientations are dependent on the density of soils, and independent of their initial stress states. The use of the non‐coaxial model also softens the shear stress evolutions, compared with the coaxial model. It is also found that the ultimate shear strengths by using the coaxial and non‐coaxial models are dependent on the plastic flow rules in the deviatoric plane. Copyright © 2006 John Wiley & Sons, Ltd.     

Hydrometeorological behaviour of pine and larch forests in eastern Siberia

Seasonal changes in the water and energy exchanges over a pine forest in eastern Siberia were investigated and compared with published data from a nearby larch forest. Continuous observations (April to August 2000) were made of the eddy‐correlation sensible heat flux and latent heat flux above the canopy. The energy balance was almost closed, although the sum of the turbulent fluxes sometimes exceeded the available energy flux (R_n ? G) when the latent heat flux was large; this was related to the wind direction. We examined the seasonal variation in energy balance components at this site. The seasonal variation and magnitude of the sensible heat flux (H) was similar to that of the latent heat flux (λE), with maximum values occurring in mid‐June. Consequently, the Bowen ratio was around 1·0 on many days during the study period. On some clear days just after rainfall, λE was very large and the sum of H and λE exceeded R_n ? G. The evapotranspiration rate above the dry canopy from May to August was 2·2 mm day^?1. The contributions of understory evapotranspiration (E_u) and overstory transpiration (E_o) to the evapotranspiration of the entire ecosystem (E_t) were both from 25 to 50% throughout the period analysed. These results suggest that E_u plays a very important role in the water cycle at this site. From snowmelt through the tree growth season (23 April to 19 August 2000), the total incoming water, comprised of the sum of precipitation and the water equivalent of the snow at the beginning of the melt season, was 228 mm. Total evapotranspiration from the forest, including interception loss and evaporation from the soil when the canopy was wet, was 208–254 mm. The difference between the incoming and outgoing amounts in the water balance was from +20 to ?26 mm. The water and energy exchanges of the pine and larch forest differed in that λE and H increased slowly in the pine forest, whereas λE increased rapidly in the larch forest and H decreased sharply after the melting season. Consequently, the shape of the Bowen ratio curves at the two sites differed over the period analysed, as a result of the differences in the species in each forest and in soil thawing. Copyright © 2003 John Wiley & Sons, Ltd.     

Comparison of infiltration models

Fourteen popular, representative infiltration models, some physically based, some semi‐empirical and some empirical, were selected for a comparative evaluation. Using the Nash and Sutcliffe efficiency criterion, the models were evaluated and compared for 243 sets of infiltration data collected from field and laboratory tests conducted in India and the USA on soils ranging from coarse sand to fine clay. Based on a relative grading scale, the semi‐empirical Singh–Yu general model, Holtan model and Horton model were graded respectively as 6·52, 5·57 and 5·48 out of 10. The empirical Huggins and Monke model, modified Kostiakov and Kostiakov model were graded as 5·57, 5·30 and 5·22, respectively. The physically based non‐linear and linear models of Smith–Parlange were graded as 5·48 and 5·22, respectively. Other models were ranked lower than these models. All the models generally performed poorly in field tests on Georgia's sandy soils, except the Robertsdale loamy sand. Copyright © 2003 John Wiley & Sons, Ltd.     

Seasonal evolution of meltwater generation,storage and discharge at a non‐temperate glacier in Svalbard

In glacierized catchments, meteorological inputs driving surface melting are translated into runoff outputs mediated by the glacier hydrological system: analysis of the relationship between meteorology and diurnal and seasonal patterns of runoff should reflect the functioning of that system, with the role of meltwater storage likely to be of particular importance. Daily meltwater storage is determined for a glacier at 78 °N in the Svalbard archipelago, by comparing inputs calculated from a surface energy balance model with measured outputs (proglacial discharge). Solar radiation, air temperature, wind speed and proglacial discharge are then analysed by regression and time‐series methods, in order to assess the meteorology–discharge relationship and its variation at diurnal and seasonal time‐scales. The recorded discharge time‐series can be divided into two contrasting intervals: up to early August, proglacial discharge was high and variable, mean hydrographs showed little indication of diurnal cycling, ARIMA models of discharge indicated a non‐seasonal, moving‐average generating process, and there was a net loss of meltwater from storage; from early August, proglacial discharge was low and relatively invariable, but with clearer diurnal cycles, regression models of discharge showed substantially improved correlations with air temperature and solar radiation, ARIMA models indicated a non‐seasonal, autoregressive generating process, and eventually a seasonal component, and there was a net gain in meltwater storage. The transition between the two periods is brief compared with the duration of the melt season. The runoff response to meteorology therefore lacks the strongly progressive element previously identified in mid‐latitude glacierized catchments. In particular, the glacier hydrological system only appears responsive to diurnal forcing following the depletion of the seasonal snowpack meltwater store. Copyright © 2001 John Wiley & Sons, Ltd.