Skip to main content

Drought severity and all-cause mortality rates among adults in the United States: 1968–2014



Little is known about the effect of drought on all-cause mortality, especially in higher income countries such as the United States. As the frequency and severity of droughts are likely to increase, understanding the connections between drought and mortality becomes increasingly important.


Our exposure variable was an annual cumulative drought severity score based on the 1-month, county-level Standardized Precipitation Evapotranspiration Index. The outcome variables of demographic subgroup-specific all-cause mortality count data per year were obtained from the National Vital Statistics System. Any counts below 10 deaths were censored in that demographic group per county. We modeled county-stratum-year mortality using interval-censored negative binomial regression with county-level random intercepts, for each combined age-race-sex stratum either with or without further stratification by climate regions. Fixed effects meta-regression was used to test the associations between age, race, sex, and region with the drought-mortality regression coefficients. Predictive margins were then calculated from the meta-regression model to estimate larger subgroup (e.g., ‘race’ or ‘sex’) associations of drought with mortality.


Most of the results were null for associations between drought severity and mortality, across joint strata of race, age, sex and region, but incidence rate ratios (IRRs) for 17 subgroups were significant after accounting for the multiple testing; ten were < 1 indicating a possible protective effect of drought on mortality for that particular subpopulation. The meta-regression indicated heterogeneity in the association of drought with mortality according to race, climate region, and age, but not by sex. Marginal means of the estimated log-incidence rate ratios differed significantly from zero for age groups 25–34, 35–44, 45–54 and 55–64; for the white race group; and for the South, West and Southwest regions, in the analysis that included wet county-years. The margin of the meta-regression model suggested a slightly negative, but not statistically significant, association of drought with same-year mortality in the overall population.


There were significant, heterogeneous-direction associations in subpopulation-stratified models, after controlling for multiple comparisons, suggesting that the impacts of drought on mortality may not be monolithic across the United States. Meta-regression identified systematic differences in the associations of drought severity with all-cause mortality according to climate region, race, and age. These findings suggest there may be important contextual differences in the effects of drought severity on mortality, motivating further work focused on local mechanisms. We speculate that some of the estimated negative associations of drought severity with same-year mortality could be consistent with either a protective effect of drought on total mortality in the same year, or with a delayed health effect of drought beyond the same year. Further research is needed to clarify associations of drought with more specific causes of death and with sublethal health outcomes, for specific subpopulations, and considering lagged effects occurring beyond the same year as the drought.

Peer Review reports


Drought generally occurs as a result of a water deficit for a given area, but the actual definition of drought differs depending on the measures considered, such as precipitation and temperature [1]. Droughts are commonly divided into four major types based on their environmental and human impacts, including meteorological drought (i.e., abnormally low precipitation), hydrological drought (i.e., precipitation shortages that impact the surface and groundwater levels), agricultural drought (i.e., decreased soil moisture that impacts crops), and socioeconomic drought (i.e., decreased water supply that affects people and supply of goods) [1]. Further, droughts can be characterized by additional factors such as duration, intensity, spatial distribution, frequency, and rate of onset [2].

In 2012, the Intergovernmental Panel on Climate Change (IPCC) reported that medium confidence exists that droughts in some areas of the world have increased in duration and intensity since the 1950s and this pattern will likely continue to over the next century, due to decreased precipitation and/or increased evapotranspiration [3]. A recent IPCC special report found that a shift from 1.5° to 2 °C change in global temperature would cause some areas to experience hotter extreme temperatures, which could increase the probability of dryness and reduced water availability [4]. Within the United States, future climate projections suggest that droughts will be more intense in the Southwest and consecutive dry days are projected to increase over much of the country [5]. At the same time, increases in the heaviest precipitation events are also projected for the U.S., suggesting an increase in both wet and dry extremes could occur [5].

Droughts can lead to widespread and multifaceted regional impacts. Recent droughts in Kenya, the Mediterranean and California have led to crises including food insecurity, political instability, and severe economic damage, respectively [6]. Numerous studies and articles have researched and constructed pathways that connect drought to health effects [7,8,9,10,11,12,13]. There are multiple mechanisms by which negative health effects could result from drought. Droughts can cause decreased freshwater availability and impact water quality [13], including a possible increased risk of microbial contamination and cyanobacteria blooms [8, 13, 14]. Other potential pathways through which drought might affect health are through effects on agricultural production and food insecurity [15, 16], increased occurrence and intensity of heat events [6, 17], effects on air quality [18] including through increased risk of wildfires [6, 19, 20] and dust events (which may impact microbial transport and transmission of coccidioidomycosis (Valley fever)) [6, 21,22,23,24], altered vector-borne disease transmission dynamics (e.g. West Nile virus) [25, 26], and by possibly acting as a catalyst for negative behavioral or mental health outcomes [27]. Some of these mechanisms can result in disease morbidity and/or mortality [28,29,30,31,32,33,34,35,36,37,38,39,40,41,42,43,44], which may result in certain populations being more at-risk.

A recent systematic review of evidence for the health effects of droughts showed that high-quality, quantitative studies on the association of drought and mortality are limited; while some studies have shown high prevalence of mortality at the time of drought, they do not prove causality [7]. A previous study of all-cause mortality and drought in the United States estimated the percentage of drought-associated health risk in adults aged 65 years and older, comparing full drought, non-drought and worsening drought periods [12]. The study found a 1.55% increased risk of all-cause mortality during high-severity worsening drought periods in western counties of the United States and increased mortality during worsening drought compared to non-drought periods in counties where drought occurred less frequently [12]. A few other studies have examined mortality related to drought among populations outside of the United States. A meta-analysis of aggregated mortality surveys of death rates in Ethiopia concluded there was no likely association between drought and death rates of children under the age of five, but did find an increase in mortality with an increase in the prevalence of global acute malnutrition [9]. A study of suicide risk between 1970 and 2007, among 30–49 year old males in rural New South Wales, Australia, found the relative risk increased 15% (95% confidence interval = 8–22%) when the drought index increased from the 1st to the 3rd quartile, after controlling for season, region and long-term trends; however the risk decreased for rural females over 30 years of age with increasing drought index [11]. Also, a recent study in Spain found statistically significant associations between drought periods and daily mortality in Galicia, north-western Spain [45]. Categorizing the death cases into natural, circulatory and respiratory causes, this article concludes that respiratory cause of mortality was the most strongly associated cause of death [45]. These studies suggest there may be an effect of drought on mortality, but a dearth of studies on the association between drought and mortality in higher income countries remains. In addition, drought health effects could differ by subpopulation, and this is another area in need of more research. Because of this, we chose to evaluate the associations between droughts and same-year all-cause mortality in adults in United States from 1968 to 2014 in the current study.



Standardized precipitation evapotranspiration index

