Search NASASearch

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 163 records · Page 9

Comparison of MERRA-2 and ECCO-V4 Ocean Surface Heat Fluxes: Consequences of Different Forcing Feedbacks on Ocean Circulation and Implications for Climate Data Assimilation

Forcing ocean models with reanalysis data is a common practice in ocean modeling. As part of this practice, prescribed atmospheric state variables and interactive ocean SST (Sea Surface Temperature) are used to calculate fluxes between the ocean and the atmosphere. When forcing an ocean model with reanalysis fields, errors in the reanalysis data, errors in the ocean model and errors in the forcing formulation will generate a different solution compared to other ocean reanalysis solutions (which also have their own errors). As a first step towards a consistent coupled ocean-atmosphere reanalysis, we compare surface heat fluxes from a state-of-the-art atmospheric reanalysis, the Modern-Era Retrospective analysis for Research and Applications, Version 2 (MERRA-2), to heat fluxes from a state-of-the-art oceanic reanalysis, the Estimating the Circulation and Climate of the Ocean Version 4, Release 2 (ECCO-v4). Then, we investigate the errors associated with the MITgcm (Massachusetts Institute of Technology general circulation model) ocean model in its ECCO-v4 ocean reanalysis configuration (1992-2011) when it is forced with MERRA- 2 atmospheric reanalysis fields instead of with the ECCO-v4 adjoint optimized ERA-interim state variables. This is done by forcing ECCO-v4 ocean with and without feedbacks from MERRA-2 related to turbulent fluxes of heat and moisture and the outgoing long wave radiation. In addition, we introduce an intermediate forcing method that includes only the feedback from the interactive outgoing long wave radiation. The resulting ocean circulation is compared with ECCO-v4 reanalysis and in-situ observations. We show that, without feedbacks, imbalances in the energy and the hydrological cycles of MERRA-2 (which are directly related to the fact it was created without interactive ocean) result in considerable SST drifts and a large reduction in sea level. The bulk formulae and interactive outgoing long wave radiation, although providing air-sea feedbacks and reducing model-data misfit, strongly relax the ocean to observed SST and may result in unwanted features such as large change in the water budget. These features have implications in a desired forcing recipe to be used. The results strongly and unambiguously argue for next generation data assimilation climate studies to involve fully coupled systems.

MERRA-2

Variational data assimilation with a semi-Lagrangian semi-implicit global shallow-water equation model and its adjoint

An adjoint model is developed for variational data assimilation using the 2D semi-Lagrangian semi-implicit (SLSI) shallow-water equation global model of Bates et al. with special attention being paid to the linearization of the interpolation routines. It is demonstrated that with larger time steps the limit of the validity of the tangent linear model will be curtailed due to the interpolations, especially in regions where sharp gradients in the interpolated variables coupled with strong advective wind occur, a synoptic situation common in the high latitudes. This effect is particularly evident near the pole in the Northern Hemisphere during the winter season. Variational data assimilation experiments of 'identical twin' type with observations available only at the end of the assimilation period perform well with this adjoint model. It is confirmed that the computational efficiency of the semi-Lagrangian scheme is preserved during the minimization process, related to the variational data assimilation procedure.

Li, Y.

VAS data assimilation into a mesoscale model

In numerical weather prediction, one of the arising problems has always been related to insufficient data, especially in cases involving scales smaller than the synoptic rawinsonde network. The Visible Infrared Spin-Scan Radiometer (VISSR) Atmospheric Sounder (VAS) instrument on the Geostationary Operational Environmental Satellite (GOES) can provide mesoscale vertical profiles of temperature and moisture several times over a twelve hour period. However, as this data is only available over a relatively small area at any one time, the VAS data is more suited to data assimilation. A variational method has been developed to assimilate VAS data into a mesoscale model. In the present paper, the assimilation method and VAS data from one case are discussed. It is found that the greatest improvements to the forecast resulted when the VAS data resolved the mesoscale structure of the temperature and relative humidity fields.

Cram, J. M.

Assessing the Impact of Pre-gpm Microwave Precipitation Observations in the Goddard WRF Ensemble Data Assimilation System

