Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “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 523 records · Page 29

I/O Parallelization for the Goddard Earth Observing System Data Assimilation System (GEOS DAS)

The National Aeronautics and Space Administration (NASA) Data Assimilation Office (DAO) at the Goddard Space Flight Center (GSFC) has developed the GEOS DAS, a data assimilation system that provides production support for NASA missions and will support NASA's Earth Observing System (EOS) in the coming years. The GEOS DAS will be used to provide background fields of meteorological quantities to EOS satellite instrument teams for use in their data algorithms as well as providing assimilated data sets for climate studies on decadal time scales. The DAO has been involved in prototyping parallel implementations of the GEOS DAS for a number of years and is now embarking on an effort to convert the production version from shared-memory parallelism to distributed-memory parallelism using the portable Message-Passing Interface (MPI). The GEOS DAS consists of two main components, an atmospheric General Circulation Model (GCM) and a Physical-space Statistical Analysis System (PSAS). The GCM operates on data that are stored on a regular grid while PSAS works with observational data that are scattered irregularly throughout the atmosphere. As a result, the two components have different data decompositions. The GCM is decomposed horizontally as a checkerboard with all vertical levels of each box existing on the same processing element(PE). The dynamical core of the GCM can also operate on a rotated grid, which requires communication-intensive grid transformations during GCM integration. PSAS groups observations on PEs in a more irregular and dynamic fashion.

Lucchesi, Rob↗

Implementation of a Parallel Kalman Filter for Stratospheric Chemical Tracer Assimilation

A Kalman filter for the assimilation of long-lived atmospheric chemical constituents has been developed for two-dimensional transport models on isentropic surfaces over the globe. An important attribute of the Kalman filter is that it calculates error covariances of the constituent fields using the tracer dynamics. Consequently, the current Kalman-filter assimilation is a five-dimensional problem (coordinates of two points and time), and it can only be handled on computers with large memory and high floating point speed. In this paper, an implementation of the Kalman filter for distributed-memory, message-passing parallel computers is discussed. Two approaches were studied: an operator decomposition and a covariance decomposition. The latter was found to be more scalable than the former, and it possesses the property that the dynamical model does not need to be parallelized, which is of considerable practical advantage. This code is currently used to assimilate constituent data retrieved by limb sounders on the Upper Atmosphere Research Satellite. Tests of the code examined the variance transport and observability properties. Aspects of the parallel implementation, some timing results, and a brief discussion of the physical results will be presented.

Chang, Lang-Ping↗

I/O Parallelization for the Goddard Earth Observing System Data Assimilation System (GEOS DAS)

The National Aeronautics and Space Administration (NASA) Data Assimilation Office (DAO) at the Goddard Space Flight Center (GSFC) has developed the GEOS DAS, a data assimilation system that provides production support for NASA missions and will support NASA's Earth Observing System (EOS) in the coming years. The DAO's support of the EOS project along with the requirement of producing long-term reanalysis datasets with an unvarying system levy a large I/O burden on the future system. The DAO has been involved in prototyping parallel implementations of the GEOS DAS for a number of years and is now converting the production version from shared-memory parallelism to distributed-memory parallelism using the portable Message-Passing Interface (MPI). If the MPI-based GEOS DAS is to meet these production requirements, we must make I/O from the parallel system efficient. We have designed a scheme that allows efficient I/O processing while retaining portability, reducing the need for post-processing, and producing data formats that are required by our users, both internal and external. The first phase of the GEOS DAS Parallel I/O System (GPIOS) will expand upon the common method of gathering global data to a Single PE for output. Instead of using a PE also tasked with primary computation, a number of PEs will be dedicated to I/O and its related tasks. This allows the data transformations and formatting required prior to output to take place asynchronously with respect to the GEOS DAS assimilation cycle, improving performance and generating output data sets in a format convenient for our users. I/O PEs can be added as needed to handle larger data volumes or to meet user file specifications. We will show I/O performance results from a prototype MPI GCM integrated with GPIOS. Phase two of GPIOS development will examine ways of integrating new software technologies to further improve performance and build scalability into the system. The maturing of MPI-IO implementations and other supporting libraries such as parallel HDF should provide performance gains while retaining portability.

Lucchesi, R.↗

Impact of Assimilating GOES-Derived Land Surface Variables into the PSU/NCAR MM5

