Search

Minimizing loss of injected hydraulic fracturing fluids into shale along fracture-matrix boundaries is desired because imbibed water restricts gas production and wastes valuable water resources. This problem has motivated the addition of surfactants into water-based hydraulic fracturing fluids in order to reduce the capillary driving force for imbibition. Here, we show that reduction in interfacial tension and wettability alteration has negligible ability to reduce imbibition in deep gas reservoirs. The effectiveness of altering capillary forces acting at the wetting front also depends on the injection pressure acting at the fracture-matrix boundary. The pressure at the interface between the fracture and the shale matrix is constrained between the reservoir pore pressure and formation pressure (rock fracture pressure, also known as breakdown pressure of the rock) and increases with depth to magnitudes that greatly exceed that of capillary pressures. The analyses presented here show that even maximum alteration of interfacial properties that result in strongly hydrophobic interactions between the fracturing fluid and reservoir rock is incapable of significantly reducing imbibition in deep reservoirs. Instead of using surfactants, this analysis points to decreases in wellbore shut-in pressures and shut-in times as practical options for reducing imbibition losses of water-based fluids.

Atmospheric nitrous oxide contributes directly to global warming, yet models of the nitrogen cycle do not account for bedrock, the largest pool of terrestrial nitrogen, as a source of nitrous oxide. Although it is known that release rates of nitrogen from bedrock are large, there is an incomplete understanding of the connection between bedrock-hosted nitrogen and atmospheric nitrous oxide. Here, we quantify nitrogen fluxes and mass balances at a hillslope underlain by marine shale. We found that, at this site, bedrock weathering contributes 78% of the subsurface reactive nitrogen, while atmospheric sources (commonly regarded as the sole sources of reactive nitrogen in pristine environments) account for only the remaining 22%. About 56% of the total subsurface reactive nitrogen denitrifies, including 14% emitted as nitrous oxide. The remaining reactive nitrogen discharges in porewaters to a floodplain where additional denitrification probably occurs. We also found that the release of bedrock nitrogen occurs primarily within the zone of the seasonally fluctuating water table and suggest that the accumulation of nitrate in the vadose zone, often attributed to fertilization and soil leaching, may also include contributions from weathered nitrogen-rich bedrock. Our hillslope study suggests that, under oxygenated and moisture-rich conditions, weathering of deep, nitrogen-rich bedrock makes an important contribution to the nitrogen cycle. Weathering of deep bedrock releases reactive nitrogen into the subsurface, which contributes to the flux of nitrous oxide to the atmosphere, according to a field study that combines soil, rock and groundwater data within a river catchment.

The mixed-wet nature of reservoir formations imposes a wide range of rock wettability from strong resident-fluid wetting to strong invading-fluid wetting. The characteristics of two-phase flow in porous media composed of mixed-wetting surfaces remain poorly understood. In this study, we investigated the displacement of resident ethylene glycol (EG) by hexane in two mixed-wet micromodels of identical 2.5-D geometry heterogeneity, with uniformly or heterogeneously distributed patches strongly wetting to hexane. These patches are mixed among pores with unaltered EG-wetting surfaces. Along with control tests in the originally EG-wet micromodel, we show the classic fingering and transitions in flow regimes at (Formula presented.) (capillary number) from −7.2 to −3.9. Moreover, pore-scale distributions of wettability and their spatial correlation influence displacement efficiency. In the two mixed-wet micromodels, we found (a) an increase of steady-state hexane saturation at the end of experiments by up to 0.12 in the capillary fingering regime and a decrease of at most by 0.06 in the viscous fingering regime, compared to the EG-wet micromodel, and (b) dispersed and fragmented hexane distribution after displacement. Brine drainage during supercritical CO2 (scCO2) injections in these micromodels occurs with lower wettability contrasts, and under similar viscosity ratios and interfacial tensions resulted in higher displacement efficiency relative to displacement of EG by hexane. While mixed-wettability can enhance displacement efficiency compared to uniform wettability, the dynamics of immiscible fluids in strong mixed-wet reservoirs are expected to be less pronounced in contributing to the efficiency of geological CO2 sequestration, oil recovery, and remediation of hydrocarbon-contaminated aquifers.

The ability to predict multiphase fluid transport in nanoporous rocks such as shales is critical for many geoscience applications, for example unconventional hydrocarbon production, geologic carbon sequestration, and nuclear waste disposal. When the pore sizes approach nanoscales, the impact of the molecular interaction forces between fluids and solids becomes increasingly important. These forces can alter macroscopic fluid phase behavior and control transport. Recent experimental studies have shown that capillary condensation and subsequent imbibition of liquid water can occur in hydrophilic nanoporous media even if the vapor phase is at a critical relative humidity (rhcrit) well below vapor saturation. This study presents a theoretical investigation of the processes controlling adsorption, capillary condensation and imbibition in nanoporous media, using the square-gradient classical density functional theory. The proposed theoretical model explicitly includes the relevant interaction forces among fluids and solids in macroscopic porous media. Application of the model to a relative-humidity-controlled water adsorption experiment is presented to demonstrate the impact of water-pore wall attractive forces on multiphase water behavior in a hydrophilic silicon nanoporous medium. The model represents well the measured time-dependent evolution of the water imbibition front inside the nanoporous medium and also explains the diffusion-like water transport regimes observed at rhrhcrit. The study furthermore gives an insight on hysteresis phenomenon in adsorption and desorption isotherms.

Author(s): Tokunaga, TK; Wan, J; Williams, KH; Brown, W; Henderson, A; Kim, Y; Tran, AP; Conrad, ME; Bill, M; Carroll, RWH; Dong, W; Xu, Z; Lavy, A; Gilbert, B; Carrero, S; Christensen, JN; Faybishenko, B; Arora, B; Siirila-Woodburn, ER; Versteeg, R; Raberg, JH; Peterson, JE; Hubbard, SS | Abstract: Major components of hydrologic and elemental cycles reside underground, where their complex dynamics and linkages to surface waters are obscure. We delineated seasonal subsurface flow and transport dynamics along a hillslope in the Rocky Mountains (USA), where precipitation occurs primarily as winter snow and drainage discharges into the East River, a tributary of the Gunnison River. Hydraulic and geochemical measurements down to 10 m below ground surface supported application of transmissivity feedback of snowmelt to describe subsurface flow and transport through three zones: soil, weathering shale, and saturated fractured shale. Groundwater flow is predicted to depths of at least 176 m, although a shallower limit exists if hillslope-scale hydraulic conductivities are higher than our local measurements. Snowmelt during the high snowpack water year 2017 sustained flow along the weathering zone and downslope within the soil, while negligible downslope flow occurred along the soil during the low snowpack water year 2018. We introduce subsurface concentration-discharge (C-Q) relations for explaining hillslope contributions to C-Q observed in rivers and demonstrate their calculations based on transmissivity fluxes and measured pore water specific conductance and dissolved organic carbon. The specific conductance data show that major ions in the hillslope pore waters, primarily from the weathering and fractured shale, are about six times more concentrated than in the river, indicating hillslope solute loads are disproportionately high, while flow from this site and similar regions are relatively smaller. This methodology is applicable in different representative environments within snow-dominated watersheds for linking their subsurface exports to surface waters.

