Berkeley Lab

Anoxia stimulates microbially catalyzed metal release from Animas River sediments

Terrain/Satellite view of the study site. Stars indicate sampling locations and the red circle indicates the location of the Gold King Mine.

Published in Environmental Science: Processes & Impacts, this study investigated the fate of heavy metals adsorbed onto riverbed sediments following the August 2015 Gold King Mine Spill in Colorado’s San Juan Mountains. It represents the first biogeochemical study of spill-impacted sediments from the Animas River and revealed mobilization of zinc, arsenic, and molybdenum species accompanying microbe-catalyzed dissolution of metal oxide sorbents.

Concentration changes in aqueous metal cations and anions from the three sediment types across 28-day microcosm incubations. The dashed line in the SO42- panel indicates the time point where exogenous SO42- was added to microcosms to stimulate additional SO42- reduction.


Following the Gold King Mine waste spill, metal contaminants adsorbed onto riverbed sediments along the spill flow path. Differences in sediment mineralogy and adsorbed metals were strongly linked to sampling locations and proximity to the mine. Results suggest that anaerobic microbial metabolisms, stimulated by natural organic carbon pools, will play a significant role in mobilizing adsorbed metal pools following the onset of anoxia in buried riverbed sediments. The site-specific nature of metal release may be linked to different reductive metabolisms, with microbial iron reduction driving dissolution of grain coatings, and alkalinity increases during sulfate reduction offering another mechanism for metal desorption from fluvial sediments. Given the iron and sulfur-rich nature of the Colorado Basin, these complex processes represent a challenge for the tracking of mining-impacted biogeochemistry and associated water quality issues, and emphasize the need for monitoring efforts that account for the dynamic nature of fluvial systems, and their ability to moderate strong spatial and temporal gradients in redox status.

This study provides valuable insight into metal mobility, particularly in mining-impacted watersheds. These results highlight the importance of long-term river water quality monitoring as river sediments undergo sedimentation and burial processes, driving the onset of anoxic conditions which favor metal (re)mobilization.


Saup, C.M., K.H. Williams, L. Rodriguez-Freire, J. Manuel Cerrato, M.D. Johnston, M.J. Wilkins (2017). Anoxia stimulates microbially catalyzed metal release from Animas River sediments. Environmental Science: Processes & Impacts. doi:10.1039/C7EM00036G

Improved modeling of floodplain nutrient and metal cycling using new multi-omics and isotope fractionation information

Water table peaking event mixes oxygen and nitrate into the anoxic Rifle floodplain aquifer. Naturally reduced zones containing sediments higher in organic matter, iron sulfides, and U(IV) rapidly consume DO and nitrate to maintain anoxic conditions, yielding Fe(II) from FeS oxidation, nitrite from denitrification, and U(VI) from nitrite-promoted U(IV) oxidation. Redox cycling is facilitated by coupled geochemistry, heterotrophy, and chemolithoautotrophy.

The study, published in Environmental Science & Technology, is one of the most comprehensive syntheses of processes and datatypes published to date and links BER’s investments in advanced multi-omics techniques and high-performance computing, representing the first critical stage in building a predictive understanding of natural hydrobiogeochemical processes at floodplain and larger scales.


Three-dimensional variably saturated flow and multicomponent biogeochemical reactive transport modeling 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 of the Colorado River in Rifle, Colorado. In this river-aquifer-vadose zone 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.


Yabusaki, S., M. Wilkins, Y. Fang, K. Williams, B. Arora, J. Bargar, H. Beller, N. Bouskill, E. Brodie, J. Christensen, M. Conrad, R. Danczak, E. King, M. Soltanian, N. Spycher, C. Steefel, T. Tokunaga, R. Versteeg, S. Waichler, H. Wainwright (2017).  Water Table Dynamics and Biogeochemical Cycling in a Shallow, Variably-Saturated Floodplain.  Environmental Science and Technology, 51 (6), 3307-3317,  DOI: 10.1021/acs.est.6b04873

Diverse Microbial Metabolism in Aquifer BGC Hot Spot


Summary of key metabolic pathways expressed by a prominent bacterium (Hydrogenophaga b174) in an NRZ biogeochemical hot spot in the Rifle aquifer. Surprisingly, this bacterium actively catalyzed both heterotrophic and chemolithoautotrophic processes and influenced biogeochemical cycling of several elements, including C, N, and S. Unexpectedly, denitrification played an important role in this metabolism.

Organic matter deposits in alluvial aquifers have been shown to result in the formation of NRZs, which can modulate aquifer redox status and influence the speciation and mobility of metals, significantly affecting groundwater geochemistry. This study (Jewell et al., Frontiers in Microbiology) sought to better understand how natural organic matter fuels microbial communities within anoxic biogeochemical hot spots (NRZs) in a shallow alluvial aquifer at the Rifle (CO) site. Overall, the results highlighted the complex nature of organic matter transformation in NRZs and the microbial metabolic pathways that interact to mediate redox status and elemental cycling.


