Search NASA⌕ Search

SEARCH · Search NASA

Results for “Covariance Analysis”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

DESI DR1 Ly α 1D power spectrum: Validation of estimators

The Data Release 1 (DR1) of the Dark Energy Spectroscopic Instrument (DESI) is the largest sample to date for small-scale Lyα forest cosmology, accessed through its one-dimensional power spectrum (P 1D ). The Lyα forest P 1D is extracted from quasar spectra that are highly inhomogeneous (both in wavelength and between quasars) in noise properties due to intrinsic properties of the quasar, atmospheric and astrophysical contamination, and also sensitive to low-level details of the spectral extraction pipeline. We employ two estimators in DR1 analysis to measure P 1D : the optimal estimator and the fast Fourier transform (FFT) estimator. To ensure robustness of our DR1 measurements, we validate these two power spectrum and covariance matrix estimation methodologies against the challenging aspects of the data. First, using a set of 20 synthetic 1D realizations of DR1, we derive the masking bias corrections needed for the FFT estimator and the continuum fitting bias needed for both estimators. We demonstrate that both estimators, including their covariances, are unbiased with these corrections using the Kolmogorov-Smirnov test. Second, we substantially extend our previous suite of CCD image simulations to include 675,000 quasars, allowing us to accurately quantify the pipeline's performance. This set of simulations reveals biases at the highest k values, corresponding to a resolution error of a few percent. We base the resolution systematics error budget of DR1 P 1D on these values, but do not derive corrections from them since the simulation fidelity is insufficient for precise corrections.

Lyman alpha forest↗
DOE OSTI2020

The Modulation of Gulf Stream Influence on the Troposphere by the Eddy-Driven Jet

This study suggests that the Gulf Stream influence on the wintertime North Atlantic troposphere is most pronounced when the eddy-driven jet (EDJ) is farthest south and better collocated with the Gulf Stream. Using the reanalysis dataset NCEP-CFSR for December–February 1979–2009, the daily EDJ latitude is separated into three regimes (northern, central, and southern). It is found that the average trajectory of atmospheric fronts covaries with EDJ latitude. In the southern EDJ regime (~19% of the time), the frequency of near-surface atmospheric fronts that pass across the Gulf Stream is maximized. Analysis suggests that this leads to significant strengthening in near-surface atmospheric frontal convergence resulting from strong air–sea sensible heat flux gradients (due to strong temperature gradients in the atmosphere and ocean). In recent studies, it was shown that the pronounced band of time-mean near-surface wind convergence across the Gulf Stream is set by atmospheric fronts. Here, it is shown that an even smaller subset of atmospheric fronts—those associated with a southern EDJ—primarily sets the time mean, due to enhanced Gulf Stream air–sea interaction. Furthermore, statistically significant anomalies in vertical velocity extending well above the boundary layer are identified in association with changes in EDJ latitude. These anomalies are particularly strong for a southern EDJ and are spatially consistent with increases in near-surface atmospheric frontal convergence over the Gulf Stream. These results imply that much of the Gulf Stream influence on the time-mean atmosphere is modulated on synoptic time scales, and enhanced when the EDJ is farthest south.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Sex Differences in Odds of Brain Metastasis and Outcomes by Brain Metastasis Status after Advanced Melanoma Diagnosis

