details of all known meningitis epidemics occurring before 2000 in countries comprising continental africa were compiled from information documented in the published literature and unpublished institutional reports at the end of june 2001 . all epidemics reported in the medical literature were identified in pubmed s online database of medical literature ( united states national library of medicine , available from : url : http://www.ncbi.nlm.nih.gov/pubmed/ ) through manual searches and by cross - referencing publications . we obtained unpublished information from web searches or directly from international and national organizations involved in disease control and humanitarian aid ( 1 ) . epidemics reported at the provincial ( second administrative level ) and district level were located by using administrative boundaries available from the u.s . geological survey eros data centre africa data dissemination service , sioux falls , sd ( available from : url : http://edcintl.cr.usgs.gov/adds/ ) . locations of villages or towns were verified by using the national imagery and mapping agency s geonet gazetteer ( available from : url : http://www.nima.mil/ ) , and historical reference atlases and maps contained in the original reports were mapped to the current administrative boundaries by using arcview 3.1 geographic information system ( gis ) ( esri , redlands , ca ) . we assumed that events reported at the provincial level affected all constituent districts and excluded epidemics reported only at the national level . events affecting parts of any one country in one epidemic season were considered as one epidemic . we then classified the 3,281 districts of africa as ever or never having experienced a documented epidemic of meningitis . in this analysis , environmental information for the african continent was obtained from a variety of sources . in this analysis , we restricted data to variables available in the public domain with digital grid  based uniform continental coverage , which we considered important ( 7 ) . variables included monthly means ( long - term averages ) of absolute humidity , absorbing aerosols ( dust ) and rainfall , and land - surface maps of land - cover type and population density ( table 1 ) . mean monthly climate averages ( 19611990 ) for absolute humidity and rainfall ( the former derived from vapor pressure and mean temperature ) were obtained from the climate research unit , university of east anglia , uk ( available from : url : http://ipcc-ddc.cru.uea.ac.uk/ ) . national aeronautics and space administration ( nasa ) , goddard space flight center , maryland ( available from : url : http://www.gsfc.nasa.gov/ ) , excluding the period may 1993 to june 1996 for which data are not available . land - cover type was obtained from the usgs nasa pathfinder 1 km project ( available from : url : http://edcwww.cr.usgs/gov/landdaac/1km/ ) . u.s . geologic survey 1990 population density forecasts ( available from : url : http://grid2.cr.usgs.gov/globalpop/africa/ ) . to collect the seasonal variation in climate while reducing the number of explanatory variables without loss of information , we reprocessed the monthly means of each variable to create a single surface comprising categories representing unique seasonal profiles . this reprocessing involved submitting the monthly means to a principal components analysis followed by a clustering procedure in which we grouped regions with a similar seasonal profile using addapix software ( version 2.05 , s. griguolo , university of venice , italy , available from : url : http://metart.fao.org/t_i/gbr/tools/eaddapix.htm ) . this software performs spatial and temporal analyses of time - series data in continuous grid - based surfaces ( 8) . the use of seasonal profiles for describing the climate of an area is widely used in crop monitoring in agriculture and has been used in modeling malaria prevalence in gambia ( 9 ) . a ) spatial variation in profile class ; b ) representative profile class . for every district , a value was extracted in gis , representing those grid cells of each profile surface contained within the district boundary . these values comprised the most common seasonal profile class for each variable ( absolute humidity , dust , and rainfall ) , the most common land - cover type , and the geometric mean population density . a logistic regression analysis was used to identify associations between a district ever or never having experienced an epidemic and the environment by using spss 11.0 software ( spss inc . , explanatory variables were first investigated individually and then entered stepwise into a forward conditional multiple regression analysis . because of the need to reduce colinearity , we did not analyze environmental variables with similar geographic distributions together in the same model . the final model was based on the simplest approach and a combination of variables that best predicted the distribution of epidemics . this model was created with and without weighting for the inverse size of the district to assign less weight to larger districts , which may have been more prone to ecologic variation and therefore inadequately represented by a single value . the probability of each district ever having had an epidemic was predicted by using the model . these probabilities were grouped into risk categories and mapped , and the estimated total population was extracted in gis derived from population density forecasts ( u.s . geological survey , 1990 ) . the sensitivity and specificity of the model were assessed by examining the agreement between predicted and observed epidemic experience by using a receiver - operator characteristics ( roc ) curve to select the optimal probability cutoff values on which predictions are based ( 10 ) . the dataset was then split at random into two parts containing approximately 60% and 40% of the 3,281 districts . the model was recreated with the 60% dataset , using the same variables as above and used to predict the risk for epidemics in the remaining 40% ( the validation set ) . we repeated this process 10 times and compared the mean sensitivity and specificity of the validation set with the model derived from the entire dataset . residuals resulting from the differences between the observed and predicted risks were also calculated and mapped to establish whether errors in the model were randomly distributed , thereby supporting the validity of the model throughout the continent . details of all known meningitis epidemics occurring before 2000 in countries comprising continental africa were compiled from information documented in the published literature and unpublished institutional reports at the end of june 2001 . all epidemics reported in the medical literature were identified in pubmed s online database of medical literature ( united states national library of medicine , available from : url : http://www.ncbi.nlm.nih.gov/pubmed/ ) through manual searches and by cross - referencing publications . we obtained unpublished information from web searches or directly from international and national organizations involved in disease control and humanitarian aid ( 1 ) . epidemics reported at the provincial ( second administrative level ) and district level were located by using administrative boundaries available from the u.s . geological survey eros data centre africa data dissemination service , sioux falls , sd ( available from : url : http://edcintl.cr.usgs.gov/adds/ ) . locations of villages or towns were verified by using the national imagery and mapping agency s geonet gazetteer ( available from : url : http://www.nima.mil/ ) , and historical reference atlases and maps contained in the original reports were mapped to the current administrative boundaries by using arcview 3.1 geographic information system ( gis ) ( esri , redlands , ca ) . we assumed that events reported at the provincial level affected all constituent districts and excluded epidemics reported only at the national level . events affecting parts of any one country in one epidemic season were considered as one epidemic . we then classified the 3,281 districts of africa as ever or never having experienced a documented epidemic of meningitis . in this analysis , environmental information for the african continent was obtained from a variety of sources . in this analysis , we restricted data to variables available in the public domain with digital grid  based uniform continental coverage , which we considered important ( 7 ) . variables included monthly means ( long - term averages ) of absolute humidity , absorbing aerosols ( dust ) and rainfall , and land - surface maps of land - cover type and population density ( table 1 ) . all data grids were incorporated into the gis . mean monthly climate averages ( 19611990 ) for absolute humidity and rainfall ( the former derived from vapor pressure and mean temperature ) were obtained from the climate research unit , university of east anglia , uk ( available from : url : http://ipcc-ddc.cru.uea.ac.uk/ ) . national aeronautics and space administration ( nasa ) , goddard space flight center , maryland ( available from : url : http://www.gsfc.nasa.gov/ ) , excluding the period may 1993 to june 1996 for which data are not available . land - cover type was obtained from the usgs nasa pathfinder 1 km project ( available from : url : http://edcwww.cr.usgs/gov/landdaac/1km/ ) . u.s . geologic survey 1990 population density forecasts ( available from : url : http://grid2.cr.usgs.gov/globalpop/africa/ ) . to collect the seasonal variation in climate while reducing the number of explanatory variables without loss of information , we reprocessed the monthly means of each variable to create a single surface comprising categories representing unique seasonal profiles . this reprocessing involved submitting the monthly means to a principal components analysis followed by a clustering procedure in which we grouped regions with a similar seasonal profile using addapix software ( version 2.05 , s. griguolo , university of venice , italy , available from : url : http://metart.fao.org/t_i/gbr/tools/eaddapix.htm ) . this software performs spatial and temporal analyses of time - series data in continuous grid - based surfaces ( 8) . the use of seasonal profiles for describing the climate of an area is widely used in crop monitoring in agriculture and has been used in modeling malaria prevalence in gambia ( 9 ) . a ) spatial variation in profile class ; b ) representative profile class . for every district , a value was extracted in gis , representing those grid cells of each profile surface contained within the district boundary . these values comprised the most common seasonal profile class for each variable ( absolute humidity , dust , and rainfall ) , the most common land - cover type , and the geometric mean population density . a logistic regression analysis was used to identify associations between a district ever or never having experienced an epidemic and the environment by using spss 11.0 software ( spss inc . , explanatory variables were first investigated individually and then entered stepwise into a forward conditional multiple regression analysis . because of the need to reduce colinearity , we did not analyze environmental variables with similar geographic distributions together in the same model . the final model was based on the simplest approach and a combination of variables that best predicted the distribution of epidemics . this model was created with and without weighting for the inverse size of the district to assign less weight to larger districts , which may have been more prone to ecologic variation and therefore inadequately represented by a single value . the probability of each district ever having had an epidemic was predicted by using the model . these probabilities were grouped into risk categories and mapped , and the estimated total population was extracted in gis derived from population density forecasts ( u.s . geological survey , 1990 ) . the sensitivity and specificity of the model were assessed by examining the agreement between predicted and observed epidemic experience by using a receiver - operator characteristics ( roc ) curve to select the optimal probability cutoff values on which predictions are based ( 10 ) . the dataset was then split at random into two parts containing approximately 60% and 40% of the 3,281 districts . the model was recreated with the 60% dataset , using the same variables as above and used to predict the risk for epidemics in the remaining 40% ( the validation set ) . we repeated this process 10 times and compared the mean sensitivity and specificity of the validation set with the model derived from the entire dataset . residuals resulting from the differences between the observed and predicted risks were also calculated and mapped to establish whether errors in the model were randomly distributed , thereby supporting the validity of the model throughout the continent . the earliest documented meningitis outbreak in africa affected a french garrison in 1841 in algiers ( 11 ) , and > 425 epidemics were documented at the subnational level for the next 158 years . these epidemics affected at least 1,231 ( 38% ) of the 3,281 continental districts . supporting our earlier findings , epidemics were not evenly distributed across the continent , instead affecting mainly districts in the sahel and south of this region and extending from northern uganda and the eastern part of democratic republic of congo , through the great lakes and the rift valley regions to malawi and northeastern mozambique , and from northeastern mozambique west and south to include other parts of southern africa ( figure 2a ) ( 1 ) . absolute humidity , dust and rainfall profiles , land - cover type , and population densities were independently associated with the location of epidemics . however , we found that absolute humidity profile and land - cover type were the best predictors in the final multivariate model . certain dust profiles and population density made only marginal difference to the performance of the model and were excluded to maintain simplicity . in addition , absolute humidity and rainfall profiles predicted similar risks in geographic locations in north and west africa , but a model including the former performed better for the entire continent than one including rainfall ; for this reason , we kept absolute humidity as a variable . weighting for district size did not improve the performance of the model and was discarded . the model , which is based on the absolute humidity profile and land - cover type , is described in terms of its baseline characteristics , estimated coefficients , standard errors , and contribution of the variables in table 2 . the risk map for epidemic experience in africa derived from this model is presented in figure 2b . the most important factor associated with areas without a marked distinction between wet and dry seasons were less likely to have had epidemics than those with contrasting seasons . the areas without distinction between wet and dry seasons include deserts and the humid parts of coastal and central africa , much of which are forested , and the areas with contrasting seasons comprise the semiarid savannah and grasslands of the sahel and east and southern africa . surface maps of africa demonstrated a close correspondence between humidity and land - cover types in these regions . the model also showed that , having accounted for the effects of humidity , sparsely vegetated and barren regions , areas of woodland mosaic , and shrub land were less likely than other regions to have ever had an epidemic . the sahel , which has a prolonged dry season with low humidity was identified as the area with the greatest risk ( p>0.6 ) . peripheral regions along its southern borders , where the dry season is shorter and less extreme , carry a moderate risk ( p>0.4 ) . the peripheral region extends from southern sudan and ethiopia to the great lakes and rift valley regions and parts of southern africa peripheral to desert areas . excludes five districts for which environmental data were unavailable ( in sinai , egypt , and pemba , tanzania ) . the roc curve used to describe the performance of the model in terms of its sensitivity and specificity at various cutoffs ; its overall accuracy independent of a single probability cutoff is shown in figure 3 . the model can discriminate between districts that have experienced epidemics and those that have never been affected . when the values for randomly selected districts were entered into the model , the epidemic risk assigned was higher for districts with epidemics than for those without in 82% of the cases . the areas identified by the model coincide to a large extent with the areas observed to have experienced epidemics . by using a probability cutoff value of > 0.4 for predicting epidemic experience , the model had a sensitivity and specificity of 83% and 67% , respectively ; these statistics were confirmed in the validation process ( table 3 ) . the map of residual risk unaccounted for by the model had a random distribution , as expected for a model that worked well across the continent ( not shown ) . according to the model , 7% of the population of africa live in very high - risk areas , 17% in high , and 27% in moderate ( based on 1990 estimates ) , corresponding to 44 , 102 , and 162 million people , respectively ( table 4 ) . this study represents a comprehensive and detailed description of the spatial distribution of meningitis epidemics at the district level in africa and the first attempt to develop a spatial forecasting model for meningitis epidemics on the basis of the environmental characteristics of the continent considered a priori to be related to the spatial distribution of epidemics . the data have limitations that need to be considered for the proper interpretation of the models . for example , census data for africa have limited accuracy ( 12 ) , variables derived through remote sensing may only partially capture surface conditions , and data from meteorologic stations in africa are often incomplete ( 1315 ) . the epidemiologic data span more than a century whereas the environmental and population data are relatively recent . in addition , some epidemics were likely never reported , and small outbreaks and those occurring in the 19th and early 20th centuries are likely to be disproportionately underrepresented . problems in defining epidemics exist as well , since most reports lack clear or internationally recognized criteria ; we had to accept the perception of an increased incidence that prompted outbreak reports . moreover , the aggregation of local level statistics to an often larger and somewhat arbitrary district level , discrepancies between where people become ill and the location of notifying health facilities , and changes to district boundaries over time may each have resulted in potential loss of specificity . despite these limitations , major outbreaks were unlikely to have gone completely unreported , and the long - term cumulative distribution of events is likely not misrepresented on a pancontinental scale . while population densities in africa have increased greatly during the last 150 years and substantial land - use change ( particularly in west africa ) is known to have occurred , the model was still able to identify the meningitis belt and areas previously described at risk beyond the sahel ( 25,16 ) ; reports of epidemics occurring since 1999 have coincided with this description ( available from : url : http://www.who.int/disease-outbreak-news/ ) . we are aware of those imposed on the analysis by spatial relationships in the data , which indicate that the importance of regression effects may be overstated ( 17 ) . no perfect model exists , but one can be developed that is based on the simplest approach and combination of variables that can be used to distinguish between areas with high and low risk of epidemic experience . we were restricted by the limited availability of suitable datasets and chose the simplest combination of variables that best predicted outcome . while this choice is likely to have oversimplified the association between meningitis epidemics and the environment , the model has the advantage of using data that are available in the public domain and , being based on a simple combination of variables , are an important basis on which to develop research and future operational applications in resource - limited settings . the model indicates not only that absolute humidity profiles and land - cover types can be used to distinguish between areas with high and low risk of epidemics but also that population density and dust may also be implicated . the incidence of meningococcal disease has previously been correlated with dry and dusty conditions in tropical and temperate climates ( 1821 ) . humidity and land cover were included in the final model for statistical reasons , but dust and population density still have an independent effect and may be important in determining epidemic occurrence ( 22 ) . the potential role of dust in precipitating epidemics is particularly interesting since dustiness in the meningitis belt has increased dramatically since the sahelian droughts of the 1970s and 1980s . however , how environmental variables interact is unclear and remains the subject of extensive climatologic research . furthermore , we did not take into account the effect of other nonenvironmental factors likely to be related to epidemics , such as population movement , vaccination coverage , and recent epidemics in the area . a combination of conditions is likely to be necessary for an epidemic to occur , and these nonenvironmental variables are likely to have additional predictive potential and should be considered in further studies . risk maps of vector - borne diseases in africa based on environmental data have received considerable attention in recent years and are tools with public health potential ( 10,15 ) . the mechanisms by which environmental factors influence meningitis epidemics in africa are unclear ( 3 ) . areas within the traditional meningitis belt and beyond , however , are environmentally susceptible to epidemics with potentially large populations at risk ; markers such as absolute humidity and dust profiles , land - cover type , and population density are independent predictors of these areas . in addition , rainfall and dust are predictors in some , but not all regions , and the potential to develop region - specific models that could be more sensitive within given ecologic zones warrants further study . our findings could facilitate the development of models to identify regions with increased vulnerability to epidemics in the future and provide a basis for monitoring the impact of climate variability and environmental change on epidemic occurrence in africa .
epidemics of meningococcal meningitis occur in areas with particular environmental characteristics . 
 we present evidence that the relationship between the environment and the location of these epidemics is quantifiable and propose a model based on environmental variables to identify regions at risk for meningitis epidemics . 
 these findings , which have substantial implications for directing surveillance activities and health policy , provide a basis for monitoring the impact of climate variability and environmental change on epidemic occurrence in africa .