Experimental and Numerical Study of Radioiodine Sorption and Transport in Hanford Sediments.

Radioiodine (129I) poses a potential risk to human health and the environment at several U.S. Department of Energy sites, including the Hanford Site, located in southeastern Washington State. Experimental studies and numerical modeling were performed to provide a technical basis for field-scale modeling of iodine sorption and transport behavior. The experiments were carried out using six columns of repacked contaminated sediments from the Hanford Site. Although iodate has been determined to be the dominant iodine species at the Hanford Site, the sorption and transport behaviors of different iodine species were investigated in a series of column experiments by first leaching sediments with artificial groundwater (AGW) followed by AGW containing iodate (IO3-), iodide (I-), or organo-iodine (2-iodo-5-methoxyphenol, C7H7IO2). Ferrihydrite amendments were added to the sediments for three of the columns to evaluate the impact of ferrihydrite on 129I attenuation. The results showed that ferrihydrite enhanced the iodate sorption capacity of the sediment and retarded the transport but had little effect on iodide or organo-I, providing a technical basis for developing a ferrihydrite-based remedial strategy for iodate under oxidizing conditions. Data from the column transport experiments were modeled using the linear equilibrium Freundlich isotherm model, the kinetic Langmuir adsorption model, and a distributed rate model. Comparisons of the experimental data and modeling results indicated that sorption was best represented with the distributed rate model with rates and maximum sorption extents varying by iodine species and ferrihydrite treatment. However, the linear Freundlich isotherm (Kd) model was also found to fit the laboratory experimental data relatively well, suggesting that the Kd model could also be used to represent iodine transport at the field scale.

A review of the behavior of radioiodine in the subsurface at two DOE sites.

Multiple processes affect the fate of the radioactive isotope 129I in the environment. Primary categories of these processes include electron transfer reactions mediated by minerals and microbes, adsorption to sediments, interactions with organic matter, co-precipitation, and volatilization. A description of dominant biogeochemical processes is provided to describe the interrelationship of these processes and the associated iodine chemical species. The majority of the subsurface iodine fate and transport studies in the United States have been conducted at U.S. Department of Energy (DOE) sites where radioisotopes of iodine are present in the environment and stored waste. The DOE Hanford Site and Savannah River Site (SRS) are used to illustrate how the iodine species and dominant processes at a site are controlled by the prevailing site biogeochemical conditions. These sites differ in terms of climate (arid vs. sub-tropical), major geochemical parameters (e.g., pH ~7.5 vs. 4), and mineralogy (carbonate vs. Fe/Al oxide dominated). The iodine speciation and dominant processes at a site also have implications for selection and implementation of suitable remedy approaches for 129I.

Relative permeability for water and gas through fractures in cement.

Relative permeability is an important attribute influencing subsurface multiphase flow. Characterization of relative permeability is necessary to support activities such as carbon sequestration, geothermal energy production, and oil and gas exploration. Previous research efforts have largely neglected the relative permeability of wellbore cement used to seal well bores where risks of leak are significant. Therefore this study was performed to evaluate fracturing on permeability and relative permeability of wellbore cement. Studies of relative permeability of water and air were conducted using ordinary Portland cement paste cylinders having fracture networks that exhibited a range of permeability values. The measured relative permeability was compared with three models, 1) Corey-curve, often used for modeling relative permeability in porous media, 2) X-curve, commonly used to represent relative permeability of fractures, and 3) Burdine model based on fitting the Brooks-Corey function to fracture saturation-pressure data inferred from x-ray computed tomography (XCT) derived aperture distribution results. Experimentally-determined aqueous relative permeability was best described by the Burdine model. Though water phase tended to follow the Corey-curve for the simple fracture system while air relative permeability was best described by the X-curve.

Quantifying community assembly processes and identifying features that impose them.

