Ancient bears provide insights into Pleistocene ice age refugia in Southeast Alaska

During the Late Pleistocene, major parts of North America were periodically covered by ice sheets. However, there are still questions about whether ice-free refugia were present in the Alexander Archipelago along the Southeast (SE) Alaska coast during the last glacial maximum (LGM). Numerous subfossils have been recovered from caves in SE Alaska, including American black (Ursus americanus) and brown (U. arctos) bears, which today are found in the Alexander Archipelago but are genetically distinct from mainland bear populations. Hence, these bear species offer an ideal system to investigate long-term occupation, potential refugial survival and lineage turnover. Here, we present genetic analyses based on 99 new complete mitochondrial genomes from ancient and modern brown and black bears spanning the last ~45,000 years. Black bears form two SE Alaskan subclades, one preglacial and another postglacial, that diverged >100,000 years ago. All postglacial ancient brown bears are closely related to modern brown bears in the archipelago, while a single preglacial brown bear is found in a distantly related clade. A hiatus in the bear subfossil record around the LGM and the deep split of their pre- and postglacial subclades fail to support a hypothesis of continuous occupancy in SE Alaska throughout the LGM for either species. Our results are consistent with an absence of refugia along the SE Alaska coast, but indicate that vegetation quickly expanded after deglaciation, allowing bears to recolonize the area after a short-lived LGM peak.

Late-Glacial and Holocene Lake-Level Fluctuations on the Kenai Lowland, Reconstructed from Satellite-Fen Peat Deposits and Ice-Shoved Ramparts, Kenai Peninsula, Alaska

Recent decades of warmer climate have brought drying wetlands and falling lake levels to southern Alaska. These recent changes can be placed into a longer-term context of postglacial lake-level fluctuations that include low stands that were as much as 7 m lower than present at eight lakes on the Kenai Lowland. Closed-basin lakes on the Kenai Lowland are typically ringed with old shorelines, usually as wave-cut scarps, cut several meters above modern lake levels; the scarps formed during deglaciation at 25–19 ka in a kettle moraine topography on the western Kenai Lowland. These high-water stands were followed by millennia of low stands, when closed-basin lake levels were drawn down by 5–10 m or more. Peat cores from satellite fens near or adjoining the eight closed-basin lakes show that a regional lake level rise was underway by at least 13.4 ka. At Jigsaw Lake, a detailed study of 23 pairs of overlapping sediment cores, seismic profiling, macrofossil analysis, and 58 AMS radiocarbon dates reveal rapidly rising water levels at 9–8 ka that caused large slabs of peat to slough off and sink to the lake bottom. These slabs preserve an archive of vegetation that had accumulated on a lakeshore apron exposed during the preceding drawdown period. They also preserve evidence of a brief period of lake level rise at 4.7–4.5 ka. We examined plant succession using in situ peat sequences in nine satellite fens around Jigsaw Lake that indicated increased effective moisture between 4.6 and 2.5 ka synchronous with the lake level rise. Mid- to late-Holocene lake high stands in this area are recorded by numerous ice-shoved ramparts (ISRs) along the shores. ISRs at 15 lakes show that individual ramparts typically record several shove events, separated by hundreds or thousands of years. Most ISRs date to within the last 5200 years and it is likely that older ISRs were erased by rising lake levels during the mid- to late Holocene. This study illustrates how data on vegetation changes in hydrologically coupled satellite-fen peat records can be used to constrain the water level histories in larger adjacent lakes. We suggest that this method could be more widely utilized for paleo-lake level reconstruction.

Steering iceberg armadas

The Asian-Pacific tropics likely instigated millennial-scale climate changes

High sensitivity of Bering Sea winter sea ice to winter insolation and carbon dioxide over the last 5500 years

