Retreating Laurentide Ice Sheet Research Articles

Imprint of relative sea level histories on Last Interglacial coral preservation

SUMMARY Fossil corals are commonly used to reconstruct Last Interglacial (∼125 ka, LIG) sea level. Sea level reconstructions assume the water depth at which the coral lived, called the ‘relative water depth’. However, relative water depth varies in time and space due to coral reef growth in response to relative sea level (RSL) changes. RSL changes can also erode coral reefs, exposing older reef surfaces with different relative water depths. We use a simplified numerical model of coral evolution to investigate how sea level history systematically influences the preservation of corals in the Bahamas and western Australia, regions which house >100 LIG coral fossils. We construct global ice histories spanning the uncertainty of LIG global mean sea level (GMSL) and predict RSL with a glacial isostatic adjustment model. We then simulate coral evolution since 132 ka. We show that preserved elevations and relative water depths of modelled LIG corals are sensitive to the magnitude, timing and number of GMSL highstand(s). In our simulations, the influence of coral growth and erosion (i.e. the ‘growth effect’) can have an impact on RSL reconstructions that is comparable to glacial isostatic adjustment. Thus, without explicitly accounting for the growth effect, additional uncertainty is introduced into sea level reconstructions. Our results suggest the growth effect is most pronounced in western Australia due to Holocene erosion, but also plays a role in the Bahamas, where LIG RSL rose rapidly due to the collapsing peripheral bulge associated with Laurentide Ice Sheet retreat. Despite the coral model's simplicity, our study highlights the utility of process-based RSL reconstructions.

Periglacial resurfacing of hillslopes and channels with large boulders in the Virginia Appalachians

AbstractLarge, resistant, quartz‐rich boulders deposited on hillslopes and in channels armour the landscape, trap sediment and influence hillslope angle and erodibility. In the Virginia Appalachians, such boulders are a significant component of hillslopes and channels. Establishing the timing of and processes responsible for bedrock fracture and boulder deposition is a critical piece of understanding the landscape as a system. In this study, we use cosmogenic 10Be exposure dating to resolve the timing of boulder deposition at three sites in the Virginia Valley and Ridge province: Gap Mountain, Brush Mountain and Little Stony Creek, and at one site in the Virginia Blue Ridge: Devil's Marbleyard. The correlation between measured boulder exposure ages (101.7 ± 6.9 ka to 10.8 ± 0.8 ka; n = 23) and the Wisconsin Glacial Stage and subsequent Laurentide Ice Sheet (LIS) deglaciation (~115–11.7 ka) suggests a periglacial origin for deposition of large hillslope and channel boulders in the Virginia Appalachians. The lack of boulder exposure ages corresponding to the Last Interglacial Stage or following Wisconsin LIS retreat suggests interglacial non‐deposition and stability. The absence of exposure ages from the penultimate Illinoian or older Quaternary Glacial Stages suggests that periglacial hillslope processes allow the landscape to be resurfaced with large boulders during each return to cold climate conditions. This cyclic resurfacing of hillslopes and channels is an example of how climatic oscillations insert disequilibrium into the landscape cycle and contributes to our appreciation of the timescales over which contemporary climate change may impact boulder dominated landscapes in rapidly warming alpine and arctic environments.

In Situ Cosmogenic 10Be Dating of Laurentide Ice Sheet Retreat from Central New England, USA

Constraining the timing and rate of Laurentide Ice Sheet (LIS) retreat through the northeastern United States is important for understanding the co-evolution of complex climatic and glaciologic events that characterized the end of the Pleistocene epoch. However, no in situ cosmogenic 10Be exposure age estimates for LIS retreat exist through large parts of Connecticut or Massachusetts. Due to the large disagreement between radiocarbon and 10Be ages constraining LIS retreat at the maximum southern margin and the paucity of data in central New England, the timing of LIS retreat through this region is uncertain. Here, we date LIS retreat through south-central New England using 14 new in situ cosmogenic 10Be exposure ages measured in samples collected from bedrock and boulders. Our results suggest ice retreated entirely from Connecticut by 18.3 ± 0.3 ka (n = 3). In Massachusetts, exposure ages from similar latitudes suggest ice may have occupied the Hudson River Valley up to 2 kyr longer (15.2 ± 0.3 ka, average, n = 2) than the Connecticut River Valley (17.4 ± 1.0 ka, average, n = 5). We use these new ages to provide insight about LIS retreat timing during the early deglacial period and to explore the mismatch between radiocarbon and cosmogenic deglacial age chronologies in this region.

