Characterizing human mobility patterns in rural settings of sub-Saharan Africa

  1. Hannah R Meredith  Is a corresponding author
  2. John R Giles
  3. Javier Perez-Saez
  4. Théophile Mande
  5. Andrea Rinaldo
  6. Simon Mutembo
  7. Elliot N Kabalo
  8. Kabondo Makungo
  9. Caroline O Buckee
  10. Andrew J Tatem
  11. C Jessica E Metcalf
  12. Amy Wesolowski  Is a corresponding author
  1. Department of Epidemiology, Johns Hopkins Bloomberg School of Public Health, United States
  2. Bureau d'Etudes Scientifiques et Techniques - Eau, Energie, Environnement (BEST-3E), Burkina Faso
  3. Dipartimento di Ingegneria Civile Edile ed Ambientale, Università di Padova, Italy
  4. Laboratory of Ecohydrology, School of Architecture, Civil and Environmental Engineering, École Polytechnique Fédérale de Lausanne, Switzerland
  5. Department of International Health, International Vaccine Access Center, Johns Hopkins Bloomberg School of Public Health, United States
  6. Macha Research Trust, Zambia
  7. Zambia Information and Communications Technology Authority, Zambia
  8. Zamtel, Zambia
  9. Department of Epidemiology and the Center for Communicable Disease Dynamics, Harvard TH Chan School of Public Health, United States
  10. WorldPop, School of Geography and Environmental Science, University of Southampton, United Kingdom
  11. Department of Ecology and Evolutionary Biology and the Princeton School of Public and International Affairs, Princeton University, United States

Abstract

Human mobility is a core component of human behavior and its quantification is critical for understanding its impact on infectious disease transmission, traffic forecasting, access to resources and care, intervention strategies, and migratory flows. When mobility data are limited, spatial interaction models have been widely used to estimate human travel, but have not been extensively validated in low- and middle-income settings. Geographic, sociodemographic, and infrastructure differences may impact the ability for models to capture these patterns, particularly in rural settings. Here, we analyzed mobility patterns inferred from mobile phone data in four Sub-Saharan African countries to investigate the ability for variants on gravity and radiation models to estimate travel. Adjusting the gravity model such that parameters were fit to different trip types, including travel between more or less populated areas and/or different regions, improved model fit in all four countries. This suggests that alternative models may be more useful in these settings and better able to capture the range of mobility patterns observed.

Introduction

Human mobility patterns are a reflection of behaviors, ranging from routine (e.g., commuting daily for work or school, traveling for holidays and religious gatherings, or seeking seasonal work opportunities) to irregular (e.g., relocating due to environmental changes or crises or social distancing due to a pandemic) (Charaudeau et al., 2014; Haberfeld et al., 1999; International Organization for Migration, 2019; Lessler et al., 2014; Lu et al., 2012; Pullano et al., 2020). Consequently, characterizing human mobility patterns is important for a wide range of applications, including predicting the spread of infectious diseases, traffic forecasting, designing intervention strategies, assessing health service accessibility, planning natural disaster relief efforts, and estimating migratory flows to estimate changes in population demographics (Charaudeau et al., 2014; Dotse-Gborgbortsi et al., 2020; Findlater and Bogoch, 2018; Finger et al., 2016; Gilbert et al., 2020; Lu et al., 2012; Mari et al., 2017; Palchykov et al., 2014; Stoddard et al., 2009; Wesolowski et al., 2012). Depending on the spatial and temporal resolution of travel relevant to the question of interest, there are many data sources that can quantify travel, such as national censuses, traffic or commuting data, travel surveys, or mobile phone data (Tatem, 2014). However, when mobility data are limited or unavailable, spatial interaction models are often used to estimate mobility patterns. One of the most commonly used models is the gravity model (Wesolowski et al., 2015b), which assumes the number of trips between an origin and destination over a fixed time period will increase as a function of the destination and origin population sizes and decrease with distance. While the gravity model was originally developed to describe commuting in a high-income setting (Zipf, 1946), it has been used to model travel in low- and middle-income countries (LMICs) when data are limited or unavailable (Stone et al., 2019; Wells et al., 2019; Wesolowski et al., 2015b). In high-income settings, the standard gravity model has been shown to perform well when predicting commuter movement between cities (Masucci et al., 2013) and perform poorly when predicting movement across areas with heterogeneity in demographics and population density or in rural areas (Truscott and Ferguson, 2012; Xia et al., 2004). However, the degree to which geographic, sociodemographic, and infrastructure differences may impact the ability for models to capture travel patterns in LMICs, particularly in more rural settings, remains to be determined.

