ABSTRACT
In African nations with complex topographies, alternative rainfall estimation methods such as satellites are crucial. This study aimed at predicting the spatial and temporal distribution of rainfall in the lake Tana sub-basin from 1990 to 2020. A satellite-based rainfall estimate of Climate Hazards Group Infrared Precipitation with Station (CHIRPS) data was used with the same spanning period. The validation process employs point-to-pixel analysis, comparing CHIRPS estimates with observed data at specific gauge stations. The findings showed that CHIRPS had well estimated rainfall incidence in the highland areas and significantly overestimated it in the lowland areas. The Mann–Kendall trends for January, June, and August indicate decreasing trends, while the winter and spring seasons show notable declines. Regression analysis reveals a non-significant decrease in annual rainfall with the highest in the summer and relatively dry winters. In addition, the coefficient of variation value of 26.37% suggests a moderate level of variability around the mean annual rainfall. In conclusion, the CHIRPS satellite exhibited varied performance across the Tana sub-basin, with site-specific discrepancies and notable inaccuracies at certain stations. The study underscores the importance of considering local factors and topography in satellite-based rainfall assessments, providing valuable insights for agricultural planning in the region.
HIGHLIGHTS
The variability of rainfall was investigated using the standardized anomaly index, coefficient of variation, and precipitation concentration index.
The percentage of daily rainfall events with high intensity was overestimated while the number of daily rainfall events with light precipitation was underestimated by Climate Hazards Group Infrared Precipitation with Station data.
INTRODUCTION
Rainfall is the most important climate input for agriculture, providing the water resources that sustain crops and livestock (Namara et al. 2010; Xu et al. 2022; Chen et al. 2024). Accurate estimation of rainfall amount and variability is critical for agricultural planning and food security (Kogo et al. 2021; Wang et al. 2022; Zhou et al. 2023). Traditionally, rainfall has been measured by ground-based rain gauge networks (Cimini et al. 2013; He et al. 2021; Guan et al. 2023). However, rain gauges provide only point data and are sparsely distributed in many parts of the world, hampering accurate spatial rainfall estimation (Hasan et al. 2016; Li et al. 2020; Dong et al. 2023). In recent decades, satellite-based rainfall products have emerged as an alternative to augment ground observations and provide better spatial coverage (Brunetti et al. 2021; M. Wu et al. 2022; X. Wu et al. 2022; Dong et al. 2023). One such product is the Climate Hazards Group Infrared Precipitation with Station (CHIRPS) data, which combines thermal infrared (TIR) satellite imagery with station data to create high-resolution rainfall time series over land (Shen et al. 2020; Cheng et al. 2023; Yuan et al. 2023). CHIRPS was selected for this study because it has been widely validated across Africa and shown to perform reasonably well compared with ground observations in capturing rainfall variability and trends (Ayehu et al. 2018; Dinku et al. 2018; Esayas et al. 2019).
Several studies have validated the performance of CHIRPS rainfall estimates across Africa. Dinku et al. (2018) compared CHIRPS with ground observations across all of mainland Africa and found a good correlation (r = 0.70) on average, with higher performance in tropical regions compared with arid regions. However, CHIRPS had a dry bias and underestimated rainfall variability. Toté et al. (2015) and Zhou et al. (2022a) evaluated CHIRPS in three climatologically diverse regions in Africa and found it represented the various rainfall regimes well, with correlation coefficients ranging from 0.51 to 0.81 against stations. However, CHIRPS again showed limitations in capturing heavy rainfall events. In Ethiopia specifically, Dinku et al. (2018) found CHIRPS had a high correlation (r = 0.81) with ground data, but underestimated rainfall, especially in the Ethiopian highlands. By contrast, Worqlul et al. (2018) found CHIRPS overestimated annual rainfall in the Upper Blue Nile basin in Ethiopia compared with stations. Both studies noted challenges with properly capturing topography-induced rainfall in the mountainous regions of Ethiopia. Regarding the Tana sub-basin in particular, no studies were found that specifically validated CHIRPS performance in this region.
This study focuses on the Tana sub-basin in the Ethiopian highlands. The region is important for agriculture, but suffers from lack of ground rainfall data. However, no studies have specifically evaluated the capability of CHIRPS to estimate rainfall in the study area. Therefore, there remain knowledge gaps regarding the accuracy of CHIRPS in this region and its potential to augment sparse ground observations for agricultural planning. The objective of this study is to evaluate the capability of CHIRPS satellite rainfall estimates to accurately estimate and predict rainfall variability, compared with ground-based measurements.
MATERIALS AND METHODS
Description of the study area
Locational map of the study area: (a) Ethiopian regional states, (b) Amhara Regional State and (c) lake Tana sub-basin.
Locational map of the study area: (a) Ethiopian regional states, (b) Amhara Regional State and (c) lake Tana sub-basin.
Methods
Data source
Gauge station rainfall data
The study utilized rainfall data from ten meteorological stations covering the period from 1990 to 2020, sourced from the Ethiopian Meteorological Institute (EMI), namely Delgi, Gorgora, Enfranz, Maksegnit, Addis Zemen, Woreta, Dekestifanos, Zegie, Merawi, and Dangila. They were selected as reference points to assess the precision of the CHIRPS rainfall product across the study region. By employing the rain gauge measurements from these stations, the research aimed to evaluate how accurately the CHIRPS rainfall data reflected the actual observed rainfall patterns in the specified area. The chosen time period spans three decades, allowing for the analysis of long-term trends and variability in rainfall patterns. This duration is significant for capturing potential climate changes or shifts in precipitation that may occur over extended periods.
CHIRPS satellite rainfall estimate
The study employed the CHIRPS satellite rainfall estimate as a key data source. CHIRPS is a comprehensive database developed collaboratively by the US Geological Survey (USGS) and the University of California, Santa Barbara (UCSB). CHIRPS offers a 0.05° resolution gridded rainfall time series on a daily, monthly, and annual basis at a global scale. To assess the predictive accuracy of the satellite rainfall estimate over the study area, monthly and annual CHIRPS data for the past 30 years was downloaded from https://earlywarning.usgs.gov/fews. This extensive temporal coverage is valuable for understanding changes in precipitation patterns over the past three decades.
Methods of data processing and analysis
Validation techniques
The primary focus of the validation work was to assess the performance of CHIRPS in predicting rainfall for data-scarce areas within the lake Tana sub-basin, aiming to enable its utilization in spatial and temporal analyses of rainfall variability within the basin. Conducted over the 1990–2020 period, which featured reliable ground-based observations from the EMI, the validation compared CHIRPS rainfall products with observed data. Given the sparse distribution of observation stations in the lake Tana sub-basin, traditional interpolation of ground-based gauge measurements to a gridded dataset was deemed impractical due to significant uncertainties. Consequently, the evaluation adopted a point-to-pixel analysis, overlooking station locations within pixels and anticipating the introduction of random errors. This approach acknowledged uncertainties arising from rain gauge measurements and discrepancies in the spatial alignment of reference gauges and pixel centers. The assessment differentiated between lowland and highland regions, recognizing potential variations in the performance of CHIRPS estimates based on distinct rainfall characteristics.
To evaluate CHIRPS' ability to detect rainfall events on a daily timescale, basic categorical validation statistics were employed, derived from a contingency table (Table 1) encompassing four event combinations: hits (A), false alarms (B), misses (C), and correct negatives (D). The assessment utilized a rainfall threshold of 1 mm to determine rain or no rain events, a value commonly used in similar regional studies (Ayehu et al. 2018; Dinku et al. 2018; Liu et al. 2023) and representing the smallest nonzero recorded value at gauging stations. The categorical validation metrics included the probability of detection (POD), false alarm ratio (FAR), frequency bias index (FBI), and Heidke skill score (HSS). POD gauges the satellite product's skill in detecting rainfall occurrences, ranging from 0 to a perfect score of 1. FAR assesses false detections, with a perfect score of 0. FBI compares the satellite product's rainfall day detection frequency with gauge measurements, indicating underestimation (<1) or overestimation (>1) by CHIRPS. HSS measures estimation accuracy while considering random chance, ranging from minus infinity to 1, where HSS < 0 signifies random chance superiority, HSS = 0 denotes no skill, and HSS = 1 indicates perfect precipitation estimation by the satellite product (Table 2).
Contingency table for comparing CHIRPS satellite-based rainfall estimates and rain gauge observations (daily rainfall threshold used is 1 mm)
. | CHIRPS ≥ 1 mm . | CHIRPS < 1 mm . |
---|---|---|
Gauge ≥ 1 mm | A | C |
Gauge < 1 mm | B | D |
. | CHIRPS ≥ 1 mm . | CHIRPS < 1 mm . |
---|---|---|
Gauge ≥ 1 mm | A | C |
Gauge < 1 mm | B | D |
Statistical formulas for rainfall detection capabilities of CHIRPS based on categorical metrics
Statistics . | Equation . | Range . | Best value . |
---|---|---|---|
POD | ![]() | 0 to 1 | 1 |
FAR | ![]() | 0 to 1 | 0 |
FBI | ![]() | 0 to ∞ | 1 |
HSS | ![]() | −∞ to 1 | 1 |
Statistics . | Equation . | Range . | Best value . |
---|---|---|---|
POD | ![]() | 0 to 1 | 1 |
FAR | ![]() | 0 to 1 | 0 |
FBI | ![]() | 0 to ∞ | 1 |
HSS | ![]() | −∞ to 1 | 1 |
In assessing the accuracy of the CHIRPS rainfall product in quantifying precipitation, the authors employed four widely recognized error metrics: mean error (ME), mean absolute error (MAE), the Nash–Sutcliffe efficiency (NSE) coefficient, and bias, examining both daily and monthly timeframes (Table 3). Correlation coefficients were omitted from the evaluation due to their sensitivity to extreme values (outliers) and insensitivity to additive and proportional differences between satellite products and gauge measurements, particularly in daily estimations. ME and MAE, with dimensions in mm, provide insights into the average magnitude of estimation errors, aiming for a perfect score of 0. The root mean squared error was replaced by MAE to mitigate the impact of extremely high rainfall values or outliers. NSE, a dimensionless metric ranging from −∞ to 1, gauges the relative magnitude of residuals compared with observed precipitation values, with higher values indicating better agreement between satellite estimates and gauges. Negative values signify the reference mean as a superior estimate, while 0 implies the reference mean equals the satellite estimates. Bias, ranging from 0 to 1, measures the correspondence between estimated and gauge mean values, with 1 denoting a perfect score. Instances of bias > 1 and positive ME values suggest overestimation.
Statistical measures of performance used for analysis based on continuous metrics
Statistics . | Equation . | Range . | Best value . |
---|---|---|---|
ME | ![]() | −∞ to ∞ | 0 |
MAE | ![]() | 0 to ∞ | 0 |
NSE | ![]() | −∞ to 1 | 1 |
Bias | ![]() | 0 to ∞ | 1 |
Statistics . | Equation . | Range . | Best value . |
---|---|---|---|
ME | ![]() | −∞ to ∞ | 0 |
MAE | ![]() | 0 to ∞ | 0 |
NSE | ![]() | −∞ to 1 | 1 |
Bias | ![]() | 0 to ∞ | 1 |
Assessment of spatial–temporal variability and trends of rainfall
Following the validation of the CHIRPS product within the lake Tana sub-basin, the study proceeds to evaluate the spatial and temporal variations in rainfall using the annual and monthly CHIRPS rainfall datasets over a 30-year time period. The analysis was conducted by utilizing ArcMap 10.8.1 and R software, allowing for a comprehensive examination of how rainfall patterns evolved over the specified period. This paper delves into both temporal and spatial aspects, exploring the variability and trends in rainfall within the basin. This helps to provide a nuanced understanding of how precipitation changed across different locations within the lake Tana sub-basin over the three-decade timeframe.
To investigate the long-term evaluation of climatic characteristics through time, several experiments are available. The standardized rainfall anomaly index (SRAI), which examines the frequency and intensity of precipitation, was used to measure the changes in observed climatic trends and variability across time. The Mann–Kendall (MK) trend test was used to detect climate trends in time series data, the coefficient of variation (CV) was used to estimate the seasonal and annual variation of the rainfall patterns, and the standardized precipitation index (SPI) was used to quantify the rainfall deficit for the observed time scale in the study area (Asfaw et al. 2018; Shang et al. 2023).





