Determinants of airborne benzene concentrations in rural areas of western Canada

This study estimated the level and determinants of airborne benzene concentrations in rural western Canada. A multi-site, multi-month unbalanced two-factorial design was used to collect air samples at 1206 fixed sites across a geographic area associated with primary oil and gas industry in Canadian provinces of Alberta, north-eastern British Columbia, and central and southern Saskatchewan from April 2001 to December 2002. Benzene concentrations integrated over 1 calendar month were determined using passive organic vapour monitors. Linear mixed effects models were applied to identify the determinants of airborne benzene concentrations, in particular the proximity to oil and gas facilities. The observed geometric mean of benzene concentrations was 158 ng m⁻³, with large geometric standard deviation: 4.9. Benzene concentrations showed a seasonal variation with maxima in winter and minima in summer. Emissions from oil well (within 2 km) and compressor influenced monthly airborne benzene concentrations. However, in our study, being located in the general area of a gas plant seems to be the most important in determining monthly airborne benzene concentrations. These findings support the need for investigation of the impact of oil and gas industry on quality of rural air.     

A characterization of solution gas flaring in Alberta

Information reported here is the result of a detailed analysis of data on flared and vented solution gas in the Province of Alberta in 1999. A goal of characterizing these flares was to aid in the improved management of solution gas flaring. In total, 4499 oil and bitumen batteries reported flaring or venting with a combined gas volume of 1.42 billion m3. There was significant site-to-site variation in volumes of gas flared or vented, gas composition, and flare design. Approximately 5% of physical batteries generate 35.7% of the gas flared and vented from oil and bitumen batteries. Therefore, if one were to attempt to mitigate flaring, significant progress could be made by starting with only the largest sites. The monthly variability of gas volumes was considered because high variability could affect implementation of alternative technologies. It was found that slightly more than 40% of the sites were reasonably steady and had monthly deviations of 100% or less from the average flared volume. The variability in monthly volumes was less for the larger batteries. Data from individual well sites show significant variability in the relative concentrations of each of the major species contained in solution gas.     

The recovery of ammonia and hydrogen sulfide from ground-level area sources using dynamic isolation flux chambers: bench-scale studies

Controlled bench-scale laboratory experiments were conducted to evaluate the recovery of ammonia (NH3) and hydrogen sulfide (H2S) from dynamic isolation flux chambers. H2S (80-4000 ppb) and NH3 (5000-40,000 ppb) samples were diffused through the flux chamber to simulate ground level area source emissions while measuring the inlet and outlet flux chamber concentrations simultaneously. Results showed that the recovery of H2S during a 30-min sampling time was almost complete for concentrations >2000 ppb. At the lowest concentration of 80 ppb, 92.55% of the H2S could be recovered during the given sampling period. NH3 emissions exhibited similar behavior between concentrations of 5000-40,000 ppb. Within the 30-min sampling period, 92.62% of the 5000-ppb NH3 sample could be recovered. Complete recovery was achieved for concentrations >40,000 ppb. Predictive equations were developed for gas adsorption. From these equations, the maximum difference between chamber inlet and outlet concentrations of NH3 or H2S was predicted to be 7.5% at the lowest concentration used for either gas. In the calculation of emission factors for NH3 and H2S, no adsorption correction factor is recommended for concentrations >37,500 ppb and 2100 ppb for NH3 and H2S, respectively. The reported differences in outlet and inlet concentration above these ranges are outside the fullscale sensitivity of the gas sensing equipment. The use of 46-90 m of Teflon tubing with the flux chambers has apparently no effect on gas adsorption, because recovery was completed almost instantaneously at the beginning of the tests.     

A Characterization of Solution Gas Flaring in Alberta

ABSTRACT

