The sensitivity of the tropical circulation and Maritime Continent precipitation to climate model resolution

The dependence of the annual mean tropical precipitation on horizontal resolution is investigated in the atmospheric version of the Hadley Centre General Environment Model. Reducing the grid spacing from about 350 km to about 110 km improves the precipitation distribution in most of the tropics. In particular, characteristic dry biases over South and Southeast Asia including the Maritime Continent as well as wet biases over the western tropical oceans are reduced. The annual-mean precipitation bias is reduced by about one third over the Maritime Continent and the neighbouring ocean basins associated with it via the Walker circulation. Sensitivity experiments show that much of the improvement with resolution in the Maritime Continent region is due to the specification of better resolved surface boundary conditions (land fraction, soil and vegetation parameters) at the higher resolution. It is shown that in particular the formulation of the coastal tiling scheme may cause resolution sensitivity of the mean simulated climate. The improvement in the tropical mean precipitation in this region is not primarily associated with the better representation of orography at the higher resolution, nor with changes in the eddy transport of moisture. Sizeable sensitivity to changes in the surface fields may be one of the reasons for the large variation of the mean tropical precipitation distribution seen across climate models.


Introduction
It is widely held that scale interactions (between short/ small and long/large scales in time and space) are an important feature of the climate system, in particular in the tropics (e.g., Slingo et al. 2003Slingo et al. , 2009)).Consequently, a change in the resolution of a climate model may be expected to impact on the simulated climate, not only at the truncation scale but also at much larger scales.Indeed, systematic resolution dependencies have been described for a number of models and phenomena (e.g., Manabe et al. 1970;Boville 1991;Senior 1995;Williamson et al 1995;Stephenson et al. 1998;Pope and Stratton 2002).While it is well known that the increase of resolution is not a panacea for removing climate-model biases (e.g., Slingo et al. 2003), beneficial impacts have been documented in a range of more recent studies both for the atmosphere and the ocean.Roberts et al. (2004) show that increasing the ocean resolution in a coupled global circulation model to 1/3°allows the ocean model to represent eddies and is associated with an improved mean ocean state.In the southeastern Pacific, the HiGEM model at 1/3°oceanic resolution represents small-scale (\450 km) transient features and associated heat transports that agree better with available observations than in the coarser-resolution (*1°9 1°) Hadley Centre General Environment Model (HadGEM1) (Toniazzo et al. 2009).At sufficient resolution in the atmosphere and the ocean, Roberts et al (2009) and Shaffrey et al (2009) find for the HiGEM model that small-scale interactions between the atmophere and instability waves in the tropical Pacific can be captured realistically improving the simulation of the mean state of the tropical Pacfic, the Walker circulation, and of all aspects of ENSO.Strachan et al (2012) report improvements in the number and geographical distribution of tropical cyclones in a series of atmosphere-only simulations with the HadGEM1 model; and Demory et al (2013) find, using the same model simulations, that the global water and energy cycles are represented more realistically at higher resolution.
Despite regular improvements, the representation of the (tropical) mean climate continues to be a challenge for state-of-the-art general circulation models (GCMs).Notably the tropical precipitation distribution has been difficult to simulate [e.g., Randall et al (2007), section 8.3.1.2,Martin et al (2010), Gent et al (2011)].Characteristic precipitation biases occur in atmosphere-only simulations and may be exacerbated in coupled models (Hack et al. 2006).Precipitation is a key predictive variable from an environmental and socioeconomic point of view.Moreover, precipitation biases are an indicator of tropical heating errors and limit the confidence in GCMs in applications such as seasonal forecasting (e.g., Wang et al. 2007) or the identification of remote (tropical) forcing of events in the extratropics (e.g., Rodwell and Hoskins 1996).These biases are very pronounced over the Maritime Continent and the neighbouring ocean basins (Neale and Slingo 2003;Qian 2008).This is the region of largest precipitation and tropical heating, and the correct representation of its climate is crucial for both the tropical and extratropical mean state.For example, the ability of a model to capture atmospheric moisture transport and precipitation over the Maritime Continent is tightly linked to its ability to represent the zonal tropical overturning circulation, and may be a prerequisite for correctly simulating the present and future intensity and location of the Walker circulation (Held and Soden 2006;Vecchi et al. 2006).
In this study, we investigate the sensitivity of the tropical mean climate to horizontal resolution in a recent version of the Met Office Hadley Centre climate model.More specifically, we aim to characterize the changes in the mean tropical precipitation due to a decrease of the atmospheric grid spacing from about 350 to 110 km, and to advance our understanding of why and how the model exhibits resolution dependence.We use a series of model experiments designed to this end; and focus on the Maritime Continent region where model resolution might be of particular relevance due to the complex arrangement of land and sea (Neale and Slingo 2003).
The outline of the paper is as follows.Section 2 describes the GCM used here, the setup of the model experiments, and observations and reanalysis data.In Sect.3, we give an overview of how the mean tropical climate depends on resolution and raise hypotheses of what may cause the resolution dependence.In Sect.4, sensitivity experiments are described that assesses the dependence of the simulated climate on the resolution of surface boundary fields over the Maritime Continent.The paper is concluded in Sect. 5.

