Search NASA⌕ Search

SEARCH · Search NASA

Results for “DATA ASSIMILATION”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
NASA NTRS2018

Benefits and Pitfalls of GRACE Terrestrial Water Storage Data Assimilation

Satellite observations of terrestrial water storage (TWS) from the Gravity Recovery and Climate Experiment (GRACE) mission have a coarse resolution in time (monthly) and space (roughly 150,000 sq km at midlatitudes) and vertically integrate all water storage components over land, including soil moisture and groundwater. Nonetheless, data assimilation can be used to horizontally downscale and vertically partition GRACE-TWS observations. This presentation illustrates some of the benefits and drawbacks of assimilating TWS observations from GRACE into a land surface model over the continental United States and India. The assimilation scheme yields improved skill metrics for groundwater compared to the no-assimilation simulations. A smaller impact is seen for surface and root-zone soil moisture. Further, GRACE observes TWS depletion associated with anthropogenic groundwater extraction. Results from the assimilation emphasize the importance of representing anthropogenic processes in land surface modeling and data assimilation systems.

TWS↗
DOE OSTI2025

Data Assimilation with Machine Learning Surrogate Models: A Case Study with FourCastNet

Modern data-driven surrogate models for weather forecasting provide accurate short-term predictions but inaccurate and nonphysical long-term forecasts. This paper investigates online weather prediction using machine learning surrogates supplemented with partial and noisy observations. We empirically demonstrate and theoretically justify that, despite the long-time instability of the surrogates and the sparsity of the observations, filtering estimates can remain accurate in the long-time horizon. As a case study, we integrate FourCastNet, a weather surrogate model, within a variational data assimilation framework using partial, noisy ERA5 data. Our results show that filtering estimates remain accurate over a year-long assimilation window and provide effective initial conditions for forecasting tasks, including extreme event prediction.