Our study’s definition of drought is based on 1-month, county-level Standardized Precipitation Evapotranspiration Index (SPEI) data for the study period of record for the contiguous United States. The SPEI is a climatic drought index based on precipitation and temperature data [46]. Drought indices are quantitative measures derived from the integration of relevant drought indicators (i.e. precipitation) into a single numerical value for drought characterization [2]. Different indices use different variables as indicators, and therefore may reflect different characteristics of drought (e.g. meteorological vs. hydrological). The SPEI is similar to the Standardized Precipitation Index (SPI), but also accounts for the effects of temperature variability on drought assessments through potential evapotranspiration [46]. The SPEI uses meteorological data to determine the summation of deficits and surpluses of available moisture across time, and is based off of historical climatology data for the region [46]. For an SPEI of a given time period, say t months, a time series is constructed by summing the potential evapotranspiration (PET) values from the preceding t-1 months, then fitted to a log-logistic density function to standardize the values (mean 0, SD 1) and make them comparable across time, space, and time scales [47]. For the SPEI used in this analysis, we only used a one-month SPEI. The SPEI ranges from − 3 to 3. Values above 0 indicate relatively wet conditions, and values below 0 indicate relatively dry conditions, with numbers further from zero in either direction representing more extreme conditions [46].

Drought severity score

We defined drought to be when two or more contiguous months in a year have each sustained at least at a “moderately dry” category with SPEI ≤ − 1, with the drought accruing a cumulative SPEI ≤ − 5. Each month is associated with a binary indicator of 0 or 1, depending on being part of a drought event or not. We multiplied the indicator for a month being contained in a drought by the SPEI of that month providing a monthly drought score for each county. The annual drought severity score was then calculated by summing the monthly drought scores within each year, for each county. The code for calculating this score is provided in Additional file 1: Appendix I.

We also created an analogous wetness severity measure. An abnormally wet period was defined as when 2 or more contiguous months have sustained at least at a “moderately wet” category with SPEI ≥1, with the wet period accruing cumulative SPEI > 5. We multiplied the wet period indicator by the SPEI value (positive during wet period) to obtain a monthly wetness score, which was likewise summed over each year to obtain an annual wetness severity score.

All-cause mortality

Mortality data and population counts were extracted from the Mortality Information and Research Analytics (MOIRA) system (formerly known as the Mortality and Population Data System (MPDS) [48, 49], which is a repository and retrieval system for detailed mortality data obtained from the National Center for Health Statistics (NCHS) and the US Census. The database contains death counts and populations within the contiguous United States over 1968–2014, within demographic strata (i.e., age group, race and sex joint categories). The counts for any sub-national cell with fewer than 10 deaths are suppressed per NCHS policy. Our analysis included data from adults ages 25 and older with age categories 25–34, 35–44, 45–54, 55–64, 65–74, 75–84, and 85 years and older.

Study sample

This was a complete case analysis that included a total of 5,099,479 observations from 42 demographic strata in 3004 counties and county-equivalents with complete data on stratum-year mortality counts, stratum-year population, and county-year drought scores. Of the observations, 1,610,899 (31.6%) were censored. Of the 6,121,164 observations with mortality counts, county-year-stratum population size was recorded as 0 for 1,013,237 of the observations. A majority of these were from the black and “other” race groups: 596 (0.06%) were from the white race group, 521,091 (51.43%) were from the black race group, and 491,550 (48.51%) were from the “other” race group. Of the observations with population recorded as 0, 1,001,345 (98.85%) had mortality counts recorded as “null” (0). We regarded observations with population equal to 0 as indicating missing data for the actual population size in that county-demographic stratum-year, although it is possible that in some instances a demographic stratum was not represented within a county within a year. These records were therefore excluded from the complete case analysis.

SPEI data was missing for some months or years of some counties, and therefore the calculated drought severity score used in this study was missing. This was the case for 244 county-years (10,248 stratum-county-years). Nine counties were missing SPEI data for all years (0.17%) and were thus excluded entirely from our study. There were 16,380 stratum-county-years without mortality count data. Missing mortality data usually occurred over a consecutive span of several years. For example, La Paz County, Arizona was missing 26 county-years (1092 stratum-county-years) from 1968 to 1993, but had complete data from the remaining years, 1994–2014. Furthermore, La Paz County, Arizona did not form officially until 1983 after separating from Yuma County [50]. This highlights the fact that some of the missing data may be related to changing counties, or counties that formed after 1968, and associated Federal Information Processing System (FIPS) codes for the given spatial location. FIPS codes are unique numerical identifiers for each county or county-equivalent, but the geographic boundaries of counties can change over time, and new counties can form or merge with other counties [50]. As a result, changes to the size and location of a given county in our dataset could occur across time. For example, Yuma County would have decreased in size with the formation of La Paz County.

Table 1 Climate regions and included states of the United States

Overall, out of 6,137,544 potential stratum-county-years that could have been included from the combined datasets, 1,038,065 observations were missing data on drought conditions, mortality, and/or their population had been recorded as 0, resulting in 16.9% of potential county-years being excluded from analysis. However, some of these missing data may be related to changing county boundaries (e.g., counties that formed after 1968), and so spatially those areas may have still been covered even if they appeared to be incomplete records in our dataset. The climate regions, based on those determined by National Oceanic and Atmospheric Association (NOAA), (Table 1) for the United States were used to explore geographic heterogeneity in the association between droughts and mortality [51]. Washington, DC was included with the Northeast climate region because of its geographical location. 

Statistical analysis

We modeled county-stratum-year mortality using an interval-censored (ranging from 1 to 9) negative binomial regression model with random intercepts. This model is similar to censored mixed negative binomial regression models used in other studies, except it accounts for interval censoring rather than left or right censoring [52, 53]. One other similar-design epidemiological study used a Bayesian censored mixed-effects Poisson regression model to estimate associations with censored autism counts [54].

The interval censoring accounts for the censored death counts, while the random intercepts allow different county-level baseline mortality rates and account for correlations due to repeated observations of the same county. Our negative binomial regression model assumed a linear predictor for the log of the mean number of deaths of the following form:

$$ \ln \left({\upmu}_{\mathrm{ij}}\right)=\left({\upbeta}_0+{\mathrm{b}}_{0\mathrm{i}}\right)+{\upbeta}_1{\mathrm{X}}_{\mathrm{ij}}+{\upbeta}_2\mathrm{Z}+{\mathrm{f}}_{\mathrm{ij}.} $$

where Xij is the drought score for county i in year j, Z is the potential confounder (year, centered at 1991), and (b0i) the random intercept. The offset, fij, is the ln (population) for each county-stratum-year. When the number of deaths (yij) is 0 or > 10 and not censored, the likelihood contribution for county i in year j, conditional on the random intercept, is given by

$$ \Pr \left({\mathrm{Y}}_{\mathrm{ij}}={\mathrm{y}}_{\mathrm{ij}}|{\mathrm{X}}_{\mathrm{ij}}={\mathrm{x}}_{\mathrm{ij}},\mathrm{Z}=\mathrm{z}\right)=\frac{\Gamma \left({\upalpha}^{-1}+{\mathrm{y}}_{\mathrm{ij}}\right)}{\Gamma \left({\upalpha}^{-1}\right){\mathrm{y}}_{\mathrm{ij}}!}{\left(\frac{\upalpha {\upmu}_{\mathrm{ij}}}{1+\upalpha {\upmu}_{\mathrm{ij}}}\right)}^{{\mathrm{y}}_{\mathrm{ij}}}{\left(\frac{1}{1+\upalpha {\upmu}_{\mathrm{ij}}}\right)}^{\frac{1}{\upalpha}} $$

where α is the negative binomial dispersion parameter. When the number of deaths is censored over the interval [1, 9], the conditional likelihood contribution, based on the negative binomial model, becomes:

$$ \Pr \left(1\le {\mathrm{Y}}_{\mathrm{ij}}\le 9|{\mathrm{X}}_{\mathrm{ij}}={\mathrm{x}}_{\mathrm{ij}},\mathrm{Z}=\mathrm{z}\right) $$

We implemented maximum likelihood (ML) estimation for this model using the general likelihood facility available in the SAS NLMIXED procedure using SAS 9.3 software (SAS Institute Inc., Cary, NC), which permits user specification of the log-likelihood conditional on the random effects. We specified the contributions for censored observations by taking advantage of recursive properties of the gamma function.

We ran models separately for each of the 42 age-race-sex joint strata. In a second analysis, we further stratified by the nine NOAA climate regions, for a total of 378 strata, to account for potential geographic heterogeneity in the association of drought score with mortality. We repeated the analysis after excluding abnormally wet county-years (wetness severity score > 0) to obtain a comparison of “drought” vs. “normal” years, rather than “drought” vs. “non-drought” years. Forest plots of the results were created using SAS University Edition/SAS 9.4 M6.

The focus for each model was the estimate of β1, i.e., the natural log of the incidence rate ratio (IRR) corresponding to a 1-unit increase in the yearly drought score. Using SAS University Edition/SAS 9.4 M6 (SAS Institute Inc., Cary, NC), we accounted for multiple testing of the ln(IRR)s against the null of 0 using the SAS PROC MULTTEST procedure using the FDR option, which estimates false discovery rate (FDR)-adjusted p-values. We adjusted all the p-values from the four sets of analyses: with and without stratification by climate regions, and with and without inclusion of abnormally wet years. We used the estimates from the models that successfully converged. Only models that converged were included in the multiple testing correction as specified in the Results section below.

For both the analyses with and without abnormally wet county-years, we calculated estimated mean deaths for the minimum (0) and maximum (15.01) drought score and the difference between the two, to estimate attributable deaths, for the West U.S. region. We then calculated the total estimated deaths by summing the stratum-specific estimated deaths. We chose the West region because the other regions had models that did not successfully converge. For this region, the drought score was 0 for 96.7% of the observations in the main analysis. We used SAS University Edition/SAS 9.4 M6 (SAS Institute Inc., Cary, NC) for this calculation.

To calculate estimated mean deaths, we exponentiated the right side of the negative binomial equation ((β0 + b0i) + β1 Xij + β2 Z + fij)), setting b0i = 0 because the mean and median are 0 by assumption; Xij = 0 (no drought) or the maximum of the drought score for years 1968–2014; Z = 0 (in order to use the mean year, 1991, because we centered the years at 1991); and fij = ln (mean population for each demographic subgroup across 1968–2014, rounded to nearest integer, among observations without missing values). In calculating the mean population and minimum and maximum drought score, we excluded observations where the population was 0, where the death rate was missing, or where the drought score was missing, because these were not included in the analyses. We rounded estimated deaths to the nearest integer.

