Printer Friendly

Evaluation of coupled model forecasts of Ethiopian highlands summer climate.

1. Introduction

Agricultural production is typically planned around a range of climatic conditions that take into account the possibility of flood or drought every decade. Commercial farmers have access to technology and finance, while subsistence farmers get by on local resources [1-4]. According to current FAO statistics, 76% of Ethiopia's 88 million people are engaged in farming on 15% of the land. There are orographic rains and a rich vegetation cover (Figure 1(a)), but a population density > 100 people/[km.sup.2] and crop yields <2 T/ha put a strain food supplies. Drought, for example in 1992-1993 and 20022003, caused malnutrition that required state assistance to 710 million people [5].

The climate observing network across northeast Africa, although adequate (Figure 1(b)), is declining from civil instability and financial constraints. While satellites supplement the atmospheric observations, seasonal forecast models need measurements in the upper ocean. These have been enhanced through the Global Ocean Observing System [6, 7] by extension of buoys and drifters near Africa ([8], Figure 1(c)). New surface flux data for coupled modelling have become available via satellite estimation of transpiration, soil water, and marine winds.

The predictability of climate is partly attributable to the Pacific El Nino Southern Oscillation (ENSO) and its overlying zonal circulation [9]. During ENSO warm phase, convection over most of Africa is suppressed [10, 11] and warm water spreads westward across the Indian Ocean [1214]. The Atlantic Ocean has a meridional sea temperature dipole that affects West African convection [15]. The three tropical ocean basins each have distinct rhythms that may interfere or reinforce each other [16]. In addition to large scale ocean-atmosphere coupling, there are sources of local forcing that can alter the incoming climate signals. These derive from rapid feedback between land surface fluxes and the lower atmosphere [17]. Decadal rainfall cycles have been found in Ethiopia [18] driven by slow eastward moving air pressure waves [16] that afford opportunities for seasonal forecasts even during ENSO transitions. The Ethiopian National Meteorological Agency (NMA) and the regional Climate Outlook Forum (GHACOF) have made seasonal forecasts since the 1990s using statistical methods [19]. These can be aimed at integrating targets such as crop yield [20] or river flow but have the disadvantage of assuming historical replication. More recently coupled ensemble models have been tested and used in seasonal forecasts [21].

In addition to the zonal circulation, the Hadley cell induces a large annual cycle over Ethiopia: soil moisture is depleted after January and replenished after June, when the equatorial trough reaches its northernmost limit. Given that highland crops are planted in June and harvested in October, users need seasonal forecasts during the warm spring season (March-April). This study evaluates coupled ensemble model forecasts (as in [22, 23]), considers factors driving Ethiopian climate fluctuations, and discusses potential applications.

2. Data and Methods

