Contrasting contributions of surface hydrological pathways in convection permitting and parameterised climate simulations over Africa and their feedbacks on the atmosphere

The partitioning of rainfall at the land surface into interception, infiltration and surface runoff plays an important role in the water cycle as it controls the time scale at which water returns to the atmosphere. Rainfall intensity is of crucial importance to this partition. High resolution convection permitting models significantly improve simulated sub-daily rainfall intensity distributions, in particular those associated with convective rainfall in the tropics. Here we compare the land surface hydrological response in a pair of 10-year simulations over an African domain performed using the Met Office Unified Model: a typical configuration using parameterised convection operating at 25 km and the second a high resolution convection permitting simulation at 4.5 km with the parametrized convection switched off. Overall pan-African interception in the convection permitting scheme is 70% lower, whilst surface runoff is 43% higher than the parameterized convection model. These changes are driven by less frequent, but more intense rainfall with a 25% increase in rainfall above 20 mm h−1 in the 4.5 km model. The parameterised scheme has a ~ 50% canopy water contribution to evaporative fraction which is negligible in the convection permitting scheme. Conversely, the convection permitting scheme has higher throughfall and infiltration leading to higher soil moisture in the weeks following rain resulting in a 30–50% decrease in the daytime sensible heat flux. We examine how important the sub-grid rainfall parameterisation in the model is for the differences between the two configurations. We show how, switching a convective parameterisation off can substantially impact land surface behaviour.


Introduction
When rainfall reaches the land surface, it follows different hydrological pathways interception by vegetation, infiltration and percolation through the soil column and runoff. Each pathway has distinct characteristics within the overall water cycle. The vegetated surface accounts for a relatively small reservoir of surface water storage, filling and emptying within hours, yet evaporation of this intercepted water can reach up to ~ 15% of annual rainfall in tropical rainforests (Miralles et al. 2011). The water that reaches the ground, either directly, or falling from the vegetation canopy store, can then infiltrate into the larger soil store. This soil water store is depleted relatively slowly, via transpiration, direct evaporation from the soil surface, and drainage. The third pathway is surface runoff, which occurs when rainfall is intense (infiltration excess runoff), or the soil is already saturated (saturation excess runoff). The physical properties of the surface e.g. vegetation cover and surface infiltrability are important for determining the partitioning between these pathways, as is the rainfall intensity. Light rainfall may be sufficient to wet a vegetated surface but insufficient to contribute to infiltration and surface runoff, whilst higher rainfall rates will quickly exceed canopy capacities and activate the soil surface, with water being partitioned into infiltration (increased soil moisture) and runoff. How the rainfall is partitioned is critical for plant productivity, agriculture and water resources. Equally, the different evaporative time scales associated with these pathways provides an important feedback on the atmosphere which climate models need to capture (Lawrence et al. 2007;Scott et al. 1995). The representation of this partitioning at the land surface is challenging for global circulation models (GCMs) (Lian et al. 2018) as the key processes are sub-grid, and can occur in a matter of minutes. This is further compounded by the fact that the rainfall produced by GCMs have significant biases in the frequency and intensity of simulated rainfall (Fiedler et al. 2020;Stephens et al. 2010). These models operate on spatial scales of the order of 100 km whilst the very intense rainfall-producing storms are an order of magnitude smaller, requiring parameterisations that relate temperature and humidity profiles to convective rainfall. As a result GCMs tend to have incorrect statistics of daily and sub-daily precipitation, producing too many low intensity rain events. Within the framework of parameterised convection schemes, several land surface models include statistically-based parameterisations that attempt to account for the sub-grid variability of the incoming rainfall within a single GCM grid box (Best et al. 2011;Wang et al. 2009;Decharme and Douville 2007;Schmidt et al. 2006). Furthermore steps to improve the runoff generating characterises include additional parameterisation of surface heterogeneity which also impact soil moisture (Beven and Kirkby 1979;Moore 1985;Wood et al. 1992).
Experiments with GCMs have shown sensitivity of canopy evaporation to near surface air temperature and humidity and in turn the development of clouds (Lawrence et al. 2007;Wang et al. 2009;Davies-Barnard et al. 2014). A recent study of evaporation partitioning into plant transpiration and bare soil and canopy evaporation in CMIP5 models (Lian et al. 2018) found interception losses contributed to incorrect partitioning of total evaporation, with rainfall frequency a likely contributor to this overestimate.
Since problems with the partitioning of water in climate models is attributed to biases in rainfall intensity as well as the treatment of rainfall by the land surface scheme, recent advances to improve rainfall characteristics have impacted the hydrological partitioning. Super-parameterizations in GCMs (Demott et al. 2007) improve the depiction of rainfall intensity and frequency, and subsequent canopy interception, which also led to improved coupling strength between the land and atmosphere over tropical forests (Sun and Pritchard 2016). Renewed interest in the hydrological partitioning of rainfall from the climate modelling community has also been triggered by the increasingly widespread use of models with much lower horizontal grid spacings (< 5 km), allowing the deep convective parameterisation to be switched off (Prein et al. 2015) and more recently in ensemble simulations (Ban et al. 2021;Pichelli et al. 2021). These convection-permitting models (CPMs) simulate the vertical updrafts and down draughts within single storms, and can resolve the broader scale organization of convection into mesoscale convective systems (Clark et al. 2016). They tend to improve the diurnal timing of rainfall (Clark et al. 2016), the longevity of storms , the distribution of sub-daily rainfall intensities (Senior et al. 2021), and the depiction of land-atmosphere feedbacks (Taylor et al. 2013;Hohenegger et al. 2009). Such models have been routinely used for short-term operational forecasts for over a decade e.g. the UK Met Office (MetUM) configuration (Lean et al. 2008), but advances in computing power and more efficient codes allow multi-year regional simulations, and even month-long global simulations (Stevens et al. 2019).
The ability to run CPM simulations longer than a few days means that any changes in partitioning of rainfall into its different pathways due to improved spatial and temporal model resolution can have potentially important consequences for the atmosphere via feedbacks in evaporative response. For example, Hu et al. (2020) have shown how the intense rain produced by organised convective systems in the United States can percolate deeper, and affect the time scale of feedback with subsequent storms. Several authors have noted a tendency towards drier soil moisture stores when comparing CPM runs with parameterised models. For example (Ban et al. 2014) speculated on the role of more frequent light rain simulated by their convective parameterisation for enhancing soil moisture deficit in a 10-year European simulation. In a seasonal CPM simulation over China, (Li et al. 2020) also found increased sensible heat flux at the expense of latent heat which they argued was contributing to biases in the diurnal cycle of rainfall. Over India, Willetts et al. (2017) attributed increased sensible heat in a CPM simulation to increased insolation coupled with soil moisture deficit from reduced interception loss and increased runoff. They found that this enhanced the land-sea pressure difference (compared to a parameterised simulation) and increased onshore flow.
In this study we provide a detailed examination of how switching from parameterised convection to convectionpermitting in a regional climate model affects the partition of rainfall into different hydrological pathways, and how this partitioning feeds back into the simulation of the atmosphere. We focus on a pair of 10-year simulations over a pan-African domain using the UK Met Office model ). This model is increasingly being used in long convection-permitting simulations, and previous studies with it have noted important changes in surface hydrological behaviour when switching to convection-permitting mode e.g. Li et al. (2020), Kendon et al. (2012) and Willetts et al. (2017). By looking in detail at the surface hydrology parameterisation in this model, our analysis sheds new light on those studies. However, our focus on changes in surface hydrological fluxes and stores is also relevant for long integrations using other CPMs which have reported enhanced dry soil biases e.g. (Ban et al. 2014) compared to coarser resolution simulations.
In Sect. 2 we describe the model and experimental design, with a particular focus on the surface hydrology. This is followed by analysis of the continental scale differences in rainfall and its partitioning in Sect. 3.1, the impact of these differences on the diurnal cycle via interception feedbacks (3.2), and the feedback on the atmosphere via soil moisture on seasonal time scales (3.3). The results are discussed in Sect. 4.

