We document the development of the first version of the U.K. Earth System Model UKESM1.The model represents a major advance on its predecessor HadGEM2-ES, with enhancements to all component models and new feedback mechanisms. These include a new core physical model with a well-resolved stratosphere; terrestrial biogeochemistry with coupled carbon and nitrogen cycles and enhanced land management; tropospheric-stratospheric chemistry allowing the holistic simulation of radiative forcing from ozone, methane, and nitrous oxide; two-moment, five-species, modal aerosol; and ocean biogeochemistry with two-way coupling to the carbon cycle and atmospheric aerosols. The complexity of coupling between the ocean, land, and atmosphere physical climate and biogeochemical cycles in UKESM1 is unprecedented for an Earth system model. We describe in detail the process by which the coupled model was developed and tuned to achieve acceptable performance in key physical and Earth system quantities and discuss the challenges involved in mitigating biases in a model with complex connections between its components. Overall, the model performs well, with a stable pre-industrial state and good agreement with observations in the latter period of its historical simulations. However, global mean surface temperature exhibits stronger-than-observed cooling from 1950 to 1970, followed by rapid warming from 1980 to 2014. Metrics from idealized simulations show a high climate sensitivity relative to previous generations of models: Equilibrium climate sensitivity is 5.4 K, transient climate response ranges from 2.68 to 2.85 K, and transient climate response to cumulative emissions is 2.49 to 2.66 K TtC −1 . Plain Language SummaryWe describe the development and behavior of UKESM1, a novel climate model that includes improved representations of processes in the atmosphere, ocean, and on land. These processes are inter-related: For example, dust is produced on the land and blown up into the atmosphere where it affects the amount of sunlight falling on Earth. Dust can also be dissolved in the ocean, where it affects marine life. This in turn changes both the amount of carbon dioxide absorbed by the ocean and the material emitted from the surface into the atmosphere, which has an affect on the formation of clouds. UKESM1 includes many processes and interactions such as these, giving it a high level of complexity. Ensuring realistic process behavior is a major challenge in the development of our model, and we have carefully tested this. UKESM1 performs well, correctly exhibiting stable results from a continuous pre-industrial simulation (used to provide a reference for future experiments) and showing good agreement

Abstract. We describe Global Atmosphere 7.0 and Global Land 7.0 (GA7.0/GL7.0), the latest science configurations of the Met Office Unified Model (UM) and the Joint UK Land Environment Simulator (JULES) land surface model developed for use across weather and climate timescales. GA7.0 and GL7.0 include incremental developments and targeted improvements that, between them, address four critical errors identified in previous configurations: excessive precipitation biases over India, warm and moist biases in the tropical tropopause layer (TTL), a source of energy non-conservation in the advection scheme and excessive surface radiation biases over the Southern Ocean. They also include two new parametrisations, namely the UK Chemistry and Aerosol (UKCA) GLOMAP-mode (Global Model of Aerosol Processes) aerosol scheme and the JULES multi-layer snow scheme, which improve the fidelity of the simulation and were required for inclusion in the Global Atmosphere/Global Land configurations ahead of the 6th Coupled Model Intercomparison Project (CMIP6). In addition, we describe the GA7.1 branch configuration, which reduces an overly negative anthropogenic aerosol effective radiative forcing (ERF) in GA7.0 whilst maintaining the quality of simulations of the present-day climate. GA7.1/GL7.0 will form the physical atmosphere/land component in the HadGEM3–GC3.1 and UKESM1 climate model submissions to the CMIP6.

