Spatiotemporal dynamic temperature variation dominated by ion behaviors during groundwater remediation using direct current.

Direct current (DC) electric field has shown promising performance in contaminated site remediation, in which the Joule heating effect plays an important role but has been previously underappreciated. This study focuses on the spatiotemporal characteristics and mechanism of temperature change in heterogeneous porous media with applied DC. The heating process can be divided into four phases: preferential heating of the low permeability zone (LPZ), rapid heating in the middle region, temperature drop and hot zone shift, and reheating. The dynamic ion behaviors with complex interplays among reactions, electrokinetic-driven migration, and mixed convection induced an uneven redistribution of ions and dominated the heating rate and temperature distribution. The concentration of major ions near the pH jump decreased to 1% of the initial value, even though ions were continuously pumped into the heating zone. This ion depletion caused a drop in current, heating rate, and temperature. Here ions cannot be delivered rapidly into the ion-depleted zone by electromigration due to the potential flattening in the surrounding region. The presence of LPZ intensified the nonuniformity of ion redistribution, where a regional focusing of water-soluble ions was observed, and weakened the temperature rebound compared with that using homogeneous sand. These results provide a new perspective on the regulation of DC heating in site remediation.

Simulating field-scale thermal conductive heating with the potential for the migration and condensation of vapors.

Thermal conductive heating (TCH) is an in-situ thermal treatment (ISTT) technology for treating non-aqueous phase liquid (NAPL) source zones. Numerical models can be useful tools for improving remedial performance, but traditional multiphase flow models are rarely used to simulate mass recovery during ISTT applications at the field scale due to their computational expense. This study developed a 3D model based on macroscopic invasion percolation to simulate the vaporization of NAPL, and the subsequent vapor migration and potential condensation at the field scale. The model was used to simulate the mass recovery of trichloroethene (TCE) from a NAPL source zone under seven scenarios of different heater placements, including three scenarios with an undersized target treatment zone (TTZ). Simulation results showed that TCH was effective in removing NAPL within the TTZ, but the treatment zone did not extend far from the perimeter heaters. In addition, during heating, NAPL condensation outside the TTZ due to the escaping vapor was observed in all scenarios. Overall, the resulting mass recovery was lower in the three scenarios with an undersized TTZ (91-95%) than in the other four scenarios (≈ 99%). Moreover, the locations of unrecovered/condensed NAPL could be inferred by monitoring mass recovery tailing at individual extraction wells.

Modelling gas-phase recovery of volatile organic compounds during in situ thermal treatment.

In situ thermal treatment (ISTT) technologies can be used to remove mass from non-aqueous phase liquid (NAPL) source zones. Ensuring the vaporization of NAPL and the capture of vapors are crucial, and numerical models are useful for understanding the processes that affect performance to help improve design and operation. In this paper, a two-dimensional model that combines a continuum approach based on finite difference for heat transfer with a macroscopic invasion percolation (macro-IP) approach for gas migration was developed to simulate thermal conductive heating (TCH) applications at the field-scale. This approach simulates heat transport and gas migration, but is different than a traditional continuum multiphase approach. Mass recovery for 60 randomly generated realizations under three degrees of heterogeneity of the permeability field were simulated. The mass recovery curves had an overall similar shape for the various permeability fields. However, a wider range of completion times was observed for domains with a higher permeability variance. Results also showed that NAPL pools that were highly saturated, deep, and away from the heaters needed more heating time to be depleted, and that total NAPL mass was not a good indicator of completion time. The completion time was positively correlated with the maximum value of the mixed spatial moment of NAPL saturation about the heaters in the lateral and vertical direction, and the NAPL pool with the highest moment could increase the heating time by as much as 35%. This effect was most notable in simulations with a high permeability variance and suggests the potential to reduce heating time by locating the largest NAPL pools and placing TCH heaters accordingly.

A numerical model for estimating the removal of volatile organic compounds in laboratory-scale treatability tests for thermal treatment of NAPL-impacted soils.

Treatability tests can be carried out to assess the potential effectiveness of thermal treatment technologies under different site conditions and are important for specific technology selection and design. In order to reduce the costs for laboratory tests and expand the insights from previous treatability studies, a one-dimensional (1D) radial finite difference model was developed to simulate the removal of volatile organic compounds (VOCs) in laboratory thermal treatability tests. The processes considered in the model include heat conduction, co-boiling of single-component or multi-component NAPLs with water, and water boiling. An explicit approach is used to simulate the evolution of NAPL composition for multi-component NAPLs during heating. The developed model adopts only two fitting parameters and was calibrated and validated using previous laboratory experiments. In this paper, the developed model was first calibrated to three laboratory experiments using temperature measurements, which resulted in matches to the NAPL and gas saturations. After calibration, the model was able to predict the temperature, NAPL and gas saturations for the remaining seven experiments, including those with single and multi-component NAPLs, using the average value of each fitting parameter.

Arsenic concentration and speciation in five freshwater fish species from Back Bay near Yellowknife, NT, CANADA.