Anomalously low winter sea ice extent and early retreat in CE 2018 and 2019 challenge previous notions that winter sea ice in the Bering Sea has been stable over the instrumental record, although long-term records remain limited. Here, we use a record of peat cellulose oxygen isotopes from St. Matthew Island along with isotope-enabled general circulation model (IsoGSM) simulations to generate a 5500-year record of Bering Sea winter sea ice extent. Results show that over the last 5500 years, sea ice in the Bering Sea decreased in response to increasing winter insolation and atmospheric CO₂, suggesting that the North Pacific is highly sensitive to small changes in radiative forcing. We find that CE 2018 sea ice conditions were the lowest of the last 5500 years, and results suggest that sea ice loss may lag changes in CO₂ concentrations by several decades.

The role of Northeast Pacific meltwater events in deglacial climate change

Columbia River megafloods occurred repeatedly during the last deglaciation, but the impacts of this fresh water on Pacific hydrography are largely unknown. To reconstruct changes in ocean circulation during this period, we used a numerical model to simulate the flow trajectory of Columbia River megafloods and compiled records of sea surface temperature, paleo-salinity, and deep-water radiocarbon from marine sediment cores in the Northeast Pacific. The North Pacific sea surface cooled and freshened during the early deglacial (19.0-16.5 ka) and Younger Dryas (12.9-11.7 ka) intervals, coincident with the appearance of subsurface water masses depleted in radiocarbon relative to the sea surface. We infer that Pacific meltwater fluxes contributed to net Northern Hemisphere cooling prior to North Atlantic Heinrich Events, and again during the Younger Dryas stadial. Abrupt warming in the Northeast Pacific similarly contributed to hemispheric warming during the Bølling and Holocene transitions. These findings underscore the importance of changes in North Pacific freshwater fluxes and circulation in deglacial climate events.

North Pacific freshwater events linked to changes in glacial ocean circulation

There is compelling evidence that episodic deposition of large volumes of freshwater into the oceans strongly influenced global ocean circulation and climate variability during glacial periods^1,2. In the North Atlantic region, episodes of massive freshwater discharge to the North Atlantic Ocean were related to distinct cold periods known as Heinrich Stadials^1-3. By contrast, the freshwater history of the North Pacific region remains unclear, giving rise to persistent debates about the existence and possible magnitude of climate links between the North Pacific and North Atlantic oceans during Heinrich Stadials^4,5. Here we find that there was a strong connection between changes in North Atlantic circulation during Heinrich Stadials and injections of freshwater from the North American Cordilleran Ice Sheet to the northeastern North Pacific. Our record of diatom δ¹⁸O (a measure of the ratio of the stable oxygen isotopes ¹⁸O and ¹⁶O) over the past 50,000 years shows a decrease in surface seawater δ¹⁸O of two to three per thousand, corresponding to a decline in salinity of roughly two to four practical salinity units. This coincided with enhanced deposition of ice-rafted debris and a slight cooling of the sea surface in the northeastern North Pacific during Heinrich Stadials 1 and 4, but not during Heinrich Stadial 3. Furthermore, results from our isotope-enabled model⁶ suggest that warming of the eastern Equatorial Pacific during Heinrich Stadials was crucial for transmitting the North Atlantic signal to the northeastern North Pacific, where the associated subsurface warming resulted in a discernible freshwater discharge from the Cordilleran Ice Sheet during Heinrich Stadials 1 and 4. However, enhanced background cooling across the northern high latitudes during Heinrich Stadial 3-the coldest period in the past 50,000 years⁷-prevented subsurface warming of the northeastern North Pacific and thus increased freshwater discharge from the Cordilleran Ice Sheet. In combination, our results show that nonlinear ocean-atmosphere background interactions played a complex role in the dynamics linking the freshwater discharge responses of the North Atlantic and North Pacific during glacial periods.

Deglaciation of the Pacific coastal corridor directly preceded the human colonization of the Americas