We applied fixed effects meta-regression to the age-race-sex-region stratified estimated ln(IRR)s (from analyses with wet years included, or without wet years included) using vwls in Stata/SE 14.2. The objective was to assess the association between the regression coefficients from the main analysis (i.e., differences in ln-rates of mortality given drought, within each stratum, adjusted for year and conditional on fixed effects) and potential modifiers of the drought-to-mortality association: age, race, sex, and region. Significance of the association (alpha = 0.05) of each predictor, adjusted for the others, with the drought-mortality regression coefficient was assessed by a Wald test [55].

We then extracted the marginal effects using margins in Stata/SE 14.2 for each subgroup and overall. The marginal effects are the adjusted mean value of the estimated ln (IRR) for each demographic and regional subgroup, and overall across all subgroups.

We also completed a logistic regression of missing population with drought severity score as a predictor to check for possible selection bias using SAS University Edition/SAS 9.4 M6 (SAS Institute Inc., Cary, NC).


Most models in the analyses successfully converged, but there were some estimation problems, where models failed to converge. This issue occurred among the “other” race female subgroup for age categories (years) 25–34, 35–44 and 45–54 in age-race-sex stratified analysis. When stratified by a fourth variable, climate region, 36 of the 378 models failed to converge reliably in the analysis that included abnormally wet county-years. We excluded two additional models that converged, but had estimates suggesting additional convergence issues, i.e. with incident rate ratios (IRRs) of 0.119 (95% confidence interval (CI) = 5.001E− 8, 283,832.31) and 0.14062 (95% CI = 25609E-14, 772,105,431,629). For the analysis stratified by climate regions which excluded abnormally wet county-years, 31 of 378 models failed to converge. We excluded one additional model that resulted in unrealistic IRR of 0.119 (95% CI = 4.797E-8, 297,567.77). The majority of the models that failed to converge were in the “other” race strata.

Overall, most effect estimates were non-significant across categories of race, age, sex, and region, with IRRs for all-cause mortality close to 1. IRRs less than 1 suggest decreasing mortality rates with increasing drought severity, and IRRs greater than 1 suggest increasing mortality rates with increasing drought severity, for a given stratum. When further stratified by climate regions, the majority of results were also null. Exclusion of wet county-years resulted in little to no change in the effect estimates.

A small number of the IRRs were significant after accounting for the multiple testing. These results are presented in Table 2 for the analysis that included abnormally wet county-years, and Additional file 3: Appendix III, Table A1 for the sensitivity analysis that excluded abnormally wet county-years. Among the age-race-sex stratified results, from both the analyses with and without inclusion of abnormally wet county-years, there were significant IRRs, all less than 1, for 5 joint strata: white males 25–34, 35–44, 75–84, and 85+ years of age, and white females 25–34 years of age.

Table 2 Incidence Rate Ratio (IRR) of all-cause mortality per increasing drought severity by demographic subgroup, with 95% confidence intervals (LCL, UCL), raw p-values and false discovery rate-adjusted p-values, for IRRs with adjusted p values < 0.05. Abnormally wet years included in analysis

When further stratified by climate regions, there were 17 significant IRRs from the analysis that included wet county-years (Table 3), and 18 from the analysis that excluded wet county-years (Additional file 3: Appendix III, Table A2). Most of these significant estimates were seen in the white, male subgroups, across a range of ages and regions.

Table 3 Incidence Rate Ratio (IRR) of all-cause mortality per increasing drought severity by demographic and climate region subgroup, with 95% confidence intervals (LCL, UCL), raw p-values and false discovery rate adjusted p-values, for IRRs with adjusted p values < 0.05. Abnormally wet years included in analysis

For both analyses, ten of the significant IRRs were less than 1. Seven significant IRRs were greater than 1 from the analysis that included abnormally wet county-years, and 8 were greater than 1 from the analysis that excluded abnormally wet county-years. The additional significant IRR from the analysis without abnormally wet county-years was for the 65 to 75-year-old “other” race males in the West North Central climate region, with an IRR of 1.066 (95% CI = 1.033, 1.099). This model did not converge in the analysis that included the abnormally wet county-years.

The age-race-sex-region joint strata with IRRs greater than 1 included: 65–74 and 75–84 year-old white females in the Southwest region; 75–84 year-old black males and females, and 85 and older black females in the Northeast region; 45–54 year-old white males in the Central region; and 65–74 year-old black males in the East North Central region. The results for the black and “other” race subgroups usually had wider confidence intervals due to smaller at-risk populations than the white subgroup in some climate regions. Results are also displayed in forest plots in Additional file 4: Appendix IV.

