SEARCH · Engineering Papers

Results for “global optimization”

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.

SourceResearch connectionFrom yearThrough year

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

At least 1,657 records · Page 92

All sources NASA NTRS DOE OSTI

NASA NTRS1994

TOPEX/POSEIDON mission overview

TOPEX/POSEIDON is the first space mission specifically designed and conducted for studying the circulation of the world's oceans. A state-of-the-art radar altimetry system is used to measure the precise height of sea level, from which information on the ocean circulation is obtained. The satellite, launched on August 10, 1992, has been making observations of the global oceans with unprecedented accuracy since late September 1992. To meet the stringent measurement accuracy required for ocean circulation studies, a number of innovative improvements have been made to the mission design, including the first dual-frequency space-borne radar altimeter capable of retrieving the ionospheric delay of the radar signal, a three-frequency microwave radiometer for retrieving the signal delay caused by the water vapor in the troposphere, an optimal model of the Earth's gravity field and multiple satellite tracking systems for precision orbit determination. Additionally, the satellite also carries two experimental instruments to demonstrate new technologies: a single-frequency solid-state altimeter for the technology of low-power, low-weight altimeter and a Global Positioning System receiver for continuous,precise satellite tracking. The performance of the mission's measurement system has been tested by numerous verification studies. The results indicate that the root-sum-square accuracy of a single-pass sea level measurement is 4.7 cm for the TOPEX system and 5.1 cm for the POSEIDON system; both are more than a factor of 2 better than the requirement of 13.7 cm. This global data set is being analyzed to improve understanding of the global ocean circulation as well as the ocean tides, geodesy, and geodynamics, and ocean wind and waves. The mission is designed to last for at least 3 years with a possible extension to 6 years. The multiyear global data set will go a long way toward understanding the ocean circulation and its variability in relation to climate change. A summary of the mission's systems and their performance as well as the mission's science team is presented.

Fu, Lee-Lueng↗

DOE OSTI2016

The CU 2-D-MAX-DOAS instrument – Part 2: Raman scattering probability measurements and retrieval of aerosol optical properties

Abstract. The multiannual global mean of aerosol optical depth at 550 nm (AOD550) over land is ∼ 0.19, and that over oceans is ∼ 0.13. About 45 % of the Earth surface shows AOD550 smaller than 0.1. There is a need for measurement techniques that are optimized to measure aerosol optical properties under low AOD conditions. We present an inherently calibrated retrieval (i.e., no need for radiance calibration) to simultaneously measure AOD and the aerosol phase function parameter, g, based on measurements of azimuth distributions of the Raman scattering probability (RSP), the near-absolute rotational Raman scattering (RRS) intensity. We employ radiative transfer model simulations to show that for solar azimuth RSP measurements at solar elevation and solar zenith angle (SZA) smaller than 80°, RSP is insensitive to the vertical distribution of aerosols and maximally sensitive to changes in AOD and g under near-molecular scattering conditions. The University of Colorado two-dimensional Multi-AXis Differential Optical Absorption Spectroscopy (CU 2-D-MAX-DOAS) instrument was deployed as part of the Two Column Aerosol Project (TCAP) at Cape Cod, MA, during the summer of 2012 to measure direct sun spectra and RSP from scattered light spectra at solar relative azimuth angles (SRAAs) between 5 and 170°. During two case study days with (1) high aerosol load (17 July, 0.3 < AOD430 < 0.6) and (2) near-molecular scattering conditions (22 July, AOD430 < 0.13) we compare RSP-based retrievals of AOD430 and g with data from a co-located CIMEL sun photometer, Multi-Filter Rotating Shadowband Radiometer (MFRSR), and an airborne High Spectral Resolution Lidar (HSRL-2). The average difference (relative to DOAS) for AOD430 is +0.012 ± 0.023 (CIMEL), −0.012 ± 0.024 (MFRSR), −0.011 ± 0.014 (HSRL-2), and +0.023 ± 0.013 (CIMELAOD − MFRSRAOD) and yields the following expressions for correlations between different instruments: DOASAOD = −(0.019 ± 0.006) + (1.03 ± 0.02) × CIMELAOD (R2 = 0.98), DOASAOD = −(0.006 ± 0.005) + (1.08 ± 0.02) × MFRSRAOD (R2 = 0.98), and CIMELAOD = (0.013 ± 0.004) + (1.05 ± 0.01) × MFRSRAOD (R2 = 0.99). The average g measured by DOAS on both days was 0.66 ± 0.03, with a difference of 0.014 ± 0.05 compared to CIMEL. Active steps to minimize the error in the RSP help to reduce the uncertainty in retrievals of AOD and g. As AOD decreases and SZA increases, the RSP signal-to-noise ratio increases. At AOD430 ∼ 0.4 and 0.10 the absolute AOD errors are ∼ 0.014 and 0.003 at 70° SZA and 0.02 and 0.004 at 35° SZA. Inherently calibrated, precise AOD and g measurements are useful to better characterize the aerosol direct effect in urban polluted and remote pristine environments.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Valence Disproportionation of GeS in the PbS Matrix Forms Pb 5 Ge 5 S 12 Inclusions with Conduction Band Alignment Leading to High n-Type Thermoelectric Performance