The authors used an anaerobic microcosm experiment in which NRZ sediments served as the sole source of electron donors and microorganisms. Biogeochemical data indicated that the decomposition of native organic matter occurred in different phases, beginning with mineralization of dissolved organic matter (DOM) to CO2 during the first week of incubation, followed by a pulse of acetogenesis that dominated carbon flux after two weeks. The depletion of DOM over time was strongly correlated with increases in expression of many genes associated with heterotrophy (e.g., amino acid, fatty acid, and carbohydrate metabolism) belonging to a Hydrogenophaga strain that accounted for a relatively large percentage (~8%) of the metatranscriptome. This Hydrogenophaga strain also expressed genes indicative of chemolithoautotrophy, including CO2 fixation, H2 oxidation, S-compound oxidation, and denitrification. The pulse of acetogenesis appears to have been collectively catalyzed by a number of different organisms and metabolisms, most prominently pyruvate:ferredoxin oxidoreductase. Unexpected genes were identified among the most highly expressed (>98th percentile) transcripts, including acetone carboxylase and cell-wall-associated hydrolases with unknown substrates. Many of the most highly expressed hydrolases belonged to a Ca. Bathyarchaeota strain and may have been associated with recycling of bacterial biomass.

New Bacteria Groups, and Stunning Diversity, Discovered Underground

Image: Bacterial tree of life (credit: Banfield group)

In a new publication in Nature Communications (Anantharaman et al. 2016), research borne from the Watershed Function SFA provides new clues about the roles of subsurface microbes in globally important cycles. Read More »

Vadose zone respiration contributes significantly to CO2 fluxes


First column: Rifle (Colorado) floodplain vadose zone profile.
Second column: Instrumentation for monitoring pore water and gas profiles down to 3.5 m depth.
Third column: Respiration profiles are sustained by organic carbon carried in infiltration water.

In this SFA study (Tokunaga et al., Vadose Zone Journal), vertical profiles of carbon dioxide concentrations in pore gases, measured from the soil surface down to the water table in a semi-arid floodplain, were shown to have significant contributions from microbial respiration in the deeper vadose zone, well below the rooting depth and above the water table. Approximately 17% of the surface CO2 flux originates from depths between 2.0 and 3.5 meters. Respiration in the vadose zone and low net infiltration limit organic carbon transport to the aquifer. These contributions are not typically accounted for in Earth System Models (ESMs).


Although CO2 fluxes from soils are often assumed to originate within shallow soil horizons (< 1 m depth), relatively little is known about respiration rates at greater depths. The authors compared measured and calculated CO2 fluxes at the Rifle floodplain along the Colorado River, and measured CO2 production rates of floodplain sediments in order to determine the relative importance of vadose zone respiration. Calculations based on measured CO2 gradients and estimated effective diffusion coefficients yielded fluxes that are generally consistent with measurements obtained at the soil surface (326 g C m-2 yr-1). CO2 production from the 2.0 to 3.5 m depth interval was calculated to contribute 17% of the total floodplain respiration, with rates that were larger than some parts of the shallower vadose zone and underlying aquifer. Microbial respiration rates determined from laboratory incubation tests of the sediments support this conclusion. The deeper unsaturated zone typically maintains intermediate water and air saturations, lacks extreme temperatures and salinities, and is annually resupplied with organic carbon from snowmelt-driven recharge and by water table decline. This combination of favorable conditions supports deeper unsaturated zone microbial respiration throughout the year.

Innovative Geophysical Approach for Mapping Biogeochemical Hotspots in Floodplains

The new method (Wainwright et al., Water Resources Research) can be used to map biogeochemical hotspots in a spatially extensive and non-invasive manner, and is particularly useful for developing and parameterizing 3D biogeochemical models. Being able to map such hotspots will enable scientists to quantify the impact of floodplain subsurface processes on biogeochemical cycling through the terrestrial system and associated water quality exports to rivers.

3D estimation of subsurface geological layers

A three-dimensional estimation of subsurface geological layers (green and purple layers) and the naturally reduced zones (red) obtained using the new method at the Rifle Colorado study site. The region shown is the modeling domain, and the estimates obtained were used to parameterize a reactive transport model.


In floodplain environments, naturally reduced zones (NRZ) are considered to be common biogeochemical hotspots that have distinct microbial and geochemical characteristics. Although important for understanding their role in mediating floodplain biogeochemical processes, mapping the subsurface distribution of NRZs over the dimensions of a floodplain or larger is challenging, as conventional wellbore data are typically spatially limited and the distribution of NRZs is heterogeneous. In this study, the authors present an innovative methodology for the probabilistic mapping of NRZs within a three-dimensional (3D) subsurface domain using induced polarization imaging, which is a non-invasive geophysical technique. Measurements consisted of surface geophysical surveys and drilling-recovered sediments, both acquired at the U.S. Department of Energy field site near Rifle, CO. Inversion of surface time-domain induced polarization (TDIP) data yielded 3D images of the complex electrical resistivity, in terms of magnitude and phase, which are associated with mineral precipitation and other lithological properties. To estimate the spatial distribution of NRZs, the team developed a Bayesian hierarchical model to integrate the geophysical and wellbore data. Results showed that the integration of geophysical imaging and wellbore data using this proposed approach was capable of mapping spatially heterogeneous interfaces and NRZ distributions. As described by a follow-on paper, the geophysical estimates were subsequently used to parameterize a reactive transport model, to predict floodplain biogeochemical cycles and river export.

A New View of the Tree of life

Tree of Life

An artistic representation of the tree of life, with the many groups of bacteria on the left, the uncultivable bacteria at upper right (purple), and the Archaea and eukaryotes (green) – which includes humans – at the lower right. (Graphic by Zosia Rostomian, Lawrence Berkeley National Laboratory)


Supported primarily by the Watershed Function SFA, UC Berkeley researchers have reconfigured the tree of life to account for newly discovered microscopic life forms. The new tree, published online in the journal Nature Microbiology, reinforces once again that the life we see around us – plants, animals, humans and other so-called eukaryotes – represent a tiny percentage of the world’s biodiversity. The high-impact paper has generated a great deal of press coverage.

Read more about the research in this press release from UC Berkeley.