Sex differences in cancer are well-established. However, less is known about sex differences in diagnosis of brain metastasis and outcomes among patients with advanced melanoma. Using a United States nationwide electronic health record-derived de-identified database, we evaluated patients diagnosed with advanced melanoma from 1 January 2011–30 July 2022 who received an oncologist-defined rule-based first line of therapy (n = 7969, 33% female according to EHR, 35% w/documentation of brain metastases). The odds of documented brain metastasis diagnosis were calculated using multivariable logistic regression adjusted for age, practice type, diagnosis period (pre/post-2017), ECOG performance status, anatomic site of melanoma, group stage, documentation of non-brain metastases prior to first-line of treatment, and BRAF positive status. Real-world overall survival (rwOS) and progression-free survival (rwPFS) starting from first-line initiation were assessed by sex, accounting for brain metastasis diagnosis as a time-varying covariate using the Cox proportional hazards model, with the same adjustments as the logistic model, excluding group stage, while also adjusting for race, socioeconomic status, and insurance status. Adjusted analysis revealed males with advanced melanoma were 22% more likely to receive a brain metastasis diagnosis compared to females (adjusted odds ratio [aOR]: 1.22, 95% confidence interval [CI]: 1.09, 1.36). Males with brain metastases had worse rwOS (aHR: 1.15, 95% CI: 1.04, 1.28) but not worse rwPFS (adjusted hazard ratio [aHR]: 1.04, 95% CI: 0.95, 1.14) following first-line treatment initiation. Among patients with advanced melanoma who were not diagnosed with brain metastases, survival was not different by sex (rwOS aHR: 1.06 [95% CI: 0.97, 1.16], rwPFS aHR: 1.02 [95% CI: 0.94, 1.1]). This study showed that males had greater odds of brain metastasis and, among those with brain metastasis, poorer rwOS compared to females, while there were no sex differences in clinical outcomes for those with advanced melanoma without brain metastasis.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Within and Among Fish Species Differences in Simulated Turbine Blade Strike Mortality: Limits on the Use of Surrogacy for Untested Species

Use of surrogacy remains a useful method for prioritizing research on representatives of at-risk groups of fishes, yet quantifiable evidence in support of its use is generally not available. Blade strike impact represents one of the most traumatic stressors experienced by fish during non-volitional movements through hydropower turbines. Here, we use data generated from laboratory trials on blade strike impact experiments to directly test use of surrogacy for salmonid and clupeid fishes. Results of logistic regression indicated that a -taxonomic (genus) variable was not a significant predictor of mortality among large rainbow trout and brook trout. Similar results were found for young-of-the-year shad species, but genus-level taxonomy was a significant predictor of mortality while species was not. Multivariate analysis of morphometric data showed that shad clustered together based on similarities in fish shape which was also closely associated with genus. Logistic regression including size as a major covariate suggested total fish length was not a significant predictor of mortality, yet dose–response data suggest differential susceptibility to lower strike velocities. We suggest that use of surrogacy among species is justifiable but should be avoided within a species since the effects of size remain unclear.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Propagating neutronic uncertainties for FFTF LOFWOS Test #13

The safety evaluation conducted for licensing a Sodium-cooled Fast Reactor (SFR) may require modeling transients with best-estimate calculation tools that must first be validated against real-world measurements. To provide the community with a valuable benchmarking opportunity for validating SFR analysis tools and methods, the International Atomic Energy Agency (IAEA) initiated a coordinated research project (CRP) in 2018 for the analysis of the Fast Flux Test Facility (FFTF) Loss of Flow Without Scram (LOFWOS) Test #13. The impact of nuclear data uncertainties on neutronics parameters was previously investigated based on the COMMARA-2.0 covariance matrix. Since the transient simulation results are very sensitive to certain reactivity coefficients, it was decided to employ rigorous uncertainty propagation methods to quantify the impact of nuclear data uncertainties on the best-estimate predication of FFTF LOFWOS Test #13. The DAKOTA code is used to propagate neutronic uncertainties through SAS4A/SASSYS-1 transient simulations, while taking into account spatial and reaction-wise correlations within these uncertainties. This study shows that the remaining discrepancies observed between the Argonne National Laboratory (ANL) best-estimate results and the experimental measurements can be partly explained by the uncertainty associated with Gas Expansion Module (GEM) worth, which contributes the majority of the overall nuclear data uncertainty on the output from the FFTF LOFWOS Test #13 transient simulation. This study also confirmed the importance of including spatial and reaction-wise correlations of nuclear data uncertainties on feedback coefficients in the uncertainty propagation to avoid under-estimating their impact during the transient simulations.

42 ENGINEERING↗
DOE OSTI2020

Online Control of Process Variance Using Feedback

For many stochastic control problem formulations, design objectives typically include constraining the covariance of the outputs of the system. Recent engineering challenges have demonstrated an interest in controlling system output covariances to a setpoint. However, open-loop control of system output covariance can be exceptionally challenging when the system has significant plant uncertainties. This paper explores closed-loop control of the system’s output covariance using a windowed variance function. The stochastic nature of the windowed variance feedback means that root locus and other classical control design techniques cannot be immediately used. Expected value analysis is used to generate transfer functions for the average response of such systems, and the results are used to design controllers. Stochastic simulations are presented to verify the analytical results and make other observations concerning closed-loop variance control.