The forthcoming Global Precipitation Measurement (GPM) Mission will provide next generation precipitation observations from a constellation of satellites. Since precipitation by nature has large variability and low predictability at cloud-resolving scales, the impact of precipitation data on the skills of mesoscale numerical weather prediction (NWP) is largely affected by the characterization of background and observation errors and the representation of nonlinear cloud/precipitation physics in an NWP data assimilation system. We present a data impact study on the assimilation of precipitation-affected microwave (MW) radiances from a pre-GPM satellite constellation using the Goddard WRF Ensemble Data Assimilation System (Goddard WRF-EDAS). A series of assimilation experiments are carried out in a Weather Research Forecast (WRF) model domain of 9 km resolution in western Europe. Sensitivities to observation error specifications, background error covariance estimated from ensemble forecasts with different ensemble sizes, and MW channel selections are examined through single-observation assimilation experiments. An empirical bias correction for precipitation-affected MW radiances is developed based on the statistics of radiance innovations in rainy areas. The data impact is assessed by full data assimilation cycling experiments for a storm event that occurred in France in September 2010. Results show that the assimilation of MW precipitation observations from a satellite constellation mimicking GPM has a positive impact on the accumulated rain forecasts verified with surface radar rain estimates. The case-study on a convective storm also reveals that the accuracy of ensemble-based background error covariance is limited by sampling errors and model errors such as precipitation displacement and unresolved convective scale instability.

microwave observations

Recent Advances in Ozone Data Assimilation at the GMAO - Towards a New Reanalysis

This presentation summarized ongoing work on improving the representation of ozone in the GEOS Data Assimilation Systems. Data from two EOS Aura sensors was used: the total column ozone from the Ozone Monitoring Instrument (OMI) and high vertical resolution stratospheric profiles from Microwave Limb Sounder (MLS, version 3.3). As several previous studies have demonstrated, assimilation of this data can constrain the stratospheric and tropospheric ozone columns with relatively good accuracy. However, the representation of the vertical structures in the troposphere and near tropopause region is often deficient. Since both these layers of the atmosphere are critical to the understanding of the radiative forcing as well as the ozone budget in the troposphere, current work will focus on improving the assimilated product between the surface and the 50 hPa pressure level. The discussion included recent steps that have been taken towards refining the treatment of ozone in GEOS-5. Impacts of improved tropospheric chemistry model were discussed including the introduction of efficiency factors ("averaging kernels") for OMI total ozone, and direct assimilation of radiances from the MLS instrument. In particular, advantages and challenges involved in assimilating limb radiances rather than retrieved product were discussed. This work is, in part, a preparation for a planned reanalysis of the EOS Aura data from 2005 to present.

Krzysztof, Wargan

Results from NASA's Near Real Time Ozone Data Assimilation Systems

Organizations in Europe, Australia, and the United States have recently broadened constituent assimilation activities beyond water vapor, which has been assimilated for years in the numerical weather prediction applications. Many of these activities have focused on ozone, with some efforts focused on the entire suite of reactive constituents that control the ozone distribution. This talk will draw from results from the near real-time ozone data assimilation system being run by NASA's Data Assimilation Office. This system utilizes ozone observations from both the TOMS and the SBUV instrument to generate global synoptic maps of ozone. The initial application of this product is to provide ozone fields to assist in the atmospheric corrections' that are necessary for the retrieval of information from other NASA instruments. The validation of the ozone assimilation system shows that the assimilated product agrees well with independent HALOE and ozonesonde observations. Aside from providing a global synoptic map, there is verifiable geophysical information at higher vertical resolution than either of the date types input into the system. This talk will establish the validation results and enumerate applications of the ozone data assimilation system. Results from exploratory research will be presented. The applications being considered include estimates of tropospheric ozone, provision of ozone fields for interactive retrievals, use of analysis increments from the assimilation to evaluate model performance, and development of long-term consistent three-dimensional global ozone fields. The results from the exploratory studies are promising, and help demonstrate how assumptions made in the development of the ozone assimilation impact the other applications. For instance, RMS errors in the current product are large near the tropopause, which is sensitive to the specification of vertical correlation functions, which in turns impacts the amount of ozone analyzed to be in the troposphere. How these sensitivities impact the different applications will also be discussed.

Rood, Richard B.

Global 3-D Ionospheric Electron Density Reanalysis Based on Multisource Data Assimilation

