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 1,567 records · Page 87

One-to-one aeroservoelastic validation of operational loads and performance of a 2.8 MW wind turbine model in OpenFAST

Abstract. This article presents a validation study of the popular aeroservoelastic code suite OpenFAST leveraging weeks of measurements obtained during normal operation of a 2.8 MW land-based wind turbine. Measured wind conditions were used to generate one-to-one turbulent flow fields (i.e., comparing simulation to measurement in 10 min increments, or bins) through unconstrained and constrained assimilation methods using the kinematic turbulence generators TurbSim and PyConTurb. A total of 253 bins of 10 min of normal turbine operation were selected for analysis, and a statistical comparison in terms of performance and loads is presented. We show that successful validation of the model was not strongly dependent on the type of inflow assimilation method used for mean quantities of interest, which had median modeling errors per wind-speed interval generally within 5 %–10 % of the measurement. The type of inflow assimilation method did have a larger effect on the fatigue predictions for blade-root flapwise and tower-base fore–aft quantities, which surprisingly saw larger errors from the assumed higher-fidelity assimilation methods. Avenues for further work are discussed and include possible improvements to the aerodynamic, structural, and controller modeling that may offer insight on the origin of the up to ∼ 40 % median overprediction of fatigue for these quantities.

17 WIND ENERGY↗

Analysis and forecast experiments incorporating satellite soundings and cloud and water vapor drift wind information

A system for assimilating conventional meteorological data and satellite-derived data in order to produce four-dimensional gridded data sets of the primary atmospheric variables used for updating limited area forecast models is described. The basic principles of a data assimilation scheme as proposed by Lorenc (1984) are discussed. The design of the system and its incremental assimilation cycles are schematically presented. The assimilation system was tested using radiosonde, buoy, VAS temperature, dew point, gradient wind data, cloud drift, and water vapor motion data. The rms vector errors for the data are analyzed.

Goodman, Brian M.↗

Evaluation of crustal recycling during the evolution of Archean-age Matachewan basaltic magmas

The simplest model for the Matachewan-Hearst Dike (MHD) magmas is assimilation-fractional crystallization (AFC), presumably occurring at the base of the crust during underplating. Subduction zone enriched mantle sources are not required. Trace elements suggest that the mantle sources for the MHD were depleted, but possessed a degree of heterogeneity. Rates of assimilation were approximately 0.5 (= Ma/Mc); the contaminant mass was less than 20 percent. The contaminant was dominated by tonalites-randodiorites, similar to xenoliths and rocks in the Kapuskasing Structural Zone (KSZ). Assimilation of partial melts of light-rare earth and garnet-bearing basaltic precursors may have produced some the MHD magmas. Apparently, previous underplating-AFC processes had already produced a thick crust. The silicic granitoid assimilant for the MHD magmas was probably produced by earlier processing of underplated mafic crust (4, 5, 10, 21 and 30). Calculations suggest that the derived silicic rocks possess negative Ta and Ti anomalies even though they were not the product of subduction.

Nelson, Dennis O.↗

The behavior of forecast error covariances for a Kalman filter in two dimensions

The behavior of forecast error covariances in a fairly realistic setting is demonstrated via a Kalman filter algorithm. It is used to assimilate simulated data from the existing radiosonde network, from the demonstration network of 31 Doppler wind profilers in the central U.S., and from hypothetical radiometers located at five of the profiler sites. Some theoretical justification of the hypothesis advanced by Phillips (1982), and the hypothesis is used to formulate the model error covariance matrix required by the Kalman filter. The results show that assimilating the profiler wind data leads to a large reduction of forecast/analysis error in heights as well as in winds, over the profiler region and also downstream, when compared with the results of assimilating the radiosonde data alone. The forecast error covariance matrices that the Kalman filter calculates to obtain this error reduction differ considerably from those prescribed by the optimal interpolation schemes that are employed for data assimilation at operational centers.

Cohn, Stephen E.↗

A Variational Formalism for the Radiative Transfer Equation: Prelude to Model 3

