Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ensemble 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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS2012

Inferring Land Surface Model Parameters for the Assimilation of Satellite-Based L-Band Brightness Temperature Observations into a Soil Moisture Analysis System

The Soil Moisture and Ocean Salinity (SMOS) satellite mission provides global measurements of L-band brightness temperatures at horizontal and vertical polarization and a variety of incidence angles that are sensitive to moisture and temperature conditions in the top few centimeters of the soil. These L-band observations can therefore be assimilated into a land surface model to obtain surface and root zone soil moisture estimates. As part of the observation operator, such an assimilation system requires a radiative transfer model (RTM) that converts geophysical fields (including soil moisture and soil temperature) into modeled L-band brightness temperatures. At the global scale, the RTM parameters and the climatological soil moisture conditions are still poorly known. Using look-up tables from the literature to estimate the RTM parameters usually results in modeled L-band brightness temperatures that are strongly biased against the SMOS observations, with biases varying regionally and seasonally. Such biases must be addressed within the land data assimilation system. In this presentation, the estimation of the RTM parameters is discussed for the NASA GEOS-5 land data assimilation system, which is based on the ensemble Kalman filter (EnKF) and the Catchment land surface model. In the GEOS-5 land data assimilation system, soil moisture and brightness temperature biases are addressed in three stages. First, the global soil properties and soil hydraulic parameters that are used in the Catchment model were revised to minimize the bias in the modeled soil moisture, as verified against available in situ soil moisture measurements. Second, key parameters of the "tau-omega" RTM were calibrated prior to data assimilation using an objective function that minimizes the climatological differences between the modeled L-band brightness temperatures and the corresponding SMOS observations. Calibrated parameters include soil roughness parameters, vegetation structure parameters, and the single scattering albedo. After this climatological calibration, the modeling system can provide L-band brightness temperatures with a global mean absolute bias of less than 10K against SMOS observations, across multiple incidence angles and for horizontal and vertical polarization. Third, seasonal and regional variations in the residual biases are addressed by estimating the vegetation optical depth through state augmentation during the assimilation of the L-band brightness temperatures. This strategy, tested here with SMOS data, is part of the baseline approach for the Level 4 Surface and Root Zone Soil Moisture data product from the planned Soil Moisture Active Passive (SMAP) satellite mission.

Reichle, Rolf H.↗
NASA NTRS2016

SMAP Data Assimilation at the GMAO

The NASA Soil Moisture Active Passive (SMAP) mission has been providing L-band (1.4 GHz) passive microwave brightness temperature (Tb) observations since April 2015. These observations are sensitive to surface(0-5 cm) soil moisture. Several of the key applications targeted by SMAP, however, require knowledge of deeper-layer, root zone (0-100 cm) soil moisture, which is not directly measured by SMAP. The NASA Global Modeling and Assimilation Office (GMAO) contributes to SMAP by providing Level 4 data, including the Level 4 Surface and Root Zone Soil Moisture(L4_SM) product, which is based on the assimilation of SMAP Tb observations in the ensemble-based NASA GEOS-5 land surface data assimilation system. The L4_SM product offers global data every three hours at 9 km resolution, thereby interpolating and extrapolating the coarser- scale (40 km) SMAP observations in time and in space (both horizontally and vertically). Since October 31, 2015, beta-version L4_SM data have been available to the public from the National Snow and Ice Data Center for the period March 31, 2015, to near present, with a mean latency of approx. 2.5 days.

Land Surface Model↗
NASA NTRS2010

Evaluating the Utility of Remotely-Sensed Soil Moisture Retrievals for Operational Agricultural Drought Monitoring

