John E. McCray

Bark beetle infestation impacts on nutrient cycling, water quality and interdependent hydrological effects

Bark beetle populations have drastically increased in magnitude over the last several decades leading to the largest-scale tree mortality ever recorded from an insect infestation on multiple wooded continents. When the trees die, the loss of canopy and changes in water and nutrient uptake lead to observable changes in hydrology and biogeochemical cycling. This review aims to synthesize the current research on the effects of the bark beetle epidemic on nutrient cycling and water quality while integrating recent and relevant hydrological findings, along with suggesting necessary future research avenues. Studies generally agree that snow depth will increase in infested forests, though the magnitude is uncertain. Changes in evapotranspiration are more variable as decreased transpiration from tree death may be offset by increased understory evapotranspiration and ground evaporation. As a result of such competing hydrologic processes that can affect watershed biogeochemistry along with the inherent variability of natural watershed characteristics, water quality changes related to beetle infestation are difficult to predict and may be regionally distinct. However, tree-scale changes to soil–water chemistry (N, P, DOC and base cation concentrations and composition) are being observed in association with beetle outbreaks which ultimately could lead to larger-scale responses. The different temporal and spatial patterns of bark beetle infestations due to different beetle and tree species lead to inconsistent infestation impacts. Climatic variations and large-scale watershed responses provide a further challenge for predictions due to spatial heterogeneities within a single watershed; conflicting reports from different regions suggest that hydrologic and water quality impacts of the beetle on watersheds cannot be generalized. Research regarding the subsurface water and chemical flow-paths and residence times after a bark beetle epidemic is lacking and needs to be rigorously addressed to best predict watershed or regional-scale changes to soil–water, groundwater, and stream water chemistry.

Biosurfactant-enhanced solubilization of NAPL mixtures

Remediation of nonaqueous phase liquids (NAPLs) by conventional pump-and-treat methods (i.e., water flushing) is generally considered to be ineffective due to low water solubilities of NAPLs and to mass-transfer constraints. Chemical flushing techniques, such as surfactant flushing, can greatly improve NAPL remediation primarily by increasing the apparent solubility of NAPL contaminants. NAPLs at hazardous waste sites are often complex mixtures. However, the equilibrium and nonequilibrium mass-transfer characteristics between NAPL mixtures and aqueous surfactant solutions are not well understood. This research investigates the equilibrium solubilization behavior of two- and three-component NAPL mixtures (containing akylbenzenes) in biosurfactant solutions. NAPL solubilization is found to be ideal in water (i.e., obeys Raoults Law), while solubilization in biosurfactant solutions was observed to be nonideal. Specifically, the relatively hydrophobic compounds in the mixture experienced solubility enhancements that were greater than those predicted by ideal enhanced solubilization theory, while the solubility enhancements for the relatively hydrophilic compounds were less than predicted. The degree of nonideality is shown to be a nonlinear function of the NAPL-phase mole fraction. Empirical relationships based on the NAPL-phase mole fraction and/or micelle-aqueous partition coefficients measured in single-component NAPL systems are developed to estimate values for the multicomponent partition coefficients. Empirical relationships that incorporate both the NAPL-phase mole fraction and single-component partition coefficients yield much improved estimates for the multicomponent partition coefficient.

Metal release from sandstones under experimentally and numerically simulated CO2 leakage conditions.

Leakage of CO2 from a deep storage formation into an overlying potable aquifer may adversely impact water quality and human health. Understanding CO2-water-rock interactions is therefore an important step toward the safe implementation of geologic carbon sequestration. This study targeted the geochemical response of siliclastic rock, specifically three sandstones of the Mesaverde Group in northwestern Colorado. To test the hypothesis that carbonate minerals, even when present in very low levels, would be the primary source of metals released into a CO2-impacted aquifer, two batch experiments were conducted. Samples were reacted for 27 days with water and CO2 at partial pressures of 0.01 and 1 bar, representing natural background levels and levels expected in an aquifer impacted by a small leakage, respectively. Concentrations of major (e.g., Ca, Mg) and trace (e.g., As, Ba, Cd, Fe, Mn, Pb, Sr, U) elements increased rapidly after CO2 was introduced into the system, but did not exceed primary Maximum Contaminant Levels set by the U.S. Environmental Protection Agency. Results of sequential extraction suggest that carbonate minerals, although volumetrically insignificant in the sandstone samples, are the dominant source of mobile metals. This interpretation is supported by a simple geochemical model, which could simulate observed changes in fluid composition through CO2-induced calcite and dolomite dissolution.

