Research

8 March 2023

Short-Term Association between Sulfur Dioxide and Mortality: A Multicountry Analysis in 399 Cities

Authors: Edward O’Brien, Pierre Masselot, Francesco Sera, Dominic Roye, Susanne Breitner, Chris Fook Sheng Ng, Micheline de Sousa Zanotti Stagliorio Coelho, … Show All … , Joana Madureira, Aurelio Tobias, Ana Maria Vicedo-Cabrera, Michelle L. Bell, Eric Lavigne, Haidong Kan, Antonio Gasparrini https://orcid.org/0000-0002-2271-3568, and MCC Collaborative Research NetworkAuthors Info & Affiliations

Publication: Environmental Health Perspectives

Volume 131, Issue 3

CID: 037002

https://doi.org/10.1289/EHP11112

PDF

Abstract

Background:

Epidemiological evidence on the health risks of sulfur dioxide (

{SO}_{2}

) is more limited compared with other pollutants, and doubts remain on several aspects, such as the form of the exposure–response relationship, the potential role of copollutants, as well as the actual risk at low concentrations and possible temporal variation in risks.

Objectives:

Our aim was to assess the short-term association between exposure to

{SO}_{2}

and daily mortality in a large multilocation data set, using advanced study designs and statistical techniques.

Methods:

The analysis included 43,729,018 deaths that occurred in 399 cities within 23 countries between 1980 and 2018. A two-stage design was applied to assess the association between the daily concentration of

{SO}_{2}

and mortality counts, including first-stage time-series regressions and second-stage multilevel random-effect meta-analyses. Secondary analyses assessed the exposure–response shape and the lag structure using spline terms and distributed lag models, respectively, and temporal variations in risk using a longitudinal meta-regression. Bi-pollutant models were applied to examine confounding effects of particulate matter with an aerodynamic diameter of

less than or equal to 10 micrometers

(

particulate matter begin subscript 10 end subscript

) and

2.5 micrometers

(

particulate matter begin subscript 2.5 end subscript

), ozone, nitrogen dioxide, and carbon monoxide. Associations were reported as relative risks (RRs) and fractions of excess deaths.

Results:

The average daily concentration of

{SO}_{2}

across the 399 cities was

11.7 micrograms per meter cubed

, with 4.7% of days above the World Health Organization (WHO) guideline limit (

40 micrograms per meter cubed

, 24-h average), although the exceedances occurred predominantly in specific locations. Exposure levels decreased considerably during the study period, from an average concentration of

19.0 micrograms per meter cubed

in 1980–1989 to

6.3 micrograms per meter cubed

in 2010–2018. For all locations combined, a

10 microgram per meter cubed

increase in daily

{SO}_{2}

was associated with an RR of mortality of 1.0045 [95% confidence interval (CI): 1.0019, 1.0070], with the risk being stable over time but with substantial between-country heterogeneity. Short-term exposure to

{SO}_{2}

was associated with an excess mortality fraction of 0.50% [95% empirical CI (eCI): 0.42%, 0.57%] in the 399 cities, although decreasing from 0.74% (0.61%, 0.85%) in 1980–1989 to 0.37% (0.27%, 0.47%) in 2010–2018. There was some evidence of nonlinearity, with a steep exposure–response relationship at low concentrations and the risk attenuating at higher levels. The relevant lag window was 0–3 d. Significant positive associations remained after controlling for other pollutants.

Discussion:

The analysis revealed independent mortality risks associated with short-term exposure to

{SO}_{2}

, with no evidence of a threshold. Levels below the current WHO guidelines for 24-h averages were still associated with substantial excess mortality, indicating the potential benefits of stricter air quality standards. https://doi.org/10.1289/EHP11112

Introduction

Sulfur dioxide (

{SO}_{2}

) is an important air pollutant linked with increased health risks.^1,2 It originates largely from the combustion of fossil fuels to generate electricity and transportation.^1,3,4

{SO}_{2}

is also released during industrial processes, mainly the production of metals, such as copper, and other chemical plants, or from emissions from fuel combustion in shipping, whereas a very small percentage occurs naturally from volcanoes and fissures. In many developed nations, the desulfurization of cars and power plants has dramatically reduced emissions and the related population exposure, but the same cannot be said of many developing countries, where levels remain relatively very high.² In 2005, the World Health Organization (WHO) published their recommendation for air quality standards, stating that the average

{SO}_{2}

concentration over a 24-h period should not exceed

20 micrograms per meter cubed

500 micrograms per meter cubed

over 10 min.⁵ The WHO guidelines were revised in 2021, with the 24-h limit increased to

40 micrograms per meter cubed

following a new criterion based on the distribution of daily

{SO}_{2}

concentrations and the corresponding limit of annual averages.⁶

Short-term mortality risks of

{SO}_{2}

have been assessed in several ecological studies primarily based on time-series data. Early multicity studies were conducted in Europe,^7,8 the United States,⁹ and East Asia.¹⁰ More recently, large investigations were performed in China, where the

{SO}_{2}

concentrations far exceed those of most high-income countries.^11,12 A recent meta-analysis of 67 eligible studies systematically reviewed the evidence and provided pooled estimates of the association.¹³ When restricting the unit of analysis to 24-h averages of

{SO}_{2}

, a

10 microgram per meter cubed

increment of

{SO}_{2}

was associated with an increase of 0.59% (95% CI: 0.46%, 0.71%) in all-cause mortality. The association remained when controlling for particulate matter (PM), but not for nitrogen dioxide (

{NO}_{2}

) or ozone (

O_{3}

). Moreover, there was no evidence of an exposure threshold below which no risk can be assumed.¹³

However, several gaps in knowledge still exist, as discussed in a comprehensive report from the U.S. Environmental Protection Agency (EPA).¹⁴ In addition to the uncertainty regarding the potential confounding effects from copollutants mentioned above, limited evidence is available on other aspects of the short-term association between exposure to

{SO}_{2}

and mortality risks. For instance, there is no conclusive evidence on the shape of the exposure–response relationships and possible nonlinearities or about the presence of more complex temporal dependencies and lagged associations. Further, it is still unclear whether results from studies from China are generalizable elsewhere, or if the risk shows geographical heterogeneity. More importantly, the published analyses assessed the association at relatively high exposure ranges. It is unclear to what extent the risk can be extrapolated at lower concentrations, for instance, below current air quality guidelines. This information is critical for revising air quality limits using evidence-based processes.

In this contribution, we address these limitations and present results from an investigation of the short-term mortality risks associated with

{SO}_{2}

exposure using data from 399 cities in 23 countries across the globe. The analysis used advanced study designs and statistical methods to characterize the associations of interest, whereas the large database and broad exposure contrasts offered enough statistical power to assess geographical variations and risks at low exposure levels.

Methods

Data Collection

The data were collected within the Multi-Country Multi-City (MCC) Collaborative Research Network, an international collaboration investigating environmental stressors and their impacts on human health.¹⁵ The MCC database has been used in previous publications that evaluated associations between air pollutants and mortality.^15–20 Mortality data were gathered from the local health authorities and were represented as daily counts of all causes [International Classification of Diseases, Manual of the International Statistical Classification of Diseases, Injuries and Causes of Death (ICD-9),²¹ codes 0–799], if available, or nonexternal [International Statistical Classification of Diseases and Related Health Problems, Tenth Revision (ICD-10),²² codes A00–R99] deaths. Nonexternal causes of death exclude intentional and unintentional injury, poisoning (including drug overdose), and complications of medical or surgical care.

The data set also contains 24-h average city-level concentrations of

{SO}_{2}