The route and timing of early human migration to the Americas have been a contentious topic for decades. Recent paleogenetic analyses suggest that the initial colonization from Beringia took place as early as 16 thousand years (ka) ago via a deglaciated corridor along the North Pacific coast. However, the feasibility of such a migration depends on the extent of the western Cordilleran Ice Sheet (CIS) and the available resources along the hypothesized coastal route during this timeframe. We date the culmination of maximum CIS conditions in southeastern Alaska, a potential bottleneck region for human migration, to ~20 to 17 ka ago with cosmogenic ¹⁰Be exposure dating and ¹⁴C dating of bones from an ice-overrun cave. We also show that productive marine and terrestrial ecosystems were established almost immediately following deglaciation. We conclude that CIS retreat ensured that an open and ecologically viable pathway through southeastern Alaska was available after 17 ka ago, which may have been traversed by early humans as they colonized the Americas.

Cordilleran Ice Sheet mass loss preceded climate reversals near the Pleistocene Termination

The Cordilleran Ice Sheet (CIS) once covered an area comparable to that of Greenland. Previous geologic evidence and numerical models indicate that the ice sheet covered much of westernmost Canada as late as 12.5 thousand years ago (ka). New data indicate that substantial areas throughout westernmost Canada were ice free prior to 12.5 ka and some as early as 14.0 ka, with implications for climate dynamics and the timing of meltwater discharge to the Pacific and Arctic oceans. Early Bølling-Allerød warmth halved the mass of the CIS in as little as 500 years, causing 2.5 to 3.0 meters of sea-level rise. Dozens of cirque and valley glaciers, along with the southern margin of the CIS, advanced into recently deglaciated regions during the Bølling-Allerød and Younger Dryas.

Synchronous deglacial thermocline and deep-water ventilation in the eastern equatorial Pacific

The deep ocean is most likely the primary source of the radiocarbon-depleted CO₂ released to the atmosphere during the last deglaciation. While there are well-documented millennial scale Δ¹⁴C changes during the most recent deglaciation, most marine records lack the resolution needed to identify more rapid ventilation events. Furthermore, potential age model problems with marine Δ¹⁴C records may obscure our understanding of the phase relationship between inter-ocean ventilation changes. Here we reconstruct changes in deep water and thermocline radiocarbon content over the last deglaciation in the eastern equatorial Pacific (EEP) using benthic and planktonic foraminiferal ¹⁴C. Our records demonstrate that ventilation of EEP thermocline and deep waters occurred synchronously during the last deglaciation. In addition, both gradual and rapid deglacial radiocarbon changes in these Pacific records are coeval with changes in the Atlantic records. This in-phase behaviour suggests that the Southern Ocean overturning was the dominant driver of changes in the Atlantic and Pacific ventilation during deglaciation.

Potential age model problems with marine Δ¹⁴C records have obscured our understanding of the role of the deep-ocean in deglacial atmospheric CO₂ rise. Here, the authors show that deglacial ventilation of EEP thermocline and deep waters occurred synchronously and was coeval with changes in Atlantic records.

Estimating the anomalous diffusion exponent for single particle tracking data with measurement errors - An alternative approach

Accurately characterizing the anomalous diffusion of a tracer particle has become a central issue in biophysics. However, measurement errors raise difficulty in the characterization of single trajectories, which is usually performed through the time-averaged mean square displacement (TAMSD). In this paper, we study a fractionally integrated moving average (FIMA) process as an appropriate model for anomalous diffusion data with measurement errors. We compare FIMA and traditional TAMSD estimators for the anomalous diffusion exponent. The ability of the FIMA framework to characterize dynamics in a wide range of anomalous exponents and noise levels through the simulation of a toy model (fractional Brownian motion disturbed by Gaussian white noise) is discussed. Comparison to the TAMSD technique, shows that FIMA estimation is superior in many scenarios. This is expected to enable new measurement regimes for single particle tracking (SPT) experiments even in the presence of high measurement errors.