Effect of soil moisture dynamics on dense nonaqueous phase liquid (DNAPL) spill zone architecture in heterogeneous porous media

The amount, location, and form of NAPL in contaminated vadose zones are controlled by the spatial distribution of water saturation and soil permeability, the NAPL spill scenario, water infiltration events, and vapor transport. To evaluate the effects of these processes, we used the three-phase flow simulator STOMP, which includes a new permeability-liquid saturation-capillary pressure (k-S-P) constitutive model. This new constitutive model considers three NAPL forms: free, residual, and trapped. A 2-D vertical cross-section with five stratigraphic layers was assumed, and simulations were performed for seven cases. The conceptual model of the soil heterogeneity was based upon the stratigraphy at the Hanford carbon tetrachloride (CT) spill site. Some cases considered co-disposal of NAPL with large volumes of wastewater, as also occurred at the Hanford CT site. In these cases, the form and location of NAPL were most strongly influenced by high water discharge rates and NAPL evaporation to the atmosphere. In order to investigate the impact of heterogeneity, the hydraulic conductivity within the lower permeability layer was modeled as a realization of a random field having three different classes. For six extreme cases of 100 realizations, the CT mass that reached the water table varied by a factor of two, and was primarily controlled by the degree of lateral connectivity of the low conductivity class within the lowest permeability layer. The grid size at the top boundary had a dramatic impact on NAPL diffusive flux just after the spill event when the NAPL was present near the ground surface. NAPL evaporation with a fine grid spacing at the top boundary decreased CT mass that reached the water table by 74%, compared to the case with a coarse grid spacing, while barometric pumping had a marginal effect for the case of a continuous NAPL spill scenario considered in this work. For low water infiltration rate scenarios, the distribution of water content prior to a NAPL spill event decreased CT mass that reached the water table by 98% and had a significant impact on the formation of trapped NAPL. For all cases simulated, use of the new constitutive model that allows the formation of residual NAPL increased the amount of NAPL retained in the vadose zone. Density-driven advective gas flow from the ground surface controlled vapor migration in strongly anisotropic layers, causing NAPL mass flux to the lower layer to be reduced. These simulations indicate that consideration of the formation of residual and trapped NAPLs and dynamic boundary conditions (e.g., areas, rates, and periods of different NAPL and water discharge and fluctuations of atmospheric pressure) in the context of full three-phase flow are needed, especially for NAPL spill events at the ground surface. In addition, NAPL evaporation, density-driven gas advection, and NAPL vertical movement enhanced by water flow must be considered in order to predict NAPL distribution and migration in the vadose zone.     

Modeling the effects of concentration history on the slow desorption of trichloroethene from a soil at 100% relative humidity

To determine the effects of concentration history on slow sorption processes, desorption kinetic profiles for trichloroethene (TCE) were measured for a soil at 100% relative humidity subject to different exposure concentrations and exposure times. Exposure concentrations ranged from 1% to 80% of the saturation vapor pressure (Ps) for TCE, and exposure times (i.e., time allowed for sorption before desorption begins) ranged from 1 to 96 days. A spherical diffusion model based on a gamma distribution of sorption rates and a gamma distribution of desorption rates was developed and applied to the data. At 80% P/Ps, the entire gamma distributions of sorption and desorption rates were available for TCE. In accordance with a micropore filling mechanism, the fraction of these distributions available for TCE sorption decreased with decreasing P/Ps. Experimental results are consistent with a micropore-filling mechanism, where the amount of slow desorbing mass decreased with decreasing exposure time, and the fraction of slow desorbing sites filled decreased with decreasing exposure concentration. Simulation results suggest that diffusion limits the rates that micropores fill, and that rates of sorption and desorption for soil contaminated at smaller values of P/Ps are, on average, less than those at larger values of P/Ps (i.e., slow desorption rates are a function of exposure concentration). Simulation results also suggest that the model adequately describes the effects of exposure concentration and exposure time on the rates of sorption and desorption, but not on the capacity of the slow sites for TCE. This work is important because contaminant concentrations in the subsurface vary in space and time, and the proposed model represents a new and mechanistically based approach to capture the effects of this heterogeneity on slow desorption.     

Quantifying dispersal and establishment limitation in a population of an epiphytic lichen

