The goal of this study is to evaluate the effects of anthropogenic climate change on air quality, in particular on ozone, during the summer in the U.S. mid-Atlantic region. First, we establish a connection between high-ozone (HO) days, defined as those with observed 8-h average ozone concentration greater than 70 parts per billion (ppb), and certain weather patterns, called synoptic types. We identify four summer synoptic types that most often are associated with HO days based on a 30-yr historical period (1986–2015) using NCEP–NCAR reanalysis. Second, we define thresholds for mean near-surface temperature and precipitation that characterize HO days during the four HO synoptic types. Next, we look at climate projections from five models from phase 5 of the Coupled Model Intercomparison Project (CMIP5) for the early and late midcentury (2025–34 and 2045–54) and analyze the frequency of HO days. We find a general increasing trend, weaker in the early midcentury and stronger in the late midcentury, with 2 and 5 extra HO days per year, respectively, from 16 in 2015. These 5 extra days are the result of two processes. On one hand, the four HO synoptic types will increase in frequency, which explains about 1.5–2 extra HO days. The remaining 3–3.5 extra days are explained by the increase in near-surface temperatures during the HO synoptic types. Future air quality regulations, which have been successful in the historical period at reducing ozone concentrations in the mid-Atlantic, may need to become stricter to compensate for the underlying increasing trends from global warming.
The most recent report of the Intergovernmental Panel for Climate Change (IPCC) indicates that an average global warming of about 2°C is likely by the end of the century, based on medium- to high-emission scenarios (IPCC 2013). Along the U.S. East Coast, the warming is projected to be even stronger with an increase of 2.5°–5.5°C (Maloney et al. 2014). Heat waves are expected to be more frequent and more intense and to last longer. For example, the number of days above 90°F (32°C) may increase by 50 days yr−1 in the mid-Atlantic region by the middle of the century if CO2 emissions continue on their current trajectory (Melillo et al. 2014). Such increases in extreme temperatures are likely to have a negative impact on public health all along the mid-Atlantic, as high temperatures are associated with increased mortality. This temperature–mortality link appears to be strongest in the eastern United States compared to other parts of the country (Curriero et al. 2002; Anderson and Bell 2009).
Global warming may exacerbate indirect health impacts via near-surface ozone (O3) pollution (Monks et al. 2015). Ozone is a summer pollutant, regulated by the U.S. Environmental Protection Agency (EPA), that forms in the atmosphere as a result of photochemical reactions between nitrogen oxides (NOx) and volatile organic compounds. Ozone has large negative impacts on health, especially affecting the cardiopulmonary and respiratory systems (EPA 2018b). In Delaware, ozone is the only regulated pollutant that exceeds national and state standards regularly [Delaware Department of Natural Resources and Environmental Control, Division of Air Quality (DNREC-DAQ) 2016], as will be discussed in detail later (section 2a). We focus on Delaware as representative of the mid-Atlantic region because of its central location. Further, since Delaware is a small state, even a few stations can represent ozone conditions over the state well.
Ozone is expected to reach higher concentrations (i.e., increase in intensity) and/or exceed national and state standards more often (i.e., increase in frequency) because of climate change via several mechanisms. First, higher temperatures enhance photochemical reactions that form ozone. Second, high-ozone (HO) episodes tend to occur during certain synoptic conditions, such as a high pressure system (e.g., the Bermuda high). These synoptic conditions may occur more often in the future, leading to more frequent HO events (Shen et al. 2015; Leibensperger et al. 2008). This paper will focus on this second indirect mechanism, that is, changes in the frequency of HO events in the future caused by changes in synoptic conditions. Since Delaware is a state with low local emissions and its ozone pollution is dominated by transport from upwind states (Moghani et al. 2018), it is a region where ozone tends to be dominated by synoptic, more than local, weather conditions.
How can we make projections of whether or not HO events may occur more frequently in the future? The most common tool used to create projections of future climate is coupled atmosphere–ocean models, such as the models included in phase 5 of the Coupled Model Intercomparison Project (CMIP5; Taylor et al. 2012). The results of these model simulations form the basis of the papers that are used in the IPCC Assessments, for example, IPCC (2013). However, most models (unless specifically outfitted with complex chemistry models) do not include ozone as a prognostic variable. Therefore, to make projections of how HO days may change in the future, we must relate ozone concentrations to variables that all coupled atmosphere–ocean climate models simulate, like pressure and temperature. Then we can determine how those variables are changing and the implications of their changes on the frequency of HO days.
Here we make use of CMIP5 simulations to project how the number of HO days may change over the mid-Atlantic region (and more specifically, Delaware), for two midcentury periods, 2025–34 and 2045–54, and compare this to the frequency of ozone days for the present day (1986–2015). Section 2 gives an overview of the observed and model data sources, while section 3 details the synoptic typing and other criteria used to classify days as “high ozone.” Section 4 discusses the results, while section 5 offers conclusions and directions for future research.
2. Data sources
a. Delaware ozone data
The period of interest for this study is 1986–2015. The DNREC DAQ operates a network of air quality stations throughout the state, recording a number of parameters, including ozone. The first three ozone monitoring sites (Claymont, Lums Pond, and Bellefonte) became operational in 1982, with additional stations being added, removed, or replaced since then. Since 2011, there are seven stations that record ozone with 1-h frequency (DNREC-DAQ 2018). The locations of these seven ozone monitoring stations are depicted in Fig. 1 and listed in Table 1. Since Delaware is a small state composed of only three counties, these seven stations provide a reasonable spatial coverage, with two stations in southernmost Sussex County, one station in central Kent County, and four stations in the more heavily populated northernmost New Castle County. These data are archived by DNREC and made available via the U.S. EPA Air Data portal (EPA 2018a).
Even though hourly values are recorded, the daily maximum 8-h average ozone concentration is used in this analysis to represent the daily value of ozone at each station. These 8-h averages have been error checked by DNREC and EPA. For example, an 8-h average is considered valid if at least six of the 8-hourly concentrations for the 8-h period are available and a daily maximum 8-h average ozone concentration is considered valid if valid 8-h averages are available for at least 13 of the 17 consecutive 8-h periods starting from 0700 to 2300 LT (EPA 2015). The daily maximum 8-h average ozone value is used to calculate the so-called design value of each station, which is the 3-yr average of the annual fourth-highest daily maximum 8-h average ozone concentration at each station. The National Ambient Air Quality Standard (NAAQS) for ozone is met when the design value is less than or equal to 0.070 ppm and it is violated when the design value is greater than 0.070 ppm (EPA 2015). In Delaware, three of the seven ozone monitoring stations were in violation of the NAAQS in 2018 (Table 1).
b. NCEP–NCAR reanalyses
Reanalysis data used for the synoptic typing come from the NCEP–NCAR reanalysis, version 1 (Kalnay et al. 1996). Reanalysis datasets are produced by assimilating various quality-controlled observations, including rawinsondes, surface marine data, aircraft data, surface land synoptic data, satellite sounder data, special sensing microwave imager, and satellite cloud drift winds, into a numerical weather prediction model (Parker 2016). The result is a gridded product that provides four-times daily “snapshots” of the atmosphere. The horizontal resolution of the NCEP–NCAR reanalysis products is about 2.5° latitude × 2.5° longitude. Here we use the variables sea level pressure (SLP), surface air temperature, 500-hPa height, and precipitation rate. Data are available for the period from 1948 to the present, but again we focus on the period of 1986–2015.
c. CMIP5 climate model data
The CMIP5 archive contains model output from a plethora of models worldwide, run for dozens of experiments, as discussed in section 1. Here we focus on the historical experiment simulations for the period 1986–2005 and the RCP8.5 scenario for 2006–15 (1986–2015 will be referred to as the “historical” period hereafter). Several different scenarios based on the radiative forcing by the end of the twenty-first century were used to generate projections of future climate in the CMIP5 project. The scenarios range from a low of 2.6 W m−2 to a high of 8.5 W m−2. The RCP8.5 scenario explores a world in which population increases to 12 billion by 2100 along with income inequality. Meanwhile, limited technological advances mean that there is no concomitant increases in energy efficiency (Riahi et al. 2011). This scenario thus represents a “business as usual” trajectory for the future (van Vuuren et al. 2011; Taylor et al. 2012) that is realistic given the failure of countries to so far meet the goals of the Paris Agreement (Masson-Delmotte et al. 2018). Regardless of the choice of scenario, the sign of most changes in the future projections is consistent between the different scenarios, meaning if temperature increases in one scenario with more radiative forcing, it increases in another simulation with less radiative forcing, but to a lesser extent (Giorgi and Mearns 2002). Two future RCP8.5 10-yr periods are considered: the early (2025–34) and late (2045–54) midcentury. Both periods are close enough in time to be directly affected by policy decisions made today. In the early midcentury period, however, the underlying climate is not too different from today’s, while the late midcentury period is likely to show stronger effects from global warming.
We used output from five different coupled atmosphere–ocean global circulation models in the CMIP5 archive (listed in Table 2) based on their performance over eastern North America, horizontal grid spacing, and availability of the relevant fields at daily and subdaily time scales. The horizontal grid spacing of the five models (degrees latitude by degrees longitude) is reported in Table 2. Pierce et al. (2009) have shown that using five models is sufficient to represent the behavior of the ensemble as a whole. Daily gridded data for the following variables were used for the analysis: SLP, surface air temperature, and precipitation. Because of Delaware’s small size, all or most of the state is contained in a single grid cell of the CMIP5 models used, and thus the values for this single cell are used when local data are needed for temperature and precipitation. Additionally, the full gridded data for SLP over the eastern United States is used for spatial correlations with the NCEP reanalyses.
a. High-ozone days from observations
The NAAQS for ozone would be the most obvious choice to assess observed high-ozone days, but this metric has not been constant over the years. For example, the 8-h average standard for ozone was changed from 0.080 ppm in 1997 (EPA 1997) to 0.075 ppm in 2008 (EPA 2008) and then further reduced to 0.070 ppm in 2015 (EPA 2015). In addition, the NAAQS for ozone is based on a rather complicated formula, including the calculation of the fourth-highest percentile each year and its average over three consecutive years (EPA 2015). Therefore, following the Air Quality Index (AQI; EPA 2017), a day is classified here as an HO day when the 8-h average ozone concentration at any station in Delaware exceeds the threshold of 70 ppb. In contrast to a large state, which would require multiple stations to identify observed high-ozone events, even just a single station is reasonable for a small state like Delaware.
Based on these criteria, we counted the frequency of observed HO events during the 30-yr observational record of ozone from 1986 to 2015. As shown in Fig. 2, the number of observed HO days decreased over the 30-yr period. For example, during the 1980s the average number of HO days was approximately 75, whereas by 2015 it had decreased to less than 20. On average, each year during the study period had about 2 fewer days of high ozone per year (as the trend line in Fig. 2 has a slope of −1.82 HO days yr−1). This steady decline resulted from stricter regulations (such as federal and state air quality standards) and decreases in the emissions of ozone-forming pollutants.
b. Synoptic typing
Episodes of high ozone are associated with certain synoptic weather conditions, such as the passage of midlatitude cyclones and the western expansion of the Bermuda high (Shen et al. 2015). Because of climate change, the synoptic weather conditions conducive to high ozone concentrations may occur with greater frequency, therefore increasing the potential for more HO episodes.
The synoptic weather conditions present on any given day can be classified by the use of synoptic typing. Synoptic typing uses an assortment of meteorological variables, including temperature, humidity, sea level pressure, precipitation, cloud cover, and wind velocity, to investigate common trends in the synoptic conditions, in order to identify a limited number of types that explain the underlying synoptic-scale variability. The origins of these techniques can be traced to visual classifications by Lamb (1972) and computer-assisted statistical relationships of sea level pressure by Lund (1963).
The synoptic typing used in this study was originally developed in Siegert et al. (2017), using an eigenvector-based approach (Kalkstein and Corrigan 1986; Yarnal 1993). Siegert et al. (2017) employed synoptic typing to assess hydroclimate in the mid-Atlantic states, with the same dataset expanded and used by Suriano and Leathers (2017) in evaluating lake-effect snowfall in the Great Lakes region. Siegert et al. (2017) describe the procedure used to create the synoptic types as follows. Weather data from the period from 1948 to 2015 (temperature, dewpoint temperature, u and υ wind components, atmospheric pressure, and cloud cover) taken four times daily (0900, 1500, 2300, 0300 UTC) at the Philadelphia International Airport Station were input into an unrotated R-mode principal components analysis. They then identified the five principal components explaining at least 75% of the variance in the mid-Atlantic region of the United States (i.e., Pennsylvania, New Jersey, Maryland, and Delaware) for winter, spring, summer, and autumn. Within-group average linkage clustering analysis was used on the seasonal component loadings to produce clusters, or synoptic types. Initially 20 clusters were retained and composite maps were developed for each cluster representing SLP, 500-hPa geopotential height, surface air temperature, and surface precipitation rate using the NCEP–NCAR reanalysis (Kalnay et al. 1996). Because some of the clusters showed similar characteristics in the composite maps, they were combined based on expert assessment. This resulted in a variable number of synoptic types for each season: 12 in winter, 13 in spring, 9 in summer, and 11 in autumn. This research focuses on the nine synoptic types from the summer identified in Siegert et al. (2017), listed in Table 3, since the vast majority of high-ozone events in Delaware take place during the summer months.
Each summer day of the 30-yr record evaluated here (1986–2015) was assigned to one of the nine summer synoptic types, as defined by (Siegert et al. 2017), and to an ozone category (HO or non-HO), as described in section 3a. In this way, the frequency distribution of synoptic types associated with HO days could be determined (Fig. 3). The nine types and their frequencies are shown in yellow in Fig. 3; the distribution of HO days by synoptic type is shown in blue.
The four most common synoptic types for the entire 30-yr period are southwesterly flow (SW), New England low (NELOW), overhead high (OHHI), and weak pattern (WEAK); the two synoptic types with the highest number of HO days are SW and NELOW (Table 3). Several synoptic types have spatial patterns that are highly correlated with each other, as shown in Table 4. The spatial pattern correlation between two SLP maps A and B representing either two synoptic types or an individual day’s map and a synoptic type was calculated as follows:
where and represent the mean SLP for maps A and B, respectively; and represent the SLP at grid point i in maps A and B, respectively; and n represents the number of grid points. For example, both types HBLOW and WEAKSW are highly correlated (above 0.93, Table 4) with type SW, since all share a generally southwesterly flow associated with the presence of the Bermuda high located around 30°N, 50°W, which advects high temperatures and ozone precursors over the mid-Atlantic (Fig. 4). In contrast, the NELOW type is characterized by northwesterly flow that would advect precursors from the upper Midwest. The near-surface air temperature patterns between all four of these synoptic types are extremely similar (Fig. 5), with mean temperatures exceeding 25°C over the mid-Atlantic, except for the NELOW pattern, which has slightly lower temperatures. Additionally, many of the OHHI HO days occurred in the first decade (50 days in 1986–95, not shown), and both the frequency of HO days that occurred during OHHI and the overall number of days of the OHHI synoptic types trended downward during the two subsequent decades (from 50 days in 1986–95 to 34 days in 1996–2005 to 15 days in 2006–15), reducing the influence of this type in more recent years. Therefore, we classified Hudson Bay low (HBLOW) and weak southwesterly flow (WEAKSW) as high-ozone types as well. Thus, our four HO synoptic types are: NELOW, SW, HBLOW, and WEAKSW.
c. High-ozone days from the CMIP5 climate models
Once the synoptic types associated with observed HO days have been determined, this information is used to classify which climate model days (historical and future simulations) have the potential to be HO days.
First, for each model day t we calculated its spatial correlation with each of the nine summer synoptic types using SLP [Eq. (1)]. The day t was assigned the type j that had the highest correlation, as long as it was greater than a minimum threshold of 0.5. Furthermore, if the model day was classified as NELOW, SW, HBLOW, or WEAKSW, then it was a possible HO day.
However, classifying a climate model day as one of the four HO synoptic types is necessary but not sufficient for a HO day to occur. Based on our analysis of observed HO days in observations and reanalysis data, there are two additional conditions that must be met for a high-ozone day to occur.
First, high temperatures are required within those four HO synoptic types. We determined that 90% of the HO events that occurred during synoptic type SW (the warmest type) had a daily mean temperature of 25.5°C or higher. Although the other three HO synoptic types had lower temperature thresholds, we selected 25.5°C in order to capture the more “extreme” cases of the other three types and to minimize the probability of false HO days. Note that this threshold is for the mean daily temperature; the daily highest temperature on HO days consistently exceeded 32°C (90°F), but mean temperature is used since that variable is more readily available in the CMIP5 archive. Also, the choice of 25.5°C allowed us to better capture the observed distribution of HO days over the four synoptic types during the last decade (2006–15). Since ultimately we want to compare future climate to the present, it is more reasonable to calibrate our criteria based on the last decade of the historical period.
Second, only days when there was no measurable precipitation in the models (defined as less than 1 mm day−1) were considered to be favorable to high ozone. Precipitation prevents ozone formation because of both reduction of sunlight, which drives the photochemical reactions, and washout of the available precursor pollutants.
In summary, any model day that met all three criteria (HO synoptic type, mean daily temperature above 25.5°C, and precipitation below 1 mm day−1) was classified as an HO day.
This section describes the results of the application of the three criteria described previously in section 3c to the five selected CMIP5 models, for the 30-yr historical record (1986–2015) and two future decades: 2025–34 and 2045–54. Before presenting the results, we would like to point out a few important properties of climate simulations that will help interpreting the findings of this study.
First, we do not expect an HO day in the model to exactly match an HO day in the real world. Weather evolves in a climate model the same way as it does in the real world, but not on a day-to-day or year-to-year basis. For example, if temperature was above 35°C on 1 September 2015, it is not expected to be above 35°C on the same day or year in the climate model. Rather, the climate model should have the same frequency of days above 35°C as the real world over the course of a decade. With respect to ozone, climate models are not expected to perfectly capture the number of HO days year by year, but rather the characteristics of the underlying variability and frequency over time scales of a decade or longer. In other words, the model results show how the number of HO days would have changed over the past 30 years if the regulatory and societal factors were held constant at 2015 levels (i.e., a standard of 70 ppb for 8-h average ozone, along with 2015 emission standards and industry pollution controls). This allows the climate models to provide estimates of how the number of HO days would fluctuate solely because of the climate, and not to societal and regulatory changes.
Second, since climate models do not account for changes in regulations or emissions that impact air pollution but do account for the positive correlation between ozone and temperature, the historical climate simulations will show increasing ozone days because of warming (as discussed shortly in section 4a), while observations previously indicated decreasing HO days because of air quality protection legislation and emission regulations (section 3a). Nonetheless, there is still value in analyzing trends in climate models over the historical period because 1) they explain how HO frequency is impacted by climate change alone and 2) they provide a baseline for comparison against future projections of HO days.
a. Trends in high-ozone days from the CMIP5 climate models
The selected CMIP5 models were first evaluated for the 30-yr historical period from 1986 to 2015. This allows for assessment of the performance of the detection algorithm for a period of time for which ozone observation data were readily available. Recall that we calibrated our ozone day detection algorithm using the period 2006–15 (section 3c). Using this decade, we see that the average number of ozone days is about 16 yr−1 by 2015 (Fig. 6a), which slightly underestimates the observed number of days by three according to the trend line in Fig. 2. The most noticeable feature in Fig. 6a is that the number of HO days simulated by the CMIP5 models increases, which may appear to contradict the observations shown in Fig. 2. However, regulations introduced during the 30-yr period, accompanied by advances in emissions reducing technologies (e.g., improved automobile fuel economy), resulted in a decline in HO days. The effects of these changes are not included in climate models. Therefore, the models show increasing HO days because of the higher temperatures simulated in response to elevated greenhouse gases. Effectively, the models show an approximation of the number of HO days that would have been expected earlier in the time period if emissions standards were consistent over the full 30-yr period. The models that simulate the highest number of HO days are the HadGEM2 and MIROC5 models, with between 20 and 23 HO days (Fig. 6b). All the models show a positive trend over 1986–2015, although the values vary from +0.04 to +0.37 HO days yr−1. The slope of the model median over the entire 30-yr period is +0.225 HO days yr−1, which means that it takes about 5 years to add another HO day per year. Since the criteria were tuned to the most recent decade (2006–15), the model results match the observations very well during those 10 years, but underestimate the number of HO days in the two earlier decades.
Next we look at two future decades (2025–34 and 2045–54) in Figs. 6. The slope of the median during the early midcentury is similar to the historical period (+0.230 HO days yr−1). Therefore, it takes about four years to increase the number of HO days by one. However, the trend accelerates during the late midcentury, increasing to +0.582 HO days yr−1, corresponding to an additional ozone day every 2 years. The variability does not change between the three periods, as the coefficient of variation is approximately the same (0.53, 0.52, 0.51, respectively, for 2006–15, 2025–34, and 2045–54).
All models show an increase in the number of HO days in the two midcentury periods compared to the historical period (Fig. 6b). There is a fairly wide spread in the number of HO days projected by the models with a range of 10–35 HO days by 2055. The Meteorological Research Institute (MRI) model projects the lowest number of HO days in all three periods, although notably the number of HO days increases in the later two periods. MIROC5 shows consistently high numbers of HO days in all three periods, with 25 HO days in the historical period and 40 HO days by the end of 2055. CCSM4 is similar to the model median for all periods. Some models show variability between the two midcentury periods, with more HO days during 2025–35 than during 2045–55. This is not unexpected because of the use of two fairly short periods in the midcentury, since the different representative concentration pathways do not diverge until the midcentury (IPCC 2013) and the radiative forcing remains below 5 W m−2.
b. Climatic factors driving high-ozone-day changes
The increases in HO days projected by the models may be the result of three factors: an increase in the frequency of synoptic types conducive to HO days (i.e., the four HO synoptic types SW, NELOW, HBLOW, and WEAKSW); an increase in daily mean temperature; and a decrease in the number of days with precipitation.
Starting with the synoptic types first, Fig. 7 shows the median predictions of number of days in each of the four HO synoptic types (without using the temperature or precipitation criteria). For two of the types (WEAKSW and NELOW), the frequency does not change substantially between the three periods, although it never decreases from the historical to the late midcentury period. However, there is a noticeable increase in the number of SW synoptic type days (by +1.5 days yr−1, from 7 to 8.5 days yr−1) in 2045–54 relative to the historical period. Also synoptic type HBLOW is projected to occur with higher frequency in 2045–54 than in the historical period (by +1 day yr−1, from 10 to 11 days yr−1), although less significantly that the SW type.
The increase in the number of SW and HBLOW synoptic type events (+1.5 and +1 days yr−1, respectively, for a total of +2.5 days yr−1) is smaller than the overall increase in HO days (+5 days yr−1, Fig. 6). Because HO days in the model simulations must be one of these four synoptic types (i.e., changes in other synoptic types are not considered), increases in temperature or increases in the frequency of dry days (or decreases in the frequency of days with precipitation) must be contributing to the increase in HO days. To evaluate this, we examined what changes the models showed over the summer months in the average daily mean temperature and in the number of dry days predicted. Overall, the models did not show any statistically significant trends in the number of dry days (not shown). However, the models did predict an increase in the average daily mean temperature (Fig. 8); the overall model mean increased from 24°C in 1986 to nearly 27°C in 2054, a change in mean summer temperature of +0.38°C decade−1. By the decade 2045–54, in fact, the five-model daily mean temperature exceeded the 25.5°C threshold that we adopted for HO days every year. Furthermore, it can be seen in Fig. 8 that the model with the lowest mean summer temperature, MRI, also predicted the fewest HO days (Fig. 6b), while one of the warmest climate models, MIROC5, also predicted the most HO days. This confirms qualitatively the importance of underlying temperature changes, not only synoptic type frequency changes, to the projected increases in the number of HO days by midcentury.
5. Discussion and conclusions
In this study, we examined climate model projections for the mid-twenty-first century to determine how the number of days with ozone levels high enough to be unhealthy for sensitive populations—HO days—may change over the mid-Atlantic. We used a two-step approach. First we identified the atmospheric circulation patterns (“synoptic types”) that were most likely to occur during observed HO days in the past. Four synoptic patterns are associated with high-ozone days over the mid-Atlantic: New England low (NELOW), southwesterly flow (SW), weak southwesterly flow (WEAKSW), and Hudson Bay low (HBLOW). Second, we imposed two other criteria: no precipitation and high daily mean temperatures (exceeding 25.5°C).
Next we applied the three criteria above to five selected CMIP5 models (HadGEM2, MIROC5, CCSM4, CNRM-CM5, and MRI) to calculate the frequency of HO days in the future. Before looking at the RCP8.5 midcentury projections, however, we prepared model predictions for the historical period of 1986–2015, to establish a baseline for comparison. Midcentury projections from the five CMIP5 models show that HO days are expected to increase, from about one extra HO day every 5 years to 1 extra HO day every 2 years. This is due to an increase in the frequency of HO synoptic types (in particular, the dominant SW type associated with southwesterly flow) and to increases in temperature, without noticeable changes in the frequency of dry days.
Historically, the number of high-ozone days in Delaware decreased by 2 days every year because of stricter air quality regulations, despite the underlying increases in temperature due to increasing greenhouse gases. Our results indicate that this progress will be partially negated as temperatures continue to increase and southwesterly-flow synoptic types become more frequent. Therefore, maintaining current ozone levels will require even stricter regulations than currently exist. Our study focused on the frequency of HO days. However, the ozone concentration during HO days itself may also increase (i.e., more intense HO days).
Future work should focus on quantifying these potential changes and their direct health impacts. In addition, other synoptic types, not currently associated with HO days, may become HO types in the future; some synoptic types may disappear; and entirely new HO synoptic types may manifest themselves. Thus we may be underestimating the changes in future HO days. Future studies should also examine spring and fall, as these transition seasons are expected to become warmer in the future, and look for potential changes in the duration of HO events.
This research was sponsored by the Delaware Department of Natural Resources and Environmental Control (DNREC), Division of Energy and Climate (DEC) and Division of Air Quality (DAQ), Grant 16A01486. In particular, the authors would like to acknowledge the support of Ali Mirzakhalili, former Director of DAQ at DNREC. The authors would also like to thank Daniel Leathers and Zachary Suriano for providing synoptic typing data. We thank Emily Aiken, who assisted with downloading CMIP5 data. NCEP reanalysis data were provided by the NOAA/OAR/ESRL PSD, Boulder, Colorado, and are available from their website at http://www.esrl.noaa.gov/psd/. This research was conducted in part on the Mills and Farber high-performance computer clusters of the University of Delaware.
Current affiliation: Department of Marine and Coastal Sciences, Rutgers, The State University of New Jersey, New Brunswick, New Jersey.