Next Article in Journal
The Impact Evaluation of Acid Mine Drainage on Zebrafish (Danio rerio) and Water Fleas (Daphnia magna) in the Vicinity of the Geum River Basin in Korea
Next Article in Special Issue
Hydrological Analysis of Batu Dam, Malaysia in the Urban Area: Flood and Failure Analysis Preparing for Climate Change
Previous Article in Journal
The Big Five Personality Traits as Predictors of Negative Emotional States in University Students in Taiwan
Previous Article in Special Issue
Vegetation Changing Patterns and Its Sensitivity to Climate Variability across Seven Major Watersheds in China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Meteorological Drought Variability and Its Impact on Wheat Yields across South Africa

1
Institute of Water and Environmental Management, Faculty of Agricultural and Food Sciences and Environmental Management, University of Debrecen, Böszörményi 138, 4032 Debrecen, Hungary
2
Department of Civil Engineering, College of Engineering, King Saud University, P.O. Box 800, Riyadh 11421, Saudi Arabia
3
School of Geographical Sciences, Nanjing University of Information Science and Technology, Nanjing 210044, China
4
Department of Civil, Geo and Environmental Engineering, Technical University of Munich, 80333 Munich, Germany
5
Institutes for Agricultural Research and Educational Farm, University of Debrecen, Böszörményi 138, 4032 Debrecen, Hungary
6
Institute of Land Use, Technical and Precision Technology, Faculty of Agricultural and Food Sciences and Environmental Management, University of Debrecen, Böszörményi 138, 4032 Debrecen, Hungary
7
Department of Geography, The Islamia University of Bahawalpur, Bahawalpur 63100, Pakistan
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2022, 19(24), 16469; https://doi.org/10.3390/ijerph192416469
Submission received: 30 October 2022 / Revised: 28 November 2022 / Accepted: 6 December 2022 / Published: 8 December 2022
(This article belongs to the Special Issue Hydrological Responses to Climate Change)

Abstract

:
Drought is one of the natural hazards that have negatively affected the agricultural sector worldwide. The aims of this study were to track drought characteristics (duration (DD), severity (DS), and frequency (DF)) in South Africa between 2002 and 2021 and to evaluate its impact on wheat production. Climate data were collected from the South African Weather Service (SAWS) along with wheat yield data from the Department of Agriculture, Forestry and Fisheries (2002–2021). The standard precipitation index (SPI) was calculated on 3-, 6-, 9-, and 12-month time scales, and the trend was then tracked using the Mann–Kendall (MK) test. To signify the climatic effects on crop yield, the standardized yield residual series (SYRS) was computed along with the crop-drought resilience factor (CR) on a provincial scale (2002–2021). The output of the SPI analysis for 32 stations covering all of South Africa indicates a drought tendency across the country. On a regional scale, western coastal provinces (WES-C and NR-C) have been more vulnerable to meteorological droughts over the past 20 years. Positive correlation results between SYRS and wheat yield indicate that the WES-C province was highly influenced by drought during all stages of wheat growth (Apr–Nov). Historical drought spells in 2003, 2009, and 2010 with low CR = 0.64 caused the province to be highly impacted by the negative impacts of droughts on yield loss. Overall, drought events have historically impacted the western part of the country and dominated in the coastal area. Thus, mitigation plans should be commenced, and priority should be given to this region. These findings can assist policymakers in budgeting for irrigation demand in rainfed agricultural regions.

1. Introduction

Climate change has been observed globally [1]. According to the Intergovernmental Panel on Climate Change (IPCC), the universal mean surface temperature is predicted to rise by 1.5–2 °C by 2100 [2,3]. The atmospheric carbon dioxide (CO2) concentration is expected to rise to the range of 540 to 970 parts per million around the 2100 period [4]. When the concentration levels of greenhouse gases such as carbon dioxide, methane, water vapor, and nitrogen oxide rise in the atmosphere, the Earth’s surface temperature increases [5].
This rapid increase in the Earth’s surface temperature has altered the hydrological cycle, leading to extreme weather events, such as drought, floods, heat weaves, and many others. Among them, drought has been ranked as the most damaging weather-induced disaster [6] and is considered very costly [7]. Recently, drought events have rapidly increased worldwide [8], with a catastrophic impact on ecosystems and the welfare of human beings [6]. For instance, drought killed nearly 1.3 million people between 1967 and 1991. Between 1960 and 2016, drought events killed approximately 2.2 million people and affected over 2.6 billion people worldwide. During that period (1960–2016), 669 drought events were recorded, which created USD 146 billion in economic damage [9].
Drought can be grouped into meteorological, streamflow or groundwater, and socio-economic perspectives [10]. To date, more than 100 drought indices have been suggested to analyze and characterize drought events (i.e., intensity, duration, and frequency) [8]. These indices vary in structure, computation, input, and complexity.
The most prominent of these include the standardized precipitation index (SPI) [11], reconnaissance drought index (RDI) [12], standardized precipitation evapotranspiration index (SPEI) [13], Palmer drought severity index (PDSI) [14], crop moisture index (CMI) [15], etc. The PDSI is a landmark of hydrological drought monitoring and incorporates precipitation, evaporation, runoff, and soil water recharge into the soil water balance equation. However, its complex computation and large number of input variables makes it less readily available for use [14,15,16]. The SPI is a more generalized and widely applied index for meteorological drought monitoring based on only one input, i.e., precipitation [16,17]. The performance of the SPI was further enhanced by developing the SPEI and RDI, which consider the temperature and evapotranspiration and are used for agricultural drought monitoring [18].
The World Meteorological Organization (WMO) recommended the use of the SPI in the monitoring of meteorological droughts [19]. Some other studies also used the SPI in combination with other indices to monitor droughts on global and regional scales [20,21], and it has been found to be highly correlated with the SPEI as compared to other Palmer’s indices [16]. Moreover, the versatility of the SPI to compare drought events in multiple climatic conditions at varying temporal scales makes it advantageous to use over other meteorological indices, and it is considered to be well suited for drought-related decision-making processes [22,23]. Hence, the SPI is designed to measure the rainfall deficit from the short term, i.e., SPI-1, SPI-3, and SPI-6, to the long term, i.e., SPI-12, SPI-36, and SPI-48, with negative values depicting drought and positive ones depicting moist conditions. However, the short-term SPI has been found to be well suited to address meteorological and agricultural droughts, while the long-term SPI better examines hydrological stress periods [24,25]. Thus, it has been reportedly applied to detect multiscale drought events in Eswatini [26], Cameroon [27], Hungary [8,28], Syria [19], and China [29].
Africa is prone to the adverse impacts of climatic extremes such as droughts because the continent is composed of developing economies with no proper infrastructure or adequate mitigation plans. Drought makes up the lion’s share of recurrent climate-related crises occurring across many African countries, usually with devastating consequences for agriculture, households, food, and energy security [30]. The World Meteorological Organization (WMO) report highlights that 1.3 billion African people are at risk of climatic extremes, as the continent’s average temperature warms faster than the rest of the world, even though the continent only contributes 4% of global greenhouse emissions [31]. Sub-Saharan Africa is projected to be particularly impacted by climate change since it already experiences hot weather with decreased and variable precipitation. African dryland crop and livestock farming has adapted to domestic environmental conditions, but farming income will decrease if more hot or dry spells persist [32].
Studies have been carried out to monitor meteorological and agricultural droughts at the global scale [33,34,35] and also at the regional scale across Africa [36,37]. For instance, Ayugi et al. [36] reviewed meteorological drought events across Africa, along with historical trends, impacts, and mitigation practices. Crop yield is a prominent factor in agroecology that is highly affected by meteorological droughts in drought-prone areas of Africa. Abebe et al. [38] reported that more than 50% variability in crop yield was explained by the SPI on Ethiopian agricultural land. Similar studies assessing drought affecting multiple crop yields have been conducted in Central Malawi and Ethiopia using meteorological indices such as the SPI, SPEI, and PDSI [39,40]. Other than meteorological drought indices, composite drought indices such as the crop drought vulnerability index (CDVI) derived from the SPI are also being utilized to measure the vulnerability of crop yield towards droughts [41].
South Africa is one of the drought-prone African countries and has experienced several interconnected climatic extremes, including drought [42], floods [43], and heat stress [44]. This has resulted in infrastructure destruction, agricultural damage, and the loss of human lives [45,46]. It is also the second-largest wheat-producing region in Sub-Saharan Africa, with a mix of irrigated and dryland production. However, the low relative production of wheat in this region is attributed to droughts and heat stress, along with biotic diseases [47]. The recent drought in 2015–2017 in the Western Cape (i.e., the largest wheat-producing region of South Africa) strongly affected wheat production and reduced its exports [48,49]. Thus, drought is one of the most catastrophic weather events in South Africa, with a drastic impact on the agricultural and environmental sectors [50]. Furthermore, Mpandeli et al. [51] reported a rise in the drought intensity and frequency in South Africa due to increasing temperatures and a rainfall deficit (1960–2015). In addition to this, another study conducted by Lottering et al. [52] also indicated that local areas such as “uMsinga” in the “KwaZulu-Natal” province of South Africa experienced drought-threatened agricultural productivity on small-scale farms. Thus, climatic pressure induced by droughts is largely affecting the agricultural production system in this region due to inadequate irrigation water management [53]. Future climatic models of the region also reveal an approximate 20% decline in rainfall in the next 3 decades, which will ultimately be responsible for below-average crop yields [54].
Hence, keeping in view the background of the study and region, the rainfall-based meteorological assessment of drought and its interrelationship with wheat yield is much-needed research of this time. Although some studies have been conducted in this region on a smaller scale using multiple drought indices to monitor agricultural droughts in the region [55,56], to the best of our knowledge, no specific studies have been conducted to address long-term meteorological droughts and their impact on wheat yields on a broader scale in this region. Thus, the aim of this study is two-fold, with a focus on the regional scale covering 31 meteorological stations:
  • To examine and analyze the variability in short-term drought (SPI-3) occurrence and trends at all meteorological stations in the region over a time period of the past 20 years (2002–2021);
  • To explore the impacts of short-term meteorological droughts (SPI-3) on wheat yield loss and its resistance using a standardized yield residual series (SYRS) and the crop drought resistance factor (CR) in all provinces in the region.
This study provides fruitful findings on a historical baseline to address this issue for better drought hazard management in the future at the regional level.

2. Materials and Methods

2.1. Study Area and Data Collection

