The 2019/20 Black Summer bushfire disaster in southeast Australia was unprecedented: the extensive area of forest burnt, the radiative power of the fires, and the extraordinary number of fires that developed into extreme pyroconvective events were all unmatched in the historical record. Australia’s hottest and driest year on record, 2019, was characterised by exceptionally dry fuel loads that primed the landscape to burn when exposed to dangerous fire weather and ignition. The combination of climate variability and long-term climate trends generated the climate extremes experienced in 2019, and the compounding effects of two or more modes of climate variability in their fire-promoting phases (as occurred in 2019) has historically increased the chances of large forest fires occurring in southeast Australia. Palaeoclimate evidence also demonstrates that fire-promoting phases of tropical Pacific and Indian ocean variability are now unusually frequent compared with natural variability in pre-industrial times. Indicators of forest fire danger in southeast Australia have already emerged outside of the range of historical experience, suggesting that projections made more than a decade ago that increases in climate-driven fire risk would be detectable by 2020, have indeed eventuated. The multiple climate change contributors to fire risk in southeast Australia, as well as the observed non-linear escalation of fire extent and intensity, raise the likelihood that fire events may continue to rapidly intensify in the future. Improving local and national adaptation measures while also pursuing ambitious global climate change mitigation efforts would provide the best strategy for limiting further increases in fire risk in southeast Australia.

The modeling of paleoclimate, using physically based tools, is increasingly seen as a strong out-of-sample test of the models that are used for the projection of future climate changes. New to the Coupled Model Intercomparison Project (CMIP6) is the Tier 1 Last Interglacial experiment for 127 000 years ago (lig127k), designed to address the climate responses to stronger orbital forcing than the midHolocene experiment, using the same state-of-the-art models as for the future and following a common experimental protocol. Here we present a first analysis of a multi-model ensemble of 17 climate models, all of which have completed the CMIP6 DECK (Diagnostic, Evaluation and Characterization of Klima) experiments. The equilibrium climate sensitivity (ECS) of these models varies from 1.8 to 5.6 ^∘C. The seasonal character of the insolation anomalies results in strong summer warming over the Northern Hemisphere continents in the lig127k ensemble as compared to the CMIP6 piControl and much-reduced minimum sea ice in the Arctic. The multi-model results indicate enhanced summer monsoonal precipitation in the Northern Hemisphere and reductions in the Southern Hemisphere. These responses are greater in the lig127k than the CMIP6 midHolocene simulations as expected from the larger insolation anomalies at 127 than 6 ka. New synthesis for surface temperature and precipitation, targeted for 127 ka, have been developed for comparison to the multi-model ensemble. The lig127k model ensemble and data reconstructions are in good agreement for summer temperature anomalies over Canada, Scandinavia, and the North Atlantic and for precipitation over the Northern Hemisphere continents. The model–data comparisons and mismatches point to further study of the sensitivity of the simulations to uncertainties in the boundary conditions and of the uncertainties and sparse coverage in current proxy reconstructions. The CMIP6–Paleoclimate Modeling Intercomparison Project (PMIP4) lig127k simulations, in combination with the proxy record, improve our confidence in future projections of monsoons, surface temperature, and Arctic sea ice, thus providing a key target for model evaluation and optimization.

The early part of the last deglaciation is characterised by a ~40 ppm atmospheric CO₂ rise occurring in two abrupt phases. The underlying mechanisms driving these increases remain a subject of intense debate. Here, we successfully reproduce changes in CO₂, δ¹³C and Δ¹⁴C as recorded by paleo-records during Heinrich stadial 1 (HS1). We show that HS1 CO₂ increase can be explained by enhanced Southern Ocean upwelling of carbon-rich Pacific deep and intermediate waters, resulting from intensified Southern Ocean convection and Southern Hemisphere (SH) westerlies. While enhanced Antarctic Bottom Water formation leads to a millennial CO₂ outgassing, intensified SH westerlies induce a multi-decadal atmospheric CO₂ rise. A strengthening of SH westerlies in a global eddy-permitting ocean model further supports a multi-decadal CO₂ outgassing from the Southern Ocean. Our results highlight the crucial role of SH westerlies in the global climate and carbon cycle system with important implications for future climate projections.

Atmospheric CO₂ was ∼90 ppmv lower at the Last Glacial Maximum (LGM) compared to the late Holocene, but the mechanisms responsible for this change remain elusive. Here we employ a carbon isotope‐enabled Earth System Model to investigate the role of ocean circulation in setting the LGM oceanic δ¹³C distribution, thereby improving our understanding of glacial/interglacial atmospheric CO₂ variations. We find that the mean ocean δ¹³C change can be explained by a 378 ± 88 Gt C(2σ) smaller LGM terrestrial carbon reservoir compared to the Holocene. Critically, in this model, differences in the oceanic δ¹³C spatial pattern can only be reconciled with a LGM ocean circulation state characterized by a weak (10–15 Sv) and relatively shallow (2000–2500 m) North Atlantic Deep Water cell, reduced Antarctic Bottom Water transport (≤10 Sv globally integrated), and relatively weak (6–8 Sv) and shallow (1000–1500 m) North Pacific Intermediate Water formation. This oceanic circulation state is corroborated by results from the isotope‐enabled Bern3D ocean model and further confirmed by high LGM ventilation ages in the deep ocean, particularly in the deep South Atlantic and South Pacific. This suggests a poorly ventilated glacial deep ocean which would have facilitated the sequestration of carbon lost from the terrestrial biosphere and atmosphere.

An increasing number of Southern Ocean models now include Antarctic ice-shelf cavities, and simulate thermodynamics at the ice-shelf/ocean interface. This adds another level of complexity to Southern Ocean simulations, as ice shelves interact directly with the ocean and indirectly with sea ice. Here, we present the first model intercomparison and evaluation of present-day ocean/sea-ice/ice-shelf interactions, as simulated by two models: a circumpolar Antarctic configuration of MetROMS (ROMS: Regional Ocean Modelling System coupled to CICE: Community Ice CodE) and the global model FESOM (Finite Element Sea-ice Ocean Model), where the latter is run at two different levels of horizontal resolution. From a circumpolar Antarctic perspective, we compare and evaluate simulated ice-shelf basal melting and sub-ice-shelf circulation, as well as sea-ice properties and Southern Ocean water mass characteristics as they influence the sub-ice-shelf processes. Despite their differing numerical methods, the two models produce broadly similar results and share similar biases in many cases. Both models reproduce many key features of observations but struggle to reproduce others, such as the high melt rates observed in the small warm-cavity ice shelves of the Amundsen and Bellingshausen seas. Several differences in model design show a particular influence on the simulations. For example, FESOM's greater topographic smoothing can alter the geometry of some ice-shelf cavities enough to affect their melt rates; this improves at higher resolution, since less smoothing is required. In the interior Southern Ocean, the vertical coordinate system affects the degree of water mass erosion due to spurious diapycnal mixing, with MetROMS' terrain-following coordinate leading to more erosion than FESOM's z coordinate. Finally, increased horizontal resolution in FESOM leads to higher basal melt rates for small ice shelves, through a combination of stronger circulation and small-scale intrusions of warm water from offshore.