Increasing mobile phone ownership (over 67% of the global population in 2019) and the mobility patterns that can be extracted from these data provide a valuable resource for evaluating how well spatial interaction models can capture human mobility patterns in low- and middle-income settings (GSM Association, 2020; Wesolowski et al., 2016). However, the challenge of procuring mobile phone datasets has resulted in relatively few studies that fit mobility models to relevant travel data from LMICs. Of these, a few studies have explored model adjustments and show that predictions can be improved by adjusting the gravity model to account for factors such as education levels, economic opportunities, gender, environmental factors, trip duration, contiguity of origin and destination, and proportion of population living in urban areas (Garcia et al., 2015; Henry et al., 2003; Wesolowski et al., 2015b). Typically, studies assume that these factors impact all trips homogeneously, thus fitting a single set of parameters for all trips is a sufficient adjustment. Yet, even with such adjustments, gravity models may still fail to accurately capture mobility patterns in LMICs, particularly in rural areas (Henry et al., 2003; Wesolowski et al., 2015b), suggesting that models which account for regional heterogeneity in travel may improve their estimates in these settings.

Here, we examine the impact that regional heterogeneity and urbanicity has on travel patterns as well as how well spatial interaction models can reproduce travel patterns in LMICs by examining four countries in Sub-Saharan Africa: Namibia, Kenya, Burkina Faso, and Zambia (Figure 1A–D). These countries capture a range of levels of urbanization, population distribution, and importantly include large rural areas common in many LMICs. Using mobile phone data, we first characterized human mobility patterns from these four countries and determined which features were not well captured by basic gravity models. Next, we conducted an analysis of six variations of the basic gravity model, including variations that account for regional travel and urbanicity, as well as a basic radiation model to determine which provided the best trip estimates for each country (Simini et al., 2012). Finally, we compared the different model fits across the four countries to evaluate if the adjustments produced similar improvements for all countries. To our knowledge, no other models have captured these mobility features, been tested consistently against mobility data from multiple LMICs, and provided clear guidance on which model to use and when. This study provides further insight on the mobility patterns in LMICs and highlights where mobility model estimates may deviate when applied in other similar settings, ultimately improving our understanding of topics such as disease spread, migratory flows, and intervention efficacy in LMICs.

Figure 1 with 9 supplements see all
The mobility patterns extracted from mobile phone data from four countries in Sub-Saharan Africa.

(A–D) Data from four Sub-Saharan African countries were selected to characterize human mobility patterns: Namibia, Kenya, Burkina Faso, and Zambia. Travel between districts (administrative level 2) was estimated via mobile phone data from each country. A basic gravity model was fit to trip data from each country which assumes that the number of trips decreases with distance and increases with population size (E–H). Here, one rural (left panel) and one urban (right panel) destination were selected from each country to show that, while the observed trips (black) from different origins do generally follow the assumptions of the gravity model (red), the gravity model is not fully capturing the observed trip patterns. See Figure 1—figure supplements 18 for this comparison made for all districts in each country. (I–L) Comparisons of origin-destination matrices colored by trip proportions estimated by the basic gravity model and the mobile phone data (observed) highlight how the basic gravity model tends to overestimate many trips, particularly those that are off-diagonal (e.g., inter-regional trips). The columns and rows of the OD matrix are ordered by district ID, which were assigned such that districts within the same region (adminstrative level 1) were clustered together. The capital district is indicated by the black arrow on the x- and y-axes. The colors indicate the proportion of an origin’s trips made to each destination (with light blue representing destinations visited infrequently and dark blue representing destinations visited most frequently). (See Supplementary file 1B for the key to the origin and destination numbers and Figure 1—figure supplement 9 for district level maps).

Results

We analyzed the average number of daily trips, defined here as the number of subscribers moving from one district (administrative level 2) to another, per month extracted from each mobile phone dataset (see Materials and methods). Similar to other studies, we found that the average number of monthly trips between districts from Namibia, Kenya, Burkina Faso, and Zambia generally decreased with trip distance and there were more trips from and to more populated areas (Figure 1E–H, Figure 1—figure supplements 18). Trips were concentrated between districts within the same region (administrative level 1) to varying degrees (intra-regional trips made up 30% of all trips in Burkina Faso, 45% in Kenya, 62% in Namibia, and 72% in Zambia) and to a few common destinations, including the district where the capital was located (Figure 1I–L, Supplementary file 1A). Although Namibia, Burkina Faso, and Zambia each consisted of ~95% predominantly rural districts, the distribution of monthly trips between urban and rural districts varied across countries. The majority of Namibia’s and Burkina Faso’s trips were between rural locations (62% and 70.5% of all trips, respectively), while Zambia’s trips were split between rural locations (53%) or rural and urban locations (46%). Kenya, with 56% predominantly urban districts, had the largest proportion of monthly trips between urban locations (70%). As a basic model, we first estimated trip counts using a basic gravity model (single parameter set fitted for all trips), which is based on the population sizes of the origin and destination and the distance between locations, with two variants on the distance kernel (power or exponential decay) (Equation 1). Estimates from this basic model overestimated trip counts and missed important features of the data, such as higher trip counts in short-distance and within-region trips, relative to long-distance and between-region trips (Figure 1E–L, Figure 1—figure supplements 18). Travel involving predominantly rural locations tended to be overestimated and these results were largely observed for both distance kernels and for all countries.

