Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Variational data 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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Evaluating the Global Water and Energy Cycles in Reanalyses Using Ensemble Spread

Reanalyses provide continuous observation based global data over weather and climate time scales through observational assimilation. Ensemble forecasts are a component of modern data assimilation systems, providing information to the analysis process. Storing this spread over the duration of a climate reanalysis will provide insight to the variations of the model that underlies the reanalysis data. Global water and energy cycles (WECs) are a critical aspect of both weather and climate, yet there are some large gaps in the observation of key components as well as their representation in models and reanalyses. NASA’s Water and Energy cycle Studies (NEWS) program is developing a corrected and balanced budget data set that is closed at monthly and regional scales during the 21st century, so far. In this paper, we will evaluate reanalysis WEC components against the initial NEWS collection of observation, but integrate the analysis spread from the reanalysis data assimilation as a new perspective on the evaluation. The initial evaluation will cover a few more recent years using ERA5 10 member spread. As the GMAO’s production of the Modern-Era Retrospective analysis for Research and Applications 21st Century (MERRA-21C) progresses, the effort will include the MERRA-21C 32 member spread. Ultimately, we will explore the potential of the ensemble spread as a measure of uncertainty in the reanalysis.

Michael Bosilovich↗
NASA NTRS1987

Assimilation of scalar versus horizontal gradient information from the VAS into a mesoscale model

Comparisons are made between analyses and forecasts which incorporate VAS geopotential data as either scalar or horizontal gradient information for a case study on the AVE/VAS day of Mar. 6, 1982. On this day, incorporating the VAS information in analysis as a variational constraint on horizontal geopotential gradients significantly mitigated the effects of large data biases which made VAS assimilation by standard scalar methods very difficult. A subsequent forecast made from the gradient assimilation was superior to one made from the standard analysis and of comparable quality in geopotentials to a control forecast from synoptic data. Most impact was noted in the forecasts of vertical motion and precipitation in the gradient vs this control simulation.

Diak, George↗
NASA NTRS1992

GEWEX

The GEWEX (Global Energy and Water cycle EXperiment) objectives are as follows: to determine the hydrologic cycle and energy fluxes by means of global measurements of observable atmospheric and surface properties; to model the global hydrologic cycle and its effects on the atmosphere and oceans; to foster the development of observing techniques and data management and assimilation systems; and to develop the ability to predict the variations of global and regional hydrologic processes and water resources, and their response to environmental change.

Source record↗
NASA NTRS2008

Assimilation of Atmospheric InfraRed Sounder (AIRS) Profiles using WRF-Var

The Weather Research and Forecasting (WRF) model contains a three-dimensional variational (3DVAR) assimilation system (WRF-Var), which allows a user to join data from multiple sources into one coherent analysis. WRF-Var combines observations with a background field traditionally generated using a previous model forecast through minimization of a cost function. In data sparse regions, remotely-sensed observations may be able to improve analyses and produce improved forecasts. One such source comes from the Atmospheric Infrared Sounder (AIRS), which together with the Advanced Microwave Sounding Unit (AMSU), represents one of the most advanced space-based atmospheric sounding systems. The combined AIRS/AMSU system provides radiance measurements used as input to a sophisticated retrieval scheme which has been shown to produce temperature profiles with an accuracy of 1 K over 1 km layers and humidity profiles with accuracy of 15% in 2 km layers in both clear and partly cloudy conditions. The retrieval algorithm also provides estimates of the accuracy of the retrieved values at each pressure level, allowing the user to select profiles based on the required error tolerances of the application. The purpose of this paper is to describe a procedure to optimally assimilate high-resolution AIRS profile data into a regional configuration of the Advanced Research WRF (ARW) version 2.2 using WRF-Var. The paper focuses on development of background error covariances for the regional domain and background field type using gen_be and an optimal methodology for ingesting AIRS temperature and moisture profiles as separate overland and overwater retrievals with different error characteristics in the WRF-Var. The AIRS thermodynamic profiles are obtained from the version 5.0 Earth Observing System (EOS) science team retrieval algorithm and contain information about the quality of each temperature layer. The quality indicators are used to select the highest quality temperature and moisture data for each profile location and pressure level. Analyses are run to produce quasi-real-time regional weather forecasts over the continental U.S. The preliminary assessment of the impact of the AIRS profiles will focus on intelligent use of the quality indicators, optimized tuning of the WRF-Var, and comparison of analysis soundings to radiosondes.