Models, methods, and data
The atmospheric version of the Hadley Centre General Environment Model (adGEM1; Johns et al. 2006) is used at different horizontal resolutions (N48-N216; see Table 1) and otherwise very similar model formulation.The model configurations at N48 and N96 resolution have a time step of 30 min and are the same as in Strachan et al (2012).
The N144 model was built for this study.It uses the HadGAM1 model formulation and the same surface boundary condition data as the atmospheric part of the High-resolution General Environmental Model [HiGEM; Shaffrey et al (2009)].The time step for this model is 20 min.Small adjustments (of the settings in the N48 and N96 model) are made to the dynamics of this model, including the vertical velocity threshold for targeted diffusion and polar filtering, to permit it to run multi-decadal climate simulations.Our sensitivity studies indicate that these adjustments have minimal impact on the simulated climate compared to the impact of resolution changes.It should also be noted that, for historical reasons, the land fraction in the N144 model is based on a slightly different dataset than in the N48 and N96 models.1This is reflected in a somewhat higher land fraction in the N144 model (Table 2).
We have also conducted sensitivity experiments with the N96 and N144 models where we regionally replace the original (N96 or N144) surface boundary fields with coarse (N48) surface data.These experiments are described in detail in Sect. 4.
In addition to the simulations at N48-N144 resolutions, simulations at N216 resolution and with a time step of 15 min are available.These are based on a model configuration that is still similar to the N48-N144 models, yet some more modifications were necessary to ensure model stability at the even higher resolution.These modifications include the reduction of the Convective Available Potential Energy (CAPE) timescale from 60 min at N48-N144 to 30 min at N216, enabling horizontal diffusion, and a reduction of the autoconversion parameter.Our sensitivity experiments show that the impact of these modifications on the simulated climate might not be negligible.We therefore briefly include a result of the N216 model in Sect.3, but our conclusions regarding the resolution-dependence of the HadGAM1 model climate are based on the N48-N144 models.
We analyse AMIPII-style [Atmospheric Model Intercomparison Project; Taylor et al (2000)] simulations covering the period 1979-2002.Three such simulations, i.e. ensemble members, generated by perturbing the least significant bit of the numerical values of the initial potential temperature field, are available for the N48, N96, and N216 models.For the N144 model that was specifically built and run for this study, one simulation is available.One year of spin-up time is discarded in all simulations so that the reference period for all analyses in this study is 1980-2002.We use version 2.1 of the Global Precipitation Climatology Project (GPCP) satellite-gauge product as the observational reference data for precipitation, HadSST2 data for sea surface temperature, CRUTEM3 data for land temperature, and evaluate atmospheric fields against ERA-Interim data (see also Table 1).