{NO}_{2}

O_{3}

, carbon monoxide (CO), and PM with an aerodynamic diameter of

less than or equal to 10 micrometers

(

particulate matter begin subscript 10 end subscript

) and

2.5 micrometers

(

particulate matter begin subscript 2.5 end subscript

), in addition to daily temperature. All pollutants were harmonized using the unit of micrograms per meter cubed, except CO, which was harmonized with milligrams per meter cubed. The country-specific data sets generally include all the major cities; however, this analysis was restricted to the 399 locations within 23 countries with

{SO}_{2}

data available and at least 365 d of measurement. The United States of America (USA) was divided into nine regions to account for the large heterogeneity in

{SO}_{2}

values, for a total of 31 areas. The geographical location of each city and the mean

{SO}_{2}

concentrations are displayed in Figure 1. Detailed information on data collection is reported in the Supplemental Material, “Information on country-specific datasets.”

Figure 1. Geographical location of the 399 urban areas and related average annual concentrations of ${SO}_{2}$ (in ${μ g / m}^{3}$ ) within the study period 1980–2018. Data can be found in Table S1. Note: ${SO}_{2}$ , sulfur dioxide.

Statistical Methods

Main Model

We applied a two-stage procedure to analyze the short-term association between

{SO}_{2}

and mortality. In the first stage, we performed city-specific time-series analyses using a quasi-Poisson generalized linear model with distributed lag terms.^23,24 The city-specific model included an indicator for the day of the week to account for within-week variation and a natural spline function with 7 degrees of freedom per year to control for long-term trends and seasonal variations. Air temperature was modeled with a distributed lag nonlinear model (DLNM), composed of a quadratic B-spline with three knots placed at the 10th, 75th, and 90th percentiles for the exposure–response and a step function with strata lags of 0 and 1–3 d for the lag–response. In the main model,

{SO}_{2}

was modeled assuming a linear exposure–response relationship of the moving average computed over lag 0–3 d.

In the second stage, we combined the city-specific estimates using a multilevel random-effect meta-analysis fitted with restricted maximum likelihood (REML) and nested random effects defined by city and country.²⁵ The pooled estimate represents the global average

{SO}_{2} -mortality association

, whereas city and country-specific estimates were derived as best linear unbiased predictions (BLUPs) at the corresponding aggregation level. The BLUPs use information from pooled associations to make more accurate location-specific estimates by borrowing information from the whole sample, especially for cities/countries with higher uncertainty, while at the same time accounting for heterogeneity in risks.^25,26 All estimates are reported as the relative risk (RR) for a

10 microgram per meter cubed

increase in

{SO}_{2}

, with corresponding 95% confidence intervals (CIs). Heterogeneity was reported as

I^{2}

statistics and tested with the Cochran’s

Q

test.^25–27

Secondary Analyses

We performed a series of secondary analyses. First, we explored potential nonlinear exposure–response shapes and more complex lag structures of the

{SO}_{2} - mortality relationship

, extending the model first using a quintic polynomial, and then using a distributed lag model (DLM) with a natural spline with knots at lag 1 and 3 plus an intercept over lag 0–7 d. The polynomial parameterization was adopted to decrease the sensitivity of the estimates to different ranges of

{SO}_{2}

concentrations because polynomial terms are not local and are, instead, defined across the whole variable range. In both extensions, city-specific estimates of the multiparameter associations were pooled using a multivariate multilevel meta-analysis.²⁸

Second, we evaluated possible changes in risk over time by subsetting the city-specific data and performing the first-stage model in multiple subperiods, splitting the time series into

approximately 5 year

intervals. City and period-specific estimates were pooled in a longitudinal multilevel meta-regression using time (defined as the midyear of each subperiod) as a continuous fixed-effect term.²⁸

Third, we assessed the potential confounding effect of other pollutants in bi-pollutant models, where

particulate matter begin subscript 10 end subscript

particulate matter begin subscript 2.5 end subscript

O_{3}

(8-h daily maximum),

{NO}_{2}

, and CO entered the model linearly, using a moving average of lag 0–3 d. These choices were informed by previous studies.^12,29,30 Only one other pollutant was controlled for at a time because of the high correlation between pollutants.^14,31 These models were fitted both with and without adjustment in the subset of cities, providing measurements for both pollutants. Information on the levels of copollutants across cities and countries has been provided in previous articles.^16–20

Computation of Excess Mortality

Finally, using the main model, we computed the excess mortality associated in the short term with exposure to

{SO}_{2}

in each city, using a previously described method.³² Briefly, the cumulative RR within lag 0–3 was applied to compute the excess daily deaths, adopting a forward perspective using the standard formula

(1 - e x p (- β_{j} {(x_{j t} - c)}_{+})) \times d_{j t}

for continuous exposure, as in previous analyses.¹⁷ In the formula,

β_{j}

represents the log-RR for a unit increase in

{SO}_{2}

, defined as the country-specific BLUP for city

j

, and

x_{j t}

and

d_{j t}

are the corresponding

{SO}_{2}

levels at day

t

and the average daily mortality in the same and next 3 d, respectively. The term

{(x_{j t} - c)}_{+}

represents the exceedance in

{SO}_{2}

concentration above a limit

c

. We used

lowercase italic c equals 0

and

lowercase italic c equals 40

to compute the burden attributable to short-term exposure to

{SO}_{2}

in total and above the WHO guideline, respectively. The results are reported as fraction of excess deaths, both in total and for levels below the WHO guideline, together with 95% empirical CIs (eCIs).

