Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “cumulative”

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 181 records · Page 10

Deep learning model for fast, science-based forecasting of fluid migration along faults in geologic carbon storage scenarios

Effective long-term geologic storage depends on robust site selection and credible, science-based forecasting of subsurface behavior to ensure storage integrity. For this work, we develop a deep learning–based reduced-order model (ROM) to quantify potential carbon dioxide (CO₂) and brine migration through geological faults. The ROM combines a Transformer model for binary classification and a Stacked Ensemble for regression, trained on a comprehensive dataset generated from 1400 physics-based reservoir simulations. Key geologic and operational parameters—including fault geometry, reservoir structure, and injection conditions—were systematically varied to capture a wide range of fluid migration scenarios. The ROM accurately predicts the onset of migration, cumulative migration volumes of both CO₂ and brine, and associated migration rates, as compared to an independent set of validation simulations, while significantly reducing computational cost compared to traditional simulation methods. Model performance was evaluated across diverse fault configurations, revealing that shallow reservoir geometry and fault angle are among the most influential factors governing migration behavior. Sensitivity analysis using SHapley Additive exPlanations (SHAP) provided interpretability, revealing distinct patterns in how geological and operational features drive transient versus cumulative migration outcomes. The ROM’s ability to rapidly simulate fault migration scenarios enables efficient sensitivity analyses, scenario evaluations, and decision support for site selection and monitoring design. This approach enhances the safety, scalability, and long-term operational performance of geologic carbon storage (GCS) systems by providing a robust, interpretable tool for predicting subsurface fluid migration and assessing fault-related migration potential.

42 ENGINEERING↗

Spinal Cord Dose Tolerance to Stereotactic Body Radiation Therapy

Spinal cord tolerance data for stereotactic body radiation therapy (SBRT) were extracted from published reports, reviewed, and modelled. For de novo SBRT delivered in 1 to 5 fractions, the following spinal cord point maximum doses (D{sub max}) are estimated to be associated with a 1% to 5% risk of radiation myelopathy (RM): 12.4 to 14.0 Gy in 1 fraction, 17.0 Gy in 2 fractions, 20.3 Gy in 3 fractions, 23.0 Gy in 4 fractions, and 25.3 Gy in 5 fractions. For reirradiation SBRT delivered in 1 to 5 fractions, reported factors associated with a lower risk of RM include cumulative thecal sac equivalent dose in 2 Gy fractions with an alpha/beta of 2 (EQD2{sub 2}) D{sub max} ≤70 Gy; SBRT thecal sac EQD2{sub 2} D{sub max} ≤25 Gy, thecal sac SBRT EQD2{sub 2} D{sub max} to cumulative EQD2{sub 2} D{sub max} ratio ≤0.5, and a minimum time interval to reirradiation of ≥5 months. Larger studies containing complete institutional cohorts with dosimetric data of patients treated with spine SBRT, with and without RM, are required to refine RM risk estimates.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Dose Summation Strategies for External Beam Radiation Therapy and Brachytherapy in Gynecologic Malignancy: A Review from the NRG Oncology and NCTN Medical Physics Subcommittees

Definitive, nonsurgical management of gynecologic malignancies involves external beam radiation therapy (EBRT) and/or brachytherapy (BT). Summation of the cumulative dose is critical to assess the total biologic effective dose to targets and organs at risk. Cumulative dose calculation from EBRT and BT can be performed with or without image registration (IR) and biologic dose summation. Among these dose summation strategies, linear addition of dose-volume histogram (DVH) parameters without IR is the global standard for composite dose reporting. This approach stems from an era without image guidance and simple external beam and brachytherapy treatment approaches. With technological advances, EBRT and high-dose-rate BT have evolved to allow for volume-based treatment planning and delivery. Modern conformal therapeutic radiation involves volumetric or intensity modulated EBRT, capable of simultaneously treating multiple targets at different specified dose levels. Therefore, given the complexity of modern radiation treatment, the linear addition of DVH parameters from EBRT and high-dose-rate BT is challenging to represent the combined dose distribution. Deformable image registration (DIR) between EBRT and image guided brachytherapy (IGBT) data sets may provide a more nuanced calculation of multimodal dose accumulation. However, DIR is still nascent in this regard, and needs further development for accuracy and efficiency for clinical use. Biologic dose summation can combine physical dose maps from EBRT and each IGBT fraction, thereby generating a composite DVH from the biologic effective dose. However, accurate radiobiologic parameters are tissue-dependent and not well characterized. A combination of voxel-based DIR and biologic weighted dose maps may be the best approximation of dose accumulation but remains invalidated. The purpose of this report is to review dose summation strategies for EBRT and BT, including conventional equivalent dose in 2-Gy fractions dose summation without image registration, physical dose summation using 3-dimensional rigid IR and DIR, and biologic dose summation. We also provide general clinical workflows for IGBT with a focus on cervical cancer.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