South Africa is situated between latitudes 22 and 35 °S and longitudes 17 and 33 °E and is neighbored by two oceans, the Indian Ocean in the east and the Atlantic Ocean in the west. The country spans a land area of 1,220,813 km2, partitions political boundaries with Mozambique, Zimbabwe, Botswana, Namibia, and Eswatini (Swaziland), and totally landlocks Lesotho [57] (Figure 1). It is described as a country with a semiarid climate that lies in the subtropics and the mid-latitudes [58].
South Africa is populated by 59.62 million people, and the major industrial activities are manufacturing, financial services, mining, tourism and trade, agriculture, and telecommunications [59]. On a regional scale, the country is divided into nine large provinces. Wheat is the main staple crop of the region after maize, with dryland and irrigated production. Climatic variations in the region make it vulnerable to extreme weather events such as drought and heat stress [47]. Rainfall in the region is quite variable from the eastern to western coastline based on the movement of oceanic currents [60]. The unpredicted rainfall variability in different seasons exposes wheat production to more climatic extremes [56].

2.2. Data Collection

For the drought analysis, the available historical climatic data were commissioned from the South African Weather Service (SAWS). Rainfall and temperature data spanning the 2000–2021 period was collected from 31 meteorological stations covering the whole region of South Africa (Table 1). Data quality and homogeneity were ensured by the South African Weather Service. The average temperature was 30.22 °C in January and 23.29 °C in June, while the rainfall ranged between 108.44 mm in January and 2.5 mm in August (Figure 2).
Wheat yield data at the provincial level (Table 1) were collated from the Crop Estimates Committee of the Department of Agriculture, Forestry and Fisheries for the period from 2000 to 2021. Wheat plays a crucial role in the agricultural economy of the country [47]. The wheat growth cycle in the region prevails from April to November, with slight variations in its summer and winter characteristics [61]. Irrigated wheat is planted in summer rainfall regions in eastern states from mid-May to the end of July [62], while in winter rain areas such as WES-C, it is planted from mid-April to mid-June [63]. The harvesting period of wheat in the region runs from October to November [61].
Between 500,000 and 900,000 hectares of cereal crop is cultivated per annum, with a mean yearly output of 1.3 to 2.4 million tons between 2000 and 2021. The area under irrigation produces approximately 5 tons per hectare on an annual basis, whereas the dryland produces 2–2.5 tons per hectare [64].

2.3. Standard Precipitation Index (SPI)

McKee et al. [11] prescribe the standard precipitation index (SPI) to capture the spatiotemporal variation in drought properties [65]. The SPI uses long-term monthly rainfall as input data and computes the divergence in rainfall from the average numerical parameter in a certain region during a specific time span [66]. The probability density function, such as the gamma statistical distribution function, is fitted to the rainfall data, which, according to Lloyd-Hughes and Saunders in [67], fits very well. The normalization of the gamma cumulative distribution function then follows [68]. The index can assimilate the drought duration, amount, and intensity on different time scales (3, 6, 9, 12, and 24 months). More details about the SPI calculation and classification can be found in McKee et al. [11]. Table 2 shows the categories of SPI index values [11]. According to Table 2, drought can be categorized according to specific SPI values, with extreme drought corresponding to less than −2 [69,70].

2.4. Drought Analysis

2.4.1. Drought Trend and Characteristics

The Mann–Kendall (MK) test [71] is applied to assess the monotonic tendencies of studied variables. This statistical test is nonparametric and assumes no normality but independent data [72]. In this case, the null (H0) assumes that there is no tendency in the studied variable, whereas the alternative hypothesis (Ha) presumes there is a tendency [73]. In addition, the Sen slope estimator was implemented to capture the values of changes during the study period [74].
Moreover, the drought characteristics computed in this study include the drought duration (DD) in number of months and the drought sum (DS), which defines the sum of all SPI values during a particular drought spell in months or years. The frequency of drought events is calculated by employing the equation [75,76]
DF = n s N s × 100
where n s represents the number of drought events during the selected time, and N s represents the total number of months, i.e., 240 in the currently selected period of 20 years (2002–2021).

2.4.2. Drought Impact on the Agricultural Sector

South African cereal production, particularly wheat (Figure 3), has been steadily increasing, mainly due to the expansion of cultivated land, an increase in irrigated areas, and the adoption of advanced agricultural inputs, such as improved seed cultivars [77].
Evolution of wheat yield on provincial scale over the selected time of 20 years is shown in Figure 3 where red line presents the yield trend and black dotted line presents the significance level.
To investigate the intercorrelation between agricultural drought (SPI-3) and wheat yield, the standardized yield residual series (SYRS) was calculated on a regional scale. A polynomial regression model was deployed to offset climatic, economic, and technological factors [78,79]. Yield variations due to nonclimatic factors (e.g., fertilizers, improved seed breeds, and irrigation operations) were separated using the detrending technique, and the remaining detrended yield was employed [80]. To signify the climatic effects on crop yield, the standardized yield residual series (SYRS) was computed employing the following equation:
SYRS = X ^ i μ σ
where X ^ i represents the residuals of the detrended yield; μ is the mean of the residuals of the detrended yield; and σ denotes the standard deviation (i: year). Table 3 shows the categories of the SYRS.
To showcase the impact of seasonal agrarian drought (SPI-3) on wheat production on a regional scale, the crop drought resilience factor (CR) was computed per province across South Africa. Mohammed et al. [81] define the CR as the ability of the cultivated crop to resist environmental pressures (such as drought) while maintaining its physiological, biochemical, and morphological duties. The CR [82] was calculated by the following equation:
CR = D d x D d t y
where D d x denotes values of yield in drought production seasons at the provincial level, and D d t represents values of detrended yield in a similar production season. Table 3 [82] depicts the CR categorization.

2.4.3. Correlation Analysis between Crop Yields and Agricultural Drought Indices

A nonlinear regression model was employed to explore the temporal relationship between drought indices and the SYRS. This relationship was studied to determine the impact of seasonal agrarian drought (SPI-3) on wheat crop yields on a regional scale. The Pearson correlation coefficient values between the SYRS and SPI-3 were calculated to determine the provinces that were severely affected by drought in South Africa (2002–2021).

3. Results

3.1. SPI Trend and Frequency on Regional and Provincial Scales

The output of the SPI analysis reveals that the northern and eastern parts of the country were less prone to drought compared with the southern and western parts (Figure 4 and Table 4). Both SPI-3 and SPI-6 showed that the majority of the stations experienced a negative MK trend (an increase in drought events) (Table 4) and shown by inverted red triangles in Figure 4. For SPI-3, only five stations (15.6%) experienced a positive MK trend (decrease in drought events), while eight stations (28.12%) witnessed a significant (p < 0.05) increase in drought events (Table 4) and shown in circular inverted red triangles in Figure 4. For SPI-6, only three stations (9.3%) were less vulnerable to drought (p < 0.05), while the rest (90%) witnessed an increase in the drought trend. In terms of SPI-9, fourteen stations (43%) were more susceptible to drought (p < 0.05), while four stations depicted a significant positive (p < 0.05) MK trend. Like SPI-9, the MK test for SPI-12 revealed that thirteen (40.6%) stations witnessed a significant (p < 0.05) increase in drought events, while only two stations were significantly (p < 0.05) less vulnerable to drought (Table 4).
Three stations, POR, IRE, and JHB BT, exhibited a significant positive (p < 0.05) trend (i.e., less affected by drought) (2000–2020) on the four SPI time scales (SPI-3, 6, 9, and 12) (Figure 4 and Table 4). Interestingly, these four stations were in the northern mountains and received an average monthly rainfall of 43.94–55 mm (Table 1 and Figure 1). In contrast, nine stations, namely, RSB (northern part), SKZ (eastern part), SPB (western part), PTA (southern part), PE (southern part), PTMBG (eastern part), CAL (southern part), CCE (southern part), and CPT (southern part), showed a fixed significant negative (p < 0.05) trend (Figure 1 and Figure 4; Table 4). Overall, drought events historically impacted the western part of the country and dominated in the coastal area.
On a regional scale, six provinces exhibited a negative SPI-3 trend (increase in drought) between 2002 and 2021. The highest decrease was recorded in WES-C (western part of the country, p < 0.05), followed by NR-C and NW (Figure 5). Two provinces experienced a positive SPI-3 value, namely, LPP and GG; however, this trend was not significant (Figure 5). Looking in depth at the SPI-3 values across all the provinces, most of the provinces experienced at least one event with less than −1.5, categorized as a severe drought event. However, the lowest value (−2.4) was recorded in WES-C (Figure 5).
Figure 6 presents the percentages of drought frequencies categorized from “no drought” to “extreme drought” based on SPI-3 ranges (Table 2) in all provinces of South Africa over a period of 20 years. A total of 1.3% of months had extreme drought events in the FS province, followed by 0.8% of months in the GG province and 0.4 % in NW. Similarly, the highest percentage of moderate to severe drought months was also experienced in the FS province, i.e., 11.3 and 3.4%, followed by 8.8% of months with moderate droughts in the LPP province of the region. The WES-C province of the region experienced 41.6% of months with mild droughts, 5.9% of months with moderate droughts, and only 0.8% of months with extreme drought events over a period of 20 years. Overall, the highest percentage of drought months was experienced by the FS province in the region, which makes it vulnerable to negative impacts. Overall, the highest percentage, i.e., 53.7%, of all droughts (mild to extreme) was experienced by ES-C, followed by 51.2% in NR-C, 50.8% in MP, and 50.4% in LPP. The lowest of all drought percentages, i.e., 46.2%, was experienced in the GG province of the region.

3.2. Impact of Drought on Wheat Production (SYRS)

The main idea of implementing the SYRS is to isolate the impact of agricultural development on crop yield due to climate conditions. The output of the SYRS could provide an overview of the direct impact of drought on wheat production in South Africa. In this research, SPI-3 was chosen as a representative of agricultural drought, and then the SYRS was calculated on a regional scale (i.e., the nine provinces).
In the NW province, the lowest SYRS value was recorded in 2006, where SYRSSPI-3 = −1.95, indicating the high impact of drought on wheat yield. In 2015, 2016, and 2019, the SYRSSPI-3 values were −1.21 (moderate losses), −0.89 (acceptable losses), and −1.00 (acceptable losses), respectively. However, the impact of drought (SPI-3) in the other years could be neglected (Table 5). For the GG province in 2003, 2006, and 2012, there was a negative impact on wheat yield, where the SYRSSPI-3 values were −1.89 (moderate losses), −1.59 (moderate losses), and −1.63 (moderate losses), respectively (Table 5). In the LPP province, the lowest SYRSSPI-3 value was recorded in 2012 (−1.82). The drought impacts were extreme in MP, where the SYRSSPI-3 value reached −2.37 in 2003. For both KZN and EsC provinces, drought had extreme impacts on wheat production, where SYRSSPI-3 = −1.79 and SYRSSPI-3 = −1.72, respectively. The highest value of SYRSSPI-3 was recorded in 2009, revealing extreme wheat loss in the FS province. In 2012, the SYRSSPI-3 reached the lowest value (−1.78) in the NR-C province (Table 5). However, the WES-C province was affected by drought in 2017 and 2019, where the SYRSSPI-3 analysis showed extreme drought impacts (SYRSSPI-3 = −2 (2017), SYRSSPI-3 = −1.74 (2019)).
By tracking the impact of drought (SPI-3) on wheat yield across the country, three years can be distinguished based on the values of SYRSSPI-3 (≤−1.5). In this sense, 2003 had a negative impact on wheat production in MP, LPP, and GG provinces. Three provinces, ES-C, KZN, and NW, were affected in 2005, while a recent drought in 2019 had a direct impact on FS, WES-C, and NR-C (Table 5).