All analyses were conducted in R (version 4.2.2; R Development Core Team), using the dlnm and mixmeta packages. The R code for the original analysis and for performing a reproducible example using simulated data is available in a GitHub repository (https://github.com/gasparrini/MCC-SO2).

Results

Descriptive Analysis

The analysis included 43,729,018 deaths across 399 cities, in 23 countries (separating the USA into nine regions), with an average period of 14.5 y. Table 1 shows the total deaths, number of cities, period of analysis, and levels of

{SO}_{2}

across cities in each country. The list of cities, together with basic information on study periods and

{SO}_{2}

exposure levels, is provided in Table S1. The

{SO}_{2}

exposure was widely heterogeneous both within and between countries, with an average across cities of

11.7 micrograms per meter cubed

and an average 5th–95th percentile range of

2.6 to 29.3 micrograms per meter cubed

(see the map in Figure 1 for the geographical distribution). The highest mean concentrations were found in cities of China and Iran (50.1 and

98.9 micrograms per meter cubed

, respectively), and the lowest concentrations were in Portugal and Estonia (2.6 and

3.5 micrograms per meter cubed

, respectively). Across the whole set of cities and periods, 4.7% of days registered concentrations higher than the WHO limits of

40 micrograms per meter cubed

.⁶

Table 1 Descriptive statistics reported by country (with the United States separated into nine regions): number of cities, total days, and deaths in the time series, distribution of sulfur dioxide (

{SO}_{2}

) across cities, and percentage of days with

{SO}_{2}

concentrations

greater than 40 micrograms per meter cubed

Country	Cities ( $n$ )	Period	City-days ( $n$ )	Deaths ( $n$ )	Average mean and 5th–95th percentile range of ${SO}_{2}$ ( ${μ g / m}^{3}$ )	Percentage of days $greater than 40 micrograms per meter cubed$ (%)^a
Canada	24	2000–2015	109,953	1,741,439	5.6 (0.3–17.7)	1.1
USA-Central	30	1985–2006	216,457	3,724,610	19.6 (3.5–49.0)	9.9
USA-NE Central	14	1985–2006	76,137	1,518,119	10.2 (1.0–29.2)	2.4
USA-Northeast	42	1985–2006	291,202	6,912,465	19.8 (3.6–49.8)	10.2
USA-Northwest	5	1985–2006	19,994	394,004	14.6 (3.8–32.3)	2.3
USA-NW Central	2	1985–2006	10,894	69,007	2.8 (0.2–9.7)	0.4
USA-South	17	1985–2006	105,338	1,513,772	8.7 (0.9–26.1)	2.7
USA-Southeast	29	1985–2006	173,040	2,945,570	9.5 (1.5–25.9)	2.2
USA-Southwest	9	1985–2006	50,766	794,800	8.0 (0.8–22.6)	1.6
USA-West	15	1985–2006	81,604	3,367,343	5.0 (0.2–16.9)	0.2
Puerto Rico	1	2009–2016	2,877	26,161	4.6 (0.0–9.4)	0.7
Brazil	1	1997–2011	5,076	909,305	12.3 (4.0–26.3)	0.6
Colombia	1	1998–2013	5,800	423,344	20.9 (5.3–41.6)	6.4
Ecuador	1	2014–2018	1,819	44,369	16.7 (8.0–26.0)	0.2
Peru	1	2010–2014	1,489	148,775	12.8 (3.4–27.8)	0.1
Estonia	4	2002–2018	18,239	96,455	2.6 (0.2–10.1)	0.6
Finland	1	1994–2014	7,663	153,166	9.1 (2.1–25.9)	1.0
UK	31	1990–2016	157,489	3,823,644	8.6 (0.5–23.6)	1.9
Czech Republic	1	1994–2009	5,835	213,706	14.6 (2.0–55.7)	8.6
Germany	12	1993–2015	70,542	2,098,705	7.9 (1.9–23.8)	2.0
Romania	8	2008–2016	14,964	221,816	8.4 (3.2–17.7)	0.4
Switzerland	4	1995–2013	24,229	151,898	5.6 (0.7–16.9)	0.2
Portugal	6	1995–2018	28,860	847,377	3.5 (0.5–9.9)	0.3
Spain	48	2002–2014	197,008	1,480,869	5.5 (2.6–10.4)	0.1
Iran	1	2002–2015	5,025	683,739	98.9 (17.0–279.1)	70.7
China	15	1996–2015	23,139	1,181,405	50.1 (14.8–124.2)	45.8
Japan	46	1980–2015	133,940	4,382,591	6.0 (2.2–12.6)	0.3
South Korea	7	1999–2015	42,979	1,658,788	15.5 (7.5–29.1)	1.3
Thailand	18	1999–2008	52,980	722,911	10.1 (3.0–20.7)	2.3
Taiwan	3	1994–2014	22,981	1,208,118	15.4 (6.1–32.1)	3.6
Australia	2	2000–2009	6,609	270,751	5.1 (0.8–18.8)	0.8
All MCC countries^b	399	1980–2018	1,964,928	43,729,018	11.7 (2.6–29.3)	4.7

Note: The analysis includes data from 399 cities within the study period 1980–2018 from the Multi-Country Multi-City (MCC) Collaborative Research Network.

Current limit of daily concentration of

{SO}_{2}

in the World Health Organization guidelines.⁶

Data shown represents the statistics for all the areas included above.

Figure 2 displays the distribution of city-specific average

{SO}_{2}

concentrations over the years, revealing a strong attenuation in exposure levels within the study period, although the comparison should account for the different sample of cities/countries contributing to each interval. Actual figures of

{SO}_{2}

levels and percentage exceedances of WHO limits aggregated by decade are reported in Table S2. Results show that the average daily concentration decreased from

19.0 micrograms per meter cubed

in 1980–1989 to

6.3 micrograms per meter cubed

in 2010–2018, and the corresponding percentage of days

greater than 40 micrograms per meter cubed

from 11.7% to 1.3%.

Figure 2. Box plot of the distribution of the average concentration of ${SO}_{2}$ (in ${μ g / m}^{3}$ ) across cities for each year. The horizontal line identifies the current limit of daily concentration of ${SO}_{2}$ in the WHO guidelines ( $40 micrograms per meter cubed$ ). The y-axis is represented in a logarithmic scale. The analysis includes data from 399 cities within the study period 1980–2018. Note that a different set of countries contributes to each study period (see Table 1 for details). Note: ${SO}_{2}$ , sulfur dioxide; WHO, World Health Organization.

Association between $S O_{2}$ and Mortality

The risk associations estimated from the main model using a linear exposure–response relationship and a moving average of lag 0–3 are illustrated in Figure 3, with the pooled RR and related country-specific BLUPs (see Table S3 for the numeric data). On average across all cities and countries, each

10 microgram per meter cubed

increase of

{SO}_{2}

was associated with an RR of mortality of 1.0045 (95% CI: 1.0019, 1.0070). Although country-specific estimates were less precise, there is evidence of substantial heterogeneity (

uppercase italic i squared equals 45.8 percent

, Cochran

uppercase italic q lowercase italic p less than 0.001

), with the RRs ranging from 0.9968 (95% CI: 0.9883, 1.0054) in Finland to 1.0234 (95% CI: 1.0168, 1.0300) in Brazil, and a few nonsignificant negative estimates.

Figure 3. Country-specific and pooled relative risks (RRs, with 95% CIs) for mortality corresponding to a $10 microgram per meter cubed$ increase in ${SO}_{2}$ over lag 0–3 d. The analysis includes data from 399 cities within the study period 1980–2018. Data can be found in Table S3. Note: CI, confidence interval; ${SO}_{2}$ , sulfur dioxide.

Exposure and Lag–Response Relationships

Results from secondary analyses are displayed in Figure 4. Figure 4A shows the estimate of the pooled

{SO}_{2} -mortality relationship

from the flexible model allowing nonlinear exposure–response associations. The graph is augmented with the estimated log-linear relationship from the main model and a bar representing the number of cities with

{SO}_{2}

measurements within the related exposure range. The figure indicates some evidence of nonlinearity, with a steep increase in risk and an attenuation at high concentrations. It is unclear if this supralinear shape resulted from lower risks at high

{SO}_{2}

exposures or if it is attributable to a different sample of countries contributing at various ranges, given that only 23.3% of cities were exposed to levels

greater than 150 micrograms per meter cubed

. In any case, the nonlinear parameterization confirms the evidence of mortality risks for exposures below the WHO limit of

40 micrograms per meter cubed

. A comparison of country-specific exposure–response relationships (as BLUPs) estimated using the linear and nonlinear models is presented in Figure S1. The results of the second modeling extension with the application of a DLM over lag 0–7 to assess the lag structure are displayed in Figure 4B. The graph suggests no risk of same-day exposure to

{SO}_{2}

, with the risk then increasing in the next 3–4 d. This analysis indicated that the main model with the moving average of lag 0–3 can capture these lagged associations and can provide valid estimates of the association.

Figure 4. Secondary analysis on the short-term association between ${SO}_{2}$ (per $10 microgram per meter cubed$ increase) and mortality. (A) Pooled exposure–response curve obtained using a linear term (dashed line) and a quintic polynomial (continuous line, with 95% confidence intervals), with a bar representing the percentage of studies contributing to the specific exposure range. (B) Pooled lag–response curve obtained using a natural spline with knots at lags 1 and 3, plus intercept. (C) Analysis of temporal variation of the pooled relative risk (RR) associated with a $10 microgram per meter cubed$ increase in ${SO}_{2}$ over lag 0–3 d. Shaded areas represent the 95% confidence intervals. The analysis includes data from 399 cities within the study period 1980–2018. Summary data on city-specific exposure distributions can be found in Table S1. Note: ${SO}_{2}$ , sulfur dioxide.

Analysis of Temporal Variation in Risk

The extension of the two-stage design for assessing potential temporal changes in risk involved the analysis by subperiod and the pooling of estimates using time as a continuous meta-predictor. The results are reported in Figure 4C, which shows the pooled mortality RR for a

10 microgram per meter cubed

increase in

{SO}_{2}

along the years predicted from the longitudinal meta-analytical model. The graph suggests little evidence of variation in time in the short-term association, with a

p

-value of the time term equal to 0.67.

Bi-Pollutant Analysis

Results from bi-pollutant analyses are provided in Table 2, with information on the number of cities contributing to the analysis of each copollutant provided in Table S4. The comparison suggests that an independent risk associated with

{SO}_{2}

remained even after adjustment for each of the five other pollutants, although with some variations indicative of partial confounding. Specifically, the estimated risks seemed to attenuate after controlling for

particulate matter begin subscript 10 end subscript

and

{NO}_{2}

, whereas they increased when including

particulate matter begin subscript 2.5 end subscript

in the model. Results were only negligibly affected by control for

O_{3}

and CO.

Table 2 Relative risk with

95% confidence interval

) [

relative risk 95% confidence interval

)] associated with a