Multiphase flow phenomena in nanoporous media are encountered in many science and engineering applications. Shales, for example, possessing complex nanopore networks, have considerable importance as source rocks for unconventional oil and gas production and as low-permeability seals for geologic carbon sequestration or nuclear waste disposal. This study presents a theoretical investigation of the processes controlling adsorption, capillary condensation, and imbibition in such nanoporous media, with a particular focus on understanding the effects of fluid-fluid and fluid-pore wall interaction forces in the interconnected nanopore space. Building on a new theoretical framework, we developed a numerical model for the multiphase nanoporous flow and tested it against water vapor uptake measurements conducted on a shale core sample. The model, which is based on the density functional approach, explicitly includes the relevant interaction forces among fluids and solids while allowing for a continuum representation of the porous medium. The experimental data include gravimetrically measured mass changes in an initially dry core sample exposed to varying levels of relative humidity, starting with a low relative humidity (rh = 0.31) followed by a period of a higher relative humidity (rh = 0.81). During this process, water vapor uptake in the dry core is recorded as a function of time. Our model suggests that, under low rh conditions, the flow within the shale sample is controlled by adsorption- and diffusion-type processes. After increasing the rh to 0.81, the uptake of water vapor becomes more significant, and according to our model, this can be explained by capillary condensation followed by immiscible displacement in the core sample. It appears that strong fluid-pore wall attractive forces cause condensation near the inlet, which then induces water imbibition further into sample.

Author(s): Chang, Chun; Kneafsey, Timothy J; Wan, Jiamin; Tokunaga, Tetsu K; Nakagawa, Seiji

Author(s): Tokunaga, TK | Abstract: Predicting water imbibition into porous materials is important in a wide variety of fields, and hydraulic fracturing of low permeability hydrocarbon reservoirs has emerged as an application that is imposing a large water footprint. Reliable predictions of imbibition are needed to better manage water use, yet are challenging because of uncertainties in both the permeability and capillary pressure driving force. Here, this uncertainty is reduced through evaluating correlations between the permeability and the effective capillary pressure associated with the wetting front, Pc,f. These correlations allow elimination of Pc,f from the Green and Ampt equation and concentrate all uncertainties in fluxes on the effective permeability k. Over a wide range of k and porosities n, imbibition scales approximately with k1/3. Although Leverett k1/4 scaling for predicting Pc,f is shown to be inferior when tested with data spanning a wide range of n, it nevertheless predicted imbibition fairly well. From simple imbibition measurements, both the empirical and Leverett scaling approaches allow estimates of k that have root-mean-square deviations of about one order of magnitude relative to measurements that ranged over 10 orders of magnitude in k.

Author(s): Wan, J; Tokunaga, TK; Dong, W; Williams, KH; Kim, Y; Conrad, ME; Bill, M; Riley, WJ; Hubbard, SS | Abstract: Understanding terrestrial carbon cycling has relied primarily on studies of topsoils that are typically characterized to depths shallower than 0.5nm. At a semiarid site instrumented down to 7nm, we measured seasonal- and depth-resolved carbon inventories and fluxes and groundwater and unsaturated zone flow rates. Measurements showed that ~30% of the CO2 efflux to the atmosphere (60% in winter) originates from below 1nm, contrary to predictions of less than 1% by Earth System Model land modules. Respiration from deeper roots and deeper microbial communities is supported by favorable subsurface temperatures, moisture, and oxygen availability. Below 1nm, dissolved organic carbon fluxes from the overlying soil and C from deep roots and exudates are expected to be important in sustaining microbial respiration. Because these conditions are characteristic of semiarid climate regions, we contend that Earth System Model land modules should incorporate such deeper soil processes to improve CO2 flux predictions.

Author(s): Wang, L; Wan, J; Tokunaga, TK; Kim, Y; Yu, Q | Abstract: Shale gas equilibrates through gas-liquid-solid interactions in reservoirs, but the role of moisture is rarely investigated. To determine how adsorbed water influences methane behavior, three carboniferous shale samples from the Qaidam Basin, China, were humidified at five levels up to a relative humidity of 89%, and their methane capacities at pressures up to 12nMPa were studied. The experimental results indicate that two water-related mechanisms, “water blocking for methane transport” and “surface competition for gas-solid interaction,” are primarily responsible for the methane capacity variations. A compositional comparison suggests that a high abundance of clay minerals plays a favorable role in methane migration by retaining water in interlayer pores. Based on the experimental data, an optimized method for calculating the adsorption amount based on an approximation of density distribution is proposed. The model predicts the average thickness of the adsorption layer and the adsorbed methane density distribution on the surface at a given pressure. The methane adsorption layer “thins” in a stepped pattern by up to 45% in the presence of water, with little further change observed at relative humidities greater than 75% in the studied samples.

Interactions between supercritical (sc) CO2 and minerals are important when CO2 is injected into geologic formations for storage and as working fluids for enhanced oil recovery, hydraulic fracturing, and geothermal energy extraction. It has previously been shown that at the elevated pressures and temperatures of the deep subsurface, scCO2 alters smectites (typical swelling phyllosilicates). However, less is known about the effects of scCO2 on nonswelling phyllosilicates (illite and muscovite), despite the fact that the latter are the dominant clay minerals in deep subsurface shales and mudstones. Our studies conducted by using single crystals, combining reaction (incubation with scCO2), visualization [atomic force microscopy (AFM)], and quantifications (AFM, X-ray photoelectron spectroscopy, X-ray diffraction, and off-gassing measurements) revealed unexpectedly high CO2 uptake that far exceeded its macroscopic surface area. Results from different methods collectively suggest that CO2 partially entered the muscovite interlayers, although the pathways remain to be determined. We hypothesize that preferential dissolution at weaker surface defects and frayed edges allows CO2 to enter the interlayers under elevated pressure and temperature, rather than by diffusing solely from edges deeply into interlayers. This unexpected uptake of CO2, can increase CO2 storage capacity by up to ∼30% relative to the capacity associated with residual trapping in a 0.2-porosity sandstone reservoir containing up to 18 mass % of illite/muscovite. This excess CO2 uptake constitutes a previously unrecognized potential trapping mechanism.

