the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
The Impact of Convection-Permitting Rainfall on the Dryland Water Balance
Abstract. In drylands, rainfall is typically delivered during short-lived and localised convective storms, the characteristics of which determine how water is partitioned into different terrestrial stores. However, rainfall datasets used in hydrological modelling and assessments of water resources are typically derived from climate models that are too coarse to represent convective processes occurring at scales smaller than the model. In this paper, we quantify the impact of climate model representation of convection on the simulated water balance at four locations in the Horn of Africa: a humid site in the Ethiopian Highlands, a semi-arid site in southern Kenya, an arid site in eastern Ethiopia, and a hyper-arid site in northern Somalia. We benchmark the novel pan-Africa convection-permitting climate model (CP4A) and its parameterised counterpart (P25) against high-resolution satellite-derived gridded datasets of rainfall (IMERG) and PET (hPET). The comparison shows that explicitly resolving convection improves the characterisation of rainfall frequency, intensity, and the relative contribution of low vs high-intensity rainfall to annual totals. We also demonstrate that the representation of convection can impact model PET. However, differences in PET between CP4A and P25 are more muted relative to rainfall, and both can capture seasonal and diurnal PET dynamics. To establish how climate model representation of convection can impact hydrology, we then ran a series of one-dimensional hydrological model experiments along an aridity gradient across the Horn of Africa using Hydrus 1-D, where at each of our four sites, Hydrus was driven by rainfall and PET from CP4A and P25 (and hPET). The ‘drizzle’ bias in P25 means that when rainfall is propagated through Hydrus, wetting fronts are more restricted to upper soil layers, resulting in higher evaporative losses, lower soil moisture, and bottom drainage in drylands. While at our humid site in the Ethiopian Highlands, there are minimal differences in hydrological outcomes; in drylands, the more intense and intermittent rainfall in CP4A means surface runoff is up to ten times higher and bottom drainage up to 25 times higher. We conclude that dryland hydrology is highly sensitive to climate model representation of convection and that forcing hydrological model projections with convectional climate models that parameterise the average effects of convection risks underestimating future crop health, groundwater availability, or flood risk.
- Preprint
(3344 KB) - Metadata XML
- BibTeX
- EndNote
Status: final response (author comments only)
-
RC1: 'Comment on egusphere-2025-1154', Bo Huang, 16 May 2025
Reviews of Manuscript No.: egusphere-2025-1154
Title: The Impact of Convection-Permitting Rainfall on the Dryland Water Balance
Author(s): George Blake et al.
Overall conclusions and recommendations:
This study investigates how convection-permitting climate models (CPMs) influence the representation of rainfall characteristics and subsequent hydrological partitioning in the Horn of Africa (HOA), compared to parameterized models. Using the CP4A (4.5 km resolution) and P25 (25 km resolution) models, the authors analyze rainfall intensity, dry spell duration, and extreme events, then propagate these outputs through Hydrus 1-D to assess impacts on soil moisture, runoff, and drainage. Key findings include improved rainfall characterization in CP4A, higher surface runoff and drainage in drylands, and sensitivity of hydrological outcomes to convection representation.
This study falls in the scope of the journal ‘HESS’. This manuscript is well organized and the topic and conclusions are likely to attract the interest of readers working in areas of climate and hydrology scientists focus on regional climate and water resource. It has the potential to be journal publish. However, there are several critical and minor issues that need to be addressed.
Specific comments:
- The authors use the FAO Penman-Monteith method (Section 2.2) to calculate potential evapotranspiration (PET) and list seven atmospheric variables. However, the equation itself and the role of these variables in its application are not explicitly provided. Additionally, Penman-Monteith method uses wind speed at 2m height but this manuscript does not clarify how wind speed measured at 10m height is adjusted to 2m. Please include the equation, explain the variables, and describe the methodology for converting 10m wind speed to 2m (e.g., logarithmic wind profile adjustment or FAO-recommended constants).
- The climatological analysis appears to aggregate data across the entire study period. Are there any differences between wet and dry seasons? Given the region’s likely seasonal contrasts (Figure C1), will the results/signal change in wet and dry seasons?
- The aridity index (AI) is referenced repeatedly (Lines 37, 226, 527), but its definition (e.g., ratio of precipitation to PET or another formula) is not provided. Please clarify the specific equation or source used for calculating AI to ensure reproducibility and reader comprehension. Or directly cite the reference to indicate the four regions other than use value of aridity index.
- Figures 7 and 9 present soil moisture distributions, but the frequency ranges differ significantly between the two. What explains this discrepancy? Additionally, Figure 9 is described as including green and red dashed lines in the caption, but these are not visible in the figure as provided. Please resolve this inconsistency.
- Figure 10 is labeled as depicting rainfall and PET but does not include these variables in the presented panels. Please revise the caption to accurately reflect the displayed data.
- The symbol “θs” is used extensively without definition, which may confuse readers unfamiliar with soil moisture notation. Please explicitly define “θs” or use a reader friendly abbreviation.
- Duplicate scenario name in Table A1, consider to re-name these scenarios.
Citation: https://doi.org/10.5194/egusphere-2025-1154-RC1 -
AC1: 'Reply on RC1', George Blake, 06 Jun 2025
We would like to thank Dr Bo Huang for his helpful comments on our manuscript. Below please find responses to each comment, which we hope help clarify some of the points raised.
Here we will provide responses to points 4-7 raised by Dr Huang, points 1-3 are addressed in the attached supplement as the internal HESS response platform does not allow for easy presentation of equations needed to address the issues raised.
4. Figures 7 and 9 present soil moisture distributions, but the frequency ranges differ significantly between the two. What explains this discrepancy? Additionally, Figure 9 is described as including green and red dashed lines in the caption, but these are not visible in the figure as provided. Please resolve this inconsistency.
Discrepancies between soil moisture distributions reflect soil properties at each site, total rainfall, rainfall characteristics, and PET. For example, at our more arid sites (Sites A and HA) soil moisture is limited to 30-50% saturation, whereas Sites SA and HU see soil moisture saturations reach greater than 80%. This is likely due to the higher total rainfall and lower PET at these sites compared to our more arid locations. However, there are also differences in soil moisture distributions between Site A and Site HA, with soil moisture on average being higher at our hyper-arid site. This may appear counter-intuitive given these is our most arid location, but soil moisture is also sensitive to root water uptake. At Site A shrubs are present throughout the profile, and take up 12-21% of infiltration, while Site HA is bare soil. We have also shown that soil moisture is sensitive to rainfall characteristics, at Site HA more rainfall is delivered via heavy events (30-82% vs 12-49%) and the magnitude of extreme (99th percentile) rainfall is also higher (17.9 vs 11.2 mm/h). You can see the deeper penetration of wetting fronts at Site HA in fig 8.
Apologies - the caption of Figure 9 needs updating. The correct caption should read: Soil Moisture Distributions with Wilting Points. Modelled distribution of soil moisture at 1.2 mbgl at Site SA (a) and A (b) using P25 (blue) and CP4A (red) rainfall and PET. The dashed orange line represents the wilting point range for Acacia shrubs, based on taking the upper and lower Feddes’ parameters given in Table A2 (wilting point = P2H) (Sela et al, 2015).
5. Figure 10 is labeled as depicting rainfall and PET but does not include these variables in the presented panels. Please revise the caption to accurately reflect the displayed data.Figure 10 is labeled as depicting rainfall and PET but does not include these variables in the presented panels. Please revise the caption to accurately reflect the displayed data.
We apologise for this error. The plot previously included rainfall and PET but they were removed as it made visualising the other variables more challenging. Cumulative rainfall and PET are given in Table 1. Additionally, infiltration largely follows cumulative rainfall, with differences between P25 and CP4A Hydrus simulations largely due to the infiltration partitioning between different stores. If it would improve clarity, we can add rainfall and PET to Figure 10.
If rainfall and PET is removed, the Figure 10 caption will read: Cumulative Infiltration Runoff, Evaporation, Transpiration, and Deep Infiltration. Modelled components of the water balance using CP4A (solid lines) and P25 (dashed lines) rainfall/PET as input for Hydrus 1-D. Plots show the 500 infiltration (a), surface run-off (b), bottom drainage (c), evapotranspiration (d), evaporation (e), and transpiration (f) at our semi-arid, arid, and hyper-arid locations.
6. The symbol “θs” is used extensively without definition, which may confuse readers unfamiliar with soil moisture notation. Please explicitly define “θs” or use a reader friendly abbreviation.
On first use in the manuscript θs is defined – please see line 419: Fig. 7a-d show histograms of depth-integrated θs (% soil saturation in entire 300 cm of the soil profile). But to increase clarity we can more explicitly define θs in 2.4.3 (Hydrological Model Set Up, Data, and Sensitivity). The revised text could read: Throughout this manuscript we will express soil moisture as a saturation percentage (given the symbol θs), which reflects the proportion of pore spaces that are filled with water relative to if all pore space is saturated (eg 100% means all pore space is filled with water, 0% mean all pores are filled with air). Pore space differs by soil properties hence a relative metric is used to compare our hydrological sites.
7. Duplicate scenario name in Table A1, consider to re-name these scenarios.
Table A1 details the soil parameters used in our Hydrus simulations, to establish whether relative biases in hydrological outcomes when driving Hydrus using CP4A/P25 remained regardless of soil parameters, we utilised three soil parameters sets (as stated in section 2.4.3 and Appendix A).
These parameters were based calculated using the Genuchten-Mualem (Van Genuchten, 1980) equations with soil texture values taken from the iSDAsoil database (Hengl et al, 2021). iSDA provides a lower and upper bound of sand, silt, and clay percentages, our default (‘def’) parameter set was estimated using the mid-point of these lower and upper bounds. To create our low (‘lowK’) and high hydraulic conductivity (‘highK’) soil parameters, we used the lower and upper bound of the sand percentage respectively. We then proportionally adjusted the silt and clay percentage to ensure values equalled 100. We used the same labelling for each site as they follow the same methodology and are designed to be comparable across sites. Of course, soil parameters will differ based on the relative proportion of sand, silt, and clay at each site.
References
Allen, R. G., Pereira, L. S., Raes, D. & Smith, M. Crop evapotranspiration-Guidelines for computing crop water requirements-FAO Irrigation and drainage paper 56. Fao, Rome 300, D05109 (1998).
Tetens, O. Uber einige meteorologische Begriffe. Z. geophys 6, 297–309 (1930).
Singer, M.B., Asfaw, D.T., Rosolem, R., Cuthbert, M.O., Miralles, D.G., MacLeod, D., Quichimbo, E.A. and Michaelides, K., 2021. Hourly potential evapotranspiration at 0.1 resolution for the global land surface from 1981-present. Scientific Data, 8(1), p.224.
Zomer, R., Trabucco, A., van Straaten, O., and Bossio, D.: Carbon, land and water: A global analysis of the hydrologic dimensions of climate change mitigation through afforestation/reforestation, IWMI, Vol. 101, 2007.
Mirzabaev, A., et al.: Desertification, in: Climate Change and Land: an IPCC special report on climate change, desertification, land 880 degradation, sustainable land management, food security, and greenhouse gas fluxes in terrestrial ecosystems, IPCC, 2019.
Hijmans, R.J., Cameron, S.E., Parra, J.L., Jones, P.G. and Jarvis, A., 2005. Very high resolution interpolated climate surfaces for global land areas. International Journal of Climatology, 25(15), pp.1965-1978.
-
RC2: 'Comment on egusphere-2025-1154', Federico Gómez-Delgado, 05 Jun 2025
I have included my comments and suggestions for revision in my Referee Report.
-
AC2: 'Reply on RC2', George Blake, 21 Jun 2025
We would like to thank Dr Gómez-Delgado for his very detailed and thoughtful review. His comments are very helpful and we are glad to have the opportunity to elaborate on some of the issues raised. Please find attached our response to his very detailed review. We are more than happy to provide any further clarifications or details where required.
-
AC2: 'Reply on RC2', George Blake, 21 Jun 2025
Viewed
HTML | XML | Total | BibTeX | EndNote | |
---|---|---|---|---|---|
217 | 49 | 15 | 281 | 9 | 15 |
- HTML: 217
- PDF: 49
- XML: 15
- Total: 281
- BibTeX: 9
- EndNote: 15
Viewed (geographical distribution)
Country | # | Views | % |
---|
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1