Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Confidence interval”

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

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

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Comparing traditional and Bayesian approaches to ecological meta‐analysis

Abstract Despite the wide application of meta‐analysis in ecology, some of the traditional methods used for meta‐analysis may not perform well given the type of data characteristic of ecological meta‐analyses. We reviewed published meta‐analyses on the ecological impacts of global climate change, evaluating the number of replicates used in the primary studies ( n i ) and the number of studies or records ( k ) that were aggregated to calculate a mean effect size. We used the results of the review in a simulation experiment to assess the performance of conventional frequentist and Bayesian meta‐analysis methods for estimating a mean effect size and its uncertainty interval. Our literature review showed that n i and k were highly variable, distributions were right‐skewed and were generally small (median n i = 5, median k = 44). Our simulations show that the choice of method for calculating uncertainty intervals was critical for obtaining appropriate coverage (close to the nominal value of 0.95). When k was low (<40), 95% coverage was achieved by a confidence interval (CI) based on the t distribution that uses an adjusted standard error (the Hartung–Knapp–Sidik–Jonkman, HKSJ), or by a Bayesian credible interval, whereas bootstrap or z distribution CIs had lower coverage. Despite the importance of the method to calculate the uncertainty interval, 39% of the meta‐analyses reviewed did not report the method used, and of the 61% that did, 94% used a potentially problematic method, which may be a consequence of software defaults. In general, for a simple random‐effects meta‐analysis, the performance of the best frequentist and Bayesian methods was similar for the same combinations of factors ( k and mean replication), though the Bayesian approach had higher than nominal (>95%) coverage for the mean effect when k was very low ( k < 15). Our literature review suggests that many meta‐analyses that used z distribution or bootstrapping CIs may have overestimated the statistical significance of their results when the number of studies was low; more appropriate methods need to be adopted in ecological meta‐analyses.

Pappalardo, Paula↗
DOE OSTI2023

Experimental Validation of a Command and Control Traffic Detection Model

Network intrusion detection systems (NIDS) are commonly used to detect malware communications, including command-and-control (C2) traffic from botnets. NIDS performance assessments have been studied for decades, but mathematical modeling has rarely been used to explore NIDS performance. This paper details a mathematical model that describes a NIDS performing packet inspection and its detection of malware's C2 traffic. Here, the paper further describes an emulation testbed and a set of cyber experiments that used the testbed to validate the model. These experiments included a commonly used NIDS (Snort) and traffic with contents from a pervasive malware (Emotet). Results are presented for two scenarios: a nominal scenario and a “stressed” scenario in which the NIDS cannot process all incoming packets. Model and experiment results match well, with model estimates mostly falling within 95 % confidence intervals on the experiment means. Model results were produced 70-3000 times faster than the experimental results. Consequently, the model's predictive capability could potentially be used to support decisions about NIDS configuration and effectiveness that require high confidence results, quantification of uncertainty, and exploration of large parameter spaces. Furthermore, the experiments provide an example for how emulation testbeds can be used to validate cyber models that include stochastic variability.

mathematical model↗
DOE OSTI2021

Defense Waste Processing Facility Nitric-Glycolic Flowsheet Chemical Process Cell Chemistry: Part 2