We report preliminary results of a global 3-D ionospheric electron density reanalysis demonstration study during 2002-2011 based on multisource data assimilation. The monthly global ionospheric electron density reanalysis has been done by assimilating the quiet days ionospheric data into a data assimilation model constructed using the International Reference Ionosphere (IRI) 2007 model and a Kalman filter technique. These data include global navigation satellite system (GNSS) observations of ionospheric total electron content (TEC) from ground-based stations, ionospheric radio occultations by CHAMP, GRACE, COSMIC, SAC-C, Metop-A, and the TerraSAR-X satellites, and Jason-1 and 2 altimeter TEC measurements. The output of the reanalysis are 3-D gridded ionospheric electron densities with temporal and spatial resolutions of 1 h in universal time, 5deg in latitude, 10deg in longitude, and approx.30 km in altitude. The climatological features of the reanalysis results, such as solar activity dependence, seasonal variations, and the global morphology of the ionosphere, agree well with those in the empirical models and observations. The global electron content derived from the international GNSS service global ionospheric maps, the observed electron density profiles from the Poker Flat Incoherent Scatter Radar during 2007-2010, and foF2 observed by the global ionosonde network during 2002-2011 are used to validate the reanalysis method. All comparisons show that the reanalysis have smaller deviations and biases than the IRI-2007 predictions. Especially after April 2006 when the six COSMIC satellites were launched, the reanalysis shows significant improvement over the IRI predictions. The obvious overestimation of the low-latitude ionospheric F region densities by the IRI model during the 23/24 solar minimum is corrected well by the reanalysis. The potential application and improvements of the reanalysis are also discussed.

Xinan Yue

Sequential estimation and satellite data assimilation in meteorology and oceanography

The central theme of this review article is the role that dynamics plays in estimating the state of the atmosphere and of the ocean from incomplete and noisy data. Objective analysis and inverse methods represent an attempt at relying mostly on the data and minimizing the role of dynamics in the estimation. Four-dimensional data assimilation tries to balance properly the roles of dynamical and observational information. Sequential estimation is presented as the proper framework for understanding this balance, and the Kalman filter as the ideal, optimal procedure for data assimilation. The optimal filter computes forecast error covariances of a given atmospheric or oceanic model exactly, and hence data assimilation should be closely connected with predictability studies. This connection is described, and consequences drawn for currently active areas of the atmospheric and oceanic sciences, namely, mesoscale meteorology, medium and long-range forecasting, and upper-ocean dynamics.

Ghil, M.

Data Assimilation and Reanalysis

This presentation introduces the Data Assimilation and Reanalysis principle, then details the NASA Modern-Era Retrospective analysis for Research and Applications Version 2 (MERRA-2). The MERRA-2 is atmospheric reanalysis data spanning 1980 to the present. It has been produced by the NASA Global Modeling and Assimilation Office (GMAO) and is distributed by the NASA Goddard Earth Sciences Data and Information Services Center (GES DISC). In this presentation, I will introduce the MERRA-2 datasets associated with aerosol and air quality studies and use case studies to demonstrate the data tools developed at GES DISC to analyze and visualize MERRA-2 data, such as Giovanni and Jupyter Python notebook.

Xiaohua Pan

Observation impacts in the lower troposphere and challenges of Planetary Boundary Layer data assimilation

The Goddard Earth Observing System (GEOS) developed by the NASA Global Modeling and Assimilation Office assimilates a wide range of observations to support various NASA Earth Science missions. To set the stage for follow-on Planetary Boundary Layer (PBL) science and prepare for future observing systems of the next decade, we have assessed the effectiveness of the use of existing observing systems in the lower troposphere in GEOS, and analyzed model responses to the incremental analysis update (IAU) forcing. With a better understanding of the GEOS data assimilation algorithms in the PBL, we have developed strategies for improved PBL data assimilation in GEOS. The strategies to enhance data usages in both the data assimilation system and forecast model will be presented, and the utilization of PBL height data from multiple observing systems will be discussed as well.

Yanqiu Zhu

Advances in Data Assimilation and Weather Prediction Using TRMM Observations

Understanding the Earth's climate and how it responds to climate perturbations requires knowledge of how atmospheric moisture, clouds, latent heating, the large-scale circulation and energy fluxes vary with changing climatic conditions. The physical process linking these climate elements is precipitation. Accurate knowledge of how precipitation varies in space and time and how it couples with other atmospheric variables is essential for understanding the global water and energy cycle. In recent years, TRMM data products have played a key role in advancing the field of data assimilation to provide better global analyses for climate research and numerical weather prediction. TRMM research has demonstrated the effectiveness of microwave-based rainfall and total precipitable water (TPW) observations in improving the quality of assimilated datasets and upgrading forecast skills. TRMM latent heating products have also stimulated experimentation with innovative techniques to use this type of information to improve global analyses. We discuss strategies of assimilating TRMM observations at NASA s Data Assimilation Office and present results on the impact assimilating TRMM data on the Goddard Earth Observing System (GEOS) analyses and forecast capabilities.