10 microgram per meter cubed

increase in sulfur dioxide (

{SO}_{2}

, in

{μ g / m}^{3}

) over lag 0–3 d with and without adjustment for each copollutant in selected cities with both measurements within the study period 1980–2018.

Copollutant	Cities ( $n$ )	$RR \pm 95 % CI$
Copollutant	Cities ( $n$ )	Without adjustment	With adjustment
$particulate matter begin subscript 10 end subscript$	265	1.0045 (1.0017, 1.0074)	1.0028 (1.0004, 1.0052)
$particulate matter begin subscript 2.5 end subscript$	217	1.0028 (1.0005, 1.0050)	1.0056 (1.0028, 1.0084)
$O_{3}$ 8h	309	1.0044 (1.0023, 1.0065)	1.0040 (1.0017, 1.0062)
${NO}_{2}$	358	1.0046 (1.0019, 1.0073)	1.0032 (1.0012, 1.0052)
CO	302	1.0038 (1.0016, 1.0060)	1.0036 (1.0023, 1.0048)

Note: Number of cities included in the bi-pollutant models are reported in Table S4. CO, carbon monoxide;

{NO}_{2}

, nitrogen dioxide;

O_{3}

, ozone;

particulate matter begin subscript 2.5 end subscript

, particulate matter with an aerodynamic diameter of

2.5 micrometers

;

particulate matter begin subscript 10 end subscript

, particulate matter with an aerodynamic diameter of

less than or equal to 10 micrometers

Excess Mortality

Finally, Table 3 depicts the excess mortality fraction associated with

{SO}_{2}

exposure computed from the main model, assuming a linear relationship both by decade and within the whole study period of each area. Overall across the 399 cities,

{SO}_{2}

was found to be associated with an excess increase of 0.50% (95% eCI: 0.42%, 0.57%) of the total deaths. This fraction showed a strong decrease in time and variation across regions, consistent with the reduction in

{SO}_{2}

levels and the geographical variation in risk described above. Specifically, the excess mortality overall decreased from 0.74% (0.61%, 0.85%) in 1980–1989 to 0.37% (0.27%, 0.47%) in 2010–2018. The corresponding quota of excess deaths attributable to levels

below 40 micrograms per meter cubed

is shown in Table S5, indicating that 93.8% of the excess was due to exposures below the WHO limits on average across the study period.

Table 3 Excess mortality fraction with 95% empirical confidence interval [% (95% eCI)] attributable to short-term exposure to sulfur dioxide (

{SO}_{2}

, per

10 micrograms per meter cubed

) by country (separating the USA into nine regions) and decade.

Country/region	1980–1989	1990–1999	2000–2009	2010–2019	Full period
Canada	—	—	0.24 (0.12, 0.37)	0.11 (0.05, 0.17)	0.19 (0.10, 0.29)
USA-Central	0.82 (0.56, 1.05)	0.57 (0.40, 0.73)	0.40 (0.28, 0.51)	—	0.58 (0.40, 0.73)
USA-NE Central	0.71 (0.28, 1.16)	0.50 (0.19, 0.82)	0.35 (0.12, 0.58)	—	0.50 (0.19, 0.82)
USA-Northeast	1.09 (0.79, 1.36)	0.77 (0.56, 0.97)	0.55 (0.40, 0.70)	—	0.78 (0.57, 0.98)
USA-Northwest	0.15 ( $negative 0.52$ , 0.80)	0.15 ( $negative 0.48$ , 0.73)	0.06 ( $negative 0.32$ , 0.39)	—	0.13 ( $negative 0.46$ , 0.67)
USA-NW Central	0.20 ( $negative 0.28$ , 0.65)	0.22 ( $negative 0.31$ , 0.72)	0.08 ( $negative 0.11$ , 0.25)	—	0.16 ( $negative 0.22$ , 0.53)
USA-South	0.63 (0.45, 0.83)	0.49 (0.35, 0.64)	0.33 (0.22, 0.44)	—	0.47 (0.33, 0.61)
USA-Southeast	0.11 ( $negative 0.06$ , 0.26)	0.07 ( $negative 0.02$ , 0.17)	0.05 ( $negative 0.01$ , 0.11)	—	0.07 ( $negative 0.02$ , 0.15)
USA-Southwest	0.55 (0.21, 0.88)	0.47 (0.18, 0.74)	0.26 (0.09, 0.43)	—	0.40 (0.16, 0.63)
USA-West	1.07 (0.86, 1.29)	0.80 (0.66, 0.94)	0.76 (0.64, 0.88)	—	0.84 (0.70, 0.99)
Puerto Rico	—	—	0.04 ( $negative 0.32$ , 0.40)	0.05 ( $negative 0.41$ , 0.50)	0.05 ( $negative 0.40$ , 0.49)
Brazil	—	3.75 (2.69, 4.70)	2.85 (2.04, 3.57)	1.69 (1.21, 2.12)	2.87 (2.05, 3.60)
Colombia	—	0.14 ( $negative 2.02$ , 2.26)	0.10 ( $negative 1.55$ , 1.73)	0.03 ( $negative 0.49$ , 0.56)	0.09 ( $negative 1.32$ , 1.47)
Ecuador	—	—	—	0.76 ( $negative 1.07$ , 2.55)	0.76 ( $negative 1.07$ , 2.55)
Peru	—	—	—	0.83 ( $negative 0.41$ , 2.17)	0.83 ( $negative 0.41$ , 2.17)
Estonia	—	—	0.08 ( $negative 0.12$ , 0.28)	0.05 ( $negative 0.05$ , 0.14)	0.06 ( $negative 0.07$ , 0.19)
Finland	—	$negative 0.35$ ( $negative 1.25$ , 0.59)	$negative 0.28$ ( $negative 1.00$ , 0.47)	$negative 0.27$ ( $negative 0.98$ , 0.46)	$negative 0.29$ ( $negative 1.07$ , 0.50)
UK	—	$negative 0.49$ ( $negative 0.79$ , $negative 0.14$ )	$negative 0.14$ ( $negative 0.21$ , $negative 0.08$ )	$negative 0.06$ ( $negative 0.10$ , $negative 0.02$ )	$negative 0.24$ ( $negative 0.39$ , $negative 0.09$ )
Czech Republic	—	$negative 0.03$ ( $negative 1.65$ , 1.69)	$negative 0.01$ ( $negative 0.36$ , 0.38)	—	$negative 0.01$ ( $negative 0.88$ , 0.91)
Germany	—	0.38 (0.16, 0.59)	0.15 (0.07, 0.22)	0.09 (0.02, 0.15)	0.24 (0.11, 0.35)
Romania	—	—	0.54 ( $negative 0.27$ , 1.31)	0.39 ( $negative 0.13$ , 0.89)	0.42 ( $negative 0.16$ , 0.98)
Switzerland	—	0.82 (0.15, 1.57)	0.41 (0.14, 0.71)	0.21 (0.07, 0.36)	0.48 (0.13, 0.86)
Portugal	—	0.42 (0.04, 0.81)	0.27 (0.07, 0.46)	0.11 (0.04, 0.18)	0.24 (0.06, 0.40)
Spain	—	—	0.82 (0.60, 1.04)	0.60 (0.46, 0.74)	0.73 (0.55, 0.92)
Iran	—	—	$negative 0.37$ ( $negative 4.30$ , 3.31)	$negative 0.13$ ( $negative 1.45$ , 1.16)	$negative 0.26$ ( $negative 3.02$ , 2.35)
China	—	0.58 (0.11, 1.07)	1.55 (1.12, 1.98)	1.68 (0.79, 2.55)	1.48 (1.04, 1.95)
Japan	0.17 ( $negative 0.18$ , 0.51)	0.11 ( $negative 0.10$ , 0.31)	0.06 ( $negative 0.04$ , 0.16)	0.11 (0.02, 0.19)	0.11 ( $negative 0.05$ , 0.25)
South Korea	—	2.32 (1.72, 2.91)	1.49 (1.09, 1.88)	1.37 (1.00, 1.74)	1.49 (1.09, 1.88)
Thailand	—	1.21 (0.80, 1.61)	1.07 (0.74, 1.40)	—	1.08 (0.75, 1.41)
Taiwan	—	$negative 0.13$ ( $negative 0.87$ , 0.53)	$negative 0.08$ ( $negative 0.54$ , 0.35)	$negative 0.06$ ( $negative 0.42$ , 0.27)	$negative 0.09$ ( $negative 0.59$ , 0.37)
Australia	—	—	0.20 ( $negative 0.10$ , 0.50)	—	0.20 ( $negative 0.10$ , 0.50)
All MCC countries	0.74 (0.61, 0.85)	0.46 (0.38, 0.53)	0.50 (0.40, 0.60)	0.37 (0.27, 0.47)	0.50 (0.42, 0.57)