The conversions of nitrite to nitrate, the destruction of glycolate, and the conversion of glycolate to formate and oxalate were modeled for the Nitric-Glycolic flowsheet using data from Chemical Process Cell (CPC) simulant runs conducted by Savannah River National Laboratory (SRNL) from 2011 to 2016. The goal of this work was to develop empirical correlation models to predict these values from measurable variables from the chemical process so that these quantities could be predicted a-priori from the sludge or simulant composition and measurable processing variables. The need for these predictions arises from the need to predict the REDuction/OXidation (REDOX) state of the glass from the Defense Waste Processing Facility (DWPF) melter. This report summarizes the work on these correlations based on the aforementioned data. Previous work on these correlations was documented in a technical report covering data from 2011-2015. This current report supersedes this previous report. Further refinement of the models as additional data are collected is recommended. The glass REDOX depends on the concentrations of nitrate and manganese (oxidants), and of glycolate, formate, oxalate, carbon, and antifoam (reductants) in the melter feed. The waste sludge contains nitrite, nitrate, manganese (Mn), and oxalate. Virtually all of the nitrite is converted to nitrate or NO+NO 2 +N 2 O gases in the CPC. The portion of the nitrite converted to nitrate increases the amount of nitrate in the sludge. The amount of glycolate in the final melter feed depends on the amount of the glycolic acid feed that is destroyed. Similarly, the amounts of formate and oxalate formed during the decomposition of glycolic acid are required. The material balance on carbon was found to not close in most cases. Generally, there was less carbon at the end of testing compared to the inputs. The most uncertain product variable was glycolate, so material balances were performed where the glycolate concentration was adjusted, usually upward, to close the balance. Correlation versus the original, as-measured, data was generally poor, but correlation against the material balance adjusted values was greatly improved. It was also shown that the correlation of the measured REDOX versus the predicted REDOX was much better when the material balance adjusted glycolate values were used. Three data series were primarily used during the regressions of the data; these series were 1) Sludge Batch 9 NG flowsheet simulant runs NG51-62 (SB9-NG); 2) Scaled Runs + Bounding Hydrogen Runs (SR+BH); and 3) Runs GN43-50 and 57 (43-50,57). The glycolate destruction was found to correlate with acid stoichiometry (AS), percent reducing acid (PRA), and for some data series, headspace to simulant volume ratio (HSV), mercury (Hg), and nitrate. Although glycolate destruction for pairs of data series (e.g., [SB9-NG] and [SR+BH]) were found to depend on HSV, the combination of all three data series was not found to have significant dependence on this variable. The best model for glycolate destruction depended on AS, nitrate, and Hg. This model predicted the product glycolate compositions of the data to within 92-106%. The conversion of glycolate to formate was high when noble metals and Hg were not present, with values up to 100%. When noble metals and Hg were present, this conversion ranged from zero to 7%, and was dependent on AS. Lower AS gave higher conversions to formate. The conversion to oxalate was found to depend on the AS and the initial concentration of nitrite. An alternative fit versus AS and the form of ruthenium (Ru) used is a possible alternative. This fit was somewhat less statistically significant. This second model predicts that more oxalate is formed when Ru-nitrosyl nitrate is used rather than Ru chloride. The conversion of glycolate to oxalate ranged from zero to 6%. The conversion of nitrite to nitrate depended primarily on AS and PRA, with HSV and Hg being significant when these variables were varied. For multiple series of data, nitrite was also needed to SRNL-STI-2017-00172 5HYLVLRQ viL distinguish between data series, and the effect of HSV became insignificant. The best model for nitrite to nitrate conversion depended on AS, PRA, nitrite, and Hg. The 95% confidence intervals on the predicted values of glycolate destruction, glycolate to oxalate conversion, and nitrite to nitrate conversion were used to determine the uncertainty in the predicted REDOX when starting with only the composition of the sludge, AS, and PRA. Using the 95% confidences on an individual value (that is the confidence in getting a particular value for one single test as opposed to what the mean would be for multiple tests), the uncertainty in the predicted REDOX was calculated. The uncertainty in the actual product composition glycolate, oxalate, formate, and nitrate concentrations translated to an uncertainty in the REDOX value of ±0.1,which is approximately the uncertainty claimed in the REDOX model itself.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Probing the consistency of cosmological contours for supernova cosmology

As the scale of cosmological surveys increases, so does the complexity in the analyses. This complexity can often make it difficult to derive the underlying principles, necessitating statistically rigorous testing to ensure the results of an analysis are consistent and reasonable. This is particularly important in multi-probe cosmological analyses like those used in the Dark Energy Survey (DES) and the upcoming Legacy Survey of Space and Time, where accurate uncertainties are vital. In this paper, we present a statistically rigorous method to test the consistency of contours produced in these analyses and apply this method to the Pippin cosmological pipeline used for type Ia supernova cosmology with the DES. We make use of the Neyman construction, a frequentist methodology that leverages extensive simulations to calculate confidence intervals, to perform this consistency check. A true Neyman construction is too computationally expensive for supernova cosmology, so we develop a method for approximating a Neyman construction with far fewer simulations. We find that for a simulated dataset, the 68% contour reported by the Pippin pipeline and the 68% confidence region produced by our approximate Neyman construction differ by less than a percent near the input cosmology; however, they show more significant differences far from the input cosmology, with a maximal difference of 0.05 in $Ω$ M and 0.07 in w. In conclusion, this divergence is most impactful for analyses of cosmological tensions, but its impact is mitigated when combining supernovae with other cross-cutting cosmological probes, such as the cosmic microwave background.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