RESULTS AND DISCUSSION
CHIRPS rainfall detection potential
Statistical indicators such as POD, FAR, FBI, and HSS were calculated for the years 1990–2020 across the Tana sub-basin to evaluate the performance of the CHIRPS satellite rainfall estimate. Ten meteorological stations in the Tana sub-basin were quantified for the pixel-to-pixel comparison of gauge and satellite rainfall (Table 4). The POD, which runs from 0 to a flawless 1, assesses how well the satellite product can identify when it rains. With a POD value score of 1, the CHIRPS satellite can accurately predict the amount of rainfall over the Dangila station. This demonstrates that the rainfall measurement at this station was entirely consistent with the rainfall estimate from the CHIRPS satellite. However, compared with rain gauge measurements of rainfall, the CHIRPS satellite's estimate of rainfall above the Dekestifanos station was 15% off. A miss estimation of rainfall with less than 10% from gauge station measurements was found in CHIRPS satellite rainfall estimates across the remaining stations (Table 4). The FAR statistic, which rates false detections from 0 to 1 with a perfect score of 0, evaluates the effectiveness of the CHIRPS satellite rainfall estimate. The findings showed that CHIRPS satellite rainfall estimates over the Addis Zemen station were inaccurate by 25%, according to FAR values. This demonstrates that the amount of rainfall determined by the CHIRPS satellite has probably detected a false alert with a 25% variation at this site. As a result, it was noted that the CHIRPS satellite's performance over the Addis Zemen station was much worse as compared with that of the other similar stations. With FAR values of 5% and 4%, respectively, the CHIRPS satellite conducts rainfall estimation quite effectively over the Woreta and Gorgora stations. Over the study area of the Tana sub-basin, CHIRPS satellite estimates rainfall over the Zegie and Delgi stations with an equivalent FAR of 9%. A FAR of 11% was experienced by CHIRPS rainfall estimations across the Merawi and Enfranz stations at the same time. This is related to the stations' similar patterns in agro-climatic conditions. This indicates that the CHIRPS satellite estimate false alarm rates are similar at stations in roughly similar agro-climatic zones.
Performance of the CHIRPS satellite rainfall estimate over the Tana sub-basin meteorological stations
Station . | Statistics . | |||
---|---|---|---|---|
POD . | FAR . | FBI . | HSS . | |
Dangila | 1.00 | 0.19 | 1.24 | 0.40 |
Addis Zemen | 0.95 | 0.25 | 1.27 | 0.28 |
Dekestifanos | 0.85 | 0.10 | 0.95 | 0.62 |
Enfranz | 0.93 | 0.11 | 1.05 | 0.73 |
Maksegnit | 0.90 | 0.10 | 1.00 | 0.41 |
Merawi | 0.89 | 0.11 | 1.10 | 0.35 |
Woreta | 0.95 | 0.05 | 0.85 | 0.39 |
Zegie | 0.91 | 0.09 | 1.05 | 0.36 |
Delgi | 0.91 | 0.09 | 1.02 | 0.39 |
Gorgora | 0.96 | 0.04 | 0.86 | 0.40 |
Station . | Statistics . | |||
---|---|---|---|---|
POD . | FAR . | FBI . | HSS . | |
Dangila | 1.00 | 0.19 | 1.24 | 0.40 |
Addis Zemen | 0.95 | 0.25 | 1.27 | 0.28 |
Dekestifanos | 0.85 | 0.10 | 0.95 | 0.62 |
Enfranz | 0.93 | 0.11 | 1.05 | 0.73 |
Maksegnit | 0.90 | 0.10 | 1.00 | 0.41 |
Merawi | 0.89 | 0.11 | 1.10 | 0.35 |
Woreta | 0.95 | 0.05 | 0.85 | 0.39 |
Zegie | 0.91 | 0.09 | 1.05 | 0.36 |
Delgi | 0.91 | 0.09 | 1.02 | 0.39 |
Gorgora | 0.96 | 0.04 | 0.86 | 0.40 |
The FBI statistic, which ranges from 0 to infinity with a perfect score of 1, compares the rainfall detection frequency of the satellite product with that of the gauge measurements and is another statistical method used to evaluate the performance of the CHIRPS satellite. An FBI of less than 1 suggests that CHIRPS has underestimated the number of rainy days, whereas an FBI of more than 1 indicates the opposite. Based on this stance, the FBI gave a perfect score of 1 for the CHIRPS satellite rainfall estimate's detection capabilities over the Maksegnit station. This indicates that gauge measurements and the frequency of rainfall detection by the CHIRPS satellite were in agreement. The CHIRPS satellite has estimated rainfall day frequencies of 1.24 and 1.27, respectively, for the Dangila and Addis Zemen stations on a wide scale (Table 4). With FBI values of 0.85, 0.86, and 0.95, respectively, the CHIRPS satellite underestimated the frequency of rainy days over the stations in Woreta, Gorgora, and Dekestifanos.
While accounting for matches owing to random chances (after rain events discovered by chance have been excluded), the HSS statistic was also calculated to determine the accuracy of the estimates. Its range is from minus infinity to 1. HSS 0 denotes that chance is better than the satellite product, HSS = 0 denotes that the product has no skill and HSS = 1 denotes that the satellite estimate product has accurately predicted the amount of precipitation. This claim is supported by the fact that the CHIRPS satellite has the ability to predict rainfall over the Enfranz station with HSS values of 0.73, or the satellite has the opportunity to estimate rainfall at this station at a random percentage of 73%.
Accuracy of the CHIRPS satellite rainfall estimate over the Tana sub-basin
In comparison with rainfall recorded at the ten gauge stations, the accuracy of CHIRPS satellite rainfall estimates across the Tana sub-basin was assessed. At each gauge station over the Tana sub-basin, statistical metrics including ME, MAE, NSE coefficient, and bias error were used to assess the accuracy of CHIRPS satellite rainfall estimates. The statistical metric of MAE was used to assess the variation of the CHIRPS satellite rainfall estimate from gauge station measurements. To examine the fluctuation of the CHIRPS satellite rainfall estimate over each individual site, gauge-specific validation was carried out, and for each of the validation sites, MAE was calculated. As a result, the rainfall estimates from the CHIRPS satellite varied over each validation site.
The MAE value for each of the ten validation locations across the research region is displayed in Table 5. The best value for MAE between two measurements is zero, and its range is zero to positive infinity. The variance of the CHIRPS satellite rainfall estimate over each validation site increases with increasing MAE from zero. With MAE values of 0.95 and 0.85, respectively, the CHIRPS satellite performs better than other gauge station measurements in its estimation of rainfall over the Woreta and Zegie gauge stations. The CHIRPS satellite rainfall estimate across four places differs significantly from the ground measurements of rainfall obtained from gauge stations. Therefore, CHIRPS satellite rainfall estimates over the stations at Dangila, Maksegnit, Delgi, and Gorgora had high MAE values of 1.43, 1.47, 1.57, and 1.45 mm, respectively, compared with ground measurements.
Statistical validation results of the CHIRPS satellite rainfall estimate
Station . | Statistics . | |||
---|---|---|---|---|
ME . | MAE . | NSE . | BIAS . | |
Dangila | 0.52 | 1.43 | −2.54 | 1.17 |
Addis Zemen | 1.00 | 1.30 | −3.06 | 1.23 |
Dekestifanos | 0.80 | 1.20 | −1.57 | 1.20 |
Enfranz | 0.68 | 1.29 | −4.39 | 1.40 |
Maksegnit | 0.87 | 1.47 | −5.21 | 1.46 |
Merawi | 0.35 | 1.16 | −0.98 | 1.20 |
Woreta | 0.70 | 0.95 | −1.33 | 1.13 |
Zegie | 0.35 | 0.85 | −0.26 | 1.13 |
Delgi | 1.11 | 1.57 | −4.73 | 1.42 |
Gorgora | 1.30 | 1.45 | −5.15 | 1.41 |
Station . | Statistics . | |||
---|---|---|---|---|
ME . | MAE . | NSE . | BIAS . | |
Dangila | 0.52 | 1.43 | −2.54 | 1.17 |
Addis Zemen | 1.00 | 1.30 | −3.06 | 1.23 |
Dekestifanos | 0.80 | 1.20 | −1.57 | 1.20 |
Enfranz | 0.68 | 1.29 | −4.39 | 1.40 |
Maksegnit | 0.87 | 1.47 | −5.21 | 1.46 |
Merawi | 0.35 | 1.16 | −0.98 | 1.20 |
Woreta | 0.70 | 0.95 | −1.33 | 1.13 |
Zegie | 0.35 | 0.85 | −0.26 | 1.13 |
Delgi | 1.11 | 1.57 | −4.73 | 1.42 |
Gorgora | 1.30 | 1.45 | −5.15 | 1.41 |
The NSE is a metric used to assess the performance of rainfall models, with a value of 1 indicating a perfect match between the modeled and observed data. As per the analysis in Table 5, the Merawi and Zegie stations exhibit NSE values greater than −1, specifically −0.98 and −0.26, respectively. This discrepancy can be attributed due to several factors. In the case of the Merawi and Zegie stations, the CHIRPS satellite rainfall estimates might align more closely with the observed data, resulting in lower error variances and higher efficiency scores. This can be attributed to specific local conditions, topography, and microclimates unique to these stations, making the CHIRPS satellite estimates relatively more accurate (Ritter & Muñoz-Carpena 2013; Seong et al. 2022). Conversely, for the Maksegnit and Gorgora stations, the NSE values are notably below 1, indicating poor model performance. The unfavorable performance at these stations could be influenced by diverse local climatic patterns, geographical features, or factors affecting rainfall distribution that are not adequately captured by the CHIRPS satellite model.
The overall unreliability of the CHIRPS satellite rainfall estimation across the Tana sub-basin, as indicated by NSE values, underscores the need for considering site-specific characteristics and potential limitations of the satellite model in capturing complex local variations. The favorable ME of 0.35 for the Merawi and Zegie stations suggests that the CHIRPS satellite tends to slightly overestimate rainfall at these locations, possibly due to specific regional dynamics not fully accounted for in the model. In contrast, the inconsistent performance observed in the Maksegnit and Gorgora stations can be attributed to their distinct climatic conditions, leading to significant discrepancies between satellite estimates and ground measurements. This nuanced understanding provides a scientifically grounded basis for the observed statistical outcomes, emphasizing the importance of considering local context and environmental factors in satellite-based rainfall assessments.
A systematic inaccuracy known as bias results in an estimate of an effect or association that is erroneous. To calculate the bias difference between gauge station measurements of rainfall and the CHIRPS satellite rainfall estimate, the bias error of each validation site was also assessed. All of the validation sites were determined to have positive bias errors. This demonstrates that the rainfall measurement acquired from the gauge station has a lower value than the rainfall estimates from the CHIRPS satellite. The overestimation of rainfall events by the CHIRPS satellite over the Tana sub-basin was blamed for the higher values of the CHIRPS satellite rainfall estimate at each validation site. To ascertain the discrepancy between the gauge station measurement and the rainfall estimate provided by the CHIRPS satellite, the ME of rainfall measurement was also assessed. The discrepancy between the ground observational rainfall estimates and the CHIRPS satellite estimate represents equivalent measurement values, according to an ME value that is near zero, while certain validation sites in the Tana sub-basin showed positive and negative ME values that were overestimated and underestimated by the CHIRPS satellite, Because of this, ME values found in almost every location showed that the ground observation and the CHIRPS satellite estimate of rainfall are very similar. With ME values of 0.35, which are substantially closer to zero than the other station measurements, the CHIRPS satellite specifically delivers better rainfall estimation over the Merawi and Zegie stations.
Cumulative distribution function of CHIRPS rainfall estimates and gauge measurements over the Tana sub-basin
The cumulative distribution functions (CDFs) were created using the daily rainfall data for each of the ten rainfall gauge stations in the research area in order to determine the probabilities of the daily rainfall volumes. However, four target rainfall stations were chosen from the total of ten rainfall stations in order to describe the results of CDFs and their related P–P plots due to the vast quantity of data gathered, in an effort to better comprehend the approach utilized. The CHIRPS satellite rainfall estimate's cumulative probability over a given time period is calculated using the CDF. The probability that a chance observation of a CHIRPS satellite rainfall estimate will be less than or equal to gauge station rainfall observations can be calculated using CDF values. The likelihood that a CHIRPS satellite rainfall estimate is greater than actual rainfall measurement data received from gauge stations, or between two values, is also calculated using CDF information.
CDF distribution of CHIRPS rainfall events over Addis Zemen station of the Tana sub-basin.
CDF distribution of CHIRPS rainfall events over Addis Zemen station of the Tana sub-basin.
CDF distribution of CHIRPS rainfall events over Maksegnit station of the Tana sub-basin.
CDF distribution of CHIRPS rainfall events over Maksegnit station of the Tana sub-basin.
CDF distribution of CHIRPS rainfall events over Dekestifanos station of the Tana sub-basin.
CDF distribution of CHIRPS rainfall events over Dekestifanos station of the Tana sub-basin.
Spatiotemporal variability of rainfall over the Tana sub-basin
Spatial variability of rainfall
Spatial variability of rainfall across each month over the lake Tana sub-basin (1990–2020).
Spatial variability of rainfall across each month over the lake Tana sub-basin (1990–2020).
Examining seasonal variations in rainfall patterns is pivotal for understanding agricultural production dynamics in the Tana sub-basin, celebrated for its substantial agricultural productivity (Fazzini et al. 2015; Behailu et al. 2021). Following the Ethiopian agro-climatology classification, which recognizes winter, autumn, summer, and spring, each spanning three specific months, this study delves into the unique characteristics of these seasons. The winter season, encompassing January, February, and March, is identified as the driest among the three main seasons (Behailu et al. 2021). Autumn, comprising April, May, and June, aligns with the onset of the rainy season and is a key cropping period. The summer months of July, August, and September are anticipated to receive the highest rainfall, crucial for the primary cropping season, during which over 80% of agricultural crop production depends on rainfall (Gebremichael et al. 2014; Wakjira et al. 2021). In contrast, the months of October, November, and December, categorized as spring following the summer season, experience less precipitation compared with the summer season. The northern half of the study region encounters minimal winter rainfall, except on the southern edge of the basin. Conversely, the southern part experiences peak rainfall during the autumn season. Surprisingly, the northern region receives the least rainfall in autumn, potentially linked to the basin's topographic features. Higher elevations correlate with increased rainfall, with the southernmost region of the basin exhibiting both elevated terrain and the highest rainfall levels (Aher et al. 2019).
Seasonal rainfall variability of the CHIRPS satellite rainfall estimate over the Tana sub-basin: (a) winter season, (b) autumn season, (c) summer season, and (d) spring season.
Seasonal rainfall variability of the CHIRPS satellite rainfall estimate over the Tana sub-basin: (a) winter season, (b) autumn season, (c) summer season, and (d) spring season.
Annual rainfall variability of the CHIRPS satellite rainfall estimate over the past 30 years in the case of the lake Tana sub-basin.
Annual rainfall variability of the CHIRPS satellite rainfall estimate over the past 30 years in the case of the lake Tana sub-basin.
Exploring the relationship between CHIRPS rainfall and elevation over the Tana sub-basin
Comparison of the (a) CHIRPS satellite rainfall estimate and (b) elevation over the Tana sub-basin.
Comparison of the (a) CHIRPS satellite rainfall estimate and (b) elevation over the Tana sub-basin.
Scatter plot depicting the relationship between the CHIRPS satellite rainfall estimate and elevation over the study region.
Scatter plot depicting the relationship between the CHIRPS satellite rainfall estimate and elevation over the study region.
SD of the CHIRPS satellite rainfall estimate over the Tana sub-basin (1990–2020).
SD of the CHIRPS satellite rainfall estimate over the Tana sub-basin (1990–2020).
Rainfall trend analysis
Rainfall summary statistic
The CV in the analysis of mean annual rainfall over the last 30 years holds significance as it reflects the relative variability or stability in precipitation patterns within the study area during this period. A CV value of 26.37% suggests a moderate level of variability around the mean annual rainfall. This information is crucial for understanding the consistency or fluctuations in precipitation, which is pivotal for various applications, including agriculture, water resource management, and climate studies (Ogunrinde et al. 2019; Zhu et al. 2022). Moderate variance was reported in the months of February and December, with CV values of 167.87% and 222.89%, respectively. The observed highest variances in rainfall, particularly in February and December, can be attributed to several factors influencing precipitation patterns in the study area. In February, factors such as atmospheric circulation patterns, local topography, and climate anomalies can contribute to the pronounced variability. December, being the beginning of the dry season, experiences irregularities such as sporadic rainfall events and the influence of climate phenomena such as El Niño or La Niña, which can significantly impact precipitation patterns (Supari et al. 2018).
Mann–Kendall test
The MK trend test is a popular approach for detecting rainfall trends in time series data. The MK test is used to find monotonic (increasing or declining) patterns in yearly and seasonal rainfall parameter bases. The MK trend test and Sen's estimator were used to explore the existence of long-term changes in rainfall indices. The MK test detects annual and seasonal trend changes with decreased sensitivity to climate outliers. However, if autocorrelation occurs in the time series data, the MK test result may contain some mistake. To address this issue, the pre-whitening process was used, and there was no significant amount of serial autocorrelation at all lags and with no modifications. The MK test from the Z value and trend from Sen's slope (s) calculation were computed based on monthly and annual rainfall data from 1990 to 2020 in the Tana sub-basin, and the results are shown in Table 6.
Summary statistical and monthly MK test results for Tana sub-basin (1990–2020)
Month . | Mean (mm) . | SD (mm) . | CV (%) . | Sen's (s) . | MK test (P-value) . | MK test (Z) . | Trend . |
---|---|---|---|---|---|---|---|
Jan | 7.589 | 30.886 | 407.008 | −84 | 0.019* | −1.293 | D |
Feb | 12.09 | 20.295 | 167.872 | −39 | 0.648 | −0.456 | D |
March | 36.528 | 43.54 | 119.197 | −313 | 0.0004*** | −3.523 | D |
April | 44.459 | 30.529 | 68.667 | −133 | 0.137 | −1.486 | D |
May | 58.578 | 40.586 | 69.286 | 111 | 0.215 | 1.238 | I |
June | 100.46 | 54.694 | 54.444 | −213 | 0.017** | −2.389 | D |
July | 325.293 | 122.713 | 37.724 | 2.9 | 0.753 | 0.314 | I |
Aug | 334.986 | 93.691 | 27.969 | −1.7 | 0.087** | −1.708 | D |
Sep | 114.025 | 50.833 | 44.58 | −86 | 0.338 | −0.956 | D |
Oct | 34.563 | 39.283 | 113.658 | 3.3 | 0.717 | 0.361 | I |
Nov | 9.546 | 15.068 | 157.837 | 112 | 0.017* | 1.361 | I |
Dec | 10.028 | 22.351 | 222.899 | −79 | 0.297 | −1.041 | D |
Annual | 1,088.147 | 286.971 | 26.372 | −133 | 0.138 | −1.482 | D |
Month . | Mean (mm) . | SD (mm) . | CV (%) . | Sen's (s) . | MK test (P-value) . | MK test (Z) . | Trend . |
---|---|---|---|---|---|---|---|
Jan | 7.589 | 30.886 | 407.008 | −84 | 0.019* | −1.293 | D |
Feb | 12.09 | 20.295 | 167.872 | −39 | 0.648 | −0.456 | D |
March | 36.528 | 43.54 | 119.197 | −313 | 0.0004*** | −3.523 | D |
April | 44.459 | 30.529 | 68.667 | −133 | 0.137 | −1.486 | D |
May | 58.578 | 40.586 | 69.286 | 111 | 0.215 | 1.238 | I |
June | 100.46 | 54.694 | 54.444 | −213 | 0.017** | −2.389 | D |
July | 325.293 | 122.713 | 37.724 | 2.9 | 0.753 | 0.314 | I |
Aug | 334.986 | 93.691 | 27.969 | −1.7 | 0.087** | −1.708 | D |
Sep | 114.025 | 50.833 | 44.58 | −86 | 0.338 | −0.956 | D |
Oct | 34.563 | 39.283 | 113.658 | 3.3 | 0.717 | 0.361 | I |
Nov | 9.546 | 15.068 | 157.837 | 112 | 0.017* | 1.361 | I |
Dec | 10.028 | 22.351 | 222.899 | −79 | 0.297 | −1.041 | D |
Annual | 1,088.147 | 286.971 | 26.372 | −133 | 0.138 | −1.482 | D |
Note: *, **, *** significant at 0.1, 0.05, 0.01 levels, respectively; I = increasing trend, D = decreasing trend.
The MK trend test revealed a statistically significant decreasing trend for the months of January, March, June, and August at 10%, 1%, 5%, and 5% levels of significance, respectively (Table 6), while there was a considerable upward trend in November. In October and July, there was no significant increase in rainfall, whereas there was no significant decrease in February, April, September, and December. In terms of rainfall volumes, the monthly share of rainfall showed a 61.54% negative trend and a 39.33% positive trend. Winter and spring, on the other hand, showed a decreasing tendency at a 5% and 1% level of significance, respectively. The geographical trend analysis for the study region likewise reveals a consistent pattern with the area-averaged trend analysis, which is a lowering tendency of the winter and spring rainfall in the region. Annual rainfall in summer and autumn has decreased in a non-significant way.
Seasonal rainfall variability is illustrated in Tables 7and 8 for the winter, autumn, summer, and spring seasons. The seasonal rainfall variability trend for the winter, autumn, summer, and spring seasons was −0.431, −0.676, −1.123, and 0.043 mm/year, respectively. The negative and positive numbers indicate, respectively, a declining and increasing trend of rainfall distribution over the measured period. The autumn rainfall pattern shows a decreasing trend but is not statistically significant with p-value 0.1567, which is greater than the significance level (0.05). The autumn season's rainfall is critical for farmers since it impacts their preparation and planting operations in the studied area.
Seasonal MK test results for Tana sub-basin (1990–2020)
Seasons . | Mean (mm) . | SD (mm) . | CV (%) . | Sen's (s) . | MK test (P-value) . | MK test (Z) . |
---|---|---|---|---|---|---|
Winter | 9.902 | 2.253 | 22.754 | −55 | 0.5521 | −0.0059*** |
Autumn | 46.522 | 11.169 | 24.008 | −127 | 0.1567 | −1.4164 |
Summer | 253.580 | 132.694 | 52.328 | −106 | 0.238 | −1.1797 |
Spring | 52.711 | 54.553 | 103.493 | −3 | 0.982 | −0.0224* |
Seasons . | Mean (mm) . | SD (mm) . | CV (%) . | Sen's (s) . | MK test (P-value) . | MK test (Z) . |
---|---|---|---|---|---|---|
Winter | 9.902 | 2.253 | 22.754 | −55 | 0.5521 | −0.0059*** |
Autumn | 46.522 | 11.169 | 24.008 | −127 | 0.1567 | −1.4164 |
Summer | 253.580 | 132.694 | 52.328 | −106 | 0.238 | −1.1797 |
Spring | 52.711 | 54.553 | 103.493 | −3 | 0.982 | −0.0224* |
Note: *, *** significant at 0.1, 0.01 levels, respectively.
Linear regression results of seasonal rainfall over Tana sub-basin (1990–2020)
Seasons . | Mean (mm) . | SD (mm) . | CV (%) . | R2 . | MK test (P-value) . | Change in rainfall (mm/year) . |
---|---|---|---|---|---|---|
Winter | 9.902 | 2.253 | 22.754 | 0.088 | 0.5521 | −0.431 |
Autumn | 46.522 | 11.169 | 24.008 | 0.083 | 0.1567 | −0.676 |
Summer | 253.580 | 132.694 | 52.328 | 0.037 | 0.238 | −1.123 |
Spring | 52.711 | 54.553 | 103.493 | 0.005 | 0.982 | 0.043 |
Seasons . | Mean (mm) . | SD (mm) . | CV (%) . | R2 . | MK test (P-value) . | Change in rainfall (mm/year) . |
---|---|---|---|---|---|---|
Winter | 9.902 | 2.253 | 22.754 | 0.088 | 0.5521 | −0.431 |
Autumn | 46.522 | 11.169 | 24.008 | 0.083 | 0.1567 | −0.676 |
Summer | 253.580 | 132.694 | 52.328 | 0.037 | 0.238 | −1.123 |
Spring | 52.711 | 54.553 | 103.493 | 0.005 | 0.982 | 0.043 |
Seasonal rainfall variability trend over the study region (1990–2020).
Rainfall anomaly index
Estimated RAI with their respective frequency of years over the past 30 years (1990–2020)
RAI class description . | Frequency . | Years . |
---|---|---|
Near-normal | 7 | 1990, 1991, 1992, 1996, 1997, 2007 and 2009 |
Normal | 2 | 2008 and 2020 |
Moderately wet | 3 | 1993, 1998 and 2006 |
Moderately dry | 8 | 2000, 2001, 2002, 2003, 2004, 2005, 2006 and 2011 |
Slightly dry | 3 | 2012, 2015 and 2017 |
Slightly wet | 7 | 1994, 1995, 1999, 2010, 2013, 2014 and 2016 |
RAI class description . | Frequency . | Years . |
---|---|---|
Near-normal | 7 | 1990, 1991, 1992, 1996, 1997, 2007 and 2009 |
Normal | 2 | 2008 and 2020 |
Moderately wet | 3 | 1993, 1998 and 2006 |
Moderately dry | 8 | 2000, 2001, 2002, 2003, 2004, 2005, 2006 and 2011 |
Slightly dry | 3 | 2012, 2015 and 2017 |
Slightly wet | 7 | 1994, 1995, 1999, 2010, 2013, 2014 and 2016 |
Frequency of RAI distribution with index class description in the study region (1990–2020).
Frequency of RAI distribution with index class description in the study region (1990–2020).
Monthly rainfall anomaly index calculated for the study area (1990–2020).
Precipitation concentration index
The study area PCI was also generated to assess rainfall variability across the study period (1990–2020). The study period was divided into three classes to allow for a clearer evaluation of the PCI trend over the study region. The first designation was an almost ten-year study period examination spanning the years 1990 to 2000. The second group was a ten-year study period description beginning in 2000 and ending in 2010. The final designations were a ten-year study period description beginning in 2010 and ending in 2020. As shown in Figure 17, the PCI values showed an up-and-down trend for ten consecutive years from 1990 to 2000. It dropped from 1990 to 1997 with PCI values of 21.31–15.19 within seven consecutive years.
Meanwhile, the lower the PCI values of the rainfall pattern, the lower the annual rainfall fluctuations will be. This argument is supported by the fact that rainfall fluctuation has steadily decreased from 1990 to 2020, reflecting the lowering tendencies of yearly rainfall fluctuation over 30 years. This rainfall fluctuation significantly affected the vegetation condition in the study region and this situation was confirmed with field observations from the ground. The rainfall distribution pattern was somewhat lowered with an up-and-down trend from 1990 to 2020, which was mainly linked to the state of vegetation condition. In summary, the PCI values of the research region demonstrated that the rainfall pattern distribution was aberrant. The annual rainfall fluctuation in the study area can be said to be particularly high throughout the past 30 years of investigation.
DISCUSSIONS
The robust capability of satellite rainfall estimates to detect precipitation is evident through a synergistic combination of high POD, FBI, and HSS, along with a low FAR. In the Tana sub-basin, the CHIRPS rainfall product tends to overestimate the frequency of rainy days at the Dangila station (POD = 1.00) and underestimate the frequency at the Addis Zemen gauge station (POD = 0.95, FAR = 0.25, and FBI = 1.27). Our investigation demonstrates the commendable performance of the CHIRPS product, especially when compared with a previous study in the surrounding area (Fenta et al. 2018), where POD ranged between 0.54 and 0.6 and FBI ranged between 0.70 and 0.87.
CHIRPS utilizes a fixed TIR cloud coldness duration threshold of 235 K, without temporal or spatial fluctuations. This fixed threshold may contribute to the slight underestimation of rainy days in highland regions, where warm-rain mechanisms related to the warm orographic rain process dominate, rendering the threshold relatively cold. Conversely, overestimation of rainy days in lowland regions may be attributed to sub-cloud evaporation, where the dry atmosphere beneath the clouds evaporates before reaching the surface, leading to an overestimation of rainfall (Dinku et al. 2018; Areerachakul et al. 2022).
In the study region spanning from 1990 to 2020, the annual rainfall exhibited a smaller CV (26.37%) compared with the interannual variability in rainfall. Likewise, there was reduced variability in rainfall amounts during the months of July, August, and September when contrasted with other months. This concurs with recent scholarly investigations highlighting the notable diversity and unpredictability in Ethiopia's seasonal and annual climate patterns (Gebremichael et al. 2014; Gizaw & Gan 2017; Mesfin et al. 2020). Standard rainfall anomalies, reflecting both positive and negative deviations, signify the interannual variability in rainfall over the observed time series. The year 2001 recorded the highest positive anomaly (+2.74), while 2015 marked the most substantial negative anomaly (−3.62). The adverse anomaly in 2015 was predominantly attributed to El Niño, causing detrimental impacts on the primary means of subsistence for rural communities in various regions of Ethiopia (Haile et al. 2021).
Recent studies underscore that Ethiopia has undergone approximately 12 significant historical droughts, influencing the country's economic development (Liou & Mulualem 2019; Kassaye et al. 2021). The observed negative influence of Ethiopian drought events on the livelihoods of millions of people aligns with the findings in our study, particularly in the Tana sub-basin region, from 1990 to 2020. The documented severe droughts in Ethiopia over the past three decades, connected to El Niño–Southern Oscillation occurrences (Zeleke et al. 2017; M. Wu et al. 2022; X. Wu et al. 2022), resonate with our examination of annual rainfall anomalies, revealing significant negative anomalies in 2015, a year associated with El Niño. Our study identifies drought years coinciding with El Niño episodes, including 2015, which had the largest unfavorable anomaly. This concurs with recent studies (Gleixner et al. 2017) that emphasize the effects of El Niño events on rainfall distribution in Ethiopia.
CONCLUSION
This study focuses on a comprehensive evaluation of CHIRPS satellite rainfall estimates in the Tana sub-basin by employing statistical indicators such as POD, FAR, FBI, and HSS. The analysis of CHIRPS estimates against gauge station measurements revealed fluctuations in ME, MAE, NSE, and bias error at different stations. This study emphasized the importance of understanding monthly and seasonal rainfall patterns, including the impact of elevation on precipitation distribution. Trend analyses via the MK test uncovered significant decreasing trends in specific months and seasons. RAI and PCI assessments indicated distinct patterns of moderately dry conditions and changing trends in rainfall concentration over the past three decades. Overall, the study provides crucial insights into the performance and implications of CHIRPS estimates, spatial–temporal rainfall dynamics, and long-term trends in the Tana sub-basin, influencing various sectors such as agriculture, water resource management, and climate studies.
ACKNOWLEDGEMENTS
The authors acknowledge Department of Earth Science, College of Natural and Computational Science, Wollega University Nekemte Campus; Department of Geography and Environmental Studies, College of Social Science and Humanities, Arsi University; Department of Forestry, College of Natural Resource and Agricultural Economics, Mettu University Bedele Campus, Bedele, Ethiopia; Department of Geography and Environmental Studies, Raya University, College of Natural Resource and Agricultural Economics; Department of Forestry and Jimma University College of Agriculture and Veterinary Medicine for the existing facilities to conduct this research.
CONTRIBUTIONS OF AUTHORS
S.S.D. participated in research design, data collection, and data analysis and draft manuscript. Z.R.R. and M.B.M. participated in methodology. K.T.D., H.H.H., and D.O.G. worked on literature and edited the manuscript to the journal style. All authors read and approved the final manuscript.
DATA AVAILABILITY STATEMENT
All relevant data are included in the paper or its Supplementary Information.
CONFLICT OF INTEREST
The authors declare there is no conflict.