The velocity of DNAPL fingering in water-saturated porous media: laboratory experiments and a mobile–immobile–zone model

Dense nonaqueous phase liquids (DNAPLs) are immiscible with water and can give rise to highly fingered fluid distributions when infiltrating through water-saturated porous media. In this paper, a conceptual mobile–immobile–zone (MIZ) model is presented to describe the structure of a DNAPL finger in water-saturated porous media and the velocity of finger propagation. A finger is composed of a finger body and a tip. The finger body has a mobile core and an immobile sheath. All the DNAPL within the tip of a finger is mobile. Lab experiments utilizing image analyses of a DNAPL (PCE) penetrating into water-saturated homogeneous glass beads were carried out in a two-dimensional transparent chamber. The results show that the fingers elongated almost linearly with time. The fingers did not grow laterally after the tip of the finger had passed. The average finger diameters were between 3.9 and 5.4 mm for PCE propagation in water-saturated glass bead porous media with mean particle diameters from 0.32 to 1.36 mm. The estimated mobile core diameters were 51–60% of the average finger diameters.     

Contaminant transport in dual-porosity media with dissolved organic matter and bacteria present as mobile colloids

In riverbank filtration, contaminant transport is affected by colloidal particles such as dissolved organic matter (DOM) and bacterial particles. In addition, the subsurface heterogeneity influences the behavior of contaminant transport in riverbank filtration. A mathematical model is developed to describe the contaminant transport in dual-porosity media in the presence of DOM and bacteria as mobile colloids. In the model development, a porous medium is divided into the mobile and immobile regions to consider the presence of ineffective micropores in physically heterogeneous riverbanks. We assume that the contaminant transport in the mobile region is controlled by the advection and dispersion while the contaminant transport in the immobile region occurs due to the molecular diffusion. The contaminant transfer between the mobile and immobile regions takes place by diffusive mass transfer. The mobile region is conceptualized as a four-phase system: two mobile colloidal phases, an aqueous phase, and a solid matrix. The complete set of governing equations is solved numerically with a fully implicit finite difference method. The model results show that in riverbank filtration, the contaminant can migrate further than expected due to the presence of DOM and bacteria. In addition, the contaminant mobility increases further in the presence of the immobile region in aquifers. A sensitivity analysis shows that in dual-porosity media, earlier breakthrough of the contaminant takes place as the volumetric fraction of the mobile region decreases. It is also demonstrated that as the contaminant mass transfer rate coefficient between the mobile and immobile regions increases, the contaminant concentration gradient between the two regions reverses at earlier pore volumes. The contaminant mass transfer coefficient between the mobile and immobile regions mainly controls the tailing effect of the contaminant breakthrough. The contaminant breakthrough curves are sensitive to changes in contaminant adsorption and desorption rate coefficients on DOM and bacteria. In situations where the contaminant is released in the presence of DOM and bacteria in dual-porosity media, the early breakthrough and tailing occur due to the colloidal facilitation and presence of immobile regions.     

Modelling the migration of contaminants through variably saturated dual-porosity, dual-permeability chalk

In the Hesbaye region in Belgium, tracer tests performed in variably saturated fissured chalk rocks presented very contrasting results in terms of transit times, according to artificially controlled water recharge conditions prevailing during the experiments. Under intense recharge conditions, tracers migrated across the partially or fully saturated fissure network, at high velocity in accordance with the high hydraulic conductivity and low effective porosity (fracture porosity). At the same time, a portion of the tracer was temporarily retarded in the almost immobile water located in the matrix. Under natural infiltration conditions, the fissure network remained inactive. Tracers migrated downward through the matrix, at low velocity in relation with the low hydraulic conductivity and the large porosity of the matrix. Based on these observations, Brouyère et al. (2004a) [Brouyère, S., Dassargues, A., Hallet, V., 2004a. Migration of contaminants through the unsaturated zone overlying the Hesbaye chalky aquifer in Belgium: a field investigation, J. Contam. Hydrol., 72 (1-4), 135-164, doi: 10.1016/j.conhyd.2003.10.009] proposed a conceptual model in order to explain the migration of solutes in variably saturated, dual-porosity, dual-permeability chalk. Here, mathematical and numerical modelling of tracer and contaminant migration in variably saturated fissured chalk is presented, considering the aforementioned conceptual model. A new mathematical formulation is proposed to represent the unsaturated properties of the fissured chalk in a more dynamic and appropriate way. At the same time, the rock water content is partitioned between mobile and immobile water phases, as a function of the water saturation of the chalk rock. The groundwater flow and contaminant transport in the variably saturated chalk is solved using the control volume finite element method. Modelling the field tracer experiments performed in the variably saturated chalk shows the adequacy and usefulness of the new conceptual, mathematical and numerical model.     

