Skip to main content
Download PDF

Associations of inter-annual rainfall decreases with subsequent HIV outcomes for persons with HIV on antiretroviral therapy in Southern Africa: a collaborative analysis of cohort studies

BMC Infectious Diseases volume 23, Article number: 889 (2023) Cite this article

Abstract

Background

Periods of droughts can lead to decreased food security, and altered behaviours, potentially affecting outcomes on antiretroviral therapy (ART) among persons with HIV (PWH). We investigated whether decreased rainfall is associated with adverse outcomes among PWH on ART in Southern Africa.

Methods

Data were combined from 11 clinical cohorts of PWH in Lesotho, Malawi, Mozambique, South Africa, Zambia, and Zimbabwe, participating in the International epidemiology Databases to Evaluate AIDS Southern Africa (IeDEA-SA) collaboration. Adult PWH who had started ART prior to 01/06/2016 and were in follow-up in the year prior to 01/06/2016 were included. Two-year rainfall from June 2014 to May 2016 at the location of each HIV centre was summed and ranked against historical 2-year rainfall amounts (1981–2016) to give an empirical relative percentile rainfall estimate. The IeDEA-SA and rainfall data were combined using each HIV centre’s latitude/longitude. In individual-level analyses, multivariable Cox or generalized estimating equation regression models (GEEs) assessed associations between decreased rainfall versus historical levels and four separate outcomes (mortality, CD4 counts < 200 cells/mm3, viral loads > 400 copies/mL, and > 12-month gaps in follow-up) in the two years following the rainfall period. GEEs were used to investigate the association between relative rainfall and monthly numbers of unique visitors per HIV centre.

Results

Among 270,708 PWH across 386 HIV centres (67% female, median age 39 [IQR: 32–46]), lower rainfall than usual was associated with higher mortality (adjusted Hazard Ratio: 1.18 [95%CI: 1.07–1.32] per 10 percentile rainfall rank decrease) and unsuppressed viral loads (adjusted Odds Ratio: 1.05 [1.01–1.09]). Levels of rainfall were not strongly associated with CD4 counts < 200 cell/mm3 or > 12-month gaps in care. HIV centres in areas with less rainfall than usual had lower numbers of PWH visiting them (adjusted Rate Ratio: 0.80 [0.66–0.98] per 10 percentile rainfall rank decrease).

Conclusions

Decreased rainfall could negatively impact on HIV treatment behaviours and outcomes. Further research is needed to explore the reasons for these effects. Interventions to mitigate the health impact of severe weather events are required.

Peer Review reports

Introduction

Climate change is increasing the frequency and severity of extreme weather events, including heatwaves, flooding, cyclones, and meteorological droughts [1], defined as an exceptional lack of rainfall compared to normal circumstances [2]. These severe weather anomalies impact human health through various mechanisms, including altered patterns of vector-borne and water-borne diseases linked to drought [3], and heat-related illnesses [4]. Southern Africa is one of the regions most affected by climate change [1] due to increasing risks of drought caused by changes in precipitation and limited water storage, as well as limited capacity and resources to support adaptation [1]. Drought is an ongoing and worsening trend [5]. The fraction of sub-Saharan Africa experiencing severe drought increased from < 5% in 1901 to ~ 15% in 2013 [6], and is projected to further increase in this century [7].

Climate change has been hypothesised by UNAIDS to impact the HIV epidemic in settings where HIV prevalence is highest [8], such as Southern Africa [9]. A lack of rainfall impacts vegetation, including agricultural production [10]. This can lead to increased food insecurity and poverty, which in turn can lead to behaviours that result in sub-optimal HIV treatment outcomes through a variety of mechanisms [11, 12]. For instance, evidence suggests that the resulting income shocks caused by drought and food insecurity can disrupt antiretroviral therapy (ART) treatment schedules [13]. This may occur through reduced adherence [14] due to people finding it difficult to travel to clinics because they have to prioritise food and money over attending clinics for ART [13], potentially leading to decreases in CD4 cell counts [15] and increases in viral load [16]. Additionally, some ART medications are more effective when taken with food [17], while some can be taken with or without food. Unfortunately, research into this field has so far been limited because of a lack of cross-disciplinary expertise, difficulties in identifying appropriate data, and complexity in the causal pathway between climate exposures and HIV-related outcomes.

Using a large, established multi-country HIV clinical cohort in Southern Africa, linked with spatial and temporal data on levels of rainfall, we investigated whether living somewhere that has recently had levels of rainfall lower than historical averages is associated with subsequent adverse treatment outcomes among persons with HIV (PWH) on ART.

Methods

Cohort data

