Spatiotemporal Patterns of Ozone and Cardiovascular and Respiratory Disease Mortalities Due to Ozone in Shenzhen

In order to explore the temporal-spatial patterns and possible health effects of ozone in Shenzhen, daily concentrations of ozone and the daily mortality caused by cardiovascular and respiratory diseases were collected. Using Geographic Information System (GIS) and SPSS, the spatial and temporal patterns of ozone in Shenzhen were illustrated. Using a generalized additive model (GAM), the associations between ozone and cardiovascular and respiratory diseases causing mortality were analyzed, adjusted for meteorological factors and other major air pollutants including fine particulate matter (PM2.5), nitrogen dioxide (NO2) and carbon monoxide (CO), and stratified by lag, age, and gender. The results showed that, during 2013, ozone was the second main air pollutant in Shenzhen apart from PM2.5, with average daily concentrations of 95.9 μg/m3 and 76.8 μg/m3 for the ozone 1-h mean and the daily ozone 8-h maximum concentration, respectively. The daily level of ozone had a higher concentration from September to October, and relatively low concentration from May to June. Obviously, a higher concentration was found in central parts of Shenzhen with the largest population, indicating higher risks. The excess risk (ER) percentage of the cardio-respiratory mortality rate showed a clearly accumulative effect at L03, with the highest ER percentage of 1.08 (0.88–1.27) per 10 μg/m3 increase in the ozone 8-h maximum concentration for all the population. Males were found to be more sensitive to ozone compared with females, and the elderly were more susceptible to ozone exposure than younger people.


Introduction
Ozone is a secondary pollutant which is not directly emitted into the atmosphere but is created and destroyed by chemical reactions of other emitted species.Ozone is formed in the atmosphere by photo-chemical reactions in the presence of sunlight and precursor pollutants, such as the oxides of nitrogen (NO x ) and volatile organic compounds (VOCs).It is destroyed by reactions and then is deposited to the ground.Various other toxic photochemical oxidants arising from similar sources, including nitric acid and hydrogen peroxide, could also influence ozone formation [1].
Through epidemiologic and experimental studies, ground level ozone exposure has been solidly proved to be associated with total mortality, hospital admissions, respiratory symptoms, cough, chest pain, wheezing, asthma, airways inflammation, bronchitis, and cardiopulmonary disease [2][3][4][5].For the association of ozone exposure and cardiovascular disease, the evidence was less conclusive, more heterogeneous, and subject to substantial uncertainty according to a recent review [6], and the mechanisms of ozone associated with cardiovascular diseases were also less clear because ozone reacts with respiratory tissues firstly.However, there have been epidemiological studies continually identifying that ozone could enhance cardiovascular related morbidity, probably through its pro-inflammatory effects on the lung [7][8][9][10].Many factors could influence health effects caused by ozone, including the concentration and duration of the exposure [11], and the susceptibility of the subject itself.Other confounders like the meteorological conditions, other air pollutants and even individual behavior could also influence the association between health impact and ozone [12,13].
According to former research, ozone is more likely to reach unhealthy levels under hot and dry weather [14].With global warming, ozone pollution, particularly in hot coastal regions, will become riskier for public health.The association between daily mortality and ozone levels has been proved to be independent of the effects of particulate matter, but not for NO x and CO, which could accelerate the formation or decomposition of ozone and thus influence the concentration [15].Controlling these confounders in ozone related epidemiological studies is crucial to identifying its health risks.
Ground level ozone has become the secondary ambient air pollutant, following fine particulate matter (PM 2.5 ), in many southern cities in China.However, recent studies on air pollution and health have mainly focused on the primary pollutants such as nitrogen oxides (NO x ) and fine particulate matter [16][17][18][19][20][21][22]; less focus has been put on ozone research and related health risks assessment.With the increasing emphasis on primary air pollutants control (mainly NO x and SO 2 ), it is possible to experience an increase in the ozone concentration in China [23].Hence, it is necessary to explore the ozone distribution pattern and its health risks to provide scientific evidence on a strategic approach to control ozone pollution, and more importantly to facilitate health precaution strategies.
In this study, we carried out a systemic analysis on temporal-spatial patterns of ozone in Shenzhen and explored the associations between ozone and cardiovascular-respiratory disease-caused mortality, adjusted for meteorological factors and other air pollutants (PM 2.5 , NO 2 and CO), and stratified by lag days, gender and age.The objectives were to identify the spatial and temporal patterns of ozone and to elucidate the health risks caused by ozone pollution.

