Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Unsaturated soil”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3

A Rapid Decision Support Tool for Estimating Impacts of a Vadose Zone Volatile Organic Compound Source on Groundwater and Soil Gas

Diminishing rates of subsurface volatile contaminant removal by soil vapor extraction (SVE) oftentimes warrants an in‐depth performance assessment to guide remedy decision‐making processes. Such a performance assessment must include quantitative approaches to better understand the impact of remaining vadose zone contamination on soil gas and groundwater concentrations. The spreadsheet‐based Soil Vapor Extraction Endstate Tool (SVEET) software functionality has recently been expanded to facilitate quantitative performance assessments. The updated version, referred to as SVEET2, includes expansion of the input parameter ranges for describing a site (site geometry, source characteristics, etc.), an expanded list of contaminants, and incorporation of elements of the Vapor Intrusion Estimation Tool for Unsaturated‐zone Sources software to provide soil gas concentration estimates for use in vapor intrusion evaluation. As part of the update, SVEET2 was used to estimate the impact of a tetrachloroethene (PCE) vadose zone source on groundwater concentrations, comparing SVEET2 results to field‐observed values at an undisclosed site where SVE was recently terminated. PCE concentrations from three separate monitoring wells were estimated by SVEET2 to be within the range of 6.0–6.7 μg/L, as compared to actual field concentrations that ranged from 3 to 11 μg/L PCE. These data demonstrate that SVEET2 can rapidly provide representative quantitative estimates of impacts from a vadose zone contaminant source at field sites. Finally, in the context of the SVE performance assessment, such quantitative estimates provide a basis to support remedial and/or regulatory decisions regarding the continued need for vadose zone volatile organic compound remediation or technical justification for SVE termination, which can significantly reduce the cost to complete for a site.

54 ENVIRONMENTAL SCIENCES↗

Large-scale paleo water-table rise in a deep desert aquifer recorded by dissolved noble gases

Illuminating past hydrological and climatological conditions in arid regions may provide insights into future trends in groundwater availability. The main objective of this study is to explore paleorecharge processes in the deep regional Nubian Sandstone Aquifer (NSA), which stretches below the hyperarid deserts of the Sinai Peninsula (Egypt) and the Negev (Israel), using dissolved noble gas data. Extremely high amounts of dissolved excess air (Ne concentrations up to 4 × solubility equilibrium) were observed in the ancient ( 81 Kr-depleted) groundwater that was recharged during past pluvial epochs in the basin. Here, the observed unique excess air signal and the clear spatiotemporal structure of the noble gas compositions in the aquifer are hypothesized to reflect two major characteristics of the groundwater system: (a) large-scale, long-term rises in the water table that facilitate the entrapment and dissolution of air bubbles under increased hydrostatic pressure, and (b) the geological settings in the southern recharge area in Sinai, including an almost horizontal position of the layers and intercalation of low-permeability rock formations within the aquifer, which inefficiently reject air bubbles from groundwater, thereby allowing for substantial entrapment of bubbles that dissolve as the water table rises. Enhanced dissolution of entrapped air was also observed in other paleo-groundwater-containing regional sandstone aquifers across North Africa, which, together with the recent findings from the NSA, suggest that largescale water table fluctuations have likely extensively occurred during past pluvial periods over these (and possibly other) arid regions. A preliminary assessment of noble gas recharge temperatures (NGTs) indicates an apparent decoupling of surface and water table temperatures in the case of deep aquifers with a thick (hundreds of m) unsaturated zone. This observed decoupling calls for a re-evaluation of previously obtained NGT records and a need for future work to consider the modification of NGTs due to soil air fractionation and geothermal heating of the underlying deep unsaturated zones.

54 ENVIRONMENTAL SCIENCES↗

Decoupling of redox processes from soil saturation in Arctic tundra

Permafrost thaw in warming Arctic landscapes alters hydrology and saturation-driven biogeochemical processes. Models assume that aerobic respiration occurs in drained soils while saturated soils support methanogenesis; however, saturated soils maintain redox gradients that host a range of anaerobic metabolisms. We evaluated how redox potential and redox-active solutes vary with soil moisture in the active layer of permafrost-affected acidic and non-acidic tundra hillslopes. Oxidizing conditions persisted in highly permeable organic horizons of both unsaturated tussock tundra and saturated wet sedge meadows. Redox potential decreased with depth in all soils as increasing soil bulk density restricted groundwater flow and oxygen diffusion. High concentrations of dissolved iron, phosphate, and organic carbon coincided with redox boundaries below the soil surface in acidic tundra, indicating active iron redox cycling and potential release of adsorbed phosphate during iron (oxyhydr)oxide dissolution. In non-acidic tundra, weatherable minerals affected nutrient dynamics more than redox-driven iron cycling, especially in low-lying, saturated areas where thaw reached mineral soils. The role of thaw depth and the ability of saturated soils to maintain oxidizing conditions in organic surface layers highlight the importance of soil physical properties and hydrology in predicting biogeochemical processes and greenhouse gas emissions.