The International epidemiology Databases to Evaluate AIDS (IeDEA) (https://www.iedea.org/) is an international research consortium that collects deidentified patient-level data from approximately two million PWH across 46 countries [18]. IeDEA Southern Africa [IeDEA-SA] (https://www.iedea.org/regions/southern-africa/), one of four African IeDEA regions, comprises ART programs that collect data from facilities across Lesotho, Malawi, Mozambique, South Africa, Zambia and Zimbabwe [19]. Local review boards and ethics committees approved the use of IeDEA data for research within the IeDEA collaboration. The Ethics Committee of the Canton of Bern (150/14, PB 2016–00273), Switzerland, approved data merging and collaborative analyses of the IeDEA-SA data (the methods of individual IeDEA-SA concept sheets are not reviewed by an ethics committee). Informed consent for the use of IeDEA routinely collected data has been obtained or waived according to local requirements of each cohort. For this analysis, data were available for 11 IeDEA-SA HIV cohorts from all included countries that had longitude and latitude coordinates available for their participating HIV centres.

Rainfall data

We used longitudinal, gridded rainfall data to develop a measure to capture weather shocks compared with a historical norm for each location in our study. Data on rainfall estimates from the Climate Hazards Group InfraRed Precipitation with Station Data (CHIRPS) at 0.05° resolution (roughly 30km2) were used to quantify rainfall [20]. This gridded dataset was prepared by the Vulnerability Analysis and Mapping (VAM) Geospatial Analysis Team at the Analysis and Trends Service of the World Food Programme (WFP). The 2-year total rainfall from June 2014 to May 2016 for each grid location was summed and then ranked compared to all 2-year rainfall amounts within the 1981–2016 period for that grid location and converted to an empirical relative percentile rank. This produced a variable with values from 1–100, with 50 indicating median rainfall versus historical levels, and lower values indicating less rainfall than usual [21]. We combined the IeDEA-SA HIV cohort data with the gridded data on rainfall using latitude and longitude data of each HIV centre and overlaying them with the gridded rainfall dataset. Further information on the climate contexts in the region during this time-period are given in the supplementary materials. The 2014–2016 period was chosen as Southern Africa experienced a severe drought during this time-period [21].

Individual-level inclusion criteria

PWH on ART were eligible for inclusion if they had started ART prior to 1st June 2016, were aged ≥ 16 years on this date, had a recorded follow-up visit between 1st June 2015 and the 31st May 2016, and had not been recorded by the cohorts as lost-to-follow-up or dead before the 1st June 2016. PWH were dropped from the analyses if they were missing data on which HIV centre they had attended (N = 43,531). A small number of PWH were dropped who attended HIV centres that did not have longitude or latitude coordinates available (N = 169), were recorded as having closed prior to 2016 (N = 32), had missing data on age (N = 28) or sex at birth (N = 24), or started ART before 2000 (N = 11). For each separate analysis, further cohort-level exclusions were made to ensure the included cohorts captured the outcomes of interest. For the analysis of all-cause mortality, one cohort (CIDRZ) was dropped due to suspected under-ascertainment of mortality, determined by a very low percentage of deaths compared to the other cohorts, whilst another (Kheth’Impilo) was dropped due to a lack of data on dates of death. For the analyses of CD4 counts < 200 cells/mm3 or viral load > 400 copies/mL, cohorts were dropped if > 50% of eligible PWH were missing data on the respective outcomes during follow-up. Additionally, for these two analyses, cohorts were dropped if > 50% of eligible PWH were missing a CD4 count or viral load measurement, respectively, at follow-up start. Eight cohorts (Lighthouse Trust, SMART Zimbabwe, Hlabisa, Gugulethu Community Health Centre, Themba Lethu, Khayelitsha, CIDRZ, and SMART Lesotho) were dropped from the CD4 analysis and five were dropped from the viral load analysis (CIDRZ, SMART Lesotho, SMART Mozambique, SMART Zimbabwe, and Lighthouse Trust).

Individual-level analyses

Our analyses investigated whether changes in rainfall between June 2014 – May 2016 (compared with historical values) affected HIV outcomes up to 2 years afterwards. Four outcomes were investigated in individual-level analyses. For each outcome of interest, models adjusted for cohort and evaluated the association of the outcome per 10 percentile relative decreases in rainfall rank (the scale of the rainfall variable was reversed for interpretability). Although we considered using quadratic terms to account for both rainfall extremes (drought and flooding) only one centre in the dataset experienced a relative rainfall value that could be considered extremely high as compared to historical values (≥ 85), so we instead included just a linear term.

The associations between the outcomes of all-cause mortality and gaps in follow-up of ≥ 12 months with the variable on the relative change in rainfall were analysed using Cox proportional hazards models. For analyses where CD4 counts < 200 cells/mm3 or viral load > 400 copies/mL were the (binary) outcomes, generalized estimating equation (GEE) regression models were fitted (binomial, with a logit link) accounting for each recorded CD4 or viral load value (up to the 1st June 2018) and additionally adjusting for the time between follow-up start and the date of each recorded value.

Follow-up started on 1st June 2016 and ended at the earliest of either the cohort/centre-specific administrative censoring date, 1st June 2018, or, for the analyses of all-cause mortality and gaps in follow-up of > 1 year, the date of this outcome. The end date 1st June 2018 was chosen to limit the period in which the outcomes could be recorded to 2 years following the drought period. A shorter follow-up period was investigated in sensitivity analyses. If a person had a gap of over a year between their last recorded visit date and the cohort/centre-specific administrative censoring date, then they were considered as lost-to-follow-up one year after their last recorded visit date.

Models additionally adjusted for sex at birth, and data recorded at follow-up start on age (16–25, 26–35, 36–45, 46–55, 56–65, ≥ 66), CD4 count cells/mm3 (0–99, 100–199, 200–349, 350–499, ≥ 500, missing), viral load in copies/mL (< 400, ≥ 400, missing), AIDS status (no AIDS, AIDS, unknown – using WHO staging), and the time since starting ART (< 6 months, 6 months to 1 year, 1–3 years, 3–6 years, 6–10 years, and ≥ 10 years). The models were also adjusted for a binary variable to indicate locations with a Normalized Difference Vegetation Index (NDVI) value ≥ 0.3. NDVI assesses whether a pixel contains live green vegetation and is both an indicator of plant health and land use (https://gisgeography.com/ndvi-normalized-difference-vegetation-index/). NDVI values were taken on 1st Jun 2016 (when analyses started). Low values indicate a place has less vegetation. These data were linked to the IeDEA-SA data using the longitude and latitude of the HIV centres. To capture the values for CD4 and viral loads at follow-up start, a window from 1st January 2015 to 31st May 2016 was used, taking the value closest to follow-up start. Models were run including PWH from all HIV centres, as well as stratifying by rural/urban status of each HIV centre.

Sensitivity analyses

For analyses of each outcome, sensitivity analyses were performed removing one cohort at a time and limiting follow-up time at 1 year instead of 2 years. For the analysis of mortality, two additional sensitivity analyses were performed 1) including the Kheth’Impilo cohort that did not have dates of death available and using logistic regression rather than Cox proportional hazards models because of this lack of data on death dates, and 2) not including the vital registration linkage that was available only for the South African cohorts.

Analysis of visitors per HIV centre

We used generalized estimating equation regression models (Negative binomial, with an identity link) to investigate whether relative rainfall impacted the number of PWH attending HIV centres between June 2014 and May 2016. HIV centres from participating cohorts were included if they had at least 20 unique visitors in June 2014 and were still open by May 2016. The variables included in the models were the rainfall change variable, the NDVI ≥ 0.3 variable, urban/rural region, the number of months since June 2014, and the cohort. Models were re-run stratifying by rural and urban region. There were no missing data for any covariates included in this analysis.