Zavodsky, Brad↗
NASA NTRS2003

Lightning and Other Influences On Tropical Tropospheric Ozone: Empirical Studies of Covariation

Tropical and subtropical tropospheric ozone are important radiatively active species, with particularly large effects in the upper third of the troposphere. Temporal variability of O3 has proved difficult to simulate day by day in process models. Thus, individual roles of lightning, biomass burning, and other pollution in providing precursor NO(x), radicals, and chain carriers (CO, hydrocarbons) remain unquantified by simulation, and it is theoretically reasonable that individual roles are magnified by a joint synergy. We use wavelet analysis and Burg-algorithm maximum entropy spectral analyses to describe time-scales and correlation of ozone with proxies for processes controlling its concentration. Our empirical studies link time variations apparent in several datasets: the SHADOZ (Southern Hemisphere Additional Ozonesondes) network stations (Nairobi, Fiji), and auxiliary series with power to explain ozone-determining processes, with some interpretation based on the TTO (Tropical Tropospheric Ozone) product derived from TOMS (the Total Ozone Mapping Spectrometer). The auxiliary series are The OTD/LIS(Optical Transient Detector/Lightning Imaging Sensor) measurements of the lightning NO(x) source, the OLR (Outgoing Longwave Radiation)measurement of high-topped clouds, and standard meteorological variables from the United States NCEP (National Centers for Environmental Prediction) and Data Assimilation Office analyses. Concentrating on equatorial ozone, we compare the statistical evidence on the variability of tropospheric ozone. Important variations occur on approximately two-week, two-month (Madden-Julian Oscillation) and annual scales, and relations with OLR suggest controls associated with continental clouds. Hence we are now using the Lightning Imaging Sensor data set to indicate NO(x) sources. We report initial results defining relative roles of the process mentioned affecting O3 using their covariance properties.

Chatfield, Robert B.↗
NASA NTRS1987

A VAS-numerical model impact study using the Gal-Chen variational approach

Numerical experiments to assess the impact of incorporating temperature data from the VISSR Atmospheric Sounder (VAS) using the assimilation technique developed by Gal-Chen (1986) modified for use in the Mesoscale Atmospheric Simulation System (MASS) model were conducted. The scheme is designed to utilize the high temporal and horizontal resolution of satellite retrievals while maintaining the fine vertical structure generated by the model. This is accomplished by adjusting the model lapse rates to reflect thicknesses retrieved from VAS and applying a three-dimensional variational that preserves the distribution of the geopotential fields in the model. A nudging technique whereby the model temperature fields are gradually adjusted toward the updated temperature fields during model integration is also tested. An adiabatic version of MASS is used in all experiments to better isolate mass-momentum imbalances. The method has a sustained impact over an 18 hr model simulation.

Aune, Robert M.↗
NASA NTRS2003

Prospect of Using Numerical Dynamo Model for Prediction of Geomagnetic Secular Variation

Modeling of the Earth's core has reached a level of maturity to where the incorporation of observations into the simulations through data assimilation has become feasible. Data assimilation is a method by which observations of a system are combined with a model output (or forecast) to obtain a best guess of the state of the system, called the analysis. The analysis is then used as an initial condition for the next forecast. By doing assimilation, not only we shall be able to predict partially secular variation of the core field, we could also use observations to further our understanding of dynamical states in the Earth's core. One of the first steps in the development of an assimilation system is a comparison between the observations and the model solution. The highly turbulent nature of core dynamics, along with the absence of any regular external forcing and constraint (which occurs in atmospheric dynamics, for example) means that short time comparisons (approx. 1000 years) cannot be made between model and observations. In order to make sensible comparisons, a direct insertion assimilation method has been implemented. In this approach, magnetic field observations at the Earth's surface have been substituted into the numerical model, such that the ratio of the multiple components and the dipole component from observation is adjusted at the core-mantle boundary and extended to the interior of the core, while the total magnetic energy remains unchanged. This adjusted magnetic field is then used as the initial field for a new simulation. In this way, a time tugged simulation is created which can then be compared directly with observations. We present numerical solutions with and without data insertion and discuss their implications for the development of a more rigorous assimilation system.

Kuang, Weijia↗
NASA NTRS1991

Marine research using data from TOPEX/Poseidon

