Search NASA⌕ Search

SEARCH · Search NASA

Results for “Residence Times”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Hyporheic Exchange Driven by Submerged Rigid Vegetation: A Modeling Study

Flow-vegetation interaction affects fluid flow hydraulics and associated material transport in river corridors. Concomitant changes in pressure within the flow field due to the presence of vegetation may act as a driver for the formation of hyporheic flow across the sediment-water interface. This potentially important process, however, has yet to be studied. In order to investigate vegetation-induced hyporheic exchange, a series of numerical models of interlinked surface-subsurface flow modified by plant stems was conducted. Periodically staggered plant stem arrays on a flat sediment bed were considered within a coupled multiphysics computational fluid dynamics approach. Plants were idealized as rigid cylinders and arranged in different streamwise and spanwise spacing distances. Each vegetation array was then subjected to a broad range of flow Reynolds Numbers (Re). The results showed that hyporheic flow occurs in all conditions with the presence of vegetation. The vegetation-induced hyporheic flux is found to be a function of Re via a power law. The flux increases with interstem space until the space reaches the distance that rigid stems no longer affect the flow structures in the vicinity of each other. Larger intervegetation distances lead to a larger hyporheic zone. A direct comparison with bedforminduced hyporheic flow showed that vegetation can induce higher hyporheic flux through relatively shallower exchange zones. The results of all the simulations were synthesized into predictive models for hyporheic flux, bulk residence time and exchange depth based on drag coefficient, vegetation density, and Reynolds Number.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Atmospheric Transport of North African Dust-Bearing Supermicron Freshwater Diatoms to South America: Implications for Iron Transport to the Equatorial North Atlantic Ocean

The equatorial North Atlantic Ocean (NAO) is a nutrient-limited ecosystem that relies on the deposition of long-range transported iron (Fe)-containing aerosols to stimulate primary productivity. Using microscopy, we characterized supermicron and supercoarse mode African aerosols transported to the western NAO in boreal winter/spring. In this work, we detected three particle types including African dust, primary biological aerosol particles, and freshwater diatoms (FDs). FDs contained 4% Fe by weight due to surficial dust inclusions that may be susceptible to chemical processing and dissolution. FDs were typically larger than dust particles and comprised 38% of particles between 10 and 18 μm in diameter. The low density, high surface-area-to-volume ratio, and large aspect ratios of FD particles suggest a mechanism by which they can be carried great distances aloft. These same properties likely increase the residence time of FDs in surface waters thereby increasing the time for Fe dissolution and their potential impact on marine biogeochemical cycles.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Effect of Glacial/Interglacial Recharge Conditions on Flow of Meteoric Water Through Deep Orogenic Faults: Insights Into the Geothermal System at Grimsel Pass, Switzerland

Many meteoric-recharged, fault-hosted geothermal systems in amagmatic orogenic belts have been active through the Pleistocene glacial/interglacial climate fluctuations. The effects of such climate-induced recharge variations on fluid flow patterns and residence times of the thermal waters are complex and may influence how the geothermal and mineralization potential of the systems are evaluated. We report systematic thermal–hydraulic simulations designed to reveal the effects of recharge variations, using a model patterned on the orogenic geothermal system at Grimsel Pass in the Swiss Alps. Fault-bound circulation of meteoric water through the Grimsel fault is driven to depths of ~10 km by the high alpine topography. Simulations suggest that the current single pass flow is typical of interglacial periods, during which 1) meteoric recharge into the fault is high (above tens of centimeters per year), 2) conditions are at or somewhat below the critical Rayleigh number, and 3) the hydraulic connectivity along the fault plane is extensive (an extent of at least 10 km into increasingly higher terrain is required to explain the 10 km penetration depth). The subcritical condition constrains the bulk fault permeability to < 1e-14 m2. In contrast, the limited recharge during the numerous Pleistocene glaciation events likely induced a layered flow system with single pass flow confined to shallow depths while non-Rayleigh convection occurred in the deep fault. The same layering can be observed at low aspect ratios (length/depth) of the fault plane, when the available recharge area limits the flux through the fault.