Results

Data on 843,289 PWH of all ages, on and off ART were available in the IeDEA-SA dataset, which spanned the period 1997–2021. Of these, 43,531 PWH were excluded due to having missing data on which HIV centre they had attended. Of the remainder, 270,708 adult PWH had started ART prior to 01/06/2016 and were still alive at this point and in follow-up in the year prior to this, so were eligible for inclusion in this study. These 270,708 PWH were included from 386 HIV centres (Fig. 1; Table 1), out of 393 in the overall dataset (5 centres had missing geocode coordinates and 2 were not operating at the start of follow-up). 181,722 (67.1%) of the 270,708 PWH were female and the median age was 39 years (interquartile range [IQR]: 32–46) (Table 2). The median empirical percentile of average rainfall for June 2014 to May 2016 versus historical values was 22 (interquartile range [IQR]: 4–31; range: 1–89), where 50 would indicate average rainfall. Of the 11 cohorts included, 6 and 3 were based entirely in urban and rural areas, respectively, with the other 2 containing a mix of urban and rural areas. For the 222,865 (82.3%) PWH that attended HIV centres in urban areas, the median rainfall percentile was 25 (IQR: 4–31; range: 1–54), whilst it was 10 (IQR: 7–25; range: 1–89) for the 47,843 (17.7%) PWH that attended HIV centres in rural areas. The number of PWH attending HIV centres in regions with an NDVI value ≥ 0.3 was 68,208 (25.2%). At the start of follow-up, the median length of time that the included persons had been on ART was 3.5 (IQR: 1.4, 6.3) years and 48,071 (17.8%) had experienced a prior AIDS event. Among the 154,091 (56.9%) PWH with data on CD4 cell count at follow-up start, the median value was 425 (IQR: 270–600) cells/mm3.

Fig. 1
figure 1

Locations of the HIV centres (red dots) and relative rainfall rank for June 2014 to May 2016 versus historical levels in an observational cohort study of the effect of rainfall on HIV outcomes in Southern Africa

Full size image
Table 1 Characteristics of each HIV cohort in an observational cohort study of the effect of rainfall on HIV outcomes in Southern Africa
Full size table
Table 2 Time-updated characteristics at follow-up start (1st June 2016) of included PWH, stratified by rural or urban HIV centre location in an observational cohort study of the effect of rainfall on HIV outcomes in Southern Africa
Full size table

Mortality

For the analyses examining mortality, 104,138 PWH were included from 9 cohorts, among whom there were 2,951 deaths (2.8%) in 327,414 person-years. In the overall analysis (Table 3; Supplementary Table 1), decreases in relative rainfall levels versus historical values were associated with higher mortality, adjusted hazard ratio (aHR) 1.18 (95% confidence interval [95%CI]: 1.07–1.32) per 10 percentile decrease in rainfall. In analyses of PWH in rural sites, these results persisted. For PWH in urban sites, there was no association between mortality and changes in rainfall due to a lack of within-cohort variation in rainfall for the subgroup included in the analysis. The association persisted across 12/13 sensitivity analyses (Supplementary Table 2), although the confidence interval was much wider when dropping the SMART Mozambique cohort, aHR 1.08 (95%CI: 0.78–1.49). The association was stronger when using 1-year follow-up after the drought period, rather than 2 years, aHR 1.36 (95%CI 1.14–1.61).

Table 3 Adjusteda hazard ratios and odds ratios for each outcome measure per decrease of 10 percentiles in rainfall versus historical values in an observational cohort study of the effect of rainfall on HIV outcomes in Southern Africa
Full size table

CD4 counts < 200 cells/mm.3

When investigating CD4 counts < 200 cells/mm3 as the outcome, 27,580 PWH were included from 3 cohorts, among whom 3,990 (14.5%) had a CD4 count < 200 cells/mm3 recorded during follow-up (Table 3; Supplementary Table 3). In the overall analysis, there was weak evidence of a negative association between decreases in relative rainfall levels versus historical values and having a CD4 count < 200 cells/mm3, adjusted odds ratio (aOR) 0.94 (95%CI: 0.89–1.00) per 10% percentile decrease in the rainfall level. When stratifying by rural/urban location, the association was very similar in urban areas and was much attenuated in rural areas, aOR 0.99 (0.89–1.10). Results for 5/5 of the sensitivity analyses were similar to those in the main analysis (Supplementary Table 4).

HIV-1 viral loads > 400 copies/mL

For the analysis with HIV-1 viral load > 400 copies/mL as the outcome, 6 cohorts were included containing 82,260 PWH, among whom 13,788 (16.8%) had viral loads > 400 copies/mL recorded during follow-up (Table 3; Supplementary Table 5). In the overall analysis, lower relative rainfall levels were positively associated with having a viral load > 400 copies/mL recorded, aOR 1.05 (95%CI: 1.01–1.09) per 10 percentile decrease in rainfall. This analysis could not be repeated in rural sites due to a lack of data on this outcome, but the results persisted for PWH in urban areas. In 7/8 of the sensitivity analyses, similar odds ratios and p-values were seen (Supplementary Table 6), but the results attenuated substantially when removing the Kheth’Impilo cohort, which is driving the overall association (there was no within-cohort variation in rainfall among the remaining cohorts).

 ≥ 12-month gaps in care

All 11 cohorts and 270,708 PWH were included in the analyses of ≥ 12-month gaps in care (Table 3; Supplementary Table 7). There were 9,426 PWH who had a ≥ 12-month gap in 567,272 person-years. Overall, there was no evidence of an association between relative rainfall levels versus historical values and ≥ 12-month gaps in care, aHR 0.98 (95%CI: 0.94–1.01), per 10% decrease. However, when stratifying by urban and rural location, there was evidence of a protective association between lower relative rainfall and ≥ 12-month gaps in care in both urban and rural areas, aORs 0.96 (0.92–1.00) and 0.86 (0.80–0.93), respectively, per 10 percentile decrease in rainfall rank. In the 12 sensitivity analyses, the results remained consistent with the overall analysis (supplementary Table 8).

