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.     

The evolution of the terrestrial-terminating Irish Sea glacier during the last glaciation

Here we reconstruct the last advance to maximum limits and retreat of the Irish Sea Glacier (ISG), the only land-terminating ice lobe of the western British Irish Ice Sheet. A series of reverse bedrock slopes rendered proglacial lakes endemic, forming time-transgressive moraine- and bedrock-dammed basins that evolved with ice marginal retreat. Combining, for the first time on glacial sediments, optically stimulated luminescence (OSL) bleaching profiles for cobbles with single grain and small aliquot OSL measurements on sands, has produced a coherent chronology from these heterogeneously bleached samples. This chronology constrains what is globally an early build-up of ice during late Marine Isotope Stage 3 and Greenland Stadial (GS) 5, with ice margins reaching south Lancashire by 30 ± 1.2 ka, followed by a 120-km advance at 28.3 ± 1.4 ka reaching its 26.5 ± 1.1 ka maximum extent during GS-3. Early retreat during GS-3 reflects piracy of ice sources shared with the Irish-Sea Ice Stream (ISIS), starving the ISG. With ISG retreat, an opportunistic readvance of Welsh ice during GS-2 rode over the ISG moraines occupying the space vacated, with ice margins oscillating within a substantial glacial over-deepening. Our geomorphological chronosequence shows a glacial system forced by climate but mediated by piracy of ice sources shared with the ISIS, changing flow regimes and fronting environments.     

Landscape and prehistoric chronology of West-Central Jordan

Remnant lake and stream terraces of the Wadi el Hasa (west-central Jordan) are associated with in situ prehistoric sites spanning > 100,000 years. Eighteen radiocarbon dates from cultural and geological deposits on the terraces facilitate the first comprehensive prehistoric landscape chronology for the southern Levant east of the Jordan Rift. In the eastern Hasa basin, the uppermost of three cut and fill surfaces (>20 m) is linked to massive fossil spring deposits and an early Middle Paleolithic occupation (100,000–70,000 B. P.), suggestive of considerably wetter climates. A later Middle Paleolithic occupation may be synchronous with the emergence of Pleistocene Lake Hasa (ca. 70,000 B. P.). Peak lake levels were attained 40,000 years ago. Dates proliferate after 25,000 B. P. and register recession of Lake Hasa (ca. 20,000 B. P.), an intervening erosional phase, and the initiation of complex humid-desiccation cycles for the terminal Pleistocene—Holocene (17,000–9,000 B. P.). The contemporary Wadi el Hasa channel began aggrading its floodplain after 8000 B. P. and was incised to its present depths 1000–500 years ago. The prehistoric landscape history of the Hasa drainage is broadly synchronous with sequences in the Rift Valley and Negev desert and offers baseline chronologies for the Late Quaternary of eastern Jordan and the Arabian peninsula. © 1994 John Wiley & Sons, Inc.     

Interactions between submarine channel systems and deformation in deepwater fold belts: Examples from the Levant Basin, Eastern Mediterranean sea

Submarine channel levee systems form important hydrocarbon reservoirs in many deep marine settings and are often deposited within a structurally active setting. This study focuses on recent submarine channels that developed within a deepwater fold and thrust belt setting from the Levant Basin, eastern Mediterranean Sea. Compressional deformation within the study area is driven by the up-dip collapse of the Nile cone above the ductile Messinian Evaporites. Structures such as folds and strike slip faults exert a strong control on channel location and development over time. From this study four end-member submarine channel–structure interactions can be defined: Confinement, diversion, deflection and blocking. Each of these channel–structure interactions results in a distinct submarine channel morphology and pattern of development compared to unconfined channel levee systems. Each interaction can also be used to assess timing relationships between submarine channel development and deformation.     

Revised recommended methods for analyzing crater size-frequency distributions

Impact crater populations help us to understand solar system dynamics, planetary surface histories, and surface modification processes. A single previous effort to standardize how crater data are displayed in graphs, tables, and archives was in a 1978 NASA report by the Crater Analysis Techniques Working Group, published in 1979 in Icarus. The report had a significant lasting effect, but later decades brought major advances in statistical and computer sciences while the crater field has remained fairly stagnant. In this new work, we revisit the fundamental techniques for displaying and analyzing crater population data and demonstrate better statistical methods that can be used. Specifically, we address (1) how crater size-frequency distributions (SFDs) are constructed, (2) how error bars are assigned to SFDs, and (3) how SFDs are fit to power-laws and other models. We show how the new methods yield results similar to those of previous techniques in that the SFDs have familiar shapes but better account for multiple sources of uncertainty. We also recommend graphic, display, and archiving methods that reflect computers’ capabilities and fulfill NASA's current requirements for Data Management Plans.     

Eocene tectono-thermal rejuvenation of an upper Paleozoic-lower Mesozoic terrane in the Cordillera de Carabaya, Puno, southeastern Peru, revealed by K---Ar and Ar/Ar dating

