Abstract
A dry deposition model was employed to estimate air pollutant uptake by Sacramento’s urban forest. Assuming 1990 air pollutant concentrations, model simulations estimated that approximately 1,457 metric tons of air pollutant are absorbed annually, at an implied value of US$28.7 million. The growing season daily uptake for ozone was approximately 2.4 metric tons per day, while particulate matter (< 10 μ diameter, PM10) uptake was slightly greater, at 2.7 metric tons per day. Daily uptake of NO2 and particulate matter represented 1 % to 2% of anthropogenic emissions for the county. Estimated growing-season annual air pollutant uptake rates averaged 10.9 kg/(ha land area per yr) for the entire study area, 13.9 kg/(ha land area per yr) for urban areas and 4.2 kg/(ha land area per yr) for rural areas. Pollutant uptake rates decreased with decreasing tree canopy cover, along an urban-to-rural gradient.
Tree planting in Sacramento, California, and in other urban areas provides energy and air-quality benefits by direct shading of buildings (reducing energy demand for cooling), by cooling the atmosphere through transpiration of water from leaves, and by the direct absorption of air pollutants by leaf surfaces (Landsberg 1981; Akbari et al. 1992; Rosenfeld et al. 1995; Sailor 1995; McPherson et al. 1994, 1997a, 1997b; Simpson and McPherson 1996; Taha 1996; Simpson 1998). As urban forestry has expanded from the notion of municipal street tree management to urban ecosystem management, new partnerships among local government, electric utilities, and volunteer associations have formed. Local air-quality management districts represent a potential partner to the extent that healthy urban forests provide air-quality benefits. The aim of this research was to produce estimates of both annual air pollutant uptake by Sacramento’s urban forest and its economic value. This study was one component of the Sacramento Urban Forest Ecosystem Study (SUFES), whose goal is to determine relationships between urban forest structure and function and associated benefits and costs. These and other SUFES results are being used as inputs for a geo-referenced cost-benefit analysis of the region’s urban forest.
Methods
Overview
The Sacramento metropolitan area is a rapidly urbanizing region approximately 120 km (75 mi) northeast of San Francisco, located in California’s Central Valley. The air basin is bounded on the east by the Sierra Nevada and on the west by the Coast Range mountains. Topography and meteorological conditions, combined with growing emissions, contribute to episodes of poor air quality. Summer months are warm and virtually cloud-free, with light winds, relative humidities less than 20%, and days that can exceed 38°C (100°F) (NOAA 1990). Prevailing wind direction is southerly, except in spring and autumn, when it is northerly. During warm months, cool oceanic air often flows into the valley through a sealevel gap in the Coast Range, resulting in rapid cooling of the Sacramento area during late afternoon and nighttime hours.
Sacramento County air-quality monitors registered violations of federal (> 0.12 ppm) and California (> 0.09 ppm) standards for ozone (O3) an annual average of 9 and 46 days, respectively, for the period 1990 through 1992 (ARB 1991, 1992, 1993). Federal and state laws allow no more than one event above the hourly standard (“exceedance”) at any air monitoring site per year, averaged over 3 years. As a result, the Sacramento area has been designated a “nonattainment area” for both federal and state health-based standards for ozone. Ozone is formed in the atmosphere through photochemical reactions involving precursors such as oxides of nitrogen (NOx), hydrocarbon compounds, and oxygen. In 1990 estimated daily emission of NOx and hydrocarbons in Sacramento County was 90.9 and 129.0 t, respectively (SMAQMD 1994).
Particulate matter less than 10 μ in diameter (PM10) is generated by such processes as agricultural tillage, construction and demolition, road or vehicle wear, fuel combustion, and atmospheric photochemical reactions involving hydrocarbons, NOx, and oxides of sulfur. Because PM10 can readily enter respiratory airways, adverse health effects can occur from acute and chronic exposures. Federal and California standards for average daily and annual PM10 concentrations were exceeded in each year from 1990 to 1992. In 1990, estimated daily emission of PM10 in Sacramento County was 119 t per day (SMAQMD 1995).
Sulfur oxides (SOx), such as sulfur dioxide (SO2), are generated by combustion of sulfur-containing fuels, such as coal or fuel oil, and emitted principally from stationary sources such as electric utilities, smelters, and pulp and paper mills. Sulfur dioxide can be converted to sulfate aerosols such as sulfuric acid (H2SO4), a component of acid deposition. Unlike many other urban regions, the Sacramento area has few significant stationary SOx sources. Home heating-oil combustion and heavy-duty truck operation comprise the main SO2 sources for the county. In 1990, Sacramento County SOx emissions were 4.2 t per day (SMAQMD, personal communication 1996) and ambient concentrations did not exceed either federal or state SO2 standards.
Trees absorb gaseous pollutants through leaf stomata and bind pollutant particles onto leaf surfaces. When trees absorb gaseous or intercept particulate pollutants without aid of precipitation, it is called “dry deposition.” Deposited pollutant gases and particles can be chemically altered by plant tissues and may be metabolized or cause foliar injury (Smith 1978, 1981). Particles can be resuspended by turbulence or other mechanical action. Absorbed pollutants can be deposited to the ground surface as litter or leaf fall.
Air pollutant uptake by Sacramento’s urban forest was estimated using canopy cover information from the SUFES project and dry deposition algorithms developed for regional air-quality models (Killus et al. 1984; Wesley 1989; Nowak 1994). Air pollutant deposition to buildings, streets, or other surfaces was neglected. To estimate pollutant uptake by trees, dry deposition model calculations incorporated information about air pollutant concentrations as well as meteorological and urban forest conditions. Hourly meteorological data and air pollutant concentrations, together with canopy cover areas from SUFES, were used as inputs. Model simulations were run for a season corresponding to the period when local deciduous trees are in leaf (March 15 through November 15). Hourly pollutant uptake for ozone (O3), nitrogen dioxide (NO2), sulfur dioxide (SO2), and particulate matter (PM10) were summed to yield total monthly and annual estimates of air pollutant removal by trees in the Sacramento area.
(For a complete description of the study area, sampling units, and land-use classes, see pages 175-177 of this issue.)
Emission-control cost factors (dollars per metric ton) were used to estimate economic value of air pollutant uptake by trees. Metric tons of pollutant removed by trees were multiplied by the control-cost factor for the particular pollutant species (Table 1). Control-cost factors were developed by the local air-quality management district (SMAQMD 1993) to determine cost effectiveness of pollutant emission-control technologies. Control-cost factors reflect what local regulators, industry, and other parties have determined to be society’s willingness to pay for controlling precursor and pollutant emissions. We assigned the control-cost factor for NO2 to the value of ozone uptake by trees because direct control costs for ozone have not been defined (e.g., ozone is not directly emitted into the atmosphere).
Deposition
Dry deposition of an air pollutant to a forest can be expressed as the product of a deposition velocity Vd (cm/sor m/s), a pollutant concentration C (g/m3), the tree canopy coverage for a given land area (m2), and a time step. When the canopy is treated as a single layer, the formulation is sometimes called a “big leaf” model. The deposition velocity is represented as an inverse sum of resistances, Vd = 1/(ra + rb + rc), where ra and rb are aerodynamic resistances, and rc is canopy resistance. Hourly deposition velocities for each pollutant were calculated for the Sacramento area for a base year (1990) growing season (March 15 through November 15) using formulations described in the appendix. Hourly meteorological data for wind speed, solar radiation, and precipitation from a California Department of Water Resources (CDWR) monitoring site in Sacramento County, together with nighttime 3-hourly cloud cover data from a local airport, were used as input data. Hourly pollutant deposition to SubRAD (Sub-Regional Assessment District) canopy areas were then summed to obtain monthly and growing-season annual pollutant deposition.
Pollutant concentrations
Hourly concentrations for NO2, O3, and SO2 (ppm), and 24-hour average PM10 concentrations (μg/m3) were obtained from the U.S. EPA Aerometric Information Retrieval System (AIRS) database for 1990 for 10 monitoring stations throughout Sacramento County. One monitoring station was located in a rural area, whereas others were in urban areas of the county. Seven stations monitored for ozone, 6 for NO2, 7 for PM10, and 2 for SO2. Hourly ozone and NO2 concentrations averaged among monitoring stations in urban areas were used to estimate deposition to the city and suburban sectors. Ozone and NO2 concentrations from the rural monitoring station were used to estimate pollutant deposition to the rural sector. PM10 was monitored episodically (approximately one 24-hour period every 6 days) and on different days throughout the study area. Average monthly concentrations for the urban stations were used to estimate PM10 deposition to the city and suburban sectors, while average monthly concentrations from the single rural site were used to estimate particulate deposition to the rural sector. Sulfur dioxide was monitored in only 2 urban sites. Having no other option, we assumed SO2 to be homogeneously distributed throughout the study area.
Pollutant concentrations in Sacramento in 1990 were typical of concentrations for the period 1989 to 1992. For the growing season, O3 concentration between 1500 and 1600 hours averaged 0.061 ppm, with higher average hourly concentrations occurring in July (0.079 ppm), August (0.771 ppm), and September (0.071 ppm). Lowest average hourly afternoon concentrations occurred in March (0.048 ppm) and November (0.045 ppm). In Sacramento County, the 0.12 ppm National Ambient Air Quality Standard (NAAQS) for ozone was exceeded 15 times between June and September (ARB 1991). The growing-season NO2 concentrations averaged 0.017 ppm in the morning hours and 0.023 ppm in the afternoon. The average hourly concentrations were highest in October (0.030, 0.045 ppm) and November (0.029, 0.044 ppm), occurring from 700 to 900 hours and from 1700 to 1900 hours. Monthly PM10 concentrations averaged 33 μg/m3, while peak concentrations during the growing season occurred in September (39 μg/m3), October (44 μg/m3), and November (60 μg/m3). The growing season average hourly SO2 concentration was low—0.001 ppm.
Results and Discussion
Total estimated pollutant uptake by Sacramento’s urban forest was 1,457 metric tons (t) (1 t = 1,000 kg). Ozone and PM10 had the highest removal rates (Table 1). Peak monthly removal rates for ozone occurred in April and May (83.2 and 83.8 t) and for PM10 in September and October (119.7 and 130.2 t). Average growing-season daily uptake for all pollutants combined was approximately 5.91 per day, with ozone and PM10 at approximately 2.4 and 2.7 t per day. Average daily uptake of NO2 was approximately 0.61 per day, which is approximately 1% of the daily anthropogenic oxides of nitrogen emissions in Sacramento County (assuming NOx is converted to NO2) (ARB 1995). Daily removal of SO2 (0.11 t/day) and particulate matter by trees represented approximately 3% and 2% of the daily anthropogenic SO2 and PM10 emissions, respectively, for Sacramento County (ARB 1995).
Because daytime deposition velocities were greater in April than in August, maximum ozone uptake occurred in April, even though ozone concentrations were greater in August. Moderate daytime air temperatures in April (~23°C [~73°F] compared to 32°C [90°F] in August) resulted in bulk canopy stomatal resistances (rs) of approximately 2.5 s/cm from 1200 to 1500 hours, whereas the August daytime rs averaged 20.3 s/cm. Combined with other resistance components, bulk canopy resistance (rc) for April averaged 157.8 s/cm from 1200 to 1500 hours, half that of the August rc (248.3 s/cm). Resulting deposition velocities for ozone during the hours 1200 to 1500 ranged from 0.60 to 0.54 cm/s in April (averaging 0.57 cm/s), a factor of 1.27 greater than deposition velocities in August (0.48 to 0.43 cm/s). For the same daytime period, April ozone concentrations (~0.049 ppm) were 16% lower than August ozone concentrations (~0.058 ppm). Taken together, these two factors account for much of the higher April uptake, compared to August. For example, ozone deposition to the city sector in April was approximately 15.91, whereas deposition in August was 14.7 t (14.7 t × 1.27 × 0.84 ≈ 15.7 t).
Ozone and NO2 concentrations at the rural monitoring station were, on average, 19% and 29% lower, respectively, than concentrations in urban areas. As a result, deposition to the rural sector was slightly less than deposition to combined city and suburban sectors (Table 1), even though total canopy cover area was greater in the rural sector than combined city and suburban sectors (Table 2).
Annual PM10 deposition to combined city and suburban sectors (3411) was nearly equivalent to deposition for the rural sector (338.1 t) (Table 1). Differences between monthly urban and rural PM10 concentrations used in deposition calculations were small for most of the growing season. However, in March the average hourly PM10 concentration for combined city and suburban sectors was 29 μg/m3 and in October the average was 52 μg/m3. Rural sector concentrations for these 2 months were 21 and 43 μg/m3, respectively. Because PM10 concentrations were roughly equivalent for much of the growing season, deposition to the rural sector could have been greater due to greater total canopy cover. However, lower surface roughness zo (Table A1; see appendix) in the rural area contributed to greater aerodynamic resistance, thereby reducing deposition velocities. For example, the October average hourly deposition velocity for particulates in the city sector between 1100 and 1700 hours was 1.3 cm/s (suburban sector average was 1.2 cm/s), while the rural Vd average was 1.1 cm/s. Hence, the rural sector canopy cover advantage (factor = 1.16) is nearly canceled out by decreased deposition velocities (factor = 1.3/1.1 = 1.18). Slightly greater PM10 deposition to combined city and suburban sectors may then be due to higher March and October PM10 concentrations.
Sulfur dioxide deposition to the rural sector was approximately 12% greater than deposition to urban sectors. Because SO2 was assumed to be homogeneously distributed throughout the study area, the difference in total deposition may be attributed to the greater overall canopy cover area in the rural sector (10,261 ha, Table 2), offset by smaller deposition velocities due to greater aerodynamic resistance for the rural sector. Average hourly deposition velocities between 1100 and 1700 hours in April, when monthly SO2 uptake was greatest, for the city, suburban, and rural sectors were 0.46, 0.44, and 0.43 cm/s, respectively.
Total monetary value of pollutant removal (Table 1) was estimated to be US$28.7 million ($1,500/[ha tree cover per yr]). The highest value was for ozone removal (57% of the total value) followed by PM10 (27%). With an estimated 6 million trees in the Sacramento area (McPherson 1998), the annual air pollutant removal benefit is approximately $5 per tree per year. Distribution of economic benefit by sector follows the distribution of pollutant deposition by canopy cover (Table 2). For example, due to higher pollutant concentrations in urban areas, monetary value of NO2 and ozone uptake is greater in combined city and suburban sectors than in the rural sector. On the other hand, economic values of PM10 uptake in the urban and rural sectors are comparable, while SO2 uptake is slightly greater in the rural sector due to greater canopy cover and assumed homogeneous distribution of SO2.
Pollutant deposition by land use
Over half of the tree cover in city and suburban sectors is located in the residential low-density land-use category (Table 2). Another 27% to 30% of the tree cover is in the institutional and vacant/wild categories (McPherson 1998). Consequently, it is estimated that over 80% of pollutant deposition to trees in city and suburban sectors occurs in these 3 land-use categories. Over half the tree cover in the rural sector is located in the vacant/wild category, with another 25% in agriculture. Similarly, the bulk of pollutant deposition to trees in the rural sector occurs in these land-use categories. It is not surprising that low-density residential areas in city and suburban sectors contain a large fraction of the sector’s canopy cover. It is surprising that natural and unmanaged areas (e.g., vacant/wild) comprise another major portion of “urban” forest cover, nearly equal to that of institutional lands. This is because vacant/wild lands include riparian corridors, which are present along the American River and numerous creeks, streams, sloughs, and canals throughout the urbanized Sacramento area.
Normalizing total annual deposition to trees by land area results in deposition rates (kg/[ha per yr]) by SubRAD, sector, or land use. The deposition rate reflects the distribution of canopy cover across land uses, as well as surface roughness and pollutant concentrations. Annual pollutant deposition rates by SubRAD and sector land use are shown in Figure 1 and Table 3.
Deposition rates are greatest in a transect from residential Sacramento city neighborhoods and along a northeast corridor between highways l-80 and US 50. Estimated deposition rates in the city and suburban sectors are greatest in the residential low-density, residential high-density, and institutional landuse categories (Table 3). Although vacant/wild lands are areally extensive and contain between 13% and 17% of city and suburban canopy cover (Table 2), they are characterized by relatively open “oak savanna” woodlands (with the exception of riparian habitat). As a result, deposition rates for city and suburban sectors’ vacant/wild lands are less than those of the more densely tree-covered residential and institutional land uses.
These estimated annual pollutant deposition rates are in agreement with simulations for other urban areas. For example, pollutant deposition rates for the Chicago, Illinois, area were estimated to range between 9.7 and 19.4 kg/(ha land area per yr) (Nowak 1994). For the Sacramento study area, deposition rates averaged 10.9 kg/(ha land area per yr). Estimated annual ozone deposition rate per ha of tree cover for the Sacramento area (40.1 kg/[ha tree cover per yr]) was greater than the rate estimated for the Chicago study area (30.7 kg/[ha tree cover per yr]) by a factor of 1.3. Longer local foliation period and greater pollutant concentrations may contribute to the higher estimate for Sacramento. The growing season for Sacramento (246 days) exceeds Chicago’s (May 1 through October 31, 184 days) by a factor of 1.3. Summer daytime peak ozone concentrations for Sacramento averaged 0.061 ppm, a factor 1.4 greater than the Chicago average (~0.043 ppm). Taken together, pollutant uptake rates for the Sacramento area might be expected to exceed rates for the Chicago area by a factor of 1.8 (1.3 × 1.4). However, the hourly-varying daytime canopy resistances of the present model were somewhat greater than the constant values used in the Chicago study. Therefore, resulting deposition velocities with this model may be less than those computed using constant rc values. Extended foliation period and higher pollutant concentrations could presumably result in greater deposition rates for urban forests in other regions, such as the Los Angeles, California, basin.
Due to differences in total canopy cover area, estimated total annual pollutant uptake for Sacramento (1,457 t) is approximately 25% of that estimated for the Chicago study area (5,575 t). Canopy cover for the Chicago study area was estimated at 19% of the 3,350 km2 land area (e.g., 65,000 ha), which is approximately a factor 3.4 greater than the estimated canopy cover for the Sacramento area (19,000 ha). Although canopy cover and estimated total pollutant uptake are greater for the Chicago study area, the estimated monetary value of air pollutant uptake by trees in the Sacramento area (US$28.7 million) is 3 times greater than Chicago’s due to differences in control-cost pricing for the individual pollutants. For example, PM10 uptake by trees in the Chicago and suburban Cook County area (2,479 km2 [953 mi2], similar to the SUFES study area) was estimated at 1,391 t for a value of $2 million ($1,441/t). By contrast, total PM10 uptake by trees in the SUFES study area (679.1 t) is approximately half that for the Chicago and suburban Cook County area, but local control cost ($11,681 t) is 8.1 times greater. Taken together, the value of PM10 uptake ($7.9 million) by Sacramento trees is 4 times greater than the estimate for the Chicago and suburban Cook County area.
Sources of uncertainty and model limitations
Pollutant deposition estimates are sensitive to uncertainties associated with canopy resistance, particle resuspension rates, and inhomogeneities in the spatial distribution of air pollutants. Our use of canopy resistance for deciduous forests may overpredict daytime canopy resistance for well-watered urban trees. In the present model, bulk canopy stomatal resistance (rs) increases with temperature. However, leaf stomata of urban trees with abundant water sources may remain open, in contrast to natural forest stands, which may become water stressed. As a result, deposition to urban forests during warm periods might be underestimated. On the other hand, Massman et al. (1994) found that the Wesley algorithms underpredicted ozone canopy resistance for California Central Valley crops by a factor of 1.4 to 2, which suggests possible overestimation of ozone deposition. An additional source of uncertainty associated with canopy resistance is local climate effects. In the present model, we have employed hourly meteorological data from a single station, which may not be representative of conditions in local atmospheric surface layers.
Other sources of uncertainty are canopy resistance parameter inputs based upon land use (Table A2; see appendix). By using parameters for deciduous forests, we have assumed that the urban forest has a closed canopy and is horizontally homogeneous at scales of several kilometers. Alternatively, urban forests may be better characterized as well-watered “sparse canopies,” comprising patches or stands of trees, impervious surfaces (streets and buildings), turf, and dispersed, open-grown individual trees. Bulk canopy or surface resistances for urban forests are currently unavailable.
Factors that influence particulate dry deposition include atmospheric characteristics (e.g., turbulence), surface properties (e.g., canopy architecture, roughness, albedo, wetness, chemical reactivity), and properties of the depositing particle species (size, mass, chemical composition). Particle resuspension is also influenced by these factors (Davidson and Wu 1990). Based on limited literature for open-grown urban tree canopies (Dochinger 1980; Nowak 1994) we assumed a 50% resuspension rate as a base case. Resuspension rates may range from 20% (resulting in an underestimate by a factor of 1.6 from the base case) to 80% (causing an overestimate by a factor of 2.5 from the base case). In this model, we have assumed that PM10 deposition in the Sacramento area is comprised of relatively large particles, characteristic of fugitive dust from roads, construction, and agricultural activities. Smaller secondary aerosols (PM2.5) are produced by photochemical reactions involving NOx, SOx, and hydrocarbons and comprise a significant portion of the overall mass of particulate pollution in urban areas. So-called process-oriented models, which treat particle removal processes as a function of particle size, will be required to model deposition of both PM10 and PM2.5 to urban forests, rather than the bulk-resistance-oriented model employed in this work (Ruijgrok et al. 1995).
Another source of uncertainty is the inhomogeneity of the pollutant field. Pollutant concentrations are influenced by such factors as prevailing winds and transport from nearby pollutant sources. Monitoring stations are generally located in populated areas to detect elevated concentrations or violations of ambient air-quality standards (ARB 1994). Extrapolation of monitoring station concentrations to distant canopy areas may then result in over estimations of pollutant dry deposition.
In this analysis, we applied the NO2 control cost to the value of ozone uptake by the urban forest. This method is problematic because it treats ozone as a primary pollutant and neglects photochemical reactions involved in the formation of ozone. In addition, biogenic hydrocarbon emissions from trees (BVOCs), not considered here, play a role in ozone formation (Corchnoy et al. 1992; Winer et al. 1995; Benjamin et al. 1996; Guenther et al. 1996). An alternative valuation scheme could involve estimating urban forest BVOC emissions and adjusting associated BVOC emission cost with a correction factor representing a ratio of ozone uptake to ozone produced as a result of BVOC emissions. The interaction of urban forests with climate and photochemical air pollution has been investigated for the Atlanta, Georgia, and Los Angeles, California, metropolitan areas using photochemical and mesoscale meteorological models (Cardelino and Chameides 1990; Rosenfeld et al. 1995; Sailor 1995; Taha 1996). These efforts will benefit from improved algorithms for the representation of dry deposition to urban vegetation and the emission of BVOCs. Modeled BVOC emissions will also be influenced by a better understanding of urban forest canopy structure and the distribution of foliar biomass. By altering urban heat islands, urban forests may affect wind fields and convergence/divergence zones and thereby the distribution of both air pollutants and precipitation. In turn, by intercepting precipitation, urban forest canopies may alter stormwater runoff and water quality, as deposited pollutants drip or flow to the ground, in ways not explicitly treated in current urban hydrologic models (Xiao et al. 1998). Urban forests may then provide benefits for nonpoint-source water pollution and urban stormwater management.
Summary
A dry deposition model has been applied to estimate air pollutant deposition to Sacramento’s urban forest. Assuming 1990 pollutant concentrations, approximately 1,4571 of air pollutants were absorbed by trees throughout the study area, at an implied value of US$28.7 million. Pollutant uptake by trees represented 1% to 2% of anthropogenic emissions for such pollutants as NOx and PM10. Urban land uses with the highest rates of pollutant uptake included residential areas, institutions (e.g., parks, campuses), and vacant, unmanaged, or natural areas. Furthermore, the institutional, commercial/industrial, and vacant/wild land uses for the city and suburban sectors also have high potential for additional tree planting (McPherson 1998). Unlike urban areas in the midwest or eastern United States, canopy cover in the Sacramento area decreases along an urban-to-rural gradient. Therefore, estimated pollutant uptake rates per unit land area were highest for residential and institutional land uses, compared to natural or unmanaged lands. Possible management implications of these estimates are that air pollutant uptake benefits from tree planting may be optimized by planting in areas where air pollutant concentrations are elevated and where relatively high planting densities can be achieved. Trees are an integral component of the complex web of roads, buildings, and waterways that comprise the urban ecosystem. Knowledge about their role in urban ecosystems will help to inform ecosystem management and enhance the health of urban dwellers.
Acknowledgments
The authors thank Coe Owen (US Environmental Protection Agency, Region 9), Mark Ott (Sacramento Metropolitan Air Quality Management District), and Jim Goodridge (California Department of Water Resources), who provided air quality and meteorological data. Thanks also to William Massman (USDA-Forest Service) and Bruce Jackson (California Air Resources Board), who reviewed earlier versions of this manuscript.
Appendix: Summary of the Dry Deposition Model
Dry deposition is defined as the product of a deposition velocity Vd (cm/sor m/s), a pollutant concentration C (g/m3), canopy area (m2), and a time step (Killus et al. 1984; Nowak 1994). By analogy with Ohm’s Law, deposition velocity is defined as the inverse sum of a series of resistances, Vd = 1/(ra + rb + rc), where ra is aerodynamic resistance, rb is quasi-laminar boundary-layer resistance, and rc is the canopy resistance (Baldocchi et al. 1987). Aerodynamic resistance, ra, is expressed as where u(z) is the wind speed (m/s) at height z (defined as 10 m above a zero-plane displacement length, e.g., z–d= 10), and u* is the frictional velocity (m/s). Resistance rb is expressed as where The frictional velocity u*(m/s) is where Ψm is the dimensionless stability function for momentum, k is von Karman’s constant (0.4), d is a zeroplane displacement length, zo is an aerodynamic roughness length for the 3 sectors (city, suburban, rural), and L is the Monin-Obukhov length (m) (van Ulden and Holtslag 1985). L was estimated by classifying hourly meteorological data (solar radiation, wind speed, nighttime cloud cover) using an automated look-up table to estimate 1/L as a function of Pasquill classes (Golder 1972; Gifford 1976). Nine stability classes were used. When atmospheric conditions are unstable (L ≤ 0), the function Ψm (van Ulden and Holtslag 1985) is where the dimensionless X (Dyer and Bradley 1982) is When conditions are stable (L > 0), the function Ψm (van Ulden and Holtslag 1985) is
Surface morphology from SUFES was used to estimate surface roughness (zo) for each sector by land use. Land use was determined through aerial photo interpretation and ground sampling. Tree dimensions were determined from SUFES measurements of trees within 675 randomly distributed 100 m2 plots (McPherson 1998). The formulation of Lettau (1969) was used to estimate surface roughness due to trees: where c2 is a constant (≈0.5), h is an average height of roughness elements (trees) in the area of interest (e.g., sector land use), and λ is a roughness density, defined as a ratio of average silhouette area of obstacles in the area of interest to average lot area. Average lot area is defined as the quotient of total area over which the elements are being considered (e.g., sector land use), divided by the number of elements within the area of interest. Tree height and crown silhouette area were determined from measurements taken of trees located in the plots. Tree silhouette area was calculated using formulas for cylindrical, paraboloid, and ellipsoid crown geometries. Building roughness values for Sacramento land-use categories were adapted from Myrup and Morgan (1972). Sector composite roughness values were derived by selecting the greatest of 2 values of roughness length zo (buildings or trees) and weighting the selected value by the corresponding land use area percentage (Table A1).
Canopy or surface resistance (rc) was parameterized following formulations of Wesley (1989) as
Bulk canopy stomatal resistance rs is represented as where ri is a minimum stomatal resistance for water vapor, G is solar radiation (W/m2), and Ts is the air temperature near the surface (°C) between 0°C and 40°C. Bulk canopy stomatal resistance is set to a large value when temperatures are outside this range, based on the assumption that stomatal uptake is reduced. At nighttime, bulk canopy stomatal resistance is also set to a large value.
A combined minimum stomatal and mesophyll resistance for a pollutant gas x is represented as where Dx is the molecular diffusivity of pollutant gas x in air, DH2O is the molecular diffusivity of water vapor, and rmx is tine mesophyll resistance for the gas of interest.
The mesophyll resistance for gas x is computed as where the parameter H* is an effective Henry’s Law constant, and f0 is a reactivity factor. The 2 parameters represent parallel diffusion pathways to extracellular water within leaf stomata and are used to scale the uptake rates. Parameters H* and f0 are also used in expressions for the resistances rlux, rclx, and rgsx. The resistance rlux represents outer surfaces of the upper canopy. Deposition pathways to leaves and stems in the subcanopy are represented by rclx. Deposition to the “ground” surface (e.g., leaf litter, soil) is represented by rac and rgsx. Resistances rlux rclx and rgsx are computed as where rclO, rclS, rgsS, and rgsS are reference resistances corresponding to those for O3 and SO2.
Buoyant convection and slope effects in plant canopies are represented by the resistance rdc given by where θ (radians) is the slope of the local terrain. In the present model, we let θ = 1.66 × 10–3, representing a transect extending from a National Weather Service observing station just south of the downtown area to an air pollutant monitoring station 28.2 km (10.8 mi) to the northeast. Input parameters and resistances are listed in Table A2.
Canopy resistance for particles (PM10) was estimated based upon an average deposition velocity minus an average ra and rb. In a field experiment of particle interception by open-grown tree stands, Dochinger (1980) reported resuspension rates ranging from 20% to 80%. For trees in Sacramento, a 50% resuspension rate was chosen as a midvalue. Particle interception by evergreen trees and deciduous tree stem surfaces during fall and winter months was neglected. Hourly deposition velocities for each pollutant were calculated using estimates for ra, rb, and rc. Deposition velocities were set equal to zero during periods of precipitation (the CDWR monitoring station registered approximately 137.2 mm (5.4 in.) of precipitation during the growing season).
- © 1998, International Society of Arboriculture. All rights reserved.