Search

Founding populations of the first Americans likely occupied parts of Beringia during the Last Glacial Maximum (LGM). The timing, pathways, and modes of their southward transit remain unknown, but blockage of the interior route by North American ice sheets between ~26 and 14 cal kyr BP (ka) favors a coastal route during this period. Using models and paleoceanographic data from the North Pacific, we identify climatically favorable intervals when humans could have plausibly traversed the Cordilleran coastal corridor during the terminal Pleistocene. Model simulations suggest that northward coastal currents strengthened during the LGM and at times of enhanced freshwater input, making southward transit by boat more difficult. Repeated Cordilleran glacial-calving events would have further challenged coastal transit on land and at sea. Following these events, ice-free coastal areas opened and seasonal sea ice was present along the Alaskan margin until at least 15 ka. Given evidence for humans south of the ice sheets by 16 ka and possibly earlier, we posit that early people may have taken advantage of winter sea ice that connected islands and coastal refugia. Marine ice-edge habitats offer a rich food supply and traversing coastal sea ice could have mitigated the difficulty of traveling southward in watercraft or on land over glaciers. We identify 24.5 to 22 ka and 16.4 to 14.8 ka as environmentally favorable time periods for coastal migration, when climate conditions provided both winter sea ice and ice-free summer conditions that facilitated year-round marine resource diversity and multiple modes of mobility along the North Pacific coast.

During the last deglaciation, dozens of glacial outburst floods-among the largest known floods on Earth-scoured the Channeled Scabland landscape of eastern Washington. Over this same period, deformation of the Earth's crust in response to the growth and decay of ice sheets changed the topography by hundreds of meters. Here, we investigated whether glacial isostatic adjustment affected routing of the Missoula floods and incision of the Channeled Scabland from an impounded, glacial Lake Columbia. We used modern topography corrected for glacial isostatic adjustment as an input to flood models that solved the depth-averaged, shallow water equations and compared the results to erosion constraints. Results showed that floods could have traversed and eroded parts of two major tracts of the Channeled Scabland-Telford-Crab Creek and Cheney-Palouse-near 18 ka, whereas glacial isostatic adjustment limited flow into the Cheney-Palouse tract at 15.5 ka. Partitioning of flow between tracts was governed by tilting of the landscape, which affected the filling and overspill of glacial Lake Columbia directly upstream of the tracts. These results highlight the impact of glacial isostatic adjustment on megaflood routing and landscape evolution., Competing Interests: The authors declare no competing interest., (Copyright © 2022 the Author(s). Published by PNAS.)

An understanding of the mechanisms that control CO2 change during glacial-interglacial cycles remains elusive. Here we help to constrain changing sources with a high-precision, high-resolution deglacial record of the stable isotopic composition of carbon in CO2(δ(13)C-CO2) in air extracted from ice samples from Taylor Glacier, Antarctica. During the initial rise in atmospheric CO2 from 17.6 to 15.5 ka, these data demarcate a decrease in δ(13)C-CO2, likely due to a weakened oceanic biological pump. From 15.5 to 11.5 ka, the continued atmospheric CO2 rise of 40 ppm is associated with small changes in δ(13)C-CO2, consistent with a nearly equal contribution from a further weakening of the biological pump and rising ocean temperature. These two trends, related to marine sources, are punctuated at 16.3 and 12.9 ka with abrupt, century-scale perturbations in δ(13)C-CO2 that suggest rapid oxidation of organic land carbon or enhanced air-sea gas exchange in the Southern Ocean. Additional century-scale increases in atmospheric CO2 coincident with increases in atmospheric CH4 and Northern Hemisphere temperature at the onset of the Bølling (14.6-14.3 ka) and Holocene (11.6-11.4 ka) intervals are associated with small changes in δ(13)C-CO2, suggesting a combination of sources that included rising surface ocean temperature.

Erosion, sediment production, and routing on a tectonically active continental margin reflect both tectonic and climatic processes; partitioning the relative importance of these processes remains controversial. Gulf of Alaska contains a preserved sedimentary record of the Yakutat Terrane collision with North America. Because tectonic convergence in the coastal St. Elias orogen has been roughly constant for 6 My, variations in its eroded sediments preserved in the offshore Surveyor Fan constrain a budget of tectonic material influx, erosion, and sediment output. Seismically imaged sediment volumes calibrated with chronologies derived from Integrated Ocean Drilling Program boreholes show that erosion accelerated in response to Northern Hemisphere glacial intensification (∼ 2.7 Ma) and that the 900-km-long Surveyor Channel inception appears to correlate with this event. However, tectonic influx exceeded integrated sediment efflux over the interval 2.8-1.2 Ma. Volumetric erosion accelerated following the onset of quasi-periodic (∼ 100-ky) glacial cycles in the mid-Pleistocene climate transition (1.2-0.7 Ma). Since then, erosion and transport of material out of the orogen has outpaced tectonic influx by 50-80%. Such a rapid net mass loss explains apparent increases in exhumation rates inferred onshore from exposure dates and mapped out-of-sequence fault patterns. The 1.2-My mass budget imbalance must relax back toward equilibrium in balance with tectonic influx over the timescale of orogenic wedge response (millions of years). The St. Elias Range provides a key example of how active orogenic systems respond to transient mass fluxes, and of the possible influence of climate-driven erosive processes that diverge from equilibrium on the million-year scale.

Understanding responses of oceanic primary productivity, carbon export, and burial to climate change is essential for model-based projection of biological feedbacks in a high-CO2 world. Here we compare estimates of productivity based on the composition of fossil diatom floras with organic carbon burial off Oregon in the Northeast Pacific across a large climatic transition at the last glacial termination. Although estimated primary productivity was highest during the Last Glacial Maximum, carbon burial was lowest, reflecting reduced preservation linked to low sedimentation rates. A diatom size index further points to a glacial decrease (and deglacial increase) in the fraction of fixed carbon that was exported, inferred to reflect expansion, and contraction, of subpolar ecosystems that today favor smaller plankton. Thus, in contrast to models that link remineralization of carbon to temperature, in the Northeast Pacific, we find dominant ecosystem and sea floor control such that intervals of warming climate had more efficient carbon export and higher carbon burial despite falling primary productivity.

Deciphering the evolution of global climate from the end of the Last Glacial Maximum approximately 19 ka to the early Holocene 11 ka presents an outstanding opportunity for understanding the transient response of Earth's climate system to external and internal forcings. During this interval of global warming, the decay of ice sheets caused global mean sea level to rise by approximately 80 m; terrestrial and marine ecosystems experienced large disturbances and range shifts; perturbations to the carbon cycle resulted in a net release of the greenhouse gases CO(2) and CH(4) to the atmosphere; and changes in atmosphere and ocean circulation affected the global distribution and fluxes of water and heat. Here we summarize a major effort by the paleoclimate research community to characterize these changes through the development of well-dated, high-resolution records of the deep and intermediate ocean as well as surface climate. Our synthesis indicates that the superposition of two modes explains much of the variability in regional and global climate during the last deglaciation, with a strong association between the first mode and variations in greenhouse gases, and between the second mode and variations in the Atlantic meridional overturning circulation.