54 ENVIRONMENTAL SCIENCES↗

Electrical Resistivity Tomography Monitoring of In Situ Soil Flushing at the Hanford 100-K East Area: 100KE Soil Flushing Monitoring

Former operations in the 100-KE Area at the Hanford Site resulted in vadose zone hexavalent chromium that provides a source of groundwater contamination. Previous efforts to remediate vadose zone chromium involved excavation and offsite treatment of contaminated soils. Although much of the chromium was removed, contamination still exists in native soil between the water table and the bottom of the pit, which has since been backfilled. In-situ soil flushing was tested at the 100-KE Area to accelerate the removal of the remaining hexavalent chromium in the vadose zone. Soil flushing works by applying clean flush water at the surface, which mobilizes and transports chromium as it migrates downward to the water table. Once the chromium reaches the water table, it is hydraulically contained, extracted, and treated through pump and treat operations. The efficacy of soil flushing is directly related to the volume of clean water that infiltrates through contaminated soils. Therefore, performance monitoring requires and understanding of flush water migration through the vadose zone. It is challenging to comprehensively assess soil flushing performance via borehole access due to the limited volume of investigation afforded by a given wellbore. As an alternative, 3D time-lapse electrical resistivity tomography (ERT) was tested as a method of monitoring the distribution of flush water over time within the vadose zone. ERT is a method of imaging the bulk electrical conductivity of the subsurface, which is highly dependent on water saturation levels in unconsolidated and unsaturated sediments. Thus, the timing and location of changes in bulk conductivity caused by soil flushing can be used to infer flush water migration pathways. In this report we present results of 3D time-lapse ERT imaging during two separate soil flushing campaigns conducted during the spring and summer months of 2022 and 2023. Results show generally that: 1) Pit backfill materials appear to nominally have larger permeability than native soils. Consequently, the boundary between pit backfill and native soil had a significant impact on flush water migration, causing some flush water to migrate along the pit boundary to the bottom of the pit. 2) Non-uniform flows, likely caused by variations in hydrogeologic properties, developed in the pit backfill materials, resulting in uneven flush water distribution on the southern margin of the soil flushing zone. 3) Redistribution of water at the interface between backfill materials and Hanford formation materials likely facilitated elevated and uniform distribution of flush water within native soils beneath the deeper parts of the pit boundary beneath the center of the flush zone, which presumably overlie soils with elevated chromium contamination. These areas appear to have been infiltrated by higher volumes of flush water than the northern and southern margins of the flush zone. 4) In comparison to 2022, high flush water application rates significantly improved flush water distribution throughout the target flushing zone. 5) Imaging resolution was limited to a depth of approximate 20 meters below the land surface, due primarily to limitations on the lateral extent of the surface ERT array. The ~10 m region of the vadose zone between approximately 130 meters elevation and the water table at approximately 120 meters elevation was unresolved.

100 Area↗

Biogeochemical controls on iron speciation and cycling across upland to shoreline gradients in freshwater and estuarine coastal soils (Lake Erie and Chesapeake Bay, United States)