Gas production from shale gas reservoirs plays a significant role in satisfying increasing energy demands. Compared with conventional sandstone and carbonate reservoirs, shale gas reservoirs are characterized by extremely low porosity, ultra-low permeability and high clay content. Slip flow, diffusion, adsorption and desorption are the primary gas transport processes in shale matrix, while Darcy flow is restricted to fractures. Understanding methane diffusion and adsorption, and gas flow and equilibrium in the low-permeability matrix of shale is crucial for shale formation evaluation and for predicting gas production. Modeling of diffusion in low-permeability shale rocks requires use of the Dusty gas model (DGM) rather than Fick’s law. The DGM is incorporated in the TOUGH2 module EOS7C-ECBM, a modified version of EOS7C that simulates multicomponent gas mixture transport in porous media. Also included in EOS7C-ECBM is the extended Langmuir model for adsorption and desorption of gases. In this study, a column shale model was constructed to simulate methane diffusion and adsorption through shale rocks. The process of binary $$\hbox {CH}_{4}{-}\hbox {N}_{2}$$ diffusion and adsorption was analyzed. A sensitivity study was performed to investigate the effects of pressure, temperature and permeability on diffusion and adsorption in shale rocks. The results show that methane gas diffusion and adsorption in shale is a slow process of dynamic equilibrium, which can be illustrated by the slope of a curve in $$\hbox {CH}_{4}$$ mass variation. The amount of adsorption increases with the pressure increase at the low pressure, and the mass change by gas diffusion will decrease due to the decrease in the compressibility factor of the gas. With the elevated temperature, the gas molecules move faster and then the greater gas diffusion rates make the process duration shorter. The gas diffusion rate decreases with the permeability decrease, and there is a limit of gas diffusion if the permeability is less than $$1.0\,\times \,10^{-15}\, \hbox { m}^{2}$$ . The results can provide insights for a better understanding of methane diffusion and adsorption in the shale rocks so as to optimize gas production performance of shale gas reservoirs.

Author(s): Hubbard, SS; Williams, KH; Agarwal, D; Banfield, J; Beller, H; Bouskill, N; Brodie, E; Carroll, R; Dafflon, B; Dwivedi, D; Falco, N; Faybishenko, B; Maxwell, R; Nico, P; Steefel, C; Steltzer, H; Tokunaga, T; Tran, PA; Wainwright, H; Varadharajan, C | Abstract: Extreme weather, fires, and land use and climate change are significantly reshaping interactions within watersheds throughout the world. Although hydrological–biogeochemical interactions within watersheds can impact many services valued by society, uncertainty associated with predicting hydrologydriven biogeochemical watershed dynamics remains high. With an aim to reduce this uncertainty, an approximately 300-km2 mountainous headwater observatory has been developed at the East River, CO, watershed of the Upper Colorado River Basin. The site is being used as a testbed for the Department of Energy supported Watershed Function Project and collaborative efforts. Building on insights gained from research at the “sister” Rifle, CO, site, coordinated studies are underway at the East River site to gain a predictive understanding of how the mountainous watershed retains and releases water, nutrients, carbon, and metals. In particular, the project is exploring how early snowmelt, drought, and other disturbances influence hydrological–biogeochemical watershed dynamics at seasonal to decadal timescales. A system-of-systems perspective and a scale-adaptive simulation approach, involving the combined use of archetypal watershed subsystem “intensive sites” are being tested at the site to inform aggregated watershed predictions of downgradient exports. Complementing intensive site hydrological, geochemical, geophysical, microbiological, geological, and vegetation datasets are long-term, distributed measurement stations and specialized experimental and observational campaigns. Several recent research advances provide insights about the intensive sites as well as aggregated watershed behavior. The East River “community testbed” is currently hosting scientists from more than 30 institutions to advance mountainous watershed methods and understanding.

How the wettability of pore surfaces affects supercritical (sc) CO2 capillary trapping in geologic carbon sequestration (GCS) is not well understood, and available evidence appears inconsistent. Using a high-pressure micromodel-microscopy system with image analysis, we studied the impact of wettability on scCO2 capillary trapping during short-term brine flooding (80 seconds, 8 to 667 pore volumes). Experiments on brine displacing scCO2 were conducted at 8.5 MPa and 45°C in water-wet (static contact angle θ = 20° ± 8°) and intermediate-wet (θ = 94° ± 13°) homogeneous micromodels under four different flow rates (capillary number Ca ranging from 9 × 10−6 to 8 × 10−4) with a total of eight conditions (four replicates for each). Brine invasion processes were recorded and statistical analysis was performed for over two thousand images of scCO2 saturations, and scCO2 cluster characteristics. The trapped scCO2 saturation under intermediate-wet conditions is 15% higher than under water-wet conditions under the slowest flow rate (Ca ∼ 9 × 10−6). Based on the visualization and scCO2 cluster analyses, we show that the scCO2 trapping process in our micromodels is governed by bypass trapping that is enhanced by the larger contact angle. Smaller contact angles enhance cooperative pore filling and widen brine fingers (or channels), leading to smaller volumes of scCO2 being bypassed. Increased flow rates suppress this wettability effect.

Wettability is an important factor controlling the displacement of immiscible fluids in porous media and therefore affects the flow and transport of supercritical (sc) CO2 in geologic carbon sequestration. Because few studies have focused on the wettability effect in scCO2-brine flow systems, we experimentally and numerically tested influences of wettability on drainage at the pore network scale. Using a high-pressure micromodel-microscopy system, we performed experiments of scCO2 invasion into brine-saturated water-wet and intermediate-wet micromodels, and recorded the scCO2 invasion morphology under reservoir relevant conditions. We also performed pore-scale numerical simulations to infer 3D details of fluid–fluid displacement processes. During drainage under intermediate-wet conditions, we found higher scCO2 saturation, wider scCO2 fingering, and more compact displacement patterns. The simulation results are qualitatively consistent with the experiments. Through quantitative analyses of the experiments, we found that the reduced wettability decreases the displacement front velocity, promotes non-wetting phase pore-filling events in the longitudinal direction, delays the breakthrough time of invading fluid, and increases the displacement efficiency. Simulated results also show that the fluid–fluid interface area follows a unified power-law relation with scCO2 saturation, and show smaller interface area in intermediate-wet case which suppresses the mass transfer between the phases. These pore-scale results provide insights for the wettability effects on CO2–brine immiscible displacement in geologic carbon sequestration.

© 2017 American Chemical Society. Three-dimensional variably saturated flow and multicomponent biogeochemical reactive transport modeling, based on published and newly generated data, is used to better understand the interplay of hydrology, geochemistry, and biology controlling the cycling of carbon, nitrogen, oxygen, iron, sulfur, and uranium in a shallow floodplain. In this system, aerobic respiration generally maintains anoxic groundwater below an oxic vadose zone until seasonal snowmelt-driven water table peaking transports dissolved oxygen (DO) and nitrate from the vadose zone into the alluvial aquifer. The response to this perturbation is localized due to distinct physico-biogeochemical environments and relatively long time scales for transport through the floodplain aquifer and vadose zone. Naturally reduced zones (NRZs) containing sediments higher in organic matter, iron sulfides, and non-crystalline U(IV) rapidly consume DO and nitrate to maintain anoxic conditions, yielding Fe(II) from FeS oxidative dissolution, nitrite from denitrification, and U(VI) from nitrite-promoted U(IV) oxidation. Redox cycling is a key factor for sustaining the observed aquifer behaviors despite continuous oxygen influx and the annual hydrologically induced oxidation event. Depth-dependent activity of fermenters, aerobes, nitrate reducers, sulfate reducers, and chemolithoautotrophs (e.g., oxidizing Fe(II), S compounds, and ammonium) is linked to the presence of DO, which has higher concentrations near the water table. (Figure Presented).