The history of ice-sheet retreat on North America during Termination 5: Implications for the origin of the sea-level highstand during interglacial stage 11

Termination (T) 5, ∼424 ka, involved the biggest deglaciation of land-ice mass during the Quaternary. Warming and ice-sheet retreat during T5 led to an exceptionally long period of interglacial warmth known as Marine Isotope Stage (MIS) 11, ∼424–395 ka. A detailed understanding of the history of continental ice-sheet decay during T5 is required to disentangle regional contributions of ice-sheet retreat to sea-level rise (that range between ∼1 and 13 m above present day) and to correct it for glacio-isostatic adjustments (GIA). Yet little is known about the timing and magnitude of retreat during this time of the volumetrically most important continental ice sheet in the Northern Hemisphere, the Laurentide Ice Sheet (LIS). Here we present new authigenic Fe-Mn oxyhydroxide-derived high-resolution records of Pb isotope data and associated rare earth element profiles for samples spanning T5 from Labrador Sea IODP Site U1302/3. These records feature astronomically-paced radiogenic Pb isotope excursions that track increases in chemical weathering of North American bedrock and freshwater routing to the Labrador Sea via Hudson Straits associated with LIS retreat. Our records show that LIS retreat during T5 began 429. 2 ± 7.9 ka (2σ) and likely occurred over a longer timescale (by ∼10 to 5 kyr) than that observed for T2 and T1. They also show that Hudson Bay Ice Saddle collapse (and therefore LIS break-up) occurred ∼419 ± 4.7 ka (2σ), around the same time as best estimates of southern Greenland deglaciation, but ∼12 kyr before LIS deglaciation and the sea-level high-stand associated with the latter half of MIS 11 likely occurred. Our findings therefore highlight that ice-mass loss on North America likely played an important role in the seemingly protracted nature of T5 sea-level rise. A comparison of the deglaciation histories of the LIS and the southern Greenland Ice Sheet during T5, T2 and T1 also demonstrates that the well-constrained history of regional ice-sheet retreat during T1 is not always applicable as a template for older late Pleistocene terminations in GIA modelling.

Seismic reflection stratigraphy of Lac Simard, Quebec, Canada: mass flow sedimentation in glacial Lake Barlow-Ojibway

Lac Simard is a glaciated Shield basin situated in the Western Quebec Seismic Zone, an area of recurrent intraplate seismicity encompassing several densely populated areas. Previous works nearby have demonstrated that sediments preserved in lake basins hold the potential to serve as valuable archives of paleoseismicity. Correspondingly, a high-resolution chirp seismic reflection survey was conducted in Lac Simard, with results showing four acoustic facies (Af-1 to Af-4) within two stratigraphic successions (SS-1 and SS-2). The lowermost SS-1 is dominated by high-amplitude, laminated Af-2, interpreted as lateglacial varves deposited by suspension from primarily overflows and interflows during glacial lake ponding following deglaciation. Suspension deposition of Af-2 was episodically interrupted by stacked, chaotically bedded, weakly graded to opaque, channelled Af-3 interpreted as either debris flows or high-density turbidity currents. Presence of silt-clast breccias in Af-3 also suggests downslope collapse and reworking of varves from basin sidewalls. The overlying SS-2 is dominated by weakly reflective, laminated, high-frequency varved Af-4, suggesting a continuing seasonal control on postglacial sedimentation but with relative sediment starvation. A basin-wide erosional unconformity separating SS-1 and postglacial SS-2 records the abrupt drainage of glacial Lake Barlow-Ojibway around 8 000 years ago. Considering the history of lake floor disturbance in the region, mass flow facies Af-3 may be seismogenic, and thus, have regional tectonic significance. However, we cannot rule out the possibility of mass flows triggered by abrupt inflows of meltwater and sediment from a highly dynamic retreating Laurentide Ice Sheet prone to surging.