Retrospective forecasts from coupled ensemble models are compared with June to August (JJA) observations in the period 1981-2006 using data from the Climate Explorer website ( The models include CFSv2 (Coupled Forecast System, [24]) and ECMWFv3 (European Community Medium-range Weather Forecast, [25, 26]). Three observed products serve as reference: NMA interpolated rain gauge [27], GPCCv5 rainfall (Global Precipitation Climatology Center, [28, 29]), and CRUv3 (Climate Research Unit, [30]) maximum temperature, based on monthly 50 km gridded station data. Supplementary observations include low-resolution satellite-interpolated rainfall GPCPv2 (Global Precipitation Climatology Project, [31, 32]) and high-resolution reanalyses from CFS [33] and ECMWF [34]. Intercomparison of highlands area-average (7-14N, 3640E)NMA interpolated rainfall yields GPCC r = 0.92,GPCP 0.92, CRU3 0.85, CFS 0.71, and ECMWF 0.54, for continuous monthly data in the period 1981-2006. To form the observed highlands JJA rainfall dataset, NMA, and GPCC anomalies were averaged.

To evaluate the coupled ensemble models, correlations between forecast and observed JJA seasonal rainfall (R) and maximum temperature (Tx) anomalies were analyzed at various lead times from January to June and for the first and second half of sample (breakpoint 1993). Correlation maps were calculated to determine the spatial patterns of model performance, while temporal trends were evaluated using averages over the highlands 7-14N and 36-40E. This target is big enough to capture climate signals at long lead time yet falls within a singular climate regime [35]. In the seasonal model evaluation, statistical significance (at 90% confidence) is reached with r > 0.32 for 26 degrees of freedom.

Scatterplots of forecast and observed JJA anomalies were evaluated for slope, range, and outliers. Using anomalies helps to offset the mean bias and is consistent with operational forecasts of departures. Global signals driving local climate fluctuations were studied by correlation of NCEP [36] zonal winds, temperature, humidity, and vertical motion in a north-south vertical slice over the highlands. Satellite vegetation (NDVI) data [37] were analyzed for standard deviation, as a measure of agricultural vulnerability to annual and year-to-year fluctuations.

3. Results

3.1. JJA Pattern andAnnual Cycle. The 1981-2006 climatology of rainfall and maximum temperature as simulated and observed are illustrated in Figures 2(a) and 2(b). Both the CFS and the ECMWF models exhibit a cool wet bias compared with CRU3 and GPCC/GPCP observations in JJA season (e.g., Tx is ~2[degrees]C below observed, ~ 2 mm/day above). The model outputs could be "real" given that most observations are taken at urban airports located in warm dry valleys.

The ECMWF JJA rain pattern is shifted too far east, suggesting that orographic uplift on the escarpment at 40E is overplayed. The CFS JJA maximum temperature pattern is well-located but cool. The annual cycle of model outputs follows the R and Tx observations (Figures 2(a) and 2(b)). ECMWF has too much rain in JJA (consistent with [23]), while the CFS reflects an earlier onset of rains that conforms to the observed seasonal shape. Both models follow the annual cycle of maximum temperature with a cool offset. ECMWF and CFS Tx climatologies match in the rainy season, but ECMWF is warmer in the dry season and close to CRU3 observations then.

3.2. Spatial Performance. The correlation maps of CFS and ECMWF model outputs with respect to GPCC and CRU3 observations in the period 1981-2006 are given in Figure 3 for March and April forecasts. For CFS rainfall, the forecasts reach statistical significance in the northeast and southern highlands but not in the east and west escarpments. For ECMWF the rain forecasts are only significant in the northeast and actually negative in the southwest highlands particularly in March. Maximum temperature forecasts are weaker than rainfall in both models. The CFS obtains significant values in the central highlands in March which fade out in April, while ECMWF Tx forecasts are only significant in the southern highlands and weak or negative in the north. Forecasts are better for R than Tx, suggesting that model resolution (~1[degrees]) is not an issue.

3.3. Temporal Performance. Bar charts of highland area-averaged correlations between predicted and observed R and Tx for forecasts issued from January to June are given in Figure 4. CFS JJA rainfall forecasts are modest and comparable to ECMWF (Figure 4(a)) with r > 0.40 after April, consistent with Ndiaye et al. [21] for the Sahel. ECMWF rainfall forecasts (Figure 4(b)) dip in March and rise thereafter. CFS area-averaged JJA maximum temperature forecasts are weak and rise in March (Figure 4(d)), while ECMWF Tx forecasts start well in January and slump from March to April (Figure 4(e)). Forecasts show slight improvement over time (Figure 4(c)) mainly for CFS maximum temperature, possibly owing to improved satellite and ocean measurements (model technology is fixed). The steadiness of early season forecasts is assessed in Figure 4(f), wherein it is seen that January forecasts remain steady for CFS outputs of R and Tx. However ECMWF January forecasts slump in February and suggest better value after March. It is thought that the seasonal weakening of Pacific ENSO [38] and ambiguous coupling with Atlantic and Indian Oceans [39] are the cause of instability.

Scatterplots of area-averaged CFS March and ECMWF April forecasts versus observed JJA seasonal anomalies are illustrated in Figures 5(a)-5(d). It is evident that ECMWF forecasts have a narrower range for both R and Tx and thus tend toward the mean more than CFS. ECMWF rainfall forecasts exhibit a suitable 1: 1 slope with the 1984 drought as an outlier. CFS rainfall forecasts are well distributed and show highest fit (27%), but the flat 0.29 slope indicates over-prediction. Surprizingly, Tx forecasts are too dispersed in both models and consequently have insignificant fit (10%). Considering the outliers: CFS forecast Tx is too warm in 1986 while ECMWF forecast Tx is too warm in 1994. Both models under-predict 2002, anticipating neutral conditions instead of drought. The scatterplot equations-of-fit suggest potential adjustments to operational model outputs.

3.4. Climate Signals. Skillful forecasts depend on model ability to simulate ocean-atmosphere coupling and transmit global circulation anomalies to northeast Africa [9, 40-43]. For model simulations in the 1981-2006 period, inter-annual forcing deserves attention. Quasi-Biennial Oscillation (QBO, 30 mb tropical zonal wind) and ENSO (Pacific SST EOF1) indices are correlated with key variables in a north-south slice over the study area (Figures 6(a) and 6(b)). Both exhibit quite similar patterns. In west phase QBO and warm phase ENSO, easterly winds (r < 0) accelerate below 500 mb, driving away Congo moisture. Sinking motions warm the lower atmosphere, while westerly winds above 300 mb (r > 0) shear the convection and inhibit Indian monsoon outflow. Given the adequate performance by CFS and ECMWF models noted above, it is likely that these signals are initialized and transmitted. Yet forecast skill is marginalized by confounding influences from the Atlantic and Indian Oceans, and an opposing ENSO response in southern Ethiopia where the equatorial trough shifts rapidly in spring [44]. It is beyond the scope of this paper to evaluate model diagnostics.

3.5. Vulnerability. Considering the amplitude of vegetation (NDVI) response to climate impacts in the period 19812006, the standard deviation is calculated on monthly fields and departures (Figures 7(a) and 7(b)). The former identifies annual range, the latter inter-annual fluctuations. Annual range is greatest near the Sudan border west of Lake Tana, where maximum temperatures exceed 35 C (cf. Figure 2(b)). Annual range is low next to the large lakes and in the eastern lowlands where it is always warm. Interannual fluctuations are greatest in the southern highlands and along the eastern escarpment on 40E. Year-to-year changes of NDVI are low across the northern highlands (Tigray, Amhara). Hence vulnerability to climate in the southern highlands, eastern escarpment and western lowlands, makes the uptake of forecasts there critical to food security.

3.6. Application. The Ethiopian Institute for Agriculture Research (EIAR) uses seasonal forecasts from the NMA/ GHACOF and CFS/ECMWF modeling centers to develop an initial outlook for the season that guides farmers on how much area is planted and which hybrid seeds are used. The initial outlook is followed with bimonthly updates and advisories to optimize farming activities. The commercial sector carries bigger risks and is more responsive to technological inputs than the subsistence sector. As the season progresses, NDVI anomalies in cropped areas are analyzed (, cf. Figure 1(a)) and the EIAR obtains direct feedback from farm liaison officers. Agroclimate information networks are utilized to help the rural population avoid risks in dryyears and secure resources in wet years. Interventions are made at planting time, in the event of flood or drought and to collect harvest data. At the EIAR experimental farm in Melkassa, staff monitor crops and develop ways to improve yields. This is critical, because an upward trend of ~0.1 T [ha.sup.-1]/yr is needed to keep pace with Ethiopia's growing population.

4. Summary

This study has evaluated summer rainfall and maximum temperature forecasts by ECMWFv3 and CFSv2 models (cf. [22, 23]) via spatial correlation maps and area-averaged temporal analyses. Reference data were comprised of gridded NMA + GPCC R and CRU3 Tx observations in the Ethiopian highlands 7-14N 36-40E. Both models simulate the JJA mean spatial pattern with a 10% cool wet bias, and their March-April forecasts correlate positively with summer observations from 1981 to 2006. Considering that a cost benefit is possible with model "fit" above half [45], such skill is reached for CFS and ECMWF April forecasts of JJA rainfall over most of the highlands (cf. Figures 3, 4(a), and 4(b)). Yet the limit of predictability is evident in Figure 5, 25-27% of variance for rainfall and 9-10% for maximum temperature. Further work is recommended to understand causes of instability in early season forecasts, determine why maximum temperature forecasts are weak, employ more robust reference data [46], and develop bias corrections for improved model skill. At the EIAR experimental farm, numerical and statistical forecasts will be compared and utilized to develop mitigating strategies that boost crop yields.

Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.


This study is part of a Rockefeller Foundation project with the Ethiopian Institute for Agriculture Research, Melkasa.


[1] P. C. Stern and W. E. Easterling, Eds., Making Climate Forecasts Matter, National Academy Press, Washington, D.C., USA, 1999.

[2] R. Blench, "Seasonal climatic forecasting: who can use it and how should it be disseminated?" Natural Resource Perspectives 47, Overseas Development Institute, London, UK, 1999.

[3] J. W. Jones, J. W. Hansen, F S. Royce, and C. D. Messina, "Potential benefits of climate forecasting to agriculture," Agriculture, Ecosystems and Environment, vol. 82, no. 1-3, pp. 169-184, 2000.

[4] G. L. Hammer, J. W. Hansen, J. G. Phillips et al., "Advances in application of climate prediction in agriculture," Agricultural Systems, vol. 70, no. 2-3, pp. 515-553, 2001.

[5] A. S. Taffessee, "Decomposition of growth in cereal production in Ethiopia: the role of agriculture," DFID Report, Economics Assoc., Addis Ababa, Ethiopia, 2008.

[6] M. R. Jury and M. S. J. Harrison, "A strategic plan for climate research in Africa," WCRP Internal Report, 1999.

[7] D. Roemmich and W. Brechner Owens, "The Argo project: global ocean observations for understanding and prediction of climate variability," Oceanography, vol. 13, no. 2, pp. 45-50, 2000.

[8] M. J. McPhaden, G. Meyers, K. Ando et al., "RAMA: the research moored array for African-Asian-Australian monsoon analysis and prediction," Bulletin of the American Meteorological Society, vol. 90, no. 4, pp. 459-480, 2009.

[9] A. Yeshanew and M. R. Jury, "North African climate variability. Part 2. Tropical circulation systems," Theoretical and Applied Climatology, vol. 89, no. 1-2, pp. 37-49, 2007

[10] C. F Ropelewski and M. S. Halpert, "Quantifying southern oscillation-precipitation relationships," Journal of Climate, vol. 9, no. 5, pp. 1043-1059, 1996.

[11] S. J. Mason and L. Goddard, "Probabilistic precipitation anomalies associated with ENSO," Bulletin of the American Meteorological Society, vol. 82, no. 4, pp. 619-638, 2001.

[12] W. B. White, "Coupled Rossby Waves in the Indian Ocean on interannual timescales," Journal of Physical Oceanography, vol. 30, no. 11, pp. 2972-2988, 2000.

[13] M. R. Juryand B. Huang, "The Rossbywave asakeymechanism of Indian Ocean climate variability," Deep-Sea Research I: Oceanographic Research Papers, vol. 51, no. 12, pp. 2123-2136, 2004.

[14] T Yamagata, S. K. Behera, J.-J. Luo, S. Masson, M. R. Jury, and S. A. Rao, "Coupled ocean-atmosphere variability in the tropical Indian Ocean," in Earth Climate: The Ocean-Atmosphere Interaction, vol. 147 of Geophysical Monograph Series, pp. 189212, The American Geophysical Union, Washington, D.C., USA, 2004.

[15] P. Chang, T Yamagata, P. Schopf et al., "Climate fluctuations of tropical coupled systems: the role of ocean dynamics," Journal of Climate, vol. 19, no. 20, pp. 5122-5174, 2006.

[16] W. B. White and Y. M. Tourre, "Global SST/SLP waves during the 20th century," Geophysical Research Letters, vol. 30, no. 12, pp. 53-57, 2003.

[17] H. Chikoore and M. R. Jury, "Intraseasonal variability of satellite-derived rainfall and vegetation over Southern Africa," Earth Interactions, vol. 14, pp. 1-26, 2010.

[18] M. R. Jury, "Ethiopian decadal climate variability," Theoretical and Applied Climatology, vol. 101, no. 1, pp. 29-40, 2010.

[19] D. Korecha and A. G. Barnston, "Predictability of June-September rainfall in Ethiopia," Monthly Weather Review, vol. 135, no. 2, pp. 628-650, 2007.

[20] M. R. Jury, "Ethiopian highlands crop-Climate prediction: 1979-2009," Journal of Applied Meteorology and Climatology, vol. 52, no. 5, pp. 1116-1126, 2013.

[21] O. Ndiaye, M. N. Ward, and W. M. Thiaw, "Predictability of seasonal Sahel rainfall using GCMs and lead-time improvements through the use of a coupled model," Journal of Climate, vol. 24, no. 7, pp. 1931-1949, 2011.

[22] H.-M. Kim, P. J. Webster, and J. A. Curry, "Seasonal prediction skill of ECMWF System 4 and NCEP CFSv2 retrospective forecast for the Northern Hemisphere Winter," Climate Dynamics, vol. 39, no. 12, pp. 2957-2973, 2012.

[23] E. Dutra, F. di Giuseppe, F Wetterhall, and F. Pappenberger, "Seasonal forecasts of droughts in African basins using the standardized precipitation index," Hydrology and Earth System Sciences, vol. 17, no. 6, pp. 2359-2373, 2013.

[24] S. Saha, S. Moorthi, H.-L. Pan et al., "The NCEP climate forecast system reanalysis," Bulletin of the American Meteorological Society, vol. 91, pp. 1015-1057, 2010.

[25] J. Vialard, F. Vitart, M. A. Balmaseda, T. N. Stockdale, and D. L. T. Anderson, "An ensemble generation method for seasonal forecasting with an ocean-atmosphere coupled model," Monthly Weather Review, vol. 133, no. 2, pp. 441-453, 2005.

[26] T. N. Stockdale, D. L. T. Anderson, M. A. Balmaseda et al., "ECMWF seasonal forecast system 3 and its prediction of sea surface temperature," Climate Dynamics, vol. 37, no. 3, pp. 455471, 2011.

[27] M. M. Bitew and M. Gebremichael, "Assessment of satellite rainfall products for streamflow simulation in medium watersheds of the Ethiopian highlands," Hydrology and Earth System Sciences, vol. 15, no. 4, pp. 1147-1155, 2011.

[28] U. Schneider, T. A. Fuchs, A. Meyer-Christoffer, and B. Rudolf, "Global precipitation analysis products of the GPCC," Global Precipitation Climatology Centre Rep, 2008.

[29] U. Schneider, A. Becker, P. Finger, A. Meyer-Christoffer, M. Ziese, and B. Rudolf, "GPCC's new land surface precipitation climatology based on quality-controlled in situ data and its role in quantifying the global water cycle," Theoretical and Applied Climatology, vol. 115, no. 1-2, pp. 15-40, 2014.

[30] T. D. Mitchell and P. D. Jones, "An improved method of constructing a database of monthly climate observations and associated high-resolution grids," International Journal of Climatology, vol. 25, no. 6, pp. 693-712, 2005.

[31] R. F. Adler, G. J. Huffman, A. Chang et al., "The version-2 global precipitation climatology project (GPCP) monthly precipitation analysis (1979-present)," Journal of Hydrometeorology, vol. 4, no. 6, pp. 1147-1167, 2003.

[32] G. J. Huffman, R. F. Adler, D. T Bolvin, and G. Gu, "Improving the global precipitation record: GPCP Version 2.1," Geophysical Research Letters, vol. 36, no. 17, Article ID L17808, 2009.

[33] S. Saha, S. Nadiga, C. Thiaw et al., "The NCEP climate forecast system," Journal of Climate, vol. 19, no. 15, pp. 3483-3517, 2006.

[34] D. P. Dee, S. M. Uppala, A. J. Simmons et al., "The ERA-Interim reanalysis: configuration and performance of the data assimilation system," Quarterly Journal of the Royal Meteorological Society, vol. 137, no. 656, pp. 553-597, 2011.

[35] M. R. Jury, H. Rautenbach, M. Tadross, and A. Philipp, "Evaluating spatial scales of climate variability in sub-Saharan Africa," Theoretical and Applied Climatology, vol. 88, no. 3-4, pp. 169-177, 2007

[36] E. Kalnay, M. Kanamitsu, R. Kistler et al., "The NCEP/NCAR 40-year reanalysis project," Bulletin of the American Meteorological Society, vol. 77, no. 3, pp. 437-471, 1996.

[37] C. J. Tucker, J. E. Pinzon, M. E. Brown et al., "An extended AVHRR 8-km NDVI dataset compatible with MODIS and SPOT vegetation NDVI data," International Journal of Remote Sensing, vol. 26, no. 20, pp. 4485-4498, 2005.

[38] R. Wu, B. P Kirtman, and H. van den Dool, "An analysis of ENSO prediction skill in the CFS retrospective forecasts," Journal of Climate, vol. 22, no. 7, pp. 1801-1818, 2009.

[39] C.-P Chang, P Harr, and J. Ju, "Possible roles of atlantic circulation on the weakening Indian monsoon rainfall-ENSO relationship," Journal of Climate, vol. 14, no. 11, pp. 2376-2380, 2001.

[40] D. Conway, M. Krol, J. Alcamo, and M. Hulme, "Future availability of water in Egypt: the interaction of global, regional, and basin scale driving forces in the Nile Basin," Ambio, vol. 25, no. 5, pp. 336-342, 1996.

[41] P Camberlin, "Rainfall anomalies in the source region of the Nile and their connection with the Indian summer monsoon," Journal of Climate, vol. 10, no. 6, pp. 1380-1392, 1997

[42] R. Washington and T. E. Downing, "Seasonal forecasting of African rainfall: prediction, responses and household food security," Geographical Journal, vol. 165, no. 3, pp. 255-274, 1999.

[43] Z. T. Segele, P J. Lamb, and L. M. Leslie, "Seasonal-to-interannual variability of Ethiopia/horn of Africa monsoon. Part I: associations of wavelet-filtered large-scale atmospheric circulation and global sea surface temperature," Journal of Climate, vol. 22, no. 12, pp. 3396-3421, 2009.

[44] E. E. Riddle and K. H. Cook, "Abrupt rainfall transitions over the Greater Horn of Africa: observations and regional model simulations," Journal of Geophysical Research D: Atmospheres, vol. 113, no. 15, Article ID D15109, 2008.

[45] A. H. Murphy, "Forecast verification," in The Economic Value of Weather and Climate Forecasts, R. W. Katz and A. H. Murphy, Eds., Cambridge University Press, Cambridge, Mass, USA, 1997.

[46] G. M. Tsidu, "High-resolution monthly rainfall database for Ethiopia: homogenization, reconstruction, and gridding," Journal of Climate, vol. 25, no. 24, pp. 8422-8443, 2012.

Mark R. Jury (1,2)

(1) University of Zululand, KwaDlangezwa 3886, South Africa

(2) Physics Department, University of Puerto Rico, Mayaguez, PR 00681, USA

Correspondence should be addressed to Mark R. Jury;

Received 1 April 2014; Revised 19 August 2014; Accepted 18 September 2014; Published 14 October 2014

Academic Editor: George Kallos
COPYRIGHT 2014 Hindawi Limited
No portion of this article can be reproduced without the express written permission from the copyright holder.
Copyright 2014 Gale, Cengage Learning. All rights reserved.

Article Details
Printer friendly Cite/link Email Feedback
Title Annotation:Research Article
Author:Jury, Mark R.
Publication:Advances in Meteorology
Article Type:Report
Date:Jan 1, 2014
Previous Article:Variations of near surface energy balance caused by land cover changes in the semiarid grassland area of China.
Next Article:Dynamical downscaling of climate change impacts on wind energy resources in the contiguous united states by using a limited-area model with...

Terms of use | Privacy policy | Copyright © 2021 Farlex, Inc. | Feedback | For webmasters |