SEARCH · Engineering Papers

Results for “likelihood ratio”

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 127 records · Page 7

All sources NASA NTRS DOE OSTI

DOE OSTI2020

The clustering of DESI-like luminous red galaxies using photometric redshifts

ABSTRACT We present measurements of the redshift-dependent clustering of a DESI-like luminous red galaxy (LRG) sample selected from the Legacy Survey imaging data set, and use the halo occupation distribution (HOD) framework to fit the clustering signal. The photometric LRG sample in this study contains 2.7 million objects over the redshift range of 0.4 < z < 0.9 over 5655 deg2. We have developed new photometric redshift (photo-z) estimates using the Legacy Survey DECam and WISE photometry, with σNMAD = 0.02 precision for LRGs. We compute the projected correlation function using new methods that maximize signal-to-noise ratio while incorporating redshift uncertainties. We present a novel algorithm for dividing irregular survey geometries into equal-area patches for jackknife resampling. For a five-parameter HOD model fit using the MultiDark halo catalogue, we find that there is little evolution in HOD parameters except at the highest redshifts. The inferred large-scale structure bias is largely consistent with constant clustering amplitude over time. In an appendix, we explore limitations of Markov chain Monte Carlo fitting using stochastic likelihood estimates resulting from applying HOD methods to N-body catalogues, and present a new technique for finding best-fitting parameters in this situation. Accompanying this paper, we have released the Photometric Redshifts for the Legacy Surveys catalogue of photo-z’s obtained by applying the methods used in this work to the full Legacy Survey Data Release 8 data set. This catalogue provides accurate photometric redshifts for objects with z < 21 over more than 16 000 deg2 of sky.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2020

Massive peatland carbon banks vulnerable to rising temperatures

Peatlands contain one-third of the world’s soil carbon (C). If destabilized, decomposition of this vast C bank could accelerate climate warming; however, the likelihood of this outcome remains unknown. Here, we examine peatland C stability through five years of whole-ecosystem warming and two years of elevated atmospheric carbon dioxide concentrations (eCO 2 ). Warming exponentially increased methane (CH 4 ) emissions and enhanced CH 4 production rates throughout the entire soil profile; although surface CH 4 production rates remain much greater than those at depth. Additionally, older deeper C sources played a larger role in decomposition following prolonged warming. Most troubling, decreases in CO 2 :CH 4 ratios in gas production, porewater concentrations, and emissions, indicate that the peatland is becoming more methanogenic with warming. We observed limited evidence of eCO 2 effects. Our results suggest that ecosystem responses are largely driven by surface peat, but that the vast C bank at depth in peatlands is responsive to prolonged warming.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Value-of-Information Analysis of a Fracture Prediction Method

A core-based fracture prediction method is used to illustrate a Value of Information (VOI) decision-analysis protocol to inform completion decisions in tight-gas sandstones. The ratio of late host-rock cement to available pore volume or degradation index uses petrographic observations of cement distributions in core (including sidewall cores) to predict whether nearby but unsampled fractures (widths > 0.5 to 1 mm) are sealed (non-conductive) or open (conductive). Measurements from four sandstone plays suggest that the index correctly predicts open vs. sealed fractures with an accuracy in excess of 80%. The value added is calculated using Bayesian inference in which accuracy of the index serves as the likelihood of the prior distribution of open fractures to assess the posterior probability that data represent a useful predictor of producibility. VOI of the prediction method is more than three times the cost to acquire the data. VOI is most sensitive to play-specific geologic and cost parameters including cost to drill, expected revenue from a successful well, cost of completion, cost of acquiring data for the index, and fracture probability distributions. The approach provides a way to value acquiring fracture data and points to a need for zone-specific production data in tight-gas sandstones.

02 PETROLEUM↗

DOE OSTI2023

Measurement of the production of a $W$ boson in association with a charmed hadron in $pp$ collisions at $\sqrt{s}$ = $\mathrm{13 TeV}$ with the ATLAS detector