The Proudman Oceanographic Laboratory (POL), the Institute of Oceanographic Sciences Deacon Laboratory (IOSDL), and the British National Space Centre (BNSC) plan to use radar altimeter data from the TOPEX/POSEIDON mission in five main areas of research. These are: (1) investigation of variations in Southern Ocean sea levels and circulation; (2) related numerical modeling and data assimilation; (3) global and regional ocean tides; (4) North Sea and northeast Atlantic studies; and (5) winds and waves. Each one of these topics is described briefly below while more information can be extracted from the United Kingdom response to the TOPEX/POSEIDON Announcement of Opportunity.

Woodworth, Philip L.↗
NASA NTRS2001

Establishing an IERS Sub-Center for Ocean Angular Momentum

With the objective of establishing the Special Bureau for the Oceans (SBO), a new archival center for ocean angular momentum (OAM) products, we have computed and analyzed a number of OAM products from several ocean models, with and without data assimilation. All three components of OAM (axial term related to length of day variations and equatorial terms related to polar motion) have been examined in detail, in comparison to the respective Earth rotation parameters. An 11+ year time series of OAM given at 5-day intervals has been made publicly available. Other OAM products spanning longer periods and with higher temporal resolution, as well as products calculated from ocean model/data assimilation systems, have been prepared and should become part of the SBO archives in the near future.

Ponte, Rui M.↗
NASA NTRS2015

RACORO Continental Boundary Layer Cloud Investigations: 1. Case Study Development and Ensemble Large-Scale Forcings

Observation-based modeling case studies of continental boundary layer clouds have been developed to study cloudy boundary layers, aerosol influences upon them, and their representation in cloud- and global-scale models. Three 60 h case study periods span the temporal evolution of cumulus, stratiform, and drizzling boundary layer cloud systems, representing mixed and transitional states rather than idealized or canonical cases. Based on in situ measurements from the Routine AAF (Atmospheric Radiation Measurement (ARM) Aerial Facility) CLOWD (Clouds with Low Optical Water Depth) Optical Radiative Observations (RACORO) field campaign and remote sensing observations, the cases are designed with a modular configuration to simplify use in large-eddy simulations (LES) and single-column models. Aircraft measurements of aerosol number size distribution are fit to lognormal functions for concise representation in models. Values of the aerosol hygroscopicity parameter, kappa, are derived from observations to be approximately 0.10, which are lower than the 0.3 typical over continents and suggestive of a large aerosol organic fraction. Ensemble large-scale forcing data sets are derived from the ARM variational analysis, European Centre for Medium-Range Weather Forecasts, and a multiscale data assimilation system. The forcings are assessed through comparison of measured bulk atmospheric and cloud properties to those computed in "trial" large-eddy simulations, where more efficient run times are enabled through modest reductions in grid resolution and domain size compared to the full-sized LES grid. Simulations capture many of the general features observed, but the state-of-the-art forcings were limited at representing details of cloud onset, and tight gradients and high-resolution transients of importance. Methods for improving the initial conditions and forcings are discussed. The cases developed are available to the general modeling community for studying continental boundary clouds.

Boundary layers↗
NASA NTRS2010

Impact of Atmospheric Infrared Sounder (AIRS) Thermodynamic Profiles on Regional Weather Forecasting

In data sparse regions, remotely-sensed observations can be used to improve analyses and lead to better forecasts. One such source comes from the Atmospheric Infrared Sounder (AIRS), which together with the Advanced Microwave Sounding Unit (AMSU), provides temperature and moisture profiles with accuracy comparable to that of radiosondes. The purpose of this paper is to describe a procedure to assimilate AIRS thermodynamic profile data into a regional configuration of the Advanced Research Weather Research and Forecasting (WRF-ARW) model using its three-dimension variational (3DVAR) analysis component (WRF-Var). Quality indicators are used to select only the highest quality temperature and moisture profiles for assimilation in both clear and partly cloudy regions. Separate error characteristics for land and water profiles are also used in the assimilation process. Assimilation results indicate that AIRS profiles produce an analysis closer to in situ observations than the background field. Forecasts from a 37-day case study period in the winter of 2007 show that AIRS profile data can lead to improvements in 6-h cumulative precipitation forecasts due to instability added in the forecast soundings by the AIRS profiles. Additionally, in a convective heavy rainfall event from February 2007, assimilation of AIRS profiles produces a more unstable boundary layer resulting in enhanced updrafts in the model. These updrafts produce a squall line and precipitation totals that more closely reflect ground-based observations than a no AIRS control forecast. The location of available high-quality AIRS profiles ahead of approaching storm systems is found to be of paramount importance to the amount of impact the observations will have on the resulting forecasts.