Given how the basic model’s estimates of urbanicity and regional movement deviated from the call data records, we tested six additional variations of a gravity model that allowed for parameters to capture increasingly complex features in the data (Equations 2–9). These included fitting parameters for the origin and destination population sizes and trip distance to: (1) trips defined by origin and destination population density (urbanicity model: higher population density = predominantly urban, lower population density = predominantly rural), (2) trips within- and between-regions (regional model), and (3) trips defined by both region and urbanicity (regional-urbanicity model) (Figure 2). Each model variation was tested with two distance kernels (power and exponential decay). We also evaluated a basic radiation model (Equation 11; Masucci et al., 2013; Simini et al., 2012).

Figure 2 with 2 supplements see all
Variations of the gravity model were fit to data to capture different types of trips.

Here, Kenya is used to demonstrate the trip types that could be defined by the region (represented by color) and/or urbanicity (solid = predominantly rural, dotted = predominantly urban) of the trip’s origin and destination to outline the various models fit. See Figure 2—figure supplements 12 for the fitted parameters.

Figure 2—source data 1

Table of model parameter values fit for each country, both distance kernels, and all trip types.

https://cdn.elifesciences.org/articles/68441/elife-68441-fig2-data1-v1.docx

All parameterized gravity models outperformed the basic model. Allowing for parameters to be fit by urbanicity or region improved model fit to varying degrees, depending on the country (Table 1, Supplementary file 2A). Accounting for urbanicity provided a larger improvement in model fit than region for all countries; however, this improvement was larger in Kenya and Burkina Faso than Namibia and was only observed for Zambia if the lower urbanicity threshold (10%) was implemented. In these models, the degree to which the parameters varied increased with model complexity, with the distance parameter varying the most, allowing for different weights to be applied based on the trip type (Figure 2—figure supplements 12). Overall, the most complex model, the regional-urbanicity model, had the best fit for all countries (the power variant reduced the basic model’s deviance information criterion (DIC) by 41% for Namibia, 30% for Kenya, 28% for Burkina Faso, and 16% for Zambia) (Table 1, Supplementary file 2A). Its flexibility allowed for the importance of origin and destination population sizes to vary by trip type, distinguished differences in the relative importance of distance by trip type, was able to better distinguish between inter- and intra-regional trip estimates, and adjusted for the lower trip counts between rural locations, relative to other trip types (Figure 3A, Figure 3—figure supplements 14). This was true regardless of the distance decay function assumed, although the best fitting decay function did vary by country (Supplementary file 2A). Interestingly, the radiation model had the most variable performance of all models, where it was the poorest performing in Namibia, Burkina Faso, and Zambia, but outperformed the basic gravity model (exponential form) in Kenya. Differences in the population distribution (more heterogenous in Namibia, Burkina Faso, and Zambia than Kenya) and locations of more populated areas may be the cause of these conflicting findings across countries (Figure 1—figure supplement 9; Linard et al., 2012; Masucci et al., 2013). These results were consistent across a range of administrative levels (1-3), suggesting that including locations’ urbanicity and/or region in a model will generally improve model fit across spatial scales (Supplementary file 2B).

Table 1
Gravity model variations (using power distance kernel) and radiation model, ranked for each country based on Deviance Information Criterion (DIC) and percent change (%Δ) from the basic gravity model.
RankNamibiaKenyaBurkina FasoZambia
ModelDIC (%Δ)ModelDIC (%Δ)ModelDIC (%Δ)ModelDIC (%Δ)
1Reg - Urb.3.62E + 06 (41.0)Reg - Urb.2.43E + 08 (30.1)Reg - Urb.1.93E + 05 (27.7)Reg - Urb.2.01E + 06 (16.3)
2Urbanicity4.56E + 06 (25.6)Urbanicity2.53E + 08 (27.2)Urbanicity2.05E + 05 (23.1)Regional2.08E + 06 (13.4)
3Regional4.61E + 06 (24.8)Regional3.40E + 08 (2.1)Regional2.52E + 05 (5.7)Urbanicity2.38E + 06 (1.1)
4Basic6.12E + 06 (0.0)Basic3.48E + 08 (0.0)Basic2.67E + 05 (0.0)Basic2.40E + 06(0.0)
5Radiation8.68E + 06 (-41.7)Radiation4.26E + 08 (-22.4)Radiation3.39E + 05 (-27.2)Radiation4.3E + 06 (-79.9)
Figure 3 with 7 supplements see all
The modeled estimates of trip counties and fits by type of trip.

(A) The proportion of trips from each origin location in Namibia was estimated by five different spatial models (power distance kernel displayed here) and ordered by the origin and destination ID. Regional clustering was more pronounced and there were fewer inter-regional trips in the adjusted models and radiation model, relative to the basic model. See Figure 3—figure supplements 14 for all countries and both distance kernel functions. The columns and rows of the OD matrix are ordered by district ID. (B) The ratio of predicted to observed trip counts in Namibia was calculated to determine the distribution of trips that were over- (ratio >1) or underestimated (ratio <1) in each trip type by model. The median ratio (solid black vertical line) for each trip type is compared with the equity line (white vertical line) for each model (shown as different colors). The proportion of trips that fall within the selected interval (± 10% of the observed trip count) was also used to assess a model’s ability to capture trips in that category. Generally, the basic model captured urban-to-urban trips the best and overestimated the other trip types. See Figure 3—figure supplements 57 for all countries, trip types, and distance kernel functions.