The MODEL III variational data assimilation model is the third of four general assimilation models designed to blend weather data measured from space based platforms in the meteorological data mainstream in a way that maximizes the information content of the satellite data. Because there are many different observation locations and there are many instruments with different measurement error characteristics, it is also necessary to require that the blending be done to maximize the information content of the data and simultaneously to retain a dynamically consistent and reasonably accurate description of the state of the atmosphere. This is ideally a variational problem for which the data receive relative weights that are inversely proportional to measurement error and are adjusted to satisfy a set of dynamical equations that govern atmospheric processes. The advantage of MODEL III over the previous two models is that radiance, the atmospheric variable measured by satellite, becomes a dependent variable. In the previous versions, mean layer temperatures that had been retrieved from the radiances by some method, were included in the assimilation by substituting them in place of the rawinsonde temperatures. Now both rawinsonde temperatures and satellite radiances are included independently in the assimilation.

Achtemeier, Gary L.↗

Retrieval of wind temperature and pressure from single Doppler radar and a numerical model

A 4D data assimilation algorithm to obtain 3D wind and thermodynamic fields from radial velocities and a numerical model is derived. It is a blend of the continuous updating technique of Charney et al. (1969) and a diagnostic pressure and temperature retrieval technique suggested by Gal-Chen (1978). This model is tested only against model-generated data which are viewed as 'data observed from a real atmosphere'. A thermal bubble developing in a dry, neutral environment is simulated. The thermal recovery technique is found to be a necessary procedure to obtain successful data assimilation results. The present assimilation method is found to be capable of reducing the 'observational error' and making the assimilation run converge toward the control run. Without data insertion, errors in the nonobserved wind component are large.

Liou, Yu-Chieng↗

Searching for neuKREEP: An EMP study of Apollo 11 Group A basalts

The Apollo 11 and 17 landing sites are characterized by the presence of high-Ti basalts (TiO2 greater than 6 percent). The Group A basalts of Apollo 11 have elevated K compositions (greater than 2000 ppm); and are enriched in incompatible trace elements relative to the other types of high-Ti basalt found in the region. These unique basalts also are the youngest of all high-Ti basalts, with an age of 3.56 +/- 0.02 Ga. Recent modelling of the Apollo 11 Group A basalts by Jerde et al. has demonstrated that this unique variety of high-Ti basalt may have formed through fractionation of a liquid with the composition of the Apollo 11 orange glass, coupled with assimilation of evolved material (dubbed neuKREEP and having similarities to lunar quartz monzodiorite). Assimilation of this material would impart its REE signature on the liquid, resulting in the elevated REE abundances observed. Minerals such as whitlockite which contain a large portion of the REE budget can be expected to reflect the REE characteristics of the assimilant. To this end, an examination of the whitlockite present in the Apollo 11 Group A basalts was undertaken to search for evidence of the neuKREEP material assimilated.

Jerde, Eric A.↗

Optimal gains for a single polar orbiting satellite

Gains are the spatial weighting of an observation in its neighborhood versus the local values of a model prediction. They are the key to data assimilation, as they are the direct measure of how the data are used to guide the model. As derived in the broad context of data assimilation by Kalman and in the context of meteorology, for example, by Rutherford, the optimal gains are functions of the prediction error covariances between the observation and analysis points. Kalman introduced a very powerful technique that allows one to calculate these optimal gains at the time of each observation. Unfortunately, this technique is both computationally expensive and often numerically unstable for dynamical systems of the magnitude of meteorological models, and thus is unsuited for use in PMIRR data assimilation. However, the optimal gains as calculated by a Kalman filter do reach a steady state for regular observing patterns like that of a satellite. In this steady state, the gains are constants in time, and thus could conceivably be computed off-line. These steady-state Kalman gains (i.e., Wiener gains) would yield optimal performance without the computational burden of true Kalman filtering. We proposed to use this type of constant-in-time Wiener gain for the assimilation of data from PMIRR and Mars Observer.

Banfield, Don↗

Variations in Atmospheric Water Vapor as Seen in Satellite Data and Model Reanalysis Fields

