Assessing seasonality and the role of its potential drivers in environmental epidemiology: a tutorial

The Multi-Country Multi-City Collaborative Research Network: An international research consortium investigating environment, climate, and health

Research on the health risks of environmental factors and climate change requires epidemiological evidence on associated health risks at a global scale. Multi-center studies offer an excellent framework for this purpose, but they present various methodological and logistical problems. This contribution illustrates the experience of the Multi-Country Multi-City Collaborative Research Network, an international collaboration working on a global research program on the associations between environmental stressors, climate, and health in a multi-center setting. The article illustrates the collaborative scheme based on mutual contribution and data and method sharing, describes the collection of a huge multi-location database, summarizes published research findings and future plans, and discusses advantages and limitations. The Multi-Country Multi-City represents an example of a collaborative research framework that has greatly contributed to advance knowledge on the health impacts of climate change and other environmental factors and can be replicated to address other research questions across various research fields.

Incidence and seasonality of Kawasaki disease in children in the Philippines, and its association with ambient air temperature

Background

Despite an unknown cause, Kawasaki disease (KD) is currently the primary leading cause of acquired heart disease in developed countries in children and has been increasing in recent years. Research efforts have explored environmental factors related to KD, but they are still unclear especially in the tropics. We aimed to describe the incidence of KD in children, assess its seasonality, and determine its association with ambient air temperature in the National Capital Region (NCR), Philippines from January 2009 to December 2019.

Methods

Monthly number of KD cases from the Philippine Pediatric Society (PPS) disease registry was collected to determine the incidence of KD. A generalized linear model (GLM) with quasi-Poisson regression was utilized to assess the seasonality of KD and determine its association with ambient air temperature after adjusting for the relevant confounders.

Results

The majority of KD cases (68.52%) occurred in children less than five years old, with incidence rates ranging from 14.98 to 23.20 cases per 100,000 population, and a male-to-female ratio of 1.43:1. Seasonal variation followed a unimodal shape with a rate ratio of 1.13 from the average, peaking in March and reaching the lowest in September. After adjusting for seasonality and long-term trend, every one-degree Celsius increase in the monthly mean temperature significantly increased the risk of developing KD by 8.28% (95% CI: 2.12%, 14.80%). Season-specific analysis revealed a positive association during the dry season (RR: 1.06, 95% CI: 1.01, 1.11), whereas no evidence of association was found during the wet season (RR: 1.10, 95% CI: 0.95, 1.27).

Conclusion

We have presented the incidence of KD in the Philippines which is relatively varied from its neighboring countries. The unimodal seasonality of KD and its linear association with temperature, independent of season and secular trend, especially during dry season, may provide insights into its etiology and may support enhanced KD detection efforts in the country.

Heat-related mortality and ambulance transport after a power outage in the Tokyo metropolitan area

Background

Air conditioners can prevent heat-related illness and mortality, but the increased use of air conditioners may enhance susceptibility to heat-related illnesses during large-scale power failures. Here, we examined the risks of heat-related illness ambulance transport (HIAT) and mortality associated with typhoon-related electricity reduction (ER) in the summer months in the Tokyo metropolitan area.

Methods

We conducted event study analyses to compare temperature-HIAT and mortality associations before and after the power outage (July to September 2019). To better understand the role of temperature during the power outage, we then examined whether the temperature-HIAT and mortality associations were modified by different power outage levels (0%, 10%, and 20% ER). We computed the ratios of relative risks to compare the risks associated with various ER values to the risks associated without ER.

Results

We analyzed the data of 14,912 HIAT cases and 74,064 deaths. Overall, 93,200 power outage cases were observed when the typhoon hit. Event study results showed that the incidence rate ratio was 2.01 (95% confidence interval [CI] = 1.42, 2.84) with effects enduring up to 6 days, and 1.11 (95% CI = 1.02, 1.22) for mortality on the first 3 days after the typhoon hit. Comparing 20% to 0% ER, the ratios of relative risks of heat exposure were 2.32 (95% CI = 1.41, 3.82) for HIAT and 0.95 (95% CI = 0.75, 1.22) for mortality.

Conclusions

A 20% ER was associated with a two-fold greater risk of HIAT because of summer heat during the power outage, but there was little evidence for the association with all-cause mortality.

Seasonal variation of Covid-19 incidence and role of land surface and air temperatures: a case study in the west of Iran

.In this study, we assessed the impact of satellite-based Land Surface Temperature (LST) and Air Temperature (AT) on covid-19. First, we spatio-temporally kriged the LST and applied bias correction. The epidemic shape, timing, and size were compared after and before adjusting for the predictors. Given the non-linear behavior of a pandemic, a semi-parametric regression model was used. In addition, the interaction effect between the predictors and season was assessed. Before adjusting for the predictors, the peak happened at the end of hot season. After adjusting, it was attenuated and slightly moved forward. Moreover, the Attributable Fraction (AF) and Peak to Trough Relative (PTR) were % 23 (95% CI; 15, 32) and 1.62 (95%CI; 1.34, 1.97), respectively. We found that temperature might have changed the seasonal variation of covid-19. However, given the large uncertainty after adjusting for the variables, it was hard to provide conclusive evidence in the region we studied.

Seasonality of mortality under climate change: a multicountry projection study

