SEARCH · Search NASA

Results for “Error Covariance”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11

All sources NASA NTRS DOE OSTI

NASA NTRS2002

Evaluation of the DAO Retrospective Data Assimilation System

We have developed and implemented a retrospective data assimilation system (RDAS) as an upgrade to the operational DAO/Terra data assimilation system. This formulation aims at improving analysis over filter analysis by the dynamically consistent incorporation of observation information past a given analysis time. The current implementation of the RDAS uses the adjoint of the tangent linear model of a simplified version of the Terra general circulation model and extensions to the physical-space statistical analysis system to propagate observation information back in time. The RDAS adopts the same assumptions of the regular data assimilation system, particularly, no explicit propagation of error covariances are involved therefore rendering a procedure that is computationally affordable. In this study, we show results of experiments conducted to investigate the performance of the 6-hour (lag-1) RDAS. Statistical results obtained over one month during a winter season indicate that the RDAS represents considerable improvement over the regular assimilation. Plans for implementation of the RDAS capability in our new finite-volume data assimilation system will also be presented at the time of the conference.

Zhu, Yanqiu↗

NASA NTRS2003

Collaborative Research: Atlantic Ocean Tropical/Subtropical Processes from Seasonal to Decadal Time Scales: Model/Data, Model/Model Comparison and Model/Data Synthesis Through Assimilation

The effort of this first year of research has been focused on the assimilation of TOPEX/Poseidon altimetric data into a primitive equation model of the Atlantic tropical/subtropical circulation. A reduced-rank, stationary Kalman filter has been constructed to assimilate the altimetric sea surface height anomaly (SHA) into the model. The goal is to assess how the inter-hemispheric transports between the Atlantic subtropics and tropics are affected by the assimilation and how the subsurface thermocline structure , and its variability ,is dynamically constrained by the SHA. The model is a reduced-gravity primitive equation GCM of the upper Atlantic Ocean between 30 S and 30 N. The assimilation scheme is an approximation to the extended Kalman filter in which the error covariances of the state estimates are calculated only in a reduced- dimension subspace. The subspace is defined by the leading empirical orthogonal functions calculated from an unconstrained model calculation. Both an identical twin experiment using simulated SHA observations and assimilation of the real TOPEX data were performed. Results from the twin experiments demonstrate the ability of the method to constrain the ocean circulation and the subsurface temperature structure. The impact on the subsurface temperature structure of TOPEX assimilation was assessed using data from expandable bathythermographs. This showed a substantial improvement in the estimated temperature variability only within 13 degrees in latitude around the equator. The impact of TOPEX SHA assimilation on zonally integrated meridional transport across different latitudes was also estimated. Again within 13 degrees from the equator both the mean amplitude and interannual variability of the surface and subsurface transports were significantly enhanced, while the transports were insensitive to the assimilation in the subtropics.

Malannotte-Rizzoli, Paola↗

NASA NTRS2003

Analysis and Forecast of Two Storms Characterized by Extreme Deepening Rates

Between 25 and 27 December 1999 two very intense cyclones, named Lothar and Martin, swept across northern and western France causing substantial life and property loss. In this work, the finite volume general circulation model and data assimilation system (fvDAS) developed at the Data Assimilation Office of the NASA Goddard Space and Flight Center is being used to investigate these storms. In the first part of this article the dynamics of the storms is analyzed, and some important mechanisms are unveiled. The second part describes a set of eleven data assimilation experiments to study the impact of different data types on the automated analyses. Cloud-track winds provided by EUMETSAT and surface winds from QuikSCAT are being used. These data are assimilated with a range of different parameter settings of the forecast error covariance model. The results show that generally the additional wind data set have positive impacts on the analyses: particularly, the analysis of Lothar can be slightly improved by using the Eumetsat winds, and the analysis of Martin can be strongly improved by using the full-resolution QuikSCAT winds with a more localized influence. The third part of this article is focused on the forecast of Lothar which is very well predicted in the 1-5 day range by the fvDAS system.