Soil moisture is a fundamental data source used by the United States Department of Agriculture (USDA) International Production Assessment Division (IPAD) to monitor crop growth stage and condition and subsequently, globally forecast agricultural yields. Currently, the USDA IPAD estimates surface and root-zone soil moisture using a two-layer modified Palmer soil moisture model forced by global precipitation and temperature measurements. However, this approach suffers from well-known errors arising from uncertainty in model forcing data and highly simplified model physics. Here we attempt to correct for these errors by designing and applying an Ensemble Kalman filter (EnKF) data assimilation system to integrate surface soil moisture retrievals from the NASA Advanced Microwave Scanning Radiometer (AMSR-E) into the USDA modified Palmer soil moisture model. An assessment of soil moisture analysis products produced from this assimilation has been completed for a five-year (2002 to 2007) period over the North American continent between 23degN - 50degN and 128degW - 65degW. In particular, a data denial experimental approach is utilized to isolate the added utility of integrating remotely-sensed soil moisture by comparing EnKF soil moisture results obtained using (relatively) low-quality precipitation products obtained from real-time satellite imagery to baseline Palmer model runs forced with higher quality rainfall. An analysis of root-zone anomalies for each model simulation suggests that the assimilation of AMSR-E surface soil moisture retrievals can add significant value to USDA root-zone predictions derived from real-time satellite precipitation products.

Bolten, John D.↗
NASA NTRS1999

The Behavior of Filters and Smoothers for Strongly Nonlinear Dynamics

The Kalman filter is the optimal filter in the presence of known Gaussian error statistics and linear dynamics. Filter extension to nonlinear dynamics is non trivial in the sense of appropriately representing high order moments of the statistics. Monte Carlo, ensemble-based, methods have been advocated as the methodology for representing high order moments without any questionable closure assumptions (e.g., Miller 1994). Investigation along these lines has been conducted for highly idealized dynamics such as the strongly nonlinear Lorenz (1963) model as well as more realistic models of the oceans (Evensen and van Leeuwen 1996) and atmosphere (Houtekamer and Mitchell 1998). A few relevant issues in this context are related to the necessary number of ensemble members to properly represent the error statistics and, the necessary modifications in the usual filter equations to allow for correct update of the ensemble members (Burgers 1998). The ensemble technique has also been applied to the problem of smoothing for which similar questions apply. Ensemble smoother examples, however, seem to quite puzzling in that results of state estimate are worse than for their filter analogue (Evensen 1997). In this study, we use concepts in probability theory to revisit the ensemble methodology for filtering and smoothing in data assimilation. We use Lorenz (1963) model to test and compare the behavior of a variety implementations of ensemble filters. We also implement ensemble smoothers that are able to perform better than their filter counterparts. A discussion of feasibility of these techniques to large data assimilation problems will be given at the time of the conference.

Zhu, Yanqiu↗
NASA NTRS1999

Behavior of Filters and Smoothers for Strongly Nonlinear Dynamics

The Kalman filter is the optimal filter in the presence of known gaussian error statistics and linear dynamics. Filter extension to nonlinear dynamics is non trivial in the sense of appropriately representing high order moments of the statistics. Monte Carlo, ensemble-based, methods have been advocated as the methodology for representing high order moments without any questionable closure assumptions. Investigation along these lines has been conducted for highly idealized dynamics such as the strongly nonlinear Lorenz model as well as more realistic models of the means and atmosphere. A few relevant issues in this context are related to the necessary number of ensemble members to properly represent the error statistics and, the necessary modifications in the usual filter situations to allow for correct update of the ensemble members. The ensemble technique has also been applied to the problem of smoothing for which similar questions apply. Ensemble smoother examples, however, seem to be quite puzzling in that results state estimates are worse than for their filter analogue. In this study, we use concepts in probability theory to revisit the ensemble methodology for filtering and smoothing in data assimilation. We use the Lorenz model to test and compare the behavior of a variety of implementations of ensemble filters. We also implement ensemble smoothers that are able to perform better than their filter counterparts. A discussion of feasibility of these techniques to large data assimilation problems will be given at the time of the conference.

Zhu, Yanqui↗
NASA NTRS2016

Reducing Multisensor Satellite Monthly Mean Aerosol Optical Depth Uncertainty: 1. Objective Assessment of Current AERONET Locations