Analysis of visitors per HIV centre

For the analysis of numbers of unique PWH visiting HIV centres per month, 187 HIV centres were included. The median number of unique PWH visiting each HIV centre per month increased from 196 in June 2014 to 220 in May 2016 (Supplementary Fig. 1). Each December saw a drop in urban areas. The numbers of visitors were higher in urban than in rural areas. Table 4 shows the associations of relative rainfall levels with monthly numbers of unique visitors per HIV centre. In fully adjusted analyses, HIV centres in areas with lower levels of rainfall than usual had lower numbers of PWH visiting them, with the number of unique visitors reducing by a factor of 0.80 (95%CI: 0.66–0.98) per 10% decrease in rainfall. There was a positive association between the number of visitors and the number of months since January 2014, i.e. the clinics became busier over time. In both the rural and urban regions, the point estimate for the association between relative rainfall levels and monthly clinic visitors was similar to the overall analysis, but confidence intervals were wider: numbers of unique visitors decreased by factors of 0.80 (0.53–1.19) and 0.83 (0.66–1.05) for rural and urban areas, respectively, per 10 percentile decrease in rainfall rank.

Table 4 Adjusted rate ratios for the numbers of unique visitors per HIV centre for each month between June 2014 and May 2016 in an observational cohort study of the effect of rainfall on HIV outcomes in Southern Africa
Full size table

Discussion

We found higher mortality among PWH on ART living in regions with lower rainfall than usual, both overall and when including just rural areas. We also found higher odds of unsuppressed viral loads among PWH on ART in regions with lower-than-normal rainfall, overall and in urban areas. However, we did not find convincing evidence of relationships between decreases in rainfall and recording low CD4 counts or ≥ 12-month gaps in care among PWH on ART. Additionally, we found that HIV centres in areas with less rainfall than usual had lower numbers of PWH visiting them. Overall, among PWH on ART we saw some evidence of negative HIV treatment outcomes in areas with lower rainfall than historically. However, the evidence was not consistent across outcomes or when stratifying by rural/urban location, and, therefore, the pathways between decreases in rainfall and HIV outcomes require further elucidation.

Comparisons with other literature

Climate change is hypothesised to negatively impact on HIV treatment outcomes via increases in food insecurity and poverty [8], although evidence on this is limited. There have been various studies that have established that food insecurity, which is linked to poverty, can adversely impact on HIV treatment outcomes [22,23,24]. This includes systematic reviews on the relationships between food insecurity and adherence to ART [13], HIV-1 viral suppression [16], and CD4 count [15]. However, to our knowledge, this is the first study to look at the effect of drought or rainfall levels on ART outcomes including mortality. Otherwise, previous cross-sectional studies of weather and HIV in sub-Saharan Africa have found associations between local rainfall shocks and heightened HIV prevalence [25], as well as between drought and increased HIV prevalence and riskier sexual behaviour in young women in rural areas [26]. A recent study in South Africa found that adherence and retention in care decreased during years of drought [14]. Another study found evidence that unusually heavy rainfall was associated with a higher HIV burden across 21 sub-Saharan African countries [27]. Other analyses have also found that droughts increased transactional sex among women employed in agriculture in Malawi, as well as HIV prevalence among both men and women [28]. Finally, a study across 10 sub-Saharan African countries found that droughts were associated with lower odds of HIV testing and higher odds of condomless sex [29]. That much of this literature found associations between droughts and HIV-related outcomes among women living in rural areas show both the gender and urban/rural divides in the pathways linking drought and HIV. Studies suggest that the impacts of climate on HIV are seen particularly among women because of greater dependence on transactional sex [30], and that the effects are concentrated in rural areas, due to changes in climate having very direct impacts on the livelihoods and behaviours of farmers, whilst any affects in urban areas would likely be less direct, potentially occurring through drought affecting the general economy [14]. Other literature has looked at the effects of drought in sub-Saharan Africa on other health outcomes that, like HIV, would be more indirectly affected by climate. This includes studies suggesting that childhood vaccination levels are reduced in the presence of drought, potentially due to drought creating extra barriers to clinic access [31].

Strengths and limitations

A strength of this analysis is the use of a very large, longitudinal cohort of PWH on ART spanning multiple Southern African countries combined with a measure of change in rainfall. This has allowed us to analyse outcomes that occurred after these exposures were measured, removing the difficulties of interpreting observed associations in cross-sectional data. However, there are several limitations. Although in each analysis we adjusted for cohort, the associations seen between lower rainfall and HIV treatment outcomes could instead be capturing other differences between regions. The cohorts capture different populations and some of the associations seen could be reflecting cross-cohort and cross-country differences in unmeasured confounders, such as types of HIV care provider, treatment practices, income, and food insecurity. Missing data on potential confounders is a common issue with routinely collected data. Furthermore, there is differential reporting and recording of variables between cohorts, including of CD4 cell counts and viral loads, and in the analyses of mortality, linkage to vital registration was only possible for some of the South African cohorts. In particular, few cohorts collected robust data on CD4 count monitoring due to a transition to viral load monitoring around this time [19]. Therefore, the analysis of CD4 counts only contained data from 3 cohorts. The inconsistency of the results across the various outcome measures (CD4 counts, viral loads, gaps in care, and mortality) could be due to the differential recording of these variables, particularly in the South African cohorts where there was linkage to death registry data, so this will have been recorded more fully than the HIV biomarkers. To account for these various limitations, we performed many sensitivity analyses, including removing the vital registration linkage for the South African cohorts, and investigating how leaving out each cohort from each analysis affected the results. Additionally, several cohorts only had one or two centres, meaning that there was no within-cohort variation in rainfall. This meant that the cohorts with more centres were more influential in the analyses and further limited the number of cohorts included and variation in rainfall that was captured when stratifying analyses by rural and urban areas. Potential mechanisms between rainfall and HIV epidemiology are less clear in areas where there is little subsistence farming, particularly urban areas (which often draw their water from dams located elsewhere), so, any associations in analyses specific to urban areas should be interpreted with caution. We also do not know the timespan in which rainfall levels will affect treatment outcomes – this requires further research. However, we investigated using a shorter 1-year follow-up duration in sensitivity analyses and found similar results for low CD4 counts, unsuppressed viral loads, and 12-month gaps in care, and a stronger association between decreased rainfall and increased mortality than when using 2 years of follow-up. For the analysis of unique visitor numbers, this could have been affected by lengthening ART dispensing periods over time, although this would likely have been happening across all HIV cohorts and centres. As data were only available on the locations of the HIV centres, rather than on place of residence, we had to assume that their place of residence had comparable weather to their HIV centre. When analysing mortality, we were only able to look at all-cause mortality rather than cause-specific mortality, so were unable to understand whether the increased mortality in regions with less rainfall than usual was AIDS-related or due to other causes.