3.3. Correlation between SYRS and SPI-3 on a Monthly Time Scale

The impact of drought on crop yield varied between the provinces. Table 6 depicts the correlation between SYRS and SPI-3 on a monthly scale in each province across South Africa. For WES-C, the SYRS had a positive correlation for all months. However, the highest correlations were in March (r SPI-3 vs. SYRS = 0.53), July (r SPI-3 vs. SYRS = 0.55), and August (r SPI-3 vs. SYRS = 0.55) (Table 6). In NR-C, the highest correlation was recorded in the growing cycle (Table 6). Like NR-C, the highest correlation in the FS province was recorded between October (r SPI-3 vs. SYRS = 0.54) and December (r SPI-3 vs. SYRS = 0.54) (Table 6). For both NS-C and NKZ, there was a low correlation between SPI-3 and SYRS (Table 6). For MP, the highest correlation was recorded in May (r SPI-3 vs. SYRS = 0.41) (Table 6). The correlation was weak in LPP from January to May and had the lowest value (r SPI-3 vs. SYRS = 0.2, April–May) (Table 6). In the GG province, the growing cycle reflects a good correlation with SPI-3 in May (r SPI-3 vs. SYRS = 0.35), June (r SPI-3 vs. SYRS = 0.51), and July (r SPI-3 vs. SYRS = 0.53) (Table 6). The growing cycle in the NW province did not reveal any notable correlation between SYRS values and drought (SPI-3) (Table 6).

3.4. Drought Resilience (CR) of Wheat on a Regional Scale

The CR analysis provides an overview of crop resistance to drought. In this research, the CR was calculated for the nine provinces across South Africa. All results were compared with the threshold CR = 0.8 (Table 3) to indicate whether wheat yield was moderately resilient to drought events or not. The most affected province by drought was WES-C (western part of the country), where the CR value was 0.65, followed by FS with CR = 0.65 (severely nonresilient) (Table 7). Interestingly, the SYRS analysis also revealed that WES-C had a positive correlation with SPI-3 in all months, while FS also had a positive correlation with SPI-3 in the harvesting months of Oct–Nov (Table 6), thus designating these provinces as the most affected by agricultural drought. Regarding the rest of the provinces, the CR values were above the threshold, i.e., ranging from 0.85 to 0.97, revealing a good deal of drought resilience.
Table 7 presents an extensive analysis of the yield loss percentage (YL %) in different stages of wheat growth and during different drought events. Consistent with the correlation evaluation, WES-C revealed the highest YL of 35% during the longest SPI-3 drought duration (DD) of 20 months, i.e., from August 2016 to March 2018, in the whole growing cycle (GC) of wheat, with a drought sum (DS) of 13.3. This was followed by a yield loss of 29.5% during another longer DD of 18 months from August 2018 to January 2020 in the GC of wheat with a DS of 14. WES-C also experienced a YL of 29.3% in 3 months of the sowing period (SP) in 2004 and 18.6% in 6 months from the SP to the growing period (GP) and 4 months from the GP to the harvesting period (HP) in 2003 and 2010, respectively.
Like WES-C, another province in the region i.e., FS, also experienced the highest SPI-3 drought-associated YL of 34.5% in 7 months from the growing period to the harvesting period (GP-HP) from July 2019 to Jan 2020. Other significant YLs of 20.4% and 17% in the same (GP-HP) stages of wheat growth were observed in 2010 and 2004, with short DDs of 4 and 5 months. The results also revealed GP to be most closely linked to YL, followed by HP, GC, and SP collectively in both provinces in the region (Table 7).
Similar results are found for other provinces in the region during different growth stages but with a good drought resilience of above 0.8. Another interesting finding is that a minimum of 7% to a maximum of 17.8% YL was found in NW, KZN, NR-C, and ES-C provinces of the region during different wet years when there was no drought identified (Table 7).

4. Discussion

4.1. Current and Future Drought across South Africa

In South Africa, a few studies have examined the drought trend and its intercorrelation with crop production. In this sense, this research was designed to bridge the gap regarding agricultural drought and its impact in South Africa and to highlight the region’s most vulnerable provinces to drought. Thus, local governments and policymakers can take action to minimize drought impacts and ensure food security through adaptation and mitigation plans. The output of the SPI analysis for 32 stations covering the whole of South Africa indicates a drought tendency across the country at different drought levels, e.g., SPI-3, SPI-6, SPI-9, and SPI-12 (Figure 4 and Table 4). Notably, most of the provinces showed a negative trend based on the SPI-3 analysis (Figure 4 and Figure 5). This drought trend could be explained by the El Niño–Southern Oscillation. In this sense, the recent 2016 El Niño indicated that the country is prone to drought trends [83,84]. The Pacific El Niño has a history of causing meteorological variability [85,86]. The El Niño–Southern Oscillation is reportedly causing drought periods in the northern region of the country [87]. In this context, the arid and semiarid climates of the study area, along with fluctuating rainfall, has accelerated the evolution of drought in the country [88]. However, a drought trend has been reported previously across Africa and in South Africa (Table 8).
Future climate projection by GCM models (from the 1960–2000 baseline timeframe) highlighted in the IPCC 4th Assessment Report [100] shows that South African rainfall is predicted to decrease by 2030–2060. The report mentions that the temperature is expected to increase by 1 and 3 degrees Celsius in most inland areas by 2060, while coastal areas will experience lower increases compared with the interior. This indicates that the surface wind direction and speed are predicted to change with high pressures from anticyclones in both the Pacific and Indian Oceans.

4.2. Drought Impacts on Wheat and Its Resilience

Since most droughts in Africa occur in the temporal and permanent domains, crop prediction models that are not based on these variable rates may give false positives. Drought has severe agricultural impacts, from the loss of income for farmers due to crop yield and livestock losses [101] to regional food security shortages [102]. Based on the SPI-3 versus SYRS analysis (Table 5 and Table 6), drought has had an impact on wheat production, especially in the western parts of the country. Pearson correlations between SYRS and SPI-3 on the monthly scale clearly show the strength of positive and negative relationships between yield loss and drought during stages of wheat growth. A positive correlation between them is found in WES-C during all stages of wheat growth, i.e., from Apr to Nov, over a period of 20 years. Severe droughts in 2003, also studied by Rouault and Richard [103], followed by low to moderate droughts in 2009 and 2010, became a cause of drought-associated yield loss here. Furthermore, after 2015, WES-C faced a continuous long-term drought condition that significantly impacted the whole wheat growing cycle and caused the highest YL of 29–35% in the whole of South Africa (Table 7) [56]. Similarly, LPP, GG, and NW provinces are also impacted by drought during different stages of wheat growth, with positive correlations between SYRS and SPI-3 from June to November, April to July, and January to April, respectively. In contrast, three provinces in the region, i.e., FS, ES-C, and KZN, revealed negative correlations between the SYRS and SPI-3 from Mar to Aug and April to December, respectively. A study by Shew et al. [47] revealed wheat yield loss in the dryland cropping system of FS due to heat stress, which is consistent with our results of significant SPI-3-associated YLs of 10.8 to 34.5% from the sowing to harvesting period (SP-HP) of wheat growth (Table 7).
Drought impacts crop risk profiling and seasonal crop-water requirements, which are fundamental to crop life. [104]. Unfortunately, only 25% of the total cropped area in the region is under irrigation [105]; therefore, the SPI-associated yield loss is quite evident in the results of our study.
A few studies have been carried out in South Africa to highlight the impact of drought on crop production. For instance, Unganai et al. [106] focused on the interaction between drought monitoring and corn yield predictions using remote sensing, with the findings highlighting the need to further investigate the severity of vegetation stress. Since most African countries do not possess long-term instrumental climate records, drought prediction using statistical methods is not prevalent. Senay and Verdin [107] allude to using modern GIS water balance algorithms to forecast seasonal drought to offset the negative impact of drought in agropastoral systems and ensure water allocation.

4.3. Strategies for Drought Mitigation in South Africa and Future Steps

The way to combat agricultural drought is to integrate robust sectoral strategies [108]. Adaptation strategies surely depend on the agropastoral practices that farmers are using. For example, cereal farmers may need to grow fall-sown crops and use better cultivars to offset the impact of drought. Other strategies for cereals sown in the winter include cultivars with reduced vernalization periods. Land management practices such as shifting from rainfed agriculture to irrigated agriculture can alleviate water stress in crops during drought seasons [109]. Reducing greenhouse gases by adopting precision agriculture practices will help alleviate carbon escape from most croplands [110]. Crop diversity, stubble residue management, and nutrient recycling procedures may help reduce the impact of drought [111].
Given that agricultural drought results in a persistent deficit in soil moisture content, which is associated with wilting crops [112], it is urgent for farmers and policymakers to adopt measures such as irrigation supply, water demand, and aftermath mitigation measures. The first two measures deal with water shortages, while the last addresses the socio-environmental impact of drought. Earthwork can improve water availability through the abstraction of groundwater and the construction of artificial dams, reservoirs, canals, or rivers. Non-earthwork systems can sound early warning alarms for drought detection and thus enable emergency responses, such as insurance aid, rehabilitation, and recovery measures. These types of measures strengthen institutional structures and their capacity to better prepare for drought [113].
Even though this research was based on input from 32 stations, the distribution of the climate stations covers the whole region and, thus, could be used to reach a robust result. Forthcoming studies will include engaging other climate indices, such as the SPEI and the Palmer drought severity index (PDSI), to identify the drought frequency and intensity based on different inputs. Furthermore, the ecosystem response to drought will be analyzed using satellite images.

4.4. Limitations of Study

Our study has some limitations that open new research questions for the future as well. Firstly, the study was limited to analyzing the meteorological drought impact on wheat yield using a single drought index, i.e., the SPI, which is comparable to other indices, such as the SPEI, PDSI, and scPDSI, for incorporating the impacts of temperature and evapotranspiration. A similar type of study for drought assessment in east Africa utilized the SPI and SPEI and presented quite similar results, with some uncertainties in the over- and underestimation of drought events [98]. However, the results of multiple indices may vary from region to region and can be utilized in the future to differentiate drought impacts in rainfed and irrigated wheat regions. In the next step, multiple indices, such as the SPEI, PDSI, and composite drought index (CDI) [114], will be used to capture drought events in South Africa and to assess their impacts on the agricultural sector.