In the attributable deaths analysis, for the West climate region, in the analysis with wet county-years included, we estimated that 8837 fewer deaths would have occurred at a drought score of 15.01 (most severe) compared to a drought score of 0 (no drought) in 1991 (Table 4). Excluding the wet county-years, the total estimated attributable deaths for the West climate region was − 8638 (Additional file 3: Table A4). The attributable deaths ranged from − 2092 for 75–84 year old white females to 121 for 65–74 year old black males in the analysis including wet county-years, and from − 2061 for 75–84 year old white females to 124 for 75–84 year old black males for the analysis without wet county-years. Most of the estimated increases in deaths were among minority subpopulations and decreases in deaths were among with white subpopulations.

Table 4 Expected same-year deaths in Nevada and California in 1991 for counterfactual scenarios of no drought or most severe drought (from years 1968–2014), and attributable deaths for each age-race-sex stratum. Abnormally wet county-years included in analysis

The tests for the significance of associations of the categories of race, age, and climate region, with the drought-mortality regression coefficient from the meta-regression, were significant at alpha = 0.05 in the analyses with and without inclusion of wet county-years (Table 5).

Table 5 Meta-regression p values by covariate, with and without abnormally wet years included in analysis

For the marginal effect estimates, besides a couple of differences, the results were mostly similar between the analyses with and without inclusion of wet county-years (Table 6 with wet county-years; Additional file 3: Appendix III, Table A3 without wet county-years)

Table 6 Incident Rate Ratios (IRRs) by demographic and NOAA climate region subgroup and overall across all groups of margins from fixed effects meta-regression controlling for demographic and climate region variables with 95% confidence intervals (LCL, UCL) and p-values. Abnormally wet years included in analysis

For the marginal effects, by region, IRRs were significant for the Southwest, South, and West regions for both the analyses with and without wet county-years. The IRR for the Southwest region was greater than 1, while the IRRs for the South and West regions were less than 1. When stratified by race, the IRR for the white subgroup was significant, with a pooled IRR 0.999 (95% CI = 0.998–0.999) for both analyses. The age group stratified meta-regression margins for the analyses that included wet county-years resulted in significant IRRs for 4 age groups: 25–34, 35–44, 45–54, and 55–64 year-olds. The IRRs for the 25–34 and 35–44 year-olds were less than 1, while the IRRs for 45–54 and 55–64 year-olds were greater than 1. Differing from the analysis with wet county-years, results for the 85 year and older age group were significant (IRR > 1), and results for the 25–34 year-old age group were not significant, when wet county-years were excluded. The overall marginal effect estimate without stratification, after adjusting for drought-mortality heterogeneity by demographic and climate variables, was not significant.


The majority of demographic and climate region joint strata showed no significant effect of drought severity on mortality rates within the same year in the contiguous United States from 1968 to 2014. Nonetheless, after adjusting for the false discovery rate, we found significant associations remained for some of the stratified analyses, which differed by subpopulation, suggesting possible health effects for certain groups. Based on the estimated attributable deaths, the net effect of drought severity on mortality was a reduction in deaths in the same year. However as shown in Table 4, a majority of increases in deaths occurred with minority populations.

The tests for significance of the fixed effects meta-regression suggest a linear relationship may exist between the regression coefficients from the main analyses and race, age, and climate region, but not sex, after adjusting for year and other demographic and region variables. The marginal effect estimates showed significant results for specific subgroups of these covariates. In the analysis that included wet county-years, the results of the marginal effect estimates suggest that, after adjustment, average mortality rates decreased slightly with increasing drought severity for 25–34 and 35–44 year-olds, while mortality rates increased slightly with increasing drought severity for 45–54 and 55–64 year-olds. In the analysis with wet county-years, the results did not reflect an increase in mortality rate with increasing drought severity in any of the age groups 65 years or older, but when wet county-years were excluded, the results suggested an increase in mortality with increasing drought severity for the 85 year and older age group. The significant IRRs less than 1 for the 25–34 and 35–44 year-olds could suggest that drought is protective against all-cause mortality in these age strata. If there exist both protective and harmful effects of drought on mortality rates, the protective effects, such as potentially decreased flooding, may drive the IRR below 1 for age groups more likely to be exposed to such events, even though harmful effects might still occur. Older populations may be especially vulnerable to a broader range of the potential effects of drought, such as changes in air quality and heat waves [30, 56], and exposures might differ by age group based on factors such as level of outdoor activity or occupational status, which could explain differential effects of drought on mortality rates by age.

A possible protective effect of drought severity on mortality rates was also seen with the significant marginal effects estimate for the white race group. Although no significant marginal effect estimates were found for the black and “other” race groups, smaller sample sizes and non-converging models could result in bias. If these groups are more sensitive to the effects of drought, the results could be conservative and towards the null. Heterogeneity by race subgroup could exist because race is related to socioeconomic status, which is in turn associated with health disparities [57]. Even after controlling for socioeconomic status, disparities exist which may be attributable to factors such as psychosocial stressors, for example [57].

Finally, for the climate regions, the significant marginal effects estimate for the Southwest region (IRR > 1) suggests average increasing mortality rates with increasing drought severity, while the significant marginal effect estimates for the South and West regions (IRRs < 1) suggest average decreasing mortality rates with increasing drought severity. Particulate matter has been shown to increase in the Southwest during droughts, with one study projecting a 20–130% increase in premature attributable deaths due to particulate matter in this region under a number of climate change projections [58]. This may be a potential mechanism for increased mortality rates with increased drought severity, although this mechanism could exist in other climate regions as well [18]. The results of the West and South regions were unexpected, e.g., because harmful mechanisms such as wildfires in the West are expected to exist, and further research would be needed to see if the results could be replicated.

Overall, based on these results, to the extent to which heterogeneity within different combined demographic groups is negligible, there seems to be a region, age and race effect. Since the marginal effect estimates represent an average of the effects across the joint age-race-sex-region strata from the analyses, certain results could be driving the significance of the groups in the meta-regression, and heterogeneity may still exist. This possible heterogeneity is supported by the significant IRRs found for certain subpopulations.

Regarding age, heterogeneity of effects is suggested by the significant results for some of the age-race-sex-region stratified estimates. Within the 25–34 and 35–44 year-old age groups, results suggested decreasing mortality with increasing drought severity for white males in the South, however there were no other significant race-sex-region subgroups. Among 65–74 year-olds, the IRR was greater than 1 for white females in the Southwest, but the remaining significant stratum-specific IRRs for this age group were less than 1. This heterogeneity within age groups is further supported by IRRs both greater than 1 and less than 1 for strata of age groups 75–84 and 85 years and older.