58 GEOSCIENCES↗
DOE OSTI2020

Predictive Inverse Model for Advective Heat Transfer in a Short–Circuited Fracture: Dimensional Analysis, Machine Learning, and Field Demonstration

Identifying fluid flow maldistribution in planar geometries is a well–established problem in subsurface science/engineering. Of particular importance to the thermal performance of enhanced (or “engineered”) geothermal systems is identifying the existence of nonuniform (i.e., heterogeneous) permeability and subsequently predicting advective heat transfer. Here, machine learning via a genetic algorithm (GA) identifies the spatial distribution of an unknown permeability field in a two–dimensional Hele–Shaw geometry (i.e., parallel plates). The inverse problem is solved by minimizing the L2 norm between simulated residence time distribution (RTD) and measurements of an inert tracer breakthrough curve (BTC) (C–Dot nanoparticle). Principal component analysis (PCA) of spatially correlated permeability fields enabled reduction of the parameter space by more than a factor of 10 and restricted the inverse search to reservoir–scale permeability variations. Thermal experiments and tracer tests conducted at the mesoscale Altona Field Laboratory (AFL) demonstrate that the method accurately predicts the effects of extreme flow channeling on heat transfer in a single bedding–plane rock fracture. However, this is only true when the permeability distributions provide adequate matches to both tracer RTD and frictional pressure loss. Without good agreement to frictional pressure loss, it is still possible to match a simulated RTD to measurements, but subsequent predictions of heat transfer are grossly inaccurate. Here, the results of this study suggest that it is possible to anticipate the thermal effects of flow maldistribution, but only if both simulated RTDs and frictional pressure loss between fluid inlets and outlets are in good agreement with measurements.

42 ENGINEERING↗
DOE OSTI2022

The Relation of Environmental Conditions With Charge Structure in Central Argentina Thunderstorms

In this study we explored the environmental conditions hypothesized to induce a dominant charge structure in thunderstorms in the province of Cordoba, Argentina, during the RELAMPAGO-CACTI (Remote sensing of Electrification, Lightning, And Mesoscale/microscale Processes with Adaptive Ground Observations-Clouds, Aerosols, Complex Terrain Interactions) field campaigns. Hypothesized environmental conditions are thought to be related to small warm cloud residence time and warm rain growth suppression, which lead to high cloud liquid water contents in the mixed-phase zone, contributing to positive charging of graupel and anomalous charge structure storms. Data from radiosondes, a cloud condensation nuclei (CCN) ground-based instrument and reanalysis were used to characterize the proximity inflow air of storms with anomalous and normal charge structures. Consistent with the initial hypothesis, anomalous storms had small warm cloud depth caused by dry low-level humidity and low 0°C height. Anomalous storms were associated with lower CCN concentrations than normal storms, an opposite result to the initial expectation. High CAPE is not an important condition for the development of anomalous storms in Argentina, as no clear pattern could be found among the different parameters calculated for updraft proxy that would be consistent with the initial hypothesis.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

How Does LCL Height Influence Deep Convective Updraft Width?

Previous studies have hypothesized that the width of deep convection should positively scale with lifting condensation level (LCL) height. To evaluate this hypothesis, we analyzed idealized large-eddy simulations with varying LCL heights and initial warm bubble widths in unsheared environments with comparable convective available potential energy. For a given initial warm bubble width, simulations with higher LCLs result in wider, deeper, and stronger cloudy updrafts compared to simulations with lower LCLs. Rising dry thermals in higher LCL simulations experience longer residence times within the sub-cloud layer, and consequently entrain more conditionally unstable air and grow wider before reaching the LCL. The resulting cloudy updrafts are wider, deeper, and have faster vertical velocities because of a reduction in entrainment-driven dilution of buoyancy, relative to lower LCL simulations. These results confirm the hypothesized positive relationship between LCL height and deep convective updraft width, and provide a physical explanation for this relationship.