Reale, Oreste↗

NASA NTRS2002

Sensitivity Studies for Assimilated Ozone Profiles

An ozone data assimilation system at the NASA/Goddard Data Assimilation Office (DAO) produces three-dimensional global ozone fields. They are obtained by assimilating ozone retrieved from the Solar Backscatter UltraViolet/2 (SBUV/2) instrument and the Earth Probe Total Ozone Mapping Spectrometer (EP TOMS) measurements into an off-line parameterized chemistry and transport model. In this talk we focus on the quality of lower stratospheric assimilated ozone profiles. Ozone in the lower stratosphere plays a key role in the forcing of climate. A biased ozone field in this region will adversely impact calculations of the stratosphere-troposphere exchange and, when used as a first guess in retrievals, the values determined from satellite observations. The SBUV/2 ozone data have a coarse vertical resolution with increased uncertainty below the ozone maximum, and TOMS provides only total ozone columns. Thus, the assimilated ozone profiles in the lower stratosphere are only weakly constrained by the incoming SBUV and TOMS data. Consequently, the assimilated ozone distribution should be sensitive to changes in inputs to the statistical analysis scheme. We investigate the sensitivity of assimilated ozone profiles to changes in a variety of system inputs: TOMS and SBUV/2 data selection, forecast and observations error covariance models, inclusion or omission of a parameterized chemistry model, and different versions of DAO assimilated wind fields used to drive the transport model. Comparisons of assimilated ozone fields with independent observations, primarily ozone sondes, are used to determine the impact of each of these changes.

Stajner, Ivanka↗

NASA NTRS2002

Monitoring and Assimilation of MLS Measurements in the DAO Ozone Data Assimilation System

Since 1999 the DAO (Data Assimilation Office) at NASA Goddard has operationally assimilated ozone measurements in near real time. Currently, the assimilation system analyzes SBUV profile and total column measurements using an off line CTM with parameterized chemistry within the 3D-PSAS algorithm. During the last year the assimilation system was modified to either monitor or actively assimilate MLS (Microwave Limb Sounder) measurements in conjunction with the active assimilation of TOMS total column and SBUV profiles. It is expected that the active assimilation of MLS profiles will improve analysis results in two ways. First, there should be an improvement in the vertical resolution. Second, there should be an improvement in regions where SBUV measurements do not exist (such as in the polar night). A series of experiments using UARS (Upper Atmospheric Research Satellite) MLS data from December 1991 to March 1992 were run. In these experiments MLS data was either monitored in conjunction with active assimilation of the SBUV profile and TOMS total column measurements, or some combination of MLS, TOMS, and SBUV observations was actively assimilated. Monitoring of MLS data indicated that the MLS observations contain verifiable information that should improve the vertical structure of the analysis results, especially below the ozone peak and above the tropopause. The monitoring also substantiated the potential to improve the assimilation in the polar night. Active assimilation of MLS data does indeed improve the analysis results in these two ways, although the quality of the improvements is not uniform. This suggests that refinement of the specification of the error covariances might be needed to optimize the system. In addition it may be necessary to account for biases between the different sources of ozone information.

Atlas, Robert↗

NASA NTRS2003

Assimilation of Satellite Ozone Observations

This talk will discuss assimilation of ozone data from satellite-borne instruments. Satellite observations of ozone total columns and profiles have been measured by a series of Total Ozone Mapping Spectrometer (TOMS), Solar Backscatter Ultraviolet (SBUV) instruments, and more recently by the Global Ozone Monitoring Experiment. Additional profile data are provided by instruments on NASA's Upper Atmosphere Research Satellite and by occultation instruments on other platforms. Instruments on Envisat' and future EOS Aura satellite will supply even more comprehensive data about the ozone distribution. Satellite data contain a wealth of information, but they do not provide synoptic global maps of ozone fields. These maps can be obtained through assimilation of satellite data into global chemistry and transport models. In the ozone system at NASA's Data Assimilation Office (DAO) any combination of TOMS, SBUV, and Microwave Limb sounder (MLS) data can be assimilated. We found that the addition of MLS to SBUV and TOMS data in the system helps to constrain the ozone distribution, especially in the polar night region and in the tropics. The assimilated ozone distribution in the troposphere and lower stratosphere is sensitive also to finer changes in the SBUV and TOMS data selection and to changes in error covariance models. All results are established by comparisons of assimilated ozone with independent profiles from ozone sondes and occultation instruments.