Land surface processes are known to have a profound impact on the overlying atmosphere over a wide range of spatial and temporal scales. Many atmospheric numerical models include special parameterizations to improve the specification and partitioning of surface fluxes which are critical to the accurate prediction of warm season boundary layer behavior, organized mesoscale circulations, and convective precipitation. However, the added degrees of freedom resulting from the inclusion of vegetation and soil schemes require the specification of additional surface parameters such as vegetative resistances, green vegetation fraction, leaf area index, soil physical and hydraulic characteristics, and the vertical distribution of soil moisture. As satellite data have become more readily available in recent years, many investigations have attempted to use these new measurements to infer missing components of the surface energy budget. Sensitivity studies have shown land-skin temperature (LST) tendencies during the mid-morning hours are strongly sensitive to the surface moisture availability (a function of soil wetness and vegetation) and less sensitive to other parameters such as surface roughness. Based upon results from these and other studies, developed a simple technique that dynamically assimilates Geostationary Operational Environmental Satellite (GOES) derived land-surface products into the surface energy budget of a mesoscale model. The purpose of this paper is to demonstrate that assimilating the GOES satellite data has the potential to improve the representation of land surface characteristics within the model without prior knowledge of the land surface characteristics. The assimilation technique is presented in Section 2 and the numerical experiments are detailed in Section 3. Preliminary results and conclusions are presented in Sections 4 and 5, respectively.

Lapenta, William M.↗

Some Continuum Aspects of Data Assimilation

A long-sought goal in data assimilation is to build more fully the time dimension into the data assimilation process. Truly four-dimensional data assimilation requires evolving second-moment information, namely the estimation error covariance, along with the state estimate itself. The continuum evolution equations for the estimation error covariance constitute a system of partial differential equations in six space dimensions, forced by a model error covariance and by the observations themselves. The high dimensionality of this system poses a difficult computational problem. Numerous methods have been proposed for approximate, discrete solution of this system of equations. Four-dimensional variational schemes solve these equations indirectly, while schemes based on Kalman filtering solve the equations more directly, usually on a small subspace of the full six-dimensional space. In both cases, most proposed solution methods are not derived from the continuum covariance evolution equations themselves. Instead, most methods simply inherit the discretization used for evolving the state estimate, either in the form of a discrete tangent linear model or, in the case of ensemble schemes, in the form of a fully nonlinear discrete model. In this lecture we show that solutions of the continuum covariance evolution equations possess simple properties that are not readily satisfied upon discretization of the equations. we give numerical examples illustrating that "inherited" discretizations sometimes fail to produce meaningfully accurate solutions. Finally, we suggest alternative discretization methods that may yield more faithful approximate solutions of the continuum problem.

Cohn, Stephen E.↗

Satellite Derived Land Surface Temperature for Model Assimilation

Studies have shown that land surface temperature (LST) tendencies are sensitive to the surface moisture availability which is a function of soil moisture and vegetation. The assimilation of satellite derived LST tendencies into the surface energy budget of mesoscale models has shown promise in improving the representation of the complex effects of both soil moisture and vegetation within the models for short term simulations. LST derived from geostationary satellites has the potential of providing the temporal and spatial resolution needed for an LST assimilation process. This paper presents an analysis comparing the LST derived from GOES-8 infrared measurements with LST calculated by the MM5 numerical model. The satellite derived LSTs are calculated using a physical split window approach using channels 4 and 5 of GOES-8. The differences in the LST data sets, especially the tendencies, are presented and examined. Quantifying the differences between the data sets provide insight of possible weaknesses in the model parameterizations affecting the surface energy budget calculations and an indication of the potential effectiveness o f assimilating LST into the models.

Suggs, Ronnie J.↗

Middle Atmosphere Transport Properties of Assimilated Datasets

One of the most compelling reasons for performing data assimilation in the middle atmosphere is to obtain global, balanced datasets for studies of trace gas transport and chemistry. This is a major motivation behind the Goddard Earth observation System-Data Assimilation System (GEOS-DAS). Previous studies have shown that while this and other data assimilation systems can generally obtain good estimates of the extratropical rotational velocity field, the divergent part of the dynamical field is deficient; this impacts the "residual circulation" and leads to spurious trace gas transport on seasonal and interannual timescales. These problems are impacted by the quality and the method of use of the observational data and by deficiencies in the atmospheric general circulation model. Whichever the cause at any place and time, the "solution" is to introduce non-physical forcing terms into the system (the so-called incremental analysis updates); these can directly (thermal) or indirectly (mechanical) affect the residual circulation. This paper will illustrate how the divergent circulation is affected by deficiencies in both observations and models. Theoretical considerations will be illustrated with examples from the GEOS-DAS and from simplified numerical experiments. These are designed to isolate known problems, such as the inability of models to sustain a quasi-biennial oscillation and sparse observational constraints on tropical dynamics, or radiative inconsistencies in the presence of volcanic aerosols.

