Search NASA⌕ Search

SEARCH · Search NASA

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 199 records · Page 11

Assimilation of satellite color observations in a coupled ocean GCM-ecosystem model

Monthly average coastal zone color scanner (CZCS) estimates of chlorophyll concentration were assimilated into an ocean global circulation model(GCM) containing a simple model of the pelagic ecosystem. The assimilation was performed in the simplest possible manner, to allow the assessment of whether there were major problems with the ecosystem model or with the assimilation procedure. The current ecosystem model performed well in some regions, but failed in others to assimilate chlorophyll estimates without disrupting important ecosystem properties. This experiment gave insight into those properties of the ecosystem model that must be changed to allow data assimilation to be generally successful, while raising other important issues about the assimilation procedure.

Sarmiento, Jorge L.↗

Assimilation of solids during ascent of magmas from the Bartoy Field of the Baikal Region, Siberia

Most investigators ascribe mare basalt magma genesis to partial melting at depths of approximately 130 to greater than 400 km within the cumulate pile deposited from a lunar magma ocean. Mare basalts share with mid-ocean ridge basalts the characteristic of relative depletion in LREE and other incompatible trace elements that arises from melting within 'used' mantle, from which crust-forming elements have already been separated. Some mare basalt types do not show the classical, La-Nd depleted mare basalt REE distributions; however, some types are isotopically heterogeneous. These differences have been ascribed to assimilation, mainly AFC-style, of KREEPy highland material overlying the source region. Might such assimilation occur during magma ascent through the KREEPy material? To gain information from a terrestrial setting on possible assimilation during ascent, we have studied a suite of Quaternary nepheline-hawalites and nepheline-mugearites from the Bartoy cinder cone complex of the Baikal Rift, Siberia. The Bartoy magmas originated from greater than 80 km deep, and erupted through thick Archean crust. We find evidence for assimilation of approximately 31 wt. percent xenocrysts of garnet, aluminous clinopyroxene, kaersutite, and olivine, all presumably from the basalt source region, but no appreciable assimilation of overlying crust, consistent with isotopic constraints. Magmatic superheat made available by rapid ascent and decomposition accounts adequately for the energy of assimilation; no accompanying fractional crystallization is required or evident.

Haas, Johnson R.↗

Improving Assimilated Global Data Sets using TMI Rainfall and Columnar Moisture Observations

A global analysis that optimally combine observations from diverse sources with physical models of atmospheric and land processes can provide a comprehensive description of the climate systems. Currently, such data products contain significant errors in primary hydrological fields such as precipitation and evaporation, especially in the tropics. In this study, we show that assimilating precipitation and total precipitable water (TPW) retrievals derived from the TRMM Microwave Imager (TMI) improves not only the hydrological cycle but also key climate parameters such as clouds, radiation, and the large-scale circulation produced by the Goddard Earth Observing System (GEOS) data assimilation system (DAS). In particular, assimilating TMI rain improves clouds and radiation in areas of active convection, as well as the latent heating distribution and the large-scale motion field in the tropics, while assimilating TMI TPW heating distribution and the large-scale motion field in the tropics, while assimilating TMI TPW retrievals leads to reduced moisture biases and improved radiative fluxes in clear-sky regions. The improved analysis also improves short-range forecasts in the tropics. Ensemble forecasts initialized with the GEOS analysis incorporating TMI rain rates and TPW yield smaller biases in tropical precipitation forecasts beyond 1 day and better 500 hPa geopotential height forecasts up to 5 days. Results of this study demonstrate the potential of using high-quality space-borne rainfall and moisture observations to improve the quality of assimilated global data for climate analysis and weather forecasting applications

Hou, Arthur Y.↗

A Technique for Assimilating GOES-Derived Land Surface Products into Regional Models to Improve the Representation of Land Surface Forcing

As the parameterizations of surface energy budgets in regional models have become more complete physically, models have the potential to be much more realistic in simulations of coupling between surface radiation, hydrology, and surface energy transfer. Realizing the importance of properly specifying the surface energy budget, many institutions are using land-surface models to represent the lower boundary forcing associated with biophysical processes and soil hydrology. However, the added degrees of freedom due to inclusion of such land-surface schemes require the specification of additional parameters within the model system such as vegetative resistances, green vegetation fraction, leaf area index, soil physical and hydraulic characteristics, stream flow, runoff, and the vertical distribution of soil moisture. A technique has been developed for assimilating GOES-IR skin temperature tendencies 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 is that it does not require a complex land-surface formulation within the atmospheric model. As a result, the need to specify poorly known soil and vegetative characteristics is eliminated. The GOES assimilation technique has been incorporated into the PSU/NCAR MM5. Results will be presented to demonstrate the ability of the assimilation scheme to improve short-term simulations of near-surface air temperature and mixing ratio during the warm season for several selected cases during June of 1995 and June 1997 which exhibit a variety of atmospheric and land-surface conditions. In addition, the simulations produced by the assimilation technique will be compared with those produced by the MM5/BATS coupled model initialized with soil water content derived from a soil hydrology model.