Coastal environments are dynamic interfaces that mediate carbon and nutrient exchanges between terrestrial landscapes and open waters, and understanding the biogeochemical factors controlling these exchanges, particularly iron (Fe) redox transformations, is crucial for predicting coastal ecosystem functions. Here, we investigated the mechanisms controlling Fe speciation changes across upland-to-shoreline gradients in freshwater and estuarine soils using Fe K-edge X-ray absorption spectroscopy, solid and porewater composition analysis, and 16S rRNA sequencing analysis. We show that Fe transformations depend primarily on inundation patterns. In unsaturated uplands, Fe occurs as Fe(III) oxyhydroxides, mainly goethite (9–35 %), Fe(II,III)-phyllosilicates (39–89 %), and Fe(III)-organic species (0–61 %). Soils influenced by estuarine waters exhibit porewater sulfide concentrations reaching up to 221 μM, Fe- and S-cycling bacteria, and up to 81 % pyrite (FeS 2 ), indicating that sulfur-driven redox dynamics control Fe transformations. In lacustrine wetlands, Fe(III) reduction is indicated by porewater Fe(II) concentrations increasing to 1.0–2.1 mM, and ~10–15 % of Fe as Fe(II,III)-(hydr)oxides (green rust), vivianite (Fe 3 (PO 4 ) 2 ·8H 2 O), and/or adsorbed Fe(II) species. EXAFS data also indicate reduction of structural Fe(III) to Fe(II) in phyllosilicates. The presence of Fe- and S-cycling bacteria, as well as sulfide (0–10 μM), suggests that Fe-cycling is microbially driven and potentially coupled with cryptic S-cycling. Fe(II) oxidation was indicated above/near the water table by the presence of Fe(III) oxyhydroxides (ferrihydrite, lepidocrocite). Furthermore, negligible Fe(III) or sulfate reduction was observed at some water-saturated sites located at the upland-wetland transition, likely due to oxic (sub-)surface water inputs. Overall, our results highlight the importance of considering both Fe-speciation and hydro-biogeochemical dynamics when predicting Fe-cycling at coastal interfaces.

54 ENVIRONMENTAL SCIENCES↗

Benzene homologues contaminants in a former herbicide factory site: distribution, attenuation, risk, and remediation implication

Benzene homologues often used as organic raw materials or as detergents in chemical industry are prone to accidental release into the environment which can cause serious long-term soil pollutions. In a large former herbicide factory site, we investigated 43 locations for benzene homologues contaminations in soil, soil gas, and groundwater; and studied the hydrogeological conditions. An inverse distance weighted interpolation method was employed to determine the pollutants three-dimensional spatial distribution in the soils. Results showed that benzene homologues residues were mainly originated from the herbicide production workshop; and that the pollution had horizontally expanded at the deeper soil layer. Contaminants had already migrated 15m downward from ground surface. Contaminant phase distribution study showed that NAPL was the primary phase (>99%) for the pollutants accumulated in the unsaturated zone; while it had not migrated to groundwater. The primary mechanism for contaminant transport and attenuation included dissolution of “occluded” NAPL into pore water and pollutant volatilization into soil pore space. Risk assessment revealed that the pollutants brought unacceptable high carcinogenic and non-carcinogenic risks to public health. In order to convert this former chemical processing factory site into a residential area, a remediation to the polluted production workshop sites is urgently required.

Yang, Shuo↗

Organic amendments change soil organic C structure and microbial community but not total organic matter on sub-decadal scales

Organic C has many benefits for soil, but it is depleted by tillage and crop harvest, and especially so for biofuel crops. Accordingly, strategies such as partially retaining stover or planting a cover crop can help ameliorate the negative effect of C removal. We used a long-term field experiment to study the impacts of stover retention and planting cover crop on soil organic matter (SOM), its extractable components and the soil microbial community. SOM chemical composition characterization was determined by electrospray ionization (ESI) coupled with Fourier transform ion cyclotron resonance mass spectrometry (FT-ICR-MS) in sequential water, methanol (MeOH), and chloroform (CHCl3) extracts. The characteristics of the soil bacterial community were measured by phospholipid fatty acid (PLFA), real-time quantitative PCR, and 16S rRNA gene sequence. The variations in total SOM content, total microbial biomass, and bacterial population were slight among treatments, but SOM chemical compounds, arbuscular mycorrhizal fungi (AMF) biomass, and bacterial structure changed significantly, and especially so in the coupled application of stover retention and cover crop. Specifically, stover retention enriched more lignin-like compounds in soil, whereas cover crop enriched more condensed hydrocarbons, and had more compounds with an aromaticity index (AI) >0.5. The bacterial community was not altered by the cover crop, but the corn stover retention increased the relative abundances of Myxococcales (Deltaproteobacteria) and decreased that of Actinobacteria. Redundancy analysis (RDA) further unveiled that the bacterial community in the stover treatments had significantly positive association with CHCl3-extracted chemical classes, i.e. unsaturated hydrocarbon and lipid with the coupled application (stover and cover crop), and lignin and protein with the single corn stover treatment. Taken together, our study shows how different C addition practices influence the molecular composition of SOM and the structure of soil microbial communities.

Xu, Jiangbin↗

Changes in water vapor adsorption and water film thickness in clayey materials as a function of relative humidity

