Disentangling the complexity of permafrost soil by using high resolution profiling of microbial community composition, key functions and respiration rates.

Thawing permafrost can stimulate microbial activity, leading to faster decomposition of formerly preserved organic matter and CO2 release. Detailed knowledge about the vertical distribution of the responsible microbial community that is changing with increasing soil depth is limited. In this study, we determined the microbial community composition from cores sampled in a high Arctic heath at Svalbard, Norway; spanning from the active layer (AL) into the permafrost layer (PL). A special aim has been on identifying a layer of recently thawed soil, the transition zone (TZ), which might provide new insights into the fate of thawing permafrost. A unique sampling strategy allowed us to observe a diverse and gradually shifting microbial community in the AL, a Bacteroidetes dominated community in the TZ and throughout the PL, a community strongly dominated by a single Actinobacteria family (Intrasporangiaceae). The contrasting abundances of these two taxa caused a community difference of about 60%, just within 3 cm from TZ to PL. We incubated subsamples at about 5°C and measured highest CO2 production rates under aerobic incubations, yet contrasting for five different layers and correlating to the microbial community composition. This high resolution strategy provides new insights on how microbial communities are structured in permafrost and a better understanding of how they respond to thaw.

Wollastonite Carbonation in Water-Bearing Supercritical CO2: Effects of Particle Size.

The performance of geologic CO2 sequestration (GCS) can be affected by CO2 mineralization and changes in the permeability of geologic formations resulting from interactions between water-bearing supercritical CO2 (scCO2) and silicates in reservoir rocks. However, without an understanding of the size effects, the findings in previous studies using nanometer- or micrometer-size particles cannot be applied to the bulk rock in field sites. In this study, we report the effects of particle sizes on the carbonation of wollastonite (CaSiO3) at 60 °C and 100 bar in water-bearing scCO2. After normalization by the surface area, the thickness of the reacted wollastonite layer on the surfaces was independent of particle sizes. After 20 h, the reaction was not controlled by the kinetics of surface reactions but by the diffusion of water-bearing scCO2 across the product layer on wollastonite surfaces. Among the products of reaction, amorphous silica, rather than calcite, covered the wollastonite surface and acted as a diffusion barrier to water-bearing scCO2. The product layer was not highly porous, with a specific surface area 10 times smaller than that of the altered amorphous silica formed at the wollastonite surface in aqueous solution. These findings can help us evaluate the impacts of mineral carbonation in water-bearing scCO2.

Gravity-driven controls on fluid and carbonate precipitation distributions in fractures.

Many challenges related to carbon-dioxide ([Formula: see text]) sequestration in subsurface rock are linked to the injection of fluids through induced or existing fracture networks and how these fluids are altered through geochemical interactions. Here, we demonstrate that fluid mixing and carbonate mineral distributions in fractures are controlled by gravity-driven chemical dynamics. Using optical imaging and numerical simulations, we show that a density contrast between two miscible fluids causes the formation of a low-density fluid runlet that increases in areal extent as the fracture inclination decreases from 90[Formula: see text] (vertical fracture plane) to 30[Formula: see text]. The runlet is sustained over time and the stability of the runlet is controlled by the gravity-driven formation of 3D vortices that arise in a laminar flow regime. When homogeneous precipitation was induced, calcium carbonate covered the entire surface for horizontal fractures (0[Formula: see text]). However, for fracture inclinations greater than 10[Formula: see text], the runlet formation limited the areal extent of the precipitation to less than 15% of the fracture surface. These insights suggest that the ability to sequester [Formula: see text] through mineralization along fractures will depend on the fracture orientation relative to gravity, with horizontal fractures more likely to seal uniformly.

Dewetting of silica surfaces upon reactions with supercritical CO2 and brine: pore-scale studies in micromodels.

Wettability of reservoir minerals and rocks is a critical factor controlling CO(2) mobility, residual trapping, and safe-storage in geologic carbon sequestration, and currently is the factor imparting the greatest uncertainty in predicting capillary behavior in porous media. Very little information on wettability in supercritical CO(2) (scCO(2))-mineral-brine systems is available. We studied pore-scale wettability and wettability alteration in scCO(2)-silica-brine systems using engineered micromodels (transparent pore networks), at 8.5 MPa and 45 °C, over a wide range of NaCl concentrations up to 5.0 M. Dewetting of silica surfaces upon reactions with scCO(2) was observed through water film thinning, water droplet formation, and contact angle increases within single pores. The brine contact angles increased from initial values near 0° up to 80° with larger increases under higher ionic strength conditions. Given the abundance of silica surfaces in reservoirs and caprocks, these results indicate that CO(2) induced dewetting may have important consequences on CO(2) sequestration including reducing capillary entry pressure, and altering quantities of CO(2) residual trapping, relative permeability, and caprock integrity.