Data and models
This study takes advantage of a pair of 10-year pan-African simulations using the UK Met Office Unified Model designed to understand how the representation of convection affects African climate and climate change (Senior et al. 2021). One simulation (R25) is run with a conventional convective parameterization on a horizontal grid of ~ 26 × 39 km over a domain from 39 N to 45S, 25 W to 56 E (Fig. 1). A second convection-permitting simulation (CP4), is run on the same domain, with the same boundary conditions, but with a regular horizontal grid of 4.5 × 4.5 km at the equator, and parameterized convection switched off. The cloud schemes used in the two models are different. CP4 has 80 model levels and uses the diagnostic cloud scheme (Smith 1990) whereas R25 has 63 model levels and uses the prognostic PC2 cloud scheme (Wilson et al. 2008). Rainfall in CP4 is determined only by the large-scale precipitation scheme. Both models are driven by 3-hourly lateral boundary conditions from a global 25 km simulation, N512L85 protype version of the Global Atmosphere 7 (GA7) and Global Land 7 (GL7) (Walters et al. 2019), which is forced with present day sea surface and lake temperatures for the years 1997 to 2006, see Stratton et al. (2018) for more details. The CP4 simulations are computationally demanding, but have provided a unique continental-scale depiction of climate change over tropical land with a convection-permitting model. Many analyses have already been published on these simulations (Senior et al. 2021), which have demonstrated advances in the representation and understanding of tropical land-based convection and how this might change over the 21st Century. Comparisons with observations have shown that CP4 exhibits more realistic rainfall properties than R25, including improvements in the diurnal cycle , intensity distribution (Senior et al. 2021;Berthou et al. 2019) and storm size and duration . None of the studies to date have focused on how the representation of land-atmosphere interactions is affected by the representation of convection, which is the topic of this work. We analyse only the two simulations under current climate conditions, and leave analysis of how land-atmosphere coupling within a convection-permitting framework may influence climate change projections to future work. Whilst more general details of the model set up are described elsewhere ), below we outline aspects of the land surface model which are key to understanding differences in land-atmosphere coupling between the two simulations. The Joint UK Land Environment Simulator (JULES) land surface model (Best et al. 2011) provides the lower boundary of energy and mass fluxes to the atmosphere. Fluxes are computed on 9 land tiles according to 5 plant functional types (broadleaved tree, needle leaved tree, shrub C4 and C3 grasses) and 4 non-plant functional types (urban, water, bare soil and ice), which share a single four level soil moisture store. The configuration of JULES is based on the "GL7" configuration (Wiltshire et al. 2020) and is used in both CP4 and R25.
The vegetation fractions used here were derived from European Space Agency-Climate Change Initiative land cover datasets (Poulter et al. 2015), providing mapping from biome to model plant functional type, and give a more realistic distribution of bare soil. A climatology of leaf area index derived from MODIS collection 5 dataset provide seasonal vegetation changes. Previous convectionpermitting studies over West Africa found the mapping of soil properties introduce artificial boundaries that can unduly influence the sensitive soil moisture -precipitation feedbacks (Taylor et al. 2013). Therefore in these simulations spatially uniform properties of sand were adopted. Soil moisture initialisations consistent with these uniform soil properties were generated from offline simulations driven by WATCH Era Interim Forcing data (Weedon et al. 2011). The 0.5 degree soil moisture climatology was computed from 10 years of simulation and downscaled to the CP4 and R25 model resolutions.
The JULES model simulates a canopy water balance above four soil layers on which soil moisture and temperature are computed. The canopy water balance is solved for each tile, with each tile having its own maximum canopy capacity (up to 1.05 mm for broadleaf trees), and which depends on leaf area index. The canopy store is replenished by rainfall and depleted by evaporation. Canopy evaporation is rapid as there is no evaporative resistance. Rainfall in excess of the canopy capacity reaches the soil surface as throughfall. Through fall is partitioned between the top level soil moisture and surface runoff. The latter includes both infiltration-excess runoff (when throughfall rates exceed the infiltration capacity) and saturation excess runoff. In this configuration saturation excess runoff (that is surface runoff produced when rain falls onto saturated ground) is modelled by the probability distributed moisture scheme (Moore, 1985), and related to the soil moisture in the top 1 m (3 levels). The infiltration capacity is based on the product of the saturated soil hydraulic conductivity and a vegetation enhancement factor, which takes values of 2 for grass and 4 for trees. Each grid box has soil moisture reservoir which is shared cross all the tiles, and divided into four soil layers of depths 0.1, 0.25, 0.65 and 2.0 m depth. Flow between the soil layers is computed following Richards equation for unsaturated flow, with water draining from the lowest level being sub-surface runoff. The relationship between soil matric potential and soil moisture and unsaturated hydraulic conductivity is defined by van Genuchten (1980). Transpiration is extracted from the soil profile based on an exponential root distribution. Bare soil evaporation occurs only from the top level, and sub-surface runoff is made up of water that drains out from the lowest soil level.
To address the problem of modelling sub-grid variability in surface hydrological fluxes over a GCM grid box described in the Introduction, the Unified Model assumes an exponential sub-grid rainfall distribution when calculating the canopy water balance and surface runoff (Dolman and Gregory 1992).
where P l is the distribution of rainfall intensities for a given grid box mean rainfall, P and the rainfall area fraction, µ controls the shape of the distribution of rainfall intensities. For rainfall from the convective parameterisation, µ varies with rainfall rate, typically taking a value ~ 0.3, whilst for rainfall from the large-scale scheme, µ = 1. Whilst the land surface parameters specified in the two model configurations are identical, they differ significantly in their assumed subgrid rainfall distribution, because the overwhelming majority of rain in R25 is produced by the convective scheme (so µ ~ 0.3), whilst CP4 rain is all large-scale (µ = 1).
The exponential distribution is used to derive throughfall, where C is the canopy water, C m the maximum canopy capacity, and Δt is the model time step. Infiltration excess runoff is given by (Best et al. 2011;Dolman and Gregory 1992), where K is the surface infiltration rate, and is computed from the saturated hydraulic conductivity multiplied by an infiltration enhance factor, which in turn is related to the land cover type. Eqs. 2 and 3 are applied to each tile, and grid box mean fluxes and canopy water computed by area-weighting the tiles. The difference between grid box mean values of throughfall and infiltration excess runoff enter the top model soil level. During the CP4 simulations, the soil moisture was inadvertently reset to saturation when the model was restarted on 1/7/1997 i.e. 6 months into the simulation. A repeat 12-month simulation initialised with the correct soil moisture on 1/7/1997 was performed, and the output of the model patched into the 10-year time series. This approach corrected the large but short-lived impact on surface fluxes (see Appendix A of Stratton et al. 2018). However the patching in of the data did result in some discontinuities in the time series of deep (i.e. below the root zone) soil moisture, and sub-surface runoff in arid and sparsely-vegetated regions (notably the Sahara). This patching process didn't have any noticeable impact on the atmosphere at the monthly time scale, and its impact on the terrestrial water budget is discussed later.
The simulations output hourly-mean land surface fluxes, with accumulated rainfall available hourly in R25 and every 15 min in CP4. In our analysis of the diurnal cycle, these variables have been adjusted to local solar time. The results presented in this paper are computed using 9 years of simulation from 1/1/1998 to 30/12/2006, and hourly and daily diagnostics.