An arbitrarily high-order three-dimensional Cartesian-grid method for reconstructing interfaces from volume fraction fields

Here Tthis work describes a newly developed, arbitrarily high-order Cartesian-grid method for reconstructing material interfaces from a volume fraction field. The method begins by identifying all of the grid cells in the volume fraction field that are intersected by the interface and need to be approximated by the reconstruction scheme. Finite-differences are used to calculate the gradient of the volume fraction field and provide an estimate of the surface normal in all of the interfacial grid cells. Groups of connected grid cells are then identified which all have the same dominant component of the normal vector. This grouping by orientation determines the proper dependent variable to use in the surface reconstruction (e.g. for a 2D curve, this step determines if the surface will be approximated by a function of x or y). A cumulative integral over the surface is constructed and fit using b-splines for two-dimensional problems or tensor-product b-splines for three-dimensional problems. This construction allows for the interface to be recovered through application of the second fundamental theorem of calculus. Fitting the cumulative integral with $\mathscr{N}$ th-order b-splines (or tensor-product b-splines) yields an ($\mathscr{N}$-1) th-order convergence rate of the interface shape. Differentiation of the b-spline interface function(s) allows for the high-order approximation of the normal vector and curvature to be obtained directly anywhere along b-spline. Together, the proposed reconstruction technique can achieve arbitrarily high mesh convergence rates. Validation tests are presented with mesh convergence rates ranging from fourth- to tenth-order.

97 MATHEMATICS AND COMPUTING↗

Effects of residue removal and tillage on greenhouse gas emissions in continuous corn systems as simulated with RZWQM2

Agricultural production is a major source of carbon dioxide (CO 2 ) and nitrous oxide (N 2 O) globally. The effects of conservation practices on soil CO 2 and N 2 O emissions remain a high degree of uncertainty. In this study, soil CO 2 and N 2 O emissions under different residue and tillage practices in an irrigated, continuous corn system, were investigated using the Root Zone Water Quality Model (RZWQM2). Combinations of no/high stover removal (NR and HR, respectively) and no-till/conventional tillage (NT and CT, respectively) field experiments were tested over the four crop-years (Apr. 2011–Apr. 2015). The model was calibrated using the NRCT, and validated with other treatments. The simulation results showed that soil volumetric water content (VWC) in the NR treatments (i.e., NRCT and NRNT) was 1.3%–1.9% higher than that in the HR treatments (i.e., HRCT and HRNT) averaged across the four years. A higher amount of CO 2 and N 2 O emissions were simulated in the NRCT across the four years (annual average: 7034 kg C/ha/yr for CO 2 and 3.8 kg N/ha/yr for N 2 O), and lower emissions were in the HRNT (annual average: 6329 kg C/ha/yr and 3.7 kg N/ha/yr for N 2 O). A long-term simulation (2001–2015) suggested that the CO 2 and N 2 O emissions were closely correlated with the stover removal degree (SRD), tillage, VWC, soil temperature (ST), years in management (Y), and fertilizer application. Stover and tillage practices had cumulative effects on CO 2 emissions. The simulated annual CO 2 emissions in 1st year from NRCT, NRNT, and HRCT were 7.8%, 0.0%, and 7.7% higher than that from HRNT, respectively; then the emissions in 15th year were 63.6%, 47.7%, and 29.1% higher, respectively. Meanwhile, there were no cumulative effects on N 2 O emissions. Furthermore, the results also demonstrated that the RZWQM2 is a promising tool for evaluating the long-term effects of CO 2 and N 2 O emissions on different conservation practices.

54 ENVIRONMENTAL SCIENCES↗

A robust deep learning workflow to predict multiphase flow behavior during geological C O 2 sequestration injection and Post-Injection periods