Study Area and Population
Shenzhen is a coastal city in the south of the Guangdong province in China, located in the Pearl River Delta.It has a subtropical oceanic climate, with warm weather and abundant rainfall.The annual average temperature is 22.4 • C, the average monthly temperature is 15.4 • C in January and 28.9 • C in July.Shenzhen has become China's most crowded city and is the fifth most densely populated city in the world, with a population density of 17,150 per square km.According to the Ministry of Environmental Protection of China, air quality is relatively good in Shenzhen and it is ranked seventh among the 74 first-stage cities implementing the National Ambient Air Quality Standard (NAAQS, GB3095-2012) in 2013.However, Shenzhen has been experiencing elevated levels of air pollution in recent years because of rapid economic development, particularly concerning ozone pollution.

Mortality Data
Cardiovascular and respiratory disease-caused mortality data from 1 January to 31 December 2013, was obtained from death certificates recorded at the Shenzhen Center for Disease Control and Prevention.The causes of death are coded according to the International Classification of Disease revision 10 (ICD-10).In this study, cardiovascular disease-caused death (CVD) (I00-I99) and respiratory disease-caused death (RD) (J00-J98) were identified.

Air Pollutants and Meteorological Data
In 2013, Shenzhen introduced the NAAQS.Air quality data was provided by the Shenzhen Environmental Monitoring Center and the China National Environmental Monitoring Center.The daily ambient air concentration of ozone was provided as the daily 1-h mean value and the daily 8-h maximum value, measured from eleven state-controlled monitoring stations in Shenzhen.Daily mean concentrations of PM 2.5 , nitrogen dioxide (NO 2 ), and carbon monoxide (CO) were also obtained from Shenzhen monitoring stations.
The location of these monitoring stations is presented in Figure 1.According to the technical guidelines of the Chinese government, the location of these monitoring stations must not be in the direct vicinity of traffic intersections or of major industrial polluters, and should also be located at a sufficient distance from any other emitting source.Thus, the monitoring data reflect the general background urban air pollution level in the urban area of Shenzhen.It is significant to note that the Nan'ao (NA) station in the east part is in the tourist area with a sparse population and the Overseas Chinese Town (Huaqiaocheng, HQC) station is in the downtown area with the densest population.

Air Pollutants and Meteorological Data
In 2013, Shenzhen introduced the NAAQS.Air quality data was provided by the Shenzhen Environmental Monitoring Center and the China National Environmental Monitoring Center.The daily ambient air concentration of ozone was provided as the daily 1-h mean value and the daily 8-h maximum value, measured from eleven state-controlled monitoring stations in Shenzhen.Daily mean concentrations of PM2.5, nitrogen dioxide (NO2), and carbon monoxide (CO) were also obtained from Shenzhen monitoring stations.
The location of these monitoring stations is presented in Figure 1.According to the technical guidelines of the Chinese government, the location of these monitoring stations must not be in the direct vicinity of traffic intersections or of major industrial polluters, and should also be located at a sufficient distance from any other emitting source.Thus, the monitoring data reflect the general background urban air pollution level in the urban area of Shenzhen.It is significant to note that the Nan'ao (NA) station in the east part is in the tourist area with a sparse population and the Overseas Chinese Town (Huaqiaocheng, HQC) station is in the downtown area with the densest population.To control for the effects of weather on mortality, meteorological data (daily mean temperature, relative humidity, air/barometric pressure and wind speed, etc.) was obtained from the Meteorological Bureau of Shenzhen Municipality.The meteorological data were measured at a fixed-site station located in the study area.This station belongs to the Meteorological Bureau of Shenzhen Municipality, the monitoring standard is consistent with international WMO (World Meteorological Organization) standard, and no data was missing.