Stajner, I.↗

NASA NTRS2003

Sensitivity of the Assimilated Ozone in the UTLS to Model and Data Selection Changes

This presentation will discuss the sensitivity of assimilated ozone fields in the upper troposphere and lower stratosphere (UTLS) to a number of factors, focusing mainly on aspects of data selection and the prediction model. This is important, because assimilation represents an attempt to construct our best estimates of the true ozone field; however, inaccuracies in the UTLS ozone distribution translate into an uncertainty in factors such as the calculated radiative forcing of climate or the inferred stratosphere-troposphere exchange (STE) of ozone. The 3D ozone data assimilation system, from NASA's Global Modeling and Assimilation Office (GMAO), combines observations of total ozone column and stratospheric profiles with predictions from an off-line, parameterized chemistry and transport model (pCTM) to produce six-hourly, global analyses. The first experiments discussed assimilate ozone retrievals from the Earth-Probe Total Ozone Mapping Spectrometer (EPTOMS) and stratospheric profiles from the Solar Backscatter UltraViolet/2 (SBUV/2) instrument. The SBUV/2 ozone data have a coarse vertical resolution, with increased uncertainty below the ozone maximum, and TOMS provides only total ozone columns. Thus, the assimilated ozone profiles in the UTLS region are only weakly constrained by the incoming SBUV and TOMS data. Consequently, the assimilated ozone distribution should be sensitive to changes in inputs to the statistical analysis scheme. Sensitivity studies have been conducted to examine the responses to TOMS and SBUV/2 data selection, modifications of the forecast and observation error covariance models, and the model formulation (turning off chemistry or using different wind analyses in the pCTM). The second set of experiments includes an additional data type: ozone retrieved from infrared limb-emission by MIPAS on Envisat. These data offer not only improved vertical resolution in the stratosphere, but also give measurements in the polar night. Comparisons of the assimilated ozone fields from both sets of experiments with independent observations, primarily ozone sondes, are used to determine the impact of each of these changes. It is shown that many of the changes have a significant impact on the UTLS ozone estimates. Implications for interpretation of STE and radiative forcing of climate are discussed.

Pawson, Steven↗

NASA NTRS2004

Can Assimilation of Satellite Ozone Data Contribute to the Understanding of the Lower Stratospheric Ozone?

We study the quality of lower stratospheric ozone fields from a three- dimensional global ozone assimilation system. Ozone in this region is important for the forcing of climate, but its global distribution is not fully known because of its large temporal and vertical variability. Modeled fields often have biases due to the inaccurate representation of transport processes in this region with strong gradients. Accurate ozonesonde or satellite occultation measurements have very limited coverage. Nadir measurements, such as those from the Solar Backscatter Ultraviolet/2 (SBUV/2) instrument that provide wide latitudinal coverage, lack the vertical resolution needed to represent sharp vertical features. Limb measurements, such as those from the Michelson Interferometer for Passive Atmospheric Sounding (MIPAS), provide a finer vertical resolution. We show that assimilation of MIPAS data in addition to SBUV/2 data leads to better estimates of ozone in comparison with independent high quality satellite, aircraft, and ozone sonde measurements. Other modifications to the statistical analysis that have an impact on the lower stratospheric ozone will be mentioned: error covariance modeling and data selection. Direct and indirect impacts of transport and chemistry models will be discussed. Implications for multi-year analyses and short-tern prediction will be addressed.

Stajner, I.↗

NASA NTRS2004

ENSO Bred Vectors in Coupled Ocean-Atmosphere General Circulation Models