Resolution dependence of the tropical mean climate
In this section, we give an overview of the mean-state biases of HadGAM1 and of how these biases depend on resolution.At all resolutions, the model qualitatively captures many of the large-scale features of the global precipitation distribution (Fig. 1a-c).This includes the Intertropical Convergence Zone (ITCZ), the South Pacific Convergence Zone (SPCZ) and the midlatitude stormtracks, as well as the subtropical subsidence regions.At the same time, biases can be identified when comparing to the observed precipitation (Fig. 1d-f).There is too much precipitation over the oceans, in particular in the western part of the basins in the tropics.HadGAM1 has pronounced dry biases over the Indian Subcontinent and Southeast Asia including the Maritime Continent.These are long-standing errors in the HadGEM1 family of climate models (e.g., Martin et al. 2006;Shaffrey et al. 2009) and they are of the same order of magnitude as the total model precipitation in some regions (for example, in the Maritime Continent box shown in Fig. 1, the annual mean precipitation is 2.64 mm day -1 in the N48 model and 6.23 mm day -1 in GPCP).
The increase in resolution reduces the magnitude of the biases in the western tropical oceans, the Maritime Continent, and, to a smaller extent, over India.This can be seen by comparing the differences between the N48 and the higher-resolution models (Fig. 1g, h) with the biases shown in Fig. 1d-f.The effect of increased resolution and the biases are of opposite sign in many regions.The difference between the N96 and N144 models is much smaller than the difference between each of these models and the N48 model.There are two reasons that may explain this.Firstly, the increase in resolution from N48 to N96 is a larger step (about 350 to 170 km grid spacing) than the increase from N96 to N144 (about 170 to 110 km).Secondly, the regional mean land fraction is higher in the N144 model (see Sect. 2 and Table 2).Since the warm ocean is the main source of evaporation in the region, the increase in latent heat flux due to the increase in resolution appears to be partly offset by the fact that there is slightly less evaporation from the smaller ocean surface in the N144 model.In the N216 model (Fig. 2), the sign and geographical distribution of The columns show the name, resolution, grid spacings (in degrees longitude, latitude, and in km at the equator), the time step in minutes, and references.Numbers in parenthesis refer to the spatial grid on which interpolated observational datasets are made available the precipitation bias are similar to those in the N48-N144 models, while the magnitude of the bias tends to be somewhat smaller.The box over the central part of the Maritime Continent (95.625-140.625°E,-8.75 to 3.75°N; Fig. 1) corresponds to a region where there is a dry precipitation bias at all resolutions, and where the impact of higher resolution is such as to increase precipitation throughout the region.This domain will be used for the analyses discussed later in the paper.
Figure 3 shows the land fraction, precipitation P, evaporation E, P-E, and midtropospheric vertical velocity meridionally averaged in a tropical channel (covering the same latitudes as the domain highlighted by the box in Fig. 1).The land fraction shows the approximate location of the African, Maritime, and American continents.At all resolutions, the contrast between a comparatively dry Maritime Continent and wet Indian and West Pacific oceans is larger in HadGAM1 than in GPCP or ERA-Interim.The evaporation in the model is larger than in ERA-Interim over the oceans, and about the same over land.Accordingly, the bias in P-E is smaller than the precipitation bias in regions where the precipitation bias is positive (such as the West Pacific), and larger where the precipitation bias is negative (such as the East Pacific).In general, P-E biases follow the precipitation biases, in particular over land.P-E, i.e. regions of climatological moisture convergence (divergence), are in turn closely related to regions of ascent (descent) as shown by the midtroposheric vertical velocity.It is notable that the Maritime Continent is a region of time-mean descent in the models, in contrast to ERA-Interim.
The effect of resolution is to reduce the dry bias over the Maritime continent by about one third.There are reductions of similar magnitude of the wet biases over the We formulate two hypotheses that may explain the observed resolution dependence.First, we test if higher resolution allows for more moisture to be transported in atmospheric ''eddies''.The second hypothesis concerns the surface boundary conditions prescribed to the atmospheric model.These surface fields include orographic fields, soil and vegetation parameters, and the fractional land-sea mask.As the resolution is increased, surface fields are prescribed with greater detail.It is clear that these changes have the potential to modify many aspects of the simulated climate such as the atmospheric circulation, the surface fluxes, and also the precipitation (e.g., Neale and Slingo 2003).We test the first hypothesis in the remainder of this section; the second hypothesis is addressed in terms of sensitivity experiments described in Sect. 4.