Data Analysis
In order to quantify the association between daily levels of ozone and daily mortality, adjusted by confounders in the multivariable modeling, consistent with other time-series studies [21,22,24], we used the generalized additive model (GAM) with penalized splines to analyze the risk of cardiorespiratory mortality attributable to ozone pollution.Because the daily mortality was small and To control for the effects of weather on mortality, meteorological data (daily mean temperature, relative humidity, air/barometric pressure and wind speed, etc.) was obtained from the Meteorological Bureau of Shenzhen Municipality.The meteorological data were measured at a fixed-site station located in the study area.This station belongs to the Meteorological Bureau of Shenzhen Municipality, the monitoring standard is consistent with international WMO (World Meteorological Organization) standard, and no data was missing.

Data Analysis
In order to quantify the association between daily levels of ozone and daily mortality, adjusted by confounders in the multivariable modeling, consistent with other time-series studies [21,22,24], we used the generalized additive model (GAM) with penalized splines to analyze the risk of cardio-respiratory mortality attributable to ozone pollution.Because the daily mortality was small and typically followed a Poisson distribution, the core analysis was a GAM with log link and Poisson error that accounted for smooth fluctuations in the daily mortality.
Before conducting the model analyses, there were two steps in the procedure of the model building and model fit: development of the best base model (without a pollutant) and development of the main model (with a pollutant).The latter was achieved by adding the ozone to the final cause-specific best base model, assuming a linear relationship between the logarithmic mortality number and the ozone concentration.
First, we constructed the basic pattern of mortality number excluding ozone.We incorporated smoothed spline functions of time, weather conditions and the other air pollutants, which can include non-linear and non-monotonic links between mortality and time/weather conditions, offering a flexible modelling tool.The day of the week was also included in the basic models.
After we established the basic models, we introduced the ozone and analyzed its association with mortality.To compare the relative quality of the mortality predictions across these non-nested models, Akaike's Information Criterion (AIC) was used as a measure of how well the model fit the data.Smaller AIC values indicated the preferred model.Briefly, we fit the following log-linear generalized additive models to obtain the estimated pollution log-relative rate β in the study district: Here E(Yt) represents the expected number of mortality at day t; β represents the log-relative rate of mortality associated with a unit increase of ozone; Xi indicates the concentrations of pollutants at day t; and Wt(week) is the dummy variable for the day of the week.∑ p j=1 fj(Zj, df) is the non-parametric spline function of calendar time, temperature, barometric pressure, wind speed and humidity.A detailed introduction of the GAM has been previously described in Wood's book [25].We initialized the df as 9 df/year for time; 6 df for temperature, barometric pressure, and PM 2.5 ; 7 df for humidity, NO 2 and CO; and 8 df for wind speed.We also discussed lag effects of the air pollutants, and for the lag effects model we examined the effect of air pollutants with different lag (L) structures of single-day lag (distributed lag; from L0 to L3, L0 corresponds to the current-day pollution, and L1 refers to the previous-day concentration) and multi-day lag (moving average lag; L01 to L03, L03 corresponds to the four-day moving average of the pollutant concentration of the current and the previous three days).The meteorological factors and other confounding air pollutants used in the lag models (distributed lag model, moving average model) were the current day data.
Pearson correlation coefficients among the variables were analyzed using SPSS 22.0 (IBM Company 2013, North Castle, NY, USA).Temporal changes of ozone were summarized by Origin 9.0 (Origin Lab, Northampton, MA, USA).In the Macroscopic regional scale, the spatial distribution of the ozone concentration followed the basic assumption of "the first law of geography", namely, the regional concentrations in nearby areas are more similar than in areas farther away.Therefore, inverse distance weighted (IDW) model interpolation analysis was used to analyze spatial distributions of ozone.Spatial differences of ozone were presented by ArcGIS 10.2 (Esri, Redlands, CA, USA) using IDW.
All other statistical analyses were conducted in R3.1.0(R Foundation for Statistical Computing, Vienna, Austria) using the MGCV package (R Development Core Team, 2014).The results obtained were expressed as the excess risk (ER) percentage change in the number of mortalities caused by cardiovascular and respiratory diseases per 10 µg/m 3 increases of ozone concentrations (ER = (e δC − 1) × 100, where δC is the increased amount of air pollutants, in this study we used 10 µg/m 3 ).