Large scale water balance and partitioning
To begin, we compare continent-wide terrestrial water budgets between the two simulations. Mean annual totals for mainland Africa (excluding Madagascar) are presented in Table 1. Averaged across Africa, mean rainfall is 1.86 mm day −1 for both models, (Table 1), though there are notable rainfall differences within the continent (Fig. 1). Despite comparable rainfall inputs, the two simulations partition the surface water very differently. At the pan-Africa scale (Table 1) total evaporation (that is transpiration, bare soil and canopy evaporation) in R25 is higher and represents a higher proportion of the surface water budget (75% of rainfall) compared to CP4 (62% of rainfall). Crucially, it is the difference in modelled interception (0.23 mm day −1 ) that accounts for the vast majority of the difference in total evaporation (0.24 mm day −1 ) between the two models. Evaporation of water intercepted by the canopy in R25 is a remarkable 3.6 times larger than in CP4. The two simulations also differ substantially in total runoff. Annual runoff Table 1 Long term water balance components in mm day −1 (percentage of rainfall) over mainland Africa  (2002); (e) Rodell et al. (2015)  is higher in CP4 (40% of rainfall) compared to R25 (25%). This difference is made up of contributions from both surface runoff and subsurface runoff. We note that because of the issues discussed above concerning the restart of CP4 and its erroneous wet soil re-initialisation, the deep soil moisture was not fully spun-up in drier parts of the domain, even after the remaining 8.5 years of simulation. This results in a small increase (0.04 mm/day) in sub-surface runoff in unvegetated arid areas (e.g. the Sahara) where losses are by drainage only, but overall the differences in sub-surface runoff between the simulations are small. On the other hand, we found that the impact of the soil moisture re-initialisation on evaporation was negligible after 12 months. This is because spin-up issues in level 4 soil moisture were only evident in dry regions (notably the Sahara). In these areas there is no deep-rooted vegetation to extract soil moisture, and hence no ability to influence transpiration.
In Table 1 we also include observational estimates of the different terms in the pan-African water balance, where available. Comparing with mean rainfall from the CHIRPS v2 dataset (Funk et al. 2015), the simulations are close to the observed over the same years (1998 to 2006), though at sub-continental scale, rainfall biases in one or both simulations can be larger, as previously presented by Stratton et al. (2018) and Jackson et al. (2020). For the remaining water fluxes, the observational estimates come with large uncertainties and are included here only to provide a "ballpark figure" to compare with the two simulations. Considering total evaporation, R25 sits within the estimated range, whereas CP4 evaporation lies 12% below the minimum estimate. The only available global estimate of canopy evaporation (Martens et al. 2017) suggests a 320% over-prediction by R25, but very similar interception loss to CP4. Finally, we note that pan-African runoff, R25 again is within the estimated range whereas CP4 over-predicts by up to 85%, though the observational datasets include losses due to lakes, reservoirs, abstractions and within river losses which are not represented by these simulations. Comparisons with observational estimates of evaporation and runoff over a region as sparely-observed as Africa are fraught with difficulty. However, we do note that the estimate of continent-wide interception losses as a percentage of rainfall from R25 does also appear to be excessive in comparison with in situ observations from tropical broadleaved evergreen forests; 8.9% (Lloyd et al. 1988); 11.6%; (Urabana 1997) 13.3% (Cuartas et al. 2007); 13.5% (Tobon Marin et al. 2000); 11.6% (Czikowsky and Fitzjarrald 2009); 9.1% (Shuttleworth 1988) in Amazonia and 9.2% (Hutjes et al. 1990) in West Africa. As interception loss should decrease strongly away from the densely-vegetated forests of Central and Western Africa, we conclude that the continent-wide interception loss of 18% in R25 is not realistic. In summary, at the continental scale the two models have similar total rainfall but quite distinct partitioning of that rainfall. The CP4 simulation produces much less evaporation than R25, and this is primarily due to reduced canopy interception. The reduction in evaporation in CP4 is offset by an increase in both surface and sub-surface runoff. Observation-based estimates of terrestrial water fluxes over Africa are rather uncertain, but the high rates of interception loss in R25 appear unrealistic in comparison to both pan-African and local estimates.
The sub-grid rainfall assumption applied to the model's canopy hydrology and infiltration excess runoff calculations (described in Sect. 2) is expected to reduce changes in the partition of rainfall between hydrological pathways at different model resolutions. It is clear from the above discussion that that goal is not achieved in this pair of simulations. We can understand the differences between simulations by considering the differences in rainfall intensities, as shown for a small region in Central Africa, 27.8 to 28.8°E, 4.5 to 5.5°N (Fig. 2a, and repeated in Fig. 2b). This region is selected for having similar mean hourly rainfall in both simulations. Here the hourly rainfall for 9 years in the two simulations have been binned following the logarithm-based formula presented by Klingaman et al. (2017), and shows the contribution of each bin to the hourly mean rainfall. It should be noted that for this analysis of grid box scale processes, we are sampling the models on their native grids. As highlighted in previous analyses of these simulations using a common 25 km grid (Senior et al. 2021), ) CP4 produces a larger fraction of its total rainfall in the higher intensity bins than R25, bringing CP4 much more into line with observations (Senior et al. 2021). The impacts on canopy interception and infiltrationexcess runoff of simply switching from a rainy area fraction of µ = 0.3 (a typical value for R25) to µ = 1 (as prescribed in CP4) are illustrated in Fig. 2c and d. These plots show the functional relationships between the grid box mean rainfall, and interception (rainfall minus throughfall) and infiltration excess runoff, as described in Eqs. 2 and 3. Here we have assumed values of canopy capacity and vegetation infiltration enhancement factor of 0.55 mm and 2 (as used for the grass tiles), a dry canopy, and saturated hydraulic conductivity of 0.019 mm s −1 (as prescribed uniformly in both models). For each bin and each simulations we have then multiplied the simulated rainfall spectra (Figs. 2a, b) by their respective sub-grid rainfall distributions (Figs. 2c, d) to compute the combined impact of differences in rainfall distribution and differences in µ on interception (Fig. 2e) and infiltration-excess runoff (Fig. 2f). Comparing the sum of the bins we estimate the interception to be 0.026 mm/h and 0.017 mm/h for R25 and CP4 respectively for this region. Figure 2e illustrates clearly why both interception and percentage interception loss (interception/rainfall) is higher in R25 (15%) than CP4 (9%). In R25, the high frequency of light rainfall (less than 1 mm hr −1 ) contributes to the larger interception losses, in spite of the assumption of rainfall only occurring over 30% of the grid box. At the larger rainfall rates typical of CP4, the vast majority of rainfall falls through the canopy to the ground. At the pan-African scale, this leads to canopy interception in R25 almost 4 times as large as in CP4 (Table 1). The impact of contrasting rainfall intensities on infiltration excess runoff between the two simulations accounts for a smaller proportion of the water budget than interception. As expected, higher rainfall rates in CP4 yield a higher infiltration excess runoff i.e. 0.004 mm/h compared to 0.001 mm/h in R25. However, the combined impact of interception and infiltration excess runoff means that, CP4 delivers a similar proportion (89%) of its rainfall as infiltration as R25 (85%). In this particular example, infiltration excess runoff is larger in CP4 than R25, but note that this ignores the role of the large scale hydrology scheme that also contributes to the surface runoff reported in Table 1.
Finally we note that CP4 has lower insolation than R25 due to differences in the cloud scheme, and higher, brighter clouds in CP4. This difference is particularly true in the presence of convective rainfall ). However in the next section we do not see a strong impact of this on the role of canopy evaporation in the hours after rain.