5. Conclusions

Our study examined the temporal interaction of meteorological droughts and their impacts on wheat yield across the nine provinces of South Africa. A widely utilized meteorological drought index, i.e., the SPI, was used to examine the duration, frequency, and trend of meteorological droughts in dryland and irrigated provinces in the region. The SYRS and CR computed from trended and detrended wheat yields across all provinces were revealed to be significantly impacted by meteorological drought over a period of 20 years. The following major conclusions can be drawn from the results of this study:
  • The frequency of drought events revealed that ES-C experienced the highest percentage of drought, i.e., 53.7%, followed by NR-C, MP, and LPP provinces of the region.
  • SPI-3 trend analysis reveals a significant negative trend across many provinces in the region. Specifically, the western coastal provinces WES-C and NR-C have been more vulnerable to meteorological droughts over the past 20 years.
  • Wheat yield loss analysis reveals that the highest SYRS = −2.52 was found in FS in 2019, followed by −2.37 in MP in 2003, −1.95 in NW in 2006, and −1.899 in GG and ES-C in 2003 and 2014, respectively. The most dreadful drought impacts on wheat yield were observed in the years 2015–2016, when all provinces experienced significant yield losses.
  • Positive correlation results between the SYRS and wheat yield indicate that the WES-C province was highly influenced by drought during all stages of wheat growth, i.e., Apr–Nov. Historical drought spells in 2003, 2009, and 2010 and a low CR = 0.64 caused the province to be highly impacted by the negative impacts of droughts on yield loss.
  • Some provinces in the region, including FS, ES-C, and KZN, were not found to be highly impacted by droughts, with negative correlations between the SYRS and SPI-3 during the wheat growth cycle from Apr to Nov.
  • The WES-C and FS provinces of the region experienced the highest yield loss % during the SP-GP-HP of wheat growth stages with a CR of 0.65, indicating extremely low resilience. Overall, the growing period of wheat was found to be the most associated with yield loss, followed by the harvesting and sowing periods. Yield loss in the WES-C province is linked to the whole growing cycle in all months of wheat growth (Apr–Nov).
However, the current study efficiently examined the meteorological drought variations and related wheat yield losses in both rainfed and irrigated dryland provinces of South Africa, but still, the results of the study can be further enhanced by incorporating other indices, such as the SPEI and PDSI. Other than these, the incorporation of remote-sensing-based indices such as NDVI can also be utilized to examine drought-associated stress and yield loss. This research can be helpful in robust yield prediction modeling associated with climatic changes in the region.
Other than this, the major findings of the study also suggest adapting and focusing on drought-resistant agricultural practices in the western coastal parts of South Africa to prevent future yield loss risk. Our study recommends an immediate climate adaptation and mitigation plan to support farmers and stakeholders in combating climate change. A regional plan for climate awareness should also be formulated to protect the country’s food security.

Author Contributions

Conceptualization: S.M. and G.N.; writing—original draft: G.N and S.M.; methodology: S.M. and K.A.; writing—review and editing: B.B., H.B., K.A., E.H. and S.A. All authors have read and agreed to the published version of the manuscript.

Funding