Atlas, Robert

TRMM Data Assimilation at NASA

Current global analyses contain significant errors in primary hydrological fields such as precipitation, evaporation, and related cloud and moisture in the tropics. Work has been underway at NASA's Data Assimilation Office to explore the use of Tropical Rainfall Measuring Mission (TRMM) and Special Sensor Microwave Imager (SSM/I)-derived rainfall and total precipitable water (TPW) data in global data assimilation to directly constrain these hydrological parameters. We found that assimilating these data types improves not only the precipitation and moisture estimates but also key climate parameters directly linked to convection such as the outgoing longwave radiation, clouds, and the large-scale circulation in the tropics. We will present results showing that assimilating TRMM and SSM/I 6-hour averaged rain rates and TPW estimates significantly reduces the state-dependent systematic errors in assimilated products. Specifically, rainfall assimilation improves cloud and latent heating distributions, which, in turn, improves the cloudy-sky radiation and the large-scale circulation, while TPW assimilation reduces moisture biases to improve radiation in clear-sky regions. Rainfall and TPW assimilation also improves tropical forecasts beyond I day.

Hou, Arthur Y.

Fundamentals of Modeling, Data Assimilation, and High-performance Computing

This lecture will introduce the concepts of modeling, data assimilation and high- performance computing as it relates to the study of atmospheric composition. The lecture will work from basic definitions and will strive to provide a framework for thinking about development and application of models and data assimilation systems. It will not provide technical or algorithmic information, leaving that to textbooks, technical reports, and ultimately scientific journals. References to a number of textbooks and papers will be provided as a gateway to the literature.

Rood, Richard B.

Using Data Assimilation Methods of Prediction of Solar Activity

The variable solar magnetic activity known as the 11-year solar cycle has the longest history of solar observations. These cycles dramatically affect conditions in the heliosphere and the Earth's space environment. Our current understanding of the physical processes that make up global solar dynamics and the dynamo that generates the magnetic fields is sketchy, resulting in unrealistic descriptions in theoretical and numerical models of the solar cycles. The absence of long-term observations of solar interior dynamics and photospheric magnetic fields hinders development of accurate dynamo models and their calibration. In such situations, mathematical data assimilation methods provide an optimal approach for combining the available observational data and their uncertainties with theoretical models in order to estimate the state of the solar dynamo and predict future cycles. In this presentation, we will discuss the implementation and performance of an Ensemble Kalman Filter data assimilation method based on the Parker migratory dynamo model, complemented by the equation of magnetic helicity conservation and long-term sunspot data series. This approach has allowed us to reproduce the general properties of solar cycles and has already demonstrated a good predictive capability for the current cycle, 24. We will discuss further development of this approach, which includes a more sophisticated dynamo model, synoptic magnetogram data, and employs the DART Data Assimilation Research Testbed.

Prediction

Using Data Assimilation Methods for Prediction of Solar Activity

The variability of solar magnetic activity known as the 11-­‐year solar cycles has the longest history of observations. These solar cycles dramatically affect conditions in the heliosphere and the Earth’s space environment. Our current understanding of the physical processes that make up global solar dynamics and the dynamo that generates the magnetic fields is sketchy, resulting in unrealistic descriptions in theoretical and numerical models of the solar cycles. The absence of long-­‐term observations of solar interior dynamics and photospheric magnetic fields hinders development of accurate dynamo models and their calibration. In such situations, mathematical data assimilation methods provide an optimal approach for combining the available observational data and their uncertainties with theoretical models in order to estimate the state of the solar dynamo and predict future cycles. In this presentation, we will discuss the implementation and performance of an Ensemble Kalman Filter data assimilation method based on the Parker migratory dynamo model complemented by the equation of magnetic helicity conservation and long-­‐term sunspot data series. This approach has allowed us to reproduce the general properties of the solar cycles and has already demonstrated a good predictive capability for the current cycle, 24. We will discuss further development of this approach, which includes a more sophisticated dynamo model, synoptic magnetogram data, and employs the DART Data Assimilation Research Testbed.

Solar Activity

Data Assimilation Experiments Using Quality Controlled AIRS Version 5 Temperature Soundings

