Search NASA⌕ Search

SEARCH · Search NASA

Results for “modeling discrepancy”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
NASA NTRS1996

Responses of the Tropical Pacific to Wind Forcing as Observed by Spaceborne Sensors and Simulated by an Ocean General Circulation Model

In this study, satellite observations, in situ measurements, and model simulations are combined to assess the oceanic response to surface wind forcing in the equatorial Pacific. The surface wind fields derived from observations by the spaceborne special sensor microwave imager (SSM/I) and from the operational products of the European Centre for Medium-Range Weather Forecasts (ECMWF) are compared. When SSM/I winds are used to force a primitive-equation ocean general circulation model (OGCM), they produce 3 C more surface cooling than ECMWF winds for the eastern equatorial Pacific during the cool phase of an El Nino-Southern Oscillation event. The stronger cooling by SSM/I winds is in good agreement with measurements at the moored buoys and observations by the advanced very high resolution radiometer, indicating that SSM/I winds are superior to ECMWF winds in forcing the tropical ocean. In comparison with measurements from buoys, tide gauges, and the Geosat altimeter, the OGCM simulates the temporal variations of temperature, steric, and sea level changes with reasonable realism when forced with the satellite winds. There are discrepancies between model simulations and observations that are common to both wind forcing fields, one of which is the simulation of zonal currents; they could be attributed to model deficiencies. By examining model simulations under two winds, vertical heat advection and uplifting of the thermocline are found to be the dominant factors in the anomalous cooling of the ocean mixed layer.

Liu, W. Timothy↗
DOE OSTI2022

Peatland Loss in Southeast Asia Contributing to U.S. Biofuel’s Greenhouse Gas Emissions

Land use change (LUC) induced by biofuel production could lead to greenhouse gas (GHG) emissions, which potentially increase biofuel’s carbon intensity. Among the sources of LUC-related emissions for soy biodiesel, the contribution from peatland loss to agricultural plantations in Southeast Asia remains uncertain. Here, in this study, we analyzed LUC in Malaysia and Indonesia and modeled its impacts on the GHG emissions of soy biodiesel produced in the United States. It shows that oil palm plantations have more than doubled over 2001–2016 and the area of palm-on-peatlands (PoP) has expanded 3.7 times. Over new palm plantations, the share of PoP is about 19% regardless of time and location and the emission factor (EF) for peatland-to-palm conversion is estimated to be 41.5 Mg CO 2 ha –1 yr –1 . With these updates on PoP and EF, the contribution of peatland loss (0.7–5.1 g CO 2 e MJ –1 ) to biodiesel emissions is only 40–65% of previous estimates, which reduces discrepancies among model simulations used by different agencies. Based on emerging evidence on LUC and related carbon changes, our analysis reexamines regional peatland loss and its impacts on LUC emissions modeling and provides new insights into the estimation of LUC impacts on biofuels’ carbon intensity.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2012

Longwave Band-by-band Cloud Radiative Effect and its Application in GCM Evaluation

The cloud radiative effect (CRE) of each longwave (LW) absorption band of a GCM fs radiation code is uniquely valuable for GCM evaluation because (1) comparing band-by-band CRE avoids the compensating biases in the broadband CRE comparison and (2) the fractional contribution of each band to the LW broadband CRE (f(sub CRE)) is sensitive to cloud top height but largely insensitive to cloud fraction, presenting thus a diagnostic metric to separate the two macroscopic properties of clouds. Recent studies led by the first author have established methods to derive such band ]by ]band quantities from collocated AIRS and CERES observations. We present here a study that compares the observed band-by-band CRE over the tropical oceans with those simulated by three different atmospheric GCMs (GFDL AM2, NASA GEOS-5, and CCCma CanAM4) forced by observed SST. The models agree with observation on the annual ]mean LW broadband CRE over the tropical oceans within +/-1W/sq m. However, the differences among these three GCMs in some bands can be as large as or even larger than +/-1W/sq m. Observed seasonal cycles of f(sub CRE) in major bands are shown to be consistent with the seasonal cycle of cloud top pressure for both the amplitude and the phase. However, while the three simulated seasonal cycles of f(sub CRE) agree with observations on the phase, the amplitudes are underestimated. Simulated interannual anomalies from GFDL AM2 and CCCma CanAM4 are in phase with observed anomalies. The spatial distribution of f(sub CRE) highlights the discrepancies between models and observation over the low-cloud regions and the compensating biases from different bands.