Early calendar life and health prediction of silicon batteries via machine learning with uncertainty quantification

Lithium-ion batteries with silicon anodes promise high energy density but are limited by calendar lifetime. Reducing the long iteration time to obtain experimental results requires predicting calendar lifetime early in a cell's life. In this study, we demonstrate that lightweight machine learning models with feature engineering can provide calendar lifetime estimates from early electrochemical signals. After 1 month of electrochemical aging, the best models achieve 10% error in calendar-life prediction and can separate "bad" from "good" lifetime cells with a mean F1 score of 0.857. As battery systems exhibit inherent variability, four methods for uncertainty quantification are compared, and confidence intervals are demonstrated with an uncertainty of +-3.6 months in lifetime prediction. A feature importance analysis indicates that early patterns in voltage decay are the strongest indicators of calendar lifetime. Finally, this modeling approach has high error when generalizing to new electrode chemistries or testing conditions but with appropriately low confidence.

25 ENERGY STORAGE↗
DOE OSTI2023

A stochastic model of future extreme temperature events for infrastructure analysis

Applying extreme temperature events for future conditions is not straightforward for infrastructure resilience analyses. This work introduces a stochastic model that fills this gap. The model uses at least 50 years of daily extreme temperature records, climate normals with 10%-90% confidence intervals, and shifts/offsets for increased frequency and intensity of heat wave events. Intensity and frequency are shifted based on surface temperature anomaly from 1850-1900 for 32 models from CMIP6. A case study for Worcester, Massachusetts passed 85% of cases using the two-sided Kolmogorov-Smirnov -value test with 95% confidence for both temperature and duration. Future shifts for several climate scenarios to 2020, 2040, 2060, and 2080 had acceptable errors between the shifted model and 10- and 50-year extreme temperature event thresholds with the largest error being 2.67 degrees C. The model is likely to be flexible enough for other patterns of extreme weather such as extreme precipitation and hurricanes.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2025

Personalized and uncertainty-aware coronary hemodynamics simulations: From Bayesian estimation to improved multi-fidelity uncertainty quantification

Non-invasive simulations of coronary hemodynamics have improved clinical risk stratification and treatment outcomes for coronary artery disease, compared to relying on anatomical imaging alone. However, simulations typically use empirical approaches to distribute total coronary flow amongst the arteries in the coronary tree, which ignores patient variability, the presence of disease, and other clinical factors. Further, uncertainty in the clinical data often remains unaccounted for in the modeling pipeline. We present an end-to-end uncertainty-aware pipeline to (1) personalize coronary flow simulations by incorporating vessel-specific coronary flows as well as cardiac function; and (2) predict clinical and biomechanical quantities of interest with improved precision, while accounting for uncertainty in the clinical data. We assimilate patient-specific measurements of myocardial blood flow from clinical CT myocardial perfusion imaging to estimate branch-specific coronary artery flows. Simulated noise in the clinical data is used to estimate the joint posterior distributions of the model parameters using adaptive Markov Chain Monte Carlo sampling. Additionally, the posterior predictive distribution for the relevant quantities of interest is determined using a new approach combining multi-fidelity Monte Carlo estimation with non-linear, data-driven dimensionality reduction. This leads to improved correlations between high- and low-fidelity model outputs. Our framework accurately recapitulates clinically measured cardiac function as well as branch-specific coronary flows under measurement noise uncertainty. We observe substantial reductions in confidence intervals for estimated quantities of interest compared to single-fidelity Monte Carlo estimation and state-of-the-art multi-fidelity Monte Carlo methods. This holds especially true for quantities of interest that showed limited correlation between the low- and high-fidelity model predictions. In addition, the proposed multi-fidelity Monte Carlo estimators are significantly cheaper to compute than traditional estimators, under a specified confidence level or variance. The proposed pipeline for personalized and uncertainty-aware predictions of coronary hemodynamics is based on routine clinical measurements and recently developed techniques for CT myocardial perfusion imaging. The proposed pipeline offers significant improvements in precision and reduction in computational cost.

Bayesian parameter estimation↗
DOE OSTI2024