The effect of system variables on in situ sweep-efficiency improvements via viscosity modification

Laboratory experiments and numerical simulations were performed to critically evaluate the utility of viscosity modification as a technique to improve injected fluid sweep efficiencies within texturally heterogeneous geomedia. The objective of this technique is to improve the subsurface distribution of fluids by mitigating the potential for preferential flow and bypassing of lower permeability media that can limit the effectiveness of in situ remediation applications. The results of two-dimensional sand tank experiments and numerical simulations demonstrate that viscosity modification, via polymer amendment, can improve sweep efficiencies within layered heterogeneous structures by up to 90%, relative to the no-polymer case. The amount of sweep efficiency improvement depended on a number of system variables, including: the degree of layering, the relative positioning of layers within the system, the permeability contrast between layers, fluid viscosity, and the rheological character of the fluid utilized. Although significant sweep-efficiency improvement was observed, achieving 100% sweep in one pore volume was only possible when the permeability contrast was less than a factor of four, regardless of the viscosity and the rheological character of the fluid.

Bioaugmentation for treatment of dense non-aqueous phase liquid in fractured sandstone blocks.

Laboratory experiments were performed in discretely fractured sandstone blocks to evaluate the use of bioaugmentation to treat residual dense non-aqueous phase liquid (DNAPL) tetrachloroethene (PCE). Significant dechlorination of PCE and growth of Dehalococcoides spp. (DHC) occurred within the fractures. DNAPL dissolution was enhanced during bioaugmentation by up to a factor of approximately 3.5, with dissolved PCE concentrations at or near aqueous solubility. The extent of dechlorination and DNAPL dissolution enhancement were dependent upon the fracture characteristics, residence time in the fractures, and dissolved concentration of PCE. No relationship was observed between planktonic DHC concentrations exiting the fracture and the observed extents of PCE dechlorination and DNAPL dissolution. Measured planktonic DHC concentrations exiting the fracture increased with increasing flow rate and bioaugmentation dosage, suggesting that these parameters may be important for distribution of DHC to treat dissolved chlorinated ethenes migrating downgradient of the DNAPL source. Bioaugmentation dosage, for the DHC dosages and conditions studied, did not have a measurable impact on DNAPL dissolution or dechlorination within the fractures themselves. Overall, these results indicate that bioaugmentation may be a viable remedial option for treating DNAPL sources in bedrock.

Solubility enhancement of seven metal contaminants using carboxymethyl-β-cyclodextrin (CMCD)

Carboxymethyl-beta-cyclodextrin (CMCD) has been suggested as a complexing agent for remediation of sites co-contaminated with metals and organic pollutants. As part of an attempt to construct a geochemical complexation model for metal-CMCD interactions, conditional formation constants for the complexes between CMCD and 7 metal ions (Ba, Ca, Cd, Ni, Pb, Sr, and Zn) are estimated from experimental data. Stable metal concentrations were reached after approximately 1 day and estimated logarithmic conditional formation constants range from -3.2 to -5.1 with confidence intervals within +/-0.08 log units. Experiments performed at 10 degrees C and 25 degrees C show that temperature affects the solubility of the metal salts but the strength of CMCD-metal complexes are not affected by this temperature variation. The conditional stability constants and complexation model presented in this work can be used to screen CMCD as a potential remediation agent for clean-up of contaminated soil and groundwater.

Kinetic metal release from competing processes in aquifers.

Understanding groundwater time scales wherein kinetic metal-desorption and mineral-dissolution are important mechanisms is essential for realistic modeling of metal release. In this study, release rate constants were compiled and the Damköhler number was applied to calculate residence times where kinetic formulations are relevant. Desorption rate constants were compiled for arsenic, barium, cadmium, copper, lead, mercury, nickel, and zinc, and span 6 orders of magnitude, while mineral-dissolution rate constants compiled for calcite, kaolinite, smectite, anorthite, albite, K-feldspar, muscovite, quartz, goethite, and galena ranged over 13 orders of magnitude. This Damköhler analysis demonstrated that metal-desorption kinetics are potentially influential at residence times up to about two years, depending on the metal and groundwater conditions. Kinetic mineral-dissolution should be considered for nearly all residence times relevant to groundwater modeling, provided the rate, solubility, and availability of the mineral generates a non-negligible concentration. Geochemical models of competitive desorption and dissolution for an illustrative metal demonstrate total metal concentrations may be sensitive to dissolution rate variations despite the predominance of release from desorption. Ultimately, this analysis provides constraints on relevant processes for incorporation into transport models.