Statistical Results
Table 1 summarizes the distribution of the annual mean, percentage value on the daily mortality caused by cardiovascular and respiratory diseases, concentrations of ozone 1-h mean and ozone 8-h maximum, confounding air pollutants, and meteorological factors of Shenzhen during the study period.During the study period, the mean daily temperature and humidity were 23.1 • C and 74.8%, respectively.The mean daily temperature ranged from 9.8 • C to 31.2 • C, and the mean daily humidity ranged from 24% to 100%, reflecting the subtropical oceanic climate of Shenzhen.For the confounding air pollutants, the daily NO 2 concentration ranged from 14 µg/m 3 to 111 µg/m 3 , the daily CO concentration ranged from 700 µg/m 3 to 1966 µg/m 3 , and the daily PM 2.5 concentration ranged from 20 µg/m 3 to 135 µg/m 3 with the average being 39.6 µg/m 3 .Average concentrations of ozone were 95.9 µg/m 3 for the 1-h mean and 76.8 µg/m 3 for the daily 8-h maximum ozone concentration, and the average values were higher than the median values for both the concentration of ozone 1-h mean and the daily 8-h maximum ozone concentration, indicating a larger contribution of high ozone concentration days to the annual average concentration of ozone.4646 residents died from cardiovascular diseases and 626 died from respiratory diseases in Shenzhen in 2013.The daily mortality for cardiovascular-respiratory diseases ranged from 4 to 30.

Correlation Coefficients
Meteorological conditions can alter the formation and removal of ozone and thus influence the concentration [12].Pearson correlation coefficients among ozone and meteorological factors are presented in Table 2.A significant positive correlation was found between ozone and barometric pressure, which is in line with the former research [26].Significant negative correlations were found between ozone and other meteorological factors (temperature, humidity and wind speed).
The health impact of ozone has also been proved to be influenced by meteorological factors.High temperature and ozone might have a synergistic effect on health while the opposite effect occurs for mild temperature when exposed to a high ozone concentration [14].Therefore, the health effect analysis in this study was adjusted by those meteorological confounders.Ambient NO 2 and CO could influence the formation and decomposition of ozone, and significant positive correlations between those two air pollutants and ozone are found in Table 2.Although the health effect of ozone exposure was independent from particulate matter [15], we included PM 2.5 as a confounder because it was the major air pollutant in Shenzhen and was contributing to significant air pollution-related mortalities [27].

Temporal Changes
During 2013, the average daily concentrations were 95.9 µg/m 3 ranging from 26 to 234 µg/m 3 for the ozone 1-h mean and 76.8 µg/m 3 ranging from 17 to 186 µg/m 3 for the daily 8-h maximum ozone concentration.According to the NAAQS, the Technical Regulation on Ambient Air Quality Index (on trial) (HJ633-2012) (TRAAQI), and the Technical Regulation for Ambient Air Quality Assessment (on trial) (HJ663-2013) (TRAAQA), ozone was the second main air pollutant in Shenzhen, with 30 days as a primary pollutant and 4 days as a non-attainment pollutant in the year.
Figure 2 shows the daily characteristics of the 1-h ozone concentration and the daily 8-h maximum ozone concentration in Shenzhen.The daily level of ozone had a higher value during September to October, while there was a relative lower daily concentration from May to June.The health impact of ozone has also been proved to be influenced by meteorological factors.High temperature and ozone might have a synergistic effect on health while the opposite effect occurs for mild temperature when exposed to a high ozone concentration [14].Therefore, the health effect analysis in this study was adjusted by those meteorological confounders.Ambient NO2 and CO could influence the formation and decomposition of ozone, and significant positive correlations between those two air pollutants and ozone are found in Table 2.Although the health effect of ozone exposure was independent from particulate matter [15], we included PM2.5 as a confounder because it was the major air pollutant in Shenzhen and was contributing to significant air pollution-related mortalities [27].

