EPA/600/A-97/090 Ecological Effects of Tropospheric Ozone: A U.S. Perspective - Past, Present and Future W.E. Hogsett and Christian P. Andersen U.S. EPA National Health and Environmental Effects Research Laboratory Western Ecology Division Corvallis, Oregon ABSTRACT Understanding the effects of tropospheric ozone on vegetation, as called for in the U.S. Clean Air Act, has involved collection of experimental data at the species level and, in particular, at the level of the individual and populations. Frequently the studies have been regression designs involving single species and single pollutants resulting in quantitative • exposure-response functions that characterize the effects on biomass or reproduction (crop yield). Recently, the 1996 EPA Oxidant Criteria Document reviewed the published research on oxidant effects on crops, forests, and ecological resources, and concluded that the current secondary National Ambient Air Quality Standard (NAAQS) for ozone is neither protective or appropriate. The subsequent discussion and decision on what is an appropriate form and level of the secondary NAAQS utilized results from single pollutant and single species studies. To estimate concentrations causing effects composite response-functions for crops and trees from a large number of species were developed, predicting crop yield loss or annual biomass loss in seedlings. This approach is useful based on available quantitative data on biological effects; however it assumes that individual plant response does not change in the presence of other stresses or in natural systems which are considerably more complex (e.g. species' assemblages and competition for resources). For example, the importance of biological complexity is illustrated when we examine the below-ground ecosystem, an often overlooked portion of the ecosystem. We have found that incorporating natural biological complexity into potting soils can result in carbon fluxes opposite those predicted from individual plant studies using artificial media lacking natural soil foodwebs. A future research approach to understanding 03 effects on ecosystems is required that will develop the necessary linkages to extrapolate experimental data taken at the individual level, often in artificial conditions, to predict changes on individuals or populations in more complex native environments. We will present experimental and modelling activities from our laboratory that' show how we are beginning to address the problems of scale, complexity, and multiple stresses in forested ecosystems exposed to ozone stress. Questions frequently posed in ------- ecological risk assessments. I. INTRODUCTION The Clean Air Act mandates the protection of ecological resources from adverse effects of air pollutants through setting a secondary National Ambient Air Quality Standard (NAAQS). Ecological resources are part of public welfare and include crops, forests, soils, ecosystems. As part of the standard setting process, EPA critically reviews all the pertinent literature every 5 years, and presents the conclusions in the Ozone and other Photochemical Qxidants Criteria Document (U.S. EPA, 1996). This document is scientifically peer-reviewed and reviewed by the Clean Air Science and Advisory Committee (CASAC). The revised and approved document then provides the scientific basis for EPA's decisions regarding the NAAQS for ozone. The most recent Oxidant Criteria Document was completed in 1996 (U.S. EPA, 1996), and concluded that the general conclusions of the previous 1986 Criteria Document and the 1992 Supplement were not altered. The general conclusions were (a) current ambient ozone concentrations in many areas of the country were sufficient to impair growth and yield of plants; (b) effects occur with only a few hourly occurrences above 0.08 ppm; (c) several species exhibited growth and yield effects when the mean ozone concentration exceeded 0.05 ppm for 4-6 hr/day for at least 2 weeks; (d) regression analyses of the National Crop Loss Assessment Network (NCLAN) data developing exposure-response functions for yield loss indicated that at least 50% of the crops were predicted to exhibit a 10% yield loss at 7-hr seasonal mean ozone concentrations of 0.05 ppm or less; (e) European crop yield loss studies substantiated the effects observed in this country; and (f) studies of forest tree seedlings substantiated pre-1986 studies indicating the sensitivity of a number of species, at least as seedlings.The experimental studies showed that tree seedling growth is altered at ozone concentrations observed in many areas of the U.S. The 1996 Oxidant Criteria Document also reported that the conclusions of the 1992 Supplement were still valid: (a) The current 1 hour form of the standard (120 ppb) is not protective nor is it an appropriate air quality indicator for plants; (b) the 7- and 12-hr seasonal mean was also not appropriate because of it treats all concentrations equally, and does not consider exposure duration; and (c) the experimental studies suggest the indicator of air quality for a secondary NAAQS should cumulate all hourly concentrations during the daylight hours and weight the higher concentrations (the 'best' weighting scheme or the relative importance of concentration ranges is not known; however, concentrations above 60 ppb are generally considered important in causing effects) As a result of the scientific data, EPA has concluded that the current secondary NAAQS is not protective of ecological resources, in particular protecting crops and forests in the U.S. The proposal for an alternative secondary NAAQS consists of 2 options for consideration. (1) Set the secondary NAAQS equal to the newly proposed alternative primary NAAQS (8 hr maximum of 80-90 ppb with 1-5 exceedances). While this proposal is based on air quality analyses, and is not biologically related to plant response, the attainment of the alternative proposed primary NAAQS would result in a reduction in ozone ------- concentrations harmful to plant growth.; or (2) set the secondary as a 3 month, 12 hr SUM06 of 26,4 ppm hr. This air quality indicator is based on current biological response data and would prevent tree seedling biomass losses and crop yield losses greater that 10% in 50% of the species. The decision regarding which proposal to adopt is still pending. The process of assembling the Criteria Document also reveals the scientific gaps in our knowledge of ozone effects on plants. Primary among these are (a) the problems associated with relating 2-4 year-old tree seedling response to ozone over 1-3 seasons to the response of large, mature trees over a 40-200 year life cycle; (b) the problems associated with extrapolating responses developed on individuals and relating it to the response of a forest stand; (c) how little we know about natural system complexity, wherein the interaction of species, genotypes, and the multitude of past and present environmental influences dictate the eventual response of the species or community in question. We will provide a brief overview of approaches we are taking to address these major areas of uncertainty, 2. PROBLEM: Predicting Effects of Ozone on Long-lived Species overtime, across varied environments, and in assemblage with other species. The knowledge gaps pointed out in the 1996 Criteria Document represent the source of much of the uncertainty in any risk characterization for any atmospheric pollutant, i.e. the ability to predict effects over long time periods, across diverse regions, and across scales of biological organizational. Fig. 1 illustrates three important considerations when extrapolating data: (1) the influential role of size and ontogeny : seedlings to mature trees: (2)complexity of the environment, including belowground processes and trophic groups, species diversity and stand structure: chambers to forests: and (3) the role of time : single season to multiple years. The whole tree in Fig. 1, including the aboveground and belowground tissues, is the experimental subject. Experimental studies have developed quantitative exposure-response functions to a single pollutant, for a single species at the individual to population level, of a given age, size and phenotypic stage, over a single season, and growing in a defined environment, either artificial or natural (lower left hand corner of the triangle). The information gained from each of these exposure experime'nts is then invariably used in risk characterizations for prediction of the effect along one or more of the lines of the triangle, i.e. mature trees, forests, and multiple years. Areas where the experimental data is limited. The prediction of risk to forests over multiple years requires an understanding of all the components of a forest, including the belowground component, species assemblages and age structure, as well as the spatial and temporal factors, in order to make inferences as to how a forest will respond to a given level of pollutant over time. The first concern is the significant limitations in extrapolating effects at the individual level. Changes in physical size as the tree ages, as well as changes in processes that occur with maturation need to be accounted for in making inferences beyond the experimental. For example, forest tree seedlings that range 1-10 ft in height are studied and potential effects are inferred on mature trees that will attain heights of 150 ft. Because of differences in resource utilization in trees of different size and age (Wareing, 1993; Hinckley et al., 1997), we anticipate their being differences in both the ------- mature trees seedlings multiple years Single season chambers COMPLEXITY forests Figure 1. Extrapolation of stressor effects data (e.g. exposure response functions) across: (1) SIZE AND ONTOGENY of species, e.g.from seedlings to mature trees; (2) TIME, including the duration of the stress exposure, season of exposure, year-to-year variation of exposure, phenology and life span of species; and (3) COMPLEXITY, from experimental chamber studies to 'real-world' conditions, artificial soils to native soils with inclusions of intact soil food webs, or single individuals to communities to landscapes uptake of ozone, as well as the ultimate effect on growth, to result from age and size differences. A few studies have indicated such differences in response and attributed those differences to uptake of ozone (Grulke et. al., 1989; Samuelson and Edwards, 1993; Fredericksen et al., 1996). In addition to these age/size limitations, the spatial and temporal distribution of environmental and exposure factors influencing the response to a pollutant need to be considered in the extrapolation. Water and nitrogen availability is of particular importance in understanding response to a pollutant over the species' range. The ozone exposure varies from year to year with weather and emissions. The year-to-year variation of exposure is not accounted for in most experimental data. How to appropriately use available data in risk characterizations is a question for future activity, as well as conducting new research to address questions of the influence biological scale and multiple environmental factors on the effect of ozone. ------- 3. RESEARCH APPROACH - How to address the shortcomings in existing data and reduce uncertainty in extrapolations needed for risk characterization A multifaceted, interactive research approach is required, including experimental and modelling components, with each informing the other. The objective is a mechanistic understanding of each of the lines of extrapolation illustrated in Fig. 1. Plants are exposed to a wide range of spatially distributed environmental extremes, including anthropogenic stressors (Taylor et al., 1994). Plants utilize available resources in order to optimize growth and reproductive outputs under constantly changing conditions. Air pollution represents a stress similar to other stresses in that plant response is the result of both avoidance and tolerance mechanisms (Aliens et al., 1976; Tingey and Andersen 1991). Avoidance factors include factors that restrict the pollutant from reaching the target site, e.g., stomatal closure. Tolerance factors include physiological shifts that result when the stress reaches the target site and causes a metabolic change. Despite the nature of tolerance mechanisms, plant response to a stress involves a shift from one metabolic condition to another, and therefore represents a stress response. BIOLOGICAL logeny •morphology -phenology EXPOSURE \ -biochemistry ^mcencentration -duration -frequency -temporal pattern -tpabal pattern ENVIRONMENTAL -*bfctic(e4.T,ppt) -btotkleg. -outrien ••nib ECOLOGICAL *ll from "Community ••pecies' range •migration -biofeodiemkal eycks EXPOSURE " •coBceDcentratioo -duration -frequency -temporal pattern ••Banal pattern ENVIRONMENTAL •al from "Communilies' •feomorphology -regional hydrotogy -landcover pattern •iymopck dunale ECOLOGICAL -all from "Individ -^accession -competition prtdation -diversity physiognomy EXPOSURE -concencentration -duration -frequency -temporal pattern -cpatial pattern COMMUNITY RESPONSE ENVIRONMENTAL -aO from "Individual" •dope -•spec! •elevation -microdimate = LANDSCAPE RESPONSE C. •= INDIVIDUAL RESPONSE A. Figure 2. Factors influencing plant response to stress at the Individual/Population level (A) or Community/Stand level (B) or the Landscape level (C) ------- There are several factors that affect the magnitude of plant response to ozone stress (Hogsett et al., 1988), and Figure 2 shows a list of important factors at each spatial scale of interest. Figure 2a provides a list of factors influencing individual plant response to ozone, including (1) the past and present environmental milieu, including daily and seasonal weather conditions, soil fertility and water holding capacity, soil food web, and disturbances of various kinds (pest infestations, wildfire, floods, land-use, and so on); (2) biological factors, including genetic variability and constraints, seasonal phenology, limits to biochemical responses and (3) components of stressor exposure, including concentration of pollutants, frequency and duration of exposure, and daily, seasonal and annual patterns of exposure. All of these factors need to be considered to extrapolate the individual plant's response across temporal and spatial scales. Understanding the factors that influence response to ozone stress is considerably more complex at the level of communities (Fig 2B) than at the level of the individual (Fig 2A). Predicting behavior of individuals growing in mixed age, mixed species stands includes all of the factors that influence individual responses, plus such factors as variations in elevation, steepness of slopes, and aspect, and the predictable characteristics of climate near the ground (microclimate). In addition, ecological processes (plant succession from least to most shade tolerant species; competition for sunlight, nutrients, and soil moisture; effects of biodiversity and the physiognomy of vegetation) cannot be derived directly from knowledge of the biological attributes inherent in individual plants. Yet, these processes play a critical rcrle in determining individual plant response to stressors as they form groups of species (populations). Indeed, exposure to stressors must now include a careful review of plant spatial distribution as controlled by the environmental variables included in Figure 2B. The shift from consideration of relatively uniform stands of vegetation gathered into communities, to highly diverse ecosystems gathered into watersheds and landscape units involves additional processes that modify response to stressors (Fig 2C). Environmental effects of geomorphology (landforms, geology, drainage patterns) and regional hydrology (seasonal movement of moisture in soils, into streams and the atmosphere) Combine with synoptic climate variables (air mass frequencies, teleconnected patterns such as el nifio) to generate a regional response which will not equal the average response of populations in the region. Land uses create varying sizes and connectivities of wildlands, each of which functions differently. Related ecosystem processes become significant at this dimension. For example, they define the geographic ranges of animal and plant species distribution, and the daily and seasonal movements (animal migrations) or long-term changes in their boundaries (plant migration). Although the exposure variables at the landscape level are mostly those useful at the level of individual plants, the converse is not true: the landscape-level phenomena can be predicted without much consideration of individual plant characteristics or their environmental constraints. The extrapolation from individual to community to landscape and on to the biome level requires multiple steps (Hinckley et al., 1997). Typically, one works with empirical information generated at one level of organization, using it to model and validate at the next level of organization, and subsequently using that information (model) to predict (implications) at a third level of organization. The complexity characteristic of each level of ------- hierarchy precludes extrapolation beyond 2 or 3 levels of scale. For the purposes of studying ozone effects, it is necessary to bridge scales across each of several levels of complexity: cells-organs-individuals; organs-individuals-stands; individuals-stands-ecosystems; stands- ecosystems-landscapes. Ideally, the complexities described above for each level of scale could be addressed experimentally using experimental designs that are quantitative and enable prediction of effects over a range of concentrations, including the variation in that response (Hogsett et al. 1988). Due to logistical problems, exposure studies over several years are not easily accomplished on whole forests, or even whole trees. Therefore, extrapolation from smaller scaled exposure studies is necessary. To accomplish the extrapolations, we use a combination of process-based and stand-level model simulations. Models can be parametrized for large trees or stands that simulate the interaction of multiple factors in response to ozone based on the experimental studies of the multiple interactions, such as those illustrated in Figure 2. Models are essential components of any risk characterization of ecosystems response to natural and anthropogenic stressors (Rastetter, 1996). Driving forces such as temperature, precipitation, nutrient inputs, topography, soil moisture, and soil biota make it impossible to predict future response to ozone based on single-factor experiments alone (Rastetter et al., 1996). Process- based models at the individual or the stand level can provide a self-consistent synthesis of the results of many experiments. The synthesis provided by these models includes the interactions among physiological and ecological processes that give rise to fhe synergistic responses to multiple environmental factors and stresses (Rastetter et al., 1996). Coupled with modelling efforts, correlative, observational studies can be carried out in forested areas along pollution gradients where monitoring is concurrent and other environmental gradients (e.g. water availability) are present. These studies are possible in some regions of the U.S. (Miller et al., 1982;1989), and will be useful in helping us understand population-level responses as well as validating model simulations. As a way to illustrate how we are combining experimental and modelling efforts, several examples of our current research activity in Corvallis are discussed. We have associated each example with a particular level of biological organization, i.e., individual/population, community, and landscape; however, there is overlap among the levels both in the experimental information obtained and its subsequent use. 3.1. Individual/Population 3.1.1. Extrapolation of effects on a few species to reflect the effect on all species A major portion of the research activity over the past decade has been directed at developing exposure-response functions for both crops and forest tree species (Heagle et al., 1988; Fuhrer et al., 1989; Hogsett et al,, 1989; Tingey et al. 1991; Neufeld et al., 1995; Kamosky et al., 1995.). Many of the studies were designed to predict growth and yield effects at the individual and population level at a range of ozone concentrations. The primary means of exposing plants during the last 15 years has included the use of the open top chamber methodology for developing the empirical database of ozone effects on crop and ------- tree species (Heagle et al., 1988; Taylor et al., 1994). Open-top chambers provide a sufficient range of ozone treatments for quantifying effects (i.e. exposure-response functions), and still provide growing conditions that closely match those in the plants' natural growing environment. In addition, the multitude of chamber sites has permitted adequate replication of studies from year to year. Examples of these studies include the National Crop Loss Assessment Network (NCLAN) (Heck et al, 1991), the European Open-Top Chamber studies (Weigel et al., 1987), and several studies of tree species (Hogsett et al.,1989; Matyssek et al., 1993; Neufeld et al., 1995; Karnosky et al., 1995; Rebbeck, 1996). The studies included a wide range of response functions illustrating that ozone response is species dependent. For the first few years of the our project studying forest tree species, we developed exposure response functions for 11 forest tree species during 1-2 year exposure seasons using open-top chambers at locations across the U.S. (Table 1). Sites were selected in different regions of the U.S. where the species of interest were indigenous. A common research protocol was used at each site, and several of the species were examined at more than one site. Table 1. Tree Species: Locations and exposure durations used to construct composite exposure response function. Species/family Aspen — wild Douglas fir Ponderosa pine Red alder Black cherry Red maple Tulip popular Virginia Pine Loblolly Sugar Maple E. white pine Site Oregon Michigan Oregon Oregon Oregon Oregon Oregon Oregon Smokey Mm Nat'l Park Smokey Mm Nat'l Park Smokey Mm Nat'l Park Smokey Mrn Nat'l Park Ohio Smokey Mm Nat'l Park Alabama Michigan Ohio . Michigan Ohio Exposure duration (days) 112 98 234 234 230 280 118 112 76 140 55 184 222 98 555 180 222 180 222 ------- Figure 3 shows examples of exposure-response functions for total biomass (foliage, stem and root) after one growing season. Douglas-fir (Psuedotsuga menziesi (Mirb.) Franco) seedlings were relatively insensitive to ozone (Fig. 3a), while ponderosa pine (Pinus ponderosa Laws.) (Fig 3b) and quaking aspen (Populus tremuloides Michx.) (Fig. 3c) were relatively sensitive to ozone. The response of the other 7 species varied from insensitive (e.g. Virginia pine, Pinus virginiana Mill.) to moderately sensitive (e.g. Sugar Maple, Acer saccharum Marsh.) to sensitive (e.g. Black cherry, Prunus serotina Ehrh.). Given the range in response among species, it is difficult to generalize or extrapolate an effect across all species in order to provide information for decision making regarding risk characterization and standard setting. However, one approach to accomplishing such a task is to develop a composite exposure-response function for tree species and one for crops. Figure 4 shows the composite functions for crops and tree species. The percentile distributions illustrate the range of predicted relative biomass loss as a result of ozone exposure. At 30 ppm-hr ozone, the predicted biomass loss of tree seedlings ranges from 0-20%. Douglas Fir Ponderosa Pine Aspen i it * * to SUMMIppm*) Figure 3. Exposure-response functions for Douglas-fir, Ponderosa Pine and Aspen. Total biomass response after one season of exposure. Exposure is expressed as SUM06 (sum of all hourly values equal to or greater than 60 ppb during the daylight hours (12 hr) over the growing season) The Weibull or linear exposure-response model was used to relate final harvest biomass and exposure, expressed as the 12-h SUM06 index for 31 NCLAN crop studies (Fig 4A) and 26 tree seedling studies (Fig 4B). Separate regressions were calculated for studies with multiple harvests or cultivars, resulting in a total of 54 individual equations from the 31 NCLAN studies and 56 equations from the 26 seedling studies. For crops, each equation was used to calculate the predicted relative yield loss at 10,20, 30, 40, 50, and 60 ppm-h, and the distributions of the resulting loss were plotted (Fig 4A). For seedlings having been exposed to ozone for one or multiple exposure seasons of varying durations, each equation was used to calculate the predicted biomass loss at 10, 20, 30, 40, 50, 60, and 70 ppm-h per 92 days a year, and the distributions of the resulting loss were plotted (Fig 4B). ------- 0.7 ,0,6 I O.S 10.4 ! 0.3 A. CROPS PRYL. 1 -oqp|-(SUM06/87.42449)"l. 8209501 10 30 40 12 kr SUMO6 (ppm-hr) 80 60 B. TREE 8EEDLIHQ8 PRBL - 1 • exp[-tSUM06/!4<3-S«2&3**1.30G??9j 10 30 4O SO 12 hi 8UMO6 (ppm-hr) 60 Figure 4. Composite exposure response functions for crops (A) and tree seedlings (B). Median response (50th percentile) is illustrated with solid line and weibull model given in each figure, 25th and 75th percentile shown with shadowed bars at each exposure value. At each SUM06 value, the median (50th percentile) response was used to characterize the typical plant response to ozone across studies. The composite curve for crops or seedlings is the calculated Weibull model fit to the median response points. Using these composite functions which composit species response, one can select an ozone concentration for a standard that is protective of trees or crops, depending on the level ------- of protection sought. The line represents the Weibull curve that describes the response in 50% of the seedlings or crops. Using this as an estimation of the mean response in tree species, concentrations can be determined that are protective against 10,20, or 30% biomass losses (or crop yield losses), in 50% of the species. Although this approach provides some means to estimate the impact of ozone on crop and forest species, it includes large uncertainties because of the species-dependence of the response, and it does not account for the variation in the response resulting from different growing environments, exposure or other biological/ecological considerations (Fig. 2). 3,1,2, Extrapolation of effects along the line of increasing complexity As controlled experimental studies begin to incorporate the multitude of factors typical of natural systems, (e.g. temporal and seasonal variation in exposure concentrations, variation in availability of water and nitrogen, presence of intact soil food web and rhizosphere biota, presence of competing species, etc), an understanding of both the direction and the magnitude of the modification by these influences on the exposure-response will result. The information can be used to reduce uncertainty in risk assessments, as well as improving process-based models simulating long-term effects of ozone on forest communities. Many earlier studies have examined the influence of these factors, including water stress on crops(e.g.Tingey & Hogsett, 1985; Miller et al. 1989; Temple et al 1988.), relative humidity (McLaughlin & Taylor, 1980), nutrient status (e.g. Tjoelker & Luxmore 1991; Bytnerowicz et al, 1990), exposure components (Walmsley et al., 1980; Musselman et al., 1983;1986; Hogsett et al., 1985; 1988; Hogsett and Tingey, 1991). Three examples of ongoing studies at our lab looking at this question from the perspective of the individual, community (Sec. 3.2.1), and landscape (Sec. 3.3) (Fig 2) are presented, Nitrogen Status: In a study of the influence of nitrogen supply on the response of ponderosa pine to ozone, we have measured both an individual response and a 'system' level response (Andersen and Scagel, 1997). Three levels of nitrogen- full, 2/3 or 1/3 rate was provided the year prior to exposure and continued through the year of exposure. The exposures were continued for 2 years. Figure 5 shows the response in biomass and foliar nitrogen content after one year of ozone exposure as not obviously different with varying nitrogen. The observed response suggests there is not a significant influence of nitrogen status on the growth response to ozone, and thus it may be appropriate to extrapolate exposure-response functions describing biomass at the level of the individual across regions differing in nitrogen availability. One measure of system level response, however, CO2 flux from the soil, did indicate an influence of nitrogen status. Figure 5 shows a more pronounced response to ozone at the low nitrogen level, suggesting a possible interaction. The low nitrogen-status may influence the ability to maintain roots and increased root mortality may be contributing to increased metabolic activity of soil biota and thus account for the increased CC>2 flux from soil. These studies are not completed and are not thoroughly analyzed, but do suggest ------- possible implications in our ability to extrapolate exposure response Functions spatially as a result of available nitrogen. The measured response, biomass loss or CO2 flux from soil, are both valid and provide critical data for risk characterization. However, from this study two different conclusions might be drawn on the role of nitrogen depending on whether the endpoint for assessment was at the level of the individual, or at the system level. A. BIOMASS - LEAF, STEM & ROOT B. FOLIAR NITROGEN 100 80 3 g 60 !«• Q 20 100% 67% 33% 1.6 1.4 Leaf z Stem | Root 1.2 1.0 0.8 0.6 0.4 CF •-. Ep-23 Ep-31 0 90120 0 90120 0 90120 Ozone Exposure (ppm-hr) 100% 67% • 33% C. SOIL CO2 FLUX 8 -7-3 £5^ o 4. IM 8,. 0 1/3 Nitrogen 2/3 Nitrogen Full Nitrogen Control EP90 EP120 Ozone Treatment Level Figure 5. Response of ponderosa pine to 3 levels of ozone (Episodic [Ep] regimes having a SUMO of 90 & 120 ppm hr and charcoal filtered[CF]) under 3 levels of nitrogen (100%, 67% & 33% of precribed application) in native soils. (A) biomass response after one season. (B) Total foliar nitrogen after one season. (C) Soil CC>2 flux after one season. ------- 3.1.3. Extrapolation of effects from seedlings to mature trees: Understanding the role of size and age in determining tree response to ozone is very important. As discussed above, much of the data we have today was obtained on various species of trees in the seedling stage, usually 1-5 years of age. As a result, the uncertainty in using this data to extrapolate the effects on large mature trees is very significant for any risk assessment. How maturation (ontogeny) effects tree response to air pollution is not known. It is thought that changes in pollutant uptake may occur with size, as well as differences in carbon sink sizes, the size of respiratory pools, and priorities (such as growth versus maintenance respiration) may affect the way the tree responds to ozone. Very few studies have addressed this question due to the technical difficulties of exposing large trees to ozone (Samuelson and Edwards, 1993; Grulke et al., 1989; Fredericksen et al., 1995). Studies in Corvallis are planned focusing on resource utilization (C, N, H20) by trees (e.g. Ponderosa pine) as a function of size and age. The research will develop a mechanistic basis for scaling (Wareing, 1993), and will provide data to parametrize models that will incorporate age differences in resource utilization (Hinckley et al., 1997) affected by ozone, as well as incorporating gas exchange differences which effect ozone uptake. 3.2. Community 3.2.1. Extrapolation of effects along the lines of complexity: Presence of other species and the soil food web: In the second example, studies at Corvallis have indicated the influence of intact soil food webs and the importance of including belowground complexity in when studying tree response to ozone. Mesocosms in open-top chambers were planted with varying densities of blue wild rye (Elymus glaucus Buckl.) (a natural competitor of ponderosa pine seedlings) and a single ponderosa pine seedling in soil from a ponderosa pine forest (Andersen, et al. 1997). The response to ozone exposure was measured as changes in carbon dynamics, including allocation in individuals (species) and soil respiration as C02 flux, over 3 years. This was not a true competition study design because it lacks the necessary permutations, however, it does indicate the influence of increasing complexity of the environment on the exposure response. The preliminary data would suggest that there is not an interaction between plant density and ozone. Figure 6 shows an increase in C02 flux and soil organic matter as grass density and ozone concentration increases. However, what we are learning from this study is the importance of an intact soil food-web, which was included here with the native ponderosa pine forest soil. Earlier exposure studies of carbon allocation belowground, both from our lab and results from other laboratories, showed a decrease in allocation of carbon belowground with ozone exposure (Gorissen and van Veen, 1988; Spence et al. 1990; Andersen and Rygiewicz,1995). ------- Grass 0 Grass 1 Grass 2 Grass 0 Grass 1 Grass 2 Level of Grass Competition Figure 6. Soil CO2 flux (A) and Soil organic matter (B) following exposure to ozone for 1 year in mesocosms containing pine seedling and blue stem rye grass at 3 increasing densities ("Grass 0,1,2"), Bars indicate 95% confidence interval. Differences were significant at the 5% level From these individual plant studies we would have predicted that ozone exposure would decrease C02 flux from soil due to less allocation of carbon belowground, rather than increase it as shown in Figure 6. The increase in both soil organic matter and soil COi flux may be the result of increased root turnover and subsequent degradation by soil organisms in ozone exposed systems, which could increase total soil COi flux. If this is the case, perhaps with' time a decrease in flux will be observed as the seedlings get older and the total root system size decreases. Exactly how this information can be included in extrapolation of results for risk characterization is not known, however it does indicate that when studies are accomplished in simple systems and over short time frames (1-3 years), some degree of uncertainty should be included. ------- 3.3 Landscape level We have employed an approach for extrapolating tree biomass and crop yield response functions across each species' range to determine the extent and magnitude of possible ozone risk to ecological resources in the U.S. (Hogsett et al., 1997). The approach uses a geographical information system (GIS) to integrate all the necessary data in a spatial context., including (1) estimated exposure over a growing season, (2) characteristics of the growing environment (e.g. precipitation, temperature, soil water, soil nitrogen levels, etc), (3) species' distribution, (4) species' inventory data (e.g. standing biomass, productivity) and (5) landscape features (e.g. elevation) (Figure 7). ------- Potential Exposure Estimated Exposure Exposure Effects Risk Characterization J' I- »0 MO 1M 300 310 >K *30 MO 3M MO n>t«ntial Expamir* VWu* PWBl- 1 ««p( h»IJB/BI I » w •» to m Expocun (SUM06 ppm-hl Figure 7. The GIS data layers for characterizing risk of ozone to tree species, including GIS-based estimation of ozone exposure in non-monitored areas for a given year, seedling or model-simulated tree exposure exposure response functions for each species, and the spatial data on environmental factors such as precipitation, and species density. The combination of these layers results in a new map of risk. ------- The approach is useful for distributing experimentally-derived exposure-response functions from previous crops or tree seedling exposure studies, as well as extrapolating process-based whole-tree model simulations of large trees using TREGRO (Weinstein et al. 1991) over a 3 year time period incorporating the role of environmental influence such as water or nitrogen, or stand-level model simulations of community structure (ZELIG) incorporating very long-term effects (100 years) as a result of ozone, climate, soil, water and nutrient availability. Theses extrapolation are based on the spatial and temporal distribution of ozone exposure and climate variables, e.g., precipitation and temperature that are present across the species' range or .the region of the forest-type. The results of spatially distributing seedling exposure-response functions and then area-weighting that response is given in Figure 8. A sensitivity-ranking based on ozone sensitivity of each species and the potential area of the species impacted can be made (Hogsett et al 1997). 1988 T . ¥ 7 CD f 1989 •i Fl OWTV %• U M SUH° «. Sptclti Spcelt* Figure 8. Area-weighted annual biomass loss in seedlings from 1988 and 1989 ozone exposure. Percentile distributions of species' area is shown for 10th, 25th, 50th, 75th and 90th percentile. The response for 50% of the species' area is shown as the midline in the shadowed box. The 25 and 75% of the area are the ends of the shadowed box. To illustrate the process-based model simulation approach, Figure 9 shows the response in loblolly pine biomass after 3 years of exposure. We compared a hot, dry year having relatively high ozone concentrations (1988) to a cooler, wetter year having relatively low ozone concentrations (1989). TREGRO was parametrized from available literature on ------- 8 Loblolly Pin* 1988 HM(Jx) 8UM08 (ppn-hr) P«rt»nt at Total Growing Ana In Each Map Pass I" 1* 1" I" 0-4 10-19 >-2« ttmt Yttf ProduclMy Induction (pd.) LoNo«yP(fw19B9 tHI(Sx) 8UM06 (ppm-hf) Ptrant of ToM Growing Area hi Etch Map Class to TO 8* • N 1- "5 30 I W ThfM mm O-J lYwr • Produdrtfty Riductton (pd.) Figure 9. Regional extrapolation of TREGRO-generated aboveground biomass response functions. Two ozone exposure years were compared to illustrate response to a year of high ozone exposure (1988) and a more typical year of ozone exposure (1989) (C). Simulations were for 3 years with the model parametrized for a 30 year-old loblolly pine. The range of response possible across the region of loblolly was accomplished using.climate data from 3 sites: a relatively wet area (Cumberland Co., NC), a relatively dry site (Cado Parish, LA), and an average site (Shelby Co., AL). From the response functions generated for each of these 3 sites (A), biomass loss functions were derived (B) and combined with estimated ozone exposures for 1988 or 1989 (C) and the loss interpolated across a region expected to have the same climate (D). Histograms (E) illustrate the area of loblolly experiencing 0-s»2% loss over the 3 year simulation. ------- large mature loblolly pines and used meteorological and ozone exposure data from 3 sites across the range of loblolly. The sites represented a relatively drier region (Cado Parish, Louisiana), a relatively wetter region (Cumberland Co. North Carolina) and an average site (Shelby Co., Alabama) based on a 10 year average of all met sites in this region. Using 3 sites for the simulations allows a means to demonstrate the range of response to ozone across the species' range based on the meteorological variables of precipiation and temperature. The biomass response functions resulting from the 3 year runs (1988-1988-1988 and 1989-1989- 1989) for each met site are interpolated over the region that is thought to be experiencing the same climate. Cado Parish response functions for the two 3-year scenarios are shown in Figure 9A. The biomass loss function (Fig 9B) is derived from the shoot biomass response function (Fig 9A) and is applied across the appropriate region with the estimated ozone exposure (Fig 9C) and the resulting biomass loss maps are generated (Fig 9D). The histograms (Fig 9E) show the area of loblolly predicted to experience various reductions in growth over the 3 year period. Clear differences are apparent comparing the higher ozone exposures of 1988 to the 1989 scenario of relatively low ozone. Sixty percent of the loblolly area is predicted to experience 1-2% reductions in productivity over the three years of 1988 ozone, compared to 20% of the area when 1989 exposure scenarios are used. The TREGRO response function generated using the climate data for the relatively wet year of 1989 compared to the function generated using the relatively dry year predicts changes in growth at less ozone concentration (150-200 ppmhr for 3 year exposure) than that from the dry year of 1988 (250-300 ppm-hr for 3 years) (Fig 9A) suggesting the role of precipiation in the response to ozone. In the drier year, the reduced conductivity would result in less ozone being taken up and thus higher exposures are needed to cause effects. To illustrate the effect of precipitation on the biomass loss across the region of loblolly pine, the response function from 1989 is used with the estimated ozone exposures of 1988 and 1989 (Figure 10). With the higher ozone exposure of 1988 and the greater ozone response simulated using the higher precipitation there is a greater loss exhibited than when this same response function is used with the 1989 estimated ozone exposure (Fig 10D). If the same ozone exposure year (1988) is compared using the 3 year 1988 and 3 year 1989 generated response function, there is substantially more area experiencing greaeter than 2% productivity loss with the wet year (1989)-generated response function (Fig 9D and 10D). This simulation is in keeping with experimental data demonstrating the reduced effect of ozone under drought conditions (Tingey and Hogsett 1985). Although there are higher exposures in 1988, because of the higher temperatures, and less precipitation, there is presumably less ozone uptake due to stomatal closure. The model simulated response function demonstrated this. These results also suggest a possible weighting scheme as year to year variation in exposure is taken into consideration for exposure-response functions and developing air quality indices. The result of such an analysis in a risk characterization is that less weight would be given to warm years where precipitation is low compared to cooler, wetter years - (even though less ozone exposure would be anticipated). ------- Ozont Expown 1988 D LoWolly Pint 1988 P«rwnt of Total growing Ana In Each Map Class 0-j ian >-ao T1li»» Yttr PredudlTtty Reduction put.) Ptmnt of Totil Growing Ana In Each Map Class I • M 1- ~5 X £ » O-J 10-lt >-lt Thin Y«» Produdtttty Rtducttcn (pa.) Figure 10. The role of Precipitation in aboveground biomass response to ozone. TREGRO generated response function of shoot biomass response for a relatively wet year of 1989 (A) was used to derive a biomass loss function (B) and that was applied to 1988 and 1989 estimated ozone exposure (C) across the region of loblolly using the GIS. The loss was interpolated across the region expected to have the same climate to produce maps (D) and from the maps, histograms reflect the area of loblolly predicted to experience various reductions in productivity over 3 years (E), ------- Using 3 different sites having extremes and average climates was our attempt at addressing the influence the different climates/meterorology have on the response to ozone across a species' range. As expected, the result of such simulations illustrate a range of responses based on climate and ozone, but an effort is needed to formally incorporate the uncertainty in the extrapolation of the effect of ozone on each species. Using the GIS as a framework for integrating layers of interacting factors in the exposure-response (i.e. exposure, environment, biological/ecological) may provide a means to accomplish these extrapolations that are needed in characterizing the risk of a spatially and temporally distributed atmospheric pollutant like ozone. Each pixel in the area or region of interest in a risk assessment can be characterized by including each of the factors contributing to the exposure-response at the appropriate spatial and temporal scale, and then the response generated with simulation modeling at either the individual level (species) or the stand-level (species assemblages). 4. CONCLUSIONS Both controlled laboratory and field studies are planned and are essential for providing the mechanistic basis underlying the modelling and extrapolation approach. In order to understand the complexities of spatial and temporal components of ecosystem response, both field and lab studies must progress in parallel. Collaborative research efforts among modelers and experimentalist are essential to improve models and guide data collection. As models improve, our ability to extrapolate should improve and thus reduce the uncertainties in risk characterization. The models need to be accurate and sensitive enough to capture the important factors that modify the response to ozone. The experimental data are essential for providing the mechanistic underpinnings of the response models at all scales. REFERENCES: Andersen, C.P. and P.T. Rygiewicz, 1995. Effects of ozone on temporal allocation of carbon in mycorrhizal Pinus ponderosa seedlings. NewPhytol. 131:471-480. Andersen, C.P., and C.F. Scagel, 1997. Nutrient availability alters below-ground respiration of ozone exposed ponderosa pine. Tree Physio!., in press. Ariens, E.J., A.M. Simonis, and J. Offermerier. 1976. Introduction to general toxicology. Academic Press, NY. Fredericksen, T.S., J.M. Skelly, K.C. Steiner, T.E. Kolb, and K.B. Kouterick. 1996. Size- mediated foliar response to ozone in black cherry trees. Environ. Pollut. 91:53-63. Fuhrer, J.D., A. Egger, B. Lehnherr, A. Grandjean and W. Taschannen. 1989. Effects of ozone on the yield of spring wheat (Triticum aestivum L., cv. Albis) grown in open-top field chambers. Environ. Pollut. 60:273-289. Fuhrer, J.D., A.G. Grimm, W. Tschannen and H. Sariat-Madari. 1992. The response of 24 ------- spring wheat (Triticum aestivum L.) to ozone at higher elevations: II. Changes in yield, yield components and grain quality in response to ozone flux. New Phytol. 121:211-219. Gorissen, A. and J.A. van Veen. 1988. Temporary disturbance of translocation of assimilates in Douglas-fir caused by low levels of ozone and sulfur dioxide. Plant Physiol. 88:559-563. Grulke, N.E., P.R. Miller, R.D. Wilborn, and S. Hahn. 1989. Photosynthetic response of giant sequoia seedlings and rooted branchlets of mature foliage to ozone fumigation. In Effects of Air Pollution of Western Forests, R.K. Olson and A.S. Lefohn, eds., AWMA, Anaheim, CA. Heagle, A.S. L.W. Kress, P.J. Temple, RJ. Kohut, J.E. Miller, and H.E. Heggestad. 1988. Factors influencing ozone dose-yield response relationships in open-top field chambers studies. IN: Heck, W.W., O.C. Taylor and D.T. Tingey, eds. Assessment of Crop Loss from Air Pollutants: Proceedings of an International Conference, pp 107-138. Elsevier Applied Science. New York. Heck, W.W., A.S. Heagle, J.E. Miller and J.O. Rawlings. 1991. A national program (NCLAN) to assess the impact of ozone on agricultural resources. IN: Berglund, R.L., D.R. Lawson and D.J. McKee, eds. Tropospheric Ozone and the Environment: Papers from an Interantional Conference, pp 225-254. Air and Waste Management Association, Pittsburg, Pa. (A&WMA Transaction Sens no. TR-19). Hinckley, T.M., D.G. Spruegel, J.R. Brookds, K.J. Brown, T.A. Martin, D.A. Roberts, W. Schaap and D. Wang. 1997. Scaling and integration in trees. IN: Peterson, D.L. and V. T. Parker, Eds. Ecological Scale: Theory and Application. Columbia University Press, New York. Hogsett, W.E., D.T. Tingey and S.R. Holman. 1985. A programmable exposure control system for determination of the effects of pollutant exposure regimes on plant growth. Atmos, Environ. 19: 1135-1145. Hogsett, W.E., D.T. Tingey and E. H. Lee. 1988. Ozone exposure indices: Concepts for development and evaluation of their use. IN: Heck, W.W., O.C. Taylor and D.T. Tingey, eds. Assessment of Crop Loss from Air Pollutants: Proceedings of an International Conference, pp 107-138. Elsevier Applied Science. New York. Hogsett, W.E., D.T. Tingey, C. Hendricks and D. Rossi. 1989. Sensitivity of western conifers to SO2 and seasonal interaction of acid fog and ozone. IN: Olson, R.K. and A.S. Lefohn, eds. Effects of Air Pollution on Western Forests. Transaction Series no 16. pp 469-489. Air and Waste Management Association, Pittsburg, PA.IN: Olson, R.K. and A.S. Lefohn, eds. Effects of Air Pollution on Western Forests. Transaction Series no 16. pp 309-323. Air and Waste Management Association, Pittsburg, PA. Karnosky, D., Z.E. Gagnon, R.E. Dickson, M.D. Coleman, E.H. Lee, and J.G. Isebrands. 1995. Changes in growth, leaf abscission and biomass associated with seasonal ozone exposures of Populus'tremuloides clones and seedlings. Canadian Journal Forest Science 25 ------- Matyssek, R., MS. Gunthardt-Goerg , W. Landolt, and T. Keller. 1993. Whole plant growth and leaf formation in ozonated hybrid poplar (Populus x euramericana). Environ, Pollul. 79:1-7. Miller J.E., A.S. Heagle, S.F. Vozzo, R.B. Philbeck and W.W. Heck, 1989. Effects of ozone and water stress, separately and in combination on soybean yield. J. Environ. Qual. 18: 330-336 Miller, P.R., O.C. Taylor and R.G. Wilhour. 1982. Oxidant air pollutant effects on a western coniferous forest ecosystem. U.S. EPA, Environmental Research Laboratory, Corvallis, Oregon EPA-600/D-82-276. Miller, P.R., J.R. McBride, S.L. Schilling and A.P. Gomex. 1989. Trend of ozone damage to conifer forests between 1974 and 1988 in the San Bernardino Mountains of southern California. IN: Olson, R.K. and A.S. Lefohn, eds. Effects of Air Pollution on Western Forests. Transaction Series no 16. pp 309-323. Air and Waste Management Association, Pittsburg, PA. Musselman, R.C., R.J. Oshima, and R.E. Gallavan. 1983. Significance of pollutant concentration distribution in rthe repsonse of'red kidney' beans to ozone. J. Am. Soc. Hort.Scl 108:347-351. Musselman, R.C., A.J. Huerta, P.M. McCool, and R.J. Oshima. 1986. Response of beans to simulated ambient and uniform ozone distributions with equal peak concentration. J, Am. Soc. Hortic. Sci. 111:4 70-473. McLaughlin, S.B., and G.E. Taylor, Jr. 1980. Relative humidity: Important modifier of pollutant uptake by plants. Science 211:167-169. Neufeld, H.S., E.H. Lee, J.R. Renfro, W.D. Hacker, and B.H. Yu. 1995. Sensitivity of seedlings of black cherry (Prunus serotina Ehrh.) to ozone in Great Smokey Mountains National Park: 1. Exposure-response curves for biomass. New Phytologist 130:447-459. Rebbeck, Joanne. 1996. Chronic ozone effect on three notheastern hardwood specie: growth and biomass. Canadian Journal Forest Science. 26:1788-1798. Samuelson, L.J., and G.S. Edwards. 1993. A comparison of sensitivity to ozone in seedlings and trees of Quercus rubra L. NewPhytoI. 125:373-379. Spence, R.D., E.J. Rykiel Jr., and P.J.H. Sharpe. 1990. Ozone alters carbon allocation in loblolly pine: Assessment with carbon-11 labeling. Environ. Pollut. 64:93-106. Taylor, G.E., D.W. Johnson, and C.P. Andersen, 1994. Air Pollution and Forest Ecosystems: A Regional to Global Perspective. Ecological Applications 4:662-689. Temple, P.J., L.F. Benoit, R.W. Lennox, C.A. Reagan and O.C. Taylor. 1988. Combined effects of ozone and water stress on alafalfa growth and yield. J. Environ. Qual. 17: 108-113. Tingey, D,T. and W.E. Hogsett. 1985. Water stress reduces ozone injury via a stomatal mechanism. Plant Physiol. 77:944-947. Tingey, D.T. And C.P, Andersen. 1991. The physiological basis of differential plant sensitivity to changes in atmospheric quality. In Ecological Genetics and Air Pollution, G.E. Taylor Jr., L.F. Pitelka and M.T. Clegg, eds., Springer-Verlag, NY.' 26 ------- Tingey, D.T., W.E. Hogsett, E.H. Lee, A.A. Herstrom and S.H. Azevedo. 1991. An evaluation of various alternative ambient ozone standards based on crop yield loss data. IN: Berglund, R.L., D.R. Lawson and D.J. McKee, eds. Tropospheric Ozone and the Environment: Papers from an International Conference, pp 272-288. Air and Waste Management Association, Pittsburg, Pa. (A&WMA Transaction Series no. TR-19). Tjoelker, M.G. and R.J. Luxmoore. 1991. Soil nitrogen and chronic ozone stress influence physiology, growth and nutrient status ofPinus taeda L. and Liriodendron tulipifera L. seedlings. New Phytol. 119:69-81. U.S. EPA. 1986 Air Quality Criteria for Ozone and Related Photochemical Document, Vol III.EPA/600/8-84-020cF. Washington, D.C. U.S. EPA. 1992. Summary of selected new information on effects of ozone on health and vegetation: Supplement to 1986 Air Quality Criteria for Ozone and Related Photochemical Oxidants. EPA/600/8-88/105F. U.S. EPA. 1996.Air Quality Criteria for Ozone and Related Photochemical Documents, Vol II. EPA/600/P-93/004bf. Washington, D.C.. Walmsley, L., M.R. Ashmore, J.N.B. Bell 1980. Adaptation of radish Raphanus sativus L. in response to continuous exposure to ozone. Environ. Pollui. 23: 165-177. Waring, R.H. 1993. How ecophsyiologists can help scale from leaves to landscapes. IN: Ehleringer, J.R. and C.B. Field, Eds. Scaling Physiological Processes: Leaf to Globe, pp 159-166. Academic Press, New York. Weigel, H.J., G. Adaros and H.-J., Jaeger. 1987. An open-top chamber study with filtered and non-filtered air to evaluate the effects of air pollutants on crops. Environ. Pollut. 47: 231-244. Weinstein, D.A., R.M. Beloin, R.D. Yanai. 1991. Modeling changes in red spruce carbon balance and allocation in response to interacting ozone and nutrient stresses. Tree Physiol. 9:127-146. 27 ------- The information in this article has been funded by the U.S. Environmental Protection Agency. It has been subjected to the Agency's peer and administrative review, and it has been approved for publication as an EPA document Mention of trade names or commercial products does not constitute endorsement or recommendation for use. ------- NHEERL-COR-2166A TECHNICAL REPORT DATA (Please read instructions on the reverse before com/ 1. REPORT NO. EPA/600/A-97/090 2. 4. TITLE AND SUBTITLE Ecological effects on tropospheric ozone: a U.S. perspective - past, present and future. 7. AUTHOR(S) W.E. Hogsett, Christian P. Andersen 9. PERFORMING ORGANIZATION NAME AND ADDRESS U.S. EPA National Health and Environmental Effects Research Laboratory Western Ecology Division 200 SW 35th Street Corvallis, Oregon 97333 12. SPONSORING AGENCY NAME AND At US EPA ENVIRONMENTAL RESE/ 200 SW 35th Street Corvallis, OR 97333 DDRESS ^RCH LABORATORY 3. F 5. REPORT DATE 6. PERFORMING ORGANIZATION CODE 8. PERFORMING ORGANIZATION REPORT NO. 10. PROGRAM ELEMENT NO. 1 1 . CONTRACT/GRANT NO. 13. TYPE OF REPORT AND PERIOD COVERED 14. SPONSORING AGENCY CODE EPA/600/02 15. SUPPLEMENTARY NOTES: 5th US-Dutch International Symposiuml 1 6. Abstract: Understanding the effects of tropospheric ozone on vegetation, as called for in the U.S. Clean Air Act, has involved collection of experimental data at the species level and, in particular, at the level of the individual populations. Frequently the studies have been regression designs involving single species and single pollutants resulting in quantitative exposure-response functions that characterize the effects on biomass or reproduction (crop yield). Recently, the 1 996 EPA Oxidant Criteria Document reviewed the published research on oxidant effects on crops, forests, and ecological resources, and concluded that the current secondary National Ambient Air Quality Standard (NAAQS) for ozone is neither protective or appropriate. The subsequent discussion and decision on what is an appropriate form and level of the secondary NAAQS utilized results from single pollutant and single species studies. To estimate concentrations causing effects composite response-functions for crops and trees from a large number of species were developed, predicting crop yield loss or annual biomass loss in seedlings. This approach is useful based on available quantitative data on biological effects; however it assumes that individual plant response does not change in the presence of other stresses or in natural systems which are considerably more complex (e.g. species' assemblages and competition for resources). For example, the importance of biological complexity is illustrated when we examine the below-ground ecosystem, an often overlooked portion of the ecosystem. We have found that incorporating natural biological complexity into potting soils can result in carbon fluxes opposite those predicted from individual plant studies using artificial media lacking natural soil foodwebs. A future research approach to understanding 03 effects on ecosystems is required that will develop the necessary linkages to extrapolate experimental data taken at the individual level, often in artificial conditions, to predict changes on individuals or populations in more complex native environments. We will present experimental and modelling activities from our laboratory that show how we are beginning to address the problems of scale, complexity, and multiple stresses in forested ecosystems exposed to ozone stress. 17. a. DESCRIPTORS Soil foodwebs, below-ground ecosystem, carbon fluxes, ozone stress 18. DISTRIBUTION STATEMENT KEY WORDS AND DOCUMENT ANALYSIS b. IDENTIFIERS/OPEN ENDED TERMS 1 9. SECURITY CLASS (This Report) 20. SECURITY CLASS (This page) c. COSATI Field/Group 21. NO. OF PAGES: 25 22. PRICE EPA Form 2220-1 (Rev. 4-77) PREVIOUS EDITION IS OBSOLETE ------- |