Overall model fit may not accurately describe how well each model can estimate certain trip types which may be relevant for particular questions. We further evaluated how close model estimates were to the observed trip count using the ratio of predicted to observed trip counts (Figure 3B, Figure 3—figure supplements 57). Unsurprisingly, the regional-urbanicity model (with either the exponential or power decay) produced some of the most accurate model estimates ( ± 10% of the observed trip counts) for most trip types in Namibia (11/14 trip categories), Kenya (7/14 trip categories), and Burkina Faso (6/9 trip categories) (Supplementary file 3A); however, the regional model produced the most accurate estimates for a range of trips in Zambia (5/13). The proportion of trips estimated within the 10% margin of error was lowest for Kenya, with the most accurate models only well estimating 3–7% of trips, and highest for Namibia, with some models estimating >30% of certain trip types. Depending on the type of trip and country, simpler models may provide more accurate trip estimates than the more complex regional-urbanicity model. We found that just accounting for region may be sufficient for best capturing rural-to-rural as well as general or specific inter-regional travel in Namibia. Further, a wide range of trip types was well estimated by the radiation model in Kenya and by a basic model (exponential decay) in Zambia.

Discussion

To date, there have been limited evaluations of mobility patterns in Sub-Saharan African countries and how the geographic, demographic, and economic differences in many LMICs impact the validity of modeling assumptions to estimate travel. By comparing and evaluating a range of gravity models in four Sub-Saharan African countries, we identified clear patterns of travel that were not well approximated using a basic gravity model. By allowing these models to accommodate differences in travel between predominantly rural and/or urban areas and within- versus between-regions, we found that the best fitting model allowed for the greatest flexibility for estimating different trip types. However, there were differences in which model was best able to estimate particular types of travel, the importance of distance and population, and the overall fit of all models to the data by country. Therefore, selecting a model that adjusts for the trip type and context of interest is important for improving estimates.

If a given application is focused on estimating specific trip types (e.g., urban residents traveling to rural areas for work), these findings may be informative for model selection in similar settings. If a country is sparsely populated with a single predominantly urban district (e.g., capital district), like Burkina Faso, then the regional-urbanicity model would be recommended, as it maximizes the proportion of trips reasonably estimated for most trip types. If specifically considering trips between predominantly rural and urban locations in a sparsely populated country with a few predominantly urban districts, like Namibia, then the regional model would provide the best estimates in these settings. Alternatively, if a country is more homogeneously populated, like Kenya, the radiation model may provide the best estimates for most specific trip types. If data are limited or unavailable, then the parameters fit to these four exemplar countries may serve as a proxy.

Moving forward, the need to select a single model could be mitigated and estimates of mobility patterns could be improved by developing an ensemble model. Recently, an ensemble model outperformed individual models in estimating human mobility patterns in Australia, combining different mobility models as well as data types to optimize mobility estimates across different spatial scales (McCulloch et al., 2021). Using this approach, additional individual level information (e.g., gender, age, occupation) collected from other data sources could be incorporated to study their impact on model fit. Future work could also evaluate the models tested here in other countries, both within and outside of Sub-Saharan Africa. Model estimates of travel patterns in high-income countries may also benefit from accounting for urbanicity and regionality (Truscott and Ferguson, 2012; Xia et al., 2004). The increasing availability of mobility data has enabled the comparison of global human mobility patterns and revealed that (a) longer distance trip ( > 20 km) patterns are similar across low- and high-income settings (Kraemer et al., 2020) and (b) that mobility patterns in rural (sparsely populated) areas differ from those in urban (densely populated) areas (Liu et al., 2015). This suggests that our findings based on trips aggregated to the region or district level will be generalizable to other countries from a range of income settings. However, mobility patterns at smaller spatial scales appear to be different for low- and high-income settings (Kraemer et al., 2020), and requires further investigation.