New data sets used to describe components of the hydrologic and energy cycles of the Earth system are currently being produced and disseminated through the NASA EOS DAACs and several data assimilation centers supported by such organizations as NASA and NOAA. These data sets incorporate hybrid data analysis schemes and portray satellite and radiosonde data combined in a diagnostic sense and in a forecast mode. There is a need to develop a better understanding of the accuracy and utility of these global, relatively long-term, datasets to describe components of the hydrologic cycle and to understand atmospheric moisture variability and its relation to climatological significant events Much progress has been made in the last ten years in the development of global atmospheric models and analysis of satellite data for global studies. The atmospheric models have improved in their ability to predict both short term and longer term weather events. This has been possible through better understanding of atmosphere dynamics and very rapid advances in computer technology. Over the same period, data assimilation methods have advanced and unconventional data sources such as aircraft and satellite data, drifting buoys, etc., can be assimilated at non-synaptic times. As a result of these improvements, comprehensive global atmospheric fields based on four-dimensional data assimilation methods now provide one of the most reliable methods for studying dynamical and physical behavior in the atmosphere.

Jedlovec, Gary J.↗

The Impact of Anisotropic Error Correlation Modelling

Most data assimilation systems assume isotropic forecast error correlations, but results from two dimensional Kalman Filter experiments indicate that the correlations can be far from isotropic. In this paper we use a simple two-dimensional data assimilation system, which analyses trace chemical species such as ozone, to assess different approaches to modelling the error correlations. We compare assimilation results using isotropic correlations with results obtained using different approaches to modelling anisotropic correlations: first, using correlations based on the concentrations of the trace chemicals, and secondly using an advective correlation model. We show that these relatively cheap ways of modelling anisotropic correlations give objectively better results than using isotropic correlations. We discuss the possible extension of these approaches to a full 3-D meteorological data assimilation system.

Swinbank, R.↗

Improved Stratospheric Temperature Retrievals for Climate Reanalysis

The Data Assimilation Office (DAO) is embarking on plans to generate a twenty year reanalysis data set of climatic atmospheric variables. One of the focus points will be in the evaluation of the dynamics of the stratosphere. The Stratospheric Sounding Unit (SSU), flown as part of the TIROS Operational Vertical Sounder (TOVS), is one of the primary stratospheric temperature sensors flown consistently throughout the reanalysis period. Seven unique sensors made the measurements over time, with individual instrument characteristics that need to be addressed. The stratospheric temperatures being assimilated across satellite platforms will profoundly impact the reanalysis dynamical fields. To attempt to quantify aspects of instrument and retrieval bias we are carefully collecting and analyzing all available information on the sensors, their instrument anomalies, forward model errors and retrieval biases. For the retrieval of stratospheric temperatures, we adapted the minimum variance approach of Jazwinski (1970) and Rodgers (1976) and applied it to the SSU soundings. In our algorithm, the state vector contains an initial guess of temperature from a model six hour forecast provided by the Goddard EOS Data Assimilation System (GEOS/DAS). This is combined with an a priori covariance matrix, a forward model parameterization, and specifications of instrument noise characteristics. A quasi-Newtonian iteration is used to obtain convergence of the retrieved state to the measurement vector. This algorithm also enables us to analyze and address the systematic errors associated with the unique characteristics of the cell pressures on the individual SSU instruments and the resolving power of the instruments to vertical gradients in the stratosphere. The preliminary results of the improved retrievals and their assimilation as well as baseline calculations of bias and rms error between the NESDIS operational product and col-located ground measurements will be presented.

Rokke, L.↗

Effect of Combined Spaceborne Microwave and Continuous Lightning Measurements on Precipitation Forecasts of the 1998 Ground-Hog Day Storm

