skip to main content


Title: Landscape evolution under the southern Laurentide Ice Sheet

Ancient glaciated landscapes reveal interactions among ice dynamics, bed composition, and preglacial river networks.

Subglacial landscapes, revealed in regions of recent ice-sheet retreat, provide a window into ice-sheet dynamics and interactions with evolving subglacial topography. Here, we document landscape evolution beneath the southern Laurentide Ice Sheet of North America since the end of the Pliocene, 2.6 million years (Ma) ago, by reconstructing the isostatically adjusted preglacial surface and modern bedrock topography at 250 m horizontal resolution. We use flow routing to reconstruct drainage networks and river longitudinal profiles, revealing the pattern and extent of their glacially forced reorganization. The overall mean Quaternary (2.6 Ma ago to present) erosion rate is 27 m/Ma, rising within ice-streaming corridors to 35 m/Ma (and locally reaching 400 m/Ma) and falling to 22 m/Ma in non–ice-streaming regions. Our results suggest that subglacial erosion was sufficient to lower the southern Laurentide Ice Sheet into warmer environments, thereby enhancing ablation and reducing ice-sheet extent over time. 
more » « less
Award ID(s):
1903606
NSF-PAR ID:
10471349
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
AAAS
Date Published:
Journal Name:
Science Advances
Volume:
7
Issue:
48
ISSN:
2375-2548
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Using offshore detrital apatite (U‐Th)/He thermochronometry and 3D thermo‐kinematic modeling of the catchment topography, we constrain the timing of major topographic change at Bourgeois Fjord, Antarctic Peninsula (AP). While many mid‐latitude glacial landscapes developed primarily in response to global cooling over the last ~2.6 Ma, we find that kilometer‐scale landscape evolution at Bourgeois Fjord began ~30–12 Ma ago and <2 km of valley incision has occurred since ~16 Ma. This early onset of major topographic change occurred following the initiation of alpine glaciation at this location and prior to the development of a regional polythermal ice sheet inferred from sedimentary evidence offshore of the AP. We hypothesize that topographic change relates to (i) feedbacks between an evolving topography and glacial erosion processes, (ii) effects of glacial‐interglacial variability, and (iii) the prevalence of subglacial meltwater. The timing and inferred spatial patterns of long‐term exhumation at Bourgeois Fjord are consistent with a hypothesis that glacial erosion processes were suppressed at the AP during global Plio‐Pleistocene cooling, rather than enhanced. Our study examines the long‐term consequences of glacial processes on catchment‐wide erosion as the local climate cooled. Our findings support the hypothesis that landscapes at different latitudes had different responses to global cooling. Our results also suggest that erosion is enhanced along the plateau flanks of Bourgeois Fjord today, which may be due to periglacial processes or mantling via subglacial till. If regional warming persists and meltwater becomes more pronounced, we predict that enhanced erosion along the plateau flank will accelerate topographic change.

     
    more » « less
  2. Richard B. Waitt ; Glenn D. Thackray ; Alan R. Gillespie (Ed.)
    The northward retreat history of the Laurentide ice sheet through the lowlands of the northeastern United States during the last deglaciation is well constrained, but its vertical thinning history is less well known because of the lack of direct constraints on ice thickness through time and space. In addition, the highest elevations in New England are characterized by gently sloping upland surfaces and weathered block fields, features with an uncertain history. To better constrain ice-sheet history in this area and its relationship to alpine geomorphology, we present 20 new 10Be and seven in situ 14C cosmogenic nuclide measurements along an elevation transect at Mount Washington, New Hampshire, the highest mountain in the northeastern United States (1917 m above sea level [a.s.l.]). Our results suggest substantially different exposure and erosion histories on the upper and lower parts of the mountain. Above 1600 m a.s.l., 10Be and in situ 14C measurements are consistent with upper reaches of the mountain deglaciating by 18 ka. However, some 10Be ages are up to several times greater than the age of the last deglaciation, consistent with weakly erosive, cold-based ice that did not deeply erode preglacial surfaces. Below 1600 m a.s.l., 10Be ages are indistinguishable over a nearly 900 m range in elevation and imply rapid ice-surface lowering ca. 14.1 ± 1.1 ka (1 standard deviation; n = 9). This shift from slow thinning early in the deglaciation on the upper part of the mountain to abrupt thinning across the lower elevations coincided with accelerated ice-margin retreat through the region recorded by Connecticut River valley varve records during the Bølling interstadial. The Mount Washington cosmogenic nuclide vertical transect and the Connecticut River valley varve record, along with other New England cosmogenic nuclide records, suggest rapid ice-volume loss in the interior northeastern United States in response to Bølling warming. 
    more » « less
  3. Bedrock topography is a key boundary condition for ice sheet modeling, and determining changes in subglacial topography through time can provide insight into the timing of ice sheet development, the magnitude of glacial erosion, and the co-development of glaciers and glacial topography. West Antarctica hosts an unusually high geothermal gradient supported by hot, low-viscosity mantle which likely enhanced the lithospheric response to West Antarctic Ice Sheet (WAIS) cycles of growth and increased the sensitivity of thermochronometers to landscape evolution on million-year timescales. Thus, a valuable record of glacial landscape change might be recovered from apatite fission track [AFT 80-130°C range] and (U-Th)/He [AHe; 50-90°C] dating, provided that landscape evolution can be distinguished from tectonic signals, including the effects of faults. This study utilizes AFT-AHe thermochronology and thermo-kinematic Pecube modeling to investigate interactions between the hot geotherm, glacial erosion, and inferred crustal structures in the Ford Ranges and the DeVicq Glacier trough in western and central Marie Byrd Land (MBL), respectively. The Ford Ranges host glacial troughs (up to 3km relief) dissecting a low-relief erosional surface. Previous work suggests a majority of bedrock exhumation and cooling occurred at/by 80 Ma. However, new data hint at renewed exhumation linked to glacial incision since WAIS formation at 34 or 20 Ma. Prior (U-Th)/He zircon dates from exposures of crystalline bedrock span 90 – 67 Ma. New AHe bedrock dates are 41 to 26 Ma, while two glacial erratics (presumed to be eroded from walls or floor of glacial troughs) yielded AHe dates of 37 Ma and 16 Ma. Initial modeling results suggest a tectonic boundary between the Ford Ranges and Edward VII Peninsula separating regions with distinct exhumation histories. The boundary may cause differential WAIS incision at 34 or 20 Ma, a possibility being investigated with new models. The DeVicq Glacier trough (>3.5km relief) coincides with a prominent crustal lineament but lacks temperature-time information compared to other regions. The crustal structure may have accommodated motion between elevated central MBL and the subdued crust of the Ford Ranges. Here, owing to the lack of onshore non-volcanic bedrock exposure, we have employed AHe and AFT dating of glacial sediment marine core samples offshore of the DeVicq Glacier to investigate the timing and rates of exhumation of the bedrock carved by the DeVicq trough, with initial results revealing detrital AHe ages as young as 24 Ma. Our new Pecube models test a series of thermal, tectonic, and landscape evolution scenarios against a suite of thermochronologic data, allowing us to assess the timing of glacial incision and WAIS initiation in the Ford Ranges, and to seek evidence of an inferred tectonic boundary at DeVicq Trough. Modeling efforts will be aided by new AHe and AFT analyses from ongoing work. These models combine topographic, tectonic, thermal, and key thermochronologic datasets to produce new insight into the unique cryosphere-lithosphere interactions affecting landscape change in West Antarctica. 
    more » « less
  4. null (Ed.)
    The Indian (southwest) summer monsoon is one of the most intense climatic phenomena on Earth. Its long-term development has been linked to the growth of high topography in South and Central Asia. The Indian continental margin, adjoining the Arabian Sea, offers a unique opportunity to investigate tectonic–climatic interactions and the net impact of these processes on weathering and erosion of the western Himalaya. During International Ocean Discovery Program Expedition 355, two sites (U1456 and U1457) were drilled in Laxmi Basin in the eastern Arabian Sea to document the coevolution of mountain building, weathering, erosion, and climate over a range of timescales. In addition, recovering basement from the eastern Arabian Sea provides constraints on the early rifting history of the western continental margin of India with special emphasis on continental breakup between India and the Seychelles and its relationship to the plume-related volcanism of the Deccan Plateau. Drilling and coring operations during Expedition 355 recovered sediment from Sites U1456 and U1457 in the Laxmi Basin, penetrating 1109.4 and 1108.6 m below seafloor (mbsf), respectively. Drilling reached sediment dated to 13.5–17.7 Ma (late early to early middle Miocene) at Site U1456, although with a large hiatus between the lowermost sediment and overlying deposits dated to <10.9 Ma. At Site U1457, a much longer hiatus occurs near the base of the cored section, spanning from 10.9 to ~62 Ma. At both sites, hiatuses span ~8.2–9.2 and ~3.6–5.6 Ma, with a possible condensed section spanning ~2.0–2.6 Ma, although the total duration for each hiatus is slightly different between the two sites. A major submarine fan draining the western Himalaya and Karakoram must have been supplying sediment to the eastern Arabian Sea since at least ~17 Ma. Sand mineral assemblages indicate that the Greater Himalayan Crystalline Sequence was fully exposed to the surface by this time. Most of the recovered sediment appears to be derived from the Indus River and includes minerals that are unique to the Indus Suture Zone, in particular glaucophane and hypersthene, most likely originating from the structural base of the Kohistan arc. Pliocene sandy intervals at Site U1456 were deposited in lower fan “sheet lobe” settings, with intervals of basin plain turbidites separated by hemipelagic muddy sections deposited during the Miocene. Site U1457 is more distal in facies, reflecting its more marginal setting. No major active lobe appears to have affected the Laxmi Basin since the Middle Pleistocene (~1.2 Ma). We succeeded in recovering sections spanning the 8 Ma climatic transition, when monsoon intensity is believed to have changed strongly, although the nature of this change awaits postcruise analysis. We also recovered sediment from a large mass transport deposit measuring ~330 and ~190 m thick at Sites U1456 and U1457, respectively. This section includes an upper sequence of slump-folded muddy and silty rocks, as well as underlying calcarenites and limestone breccias, together with smaller amounts of volcanic clasts, all of which are likely derived from the western Indian continental shelf. Identification of similar facies on the regional seismic lines in Laxmi Basin suggests that these deposits form parts of one of the world’s largest mass transport deposits. Coring of igneous basement was successful at Site U1457. Recovery of massive basalt and associated volcaniclastic sediment at this site should address the key questions related to rifting and volcanism associated with formation of Laxmi Basin. Geochemical analysis is required to understand the petrogenesis and thus the tectonic setting of volcanism that will reveal whether it is oceanic basalt or volcanic rock contaminated by underlying continental crust or continental flood basalt. However, the fact that the lavas are massive and have few vesicles implies water depths of eruption likely deeper than 2000 m. This precludes opening of the basin in the presence of a major mantle thermal anomaly, such as that associated with the Deccan Large Igneous Province. Other observations made at the two sites during Expedition 355 provide vital constraints on the rift history of this margin. Heat flow measurements at the two drill sites were calculated to be ~57 and ~60 mW/m2. Such heat flow values are compatible with those observed in average oceanic crust of 63–84 Ma age, as well as with the presence of highly extended continental crust. Postcruise analyses of the more than ~1722 m of core will provide further information about the nature of tectonic–climatic interactions in this global type area for such studies. 
    more » « less
  5. Abstract

    The Pliocene Epoch (∼5.3–2.6 million years ago, Ma) was characterized by a warmer than present climate with smaller Northern Hemisphere ice sheets, and offers an example of a climate system in long‐term equilibrium with current or predicted near‐future atmospheric CO2concentrations (pCO2). A long‐term trend of ice‐sheet expansion led to more pronounced glacial (cold) stages by the end of the Pliocene (∼2.6 Ma), known as the “intensification of Northern Hemisphere Glaciation” (iNHG). We assessed the spatial and temporal variability of ocean temperatures and ice‐volume indicators through the late Pliocene and early Pleistocene (from 3.3 to 2.4 Ma) to determine the character of this climate transition. We identified asynchronous shifts in long‐term means and the pacing and amplitude of shorter‐term climate variability, between regions and between climate proxies. Early changes in Antarctic glaciation and Southern Hemisphere ocean properties occurred even during the mid‐Piacenzian warm period (∼3.264–3.025 Ma) which has been used as an analog for future warming. Increased climate variability subsequently developed alongside signatures of larger Northern Hemisphere ice sheets (iNHG). Yet, some regions of the ocean felt no impact of iNHG, particularly in lower latitudes. Our analysis has demonstrated the complex, non‐uniform and globally asynchronous nature of climate changes associated with the iNHG. Shifting ocean gateways and ocean circulation changes may have pre‐conditioned the later evolution of ice sheets with falling atmosphericpCO2. Further development of high‐resolution, multi‐proxy reconstructions of climate is required so that the full potential of the rich and detailed geological records can be realized.

     
    more » « less