Ozone source apportionment during peak summer events over southwestern Europe

Geographic sources of ozone air pollution and mortality burden in Europe

Ground-level ozone (O3) is a harmful air pollutant formed in the atmosphere by the interaction between sunlight and precursor gases. Exposure to current O3 levels in Europe is a major source of premature mortality from air pollution. However, mitigation actions have been mainly designed and implemented at the national and regional scales, lacking a comprehensive assessment of the geographic sources of O3 pollution and its associated health impacts. Here we quantify both national and imported contributions to O3 and their related mortality burden across 813 contiguous regions in 35 European countries, representing about 530 million people. Imported O3 contributed to 88.3% of all O3-attributable deaths (intercountry range 83-100%). The greatest share of imported O3 had its origins outside the study domain (that is, hemispheric sources), which was responsible for 56.7% of total O3-attributable mortality (range 42.5-87.2%). It was concluded that achieving the air-quality guidelines set out by the World Health Organization and avoiding the health impacts of O3 require not only the implementation of national or coordinated pan-European actions but also global strategies.

Dynamic evaluation of modeled ozone concentrations in Germany with four chemistry transport models

Simulating the ozone variability at regional scales using chemistry transport models (CTMs) remains a challenge. We designed a multi-model intercomparison to evaluate, for the first time, four regional CTMs on a national scale for Germany. Simulations were conducted with LOTOS-EUROS, REM-CALGRID, COSMO-MUSCAT and WRF-Chem for January 1st to December 31st, 2019, using prescribed emission information. In general, all models show good performance in the operational evaluation with average temporal correlations of MDA8 O3 in the range of 0.77-0.87 and RMSE values between 16.3 μg m-3 and 20.6 μg m-3. On average, better models' skill has been observed for rural background stations than for the urban background stations as well as for springtime compared to summertime. Our study confirms that the ensemble mean provides a better model-measurement agreement than individual models. All models capture the larger local photochemical production in summer compared to springtime and observed differences between the urban and the rural background. We introduce a new indicator to evaluate the dynamic response of ozone to temperature. During summertime a large ensemble spread in the ozone sensitivities to temperature is found with (on average) an underestimation of the ozone sensitivity to temperature, which can be linked to a systematic underestimation of mid-level ozone concentrations. During springtime we observed an ozone episode that is not covered by the models which is likely due to deficiencies in the representation of background ozone in the models. We recommend to focus on a diagnostic evaluation aimed at the model descriptions for biogenic emissions and dry deposition as a follow up and to repeat the operational and dynamic analysis for longer timeframes.

Assessing ozone abatement scenarios in the framework of the Spanish ozone mitigation plan

Tropospheric ozone (O3) is a secondary air pollutant that affects human health, vegetation and climate, especially in Mediterranean countries such as Spain. In order to tackle this long-standing issue, the Spanish government recently started to design the Spanish O3 Mitigation Plan. To support this initiative and ultimately provide recommendations, we performed a first ambitious emission and air quality modeling exercise. This study presents the development of different emission scenarios - aligned with or beyond the measures planned for 2030 in Spain - and the modeling of their respective impact on the O3 pollution across Spain (in July 2019) with both MONARCH and WRF-CMAQ air quality models. The modeling experiments include a base case scenario, a so-called planned emission (PE) scenario integrating the expected emission changes related to 2030, and a set of specific emission scenarios in which additional emission changes are applied to specific sectors (on e.g., road transport, maritime traffic) on top of the PE scenario. The planned emission scenario considerably reduces daily 8-h maximum O3 concentrations (-4 μg/m3 on average), with strongest reductions in Madrid region, north of Catalonia, Valencia region, Galicia and Andalusia. The frequency of observed daily exceedances of the 120 μg/m3 daily 8-h maximum target value and 180 μg/m3 hourly information threshold could be reduced by -37 and -77 %, respectively. The results of the specific scenarios highlight road transport and maritime traffic as two key emission sectors contributing to O3 pollution, over the entire country and the Mediterranean coast, respectively, while solvent use and industry emissions have a more limited and localized impact on O3. In any case, even with the implementation of all the emission scenarios, daily exceedances of the aforementioned thresholds will still be recorded over the country.

Summertime tropospheric ozone source apportionment study in the Madrid region (Spain)