Among race, heterogeneity was also seen. The white race subgroup had the greatest number of significant IRRs after adjusting for multiple testing. The majority of these significant IRRs suggested decreasing mortality rates with increasing drought severity, across a range of age, sex, and climate region subgroups. One possible reason for the greatest proportion of significant findings from the white race subgroup is that this subgroup frequently had the largest sample size, and more precise estimates. The fact that all significant IRRs for the black race group were greater than 1 could indicate this group is more sensitive to the effects of drought severity on mortality rate, for a given age, sex and region combination. For example, among 75–84 year-old black males and females, and 85 and older black females, in the Northeast, all IRRs were greater than 1. Additionally, the only significant result for “other” race group was also greater than 1 and resulted in the largest effect estimate of all the strata-specific results (IRR:1.066 (95% CI = 1.033, 1.099)); this was for 65–74 year-old males in the West North Central region, in the analysis that excluded abnormally wet years. Although the meta-regression and margins did not indicate sex-specific differences in the average effects of drought severity on mortality rates, there are again certain significant age-race-region subgroups, further supporting heterogeneity by combined strata.

Finally, by region, the results align with the marginal effect estimates in that all the significant results for the Southwest region had IRRs greater than 1 and the South and West regions had IRRs less than 1. These results again suggest heterogeneity because they were only found for certain joint strata. It is possible that any effects of drought severity on mortality rates seen in the Southwest region could differ by exposure and vulnerability among certain demographic groups. There may also be specific health effects related to drought in the Northeast that leave vulnerable populations susceptible to exposure or harmful effects, as seen by the significant results of certain combined age-race-sex strata in this region. The observed heterogeneity indicates differential effects and suggests the effects of drought are contextual, based on specific characteristics of the drought and the vulnerability of the populations. However, we cannot exclude the possibility of false positive results. Further research could determine if the results can be replicated in other settings.

Differences in exposure and vulnerability to climate-related events, potential key determinants of disaster risk and impacts, often leave marginalized groups disproportionately impacted [59]. These can often be attributed to intersecting social processes such as discrimination based on class, gender, or ethnicity, for example, rather than to a single cause, and the effects of the individual’s experiences [59]. Stratifying the analysis by age, race, sex and region may be an oversimplification of this potential intersectionality, because the effects of the individual’s experiences based on these categories may not be additive [60]. Future analyses may examine differences in risk at a local level to examine population and drought heterogeneity on a finer spatial scale. Additionally, other variables not considered in the analysis, such as socioeconomic status or rural vs. urban residence, may be important modifiers to examine in future studies.

Although we did not detect significant associations of drought severity with same-year mortality for most subgroups, it is possible that drought conditions could increase some cause-specific mortality rates while decreasing others, resulting in an overall null result. Our analysis only considered all-cause mortality, and specific types of mortality, such as respiratory mortality or cardiovascular mortality, may still have increased rates due to drought and merit further investigation. If there are cause-specific mortalities that increase due to drought, it is important to understand the causes for targeted interventions. The net effect of a reduction in deaths in the same year could be consistent with a lagged effect of drought severity on mortality beyond the same-year time interval considered in the study.

A recent study of drought and the risk of all-cause mortality among people ages 65 years and older in the western United States found no significant association during full-drought periods or low-severity worsening drought periods, but did find an increased risk of mortality during high-severity worsening drought periods [12]. They also found increased mortality during worsening drought compared to non-drought periods in counties where drought occurred less frequently [12]. This study differed from ours in that it focused on counties west of the Mississippi river only, and was restricted to ages 65 and older [12]. Additionally, the study compared days in drought versus non-drought periods, using the weekly U.S. Drought Monitor, so it used a different index for classifying drought and a different temporal scale [12]. This study further stratified by drought severity, instead of using a continuous variable, and considered factors such as worsening drought, or frequency of droughts in the county, which were not included in our study [12]. For these reasons, the results may not be directly comparable. Considering worsening droughts or frequency of droughts in our models could be an option for further analysis.

A strength of the study was the creation of a novel drought severity score designed to capture drought intensity and duration over the year with a continuous scale from the available one-month SPEI data. This scale allowed us to look for effects related to drought severity, and not simply classify exposure as drought vs. non-drought in the analysis. Another strength of the study is that we developed a way to deal with interval-censored CDC mortality data through an interval-censored negative binomial regression model with random county effects. We then supplemented the subgroup-specific model results by performing the meta-regression and marginal effects analyses. This provides a principled way to incorporate all of the available information about county/month-specific mortality counts, while accounting for the differences between the combined age-race-sex-region strata.

This study has some limitations, including the ecological study design; future analyses might consider alternative designs. There was censoring of the number of deaths over the interval 1 to 9 in the dataset. We used interval censored negative binomial regression to handle the missing data, but this necessarily reduces estimation precision relative to the unavailable complete data, particularly for strata with small numbers at risk. Additionally, there were missing data on population, mortality, and SPEI, further decreasing sample size, and potentially leading to bias in our estimated drought effects if the missing observations differed in terms of drought exposure or mortality from those included. A logistic regression of missing population frequency with drought score as the predictor was completed for the overall dataset, and for the significant strata (Additional file 3: Table A5). Several of the analyses found the drought score estimate significant in the model. As there was an association between exposure and inclusion of a county-year, there is a formal possibility that there could be selection bias if inclusion of a county-year was also differential by the mortality rate in that county-year.

Small strata resulted in non-convergence and wide confidence intervals for some models. This especially affected the “other” race category, which had the greatest number of non-converging models. If the more frequently excluded demographic groups due to smaller sample size are more vulnerable than the included groups, and if they also differ in terms of inclusion based on drought exposure status, selection bias could result for the overall effect estimate, and we may be unable to detect impacts of drought on mortality among particularly vulnerable groups.

Also, despite the strengths of the SPEI, we only used one index for creating the annual drought severity score, whereas another index might have resulted in different classification of drought periods and severity. Additionally, the analysis did not specifically account for differences in the types of droughts, such as agricultural vs. meteorological, which could potentially affect health outcomes differently since they focus on different issues (e.g. effects on crops for agricultural drought) [1]. Further, we did not account for the effect of previous years of drought, except when classifying drought vs. non-drought months for the derivation of the drought severity score. Differences in duration or frequency could theoretically result in different exposure durations to potentially harmful or preventative drought-related conditions or impact a community’s adaptive capacity to droughts.

This analysis did not account for potential lagged effects of drought, or interactions between abnormally wet periods following drought, or drought following abnormally wet periods. Differences could hypothetically exist from mechanisms resulting from increased rainfall after a drought such as landslides, or indirect factors such as changes in water quality [61] and disease transmission (e.g. increased incidence of Valley fever after a period of rain following drought) [7, 62]. Drought and wildfires (which might occur more frequently during droughts) [20] could lead to landslides, in part due to removal of protective vegetation [63], and landslides may be triggered by periods of intense rainfall [63]. Interaction mechanisms could also occur from drought following wet periods. Increased plant growth from a wet period followed by drought might lead to build-up of fuel for wildfires in the form of vegetation, which is related to fire risk and spread [64]. The mechanisms that could impact health related to dry or wet conditions could also differ by geographic location, as suggested by a study that found differences by location (i.e. California and Arizona) in direction of the association between coccidioidomycosis incidence and soil moisture [65]. These types of interactions may be useful avenues for further investigation. Future investigations might also consider a larger suite of potential confounders.

While we also did not directly include a variable for temperature in the analysis, SPEI accounts for temperature when calculating evapotranspiration [46]. Temperature has a complex relationship with drought. Since heat can impact heat-related mortality [37] and influence drought severity through evapotranspiration [47, 66, 67] it could be considered a confounder. Multiple variables are considered in the calculation of evapotranspiration [46] and two droughts of the same SPEI could have different extreme or average high temperatures, and therefore differences in heat-related deaths; however, because heat is a component of drought [46] it may not be logical to separate the two. Heat could also be a mediator on the causal pathway between drought and mortality, since drought conditions can increase heat events [17], which can in turn impact heat-related mortality. Interactions could also possibly occur between the effects of heat and drought on all-cause mortality.