58 GEOSCIENCES↗
DOE OSTI2021

Hydrate Formation on Marine Seep Bubbles and the Implications for Water Column Methane Dissolution

Abstract Methane released from seafloor seeps contributes to a number of benthic, water column, and atmospheric processes. At seafloor seeps within the methane hydrate stability zone, crystalline gas hydrate shells can form on methane bubbles while the bubbles are still in contact with the seafloor or as the bubbles begin ascending through the water column. These shells reduce methane dissolution rates, allowing hydrate‐coated bubbles to deliver methane to shallower depths in the water column than hydrate‐free bubbles. Here, we analyze seafloor videos from six deepwater seep sites associated with a diverse range of bubble‐release processes involving hydrate formation. Bubbles that grow rapidly are often hydrate‐free when released from the seafloor. As bubble growth slows and seafloor residence time increases, a hydrate coating can form on the bubble's gas‐water interface, fully coating most bubbles within ∼10 s of the onset of hydrate formation at the seafloor. This finding agrees with water‐column observations that most bubbles become hydrate‐coated after their initial ∼150 cm of rise, which takes about 10 s. Whether a bubble is coated or not at the seafloor affects how much methane a bubble contains and how quickly that methane dissolves during the bubble's rise through the water column. A simplified model shows that, after rising 150 cm above the seafloor, a bubble that grew a hydrate shell before releasing from the seafloor will have ∼5% more methane than a bubble of initial equal volume that did not grow a hydrate shell after it traveled to the same height.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Examining conditions supporting the development of anomalous charge structures in supercell thunderstorms in the Southeastern United States

Hypotheses regarding favorable conditions for anomalous charging have primarily resulted from studies within the Great Plains region of the United States, where the efficiency of warm precipitation processes is thought to be fundamental. Rare observations of anomalous charge structures in the Southeastern region challenge existing conceptual models used to explain anomalous charging. As a rigorous evaluation of conditions that support anomalous charge structures, environmental characteristics and bulk kinematic and microphysical properties of two normal and two anomalous supercell thunderstorms observed in the Southeast were compared in this paper. Within the anomalous supercells, greater quantities of precipitation ice were identified at higher altitudes and colder temperatures, suggesting a greater depth of riming growth and increased vertical transport of rimed hydrometeors. Deeper anomalous supercell updrafts were larger and stronger in the upper mixed-phase and glaciated regions. However, normal supercells were characterized by more robust low-level updrafts, resulting in comparable warm cloud residence times that suggested warm precipitation processes were not necessarily less efficient in the anomalous supercells. Indications of enhanced mixed-phase liquid water content in favor of anomalous charging were observed in the anomalous supercells, though contrasts in related environmental parameters were not as large as observed in other comparative studies. Anomalous supercell environments were characterized by increased instability, shallower warm cloud depth, as well as lower relative humidity in the 700 mb to 500 mb layer. Evidence of impacts from dry air in anomalous storm structures suggested that water vapor content may have affected particle-scale charge transfer in support of anomalous charge structure development.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Experimental Determination of the Relationship Between Organic Aerosol Viscosity and Ice Nucleation at Upper Free Tropospheric Conditions