Project no. TKP2021-NKTA-32 was implemented with support provided by the National Research, Development, and Innovation Fund of Hungary, financed under the TKP2021-NKTA funding scheme. This research was supported by the RESEARCHERS SUPPORTING PROJECT, grant number RSP-2021/296, King Saud University, Riyadh, Saudi Arabia.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank Debrecen University for their unlimited support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Houghton, J.T.; Ding, Y.; Griggs, D.J.; Noguer, M.; van der Linden, P.J.; Dai, X.; Maskell, K.; Johnson, C. Climate Change 2001: The Scientific Basis: Contribution of Working Group I to the Third Assessment Report of the Intergovernmental Panel on Climate Change; Cambridge University Press: Cambridge, UK, 2001. [Google Scholar]
  2. Pörtner, H.-O.; Roberts, D.C.; Adams, H.; Adler, C.; Aldunce, P.; Ali, E.; Begum, R.A.; Betts, R.; Kerr, R.B.; Biesbroek, R. Climate change 2022: Impacts, adaptation and vulnerability. In IPCC Sixth Assessment Report; IPCC: Geneva, Switzerland, 2022. [Google Scholar]
  3. Zhongming, Z.; Linong, L.; Xiaona, Y.; Wangqiang, Z.; Wei, L. AR6 Climate Change 2021: The Physical Science Basis; IPCC: Geneva, Switzerland, 2021. [Google Scholar]
  4. Prentice, I.C.; Farquhar, G.; Fasham, M.; Goulden, M.L.; Heimann, M.; Jaramillo, V.; Kheshgi, H.; Le Quéré, C.; Scholes, R.J.; Wallace, D.W. The carbon cycle and atmospheric carbon dioxide. In Climate Change 2001: The Scientific Basis, Intergovernmental Panel on Climate Change; IPCC: Geneva, Switzerland, 2001. [Google Scholar]
  5. Gillett, N.P.; Graf, H.F.; Osborn, T.J. Climate change and the North Atlantic oscillation. Geophys. Monogr.-Am. Geophys. Union 2003, 134, 193–210. [Google Scholar]
  6. Kogan, F.N. Global drought watch from space. Bull. Am. Meteorol. Soc. 1997, 78, 621–636. [Google Scholar] [CrossRef]
  7. Kogan, F.; Guo, W.; Yang, W. Drought and food security prediction from NOAA new generation of operational satellites. Geomat. Nat. Hazards Risk 2019, 10, 651–666. [Google Scholar] [CrossRef] [Green Version]
  8. Alsafadi, K.; Mohammed, S.; Ayugi, B.; Sharaf, M.; Harsányi, E. Spatial–temporal evolution of drought characteristics over Hungary between 1961 and 2010. Pure Appl. Geophys. 2020, 177, 3961–3978. [Google Scholar] [CrossRef] [Green Version]
  9. Mekonen, A.A.; Berlie, A.B.; Ferede, M.B. Spatial and temporal drought incidence analysis in the northeastern highlands of Ethiopia. Geoenviron. Disasters 2020, 7, 10. [Google Scholar] [CrossRef] [Green Version]
  10. Zhao, L.; Lyu, A.; Wu, J.; Hayes, M.; Tang, Z.; He, B.; Liu, J.; Liu, M. Impact of meteorological drought on streamflow drought in Jinghe River Basin of China. Chin. Geogr. Sci. 2014, 24, 694–705. [Google Scholar] [CrossRef]
  11. Mckee, T.; Doesken, N.; Kleist, J. The relationship of drought frequency and duration to time scales. Appl. Climatol. 1993, 17, 179–183. [Google Scholar]
  12. Tsakiris, G.; Pangalou, D.; Vangelis, H. Regional Drought Assessment Based on the Reconnaissance Drought Index (RDI). Water Resour. Manag. 2007, 21, 821–833. [Google Scholar] [CrossRef]
  13. Vicente-Serrano, S.M.; Beguería, S.; López-Moreno, J.I. A Multiscalar Drought Index Sensitive to Global Warming: The Standardized Precipitation Evapotranspiration Index. J. Clim. 2010, 23, 1696–1718. [Google Scholar] [CrossRef] [Green Version]
  14. Alley, W.M. The palmer drought severity index as a measure of hydrologic drought. J. Am. Water Resour. Assoc. 1985, 21, 105–114. [Google Scholar] [CrossRef]
  15. Palmer, W.C. Keeping Track of Crop Moisture Conditions, Nationwide: The New Crop Moisture Index. Weatherwise 1968, 21, 156–161. [Google Scholar] [CrossRef]
  16. Vicente-Serrano, S.M.; Beguería, S.; Lorenzo-Lacruz, J.; Camarero, J.J.; López-Moreno, J.I.; Azorin-Molina, C.; Revuelto, J.; Morán-Tejeda, E.; Sanchez-Lorenzo, A. Performance of Drought Indices for Ecological, Agricultural, and Hydrological Applications. Earth Interact. 2012, 16, 1–27. [Google Scholar] [CrossRef]
  17. Du, J.; Fang, J.; Xu, W.; Shi, P. Analysis of dry/wet conditions using the standardized precipitation index and its potential usefulness for drought/flood monitoring in Hunan Province, China. Stoch. Environ. Res. Risk Assess. 2013, 27, 377–387. [Google Scholar] [CrossRef]
  18. Xu, X.; Gao, P.; Zhu, X.; Guo, W.; Ding, J.; Li, C. Estimating the responses of winter wheat yields to moisture variations in the past 35 years in Jiangsu Province of China. PLoS ONE 2018, 13, e0191217. [Google Scholar] [CrossRef] [PubMed]
  19. Mohammed, S.; Alsafadi, K.; Al-Awadhi, T.; Sherief, Y.; Harsanyie, E.; El Kenawy, A.M. Space and time variability of meteorological drought in Syria. Acta Geophys. 2020, 68, 1877–1898. [Google Scholar] [CrossRef]
  20. Asadi Zarch, M.A.; Sivakumar, B.; Sharma, A. Droughts in a warming climate: A global assessment of Standardized precipitation index (SPI) and Reconnaissance drought index (RDI). J. Hydrol. 2015, 526, 183–195. [Google Scholar] [CrossRef]
  21. Tirivarombo, S.; Osupile, D.; Eliasson, P. Drought monitoring and analysis: Standardised Precipitation Evapotranspiration Index (SPEI) and Standardised Precipitation Index (SPI). Phys. Chem. Earth Parts A/B/C 2018, 106, 1–10. [Google Scholar] [CrossRef]
  22. Blain, G.C. Revisiting the probabilistic definition of drought: Strengths, limitations and an agrometeorological adaptation. Bragantia 2012, 71, 132–141. [Google Scholar] [CrossRef] [Green Version]
  23. Okpara, J.N.; Afiesimama, E.A.; Anuforom, A.C.; Owino, A.; Ogunjobi, K.O. The applicability of Standardized Precipitation Index: Drought characterization for early warning system and weather index insurance in West Africa. Nat. Hazards 2017, 89, 555–583. [Google Scholar] [CrossRef]
  24. Guttman, N.B. Accepting the standardized precipitation index: A calculation algorithm 1. J. Am. Water Resour. Assoc. 1999, 35, 311–322. [Google Scholar] [CrossRef]
  25. Yerdelen, C.; Abdelkader, M.; Eris, E. Assessment of drought in SPI series using continuous wavelet analysis for Gediz Basin, Turkey. Atmos. Res. 2021, 260, 105687. [Google Scholar] [CrossRef]
  26. Jordaan, A.J.; Mlenga, D.H.; Mandebvu, B. Monitoring droughts in Eswatini: A spatiotemporal variability analysis using the Standard Precipitation Index. Jàmbá J. Disaster Risk Stud. 2019, 11, 1–11. [Google Scholar]
  27. Guenang, G.M.; Kamga, F.M. Computation of the standardized precipitation index (SPI) and its use to assess drought occurrences in Cameroon over recent decades. J. Appl. Meteorol. Climatol. 2014, 53, 2310–2324. [Google Scholar] [CrossRef]
  28. Harsányi, E.; Bashir, B.; Alsilibe, F.; Alsafadi, K.; Alsalman, A.; Széles, A.; Rahman, M.H.u.; Bácskai, I.; Juhász, C.; Ratonyi, T. Impact of agricultural drought on sunflower production across Hungary. Atmosphere 2021, 12, 1339. [Google Scholar] [CrossRef]
  29. Jiang, R.; Xie, J.; He, H.; Luo, J.; Zhu, J. Use of four drought indices for evaluating drought characteristics under climate change in Shaanxi, China: 1951–2012. Nat. Hazards 2015, 75, 2885–2903. [Google Scholar] [CrossRef]
  30. Araneda-Cabrera, R.J.; Bermúdez, M.; Puertas, J. Assessment of the performance of drought indices for explaining crop yield variability at the national scale: Methodological framework and application to Mozambique. Agric. Water Manag. 2021, 246, 106692. [Google Scholar] [CrossRef]
  31. WMO. State of the Climate in Africa. Available online: https://public.wmo.int/en/our-mandate/climate/wmo-statement-state-of-global-climate/Africa (accessed on 25 March 2022).
  32. Kurukulasuriya, P.; Mendelsohn, R.; Hassan, R.; Benhin, J.; Deressa, T.; Diop, M.; Eid, H.M.; Fosu, K.Y.; Gbetibouo, G.; Jain, S. Will African agriculture survive climate change? World Bank Econ. Rev. 2006, 20, 367–388. [Google Scholar] [CrossRef] [Green Version]
  33. Afshar, M.H.; Bulut, B.; Duzenli, E.; Amjad, M.; Yilmaz, M.T. Global spatiotemporal consistency between meteorological and soil moisture drought indices. Agric. For. Meteorol. 2022, 316, 108848. [Google Scholar] [CrossRef]
  34. Turco, M.; Jerez, S.; Donat, M.G.; Toreti, A.; Vicente-Serrano, S.M.; Doblas-Reyes, F.J. A Global Probabilistic Dataset for Monitoring Meteorological Droughts. Bull. Am. Meteorol. Soc. 2020, 101, E1628–E1644. [Google Scholar] [CrossRef] [Green Version]
  35. Sánchez, N.; González-Zamora, Á.; Martínez-Fernández, J.; Piles, M.; Pablos, M. Integrated remote sensing approach to global agricultural drought monitoring. Agric. For. Meteorol. 2018, 259, 141–153. [Google Scholar] [CrossRef]
  36. Ayugi, B.; Eresanya, E.O.; Onyango, A.O.; Ogou, F.K.; Okoro, E.C.; Okoye, C.O.; Anoruo, C.M.; Dike, V.N.; Ashiru, O.R.; Daramola, M.T. Review of meteorological drought in Africa: Historical trends, impacts, mitigation measures, and prospects. Pure Appl. Geophys. 2022, 179, 1365–1386. [Google Scholar] [CrossRef] [PubMed]
  37. Ngcamu, B.S.; Chari, F. Drought influences on food insecurity in Africa: A Systematic literature review. Int. J. Environ. Res. Public Health 2020, 17, 5897. [Google Scholar] [CrossRef] [PubMed]
  38. Abebe, S.; Solomon, A.; Suryabhagavan, K.V. Mapping the spatial and temporal variation of agricultural and meteorological drought using geospatial techniques, Ethiopia. Environ. Syst. Res. 2021, 10, 15. [Google Scholar] [CrossRef]
  39. Chikabvumbwa, S.R.; Salehnia, N.; Manzanas, R.; Abdelbaki, C.; Zerga, A. Assessing the effect of spatial–temporal droughts on dominant crop yield changes in Central Malawi. Environ. Monit. Assess. 2022, 194, 63. [Google Scholar] [CrossRef]
  40. Temam, D.; Uddameri, V.; Mohammadi, G.; Hernandez, E.A.; Ekwaro-Osire, S. Long-Term Drought Trends in Ethiopia with Implications for Dryland Agriculture. Water 2019, 11, 2571. [Google Scholar] [CrossRef] [Green Version]
  41. Kamali, B.; Abbaspour, K.C.; Wehrli, B.; Yang, H. Drought vulnerability assessment of maize in Sub-Saharan Africa: Insights from physical and social perspectives. Glob. Planet. Chang. 2018, 162, 266–274. [Google Scholar] [CrossRef] [Green Version]
  42. Ziervogel, G.; New, M.; Archer van Garderen, E.; Midgley, G.; Taylor, A.; Hamann, R.; Stuart-Hill, S.; Myers, J.; Warburton, M. Climate change impacts and adaptation in South Africa. Wiley Interdiscip. Rev. Clim. Chang. 2014, 5, 605–620. [Google Scholar] [CrossRef]
  43. Olorunfemi, F. Managing flood disasters under a changing climate: Lessons from Nigeria and South Africa. In NISER Research Seminar Series; NISER: Ibadan, Nigeria, 2011; pp. 1–44. [Google Scholar]
  44. Nesamvuni, E.; Lekalakala, R.; Norris, D.; Ngambi, J. Effects of climate change on dairy cattle, South Africa. Afr. J. Agric. Res. 2012, 7, 3867–3872. [Google Scholar] [CrossRef]
  45. Calzadilla, A.; Zhu, T.; Rehdanz, K.; Tol, R.S.; Ringler, C. Climate change and agriculture: Impacts and adaptation options in South Africa. Water Resour. Econ. 2014, 5, 24–48. [Google Scholar] [CrossRef]
  46. Wright, C.Y.; Kapwata, T.; Du Preez, D.J.; Wernecke, B.; Garland, R.M.; Nkosi, V.; Landman, W.A.; Dyson, L.; Norval, M. Major climate change-induced risks to human health in South Africa. Environ. Res. 2021, 196, 110973. [Google Scholar] [CrossRef]
  47. Shew, A.M.; Tack, J.B.; Nalley, L.L.; Chaminuka, P. Yield reduction under climate warming varies among wheat cultivars in South Africa. Nat. Commun. 2020, 11, 4408. [Google Scholar] [CrossRef] [PubMed]
  48. Otto, F.E.; Wolski, P.; Lehner, F.; Tebaldi, C.; Van Oldenborgh, G.J.; Hogesteeger, S.; Singh, R.; Holden, P.; Fučkar, N.S.; Odoulami, R.C. Anthropogenic influence on the drivers of the Western Cape drought 2015–2017. Environ. Res. Lett. 2018, 13, 124010. [Google Scholar] [CrossRef]
  49. Cullis, J.; Alton, T.; Arndt, C.; Cartwright, A.; Chang, A.; Gabriel, S.; Gebretsadik, Y.; Hartley, F.; de Jager, G.; Makrelov, K. An Uncertainty Approach to Modelling Climate Change Risk in South Africa; 929230934X; WIDER Working Paper: Helsinki, Finland, 2015. [Google Scholar]
  50. Meza, I.; Rezaei, E.E.; Siebert, S.; Ghazaryan, G.; Nouri, H.; Dubovyk, O.; Gerdener, H.; Herbert, C.; Kusche, J.; Popat, E. Drought risk for agricultural systems in South Africa: Drivers, spatial patterns, and implications for drought risk management. Sci. Total Environ. 2021, 799, 149505. [Google Scholar] [CrossRef] [PubMed]
  51. Mpandeli, S.; Nhamo, L.; Moeletsi, M.; Masupha, T.; Magidi, J.; Tshikolomo, K.; Liphadzi, S.; Naidoo, D.; Mabhaudhi, T. Assessing climate change and adaptive capacity at local scale using observed and remotely sensed data. Weather. Clim. Extrem. 2019, 26, 100240. [Google Scholar] [CrossRef]
  52. Lottering, S.J.; Mafongoya, P.; Lottering, R. The impacts of drought and the adaptive strategies of small-scale farmers in uMsinga, KwaZulu-Natal, South Africa. J. Asian Afr. Stud. 2021, 56, 267–289. [Google Scholar] [CrossRef]
  53. Dube, E.; Tsilo, T.J.; Sosibo, N.Z.; Fanadzo, M. Irrigation wheat production constraints and opportunities in South Africa. South Afr. J. Sci. 2020, 116, 1–6. [Google Scholar] [CrossRef]
  54. Conway, D.; van Garderen, E.A.; Deryng, D.; Dorling, S.; Krueger, T.; Landman, W.; Lankford, B.; Lebek, K.; Osborn, T.; Ringler, C.; et al. Climate and southern Africa’s water–energy–food nexus. Nat. Clim. Chang. 2015, 5, 837–846. [Google Scholar] [CrossRef] [Green Version]
  55. Kganvago, M.; Mukhawana, M.B.; Mashalane, M.; Mgabisa, A.; Moloele, S. Recent Trends of Drought Using Remotely Sensed and In-Situ Indices: Towards an Integrated Drought Monitoring System for South Africa. In Proceedings of the 2021 IEEE International Geoscience and Remote Sensing Symposium IGARSS, Brussels, Belgium, 11–16 July 2021; pp. 6225–6228. [Google Scholar]
  56. Theron, S.; Archer, E.; Midgley, S.; Walker, S. Agricultural perspectives on the 2015–2018 western cape drought, South Africa: Characteristics and spatial variability in the core wheat growing regions. Agric. For. Meteorol. 2021, 304, 108405. [Google Scholar] [CrossRef]
  57. Knight, J.; Rogerson, C.M. The Geography of South Africa; Springer: Berlin/Heidelberg, Germany, 2019. [Google Scholar]
  58. Engelbrecht, C.J.; Landman, W.A.; Engelbrecht, F.A.; Malherbe, J. A synoptic decomposition of rainfall over the Cape south coast of South Africa. Clim. Dyn. 2015, 44, 2589–2607. [Google Scholar] [CrossRef]
  59. SAG. South Africa at a Glance. Available online: https://www.gov.za/about-sa/south-africa-glance#:~:text=Mining%2C%20transport%2C%20energy%2C%20manufacturing%2C%20tourism%20and%20agriculture (accessed on 25 April 2022).
  60. McBride, C.M.; Kruger, A.C.; Dyson, L. Changes in extreme daily rainfall characteristics in South Africa: 1921–2020. Weather. Clim. Extrem. 2022, 38, 100517. [Google Scholar] [CrossRef]
  61. Franch, B.; Cintas, J.; Becker-Reshef, I.; Sanchez-Torres, M.J.; Roger, J.; Skakun, S.; Sobrino, J.A.; Van Tricht, K.; Degerickx, J.; Gilliams, S.; et al. Global crop calendars of maize and wheat in the framework of the WorldCereal project. GISci. Remote Sens. 2022, 59, 885–913. [Google Scholar] [CrossRef]
  62. Tadesse, W.; Bishaw, Z.; Assefa, S. Wheat production and breeding in Sub-Saharan Africa: Challenges and opportunities in the face of climate change. Int. J. Clim. Chang. Strateg. Manag. 2018, 11, 696–715. [Google Scholar] [CrossRef] [Green Version]
  63. RSA. Production Guideline for Wheat; agriculture, forestry and fisheries; Republic of South Africa: Bloemfontein, South Africa, 2016.
  64. Crop Estimate Committee. Available online: https://www.sagis.org.za/cec_reports.html (accessed on 23 April 2022).
  65. Laddimath, R.S.; Patil, N.S.; Rao, P.; Nagendra, N. Assessing the impacts of climate change on drought-prone regions in Bhima sub-basin (India) using the Standard Precipitation Index. J. Water Clim. Chang. 2022, 13, 817–838. [Google Scholar] [CrossRef]
  66. Zuo, D.-D.; Hou, W.; Zhang, Q.; Yan, P.-C. Sensitivity analysis of standardized precipitation index to climate state selection in China. Adv. Clim. Chang. Res. 2022, 13, 42–50. [Google Scholar] [CrossRef]
  67. Lloyd-Hughes, B.; Saunders, M.A. A drought climatology for Europe. Int. J. Climatol. A J. R. Meteorol. Soc. 2002, 22, 1571–1592. [Google Scholar] [CrossRef]
  68. Mathier, L.; Perreault, L.; Bobée, B.; Ashkar, F. The use of geometric and gamma-related distributions for frequency analysis of water deficit. Stoch. Hydrol. Hydraul. 1992, 6, 239–254. [Google Scholar] [CrossRef]
  69. Svoboda, M.; Hayes, M.; Wood, D. Standardized Precipitation Index: User Guide; World Meteorological Organization: Geneva, Switzerland, 2012. [Google Scholar]
  70. Kim, D.-W.; Byun, H.-R.; Choi, K.-S. Evaluation, modification, and application of the Effective Drought Index to 200-Year drought climatology of Seoul, Korea. J. Hydrol. 2009, 378, 1–12. [Google Scholar] [CrossRef]
  71. Mann, H.B. Nonparametric tests against trend. Econom. J. Econom. Soc. 1945, 13, 245–259. [Google Scholar] [CrossRef]
  72. Wasserstein, R.L.; Schirm, A.L.; Lazar, N.A. Moving to a world beyond “p < 0.05”. Am. Stat. 2019, 73, 1–19. [Google Scholar]
  73. Koudahe, K.; Koffi, D.; Kayode, J.; Awokola, S.; Adebola, A. Impact of climate variability on crop yields in southern Togo. Environ. Pollut. Clim. Chang. 2018, 2, 1–9. [Google Scholar]
  74. Li, Z.; Ali, Z.; Cui, T.; Qamar, S.; Ismail, M.; Nazeer, A.; Faisal, M. A comparative analysis of pre- and post-industrial spatiotemporal drought trends and patterns of Tibet Plateau using Sen slope estimator and steady-state probabilities of Markov Chain. Nat. Hazards 2022, 113, 547–576. [Google Scholar] [CrossRef]
  75. Zambreski, Z.T.; Lin, X.; Aiken, R.M.; Kluitenberg, G.J.; Pielke Sr, R.A. Identification of hydroclimate subregions for seasonal drought monitoring in the U.S. Great Plains. J. Hydrol. 2018, 567, 370–381. [Google Scholar] [CrossRef]
  76. Gebremichael, H.B.; Raba, G.A.; Beketie, K.T.; Feyisa, G.L. Temporal and spatial characteristics of drought, future changes and possible drivers over Upper Awash Basin, Ethiopia, using SPI and SPEI. Environ. Dev. Sustain. 2022, 1–39. [Google Scholar] [CrossRef]
  77. Nhamo, L.; Matchaya, G.; Mabhaudhi, T.; Nhlengethwa, S.; Nhemachena, C.; Mpandeli, S. Cereal production trends under climate change: Impacts and adaptation strategies in southern Africa. Agriculture 2019, 9, 30. [Google Scholar] [CrossRef] [Green Version]
  78. Liu, W.; Kogan, F. Monitoring regional drought using the vegetation condition index. Int. J. Remote Sens. 1996, 17, 2761–2782. [Google Scholar] [CrossRef]
  79. Mavromatis, T. Drought index evaluation for assessing future wheat production in Greece. Int. J. Climatol. A J. R. Meteorol. Soc. 2007, 27, 911–924. [Google Scholar] [CrossRef]
  80. Wu, H.; Hubbard, K.G.; Wilhite, D.A. An agricultural drought risk-assessment model for corn and soybeans. Int. J. Climatol. A J. R. Meteorol. Soc. 2004, 24, 723–741. [Google Scholar] [CrossRef] [Green Version]
  81. Mohammed, S.; Alsafadi, K.; Enaruvbe, G.O.; Bashir, B.; Elbeltagi, A.; Széles, A.; Alsalman, A.; Harsanyi, E. Assessing the impacts of agricultural drought (SPI/SPEI) on maize and wheat yields across Hungary. Sci. Rep. 2022, 12, 8838. [Google Scholar] [CrossRef]
  82. Sharma, A.; Goyal, M.K. Assessment of ecosystem resilience to hydroclimatic disturbances in India. Glob. Chang. Biol. 2018, 24, e432–e441. [Google Scholar] [CrossRef]
  83. Baudoin, M.-A.; Vogel, C.; Nortje, K.; Naik, M. Living with drought in South Africa: Lessons learnt from the recent El Niño drought period. Int. J. Disaster Risk Reduct. 2017, 23, 128–137. [Google Scholar] [CrossRef]
  84. Nicholson, S.E.; Kim, J. The relationship of the El Niño–Southern oscillation to African rainfall. Int. J. Climatol. A J. R. Meteorol. Soc. 1997, 17, 117–135. [Google Scholar] [CrossRef]
  85. Chikoore, H.; Jury, M.R. Intraseasonal variability of satellite-derived rainfall and vegetation over Southern Africa. Earth Interact. 2010, 14, 1–26. [Google Scholar] [CrossRef]
  86. Mahlalela, P.; Blamey, R.; Hart, N.; Reason, C. Drought in the Eastern Cape region of South Africa and trends in rainfall characteristics. Clim. Dyn. 2020, 55, 2743–2759. [Google Scholar] [CrossRef]
  87. Lakhraj-Govender, R.; Grab, S.W. Assessing the impact of El Niño–Southern Oscillation on South African temperatures during austral summer. Int. J. Climatol. 2019, 39, 143–156. [Google Scholar] [CrossRef] [Green Version]
  88. Botai, C.M.; Botai, J.O.; Adeola, A.M. Spatial distribution of temporal precipitation contrasts in South Africa. S. Afr. J. Sci. 2018, 114, 70–78. [Google Scholar] [CrossRef]
  89. Hoffman, M.T.; Carrick, P.; Gillson, L.; West, A. Drought, climate change and vegetation response in the succulent karoo, South Africa. S. Afr. J. Sci. 2009, 105, 54–60. [Google Scholar] [CrossRef] [Green Version]
  90. Rojas, O.; Vrieling, A.; Rembold, F. Assessing drought probability for agricultural areas in Africa with coarse resolution remote sensing imagery. Remote Sens. Environ. 2011, 115, 343–352. [Google Scholar] [CrossRef]
  91. Haile, G.G.; Tang, Q.; Leng, G.; Jia, G.; Wang, J.; Cai, D.; Sun, S.; Baniya, B.; Zhang, Q. Long-term spatiotemporal variation of drought patterns over the Greater Horn of Africa. Sci. Total Environ. 2020, 704, 135299. [Google Scholar] [CrossRef]
  92. Kasei, R.; Diekkrüger, B.; Leemhuis, C. Drought frequency in the Volta basin of West Africa. Sustain. Sci. 2010, 5, 89–97. [Google Scholar] [CrossRef]
  93. Amri, R.; Zribi, M.; Lili-Chabaane, Z.; Duchemin, B.; Gruhier, C.; Chehbouni, A. Analysis of vegetation behavior in a North African semi-arid region, using SPOT-VEGETATION NDVI data. Remote Sens. 2011, 3, 2568–2590. [Google Scholar] [CrossRef] [Green Version]
  94. Mbatha, N.; Xulu, S. Time Series Analysis of MODIS-Derived NDVI for the Hluhluwe-Imfolozi Park, South Africa: Impact of Recent Intense Drought. Climate 2018, 6, 95. [Google Scholar] [CrossRef] [Green Version]
  95. Musonda, B.; Jing, Y.; Iyakaremye, V.; Ojara, M. Analysis of long-term variations of drought characteristics using standardized precipitation index over Zambia. Atmosphere 2020, 11, 1268. [Google Scholar] [CrossRef]
  96. Hadri, A.; Saidi, M.E.M.; Boudhar, A. Multiscale drought monitoring and comparison using remote sensing in a Mediterranean arid region: A case study from west-central Morocco. Arab. J. Geosci. 2021, 14, 118. [Google Scholar] [CrossRef]
  97. Desbureaux, S.; Damania, R. Rain, forests and farmers: Evidence of drought induced deforestation in Madagascar and its consequences for biodiversity conservation. Biol. Conserv. 2018, 221, 357–364. [Google Scholar] [CrossRef]
  98. Uwimbabazi, J.; Jing, Y.; Iyakaremye, V.; Ullah, I.; Ayugi, B. Observed Changes in Meteorological Drought Events during 1981–2020 over Rwanda, East Africa. Sustainability 2022, 14, 1519. [Google Scholar] [CrossRef]
  99. Black, E.; Tarnavsky, E.; Maidment, R.; Greatrex, H.; Mookerjee, A.; Quaife, T.; Brown, M. The use of remotely sensed rainfall for managing drought risk: A case study of weather index insurance in Zambia. Remote Sens. 2016, 8, 342. [Google Scholar] [CrossRef] [Green Version]
  100. Reddy, S.S. IPCC Special Report Emission Scenarios: Summary for Plicy Makers: A Speical Report of Working Group III of Intergovernmental Panel of Climatic Change; Cambridge University: Cambridge, UK, 2000. [Google Scholar]
  101. Ding, Y.; Hayes, M.J.; Widhalm, M. Measuring economic impacts of drought: A review and discussion. Disaster Prev. Manag. Int. J. 2011, 20, 434–446. [Google Scholar] [CrossRef] [Green Version]
  102. Ziolkowska, J.R. Socio-Economic Implications of Drought in the Agricultural Sector and the State Economy. Economies 2016, 4, 19. [Google Scholar] [CrossRef] [Green Version]
  103. Rouault, M.; Richard, Y. Intensity and spatial extent of droughts in southern Africa. Geophys. Res. Lett. 2005, 32, L15702. [Google Scholar] [CrossRef]
  104. Rulinda, C.M.; Dilo, A.; Bijker, W.; Stein, A. Characterising and quantifying vegetative drought in East Africa using fuzzy modelling and NDVI data. J. Arid. Environ. 2012, 78, 169–178. [Google Scholar] [CrossRef]
  105. USDA. Wheat Production in South Africa. Available online: https://ipad.fas.usda.gov/highlights/2015/09/RSA/index.htm (accessed on 27 April 2022).
  106. Unganai, L.S.; Kogan, F.N. Drought monitoring and corn yield estimation in Southern Africa from AVHRR data. Remote Sens. Environ. 1998, 63, 219–232. [Google Scholar] [CrossRef]
  107. Senay, G.B.; Verdin, J. Characterization of yield reduction in Ethiopia using a GIS-based crop water balance model. Can. J. Remote Sens. 2003, 29, 687–692. [Google Scholar] [CrossRef]
  108. Rosenzweig, C.; Tubiello, F.N. Adaptation and mitigation strategies in agriculture: An analysis of potential synergies. Mitig. Adapt. Strateg. Glob. Chang. 2007, 12, 855–873. [Google Scholar] [CrossRef] [Green Version]
  109. Rosenzweig, C.; Karoly, D.; Vicarelli, M.; Neofotis, P.; Wu, Q.; Casassa, G.; Menzel, A.; Root, T.L.; Estrella, N.; Seguin, B. Attributing physical and biological impacts to anthropogenic climate change. Nature 2008, 453, 353–357. [Google Scholar] [CrossRef]
  110. Zhang, W.-f.; Dou, Z.-x.; He, P.; Ju, X.-T.; Powlson, D.; Chadwick, D.; Norse, D.; Lu, Y.-L.; Zhang, Y.; Wu, L. New technologies reduce greenhouse gas emissions from nitrogenous fertilizer in China. Proc. Natl. Acad. Sci. USA 2013, 110, 8375–8380. [Google Scholar] [CrossRef]
  111. Holland, J.M. The environmental consequences of adopting conservation tillage in Europe: Reviewing the evidence. Agric. Ecosyst. Environ. 2004, 103, 1–25. [Google Scholar] [CrossRef]
  112. Mishra, A.K.; Singh, V.P. A review of drought concepts. J. Hydrol. 2010, 391, 202–216. [Google Scholar] [CrossRef]
  113. Wilhite, D.A. National drought management policy guidelines: A template for action. In Integrated Drought Management Programme (IDMP) Tools and Guidelines Series; WMO: Geneva, Switzerland, 2014; Volume 1. [Google Scholar]
  114. Faiz, M.A.; Zhang, Y.; Zhang, X.; Ma, N.; Aryal, S.K.; Ha, T.T.V.; Baig, F.; Naz, F. A composite drought index developed for detecting large-scale drought characteristics. J. Hydrol. 2022, 605, 127308. [Google Scholar] [CrossRef]