Deglacial release of petrogenic and permafrost carbon from the Canadian Arctic impacting the carbon cycle

The changes in atmospheric pCO2 provide evidence for the release of large amounts of ancient carbon during the last deglaciation. However, the sources and mechanisms that contributed to this process remain unresolved. Here, we present evidence for substantial ancient terrestrial carbon remobilization in the Canadian Arctic following the Laurentide Ice Sheet retreat. Glacial-retreat-induced physical erosion of bedrock has mobilized petrogenic carbon, as revealed by sedimentary records of radiocarbon dates and thermal maturity of organic carbon from the Canadian Beaufort Sea. Additionally, coastal erosion during the meltwater pulses 1a and 1b has remobilized pre-aged carbon from permafrost. Assuming extensive petrogenic organic carbon oxidation during the glacial retreat, a model-based assessment suggests that the combined processes have contributed 12 ppm to the deglacial CO2 rise. Our findings suggest potentially positive climate feedback of ice-sheet retreat by accelerating terrestrial organic carbon remobilization and subsequent oxidation during the glacial-interglacial transition.

Evidence of An Ice-Dammed Lake and Laurentide Readvance Upper Susquehanna Valley, New York State

Landforms and well logs document a system of ice-contact and proglacial lakes in the upper Susquehanna valley during Laurentide Ice Sheet retreat from the Appalachian Plateau, central New York State.  Recessional moraines formed dams for all lakes, except a newly revealed “Ancestral Goodyear Lake” retained behind an ephemeral ice dam stranded at Colliersville.  A prominent dead-ice sink currently occupies the valley floor at the dam site
 
 Ancestral Goodyear Lake held a stable lake level at 1360 feet as represented by thick lake sediments perched in water well logs on the valley wall above Goodyear Lake.  A deltaic terrace at 1250 feet in the same vicinity marks a second, lower lake strand.   
 
 In addition, water well logs on the adjacent Portlandville Moraine contain lake sediments bound above and below by ice-contact material deposits thus demonstrating a Laurentide readvance that subsequently dammed the valley to form Glacial Lake Milford as part of the Susquehanna Lake System.

Laurentide Ice Sheet extent over the last 130 thousand years traced by the Pb isotope signature of weathering inputs to the Labrador Sea

Understanding the history of continental ice-sheet growth on North America, and in particular that of the Laurentide Ice Sheet (LIS), is important for palaeoclimate and sea-level reconstructions. Information on ice-sheet extent pre-dating the Last Glacial Maximum (LGM) is heavily reliant, though, on the outputs of numerical models underpinned by scant geological data. Important aspects of LIS history that remain unresolved include the timing of its collapse during Termination 2, the first time that it expanded significantly during the Last Glacial Cycle, and whether its volume was significantly reduced during marine isotope stage (MIS) 3. To address these issues and more, we present authigenic iron-manganese (Fe–Mn) oxyhydroxide-derived high-resolution records of Pb isotope data and associated rare earth element profiles for samples spanning the past ∼130 kyr from northwest North Atlantic Labrador Sea, IODP Site U1302/3. We use these new data to track chemical weathering intensity and solute flux to the Labrador Sea associated with LIS extent on the adjacent highly radiogenic (high Pb isotope composition) North American Superior Province (SP) craton since the Penultimate Glacial Maximum (PGM). Our new records show that relatively high (radiogenic) values characterise warm marine isotope stages (MIS) 5, 3 and 1 and the lowest (most unradiogenic) values occurred during cold stages MIS 6, 4 and 2. The radiogenic Pb isotope excursion associated with Termination 2 is short-lived relative to the one documented for Termination 1, suggesting that LIS retreat during the PGM was relatively fast compared to the LGM and that its collapse during the last interglacial occurred ∼125 ka. Highly radiogenic inputs to the Labrador Sea during MIS 5d-a, ∼116–71 ka, most likely reflect a spin-up in Labrador Current vigour, incipient glaciation and renewed glacial erosion of high grounds of the eastern SP craton by localised wet-based ice-caps. A large decrease in Pb isotope values towards unradiogenic LGM-like compositions between ∼75–65 ka across the MIS 5/4 transition likely reflects a slow-down in Labrador Current vigour, an increase in subaerial deposition of aeolian dust and a significant advance of the LIS across Hudson Bay caused a strong reduction or even abandonment of Pb sourcing from the SP. The relatively radiogenic Pb isotope composition of bottom-waters bathing our study site during MIS 3, 57–29 ka, is unlikely to support a recently proposed major reduction in LIS extent for this time. Instead, we argue these values are better explained by southern Greenland Ice Sheet retreat, increased chemical weathering of the Ketelidian Mobile Belt and subsequent Pb runoff from Greenland.

