Fluorescent Tracer Studies of Pollutant Transport in the San Francisco Bay Area

To simulate the transport and diffusion of airborne contaminants across a metropolitan region, point-source releases of fluorescent tracer material were made near various urban centers and some 50 samplers were arrayed in expected downwind directions. The effects of land-water, hill-valley, and urban-rural differences on airflow and diffusion were observed in their existing interrelationships during these experiments. Since the tracer could be assessed with high sensitivity over great distances, tracer results provided a quantitative indicator of pollutant dispersion across an extensive metropolitan complex.

From July 1967 through June 1968, the test series included typical seasonal weather patterns, with emphasis on those conducive to the travel and accumulation of pollutants. In each test about 15 kilograms of tracer material were released during two-hour periods, and significant dosages were found at downwind distances up to 80 kilometers. All tests were conducted during daylight hours, to coincide better with the oxidant-type pollution important in this region.

Dispersion characteristics showed much greater complexity than predictable from classical models, thus limiting the applicability of such models in this region. Built-up urban areas increased the initial dispersion rates of tracer clouds, and travel over water tended to decrease them. Hilly terrain resulted in increased dispersion, but channeling associated with such terrain caused locally higher concentrations. The complex horizontal dosage patterns obtained did confirm previously observed airflow patterns as aids in predicting pollutant distributions.     

Tracer Study In an Urban Valley

Results of a series of nighttime tracer experiments conducted during the Autumn of 1966 in the industrialized valley of Johnstown, Pa., are discussed. Quite atypical meteorology and dispersion occur within a classical drainage flow framework. An urban heat island effect is observed creating uniform temperature and wind structures within a layer of air flowing through the valley. Dispersion in the valley at night is comparable to that of neutral conditions over open country.     

Tracer Study In an Urban Valley

Results of a series of nighttime tracer experiments conducted during the Autumn of 1966 in the industrialized valley of Johnstown, Pa., are discussed. Quite atypical meteorology and dispersion occur within a classical drainage flow framework. An urban heat island effect is observed creating uniform temperature and wind structures within a layer of air flowing through the valley. Dispersion in the valley at night is comparable to that of neutral conditions over open country.     

Experiments on vertical transverse mixing in a large-scale heterogeneous model aquifer

Vertical transverse mixing is known to be a controlling factor in natural attenuation of extended biodegradable plumes originating from continuously emitting sources. We perform conservative and reactive tracer tests in a quasi two-dimensional 14 m long sand box in order to quantify vertical mixing in heterogeneous media. The filling mimics natural sediments including a distribution of different hydro-facies, made of different sand mixtures, and micro-structures within the sand lenses. We quantify the concentration distribution of the conservative tracer by the analysis of digital images taken at steady state during the tracer-dye experiment. Heterogeneity causes plume meandering, leading to distorted concentration profiles. Without knowledge about the velocity distribution, it is not possible to determine meaningful vertical dispersion coefficients from the concentration profiles. Using the stream-line pattern resulting from an inverse model of previous experiments in the sand box, we can correct for the plume meandering. The resulting vertical dispersion coefficient is approximately approximately 4 x 10(-)(9) m(2)/s. We observe no distinct increase in the vertical dispersion coefficient with increasing travel distance, indicating that heterogeneity has hardly any impact on vertical transverse mixing. In the reactive tracer test, we continuously inject an alkaline solution over a certain height into the domain that is occupied otherwise by an acidic solution. The outline of the alkaline plume is visualized by adding a pH indicator into both solutions. From the height and length of the reactive plume, we estimate a transverse dispersion coefficient of approximately 3 x 10(-)(9) m(2)/s. Overall, the vertical transverse dispersion coefficients are less than an order of magnitude larger than pore diffusion coefficients and hardly increase due to heterogeneity. Thus, we conclude for the assessment of natural attenuation that reactive plumes might become very large if they are controlled by vertical dispersive mixing.     

A Microcomputer-based Forecasting Model: Potential Applications for Emergency Response Plans and Air Quality Studies