Abstract Soil water, with adsorbed water being an important component, plays an important role in fluid flow and chemical movement in the unsaturated zone. The adsorbed water consists of water film and interlayer water, and its content is related to the relative humidity (RH). This study theoretically and experimentally focuses on the variation in adsorbed water content in clayey materials as the RH changes. Based on a slit‐pore model, three types of water (water film, interlayer water, and capillary water) could be present in the water vapor adsorption process. The variation of water film is obtained from the analysis of interfacial forces, and then the total water content can be quantitatively subdivided into these three types of water for different RH values. The slit‐pore adsorption model was used with experimental measurements (water vapor adsorption and N 2 physisorption) to quantify the amounts and interrelationships of these three different types of water to six clayey materials, namely three clay minerals (kaolinite, montmorillonite, and illite–smectite mixed layer) and three clay‐rich sediments from the Jianghan Plain in China. The main results are that (a) the volumes of water vapor adsorption are much greater than those of N 2 adsorption except for kaolinite; (b) interlayer water largely dominates the growth of total water content in montmorillonite and results in a concentration of 0.090 ml g −1 at 95% RH; and (c) the minimum thickness of water film is calculated to be 0.23 nm, and the maximum value is one‐third of the pore width by considering the interfacial forces.

54 ENVIRONMENTAL SCIENCES↗

AmeriFlux FLUXNET-1F US-ADR Amargosa Desert Research Site (ADRS)

This is the AmeriFlux Management Project (AMP) created FLUXNET-1F version of the carbon flux data for the site US-ADR Amargosa Desert Research Site (ADRS). This is the FLUXNET version of the carbon flux data for the site US-ADR Amargosa Desert Research Site (ADRS) produced by applying the standard ONEFlux (1F) software. Site Description - This tower is located at the Amargosa Desert Research Site (ADRS). The U.S. Geological Survey (USGS) began studies of unsaturated zone hydrology at ADRS in 1976. Over the years, USGS investigations at ADRS have provided long-term "benchmark" information about the hydraulic characteristics and soil-water movement for both natural-site conditions and simulated waste-site conditions in an arid environment. The ADRS is located in a creosote-bush community adjacent to disposal trenches for low-level radioactive waste.

Moreo, Michael [U.S. Geological Survey]↗

Field Test Bed for Vadose-Zone Monitoring Approaches - 20404

In the Central Plateau at the U.S. Department of Energy Hanford Site, a large inventory of contaminants resides in unsaturated sediments within the approximately 100-meter-thick vadose zone, posing a potential continuing risk to groundwater. Vadose zone remedies used to address these contaminants will require performance monitoring to provide feedback during implementation and for long-term verification that remedial action objectives have been met. Passive approaches may also need long-term monitoring to demonstrate that the flux of contaminants from the vadose zone to the groundwater are below thresholds established to meet groundwater protection goals. Collection of physical (e.g., groundwater or sediment) samples is a common method for identifying contaminant concentration distributions, but this approach is limited by the number of locations and the frequency with which data can be collected. In situ vadose zone measurements have evolved over the past few years to include key measurements of water content, soil water pressure, temperature, and chemical concentration. However, the current generation of sensors is designed for relatively short-term use in near-surface soils or sediments. Geophysical methods have been evolving but are also limited in that they have not been designed for the specific long-term vadose zone monitoring needs at the Hanford Site. Overall, monitoring under unsaturated conditions can be difficult due to the need to install and maintain instrumentation over a large area and depth and the need to identify preferential flow pathways due to geologic and chemical heterogeneities over long periods. Thus, a vadose zone monitoring test bed was initiated to address these challenges and identify cost-effective approaches for implementation and postclosure monitoring of the deep vadose zone. The monitoring test bed is expected to provide valuable field-scale information for the design of vadose zone monitoring systems. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Collaborative Research: Natural Organic Matter and Microbial Controls on Mobilization/Immobilization of I and Pu in Soils and Waters Affected by Radionuclide Releases in USA and Japan