The lack of significant association between drought severity and all-cause mortality overall, and for many, but not all, subgroups after adjusting for multiple testing, could be, in part, due to the heterogeneity of effects based on categories of race, age, climate region, or other unmeasured variables. Other potential explanations include: contextual heterogeneity in the effects of drought, an observed null effect of subgroups due to cancellation of the protective and harmful effects, effects of drought on cause-specific but not all-cause mortality, limitations in the methods of our analysis, lagged effects of drought on mortality, or some combination of these factors. Our results suggest that the potential associations between drought and all-cause mortality could differ by race, age, and climate region and might be more nuanced than the situations considered here. Therefore, the associations may merit additional research based on more detailed data.

Availability of data and materials

The datasets generated during and/or analyzed during the current study are available in the CDC Vital Statistics Online Data Portal, The drought/climate datasets during and/or analyzed during the current study are available from the authors on reasonable request.


  1. Wilhite DA, Glantz MH. Understanding the drought phenomenon: the role of definitions. Water Int. 1985;10(3):111–20..

    Article  Google Scholar 

  2. Zargar A, Sadiq R, Naser B, Khan FI. A review of drought indices. Environ Rev. 2011;19(NA):333–49..

    Article  Google Scholar 

  3. Intergovernmental Panel on Climate Change. Summary for Policymakers. A Special Report of Working Groups I and II of the Intergovernmental Panel on Climate Change. In: Field CB, Barros V, Stocker TF, Qin D, Dokken DJ, Ebi KL, Mastrandrea MD, Mach KJ, Plattner G-K, Allen SK, Tignor M, Midgley PM, editors. Cambridge: Cambridge University Press; 2012. p. 3–21.

  4. Hoegh-Guldberg O, D. Jacob, M. Taylor, M. Bindi, S. Brown, I. Camilloni, A. Diedhiou, R. Djalante, K.L. Ebi, F. Engelbrecht, J. Guiot, Y. Hijioka, S. Mehrotra, A. Payne, S.I. Seneviratne, A. Thomas, R. Warren, and G. Zhou. IPCC: Impacts of 1.5°C Global Warming on Natural and Human Systems. In: Global Warming of 15°C An IPCC Special Report on the impacts of global warming of 15°C above pre-industrial levels and related global greenhouse gas emission pathways, in the context of strengthening the global response to the threat of climate change, sustainable development, and efforts to eradicate poverty Edited by Masson-Delmotte V, P. Zhai, H.-O. Pörtner, D. Roberts, J. Skea, P.R. Shukla, A. Pirani, W. Moufouma-Okia, C. Péan, R. Pidcock, S. Connors, J.B.R. Matthews, Y. Chen, X. Zhou, M.I. Gomis, E. Lonnoy, T. Maycock, M. Tignor, and T. Waterfield; 2018.

  5. Walsh J, Wuebbles D, Hayhoe K, Kossin J, Kunkel K, Stephens G, Thorne P, Vose R, Wehner M, Willis J, Anderson D, Doney S, Feely R, Hennon P, Kharin V, Knutson T, Landerer F, Lenton T, Kennedy J, Somerville R. Ch. 2: Our Changing Climate. In: JMM, Richmond T(TC), Yohe GW, editors. Climate Change Impacts in the United States: The Third National Climate Assessment, U.S. Global Change Research Program, United States. 2014; p. 19–67.

  6. Bell JE, Brown CL, Conlon K, Herring S, Kunkel KE, Lawrimore J, Luber G, Schreck C, Smith A, Uejio C. Changes in extreme events and the potential impacts on human health. J Air Waste Manage Assoc. 2017;68(4):265–87.

  7. Stanke C, Kerac M, Prudhomme C, Medlock J, Murray V. Health effects of drought: a systematic review of the evidence. PLoS Curr. 2013;5(1).

  8. Walter JM, Lopes FAC, Lopes-Ferreira M, Vidal LM, Leomil L, Melo F, de Azevedo GS, Oliveira RMS, Medeiros AJ, Melo ASO, et al. Occurrence of harmful cyanobacteria in drinking water from a severely drought-impacted semi-arid region. Front Microbiol. 2018;9:176.

    Article  Google Scholar 

  9. Delbiso TD, Altare C, Rodriguez-Llanes JM, Doocy S, Guha-Sapir D. Drought and child mortality: a meta-analysis of small-scale surveys from Ethiopia. Sci Rep. 2017;7(1):2212.

    Article  CAS  Google Scholar 

  10. Yusa A, Berry P, Cheng JJ, Ogden N, Bonsal B, Stewart R, Waldick R. Climate change, drought and human health in Canada. Int J Environ Res Public Health. 2015;12(7):8359–412.

    Article  CAS  Google Scholar 

  11. Hanigan IC, Butler CD, Kokic PN, Hutchinson MF. Suicide and drought in New South Wales, Australia, 1970-2007. Proc Natl Acad Sci U S A. 2012;109(35):13950–5.

    Article  CAS  Google Scholar 

  12. Berman JD, Ebisu K, Peng RD, Dominici F, Bell ML. Drought and the risk of hospital admissions and mortality in older adults in western USA from 2000 to 2013: a retrospective study. Lancet Planetary Health. 2017;1(1):e17–25.

    Article  Google Scholar 

  13. Bell JE, Herring SC, Jantarasami L, Adrianopoli C, Benedict K, Conlon K, Escobar V, Hess J, Luvall J, Garcia-Pando CP, Quattrochi D, Runkle J, Schreck CJ III. Ch. 4: Impacts of Extreme Events on Human Health. In: The Impacts of Climate Change on Human Health in the United States: A Scientific Assessment. Washington, DC: U.S. Global Change Research Program; 2016. p. 99–128.

    Chapter  Google Scholar 

  14. Trtanj J, Jantarasami L, Brunkard J, Collier T, Jacobs J, Lipp E, McLellan S, Moore S, Paerl H, Ravenscroft J, Sengco M, Thurston J. Ch. 6: Climate Impacts on Water-Related Illness. In: The Impacts of Climate Change on Human Health in the United States: A Scientific Assessment. Washington, DC: U.S. Global Change Research Program; 2016. p. 157–88.

    Chapter  Google Scholar 

  15. Leng G, Hall J. Crop yield sensitivity of global major agricultural countries to droughts and the projected changes in the future. Sci Total Environ. 2019;654:811–21.

    Article  CAS  Google Scholar 

  16. Porter JR, Xie L, Challinor AJ, Cochrane K, Howden SM, Iqbal MM, Lobell DB, Travasso MI. Food security and food production systems. In: Field CB, Barros VR, Dokken DJ, Mach KJ, Mastrandrea MD, Bilir TE, Chatterjee M, Ebi KL, Estrada YO, Genova RC, Girma B, Kissel ES, Levy AN, MacCracken S, Mastrandrea PR, White LL, editors. Climate Change 2014: Impacts, Adaptation, and Vulnerability Part A: Global and Sectoral Aspects Contribution of Working Group II to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. Cambridge: Cambridge University Press; 2014. p. 485–533.

  17. Mueller B, Seneviratne SI. Hot days induced by precipitation deficits at the global scale. Proc Natl Acad Sci. 2012;109(31):12398.

    Article  CAS  Google Scholar 

  18. Wang Y, Xie Y, Dong W, Ming Y, Wang J, Shen L. Adverse effects of increasing drought on air quality via natural processes. Atmos Chem Phys. 2017;17(20):12827–43.

    Article  CAS  Google Scholar 

  19. Mallia DV, Lin JC, Urbanski S, Ehleringer J, Nehrkorn T. Impacts of upwind wildfire emissions on CO, CO2, and PM2.5 concentrations in Salt Lake City, Utah. J Geophysical Res. 2015;120(1):147–66.

    CAS  Google Scholar 

  20. Vose JM, Clark JS, Luce CH, Patel-Weynand T, eds. Effects of drought on forests and rangelands in the United States: A comprehensive science synthesis. Gen. Tech. Rep. WO-93b. Washington: U.S. Department of Agriculture, Forest Service, Washington Office; 2016. p. 1-289.

  21. Rosselli R, Fiamma M, Deligios M, Pintus G, Pellizzaro G, Canu A, Duce P, Squartini A, Muresu R, Cappuccinelli P. Microbial immigration across the Mediterranean via airborne dust. Sci Rep. 2015;5:16306.

    Article  CAS  Google Scholar 

  22. Centers for Disease Control and Prevention. Increase in coccidioidomycosis--Arizona, 1998-2001. MMWR Morb Mortal Wkly Rep. 2003;52(6):109–12.

  23. Pu B, Ginoux P. Projection of American dustiness in the late 21(st) century due to climate change. Sci Rep. 2017;7(1):5553.

    Article  CAS  Google Scholar 

  24. Tong DQ, Wang JXL, Gill TE, Lei H, Wang B. Intensified dust storm activity and valley fever infection in the southwestern United States. Geophys Res Lett. 2017;44(9):4304–12.

    Article  Google Scholar 

  25. Paull SH, Horton DE, Ashfaq M, Rastogi D, Kramer LD, Diffenbaugh NS, Kilpatrick AM. Drought and immunity determine the intensity of West Nile virus epidemics and climate change impacts. Proc Biol Sci. 2017;284(1848):20162078.

  26. Johnson BJ, Sukhdeo MV. Drought-induced amplification of local and regional West Nile virus infection rates in New Jersey. J Med Entomol. 2013;50(1):195–204.

    Article  CAS  Google Scholar 

  27. Vins H, Bell J, Saha S, Hess JJ. The mental health outcomes of drought: a systematic review and causal process diagram. Int J Environ Res Public Health. 2015;12(10):13251–75.

    Article  Google Scholar 

  28. Dockery DW, Pope CA, Xu X, Spengler JD, Ware JH, Fay ME, Ferris BG, Speizer FE. An association between air pollution and mortality in six U.S. cities. N Engl J Med. 1993;329(24):1753–9.

    Article  CAS  Google Scholar 

  29. Qiu H, Pun VC, Tian L. Short-term effects of fine and coarse particles on deaths in Hong Kong elderly population: An analysis of mortality displacement. Environ Poll (Barking, Essex : 1987). 2018;241:148–54.

    Article  CAS  Google Scholar 

  30. Zhang L, Zhang Z, Ye T, Zhou M, Wang C, Yin P, Hou B. Mortality effects of heat waves vary by age and area: a multi-area study in China. Environ Health. 2018;17(1):54.

    Article  Google Scholar 

  31. Zhang K, Chen TH, Begley CE. Impact of the 2011 heat wave on mortality and emergency department visits in Houston, Texas. Environ Health. 2015;14:11.

    Article  Google Scholar 

  32. Faustini A, Alessandrini ER, Pey J, Perez N, Samoli E, Querol X, Cadum E, Perrino C, Ostro B, Ranzi A, et al. Short-term effects of particulate matter on mortality during forest fires in southern Europe: results of the MED-PARTICLES project. Occup Environ Med. 2015;72(5):323–9.

    Article  Google Scholar 

  33. Fann N, Alman B, Broome RA, Morgan GG, Johnston FH, Pouliot G, Rappold AG. The health impacts and economic value of wildland fire episodes in the U.S.: 2008–2012. Sci Total Environ. 2018;610-611:802–9.

    Article  CAS  Google Scholar 

  34. Crooks JL, Cascio WE, Percy MS, Reyes J, Neas LM, Hilborn ED. The association between dust storms and daily non-accidental mortality in the United States, 1993-2005. Environ Health Perspect. 2016;124(11):1735–43.

    Article  Google Scholar 

  35. Carmichael WW, Azevedo SM, An JS, Molica RJ, Jochimsen EM, Lau S, Rinehart KL, Shaw GR, Eaglesham GK. Human fatalities from cyanobacteria: chemical and biological evidence for cyanotoxins. Environ Health Perspect. 2001;109(7):663–8.

    Article  CAS  Google Scholar 

  36. Petersen LR, Brault AC, Nasci RS. West Nile virus: review of the literature. Jama. 2013;310(3):308–15.

    Article  CAS  Google Scholar 

  37. Cheng J, Xu Z, Bambrick H, Su H, Tong S, Hu W. Heatwave and elderly mortality: An evaluation of death burden and health costs considering short-term mortality displacement. Environ Int. 2018;115:334–42.

    Article  Google Scholar 

  38. Nguyen C, Barker BM, Hoover S, Nix DE, Ampel NM, Frelinger JA, Orbach MJ, Galgiani JN. Recent advances in our understanding of the environmental, epidemiological, immunological, and clinical dimensions of coccidioidomycosis. Clin Microbiol Rev. 2013;26(3):505–25.

    Article  CAS  Google Scholar 

  39. Youssouf H, Liousse C, Roblou L, Assamoi EM, Salonen RO, Maesano C, Banerjee S, Annesi-Maesano I. Non-accidental health impacts of wildfire smoke. Int J Environ Res Public Health. 2014;11(11):11772–804.

    Article  Google Scholar 

  40. Horne BD, Joy EA, Hofmann MG, Gesteland PH, Cannon JB, Lefler JS, Blagev DP, Korgenski EK, Torosyan N, Hansen GI, et al. Short-term elevation of fine particulate matter air pollution and acute lower respiratory infection. Am J Respir Crit Care Med. 2018;198(6):759–66.

    Article  Google Scholar 

  41. Analitis A, Katsouyanni K, Dimakopoulou K, Samoli E, Nikoloulopoulos AK, Petasakis Y, Touloumi G, Schwartz J, Anderson HR, Cambra K, et al. Short-term effects of ambient particles on cardiovascular and respiratory mortality. Epidemiol (Cambridge, Mass). 2006;17(2):230–3.

    Article  Google Scholar 

  42. Reid CE, Brauer M, Johnston FH, Jerrett M, Balmes JR, Elliott CT. Critical review of health impacts of wildfire smoke exposure. Environ Health Perspect. 2016;124(9):1334–43.

    Article  Google Scholar 

  43. Morman SA, Plumlee GS. The role of airborne mineral dusts in human disease. Aeolian Res. 2013;9:203–12.

    Article  Google Scholar 

  44. Dimsdale JE. Psychological stress and cardiovascular disease. J Am Coll Cardiol. 2008;51(13):1237–46.

    Article  Google Scholar 

  45. Salvador C, Nieto R, Linares C, Diaz J, Gimeno L. Effects on daily mortality of droughts in Galicia (NW Spain) from 1983 to 2013. Sci Total Environ. 2019;662:121–33.

    Article  CAS  Google Scholar 

  46. Vicente-Serrano SM, Beguería S, López-Moreno JI. A multiscalar drought index sensitive to global warming: the standardized precipitation evapotranspiration index. J Clim. 2009;23(7):1696–718.

    Article  Google Scholar 

  47. Vicente-Serrano SM and National Center for Atmospheric Research Staff (Eds). Last modified 18 Jul 2015. "The Climate Data Guide: Standardized Precipitation Evapotranspiration Index (SPEI)." Retrieved from

  48. Marsh GM, Youk AO, Stone RA, Sefcik S, Alcorn C. OCMAP-PLUS: a program for the comprehensive analysis of occupational cohort data. J Occup Environ Med. 1998;40(4):351–62.

    Article  CAS  Google Scholar 

  49. Buchanich JM, Doerfler SM, Lann MF, Marsh GM, Burke DS. Improvement in racial disparities in years of life lost in the USA since 1990. PLoS One. 2018;13(4):e0194308.

    Article  CAS  Google Scholar 

  50. Substantial Changes to Counties and County Equivalent Entities: 1970-Present [].

  51. Karl T, Koss WJ. Regional and national monthly, seasonal, and annual temperature weighted by area, 1895–1983. In: Historical climatology series ; 4–3. Asheville: National Climatic Data Center; 1984.

    Google Scholar 

  52. Quiroz J, Wilson JR, Roychoudhury S. Statistical analysis of data from dilution assays with censored correlated counts. Pharm Stat. 2012;11(1):63–73.

    Article  Google Scholar 

  53. Lee WY, Lee K, Chun J, Choe J, Jabłoński P, Lee S-I. Comparison of a culture-based and a PCR-based methods for estimating bacterial abundance on eggshells, with comments on statistical analyses, vol. 84; 2013.

    Google Scholar 

  54. Bartell SM, Lewandowski TA. Administrative censoring in ecological analyses of autism and a Bayesian solution. J Environ Public Health. 2011;2011:202783.

    Article  Google Scholar 

  55. Harbord RM, Higgins JPT. Metaregresion in Stata. Stata J. 2008;8(4):493–519.

    Article  Google Scholar 

  56. Simoni M, Baldacci S, Maio S, Cerrai S, Sarno G, Viegi G. Adverse effects of outdoor pollution in the elderly. J Thorac Dis. 2015;7(1):34–45.

    Google Scholar 

  57. Williams DR, Priest N, Anderson NB. Understanding associations among race, socioeconomic status, and health: patterns and prospects. Health Psychol. 2016;35(4):407–11.

    Article  Google Scholar 

  58. Achakulwisut P, Mickley LJ, Anenberg SC. Drought-sensitivity of fine dust in the US southwest: implications for air quality and public health under future climate change. Environ Res Lett. 2018;13(5):054025.

    Article  CAS  Google Scholar 

  59. Intergovernmetal Panel on Climate Change. Part A: Global and Sectoral Aspects. Contribution of Working Group II to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. In: Field CB, Barros VR, Dokken DJ, Mach KJ, Mastrandrea MD, Bilir TE, Chatterjee M, Ebi KL, Estrada YO, Genova RC, Girma B, Kissel ES, Levy AN, MacCracken S, Mastrandrea PR, White LL, editors. Climate Change 2014: Impacts, Adaptation, and Vulnerability.  Cambridge: Cambridge University Press; 2014. p. 1132.

  60. Bauer GR. Incorporating intersectionality theory into population health research methodology: challenges and the potential to advance health equity. Soc Sci Med. 2014;110:10–7.

    Article  Google Scholar 

  61. Bi W, Weng B, Yuan Z, Yang Y, Xu T, Yan D, Ma J. Evolution of Drought-flood abrupt alternation and its impacts on surface water quality from 2020 to 2050 in the Luanhe River Basin. Int J Environ Res Public Health. 2019;16(5):691.

  62. Pappagianis D. Marked increase in cases of coccidioidomycosis in California: 1991, 1992, and 1993. Clin Infect Dis. 1994;19(Suppl 1):S14–8.

    Article  Google Scholar 

  63. Highland LM, Bobrowsky P. The landslide handbook—A guide to understanding landslides. Reston: U.S. Geological Survey Circular 1325; 2008:129.

  64. Whitlock C, Shafer SL, Marlon J. The role of climate and vegetation change in shaping past and future fire regimes in the northwestern US and the implications for ecosystem management. For Ecol Manag. 2003;178(1):5–21.

    Article  Google Scholar 

  65. Coopersmith EJ, Bell JE, Benedict K, Shriber J, McCotter O, Cosh MH. Relating coccidioidomycosis (valley fever) incidence to soil moisture conditions. GeoHealth. 2017;1:51–63.

    Article  CAS  Google Scholar 

  66. Cook BI, Miller RL, Seager R. Amplification of the north American "dust bowl" drought through human-induced land degradation. Proc Natl Acad Sci U S A. 2009;106(13):4997–5001.

    Article  CAS  Google Scholar 

  67. Huang S, Leng G, Huang Q, Xie Y, Liu S, Meng E, Li P. The asymmetric impact of global warming on US drought types and distributions in a large ensemble of 97 hydro-climatic simulations. Sci Rep. 2017;7:5891.

    Article  CAS  Google Scholar 

