Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “bias correction”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Comparative Analysis of TRGBs (CATs) from Unsupervised, Multi-halo-field Measurements: Contrast is Key

The tip of the red giant branch (TRGB) is an apparent discontinuity of the luminosity function (LF) due to the end of the red giant evolutionary phase and is used to measure distances in the local universe. In practice, tip localization via edge detection response (EDR) relies on several methods applied on a case-by-case basis. It is hard to evaluate how individual choices affect a distance estimation using only a single host field while also avoiding confirmation bias. To devise a standardized approach, we compare unsupervised, algorithmic analyses of the TRGB in multiple halo fields per galaxy. We first optimize methods for the lowest field-to-field dispersion, including spatial filtering, smoothing, and weighting of LF, color band selection, and tip selection based on the number of likely RGB stars and the ratio of stars below versus above the tip (R). We find R, which we call the tip contrast, to be the most important indicator of the quality of EDR measurements; higher R selection can decrease field-to-field dispersion. Further, since R is found to correlate with the age or metallicity of the stellar population based on theoretical modeling, it might result in a displacement of the detected tip magnitude. We find a tip-contrast relation with a slope of -0.023 ± 0.0046 mag/ratio, an ~5σ result that can be used to correct these variations in the detections. When using TRGB to establish a distance ladder, consistent TRGB standardization using tip-contrast relation across rungs is vital to make robust cosmological measurements.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Correcting correlation functions for redshift-dependent interloper contamination

ABSTRACT The construction of catalogues of a particular type of galaxy can be complicated by interlopers contaminating the sample. In spectroscopic galaxy surveys this can be due to the misclassification of an emission line; for example in the Hobby–Eberly Telescope Dark Energy Experiment (HETDEX) low-redshift [O ii] emitters may make up a few per cent of the observed Ly α emitter (LAE) sample. The presence of contaminants affects the measured correlation functions and power spectra. Previous attempts to deal with this using the cross-correlation function have assumed sources at a fixed redshift, or not modelled evolution within the adopted redshift bins. However, in spectroscopic surveys like HETDEX, where the contamination fraction is likely to be redshift dependent, the observed clustering of misclassified sources will appear to evolve strongly due to projection effects, even if their true clustering does not. We present a practical method for accounting for the presence of contaminants with redshift-dependent contamination fractions and projected clustering. We show using mock catalogues that our method, unlike existing approaches, yields unbiased clustering measurements from the upcoming HETDEX survey in scenarios with redshift-dependent contamination fractions within the redshift bins used. We show our method returns autocorrelation functions with systematic biases much smaller than the statistical noise for samples with at least as high as 7 per cent contamination. We also present and test a method for fitting for the redshift-dependent interloper fraction using the LAE–[O ii] galaxy cross-correlation function, which gives less biased results than assuming a single interloper fraction for the whole sample.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Forward modeling fluctuations in the DESI LRGs target sample using image simulations

We use the forward modeling pipeline, Obiwan, to study the imaging systematics of the Luminous Red Galaxies (LRGs) targeted by the Dark Energy Spectroscopic Instrument (DESI). Imaging systematics refers to the false fluctuation of galaxy densities due to varying observing conditions and astrophysical foregrounds corresponding to the imaging surveys from which DESI LRG target galaxies are selected. We update the Obiwan pipeline, which we previously developed to simulate the optical images used to target DESI data, to further simulate WISE images in the infrared. This addition allows simulating the DESI LRGs sample, which utilizes WISE data in the target selection. Deep DESI imaging data combined with a method to account for biases in their shapes is used to define a truth sample of potential LRG targets. We inject these data evenly throughout the DESI Legacy Imaging Survey footprint at declinations between -30 and 32.375 degrees. We simulate a total of 15 million galaxies to obtain a simulated LRG sample (Obiwan LRGs) that predicts the variations in target density due to imaging properties. We find that the simulations predict the trends with depth observed in the data, including how they depend on the intrinsic brightness of the galaxies. We observe that faint LRGs are the main contributing source of the imaging systematics trend induced by depth. We also find significant trends in the data against Galactic extinction that are not predicted by Obiwan. These trends depend strongly on the particular map of Galactic extinction chosen to test against, implying systematic contamination in the Galactic extinction maps is a likely root cause (e.g., Cosmic-Infrared Background, dust temperature correction). We additionally observe a morphological change of the DESI LRGs population evidenced by a correlation between OII emission line average intensity and the size of the z-band PSF. This effect most likely results from uncertainties in background subtraction. The detailed findings we present should be used to guide any observational systematics mitigation treatment for the clustering of the DESI LRGs sample.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Reassessing long-standing meteorological records: an example using the national hottest day in Ireland