Temporal Changes
During 2013, the average daily concentrations were 95.9 µg/m 3 ranging from 26 to 234 µg/m 3 for the ozone 1-h mean and 76.8 µg/m 3 ranging from 17 to 186 µg/m 3 for the daily 8-h maximum ozone concentration.According to the NAAQS, the Technical Regulation on Ambient Air Quality Index (on trial) (HJ633-2012) (TRAAQI), and the Technical Regulation for Ambient Air Quality Assessment (on trial) (HJ663-2013) (TRAAQA), ozone was the second main air pollutant in Shenzhen, with 30 days as a primary pollutant and 4 days as a non-attainment pollutant in the year.
Figure 2 shows the daily characteristics of the 1-h ozone concentration and the daily 8-h maximum ozone concentration in Shenzhen.The daily level of ozone had a higher value during September to October, while there was a relative lower daily concentration from May to June.According to the daily mortality in Figure 3, males had higher rate cardiovascular and respiratory disease-caused mortality than females, and the elderly group (65+) had significantly higher mortality rate compared to the younger.For the total cardiovascular and respiratory disease-caused mortality, a slightly higher number was found between December to January, while no obvious seasonal patterns of the mortality were found among the males, the females, the younger or the elderly group.
higher mortality rate compared to the younger.For the total cardiovascular and respiratory diseasecaused mortality, a slightly higher number was found between December to January, while no obvious seasonal patterns of the mortality were found among the males, the females, the younger or the elderly group.

Spatial Differences
To represent spatial differences of ozone in Shenzhen more directly, we analyzed monthly and hourly ozone concentrations at HQC and NA sites, in the downtown and tourist areas of Shenzhen, respectively.These two parts of the city serve distinct and quite different urban functions, and differences in air quality might indicate that they were affected by different pollutant emission sources.
Figure 4 presents the spatial distribution of ozone in Shenzhen during 2013.The highest concentration could be found in Yantian, which might be attributable to the high emission from the major cargo port in Shenzhen (Yantian port).The Luohu district where Honghu and Liyuan are located is the city center with intense traffic and human activities, which could explain the high ozone concentration found there.Meanwhile, an obvious higher concentration was also found in the downtown of Shenzhen (HQC) and tourist area (NA).According to NAAQS, TRAAQI, and TRAAQA, annual concentrations of ozone at all the monitoring stations met grade II national air quality standard.

Spatial Differences
To represent spatial differences of ozone in Shenzhen more directly, we analyzed monthly and hourly ozone concentrations at HQC and NA sites, in the downtown and tourist areas of Shenzhen, respectively.These two parts of the city serve distinct and quite different urban functions, and differences in air quality might indicate that they were affected by different pollutant emission sources.
Figure 4 presents the spatial distribution of ozone in Shenzhen during 2013.The highest concentration could be found in Yantian, which might be attributable to the high emission from the major cargo port in Shenzhen (Yantian port).The Luohu district where Honghu and Liyuan are located is the city center with intense traffic and human activities, which could explain the high ozone concentration found there.Meanwhile, an obvious higher concentration was also found in the downtown of Shenzhen (HQC) and tourist area (NA).According to NAAQS, TRAAQI, and TRAAQA, annual concentrations of ozone at all the monitoring stations met grade II national air quality standard.