Mean and eddy transport of moisture
We have seen that the resolution dependence of P-E is very similar to the resolution dependence of precipitation.Since r Á ðqvÞ ¼ PÀE, moisture convergence is a suitable quantity to investigate the role of eddies for the observed resolution dependence.To this end, we decompose the moisture flux convergence into mean and eddy component where q ¼ " q þ q 0 and v ¼ " v þ v 0 is the decomposition of the specific humidity and the horizontal velocity into the monthly-mean and eddy (deviation from the monthly mean at each time step) part.The " qv 0 and q 0 " v terms in Eq. ( 1) are zero by construction.The integration is carried out over the whole mass of air in the atmospheric column.The result of this decomposition is shown in Fig. 4 for the N48 and N96 models.In the tropics, the dominant contribution to moisture convergence is in the mean term.In agreement with the previous discussion, the effect of increasing the resolution is to increase mean moisture convergence over the Maritime Continent, and to decrease it over the neighbouring ocean areas (Fig. 4c).The resolution dependence of mean moisture convergence is very similar to the resolution dependence of the precipitation shown in Fig. 1 and it is larger than the differences in the eddy part shown in Fig. 4f.The same results are obtained for the N144 model (not shown).This suggests that moisture transported in eddies does not constitute the main part of the climatological resolution dependence, and so we reject the first hypothesis formulated above.

Role of surface boundary conditions
We now turn to the second hypothesis formulated in the preceeding section and investigate the role of prescribed (or seasonally varying) surface boundary conditions at different resolutions.As discussed above, some sensitivity to the surface fields can be expected, yet it is not straightforward to estimate a priori what the net consequences for the simulated climate should be.We therefore conduct the sensitivity experiments described below.

Methodology
As the resolution of a model is increased, boundary data are prescribed at a higher resolution, and a part of the subgrid-scale variability in the boundary fields at the lower resolution will be explicitly represented at the higher resolution.We aim to quantify the impact this has on the simulated climate, and to separate this effect from other changes due to increased resolution.To this end, we use the high-resolution (N96 or N144) atmospheric models, but regionally replace surface boundary conditions by their coarse-resolution (N48) counterpart.This is illustrated in Fig. 5 for the fractional land-sea mask.The top panels (Fig. 5a, c) show the original land-sea masks of the N96 and N144 models, and the bottom panels (Fig. 5b, d) show the land-sea masks obtained by replacing the high-resolution mask by the coarse-resolution mask over the Maritime Continent (in the region outlined by the black box).
To retain the coarse-resolution structure of the N48 boundary data, the value of a single N48 grid cell is indentically assigned to all corresponding grid cells at the higher resolution.For the N144 model, this is straightforward since nine grid cells of the N144 model exactly overlap one grid cell of the N48 model as shown in Fig. 6a.
For the N96 model, there are a number of ways to match an N48 cell to four N96 cells, and we choose to shift the N48 values to the southwest by one half of the grid spacing of the N96 model (Fig. 6b).
Coarse surface boundary fields are introduced for the land fraction, for orographic fields, for soil moisture, soil temperature, and a number of soil parameters; as well as for leaf-area index, canopy height, and fractional coverage by vegetation types.We refer to the experiments where all these surface boundary fields have been replaced by the N48 fields as 'N96DA' and 'N144DA'.Additionally, we test the effect of introducing coarse-resolution orographic fields only (and leaving all other surface boundary fields unchanged).We conduct this experiment with the N96 model only and refer to it as 'N96DO'.