This analysis highlights the potential value in reanalysing early national meteorological records from around the world. These were oftentimes measured via techniques that preceded standardisation of instrumentation and methods of observation and thus could be subject to considerable biases and uncertainties. This analysis uses the techniques pioneered by WMO record assessment teams. The highest currently recognised air temperature (33.3 °C) ever recorded in the Republic of Ireland was logged at Kilkenny Castle in 1887. The original observational record however no longer exists. Given that Ireland is now the only country in Europe to have a national heat record that was set in the 19th century, a reassessment of the verity of this record is both timely and valuable. The present analysis undertakes a fundamental reassessment of the plausibility of the 1887 temperature record using methods similar to those used to assess various weather extremes under WMO auspices over recent years. Specifically, we undertake an inter-station reassessment using sparse available records and make recourse to the new and improved 20CRv3 sparse-input reanalysis product. Neither surrounding available stations nor the reanalysis offer substantive support for the Kilkenny record of 33.3 °C being correct. Moreover, recent data rescue efforts have uncovered several earlier extreme values, one of which exceeds the Kilkenny value (33.5 °C on 16 July 1876 recorded at the Phoenix Park). However, the sparsity of early observational networks, a distinct lack of synoptic support from 20CRv3 for many of the extreme heat values, and the fact that these measurements were obtained using non-standard exposures lead us to conclude that there is grossly insufficient evidence to support any of these 19th century extremes as robust national heat record candidates. Data from the early 20th century onwards benefit from a denser network of stations undertaking measurements in a more standardised manner, many under the direct auspices of Met Éireann and its predecessors, adhering to WMO guidance and protocols. This enables more robust cross-checking of records. We argue that the Met Éireann-recognised 20th century heat record from Boora in 1976 is verified as the most plausible robust national temperature record based upon the synoptic situation and comparisons with nearby neighbouring stations. This measurement of 32.5°C thus likely constitutes the highest reliably recorded temperature measurement in the Republic of Ireland. Ultimately, the formal decision on any reassessment and reassignment of the national record rests with the national meteorological service, Met Éireann.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Fast External Quantum Efficiency Measurements of Photovoltaic Modules

The external quantum efficiency (EQE) or associated spectral response is widely used to evaluate the performance of photovoltaics, light emitting diodes, photodetectors, photodiodes, semiconductor lasers, laser-induced refrigeration of solids etc. For photovoltaic (PV) devices, The EQE is used to determine the spectral mismatch correction that is used to convert a performance measurement from a non-ideal simulator spectrum to a standard spectrum in PV device calibrations. Moreover, EQE can provide important insight into the physics of a PV device. For instance, the EQE shape can be used to infer mechanisms limiting the performance of the device. However, EQE measurements of full-size PV modules remain a challenge in the PV community due to rare availability of reliable module QE measurement tools worldwide and the size and complexity of module architecture that comprises several cells connected in strings (series or/and parallel). Here, we use a customized nondestructive tool made of high-power light emitting diodes to perform high-throughput EQE measurements on PV modules. Circuit simulations show that the measurement system can safely be used to measure both crystalline silicon and thin film-based PV modules with series connected cells without the risk of damage through reverse bias. In addition, EQE mapping combined with electroluminescence show a strong correlation and provide a way to determine both spatial non-uniformity of the spectral response and defects/degradation effects across the module.

PV module↗
DOE OSTI2022

Fabrication of Black Body Grids by Thick Film Printing for Quantitative Neutron Imaging