There are a number of caveats to be noted in this study. Mobile phone data has inherent owner, user, and coverage biases, especially in LMICs where mobility data tend to be concentrated around urban areas and roads (Kraemer et al., 2020; Wesolowski et al., 2013). Furthermore, the dataset from Burkina Faso was limited to 100,000 randomly selected subscribers (1.4% of Burkina Faso’s Telecel Faso subscribers). This smaller sample size may have resulted in the models overestimating routes that were missed in the actual dataset due to relatively small trip counts. Regardless, mobile phone data remain one of the most direct ways to gather information on a large population in LMICs. Alternative human mobility data sources that have recently become more available, such as Google Mobility or Facebook Data for Good (Kissler et al., 2020; Ojal, 2020), could verify these observations in future studies; however, smart phones are also associated with ownership biases and are currently less pervasive than general mobile phone ownership (GSM Association, 2020). It should be noted that models parameterized by different mobility data sources have produced different outcomes (Panigutti et al., 2017). Thus, the generalizability of the results ascertained from these four LMICs should continue to be evaluated and expanded upon as more mobility datasets from other settings become available. Similarly, the ability of other spatial interaction models to better estimate human mobility should be explored. We opted to focus on the gravity model because it has been the most commonly used spatial interactions model; however, one of its shortfalls is that it does not address the competition or synergism that often occurs between potential destinations (Bjørnstad et al., 2019). While this is addressed by the radiation model, the degree to which accounting for region and urbanicity can improve its estimates remains to be considered in the future (Bjørnstad et al., 2019). Another shortfall of the gravity and radiation models is their inability to capture the temporally dynamic nature of human mobility. Here, we focused on one definition and measurement of human mobility patterns; however, seasonal mobility patterns, such as those driven by seasonal work or societal factors (Buckee et al., 2017; Wesolowski et al., 2015a), and irregular migration, such as those driven by political or natural crises and environmental changes (International Organization for Migration, 2019), are ubiquitous. Mobile phone data could also be used to quantify these movements and improve estimates of population distribution and movement over time (Facebook data for good, 2021; Finger et al., 2016; Wesolowski et al., 2015a). Finally, the modifiable areal unit problem (MAUP) is a source of bias in the spatial distribution and aggregation of both CDRs (and other movement data) and population that impacts the results of models dependent on these inputs (Fotheringham and Wong, 2016). Although the general trends in model fits were preserved across a range of administrative levels, the generalization of the results to movement aggregated at finer spatial scales or for different geographical boundaries remains to be determined.

Ultimately, gaining a more complete understanding of travel in diverse geographies will help inform applications in the health, economic, social, and transportation sciences. While incorporating urbanicity and region did improve the gravity model fit to various degrees for different countries, this study highlights the need to continue honing a model framework that can better capture mobility patterns and behavioral nuances in LMICs.

Materials and methods

Population, urbanicity, and geolocation data

Request a detailed protocol