Arctic drainage of Laurentide Ice Sheet meltwater throughout the past 14,700 years

During the last deglaciation substantial volumes of meltwater from the decaying Laurentide Ice Sheet were supplied to the Arctic, Gulf of Mexico and North Atlantic along different drainage routes, sometimes as catastrophic flood events. These events are suggested to have impacted global climate, for example initiating the Younger Dryas cold period. Here we analyze the authigenic Pb isotopic composition of sediments in front of the Arctic Mackenzie Delta, a sensitive tracer for elevated freshwater runoff of the retreating Laurentide Ice Sheet. Our data reveal continuous meltwater supply to the Arctic along the Mackenzie River since the onset of the Bølling–Allerød. The strongest Lake Agassiz outflow event is observed at the end of the Bølling–Allerød close to the onset of the Younger Dryas. In context of deglacial North American runoff records from the southern and eastern outlets, our findings provide a detailed reconstruction of the deglacial drainage chronology of the disintegrating Laurentide Ice Sheet.

Reconstruction of isostatically adjusted paleo-strandlines along the southern margin of the Laurentide Ice Sheet in the Great Lakes, Lake Agassiz, and Champlain Sea basins

Strandlines document the former presence of lakes and a sea in east-central North America along the southern margin of the retreating Laurentide Ice Sheet (LIS). The strandlines of these formerly level water bodies are uplifted to the north and provide evidence of glacial isostatic adjustment (GIA) of the Earth’s crust to the former ice load. We compile published ages and measurements of the present elevation and location of shore features in the strandlines of eight major paleo-water bodies from the St. Lawrence Valley to the northern Great Plains in digital format as an aid for the numerical modelling of GIA. Data for eastern water bodies were extracted and digitized from publications during the past 120 years. Digital position coordinates were scaled from published maps of survey sites or were determined using Google Earth Pro software. Published data for paleo-lakes Duluth and Agassiz were mainly obtained from field measurements and digital elevation models. Two-sigma or 95% probability values are provided for the strandline ages and for isobase (contour) positions representing the deformed water surfaces. Peak strandline gradients reported here were largest at about ca. 13 000 years ago. Lower strandline gradients for older shores may reveal areas closer to the peripheral bulge and areas of thinner ice (lighter crustal loads). Concave upward strandline profiles characterize most paleo-basins, whereas a linear uplift profile characterizes the Champlain Sea strandline. Directions of strandline maximum uplift within the former water body basins point toward the thickest part of the LIS near the Quebec–Labrador ice dome.

The potential of Hudson Valley glacial floods to drive abrupt climate change