Abstract This work presents measurements of the ice nucleating ability of secondary organic material recorded between −40 and −70°C and relative humidity with respect to ice (RH ice ) between 150% and 220%. For a subset of systems, temperature and humidity dependence of particle viscosity as well as dry glass transition temperature were characterized using the dual tandem differential mobility analyzer method. Eleven unique monoterpene, sesquiterpene and aromatic precursors were used to generate secondary organic material (SOM) using either an oxidation flow reactor (OFR) or an environmental chamber (EC). For the SOM for which viscosity was measured, the particle glass transition temperatures varied between 6 and 23°C ( n = 8). Measurements were performed to verify that increased relative humidity did not plasticize the particles below −10°C at residence times similar to those in the ice nucleation instrument. No heterogeneous ice nucleation was observed at the ∼0.5% onset threshold for any of the materials generated. The ice nucleation occurs by the freezing of SOM solution droplets consistent with homogeneous freezing indicating that they form an aqueous solution, or the SOM particles required water saturation to freeze, indicating that they were hydrophobic. Experiments exploring the influence of functional groups and mass loading did not reveal any obvious influence of particle chemistry or generation conditions on the results. Close structural matches between known organic ice nucleating particles as precursor or formed products did not yield materials that promoted freezing. These experiments suggest that heterogeneous ice formation of glassy secondary organic materials is likely uncommon under upper free tropospheric conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Spatial Microbial Respiration Variations in the Hyporheic Zones Within the Columbia River Basin

Abstract While the hyporheic zone (HZ) accounts for a significant portion of whole stream CO 2 concentrations, HZ respiration modeling studies are lacking in quantifying their contributions to the total CO 2 at large watershed/basin scales. Quantifying the contribution of anaerobic respiration is also underappreciated. This study used a carbon‐nitrogen‐coupled river corridor model to quantify HZ aerobic and anaerobic respiration and determined the key factors controlling their spatial variability within the Columbia River Basin (CRB). The modeled respiration patterns showed high spatial variability. Among the nine sub‐basins composing the CRB, the Lower Columbia and the Willamette, which receive higher precipitation, had higher respiration. Medium‐sized rivers (fourth to sixth orders) produced the highest aerobic and anaerobic respiration among reaches of different sizes. At the basin scale, aerobic respiration is dominant, representing approximately 98.7% of the total respiration across the CRB. While most of the reaches were dominant with aerobic respiration, reaches in agricultural land showed a relatively higher anaerobic respiration (18%) ratio. A variable importance analysis showed that hyporheic exchange flux controlled most of the spatial variability of HZ respiration, dominating over other physical variables such as residence time, stream dissolved organic carbon (DOC), nitrate, and dissolved oxygen (DO). The influence of substrate concentration (DOC and DO) is larger in modeling anaerobic respiration than aerobic respiration. Future efforts will focus on improving the estimation of the HZ exchange flux and the implementation of spatially explicit parameterizations for the reactions of interest to reduce model uncertainty.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Local and Global Sensitivity Analysis of a Reactive Transport Model Simulating Floodplain Redox Cycling

Reactive transport models (RTMs) are essential tools that simulate the coupling of advective, diffusive, and reactive processes in the subsurface, but their complexity makes them difficult to understand, develop and improve without accompanying statistical analyses. Although global sensitivity analysis (SA) can address these issues, the computational cost associated with most global SA techniques limits their use with RTMs. In this study, we apply distance-based generalized sensitivity analysis (DGSA), a novel and computationally efficient method of global SA, to a floodplain-scale RTM and compare DGSA results to those from local SA. Our test case focuses on the impact of 17 uncertain environmental parameters on spatially and temporally variable redox conditions within a floodplain aquifer. The input parameters considered include flow and diffusion rates, geochemical reaction rates, and the spatial distribution of sediment facies. Sensitivity was evaluated for three distinct components of the model response, encompassing both multidimensional and categorical output. Parameter rankings differ between local SA and DGSA, due to nonlinear effects of individual parameters and interaction effects between parameters. DGSA results show that fluid residence time, which is controlled by aquifer permeability, generally exerts a stronger control on redox conditions than do geochemical reaction rates. Sensitivity indices also demonstrate that sulfate reduction is key for establishing and maintaining reducing conditions throughout the aquifer. Furthermore, these results provide insights into the key drivers of heterogeneous redox processes within floodplain aquifers, as well as the main sources of uncertainty when modeling complex subsurface systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Organizational Principles of Hyporheic Exchange Flow and Biogeochemical Cycling in River Networks Across Scales