Simulation of multiphase flow in porous media is essential to manage the geologic CO 2 sequestration (GCS) process, and physics-based simulation approaches usually take prohibitively high computational cost due to the nonlinearity of the coupled physics. This paper contributes to the development and evaluation of a deep learning workflow that accurately and efficiently predicts the temporal-spatial evolution of pressure and CO 2 plumes during injection and post-injection periods of GCS operations. Based on a Fourier Neural Operator, the deep learning workflow takes input variables or features including rock properties, well operational controls and time steps, and predicts the state variables of pressure and CO 2 saturation. To further improve the predictive fidelity, separate deep learning models are trained for CO 2 injection and post-injection periods due to the difference in primary driving force of fluid flow and transport during these two phases. We also explore different combinations of features to predict the state variables. We use a realistic example of CO 2 injection and storage in a 3D heterogeneous saline aquifer, and apply the deep learning workflow that is trained from physics-based simulation data and emulate the physics process. Through this numerical experiment, we demonstrate that using two separate deep learning models to distinguish post-injection from injection period generates the most accurate prediction of pressure, and a single deep learning model of the whole GCS process including the cumulative injection volume of CO 2 as a deep learning feature, leads to the most accurate prediction of CO 2 saturation. For the post-injection period, it is key to use cumulative CO 2 injection volume to inform the deep learning models about the total carbon storage when predicting either pressure or saturation. The deep learning workflow not only provides high predictive fidelity across temporal and spatial scales, but also offers a speedup of 250 times compared to full physics reservoir simulation, and thus will be a significant predictive tool for engineers to manage the long-term process of GCS.

58 GEOSCIENCES↗

A model for trapping and re-solution regarding intra-granular bubbles in UO 2 , linked to atomic-scale simulations

In the literature, a clear definition of the irradiation re-solution frequency of gas from bubbles in the UO 2 fuel is absent. Moreover, for intra-granular bubbles, a detailed calculation of the cumulated displaced gas quantities in function of the distance from the radius of the bubble after a re-solution event has never been published. The assessment of these two elements is very useful if we want to increase the adherence of fission gas release codes to our present knowledge of the behavior of fission gases. Hence, we suggest to link the definition of the re-solution frequency to atomic-scale simulations. Furthermore, we present the cumulated displaced gas quantities obtained from Molecular Dynamics calculations, from which we have derived a re-solution profile that can be exploited to better consider the irradiation re-solution phenomenon inside Fission Gas Release codes. On top of that, we have built a new trapping/re-solution model for intra-granular bubbles linked to Molecular Dynamics simulations that can be easily incorporated into Fission Gas Release codes. In conclusion, we also check that the model is properly built through the comparison of the new model against a reference.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Aperture size distribution, length, and preferential location of bed-parallel veins in shale

Bed-parallel, calcite-filled veins (BPVs) are common in shale formations, and although they have been widely described in other studies, little is known about their population aperture size distribution. To address this knowledge gap, we analyzed BPV sizes in outcrops and cores from the Vaca Muerta Formation, Neuquén Basin, Argentina; in two cores from the Marcellus Formation, Appalachian Basin, northeast Pennsylvania; and in one core from the Wolfcamp Shale, Delaware Basin, West Texas. Nine out of ten aperture size populations follow a negative exponential distribution, with one following a weak power law. Bed-parallel vein size distribution and intensity vary among formations and within the same shale. We define three groups of distributions: (1) Vaca Muerta outcrops, with the highest BPV intensity and the largest BPVs (cumulative frequency of 4.9 BPVs per meter [BPVs/m] for apertures 0.265 mm to 8.7 cm); (2) Vaca Muerta cores with a similar BPV intensity overall but with no apertures wider than 1.2 cm; and (3) Vaca Muerta, Wolfcamp, and Marcellus cores with the fewest BPVs (cumulative frequency up to 0.63 BPVs/m) and very few wider than 1 cm. Aperture and length in two outcrop data sets are weakly positively correlated and follow power laws with exponents of 0.44 and 0.49. Mechanical interfaces at boundaries between different lithologies exert a strong control on BPV location, with 65–75% of observed interfaces having BPVs along them. Only 25–30% of the BPVs occur at observed material interfaces, however, and unless subtle, unobserved mechanical layering is present, other factors must also control location. BPV intensity and organic richness (TOC) from Vaca Muerta well logs are correlated in some instances but not in others, indicating TOC is not always a good proxy for BPV location or intensity. Furthermore, these findings provide useful information for modeling of hydraulic fracture treatments where BPVs may influence development of the stimulated fracture network, for example by limiting height growth.