Information reported here is the result of a detailed analysis of data on flared and vented solution gas in the Province of Alberta in 1999. A goal of characterizing these flares was to aid in the improved management of solution gas flaring. In total, 4499 oil and bitumen batteries reported flaring or venting with a combined gas volume of 1.42 billion m³. There was significant site-to-site variation in volumes of gas flared or vented, gas composition, and flare design.

Approximately 5% of physical batteries generate 35.7% of the gas flared and vented from oil and bitumen batteries. Therefore, if one were to attempt to mitigate flaring, significant progress could be made by starting with only the largest sites. The monthly variability of gas volumes was considered because high variability could affect implementation of alternative technologies. It was found that slightly more than 40% of the sites were reasonably steady and had monthly deviations of 100% or less from the average flared volume. The variability in monthly volumes was less for the larger batteries. Data from individual well sites show significant variability in the relative concentrations of each of the major species contained in solution gas.     

Cumulative impact of 40 years of industrial sulfur emissions on a forest soil in west-central Alberta (Canada)

The impact of 40 years of sulfur (S) emissions from a sour gas processing plant in Alberta (Canada) on soil development, soil S pools, soil acidification, and stand nutrition at a pine (Pinus contorta x Pinus banksiana) ecosystem was assessed by comparing ecologically analogous areas subjected to different S deposition levels. Sulfur isotope ratios showed that most deposited S was derived from the sour gas processing plant. The soil subjected to the highest S deposition contained 25.9 kmol S ha(-1) (uppermost 60 cm) compared to 12.5 kmol S ha(-1) or less at the analogues receiving low S deposition. The increase in soil S pools was caused by accumulation of organic S in the forest floor and accumulation of inorganic sulfate in the mineral soil. High S inputs resulted in topsoil acidification, depletion of exchangeable soil Ca2+ and Mg2+ pools by 50%, podzolization, and deterioration of N nutrition of the pine trees.     

Measurements of sulfur dioxide and formaldehyde in Taipei using a differential optical absorption spectrometer

Taipei, the capital city of Taiwan, lies in a basin, and its topography prevents the dispersion of pollutants in the city. As a continuation of our air quality study, from February 1999 through June 1999, we measured the concentrations of SO2 at six different locations and of formaldehyde at five locations using a differential optical absorption spectrometer (DOAS). The average concentration of SO2 varied from 3.5 to 6.6 ppb. The average concentration was highest at Toucheng because of its proximity to point sources. The level in Hsientien was close to that in Toucheng, with Hsinyu showing the lowest concentrations. The DOAS and the Taiwan Air Quality Monitoring Network (TAQMN) measurements for SO2 were highly correlated (r2 > 0.9) for Toucheng, Panchiao, and Hsientien. However, DOAS SO2 concentrations were 2 times higher for Hsientien and slightly lower for Panchiao than the TAQMN concentrations were. The average concentration of formaldehyde varied from 7 to 10 ppb. Diurnal variation of formaldehyde closely followed the variation of ozone, especially when the 1-hr peak ozone concentration was > 60 ppb. Photochemical formation accounted for the ambient levels of formaldehyde in Taipei. Concentration of formaldehyde became significant on days when O3 concentration was high. Our results indicate that DOAS can replace conventional measurement techniques.     

An analysis of flaring and venting activity in the Alberta upstream oil and gas industry