In this project, the relationship between natural organic matter (NOM) and two radioactive elements that are relevant to nuclear waste disposal were studied: Plutonium (Pu) and Iodine. The human and environmental risks associated with Pu stem mainly from the very long half-lives of several of its isotopes ( 238 Pu, 88 yr; 239 Pu, 24,100 yr; 240 Pu, 6560 yr) and its radiotoxicity. Understanding Pu biogeochemical behavior in both near-field (>10-11M) and far-field scales (<10-11M) is imperative to the development of approaches for reprocessing Pu, remediation of Pu contamination and accurate assessment of risks posed by disposal practices for Pu-bearing wastes. The environmental mobility of Pu can be affected by redox potential, pH, adsorption, precipitation, complexation, colloid formation, and microbial activity, of which the first characteristic has the most profound influence. Numerous studies have shown high affinity of Pu towards NOM, as well as to mineral phases. NOM is ubiquitous in the environment, e.g., both fulvic and humic acids are able to reduce Pu(V,VI) to Pu(IV) and the redox potential of NOM is positively related to the abundance of phenolic/acidic OH groups. NOM can either facilitate or limit actinide migration, depending on specific biogeochemical conditions including pH, mineral and organic matter characteristics, etc. The other radionuclide of interest is radioiodine ( 129 I). 129I is a major by-product of nuclear fission and of serious concern to the Department of Energy (DOE) as it is among the top risk drivers at existing and potential radiowaste-contaminated sites. The risk of 129 I stems largely from its high bioconcentration factor (90% of the body’s iodine is accumulated in the thyroid), a high inventory at source terms, a very long-half life (16M years), and rapid mobility in the subsurface environment. As a consequence, 129 I has the lowest drinking water standard (1 pCi/L) among all radionuclides in the Federal Register. With a novel and sensitive gas chromatography-mass spectrometry (GC-MS) method developed in our lab, it is possible to quickly and simultaneously determine the distribution of 129 I and stable 127 I forms in environments, as low as 2 pCi/L for 129 I. This method was subsequently validated using accelerator mass spectroscopy, AMS. IO 3 - and organo-I were determined as major species in the groundwater of SRS and the Hanford Site, contrary to thermodynamic predictions that I- should be the dominant species at these sites. Mobility of 129 I was also demonstrated to depend greatly on the I species and its concentration, sediment pH, and redox state, with times to achieve equilibrium taking up to 12 weeks. Along the groundwater pathway in the F-Area of SRS, 129 I- supplied from the seepage basins was transformed to 129 IO 3 - and organo- 129 I with increasing iodine sediment sorption, causing the lower total 127 I and 129 I concentrations along the gradient transect of the waste plume. By contrast, groundwater 129 I concentrations in the wetlands (as high as 1617.3 pCi/L) were greatly elevated with respect to the source term (159.3 pCi/L). While the NOM promoted the uptake of 129 I to the wetland sediment, it also promoted the formation of soluble organic fraction. A small fraction of NOM that is bound to iodine can behave as a mobile organo-I source. Iodide was enzymatically incorporated into NOM, whereas both iodide and iodate were abiotically bound to NOM, under certain conditions. Iodate removal from the mobile aqueous phase can also occur through incorporation into carbonate (e.g., at the Hanford Site, USA). Thus immobilization and re-mobilization of iodine species were influenced by pH, Eh and the presence of NOM and metal oxides, which adds to the complexity of site remediation action. A ground-breaking result was to elucidate the products (i.e. organo-iodine moieties formed via enzymatic and non-enzymatic processes) at the molecular level by nuclear magnetic resonance (NMR) and electrospray ionization Fourier transform ion cyclotron resonance mass spectrometry (ESI-FTICRMS). We found that iodine-NOM interactions may be influenced by NOM hydrophobic aliphatic moieties. From the perspective of ESI-FTICRMS, organo-iodine formulas were ascribed to the groups of unsaturated hydrocarbons, lignins and proteins. Iodate is likely abiotically reduced to reactive iodine species by lignin- and tannin-like compounds or carboxylic-rich alicyclic molecules (CRAM). We also investigated microbial mechanisms in iodine incorporation into NOM. We established that soil bacteria isolated from F-Area of SRS did not accumulate significant amounts of I- (0.2-2%). Intracellular uptake of I- decreases with increasing pH when pH ranged from 4 to 6. In contrast, 44 out of 84 strains isolated from the F-Area of SRS can transform I- to IO 3 - and organo-iodine. In some cases, oxidation was facilitated in the presence of H 2 O 2 . Microbes can also excrete organic acids that enhance I- oxidation by lowering the ambient pH and reacting with H 2 O 2 to form peroxy carboxylic acids. At lower pH values (≤5), H 2 O 2 hydrolysis was the driving force for iodide-oxidation; whereas, at pH ≥ 6, spontaneous decomposition of peroxy carboxylic acids, originating from H 2 O 2 and organic acids were the primary cause of iodide oxidation. Lastly, it was determined that microbial processes involved in Mn (II) are capable of directly oxidizing I- via enzymatic catalysis (i.e., multicopper oxidases), or indirectly through the formation of reactive oxygen species (ROS) and/or biogenic manganese oxides. ROS-mediated oxidation of I- was found to predominate at pH >5, whereas the enzymatic and Mn oxide pathways were more active at pH < 5. Together, this project has resulted in 9 publications in high-impact journals, and the training of 1 Ph.D and 4 undergraduate students.