Cross-Scale Molecular Analysis of Chemical Heterogeneity in Shale Rocks.

The organic and mineralogical heterogeneity in shale at micrometer and nanometer spatial scales contributes to the quality of gas reserves, gas flow mechanisms and gas production. Here, we demonstrate two molecular imaging approaches based on infrared spectroscopy to obtain mineral and kerogen information at these mesoscale spatial resolutions in large-sized shale rock samples. The first method is a modified microscopic attenuated total reflectance measurement that utilizes a large germanium hemisphere combined with a focal plane array detector to rapidly capture chemical images of shale rock surfaces spanning hundreds of micrometers with micrometer spatial resolution. The second method, synchrotron infrared nano-spectroscopy, utilizes a metallic atomic force microscope tip to obtain chemical images of micrometer dimensions but with nanometer spatial resolution. This chemically "deconvoluted" imaging at the nano-pore scale is then used to build a machine learning model to generate a molecular distribution map across scales with a spatial span of 1000 times, which enables high-throughput geochemical characterization in greater details across the nano-pore and micro-grain scales and allows us to identify co-localization of mineral phases with chemically distinct organics and even with gas phase sorbents. This characterization is fundamental to understand mineral and organic compositions affecting the behavior of shales.

Landscape topography structures the soil microbiome in arctic polygonal tundra.

In the Arctic, environmental factors governing microbial degradation of soil carbon (C) in active layer and permafrost are poorly understood. Here we determined the functional potential of soil microbiomes horizontally and vertically across a cryoperturbed polygonal landscape in Alaska. With comparative metagenomics, genome binning of novel microbes, and gas flux measurements we show that microbial greenhouse gas (GHG) production is strongly correlated to landscape topography. Active layer and permafrost harbor contrasting microbiomes, with increasing amounts of Actinobacteria correlating with decreasing soil C in permafrost. While microbial functions such as fermentation and methanogenesis were dominant in wetter polygons, in drier polygons genes for C mineralization and CH4 oxidation were abundant. The active layer microbiome was poised to assimilate N and not to release N2O, reflecting low N2O flux measurements. These results provide mechanistic links of microbial metabolism to GHG fluxes that are needed for the refinement of model predictions.

Composite thermal conductivity in a large heterogeneous porous methane hydrate sample.

By employing inverse modeling to analyze the laboratory data, we determined the composite thermal conductivity (k(theta), W/m/K) of a porous methane hydrate sample ranged between 0.25 and 0.58 W/m/K as a function of density. The calculated composite thermal diffusivities of porous hydrate sample ranged between 2.59 x 10(-7) m(2)/s and 3.71 x 10(-7) m(2)/s. The laboratory study involved a large heterogeneous sample (composed of hydrate, water, and methane gas). The measurements were conducted isobarically at 4.98 MPa over a temperature range of 277.3-279.1 K. Pressure and temperature were monitored at multiple locations in the sample. X-ray computed tomography (CT) was used to visualize and quantify the density changes that occurred during hydrate formation from granular ice. CT images showed that methane hydrate formed from granular ice was heterogeneous and provided an estimate of the sample density variation in the radial direction. This facilitated quantifying the density effect on composite thermal conductivity. This study showed that the sample heterogeneity should be considered in thermal conductivity measurements of hydrate systems. Mixing models (i.e., arithmetic, harmonic, geometric mean, and square root models) were compared to the estimated composite thermal conductivity determined by inverse modeling. The results of the arithmetic mean model showed the best agreement with the estimated composite thermal conductivity.

Non-aqueous phase liquid spreading during soil vapor extraction.

