Search NASA⌕ Search

SEARCH · Search NASA

Results for “estimation and forecasting”

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

Hourly PM 2.5 Estimates across California from 2018 to 2023

This study presents a new data set of hourly PM 2.5 concentrations across California from 2018 to 2023 at a three-kilometer resolution. This data set was developed by assimilating observations from PurpleAir and the U.S. EPA Air Quality System monitors into wildfire smoke forecasts from the High-Resolution Rapid Refresh Smoke (HRRR-Smoke) model using the Gridpoint Statistical Interpolation (GSI) three-dimensional variational data assimilation framework. Archived forecasts of modeled wildfire smoke PM 2.5 from HRRR-Smoke create the background field for assimilation, which is then corrected using surface observations of total PM 2.5 . The resulting reanalysis from GSI provides an estimate of total PM 2.5 that minimizes error from both the observational and the model data. Validation results indicate strong performance, with monthly R 2 values ranging from 0.73 to 0.91 across the six-year data set, comparable to other PM 2.5 data sets. Case studies are presented for three major fire events, the 2018 Camp Fire, 2019 Kincade Fire, and 2020 Lightning Complex Fires to demonstrate the data set’s fidelity in resolving plume dynamics and local exposure patterns. Root-mean-squared error averaged over each month scales with average PM 2.5 concentrations, resulting in a low error under typical conditions but higher absolute errors during extreme smoke events. This is the first long-term, hourly PM 2.5 data set of its kind for California and enables the generation of subdaily exposure metrics, such as peak hourly concentrations, exceedance durations, and time-of-day exposure peaks. The novelty and strong validation of this data set make it a compelling resource for future studies on the impact and significance of subdaily PM 2.5 exposure.

PM2.5↗
DOE OSTI2024

Evaluating the Accuracy of Machine Learning Forecasts

To improve the accuracy of forecasting in machine learning, we must investigate multiple machine learning models and see how accurately they can predict values after training. We used seven machine learning models to try and get more accurate predictions. The models that were used were ARIMA, SES, MLP, CART, LightGBM, and XGBoost. We used a processed dataset from a Terminal at LAX that had the number of people traveling through terminal X every hour in March from 2015-2019. We trained our models with the dates March 6 - March 19 to predict the value for March 20th and the hours 6:00 am to 6:00 pm since those are the most popular traveling hours. By using the different models, we had varying results of accuracy when estimating the amount of people traveling through terminal X on March 20th. We know that machine learning models are helpful for forecasting and by seeing how accurately these models can predict, we can see how forecasting can be helpful for other issues. Using these methods, airports can use forecasting to predict the amount of people coming in and out and can use these predictions to prepare their resource management, operational efficiency, and overall passenger experience.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

A modeling examination of cloud seeding conditions under the warmer climate in Utah, USA

Assessment of how the wintertime cloud seeding conditions over Utah's mountains may change provides essential information for the state to work on increasing its snowpack in the warmer climate. Here in this study, data from the National Center for Atmospheric Research's high-resolution, convection-permitting Weather Research and Forecasting model are used to evaluate suitable conditions for wintertime cloud seeding over Utah's mountains, and to estimate future changes to these conditions under a pseudo-global warming scenario. While previous studies have examined historical cloud seeding conditions, this is the first study that explores the impact of future climate change on these conditions. Based on two atmospheric variables commonly used in cloud-seeding operations, i.e. average temperature from the surface to 1 km above ground level (-20 °C≤T≤-6 °C) and vertically integrated supercooled liquid water (SLW > 0.01 mm), Utah's mountains are suitable for seeding more than 20% of the time during winter. A greater rate of suitability exists in the northern and Uinta Mountains. In the warmer climate under a high-emissions scenario, these seeding conditions may wane as the percentage of precipitating clouds suitable for seeding, about 60% under current conditions, would decline to about 40% across the state with significance. This projected decrease is due to rising temperatures and a decreased frequency of precipitation events. These findings imply that climate warming will narrow the window of opportunity for winter cloud seeding operation in Utah.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Boundary layer moisture variability at the Atmospheric Radiation Measurement (ARM) Eastern North Atlantic observatory during marine conditions

