Results are presented of the GASS/EUCLIPSE single‐column model intercomparison study on the subtropical marine low‐level cloud transition. A central goal is to establish the performance of ...state‐of‐the‐art boundary‐layer schemes for weather and climate models for this cloud regime, using large‐eddy simulations of the same scenes as a reference. A novelty is that the comparison covers four different cases instead of one, in order to broaden the covered parameter space. Three cases are situated in the North‐Eastern Pacific, while one reflects conditions in the North‐Eastern Atlantic. A set of variables is considered that reflects key aspects of the transition process, making use of simple metrics to establish the model performance. Using this method, some longstanding problems in low‐level cloud representation are identified. Considerable spread exists among models concerning the cloud amount, its vertical structure, and the associated impact on radiative transfer. The sign and amplitude of these biases differ somewhat per case, depending on how far the transition has progressed. After cloud breakup the ensemble median exhibits the well‐known “too few too bright” problem. The boundary‐layer deepening rate and its state of decoupling are both underestimated, while the representation of the thin capping cloud layer appears complicated by a lack of vertical resolution. Encouragingly, some models are successful in representing the full set of variables, in particular, the vertical structure and diurnal cycle of the cloud layer in transition. An intriguing result is that the median of the model ensemble performs best, inspiring a new approach in subgrid parameterization.
Key Points
SCM simulations of low‐level cloud transitions are confronted with LES results
Longstanding problems are identified along with encouraging progress
The model‐ensemble median outperforms the individual models
ABSTRACT
We present results from 3D radiative-hydrodynamical simulations of HD 209458b with a fully coupled treatment of clouds using the EddySed code, critically, including cloud radiative feedback ...via absorption and scattering. We demonstrate that the thermal and optical structure of the simulated atmosphere is markedly different, for the majority of our simulations, when including cloud radiative effects, suggesting this important mechanism cannot be neglected. Additionally, we further demonstrate that the cloud structure is sensitive to not only the cloud sedimentation efficiency (termed fsed in EddySed), but also the temperature–pressure profile of the deeper atmosphere. We briefly discuss the large difference between the resolved cloud structures of this work, adopting a phase-equilibrium and parametrized cloud model, and our previous work incorporating a cloud microphysical model, although a fairer comparison where, for example, the same list of constituent condensates is included in both treatments is reserved for a future work. Our results underline the importance of further study into the potential condensate size distributions and vertical structures, as both strongly influence the radiative impact of clouds on the atmosphere. Finally, we present synthetic observations from our simulations reporting an improved match, over our previous cloud-free simulations, to the observed transmission, HST WFC3 emission, and 4.5 μm Spitzer phase curve of HD 209458b. Additionally, we find all our cloudy simulations have an apparent albedo consistent with observations.
Aims. To understand and compare the 3D atmospheric structure of HD 209458 b and HD 189733 b, focusing on the formation and distribution of cloud particles, as well as their feedback on the dynamics ...and thermal profile. Methods. We coupled the 3D Met Office Unified Model (UM), including detailed treatments of atmospheric radiative transfer and dynamics, to a kinetic cloud formation scheme. The resulting model self-consistently solves for the formation of condensation seeds, surface growth and evaporation, gravitational settling and advection, cloud radiative feedback via absorption, and crucially, scattering. We used fluxes directly obtained from the UM to produce synthetic spectral energy distributions and phase curves. Results. Our simulations show extensive cloud formation in both HD 209458 b and HD 189733 b. However, cooler temperatures in the latter result in higher cloud particle number densities. Large particles, reaching 1 μm in diameter, can form due to high particle growth velocities, and sub-μm particles are suspended by vertical flows leading to extensive upper-atmosphere cloud cover. A combination of meridional advection and efficient cloud formation in cooler high latitude regions, results in enhanced cloud coverage for latitudes above 30° and leads to a zonally banded structure for all our simulations. The cloud bands extend around the entire planet, for HD 209458 b and HD 189733 b, as the temperatures, even on the day side, remain below the condensation temperature of silicates and oxides. Therefore, the simulated optical phase curve for HD 209458 b shows no “offset”, in contrast to observations. Efficient scattering of stellar irradiation by cloud particles results in a local maximum cooling of up to 250 K in the upper atmosphere, and an advection-driven fluctuating cloud opacity causes temporal variability in the thermal emission. The inclusion of this fundamental cloud-atmosphere radiative feedback leads to significant differences with approaches neglecting these physical elements, which have been employed to interpret observations and determine thermal profiles for these planets. This suggests that readers should be cautious of interpretations neglecting such cloud feedback and scattering, and that the subject merits further study.
ABSTRACT
We present results of 3D hydrodynamical simulations of HD209458b including a coupled, radiatively active cloud model (eddysed). We investigate the role of the mixing by replacing the default ...convective treatment used in previous works with a more physically relevant mixing treatment (Kzz) based on global circulation. We find that uncertainty in the efficiency of sedimentation through the sedimentation factor fsed plays a larger role in shaping cloud thickness and its radiative feedback on the local gas temperatures – e.g. hotspot shift and day-to-night side temperature gradient – than the switch in mixing treatment. We demonstrate using our new mixing treatments that simulations with cloud scales that are a fraction of the pressure scale height improve agreement with the observed transmission spectra, the emission spectra, and the Spitzer 4.5 µm phase curve, although our models are still unable to reproduce the optical and ultraviolet transmission spectra. We also find that the inclusion of cloud increases the transit asymmetry in the optical between the east and west limbs, although the difference remains small ($\lesssim 1{{\ \rm per\ cent}}$).
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).
The main aim of this two-part study is to use a perturbed parameter ensemble (PPE) to select plausible and diverse variants of a relatively expensive climate model for use in climate projections. In ...this first part, the extent to which climate biases develop at weather forecast timescales is assessed with two PPEs, which are based on 5-day forecasts and 10-year simulations with a relatively coarse resolution (N96) atmosphere-only model. Both ensembles share common parameter combinations and strong emergent relationships are found for a wide range of variables between the errors on two timescales. These relationships between the PPEs are demonstrated at several spatial scales from global (using mean square errors), to regional (using pattern correlations), and to individual grid boxes where a large fraction of them show positive correlations. The study confirms more robustly than in previous studies that investigating the errors on weather timescales provides an affordable way to identify and filter out model variants that perform poorly at short timescales and are likely to perform poorly at longer timescales too. The use of PPEs also provides additional information for model development, by identifying parameters and processes responsible for model errors at the two different timescales, and systematic errors that cannot be removed by any combination of parameter values.
The numerical weather prediction (NWP) of fog remains a challenge, with accurate forecasts relying on the representation of many interacting physical processes. The recent Local And Non‐local Fog ...EXperiment (LANFEX) has generated a detailed observational dataset, creating a unique opportunity to assess the NWP of fog events. We evaluate the performance of operational and research configurations of the Met Office Unified Model (MetUM) with three horizontal grid lengths, 1.5 km and 333 and 100 m, in simulating four LANFEX case studies. In general, the subkilometre (sub‐km) scale versions of MetUM are in better agreement with the observations; however, there are a number of systematic model deficiencies. MetUM produces valleys that are too warm and hills that are too cold, leading to valleys that do not have enough fog and hills that have too much. A large sensitivity to soil temperature was identified from a set of parametrisation sensitivity experiments. In all the case studies, the model erroneously transfers heat too readily through the soil to the surface, preventing fog formation. Sensitivity tests show that the specification of the soil thermal conductivity parametrisation can lead to up to a 5‐hr change in fog onset time. Overall, the sub‐km models demonstrate promise, but they have a high sensitivity to surface properties.
Fog remains a challenge to forecast accurately using numerical weather prediction. We evaluate the performance of the Met Office Unified Model (MetUM) at both kilometre and subkilometre (sub‐km) grid lengths. MetUM produces valleys that are too warm and hills that are too cold, leading to valleys that do not have enough fog and hills that have too much. The sub‐km scale configurations generally outperform the km scale, but they are highly sensitive to the soil thermal conductivity.
In this paper we define the first Regional Atmosphere and Land (RAL) science configuration for kilometre-scale modelling using the Unified Model (UM) as the basis for the atmosphere and the Joint UK ...Land
Environment Simulator (JULES) for the land. RAL1 defines the science configuration of the dynamics and physics schemes of the atmosphere and
land. This configuration will provide a model baseline for any future weather or climate model developments to be described
against, and it is the intention that from this point forward significant changes to the system will be
documented in the literature. This reproduces the process used for global configurations of the UM,
which was first documented as a science configuration in 2011. While it is our goal to have a single
defined configuration of the model that performs effectively in all regions, this has not yet been
possible. Currently we define two sub-releases, one for mid-latitudes (RAL1-M) and one for tropical
regions (RAL1-T). The differences between RAL1-M and RAL1-T are documented, and where appropriate we define how the model
configuration relates to the corresponding configuration of the global forecasting model.
We present significant differences in the simulated atmospheric flow for warm, tidally locked small Neptunes and super Earths (based on a nominal GJ 1214b) when solving the simplified, and commonly ...used, primitive dynamical equations or the full Navier-Stokes equations. The dominant prograde, superrotating zonal jet is markedly different between the simulations, which are performed using practically identical numerical setups, within the same model. The differences arise due to the breakdown of the so-called "shallow-fluid" and traditional approximations, which worsens when rotation rates are slowed, and day-night temperature contrasts are increased. The changes in the zonal advection between simulations solving the full and simplified equations, give rise to significant differences in the atmospheric redistribution of heat, altering the position of the hottest part of the atmosphere and temperature contrast between the daysides and nightsides. The implications for the atmospheric chemistry, and therefore, observations need to be studied with a model including a more detailed treatment of the radiative transfer and chemistry. Small Neptunes and super Earths are extremely abundant and important, potentially bridging the structural properties (mass, radius, and composition) of terrestrial and gas giant planets. Our results indicate care is required when interpreting the output of models solving the primitive equations of motion for such planets.
Aerosol processes and, in particular, aerosol‐cloud interactions cut across the traditional physical‐Earth system boundary of coupled Earth system models and remain one of the key uncertainties in ...estimating anthropogenic radiative forcing of climate. Here we calculate the historical aerosol effective radiative forcing (ERF) in the HadGEM3‐GA7 climate model in order to assess the suitability of this model for inclusion in the UK Earth system model, UKESM1. The aerosol ERF, calculated for the year 2000 relative to 1850, is large and negative in the standard GA7 model leading to an unrealistic negative total anthropogenic forcing over the twentieth century. We show how underlying assumptions and missing processes in both the physical model and aerosol parameterizations lead to this large aerosol ERF. A number of model improvements are investigated to assess their impact on the aerosol ERF. These include an improved representation of cloud droplet spectral dispersion, updates to the aerosol activation scheme, and black carbon optical properties. One of the largest contributors to the aerosol forcing uncertainty is insufficient knowledge of the preindustrial aerosol climate. We evaluate the contribution of uncertainties in the natural marine emissions of dimethyl sulfide and organic aerosol to the ERF. The combination of model improvements derived from these studies weakens the aerosol ERF by up to 50% of the original value and leads to a total anthropogenic historical forcing more in line with assessed values.
Key Points
The HadGEM3‐GA7 climate model has a large, negative aerosol ERF resulting in an unrealistic negative total anthropogenic forcing of climate
The aerosol ERF is shown to be highly sensitive to the underlying physical and aerosol model processes and parameterizations
Through a combination of scientific model improvements the aerosol ERF is reduced by up to 50% from –2.75 to –1.45 W/m2