Although bedrock weathering strongly influences water quality and global carbon and nitrogen budgets, the weathering depths and rates within subsurface are not well understood nor predictable. Determination of both porewater chemistry and subsurface water flow are needed in order to develop more complete understanding and obtain weathering rates. In a long-term field study, we applied a multiphase approach along a mountainous watershed hillslope transect underlain by marine shale. Here we report three findings. First, the deepest extent of the water table determines the weathering front, and the range of annually water table oscillations determines the thickness of the weathering zone. Below the lowest water table, permanently water-saturated bedrock remains reducing, preventing deeper pyrite oxidation. Secondly, carbonate minerals and potentially rock organic matter share the same weathering front depth with pyrite, contrary to models where weathering fronts are stratified. Thirdly, the measurements-based weathering rates from subsurface shale are high, amounting to base cation exports of about 70 kmolc ha−1 y−1, yet consistent with weathering of marine shale. Finally, by integrating geochemical and hydrological data we present a new conceptual model that can be applied in other settings to predict weathering and water quality responses to climate change.

SummaryWatersheds are important for supplying fresh water, the quality of which depends on complex interplay involving physical, chemical and biological processes. As water percolates through the soil and underlying weathering rock en route to the river corridor, microorganisms mediate key geochemical transformations, yet the distribution and functional capacities of subsurface microbial communities remain little understood. We have studied metabolic capacities of microbial communities along a meadow to floodplain hillslope transect within the East-River watershed, Colorado, using genome resolved metagenomics and carbon and hydrogen stable isotopes. Very limited strain/species overlap was found at different depths below the ground surface and at different distances along the hillslope, possibly due to restricted hydraulic connectivity after early stages of snowmelt. Functions such as carbon fixation and selenate reduction were prevalent at multiple sites, although the lineages of organisms responsible tend to be location-specific. Based on its abundance, sulfur is significantly more important for microbial metabolism at the floodplain compared to on the hillslope. Nitrification and methylamine oxidation are likely only occurring within the floodplain, with nitrification capacity in shallow soil, and methylamine oxidation in deeper unsaturated sediment. Biogenic methane was detected in deep surface samples, but methanogenic organisms were not identified.Originality-Significance StatementIn a previous study within a hillslope to riparian zone transect of a sub-alpine watershed, the community structure was explored using ribosomal protein S3 genes, and the metabolic potential was hypothesized based on the presence of metabolism related genes. However, tying specific strains and species to metabolic functioning was not discussed as resolved genomes were not available.In the current study, we use genome-resolved metagenomics along with carbon and hydrogen stable isotopes to explore the spatial distribution of biogeochemical processes. By linking taxonomy and function, using multiple functional genes indicative of full metabolic pathways, we detect heterogeneity in the distribution of metabolic potential and the organisms involved with depth and landscape position. Thus, we infer how microbiome genomic variation impacts biogeochemical cycling across the watershed.We found very limited strain/species overlap at different depths below the surface and along the hillslope, possibly due to the restricted site to site hydraulic connectivity, and show that communities are largely distinct in their metabolic capacities. Both proximity to the river and the underlying Mancos shale apparently control species distribution and metabolic potential.Functions such as carbon fixation and selenate reduction were prevalent at multiple sites, although the lineages of organisms responsible tend to be location-specific. Arsenate detoxification was found to be prevalent in the riparian zone whereas selenate reduction was detected within weathered Mancos shale. We conclude that important ecosystem functions are strongly associated with the riparian zone, some of which may have crucial implications as to water quality and human health.

Capillary pressure (Pc) – saturation (Sw) relations are essential for predicting equilibrium and flow of immiscible fluid pairs in soils and deeper geologic formations. In systems that are difficult to measure, behavior is often estimated based on capillary scaling of easily measured Pc–Sw relations (e.g., air-water, and oil-water), yet the reliability of such approximations needs to be examined. In this study, seventeen sets of brine drainage and imbibition curves were measured with air-brine, decane-brine, and supercritical (sc) CO2-brine in homogeneous quartz and carbonate sands, using porous plate systems under ambient (0.1 MPa, 23 ˚C) and reservoir (12.0 MPa, 45 ˚C) conditions. Comparisons between these measurements showed significant differences in residual nonwetting phase saturation, Snw,r. Through applying capillary scaling, changes in interfacial properties were indicated, particularly wettability. With respect to the residual trapping of the nonwetting phases, Snwr, CO2 > Snwr, decane > Snwr, air. Decane-brine and scCO2-brine Pc–Sw curves deviated significantly from predictions assuming hydrophilic interactions. Moreover, neither the scaled capillary behavior nor Snw,r for scCO2-brine were well represented by decane-brine, apparently because of differences in wettability and viscosities, indicating limitations for using decane (and other organic liquids) as a surrogate fluid in studies intended to apply to geological carbon sequestration. Thus, challenges remain in applying scaling for predicting capillary trapping and multiphase displacement processes across such diverse fields as vadose zone hydrology, enhanced oil recovery, and geologic carbon sequestration. This article is protected by copyright. All rights reserved.

CO 2 flooding is a proven enhanced oil recovery (EOR) technique and is also considered as a potential method for CO 2 sequestration. Despite having successful field trials on CO 2 EOR, the effects of reservoir wettability and permeability heterogeneity on the efficiency of miscible CO 2 flooding are not well understood. In this work, laboratory investigations have been carried out to evaluate the influence of these properties on the miscible CO 2 EOR performance. The wettability of hydrophilic Berea core samples was altered to be oil-wet by vacuum saturation of the clean and dry core samples with n-hexadecane. The permeability heterogeneity was obtained by combining two half pieces of axially split water-wet core samples of different permeabilities. Core flooding experiments were conducted for n-hexadecane – synthetic brine – CO 2 systems at 1400 psig backpressure to achieve minimum miscibility pressure (MMP) of CO 2 in n-hexadecane at the test temperature (24 ± 1 °C). It was found that wettability strongly influences CO 2 EOR. For the alternate cases of previously brine flooded (to remaining oil saturation) oil-wet and water-wet core samples, five pore volumes (PVs) of CO 2 recovered 100% and only 43% of remaining oil in place (ROIP) respectively. Three PVs of CO 2 could recover only about 0–5% ROIP from the split core samples. The mechanisms underlying these results are discussed. This study sheds light on the significant influence of reservoir wettability and permeability heterogeneity on the performance of miscible CO 2 EOR.

Recharge of alluvial aquifers is a key component in understanding the interaction between floodplain vadose zone biogeochemistry and groundwater quality. The Rifle Site (a former U-mill tailings site) adjacent to the Colorado River is a well-established field laboratory that has been used for over a decade for the study of biogeochemical processes in the vadose zone and aquifer. This site is considered an exemplar of both a riparian floodplain in a semiarid region and a post-remediation U-tailings site. In this paper we present Sr isotopic data for groundwater and vadose zone porewater samples collected in May and July 2013 to build a mixing model for the fractional contribution of vadose zone porewater (i.e. recharge) to the aquifer and its variation across the site. The vadose zone porewater contribution to the aquifer ranged systematically from 0% to 38% and appears to be controlled largely by the microtopography of the site. The area-weighted average contribution across the site was 8% corresponding to a net recharge of 7.5 cm. Given a groundwater transport time across the site of ~1.5 to 3 years, this translates to a recharge rate between 5 and 2.5 cm/yr, and with the average precipitation to the site implies a loss from the vadose zone due to evapotranspiration of 83% to 92%, both ranges are in good agreement with previously published results by independent methods. A uranium isotopic (234U/238U activity ratios) mixing model for groundwater and surface water samples indicates that a ditch across the site is hydraulically connected to the aquifer and locally significantly affects groundwater. Groundwater samples with high U concentrations attributed to natural bio-reduced zones have 234U/238U activity ratios near 1, suggesting that the U currently being released to the aquifer originated from the former U-mill tailings.