The design of emission abatement measures to effectively reduce high ground-level ozone (O3) concentrations in urban areas is very complex. In addition to the strongly non-linear chemistry of this secondary pollutant, precursors can be released by a variety of sources in different regions, and locally produced O3 is mixed with that transported from the regional or continental scales. All of these processes depend also on the specific meteorological conditions and topography of the study area. Consequently, high-resolution comprehensive modeling tools are needed to understand the drivers of photochemical pollution and to assess the potential of local strategies to reduce adverse impacts from high tropospheric O3 levels. In this study, we apply the Integrated Source Apportionment Method (ISAM) implemented in the Community Multiscale Air Quality (CMAQ v5.3.2) model to investigate the origin of summertime O3 in the Madrid region (Spain). Consistent with previous studies, our results confirm that O3 levels are dominated by non-local contributions, representing around 70 % of mean values across the region. Nonetheless, precursors emitted by local sources, mainly road traffic, play a more important role during O3 peaks, with contributions as high as 25 ppb. The potential impact of local measures is higher under unfavorable meteorological conditions associated with regional accumulation patterns. These findings suggest that this modeling system may be used in the future to simulate the potential outcomes of specific emission abatement measures to prevent high-O3 episodes in the Madrid metropolitan area.

Source-Receptor Relationships Between Precursor Emissions and O3 and PM2.5 Air Pollution Impacts

Reduced complexity tools that provide a representation of both primarily emitted particulate matter with an aerodynamic diameter less than 2.5 μm (PM2.5), secondarily formed PM2.5, and ozone (O3) allow for a quick assessment of many iterations of pollution control scenarios. Here, a new reduced complexity tool, Pattern Constructed Air Pollution Surfaces (PCAPS), that estimates annual average PM2.5 and seasonal average maximum daily average 8 h (MDA8) O3 for any source location in the United States is described and evaluated. Typically, reduced complexity tools are not evaluated for skill in predicting change in air pollution by comparison with more sophisticated modeling systems. Here, PCAPS was compared against multiple types of emission control scenarios predicted with state-of-the-science photochemical grid models to provide confidence that the model is realistically capturing the change in air pollution due to changing emissions. PCAPS was also applied with all anthropogenic emissions sources for multiple retrospective years to predict PM2.5 chemical components for comparison against routine surface measurements. PCAPS predicted similar magnitudes and regional variations in spatial gradients of measured chemical components of PM2.5. Model performance for capturing ambient measurements was consistent with other reduced complexity tools. PCAPS also did well at capturing the magnitude and spatial features of changes predicted by photochemical transport models for multiple emissions scenarios for both O3 and PM2.5. PCAPS is a flexible tool that provides source-receptor relationships using patterns of air quality gradients from a training data set of generic modeled sources to create interpolated air pollution gradients for new locations not part of the training database. The flexibility provided for both sources and receptors makes this tool ideal for integration into larger frameworks that provide emissions changes and need estimates of air quality to inform downstream analytics, which often includes an estimate of monetized health effects.

Spatiotemporal variations of tropospheric ozone in Spain (2008-2019)

This study aims to support the development of Spain's Ozone Mitigation Plan by evaluating the present-day spatial variation (2015-2019) and trends (2008-2019) for seven ground-level ozone (O₃) metrics relevant for human/ecosystems exposure and regulatory purposes. Results indicate that the spatial variation of O₃ depends on the part of the O₃ distribution being analyzed. Metrics associated with moderate O₃ concentrations depict an increasing O₃ gradient between the northern and Mediterranean coasts due to climatic factors, while for metrics considering the upper end of the O₃ distribution, this climatic gradient tends to attenuate in favor of hotspot regions pointing to relevant local/regional O₃ formation. A classification of atmospheric regions in Spain is proposed based on their O₃ pollution patterns, to identify priority areas (or O₃ hotspots) where local/regional precursor abatement might significantly reduce O₃ during pollution episodes. The trends assessment reveals a narrowing of the O₃ distribution at the national level, with metrics influenced by lower concentrations tending to increase over time, and those reflecting the higher end of the O₃ distribution tending to decrease. While most stations show no statistically significant variations, contrasting O₃ trends are evident among the O₃ hotspots. The Madrid area exhibits the majority of upward trends across all metrics, frequently with the highest increasing rates, implying increasing O₃ associated with both chronic and episodic exposure. The Valencian Community area exhibits a mixed variation pattern, with moderate to high O₃ metrics increasing and peak metrics decreasing, while O₃ in areas downwind of Barcelona, the Guadalquivir Valley and Puertollano shows no variations. Sevilla is the only large Spanish city with generalized O₃ decreasing trends. The different O₃ trends among hotspots highlight the need for mitigation measures to be designed at a local/regional scale to be effective. This approach may offer valuable insights for other countries developing O₃ mitigation plans.