Chou, Shih-Hung↗
DOE OSTI2024

Earth System Reanalysis in Support of Climate Model Improvements

Recent climate model developments, established through increased model resolution, have led to substantial improvements in model simulations of the time-evolving, coupled Earth system and its subcomponents. However, regardless of resolution, climate models will always produce climate features and variability that differ from the real world and will be prone to biases. This is due to many remaining uncertainties, such as in parametric and structural model uncertainty, in the initial conditions prescribed, and in the prescribed (scenario) forcing which varies on decadal to centennial timescales. Further model improvements are expected to arise specifically from improved representation of physical processes realized through model-data fusion. This will create an unprecedented opportunity to better exploit a large array of Earth observations, from in situ measurements to weather radars and satellite observations, as the resolved scales of the models approach those of the observations. For this, climate DA will be the central tool to bring models and observations into consistency, by improving initial conditions, inferring uncertain model parameters and structure, and quantifying uncertainty. Generally, there will be advantages and complementarities of adjoint-based smoother approaches, ensemble-based filter approaches, or new ML-inspired approaches. Yet, the ever-increasing model resolution will present growing challenges arising from computational cost, calling for new ways of performing data assimilation and model optimization. Using the complementarity in a hybrid approach, blending tools and concepts from variational, ensemble and ML methods might be what is required in the future. In this context ML could be important to handle non-linear responses, and to better approximate non-Gaussian distributions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Development and evaluation of a new 4DEnVar-based weakly coupled ocean data assimilation system in E3SMv2

The development, implementation, and evaluation of a new weakly coupled ocean data assimilation (WCODA) system for the fully coupled Energy Exascale Earth System Model version 2 (E3SMv2) utilizing the four-dimensional ensemble variational (4DEnVar) method are presented in this study. The 4DEnVar method, based on the dimension-reduced projection four-dimensional variational (DRP-4DVar) approach, replaces the adjoint model with the ensemble technique, thereby reducing computational demands. Monthly mean ocean temperature and salinity data from the EN4.2.1 reanalysis are integrated into the ocean component of E3SMv2 from 1950 to 2021 with the goal of providing realistic initial conditions for decadal predictions and predictability studies. The performance of the WCODA system is assessed using various metrics, including the reduction rate of the cost function, root mean square error (RMSE) differences, correlation differences, and model biases. Results indicate that the WCODA system effectively assimilates the reanalysis data into the climate model, consistently achieving negative reduction rates of the cost function and notable improvements in RMSE and correlation across various ocean layers and regions. Significant enhancements are observed in the upper ocean layers across the majority of global ocean regions, particularly in the north Atlantic, north Pacific, and Indian Ocean. Model biases in sea surface temperature and salinity are also substantially reduced. For sea surface temperature, cold biases in the north Pacific and north Atlantic are diminished by about 1–2 °C, and warm biases in the Southern Ocean are corrected by approximately 1.5–2.5 °C. In terms of salinity, improvements are observed with bias reductions of about 0.5–1 psu in the north Atlantic and north Pacific and up to 1.5 psu in parts of the Southern Ocean. The ultimate goal of the WCODA system is to advance the predictive capabilities of E3SM for subseasonal to decadal climate predictions, thereby supporting research on strategic energy-sector policies and planning.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2002

Assimilation Experiments of One-Dimensional Variational Analyses with GPS/MET Refractivity

The GPS/MET 1995 experiment demonstrated the feasibility of remote soundings of the Earth's atmosphere using the GPS radio occultation technique. Subsequent studies assessed the information content of such measurements. Several groups have attempted to assimilate the GPS radio occultation measurements into a global Data Assimilation System (DAS). These attempts followed in general the methodology proposed by Eyre, and range from die act four-dimensional variational assimilation of bending angles to assimilation of inverted profiles of temperature and/or humidity. We present here an hybrid approach. We attempted to take advantage of the accuracy of current, weather models to constrain one-dimensional variational analyses using GPS refractivity. Retrieved profiles of atmospheric parameters were then assimilated like other types of observations in a global DAS for issuing the next forecast period.

Poli, Paul↗
DOE OSTI2020