54 ENVIRONMENTAL SCIENCES↗

Laboratory Evaluation to Increase Effectiveness of Field-Scale Soil Flushing in the Hanford 100 Areas

This laboratory study was initiated to develop an improved technical understanding of the key controlling geochemical and physical processes of mobilizing residual Cr(VI) in Hanford sediments to provide increased effectiveness in soil flushing activities planned for Hanford 100 Areas. This was accomplished by (1) quantifying Cr mass, release rate, and surface phase changes as Cr is leached from different Hanford sediments; (2) maximizing Cr leaching during soil flushing by evaluating different leach solutions and leach strategies; and (3) minimizing leaching of remaining residual Cr after soil flushing ends. To quantify geochemical controls on Cr leaching from sediments, sequential liquid extractions were used to identify aqueous, adsorbed, and solid Cr surface phases (precipitates or Cr incorporated into other phases) before and after water-saturated and unsaturated leach experiments with pH 8 artificial groundwater and with other amendment solutions. Changes in the Cr release rate from sediments were also correlated to changes in Cr surface phases. During unsaturated infiltration, water advection occurs primarily through larger pores, with less movement through smaller grain layers leading to a slower release of Cr trapped in smaller pores next to smaller grains. To quantify these physical controls on Cr leaching from sediments, 10 ft high 1-D infiltration experiments were conducted at differing leach solution application rates, and Cr leaching during application and subsequent residual water flow was quantified. Qualitative 2-D experiments were also conducted to evaluate surfactant addition to increase flow in low-permeability zones. Cr-contaminated Hanford sediments from the 100-H, 100-K, 100-D, and 200 East Areas exhibited a wide variety of leach behavior from rapid release of a small amount of Cr to fast (i.e., minutes) and slow (i.e., hours or more) release of high Cr. Sediments with low labile Cr (< 0.2 µg/g) exhibited fast Cr release but less overall release, and most (>95%) leached from the sediments within a few pore volumes in water-saturated columns. With nearly no sorption of chromate, aqueous and sorbed Cr (as chromate) were released quickly from sediment within the first few pore volumes, resulting in high leach concentrations, fast initial Cr release rates, and changes in pre- and post-leach extractions. In contrast, sediments with high labile Cr (0.2 to 20 µg/g) exhibited a combination of fast and slow Cr release due to multiple Cr surface phases contributing to Cr release from the sediment at different rates with greater overall release. This resulted in elevated Cr effluent concentrations even after dozens of pore volumes were leached. In addition, the Cr release rate decreased with increasing time during leaching. A decrease in labile Cr (aqueous Cr, adsorbed Cr, pH 5 acetate dissolved precipitates) in post-leach extractions compared to pre-leach extractions showed initial aqueous and adsorbed Cr release, followed by dissolution of high-solubility precipitates such as CaCrO 4 (if present), then slow dissolution of calcite that may contain some chromate and possibly BaCrO 4 . The labile Cr in all sediments correlated well with the Cr release rate (r 2 = 0.81). Labile Cr also correlated well with the leached mass in water-saturated columns (r 2 = 0.98), and therefore can be used for prediction of potential removal with flushing applications.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Vadose Zone Soil Flushing for Chromium Remediation: A Laboratory Investigation to Support Field‐scale Application

Cr(VI) flushing from the vadose zone to the groundwater (with subsequent Cr(VI) removal in groundwater by pump-and-treat system) is a promising remedial technique that has recently been used at field scale. This laboratory study was conducted to provide the technical basis to design a field soil flushing strategy. The objectives were to (1) quantify the relationship between sediment Cr(VI) and Cr(III) mass and release rates and subsequent Cr(VI) leaching; (2) investigate different methodologies to maximize Cr(VI) leaching, and (3) investigate methods to minimize leaching of remaining residual Cr. Characterization of Cr-contaminated sediments (Hanford Site, WA) exhibited Cr(VI) showed that leach rates that were correlated to different Cr surface phases. Sediments with low leachable Cr(VI) (<2 μg/g) leached Cr rapidly, so slow infiltration of water in a single pulse was sufficient to leach most Cr. In contrast, sediments with high Cr (2 to 200 μg/g) released some Cr(VI) quickly but 10 to 50% Cr(VI) slowly (tens to hundreds of hours). Efficient unsaturated leaching of these sediments required a different infiltration strategy that includes: multiple slow leach pulses with time between flushing cycles; the use of a surfactant to increase Cr leaching from low-permeability zones, and the use of a reductant (Na-dithionite or Ca-polysulfide) in the final leach water was highly effective at decreasing residual Cr leaching. This study clearly demonstrated that the methodology of basing laboratory Cr flushing on parameters such as Cr release mass and rates could be used to improve the efficiency of soil flushing at field scale.