Spatial turnover in the composition of biological communities is governed by (ecological) Drift, Selection and Dispersal. Commonly applied statistical tools cannot quantitatively estimate these processes, nor identify abiotic features that impose these processes. For interrogation of subsurface microbial communities distributed across two geologically distinct formations of the unconfined aquifer underlying the Hanford Site in southeastern Washington State, we developed an analytical framework that advances ecological understanding in two primary ways. First, we quantitatively estimate influences of Drift, Selection and Dispersal. Second, ecological patterns are used to characterize measured and unmeasured abiotic variables that impose Selection or that result in low levels of Dispersal. We find that (i) Drift alone consistently governs ∼25% of spatial turnover in community composition; (ii) in deeper, finer-grained sediments, Selection is strong (governing ∼60% of turnover), being imposed by an unmeasured but spatially structured environmental variable; (iii) in shallower, coarser-grained sediments, Selection is weaker (governing ∼30% of turnover), being imposed by vertically and horizontally structured hydrological factors;(iv) low levels of Dispersal can govern nearly 30% of turnover and be caused primarily by spatial isolation resulting from limited exchange between finer and coarser-grain sediments; and (v) highly permeable sediments are associated with high levels of Dispersal that homogenize community composition and govern over 20% of turnover. We further show that our framework provides inferences that cannot be achieved using preexisting approaches, and suggest that their broad application will facilitate a unified understanding of microbial communities.

Persistence of uranium groundwater plumes: contrasting mechanisms at two DOE sites in the groundwater-river interaction zone.

We examine subsurface uranium (U) plumes at two U.S. Department of Energy sites that are located near large river systems and are influenced by groundwater-river hydrologic interaction. Following surface excavation of contaminated materials, both sites were projected to naturally flush remnant uranium contamination to levels below regulatory limits (e.g., 30 µg/L or 0.126 µmol/L; U.S. EPA drinking water standard), with 10 years projected for the Hanford 300 Area (Columbia River) and 12 years for the Rifle site (Colorado River). The rate of observed uranium decrease was much lower than expected at both sites. While uncertainty remains, a comparison of current understanding suggests that the two sites have common, but also different mechanisms controlling plume persistence. At the Hanford 300 A, the persistent source is adsorbed U(VI) in the vadose zone that is released to the aquifer during spring water table excursions. The release of U(VI) from the vadose zone and its transport within the oxic, coarse-textured aquifer sediments is dominated by kinetically-limited surface complexation. Modeling implies that annual plume discharge volumes to the Columbia River are small (

Stochastic simulation of uranium migration at the Hanford 300 Area.

This work focuses on the quantification of groundwater flow and subsequent U(VI) transport uncertainty due to heterogeneity in the sediment permeability at the Hanford 300 Area. U(VI) migration at the site is simulated with multiple realizations of stochastically-generated high resolution permeability fields and comparisons are made of cumulative water and U(VI) flux to the Columbia River. The massively parallel reactive flow and transport code PFLOTRAN is employed utilizing 40,960 processor cores on DOE's petascale Jaguar supercomputer to simultaneously execute 10 transient, variably-saturated groundwater flow and U(VI) transport simulations within 3D heterogeneous permeability fields using the code's multi-realization simulation capability. Simulation results demonstrate that the cumulative U(VI) flux to the Columbia River is less responsive to fine scale heterogeneity in permeability and more sensitive to the distribution of permeability within the river hyporheic zone and mean permeability of larger-scale geologic structures at the site.

Migration of saline solutions in variably saturated porous media.

Migration of concentrated NaNO3 solutions in homogeneous packs of pre-wetted silica sands was investigated using a light transmission system. Solutions of 5 molal NaNO3 were found to migrate downward 24-62% faster than pure water, in an unstable, fingered manner. This behavior was attributed primarily to a surface tension induced, non-zero apparent contact angle between the imbibing and the resident fluids. For saline solutions of similar surface tension to that of pure water (achieved by the addition of 2% methanol), the migration rates and plume shapes were comparable to that of water, demonstrating that density was not the primary source of the observed differences in migration patterns. At depths where resident saturation increased above residual, the migration process appeared to occur via film flow with slight changes in saturation (<4%), rather than in a series of abrupt jumps, as observed at shallower depths. A method for contact angle scaling was used to illustrate the effects of non-zero contact angles on capillary pressure-saturation curves.