Greenhouse gases and aerosols play a major role in controlling global climate change. Greenhouse gases drive a radiative imbalance which warms the ocean, while aerosols cool the ocean. Since 1980, the effective radiation felt by the planet due to anthropogenic aerosols has levelled off, global ocean cooling due to aerosols has decelerated, and greenhouse gas-driven ocean warming has accelerated. We explore the deceleration of aerosol-driven ocean cooling by quantifying a time- and spatially-varying ocean heat uptake efficiency, defined as the change in the rate of global ocean heat storage per degree of cooling surface temperature. In aerosol-only simulations, ocean heat uptake efficiency has decreased by 69% since the 1900s. The tropics and sub-tropics have driven this decrease, while the coldest fraction of the ocean continues to sustain cooling and high ocean heat uptake efficiency. Our results identify a growing trend towards less efficient ocean cooling due to aerosols.

The ocean has absorbed approximately 90% of the accumulated heat in the climate system since 1970. As global warming accelerates, understanding ocean heat content changes and tracing these to surface heat input is becoming increasingly important. We introduce a novel tracer-percentile framework in which we organise the ocean into temperature percentiles from warmest to coldest, allowing us to trace changes in ocean temperature to changes in air-sea heat fluxes and mixing. Applying this framework to observations and historical CMIP6 simulations, we find that 40% of heat uptake between 1970 and 2014 occurs in the warmest 10% ocean volume. However, the coolest 90% ocean volume outcrops over 23% of the ocean’s surface area, implying a disproportionately large heat input per unit area. Additionally, a cold bias in the CMIP6 models is traced to inaccurate sea surface temperatures and surface heat fluxes into the warmest 5 – 20% ocean volume.

The largest contributor to the planetary energy imbalance is well-mixed greenhouse gases (GHGs), which are partially offset by poorly-mixed (and thus northern mid-latitude dominated) anthropogenic aerosols (AAs). To isolate the effects of GHGs and AAs, we analyze data from the CMIP5 historical (i.e. all natural and anthropogenic forcing) and single forcing (GHG-only and AA-only) experiments. Over the duration of the historical experiment (1861-2005) excess heat uptake at the top of the atmosphere and ocean surface occurs almost exclusively in the Southern Hemisphere, with AAs canceling the influence of GHGs in the Northern Hemisphere. This interhemispheric asymmetry in surface heat uptake is eliminated by a northward oceanic transport of excess heat, as there is little hemispheric difference in historical ocean heat storage after accounting for ocean volume. Data from the 1pctCO2 and RCP 8.5 experiments suggests that the future storage of excess heat will be skewed towards the Northern Hemisphere oceans. PLAIN LANGUAGE SUMMARY Climate change is fundamentally an energy balance problem. Due to the influence of greenhouse gas (GHG) emissions, the amount of solar energy absorbed by Earth is currently greater than the amount of energy radiated to space. This energy imbalance is partially offset by particulate matter released into the atmosphere from burning fossil fuels (anthropogenic aerosols; AAs), which is most concentrated in the northern mid-latitudes. In this study, model simulations of the historical and future climate are compared to single forcing simulations that apply only GHG or AA emissions. We find that the historical uptake of excess heat is strongly skewed towards the Southern Hemisphere because AAs cancel the influence of GHGs in the Northern Hemisphere. The oceanic storage of that heat shows little difference between the hemispheres due to a strong northward transport of excess heat. In future, the models suggest excess heat storage will skew towards the Northern Hemisphere. KEY POINTS - Single forcing simulations of the historical (1861-2005) climate suggest anthropogenic aerosols cause most excess heat uptake to occur in the Southern Hemisphere - This interhemispheric asymmetry in heat uptake is eliminated (after accounting for ocean volume) by a northward oceanic transport of excess heat - In future, the storage of excess heat will be skewed towards the Northern Hemisphere