Accelerated Cr Flushing↗

Multiscale simulations for multi-continuum Richards equations

In this paper, we study a multiscale method for simulating a dual-continuum unsaturated flow problem within complex heterogeneous fractured porous media. Mathematically, each of the dual continua is modeled by a multiscale Richards equation (for pressure head), and these equations are coupled to one another by transfer terms. On its own, Richards equation is already a nonlinear partial differential equation, and it is exceedingly difficult to solve numerically due to the extra nonlinear dependencies involving the soil water. To deal with multiple scales, our strategy is that starting from a microscopic scale, we upscale the coupled system of dual-continuum Richards equations via homogenization by the two-scale asymptotic expansion, to obtain a homogenized system, at an intermediate scale (level). Based on a hierarchical approach, the homogenization’s effective coefficients are computed through solving the arising cell problems. Furthermore, to tackle the nonlinearity, after time discretization, we use Picard iteration procedure for linearization of the homogenized Richards equations. At each Picard iteration, some degree of multiscale still remains from the intermediate level, so we utilize the generalized multiscale finite element method (GMsFEM) combining with a multi-continuum approach, to upscale the homogenized system to a macroscopic (coarse-grid) level. This scheme involves building uncoupled and coupled multiscale basis functions, which are used not only to construct coarse-grid solution approximation with high accuracy but also (with the coupled multiscale basis) to capture the interactions among continua. These prospects and convergence are demonstrated by several numerical results for the proposed method.

97 MATHEMATICS AND COMPUTING↗

Carbon Dioxide Production in Bedrock Beneath Soils Substantially Contributes to Forest Carbon Cycling

Abstract Soils are widely considered the primary terrestrial organic matter pool mediating carbon transactions with the atmosphere and groundwater. Because soils are both a host and a product of rhizosphere activity, they are thought to mark the location where photosynthetic fixation of carbon dioxide (CO 2 ) is balanced by the oxidation of organic matter. However, in many terrestrial environments, the rhizosphere extends below soils and into fractured bedrock, and it is unknown if the resulting biological and hydrologic dynamics in bedrock have a significant impact on carbon cycling. Here we show substantial production of CO 2 in weathered bedrock at 4–8 m below the thin soils (<0.5 m thick) of a Northern California forest using innovative monitoring technology for sampling gases and water in fractured rock. The deep CO 2 production supports a persistent upward flux of CO 2(g) year‐round from bedrock to soil, constituting between 2% and 29% of the average daily CO 2 efflux from the land surface. When water is rapidly transported across the fractured bedrock vadose zone, nearly 50% of the CO 2 produced in the bedrock dissolves into water, promoting water‐rock interaction and export of dissolved inorganic carbon (DIC) from the unsaturated zone to groundwater, constituting as much as 80% of the DIC exiting the hillslope. Such CO 2 production in weathered bedrock is subject to unique moisture, temperature, biological, and mineralogical conditions which are currently missing from terrestrial carbon models.

54 ENVIRONMENTAL SCIENCES↗

Characterization of Sequentially Extracted Soil Organic Matter by Electrospray Ionization and Atmospheric Pressure Photoionization Fourier Transform Ion Cyclotron Resonance Mass Spectrometry