Dispersal is a process critical for the dynamics and persistence of metapopulations, but it is difficult to quantify. It has been suggested that the old-forest lichen Lobaria pulmonaria is limited by insufficient dispersal ability. We analyzed 240 DNA extracts derived from snow samples by a L. pulmonaria-specific real-time PCR (polymerase chain reaction) assay of the ITS (internal transcribed spacer) region allowing for the discrimination among propagules originating from a single, isolated source tree or propagules originating from other locations. Samples that were detected as positives by real-time PCR were additionally genotyped for five L. pulmonaria microsatellite loci. Both molecular approaches demonstrated substantial dispersal from other than local sources. In a landscape approach, we additionally analyzed 240 snow samples with real-time PCR of ITS and detected propagules not only in forests where L. pulmonaria was present, but also in large unforested pasture areas and in forest patches where L. pulmonaria was not found. Monitoring of soredia of L. pulmonaria transplanted to maple bark after two vegetation periods showed high variance in growth among forest stands, but no significant differences among different transplantation treatments. Hence, it is probably not dispersal limitation that hinders colonization in the old-forest lichen L. pulmonaria, but ecological constraints at the stand level that can result in establishment limitation. Our study exemplifies that care has to be taken to adequately separate the effects of dispersal limitation from a limitation of establishment.     

Impact of climate change and seasonal trends on the fate of Arctic oil spills

We investigated the effects of a warmer climate, and seasonal trends, on the fate of oil spilled in the Arctic. Three well blowout scenarios, two shipping accidents and a pipeline rupture were considered. We used ensembles of numerical simulations, using the OSCAR oil spill model, with environmental data for the periods 2009–2012 and 2050–2053 (representing a warmer future) as inputs to the model. Future atmospheric forcing was based on the IPCC’s A1B scenario, with the ocean data generated by the hydrodynamic model SINMOD. We found differences in “typical” outcome of a spill in a warmer future compared to the present, mainly due to a longer season of open water. We have demonstrated that ice cover is extremely important for predicting the fate of an Arctic oil spill, and find that oil spills in a warming climate will in some cases result in greater areal coverage and shoreline exposure.     

Investigation of surfactant-enhanced mass removal and flux reduction in 3D correlated permeability fields using magnetic resonance imaging

Magnetic resonance imaging (MRI) was used to visualize the NAPL source zone architecture before and after surfactant-enhanced NAPL dissolution in three-dimensional (3D) heterogeneously packed flowcells characterized by different longitudinal correlation lengths: 2.1 cm (aquifer 1) and 1.1 cm (aquifer 2). Surfactant flowpaths were determined by imaging the breakthrough of a paramagnetic tracer (MnCl(2)) analyzed by the method of moments. In both experimental aquifers, preferential flow occurred in high permeability materials with low NAPL saturations, and NAPL was preferentially removed from the top of the aquifers with low saturation. Alternate flushing with water and two surfactant pulses (5-6 pore volumes each) resulted in approximately 63% of NAPL mass removal from both aquifers. However, overall reduction in mass flux (Mass Flux 1) exiting the flowcell was lower in aquifer 2 (68%) than in aquifer 1 (81%), and local effluent concentrations were found to increase by as high as 120 times at local sampling ports from aquifer 2 after surfactant flushing. 3D MRI images of NAPL revealed that NAPL migrated downward and created additional NAPL source zones in previously uncontaminated areas at the bottom of the aquifers. The additional NAPL source zones were created in the direction transverse to flow in aquifer 2, which explains the higher mass flux relative to aquifer 1. Analysis using a total trapping number indicates that mobilization of NAPL trapped in the two coarsest sand fractions is possible when saturation is below 0.5 and 0.4, respectively. Results from this study highlight the potential impacts of porous media heterogeneity and NAPL source zone architecture on advanced in-situ flushing technologies.     

Evaluation of simplified mass transfer models to simulate the impacts of source zone architecture on nonaqueous phase liquid dissolution in heterogeneous porous media