Bieniek, Mitchell↗
DOE OSTI2021

Spectral Denoising for Accelerated Analysis of Correlated Ionic Transport

Computation of correlated ionic transport properties from molecular dynamics in the Green-Kubo formalism is expensive, as one cannot rely on the affordable mean square displacement approach. We use spectral decomposition of the short-time ionic displacement covariance to learn a set of diffusion eigenmodes that encode the correlation structure and form a basis for analyzing the ionic trajectories. This allows systematic reduction of the uncertainty and accelerate computations of ionic conductivity in systems with a steady-state correlation structure. We provide mathematical and numerical proofs of the method’s robustness and demonstrate it on realistic electrolyte materials.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Nuclear data uncertainty propagation applied to the versatile test reactor conceptual design

We report the Versatile Test Reactor (VTR) currently under development is a 300 MWth sodium-cooled fast reactor (SFR) fueled with ternary metal alloy fuel, which aims to accelerate the testing of advanced nuclear fuels, materials, instrumentation, and sensors in high flux environments that are necessary to license the next generation of advanced reactor concepts. To support the VTR design process, uncertainties associated with the nuclear data has been propagated through the reactor core neutronics calculation to global parameters of interest, such as the core multiplication factor, kinetic parameters, and various reactivity feedback coefficients, following the sensitivity based uncertainty propagation approach. By folding the sensitivity coefficients, separately computed by the generalized perturbation theory code PERSENT and Monte Carlo code Serpent 2, with the variance-covariance matrices from COMMARA-2.0, we obtain the reaction-wise, isotope-wise, and overall uncertainties for each response of interest due to nuclear data uncertainty. With Serpent 2, the statistical error of the uncertainty is obtained by propagating the statistical error of the sensitivity coefficients through the same process using a newly developed uncertainty propagation method. From both codes, the overall top uncertainty contributors are found to be the cross section of Fe-56 elastic scattering, Na-23 elastic scattering, and U 238 inelastic scattering. The large contributions of the Fe-56 elastic scattering cross sections to global parameters are due to its relatively large relative uncertainty of 5–10% in nuclear data and the large volume of Fe-containing reflector assemblies in the fairly compact VTR core design. Both codes agreed well for the overall uncertainty estimates of all responses of interest, except the delayed neutron fraction, prompt neutron generation time, and the coolant density feedback coefficient, where Serpent 2 yielded a much larger value than PERSENT due to the large statistical error of sensitivity coefficients. The calculated uncertainties are also compared to those associated with other SFR cores. Another outcome of this study is a variance-covariance matrix of reactivity coefficients, which can be used in the subsequent uncertainty propagation to the system level to investigate the impact of identified uncertainties on system responses in the safety analysis.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Understanding Hailstone Temporal Variability and Contributing Factors over the U.S. Southern Great Plains

Abstract Hailstones are a natural hazard that pose a significant threat to property and are responsible for significant economic losses each year in the United States. Detailed understanding of their characteristics is essential to mitigate their impact. Identifying the dynamic and physical factors contributing to hail formation and hailstone sizes is of great importance to weather and climate prediction and policymakers. In this study, we have analyzed the temporal and spatial variabilities of severe hail occurrences over the U.S. southern Great Plains (SGP) states from 2004 to 2016 using two hail datasets: hail reports from the Storm Prediction Center and the newly developed radar-retrieved maximum expected size of hail (MESH). It is found that severe and significant severe hail occurrences have a considerable year-to-year temporal variability in the SGP region. The interannual variabilities have a strong correspondence with sea surface temperature anomalies over the northern Gulf of Mexico and there is no outlier. The year 2016 is identified as an outlier for the correlations with both El Niño–Southern Oscillation (ENSO) and aerosol loading. The correlations with ENSO and aerosol loading are not statistically robust to inclusion of the outlier 2016. Statistical analysis without the outlier 2016 shows that 1) aerosols that may be mainly from northern Mexico have the largest correlation with hail interannual variability among the three factors and 2) meteorological covariation does not significantly contribute to the high correlation. These analyses warrant further investigations of aerosol impacts on hail occurrence.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

On the statistical theory of self-gravitating collisionless dark matter flow