Lapenta, William M.↗

Results from Global Land-Surface Data Assimilation Methods

Realistic representation of the land surface is crucial in global climate modeling (GCM). Recently, the Mosaic land-surface Model (LSM) has been driven off-line using GEOS DAS (Goddard Earth Observing System Data Assimilation System) atmospheric forcing, forming the Off-line Land-surface Global Assimilation (OLGA) system. This system provides a computationally efficient test bed for land surface data assimilation. Here, we validate the OLGA simulation of surface processes and the assimilation of ISCCP surface temperatures. Another component of this study as the incorporation of the Physical-space Statistical Analysis System (PSAS) into OLGA, in order to assimilate surface temperature observations from the International Satellite Cloud Climatology Project (ISCCP). To counteract the subsequent forcing of the analyzed skin temperature back to the initial state following the analysis. incremental bias correction (IBC) was included in the assimilation. The IBC scheme effectively removed the time mean bias, but did not remove him in the mean diurnal cycle. Therefore, a diurnal him correction (DBC) scheme was developed, where the time-dependent bias was modeled with a sine wave parameterization. In addition, quality control of the ISCCP data and anisotropic temperature correction were implemented in PSAS. Preliminary results showed a substantial impact from the inclusion of PSAS and DBC that was visible in the surface meteorology fields and energy budget. Also, the monthly mean diurnal cycle from the experiment closely matched the diurnal cycle from the observations.

Radakovich, Jon D.↗

The Evolution of Stratospheric Data Assimilation Systems

The use of model-assimilated meteorological observations for stratospheric research has become routine since the late 1980's. The first stratospheric assimilation systems were straightforward extensions of systems developed for tropospheric weather forecasting. During the 1990's systems were developed that more directly addressed the specifics of stratospheric applications. These developments include better treatment of the satellite observations and improved models that better represent the residual circulation in the assimilated data sets. This talk will review the evolution of stratospheric data assimilation and its application, especially to problems of tracer transport. The new data assimilation currently under validation at NASA will be described in some detail, and results from the validation exercise will be presented. This data assimilation system sits at the foundation of a proposed stratospheric reanalysis that covers the era of the Upper Atmosphere Research Satellite (UARS).

Rood, Richard B.↗

Plans for Assimilating AQUA data at NASA's DAO

NASA's Data Assimilation Office (DAO) is expanding its work with TIROS Operational Vertical Sounder (TOVS) to assimilate data from the advanced instruments which will fly on NASA's AQUA satellite in early 2002. The Atmospheric Infrared Sounder (AIRS), which has over 2000 channels, together with the Advanced Microwave Sounding Unit-A (AMSU-A) and the Humidity Sounder Brazil (HSB) will provide many technical challenges for data assimilation centers. One of the primary concerns is how best to subset the data in order to efficiently extract information about the Earth's atmosphere and surface. This includes static and dynamic channel selection as well as pixel thinning. The DAO is currently experimenting with simulated AIRS/AMSU/HSB radiances within the framework of our finite volume data assimilation system (NDAS) using the OPTRAN radiative transfer code developed as part of the NOAA/NASA Joint Center for Satellite Data Assimilation. The short-term goals include the assessment of the cost of processing various data subsets and preparation for near-real-time assimilation within a few months of launch. We will also discuss plans and tools for evaluating the quality of AIRS data including radiances and level 2 products from the AIRS science team.

Frank, Donald↗

The Computational Complexity, Parallel Scalability, and Performance of Atmospheric Data Assimilation Algorithms

The computational complexity of algorithms for Four Dimensional Data Assimilation (4DDA) at NASA's Data Assimilation Office (DAO) is discussed. In 4DDA, observations are assimilated with the output of a dynamical model to generate best-estimates of the states of the system. It is thus a mapping problem, whereby scattered observations are converted into regular accurate maps of wind, temperature, moisture and other variables. The DAO is developing and using 4DDA algorithms that provide these datasets, or analyses, in support of Earth System Science research. Two large-scale algorithms are discussed. The first approach, the Goddard Earth Observing System Data Assimilation System (GEOS DAS), uses an atmospheric general circulation model (GCM) and an observation-space based analysis system, the Physical-space Statistical Analysis System (PSAS). GEOS DAS is very similar to global meteorological weather forecasting data assimilation systems, but is used at NASA for climate research. Systems of this size typically run at between 1 and 20 gigaflop/s. The second approach, the Kalman filter, uses a more consistent algorithm to determine the forecast error covariance matrix than does GEOS DAS. For atmospheric assimilation, the gridded dynamical fields typically have More than 10(exp 6) variables, therefore the full error covariance matrix may be in excess of a teraword. For the Kalman filter this problem can easily scale to petaflop/s proportions. We discuss the computational complexity of GEOS DAS and our implementation of the Kalman filter. We also discuss and quantify some of the technical issues and limitations in developing efficient, in terms of wall clock time, and scalable parallel implementations of the algorithms.

