Eagle Hill Masthead



Northeastern Naturalist
    NENA Home
    Range and Scope
    Board of Editors
    Staff
    Editorial Workflow
    Publication Charges
    Subscriptions

Other Eagle Hill Science Journals
    Southeastern Naturalist
    Caribbean Naturalist
    Neotropical Naturalist
    Urban Naturalist
    Prairie Naturalist
    Eastern Paleontologist
    Journal of the North Atlantic
    eBio

Eagle Hill Institute Home

  Help

About Northeastern Naturalist

 

Ozone Injury to Plants within the Seney National Wildlife Refuge in Northern Michigan
Donald D. Davis

Northeastern Naturalist, Volume 14, Issue 3 (2007): 415–424

Full-text pdf (Accessible only to subscribers.To subscribe click here.)

 

Access Journal Content

Open access browsing of table of contents and abstract pages. Full text pdfs available for download for subscribers.



Current Issue: Vol. 30 (3)
NENA 30(3)

Check out NENA's latest Monograph:

Monograph 22
NENA monograph 22

All Regular Issues

Monographs

Special Issues

 

submit

 

subscribe

 

JSTOR logoClarivate logoWeb of science logoBioOne logo EbscoHOST logoProQuest logo

2007 NORTHEASTERN NATURALIST 14(3):415–424 Ozone Injury to Plants within the Seney National Wildlife Refuge in Northern Michigan Donald D. Davis* Abstract - Annual field surveys were conducted from 1999–2004 within the Seney National Wildlife Refuge in northern Michigan to determine if ambient ozone levels at this remote location were great enough to injure refuge vegetation. Ozone injury was observed on sensitive bioindicator plants during each survey year; however, the incidence (percentage) of plants exhibiting symptoms was low and varied among species and years. Ozone-induced symptoms occurred on Sambucus canadensis (American elder), Prunus serotina (black cherry), Asclepias syriaca (common milkweed), and Apocynum androsaemifolium (spreading dogbane). The most sensitive species was spreading dogbane. In addition, ozone injury was observed on a viburnum species, tentatively identified as Viburnum nudum var. cassinoides (withe-rod). Ambient ozone has been monitored since 2002 at an EPA monitoring site within the refuge. Cumulative SUM60 ozone levels (ppb-hrs) by the end of August for each survey year were greatest in 2003, followed by 2002, and least in 2004. The annual incidence of ozone injury for the 3 years was not directly related to level of ambient ozone, but was likely confounded by environmental factors such as drought. Based on the 2004 survey, the threshold level of SUM60 ozone needed to induce visible symptoms on sensitive vegetation in this remote refuge is close to 5000 ppb-hrs. Introduction The Seney National Wildlife Refuge (Seney NWR), located in the Upper Peninsula (UP) of northern Michigan (Fig. 1), was established in 1935 as a refuge and breeding ground for migratory birds and wildlife. The refuge is located near the villages of Seney and Germfask, approximately 130 km northwest of the Mackinac Bridge that joins the UP with lower Michigan. The refuge contains 38,600 ha that are managed to provide quality habitat for a diversity of waterfowl and wildlife species. In 1970, Congress conferred wilderness status on part of the Seney NWR, naming it the “Seney Wilderness,” and designating it as a Class-I air-quality area, affording the refuge further protection under the Clean Air Act as amended in 1977 (US Congress 1977). In that act, the US Congress gave federal land managers of Class-I areas responsibility to protect all air-quality related values (AQRVs) in these areas. AQRVs include vegetation, wildlife, water, soils, visibility, and cultural resources. By federal law, AQRVs in Class-I areas must be protected from deterioration. However, despite this special protection, significant levels of ambient ozone impinge on many Class-I air-quality areas in the eastern US, adversely affecting AQRVs (Lefohn and Manning 1995, *Department of Plant Pathology and Penn State Institutes of the Environment, 211 Buckhout Laboratory, The Pennsylvania State University, University Park, PA 16802; ddd2@psu.edu. 416 Northeastern Naturalist Vol. 14, No. 3 Manning et al. 1996). This paper is the third in a series dealing with ozone injury on vegetation within our national wildlife refuges (Davis 2007, Davis and Orendovici 2006). Refuge characteristics (adapted from refuge brochures) The Seney NWR, located in the Great Manistique Swamp, contains natural wetlands, grasslands, and forests. Pinus strobus L. (white pine) once dominated the forested uplands of the refuge, but these were logged during the late 1800s and the area burned many times. Controlled burning, a habitat management treatment to encourage those plant species and communities favored by fire, is practiced today within the refuge. The upland portions of the refuge are now covered with second- and third-growth stands of Populus spp. (aspen), Pinus banksiana Lamb. (jack pine), Pinus resinosa Ait. (red pine), and white pine. Wet areas are dominated by thick stands of Alnus rugosa (Du Roi) Spren. (speckled alder) and Picea mariana (Miller) BSP (black spruce), with Typha latifolia L. (common cattail) growing in the open areas. Strangmoor bogs, or “string bogs,” sub-arctic formations rare at this latitude, characterize more than half of the Wilderness Area. The string bogs are dominated by sparse stands of aspen, jack pine, and red pine, fringed with thick clumps of speckled alder, black spruce, and occasional Larix laricina (Du Roi) K. Koch (tamarack). A series of manmade dikes, resulting in 21 artificial pools and 750 islands, were built within the refuge to encourage use by waterfowl. Because of this varied habitat, the refuge is now home to more than 400 species of plants. The most predominant, ozone-sensitive bioindicator species (US DOI 2003) found within the Seney NWR is common Asclepias syriaca L. (milkweed); many thousands grow on the open, exposed, sandy dikes of the refuge. Populus species were also common in the refuge, but those growing near the survey sites were often infected with leafspots, confounding symptom evaluation. Other bioindicators in the refuge include Sambucus canadensis L. (American elder), Prunus serotina Ehrh. (black cherry), and Apocynum androsaemifolium L. (spreading dogbane). These sensitive species generally respond to ozone by producing a characteristic adaxial “stipple,” a symptom first described by Richards et al. (1958). Stipple on the leaves of ozone-sensitive plants is a reliable diagnostic symptom that can be used to evaluate the incidence of ozone injury during field surveys within wildlife refuges (Davis 2007, Davis and Orendovici 2006). Ambient ozone levels Ground-level ozone is the most important plant-damaging air pollutant in eastern North America, where elevated ozone concentrations occur annually (Comrie 1994, Coulston et al. 2003). These elevated ozone levels can travel long distances and are capable of injuring native plants in remote locations such as wilderness areas (Lefohn and Manning 1995, Manning et al. 1996) and wildlife refuges (Davis 2007, Davis and Orendovici 2006). Prior to the first survey in 1999, the 1997–1998 ozone data from the nearest EPA 2007 D.D. Davis 417 monitoring site (EPA AIRS site #55-029-0004) were examined. Data from this monitor, located approximately 120 km southwest of the refuge at Ellison Bay, WI, revealed that elevated levels of ambient ozone had occurred in both 1997 and 1998. Assuming ozone levels monitored at Ellison Bay were comparable to those occurring at the Seney NWR, located 120 km from the monitor, ozone injury was likely to occur to sensitive vegetation within the refuge. The ambient ozone levels monitored at Ellison Bay, along with the presence of ozone-sensitive bioindicators (US DOI 2003) within the refuge, provided impetus for conducting this survey. In 2002, an ozone monitor (EPA AIRS site #26-153-0001) was established within the boundaries of the Seney NWR. Comparison of 2002– 2004 ambient ozone data from the EPA site at Ellison Bay with data monitored within the Seney NWR revealed that the levels and patterns of ozone were not similar between the two monitoring sites (data not shown). Therefore, only ozone data monitored within the Seney NWR are presented in this paper. The main objective of this study was to determine if ozone injury occurred on vegetation growing in the remote Seney NWR, and if so, to determine the incidence (percentage) of the various species of bioindicator plants exhibiting stipple. A secondary objective was to determine if the incidence of ozone injury was related to ambient ozone levels and soil-moisture stress. This second objective was limited since only 3 years of monitored ozone data were available. To meet these objectives, annual surveys were conducted in the refuge during the 6-year period from 1999–2004. Methods General survey methods were similar to those used within the Brigantine NWR in New Jersey (Davis and Orendovici 2006) and Moosehorn NWR in Maine (Davis 2007). All fieldwork was conducted solely by the author. In 1998, maps of the Seney NWR were examined to select tentative survey sites that had unrestricted air movement and sunlight, criteria for suitable sampling sites (Anderson et al. 1989). Based on the maps, 25 tentative survey sites were selected and visited. Eleven sites did not contain suitable bioindicators and were eliminated, resulting in 14 final sampling sites (Fig. 1). With some modifications, the general location of these 14 sites formed the basis for the 1999–2004 field surveys. Data were not taken each year at each site, depending on the foliar condition of the bioindicator plants. The refuge was surveyed once a year during August of 1999–2004 on the dates shown in Table 1. American elder, black cherry, common milkweed, and spreading dogbane were examined at each site (Table 1). The number of plants examined within a species varied from year to year, mainly due to insect infestations that rendered foliage unusable in some years. The total number of individual plants within each bioindicator species, as well as the number of plants exhibiting stipple, were counted. Ozone-induced stipple 418 Northeastern Naturalist Vol. 14, No. 3 was simply recorded as present or absent on individual bioindicator plants; severity of injury was not evaluated. Incidence (percentage) was calculated as the number of symptomatic plants/number of plants examined for each bioindicator species. A general linear model (GLM) was used to determine if significant (p = 0.05) differences in incidence of ozone injury occurred among species across all 6 years (Minitab 2003). Since species (US DOI 2003), ozone level (Hildebrand et al. 1996), and drought stress (Showman 1991, Yuska et al. 2003) influence the level of ozone injury, we investigated the relationship between these factors using binomial logistic regression (Davis and Orendovici 2006). Annual ambient ozone levels were expressed as SUM60 (ppb-hrs), the cumulative sum of hourly ozone concentrations equaling or exceeding 60 ppb during the growing-season. Annual soilmoisture stress as of July 31 of each survey year was expressed as the Palmer Drought Severity Index (PDSI) (Palmer 1965), using data from the web at http://climvis.ncdc.noaa.gov/cgi-bin/ginterface. Predictive abilities of binomial regression analyses were severely limited since only 3 years of ambient ozone data were available. To illustrate the pattern of ambient ozone levels that impinge upon the refuge during the growing season, SUM60 ozone levels for May through September of each survey year were graphed (Fig. 2). Figure 1. Location of 14 survey sites (circles with site numbers) in the Seney National Wildlife Refuge within the Upper Peninsula of northern Michigan (map courtesy of the US Fish and Wildlife Service). 2007 D.D. Davis 419 Results and Discussion Symptom description and incidence Although not all bioindicator species were present at all survey sites, American elder, black cherry, common milkweed, and spreading dogbane were the most common bioindicators in the Seney NWR, and all exhibited classic dark, adaxial stipple as first described by Richards et al. (1958) and illustrated by Skelly (2000). These symptoms were similar to those observed on ozone-sensitive bioindicator plants in wildlife refuges in New Jersey and Maine (Davis 2007, Davis and Orendovici 2006). Chlorosis, reddening, and premature defoliation also occurred occasionally on bioindicators, but were not recorded since such symptoms can be elicited by other stresses such as heat, moisture stress, and early onset of autumnal coloration/senescence (Orendovici et al. 2003). Overall, the incidence of stipple within this remote refuge in northern Michigan was very light, except on spreading dogbane. The incidence of ozone injury on black cherry and common milkweed was much less than that reported in the more polluted areas of southwestern Michigan (Bennett et al. Table 1. Summary of observations made during the 1999–2004 surveys at the Seney National Wildlife Refuge. Numbers in table refer to number of plants exhibiting ozone-induced stipple as compared to the total number of plants evaluated for that bioindicator; data also expressed as percentages. Black American Common Spreading SUM and Year (survey date) cherry elder milkweed dogbane average 1999 (Aug 11–15) Number plants examined 49 20 7068 47 7184 Number plants injured 1 0 18 14 33 Percentage 2.04% 0.00% 0.25% 29.79% 0.46% 2000 (Aug 24–27) Number plants examined 22 35 7060 325 7442 Number plants injured 0 0 39 99 138 Percentage 0.00% 0.00% 0.55% 30.46% 1.85% 2001 (Aug 20–23) Number plants examined 62 42 5416 480 6000 Number plants injured 1 5 64 111 181 Percentage 1.61% 11.90% 1.18% 23.13% 3.02% 2002 (Aug 18–22) Number plants examined 36 20 3514 108 3678 Number plants injured 0 0 51 15 66 Percentage 0.00% 0.00% 1.45% 13.89% 1.79% 2003 (Aug 20–23) Number plants examined 90 48 2228 50 2416 Number plants injured 1 2 44 3 50 Percentage 1.11% 4.17% 1.97% 6.00% 2.07% 2004 (Aug 15–17) Number plants examined 50 21 2055 51 2177 Number plants injured 0 0 33 3 0 Percentage 0.00% 0.00% 1.61% 5.88% 0.00% Average 0.97% 3.76% 0.91% 23.09% 1.62% 420 Northeastern Naturalist Vol. 14, No. 3 2006). The overall ozone-sensitivity ranking of the bioindicators, based on mean percentage of individuals exhibiting stipple (in parentheses) across all years was: spreading dogbane (23.1%) > American elder (3.8%) > black cherry (1.0%) > common milkweed (0.9%) (Table 1). However, only spreading dogbane exhibited significantly (p = 0.05) greater incidence values, as compared to the other species, which had statistically similar incidence ratings. The incidence of ozone injury on spreading dogbane ranged from approximately 6% in 2004 to 30% in 1999 and 2000. However, these high values were mainly due to one very sensitive clone growing at one site. By 2003 and 2004, this clone was being overgrown by competing vegetation, complicating evaluation of temporal patterns of ozone injury on this species. Nevertheless, the high sensitivity of spreading dogbane also has been reported in Maine within both the Acadia National Park (Kohut et al. 2000) and the Moosehorn NWR (Davis 2007). In addition, Bergweiler and Manning (1999) reported that flowering and reproduction of spreading dogbane was inhibited by ambient ozone in Massachusetts. These observations indicate that spreading dogbane is a sensitive bioindicator and useful for detecting elevated levels of ozone. However, both Davis (2007) and Kohut et al. (2000) observed that the leaves of spreading dogbane became chlorotic and spotted, and began to senesce by late summer. These foliar disorders limit the usefulness of spreading dogbane as a bioindicator. Figure 2. Sum of hourly ozone concentrations equaling or exceeding 60 ppb (SUM60, ppb-hrs) recorded from May 1 to September 31, 2002–2004, at a monitoring location (EPA AIRS site #26-153-0001) within the Seney National Wildlife Refuge. 2007 D.D. Davis 421 American elder, black cherry, and common milkweed were also injured by ambient ozone within the Seney NRW. These species are known ozone bioindicators (US DOI 2003), but American elder and black cherry plants were found only at widely scattered locations in the refuge, making them less useful. Also, although black cherry has been reported to be sensitive to ozone (Bennett et al. 2006, Davis and Skelly 1992, Davis et al. 1981), this species exhibited a very low incidence (less than 1%) of ozone injury within the Seney refuge. This low level of ozone injury was similar to the incidence of ozone injury observed on black cherry reported for wildlife refuges in New Jersey (Davis and Orendovici 2006) and Maine (Davis 2007). The reason for low level of injury on black cherry observed during these surveys is unknown, but may be due to genetic differences in ozone sensitivity among different populations of black cherry. Common milkweed is known to be very sensitive to ozone (US DOI 2003) and occurred in very high numbers on dikes throughout the Seney NWR. More than 27,000 milkweed plants were examined during these surveys, illustrating that milkweed was the most common bioindicator in the refuge. Unlike spreading dogbane, injured milkweed plants appeared to be widely scattered among the population. One species of viburnum, tentatively identified as Viburnum nudum L. var. cassinoides (L.) Torr. & Gray (witherod), also exhibited classic stipple, but occurred infrequently. This species also exhibited classic ozone stipple at the Moosehorn NWR in Maine (Davis 2007). If stipple can be induced on withe-rod under controlled, ambient ozone dosages to prove cause and effect, withe-rod may prove useful as an ozone bioindicator in eastern North America. Ambient ozone levels The cumulative SUM60 ozone levels (ppb-hrs) monitored at the EPA site within the Seney NWR were: 22,200 (2003) > 12,685 (2002) > 5171 (2004). Ozone accumulation during the growing season followed fairly similar patterns from year to year, gradually increasing from May to the end of the summer, and becoming fairly constant in September (Fig. 2). However, the pattern was slightly different in 2002, when the ozone levels rose quickly early in the growing season. In comparison to other wildlife refuges in eastern and central US, ambient ozone levels monitored within the remote, northern Seney refuge were very low. There are few industrial or urban areas upwind from the Seney NWR, and these areas are relatively unpolluted and produce few ozone precursors. Ozone levels at Seney were less than at the Moosehorn NWR in northeastern Maine, as estimated by the nearest EPA monitor (EPA AIRS site #23-009- 0102) at Acadia National Park, ME. At Acadia, the cumulative SUM60 ozone levels ranged from approximately 26,000–33,000 ppb-hrs by late summer in 1999–2003 (Davis 2007). It is likely that the higher ozone levels at Acadia were related to ozone precursors forming in the megalopolis along the eastern seaboard southwest of Maine. In New Jersey, the cumulative SUM60 ozone levels within the Forsythe NWR (EPA AIRS site #34-001-0005), often 422 Northeastern Naturalist Vol. 14, No. 3 exceeded 40,000 ppb-hrs by the end of the summer, and actually exceeded 70,000 ppb-hrs during a high-ozone year (Davis and Orendovici 2006). These high levels were also likely influenced by ozone precursors forming in the megalopolis along the eastern seaboard. High SUM60 ozone levels, also exceeding 70,000 ppb-hrs by late summer, have been reported as well near the Mingo NWR in southeastern Missouri, as estimated from the nearest ozone monitor (EPA AIRS Site #29-186-0005). Relationship of incidence to ambient ozone and drought levels Although the incidence of ozone injury was not related to ambient ozone levels for 2002–2004, the analyses were severely limited by having only 3 years of ambient ozone data. However, Eckert et al. (1999) also reported the lack of a relationship between ambient ozone levels in Acadia National Park (Maine) and ozone injury, and attributed the lack of correlation to the confounding effects of moisture stress on stomatal functioning and resultant gas uptake. In the current study, the incidence of ozone injury on spreading dogbane, but not other species, was weakly, but significantly, related to the drought index (PDSI) according to the Pearson goodness-of-fit statistic (Minitab 2003). However, this relationship was too weak to be used for predictive purposes. Nevertheless, the threshold SUM60 ozone level needed to induce stipple on sensitive plants within the Seney refuge is likely 5000 ppb-hrs under the environmental conditions of these surveys. If this threshold level is correct, undetected ozone injury might have occurred within the refuge on plants in a sensitive stage of development as early as June in 2002 and 2003 (Fig. 2). Since most ozone injury surveys are conducted in mid- to late August, such early-season injury might no longer be apparent at survey time, due to abscission of symptomatic leaves and reduced visibility resulting from new growth. Also, high ozone levels in June may cause injury on those plant species that emerge and complete their life cycles early in the growing season, such as late-spring or early-summer ephemerals (Davis and Orendovici 2006). Ozone-induced stipple occurred on at least one bioindicator species within the refuge during each survey year, revealing that bioindicators in this remote area can be used as predictive tools as indicators of environmental deterioration. Ambient ozone levels at this remote refuge are likely great enough during most years to cause injury on refuge plants, including those growing in the Seney Class-I wilderness area. The US Fish and Wildlife Service can utilize the results of these surveys when making air-quality management decisions, including those related to review of Prevention of Significant Deterioration permits, and such data can be used to strengthen our National Ambient Air Quality Standards for ozone (US EPA 1996). Acknowledgments The author gratefully acknowledges receipt of financial support and ambient-ozone datasets from the US Fish and Wildlife Service, Air Quality Branch, Denver, CO. 2007 D.D. Davis 423 Literature Cited Anderson, R.L., C.M. Huber, R.P. Belanger, J. Knighten, T. McCartney, and B. Book. 1989. Recommended survey procedures for assessing on bioindicator plants in Region-8 Class-1 Wilderness areas. USDA Forest Service, Forest Pest Management, Asheville Field Office, Asheville, NC. Report 89-1-36. 6 pp. Bennett, J.P., E.A. Jepson, and J.A. Roth. 2006. Field responses of Prunus serotina and Asclepias syriaca to ozone around southern Lake Michigan. Environmental Pollution 142:354–366. Bergweiler, C.J., and W.J. Manning. 1999. Inhibition of flowering and reproductive success in spreading dogbane, Apocynum androsaemifolium (L), by exposure to ambient ozone. Environmental Pollution 105:333–339. Comrie, A.C. 1994. A synoptic climatology of surface ozone in rural ozone pollution at three forest sites in Pennsylvania. Atmospheric Environment 28:1601–1614. Coulston, J.W., G.C. Smith, and W.D. Smith. 2003. Regional assessment of ozone sensitive tree species using bioindicator plants. Environmental Monitoring and Assessment 83:117–127. Davis, D.D. 2007. Ozone-induced symptoms on vegetation within the Moosehorn National Wildlife Refuge in Maine. Northeastern Naturalist 14:403–414. Davis, D.D., and T. Orendovici. 2006. Incidence of ozone symptoms on vegetation within a National Wildlife Refuge in New Jersey, USA. Environmental Pollution 143:555–564. Davis, D.D., and J.M. Skelly. 1992. Foliar sensitivity of eight eastern hardwood tree species to ozone. Water, Air, and Soil Pollution 62:269–277. Davis, D.D., D.M. Umbach, and J.B. Coppolino. 1981. Susceptibility of tree and shrub species and response of black cherry foliage to ozone. Plant Disease 65:904–907. Eckert, R., R. Kohut, T. Lee, and K. Stapelfeldt. 1999. Foliar ozone injury on native vegetation at Acadia National Park: Results from a six-year (1992–1997) field survey. National Park Service , Boston, MA. 41 pp. Technical Report NPS/BSO– RNR/NRTR/00–12. Hildebrand, E., J.M. Skelly, and T.S. Fredericksen. 1996. Foliar response of ozone sensitive hardwood tree species from 1991 to 1993 in the Shenandoah National Park, VA. Canadian Journal for Forestry Research 26:658–669. Kohut, R., J. Laurence, P. King, and R. Raba. 2000. Identification of bioindicator species for ozone and assessment of the responses to ozone of native vegetation at Acadia National Park. National Park Service , Boston, MA. Technical Report NPS/BSO–RNR/NRTR/00–13. 126 pp. Lefohn, A.S., and W.J. Manning. 1995. Ozone exposures near Class-I wilderness areas in New Hampshire and Vermont. Atmospheric Environment 29:601–606. Manning, W.J., S.V. Krupa, C.J. Bergweiler, and K.I. Nelson. 1996. Ambient ozone (O3) in three Class-I wilderness areas in northeastern USA: Measurements with Ogawa passive samplers. Environmental Pollution 91:399–403. Minitab Inc. 2003. Quality Plaza, 1829 Pine Hall Road, State College, PA. Orendovici, T., J.M. Skelly, J.A. Ferdinand, J.E. Savage, M.-J. Sanz, and G.C. Smith. 2003. Response of native plants of northeastern United States and southern Spain to ozone exposures: Determining exposure/response relationships. Environmental Pollution 125:31–40 Palmer, W.C. 1965. Meteorological Drought. Research Paper Number 45. US Department of Commerce, Washington, DC. 58 pp. 424 Northeastern Naturalist Vol. 14, No. 3 Richards, B.L., J.T. Middleton, and W.B. Hewitt. 1958. Air pollution with relation to agronomic crops. V. Oxidant stipple of grape. Agronomy Journal 50:559–561. Showman, R.E. 1991. A comparison of ozone injury to vegetation during moist and drought years. Journal Air and Waste Management Association 41:63–64. Skelly, J.M. 2000. Tropospheric ozone and its importance to forests and natural plant communities of the northeastern United States. Northeastern Naturalist 7:221–236. United States Congress (US Congress). 1977. The Clean Air Act as amended August 1977. P.L. 95–95. US Government Printing Office, Washington, DC. United States Department of the Interior (US DOI). 2003. Ozone sensitive plant species on National Park Service and US Fish and Wildlife Service lands: Results of a June 24–25, 2003 workshop, Baltimore, MD. US DOI, Denver, CO. Natural Resource Report NPS/NRARD/NRR-2003/01. 21 pp. United States Environmental Protection Agency (US EPA). 1996. Air-quality criteria for ozone and related photochemical oxidants. Vol. 1 of 3 reports. EPA/600/P- 93/004aF. Research Triangle Park, NC. Yuska, D.E., J.M. Skelly, J.A. Ferdinand, R.E. Stevenson, J.E. Savage, J.D. Mulik, and A. Hines. 2003. Use of bioindicators and passive sampling devices to evaluate ambient ozone concentrations in north central Pennsylvania. Environmental Pollution 125:71–80.