Dark matter, if it exists, accounts for five times as much as the ordinary baryonic matter. Compared to hydrodynamic turbulence, the flow of dark matter might possess the widest presence in our universe. This paper presents a statistical theory for the flow of dark matter that is compared with N-body simulations. By contrast to hydrodynamics of normal fluids, dark matter flow is self-gravitating, long-range, and collisionless with a scale-dependent flow behavior. The peculiar velocity field is of constant divergence nature on small scale and irrotational on large scale. The statistical measures, i.e., correlation, structure, dispersion, and spectrum functions, are modeled on both small and large scales, respectively. Kinematic relations between statistical measures are fully developed for incompressible, constant divergence, and irrotational flow. Incompressible and constant divergence flow share the same kinematic relations for even order correlations. The limiting correlation of velocity $\mathrm{ρ_{L}=1/2}$ on the smallest scale ( r = 0) is a unique feature of collisionless flow (⁠$\mathrm{ρ_{L}=1}$ for incompressible flow). On large scale, transverse velocity correlation has an exponential form $T_{2}∝e^{–r/r_2}$ with a constant comoving scale r 2 =21.3 Mpc/h that may be related to the horizon size at matter–radiation equality. All other correlation, structure, dispersion, and spectrum functions for velocity, density, and potential fields are derived analytically from kinematic relations for irrotational flow. On small scale, longitudinal structure function follows one-fourth law of ${S}_{2}^{1}∝r^{1/4}$. All other statistical measures can be obtained from kinematic relations for constant divergence flow. Vorticity is negatively correlated for scale r between 1 and 7 Mpc/h. Divergence is negatively correlated for r > 30 Mpc/h that leads to a negative density correlation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

A unified catalogue-level reanalysis of stage-III cosmic shear surveys

ABSTRACT Cosmological parameter constraints from recent galaxy imaging surveys are reaching percent-level accuracy on the effective amplitude of the lensing signal, S8. The upcoming Legacy Survey of Space and Time (LSST) of the Vera C. Rubin Observatory will produce subpercent level measurements of cosmological parameters, providing a milestone test of the ΛCDM model. To supply guidance to the upcoming LSST analysis, it is important to understand thoroughly the results from different recent galaxy imaging surveys and assess their consistencies. In this work, we perform a unified catalogue-level reanalysis of three cosmic shear data sets: the first year data from the Dark Energy Survey (DES-Y1), the 1000 deg2 data set from the Kilo-Degree Survey (KiDS-1000), and the first year data from the Hyper Suprime-Cam Subaru Strategic Program (HSC-Y1). We utilize a pipeline developed and rigorously tested by the LSST Dark Energy Science Collaboration to perform the reanalysis and assess the robustness of the results to analysis choices. We find the S8 constraint to be robust to two different small-scale modelling approaches, and varying choices of cosmological priors. Our unified analysis allows the consistency of the surveys to be rigorously tested, and we find the three surveys to be statistically consistent. Due to the partially overlapping footprint, we model the cross-covariance between KiDS-1000 and HSC-Y1 approximately when combining all three data sets, resulting in a 1.6–1.9 per cent constraint on S8 given different assumptions on the cross-covariance.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Quantification of the Radiative Effect of Aerosol–Cloud Interactions in Shallow Continental Cumulus Clouds

The indirect radiative effect of aerosol variability on shallow cumulus clouds is realized in nature with considerable concurrent meteorological variability. Large-eddy simulations constrained by observations at a continental site in Oklahoma are performed to represent the variability of different meteorological states on days with different aerosol conditions. The total radiative effect of this natural covariation between aerosol and other meteorological drivers of total cloud amount and albedo is quantified. The changes to these bulk quantities are used to understand the response of the cloud radiative effect to aerosol–cloud interactions (ACI) in the context of concurrent processes, as opposed to attempting to untangle the effect of individual processes on a case-by-case basis. Mutual information (MI) analysis suggests that meteorological variability masks the strength of the relationship between cloud drop number concentration and the cloud radiative effect. This is shown to be mostly due to variation in solar zenith angle and cloud field horizontal heterogeneity masking the relationship between cloud drop number and cloud albedo. By combining MI and more traditional differential analyses, a framework to identify important modes of covariation between aerosol, clouds, and meteorological conditions is developed. This shows that accounting for solar zenith angle variation and implementing an albedo bias correction increases the detectability of the radiative effects of ACI in simulations of shallow cumulus.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Nuclear data uncertainty propagation and modeling uncertainty impact evaluation in neutronics core simulation