Author(s): Tokunaga, TK; Kim, Y; Wan, J; Bill, M; Conrad, M; Dong, W | Abstract: Transport from the soil surface to groundwater is commonly mediated through deeper portions of the vadose zone and capillary fringe, where variations in temperature and water saturation strongly influence biogeochemical processes. This technical note describes a sediment column design that allows laboratory simulation of thermal and hydrologic conditions found in many field settings. Temperature control is particularly important because room temperature is not representative of most subsurface environments. A 2.0-m-tall column was capable of simulating profiles with temperatures ranging from 3 to 22°C, encompassing the full range of seasonal temperature variation observed in the deep vadose zone and capillary fringe of a semiarid floodplain in western Colorado. The water table was varied within the lower 0.8-m section of the column, and profiles of water content and matric potential were measured. Vadose zone CO2 collected from depth-distributed gas samplers under representative seasonal conditions reflected the influences of temperature and water table depth on microbial respiration. Thus, realistic subsurface biogeochemical dynamics can be simulated in the laboratory through establishing column profiles that represent seasonal thermal and hydrologic conditions.

Author(s): Cihan, A; Wang, S; Tokunaga, TK; Birkholzer, JT | Abstract: Understanding the main mechanisms affecting long-term migration and redistribution of injected CO2 in geological carbon storage is needed for developing predictive models to assess environmental risks and designing monitoring schemes. Preparation of a postinjection site care plan is required for CO2 injection wells, including monitoring of pressure changes and injected CO2 plume. Knowledge gaps exist regarding assessment of postinjection monitoring timeframes for the CO2 plume because the processes driving long-term CO2 plume migration and trapping are not fully understood. In the postinjection stage of geological carbon storage, redistribution of CO2 occurs mainly due to buoyancy and capillary forces. This work presents experimental and modeling studies to investigate processes contributing to postinjection plume distribution and stabilization. We conducted a flow cell experiments (0.5nmn×n0.05nmn×n0.01nm) with two immiscible fluid phases in a glass bead porous medium to study postinjection plume behavior. We employed a hysteretic macroscopic two-phase flow model to interpret the experimental results and to understand main processes leading to plume stabilization. Our findings show that capillary pressure hysteresis explains the experimentally observed plume shape and redistribution at early postinjection stages; however, the long-term plume migration and eventual plume stabilization can only be represented when in addition microscale heterogeneity is accounted for. Results also show that plume stabilization can be extremely slow and that the migration of the plume front can occur through multiple intermittent bursts over long times. Further studies are needed to understand implications of the results for more realistic porous media and large-scale storage reservoirs.

Author(s): Hu, R; Wan, J; Yang, Z; Chen, YF; Tokunaga, T | Abstract: When a more viscous fluid displaces a less viscous one in porous media, viscous pressure drop stabilizes the displacement front against capillary pressure fluctuation. For this favorable viscous ratio conditions, previous studies focused on the front instability under slow flow conditions but did not address competing effects of wettability and flow rate. Here we study how this competition controls displacement patterns. We propose a theoretical model that describes the crossover from fingering to stable flow as a function of invading fluid contact angle θ and capillary number Ca. The phase diagram predicted by the model shows that decreasing θ stabilizes the displacement for θ≥45° and the critical contact angle θc increases with Ca. The boundary between corner flow and cooperative filling for θ l 45° is also described. This work extends the classic phase diagram and has potential applications in predicting CO2 capillary trapping and manipulating wettability to enhance gas/oil displacement efficiency.

Significance Reliable estimates of geologic carbon storage capacities (needed for policymaking) in both saline aquifers and unconventional gas/oil shales rely on understanding trapping mechanisms. We found that CO2 uptake by muscovite (a common mineral and a conservative proxy for illite) far exceeds the maximum adsorption capacity of its external surface area. Our measurements using different methods collectively suggest that CO2 enters muscovite interlayers without bulk interlayer expansion, contrary to the conventional wisdom that only swelling clays take up CO2 into interlayers. Because the nonswelling illitic clay is the major clay mineral in deep subsurface tight rocks, their excess uptake of CO2 may significantly contribute to CO2 storage capacity and warrants further in-depth studies., Interactions between supercritical (sc) CO2 and minerals are important when CO2 is injected into geologic formations for storage and as working fluids for enhanced oil recovery, hydraulic fracturing, and geothermal energy extraction. It has previously been shown that at the elevated pressures and temperatures of the deep subsurface, scCO2 alters smectites (typical swelling phyllosilicates). However, less is known about the effects of scCO2 on nonswelling phyllosilicates (illite and muscovite), despite the fact that the latter are the dominant clay minerals in deep subsurface shales and mudstones. Our studies conducted by using single crystals, combining reaction (incubation with scCO2), visualization [atomic force microscopy (AFM)], and quantifications (AFM, X-ray photoelectron spectroscopy, X-ray diffraction, and off-gassing measurements) revealed unexpectedly high CO2 uptake that far exceeded its macroscopic surface area. Results from different methods collectively suggest that CO2 partially entered the muscovite interlayers, although the pathways remain to be determined. We hypothesize that preferential dissolution at weaker surface defects and frayed edges allows CO2 to enter the interlayers under elevated pressure and temperature, rather than by diffusing solely from edges deeply into interlayers. This unexpected uptake of CO2, can increase CO2 storage capacity by up to ∼30% relative to the capacity associated with residual trapping in a 0.2-porosity sandstone reservoir containing up to 18 mass % of illite/muscovite. This excess CO2 uptake constitutes a previously unrecognized potential trapping mechanism.

It is commonly true that a diluted colloidal suspension is more stable over time than a concentrated one because dilution reduces collision rates of the particles and therefore delays the formation of aggregates. However, this generalization does not apply for some engineered ligand-coated nanoparticles (NPs). We observed the opposite relationship between stability and concentration of NPs. We tested 4 different types of NPs: CdSe-11-mercaptoundecanoic acid, CdTe-polyelectrolytes, Ag-citrate, and Ag-polyvinylpyrrolidone. The results showed that dilution alone induced aggregation and subsequent sedimentation of the NPs that were originally monodispersed at very high concentrations. Increased dilution caused NPs to progressively become unstable in the suspensions. The extent of the dilution impact on the stability of NPs is different for different types of NPs. We hypothesize that the unavoidable decrease in free ligand concentration in the aqueous phase following dilution causes detachment of ligands from the suspended NP cores. The ligands attached to NP core surfaces must generally approach exchange equilibrium with free ligands in the aqueous phase; therefore, ligand detachment and destabilization are expected consequences of dilution. More studies are necessary to test this hypothesis. Because the stability of NPs determines their physicochemical and kinetic behavior including toxicity, dilution-induced instability needs to be understood to realistically predict the behavior of engineered ligand-coated NPs in aqueous systems. Environ Toxicol Chem 2018;37:1301-1308. © 2018 SETAC.