Impact of varying soil structure on transport processes in different diagnostic horizons of three soil types

When soil structure varies in different soil types and the horizons of these soil types, it has a significant impact on water flow and contaminant transport in soils. This paper focuses on the effect of soil structure variations on the transport of pesticides in the soil above the water table. Transport of a pesticide (chlorotoluron) initially applied on soil columns taken from various horizons of three different soil types (Haplic Luvisol, Greyic Phaeozem and Haplic Cambisol) was studied using two scenarios of ponding infiltration. The highest infiltration rate and pesticide mobility were observed for the Bt₁ horizon of Haplic Luvisol that exhibited a well-developed prismatic structure. The lowest infiltration rate was measured for the Bw horizon of Haplic Cambisol, which had a poorly developed soil structure and a low fraction of large capillary pores and gravitational pores. Water infiltration rates were reduced during the experiments by a soil structure breakdown, swelling of clay and/or air entrapped in soil samples. The largest soil structure breakdown and infiltration decrease was observed for the Ap horizon of Haplic Luvisol due to the low aggregate stability of the initially well-aggregated soil. Single-porosity and dual-permeability (with matrix and macropore domains) flow models in HYDRUS-1D were used to estimate soil hydraulic parameters via numerical inversion using data from the first infiltration experiment. A fraction of the macropore domain in the dual-permeability model was estimated using the micro-morphological images. Final soil hydraulic parameters determined using the single-porosity and dual-permeability models were subsequently used to optimize solute transport parameters. To improve numerical inversion results, the two-site sorption model was also applied. Although structural changes observed during the experiment affected water flow and solute transport, the dual-permeability model together with the two-site sorption model proved to be able to approximate experimental data.     

Simulating the fate and transport of TCE from groundwater to indoor air

This work provides an exploratory analysis on the relative importance of various factors controlling the fate and transport of volatile organic contaminants (in this case, TCE) from a DNAPL source zone located below the water table and into the indoor air. The analysis is conducted using the multi-phase compositional model CompFlow Bio, with the base scenario problem geometry reminiscent of a field experiment conducted by Rivett [Rivett, M.O., (1995), Soil–gas signatures from volatile chlorinated solvents: Borden field experiments. Groundwater, 33(1), 84–98.] at the Borden aquifer where groundwater was observed to transport a contaminant plume a substantial distance without vertical mass transport of the contaminant across the capillary fringe and into the vadose zone. Results for the base scenario model indicate that the structure of the permeability field was largely responsible for deflecting the groundwater plume upward towards the capillary fringe, permitting aqueous phase diffusion to transport the TCE into the vadose zone. Alternative permeability realizations, generated as part of a Monte Carlo simulation process, at times deflected the groundwater plume downwards causing the extended thickness of the saturated zone to insulate the vadose zone from exposure to the TCE by upward diffusive transport. Comparison of attenuation coefficients calculated using the CompFlow Bio and Johnson and Ettinger [Johnson, P.C. and Ettinger, R.A., (1991), Heuristic model for predicting the intrusion rate of contaminant vapors into buildings. Environmental Science and Technology, 25, 1445–1452.] heuristic model exhibited fortuitous agreement for the base scenario problem geometry, with this agreement diverging for the alternative permeability realizations as well as when parameters such as the foundation slab fracture aperture, the indoor air pressure drop, the capillary fringe thickness, and the infiltration rate were varied over typical ranges.     

Implementation of a windblown dust parameterization into MODELS-3/CMAQ: Application to episodic PM events in the US/Mexico border