Results of the sensitivity experiments
The impact of introducing coarse surface boundary conditions on the mean precipitation is shown in Fig. 7.Over the Maritime Continent, where the modified surface fields have been applied, coarse boundary conditions lead to a reduction of precipitation, both in the N96DA and N144DA experiments.
In the Indian Ocean, over the Maritime Continent, and the West Pacific, the response to coarsening the surface fields is similar to the impact of resolution (compare Figs. 1g,h,7a,b).This suggests that the change in the model-simulated climate with resolution is largely due to the associated change in the surface boundary conditions.The sensitivity is not dominated by the orographic fields as shown by the smaller and more localized response in the N96DO experiment shown in Fig. 7c.This implies that the response is largely due to the different representation of the coastline and of land-surface parameters at different resolutions.
We proceed with an analysis of the surface climate in the central Maritime Continent region introduced in Fig. 1. Figure 8 shows the precipitation, temperature, net radiation, and latent heat flux at the surface.When compared with the observed precipitation and the re-analysis data, all model experiments show less precipitation, higher landsurface temperature, and larger net surface radiation.The latent heat flux over land is smaller than in ERA-Interim, whereas the bias is relatively small over the ocean where the surface temperature is prescribed.The biases in all quantities are consistent with a drier and warmer climate in HadGAM1 than seen in ERA-Interim and observations.The effect of higher resolution is to decrease the model biases (compare the N96 and N144 model with the N48 model for the following discussion).The precipitation bias with respect to GPCP in the N96 model is about 40 % smaller than in the N48 model [(4.06 -2.64)/(6.23 -2.64)].The reduction of the bias in the N144 model is somewhat smaller, about 30 %.In both the N96 and the N144 model, the dry precipitation bias is reduced by a larger amount over land than over the ocean.
The biases in the other variables undergo a consistent reduction with resolution.The land-surface temperature decreases below the value of the sourrounding oceans bringing the N96 and N144 models into better agreement with ERA-Interim than N48.Net radiation descreases over land and changes only little over the ocean.Latent heating increases with resolution, albeit only by a very small amount in the N144 model over the ocean.The overall effect of higher resolution is a cooler and moister model climate in better agreement with the reference data.
Much of the resolution dependence described above can also be seen in the sensitivity experiments using coarse surface boundary fields.For example, the precipitation both over the land and the ocean is only slightly larger in the N96DA and N144DA experiments than in the N48 model.Net radiation over land increases due to the introduction of coarse surface fields, whereas the latent heat flux decreases both over land and ocean.This further suggests that the improvement with higher resolution of the simulated climate of the Maritime Continent is predominantly Conversion between mm day -1 and W m -2 is based on a latent heat of condensation of 2.5 9 10 6 J kg -1 .Ordinates do not start with zero due to the fact that surface boundary conditions are prescribed in greater detail as the resolution is increased.