Abstract. Boundary layer moisture variability at the Eastern North Atlantic (ENA) site during marine conditions is examined at monthly and daily timescales using 5 years of ground-based observations and output from the European Center for Medium range Weather Forecast (ECMWF) reanalysis model. The annual cycle of the mixed-layer total water budgets is presented to estimate the relative contribution of large-scale advection, local moisture tendency, entrainment, and precipitation to balance the moistening due to surface latent heat flux on monthly timescales. When marine conditions prevail, advection of colder and dry air from the north acts as an important moisture sink (∼ 50 % of the overall budget) during fall and winter driving the seasonality of the budget. Entrainment and precipitation contribute to the drying of the boundary layer (∼ 25 % and ∼ 15 % respectively), and the local change in moisture contributes to a small residual part. On a daily temporal scale, moist and dry mesoscale columns of vapor (∼ 10 km) are analyzed during 10 selected days of precipitating stratocumulus clouds. Adjacent moist and dry columns present distinct mesoscale features that are strongly correlated with clouds and precipitation. Dry columns adjacent to moist columns have more frequent and stronger downdrafts immediately below the cloud base. Moist columns have more frequent updrafts, stronger cloud-top cooling, and higher liquid water path and precipitation compared to the dry columns. This study highlights the complex interaction between large-scale and local processes controlling the boundary layer moisture and the importance of spatial distribution of vapor to support convection and precipitation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Real-Time Fault Tracking and Ground Motion Prediction for Large Earthquakes With HR-GNSS and Deep Learning

Earthquake early warning (EEW) systems aim to forecast the shaking intensity rapidly after an earthquake occurs and send warnings to affected areas before the onset of strong shaking. The system relies on rapid and accurate estimation of earthquake source parameters. However, it is known that source estimation for large ruptures in real-time is challenging, and it often leads to magnitude underestimation. In a previous study, we showed that machine learning, HR-GNSS, and realistic rupture synthetics can be used to reliably predict earthquake magnitude. This model, called Machine-Learning Assessed Rapid Geodetic Earthquake model (M-LARGE), can rapidly forecast large earthquake magnitudes with an accuracy of 99%. Here, we expand M-LARGE to predict centroid location and fault size, enabling the construction of the fault rupture extent for forecasting shaking intensity using existing ground motion models. We test our model in the Chilean Subduction Zone with thousands of simulated and five real large earthquakes. The result achieves an average warning time of 40.5 s for shaking intensity MMI4+, surpassing the 34 s obtained by a similar GNSS EEW model. Our approach addresses a critical gap in existing EEW systems for large earthquakes by demonstrating real-time fault tracking feasibility without saturation issues. This capability leads to timely and accurate ground motion forecasts and can support other methods, enhancing the overall effectiveness of EEW systems. Additionally, the ability to predict source parameters for real Chilean earthquakes implies that synthetic data, governed by our understanding of earthquake scaling, is consistent with the actual rupture processes.

58 GEOSCIENCES↗
DOE OSTI2021

Freight Analysis Framework Version 5 (FAF5) Base Year 2017 Data Development Technical Report

The Freight Analysis Framework (FAF) integrates data from a variety of sources to create a comprehensive national picture of freight movements among states and major metropolitan areas by all modes of transportation. The latest of this data series is FAF5, which is the fifth generation FAF and is benchmarked on Commodity Flow Survey (CFS) 2017. Except for FAF1 that provided estimates for truck, rail, and water tonnage for calendar year 1998, later generations of FAF (FAF2 through FAF5) were built based on their benchmark year CFS data, for 2002, 2007, 2012, and 2017 respectively. The FAF is produced under a partnership between Bureau of Transportation Statistics (BTS) and Federal Highway Administration (FHWA). As a major data product of the FAF program, the FAF regional database provides a national picture of freight flows to, from, and within the United States (among regions and states), by commodity and mode for the base year, as well as for forecasts up to 30 years into the future in a 5-year interval. Additional FAF data products also include FAF network flows database, where truck movements are routed onto the national highway network, estimates of annual projections, and synchronized historical data series. This report is a technical document prepared to describe the data sources and methodologies applied in the process of building the FAF5 base-year 2017 regional database, released as FAF5.0 in February 2021. This report offers a description of the diverse data sources and modeling methods used in constructing the base year FAF5 regional database. The FAF5 base-year database is used as the base for development of forecasts and for assignment of truck flows on highway network. Similarly, the FAF5 base-year database will be used as the base to generate FAF5 annual estimates. In addition to this report, users are encouraged to refer to the FAF5 User’s Guide, which provides basic information of the data, including definitions of the data attributes, information on how to access the data and tool, as well as detailed data dictionary and code tables.