Uncertainty analysis is a critical requirement in reactor simulation as it is used to quantify the reliability of best-estimate calculation. A comprehensive uncertainty analysis should characterize all sources of uncertainties in a computationally-feasible and scientifically-defendable manner. Here we employ a well-established reduced order modeling (ROM) based uncertainty quantification methodology to propagate uncertainties throughout neutronic calculations. ROM relies on recent advances in randomized data mining techniques applied to large data streams. In our proposed implementation, the nuclear data uncertainties are first propagated from multi-group level through lattice physics calculation to generate few-group parameter uncertainties, described using a vector of mean values and a covariance matrix. Employing an ROM-based compression of the covariance matrix, the few-group uncertainties are then propagated through downstream core simulation in a computationally efficient manner. This straightforward approach, albeit efficient as compared to brute force forward and/or adjoint-based methods, often employs a number of assumptions that have been unquestioned in the literature of neutronic uncertainty analysis. This manuscript argues that these assumptions could introduce another source of uncertainty referred to as modeling uncertainties, whose magnitude needs to be quantified in tandem with nuclear data uncertainties. Thus, our primary goal is to explore the interactions between these two uncertainty sources in order to assess whether modeling uncertainties have an impact on parameter uncertainties. To explore this endeavor, the impact of a number of modeling assumptions on core attributes uncertainties is quantified. The study employs a CANDU reactor model, with Serpent and NEWT as lattice physics solvers and NESTLE-C as core simulator. The modeling assumptions investigated include those related with the uncertainty propagation method employed, e.g., deterministic vs. stochastic, the few-group energy structure employed to represent the cross-sections, the resonance treatment in lattice physics calculation, the reference values for the cross-section, and the number of samples employed to render ROM compression. Results indicate that some of the modeling assumptions could have a non-negligible impact on the core responses propagated uncertainties, highlighting the need for a more comprehensive approach to combine parameter and modeling uncertainties.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Demonstration of TOFFEE: A Response Uncertainty Quantification Tool

A key characteristic in neutron transport is nuclear data. Cross-section uncertainty is not used in MCNP6.3 to propagate response uncertainty without external analysis. Here, the TOol For Fast Error Estimation (TOFFEE) is a Python-based code developed to automate the propagation of cross-section uncertainty for MCNP evaluations. TOFFEE implements the sandwich rule to calculate the uncertainty from cross sections with sensitivity coefficients from MCNP6.3 and ENDF/B covariance data. In this paper, TOFFEE has been tested with benchmark experiments, and it has been compared to the uncertainty quantification capabilities of Sampler and TSUNAMI, within SCALE, to verify the application’s capabilities.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Vegetation classification map and covariates associated with NEON AOP survey, East River, CO 2018

This package includes geospatial data layers developed to investigate how environmental gradients—specifically topography and near-surface soil properties—drive the spatial arrangement of dominant plant communities in mountainous watersheds. The geospatial products, which support the analysis of these ecological relationships, are derived from airborne hyperspectral and LiDAR datasets acquired by the National Ecological Observatory Network (NEON) Airborne Observation Platform (AOP), in conjunction with an extensive ground field campaign conducted in summer 2018. This work is part of the DOE Watershed Function Science Focus Area (SFA) and features geospatial datasets developed based on observations and ground data collected at East River, Colorado, in collaboration with the National Ecological Observatory Network (NEON) Airborne Observation Platform (AOP) survey in June 2018. Classification Map: - Classification Map (PNG, GeoTIFF): Derived from hyperspectral and LiDAR airborne data using a machine learning approach. - Class Code Mapper (CSV): Associates pixel values with corresponding vegetation/non-vegetation classes. - Classification Reference Data (CSV): Reference data used in the machine learning procedure. LiDAR-Derived Products: - Topographical Metrics (GeoTIFFs): Elevation, slope, curvature, TWI, TPI, solar insolation, and canopy height model (CHM), smoothed with a 5x5 pixel window. Vegetation Indices: - GeoTIFFs of NDVI, NDNI, NDWI: Vegetation indices derived from hyperspectral data. Urban Masks: - Urban Mask (GeoTIFF): Applied to the mapping to convert bare soil classes to urban classes. Software Compatibility: GeoTIFFs: Can be visualized with GIS software or libraries that support GeoTIFF images. CSV Files: Can be opened with any software that handles comma-separated values. The FLMD file provides details and links to the source datasets used to derive the products. The manuscript (in the Method session) provides details on how each product was derived. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231. Update on 2026-03-25: Since the original dataset publication date of 02/28/2020, this package has a new classification map derived by an improved methodology. This update also includes additional ground data that improved the representation of some of the communities. See the methods for further details on what has changed between versions.