Alberta, Canada, is an important global producer of petroleum resources. In association with this production, large amounts of gas (1.14 billion m3 in 2008) are flared or vented. Although the amount of flaring and venting has been measurably reduced since 2002, data from 2005 reveal sharp increases in venting, which have important implications in terms of resource conservation and greenhouse gas emissions (which exceeded 8 million tonnes of carbon dioxide equivalent in 2008). With use of extensive monthly production data for 18,203 active batteries spanning the years 2002-2008 obtained in close cooperation with the Alberta Energy Resources Conservation Board, a detailed analysis has been completed to examine activity patterns of flaring and venting and reasons behind these trends in the Alberta upstream oil and gas industry. In any given year, approximately 6000 batteries reported flaring and/or venting, but the distribution of volumes flared and vented at individual sites was highly skewed, such that small numbers of sites handled large fractions of the total gas flaring and venting in the Province. Examination of month-to-month volume variability at individual sites, cast in terms of a nominal turndown ratio that would be required for a compressor to capture that gas and direct it into a pipeline, further revealed that volumes at a majority of sites were reasonably stable and there was no evidence that larger or more stable sites had been preferentially reduced, leaving potential barriers to future mitigation. Through linking of geospatial data with production data coupled with additional statistical analysis, the 31.2% increase in venting volumes since 2005 was revealed to be predominantly associated with increased production of heavier oils and bitumen in the Lloydminster region of the Province. Overall, the data suggest that quite significant reductions in flaring and venting could be realized by seeking mitigation solutions for only the largest batteries in the Province.     

Long Term Monitoring of Hydrocarbon Contamination Using Multi-Level Vapor Phase Piezometers

This study evaluated the feasibility of supplementing groundwater monitoring protocols by assessing the vadose zone for the extent of residual subsurface contamination. The study also characterized the response of the soil gas signatures with respect to different soil types and degrees of contamination. A field study was conducted at a former gasoline vending station located in Ottawa, Canada. The current state of contamination was determined by analysis of soil samples taken from boreholes. A series of 10 nested soil gas wells with monitoring depths of 0.75, 1.5, 2.25 and 3.0 m were then installed. Using these wells, soil gas surveys were performed at regular intervals over an extended period to quantify Gaseous TPH (TPH g ), oxygen and carbon dioxide concentrations in the soil gas. Results indicate that soil gas wells located near the source term exhibited characteristic soil gas signatures and significant fluctuations in TPH g , oxygen, and carbon dioxide concentrations with time. Soil gas wells located beyond the soil contamination demonstrated limited correlation between TPH g , oxygen and carbon dioxide concentrations and decreased seasonal variability.     

Predictions of Maximum Ground-Level H2S Concentrations Resulting from Two Sour Gas Well Blowouts

Alberta has recently experienced two sour gas well blowouts: Lodgepole and Claresholm. Sulphur emissions associated with the blowouts were about 1400 and 2 tonnes per day, respectively. The Lodgepole blowout was not only of much greater magnitude but also lasted significantly longer than the Claresholm blowout (67 vs. 4 days). Special air quality monitoring with respect to H₂S was conducted to assess impacts of the blowouts. Monitoring was especially extensive for the Lodgepole incident. Maximum observed ground-level H₂S concentrations were compared to predictions obtained using a Gaussian model which makes allowance for the effects of sonic exit velocity on plume spread and the effects of wind shear on plume transport. There was appreciable agreement between predicted and observed values.     

Long Term Monitoring of Hydrocarbon Contamination Using Multi-Level Vapor Phase Piezometers

This study evaluated the feasibility of supplementing groundwater monitoring protocols by assessing the vadose zone for the extent of residual subsurface contamination. The study also characterized the response of the soil gas signatures with respect to different soil types and degrees of contamination.A field study was conducted at a former gasoline vending station located in Ottawa, Canada. The current state of contamination was determined by analysis of soil samples taken from boreholes. A series of 10 nested soil gas wells with monitoring depths of 0.75, 1.5, 2.25 and 3.0 m were then installed. Using these wells, soil gas surveys were performed at regular intervals over an extended period to quantify Gaseous TPH (TPH_g), oxygen and carbon dioxide concentrations in the soil gas.Results indicate that soil gas wells located near the source term exhibited characteristic soil gas signatures and significant fluctuations in TPH_g, oxygen, and carbon dioxide concentrations with time. Soil gas wells located beyond the soil contamination demonstrated limited correlation between TPH_g, oxygen and carbon dioxide concentrations and decreased seasonal variability.     

Recharge processes drive sulfate reduction in an alluvial aquifer contaminated with landfill leachate