42 ENGINEERING↗
DOE OSTI2024

Deep-learning-derived planetary boundary layer height from conventional meteorological measurements

Abstract. The planetary boundary layer (PBL) height (PBLH) is an important parameter for various meteorological and climate studies. This study presents a multi-structure deep neural network (DNN) model, which can estimate PBLH by integrating the morning temperature profiles and surface meteorological observations. The DNN model is developed by leveraging a rich dataset of PBLH derived from long-standing radiosonde records augmented with high-resolution micro-pulse lidar and Doppler lidar observations. We access the performance of the DNN with an ensemble of 10 members, each featuring distinct hidden-layer structures, which collectively yield a robust 27-year PBLH dataset over the southern Great Plains from 1994 to 2020. The influence of various meteorological factors on PBLH is rigorously analyzed through the importance test. Moreover, the DNN model's accuracy is evaluated against radiosonde observations and juxtaposed with conventional remote sensing methodologies, including Doppler lidar, ceilometer, Raman lidar, and micro-pulse lidar. The DNN model exhibits reliable performance across diverse conditions and demonstrates lower biases relative to remote sensing methods. In addition, the DNN model, originally trained over a plain region, demonstrates remarkable adaptability when applied to the heterogeneous terrains and climates encountered during the GoAmazon (Green Ocean Amazon; tropical rainforest) and CACTI (Cloud, Aerosol, and Complex Terrain Interactions; middle-latitude mountain) campaigns. These findings demonstrate the effectiveness of deep learning models in estimating PBLH, enhancing our understanding of boundary layer processes with implications for improving the representation of PBL in weather forecasting and climate modeling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

An Adaptive-Importance-Sampling-Enhanced Bayesian Approach for Topology Estimation in an Unbalanced Power Distribution System

The reliable operation of a power distribution system relies on a good prior knowledge of its topology and its system state. Although crucial, due to the lack of direct monitoring devices on the switch statuses, the topology information is often unavailable or outdated for the distribution system operators for real-time applications. Apart from the limited observability of the power distribution system, other challenges are the nonlinearity of the model, the complicated, unbalanced structure of the distribution system, and the scale of the system. To overcome the above challenges, we, in this paper, propose a Bayesian-inference framework that allows us to simultaneously estimate the topology and the state of a three-phase, unbalanced power distribution system. Specifically, by using the very limited number of measurements available that are associated with the forecast load data, we efficiently recover the full Bayesian posterior distributions of the system topology under both normal and outage operation conditions. This is performed through an adaptive importance sampling procedure that greatly alleviates the computational burden of the traditional Monte-Carlo (MC)-sampling-based approach while maintaining a good estimation accuracy. The simulations conducted on the IEEE 123-bus test system and an unbalanced 1282-bus system reveal the excellent performances of the proposed method.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Simulations of axion minihalos

The axion, motivated as a solution to the strong CP problem, is also a viable dark matter candidate. The axion field takes random values in causally disconnected regions if the symmetry breaking that establishes the particle occurs after inflation, leading to white-noise density fluctuations at low wave numbers and forming dense minihalos with sub-planetary masses subsequently. There have been two recent proposals that appear capable of testing this scenario, namely using pulsar timing arrays and studying cosmological microlensing caustics. Motivated by these proposals, we use N-body simulations to study the formation of substructures from white-noise density fluctuations. The density profiles of our relaxed axion minihalos can be described by the Navarro–Frenk–White profile, and the minihalos’ concentration number agrees well with a simple, physically-motivated model. We develop a semi-analytic formula to fit the mass function from our simulation, which agrees broadly at different redshifts and only differs at factor of two level from classic halo mass functions. This analytic mass function allows us to consider uncertainties in the post-inflation axion scenario, as well as extrapolate our high-redshift simulations results to the present. Our work estimates the present-day abundance of axion substructures, as is necessary for predicting their effect on cosmological microlensing caustics and pulsar timing. Our calculations suggest that if pulsar timing and microlensing probes can reach recent sensitivity forecasts, they may be sensitive to the post-inflation axion dark matter scenario, even when accounting for uncertainties pertaining to axion strings. For pulsar timing, the most significant caveat is whether axion minihalos are disrupted by stars, which our estimates show is mildly important at the most relevant masses. In conclusion, as our gravitational simulations are scale invariant, the results can be extended to models where the dark matter is comprised of other axion-like particles and even clusters of primordial black holes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Lake-Effect Snowstorm Events and Associated Snowfall Totals Integrated from NOAA Storm Reports, ERA5, and HRRR for the Laurentian Great Lakes (1997–2024)