Most atmospheric transport and diffusion models within emergency response systems have very limited physics and are forced to rely on the assumption that wind and turbulence conditions at the time of the release will be representative over the period for which dispersion must be predicted. For releases where the principal concern is about the first few kilometers of travel, such an assumption is appropriate. However, for large accidental releases during stable conditions, the plume may travel for several hours before it is diluted to safe levels and the assumption of persistence may be inappropriate, particularly for transport in complex terrain. Under these circumstances, a model that can forecast changes in wind and turbulence conditions is required. We have installed such a model on microcomputers and tested it in complex terrain near Salt Lake City. One-hour tracer releases produced surface concentrations that remain high for much longer times than that expected based on one hour’ travel time with the mean wind at the source height. Furthermore, relatively large concentrations were found at distances of over 40 km from the source. The model was generally able to reproduce the principal features described by the measurements, although some effects of subgrid scale terrain were missed.     

Description of pollutant dispersion in an urban street canyon using a two-dimensional lattice model

The pollutant dispersion in a street canyon has been described in this work by using an isothermal two-dimensional lattice model coupled to the Smagorinsky sub-grid scale model. The influence of the ratio between the height of the upstream and downstream canyon walls, as well as the gap distance between them on the flow pattern, was analyzed considering the situations of ‘open country’ or isolated street canyon and ‘urban roughness’ in which the influence of an urban fabric was considered. The model determined the trajectories of a large number of passive tracer particles released in the computational domain, making it easy to visualize the flow regimes established in each case. The results agreed with the observations reported from the experiments showing a strong influence on the flow inside the canyon exerted by the upstream landscape configuration.     

Atmospheric Tracer Techniques and Gas Transport in the Primary Aluminum Industry

The results of 35 Individual SF₆ tracer tests conducted in Norway during 1978 demonstrate the applicability of tracer techniques to the study of a wide variety of pollutant transport problems found in the primary aluminum industry. Tracer methods were employed to determine the efficiency of the pollutant control system over a single reduction cell under a variety of operating conditions. Two tests conducted during normal operation gave efficiencies equal to 100 ±19% and 79 ± 12%, while a test performed during the occurrence of an anode effect yielded an efficiency equal to 66 ± 22%.

Tracer investigations of flow in the wake of a smelter hall indicated that between 1 % and 11 % of secondary, roof-top emissions can become entrained in the recirculation cavity and reenter the hall through the ventilation fresh air supply. These reentry rates were observed for release heights as high as 8 m above the existing roof exhaust duct. Tracer dispersion data collected within 20 building heights of the smelter agreed very well with extrapolations of McEIroy- Pooler dispersion curves for an urban area. Dispersion curves determined from a previous wind tunnel study of flow downwind of an isolated building underestimated dispersion downwind of the vs.melter complex.

The total fluoride mass flow rate measured downwind of a smelter during wet, foggy conditions indicated that wet removal rates of fluorides are in the range 3.2 × 10^?4/s to 6.4 × 10^?4/s. Simulation of the source with several tracer point releases and simultaneous measurement of fluoride and tracer ground-level concentrations downwind of the smelter eliminated the need for measurements of vertical profiles of wind speed and fluoride concentration during the experiment.     

A simple urban dispersion model tested with tracer data from Oklahoma City and Manhattan

A simple urban dispersion model is tested that is based on the Gaussian plume model and modifications to the Briggs urban dispersion curves. An initial dispersion coefficient (σ_o) of 40 m is assumed to apply in built-up downtown areas, and the stability is assumed to be slightly unstable during the day and slightly stable during the night. Observations from tracer experiments during the Joint Urban 2003 (JU2003) field study in Oklahoma City and the Madison Square Garden 2005 (MSG05) field study in Manhattan are used for model testing. The tracer SF₆ was released during JU2003 near ground level in the downtown area and concentrations were observed at over 100 locations within 4 km from the source. Six perfluorocarbon tracer (PFT) gases were released near ground level during MSG05 and sampled by about 20 samplers at the surface and on building roofs. The evaluations compare concentrations normalized by source release rate, C/Q, for each sampler location and each tracer release, where data were used only if both the observed and predicted concentrations exceeded threshold levels. At JU2003, for all samplers and release times, the fractional mean bias (FB) is about 0.2 during the day (20% mean underprediction) and 0.0 during the night. About 45 –50% of the predictions are within a factor of two (FAC2) of the observations day and night at JU2003. The maximum observed C/Q is about two times the maximum predicted C/Q both day and night. At MSG05, for all PFTs, surface samplers, and release times, FB is 0.14 and FAC2 is about 45%. The overall 60 min-averaged maximum C/Q is underpredicted by about 40% for the surface samplers and is overpredicted by about 25% for the building-roof samplers.     

Near-field dispersion modeling for regulatory applications

This paper evaluates the application of dispersion models to estimate near-field pollutant concentrations in two case studies. The Industrial Source Complex Short-Term Model (ISCST3) was evaluated with hexavalent chromium measurements collected within 100 m of two facilities in Barrio Logan, San Diego, CA. ISCST3 provided reasonable estimates for higher pollutant concentrations but underestimated lower concentrations. To understand the observed distribution of concentrations in Barrio Logan, a recently conducted tracer experiment was analyzed. The tracer, sulfur hexafluoride, was released at ambient temperature from an urban facility at the University of California at Riverside, and concentrations were measured within 20 m of the source. Modeling results indicated that Industrial Source Complex-Plume Rise Model Enhancement and American Meteorological Society/U.S. Environmental Protection Agency Regulatory Model-Plume Rise Model Enhancement overestimated high concentrations and underestimated low concentrations. A diagnostic study with a simple Gaussian dispersion model that incorporated site-specific meteorology was used to evaluate model results. This study found that incorporating lateral meandering for nonbuoyant urban plumes in Gaussian dispersion models could improve concentration estimates even when downwash is not considered. Incorporating a meandering component in ISCST3 resulted in improvements in estimating hexavalent chromium concentrations in Barrio Logan. Credible near-source concentration estimates depend on accurate characterization of emissions, onsite micrometeorology, and a method to account for lateral meandering in the near field.     

On the influence of the urban roughness sublayer on turbulence and dispersion