A new data-driven map predicts substantial undocumented peatland areas in Amazonia

Tropical peatlands are among the most carbon-dense terrestrial ecosystems yet recorded. Collectively, they comprise a large but highly uncertain reservoir of the global carbon cycle, with wide-ranging estimates of their global area (441 025–1700 000 km 2 ) and below-ground carbon storage (105–288 Pg C). Substantial gaps remain in our understanding of peatland distribution in some key regions, including most of tropical South America. Here we compile 2413 ground reference points in and around Amazonian peatlands and use them alongside a stack of remote sensing products in a random forest model to generate the first field-data-driven model of peatland distribution across the Amazon basin. Our model predicts a total Amazonian peatland extent of 251 015 km 2 (95th percentile confidence interval: 128 671–373 359), greater than that of the Congo basin, but around 30% smaller than a recent model-derived estimate of peatland area across Amazonia. The model performs relatively well against point observations but spatial gaps in the ground reference dataset mean that model uncertainty remains high, particularly in parts of Brazil and Bolivia. For example, we predict significant peatland areas in northern Peru with relatively high confidence, while peatland areas in the Rio Negro basin and adjacent south-western Orinoco basin which have previously been predicted to hold Campinarana or white sand forests, are predicted with greater uncertainty. Similarly, we predict large areas of peatlands in Bolivia, surprisingly given the strong climatic seasonality found over most of the country. Very little field data exists with which to quantitatively assess the accuracy of our map in these regions. Data gaps such as these should be a high priority for new field sampling. This new map can facilitate future research into the vulnerability of peatlands to climate change and anthropogenic impacts, which is likely to vary spatially across the Amazon basin.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Calculating Interval Uncertainties for Calibration Standards That Drift with Time

Calibrated values of many devices exhibit predictable drift over time. To provide an uncertainty statement valid over the entire calibration interval, one must account for drift. In this article, a method of accounting for drift is proposed based on guidance in the Guide to Expression of Uncertainty in Measurement. An additional uncertainty term is computed using a linear regression of historical measurement data, which is included along with the time-of-test uncertainty. This method is evaluated by analyzing its average out-of-tolerance (OOT) rate using a Monte Carlo simulation, which results in the desired 5% average OOT rate when the total uncertainty is expanded to a 95% confidence interval.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Optimizing cloud motion estimation on the edge with phase correlation and optical flow

Abstract. Phase correlation (PC) is a well-known method for estimating cloud motion vectors (CMVs) from infrared and visible spectrum images. Commonly, phase shift is computed in the small blocks of the images using the fast Fourier transform. In this study, we investigate the performance and the stability of the blockwise PC method by changing the block size, the frame interval, and combinations of red, green, and blue (RGB) channels from the total sky imager (TSI) at the United States Atmospheric Radiation Measurement user facility's Southern Great Plains site. We find that shorter frame intervals, followed by larger block sizes, are responsible for stable estimates of the CMV, as suggested by the higher autocorrelations. The choice of RGB channels has a limited effect on the quality of CMVs, and the red and the grayscale images are marginally more reliable than the other combinations during rapidly evolving low-level clouds. The stability of CMVs was tested at different image resolutions with an implementation of the optimized algorithm on the Sage cyberinfrastructure test bed. We find that doubling the frame rate outperforms quadrupling the image resolution in achieving CMV stability. The correlations of CMVs with the wind data are significant in the range of 0.38–0.59 with a 95 % confidence interval, despite the uncertainties and limitations of both datasets. A comparison of the PC method with constructed data and the optical flow method suggests that the post-processing of the vector field has a significant effect on the quality of the CMV. The raindrop-contaminated images can be identified by the rotation of the TSI mirror in the motion field. The results of this study are critical to optimizing algorithms for edge-computing sensor systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Metabolic Costs of Walking with Weighted Vests