58 GEOSCIENCES↗

IRDFF-II: A New Neutron Metrology Library

High quality nuclear data is the most fundamental underpinning for all neutron metrology applications. This paper describes the release of version II of the International Reactor Dosimetry and Fusion File (IRDFF-II) that contains a consistent set of nuclear data for fission and fusion neutron metrology applications up to 60 MeV neutron energy. The library is intended to support: a) applications in research reactors; b) safety and regulatory applications in the nuclear power generation in commercial fission reactors; and c) material damage studies in support of the research and development of advanced fusion concepts. The paper describes the contents of the library, documents the thorough verification process used in its preparation, and provides an extensive set of validation data gathered from a wide range of neutron benchmark fields. The new IRDFF-II library includes 119 metrology reactions, four cover material reactions to support self-shielding corrections, five metrology metrics used by the dosimetry community, and cumulative fission products yields for seven fission products in three different neutron energy regions. In support of characterizing the measurement of the residual nuclei from the dosimetry reactions and the fission product decay modes, the present document lists the recommended decay data, particle emission energies and probabilities for 68 activation products. It also includes neutron spectral characterization data for 29 neutron benchmark fields for the validation of the library contents. Additional six reference fields were assessed (four from plutonium critical assemblies, two measured fields for thermal-neutron induced fission on 233U and 239Pu targets) but not used for validation due to systematic discrepancies in C/E reaction rate values or lack of reaction-rate experimental data. Another ten analytical functions are included that can be useful for calculating average cross sections, average energy, thermal spectrum average cross sections and resonance integrals. The IRDFF-II library and comprehensive documentation is available online at www-nds.iaea.org/IRDFF/. Evaluated cross sections can be compared with experimental data and other evaluations at www-nds.iaea.org/exfor/endf.htm. The new library is expected to become the international reference in neutron metrology for multiple applications.

IAEA IRDFF↗

Survey of tungsten gross erosion from main plasma facing components in WEST during a L-mode high fluence campaign

An initial high fluence campaign was performed in WEST, in 2023, on the newly installed actively cooled tungsten divertor composed of ITER-grade monoblocks. The campaign consisted in the repetition of a 60 s long Deuterium L-mode pulse in attached divertor conditions, cumulating over 10000s of plasma exposure. A maximum deuterium fluence of approximately 5⋅1⁢026 m−2 was reached in the outer strike point region, representative of a few high performance ITER pulses. Gross tungsten erosion inferred from visible spectroscopy shows that the most eroded plasma facing component is the inner divertor target with rates ten times larger than on the outer divertor target. The outer midplane tungsten bumpers, located a few centimeters from the plasma, show gross erosion rates two times lower than at the outer divertor. We conclude that the outer midplane bumpers have a negligible contribution to the long range tungsten migration and deposition onto the lower divertor. The cumulated gross erosion rate on the inner divertor translates in an effective gross erosion thickness of about 20μ⁢m, while it is about 2μ⁢m for the outer divertor. Strikingly, these orderings coincide with the thickness of deposits found locally on the divertor: the exposed surfaces of high field side monoblocks are covered with several tens of μ⁢m tungsten deposits, while on the lower field side, few μ⁢m thin tungsten deposits are only found on the magnetically shadowed parts of monoblocks. The strong impact of those deposits on WEST operation, namely perturbation of surface temperature measurement with infra-red thermography, and the emission of flakes causing radiative perturbation of the confined plasma, calls for anticipating similar issues in ITER. In particular, the start of research operation shall consider the definition of a divertor erosion budget in order to anticipate the formation of deleterious deposits.

Fedorczak, N.↗

Pseudorapidity dependence of anisotropic flow and its decorrelations using long-range multiparticle correlations in Pb–Pb and Xe–Xe collisions