Converting waste heat into useful electricity using solid-state thermoelectrics has a potential for enormous global energy savings. Lead chalcogenides are among the most prominent thermoelectric materials, whose performance decreases with an increase in chalcogen amounts (e.g., PbTe>PbSe>PbS). Herein, we demonstrate the simultaneous optimization of the electrical and thermal transport properties of PbS-based compounds by alloying with GeS. The addition of GeS triggers a complex cascade of beneficial events as follows: Ge 2+ substitution in Pb 2+ and discordant off-center behavior; formation of Pb 5 Ge 5 S 12 as stable second phase inclusions through valence disproportionation of Ge 2+ to Ge0 and Ge 4+ . PbS and Pb 5 Ge 5 S 12 exhibit good conduction band energy alignment that preserves the high electron mobility; the formation of Pb 5 Ge 5 S 12 increases the electron carrier concentration by introducing S vacancies. Sb doping as the electron donor produces a large power factor and low lattice thermal conductivity (κ lat ) of ~0.61 Wm -1 K -1 . The highest performance was obtained for the 14% GeS-alloyed samples, which exhibited an increased room temperature electron mobility of ~121 cm 2 V -1 s -1 for 3 × 10 19 cm -3 carrier density, and a ZT, of 1.32 at 923 K. This is ~ 55% greater that the corresponding Sb-doped PbS sample and is one of the highest reported for the n-type PbS system. Moreover, the average ZT (ZT avg ) of ~0.76 from 400 to 923 K is the highest for PbS-based systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2022

Constraining the nuclear gluon PDF with inclusive hadron production data

The nuclear parton distribution functions (nPDFs) of gluons are known to be difficult to determine with fits of deep inelastic scattering (DIS) and Drell-Yan (DY) data alone. Therefore, the nCTEQ15 analysis of nuclear PDFs added inclusive neutral pion production data from RHIC to help in constraining the gluon. In this analysis, we present a new global analysis of nuclear PDFs based on a much larger set of single inclusive light hadron data from RHIC and the LHC. Using our new nCTEQ code (nCTEQ++) with an optimized version of INCNLO we study systematically the limitations of the theory and the impact of the fragmentation function uncertainty.

Duwentäster, Pit↗

DOE OSTI2022

Co-Optimization of Fuels & Engines Findings and Impact (FY2015-FY2021 Report)

Breakthroughs from the U.S. Department of Energy Co-Optimization of Fuels & Engines (Co-Optima) initiative could make it possible to more rapidly cut emissions, reduce dependence on international petroleum, and contribute to ambitious national goals to slow global warming on the land, in the air, and across the water. Findings from the 6-year collaborative undertaking are now available in this report.

33 ADVANCED PROPULSION SYSTEMS↗

NASA NTRS1987

Numerical analysis of discrete fractional integrodifferential structural dampers

This paper develops solution algorithms enabling the handling of the dynamic response of nonlinear structures contained discretely attached dampers modeled by fractional integrodifferential operators of the Grunwald-Liouville-Riemann type. The development consists of two levels of formulation, namely: (1) numerical approximations of fractional operators and, (2) the establishment of global level implicit schemes enabling the solution to nonlinear structural formulations. To generalize the overall results, error estimates are derived for the fractional operator approximation algorithm. These enable an ongoing optimization of solution efficiency for a given error tolerance. To benchmark the scheme, the results of several numerical experiments are presented. These illustrate the numerical characteristics of the overall formulation.

Padovan, Joseph↗

NASA NTRS1992

Structural design optimization with survivability dependent constraints application: Primary wing box of a multi-role fighter