Huang, Xianglei↗
NASA NTRS2021

Mid-Pliocene Atlantic Meridional Overturning Circulation simulated in PlioMIP2

In the Pliocene Model Intercomparison Project Phase 2 (PlioMIP2), coupled climate models have been used to simulate an interglacial climate during the mid-Piacenzian warm period (mPWP; 3.264 to 3.025 Ma). Here, we compare the Atlantic Meridional Overturning Circulation (AMOC), poleward ocean heat transport and sea surface warming in the Atlantic simulated with these models. In PlioMIP2, all models simulate an intensified mid-Pliocene AMOC. However, there is no consistent response in the simulated Atlantic ocean heat transport nor in the depth of the Atlantic overturning cell. The models show a large spread in the simulated AMOC maximum, the Atlantic ocean heat transport and the surface warming in the North Atlantic. Although a few models simulate a surface warming of ~8–12°C in the North Atlantic, similar to the reconstruction from Pliocene Research, Interpretation and Synoptic Mapping (PRISM) version 4, most models appear to underestimate this warming. The large model spread and model–data discrepancies in the PlioMIP2 ensemble do not support the hypothesis that an intensification of the AMOC, together with an increase in northward ocean heat transport, is the dominant mechanism for the mid-Pliocene warm climate over the North Atlantic.

