Mineral dust plays an important role in Earth system models and is linked to many components, including atmospheric wind speed, precipitation and radiation, surface vegetation cover and soil properties and oceanic biogeochemical systems. In this paper, the dust scheme in the first configuration of the United Kingdom Earth System Model UKESM1 is described, and simulations of dust and its radiative effects are presented and compared with results from the parallel coupled atmosphere–ocean general circulation model (GCM) HadGEM3-GC3.1. Not only changes in the driving model fields but also changes in the dust size distribution are shown to lead to considerable differences to the present-day dust simulations and to projected future changes. UKESM1 simulations produce a present-day, top-of-the-atmosphere (ToA) dust direct radiative effect (DRE – defined as the change in downward net flux directly due to the presence of dust) of 0.086 W m
Dust is an important component of the Earth system that is linked to the atmosphere, surface, aerosols and terrestrial and oceanic biospheres (Carslaw et al., 2010). Many factors control the dust life cycle, including vegetation cover, soil moisture, wind speed and precipitation; hence, dust is potentially sensitive to changes in climate through multiple mechanisms. Dust, in turn, affects climate via radiative processes in the atmosphere, via interaction with clouds and other aerosol species, by changing surface albedo and as a source of nutrients for oceanic and terrestrial biogeochemical cycles.
Climate models have long included dust schemes (e.g. Tegen and Miller, 1998; Woodward, 2001; Miller et al., 2006; Kok et al., 2014), which typically allowed dust emission, transport and deposition to be modelled together with the direct radiative effects of atmospheric dust. Some attempts were made to include the impact of the vegetation response on dust (Tegen et al., 2002; Mahowald et al., 2006), though, as the vegetation was not simulated interactively, feedbacks between vegetation and climate were excluded. The advent of Earth system models (Collins et al., 2011; Watanabe et al., 2011), which include interactive vegetation and ocean biogeochemical schemes together with other components, has extended the possibilities for the investigation of interactions between dust and other systems. This also allows for the impact of the responses of Earth system components on the dust to be included in future climate projections.
Here we employ the state-of-the-art UKESM1 model (Sellar et al., 2019) to explore the behaviour of dust within the Earth system in the present day and under possible future development pathways. UKESM1 is built around the HadGEM3-GC3.1 coupled ocean–atmosphere model (Kulbrodt et al., 2018; Williams et al., 2017), and we have made use of this to run parallel experiments with the coupled and Earth system models, in order to investigate the importance of Earth system processes to the dust simulation. Both models use the same dust emission scheme, though the tuning parameters were adjusted for UKESM1. We include particle sizes up to 63
In this paper, we give a brief introduction to UKESM1 and describe the dust scheme and the experiments used before assessing the present-day simulation of dust load, size distribution, aerosol optical depth (AOD) and deposition against observations, and we also present the simulated dust direct radiative effect (DRE) results. We then compare the UKESM1 simulations with parallel results from HadGEM3-GC3.1 and investigate the causes of the differences. In particular, we examine the consequences of including interactive vegetation and other Earth system processes and the effects of the changes to the dust scheme settings, including the impacts of retuning on load and size distribution and of excluding emissions from seasonal sources. In the next section, we investigate the response of dust to potential climate change as driven by two Shared Socioeconomic Pathways and consider the mechanisms involved. Finally, we present a discussion of the results and summarise our main findings.
UKESM1 is the latest-generation UK Earth system model built around the HadGEM3-GC3.1 coupled atmosphere–ocean general circulation model (Kulhbrodt et al., 2018; Williams et al., 2017) combined with the Model of Ecosystem Dynamics, nutrient Utilisation, Sequestration and Acidification (MEDUSA) biogeochemical model (Yool et al., 2013). The interface between these is provided by the OASIS coupler (Craig et al., 2017). The component models within HadGEM3 are the Met Office Unified Model (UM) atmospheric model (Williams et al., 2017), which contains the United Kingdom Chemistry and Aerosols (UKCA) stratospheric–tropospheric chemistry model (Archibald et al., 2020) with the GLOMAP-mode version of the Global Model of Aerosol Processes scheme (Mulcahy et al., 2020), together with the Joint UK Land Environment Simulator (JULES) model (Walters et al., 2019), the Nucleus for European Modelling of the Ocean (NEMO) model (Storkey et al., 2018) and the CICE sea ice model (Ridley et al., 2018). In UKESM1, these components are used for the interactive simulation of Earth system processes, such as the full atmospheric chemistry from UKCA. The modelling of vegetation and surface properties is particularly important for dust. In UKESM1 the TRIFFID (Top-down Representation of Interactive Foliage and Flora Including Dynamics) scheme within JULES simulates interactive vegetation, while, in HadGEM3-GC3.1, data from the International Geosphere–Biosphere Programme (IGBP) climatology (IGBP, 2000) are used. The soil parameters also differ. Those in HadGEM3-GC3.1 are based on Van Genuchten (Loveland et al., 2000), while those in UKESM1 are from Brooks and Corey (1964). A full description of UKESM1 is available in Sellar et al. (2019).
The mineral dust is simulated within HadGEM3 by the fully interactive dust scheme described below, which is called each atmospheric model time step. The driving fields are calculated directly by the UM and JULES, and dust impacts the rest of the model through radiative interactions with the UM atmosphere and through input to the ocean biogeochemistry in MEDUSA. In the current configuration, the dust is externally mixed with other aerosols.
The dust scheme is a development of that described in Woodward (2001), with
significant improvements to the emission scheme and newer refractive index
data. Dust emission in six size bins, with boundaries at 0.06324, 0.2,
0.6324, 2.0, 6.324, 20.0 and 63.24
Flux calculations are based on the method of Marticorena and Bergametti (1995). Horizonal flux in bin
In UKESM1 dust is emitted only from the bare soil fraction of a grid box, though there is also an option in the code to allow emissions from seasonally bare sources, based on the leaf area index (LAI). Dust emissions are prevented if snow is present, if the ground is frozen, on steep slopes, if soil moisture exceeds a threshold (see below), and at coastal points where the lowest-level wind speed over land may be anomalously high. No preferential source terms are used.
The driving fields for the scheme are the model's grid box mean time step
mean fields, but Eq. (1) was derived from shorter timescale measurements at single locations. Corrections are therefore needed to account for the effect of spatial and temporal averaging. Here, model friction velocity
Dry threshold friction velocity (
The Fécan et al. (1999) treatment of soil moisture was designed for use in arid or semi-arid areas, and in order to apply it to the whole Earth, a further constraint on dust production from moist soil is required. Dust production is inhibited when soil moisture exceeds a particle-size-dependent threshold (
The scheme described here uses the same code as the UM global atmosphere
GA7.1 and global coupled GC3.1 configurations but with the following two changes in the settings: (1) the three tuning terms (
The model provides a self-consistent set of diagnostics of dust emissions,
depositions and mixing ratios (concentrations). While the concentration,
wet deposition and dry deposition over water are directly comparable to
results from other models and to observations, the emission and dry
deposition in source regions are not. Atmospheric lifetime, being derived
from these terms, is also not comparable. The emission and dry deposition
diagnostics include all particles released from the surface, whether or not
they are added to the atmospheric load or interact with the model atmosphere
in any way. In the UM the flux of particles released from the surface is calculated first, and then another section of code performs both the mixing into the atmosphere and the dry deposition back to the surface in a single loop. As a result, a fraction of the particles released from the surface is immediately dry deposited and never added to the atmospheric burden. In physical terms the redeposited fraction may be considered to be those particles which fall back to the surface within the model time step. The magnitude of this effect is greater for larger particles due to the size-dependence of sedimentation and, because the particle size range extends to 63
The parameters
The aim of the tuning process was to obtain the optimum simulation of emissions by correcting for terms and processes not included in the climate model but which impact dust emissions, such as the effect of gustiness, the relation between moisture in the model's top soil level and at the soil surface, etc. In the absence of emission measurements, observations of other dust properties at some distance from the sources must be used for model evaluation. These are affected not only by the emissions but also by transport and deposition processes and, in some cases, by the dust radiative properties. Thus the tuning is inevitably influenced by, and compensates for, any biases in these processes and properties in the dust scheme, in addition to biases in the driving model fields, such as in bare soil fraction, which affected the dust. In an attempt to reduce the effects of these biases, higher importance was given to comparisons with direct observations of dust concentration and size than with observations of properties such as AOD and with observations nearer sources than with observations at remote locations. Correct simulation of near-source concentration and size distribution is also important because the impact on the dust radiative forcing and potential feedbacks will be largest in source areas.
Initial sets of tuning experiments were used to explore the parameter space, finding how the model responded to changes in each of the variables and constraining their range; later experiments focussed on the finer tuning of the variables. Due to time constraints in the development of UKESM1, dust scheme tuning was carried out in parallel with the later stages of the Earth system model (ESM) development. As a consequence of this, tuning terms had to be adjusted on an ongoing basis in response to changes in the driving model, and there was insufficient time for a highly detailed tuning process. Though otherwise not ideal, this parallel development had the advantage that minor changes could be made to other parts of the ESM to improve the dust simulation. In particular, the lai_min term for grasses in JULES was reduced to 0.3, allowing the grasses to spread more readily and improving the simulation of bare soil fraction (Sellar et al., 2019).
The tuning of UKESM1 dust was carried out primarily with UKESM1-CN, the offline chemistry version of UKESM1 (Sellar et al., 2019), as this ran more quickly than the full ESM, and the interactive chemistry was expected to have little impact on the dust. Comparisons of dust simulations from the two versions of the ESM showed only small differences, with the various tuning settings tending to perform similarly in each case. Preindustrial (PI) simulations were used for the tuning as a fully spun-up present-day experiment was not available at that stage. While this meant that the model climate did not match the period of the observations, this discrepancy was likely to have had a relatively small effect compared with the cumulative effect of the many uncertainties in the dust scheme, such as biases in the simulated bare soil fraction, soil moisture and wind speed, resolution-related limitations in capturing some emissions mechanisms and the omission of surface-crusting effects and other poorly understood processes. Tuning tests were also carried out with “pseudo-present-day” experiments, initialised from PI runs but then run to equilibrium with present-day (PD) atmosphere and land use settings. These produced similar results to the main tuning experiments and gave confidence in the use of PI simulations to establish tuning parameters. However, the climate in the fully spun-up present-day UKESM1 simulation reported here was slightly different from that in the tuning runs.
The dust in HadGEM3-GC3.1 was not retuned for this work. HadGEM3 is a widely used model (Williams et al., 2018) which has been regularly updated, with some dust retuning being carried out before the release of a new version if the update had a major effect on the dust. Here we make use of the existing HadGEM3-GC3.1 configuration (Kuhlbrodt et al., 2018), as used for the CMIP6 simulations (Eyring et al., 2016), without any changes, except in experiments H3_TUK_EXSS and H3_TUK_INSS, as described in Sect. 4. The dust tuning parameters had been set to
The results presented here were primarily obtained from UKESM1 and HadGEM3-GC3.1 historical experiments parallel to those performed as part of CMIP6 (Eyring et al., 2016) but with extra diagnostics. In order to obtain the diagnostics needed for calculating DRE, a copy of each relevant CMIP6 experiment was run for 20 years with the “double call” method. This involved calling the radiation scheme twice each time step, with the radiative effect of dust excluded from the first call but included in the second which is used to progress the model. The dust DREs are then calculated as the differences between the fluxes from each call. In addition, two parallel HadGEM3-GC3.1 experiments were run. In both of these, the tuning terms were set to UKESM1 values, and in one the seasonal sources were also deactivated. These allowed us to investigate the relative importance of the various differences between the main simulations. Two UKESM1 Atmospheric Model Intercomparison Project (AMIP) experiments from CMIP6 were also used. The experiments are summarised in Table 1.
Summary of experiments. The UKESM1 tuning settings for
Figure 1 shows dust load from the UKESM1 and HadGEM3-GC3.1 present-day
simulations (UK_PD and H3_PD), together with the differences and fractional differences. Both models capture the expected global dust distributions qualitatively well. The UK_PD global load of 19.5 Tg is 30 % higher than the 15.0 Tg H3_PD value. These values span the AeroCom phase III mean of 16.6 Tg and are within the range of 5.7 to 22.3 Tg shown by participating models (Gliß et al., 2021). They are also consistent with the more observationally constrained estimate by Kok et al. (2017) of a PM
Dust load (Tg) from H3_PD, UK_PD, difference and fractional difference.
Scatterplots of model versus observational multi-annual seasonal
mean concentrations at remote marine locations and AODs in dusty areas for
Figure 2 shows scatterplots of seasonal mean modelled dust concentrations versus observations from remote sites from the University of Miami network and AERONET sites (Holben et al., 2001) chosen for data availability and dustiness, based on the Ångström exponent (see Sect. 5.3). At least 4 years of monthly data were available for the chosen sites. Both UK_PD and H3_PD show good agreement with the observations, with the models being within a factor of 3 of the observations at most stations for most seasons. The correlation coefficients for the concentration and AOD data shown in Fig. 2 are 0.89 for UKESM1 and 0.87 for HadGEM3-GC3.1. The slight high bias at a few North Pacific stations in spring and summer is increased in UKESM1 due to larger bare soil areas in Asia simulated by TRIFFID, with associated increased wind speed and reduced soil moisture. Similarly, increased bare soil in Australia results in an increased high bias at the nearest South Pacific station. The low biases almost all occur at two stations on the Antarctic Peninsula, where the dust from Patagonia dominates. The concentrations here are likely to be very sensitive to the westerly winds in that region. The bias is worsened in UKESM1 due to a low bias in the simulated bare soil fraction. Overall, the UKESM1 surface concentrations show good agreement with observations, compared with other ESMs (Checa-Garcia et al., 2021). The similarity in the level of performance between UKESM1 and HadGEM3-GC3.1 is noteworthy, given the many extra processes and feedbacks within the ESM and, in particular, its use of interactive vegetation.
UKESM1 dust was tuned to give better agreement with the size distribution data from the FENNEC campaign (Ryder et al., 2013) than was shown in HadGEM3-GC3.1. These data were particularly useful as they provided measurements near sources where the effects of deposition and transport on the size distribution would be relatively small and, unusually, included data from several instruments for measuring larger particles with diameters above 3
Dust volume size distributions, including the
Figure 3 presents normalised volume size distributions from 20-year June means from UK_PD and H3_PD compared with a fit to the FENNEC observations. UK_PD shows good agreement with the observations throughout the size range. The lower-level data have a somewhat greater coarse particle fraction than the higher-level data, as might be expected in a source region, given the short lifetimes of the largest particles. The peak diameter is slightly smaller than in the observations. This would be consistent with the 20-year June mean containing a smaller fraction of freshly emitted particles than were measured in observations which were intended to sample dusty conditions.
H3_PD concentrations overlap with the FENNEC observations, though the overall agreement is poorer, as the model is outside the range of the observations for the two smallest bins, and the peak of the distribution is in bin 4 (2–6
Comparison of monthly mean model AODs with data from AERONET (Holben et
al., 2001) sites in dusty areas (Fig. 4) shows that the modest low bias in
H3_PD is slightly worsened in UK_PD, though the model mean is within 2 standard deviations of the observations in almost all cases. The stations were selected to obtain as realistic observational climate means as possible, given the considerable variability in dust and the short data record, and to minimise the effects of other aerosols. Sites were chosen from those in potentially dusty areas as having a minimum of 4 years of monthly data, with at least 10 daily means per month, and a monthly mean Ångström exponent (870–440) below 0.5 for at least 10 months of the year. The only sites to fulfil these criteria were Tamanrasset_INM in the Sahara, Cabo Verde off the West African coast and six Sahelian stations all in a narrow band between 12 and 16
Annual cycle of AOD (440 nm) at dust-dominated AERONET stations from UK_PD and H3_PD. Site locations are Agoufou (1.48
The observed multi-annual mean AOD at 440 nm averaged over these AERONET sites is 0.45, while UK_PD simulates 0.28 and H3_PD simulates 0.37. Comparison of the simulations shows that the lower AOD in UK_PD is due to a combination of lower dust optical depth (DOD; 0.13 in UK_PD and 0.19 in H3_PD) with lower optical depth due to other species (0.15 in UK_PD and
0.18 in H3_PD). The lower DOD in UK_PD is caused by the larger particle sizes in that model, as evidenced by the Ångström exponent for the DOD (440–870), which is
Studies by Marticorena et al. (2010, 2017) have shown that the annual cycle
of dust in the western Sahel is related to the timing of the West African
Monsoon (WAM) and the annual north–south shift in the Intertropical Convergence Zone (ITCZ). Through winter and early spring, the Harmattan flow brings dust from the Sahara to this region, and the increasing AODs over the first few months of the year are consistent with the annual cycle of Saharan dust, as seen at Tamanrasset INM. From May to October the WAM dominates the region, bringing precipitation and wind from the southwest. During the early part of the monsoon season, strong but sporadic local emission events driven by mesoscale convective systems (MCSs; Caton-Harrison et al., 2019, 2020) produce similar levels of dust to those seen in the dry season. The UK_PD simulations slightly underestimate the AODs due to Saharan dust in the dry season, but the main sources of error appear to be in the timing of the arrival of the WAM, with AODs already decreasing in April, and in particular a failure to simulate the dust production from local sources. These failings are unsurprising, given the difficulties in simulating the WAM and MCSs in an N96 climate model with a resolution of 1.85
AOD (550 nm) from
Figure 5, which compares the total aerosol optical depth at 550 nm in UK_PD with MODIS data, shows a low bias in the Sahel, as mentioned above, and also in dust from the Bodélé Depression. The Bodélé dust source is very difficult to simulate in global climate models because a resolution of a few tens of kilometres is needed to represent the Bodélé low-level jet, which is responsible for much of the dust emissions (Todd et al., 2007). The dust-dominated AOD over the northern Sahara and the Arabian Peninsula is also low, while over Australia it is too high. The high bias on the southwestern side of the Himalayas and low bias on the northeastern side either suggest that the model may be failing to transport aerosol over the steep orography there or may be associated with other aerosol species. These biases reflect model weaknesses and also the difficulty of finding a single set of tuning terms which gives a good performance against all metrics. The settings which gave improvement in AOD over North Africa resulted in the worsening of biases associated with Asian and Australian dust.
The effects of non-sphericity on the dust particle optical properties were not included in the models. It has been estimated that the extinction efficiency could increase by as much as 29 % (Kok et al., 2017). An indication of the potential effect this might have was estimated by multiplying the UK_PD dust optical depth by a factor of 1.29, which resulted in an improved agreement with the satellite data (Fig. 5).
Though the model deposition diagnostics will not represent observed deposition in source regions (see above), they are comparable to measurements in locations where there are no local emissions. Deposition rates have been evaluated against annual mean measured fluxes from a range of data sources at sites remote from emission areas, as selected by Huneeus et al. (2011), for the evaluation of dust in AeroCom models. These include deposition observations from Ginoux et al. (2001), dust and iron deposition fluxes reported in Mahowald et al. (2009), ice core data from Mahowald et al. (1999) and sediment trap data from the DIRTMAP database (Tegen et al., 2002; Kohfeld and Harrison, 2001). Though this collection of deposition data from 84 sites is one of the most comprehensive available, there are considerable uncertainties associated with the observations, most notably that, although sites with less than 50 d of data were excluded, some of the records are not long enough to be considered climatological (Huneeus et al., 2011). This is particularly problematic for dust fields which show very strong variability, as a large fraction of annual deposition may occur over a few days of the year (Prospero et al., 2010).
Scatterplot of deposition rates (g m
A scatterplot of UK_PD mean deposition rates against these observations (Fig. 6) shows that the model agrees reasonably well with the observations, and model results are within a factor of 10 of the observations at most locations. The only area with a noticeable bias is in Antarctica and the Southern Ocean. The four Antarctic stations where the model significantly overestimates deposition rates are very close together near the Dumont station, between 64.60–64.97
Figure 7 and Table 2 show the dust direct radiative effect (DRE), which is the change in flux directly due to the presence of dust, excluding any secondary effects. This is calculated using the “double call” method, as described in Sect. 4.
Dust direct radiative effects (W m
Global mean load and all-sky direct radiative effects, due to dust in a present-day climate, from various experiments. Alternate rows show experimental results and the absolute differences, followed by percentage differences (in parentheses) between the experiment below and the experiment above.
In UKESM1 dust is generally more reflective than the surface in the
shortwave (SW), except over ice and the brightest deserts, with the result
that the net downward shortwave at the top of the atmosphere (ToA) is reduced everywhere but over these very light surfaces. The global mean ToA SW dust DRE is
In HaGEM3-GC3.1 the global mean net DRE at ToA is
The relatively small global mean net DREs are the residuals of the partial cancellation of areas of larger positive and negative net DREs, in addition to the partial cancellation of component DREs of different sign, each of which are sensitive to changes in the dust load, spatial distribution, size distribution and radiative properties. As a result, there are large uncertainties associated with estimates of global mean dust DREs.
Running simulations with almost identical dust schemes in two similar models allows us to explore the impacts of a number of factors on the dust simulation. The causes of the differences between the HadGEM3-GC3.1 and UKESM1 present-day dust load and DREs may be divided into the following four groups: (1) the change in size distribution with retuning, (2) the change in total load due to retuning, (3) the (de-)activation of seasonal sources and (4) the change in driving model. These are investigated using results from H3_PD, UK_PD and two additional experiments representing intermediate configurations. H3_TUK_EXSS is a HadGEM3-GC3.1 experiment with UKESM1 dust tuning and excluding seasonal sources. Comparison of this with UK_PD helps to isolate the effects of the change in model. H3_TUK_INSS is parallel to this but includes seasonal sources to allow investigation of the effect of these sources (see Table 1). Another set of dust DRE results H3_TUK_INSS_NL has been generated by normalising the DREs from H3_TUK_INSS by the ratio of H3_PD load to H3_TUK_INSS load. This allows us to separate the effects of the change in size distribution (from comparing H3_TUK_INSS_NL and H3_PD) from the effects of the load change (from comparing H3_TUK_INSS_NL and H3_TUK_INSS). Global mean volume size distributions from UK_PD, H3_PD and H3_TUK_INSS_NL are shown in Fig. 3c. The implicit assumption that the DREs are proportional to total load for a given size distribution is a reasonable first approximation, although it ignores the feedback between dust DREs and emissions (Miller et al., 2004; Woodage and Woodward, 2014; Kok et al., 2018). Load and DREs from each experiment, together with the differences due to each factor, are listed in Table 2.
Soil moisture in the top model layer, dust source area fraction, 10 m wind speed, dust emission diagnostic and dust load from H3_TUK_EXSS and UK_PD and the differences and fractional differences between experiments.
Overall, the changes from H3_PD to UK_PD result in an increase in total load of 30 %. The load change due to retuning (factor 2) is responsible for a reduction of 31 %, while the disabling of seasonal sources (factor 3) only decreases global load by a further 9 %. The change in driving model (factor 4) more than doubles the load and produces much the largest individual impact on the atmospheric burden. The most important element in this change is the vegetation. The disparity between the bare soil simulated by TRIFFID in UKESM1 and the IGBP climatology used in HadGEM3-GC3.1 accounts for 70 % of the extra emissions in UK_PD compared with H3_TUK_EXSS. The global average bare soil fraction is 0.26 in UK_PD and 0.24 in H3_UK_EXSS, but regional variations are much larger, and their geographic distribution promotes extra dust production in UK_PD. The 29 % of the UK_PD bare soil area that is vegetated in H3_TUK_EXSS is mostly in arid areas at the margins of existing deserts where conditions favour dust emission, while the 25 % of the H3_TUK_EXSS bare soil area that is vegetated in UK_PD is mainly in regions where moisture limits dust emission (Fig. 8g). Areas where there is less vegetation, and thus reduced roughness, might be expected to be associated with higher near-surface wind speed and hence increased evaporation and reduced soil moisture, providing conditions particularly favourable to dust production. This can be seen in the regions of greatest difference between UK_PD and H3_TUK_EXSS such as the Sahel, India and the Kazakh Steppe (Fig. 8g, d, k and o).
Global mean dust net DREs in UKESM1 are smaller than in HadGEM3_GC3.1, despite the larger load (Table 2, Fig. 7). The net surface DRE of
The global mean preindustrial to present-day dust forcing (defined here as
the change in DRE) is small in both models. It is calculated as the
difference between dust DREs in the 20-year mean UK_PD (and H3_PD) runs and a 20-year mean from the respective PI control runs. UKESM1 (HadGEM3-GC3.1) simulates a forcing of
The role of the vegetation response in the historical dust changes is investigated using the results of atmosphere-only AMIP experiments in which the driving fields, except vegetation, were taken from a UKESM1 preindustrial simulation. In the control experiment (A_UK_PI), the vegetation fields – LAI, canopy height and vegetation fraction – were also taken from the preindustrial simulation; in the parallel A_UK_PI_V2014NOLU experiment, the 2014 vegetation values from a UKESM1 historical experiment which excluded anthropogenic land use changes were used. In both experiments, the dust settings were as in UKESM1. The difference between these simulations gives an estimate of the effect of the change in vegetation alone on dust. The effect of climate change excluding the vegetation response (but including land use change) is estimated from the difference between HadGEM3-GC3.1 simulations H3_PD and H3_PI.
Difference in dust load (mg m
Figure 9 shows the differences in dust load, due to vegetation response and
to climate estimated in this way, together with the sum of these two changes, and the equivalent change in UKESM1 from the preindustrial UK_PI to the present-day UK_PD. The similarities between the patterns of load change due to the combined vegetation and climate changes and due to the changes in UKESM1 is notable, giving confidence that dust changes simulated by the models are comparable even though the present-day dust is somewhat different. From these experiments the global totals of the PD to PI differences in dust load are
The dust changes due to possible future changes in climate are explored using a set of scenario experiments representing future Shared Socioeconomic Pathways (Sellar et al., 2020; Riahi et al., 2017). Means of the last 20 years (2081–2100) of SSP5–8.5 and SSP2–4.5 experiments UK_S5, H3_S5, UK_S2 and H3_S2 are compared with the present-day (1995–2014) period of the historical experiments UK_PD and H3_PD. The Shared Socioeconomic Pathway SSP5–8.5 represents high-end projections of fossil fuel and energy use, food demand and greenhouse gas emissions, assuming fossil-fuelled development; it has a radiative forcing pathway similar to the highest Representative Concentration Pathway, RCP8.5 (Moss et al., 2010; van Vuuren et al., 2011). SSP2–4.5 represents a middle-of-the-road future, with social, economic and technological developments broadly following historical patterns, giving a radiative forcing pathway similar to RCP4.5. The scenario experiments are initialised from the end of the respective historical runs.
Dust source area fraction, atmospheric load and DRE from simulations of present-day and future climates. (Note that DRE diagnostics were not available for the H3_S2 and UK_S2 experiments.)
Results are summarised in Table 3 and Figs. 11 and 12. UKESM1 simulates a reduction of 23 % in the total dust load from the present day to the end of the SSP5 scenario, with the greatest reductions from the major deserts and only small areas experiencing increases (Fig. 11). In contrast, the global load in HadGEM3-GC3.1 decreases by only 4 %, with smaller local changes than in UKESM1 and reduced emissions from some deserts being almost completely balanced by increases elsewhere (Fig. 12). The equivalent results for SSP2–4.5 are reductions of 19 % for UKESM1 and 7 % for HadGEM3-GC3.1. In each model the pattern of load changes from the present to the end of SSP2–4.5 has a very similar geographical distribution to the equivalent SSP5–8.5 results. This is also true of the main drivers of dust emissions, namely soil moisture, source areas and wind speed. The slightly larger decrease in dust load in H3_S2 compared with H3_S5 is due to minor differences in the residuals from the cancellation of areas of positive and negative change. Given the similarities in the patterns of change in dust and its drivers in both pathways, the following analysis will focus only on the SSP5–8.5 experiments.
Grass fraction, tree and shrub fraction, bare soil fraction and dust source area fraction from UK_PD, UK_SSP5, H3_PD and H3_SSP5. (The H3_PD and H3_SSP5 source areas include seasonal sources.)
Soil moisture in the top layer, dust source area fraction, 10 m wind speed, dust emission diagnostic and dust load from UK_PD and UK_SSP5, together with differences and fractional differences between experiments.
As in Fig. 11 but for H3_PD and H3_SSP5.
Comparison of the changes in the drivers of dust emissions helps reveal the
causes of the differences between the dust changes in the two models. Source
area changes and associated vegetation changes are shown in Fig. 10. In UKESM1, the dynamic vegetation responds to the warming climate and enhanced
CO
Figures 11 and 12 show the changes in dust fields and the drivers of emissions, simulated by UKESM1 and HadGEM3-GC3.1 respectively. The top-level soil moisture in UKESM1 is seen to be reduced in all areas except deserts, giving a global mean loss of 12 %. A similar reduction of 15 % occurs in HadGEM3-GC3.1, suggesting that this effect is mostly a response to the changing climate, mediated by evaporation and somewhat mitigated by the reduction in bare soil fraction in the Earth system model.
In UKESM1 the lowest-level wind speed over land is reduced in northern mid- and high-latitudes, and also in Australia and southern South America, but is increased in parts of the tropics, particularly tropical South America, though not in desert regions. HadGEM3-GC3.1 exhibits similar increases in the tropics, but the mid- and high-latitude reductions are much smaller. This indicates that the tropical increases in wind speed are likely to be primarily a response to the changing climate, while the reductions likely are due to increased roughness produced by vegetation growth in the corresponding areas in UKESM1.
The increase in dust emissions at high latitudes in HadGEM3 is driven predominantly by the reduction in soil moisture and at low latitudes by the increases in wind speed, while the reduced emissions from arid areas are caused by the slightly reduced wind speeds and moister soil in those regions. These processes also occur in UKESM1, but the vegetation-driven loss of bare soil and larger reduction in wind speed have a greater impact, resulting in an enhanced reduction in emissions, with emissions increases occurring in only a few small regions. While there is a global near balance of dust load increases and decreases in HadGEM3-GC3, UKESM1 simulates a global load reduction of 23 % by the end of SSP5–8.5.
The difference in dust responses between the models is much larger than the difference in dust responses between the pathways. HadGEM3-GC3.1 results suggest that the global dust burden dust will remain largely unchanged, whatever socioeconomic pathway is followed, with increases mostly in North and South America and Australia balanced by decreases mainly in North Africa, Asia and Europe. In contrast, in UKESM1 the addition of the extra Earth system process, and particularly the interactive vegetation, results in projections of reduced global total load, with reductions from most of the main desert regions and only a few small areas of increase in the tropics. The global load reduction of 23 %, associated with fossil-fuelled development, is somewhat larger than the 19 % reduction of the middle-of-the-road pathway. We note the non-linearity of this response to forcing, despite the linear response of the source area, and speculate that this may be due to the dust emission process involving non-linear dependence on various factors, which themselves may respond non-linearly to radiative forcing, and that dust feedbacks, which may enhance or limit emissions (e.g. Miller et al., 2004; Woodage and Woodward, 2014), could also introduce non-linearity. Our simulations suggest that the impact of the vegetation response on dust is larger than the direct impact of future climate change, and the differences due to including Earth system processes in the simulations are larger than the differences between pathways.
The reduction in load in UKESM1 results in a decrease in ToA DRE from
UK_PD to UK_SP5 of 44 %, though the absolute values are small (
The simulation of present-day dust in UKESM1 shows good agreement with most of the observations, though the AOD is somewhat low. The performance is generally comparable with that of HadGEM3_GC3.1, despite the extra complexity of the additional processes and feedbacks in the Earth system model. The tuning of the lai_min term for grasses in JULES, which limits the rate at which grasses can spread, was an important factor in this, as it improved the simulation of bare soil fraction in source regions. Some degradation in model performance is to be expected as a cost of the increased complexity, which is included to allow investigation of the feedback processes (Jones et al., 2011). Based on the evaluations performed, it was not found necessary to use preferential sources to add further restrictions to the dust emissions in UKESM1 compared to HadGEM3-GC3.1, so the dust was able to respond to climate-driven changes without any potentially unrealistic constraints. It is, of course, possible that other, as yet unevaluated, properties of the present-day simulations, such as frequency of occurrence of high dust optical depth, might be improved by the use of preferential sources.
Tuning of the dust scheme in the final UKESM1 configuration was limited by time constraints, as is often the case with the development of a model which includes dust. Further work might possibly have resulted in an alternative choice of settings. For example, in the UKESM1.1 model configuration dust was retuned to give better agreement with observations of AODs, though at the expense of agreement with other observations of size distribution and of concentrations in some areas (Mulcahy et al., 2022).
The low bias in AOD in the Saharan plume is the main weakness of the UKESM1 tuning settings, though as AOD depends on dust optical properties in addition to size distribution and concentrations, it is unclear quite how much of this bias is due to tuning. While optical properties are based on observations, they may well not be representative of all Saharan dust. In particular, dust from the Bodélé is made up predominantly of diatomite, which has different physical properties compared with dust from other Saharan sources (Todd et al., 2007). The size distribution was shown to agree well with observations.
The simulation of the Saharan plume is also affected by model resolution. Up to half of the dust emissions from West Africa may be associated with cold pool outflows from moist convection (haboobs), particularly in summer when this fraction is most important (Caton-Harrison et al., 2019, 2020), but the parameterised rather than explicit convection used in the N96 model renders the simulation of these events very difficult (Marsham et al., 2011; Knippertz and Todd, 2012), likely leading to underestimation of emissions in that area. Emissions from the Sahel are probably also underestimated for similar reasons. The omission of dust ageing through chemical processing has the effect of increasing lifetime and therefore increasing dust concentrations remotely from sources, and though this is not a very large effect, the inclusion of ageing could impact the choice of tuning settings, allowing higher Saharan AODs without too much overestimation of the concentrations in other areas. A higher present-day AOD over the Sahara would likely be associated with larger DREs in both the SW and LW, assuming that the size distribution was unchanged. These would still largely cancel, but the net DRE would probably be larger, as would the projected forcings for 2100.
The results demonstrate the crucial role of the particle size distribution
in the simulation of the 3D dust distribution and particularly the radiative
effects. Comparison of H3_PD with H3_TUK_INSS_NL (and H3_TUK_INSS) shows that quite small shifts in the size distribution towards coarser particles can change the balance of global mean SW and LW DREs so that they almost completely balance, giving net DREs close to zero. The particle size range is the same in all cases, with a maximum diameter of 63
The DRE results from UKESM1 are consistent with the estimates of Kok et al. (2017), who find that most models underestimate the size of atmospheric dust compared with measurements and use analysis of size and other observations to constrain model estimates of dust ToA DRE to a range between
The response of dust to climate change in each of the models is very different, predominantly due to the direct and indirect impacts of the vegetation changes in UKESM1, which are not included in HadGEM3-GC3.1. This suggests that, in the future, the main anthropogenic impact on dust may be via the change in vegetation consequent on fossil-fuel-driven emissions rather than through the changes in climate variables, while changes in land use have a smaller effect. The differences between the models are greater than the differences between the SSP2 and SSP5 simulations from either model. These results illustrate the importance of including Earth system processes when simulating the response of dust to climate change. Indeed, a realistic future dust simulation may not be possible without including vegetation changes and its effects.
Dust is particularly sensitive to the interactively simulated vegetation in
an Earth system model. Marginal areas for vegetation at the edges of deserts
may be difficult for a vegetation scheme to simulate realistically but are
potentially strong dust sources. Also, where vegetation is lost, not only do
source areas increase but the consequent reduced roughness may lead to increased wind speed and, hence, via increased evaporation to reduced soil moisture changes which tend to enhance dust emissions. However, this effect
is small in our climate change studies, as vegetation tends to grow rather
than die back in response to increased CO
A feature of the pattern of emission changes, particularly in UKESM1, is the location of areas of large increase very close to areas of large decrease, with a sharp boundary between them. This demonstrates the sensitivity of the dust to the balance of the various driving fields in addition to the high sensitivity to wind speed. Even modest biases in driving model fields could produce shifts in such balances, and the areas of sharp gradients in dust changes are unlikely to be captured correctly by any current schemes and present a challenge for future dust modelling. There are, inevitably, considerable uncertainties associated with such simulations as these. Biases in the driving model can have a significant impact on the dust simulation. Dust production is highly sensitive to friction velocity, but emission only occurs when this exceeds a threshold which depends on soil moisture; hence, biases in model wind speed, precipitation or soil hydrology can have a large impact on dust. Feedbacks, such as those between soil moisture, vegetation cover and wind speed, can exacerbate the effects of biases. Precipitation biases have been noted in both models, particularly along the ITCZ (Williams et al., 2017; Sellar et al., 2019), and these impact soil moisture. In UKESM1 the precipitation biases also affect vegetation, with a lack of monsoon rain resulting in lack of C3 grass and excessive bare soil in India and the Sahel (Sellar et al., 2019). The model resolution limits the realistic representation of some processes, such as mesoscale convective systems, which results in biases in precipitation and dust mobilisation in the Sahel. These biases could be exacerbated by dust–precipitation feedbacks, as demonstrated by Yoshioka et al. (2007). Biases in the bare soil fraction also impact emissions, an effect which will be larger in UKESM1 due to the use of the interactive vegetation scheme. Transport and deposition are also affected by driving model biases, particularly in the circulation and precipitation.
Simplifications, uncertainties and missing processes in the dust scheme will all produce biases and contribute to uncertainty in the results. Seasonal sources are excluded from the scheme because they depend on the distributions of various plant types, and it was not clear whether these could be simulated sufficiently accurately in the ESM. Comparison of experiments H3_TUK_INSS and H3_TUK_EXSS indicates that, without other changes, the omission of seasonal sources would have reduced the load by approximately 10 %, though we have not assessed the realism of dust production from such sources in HadGEM3-GC3.1. The emission scheme ignores some factors, including re-entrainment and the effects of surface crusting and surface geomorphology. The global tuning terms can only compensate for the effects of temporal and spatial averaging very approximately, and the impact of heterogeneous, short-term and sub-grid-scale phenomena such as gustiness cannot be well represented. There are many uncertainties associated with deposition, as the role of such factors as electrical charging, electrophoresis, diffusiophoresis and rear capture in below-cloud scavenging is not yet fully understood. The treatment of deposition here is relatively simple, and in-cloud scavenging is ignored as it is a relatively small effect compared to below-cloud scavenging for insoluble dust. Dust ageing through chemical processing is also ignored. The inclusion of ageing and in-cloud scavenging would be likely to result in lower concentrations remotely from sources. The use of a uniform set of refractive index data for all dust will result in biases in the DRE (Scanza et al., 2015). Lack of observational input data, such as global surface geomorphology, limits the representation of dust. There is also a paucity of data for the validation of dust on climate timescales.
The accuracy of the estimates of the response of dust to climate change is difficult to assess. Any biases, uncertainties and feedbacks in the present-day simulations could well change with changing climate and will inevitably impact the estimates of the dust responses and feedbacks. Factors associated with model resolution, driving fields, limitations of the dust scheme, missing processes or missing feedbacks could all be important. There is some evidence that a climate model may be unable to simulate the observed variability in dust over decadal timescales (Yoshioka et al., 2007; Mahowald et al., 2010), which would be indicative of missing feedbacks and could result in the underestimation of changes in future projections. We have not analysed the variability in the modelled dust compared with observations, though UKESM1 does exhibit higher dust variability than HadGEM3-GC3.1 probably due to the additional processes and feedbacks in the ESM.
Despite the inherent limitations and uncertainties, these results can provide much useful information, provided that the caveats are borne in mind. Further work is required to provide the missing observations, to study the processes affecting dust, and to improve dust schemes and the models driving them, in order to enhance understanding and improve the simulation of the effects of climate change on dust.
The dust scheme used in UKESM1 and HadGEM3-GC3.1 is a development of the Woodward (2001) scheme with an improved emissions parameterisation. It was initially used in HadGEM3 and then retuned for UKESM1. Seasonal sources were also deactivated in UKESM1, reducing dependence on multiple plant types in the interactive vegetation scheme, while having only a small effect on load (less than 10 % in HadGEM3-GC3.1). Evaluation of the UKESM1 present-day dust simulation showed good agreement with observations, comparable with that of HadGEM3-GC3.1. This is particularly encouraging, given the additional uncertainties produced by the extra processes and feedbacks within the Earth system model and the fact that no extra constraints, such as preferential sources, were applied to limit the dust emissions in UKESM1.
The differences between global mean UKESM1 and HadGEM3-GC3.1 present-day DREs have been shown to depend on the change in size distribution consequent on the retuning as much as on the change in driving model. The change in load due to retuning had a lesser impact, and the activation of seasonal sources produced only a small effect. This demonstrates the importance of the simulation of size distribution, which impacts the global dust concentrations through the size-dependent deposition processes and, additionally, impacts DREs through the size-dependent radiative properties. The magnitude of this effect will depend on the choice of refractive index data.
The response of the dust under future socioeconomic pathways is highly model dependent. In HadGEM3-GC3.1 the climate response produces drying of moist soils, which tends to increase emissions slightly, but this is balanced by a small reduction in wind speeds in major source regions leading to the total emissions and load remaining almost unchanged between the present day and 2100. While these processes also occur in UKESM1, a greater impact on dust comes from the vegetation response, which is simulated interactively in this model. Enhanced vegetation growth produces a decrease in bare soil, leading to further reductions in wind speed, both of which result in lower emissions and load. The differences between the models are greater than the differences between the SSP2–4.5 and SSP5–8.5 pathways, and though there are some considerable uncertainties associated with these results, they indicate the importance of including the vegetation response in projections of dust in future climates.
These results provide useful new information about the interactions between dust and climate. They highlight the need to represent the full dust size distribution as realistically as possible and indicate that the effect of doing this may be to reduce the global mean net dust DREs at ToA to a value close to zero (though impacts on climate may be larger, particularly on a regional scale). They also show the importance of including Earth system interactions in dust simulations over climate timescales, as the greatest driver of dust change may be vegetation changes and the consequent modifications to source areas, soil moisture and wind speed.
Information on the UKESM1 configuration and its components and the prerequisites for using it may be found at
The simulations used in this work are based on CMIP6 simulations which are archived on the Earth System Grid Federation (ESGF) node
SW was primarily responsible for the development of the dust code and its application in UKESM1 and ran the majority of the experiments and performed the analysis. AAS was scientific manager for UKESM1, coordinated the implementation of the dust code in the ESM and provided invaluable suggestions and comments throughout. YT implemented the dust settings into the UKESM1 model. MS wrote the code for dust deposition into the ocean. AY provided data and advice on deposition of dust to the ocean. ER ran the UK AMIP experiments and advised on their use. AW advised on the lai_min vegetation changes in UKESM1.
The contact author has declared that none of the authors has any competing interests.
Publisher’s note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
This study would not have been possible without the work of many staff in the UKESM1 Core Group and the Met Office, who developed the UKESM1 model and set up and ran the CMIP6 experiments. We thank Claire Ryder, for providing information about the FENNEC campaign. We are also grateful to all the researchers who generated the observational data we have used, including those of the University of Miami Aerosol Network and the AERONET community. We thank Didier Tanre, Jean-Louis Rajot, Philippe Goloub, Sidi Baika, Africa Barreto and Mahamadou M. Zakari, for establishing and maintaining the AERONET sites at Agoufou, Ouagadougou, Cabo Verde, Banizoumbou, Dakar, IER Cinzana, Tamanrasset INM and Zinder Airport.
This research has been supported by the UK Government Department for Business, Energy and Industrial Strategy and Department for Food, Environment and Rural Affairs (Met Office Hadley Centre Climate Programme). Marc Stringer has been supported by the Natural Environment Research Council (grant no. NE/N017978/1). Andrew Yool has been supported by EU Horizon 2020 project CRESCENDO (grant no. 641816).
This paper was edited by Susannah Burrows and reviewed by Ron L. Miller and one anonymous referee.