The superior survivability of a multirole fighter is dependent upon balanced integration of technologies for reduced vulnerability and susceptability. The objective is to develop a methodology for structural design optimization with survivability dependent constraints. The design criteria for optimization will be survivability in a tactical laser environment. The following analyses are studied to establish a dependent design relationship between structural weight and survivability: (1) develop a physically linked global design model of survivability variables; and (2) apply conventional constraints to quantify survivability dependent design. It was not possible to develop an exact approach which would include all aspects of survivability dependent design, therefore guidelines are offered for solving similar problems.

Dolvin, Douglas J.↗

NASA NTRS2005

Evolution of Galaxies and the Star Formation Rate in the Infrared

A central goal of extragalactic observational astronomy is to understand how normal galaxies evolve with redshift, and particularly when galaxies formed their stars. While optical and rest-frame UV observations have begun to address these issues, the interpretation of such data is particularly challenging because of the sensitivity to dust obscuration (at optical and UV wavelengths). The absorbed light is re-radiated at IR wavelengths, hence the optimal indicators of the star formation rate (SFR) is at a rest-frame wavelength of approx. 60 microns. The Spitzer Space Telescope mission is revolutionizing the study of the global properties and evolution of galaxies. Spitzer reaches nearly two orders of magnitude more sensitivity than previous IR space missions. This research program is to study the SFR using statistical samples of galaxies in the local universe, at intermediate redshifts, and set the stage for continuing studies up to z=5. The overall research program is divided into three main investigations: A Mid-IR Hubble Atlas and SFR estimators in the local universe, Evolution of the SFR at 0 < z < 1 using pencil beam redshift surveys, and Galaxy formation and evolution at 1 < z < 5. The first papers from Spitzer were published during the last year, including ten refereed journal papers where the PI was first or co-author.

Pahre, Michael A.↗

NASA NTRS2005

Constraining Oceanic Sources and Sinks of CO2 by Inverse Modeling: A First Step Toward Assimilation of In Situ and Remote Data

The primary goal of this collaborative project, lead by Nicholas Gruber at UCLA, is to better quantify the air-sea fluxes of pre-industrial and anthropogenic carbon dioxide CO2 for the global ocean with a newly developed inverse method. The method involves using ocean general circulation models (OGCMs) to characterize the subsurface distribution of tracers injected for different regions of the ocean surface. By linearly combining a series of these basis (or Greens) functions, one should be able to recreate any arbitrary ocean tracer field. In our case, we are using the global database of ocean dissolved inorganic carbon (DIC) from the WOCE/JGOFS global CO2 survey from the 1990s. The observed DIC field has been partitioned into pre-industrial and anthropogenic components using the empirical C* technique. We then solve for an optimal combination of the regional basis functions that best match the two sets observations in a least-squares sense with singular value decomposition. This in turn provides solutions for the regional air-sea CO2 fluxes (pre-industrial and anthropogenic) that are independent of more traditional estimates derived from surface water disequilibria in the partial pressure of CO2 and wind speed gas exchange parameterizations. We are apply the method to a suite of nine different OGCMs to test the robustness of the derived surface fluxes fields to differences in ocean circulation. Together with a collegue Keith Lindsay at NCAR, S. Doney's role in this project was to create the appropriate steady-state and temporally evolving basis functions from the NCAR CSM-1 ocean circulation model. Because of a slow start across the project in defining a common set of experimental protocols, we requested and were given a one year, no cost extension. Both the steady-state and temporally evolving basis functions were delivered to the UCLA group for further analysis. Additional information is included in the original extended abstract.

Doney, Scott↗

DOE OSTI2021

Global study of separable pairing interaction in covariant density functional theory

A systematic global investigation of pairing properties based on all available experimental data on pairing indicators has been performed for the first time in the framework of covariant density functional theory. It is based on the separable pairing interaction. The optimization of the scaling factors of this interaction to experimental data clearly reveals its isospin dependence in the neutron subsystem. However, the situation is less certain in the proton subsystem since similar accuracy of the description of pairing indicators can be achieved both with isospin-dependent and mass-dependent scaling factors. The differences in the functional dependencies of scaling factors lead to uncertainties in the prediction of proton and neutron pairing properties which are especially pronounced at high isospin and could have a significant impact on some physical observables. For a given part of the nuclear chart the scaling factors for spherical nuclei are smaller than those for deformed ones; this feature exists also in nonrelativistic density functional theories. Furthermore, its origin is traced back to particle-vibration coupling in odd-A nuclei which is missing in all existing global studies of pairing. Although the present investigation is based on the NL5(E) covariant energy density functional (CEDF), its general conclusions are expected to be valid also for other CEDFs built at the Hartree level.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2014