Temporal Differences between Different Urban Function Areas
In the temporal differences, we considered monthly differences and hourly differences of ozone for the above two different city function areas.
Figure 5a shows the monthly average concentration of ozone at HQC and NA from January to November in 2013.Both the downtown and tourist areas had the highest monthly average 1-h ozone level in January and the lowest concentration in July.The daily 8-h maximum ozone concentration showed the highest values in January and October, and the lowest level in July in the tourist area, while the downtown area had the peak value in October and lowest level in July.The ozone concentration in the tourist area was higher than that that in downtown area.
Figure 5b shows the average hourly ozone concentrations at HQC and NA in the study period.The 1-h ozone level showed the same hourly trend in the downtown and tourist areas, with the highest concentration at 16:00 and lowest level at 8:00.The daily 8-h maximum ozone concentration also presented similar hourly patterns in the two different areas, which had peak values at 19:00 and the lowest levels at 11:00.

Temporal Differences between Different Urban Function Areas
In the temporal differences, we considered monthly differences and hourly differences of ozone for the above two different city function areas.
Figure 5a shows the monthly average concentration of ozone at HQC and NA from January to November in 2013.Both the downtown and tourist areas had the highest monthly average 1-h ozone level in January and the lowest concentration in July.The daily 8-h maximum ozone concentration showed the highest values in January and October, and the lowest level in July in the tourist area, while the downtown area had the peak value in October and lowest level in July.The ozone concentration in the tourist area was higher than that that in downtown area.
Figure 5b shows the average hourly ozone concentrations at HQC and NA in the study period.The 1-h ozone level showed the same hourly trend in the downtown and tourist areas, with the

Health Effects
Table 3 presents the excess risk percentage (ER, 95% confidence interval) of the daily cardio-respiratory mortality with every 10 µg/m 3 increase in ozone concentration.To identify a possible time-delay of ozone exposure-health effects on mortality, we analyzed the lag effects of ozone.Changes on ER in cardiovascular and respiratory disease-caused mortality with a 10 µg/m 3 increase of ozone for single-day measures, 1-3 days prior to the mortality (L0-L3), and moving average measures from day 0 and day 1 to day 3 prior to mortality (L01-L03) are also listed in Table 3.During the modelling, lag effects of more than 3 days for ozone were also taken into consideration, but little relationship was found, so we only included the lag till lag 3 (L0-L3) and the cumulative lag 0-3 (L01-L03).In this study, we took the daily mortality caused by cardiovascular disease and respiratory disease as one group, and gender and age differences were also taken into consideration.Other major air pollutants and the meteorological conditions were adjusted in the model because these either can influence the health impact of ozone or can simultaneously pose significant health risks [14].
Sustainability 2017, 9, 559 9 of 13 highest concentration at 16:00 and lowest level at 8:00.The daily 8-h maximum ozone concentration also presented similar hourly patterns in the two different areas, which had peak values at 19:00 and the lowest levels at 11:00.

Health Effects
Table 3 presents the excess risk percentage (ER, 95% confidence interval) of the daily cardiorespiratory mortality with every 10 µg/m 3 increase in ozone concentration.To identify a possible time-delay of ozone exposure-health effects on mortality, we analyzed the lag effects of ozone.Changes on ER in cardiovascular and respiratory disease-caused mortality with a 10 µg/m 3 increase of ozone for single-day measures, 1-3 days prior to the mortality (L0-L3), and moving average measures from day 0 and day 1 to day 3 prior to mortality (L01-L03) are also listed in Table 3.During the modelling, lag effects of more than 3 days for ozone were also taken into consideration, but little relationship was found, so we only included the lag till lag 3 (L0-L3) and the cumulative lag 0-3 (L01-L03).In this study, we took the daily mortality caused by cardiovascular disease and respiratory disease as one group, and gender and age differences were also taken into consideration.Other major air pollutants and the meteorological conditions were adjusted in the model because these either can influence the health impact of ozone or can simultaneously pose significant health risks [14].
Generally speaking, the association between cardio-respiratory mortality and the daily 8-h maximum ozone exposure was more obvious than the ozone 1-h concentration with higher excess risk (ER) percentages; the ER percentage was higher at L2 for single day lag for both 8-h maximum  Generally speaking, the association between cardio-respiratory mortality and the daily 8-h maximum ozone exposure was more obvious than the ozone 1-h concentration with higher excess risk (ER) percentages; the ER percentage was higher at L2 for single day lag for both 8-h maximum and ozone 1-h concentration, and most of the ER percentages reached the highest points at L03, indicating the cumulative effect.In terms of the gender and age stratification, we found significant higher risks for the male group and the elderly group (65 years and older), and the association between cardio-respiratory mortality and ozone exposure was less evidenced for the females and particularly less for the younger group (0-65 years).The highest risk was found at L03 with the ER percentage of 1.08 (0.88-1.27) per 10 µg/m 3 increase in the 8-h maximum ozone concentration for cardio-respiratory mortality of all population, at L03 with the ER percentage of 1.14 (0.90-1.38) for the male group, at L2 with the ER percentage of 1.28 (1.00-1.48)for the female group, and at L03 with the ER percentage of 1.59 (1.35-1.83)for the elderly group.

