Ocean Heat Content Responses to Changing Anthropogenic Aerosol Forcing Strength: Regional and Multi‐Decadal Variability

The causes of decadal variations in global warming are poorly understood, however it is widely understood that variations in ocean heat content (OHC) are linked with variations in surface warming. To investigate the forced response of OHC to anthropogenic aerosols (AA), we use an ensemble of historical simulations, which were carried out using a range of anthropogenic aerosol forcing magnitudes in a CMIP6-era global circulation model. We find that the centennial scale linear trends in historical OHC are significantly sensitive to AA forcing magnitude (−3.0 ± 0.1 × 10 5 (J m −3 century −1 )/(W m −2 ), R 2 = 0.99), but interannual to multi-decadal variability in global OHC appear largely independent of AA forcing magnitude. Comparison with observations find consistencies in different depth ranges and at different time scales with all but the strongest aerosol forcing magnitude, at least partly due to limited observational accuracy. We find that OHC is sensitive to aerosol forcing magnitude across much of the tropics and sub-tropics, and stronger negative forcing induces more ocean cooling. The polar regions and North Atlantic show the strongest heat content trends, and also show the strongest dependence on aerosol forcing magnitude. However, the OHC response to increasing aerosol forcing magnitude in the North Atlantic and Southern Ocean is either dominated by internal variability, or strongly state dependent, showing different behavior in different time periods. Our results suggest the response to aerosols in these regions is a complex combination of influences from ocean transport, atmospheric forcings, and sea ice.