Nonaqueous phase liquid (NAPL) dissolution was studied in three-dimensional (3D) heterogeneous experimental aquifers (25.5 cm x 9 cm x 8.5 cm) with two different longitudinal correlation lengths (2.1 cm and 1.1 cm) and initial spill volumes (22.5 ml and 10.5 ml). Spatial and temporal distributions of NAPL during dissolution were measured using magnetic resonance imaging (MRI). At high NAPL spill volume, average effluent concentrations initially increased during dissolution, as NAPL pools transitioned to NAPL ganglia, and then decreased as the total NAPL-water interfacial area decreased over time. Experimental results were used to test six dissolution models: (i and ii) a one-dimensional (1D) model using either specific NAPL-water interfacial area values estimated from MR images at each time step (i.e., 1D quasi-steady state model), or an empirical mass transfer (Sh') correlation (i.e., 1D transient model), (iii and iv) a multiple analytical source superposition technique (MASST) using either the NAPL distribution determined from MR images at each time step (i.e., MASST steady state model), or the NAPL distribution determined from mass balance calculations (i.e., MASST transient model), (v) an equilibrium streamtube model, and (vi) a 3D grid-scale pool dissolution model (PDM) with a dispersive mass flux term. The 1D quasi-steady state model and 3D PDM captured effluent concentration values most closely, including some concentration fluctuations due to changes in the extent of flow reduction. The 1D transient, MASST steady state and transient, and streamtube models all showed a monotonic decrease in effluent concentration values over time, and the streamtube model was the most computationally efficient. Changes during dissolution of the effective NAPL-water interfacial area estimated from imaging data are similar to changes in effluent concentration values. The 1D steady state model incorporates estimates of the effective NAPL-water interfacial area directly at each time point; the 3D PDM does so indirectly through mass balance and a relative permeability function, which causes reduced water flow through high saturation NAPL regions. Hence, when model accuracy is required, the results indicate that a surrogate of this effective interfacial area is required. Approaches to include this surrogate in the MASST and streamtube models are recommended.     

Environmental assessment of ester-based lubricants after application

GOAL, SCOPE AND BACKGROUND: Lubricants based on renewable resources are an environmentally friendly alternative to petrochemical products due to their better ecotoxicological performance and their excellent biodegradability. To improve the technical performance of lubricants, and to reduce friction and wear, the use of additives is nowadays obligatory. The collaborative research center SFB 442 aims at developing environmentally acceptable lubricants that facilitate the avoidance of these additives by transferring their function to suitable coatings. For a complete assessment of the ecological performance of these newly developed lubricants, the whole life cycle including production, application as well as disposal and fate in the environment is studied. In the following study the focus was on the application and its influence on the environmental behavior of the lubricant. The application of lubricants leads, among other things, to the intake of metals due to abrasion from tools, work pieces or mechanical components. Previous examinations indicated a possible influence of metals on the toxicity of eluates prepared from used lubricants (Erlenkaemper et al. 2005). To clarify if the apparent toxicity of used lubricants solely results from the intake of metals, the extractability of these metals from the oil matrix is determined. By combining chemical analyses with bioassays, the bioavailability of metals that are present in the extract is estimated. To further investigate the relevance of metals on toxicity, toxic units (TU) were calculated and related to the results of the bioassays. Interactions between the metals were investigated with aqueous mixtures of metal chlorides and calculations based on the concept of concentration addition and independent action. METHODS: A lubricant mixture was applied to a tribological test bench that simulates real conditions of use and extremely short time load, respectively. Samples were taken at particular times, water soluble fractions (WSF) of these fluids were prepared and dilution series were investigated in several bioassays. Concentration of metals and total organic carbon (TOC) were determined in the eluates. TUs were calculated according to Sprague (1970) and mixture toxicity was calculated according to the concept of concentration addition (Loewe and Muischnek 1926) and independent action (Bliss 1939). RESULTS: Analyses of the metal content of the lubricant and the eluates clearly revealed the availability of the metals in the aqueous extracts. Especially copper, zinc, nickel and chromium were found and their concentrations increased during the time of use. The water extractable fraction, e.g., of copper, rose from 8.8% to 45.3% of the total content in the lubricant after 33.5 hours of use. Tests performed with the algal growth inhibition assay and the luminescence inhibition assay revealed the uptake or absorption by the organisms and, thus, the bioavailability of the metals. The calculation of TUs partly indicated a possible influence of the metals on ecotoxicity of the eluates. Copper reached concentrations equal to or higher than the EC50 value of copper chloride in the growth inhibition assays with algae and Ps. putida as well in the immobilization assay with daphnids. TUs for copper are also larger than 1 for the algal growth inhibition assay. The EL50 values indicated that the luminescence inhibition assay was the most sensitive test system, with values between 4.7% and 9.6%. While the toxicity towards algae and V. fischeri in the growth inhibition assay decreased until both organisms were no longer influenced by the exposure, the EL50 values for the D. magna immobilization assay and the Ps. putida growth inhibition assay decreased with the progressing use of the lubricant. The tested metal salt mixtures showed that Ps. putida, algae and daphnids are the most sensitive organisms with EC50 values below 1 mg/l. DISCUSSION: Although the intake of metals mainly occurred via abrasion of particles, the results revealed the availability of these metals in water. The availability varied for each of the four metals. For both the algal growth inhibition assay and the luminescence inhibition assay, an uptake or absorption of the metals could be demonstrated. The calculated TUs indicated an effect in some bioassays that was not verified in the test itself. The influence of copper on V. fischeri, for example, was not confirmed while the EL50 values for the daphnid bioassay decreased, meaning that the eluates became more toxic with progressing use of the lubricant. The calculations of mixture toxicity based on the concept of concentration addition demonstrated good correlations with the tested metal mixtures, but also a different sensitivity of the organisms. CONCLUSIONS: The results presented here reveal the availability of those metals in water that were taken in during the use of the lubricant in a tribological test bench and, thus, have the possibility of interacting with the organisms. The availability of the metals in the bioassays was proven by chemical analyses. The calculation of TUs and the corresponding EL50 values, however, indicate different availabilities of the metals. The results of the metal salt mixtures show good correlations with calculations of mixture toxicity based on concentration addition. Moreover, the varying sensitivity of the organisms when exposed to eluates or metal mixtures indicates a different bioavailability of the metals and/or the presence of other compounds that exert toxic action. RECOMMENDATIONS AND PERSPECTIVES: For further investigations, the organic oil matrix and its influence on the toxicity have to be taken into account. The toxicity of the eluates may not only be due to metals; additional effects could arise from changes in the lubricant itself.     

Magnetic resonance imaging of nonaqueous phase liquid during soil vapor extraction in heterogeneous porous media

Soil vapor extraction (SVE) is commonly used to remediate nonaqueous phase liquids (NAPLs) from the vadose zone. This paper aims to determine the effect of grain size heterogeneity on the removal of NAPL in porous media during SVE. Magnetic resonance imaging (MRI) was used to observe and quantify the amount and location of NAPL in flow-through columns filled with silica gel grains. MRI is unique because it is nondestructive, allowing three-dimensional images to be taken of the phases as a function of space and time. Columns were packed with silica gel in three ways: coarse grains (250-550 microm) only, fine grains (32-63 microm) only, and a core of fine grains surrounded by a shell of coarse grains. Columns saturated with water were drained under a constant suction head, contaminated with decane, and then drained to different decane saturations. Each column was then continuously purged with water-saturated nitrogen gas and images were taken intermittently. Results showed that at residual saturation, a sharp volatilization front moved through the columns filled with either coarse-grain or fine-grain silica gel. In the heterogeneous columns, the volatilization front in the core lagged just behind the shell because gas flow was greater through the shell and decane in the core diffused outward to the shell. When decane saturation in the core was above residual saturation, decane volatilization occurred near the inlet, the relative decane saturation throughout the core dropped uniformly, and decane in the core flowed in the liquid phase to the shell to replenish volatilized decane. These results indicate that NAPL trapped in low-permeability zones can flow to replenish areas where NAPL is lost due to SVE. However, when residual NAPL saturation is reached, NAPL flow no longer occurs and diffusion limits removal from low-permeability zones.     

Oil spill response capabilities and technologies for ice-covered Arctic marine waters: A review of recent developments and established practices

Renewed political and commercial interest in the resources of the Arctic, the reduction in the extent and thickness of sea ice, and the recent failings that led to the Deepwater Horizon oil spill, have prompted industry and its regulatory agencies, governments, local communities and NGOs to look at all aspects of Arctic oil spill countermeasures with fresh eyes. This paper provides an overview of present oil spill response capabilities and technologies for ice-covered waters, as well as under potential future conditions driven by a changing climate. Though not an exhaustive review, we provide the key research results for oil spill response from knowledge accumulated over many decades, including significant review papers that have been prepared as well as results from recent laboratory tests, field programmes and modelling work. The three main areas covered by the review are as follows: oil weathering and modelling; oil detection and monitoring; and oil spill response techniques.