Summer O3 pollution cycle characteristics and VOCs sources in a central city of Beijing-Tianjin-Hebei area, China

One of the major pollutants influencing urban air quality in China is O₃. O₃ is the second most important pollutant affecting air quality in Shijiazhuang, which is the third largest city in the Beijing-Tianjin-Hebei area and the provincial capital of Hebei province. To fully understand the characteristics of O₃ and volatile organic compounds (VOCs), which are O₃ precursors, and the role of VOCs to ozone formation, we measured the hourly concentrations of O₃ and 85 VOCs in Shijiazhuang continuously from January to November 2020, and the concentration characteristics of both together with the chemical reactivity and sources of VOCs were analyzed from a seasonal perspective. The O₃ concentration in Shijiazhuang showed a phenomenon of high summer and low winter, and the VOCs showed a phenomenon of high winter and low spring. In the summer when the O₃ exceedance rate is the highest, the time-domain variation characteristics of O₃ were analyzed by wavelet analysis model, and the main periods controlling the O₃ concentration variation in Shijiazhuang in summer 2020 were 52 days, 32 days, 19 days and 12 days. The maximum incremental reactivity (MIR) and propylene equivalence method indicated ethene, propylene and 1-pentene were common substances in the top five species of each season. The T/B, Iso-p/N-p, Iso-p/E, N-p/E, and positive matrix factorization (PMF) model showed that industrial source (18.62%-22.03%) and vehicle emission (13.20%-17.69%) were the major VOCs sources in Shijiazhuang. Therefore, to control the O₃ concentration in Shijiazhuang, it is necessary to decrease alkenes emissions as well as VOCs from industrial source and vehicle emission.

A new NMVOC speciated inventory for a reactivity-based approach to support ozone control strategies in Spain

Ozone (O₃) pollution is a persistent problem in many regions of Spain, so understanding O₃ precursor emissions and trends is essential to design effective control strategies. We estimated the impact of Non-Methane Volatile Organic Compounds (NMVOC) species upon O₃ formation potential (OFP) using the maximum incremental reactivity approach. For this, we developed a speciated NMVOC emission inventory for Spain from 2010 to 2019 combining national reported emissions with state-of-the-art speciation profiles, which resulted in a database of emissions for over 900 individual NMVOC species and 153 individual sectors. Additionally, we analysed 2030 emission projections to quantify the expected impact of planned measures on future OFP levels. Overall, the main activities contributing to OFP in Spain are paint manufacturing and applications (20 %), manure management (16 %), and domestic solvent use (6 %). These activities contribute unevenly across regions. The more urbanised areas report a larger contribution from the solvent sector (64 % in Madrid), while in rural areas, manure management and agricultural waste burning gain importance (24 % in Extremadura), indicating that local control measures should be implemented. The top 10 NMVOC species contributing to OFP are ethanol, ethene, xylenes, propene, toluene, formaldehyde, 1,3-butadiene, styrene, n-butane, and cyclopentane, which together are responsible for 54 % of the total OFP. Our trend analysis indicates a reduction of NMVOC emissions and OFP of -5 % and -10 % between 2010 and 2019, respectively. The larger decrease in OFP is driven by a bigger reduction in xylenes (-29 %) and toluene (-28 %) from paint application industries and the road transport sector. By 2030 a significant increase (+37 %) in the OFP from the public electricity sector is expected due to the planned increase in biomass use for power generation. Our results indicate that policies should focus on paint reformulation, limiting aerosol products, and implementing NMVOC control devices in future biomass power plants.

Impact of traffic on air pollution in a mid-sized urban city during COVID-19 lockdowns

