Pattern,style and timing of British–Irish Ice Sheet advance and retreat over the last 45 000 years: evidence from NW Scotland and the adjacent continental shelf

Predicting the future response of ice sheets to climate warming and rising global sea level is important but difficult. This is especially so when fast-flowing glaciers or ice streams, buffered by ice shelves, are grounded on beds below sea level. What happens when these ice shelves are removed? And how do the ice stream and the surrounding ice sheet respond to the abruptly altered boundary conditions? To address these questions and others we present new geological, geomorphological, geophysical and geochronological data from the ice-stream-dominated NW sector of the last British–Irish Ice Sheet (BIIS). The study area covers around 45 000 km² of NW Scotland and the surrounding continental shelf. Alongside seabed geomorphological mapping and Quaternary sediment analysis, we use a suite of over 100 new absolute ages (including cosmogenic-nuclide exposure ages, optically stimulated luminescence ages and radiocarbon dates) collected from onshore and offshore, to build a sector-wide ice-sheet reconstruction combining all available evidence with Bayesian chronosequence modelling. Using this information we present a detailed assessment of ice-sheet advance/retreat history, and the glaciological connections between different areas of the NW BIIS sector, at different times during the last glacial cycle. The results show a highly dynamic, partly marine, partly terrestrial, ice-sheet sector undergoing large size variations in response to sub-millennial-scale climatic (Dansgaard–Oeschger) cycles over the last 45 000 years. Superimposed on these trends we identify internally driven instabilities, operating at higher frequency, conditioned by local topographic factors, tidewater dynamics and glaciological feedbacks during deglaciation. Specifically, our new evidence indicates extensive marine-terminating ice-sheet glaciation of the NW BIIS sector during Greenland Stadials 12 to 9 – prior to the main ‘Late Weichselian’ ice-sheet glaciation. After a period of restricted glaciation, in Greenland Interstadials 8 to 6, we find good evidence for rapid renewed ice-sheet build-up in NW Scotland, with the Minch ice-stream terminus reaching the continental shelf edge in Greenland Stadial 5, perhaps only briefly. Deglaciation of the NW sector took place in numerous stages. Several grounding-zone wedges and moraines on the mid- and inner continental shelf attest to significant stabilizations of the ice-sheet grounding line, or ice margin, during overall retreat in Greenland Stadials 3 and 2, and to the development of ice shelves. NW Lewis was the first substantial present-day land area to deglaciate, in the first half of Greenland Stadial 3 at a time of globally reduced sea-level c. 26 kabp , followed by Cape Wrath at c. 24 kabp. The topographic confinement of the Minch straits probably promoted ice-shelf development in early Greenland Stadial 2, providing the ice stream with additional support and buffering it somewhat from external drivers. However, c. 20–19 kabp , as the grounding-line migrated into shoreward deepening water, coinciding with a marked change in marine geology and bed strength, the ice stream became unstable. We find that, once underway, grounding-line retreat proceeded in an uninterrupted fashion with the rapid loss of fronting ice shelves – first in the west, then the east troughs – before eventual glacier stabilization at fjord mouths in NW Scotland by ~17 kabp. Around the same time, ~19–17 kabp , ice-sheet lobes readvanced into the East Minch – possibly a glaciological response to the marine-instability-triggered loss of adjacent ice stream (and/or ice shelf) support in the Minch trough. An independent ice cap on Lewis also experienced margin oscillations during mid-Greenland Stadial 2, with an ice-accumulation centre in West Lewis existing into the latter part of Heinrich Stadial 1. Final ice-sheet deglaciation of NW mainland Scotland was punctuated by at least one other coherent readvance at c. 15.5 kabp , before significant ice-mass losses thereafter. At the glacial termination, c. 14.5 kabp , glaciers fed outwash sediment to now-abandoned coastal deltas in NW mainland Scotland around the time of global Meltwater Pulse 1A. Overall, this work on the BIIS NW sector reconstructs a highly dynamic ice-sheet oscillating in extent and volume for much of the last 45 000 years. Periods of expansive ice-sheet glaciation dominated by ice-streaming were interspersed with periods of much more restricted ice-cap or tidewater/fjordic glaciation. Finally, this work indicates that the role of ice streams in ice-sheet evolution is complex but mechanistically important throughout the lifetime of an ice sheet – with ice streams contributing to the regulation of ice-sheet health but also to the acceleration of ice-sheet demise via marine ice-sheet instabilities.     