ABSTRACT Introduction The US Army Load Carriage Decision Aid (LCDA) metabolic model is used by militaries across the globe and is intended to predict physiological responses, specifically metabolic costs, in a wide range of dismounted warfighter operations. However, the LCDA has yet to be adapted for vest-borne load carriage, which is commonplace in tactical populations, and differs in energetic costs to backpacking and other forms of load carriage. Purpose The purpose of this study is to develop and validate a metabolic model term that accurately estimates the effect of weighted vest loads on standing and walking metabolic rate for military mission-planning and general applications. Methods Twenty healthy, physically active military-age adults (4 women, 16 men; age, 26 ± 8 yr old; height, 1.74 ± 0.09 m; body mass, 81 ± 16 kg) walked for 6 to 21 min with four levels of weighted vest loading (0 to 66% body mass) at up to 11 treadmill speeds (0.45 to 1.97 m·s −1 ). Using indirect calorimetry measurements, we derived a new model term for estimating metabolic rate when carrying vest-borne loads. Model estimates were evaluated internally byk-fold cross-validation and externally against 12 reference datasets (264 total participants). We tested if the 90% confidence interval of the mean paired difference was within equivalence limits equal to 10% of the measured walking metabolic rate. Estimation accuracy, precision, and level of agreement were also evaluated by the bias, standard deviation of paired differences, and concordance correlation coefficient (CCC), respectively. Results Metabolic rate estimates using the new weighted vest term were statistically equivalent (P< 0.01) to measured values in the current study (bias, −0.01 ± 0.54 W·kg −1 ; CCC, 0.973) as well as from the 12 reference datasets (bias, −0.16 ± 0.59 W·kg −1 ; CCC, 0.963). Conclusions The updated LCDA metabolic model calculates accurate predictions of metabolic rate when carrying heavy backpack and vest-borne loads. Tactical populations and recreational athletes that train with weighted vests can confidently use the simplified LCDA metabolic calculator provided as Supplemental Digital Content to estimate metabolic rates for work/rest guidance, training periodization, and nutritional interventions.

Sport Sciences↗
DOE OSTI2020

Bayesian inference of nuclear symmetry energy from measured and imagined neutron skin thickness in Sn 116 , 118 , 120 , 122 , 124 , 130 , 132 , Pb 208 , and Ca 48

The neutron skin thickness Δr np in heavy nuclei has been known as one of the most sensitive terrestrial probes of the nuclear symmetry energy E sym (ρ) around $\frac{2}{3}$ of the saturation density ρ 0 of nuclear matter. Existing neutron skin data mostly from hadronic observables suffer from large uncertainties and their extraction from experiments are often strongly model dependent. While waiting eagerly for the promised model-independent and high-precision neutron skin data for 208 Pb and 48 Ca from the parity-violating electron scattering experiments (PREX-II and CREX at JLab as well as MREX at MESA), within the Bayesian statistical framework using the Skyrme-Hartree-Fock model we infer the posterior probability distribution functions (PDFs) of the slope parameter L of the nuclear symmetry energy at ρ 0 from imagined Δr np ( 208 Pb)=0.15, 0.20, and 0.30 fm with a 1σ error bar of 0.02, 0.04, and 0.06 fm, respectively, as well as Δr np ( 48 Ca)=0.12, 0.15, and 0.25 fm, with different 1σ error bar of 0.01 and 0.02 fm, respectively. The results are compared with the PDFs of L inferred using the same approach from the available Δr np data for 116, 118, 120, 122, 124, 130, 132 Sn from hadronic probes. They are also compared with results from a recent Bayesian analysis of the radius and tidal deformability data of canonical neutron stars from GW170817 and NICER. The neutron skin data for Sn isotopes gives L = 45.5 $^{+ 26.5}_{-21.6}$ MeV surrounding its mean value or L = 53 . 4 $^{+ 18.6}_{ -29.5}$ MeV surrounding its maximum a posteriori value, respectively, with the latter smaller than but consistent with the L = 66 $^{+ 12}_{-20}$ MeV from the neutron star data within their 68% confidence intervals. We found that Δr np = 0.17 –0.18 fm in 208 Pb with an error bar of about 0.02 fm leads to a PDF of L compatible with that from analyzing the Sn data. To provide additionally useful information on L extracted from the Δr np of Sn isotopes, the experimental error bar of Δr np in 208 Pb should be at least smaller than 0.06 fm aimed by some current experiments. In addition, the Δr np ( 48 Ca) needs to be larger than 0.15 fm but smaller than 0.25 fm to be compatible with the Sn and/or neutron star results. To further improve our current knowledge about L and distinguish its PDFs in the examples considered, even higher precisions leading to significantly less than ±20 MeV error bars for L at 68% confidence level are necessary.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Cosmological neutrino mass: a frequentist overview in light of DESI