Natural attenuation of contaminants in groundwater depends on an adequate supply of electron acceptors to stimulate biodegradation. In an alluvial aquifer contaminated with leachate from an unlined municipal landfill, the mechanism of recharge infiltration was investigated as a source of electron acceptors. Water samples were collected monthly at closely spaced intervals in the top 2 m of the saturated zone from a leachate-contaminated well and an uncontaminated well, and analyzed for delta(18)O, delta(2)H, non-volatile dissolved organic carbon (NVDOC), SO(4)(2-), NO(3)(-) and Cl(-). Monthly recharge amounts were quantified using the offset of the delta(18)O or delta(2)H from the local meteoric water line as a parameter to distinguish water types, as evaporation and methanogenesis caused isotopic enrichment in waters from different sources. Presence of dissolved SO(4)(2-) in the top 1 to 2 m of the saturated zone was associated with recharge; SO(4)(2-) averaged 2.2 mM, with maximum concentrations of 15 mM. Nitrate was observed near the water table at the contaminated site at concentrations up to 4.6 mM. Temporal monitoring of delta(2)H and SO(4)(2-) showed that vertical transport of recharge carried SO(4)(2-) to depths up to 1.75 m below the water table, supplying an additional electron acceptor to the predominantly methanogenic leachate plume. Measurements of delta(34)S in SO(4)(2-) indicated both SO(4)(2-) reduction and sulfide oxidation were occurring in the aquifer. Depth-integrated net SO(4)(2-) reduction rates, calculated using the natural Cl(-) gradient as a conservative tracer, ranged from 7.5x10(-3) to 0.61 mM.d(-1) (over various depth intervals from 0.45 to 1.75 m). Sulfate reduction occurred at both the contaminated and uncontaminated sites; however, median SO(4)(2-) reduction rates were higher at the contaminated site. Although estimated SO(4)(2-) reduction rates are relatively high, significant decreases in NVDOC were not observed at the contaminated site. Organic compounds more labile than the leachate NVDOC may be present in the root zone, and SO(4)(2-) reduction may be coupled to methane oxidation. The results show that sulfur (and possibly nitrogen) redox processes within the top 2 m of the aquifer are directly related to recharge timing and seasonal water level changes in the aquifer. The results suggest that SO(4)(2-) reduction associated with the infiltration of recharge may be a significant factor affecting natural attenuation of contaminants in alluvial aquifers.     

Alfalfa nutritive quality for ruminant livestock as influenced by ambient air quality in west-central Alberta

Alfalfa (Medicago sativa) nutritive quality response to ambient ozone (O(3)), sulfur dioxide (SO(2)) and oxides of nitrogen (NO(x)) were assessed at three locations in west-central Alberta, Canada (1998-2002). Yield data were segregated into high and low relative to overall median yield. Ozone concentrations (hourly median and 95th-percentile) and precipitation (P) contributed 69 and 29%, respectively, to the variability in crude protein (CP) concentration in low-yielding alfalfa, whereas mean temperature (T) and relative humidity (RH) collectively influenced 98% of the variation in CP in high-yielding alfalfa. Three-fourths of the accounted variation in relative feed value (RFV) of low-yielding alfalfa was attributable to P, T and RH, whereas median and 95th-percentile hourly O(3) concentrations and SO(2) and NO(x) exposure integrals contributed 25%. In contrast, air quality, (mainly O(3)) influenced 86% of the accounted variation in RFV of high-yielding alfalfa, and T and P collectively contributed 14%.     

The Recovery of Ammonia and Hydrogen Sulflde from Ground-Level Area Sources Using Dynamic Isolation Flux Chambers: Bench-Scale Studies

Abstract