We evaluated the impact of several newly available sources of meteorological data on mesoscale model forecasts of precipitation produced by the extra-tropical cyclone that struck Florida on February 2, 1998. Precipitation distributions of convective rainfall events were derived from Special Sensor Microwave Imager (SSM/I) and Multi-Channel Passive Microwave Sensor (TMI) microwave radiometric data by means of the Goddard PROFiling (GPROF) algorithm. Continuous lightning distributions were obtained from sferics measurements obtained from a network of VLF radio receivers. Histograms of coincident sferics frequency distributions were matched to those of precipitation to derive bogus convective rainfall rates from the continuously available sferics measurements. SSM/I and TMI microwave data were used to derive Integrated Precipitable Water (IPW) distributions. The TMI also provided sea surface temperatures (SSTS) of the Loop Current and Gulf Stream with improved structural detail. A series of experiments assimilated IPW and latent heating from the bogus convective rainfall for six-hours in the MM5 mesoscale forecast model to produce nine-hour forecasts of all rainfall as well as other weather parameters. Although continuously assimilating latent heating only slightly improved the surface pressure distribution forecast, it significantly improved the precipitation forecasts. Correctly locating convective rainfall was found critical for assimilating latent heating in the forecast model, but measurement of the rainfall intensity proved to be less important. The improved SSTs also had a positive impact on rainfall forecasts for this case. Assimilating bogus rainfall in the model produced nine-hour forecasts of radar reflectivity distributions that agreed well with coincident observations from the TRMM spaceborne precipitation radar, ground based radar and spaceborne microwave measurements.

Weinman, James A.↗

Adaptive Error Estimation in Linearized Ocean General Circulation Models

Data assimilation methods are routinely used in oceanography. The statistics of the model and measurement errors need to be specified a priori. This study addresses the problem of estimating model and measurement error statistics from observations. We start by testing innovation based methods of adaptive error estimation with low-dimensional models in the North Pacific (5-60 deg N, 132-252 deg E) to TOPEX/POSEIDON (TIP) sea level anomaly data, acoustic tomography data from the ATOC project, and the MIT General Circulation Model (GCM). A reduced state linear model that describes large scale internal (baroclinic) error dynamics is used. The methods are shown to be sensitive to the initial guess for the error statistics and the type of observations. A new off-line approach is developed, the covariance matching approach (CMA), where covariance matrices of model-data residuals are "matched" to their theoretical expectations using familiar least squares methods. This method uses observations directly instead of the innovations sequence and is shown to be related to the MT method and the method of Fu et al. (1993). Twin experiments using the same linearized MIT GCM suggest that altimetric data are ill-suited to the estimation of internal GCM errors, but that such estimates can in theory be obtained using acoustic data. The CMA is then applied to T/P sea level anomaly data and a linearization of a global GFDL GCM which uses two vertical modes. We show that the CMA method can be used with a global model and a global data set, and that the estimates of the error statistics are robust. We show that the fraction of the GCM-T/P residual variance explained by the model error is larger than that derived in Fukumori et al.(1999) with the method of Fu et al.(1993). Most of the model error is explained by the barotropic mode. However, we find that impact of the change in the error statistics on the data assimilation estimates is very small. This is explained by the large representation error, i.e. the dominance of the mesoscale eddies in the T/P signal, which are not part of the 21 by 1" GCM. Therefore, the impact of the observations on the assimilation is very small even after the adjustment of the error statistics. This work demonstrates that simult&neous estimation of the model and measurement error statistics for data assimilation with global ocean data sets and linearized GCMs is possible. However, the error covariance estimation problem is in general highly underdetermined, much more so than the state estimation problem. In other words there exist a very large number of statistical models that can be made consistent with the available data. Therefore, methods for obtaining quantitative error estimates, powerful though they may be, cannot replace physical insight. Used in the right context, as a tool for guiding the choice of a small number of model error parameters, covariance matching can be a useful addition to the repertory of tools available to oceanographers.

Chechelnitsky, Michael Y.↗

A Systematic Error Correction Method for TOVS Radiances

Treatment of systematic errors is crucial for the successful use of satellite data in a data assimilation system. Systematic errors in TOVS radiance measurements and radiative transfer calculations can be as large or larger than random instrument errors. The usual assumption in data assimilation is that observational errors are unbiased. If biases are not effectively removed prior to assimilation, the impact of satellite data will be lessened and can even be detrimental. Treatment of systematic errors is important for short-term forecast skill as well as the creation of climate data sets. A systematic error correction algorithm has been developed as part of a 1D radiance assimilation. This scheme corrects for spectroscopic errors, errors in the instrument response function, and other biases in the forward radiance calculation for TOVS. Such algorithms are often referred to as tuning of the radiances. The scheme is able to account for the complex, air-mass dependent biases that are seen in the differences between TOVS radiance observations and forward model calculations. We will show results of systematic error correction applied to the NOAA 15 Advanced TOVS as well as its predecessors. We will also discuss the ramifications of inter-instrument bias with a focus on stratospheric measurements.