The periodic input of meltwater into the ocean from a retreating Laurentide Ice Sheet is often hypothesized to have weakened the Atlantic meridional overturning circulation (AMOC) and triggered several cold periods during the last deglaciation (21,000 to 8,000 years before present). Here, we use a numerical model to investigate whether the Intra-Allerød Cold Period was triggered by the drainage of Glacial Lake Iroquois, ~13,300 years ago. Performing a large suite of experiments with various combinations of single and successive, short (1 month) and long (1 year) duration flood events, we were unable to find any significant weakening of the AMOC. This result suggests that although the Hudson Valley floods occurred close to the beginning of the Intra-Allerød Cold Period, they were unlikely the sole cause. Our results have implications for re-evaluating the relationship of meltwater flood events (past and future) to periods of climatic cooling, particularly with regards to flood input location, volume, frequency, and duration.

Long-term and recent ecohydrological dynamics of patterned peatlands in north-central Quebec (Canada)

Peatlands are natural ecosystems that provide archives of the hydrological cycle, ecological processes and terrestrial carbon dynamics. In the north-central region of Quebec (eastern Canada), patterned peatlands developed in topographic depressions of the Precambrian Shield following the Laurentide Ice Sheet retreat. These peatlands display characteristics similar to appa mires and other peatlands that developed at the ecotone between the open (taiga) and closed boreal forest biomes of the Northern Hemisphere, and also correspond to the biogeographic limit between ombrotrophic and minerotrophic peatlands. During the Neoglacial cooling period in northeastern Canada, patterned peatlands, mainly oligotrophic fens, registered a hydrological disequilibrium expressed by an increase in surface wetness as aquatic microforms expanded to the detriment of terrestrial surfaces. Ecohydrological trajectories were reconstructed from a detailed study of two patterned peatlands in order to document their sensitivity to climate variations. To do this, plant macrofossil and testate amoeba data were combined with peat carbon accumulation rates, C:N ratios,210Pb and14C chronologies. Data show that peatlands initiated ca 6500 cal. y BP as ombrotrophic or minerotrophic systems depending on site-specific conditions, followed by a general increase in surface wetness during the Neoglacial cooling until the end of the Little Ice Age. A relatively synchronous ecosystem state shift from oligotrophic to more ombrotrophic conditions was registered at the beginning of the 20th century in central and lateral cores of both study sites, evoking the likely influence of recent warming on peat accumulation. These results suggest a potential northward migration of the biogeographic limit of the ombrotrophic peatland distribution during the 20th century, which could have implications for the role of these ecosystems as C sinks at the continental scale. Overall, these peatlands have stored a mean carbon mass of ca 100 kg m−2.

Atmospherically produced beryllium-10 in annually laminated late-glacial sediments of the North American Varve Chronology

Abstract. We attempt to synchronize the North American Varve Chronology (NAVC) with ice core and calendar year timescales by comparing records of atmospherically produced 10Be fallout in the NAVC and in ice cores. The North American Varve Chronology (NAVC) is a sequence of 5659 varves deposited in a series of proglacial lakes adjacent to the southeast margin of the retreating Laurentide Ice Sheet between approximately 18 200 and 12 500 years before present. Because properties of NAVC varves are related to climate, the NAVC is also a climate proxy record with annual resolution, and our overall goal is to place the NAVC and ice core records on the same timescale to facilitate high-resolution correlation of climate proxy variations in both. Total 10Be concentrations in NAVC sediments are within the range of those observed in other lacustrine records of 10Be fallout, but 9Be and 10Be concentrations considered together show that the majority of 10Be is present in glacial sediment when it enters the lake, and only a minority of total 10Be derives from atmospheric fallout at the time of sediment deposition. Because of this, an initial experiment to determine whether or not 10Be fallout variations were recorded in NAVC sediments by attempting to observe the characteristic 11-year solar cycle in short varve sections sampled at high resolution was inconclusive: short-period variations at the expected magnitude of this cycle were not distinguishable from measurement scatter. On the other hand, longer varve sequences sampled at decadal resolution display centennial-period variations in reconstructed 10Be fallout that have similar properties as coeval 10Be fallout variations recorded in ice core records. These are most prominent in glacial sections of the NAVC that were deposited in proglacial lakes and are suppressed in paraglacial sections of the NAVC that were deposited in lakes lacking direct glacial sediment input. We attribute this difference to the fact that buffering of 10Be fallout by soil adsorption can filter out short-period variations in an entirely deglaciated watershed, but such buffering cannot occur in the ablation zone of an ice sheet. This implies that proglacial lakes whose watershed is mostly glacial may effectively record 10Be fallout variations. We attempted to match centennial-period variations in reconstructed 10Be fallout flux from two segments of the NAVC with ice core fallout records. For both records, it is possible to obtain matches that result in acceptable correlation between NAVC and ice core 10Be fallout records, but the best-fitting matches for the two segments disagree, and only one of them is consistent with independent calendar year calibrations of the NAVC and therefore potentially valid. This leaves several remaining ambiguities in whether or not 10Be fallout variations can, in fact, be used for synchronizing NAVC and ice core timescales, but these could most likely be resolved by higher-resolution and replicate 10Be measurements on targeted sections of the NAVC.