Conclusion

Our results add to the evidence that changes in rainfall can be associated with the epidemiology of HIV in high HIV prevalence areas such as Southern Africa through changing behaviours [28, 29]. Previous research has highlighted that food insecurity and income are key links in the chain between climate change and HIV [12]. Additional research is required using longitudinal cohort data where information on these links can be accounted for to further our understanding. Improvements in HIV care could mitigate the potential future negative effects of reduced rainfall on HIV treatment outcomes, particularly in rural areas. However, the frequency and severity of droughts are projected to increase in Southern Africa [7], which could blunt the impact of improvements in HIV care on treatment outcomes, hindering progress towards UNAIDS’ goals of controlling the HIV epidemic by 2030 [32]. Additionally, climate change could cause increased conflict over land and resources, with such conflicts also likely to negatively impact on the HIV care cascade [33]. Our results regarding the effect of drought on HIV treatment outcomes, also have implications for the monitoring and treatment of other chronic conditions that may worsen in drought situations, although reviews have noted a surprising lack of research into this topic [34, 35]. Whilst human-driven climate change is a global issue that will require global efforts and investment to deal with [1], environmental and HIV interventions could be considered in high HIV burden areas that are likely to experience droughts. Localised interventions in Southern Africa that could mitigate the effects of climate change on HIV outcomes could include the use of newly developed long-lasting antiretrovirals that do not need to be taken daily [36], multi-month dispensing [37], microcredit interventions to reduce reliance on subsistence agriculture [38], reducing the travel costs and time associated with obtaining ART [39, 40], and building early warning systems to help health systems anticipate extreme weather conditions [41].

Availability of data and materials

The data file designed by the World Food Programme's Vulnerability Analysis and Mapping Geospatial Analysis Team that compares historical rainfall patterns with rainfall between June 2014 to May 2016 is available on request—please contact the corresponding author for access as the datasets are very large. The HIV cohort data are not publicly available. Please contact IeDEA-SA for further information (https://www.iedea-sa.org/).