Lake-effect snowstorms are localized, impactful winter weather phenomena that can generate substantial snowfall totals and pose significant challenges for forecasting, transportation, and regional infrastructure. To support the analysis and modeling of these events, this dataset compiles observational reports of lake-effect snowstorms alongside corresponding snowfall estimates derived from gridded atmospheric datasets. The observational component of the data originates from the National Weather Service (NWS) winter storm report, subset to lake-effect snow event type, covering 1997–2024. For each lake-effect snow event, this data provides the impacted county, event start and end datetimes at an hourly resolution, as well as relevant storm narratives. The complementary reanalysis-derived data is sourced from European Centre for Medium-Range Weather Forecasts (ECMWF) Reanalysis 5 (ERA5) and High-Resolution Rapid Refresh (HRRR) gridded data. For both gridded datasets, the maximum total snowfall (in units mm) was extracted, constrained by the county and datetimes specified by the observational report. ERA5 data covers the entire observational period (1997–2024), whereas HRRR data is only available from November 2016 – December 2024. Three CSV files are provided here: (1) the observational lake-effect snow event report, (2) ERA5 maximum snowfall detections for each event, and (3) HRRR maximum snowfall detections for each event. Relevant data from the observational files, such as impacted state and county, event datetimes, and event IDs, were included for convenience. Users can inspect and visualize the data using tools such as Microsoft Excel and Python pandas/matplotlib packages. This dataset may support a variety of applications, including climatological analyses of lake-effect snowfall, evaluation of snowfall representation in atmospheric datasets and numerical weather prediction models, and the development of machine learning approaches for detecting or predicting lake-effect snowfall events.

EARTH SCIENCE > ATMOSPHERE > PRECIPITATION > SOLID↗
DOE OSTI2025

Preserving nonlinear constraints in variational flow filtering data assimilation

Data assimilation aims to estimate the states of a dynamical system by optimally combining sparse and noisy observations of the physical system with uncertain forecasts produced by a computational model. The states of many dynamical systems of interest obey nonlinear physical constraints, and the corresponding dynamics is confined to a certain sub-manifold of the state space. Standard data assimilation techniques applied to such systems yield posterior states lying outside the manifold, violating the physical constraints. This work focuses on particle flow filters which use stochastic differential equations to evolve state samples from a prior distribution to samples from an observation-informed posterior distribution. The variational Fokker-Planck (VFP)—a generic particle flow filtering framework—is extended to incorporate non-linear, equality state constraints in the analysis. To this end, two algorithmic approaches that modify the VFP stochastic differential equation are discussed: (i) VFPSTAB, to inexactly preserve constraints with the addition of a stabilizing drift term, and (ii) VFPDAE, to exactly preserve constraints by treating the VFP dynamics as a stochastic differential-algebraic equation (SDAE). Additionally, an implicit-explicit time integrator is developed to evolve the VFPDAE dynamics. The strength of the proposed approach for constraint preservation in data assimilation is demonstrated on three test problems: the double pendulum, Korteweg-de-Vries, and the incompressible Navier-Stokes equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Joint impedance and facies inversion of time-lapse seismic data for improving monitoring of CO 2 incidentally stored from CO 2 EOR

