Engineering topics

Rienecker, Michele M.

Publications and source records attributed to Rienecker, Michele M..

At least 19 records

NASA NTRS2014

Background Error Covariance Estimation using Information from a Single Model Trajectory with Application to Ocean Data Assimilation into the GEOS-5 Coupled Model

An attractive property of ensemble data assimilation methods is that they provide flow dependent background error covariance estimates which can be used to update fields of observed variables as well as fields of unobserved model variables. Two methods to estimate background error covariances are introduced which share the above property with ensemble data assimilation methods but do not involve the integration of multiple model trajectories. Instead, all the necessary covariance information is obtained from a single model integration. The Space Adaptive Forecast error Estimation (SAFE) algorithm estimates error covariances from the spatial distribution of model variables within a single state vector. The Flow Adaptive error Statistics from a Time series (FAST) method constructs an ensemble sampled from a moving window along a model trajectory. SAFE and FAST are applied to the assimilation of Argo temperature profiles into version 4.1 of the Modular Ocean Model (MOM4.1) coupled to the GEOS-5 atmospheric model and to the CICE sea ice model. The results are validated against unassimilated Argo salinity data. They show that SAFE and FAST are competitive with the ensemble optimal interpolation (EnOI) used by the Global Modeling and Assimilation Office (GMAO) to produce its ocean analysis. Because of their reduced cost, SAFE and FAST hold promise for high-resolution data assimilation applications.

Data Assimilation↗

NASA NTRS2013

Decadal Prediction Skill in the GEOS-5 Forecast System

A suite of decadal predictions has been conducted with the NASA Global Modeling and Assimilation Office's (GMAO's) GEOS-5 Atmosphere-Ocean general circulation model. The hind casts are initialized every December 1st from 1959 to 2010, following the CMIP5 experimental protocol for decadal predictions. The initial conditions are from a multivariate ensemble optimal interpolation ocean and sea-ice reanalysis, and from GMAO's atmospheric reanalysis, the modern-era retrospective analysis for research and applications. The mean forecast skill of a three-member-ensemble is compared to that of an experiment without initialization but also forced with observed greenhouse gases. The results show that initialization increases the forecast skill of North Atlantic sea surface temperature compared to the uninitialized runs, with the increase in skill maintained for almost a decade over the subtropical and mid-latitude Atlantic. On the other hand, the initialization reduces the skill in predicting the warming trend over some regions outside the Atlantic. The annual-mean Atlantic meridional overturning circulation index, which is defined here as the maximum of the zonally-integrated overturning stream function at mid-latitude, is predictable up to a 4-year lead time, consistent with the predictable signal in upper ocean heat content over the North Atlantic. While the 6- to 9-year forecast skill measured by mean squared skill score shows 50 percent improvement in the upper ocean heat content over the subtropical and mid-latitude Atlantic, prediction skill is relatively low in the sub-polar gyre. This low skill is due in part to features in the spatial pattern of the dominant simulated decadal mode in upper ocean heat content over this region that differ from observations. An analysis of the large-scale temperature budget shows that this is the result of a model bias, implying that realistic simulation of the climatological fields is crucial for skillful decadal forecasts.

Decadal Prediction↗

NASA NTRS2013

An Assessment of the Skill of GEOS-5 Seasonal Forecasts

The seasonal forecast skill of the NASA Global Modeling and Assimilation Office coupled global climate model (CGCM) is evaluated based on an ensemble of 9-month lead forecasts for the period 1993 to 2010. The results from the current version (V2) of the CGCM consisting of the GEOS-5 AGM coupled to the MOM4 ocean model are compared with those from an earlier version (V1) in which the AGCM (the NSIPP model) was coupled to the Poseidon Ocean Model. It was found that the correlation skill of the Sea Surface Temperature (SST) forecasts is generally better in V2, especially over the sub-tropical and tropical central and eastern Pacific, Atlantic, and Indian Ocean. Furthermore, the improvement in skill in V2 mainly comes from better forecasts of the developing phase of ENSO from boreal spring to summer. The skill of ENSO forecasts initiated during the boreal winter season, however, shows no improvement in terms of correlation skill, and is in fact slightly worse in terms of root mean square error (RMSE). The degradation of skill is found to be due to an excessive ENSO amplitude. For V1, the ENSO amplitude is too strong in forecasts starting in boreal spring and summer, which causes large RMSE in the forecast. For V2, the ENSO amplitude is slightly stronger than that in observations and V1 for forecasts starting in boreal winter season. An analysis of the terms in the SST tendency equation, shows that this is mainly due to an excessive zonal advective feedback. In addition, V2 forecasts that are initiated during boreal winter season, exhibit a slower phase transition of El Nino, which is consistent with larger amplitude of ENSO after the ENSO peak season. It is found that this is due to weak discharge of equatorial Warm Water Volume (WWV). In both observations and V1, the discharge of equatorial WWV leads the equatorial geostrophic easterly current so as to damp the El Nino starting in January. This process is delayed by about 2 months in V2 due to the slower phase transition of the equatorial zonal current from westerly to easterly.