References

  1. Trisos CH, Adelekan IO, Totin E, Ayanlade A, Efitre J, Gemeda A, et al. 2022: Africa. Climate Change 2022: Impacts, Adaption and Vulnerability. Cambridge, UK and New York, NY, USA: Cambridge University Press; 2022. p. 1285–455.

    Google Scholar 

  2. Van Loon AF, Gleeson T, Clark J, Van Dijk AIJM, Stahl K, Hannaford J, et al. Drought in the Anthropocene. Nat Geosci. 2016;9(2):89–91.

    Article  Google Scholar 

  3. Campbell-Lendrum D, Manga L, Bagayoko M, Sommerfeld J. Climate change and vector-borne diseases: what are the implications for public health research and policy? Philos T R Soc B. 2015;370(1665):20130552.

    Article  Google Scholar 

  4. Morefield PE, Fann N, Grambsch A, Raich W, Weaver CP. Heat-Related Health Impacts under Scenarios of Climate and Population Change. Int J Env Res Pub He. 2018;15(11):2438.

  5. NASA. Earth Observatory. Drought Threatens Millions in Southern Africa. 2019.

    Google Scholar 

  6. Masih I, Maskey S, Mussa FEF, Trambauer P. A review of droughts on the African continent: a geospatial and long-term perspective. Hydrol Earth Syst Sc. 2014;18(9):3635–49.

    Article  Google Scholar 

  7. Ahmadalipour A, Moradkhani H, Castelletti A, Magliocca N. Future drought risk in Africa: Integrating vulnerability, climate change, and population growth. Sci Total Environ. 2019;662:672–86.

    Article  PubMed  CAS  Google Scholar 

  8. UNAIDS. Climate Change and AIDS: A Joint Working Paper. 2008.

    Google Scholar 

  9. UNAIDS. Seizing the moment: Global AIDS update 2020. Geneva, Switzerland: UNAIDS; 2020.

  10. Ngcamu BS, Chari F. Drought Influences on Food Insecurity in Africa: A Systematic Literature Review. Int J Environ Res Public Health. 2020;17(16):5897.

  11. Chop E, Duggaraju A, Malley A, Burke V, Caldas S, Yeh PT, et al. Food insecurity, sexual risk behavior, and adherence to antiretroviral therapy among women living with HIV: A systematic review. Health Care Women In. 2017;38(9):927–44.

    Article  Google Scholar 

  12. Orievulu KS, Ayeb-Karlsson S, Ngema S, Baisley K, Tanser F, Ngwenya N, et al. Exploring linkages between drought and HIV treatment adherence in Africa: a systematic review. Lancet Planet Health. 2022;6(4):359–70.

    Article  Google Scholar 

  13. Singer AW, Weiser SD, Mccoy SI. Does Food Insecurity Undermine Adherence to Antiretroviral Therapy? A Systematic Review Aids Behav. 2015;19(8):1510–26.

    Article  PubMed  Google Scholar 

  14. Iwuji CC, Baisley K, Maoyi ML, Orievulu K, Mazibuko L, Ayeb-Karlsson S, et al. The Impact of Drought on HIV Care in Rural South Africa: An Interrupted Time Series Analysis. Ecohealth. 2023;20(2):178–93.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Aibibula W, Cox J, Hamelin AM, Mamiya H, Klein MB, Brassard P. Food insecurity and low CD4 count among HIV-infected people: a systematic review and meta-analysis. AIDS Care. 2016;28(12):1577–85.

    Article  PubMed  Google Scholar 

  16. Aibibula W, Cox J, Hamelin AM, McLinden T, Klein MB, Brassard P. Association Between Food Insecurity and HIV Viral Suppression: A Systematic Review and Meta-Analysis. Aids Behav. 2017;21(3):754–65.

    Article  PubMed  Google Scholar 

  17. Young S, Wheeler AC, McCoy SI, Weiser SD. A Review of the Role of Food Insecurity in Adherence to Care and Treatment Among Adult and Pediatric Populations Living with HIV and AIDS. Aids Behav. 2014;18:S505–15.

    Article  PubMed  Google Scholar 

  18. Chammartin F, Ostinelli CH, Anastos K, Jaquet A, Brazier E, Brown S, et al. International epidemiology databases to evaluate AIDS (IeDEA) in sub-Saharan Africa, 2012–2019. Bmj Open. 2020;10(5):e035246.

  19. Zaniewski E, Brazier E, Ostinelli CH, Wood R, Osler M, Technau KG, et al. Regression discontinuity analysis demonstrated varied effect of Treat-All on CD4 testing among Southern African countries. J Clin Epidemiol. 2021;140:101–10.

    Article  PubMed  PubMed Central  Google Scholar 

  20. Funk C, Peterson P, Landsfeld M, Pedreros D, Verdin J, Shukla S, et al. The climate hazards infrared precipitation with stations-a new environmental record for monitoring extremes. Sci Data. 2015;2:150066.

  21. World Food Programme. Southern Africa growing season 2016–2017: Recovery after two years of drought?. Rome; 2016. https://documents.wfp.org/stellent/groups/public/documents/ena/wfp290431.pdf.

  22. Weiser SD, Tsai AC, Gupta R, Frongillo EA, Kawuma A, Senkungu J, et al. Food insecurity is associated with morbidity and patterns of healthcare utilization among HIV-infected individuals in a resource-poor setting. AIDS. 2012;26(1):67–75.

    Article  PubMed  Google Scholar 

  23. Weiser SD, Palar K, Frongillo EA, Tsai AC, Kumbakumba E, dePee S, et al. Longitudinal assessment of associations between food insecurity, antiretroviral adherence and HIV treatment outcomes in rural Uganda. AIDS. 2014;28(1):115–20.

    Article  PubMed  CAS  Google Scholar 

  24. Haacker M, Birungi C. Poverty as a barrier to antiretroviral therapy access for people living with HIV/AIDS in Kenya. Ajar-Afr J Aids Res. 2018;17(2):145–52.

    Article  Google Scholar 

  25. Burke M, Gong E, Jones K. Income Shocks and HIV in Africa. Econ J. 2015;125(585):1157–89.

    Article  Google Scholar 

  26. Low AJ, Frederix K, McCracken S, Manyau S, Gummerson E, Radin E, et al. Association between severe drought and HIV prevention and care behaviors in Lesotho: A population-based survey 2016–2017. Plos Med. 2019;16(1):e1002727

  27. Nagata JM, Hampshire K, Epstein A, Lin F, Zakaras J, Murnane P, et al. Analysis of Heavy Rainfall in Sub-Saharan Africa and HIV Transmission Risk, HIV Prevalence, and Sexually Transmitted Infections, 2005–2017. Jama Netw Open. 2022;5(9):e2230282.

  28. Treibich C, Bell E, Blanc E, Lepine A. From a drought to HIV: An analysis of the effect of droughts on transactional sex and sexually transmitted infections in Malawi. Ssm-Popul Hlth. 2022;19:101221.

  29. Epstein A, Nagata JM, Ganson KT, Nash D, Saberi P, Tsai AC, et al. Drought, HIV Testing, and HIV Transmission Risk Behaviors: A Population-Based Study in 10 High HIV Prevalence Countries in Sub-Saharan Africa. Aids Behav. 2022;27(3):855–63.

    Article  PubMed  Google Scholar 

  30. Austin KF, Noble MD, Berndt VK. Drying Climates and Gendered Suffering: Links Between Drought, Food Insecurity, and Women’s HIV in Less-Developed Countries. Soc Indic Res. 2021;154(1):313–34.

  31. Nagata JM, Epstein A, Ganson KT, Benmarhnia T, Weiser SD. Drought and child vaccination coverage in 22 countries in sub-Saharan Africa: A retrospective analysis of national survey data from 2011 to 2019. Plos Med. 2021;18(9):e1003678.

  32. UNAIDS. Fast-track: Ending the AIDS epidemic by 2030. Geneva, Switzerland: UNAIDS; 2014.

    Google Scholar 

  33. Mann M, Diero L, Kemboi E, Mambo F, Rono M, Injera W, et al. Antiretroviral Treatment Interruptions Induced by the Kenyan Postelection Crisis Are Associated With Virological Failure. JAIDS. 2013;64(2):220–4.

  34. Stanke C, Kerac M, Prudhomme C, Medlock J, Murray V. Health effects of drought: a systematic review of the evidence. PLoS Curr. 2013;5:5:ecurrents.dis.7a2cee9e980f91ad7697b570bcc4b004.

  35. Asmall T, Abrams A, Roosli M, Cisse G, Carden K, Dalvie MA. The adverse health effects associated with drought in Africa. Sci Total Environ. 2021;793:148500.

  36. Bares SH, Scarsi KK. A new paradigm for antiretroviral delivery: long-acting cabotegravir and rilpivirine for the treatment and prevention of HIV. Curr Opin Hiv Aids. 2022;17(1):22–31.

    Article  PubMed  CAS  Google Scholar 

  37. Lopes J, Grimwood A, Ngorima-Mabhena N, Tiam A, Tukei BB, Kasu T, et al. Out-of-Facility Multimonth Dispensing of Antiretroviral Treatment: A Pooled Analysis Using Individual Patient Data From Cluster-Randomized Trials in Southern Africa. JAIDS. 2021;88(5):477–86.

  38. Linnemayr S, Buzaalirwa L, Balya J, Wagner G. A Microfinance Program Targeting People Living with HIV in Uganda: Client Characteristics and Program Impact. J Int Assoc Provid AIDS Care. 2017;16(3):254–60.

    Article  PubMed  Google Scholar 

  39. Lankowski AJ, Siedner MJ, Bangsberg DR, Tsai AC. Impact of Geographic and Transportation-Related Barriers on HIV Outcomes in Sub-Saharan Africa: A Systematic Review. Aids Behav. 2014;18(7):1199–223.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Tuller DM, Bangsberg DR, Senkungu J, Ware NC, Emenyonu N, Weiser SD. Transportation Costs Impede Sustained Adherence and Access to HAART in a Clinic Population in Southwestern Uganda: A Qualitative Study. Aids Behav. 2010;14(4):778–84.

    Article  PubMed  Google Scholar 

  41. Ebi KL, Vanos J, Baldwin JW, Bell JE, Hondula DM, Errett NA, et al. Extreme Weather and Climate Change: Population Health and Health System Implications. Annu Rev Publ Health. 2021;42:293–315.

    Article  Google Scholar 