© 2017 Elsevier Inc. When a wetting liquid is displaced by air in a capillary tube, a wetting film develops between the tube wall and the air that is responsible for the snap-off mechanism of the gas phase. By dissolving a dye in the wetting phase it is possible to relate a measure of the absorbance in the capillary to the thickness of liquid films. These data could be used to compare with cutting edge numerical simulations of the dynamics of snap-off for which experimental and numerical data are lacking. Drainage experiments in constricted capillary tubes were performed where a dyed wetting liquid is displaced by air for varying flow rates. We developed an optical method to measure liquid film thicknesses that range from 3 to 1000 μm. The optical measures are validated by comparison with both theory and direct numerical simulations. In a constricted capillary tube we observed, both experimentally and numerically, a phenomenon of snap-off coalescence events in the vicinity of the constriction that bring new insights into our understanding and modeling of two-phase flows. In addition, the good agreement between experiments and numerical simulations gives confidence to use the numerical method for more complex geometries in the future.

In geologic carbon sequestration, capillary pressure (Pc)-saturation (Sw) relations are needed to predict reservoir processes. Capillarity and its hysteresis have been extensively studied in oil-water and gas-water systems, but few measurements have been reported for supercritical (sc) CO2-water. Here, Pc-Sw relations of scCO2 displacing brine (drainage), and brine rewetting (imbibition) were studied to understand CO2 transport and trapping behavior under reservoir conditions. Hysteretic drainage and imbibition Pc-Sw curves were measured in limestone sands at 45 °C under elevated pressures (8.5 and 12.0 MPa) for scCO2-brine, and in limestone and dolomite sands at 23 °C (0.1 MPa) for air-brine using a new computer programmed porous plate apparatus. scCO2-brine drainage and imbibition curves shifted to lower Pc relative to predictions based on interfacial tension, and therefore deviated from capillary scaling predictions for hydrophilic interactions. Fitting universal scaled drainage and imbibition curves show that wettability alteration resulted from scCO2 exposure over the course of months-long experiments. Residual trapping of the nonwetting phases was determined at Pc = 0 during imbibition. Amounts of trapped scCO2 were significantly larger than for those for air, and increased with pressure (depth), initial scCO2 saturation, and time. These results have important implications for scCO2 distribution, trapping, and leakage potential.

Renewed interest in managing C balance in soils is motivated by increasing atmospheric concentrations of CO2 and consequent climate change. Here, experiments were conducted in soil columns to determine C mass balances with and without addition of CaSO4-minerals (anhydrite and gypsum), which were hypothesized to promote soil organic carbon (SOC) retention and soil inorganic carbon (SIC) precipitation as calcite under slightly alkaline conditions. Changes in C contents in three phases (gas, liquid and solid) were measured in unsaturated soil columns tested for one year and comprehensive C mass balances were determined. The tested soil columns had no C inputs, and only C utilization by microbial activity and C transformations were assumed in the C chemistry. The measurements showed that changes in C inventories occurred through two processes, SOC loss and SIC gain. However, the measured SOC losses in the treated columns were lower than their corresponding control columns, indicating that the amendments promoted SOC retention. The SOC losses resulted mostly from microbial respiration and loss of CO2 to the atmosphere rather than from chemical leaching. Microbial oxidation of SOC appears to have been suppressed by increased Ca(2+) and SO4(2)(-) from dissolution of CaSO4 minerals. For the conditions tested, SIC accumulation per m(2) soil area under CaSO4-treatment ranged from 130 to 260 g C m(-1) infiltrated water (20-120 g C m(-1) infiltrated water as net C benefit). These results demonstrate the potential for increasing C sequestration in slightly alkaline soils via CaSO4-treatment.

The current ambiguity on wettability of minerals in CO2–brine systems under the geological CO2 sequestration (GCS) reservoir conditions imparts the greatest uncertainty in predicting capillary behavior controlling safe-storage of CO2. To address this issue we conducted a series of experiments using muscovite as a representative of common aluminosilicate minerals. Based on new experimental results we identified several possible causes of the ambiguity problem in contact angle (CA) measurements. We also found that reaction with water-saturated supercritical (sc) CO2 (but not with scN2) phase severely roughened the muscovite surfaces, largely increased CA hysteresis and CO2 adhesion. Although some methodological influences on contact angle uncertainty can be reduced, the high surface-energy of clean and pristine aluminosilicate minerals have strong tendency to adsorb oppositely changed molecules and particles to reduce their surface energy, resulting in less reproducible CA values. Giving the fact that such clean and pristine mineral surfaces do not exist in real reservoirs, our future investigations shall focus on improving understanding of the effect of long-term CO2–mineral–brine reactions on reservoir wettability under realistic reservoir geochemical conditions.

© 2017 American Chemical Society. Environmentally benign, economical, and effective surfactants and additives are needed in engineering subsurface energy recovery processes. Biosurfactants have some advantages over chemically synthesized surfactants, but the high costs of microbial-biosynthesis limit their applications in subsurface engineering. Here we propose to use naturally occurring biosurfactants contained within Earth's readily available and inexpensive humus deposits for subsurface engineering applications. We collected humus samples of different types from four different regions, and developed a simple method for extracting natural biosurfactant (NBS) using only four common chemicals. The average NBS extraction yields are 16 ± 3% of the raw humus. No significant differences in elemental composition and functional group chemistry were found among the NBS extracted from humus of different origins, suggesting that any humus deposit can be used for NBS extraction. Measurements of interfacial tensions between air-water and supercritical (sc) CO2-water interfaces indicate that the NBS is a highly effective surfactant. NBS has good foaming ability. Preliminary tests with only 0.5 mass % NBS in the aqueous phase (no other additives) yielded scCO2-in-water foams of 83% foam quality. The apparent viscosity of 13 cP measured at a shear rate of 2320 s-1indicates that much higher viscosities are achievable at lower shear rates. These results suggest that NBS merits further research and development for potential applications in subsurface energy production.

© 2017 American Chemical Society. Diffusion is important in controlling local solute transport and reactions in unsaturated soils and geologic formations. Although it is commonly assumed that thinning of water films controls solute diffusion at low water contents, transport under these conditions is not well understood. We conducted experiments in quartz sands at low volumetric water contents (θ) to quantify ion diffusion within adsorbed films. At the lowest water contents, we employed fixed relative humidities to control water films at nm thicknesses. Diffusion profiles for Rb+ and Br- in unsaturated sand packs were measured with a synchrotron X-ray microprobe, and inverse modeling was used to determine effective diffusion coefficients, De, as low as ∼9 × 10-15 m2 s-1 at θ = 1.0 × 10-4 m3 m-3, where the film thickness = 0.9 nm. Given that the diffusion coefficients (Do) of Rb+ and Br- in bulk water (30 °C) are both ∼2.4 × 10-9 m2 s-1, we found the impedance factor f = De/(θDo) is equal to 0.03 ± 0.02 at this very low saturation, in agreement with the predicted influence of interface tortuosity (τa) for diffusion along grain surfaces. Thus, reduced cross-sectional area (θ) and tortuosity largely accounted for the more than 5 orders of magnitude decrease in De relative to Do as desaturation progressed down to nanoscale films.