Controlled bench-scale laboratory experiments were conducted to evaluate the recovery of ammonia (NH₃) and hydrogen sulflde (H₂S) from dynamic isolation flux chambers. H₂S (80–4000 ppb) and NH₃ (5000–40,000 ppb) samples were diffused through the flux chamber to simulate ground level area source emissions while measuring the inlet and outlet flux chamber concentrations simultaneously. Results showed that the recovery of H₂S during a 30-min sampling time was almost complete for concentrations >2000 ppb. At the lowest concentration of 80 ppb, 92.55% of the H₂S could be recovered during the given sampling period. NH₃ emissions exhibited similar behavior between concentrations of 5000–40,000 ppb. Within the 30-min sampling period, 92.62% of the 5000-ppb NH₃ sample could be recovered. Complete recovery was achieved for concentrations >40,000 ppb. Predictive equations were developed for gas adsorption. From these equations, the maximum difference between chamber inlet and outlet concentrations of NH₃ or H₂S was predicted to be 7.5% at the lowest concentration used for either gas. In the calculation of emission factors for NH₃ and H₂S, no adsorption correction factor is recommended for concentrations >37,500 ppb and 2100 ppb for NH₃ and H₂S, respectively. The reported differences in outlet and inlet concentration above these ranges are outside the full-scale sensitivity of the gas sensing equipment. The use of 46–90 m of Teflon tubing with the flux chambers has apparently no effect on gas adsorption, because recovery was completed almost instantaneously at the beginning of the tests.     

Foliar sulphur species in pine: a new indicator of a forest ecosystem under air pollution stress

An eleven-year foliar sulphur (S) monitoring program was carried out from 1976 to 1986 near a sulphur recovery-gas plant in west-central Alberta, Canada, as part of a case study designed to determine the effects of chronic, low concentration sulphur gas emissions on the forest ecosystem surrounding the gas plant. Measurements of both foliar total sulphur (ST) and foliar inorganic sulphur (SO4-S) concentration in lodgepole x jack pine trees at the end of each of the 11 growing seasons were taken to provide an indication of S loading of the forest from industrial sulphur emissions. To measure the state of the forest ecosystem, foliar ST was separated into foliar accumulated sulphur (inorganic sulphur or SO4-S) and foliar assimilated sulphur (organic sulphur or S0) and the ratio of SO4-S/S0 taken. Foliar S0 was calculated as the difference between foliar ST and foliar SO4-S. The median SO4-S/S0 ratio, with all three years of needles considered, varied from 0.29 at a reference location (AV) to 0.88 at the location with the highest stress (AI). The corresponding mean values ranged from 0.3 at the reference location to 2.2 at the location of highest stress. The mean seasonal photosynthetic rate of current year's foliage of the pine trees and soil pH were reduced at a stressed location (AI) compared to the reference location (AV), between 1976 and 1981. Over this same time period the mean foliar SO4-S/S0 ratio increased from 0.4 +/- 0.1 to 1.0 +/- 0.3 at the stressed location (AI) and remained nearly the same at the reference location (AV) at 0.3 +/- 0.1. This research suggests that the foliar SO4-S/S0 ratio is a useful indicator of the state of forest ecosystems under S air pollution stress. It is concluded that foliar S separated into various fractions has potential as an early warning environmental management tool.     

Compositions and greenhouse gas emission factors of flared and vented gas in the Western Canadian Sedimentary Basin

A significant obstacle in evaluating mitigation strategies for flaring and venting in the upstream oil and gas industry is the lack of publicly available data on the chemical composition of the gas. This information is required to determine the economic value of the gas, infrastructure and processing requirements, and potential emissions or emissions credits, all of which have significant impact on the economics of such strategies. This paper describes a method for estimating the composition of solution gas being flared and vented at individual facilities, and presents results derived for Alberta, Canada, which sits at the heart of the Western Canadian Sedimentary Basin. Using large amounts of raw data obtained through the Alberta Energy Resources Conservation Board, a relational database was created and specialized queries were developed to link production stream data, raw gas samples, and geography to create production-linked gas composition profiles for approximately half of the currently active facilities. These were used to create composition maps for the entire region, to which the remaining facilities with unknown compositions were geographically linked. The derived data were used to compute a range of solution gas composition profiles and greenhouse gas emission factors, providing new insight into flaring and venting in the region and enabling informed analysis of future management and mitigation strategies.

