Search NASA⌕ Search

SEARCH · Search NASA

Results for “modeling and filtering”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Earth-Centered, Earth-Fixed Inertial Navigation System & Error-State Kalman Filter Reference Manual

This is a self-contained reference document that derives the equations necessary to build a combined inertial navigation system and error-state Kalman filter. Coordinate transform, linear time invariant system, inertial sensing, and error-state Kalman filtering theory is built up from first principles. This theory is then leveraged to derive the system equations for two combined inertial navigation system and error-state Kalman filters: (1) a 15-state system modeling white-noise-integrating accelerometer and gyroscope biases, and (2) a 39-state system modeling static and first-order Gauss-Markov accelerometer and gyroscope biases, scale factor errors, and cross-axis sensitivity errors.

42 ENGINEERING↗
DOE OSTI2025

Studies of optical and chemical properties of aged and fresh biomass burning absorbing aerosols for climate models

Particulate emissions from burning fuels from sub-Saharan Africa is currently poorly characterized. With increasing population and industrialization, emissions are expected to grow. Optical methods for measuring light absorption in the ultraviolet (UV) and visible (vis) regions of the spectrum are crucial for understanding the optical properties of BB aerosols (Bond and Bergstrom, 2006; Lack and Cappa, 2010; Moschos et al., 2021). One such instrument is the dual-spot AE33 model aethalometer is a common filter-based instrument that can be used to measure aerosol absorption properties (Hansen et al., 1984). This instrument has many figures of merit; however, the method's limitations include potential biases from scattering within the filter material, multiple scattering effects of the particles, filter loading (i.e. particle shadowing), and the need for wavelength-specific correction factors (Cλ), which depend on accurate intercomparison with other in situ techniques (Bond et al., 1999; Collaud Coen et al., 2010). This work aims to: 1. Optically characterize BB aerosol generated from African fuel 2. Correct the aethalometer at multiple wavelengths 3. Parametrize these correction factors. 4. Combusting the same fuels with two difference tube furnace systems assess the general applicability of this flexible means of producing BB aerosol. This was accomplished by conducting laboratory studies at Los Alamos National Laboratory (LANL) and North Carolina A&T State University (NCA&T), where aethalometer measurements were compared to in situ optical measurements. Mass cross sections of BB aerosol from these under-studied fuels will be determined along with their wavelength dependence.

09 BIOMASS FUELS↗
DOE OSTI2023

Novel Calibration Approach for Monitoring Aerosol Hydrogen Isotopes Using Laser-Induced Breakdown Spectroscopy for Molten Salt Reactor Off-Gas Streams

Online monitoring is a key challenge for the continued development of molten salt reactor (MSR) technology. Laser-induced breakdown spectroscopy (LIBS) has previously been demonstrated to be a viable tool for monitoring aerosolized species and noble gases in real time, but the ability to discern varying isotopes in these streams has not yet been investigated for MSR applications. Tritium will form in MSRs from ternary fission and from (n,α)-reactions occurring in lithium-containing salts. This study compares three spectrometers of varying resolutions and types for measuring hydrogen isotope shifts in LIBS spectra of wetted filter paper. For each spectrometer, multivariate models were built (i.e., principal component regression, partial least squares regression, and multivariate curve resolution) to quantify the isotope ratio. The top models were then modified and corrected to apply the models to aerosol samples with varying isotope ratios. This novel calibration strategy offers an 82% reduction in volume of the calibration samples needed and is a more viable pathway for calibrating deployable LIBS systems. Lastly, this calibration model was compared with an all-aerosol trained model for monitoring hydrogen isotopes during a real-time test where the protium/deuterium ratio, along with representative salt species (i.e., lithium, sodium, and potassium) were adjusted dynamically. Results of this test validated the predictive capabilities of the transferred model and highlighted the capabilities of LIBS for real-time monitoring of MSR effluent streams.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Calculating Potential Radiological Emissions for Waste Management Activities at INL - 20068