Neutron imaging offers deep penetration through many high-Z materials while also having high sensitivity to certain low-Z isotopes such as 1H, 6Li, and 10B. This unique combination of properties has made neutron imaging an attractive tool for a wide range of material science and engineering applications. However, measurements made by neutron imaging or tomography are generally qualitative in nature due to the inability of detectors to discriminate between neutrons which have been transmitted through the sample and neutrons which are scattered by the sample or within the detector. Recent works have demonstrated that deploying a grid of small black bodies (BBs) in front of the sample can allow for the scattered neutrons to be measured at the BB locations and subsequently subtracted from the total measured intensity to yield a quantitative transmission measurement. While this method can be very effective, factors such as the scale and composition of the sample, the beam divergence, and the resolution and construction of the detector may require optimization of the grid design to remove all measurement biases within a given experimental setup. Therefore, it is desirable to have a method by which BB grids may be rapidly and inexpensively produced such that they can easily be tailored to specific applications. In this work, we present a method for fabricating BB patterns by thick film printing of Gd2O3 and evaluate the performance with variation in feature size and number of print layers with cold and thermal neutrons.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Parameter inference with non-linear galaxy clustering: accounting for theoretical uncertainties

ABSTRACT We implement euclidemulator (version 1), an emulator for the non-linear correction of the matter power spectrum, into the Markov chain Monte Carlo forecasting code montepython. We compare the performance of halofit, hmcode, and euclidemulator1, both at the level of power spectrum prediction and at the level of posterior probability distributions of the cosmological parameters, for different cosmological models and different galaxy power spectrum wavenumber cut-offs. We confirm that the choice of the power spectrum predictor has a non-negligible effect on the computed sensitivities when doing cosmological parameter forecasting, even for a conservative wavenumber cut-off of $0.2\, h\, {\rm Mpc}^{-1}$. We find that euclidemulator1 is on average up to 17 per cent more sensitive to the cosmological parameters than the other two codes, with the most significant improvements being for the Hubble parameter of up to 42 per cent and the equation of state of dark energy of up to 26 per cent, depending on the case. In addition, we point out that the choice of the power spectrum predictor contributes to the risk of computing a significantly biased mean cosmology when doing parameter estimations. For the four tested scenarios we find biases, averaged over the cosmological parameters, of between 0.5σ and 2σ (from below 1σ up to 6σ for individual parameters). This paper provides a proof of concept that this risk can be mitigated by taking a well-tailored theoretical uncertainty into account as this allows to reduce the bias by a factor of 2 to 5, depending on the case under consideration, while keeping posterior credibility contours small: the standard deviations are amplified by a factor of ≤1.4 in all cases.

Astronomy & Astrophysics↗
DOE OSTI2022

BASS. XXV. DR2 Broad-line-based Black Hole Mass Estimates and Biases from Obscuration

We present measurements of broad emission lines and virial estimates of supermassive black hole masses (M BH ) for a large sample of ultrahard X-ray-selected active galactic nuclei (AGNs) as part of the second data release of the BAT AGN Spectroscopic Survey (BASS/DR2). Our catalog includes M BH estimates for a total of 689 AGNs, determined from the Hα, Hβ, Mg II λ2798, and/or C IV λ1549 broad emission lines. The core sample includes a total of 512 AGNs drawn from the 70 month Swift/BAT all-sky catalog. We also provide measurements for 177 additional AGNs that are drawn from deeper Swift/BAT survey data. We study the links between M BH estimates and line-of-sight obscuration measured from X-ray spectral analysis. We find that broad Hα emission lines in obscured AGNs ($\mathrm{log}({N}_{{\rm{H}}}/{\mathrm{cm}}^{-2})\gt 22.0$) are on average a factor of ${8.0}_{-2.4}^{+4.1}$ weaker relative to ultrahard X-ray emission and about ${35}_{-12}^{\,+7}$% narrower than those in unobscured sources (i.e., $\mathrm{log}({N}_{{\rm{H}}}/{\mathrm{cm}}^{-2})\lt 21.5$). This indicates that the innermost part of the broad-line region is preferentially absorbed. Consequently, current single-epoch M BH prescriptions result in severely underestimated (>1 dex) masses for Type 1.9 sources (AGNs with broad Hα but no broad Hβ) and/or sources with $\mathrm{log}({N}_{{\rm{H}}}/{\mathrm{cm}}^{-2})\gtrsim 22.0$. We provide simple multiplicative corrections for the observed luminosity and width of the broad Hα component (L[bHα] and FWHM[bHα]) in such sources to account for this effect and to (partially) remedy M BH estimates for Type 1.9 objects. As a key ingredient of BASS/DR2, our work provides the community with the data needed to further study powerful AGNs in the low-redshift universe.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

A Comparison of Pre‐Construction and Operational Wake Loss Estimates for Land‐Based Wind Plants