Download references

Acknowledgements

We would like to thank all those who contributed data or were involved in data collection.

AT acknowledges funding from Wellcome Trust (222770/Z/21/Z). PV acknowledges funding from NIAID and NIDA (R01AI147490, R01DA033679 and R21DA047902), the Wellcome Trust (WT 226619/Z/22/Z), and support from the UK National Institute for Health and Care Research (NIHR) Health Protection Research Unit (HPRU) in Behavioural Science. TW acknowledges funding from the Alexander von Humboldt Foundation in the framework of the Alexander von Humboldt Professorship endowed by the German Federal Ministry of Education and Research.

Research reported in this publication was supported by the U.S. National Institutes of Health’s National Institute of Allergy and Infectious Diseases (NIAID), the Eunice Kennedy Shriver National Institute of Child Health and Human Development (NICHD), the National Cancer Institute (NCI), the National Institute on Drug Abuse (NIDA), the National Heart, Lung, and Blood Institute (NHLBI), the National Institute on Alcohol Abuse and Alcoholism (NIAAA), the National Institute of Diabetes and Digestive and Kidney Diseases (NIDDK) and the Fogarty International Center (FIC) under Award Number U01AI069924. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.

Site investigators and cohorts:

Gary Maartens, Aid for AIDS, South Africa; Carolyn Bolton, Center for Infectious Disease Research in Zambia (CIDRZ), Zambia; Robin Wood, Gugulethu (Desmond Tutu HIV Centre), South Africa; Nosisa Sipambo, Harriet Shezi Children’s Clinic, South Africa; Frank Tanser, Collins Iwuji, Hlabisa (Africa Health Research Institute), South Africa; Jonathan Euvrard, Khayelitsha ART Programme, South Africa; Geoffrey Fatti, Kheth’Impilo AIDS Free Living, South Africa; Safari Mbewe, Lighthouse Trust, Malawi; Mazvita Muchengeti, National Cancer Registry (National Health Laboratory Service), South Africa; Cleophas Chimbetete, Newlands Clinic (Ruedi Luethy Foundation Zimbabwe), Zimbabwe; Karl Technau, Rahima Moosa Mother and Child Hospital, South Africa; Brian Eley, Red Cross War Memorial Children’s Hospital, South Africa; Irene Ayakaka, SolidarMed Lesotho; Idivino Rafael, SolidarMed Mozambique; Cordelia Kunzekwenyika, SolidarMed Zimbabwe, Matthew P Fox, Themba Lethu Clinic, South Africa; Hans Prozesky, Tygerberg Hospital, South Africa; Andrew Boulle, Western Cape Provincial Health Data Centre.

Data centers:

John Andoh, Marie Ballif, Nicolas Banholzer, Cam Ha Dao Ostinelli, Matthias Egger, Lukas Fenner, Nathalie Fernandez, Andreas Haas, Eliane Rohner, Julien Riou, Tiana Schwab, Veronika W Skrivankova, Lilian Smith, Katayoun Taghavi, Per von Groote, Gilles Wandeler, Elizabeth Zaniewski, Institute of Social and Preventive Medicine, University of Bern, Switzerland; Nanina Anderegg, Kim Anderson, Andrew Boulle, Chido Chinogurei, Morna Cornell, Mary-Ann Davies, Leigh Johnson, Reshma Kassanjee, Amohelang Lehloa Nicola Maxwell, Haroon Moolla, Carl Morrow, Patience Nyakato, Gem Patten, Mpho Tlali, Renee de Waal, Wendy Wiemers, Center for Infectious Disease Epidemiology and Research, School of Public Health, University of Cape Town, South Africa.

Funding

This research was funded in whole, or in part, by the Wellcome Trust (222770/Z/21/Z). For the purpose of Open Access, the author has applied a CC BY public copyright licence to any Author Accepted Manuscript version arising from this submission.

Author information