At Idaho National Laboratory (INL), work involving radionuclides is evaluated for potential emissions from a project in order to comply with the National Emission Standards for Hazardous Air Pollutants (NESHAP) regulations, 40 CFR 61 Subpart H. Emission calculations are documented in an Air Permitting Applicability Determination (APAD) to analyze unmitigated and mitigated emissions and determine if an Application to Construct (ATC) or continuous monitoring is required. To calculate the unmitigated and mitigated emissions, a spreadsheet was developed to provide ease in determining potential emissions by providing the maximum operating temperature and the material being used. The spreadsheet aids in determining the potential emissions for research projects and waste management activities at Materials and Fuels Complex (MFC) and other locations across the INL site. Furthermore, it can also be used for periodic confirmatory measurements (PCM) to justify low emissions. Elements that factor into the unmitigated and mitigated calculations include the amount of each radionuclide used (in curies or grams), specific activity (if amount is given in grams), the temperature the material is heated to in Celsius, the dose conversion factor which is derived from Clean Air Act Assessment Package - 1988 (CAP-88) modeling, and the number of HEPA filters used for mitigated measures. The main drivers for calculating the unmitigated emissions for a project are the amount used per radionuclide, the maximum operating temperature, and the location of the work. The maximum operating temperature determines the airborne release factor which is dependent on the physical state of the radionuclide. Prior to October 2017, if the radionuclide was heated to greater than 100 deg. C, the radionuclide was assumed to be a gas, which has the highest airborne release factor. This assumption would be overly conservative for radionuclides with high melting and boiling points, which provided a challenge to demonstrate low emissions. In October 2017, the Environmental Protection Agency (EPA) approved an alternative method for INL. This method allows the airborne release factor to be determined by using the melting point and 90% of the boiling point of the radionuclide. This methodology was included in the spreadsheet to allow unmitigated emission calculations for APADs to be completed more efficiently and effectively. Results show a reduction in time completing air emission calculations as well as lower total emissions across all facilities at INL. MFC annual emissions were reduced by 62% from the previous year and Research and Education Campus (REC) facilities were reduced by 38% due to implementation of the approved alternative method. Time spent on APADs, PCMs, and documentation for the annual NESHAP report was also reduced significantly. The spreadsheet provided in Table I provides the potential emission calculations for the 'Advanced Retrieval and Disposition Techniques for Remote Handled Mixed Low Level Waste (RH MLLW) at the Radioactive Scrap and Waste Facility (RSWF)' project. Calculations show the Potential Effective Dose Equivalent (PEDE) at RSWF to be 7.27 E-04 mrem/yr (7.27 E-09 Sv/yr) which is well below the 0.1 mrem/yr threshold. (authors)

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

Neural networks for large eddy simulations of wall-bounded turbulence: numerical experiments and challenges

Here, we examine the application of neural network-based methods to improve the accuracy of large eddy simulations of incompressible turbulent flows. The networks are trained to learn a mapping between flow features and the subgrid scales, and applied locally and instantaneously—in the same way as traditional physics-based subgrid closures. Models that use only the local resolved strain rate are poorly correlated with the actual subgrid forces obtained from filtering direct numerical simulation data. We see that highly accurate models in a priori testing are inaccurate in forward calculations, owing to the preponderance of numerical errors in implicitly filtered large eddy simulations. A network that accounts for the discretization errors is trained and found to be unstable in a posteriori testing. We identify a number of challenges that the approach faces, including a distribution shift that affects networks that fail to account for numerical errors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Effect of particle size on the capture of uranium oxide colloidal particles from aqueous suspensions via high-gradient magnetic filtration

The effectiveness of High Gradient Magnetic Filtration (HGMF) in capturing uranium oxide particles from suspensions was investigated in this study. Two sets of experiments were performed to evaluate the importance of size on the capture of uranium oxide particles. The first considered two batches sieved into size bins of< 5, 5–10, 10–15, and 15–20 µm, while the second was performed using two suspensions with diameters smaller than 1.0 µm and between 1.0 and 1.5 µm. Iron oxide experiments, with particles between 0.3 and 0.8 µm, were performed for calibration purposes. In all experiments, a surfactant (Triton-X100 or sodium dodecyl sulfate) was used to prevent particle aggregation and limit the influence of non-magnetic capture mechanisms. A magnetic field of approximately 1.1 Tesla was generated using a water cooled electromagnet. HGMF was performed using tubular filters packed with ferromagnetic stainless-steel wool. Of the initial four uranium oxide particle sizes, magnetic capture was only observed for particles with a diameter of less than 5 µm, while larger particles experienced no magnetic and minimal total capture. For particles with diameters smaller than 1.0 µm and between 1.0 and 1.5 µm, capture efficiencies increased by 39 ± 9% and 34 ± 6% respectively, solely due to the magnetic field. Although the magnetic force is proportional to particle diameter, the capture efficiency decreased as diameter increased. So these results suggest that Brownian diffusion, which is influential for micron sized particles and increases with decreasing particle size, is acting in conjunction with the magnetic force to influence the efficacy of HGMF for uranium oxide. This important finding underscores the effectiveness of Brownian diffusion in increasing the rate of collision between particles and collector fibers. A stochastic trajectory model was developed to incorporate the influence of Brownian motion on particle behavior and filter removal efficiency. Modeling results are discussed and compared for uranium and iron oxide particles.