Changes in timber haul emissions in the context of shifting forest management and infrastructure

Background  

Although significant amounts of carbon may be stored in harvested wood products, the extraction of that carbon from the forest generally entails combustion of fossil fuels. The transport of timber from the forest to primary milling facilities may in particular create emissions that reduce the net sequestration value of product carbon storage. However, attempts to quantify the effects of transport on the net effects of forest management typically use relatively sparse survey data to determine transportation emission factors. We developed an approach for systematically determining transport emissions using: 1) -remotely sensed maps to estimate the spatial distribution of harvests, and 2) - industry data to determine landscape-level harvest volumes as well as the location and processing totals of individual mills. These data support spatial network analysis that can produce estimates of fossil carbon released in timber transport.     

How accurately may we project tropical forest-cover change? A validation of a forward-looking baseline for REDD

The Reduced Emissions from Deforestation and forest Degradation (REDD+) mechanism of a future post-2012 global climate-change treaty would aim to give incentive to tropical countries to reduce deforestation and thus forest-carbon emissions. It would do so by crediting tropical countries for reducing deforestation relative to a baseline scenario describing carbon emissions and removals from forest-cover change expected in the absence of REDD+. Defining a credible and accurate baseline is both critical and challenging. One approach considered promising is spatial modelling to project forest-cover change on the basis of historical trends; yet few such projections have been validated at a national scale. We develop and validate a novel GEOMOD projection of forest-cover change in Panama over 2000–2008, based on trends over 1990–2000 and 25 drivers of forest-cover change. Compared with the actual landscape of 2008, our projection is 85.2% accurate at a 100-m pixel resolution. More error is attributable to the location of projected forest (8.6%) than to its area (6.2%). Accuracy was least where forest regeneration predominated (80%), and greatest where deforestation predominated (90%). Despite the sophistication of our projection, it is slightly less accurate than if we had assumed no forest-cover change over 2000–2008. We identify factors limiting projection accuracy, including the complexity of forest-cover change, the spatial variability of forest-carbon density, and the relatively small area of change at the national scale. We conclude that, with the exception of contexts where forest-cover change is significant and straightforward and where forest-carbon density relatively uniform (e.g., agricultural frontiers), spatially projected baselines are of limited value for REDD+ – their accuracy is too limited given their relative lack of transparency. Simpler, relatively coarse scale, retrospective baselines are recommended instead.     

Constitutive approach for rate‐sensitive anisotropic structured clays

This paper describes a constitutive approach to model the behavior of rate‐dependent anisotropic structured clay. Rate‐sensitivity is modeled using overstress viscoplasticity. Clay structure is treated as a viscous phenomenon whereby the viscosity of the undisturbed structured clay is initially very high and the viscosity degrades or decreases with plastic straining until the intrinsic or residual viscosity is reached. A microstructure tensor approach is used to make the structured viscosity anisotropic; whereas, the intrinsic viscosity is assumed to be isotropic. The behavior of the constitutive model is compared with the measured response of two clays (Gloucester and St. Vallier clay) from Eastern Canada during triaxial compression tests on specimens trimmed at different orientations to the vertical. The comparisons show that the constitutive framework is able to describe the anisotropic and rate‐sensitive response of both clays. The response of the model is also examined for the more general case of anisotropic consolidated triaxial compression and extension. Copyright © 2010 John Wiley & Sons, Ltd.     

Silurian and Devonian systems