Authors and Affiliations

  1. Population Health Sciences, University of Bristol, Bristol, UK

    Adam Trickey & Peter Vickerman

  2. Centre for Infectious Disease Epidemiology and Research, School of Public Health and Family Medicine, University of Cape Town, Cape Town, South Africa

    Leigh F. Johnson & Jonathan Euvrard

  3. Department of Civil Engineering, University of Bristol, Bristol, UK

    Fai Fung

  4. UK Meteorological Office, Exeter, UK

    Fai Fung

  5. Climate and Earth Observation Unit, Research Assessment and Monitoring Division, World Food Programme HQ, Rome, Italy

    Rogerio Bonifacio

  6. Africa Health Research Institute, KwaZulu-Natal, South Africa

    Collins Iwuji

  7. Department of Global Health Infection, Brighton and Sussex Medical School, University of Sussex, Brighton, UK

    Collins Iwuji

  8. ICAP at Columbia University, Nakasero, Kampala, Uganda

    Samuel Biraro

  9. Centre for Infectious Disease Research in Zambia, Lusaka, Zambia

    Samuel Bosomprah & Guy Muula

  10. Department of Biostatistics, School of Public Health, University of Ghana, Legon, Accra, Ghana

    Samuel Bosomprah

  11. Newlands Clinic, Harare, Zimbabwe

    Linda Chirimuta

  12. Kheth’Impilo AIDS Free Living, Cape Town, South Africa

    Geoffrey Fatti

  13. Division of Epidemiology and Biostatistics, Department of Global Health, Faculty of Medicine and Health Sciences, Stellenbosch University, Cape Town, South Africa

    Geoffrey Fatti

  14. Health Economics and Epidemiology Research Office, Faculty of Health Sciences, University of the Witwatersrand, Johannesburg, South Africa

    Matthew P. Fox

  15. Department of Global Health and Department of Epidemiology, Boston University School of Public Health, Boston, MA, USA

    Matthew P. Fox

  16. Institute of Social and Preventive Medicine, University of Bern, Bern, Switzerland

    Per Von Groote

  17. Lighthouse Trust, Mzimba, Malawi

    Joe Gumulira

  18. Department of Civil Engineering and Cabot Institute of the Environment, University of Bristol, Bristol, UK

    Guy Howard

  19. Desmond Tutu Health Foundation, Institute of Infectious Diseases and Molecular Medicine, Department of Medicine, University of Cape Town, Cape Town, South Africa

    Lauren Jennings

  20. Research Division, African Population and Health Research Center, Nairobi, Kenya

    Agnes Kiragga

  21. Centre for Epidemic Response and Innovation, School of Data Science and Computational Thinking, Stellenbosch University, Stellenbosch, South Africa

    Frank Tanser

  22. School of Nursing and Public Health, University of KwaZulu-Natal, Durban, South Africa

    Frank Tanser

  23. Institute of Environmental Science and Geography, University of Potsdam, Potsdam, Germany

    Thorsten Wagener

  24. Department of Epidemiology, Mailman School of Public Health, Columbia University, New York, NY, USA

    Andrea Low

  25. NIHR Health Protection Research Unit in Behavioural Science and Evaluation at University of Bristol, Bristol, UK

    Peter Vickerman

Authors
  1. Adam Trickey
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Leigh F. Johnson
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Fai Fung
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Rogerio Bonifacio
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Collins Iwuji
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Samuel Biraro
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. Samuel Bosomprah
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Linda Chirimuta
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Jonathan Euvrard
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. Geoffrey Fatti
    View author publications

    You can also search for this author in PubMed Google Scholar

  11. Matthew P. Fox
    View author publications

    You can also search for this author in PubMed Google Scholar

  12. Per Von Groote
    View author publications

    You can also search for this author in PubMed Google Scholar

  13. Joe Gumulira
    View author publications

    You can also search for this author in PubMed Google Scholar

  14. Guy Howard
    View author publications

    You can also search for this author in PubMed Google Scholar

  15. Lauren Jennings
    View author publications

    You can also search for this author in PubMed Google Scholar

  16. Agnes Kiragga
    View author publications

    You can also search for this author in PubMed Google Scholar

  17. Guy Muula
    View author publications

    You can also search for this author in PubMed Google Scholar

  18. Frank Tanser
    View author publications

    You can also search for this author in PubMed Google Scholar

  19. Thorsten Wagener
    View author publications

    You can also search for this author in PubMed Google Scholar

  20. Andrea Low
    View author publications

    You can also search for this author in PubMed Google Scholar

  21. Peter Vickerman
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

AT conceived and designed the study, with assistance from LJo, PV, and AL. AT wrote the original draft and conducted all statistical analyses. FF, RB, GH, and TW assisted with provision and interpretation of the rainfall and NDVI data. IeDEA-SA cohort representatives (CI, SBo, LC, JE, GF, MF, PVG, LJ, GM, FT, and LJo) contributed to the provision of cohort data. All authors contributed to the interpretation of data and critical revisions of the manuscript for important intellectual content. AT had final responsibility for the decision to submit for publication.

Corresponding author

Correspondence to Adam Trickey.

Ethics declarations

Ethics approval and consent to participate

Local review boards and ethics committees approved the use of IeDEA data for research within the IeDEA collaboration. The Ethics Committee of the Canton of Bern (150/14, PB 2016–00273), Switzerland, approved data merging and collaborative analyses. Informed consent for the use of IeDEA routinely collected data has been obtained or waived according to local requirements.

Consent for publication

Not applicable.

Competing interests

PV reports research grants from Gilead unrelated to this work. CI reports research grants and conference attendance support from Gilead Sciences unrelated to this work. All other authors report no conflicts of interest.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Supplementary table 1.

Full results of the analysis with mortality as the outcome. Supplementary table 2. Sensitivity analysis results for the analysis with mortality as the outcome. Supplementary table 3. Full results for the analysis with CD4 counts<200 cells/mm3 as the outcome. Supplementary table 4. Sensitivity analysis results for the analysis with CD4 counts<200 cells/mm3 as the outcome. Supplementary table 5. Full results for the analysis with viral loads≥400 copies/mL as the outcome. Supplementary table 6. Sensitivity analysis results for the analysis with viral loads≥400 copies/mL as the outcome. Supplementary table 7. Full results for the analysis with 12-month gaps in care as the outcome. Supplementary table 8. Sensitivity analyses results for the analysis with 12-month gaps in care as the outcome. Supplementary table 9. Strengthening the Reporting of Observational studies in Epidemiology (STROBE) checklist. Supplementary Figure 1. Median unique PWH visiting each HIV centre per month.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Trickey, A., Johnson, L.F., Fung, F. et al. Associations of inter-annual rainfall decreases with subsequent HIV outcomes for persons with HIV on antiretroviral therapy in Southern Africa: a collaborative analysis of cohort studies. BMC Infect Dis 23, 889 (2023). https://doi.org/10.1186/s12879-023-08902-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12879-023-08902-9

Keywords

BMC Infectious Diseases

ISSN: 1471-2334

Contact us