The breeding method has been implemented in the NASA Seasonal-to-Interannual Prediction Project (NSIPP) Coupled General Circulation Model (CGCM) with the goal of improving operational seasonal to interannual climate predictions through ensemble forecasting and data assimilation. The coupled instability as cap'tured by the breeding method is the first attempt to isolate the evolving ENSO instability and its corresponding global atmospheric response in a fully coupled ocean-atmosphere GCM. Our results show that the growth rate of the coupled bred vectors (BV) peaks at about 3 months before a background ENSO event. The dominant growing BV modes are reminiscent of the background ENSO anomalies and show a strong tropical response with wind/SST/thermocline interrelated in a manner similar to the background ENSO mode. They exhibit larger amplitudes in the eastern tropical Pacific, reflecting the natural dynamical sensitivity associated with the presence of the shallow thermocline. Moreover, the extratropical perturbations associated with these coupled BV modes reveal the variations related to the atmospheric teleconnection patterns associated with background ENSO variability, e.g. over the North Pacific and North America. A similar experiment was carried out with the NCEP/CFS03 CGCM. Comparisons between bred vectors from the NSIPP CGCM and NCEP/CFS03 CGCM demonstrate the robustness of the results. Our results strongly suggest that the breeding method can serve as a natural filter to identify the slowly varying, coupled instabilities in a coupled GCM, which can be used to construct ensemble perturbations for ensemble forecasts and to estimate the coupled background error covariance for coupled data assimilation.

Yang, S. C.↗

NASA NTRS2003

An optimal modification of a Kalman filter for time scales

The Kalman filter in question, which was implemented in the time scale algorithm TA(NIST), produces time scales with poor short-term stability. A simple modification of the error covariance matrix allows the filter to produce time scales with good stability at all averaging times, as verified by simulations of clock ensembles.

time scale Kalman filter↗

NASA NTRS2007

Non-Parametric Collision Probability for Low-Velocity Encounters

An implicit, but not necessarily obvious, assumption in all of the current techniques for assessing satellite collision probability is that the relative position uncertainty is perfectly correlated in time. If there is any mis-modeling of the dynamics in the propagation of the relative position error covariance matrix, time-wise de-correlation of the uncertainty will increase the probability of collision over a given time interval. The paper gives some examples that illustrate this point. This paper argues that, for the present, Monte Carlo analysis is the best available tool for handling low-velocity encounters, and suggests some techniques for addressing the issues just described. One proposal is for the use of a non-parametric technique that is widely used in actuarial and medical studies. The other suggestion is that accurate process noise models be used in the Monte Carlo trials to which the non-parametric estimate is applied. A further contribution of this paper is a description of how the time-wise decorrelation of uncertainty increases the probability of collision.

Carpenter, J. Russell↗

NASA NTRS2001

Stable and verifiable state estimation methods and systems with spacecraft applications

