Heterogeneity of use, access, and retention of insecticide-treated nets: Implications for subnational tailoring to maximise malaria control

  1. Andrew C Glover  Is a corresponding author
  2. Hannah Koenker
  3. El Hadji Amadou Niang
  4. Kate Kolaczinski
  5. Thomas S Churcher
  1. Imperial College London, United Kingdom
  2. Program for Appropriate Technology in Health, United States
  3. Laboratoire d’Écologie Vectorielle et Parasitaire, Département de Biologie Animale, Université Cheikh Anta Diop, Senegal
  4. The Global Fund to Fight AIDS, Tuberculosis and Malaria, Switzerland
10 figures, 5 tables and 1 additional file

Figures

Population-weighted distribution of expected ITN retention times and duration of use by country.

Median estimates of the mean ITN retention time and duration of use for each subnational region were weighted by population to produce country-level distributions. These histograms represent probability distributions of retention time and duration of use for an individual randomly selected from each country, assuming model-derived median estimates for all subnational regions. Subnational variability in retention times and the duration of use reflects differences in model-derived estimates between regions, since the model assumes these values are constant within regions. Vertical dashed lines indicate the overall population-weighted mean (black) and the population-weighted mean for each country (coloured).

Figure 2 with 2 supplements
ITN use, access, and PfPR over time in rural Burkina Faso.

(A) Survey estimates of ITN use (probability someone slept under an ITN the previous night, points sized proportionally to sample size) and model mean estimates of the probability of use over time (solid lines) for rural portions of subnational regions in Burkina Faso. Mean estimates are shown for the probability of access (dotted lines) and the probability of using an ITN sourced from a continuous channel (dashed lines). The model accounts for differing probabilities of use between individuals within a region; this results in notable uncertainty around the probability an individual selected at random uses an ITN. This uncertainty is illustrated in the shaded region, with the 50%, 80%, and 95% credible intervals as indicated by progressively lighter shaded regions. (B) PfPR2−10 estimates from the transmission dynamics model (darker lines) were fitted to annual Malaria Atlas Project estimates of PfPR2−10 compiled by site (Winskill, 2024b) (hollow points); model estimates of PfPR6−59mo are also shown (lighter lines), in addition to observed DHS PfPR6−59mo and associated 95% credible intervals due to measurement uncertainty (solid points and vertical lines).

Figure 2—figure supplement 1
Use, access, and PfPR over time in rural Senegal.

Model estimates in a lower transmission context in Senegal are shown. All features remain unchanged from Figure 2. Note, in B, some credible intervals extend beyond the y-axis limits, which have been restricted to 30% to aid visibility.

Figure 2—figure supplement 2
Use of ITNs from continuous channels.

For each subnational region, coloured points indicate median estimates of the mean use of any ITN and of those sourced from continuous channels, immediately following a mass campaign (A) and over the subsequent 3 years (B) with associated 95% credible intervals. For each country, the proportion of ITNs that were used which were sourced from continuous channels immediately following a campaign (A) and over the subsequent 3 years (B) can be estimated from the slopes of the coloured population-weighted linear regression lines, which assume a zero intercept. Estimates over all countries are shown by the population-weighted linear regression line in black.

Figure 3 with 3 supplements
Mean use and access with 3-year campaigns in subnational rural areas.

Central estimates of mean overall proportion of people using an ITN the previous night, access to an ITN, and use given access for 3-year mass campaign intervals are shown in the top three rows. The number of months following a mass campaign where overall ITN use exceeds 80% are shown in the bottom row.

Figure 3—figure supplement 1
Mean use and access with 3-year campaigns in subnational urban areas.

Central estimates of mean overall proportion of people using an ITN the previous night, access to an ITN, and use given access for 3-year mass campaign intervals are shown in the top three rows. The number of months following a mass campaign where overall ITN use exceeds 80% are shown in the bottom row.

Figure 3—figure supplement 2
Mean rural use and access with 2-year campaigns.

Central estimates of mean overall subnational ITN use, access, and use given access for 2-year mass campaign intervals in rural areas are shown in the top three rows. The number of months following a mass campaign where overall ITN use exceeds 80% are shown in the bottom row.

Figure 3—figure supplement 3
Mean urban use and access with 2-year campaigns.

Central estimates of mean overall subnational ITN use, access, and use given access for 2-year mass campaign intervals in urban areas are shown in the top three rows. The number of months following a mass campaign where overall ITN use exceeds 80% are shown in the bottom row.

Figure 4 with 2 supplements
Equity of use and access in rural Thiès and Ziguinchor, Senegal.