We derive constraints on the neutrino mass using a variety of recent cosmological datasets, including DESI BAO, the full-shape analysis of the DESI matter power spectrum and the one-dimensional power spectrum of the Lyman-α forest (P1D) from eBOSS quasars as well as the cosmic microwave background (CMB). The constraints are obtained in the frequentist formalism by constructing profile likelihoods and applying the Feldman-Cousins prescription to compute confidence intervals. This method avoids potential prior and volume effects that may arise in a comparable Bayesian analysis. Parabolic fits to the profiles allow one to distinguish changes in the upper limits from variations in the constraining power σ of the different data combinations. We find that all profiles in the ΛCDM model are cut off by the ∑m ν ≥ 0 bound, meaning that the corresponding parabolas reach their minimum in the unphysical sector. The most stringent 95% C.L. upper limit is obtained by the combination of DESI DR2 BAO, Planck PR4 and CMB lensing at 53 meV, below the minimum of 59 meV set by the normal ordering. The corresponding constraining power σ is 43 meV, which highlights the importance of the cut-off by negative values in the determination of the upper limit. Extending ΛCDM to non-zero curvature and w 0 w a CDM relaxes the constraints past 59 meV again, but only w 0 w a CDM exhibits profiles with a minimum at a positive value. Additionally, we extend the formalism to constrain the lightest neutrino mass. For DESI DR2 BAO, Planck PR4 and CMB lensing, we find confidence limits at 20 and 19 meV for normal and inverted ordering, respectively. Using a combination of DESI DR1 full-shape, BBN and eBOSS Lyman-α P1D, we successfully constrain the neutrino mass independently of the CMB. This combination yields m l ≤ 97 and 98 meV in the normal and inverted orderings, and total neutrino mass ∑m ν ≤ 285 meV (95% C.L.). The addition of DESI full-shape or Lyman-α P1D to CMB and DESI BAO results in small but noticeable improvement of the constraining power of the data. Lyman-α free-streaming measurements especially improve the constraint. Since they are based on eBOSS data, this sets a promising precedent for upcoming DESI data.

Frequentist statistics↗
DOE OSTI2021

A comprehensive and synthetic dataset for global, regional, and national greenhouse gas emissions by sector 1970–2018 with an extension to 2019

To track progress towards keeping global warming well below 2 °C or even 1.5 °C, as agreed in the Paris Agreement, comprehensive up-to-date and reliable information on anthropogenic emissions and removals of greenhouse gas (GHG) emissions is required. Here we compile a new synthetic dataset on anthropogenic GHG emissions for 1970–2018 with a fast-track extension to 2019. Our dataset is global in coverage and includes CO 2 emissions, CH 4 emissions, N 2 O emissions, as well as those from fluorinated gases (F-gases: HFCs, PFCs, SF 6 , NF 3 ) and provides country and sector details. We build this dataset from the version 6 release of the Emissions Database for Global Atmospheric Research (EDGAR v6) and three bookkeeping models for CO 2 emissions from land use, land-use change, and forestry (LULUCF). We assess the uncertainties of global greenhouse gases at the 90 % confidence interval (5th–95th percentile range) by combining statistical analysis and comparisons of global emissions inventories and top-down atmospheric measurements with an expert judgement informed by the relevant scientific literature. We identify important data gaps for F-gas emissions. The agreement between our bottom-up inventory estimates and top-down atmospheric-based emissions estimates is relatively close for some F-gas species (~ 10 % or less), but estimates can differ by an order of magnitude or more for others. Our aggregated F-gas estimate is about 10 % lower than top-down estimates in recent years. However, emissions from excluded F-gas species such as chlorofluorocarbons (CFCs) or hydrochlorofluorocarbons (HCFCs) are cumulatively larger than the sum of the reported species. Using global warming potential values with a 100-year time horizon from the Sixth Assessment Report by the Intergovernmental Panel on Climate Change (IPCC), global GHG emissions in 2018 amounted to 58 ± 6.1 GtCO 2 eq. consisting of CO 2 from fossil fuel combustion and industry (FFI) 38 ± 3.0 GtCO 2 , CO 2 -LULUCF 5.7 ± 4.0 GtCO 2 , CH 4 10 ± 3.1 GtCO 2 eq., N2O 2.6 ± 1.6 GtCO 2 eq., and F-gases 1.3 ± 0.40 GtCO 2 eq. Initial estimates suggest further growth of 1.3 GtCO 2 eq. in GHG emissions to reach 59 ± 6.6 GtCO 2 eq. by 2019. Our analysis of global trends in anthropogenic GHG emissions over the past 5 decades (1970–2018) highlights a pattern of varied but sustained emissions growth. There is high confidence that global anthropogenic GHG emissions have increased every decade, and emissions growth has been persistent across the different (groups of) gases. There is also high confidence that global anthropogenic GHG emissions levels were higher in 2009–2018 than in any previous decade and that GHG emissions levels grew throughout the most recent decade. While the average annual GHG emissions growth rate slowed between 2009 and 2018 (1.2 % yr –1 ) compared to 2000–2009 (2.4 % yr –1 ), the absolute increase in average annual GHG emissions by decade was never larger than between 2000–2009 and 2009–2018. Our analysis further reveals that there are no global sectors that show sustained reductions in GHG emissions. There are a number of countries that have reduced GHG emissions over the past decade, but these reductions are comparatively modest and outgrown by much larger emissions growth in some developing countries such as China, India, and Indonesia. There is a need to further develop independent, robust, and timely emissions estimates across all gases. As such, tracking progress in climate policy requires substantial investments in independent GHG emissions accounting and monitoring as well as in national and international statistical infrastructures. The data associated with this article (Minx et al., 2021) can be found at https://doi.org/10.5281/zenodo.5566761.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Iron Heterogeneity in Early Active Multiple Sclerosis Lesions