Feedbacks on the diurnal cycle via canopy interception
We next consider how the differences in canopy interception between the two simulations affect the atmosphere over the Fig. 2 a and b show the contribution of hourly rainfall in CP4 (blue) and R25 (orange) to average grid box rainfall over a Central Africa sub-region (27.8-28.8° E, 4.5-5.5° N). The same distributions are shown on a and b but the x-axis range is reduced for (a, c and e). The sub-grid rainfall relationships from Eqs. 2 and 3 for c interception and d infiltration excess runoff respectively, under different assumptions of rainfall area fraction corresponding those of to CP4 (µ = 1) and R25 (µ = 0.3). The contribution of each rainfall bin to e interception and f infiltration excess runoff, constructed by weighting the interception and infiltration excess runoff curves by the corresponding model rainfall spectra diurnal cycle. This feedback is rapid as the canopy store is small (typically 0.5 mm), and when full, evaporates at potential rates ~ 3-5 mm day −1 , dominating the energy partition in the hours after rain. Moreover, the diurnal cycles of rainfall differ between CP4 and R25 (Fig. 3). As shown in other studies e.g. Dai (2006) and Guichard et al. (2004) the parameterised convective schemes tend to trigger convection too easily, leading to a model that starts raining earlier in the day. This is true for R25 in which rainfall peaks at 1300 h compared with 1600 h in CP4. The evaporative fraction (EF) is a simple measure of the surface flux partitioning of available energy (total energy at the surface) into latent and sensible heat, defined as the latent flux divided by the sum of sensible and latent fluxes. In Fig. 4 we show how the contrasting diurnal cycles of rainfall between the two simulations, in combination with their distinct canopy interception behaviours, affects the diurnal cycle of EF. Maps of the afternoon minus the morning EF differences for R25 and CP4 for March to May are shown Figs. 4 a and b respectively, though similar features are found in other seasons. In R25, afternoon EF is typically 0.05-0.1 higher than in the morning right across sub-Saharan Africa, with a particularly strong signal over West Africa. By contrast, CP4 has only a weak morning to afternoon increase in EF. The difference between these two fields ( Fig. 4c) demonstrates across tropical North Africa a much greater increase in latent heat flux during the afternoon in R25 at the expense of the sensible heat flux.
We now use the southern West Africa region (− 8 to 6° E, 5 to 10° N, land points only, indicated in Fig. 6) to look at the different model responses over the diurnal cycle, again averaged over the March-May season (Fig. 5). Here EF is computed for daytime hours when total available energy is greater than 100 Wm −2 . The EF is computed on the native model grid and regridded to the 25 km model grid for comparison. Considering the parameterised model first, the rainfall (Fig. 5a) is near zero at 0600 LT but increases through the morning, peaking at 1300 LT. Meanwhile the partitioning of energy (Fig. 5b) switches from being equally apportioned to sensible and latent heat (EF of 0.5 at 1200 LT) to one dominated by latent heat (EF ~ 0.67 at 1500 LT). The dashed line indicates the canopy evaporation component of evaporative fraction, computed as the ratio of canopy latent heat flux to the available energy. Canopy EF rises from 0.11 at 1000 LT to 0.27 at 1500 LT in R25, and is clearly responsible for the overall daytime rise in EF seen in Fig. 5b. In CP4 on the other hand, EF is relatively constant during the day, but rises in late afternoon in response to the peak in rainfall and associated canopy evaporation. In Fig. 5c we see that up to mid-morning, with similar levels of available energy to CP4, R25 produces slightly stronger sensible heat (when the canopy is fairly dry), and has a more rapid decline in the afternoon. That rapid afternoon decline is due to a combination of a wetter canopy (which increases EF) and reduced available energy. The feedback in R25, whereby excessive canopy evaporation and changes in surface insolation suppress afternoon heating and thereby reduces convergence of moisture at low levels, contributes to the reduced amplitude of the diurnal cycle of column moisture flux convergence in R25 (Fig. 5d). Effectively, canopy evaporation in R25 is bringing forward the phase of diurnal heating by 1-2 h, enhanced by additional cloud shading.
We now consider how the rapid afternoon decrease in sensible heat in R25 impacts the atmosphere. In this coastal region, we see differences in the sea breeze circulation between the two simulations (Fig. 6). Suppressed afternoon surface heating in the rainiest coastal areas in R25 reduces the land-sea pressure difference, in turn weakening the sea breeze. Though weaker than at the coast, this pattern of reduced afternoon convergent flow in R25 is evident across the interior of the region. The feedback in R25, whereby excessive canopy evaporation suppresses afternoon heating and thereby reduces convergence of moisture at low levels, contributes to the reduced amplitude of the diurnal cycle of column moisture flux convergence in R25 (Fig. 5d). Better resolution in CP4 of purely atmospheric processes presumably also plays a role in the moisture flux convergence differences between the two simulations (Finney et al. 2019). Overall, the more sustained moisture convergence evident during late afternoon in CP4 will help to maintain rainfall well into the evening in that model (Fig. 5a,d), in better agreement with observations.