Discussion
This study illustrated the spatial and temporal patterns of ozone pollution in Shenzhen, and analyzed its health risks focusing on cardiovascular and respiratory disease-caused mortality, adjusted for meteorological factors and major air pollutants.Age and gender were also stratified in the risk study to identify the susceptible group.
The ozone concentration was obviously higher in September and October, and lower from May to June.Anthropogenic activities can largely explain the highest ozone concentrations in the harbor and the city center, while a higher ozone concentration was also identified in the tourist area (NA) with less population.This was not the only case in Shenzhen.With the increasing emission of the air pollutants in the city center where the emissions were high, ozone could be decomposed together with its precursors, transferred to downwind areas and then formed into ozone again.
Meteorological factors were found to relate to the ozone concentration through its formation and decomposition [26].Barometric pressure was found to positively correlate with both O 3 1-h and O 3 8-h in this study, while, in opposition to former research, temperature was negatively associated with ozone concentration [28].This might be because the months with high temperature were accompanied with high humidity (the rainy season from May to August in Shenzhen), which can influence ultraviolet radiation and thus slow down the formation of ozone.The relationships between meteorological factors and ozone were complex and highly dependent on the range of those factors and the type of the weather clusters, and even the weather frequency change [26].With the increase of climate change, the frequency of the weather clusters and extreme weather are predicted to increase and thus the health precaution of ozone pollution could become even more challenging.
In terms of health damage, cardiovascular and respiratory mortality risk of ozone exposure was found among all the population at cumulative lags both for 1-h and daily 8-h maximum ozone concentrations, with the highest ER percentage of 1.08 (0.88-1.27).Ozone is highly oxidant and toxic.Numerous epidemiology studies have shown the association between ozone exposure and respiratory and cardiovascular diseases [2,5,29].Inhaled ozone could cause decreased lung function and increased airway hyper-reactivity, particularly for those with preexisting conditions [11,30].A recent study using the rat experiment model suggested that ozone exposure could lead to significant perturbations of coronary vascular functions possibly related to intracellular redox disturbances [31].Some other studies, using human volunteers, identified that exposure to ozone could cause adverse systemic inflammation and cardiac autonomic effects through changes in heart rate variability and C-reactive protein, and thus contribute to cardiovascular mortality [7,10].Though the mechanism is not conclusive yet, elevated ozone concentration could lead to respiratory and cardiovascular diseases according to many epidemiological studies [2][3][4][32][33][34][35]. Our study came to a similar conclusion, with identified susceptible groups being the elderly and males.
Because of largely existing pre-conditions or deprived immune systems, elderly people are more vulnerable to ozone exposure than the younger generation [36].The elderly group (65 years and older) showed the highest cardio-respiratory mortality risk with an ER percentage of 1.59 (1.35-1.83)per daily 8-h maximum ozone concentration increasing by 10 µg/m 3 in this study.From a gender perspective, animal experimental research and meta-analysis have shown that women were more susceptible to inflammatory lung disease induced by air pollution and show worse adverse pulmonary health outcomes than men [36,37].By contrast, our research found a more obvious negative association between ozone exposure and cardiovascular and respiratory mortality among males than females.This could be explained by the fact that men are more apt to be exposed to ozone because of their work outside in China.Elevated ozone exposure could affect the health of any individual, but in our research, no obvious exposure-effect correlations between ozone exposure and cardiovascular and respiratory disease-caused mortality among the younger group (0-65 years) were found.We speculated that the excess risk (ER) was estimated with less precision because of the smaller numbers of deaths in this age group compared with the elderly group.For all the sub groups, the associations of elevated ozone exposure and cardiovascular and respiratory mortality were much better proved among the cumulative lags for 8-h exposure.Susceptible individuals for ozone exposure were males and the elderly in Shenzhen.
The health risk caused by ozone exposure is influenced by the emissions of precursor chemical species, meteorology and population themselves [12,14,36,37].Precursor chemicals mainly are methane (CH 4 ) and carbon monoxide (CO), which, with the emissions of nitrogen oxides (NO x = NO + NO 2 ), contribute to a general hemispheric 'background' of ozone.Temperature and humidity are the two major meteorological factors influencing the formation and concentration of ozone, and more directly, intervening the formation of clot through the activation or suppression of the fibrinolytic pathway [12,14,38].In this research, we analyzed the association between ozone concentration and cardiovascular and respiratory disease-caused mortality confounded by those air pollutants.However, other confounders including socio-economic position and individual behavior could also influence the health effects of ozone.Because of the data limitation, we did not include those confounders.