42 ENGINEERING↗
DOE OSTI2023

Combining Remotely Sensed Evapotranspiration and an Agroecosystem Model to Estimate Center‐Pivot Irrigation Water Use at High Spatio‐Temporal Resolution

Abstract Estimating irrigation water use accurately is critical for sustainable irrigation and studying terrestrial water cycle in irrigated croplands. However, irrigation is not monitored in most places, and current estimations of irrigation water use has coarse spatial and/or temporal resolutions. This study aims to estimate irrigation water use at the daily and field scale through the proposed model‐data fusion framework, which is achieved by particle filtering with two configurations (concurrent, CON, and sequential, SEQ) by assimilating satellite‐based evapotranspiration (ET) observations into an advanced agroecosystem model, ecosys . Two types of experiments using synthetic and real ET observations were conducted to study the efficacy of the proposed framework for estimating irrigation water use at the irrigated fields in eastern and western Nebraska, United States. The experiments using synthetic ET observations indicated that, for two major sources of uncertainties of ET difference between observations and model simulations, which are bias and noise, noise had larger impacts on degrading the estimation performance of irrigation water use than bias. For the experiments using real ET observations, monthly and annual estimations of irrigation water use matched well with farmer irrigation records, with Pearson correlation coefficient ( r ) around 0.80 and 0.50, respectively. Although detecting daily irrigation records was very challenging, our method still gave a good performance with RMSE, BIAS, and r around 2.90, 0.03, and 0.4 mm/d, respectively. Our proposed model‐data fusion framework for estimating irrigation water use at high spatio‐temporal resolution could contribute to regional water management, sustainable irrigation, and better tracking terrestrial water cycle.

Environmental Sciences & Ecology↗
DOE OSTI2023

Assessing the nexus between groundwater and solar energy plants in a desert basin with a dual-model approach under uncertainty

Globally, many solar power plants and other types of renewable energy are being located in water-scarce regions. Many projects rely on groundwater resources whose sustainability is uncertain. In the Chuckwalla Basin in California, quantification of recharge and trans-valley underflow is needed to estimate the impacts of solar project withdrawals on the water table. However, such estimates are highly challenging due to data scarcity, heterogeneous soils and long residence times. Conventional assessment employs isolated groundwater models configured with crude and uniform estimates of recharge. Here, we employ a data-constrained surface subsurface processes model, PAWS+CLM, to provide an ensemble of recharges and underflows with perturbed parameters. Then, the Parameter Estimation (PEST) package is used to calibrate MODFLOW aquifer conductivity and filter out implausible recharges. The novel dual-model approach, potentially applicable in other arid regions, can effectively assimilate groundwater head observations, reject unrealistic parameters, and narrow the range of estimated drawdowns. Simulated recharge concentrates along alluvial fans at the mountain fronts and ephemeral washes where run-off water infiltrates. If an evenly distributed recharge was assumed, it resulted in under-estimated drawdown and larger uncertainty bounds. The withdrawals are approaching total inflow, suggesting the system will be nearing, if not exceeding, its sustainable groundwater production capacity, and a boom of such projects will not be sustainable. Especially, the cost/benefit of pumped-storage projects is called into question as the initial-fill phase depletes entire area’s recharge. Our study highlights the stress on groundwater resources of solar development, and that the speed of groundwater recovery does not indicate sustainability. Main point 1: A novel dual model approach, involving an integrated surface/subsurface model and a groundwater parameter-estimation model, was able to better constrain the model. Main point 2: The groundwater system may be nearing, if not exceeding, its sustainable groundwater production capacity and the speed of recovery is not indicative of sustainability. Main point 3: Results from using conventionally-assumed uniform recharge distort calibrated K fields and impacts assessment