The pseudorapidity dependence of elliptic (v 2 ), triangular (v 3 ), and quadrangular (v 4 ) flow coefficients of charged particles measured in Pb–Pb collisions at a centre-of-mass energy per nucleon pair of $\sqrt{S_{NN}}$ = 5.02 TeV and in Xe–Xe collisions at $\sqrt{S_{NN}}$ = 5.44 TeV with ALICE at the LHC are presented. The measurements are performed in the pseudorapidity range − 3.5 < η < 5 for various centrality intervals using two- and multi-particle cumulants with the subevent method. The flow probability density function (p.d.f.) is studied with the ratio of flow coefficient v 2 calculated with four- and two-particle cumulant, and suggests that the variance of flow p.d.f. is independent of pseudorapidity. The decorrelation of the flow vector in the longitudinal direction is probed using two-particle correlations. The results measured with respect to different reference regions in pseudorapidity exhibit differences, argued to be a result of saturating decorrelation effect above a certain pseudorapidity separation, in contrast to previous publications which assign this observation to non-flow effects. The results are compared to 3 + 1 dimensional hydrodynamic and the AMPT transport model calculations. Neither of the models is able to simultaneously describe the pseudorapidity dependence of measurements of anisotropic flow and its fluctuations. The results presented in this work highlight shortcomings in our current understanding of initial conditions and subsequent system expansion in the longitudinal direction. Therefore, they provide input for its improvement.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Computational studies of hydraulic stressors for biological performance assessment in a hydropower plant with Kaplan turbine

We report Hydropower is currently one of the preeminent sources of renewable energy in the United States and globally. Hydropower plants also have detrimental impacts on the environment and ecology, including direct impacts to anadromous fish populations. The computational fluid dynamics (CFD) – based Biological Performance Assessment (BioPA) toolset is used for biological evaluations of fish passage through hydropower plants. The hydraulic environment of a hydropower plant was evaluated using CFD coupled with discrete element method (DEM) simulations. The predicted flow field and particle collision rate were validated against the experimental data in a water flume that has an idealized hydroturbine distributor geometry. Flow simulations were conducted to evaluate the hydraulic stressors, such as nadir pressure, fluid shear, runner collision, in a physical scale in a hydropower plant with Kaplan turbine which are responsible for injury and mortality of fish in a downstream migration. The cumulative exposure probability for the nadir pressure and collision with turbine runner was found to decrease with increased discharge rate. The lowest discharge rate shows the higher value of cumulative shear exposure probability. The maximum value of collision velocity increases with increased discharge rate. We offer the conclusions that will help in understanding various hydraulic stressors for biological assessment for environmentally sustainable hydroturbine passage.

13 HYDRO ENERGY↗

Water level changes in Lake Erie drive 21st century CO 2 and CH 4 fluxes from a coastal temperate wetland

Wetland water depth influences microbial and plant communities, which can alter the above- and below-ground carbon cycling of a wetland. Wetland water depths are likely to change due to shifting precipitation patterns, which will affect projections of greenhouse gas emissions; however, these effects are rarely incorporated into wetland greenhouse gas models. Seeking to address this gap, we used a mechanistic model, ecosys, to simulate a range of water depth scenarios in a temperate wetland, and analyzed simulated predictions of carbon dioxide (CO 2 ) andmethane (CH 4 ) fluxes over the 21st century. We tested our model using eddy covariance measurements of CO 2 and CH 4 fluxes collected at the Old Woman Creek National Estuarine Research Reserve (OWC) during 2015 and 2016. OWC is a lacustrine, estuarine, freshwater, mineral-soil marsh. An empirical model found that the wetland water depth is highly dependent on the water depth of the nearby Lake Erie. Future wetland surface water depths were modeled based on projection of Lake Erie's water depth using four separate NOAA projections, resulting in four wetland water-depth scenarios. Two of the four 21st century projections for Lake Erie water depths used in this study indicated that the water depth of the wetland would remain nearly steady; however, the other two indicated decreases in the wetland water depth. In our scenario where the wetland dries out, we project the wetland's climatological warming effect will decrease due to smaller CH 4 fluxes to the atmosphere and larger CO 2 uptake by the wetland. We also found that increased water level can lower emissions by shifting the site towards more open water areas, which have lower CH 4 emissions. We found that decreased water depths would cause more widespread colonization of the wetland by macrophyte vegetation. Using an empirical relationship, we also found that further drying could result in other, non-wetland vegetation to emerge, dramatically altering soil carbon cycling. In three of our four projections, we found that in general the magnitude of CO 2 and CH 4 fluxes steadily increase over the next 100 years in response to higher temperatures. However, in our driest simulations, we projected a different response due to increased oxidation of soil carbon, with CH 4 emissions decreasing substantially from an annual cumulative peak of 224.6 to a minimum of 104.7 gC m –2 year –1 . In that same simulation, net cumulative flux of CO 2 changed from being a sink of 56.5 gC m –2 year –1 to a source of 369.6 gCm –2 year –1 over the same period, despite a temperature increase from 13.7 °C to 14.2 °C. This temperature shift in our other three cases with greater water depths increased the source strength of CH 4 and the sink strength of CO 2 . Here, we conclude that the magnitude of wetland greenhouse-gas fluxes depended on the water depth primarily as it affected the areal percentage of the wetland available for plant colonization, but dramatic decreases in water depths could cause significant reductions in the wetland CH 4 fluxes, while simultaneously altering the wetland vegetation.