The overall bias between pre‐construction energy yield assessment (EYA) estimates of wind plant energy production and the achieved operational production is improving in the wind industry, but uncertainty remains high for individual wind plants. Wake effects within wind plants are one of the largest sources of energy loss considered in the EYA process, and previous work shows wake loss estimates to be a major source of disagreement among wind energy consultants who perform EYAs. To better understand the accuracy of wake loss predictions, we compare overall operational wake loss estimates based on supervisory control and data acquisition data to pre‐construction estimates provided by six wind energy consultants for five land‐based wind plants in North America. By augmenting existing approaches for quantifying operational wake losses, we estimate wake losses during the period of record for which operational data are available as well as the expected long‐term wake losses, based on historical reanalysis weather data, to which the EYA estimates are compared. To account for power variations at different turbine locations caused by terrain‐induced wind resource heterogeneity, we correct the operational wake loss estimates using predicted freestream wind speed variations from the Wind Systems Engineering Reynolds‐averaged Navier–Stokes (RANS) tool. We identify long‐term corrected operational wake losses between 1.9% and 6.4% for the five plants, with a mean loss of 4%. For the project deemed most acceptable for operational wake loss assessment, which is located in the simplest terrain and isolated from neighboring plants, the mean EYA wake loss estimate is within 0.7 percentage points of the operational value of 6.4%. For most of the remaining plants, results suggest that wake losses are generally overpredicted by 2.6–6.3 percentage points. However, operational wake losses may be underestimated for many of these projects because of spatial wind resource variations not captured by the RANS model, external wake effects that are unaccounted for in the estimation process, and wind plant blockage effects. To better understand factors that contribute to the observed wake losses, we investigate operational wake losses as a function of wind direction and wind speed. As expected, wake losses are generally concentrated near wind directions that are aligned with rows of closely spaced turbines and at below‐rated wind speeds; however, for some projects, the energy produced by the wind plant exceeds the estimated potential energy of the plant without wake interactions for certain wind directions and wind speeds, suggesting inaccurate assumptions in the wake loss estimation method for those plants. Lastly, we compare predicted and operational wake losses for individual wind turbines, finding that even when overall wake losses are predicted accurately, large uncertainty exists at the turbine level.

17 WIND ENERGY↗
DOE OSTI2023

Measurement of Plutonium Ratios for CRM 137A Recertification

This report documents the analytical methodologies used by the Chemical and Isotopics Mass Spectrometry (CIMS) Group at the US Department of Energy’s Oak Ridge National Laboratory (ORNL) for plutonium isotopic ratio characterizations on purified ~1 mg units prepared by Lawrence Livermore National Laboratory (LLNL) from a mother solution created from an existing 250 mg unit of certified reference material (CRM) 137. The analytical processes described are designed to maintain traceability of the final certified attributes and to eliminate sources of systematic bias for the measurements described. The measurands for certification are the plutonium isotope ratios: 238 Pu/ 239 Pu, 240 Pu/ 239 Pu, 241 Pu/ 239 Pu, and 242 Pu/ 239 Pu. The preparation and analytical process included the following activities: Preparation of reference materials to be used as calibrants to establish traceability to the SI unit of mass for the isotopic analyses and quality controls. For the plutonium isotopic analyses, CRM 128 (a CRM comprising equal atoms of 239 Pu and 242 Pu) was used as the primary calibrant; Each CRM 137A unit was reconstituted using dilute high-purity HNO3. From each, seven replicates were distributed for the isotopic characterizations by thermal ionization mass spectrometry (TIMS) using the total evaporation (TE) method for isotope ratio measurements and MC-ICP-MS using a dynamic peak hooping method for secondary isotope ratio measurements; Due to time delay in the analyses from the last purification of the original CRM 137 mother solution by LLNL, an aliquot of the three provided plutonium units was purified using Eichrom TEVA separation methods to remove americium in-growth and uranium impurities; Minimally corrected data were reported from each analysis using reporting templates provided by the New Brunswick Laboratory Program Office (NBL PO) for the specific methods used for isotopic measurements. In addition, a copy of the instrument exports, Excel templates for isotopic data corrections, and the uncertainty budgets completed with the GUM Workbench software package were submitted to NBL PO.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2024

Tunneling magnetoresistance in MgO tunnel junctions with Fe-based leads in empirically corrected density functional theory