Discussion of mechanisms
Further insight into the increase in the Maritime Continent precipitation and surface latent heat flux with resolution can be gained from a more detailed consideration of the land surface scheme.In common with many other climate models, HadGEM1 uses a coastal-tiling surface scheme (Martin et al. 2006), where coastal grid-boxes are partitioned into land and ocean fractional tiles at the surface level.Diagnostic fields used in the calculation of surface fluxes have different values for the land and sea parts at the surface of the coastal grid-box, whereas no distinction between land and sea is made at the lowest atmospheric level (Fig. 9).The vertical exchange coefficients governing the fluxes depend on quantities both at the surface and at the lowest atmospheric level, and therefore also have different values over the land and sea parts of the coastal grid-box.
The formulation of the coastal tiling scheme can result in unrealistic surface fluxes over coastal grid-boxes (Strachan 2007).We illustrate this here for the surface latent heat flux in the N144 model (Fig. 10).Figure 10a shows the land fraction; the colour bar is such that the lowest (highest) level corresponds to pure sea (land) grid-boxes, and the intermediate three levels correspond to coastal grid-boxes.The surface latent heat flux is shown in Fig. 10b.Over the coastal grid-boxes surrounding the major Maritime Continent islands the latent heat flux is smaller than over nearby grid-boxes over the open sea and in the interior of the islands.
Over coastal grid-boxes, the latent heat flux is the mean of the flux over the sea and land parts of the grid-box weighted by the land fraction.Figure 10c shows the latent heat flux from pure sea grid-boxes and from the sea part of coastal grid-boxes, and Fig. 10d shows the complementary latent heat flux from pure land grid-boxes and the land part of coastal grid-boxes.Over most of the coastal grid-boxes, surface lowest atmospheric level T S , q S , CH S T L , q L , CH L T 1 , q 1 , u 1 , v 1 , ... ), e wind speed at 10 m over pure sea and coastal grid-boxes (m s -1 ).Pure land or sea gridboxes where the respective diagnostic is not defined are shown in grey the latent heat flux is substantially smaller over the sea part of the grid-box than over the land part.This shows that the reduction of the latent heat flux over coastal grid-boxes (Fig. 10b) is largely due to the small latent heat flux over the sea part of coastal grid-boxes (Fig. 10c).
To investigate further, the 10-m wind speed over pure sea and coastal grid-boxes is shown in Fig. 10e.The 10-m wind speed is smaller over coastal grid-boxes than over pure sea grid-boxes.Comparison of Fig. 10b, c, e shows that the spatial variation of the wind is fairly closely associated with that of the latent heat flux over both sea and coastal grid-boxes.The reduction of the latent heat flux over the coastal grid-boxes is consistent with the formulation of the coastal tiling scheme: It is plausible that over coastal grid-boxes the increased surface roughness over the land part of the grid-box will reduce the wind speed at the lowest model level.The reduced lowest model-level wind is then used to calculate a reduced surface latent heat flux over the sea part of the coastal grid-box.
It is clear that the suppression of the latent heat flux over coastal points can lead to resolution dependencies.To show this explicitly, the differences in surface latent heat flux averaged over N48 grid-boxes were compared between the N144 control and the degraded land surface (N144DA) experiments over the Maritime Continent.The differences in the surface latent heat flux as a function of land fraction can be seen in Fig. 11.When the land fraction within a N48 gridbox is close to zero or one (i.e.all sea or land) the surface latent heat fluxes in the N144 control and N144DA experiment are similar on average.For an N48 grid-box with a land fraction of 0.5 (i.e.50 % land and 50 % ocean) the latent heat flux in the N144 control experiment is roughly 15-20 W m -2 larger than that of the N144DA experiment.We have verified that this difference is statistically significant.
Since the Maritime Continent is a region of complex landsea contrasts, a large number of grid-boxes are prescribed as being coastal.As the resolution increases then naturally the fraction of the prescribed coastal grid-boxes decreases (Table 3).Subsequently this leads to an increase in the surface latent heat flux averaged over the Maritime Continent.We argue that the response seen in the precipitation is due to the direct forcing of the surface wind and turbulent fluxes by the modified surface boundary fields, and a feedback associated with an enhancement of moisture convergence and the Walker circulation.We note in passing that the fraction of coastal grid-boxes is comparatively similar in the N96 and N144 models (Table 3), which provides an additional possible explanation of the rather small difference in the mean state between these two models discussed in Sect.3.
These results suggest two possible conclusions.Firstly, by reducing the need to parametrise coastal regions, increasing the resolution of climate models enables land surface schemes to better represent regions such as the Maritime Continent and to reduce associated long-standing biases in tropical precipitation.Secondly, these results also highlight a need to improve land surface schemes so that coastal tiling parametrisations do not exhibit such large resolution dependencies when used in coarse resolution climate models.Indeed, a so-called ''buddy'' scheme for coastal grid-boxes has been introduced in a more recent model version (HadGEM3-GA3.0,Walters et al. 2011), which uses an average wind speed over neighbouring gridboxes to split the wind speed at the lowest atmospheric level into separate land and sea contributions.This enhances the wind speed over the sea part of the grid-box giving improved scalar fluxes there and a better mean state over the Maritime Continent region.