GMAO↗

NASA NTRS2013

Ensemble Data Assimilation Without Ensembles: Methodology and Application to Ocean Data Assimilation

Two methods to estimate background error covariances for data assimilation are introduced. While both share properties with the ensemble Kalman filter (EnKF), they differ from it in that they do not require the integration of multiple model trajectories. Instead, all the necessary covariance information is obtained from a single model integration. The first method is referred-to as SAFE (Space Adaptive Forecast error Estimation) because it estimates error covariances from the spatial distribution of model variables within a single state vector. It can thus be thought of as sampling an ensemble in space. The second method, named FAST (Flow Adaptive error Statistics from a Time series), constructs an ensemble sampled from a moving window along a model trajectory. The underlying assumption in these methods is that forecast errors in data assimilation are primarily phase errors in space and/or time.

Data Assimilation↗

NASA NTRS2012

The GEOS-iODAS: Description and Evaluation

This report documents the GMAO's Goddard Earth Observing System sea ice and ocean data assimilation systems (GEOS iODAS) and their evolution from the first reanalysis test, through the implementation that was used to initialize the GMAO decadal forecasts, and to the current system that is used to initialize the GMAO seasonal forecasts. The iODAS assimilates a wide range of observations into the ocean and sea ice components: in-situ temperature and salinity profiles, sea level anomalies from satellite altimetry, analyzed SST, and sea-ice concentration. The climatological sea surface salinity is used to constrain the surface salinity prior to the Argo years. Climatological temperature and salinity gridded data sets from the 2009 version of the World Ocean Atlas (WOA09) are used to help constrain the analysis in data sparse areas. The latest analysis, GEOS ODAS5.2, is diagnosed through detailed studies of the statistics of the innovations and analysis departures, comparisons with independent data, and integrated values such as volume transport. Finally, the climatologies of temperature and salinity fields from the Argo era, 2002-2011, are presented and compared with the WOA09.

Climate↗

NASA NTRS2012

Global Modeling and Assimilation Office Annual Report and Research Highlights 2011-2012

Over the last year, the Global Modeling and Assimilation Office (GMAO) has continued to advance our GEOS-5-based systems, updating products for both weather and climate applications. We contributed hindcasts and forecasts to the National Multi-Model Ensemble (NMME) of seasonal forecasts and the suite of decadal predictions to the Coupled Model Intercomparison Project (CMIP5).

Rienecker, Michele M.↗

NASA NTRS2012

Atmospheric Reanalyses-Recent Progress and Prospects for the Future. A Report from a Technical Workshop, April 2010

In April 2010, developers representing each of the major reanalysis centers met at Goddard Space Flight Center to discuss technical issues - system advances and lessons learned - associated with recent and ongoing atmospheric reanalyses and plans for the future. The meeting included overviews of each center s development efforts, a discussion of the issues in observations, models and data assimilation, and, finally, identification of priorities for future directions and potential areas of collaboration. This report summarizes the deliberations and recommendations from the meeting as well as some advances since the workshop.

Rienecker, Michele M.↗

NASA NTRS2010

Decadal Prediction Efforts in GMAO (Global Modeling and Assimilation Office)

The Global Modeling and Assimilation Office (GMAO) plans to use our GEOS-5 atmosphere-ocean general circulation model (AOGCM) to explore issues associated with predictability on decadal time scales and to contribute to the decadal prediction project that is part ofCMIP5. The GEOS-5 AOGCM is comprised of the GEOS-5 AGCM with the Catchment Land Surface Model, coupled to GFDL's MOM, version 4. We have assimilation systems for both the atmosphere and ocean. For our climate prediction efforts, the atmosphere will be initialized from the GEOS-5 Modem Era Retrospective-analysis for Research and Applications (MERRA), available from 1979 to present at 112 resolution, and from 1948 to present at 2 resolution. The ocean assimilation is conducted within the coupled model framework, using the MERRA as a constraint for both the atmosphere and the ocean. The decadal prediction experiments will be conducted with a 1 atmosphere and a 112 ocean. Some initial results will be presented, focusing on initialization aspects of the GEOS-5 system.