Pliocene Model Intercomparison Project Phase 2 (Pl↗
NASA NTRS2012

Carbon Consequences of Forest Disturbance and Recovery Across the Conterminous United States

Forests of North America are thought to constitute a significant long term sink for atmospheric carbon. The United States Forest Service Forest Inventory and Analysis (FIA) program has developed a large data base of stock changes derived from consecutive estimates of growing stock volume in the US. These data reveal a large and relatively stable increase in forest carbon stocks over the last two decades or more. The mechanisms underlying this national increase in forest stocks may include recovery of forests from past disturbances, net increases in forest area, and growth enhancement driven by climate or fertilization by CO2 and Nitrogen. Here we estimate the forest recovery component of the observed stock changes using FIA data on the age structure of US forests and carbon stocks as a function of age. The latter are used to parameterize forest disturbance and recovery processes in a carbon cycle model. We then apply resulting disturbance/recovery dynamics to landscapes and regions based on the forest age distributions. The analysis centers on 28 representative climate settings spread about forested regions of the conterminous US. We estimate carbon fluxes for each region and propagate uncertainties in calibration data through to the predicted fluxes. The largest recovery-driven carbon sinks are found in the South central, Pacific Northwest, and Pacific Southwest regions, with spatially averaged net ecosystem productivity (NEP) of about 100 g C / square m / a driven by forest age structure. Carbon sinks from recovery in the Northeast and Northern Lake States remain moderate to large owing to the legacy of historical clearing and relatively low modern disturbance rates from harvest and fire. At the continental scale, we find a conterminous U.S. forest NEP of only 0.16 Pg C/a from age structure in 2005, or only 0.047 Pg C/a of forest stock change after accounting for fire emissions and harvest transfers. Recent estimates of NEP derived from inventory stock change, harvest, and fire data show twice the NEP sink we derive from forest age distributions. We discuss possible reasons for the discrepancies including modeling errors and the possibility of climate and/or fertilization (CO2 or N) growth enhancements.

Williams, Christopher A.↗
DOE OSTI2023

Evaluation of liquid cloud albedo susceptibility in E3SM using coupled eastern North Atlantic surface and satellite retrievals

The impact of aerosol number concentration on cloud albedo is a persistent source of spread in global climate predictions due to multi-scale, interactive atmospheric processes that remain difficult to quantify. We use 5 years of geostationary satellite and surface retrievals at the US Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) eastern North Atlantic (ENA) site in the Azores to evaluate the representation of liquid cloud albedo susceptibility for overcast cloud scenes in the DOE Energy Exascale Earth System Model version 1 (E3SMv1) and provide possible reasons for model–observation discrepancies. The overall distribution of surface 0.2 % CCN concentration values is reasonably simulated, but simulated liquid water path (LWP) is lower than observed and layer mean droplet concentration (N d ) comparisons are highly variable depending on the N d retrieval technique. E3SMv1's cloud albedo is greater than observed for given LWP and N d values due to a lower cloud effective radius than observed. However, the simulated albedo response to N d is suppressed due to a correlation between the solar zenith angle (SZA) and N d created by the seasonal cycle that is not observed. Controlling for this effect by examining the cloud optical depth (COD) shows that E3SMv1's COD response to CCN concentration is greater than observed. For surface-based retrievals, this is only true after controlling for cloud adiabaticity because E3SMv1's adiabaticities are much lower than observed. Assuming a constant adiabaticity in surface retrievals as done in top-of-atmosphere (TOA) retrievals narrows the retrieved ln N d distribution, which increases the cloud albedo sensitivity to ln N d to match the TOA sensitivity. The greater sensitivity of COD to CCN is caused by a greater Twomey effect in which the sensitivity of N d to CCN is greater than observed for TOA-retrieved N d , and once model–observation cloud adiabaticity differences are removed, this is also true for surface-retrieved N d . The LWP response to N d in E3SMv1 is overall negative as observed. Despite reproducing the observed LWP–N d relationship, observed clouds become much more adiabatic as N d increases, while E3SMv1 clouds do not, associated with more heavily precipitating clouds that are partially but not completely caused by deeper clouds and weaker inversions in E3SMv1. These cloud property differences indicate that the negative LWP–N d relationship is likely not caused by the same mechanisms in E3SMv1 and observations. The negative simulated LWP response also fails to mute the excessively strong Twomey effect, highlighting potentially important confounding factor effects that likely render the LWP–N d relationship non-causal. N d retrieval scales and assumptions, particularly related to cloud adiabaticity, contribute to substantial spreads in the model–observation comparisons, though enough consistency exists to suggest that aerosol activation, drizzle, and entrainment processes are critical areas to focus E3SMv1 development for improving the fidelity of aerosol–cloud interactions in E3SM.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2006

Radiative Effect of Clouds on Tropospheric Chemistry in a Global Three-Dimensional Chemical Transport Model

Clouds exert an important influence on tropospheric photochemistry through modification of solar radiation that determines photolysis frequencies (J-values). We assess the radiative effect of clouds on photolysis frequencies and key oxidants in the troposphere with a global three-dimensional (3-D) chemical transport model (GEOS-CHEM) driven by assimilated meteorological observations from the Goddard Earth Observing System data assimilation system (GEOS DAS) at the NASA Global Modeling and Assimilation Office (GMAO). We focus on the year of 2001 with the GEOS-3 meteorological observations. Photolysis frequencies are calculated using the Fast-J radiative transfer algorithm. The GEOS-3 global cloud optical depth and cloud fraction are evaluated and generally consistent with the satellite retrieval products from the Moderate Resolution Imaging Spectroradiometer (MODIS) and the International Satellite Cloud Climatology Project (ISCCP). Results using the linear assumption, which assumes linear scaling of cloud optical depth with cloud fraction in a grid box, show global mean OH concentrations generally increase by less than 6% because of the radiative effect of clouds. The OH distribution shows much larger changes (with maximum decrease of approx.20% near the surface), reflecting the opposite effects of enhanced (weakened) photochemistry above (below) clouds. The global mean photolysis frequencies for J[O1D] and J[NO2] in the troposphere change by less than 5% because of clouds; global mean O3 concentrations in the troposphere increase by less than 5%. This study shows tropical upper tropospheric O3 to be less sensitive to the radiative effect of clouds than previously reported (approx.5% versus approx.20-30%). These results emphasize that the dominant effect of clouds is to influence the vertical redistribution of the intensity of photochemical activity while global average effects remain modest, again contrasting with previous studies. Differing vertical distributions of clouds may explain part, but not the majority, of these discrepancies between models. Using an approximate random overlap or a maximum-random overlap scheme to take account of the effect of cloud overlap in the vertical reduces the impact of clouds on photochemistry but does not significantly change our results with respect to the modest global average effect.