Windblown dust is known to impede visibility, deteriorate air quality and modify the radiation budget. Arid and semiarid areas with unpaved and unvegetated land cover are particularly prone to windblown dust, which is often attributed to high particulate matter (PM) pollution in such areas. Yet, windblown dust is poorly represented in existing regulatory air quality models. In a study by the authors on modeling episodic high PM events along the US/Mexico border using the state-of-the-art CMAQ/MM5/SMOKE air quality modeling system [Choi, Y.-J., Hyde, P., Fernando, H.J.S., 2006. Modeling of episodic particulate matter events using a 3D air quality model with fine grid: applications to a pair of cities in the US/Mexico border. Atmospheric Environment 40, 5181–5201], some of the observed PM₁₀ NAAQS exceedances were inferred as due to windblown dust, but the modeling system was incapable of dealing with time-dependent episodic dust entrainment during high wind periods. In this paper, a time-dependent entrainment parameterization for windblown dust is implemented in the CMAQ/MM5/SMOKE modeling system with the hope of improving PM predictions. An approach for realizing windblown dust emission flux for each grid cell over the study domain on an hourly basis, which accounts for the influence of factors such as soil moisture content, atmospheric stability and wind speed, is presented in detail. Comparison of model predictions with observational data taken at a pair of US/Mexico border towns shows a clear improvement of model performance upon implementation of the dust emission flux parameterization.     

The impacts of different kinds of dust events on PM10 pollution in northern China

In this study the frequencies of PM₁₀ (as key urban pollutant) in 14 key environmental protection cities in northern China were analyzed. It follows that the PM₁₀ concentration in the high-frequency period is higher with an extent 0.009–0.066 mg m⁻³ than in the low-frequency period of 2001–2002. Further the impacts of three kinds of dust events on the PM₁₀ concentration in four cities (Beijing, Hohhot, Xi’an and Lanzhou) were explored. The results showed that different kinds of dust events have different influences on variation of PM₁₀ concentration in these four cities. In Lanzhou and Hohhot, which are near the source areas of dust events, the contribution degree of these three dust events to the PM₁₀ is: floating dust>dust storm>blowing dust. Whereas, in Beijing and Xi’an situated in dust event passing areas, the mean value of PM₁₀ concentration is higher in blowing dust than in floating dust (no dust storm). In addition, the influences of dust events on PM₁₀ concentration are different in the cities on different dust event paths. In Beijing and Hohhot (on the northern path), the high PM₁₀ concentration is usually caused by blowing dust. But in both Lanzhou and Xi’an (on the western/northwestern path) the high PM₁₀ pollution concentration is usually caused by floating dust.     

Coupled infrared extinction spectra and size distribution measurements for several non-clay components of mineral dust aerosol (quartz, calcite, and dolomite)

Simultaneous size distributions and Fourier transform infrared (FTIR) extinction spectra have been measured for several representative components of mineral dust aerosol (quartz, calcite, and dolomite) in the fine particle size mode (D=0.1–1 μm). Optical constants drawn from the published literature have been used in combination with the experimentally determined size distributions to simulate the extinction spectra. In general, Mie theory does not accurately reproduce the peak position or band shape for the prominent IR resonance features in the 800–1600 cm⁻¹ spectral range. The resonance peaks in the Mie simulation are consistently blue shifted relative to the experimental spectra by 20–50 cm⁻¹. Spectral simulations, derived from a simple Rayleigh-based analytic theory for a “continuous distribution of ellipsoids” particle shape model, better reproduce the experimental spectra, despite the fact that the Rayleigh approximation is not strictly satisfied in these experiments. These results differ from our previous studies of particle shape effects in silicate clay mineral dust aerosols where a disk-shaped model for the particles was found to be more appropriate.     

Detection of Asian dust aerosols using meteorological satellite data and suspended particulate matter concentrations

The advection and dispersion of Asian dust events from China to the Pacific Ocean around Japan during 2000–2002 were investigated using the meteorological satellite data of NOAA/AVHRR and GMS-5/VISSR. Aerosol vapour index images, taking the brightness temperature difference between 11 and 12 μm, are very effective for monitoring the Asian dust phenomenon in the East Asia region, with their capacity for detection during the day or night. We discuss the dust events, focusing on the advection patterns shown in satellite images, which are classified into three types as ‘dry slot’, ‘high-pressure wedge’ and ‘travelling high’, based on synoptic patterns. The results are compared with suspended particulate matter concentrations measured at Japanese surface stations and with ground-based observations of Sakurajima volcano by a web camera system at Kagoshima in Kyushu, Japan. We found that the passage of cold fronts caused a rapid increase of suspended particulate matter (SPM) concentrations, which exceeded 100 μg m⁻³, and that deep low-pressure complexes strengthened the dust phenomenon. The ‘high-pressure wedge’ type is seen much more clearly in satellite images than the ‘travelling high’ type, but SPM concentrations and visibility were similar in both owing to the differences in the vertical distribution of the dust and in viewing conditions.     