Various space-based sensors have been designed and corresponding algorithms developed to retrieve aerosol optical depth (AOD), the very basic aerosol optical property, yet considerable disagreement still exists across these different satellite data sets. Surface-based observations aim to provide ground truth for validating satellite data; hence, their deployment locations should preferably contain as much spatial information as possible, i.e., high spatial representativeness. Using a novel Ensemble Kalman Filter (EnKF)- based approach, we objectively evaluate the spatial representativeness of current Aerosol Robotic Network (AERONET) sites. Multisensor monthly mean AOD data sets from Moderate Resolution Imaging Spectroradiometer, Multiangle Imaging Spectroradiometer, Sea-viewing Wide Field-of-view Sensor, Ozone Monitoring Instrument, and Polarization and Anisotropy of Reflectances for Atmospheric Sciences coupled with Observations from a Lidar are combined into a 605-member ensemble, and AERONET data are considered as the observations to be assimilated into this ensemble using the EnKF. The assessment is made by comparing the analysis error variance (that has been constrained by ground-based measurements), with the background error variance (based on satellite data alone). Results show that the total uncertainty is reduced by approximately 27% on average and could reach above 50% over certain places. The uncertainty reduction pattern also has distinct seasonal patterns, corresponding to the spatial distribution of seasonally varying aerosol types, such as dust in the spring for Northern Hemisphere and biomass burning in the fall for Southern Hemisphere. Dust and biomass burning sites have the highest spatial representativeness, rural and oceanic sites can also represent moderate spatial information, whereas the representativeness of urban sites is relatively localized. A spatial score ranging from 1 to 3 is assigned to each AERONET site based on the uncertainty reduction, indicating its representativeness level.

aerosol optical depth↗
NASA NTRS2014

Assimilation of SMOS Retrieved Soil Moisture into the Land Information System

Soil moisture retrievals from the Soil Moisture and Ocean Salinity (SMOS) instrument are assimilated into the Noah land surface model (LSM) within the NASA Land Information System (LIS). Before assimilation, SMOS retrievals are bias-corrected to match the model climatological distribution using a Cumulative Distribution Function (CDF) matching approach. Data assimilation is done via the Ensemble Kalman Filter. The goal is to improve the representation of soil moisture within the LSM, and ultimately to improve numerical weather forecasts through better land surface initialization. We present a case study showing a large area of irrigation in the lower Mississippi River Valley, in an area with extensive rice agriculture. High soil moisture value in this region are observed by SMOS, but not captured in the forcing data. After assimilation, the model fields reflect the observed geographic patterns of soil moisture. Plans for a modeling experiment and operational use of the data are given. This work helps prepare for the assimilation of Soil Moisture Active/Passive (SMAP) retrievals in the near future.

Data Assimilation↗
NASA NTRS2018

Synergy of Observations and Dynamo Models to Understand and Predict Solar Activity Cycles

The long-standing problem of understanding the evolution of the global magnetic fields that drive solar activity through different temporal scales is becoming more tractable because, in addition to 400 years of sunspot records, we now have almost 4 solar cycles of magnetic field observations. These observations allow us to discern physical connections between dynamo model variables and observations using data assimilation analysis. In particular, the Ensemble Kalman Filter approach takes into account uncertainties in both observations and modeling and allows us to make reliable forecasts of solar cycle activity by using a relatively simple non-linear dynamical model of the solar dynamo. To expand this approach for more complex 2D and 3D dynamo modeling, it is necessary to decompose the observed synoptic magnetograms into poloidal and toroidal field components. In this presentation I will present initial results on magnetogram decomposition and assimilation of magnetogram data into dynamo modeling.

Kitiashvili, Irina N.↗
NASA NTRS2018

Synergy of Observations and Dynamo Models to Understand and Predict Solar Activity Cycles

The long-standing problem of understanding the evolution of the global magnetic fields that drive solar activity through different temporal scales is becoming more tractable because, in addition to 400 years of sunspot records, we now have almost 4 solar cycles of magnetic field observations. These observations allow us to discern physical connections between dynamo model variables and observations using data assimilation analysis. In particular, the Ensemble Kalman Filter approach takes into account uncertainties in both observations and modeling and allows us to make reliable forecasts of solar cycle activity cycles by using a relatively simple non-linear dynamical model of the solar dynamo. To expand this approach for more complex 2D and 3D dynamo modeling, it is necessary to decompose the observed synoptic magnetograms into poloidal and toroidal field components. In this presentation I will present initial results on magnetogram decomposition and assimilation of magnetogram data into dynamo modeling.