Liu, Hongyu↗
DOE OSTI2026

Scaling dynamics in low-salt-rejection reverse osmosis for high-salinity produced water desalination: Mechanistic modeling and membrane autopsy

Membrane scaling remains a critical barrier to the reliable operation of desalination systems, particularly for hypersaline produced water (PW) treatment. This study fills the knowledge gap of autopsy-based model validation for PW desalination by elucidating scaling mechanisms in a Low-Salt-Rejection Reverse Osmosis (LSRRO) system through the integration of pilot-scale experimentation and complementary modeling approaches. A semi-empirical modeling framework was developed and applied to a multistage pilot LSRRO system equipped with nanofiltration and RO membranes treating high-salinity PW from the Permian Basin. Water quality analysis showed that total dissolved solids decreased from ~130,000 mg/L to ~1900 mg/L in the permeate, then further reduced to ~300 mg/L by a second-pass RO. Two different thermodynamic modeling approaches were evaluated: the first extends the LSRRO framework by incorporating system complexity and scaling phenomena, whereas the second method explicitly captures concentration polarization in localized supersaturation. Both methods illustrate the tendency for carbonate and sulfate scaling throughout the stages. Membrane autopsies revealed a silica-dominated deposit matrix, localized CaSO 4 at Stage 2, and minor barite/celestite despite their prominence in model predictions. Quantum-chemical calculations indicated silica scaling can be rationalized by favorable adsorption of H 4 SiO 4 on Fe-oxide surfaces (ΔG ≈ −44 kJ/mol), providing a kinetic pathway for interfacial inorganic polymerization even when bulk equilibrium predictions are conservative. Overall, the thermodynamic scaling modeling and membrane autopsy revealed heterogeneous, localized deposits with limited impact on LSRRO performance, while quantum analysis rationalized the thermodynamically unfavorable precipitation formation under bulk equilibrium, reconciling model–autopsy discrepancies. These insights support targeted pretreatment and silica-specific antiscalants to extend membrane lifetime and increase recovery, providing a transferable framework for hypersaline water desalination systems. The combined experimental–computational approach provides new mechanistic insight into scaling in hypersaline membrane systems and establishes a transferable framework for predicting and mitigating scaling in next-generation desalination technologies.

Low-salt-rejection reverse osmosis↗
DOE OSTI2026

First study of the nuclear response to fast hadrons via angular correlations between pions and slow protons in electron-nucleus scattering

We report on the first measurement of angular correlations between high-energy pions and slow protons in electron-nucleus (eA) scattering, providing a new probe of how a nucleus responds to a fast-moving quark. The experiment employed the CLAS detector with a 5-GeV electron beam incident on deuterium, carbon, iron, and lead targets. For heavier nuclei, the pion-proton correlation function is more spread-out in azimuth than for lighter ones, and this effect is more pronounced in the πp channel than in earlier ππ studies. The proton-to-pion yield ratio likewise rises with nuclear mass, although the increase appears to saturate for the heaviest targets. These trends are qualitatively reproduced by state-of-the-art eA event generators, including BeAGLE, eHIJING, and GiBUU, indicating that current descriptions of target fragmentation rest on sound theoretical footing. At the same time, the precision of our data exposes model-dependent discrepancies, delineating a clear path for future improvements in the treatment of cold-nuclear matter effects in eA scattering.

Correlations↗
DOE OSTI2024

Synthesis and characterization of isotopically barcoded nickel, molybdenum, and tungsten taggants for intentional nuclear forensics