The production of a $W$ boson in association with a single charm quark is studied using 140 fb –1 of $\sqrt{s}$ = $13$ $\mathrm{TeV}$ proton-proton collision data collected with the ATLAS detector at the Large Hadron Collider. The charm quark is tagged by the presence of a charmed hadron reconstructed with a secondary-vertex fit. The $W$ boson is reconstructed from the decay to either an electron or a muon and the missing transverse momentum present in the event. The charmed mesons reconstructed are D + → K – π + π + and D *+ → D 0 π + → (K – π + ) π + and the charge conjugate decays in the fiducial regions where $p$ T ($e,μ$) > 30 GeV, |$η(e,μ)$| < 2.5, $p$ T (D (*) ) > 8 GeV, and |$η$(D (*) )| < 2.2. The integrated and normalized differential cross sections as a function of the pseudorapidity of the lepton from the $W$ boson decay, and of the transverse momentum of the charmed hadron, are extracted from the data using a profile likelihood fit. The measured total fiducial cross sections are ${σ}_{fid}^{OS – SS}$ (W – + D + ) = 50.2 ± 0.2 ${(stat)}_{–2.3}^{+2.4}$(syst) pb, ${σ}_{fid}^{OS – SS}$ (W + + D – ) = 48.5 ± 0.2 ${(stat)}_{–2.2}^{+2.3}$(syst) pb, ${σ}_{fid}^{OS – SS}$ (W – + D *+ ) = 51.1 ± 0.4 ${(stat)}_{–1.8}^{+1.9}$(syst) pb, ${σ}_{fid}^{OS – SS}$ (W + + D *– ) = 50.0 ± 0.4 ${(stat)}_{–1.8}^{+1.9}$(syst) pb. Results are compared with the predictions of next-to-leading-order quantum chromodynamics calculations performed using state-of-the-art parton distribution functions. Additionally, the ratio of charm to anticharm production cross sections is studied to probe the $s$-$\bar{s}$ quark asymmetry. The ratio is found to be ${R}_{c}^{±}$= 0.971 ± 0.006 (stat) ± 0.011 (syst). The ratio and cross-section measurements are consistent with the predictions obtained with parton distribution function sets that have a symmetric $s$-$\bar{s}$ sea, indicating that any $s$-$\bar{s}$ asymmetry in the Bjorken-x region relevant for this measurement is small.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2022

A full-ring variable-aperture cadmium zinc telluride system for whole-body single photon emission computed tomography: realistic simulations with phantoms

Single photon emission computed tomography (SPECT) is an imaging modality that has demonstrated its utility in a number of clinical indications. Despite this progress, a high sensitivity, high spatial resolution, multi-tracer SPECT with a large field of view suitable for whole-body imaging of a broad range of radiotracers for theranostics is not available. Purpose We have designed a cadmium zinc telluride (CZT) variable-aperture full-ring SPECT scanner instrumented with a broad-energy tungsten collimator intended to fill this technological gap. The final purpose is to provide a multi-tracer solution for brain and whole-body imaging. Our static SPECT scanner breaks the paradigm of the standard dual- and triple-head rotational SPECT systems, utilizing a larger detector area in each scan increasing the sensitivity. We provide a demonstration of the performance of our design using a realistic model of our detector with simulated body-sized 99mTc phantoms. Methods We developed a realistic model of our detector by using a combination of a Geant4 Monte Carlo simulation and a CZT detector response model based on a finite element model. Our approach models the characteristic low-energy tail effect in CZT that noticeably affects the sensitivity and the quality of the scatter correction in CZT detectors. We implement a modified dual energy window scatter correction adapted to include the CZT low-energy tail effect. A dedicated correction is also developed to eliminate the undesirable truncation observed in images given the presence of detector edges and gaps between detectors, due to the non-rotational nature of our device. Corrections for the attenuation, detector response and the presence of collimators are also included. The images are reconstructed using the maximum-likelihood expectation-maximization algorithm implemented in the reconstruction open software STIR. Detector and reconstruction performance are characterized with a Derenzo phantom and a body-sized National Electrical Manufacturers Association (NEMA) Image Quality (IQ) phantom containing 99mTc. Results Our SPECT design can resolve 6.4mm rods in a Derenzo phantom and obtain a good image contrast with the IQ phantom. Explicit testing of the gap and edge correction is provided, showing an excellent performance in eliminating the image truncation artifacts. Our modified scatter correction shows no overestimation of the contrast-recovery ratio for our realistic CZT detector model, as opposed to the cases without correction and with a standard dual-energy window scatter correction. Conclusions In this paper, we further demonstrate the performance of our design for whole-body imaging purposes. This adds to our previous demonstration of improved qualitative and quantitative 99mTc imaging for brain perfusion and 123I imaging for dopamine transport with respect to state-of-the-art NaI dual-head cameras. We show that our design performs similarly to the VERITON SPECT from Spectrum Dynamics, a commercial full-ring CZT SPECT camera, with the potential advantage of the broader energy range of application given by our custom-design tungsten collimators. Here, our device combines high sensitivity and image resolution with a broad-energy imaging application for the purpose of clinical imaging and theranostics of emerging radionuclides.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