Time-lapse seismic monitoring is an effective and proven technology for mapping the distribution of CO 2 in a subsurface reservoir. When injected CO 2 displaces other reservoir fluids, porous-medium properties are changed and thus the seismic impedance changes, causing time-lapse seismic amplitude differences in the injection zones. The analysis and interpretation of images created from these amplitude differences can provide information about reservoir architecture and the CO 2 migration within the reservoir. Incorporating seismic inversion and rock physics into the interpretation of time-lapse seismic data can considerably improve the modeling and monitoring to detect and assess the location of CO 2 over time. The joint inversion method presented in this paper has an integral representation of the geology in the inversion algorithm using elastic facies, which provides information about the spatial distribution of the geologic heterogeneities controlling the movement of fluids in the reservoir. The method was successfully applied to time-lapse seismic data from a mature oil field undergoing CO 2 enhanced oil recovery. The estimated seismic acoustic impedances and facies reflect the characteristics of individual geologic facies and fluid conditions of the reservoir subject to CO 2 injection. Finally, the probabilities estimated by the joint impedance and facies inversion for the reservoir's litho-fluid facies can be used for forecasting CO 2 saturation and pressure changes within the target reservoir.

4D seismic↗
DOE OSTI2021

Cherenkov Telescope Array sensitivity to the putative millisecond pulsar population responsible for the Galactic Centre excess

ABSTRACT The leading explanation of the Fermi Galactic Centre γ-ray excess is the extended emission from an unresolved population of millisecond pulsars (MSPs) in the Galactic bulge. Such a population would, along with the prompt γ-rays, also inject large quantities of electrons/positrons (e±) into the interstellar medium. These e± could potentially inverse-Compton (IC) scatter ambient photons into γ-rays that fall within the sensitivity range of the upcoming Cherenkov Telescope Array (CTA). In this article, we examine the detection potential of CTA to this signature by making a realistic estimation of the systematic uncertainties on the Galactic diffuse emission model at TeV-scale γ-ray energies. We forecast that, in the event that e± injection spectra are harder than E−2, CTA has the potential to robustly discover the IC signature of a putative Galactic bulge MSP population sufficient to explain the Galactic Centre excess for e± injection efficiencies in the range of ≈2.9–74.1 per cent, or higher, depending on the level of mismodelling of the Galactic diffuse emission components. On the other hand, for spectra softer than E−2.5, a reliable CTA detection would require an unphysically large e± injection efficiency of ${\gtrsim} 158{{\ \rm per\ cent}}$. However, even this pessimistic conclusion may be avoided in the plausible event that MSP observational and/or modelling uncertainties can be reduced. We further find that, in the event that an IC signal were detected, CTA can successfully discriminate between an MSP and a dark matter origin for the radiating e±.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Prospects for detecting heavy WIMP dark matter with the Cherenkov Telescope Array: The Wino and Higgsino

TeV-scale particles that couple to the standard model through the weak force represent a compelling class of dark matter candidates. The search for such weakly interacting massive particles has already spanned multiple decades, and whilst it has yet to provide any definitive evidence for their existence, viable parameter space remains. In this paper, we show that the upcoming Cherenkov Telescope Array (CTA) has significant sensitivity to uncharted parameter space at the TeV mass scale. To do so, we focus on two prototypical dark matter candidates, the Wino and Higgsino. Additionally, sensitivity forecasts for both models are performed including the irreducible background from misidentified cosmic rays, as well as a range of estimates for the Galactic emissions at TeV energies. For each candidate, we find substantial expected improvements over existing bounds from current imaging atmospheric Cherenkov telescopes. In detail, for the Wino we find a sensitivity improvement of roughly an order of magnitude in <σ v>, whereas for the Higgsino we demonstrate that CTA has the potential to become the first experiment that has sensitivity to the thermal candidate. Taken together, these enhanced sensitivities demonstrate the discovery potential for dark matter at CTA in the 1–100 TeV mass range.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Predictability and empirical dynamics of fisheries time series in the North Pacific