14 SOLAR ENERGY↗
DOE OSTI2021

Decentralized Filtering Adaptive Neural Network Control for Uncertain Switched Interconnected Nonlinear Systems

This article presents a novel decentralized filtering adaptive neural network control framework for uncertain switched interconnected nonlinear systems. Each subsystem has its own decentralized controller based on the established decentralized state predictor. For each subsystem, the nonlinear uncertainties are approximated by a Gaussian radial basis function (GRBF) neural network incorporated with a piecewise constant adaptive law, where the adaptive law will update adaptive parameters from the error dynamics between the host system and the decentralized state predictor by discarding the unknowns, whereas a decentralized filtering control law is derived to cancel both local and mismatched uncertainties from other subsystems, as well as achieve the local objective tracking of the host system. The achievement of global objective depends on the achievement of local objective for each subsystem. The matched uncertainties are canceled directly by adopting their opposite in the control signal, whereas a dynamic inversion of the system is required to eliminate the effect of the mismatched uncertainties on the output. By exploiting the average dwell time principle, the error bounds between the real system and the virtual reference system, which defines the best performance that can be achieved by the closed-loop system, are derived. A numerical example is given to illustrate the effectiveness of the decentralized filtering adaptive neural network control architecture by comparing against the model reference adaptive control (MRAC).

Average dwell time, decentralized, filtering adapt↗
DOE OSTI2025

Moltensaltpropnet

MoltenSaltPropnet is a physics-informed machine learning framework that aims to predict the thermophysical properties of molten fluoride and chloride salt mixtures, which are crucial for the design and safety of Generation IV molten salt reactors. The code processes data from the Molten-Salt Thermal Properties Database (MSTDB-TP) and the Janz compendium, converting critically evaluated correlations into fast, differentiable surrogate models for density, viscosity, thermal conductivity, and heat capacity across 448 distinct salt systems. The implementation consists of several key components: 1. Data Curation: The code parses and cleans the raw data, normalizing elemental mole fractions and extracting relevant regression coefficients for various thermophysical properties. 2. Feature Engineering: It generates fixed-length numerical descriptors that encapsulate the composition and temperature, incorporating polynomial interaction terms and dimensionality-reduction techniques to optimize model performance. 3. Coefficient Learning: Four different machine learning architectures are employed: a deep residual network (ResNet), a Kolmogorov–Arnold network (KAN), a sparsity-inducing neural network (SNN), and classical regression models. Each model learns to predict coefficients that define the temperature-dependent correlations for the thermophysical properties. 4. Property Reconstruction: The predicted coefficients are used to compute temperature-dependent property values, ensuring positivity and monotonic trends through a composite loss function that enforces physical constraints. 5. User Interface: An open-source web application enables users to filter the database, train task-specific models, and visualize the results, allowing for rapid exploration of candidate salt mixtures. MoltenSaltPropnet bridges the gap between limited experimental data and high-fidelity reactor simulations, providing a powerful tool for researchers in the field of molten salt reactors and advanced nuclear energy systems.