Hyporheic zones increase freshwater ecosystem resilience to hydrological extremes and global environmental change. However, current conceptualizations of hyporheic exchange, residence time distributions, and the associated biogeochemical cycling in streambed sediments do not always accurately explain the hydrological and biogeochemical complexity observed in streams and rivers. Specifically, existing conceptual models insufficiently represent the coupled transport and reactivity along groundwater and surface water flow paths, the role of autochthonous organic matter in streambed biogeochemical functioning, and the feedbacks between surface-subsurface ecological processes, both within and across spatial and temporal scales. While simplified approaches to these issues are justifiable and necessary for transferability, the exclusion of important hyporheic processes from our conceptualizations can lead to erroneous conclusions and inadequate understanding and management of interconnected surface water and groundwater environments. This is particularly true at the landscape scale, where the organizational principles of spatio-temporal dynamics of hyporheic exchange flow (HEF) and biogeochemical processes remain largely uncharacterized. This article seeks to identify the most important drivers and controls of HEF and biogeochemical cycling based on a comprehensive synthesis of findings from a wide range of river systems. We use these observations to test current paradigms and conceptual models, discussing the interactions of local-to-regional hydrological, geomorphological, and ecological controls of hyporheic zone functioning. This improved conceptualization of the landscape organizational principles of drivers of HEF and biogeochemical processes from reach to catchment scales will inform future river research directions and watershed management strategies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Impact of Coastal Marsh Eco‐Geomorphologic Change on Saltwater Intrusion Under Future Sea Level Rise

Abstract Coastal saltwater intrusion (SWI) is one key factor that affects the hydrology, ecology, and biogeochemistry of coastal ecosystems. Future climate change, especially intensified sea level rise (SLR), is expected to trigger SWI to encroach on coastal freshwater aquifers more intensively. Numerous studies have investigated decadal/century scale SWI under SLR by assuming a static coastal landscape topography. However, coastal landscapes are highly dynamic in response to SLR, and the impact of coastal landscape evolution on SWI has received very little attention. Thus, this study used a coastal marsh landscape as an example and investigated how coastal marsh evolution affects future SWI with a physically‐based coastal hydro‐eco‐geomorphologic model, Advanced Terrestrial Simulator. Our numerical experiments showed that it is very likely that the marsh elevation increases with future SLR due to sediment deposition, and a depression zone is formed due to different marsh accretion rates between the ocean boundary and the inland. We found that marsh accretion may significantly reduce the surface saltwater inflow at the ocean boundary, and the evolved topographic depression zone may prolong the residence time of surface ponded saltwater, affecting subsurface salinity distribution differently. We also predicted that marshlands may become more sensitive to upland freshwater supply under future SLR, compared with previous predictions without marsh evolution. This study demonstrates the importance of coastal evolution to coastal freshwater‐saltwater interaction. The eco‐geomorphologic effect may not be ignored when evaluating coastal SWI under SLR at decadal or century scales.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Arctic Permafrost Thawing Enhances Sulfide Oxidation