In this study, we evaluated the changes in air pollutant concentrations around Milwaukee, WI, during and after lockdown due to the COVID-19 pandemic for a period of 126 days. Measurements of particulate matter (PM1, PM2.5, and PM10), NH3, H2S, and O3 + NO2, were made on a 74-km route of arterial and highway roads from April to August 2020 using a Sniffer 4D sensor mounted to a vehicle. Traffic volume during measurement periods were estimated from smartphone-based traffic data. From lockdown (March 24, 2020-June 11, 2020) to post-lockdown (June 12, 2020-August 26, 2020) median traffic volume increased roughly 30-84%, depending upon the road type. In addition, increases in mean concentrations of NH3 (277%), PM (220-307%), and O3 + NO2 (28%) were also observed. For both traffic and air pollutants, abrupt changes in the data were observed mid-June, shortly after lockdown measures were lifted in Milwaukee County. Indeed, traffic was able to explain up to 57% of PM, 47% of NH3, and 42% of O3 + NO2 variance in pollutant concentrations on arterial and highway road segments. Two arterial roads that did not have statistically significant changes in traffic patterns during the lockdown exhibited no statistically significant trends between traffic and air quality parameters. This study demonstrated that COVID-19 lockdowns in Milwaukee, WI, caused significant decreases in traffic, which in turn had a direct impact on air pollutants. It also highlights the need for traffic volume and air quality data at relevant spatial and temporal scales for accurately assessing source apportionment of combustion-based air pollutants, which cannot be captured with typical ground-based sensor systems.

16 Years (2006-2021) of Surface Ozone Measurements in Córdoba (Southern Spain): Trends and the Impact of the COVID-19 Lockdown

Surface ozone concentrations (O₃) during the period 2006-2021 are analysed at Córdoba city (southern Iberian Peninsula) in suburban and urban sampling sites. The aims are to present the levels and temporal variations, to explore trends and to quantity the variation in O₃ concentrations in the context of the COVID-19 lockdown. The O₃ means are higher in the suburban station (62 µg m^-3 and 51.3 µg m^-3), being the information level threshold only exceeded twice during this period. The daily evolution shows a maximum at about 17:00 UTC, whereas the minimum is reached at about 9:00 UTC, with higher levels in the suburban station. The seasonal evolution of this daily cycle also presents monthly differences in shape and intensity between stations. The trends are analysed by means of daily averages and daily 5th and 95th percentiles, and they show a similar increase in all of these parameters, with special emphasis on the daily P95 concentrations, with 0.27 µg m^-3 year^-1 and 0.24 µg m^-3 year^-1. Finally, the impact of the COVID-19 lockdown shows a decline in O₃ concentrations over 10%.

Reasons for the observed tropospheric ozone weakening over south-western Europe during COVID-19: Strict lockdown versus the new normal

In this work we investigate the variation in tropospheric ozone concentrations in south-western Europe in March and April 2020 in the context of COVID-19 disease, and to what extent the former situation was recovered one year after the pandemic outbreak. To carry this study, data from 15 regional background sites in Spain, from 2010 onwards, are used. Historic (2010-2019) and most recent tropospheric ozone concentrations are compared. March and April 2020 ozone concentrations declined over 15% in most cases, rising to 23-28% at sites facing the Mediterranean. Most of the decay was related to the reduction of hemispheric background concentrations, but those sites downwind continental emissions from the Iberian Peninsula and neighbouring countries experienced an additional lessening. By exploring O₃ concentrations one year after, March and April 2021, the general decline with respect to 2010-2019 persist but its magnitude was substantially lessened with respect to the strict lockdown period. The pandemic situation unveiled that air pollution is not an endemic matter but it should be tackle with adequate actions. Ozone abatement plans for Mediterranean countries should need a pan-regional covenant in order to drop precursor emissions.

Understanding temporal and spatial changes of O3 or NO2 concentrations combining multivariate data analysis methods and air quality transport models

The application of the multivariate curve resolution method to the analysis of temporal and spatial data variability of hourly measured O₃ and NO₂ concentrations at nineteen air quality monitoring stations across Catalonia, Spain, during 2015 is shown. Data analyzed included ground-based experimental measurements and predicted concentrations by the CALIOPE air quality modelling system at three horizontal resolutions (Europe at 12 × 12 km², Iberian Peninsula at 4 × 4 km² and Catalonia at 1 × 1 km²). Results obtained in the analysis of these different data sets allowed a better understanding of O₃ and NO₂ concentration changes as a sum of a small number of different contributions related to daily sunlight radiation, seasonal dynamics, traffic emission patterns, and local station environments (urban, suburban and rural). The evaluation of O₃ and NO₂ concentrations predicted by the CALIOPE system revealed some differences among data sets at different spatial resolutions. NO₂ predictions, showed in general a better performance than O₃ predictions for the three model resolutions, specially at urban stations. Our results confirmed that the application of the trilinearity constraint during the multivariate curve resolution factor analysis decomposition of the analyzed data sets is a useful tool to facilitate the understanding of the resolved variability sources.