Tritium and Cl as constraints on fast fracture flow and percolation flux in the unsaturated zone at Yucca Mountain

An analysis of tritium and ³⁶Cl data collected at Yucca Mountain, Nevada suggests that fracture flow may occur at high velocities through the thick unsaturated zone. The mechanisms and extent of this “fast flow” in fractures at Yucca Mountain are investigated with data analysis, mixing models and several one-dimensional modeling scenarios. The model results and data analysis provide evidence substantiating the weeps model [Gauthier, J.H., Wilson, M.L., Lauffer, F.C., 1992. Proceedings of the Third Annual International High-level Radioactive Waste Management Conference, vol. 1, Las Vegas, NV. American Nuclear Society, La Grange Park, IL, pp. 891–989] and suggest that fast flow in fractures with minimal fracture–matrix interaction may comprise a substantial proportion of the total infiltration through Yucca Mountain. Mixing calculations suggest that bomb-pulse tritium measurements, in general, represent the tail end of travel times for thermonuclear-test-era (bomb-pulse) infiltration. The data analysis shows that bomb-pulse tritium and ³⁶Cl measurements are correlated with discrete features such as horizontal fractures and areas where lateral flow may occur. The results presented here imply that fast flow in fractures may be ubiquitous at Yucca Mountain, occurring when transient infiltration (storms) generates flow in the connected fracture network.     

Transport and degradation of propylene glycol in the vadose zone: model development and sensitivity analysis

Transport and degradation of de-icing chemical (containing propylene glycol, PG) in the vadose zone were studied with a lysimeter experiment and a model, in which transient water flow, kinetic degradation of PG and soil chemistry were combined. The lysimeter experiment indicated that aerobic as well as anaerobic degradation occurs in the vadose zone. Therefore, the model included both types of degradation, which was made possible by assuming advection-controlled (mobile) and diffusion-controlled (immobile) zones. In the mobile zone, oxygen can be transported by diffusion in the gas phase. The immobile zone is always water-saturated, and oxygen only diffuses slowly in the water phase. Therefore, the model is designed in a way that the redox potential can decrease when PG is degraded, and thus, anaerobic degradation can occur. In our model, manganese oxide (MnO₂, which is present in the soil) and NO \(_{3}^{-}\) (applied to enhance biodegradation) can be used as electron acceptors for anaerobic degradation. The application of NO \(_{3}^{-}\) does not result in a lower leaching of PG nor in a slower depletion of MnO₂. The thickness of the snowcover influences the leached fraction of PG, as with a high infiltration rate, transport is fast, there is less time for degradation and thus more PG will leach. The model showed that, in this soil, the effect of the water flow dominates over the effect of the degradation parameters on the leaching at a 1-m depth.     

Mobilization of colloids in groundwater due to infiltration of water at a coal ash disposal site

We investigated groundwaters in the vicinity of a coal ash site near an electric generating station in the western U.S.A. The purpose of the study was to ascertain why fine particles or colloids appear in some subsurface water samples there. If such fine particles are merely introduced during bailing or pumping operations which suspend otherwise immobile soil colloids, we should exclude these particulate materials from the water samples before analysis intended to quantify what is moving through the aquifer. However, if the colloids were truly suspended and moving with the groundwater flow in situ, then we should includes their contribution to our assessment of the mobile loads.Application of very careful sampling techniques (slow pumping rates, no atmospheric exposure) did not cause the large quantities of colloids observed previously to disappear from well water in which they occured. Additionally, the same sampling procedures did not cause similar abundances of colloids to appear in waters collected from neighboring wells installed and developed in the same manner and in the same geologic strata. Thus we believe sampling artifacts do not explain the colloids' presence in the groundwater samples.On the other hand, the groundwater chemistry and the nature of the suspended colloids (size, composition) strongly suggest these fine particles were suspended and therefore moving with the groundwater flow. At wells exhibiting large amounts of suspended colloids (≈10–100 mg L⁻¹), the water was enriched in CO₂ and depleted in O₂ relative to nearby locations. The colloids were typically between 0.1 and 2 μm in size and were primarily silicates. These results suggest to us that, where infiltrating water is percolating through a site that has been mixed with coal ash, the secondary carbonate mineral in the soils are being dissolved; removal of this cementing carbonate phase may consequently release soil silicate colloids to be carried in the flowing water.Such processes may enhance contaminant transport in groundwater by augmenting the pollutant load moving in the groundwater, and increasing the permeability of the porous medium to pollutant infiltration with water water and/or rainwater.     