Intentional nuclear forensics is a concept wherein the deliberate addition of benign and persistent material signatures to nuclear material can be used to reduce the time between the discovery of material outside of regulatory control and determination of its original provenance. One concept within intentional nuclear forensics involves the use of perturbed stable isotopes to generate unique isotope ratio “barcodes” to encode information (e.g., production batch, location, etc.) and track material throughout the nuclear fuel cycle. Synthesis of taggant species of nickel (Ni), molybdenum (Mo), and tungsten (W) was undertaken via a double-spike mechanism, wherein two highly enriched isotopes of interest per elemental taggant were mixed to form an enriched “double-spike” which was subsequently isotopically diluted with bulk material having a natural isotopic composition. Two taggant species perturbing isotopic ratios, alpha (α) and beta (β), for each of Ni, Mo, and W were synthesized. Independent measurements of double spikes and alpha and beta taggant species agreed within uncertainty and are clearly resolvable from natural compositions. High-precision analyses were independently performed by MC-ICP-MS at two U.S. National Laboratories, with consensus values and uncertainties calculated for all samples. Observed isotopic perturbations in the final taggant species measured on the order of hundreds to thousands of permille (‰) with respect to natural for isotope ratios of interest (e.g., 60 Ni/ 58 Ni, 100 Mo/ 98 Mo, 186 W/ 183 W). Discrepancies between modeled and measured isotopic compositions were observed and are largely attributed to imprecise vendor assay values for starting materials. Using measured starting material compositions as inputs for the mixing model improved the level of agreement between predicted and measured α and β taggant isotope ratios. Overall, characterization of all taggant species demonstrates that this “barcode” concept could have viability for use in nuclear forensics. Finally, it is expected that for any two-isotope mixing array dozens of isotopic barcodes could be encoded into a material system and subsequently resolved utilizing modern mass spectrometric methods.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2022

Transient vibration and product formation of photoexcited CS 2 measured by time-resolved x-ray scattering

Here, we have observed details of the internal motion and dissociation channels in photoexcited carbon disulfide (CS 2 ) using time-resolved x-ray scattering (TRXS). Photoexcitation of gas-phase CS 2 with a 200 nm laser pulse launches oscillatory bending and stretching motion, leading to dissociation of atomic sulfur in under a picosecond. During the first 300 fs following excitation, we observe significant changes in the vibrational frequency as well as some dissociation of the C–S bond, leading to atomic sulfur in the both 1 D and 3 P states. Beyond 1400 fs, the dissociation is consistent with primarily 3 P atomic sulfur dissociation. This channel-resolved measurement of the dissociation time is based on our analysis of the time-windowed dissociation radial velocity distribution, which is measured using the temporal Fourier transform of the TRXS data aided by a Hough transform that extracts the slopes of linear features in an image. The relative strength of the two dissociation channels reflects both their branching ratio and differences in the spread of their dissociation times. Measuring the time-resolved dissociation radial velocity distribution aids the resolution of discrepancies between models for dissociation proposed by prior photoelectron spectroscopy work.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Cosmological implications of DESI DR2 BAO measurements in light of the latest ACT DR6 CMB data

We report cosmological results from the Dark Energy Spectroscopic Instrument (DESI) measurements of baryon acoustic oscillations (BAO) when combined with recent data from the Atacama Cosmology Telescope (ACT). By jointly analyzing ACT and Planck data and applying conservative cuts to overlapping multipole ranges, we assess how different 𝑃⁢𝑙⁢𝑎⁢𝑛⁢𝑐⁢𝑘 + ACT dataset combinations affect consistency with DESI. While ACT alone exhibits a tension with DESI exceeding 3⁢𝜎 within the Λ⁢ CDM model, this discrepancy is reduced when ACT is analyzed in combination with Planck. For our baseline DESI DR2 BAO + 𝑃⁢𝑙⁢𝑎⁢𝑛⁢𝑐⁢𝑘 P⁢R⁢4 + ACT likelihood combination, the preference for evolving dark energy over a cosmological constant is about 3⁢𝜎, increasing to over 4⁢𝜎 with the inclusion of type Ia supernova data. While the dark energy results remain quite consistent across various combinations of Planck and ACT likelihoods with those obtained by the DESI collaboration, the constraints on neutrino mass are more sensitive, ranging from ∑𝑚 𝜈 < 0.061 eV in our baseline analysis, to ∑𝑚 𝜈 < 0.077 eV (95% confidence level) in the CMB likelihood combination chosen by ACT when imposing the physical prior ∑𝑚 𝜈 > 0 eV.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Assessing Impacts of Waves on Hub-Height Winds off the U.S. West Coast Using Lidar Buoys and Coupled Modeling Approaches

