Rainfall will become more variable from day to day and year to year in global wet regions under future warming.

A physically based method for parameterizing the role of subgrid-scale turbulence in the production and maintenance of supercooled liquid water and mixed-phase clouds is presented. The approach used is to simplify the dynamics of supersaturation fluctuations to a stochastic differential equation that can be solved analytically, giving increments to the prognostic liquid cloud fraction and liquid water content fields in a general circulation model (GCM). Elsewhere, it has been demonstrated that the approach captures the properties of decameter-resolution large-eddy simulations of a turbulent mixed-phase environment. In this paper, it is shown that it can be implemented in a GCM, and the effects that this has on Southern Ocean biases and on Arctic stratus are investigated.

The Earth system models (ESMs) that participated in the sixth Coupled Model Intercomparison Project (CMIP6) tend to simulate excessive cooling in surface air temperature (TAS) between 1960 and 1990. The anomalous cooling is pronounced over the Northern Hemisphere (NH) midlatitudes, coinciding with the rapid growth of anthropogenic sulfur dioxide (SO₂) emissions, the primary precursor of atmospheric sulfate aerosols. Structural uncertainties between ESMs have a larger impact on the anomalous cooling than internal variability. Historical simulations with and without anthropogenic aerosol emissions indicate that the anomalous cooling in the ESMs is attributed to the higher aerosol burden in these models. The aerosol forcing sensitivity, estimated as the outgoing shortwave radiation (OSR) response to aerosol concentration changes, cannot well explain the diversity of pothole cooling (PHC) biases in the ESMs. The relative contributions to aerosol forcing sensitivity from aerosol–radiation interactions (ARIs) and aerosol–cloud interactions (ACIs) can be estimated from CMIP6 simulations. We show that even when the aerosol forcing sensitivity is similar between ESMs, the relative contributions of ARI and ACI may be substantially different. The ACI accounts for between 64 % and 87 % of the aerosol forcing sensitivity in the models and is the main source of the aerosol forcing sensitivity differences between the ESMs. The ACI can be further decomposed into a cloud-amount term (which depends linearly on cloud fraction) and a cloud-albedo term (which is independent of cloud fraction, to the first order), with the cloud-amount term accounting for most of the inter-model differences.

Projected changes of future precipitation extremes exhibit substantial uncertainties among climate models, posing grand challenges to climate actions and adaptation planning. Practical methods for narrowing the projection uncertainty remain elusive. Here, using large model ensembles, we show that the uncertainty in projections of future extratropical extreme precipitation is significantly correlated with the model representations of present-day precipitation variability. Models with weaker present-day precipitation variability tend to project larger increases in extreme precipitation occurrences under a given global warming increment. This relationship can be explained statistically using idealized distributions for precipitation. This emergent relationship provides a powerful constraint on future projections of extreme precipitation from observed present-day precipitation variability, which reduces projection uncertainty by 20–40% over extratropical regions. Because of the widespread impacts of extreme precipitation, this has not only provided useful insights into understanding uncertainties in current model projections, but is also expected to bring potential socio-economic benefits in climate change adaptation planning.

High-resolution simulations of a Southern Ocean cyclone are compared to satellite-derived observations of liquid water path, cloud-top properties, and top-of-atmosphere radiative fluxes. The focus is on the cold-air-outflow region, where there are contributions to the hydrological budget from the microphysical growth of ice particles by riming and vapor deposition and transport by turbulent mixing. The sensitivity of the simulation to the parameterization of these processes is tested and the relative importance of ice-nucleation temperature is identified. It is shown that ice-phase microphysics is a key factor determining the phase composition of Southern Ocean clouds and physically reasonable parameterization changes are identified that affect the liquid water content of these clouds. The information gained from the sensitivity tests is applied to global model development, where it is shown that a modification to the riming parameterization improves climate mean-state biases in the Southern Ocean region.

A key contribution to the latest generation of climate projections for the UK (UKCP18) was a perturbed parameter ensemble (PPE) of global coupled models based on HadGEM3-GC3.05. Together with 13 CMIP5 simulations, this PPE provides users with a dataset that samples modelling uncertainty and is ideal for use in impacts studies. Evaluations of global mean surface temperatures for this PPE have shown twenty-first century warming rates consistently at the top end of the CMIP5 range. Here we investigate one potential contributory factor to this lack of spread: that the methodology to select plausible members from a larger, related PPE of atmosphere-only experiments preferentially ruled out those predicted to have more negative climate feedbacks (i.e. lower climate sensitivities). We confirm that this is indeed the case. We show that performance in extratropical long-wave cloud forcing played a key role in this by constraining ice cloud parameters, which in turn constrained the feedback distribution (though causal links are not established). The relatively weak relationship driving this constraint is shown to arise from stronger relationships for the long-wave and short-wave cloud feedback components, which largely cancel out due to changes in tropical high clouds. Moreover, we show that the strength of these constraints is due to a structural bias in extratropical long-wave cloud forcing across the PPE. We discuss how choices made in the methodology to pick the plausible PPE members may result in an overly strong constraint when there is a structural bias and possible improvements to this methodology for the future.