Impact of the Saharan dust outbreaks on the ambient levels of total suspended particles (TSP) in the marine boundary layer (MBL) of the Subtropical Eastern North Atlantic Ocean

Six years (1998–2003) of measurements of ambient air concentrations of total suspended particulate (TSP) measured at a rural background monitoring station in Tenerife (Canary Islands), the El Río station (ER, 28°08′35″N, 16°39′20″W, 500 m a.s.l.) were studied. African dust outbreaks were objectively identified using a new quantitative tool, called the African Index. This index indicates the percentage of time that an air mass remained over an African region at one of three possible height intervals of the lower troposphere. After identifying these episodes, a study of the background TSP levels at the ER station and of direct and indirect (those which cause vertical deposition of dust) African air mass intrusion impacts was performed. Taking into account both direct and indirect episodes, a total of 322 days of African dust intrusion were objectively identified (a mean of 54 episodes per year) in the period 1998–2003, some of them caused by “transition episodes” or “return African air masses”. A subjective method confirmed that 256 of these days were caused by direct impacts of African dust on the ER station. A mean TSP value of 21.6 μg m⁻³ was found at the station during this period. All the episodes occurred when the TSP concentration was >28.5 μg m⁻³. The TSP background (14 μg m⁻³) can be assumed to be representative of the MBL of the Eastern North Atlantic subtropical region. The highest number of dust gravitational settlement (or indirect) episodes occurs in summer, but the highest contribution of these episodes to the TSP levels is in March with a monthly mean TSP contribution of up to 30.5 μg m⁻³.     

The ratio between nitrogen oxides and carbon monoxide total emissions as precursors of tropospheric hydroxyl content evolution

The Main Geophysical Observatory 2D channel photochemical model is used to study the behavior of tropospheric OH within the 30–60°N zonal belt in relation to changing NO_X and CO emissions. The changes of tropospheric OH as a function of the contributions by NO_X and CO emissions during the period 1850–2050 are calculated. Our estimations show that the largest annual increment of total tropospheric OH within the belt considered occurs in the 1985–1995 period, about 0.27% yr⁻¹. Based on scenarios of tropospheric pollution emissions in the first half of 21st century, the total tropospheric OH content will increase more slowly, by 0.12–0.15% yr⁻¹. The maximum growth of OH concentration occurs close to air pollution locations—in the lower troposphere during 1850–1995 but in the upper troposphere in the 21st century when the NO_X source from subsonic aircraft increases faster than the surface source.     

Phosphorus characterization in sediments impacted by septic effluent at four sites in central Canada

Characterization of phosphorus (P) enriched solids was undertaken in the sediments below four mature septic system infiltration beds, where previous monitoring of phosphate (PO₄) concentrations in the groundwater had indicated that substantial retention of P was occurring in the vadose zone. At each site, zones of sediment P enrichment were identified by an acid extraction procedure. Acid extractable sediment P concentrations were found to be 2–5 times higher than background values, within narrow discrete zones generally 10–30 cm in thickness, located within one meter of the infiltration pipes. Back scattered electron images of the P enriched zones indicated that the P solids occurred as distinct authigenic grains (up to 300 μm diameter) and as grain coatings. Microprobe analyses indicated predominantly Fe–P in calcareous sediments (Cambridge and Langton) and Al–Fe–P in non-calcareous sediments (Muskoka and Harp Lake). Porewater analyses indicated that the zones of P accumulation were closely associated with zones of redox change characterized by the conversion of effluent NH₄⁺ to NO₃⁻. The data suggests that a substantial amount of the septic derived P is being attenuated by mineral precipitation reactions that occur rapidly after the effluent encounters subsurface sediments. Reductive dissolution of ferric (oxy)hydroxide minerals as a consequence of reducing environments near the infiltrations pipes, the release of Fe²⁺ in solution and subsequent conversion of Fe²⁺ to Fe³⁺ may promote the precipitation of ferric or ferrosoferric PO₄ minerals. In sediments with limited buffering capacity (calcite deficient), the decrease in pH resulting from effluent oxidation may cause Al (oxy)hydroxide dissolution and subsequent precipitation of Al–P rich phases. These precipitation reactions appear to have the capacity to immobilize a substantial amount of septic derived P (25–99% at these sites) for a considerable period of time.     