BACKGROUND

Climate change can directly impact temperature-related excess deaths and might subsequently change the seasonal variation in mortality. In this study, we aimed to provide a systematic and comprehensive assessment of potential future changes in the seasonal variation, or seasonality, of mortality across different climate zones.

METHODS

In this modelling study, we collected daily time series of mean temperature and mortality (all causes or non-external causes only) via the Multi-Country Multi-City Collaborative (MCC) Research Network. These data were collected during overlapping periods, spanning from Jan 1, 1969 to Dec 31, 2020. We projected daily mortality from Jan 1, 2000 to Dec 31, 2099, under four climate change scenarios corresponding to increasing emissions (Shared Socioeconomic Pathways [SSP] scenarios SSP1-2.6, SSP2-4.5, SSP3-7.0, and SSP5-8.5). We compared the seasonality in projected mortality between decades by its shape, timings (the day-of-year) of minimum (trough) and maximum (peak) mortality, and sizes (peak-to-trough ratio and attributable fraction). Attributable fraction was used to measure the burden of seasonality of mortality. The results were summarised by climate zones.

FINDINGS

The MCC dataset included 126 809 537 deaths from 707 locations within 43 countries or areas. After excluding the only two polar locations (both high-altitude locations in Peru) from climatic zone assessments, we analysed 126 766 164 deaths in 705 locations aggregated in four climate zones (tropical, arid, temperate, and continental). From the 2000s to the 2090s, our projections showed an increase in mortality during the warm seasons and a decrease in mortality during the cold seasons, albeit with mortality remaining high during the cold seasons, under all four SSP scenarios in the arid, temperate, and continental zones. The magnitude of this changing pattern was more pronounced under the high-emission scenarios (SSP3-7.0 and SSP5-8.5), substantially altering the shape of seasonality of mortality and, under the highest emission scenario (SSP5-8.5), shifting the mortality peak from cold seasons to warm seasons in arid, temperate, and continental zones, and increasing the size of seasonality in all zones except the arid zone by the end of the century. In the 2090s compared with the 2000s, the change in peak-to-trough ratio (relative scale) ranged from 0·96 to 1·11, and the change in attributable fraction ranged from 0·002% to 0·06% under the SSP5-8.5 (highest emission) scenario.

INTERPRETATION

A warming climate can substantially change the seasonality of mortality in the future. Our projections suggest that health-care systems should consider preparing for a potentially increased demand during warm seasons and sustained high demand during cold seasons, particularly in regions characterised by arid, temperate, and continental climates.

FUNDING

The Environment Research and Technology Development Fund of the Environmental Restoration and Conservation Agency, provided by the Ministry of the Environment of Japan.

Ambient temperature and seasonal variation in inpatient mortality from respiratory diseases: a retrospective observational study

Background

The seasonal fluctuation in mortality and hospital admissions from respiratory diseases, with a winter peak and a summer trough, is widely recognized in extratropical countries. However, little is known about the seasonality of inpatient mortality and the role of ambient temperature remains uncertain. We aimed to analyse the association between ambient temperature and in-hospital mortality from respiratory diseases in the provinces of Madrid and Barcelona, Spain.

Methods

We used data on daily hospitalisations, weather (ie, temperature and relative humidity) and air pollutants (ie, PM2.5, PM10, NO2 and O3) for the Spanish provinces of Madrid and Barcelona during 2006-2019. We applied a daily time-series quasi-Poisson regression in combination with distributed lag non-linear models (DLNM) to assess, on the one hand, the seasonal variation in fatal hospitalisations and the contribution of ambient temperature, and on the other hand, the day-to-day association between temperature and fatal hospital admissions. The analyses were stratified by sex, age and primary diagnostic of hospitalisation.

Findings

The study analysed 1 710 012 emergency hospital admissions for respiratory diseases (mean [SD] age, 60.4 [31.0] years; 44.2% women), from which 103 845 resulted in in-hospital death (81.4 [12.3] years; 45.1%). We found a strong seasonal fluctuation in in-hospital mortality from respiratory diseases. While hospital admissions were higher during the cold season, the maximum incidence of inpatient mortality was during the summer and was strongly related to high temperatures. When analysing the day-to-day association between temperature and in-hospital mortality, we only found an effect for high temperatures. The relative risk (RR) of fatal hospitalisation at the 99th percentile of the distribution of daily temperatures vs the minimum mortality temperature (MMT) was 1.395 (95% eCI: 1.211-1.606) in Madrid and 1.612 (1.379-1.885) in Barcelona. In terms of attributable burden, summer temperatures (June-September) were responsible for 16.2% (8.8-23.3) and 22.3% (15.4-29.2) of overall fatal hospitalisations from respiratory diseases in Madrid and Barcelona, respectively. Women were more vulnerable to heat than men, whereas the results by diagnostic of admission showed heat effects for acute bronchitis and bronchiolitis, pneumonia and respiratory failure.

Interpretation

Unless effective adaptation measures are taken in hospital facilities, climate warming could exacerbate the burden of inpatient mortality from respiratory diseases during the warm season.

Funding

European Research Council Consolidator Grant EARLY-ADAPT, European Research Council Proof-of-Concept Grants HHS-EWS and FORECAST-AIR.