Feedbacks on the seasonal time scale via soil moisture
Having identified a strong impact of water budget partitioning on the atmosphere at the diurnal time scale, we now examine its impact at the seasonal time scale, via soil moisture. Soil moisture state integrates the effect of rainfall partitioning and evaporation over a range of timescales. In Fig. 7 we present the annual cycle of the key components of the water balance (rainfall (P), throughfall, evaporation (E), runoff, soil moisture and P-E) as a function of latitude. Both models show the seasonal migration of the rain band which characterises African climate. However south of 20° S CP4 has lower rainfall in all months (Fig. 7c), related to the improved representation of convection in tropical -extratropical cloud bands (Hart et al. 2018). In Equatorial and Northern Africa, CP4 has lower rainfall during months January to April but higher rainfall later in the year. The large differences in interception loss between models means that the pattern in throughfall (Fig. 7d-f) is modulated compared to the rainfall forcing. The surface beneath the canopy receives a larger fraction of rain in CP4 than R25, shifting the features in Fig. 7c towards more positive values in Fig. 7f. This shift towards enhanced water fluxes in CP4 is reinforced when considering surface runoff. Even in latitude-months where CP4 delivers much less rain than R25 (e.g. Equatorial Africa during boreal winter), this difference is not translated into less surface runoff. Compared to R25, enhanced surface runoff in CP4 is due to a combination of reduced interception loss and more efficient runoff generation because of higher rainfall intensities. Similarly, sub-surface runoff is greater in CP4 for all latitude-months. This reflects increased infiltration (throughfall minus surface runoff) and hence soil moisture in CP4. The sub-surface runoff differences (Fig. 7l) lag the surface runoff pattern (Fig. 7i) by 1 to 2 months due to the time taken for rainfall to travel down through the profile to the bottom soil layer (between 1 and 3 m depth). Figure 7m, n, o show the month on month change in soil moisture storage summed across all four model soil levels. This shows the extent to which soil water is replenished by rainfall (throughfall) and its reduction in the dry season by drainage, transpiration and bare soil evaporation. The seasonal soil moisture amplitude in the whole soil column is 50% higher for CP4, 259 mm in the equatorial latitudes, compared to 184 mm in R25. Total evaporation (the sum of evaporation from bare soil, transpiration and canopy evaporation (Fig. 7p, q, r) is markedly lower in CP4 within the climatological rain band. This again highlights the key role of rapid re-evaporation of water intercepted by the canopy in explaining the different land-atmosphere coupling in the two simulations. Only as the rain band retreats southwards in September to December is there evidence of enhanced total evaporation in CP4 over Tropical North Africa (Fig. 7r). This brief post-rainy season effect in CP4 is due principally to an extended period of transpiration, which can be maintained because of increased root zone soil moisture. Overall, the slower evaporative response of CP4 favours a stronger influence of soil moisture (rather than canopy water) on surface fluxes. Finally, Fig. 7s, t, u reveals a strong inter-model difference in P-E, which at monthly time scales is equivalent to atmospheric moisture flux convergence, because change in atmospheric moisture storage om this time scale is close to zero. In the most active phases of the climatological rain band, P-E in CP4 is dramatically different, reaching values of up to 60 mm/month, or up to 2 times higher than R25. For the atmospheric water balance, enhanced moisture convergence into the rainband in CP4 helps to offset the loss of rapidly recycled moisture from the surface. In this pair of simulations, it is more challenging to identify the impact on the atmosphere of root zone soil moisture than canopy water. That is because in R25, the coupling between canopy water and total evaporation is so strong that any analysis of the two models during rainy months will be dominated by interception. Instead we examine the impact of post-rainy season soil moisture anomalies on atmospheric circulation. Such feedbacks are also likely to be operating within the rainy season, and potentially affecting rainfall, though in the case of R25, we would expect the soil moisture pathway to be bypassed by canopy interception. As discussed above, there are strong contrasts in surface fluxes between R25 and CP4 in Tropical North Africa during the onset of the dry season in October and November. We consider the response in November across West Africa when rainfall and interception are negligible and the role of soil moisture is most clearly expressed. Focusing on West Africa, Fig. 8 shows that daytime (0500 to 1700 LT) sensible heat fluxes are typically 20-60 Wm −2 lower in CP4 than R25, equivalent to a decrease of 30-50%. This is driven by enhanced root zone soil moisture in CP4 from increased infiltration in the preceding months. Covering an area ~ 2500 × 1000 km, the negative surface heating anomaly is strong enough to induce a divergent anomaly in the regional scale circulation. This feedback suppresses moisture flux convergence into the region, and would effectively act as a negative soil

Discussion
The recent application of CPMs for climate studies raises important questions about the fate of rainfall once it reaches the ground, and its potential feedback on the atmosphere. The higher rain rates produced by CPMs would be expected to substantially increase surface runoff in comparison to conventional models with parameterised convection via higher infiltration excess runoff. We do see that effect in our analysis, though the increase in total surface runoff when comparing CP4 with R25 (amounting to 8% of total rainfall) is modest compared to the change in high intensity rainfall (with a 25% increase in rainfall above 20 mm h −1 ). The effect of enhanced surface runoff in CP4 is limited for two reasons. Firstly, the sub-grid rainfall distribution assumed by the land surface scheme in R25 allows significant runoff to be generated even at rainfall rates of 10 mm h −1 . Secondly, both simulations have high values of saturated hydraulic conductivity (from the uniformly sandy soils) and this limits infiltration-excess runoff overall. On the other hand, for the sandy soils adopted in these simulations, there is relatively rapid moisture transport down the profile and into sub-surface runoff. Indeed the increase for CP4 (compared to R25) in sub-surface runoff is of a similar magnitude to the increase due to surface runoff. Taken together, total runoff in CP4 appears very high compared to available observational estimates (Table 1).
Relatively speaking, there is a much larger difference in interception loss than surface runoff between models. In this case, the sub-grid rainfall distribution cannot offset the impact of strong biases in rainfall forcing in R25 (i.e. overly frequent light rain). As a result, interception loss is unrealistically high in R25. The additional throughfall simulated by CP4 helps to offset soil moisture losses due to runoff, and allows CP4 to maintain higher rates of transpiration, which extend further into the dry season. Importantly, the effect of reduced wet season evapotranspiration in CP4 does not simply translate to increased sensible heat flux, as CP4 also receives less insolation due to differences in the cloud scheme . This is in contrast to the situation described by (Ban et al. 2014) over Europe, where drier soils in a CPM combined with increased insolation to amplify warm biases found in a coarser resolution version of the model.
Although R25 is capable of producing high intensity rainfall on the model time-step , our study shows that the frequency of low intensity rainfall (< 1 mm day −1 ) remains a dominant control on the overall surface hydrological response. We see a strong impact of interception errors in R25 on the diurnal cycle. The switch to a latent-dominated energy balance in the afternoon means that the model effectively brings forward the phase of maximum surface heating by 1-2 h. This reduces moisture flux convergence in the afternoon, limiting the ability of that model to sustain rainfall into the evening. Similar impacts of surface heating on the diurnal moisture flux over West Africa have been discussed previously by Birch et al. (2014), though there the differences were in association with cloud cover, affecting total available energy, rather than its partition. On the other hand, an impact of differences in interception loss on land-sea pressure differences and convergence has been identified in sub-seasonal simulations over India with an earlier version of this model (Willetts et al. 2017). The change in rainfall frequency and diurnal timing which is often reported when switching off parameterisations of deep convection is expected to affect the surface radiation budget through changes in cloud cover. We found differences in surface insolation here, though it was not possible to isolate the impact of the representation of convection from the effect of using different cloud schemes in the two simulations. Along with the potential for land feedbacks via hydrological pathways that was our focus, CP modellers should also consider radiative land feedbacks, particularly in regions (and models) where evapotranspiration is less sensitive to the frequency and intensity of rainfall.
This study suggests that the assumptions made in the exponential distribution of sub-grid rainfall need to be revisited for both convection-permitting and parameterised configurations of this model. The inability of the sub-grid scheme to overcome the overly frequent light rain in R25 (and presumably also other coarse resolution configurations of the model) suggests a need for further model development. We note that as well as the description of convection, differences in horizontal resolution between the two simulations could also impact the simulated hydrological pathways. Given the fundamental differences in both the rainfall characteristics and the representation of sub-grid hydrology however, we expect resolution effects to be minor. Moreover, our comparisons of the surface and atmospheric water budgets (for example in Fig. 7) are performed over large regions rather than at the grid point scale. The excessive amount of time during the day when the vegetation canopy is wet also has implications for other aspects of the model not explored here. For example, the transpiration pathway is effectively bypassed whenever the canopy is wet. A canopy that remains wetted through overly frequent rainfall can also reduce the sensitivity of interception loss to canopy capacity differences from different land covers. Together these factors may impact simulations of deforestation, and deleteriously impact on plant productivity, which would compromise the carbon cycle simulated in future Earth System type configurations of the model. Our results also raise the question of whether there is still a need for a sub-grid rainfall assumption in a CPM given that, unlike in GCMs, grid box rainfall rates are already large enough to trigger infiltration-excess runoff.

Conclusion
The pioneering 10 year convection permitting simulation over Africa has provided the opportunity to examine the role of the land surface under a wholly different rainfall regime to that in which these models have been developed. Convection-permitting models offer significant improvement on the depiction of rainfall intensities. Our results show this has a significant impact on the surface hydrology, and in turn on how that can feed back on the atmosphere.
We find that averaged across the continent, the CP4 model produces very similar annual rainfall to R25, and to observed values, but with large contrasts in the partitioning of that rain into the different hydrological pathways. Annual ET falls sharply in CP4, amounting to a decrease of 13% of total rain, and taking it below the range of estimates in the literature. The drop in ET is almost entirely due to a reduction in the re-evaporation of rainfall intercepted by the canopy, which in R25 appears unrealistically high when considering observational evidence. With ET suppressed, there is a much greater contribution to rainfall in the atmospheric water budget from moisture flux convergence. This is particularly clear within the rain band that migrates seasonally across the continent.
The drop in ET is almost entirely due to a reduction in the re-evaporation of rainfall intercepted by the canopy, which in R25 appears unrealistically high. Excessive canopy interception in R25 is due to overly-frequent light rain typical of coarse resolution models. Even with a sub-grid rainfall distribution assumed, the surface scheme is unable to adequately address this bias. As a result of excessive canopy interception in R25, in combination with a typical rainfall diurnal cycle bias, sensible heat in R25 weakens rapidly in the afternoon, an effect which feeds back on sea breeze circulations and moisture flux convergence.
Relative to R25, throughfall is increased in the CP4 simulation, and occurs less frequently but at higher rates. This shift in intensity distribution triggers an increase in surface runoff in CP4 equivalent to 9% of total rain. More rain also infiltrates into the soil in CP4, producing a larger amplitude seasonal cycle of soil moisture and increased sub-surface runoff. Evidence indicates that where increased soil moisture storage allows enhanced transpiration during the dry season, the impact of accompanying suppressed sensible heat in CP4 can induce a regional scale divergent flow.
This study shows how, even when total continental rain remains unchanged, switching a convective parameterisation off can substantially impact land surface behaviour, via changes in timing and intensity of rainfall, which in turn can feedback on the atmosphere. With the increasing use of CPMs for climate simulations, we suggest that detailed diagnosis is needed of how existing land parameterisations, usually developed for GCM studies, respond to the more realistic rainfall forcing that CPMs produce. Fig. 7 Zonal mean of modelled water balance components for R25 (left column), CP4 (middle) and model difference (CP4-R25; right). With each component shown in rows: rainfall (P; top row), throughfall (row 2), surface runoff (row 3), sub surface runoff (row 4), change in total column soil moisture (dSM/dt; row 5), total evaporation from the surface (row 6) and P-E (bottom), all units in mm/month. The contours on the difference plots (c, f, i, l, o) indicate model rainfall differences (CP4 minus R25) of ± 5 mm. The contours on the difference plots (r, u) indicate mean modelled rainfall in CP4 of 100 and 200 mm/month ◂ Fig. 8 Model differences (CP4 minus R25) in November mean daytime (0500-1700 LT) sensible heat, vectors show the difference in daytime (0500-1700 LT) 10 m winds