Soil organic matter (SOM) is a complex mixture of small molecules and biopolymers that are active in various biogeochemical processes. However, the chemical diversity of biopolymer-derived SOM remains poorly explored. Identifying this diversity is important because global environmental changes may well alter SOM chemistry, as field experiments are beginning to show. Here, organic solvent-extractable (DcMe-SOM), base-hydrolyzable (KOH Hy -SOM), and CuO-oxidizable (CuO Ox -SOM) SOM fractions from a forest with a long-term nitrogen addition experiment were sequentially extracted and characterized by Fourier transform ion cyclotron resonance mass spectrometry (FT-ICR MS) coupled with negative-ion electrospray ionization (ESI) or atmospheric pressure photoionization (APPI). From DcMe-SOM to CuO Ox -SOM, the total number of assigned formulas, average O/C ratio, aromaticity, and unsaturation degree of SOM continuously increased, while the average m/z and H/C ratio decreased. Moreover, the dominant chemical category shifted from lipid-like components to phytochemical- and protein-like components. Complementary to ESI, APPI effectively facilitated detection of additional compounds with low polarity. With long-term nitrogen addition, the average m/z, unsaturation degree, aromaticity, and oxidation state of SOM increased, and more aromatic nitrogen-containing formulas were detected in CuOOx-SOM. Our study demonstrates that chronic nitrogen deposition in forests alters both the small molecules and biopolymers of SOM fundamentally.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Development of inter-grid-cell lateral unsaturated and saturated flow model in the E3SM Land Model (v2.0)

Abstract. The lateral transport of water in the subsurface is important in modulating terrestrial water energy distribution. Although a few land surface models have recently included lateral saturated flow within and across grid cells, it is not a default configuration in the Climate Model Intercomparison Project version 6 experiments. In this work, we developed the lateral subsurface flow model within both unsaturated and saturated zones in the Energy Exascale Earth System Model (E3SM) Land Model version 2 (ELMv2.0). The new model, called ELMlat, was benchmarked against PFLOTRAN, a 3D subsurface flow and transport model, for three idealized hillslopes that included a convergent hillslope, divergent hillslope, and tilted V-shaped hillslope with variably saturated initial conditions. ELMlat showed comparable performance against PFLOTRAN in terms of capturing the dynamics of soil moisture and groundwater table for the three benchmark hillslope problems. Specifically, the mean absolute errors (MAEs) of the soil moisture in the top 10 layers between ELMlat and PFLOTRAN were within 1 %±3 %, and the MAEs of water table depth were within ±0.2 m. Next, ELMlat was applied to the Little Washita experimental watershed to assess its prediction of groundwater table, soil moisture, and soil temperature. The spatial pattern of simulated groundwater table depth agreed well with the global groundwater table benchmark dataset generated from a global model calibrated with long-term observations. The effects of lateral groundwater flow on the energy flux partitioning were more prominent in lowland areas with shallower groundwater tables, where the difference in simulated annual surface soil temperature could reach 0.3–0.4 ∘C between ELMv2.0 and ELMlat. Incorporating lateral subsurface flow in ELM improves the representation of the subsurface hydrology, which will provide a good basis for future large-scale applications.

58 GEOSCIENCES↗

Plant organic matter inputs exert a strong control on soil organic matter decomposition in a thawing permafrost peatland

Peatlands are climate critical carbon (C) reservoirs that could become a C source under continued warming. A strong relationship between plant tissue chemistry and the soil organic matter (SOM) that fuels C gas emissions is inferred, but rarely examined at the molecular level. Here we compared Fourier transform infrared (FT-IR) spectroscopy measurements of solid phase functionalities in plants and SOM to ultra-high-resolution mass spectrometric analyses of plant and SOM water extracts across a palsa-bog-fen thaw and moisture gradient in an Arctic peatland. From these analyses we calculated the C oxidation state (NOSC), a measure which can be used to assess organic matter quality. Palsa plant extracts had the highest NOSC, indicating high quality, whereas extracts of Sphagnum, which dominated the bog, had the lowest NOSC. The percentage of plant compounds that are less bioavailable and accumulate in the peat, increases from palsa (25%) to fen (41%) to bog (47%), reflecting the pattern of percent Sphagnum cover. The pattern of NOSC in the plant extracts was consistent with the high number of consumed compounds in the palsa and low number of consumed compounds in the bog. However, in the FT-IR analysis of the solid phase bog peat, carbohydrate content was high implying high quality SOM. Here, we explain this discrepancy as the result of low solubilization of bog SOM facilitated by the low pH in the bog which makes the solid phase carbohydrates less available to microbial decomposition. Plant-associated condensed aromatics, tannins, and lignin-like compounds declined in the unsaturated palsa peat indicating decomposition, but lignin-like compounds accumulated in the bog and fen peat where decomposition was presumably inhibited by the anaerobic conditions. A molecular-level comparison of the aboveground C sources and peat SOM demonstrates that climate-associated vegetation shifts in peatlands are important controls on the mechanisms underlying changing C gas emissions.

54 ENVIRONMENTAL SCIENCES↗