WorldPop population data for each country were analyzed as people per pixel for each district (https://www.worldpop.org/) (Supplementary file 1A). WorldPop gridded building pattern datasets were used to categorize grid cells of each district as urban or rural as described elsewhere ( Dooley et al., 2020) using QGIS v3.6. Districts with more or less than 50% urban grid cells were categorized as urban and rural, respectively. A sensitivity analysis in which urban thresholds of 10% and 50% urban grid cells were compared showed that, while the general trends in model fits remained the same, the overall model fits were worse for the lower threshold in all countries but Zambia (Supplementary file 2C). Trip distances were defined as the haversine distance between centroids of districts. Shapefiles for the different countries were downloaded from DIVA-GIS (https://www.diva-gis.org/). Generally, the district identification numbers were assigned by the map source, with district IDs being clustered within their respective region.

Mobile phone data

Request a detailed protocol

Anonymized call data records (CDRs) were provided by the leading mobile phone provider in each country (Supplementary file 1A). Two districts in Namibia did not have data (Oshakati and Uuvudhiya) and were excluded from analysis. The Burkina Faso provider shared CDRs from a subset of randomly selected subscribers (100,000, ~ 1.4% of subscribers), as opposed to the other countries’ providers that shared CDRs from all of their subscribers. The duration of and year(s) covered by the CDR datasets varied by country, defined by different data sharing agreements: Namibia’s ran October 2, 2010 – April 30, 2014; Kenya’s ran June 1, 2008 – July 3, 2009 (excluding the month of February, 2009); Burkina Faso’s ran January 1, 2016 – December 31, 2016; and Zambia’s ran August 1, 2020 – December 30, 2020. Briefly, CDRs for each country were first aggregated to tower locations and then to districts for each country. A similar method described elsewhere (Zu Erbach-Schoenberg et al., 2016) was used to assign cell towers to districts. Briefly, if a cell tower’s coverage zone fell squarely within one district, all CDRs associated with that tower were assigned to that district. If the coverage zone spanned more than one district, the number of CDRs assigned to each district was split according to the area of overlap between the coverage zone and districts. We only considered travel that crossed district boundaries, not local movement within the district. The average number of total monthly trips taken between each origin and destination was calculated and the proportion of trips was calculated for each origin by normalizing the trip counts to a given destination by the total trips made from that origin. Trip types were defined by origin and destination, either by urbanicity (urban or rural) or by region (intra- or inter-regional). Statistical and spatial analysis was done in R v3.6.3.

Mobility models

Request a detailed protocol

We compared the ability of eight variations of the gravity model and a basic radiation model to capture the heterogeneity in trip counts (Tij) between each origin (i) and destination (j).

Gravity models

Request a detailed protocol

The gravity model estimates the trip counts, T^i,j, as a function of the population sizes at the origin (Pi) and destination (Pj) and deterrence function that depends on the distance between the two locations (di,j) (Equation 1).

(1) T^i,j=θPiαPjβf(di,j)

Here, α and β are non-negative parameters that scale the strength of association between i and j; θ acts as a proportionality constant, and fdij is the penalty associated with a trip distance (d, in kilometers). Both the power fdij=dijγ and exponential fdij=expdi,jD forms of the deterrence function were tested, where γ is a non-negative parameter that determines the rate at which the number of trips decays with trip distance and D is a non-negative parameter that captures the deterrence distance (Chen, 2015). The number of trips increases with larger values of α and β and smaller values of γ or D. We tested eight model variations of Equation 1 in which these parameters were allowed to vary according to aspects of the origin and destination.

Power variants:

Basic: parameters are fitted to the full set of trips,

(2) T^i,j=θPiαPjβdi,jγ

Urbanicity: parameters are fitted to trips categorized by the urbanicity of the origin and destination (k = 1 : 4 for rural-rural, rural-urban, urban-rural, and urban-urban),

(3) T^i,j=θPiαkPjβkdi,jγk{k=1ifurbanicityi=ruralurbanicityj=ruralk=2ifurbanicityi=ruralurbanicityj=urbank=3ifurbanicityi=urbanurbanicityj=ruralk=4ifurbanicityi=urbanurbanicityj=urban

Regional: parameters are fitted based on trips categorized by whether the origin and destination of a trip were both in the same region (intra-regional) or in different regions (inter-regional) (m = 1 : 2),

(4) T^i,j=θPiαmPjβmdi,jγm{m=1ifregioni=regionjm=2ifregioniregionj

Regional-Urbanicity: parameters are fitted to trips categorized by both the region and urbanicity of the origin and destination (n = 1 : 8 for intra-regional- rural-to-rural, inter-regional-rural-to-rural, etc.).

(5) T^i,j=θPiαnPjβndi,jγn{n=1ifregioni=regionjurbanicityi=ruralurbanicityj=ruraln=2ifregioniregionjurbanicityi=ruralurbanicityj=urbann=3ifregioni=regionjurbanicityi=urbanurbanicityj=ruraln=4ifregioniregionjurbanicityi=urbanurbanicityj=urbann=5ifregioni=regionjurbanicityi=ruralurbanicityj=ruraln=6ifregioniregionjurbanicityi=ruralurbanicityj=urbann=7ifregioni=regionjurbanicityi=urbanurbanicityj=ruraln=8ifregioniregionjurbanicityi=urbanurbanicityj=urban

Exponential variants:

Basic: parameters are fitted to the full set of trips,

(6) T^i,j=θPiαPjβexp(di,jD)

Urbanicity: parameters are fitted to trips categorized by the urbanicity of the origin and destination (k = 1 : 4 for rural – rural, rural – urban, urban – rural, and urban – urban. See definitions in Equation 3),

(7) T^i,j=θPiαkPjβkexp(di,jDk)

Regional: parameters are fitted based on trips categorized by whether the origin and destination of a trip were both in the same region (intra-regional) or in different regions (inter-regional) (m = 1 : 2. See definitions in Equation 4),

(8) T^i,j=θPiαmPjβmexp(di,jDm)

Regional-Urbanicity: parameters are fitted to trips categorized by both the region and urbanicity of the origin and destination (n = 1 : 8 for intra-regional- rural-to-rural, inter-regional-rural-to-rural, etc. See definitions in Equation 5),

(9) T^i,j=θPiαnPjβnexp(di,jDn)

We fit the gravity model parameters θ,α,β,γ,andD to observed trip counts extracted from mobile phone data (mi,j) using Bayesian inference, where the model likelihood was assumed to have Poisson error structure (Equation 10) and parameters were given uninformative Gamma priors. The gravity models were fitted to call data records using the R package ‘mobility’, which employs the JAGS (Just Another Gibbs Sampler) Bayesian MCMC algorithm and ‘rjags’ R package (found at https://github.com/COVID-19-Mobility-Data-Network/mobility, John, 2021). The posterior parameter estimates were then used to simulate human mobility patterns.

(10) mi,j Pois(T^i,j)

Radiation model

Request a detailed protocol

Like the gravity model, the radiation model estimates a trip count, T^i,j, as a function of origin and destination population size (Equation 11); however, it differs in the way that it assumes that the probability of making a trip is also influenced by nearby potential destinations (Simini et al., 2012). Thus, T^i,j is also dependent on the total population in the circle (si,j) centered at i with a radius equal to di,j, excluding populations in i and j, and we defined this value by summing the population sizes of districts that fell completely or partially within the radius. The number of trips emanating from origin i is Ti=σPi, where σ is the proportion of the entire country’s population that traveled over a given time period (P=iPi). We fit the parameters σ associated with each Ti,j to trips calculated from mobile phone data (mi,j) using a Poisson error structure (Equation 10). While additional forms of the radiation model have been explored elsewhere (Bjørnstad et al., 2019), we focused on a form that normalizes Ti,j for a finite system (Masucci et al., 2013).

(11) T^i,j=σPi1PiPPiPj(Pi+si,j)(Pi+Pj+si,j)

Model comparisons

Request a detailed protocol

Models were compared using the Deviance Information Criterion (DIC), a criterion designed for MCMC outputs that assesses a model’s trade-off between goodness of fit and complexity (Shriner and Yi, 2009; Spiegelhalter et al., 2002). Models fit to the same datasets (e.g., from the same country) were compared and those with the lowest DIC were selected as the best model. To determine the distribution of trips that were over- or underestimated for a given model, the ratio of estimated to observed trip counts for each route was calculated. Given that the distribution of ratios ranged nine orders of magnitude, the general accuracy of model estimates for specific trip types was evaluated by comparing the proportion of trips with model estimates that fell within ±10% of the observed trips.

Modifiable areal unit problem

Request a detailed protocol

To explore the impact of modifiable areal unit problem (MAUP) (e.g., the effect of the arbitrary definition of administrative units the spatial distribution of CDRs and population), we fit and compared models for a range of administrative units. We reran the gravity models with a power decay function at the administrative one unit (region) level for all countries and administrative three unit level for Burkina Faso, the only country whose dataset was supplied at the administrative three unit level. Note that the models involving regionality could not be run at the administrative one unit. Regardless of the administrative unit level used (e.g., smallest or largest district sizes), the general trend in model ranking was preserved (Supplementary file 2B).

Data availability

Request a detailed protocol

A different form of the datasets from Kenya and Namibia that were negotiated in a prior negation are available as supplements of (Ruktanonchai et al., 2016 and Wesolowski et al., 2015b). Individuals interested in the dataset from Zambia may contact the authors with requests.

Data availability

Due to data sharing agreements with the mobile phone companies, the call data records used in this study are not directly available. However, the Burkina Faso data sharing agreement allows for a jittered dataset of monthly aggregated trips between administrative 2 units, which can be found on Dryad as "Burkina Faso mobility data with some noise" (https://doi.org/10.5061/dryad.fn2z34tt6). A different form of the datasets from Kenya and Namibia that were negotiated in a prior negation are available as supplements of (Ruktanonchai et al., 2016) and (Wesolowski et al., 2015b). Individuals interested in the dataset from Zambia may contact the authors with requests. The code used to analyze the mobile phone data and run the models can be found on github at https://github.com/hrmeredith12/Rural-mobility-models.git (copy archived at https://archive.softwareheritage.org/swh:1:rev:cfc77221c574dad23c3204cd6c5d5fadcb1ce385).

The following data sets were generated
    1. Meredith HR
    (2021) Dryad Digital Repository
    Burkina Faso mobility datawith some noise.
    https://doi.org/10.5061/dryad.fn2z34tt6

References

Article and author information

Author details

  1. Hannah R Meredith

    Department of Epidemiology, Johns Hopkins Bloomberg School of Public Health, Baltimore, United States
    Contribution
    Conceptualization, Formal analysis, Investigation, Methodology, Writing - original draft, Writing – review and editing
    For correspondence
    hmeredi4@jhu.edu
    Competing interests
    none
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-5315-7568
  2. John R Giles

    Department of Epidemiology, Johns Hopkins Bloomberg School of Public Health, Baltimore, United States
    Contribution
    Software, Writing – review and editing
    Competing interests
    none
  3. Javier Perez-Saez

    Department of Epidemiology, Johns Hopkins Bloomberg School of Public Health, Baltimore, United States
    Contribution
    Data curation, Validation, Writing – review and editing
    Competing interests
    none
  4. Théophile Mande

    Bureau d'Etudes Scientifiques et Techniques - Eau, Energie, Environnement (BEST-3E), Ouagadougou, Burkina Faso
    Contribution
    Resources
    Competing interests
    none
  5. Andrea Rinaldo

    1. Dipartimento di Ingegneria Civile Edile ed Ambientale, Università di Padova, Padova, Italy
    2. Laboratory of Ecohydrology, School of Architecture, Civil and Environmental Engineering, École Polytechnique Fédérale de Lausanne, Lausanne, Switzerland
    Contribution
    Data curation, Funding acquisition, Resources, Writing – review and editing
    Competing interests
    none
  6. Simon Mutembo

    1. Department of International Health, International Vaccine Access Center, Johns Hopkins Bloomberg School of Public Health, Baltimore, United States
    2. Macha Research Trust, Choma, Zambia
    Contribution
    Data curation
    Competing interests
    none
  7. Elliot N Kabalo

    Zambia Information and Communications Technology Authority, Lusaka, Zambia
    Contribution
    Data curation
    Competing interests
    none
  8. Kabondo Makungo

    Zamtel, Lusaka, Zambia
    Contribution
    Data curation
    Competing interests
    none
  9. Caroline O Buckee

    Department of Epidemiology and the Center for Communicable Disease Dynamics, Harvard TH Chan School of Public Health, Boston, United States
    Contribution
    Resources
    Competing interests
    None
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-8386-5899
  10. Andrew J Tatem

    WorldPop, School of Geography and Environmental Science, University of Southampton, Southampton, United Kingdom
    Contribution
    Resources, Writing – review and editing
    Competing interests
    none
  11. C Jessica E Metcalf

    Department of Ecology and Evolutionary Biology and the Princeton School of Public and International Affairs, Princeton University, Princeton, United States
    Contribution
    Resources, Writing – review and editing
    Competing interests
    none
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-3166-7521
  12. Amy Wesolowski

    Department of Epidemiology, Johns Hopkins Bloomberg School of Public Health, Baltimore, United States
    Contribution
    Conceptualization, Data curation, Funding acquisition, Resources, Supervision, Writing – review and editing
    For correspondence
    awesolowski@jhu.edu
    Competing interests
    Reviewing editor, eLife
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-6320-3575

Funding

National Institutes of Health (DP2LM013102)

  • Amy Wesolowski

Burroughs Wellcome Fund

  • Amy Wesolowski

Swiss Agency for Development and Cooperation

  • Andrea Rinaldo

Swiss National Science Foundation (200021-172578)

  • Andrea Rinaldo

National Institutes of Health (1R01Al160780-01)

  • Amy Wesolowski

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

Research reported in this publication was supported in part by the National Library Of Medicine of the National Institutes of Health under Award Number DP2LM013102 (HRM, JRG, APW) and 1R01Al160780-01 (APW), a Career Award at the Scientific Interface by the Burroughs Wellcome Fund (HRM, JRG, APW), the Swiss Agency for Development and Cooperation - 2iE partie scientifique – through “Projet 3E Afrique, Burkina Faso” (JSP, TM, AR), and a Swiss National Science Foundation grant under Award Number 200021–172578 (JSP, AR). The content is solely the responsibility of the authors and does not necessarily represent the official views of the funders.

Copyright

© 2021, Meredith et al.

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 1,896
    views
  • 207
    downloads
  • 11
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

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. Hannah R Meredith
  2. John R Giles
  3. Javier Perez-Saez
  4. Théophile Mande
  5. Andrea Rinaldo
  6. Simon Mutembo
  7. Elliot N Kabalo
  8. Kabondo Makungo
  9. Caroline O Buckee
  10. Andrew J Tatem
  11. C Jessica E Metcalf
  12. Amy Wesolowski
(2021)
Characterizing human mobility patterns in rural settings of sub-Saharan Africa
eLife 10:e68441.
https://doi.org/10.7554/eLife.68441

Share this article

https://doi.org/10.7554/eLife.68441

Further reading

    1. Computational and Systems Biology
    2. Microbiology and Infectious Disease
    Gaetan De Waele, Gerben Menschaert, Willem Waegeman
    Research Article

    Timely and effective use of antimicrobial drugs can improve patient outcomes, as well as help safeguard against resistance development. Matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS) is currently routinely used in clinical diagnostics for rapid species identification. Mining additional data from said spectra in the form of antimicrobial resistance (AMR) profiles is, therefore, highly promising. Such AMR profiles could serve as a drop-in solution for drastically improving treatment efficiency, effectiveness, and costs. This study endeavors to develop the first machine learning models capable of predicting AMR profiles for the whole repertoire of species and drugs encountered in clinical microbiology. The resulting models can be interpreted as drug recommender systems for infectious diseases. We find that our dual-branch method delivers considerably higher performance compared to previous approaches. In addition, experiments show that the models can be efficiently fine-tuned to data from other clinical laboratories. MALDI-TOF-based AMR recommender systems can, hence, greatly extend the value of MALDI-TOF MS for clinical diagnostics. All code supporting this study is distributed on PyPI and is packaged at https://github.com/gdewael/maldi-nn.

    1. Computational and Systems Biology
    2. Genetics and Genomics
    Sanjarbek Hudaiberdiev, Ivan Ovcharenko
    Research Article

    Enhancers and promoters are classically considered to be bound by a small set of transcription factors (TFs) in a sequence-specific manner. This assumption has come under increasing skepticism as the datasets of ChIP-seq assays of TFs have expanded. In particular, high-occupancy target (HOT) loci attract hundreds of TFs with often no detectable correlation between ChIP-seq peaks and DNA-binding motif presence. Here, we used a set of 1003 TF ChIP-seq datasets (HepG2, K562, H1) to analyze the patterns of ChIP-seq peak co-occurrence in combination with functional genomics datasets. We identified 43,891 HOT loci forming at the promoter (53%) and enhancer (47%) regions. HOT promoters regulate housekeeping genes, whereas HOT enhancers are involved in tissue-specific process regulation. HOT loci form the foundation of human super-enhancers and evolve under strong negative selection, with some of these loci being located in ultraconserved regions. Sequence-based classification analysis of HOT loci suggested that their formation is driven by the sequence features, and the density of mapped ChIP-seq peaks across TF-bound loci correlates with sequence features and the expression level of flanking genes. Based on the affinities to bind to promoters and enhancers we detected five distinct clusters of TFs that form the core of the HOT loci. We report an abundance of HOT loci in the human genome and a commitment of 51% of all TF ChIP-seq binding events to HOT locus formation thus challenging the classical model of enhancer activity and propose a model of HOT locus formation based on the existence of large transcriptional condensates.