Solid vertical lines denote overall ITN use and access averaged over 3 years following the last mass campaign. Coloured bars indicate the proportion of the population with different probabilities of using or having access to an ITN. For example, 10% of the population have a probability of access between 70% and 80% in rural Ziguinchor, in comparison to 20% in rural Thiès. Vertical dashed lines denote 50% credible intervals for the probability of an individual using or having access to an ITN, indicating that half of the population are expected to have a probability of use or access within this range.

Figure 4—figure supplement 1
Subnational equity of use in Senegal.

All features remain unchanged from Figure 4.

Figure 4—figure supplement 2
Subnational equity of access in Senegal.

All features remain unchanged from Figure 4.

Figure 5 with 5 supplements
Retention time, use given access, and changes in cases for Mali.

All coloured points indicate median estimates and associated 95% credible intervals for subnational regions in Mali. (A) Mean durations of ITN use (hollow) and retention (filled) are shown in addition to national median estimates (black solid and dashed lines, respectively). (B) Mean ITN retention and use given access under a triennial campaign strategy are shown, with black lines indicating national median estimates. (C) Points indicate mean annual clinical cases under a triennial pyrethroid-PBO campaign strategy, and the projected change in cases for alternative intervention strategies with pyrethroid-only (C.i–iii), pyrethroid-PBO (C.iv–vi), or pyrethroid–chlorfenapyr (C.vii–ix) ITNs, with continuous-only distribution (C.i, iv, vii), or in conjunction with triennial (C.ii, v, viii) or biennial (C.iii, vi, ix) campaigns. The change in clinical cases is equal to zero for equivalent comparator and intervention strategies (C.v). Labelled diagonal reference lines with positive gradients indicate percentage increases in clinical cases relative to mean annual estimates under the comparator strategy; lines of the same style with negative gradients indicate equivalent percentage decreases. Urban areas were estimated to broadly have lower use given access than in rural settings, while the capital city, Bamako, was estimated to have the lowest ITN retention time of all regions.

Figure 5—figure supplement 1
Retention time, use given access, and changes in cases for Burkina Faso.

All features remain unchanged from Figure 5. Note the ordering of colours differs in panel B where numerical labels should be referred to.

Figure 5—figure supplement 2
Retention time, use given access, and changes in cases for Ghana.

All features remain unchanged from Figure 5. In 2019, several administrative regions in Ghana were subdivided: Brong-Ahafo became Bono, Bono East, and Ahafo; the Northern Region was split into Northern, Savannah, and North East; and Volta was divided into Volta and Oti. Retrospective analyses (A, B) were conducted using the pre-2019 regional boundaries, while future projections (C) used the post-2019 subdivisions. In all subnational regions, ITN retention times and use given access were estimated to be higher in rural than in urban settings.

Figure 5—figure supplement 3
Retention time and use given access for Malawi.

All features remain unchanged from Figure 5A. Clinical case projections were not conducted for Malawi due to the lack of intervention data stratified to the administrative-one level, despite DHS surveys being conducted at that scale.

Figure 5—figure supplement 4
Retention time, use given access, and changes in cases for Mozambique.

All features remain unchanged from Figure 5.

Figure 5—figure supplement 5
Retention time, use given access, and changes in cases for Senegal.

All features remain unchanged from Figure 5. While use given access was estimated to be broadly similar across urban and rural settings, ITNs were generally estimated to be retained for longer in rural areas.

Figure 6 with 9 supplements
Change in cases vs cases with triennial pyrethroid-PBO distribution.

Points represent urban and rural areas within each subnational region, are coloured by country, and are sized proportionally to the mean annual number of ITNs distributed under each strategy. Each point shows median estimates of the change in clinical cases following a switch from triennial pyrethroid-PBO distribution to alternative strategies against the projected clinical cases under the triennial pyrethroid-PBO strategy. Concurrent distribution of continuous ITNs of the same class used in mass campaigns in assumed throughout. Vertical and horizontal lines indicate 95% credible intervals. Labelled diagonal reference lines with positive gradients indicate percentage increases in clinical cases relative to mean annual estimates under the comparator strategy; corresponding lines with negative gradients represent equivalent percentage decreases and share the same line style.

Figure 6—figure supplement 1
Change in cases vs cases with triennial pyrethroid-chlorfenapyr distribution.

Points show median estimates of the change in clinical cases following a switch from triennial pyrethroid-chlorfenapyr distribution to alternative strategies against projected clinical cases under the triennial pyrethroid-chlorfenapyr strategy. Figure features otherwise remain unchanged from Figure 6.

Figure 6—figure supplement 2
Projected cases averted under different ITN distribution strategies.