Refining the glacial lake coverage of the southern Laurentide ice margin using Lidar-DEM based reconstructions: The case of Lake Obedjiwan in south-central Quebec, Canada

The existence of former glacial lakes remains poorly documented in the remote regions of the Canadian Shield where the dense forest cover and poor access complicate classical mapping methods. Here we use a high-resolution Digital Elevation Model (DEM; 1 m/pixel) based on LiDAR data to map Quaternary deposits and landforms in the vicinity of the Gouin reservoir in the Laurentian Highlands of Quebec (Canada). A total of 598 raised shorelines were identified, providing new evidence for the presence of a large proglacial lake – Lake Obedjiwan – that formed in front of the northward retreating Laurentide ice sheet margin during the last deglaciation, a few kilometers east of the larger Lake Ojibway that covered extensive areas in northern Ontario and Quebec. The elevation of shorelines outlines a tilted water plane that depicts the influence of the glacio-isostatic adjustment across the basin, with shorelines ranging from 409 m in the south to 439 m in the north. The calculated uplift gradient of 0.78 m/km indicates that the lake developed during the early stages of the deglaciation. Our reconstruction constrains the full areal extent and maximum elevation of the lake, which covered an area of 1800 km2 and reached an elevation of 450 m. It also indicates that Lake Obedjiwan developed independently from Lake Ojibway. The lake-surface elevation was controlled by an outlet located in the southeastern part of the basin where a topographic depression lying at an elevation of 390 m routed meltwater overflow southward through the Saint-Maurice River and into the St. Lawrence River. Another important outlet is located in the northernmost part of the basin and ice retreat beyond this topographic depression (427 m) caused the drawdown of the Lake Obedjiwan, which drained westward into Lake Ojibway. This discharge added around 19 km3 of meltwater, an event that likely impacted the hydrological budget and sedimentology of Lake Ojibway. The results suggest that other glacial lakes probably developed along the southern LIS ice margin and that the approach developed in this study may contribute to the refinement of their areal coverage and improve paleogeographic reconstructions.

Sedimentary and rock magnetic signatures and event scenarios of deglacial outburst floods from the Laurentian Channel Ice Stream

Eastern Canadian margin sediments bear testimony to several catastrophic deglacial meltwater discharges from the retreating Laurentide Ice Sheet. The reddish-brown plumite layers deposited on the levees of the Laurentian Fan valleys have been recognized as indications of multiple outburst floods between Heinrich events 2 and 1. Five event layers have been consistently recorded in three new gravity cores retrieved on the SW Grand Banks slope and comply with the previously published Laurentian Fan core MD95-2029. The apparently huge extent of these outburst plumes around the Laurentian Fan as well as their causes and consequences are investigated in this study using physical properties, rock magnetic and grain-size analyses, together with seismoacoustic profiling. We provide the first detailed 14C ages of the outburst event sequence and discuss their recurrence intervals in the context of regional ice retreat. Compared to the hemipelagic interlayers, event layers have overall uniform and systematic changes of rock-magnetic properties. Hematite contents increase over time and proximally while magnetite grain sizes fine upwards and spatially away from the fan. Based on the sediment composition and load, we argue that these plumites were formed by recurrent erosion of glacial mud deposits in the Laurentian Channel by meltwater outbursts. Three alternative glaciological scenarios are evaluated: in each case, the provenance of the transported sediment is not an indicator of the precise source of the meltwater.