DOE OSTI2021

Data from: "Ecophysiological variation in two provenances of Pinus flexilis seedlings across an elevation gradient from forest to alpine"

This archive contains data used to support conclusions drawn in “Ecophysiological variation in two provenances of Pinus flexilis seedlings across an elevation gradient from forest to alpine”, by Reinhardt et al., 2011. Data were collected over one summer season in plots within the Alpine Treeline Warming Experiment (ATWE), before climate manipulations began. The experiment was located on Niwot Ridge, in the Front Range of the Colorado Rocky Mountains. This data package includes five comma-separated-values (.csv) files, five Microsoft Excel (.xlsx) files, one .pdf file, and two types of geospatial files: keyhole markup language (.kml), and ESRI shapefiles (.shp). .csv files can be opened using any simple text-editing software (such as Notepad and TextEdit), R, and Microsoft Excel. .xlsx files can only be opened using Microsoft Excel. The .pdf file can be opened using Adobe Acrobat Reader or any other compatible file viewing software. The .kml file can be opened using Google Earth and Google Maps, and shapefiles can be opened using any software compatible with the file type, such as ESRI’s ArcGIS suite and QGIS.Data archived contain gas exchange and plant physiology measurements, non-structural carbohydrate data, among others. Geospatial files are also provided for additional locational context. The files and their contents in this data package are summarized under "Data Summary" in the included Data User's Guide. All files (excluding geospatial) are available in both Microsoft Excel and in .csv format, and are indicated in the Data Summary list as well.-----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------Climate change is predicted to cause upward shifts in forest tree distributions, which will require seedling recruitment beyond current forest boundaries. However, predicting the likelihood of successful plant establishment beyond current species’ ranges under changing climate is complicated by the interaction of genetic and environmental controls on seedling establishment. To determine how genetics and climate may interact to affect seedling establishment, we transplanted recently germinated seedlings from high- and low-elevation provenances (HI and LO, respectively) of Pinus flexilis in common gardens arrayed along an elevation and canopy gradient from subalpine forest into the alpine zone and examined differences in physiology and morphology between provenances and among sites. Plant dry mass, projected leaf area and shoot:root ratios were 12–40% greater in LO compared with HI seedlings at each elevation. There were no significant changes in these variables among sites except for decreased dry mass of LO seedlings in the alpine site. Photosynthesis, carbon balance (photosynthesis/respiration) and conductance increased >2× with elevation for both provenances, and were 35–77% greater in LO seedlings compared with HI seedlings. There were no differences in dark-adapted chlorophyll fluorescence (Fv/Fm) among sites or between provenances. Our results suggest that for P. flexilis seedlings, provenances selected for above-ground growth may outperform those selected for stress resistance in the absence of harsh climatic conditions, even well above the species’ range limits in the alpine zone. This indicates that forest genetics may be important to understanding and managing species’ range adjustments due to climate change.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

Effects of Turbulence and Temperature Fluctuations on Knock Development in an Ethanol/Air Mixture

The effects of turbulence on knock development and intensity for a thermally inhomogeneous stoichiometric ethanol/air mixture at a representative end-gas autoignition condition in internal combustion engines are investigated using direct numerical simulations with a skeletal reaction mechanism. Two- and three-dimensional simulations are performed by varying the most energetic length scale of temperature, $l_T$ , and its relative ratio with the most energetic length scale of turbulence, $l_T/l_e$ , together with two different levels of the turbulent velocity fluctuation, $u'$ . It is found that $l_T$ / $l_e$ and the ratio of ignition delay time to eddy-turnover time, $\tau _{ig}/\tau _t$ , are the key parameters that control the detonation development. An increase in either $l_T$ or $l_e$ enhances the detonation propensity by allowing a longer run-up distance for the detonation development. The characteristic length scale of the temperature field, $l_T$ , is significantly modified by high turbulence intensity achieved by a large $l_e$ and $u'$ . The intense turbulence mixing effectively distributes the initial temperature field to broader scales to support the developing detonation waves, thereby increasing the likelihood of the detonation formation. On the contrary, high turbulence intensity with a short mixing time scale, achieved by a small $l_e$ and a large $u'$ , reduces the super-knock intensity attributed to the finer broken-up structures of detonation waves. Either $\tau _{ig}/\tau _t$ less than unity or $l_e = l_T$ even with a large $u'$ is found to have no significant effect on super-knock mitigation. Finally, high turbulent intensity may induce high-pressure spikes comparable to the von Neumann spike. Overall, increased temperature and pressure by combustion heating, noticeably after the peak of heat release rate, significantly enhance the collision and interaction of multiple emerging autoignition fronts near the ending combustion process, resulting in localized high-pressure spikes.

