Rainfall events and daily mortality across 645 global locations: two stage time series analysis
BMJ 2024; 387 doi: https://doi.org/10.1136/bmj-2024-080944 (Published 09 October 2024) Cite this as: BMJ 2024;387:e080944Linked Editorial
Rainfall events and adverse health outcomes
- Cheng He, postdoctoral researcher1,
- Susanne Breitner-Busch, senior scientist12,
- Veronika Huber, senior scientist12,
- Kai Chen, associate professor34,
- Siqi Zhang, associate research scientist34,
- Antonio Gasparrini, professor5,
- Michelle Bell, professor67,
- Haidong Kan, professor8,
- Dominic Royé, senior researcher9,
- Ben Armstrong, professor10,
- Joel Schwartz, professor11,
- Francesco Sera, associate professor12,
- Ana Maria Vicedo-Cabrera, assistant professor1314,
- Yasushi Honda, professor emeritus15,
- Jouni J K Jaakkola, professor1617,
- Niilo Ryti, lecturer16,
- Jan Kyselý, associate professor1819,
- Yuming Guo, professor2021,
- Shilu Tong, professor2223,
- Francesca de’Donato, senior researcher24,
- Paola Michelozzi, senior scientist24,
- Micheline de Sousa Zanotti Stagliorio Coelho, researcher25,
- Paulo Hilario Nascimento Saldiva, professor26,
- Eric Lavigne, research scientist2728,
- Hans Orru, professor29,
- Ene Indermitte, lecturer29,
- Mathilde Pascal, researcher30,
- Patrick Goodman, professor31,
- Ariana Zeka, associate professor32,
- Yoonhee Kim, associate professor33,
- Magali Hurtado Diaz, professor34,
- Eunice Elizabeth Félix Arellano, researcher34,
- Ala Overcenco, senior scientist35,
- Jochem Klompmaker, researcher36,
- Shilpa Rao, senior scientist37,
- Alfonso Diz-Lois Palomares, researcher37,
- Gabriel Carrasco, assistant professor3839,
- Xerxes Seposo, associate professor40,
- Susana das Neves Pereira da Silva, data scientist41,
- Joana Madureira, assistant researcher414243,
- Iulian-Horia Holobaca, professor44,
- Noah Scovronick, assistant professor45,
- Fiorella Acquaotta, researcher46,
- Ho Kim, professor47,
- Whanhee Lee, assistant professor48,
- Masahiro Hashizume, professor49,
- Aurelio Tobias, associate research professor5051,
- Carmen Íñiguez, professor5253,
- Bertil Forsberg, professor54,
- Martina S Ragettli, project leader5556,
- Yue Leon Guo, professor575859,
- Shih-Chun Pan, postdoctoral researcher60,
- Samuel Osorio, doctoral student61,
- Shanshan Li, professor2021,
- Antonella Zanobetti, principal research scientist11,
- Tran Ngoc Dang, lecturer62,
- Do Van Dung, associate professor63,
- Alexandra Schneider, senior scientist1
- 1Institute of Epidemiology, Helmholtz Zentrum München – German Research Center for Environmental Health, Neuherberg, Germany
- 2Institute for Medical Information Processing, Biometry, and Epidemiology, Medical Faculty, Ludwig-Maximilians-Universität München, Munich, Germany
- 3Department of Environmental Health Sciences, Yale School of Public Health, New Haven, CT, USA
- 4Yale Center on Climate Change and Health, Yale School of Public Health, New Haven, CT, USA
- 5Environment and Health Modelling Lab, Department of Public Health, Environments and Society, London School of Hygiene and Tropical Medicine, London, UK
- 6School of the Environment, Yale University, New Haven, CT, USA
- 7School of Health Policy and Management, College of Health Sciences, Korea University, Seoul, Republic of Korea
- 8Department of Environmental Health, School of Public Health, Fudan University, Shanghai, China
- 9Climate Research Foundation, CIBER of Epidemiology and Public Health, Madrid, Spain
- 10Department of Public Health Environments and Society, London School of Hygiene and Tropical Medicine, London, UK
- 11Department of Environmental Health, Harvard T.H. Chan School of Public Health, Boston, MA, USA
- 12Department of Statistics, Computer Science and Applications “G Parenti,” University of Florence, Florence, Italy
- 13Institute of Social and Preventive Medicine, University of Bern, Bern, Switzerland
- 14Oeschger Center for Climate Change Research, University of Bern, Bern, Switzerland
- 15Center for Climate Change Adaptation, National Institute for Environmental Studies, Tsukuba, Japan
- 16Center for Environmental and Respiratory Health Research, University of Oulu, Oulu, Finland
- 17Finnish Meteorological Institute, Helsinki, Finland
- 18Institute of Atmospheric Physics, Czech Academy of Sciences, Prague, Czech Republic
- 19Faculty of Environmental Sciences, Czech University of Life Sciences, Prague, Czech Republic
- 20School of Public Health and Preventive Medicine, Monash University, Melbourne, VIC, Australia
- 21Climate, Air Quality Research Unit, School of Public Health and Preventive Medicine, Monash University, Melbourne, VIC, Australia
- 22National Institute of Environmental Health, Chinese Center for Disease Control and Prevention, Beijing, China
- 23School of Public Health and Social Work, Queensland University of Technology, Brisbane, QLD, Australia
- 24Department of Epidemiology, Lazio Regional Health Service, ASL ROMA 1, Rome, Italy
- 25Department of Pathology, Faculty of Medicine, University of Sao Paulo, Sao Paulo, Brazil
- 26INSPER, Sao Paulo, Brazil
- 27School of Epidemiology and Public Health, Faculty of Medicine, University of Ottawa, Ottawa, ON, Canada
- 28Environmental Health Science and Research Bureau, Health Canada, Ottawa, ON, Canada
- 29Department of Family Medicine and Public Health, University of Tartu, Tartu, Estonia
- 30Department of Environmental and Occupational Health, French National Public Health Agency, Saint Maurice, France
- 31Technological University Dublin, Dublin, Ireland
- 32Institute for Global Health, University College London, UK
- 33Department of Global Environmental Health, Graduate School of Medicine, University of Tokyo, Tokyo, Japan
- 34Department of Environmental Health, National Institute of Public Health, Cuernavaca, Morelos, Mexico
- 35National Agency for Public Health of Ministry of Health, Labour and Social Protection of the Republic of Moldova, Chisinau, Republic of Moldova
- 36National Institute for Public Health and the Environment (RIVM), Centre for Sustainability and Environmental Health, Bilthoven, Netherlands
- 37Norwegian institute of Public Health, Oslo, Norway
- 38Institute of Tropical Medicine, Universidad Peruana Cayetano Heredia, Lima, Peru
- 39Division of Infectious Diseases, Department of Medicine, University of California, San Diego, CA, USA
- 40School of Tropical Medicine and Global Health, Nagasaki University, Nagasaki, Japan
- 41National Institute of Health Dr Ricardo Jorge, Portugal
- 42EPIUnit - Instituto de Saúde Pública, Universidade do Porto, Porto, Portugal
- 43Laboratório para a Investigação Integrativa e Translacional em Saúde Pública (ITR), Porto, Portugal
- 44Faculty of Geography, Babeș-Bolyai University, Cluj-Napoca, Romania
- 45Gangarosa Department of Environmental Health, Rollins School of Public Health, Emory University, Atlanta, GA, USA
- 46Department of Earth Sciences, University of Turin, Turin, Italy
- 47Graduate School of Public Health, Seoul National University, Seoul, Republic of Korea
- 48School of Biomedical Convergence Engineering, College of Information and Biomedical Engineering, Pusan National University, Yangsan, Republic of Korea
- 49Department of Global Health Policy, Graduate School of Medicine, University of Tokyo, Tokyo, Japan
- 50Institute of Environmental Assessment and Water Research, Spanish Council for Scientific Research, Barcelona, Spain
- 51School of Tropical Medicine and Global Health, Nagasaki University, Nagasaki, Japan
- 52Department of Statistics and Computational Research, University of Valencia, Valencia, Spain
- 53Ciberesp, Madrid, Spain
- 54Department of Public Health and Clinical Medicine, Umea University, Umea, Sweden
- 55Swiss Tropical and Public Health Institute, Allschwil, Switzerland
- 56University of Basel, Basel, Switzerland
- 57Environmental and Occupational Medicine, National Taiwan University (NTU) College of Medicine and NTU Hospital, Taipei, Taiwan
- 58National Institute of Environmental Health Science, National Health Research Institutes, Zhunan, Taiwan
- 59Graduate Institute of Environmental and Occupational Health Sciences, NTU College of Public Health, Taipei, Taiwan
- 60National Institute of Environmental Health Science, National Health Research Institutes, Zhunan, Taiwan
- 61Department of Environmental Health, University of Sao Paulo, Sao Paulo, Brazil
- 62Institute of Research and Development, Duy Tan University, Da Nang, Vietnam
- 63Department of Environmental Health, Faculty of Public Health, University of Medicine and Pharmacy at Ho Chi Minh City, Ho Chi Minh City, Vietnam
- Correspondence to: C He cheng.he{at}helmholtz-munich.de (or @cheng_he_1130 on X)
- Accepted 17 September 2024
Abstract
Objective To examine the associations between characteristics of daily rainfall (intensity, duration, and frequency) and all cause, cardiovascular, and respiratory mortality.
Design Two stage time series analysis.
Setting 645 locations across 34 countries or regions.
Population Daily mortality data, comprising a total of 109 954 744 all cause, 31 164 161 cardiovascular, and 11 817 278 respiratory deaths from 1980 to 2020.
Main outcome measure Association between daily mortality and rainfall events with return periods (the expected average time between occurrences of an extreme event of a certain magnitude) of one year, two years, and five years, with a 14 day lag period. A continuous relative intensity index was used to generate intensity-response curves to estimate mortality risks at a global scale.
Results During the study period, a total of 50 913 rainfall events with a one year return period, 8362 events with a two year return period, and 3301 events with a five year return period were identified. A day of extreme rainfall with a five year return period was significantly associated with increased daily all cause, cardiovascular, and respiratory mortality, with cumulative relative risks across 0-14 lag days of 1.08 (95% confidence interval 1.05 to 1.11), 1.05 (1.02 to 1.08), and 1.29 (1.19 to 1.39), respectively. Rainfall events with a two year return period were associated with respiratory mortality only, whereas no significant associations were found for events with a one year return period. Non-linear analysis revealed protective effects (relative risk <1) with moderate-heavy rainfall events, shifting to adverse effects (relative risk >1) with extreme intensities. Additionally, mortality risks from extreme rainfall events appeared to be modified by climate type, baseline variability in rainfall, and vegetation coverage, whereas the moderating effects of population density and income level were not significant. Locations with lower variability of baseline rainfall or scarce vegetation coverage showed higher risks.
Conclusion Daily rainfall intensity is associated with varying health effects, with extreme events linked to an increasing relative risk for all cause, cardiovascular, and respiratory mortality. The observed associations varied with local climate and urban infrastructure.
Introduction
Emerging evidence from epidemiological studies suggests a compelling association between rainfall events and adverse health outcomes,12 particularly the transmission of infectious diseases.3 The influence of rainfall events on non-communicable diseases, however, remains understudied.45 Current evidence, derived mainly from specific locations,67 has failed to capture the diverse health effects of regional variations in rainfall characteristics. Extreme rainfall events can result in rapid meteorological changes that potentially lead to an interaction with harmful aerosols,89 which can contribute to an increased risk for various cardiovascular610 and respiratory conditions.10
Several challenges hinder comprehensive understanding of the association between rainfall events and health outcomes. Firstly, most studies defined rainfall events using traditional centile methods,111213 and the variation in centile thresholds used across studies has led to mixed results. This problem, compounded by limitations such as restricted geographical locations and short study durations,6 has prevented consensus from being reached. Secondly, these traditional centile methods failed to reflect local rainfall patterns and did not follow the World Meteorological Organization criteria for identifying important events14: High intensity (substantial rainfall within a brief period), low frequency, and short duration.14 Thirdly, the association between rainfall intensity and health outcomes was likely non-linear. Although moderate rainfall can mitigate summer heat and help reduce air pollution, which in turn might reduce some environmental health risks,15 high intensity rainfall events of low frequency and short duration may have a particularly harmful effect on health, as such events can result in rapid meteorological changes, increased run-off, overwhelmed infrastructures,16 increased breeding of pathogens, and increased risk of different kinds of pollutants, including particulate matter and chemicals,1317 which can potentially exacerbate health conditions.
Climate change is intensifying the variability in precipitation1819 and extreme rainfall events both daily and overall.2021 An awareness of the effects of these extreme events is crucial for understanding the complex health consequences of climate change. Using an intensity-duration-frequency model to combine the three rainfall indices (high intensity, low frequency, short duration) with mortality data from 34 countries or regions, we estimated associations between mortality (all cause, cardiovascular, and respiratory) and rainfall events with different return periods (the expected average time between occurrences of an extreme event of a certain magnitude) and crucial effect modifiers, including climatic, socioeconomic, and urban environmental conditions. The findings should help improve our understanding of the effect of rainfall on health and guide health services to better respond to rainfall related health risks in the context of climate change.
Methods
Study population
We obtained data on mortality from the Multi-Country Multi-City (MCC) Collaborative Research Network database (https://mccstudy.lshtm.ac.uk/), which is described elsewhere.222324 Daily mortality records, provided by local statistical authorities within each participating country or region, were categorised using ICD-9 and ICD-10 (international classification of diseases, ninth and 10th revisions, respectively) codes. We obtained records for non-external mortality (ICD-10 codes: A00-R99; ICD-9 codes: 001-799) or all cause mortality depending on data availability, as well as specific data for cardiovascular diseases (ICD-10 codes: I00-I99; ICD-9 codes: 390-459) and respiratory diseases (ICD-10 codes: J00-J99; ICD-9 codes: 460-519) across a total of 743 locations from 1980 to 2020. These locations were selected based on the representativeness of the mortality data in relation to demographic characteristics of populations across diverse conditions in different countries and regions.2526 The appendix in our previous study provides a list of these locations and details on the sources of mortality data for each.24
Exposure data
Our primary climate data consisted of surface precipitation rates on an hourly scale (hourly rainfall amount) sourced from the land component of the fifth generation of the European reanalysis (ERA5-Land) dataset.27 This dataset integrates satellite and ground observations through advanced assimilation and modelling to estimate climate variables globally on land. We obtained data on a regular 0.1°×0.1° grid for 1980 to 2020 for each location, although the available years differed by location. The data were aggregated to daily scales, accounting for different time zones. Numerous studies have validated the precision of this precipitation data across various regions.2829 Furthermore, considering the distinct effects of rainfall and snowfall,30 we also collected data on snow depth. To differentiate between rain and snow in subsequent analyses, we treated rainfall days with a daily mean temperature <0°C and a snow depth >0 as non-precipitation days, assuming such conditions likely indicated snowfall or mixed precipitation rather than pure rainfall.30
In addition, we sourced precipitation data from two additional datasets for sensitivity tests. The MSWEP (multi-source weighted-ensemble precipitation) dataset31 provided data with a three hour temporal resolution, matching ERA5-Land’s spatial resolution for 1980-2020. We also utilised the IMERG (integrated multi-satellite retrievals for global precipitation measurement) dataset, which offers half hourly precipitation measurements from 2000 to 2021.32 IMERG uses >10 satellites with passive microwave radiometers to deliver calibrated global estimates of precipitation, a method widely recognised for its reliability.33 From these datasets, for each city in our study, we extracted data from the grid point closest to the city centre’s coordinates. We also obtained daily mean temperature and relative humidity data from the MCC database, originally sourced from local monitoring stations. Further details on data collection and assessment procedures are provided elsewhere.26 Finally, we collected daily (24 hour average) pollutant concentrations, specifically particulate matter with aerodynamic diameters ≤2.5 μm (PM2.5) and eight hour maximum ozone, from local monitoring stations. These air pollution data were available only for 372 locations across 19 countries or regions, as described previously.23
Definitions of rainfall events
Rainfall patterns often do not follow a normal distribution owing to periods of no rain leading to many zero values and extreme weather events causing skewness.3435 Geographical factors also contribute to uneven distributions across regions.36 Therefore, centile methods that are suitable for defining heat waves or cold spells based on daily temperature data may be inadequate for rainfall analysis. Impactful events are usually intense and rare, necessitating a definition that accounts for both intensity and frequency.
Following the World Meteorological Organization’s guideline,14 we used the extreme value theorem to assess the threshold of extreme rainfall events at daily scale with a specific return period. Return periods are the expected average times between occurrences of an extreme event of a certain magnitude. This approach involves considering the frequency, duration, and intensity of extreme events at the same time.37 Specifically, we used the regional frequency analysis framework38 with the generalised extreme value function to model the distribution of annual maximum daily cumulative rainfall across all locations.39 Using this fitted generalised extreme value model, we established thresholds for rainfall events with return periods of one year, two years, and five years for each group of cities with similar climatic and precipitation conditions. For example, if the cumulative probability of daily cumulative rainfall reached 50% based on the fitted distribution, this level of rainfall was then associated with a two year return period. This does not imply that such events happen regularly every two years but rather describes their probability of occurrence over a long term horizon. The supplementary methods provide details on the methodological and validation approaches. The methodology is commonly employed in research examining the health effects of extreme climatic conditions40 and is extensively used in hydrological studies.37414243
Furthermore, to illustrate the gradation of relative risk associated with increasing rainfall intensity from relatively moderate to extreme, we calculated the percentage of daily accumulated rainfall exceeding the two year return period threshold as a continuous variable, denoted as PE2 (fig 1). This measure provides a continuous index for quantifying the relative intensity of extreme rainfall events.
As an example, if the threshold for a two year return period was 50 mm and a city experienced 70 mm in a day, the index value would be calculated as 40%, meaning that on that day the city experienced rainfall that was 40% greater than expected for a two year return period. We selected this threshold based on our preliminary research findings, where no significant effect was observed for rainfall events with a one year and two year return period.43
Statistical analysis
To investigate the association between extreme rainfall events and daily mortality, we employed an established two stage analytical method.264445 This approach follows the methods detailed in previous studies,4345 ensuring consistency and comparability in our analysis.
In the first stage of our analysis, we utilised quasi-Poisson regression models to assess the city specific associations between extreme rainfall events and daily mortality. The quasi-Poisson regression model was adopted as it has been found to fit data for daily mortality count well in many related studies.4546 For the main model, we employed distributed lag models to evaluate the cumulative impact of exposure to extreme rainfall as a binary variable. Based on return period analysis, we identified a rainfall threshold for three different return periods (one year, two years, and five years). Days with rainfall below certain threshold were assigned a value of 0, indicating either no rainfall or precipitation not meeting the criteria. Conversely, days exceeding this threshold were assigned a value of 1. We modelled the lag structure using a natural cubic B spline with four degrees of freedom and two internal knots, which were evenly spaced in the log scale.47 The lag period was set to 14 days, which was informed by our preliminary findings in east Asia.43 To control the potential confounding effects of non-optimal temperatures, we used a distributed lag non-linear model to introduce a cross basis function of daily mean temperature in the main model. This distributed lag non-linear model comprised a quadratic B spline with three internal knots placed at the 10th, 75th, and 90th centiles of location specific temperature distributions.26 Additionally, the model included a lag response curve characterised by a natural cubic spline with three internal knots evenly spaced on the log scale with the same 14 day lag setting.26 We also included a natural cubic B spline of time with seven degrees of freedom per year26 and an indicator of the day of the week to control for the long term and seasonal trends and weekly variations.26 In the second stage, we pooled the city specific estimates to obtain overall relative risk based on coefficients and covariance matrices of the results from each city using random effects meta-analysis with restricted maximum likelihood estimation.48 Specifically, these associations were pooled to generate an overall level and were based on several key factors related to the rainfall effect to explore potential effect modifications that have not been adjusted for, including local climate types (according to the Köppen climate classification, a widely used global climate zoning method49), annual cumulative precipitation, the average annual standard deviation of daily cumulative precipitation for all rainfall days, economic conditions (reflected in the country’s income class, 2010 version), population density, and the average percentage of vegetation coverage (2010 version). The supplementary file provides detailed information on the data sources for each factor. We examined heterogeneity using Cochran’s Q test and the I2 statistic.50
For the analysis of PE2, we employed a similar process, with one key distinction. In the first step, we replaced the distributed lag model function of extreme rainfall events (treated as a binary variable) with a distributed lag non-linear model cross basis function of PE2, treating it as a continuous variable. For this function, we utilised the same distributed lag non-linear model settings as those applied for daily mean temperature. This modification in the model allowed us to capture the cumulative effects of varying rainfall intensities while maintaining methodological consistency.
To strengthen the reliability of our analysis, we focused on locations that had matched a minimum of two days of rainfall events with a five year return period during the study period to exclude locations lacking adequate exposure events. Additionally, we excluded rainfall events if during event days no hourly total rainfall exceeded 4 mm, according to the World Meteorological Organization’s criteria.51 Finally, to ensure the robustness of our findings, we conducted a series of sensitivity analyses (see supplementary methods), including adjustments to some key model settings, the incorporation of alternative precipitation datasets, and controlling for air pollution and flood events.
Effect modifications
To investigate how various factors might modify the health effects of extreme rainfall events across different populations, we applied random effects meta-regressions with maximum likelihood estimation. This approach allowed us to compare the cumulative relative risks across different levels of potential effect modifiers.45 The potential effect modifiers examined include the classification of local climate types, the amount of total annual precipitation, the yearly average standard deviation of precipitation on rainfall days, economic status, population density, and average vegetation coverage.
Patient and public involvement
This study utilised data from worldwide death registries, which inherently restricts direct involvement of patients or their relatives owing to data privacy and security regulations across different counties. In addition, no funding was allocated for patient or public involvement, and given the nature of the data, such interaction was not feasible. At the same time, given the nature of this study it was not appropriate or feasible to involve members of the public directly. Moreover, the research was primarily data driven, with no clear mechanism for public participation in the analysis process. Importantly, this study was partially inspired by observed respiratory health effects after extreme rainfall events, which reflect the issues patients face.
Results
After excluding 98 locations owing to insufficient exposure events, where the number of matched extreme rainfall days with a five year return period was less than two days, this study encompassed 645 locations across 34 countries or regions. Results of the Kolmogorov-Smirnov test, a method used to verify the reliability of the results (see supplementary methods for details), showed that the generalised extreme value distributions were accepted at a 5% significance level for all included locations (see supplementary figure S1). Our analysis included 109 954 744 all cause, 31 164 161 cardiovascular, and 11 817 278 respiratory deaths. During the study period, we identified a total of 50 913 rainfall events with a one year return period, 8362 events with a two year return period, and 3301 events with a five year return period. Figure 2 illustrates the distribution patterns of the threshold of rainfall events with the three levels of return periods. Locations that experienced higher daily cumulative rainfall within the same return period were predominantly in South America, east Asia, and South East Asia. Supplementary table S1 provides a statistical overview of the included locations, detailing the economic classification, number of included locations, and number of matched extreme rainfall events during the study period for each country or region.
We identified a global positive association between all cause mortality and extreme rainfall events with a five year return period (fig 3). However, shorter return periods of one year or two years were not associated with increased risks of all cause mortality. Overall, across all locations, a day of extreme rainfall with a five year return period was associated with a cumulative relative risk across 1-14 lag days of 1.08 (95% confidence interval (CI) 1.05 to 1.11) for daily all cause mortality. Cause specific analyses showed that rainfall events with a two year return period were associated with increased daily respiratory mortality (1.14, 1.05 to 1.23), whereas no significant effect was observed for the same period for cardiovascular mortality (1.01, 0.98 to 1.04). Furthermore, five year return period events were associated with increased risk for both cardiovascular (1.05, 1.02 to 1.08) and respiratory mortality (1.29, 1.19 to 1.39), with the risk for respiratory mortality notably higher than for cardiovascular mortality (P<0.05). Conversely, rainfall events with a one year return period showed no effect on either cardiovascular or respiratory mortality. The lag-response function results (see supplementary figure S2) indicate that the peak relative risk occurred on the day of the event (lag 0), with a substantial decrease observed by lag 3 to 4 days. Sensitivity analyses (see supplementary file) further confirmed the robustness of these findings.
Figure 4 illustrates the distinct intensity-response associations between PE2 and all cause mortality. The association displayed a negative trend when PE2 was below a critical threshold. However, the curve showed a consistent upward trajectory once PE2 surpassed 50%, with a substantial positive effect emerging as PE2 surpassed about 100%, indicative of rainfall intensities more than twice the threshold for the two year return period. The pattern for cardiovascular mortality closely mirrored that of all cause mortality, which initially presented a negative association but shifted to a positive association at higher PE2 values. The association for respiratory mortality escalated steeply and more pronouncedly with increasing PE2 values.
Regional specific analyses revealed variations in the relative risk of daily mortality associated with different rainfall events (fig 3). For events with a two year return period, the relative risk ranged from 0.85 (95% CI 0.74 to 0.99) in Australia to 1.23 (1.08 to 1.38) in South America. This variation became more pronounced for events with a five year return period, with the relative risk in South America reaching 1.43 (1.19 to 1.67). When regions were compared, the estimated impact of extreme rainfall events on mortality in South America was notably higher than the global average for both two year and five year return periods. North America’s estimates were close to the global average. Additionally, increases in mortality risk (P<0.05) were observed with five year return period events in east Asia, North America, and South America. When health outcomes were compared, risk ratios for respiratory mortality exceeded those for cardiovascular mortality across most regions.
The analyses of effect modification (see table 1 and supplementary figures S3-S7) revealed that climate types modified the mortality risk associated with extreme rainfall events with a five year return period. The relative risk of all cause mortality was found to be higher in locations characterised by scarce vegetation coverage or lower variability in daily precipitation (P<0.05). Additionally, a statistically significant association with all cause mortality was observed predominantly in lower income countries. The association with respiratory mortality was also more pronounced in locations with scarce vegetation coverage and lower variability in daily precipitation. Evidence that average annual precipitation and population density modified these associations is, however, limited.
Discussion
In this study, we performed an analysis of mortality risks associated with extreme rainfall using a large dataset covering 34 countries or regions and employing a uniform methodology to define rainfall events based on intensity, frequency, and duration. Our analysis identified associations between extreme rainfall events and daily mortality from all causes, as well as specifically from cardiovascular and respiratory diseases. Furthermore, the intensity-response functions revealed a protective effect of daily rainfall at moderate to heavy intensities, which transitioned to a marked increase in mortality risk at extreme rainfall intensities. Moreover, the substantial health effects of these extreme rainfall events were found to vary by climate types. Extreme rainfall events were more severe in locations with lower variability of precipitation or scarce vegetation coverage. Our study utilised a comprehensive database encompassing information from multiple countries, territories, and continents. The analysis incorporated death records from numerous locations from 1980 to 2020, providing a global perspective on the effect of extreme rainfall events on health.
Only a few studies have explored the association between exposure to rainfall events and risk of mortality from non-external causes, with results showing notable inconsistencies.525354 This variation might be attributed to the lack of a consistent definition of rainfall events and the limited geographical scope of most previous investigations. Similar to our results, many studies have found harmful effects of extreme precipitation. A study in rural Bangladesh observed a 0.6% increase in daily mortality for every additional millimetre of rainfall beyond 100 mm per day in the 5-19 age group.53 A study in Burkina Faso found an increased risk of malaria related death associated with precipitation levels ≥453 mm.1 Additionally, a notable connection has been observed between thunderstorms accompanying extreme rainfall and an increase in respiratory morbidity and mortality, especially in people with asthma.5556 However, limited studies have shown either insignificant effects or protective effects of rainfall, aligning with certain aspects of our findings. A recent study in tropical regions of Africa found no significant association between extreme rainfalls and incidence of respiratory diseases among children.11 A study in India observed a significant reduction in risk of all cause mortality during the first week after extreme rainfall.2 Similarly, research conducted in rural areas of Ecuador found that extreme rainfall events were associated with a decreased incidence of diarrhoea after wet periods.57 These findings align with our observations that some levels of precipitation may not only be non-detrimental but could also offer certain protective effects.
Biological mechanisms
Although the biological mechanisms linking rainfall intensity with health outcomes are not fully elucidated, several plausible explanations may be helpful to clarify the complex association observed in our study. Firstly, moderate to high intensity rainfall may exert protective effects through two primary mechanisms: Improvement in air quality (rainfall can decrease the concentration of PM2.5 particles in the atmosphere,58 potentially mitigating its adverse effects on all cause, respiratory, and cardiovascular mortality), and behavioural changes (rainfall may alter people’s daily patterns, leading to more time spent indoors).59 This reduction in direct exposure to outdoor air pollution and non-optimal temperatures could explain the protective effect observed on various health outcomes when intensity is not extreme. Secondly, as rainfall intensity increases, the initial protective effects may be affected by a cascade of negative impacts: Critical disruptions to resources, physiological impacts, and indirect effects.
Critical resource disruptions—Intense rainfall can cause major disruptions to healthcare access. It might damage infrastructure, leading to power outages and hindering essential medical services, both of which can particularly affect older people and those with chronic illnesses, who often require consistent access to drugs and healthcare.60 Intense rainfall can compromise water and food quality by allowing pathogenic microorganisms to reproduce and spread,61 leading to diarrhoeal diseases, infections, and associated mortality.6263
Physiological effects—Increased humidity levels facilitate the growth of airborne pathogens,64 potentially triggering allergic reactions and respiratory problems, particularly in vulnerable individuals.65 Rapid shifts in atmospheric pressure can reduce the partial pressure of oxygen in the body,6667 leading to cardiovascular and respiratory complications.60 Furthermore, major temperature fluctuations often accompanying rainfall days can affect both respiratory and cardiovascular systems.68
Indirect effects—Extreme rainfall can have profound effects on mental health, inducing stress and anxiety that may exacerbate pre-existing mental health conditions69 and indirectly contributing to increased total mortality from non-external causes. Overall, these pathways elucidate the complex associations between rainfall intensity and health outcomes observed in our study. Moderate rainfall (PE2 <100%) shows protective effects, whereas intense rainfall (PE2 >100%) leads to important harmful effects as negative factors accumulate. Supplementary analysis of the five year return period threshold also corroborated these findings, showing a weaker protective effect and a rapid transition to harmful impacts beyond 50% of the threshold rather than 100% as seen with PE2 (see supplementary figure S7). Moreover, the rapid onset of these negative factors during extreme rainfall days explained the observed lag pattern (see supplementary figure S2), with health effects most pronounced on the first day and then gradually diminishing.
We found several factors that could influence the variability in the health effects of extreme rainfalls. Notably, the associations were modified by climate types and baseline variations in rainfall intensity. Many studies have identified an association between a region’s resilience to extreme climate events and its natural conditions.7071 Areas with stable rainfall patterns, such as arid regions or areas with low variability in rainfall, often have weaker adaptive capacities for heavy rainfall. A recent study suggested that the risk of rainfall associated infectious diseases is more pronounced in arid areas, likely due to the more major environmental changes caused by brief heavy downpours.72 Furthermore, many studies have highlighted a negative correlation between per capita gross domestic product and mortality resulting from extreme climate events.7374 This finding suggests that socioeconomic development enhances urban resilience and the ability to respond rapidly to extreme rainfall related disasters. Conversely, unemployment and social inequality can exacerbate health problems during such events. It has also been noted that the effect of extreme rainfall is influenced by vegetation coverage, with higher vegetation coverage correlating with greater adaptive capacity for managing the consequences of extreme rainfall.75 In contrast, we found no significant effect modifications linked to annual average precipitation and population density. Greater cumulative rainfall does not necessarily equate to a higher frequency of extreme rain events. Adaptation to extreme rainfall is most noticeable in areas with more frequent and variable rainfalls.76 Additionally, population density does not indicate socioeconomic conditions or residential environments,77 and thus it does not serve as a reliable predictor of resilience to disasters. Regional variations in the effects of extreme rainfall can be attributed to several factors, in addition to the varying adequacy of identified total events. In east Asia, higher intensity of extreme rainfall events (fig 2) coupled with rapid urbanisation and reduced green space coverage (see supplementary figure S5) could increase overall vulnerability. In North America, varying green space coverage (see supplementary figure S5) combined with ageing infrastructure, inadequate updates to drainage and water management systems, and uneven distribution of resources, are important explanatory factors for the observed effects.78 In South America, major effects likely arise from limited access to medical resources.79
Policy implications
Our findings have some practical implications. Policy makers should be aware of the implications of rainfall related risks on public health, encompassing both non-communicable and infectious diseases. Hospitals and clinicians need vigilance when monitoring and managing patients during the rainy season, and researchers need to investigate the exact mechanisms that could then be targeted by systematic protective action. Considering the health risks associated with the modification of rainfall by local climatic conditions and vegetation coverage, localised or health education initiatives centred on communities would be imperative to increase awareness about the risks to health of severe rainfalls, and increasing vegetation coverage is a reasonable adaptation strategy. A multidisciplinary approach and integrated strategies are essential to navigate public health interventions to cope with extreme rainfall. This necessitates not only the fortification of meteorological observation and early warning systems but also the implementation of robust measures across urban planning and social security sectors to bolster the resilience of communities and individuals against disasters.
Limitations of this study
This study has some limitations. Firstly, although the analysis included 34 major counties or regions on six continents, the analysed locations were mainly located in east Asia, Europe, and North America, with fewer in Latin America and Africa. This limitation restricts the global representativeness of the estimated risks. Furthermore, the study areas were all urban centres, therefore the associations between extreme rainfall and mortality might differ in rural areas, and we only used mortality statistics from nearly 700 cities globally. Although our results provide a broad, population level perspective, they do not capture granular details of individual factors such as age, sex, race, urban/rural residence, or specific clinical settings. In addition, this study could be limited by its ecological study design as with most previous epidemiological studies on the impact of climate related events.2680 As such, this study was unable to analyse specific personal factors, interventions, or clinic based practices. A more granular study focusing on individual features or clinical settings would be a valuable next step in this line of research. Secondly, the accuracy of our assessment outcomes may be influenced by potential exposure misclassification arising from our reliance on modelled output precipitation data.81 This issue might further be compounded by the limited spatial resolution inherent in the rainfall databases employed. We incorporated a diverse array of rainfall databases to sensitively test our findings. Thirdly, in a global study spanning several decades, diagnostic or coding errors in health data may have occurred. The effect of these errors on our results is likely to be minimal, however, given most data used in this study were from urban centres and from government data sources. Finally, our analysis captured only a part of the health effects of heavy rainfalls. It notably lacked data such as hospital admission records and health outcome data from rural areas and excluded data on external causes of mortality. These limitations highlight the need for further comprehensive studies to evaluate a broader range of health effects associated with extreme rainfall events.
Conclusions
Our multicountry analysis provided epidemiological evidence for associations between heavy rainfall events and mortality from all causes and cardiovascular or respiratory causes, based on a global dataset and standard statistical methods. Notably, the intensity-response curves for the health effect of heavy rainfalls showed a non-linear trend, transitioning from a protective effect at heavy rainfall levels to a risk of major harm as rainfall intensity became extreme. Furthermore, the significant effects of extreme events were modified by climate types and were most pronounced in locations characterised by low variability in precipitation or scarce vegetation coverage. These findings underscore the need for comprehensive public health strategies, developed through collaboration among meteorological, public health, and urban planning sectors. Such strategies are crucial to mitigate the broad health effects of extreme rainfall. This is especially important considering the well established trend of increasing short term rainfall intensity as a result of climate change.
What is already known on this topic
Climate change is intensifying the frequency and severity of short term rainfall events
Current research on rainfall related health risks primarily focuses on infectious diseases and severe storms, overlooking broader impacts such as those on cardiovascular and respiratory health, and how varying rainfall intensities affect these condition
Existing evidence is predominantly derived from specific locations, and regional variations in rainfall intensity, duration, frequency, and timing can lead to diverse health impacts across different geographical areas
What this study adds
This study found that local factors such as climate conditions, climate type, and vegetation coverage have a potential influence on rainfall related all cause and cardiovascular and respiratory mortality
These results can enhance understanding of the broad and complex health impacts of climate change, guiding the development of localised adaptation measures for policy makers
Ethics statements
Ethical approval
Not required.
Data availability statement
Health outcome data have been collected within the MCC (Multi-Country Multi-City) Collaborative Research Network (https://mccstudy.lshtm.ac.uk/) under a data sharing agreement and cannot be made publicly available. Researchers can refer to MCC participants listed as coauthors for information on accessing the data for each country. The time series rainfall data from ERA5-Land, MSWEP, and IMERG are publicly accessible from their respective websites (https://cds.climate.copernicus.eu; https://www.gloh2o.org/mswep/; https://disc.gsfc.nasa.gov/datasets/GPM_3IMERGHHL_06/summary)
Footnotes
Contributors: AS is the senior author. CH designed the study. CH, SB, VH. KC, and SZ coordinated the work, conducted the statistical analysis, and took the lead in drafting the manuscript and interpreting the results. AG, MB, HK, DR, BA, JS, FS, AMV-C, and FD provided substantial scientific input in interpreting the results and drafting the manuscript. YH, JJ, NR, YG, ST, PM, MC, PS, EL, JK, HO, EI, MP, PG, AZ, YK, MHD, EA, AO, JK, SR, SP, GC, XS, SS, JM, JH, LH, NS, FA, HKi, WL, MH, AT, CI, BF, MR, YLG, SCP, SO, SL, AZ, TND, and DVD provided the data and contributed to the revision of the submitted version of the manuscript. CH and AS are the guarantors. The corresponding author attests that all listed authors meet authorship criteria and that no others meeting the criteria have been omitted.
Funding: CH was supported by the Alexander von Humboldt Foundation. The funders had no role in considering the study design or in the collection, analysis, interpretation of data, writing of the report, or decision to submit the article for publication.
Competing interests: All authors have completed the ICMJE uniform disclosure form at www.icmje.org/disclosure-of-interest/ and declare: Support from the Alexander von Humboldt Foundation no financial relationships with any organizations that might have an interest in the submitted work in the previous three years; no other relationships or activities that could appear to have influenced the submitted work.
Transparency: The lead authors (CH and AS) affirm that the manuscript is an honest, accurate, and transparent account of the study being reported; that no important aspects of the study have been omitted; and that any discrepancies from the study as planned have been explained.
Dissemination to participants and related patient and public communities: We plan to disseminate the findings through multiple channels. A press release will be coordinated with our communication department, as well as with the communication departments of some key coauthors. The results will also be shared with students during lectures. Additionally, we will present this study at some well known academic conferences and will continue to engage with journalists, particularly when discussing the health impacts of climate change. We are also exploring the possibility of a blog post or other accessible formats to increase public engagement with the study.
Provenance and peer review: Not commissioned; externally peer reviewed.
Publisher’s note: Published maps are provided without any warranty of any kind, either express or implied. BMJ remains neutral with regard to jurisdictional claims in published maps.
This is an Open Access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/.