Pawson, Steven↗

Massively Parallel Assimilation of TOGA/TAO and Topex/Poseidon Measurements into a Quasi Isopycnal Ocean General Circulation Model Using an Ensemble Kalman Filter

A massively parallel ensemble Kalman filter (EnKF)is used to assimilate temperature data from the TOGA/TAO array and altimetry from TOPEX/POSEIDON into a Pacific basin version of the NASA Seasonal to Interannual Prediction Project (NSIPP)ls quasi-isopycnal ocean general circulation model. The EnKF is an approximate Kalman filter in which the error-covariance propagation step is modeled by the integration of multiple instances of a numerical model. An estimate of the true error covariances is then inferred from the distribution of the ensemble of model state vectors. This inplementation of the filter takes advantage of the inherent parallelism in the EnKF algorithm by running all the model instances concurrently. The Kalman filter update step also occurs in parallel by having each processor process the observations that occur in the region of physical space for which it is responsible. The massively parallel data assimilation system is validated by withholding some of the data and then quantifying the extent to which the withheld information can be inferred from the assimilation of the remaining data. The distributions of the forecast and analysis error covariances predicted by the ENKF are also examined.

Keppenne, Christian L.↗

Middle Atmospheric Transport Properties of Assimilated Datasets

One of the most compelling reasons for performing data assimilation in the middle atmosphere is to obtain global, balanced datasets for studies of trace gas transport and chemistry. This is a major motivation behind the Goddard Earth observation System-Data Assimilation System (GEOS-DAS). Previous studies have shown that while this and other data assimilation systems can generally obtain good estimates of the extratropical rotational velocity field, the divergent part of the dynamical field is deficient; this impacts the "residual circulation" and leads to spurious trace gas transport on seasonal and interannual timescales. These problems are impacted by the quality and the method of use of the observational data and by deficiencies in the atmospheric general circulation model. Whichever the cause at any place and time, the "solution" is to introduce non-physical forcing terms into the system (the so-called incremental analysis updates); these can directly (thermal) or indirectly (mechanical) affect the residual circulation. This paper will illustrate how the divergent circulation is affected by deficiencies in both observations and models. Theoretical considerations will be illustrated with examples from the GEOS-DAS and from simplified numerical experiments. These are designed to isolate known problems, such as the inability of models to sustain a quasi-biennial oscillation and sparse observational constraints on tropical dynamics, or radiative inconsistencies in the presence of volcanic aerosols.

Pawson, Steven↗

Simulated Surface Energy Budgets Over the Southeastern US: The GHCC Satellite Assimilation System and the NCEP Early Eta

A technique has been developed for assimilating GOES-derived skin temperature tendencies and insolation into the surface energy budget equation of a mesoscale model so that the simulated rate of temperature change closely agrees with the satellite observations. A critical assumption of the technique is that the availability of moisture (either from the soil or vegetation) is the least known term in the model's surface energy budget. Therefore, the simulated latent heat flux, which is a function of surface moisture availability, is adjusted based upon differences between the modeled and satellite-observed skin temperature tendencies. An advantage of this technique is that satellite temperature tendencies are assimilated in an energetically consistent manner that avoids energy imbalances and surface stability problems that arise from direct assimilation of surface shelter temperatures. The fact that the rate of change of the satellite skin temperature is used rather than the absolute temperature means that sensor calibration is not as critical. An advantage of this technique for short-range forecasts (0-48h) is that it does not require a complex land-surface formulation within the atmospheric model. As a result, we can avoid having to specify land surface characteristics such as vegetation resistances, green fraction, leaf area index, soil physical and hydraulic characteristics, stream flow, runoff, and the vertical and horizontal distribution of soil moisture.

Lapenta, William M.↗

Comparison of Radiative Forcing Calculations Due to Mineral Dust from a Transport Model, Satellite Measurements and an Assimilation System

