Information

  • Author Services

Initiatives

You are accessing a machine-readable page. In order to be human-readable, please install an RSS reader.

All articles published by MDPI are made immediately available worldwide under an open access license. No special permission is required to reuse all or part of the article published by MDPI, including figures and tables. For articles published under an open access Creative Common CC BY license, any part of the article may be reused without permission provided that the original article is clearly cited. For more information, please refer to https://www.mdpi.com/openaccess .

Feature papers represent the most advanced research with significant potential for high impact in the field. A Feature Paper should be a substantial original Article that involves several techniques or approaches, provides an outlook for future research directions and describes possible research applications.

Feature papers are submitted upon individual invitation or recommendation by the scientific editors and must receive positive feedback from the reviewers.

Editor’s Choice articles are based on recommendations by the scientific editors of MDPI journals from around the world. Editors select a small number of articles recently published in the journal that they believe will be particularly interesting to readers, or important in the respective research area. The aim is to provide a snapshot of some of the most exciting work published in the various research areas of the journal.

Original Submission Date Received: .

  • Active Journals
  • Find a Journal
  • Proceedings Series
  • For Authors
  • For Reviewers
  • For Editors
  • For Librarians
  • For Publishers
  • For Societies
  • For Conference Organizers
  • Open Access Policy
  • Institutional Open Access Program
  • Special Issues Guidelines
  • Editorial Process
  • Research and Publication Ethics
  • Article Processing Charges
  • Testimonials
  • Preprints.org
  • SciProfiles
  • Encyclopedia

ijerph-logo

Article Menu

botswana hiv aids case study geography

  • Subscribe SciFeed
  • Recommended Articles
  • PubMed/Medline
  • Google Scholar
  • on Google Scholar
  • Table of Contents

Find support for a specific problem in the support section of our website.

Please let us know what you think of our products and services.

Visit our dedicated information section to learn more about MDPI.

JSmol Viewer

Spatial analysis of hiv infection and associated risk factors in botswana.

botswana hiv aids case study geography

1. Introduction

2. materials and methods, 2.1. study area, 2.2. study data, 2.3. measures, 2.4. statistical analysis, 2.5. spatial analyses, 3.1. socio-demographic and behavioural characteristics of participants, 3.2. distribution of hiv cases, 3.3. identification of spatial clusters in botswana, 4. discussion, strengths and limitations, 5. conclusions, author contributions, institutional review board statement, informed consent statement, data availability statement, acknowledgments, conflicts of interest.

  • Piot, P.; Plummer, F.A.; Mhalu, F.S.; Lamboray, J.L.; Chin, J.; Mann, J.M. AIDS: An international perspective. Science 1988 , 239 , 573–579. [ Google Scholar ] [ CrossRef ] [ PubMed ] [ Green Version ]
  • Barankanira, E.; Molinari, N.; Niyongabo, T.; Laurent, C. Spatial analysis of HIV infection and associated individual characteristics in Burundi: Indications for effective prevention. BMC Public Health 2015 , 16 , 118. [ Google Scholar ] [ CrossRef ] [ PubMed ] [ Green Version ]
  • Agyei-Mensah, S. The HIV/AIDS epidemic in sub-Saharan Africa: Homogeneity or heterogeneity? Nor. Geogr. Tidsskr. Nor. J. Geogr. 2005 , 59 , 14–25. [ Google Scholar ] [ CrossRef ]
  • Zulu, L.C.; Kalipeni, E.; Johannes, E. Analyzing spatial clustering and the spatiotemporal nature and trends of HIV/AIDS prevalence using GIS: The case of Malawi, 1994–2010. BMC Infect. Dis. 2014 , 14 , 285. [ Google Scholar ] [ CrossRef ] [ PubMed ] [ Green Version ]
  • Bulstra, C.A.; Hontelez, J.A.C.; Giardina, F.; Steen, R.; Nagelkerke, N.J.D.; Bärnighausen, T.; De Vlas, S.J. Mapping and characterising areas with high levels of HIV transmission in sub-Saharan Africa: A geospatial analysis of national survey data. PLoS Med. 2020 , 17 , e1003042. [ Google Scholar ] [ CrossRef ] [ PubMed ] [ Green Version ]
  • Botswana, R. Third National Multisectoral HIV and AIDS Response Strategic Framework 2018–2023—Enhancing Efficiencies through Integrated Approach ; Wellness, M., Ed.; National AIDS Co-Ordinating Agency (NACA): Gaborone, Botswana, 2018. [ Google Scholar ]
  • Levinsohn, J.A.; Dinkelman, T.; Majelantle, R. When Knowledge is Not Enough: Hiv/Aids Information and Risky Behavior in Botswana. Natl. Bur. Econ. Res. 2006 . [ Google Scholar ] [ CrossRef ]
  • Lang, E.; Menon, V. Financing ARV Treatment in Botswana. Available online: https://www.healthpolicyproject.com/pubs/857_HPPBotswanaBriefindesign.pdf (accessed on 19 March 2021).
  • Sabone, M.; Magowe, M.; Busang, L.; Moalosi, J.; Binagwa, B.; Mwambona, J. Impediments for the Uptake of the Botswana Government’s Male Circumcision Initiative for HIV Prevention. Sci. World J. 2013 , 2013 , 1–7. [ Google Scholar ] [ CrossRef ] [ PubMed ] [ Green Version ]
  • Stover, J.; Fidzani, B.; Molomo, B.C.; Moeti, T.; Musuka, G. Estimated HIV Trends and Program Effects in Botswana. PLoS ONE 2008 , 3 , e3729. [ Google Scholar ] [ CrossRef ] [ Green Version ]
  • Republic of Botswana Ministry of Heath; National AIDS Coordinating Agency; Central Statistics Office. Botswana AIDS Impact Survey IV 2013 ; Central Statistics Office: Gaborone, Botswana, 2013. [ Google Scholar ]
  • Matlho, K.; Randell, M.; Lebelonyane, R.; Kefas, J.; Driscoll, T.; Negin, J. HIV prevalence and related behaviours of older people in Botswana—secondary analysis of the Botswana AIDS Impact Survey (BAIS) IV. Afr. J. AIDS Res. 2019 , 18 , 18–26. [ Google Scholar ] [ CrossRef ] [ Green Version ]
  • (UNAIDS) JUAP. Connecting People Faster to HIV Servicces ; Joint United Nations Programme on HIV/AIDS(UNAIDS): Geneva, Switzerland, 2013. [ Google Scholar ]
  • Ama, N.O.; Shaibu, S.; Burnette, J.D. HIV-related knowledge and practices: A cross-sectional study among adults aged 50 years and above in Botswana. South Afr. Fam. Pr. 2016 , 58 , 100–107. [ Google Scholar ] [ CrossRef ] [ Green Version ]
  • Kandala, N.-B.; Campbell, E.K.; Rakgoasi, S.D.; Madi-Segwagwe, B.C.; Fako, T.T. The geography of HIV/AIDS prevalence rates in Botswana. HIV/AIDS Res. Palliat. Care 2012 , 4 , 95–102. [ Google Scholar ] [ CrossRef ] [ Green Version ]
  • Tanser, F.; Bärnighausen, T.; Cooke, G.S.; Newell, M.-L. Localized spatial clustering of HIV infections in a widely disseminated rural South African epidemic. Int. J. Epidemiol. 2009 , 38 , 1008–1016. [ Google Scholar ] [ CrossRef ] [ PubMed ] [ Green Version ]
  • Chimoyi, L.A.; Musenge, E. Spatial analysis of factors associated with HIV infection among young people in Uganda, 2011. BMC Public Health 2014 , 14 , 555. [ Google Scholar ] [ CrossRef ] [ Green Version ]
  • Dwyer-Lindgren, L.; Cork, M.A.; Sligar, A.; Steuben, K.M.; Wilson, K.F.; Provost, N.R.; Mayala, B.K.; VanderHeide, J.D.; Collison, M.L.; Hall, J.B.; et al. Mapping HIV prevalence in sub-Saharan Africa between 2000 and 2017. Nat. Cell Biol. 2019 , 570 , 189–193. [ Google Scholar ] [ CrossRef ] [ PubMed ] [ Green Version ]
  • USAID. Health Finance Profile: Botswana ; Health Policy Project: Washington, DC, USA, 2016. [ Google Scholar ]
  • Republic of Botswana Ministry of Heath; National AIDS Coordinating Agency; Central Statistics Office. Botswana AIDS Impact Survey IV—Statistical Report ; Central Statistics Office: Gaborone, Botswana, 2013. [ Google Scholar ]
  • World Health Organization. Guidelines for Second Generation HIV Surveillance ; World Health Organisation: Geneva, Switzerland, 2000. [ Google Scholar ]
  • Kulldorff, M. A spatial scan statistic. Commun. Stat. Theory Methods 1997 , 26 , 1481–1496. [ Google Scholar ] [ CrossRef ]
  • Ozdenerol, E.; Williams, B.L.; Kang, S.Y.; Magsumbol, M.S. Comparison of spatial scan statistic and spatial filtering in estimating low birth weight clusters. Int. J. Health Geogr. 2005 , 4 , 19. [ Google Scholar ] [ CrossRef ] [ Green Version ]
  • Kim, S.; Jung, I. Optimizing the maximum reported cluster size in the spatial scan statistic for ordinal data. PLoS ONE 2017 , 12 , e0182234. [ Google Scholar ] [ CrossRef ] [ Green Version ]
  • Carrel, M.; Janko, M.; Mwandagalirwa, M.K.; Morgan, C.; Fwamba, F.; Muwonga, J.; Tshefu, A.K.; Meshnick, S.; Emch, M. Changing spatial patterns and increasing rurality of HIV prevalence in the Democratic Republic of the Congo between 2007 and 2013. Health Place 2016 , 39 , 79–85. [ Google Scholar ] [ CrossRef ] [ PubMed ] [ Green Version ]
  • Oseni, Z.; Seedat, F.; Kandala, N.-B. Hiv epidemic heterogeneity in zimbabwe: Evidence from successive demographic and health surveys. J. Biosoc. Sci. 2018 , 50 , 840–852. [ Google Scholar ] [ CrossRef ] [ PubMed ]
  • Steen, R.; Hontelez, J.A.; Mugurungi, O.; Mpofu, A.; Matthijsse, S.M.; De Vlas, S.J.; Dallabetta, G.A.; Cowan, F.M. Economy, migrant labour and sex work: Interplay of HIV epidemic drivers in Zimbabwe over three decades. AIDS 2019 , 33 , 123–131. [ Google Scholar ] [ CrossRef ]
  • Nicholas, P.K.; Mfono, N.; Corless, I.B.; Davis, S.M.; O’Brien, E.; Padua, J.; Ahmed, S.; Sheldon, L.K.; Oas, K.; Sadler, L.; et al. HIV vulnerability in migrant populations in southern Africa: Sociological, cultural, health-related, and human-rights perspectives. Int. J. Afr. Nurs. Sci. 2016 , 5 , 1–8. [ Google Scholar ] [ CrossRef ] [ Green Version ]
  • Uganda, K. Most at Risk Populations–Long Distance Truck Drivers and HIV/AIDS in Uganda: Synthesis of Information and Evidence to Inform the Response. Synth. Rep. 2014 . [ Google Scholar ]
  • Tanser, F.; Bärnighausen, T.; Dobra, A.; Sartorius, B. Identifying ‘corridors of HIV transmission’ in a severely affected rural South African population: A case for a shift toward targeted prevention strategies. Int. J. Epidemiol. 2018 , 47 , 537–549. [ Google Scholar ] [ CrossRef ] [ Green Version ]
  • Institute, K. HIV in Australia: Annual Surveillance Short Report 2018 ; Unsw Sydney, Kirby Institute: Sydney, Australia, 2018. [ Google Scholar ]
  • Burgoyne, A.D.; Drummond, P.D. Knowledge of HIV and AIDS in women in sub-Saharan Africa. Afr. J. Reprod. Health 2008 , 12 , 1–14. [ Google Scholar ]
  • Su, Y.Y. The failure of the American ABC HIV prevention model in Botswana. SURG J. 2010 , 4 , 93–100. [ Google Scholar ] [ CrossRef ]
  • Katisi, M.; Daniel, M. Safe male circumcision in Botswana: Tension between traditional practices and biomedical marketing. Glob. Public Health 2015 , 10 , 739–756. [ Google Scholar ] [ CrossRef ] [ Green Version ]
  • Osuafor, G.N.; Mturi, A.J. Attitude towards sexual control among women in conjugal union in the era of the HIV/AIDS epidemic in Mahikeng, South Africa. Afr. Popul. Stud. 2014 , 28 , 538–550. [ Google Scholar ] [ CrossRef ] [ Green Version ]
  • Mabaso, M.; Makola, L.; Naidoo, I.; Mlangeni, L.L.; Jooste, S.; Simbayi, L. HIV prevalence in South Africa through gender and racial lenses: Results from the 2012 population-based national household survey. Int. J. Equity Health 2019 , 18 , 1–11. [ Google Scholar ] [ CrossRef ] [ Green Version ]
  • UNAIDS. Ending AIDS: Progress towards 90-90-90 Targets ; UN Joint Programme on HIV/AIDS (UNAIDS): Geneva, Switzerland, 2017. [ Google Scholar ]
  • Miller, K.S.; Cham, H.J.; Taylor, E.M.; Berrier, F.L.; Duffy, M.; Vig, J.; Chipazi, L.; Chakalisa, C.; Sidibe, S.; Swart, K.; et al. Formative Work and Community Engagement Approaches for Implementing an HIV Intervention in Botswana Schools. Am. J. Public Health 2016 , 106 , 1439–1441. [ Google Scholar ] [ CrossRef ]
  • Muyunda, B.; Musonda, P.; Mee, P.; Todd, J.; Michelo, C. Educational Attainment as a Predictor of HIV Testing Uptake Among Women of Child-Bearing Age: Analysis of 2014 Demographic and Health Survey in Zambia. Front. Public Health 2018 , 6 , 192. [ Google Scholar ] [ CrossRef ] [ Green Version ]
  • Letshwenyo-Maruatona, S.B.; Madisa, M.; Boitshwarelo, T.; George-Kefilwe, B.; Kingori, C.; Ice, G.; A Bianco, J.; Marape, M.; Haile, Z.T. Association between HIV/AIDS knowledge and stigma towards people living with HIV/AIDS in Botswana. Afr. J. AIDS Res. 2019 , 18 , 58–64. [ Google Scholar ] [ CrossRef ] [ PubMed ]
  • Nutor, J.J.; Duah, H.O.; Agbadi, P.; Duodu, P.A.; Gondwe, K.W. Spatial analysis of factors associated with HIV infection in Malawi: Indicators for effective prevention. BMC Public Health 2020 , 20 , 1–14. [ Google Scholar ] [ CrossRef ] [ PubMed ]
  • Phillips, D.L.; Clancy, K.J. Some Effects of "Social Desirability" in Survey Studies. Am. J. Sociol. 1972 , 77 , 921–940. [ Google Scholar ] [ CrossRef ]
  • Andersson, N.; Cockcroft, A. Male circumcision, attitudes to HIV prevention and HIV status: A cross-sectional study in Botswana, Namibia and Swaziland. AIDS Care 2011 , 24 , 301–309. [ Google Scholar ] [ CrossRef ] [ PubMed ]