Note: The analysis includes data from the Multi-Country Multi-City (MCC) Collaborative Research Network for 399 cities within the study period 1980–2018. Estimates based on the main model assuming a linear exposure–response relationship and a moving average of lag 0–3 d. —, not applicable.

Discussion

To our knowledge, this study represents the most extensive epidemiological assessment of the short-term mortality risks associated with exposure to

{SO}_{2}

, investigating the relationship using a large database that includes almost

44 million

deaths from 399 cities in 23 countries across 5 continents. We found an overall increased short-term risk, with an

RR equal to 1.0045

(95% CI: 1.0019, 1.0070) per

10 microgram per meter cubed

increase in

{SO}_{2}

, although with evidence of heterogeneity across countries and cities. This translated to an annual excess corresponding to 0.5% of the total mortality on average across the study period. Stratified by decades, the analysis indicates a strong decrease in the health impacts, in line with the reduction in the concentration levels of

{SO}_{2}

although this result should be interpreted with caution given the differential temporal coverage across countries. The application of more flexible models suggested some evidence of nonlinearity, with steeper RR estimates at low exposure levels, and a complex lag structure with the maximal risk arising 1–3 d after exposure. Independent associations were still measurable even after controlling for copollutants, and there was no evidence that the risk associated with a given exposure level had changed over time.

This study provides an essential contribution to the literature on the mortality risks associated with short-term

{SO}_{2}

exposure. The pooled RR estimates are consistent with previous epidemiological studies, although slightly lower when compared with figures published in the multicity studies and the meta-analysis described in the introduction.^7–13 The difference can be due to the selection of locations, which in our study represent a broader sample and a wider geographical area. Our results strengthen the evidence on the association of short-term

{SO}_{2}

exposure with total mortality, which was determined as suggestive but not sufficient to infer a causal relationship in the U.S. EPA report mentioned above.¹⁴ One of the main reasons described in the report for such a conclusion was the limited knowledge about the potential confounding by other pollutants. Our analysis confirms that risk estimates were somewhat sensitive to control for copollutants but that there was still strong evidence of the association when

particulate matter begin subscript 10 end subscript

particulate matter begin subscript 2.5 end subscript

O_{3}

{NO}_{2}

, and CO were each included in the model. In addition, our contribution addresses other knowledge gaps on the association, such as the shape of the exposure–response relationship, the lag patterns, and possible temporal variations in risks.

There are several biological pathways and mechanisms through which

{SO}_{2}

can lead to higher mortality risks. The body of evidence from animal, experimental, and epidemiological studies is deemed sufficient to suggest a causal relationship between exposure and effects on the respiratory system.¹³ Specifically, short-term exposure to

{SO}_{2}

is known to induce neural reflex responses, release of inflammatory mediators, and modulation of allergic inflammation, leading to several end points that include bronchoconstriction and increased airway responsiveness.³³ The impacts are known to be more severe in susceptible individuals, such as those with asthma.³⁴ Previous studies have also assessed potential risks for cardiovascular outcomes.

{SO}_{2}

was shown to cause drops in measures of baroreflex sensitivity and cardiac vagal control, as well as increases in plasma fibrinogen, oxidative stress, and blood viscosity in young adults.³⁵ However, the body of evidence is not sufficient and consistent enough to establish a causal relationship.¹⁴

Another important result of this study is the evidence of the potential public health benefits achievable with more stringent air quality policies. As stated above, we found that short-term exposure to

{SO}_{2}

was associated with an excess risk of mortality. Although the impact has decreased in time consistent with the reduction in concentration levels, exposure to

{SO}_{2}

was still linked with a considerable number of excess deaths in recent years. More importantly, the majority of the additional deaths were linked with exposures sustained on days in which

{SO}_{2}

levels were at or below the current daily WHO threshold, which was recently increased from 20 to

40 micrograms per meter cubed

. These results, in line with previous multicountry studies on other pollutants,^16–20 indicate the presence of a considerable risk even at low exposure ranges that is also associated with a substantial health burden in countries that comply with the current WHO guidelines. These results therefore support the efforts to enforce national and international policy guidelines and to consider the opportunity to revise the limits downward.

An important aspect of this assessment is the application of state-of-the-art study design and statistical methods on a large multicountry database. Given the large sample of locations covering such a wide geographic area, we were able to obtain consistent evidence using a uniform model of a short-term association between exposure to

{SO}_{2}

and all-cause mortality. The use of advanced analytical techniques applied in the first stage, including time-series regression and distributed lag models, allows a nuanced characterization of complex exposure–lag–response relationships. Similarly, the analysis of hundreds of locations and the use of extended meta-analytical techniques provide both pooled and city- and country-specific estimates and offer a comprehensive geographical and temporal comparison across various regions of the globe.