Coupled climate models are prone to ‘drift’ (long-term unforced trends in state variables) due to incomplete spin-up and non-closure of the global mass and energy budgets. Here we assess model drift and the associated conservation of energy, mass and salt in CMIP6 and CMIP5 models. For most models, drift in globally-integrated ocean mass and heat content represents a small but non-negligible fraction of recent historical trends, while drift in atmospheric water vapor is negligible. Model drift tends to be much larger in time-integrated ocean heat and freshwater flux, net top-of-the-atmosphere radiation (netTOA) and moisture flux into the atmosphere (evaporation minus precipitation), indicating a substantial leakage of mass and energy in the simulated climate system. Most models are able to achieve approximate energy budget closure after drift is removed, but ocean mass budget closure eludes a number of models even after de-drifting and none achieve closure of the atmospheric moisture budget. The magnitude of the drift in the CMIP6 ensemble represents an improvement over CMIP5 in some cases (salinity and time-integrated netTOA) but is worse (time-integrated ocean freshwater and atmospheric moisture fluxes) or little changed (ocean heat content, ocean mass and time-integrated ocean heat flux) for others, while closure of the ocean mass and energy budgets after drift removal has improved.

Weather and climate science has undergone a computational revolution in recent decades, to the point where all modern research relies heavily on software and code. Despite this profound change in the research methods employed by weather and climate scientists, the reporting of computational results has changed very little in relevant academic journals. This lag has led to something of a reproducibility crisis, whereby it is impossible to replicate and verify most of today’s published computational results. While it is tempting to simply decry the slow response of journals and funding agencies in the face of this crisis, there are very few examples of reproducible weather and climate research upon which to base new communication standards. In an attempt to address this deficiency, this essay describes a procedure for reporting computational results that was employed in a recent _Journal of Climate_ paper. The procedure was developed to be consistent with recommended computational best practices and seeks to minimize the time burden on authors, which has been identified as the most important barrier to publishing code. It should provide a starting point for weather and climate scientists looking to publish reproducible research, and it is proposed that journals could adopt the procedure as a minimum standard.

The Pacific-South American (PSA) pattern is an important mode of climate variability in the mid-to-high southern latitudes. It is widely recognized as the primary mechanism by which the El Niño-Southern Oscillation (ENSO) influences the south-east Pacific and south-west Atlantic, and in recent years has also been suggested as a mechanism by which longer-term tropical sea surface temperature trends can influence the Antarctic climate. This study presents a novel methodology for objectively identifying the PSA pattern. By rotating the global coordinate system such that the equator (a great circle) traces the approximate path of the pattern, the identification algorithm utilizes Fourier analysis as opposed to a traditional Empirical Orthogonal Function approach. The climatology arising from the application of this method to ERA-Interim reanalysis data reveals that the PSA pattern has a strong influence on temperature and precipitation variability over West Antarctica and the Antarctic Peninsula, and on sea ice variability in the adjacent Amundsen, Bellingshausen and Weddell Seas. Identified seasonal trends towards the negative phase of the PSA pattern are consistent with warming observed over the Antarctic Peninsula during autumn, but are inconsistent with observed winter warming over West Antarctica. Only a weak relationship is identified between the PSA pattern and ENSO, which suggests that the pattern might be better conceptualized as preferred regional atmospheric response to various external (and internal) forcings.

Southern Hemisphere mid-to-upper tropospheric planetary wave activity is characterized by the superposition of two zonally-oriented, quasi-stationary waveforms: zonal wavenumber one (ZW1) and zonal wavenumber three (ZW3). Previous studies have tended to consider these waveforms in isolation and with the exception of those studies relating to sea ice, little is known about their impact on regional climate variability. We take a novel approach to quantifying the combined influence of ZW1 and ZW3, using the strength of the hemispheric meridional flow as a proxy for zonal wave activity. Our methodology adapts the wave envelope construct routinely used in the identification of synoptic-scale Rossby wave packets and improves on existing approaches by allowing for variations in both wave phase and amplitude. While ZW1 and ZW3 are both prominent features of the climatological circulation, the defining feature of highly meridional hemispheric states is an enhancement of the ZW3 component. Composites of the mean surface conditions during these highly meridional, ZW3-like anomalous states (i.e. months of strong planetary wave activity) reveal large sea ice anomalies over the Amundsen and Bellingshausen Seas during autumn and along much of the East Antarctic coastline throughout the year. Large precipitation anomalies in regions of significant topography (e.g. New Zealand, Patagonia, coastal Antarctica) and anomalously warm temperatures over much of the Antarctic continent were also associated with strong planetary wave activity. The latter has potentially important implications for the interpretation of recent warming over West Antarctica and the Antarctic Peninsula.