Synergy↗
NASA NTRS2019

Application of Synoptic Magnetograms for Prediction of Solar Activity Using Ensemble Kalman Filter

Solar activity predictions using the data assimilation approach have demonstrated great potential to build reliable long-term forecasts of solar activity. In particular, it has been shown that the Ensemble Kalman Filter (EnKF) method applied to a non-linear dynamo model is capable of predicting solar activity up to one sunspot cycle ahead in time, as well as estimating the properties of the next cycle a few years before it begins. These developments assume an empirical relationship between the mean toroidal magnetic field flux and the sunspot number. Estimated from the sunspot number series, variations of the toroidal field have been used to assimilate the data into the Parker-Kleeorin-Ruzmakin (PKR) dynamo model by applying the EnKF method. The dynamo model describes the evolution of the toroidal and poloidal components of the magnetic field and the magnetic helicity. Full-disk magnetograms provide more accurate and complete input data by constraining both the toroidal and poloidal global field components, but these data are available only for the last four solar cycles. In this presentation, using the available magnetogram data, we discuss development of the methodology and forecast quality criteria (including forecast uncertainties and sources of errors). We demonstrate the influence of limited time series observations on the accuracy of solar activity predictions. We present EnKF predictions of the upcoming Solar Cycle 25 based on both the sunspot number series and observed magnetic fields and discuss the uncertainties and potential of the data assimilation approach.

Kitiashvili, Irina N.↗
NASA NTRS2011

MoSST DAS: The First Working Geomagnetic Data Assimilation System

The Earth possesses an internal magnetic field (geomagnetic field) generated by convection in the outer core (geodynamo). Previous efforts have been focused along two distinct paths: (1) numerical geodynamo modeling to understand the origin of the geomagnetic field, and the mechanisms of geomagnetic secular variations (SV); and (2) geomagnetic field modeling to map the spatial/temporal variations of the field from geomagnetic data, and to derive core properties, e.g. inversion of core flow near the core-mantle boundary (CMB). Geomagnetic data assimilation is a new approach emerged over the past 5 years: surface observations are assimilated with geodynamo models for better understanding of the core dynamical state, and accurately prediction of SV. In collaboration with several geomagnetic research groups, we have developed the first working geomagnetic data assimilation system, Modular, Scalable, Self-consistent, and Three-dimensional (MoSST) DAS, that includes the MoSST numerical dynamo model; 7000 years of geomagnetic field maps from several field models utilizing satellite and ground observatory data, historical magnetic records and archeo/paleo magnetic data; and an ensemble based optimal interpolation (01) assimilation algorithm. With this system, we have demonstrated clearly that the assimilated core dynamical state is substantially different from those of pure geodynamo simulations. Ensemble assimilation runs also show the convergence of the assimilated solutions inside the core, suggesting that the simulation state is pulled closer to the truth via data assimilation. The forecasts from this system are also very accurate: the 5-year forecast of the geomagnetic field agrees very well with the observations; and the 5-year secular variation forecast is more accurate than the IGRF SV forecast models in the past. Using geomagnetic records up to 2009, we have made an SV forecast for the period from 2010-2015, and is a candidate SV model for IGRF-11.

Kuang, Weijia↗
NASA NTRS2010

Assimilation of Satellite-Derived Skin Temperature Observations into Land Surface Models

Land surface (or "skin") temperature (LST) lies at the heart of the surface energy balance and is a key variable in weather and climate models. Here we assimilate LST retrievals from the International Satellite Cloud Climatology Project (ISCCP) into the Noah and Catchment (CLSM) land surface models using an ensemble-based, off-line land data assimilation system. LST is described very differently in the two models. A priori scaling and dynamic bias estimation approaches are applied because satellite and model LST typically exhibit different mean values and variability. Performance is measured against 27 months of in situ measurements from the Coordinated Energy and Water Cycle Observations Project at 48 stations. LST estimates from Noah and CLSM without data assimilation ("open loop") are comparable to each other and superior to that of ISCCP retrievals. For LST, RMSE values are 4.9 K (CLSM), 5.6 K (Noah), and 7.6 K (ISCCP), and anomaly correlation coefficients (R) are 0.62 (CLSM), 0.61 (Noah), and 0.52 (ISCCP). Assimilation of ISCCP retrievals provides modest yet statistically significant improvements (over open loop) of up to 0.7 K in RMSE and 0.05 in anomaly R. The skill of surface turbulent flux estimates from the assimilation integrations is essentially identical to the corresponding open loop skill. Noah assimilation estimates of ground heat flux, however, can be significantly worse than open loop estimates. Provided the assimilation system is properly adapted to each land model, the benefits from the assimilation of LST retrievals are comparable for both models.