Rienecker, Michele M.↗

NASA NTRS2009

Synthesis and Assimilation Systems - Essential Adjuncts to the Observing System

Assimilation systems synthesize diverse in-situ and satellite data streams into full four-dimensional state estimates by combining the strengths of each data set and also of the model. The resulting analysis provides an integrated view of the information in the various observations as well as derived estimates of unobserved quantities. Assimilation systems are particularly important for the ocean where subsurface observations, even today, are sparse and intermittent compared with the scales needed to represent ocean variability and where satellites only sense the surface. Increasingly, models and assimilation systems are being used to provide information about the current observing system and to help in the design plans for new observations. Whether it is as a user of observations or a contributor to evaluation of the observing system, ocean synthesis and assimilation systems are now an integral part of the global ocean observing and information system. Major advances have been made over the last decade under the auspices of WCRP's Climate Variability and Predictability Project (CLIVAR) and the Global Ocean Data Assimilation Experiment (GODAE). In addition to advances in the assimilation systems, there have been major developments in the observing system, with satellite altimetry, the tropical moored buoy arrays in the Pacific and Atlantic, and more recently Argo. These developments have led to significant advances in our understanding and prediction of ocean variations at both mesoscale and climate scales. Many challenges remain. Some of these challenges lie in the observations themselves, some in the assimilation systems that, even in the more recent era of unprecedented observations from satellite altimetry and Argo, provide different views of climate variations. Yet there are many examples of successful applications from ocean assimilation products. Use of these systems for assessing the observing system helps identify the strengths of each observation type, and indicates that none of the current observations is redundant. Indeed, the indication is that the ocean remains under-sampled and that further improvements in the observing system are needed for both climate monitoring and prediction. Future developments will be increasingly towards consistent analyses across components of the Earth system using, e.g., coupled atmosphere-ocean models.

Rienecker, Michele M.↗

NASA NTRS2009

Synthesis and Assimilation Systems - Essential Adjuncts to the Global Ocean Observing System

Ocean assimilation systems synthesize diverse in situ and satellite data streams into four-dimensional state estimates by combining the various observations with the model. Assimilation is particularly important for the ocean where subsurface observations, even today, are sparse and intermittent compared with the scales needed to represent ocean variability and where satellites only sense the surface. Developments in assimilation and in the observing system have advanced our understanding and prediction of ocean variations at mesoscale and climate scales. Use of these systems for assessing the observing system helps identify the strengths of each observation type. Results indicate that the ocean remains under-sampled and that further improvements in the observing system are needed. Prospects for future advances lie in improved models and better estimates of error statistics for both models and observations. Future developments will be increasingly towards consistent analyses across components of the Earth system. However, even today ocean synthesis and assimilation systems are providing products that are useful for many applications and should be considered an integral part of the global ocean observing and information system.

Rienecker, Michele M.↗

NASA NTRS2004

Multivariate Error Covariance Estimates by Monte-Carlo Simulation for Assimilation Studies in the Pacific Ocean

One of the most difficult aspects of ocean state estimation is the prescription of the model forecast error covariances. The paucity of ocean observations limits our ability to estimate the covariance structures from model-observation differences. In most practical applications, simple covariances are usually prescribed. Rarely are cross-covariances between different model variables used. Here a comparison is made between a univariate Optimal Interpolation (UOI) scheme and a multivariate OI algorithm (MvOI) in the assimilation of ocean temperature. In the UOI case only temperature is updated using a Gaussian covariance function and in the MvOI salinity, zonal and meridional velocities as well as temperature, are updated using an empirically estimated multivariate covariance matrix. Earlier studies have shown that a univariate OI has a detrimental effect on the salinity and velocity fields of the model. Apparently, in a sequential framework it is important to analyze temperature and salinity together. For the MvOI an estimation of the model error statistics is made by Monte-Carlo techniques from an ensemble of model integrations. An important advantage of using an ensemble of ocean states is that it provides a natural way to estimate cross-covariances between the fields of different physical variables constituting the model state vector, at the same time incorporating the model's dynamical and thermodynamical constraints as well as the effects of physical boundaries. Only temperature observations from the Tropical Atmosphere-Ocean array have been assimilated in this study. In order to investigate the efficacy of the multivariate scheme two data assimilation experiments are validated with a large independent set of recently published subsurface observations of salinity, zonal velocity and temperature. For reference, a third control run with no data assimilation is used to check how the data assimilation affects systematic model errors. While the performance of the UOI and MvOI is similar with respect to the temperature field, the salinity and velocity fields are greatly improved when multivariate correction is used, as evident from the analyses of the rms differences of these fields and independent observations. The MvOI assimilation is found to improve upon the control run in generating the water masses with properties close to the observed, while the UOI failed to maintain the temperature and salinity structure.