Improving our ability to estimate the parameters that control water and heat fluxes in the shallow subsurface is particularly important due to their strong control on recharge, evaporation and biogeochemical processes. The objectives of this study are to develop and test a new inversion scheme to simultaneously estimate subsurface hydrological, thermal and petrophysical parameters using hydrological, thermal and electrical resistivity tomography (ERT) data. The inversion scheme – which is based on a nonisothermal, multiphase hydrological model – provides the desired subsurface property estimates in high spatiotemporal resolution. A particularly novel aspect of the inversion scheme is the explicit incorporation of the dependence of the subsurface electrical resistivity on both moisture and temperature. The scheme was applied to synthetic case studies, as well as to real datasets that were autonomously collected at a biogeochemical field study site in Rifle, Colorado. At the Rifle site, the coupled hydrological-thermal-geophysical inversion approach well predicted the matric potential, temperature and apparent resistivity with the Nash–Sutcliffe efficiency criterion greater than 0.92. Synthetic studies found that neglecting the subsurface temperature variability, and its effect on the electrical resistivity in the hydrogeophysical inversion, may lead to an incorrect estimation of the hydrological parameters. The approach is expected to be especially useful for the increasing number of studies that are taking advantage of autonomously collected ERT and soil measurements to explore complex terrestrial system dynamics.

© 2016 In order to understand the transport and humification processes of dissolved organic matter (DOM) within sediments of a semi-arid floodplain at Rifle, Colorado, fluorescence excitation–emission matrix (EEM) spectroscopy, humification index (HIX) and specific UV absorbance (SUVA) at 254 nm were applied for characterizing depth and seasonal variations of DOM composition. Results revealed that late spring snowmelt leached relatively fresh DOM from plant residue and soil organic matter down into the deeper vadose zone (VZ). More humified DOM is preferentially adsorbed by upper VZ sediments, while non- or less-humified DOM was transported into the deeper VZ. Interestingly, DOM at all depths undergoes rapid biological humification process evidenced by the products of microbial by-product-like (i.e., tyrosine-like and tryptophan-like) matter in late spring and early summer, particularly in the deeper VZ, resulting in more humified DOM (e.g., fulvic-acid-like and humic-acid-like substances) at the end of year. This indicates that DOM transport is dominated by spring snowmelt, and DOM humification is controlled by microbial degradation, with seasonal variations. It is expected that these relatively simple spectroscopic measurements (e.g., EEM spectroscopy, HIX and SUVA) applied to depth- and temporally-distributed pore-water samples can provide useful insights into transport and humification of DOM in other subsurface environments as well.

© 2016, Springer-Verlag Berlin Heidelberg. Technologies for increasing carbon storage in soils are gathering attention as a means for mitigating atmospheric CO2emissions. Carbon sequestration can be achieved by controlling the organic carbon stock in soil and by accelerating mineral carbonation. In this study, carbon sequestration capacity was measured in soil columns treated with flue gas desulfurization gypsum (FGDG), a by-product of electric power generation. The feasibility of using FGDG as an environmentally benign alternative to gypsum or anhydrite was examined using a toxicity characteristic leaching procedure and Microtox bioassay. While no toxic leachate was generated from the FGDG treatment, some toxic elements in the soil were removed through absorption reactions. Test results for carbon sequestration based on unsaturated soil column experiments suggest that the application of FGDG for soil treatment holds promise of less microbial CO2emission from soil. The net benefits of carbon sequestration from the FGDG treatment were calculated as 87 and 621 g C/m2/m of infiltrated water, for the 1 % calcite-added column and 3 % calcite-added columns, respectively. The presented test results show that the FGDG treatment for soil carbon sequestration holds a promise when it is applied to slightly alkaline soils.