Permafrost degradation is altering biogeochemical processes throughout the Arctic. Thaw-induced changes in organic matter transformations and mineral weathering reactions are impacting fluxes of inorganic carbon (IC) and alkalinity (ALK) in Arctic rivers. However, the net impact of these changing fluxes on the concentration of carbon dioxide in the atmosphere (pCO 2 ) is relatively unconstrained. Resolving this uncertainty is important as thaw-driven changes in the fluxes of IC and ALK could produce feedbacks in the global carbon cycle. Enhanced production of sulfuric acid through sulfide oxidation is particularly poorly quantified despite its potential to remove ALK from the ocean-atmosphere system and increase pCO 2 , producing a positive feedback leading to more warming and permafrost degradation. In this work, we quantified weathering in the Koyukuk River, a major tributary of the Yukon River draining discontinuous permafrost in central Alaska, based on water and sediment samples collected near the village of Huslia in summer 2018. Using measurements of major ion abundances and sulfate (SO 4 2- ) sulfur ( 34 S/ 32 S) and oxygen ( 18 O/ 16 O) isotope ratios, we employed the MEANDIR inversion model to quantify the relative importance of a suite of weathering processes and their net impact on pCO 2 . Calculations found that approximately 80% of SO 4 2- in mainstem samples derived from sulfide oxidation with the remainder from evaporite dissolution. Moreover, 34 S/ 32 S ratios, 13 C/ 12 C ratios of dissolved IC, and sulfur X-ray absorption spectra of mainstem, secondary channel, and floodplain pore fluid and sediment samples revealed modest degrees of microbial sulfate reduction within the floodplain. Weathering fluxes of ALK and IC result in lower values of pCO 2 over timescales shorter than carbonate compensation (~10 4 yr) and, for mainstem samples, higher values of pCO 2 over timescales longer than carbonate compensation but shorter than the residence time of marine (~10 7 yr). Furthermore, the absolute concentrations of and Mg 2+ in the Koyukuk River, as well as the ratios of SO 4 2- and Mg 2+ to other dissolved weathering products, have increased over the past 50 years. Through analogy to similar trends in the Yukon River, we interpret these changes as reflecting enhanced sulfide oxidation due to ongoing exposure of previously frozen sediment and changes in the contributions of shallow and deep flow paths to the active channel. Overall, these findings confirm that sulfide oxidation is a substantial outcome of permafrost degradation and that the sulfur cycle responds to permafrost thaw with a timescale-dependent feedback on warming.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Drivers of Phytoplankton Variability in and Near the Pearl River Estuary, South China Sea During Typhoon Hato (2017): A Numerical Study

The Pearl River estuary-coastal system in the Northern South China Sea is frequently affected by tropical cyclones (TCs) every year. Previous in-situ and remote sensing studies have found typhoons can enhance phytoplankton biomass and induce blooms in this region. However, the mechanistic links between phytoplankton blooms and typhoons have not been well elucidated due to the interplay of multiple processes along the land-ocean-atmosphere interfaces. Unravelling these interactions will have to rely on an integrated modelling system. In this work, we constructed a realistic, 3-dimensional, land-ocean-atmosphere modelling system with the marine ecosystem and sediment components for the China Great Bay Area. By using the integrated modelling system, we quantitatively investigated phytoplankton response to hydrological conditions variations under Typhoon Hato (2017), a strong typhoon case. Passive tracer experiments showed that with high river discharge induced by heavy rainfall, the residence time of Lingding Bay is as short as 15-day, less than half of that under the climatological discharge. The increase in freshwater pulse washed out the phytoplankton biomass within Lingding Bay. While for the offshore region, the source and sink terms analysis showed that the increase of phytoplankton biomass in the first week was because of the uplift of nutrient-rich subsurface water, while in the second week was because of the seaward propagated nearshore high phytoplankton biomass water. While riverine nutrients support phytoplankton growth in the third week, a large part of phytoplankton biomass was lost to zooplankton grazing, showing the system shifted from the bottom-up control to the top-down control.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Matrix Approach to Land Carbon Cycle Modeling