K---Ar dates for muscovites and biotites in granitoid rocks and hydrothermal ore deposits of the northeastern parts of the plutons making up the Triassic Carabaya batholith, underlying the axial Cordillera Oriental of northern Puno Department, southeastern Peru, are markedly variable and mutually discordant. Steep transverse gradients are defined in the apparent ages of both micas, which decrease systematically from SW to NE, delimiting a ca. 25-km-wide, longitudinal zone of anomalously young Mesozoic to Paleocene dates. Age minima of 37±1 Ma are attained in three of the four studied transects. ⁴⁰Ar/³⁹Ar step-heating analyses of selected micas confirm the occurrence of a thermal disturbance, and modeling of the spectra suggests that argon loss in muscovites attains at least ca. 75% in the northeastern part of the zone. A single K-feldspar spectrum yielded a minimum at 31 Ma, and apatite fission-track age cluster at ca. 31 and 18.5 Ma. The affected granitoid rocks generally display little megascopic evidence of tectonism, but microscopic deformational fabrics increase in intensity with apparent decreasing K---Ar age, paralleling a marked increase in alkali feldspar ordering. Secondary fluid inclusions trapped within the microfabrics reveal that the plutonic rocks were penetrated by a homogeneous H₂O---CO₂---CH₄---NaCl fluid at ca. 300–400°C and 0.7–2 kbar. This fluid is implicated in the degassing of the rocks. These diverse data are interpreted as evidence for a major, but moderate-temperature (400°C) and brief, tectono-thermal event at ca. 37±1 Ma (biotite closure temperature)—i.e., at the Eocene-Oligocene boundary. The K-feldspar ⁴⁰Ar/³⁹Ar data and the Oligocene fission-track dates may record the later stages in the event, whereas the Miocene fission-track dates are tentatively ascribed to a distinct Neogene episode. Essentially identical geochronological and petrological relationships have been documented in the Cordillera Real of northwestern Bolivia by McBride et al. (1987), permitting the delimination of a disturbed belt paralleling the South American plate boundary and more than 450 km long. The tectono-thermal domain, which we term the Zongo-San Gabán Zone, constituted the foreland boundary of the Andean orogen in the vicinity of the Arica Deflection during the late Eocene Incaic orogeny. This regional thermal event, which involved the basement, appears to have resulted from compressional or, in some segments, transpressional tectonics.     

Time scale of hydrothermal water-rock reactions in Yellowstone National Park based on radium isotopes and radon

We have measured ²²⁴Ra (3.4 d), ²²⁸Ra (5.7 yr), and ²²⁶Ra (1620 yr) and chloride in hot spring waters from the Norris-Mammoth Corridor, Yellowstone National Park. Two characteristic cold-water components mix with the primary hydrothermal water: one for the travertine-depositing waters related to the Mammoth Hot Springs and the other for the sinter-depositing Norris Geyser Basin springs. The Mammoth Hot Springs water is a mixture of the primary hydrothermal fluid with meteoric waters flowing through the Madison Limestone, as shown by the systematic decrease of the (²²⁸Ra/²²⁶Ra) activity ratio proceeding northward. The Norris Geyser Basin springs are mixtures of primary hydrothermal water with different amounts of cold meteoric water with no modification of the primary hydrothermal (²²⁸Ra/²²⁶Ra) activity ratio. Using a solution and recoil model for radium isotope supply to the primary hydrothermal water, a mean water-rock reaction time prior to expansion at 350°C and supply to the surface is 540 years assuming that 250 g of water are involved in the release of the radium from one gram of rock. The maximum reaction time allowed by our model is 1150 years.     

Lead isotopic evidence for evolutionary changes in magma-crust interaction,Central Andes,southern Peru

Lead isotopic measurements were made on Andean igneous rocks of Jurassic to Recent age in Moquegua and Tacna Departments, southernmost Peru, to clarify the petrogenesis of the rocks and, in particular, to investigate the effect of crustal thickness on rock composition. This location in the Cordillera Occidental is ideal for such a study because the ca. 2 Ga Precambrian basement rocks (Arequipa massif) have a distinct Pb isotopic signature which is an excellent tracer of crustal interaction, and because geomorphological research has shown that the continental crust was here thickened drastically in the later Tertiary.Seven samples of quartz diorites and granodiorites from the Ilo and Toquepala intrusive complexes, and seven samples of Toquepala Group subaerial volcanics were analyzed for Pb isotopic compositions. The plutonic rocks range in age from Jurassic to Eocene; the volcanic rocks are all Late Cretaceous to Eocene. With one exception, the Pb isotopic ratios are in the ranges ²⁰⁶Pb/²⁰⁴Pb= 18.52–18.75, ²⁰⁷Pb/²⁰⁴Pb= 15.58–15.65, and ²⁰⁸Pb/²⁰⁴Pb= 38.53–38.74. The data reflect very little or no interaction with old continental material of the Arequipa massif type.Lead from four Miocene Huaylillas Formation ash-flow tuffs, two Pliocene Capillune Formation andesites and five Quaternary Barroso Group andesites has lower ²⁰⁶Pb/²⁰⁴Pb than that in the pre-Miocene rocks, but relatively high ²⁰⁷Pb/²⁰⁴Pb and ²⁰⁸Pb/²⁰⁴Pb (²⁰⁶Pb/²⁰⁴Pb= 18.16–18.30, ²⁰⁷Pb/²⁰⁴Pb= 15.55–15.63, ²⁰⁸Pb/²⁰⁴Pb= 38.45–38.90). Tilton and Barreiro [9] have shown that contamination by Arequipa massif granulites can explain the isotopic composition of the Barosso Group lavas, and the new data demonstrate that this effect is evident, to varying degrees, in all the analysed Neogene volcanic rocks. The initial incorporation of such basement material into the magma coincided with the Early Miocene uplift of this segment of the Cordillera Occidental [32], and thus with the creation of a thick crustal root. The data strongly imply a relationship between crustal thickness and degree of crustal contamination of magmas in this area, but a rigorous relationship is not yet established.