© Soil Science Society of America. Although CO2fluxes from soils are often assumed to originate within shallow soil horizons (

© 2016, US Government. Flood plains play a potentially important role in the global carbon cycle. The accumulation of organic matter in flood plains often induces the formation of chemically reduced groundwater and sediments along riverbanks. In this study, our objective is to evaluate the cumulative impact of such reduced zones, water table fluctuations, and temperature gradients on subsurface carbon fluxes in a flood plain at Rifle, Colorado located along the Colorado River. 2-D coupled variably-saturated, non-isothermal flow and biogeochemical reactive transport modeling was applied to improve our understanding of the abiotic and microbially mediated reactions controlling carbon dynamics at the Rifle site. Model simulations considering only abiotic reactions (thus ignoring microbial reactions) underestimated CO2partial pressures observed in the unsaturated zone and severely underestimated inorganic (and overestimated organic) carbon fluxes to the river compared to simulations with biotic pathways. Both model simulations and field observations highlighted the need to include microbial contributions from chemolithoautotrophic processes (e.g., Fe+2and S−2oxidation) to match locally-observed high CO2concentrations above reduced zones. Observed seasonal variations in CO2concentrations in the unsaturated zone could not be reproduced without incorporating temperature gradients in the simulations. Incorporating temperature fluctuations resulted in an increase in the annual groundwater carbon fluxes to the river by 170 % to 3.3 g m−2d−1, while including water table variations resulted in an overall decrease in the simulated fluxes. We conclude that spatial microbial and redox zonation as well as temporal fluctuations of temperature and water table depth contribute significantly to subsurface carbon fluxes in flood plains and need to be represented appropriately in model simulations.

At the Savannah River Site's F-Area, wastewaters containing radionuclides were disposed into seepage basins for decades. After closure and capping in 1991, the U.S. Department of Energy (DOE) has being monitoring and remediating the groundwater plume. Despite numerous studies of the plume, its persistence for over 20 years has not been well understood. To better understand the plume dynamics, a limited number of deep boreholes were drilled to determine the current plume characteristics. A mixing model was developed to predict plume tritium and nitrate concentrations. We found that the plume trailing edges have emerged for some contaminants, and that contaminant recharge from the basin's vadose zone is still important. The model's estimated time-dependent basin drainage rates combined with dilution from natural recharge successfully predicted plume tritium and nitrate concentrations. This new understanding of source zone influences can help guide science-based remediation, and improve predictions of the natural attenuation timeframes.

Comprehensive understanding of the transport and deposition of engineered nanoparticles (NPs) in subsurface is required to assess their potential negative impact on the environment. We studied the deposition behavior of functionalized quantum dot (QD) NPs (CdTe) in different types of sands (Accusand, ultrapure quartz, and iron-coated sand) at various solution ionic strengths (IS). The observed transport behavior in ultrapure quartz and iron-coated sand was consistent with conventional colloid deposition theories. However, our results from the Accusand column showed that deposition was minimal at the lowest IS (1 mM) and increased significantly as the IS increased. The effluent breakthrough occurred with a delay, followed by a rapid rise to the maximum normalized concentration of unity. Negligible deposition in the column packed with ultrapure quartz sand (100 mM) and Accusand (1 mM) rules out the effect of straining and suggests the importance of surface charge heterogeneity in QD deposition in Accusand at higher IS. Data analyses further show that only a small fraction of sand surface area contributed in QD deposition even at the highest IS (100 mM) tested. The observed delay in breakthrough curves of QDs was attributed to the fast diffusive mass transfer rate of QDs from bulk solution to the sand surface and QD mass transfer on the solid phase. Scanning electron microscopy (SEM) and energy-dispersive X-ray (EDX) analysis were used to examine the morphology and elemental composition of sand grains. It was observed that there were regions on the sand covered with layers of clay particles. EDX spectra collected from these regions revealed that Si and Al were the major elements suggesting that the clay particles were kaolinite. Additional batch experiments using gold NPs and SEM analysis were performed and it was observed that the gold NPs were only deposited on clay particles originally on the Accusand surface. After removing the clays from the sand surface, we observed negligible QD deposition even at 100 mM IS. We proposed that nanoscale charge heterogeneities on clay particles on Accusand surface played a key role in QD deposition. It was shown that the value of solution IS determined the extent to which the local heterogeneities participated in particle deposition.

Bioreduction-based strategies for remediating uranium (U)-contaminated sediments face the challenge of maintaining the reduced status of U for long times. Because groundwater influxes continuously bring in oxidizing terminal electron acceptors (O2, NO3(-)), it is necessary to continue supplying organic carbon (OC) to maintain the reducing environment after U bioreduction is achieved. We tested the influence of OC supply rates on mobility of previously microbial reduced uranium U(IV) in contaminated sediments. We found that high degrees of U mobilization occurred when OC supply rates were high, and when the sediment still contained abundant Fe(III). Although 900 days with low levels of OC supply minimized U mobilization, the sediment redox potential increased with time as did extractable U(VI) fractions. Molecular analyses of total microbial activity demonstrated a positive correlation with OC supply and analyses of Geobacteraceae activity (RT-qPCR of 16S rRNA) indicated continued activity even when the effluent Fe(II) became undetectable. These data support our hypothesis on the mechanisms responsible for remobilization of U under reducing conditions; that microbial respiration caused increased (bi)carbonate concentration and formation of stable uranyl carbonate complexes, thereby shifted U(IV)/U(VI) equilibrium to more reducing potentials. The data also suggested that low OC concentrations could not sustain the reducing condition of the sediment for much longer time. Bioreduced U(IV) is not sustainable in an oxidizing environment for a very long time.

The oxidation status of uranium in sediments is important because the solubility of this toxic and radioactive element is much greater for U(VI) than for U(IV) species. Thus, redox manipulation to promote precipitation of UO2 is receiving interest as a method to remediate U-contaminated sediments. Presence of Fe and Mn oxides in sediments at much higher concentrations than U requires an understanding of their redox status as well. This study was conducted to determine changes in oxidation states of U, Fe, and Mn in U-contaminated sediments from Oak Ridge National Laboratory. Oxidation states of these elements were measured in real-time and nondestructively using X-ray absorption spectroscopy on sediment columns supplied with synthetic groundwater containing organic carbon (OC, 0, 3, 10, 30, and 100 mM OC as lactate) for over 400 days. In sediments supplied with OCor = 30 mM, 80% of the U was reduced to U(IV), with transient reoxidation at about 150 days. Mn(III,IV) oxides were completely reduced to Mn(II) in sediments infused with OCor = 3 mM. However, Fe remained largely unreduced in all sediment columns, showing that Fe(III) can persist as an electron acceptor in reducing sediments over long times. This result in combination with the complete reduction of all other potential electron acceptors supports the hypothesis that the reactive Fe(III) fraction was responsible for reoxidizing U(IV).

Soils contain the largest inventory of organic carbon on the Earth’s surface. Therefore, it is important to understand how soil organic carbon (SOC) is distributed in soils. This study directly measured SOC distributions within soil microaggregates and its associations with major soil elements from three soil groups (Phaeozem, Cambisol, and Ultisol), using scanning transmission X-ray microscopy (STXM) and near-edge X-ray absorption fine structure (NEXAFS) spectroscopy at a spatial resolution of 30 nm. Unlike previous studies, small intact soil microaggregates were examined directly in order to avoid preparatory procedures that might alter C speciation. We found that SOC exists as distinct particles (10s to 100s of nm) and as ubiquitous thin coatings on clay minerals and iron-oxides coatings. The distinct SOC particles have higher fractions of aromatic C than the coatings. NEXAFS spectra of the C coatings within individual microaggregates were relatively similar. In the Phaeozem soil, the pervasive spectral features were those of phenolic and carboxylic C, while in the Cambisol soil the most common spectral feature was the carboxyl peak. The Ultisol soil displayed a diffuse distribution of aromatic, phenolic, and carboxylic C peaks over all surfaces. In general, a wide range of C functional groups coexist within individual microaggregates. In this work we were able to, for the first time, directly quantify the major mineral elemental (Si, Al, Ca, Fe, K, Ti) compositions simultaneously with C distribution and speciation at the nm to μm scale. These direct microscale measurements will help improve understanding on SOC–mineral associations in soil environments.

Waste forms of uranium (U), such as those in the U.S. Department of Energy's Hanford Site, often contain high concentrations of sodium and phosphorus. Low solubility sodium uranyl phosphates such as sodium meta-autunite have the potential to form mobile colloids that can facilitate transport of this radionuclide. To understand the stability behavior of uranyl phosphate colloids, we synthesized sodium meta-autunite colloids and characterized their morphology, chemical composition, structure, dehydration, and surface charge. The stability of these synthetic, plate-shaped colloids was tested with respect to time and pH. We observed the highest aggregation rate to be at pH 3, with the rate decreasing as pH increases, indicating a higher stability of colloid dispersion under neutral and alkaline pH conditions. The synthetic colloids are all negatively charged, and no isoelectric points were found over a pH range of 3–9. The zeta potentials (ζ) of the colloids in the phosphate solution show a strong pH-dependence in the more acidic range over time, but are relatively constant in the neutral and alkaline pH range. Using Derjaguin–Landau–Verwey–Overbeek theory, the stability behavior of these synthetic colloids was interpreted. The high stability of sodium meta-autunite colloids under neutral–alkaline conditions suggests that formation of sodium meta-autunite colloids can enhance the transport of U in U-contaminated sediments.

Nuclear weapons and fuel production have left many soils and sediments contaminated with toxic levels of uranium (U). Although previous short-term experiments on microbially mediated U(VI) reduction have supported the prospect of immobilizing the toxic metal through formation of insoluble U(IV) minerals, our longer-term (17 months) laboratory study showed that microbial reduction of U can be transient, even under sustained reducing conditions. Uranium was reduced during the first 80 days, but later (100-500 days) reoxidized and solubilized, even though a microbial community capable of reducing U(VI) was sustained. Microbial respiration caused increases in (bi)-carbonate concentrations and formation of very stable uranyl carbonate complexes, thereby increasing the thermodynamic favorability of U(IV) oxidation. We propose that kinetic limitations including restricted mass transfer allowed Fe-(III) and possibly Mn(IV) to persist as terminal electron acceptors (TEAs) for U reoxidation. These results show that in-situ U remediation by organic carbon-based reductive precipitation can be problematic in sediments and groundwaters with neutral to alkaline pH, where uranyl carbonates are most stable.