2018 NEON and 2025 CHESS Campaigns↗
DOE OSTI2021

Identifying dominant environmental predictors of freshwater wetland methane fluxes across diurnal to seasonal time scales

While wetlands are the largest natural source of methane (CH 4 ) to the atmosphere, they represent a large source of uncertainty in the global CH 4 budget due to the complex biogeochemical controls on CH 4 dynamics. Here we present, to our knowledge, the first multi-site synthesis of how predictors of CH 4 fluxes (FCH4) in freshwater wetlands vary across wetland types at diel, multiday (synoptic), and seasonal time scales. In this work we used several statistical approaches (correlation analysis, generalized additive modeling, mutual information, and random forests) in a wavelet-based multi-resolution framework to assess the importance of environmental predictors, nonlinearities and lags on FCH4 across 23 eddy covariance sites. Seasonally, soil and air temperature were dominant predictors of FCH4 at sites with smaller seasonal variation in water table depth (WTD). In contrast, WTD was the dominant predictor for wetlands with smaller variations in temperature (e.g., seasonal tropical/subtropical wetlands). Changes in seasonal FCH4 lagged fluctuations in WTD by ~17 ± 11 days, and lagged air and soil temperature by median values of 8 ± 16 and 5 ± 15 days, respectively. Temperature and WTD were also dominant predictors at the multiday scale. Atmospheric pressure (PA) was another important multiday scale predictor for peat-dominated sites, with drops in PA coinciding with synchronous releases of CH 4 . At the diel scale, synchronous relationships with latent heat flux and vapor pressure deficit suggest that physical processes controlling evaporation and boundary layer mixing exert similar controls on CH 4 volatilization, and suggest the influence of pressurized ventilation in aerenchymatous vegetation. In addition, 1- to 4-h lagged relationships with ecosystem photosynthesis indicate recent carbon substrates, such as root exudates, may also control FCH4. By addressing issues of scale, asynchrony, and nonlinearity, this work improves understanding of the predictors and timing of wetland FCH4 that can inform future studies and models, and help constrain wetland CH 4 emissions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Airborne flux measurements of ammonia over the southern Great Plains using chemical ionization mass spectrometry

Abstract. Ammonia (NH3) is an abundant trace gas in the atmosphere and an important player in atmospheric chemistry, aerosol formation and the atmosphere–surface exchange of nitrogen. The accurate determination of NH3 emission rates remains a challenge, partly due to the propensity of NH3 to interact with instrument surfaces, leading to high detection limits and slow response times. In this paper, we present a new method for quantifying ambient NH3, using chemical ionization mass spectrometry (CIMS) with deuterated benzene cations as reagents. The setup aimed at limiting sample–surface interactions and achieved a 1σ precision of 10–20 pptv and an immediate 1/e response rate of < 0.4 s, which compares favorably to the existing state of the art. The sensitivity exhibited an inverse humidity dependence, in particular in relatively dry conditions. Background of up to 10 % of the total signal required consideration as well, as it responded on the order of a few minutes. To showcase the method's capabilities, we quantified NH3 mixing ratios from measurements obtained during deployment on a Gulfstream I aircraft during the HI-SCALE (Holistic Interactions of Shallow Clouds, Aerosols, and Land-Ecosystems) field campaign in rural Oklahoma during May 2016. Typical mixing ratios were 1–10 parts per billion by volume (ppbv) for the boundary layer and 0.1–1 ppbv in the lower free troposphere. Sharp plumes of up to tens of ppbv of NH3 were encountered as well. We identified two of their sources as a large fertilizer plant and a cattle farm, and our mixing ratio measurements yielded upper bounds of 350 ± 50 and 0.6 kg NH3 h−1 for their respective momentary source rates. The fast response of the CIMS also allowed us to derive vertical NH3 fluxes within the turbulent boundary layer via eddy covariance, for which we chiefly used the continuous wavelet transform technique. As expected for a region dominated by agriculture, we observed predominantly upward fluxes, implying net NH3 emissions from the surface. The corresponding analysis focused on the most suitable flight, which contained two straight-and-level legs at ∼ 300 m above ground. We derived NH3 fluxes between 1 and 11 mol km−2 h−1 for these legs, at an effective spatial resolution of 1–2 km. The analysis demonstrated how flux measurements benefit from suitably arranged flight tracks with sufficiently long straight-and-level legs, and it explores the detrimental effect of measurement discontinuities. Following flux footprint estimations, comparison to the NH3 area emissions inventory provided by the U.S. Environmental Protection Agency indicated overall agreement but also the absence of some sources, for instance the identified cattle farm. Our study concludes that high-precision CIMS measurements are a powerful tool for in situ measurements of ambient NH3 mixing ratios, and even allow for the airborne mapping of the air–surface exchange of NH3.