Some limitations must be acknowledged. First, although we were able to provide risk summaries across four inhabited continents, our results should not be considered truly global estimates given that some areas, such as Africa, South America, and the Middle East, were underrepresented or not assessed. Moreover, the study was restricted to urban populations, with several countries represented by a small number of cities, and therefore it cannot be entirely representative of the risks across whole populations. Notably, although the assessment offers a consistent comparison of risks across locations and found important differences in risks and impacts across countries, we did not attempt to characterize such heterogeneity. Part of this variability can be due to systemic differences concerning measurements from atmospheric monitors (type of station, proximity to the study area), study area boundaries, temporal coverage, and data processing, whereas the other part can be related to actual differences in susceptibility. This will be a topic for future research. In addition, the extension using polynomial terms indicates a degree of nonlinearity in the exposure–response relationship. However, it is challenging to disentangle to what extent this is due to heterogeneous risks in areas with high exposure ranges, and the main results are therefore reported from a model assuming a linear association. Further research on the shape of the exposure–response relationship and the risk at high

{SO}_{2}

levels is needed. Furthermore, although we assessed potential confounding effects in bi-pollutant models, we did not extend the analysis to evaluate possible synergistic effects between

{SO}_{2}

and other pollutants. Finally, the study was conducted using aggregated time-series data, preventing a more refined analysis of potential biological mechanisms and differential susceptibility patterns at the individual level.

Conclusions

This large multicountry study provides evidence of an independent short-term association between exposure to

{SO}_{2}

and all-cause mortality. The assessment indicates that even if current air quality guidelines for

{SO}_{2}

were enforced, many deaths would still occur, and additional health benefits could be attained by further lowering existing limits. These findings have important implications for the design of future health and environmental policy actions. More generally, they can contribute to the design and implementation of mitigation strategies to reduce the environmental risks and impacts on health in the context of climate change.

Acknowledgments

This study was supported by the Medical Research Council UK [MR/R013349/1 (to A.G.)] and the European Union’s Horizon 2020 Project Exhaustion [820655 (to A.G.)]. Additional grants supported individual researchers: H.K. was supported by the National Natural Science Foundation of China (92043301), J.M. was supported by a fellowship of Fundação para a Ciência e a Tecnlogia (SFRH/BPD/115112/2016), and A.T. was supported by Ministerio de Ciencia e Innovacion/Agencia Estatal de Investigación (MCIN/AEI/10.13039/501100011033; CEX2018-000794-S).

The Multi-City Multi-Country Collaborative Research Network group authorship includes G. Carrasco, B.-Y. Chen, A. Entezari, Y. Guo, Y.L. Guo, M. Hashizume, I.-H. Holobaca, C. Íñiguez, J.J.K. Jaakkola, H. Kim, W. Lee, S. Li, M. Maasikmets, F. Mayvaneh, B. Nunes, H. Orru, S. Osorio, M.S. Ragettli, N. Ryti, P.H.N. Saldiva, A. Schneider, J. Schwartz, S. Tong, and A. Zanobetti.

Article Notes

These authors contributed equally to this work.

†

A full list of MCC Network authors is displayed in the “Acknowledgments” section.

The authors declare they have nothing to disclose.

Supplementary Material

File (ehp11112.smcontents.508.pdf)

Download
93.42 KB

File (ehp11112.s001.acco.pdf)

Download
584.48 KB

References

Guttikunda SK, Tang Y, Carmichael GR, Kurata G, Pan L, Streets DG, et al. 2005. Impacts of Asian megacity emissions on regional air quality during spring 2001. J Geophys Res 110(D20):D20301, https://doi.org/10.1029/2004JD004921.

Google Scholar

Greenstone M. 2004. Did the Clean Air Act cause the remarkable decline in sulfur dioxide concentrations? J Environ Econ Manage. 47(3):585–611, https://doi.org/10.1016/j.jeem.2003.12.001.

Google Scholar

WHO (World Health Organization). Air pollution. https://www.who.int/health-topics/air-pollution [accessed 8 August 2019].

Google Scholar

Bell ML, Samet JM, Dominici F. 2004. Time-series studies of particulate matter. Annu Rev Public Health 25(1):247–280. https://pubmed.ncbi.nlm.nih.gov/15015920/, https://doi.org/10.1146/annurev.publhealth.25.102802.124329.

Google Scholar

WHO. 2006. WHO Air Quality Guidelines for Particulate Matter, Ozone, Nitrogen Dioxide and Sulfur Dioxide: Global Update 2005: Summary of Risk Assessment. https://apps.who.int/iris/handle/10665/69477 [accessed 25 November 2021].

Google Scholar

WHO. 2021. WHO Global Air Quality Guidelines: Particulate Matter (PM2.5 and PM10), Ozone, Nitrogen Dioxide, Sulfur Dioxide and Carbon Monoxide. https://apps.who.int/iris/handle/10665/345329 [accessed 25 November 2021].

Google Scholar

Katsouyanni K, Touloumi G, Spix C, Schwartz J, Balducci F, Medina S, et al. 1997. Short term effects of ambient sulphur dioxide and particulate matter on mortality in 12 European cities: results from time series data from the APHEA project. BMJ 314(7095):1658–1663. https://pubmed.ncbi.nlm.nih.gov/9180068/, https://doi.org/10.1136/bmj.314.7095.1658.

Google Scholar

Biggeri A, Baccini M, Bellini P, Terracini B. 2005. Meta-analysis of the Italian studies of short-term effects of air pollution (MISA), 1990–1999. Int J Occup Environ Health 11(1):107–122. https://pubmed.ncbi.nlm.nih.gov/15859200/, https://doi.org/10.1179/oeh.2005.11.1.107.

Google Scholar

Moolgavkar SH, McClellan RO, Dewanji A, Turim J, Luebeck EG, Edwards M. 2013. Time-series analyses of air pollution and mortality in the United States: a subsampling approach. Environ Health Perspect 121(1):73–78. https://pubmed.ncbi.nlm.nih.gov/23108284/, https://doi.org/10.1289/ehp.1104507.

Google Scholar

10.

Kan H, Wong CM, Vichit-Vadakan N, Qian Z, PAPA Project Teams. 2010. Short-term association between sulfur dioxide and daily mortality: the Public Health and Air Pollution in Asia (PAPA) study. Environ Res 110(3):258–264. https://pubmed.ncbi.nlm.nih.gov/20122685/, https://doi.org/10.1016/j.envres.2010.01.006.

Google Scholar

11.

Chen R, Huang W, Wong CM, Wang Z, Thach TQ, Chen B, et al. 2012. Short-term exposure to sulfur dioxide and daily mortality in 17 Chinese cities: the China Air Pollution and Health Effects Study (CAPES). Environ Res 118:101–106. https://pubmed.ncbi.nlm.nih.gov/22831556/, https://doi.org/10.1016/j.envres.2012.07.003.

Google Scholar

12.

Wang L, Liu C, Meng X, Niu Y, Lin Z, Liu Y, et al. 2018. Associations between short-term exposure to ambient sulfur dioxide and increased cause-specific mortality in 272 Chinese cities. Environ Int 117:33–39. https://pubmed.ncbi.nlm.nih.gov/29715611/, https://doi.org/10.1016/j.envint.2018.04.019.

Google Scholar

13.

Orellano P, Reynoso J, Quaranta N. 2021. Short-term exposure to sulphur dioxide (SO₂) and all-cause and respiratory mortality: a systematic review and meta-analysis. Environ Int 150:106434. https://pubmed.ncbi.nlm.nih.gov/33601225/, https://doi.org/10.1016/j.envint.2021.106434.

Google Scholar

14.

U.S. EPA (U.S. Environmental Protection Agency). 2017. Integrated Science Assessment (ISA) for Sulfur Oxides—Health Criteria (Final Report, Dec 2017). EPA/600/R-17/451. Washington, DC: U.S. EPA.

