North Atlantic meltwater during Heinrich Stadial 1 drives wetter climate with more atmospheric rivers in western North America.

Atmospheric rivers (ARs) bring concentrated rainfall and flooding to the western United States (US) and are hypothesized to have supported sustained hydroclimatic changes in the past. However, their ephemeral nature makes it challenging to document ARs in climate models and estimate their contribution to hydroclimate changes recorded by time-averaged paleoclimate archives. We present new climate model simulations of Heinrich Stadial 1 (HS1; 16,000 years before the present), an interval characterized by widespread wetness in the western US, that demonstrate increased AR frequency and winter precipitation sourced from the southeastern North Pacific. These changes are amplified with freshwater fluxes into the North Atlantic, indicating that North Atlantic cooling associated with weakened Atlantic Meridional Overturning Circulation (AMOC) is a key driver of HS1 climate in this region. As recent observations suggest potential weakening of AMOC, our identified connection between North Atlantic climate and northeast Pacific AR activity has implications for future western US hydroclimate.

Tropical mountain ice core δ18O: A Goldilocks indicator for global temperature change.

Very high tropical alpine ice cores provide a distinct paleoclimate record for climate changes in the middle and upper troposphere. However, the climatic interpretation of a key proxy, the stable water oxygen isotopic ratio in ice cores (δ18Oice), remains an outstanding problem. Here, combining proxy records with climate models, modern satellite measurements, and radiative-convective equilibrium theory, we show that the tropical δ18Oice is an indicator of the temperature of the middle and upper troposphere, with a glacial cooling of -7.35° ± 1.1°C (66% CI). Moreover, it severs as a "Goldilocks-type" indicator of global mean surface temperature change, providing the first estimate of glacial stage cooling that is independent of marine proxies as -5.9° ± 1.2°C. Combined with all estimations available gives the maximum likelihood estimate of glacial cooling as -5.85° ± 0.51°C.

Human-induced ecological cascades: Extinction, restoration, and rewilding in the Galápagos highlands.

Oceanic islands support unique biotas but often lack ecological redundancy, so that the removal of a species can have a large effect on the ecosystem. The larger islands of the Galápagos Archipelago once had one or two species of giant tortoise that were the dominant herbivore. Using paleoecological techniques, we investigate the ecological cascade on highland ecosystems that resulted from whalers removing many thousands of tortoises from the lowlands. We hypothesize that the seasonal migration of a now-extinct tortoise species to the highlands was curtailed by decreased intraspecific competition. We find the trajectory of plant community dynamics changed within a decade of the first whaling vessels visiting the islands. Novel communities established, with a previously uncommon shrub, Miconia, replacing other shrubs of the genera Alternanthera and Acalypha. It was, however, the introduction of cattle and horses that caused the local extirpation of plant species, with the most extreme impacts being evident after c. 1930. This modified ecology is considered the natural state of the islands and has shaped subsequent conservation policy and practice. Restoration of El Junco Crater should emphasize exclusion of livestock, rewilding with tortoises, and expanding the ongoing plantings of Miconia to also include Acalypha and Alternanthera.

Description and Demonstration of the Coupled Community Earth System Model v2 - Community Ice Sheet Model v2 (CESM2-CISM2).

Earth system/ice-sheet coupling is an area of recent, major Earth System Model (ESM) development. This work occurs at the intersection of glaciology and climate science and is motivated by a need for robust projections of sea-level rise. The Community Ice Sheet Model version 2 (CISM2) is the newest component model of the Community Earth System Model version 2 (CESM2). This study describes the coupling and novel capabilities of the model, including: (1) an advanced energy-balance-based surface mass balance calculation in the land component with downscaling via elevation classes; (2) a closed freshwater budget from ice sheet to the ocean from surface runoff, basal melting, and ice discharge; (3) dynamic land surface types; and (4) dynamic atmospheric topography. The Earth system/ice-sheet coupling is demonstrated in a simulation with an evolving Greenland Ice Sheet (GrIS) under an idealized high CO2 scenario. The model simulates a large expansion of ablation areas (where surface ablation exceeds snow accumulation) and a large increase in surface runoff. This results in an elevated freshwater flux to the ocean, as well as thinning of the ice sheet and area retreat. These GrIS changes result in reduced Greenland surface albedo, changes in the sign and magnitude of sensible and latent heat fluxes, and modified surface roughness and overall ice sheet topography. Representation of these couplings between climate and ice sheets is key for the simulation of ice and climate interactions.

An Efficient Ice Sheet/Earth System Model Spin-up Procedure for CESM2-CISM2: Description, Evaluation, and Broader Applicability.

Spinning up a highly complex, coupled Earth system model (ESM) is a time consuming and computationally demanding exercise. For models with interactive ice sheet components, this becomes a major challenge, as ice sheets are sensitive to bidirectional feedback processes and equilibrate over glacial timescales of up to many millennia. This work describes and demonstrates a computationally tractable, iterative procedure for spinning up a contemporary, highly complex ESM that includes an interactive ice sheet component. The procedure alternates between a computationally expensive coupled configuration and a computationally cheaper configuration where the atmospheric component is replaced by a data model. By periodically regenerating atmospheric forcing consistent with the coupled system, the data atmosphere remains adequately constrained to ensure that the broader model state evolves realistically. The applicability of the method is demonstrated by spinning up the preindustrial climate in the Community Earth System Model Version 2 (CESM2), coupled to the Community Ice Sheet Model Version 2 (CISM2) over Greenland. The equilibrium climate state is similar to the control climate from a coupled simulation with a prescribed Greenland ice sheet, indicating that the iterative procedure is consistent with a traditional spin-up approach without interactive ice sheets. These results suggest that the iterative method presented here provides a faster and computationally cheaper method for spinning up a highly complex ESM, with or without interactive ice sheet components. The method described here has been used to develop the climate/ice sheet initial conditions for transient, ice sheet-enabled simulations with CESM2-CISM2 in the Coupled Model Intercomparison Project Phase 6 (CMIP6).