Download references


The authors are also grateful for the support provided by the National Integrated Drought Information System. The authors would also like to thank J. Rennie, W. Van Panhuis, and J. Buchanich for their assistance with data access. 

Conflicts of Interest

The funding sponsor had no role in the design of the study; in the collection, analysis, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results. The authors have no potential conflicts of interest to declare.


This work was supported by NOAA through the Cooperative Institute for Climate and Satellites - North Carolina under Cooperative Agreement NA14NES432003 with the support coming from the National Integrated Drought Information System. The effort of MOG was supported by the National Institute of Environmental Health Sciences-funded HERCULES Exposome Research Center (P30 ES019776). The funders had no role in the design of the study, nor the collection, analysis, and interpretation of data, nor the writing of the manuscript.

Author information

Authors and Affiliations



KML conceived the study idea, conducted the statistical analysis, wrote the first draft of the manuscript, and lead revisions of the manuscript. RHL developed the statistical methodology used in the manuscript and supervised the data analysis and interpretation. LAW provided guidance on the statistical analysis as well as interpretation of the results. AMA contributed to the contextualization and interpretation of the findings and support of interpretation by current literature. JEB and MOG oversaw the research at all stages; JEB contributed substantive expertise and MOG methodological expertise. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Matthew O. Gribble.

Ethics declarations

Ethics approval and consent to participate

Not applicable

Consent for publication

Not applicable

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

Supplementary information

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 The Creative Commons Public Domain Dedication waiver ( 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

Lynch, K.M., Lyles, R.H., Waller, L.A. et al. Drought severity and all-cause mortality rates among adults in the United States: 1968–2014. Environ Health 19, 52 (2020).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: