Severe Surface Ozone Pollution in China: A Global Perspective

Jul 26, 2018 - Cooperative Institute for Research in Environmental Sciences ... We also find an increase in the surface ozone level over China in 2016...
1 downloads 0 Views 1MB Size
Subscriber access provided by Kaohsiung Medical University

Characterization of Natural and Affected Environments

Severe surface ozone pollution in China: a global perspective Xiao Lu, Jiayun Hong, Lin Zhang, Owen R. Cooper, Martin Schultz, Xiaobin Xu, Tao Wang, Meng Gao, Yuanhong Zhao, and Yuanhang Zhang Environ. Sci. Technol. Lett., Just Accepted Manuscript • DOI: 10.1021/acs.estlett.8b00366 • Publication Date (Web): 26 Jul 2018 Downloaded from http://pubs.acs.org on July 26, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 20

Environmental Science & Technology Letters

1

Severe surface ozone pollution in China: a global perspective

2 3

Xiao Lu1, Jiayun Hong1, Lin Zhang1*, Owen R. Cooper2,3, Martin G. Schultz4,

4

Xiaobin Xu5, Tao Wang6, Meng Gao7, Yuanhong Zhao1, Yuanhang Zhang8*

5

(1) Laboratory for Climate and Ocean-Atmosphere Studies, Department of

6

Atmospheric and Oceanic Sciences, School of Physics, Peking University, Beijing

7

100871, China

8

(2) Cooperative Institute for Research in Environmental Sciences (CIRES), University

9

of Colorado, Boulder, CO, 80309, US

10

(3) NOAA Earth System Research Laboratory, Boulder, CO, 80305, US

11

(4) Jülich Supercomputing Centre, Forschungszentrum Jülich, 52425, Germany

12

(5) State Key Laboratory of Severe Weather & Key Laboratory for Atmospheric

13

Chemistry of China

14

Meteorological Sciences, Beijing, 100081, China

15

(6) Department of Civil and Environmental Engineering, The Hong Kong Polytechnic

16

University, Hong Kong, 99907, China

17

(7) School of Engineering and Applied Sciences, Harvard University, Cambridge, MA,

18

02138, USA

19

(8) State Key Joint Laboratory of Environmental Simulation and Pollution Control,

20

College of Environmental Sciences and Engineering, Peking University, Beijing

21

100871, China

22

Correspondence to: Lin Zhang ([email protected]) and Yuanhang Zhang

23

([email protected])

Meteorological Administration, Chinese Academy of

24 25

1

ACS Paragon Plus Environment

Environmental Science & Technology Letters

26

Abstract

27

The nationwide extent of surface ozone pollution in China and its comparison to the

28

global ozone distribution have not been recognized due to the sparsity of Chinese

29

monitoring sites before 2012. Here we address this issue by using the latest 5-year

30

(2013-2017) surface ozone measurements from the Chinese monitoring network,

31

combined with the recent Tropospheric Ozone Assessment Report (TOAR) database

32

for other industrialized regions such as Japan, South Korea, Europe, and the US

33

(JKEU). We use various human health and vegetation exposure metrics. We find that

34

although the median ozone values are comparable between Chinese and JKEU cities,

35

the magnitude and frequency of high ozone events are much larger in China. The

36

national warm-season (April-September) 4th highest daily maximum 8-h average

37

(4MDA8) ozone (86.0 ppb) and number of days with MDA8 greater than 70 ppb

38

(NDGT70, 29.7 days) in China are, respectively, 6.3-30% (range of regional mean

39

differences) and 93-575% higher than the JKEU regional averages. Health exposure

40

metrics such as warm-season mean MDA8 and annual SOMO35 (Sum of Ozone

41

Means Over 35 ppb) are 6.3-16% and 25-95% higher in China. We also find an

42

increase of surface ozone over China in 2016-2017 relative to 2013-2014. Our results

43

show that on the regional scale human and vegetation exposure to ozone is greater in

44

China than other developed regions of the world with comprehensive ozone

45

monitoring.

46 47

TOC Graphic

48 2

ACS Paragon Plus Environment

Page 2 of 20

Page 3 of 20

Environmental Science & Technology Letters

49 50

1. Introduction

51

Surface ozone is an air pollutant detrimental to human health and vegetation growth1.

52

There is evidence of risks of respiratory and cardiovascular mortality due to

53

short-term (acute) exposure to high ambient ozone, and recent studies also reveal

54

long-term (chronic) exposure effects even at low ozone levels2-4. Ozone near the

55

surface is mainly generated by photochemical oxidation of carbon monoxide (CO)

56

and volatile organic compounds (VOCs) in the presence of nitrogen oxides

57

(NOx=NO+NO2) and sunlight. Severe ozone pollution in many European and United

58

States (US) urban areas has now been extensively alleviated owing to stringent

59

emission control measures since the 1990s5-8. In the meantime, South and East Asia

60

have experienced rapid urbanization and industrialization, causing significant

61

increases in anthropogenic ozone precursor emissions9,10 and potentially shifting the

62

worldwide air pollution hotspots to these populous regions11.

63 64

Eastern China experiences severe fine particulate matter (PM2.5) pollution in winter12,

65

and this issue has been the main focus of the government’s air pollution control

66

strategy over the past 5 years13. More recently, summertime ozone pollution has

67

become an emerging concern in China. The past summer of 2017 recorded

68

particularly high surface ozone in a large number of Chinese cities, with the 90th

69

percentile of the daily maximum 8-h average (MDA8) ozone in 30 of the 74 major

70

cities exceeding 200 ug m-3 (the Grade II national air quality standard for protection

71

of residential areas)14. Previous studies from field observations and satellite retrievals

72

have indicated severe ozone pollution in China, e.g. observed hourly maximum ozone

73

concentrations in China frequently exceeded 150 ppb15,16. These values are

74

comparable to or higher than the present-day annual 4th highest MDA8 (4MDA8,

75

approximately 98th percentile) ozone values in the Los Angeles Basin (about 110

76

ppb)8, an area with the most severe ozone pollution in the US. In contrast to the

77

general decreasing ozone in the US and Europe5, available surface ozone observations

78

have shown significant positive trends in China since 199017-21. Despite the increasing 3

ACS Paragon Plus Environment

Environmental Science & Technology Letters

Page 4 of 20

79

attention, the severity of ozone pollution in China compared with other industrialized

80

regions of the world has not been quantified based on extensive ozone monitoring.

81 82

This study aims to understand the current status of observed surface ozone pollution

83

in China from a global perspective by comparison to ozone observations in other

84

industrialized regions of the Northern Hemisphere, as revealed by ozone metrics

85

relevant to human health and crop/ecosystem productivity. The data come from

86

China’s recently available nationwide surface ozone measurements and from the

87

Tropospheric

88

http://www.igacproject.org/activities/TOAR) organized by the International Global

89

Atmospheric Chemistry Project (IGAC)22. We also analyze changes of ozone

90

pollution in China over the past 5 years based on the different ozone exposure

91

metrics.

Ozone

Assessment

Report

(TOAR,

92 93

2. Materials and Methods

94

The China National Environmental Monitoring Center (CNEMC) network.

95

Surface air pollutants in mainland China are monitored by the CNEMC of the

96

Ministry of Environmental Protection in China (MEPC) and data are reported hourly

97

(http://106.37.208.233:20035/). This nationwide observational network, designed for

98

monitoring urban and suburban air pollution in mainland China, became operational

99

in 2013 at 74 major cities, and by 2017 it includes 1597 non-rural sites covering 454

100

cities (Figure S1). These measurements now document the air quality in Chinese

101

cities and have been applied in recent studies23-26. In this study we use the hourly

102

measurements of ozone, carbon monoxide (CO), nitrogen dioxide (NO2), sulfur

103

dioxide (SO2), and fine particulate matter (PM2.5) from 2013 to 2017. Quality controls

104

for all species have been applied following previous studies23 to remove data outliers

105

(Supporting Information).

106 107

TOAR surface ozone database. We use the TOAR surface ozone dataset (accessed

108

from https://doi.org/10.1594/PANGAEA.876108)27 for a global comparison of 4

ACS Paragon Plus Environment

Page 5 of 20

Environmental Science & Technology Letters

109

ground-level ozone. The TOAR surface ozone database includes ozone statistics and

110

metrics relevant to studies on climate, human health and vegetation exposure,

111

calculated from hourly ozone measurements at more than 9,000 monitoring sites

112

around the world since the 1970s22. Procedures for data harmonization, quality control,

113

and metric calculation are described by Schultz et al. (2017)22. Here we analyze the

114

latest 5-year (2010-2014) aggregated ozone metrics as well as the long-term

115

(1980-2014) yearly ozone metrics from the TOAR database. Recent reports suggest

116

that levels of ozone or its precursors in the US and Europe have remained relatively

117

stable or slightly decreased after 201428,29. We use the warm-season (April-September)

118

metrics, except for SOMO35 (Sum of Ozone Means Over 35 ppb) that is aggregated

119

annually. While the TOAR database includes thousands of observational sites in

120

Europe, US, Japan, and Korea, it only has data from 32 sites located in China (half are

121

in Hong Kong). Our analysis of the CNEMC observations (not included in the TOAR

122

surface ozone database) thus fills this gap.

123 124

Ozone metrics for human health and crop/ecosystem exposure. Table 1 lists the 11

125

ozone metrics analyzed in this study. These metrics include standard statistics such as

126

median, 98th-percentile (Perc98), MDA8, and daytime average (DTAvg), as well as

127

metrics for assessing human health and ecosystem exposure impacts30. MDA8 is a

128

widely used daily ozone metric for air quality regulation and human health impact

129

studies in many regions such as the US, Europe, and China. The response of human

130

health to ozone exposure is controlled by ozone level, exposure duration frequency,

131

and physical activity4,30. Following TOAR, ozone exposure is evaluated in our study

132

by short-term human exposure metrics such as daily maximum 1-hour (MDA1) and

133

MDA8, and total number of days with MDA8 >70 ppb (NDGT70), or by metrics that

134

consider the cumulative exposure (e.g. annual SOMO35). Risk estimates from

135

epidemiological studies based on these exposure metrics are summarized in ref (4)

136

and have been recently applied in China31. Ozone damage to vegetation is cumulative

137

and can be disproportionally larger for high ozone concentrations30. Here we use the

138

AOT40 and W126 metrics which correlate with crop yields32,33 as indicators of 5

ACS Paragon Plus Environment

Environmental Science & Technology Letters

139

vegetation exposure. We calculated all metrics based on hourly CNEMC

140

measurements following the TOAR definitions, procedures, and data completeness

141

requirements (Table S1). We also calculated the annual total numbers of days with the

142

ozone level exceeding the Chinese Grade II (Exceedance) ozone air quality standard

143

(Table 1).

144 145

3. Results and discussion

146

Spatiotemporal distribution of ozone air quality across mainland China. Figure 1

147

shows the spatial distribution of annual mean MDA8 (annual AVGMDA8) in Chinese

148

cities averaged for 2013-2017. MDA8 and MDA1 are also used to determine

149

exceedances of the national ozone air quality standard (Table 1) in China. The annual

150

AVGMDA8 and AVGMDA1 ozone values averaged for the Chinese sites are 41.2±6.3

151

ppb (87.9±13.5 ug m-3) and 48.2±7.2 ppb (103.6±15.4 ug m-3), respectively. During

152

2013-2017 there were over 27 days per year that exceeded the Grade II air quality

153

standard, averaged over the monitoring sites. Spatially, ozone hotspots (annual

154

AVGMDA8 > 50 ppb and Exceedance > 40 days) extend across eastern China

155

especially in the North China Plain (NCP) and the Yangtze River Delta (YRD),

156

mainly induced by high anthropogenic sources of ozone precursors. Severe ozone

157

pollution also occurs in some western cities, e.g., in the central Gansu Province where

158

Exceedance II is greater than 80 days, comparable to the eastern regions. This is likely

159

due to unique topography (valley basins in mountainous regions) coupled with high

160

local ozone precursor emissions from the petrochemical industry and vehicles34.

161 162

The mean MDA8 ozone over China peaks in summer due to strong photochemistry,

163

similar to other regions at northern mid-latitudes, but the patterns vary among

164

different regions (Figure S2 and S3). MDA8 ozone levels in the YRD and NCP peak

165

in May (60 ppb for the regional average) and June (68 ppb), respectively, while the

166

averaged MDA8 value in the Pearl River Delta (PRD) is highest in October (57 ppb).

167

This difference is due to the timing of the Asian summer monsoon, which brings

168

cloudy and rainy weather, marine air, and strong deep convection, all unfavorable for 6

ACS Paragon Plus Environment

Page 6 of 20

Page 7 of 20

Environmental Science & Technology Letters

169

ozone chemical production and accumulation35, 36. As a consequence, ozone pollution

170

over YRD and PRD generally decreases during the summer monsoon, and peaks

171

before its arrival or after its retreat.

172 173

Comparison with other northern mid-latitude regions. Figures 2 and S4 compare

174

different human health and vegetation exposure metrics for warm-season surface

175

ozone in China, to those in other industrialized regions including Japan and South

176

Korea (JK), Europe (EU), and the US (hereafter referred as JKEU). We use the

177

non-rural sites (Supporting Information) from the TOAR database. We show that

178

warm-season Median and mean DTAvg ozone metric values averaged for China

179

(30.2±8.1 ppb and 41.6±8.6 ppb, respectively) are comparable to other regions, e.g.,

180

Median values of 29.0 ppb over JK, 31.0 ppb over EU, and 33.5 ppb over the US

181

(Figure S4). Active local anthropogenic emissions and photochemistry lead to high

182

DTAvg ozone in East Asia and California37, 38. High Median and DTAvg values also

183

occur at high elevations such as western China, the European Alps, and the western

184

US, reflecting the typical increase of ozone with altitude39-42.

185 186

The severity of Chinese surface ozone pollution becomes distinct when comparing

187

metrics such as 4MDA8 and Perc98, which focus on the high end of the ozone

188

distribution and are likely caused by local pollution episodes. As shown in Figure 2

189

and Figure S4, 4MDA8 and Perc98 ozone values averaged over the Chinese sites are

190

86.0±14.4 ppb and 80.7±14.1 ppb, respectively, approximately 20-25% higher than

191

the averages of European and US sites. High 4MDA8 ozone values above 100 ppb are

192

widely observed in the NCP, YRD and PRD (Figure 2). High ozone days also occur

193

much more frequently in China, as indicated by the NDGT70 metric (Figure 2). The

194

NDGT70 value averages 29.7±22.0 days over the Chinese sites, approximately twice

195

the value in JK, a factor of 6 higher than EU, and a factor of 3 higher than the US.

196

NDGT70 values above 70 days are common in eastern China, while in other regions

197

only a few sites in California reach such a high frequency of exceedance.

198 7

ACS Paragon Plus Environment

Environmental Science & Technology Letters

199

In addition to NDGT70, we compare two other health exposure metrics (AVGMDA8

200

in Figure 2 and SOMO35 in Figure S4). The national warm-season mean AVGMDA8

201

value in China (50 ppb) is 6.3%-16% higher than those found in JKEU (43.1-47.0

202

ppb), while the mean annual SOMO35 value in China (4.3±1.7 ×103 ppb day) is

203

25-95% higher than the regional mean SOMO35 values of 2.2-3.4 ×103 ppb day in

204

JKEU cities. The AVGMDA8 and SOMO35 values are particularly high in populous

205

eastern China (Table S2), e.g. AVGMDA8 > 55 ppb in NCP and SOMO35 > 5.5 ×103

206

ppb day in YRD. To our knowledge, no other region of the world, with extensive

207

ozone monitoring, has such a large population exposed to such high and frequent

208

ozone pollution episodes.

209 210

Comparisons of vegetation exposure metrics also indicate the potential for greater

211

ozone-induced plant damage in China. The AOT40 (Figure 2) and W126 (Figure S4)

212

metrics in China are 2.2±1.1×104 and 2.9±1.7×104 ppb hour, which are 35-100% and

213

50-170% higher than the JKEU regional means, respectively. W126 shows a larger

214

difference due to its weighting towards higher ozone levels. High AOT40 (>4×104

215

ppb hour) and W126 (>6×104 ppb hour) values are widespread in eastern and central

216

China. Previous field surveys at limited locations, and coarse-resolution modeling

217

studies (e.g., 2 degrees) estimated that ozone pollution in China could decrease the

218

wheat yield by 6.4-14.9% and the annual net primary production by about 14%33, 43.

219

These studies focused on ozone pollution before 2010, and thus we may expect even

220

larger effects for the more recent conditions recorded by the nationwide network.

221 222

Changes in surface ozone pollution over China in 2013-2017. The 5-year CNEMC

223

measurements allow us to examine the interannual variability of surface ozone

224

pollution over China. Here we focus on the sites in 74 major Chinese cities (Figure S1)

225

where continuous observations from 2013 to 2017 are available.

226 227

Figure 3 shows the evolution of different ozone metrics for the Chinese sites over

228

2013-2017 compared with those in Japan, Europe, and the US over 1980-2014. Figure 8

ACS Paragon Plus Environment

Page 8 of 20

Page 9 of 20

Environmental Science & Technology Letters

229

S5 shows interannual variations of additional ozone metrics averaged over these

230

Chinese sites during 2013-2017. We find that all ozone metrics show continuous

231

increases over the 5 years in China. Averaged over the 74 cities, the DTAvg, 4MDA8,

232

and Perc98 values increased at rates of 3.7-6.2% year-1, while the ozone exposure

233

metrics SOMO35, AOT40, and W126 increased at higher rates (11.9-15.3% year-1).

234

The increases are statistically significant in 2016-2017 relative to 2013-2014 based on

235

analysis of variance tests (Table S4), and are consistent with previous studies of

236

long-term ozone increases at a limited number of sites across China17-21. Even though

237

the 5-year period is too short to derive statistically robust trend estimates, these results

238

indicate an increasing severity of human and crops/ecosystems ozone exposure across

239

China. We find the largest increases of ozone exposure in eastern and central China,

240

especially in the NCP and YRD, overlapping with the most populous areas (Figure

241

S6). From an historical perspective, present-day ozone in the major Chinese cities is

242

comparable to or even higher than the 1980 levels in the US when emission controls

243

were just beginning to have an impact on reducing ozone (Figure 3).

244 245

The Chinese State Council implemented the Action Plan on Air Pollution Prevention

246

and Control in September 2013, resulting in nationwide reductions of ambient SO2,

247

CO, NO2, and PM2.5 levels, as shown in Figure S7. The emerging severity of ozone

248

pollution in China raises a new challenge to current emission control actions. Previous

249

observational and modeling analyses showed that ozone chemical production in the

250

NCP, YRD, and PRD regions are likely VOC-sensitive or mixed-sensitive16,44-49.

251

Recent bottom-up emission estimates and satellite formaldehyde observations

252

indicated increasing VOCs levels in eastern China that can be attributed to

253

anthropogenic sources 20,50. For those VOC-sensitive regions, either decreasing NOx

254

levels or increasing VOCs levels could potentially enhance ozone pollution. Therefore

255

VOC controls could be explored as a possible control strategy. Reduced PM2.5 levels

256

over 2013-2017 may also cause an increase of ozone due to impacts on ozone

257

photochemistry and heterogeneous chemistry on aerosol surfaces51. In addition, the

258

interannual variability of surface ozone can be influenced by meteorological 9

ACS Paragon Plus Environment

Environmental Science & Technology Letters

259

conditions. The summer of 2017 had hotter and drier weather conditions compared to

260

previous years (Figure S8), which favored ozone production and led to higher ozone

261

levels. Further studies are needed to better quantify the contributions of emission

262

changes and weather variability on recent trends in surface ozone over China.

263 264

We conclude that China has become a global hotspot of present-day surface ozone

265

pollution, and human and vegetation exposure in China is greater than in JKEU.

266

Although statistics such as Median and mean DTAvg, which focus on the mid-range

267

of the ozone distribution, are comparable to JKEU, the magnitude and frequency of

268

high ozone events are much greater in China. While high ozone pollution also occurs

269

in some other developing regions, such as India in the pre-summer monsoon season

270

(April-May)36 and Mexico City52, the lack of accessible ozone observations prevents

271

us from comparing those regions to China and JKEU. China’s surface ozone air

272

quality deteriorated in 2016-2017 compared to 2013-2014, indicating that current

273

emission control measures have not been effective for reducing ozone air pollution.

274 275

Supporting Information

276

The CNEMC data quality control (SI 1); station type in the TOAR database (SI 2);

277

calculation of the ozone metrics (Table S1), ozone metric values over NCP, YRD, and

278

PRD (Table S2); ozone metric values over China, Japan/South Korea, Europe, and the

279

US (Table S3);ANOVA test of year on year ozone increase (Table S4); site locations

280

(Figure S1); Annual cycle of MDA8 over China (Figure S2); seasonal and spatial

281

distribution of MDA8 over China (Figure S3); additional ozone metrics over China,

282

Europe, and the US (Figure S4); interannual variability of ozone metrics (Figure S5

283

and S6); interannual variations of national mean NO2, CO, SO2, and PM2.5 (Figure S7);

284

surface temperature and relative humidity anomalies in summer 2017 (Figure S8).

285 286

Notes

287

The authors declare no competing financial interest.

288 10

ACS Paragon Plus Environment

Page 10 of 20

Page 11 of 20

Environmental Science & Technology Letters

289

Acknowledgements

290

This work was funded by the National Key Research and Development Program of

291

China

292

(2014CB441303), the National Natural Science Foundation of China (41475112), and

293

the Clean Air Research Project in China (Grant No. 201509001). We acknowledge the

294

Tropospheric Ozone Assessment Report (TOAR) initiative for providing the surface

295

ozone data used in this study.

(2017YFC0210102),

China’s

National

Basic

296 297

11

ACS Paragon Plus Environment

Research

Program

Environmental Science & Technology Letters

Page 12 of 20

298

Table 1. Description of ozone metrics used in this study. More details of the

299

calculation methods are provided in Table S1. Metrics [unit]

Definition

Aggregation period

50th-percentile of hourly concentrations. 98th-percentile of hourly concentrations. daytime average ozone is the average of hourly ozone concentrations for the 12-h period from 08:00 to 19:59 local time. MDA1 [ppb] daily maximum 1-hour average. AVGMDA1 represents mean MDA1 in the aggregation period. MDA8 [ppb] daily maximum 8-hour average. AVGMDA8 represents mean MDA8 in the aggregation period. 4MDA8 [ppb] 4th highest MDA8. SOMO35 [ppb the sum of positive differences between MDA8 day] and a cut-off concentration of 35 ppb. NDGT70 total number of days with MDA8 > 70 ppb. [day] AOT40 accumulative hourly ozone concentrations above [ppb hour] 40 ppb. Daily 12-h AOT40 is calculated using hourly values for the 12-h period from 08:00 to 19:59 local time, and warm-season total AOT40 is presented. W126 Daily W126 is calculated using the formula: [ppb hour] W126 = ∑ w × C , where C denotes hourly ozone concentration in ppb for the 12-h period from 08:00 to 19:59 local time., w is the weighting index defined as Median [ppb] Perc98 [ppb] DTAvg [ppb]



w = ∙(∙

April-September April-September April-September

annual and April-September annual and April-September April-September annual; April-September April-September

April-September

, M=4403, A = 126.

 / )

Exceedance [day]

Number of days with ozone concentration annual exceeding the Chinese Grade II national air quality standard, defined as MDA8 >160 µg m-3 or MDA1 > 200 µg m-3 over residential, industrial, and rural areas.

300 301

12

ACS Paragon Plus Environment

Page 13 of 20

Environmental Science & Technology Letters

302 303 304 305 306

Figure 1. Distribution of annual mean MDA8 (left) and ozone air quality exceedance (right) averaged over 2013-2017. Only sites with at least 3-year measurements are included. Mean values and standard deviations over the sites are shown inset.

13

ACS Paragon Plus Environment

Environmental Science & Technology Letters

307 308 309 310 311 312 313 314 315 316

Figure 2. Comparison of April-September ozone metrics (4MDA8, NDGT70, AVGMDA8, AOT40) between China (including CNEMC sites and TOAR Chinese sites) and other industrialized regions (Japan and Korea (JK), Europe (EU), and the US. Only sites with at least 3-year measurements are included. Values inset are regional mean and standard deviation averaged over the N sites. The differences between China and other industrialized regions are statistically significant (p