despite its importance , current knowledge on the nature and drivers of changing endemicity in sub - saharan africa is remarkably weak . national health records in 32 highly endemic countries ( together accounting for about 90% of the global malaria burden ) are considered inadequate to assess trends in malaria cases . this stems from low care - seeking rates ( many malaria cases are not seen at formal health facilities ) , incomplete record keeping and curation ( many recorded cases are never captured in surveillance databases ) , and historically poor access to parasitological diagnosis ( malaria cases were often diagnosed presumptively with poor specificity ) . as systems have begun to improve , for example owing to greater use of rapid diagnostic testing in health facilities , these biases have been mitigated , but that presents further challenges for comparison of data through time and evaluation of trends . for these countries , the who has previously adopted a cartographic  burden measurement approach whereby a map of climatic suitability for malaria transmission is first used to stratify likely incidence rates across the continent , with these rates then progressively downgraded as intervention coverage increases according to effect sizes measured in randomized control trials . one acknowledged limitation of this approach is its reliance on the central assumption that effects observed in a limited number of short - term trials can be extrapolated to sustained continent - wide implementation . this assumption has never been validated beyond local or national - level analyses . in reality , the individual or combined efficacy of interventions will vary by setting and be contingent on many local factors including vector ecology , health systems , and coverage levels . these studies capture a wider range of real - world settings , but yield only a single pooled estimate across diverse disease transmission settings that , again , has unknown validity when extrapolated across africa . since its first use , component parts of the cartographic burden framework have incrementally improved . climatic suitability maps have been superseded by empirical endemicity maps that use model - based geostatistics to create surfaces of risk based on thousands of geolocated cross - sectional surveys measuring infection prevalence ( termed plasmodium falciparum parasite rate , pfpr ) . improvements have also been made in the estimation of clinical incidence rates as a function of pfpr , allowing clinical incidence rates , which are notoriously difficult to measure in the field , to be estimated geographically using mapped surfaces of pfpr . all these earlier studies , however , preceded the most intense period of control effort ( from 2010 to the present ) , and none were designed to formally evaluate temporal changes in disease burden or explicitly consider the effect of interventions . here , we provide the first formal quantification , with rigorously defined uncertainty , of p. falciparum infection prevalence and disease incidence across sub - saharan africa from the year 2000 to the benchmark year of 2015 , and of the role the major control interventions have had in causing these changes . our approach evaluates not only point estimates , but also presents a full treatment of uncertainty through a bayesian hierarchical model . components contributing to uncertainty in outputs included the sample size and spatiotemporal density of pfpr surveys , uncertainty in the fitted relationships between pfpr and the suite of environmental and intervention covariates , uncertainty in the input data on observed clinical incidence rates , and uncertainty in the mechanistic model parameters defining the prevalence - incidence relationship . by linking all components together in a bayesian framework , these distinct sources of uncertainty are formally propagated through the predictive model and represented as predictive posterior distributions around all output results . data on itn use and access to acts from over one million households were combined with national malaria control programme data on itn , act , and irs provision to develop time - series models of coverage of these interventions within each country . these were combined within a spatiotemporal bayesian geostatistical model with pfpr data from 27,573 georeferenced population clusters between 1995 and 2014 , along with an optimised suite of temporally dynamic environmental and sociodemographic covariates . the model adjusted pfpr observations by age , season , and type of diagnostic used , and fitted flexible functional forms to capture the effect of each intervention on declining pfpr as a function of coverage reached and the starting ( pre - intervention ) pfpr in 2000 ( extended data fig . 2 ) . following earlier work we chose to model pfpr in the 2-up - to-10 year age range , since this is associated with a plateau in the age - prevalence relationship and thus acts as a standardised comparison . the model was used to predict a spatio - temporal  cube  of age - structured pfpr at 55 km resolution across all endemic african countries for each year from 2000 to 2015 . using the empirically observed effect of each intervention , it was possible to generate counterfactual maps estimating contemporary pfpr under hypothetical scenarios without interventions . we chose to evaluate this  no intervention counterfactual to allow estimation of the total effect of interventions . for the 32 high - burden countries of africa , an ensemble model was developed to predict incidence rates of clinical malaria as a function of community pfpr . this brought together three independently developed mathematical malaria transmission models that were re - fitted to a common data set of age - structured clinical incidence measured longitudinally at 30 sites , allowing an ensemble model to be defined to predict age - specific incidence at all locations given prevalence , seasonality , level of treatment , and probable immune status of populations . we used a definition of  clinical malaria as an attributable febrile episode ( body temperature in excess of 37.5 c ) , censored by a 30-day window ( that is , multiple bouts of symptoms occurring within the same 30-day period are counted as a single episode ) . the ensemble model was then combined with the pfpr cube and underlying population surfaces to predict clinical incidence by country and year for both the real and counterfactual scenarios . for the remaining eleven low - burden countries in africa ( accounting for around 3% of cases ) where national reporting systems are more robust , we generated clinical incidence estimates with an existing approach that uses national case reports while adjusting for care - seeking behaviour , low diagnostic testing rates , and underreporting . we found that infection prevalence in children age 2-up - to-10 across endemic africa has halved since the year 2000 ( population - weighted mean pfpr2 - 10 : year 2000 = 33% , 95% credible interval 3135% ; year 2015 = 16% , 1419% ) , with around three - quarters of this decline occurring after 2005 . across africa the rate of decline in pfpr2 - 10 rose steadily to a peak yearly decline of 9% in 2011 , after which there was a slowing between 2011 and 2013 followed by resurgence in recent years back to the current rate of 5% annual decline in pfpr . our predicted surfaces of pfpr2 - 10 demonstrate the geographical pattern of this reduction across the continent ( fig . c ) , with hyper- or holo - endemic transmission ( where pfpr exceeds 50% and 75% , respectively , see fig . 1d ) common in 2000 across large swathes of central and western africa , but limited to isolated pockets by 2015 . 1d ) , with the proportion of the endemic population exposed to hyper- or holo - endemic malaria falling from 33% ( 30 - 37% ) to just 9% ( 5 - 13% ) ( table 1 ) . crucially , for the feasibility of post-2015 elimination efforts , the population of stable endemic africa experiencing very low transmission ( pf pr2 - 10 less than 1% ) has increased six - fold since 2000 ( far outpacing the 50% underlying population growth over the period ) meaning there are now 121 ( 110 - 133 ) million people living in settings where elimination campaigns can be considered . we estimated that there were 187 ( 132259 ) million clinical cases of p. falciparum malaria in africa in 2015 . case incidence declined by 40% from 321 ( 253427 ) per 1000 persons per annum in 2000 to 192 ( 135265 ) per 1000 persons p.a . in 2015 , with all but one of the 43 mainland endemic countries meeting the mdg target of reversing incidence trends by 2015 , 19 ( 17 - 25 ) achieving a > 50% decline , and 7 ( 6 - 7 ) declining by > 75% ( extended data fig . the model has been able to predict changes in mean pfpr across africa with considerable precision , reflecting the increasing abundance of pfpr surveys , their relatively large signal - to - noise ratio , the adjustments for diagnostic type , and the informative covariate suite , all of which contributed to strong predictive performance of the geostatistical model . credible intervals around the continental clinical incidence estimates were proportionately much larger and this reflected primarily the residual uncertainty around the modelled relationship between infection prevalence and clinical incidence . changes in prevalence largely followed patterns of increasing itn coverage , and itns were by far the most important intervention across africa , accounting for an estimated 68 ( 62 - 72)% of the declines in pfpr seen by 2015 ( fig . 2a ) . we estimated act and irs contributed 19 ( 15 - 24)% and 13 ( 11 - 16)% respectively , although these interventions had larger proportional contributions where their coverage was high ( extended data fig . it is important to emphasise that these proportional contributions do not necessarily reflect the comparative effectiveness of different intervention strategies but , rather , are driven primarily by how early and at what scale the different interventions were deployed . in total , we estimated that malaria control interventions have averted 663 ( 542753 ) million clinical cases since 2000 , of which 68 ( 62 - 73)% , 22 ( 17 - 28)% and 10 ( 5 - 14)% were contributed by itns , acts , and irs , respectively ( fig . 2b ) . here , for the first time , the rapidly changing landscape of malaria risk in africa has been quantified across the 15-year span of the millennium development goals . our approach is primarily data driven , informed by empirical observations in the field rather than theoretical models or extrapolated experimental results . our modelling framework requires few prior assumptions and allows patterns of change and attribution to be identified with rigorously defined metrics of uncertainty . we have shown that remarkable and widespread reductions in infection prevalence and case incidence have occurred across africa since 2000 , and that malaria control interventions have been responsible for most of the decline even though they remain well below international targets for universal coverage . itns have had by far the largest effect , but have also been generally present for longer and at higher levels of coverage . irs and acts have both made important contributions to reducing prevalence and incidence where they have been implemented at scale ( although it is important to note that the primary role of acts is in averting severe disease and death rather than reducing transmission and uncomplicated cases ) . this analysis has focused on evaluating changes in infection prevalence and clinical incidence , and we have not addressed effects on malaria mortality . data on malaria deaths are sparse both spatially and temporally , and concerted efforts must be made to both increase data collection and improve the sensitivity and specificity of malaria death attribution . integrating the results of the current study with existing malaria mortality estimation processes , to yield improved understanding of lives saved by malaria control , is an immediate priority . the modelling framework presented here has been necessitated in part by the absence of detailed and robust surveillance data collected routinely by health systems across africa . the development of more robust surveillance systems to deliver geographically detailed and timely data on malaria incidence will be an increasingly important strand of malaria control efforts , particularly if prevalence continues to decline and identification and rapid response to individual malaria cases becomes critical to achieve elimination . the efforts of the international community over the past 15 years have reduced malaria risk levels for many millions of people , and large regions of africa are now in a position to consider elimination strategies . despite this progress , many millions of people remain at risk of malaria disease and death in africa in 2015 . this analysis demonstrates that current malaria interventions have been highly effective at reducing prevalence and incidence across the continent , and provides strong support for sustaining and increasing access to these interventions as a cornerstone of post-2015 control strategies . this will need to be coupled with a redoubling of efforts to delay the spread of drug and insecticide resistance , tools for addressing the residual transmission that persists in some regions despite high vector control coverage , and concerted local programs to systematically detect and eliminate the remaining parasites . curves illustrate the predicted effect of itns as a function of coverage ( five example coverage levels are shown , specified as mean coverage over preceding 4-year period ) and baseline transmission . the baseline pfpr is shown on the horizontal axis and the suppressed pfpr given the itn coverage level shown on the vertical axis . the diagonal line ( representing zero itn effect ) is shown in black , and parameter uncertainty around each itn effect line is illustrated by the semi - transparent envelopes . results shown are derived from a bayesian geostatistical model fitted to : n = 27,573 pfpr survey points ; n = 24,868 itn survey points ; n = 96 national survey reports of act coverage ; n = 688 country - year reports on itn , act and irs distribution by national programs ; and n = 20 environmental and socioeconomic covariate grids . estimated country - level rates of all - age clinical incidence results shown are derived from a bayesian geostatistical model fitted to : n = 27,573 pfpr survey points ; n = 24,868 itn survey points ; n = 96 national survey reports of act coverage ; n = 688 country - year reports on itn , act and irs distribution by national programs ; n = 20 environmental and socioeconomic covariate grids ; and n = 30 active - case detection studies reporting p. falciparum clinical incidence . d , each map shows absolute decline in pfpr2 - 10 between 2000 and 2015 within areas of stable transmission attributable to the combined effect of itns , acts , and irs ( a ) ; and the individual effect of itns ( b ) ; acts ( c ) ; and irs ( d ) note that the colour scaling differs between the panels . results shown in all panels are derived from a bayesian geostatistical model fitted to : n = 27,573 pfpr survey points ; n = 24,868 itn survey points ; n = 96 national survey reports of act coverage ; n = 688 country - year reports on itn , act and irs distribution by national programs ; and n = 20 environmental and socioeconomic covariate grids . maps in this figure are available from the malaria atlas project ( http://www.map.ox.ac.uk/ ) under the creative commons attribution 3.0 unported license .
since the year 2000 , a concerted campaign against malaria has led to unprecedented levels of intervention coverage across sub - saharan africa . understanding the effect of this control effort is vital to inform future control planning . 
 however , the effect of malaria interventions across the varied epidemiological settings of africa remains poorly understood owing to the absence of reliable surveillance data and the simplistic approaches underlying current disease estimates . here 
 we link a large database of malaria field surveys with detailed reconstructions of changing intervention coverage to directly evaluate trends from 2000 to 2015 and quantify the attributable effect of malaria disease control efforts . 
 we found that plasmodium falciparum infection prevalence in endemic africa halved and the incidence of clinical disease fell by 40% between 2000 and 2015 . 
 we estimate that interventions have averted 663 ( 542753 credible interval ) million clinical cases since 2000 . 
 insecticide - treated nets , the most widespread intervention , were by far the largest contributor ( 68% of cases averted ) . 
 although still below target levels , current malaria interventions have substantially reduced malaria disease incidence across the continent . increasing access to these interventions , and maintaining their effectiveness in the face of insecticide and drug resistance , should form a cornerstone of post-2015 control strategies .