Implications: Accurate and transparent determination of environmental impacts of flaring and venting of gas associated with oil production, and potential benefits of mitigation, is severely hampered by the lack of publicly available gas composition data. In jurisdictions within the Western Canadian Sedimentary Basin, frameworks exist for regulating and trading carbon offset credits but current potential for mitigation is limited by a lack of standardized methods for calculating CO₂ equivalent emissions. The composition and emission factor data derived in this paper will be useful to industry, regulators, policy researchers, and entrepreneurs seeking statistically significant and openly available data necessary to manage and mitigate upstream flaring and venting activity and estimate greenhouse gas impacts.     

A prediction-based approach to modelling temporal and spatial variability of traffic-related air pollution in Montreal, Canada

Concentrations of traffic-related air pollution can be highly variable at the local scale and can have substantial seasonal variability. This study was designed to provide estimates of intra-urban concentrations of ambient nitrogen dioxide (NO₂) in Montreal, Canada, that would be used subsequently in health studies of chronic diseases and long-term exposures to traffic-related air pollution. We measured concentrations of NO₂ at 133 locations in Montreal with passive diffusion samplers in three seasons during 2005 and 2006. We then used land use regression, a proven statistical prediction method for describing spatial patterns of air pollution, to develop separate estimates of spatial variability across the city by regressing NO₂ against available land-use variables in each of these three periods. We also developed a “pooled” model across these sampling periods to provide an estimate of an annual average. Our modelling strategy was to develop a predictive model that maximized the model R². This strategy is different from other strategies whose goal is to identify causal relationships between predictors and concentrations of NO₂.Observed concentrations of NO₂ ranged from 2.6 ppb to 31.5 ppb, with mean values of 12.6 ppb in December 2005, 14.0 ppb in May 2006, and 8.9 ppb in August 2006. The greatest variability was observed during May. Concentrations of NO₂ were highest downtown and near major highways, and they were lowest in the western part of the city. Our pooled model explained approximately 80% of the variability in concentrations of NO₂. Although there were differences in concentrations of NO₂ between the three sampling periods, we found that the spatial variability did not vary significantly across the three sampling periods and that the pooled model was representative of mean annual spatial patterns.     

Dissociation of sulfur hexafluoride tracer gas in the presence of an indoor combustion source

As an odorless, nontoxic, and inert compound, sulfur hexafluoride (SF6) is one of the most widely used tracer gases in indoor air quality studies in both controlled and uncontrolled environments. This compound may be subject to reactions with water vapor under elevated temperature to form acidic inorganic compounds such as HF and H2SO4. Thus, in the presence of unvented combustion sources such as kerosene heaters, natural gas heaters, gas log fireplaces, candles, and lamps, the SF6 dissociation may interfere with measurements of the emissions from these sources. Tests were conducted in a research house with a vent-free natural gas heater to investigate these potential interferences. It was observed that the heater operation caused about a 5% reduction of SF6 concentration, which can be an error source for the ventilation rate measurement and consequently the estimated pollutant emission rates. Further analysis indicates that this error can be much greater than the observed 5% under certain test conditions because it is a function of the ventilation flow rate. Reducing the tracer gas concentration has no effect on this error. A simple theoretical model is proposed to estimate the magnitude of this error. The second type of interference comes from the primary and secondary products of the SF6 dissociation, mainly H2SO4, SO2, HF, and fine particulate matter (PM). In the presence of approximately 5 ppm SF6, the total airborne concentrations of these species increased by a factor of 4-10. The tests were performed at relatively high SF6 concentrations, which is necessary to determine the interferences quantitatively. The second type of interference can be significantly reduced if the SF6 concentration is kept at a low ppb level.     