Borovikov, Anna↗

NASA NTRS2003

Temperature Data Assimilation with Salinity Corrections: Validation for the NSIPP Ocean Data Assimilation System in the Tropical Pacific Ocean, 1993-1998

The NASA Seasonal-to-Interannual Prediction Project (NSIPP) has developed an Ocean data assimilation system to initialize the quasi-isopycnal ocean model used in our experimental coupled-model forecast system. Initial tests of the system have focused on the assimilation of temperature profiles in an optimal interpolation framework. It is now recognized that correction of temperature only often introduces spurious water masses. The resulting density distribution can be statically unstable and also have a detrimental impact on the velocity distribution. Several simple schemes have been developed to try to correct these deficiencies. Here the salinity field is corrected by using a scheme which assumes that the temperature-salinity relationship of the model background is preserved during the assimilation. The scheme was first introduced for a zlevel model by Troccoli and Haines (1999). A large set of subsurface observations of salinity and temperature is used to cross-validate two data assimilation experiments run for the 6-year period 1993-1998. In these two experiments only subsurface temperature observations are used, but in one case the salinity field is also updated whenever temperature observations are available.

Troccoli, Alberto↗

NASA NTRS2001

Design and Implementation of a Parallel Multivariate Ensemble Kalman Filter for the Poseidon Ocean General Circulation Model

A multivariate ensemble Kalman filter (MvEnKF) implemented on a massively parallel computer architecture has been implemented for the Poseidon ocean circulation model and tested with a Pacific Basin model configuration. There are about two million prognostic state-vector variables. Parallelism for the data assimilation step is achieved by regionalization of the background-error covariances that are calculated from the phase-space distribution of the ensemble. Each processing element (PE) collects elements of a matrix measurement functional from nearby PEs. To avoid the introduction of spurious long-range covariances associated with finite ensemble sizes, the background-error covariances are given compact support by means of a Hadamard (element by element) product with a three-dimensional canonical correlation function. The methodology and the MvEnKF configuration are discussed. It is shown that the regionalization of the background covariances; has a negligible impact on the quality of the analyses. The parallel algorithm is very efficient for large numbers of observations but does not scale well beyond 100 PEs at the current model resolution. On a platform with distributed memory, memory rather than speed is the limiting factor.

Keppenne, Christian L.↗

NASA NTRS2001

The Impacts of Daily Surface Forcing in the Upper Ocean over Tropical Pacific: A Numerical Study

Tropical Pacific Ocean is an important region that affects global climate. How the ocean responds to the atmospheric surface forcing (surface radiative, heat and momentum fluxes) is a major topic in oceanographic research community. The ocean becomes warm when more heat flux puts into the ocean. The monthly mean forcing has been used in the past years since daily forcing was unavailable due to the lack of observations. The daily forcing is now available from the satellite measurements. This study investigates the response of the upper ocean over tropical Pacific to the daily atmospheric surface forcing. The ocean surface heat budgets are calculated to determine the important processes for the oceanic response. The differences of oceanic responses between the eastern and western Pacific are intensively discussed.

Sui, C.-H.↗

NASA NTRS2001

A Modeling Study of Oceanic Response to Daily and Monthly Surface Forcing