Character and origin of De Geer moraines in the Seacoast region of New Hampshire, USA

ABSTRACTPreviously unrecognized clusters of regularly spaced, low‐relief ridges were recently revealed by LiDAR data acquired in the Seacoast region of New Hampshire, USA. Based on their morphometry and geologic context, these features are interpreted as De Geer moraines that formed during regional deglaciation. While De Geer moraines are widely recognized as products of ice recession in glaciomarine settings, their formative processes remain debated. Here, we integrate geospatial mapping, morphometric analyses, stratigraphic investigations and ground‐penetrating radar (GPR) subsurface profiling to characterize De Geer moraines in New Hampshire and assess their origin. On average, the moraines are 1.2 km long, 28 m wide, 1.7 m tall and spaced 104 m apart. Trenches expose stratigraphic units consistent with subaqueous deposition along a tidewater glacier margin, including subglacial traction till indicative of proglacial pushing, and gravity‐driven flow deposits on an ice‐distal moraine slope. GPR profiles reveal variable till thicknesses beneath the ridges and stratified sediment in adjacent swales. Our findings suggest De Geer moraine formation was driven by annual oscillations of the retreating Laurentide Ice Sheet margin. In this interpretation, average moraine spacing would represent a retreat rate of ∼100 m a−1, which is in line with previous estimates of recession rates across New England.

Spatial Autocorrelation, Source Water and the Distribution of Total and Viable Microbial Abundances within a Crystalline Formation to a Depth of 800 m.

Proposed radioactive waste repositories require long residence times within deep geological settings for which we have little knowledge of local or regional subsurface dynamics that could affect the transport of hazardous species over the period of radioactive decay. Given the role of microbial processes on element speciation and transport, knowledge and understanding of local microbial ecology within geological formations being considered as host formations can aid predictions for long term safety. In this relatively unexplored environment, sampling opportunities are few and opportunistic. We combined the data collected for geochemistry and microbial abundances from multiple sampling opportunities from within a proposed host formation and performed multivariate mixing and mass balance (M3) modeling, spatial analysis and generalized linear modeling to address whether recharge can explain how subsurface communities assemble within fracture water obtained from multiple saturated fractures accessed by boreholes drilled into the crystalline formation underlying the Chalk River Laboratories site (Deep River, ON, Canada). We found that three possible source waters, each of meteoric origin, explained 97% of the samples, these are: modern recharge, recharge from the period of the Laurentide ice sheet retreat (ca. ∼12000 years before present) and a putative saline source assigned as Champlain Sea (also ca. 12000 years before present). The distributed microbial abundances and geochemistry provide a conceptual model of two distinct regions within the subsurface associated with bicarbonate – used as a proxy for modern recharge – and manganese; these regions occur at depths relevant to a proposed repository within the formation. At the scale of sampling, the associated spatial autocorrelation means that abundances linked with geochemistry were not unambiguously discerned, although fine scale Moran’s eigenvector map (MEM) coefficients were correlated with the abundance data and suggest the action of localized processes possibly associated with the manganese and sulfate content of the fracture water.

Fracturing‐induced release of radiogenic 4He and 234U into groundwater during the last deglaciation: An alternative source to crustal helium fluxes in periglacial aquifers