42 ENGINEERING↗

DOE OSTI2025

Impact of Geomagnetic Induced Current Neutral Blocking Devices on Distance Relays in Sub-transmission Networks with IBR

Geomagnetically induced currents (GICs) can flow through transmission lines during geomagnetic disturbances, such as solar flares or coronal mass ejections. These currents can cause problems like transformer saturation and equipment damage. The most common method of mitigating GICs involves installing GIC neutral blocking devices (NBDs) in transformer neutrals. However, the wide application of capacitive GIC blocking devices may have unintended adverse effects on other devices, such as distance protection relays. As the number of inverter-based resources being connected to the transmission and sub-transmission systems increases, the likelihood of a sub-transmission line protected by a distance relay connected to a transformer with NBDs is increasing. Therefore, distance relays fed by IBRs and transmission lines with GIC-NBDs must be studied. This paper studies the effect of GIC-NBDs on a 69kV sub-transmission line of various lengths fed by a 25 MVA IBR and synchronous source. This work focused on the behavior of the GIC-NBDs using the measured apparent phase-to-ground and phase-to-phase impedance calculated by the relay and the source impedance ratio (SIR) during various electrical faults.

Patel, Trupal R [Sandia National Laboratories (SNL↗

DOE OSTI2021

A demonstration of improved constraints on primordial gravitational waves with delensing

We present a constraint on the tensor-to-scalar ratio, $r$, derived from measurements of cosmic microwave background (CMB) polarization $B$-modes with "delensing,'' whereby the uncertainty on $r$ contributed by the sample variance of the gravitational lensing $B$-modes is reduced by cross-correlating against a lensing $B$-mode template. This template is constructed by combining an estimate of the polarized CMB with a tracer of the projected large-scale structure. The large-scale-structure tracer used is a map of the cosmic infrared background derived from Planck satellite data, while the polarized CMB map comes from a combination of South Pole Telescope, BICEP/Keck, and Planck data. We expand the BICEP/Keck likelihood analysis framework to accept a lensing template and apply it to the BICEP/Keck data set collected through 2014 using the same parametric foreground modelling as in the previous analysis. From simulations, we find that the uncertainty on $r$ is reduced by $\sim10\%$, from $\sigma(r)$= 0.024 to 0.022, which can be compared with a $\sim26\%$ reduction obtained when using a perfect lensing template. Applying the technique to the real data, the constraint on $r$ is improved from $r_{0.05} < 0.090$ to $r_{0.05} < 0.082$ (95% C.L.). Furthermore, this is the first demonstration of improvement in an $r$ constraint through delensing.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2026

PISCES two-detector covariance matrix fit for the NOvA Experiment

NOvA is a long-baseline neutrino oscillation experiment with two functionally identical detectors: a Near Detector (ND) at Fermilab, placed 1 km from the neutrino source, and a Far Detector (FD) located 810 km away from the ND in Minnesota. NOvA's primary physics goals are the precise measurements of neutrino oscillation parameters $\theta_{23}$ and $\Delta m^2_{32}$ , determine the neutrino mass ordering, and constrain the value of $\delta_{CP}$, via the study of muon neutrino to electron neutrino oscillation. In the standard NOvA three-flavor analysis, oscillation parameters are extracted using an extrapolation technique in which the ND data constrain the FD prediction through a ratio method. While this allows for systematic uncertainties sharing the same effects in both detectors to cancel, it remains an FD-only fit and does not fully leverage the constraining power of the high-statistics ND. This analysis proposes a simultaneous ND+FD fit using the PISCES method. PISCES (Parameter Inference with Systematic Covariance and Exact Statistics) is a framework designed to support complex configurations such as a joint ND+FD fit. This allows PISCES to take full advantage of the ND data to directly constrain systematic uncertainties across all samples. In PISCES, systematic uncertainties are encoded in a fractional covariance matrix, and statistical uncertainties are handled with a Poisson likelihood, making the approach well suited for low-statistics samples. For interpretability, we further use a Newton–Raphson + PCA method to recover per-systematic pulls from the covariance formulation. This poster presents the full PISCES joint ND+FD fit for the NOvA three-flavor analysis, describes its implementation and evaluates its performance through extensive robustness tests and fake data studies. It also provides a comparison between the PISCES joint ND+FD results and the standard NOvA extrapolation method.

Rajaoalisoa, Miriama [Cincinnati U.] (ORCID:000000↗

DOE OSTI2023

Analysis of Unpolarised p+p¿ Photoproduction with the GlueX Experiment

This thesis presents measurements of spin-density matrix elements in unpolarised p+p? photoproduction on a proton target. The dominant resonance contribution to the dipion system is the r(770) meson. Due to the large production cross section for this resonance, a valid comparison can be made between the obtained final results and r(770) spin-density matrix elements measured previously with other experiments. The measurement was performed over the 3:0 ? 11:6 GeV beam energy regime, which is a more extensive energy range than has ever been studied previously for the r(770). Results were obtained by analysing data from the GlueX experiment based at Jefferson Lab. Extended maximum likelihood fits were applied to extract three spin-density matrix elements using Markov chain Monte Carlo based parameter estimations. This was performed using various binning configurations to probe the energy, mass, and four-momentum transfer dependence of the determined physics observables. Spin-density matrix elements are shown to be consistent with the model of s-channel helicity conservation at low ?t. The effects of pomeron and f2 exchanges are clearly visible in the energy dependence of the measured observables. These observations provide valuable insights into the relative strengths of both processes as a function of the photon energy, and may enable theorists to disentangle the f2=P coupling ratio. Spin-density matrix elements are seen to be highly dependent on the reconstructed resonance mass. This observation is likely to be a result of non-resonant S-wave background processes, and emphasises the need for a more detailed model of the p+p? angular distribution that considers all of the competing angular momentum components that contribute to the measured final state. The statistical precision of measurements performed for this thesis surpass what was achievable in previous studies of the r(770) by several orders of magnitude. Studies of the energy and four-momentum transfer dependence, and insights into the effects of non r(770) background contributions provide valuable input for production models. This will help inform the choice of wave-sets used for partial wave analyses, supporting GlueX in its search for exotic hybrid meson states.

Fitches, James↗

DOE OSTI2022

Diabetes and hypertension are related to amyloid-beta burden in the population-based Rotterdam Study

Abstract Higher vascular disease burden increases the likelihood of developing dementia, including Alzheimer’s disease. Better understanding the association between vascular risk factors and Alzheimer’s disease pathology at the predementia stage is critical for developing effective strategies to delay cognitive decline. In this work, we estimated the impact of six vascular risk factors on the presence and severity of in vivo measured brain amyloid-beta (Aβ) plaques in participants from the population-based Rotterdam Study. Vascular risk factors (hypertension, hypercholesterolaemia, diabetes, obesity, physical inactivity and smoking) were assessed 13 (2004–2008) and 7 years (2009–2014) prior to 18F-florbetaben PET (2018–2021) in 635 dementia-free participants. Vascular risk factors were associated with binary amyloid PET status or continuous PET readouts (standard uptake value ratios, SUVrs) using logistic and linear regression models, respectively, adjusted for age, sex, education, APOE4 risk allele count and time between vascular risk and PET assessment. Participants’ mean age at time of amyloid PET was 69 years (range: 60–90), 325 (51.2%) were women and 190 (29.9%) carried at least one APOE4 risk allele. The adjusted prevalence estimates of an amyloid-positive PET status markedly increased with age [12.8% (95% CI 11.6; 14) in 60–69 years versus 35% (36; 40.8) in 80–89 years age groups] and APOE4 allele count [9.7% (8.8; 10.6) in non-carriers versus 38.4% (36; 40.8) to 60.4% (54; 66.8) in carriers of one or two risk allele(s)]. Diabetes 7 years prior to PET assessment was associated with a higher risk of a positive amyloid status [odds ratio (95% CI) = 3.68 (1.76; 7.61), P < 0.001] and higher standard uptake value ratios, indicating more severe Aβ pathology [standardized beta = 0.40 (0.17; 0.64), P = 0.001]. Hypertension was associated with higher SUVr values in APOE4 carriers (mean SUVr difference of 0.09), but not in non-carriers (mean SUVr difference 0.02; P = 0.005). In contrast, hypercholesterolaemia was related to lower SUVr values in APOE4 carriers (mean SUVr difference −0.06), but not in non-carriers (mean SUVr difference 0.02). Obesity, physical inactivity and smoking were not related to amyloid PET measures. The current findings suggest a contribution of diabetes, hypertension and hypercholesterolaemia to the pathophysiology of Alzheimer’s disease in a general population of older non-demented adults. As these conditions respond well to lifestyle modification and drug treatment, further research should focus on the preventative effect of early risk management on the development of Alzheimer’s disease neuropathology.

van Arendonk, Joyce (ORCID:0000000172941633)↗

DOE OSTI2024

Observation of Significant Flavor-SU(3) Breaking in the Kaon Wave Function at 12 < Q 2 < 25 GeV 2 and Discovery of the Charmless Decay ψ ( 3770 ) → K S 0 K L 0

We present cross sections for the reaction e + e − → K S 0 K L 0 at center-of-mass energies ranging from 3.51 to 4.95 GeV using data samples collected in the BESIII experiment, corresponding to a total integrated luminosity of 26.5 fb − 1 . The ratio of neutral-to-charged kaon form factors at large momentum transfers ( 12 < Q 2 < 25 GeV 2 ) is determined to be 0.21 ± 0.01 , which indicates a small but significant effect of flavor-SU(3) breaking in the kaon wave function, and, consequently, excludes the possibility that flavor-SU(3) breaking is the primary reason for the strong experimental violation of the pQCD prediction | F ( π ± ) | / | F ( K ± ) | = f π 2 / f K 2 , where F ( π ± ) and F ( K ± ) are the form factors, and f π and f K are the decay constants of charged pions and kaons, respectively. We also observe a significant signal for the charmless decay ψ ( 3770 ) → K S 0 K L 0 for the first time. Within a 1 σ contour of the likelihood value, the branching fraction for ψ ( 3770 ) → K S 0 K L 0 is determined to be B = ( 2.6 3 − 1.59 + 1.40 ) × 10 − 5 , and the relative phase between the continuum and ψ ( 3770 ) amplitudes is ϕ = ( − 0.3 9 − 0.10 + 0.05 ) π . The branching fraction is in good agreement with the S - and D -wave charmonia mixing scheme proposed in the interpretation of the “ ρ π puzzle” between J / ψ and ψ ( 3686 ) decays. Published by the American Physical Society 2024

Physics↗

DOE OSTI2024

CMB-S4: Iterative Internal Delensing and r Constraints

Abstract The tightest constraints on the tensor-to-scalar ratio r can only be obtained after removing a substantial fraction of the lensing B -mode sample variance. The planned Cosmic Microwave Background (CMB)-S4 experiment ( cmb-s4.org ) will remove the lensing B -mode signal internally by reconstructing the gravitational lenses from high-resolution observations. We document here a first lensing reconstruction pipeline able to achieve this optimally for arbitrary sky coverage. We make it part of a map-based framework to test CMB-S4 delensing performance and its constraining power on r , including inhomogeneous noise and two non-Gaussian Galactic polarized foreground models. The framework performs component separation of the high-resolution maps, followed by the construction of lensing B -mode templates, which are then included in a parametric small-aperture map cross-spectra-based likelihood for r . We find that the lensing reconstruction and framework achieve the expected performance, compatible with the target σ ( r ) ≃ 5 · 10 −4 in the absence of a tensor signal, after an effective removal of 92%–93% of the lensing B -mode variance, depending on the simulation set. The code for the lensing reconstruction can also be used for cross-correlation studies with large-scale structures, lensing spectrum reconstruction, cluster lensing, or other CMB lensing-related purposes. As part of our tests, we also demonstrate the joint optimal reconstruction of the lensing potential with the lensing curl potential mode at second order in the density fluctuations.

Astronomy & Astrophysics↗

DOE OSTI2022

Monte Carlo Simulation and Reconstruction: Assessment of Myocardial Perfusion Imaging of Tracer Dynamics With Cardiac Motion Due to Deformation and Respiration Using Gamma Camera With Continuous Acquisition

Purpose: Myocardial perfusion imaging (MPI) with single photon emission computed tomography (SPECT) is routinely used for stress testing in nuclear medicine. Recently, our group extended its potential going from 3D visual qualitative image analysis to 4D spatiotemporal reconstruction of dynamically acquired data to capture the time variation of the radiotracer concentration and the estimated myocardial blood flow (MBF) and coronary flow reserve (CFR). However, the quality of reconstructed image is compromised due to cardiac deformation and respiration. The work presented here develops an algorithm that reconstructs the dynamic sequence of separate respiratory and cardiac phases and evaluates the algorithm with data simulated with a Monte Carlo simulation for the continuous image acquisition and processing with a slowly rotating SPECT camera. Methods: A clinically realistic Monte Carlo (MC) simulation is developed using the 4D Extended Cardiac Torso (XCAT) digital phantom with respiratory and cardiac motion to model continuous data acquisition of dynamic cardiac SPECT with slowly rotating gamma cameras by incorporating deformation and displacement of the myocardium due to cardiac and respiratory motion. We extended our previously developed 4D maximum-likelihood expectation-maximization (MLEM) reconstruction algorithm for a data set binned from a continuous list mode (LM) simulation with cardiac and respiratory information. Our spatiotemporal image reconstruction uses splines to explicitly model the temporal change of the tracer for each cardiac and respiratory gate that delineates the myocardial spatial position as the tracer washes in and out. Unlike in a fully list-mode data acquisition and reconstruction the accumulated photons are binned over a specific but very short time interval corresponding to each cardiac and respiratory gate. Reconstruction results are presented showing the dynamics of the tracer in the myocardium as it continuously deforms. These results are then compared with the conventional 4D spatiotemporal reconstruction method that models only the temporal changes of the tracer activity. Mean Stabilized Activity (MSA), signal to noise ratio (SNR) and Bias for the myocardium activities for three different target-to-background ratios (TBRs) are evaluated. Dynamic quantitative indices such as wash-in (K1) and wash-out (k2) rates at each gate were also estimated. Results: The MSA and SNR are higher with higher TBRs while biases were improved with higher TBRs to less than 10%. The correlation between exhalation-inhalation sequence with the ground truth during respiratory cycle was excellent. Our reconstruction method showed better resolved myocardial walls during diastole to systole as compared to the ungated 4D image. Estimated values of K1 and k2 were also consistent with the ground truth. Conclusion: The continuous image acquisition for dynamic scan using conventional two-head gamma cameras can provide valuable information for MPI. Our study demonstrated the viability of using a continuous image acquisition method on a widely used clinical two-head SPECT system. Our reconstruction method showed better resolved myocardial walls during diastole to systole as compared to the ungated 4D image. Precise implementation of reconstruction algorithms, better segmentation techniques by generating images of different tissue types and background activity would improve the feasibility of the method in real clinical environment.

60 APPLIED LIFE SCIENCES↗

DOE OSTI2026

A dendritic strontium river isoscape for fisheries applications in the Sacramento River basin, California, USA

Objective Understanding the origins and movements of fish is fundamental to effective conservation and fisheries management. Strontium isotope ratios ( 87 Sr/ 86 Sr) in otoliths provide a powerful tracer of natal origin and migratory pathways. However, existing 87 Sr/ 86 Sr isoscapes for the Sacramento River basin, an ecosystem that supports ecologically and economically important salmon populations, rely on discrete classification approaches that overlook unsampled habitats and do not incorporate spatial uncertainty. Our objective was to develop a continuous, network-explicit 87 Sr/ 86 Sr isoscape with quantified uncertainty to fill in data gaps and enable probabilistic assignments of fish origin and movement. Methods We used river water 87 Sr/ 86 Sr data from 106 sites (1997–2021) to develop spatial stream network models that use dendritic connectivity and watershed characteristics (lithology, bedrock age, and land cover) to predict river water 87 Sr/ 86 Sr throughout the basin. Models were fitted using maximum and restricted likelihood and were evaluated via Akaike’s information criterion and leave-one-out cross validation. We produced both historical (pre-dam) and present-day (below-dam) isoscapes, delineated uncertainty-informed isotopic ranges using k -means clustering, and applied a proof-of-concept Bayesian assignment to estimate natal origins and early rearing habitats for two endangered winter-run Chinook Salmon Oncorhynchus tshawytscha. Results Cross validation indicated strong performance of the 87 Sr/ 86 Sr model (leave-one-out cross validation: R 2 = 0.91; root mean square error = 0.0005). Uncertainty-informed clustering identified 19 isotopic “suites” (reaches with indistinguishable 87 Sr/ 86 Sr values) in present-day anadromous habitats and 25 suites in the historical network. Example natal and early rearing assignments included predictions that challenged expectations for juvenile salmon migration based on predicted river 87 Sr/ 86 Sr compositions. Conclusions This study developed a continuous, network-explicit 87 Sr/ 86 Sr isoscape that integrates existing river data to predict 87 Sr/ 86 Sr in unsampled reaches and the likely achievable range and resolution of otolith-based origin and life history inference. The resulting river isoscape provides a valuable tool to predict salmon movements and identify habitats supporting their survival and growth that otherwise might remain undetected. Coupling these predictions with complementary approaches that ground-truth juvenile presence (e.g., targeted fish surveys) represents an important step toward science-informed restoration and management of critical habitats throughout the Sacramento River basin.

Environmental sciences↗

DOE OSTI2022

Combining Planck and SPT Cluster Catalogs: Cosmological Analysis and Impact on the Planck Scaling Relation Calibration

We provide the first combined cosmological analysis of the South Pole Telescope (SPT) and Planck cluster catalogs. The aim is to provide an independent calibration for Planck scaling relations, exploiting the cosmological constraining power of the SPT-SZ cluster catalog and its dedicated weak lensing (WL) and X-ray follow-up observations. We build a new version of the Planck cluster likelihood. In the νΛ CDM scenario, focusing on the mass slope and mass bias of Planck scaling relations, we find ${\alpha }_{\mathrm{SZ}}={1.49}_{-0.10}^{+0.07}$ and ${\left(1-b\right)}_{\mathrm{SZ}}={0.69}_{-0.14}^{+0.07}$, respectively. The results for the mass slope show a ~4 σ departure from the self-similar evolution, α SZ ~ 1.8. This shift is mainly driven by the matter density value preferred by SPT data, Ωm = 0.30 ± 0.03, lower than the one obtained by Planck data alone, ${{\rm{\Omega }}}_{m}={0.37}_{-0.06}^{+0.02}$. The mass bias constraints are consistent both with outcomes of hydrodynamical simulations and external WL calibrations, (1 – b) ~ 0.8, and with results required by the Planck cosmic microwave background cosmology, (1 – b) ~ 0.6. From this analysis, we obtain a new catalog of Planck cluster masses M 500 . We estimate the ratio between the published Planck M SZ masses and our derived masses M 500 , as a "measured mass bias," ${\left(1-b\right)}_{M}$. We analyze the mass, redshift, and detection noise dependence of ${\left(1-b\right)}_{M}$, finding an increasing trend toward high redshift and low mass. These results mimic the effect of departure from self-similarity in cluster evolution, showing different dependencies for the low-mass, high-mass, low-z, and high-z regimes.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Framework for Quantitative Evaluation of Resilience Solutions: An Approach to Determine the Value of Resilience for a Particular Site

The paper provides the approach to providing a benefit cost analysis of energy and water alternatives to provide resilience to extreme events. The approach estimates the costs and returns of providing greater resilience of water and energy infrastructure. Extreme events are defined as high impact, low-frequency events such as, but not limited to, hurricanes, floods, storm surges and earthquakes. The provides justification for hardening water and energy infrastructure. Resilience is defined as “the ability to prepare for and to withstand an extreme event with little or no damage, or to recover more quickly from an extreme event.” The approach can be summarized as follows. The approach requires the development of a baseline with which to compare alternatives. The baseline is used to evaluate the baseline’s resilience to hazards through the probability of the hazard(s), the likelihood of damage from that the hazard through a vulnerability analysis, and the consequence to calculate a cost of the damage. The approach then evaluates proposed mitigation alternatives that would improve the resilience of the system. Each alternative is evaluated based on probability of the hazard, probability of vulnerability and consequence to determine the reduced damage that each alternative presents. The approach includes any monetary and non-monetary benefits that can quantified for each of the alternatives. Non-quantifiable benefits are evaluated based on the relative importance of each alternative to the criteria used to determine how well the alternative meets the goals and objectives of the site/facility. Then, a life cycle cost analysis should be conducted for the baseline and alternatives. Finally, the results of the life cycle analysis should be presented in a decision matrix with cost, net present value, benefit/cost ratios, and any non-monetary criteria ranked to show how well the alternatives met the criteria, weighted with the decision maker’s weights and the results presented.

54 ENVIRONMENTAL SCIENCES↗