The concentration of total arsenic and five different arsenic species [As(III), As(V), monomethylarsonic acid (MMA), dimethylarsenic acid (DMA), and arsenobetaine (AsB)], were measured in the muscle, liver and gastrointestinal tract (GIT) of five different fish species [lake whitefish (Coregonus clupeaformis), walleye (Stizostedion vitreum), northern pike (Esox lucius), white sucker (Catostomus commersoni) and longnose sucker (Catostomus catostomus)] from Back Bay, Great Slave Lake, near the city of Yellowknife, NT, Canada. The total concentration (dry weight) of arsenic in muscle ranged from 0.57 to 1.15 mg/kg, in the liver from 0.42 to 2.52 mg/kg and in the GIT from 1.48 to 8.92 mg/kg. Among fish species, C. commersoni had significantly higher total arsenic concentrations in the GIT than S. vitreum, E. lucius and C. clupeaformis, and higher total arsenic concentrations in the liver than C. clupeaformis. The mean concentration of As(III) and As(V) in the muscle of all fish ranged from < or =0.01 to 0.05 mg/kg and < or =0.01 to 0.02 mg/kg, respectively, and together comprised < or =7.5% of the total arsenic measured in muscle. The concentrations of MMA were below detection in the muscle of all five fish species. However, AsB and DMA were measured in all fish species and nearly all fish tissues. The concentrations of AsB ranged from 0.01 to 0.13 mg/kg and the concentrations of DMA ranged from <0.02 to 0.45 mg/kg. The majority (>50%) of organic arsenic in almost all of the tissues from fish caught in Back Bay was not directly identified. Evidence from the literature suggests that most of these other organic arsenic species were likely trimethylated arsenic compounds, however, further analytical work would need to be performed to verify this hypothesis.

Effects of arsenic speciation and low dissolved oxygen condition on the toxicity of arsenic to a lotic mayfly.

The influence of site-specific conditions on contaminant bioavailability and toxicity to benthic invertebrates is a key consideration in the environmental risk assessment process. This is particularly relevant for contaminants with complex speciation chemistries, such as arsenic. The present study addressed uncertainties regarding arsenic toxicity to a mayfly (Baetis tricaudatus) under low dissolved oxygen (DO) conditions characteristic of many contaminated sites. Arsenic toxicity (arsenite, As(III); arsenate, As(V)) to mayfly nymphs was assessed under two DO scenarios (68 and 84% saturation). Arsenic speciation ratios were determined during testing to confirm the nature of arsenic exposure. The present study found that As(III) was more lethal and bioaccumulated to a greater degree in B. tricaudatus compared to As(V), but the sublethal toxicities of the two arsenic species were similar. Nymph growth and development were significantly inhibited after 12 d of exposure to both 1 mg/L of As(III) and As(V). Exposure to arsenic under low DO conditions (6.5 mg/L, 68% saturation) did not significantly affect As(III) or As(V) toxicity and bioaccumulation over 12 d. The DO level of 6.5 mg/L, however, appeared to be marginally lethal to B. tricaudatus. Results indicate that the Canadian arsenic criterion for the protection of aquatic life (5 microg/L) is protective of B. tricaudatus and is low enough to accommodate differences in arsenic toxicity because of the interconversions between As(III) and As(V). These findings provide insight regarding the toxicity and speciation of arsenic under DO conditions considered to be low for this lotic mayfly species and representative of existing conditions at mine sites in northern Canada.

Quantifying uranium complexation by groundwater dissolved organic carbon using asymmetrical flow field-flow fractionation.

The long-term mobility of actinides in groundwaters is important for siting nuclear waste facilities and managing waste-rock piles at uranium mines. Dissolved organic carbon (DOC) may influence the mobility of uranium, but few field-based studies have been undertaken to examine this in typical groundwaters. In addition, few techniques are available to isolate DOC and directly quantify the metals complexed to it. Determination of U-organic matter association constants from analysis of field-collected samples compliments laboratory measurements, and these constants are needed for accurate transport calculations. The partitioning of U to DOC in a clay-rich aquitard was investigated in 10 groundwater samples collected between 2 and 30 m depths at one test site. A positive correlation was observed between the DOC (4-132 mg/L) and U concentrations (20-603 microg/L). The association of U and DOC was examined directly using on-line coupling of Asymmetrical Flow Field-Flow Fractionation (AsFlFFF) with UV absorbance (UVA) and inductively coupled plasma-mass spectrometer (ICP-MS) detectors. This method has the advantages of utilizing very small sample volumes (20-50 microL) as well as giving molecular weight information on U-organic matter complexes. AsFlFFF-UVA results showed that 47-98% of the DOC (4-136 mg C/L) was recovered in the AsFlFFF analysis, of which 25-64% occurred in the resolvable peak. This peak corresponded to a weight-average molecular weight of about 900-1400 Daltons (Da). In all cases, AsFlFFF-ICP-MS suggested that