Summary (149 words of referenced text): 46The climate impact of aerosols is highly uncertain owing primarily to their poorly quantified 47 influence on cloud properties. During 2014-15, a fissure eruption in Holuhraun (Iceland) 48 emitted huge quantities of sulphur dioxide, resulting in significant reductions in liquid cloud 49 droplet size. Using satellite observations and detailed modelling, we estimate a global mean 50 radiative forcing from the resulting aerosol-induced cloud brightening for the time of the 51 eruption of around -0.2 W.m -2 . Changes in cloud amount or liquid water path are 52 undetectable, indicating that these aerosol-cloud indirect effects are modest. It supports the 53 idea that cloud systems are well buffered against aerosol changes as only impacts on cloud 54 effective radius appear relevant from a climate perspective, thus providing a strong constraint 55 on aerosol-cloud interactions. This result will reduce uncertainties in future climate 56 projections as we are able to reject the results from climate models with an excessive liquid 57 water path response. 58 59Main Text: (3103 words of referenced text, including concluding paragraph) 60 The 2014-15 eruption at Holuhraun (486 words of referenced text): 61Anthropogenic emissions that affect climate are not just confined to greenhouse gases. 62Sulphur dioxide and other pollutants form atmospheric aerosols that can scatter and absorb 63 sunlight and can influence the properties of clouds, modulating the Earth-atmosphere energy 64 balance. Aerosols act as cloud condensation nuclei (CCN); an increase in CCN translates into 65 a higher number of smaller, more reflective cloud droplets that scatter more sunlight back to 66 space 1 (the ÔfirstÕ indirect effect of aerosols). Smaller cloud droplets decrease the efficiency 67 of collision-coalescence processes that are pivotal in rain initiation, thus aerosol-influenced 68 clouds may retain more liquid water and extend coverage/lifetime 2,3 (the ÔsecondÕ or Ôcloud 69 lifetimeÕ indirect effect). Aerosols usually co-vary with key environmental variables making 70 it difficult to disentangle aerosol-cloud impacts from meteorological variability [4][5][6] . 71Additionally, clouds themselves are complex transient systems subject to dynamical 72 feedbacks (e.g. cloud top entrainment/evaporation, invigoration of convection) which 73 influence cloud response [7][8][9][10][11][12] . These aspects present great challenges in evaluating and 74 constraining aerosol-cloud interactions (ACI) in General Circulation Models (GCM) 13-17 , 75 with particular contentious debate surrounding the relative importance of these feedback 76 mechanisms. 77Nonetheless, anthropogenic aerosol emissions are thought to cool the Earth via indirect 78 effects 17 , but the uncertainty ranges from -1.2 to -0.0 W.m -2 (90% confidence interval) due to 79 i) a lack of characterization of the pre-industrial aerosol state 15,18,19 , and ii) model parametric 80 and structural errors in representing cloud responses to aerosol chan...

Abstract. Here we present a description of the UKCA StratTrop chemical mechanism, which is used in the UKESM1 Earth system model for CMIP6. The StratTrop chemical mechanism is a merger of previously well-evaluated tropospheric and stratospheric mechanisms, and we provide results from a series of bespoke integrations to assess the overall performance of the model. We find that the StratTrop scheme performs well when compared to a wide array of observations. The analysis we present here focuses on key components of atmospheric composition, namely the performance of the model to simulate ozone in the stratosphere and troposphere and constituents that are important for ozone in these regions. We find that the results obtained for tropospheric ozone and its budget terms from the use of the StratTrop mechanism are sensitive to the host model; simulations with the same chemical mechanism run in an earlier version of the MetUM host model show a range of sensitivity to emissions that the current model does not fall within. Whilst the general model performance is suitable for use in the UKESM1 CMIP6 integrations, we note some shortcomings in the scheme that future targeted studies will address.

Abstract. A 1200×1200 km2 area of the tropical South Atlantic Ocean near Ascension Island is studied with the HadGEM climate model at convection-permitting and global resolutions for a 10-day case study period in August 2016. During the simulation period, a plume of biomass burning smoke from Africa moves into the area and mixes into the clouds. At Ascension Island, this smoke episode was the strongest of the 2016 fire season.The region of interest is simulated at 4 km resolution, with no parameterised convection scheme. The simulations are driven by, and compared to, the global model. For the first time, the UK Chemistry and Aerosol model (UKCA) is included in a regional model with prognostic aerosol number concentrations advecting in from the global model at the boundaries of the region.Fire emissions increase the total aerosol burden by a factor of 3.7 and cloud droplet number concentrations by a factor of 3, which is consistent with MODIS observations. In the regional model, the inversion height is reduced by up to 200 m when smoke is included. The smoke also affects precipitation, to an extent which depends on the model microphysics. The microphysical and dynamical changes lead to an increase in liquid water path of 60 g m−2 relative to a simulation without smoke aerosol, when averaged over the polluted period. This increase is uncertain, and smaller in the global model. It is mostly due to radiatively driven dynamical changes rather than precipitation suppression by aerosol.Over the 5-day polluted period, the smoke has substantial direct radiative effects of +11.4 W m−2 in the regional model, a semi-direct effect of −30.5 W m−2 and an indirect effect of −10.1 W m−2. Our results show that the radiative effects are sensitive to the structure of the model (global versus regional) and the parameterization of rain autoconversion. Furthermore, we simulate a liquid water path that is biased high compared to satellite observations by 22 % on average, and this leads to high estimates of the domain-averaged aerosol direct effect and the effect of the aerosol on cloud albedo. With these caveats, we simulate a large net cooling across the region, of −27.6 W m−2.