In atmospheric models with kilometer-scale grids the resolution approaches the scale of convection. As a consequence the most energetic eddies in the atmosphere are partially resolved and partially unresolved. The modeling challenge to represent convection partially explicitly and partially as a subgrid process is called the convective gray zone problem. The gray zone issue has previously been discussed in the context of regional models, but the evolution in regional models is constrained by the lateral boundary conditions. Here we explore the convective gray zone starting from a defined global configuration of the Met Office Unified Model using initialized forecasts and comparing different model formulations to observations. The focus is on convection and turbulence, but some aspects of the model dynamics are also considered. The global model is run at nominal 5km resolution and thus contributions from both resolved and subgrid turbulent and convective fluxes are non-negligible. The main conclusion is that in the present assessment, the configurations which include scale-aware turbulence and a carefully reduced and simplified mass-flux convection scheme outperform both the configuration with fully parameterized convection as well as a configuration in which the subgrid convection parameterization is switched off completely. The results are more conclusive with regard to convective organization and tropical variability than extratropical predictability. The present study thus endorses the strategy to further develop scale-aware physics schemes and to pursue an operational implementation of the global 5km-resolution model to be used alongside other ensemble forecasts to allow researchers and forecasters to further assess these simulations.

Global studies of climate change impacts that use future climate model projections also require projections of land surface changes. Simulated land surface performance in Earth System models is often affected by climate biases of the atmospheric models, leading to errors in land surface projections. Here we run the JULES-ES land surface model with ISIMIP2b bias-corrected climate model data from 4 global climate models (GCMs). The bias correction reduces the impact of the climate biases present in individual models. We evaluate JULES-ES performance against present-day observations to demonstrate its usefulness for providing required information for impacts such as fire and river flow. We simulate a historical and two future scenarios; a mitigation scenario RCP2.6 and RCP6.0, which has very little mitigation. We include a standard JULES-ES configuration without fire as a contribution to ISIMIP2b and JULES-ES with fire as a potential future development. Simulations for gross primary productivity, evapotranspiration (ET) and albedo compare well against observations. Including fire improves the simulations, especially for ET and albedo and vegetation distribution, with some degradation in shrub cover and river flow. This configuration represents some of the most current earth system science for land surface modelling. The suite associated with this configuration provides a basis for past and future phases of ISIMIP, providing a simulation setup, postprocessing and initial evaluation using ILAMB. This suite ensures that it is as straightforward, reproducible and transparent as possible to follow the protocols and participate fully in ISIMIP using JULES.