Many non-aqueous phase liquids (NAPLs) are expected to spread at the air-water interface, particularly under non-equilibrium conditions. In the vadose zone, this spreading should increase the surface area for mass transfer and the efficiency of volatile NAPL recovery by soil vapor extraction (SVE). Observations of spreading on water wet surfaces led to a conceptual model of oil spreading vertically above a NAPL pool in the vadose zone. Analysis of this model predicts that spreading can enhance the SVE contaminant recovery compared to conditions where the liquid does not spread. Experiments were conducted with spreading volatile oils hexane and heptane in wet porous media and capillary tubes, where spreading was observed at the scale of centimeters. Within porous medium columns up to a meter in height containing stagnant gas, spreading was less than ten centimeters and did not contribute significantly to hexane volatilization. Water film thinning and oil film pinning may have prevented significant oil film spreading, and thus did not enhance SVE at the scale of a meter. The experiments performed indicate that volatile oil spreading at the field scale is unlikely to contribute significantly to the efficiency of SVE.

Experimental and numerical simulation of dissolution and precipitation: implications for fracture sealing at Yucca Mountain, Nevada.

Plugging of flow paths caused by mineral precipitation in fractures above the potential repository at Yucca Mountain, Nevada could reduce the probability of water seeping into the repository. As part of an ongoing effort to evaluate thermal-hydrological-chemical (THC) effects on flow in fractured media, we performed a laboratory experiment and numerical simulations to investigate mineral dissolution and precipitation under anticipated temperature and pressure conditions in the repository. To replicate mineral dissolution by vapor condensate in fractured tuff, water was flowed through crushed Yucca Mountain tuff at 94 degrees C. The resulting steady-state fluid composition had a total dissolved solids content of about 140 mg/l; silica was the dominant dissolved constituent. A portion of the steady-state mineralized water was flowed into a vertically oriented planar fracture in a block of welded Topopah Spring Tuff that was maintained at 80 degrees C at the top and 130 degrees C at the bottom. The fracture began to seal with amorphous silica within 5 days.A 1-D plug-flow numerical model was used to simulate mineral dissolution, and a similar model was developed to simulate the flow of mineralized water through a planar fracture, where boiling conditions led to mineral precipitation. Predicted concentrations of the major dissolved constituents for the tuff dissolution were within a factor of 2 of the measured average steady-state compositions. The mineral precipitation simulations predicted the precipitation of amorphous silica at the base of the boiling front, leading to a greater than 50-fold decrease in fracture permeability in 5 days, consistent with the laboratory experiment.These results help validate the use of a numerical model to simulate THC processes at Yucca Mountain. The experiment and simulations indicated that boiling and concomitant precipitation of amorphous silica could cause significant reductions in fracture porosity and permeability on a local scale. However, differences in fluid flow rates and thermal gradients between the experimental setup and anticipated conditions at Yucca Mountain need to be factored into scaling the results of the dissolution/precipitation experiments and associated simulations to THC models for the potential Yucca Mountain repository.

Thermal dissociation behavior and dissociation enthalpies of methane-carbon dioxide mixed hydrates.

Replacement of methane with carbon dioxide in hydrate has been proposed as a strategy for geologic sequestration of carbon dioxide (CO(2)) and/or production of methane (CH(4)) from natural hydrate deposits. This replacement strategy requires a better understanding of the thermodynamic characteristics of binary mixtures of CH(4) and CO(2) hydrate (CH(4)-CO(2) mixed hydrates), as well as thermophysical property changes during gas exchange. This study explores the thermal dissociation behavior and dissociation enthalpies of CH(4)-CO(2) mixed hydrates. We prepared CH(4)-CO(2) mixed hydrate samples from two different, well-defined gas mixtures. During thermal dissociation of a CH(4)-CO(2) mixed hydrate sample, gas samples from the head space were periodically collected and analyzed using gas chromatography. The changes in CH(4)-CO(2) compositions in both the vapor phase and hydrate phase during dissociation were estimated based on the gas chromatography measurements. It was found that the CO(2) concentration in the vapor phase became richer during dissociation because the initial hydrate composition contained relatively more CO(2) than the vapor phase. The composition change in the vapor phase during hydrate dissociation affected the dissociation pressure and temperature; the richer CO(2) in the vapor phase led to a lower dissociation pressure. Furthermore, the increase in CO(2) concentration in the vapor phase enriched the hydrate in CO(2). The dissociation enthalpy of the CH(4)-CO(2) mixed hydrate was computed by fitting the Clausius-Clapeyron equation to the pressure-temperature (PT) trace of a dissociation test. It was observed that the dissociation enthalpy of the CH(4)-CO(2) mixed hydrate lays between the limiting values of pure CH(4) hydrate and CO(2) hydrate, increasing with the CO(2) fraction in the hydrate phase.