Source characterization of fine and coarse particles at the East Mediterranean coast

Fine and coarse atmospheric particles were collected in Ashdod—a midsize industrial city on the southeastern Mediterranean coast, and in Gedera—a rural site, to characterize ambient particles and to determine their long-range transport during two major seasons—winter and summer. Manual PM2.5 and PM10 samplers, dichotomous samplers, continuous automated PM10 samplers, and denuders were used to sample particulate and gaseous pollutants.Fine and coarse concentrations in Ashdod were 21.2 and 39.6 μg m⁻³, and 23.9 and 30.5 μg m⁻³ in the fall–winter and summer campaigns, respectively. Crustal material, as calcites or dolomites mixed with silicates, dominated the coarse fraction and also the fine fraction on dusty days. In the fall–winter, S, P, and Ni were coupled with minerals. Coarse Ni was associated with crustal material during dust storms, while P originated from shipping and deposition of phosphates in the urban area around.Sulfates dominated the fine fractions in the summer season averaging 12 μg m⁻³. Multivariate analysis indicated that S was associated with As and Se, V and Ni, both associated with heavy fuel combustion, and Zn and Pb. In winter, those mixed sources were local, but in summer they were part of long-range transport. In the fall–winter, Zn and Pb were strongly associated with Mn, Ga, and Cu—elements emitted from either traffic or metal processing plants.Although the influence of crustal material on both size fractions was significant, most heavy metals were associated with PM2.5. Higher concentrations were linked to a larger number of particles in this fraction, to a larger surface area available for biochemical reaction [Harrison, R., Shi, J., Xi, S., Khan, A., Mark, D., Kinnersley, R., Yin, J., Philos, T., 2000. Measurement of number, mass and size distribution of particles in the atmosphere. Philosophical Transactions of the Royal Society 358, 2567–2579], and finally to a larger concern in regards to health effects.     

Chemical compositions and radiative properties of dust and anthropogenic air masses study in Taipei Basin, Taiwan, during spring of 2004

Asia is one of the major sources of not only mineral dust but also anthropogenic aerosols. Continental air masses associated with the East Asian winter monsoon always contain high contents of mineral dust and anthropogenic species and transported southeastward to Taiwan, which have significant influences on global atmospheric radiation transfer directly by scattering and absorbing solar radiation in each spring. However, few measurements for the long-range transported aerosol and its optical properties were announced in this area, between the Western Pacific and the southeastern coast of Mainland China. The overall objective of this work is to quantify the optical characteristics of different aerosol types in the Eastern Asian. In order to achieve this objective, meteorological parameters, concentrations of PM₁₀ and its soluble species, and optical property of atmospheric scattering coefficients were measured continuously with 1 h time-resolved from 11 February to 7 April 2004 in Taipei Basin (25°00′N, 121°32′E). In this work, the dramatic changes of meteorological parameters such as temperature and winds were used to determine the influenced period of each air mass. Continental, strong continental, marine, and stagnant air masses defined by the back-trajectory analysis and local meteorology were further characterized as long-range transport pollution, dust, clean marine, and local pollution aerosols, respectively, according to the diagnostic ratios. The aerosol mass scattering efficiency of continental pollution, dust, clean marine, and local pollution aerosols were ranged from 1.3 to 1.6, 0.7 to 1.0, 1.4 and 1.4 to 2.3 m² g⁻¹, respectively. Overall, there are two distinct populations of aerosol mass scattering efficiencies, one for an aerosol chemical composition dominated by dust (<1.0 m² g⁻¹) and the other for an aerosol chemical composition dominated by anthropogenic pollutants (1.3–2.3 m² g⁻¹), which were similar to the previous measurements with high degree of temporal resolution.