The AIRS Science Team Version 5 retrieval algorithm has been finalized and is now operational at the Goddard DAAC in the processing (and reprocessing) of all AIRS data. The AIRS Science Team Version 5 retrieval algorithm contains a number of significant improvements over Version 4. Two very significant improvements are described briefly below. 1) The AIRS Science Team Radiative Transfer Algorithm (RTA) has now been upgraded to accurately account for effects of non-local thermodynamic equilibrium on the AIRS observations. This allows for use of AIRS observations in the entire 4.3 micron CO2 absorption band in the retrieval algorithm during both day and night. Following theoretical considerations, tropospheric temperature profile information is obtained almost exclusively from clear column radiances in the 4.3 micron CO2 band in the AIRS Version 5 temperature profile retrieval step. These clear column radiances are a derived product that are indicative of radiances AIRS channels would have seen if the field of view were completely clear. Clear column radiances for all channels are determined using tropospheric sounding 15 micron CO2 observations. This approach allows for the generation of accurate values of clear column radiances and T(p) under most cloud conditions. 2) Another very significant improvement in Version 5 is the ability to generate accurate case-by-case, level-by-level error estimates for the atmospheric temperature profile, as well as for channel-by-channel clear column radiances. These error estimates are used for quality control of the retrieved products. Based on error estimate thresholds, each temperature profiles is assigned a characteristic pressure, pg, down to which the profile is characterized as good for use for data assimilation purposes. We have conducted forecast impact experiments assimilating AIRS quality controlled temperature profiles using the NASA GEOS-5 data assimilation system, consisting of the NCEP GSI analysis coupled with the NASA FVGCM, at a spatial resolution of 0.5 deg by 0.5 deg. Assimilation of Quality Controlled AIRS temperature profiles down to pg resulted in significantly improved forecast skill compared to that obtained from experiments when all data used operationally by NCEP, except for AIRS data, is assimilated. These forecasts were also significantly better than to those obtained when AIRS radiances (rather than temperature profiles) are assimilated, which is the way AIRS data is used operationally by NCEP and ECMWF.

Susskind, Joel

Randomized Preconditioned Solvers for Strong Constraint 4D-Var Data Assimilation

The Strong Constraint 4D Variational (SC-4DVAR) data assimilation method is widely used in climate and weather applications. SC-4DVAR involves solving a minimization problem to compute the maximum a posteriori estimate, which we tackle using the Gauss-Newton method. The computation of the descent direction is expensive since it involves the solution of a large-scale and potentially ill-conditioned linear system, solved using the preconditioned conjugate gradient (PCG) method. Here, to address this cost, we efficiently construct scalable preconditioners using three different randomization techniques, which all rely on a certain low-rank structure involving the Gauss-Newton Hessian. The proposed techniques come with theoretical guarantees on the condition number, and at the same time, are amenable to parallelization. We also develop an adaptive approach to estimate the sketch size and choose between the reuse or recomputation of the preconditioner. We demonstrate the performance and effectiveness of our methodology on two representative model problems—the Burgers and barotropic vorticity equation—showing a drastic reduction in both the number of PCG iterations and the number of Gauss-Newton Hessian products after including the preconditioner construction cost.

Gauss-Newton

Benchmarking a Soil Moisture Data Assimilation System for Agricultural Drought Monitoring

Despite considerable interest in the application of land surface data assimilation systems (LDAS) for agricultural drought applications, relatively little is known about the large-scale performance of such systems and, thus, the optimal methodological approach for implementing them. To address this need, this paper evaluates an LDAS for agricultural drought monitoring by benchmarking individual components of the system (i.e., a satellite soil moisture retrieval algorithm, a soil water balance model and a sequential data assimilation filter) against a series of linear models which perform the same function (i.e., have the same basic inputoutput structure) as the full system component. Benchmarking is based on the calculation of the lagged rank cross-correlation between the normalized difference vegetation index (NDVI) and soil moisture estimates acquired for various components of the system. Lagged soil moistureNDVI correlations obtained using individual LDAS components versus their linear analogs reveal the degree to which non-linearities andor complexities contained within each component actually contribute to the performance of the LDAS system as a whole. Here, a particular system based on surface soil moisture retrievals from the Land Parameter Retrieval Model (LPRM), a two-layer Palmer soil water balance model and an Ensemble Kalman filter (EnKF) is benchmarked. Results suggest significant room for improvement in each component of the system.

soil moisture