The Ediacara mineral field is situated 30 km W of Beltana on the western margins of the Flinders Ranges, South Australia, and consists of silver‐lead and copper deposits in lower Cambrian carbonate rocks that contain anomalous base‐metal contents throughout the Adelaide Geosyncline. The lower Cambrian rocks, which consist of the basal Parachilna Formation and overlying Ajax Limestone, rest disconformably on the Precambrian, and at Ediacara occupy a shallow N‐S elongate syncline near the hinge zone of the Adelaide Geosyncline. The main primary ore minerals of the silver‐lead mineralization are galena and pyrite, with very minor chalcopyrite and sphalerite, and rare tetrahedrite and pearceite. The gangue consists mainly of silica (both chalcedony and quartz), with minor dolomite and rare barite. The mineralization is stratabound and occurs in conformable zones, the lowest of which commences about 30–50 m above the base of the Cambrian sequence. The host to the silver‐lead mineralization, the Ajax Limestone, can be subdivided into three units which represent a set of lithologies, structures and organic traces indicative of a shallow near‐shore carbonate environment. The silver‐lead mineralization is mainly present in sandy and laminated dolomites which were deposited in an environment ranging from sub‐tidal to bar and channel and tidal flat, respectively. Four types of mineralization have been recognized; disseminated sulphides of syngenetic and/or diagenetic origin and epigenetic concentrations along stylolites, in veins and as breccia fillings. Post‐depositional solution activity has affected a large proportion of the carbonate sequence. The effects of this activity range from stylolites through stylobreccias to solution collapse breccias. The epigenetic concentrations of mineralizations have apparently been formed by the remobilization of the disseminated sulphides during solution activity. The ore and gangue minerals of the epigenetic mineralization display both euhedral forms and distinct colloform banding, and framboidal textures have also been observed in both pyrite and galena. There is evidence of repeated episodic precipitation and no simple paragenetic sequence can be recognized. Fluid inclusions in silica and dolomite associated with the epigenetic mineralization have homogenization temperatures of 159 to 199°C and freezing temperatures that indicate the fluids to be saline brines containing NaCl with CaCl₂ and/or MgCl₂. Sulphur isotope analyses show a range of 8³⁴S values from ‐12.5 to +8.6 per mil, with no evidence of significant differences between the four types of mineralization. The data suggest deposition of the disseminated sulphides as a result of biological reduction of seawater sulphate in a system partially open with respect to sulphate supply. Subsequent remobilization of sulphides apparently involved little or no sulphur isotope fractionation. The Ediacara silver‐lead deposits have many features in common with Mississippi Valley‐type lead‐zinc deposits and appear to have similarities in terms of genesis, in that the epigenetic mineralization has been formed as a result of post‐depositional solution activity during diagenesis in a sedimentary basin. The scale of transport of the metals deposited as the epigenetic mineralization at Ediacara appears, however, to have been very much less than that of the metals in other Mississippi Valley‐type deposits.     

Soil hydrophobic effects on infiltration and catchment runoff

After dry summers or drought, eucalypt forest soils at two sites in southeastern Australia developed hydrophobic or non-wetting surface characteristics that reduced infiltration, measured using a sprinkling infiltrometer. At one site the development of hydrophobic conditions caused the rainfall to runoff conversion efficiency of a forested catchment to increase from 5 per cent to 15 per cent. Under non-hydrophobic conditions at this site, grassland always generated more runoff than forest. However, one major rainfall-runoff was recorded at a time of highly hydrophobic forest soil conditions and this storm generated greater runoff on the forested catchment than the grassland catchment. At the second site forest soils have naturally highly conductive surface layers because of a dense network of macropores and pathways for preferential flow. Hydrophobic conditions produced by drought caused soil water movement to be confined to only a few of the larger macropores exposed to surface ponded water. Even so, infiltration rates remained relatively high so that the impacts of hydrophobic soils were not translated into increased catchment runoff as at the first site.     

Complexity in astroplasmas

Ideas of complexity theory – including self-organized criticality – are being applied to the study of turbulence and astroplasmas in the solar system. An RAS MIST meeting in London brought together researchers from different disciplines to compare established paradigms with more recent ideas. Sandra Chapman , Mervyn Freeman and Sean Oughton report.     

Distorted Froude‐scaled flume analysis of large woody debris

This paper presents the results of a movable‐boundary, distorted, Froude‐scaled hydraulic model based on Abiaca Creek, a sand‐bedded channel in northern Mississippi. The model was used to examine the geomorphic and hydraulic impact of simplified large woody debris (LWD) elements. The theory of physical scale models is discussed and the method used to construct the LWD test channel is developed. The channel model had bed and banks moulded from 0·8 mm sand, and flow conditions were just below the threshold of motion so that any sediment transport and channel adjustment were the result of the debris element. Dimensions and positions of LWD elements were determined using a debris jam classification model. Elements were attached to a dynamometer to measure element drag forces, and channel adjustment was determined through detailed topographic surveys. The fluid drag force on the elements decreased asymptotically over time as the channel boundary eroded around the elements due to locally increased boundary shear stress. Total time for geomorphic adjustment computed for the prototype channel at the Q₂ discharge (discharge occurring once every two years on average) was as short as 45 hours. The size, depth and position of scour holes, bank erosion and bars created by flow acceleration past the elements were found to be related to element length and position within the channel cross‐section. Morphologies created by each debris element in the model channel were comparable with similar jams observed in the prototype channel. Published in 2001 John Wiley & Sons, Ltd.