Lyster, Peter M.↗

Advances in Land Data Assimilation Systems

Assimilation of remotely sensed land surface observations into regional to global scale numerical models have the potential to significantly advance our ability, to assess, understand, and predict surface water, energy, and carbon cycles. This session seeks to assess the state-of-the-art in data assimilation methods for integrating land surface remote sensing and modeling, with a focus on practical applications and techniques. Assimilated land surface variables of interest include (but are not limited to, soil moisture, surface temperature, snowpack, streamflow, vegetation dynamics, and carbon storage. Contributions describing the development of practical land surface data assimilation methods, multivariate land surface data assimilation strategies, evaluation of the required accuracy and resolution of remote sensing observations, the effects of scale, process complexity, and uncertainty on data assimilation, and the optimal treatment of model and observation errors are encouraged.

Houser, Paul R.↗

The GEOS Retrospective Data Assimilation System: The 6-hour lag case

The fixed-lag Kalman smoother (FLKS) has been proposed as a framework to construct data assimilation procedures capable of producing high-quality climate research datasets. Fixed-lag Kalman smoother-based systems, referred to as retrospective data assimilation systems, are an extension to three-dimensional filtering procedures with the added capability of incorporating observations not only in the past and present time of the estimate, but also at future times. A variety of simplifications are necessary to render retrospective assimilation procedures practical. In this article, we present an FLKS-based retrospective data assimilation system implementation for the Goddard Earth Observing System (GOES) Data Assimilation System (DAS). The practicality of this implementation comes from the practicality of its underlying (filter) analysis system, i.e., the physical-space statistical analysis system (PSAS). The behavior of two schemes is studied here. The first retrospective analysis (RA) scheme is designed simply to update the regular PSAS analyses with observations available at times ahead of the regular analysis times. Although our GEOS DAS implementation is general, results are only presented for when observations 6-hours ahead of the analysis time are used to update the PSAS analyses and thereby to calculate the so-called lag-1 retrospective analyses. Consistency tests for this RA scheme show that the lag-1 retrospective analyses indeed have better 6-hour predictive skills than the predictions from the regular analyses. This motivates the introduction of the second retrospective analysis scheme which, at each analysis time, uses the 6-hour retrospective analysis to replace the first-guess normally used in the PSAS analysis, and therefore allows the calculation of a revised (filter) PSAS analysis. Since in this scheme the lag-1 retrospective analyses influence the filter results, this procedure is referred to as the retrospective-based iterative analysis (RIA) scheme. Results from the RIA scheme indicate its potential for improving the overall quality of the assimilation.

Zhu, Yan-Qiu↗

High-Frequency Planetary Waves in the Polar Middle Atmosphere as seen in a data Assimilation System

This study examines the winter southern hemisphere vortex of 1998 using four times daily output from a data assimilation system to focus on the polar 2-day, wave number 2 component of the 4-day wave. The data assimilation system products are from a test version of the finite volume data assimilation system (fvDAS) being developed at Goddard Space Flight Center (GSFC) and include an ozone assimilation system. Results show that the polar 2-day wave dominates during July 1998 at 70 degrees. The period of the quasi 2-day wave is somewhat shorter than 2 days (about 1.7 days) during July 1998 with an average perturbation temperature amplitude for the month of over 2.5 K. The 2-day wave propagates more slowly than the zonal mean zonal wind, consistent with Rossby wave theory, and has EP flux divergence regions associated with regions of negative horizontal potential vorticity gradients, as expected from linear instability theory. Results for the assimilation-produced ozone mixing ratio show that the 2-day wave represents a major source of ozone variation in this region. The ozone wave in the assimilation system is in good agreement with the wave seen in the POAM (Polar Ozone and Aerosol Measurement) ozone observations for the same time period. Some differences with linear instability theory are noted as well as spectral peaks in the ozone field, not seen in the temperature field, that may be a consequence of advection.

Coy, L.↗

Sensitivity of Assimilated Tropical Tropospheric Ozone to the Meteorological Analyses

Tropical tropospheric ozone fields from two different experiments performed with an off-line ozone assimilation system developed in NASA's Data Assimilation Office (DAO) are examined. Assimilated ozone fields from the two experiments are compared with the collocated ozone profiles from the Southern Hemispheric Additional Ozonesondes (SHADOZ) network. Results are presented for 1998. The ozone assimilation system includes a chemistry-transport model, which uses analyzed winds from the Goddard Earth Observing System (GEOS) Data Assimilation System (DAS). The two experiments use wind fields from different versions of GEOS DAS: an operational version of the GEOS-2 system and a prototype of the GEOS-4 system. While both versions of the DAS utilize the Physical-space Statistical Analysis System and use comparable observations, they use entirely different general circulation models and data insertion techniques. The shape of the annual-mean vertical profile of the assimilated ozone fields is sensitive to the meteorological analyses, with the GEOS-4-based ozone being closest to the observations. This indicates that the resolved transport in GEOS-4 is more realistic than in GEOS-2. Remaining uncertainties include quantification of the representation of sub-grid-scale processes in the transport calculations, which plays an important role in the locations and seasons where convection dominates the transport.

Hayashi, Hiroo↗

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

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↗

Forward-looking Assimilation of MODIS-derived Snow Covered Area into a Land Surface Model

Snow cover over land has a significant impact on the surface radiation budget, turbulent energy fluxes to the atmosphere, and local hydrological fluxes. For this reason, inaccuracies in the representation of snow covered area (SCA) within a land surface model (LSM) can lead to substantial errors in both offline and coupled simulations. Data assimilation algorithms have the potential to address this problem. However, the assimilation of SCA observations is complicated by an information deficit in the observation SCA indicates only the presence or absence of snow, and not snow volume and by the fact that assimilated SCA observations can introduce inconsistencies with atmospheric forcing data, leading to non-physical artifacts in the local water balance. In this paper we present a novel assimilation algorithm that introduces MODIS SCA observations to the Noah LSM in global, uncoupled simulations. The algorithm utilizes observations from up to 72 hours ahead of the model simulation in order to correct against emerging errors in the simulation of snow cover while preserving the local hydrologic balance. This is accomplished by using future snow observations to adjust air temperature and, when necessary, precipitation within the LSM. In global, offline integrations, this new assimilation algorithm provided improved simulation of SCA and snow water equivalent relative to open loop integrations and integrations that used an earlier SCA assimilation algorithm. These improvements, in turn, influenced the simulation of surface water and energy fluxes both during the snow season and, in some regions, on into the following spring.

Zaitchik, Benjamin F.↗

Assimilation of SCIAMACHY Total Column CO Observations: Regional Analysis of Data Impact

Carbon monoxide (CO) total column observations from the SCanning Imaging Absorption SpectroMeter for Atmospheric CHartography (SCIAMACHY) on board ENVISAT are assimilated into the Global Modeling and Assimilation Office (GMAO) constituent assimilation system for the period July 18-October 31, 2004. This is the first assimilation of CO observations from a near infrared sounder. The impact of the assimilation on CO distribution is evaluated using independent Measurement of Ozone and Water vapor by Airbus In-service Aircraft (MOZAIC) in-situ CO profiles. Assimilation of satellite data improves agreement with MOZAIC CO globally, especially in the upper troposphere.

Tangborn, Andrew↗

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↗

Lessons Learned from Assimilating Altimeter Data into a Coupled General Circulation Model with the GMAO Augmented Ensemble Kalman Filter

Satellite altimetry measurements have provided global, evenly distributed observations of the ocean surface since 1993. However, the difficulties introduced by the presence of model biases and the requirement that data assimilation systems extrapolate the sea surface height (SSH) information to the subsurface in order to estimate the temperature, salinity and currents make it difficult to optimally exploit these measurements. This talk investigates the potential of the altimetry data assimilation once the biases are accounted for with an ad hoc bias estimation scheme. Either steady-state or state-dependent multivariate background-error covariances from an ensemble of model integrations are used to address the problem of extrapolating the information to the sub-surface. The GMAO ocean data assimilation system applied to an ensemble of coupled model instances using the GEOS-5 AGCM coupled to MOM4 is used in the investigation. To model the background error covariances, the system relies on a hybrid ensemble approach in which a small number of dynamically evolved model trajectories is augmented on the one hand with past instances of the state vector along each trajectory and, on the other, with a steady state ensemble of error estimates from a time series of short-term model forecasts. A state-dependent adaptive error-covariance localization and inflation algorithm controls how the SSH information is extrapolated to the sub-surface. A two-step predictor corrector approach is used to assimilate future information. Independent (not-assimilated) temperature and salinity observations from Argo floats are used to validate the assimilation. A two-step projection method in which the system first calculates a SSH increment and then projects this increment vertically onto the temperature, salt and current fields is found to be most effective in reconstructing the sub-surface information. The performance of the system in reconstructing the sub-surface fields is particularly impressive for temperature, but not as satisfactory for salt.

Keppenne, Christian↗