Simulation of Continental Shallow Cumulus Populations Using an Observation-Constrained Cloud-System Resolving Model

Continental shallow cumulus (ShCu) clouds observed on 30 August 2016 during the Holistic Interactions of Shallow Clouds, Aerosols, and Land-Ecosystems (HI-SCALE) field campaign are simulated by using an observation-constrained cloud-system resolving model. On this day, ShCu forms over Oklahoma and southern Kansas and some of these clouds transition to deeper, precipitating convection during the afternoon. We apply a four-dimensional ensemble-variational (4DEnVar) hybrid technique in the Community Gridpoint Statistical Interpolation (GSI) system to assimilate operational data sets and unique boundary layer measurements including a Raman lidar, radar wind profilers, radiosondes, and surface stations collected by the U.S. Department of Energy's (DOE) Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) atmospheric observatory into the Weather Research and Forecasting (WRF) model to ascertain how improved environmental conditions can influence forecasts of ShCu populations and the transition to deeper convection. Independent observations from aircraft, satellite, as well as ARM's remote sensors are used to evaluate model performance in different aspects. Several model experiments are conducted to identify the impact of data assimilation (DA) on the prediction of clouds evolution. The analyses indicate that ShCu populations are more accurately reproduced after DA in terms of cloud initiation time and cloud base height, which can be attributed to an improved representation of the ambient meteorological conditions and the convective boundary layer. Extending the assimilation to 18 UTC (local noon) also improved the simulation of shallow-to-deep transitions of convective clouds.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Large Ensemble Particle Filter for Spatial Climate Reconstructions Using a Linear Inverse Model

Proxy records that document the last 2000 years of climate provide evidence for the wide range of the natural climate variability from inter-annual to secular timescales not captured by the short window of recent direct observations. Assessing climate models ability to reproduce such natural variations is crucial to understand climate sensitivity and impacts of future climate change. Paleoclimate data assimilation (PDA) offers a powerful way to extend the short instrumental period by optimally combining the physics described by General Circulation Climate Models (GCMs) with information from available proxy records while taking into account their uncertainties. Here we present a new PDA approach based on a sequential importance resampling (SIR) Particle filter (PF) that uses Linear Inverse Modeling (LIM) as an emulator of several CMIP-class GCMs. We examine in a perfect-model framework the skill of the various LIMs to forecast the dynamics of the surface temperatures and provide spatial field reconstructions over the last millennium in a SIR PF. Our results show that the LIMs allow for skillful ensemble forecasts at 1-year lead-time based on GCMs dynamical knowledge with best prediction in the tropics and the North Atlantic. The PDA further provides a set of physically consistent spatial fields allowing robust uncertainty quantification related to climate models biases and proxy spatial sampling. Our results indicate that the LIM yields dynamical memory improving climate variability reconstructions and support the use of the LIM as a GCM-emulator in real reconstruction to propagate large ensembles of particles at low cost in SIR PF.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2001

Temperature Humidity and Sea Level Pressure Increments Induced by 1DVAR Analysis of GPS Refractivity

The Global Positioning System (GPS) transmitted signals are affected by the atmosphere. Using the radio occultation technique, where a receiver is placed on a low-Earth-orbiting platform. it is possible to perform soundings. by measuring the bending angles of the rays. The information can be converted into atmospheric refractivity. We have developed a one dimensional variational (1DVAR) analysis that uses GPS/MET 1995 refractivity and 6-hour FVDAS (Finite Volume Data Assimilation System) forecasts as background information to constrain the retrievals. The analysis increments are defined as 1DVAR minus background temperature, humidity and sea level pressure. Before assimilating the 1DVAR profiles into the FVDAS. the increments need to be understood. First, some bias could be induced in the retrievals when confronted with actual biased data: second. bias in the back-round could create undesired bias in the retrievals. Anv bias in the analyses will ultimately change the climatology of the model the retrievals will be assimilated into. We relate the increments to the reduction of the difference between observed minus computed refractivity profiles. We also point out the difference in the mean increments using backgrounds which have assimilated either NESDIS TIROS Operational Vertical Sounder (TOVS) operational retrievals or Data Assimilation Office (DAO) TOVS interactive retrievals. The climatology of the model in terms of refractivity is significantly different and this impacts the GPS 1DVAR increments. This testifies that changing the basic load of assimilated data has an influence on the impact the GPS data may have in a DAS.

Poli, Paul↗