Points indicate median estimates of the mean annual clinical cases averted by different ITN distribution strategies, in comparison to ceasing ITN distribution entirely, against clinical cases under a triennial pyrethroid-PBO distribution strategy. Figure features otherwise remain unchanged from Figure 6.

Figure 6—figure supplement 3
Change in cases vs PfPR with triennial pyrethroid-PBO distribution.

Points show median estimates of the change in clinical cases following a switch from triennial pyrethroid-PBO distribution to alternative strategies against all-age PfPR under the triennial pyrethroid-PBO strategy. Figure features otherwise remain unchanged from Figure 6.

Figure 6—figure supplement 4
Change in cases vs pyrethroid resistance.

Points show median estimates of the change in clinical cases following a switch from triennial pyrethroid-PBO distribution to alternative strategies against pyrethroid resistance. Figure features otherwise remain unchanged from Figure 6.

Figure 6—figure supplement 5
Change in cases vs duration of use.

Points show median estimates of the change in clinical cases following a switch from triennial pyrethroid-PBO distribution to alternative strategies against mean duration of use. Figure features otherwise remain unchanged from Figure 6.

Figure 6—figure supplement 6
Change in cases vs retention time.

Points show median estimates of the change in clinical cases following a switch from triennial pyrethroid-PBO distribution to alternative strategies against mean retention time. Figure features otherwise remain unchanged from Figure 6.

Figure 6—figure supplement 7
Change in cases vs use.

Points show median estimates of the change in clinical cases following a switch from triennial pyrethroid-PBO distribution to alternative strategies against mean use over any triennial distribution strategy. Figure features otherwise remain unchanged from Figure 6.

Figure 6—figure supplement 8
Change in cases vs access.

Points show median estimates of the change in clinical cases following a switch from triennial pyrethroid-PBO distribution to alternative strategies against mean access over any triennial distribution strategy. Figure features otherwise remain unchanged from Figure 6.

Figure 6—figure supplement 9
Change in cases vs use given access.

Points show median estimates of the change in clinical cases following a switch from triennial pyrethroid-PBO distribution to alternative strategies against mean use given access over any triennial distribution strategy. Figure features otherwise remain unchanged from Figure 6.

Appendix 1—figure 1
Flowchart of methods.

Data inputs and methodological steps are, respectively, summarised in orange and blue. 1DHS data (ICF, 2025) were sourced using the RDHS R package (v0.8.4) (Watson et al., 2019). 2Annual estimates of the number of ITNs delivered by country were obtained from the (The Alliance for Malaria Prevention, 2024) Net Mapping Project. 3Region-specific characteristics were sourced from the site R package (v0.2.2) (Winskill, 2024b). 4Parameter values that describe the probability of repellency and mortality as different ITN types age were sourced from Sherrard-Smith et al., 2022; Churcher et al., 2024 conditioned on annual pyrethroid (Pyr) resistance estimates for each subnational region from site. 5Mean annual PfPR2−10 Malaria Atlas Project estimates were sourced from the site R package (v0.2.2) (Winskill, 2024b); 6these were used to calibrate baseline entomological inoculation rates (EIR) with the cali (v1.0.8) R package (Winskill, 2024a) given historical ITN and non-ITN interventions, with coverage estimates of the latter sourced from site. 7Case projections were simulated for different ITN distribution strategies (D) using the malariasimulation R package (v1.6.0) (Charles et al., 2024) following characterisation of subnational differences in ITN use, access, and retention (A–C).

Appendix 2—figure 1
An illustration of the number of nets, n(t), in a region over time with continuous replenishment of routinely distributed nets, nd, where nc0 nets are distributed every Δt years through regular mass campaigns (solid black lines).

The contribution towards the total number of nets from routine distribution and previous mass campaigns are shown by the dashed and dotted blue lines, respectively. The timing of DHS surveys are shown by the red vertical dashed lines for both regular (left) and irregular (right) surveying.

Appendix 2—figure 2
Posterior predictive distributions (red) from the hierarchical model for the age of a used net in months if one is sampled randomly over time and across either sub-Saharan Africa (A), or within an exemplar country (B) or subnational region (C).

The shaded regions indicate 95% credible intervals, while the blue lines indicate the normalised empirical densities; these were generated from the pseudo-counts as described in Equation 19. The peaks at 12, 24, and 36 months in the empirical densities are believed to be artifacts from the survey data due to individuals rounding reported ITN ages to the nearest whole year.

Appendix 3—figure 1
Estimates of mass campaign timings in the subnational Est region in Burkina Faso.