Figure 1. Location of South Africa, along with the distribution of the meteorological stations.
Figure 1. Location of South Africa, along with the distribution of the meteorological stations.
Ijerph 19 16469 g001
Figure 2. Boxplots of temperature and rainfall for all stations across South Africa (2000–2020): (a) temperature (°C); (b) rainfall (mm); (☐) mean; (____) median.
Figure 2. Boxplots of temperature and rainfall for all stations across South Africa (2000–2020): (a) temperature (°C); (b) rainfall (mm); (☐) mean; (____) median.
Ijerph 19 16469 g002
Figure 3. Evolution of wheat yield. Evolution of wheat yield (ton/hectare) between 2000 and 2020 in South Africa on a provincial level: (a) WES-C; (b) NR-C; (c) FS; (d) ES-C; (e) KZN; (f) MP; (g) LPP; (h) GG; (i) NW.
Figure 3. Evolution of wheat yield. Evolution of wheat yield (ton/hectare) between 2000 and 2020 in South Africa on a provincial level: (a) WES-C; (b) NR-C; (c) FS; (d) ES-C; (e) KZN; (f) MP; (g) LPP; (h) GG; (i) NW.
Ijerph 19 16469 g003
Figure 4. Drought (SPI) trend between 2002 and 2021 in South Africa on different time scales: 3 months; 6 months; 9 months; and 12 months.
Figure 4. Drought (SPI) trend between 2002 and 2021 in South Africa on different time scales: 3 months; 6 months; 9 months; and 12 months.
Ijerph 19 16469 g004
Figure 5. SPI-3 values on a monthly scale across South African provinces between 2002 and 2021: (a) WES-C; (b) NR-C; (c) FS; (d) ES-C; (e) KZN; (f) MP; (g) LPP; (h) GG; (i) NW (red dots indicate SPI-3 values of less than −1).
Figure 5. SPI-3 values on a monthly scale across South African provinces between 2002 and 2021: (a) WES-C; (b) NR-C; (c) FS; (d) ES-C; (e) KZN; (f) MP; (g) LPP; (h) GG; (i) NW (red dots indicate SPI-3 values of less than −1).
Ijerph 19 16469 g005
Figure 6. Frequency of SPI-3 drought occurrence at provincial scale between 2002 and 2021: (1) WES-C; (2) NR-C; (3) FS; (4) ES-C; (5) KZN; (6) MP; (7) LPP; (8) GG; (9) NW.
Figure 6. Frequency of SPI-3 drought occurrence at provincial scale between 2002 and 2021: (1) WES-C; (2) NR-C; (3) FS; (4) ES-C; (5) KZN; (6) MP; (7) LPP; (8) GG; (9) NW.
Ijerph 19 16469 g006
Table 1. Full overview of the analyzed meteorological stations in South Africa (2000–2021).
Table 1. Full overview of the analyzed meteorological stations in South Africa (2000–2021).
StationAbbreviationXYElevation (m) R * (mm)ProvinceAbbreviation
BetlehemBET−28.249628.3343168858.25Free StateFS
BloemfointeinBLM−29.120426.1874135443
CalviniaCAL−31.481919.761797514.54Western CapeWES-C
Cape ColumbineCCE−32.827817.85586719.67
Cape TownCPT−33.963118.60234238.05
UpingtonUPN−28.411121.264184819.35
De AarDAR−30.665123.9927124726.58Northern CapeNR-C
KimberleyKMBLY−28.806124.7698119833.52
SpringbokSPB−29.669417.8788100617.24
East LondonEL−33.035727.816112564.01Eastern CapeES-C
GrahamstownGHT−33.290726.502664242.11
uMthathaMTT−31.549728.673974254.68
Port AlfredPTA−33.559526.88093751.97
Port ElizabethPE−33.986425.61646149.85
RusternburgRSB−25.660727.2322115730.63North-WestNW
LichtenburgLCBG−26.13326.1644148748.5
MafikengMFG−25.803725.5428127939.57
LadysmithLYS−28.575529.7503107852.19KwaZulu NatalKZN
MntunziniMTN−28.947431.70796592.8
PietermariesburgPTMBG−29.627830.402967359.57
PotchersroomPOR−26.735927.0755135143.94
Richards BayRDB−28.737832.0934889.62
JHB Bot TuineJHB BT−26.156627.9991162648.24GautengGG
JHB Int WOJHB Int WO−26.14328.2346169458.63
IreneIRE−25.910528.2106152355.31
PolokwanePOL−23.857629.4517122854.29LimpopoLPP
ThohoyandouTHY−22.984530.458361859.11
LephalaleLPL−23.676727.705184031.92
SkukuzaSKZ−24.992631.58827142.98MpumalangaMP
OudestadODD−25.1829.3394933.79
ErmeloERM−26.497729.9838173758.96
* R (mm): average monthly rainfall for the whole period.
Table 2. Classification of standard precipitation index (SPI) for drought studies.
Table 2. Classification of standard precipitation index (SPI) for drought studies.
SPI ValueCategory
−1.0 ≤ SPI ≤ 0Mild event
−1.49 < SPI ≤ −1.0Moderate event
−2.0 < SPI ≤ −1.5Severe event
SPI ≤ −2.0Extreme event
Table 3. Classification of SYRS and CR.
Table 3. Classification of SYRS and CR.
SYRS ValuesSYRS ClassesCR ValuesCR Classes
−0.5  <  SYRS  ≤  0.5Normal conditionsCR > 1Resilient
−0.5  <  SYRS  ≤  −1.0Acceptable losses due to drought 0.9 < CR < 1Slightly nonresilient
−1.0  <  SYRS  ≤  −1.5Moderate0.8 < CR < 0.9Moderately nonresilient
−1.5  <  SYRS  <  −2.0HighCR < 0.8Severely nonresilient
SYRS  ≤  −2.0Extreme
Table 4. Detailed trend of drought at the studied stations (2002–2021) based on MK test and Sen’s slope estimator.
Table 4. Detailed trend of drought at the studied stations (2002–2021) based on MK test and Sen’s slope estimator.
StationSPI-3SPI-6SPI-9SPI-12
pSen’s *pSen’spSen’spSen’s
BET0.590.010.590.0100.0300.03
BLM0.19−0.010.19−0.010.9400.610
CAL0.03−0.020.03−0.020−0.030−0.03
CCE0−0.030−0.03<0.0001−0.04<0.0001−0.05
CPT<0.0001−0.04<0.0001−0.04<0.0001−0.05<0.0001−0.06
DAR0.31−0.010.31−0.010.01−0.020.01−0.02
EL0.14−0.010.14−0.01<0.0001−0.04<0.0001−0.04
ERM0.9200.9200.40.010.290.01
GHT0.18−0.010.18−0.010−0.030.01−0.03
IRE0.010.020.010.02<0.00010.04< 0.00010.05
JHB BT00.0300.03<0.00010.05< 0.00010.05
JHB Int WO0.08−0.020.08−0.020−0.030−0.03
KMBLY0.04−0.020.04−0.020.35−0.010.44−0.01
LYS0.31−0.010.31−0.010.9400.660
LPL0.8500.8500.7100.610
LCBG0.6500.6500.120.020.070.02
MFG0.9800.9800.52−0.010.35−0.01
MTT0.35−0.010.35−0.010.22−0.010.11−0.01
MTN0.8900.8900.38−0.010.3−0.01
ODD0.2−0.010.2−0.010.03−0.010.04−0.01
PTMBG<0.0001−0.05<0.0001−0.05<0.0001−0.08<0.0001−0.09
POL0.25−0.010.25−0.010.12−0.010.08−0.02
PTA0.05−0.020.05−0.02<0.0001−0.030−0.02
PE0.03−0.020.03−0.02<0.0001−0.04<0.0001−0.04
POR<0.00010.04<0.00010.04<0.00010.06<0.00010.06
RDB0.47−0.010.47−0.010.05−0.020.09−0.01
RSB<0.0001−0.05<0.0001−0.05<0.0001−0.07<0.0001−0.08
SKZ0.01−0.030.01−0.03<0.0001−0.05<0.0001−0.05
SPB0−0.040−0.04<0.0001−0.05<0.0001−0.06
THY0.8200.8200.310.010.090.02
UPN0.07−0.020.07−0.020.08−0.010.24−0.01
* Sen’s: Sen’s slope estimator at decadal scale; gray shade indicates a significant (p < 0.05) trend.
Table 5. Year-by-year comparison of SYRSSPI-3 changes in South African provinces from 2002 to 2021.
Table 5. Year-by-year comparison of SYRSSPI-3 changes in South African provinces from 2002 to 2021.
YearWES-CNR-CFSES-CKZNMPLPPGGNW
20020.960.080.19−0.850.51.35−0.55−0.12−0.62
2003−0.820.3−0.370.15−0.33−2.370.06−1.890.03
2004−1.280−0.580.371.770.52−1.010.740.36
2005−0.10.55−0.570.61−0.680.24−0.140.270.93
20060.49−0.110.47−1.79−1.72−0.02−0.55−1.59−1.95
20070.36−0.491.191.28−0.55−0.761.601.66−0.84
20080.16−0.01−0.210.290.861.321.092.52
2009−0.11−1.290.790.65−0.430.131.020.820.6
2010−0.98−0.77−0.910.25−0.79−1.09−0.010.59−0.29
20110.262.41−0.060.380.431.750.881−0.11
20121.39−1.780.27−0.39−1.56−1.09−1.82−1.63−0.45
20130.690.880.42−0.112.041.15−0.83−0.30.5
20140.420.351.1−1.891.620.47−1.64−0.270.27
2015−0.96−0.61−1.67−0.25−0.35−1.46−1.05−0.82−1.21
20161.250.28−0.91−0.39−0.47−0.430.53−0.77−0.89
2017−21.721.31−1.180.670.361.50−0.670.97
2018−0.080.480.292.340.590.520.740.680.32
2019−1.74−1.49−2.52−0.43−0.3−0.37−0.51−0.33−1
20200.96−0.590.490.49−0.41−0.040.230.63−0.31
20211.140.091.26−0.23−0.330.30.220.911.18
Table 6. Correlation between SYRS and SPI-3 on a monthly time scale and on a regional scale (2002–2021).
Table 6. Correlation between SYRS and SPI-3 on a monthly time scale and on a regional scale (2002–2021).
ProvinceJANFEBMARAPRMAYJUNJULAUGSEPOCTNOVDEC
WES-C0.160.30.530.40.480.380.550.550.520.320.1−0.05
NR-C0.540.320.06−0.050.150.060.01−0.22−0.23−0.13−0.03−0.24
FS0.390.38−0.11−0.33−0.6−0.3−0.240.030.130.540.50.54
ES-C−0.170.020.14−0.11−0.22−0.27−0.04−0.22−0.36−0.39−0.23−0.12
KZN−0.180.09−0.070.02−0.14−0.090.03−0.11−0.14−0.2−0.23−0.18
MP0.310.250.390.320.410.120.150.12−0.150.14−0.060.22
LPP−0.05−0.18−0.19−0.22−0.20.120.530.590.250.040.250.05
GG0.05−0.11−0.010.130.350.510.35−0.03−0.320.150.240.39
NW0.240.140.260.120.05−0.04−0.05−0.24−0.45−0.12−0.12−0.06
Table 7. Shows the percentage of wheat yield losses (YLs, %) caused by drought occurrences (drought event, DE) at the provincial level over a 3-month time scale during the winter wheat stages *.
Table 7. Shows the percentage of wheat yield losses (YLs, %) caused by drought occurrences (drought event, DE) at the provincial level over a 3-month time scale during the winter wheat stages *.
Western Cape (WES-C), CR = 0.65Northern Cape (NR-C), CR = 0.93Free State (FS), CR = 0.65
YearYL%DEGSDS DDYearYL%DEGSDS DDYearYL%DEGSDSDD
2017−352016Aug–18MarGC13.3202012−8.7-nnn2019−34.52019Jul–20JanGP, HP6.17
2019−29.52018Aug–20JanGC14182019−6.82019Jul–20MarGC8.8102015−28.72015Feb–16AprSP, HP10.511
2004−29.32004Mar–4JulSP1.432009−6.42009May–9JunSP0.622010−20.42010Jul–10OctGP, HP3.94
2003−18.62003Mar–3AugSP, GP6.66 2004−172004Jun–4DecGP, HP2.35
2010−18.62010Mar–10OctGP, HP2.94 2005−15.92005Ju1–4SeptGP1.93
2015−16.52015Feb–15decGC9.411 2003−10.82002Nov–4MarGC10.817
Eastern Cape (ES-C), CR = 0.97KwaZulu−Natal (KZN), CR = 0.96Mpumalanga (MP), CR = 0.92
YearYL%DEGSDSDDYearYL%DEGSDSDDYearYL%DEGSDSDD
2006−17.8-GP-HPww2006−8.4-nnn2003−15.52003Sept–3DecGP, HP1.94
2014−142014Jul–4DecGP, HP4.262012−6.92012May–12AugSP, GP2.342015−7.52014Jun–15AugSP, GP1215
2017−7.52016Oct–17SpetSP, GP6122010−3.42010Jul–10NovGP, HP4.752010−62010Ju1–11AugGP, HP3.15
2002−8.8-GP, HPww 2012−5.72012Jun–12AugSP, GP2.83
2019−2.52019Oct–20JanGP, HP7.46 2007−4.22007Jan–12JunSP4.56
Limpopo (LPP), CR = 0.79Gauteng (GG), CR = 0.85North−West (NW) CR = 0.97
YearYL%DEGSDSDDYearYL%DEGSDSDDYearYL%DEGSDSDD
2012−20.32012Jan–12AugSP, GP682003−14.52003Sept–4JanGC652006−7.2nnn
2014−18.32014Jun–14NovGP, HP3.762006−10.82006May–6JulSP, GP2.532015−42015Jan–15JunSP3.86
2004−17.62004Jun–05OctGC10.6182012−102012Jan–12AugSP, GP6.582019−3.22019Jul–19NovGP, HP5.85
2015−14.22015Mar–15AugSP, GP2.462015−4.62015Feb–15JulSP, GP2.46
2005−10.22004Jun–05OctGC10.618
* Winter wheat growing stages (GSs) include the sowing period (SP), which runs from May to June; the growing period (GP), which runs from June to September; the harvesting period (HP), which runs from September to November; and the entire growing cycle (GC), which runs from May to November. DD is drought duration, DS is drought severity, and gray shading is used to represent the YL (%) during wet events (w) or no drought event occurrences (n).
Table 8. Previous studies on drought trends across Africa and in South Africa.
Table 8. Previous studies on drought trends across Africa and in South Africa.
BasinCountryDrought IndicesPeriodOutput and ImpactReference
Eastern Cape province South Africa Rainfall trend 1981–2018Drought was detected in all seasons since 2015Mahlalela et al. [86]
Karroo (Northern, Western, and Eastern Cape provinces)South AfricaSPI1900–2000From 1900 to 1950, dry spell patterns were detected, whereas no visible precipitation or drought trends in 1951–2000 periodHoffman et al. [89]
AfricaAfrican countriesRemote sensing imagery (NDVI and HVI)1981–2009Impact of drought on agricultureRojas et al. [90]
Greater Horn of AfricaEthiopia, Eritrea, Kenya, Rwanda, Somalia, Sudan, Uganda, Tanzania, Burundi, Djibouti, and South SudanSPEI1964–2015The previous 52 years saw rising trends in drought with different temporal and spatial patternsHaile et al. [91]
West Africa (Volta Basin)Benin, Togo, Mali, Ghana, Burkina Faso, and Ivory CoastSPI1961–2003Frequency of droughts has increased since the 1970sKasei et al. [92]
Kairouan plain Central Tunisia SPOT VEGETATION NDVI1998–2010Reveals drought year in 2000–2001 and decadal persistent temporal fluctuations in agricultureAmri et al. [93]
Hluhluwe–iMfolozi ParkKZN, South AfricaNDVI, EVI, BAI, and NDII2002–2017The indices show the vegetation experienced water stress, especially in 2003 and 2014–2016Mbatha and Xulu [94]
ZambiaAfricaSPI1981–2017Widespread floods in DJF (summer) seasons and drought periods in 1992, 1995, and 2005Musonda et al. [95]
Chichaoua–MejjateMoroccoSPI and NDVI2008–2017Trend in temporal persistent soaring drought with exception on 3-month SPI scaleHadri et al. [96]
AfricaMadagascarSPI1900–2013Result of drought-induced deforestation and loss of biodiversityDesbureaux and Damania [97]
AfricaRwandaSPEI and SPI1981–2020Much variability in rainfall and temperature with a major decline in 2010–2017 and drought events in 2015, 2016, and 2017Uwimbabazi et al. [98]
Southern AfricaZambiaJoint UK Land Environment Simulator (JULES)1995–2009Strong relationships exist between drought classifications in all spatial ranges in south, west, and east regions of ZambiaBlack et al. [99]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Nxumalo, G.; Bashir, B.; Alsafadi, K.; Bachir, H.; Harsányi, E.; Arshad, S.; Mohammed, S. Meteorological Drought Variability and Its Impact on Wheat Yields across South Africa. Int. J. Environ. Res. Public Health 2022, 19, 16469. https://doi.org/10.3390/ijerph192416469

AMA Style

Nxumalo G, Bashir B, Alsafadi K, Bachir H, Harsányi E, Arshad S, Mohammed S. Meteorological Drought Variability and Its Impact on Wheat Yields across South Africa. International Journal of Environmental Research and Public Health. 2022; 19(24):16469. https://doi.org/10.3390/ijerph192416469

Chicago/Turabian Style

Nxumalo, Gift, Bashar Bashir, Karam Alsafadi, Hussein Bachir, Endre Harsányi, Sana Arshad, and Safwan Mohammed. 2022. "Meteorological Drought Variability and Its Impact on Wheat Yields across South Africa" International Journal of Environmental Research and Public Health 19, no. 24: 16469. https://doi.org/10.3390/ijerph192416469

APA Style

Nxumalo, G., Bashir, B., Alsafadi, K., Bachir, H., Harsányi, E., Arshad, S., & Mohammed, S. (2022). Meteorological Drought Variability and Its Impact on Wheat Yields across South Africa. International Journal of Environmental Research and Public Health, 19(24), 16469. https://doi.org/10.3390/ijerph192416469

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop