Simulation of dispersion of a power plant plume using an adaptive grid algorithm

A new dynamic adaptive grid algorithm has been developed for use in air quality modeling. This algorithm uses a higher order numerical scheme—the piecewise parabolic method (PPM)—for computing advective solution fields; a weight function capable of promoting grid node clustering by moving grid nodes; and a conservative interpolation equation using PPM for redistributing the solution field after movement of grid nodes. Applications of the algorithm to a model problem, in which emissions from a point source disperse through the atmosphere in time, reflect that the algorithm is able to capture not only the regional ozone plume distribution, but also the small-scale plume structure near the source. In contrast, the small-scale plume structure was not captured in the corresponding static grid solution. Performance achieved in model problem simulations indicates that the algorithm has the potential to provide accurate air quality modeling solutions at costs that may be significantly less than those incurred in obtaining equivalent static grid solutions.     

Fugitive emissions opacity determination using the digital opacity compliance system (DOCS)

Maintenance of Department of Defense (DoD) weapon systems, conducting battlefield training exercises as well as meeting military construction and/or demolition schedules, invariably generate fugitive air emissions, many of which are visible. Although there is no codified federal method for quantifying fugitive emissions opacity, many state and local air regulatory agencies have instituted enforceable fugitive emission opacity standards at DoD facilities. The current study focused on comparing the performance of the digital opacity compliance system (DOCS) with U.S. Environment Protection Agency Method 9 (Method 9) certified human observers in quantifying the visible opacity associated with fugitive emissions produced using a commercial fog generator. By systematically repositioning both DOCS cameras and Method 9-certified observers during field testing, differences in method performance as a function of observational locations were documented. At both the 30- and 300-ft off-set distances, opacity levels reported by the DOCS technology and Method 9-certified smoke readers were found to be statistically different at the 99% confidence level. Alternatively, at the 90- and 150-ft off-set distances, results suggested that there was an insignificant difference at the 99% confidence level between the two methods. Comparing the magnitude of the each method's standard deviation suggested that, at the 30-ft off-set distance, the DOCS technology was consistently more precise than Method 9-certified readers regardless of the observer's downwind distance. However, at the 90, 150, and 300-ft off-set distances, method precision seemed to vary as a function of both off-set and downwind distance. The primary factor affecting the consistency in opacity measurements appeared to be the impact of ground-level air turbulence on fog plume dispersion and transport. Field observations demonstrated that localized wind shear played a critical and decisive role in how and to what extent fugitive emissions opacity could be determined, regardless of the method selected.     

Field demonstration of visible opacity photographic systems

The Digital Opacity Compliance System (DOCS) is an innovative method that uses digital imaging technology to quantify visible opacity of stationary sources. DOCS, which has been demonstrated at pilot and full scale as a technically defensible and economically attractive alternative to U.S. Environmental Protection Agency (EPA) Reference Method 9 (Method 9), uses commercial-off-the-shelf (COTS) digital cameras in combination with a user-friendly computer software package to determine opacity. To date, all DOCS field testing has been conducted using two models of digital cameras, notably, Kodak Models DC265 and DC290, both of which are no longer commercially available. To ensure that field-validated digital cameras will be available to future DOCS users, a suite of new digital cameras was evaluated with the opacity determination software including the following: (1) Sony Model Cybershot Model DSC-WI, (2) Nikon Model Coolpix 5200, (3) Fuji Finepix Model E500, and (4) Kodak Model DX6490. Within the opacity range of regulatory interest, that is, 0-40%, the Sony Cybershot Model DSC-WI and Nikon Coolpix Model 5200 digital cameras were found to generate plume photographs of which the DOCS opacity analysis yielded results that were statistically equivalent to the previously field-validated Kodak Model DC290. In contrast, the Fuji Finepix Model E500 generated plume photographs of which the DOCS opacity analysis were, on average, 2.2% less than those generated by the Kodak Model DC290 photographs, a difference that was determined to be statistically significant. Over the same opacity range, photographs taken by the Kodak Model DX6490 yielded DOCS opacity readings that were found to be statistically equivalent to a Method 9-certified transmissometer. Based on the results from the current digital camera validation testing approach, EPA has developed a new camera-based visible opacity measurement method titled "Determination of Visible Emission Opacity from Stationary Sources Using Computer-Based Photographic Analysis Systems." The proposed method is expected to be promulgated after closure of the public comment period.     

Simulation of the evolution of particle size distributions in a vehicle exhaust plume with unconfined dilution by ambient air

Over the past several years, numerous studies have linked ambient concentrations of particulate matter (PM) to adverse health effects, and more recent studies have identified PM size and surface area as important factors in determining the health effects of PM. This study contributes to a better understanding of the evolution of particle size distributions in exhaust plumes with unconfined dilution by ambient air. It combines computational fluid dynamics (CFD) with an aerosol dynamics model to examine the effects of different streamlines in an exhaust plume, ambient particle size distributions, and vehicle and wind speed on the particle size distribution in an exhaust plume. CFD was used to calculate the flow field and gas mixing for unconfined dilution of a vehicle exhaust plume, and the calculated dilution ratios were then used as input to the aerosol dynamics simulation. The results of the study show that vehicle speed affected the particle size distribution of an exhaust plume because increasing vehicle speed caused more rapid dilution and inhibited coagulation. Ambient particle size distributions had an effect on the smaller sized particles (approximately 10 nm range under some conditions) and larger sized particles (>2 microm) of the particle size distribution. The ambient air particle size distribution affects the larger sizes of the exhaust plume because vehicle exhaust typically contains few particles larger than 2 microm. Finally, the location of a streamline in the exhaust plume had little effect on the particle size distribution; the particle size distribution along any streamline at a distance x differed by less than 5% from the particle size distributions along any other streamline at distance x.     

Migration and natural fate of a coal tar creosote plume: 1. Overview and plume development

A volume of sand containing coal tar creosote was emplaced below the water table at CFB Borden to investigate natural attenuation processes for complex biodegradable mixtures. Coal tar creosote is a mixture of more than 200 polycyclic aromatic hydrocarbons, heterocyclic compounds and phenolic compounds. A representative group of seven compounds was selected for detailed study: phenol, m-xylene, naphthalene, phenanthrene, 1-methylnaphthalene, dibenzofuran and carbazole. Movement of groundwater through the source led to the development of a dissolved organic plume, which was studied over a 4-year period. Qualitative plume observations and mass balance calculations indicated two key conclusions: (1) compounds from the same source can display distinctly different patterns of plume development and (2) mass transformation was a major influence on plume behaviour for all observed compounds.     

Measuring visual opacity using digital imaging technology

The U.S. Environmental Protection Agency (EPA) Reference Method 9 (Method 9) is the preferred enforcement approach for verifying facility compliance with federal visible opacity standards. Supporters of Method 9 have cited its flexibility and low cost as important technological and economic advantages of the methodology. The Digital Opacity Compliance System (DOCS), an innovative technology that employs digital imaging technology for quantifying visible opacity, has been proposed as a technically defensible and economically competitive alternative to Method 9. Results from the field application of the DOCS at EPA-approved Method 9 smoke schools located in Ogden, UT, Augusta, GA, and Columbus, OH, demonstrated that, under clear sky conditions, the DOCS consistently met the opacity error rate established under Method 9. Application of hypothesis testing on the smoke school data set confirmed that the DOCS was equivalent to Method 9 under clear sky conditions. Under overcast sky conditions, human observers seemed to be more accurate than the DOCS in measuring opacity. However, within the smoke school environment, human observers routinely employ backgrounds other than sky (e.g., trees, telephone poles, billboards) to quantify opacity on overcast days. Under conditions that compel the use of sky as plume background (e.g., emission stacks having heights above the tree line), the DOCS appears to be a more accurate methodology for quantifying opacity than are human observers.     

Stable plume rise in a shear layer

Solutions are given for plume rise assuming a power-law wind speed profile in a stably stratified layer for point and finite sources with initial vertical momentum and buoyancy. For a constant wind speed, these solutions simplify to the conventional plume rise equations in a stable atmosphere. In a shear layer, the point of maximum rise occurs further downwind and is slightly lower compared with the plume rise with a constant wind speed equal to the wind speed at the top of the stack. If the predictions with shear are compared with predictions for an equivalent average wind speed over the depth of the plume, the plume rise with shear is higher than plume rise with an equivalent average wind speed.     

Biogeochemical characterisation of a coal tar distillate plume

The distillation of acidified coal tars for up to 50 years has given rise to a phenol plume approximately 500 m long, 50 m deep and containing up to 15 g l(-1) dissolved organic carbon (DOC) in the Triassic Sandstones aquifer. A conceptual biogeochemical model based on chemical and microbiological analysis of groundwater samples has been developed as a preliminary to more detailed studies of the controls on natural attenuation. While the development of redox zones and the production of methane and carbon dioxide provide evidence of natural attenuation, it appears that degradation is slow. The existence of sulphate in the plume indicates that this electron acceptor has not been depleted and that consequently methanogenesis is probably limited. Based on a simple estimate of sulphate input concentration, a half-life of about 15 years has been estimated for sulphate reduction. Geochemical modelling predicts that increased alkalinity within the plume has not led to carbonate precipitation, and thus within the limits of accuracy of the measurement, alkalinity may reflect the degree of biodegradation. This implies a loss of around 18% of the DOC over a 30-year period. Despite limited degradation, microbial studies show that there are diverse microbial communities in the aquifer with the potential for both anaerobic and aerobic biodegradation. Microbial activity was found to be greatest at the leading edge of the plume where DOC concentrations are 60 mg l(-1) or less, but activity could still be observed in more contaminated samples even though cells could not be cultured. The study suggests that degradation may be limited by the high phenol concentrations within the core of the plume, but that once diluted by dispersion, natural attenuation may proceed. More detailed studies to confirm these initial findings are identified and form the basis of associated papers.     

The organic geochemistry of a sanitary landfill leachate plume

Leachate from the North Bay municipal landfill has contaminated an unconfined, sandy aquifer throughout the 700 m flow system from the site to a discharge zone at a creek. The major organic contaminants identified are aromatic hydrocarbons, especially substituted benzenes. The high groundwater velocity of about 75 m yr⁻¹ and the low organic sorption properties of the sand have permitted non-transformed contaminants to spread throughout the total flow system. There is considerable temporal and spatial variability in groundwater chemistry.Most of the aqueous organic carbon has a nominal molecular weight of <2000 and the general decrease in the mass of this fraction relative to Cl indicates it is being mineralized significantly during transport. IR spectra indicate a general trend of increased aromaticity and decreased OH content of organic matter along the flow system. The aqueous organic matter has a significant apparent complexing capacity and so it is somewhat surprising that toxic metal concentrations in leachate-impacted groundwaters are low. In the leachate plume, this complexing capacity is taken up by major cations and H⁺ and to a far lesser extent by toxic metals such as Pb, Cd or Zn.Dispersion is clearly responsible for considerable decrease in contaminant concentration along the flow system. Biotransformation under strictly anaerobic conditions has probably caused 1,1,1-trichloroethane and trichloroethylene to be restricted to the immediate vicinity of the landfill. A simple method of comparing the concentrations of pairs of organics at points along the flow system provides relative transformation rates for pairs of organics even with variable inputs from the landfill and dispersive dilution. Relative to ethylbenzene, o-xylene is rapidly lost from this system. O-xylene may be less persistent than m- or p-xylene; a result unexpected from previous studies of these dimethylbenzenes. In the initial, strictly anarobic segment of the flow system 1,2,4-trimethylbenzene and 1,4-dichlorobenzene are equally persistent, but in the final, less anaerobic segment, the former appears to be degraded more rapidly than the latter.Contaminant distributions in aquifers reflect the results of a number of processes integrated in a complex manner and so are difficult to interpret in terms of specific processes. However, they do provide evidence for what processes are most significant in real groundwater systems and they will also provide critical tests of how well laboratory-derived information relates to real groundwater contamination situations.     

Relative velocities of DNAPL and aqueous phase plume migration

Numerical simulation is used to examine the relative velocities of DNAPL and aqueous phase plumes in sandy aquifers where lateral spreading of DNAPL has occurred at the base of the aquifer. The scenario being modeled is one where a permeable aquifer is underlain by a sloping aquitard, which results in lateral migration of the DNAPL down the slope, in addition to lateral migration of an aqueous phase plume subject to a specified hydraulic gradient. A sensitivity analysis is presented to the impacts of both DNAPL properties and geologic properties. The most important chemical properties governing the relative velocities of the DNAPL and the shallow aqueous phase plume are the DNAPL viscosity and the aqueous component soil-water partition coefficient (Kd). The dip of the underlying aquitard was found to be relatively unimportant, at least for the range of values studied. The scenario under consideration can be important in conceptual model development and remedial design, as in certain cases DNAPL could be migrating in areas without the evidence of a well-developed aqueous phase plume. The implication of this work is that the absence of a shallow aqueous phase plume directly downgradient of a DNAPL source zone does not rule out the possibility of deep occurrences of DNAPL beyond the shallow monitoring well network. A further finding of this study is that the occurrence of a highly sorbing compound in groundwater at virtually any concentration may indicate the immediate upgradient presence of residual or pooled DNAPL.     

Biogeochemistry and isotope geochemistry of a landfill leachate plume

The biogeochemical processes were identified which improved the leachate composition in the flow direction of a landfill leachate plume (Banisveld, The Netherlands). Groundwater observation wells were placed at specific locations after delineating the leachate plume using geophysical tests to map subsurface conductivity. Redox processes were determined using the distribution of solid and soluble redox species, hydrogen concentrations, concentration of dissolved gases (N(2), Ar, and CH(4)), and stable isotopes (delta15N-NO(3), delta34S-SO(4), delta13C-CH(4), delta2H-CH(4), and delta13C of dissolved organic and inorganic carbon (DOC and DIC, respectively)). The combined application of these techniques improved the redox interpretation considerably. Dissolved organic carbon (DOC) decreased downstream in association with increasing delta13C-DOC values confirming the occurrence of degradation. Degradation of DOC was coupled to iron reduction inside the plume, while denitrification could be an important redox process at the top fringe of the plume. Stable carbon and hydrogen isotope signatures of methane indicated that methane was formed inside the landfill and not in the plume. Total gas pressure exceeded hydrostatic pressure in the plume, and methane seems subject to degassing. Quantitative proof for DOC degradation under iron-reducing conditions could only be obtained if the geochemical processes cation exchange and precipitation of carbonate minerals (siderite and calcite) were considered and incorporated in an inverse geochemical model of the plume. Simulation of delta13C-DIC confirmed that precipitation of carbonate minerals happened.     

In situ sequential treatment of a mixed contaminant plume

Groundwater plumes often contain a mixture of contaminants that cannot easily be remediated in situ using a single technology. The purpose of this research was to evaluate an in situ treatment sequence for the control of a mixed organic plume (chlorinated ethenes and petroleum hydrocarbons) within a Funnel-and-Gate. A shallow plume located in the unconfined aquifer at Alameda Point, CA, was found to contain up to 218,000 μg/l of cis-1,2 dichloroethene (cDCE), 16,000 μg/l of vinyl chloride (VC) and <1000 μg/l of 1,1 dichloroethene (1,1 DCE), trans-1,2 dichloroethene (trans-1,2 DCE) and trichloroethene (TCE). Total benzene, toluene, ethylbenzene and xylenes (BTEX) concentrations were <10,000 μg/l. Contaminated groundwater was funneled into a gate, 3.0 m wide, 4.5 m long and 6.0 m deep (keyed into the underlying aquitard) where treatment occurred. The initial gate segment consisted of granular iron, for the reductive dechlorination of the higher chlorinated ethenes. The second segment, the biosparge zone, promoted aerobic biodegradation of petroleum hydrocarbons and any remaining lesser-chlorinated compounds, stimulated by dissolved oxygen (DO) and carbon dioxide (CO₂) additions via an in situ sparge system (CO₂ was used to neutralize the high pH produced from reactions in the iron wall). Groundwater was drawn through the gate by pumping two wells located at the sealed, downgradient, end. Over a 4-month period an estimated 1350 g of cDCE flowed into the treatment gate and the iron wall removed 1230 g, or 91% of the mass. The influent mass of VC was 572 g and the iron wall removed 535 g, corresponding to 94% mass removal. The other chlorinated ethenes had significantly lower influent masses (3 to 108 g) and the iron wall removed the majority of the mass resulting in >96% mass removal for any of the compounds. In spite of these high removal percentages, laboratory column tests indicated that at these levels of chlorinated contaminants, surface saturation of the iron grains likely contributed to lower than expected reaction rates. In the biosparge zone, mass removal of cDCE appeared to occur predominantly by biodegradation (65%) with volatilization (35%) being an important secondary process. The dominant removal process for VC was volatilization (70%) although significant biodegradation was also indicated (30%). Laboratory microcosm results confirmed the potential for aerobic biodegradation of cDCE and VC. When average influent field concentrations for cDCE and VC were 220,000 and 46,000 μg/l, respectively, the sequential treatment unit removed 99.6% of the total mass and when the influent concentrations decreased to 26,000 and 19,000 μg/l for cDCE and VC, respectively, >99.9% removal within the treatment gate was attained. BTEX compounds were found to be significantly retarded in the iron treatment zone. Although they did eventually break through the granular iron, and into the gravel transition zone, none of these compounds was detected in the biosparge zone. No noticeable interferences between the anaerobic (reductive) and aerobic parts of the system occurred during testing. The results of this experiment show that in situ treatment sequences are viable, although further work is needed to optimize performance.     