Given the importance of offshore wind energy development to the U.S. clean energy targets, it is vital to be able to characterize the wind resource in that environment accurately. Toward that end, two Bureau of Ocean Energy Management buoys equipped with Doppler lidar are being maintained by Pacific Northwest National Laboratory on behalf of the Department of Energy and deployed to regions of potential offshore wind development. In addition to standard meteorological and oceanographic measurements, the buoys document the wind profile between about 40 m and 250 m above the sea surface through Doppler lidar retrievals. After a multiyear deployment of two buoys along the U.S. East Coast, the buoys were redeployed to the U.S. West coast from 2020 – 2022 to locations near the Humboldt and Morro Bay lease areas. The buoys provide nearly continuous, multiyear datasets that can be used to evaluate predictions of hub-height (~100 m) wind speed for standard atmospheric models in the region. In the absence of measurements at the study site, offshore wind developers rely on model-based data to assess site conditions. Potential sources of model error in this environment include under-resolution or misrepresentation of coastal topographically forced flows, marine boundary layer dynamics and the evolution of their associated cloud and turbulence fields, the role of upwelling and other currents on surface heat fluxes into the boundary layer, and the impact of wave fields on surface momentum fluxes and thus the wind speed profile. In particular, most predictive models of wind speed do not predict wave fields at all, relying on parameterizations to represent their effects. In thus study, we focus on evaluating the role of wind / wave interactions on modeled hub-height wind speed and error by using the Coupled Ocean–Atmosphere–Wave–Sediment–Transport Modeling System to capture two-way interactions between an atmospheric model (Weather Research and Forecasting (WRF)) and a wave model (WAVEWATCHIII (WW3)) and compare to both stand-alone WRF and one-way coupled WRF / WW3 configurations. Our approach is similar to that used in Gaudet et al. (2022) to evaluate wind / wave coupling over the U.S. East Coast, but applied to the very different environment of the U.S. West Coast. We show examples for two cases, a cold-season frontal case and a warm-season low-level jet case. We find that wind / wave coupling makes little impact on model error for these cases at the location of the lidar buoys, for which other misrepresentations of model physics seems to be responsible for model-observation discrepancies. However, domain-wide evaluations, which also make use of the National Buoy Data Center network, show that a two-way coupling approach is less prone to systematic errors in the hub-height wind field than the one-way coupled approach. WRF resolution of kilometer-scale or less is needed to properly capture the sharp wind speed gradients that can be found along the coastline, and WW3 simulations driven by the downscaled WRF produce better bulk and spectral wave fields when compared to observations. Implications of the results for wind resource characterization are then discussed.

17 WIND ENERGY↗
DOE OSTI2025

Bayesian Calibration of Nuclear Graphite Property Models Accounting for Model Inadequacy and Impacts on Component Performance