Objective: Multiple sclerosis (MS) is a heterogeneous inflammatory demyelinating disease. Iron distribution is altered in MS patients' brains, suggesting iron liberation within active lesions amplifies demyelination and neurodegeneration. Whether the amount and distribution of iron are similar or different among different MS immunopatterns is currently unknown. Methods: We used synchrotron X-ray fluorescence imaging, histology, and immunohistochemistry to compare the iron quantity and distribution between immunopattern II and III early active MS lesions. We analyzed archival autopsy and biopsy tissue from 21 MS patients. Results: Immunopattern II early active lesions contain 64% more iron (95% confidence interval [CI] = 17-127, p = 0.004) than immunopattern III lesions, and 30% more iron than the surrounding periplaque white matter (95% CI = 3–64%, p = 0.03). Iron in immunopattern III lesions is 28% lower than in the periplaque white matter (95% CI = -40 to -14%, p < 0.001). When normalizing the iron content of early active lesions to that of surrounding periplaque white matter, the ratio is significantly higher in immunopattern II (p < 0.001). Microfocused X-ray fluorescence imaging shows that iron in immunopattern II lesions localizes to macrophages, whereas macrophages in immunopattern III lesions contain little iron. Interpretation: Iron distribution and content are heterogeneous in early active MS lesions. Iron accumulates in macrophages in immunopattern II, but not immunopattern III lesions. This heterogeneity in the two most common MS immunopatterns may be explained by different macrophage polarization, origin, or different demyelination mechanisms, and paves the way for developing new or using existing iron-sensitive magnetic resonance imaging techniques to differentiate among immunopatterns in the general nonbiopsied MS patient population.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Associations between cumulative environmental quality and ten selected birth defects in Texas

Abstract Background Causes of most birth defects are largely unknown. Genetics, maternal factors (e.g., age, smoking) and environmental exposures have all been linked to some birth defects, including neural tube, oral cleft, limb reduction, and gastroschisis; however, the contribution of cumulative exposures across several environmental domains in association with these defects is not well understood. Methods The Environmental Quality Index (EQI) and its domains (air, water, land, sociodemographic, built) were used to estimate county‐level cumulative environmental exposures from 2006–2010 and matched to birth defects identified from Texas Birth Defects Registry and live birth records from births in years 2007–2010 ( N = 1,610,709). Poisson regression models estimated prevalence ratios (PR) and 95% confidence intervals (CI) for associations between 10 birth defects and the EQI. Results We observed some positive associations between worst environmental quality and neural tube, anencephaly, spina bifida, oral cleft, cleft palate, cleft lip with and without cleft palate, and gastroschisis [PR range: 1.12–1.55], but near null associations with limb reduction defects. Among domain specific results, we observed the strongest positive associations with the sociodemographic domain across birth defects but varied positive associations among the air and water domains, and negative or null associations with the land and built domains. Overall, few exposure‐response patterns were evident. Conclusions Our results highlight the complexities of cumulative, simultaneous environmental exposures in the prevalence rates of 10 selected birth defects. We were able to explore the impact of overall and domain specific environmental quality on birth defects and identify potential domain specific drivers of these associations.