Environmental building policy by the use of microalgae and decreasing of risks for Canadian oil sand sector development

Environmental building recommendations aimed towards new environmental policies and management-changing decisions which as example demonstrated in consideration of the problems of Canadian oil sands operators. For the implementation of the circular economic strategy, we use an in-depth analysis of reported environmental after-consequence on all stages of the production process. The study addressed the promotion of innovative solutions for greenhouse gas emission, waste mitigation, and risk of falling in oil prices for operators of oil sands with creating market opportunities. They include the addition of microalgae biomass in tailings ponds for improvement of the microbial balance for the water speedily cleaning, recycling, and reusing with mitigation of GHG emissions. The use of food scraps for the nutrition of microalgae will reduce greenhouse gas emission minimally, on 0.33 MtCO₂eq for Alberta and 2.63 MtCO₂eq/year for Canada. Microalgae-derived biofuel can reduce this emission for Alberta on 11.9–17.9 MtCO₂eq and for Canada on 71–106 MtCO₂eq/year, and the manufacturing of other products will adsorb up to 135.6 MtCO₂ and produce 99.2 MtO₂. The development of the Live Conserve Industry and principal step from non-efficient protection of the environment to its cultivation in a large scale with mitigation of GHG emission and waste as well as generating of O₂ and value-added products by the use of microalgae opens an important shift towards a new design and building of a biological system.     

Urban and rural ozone concentrations in Alberta,Canada

Ozone concentrations in Alberta cities typically exhibit a maximum in May (up to 35 ppb) and a minimum in November (as low as 4 ppb). This behaviour is similar to that of rural Alberta O₃ concentrations. Annual O₃ concentrations at six urban monitoring stations vary from 11 ppb to 22 ppb and are about one-half the values at rural stations. In winter, urban O₃ concentrations are always smaller than rural concentrations and the cities act as sinks for O₃. Although urban stations do not exceed Canada's maximum acceptable levels of daily (25 ppb) and annual (15 ppb) O₃ concentrations as often as rural stations, the frequency is still quite large. Canada's hourly maximum desirable level (50 ppb) is exceeded 11 times more often at the remote (rural) station than at the downtown (urban) stations.     

Foliar injury symptoms of Saskatoon serviceberry (Amelanchier alnifolia Nutt.) as a biological indicator of ambient sulfur dioxide exposures

Saskatoon serviceberry or Saskatoon (Amelanchier alnifolia Nutt. cv. Smoky) seedlings were planted at five study sites within a 35,000 km(2) airshed, that is influenced by a number of isolated stationary sources of sulfur dioxide (SO(2)), oxides of nitrogen and hydrocarbons, among others. The locations of the five sites were based on the results of a meteorological dry deposition model for the oxides of sulfur and nitrogen. Visible foliar injury responses of Saskatoon were used as a biological indicator of SO(2) exposures, through monthly field surveys. During late July 1998, unifacial, interveinal chlorosis was observed on some 12% of the seedlings at one study site. By September, the chlorosis had become more severe (necrosis) on some 70% of the plants at that site. Site specific ambient SO(2) levels were relatively low (maximum 5-min concentration of 52.8 ppb). Similar data were unavailable for all, but one other site. Therefore, foliar total S and SO(4)(2-)-S concentrations were analyzed in September at four of the five study sites. Previously soil SO(4)(2-)-S at these sites had been analyzed. There were spatial variabilities among these parameters. Based on the overall examination of these data, it is concluded that the observed visible injury symptoms were due to chronic SO(2) exposures, exacerbated by the presence of ozone (O(3)). Independent of this literature based speculation, visible foliar injury responses of Saskatoon can be used as a biological indicator for acute or chronic ambient SO(2) exposures, in the presence of other phytotoxic air pollutants.