Adrian, Melissa [Univ. of Chicago, IL (United Stat↗
NASA NTRS2001

Evaluation of Bogus Vortex Techniques Using Four-Dimensional Variational Data Assimilation

The effectiveness of a four-dimensional variational data assimilation (4-DVAR) technique for creating "bogus" vortices in numerical simulations of hurricanes is evaluated in this study. A series of numerical experiments is conducted to generate initial vortices for Hurricane Georges and Bonnie (1998) in the Atlantic Ocean by assimilating bogus sea-level pressure and wind information into a mesoscale numerical model (MM5). Several different strategies are tested for investigating the sensitivity of the initial vortex representation to the type of bogus information.

Pu, Zhao-Xia↗
DOE OSTI2022

Physics‐Informed Machine Learning Method for Large‐Scale Data Assimilation Problems

Abstract We develop a physics‐informed machine learning approach for large‐scale data assimilation and parameter estimation and apply it for estimating transmissivity and hydraulic head in the two‐dimensional steady‐state subsurface flow model of the Hanford Site given synthetic measurements of said variables. In our approach, we extend the physics‐informed conditional Karhunen‐Loéve expansion (PICKLE) method to modeling subsurface flow with unknown flux (Neumann) and varying head (time‐dependent Dirichlet) boundary conditions. We demonstrate that the PICKLE method is comparable in accuracy with the standard maximum a posteriori (MAP) method, but is significantly faster than MAP for large‐scale problems. Both methods use a mesh to discretize the computational domain. In MAP, the parameters and states are discretized on the mesh; therefore, the size of the MAP parameter estimation problem directly depends on the mesh size. In PICKLE, the mesh is used to evaluate the residuals of the governing equation, while the parameters and states are approximated by the truncated conditional Karhunen‐Loéve expansions with the number of parameters controlled by the smoothness of the parameter and state fields, and not by the mesh size. For a considered example, we demonstrate that the computational cost of PICKLE increases near linearly (as N 1.15 ) with the number of grid nodes N , while that of MAP increases much faster (as N 3.28 ). We also show that once trained for one set of Dirichlet boundary conditions (i.e., one river stage), the PICKLE method provides accurate estimates of the hydraulic head for any value of the Dirichlet boundary conditions (i.e., for any river stage).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Adjoint sensitivity analysis and data assimilation for verification of dry storage cask contents

Dry cask storage is a method for interim storage of spent fuel assemblies which contain fissile isotopes of uranium and plutonium. These can present a proliferation concern and consequently there is a need for non-destructive testing methods to verify a dry cask's contents for proliferation protection. We present an application of adjoint sensitivity analysis and data assimilation to a multigroup diffusion model of dry cask storage. Adjoint sensitivity analysis allows the efficient calculation of sensitivities for use in data assimilation to calibrate imprecisely known parameter values and data consistency tests to detect diversion scenarios. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Data Assimilative Optimization of WSA Source Surface and Interface Radii using Particle Filtering

TheWang-Sheeley-Arge (WSA) model estimates solar wind speed and interplanetary magnetic field polarity in the inner heliosphere using global photospheric magnetic field maps. WSA employs the Potential Field Source Surface (PFSS) and Schatten Current Sheet (SCS) models to determine the Sun's global coronal magnetic field. The PFSS and SCS models are connected through two radial parameters, the source surface and interface radii, which specify the overlap region between the inner SCS and outer PFSS models. Though both radii values are adjustable, they have typically been fixed to 2.5 solar radii. Our work highlights how solar wind predictions improve when the radii are allowed to vary over time. Data assimilation using particle filtering (sequential Monte Carlo) is used to infer optimal values over a fixed time window. Solar wind model predictions and satellite observations are compared with a newly developed quality-of-agreement prediction metric. The agreement metric between the model and observations is assumed to correspond to the probability of the two key WSA model parameters, the source surface and interface radii, where the highest metric value implies the optimal radii.We find that the optimal particle filter values of solar radii can perform twice as well as standard values for an exploratory period during Carrington Rotation 1901, with these values also reducing nonphysical kinking effects seen in solar magnetic field lines. Data assimilation choices of input realization and time frame have implications for variation in the solar wind over time. We present this work's theoretical context and practical applications for prediction accuracy.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2002

A Variable-Resolution Stretched-Grid General Circulation Model and Data Assimilation System with Multiple Areas of Interest: Studying the Anomalous Regional Climate Events of 1998

The new stretched-grid design with multiple (four) areas of interest, one at each global quadrant, is implemented into both a stretched-grid GCM (general circulation model) and a stretched-grid data assimilation system (DAS). The four areas of interest include: the U.S./Northern Mexico, the El Nino area/Central South America, India/China, and the Eastern Indian Ocean/Australia. Both the stretched-grid GCM and DAS annual (November 1997 through December 1998) integrations are performed with 50 km regional resolution. The efficient regional down-scaling to mesoscales is obtained for each of the four areas of interest while the consistent interactions between regional and global scales and the high quality of global circulation, are preserved. This is the advantage of the stretched-grid approach. The global variable resolution DAS incorporating the stretched-grid GCM has been developed and tested as an efficient tool for producing regional analyses and diagnostics with enhanced mesoscale resolution. The anomalous regional climate events of 1998 that occurred over the U.S., Mexico, South America, China, India, African Sahel, and Australia are investigated in both simulation and data assimilation modes. Tree assimilated products are also used, along with gauge precipitation data, for validating the simulation results. The obtained results show that the stretched-grid GCM and DAS are capable of producing realistic high quality simulated and assimilated products at mesoscale resolution for regional climate studies and applications.

Fox-Rabinovitz, Michael S.↗
NASA NTRS2008

Data Assimilation Cycling for Weather Analysis

This software package runs the atmospheric model MM5 in data assimilation cycling mode to produce an optimized weather analysis, including the ability to insert or adjust a hurricane vortex. The program runs MM5 through a cycle of short forecasts every three hours where the vortex is adjusted to match the observed hurricane location and storm intensity. This technique adjusts the surrounding environment so that the proper steering current and environmental shear are achieved. MM5cycle uses a Cressman analysis to blend observation into model fields to get a more accurate weather analysis. Quality control of observations is also done in every cycle to remove bad data that may contaminate the analysis. This technique can assimilate and propagate data in time from intermittent and infrequent observations while maintaining the atmospheric field in a dynamically balanced state. The software consists of a C-shell script (MM5cycle.driver) and three FORTRAN programs (splitMM5files.F, comRegrid.F, and insert_vortex.F), and are contained in the pre-processor component of MM5 called "Regridder." The model is first initialized with data from a global model such as the Global Forecast System (GFS), which also provides lateral boundary conditions. These data are separated into single-time files using splitMM5.F. The hurricane vortex is then bogussed in the correct location and with the correct wind field using insert_vortex.F. The modified initial and boundary conditions are then recombined into the model fields using comRegrid.F. The model then makes a three-hour forecast. The three-hour forecast data from MM5 now become the analysis for the next short forecast run, where the vortex will again be adjusted. The process repeats itself until the desired time of analysis is achieved. This code can also assimilate observations if desired.

Tran, Nam↗
NASA NTRS2018

Atmosphere-Ocean Coupled Data Assimilation Using NASA GEOS: Estimation of Air-Sea Interface State Variables

Air-sea interface variables, such as the skin Sea Surface Temperature (SST) are essential for atmosphere-ocean coupling. In the NASA GMAO Data Assimilation System (DAS), the skin SST and 3-D atmospheric state are jointly estimated [1]. This presentation is focused on the prior or background error covariance that is used in this analysis. The GEOS DAS uses an ensemble-variational assimilation strategy. In that, specification of a climatological background (CB) error covariance for SST relies on the NOAA's OI SST, with estimates of standard deviation and correlation length scales based on weekly analyses of the bulk SST at 1 degree resolution. However, present analysis system is striving to resolve SST diurnal variability with six hourly analyses and assimilates a vast number of in-situ and satellite observations. The first part of this presentation re-derives the CB error covariance using OSTIA SST analyses and illustrates the impact of this update on assimilating satellite observations. In a hybrid assimilation system the CB error covariances are appended with a flow-dependent background error covariance estimate implied by the underlying ensemble. The second part of this presentation refers to: a. treatment of the skin SST in the ensemble members, b. corresponding ensemble spread, and c. impact of these additions on the data assimilation system. [1] S. Akella, et al. (2017), doi:10.1002/qj.2988

Akella, Santha↗
NASA NTRS1999

Precipitation Recycling in the NASA GEOS Data Assimilation System

Analysis of precipitation recycling can improve the understanding of regional hydrologic anomalies, especially their evolution and maintenance. Diagnostic models of the recycling of precipitation and are applied to 15 years of the NASA Goddard Earth Observing System (GEOS) Data Assimilation System (DAS). Recycled precipitation is defined as the fraction of precipitation within a given region that originated as surface evaporation from the same region. The focus of the present work is on the interannual variability of the central United States hydrologic cycle and precipitation recycling. The extreme years of 1988 (drought) and 1993 (flood) are compared with the 15 year base period mean annual cycle. The results indicate that recycling ratio (the amount of precipitation with a local source relative to the total precipitation) is greater in 1988 than both the base period mean and the 1993 season (with 1993 recycling ratio less than the mean). On the other hand, both the summers of 1988 and 1993 show less total recycled precipitation than the mean. The results also show that precipitation recycling may have been more important in the spring of 1993, when the region was primed for flooding, than the summer, when the sever flooding occurred. The diagnostic approaches to precipitation recycling suffer from some weaknesses. Numerical simulations and assimilation using passive tracers have the potential to provide more accurate calculations of precipitation recycling and the remote sources of water. This ability is being incorporated into the latest GEOS data assimilation system, and some preliminary results will be presented.

Bosilovich, Michael G.↗
DOE OSTI2024

ML-based Data Assimilation and History Matching: Application to the IBDP CCS Project

It is crucial to monitor the CO2 plume effectively throughout the life cycle of a geologic CO2 sequestration project to ensure safety and storage efficiency. However, the computational cost of existing data assimilation methods can be prohibitively expensive due to the complex physics with multi-component non-isothermal simulation and high dimensionality of large-scale reservoir models. We address this challenge by proposing an accelerated deep learning-based workflow for model calibration and prediction of CO2 plume evolution in the reservoir.The power and efficacy of our workflow is demonstrated by application to the Illinois Basin-Decatur Project (IBDP), a large-scale CO2 storage test in saline aquifer. The data assimilation process is implemented rapidly by the proposed workflow with given field measurements including distributed pressure and temperature sensing (DTS) data at an injection and a monitoring well. CO2 plume evolution is predicted by running the simulations of the calibrated reservoir models.

Nagao, Masahiro↗
NASA NTRS1999

Data Assimilation Into a Coupled Ocean Atmosphere Model: Application to the 1997-1998 El Nino

As part of JPL's ocean data assimilation effort to study ocean circulation and seasonal-interannual climate variability, sea level anomaly observed by TOPEX altimeter, together with sea surface temperature and wind stress data, are assimilated into a simple coupled ocean atmosphere model of the tropical Pacific. Model-data consistency is examined. Impact of the assimilation (as initialization) on El Nino Southern Oscillation (ENSO) forecasts is evaluated. The coupled model consists of a shallow water component with two baroclinic modes, an Ekman shear layer, a simplified mixed-layer temperature equation, and a statistical atmosphere based on dominant correlations between historical surface temperature and wind stress anomaly data. The adjoins method is used to fit the coupled model to the data over various six-month periods from late 1996 to early 1998 by optimally adjusting the initial state, model parameters, and basis functions of the statistical atmosphere. On average, the coupled model can be fitted to the data to approximately within the data and representation errors (5 cm, 0.5 C, and 10 sq m/sq m for sea level, surface temperature, and pseudo wind stress anomalies, respectively). The estimated fields resemble observed spatio-temporal structure reasonably well. Hindcasts/forecasts of the 1997/1998 El Nino initialized from forced estimated ocean states and parameters are much more realistic than those simply initialized from ocean states (see figure below). In particular, the ability of the model to produce significant warming beyond the initial state is dramatically improved. Parameter estimation, which compensates for some model errors, is found to be important to obtaining better fits of the model to data and to improving forecasts.

Lee, Tong↗
NASA NTRS2015

Technical Report Series on Global Modeling and Data Assimilation, Volume 43: Initial Evaluation of the Climate - MERRA-2

The years since the introduction of MERRA have seen numerous advances in the GEOS-5 Data Assimilation System as well as a substantial decrease in the number of observations that can be assimilated into the MERRA system. To allow continued data processing into the future, and to take advantage of several important innovations that could improve system performance, a decision was made to produce MERRA-2, an updated retrospective analysis of the full modern satellite era. One of the many advances in MERRA-2 is a constraint on the global dry mass balance; this allows the global changes in water by the analysis increment to be near zero, thereby minimizing abrupt global interannual variations due to changes in the observing system. In addition, MERRA-2 includes the assimilation of interactive aerosols into the system, a feature of the Earth system absent from previous reanalyses. Also, in an effort to improve land surface hydrology, observations-corrected precipitation forcing is used instead of model-generated precipitation. Overall, MERRA-2 takes advantage of numerous updates to the global modeling and data assimilation system. In this document, we summarize an initial evaluation of the climate in MERRA-2, from the surface to the stratosphere and from the tropics to the poles. Strengths and weaknesses of the MERRA-2 climate are accordingly emphasized.

MERRA-2↗
NASA NTRS2005

Documentation and Validation of the Goddard Earth Observing System (GEOS) Data Assimilation System, Version 4

This document describes the structure and validation of a frozen version of the Goddard Earth Observing System Data Assimilation System (GEOS DAS): GEOS-4.0.3. Significant features of GEOS-4 include: version 3 of the Community Climate Model (CCM3) with the addition of a finite volume dynamical core; version two of the Community Land Model (CLM2); the Physical-space Statistical Analysis System (PSAS); and an interactive retrieval system (iRET) for assimilating TOVS radiance data. Upon completion of the GEOS-4 validation in December 2003, GEOS-4 became operational on 15 January 2004. Products from GEOS-4 have been used in supporting field campaigns and for reprocessing several years of data for CERES.

Suarez, Max J.↗
NASA NTRS2018

The Atmosphere-Ocean Interface Layer of NASA's Goddard Earth Observing System Model and Data Assimilation System Volume 51

The Goddard Earth Observing System (GEOS) general circulation model (GCM) includes modules for sea surface temperature (SST) diurnal warming and cool-skin layers. To support the application of a coupled atmosphere-ocean data assimilation capability, the GCM needs to be flexible enough to support both coupled atmosphere ocean general circulation model (AOGCM) and atmosphere-only (AGCM) configurations, with only minor configuration changes at the user interface. This document presents a formulation of an atmosphere-ocean interface layer (AOIL) that serves this purpose. Previous work by Akella et al. (2017) described a version of a model for near-surface temperature variations, including both both diurnal warming and cool-skin effects, that has been used since 2017 in the near-real-time GEOS FP (forward processing) weather analysis and forecasting system. The diurnal cycle of SST in that version of the GEOS atmospheric data assimilation system (ADAS) undergoes a sharp decay in the late afternoon (local time). The updated AOIL presented here includes a modification of the similarity function used in the diurnal warming model. Results from offline model runs illustrate an improvement in the near-surface (less than 0:5m depth) diurnal cycle compared to the original formulation. The new formulation requires minimal parameter tuning, and the improvements are robust across long (several month) simulation periods. This new model formulation, however, retains some deficiences from the previous module, such as a small warm bias in calm wind conditions for water depths below 1m. Our future work would include surface salinification and sea-ice into the AOIL.

Koster, Randal D.↗
NASA NTRS1995

Technical Report Series on Global Modeling and Data Assimilation. Volume 7: Proceedings of the Workshop on the GEOS-1 Five-year Assimilation

The primary objective of the three-day workshop on results from the Data Assimilation Office (DAO) five-year assimilation was to provide timely feedback from the data users concerning the strengths and weaknesses of version 1 of the Goddard Earth Observing System (GEOS-1) assimilated products. A second objective was to assess user satisfaction with the current methods of data access and retrieval. There were a total of 49 presentations, with about half (23) of the presentations from scientists from outside of Goddard. The first two days were devoted to applications of data: studies of the energy diagnostics, precipitation and diabatic heating, hydrological modeling and moisture transport, cloud forcing and validation, various aspects of intraseasonal, seasonal, and interannual variability, ocean wind stress applications, and validation of surface fluxes. The last day included talks from the National Meteorological Center (NMC), the National Center for Atmospheric Research (NCAR), the Center for Ocean-Land-Atmosphere Studies (COLA), the United States Navy, and the European Center for Medium Range Weather Forecasts (ECMWF).

Suarez, Max J.↗
NASA NTRS2015

Completing the Feedback Loop: The Impact of Chlorophyll Data Assimilation on the Ocean State

In anticipation of the integration of a full biochemical model into the next generation GMAO coupled system, an intermediate solution has been implemented to estimate the penetration depth (1Kd_PAR) of ocean radiation based on the chlorophyll concentration. The chlorophyll is modeled as a tracer with sources-sinks coming from the assimilation of MODIS chlorophyll data. Two experiments were conducted with the coupled ocean-atmosphere model. In the first, climatological values of Kpar were used. In the second, retrieved daily chlorophyll concentrations were assimilated and Kd_PAR was derived according to Morel et al (2007). No other data was assimilated to isolate the effects of the time-evolving chlorophyll field. The daily MODIS Kd_PAR product was used to validate the skill of the penetration depth estimation and the MERRA-OCEAN re-analysis was used as a benchmark to study the sensitivity of the upper ocean heat content and vertical temperature distribution to the chlorophyll input. In the experiment with daily chlorophyll data assimilation, the penetration depth was estimated more accurately, especially in the tropics. As a result, the temperature bias of the model was reduced. A notably robust albeit small (2-5 percent) improvement was found across the equatorial Pacific ocean, which is a critical region for seasonal to inter-annual prediction.

MERRA-OCEAN↗
NASA NTRS2018

A New Atmosphere-Ocean Model for Studying Air-Sea Interactions and Coupled Data Assimilation

During the last two plus decades, The Goddard Earth Observing System (GEOS) and Massachusetts Institute of Technology (MIT) modeling groups have developed, respectively, atmosphere-only and ocean-only global general circulation models. These two models (GEOS and MITgcm) have demonstrated their data assimilation capabilities with the recent releases of the Modern Era Reanalysis for Research Applications, Version 2 (MERRA-2) atmospheric reanalysis and the Estimating the Circulation and Climate of the Ocean, Version 4 (ECCO-v4) ocean (and sea ice) state estimate. Independently, the two modeling groups have also produced global atmosphere-only and ocean-only simulations with km-scale grid spacing which proved invaluable for process studies and for the development of satellite and in-situ sampling strategies.Recently, a new effort has been made to couple these two models and to leverage their data-assimilation and high resolution capabilities (i.e., eddy-permitting ocean, cloud-permitting atmosphere). The focus in the model development is put on sub-seasonal to decadal time scales. In this talk, I discuss the new coupled model and present some first coupled simulation results. This will include a high-resolution coupled GEOS-MIT simulation, whereby we have coupled a cubed-sphere-720 (~ 1/8°) configuration of the GEOS atmosphere to a lat-lon-cap-1080 (~ 1/12°) configuration of the MIT ocean. We compare near-surface diagnostics of this fully coupled ocean-atmosphere set-up to equivalent atmosphere-only and ocean-only simulations. In the comparisons we focus in particular on the differences in air-sea interactions between sea surface temperature (SST) and wind for the coupled and uncoupled simulations.

Strobach, Ehud↗