Retamales, Mauricio Eduardo Tano [Idaho National L↗
DOE OSTI2023

Quantifying airborne fraction trends and the destination of anthropogenic CO 2 by tracking carbon flows in a simple climate model

Abstract Atmospheric carbon dioxide (CO 2 ) concentrations have increased as a direct result of human activity and are at their highest level over the last 2 million years, with profound impacts on the Earth system. However, the magnitude and future dynamics of land and ocean carbon sinks are not well understood; therefore, the amount of anthropogenic fossil fuel emissions that remain in the atmosphere (the airborne fraction) is poorly constrained. This work aims to quantify the sources and controls of atmospheric CO 2 , the fate of anthropogenic CO 2 over time, and the likelihood of a trend in the airborne fraction. We use Hector v3.0, a coupled simple climate and carbon cycle model with the novel ability to explicitly track carbon as it flows through the Earth system. We use key model parameters in a Monte Carlo analysis of 15 000 model runs from 1750 to 2300. Results are filtered for physical realism against historical observations and CMIP6 projection data, and we calculate the relative importance of parameters controlling how much anthropogenic carbon ends up in the atmosphere. Modeled airborne fraction was roughly 52%, consistent with observational studies. The overwhelming majority of model runs exhibited a negative trend in the airborne fraction from 1960–2020, implying that current-day land and ocean sinks are proportionally taking up more carbon than the atmosphere. However, the percentage of atmospheric CO 2 derived from anthropogenic origins can be much higher because of Earth system feedbacks. We find it peaks at over 90% between 2010–2050. Moreover, when looking at the destination of anthropogenic fossil fuel emissions, only a quarter ends up in the atmosphere while more than half of emissions are taken up by the land sink on centennial timescales. This study evaluates the likelihood of airborne fraction trends and provides insights into the dynamics of anthropogenic CO 2 in the Earth system.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Filtering adaptive output feedback control for multivariable nonlinear systems with mismatched uncertainties and unmodeled dynamics

This article synthesizes a filtering adaptive output feedback controller for multivariable nonlinear systems with mismatched uncertainties and unmodeled dynamics. The multivariable nonlinear systems under consideration have both matched and mismatched uncertainties, which satisfy the semiglobal Lipschitz condition. The unmodeled dynamics are bounded-input bounded-output stable. By adopting an estimation/cancellation strategy, a piecewise constant adaptive law drives the estimation error to zero at every time instant, which yields the adaptive parameters; a disturbance rejection control law is designed to compensate the nonlinear uncertainties within the bandwidth of low-pass filters. The matched uncertainties are cancelled directly by adopting their opposite in the control signal, while a dynamic inversion of the system is required to eliminate the effect of the mismatched uncertainties on the output. A feedforward control law is designed to track a given command. Since the virtual reference system defines the best performance that can be achieved by the closed-loop system, the uniform performance bounds are derived for the states and control signals via comparison. Both numerical and practical examples are provided to illustrate the effectiveness of the proposed filtering adaptive output feedback control architecture, comparisons with the model reference adaptive control demonstrates the superiority of the proposed control method.

filtering adaptive control, mismatched uncertainit↗
DOE OSTI2025

Performance evaluation of automated data-driven feature extraction and selection methods for practical and scalable building energy consumption prediction models

Here, this study quantifies the impact of automated feature engineering methods (feature extraction and selection) on the quality and accuracy of machine learning models that predict building energy consumption. The case study compares model performance for three main scenarios: baseline (no feature extraction and selection), feature extraction only, and feature extraction combined with feature selection (filter and/or wrapper methods) for fully trained machine learning models for 200 metered/sub-metered energy measurements across 118 real buildings. For consistency, the same machine learning model architecture (a black box deep learning neural network with probabilistic forecast output) was used for all scenarios. Based on results, all feature engineering methods provided noticeable prediction accuracy improvements (e.g., 29%-68% median prediction improvement) compared to baseline scenarios. However, in this application, feature selection methods provide little practical value due to their limited performance gains and high computational cost. Smarter algorithm development supported by better computational environments will be needed before feature selection methods can reliably and efficiently improve predictive model performance.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Estimating mass-absorption cross-section of ambient black carbon aerosols: theoretical, empirical, and machine learning models

The mass-absorption cross-section of black carbon (MAC BC ) is an essential parameter to link the atmospheric concentration of black carbon (BC) with its radiative forcing. When a direct calculation of MAC BC based on observations of aerosol light absorption and BC mass concentration is impossible, we rely on modeling and simulations to estimate MAC BC , but currently, there is no consensus model that can be relied on for accurate predictions across all atmospheric environments when BC particles have different coating thicknesses. Here, we applied five MAC BC prediction models (including three light scattering theories, an empirical model based on observations of particle mass concentrations, and a machine learning model developed in our previous work) to aerosols from three Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) field campaigns. While many studies have found that increasing the complexity of the models helps to constrain biases of the estimated MAC BC , our effort is to evaluate the models based on the criteria of simplicity and accuracy. We find that our machine learning model (support vector machine for regression, SVM) generally performs well across all DOE ARM field campaign data, while the accuracy of core-shell Mie theory depends on the bias correction algorithm applied to filter-based light absorption data. Generally, the empirical model for internally-mixed particles that we considered tends to over-predict MAC BC , while Mie theory for externally-mixed particles tends to under-predict MAC BC . An examination of the influence of coating material on BC cores suggests that the performance of our current SVM model is degraded when the BC is thickly-coated (e.g., it has undergone aging and mixing with other materials in the atmosphere).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Virtual Dynamic Grid Impedance and Its Impacts on Harmonics and Stability of Inverter Based Resources Plant