The effects of continentality, marine nature and the recirculation of air masses on pollen concentration: Olea in a Mediterranean coastal enclave

Olea pollen concentrations have been studied in relation to the typology of air masses, pollen grain sources and marine nature during advections in a coastal enclave in the south-eastern Iberian Peninsula. Since Spain is the world's leading olive producer, and olive growing extends throughout the Mediterranean basin, this location is ideal for the study of long-distance transport events (LTD) during the main pollen season (MPS). The air masses were classified using the calculation of 48-h back trajectories at 250, 500 and 750 m above ground level using the HYSPLIT model. After that, the frequency of LDT events from Africa and Europe was found to be 8.7% of the MPS days. In contrast, regional air masses were found in 38.6% of the MPS days. This was reflected in pollen concentrations, with significantly higher concentrations (p-value <0.05) on days with regional air masses compared to days with European air masses. Regarding the source areas, the importance of nearby sources with intense olive cultivation was confirmed (i.e., Andalusia). This proximity was relevant beyond the attenuations observed when the advections acquired a marine nature as the air mass back trajectories moved over the sea (p-value <0.001). The review of air mass typologies, source areas and pollen concentrations resulted in establishing peak dates and the detection of LDT associated with these peak dates. Distortions in the typical path of each air mass explained alterations in pollen concentrations on consecutive days. The recirculation and loops of the air mass back trajectories varied the pollen load that every type of air mass could originally contain.

Origin of regional springtime ozone episodes in the Sichuan Basin, China: Role of synoptic forcing and regional transport

The Sichuan Basin (SCB) located in southwestern China has long been considered the most polluted city cluster with exposure to unhealthy levels of ozone (O₃) at times. However, the features of O₃ regional transport and source contributions in SCB are poorly understood. In this study, ambient measurements, ERA5 reanalysis dataset, IASI O₃ column, and the Weather Research and Forecasting-Community Multiscale Air Quality (WRF-CMAQ) modeling system coupled with the Integrated Source Apportionment Method (ISAM) module were used to investigate the formation mechanism and sources of a severe O₃ episode in spring 2020 over the SCB. In the first stage of the O₃ episode, a high-pressure system persisted over the western SCB and caused northeasterly wind fields, leading to enhanced regional transport from the northern boundary with the O₃ contribution from the boundary exceeding 50% across the SCB. As the synoptic pattern evolved, southeasterly winds dominated the SCB and the stagnant zone over the Chengdu Plain confined O₃ originating from the southern SCB and Chongqing city, leading to the accumulation of precursors and elevated O₃ levels. During the O₃ episode, transportation and industrial sources were major contributors to O₃ formation especially for the Chengdu Plain and Chongqing city. In addition, the O₃-rich air mass in the nocturnal residual layer that formed over Chongqing city was transported to the Chengdu Plain through southeastern corridor at 400-1600m above ground-level under the prevailing southeasterly winds. With sunrise and the development of the atmospheric boundary layer, the O₃-rich air mass in the residual layer (RL) was entrained to the ground-level via vertical mixing, which further enhanced O₃ pollution across the Chengdu Plain. Our results revealed the mechanism of regional transport via northeastern and southeastern corridors during an O₃ episode and demonstrated the need for joint emission regulation across the SCB to mitigate O₃ pollution.

Environmental effects of stratospheric ozone depletion, UV radiation and interactions with climate change: UNEP Environmental Effects Assessment Panel, update 2019

This assessment, by the United Nations Environment Programme (UNEP) Environmental Effects Assessment Panel (EEAP), one of three Panels informing the Parties to the Montreal Protocol, provides an update, since our previous extensive assessment (Photochem. Photobiol. Sci., 2019, 18, 595-828), of recent findings of current and projected interactive environmental effects of ultraviolet (UV) radiation, stratospheric ozone, and climate change. These effects include those on human health, air quality, terrestrial and aquatic ecosystems, biogeochemical cycles, and materials used in construction and other services. The present update evaluates further evidence of the consequences of human activity on climate change that are altering the exposure of organisms and ecosystems to UV radiation. This in turn reveals the interactive effects of many climate change factors with UV radiation that have implications for the atmosphere, feedbacks, contaminant fate and transport, organismal responses, and many outdoor materials including plastics, wood, and fabrics. The universal ratification of the Montreal Protocol, signed by 197 countries, has led to the regulation and phase-out of chemicals that deplete the stratospheric ozone layer. Although this treaty has had unprecedented success in protecting the ozone layer, and hence all life on Earth from damaging UV radiation, it is also making a substantial contribution to reducing climate warming because many of the chemicals under this treaty are greenhouse gases.