This study uses information on mineral aerosol from a transport model to calculate global radiative forcing values. The transport model is driven by assimilated meteorology and outputs three-dimensional dust spatial information for various size ranges. The dust fields are input to an off-line radiative transfer calculation to obtain the direct radiative forcing due to the dust fields. During June, July and August of 1988 presence of dust 1) reduces the global net incoming radiation at the top of atmosphere (TOA) by 0.3 to 0.7 W/sq m and 2) reduces net incoming radiation at the earth's surface by 1.3 to 2.0 W/sq m. Over Africa our estimates of the reduction of radiation at the top of atmosphere compare well with TOA reductions derived from ERBE and TOMS satellite data. However, our heating rates are not consistent with analysis temperature increments produced by the assimilation system over regions of high aerosol loading. These increments are based on differences between temperature observations and temperatures from the assimilation general circulation model. One explanation is that the lower tropospheric temperatures retrieved by TOVS are being contaminated by mineral aerosol.

Weaver, Clark J.↗

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.↗

Data Assimilation in the Ocean and in the Atmosphere: What Should be Next?

As part of the International Symposium on Assimilation of Observation in Meteorology and Oceanography, a panel discussion was held on the evening of 15 March 1995. The purpose of this panel discussion was focus on what the next major areas of research in data assimilation should be. The panelists had five minutes each for short presentations (Kalman filters, representers, etc.) and this was followed by an open discussion. This preprocessing will require a good understanding of the fine-scale phenomena. Least square methods such as Kalman filters and variational schemes are inefficient estimators of non-Gaussian field such as chemical traces. Regardless of the modeling technique employed (Lagrangian methods seem best), a least squares assimilation scheme will smear fine structure. The estimator of maximum likelihood must be sought, by examination of tracer probability distributions.

Kalnay, Eugenia↗

An Extended Kalman Filter to Assimilate Altimetric Data into a Non-Linear Model of the Tropical Pacific

A new implementation of the extended Kaman filter is developed for the purpose of assimilating altimetric observations into a primitive equation model of the tropical Pacific. Its specificity consists in defining the errors into a reduced basis that evolves in time with the model dynamic. Validation by twin experiments is conducted and the method is shown to be efficient in quasi real conditions. Data from the first 2 years of the Topex/Poseidon mission are assimilated into the Gent & Cane [1989] model. Assimilation results are evaluated against independent in situ data, namely TAO mooring observations.

Gourdeau, L.↗

NOAA/NASA/DOD Workshop on Satellite Data Assimilation

A workshop on the assimilation of satellite sounding information using global forecast and climate models was held at College Park, MD, 23-25 August 1999. Topics discussed included: comparisons of assimilations of satellite retrievals versus satellite-observed radiances, planning for the use of advanced infrared sounders, the use of satellite sounding data affected by land surfaces, radiative transfer issues, and error characteristics of models and observations. The workshop concluded with a number of general and specific recommendations to advance the state of the art of assimilation of satellite sounding data.

Errico, Ronald M.↗

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↗

Assimilation of IASI and AIRS Data: Information Content and Quality Control

The Infrared Atmospheric Sounding Interferometer (IASI) and Atmospheric Infrared Sounder (AIRS) instruments have two orders of magnitude more channels that the current operational infrared sounder (High Resolution Infra-Red Sounder (HIRS)). This data volume presents a technological challenge for using the data in a data assimilation system. Data reduction will be a necessary for assimilation. It is important to understand the information content of the radiance measurements for data reduction purposes. In this talk, I will discuss issues relating to information content and quality control for assimilation of the AIRS and IASI data.

Joiner, J.↗

DAO's Finite-volume/Physical-space Data Assimilation System: Stratospheric Applications

In this talk we describe the next-generation data assimilation system being developed at NASA's Data Assimilation Office (DAO), with emphasis on the applications to stratospheric forecasts and stratospheric constituent transport. This data assimilation system includes the General Circulation Model jointly developed by the DAO and the Climate and Global Dynamics Division (CGDD) at NCAR. This model is based on the finite-volume dynamical core) developed at DAO with physical parameterizations from the NCAR Climate Community Model. The Physical-space Statistical Analysis System (PSAS) is used to combine a first guess from the NASA-NCAR GCM with observational data to provide an updated estimate of the state of the atmosphere. Case studies for the (northern) Winter of 2000 will be discussed.

daSilva, Arlindo↗