AbstractExternal 4He sources have been invoked to explain 4He concentrations in groundwater greater than those expected from in situ U and Th production. In a fractured aquifer of Ordovician age located in the St. Lawrence Lowlands (Quebec, Canada), 4He concentrations of up to 4.48 × 10−5 cm3 STP were measured. Such concentrations are ∼1000 times higher than would be expected from in situ production. A concomitant increase in 4He concentration and 234U/238U activity ratio is shown, suggesting a common release process in groundwater for 234U and 4He. This process has tentatively been identified as glaciation‐induced rock fracturing following the Laurentide Ice Sheet retreat. The resulting increase in exposed grain surface facilitates 234U release by α‐recoil and that of radiogenic 4He by diffusion. Using a model of helium diffusion from a spherical grain, it is shown that rock fracturing facilitated the release of accumulated 4He at rates ranging from 4.2 × 10−10 to 1.06 × 10−8 cm3 STP yr−1. These release rates are between 1000 and 30,000 times higher than the local U and Th steady state production rate, of 3.5 ± 1.4 × 10−13 cm3 STP yr−1. Integration of 4He release rates over time yields a radiogenic 4He concentration of between 3.85 × 10−6 and 7.12 × 10−5 cm3 STP , in the range of concentrations measured in the St. Lawrence Lowlands fractured aquifers. Results support the occurrence of a local radiogenic helium source to explain the excesses measured in groundwater without requiring a significant external He crustal flux.

Cosmogenic nuclide age estimate for Laurentide Ice Sheet recession from the terminal moraine, New Jersey, USA, and constraints on latest Pleistocene ice sheet history

AbstractThe time at which the Laurentide Ice Sheet reached its maximum extent and subsequently retreated from its terminal moraine in New Jersey has been constrained by bracketing radiocarbon ages on preglacial and postglacial sediments. Here, we present measurements of in situ produced10Be and26Al in 16 quartz-bearing samples collected from bedrock outcrops and glacial erratics just north of the terminal moraine in north-central New Jersey; as such, our ages represent a minimum limit on the timing of ice recession from the moraine. The data set includes field and laboratory replicates, as well as replication of the entire data set five years after initial measurement. We find that recession of the Laurentide Ice Sheet from the terminal moraine in New Jersey began before 25.2±2.1 ka (10Be, n=16, average, 1 standard deviation). This cosmogenic nuclide exposure age is consistent with existing limiting radiocarbon ages in the study area and cosmogenic nuclide exposure ages from the terminal moraine on Martha’s Vineyard ~300 km to the northeast. The age we propose for Laurentide Ice Sheet retreat from the New Jersey terminal position is broadly consistent with regional and global climate records of the last glacial maximum termination and records of fluvial incision.

Rapid thinning of the Laurentide Ice Sheet in coastal Maine, USA, during late Heinrich Stadial 1

Few data are available to infer the thinning rate of the Laurentide Ice Sheet (LIS) through the last deglaciation, despite its importance for constraining past ice sheet response to climate warming. We measured 31 cosmogenic 10Be exposure ages in samples collected on coastal mountainsides in Acadia National Park and from the slightly inland Pineo Ridge moraine complex, a ∼100-km-long glaciomarine delta, to constrain the timing and rate of LIS thinning and subsequent retreat in coastal Maine. Samples collected along vertical transects in Acadia National Park have indistinguishable exposure ages over a 300 m range of elevation, suggesting that rapid, century-scale thinning occurred at 15.2 ± 0.7 ka, similar to the timing of abrupt thinning inferred from cosmogenic exposure ages at Mt. Katahdin in central Maine (Davis et al., 2015). This rapid ice sheet surface lowering, which likely occurred during the latter part of the cold Heinrich Stadial 1 event (19–14.6 ka), may have been due to enhanced ice-shelf melt and calving in the Gulf of Maine, perhaps related to regional oceanic warming associated with a weakened Atlantic Meridional Overturning Circulation at this time. The ice margin subsequently stabilized at the Pineo Ridge moraine complex until 14.5 ± 0.7 ka, near the onset of Bølling Interstadial warming. Our 10Be ages are substantially younger than marine radiocarbon constraints on LIS retreat in the coastal lowlands, suggesting that the deglacial marine reservoir effect in this area was ∼1,200 14C years, perhaps also related to the sluggish Atlantic Meridional Overturning Circulation during Heinrich Stadial 1.