Joiner, Joanna↗

Comparison of the Variability of Constituents from a 3D Simulation with Observations in the Tropics

Fields from chemical transport models (CTMs) driven by assimilated winds have been shown to represent many aspects of observations at middle and high latitudes. Comparisons of modeled fields with observations have been much poorer in the tropics. Noise in the tropical winds produces variability in constituents in the tropics that greatly exceeds what is observed, Constituent gradients between the tropics, subtropics and middle latitudes are weaker than observed in the lower stratosphere. Tropical transport is shown to be much more realistic using a CTM driven by winds from the assimilation system that has recently been developed by the Data Assimilation Office at Goddard Space Flight Center. The Finite Volume Community Climate Model (FV-CCM), a general circulation model model that uses the NCAR CCM physics and the Lin and Rood dynamical core, is at the core of the new assimilation system. Realistic transport in the tropics, particularly troposphere to stratosphere exchange, is an obvious requirement for realistic representation of the stratosphere to troposphere transport of ozone and other stratospheric constituents to the middle latitude upper troposphere.

Douglass, A. R.↗

Mass Conservation in a Chemical Transport Model and its Effect on CO2 and SF6 Simulations

Chemical transport models (CTMs) must conserve mass to be useful for applications involving assessment of the effect of various pollutants on the troposphere and stratosphere. Furthermore, calculations of the evolution of constituents such as SF6 are used to evaluate overall model transport, and interpretation of such simulations is clouded if mass conservation is not assured. For realistic simulations or predictions, it is crucial that constituents are not produced or lost by transport or other processes in the CTMs. Analysis of CO2 and SF6 experiments using a CTM shows that problems with mass conservation can seriously degrade the simulations. Failure to conserve mass results from inconsistency of the surface pressure tendency and the divergence of horizontal mass flux when the model is forced by assimilated meteorological data. We have developed an effective method to eliminate the inconsistency by modifying the divergent part of the wind field. The changes in the wind fields are quite small but the impact on mass conservation is large. Parameterizations of physical processes such as convection or turbulent transport can also affect mass conservation. The lack of conservation is small but accumulates when integrations are lengthy such as required for SF6. This lack of conservation is found using winds from either a GCM or from an assimilation system. A simple adjustment removes much of the inaccuracy in the convective parameterization. A CO2 simulation using assimilated winds from the most recent version of the Goddard Earth Observing System Data Assimilation System will be used to illustrate the impact of these transport improvements.

Zhu, Z.↗

Investigation and Development of Data-Driven D-Region Model for HF Systems Impacts

Space Environment Corporation (SEC) and RP Consultants (RPC) are to develop and validate a weather-capable D region model for making High Frequency (HF) absorption predictions in support of the HF communications and radar communities. The weather-capable model will assimilate solar and earth space observations from NASA satellites. The model will account for solar-induced impacts on HF absorption, including X-rays, Solar Proton Events (SPE's), and auroral precipitation. The work plan includes: I . Optimize D-region model to quickly obtain ion and electron densities for proper HF absorption calculations. 2. Develop indices-driven modules for D-region ionization sources for low, mid, & high latitudes including X-rays, cosmic rays, auroral precipitation, & solar protons. (Note: solar spectrum & auroral modules already exist). 3. Setup low-cost monitors of existing HF beacons and add one single-frequency beacon. 4. Use PENEX HF-link database with HF monitor data to validate D-region/HF absorption model using climatological ionization drivers. 5. Develop algorithms to assimilate NASA satellite data of solar, interplanetary, and auroral observations into ionization source modules. 6. Use PENEX HF-link & HF-beacon data for skill score comparison of assimilation versus climatological D-region/HF absorption model. Only some satellites are available for the PENEX time period, thus, HF-beacon data is necessary. 7. Use HF beacon monitors to develop HF-link data assimilation algorithms for regional improvement to the D-region/HF absorption model.

Eccles, J. V.↗

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↗