Towards local progression estimation of pulmonary emphysema using CT

Purpose: Whole lung densitometry on chest CT images is an accepted method for measuring tissue destruction in patients with pulmonary emphysema in clinical trials. Progression measurement is required for evaluation of change in health condition and the effect of drug treatment. Information about the location of emphysema progression within the lung may be important for the correct interpretation of drug efficacy, or for determining a treatment plan. The purpose of this study is therefore to develop and validate methods that enable the local measurement of lung density changes, which requires proper modeling of the effect of respiration on density. Methods: Four methods, all based on registration of baseline and follow-up chest CT scans, are compared. The first naïve method subtracts registered images. The second employs the so-called dry sponge model, where volume correction is performed using the determinant of the Jacobian of the transformation. The third and the fourth introduce a novel adaptation of the dry sponge model that circumvents its constant-mass assumption, which is shown to be invalid. The latter two methods require a third CT scan at a different inspiration level to estimate the patient-specific density-volume slope, where one method employs a global and the other a local slope. The methods were validated on CT scans of a phantom mimicking the lung, where mass and volume could be controlled. In addition, validation was performed on data of 21 patients with pulmonary emphysema. Results: The image registration method was optimized leaving a registration error below half the slice increment (median 1.0 mm). The phantom study showed that the locally adapted slope model most accurately measured local progression. The systematic error in estimating progression, as measured on the phantom data, was below 2 gr/l for a 70 ml (6%) volume difference, and 5 gr/l for a 210 ml (19%) difference, if volume correction was applied. On the patient data an underlying linearity assumption relating lung volume change with density change was shown to hold (fitR{sup 2} = 0.94), and globalized versions of the local models are consistent with global results (R{sup 2} of 0.865 and 0.882 for the two adapted slope models, respectively). Conclusions: In conclusion, image matching and subsequent analysis of differences according to the proposed lung models (i) has good local registration accuracy on patient data, (ii) effectively eliminates a dependency on inspiration level at acquisition time, (iii) accurately predicts progression in phantom data, and (iv) is reasonably consistent with global results in patient data. It is therefore a potential future tool for assessing local emphysema progression in drug evaluation trials and in clinical practice.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2023

Controlling circuitry underlies the growth optimization of Saccharomyces cerevisiae

Microbial growth emerges from coordinated synthesis of various cellular components from limited resources. In Saccharomyces cerevisiae, cyclic AMP (cAMP)-mediated signaling is shown to orchestrate cellular metabolism; however, it remains unclear quantitatively how the controlling circuit drives resource partition and subsequently shapes biomass growth. Here we combined experiment with mathematical modeling to dissect the signaling-mediated growth optimization of S. cerevisiae. We showed that, through cAMP-mediated control, the organism achieves maximal or nearly maximal steady-state growth during the utilization of multiple tested substrates as well as under perturbations impairing glucose uptake. However, the optimal cAMP concentration varies across cases, suggesting that different modes of resource allocation are adopted for varied conditions. Under settings with nutrient alterations, S. cerevisiae tunes its cAMP level to dynamically reprogram itself to realize rapid adaptation. Moreover, to achieve growth maximization, cells employ additional regulatory systems such as the GCN2-mediated amino acid control. This study establishes a systematic understanding of global resource allocation in S. cerevisiae, providing insights into quantitative yeast physiology as well as metabolic strain engineering for biotechnological applications.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2025

Multi-plane moment-of-fluid interface reconstruction in 3D

Moment-of-fluid (MOF) methods for interface reconstruction approximate the region occupied by material in each mesh element only through reference to its geometric moments. Here, we present a 3D MOF method that represents the material (POM) in each cell as the convex intersection of the cell and multiple half-spaces, each selected to minimize the least-squares error between computed moments of the approximated material and provided reference moments. This optimization problem is highly non-linear and non-convex, making the numerical result very sensitive to the initial guess. To create an effective initial guess in each cell, we construct an ellipsoid from 0th–2nd order reference moments such that its shape corresponds with that of the POM. Within this ellipsoid we inscribe a polyhedron, and initialize the minimization problem with the half-spaces defined by each of its faces. The inscribed polyhedron has minimally 4 faces, and using up to 3rd order moments permits optimization over up to 20 unknown values. We therefore define MOF methods that utilize 4, 5, or 6 half-spaces, correspondingly initialized with the faces of a single inscribed tetrahedron, triangular prism, or hexahedron. Stability of the non-linear optimization is further improved with a prepossessing step that normalizes the reference moments according to the axes of the reference ellipsoid. Using this approach, the non-linear least-squares solver reliably converges to a near-global minimum from a single initial guess. We demonstrate accuracy and robustness using single-cell and multi-cell examples over a wide spectrum of geometry. In particular, we demonstrate our ability to exactly reproduce several important and complex features defined by up to four half-spaces, such as corners, filaments, filament tips, and embedded material in the cell.

