Robust projections and predictions of climate variability and change, particularly at regional scales, rely on the driving processes being represented with fidelity in model simulations. The role of enhanced horizontal resolution in improved process representation in all components of the climate system is of growing interest, particularly as some recent simulations suggest both the possibility of significant changes in large-scale aspects of circulation as well as improvements in small-scale processes and extremes.
However, such high-resolution global simulations at climate timescales, with
resolutions of at least 50 km in the atmosphere and 0.25
Increases in high-performance computing (HPC) resources, as well as the revised experimental design for CMIP6, now enable a detailed investigation of the impact of increased resolution up to synoptic weather scales on the simulated mean climate and its variability.
The High Resolution Model Intercomparison Project (HighResMIP) presented in this paper applies, for the first time, a multi-model approach to the systematic investigation of the impact of horizontal resolution. A coordinated set of experiments has been designed to assess both a standard and an enhanced horizontal-resolution simulation in the atmosphere and ocean. The set of HighResMIP experiments is divided into three tiers consisting of atmosphere-only and coupled runs and spanning the period 1950–2050, with the possibility of extending to 2100, together with some additional targeted experiments. This paper describes the experimental set-up of HighResMIP, the analysis plan, the connection with the other CMIP6 endorsed MIPs, as well as the DECK and CMIP6 historical simulations. HighResMIP thereby focuses on one of the CMIP6 broad questions, “what are the origins and consequences of systematic model biases?”, but we also discuss how it addresses the World Climate Research Program (WCRP) grand challenges.
Recent studies with global high-resolution climate models have demonstrated the added value of enhanced horizontal atmospheric resolution compared to the output from models in the CMIP3 and CMIP5 archive. They showed significant improvement in the simulation of aspects of the large-scale circulation such as El Niño–Southern Oscillation (ENSO) (Shaffrey et al., 2009; Masson et al., 2012), tropical instability waves (Roberts et al., 2009), the Gulf Stream (Kirtman et al., 2012), and Kuroshio (Ma et al., 2016), and their influence on the atmosphere (Minobe et al., 2008; Chassignet and Marshall, 2008; Kuwano-Yoshida et al., 2010; Small et al., 2014b; Ma et al., 2015), the global water cycle (Demory et al., 2014), snow cover (Kapnick and Delworth, 2013), the Atlantic inter-tropical convergence zone (ITCZ) (Doi et al., 2012), the jet stream (Lu et al., 2015; Sakaguchi et al., 2015), storm tracks (Hodges et al., 2011), and Euro–Atlantic blocking (Jung et al., 2012). High horizontal resolution in the atmosphere has a positive impact in representing the non-Gaussian probability distribution associated with the climatology of quasi-persistent low-frequency variability weather regimes (Dawson et al., 2012). In addition, the increased resolution enables a more realistic simulation of small-scale phenomena with potentially severe impacts such as tropical cyclones (Shaevitz et al., 2015; Zhao et al., 2009; Bengtsson et al., 2007; Murakami et al., 2015; Walsh et al., 2012; Ohfuchi et al., 2004; Bell et al., 2013; Strachan et al., 2013; Walsh et al., 2015), tropical–extratropical interactions (Baatsen et al., 2015; Haarsma et al., 2013), and polar lows (Zappa et al., 2014). Other phenomena that are sensitive to increasing resolution are ocean mixing, sea-ice dynamics, the diurnal precipitation cycle (Sato et al., 2009; Birch et al., 2014; Vellinga et al., 2016), quasi biennial oscillation (QBO) (Hertwig et al., 2015), the Madden–Julian oscillation (MJO) representation (Peatman et al., 2015), atmospheric low-level coastal jets and their impact on sea surface temperature (SST) bias along eastern boundary upwelling regions (Patricola and Chang, 2016; Zuidema et al., 2016), and monsoons (Sperber et al., 1994; Lal et al., 1997; Martin, 1999). The improved simulation of climate also results in better representation of extreme events such as heat waves, droughts (Van Haren et al., 2015), and floods. Enhanced horizontal resolution in ocean models can also have beneficial impacts on the simulations. Such impacts include improved simulation of boundary currents, Indonesian throughflow, and water exchange through narrow straits, coastal currents such as the Kuroshio, Leeuwin Current, and Eastern Australian Current, upwelling, oceanic eddies, SST fronts (Sakamoto et al., 2012; Delworth et al., 2012; Small et al., 2015), ENSO (Masumoto et al., 2004; Smith et al., 2000; Rackow et al., 2016), and sea-ice drift and deformation (Zhang et al., 1999; Gent et al., 2010). Although enhanced resolution in atmosphere and ocean models had a beneficial impact on a wide range of modes of internal variability, the relatively short high-resolution simulations make it difficult to sort that out in detail due to large decadal fluctuations in interannual variability in for instance ENSO (Sterl et al., 2007).
The requirement for a multitude of multi-centennial simulations, due to the
slow adjustment times in the Earth system, and the inclusion of Earth system
processes and feedbacks, such as those that involve biogeochemistry, have
meant that model resolution within CMIP has progressed relatively slowly. In
CMIP3, the horizontal typical resolution was 250 km in the atmosphere and
1.5
The primary goal of HighResMIP is to determine the robust benefits of
increased horizontal model resolution based on multi-model ensemble
simulations – to make this practical, vertical resolution will not be
considered. The argument for this is that the scaling between horizontal and
vertical resolution must obey
The use of process-based assessment is crucial to HighResMIP, since we aim to better understand the dynamical and physical reasons for differences in model results induced by resolution change, in order to increase our trust in the fidelity of models. Such process understanding will either contribute to bolstering our confidence in results from lower-resolution (but with greater complexity) CMIP simulations or to enabling a better understanding of the limitations of such models. There are an increasing number of studies suggesting that, in individual models, important processes are better represented at higher resolution, indicating ways to potentially increase our confidence in climate projections (e.g. Vellinga et al., 2016). A wide variety of processes will be assessed, from global and regional drivers of climate variability, down to mesoscale eddies in atmosphere and ocean – in the atmosphere these include tropical cyclones (Zhao et al., 2009; Bell et al., 2013; Rathmann et al., 2014; Roberts et al., 2015; Walsh et al., 2015) and eddy–mean flow interactions (Novak et al., 2015; Schiemann et al., 2016), while for the ocean they are an important mechanism for mesoscale air–sea interactions (Chelton and Xie, 2010; Bryan et al., 2010; Frenger et al., 2013; Ma et al., 2015, 2016), trans-basin heat transport (e.g. Agulhas leakage) (Sein et al., 2016), convection, and oceanic fronts.
HighResMIP will coordinate the efforts in the high-resolution modelling community. Joint analysis, based on process-based assessment and seeking to attribute model performance to emerging physical climate processes (without the complications of (bio)geochemical Earth system feedbacks) and sensitivity of model physics to model resolution, will further highlight the impact of enhanced horizontal resolution on the simulated climate. As the widespread impact of horizontal resolution, in the range of a few hundred to about 10 km, on climate simulation has been demonstrated in the past, it is expected that HighResMIP will contribute to many of the grand challenges of the WCRP, and hence such analysis may begin to reveal at what resolution in this range particular processes can be robustly represented.
The remainder of this paper is structured as follows. Section 2 gives an overview of the simulations, while Sect. 3 describes the tiers of simulation in detail. Section 4 makes links between these and the CMIP6 DECK and other CMIP6 MIPs, Sect. 5 describes the data storage and sharing plans, and Sects. 6 and 7 describe the analysis and potential application plans. Conclusions and discussion are contained in Sect. 8. Several appendices contain more detail of the experimental design and forcing.
Schematic outline of Tiers 1, 2, and 3. Tier 1 is a 64-year AMIP simulation from 1950 to 2014 with historical forcings. The first part of Tier 2 (coupled ocean–atmosphere simulations) consists of a 50-year integration starting from the 1950 initial state under 1950s conditions. Thereafter this simulation will be continued by two branches of 100 years: one continuing with the 1950s forcing (control run) and the other using until 2014 historical forcings and for 2015–2050 SSPx (scenario run). Tier 3 is the extension of Tier 1 from 2014 to 2050 (obliged, solid line) and 2051–2100 (optional, dashed line) for SSPx.
The main experiments will be divided into Tiers 1, 2, and 3. They are illustrated in Fig. 1. We provide an outline of these different tiers, with more detail in Sect. 3. Each set of simulations comprises model resolutions at both a standard and a high resolution, where the standard-resolution model is expected to be used in a set of CMIP6 DECK simulations and is considered the entry card for HighResMIP.
The Tier 1 experiments will be historical forced atmosphere (ForcedAtmos) runs for the period 1950–2014. A number of centres have already performed similar high-resolution simulations and published their results (CAM5 Bacmeister et al., 2014; HadGEM3 Mizielinski et al., 2014; NICAM Satoh et al., 2014; EC-Earth Haarsma et al., 2013); hence, these runs should not present prohibitively large technical difficulties. Restricting the ForcedAtmos runs to the historical period also makes it possible for numerical weather prediction (NWP) centres to contribute to the multi-model ensemble. Nineteen international groups have expressed interest in completing these simulations as shown in Appendix A. All centres participating in HighResMIP are obliged to participate at least in Tier 1.
The coupled experiments in Tier 2 are more challenging, but provide an opportunity to understand the role of natural variability, due to the centennial scale, and to investigate the impact of high resolution on future climate. Although a few centres have previously carried out high-resolution coupled simulations such as SINTEX-F2, GFDL, Hadley, MIROC, and CESM (Masson et al., 2012; Delworth et al., 2012; Mecking et al., 2016; Sakamoto et al., 2012; Small et al., 2014a), considerable issues including mean-state biases, climate drift, and ocean spin-up remain. Due to these issues and the large amount of computer resources needed to complete both a reference and a transient simulation, fewer centres (currently six) are confirmed participants for these experiments. The period of the coupled simulations is 1950–2050.
Future atmosphere-only simulations for the period 2015–2100 will be carried out in Tier 3. Although the future period covers the entire present century, the simulations can for computational reasons be restricted to the mid-century (2050).
For a clean evaluation of the impact of horizontal resolution, additional tuning of the high-resolution version of the model should be avoided. The experimental set-up and design of the standard resolution experiments will be exactly the same as for the high-resolution runs. This enables the use of HighResMIP simulations for sensitivity studies investigating the impact of resolution. If unacceptably large physical biases emerge in the high-resolution simulations, all necessary alterations should be thoroughly documented. The requirement of no additional tuning is more relevant for the coupled runs because atmosphere-only models are strongly constrained by the prescribed SSTs.
Forcings and initialization for the Historic simulations (pre-2015).
To focus on the impact of resolution on the design of the HighResMIP, simulations should minimize the difference in forcings and model set-up that would hamper the interpretation of the results.
Most of the forcing fields are the same as those used in the CMIP6 Historical Simulation that are described separately in this Special Issue (Eyring et al., 2016) and are provided via the CMIP6 data portal. For the future time period, GHG and aerosol concentrations from a high-end emission scenario of the Shared Socioeconomic Pathways (SSPs) will be prescribed, which in the following will be denoted by SSPx. A summary of the differences in forcing between the CMIP6 AMIPII protocol and the Tier 1 and 2 simulations is given in Table 1.
A potential large source of uncertainty is the aerosol forcing – for the same aerosol emissions, different models can simulate very different aerosol concentrations, hence producing different radiative forcing. In HighResMIP, each model will use its own aerosol concentration background climatology. To this will be added an anthropogenic time-varying, albeit uniform, forcing provided via the MACv2-SP method by Stevens et al. (2016). These will be computed using a new approach to prescribe aerosols in terms of optical properties and fractional change in cloud droplet effective radius to provide a more consistent representation of aerosol forcing. This will provide an aerosol forcing field that minimizes the differences between models as well as reduces the need for model tuning. This method is also the standard method in CMIP6 DECK for models without interactive aerosols.
The land surface properties will also be different from the CMIP6 AMIPII protocol. Given the requirement to make model forcing as simple as possible to aid comparability, the land surface properties will be climatological seasonally varying conditions of leaf area index (LAI), with no dynamic vegetation and a constant land use/land cover consistent with the present-day period, centered around 2000. Consideration was given to attempting to further constrain land surface properties to be more similar between groups, but this was rejected given the complex and different ways in which remotely sensed values are mapped to model land surface properties. However, an additional targeted experiment has been included to further investigate the sensitivity to land surface representation. This is outlined in Appendix C.
As discussed in Eyring et al. (2016), the initialization of land surface and atmosphere requires several years of spin-up to reach quasi-equilibrium before the simulation proper can begin. We recommend this is done using the first few years of the forcing datasets before restarting in 1950. We further recommend that the initial condition for the atmosphere and land for 1950 (for the highresSST-present and the highres-1950 experiment) come from the ERA-20C reanalysis from January 1950. If this is not possible, then the exact procedure used should be fully documented by each group.
The target for high resolution is 25–50 km, which is significantly higher than the typical CMIP5 resolution of 150 km. These ForcedAtmos runs will also be performed with the standard resolution that is used for the DECK and historical simulations.
The 1950–2014 simulation period is longer than the DECK AMIPII that spans 1979–2014. This is motivated primarily by work in many groups interested in climate variability over multi-decadal timescales, focusing on different phases of climate modes of variability such as Atlantic meridional oscillation (AMO) and Pacific decadal oscillation (PDO), as well as improved sampling of ENSO teleconnections (Sterl et al., 2007). The longer period will also improve the robustness of assessing the difference in variability between standard- and higher-resolution simulations, as well as being important for statistics of teleconnections (e.g. Rowell, 2013). Furthermore, the longer period of integration will enable a much more robust assessment of the ability of models to simulate known modes and their phases of variability, which is a big issue for climate risk assessment and decadal predictions where the combined effect of the global warming signal and natural variability will be considered.
The recommended ensemble size for the high-resolution simulations is three, but due to their computational cost many centres will probably be able to simulate only one member. Therefore although an ensemble size of three is recommended, it is not a requirement to participate in HighResMIP. The small ensemble size or absence of it will be insufficient for a rigorous estimate of the contribution of the internal variability to the total climate signal. However, by using a strictly common protocol in the various participating centres, the effective multi-model ensemble size will be much larger, enabling a much wider sampling than -pre-HighResMIP of the multi-model robustness of resolution impacts. In addition, if models can be proven to be portable, the ensemble size could be increased if auxiliary computer resources should become available at a later stage.
Although there is a significant forcing of the ocean by the atmosphere, in
particular at the mid-latitudes (Wu and Kirtman, 2007), many recent studies
have shown that gradients in SST associated with fronts and ocean eddies can
have a significant influence on the atmosphere via changes in air–sea fluxes
(Minobe et al., 2008; Parfitt et al., 2016; Ma et al., 2015; O'Reilly et al.,
2015). Similarly, there is evidence that daily variability rather than
monthly smoothed forcing can influence model simulations (de Boisséson et
al., 2012; Woollings et al., 2010). Since the high-resolution simulations
will approach 25 km, this means there is a requirement for a daily,
The coupled simulations are also aimed at addressing questions of model bias in both mean state and variability similar to the ForcedAtmos simulations. There are many examples from previous studies (e.g. Scaife et al., 2011; Bellucci et al., 2010) where these biases become much more evident in the coupled context compared to the forced atmosphere simulations. The systematic comparison between uncoupled (Tier 1) and coupled simulations for the 1950–2050 period, under different horizontal resolutions, will stimulate novel process-oriented studies tackling the origins of well-known biases affecting climate models, such as the double-ITCZ tropical bias.
These coupled runs will be the HighResMIP equivalent of the pre-industrial
control, here being a 1950s control using fixed 1950s forcing. The forcing
consists of GHG gases, including O
This will allow an evaluation of the model drift. The initial ocean
conditions are taken from version 4 of the Met Office Hadley Centre “EN”
series of datasets of global quality controlled ocean temperature and
salinity profiles and monthly objective analyses (EN4, Good et al., 2013)
over an average period of 1950–1954. As described below, a short spin-up
with these forcings is required (
These are coupled historic runs for the period 1950–2014 using an initial condition taken from Sect. 3.2.1.
For this period the external forcings are the same as in Tier 1 (see Table 1).
These are the coupled scenario simulations 2015–2050, effectively a continuation of the Sect. 3.3.2 historic simulation into the future. For the future period the forcing fields will be based on CMIP6 SSPx. Other forcings are detailed in Table 2.
The atmospheric component of the coupled models will be the same as in the
Tier 1 simulations. The minimum resolution for the high-resolution ocean
model is 0.25
The period of the historic coupled integrations is chosen to be the same as
in the Tier 1 simulations. The future end-date is based on a compromise
between what is computationally affordable by a sufficient number of centres
(
We again emphasize our interest in model error (bias, fidelity in representation of climate processes and variability) rather than climate sensitivity or transient climate response in configuring these coupled simulations, in particular whether any changes in process representation have an influence on patterns of climate variability and change. As discussed before, the number of ensemble members that will be possible, at least initially, in HighResMIP will not be sufficient to fully address internal variability, but it will form an important baseline set of simulations from which already preliminary robust conclusions can be extracted, and should be useful for many of the other CMIP6 MIPs (e.g. DCCP, GMMIP, CORDEX, CFMIP).
The HighResMIP simulations will enable the simulation of weather systems with short timescales that can provoke strong air–sea interactions such as tropical cyclones. Hence, high-frequency coupling between ocean and atmosphere is required: a 3 or 1 h frequency is highly recommended so that the diurnal timescale can be resolved, assuming sufficient vertical model resolution in the upper ocean.
Due to the large computer resources needed, a long spin-up to (near) complete
equilibrium is not possible at high resolution (and hence for consistency
will not be used at standard resolution). We recommend an alternative
approach which will use the EN4 (Good et al., 2013) analysed ocean state
representative of 1950 as the initial condition for temperature and salinity.
To reduce the large initial drift, a spin-up of
Forcings for the future climate simulations.
The Tier 3 simulations are an extension of the Tier 1 atmosphere-only simulations to 2050, with an option to continue to 2100. To allow comparison with the coupled integrations, the same scenario forcing as for Tier 2 (SSPx) will be used. However, since all the HighResMIP models will have the same SST and sea-ice forcing (described below), comparison of the Tier 2 and Tier 3 simulations can help to disentangle the impact of a model bias from forced response. This could be useful for applications such as time of emergence (e.g. Hawkins and Sutton, 2012). The forcing fields and scenario are shown in Table 2.
The future SST and sea-ice forcing are detailed in Appendix B. It broadly follows the methodology of Mizuta et al. (2008), enabling a smooth, continuous transition from the present day into the future. The rate of future warming is derived from an ensemble mean of CMIP5 RCP8.5 simulations, while the interannual variability is derived from the historic 1950–2014 period. Using SST derived from CMIP5 RCP8.5 in conjunction with a CMIP6 SSPx GHG forcing introduces an inconsistency. However, given the wide range of climate sensitivity among the climate models and the small differences in the model response up to 2050 for different scenarios, we argue that this inconsistency is minor.
In addition to the Tier 1–3 simulations above, discussions with other CMIP6 MIP participants have suggested several targeted experiments that would enable further investigation of specific processes and forcings, as well as potentially informing future CMIP protocols. These are optional experiments, and as such the details of the experimental design will be described in Appendix C. In brief they comprise the following.
Leaf area index (LAI) experiment – Impact of SST variability on large-scale atmospheric circulation –
Idealized forcing experiments with CFMIP – Abrupt 4 Tiers 2 and 3 using RCP8.5 instead of SSPx –
For the high-resolution models, completing the full set of CMIP6 DECK simulations is too expensive in terms of computer resources. Hence, there is an assumption that groups participating in HighResMIP will complete a set of DECK simulations with the standard-resolution model. The HighResMIP simulations will in that case be considered as sensitivity experiments with respect to the standard-resolution DECK runs, which are the entry cards for HighResMIP. If groups are not able to do this, because for instance the only available configuration is with prescribed SSTs, which is often the case for NWP centres, they can still participate in HighResMIP, but their simulations will only be visible as HighResMIP and not as CMIP6 runs.
Although there will be no DECK simulations at the high resolution, the comparisons between the standard-resolution simulations within HighResMIP and the DECK simulations will be informative in themselves. The relevance of HighResMIP is that the significant step in horizontal resolution enables us to clarify some of the outstanding climate science questions remaining from CMIP3 and CMIP5 exercises.
For the Tier 1 simulations, there is a strong link with the CMIP6 AMIPII simulations – the latter are likely to provide multiple ensemble members per modelling centre, but using slightly different boundary conditions and forcings (SST, sea ice, aerosols, LAI, and land use). Hence this comparison will be informative of the impacts of these changes at the standard resolution common to both AMIPII and HighResMIP. In addition, the multiple ensemble members will provide a measure of internal variability to assess whether the high-resolution simulation lies outside this envelope.
HighResMIP, as one of the CMIP6 endorsed MIPs, has links to a number of other MIPs. Collaboration with those will enhance the synergy.
There is well-known sensitivity of monsoon flow and rainfall to model resolution in the West African monsoon, Indian monsoon, and possibly East Asian monsoon. As stated in GMMIP, the monsoon rainbands are usually at a maximum width of 200 km. Climate models with low or moderate resolutions are generally unable to realistically reproduce the mean state and variability of monsoon precipitation for the right reasons. This is partly due to the model resolution. The Tier 1 ForcedAtmos runs of HighResMIP will be used in Task-4 of GMMIP to examine the performance of high-resolution models in reproducing both the mean state and year-to-year variability of global monsoons. As tropical monsoonal rainfall is sensitive to small-scale topography, high resolution has the potential to improve this. On the other hand, there is strong evidence of the importance of coupled ocean–atmosphere interactions for the summer monsoon dynamics (Robertson and Mechoso, 2000; Robertson et al., 2003; Wang et al., 2005; Nobre et al., 2012). Consideration was given to starting the HighResMIP from 1870 to better compare with GMMIP, but it would not be affordable for many groups. In addition, the quality of observational and reanalysis datasets during the earlier period, to assess the modelled variability and processes, is questionable.
HighResMIP intends to use the MACv2.0-SP simplified aerosol forcing being partly produced and analysed in RFMIP (Stevens et al., 2016). Additionally, assessment of its impact at different resolutions will contribute to understanding this simplified forcing. The impact of different aerosols on atmospheric circulation and teleconnections in the coupled climate system has been shown before and is likely dependent on model resolution (e.g. Chuwah et al., 2016).
CORDEX regional downscaling experiments provide focused downscaling over particular regions. Comparison between these and global HighResMIP simulations can give insight into the relative importance of global-scale teleconnections and interactions, against enhanced local resolution and local processes. HighResMIP can (potentially) provide boundary conditions for downscaling and provide a stimulus to cloud resolving simulations, but data volumes are likely to be prohibitive, so this will be left to individual groups to coordinate.
There is potential to jointly examine the spin-up issues in both forced ocean (OMIP) and coupled (HighResMIP) simulations, to improve the understanding of how we might better initialize coupled climate or forced ocean simulations and minimize initialization shock and the required integration time. The targeted experiment in Appendix C2 to understand the impact of mesoscale SST variability is another joint area of research. We will also exchange diagnostic/analysis techniques to understand ocean circulation changes at different resolutions.
Within the scope of LS3MIP on understanding the land–atmosphere interactions at different horizontal resolutions, HighResMIP can provide useful datasets to evaluate the role of soil moisture in extreme events, as well as the impact of LAI forcing datasets on model variability and mean state at different resolutions via targeted experiment in Appendix C1.
An increase in horizontal resolution may also improve the stratospheric basic state through vertical propagation of small-scale gravity waves, which in turn may affect tropospheric circulation. The sensitivity of such troposphere–stratosphere dynamical interactions to horizontal resolution will be analysed by the DynVAR community, and HighResMIP has actively coordinated with the DynVAR diagnostic request to make this possible.
Targeted experiments in Appendix C3, to look at the clouds and feedback response in different resolution models, can be assessed in conjunction with CFMIP experiments using the standard-resolution model.
Coordination of a sea-ice diagnostic request with SIMIP will enable coordinated assessment of the impact of model resolution on sea-ice conditions and processes. Indeed, sea-ice drift, deformation, and leads (Zhang et al., 1999; Gent et al., 2010) have been shown to be highly sensitive to model resolution in single-model studies. The robustness of these conclusions should be assessed in a coordinated multi-model exercise such as HighResMIP.
The storage and distribution of high-resolution model data are challenging issues. Since the resolution of HighResMIP approaches the scales necessary for realistic simulation of synoptic weather phenomena, daily and sub-daily data will be stored to allow the investigation of weather phenomena such as those related to mid-latitude storms, blocking, hurricanes, and monsoon systems. However, high-frequency output of all three-dimensional fields will not be affordable to store. Careful considerations are needed to limit the high-frequency output. The considerations should take into account that the information relevant for the end users is concentrated at or near the land surface where people live, so that it is desirable to store surface and near-surface variables at high temporal and spatial resolutions. Furthermore, in order to evaluate the HighResMIP ensemble, the high-frequency output should contain variables for which high-frequency observations are available as well.
HighResMIP output data will conform to all the CMIP requirements for
standardization. The CMIP6 data and diagnostic plan (Juckes et al., 2016)
describes the diagnostic request for all the CMIP6 MIPs. This data request,
including that of HighResMIP, is available from the CMIP6 website. The data
and diagnostic plan will be finalized during the boreal summer of 2016. An
estimate of the amount of data that need to be stored is given at
The data storage is divided into three priorities. This is based on a balance
between the HighResMIP data request
(
One useful approach may be to provide spatially and/or temporally coarsened model output on the ESGF, which would enable initial analysis compared to DECK simulations, and indicate which avenues of analysis may require full model resolution output, with manageable remaining volumes. It would then also be available for any automated assessment tools on the ESGF.
The analysis will focus on the impact of increasing resolution on the simulation of different climate phenomena that are strongly biased in coarse-resolution models and that could potentially benefit from higher resolution. The robustness of the impact of increasing resolution on the simulation of weather and climate phenomena such as extreme weather events, atmospheric eddy–jet stream interactions, atmospheric blocking events, typical ocean model biases, and ocean model drift among the different HighResMIP models will be investigated and their response to global warming assessed as well as their interannual variabilities.
The increased resolution will permit evaluation of whether horizontal resolution alone can generate a better simulation of regional climates. The analysis will therefore also have a focus on regional climate and relative teleconnections. Because HighResMIP will enable a more detailed simulation of small-scale weather systems, the scale interaction between these systems and the large-scale circulation will be another focus of the analysis plan. The benefit of atmosphere–ocean coupling at these high resolutions will be addressed as well since we can compare the AMIP-style simulations with fully coupled simulations. Not all modelling centres may be able to afford eddy-resolving ocean simulations; nevertheless, where possible, it will be interesting to investigate scale interactions in the ocean as well.
Five initial foci for analyses have been identified.
Current climate risk assessments rely on output from ensembles of relatively
coarse-resolution global climate models or on their downscaled products (e.g.
CORDEX) in addition to observations. For Europe, around 15 regional modelling
groups downscaled ERA-Interim simulations at 50 km and 12.5 km resolutions
(
A recent study by Jacob et al. (2014) showed that the high-resolution Euro-CORDEX simulations provide a more realistic representation of precipitation extremes over Europe and a larger increase in extreme precipitation in future simulations compared to the global models. Generally, the regional CORDEX simulations show a more sensitive response of precipitation to changes in greenhouse gas concentrations compared to their driving global models. However, the bias in the lateral boundary conditions from coarse resolution climate models can strongly affect the simulations in the regional models, such as shown for precipitation trends over Europe by van Haren et al. (2014, 2015).
The HighResMIP simulations will provide the first ensemble of global models with a comparable resolution to the current generation regional models. This will allow for a direct comparison of user-relevant parameters in HighResMIP to the CORDEX results. The comparison will focus on statistics and physics of meteorological events such as intense rainfall, droughts, storms, and heat waves. A comparison of the simulation of extreme events in the global models (which are self-contained and include global small-scale to large-scale interactions) and in regional models (forced at the boundary by another model, and typically a one-way downscaling) will be made. Results from various studies (e.g. Scaife et al., 2011; Kirtman et al., 2012), analysing the benefits of high resolution in the ocean in one single global model, indicate that increased resolution in global models leads to an improved simulation of large-scale phenomena such as the North Atlantic Current system and related surface temperature gradients. The impact of such improvements on blocking and storm tracks and the downstream effect on European climate variability and extremes will be analysed and compared to CORDEX results. Comparing HighResMIP results, with a globally high resolution, to results from both standard resolution global models and regional CORDEX simulations with a locally high-resolution domain (but boundaries based on coarse-resolution CMIP5 models) will give us insights into the importance of realistic large-scale climate conditions for local climate variations and extremes.
Studying internal variability of and long-term change in the Northern
Hemisphere sea-ice cover in the coupled HighResMIP simulations will enable us
to explore the impact of better resolved sea-ice dynamics on Arctic and
global climate. Preliminary tests conducted at
The improved simulation of synoptic-scale systems in HighResMIP enables us to analyse multi-scale phenomena such as large-scale circulation, tropical and extratropical cyclones, MJO, tropical waves, convection, and cloud in a seamless manner. For example, tropical cyclogenesis has known links to multi-scale phenomena including monsoon, synoptic-scale disturbances, and MJO (e.g. Yoshida and Ishikawa, 2013). Even for the dynamical storm track, which may be thought satisfactorily resolved by low-resolution climate models, its bias in latitudinal position is related to the cloud amount bias in CMIP5 models (Grise and Polvani, 2014). Existing high-resolution atmosphere simulations suggest that the characteristics of the jet stream (Hodges et al., 2011) and blocking (Jung et al., 2012) will be improved by higher resolution. The MJO and diurnal precipitation cycle are also of great interest. Such analysis, requiring high-frequency data, has implications for the output diagnostics – see Sect. 5 and Juckes et al. (2016).
In addition, the role of air–sea interactions at the mesoscale, such as analysed by Chelton and Xie (2010), Bryan et al. (2010), and Ma et al. (2015), can be assessed across models to understand the impact of resolution and the potential feedbacks in the system that may change the mean state.
Regarding the ocean, multi-scale phenomena can be discussed in a similar way.
By resolving eddies and having a lower dissipation due to refined resolution,
the cold bias in the north-western corner, the pathway of the Gulf
Stream/North Atlantic Current, the Southern Ocean warm bias, as well as the
Agulhas Current have been shown to be substantially improved (Sein et al.,
2016). Even at an intermediate
Process-level assessment of the simulated climate will give us some insights to improve the physics scheme in the climate models at a range of resolutions. Satellite simulators will be applied to the HighResMIP model output to evaluate cloud and precipitation processes in detail (e.g. Hashino et al., 2013). After the launch of the EarthCare satellite (planned in 2018; Illingworth et al., 2015), a new dataset including vertical distribution of cloud, precipitation, and vertical velocity is expected to be available. The fact that the horizontal resolution of the climate model is approaching that of the satellite observations also motivates us to accelerate synergetic studies between models and observations.
Process studies will aim to pin down the reasons for potentially better capturing small-scale and consequently large-scale phenomena with increasing resolution. Such process understanding will be the basis for developing schemes or error correction methods that could potentially compensate for not capturing a range of processes in standard-resolution models.
This topic has links with RFMIP (aerosols), LS3MIP (land surface processes), CFMIP (clouds), SMIP (sea ice), and DynVar (troposphere–stratosphere processes).
Many aspects of climate extremes are associated with the hydrological cycle, together with dynamical drivers such as mid-latitude storm tracks and jets. Analysis following Demory et al. (2014) will assess the multi-model sensitivity of the global hydrological cycle to model resolution, and convergence of moisture over land and ocean. In the tropics, the hydrological extremes due to monsoon systems and interactions between land and atmosphere (Vellinga et al., 2016; Martin and Thorncroft, 2015) will be investigated in conjunction with GMMIP. On a regional scale the extremes and hydrological cycle will be analysed in collaboration with CORDEX. For extremes associated with surface processes, there are links with LS3MIP.
At mid-latitudes, the representation of storm tracks and jet streams will be assessed. Novak et al. (2015) investigated the role of meridional eddy heat flux in the tilt of the North Atlantic eddy-driven jet. This behaviour may partly explain the dominant equatorward bias of the jet stream in generations of global climate simulations with model resolutions much coarser than 50 km (Kidston and Gerber, 2010; Barnes and Polvani, 2013; Lu et al., 2015). Biases in the jet stream position have been found to correlate with the meridional shift of the jet position in a warmer climate (Kidston and Gerber, 2010).
Atmospheric rivers play a key role in the global and regional water cycle (Zhu and Newell, 1998; Ralph et al., 2006; Leung and Qian, 2009; Neiman et al., 2011; Lavers and Villarini, 2013), and hydrological extremes, and have been shown to be sensitive to model resolution (Hagos et al., 2015). In both the North Pacific and North Atlantic, uncertainty in projecting atmospheric river frequency has been linked to uncertainty in projecting the meridional shift of the jet position in the future (Gao et al., 2015, 2016; Hagos et al., 2016), with consequential impacts on robust predictions of regional hydrologic extremes in areas frequented by land falling atmospheric rivers.
With the high-resolution simulations resolving more realistic orographic features in western North and South America and western Europe (Wehner et al., 2010), this motivates more detailed analysis of regional precipitation and hydrologic extremes, including changes in the amount and phase of extreme precipitation, snowpack, soil moisture, and runoff and rain-on-snow flooding events in a warmer climate than have been attempted previously with the coarser-resolution CMIP3 and CMIP5 model outputs.
Recent studies (Walsh et al., 2012, 2015; Shaevitz et al., 2014; Scoccimarro et al., 2014; Villarini et al., 2014) have highlighted the benefits of enhanced model resolution for the representation of several aspects of tropical cyclones (TCs), including the formation patterns, genesis potential index, and the relative impact on precipitation. HighResMIP will provide an ideal framework to systematically investigate the influence of model resolution on the representation of tropical cyclones in the next generation of climate models.
It is expected that by improving the representation of the background, large-scale (oceanic and atmospheric) pre-conditioning factors affecting TC dynamics (such as wind shear and ocean stratification) via a refinement of model resolution, the overall representation of TC properties (including structure and statistics) will be affected. The potential remote influence of TCs on high-latitude processes suggested by a few authors – e.g. TC impacts on sea-ice export in the Arctic region (Scoccimarro et al., 2012), extra-tropical transition (Haarsma et al., 2013), and extreme precipitation events over Europe (Krichak et al., 2015) – is another (so far, poorly explored) topic that may benefit from the HighResMIP multi-model effort.
Finally, the 1950–2050 time window targeted in HighResMIP experiments will allow an evaluation of the stationarity of the relationship between TC frequency and intensity, and the underlying, large-scale environmental conditions (Emanuel, 2015).
Given the relatively short time period for integration and small ensemble size, and the fact that Tier 3 simulations are also limited by using atmosphere-only models, we must give careful consideration to the applications for which the HighResMIP simulations can be used.
Below is a non-exhaustive list of additional issues, not discussed in the analysis plan, that can be addressed by HighResMIP.
A better assessment and attribution of the changes in extreme events that are
already occurring and of near-future changes will provide useful information
for regional climate adaptation strategies and other users of climate model
output such as infrastructure investments that have a time horizon of up to
30 years. The benefit relates to the increased physical plausibility and
reliability of simulating the circulation-driven aspects of the weather
extremes, which are more biased in coarser-resolution climate models. The
ensemble could aid in developing scenarios of potential future weather events
to which society is vulnerable (Hazeleger et al., 2015) and be used for
impact studies such as ecosystem studies, meteo-hydrological risks, and
landslides.
Another potential use of these simulations is to give a baseline of the forced response only (using the best estimate of the SST forced response and the SSPx radiative forcing) for near-term decadal predictions. This can then be combined with coupled decadal predictions (or statistical modelling) that also include the ocean variability and its influence. See for instance Hoerling et al. (2011) as a first attempt to do this with low-resolution models.
HighResMIP will for the first time coordinate high-resolution simulations and process-based analysis at an international level and perform a robust assessment of the benefits of increased horizontal resolution for climate simulation. As such it is an important step in closing the gap between climate modelling and NWP, by approaching weather resolving scales. A better representation of multiple-scale interactions is essential for a trustworthy simulation of the climate, its variability, and its response to time-varying forcings and boundary conditions. HighResMIP thereby focuses on one of the three CMIP6 questions: “what are the origins and consequences of systematic model biases?”. Specifically it will investigate the relation of these model biases to small-scale systems in the atmosphere and ocean and how well they are represented in climate models.
Despite the importance of enhancing horizontal resolution, many processes still have to be parameterized. For processes and regions where these parameterizations are crucial, increasing horizontal resolution did not improve the model bias. The role of various parameterizations in model biases will be investigated in other MIPs, for instance in AerChemMIP, CFMIP, and RFMIP. Jointly they will address the grand challenges of the WCRP from different angles.
HighResMIP will address the grand challenges of the WCRP in the following way.
HighResMIP will address this grand challenge by investigating the sensitivity to increasing resolution of water vapour loading, cloud formation, and circulation characteristics, with analysis concentrating on the relevant processes (see Sect. 6.3).
To improve the robustness of our understanding, the multi-model ensemble at different resolutions, together with the longer AMIP integrations, will allow us to
link tropospheric circulation to changing patterns of SSTs and land surface properties, and understand the role of cloud processes in natural
variability; examine the extent and limits of our understanding of patterns of
precipitation; and examine changes in model biases (such as humidity) with resolution, since there are some indications that these may be linked to climate sensitivity.
Increasing resolution affects in particular small-scale process such as the formation of clouds. Although the formation of clouds has still to be parameterized under the typical resolution used within HighResMIP, the dynamical constraints for the formation of clouds, such as the location and magnitude of upward and downward motion associated with frontal systems and orography, as well as moisture availability, are sensitive to resolution. This also applies to the response of the circulation to cloud formation.
HighResMIP is very relevant to this grand challenge. Resolution affects the
hydrological cycle by modifying the land–sea partitioning of precipitation.
Increasing resolution in general increases the moisture convergence over land
(Demory et al., 2014), although regionally this can be reversed, such as for
instance in Europe during the winter due to changes in the position of the
storm track (Van Haren et al., 2014). In addition, simulations of extreme
precipitation events are highly sensitive to increasing resolution. How
robust are these results across the multi-model ensemble? Can
higher-resolution models help to give insight into inconsistencies between
global precipitation and energy balance datasets? How surface water
availability (
HighResMIP is strongly related to this grand challenge. Increasing resolution of climate models will bring us closer to the ultimate goal of seamless prediction of weather and climate. Extremes mostly occur and are driven by processes on small temporal and spatial scales that are not well resolved by standard CMIP6 climate models. Dynamical downscaling only partially resolves this limitation due to the non-linear interaction between large and small spatial scales and the importance of representing global teleconnection patterns. We aim to improve our understanding of the interaction between global modes of variability (e.g. ENSO, NAO, PDO) and regional climate inter-decadal variability and extremes, as well as between local topographic features and the triggering of extreme events.
Regional climate information focuses on smaller scales and extreme events, which are relevant for stakeholders and adaptation strategies. This requires high-resolution modelling to provide reliable information. Increasing resolution globally allows one to better capture not only local processes that could be captured by regional climate models, but also teleconnections with distant regions which could have a strong impact on the region of interest. Recent high-resolution modelling studies (Di Luca et al., 2012; Bacmeister et al., 2014) and comparisons of CMIP3 and CMIP5 results (Watterson et al., 2014) have demonstrated the added value of increased resolution for regional climate information. Model outputs from HighResMIP could also be used by the regional climate modelling community for comparison of dynamical downscaling and global high-resolution approaches and for further dynamical downscaling by cloud resolving regional models and statistical downscaling for impact assessments.
In the Tier 2 coupled simulations, the better representation of sea-ice deformation, drift, and leads as well as heat storage and release with increased resolution can contribute to better capturing the growth and motion of sea ice, the air–sea heat flux, and deepwater production in polar regions, processes that are strongly affected by small-scale processes. Based on HighResMIP coordinated simulations we can make a robust assessment of the effect of model resolution on Arctic sea-ice variability, including sea-ice circulation and export through the Fram and Davis straits, and possible influences on mid-latitude circulation. Analysis of the cryosphere in the Tier 1 experiments will, however, be somewhat limited due to the prescribed sea-ice distribution. Its main impact will be on the distribution of snow fall and subsequent accumulation and melting of the snowpack that affect land surface hydrology.
The simulations in HighResMIP will obviously be demanding with respect to high-performance computing capability, particularly in order to complete them in a reasonable time frame. There are ongoing efforts to acquire supra-national resources in Europe and elsewhere, and the Tianhe-2 supercomputer, one of the most powerful systems in the world, also offers huge computing resources to support HighResMIP in China.
HighResMIP has evolved from the need to harmonize existing projects of high-resolution climate modelling. European Horizon2020 project PRIMAVERA, in which major European climate centres are participating, has coordinated the initiatives for a common protocol within the CMIP6 framework. As such, the simulations conducted in PRIMAVERA will be first under the HighResMIP protocol.
It is expected that HighResMIP will be a major step forward in entering the area of weather resolving climate models and thereby opening new avenues of climate research. Fundamental new scientific knowledge is expected on weather extremes, the hydrological cycle, ocean–atmosphere interactions, and multiple-scale dynamics. As such, it will contribute more trustworthy climate projections and risk assessments.
The model output from the DECK and CMIP6 historical simulations will be
distributed through the Earth System Grid Federation (ESGF) with digital
object identifiers (DOIs) assigned. As in CMIP5, the
model output will be freely accessible through data portals after
registration. In order to document CMIP6's scientific impact and enable
ongoing support of CMIP, users are obligated to acknowledge CMIP6, the
participating modelling groups, and the ESGF centres (see details on the CMIP
Panel website at
Model details from groups expressing intention to participate in at least Tier 1 simulations, together with the potential model resolutions (if known/available, blank if not).
Discussion with the HighResMIP participants suggests that the agreed
approach is to use the RCP8.5 scenario, and use the CMIP5 models to generate
the projected future trend. Numerical code for the following calculations
will be made available in Python, as will the final dataset on the
So, following Mizuta et al. (2008) for the most part, the algorithm is described below.
For HadISST2.2.0 (Rayner et al., 2016) in the period 1950–2014:
For each year
Calculate, from the monthly mean, the time mean of the period
Calculate the linear monthly trend
Finally calculate the interannual variability
Then from at least 12 CMIP5 coupled models during the period 1950–2100
(using the Historic and RCP8.5 simulations), calculate a monthly mean trend, for each model over this period, as a
difference from several years centered at 2014, so that the change in
temperature can be smoothly applied to the HadISST2 dataset.
Regrid this trend to the HadISST2
Calculate the multi-model ensemble mean of this monthly trend.
This ensemble mean still contains a large component of both spatial and
temporal variability – since the object here is to produce a large-scale,
smoothly varying background signal to the HadISST2 variability, this
multi-model trend is spatially filtered (using a 20
Then for the future period, the temperature is
This will repeat the variability from the past period into the future, but adding the model future trend. The choice of 1950 as a start date for this section is because it has the most similar phase of some of the major modes of variability (AMO, PDO, etc.) to use for the repeat.
Concatenate this section (twice) to the end of HadISST2 at 2014:
Projecting the sea-ice into the future will be based on the following
procedure.
Using observed SST and sea-ice concentration, an empirical relationship is
constructed. HadISST2 (Rayner et al., 2016) uses the inverse method to derive
SST based on sea-ice concentration). This is done by dividing the SST into bins of 0.1 K. The SST of each data
point determines in which bin the sea-ice concentration of each data point
falls. After all data points are handled in this way the mean sea-ice
concentration for each bin is computed. The relationship is different for the
Arctic and Antarctic and seasonally dependent. Using this empirical relationship between SST and sea-ice
concentration,
the sea-ice concentrations for the constructed SST are computed.
However, a couple of alternative methods are also being investigated, such as that used in HadISST2 (Titchner and Rayner, 2014), in which the sea-ice edge is located, and then the concentration is filled in from here towards the pole.
The LAI is one of the most common vegetation indices that describe vegetation activity (Chen and Black, 1992). It closely modulates the energy balance, as well as the hydrological and carbon cycles of the coupled land–atmosphere system at different spatiotemporal scales (Mahowald et al., 2016). For atmosphere–ocean GCMs, including those of HighResMIP, the mean seasonal cycle of LAI is commonly prescribed to improve the physical and biophysical simulations of the land–atmosphere system (Taylor et al., 2011). To reduce the potential uncertainties due to inconsistent LAI inputs for different models participating in HighResMIP, we propose conducting targeted LAI experiments, with a common LAI dataset.
Various remote sensing based LAI datasets have been recently developed (Fang
et al., 2013; Zhu et al., 2013). Among them, the LAI3g data have been found
to be the best in terms of continuity, quality, and extensive applications
(Zhu et al., 2013; Mao et al., 2013). For the targeted experiments we will
provide a
The impact of mesoscale air–sea coupling on the large-scale circulation (in atmosphere and ocean) is a growing area of research interest. Ma et al. (2015) have shown that mesoscale SST variability in the Kuroshio region can exert an influence on rainfall variability along the US North Pacific coast. In order to assess this, we propose parallel simulations of the high-resolution ForcedAtmos model using spatially filtered SST forcing.
The modelling approach is to conduct twin experiments – one with high-resolution SST (the reference HighResMIP simulation) and another with spatially low-pass filtered SST. This approach appears to be quite effective in dissecting the effect of mesoscale air–sea coupling. The filter should be the LOESS filter used by Ma et al. (2015) and Chelton and Xie (2010). The parallel simulation should start in 1990 from the HighResMIP simulation and be identical apart from the SST forcing.
Period of integration: 10 years. This should be done in an ensemble multi-model approach to ensure statistically significant results.
CFMIP experiments using
Period of integration: 10 years for each
CFMIP experiments use abrupt 4
Period of integration: 20 years (in parallel with the first 20 years of control-1950, after the initial spin-up period).
This option is included for centres, such as those involved in European H2020 project PRIMAVERA, that have to start their simulations before the availability of SSPx. It is motivated by the notion that the differences between SSPx and RCP8.5 will be limited up to 2050. If in a joint analysis the SSPx and RCP8.5 ensembles appear to be significantly different, then the RCP8.5 centres are recommended to repeat their simulations with SSPx, which, due to the short integration period of 36 years, should not be prohibitive.
PRIMAVERA project members (Malcolm J. Roberts, Reindert J. Haarsma, Pier Luigi Vidale, Torben Koenigk, Virginie Guemas, Susanna Corti, Jost von Hardenberg, Jin-Song von Storch, Wilco Hazeleger, Catherine A. Senior, Matthew S. Mizielinsky, Tido Semmler, Alessio Bellucci, Enrico Scoccimarro, Neven S. Fučkar) acknowledge funding received from the European Commission under grant agreement 641727 of the Horizon 2020 research programme.
Chihiro Kodama acknowledges Y. Yamada, M. Nakano, T. Nasuno, T. Miyakawa, and H. Miura for analysis ideas.
Neven S. Fučkar acknowledges support of the Juan de la Cierva-incorporación postdoctoral fellowship from the Ministry of Economy and Competitiveness of Spain.
L. Ruby Leung and Jian Lu acknowledge support from the U.S. Department of Energy Office of Science Biological and Environmental Research as part of the Regional and Global Climate Modeling Program. The Pacific Northwest National Laboratory is operated for the DOE by Battelle Memorial Institute under contract DE-AC05-76RLO1830.
Jiafu Mao is supported by the Biogeochemistry-Climate Feedbacks Scientific Focus Area project funded through the Regional and Global Climate Modeling Program in Climate and Environmental Sciences Division (CESD) of the Biological and Environmental Research (BER) Program in the U.S. Department of Energy Office of Science. Oak Ridge National Laboratory is managed by UT-BATTELLE for the DOE under contract DE-AC05-00OR22725.
Paulo Nobre acknowledges support from CNPq grant nos. 573797/2008-0 and 490237/2011-8, and FAPESP grant no. 2008/57719-9.
Chihiro Kodama and Masaki Satoh are supported by the Program for Risk Information on Climate Change (SOSEI) and the FLAGSHIP2020 within the priority study4 (Advancement of meteorological and global environmental predictions utilizing observational “Big Data”), which are promoted by the Ministry of Education, Culture, Sports, Science and Technology (MEXT), Japan.
Ping Chang is supported by US National Science Foundation grants AGS-1462127 and AGS-1067937, and National Oceanic and Atmospheric Administration grant NA11OAR4310154, as well as by China's National Basic Research Priorities Programme (2013CB956204 and 2014CB745000).
We thank Martin Juckes and his team for all their work on the HighResMIP and CMIP6 data request.
Nick Rayner and John Kennedy for allowing early access to the HadISST2 daily,