The concept of the urban roughness sublayer is discussed and this lowest atmospheric layer over a rough surface is shown to have a non-negligible vertical extension over typical urban surfaces. The existing knowledge on the turbulence and flow structure within an urban roughness sublayer is reviewed, focusing on the height dependence of turbulent fluxes and a scaling approach for turbulence statistics, such as velocity variances, in the above-roof part of the roughness sublayer. Finally, the implication of this turbulence and flow structure upon dispersion characteristics is investigated. The most prominent difference of explicitly taking into account the roughness sublayer in a dispersion simulation (as compared to assuming a `constant flux layer') is a clearly enhanced ground level concentration far downwind from the source. For the example of a tracer release experiment over a (sub) urban surface (Copenhagen) it is shown that introducing the roughness sublayer clearly improves the model performance.     

Near-Field Dispersion Modeling for Regulatory Applications

Abstract

This paper evaluates the application of dispersion models to estimate near-field pollutant concentrations in two case studies. The Industrial Source Complex Short-Term Model (ISCST3) was evaluated with hexavalent chromium measurements collected within 100 m of two facilities in Barrio Logan, San Diego, CA. ISCST3 provided reasonable estimates for higher pollutant concentrations but underestimated lower concentrations. To understand the observed distribution of concentrations in Barrio Logan, a recently conducted tracer experiment was analyzed. The tracer, sulfur hexafluoride, was released at ambient temperature from an urban facility at the University of California at Riverside, and concentrations were measured within 20 m of the source. Modeling results indicated that Industrial Source Complex–Plume Rise Model Enhancement and American Meteorological Society/U.S. Environmental Protection Agency Regulatory Model–Plume Rise Model Enhancement overestimated high concentrations and underestimated low concentrations. A diagnostic study with a simple Gaussian dispersion model that incorporated site-specific meteorology was used to evaluate model results. This study found that incorporating lateral meandering for nonbuoyant urban plumes in Gaussian dispersion models could improve concentration estimates even when downwash is not considered. Incorporating a meandering component in ISCST3 resulted in improvements in estimating hexavalent chromium concentrations in Barrio Logan. Credible near-source concentration estimates depend on accurate characterization of emissions, onsite micrometeorology, and a method to account for lateral meandering in the near field.     

Analysis of colloid and tracer breakthrough curves

We consider the dispersion and elution of colloids and dissolved nonsorbing tracers within saturated heterogeneous porous media. Since flow path geometry in natural systems is often ill-characterized macroscopic (mean) flow rates and dispersion tensors are utilized in order to account for the sub-model scale microscopic fluctuations in media structure (and the consequent hydrodynamic profile). Even for tracer migration and dispersal this issue is far from settled.Here we consider how colloid and tracer migration phenomena can be treated consistently. Theoretical calculations for model flow geometries yield two quantitative predictions for the transport of free (not yet captured) colloids with reference to a non-sorbing dissolved tracer within the same medium: the average migration velocity of the free colloids is higher than that of the tracer; and that the ratio of the equivalent hydrodynamic dispersion rates of colloids and tracer is dependent only upon properties of the colloids and the porous medium, it is independent of pathlengths and fluid flux, once length scales are large enough.The first of these is well known, since even in simple flow paths free colloids must stay more centre stream. The second, if validated suggests how solute and colloid dispersion may be dealt with consistently in macroscopic migration models. This is crucial since dispersion is usually ill-characterized and unaddressed by the experimental literature. In this paper we present evidence based upon an existing Drigg field injection test for the validity of these predictions.We show that starting from experimental data the fitted dispersion rates of both colloids and non-sorbing tracers increase with the measured elution rates (obeying slightly different rules for tracers and colloids); and that the ratio of colloid and nonsorbing tracer elution rates, and the ratio of colloid and nonsorbing tracer dispersion rates may be dependent upon properties of the colloids and the medium (not the flow regime).It is important to realize that even for unretarded species, an earlier peak in the breakthrough curve does not necessarily correspond to a faster mean elution rate, or vice versa. But rather that a colloid may elute faster but disperse less than an equivalent tracer. Hence its peak may be retarded compared to that of the tracer, even assuming no retardation. Hence one must consider a combination of mean elution rate and mean dispersion rate, and not rely on “peak times” to corroborate chromatographic effects. The importance of this lies in the fact that these processes are not independent and yet upscale differently. Thus realistic estimates of effective colloid dispersion rates should be upscaled in a way consistent with that adopted for tracers within the same system.     

Validation of the operational emergency response model at the Swedish Meteorological and Hydrological Institute using data from etex and the chernobyl accident

The Eulerian atmospheric tracer transport model MATCH (Multiscale Atmospheric Transport and Chemistry model) has been extended with a Lagrangian particle model treating the initial dispersion of pollutants from point sources. The model has been implemented at the Swedish Meteorological and Hydrological Institute in an emergency response system for nuclear accidents and can be activated on short notice to provide forecast concentration and deposition fields.The model has been used to simulate the transport of the inert tracer released during the ETEX experiment and the transport and deposition of ¹³⁷Cs from the Chernobyl accident. Visual inspection of the results as well as statistical analysis shows that the extent, time of arrival and duration of the tracer cloud, is in good agreement with the observations for both cases, with a tendency towards over-prediction for the first ETEX release. For the Chernobyl case the simulated deposition pattern over Scandinavia and over Europe as a whole agrees with observations when observed precipitation is used in the simulation. When model calculated precipitation is used, the quality of the simulation is reduced significantly and the model fails to predict major features of the observed deposition field.     

Evaluation of risks from urban air pollutants in the southeast Chicago area

Region V of the U.S. Environmental Protection Agency has conducted a comprehensive study of cancer risks from urban exposure to air pollutants in the Southeast Chicago area. This study estimated emissions of a list of 30 air carcinogens from a broad range of nontraditional, as well as traditional, source types. Using dispersion modeling and applying the appropriate unit risk factors and population data, this study estimated the risks at each receptor location and the total number of cancer cases attributable to air pollution in the area. This analysis estimated that current concentrations would cause 77 cases of cancer over the next 70 years, an average risk of 2.0 X 10(-4). Contributions from different source types and different pollutants were estimated. The total contribution from nontraditional source types was less than 0.3 percent. Although these estimates are highly uncertain, the study does suggest the nature and general magnitude of cancer risks from air pollution in the urban area studied.     

Incorporating layer- and local-scale heterogeneities in numerical simulation of unsaturated flow and tracer transport

This study characterizes layer- and local-scale heterogeneities in hydraulic parameters (i.e., matrix permeability and porosity) and investigates the relative effect of layer- and local-scale heterogeneities on the uncertainty assessment of unsaturated flow and tracer transport in the unsaturated zone of Yucca Mountain, USA. The layer-scale heterogeneity is specific to hydrogeologic layers with layerwise properties, while the local-scale heterogeneity refers to the spatial variation of hydraulic properties within a layer. A Monte Carlo method is used to estimate mean, variance, and 5th, and 95th percentiles for the quantities of interest (e.g., matrix saturation and normalized cumulative mass arrival). Model simulations of unsaturated flow are evaluated by comparing the simulated and observed matrix saturations. Local-scale heterogeneity is examined by comparing the results of this study with those of the previous study that only considers layer-scale heterogeneity. We find that local-scale heterogeneity significantly increases predictive uncertainty in the percolation fluxes and tracer plumes, whereas the mean predictions are only slightly affected by the local-scale heterogeneity. The mean travel time of the conservative and reactive tracers to the water table in the early stage increases significantly due to the local-scale heterogeneity, while the influence of local-scale heterogeneity on travel time gradually decreases over time. Layer-scale heterogeneity is more important than local-scale heterogeneity for simulating overall tracer travel time, suggesting that it would be more cost-effective to reduce the layer-scale parameter uncertainty in order to reduce predictive uncertainty in tracer transport.     

Nitrate-enhanced bioremediation of BTEX-contaminated groundwater: parameter estimation from natural-gradient tracer experiments

Two natural-gradient pulse tracer tests were conducted in a petroleum-contaminated aquifer to evaluate the potential for benzene, toluene, ethylbenzene, and xylenes (BTEX) biodegradation under enhanced nitrate-reducing conditions. Addition of nitrate resulted in loss of toluene, ethylbenzene, and m,p-xylenes (TEX) after an initial lag period of approximately 9 days. Losses of benzene were not observed over the 60-day monitoring period. Tracer breakthrough curves (BTCs) were analyzed to derive transport and biodegradation parameters, including advective velocities, retardation factors, dispersion coefficients, biodegradation rate constants, and nitrate utilization ratios. Using the parameters derived from the BTC analysis, numerical simulations of one of the tracer experiments were conducted using BIONAPL/3D [Molson, J., BIONAPL/3D User Guide, A 3D Coupled Flow and Multi-Component Reactive transport model. University of Waterloo, Waterloo, Ontario, Canada]. Simulations using the BTC-derived transport and biodegradation parameters successfully reproduced benzene, TEX, and nitrate concentrations measured during the tracer experiment. Comparisons of observed and simulated nitrate concentrations indicate that the mass ratio of nitrate-N utilized to TEX degraded increased over time during the experiment, reaching values many times that expected based on stoichiometry of TEX oxidation coupled to nitrate reduction. Excess nitrate loss is likely due to oxidation of other organics in addition to TEX.     

Effects of biofilm growth on flow and transport through a glass parallel plate fracture

The effects of biofilm growth on flow and solute transport through a sandblasted glass parallel plate fracture was investigated. The fracture was inoculated using soil microorganisms. Glucose, oxygen and other nutrients were supplied to support growth. The biomass initially formed discrete clusters attached to the glass surfaces, but over time formed a continuous biofilm. From dye tracer tests conducted during biofilm growth, it was observed that channels and low-permeability zones dominated transport. The hydraulic conductivity of the fracture showed a sigmoidal decrease with time. The hydraulic conductivity was reduced by a factor of 0.033, from 18 to 0.6 cm/s, corresponding to a 72% decrease in the hydraulic aperture, from 500 to 140 microm. In contrast, the mass balance aperture, determined from fluoride tracer tests, remained relatively constant, indicating that the impact of biomass growth on effective fracture porosity was much less than the effect on hydraulic conductivity. Analyses of pre-biofilm tracer tests revealed that both Taylor dispersion and macrodispersion were influencing transport. During biofilm growth, only macrodispersion was dominant. The macrodispersion coefficient alpha(macro) was found to increase logarithmically with hydraulic conductivity reduction.     

A simple network approach to modelling dispersion among large groups of obstacles

A simple network approach has been developed to simulate the movement of pollutant within urban areas. The model uses estimates of pollutant exchange obtained from velocity measurements in experiments with various regular obstacle arrays. The transfer of tracer material was modelled using concepts of advection along streets, well-mixed flow properties within street segments and exchange velocities (akin to aerodynamic conductances) across side and top facets of the street segments.The results predicted both the centreline concentration and lateral dispersion of the tracer with reasonable accuracy for a range of packing densities and wind directions. The basic model's concentration predictions were accurate to better than a factor of two in all cases for the region from two obstacle rows behind a source located within the array to around eight rows behind, a range of distances that falls into the so-called “neighbourhood-scale” for dispersion problems. The results supported the use of parameterized rates of exchange between regions of flow as being useful for fast, approximate dispersion modelling. It was thought that the effects of re-entrainment of tracer back into the canopy were of significance, but modelling designed to incorporate these effects did not lead to general improvements to the modelling for these steady-state source experiments.The model's limitations were also investigated. Chief amongst these was that it worked poorly among tall buildings where the well-mixed assumption within street segments was inadequate.     

Evaluation of Risks from Urban Air Pollutants in the Southeast Chicago Area

Region V of the U.S. Environmental Protection Agency has conducted a comprehensive study of cancer risks from urban exposure to air pollutants in the Southeast Chicago area. This study estimated emissions of a list of 30 air carcinogens from a broad range of nontraditional, as well as traditional, source types. Using dispersion modeling and applying the appropriate unit risk factors and population data, this study estimated the risks at each receptor location and the total number of cancer cases attributable to air pollution in the area.

This analysis estimated that current concentrations would cause 77 cases of cancer over the next 70 years, an average risk of 2.0 × 10^?4. Contributions from different source types and different pollutants were estimated. The total contribution from nontraditional source types was less than 0.3 percent. Although these estimates are highly uncertain, the study does suggest the nature and general magnitude of cancer risks from air pollution in the urban area studied.     

A stochastic multi-channel model for solute transport--analysis of tracer tests in fractured rock

Some of the basic assumptions of the advection-dispersion model (AD-model) are revisited. This model assumes a continuous mixing along the flowpath similar to Fickian diffusion. This implies that there is a constant dispersion length irrespective of observation distance. This is contrary to most field observations. The properties of an alternative model based on the assumption that individual water packages can retain their identity over long distances are investigated. The latter model is called the multi-channel model (MCh-model). Inherent in the latter model is that if the waters in the different pathways are collected and mixed, the "dispersion length" is proportional to distance. The conditions for when non-mixing between adjacent streams can be assumed are explored. The MCh- and AD-models are found to have very similar residence time distributions (RTD) for Peclet numbers larger than 3. A generalized relation between flowrate and residence time is developed, including the so-called cubic law and constant aperture assumptions. The two models extrapolate very differently when there is strong matrix interaction. The AD-model could severely underestimate the effluent concentration of a tracer pulse and overestimate the residence time. The conditions are explored for when in-filling particles in the fracture will not be equilibrated but will act as if there was seemingly a much larger flow wetted surface (FWS). It is found that for strongly sorbing tracers, relatively small particles can act in this way for systems and conditions that are typical of many tracer tests. The assumption that the tracer residence time found by cautiously injecting a small stream of traced water represents the residence time in the whole fracture is explored. It is found that the traced stream can potentially sample a much larger fraction of the fracture than the ratio between the traced flowrate and the total pumped flowrate. The MCh-model was used to simulate some recent tracer tests in what is assumed to be a single fracture at the Asp? Hard rock laboratory in Sweden. Non-sorbing tracers, HTO and Uranin were used to determine the mean residence time and its variance. Laboratory data on diffusion and sorption properties were used to "predict" the RTD of the sorbing tracers. At least 30 times larger FWS or 1000 times larger diffusion or sorption coefficients would be needed to explain the observed BTCs. Some possible reasons for such behavior are also explored.