Land ecosystems contribute to climate change mitigation by taking up approximately 30% of anthropogenically emitted carbon. However, estimates of the amount and distribution of carbon uptake across the world's ecosystems or biomes display great uncertainty. The latter hinders a full understanding of the mechanisms and drivers of land carbon uptake, and predictions of the future fate of the land carbon sink. The latter is needed as evidence to inform climate mitigation strategies such as afforestation schemes. To advance land carbon cycle modeling, we have developed a matrix approach. Land carbon cycle models use carbon balance equations to represent carbon exchanges among pools. Our approach organizes this set of equations into a single matrix equation without altering any processes of the original model. The matrix equation enables the development of a theoretical framework for understanding the general, transient behavior of the land carbon cycle. While carbon input and residence time are used to quantify carbon storage capacity at steady state, a third quantity, carbon storage potential, integrates fluxes with time to define dynamic disequilibrium of the carbon cycle under global change. The matrix approach can help address critical contemporary issues in modeling, including pinpointing sources of model uncertainty and accelerating spin-up of land carbon cycle models by tens of times. The accelerated spin-up liberates models from the computational burden that hinders comprehensive parameter sensitivity analysis and assimilation of observational data to improve model accuracy. Such computational efficiency offered by the matrix approach enables substantial improvement of model predictions using ever-increasing data availability. Overall, the matrix approach offers a step change forward for understanding and modeling the land carbon cycle.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

The Influence of Air‐Sea CO 2 Disequilibrium on Carbon Sequestration by the Ocean's Biological Pump

Abstract The ocean's biological carbon pump (BCP) affects the Earth's climate by sequestering CO 2 away from the atmosphere for decades to millennia. One primary control on the amount of carbon sequestered by the biological pump is air‐sea CO 2 disequilibrium, which is controlled by the rate of air‐sea CO 2 exchange and the residence time of CO 2 in surface waters. Here, we use a data‐assimilated model of the soft tissue BCP to quantify carbon sequestration inventories and time scales from remineralization in the water column to equilibration with the atmosphere. We find that air‐sea CO 2 disequilibrium enhances the global biogenic carbon inventory by ∼35% and its sequestration time by ∼70 years compared to identical calculations made assuming instantaneous air‐sea CO 2 exchange. Locally, the greatest enhancement occurs in the subpolar Southern Ocean, where air‐sea disequilibrium increases sequestration times by up to 600 years and the biogenic dissolved inorganic carbon inventory by >100% in the upper ocean. Contrastingly, in deep‐water formation regions of the North Atlantic and Antarctic margins, where biological production creates undersaturated surface waters which are subducted before fully equilibrating with the atmosphere, air‐sea CO 2 disequilibrium decreases the depth‐integrated sequestration inventory by up to ∼150%. The global enhancement of carbon sequestration by air‐sea disequilibrium is particularly important for carbon respired in deep waters that upwell in the Southern Ocean. These results highlight the importance of accounting for air‐sea CO 2 disequilibrium when evaluating carbon sequestration by the biological pump and for assessing the efficacy of ocean‐based CO 2 removal methods.

Environmental Sciences & Ecology↗
DOE OSTI2023

Effects of Dead‐End Fractures on Non‐Fickian Transport in Three‐Dimensional Discrete Fracture Networks

Abstract Understanding mechanistic causes of non‐Fickian transport in fractured media is important for many hydrogeologic processes and subsurface applications. This study elucidates the effects of dead‐end fractures on non‐Fickian transport in three‐dimensional (3D) fracture networks. Although dead‐end fractures have been identified as low‐velocity regions that could delay solute transport, the direct relation between dead‐end fractures and non‐Fickian transport has been elusive. We systematically generate a large number of 3D discrete fracture networks with different fracture length distributions and fracture densities. We then identify dead‐end fractures using a novel graph‐based method. The effect of dead‐end fractures on solute residence time maximizes at the critical fracture density of the percolation threshold, leading to strong late‐time tailing. As fracture density increases beyond the percolation threshold, the network connectivity increases, and dead‐end fractures diminish. Consequently, the increase in network connectivity leads to a reduction in the degree of late‐time tailing. We also show that dead‐end fractures can inform about main transport paths, such as the mean tortuosity of particle trajectories. This study advances our mechanistic understanding of solute transport in 3D fracture networks.

58 GEOSCIENCES↗