Google Scholar

15.

Multi-Country Multi-City (MCC) Collaborative Research Network. https://mccstudy.lshtm.ac.uk/ [accessed 28 August 2022].

Google Scholar

16.

Liu C, Chen R, Sera F, Vicedo-Cabrera AM, Guo Y, Tong S, et al. 2019. Ambient particulate air pollution and daily mortality in 652 cities. N Engl J Med 381(8):705–715. https://pubmed.ncbi.nlm.nih.gov/31433918/, https://doi.org/10.1056/NEJMoa1817364.

Google Scholar

17.

Vicedo-Cabrera AM, Sera F, Liu C, Armstrong B, Milojevic A, Guo Y, et al. 2020. Short term association between ozone and mortality: global two stage time series study in 406 locations in 20 countries. BMJ 368:m108. https://pubmed.ncbi.nlm.nih.gov/32041707/, https://doi.org/10.1136/bmj.m108.

Google Scholar

18.

Meng X, Liu C, Chen R, Sera F, Vicedo-Cabrera AM, Milojevic A, et al. 2021. Short term associations of ambient nitrogen dioxide with daily total, cardiovascular, and respiratory mortality: multilocation analysis in 398 cities. BMJ 372:n534. https://pubmed.ncbi.nlm.nih.gov/33762259/, https://doi.org/10.1136/bmj.n534.

Google Scholar

19.

Zhao Q, Guo Y, Ye T, Gasparrini A, Tong S, Overcenco A, et al. 2021. Global, regional, and national burden of mortality associated with non-optimal ambient temperatures from 2000 to 2019: a three-stage modelling study. Lancet Planet Health 5(7):e415–e425. https://pubmed.ncbi.nlm.nih.gov/34245712/, https://doi.org/10.1016/S2542-5196(21)00081-4.

Google Scholar

20.

Chen K, Breitner S, Wolf K, Stafoggia M, Sera F, Vicedo-Cabrera AM, et al. 2021. Ambient carbon monoxide and daily mortality: a global time-series study in 337 cities. Lancet Planet Health 5(4):e191–e199. https://pubmed.ncbi.nlm.nih.gov/33838734/, https://doi.org/10.1016/S2542-5196(21)00026-7.

Google Scholar

21.

WHO. 1997. International Classification of Diseases, Manual of the International Statistical Classification of Diseases, Injuries and Causes of Death. Geneva, Switzerland: World Health Organization.

Google Scholar

22.

WHO. 2016. International Statistical Classification of Diseases and Related Health Problems, 10th Revision. http://apps.who.int/classifications/icd10/browse/2016/en [accessed 25 November 2021].

Google Scholar

23.

Bhaskaran K, Gasparrini A, Hajat S, Smeeth L, Armstrong B. 2013. Time series regression studies in environmental epidemiology. Int J Epidemiol 42(4):1187–1195. https://pubmed.ncbi.nlm.nih.gov/23760528/, https://doi.org/10.1093/ije/dyt092.

Google Scholar

24.

Gasparrini A, Armstrong B, Kenward MG. 2010. Distributed lag non-linear models. Stat Med 29(21):2224–2234. https://pubmed.ncbi.nlm.nih.gov/20812303/, https://doi.org/10.1002/sim.3940.

Google Scholar

25.

Sera F, Armstrong B, Blangiardo M, Gasparrini A. 2019. An extended mixed‐effects framework for meta‐analysis. Stat Med 38(29):5429–5444. https://pubmed.ncbi.nlm.nih.gov/31647135/, https://doi.org/10.1002/sim.8362.

Google Scholar

26.

Gasparrini A, Armstrong B, Kenward MG. 2012. Multivariate meta-analysis for non-linear and other multi-parameter associations. Stat Med 31(29):3821–3839. https://pubmed.ncbi.nlm.nih.gov/22807043/, https://doi.org/10.1002/sim.5471.

Google Scholar

27.

Higgins JPT, Thompson SG. 2002. Quantifying heterogeneity in a meta-analysis. Stat Med 21(11):1539–1558. https://pubmed.ncbi.nlm.nih.gov/12111919/, https://doi.org/10.1002/sim.1186.

Google Scholar

28.

Sera F, Gasparrini A. 2022. Extended two-stage designs for environmental research. Environ Health 21(1):41. https://pubmed.ncbi.nlm.nih.gov/35436963/, https://doi.org/10.1186/s12940-022-00853-z.

Google Scholar

29.

Samoli E, Analitis A, Touloumi G, Schwartz J, Anderson HR, Sunyer J, et al. 2005. Estimating the exposure–response relationships between particulate matter and mortality within the APHEA multicity project. Environ Health Perspect 113(1):88–95. https://pubmed.ncbi.nlm.nih.gov/12626653/, https://doi.org/10.1289/ehp.7387.

PubMed

Google Scholar

30.

Chen R, Yin P, Meng X, Liu C, Wang L, Xu X, et al. 2017. Fine particulate air pollution and daily mortality. A nationwide analysis in 272 Chinese cities. Am J Respir Crit Care Med 196(1):73–81. https://pubmed.ncbi.nlm.nih.gov/28248546/, https://doi.org/10.1164/rccm.201609-1862OC.

Google Scholar

31.

Schwartz J, Spix C, Touloumi G, Bachárová L, Barumamdzadeh T, le Tertre A, et al. 1996. Methodological issues in studies of air pollution and daily counts of deaths or hospital admissions. J Epidemiol Community Health 50(suppl 1):S3–S11. https://pubmed.ncbi.nlm.nih.gov/8758217/, https://doi.org/10.1136/jech.50.suppl_1.s3.

Google Scholar

32.

Gasparrini A, Leone M. 2014. Attributable risk from distributed lag models. BMC Med Res Methodol 14(1):55. https://pubmed.ncbi.nlm.nih.gov/24758509/, https://doi.org/10.1186/1471-2288-14-55.

Google Scholar

33.

Johns DO, Svendsgaard D, Linn WS. 2010. Analysis of the concentration–respiratory response among asthmatics following controlled short-term exposures to sulfur dioxide. Inhal Toxicol 22(14):1184–1193. https://pubmed.ncbi.nlm.nih.gov/21126154/, https://doi.org/10.3109/08958378.2010.535220.

Google Scholar

34.

Zheng XY, Ding H, Jiang LN, Chen SW, Zheng JP, Qiu M, et al. 2015. Association between air pollutants and asthma emergency room visits and hospital admissions in time series studies: a systematic review and meta-analysis. PLoS One 10(9):e0138146. https://pubmed.ncbi.nlm.nih.gov/26382947/, https://doi.org/10.1371/journal.pone.0138146.

Google Scholar

35.

Chuang KJ, Chan CC, Su TC, Lee CT, Tang CS. 2007. The effect of urban air pollution on inflammation, oxidative stress, coagulation, and autonomic dysfunction in young adults. Am J Respir Crit Care Med 176(4):370–376. https://pubmed.ncbi.nlm.nih.gov/17463411/, https://doi.org/10.1164/rccm.200611-1627OC.

Google Scholar

Information & Authors

Information

Published In

Environmental Health Perspectives

Volume 131 • Issue 3 • March 2023

PubMed: 36883823

License Information

EHP is an open-access journal published with support from the National Institute of Environmental Health Sciences, National Institutes of Health. All content is public domain unless otherwise noted.

History

Received: 16 February 2022

Revision received: 22 January 2023

Accepted: 30 January 2023

Published online: 8 March 2023

Authors

Affiliations

Edward O’Brien