Effects of chemical oxidants on perfluoroalkyl acid transport in one-dimensional porous media columns.

In situ chemical oxidation (ISCO) is a remediation approach that is often used to remediate soil and groundwater contaminated with fuels and chlorinated solvents. At many aqueous film-forming foam-impacted sites, perfluoroalkyl acids (PFAAs) can also be present at concentrations warranting concern. Laboratory experiments were completed using flow-through one-dimensional columns to improve our understanding of how ISCO (i.e., activated persulfate, permanganate, or catalyzed hydrogen peroxide) could affect the fate and transport of PFAAs in saturated porous media. While the resultant data suggest that standard ISCO is not a viable remediation strategy for PFAA decomposition, substantial changes in PFAA transport were observed upon and following the application of ISCO. In general, activated persulfate decreased PFAA transport, while permanganate and catalyzed hydrogen peroxide increased PFAA transport. PFAA sorption increased in the presence of increased aqueous polyvalent cation concentrations or decreased pH. The changes in contaminant mobility were greater than what would be predicted on the basis of aqueous chemistry considerations alone, suggesting that the application of ISCO results in changes to the porous media matrix (e.g., soil organic matter quality) that also influence transport. The application of ISCO is likely to result in changes in PFAA transport, where the direction (increased or decreased transport) and magnitude are dependent on PFAA characteristics, oxidant characteristics, and site-specific factors.

Engineered Streambeds for Induced Hyporheic Flow: Enhanced Removal of Nutrients, Pathogens, and Metals from Urban Streams

AbstractThe hyporheic zone (HZ) has the potential to mitigate nonpoint source pollution that threatens urban streams, but limited flows and inefficient interchange typically constrain its water quality improvements. This paper presents coordinated streambed hydraulic conductivity (K) modifications, termed Biohydrochemical Enhancement Structures for Streamwater Treatment (BEST), to efficiently drive interchange, attenuating nutrients, metals, and pathogens. Numerical models consisted of sandy HZ with high and low K streambed media components, and variable slope and in situ sediment K. BEST contaminant mitigation was estimated using hyporheic flow and residence time (RT) simulations and contaminant removal rate constants from literature for specific porous media. Hyporheic flows improved, whereas mean RT decreased approximately linearly with increasing slope and in situ K. Despite vast differences in flow and RT, contaminant removal potential was relatively consistent among BEST in which reaction time scale...

Metal fate and partitioning in soils under bark beetle-killed trees

Recent mountain pine beetle infestation in the Rocky Mountains of North America has killed an unprecedented acreage of pine forest, creating an opportunity to observe an active re-equilibration in response to widespread land cover perturbation. This work investigates metal mobility in beetle-impacted forests using parallel rainwater and acid leaches to estimate solid-liquid partitioning coefficients and a complete sequential extraction procedure to determine how metals are fractionated in soils under trees experiencing different phases of mortality. Geochemical model simulations analyzed in consideration with experimental data provide additional insight into the mechanisms controlling metal complexation. Metal and base-cation mobility consistently increased in soils under beetle-attacked trees relative to soil under healthy trees. Mobility increases were more pronounced on south facing slopes and more strongly correlated to pH under attacked trees than under healthy trees. Similarly, soil moisture was significantly higher under dead trees, related to the loss of transpiration and interception. Zinc and cadmium content increased in soils under dead trees relative to living trees. Cadmium increases occurred predominantly in the exchangeable fraction, indicating increased mobilization potential. Relative increases of zinc were greatest in the organic fraction, the only fraction where increases in copper were observed. Model results reveal that increased organic complexation, not changes in pH or base cation concentrations, can explain the observed differences in metal partitioning for zinc, nickel, cadmium, and copper. Predicted concentrations would be unlikely to impair human health or plant growth at these sites; however, higher exchangeable metals under beetle-killed trees relative to healthy trees suggest a possible decline in riverine ecosystem health and water quality in areas already approaching criteria limits and drinking water standards. Impairment of water quality in important headwater streams from the increased potential for metal mobilization and storage will continue to change as beetle-killed trees decompose and forests begin to recover.