Inverse modeling to estimate LNAPL plume release timing

Methodologies are presented for dating releases of light nonaqueous phase liquids (LNAPLs) using an inverse modeling approach with simple analytical models. Models for LNAPL plume migration are presented to predict LNAPL plume velocity in the unsaturated and saturated zones as a function of basic soil and fluid properties. A relative mobility factor is introduced for LNAPL movement at the water table that depends primarily on the van Genuchten n parameter (related to the breadth of the soil pore size distribution) and the magnitude of water table fluctuations. Estimated LNAPL plume velocities compare reasonably with more rigorous numerical models, which may be used in cases where data availability warrant the greater effort entailed.Two methods of estimating release timing and its uncertainty are investigated. A direct estimation method is described that determines travel time for a single observed travel distance based on estimated soil and fluid properties. Release date uncertainty may be determined using the first order (FO) or Monte Carlo (MC) methods. The second method for estimating release date involves nonlinear parameter estimation utilizing distance vs. time measurements and other data.A case study is presented for a field site where independent estimates of release timing were obtained from a numerical modeling analysis. Release timing estimates based on direct inversion of the analytical timing model agree well with the numerical analysis. Results for a second field site indicate that release date confidence limits estimated by the FO method, assuming log-normally distributed travel times, are close to values determined by the MC method, which makes no assumption regarding the form of the travel time probability distribution.Results for a hypothetical problem indicate that LNAPL velocity and travel time may be accurately estimated if sufficient data on travel distance vs. time are available. Incorporating prior information on relevant soil and fluid properties into the objective function reduces the uncertainty in release date if prior estimates are accurate. However, biased prior estimates may lead to over- or underestimation of release date uncertainty. Simultaneous estimation of soil and fluid properties and release date is possible if prior information is available to condition the parameter estimates.     

Cooling tower plume rise analyses by airborne lidar

As part of the Atmospheric Studies in Complex Terrain (ASCOT) program, five cooling-tower plume experiments were conducted at The Geysers geothermal area in California during August 1981. The experiments were designed to investigate plume behavior during conditions of valley nocturnal drainage flow and daytime vertical mixing on a mountain ridge. The Airborne Lidar Plume and Haze Analyzer (ALPHA-1) system was used during the experiments to observe plume geometry and aerosol structure of the boundary layer. Subvisible plume rise derived from the backscatter signatures is related to visible plume-rise results recently published. This analysis indicates that the lidar-observed plumes are about 10 times higher than visually observed plumes. The lidar-observed plume rise seemed to correspond with heights of elevated aerosol layers that typically indicate presence of temperature inversions. Pictorial displays are presented which provide information on atmospheric behavior over mountainous terrain.