Year-to-year fluctuations of Ethiopia climate are investigated to develop statistical predictions at one-season lead time. Satellite vegetation data from NASA and rainfall from ARC2 are the basis for analysis. The “target” seasons are May–July and August–October, while “predictors” are December–February and March–May, respectively. Global fields of surface temperature, sea level air pressure, and upper and lower level zonal winds are employed in point-to-field correlations. After step-wise multivariate regression, the leading predictors are: surface temperature across Europe (cold-favourable), 850 mb zonal winds over the tropical Atlantic (easterly-favourable), and surface temperature in the tropical Indian Ocean (cold-favourable). Predictive algorithms for early and late rainfall exhibit a consistent r2 fit of ~0.50, while those for vegetation reach ~0.65 in late summer, indicating that fluctuations in food resources could be forewarned.
1. Introduction
Climate variability in Ethiopia directly impacts food resources available to a rural population already experiencing a deficit of ~1013 KCAL/yr [1]. An ability to predict crop yield and improve production requires a technical capacity underpinned with knowledge on how the global circulation affects regional climate. Reliable locally tailored forecasts can help communities to avoid risks, optimize resources, and restrain the macroeconomic impacts of drought and flood [2–8]. These prospects motivate our research.
When climate departs from the expected annual cycle, it is often due to slowly varying surface boundary conditions like the Pacific El Nino Southern Oscillation (ENSO) or its decadal component (PDO). While the tropical ocean thermocline is considered a driver of climate, the atmosphere’s hydrostatic response may cause long-lived circulation patterns to slowly propagate around the world, bringing opportunities for climate prediction [9–12]. Coupled ensemble numerical models are becoming more sophisticated in simulating these processes, and their forecasts already provide useful regional guidance a few months before the rainy season [13–15]. Wang and Fan [16] show how both ensemble numerical model predictions and observed spatial patterns in “analog years” can be combined to improve the prediction of east Asian summer rainfall. Yet there is still room for statistical techniques that assume historical replication via multiple regression and more complex techniques [17–22].
The main goals of this study are to develop the ability to statistically predict Ethiopian summer rainfall and vegetation by determining the most influential climate indices, formulating multivariate algorithms and evaluating their reliability. Following this introduction, Section 2 provides the data and methods. Section 3 analyzes the results, while conclusions are given in Section 4.
2. Data and Methods
The data and methods employed to develop prediction algorithms for Ethiopia are described. “Targets” include interpolated observed summer rainfall and satellite vegetation fraction, while surface and atmospheric “predictors” are drawn from satellite-era reanalysis products in the preceding season. An earlier study found a significant correlation (r=0.56) between colocated vegetation fraction in July–September and harvested maize and sorghum yields in the period 1986–2009 [23].
2.1. Target and Predictor Data
Gridded monthly rainfall data for Ethiopia were analyzed from the African Rainfall Climatology (ARCv2, [24]). A secondary dataset was derived from the Global Precipitation Climatology Center (GPCCv6, [25, 26]). Their respective resolutions/temporal coverages are 10 km/1983–2012 and 50 km/1981–2010. Comparisons with gauges from the Ethiopian National Meteorological Agency yielded better performance by ARC2. Vegetation fraction values originate from the corrected NASA satellite dataset (NDVI, [27]) at 25 km resolution in the period 1981–2006. Instead of using objective analysis to define target areas [28], crop productivity reports from the Ethiopian Central Statistical Agency (http://www.csa.gov.et/) suggest two main zones in the northern and southern highlands: 10–14N, 36–40E~Amhara and 6–10N, 35–40E~Oromia, identified in Figure 1(a). The rainy season is divided into the first and second half (May–July, August–October) hereafter MJJ and ASO. The targets are of sufficient size for averaging to damp out local “noise”, yet they still fall within a homogeneous climate regime [29]. Considering the mean annual cycle (Figure 1(b)), rainfall in the south rises in May while vegetation in the north peaks in September. Hence early (late) summer forecasts are critical to strategic planning in Oromia (Amhara). Linear trend analysis indicates no significant drift in the target records over the period since 1981 (Figure 2). Target data were drawn from the IRI Climate Library (ARC2) and Climate Explorer websites (GPCC6, NDVI).
(a) Vegetation fraction averaged 2000–2012 from MODIS, with state boundaries, target boxes and the national experimental farm (arrow). (b) Mean annual cycle of ARC2/GPCC6 rainfall and vegetation NDVI 1981–2006.
Standardized departures of (a) rainfall and (b) vegetation fraction for two zones in two seasons (MJJ NDVI omitted in north).
The predictor fields employed here include global land-sea surface air temperature from the National Climate Data Center v3 (Ts, [30]), sea level air pressure from the Hadley Center v2 reanalysis (SLP, [31]), and 850 mb and 200 mb zonal winds from Coupled Forecast System (CFS) and NASA Modern Era (MERRA) reanalysis [32, 33]. Seasonal averages were calculated: December–February (DJF) for MJJ targets and March–May (MAM) for ASO targets. Predictor fields were available from the Climate Explorer website 〈http://climexp.knmi.nl/〉 and analyzed globally in latitudes 50S–60N.
2.2. Methods
Our search for predictors was facilitated by linear point-to-field correlation analysis followed by stepwise multivariate regression of key-area time series. Correlation maps were calculated for the specified predictor-target lead time (e.g., DJF for MJJ, MAM for ASO) and fields were masked below 80% significance. For selection as a predictor, the key area should exceed 15° latitude × 20° longitude in size. All time series were converted to standardized departures (cf. Figure 2).
Statistical prediction algorithms were formulated via backward stepwise linear regression onto the target time series. The candidate predictor pool was <10 compared with a training period of ~30 years. Initially all predictors were included and their partial correlation was evaluated. Those with lower significance (or colinearity) were screened out and the algorithm was recalculated from the remaining variables. In most cases an optimal fit was reached with three predictors, thereby limiting the chance of artificial skill [6].
A 30 year training period does not readily support independent validation, so the performance of multivariate linear algorithms was evaluated by r2 fit, adjusted for the number of predictors. Our target time series exhibit minimal persistence (Figures 2(a) and 2(b)) yielding a degrees of freedom ~30. A predictive algorithm with r2 fit >0.50 indicates that cost-effective tercile forecasts (above/normal/below) can be achieved. Forecast versus observed scatterplots were analyzed for slope, tercile hits, and outliers.
3. Results3.1. Targets and Correlation Maps
The regression of rain and vegetation time series (Figures 2(a) and 2(b)) onto global fields of Ts, SLP, and U wind is calculated for the prescribed lead time, season, and area. The targets reflect biennial and decadal oscillations (cf. [29, 34]) and a more coherent structure for rainfall than vegetation. Strangely, year-to-year fluctuations of MJJ rain correlate significantly with colocated MJJ vegetation (~0.4) but not in ASO season. Many candidate predictors emerge in the correlation maps (cf. [23]), but only a few are needed to form multivariate algorithms for northern Amhara and southern Oromia zones.
For MJJ northern rainfall, the selected DJF predictors (Figures 3(a) and 3(b); Table 1(a)) are surface temperatures over northern Europe (with a coefficient of −0.38) and the central Indian Ocean (−0.39) and 850 mb zonal winds over the central Atlantic (−0.44) and South Indian Ocean (0.42). For ASO northern rainfall (Figures 3(c)–3(e)) the selected MAM predictors are surface temperature over Asia (−0.45), 200 mb zonal winds over the northwest Pacific (0.50), and 850 mb zonal winds over western Europe (0.47) and east Africa (−0.45). For MJJ southern rainfall (Figures 4(a) and 4(b)), the selected DJF predictors are surface temperature over northern Europe (−0.70) and 850 mb zonal winds over the central Atlantic (−0.56). For ASO southern (GPCC6) rainfall (Figures 4(c) and 4(d)) the selected MAM predictors are surface temperature over the southeast Atlantic (−0.28) and 200 mb zonal winds over South America (0.21) and the southwest Atlantic (−0.56). Arrows on the maps highlight their position, while Table 1(a) lists their domain. A key feature is negative Ts correlations in northern Eurasia (cf. Figures 3(a), 3(c), and 4(a)) indicating that a cold winter-spring there precedes favourable Ethiopian rainfall. This signal is part of the Arctic Oscillation and its quasi-biennial circulation anomalies ([35–37]) that affect Ethiopia as outlined in Jury [23]. The correlation dipole over the southwest Atlantic (cf. Figure 4(c)) refers to an equatorward-displaced subtropical jet that plays a role in the ASO southern rain algorithm. The Pacific ENSO-PDO pattern of strengthened equatorial trade winds (Figure 4(b)) is not selected as a predictor.
Multivariate predictors for (a) rain and (b) vegetation.
Target
Predictor
Domain
Coefficient
MJJ northern rain
Ts n Europe
45N–65N, 40W–70E
-0.38
r2=0.51
Ts c Indian O
20S–10N, 70E–95E
-0.39
Offset = 0.01
8 U c Atlantic
10S–10N, 60W–10W
-0.44
8 U s Indian O
45S–30S, 40E–110E
0.42
ASO northern rain
Ts Asia
40N–60N, 50E–130E
-0.45
r2=0.53
2 U nw Pacific
45N–60N, 140E–160W
0.50
Offset = 0.01
8 U w Europe
35N–50N, 50W–20E
0.47
8 U e Africa
15S–0N, 35E–55E
-0.45
MJJ southern rain
Ts n Europe
45N–60N, 20W–80E
-0.70
r2=0.49, offset = 0.03
8 U c Atlantic
15S–10N, 70–20W
-0.56
ASO southern rain
Ts se Atlantic
20S–0N, 15W–15E
-0.28
r2=0.51
2 U s America
32S–15S, 120W–0W
0.21
Offset = 0.05
2 U sw Atlantic
52S–40S, 80W–0W
-0.56
Target*
Predictor
Domain
Coefficient
ASO northern Veg
Ts sw Pacific
30S–12N, 140E–170E
0.59
r2=0.67
SLP n Atlantic
48N–62N, 80W–0W
0.38
SLP s Indian O
50S–30S, 65E–110E
-0.34
2 U n Africa
22N–35N, 30E–40E
-0.31
MJJ southern Veg
Ts w Indian O
18S–18N, 5W–100E
-0.34
r2=0.43
Ts se Pacific
38S–18S, 110W–75W
-0.47
2 U se Pacific
32S–15S, 125W–70W
0.28
ASO southern Veg
Ts sw Atlantic
50S–32S, 50W–10W
-0.26
r2=0.62
SLP Europe
40N–60N, 25W–15E
-0.35
SLP s Indian O
50S–30S, 60E–90E
-0.50
2 U ne Pacific
20N–40N, 140W–110W
-0.60
*MJJ northern omitted.
Correlation of MJJ northern rain with DJF predictors: (a) Ts and (b) 850 U fields. Correlation of ASO northern rain with MAM predictors: (c) Ts, (d) 200 U, and (e) 850 U fields 1983–2012. Values below 80% significance are masked here and in Figures 4–6; arrows highlight predictors employed.
Correlation of MJJ southern rain with DJF (a) Ts and (b) 850 U fields 1983–2012. Correlation of ASO southern rain with MAM predictors: (c) Ts, and (d) 200 U fields for GPCC6 rain 1981–2010.
Vegetation in the northern zone rises slowly after May (cf. Figure 1(b)) so its early season variability is not considered here. Maps for ASO northern vegetation (Figures 5(a)–5(c)) point to the following MAM predictors are surface temperature over the southwest Pacific (0.59), sea level air pressure over the north Atlantic (0.38) and south Indian Ocean (−0.34), and the 200 mb zonal wind over North Africa (−0.31). The later predictor suggests influence from the Arctic Oscillation (AO) and North Atlantic Oscillation (NAO). For MJJ southern vegetation (Figures 6(a) and 6(b)) the DJF predictors are surface temperature in the west Indian Ocean (−0.34) and the southeast Pacific (−0.47), and 200 mb zonal winds over the southeast Pacific (0.28). For ASO southern vegetation (Figures 6(c)–6(e)) the MAM predictors are surface temperature over the southwest Atlantic (−0.26), sea level air pressure over northern Europe (−0.35) and the south Indian Ocean (−0.50), and 200 mb zonal winds over the northeast Pacific (−0.60). Table 1(b) lists their domains and arrows again identify their position. One signal which stands out is negative Ts correlations in the Indian Ocean north of 20S (cf. Figure 6(a)) that affect early summer vegetation in Oromia. Another feature is the slow westward shift of 200 mb zonal wind signals in the eastern Pacific, indicative of poleward-displaced subtropical jets in both hemispheres (cf. Figures 6(b) and 6(e)) consistent with ENSO-PDO.
Correlation of ASO northern vegetation with MAM (a) Ts, (b) SLP, and (c) 200 U fields 1981–2006.
Correlation of MJJ southern vegetation with DJF (a) Ts and (b) 200 U fields. Correlation of ASO southern vegetation with MAM (c) Ts, (d) SLP, and (e) 200 U fields 1981–2006.
3.2. Algorithm Performance
Scatterplots of the multivariate algorithms are given in Figure 7 for Rain and Figure 8 for Vegetation, and Tables 1(a) and 1(b) lists the r2 fit. For northern rainfall the predictive algorithm achieves r2 of 0.51 in early summer and 0.53 in late summer, a satisfactory result. The multivariate algorithm for southern rainfall exhibits similar values of 0.49 and 0.51 for MJJ and ASO seasons, respectively. Curiously ARC2 and GPCC6 rainfall diverge over the southern area in late summer, with the later exhibiting higher amounts (cf. Figure 1(b)) and greater predictability. The scatterplots of predicted and observed rainfall show dispersion within the normal tercile and tend to taper favourably toward the extremes. Outliers are noted more in the south, wherein the algorithm yields neutral forecasts of early summer rains in 1996 (obs = 2.1σ) and in 2003 (obs = −1.2σ). The algorithm for southern rainfall underpredicts late summer rains in 1988 (obs = 2.4σ). The results suggest that statistical predictions in March for MJJ rainfall and in June for ASO rainfall should be reliable.
Scatterplots of predicted and observed rainfall according to Table 1(a). Southern ASO represented by GPCC6 rain; others are ARC2. Linear regression lines are applied.
Scatterplots of predicted and observed NDVI vegetation according to Table 1(b), with northern MJJ season omitted.
The multivariate algorithms for vegetation exhibit low skill in early summer (r2=0.43 south) possibly due to the erratic nature of warm spells. A substantially better fit is achieved in late summer when the multivariate algorithm for northern vegetation reaches 0.67, while southern vegetation r2 is 0.62. As NDVI follows crop production [23] and grazing yields, its predictability is most fortunate. The scatterplots for vegetation forecasts display correct tercile “hits” in late summer particularly in the northern zone, with drought underpredicted in 1991 (obs = −1.9σ). The algorithm for late summer southern vegetation exhibits three neutral forecasts when the corresponding observations were in the upper or lower tercile. One outlier (false alarm) is noted for early summer southern vegetation, when drought was predicted in 2005 but neutral conditions ensued. The r2 fit of our statistical model forecasts is almost triple those reported by Wang and Fan [16] for northern hemisphere summer rainfall.
4. Concluding Discussion
In this study, statistical algorithms to predict rainfall and vegetation over Amhara and Oromia, Ethiopia were formulated. The target rainfall derive from high quality ARC2 and GPCC6 reanalysis that diverge slightly in the southern zone in late summer (cf. Figure 1(b)). Correlation maps were made with respect to Ts, SLP, and 200/850 mb zonal winds which identified predictor-target relationships. From Table 1, Ts over Eurasia was the most influential predictor of rainfall (cold-favourable), followed by 850 mb zonal winds over the Atlantic (easterly-favourable), Ts over the tropical Indian Ocean (cold-favourable), and SLP over the south Indian Ocean (low-favourable). Predictors were selected evenly from the Atlantic, Indian, and Pacific basins, but ENSO-PDO signals were less influential (cf. Figures 4(b) and 6(e); w.r.t. southern vegetation). The rain targets selected zonal wind eight times, while the vegetation targets selected air pressure five times. After step-wise regression, algorithms for early and late season rainfall exhibited a consistent fit of ~0.50, while those for vegetation reached ~0.65 in late summer. Tercile hit rates exceeded two-thirds and false alarms were below one-third. How can our results be applied?
Ethiopia currently experiences food deficits because of low resource inputs, high population density, and variable climate. Shortfalls can be overcome with scientific information and practical engineering solutions [8]. Our climate predictions will be extended using crop yield simulators and formulated into coherent advice to collective decision-makers/resource managers and individual investors/farmers. Mitigating actions will be suggested by agricultural extension services, with feedback from their network of cultivators dealing with impacts.
We have sought to enhance our predictive capacity by statistically analyzing how surface temperature and pressure anomalies affect the overlying zonal circulation and tropical convection. The predictability of targeted resources is gratifying and highlights the ascendancy of extratropical signals (AO) and weakness of ENSO influence. Seasonal forecasts of early and late summer rainfall and vegetation over Ethiopia based on DJF and MAM predictors will provide the necessary lead time for strategic decisions to mitigate adverse impacts or take advantage of favourable weather. Processes underlying the apparent skill of extratropical predictors need to be tested via coupled ensemble model simulations. Further work is recommended by extending the predictor search to the subsurface oceans and meridional winds and by consideration of targets such as air temperature, evaporation, and streamflow.
Conflict of Interests
The author declares that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This study is part of a Rockefeller Foundation Project with the Ethiopian Institute for Agriculture Research, Melkassa supported by F. Getachew, A. Shimeles, and G. Mamo.
FAO2013Rome, ItalyFood and Agriculture OrganisationGlanzM. H.KatzR. W.NichollsN.1991Cambridge, UKCambridge University PressMouraA. D.SarachikE. S.Seasonal-to-interannual climate prediction and applications: new institutions, new possibilities19974643423472-s2.0-0031455220WCRP199929World Climate Research Programme, CLIVAR. WRCP Report 16/1999ICPO Publication SeriesMeinkeH.HochmanZ.HammerG. L.NichollsN.MitchellC.Using seasonal climate forecasts to manage dryland crops in northern Australia: experiences from the 1997/98 seasons2000Dordrecht, The NetherlandsKluwer Academic Publishers149160JuryM. R.Economic impacts of climate variability in South Africa and development of resource prediction models200241146552-s2.0-0036333791SadoffC.Can water resources undermine development in Ethiopia?20061814GoddardL.AitchelloucheY.BaethgenW.DettingerM.GrahamR.HaymanP.KadiM.MartínezR.MeinkeH.ConradE.Providing Seasonal-to-interannual climate information for risk management and decision-making1Procedia Environmental SciencesSeptember 2009811012-s2.0-7865084089410.1016/j.proenv.2010.09.007WhiteW. B.TourreY. M.Global SST/SLP waves during the 20th century200330125357 2-s2.0-0346783255JuryM. R.HuangB.The Rossby wave as a key mechanism of Indian Ocean climate variability20045112212321362-s2.0-884424156710.1016/j.dsr.2004.06.005YeshanewA.JuryM. R.North African climate variability. Part 2: tropical circulation systems2007891-237492-s2.0-3424721121410.1007/s00704-006-0243-7SegeleZ. T.LambP. J.LeslieL. M.Seasonal-to-interannual variability of Ethiopia/horn of Africa monsoon. Part I: associations of wavelet-filtered large-scale atmospheric circulation and global sea surface temperature20092212339634212-s2.0-6765121408210.1175/2008JCLI2859.1NdiayeO.WardM. N.ThiawW. M.Predictability of seasonal Sahel rainfall using GCMs and lead-time improvements through the use of a coupled model2011247193119492-s2.0-7995808167210.1175/2010JCLI3557.1OettliP.SultanB.BaronC.VracM.Are regional climate models relevant for crop yield prediction in West Africa?2011612-s2.0-7995296937010.1088/1748-9326/6/1/014008014008JuryM. R.Evaluation of coupled model forecasts of Ethiopian highlands summer climatein review, 2013WangH.FanK.A new scheme for improving the seasonal prediction of summer precipitation anomalies20092425485542-s2.0-6994912819910.1175/2008WAF2222171.1BrethertonC. S.SmithC.WallaceJ. M.An intercomparison of methods for finding coupled patterns in climate data19925541560BarnstonA. G.Linear statistical short-term climate predictive skill in the Northern Hemisphere1994715131564VautardR.PlautG.WangR.BrunetG.Seasonal prediction of North American surface air temperatures using space-time principal components19991223803942-s2.0-0033079860CamberlinP.JanicotS.PoccardI.Seasonality and atmospheric dynamics of the teleconnection between African rainfall and tropical sea-surface temperature: Atlantic vs. ENSO200121897310052-s2.0-003597314310.1002/joc.673KorechaD.BarnstonA. G.Predictability of June-September rainfall in Ethiopia200713526286502-s2.0-3384738809010.1175/MWR3304.1JuryM. R.MelesseA. M.Climatic factors modulating Nile River flow2011Springer267280JuryM. R.Ethiopian highlands crop-climate prediction 1979–20092013521116112610.1175/JAMC-D-12-0139.1NovellaN. S.ThiawW. M.African rainfall climatology version 2 for famine early warning systems20135258860610.1175/JAMC-D-11-0238.1SchneiderU.FuchsT.Meyer-ChristofferA.RudolfB.2008Global Precipitation Climatology Centre Rep.SchneiderU.BeckerA.FingerP.Meyer-ChristofferA.ZieseM.RudolfB.GPCC's new land surface precipitation climatology based on quality-controlled in situ data and its role in quanti-fying the global water cycle20141151-2154010.1007/s00704-013-0860-xTuckerC. J.PinzonJ. E.BrownM. E.SlaybackD.PakE. W.MahoneyR.VermoteE.El SaleousN.An extended AVHRR 8-km NDVI data set compatible with MODIS and SPOT vegetation NDVI data2005264485559810.1080/01431160500168686JuryM.RautenbachH.TadrossM.PhilippA.Evaluating spatial scales of climate variability in sub-Saharan Africa2007883-41691772-s2.0-3384738454010.1007/s00704-006-0251-7JuryM. R.Ethiopian decadal climate variability2010101129402-s2.0-7795380926310.1007/s00704-009-0200-3SmithT. M.ReynoldsR. W.PetersonT. C.LawrimoreJ.Improvements to NOAA's historical merged land-ocean surface temperature analysis (1880–2006)20082110228322962-s2.0-4534910900810.1175/2007JCLI2100.1AllanR.AnsellT.A new globally complete monthly historical gridded mean sea level pressure dataset (HadSLP2): 1850–200420061922581658422-s2.0-3374595945910.1175/JCLI3937.1SahaS.MoorthiS.PanH.-L.The NCEP climate forecast system reanalysis2010911015105710.1175/2010BAMS3001.1RieneckerM. M.SuarezM. J.GelaroR.MERRA: NASA's modern-era retrospective analysis for research and applications2011243624364810.1175/JCLI-D-11-00015.1ChangP.JiL.LiH.A decadal climate variation in the tropical Atlantic Ocean from thermodynamic air-sea interactions199738566165165182-s2.0-003096326210.1038/385516a0CohenJ.EntekhabiD.Eurasian snow cover variability and Northern Hemisphere climate predictability19992633453482-s2.0-0033082738GongG.EntekhabiD.CohenJ.Orographic constraints on a modeled Siberian snow-tropospheric-stratospheric teleconnection pathway2004176117611892-s2.0-1842739526AllenR. J.ZenderC. S.Forcing of the arctic oscillation by eurasian snow cover20112424652865392-s2.0-8415518529710.1175/2011JCLI4157.1