The minority-spin Fe/MgO interface states are at the Fermi level in density functional theory (DFT), but experimental evidence and GW calculations place them slightly higher in energy. This small shift can strongly influence tunneling magnetoresistance (TMR) in junctions with a thin MgO barrier and its dependence on the concentration of Co in the electrodes. Here, in this study, an empirical potential correction to DFT is introduced to shift the interface states up to match the tunnel spectroscopy data. With this shift, TMR in Fe/MgO/Fe junctions exceeds 800% and 3000% at 3 and 4 monolayers (ML) of MgO, respectively. We further consider the effect of alloying of the Fe electrodes with up to 30% Co or 10% V, treating them in the coherent potential approximation (CPA). Alloying with Co broadens the interface states and brings a large incoherent minority-spin spectral weight to the Fermi level. Alloying with V brings the minority-spin resonant states close to the Fermi level. However, in both cases the minority-spin spectral weight at the Fermi level resides primarily at the periphery of the Brillouin zone, which is favorable for spin filtering. Using convolutions of k II -resolved barrier densities of states calculated in CPA, it is found that TMR is strongly reduced by alloying with Co or V but still remains above 500% at 4 ML of MgO up to 30% of Co or 5% V. At 5 ML, the TMR increases above 1000% in all systems considered. However, while TMR declines sharply with increasing bias up to 0.2 eV in the tunnel junctions with pure Fe leads, it remains almost constant up to 0.5 eV if leads are alloyed with Co.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

PV Reference Cells for Outdoor Use: An Investigation of Calibration Factors

Reference cells are widely used in the photovoltaic (PV) industry to measure irradiance. For field applications and outdoor use, a variety of products are on the market, and they are often perceived as a low-cost alternative to thermopile radiometers. But reference cell characteristics differ substantially from the latter, therefore measurements made by one category of instrument cannot be substituted for or directly compared with the other. The National Renewable Energy Laboratory's (NREL's) Solar Radiation Research Laboratory is currently in a multiyear effort to develop guidance and recommendations for the design and use of outdoor reference cells, with a view to reducing inconsistencies and measurement uncertainty. The core of this effort is the long-term deployment of 39 units of 10 distinct types made by 6 manufacturers. This report compares the calibration factors provided by manufacturers - which most users rely upon - to those measured at a high level of accuracy by NREL's Cell Lab. It also reports on the application of the Broadband Outdoor Radiometer Calibration (BORCAL) method to reference cells. As a group, the factory calibrations for the crystalline silicon cells showed a small positive bias of +0.7%, with a tight distribution around this value. This is a very good result for the manufacturers. The two World Photovoltaic Scale cells constituted outliers, with a mean deviation of -0.8%, which is still within the reported uncertainty ranges. The data collected during the BORCAL procedures show a strong and systematic fluctuation of the apparent responsivity of the reference cells over the course of a day, which was anticipated. Consequently, the standard BORCAL R@45 calibration factors underestimated the reference cell responsivities by 1%-4%. An alternate responsivity indicator was developed, R'@30, incorporating adjustments for temperature response, spectral response, and directional response. This reduced the mean deviation to -0.4% compared to the Cell Lab. Further, two independent BORCAL sessions were able to produce the same values to within ±0.2%. This could be a viable approach for verifying reference cells using existing BORCAL infrastructure, for example, to evaluate their long-term stability. The daily profiles of responsivity obtained during the BORCAL sessions (provided in Appendix A) show clear patterns of differences between models, which provides evidence that their directional responses differ from each other; thus, two sensors that are calibrated correctly might give the same reading at 30° zenith angle but systematically different readings at 60°. The impact of these inconsistencies in field operation will be quantified using the long-term observations started in 2020.

14 SOLAR ENERGY↗
DOE OSTI2021

Performance of the ATLAS RPC detector and Level-1 muon barrel trigger at √(s)=13 TeV