54 ENVIRONMENTAL SCIENCES↗

The Stability of Nd Hydroxyl Complexes at Near Neutral to Alkaline pH from 25 to 75 °C: Implications for Monazite Solubility in Hydrothermal Aqueous Fluids

Speciation of rare earth elements (REE) in aqueous fluids controls their mobilization during fluid-rock interaction. Thermodynamic modeling provides important insights into the factors controlling REE mobility and fractionation in aqueous fluids, in which aqueous REE hydroxyl complexes are particularly important across a broad pH range. However, the stability of these REE hydroxyl complexes remains poorly determined due to limited experimental data above ambient temperature. In this work, UV–vis spectrophotometry experiments were conducted to determine the cumulative hydrolysis constants for neodymium (Nd) from 25 to 75 °C in near-neutral to alkaline solutions. The color indicator m -cresol purple was used to determine in situ pH. Alkaline NaOH-bearing aqueous solutions were doped with varying initial Nd concentrations (0 to ∼0.155 mmol/kg), which resulted in the release of protons (H + ) and pH decrease during hydrolysis of Nd 3+ according to Nd 3+ + n H 2 O = Nd(OH) n 3– n + H + . The average OH – ligand number coordinated to Nd 3+ was found to increase from 1.0 to 1.6 at 25 °C and from 1.2 to 2.8 at 75 °C over a pH range from 6.3 to 9.0. The measured speciation shows an increased predominance of Nd(OH) 3 0 over the Nd(OH) 2 + and Nd(OH) 2+ species with increased temperature and pH. The increasing cumulative formation constants (log*β° n , n = 1 to 3) retrieved from 25 to 75 °C differ by 0.1 to 1.9 logarithmic units in comparison to existing thermodynamic databases. These updated thermodynamic data have important implications for geochemical modeling of the speciation of Nd hydroxyl complexes and the solubility of monazite as a function of pH.

58 GEOSCIENCES↗

Linking Hydrology and Dissolved Organic Matter Characteristics in a Subtropical Wetland: A Long-Term Study of the Florida Everglades

Dissolved organic matter (DOM) acts as an important biogeochemical component of aquatic ecosystems that controls nutrient cycling, influences water quality, and links terrestrial and oceanic carbon pools, yet long-term studies of how changing environmental drivers alter its abundance and composition are rare. Using a ten-year, spatially explicit dataset from Everglades National Park (ENP), a globally significant wetland, we investigated the relationships between DOM quality/quantity and hydrologic/climatic drivers along two contrasting marsh-estuarine transects based on generalized linear modeling and a cumulative sums analysis. Analyses revealed distinct spatial, seasonal, and interannual patterns in variability of DOC and optical properties. Landscape-scale seasonal patterns showed an enrichment in microbial-like and protein-like DOM during the dry season relative to the wet season. We found that, while some compositional constituents varied with the solar calendar, responsive to temperature and photoperiod, others varied with the hydrologic calendar. Independent water level and discharge effects indicated strong hydrologic control on DOM quality that differed between the two transects, evidencing differences in their connectivity to areas of high agricultural activity. Across all sites, a significant long-term increasing trend in the fluorescence index was observed, associated with a positive correlation with precipitation and also potential changes in agricultural inputs, with other features associated with drought and hurricanes. Lastly, the cumulative sums analysis revealed differences between the two transects in the sensitivity of DOM composition to decreased water levels associated with 30-year climate scenarios, with the less hydrologically dynamic transect exhibiting greater potential sensitivity.