Health impact assessment by the implementation of Madrid City air-quality plan in 2020

OBJECTIVES

Air pollutant concentrations in many urban areas are still above the legal and recommended limits that are set to protect the citizens' health. Madrid is one of the cities where traffic causes high NO₂ levels. In this context, Madrid City Council launched the Air Quality and Climate Change Plan for the city of Madrid (Plan A), a local strategy approved by the previous government in 2017. The aim of this study was to conduct a quantitative health impact assessment to evaluate the number of premature deaths that could potentially be prevented by the implementation of Plan A in Madrid in 2020, at both citywide and within-city level. The main purpose was to support decision-making processes in order to maximize the positive health impacts from the implementation of Plan A measures.

METHODS

The Regional Statistical Office provided information on population and daily mortality in Madrid. For exposure assessment, we estimated PM_2.5, NO₂ and O₃ concentration levels for Madrid city in 2012 (baseline air-quality scenario) and 2020 (projected air-quality scenario based on the implementation of Plan A), by means of an Eulerian chemical-transport model with a spatial resolution of 1 km × 1 km and 30 vertical levels. We used the concentration-response functions proposed by two relevant WHO projects to calculate the number of attributable annual deaths corresponding to all non-accidental causes (ICD-10: A00-R99) among all-ages and the adult population (>30 years old) for each district and for Madrid city overall. This health impact assessment was conducted dependant on health-data availability.

RESULTS

In 2020, the implementation of Plan A would imply a reduction in the Madrid citywide annual mean PM_2.5 concentration of 0.6 μg/m³ and 4.0 μg/m³ for NO₂. In contrast, an increase of 1 μg/m³ for O₃ would be expected. The annual number of all-cause deaths from long-term exposure (95% CI) that could be postponed in the adult population by the expected air-pollutant concentration reduction was 88 (57-117) for PM_2.5 and 519 (295-750) for NO₂; short-term exposure accounted for 20 (7-32) for PM_2.5 and 79 (47-111) for NO₂ in the total population. According to the spatial distribution of air pollutants, the highest mortality change estimations were for the city centre - including Madrid Central and mainly within the M-30 ring road -, as compared to peripheral districts. The positive health impacts from the reductions in PM_2.5 and NO₂ far exceeded the adverse mortality effects expected from the increase in O₃.

CONCLUSIONS

Effective implementation of Plan A measures in Madrid city would bring about an appreciable decline in traffic-related air-pollutant concentrations and, in turn, would lead to significant health-related benefits.

Characteristics of Ozone Pollution, Regional Distribution and Causes during 2014–2018 in Shandong Province, East China

The summer ozone pollution of Shandong province has become a severe problem in the period 2014–2018. Affected by the monsoon climate, the monthly average ozone concentrations in most areas were unimodal, with peaks in June, whereas in coastal areas the concentrations were bimodal, with the highest peak in May and the second highest peak in September. Using the empirical orthogonal function method, three main spatial distribution patterns were found. The most important pattern proved the influences of solar radiation, temperature, and industrial structure on ozone. Spatial clustering analysis of the ozone concentration showed Shandong divided into five units, including Peninsula Coastal area (PC), Lunan inland area (LN), Western Bohai area (WB), Luxi plain area (LX), and Luzhong mountain area (LZ). Influenced by air temperature and local circulation, coastal cities had lower daytime and higher nighttime ozone concentrations than inland. Correlation analysis suggested that ozone concentrations were significantly positively correlated with solar radiation. The VOCs from industries or other sources (e.g., traffic emission, petroleum processing, and chemical industries) had high positive correlations with ozone concentrations, whereas NOx emissions had significantly negatively correlation. This study provides a comprehensive understanding of ozone pollution and theoretical reference for regional management of ozone pollution in Shandong province.