Conclusions
In conclusion, ozone was the second main air pollutant in Shenzhen.The daily concentration of ozone was higher from September to October, while there was a relative low daily concentration from May to June.Obviously, a higher ozone concentration was found in the harbor and city center, while the tourist area, which was distant from the city center or industrial area, also had a considerably high ozone concentration, indicating the regional transmission.The ozone level showed same hourly trend in the downtown and tourist areas and the ozone 8-h maximum had peak values around 19:00 and the lowest levels around 11:00.The excess risk (ER) in the cardiovascular and respiratory disease-caused mortality would increase with the concentration level increase of ozone, particularly for the ozone 8-h maximum concentration.When the ozone concentration increased, males in the cardiovascular-respiratory disease mortality group seemed to be more sensitive than females, and people older than 65 years seemed to be affected more easily than younger people.Our findings provide additional information about the ozone distribution and health risk patterns and thus provide scientific evidence for coping strategies for health protection.

Figure 1 .
Figure 1.Distribution of 11 state-controlled air quality monitoring sites in Shenzhen.

Figure 1 .
Figure 1.Distribution of 11 state-controlled air quality monitoring sites in Shenzhen.

Figure 2 .
Figure 2. Daily concentration of ozone in Shenzhen.Figure 2. Daily concentration of ozone in Shenzhen.

Figure 2 .
Figure 2. Daily concentration of ozone in Shenzhen.Figure 2. Daily concentration of ozone in Shenzhen.

Figure 5 .
Figure 5. Average monthly (a) and hourly (b) concentrations of ozone at Overseas China Town (HQC) and Nan'ao (NA) monitoring stations.

Figure 5 .
Figure 5. Average monthly (a) and hourly (b) concentrations of ozone at Overseas China Town (HQC) and Nan'ao (NA) monitoring stations.

Table 1 .
Statistical characteristic of air pollutants, meteorological factors and daily mortality.

Table 2 .
Correlation coefficients for air pollutants and meteorological factors.correlation is significant at the 0.05 level; ** correlation is significant at the 0.01 level (two-tailed). * correlation is significant at the 0.05 level; ** correlation is significant at the 0.01 level (two-tailed). *

Table 3 .
Excess risk (ER) percentage for the daily cardiovascular and respiratory disease-caused mortality with every 10 µg/m 3 increase in ozone concentration.