The goal of this study is to investigate the effect of high-frequency surface forcing (wind stresses and heat fluxes) on upper-ocean response. We use the reduced-gravity quasi-isopycnal ocean model by Schopf and Loughe (1995) for this study. Two experiments are performed: one with daily and the other with monthly surface forcing. The two experiments are referred to as DD and MM, respectively. The daily surface wind stress is produced from the SSM/I wind data (Atlas et al. 1991) using the drag coefficient of Large and Pond (1982). The surface latent and sensible heat fluxes are estimated using the atmospheric mixed layer model by Seager et al. (1995) with the time-varying air temperature and specific humidity from the NCEP-NCAR reanalysis (Kalnay et al. 1996). The radiation is based on climatological shortwave radiation from the Earth Radiation Budget Experiment (ERBE) [Harrison et al. 1993] and the daily GEWEX SRB data. The ocean model domain is restricted to the Pacific Ocean with realistic land boundaries. At the southern boundary the model temperature and salinity are relaxed to the Levitus (1994) climatology. The time-mean SST distribution from MM is close to the observed SST climatology while the mean SST field from DD is about 1.5 C cooler. To identify the responsible processes, we examined the mean heat budgets and the heat balance during the first year (when the difference developed) in the two experiments. The analysis reveals that this is contributed by two factors. One is the difference in latent heat flux. The other is the difference in mixing processes. To further evaluate the responsible processes, we repeated the DD experiment by reducing the based vertical diffusion from 1e-4 to 0.5e-5. The resultant SST field becomes quite closer to the observed SST field. SST variability from the two experiments is generally similar, but the equatorial SST differences between the two experiments show interannual variations. We are investigating the possible mechanisms responsible for the different responses.

Sui, Chung-Hsiung↗

NASA NTRS1999

Use of Ocean Remote Sensing Data to Enhance Predictions with a Coupled General Circulation Model

Surface height, sea surface temperature and surface wind observations from satellites have given a detailed time sequence of the initiation and evolution of the 1997/98 El Nino. The data have beet complementary to the subsurface TAO moored data in their spatial resolution and extent. The impact of satellite observations on seasonal prediction in the tropical Pacific using a coupled ocean-atmosphere general circulation model will be presented.

Rienecker, Michele M.↗

NASA NTRS1999

Mechanism for Surface Warming in the Equatorial Pacific during 1994-95

Mechanisms controlling the variation in sea surface temperature warm event in the equatorial Pacific were investigated through ocean model simulations. In addition, the mechanisms of the climatological SST cycle were investigated. The dominant mechanisms governing the seasonal cycle of SST vary significantly across the basin. In the western Pacific the annual cycle of SST is primarily in response to external heat flux. In the central basin the magnitude of zonal advection is comparable to that of the external heat flux. In the eastern basin the role of zonal advection is reduced and the vertical mixing is more important. In the easternmost equatorial Pacific the vertical entrainment contribution is as large as that of vertical diffusion. The model estimate of the vertical mixing contribution to the mixed layer heat budget compared well with estimates obtained by analysis of observations using the same diagnostic vertical mixing scheme. During 1994- 1995 the largest positive SST anomaly was observed in the mid-basin and was related to reduced latent heat flux due to weak surface winds. In the western basin the initial warming was related to enhanced external heating and reduced cooling effects of both vertical mixing and horizontal advection associated with weaker than usual wind stress. In the eastern Pacific where winds were not significantly anomalous throughout 1994-1995, only a moderate warm surface anomaly was detected. This is in contrast to strong El Nino events where the SST anomaly is largest in the eastern basin and, as shown by previous studies, the anomaly is due to zonal advection rather than anomalous surface heat flux. The end of the warm event was marked by cooling in July 1995 everywhere across the equatorial Pacific.

Rienecker, Michele M.↗

NASA NTRS1991

Chemical and biological structure and transport of a cool filament associated with a jet-eddy system of northern California in July 1986 (OPTOMA21)

The distributions of nutrient, pigment, bio-optical, and physical variables were mapped in a jet-eddy system off Point Reyes and Point Arena, California, from July 7 to 19, 1986, in order to describe the 3D variability of the filament and its relation to the nutrient and phytoplankton distributions offshore, to examine the interaction between the filament and coastal water, and to estimate the transport of nutrients and phytoplankton by the jet system. Several cool filaments were distinguishable at distances of more than 35-50 km from the coast in satellite imagery during this period. The juxtaposition of these features as well as the presence of an offshore anticyclone and a cyclone south of the filament anchored to the coast at Point Arena led to complex patterns in all variables, aided by the apparent alongshore variability in the source of upwelled water. This structure has implications for the fluxes of organic material in the region and is probably significant in organizing the interactions among different trophic levels within the system.

Jones, Burton H.↗