Previous studies have documented a strong relationship between marine ecosystems and large-scale modes of sea surface height (SSH) and sea surface temperature (SST) variability in the North Pacific such as the Pacific Decadal Oscillation and the North Pacific Gyre Oscillation. In the central and western North Pacific along the Kuroshio-Oyashio Extension (KOE), the expression of these modes in SSH and SST is linked to the propagation of long oceanic Rossby waves, which extend the predictability of the climate system to ~3 years. Using a multivariate physical-biological linear inverse model (LIM) we explore the extent to which this physical predictability leads to multi-year prediction of dominant fishery indicators inferred from three datasets (i.e., estimated biomasses, landings, and catches). We find that despite the strong autocorrelation in the fish indicators, the LIM adds dynamical forecast skill beyond persistence up to 5-6 years. By performing a sensitivity analysis of the LIM forecast model, we find that two main factors are essential for extending the dynamical predictability of the fishery indicators beyond persistence. The first is the interaction of the fishery indicators with the SST/SSH of the North and tropical Pacific. The second is the empirical relationship among the fisheries time series. This latter component reflects stock-stock interactions as well as common technological and human socioeconomic factors that may influence multiple fisheries and are captured in the training of the LIM. These results suggest that empirical dynamical models and machine learning algorithms, such as the LIM, provide an alternative and promising approach for forecasting key ecological indicators beyond the skill of persistence.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Deep-Learning-derived Boundary Layer Height from Meteorological Data over the SGP, GOAMAZON, CACTI

The planetary boundary-layer (PBL) height (PBLH) is an important parameter for various meteorological and climate studies. This study presents a multi-structure deep neural network (DNN) model, designed to estimate PBLH by integrating morning temperature profiles with surface meteorological observations. The DNN model is developed by leveraging a rich data set of PBLH derived from long-standing radiosonde records and augmented with high-resolution micropulse lidar and Doppler lidar observations. We access the performance of the DNN with an ensemble of 10 members, each featuring distinct hidden layer structures, which collectively yield a robust 27-year PBLH data set over the Southern Great Plains from 1994 to 2020. The influence of various meteorological factors on PBLH is rigorously analyzed through the importance test. Moreover, the DNN model's accuracy is evaluated against radiosonde observations and juxtaposed with conventional remote-sensing methodologies, including Doppler lidar, ceilometer, Raman lidar, and micropulse lidar. The DNN model exhibits reliable performance across diverse conditions and demonstrates lower biases relative to remote-sensing methods. In addition, the DNN model, originally trained over a plain region, demonstrates remarkable adaptability when applied to the heterogeneous terrains and climates encountered during the GoAmazon (tropical rainforest) and CACTI (middle-latitude mountain) campaigns. These findings demonstrate the effectiveness of deep learning models in estimating PBLH, enhancing our understanding of boundary-layer dynamics with implications for enhancing the representation of PBL in weather forecasting and climate modeling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Considerations for Long-Term Load Forecasting in Morocco

There are many factors that determine how demand for electricity may change over time. These factors include GDP, population size, and technology diffusion and adoption. We employ a simple extrapolation of current trends in Morocco GDP to estimate how the peak demand and annual consumption may change through the year 2030. We discuss the many factors that this approach does not take into account (such as adoption of air conditioning, electric vehicles, and distributed generation).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Primordial non-Gaussianity with angular correlation function: integral constraint and validation for DES

ABSTRACT Local primordial non-Gaussianity (PNG) is a promising observable of the underlying physics of inflation, characterized by $f_{\rm NL}^{\rm loc}$. We present the methodology to measure $f_{\rm NL}^{\rm loc}$ from the Dark Energy Survey (DES) data using the two-point angular correlation function (ACF) with scale-dependent bias. One of the focuses of the work is the integral constraint (IC). This condition appears when estimating the mean number density of galaxies from the data and is key in obtaining unbiased $f_{\rm NL}^{\rm loc}$ constraints. The methods are analysed for two types of simulations: ∼246 goliat-png N-body small area simulations with fNL equal to −100 and 100, and 1952 Gaussian ice-cola mocks with fNL = 0 that follow the DES angular and redshift distribution. We use the ensemble of goliat-png mocks to show the importance of the IC when measuring PNG, where we recover the fiducial values of fNL within the 1σ when including the IC. In contrast, we found a bias of ΔfNL ∼ 100 when not including it. For a DES-like scenario, we forecast a bias of ΔfNL ∼ 23, equivalent to 1.8σ, when not using the IC for a fiducial value of fNL = 100. We use the ice-cola mocks to validate our analysis in a realistic DES-like set-up finding it robust to different analysis choices: best-fitting estimator, the effect of IC, BAO damping, covariance, and scale choices. We forecast a measurement of fNL within σ(fNL) = 31 when using the DES-Y3 BAO sample, with the ACF in the 1 deg < θ < 20 deg range.

79 ASTRONOMY AND ASTROPHYSICS↗