Letting region i be the Est region of Burkina Faso, (a) shows the empirical,Zζ(i), and smoothened,Zζ(i), annual numbers of nets delivered nationally are shown by the solid blue line (as reported by the Alliance of Malaria Prevention, AMP) and shaded region. The local minima of Zζ(i), which define the midpoints between net delivery dates are shown by the vertical black lines in (a) and (d) and indicate the assumed periods when mass campaigns are allowed to have occurred. The density proportional to the number of nets received subnationally by month, Xi, is shown in red (b). The weighted version of this density, Wi, to account for older nets being under-represented in DHS surveys is shown in yellow (c). The composite density, Vi, and its smoothened counterpart, Vi, which are constructed from the stepwise densities in (a) and (c) are shown by the green line and shaded region in (d). The densities, Yik, are subsetted from Vi by the black lines. These are treated as approximations of the probability of a mass campaign occurring in each month, and the expected value for the timing of each kth mass campaign are shown by the black triangles.

Tables

Appendix 1—table 1
Core notation used throughout the main text and all appendices.
SymbolDefinition
iSubnational region index.
jTime point index.
kMass campaign index.
lITN or individual index, depending on context.
tContinuous time variable.
tjDiscrete time point corresponding to month j.
ζ(i)Country to which region i belongs.
aijProportion of the population with ITN access in region i at time tj.
uijProportion of the population with ITN use in region i at time tj.
τikTiming of the kth mass campaign in region i.
ϕijTiming of the most recent mass campaign in region i before time tj.
mijNumber of months since the most recent mass campaign in region i at time tj.
pijProportion of the population with ITN use or access in region i at time tj.
cijComponent of pij attributable to campaign-distributed ITNs.
dijComponent of pij attributable to continuously distributed ITNs.
pij0Immediate post-campaign level of ITN use or access in region i at time tj.
qijProportion of ITN use or access attributable to campaign-distributed ITNs in region i at time tj.
qi0Proportion of immediate post-campaign ITN use or access attributable to campaign-distributed ITNs in region i.
ψiUpper asymptote for pij0 in region i.
βi0Region-specific intercept parameter governing the logistic-type trajectory in immediate post-campaign use or access.
βitRegion-specific slope parameter governing the logistic-type trajectory in immediate post-campaign use or access over calendar time.
λiRegion-specific exponential decay rate for the campaign-derived component of ITN use or access.
γi1Mean duration of ITN use or access (retention time) in region i under the assumption of non-random allocation of ITNs.
Λi1Mean duration of ITN use or access (retention time) in region i under the assumption of random allocation of ITNs.
p~ijIndividual-level probability of ITN use or access in region i at time tj.
αi0Region-specific overdispersion parameter governing heterogeneity in the individual-level probability p~ij.
αiAge distribution of ITNs providing use or access in region i.
αilAge of ITN l in region i at the point of observation.
wilDHS household survey weight for ITN l in region i, rescaled to remove the DHS factor of 106.
Appendix 1—table 2
Additional notation used in Appendix 2 beyond the core notation defined in Appendix 1—table 1.
SymbolDefinition
x{u,a}Superscript indicating whether a quantity refers to ITN use (u) or access (a).
ρilThe duration ITN l in region i continues to be used or provide access.
γi1Mean duration ITNs in region i continue to be used or provide access.
N(t)Time-dependent replenishment process through which nets providing use or access are distributed at time t.
N¯Long-run mean replenishment rate of nets providing use or access.
TTime horizon over which the long-run average replenishment rate is defined.
ΔtTime period of the periodic equilibrium replenishment process.
ΔsInterval between surveys when surveys are assumed to occur at regular intervals.
rSpecific ITN age value used in the derivation of the pooled age density.
ΔrSmall age interval width used when considering ITNs with ages in [r,r+Δr].
αiVector of observed ITN ages (αil) pooled across a finite number of surveys.
γ^Empirical estimate of the exponential decay rate obtained by fitting the pooled ITN age distribution.
γζ(i)1Country-specific mean ITN retention time for the country containing region i.
σζ(i)Country-specific standard deviation governing between-region variation in mean ITN retention times within country ζ(i).
γ01Continent-level mean ITN retention time in the hierarchical model.
σ0Between-country standard deviation in mean ITN retention time.
χ0Hyperhyperprior mean for the country-level standard deviations σζ(i).
ω0Hyperhyperprior standard deviation for the country-level standard deviations σζ(i).
αilObserved age of ITN l in region i in the hierarchical likelihood.
MirPseudo-count of the number of ITNs of age α=r in region i, after applying DHS household survey weights and rounding.
NirNumber of surveyed ITNs of age r in region i.
Appendix 1—table 3
Additional notation used in Appendix 3 beyond the core notation defined in Appendix 1—table 1.
SymbolDefinition
Zζ(t)Annual number of ITNs delivered to country ζ at time t.
Zζ(t)Smoothened distribution of annual ITN deliveries to country ζ at time t.
Xi(t)Normalised density of the month in which ITNs recorded in DHS surveys in region i were received.
hilSurvey date of ITN l in region i, measured in months.
Wi(t)Weighted density of the month in which ITNs in region i were received, after correcting for the under-representation of older ITNs in DHS surveys.
γ¯i1Posterior mean estimate from the hierarchical model of the mean duration of access (retention time) in region i under the assumption of non-random allocation of ITNs.
Vi(t)Composite density for region i constructed from Wi(t) and the national annual-delivery density Zζ(i)(t).
Vi(t)Smoothened version of the composite density Vi(t) for region i.
Yik(t)Normalised density for the timing of the kth mass campaign in region i.
τ¯ikMean estimate of the timing of the kth mass campaign in region i.
στikAssigned standard deviation for the timing of mass campaign k in region i, taking values 1, 2, or 3 months to represent low, medium, or high uncertainty.
δUnit point mass at t=.
Appendix 1—table 4
Additional notation used in Appendix 4 beyond the core notation defined in Appendix 1—table 1.
SymbolDefinition
nijTotal number of surveyed individuals in region i at time tj.
pi(t)Continuous-time analogue of pij.
ci(t)Continuous-time analogue of cij.
di(t)Continuous-time analogue of dij.
pi0(t)Continuous-time analogue of pij0.
mi(t)Continuous-time analogue of mij.
μλiPosterior mean of the mean ITN retention time in region i, used to parameterise the prior on λi1.
σλi2Posterior variance of the mean ITN retention time in region i, used to parameterise the prior on λi1.
t^jStandardised survey time corresponding to time tj.
t¯Mean survey time used to standardise tj.
σtStandard deviation of survey times used to standardise tj.
β^i0Intercept parameter on the standardised time scale.
β^itSlope parameter on the standardised time scale.
yijNumber of ITNs in region i at time tj recorded or imputed as originating from campaigns.
υilIndicator variable for whether ITN l in region i originated from a campaign (υil=1) or continuous distribution (υil=0).
υilIndicator variable denoting that the source of ITN l in region i was missing in the DHS data and was therefore imputed probabilistically.
y¯iObserved proportion of ITNs providing access in region i that were recorded as originating from campaigns.
d0Proportion of the population using continuously distributed ITNs at the start of a month in the random-allocation approximation.
d1ΔtProportion of the population using continuously distributed ITNs immediately before the start of the next month.
ΔdProportion of the population that must receive ITNs through continuous distribution at the start of the next month under random allocation.
p0Proportion of the population using ITNs immediately following an arbitrary mass campaign in the random-allocation approximation.
p^(m)Approximate exponential-decay solution for population ITN use as a function of time since the most recent campaign.
Appendix 1—table 5
Additional notation used in Appendix 5 beyond the core notation defined in Appendix 1—table 1.
SymbolDefinition
ϵiRegion-specific over-reporting parameter applied to ITN use in region i.
PijTrue value of PfPR6−59mo at time tj in region i.
γNDecay rate governing loss of ITN insecticidal efficacy with ITN age.
rNProbability that a mosquito is repelled by an ITN of age α, conditional on a feeding attempt.
dNProbability that a mosquito is killed by an ITN of age α, conditional on a feeding attempt.
sNProbability that a mosquito successfully feeds in the presence of an ITN of age α, conditional on a feeding attempt.
rN0Probability that a mosquito is repelled by a new ITN.
dN0Probability that a mosquito is killed by a new ITN.
rNMLong-term probability of repellency from the physical barrier of the ITN alone.
τLCentral estimate of the timing of the last mass campaign before 2025.
ΔτInterval between future mass campaigns.
τFTiming of the first future mass campaign.
xij+Number of individuals testing positive for malaria in the calibration dataset.
xijNumber of individuals testing negative for malaria in the calibration dataset.
nij+Total number of individuals in the calibration dataset contributing positive test outcomes.
nijTotal number of individuals in the calibration dataset contributing negative test outcomes.
wijlDHS household survey weight for child l in region i at time tj, rescaled to remove the DHS factor of 106.

Additional files

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Andrew C Glover
  2. Hannah Koenker
  3. El Hadji Amadou Niang
  4. Kate Kolaczinski
  5. Thomas S Churcher
(2026)
Heterogeneity of use, access, and retention of insecticide-treated nets: Implications for subnational tailoring to maximise malaria control
eLife 14:RP108745.
https://doi.org/10.7554/eLife.108745.4