Click here to enlarge figure

Variablen (Weight %)Proportion of HIV+p-Value
Gender
Female2653 (56.4)15.6<0.001
Male2055 (43.6)9.5
Age group <0.001
15–241534 (32.6)2.5
25–341275 (27.1)7.7
35–44902 (19.1)8.4
45–54602 (12.8)4.7
55+395 (8.4)1.8
Place of residence 0.6065
Urban3040 (64.6)16.0
Rural1668 (35.4)9.1
Education level <0.001
None465 (9.9)2.8
Primary950 (20.2)7.5
Secondary2440 (51.8)12.2
Tertiary853 (18.1)2.6
Marital status <0.001
Married2066 (43.9)13.4
Never married2642 (56.1)11.7
Religion 0.747
Christianity3966 (84.2)21.1
Other742(15.8)4.0
Employed
Yes2525 (53.6)15.7<0.001
No2183 (46.6)9.4
Alcohol use
Yes1643 (34.9)8.60.584
No3065 (65.1)16.5
Condom use
Yes1991(42.3)12.9<0.001
No2717(57.7)12.3
VariablesUnivariate RegressionMultivariable Regression
OR CI *p-ValueOR CIp-Value
Gender
MaleRef Ref
Female1.381.15–1.650.0011.421.16–1.730.001
Age group
15–24Ref Ref
25–344.763.51–6.45<0.0015.043.60–7.05<0.001
35–449.296.79–12.7<0.0019.576.61–13.9<0.001
45–546.974.96–9.81<0.0017.074.59–10.9<0.001
55+3.362.25–5.03<0.0013.582.20–5.83<0.001
Marital status
Never MarriedRef Ref
Married0.600.50–0.71<0.0011.020.83–1.270.831
Education
NoneRef Ref
Primary1.521.10–2.110.0111.340.95–1.890.094
Secondary0.800.59–1.090.1611.110.77–1.600.56
Tertiary0.430.29–0.65<0.0010.420.27–0.66<0.001
Place of Residence
RuralRef Ref
Urban0.950.80–1.40.584---
Religion
ChristianRef Ref
Other1.040.81–1.330.749---
Employment
NoRef Ref
Yes0.620.52–0.74<0.0011.050.85–1.300.664
Alcohol use
NoRef Ref
Yes0.950.79–1.150.584---
Condom use
NoRef Ref
Yes1.621.36–1.94<0.0011.561.28–1.91<0.001
VariableInside the Clusters (%)Outside the Clusters (%)p-Value
HIV positive30240.012
Gender 0.222
Female5457
Male4643
Age group 0.559
15–243033
25–342727
35–441919
45–541513
55+98
Place of residence <0.001
Urban5268
Rural4832
Education level 0.001
None139
Primary2619
Secondary5152
Tertiary1020
Marital status 0.885
Married4444
Never married5656
Religion <0.001
Christianity7686
Other2413
Employed 0.002
Yes4855
No5244
Alcohol use 0.779
Yes3535
No6565
Condom use <0.312
Yes4442
No5558
MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

Solomon, M.; Furuya-Kanamori, L.; Wangdi, K. Spatial Analysis of HIV Infection and Associated Risk Factors in Botswana. Int. J. Environ. Res. Public Health 2021 , 18 , 3424. https://doi.org/10.3390/ijerph18073424

Solomon M, Furuya-Kanamori L, Wangdi K. Spatial Analysis of HIV Infection and Associated Risk Factors in Botswana. International Journal of Environmental Research and Public Health . 2021; 18(7):3424. https://doi.org/10.3390/ijerph18073424

Solomon, Malebogo, Luis Furuya-Kanamori, and Kinley Wangdi. 2021. "Spatial Analysis of HIV Infection and Associated Risk Factors in Botswana" International Journal of Environmental Research and Public Health 18, no. 7: 3424. https://doi.org/10.3390/ijerph18073424

Article Metrics

Article access statistics, further information, mdpi initiatives, follow mdpi.

MDPI

Subscribe to receive issue release notifications and newsletters from MDPI journals

Log in using your username and password

  • Search More Search for this keyword Advanced search
  • Latest content
  • Current issue
  • BMJ Journals

You are here

  • Volume 65, Issue Suppl 2
  • The geography of HIV/AIDS infection in Botswana
  • Article Text
  • Article info
  • Citation Tools
  • Rapid Responses
  • Article metrics

Download PDF

  • NB Kandala 1 , 2* ,
  • E Campbell 2 ,
  • D Rakgoasi 2 ,
  • 1 WMS, University of Warwick, Coventry, UK
  • 2 Population Studies, University of Botswana, Gaborone, Botswana
  • 3 HIV and AIDS TB, Malaria Policy Development and Harmonisation, SADC, Gaborone, Botswana

Background Botswana's Human immunodeficiency virus (HIV)/Acquired immune deficiency syndrome (AIDS) epidemic is unprecedented in magnitude and impact. For two decades starting during the early 1990s, HIV infection rates in this country of less than two million people have grown tremendously. While the impact of HIV/AIDS is clear for all to see, for a number of years, estimating the population based HIV prevalence rate was a big challenge due to absence of data. Most estimates of HIV prevalence relied heavily on estimates derived from sentinel surveillance of pregnant women attending antenatal care.

Methods Approximately 15,000 respondents (54.0% female) were asked to give blood for syphilis and HIV testing for the 2008 Botswana AIDS Impact Survey. Samples for HIV testing were dried blood spots on a filter paper card taken from a venous blood specimen. A three-stage testing procedure was used with 10% of the negative samples retested and discordant results tested by Western Blot. A Bayesian geo-additive mixed model based on Markov Chain Monte Carlo techniques was used to map the geographic distribution of HIV/AIDS prevalence at the 26 districts, accounting for important risk factors.