58 GEOSCIENCES↗
DOE OSTI2022

Airborne flux measurements of ammonia over the Southern Great Plains using chemical ionization mass spectrometry

Ammonia (NH 3 ) is an abundant trace gas in the atmosphere and an important player in atmospheric chemistry, aerosol formation and the atmosphere-surface exchange of nitrogen. It is recognized as a major source of aerosol pollution, and it may limit the formation of cloud nuclei in remote or cold parts of the atmosphere. For soil and plants, NH 3 -mediated nitrogen can act as a harmful pollutant or as a desirable nutrient, mostly in natural and agricultural settings, respectively. Agriculture is also the main source of atmospheric NH 3 via volatilization from fertilizers and manure processing in livestock farming. The accurate determination of NH 3 emission rates remains a challenge, partly due to the propensity of NH 3 to interact with instrument surfaces leading to high detection limits and slow response times. In this paper, we present a new method for quantifying ambient NH 3 , using chemical ionization mass spectrometry (CIMS) with deuterated benzene cations as reagents. The setup aimed at limiting sample-surface interactions and achieved a 1-σ precision of 10–20 pptv and an immediate 1/e response rate < 0.4 s, which compares favorably to the existing state of the art. The sensitivity exhibited an inverse humidity dependence, in particular in relatively dry conditions. Background of up to 10 % of the total signal required consideration as well, as it responded on the order of a few minutes. To showcase the method’s capabilities, we quantified NH3 mixing ratios from measurements obtained during deployment on a Gulfstream I aircraft during the HI-SCALE (Holistic Interactions of Shallow Clouds, Aerosols and Land Ecosystems) field campaign in rural Oklahoma during May 2016. Typical mixing ratios were 1–10 parts per billion by volume (ppbv) for the boundary layer and 0.1–1 ppbv in the lower free troposphere. Sharp plumes of up to 10s of ppbv of NH3 were encountered as well. We identified two of their sources as a large fertilizer plant and a cattle farm, and our mixing ratio measurements yielded upper bounds of 350 ± 50 and 0.6 kg NH 3 h –1 for their respective momentary source rates. The fast response of the CIMS also allowed us to derive vertical NH 3 fluxes within the turbulent boundary layer via eddy covariance, for which we chiefly used the continuous wavelet transform technique. As expected for a region dominated by agriculture, we observed predominantly upward fluxes, implying net NH 3 emissions from surface. The corresponding analysis focused on the most suitable flight, which contained two straight-and-level legs at ~300 m above ground. We derived NH 3 fluxes between –4 and 18 mol km –2 h –1 for these legs, at an effective spatial resolution of 1–2 km. The analysis demonstrated how flux measurements benefit from suitably arranged flight tracks with sufficiently long straight-and-level legs, and explores the detrimental effect of measurement discontinuities. Following flux footprint estimations, comparison to the NH 3 area emissions inventory provided by the US Environmental Protection Agency indicated overall agreement, but also the absence of some sources, for instance the identified cattle farm. Our study concludes that high-precision CIMS measurements are a powerful tool for in-situ measurements of ambient NH 3 mixing ratios, and even allow for the airborne mapping of the air-surface exchange of NH 3 .

54 ENVIRONMENTAL SCIENCES↗