Reichle, Rolf H.↗
NASA NTRS2019

Near Real-Time Sub/Seasonal Prediction of Aerosol at NASA Global Modeling and Assimilation Office

A new version of the coupled modeling and analysis system used to produce near real time subseasonal to seasonal forecasts was released over a year ago by the NASA/Goddard Global Modeling and Assimilation Office. The model runs at approximately 1/2 degree globally in the atmosphere and ocean, contains a realistic description of the cryosphere, and includes an interactive aerosol model. The data assimilation used to produce initial conditions is weakly coupled, in which the atmosphere-only assimilated state is coupled to an ocean data assimilation system using a Local Ensemble Transform Kalman Filter. Here will briefly describe the new system, and show results of aerosol-derived air quality from an extensive series of retrospective forecasts. The interactive aerosol is shown to improve seasonal time scale prediction skill during some "forecasts of oppurtunity". Plans for a future version of the system with predicted biomass burning from fires will also be discussed.

Molod, Andrea M.↗
NASA NTRS2019

Near Real-Time Sub/Seasonal Prediction of Aerosol and Air Quality at the NASA Global Modeling and Assimilation Office

Version 2 of the coupled modeling and analysis system used to produce near real time subseasonal to seasonal forecasts was released almost two years ago by the NASA/Goddard Global Modeling and Assimilation Office. The model runs at approximately 1/2 degree globally in the atmosphere and ocean, contains a realistic description of the cryosphere, and includes an interactive aerosol model. The data assimilation used to produce initial conditions is weakly coupled, in which the atmosphere-only assimilated state is coupled to an ocean data assimilation system using a Local Ensemble Transform Kalman Filter. Results of aerosol-derived air quality (Particulate Matter) from an extensive series of retrospective forecasts will be shown, with particular focus on the continental United States and eastern Asia. In addition, under some circumstances, the interactive aerosol is shown to improve seasonal time scale prediction skill. Plans for a future version of the system with predicted biomass burning from fires will also be discussed.

Molod, Andrea↗
NASA NTRS2016

Martian Polar Vortices: Comparison of Reanalyses

The structure and evolution of the Martian polar vortices is examined using two recently available reanalysis systems: version 1.0 of the Mars Analysis Correction Data Assimilation (MACDA) and a preliminary version of the Ensemble Mars Atmosphere Reanalysis System (EMARS). There is quantitative agreement between the reanalyses in the lower atmosphere, where Mars Global Surveyor (MGS) Thermal Emission Spectrometer (TES) data are assimilated, but there are differences at higher altitudes reflecting differences in the free-running general circulation model simulations used in the two reanalyses. The reanalyses show similar potential vorticity (PV) structure of the vortices: There is near-uniform small PV equatorward of the core of the westerly jet, steep meridional PV gradients on the polar side of the jet core, and a maximum of PV located off of the pole. In maps of 30 sol mean PV, there is a near-continuous elliptical ring of high PV with roughly constant shape and longitudinal orientation from fall to spring. However, the shape and orientation of the vortex varies on daily time scales, and there is not a continuous ring of PV but rather a series of smaller scale coherent regions of high PV. The PV structure of the Martian polar vortices is, as has been reported before, very different from that of Earth's stratospheric polar vortices, but there are similarities with Earth's tropospheric vortices which also occur at the edge of the Hadley Cell, and have near-uniform small PV equatorward of the jet, and a large increase of PV poleward of the jet due to increased stratification.

Waugh, D. W.↗
NASA NTRS2019