The stability of a recursive estimator process (e.g., a Kalman filter is assured for long time periods by periodically resetting an error covariance P(t.sub.n) of the system to a predetermined reset value P.sub.r. The recursive process is thus repetitively forced to start from a selected covariance and continue for a time period that is short compared to the system's total operational time period. The time period in which the process must maintain its numerical stability is significantly reduced as is the demand on the system's numerical stability. The process stability for an extended operational time period T.sub.o is verified by performing the resetting step at the end of at least one reset time period T.sub.r whose duration is less than the operational time period T.sub.o and then confirming stability of the process over the reset time period T.sub.r. Because the recursive process starts from a selected covariance at the beginning of each reset time period T.sub.r, confirming stability of the process over at least one reset time period substantially confirms stability over the longer operational time period T.sub.o.

Li, Rongsheng↗

NASA NTRS2012

Multi-Scale Three-Dimensional Variational Data Assimilation System for Coastal Ocean Prediction

A multi-scale three-dimensional variational data assimilation system (MS-3DVAR) has been formulated and the associated software system has been developed for improving high-resolution coastal ocean prediction. This system helps improve coastal ocean prediction skill, and has been used in support of operational coastal ocean forecasting systems and field experiments. The system has been developed to improve the capability of data assimilation for assimilating, simultaneously and effectively, sparse vertical profiles and high-resolution remote sensing surface measurements into coastal ocean models, as well as constraining model biases. In this system, the cost function is decomposed into two separate units for the large- and small-scale components, respectively. As such, data assimilation is implemented sequentially from large to small scales, the background error covariance is constructed to be scale-dependent, and a scale-dependent dynamic balance is incorporated. This scheme then allows effective constraining large scales and model bias through assimilating sparse vertical profiles, and small scales through assimilating high-resolution surface measurements. This MS-3DVAR enhances the capability of the traditional 3DVAR for assimilating highly heterogeneously distributed observations, such as along-track satellite altimetry data, and particularly maximizing the extraction of information from limited numbers of vertical profile observations.

Li, Zhijin↗

NASA NTRS2012

Constraining Middle Atmospheric Moisture in GEOS-5 Using EOS-MLS Observations

Middle atmospheric water vapor plays an important role in climate and atmospheric chemistry. In the middle atmosphere, water vapor, after ozone and carbon dioxide, is an important radiatively active gas that impacts climate forcing and the energy balance. It is also the source of the hydroxyl radical (OH) whose abundances affect ozone and other constituents. The abundance of water vapor in the middle atmosphere is determined by upward transport of dehydrated air through the tropical tropopause layer, by the middle atmospheric circulation, production by the photolysis of methane (CH4), and other physical and chemical processes in the stratosphere and mesosphere. The Modern-Era Retrospective analysis for Research and Applications (MERRA) reanalysis with GEOS-5 did not assimilate any moisture observations in the middle atmosphere. The plan is to use such observations, available sporadically from research satellites, in future GEOS-5 reanalyses. An overview will be provided of the progress to date with assimilating the EOS-Aura Microwave Limb Sounder (MLS) moisture retrievals, alongside ozone and temperature, into GEOS-5. Initial results demonstrate that the MLS observations can significantly improve the middle atmospheric moisture field in GEOS-5, although this result depends on introducing a physically meaningful representation of background error covariances for middle atmospheric moisture into the system. High-resolution features in the new moisture field will be examined, and their relationships with ozone, in a two-year assimilation experiment with GEOS-5. Discussion will focus on how Aura MLS moisture observations benefit the analyses.

Jin, Jianjun↗

NASA NTRS2012

Quantifying Interannual Variability of the UTLS Ozone Using Assimilation of Satellite Data

An accurate representation of spatial and temporal variability of the Upper Troposphere Lower Stratosphere (UTLS) ozone is essential for understanding both the tropospheric ozone budget and ozone s contribution to radiative forcing. The complex, dynamically driven structure of trace gas fields in the UTLS presents a challenge to data-based and modelling studies. Small features are not fully resolved in data from limb-sounding instruments such as the Microwave Limb Sounder on EOS-Aura (the EOS-MLS), but are captured in assimilation of those data as vertical structure is added from the assimilated meteorology. This will be demonstrated using a multi-year assimilation of EOS-MLS observations in the Goddard Earth Observing System, Version 5 (GEOS-5) data assimilation system. The results demonstrate the realism of the seasonal and year to year variability of laminar structures in the mid-latitudinal ozone field between years 2005-2007, for which independent validation data are available from the HIRDLS instrument. The analysis is done in the context of the underlying large scale dynamics. The lifetimes of most research instruments are too short for them to be used throughout the duration of long-term (at least 3 decades) reanalyses. For example, the EOS-MLS instrument has operated since mid-2004 until present. By contrast, Solar Backscatter Ultra Violet (SBUV) measurements provide continuous data since late 1978, but their vertical resolution is insufficient to represent the profile shape in the UTLS. Assimilation of these SBUV/2 observations in the GEOS-5 data assimilation system has hitherto not captured a realistic ozone structure in the UTLS, even though transport studies using GEOS-5 wind fields do show such structures. We show that careful construction of the background error covariance structure in GEOS-5 can lead to more realistic UTLS ozone fields when assimilating SBUV/2 observations. The reasoning behind this will be discussed, emphasizing the need to retain the sharp gradient of ozone concentrations across the tropopause. We analyze the UTLS ozone distributions in multi-year SBUV/2 assimilation experiments, comparing the results against the independent HIRDLSdataset and, for a longer period, with the MLS assimilation and discuss the consequences for tropospheric ozone and radiative forcing.

Wargan, K.↗

NASA NTRS2011

Reduced Kalman Filters for Clock Ensembles

This paper summarizes the author's work ontimescales based on Kalman filters that act upon the clock comparisons. The natural Kalman timescale algorithm tends to optimize long-term timescale stability at the expense of short-term stability. By subjecting each post-measurement error covariance matrix to a non-transparent reduction operation, one obtains corrected clocks with improved short-term stability and little sacrifice of long-term stability.

Kalman filter↗

NASA NTRS2015

Assimilation of NUCAPS Retrieved Profiles in GSI for Unique Forecasting Applications

Hyperspectral IR profiles can be assimilated in GSI as a separate observation other than radiosondes with only changes to tables in the fix directory. Assimilation of profiles does produce changes to analysis fields and evidenced by: Innovations larger than +/-2.0 K are present and represent where individual profiles impact the final temperature analysis.The updated temperature analysis is colder behind the cold front and warmer in the warm sector. The updated moisture analysis is modified more in the low levels and tends to be drier than the original model background Analysis of model output shows: Differences relative to 13-km RAP analyses are smaller when profiles are assimilated with NUCAPS errors. CAPE is under-forecasted when assimilating NUCAPS profiles, which could be problematic for severe weather forecasting Refining the assimilation technique to incorporate an error covariance matrix and creating a separate GSI module to assimilate satellite profiles may improve results.

data assimilation↗

NASA NTRS2014

Carbon Monitoring System Flux Estimation and Attribution: Impact of ACOS-GOSAT X(CO2) Sampling on the Inference of Terrestrial Biospheric Sources and Sinks

Using an Observing System Simulation Experiment (OSSE), we investigate the impact of JAXA Greenhouse gases Observing SATellite 'IBUKI' (GOSAT) sampling on the estimation of terrestrial biospheric flux with the NASA Carbon Monitoring System Flux (CMS-Flux) estimation and attribution strategy. The simulated observations in the OSSE use the actual column carbon dioxide (X(CO2)) b2.9 retrieval sensitivity and quality control for the year 2010 processed through the Atmospheric CO2 Observations from Space algorithm. CMS-Flux is a variational inversion system that uses the GEOS-Chem forward and adjoint model forced by a suite of observationally constrained fluxes from ocean, land and anthropogenic models. We investigate the impact of GOSAT sampling on flux estimation in two aspects: 1) random error uncertainty reduction and 2) the global and regional bias in posterior flux resulted from the spatiotemporally biased GOSAT sampling. Based on Monte Carlo calculations, we find that global average flux uncertainty reduction ranges from 25% in September to 60% in July. When aggregated to the 11 land regions designated by the phase 3 of the Atmospheric Tracer Transport Model Intercomparison Project, the annual mean uncertainty reduction ranges from 10% over North American boreal to 38% over South American temperate, which is driven by observational coverage and the magnitude of prior flux uncertainty. The uncertainty reduction over the South American tropical region is 30%, even with sparse observation coverage. We show that this reduction results from the large prior flux uncertainty and the impact of non-local observations. Given the assumed prior error statistics, the degree of freedom for signal is approx.1132 for 1-yr of the 74 055 GOSAT X(CO2) observations, which indicates that GOSAT provides approx.1132 independent pieces of information about surface fluxes. We quantify the impact of GOSAT's spatiotemporally sampling on the posterior flux, and find that a 0.7 gigatons of carbon bias in the global annual posterior flux resulted from the seasonally and diurnally biased sampling when using a diagonal prior flux error covariance.

biased sampling↗