Nuclear-grade structural graphite is extensively utilized in the core designs of various advanced nuclear reactors. In the reactor environment, graphite is subjected to prolonged exposure to extreme conditions, including high temperatures, radiation, and potentially molten salt and oxygen. Such exposure can induce several degradation mechanisms in graphite, including nonuniform volumetric strains caused by irradiation and thermal expansion, leading to stresses that may compromise the performance of graphite components. Assessing component integrity requires accurate models of graphite's thermomechanical response. This report documents the Bayesian calibration of thermomechanical properties for nuclear-grade graphite and their application to graphite component modeling and simulation using the Grizzly code. As part of this work, uncertainty-quantified models were developed for the elastic modulus, coefficient of thermal expansion, irradiation-induced dimensional change, and irradiation-induced creep for graphite grades IG-110, NBG-18, NBG-17, PCEA, and 2114. Using a hierarchical Bayesian approach, multiple experimental data sources were combined to develop Gaussian process models for the properties. Using the Kennedy O'Hagan framework, the uncertainties due to inadequacies in the model and the inherent spread in the experimental data were quantified for three different models. These uncertainty-quantified models, with a model-form correction, were subsequently applied to a coupled-physics simulation of representative graphite components, revealing that the uncertainties have a large impact on the components' deformation.

36 - MATERIALS SCIENCE↗
DOE OSTI2026

Characterization of the Polarization Beam Response of SPT-3G Using Point Sources

Precise measurements of cosmic microwave background (CMB) polarization require rigorous control of instrumental systematics. For the South Pole Telescope's third-generation camera (SPT-3G), which observes in three bands centered near 95, 150, and 220 GHz, accurate beam characterization is critical for interpreting the polarized mm-wave sky. We present direct measurements of SPT-3G's polarized beam response from observations of 100 bright extragalactic point sources. Previous SPT-3G power spectrum analyses introduced a phenomenological parameter, $β_{\rm pol}$, to describe the polarization preserved in beam sidelobes, and found evidence for significant depolarization from the requirement of inter-frequency polarization power spectrum consistency. Our direct measurements yield $β_{\rm pol}=0.89\pm0.10$ at 95 GHz, $1.08\pm0.10$ at 150 GHz, and $0.90\pm0.22$ at 220 GHz, indicating minimal sidelobe depolarization. We validate these results with systematic tests of posterior sampling versus bootstrap resampling, real-space versus Fourier-space analysis, temperature-to-polarization leakage handling, covariance determination, and source selection. Compared to values inferred from previous cosmological analyses, our results differ by an effective $1.3σ$. This apparent discrepancy is model dependent, because the point source analysis derives much of its $β_{\rm pol}$ constraining power from higher multipoles than the power spectrum analysis. These measurements therefore admit three explanations for the frequency-dependent residuals observed in the power spectrum analysis: a statistical fluctuation, the need for more sophisticated polarized beam models, or systematics other than beam depolarization.

de Haan, T. [KEK, Tsukuba]↗
DOE OSTI2025

SpecDis: Value Added Distance Catalog for 4 Million Stars from DESI Year-1 Data

We present the SpecDis value-added stellar distance catalog accompanying DESI Data Release 1. SpecDis trains a feed-forward neural network (NN) with Gaia parallaxes and gets the distance estimates. To build up an unbiased training sample, we do not apply selections on parallax error or signal-to-noise (S/N) of the stellar spectra, and instead, we incorporate parallax error into the loss function. Moreover, we employ principal component analysis to reduce the noise and dimensionality of stellar spectra. Validated by independent external samples of member stars with precise distances from globular clusters, dwarf galaxies, stellar streams, combined with blue horizontal branch stars, we demonstrate that our distance measurements show no significant bias up to 100 kpc, and are much more precise than Gaia parallax beyond 7 kpc. The median distance uncertainties are 23%, 19%, 11%, and 7% for S/N < 20, 20 ≤ S/N < 60, 60 ≤ S/N < 100, and S/N ≥ 100. Selecting stars with ${\mathrm{log}}\,g\lt 3.8$ and distance uncertainties smaller than 25%, we have more than 74,000 giant candidates within 50 kpc of the Galactic center and 1500 candidates beyond this distance. Additionally, we develop a Gaussian mixture model to identify unresolvable equal-mass binaries by modeling the discrepancy between the NN-predicted and the geometric absolute magnitudes from Gaia parallaxes and identify 120,000 equal-mass binary candidates. Our final catalog provides distances and distance uncertainties for >4 million stars, offering a valuable resource for Galactic astronomy.

astronomy data analysis↗