Summary and discussion
We have investigated how tropical precipitation biases in an atmosphere-only GCM (AGCM) depend on the model's   ).Each point in the scatterplot corresponds to the mean of 3 9 3 N144 grid-boxes, i.e. to a box of the N48 grid (as in Fig. 6a).The black line is obtained by LOWESS smoothing with span f = 0.5 (Cleveland 1979(Cleveland , 1981) ) horizontal resolution.To this end, we have used ensembles of simulations of present-day climate with the HadGAM1 AGCM at grid spacings between about 350 and 110 km, but otherwise very similar model formulation.We find that • the increase in resolution reduces long-standing tropical precipitation biases in HadGAM1.For example, the dry bias over the Maritime Continent is reduced by about one third, and associated improvements are seen over the Western Pacific and western Indian oceans.• the increase in precipitation over the Maritime Continent at higher resolution is associated with an enhancement of the Walker circulation and increased moisture convergence.Decomposition of the moisture convergence into its mean and eddy parts shows that the resolution dependence is predominantly seen in the mean part.• the resolution dependence is largely due to the specification of better resolved surface boundary conditions (such as the land-sea mask, soil and vegetation parameters) as the resolution is increased.This finding is based on sensitivity experiments where we have replaced surface boundary conditions in higher resolution models (N96 and N144 resolutions) with those from a low (N48) resolution model over the Maritime Continent.The sensitivity is not mainly associated with orographic surface fields as seen in an experiment where only these fields are changed.
In this study we find that the surface latent heat flux depends nonlinearly on the land fraction in the Had-GAM1 models.Over the Maritime Continent, the latent heat flux over coastal points is smaller than that over pure land or ocean points.As the model resolution increases, there are relatively fewer coastal points and therefore the regionally averaged latent heat flux increases.This suggests that increasing the resolution of climate models improves the representation of complex and heterogeneous regions, such as the Maritime Continent, by reducing the need to parameterise coastal areas.This can have sizable impacts on long-standing biases in the ability of climate models to simulate tropical precipitation.
It was also found that a rather similar reduction in precipitation biases was found in the N96 and N144 models.This suggests that the sensitivity of Maritime Continent precipitation to resolution in the HadGEM1 family of models may be reduced for resolutions higher than N96.This result is also consistent with preliminary results from an N216 resolution model.One key question is whether the conclusions of this study will be found in other atmospheric climate models that have very different formulations and parameterisations.It will be particularly useful to address this question with climate model experiments that systematically test the impact of model resolution on the representation of the climate system.Future research will test the conclusions of this study in a new set of climate model experiments with the HadGEM3 climate model (Hewitt et al. 2011) currently being performed (Mizielinski et al., submitted to Geosci Model Dev) with horizontal resolutions ranging from N48 (*350 km at low latitudes) all the way up to N512 (*30 km).

Fig. 2 Fig. 3
Fig. 1 a-c Annual mean precipitation in the N48, N96, and N144 models, d-f biases with respect to GPCP, and g, h resolution depence.Here and in subsequent figures: stippling shows grid points where a

Fig. 4
Fig. 4 Moisture-flux convergence in the N48 and N96 models.a, b Mean part and d, e eddy part for the two resolutions, and c, f difference in the mean and eddy part between the resolutions

Fig. 5
Fig. 5 Modification of the fractional land-sea mask over the Maritime Continent.a Original mask of the N96 model, b the mask of the N48 model is regionally introduced into the N96 model, c original mask of the N144 model, d the mask of the N48 model is

Fig. 6 Fig. 7 Fig. 8 a
Fig. 6 Schematic of the relative alignment of grid cells in a the N48 and N144 models, and b the N48 and N96 models

Fig. 9 Fig. 10
Fig.9Schematic of the coastal tiling scheme.At surface level, the model grid-box is devided into a land and a sea part, with pertinent variables such as temperature T and specific humidity q, and also exchange coefficients CH.No distinction between land and sea is made at the lowest atmospheric level

Fig. 11
Fig.11Difference in the annual-mean latent heat flux between the N144 control simulation and the N144 simulation with coarse surface boundary fields over the Maritime Continent (N144DA).Each point in the scatterplot corresponds to the mean of 3 9 3 N144 grid-boxes, i.e. to a box of the N48 grid (as in Fig.6a).The black line is obtained by LOWESS smoothing with span f = 0.5(Cleveland 1979(Cleveland , 1981) )

Table 1
Model experiments, observations, and reanalysis data

Table 2
Land fraction (%) at different resolutions, averaged over the whole globe and over the Maritime

Table 3
Area fraction (%) of coastal grid-boxes at different resolutions, averaged over the whole globe and over the Maritime Continent (95.625 -140.625°E,-8.75 to 3.75°N)