Krajewski, Alison K.↗
DOE OSTI2024

Pan‐Cancer Survival Impact of Immune Checkpoint Inhibitors in a National Healthcare System

ABSTRACT Background The cumulative, health system‐wide survival benefit of immune checkpoint inhibitors (ICIs) is unclear, particularly among real‐world patients with limited life expectancies and among subgroups poorly represented on clinical trials. We sought to determine the health system‐wide survival impact of ICIs. Methods We identified all patients receiving PD‐1/PD‐L1 or CTLA‐4 inhibitors from 2010 to 2023 in the national Veterans Health Administration (VHA) system (ICI cohort) and all patients who received non‐ICI systemic therapy in the years before ICI approval (historical control). ICI and historical control cohorts were matched on multiple cancer‐related prognostic factors, comorbidities, and demographics. The effect of ICI on overall survival was quantified with Cox regression incorporating matching weights. Cumulative life‐years gained system‐wide were calculated from the difference in adjusted 5‐year restricted mean survival times. Results There were 27,322 patients in the ICI cohort and 69,801 patients in the historical control cohort. Among ICI patients, the most common cancer types were NSCLC (46%) and melanoma (10%). ICI demonstrated a large OS benefit in most cancer types with heterogeneity across cancer types (NSCLC: adjusted HR [aHR] 0.56, 95% confidence interval [CI] 0.54–0.58,p < 0.001; urothelial: aHR 0.91, 95% CI 0.83–1.01,p = 0.066). The relative benefit of ICI was stable across patient age, comorbidity, and self‐reported race subgroups. Across VHA, 15,859 life‐years gained were attributable to ICI within 5‐years of treatment, with NSCLC contributing the most life‐years gained. Conclusion We demonstrated substantial increase in survival due to ICIs across a national health system, including in patient subgroups poorly represented on clinical trials.

Oncology↗
DOE OSTI2022

Risk and incidence of head and neck cancers in veterans living with HIV and matched HIV‐negative veterans

Background Persons living with HIV/AIDS have a higher incidence of virus‐related and tobacco/alcohol‐related cancers. This study is the first to estimate the effect of HIV versus HIV‐negative veterans on the risk of head and neck squamous cell carcinoma incidence in a large retrospective cohort study. Methods The authors constructed a retrospective cohort study using patient data from 1999 to 2016 from the National Veterans Administration Corporate Data Warehouse and the VA Central Cancer Registry. This cohort study included 45,052 veterans living with HIV/AIDS and 162,486 HIV‐negative patients matched by age, sex, and index visit (i.e., HIV diagnosis date or clinic visit date). The age‐standardized incidence rates and estimated adjusted hazard ratios were calculated with a Cox proportional hazards regression for oropharyngeal and nonoropharyngeal head and neck cancer squamous cell carcinoma (HNSCC). The authors also abstracted human papillomavirus (HPV) status from oropharyngeal HNSCC diagnosed after 2010. Results Veterans living with HIV/AIDS (VLWH) have 1.71 (95% confidence interval [CI], 1.36, 2.14) times the risk of oropharyngeal cancer and 2.06 (95% CI, 1.76, 2.42) times the hazard of nonoropharyngeal cancer compared with HIV‐negative veterans. VLWH with oropharyngeal squamous cell carcinoma (OPSCC) were more likely to be HPV‐positive ( N = 30 [81.1%]) than the HIV‐negative veterans with OPSCC ( N = 50 [67.6%]), although this difference was not significant ( p = .135). For nonoropharyngeal cancer, the increased risk of oral cavity cancer among VLWH drove the increased risk. Conclusions The study results suggest that HIV may play a role in virally mediated and nonvirally mediated HNSCC. As the HIV prevalence rises in the United States due to better survival and the incidence of HPV‐positive oropharyngeal HNSCC increases, the interaction between HPV and HIV becomes increasingly relevant.

Mazul, Angela L.↗