Assimilation of Aerosol Observations in the NASA GEOS Model

In the GEOS near real-time system, as well as in MERRA-2 which is the latest reanalysis produced at NASA's Global Modeling Assimilation Office (GMAO), the assimilation of aerosol observations is performed by means of a so-called analysis splitting method. The prognostic model is based on the GEOS model radiatively coupled to GOCART aerosol module and includes assimilation of bias-corrected Aerosol Optical Depth (AOD) at 550 nm from various space-based remote sensing platforms. In line with the transition of the GEOS meteorological data assimilation system to a hybrid Ensemble-Variational formulation, we are updating the aerosol component of our assimilation system to a variational ensemble type of scheme. In this talk we will examine the impact of replacing the current analysis splitting scheme with this new approach. Starting with the assimilation of satellite based single-channel retrievals; we will discuss the impact of this aerosol data assimilation technique on the 3D aerosol distributions by means of innovation statistics and verification against independent datasets such as the Aerosol Robotic Network (AERONET) and surface PM2.5. We will also present preliminary results related to the introduction of new aerosol data types in GEOS, including multi-spectral AOD retrievals.

AERONET↗
NASA NTRSDate not supplied

Improving the Representation of Land Surface Processes using the Data Assimilation Research Testbed (DART)

The land surface is a critical part of the earth system as processes related to water, carbon, energy and nitrogen cycling have important implications for climate forcing, air quality, water availability and seasonal atmospheric forecasting. Despite advances in land surface modeling, land surface model performance is often limited because of errors related to initial and boundary conditions, model structure, and parameters. Data assimilation (DA) techniques combined with an expanding network of earth system observations present an opportunity to reduce these errors and improve simulations. Here we apply an Ensemble Kalman Filter DA system as part of the Data Assimilation Research Testbed to a variety of land surface simulations. First, we describe the use of remotely sensed biomass observations to provide improved simulations of plant phenology, carbon and water cycling for regions highly sensitive to climate change (Western US, China, and Arctic). We discuss approaches to account for systemic biases between models and observations, including the use of spatially-varying adaptive ensemble inflation as an alternative approach to re-scaling soil moisture observations. Finally, we discuss a strategy to incorporate complementary observations (snow water equivalent, solar-induced fluorescence) to better constrain the representation of carbon and water cycling across complex terrain.

DART↗
NASA NTRS2012

EnKF OSSE Experiments Assessing the Impact of HIRAD Wind Speed and HIWRAP Radial Velocity Data on Analysis of Hurricane Karl (2010)

Previous studies (e.g., Zhang et al. 2009, Weng et al. 2011) have shown that radial velocity data from airborne and ground-based radars can be assimilated into ensemble Kalman filter (EnKF) systems to produce accurate analyses of tropical cyclone vortices, which can reduce forecast intensity error. Recently, wind speed data from SFMR technology has also been assimilated into the same types of systems and has been shown to improve the forecast intensity of mature tropical cyclones. Two instruments that measure these properties were present during the NASA Genesis and Rapid Intensification Processes (GRIP) field experiment in 2010 which sampled Hurricane Karl, and will next be co-located on the same aircraft for the subsequent NASA HS3 experiment. The High Altitude Wind and Rain Profiling Radar (HIWRAP) is a conically scanning Doppler radar mounted upon NASAs Global Hawk unmanned aerial vehicle, and the usefulness of its radial velocity data for assimilation has not been previously examined. Since the radar scans from above with a fairly large fixed elevation angle, it observes a large component of the vertical wind, which could degrade EnKF analyses compared to analyses with data taken from lesser elevation angles. The NASA Hurricane Imaging Radiometer (HIRAD) is a passive microwave radiometer similar to SFMR, and measures emissivity and retrieves hurricane surface wind speeds and rain rates over a much wider swath. Thus, this study examines the impact of assimilating simulated HIWRAP radial velocity data into an EnKF system, simulated HIRAD wind speed, and HIWRAP+HIRAD with the Weather Research and Forecasting (WRF) model and compares the results to no data assimilation and also to the Truth from which the data was simulated for both instruments.

Albers, Cerese↗