Nowadays, a large number of inverter-based resources (IBRs) are integrated into the grid at a single connection point as an IBR plant. In this article, a virtual dynamic grid impedance concept is introduced to evaluate the harmonics and stability for grid integration of an IBR plant containing multiple IBRs. First, a detailed theoretical study is conducted to build a foundation of the virtual grid impedance concept, which is a dynamic impedance that changes with the number of IBRs added into an IBR plant. Then, based on the new virtual dynamic grid impedance concept, frequency spectrum analysis is performed to explore harmonic impact at the interconnection point for an IBR plant with IBRs having L, LC, and LCL filters, respectively. An electromagnetic transient (EMT) simulation model of a grid-connected IBR plant is developed to explore the harmonics and stability of the IBR plant connected to the grid as well as the reliable operation of IBRs within the plant from the novel virtual equivalent dynamic grid impedance point of view. Hardware experiments are conducted to validate the EMT simulation evaluation. The results show that the number of IBRs added into an IBR plant influences the grid impedance, i.e., grid strength, and the grid impedance variation has multiple impacts on the IBR plant and IBRs within the plant depending on the grid-connected filters of the IBRs.

grid impedance↗
DOE OSTI2026

A generative machine learning model for designing metal hydrides applied to hydrogen storage

Developing new metal hydrides is a critical step toward efficient hydrogen storage in carbon-neutral energy systems. However, existing materials databases, such as the Materials Project, contain a limited number of well-characterized hydrides, which constrains the discovery of optimal candidates. This work presents a framework that integrates causal discovery with a lightweight generative machine learning model to generate novel metal hydride candidates that may not exist in current databases. Using a dataset of 450 samples (270 training, 90 validation, and 90 testing), the model generates 1000 candidates. After ranking and filtering, six previously unreported chemical formulas and crystal structures are identified, four of which are validated by density functional theory simulations and show strong potential for future experimental investigation. Overall, the proposed framework provides a scalable and time-efficient approach for expanding hydrogen storage datasets and accelerating materials discovery.

generative model↗
DOE OSTI2024

Probing the evolution of fault properties during the seismic cycle with deep learning

We use seismic waves that pass through the hypocentral region of the 2016 M6.5 Norcia earthquake together with Deep Learning (DL) to distinguish between foreshocks, aftershocks and time-to-failure (TTF). Binary and N-class models defined by TTF correctly identify seismograms in test with > 90% accuracy. We use raw seismic records as input to a 7 layer CNN model to perform the classification. Here we show that DL models successfully distinguish seismic waves pre/post mainshock in accord with lab and theoretical expectations of progressive changes in crack density prior to abrupt change at failure and gradual postseismic recovery. Performance is lower for band-pass filtered seismograms (below 10 Hz) suggesting that DL models learn from the evolution of subtle changes in elastic wave attenuation. Tests to verify that our results indeed provide a proxy for fault properties included DL models trained with the wrong mainshock time and those using seismic waves far from the Norcia mainshock; both show degraded performance. Our results demonstrate that DL models have the potential to track the evolution of fault zone properties during the seismic cycle. If this result is generalizable it could improve earthquake early warning and seismic hazard analysis.

58 GEOSCIENCES↗
DOE OSTI2021

Real-time ensemble microalgae growth forecasting with data assimilation

Accurate short-range (e.g., 7-day) microalgae growth forecasts will be beneficial for both production and harvesting of microalgae. This study developed an operational microalgae growth forecasting system with ensemble data assimilation (DA). The forecasting system was validated against observed Monoraphidium minutum 26B-AM growth in two outdoor pond cultures located in Mesa, Arizona, U.S. We first examined the relative roles of uncertainty in the meteorological forecast and initial conditions (i.e., algal concentration at the time of forecast) in the microalgae 7-day forecast and found initial conditions dominated the microalgae forecasting skill, suggesting the importance of implementing DA to improve initial condition characterization. To correct the systematic bias in biomass simulations, we developed a particle filter with bias estimation (PFBE) DA method to estimate biases and correct the model forecast. We found the DA forecasting system could improve the 7-day microalgae forecasting skill by about 85% on average compared to model forecasts without DA. These results suggest the potential accuracy of biomass growth forecasts may be sufficient to inform real-time operational decisions, such as harvesting planning, for commercial-scale microalgae production.

59 BASIC BIOLOGICAL SCIENCES↗