54 ENVIRONMENTAL SCIENCES↗

Tundra Underlain By Thawing Permafrost Persistently Emits Carbon to the Atmosphere Over 15 Years of Measurements

Warming of the Arctic can stimulate microbial decomposition and release of permafrost soil carbon (C) as greenhouse gases, and thus has the potential to influence climate change. At the same time, plant growth can be stimulated and offset C release. This study presents a 15-year time series comprising chamber and eddy covariance measurements of net ecosystem C exchange in a tundra ecosystem in Alaska where permafrost has been degrading due to regional warming. The site was a carbon dioxide source to the atmosphere with a cumulative total loss of 781.6 g C m -2 over the study period. Both gross primary productivity (GPP) and ecosystem respiration (R eco ) were already likely higher than historical levels such that increases in R eco losses overwhelmed GPP gains in most years. This shift to a net C source to the atmosphere likely started in the early 1990s when permafrost was observed to warm and thaw at the site. Shifts in the plant community occur more slowly and are likely to constrain future GPP increases as compared to more rapid shifts in the microbial community that contribute to increased Reco. Observed rates suggest that cumulative net soil C loss of 4.18–10.00 kg C m -2 —8%–20% of the current active layer soil C pool—could occur from 2020 to the end of the century. This amount of permafrost C loss to the atmosphere represents a significant accelerating feedback to climate change if it were to occur at a similar magnitude across the permafrost region.

54 ENVIRONMENTAL SCIENCES↗

How Does Flow Alteration Propagate Across a Large, Highly Regulated Basin? Dam Attributes, Network Context, and Implications for Biodiversity

Abstract Large dams are a leading cause of river ecosystem degradation. Although dams have cumulative effects as water flows downstream in a river network, most flow alteration research has focused on local impacts of single dams. Here we examined the highly regulated Colorado River Basin (CRB) to understand how flow alteration propagates in river networks, as influenced by the location and characteristics of dams as well as the structure of the river network—including the presence of tributaries. We used a spatial Markov network model informed by 117 upstream‐downstream pairs of monthly flow series (2003–2017) to estimate flow alteration from 84 intermediate‐to‐large dams representing >83% of the total storage in the CRB. Using Least Absolute Shrinkage and Selection Operator regression, we then investigated how flow alteration was influenced by local dam properties (e.g., purpose, storage capacity) and network‐level attributes (e.g., position, upstream cumulative storage). Flow alteration was highly variable across the network, but tended to accumulate downstream and remained high in the main stem. Dam impacts were explained by network‐level attributes (63%) more than by local dam properties (37%), underscoring the need to consider network context when assessing dam impacts. High‐impact dams were often located in sub‐watersheds with high levels of native fish biodiversity, fish imperilment, or species requiring seasonal flows that are no longer present. These three biodiversity dimensions, as well as the amount of dam‐free downstream habitat, indicate potential to restore river ecosystems via controlled flow releases. Our methods are transferrable and could guide screening for dam reoperation in other highly regulated basins.

54 ENVIRONMENTAL SCIENCES↗

Method of Distributions for Two‐Phase Flow in Heterogeneous Porous Media

Abstract Multiscale heterogeneity and insufficient characterization data for a specific subsurface formation of interest render predictions of multi‐phase fluid flow in geologic formations highly uncertain. Quantification of the uncertainty propagation from the geomodel to the fluid‐flow response is typically done within a probabilistic framework. This task is computationally demanding due to, for example, the slow convergence of Monte Carlo simulations (MCS), especially when computing the tails of a distribution that are necessary for risk assessment and decision‐making under uncertainty. The frozen streamlines method (FROST) accelerates probabilistic predictions of immiscible two‐phase fluid flow problems; however, FROST relies on MCS to compute the travel‐time distribution, which is then used to perform the transport (phase saturation) computations. To alleviate this computational bottleneck, we replace MCS with a deterministic equation for the cumulative distribution function (CDF) of travel time. The resulting CDF‐FROST approach yields the CDF of the saturation field without resorting to sampling‐based strategies. Our numerical experiments demonstrate the high accuracy of CDF‐FROST in computing the CDFs of both saturation and travel time. For the same accuracy, it is about 5 and 10 times faster than FROST and MCS, respectively.

Yang, Hyung Jun↗