3D interface reconstruction↗

NASA NTRS1999

Systematic Differences between Satellite-Based Presipitation Climatologies over the Tropical Oceans

Since the beginning of the World Climate Research Program's Global Precipitation Climatology Project (GPCP) satellite remote sensing of precipitation has made dramatic improvements, particularly for tropical regions. Data from microwave and infrared sensors now form the most critical input to precipitation data sets and can be calibrated with surface gauges to so that the strengths of each data source can be maximized in some statistically optimal sense. It is clear however that there still remain significant uncertainties with satellite precipitation retrievals which limit their usefulness for many purposes. Systematic differences i'A tropical precipitation estimates have been brought to light in comparison activities such as the GPCP Algorithm Intercomparison Project and more recent Wetnet Precipitation Intercomparison Project 3. These uncertainties are assuming more importance because of the demands for validation associated with global climate modeling and data assimilation methodologies. The objective of the present study is to determine the physical basis for systematic differences in spatial structure of tropical precipitation as portrayed by several different satellite-based data sets. The study is limited to oceanic regions only and deals primarily with aspects of spatial variability. We are specifically interested in why MSU channel 1 and GPI precipitation differences are so striking over the Eastern Pacific ITCZ and why they both differ from other microwave emission-based precipitation estimates from SSM/I and a scattering-based deep convective ice index from MSU channel 2. Our results to date have shown that MSU channel I precipitation estimates are biased high over the Eastern Pacific ITCZ because of two factors: (1) the hypersensitivity of this frequency to cloud water in contrast to falling rain drops, and (2) unaccounted for scattering effects by precipitation-size ice which depresses the signal of the liquid water emission. Likewise, cold cloud top climatologies such as the GPI show an excess (a deficit) in estimated rainfall over the E. Pacific ITCZ (Warm Pool region). We show that these algorithms need to account for regionally varying heights (or temperatures) at which tropical convection detrains to form cirrus shields. A second objective we pursue is to identify variations in the macroscale cloud physical and thermodynamic properties of precipitation regimes" and relate these differences to tropical dynamical mechanisms of tropical heat and moisture balance. Finally, we interpret the algorithm differences and their associations with tropical dynamics in terms of WCRP GPCP goals for constructing precipitation climatologies.

Robertson, Frankin R.↗

DOE OSTI2019

Coordinating Multiple Sources for Service Restoration to Enhance Resilience of Distribution Systems

When a major outage occurs on a distribution system due to extreme events, microgrids, distributed generators, storage devices, and other local resources can be used to restore critical loads and enhance resiliency. This paper proposes a decision-making method to determine the optimal restoration strategy coordinating multiple sources to serve critical loads after a major outage. The critical load restoration problem is formulated as a mixed -integer semidefinite program. The objective is maximizing the number of loads restored, weighted by their priority. The unbalanced three-phase power flow constraint and operational constraints are considered. An iterative algorithm is proposed to deal with integer variables and can attain the global optimum of the critical load restoration problem by solving a few semidefinite programs under mild conditions. The effectiveness of the proposed models and algorithm is validated by numerical simulation with the modified IEEE 13-node test feeder and the modified IEEE123-node test feeder under a large number of scenarios. The results indicate that the optimal restoration strategy can be determined efficiently in most scenarios

Wang, Ying↗

DOE OSTI2014

Low energy isomers of (H{sub 2}O){sub 25} from a hierarchical method based on Monte Carlo temperature basin paving and molecular tailoring approaches benchmarked by MP2 calculations