The ATLAS experiment at the Large Hadron Collider (LHC) employs a trigger system consisting of a first-level hardware trigger (L1) and a software-based high-level trigger. The L1 muon trigger system selects muon candidates, assigns them to the correct LHC bunch crossing and classifies them into one of six transverse-momentum threshold classes. The L1 muon trigger system uses resistive-plate chambers (RPCs) to generate the muon-induced trigger signals in the central (barrel) region of the ATLAS detector. The ATLAS RPCs are arranged in six concentric layers and operate in a toroidal magnetic field with a bending power of 1.5 to 5.5 Tm. The RPC detector consists of about 3700 gas volumes with a total surface area of more than 4000 m2. This paper reports on the performance of the RPC detector and L1 muon barrel trigger using 60.8 fb-1 of proton-proton collision data recorded by the ATLAS experiment in 2018 at a centre-of-mass energy of 13 TeV. Detector and trigger performance are studied using Z boson decays into a muon pair. Measurements of the RPC detector response, efficiency, and time resolution are reported. Measurements of the L1 muon barrel trigger efficiencies and rates are presented, along with measurements of the properties of the selected sample of muon candidates. Measurements of the RPC currents, counting rates and mean avalanche charge are performed using zero-bias collisions. Finally, RPC detector response and efficiency are studied at different high voltage and front-end discriminator threshold settings in order to extrapolate detector response to the higher luminosity expected for the High Luminosity LHC.

47 OTHER INSTRUMENTATION↗
DOE OSTI2023

Interpreting Sunyaev–Zel’dovich observations with MillenniumTNG: mass and environment scaling relations

ABSTRACT Sunyaev–Zel’dovich (SZ) measurements can dramatically improve our understanding of the intergalactic medium and the role of feedback processes in galaxy formation, allowing us to calibrate important astrophysical systematics in cosmological constraints from weak lensing galaxy clustering surveys. However, the signal is only measured in a two-dimensional projection, and its correct interpretation relies on understanding the connection between observable quantities and the underlying intrinsic properties of the gas, in addition to the relation between the gas and the underlying matter distribution. One way to address these challenges is through the use of hydrodynamical simulations such as the high-resolution, large-volume MillenniumTNG suite. We find that measurements of the optical depth, τ, and the Compton-y parameter, Y, receive large line-of-sight contributions that can be removed effectively by applying a compensated aperture photometry filter. In contrast with other τ probes (e.g. X-rays and fast radio bursts), the kinematic SZ-inferred τ receives most of its signal from a confined cylindrical region around the halo due to the velocity decorrelation along the line of sight. Additionally, we perform fits to the Y–M and τ–M scaling relations and report best-fitting parameters adopting the smoothly broken power law formalism. We note that subgrid physics modelling can broaden the error bar on these by 30 per cent for intermediate-mass haloes (${\sim }10^{13} \, {\rm M}_{\odot }$). The scatter of the scaling relations can be captured by an intrinsic dependence on concentration and an extrinsic dependence on tidal shear. Finally, we comment on the effect of using galaxies rather than haloes in observations, which can bias the inferred profiles by ∼20 per cent for L* galaxies.

Hadzhiyska, Boryana (ORCID:0000000223123121)↗
DOE OSTI2024

Facets control charge separation during photoelectrochemical water oxidation with strontium titanate (SrTiO 3 ) single crystals

The photocatalytic overall water splitting reaction provides a pathway to hydrogen fuel from sunlight. Photocatalysts must achieve the reaction without the application of an external bias, which requires an effective charge separation mechanism. Photolabeling studies and electrostatic simulations for the well-known CoOOH/Al:SrTiO 3 /Rh/Cr 2 O 3 photocatalyst suggest that charge separation is driven by work function differences at the (100) and (110) facets of SrTiO 3 , which are electron and hole selective, respectively. Here we use hydrogen annealed SrTiO 3–x single crystals to obtain the first quantitative assessment of the charge separation ability of the (100), or (110), or (111) facets during oxygen evolution. Under UV illumination (60 mW cm –2 ), the crystals exhibit variable water oxidation photocurrents (0.34, 0.82, 1.36 mA cm –2 at 1.23 V versus RHE) and photovoltage values of 1.40, 1.52 and 1.52 V for (100), (110), and (111) SrTiO 3–x , respectively. A surface photovoltage increase in that same order (0.31 V < 0.57 V < 0.67 V) is confirmed independently with vibrating Kelvin probe surface photovoltage spectroscopy (VKP-SPV) under 375 nm (1.91 mW cm –2 ) illumination. Mott Schottky measurements in aqueous K 3/4 [Fe(CN) 6 ] reveal facet-dependent flatband positions of –0.58, –0.71, and –0.74 V RHE for the (100), (110), and (111) crystals respectively. This confirms that the photoelectrochemical water oxidation performance of SrTiO 3–x crystals is controlled by the work function of each facet, which determines the electron transfer barrier height of the respective solid–liquid junctions. After correcting for differences in electron donor concentrations, barriers are found to increase in the order (100) < (111) < (110) and differ by as much as 0.16 eV, similar to an earlier prediction. Altogehter, these results explain the charge separation mechanism in SrTiO 3 photocatalysts and highlight the need for faceted semiconductor crystals as light absorbers in particle-based photocatalysts.