Findings The overall HIV prevalence was 17.6% in 2008, higher prevalence among females and in cities and town but lower prevalence among professionals. The mean age for men was lower compared to their female counterpart (25.1 years (SD: 19.3 vs 27.2 years (SD: 20.5). We observed a U-shape association between age and the prevalence of HIV. Unadjusted/adjusted Odds Ratios (OR) indicate that the highest HIV prevalence was in Selebi-Phikwe (OR and 95% CI: 3.29 (2.17 to 4.96)), Sowa (OR and 95% CI: 2.87 (1.51 to 5.49)), Francistown (OR and 95% CI: 2.75 (1.83 to 4.12)) followed by Chobe, Northeast, Ngamiland South, Central-Serowe, Central-Tutume, Central-Bobonong, Kgalagadi South, Orapa, Central-Mahalapye, Ngamiland North, Gaborone, Lobatse, Jwaneng, Ngwaketse West, Kweneng East, Central-Boteti, Kgatleng, Southern, Barolong, Ghanzi, Southeast districts, with the lowest prevalence in Kgalagadi North, Kweneng West districts.

Policy Implication Based on large population cross-sectional household survey, this study shows a clear geographic distribution of the HIV/AIDS epidemic in Botswana with highest incidence of infection is in the east-central districts. Generally, it is most prevalent in the northern and north-eastern districts of the country. HIV prevalence is also quite high in the central districts. The prevalence rate is moderate in the southern districts. It is generally believed that the geographical distribution of the virus is largely explained by trucking route from countries north of Botswana to South Africa and mine districts.

https://doi.org/10.1136/jech.2011.143586.78

Statistics from Altmetric.com

Request permissions.

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

Read the full text or download the PDF:

Javascript is currently disabled in your browser. Several features of this site will not function whilst javascript is disabled.

  • Why Publish With Us?
  • Editorial Policies
  • Author Guidelines
  • Peer Review Guidelines
  • Open Outlook
  • Submit New Manuscript

botswana hiv aids case study geography

  • Sustainability
  • Press Center
  • Testimonials
  • Favored Author Program
  • Permissions
  • Pre-Submission

Chinese website (中文网站)

open access to scientific and medical research

Back to Journals » HIV/AIDS - Research and Palliative Care » Volume 4

botswana hiv aids case study geography

The geography of HIV/AIDS prevalence rates in Botswana

  • Get Permission
  • Cite this article

Authors Kandala N , Campbell , Rakgoasi , Madi , Fako TT

Received 4 February 2012

Accepted for publication 28 March 2012

Published 18 July 2012 Volume 2012:4 Pages 95—102

DOI https://doi.org/10.2147/HIV.S30537

Review by Single anonymous peer review

Peer reviewer comments 2

Ngianga-Bakwin Kandala, 1 Eugene K Campbell, 2 Serai Dan Rakgoasi, 2 Banyana C Madi-Segwagwe, 3 Thabo T Fako 4 1 University of Warwick, Warwick Medical School, Division of Health Sciences; Populations, Evidence and Technologies Group, Warwick Evidence, Coventry, UK; 2 Department of Population Studies, University of Botswana, 3 SADC Secretariat, Directorate of Social and Human Development and Special Programmes, 4 Vice Chancellor's Office, University of Botswana, Gaborone, Botswana Background: Botswana has the second-highest human immunodeficiency virus (HIV) infection rate in the world, with one in three adults infected. However, there is significant geographic variation at the district level and HIV prevalence is heterogeneous with the highest prevalence recorded in Selebi-Phikwe and North East. There is a lack of age-and location-adjusted prevalence maps that could be used for targeting HIV educational programs and efficient allocation of resources to higher risk groups. Methods: We used a nationally representative household survey to investigate and explain district level inequalities in HIV rates. A Bayesian geoadditive mixed model based on Markov Chain Monte Carlo techniques was applied to map the geographic distribution of HIV prevalence in the 26 districts, accounting simultaneously for individual, household, and area factors using the 2008 Botswana HIV Impact Survey. Results: Overall, HIV prevalence was 17.6%, which was higher among females (20.4%) than males (14.3%). HIV prevalence was higher in cities and towns (20.3%) than in urban villages and rural areas (16.6% and 16.9%, respectively). We also observed an inverse U-shape association between age and prevalence of HIV, which had a different pattern in males and females. HIV prevalence was lowest among those aged 24 years or less and HIV affected over a third of those aged 25–35 years, before reaching a peak among the 36–49-year age group, after which the rate of HIV infection decreased by more than half among those aged 50 years and over. In a multivariate analysis, there was a statistically significant higher likelihood of HIV among females compared with males, and in clerical workers compared with professionals. The district-specific net spatial effects of HIV indicated a significantly higher HIV rate of 66% (posterior odds ratio of 1.66) in the northeast districts (Selebi-Phikwe, Sowa, and Francistown) and a reduced rate of 27% (posterior odds ratio of 0.73) in Kgalagadi North and Kweneng West districts. Conclusion: This study showed a clear geographic distribution of the HIV epidemic, with the highest prevalence in the east-central districts. This study provides age- and location-adjusted prevalence maps that could be used for the targeting of HIV educational programs and efficient allocation of resources to higher risk groups. There is need for further research to determine the social, cultural, economic, behavioral, and other distal factors that might explain the high infection rates in some of the high-risk areas in Botswana. Keywords: Botswana, HIV prevalence, geographic location, spatial autocorrelation

Creative Commons License

Contact Us   •   Privacy Policy   •   Associations & Partners   •   Testimonials   •   Terms & Conditions   •   Recommend this site •   Cookies •   Top

Contact Us   •   Privacy Policy

  • DOI: 10.1136/jech.2011.143586.78
  • Corpus ID: 129519154

The geography of HIV/AIDS infection in Botswana

  • N. Kandala , Eugene K. Campbell , +1 author B. Mádi
  • Published in Journal of Epidemiology and… 1 September 2011
  • Geography, Medicine

Related Papers

Showing 1 through 3 of 0 Related Papers

Academia.edu no longer supports Internet Explorer.

To browse Academia.edu and the wider internet faster and more securely, please take a few seconds to  upgrade your browser .

Enter the email address you signed up with and we'll email you a reset link.

  • We're Hiring!
  • Help Center

paper cover thumbnail

The geography of HIV/AIDS prevalence rates in Botswana

Profile image of Serai D Rakgoasi

2012, HIV/AIDS - Research and Palliative Care

Related Papers

International Journal of Environmental Research and Public Health

Malebogo Solomon

Botswana has the third highest human immunodeficiency virus (HIV) prevalence globally, and the severity of the epidemic within the country varies considerably between the districts. This study aimed to identify clusters of HIV and associated factors among adults in Botswana. Data from the Botswana Acquired Immunodeficiency Syndrome (AIDS) Impact Survey IV (BIAS IV), a nationally representative household-based survey, were used for this study. Multivariable logistic regression and Kulldorf’s scan statistics were used to identify the risk factors and HIV clusters. Socio-demographic characteristics were compared within and outside the clusters. HIV prevalence among the study participants was 25.1% (95% CI 23.3–26.4). HIV infection was significantly higher among the female gender, those older than 24 years and those reporting the use of condoms, while tertiary education had a protective effect. Two significant HIV clusters were identified, one located between Selibe-Phikwe and Francisto...

botswana hiv aids case study geography

Khátia Munguambe

sciepub.com SciEP

More than 2/3 of global HIV/AIDS epidemics, as well as seven in ten of all HIV/AIDs death, are reported in Africa. The economic burden of the disease in terms of treatment and prevention for an already impoverished continent is enormous. Several international organizations are engaging Africa to alleviate the disease epidemics and its economic impacts. This study attempted to determine the geographical controls of HIV/AIDS prevalence rates and established a relationship with demographic, social and economic determinants. Accordingly, significant (I = 0.71; α = 0.5) geographical clustering of adult HIV/AIDS prevalence rates were detected, such that the higher rates are clustered in the southern African countries, while the west African countries have the clustering of the lowest rates. The demographic, social and economic factors such as literacy (female) rate, life expectancy at birth, contraceptive prevalence rate, education expenditure and percent Muslim population, significantly explained the variability of the Adult HIV/AIDS prevalence rates (R 2 = 0.81 – 0.84; α = 0.05) in Africa. This result highlights the need for geographically-resolved resource allocation, particularly in the areas of primary education and gender participation for the fight against HIV/AIDS.

This study offers an alternative method rooted in GIS techniques and spatial analysis to estimate HIV/AIDS prevalence over space from an incomplete surveillance data set and explain the variation of those estimates. The results clearly show that the HIV/AIDS epidemic is complex and that it is interconnected with other geographic, historical, economic and cultural phenomena which help explain its spatial spread and variation. The regression models which were developed in this paper illustrated that variables which measure the historical context of colonialism such as resource exploitation and labor migration, gender, culture, contemporary global forces, poverty and disease burden have all contributed variously to the rapid spread of this disease both in space and time. The policy implication is that concentrating on behavior change or therapy alone may not turn the epidemic around. The attack needs to be multifaceted and interdisciplinary taking into consideration the context and the economic and social realities at multiple spatial scales.

BMC Health Services Research

Elizabeth Lutge

Development Policy Review

Francis Owusu

Lemma Merid

Annals of Global Health

Enrique Saldarriaga

Samson Nzou

Godfrey Musuka

Loading Preview

Sorry, preview is currently unavailable. You can download the paper by clicking the button above.

RELATED PAPERS

BMC Infectious Diseases

Natashia Morris

RELATED TOPICS

  •   We're Hiring!
  •   Help Center
  • Find new research papers in:
  • Health Sciences
  • Earth Sciences
  • Cognitive Science
  • Mathematics
  • Computer Science
  • Academia ©2024

Europe PMC requires Javascript to function effectively.

Either your web browser doesn't support Javascript or it is currently turned off. In the latter case, please turn on Javascript support in your web browser and reload this page.

The Geography of HIV/AIDS Prevalence Rates in Botswana

Type Journal Article - HIV/AIDS
Title The Geography of HIV/AIDS Prevalence Rates in Botswana
Author(s)
Volume 4
Publication (Day/Month/Year) 2012
Page numbers 95-102
URL
Abstract
»
»
  • LATEST INFORMATION

Botswana

  • High contrast
  • Press centre

Botswana

Search UNICEF

Botswana hiv/aids impact survey report.

BIAS

The Fifth Botswana AIDS Impact Survey, BAIS V, was a household-based national survey among adults (defined as individuals aged 15 to 64 years) and children (defined as individuals aged 6 weeks to 14 years) conducted from March 2021 to August 2021 to measure the impact of the national HIV response. The survey offered HIV counseling and testing with return of results to the participants and collected information about the uptake of HIV care and treatment services. BAIS V was led by the National AIDS and Health Promotion Agency (NAHPA), the Ministry of Health (MOH), and Statistics Botswana. 

This BAIS V data were used to estimate national HIV incidence, national and district-level HIV prevalence, and viral load suppression (VLS), defined as HIV RNA <1,000 copies per milliliter (mL) among adults living with HIV. The previous BAIS surveys were conducted in 2001, 2005, 2008, and 2013. The results of these five surveys provide critical information on national and district-level progress toward control of the HIV epidemic.

Files available for download

Related topics, more to explore, growing up with shuga, botswana first lady launches u-report hiv testing poll.

U-Report is a messaging tool that empowers young people around the world to engage with and speak out on issues that matter to them.

Children lead #CRC30 reflections

Children take center stage at Convention on the Rights of the Child 30th Anniversary National Commemoration

U.S. flag

An official website of the United States government

The .gov means it’s official. Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

The site is secure. The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

  • Publications
  • Account settings

Preview improvements coming to the PMC website in October 2024. Learn More or Try it out now .

  • Advanced Search
  • Journal List

Logo of medrxiv

This is a preprint.

Population mobility and the development of botswana’s generalized hiv epidemic: a network analysis.

1 Center for Biomedical Modeling, Department of Psychiatry and Biobehavioral Sciences, Semel Institute for Neuroscience and Human Behavior, David Geffen School of Medicine, University of California, Los Angeles, CA 90095

Justin T. Okano

Lesego busang.

2 The African Comprehensive HIV/AIDS Partnerships (ACHAP), Gaborone, Botswana

Khumo Seipone

Eugenio valdano.

3 Sorbonne Université, INSERM, Institut Pierre Louis d’Epidémiologie et de Santé Publique, IPLESP, F75012, Paris, France

Sally Blower

Author contributions:

Methodology: JS, JTO, JP, EV

Software: JS, JP

Validation: JS, JTO

Formal analysis: JS, JTO, JP

Visualization: JTO, JP

Supervision: SB

Writing—original draft: SB

Writing—review & editing: JS, JTO, JP, LB, KS, EV, SB

Funding acquisition: SB

Associated Data

All data needed to evaluate the conclusions in the paper are presented in the paper, or freely available for registered users at the IPUMS International website: https://doi.org/10.18128/D020.V7.3 . We note that we did not collect these data, nor are they permitted to be posted to other repositories. All code needed to reproduce all parts of this analysis are available from the first author’s GitHub page: https://github.com/janetsong80/pop-mobility-botswana-hiv (copy archived at: https://zenodo.org/badge/latestdoi/594401539 ).

The majority of people with HIV live in sub-Saharan Africa, where HIV epidemics are generalized. For these epidemics to develop, populations need to be mobile. However, population-level mobility has not yet been studied in the context of the development of generalized HIV epidemics. Here we do so by studying historical migration data from Botswana which has one of the most severe generalized HIV epidemics worldwide; in 2021, HIV prevalence was 21%. The country reported its first AIDS case in 1985 when it began to rapidly urbanize. We hypothesize that, during the development of Botswana’s epidemic, the population was highly mobile and there were substantial urban-to-rural and rural-to-urban migratory flows. We test this hypothesis by conducting a network analysis using a historical time series (1981 to 2011) of micro-census data from Botswana. We found 10% of the population moved their residency annually, complex migration networks connected urban with rural areas, and there were very high rates of rural-to-urban migration. Notably, we also found mining towns were both important in-flow and out-flow migration hubs; consequently, there was a very high turnover of residents in towns. Our results support our hypothesis, and together, provide one explanation for the development of Botswana’s generalized epidemic.

INTRODUCTION:

Over 25 million people live with HIV infection in sub-Saharan Africa. All of the HIV epidemics in this continent are generalized: in these type of epidemics, the epidemic is dispersed throughout the country. Therefore, a population needs to be highly mobile in order for a generalized epidemic to have developed. However, although the epidemiology of HIV in sub-Saharan Africa has been widely studied ( Farley et al., 2022 ; Read et al., 2022 ), the role of population-level mobility patterns in the development of generalized HIV epidemics has not been assessed for any country on the continent. Botswana has one of the most severe HIV epidemics worldwide, and reported its first AIDS case in the early 1980s. Since then, the epidemic has become generalized and hyper-endemic: in 2021, HIV prevalence in adults (15 to 64) was 21% ( Mine et al., 2022 ). Simultaneous to the development of its generalized HIV epidemic, Botswana has undergone rapid urbanization. It was predominantly rural in 1981 ( Tarver, 1984 ; Tarver & Miller, 1987 ), but had become predominantly urban ( Statistics Botswana, 2014 ) by 2011; by that time Botswana’s HIV epidemic had stabilized ( CSO Botswana & NACA, 2009 ). We hypothesize that – during the development of Botswana’s epidemic – there was a high level of population mobility, and substantial urban-to-rural migratory flows. We test this hypothesis by conducting a network analysis using a historical time series of micro-census (i.e., individual-level) data collected in Botswana. These data contain information both on migration and urbanization. The time series covers the time period from when the epidemic was first apparent to when it stabilized in 2011. We use these data: (i) to estimate the annual incidence (at the national level) of internal migration over three decades (1981 to 2011), (ii) to characterize migrants on the basis of gender and age, (iii) to reconstruct internal migration networks (in order to identify large-scale population movements and connectivity patterns), (iv) to identify migration hubs, and (v) to quantify urban-to-rural and rural-to-urban migratory flows. The Government of Botswana defines internal migration as residents changing their place of permanent residence within their home country. We discuss our results in the context of understanding the development of the generalized hyper-endemic HIV epidemic in Botswana.

Between 1981 and 2011, the population of Botswana increased from just under one million to over two million ( Statistics Botswana, 2014 ) and the country became progressively urbanized. In 1981, the vast majority of the population lived in small rural villages. Only ~18% of the population were living in urban areas ( Statistics Botswana, 2014 ): either in one of the two cities (Gaborone, the capital, or Francistown) or one of the four towns (Lobatse, Selebi Phikwe, Orapa, and Jwaneng). Lobatse is an administrative center, and the other three are mining towns. Selebi Phikwe is based on copper and nickel mining, and was founded in the early 1970s. Orapa and Jwaneng are based on mining diamonds, and were founded in the late 1960s and early 1980s, respectively. By 2011, only ~36% of the population were living in rural areas ( Statistics Botswana, 2014 ). The remaining population were living in one of the three types of urban centers: a city (there were still only two cities), a town (there were now five towns following the 1991 addition of Sowa, a mining town for soda ash) or an urban village. In Botswana, an urban village is defined as a settlement with at least 5,000 individuals and 75% of the workforce engaged in non-agricultural economic activities ( Statistics Botswana, 2014 ). The urban villages developed by in situ urbanization, which is defined as rural settlements transforming into urban areas by expanding their non-agricultural activities and increasing economic linkages with neighboring areas ( Moriconi-Ebrard et al., 2020 ). Botswana consists of 28 administrative districts ( Figure 2—figure supplement 1 ) ( Okano et al., 2021 ). Each city and town are separate administrative districts, and the remaining 21 districts contain at least one urban village and many rural villages.

A great deal is known about the epidemiology of HIV in Botswana in terms of the increase in prevalence and geographic variation in prevalence ( Magosi et al., 2022 ; Novitsky et al., 2015 ; Novitsky et al., 2020 ; Okano et al., 2021 ). The first case of AIDS was reported in 1985 from the nickel mining town of Selebi Phikwe ( African Natural Resources Center, 2016 ); in the late 1980s, additional cases were reported from two diamond mining towns (Jwaneng and Orapa) and the city, Francistown. In 1990, the first HIV Sentinel Surveillance Survey of antenatal clinic (ANC) attendees (15 to 49 years old) was conducted in Gaborone and in one rural district (Boteti): HIV prevalence was found to be 6% and 4%, respectively ( UNAIDS & World Health Organization, 2004 ). In 1992, a national sentinel surveillance survey began and was conducted annually to 2011. National prevalence in pregnant women was found to have already reached a very high level (18%) by 1992, and continued to increase for the next eight years: 23% by 1993, 32% by 1995, and 39% by 2000 ( UNAIDS & World Health Organization, 2004 ). Based on the ANC data ( UNAIDS & World Health Organization, 2004 ), HIV prevalence in the cities and the mines was the highest. In Francistown, prevalence was 8% in 1991, 24% by 1992, and 44% by 2000. In Gaborone, prevalence rose from 15% in 1992 to 36% by 2000. In the mining town of Selebi Phikwe, prevalence rose from 27% in 1994 to 50% by 2000. By the late 1990s, a high AIDS morbidity and mortality rate in miners led the diamond mining company Debswana to conduct an HIV prevalence survey of its (male) employees. The survey took place in 1999: 29% of miners were found to be infected with HIV ( Barnett et al., 2002 ). The first population-level survey in Botswana to test participants for HIV (BAIS II, the Botswana AIDS Impact Survey) was conducted in 2004 ( NACA & CSO Botswana, 2005 ). It found that, at that time, prevalence in the general population was ~29% in women and ~20% in men aged 15 to 49 years old. Notably, phylogenetic studies have shown that viral lineages have dispersed widely throughout Botswana suggesting large-scale population-level movements have occurred ( Magosi et al., 2022 ; Novitsky et al., 2015 ; Novitsky et al., 2020 ).

Estimated incidence of internal migration:

We estimated the Crude Migration Intensity (CMI) ( Bell et al., 2002 ) for Botswana in the 12 months before each census was conducted (i.e., between 1980 and 1981, 1990 and 1991, 2000 and 2001, and 2010 and 2011). The CMI represents the overall incidence, or level of internal migration, per hundred residents over a specified time interval ( Bell et al., 2002 ); it is an indicator of the propensity of the population to move, and is a measure of migration both between districts and within districts. The CMI remained markedly high and constant: 9.05 per hundred persons between 1980 and 1981, 10.45 per hundred persons between 1990 and 1991, 10.73 per hundred persons between 2000 and 2001, and 10.33 per hundred persons between 2010 and 2011. We found that migrants were more likely to move between districts than to move within districts; this propensity, calculated as a ratio (the number of migrants moving between districts relative to the number of migrants moving within districts), increased over time from 1.25 (1981) to 2.40 (2001 and 2011).

Age profiles of migrants:

Gender-stratified age profiles of migrants are presented in Figure 1 . These results show that the type of individual who migrated within Botswana at the time of each census was markedly similar, with respect to gender and age. Approximately 50% of migrants were women, and the most common age to migrate (for both women and men) was between 16 and 20 years old.

An external file that holds a picture, illustration, etc.
Object name is nihpp-2023.02.01.23285339v1-f0001.jpg

(A) 1981, (B) 1991, (C) 2001, and (D) 2011.

Migratory flows:

To reconstruct the internal migration networks, we first calculated migratory flows between districts. We define a district-level migratory flow as the number of migrants who changed their residency from one district to another during the 12 months prior to the census. The migratory flows for each district (within district, in-flow, and out-flow) are listed in Table 1 (1980-1981), Table 2 (1990-1991), Table 3 (2000-2001), and Table 4 (2010-2011). Over three decades (between 1981 and 2011), the number of migrants more than doubled (as did the population of Botswana).

The table presents each district’s urban/rural classification, population size, total number of migrants (both within and between districts), within district migration intensity per 100 residents (WDMI), and population turnover. The urban/rural classes were: city (C), town (T), or predominantly rural (R). Population sizes are tabulated for residents for whom there was one-year migration data. Turnover was calculated as the net change in the annual migration rate per 100 residents.

Number of Migrants Metric
IDDistrictClassPopulationWithinInOutWDMITurnover
1GaboroneC56,3901207,7506,1100.223.00
2FrancistownC31,090504,7303,8800.172.81
3LobatseT18,11001,9402,2600.00−1.74
4Selebi PhikweT27,210504,0803,8500.190.85
5OrapaT6,33008108900.00−1.25
6JwanengT5,12002,3107100.0045.45
10SouthernR103,1907,1203,3303,3906.90−0.06
11BarolongR14,2901703407701.15−2.92
20South EastR30,8803302,0001,0001.103.35
30KwenengR109,8905,9702,0102,4505.41−0.40
40KgatlengR43,1101,7401,3301,9603.98−1.44
50Central SeroweR93,2605,2403,3504,4805.55−1.20
51Central MahalapyeR78,9103,1202,5303,2003.92−0.84
52Central BobonongR44,1603,1001,8402,7306.88−1.98
53Central BotetiR24,6201,9401,2001,0507.930.61
54Central TutumeR74,8502,7903,2303,9403.69−0.94
55Central Tuli BlockR2,7402053000.9023.98
60North EastR36,3608901,8902,4802.41−1.60
70Ngamiland SouthR58,1602,8701,2801,3604.93−0.14
71Ngamiland NorthR7,760901403101.13−2.14
72ChobeR7,3904803703806.49−0.14
73Okavango DeltaR1,02020150202.2514.61
80GhanziR17,8901,6203905308.99−0.78
90Kgalagadi SouthR14,9106904504204.640.20
91Kgalagadi NorthR8,8801704302401.962.19
Botswana916,52038,59048,41048,4104.210

The table presents each district’s urban/rural classification, population size, total number of migrants (both within and between districts), within district migration intensity per 100 residents (WDMI), and population turnover. The urban/rural classes were: city (C), town (T), predominantly urban (U), partially urban (PU) and predominantly rural (R). Population sizes are tabulated for residents for whom there was one-year migration data. Turnover was calculated as the net change in the annual migration rate per 100 residents.

Number of Migrants Metric
IDDistrictClassPopulationWithinInOutWDMITurnover
1GaboroneC137,68055015,35014,6200.400.53
2FrancistownC67,740708,1106,5100.112.42
3LobatseT26,57003,1302,6400.001.88
4Selebi PhikweT40,480304,9204,1100.082.04
5OrapaT8,12009001,0800.00−2.17
6JwanengT10,52001,4501,8000.00−3.22
7SowaT1,64007804100.0029.13
10SouthernR125,4005,2005,0206,4604.10−1.14
11BarolongR18,0504207001,4002.24−3.73
20South EastU42,3209003,2802,1502.182.74
31Kweneng EastPU140,8106,6508,0606,4804.781.13
32Kweneng WestR27,7001,3601,6301,4304.950.73
40KgatlengPU55,1802,6703,0403,2004.82−0.29
50Central SeroweR126,8107,0507,3506,9005.580.36
51Central MahalapyeR97,0305,3604,6906,3605.43−1.69
52Central BobonongR51,2202,7802,2203,3205.31−2.10
53Central BotetiR37,6302,1901,9801,7705.850.56
54Central TutumeR98,9403,5905,0006,9603.56−1.94
60North EastR41,4701,0402,8303,2902.48−1.10
70Ngamiland SouthPU60,4504,6502,4704,4007.45−3.09
71Ngamiland NorthR33,7506202,6506501.956.30
72ChobePU12,5007001,3201,1005.701.79
80GhanziR26,1101,9708208107.550.04
90Kgalagadi SouthR19,2501,0109107805.280.68
91Kgalagadi NorthR10,9704107407203.740.18
Botswana (Total)1,318,34049,22089,35089,3503.730
Number of Migrants Metric
IDDistrictClassPopulationWithinInOutWDMITurnover
1GaboroneC195,1501,36022,10025,4200.69−1.67
2FrancistownC83,9002209,0609,6200.26−0.66
3LobatseT31,5101703,6604,1200.53−1.44
4Selebi PhikweT50,4201105,4205,6300.22−0.41
5OrapaT8,770101,4501,1000.124.16
6JwanengT15,360702,8702,3300.473.64
7SowaT3,15006807900.00−3.37
10SouthernR113,6004,8506,9407,6704.24−0.64
11BarolongR45,4501,3103,3103,0702.900.53
12Ngwaketse WestU10,8403809301,0003.48−0.64
20South EastU59,6706906,0303,6801.204.10
31Kweneng EastU185,8806,01011,7609,5403.271.21
32Kweneng WestR39,1102,0402,2202,5405.17−0.81
40KgatlengPU73,4602,2704,2504,0903.100.22
50Central SerowePU151,8305,7509,83010,7803.76−0.62
51Central MahalapyePU110,0104,5606,2506,9804.12−0.66
52Central BobonongPU64,3302,7104,0903,4204.261.05
53Central BotetiPU47,1103,7803,1102,1408.192.10
54Central TutumeR121,5204,3906,6308,0303.57−1.14
60North EastR49,9301,3103,7403,6302.630.22
70Ngamiland SouthU69,3702,3504,8904,3203.420.83
71Ngamiland NorthR52,8202,5502,1302,1904.82−0.11
72ChobePU16,2507501,7009904.834.57
80GhanziR32,3703,0901,5702,0309.41−1.40
90Kgalagadi SouthR25,2601,6801,4001,4206.65−0.08
91Kgalagadi NorthR16,5109101,7001,1905.693.19
Botswana (Total)1,673,58053,320127,720127,7203.190

The table presents each district’s urban/rural classification, population size, total number of migrants (both within and between districts), within district migration intensity per 100 residents (WDMI), and population turnover. The urban/rural classes were: city (C), town (T), predominantly urban (U), partially urban (PU) and predominantly rural (R). Population sizes are tabulated for residents for whom there was one-year migration data. Turnover was calculated as the net change in the annual migration rate per 100 residents. The Central Kgalagadi Game Reserve district is denoted CKGR.

Number of Migrants Metric
IDDistrictClassPopulationWithinInOutWDMITurnover
1GaboroneC232,55085022,25029,5400.37−3.04
2FrancistownC99,5801509,42012,8800.15−3.36
3LobatseT28,990603,1503,7900.21−2.16
4Selebi PhikweT49,630604,8305,2600.12−0.86
5OrapaT9,23001,4201,8100.00−4.05
6JwanengT17,910403,2002,9700.221.30
7SowaT3,580108706600.286.23
10SouthernU128,2104,2307,9408,2103.30−0.21
11BarolongR54,3902,1203,8103,3403.900.87
12Ngwaketse WestU13,8106209908904.490.73
20South EastU83,9701,0209,2105,7001.224.36
30KwenengU299,62010,27018,18012,3803.431.97
40KgatlengU90,8402,2805,6905,0602.510.70
50Central SerowePU178,6708,25011,84010,8904.620.53
51Central MahalapyePU118,3504,2506,4507,3703.59−0.77
52Central BobonongPU71,3103,0205,1704,2904.241.25
53Central BotetiPU57,4203,4203,6603,0405.961.09
54Central TutumePU146,5405,4709,6809,4003.730.19
60North EastR59,8002,0504,7504,0603.431.17
70Ngamiland SouthU89,7102,9804,4606,2103.32−1.91
71Ngamiland NorthR59,7803,0601,7702,5205.12−1.24
72ChobeR21,8301,0002,5901,9404.583.07
73Okavango DeltaR2,080503501102.4013.04
80GhanziR42,7203,8102,5201,9508.921.35
81CKGRR23010130.00604.3543.75
90Kgalagadi SouthR29,9001,5401,3901,7205.15−1.09
91Kgalagadi NorthU20,5101,0901,6601,3005.311.79
Botswana (Total)2,011,16061,710147,380147,3803.070

In 1981, the two cities and four towns were fairly small. The capital, Gaborone, had a population of 56,860 and Francistown had a population of only 31,120. The populations of the four towns ranged in size from 5,200 to 27,430. The number of individuals who lived in the other districts ranged in size from 1,050 in the Okavango Delta to 111,820 in Kweneng. All of these districts had an in-flow, and/or out-flow, of migrants ( Table 1 ). More than half also had a fairly high within district migration intensity (WDMI): for example, in Central Bobonong, ~7% of residents moved from one rural village to another between 1980 and 1981. Migratory flows continued to increase over the next three decades ( Tables 2 – 4 ), as the population size increased and the country became progressively urbanized.

We also estimated the annual turnover rate in each district in the year before each census ( Tables 1 – 4 ). This rate is defined as the net change in the district’s annual migration rate per hundred residents. A negative rate signifies that the district’s population size decreased, a positive rate signifies that it increased. By far, the highest turnover rates were in 1981 ( Table 1 ). In that year, the highest turnover rate was in the mining town, Jwaneng.

Reconstructed internal migration networks:

The chord diagrams ( Figure 2 ) show the reconstructed internal migration networks, based on the micro-census data, between: 1980-81 ( Figure 2A ), 1990-91 ( Figure 2B ), 2000-01 ( Figure 2C ), and 2010-11 ( Figure 2D ). Each district is a node in the network. Networks are shown in terms of the magnitude of the migratory in-flows and out-flows between districts, and the effect of these flows on connecting different districts throughout the country: two districts are connected if they have a migratory flow between them. Even the earliest migration network in 1981 – when Botswana was predominantly rural – can be seen to be fairly complex and shows a high degree of connectivity amongst the districts, although the migratory flows and counter-flows are fairly small. Notably, in all four networks, many of the flows and counter-flows are similar in magnitude. Both the magnitude of the migratory flows, and counter-flows, increased with time (between 1981 and 2011) as the population increased in size.

An external file that holds a picture, illustration, etc.
Object name is nihpp-2023.02.01.23285339v1-f0002.jpg

(A) 1981, (B) 1991, (C) 2001, and (D) 2011. Each diagram shows the internal migration network of the general population in the 12 months prior to the census. Each color represents a different administrative district. The thickness of each line is proportional to the number of migrants that moved between the two connected districts. The angular width of each district is proportional to the total number of migrants who moved into, or out of, that district. For clarity, in (A)-(C) only connections with greater than 200 migrants are shown, and in (D) only connections with greater than 400 migrants are shown. Consequently, some districts are not shown in the chord diagram. The total number of migrants (in and out) of every district is listed in Tables 1 – 4 .

In-flow and out-flow migration hubs:

Out-flow and in-flow migration hubs are shown in Figure 3 . Notably, the towns (Jwaneng, Selebi Phikwe, Sowa, Lobatse, and Orapa) were amongst the top five out-flow and in-flow migration hubs at every census period.

An external file that holds a picture, illustration, etc.
Object name is nihpp-2023.02.01.23285339v1-f0003.jpg

(A) Top out-flow hubs in 1981, (B) Top in-flow hubs in 1981, (C) Top out-flow hubs in 1991, (D) Top in-flow hubs in 1991, (E) Top out-flow hubs in 2001, (F) Top in-flow hubs in 2001, (G) Top out-flow hubs in 2011, (H) Top in-flow hubs in 2011. Out-flow hubs are sorted by the rate per hundred residents of a district’s population that moves out to another district; in-flow hubs are sorted by the rate per hundred residents of a district’s population that moves in from another district.

In 1981, the top in-flow hub was the diamond mining town, Jwaneng: 45% of the town’s population consisted of migrants who had moved to the town between 1980 and 1981 ( Figure 3B ). At that time, Jwaneng was also an important out-flow hub: 14% of the town’s population moved to another district between 1980 and 1981 ( Figure 3A ). The ego networks for Jwaneng show the districts that migrants returned to ( Figure 4A ), and the districts that they came from ( Figure 4B ). It can be seen that migrants moved between the other mining towns, the two cities, and many of the rural districts, with the majority of migrants moving between Jwaneng and Southern, a district abutting the town ( Figure 2—figure supplement 1 ).

An external file that holds a picture, illustration, etc.
Object name is nihpp-2023.02.01.23285339v1-f0004.jpg

The chord diagrams show migrants flowing (A) out of, and (B) into, the diamond mining town in 1981.

In 1991 and 2001, the top out-flow (and in-flow) hub was the mining town, Sowa. In 1991, 48% of the population consisted of migrants who had moved to the town in the previous 12 months ( Figure 3D ); 25% of the town’s population moved to another district between 1990 and 1991 ( Figure 3C ). In 2001, 22% of the population consisted of migrants who had moved to the town in the previous 12 months ( Figure 3F ); 25% of the town’s population moved to another district in 2000-2001 ( Figure 3E ).

Francistown was an important in-flow and out-flow hub between 1980 and 1981. In 1981, 15% of the population consisted of migrants who had moved to the town in the previous 12 months ( Figure 3B ); 13% of the town’s population moved to another district between 1980 and 1981 ( Figure 3A ). The city was also an important in-flow hub in 1991 ( Figure 3D ) and out-flow hub in 2011 ( Figure 3G ). Gaborone was only an important out-flow hub in 1991 and 2001 ( Figure 3C and Figure 3E ); it was also an important in-flow hub in 1981 and 2001 ( Figure 3B and Figure 3F ).

Quantifying urban-to-rural and rural-to-urban migratory flows:

The internal migration networks for each of the four time periods (1981, 1991, 2001, and 2011) are presented in terms of our urban/rural classification framework, both as schemata ( Figure 5 ) and as Sankey diagrams ( Figure 6 ). Both the schemata and the Sankey diagrams show the magnitude of the migratory flows amongst the five classes in the 12 months before each census. Between 1981 and 2011 the overall percentage of the population living in urban areas increased more than 3.5 fold (from ~18% to ~64%) ( Statistics Botswana, 2014 ), and the internal migration networks changed substantially over that period.

An external file that holds a picture, illustration, etc.
Object name is nihpp-2023.02.01.23285339v1-f0005.jpg

(A) 1981, (B) 1991, (C) 2001, and (D) 2011. Circles represent the five classes based on an urban/rural classification (see Methods ). The radius of each circle is proportional to the number of residents living in the districts in that specific class. The color of each arrow indicates the size of the net migration between classes. The class designation (city, town, predominantly urban, partially urban, predominantly rural) of each district over time is provided in Tables 1 – 4 .

An external file that holds a picture, illustration, etc.
Object name is nihpp-2023.02.01.23285339v1-f0006.jpg

(A) 1981, (B) 1991, (C) 2001, and (D) 2011. The diagrams show the relative magnitude of migratory flows between the different classes over time. Cities are shown in purple, towns in blue, predominantly urban districts in gold, partially urban districts in magenta, and predominantly rural districts in green. The class designation (city, town, predominantly urban, partially urban, predominantly rural) of each district over time is provided in Tables 1 – 4 .

In 1981, there were only six districts (the two cities and four towns) that were urban areas ( Table 1 , Figure 5A , and Figure 6A ). The other districts contained several hundred small rural villages and no urbanized areas. Almost all districts had an in-flow and out-flow of migrants: notably, rural-to-urban, and counter flows (urban-to-rural) are apparent. Taken together the results show that in 1981 the population was very mobile, the majority of migratory flows were within and amongst rural districts, and rural-to-urban migrations were greater than urban-to-rural counter-flows.

By 1991, due to in situ urbanization, one previously rural district (South East) had become predominantly urban, and four districts (Kweneng East, Kgatleng, Ngamiland South, and Chobe) had become partially urban ( Table 2 ). Migratory flows appear to be fairly symmetrical (in terms of the number of migrants who moved between pairs of districts) throughout the country ( Figure 5B and Figure 6B ). The towns had increased ~1.5 fold since 1981, and the cities had more than doubled in size. As at the time of the previous census, there were considerable rural-to-urban and urban-to-rural migratory flows.

By 2001, considerably more in situ urbanization had occurred: since 1981, four districts had become predominantly urban, six districts partially urban, and only nine districts remained predominantly rural ( Table 3 ). Notably, some of the migratory flows were now asymmetric: the number of migrants moving from the cities to the predominantly urban and partially urban districts were greater than the counter-flows ( Figure 5C and Figure 6C ). The growth rate of towns remained approximately the same as in the previous decade, but the growth rate of the cities had decreased to ~1.4 fold; these rates depended upon the number of migrants, births, and deaths.

By 2011, twelve districts had become predominantly or partially urbanized ( Table 4 ). The growth rate of cities had decreased to ~1.2 fold, and the towns had not increased substantially in size. At this time, while most flows were symmetric, migrants were disproportionately leaving the cities in favor of towns and urban villages ( Figure 5D and Figure 6D ). Approximately 42% of the population were living in these urban villages; there was at least one urban village in every district.

DISCUSSION:

Our results support our hypothesis that – during the development of Botswana’s generalized HIV epidemic (i.e., between 1981 and 2011) – there was a high level of population mobility, and substantial urban-to-rural migratory flows. Using our time series of historical data, we found ~10% of the population moved their residency in the 12 months before each census in 1981, 1991, 2001, and 2011. The constancy of this value at each census suggests that the annual migration rate in Botswana was stable at ~10% for each year between 1981 and 2011. Notably, the type of migrants also remained constant over this time period: evenly split by gender, with younger people more likely to migrate. We found that migration occurred within districts, as well as between districts, but that migration between districts was more common. Even in 1981, at the beginning of the period of rapid urbanization, the migration network was highly connected and therefore large-scale population movements linked almost all of the districts in the country. Notably, we found at each census, that although the number of urbanized districts increased (as did their population size) the vast majority of migratory flows between the rural and developing urban districts were balanced by counter-flows.

Our results provide insights into how HIV could have diffused through the population of Botswana as the epidemic became generalized. The first cases of AIDS in Botswana were reported in the 1980s ( African Natural Resources Center, 2016 ); this suggests that HIV may have been circulating in Botswana in the 1970s, as HIV has an incubation period of 10 to 12 years ( Hendriks et al., 1993 ; Muñoz et al., 1989 ). One phylogenetic study indicates that HIV was introduced into Botswana in the 1960s, but exponential growth did not begin until the 1980s ( Wilkinson et al., 2015 ). Another study suggests that HIV spread southwards from Kinshasa (the epicenter of the HIV epidemic) in the Democratic Republic of the Congo via travelers on the railroad that was built in the 1800s ( Faria et al., 2014 ); notably, this railroad goes through Francistown. It is possible that HIV was introduced into the mining towns of Botswana from overseas or South Africa. In the 1970s and 1980s many young men from Botswana went to work in the mines in South Africa, returning as the mining industry became established in their home country. Our results show that, by 1981, the migration networks (for women and men) linked almost every district in Botswana; this indicates that residents from Francistown ( UNAIDS & World Health Organization, 2004 ) and the mining towns ( Barnett et al., 2002 ) could have seeded sub-epidemics in multiple rural districts throughout the country. HIV could then have spread widely within these districts due to the high rates of within-district migration. Due to the high migration rates, the sub-epidemics in rural districts were likely to have been seeded numerous times. When HIV was first introduced into rural districts, the sub-epidemics in these districts may have been maintained by source-sink dynamics ( Okano et al., 2020 ): districts where transmission was high enough to be self-sustaining (e.g., mining towns) could have maintained sub-epidemics in rural districts where transmission was too low to be self-sustaining. Our results suggest that migratory flows between high and low prevalence areas are likely to have existed, at least from the early 1980s onwards. These flows would have functioned as “transmission corridors” ( Okano et al., 2021 ), and created high risk flows ( Valdano et al., 2021 ). Taken together, our quantitative results show that the highly-connected migration networks, with high migratory flows and increasing urbanization, could have contributed to the development of Botswana’s generalized epidemic.

Our results also provide insights into how the Botswana epidemic could have become hyper-endemic. Migrants have been shown to be at high risk of acquiring HV infection ( Camlin & Charlebois, 2019 ; Dobra et al., 2017 ; Olawore et al., 2018 ); our results suggest that, due to migration alone, each year ~10% of the population of Botswana was at high risk. The migrants in Botswana were young women and men, a group generally at the highest risk of infection through sexual transmission. Mining towns and cities would have been transmission hot-spots (for both women and men) due to high levels of risky sexual behavior (as evidenced by high prevalence ( Barnett et al., 2002 )) and a high probability of encountering an HIV-infected individual as a sex partner. Due to the high prevalence ( Barnett et al., 2002 ; UNAIDS & World Health Organization, 2004 ), high turnover rates, and a high level of connectivity to other districts, the mining towns and cities could have functioned as geographically-defined core groups for sexual transmission (for both women and men). Some individuals could have been part of a high-risk core group when they were living in the cities and mining towns, and decreased their risky behavior when they returned to the rural districts. For example, FSWs in Botswana can have an extremely high number of sex partners in mining towns and urban areas (~7 per week) and remain as FSWs for ~4 years ( Ministry of Health Botswana, 2013 ). Notably we found that the mining towns and cities were both migration in-flow hubs and migration out-flow hubs. The fact that they were migration in-flow hubs would have led to a high in-flow of uninfected individuals from rural districts; the fact that they were migration out-flow hubs would have led to a high out-flow of HIV-infected individuals to rural districts. Therefore, rural-to-urban migration coupled with urban-to-rural migration could have been a very important driver in the Botswana epidemic becoming hyper-endemic. Without this driver, prevalence could have reached very high levels in the mining towns and cities, but remained fairly low in rural districts.

Our study has several limitations. First, we have examined only internal migration (i.e., migration that takes place within a country), and have not included international migration. However, for the three-decade time period that we have investigated, the number of international migrants to Botswana was much smaller than the number of internal migrants ( Statistics Botswana, 2014 ; Statistics Botswana, 2015 ). Second, all of our analyses are focused on Botswana in order to analyze the impact of migration networks on a hyper-endemic epidemic; we have not evaluated migration networks in other African countries. We recommend that studies, such as we have conducted here, are conducted in these other countries by analyzing a time series of historical micro-census data, especially in the only other two countries in sub-Saharan Africa that have generalized hyper-endemic HIV epidemics: Eswatini and Lesotho. Third, a potential limitation is that the IPUMS data are samples rather than the complete censuses; we used this approach as it is generally not possible to obtain complete census data due to privacy reasons. Fourth, while we have gained insights into understanding the development of the generalized hyper-endemic HIV epidemic in Botswana, we have not conducted an analysis of the transmission dynamics of the epidemic in order to explore these insights. This is the focus of current research in which we are analyzing a geospatial mathematical model of the epidemiological evolution of the HIV epidemic in Botswana. Fifth, we have examined only two processes that affected the evolution of Botswana’s epidemic, other factors could also have been important.

All HIV epidemics in SSA are generalized; HIV can only diffuse through a population by the movement of people. To our knowledge, our study is the first to study population-level mobility patterns and rates in any country in sub-Saharan Africa. Taken together, our results identify particular characteristics of migration in Botswana that could explain why HIV prevalence rose to very high levels in many districts throughout the country. Migratory behavior in other countries with generalized epidemics may be very different. For example, migrants in other African countries tend to have a very male-biased sex ratio; consequently, women are unlikely to have seeded multiple sub-epidemics throughout these countries. If mining towns and cities in these countries were not migratory in-flow and out-flow hubs with high turnover of their populations, as in Botswana, they would not have served as major sources of infection. Botswana has recently achieved the 95-95-95 targets that UNAIDS specified needed to be reached by 2030 in order to eliminate HIV: 95% of people living with HIV to be diagnosed, 95% of the diagnosed to be on treatment, and 95% of those on treatment to be virally suppressed ( Mine et al., 2022 ). However, transmission is continuing ( Magosi et al., 2022 ), as is urbanization and large-scale population movements ( Okano et al., 2021 ). The same processes that we have identified that may have contributed to the development of the generalized hyper-endemic HIV epidemic in Botswana may prevent the elimination of HIV.

MATERIALS AND METHODS:

To conduct our analyses, we used representative samples of micro-census data extracted from the IPUMS-International (Integrated Public Use Microdata Series-International) database ( Minnesota Population Center, 2020 ): micro-census data are anonymized individual-level data ( Ruggles et al., 2015 ). IPUMS-International currently disseminates data from 547 censuses and surveys in 103 countries worldwide ( Minnesota Population Center, 2020 ). The Botswana dataset consists of representative 10% samples from the original censuses, and includes anonymized individual-level data on age, gender, and residence (current, as well as 12 months prior). The data also includes individual survey weights that allow for population-level estimation ( Ruggles et al., 2015 ).

To examine historical trends in internal migration networks and urbanization, we analyzed data collected in the 1981, 1991, 2001, and 2011 censuses. We used data on internal migration that had occurred in the 12 months prior to each census: an individual was classified as a migrant if they had changed their permanent residency within that 1-year interval.

Estimating the incidence of internal migration:

We estimated the incidence of internal migration by calculating the CMI. This statistic represents the overall incidence, or level of internal migration (between district plus within district), per hundred residents over a year. The mathematical definition of the CMI is given in Equation 1 :

where M is the total number of internal migrants and P is the population size of Botswana in a given year. M is calculated as M = Σ i D i = Σ i O i where D i represents the in-flows to each district i , and O i represents the out-flows to each district i .

Constructing gender-stratified age structure pyramids:

We aggregated the migration data from each census by gender and age (using 5-year age groupings) to construct population pyramids; these pyramids show the age-gender demographics of all individuals internally migrating in the 12 months prior to each census.

Calculating migratory flows:

We define a district-level migratory flow as the number of migrants who change their residency from one district to another during the 12 months prior to the census. We calculated migratory flows between each pair of districts. The country consists of 28 administrative districts ( Okano et al., 2021 ). Each city and town are separate administrative districts.

Calculating annual turn-over rates:

We defined the annual turn-over rate for a district as the net change in its annual migration rate per hundred residents.

Here the turn-over rate TO i for district i is a function of the number of in-migrants D i and out-migrants O i within the past year, and its population P i at the beginning of the year.

Reconstructing migration networks:

We used the micro-census data to construct Origin-Destination (OD) matrices: the origin was the district that an individual lived in 12 months prior to the census, the destination was the district they lived in at the time of the census. Coefficients of these matrices specify the number of migrants who moved between each pair of districts in the 12 months prior to each census.

Identifying in-flow and out-flow migration hubs:

Migration hubs are those districts where recent migration has a sizeable impact on the size of the resident population, by either bringing it down (out-flow hubs) or increasing it (in-flow hubs) above the average. In-flow migration hubs were identified by calculating the total number of in-migrants to a node/district and dividing by the district’s population that year. Out-flow migration hubs were identified by calculating the total number of out-migrants to a node/district and dividing by the district’s population that year. We list the top five in-flow and out-flow hubs each census year.

To measure and visualize migratory flows during urbanization, we developed a classification framework. This framework enables us to measure the magnitude of migratory flows and counter-flows between districts (based upon their degree of urbanization), and distinguishes between urban-to-rural and rural-to-urban flows.

The framework consists of five classes. The classes (at any point in time) are defined based on the degree of urbanization of the districts at the time of the most recent census: (i) predominantly rural (< 40% of the population living in urban areas), (ii) partially urban (40%-60% of the population living in urban areas), (iii) predominantly urban (> 60% of the population living in urban areas), (iv) town (100% of the population living in urban areas), or (v) city (100% of the population living in urban areas). Classes (i), (ii), and (iii) only include urban areas that develop by in situ urbanization, i.e., rural villages transforming into urban villages. There are migratory flows and counter-flows between the five classes, and flows within each of the five classes: therefore, 25 migratory flows/counter-flows are possible.

The classification framework can be visualized in two formats: (i) a schematic figure showing the five classes, the number of residents in each class, and migratory flows (and counter-flows) amongst and within the classes, and (ii) Sankey diagrams showing the number of migrants who move within, and between, the five classes.

We used the micro-census data to parameterize the framework for Botswana. Specifically, we estimated – at the time of each census – the population size of each district, and determined how many individuals lived in rural and urban areas. We used these estimates to classify all districts into one of the five classes. We then calculated the migratory flows and counter-flows between, and within, districts in the 12 months prior to each survey. These parameter estimates are given for 1981 census ( Table 1 ), the 1991 census ( Table 2 ), the 2001 census ( Table 3 ), and the 2011 census ( Table 4 ). Between each census, districts change in size due to births, deaths and migration. Also, between each census, districts can be reclassified based upon their increase in urbanization.

Supplementary Material

Acknowledgements:.

We acknowledge the Central Statistics Office (CSO) Botswana for collecting the data archived by IPUMS. We are grateful to Nelson Freimer for discussions throughout the course of this research.

JS, JTO, JP, and SB acknowledge the financial support of the National Institute of Allergy and Infectious Diseases, National Institutes of Health grants R56 AI152759 and R01 AI167713.

Competing interests: All authors declare that they have no competing interests.

Data availability statement:

References:.

  • African Natural Resources Center. (2016). Debswana Diamond Company and the Government of Botswana: an HIV/AIDS public-private partnership programme . https://www.afdb.org/fileadmin/uploads/afdb/Documents/Publications/anrc/DEBSWANA_ENG_HR_PAGES.pdf [Accessed October 17, 2022].
  • Barnett T., Fantan T., Mbakile B., & Whiteside A. (2002). The private sector responds to the epidemic: Debswana–a global benchmark . Geneva, Switzerland: UNAIDS. [ Google Scholar ]
  • Bell M., Blake M., Boyle P., Duke-Williams O., Rees P., Stillwell J., & Hugo G. (2002). Cross-national comparison of internal migration: Issues and measures . Journal of the Royal Statistical Society: Series A (Statistics in Society) , 165 ( 3 ), 435–464. [ Google Scholar ]
  • Camlin C. S., & Charlebois E. D. (2019). Mobility and its effects on HIV acquisition and treatment engagement: Recent theoretical and empirical advances . Current HIV/AIDS Reports , 16 ( 4 ), 314–323. doi: 10.1007/s11904-019-00457-2 [ PMC free article ] [ PubMed ] [ CrossRef ] [ Google Scholar ]
  • CSO Botswana, & NACA. (2009). 2008 Botswana AIDS Impact Survey III: Statistical Report . https://microdata.statsbots.org.bw/index.php/catalog/13/download/77 [Accessed October 12, 2022].
  • Dobra A., Barnighausen T., Vandormael A., & Tanser F. (2017). Space-time migration patterns and risk of HIV acquisition in rural South Africa . AIDS , 31 ( 1 ), 137–145. doi: 10.1097/qad.0000000000001292 [ PMC free article ] [ PubMed ] [ CrossRef ] [ Google Scholar ]
  • Faria N. R., Rambaut A., Suchard M. A., Baele G., Bedford T., Ward M. J., Tatem A. J., Sousa J. D., Arinaminpathy N., Pépin J., Posada D., Peeters M., Pybus O. G., & Lemey P. (2014). The early spread and epidemic ignition of HIV-1 in human populations . Science , 346 ( 6205 ), 56–61. doi: 10.1126/science.1256739 [ PMC free article ] [ PubMed ] [ CrossRef ] [ Google Scholar ]
  • Farley S. M., Wang C., Bray R. M., Low A. J., Delgado S., Hoos D., Kakishozi A. N., Harris T. G., Nyirenda R., Wadonda N., Li M., Amuri M., Juma J., Kancheya N., Pietersen I., Mutenda N., Natanael S., Aoko A., Ngugi E. W., Asiimwe F., Lecher S., Ward J., Chikwanda P., Mugurungi O., Moyo B., Nkurunziza P., Aibo D., Kabala A., Biraro S., Ndagije F., Musuka G., Ndongmo C., Shang J., Dokubo E. K., Dimite L. E., McCullough-Sanden R., Bissek A. C., Getaneh Y., Eshetu F., Nkumbula T., Tenthani L., Kayigamba F. R., Kirungi W., Musinguzi J., Balachandra S., Kayirangwa E., Ayite A., West C. A., Bodika S., Sleeman K., Patel H. K., Brown K., Voetsch A. C., El-Sadr W. M., & Justman J. J. (2022). Progress towards the UNAIDS 90-90-90 targets among persons aged 50 and older living with HIV in 13 African countries . J Int AIDS Soc , 25 Suppl 4 , e26005. doi: 10.1002/jia2.26005 [ PMC free article ] [ PubMed ] [ CrossRef ] [ Google Scholar ]
  • Hendriks J. C., Medley G. F., van Griensven G. J., Coutinho R. A., Heisterkamp S. H., & van Druten H. A. (1993). The treatment-free incubation period of AIDS in a cohort of homosexual men . AIDS , 7 ( 2 ), 231–239. doi: 10.1097/00002030-199302000-00012 [ PubMed ] [ CrossRef ] [ Google Scholar ]
  • Magosi L. E., Zhang Y., Golubchik T., DeGruttola V., Tchetgen Tchetgen E., Novitsky V., Moore J., Bachanas P., Segolodi T., Lebelonyane R., Pretorius Holme M., Moyo S., Makhema J., Lockman S., Fraser C., Essex M. M., & Lipsitch M. (2022). Deep-sequence phylogenetics to quantify patterns of HIV transmission in the context of a universal testing and treatment trial - BCPP/Ya Tsie trial . eLife , 11 , e72657. doi: 10.7554/eLife.72657 [ PMC free article ] [ PubMed ] [ CrossRef ] [ Google Scholar ]
  • Mine M., Stafford K., Laws R., Marima R., Lekone P., Ramaabya D., Makhaola K., Mapondera P., Wray-Gordon F., & Agbakwuru C. (2022). Botswana achieved the Joint United Nations Programme on HIV/AIDS (UNAIDS) 95-95-95 targets: results from the Fifth Botswana HIV/AIDS Impact Survey (BAIS V), 2021 . Paper presented at the The 24th International AIDS Conference, Montreal, Canada. [ Google Scholar ]
  • Ministry of Health Botswana. (2013). 2012 Mapping, Size Estimation & Behavioral and Biological Surveillance Survey (BBSS) of HIV/STI Among Select High-Risk Sub-Populations in Botswana . https://www.fhi360.org/sites/default/files/media/documents/BBSS_2012_Final.pdf [Accessed October 3, 2022].
  • Minnesota Population Center. (2020). Integrated Public Use Microdata Series (IPUMS), International: Version 7.3 [dataset]. 10.18128/D020.V7.3 [Accessed June 4, 2020]. [ CrossRef ]
  • Moriconi-Ebrard F., Heinrigs P., & Trémolières M. (2020). Africa’s Urbanisation Dynamics 2020: Africapolis, Mapping a New Urban Geography: OECD, Sahel and West Africa Club Secretariat .
  • Muñoz A., Wang M. C., Bass S., Taylor J. M., Kingsley L. A., Chmiel J. S., & Polk B. F. (1989). Acquired immunodeficiency syndrome (AIDS)-free time after human immunodeficiency virus type 1 (HIV-1) seroconversion in homosexual men . Multicenter AIDS Cohort Study Group. American Journal of Epidemiology , 130 ( 3 ), 530–539. doi: 10.1093/oxfordjournals.aje.a115367 [ PubMed ] [ CrossRef ] [ Google Scholar ]
  • NACA, & CSO Botswana. (2005). Botswana AIDS Impact Survey II: Statistical Report . https://www.statsbots.org.bw/sites/default/files/stats_report%282%29.pdf [Accessed October 12, 2022].
  • Novitsky V., Kuhnert D., Moyo S., Widenfelt E., Okui L., & Essex M. (2015). Phylodynamic analysis of HIV sub-epidemics in Mochudi, Botswana . Epidemics , 13 , 44–55. doi: 10.1016/j.epidem.2015.07.002 [ PMC free article ] [ PubMed ] [ CrossRef ] [ Google Scholar ]
  • Novitsky V., Zahralban-Steele M., Moyo S., Nkhisang T., Maruapula D., McLane M. F., Leidner J., Bennett K., Wirth K. E., Gaolathe T., Kadima E., Chakalisa U., Pretorius Holme M., Lockman S., Mmalane M., Makhema J., Gaseitsiwe S., DeGruttola V., & Essex M. (2020). Mapping of HIV-1C transmission networks reveals extensive spread of viral lineages across villages in Botswana treatment-as-prevention trial . The Journal of Infectious Diseases , 222 ( 10 ), 1670–1680. doi: 10.1093/infdis/jiaa276 [ PMC free article ] [ PubMed ] [ CrossRef ] [ Google Scholar ]
  • Okano J. T., Busang L., Seipone K., Valdano E., & Blower S. (2021). The potential impact of country-level migration networks on HIV epidemics in sub-Saharan Africa: the case of Botswana . The Lancet. HIV , 8 ( 12 ), e787–e792. doi: 10.1016/s2352-3018(21)00267-8 [ PMC free article ] [ PubMed ] [ CrossRef ] [ Google Scholar ]
  • Okano J. T., Sharp K., Valdano E., Palk L., & Blower S. (2020). HIV transmission and source-sink dynamics in sub-Saharan Africa . The Lancet. HIV , 7 ( 3 ), e209–e214. doi: 10.1016/s2352-3018(19)30407-2 [ PMC free article ] [ PubMed ] [ CrossRef ] [ Google Scholar ]
  • Olawore O., Tobian A. A. R., Kagaayi J., Bazaale J. M., Nantume B., Kigozi G., Nankinga J., Nalugoda F., Nakigozi G., Kigozi G., Gray R. H., Wawer M. J., Ssekubugu R., Santelli J. S., Reynolds S. J., Chang L. W., Serwadda D., & Grabowski M. K. (2018). Migration and risk of HIV acquisition in Rakai, Uganda: a population-based cohort study . The Lancet. HIV , 5 ( 4 ), e181–e189. doi: 10.1016/s2352-3018(18)30009-2 [ PMC free article ] [ PubMed ] [ CrossRef ] [ Google Scholar ]
  • Read S. W., Kim P., Marovich M., Dieffenbach C. W., & Fauci A. S. (2022). Forty years of investment in HIV research: progress towards ending the HIV pandemic and preparation for future pandemics . J Int AIDS Soc , 25 ( 12 ), e26039. doi: 10.1002/jia2.26039 [ PMC free article ] [ PubMed ] [ CrossRef ] [ Google Scholar ]
  • Ruggles S., McCaa R., Sobek M., & Cleveland L. (2015). The IPUMS collaboration: integrating and disseminating the world’s population microdata . Journal of Demographic Economics , 81 ( 2 ), 203–216. [ PMC free article ] [ PubMed ] [ Google Scholar ]
  • Statistics Botswana. (2014). Population and Housing Census 2011: Analytical Report . https://www.statsbots.org.bw/sites/default/files/publications/Population%26HousingCensusDisseminationanalyticalreport.pdf [Accessed March 21, 2022].
  • Statistics Botswana. (2015). Population and Housing Census 2011: National Statistical Tables . https://www.statsbots.org.bw/sites/default/files/publications/national_statisticsreport.pdf [Accessed April 20, 2020].
  • Tarver J. D. (1984). Migration of the population of Botswana . Pula , 4 ( 1 ), 80–94. [ PubMed ] [ Google Scholar ]
  • Tarver J. D., & Miller H. M. (1987). Rural-urban migration in a developing country: Botswana, Africa . Africa Quarterly , 24 ( 1-2 ), 22–33. [ PubMed ] [ Google Scholar ]
  • UNAIDS, & World Health Organization. (2004). Botswana: Epidemiological Fact Sheet - 2004 Update . https://data.unaids.org/publications/fact-sheets01/botswana_en.pdf [Accessed October 11, 2022].
  • Valdano E., Okano J. T., Colizza V., Mitonga H. K., & Blower S. (2021). Using mobile phone data to reveal risk flow networks underlying the HIV epidemic in Namibia . Nature Communications , 12 ( 1 ), 2837. doi: 10.1038/s41467-021-23051-w [ PMC free article ] [ PubMed ] [ CrossRef ] [ Google Scholar ]
  • Wilkinson E., Engelbrecht S., & de Oliveira T. (2015). History and origin of the HIV-1 subtype C epidemic in South Africa and the greater southern African region . Scientific Reports , 5 , 16897. doi: 10.1038/srep16897 [ PMC free article ] [ PubMed ] [ CrossRef ] [ Google Scholar ]

U.S. flag

A .gov website belongs to an official government organization in the United States.

A lock ( ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

  • Where We Work
  • Global HIV Epidemic
  • Global TB Epidemic
  • Transforming Global HIV & TB Programs with Data
  • Our Approach
  • Stories of Hope
  • Raising Awareness

Related Topics:

  • Show All Home
  • Global Health

HIV and TB Overview: Botswana

At a glance.

CDC supports HIV testing, treatment, and prevention services in Botswana through the U.S. President's Emergency Plan for AIDS Relief (PEPFAR). CDC also collaborates with partners to support tuberculosis (TB) prevention and control programs in Botswana.

Flag of the country of Botswana

CDC partnered with Botswana in 1995, translating health research into methods to strengthen TB prevention and control. In 2000, this partnership expanded to include HIV prevention and treatment. It also included development of programs to maximize the quality, coverage, and impact of Botswana's HIV response.

With support from PEPFAR, CDC continues to collaborate with Botswana’s Ministry of Health and Wellness (MOH) in their response to HIV. This partnership includes services like HIV testing, antiretroviral treatment (ART), prevention of mother-to-child transmission of HIV, and voluntary medical male circumcision. CDC also works with Botswana on TB prevention and control programs.

Download CDC Botswana's Fact Sheet‎

Hiv and tb key data.

Estimated HIV Prevalence (Ages 15-49)

Estimated AIDS Deaths (Age≥15)

Estimated Orphans Due to AIDS

Reported Number Receiving Antiretroviral Therapy (Age≥15)

Tuberculosis (TB)

Estimated TB Incidence

TB Patients with Known HIV-Status who are HIV-Positive

TB Treatment Success Rate

Surpassing UNAIDS treatment targets

Botswana demonstrates that HIV epidemic control is achievable, as measured by the 2030 UNAIDS 95-95-95 targets .

Progress toward targets for ending HIV‎

In 2021, results from the Fifth Botswana AIDS Impact Survey (BAIS V) showed that Botswana has surpassed the UNAIDS targets at 95-98-98.

The country also became the first high-burden country to attain WHO "Silver Tier" certification .

Key activities and accomplishments

Addressing gaps in hiv services.

CDC and partners will continue to use BAIS V data to prioritize service provision, especially among adolescents and young adults. This includes improving services for adolescents and children living with HIV and key populations. CDC and partners are also working to reduce HIV-related mortality, improve TB preventive therapy coverage, and cervical cancer screening. Additional efforts include implementing a 6-month multi-month dispensing of ART. CDC and partners also emphasize equity for all persons seeking HIV services regardless of age, sex, sexual orientation, or geography.

CDC is also supporting Botswana to achieve WHO “Gold Tier” certification. This includes maintaining mother-to-child transmission of less than 5% and providing prenatal care and ART to over 95% of pregnant women. CDC and partners are working to eliminate mother-to-child transmission of syphilis and hepatitis through clinical mentorship and continuous quality improvement (CQI).

Building public health capacity

CDC collaborates with partners, including the Botswana MOH, to enhance the country's public health capacity. These efforts include:

  • Extending the Clinical Mentorship Program to all health districts.
  • Enhancing public health capacity to provide high-quality HIV services.
  • Expanding MOH research capacity through partnerships and training.
  • Maintaining CQI for targeted interventions in all health districts.
  • Strengthening public health infrastructure through FETP .
  • Supporting Botswana's Public Health Institute to sustain HIV control.

Enhancing disease surveillance

CDC-Botswana provides testing for recent HIV infection in the traditional HIV testing services program, linking results to case-based surveillance. Recency testing was expanded from 10 to over 200 sites, offering essential information about new HIV diagnoses, infections, and ongoing transmission.

Public health teams use experience gained from COVID-19 case control efforts to improve HIV case-finding. Rapid response teams use recency and case-based surveillance and geographic data to identify recent HIV infection clusters. Tailored HIV testing strategies and effective treatment and prevention interventions, including PrEP are also provided to populations at highest risk.

Strengthening health information systems

CDC continues to work with Botswana’s MOH to strengthen the country’s data and laboratory systems. Data support includes improving electronic medical records, system interoperability, and analysis capacity. Laboratory support includes quality diagnostics, emerging infection response capacity development, and assistance with accrediting testing and laboratory sites with international standards.

Support for CDC's global HIV and TB efforts‎

Our success is built on the backbone of science and partnerships.

  • Botswana: Development news, research, data | World Bank
  • AIDSinfo | UNAIDS
  • Global Tuberculosis Programme (who.int)

Global HIV and TB

CDC's Division of Global HIV and TB works to end the public health threats of HIV and TB in partnership with countries supported by PEPFAR.

IMAGES

  1. IJERPH

    botswana hiv aids case study geography

  2. HIV AIDS

    botswana hiv aids case study geography

  3. (PDF) The geography of HIV/AIDS prevalence rates in Botswana

    botswana hiv aids case study geography

  4. U1 Population

    botswana hiv aids case study geography

  5. Botswana Presentation Hiv

    botswana hiv aids case study geography

  6. [PDF] The geography of HIV/AIDS prevalence rates in Botswana

    botswana hiv aids case study geography

VIDEO

  1. Week 3 Case Study Presentation

  2. HIV/AIDS in Botswana (ka Setswana)

  3. Botswana Orphan Project

  4. 13. HIV/AIDS Case. The Rhabditid Nematode Roundworm Strongyloides Stercoralis Geohelminths Part 3

COMMENTS

  1. The geography of HIV/AIDS prevalence rates in Botswana

    Current estimates from the Department of HIV/AIDS Prevention and Care of HIV prevalence among pregnant women aged 15-49 years attending antenatal care in public health clinics in Botswana was 31.8%. 3 The national HIV prevalence amongst the women surveyed showed a decline in prevalence from 36.2% in 2001 to 31.8% in 2009.

  2. The geography of HIV/AIDS prevalence rates in Botswana

    However, there is significant geographic variation at the district level, and HIV prevalence is heterogeneous in Botswana's districts, with the highest prevalence recorded in Selebi-Phikwe and North East (41.6%). This is followed by Tutume (41.1%), Bobirwa (39.6%), and Chobe (39.3%). Hukuntsi has the lowest prevalence rate (16.1%).

  3. Spatial Analysis of HIV Infection and Associated Risk Factors in Botswana

    2.1. Study Area. The study was carried out in Botswana, a land-locked country in southern Africa that shares borders with South Africa, Zimbabwe, Zambia and Namibia and covers an area of about 582,000 km 2.Botswana is divided into 10 administrative districts, and further divided into 28 sub-districts, also referred to as census districts (as of 2011).

  4. IJERPH

    Botswana has the third highest human immunodeficiency virus (HIV) prevalence globally, and the severity of the epidemic within the country varies considerably between the districts. This study aimed to identify clusters of HIV and associated factors among adults in Botswana. Data from the Botswana Acquired Immunodeficiency Syndrome (AIDS) Impact Survey IV (BIAS IV), a nationally representative ...

  5. The geography of HIV/AIDS prevalence rates in Botswana

    Botswana has the second-highest HIV infection. rate in the world after Swaziland, with one in three adults. infected. In 2007-2009, the HIV prevalence among males. and females aged 15-49 years ...

  6. The geography of HIV/AIDS infection in Botswana

    Policy Implication Based on large population cross-sectional household survey, this study shows a clear geographic distribution of the HIV/AIDS epidemic in Botswana with highest incidence of infection is in the east-central districts. Generally, it is most prevalent in the northern and north-eastern districts of the country. HIV prevalence is also quite high in the central districts.

  7. PDF Country progress report

    Botswana has an estimated HIV prevalence of 20.68% and has consistently been among the highest in the Eastern and southern African region. The number of new HIV infections steadily decreased since 2010 from 14000 to 9500 representing a decline of 34%. A third of new HIV infections occured among young people 15-24 years, with 70% of the ...

  8. The geography of HIV/AIDS prevalence rates in Botswana

    Results: Overall, HIV prevalence was 17.6%, which was higher among females (20.4%) than males (14.3%). HIV prevalence was higher in cities and towns (20.3%) than in urban villages and rural areas (16.6% and 16.9%, respectively). We also observed an inverse U-shape association between age and prevalence of HIV, which had a different pattern in ...

  9. Botswana's HIV Response: Policies, Context and Future Directions

    Reports suggest the first case of HIV in Botswana was reported in 1985 (Hardon, ... the challenges of antiretroviral treatment: studies from Botswana, Tanzania and Uganda 2006. Geneva, Switzerland: World Health Organization. ... [Google Scholar] Joint United Nations Programme on HIV/AIDS (UNAIDS). (2014). 90-90-90. An ambitious treatment target ...

  10. The geography of HIV/AIDS infection in Botswana

    Abstract. Background Botswana's Human immunodeficiency virus (HIV)/Acquired immune deficiency syndrome (AIDS) epidemic is unprecedented in magnitude and impact. For two decades starting during the ...

  11. The geography of HIV/AIDS prevalence rates in Botswana

    This study showed a clear geographic distribution of the HIV epidemic, with the highest prevalence in the east-central districts, and observed an inverse U-shape association between age and prevalence of HIV, which had a different pattern in males and females. Background Botswana has the second-highest human immunodeficiency virus (HIV) infection rate in the world, with one in three adults ...

  12. The geography of HIV/AIDS prevalence rates in Botswana

    Search life-sciences literature (42,603,876 articles, preprints and more) Search. Advanced search

  13. The geography of HIV/AIDS infection in Botswana

    A clear geographic distribution of the HIV/AIDS epidemic in Botswana with highest incidence of infection is in the east-central districts, and a U-shape association between age and the prevalence of HIV is observed. Background Botswana's Human immunodeficiency virus (HIV)/Acquired immune deficiency syndrome (AIDS) epidemic is unprecedented in magnitude and impact. For two decades starting ...

  14. The geography of HIV/AIDS prevalence rates in Botswana

    HIV/AIDS - Research and Palliative Care Dovepress open access to scientific and medical research O riginal R esearch Open Access Full Text Article The geography of HIV/AIDS prevalence rates in Botswana This article was published in the following Dove Press journal: HIV/AIDS - Research and Palliative Care 17 July 2012 Number of times this article has been viewed Ngianga-Bakwin Kandala 1 Eugene ...

  15. The geography of HIV/AIDS prevalence rates in Botswana

    The health clinics in Botswana was 31.8%.3 The national HIV BAIS III is the most recent AIDS impact survey conducted prevalence amongst the women surveyed showed a decline in in Botswana and is a valuable resource for population-based prevalence from 36.2% in 2001 to 31.8% in 2009.

  16. The geography of HIV/AIDS prevalence rates in Botswana

    Search worldwide, life-sciences literature Search. Advanced Search Coronavirus articles and preprints Search examples: "breast cancer" "breast cancer"

  17. PDF Botswana Case Study

    HIV/AIDS and other communicable diseases causing about half the deaths7. Botswana has an established generalized HIV epidemic with an estimated prevalence of 30.4% in the 15-49-year-old antenatal population8. Purpose of the review This exercise was undertaken to assess the level of integration of select intersectoral

  18. The Geography of HIV/AIDS Prevalence Rates in Botswana

    The Geography of HIV/AIDS Prevalence Rates in Botswana: Author(s) Ngianga-Bakwin Kandala. ... household, and area factors using the 2008 Botswana HIV Impact Survey. Results: Overall, HIV prevalence was 17.6%, which was higher among females (20.4%) than males (14.3%). HIV prevalence was higher in cities and towns (20.3%) than in urban villages ...

  19. Botswana HIV/AIDS Impact Survey Report

    The Fifth Botswana AIDS Impact Survey, BAIS V, was a household-based national survey among adults (defined as individuals aged 15 to 64 years) and children (defined as individuals aged 6 weeks to 14 years) conducted from March 2021 to August 2021 to measure the impact of the national HIV response. The survey offered HIV counseling and testing with return of results to the participants and ...

  20. Botswana Hiv case study

    HIV and AIDS Botswana HIV AIDS Case Study: Facts:-320,000 people living with HIV/AIDS-21% of HIV adult prevalence.-5,800 AIDS death-9,100 new HIV infections Botswana has the third highest HIV prevalence in the world, 21%. Research shows that the population mostly affected by HIV/AIDS is the female sex workers and men who have sex with men.

  21. U1 Population

    Botswana - The impact of HIV/AIDS on the population 1. Place Specific Reference - A problem in southern Africa in the late 1980s - 1985: the first case of HIV/AIDS - 1998-2009: ¼ Botswanans between the age 15-49 was HIV-positive - 2000: 44% of pregnant women had HIV/AIDS - 2001: Peak death - 2016: 1% natural population growth rate, 22% adults ...

  22. HIV In Africa

    This video is a brief introduction to the Geography case study of HIV/AIDS in Botswana, discussing impacts & prevention!LIKE, COMMENT & SUBSCRIBE!

  23. Population mobility and the development of Botswana's generalized HIV

    The first population-level survey in Botswana to test participants for HIV (BAIS II, the Botswana AIDS Impact Survey) was conducted in 2004 (NACA & CSO Botswana, 2005). It found that, at that time, prevalence in the general population was ~29% in women and ~20% in men aged 15 to 49 years old.

  24. HIV and TB Overview: Botswana

    CDC partnered with Botswana in 1995, translating health research into methods to strengthen TB prevention and control. In 2000, this partnership expanded to include HIV prevention and treatment. It also included development of programs to maximize the quality, coverage, and impact of Botswana's HIV response.