We report new global minimum candidate structures for the (H{sub 2}O){sub 25} cluster that are lower in energy than the ones reported previously and correspond to hydrogen bonded networks with 42 hydrogen bonds and an interior, fully coordinated water molecule. These were obtained as a result of a hierarchical approach based on initial Monte Carlo Temperature Basin Paving sampling of the cluster's Potential Energy Surface with the Effective Fragment Potential, subsequent geometry optimization using the Molecular Tailoring Approach with the fragments treated at the second order Møller-Plesset (MP2) perturbation (MTA-MP2) and final refinement of the entire cluster at the MP2 level of theory. The MTA-MP2 optimized cluster geometries, constructed from the fragments, were found to be within <0.5 kcal/mol from the minimum geometries obtained from the MP2 optimization of the entire (H{sub 2}O){sub 25} cluster. In addition, the grafting of the MTA-MP2 energies yields electronic energies that are within <0.3 kcal/mol from the MP2 energies of the entire cluster while preserving their energy rank order. Finally, the MTA-MP2 approach was found to reproduce the MP2 harmonic vibrational frequencies, constructed from the fragments, quite accurately when compared to the MP2 ones of the entire cluster in both the HOH bending and the OH stretching regions of the spectra.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Incorporating the Acclimation of Photosynthesis and Leaf Respiration in the Noah‐MP Land Surface Model: Model Development and Evaluation

Abstract Realistic simulation of leaf photosynthetic and respiratory processes is needed for accurate prediction of the global carbon cycle. These two processes systematically acclimate to long‐term environmental changes by adjusting photosynthetic and respiratory traits (e.g., the maximum photosynthetic capacity at 25°C ( V cmax,25 ) and the leaf respiration rate at 25°C ( R 25 )) following increasingly well‐understood principles. While some land surface models (LSMs) now account for thermal acclimation, they do so by assigning empirical parameterizations for individual plant functional types (PFTs). Here, we have implemented an Eco‐Evolutionary Optimality (EEO)‐based scheme to represent the universal acclimation of photosynthesis and leaf respiration to multiple environmental effects, and that therefore requires no PFT‐specific parameterizations, in a standard version of the widely used LSM, Noah MP. We evaluated model performance with plant trait data from a 5‐year experiment and extensive global field measurements, and carbon flux measurements from FLUXNET2015. We show that observed R 25 and V cmax,25 vary substantially both temporally and spatially within the same PFT ( C.V. >20%). Our EEO‐based scheme captures 62% of the temporal and 70% of the spatial variations in V cmax,25 (73% and 54% of the variations in R 25 ). The standard scheme underestimates gross primary production by 10% versus 2% for the EEO‐based scheme and generates a larger spread in r (correlation coefficient) across flux sites (0.79 ± 0.16 vs. 0.84 ± 0.1, mean ± S.D.). The standard scheme greatly overestimates canopy respiration (bias: ∼200% vs. 8% for the EEO scheme), resulting in less CO 2 uptake by terrestrial ecosystems. Our approach thus simulates climate‐carbon coupling more realistically, with fewer parameters.

Ren, Yanghang [Department of Earth System Science ↗

DOE OSTI2020

Methane and nitrous oxide porewater concentrations and surface fluxes of a regulated river

Greenhouse gas (GHG) emissions from rivers are a critical missing component of current global GHG models. Their exclusion is mainly due to a lack of in-situ measurements and a poor understanding of the spatiotemporal dynamics of GHG production and emissions, which prevents optimal model parametrization. We combined simultaneous observations of porewater concentrations along different beach positions and depths, and surface fluxes of methane and nitrous oxide at a plot scale in a large regulated river during three water stages: rising, falling, and low. Our goal was to gain insights into the interactions between hydrological exchanges and GHG emissions and elucidate possible hypotheses that could guide future research on the mechanisms of GHG production, consumption, and transport in the hyporheic zone (HZ). Results indicate that the site functioned as a net source of methane. Surface fluxes of methane during river water stages at three beach positions (shallow intermediate and deep) correlated with porewater concentrations of methane. However, fluxes were significantly higher in the intermediate position during low water stage, suggesting that low residence time increased methane emissions. Vertical profiles of methane peaked at different depths, indicating an influence of the magnitude and direction of the hyporheic mixing during the different river water stages on methane production and consumption. The site acted as either a sink or a source of nitrous oxide depending on the elevation of the water column. Nitrous oxide porewater concentrations peaked at the upper layers of the sediment throughout the different water stages. River hydrological stages significantly influenced porewater concentrations and fluxes of GHG, probably by influencing heterotrophic respiration (production and consumption processes) and transport to and from the HZ. Our results highlight the importance of including dynamic hydrological exchanges when studying and modeling GHG production and consumption in the HZ of large rivers.

54 ENVIRONMENTAL SCIENCES↗