London School of Hygiene & Tropical Medicine, London, UK

View all articles by this author

Pierre Masselot

Department of Public Health Environments and Society, London School of Hygiene & Tropical Medicine, London, UK

Centre on Climate Change & Planetary Health, London School of Hygiene & Tropical Medicine, London, UK

View all articles by this author

Francesco Sera

Department of Statistics, Computer Science and Applications “G. Parenti,” University of Florence, Florence, Italy

View all articles by this author

Dominic Roye

Department of Geography, University of Santiago de Compostela, Santiago de Compostela, Spain

CIBER Epidemiología y Salud Pública, Madrid, Spain

View all articles by this author

Susanne Breitner

Institute for Medical Information Processing, Biometry, and Epidemiology, Ludwig Maximilian University of Munich, Munich, Germany

Institute of Epidemiology, Helmholtz Zentrum München–German Research Center for Environmental Health, Neuherberg, Germany

View all articles by this author

Chris Fook Sheng Ng

Department of Global Health Policy, Graduate School of Medicine, University of Tokyo, Tokyo, Japan

School of Tropical Medicine and Global Health, Nagasaki University, Nagasaki, Japan

View all articles by this author

Micheline de Sousa Zanotti Stagliorio Coelho

Department of Pathology, Faculty of Medicine, University of São Paulo, São Paulo, Brazil

View all articles by this author

Joana Madureira

Department of Environmental Health, Instituto Nacional de Saúde Dr Ricardo Jorge, Porto, Portugal

EPIUnit, Instituto de Saúde Publica, Universidade do Porto, Porto, Portugal

Laboratório para a Investigação Integrativa e Translacional em Saúde Populacional, Porto, Portugal

View all articles by this author

Aurelio Tobias

Department of Global Health Policy, Graduate School of Medicine, University of Tokyo, Tokyo, Japan

Institute of Environmental Assessment and Water Research, Spanish Council for Scientific Research, Barcelona, Spain

View all articles by this author

Ana Maria Vicedo-Cabrera

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

Oeschger Center for Climate Change Research, University of Bern, Bern, Switzerland

View all articles by this author

Michelle L. Bell

School of the Environment, Yale University, New Haven, Connecticut, USA

View all articles by this author

Eric Lavigne

School of Epidemiology & Public Health, Faculty of Medicine, University of Ottawa, Ottawa, Ontario, Canada

Air Health Science Division, Health Canada, Ottawa, Ontario, Canada

View all articles by this author

Haidong Kan^*

Department of Environmental Health, School of Public Health, Fudan University, Shanghai, China

View all articles by this author

Antonio Gasparrini^* https://orcid.org/0000-0002-2271-3568

Department of Public Health Environments and Society, London School of Hygiene & Tropical Medicine, London, UK

Centre on Climate Change & Planetary Health, London School of Hygiene & Tropical Medicine, London, UK

Centre for Statistical Methodology, London School of Hygiene & Tropical Medicine, London, UK

View all articles by this author

MCC Collaborative Research Network^†

View all articles by this author

Notes

Address correspondence to Antonio Gasparrini, London School of Hygiene & Tropical Medicine, 15–17 Tavistock Place, WC1H 9SH, London, UK. Telephone: +44 (0)20 7927 2406. Email: [email protected]

Metrics & Citations

Metrics

About Article Metrics

Citations

Download citation

If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click DOWNLOAD.

Cited by

Tavella R, da Silva Júnior F, Santos M, Miraglia S, Pereira Filho R, A Review of Air Pollution from Petroleum Refining and Petrochemical Industrial Complexes: Sources, Key Pollutants, Health Impacts, and Challenges, ChemEngineering, 10.3390/chemengineering9010013, 9, 1, (13), (2025).
Crossref
Psistaki K, Richardson D, Achilleos S, Roantree M, Paschalidou A, Assessing the Impact of Climatic Factors and Air Pollutants on Cardiovascular Mortality in the Eastern Mediterranean Using Machine Learning Models, Atmosphere, 10.3390/atmos16030325, 16, 3, (325), (2025).
Crossref
Masselot P, Gasparrini A, Modelling extensions for multi-location studies in environmental epidemiology, Statistical Methods in Medical Research, 10.1177/09622802241313284, (2025).
Crossref
Walton H, Dajnak D, Holland M, Evangelopoulos D, Wood D, Brand C, Assareh N, Stewart G, Beddows A, Lee S, Fecht D, Liu Y, Davies B, Goodman A, Vu T, Beevers S, Health and associated economic benefits of reduced air pollution and increased physical activity from climate change policies in the UK, Environment International, 10.1016/j.envint.2025.109283, 196, (109283), (2025).
Crossref
Chihana S, Mbale J, Chaamwe N, Unveiling the Nexus: Sulphur Dioxide Exposure, Proximity to Mining, and Respiratory Illnesses in Kankoyo: A Mixed-Methods Investigation, International Journal of Environmental Research and Public Health, 10.3390/ijerph21070850, 21, 7, (850), (2024).
Crossref
Abed Al Ahad M, Demšar U, Sullivan F, Kulu H, Long term exposure to ambient air pollution and hospital admission burden in Scotland: 16 year prospective population cohort study, BMJ Open, 10.1136/bmjopen-2024-084032, 14, 12, (e084032), (2024).
Crossref
Masselot P, Kan H, Kharol S, Bell M, Sera F, Lavigne E, Breitner S, das Neves Pereira da Silva S, Burnett R, Gasparrini A, Brook J, Air pollution mixture complexity and its effect on PM2.5-related mortality: A multicountry time-series study in 264 cities, Environmental Epidemiology, 10.1097/EE9.0000000000000342, 8, 6, (e342), (2024).
Crossref
Gasparrini A, Vicedo-Cabrera A, Tobias A, The Multi-Country Multi-City Collaborative Research Network: An international research consortium investigating environment, climate, and health, Environmental Epidemiology, 10.1097/EE9.0000000000000339, 8, 5, (e339), (2024).
Crossref
Qin L, Yang L, Liu L, Tong S, Liu Q, Li G, Zhang H, Zhu W, Liu G, Zheng M, Jiang G, Oxidative potential and persistent free radicals in dust storm particles and their associations with hospitalization, Nature Communications, 10.1038/s41467-024-55151-8, 15, 1, (2024).
Crossref
Zhu L, Yuan Y, Mayvaneh F, Sun H, Zhang Y, Hu C, Maternal ozone exposure lowers infant’s birthweight: A nationwide cohort of over 4 million livebirths in Iran, Ecotoxicology and Environmental Safety, 10.1016/j.ecoenv.2024.116840, 283, (116840), (2024).
Crossref
See more

View Options

View options

PDF

View PDF

Restore your content access

Note: This functionality works only for purchases done as a guest. If you already have an account, log in to access the content to which you are entitled.

Abstract

Background:

Objectives:

Methods:

Results:

Discussion:

Introduction

Methods

Data Collection

Statistical Methods

Main Model

Secondary Analyses

Computation of Excess Mortality

Results

Descriptive Analysis

Association between SO2 and Mortality

Exposure and Lag–Response Relationships

Analysis of Temporal Variation in Risk

Bi-Pollutant Analysis

Excess Mortality

Discussion

Conclusions

Acknowledgments

Article Notes

Supplementary Material

References

Information

Published In

License Information

History

Authors

Affiliations

Notes

Metrics

Citations

Download citation

Cited by

View options

PDF

Restore your content access

Share

Copy the content Link

Share on social media

SUPPORT AND RESOURCES

ABOUT EHP PUBLISHING

Association between $S O_{2}$ and Mortality