08 HYDROGEN↗
DOE OSTI2020

A Causal Inference Model Based on Random Forests to Identify the Effect of Soil Moisture on Precipitation

Soil moisture influences precipitation mainly through its impact on land–atmosphere interactions. Understanding and correctly modeling soil moisture–precipitation (SM–P) coupling is crucial for improving weather forecasting and subseasonal to seasonal climate predictions, especially when predicting the persistence and magnitude of drought. However, the sign and spatial structure of SM–P feedback are still being debated in the climate research community, mainly due to the difficulty in establishing causal relationships and the high degree of nonlinearity in land–atmosphere processes. To this end, we developed a causal inference model based on the Granger causality analysis and a nonlinear machine learning model. This model includes three steps: nonlinear anomaly decomposition, nonlinear Granger causality analysis, and evaluation of the quality of SM–P feedback, which eliminates the nonlinear response of interannual and seasonal variability and the memory effects of climatic factors and isolates the causal relationship of local SM–P feedback. We applied this model by using National Climate Assessment–Land Data Assimilation System (NCA-LDAS) datasets over the United States. Here, the results highlight the importance of nonlinear atmosphere responses in land–atmosphere interactions. In addition, the strong feedback over the southwestern United States and the Great Plains both highlight the impacts of topographic factors rather than only the sensitivity of evapotranspiration to soil moisture. Furthermore, the SM–P index defined by our framework is used to benchmark Earth system models (ESMs), which provides a new metric for efficiently identifying potential model biases in modeling local land–atmosphere interactions and may help the development of ESMs in improving simulations of water cycle variability and extremes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Offshore wind energy forecasting sensitivity to sea surface temperature input in the Mid-Atlantic

Abstract. As offshore wind farm development expands, accurate wind resource forecasting over the ocean is needed. One important yet relatively unexplored aspect of offshore wind resource assessment is the role of sea surface temperature (SST). Models are generally forced with reanalysis data sets, which employ daily SST products. Compared with observations, significant variations in SSTs that occur on finer timescales are often not captured. Consequently, shorter-lived events such as sea breezes and low-level jets (among others), which are influenced by SSTs, may not be correctly represented in model results. The use of hourly SST products may improve the forecasting of these events. In this study, we examine the sensitivity of model output from the Weather Research and Forecasting model (WRF) 4.2.1 to different SST products. We first evaluate three different data sets: the Multiscale Ultrahigh Resolution (MUR25) SST analysis, a daily, 0.25∘ × 0.25∘ resolution product; the Operational Sea Surface Temperature and Ice Analysis (OSTIA), a daily, 0.054∘ × 0.054∘ resolution product; and SSTs from the Geostationary Operational Environmental Satellite 16 (GOES-16), an hourly, 0.02∘ × 0.02∘ resolution product. GOES-16 is not processed at the same level as OSTIA and MUR25; therefore, the product requires gap-filling using an interpolation method to create a complete map with no missing data points. OSTIA and GOES-16 SSTs validate markedly better against buoy observations than MUR25, so these two products are selected for use with model simulations, while MUR25 is at this point removed from consideration. We run the model for June and July of 2020 and find that for this time period, in the Mid-Atlantic, although OSTIA SSTs overall validate better against in situ observations taken via a buoy array in the area, the two products result in comparable hub-height (140 m) wind characterization performance on monthly timescales. Additionally, during hours-long flagged events (< 30 h each) that show statistically significant wind speed deviations between the two simulations, both simulations once again demonstrate similar validation performance (differences in bias, earth mover's distance, correlation, and root mean square error on the order of 10−1 or less), with GOES-16 winds validating nominally better than OSTIA winds. With a more refined GOES-16 product, which has been not only gap-filled but also assimilated with in situ SST measurements in the region, it is likely that hub-height winds characterized by GOES-16-informed simulations would definitively validate better than those informed by OSTIA SSTs.

17 WIND ENERGY↗
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↗