Search NASA⌕ Search

SEARCH · Search NASA

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 73 records · Page 4

Stratospheric Tides and Data Assimilation

In the upper stratosphere, the atmosphere exhibits significant diurnal and semi-diurnal tidal variations, with typical amplitude of about 2K in mid-latitudes. In this paper we examine how well the tidal variations in temperature are represented by the Goddard Geodesic Earth Orbiting Satellite (GEOS-2) data assimilation system. We show that the GEOS-2 atmospheric model is quite successful at simulating the tidal temperature variations. However, the assimilation of satellite temperature soundings significantly damps the simulated tides. The reason is because the tides are not well represented by the satellite retrievals used by the assimilation system (which have a typical tidal amplitude of around 1K). As a result of this study, we suggest improvements that should be made to the treatment of satellite soundings by the assimilation system.

Swinbank, R.↗

Uncertainty Estimation for SMAP Level-1 Brightness Temperature Assimilation at Different Timescales

In this study, we assimilate Soil Moisture Active Passive (SMAP) mission brightness temperature (Tb) observations into NASA’s Catchment Land Surface Model using an Ensemble Kalman filter to update surface and root-zone soil moisture simulations. Different time series components of the Tb observations are assimilated including anomalies, inter-annual variations, and high-frequency variations. To optimize the weights that the data assimilation (DA) puts on the observations, the ratio between the un-certainties of modeled and observed Tb is approximated using modeled and observed soil moisture uncertainties estimated using triple collocation analysis. Results are compared to a benchmark experiment that mimics the operational SMAP Level-4 algorithm, which assimilates Tb observations using a spatially-constant 4 Kelvin (K) observation uncertainty. All DA experiments exhibit notable skill improvements in most regions. Improvements are greatest for the inter-annual variations in the simulations of both surface and root-zone soil moisture (mean improvements in terms of Pearson correlation (-) are 0.08 and 0.06, respectively). Anomaly simulations improve similarly (0.07), and improvements in the high-frequency variations are only observed for sur-face soil moisture simulations (0.06). Strikingly, however, no notable difference in skill—neither improvement nor deterioration—is observed between the experiments that use optimized observation uncertainty parameters and the 4 K benchmark experiment. We show, analytically, that this may be explained by the presence of large observation operator errors, which have the potential to render post-update uncertainty insensitive to inaccuracies in the Kalman gain.

Alexander Gruber↗

Hyperspectral Remote Sensing of Atmosphere and Surface Properties

Atmospheric Infrared Sounder (AIRS), Infrared Atmospheric Sounding Interferometer (IASI), and Cross-track Infrared Sounder (CrIS) are all hyper-spectral satellite sensors with thousands of spectral channels. Top of atmospheric radiance spectra measured by these sensors contain high information content on atmospheric, cloud, and surface properties. Exploring high information content contained in these high spectral resolution spectra is a challenging task due to computation e ort involved in modeling thousands of spectral channels. Usually, only very small fractions (4{10 percent) of the available channels are included in physical retrieval systems or numerical weather forecast (NWP) satellite data assimilations. We will describe a method of simultaneously retrieving atmospheric temperature, moisture, cloud, and surface properties using all available spectral channels without sacrificing computational speed. The essence of the method is to convert channel radiance spectra into super-channels by an Empirical Orthogonal Function (EOF) transformation. Because the EOFs are orthogonal to each other, about 100 super-channels are adequate to capture the information content of the radiance spectra. A Principal Component-based Radiative Transfer Model (PCRTM) developed at NASA Langley Research Center is used to calculate both the super-channel magnitudes and derivatives with respect to atmospheric profiles and other properties. There is no need to perform EOF transformations to convert super channels back to spectral space at each iteration step for a one-dimensional variational retrieval or a NWP data assimilation system. The PCRTM forward model is also capable of calculating radiative contributions due to multiple-layer clouds. The multiple scattering effects of the clouds are efficiently parameterized. A physical retrieval algorithm then performs an inversion of atmospheric, cloud, and surface properties in super channel domain directly therefore both reducing the computational need and preserving the information content of the IASI measurements. The inversion algorithm is based on a non-linear Levenberg-Marquardt method with climatology covariance matrices and a priori information as constraints. One advantage of this approach is that it uses all information content from the hyper-spectral data so that the retrieval is less sensitive to instrument noise and eliminates the need for selecting a sub-set of the channels.

Liu, Xu↗

Next generation initiation techniques

Four-dimensional data assimilation strategies can generally be classified as either current or next generation, depending upon whether they are used operationally or not. Current-generation data-assimilation techniques are those that are presently used routinely in operational-forecasting or research applications. They can be classified into the following categories: intermittent assimilation, Newtonian relaxation, and physical initialization. It should be noted that these techniques are the subject of continued research, and their improvement will parallel the development of next generation techniques described by the other speakers. Next generation assimilation techniques are those that are under development but are not yet used operationally. Most of these procedures are derived from control theory or variational methods and primarily represent continuous assimilation approaches, in which the data and model dynamics are 'fitted' to each other in an optimal way. Another 'next generation' category is the initialization of convective-scale models. Intermittent assimilation systems use an objective analysis to combine all observations within a time window that is centered on the analysis time. Continuous first-generation assimilation systems are usually based on the Newtonian-relaxation or 'nudging' techniques. Physical initialization procedures generally involve the use of standard or nonstandard data to force some physical process in the model during an assimilation period. Under the topic of next-generation assimilation techniques, variational approaches are currently being actively developed. Variational approaches seek to minimize a cost or penalty function which measures a model's fit to observations, background fields and other imposed constraints. Alternatively, the Kalman filter technique, which is also under investigation as a data assimilation procedure for numerical weather prediction, can yield acceptable initial conditions for mesoscale models. The third kind of next-generation technique involves strategies to initialize convective scale (non-hydrostatic) models.

Warner, Tom↗

Exploring the Utility of the Planned CYGNSS Mission for Investigating the Initiation and Development of the Madden-Julian Oscillation

CYGNSS is a planned constellation consisting of multiple micro-satellites that leverage the Global Positioning System (GPS) to provide rapidly updated, high resolution (approx. 15-50 km, approx. 4 h) surface wind speeds (via bi-static scatterometry) over the tropical oceans in any weather condition, including heavy rainfall. The approach of the work to be presented at this conference is to utilize a limited-domain, cloud-system resolving model (Weather Research and Forecasting or WRF) and its attendant data assimilation scheme (Three-Dimensional Variational Assimilation or 3DVAR) to investigate the utility of the CYGNSS mission for helping characterize key convectiveto- mesoscale processes - such as surface evaporation, moisture advection and convergence, and upscale development of precipitation systems - that help drive the initiation and development of the Madden-Julian Oscillation (MJO) in the equatorial Indian Ocean. The proposed work will focus on three scientific objectives. Objective 1 is to produce a high-resolution surface wind dataset resolution (approx. 0.5 h, approx. 1-4 km) for multiple MJO onsets using WRF-assimilated winds and other data from the DYNAmics of the MJO (DYNAMO) field campaign, which took place during October 2011 - March 2012. Objective 2 is to study the variability of surface winds during MJO onsets at temporal and spatial scales of finer resolution than future CYGNSS data. The goal is to understand how sub-CYGNSS-resolution processes will shape the observations made by the satellite constellation. Objective 3 is to ingest simulated CYGNSS data into the WRF model in order to perform observing system simulation experiments (OSSEs). These will be used to test and quantify the potential beneficial effects provided by CYGNSS, particularly for characterizing the physical processes driving convective organization and upscale development during the initiation and development of the MJO. The proposed research is ideal for answering important questions about the CYGNSS mission, such as the representativeness of surface wind retrievals in the context of the complex airflow processes that occur during heavy precipitation, as well as the tradeoffs in retrieval accuracy that result from finer spatial resolution of the CYGNSS winds versus increased errors/noisiness in those data. Research plans and initial progress toward these objectives will be presented.

Lang, Timothy↗

The Special Sensor Microwave Imager Wind Dataset: A Source of Quantitative Information for the Ocean-to-Land Advection

It is well recognized that advection from the North Atlantic has a profound effect on the climatic conditions in central Europe. A new dataset of the ocean-surface winds, derived from the Special Sensor Microwave Imager, SSM/1, is now available. This satellite instrument measures the wind speed, but not the direction. However, variational analysis developed at the Data Assimilation Office, NASA Goddard Space Flight Center, by combining the SSM/I measurements with wind vectors measured from ships, etc., produced global maps of the ocean surface winds suitable for climate analysis. From this SSM/I dataset, a specific index I(sub na) of the North Atlantic surface winds has been developed, which pertinently quantifies the low-level advection into central Europe. For a selected time-period, the index I(sub na) reports the average of the amplitude of the wind, averaging only the speed when the direction is from the southwest (when the wind is from another direction, the contribution counts to the average as zero speed). Strong correlations were found between February I(sub na) and the surface air temperatures in Europe 50-60 deg N. In the present study, we present the correlations between I(sub na) and temperature I(sub s), and also the sensitivity of T(sub s), to an increase in I(sub na), in various seasons and various regions. We specifically analyze the flow of maritime-air from the North Atlantic that produced two extraordinary warm periods: February 1990, and early-winter 2000/2001. The very cold December 2001 was clearly due to a northerly flow. Our conclusion is that the SSM/I dataset is very useful for providing insight to the forcing of climatic fluctuations in Europe.

Otterman, J.↗

Investigation of the Potential Saturation of Information from Global Navigation Satellite System Radio Occultation Observations with an Observing System Simulation Experiment

The potential impact of large numbers of Global Navigation Satellite System radio occultation (GNSS-RO) observations on numerical weather prediction is investigated using a global observing system simulation experiment (OSSE). The hybrid four-dimensional ensemble variational Gridpoint Statistical Interpolation (GSI) data assimilation system and Global Earth Observing System (GEOS) model are used to ingest up to 100,000 GNSS-RO soundings per day in addition to the current suite of conventional and radiance data. Analysis quality, forecast skill, and forecast sensitivity to observation impact are examined with differing quantities of additional GNSS-RO profiles. It is found that saturation of information from additional RO soundings has not been reached with 100,000 soundings per day. There are some indications of suboptimal performance of the GSI in handling GNSS-RO observations particularly in the middle and lower tropospheric extratropics.

GNSS-RO↗

Data Assimilation and Uncertainties in Early Solar Cycle Predictions

Stochastic nature of solar activity variations together with our limited knowledge of the dynamo mechanism and subsurface dynamics causes uncertainty in predictions of the solar cycle. For improving the physics-­‐based predictions we can take advantage of the mathematical data assimilation approach that allows us to take into account both, observational errors and model uncertainties, and provide estimates of the next solar cycle along with prediction uncertainties. In this study we use the Parker's migratory dynamo model together with the equation of magnetic helicity balance, which reproduces main properties of the sunspot cycles and allow us to minimize discrepancies between the observed global activity variations and the model solution. The test simulation runs show that a reliable prediction can be obtained for two phases of preceding solar cycle: 1) if the polar field reversals shortly after the solar maxima (strong toroidal field and weak poloidal field), and 2) during the solar minima (strongest poloidal and weak toroidal fields). The early estimate of Cycle 25 obtained by this method shows that this cycle will start in 2019 – 2020, reach the maximum in 2023 – 2024, and that the mean sunspot number at the maximum will be about 90 (for the v2.0 sunspot number series).

Data Assimilation↗

Observing System Variations Effect on Reanalyses

Reanalyses integrate multitudes of satellite and conventional observations data assimilation and numerical weather prediction. The result is that many disparate observation platforms, discontinuous in space and time, lead to complete and consistent representations the state of the weather. The component also provides physical fields rarely or never observed. However, the numerical model bias is continuously being corrected by the observational analysis, and this bias changes as variations in the observations occur. NASA's Modern Era Retrospective-analysis for Research and Applications (MERRA) sensitivity to variations in the observing systems are explored. Specifically, we will evaluate the water budget and transport processes as they relate to the advent of SSM/I and AMSU-A radiance assimilation, and an additional case of radiosonde station that exhibits a dramatic shift in mean water states. The MERRA input observation data, now available online, is used to explore these variations.

Bosilovich, Michael G.↗

Variational Assimilation of Global Microwave Rainfall Retrievals: Physical and Dynamical Impact on GEOS Analyses and Forecasts

Global microwave rainfall retrievals from a 5-satellite constellation, including TMI from TRMM, SSWI from DMSP F13, F14 and F15, and AMSR-E from EOS-AQUA, are assimilated into the NASA Goddard Earth Observing System (GEOS) Data Assimilation System (DAS) using a 1-D variational continuous assimilation (VCA) algorithm. The physical and dynamical impact of rainfall assimilation on GEOS analyses and forecasts is examined at various temporal and spatial scales. This study demonstrates that the 1-D VCA algorithm, which was originally developed and evaluated for rainfall assimilations over tropical oceans, can effectively assimilate satellite microwave rainfall retrievals and improve GEOS analyses over both the Tropics and the extratropics where the atmospheric processes are dominated by different large-scale dynamics and moist physics, and also over the land, where rainfall estimates from passive microwave radiometers are believed to be less accurate. Results show that rainfall assimilation renders the GEOS analysis physically and dynamically more consistent with the observed precipitation at the monthly-mean and 6-hour time scales. Over regions where the model precipitation tends to misbehave in distinctly different rainy regimes, the 1-D VCA algorithm, by compensating for errors in the model s moist time-tendency in a 6-h analysis window, is able to bring the rainfall analysis closer to the observed. The radiation and cloud fields also tend to be in better agreement with independent satellite observations in the rainfall-assimilation m especially over regions where rainfall analyses indicate large improvements. Assimilation experiments with and without rainfall data for a midlatitude frontal system clearly indicates that the GEOS analysis is improved through changes in the thermodynamic and dynamic fields that respond to the rainfall assimilation. The synoptic structures of temperature, moisture, winds, divergence, and vertical motion, as well as vorticity are more realistically captured across the front. Short-term forecasts using initial conditions assimilated with rainfall data also show slight improvements. 1

Lin, Xin↗

Geomagnetic Secular Variation Prediction with Thermal Heterogeneous Boundary Conditions

It has long been conjectured that thermal heterogeneity at the core-mantle boundary (CMB) affects the geodynamo substantially. The observed two pairs of steady and strong magnetic flux lobes near the Polar Regions and the low secular variation in the Pacific over the past 400 years (and perhaps longer) are likely the consequences of this CMB thermal heterogeneity. There are several studies on the impact of the thermal heterogeneity with numerical geodynamo simulations. However, direct correlation between the numerical results and the observations is found very difficult, except qualitative comparisons of certain features in the radial component of the magnetic field at the CMB. This makes it difficult to assess accurately the impact of thermal heterogeneity on the geodynamo and the geomagnetic secular variation. We revisit this problem with our MoSST_DAS system in which geomagnetic data are assimilated with our geodynamo model to predict geomagnetic secular variations. In this study, we implement a heterogeneous heat flux across the CMB that is chosen based on the seismic tomography of the lowermost mantle. The amplitude of the heat flux (relative to the mean heat flux across the CMB) varies in the simulation. With these assimilation studies, we will examine the influences of the heterogeneity on the forecast accuracies, e.g. the accuracies as functions of the heterogeneity amplitude. With these, we could be able to assess the model errors to the true core state, and thus the thermal heterogeneity in geodynamo modeling.

Kuang, Weijia↗

Diverse Photosynthetic Capacity of Global Ecosystems Mapped by Satellite Chlorophyll Fluorescence Measurements

Photosynthetic capacity is often quantified by the Rubisco-limited photosynthetic capacity (i.e. maximum carboxylation rate, V(sub cmax)). It is a key plant functional trait that is widely used in Earth System Models for simulation of the global carbon and water cycles. Measuring V(sub cmax) is time-consuming and laborious; therefore, the spatiotemporal distribution of V(sub cmax) is still poorly understood due to limited measurements of V(sub cmax). In this study, we used a data assimilation approach to map the spatial variation of V(sub cmax) for global terrestrial ecosystems from a 11-year-long satellite-observed solar-induced chlorophyll fluorescence (SIF) record. In this SIF-derived V(sub cmax) map, the mean V(sub cmax) value for each plant function type (PFT) is found to be comparable to a widely used N-derived V(sub cmax) dataset by Kattge et al. (2009). The gradient of V(sub cmax) along PFTs is clearly revealed even without land cover information as an input. Large seasonal and spatial variations of V(sub cmax) are found within each PFT, especially for diverse crop rotation systems. The distribution of major crop belts, characterized with high V(sub cmax) values, is highlighted in this V(sub cmax) map. Legume plants are characterized with high V(sub cmax) values. This V(sub cmax) map also clearly illustrates the emerging soybean revolution in South America where V(sub cmax) is the highest among the world. The gradient of V(sub cmax) in Amazon is found to follow the transition of soil types with different soil N and P contents. This study suggests that satellite-observed SIF is powerful in deriving the important plant functional trait, i.e. V(sub cmax), for global climate change studies.

Liming He↗

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↗

Stratospheric Water Vapor Observations From the International Space Station By the Stratospheric Aerosol and Gas Experiment III (Sage III/ISS)

The Stratospheric Aerosol and Gas Experiment III has routinely observed the upper atmosphere from the International Space Station (SAGE III/ISS) since June 2017. An important greenhouse gas and molecule affecting chemistry cycles, water vapor is among the vertical profiles of trace gases measured by SAGE III/ISS. This makes SAGE III/ISS the youngest of the current fleet of space-based instruments capable of measuring the vertical distribution of stratospheric water vapor (SWV). The SWV record is mature, a candidate for multi-instrument datasets, e.g. Stratospheric Water and Ozone Satellite Homogenized (SWOOSH) dataset, and shows promise within the GEOS Constituent Data Assimilation System. This presentation reviews the variations in SWV observed by SAGE III/ISS over the past half-decade. Variations highlighted include phase changes of the Quasi-biennial Oscillation (QBO), seasonal and monsoonal moisture patterns, and terrestrial events such as the recent January 2022 Hunga Tonga eruption that perturbed the SWV in spectacular fashion.

stratospheric water vapor↗

The Gmao Hybrid 4d-Envar Observing System Simulation Experiment Framework.

This work describes the extension of the Global Modeling and Assimilation Office (GMAO) Observing System Simulation Experiment (OSSE) framework to use a hybrid 4D-EnVar scheme instead of 3D-Var. The original 3D-Var and hybrid 4D-EnVar OSSEs use the same version of the data assimilation system (DAS) so that a direct comparison is possible in terms of the validation with respect to their corresponding real cases. Rather than quantifying the differences between the two data assimilation methodologies, a short inter-comparison of upgrading from a 3D- to a 4D-OSSE is provided to highlight aspects where this change matters to the OSSE community and to identify features of data assimilation that can only be explored in a four-dimensional OSSE framework. A short validation of the hybrid 4D-EnVar OSSE shows that conclusions from previous assessments of the 3D-Var OSSE in its ability to mimic the behavior of the real system still hold with the same caveats. Furthermore, some aspects of the ensemble configuration and behavior are discussed along with forecast sensitivity to observation impacts (FSOI). Estimates of error standard deviations are shown to be smaller in the hybrid 4D-EnVar OSSE but with little impact on the character of the error. A discussion on future work directions focuses on exploring the four-dimensional aspect such as the error distribution within the assimilation window or four-dimensional handling of high-temporal density observations.

Data assimilation↗

1DVAR Analysis of Temperature and Humidity using GPS Radio Occulation Data

The Global Positioning System (GPS) enables positioning with a very small receiver. The signals transmitted by the GPS satellites are sensitive to the atmosphere and can be used to perform soundings with the radio occultation technique (e.g., Kursinski et al., 1997). The GPS signal can be converted to refractivity N via the Abel transform. The refractivity can then be related to atmospheric pressure P, temperature T, and water vapor partial pressure P (sub w) the GPS measurement, (between 0.5 and 1.5 km), its self-calibration, and it's nearly all-weather capabilities make it a good candidate for use in data assimilation systems (DAS) and numerical weather prediction (NWP). In order to demonstrate its usefulness in a DAS or NWP system, a first step is to assess its impact oil the analysis. A one-dimensional variational off-line analysis (1DVAR), meaning the data are not assimilated 'In the 3D DAS, constitutes a starting approach to which further enhancements can be made. The chosen observable to be analyzed in this study is the refractivity. One way to extract temperature (humidity) from the refractivity, is to assume a humidity (temperature) profile. One variable may then be retrieved without any a priori information. The 1DVAR approach used here resolves the ambiguity problem raised in the interpretation of these data. It enables retrieving these two atmospheric variables at a reasonable computing cost.

Poli, Paul↗

On the Possibilities of Predicting Geomagnetic Secular Variation with Geodynamo Modeling

We use our MoSST core dynamics model and geomagnetic field at the core-mantle boundary (CMB) continued downward from surface observations to investigate possibilities of geomagnetic data assimilation, so that model results and current geomagnetic observations can be used to predict geomagnetic secular variation in future. As the first attempt, we apply data insertion technique to examine evolution of the model solution that is modified by geomagnetic input. Our study demonstrate that, with a single data insertion, large-scale poloidal magnetic field obtained from subsequent numerical simulation evolves similarly to the observed geomagnetic variation, regardless of the initial choice of the model solution (so long it is a well developed numerical solution). The model solution diverges on the time scales on the order of 60 years, similar to the time scales of the torsional oscillations in the Earth's core. Our numerical test shows that geomagnetic data assimilation is promising with our MoSST model.

Kuang, Wei-Jia↗

Assimilation of GRACE Terrestrial Water Storage Data into a Land Surface Model: Results for the Mississippi River Basin

NASA's GRACE mission has the potential to be extremely valuable for water resources applications and global water cycle research. What makes GRACE unique among Earth Science satellite systems is that it is able to monitor variations in water stored in all forms, from snow and surface water to soil moisture to groundwater in the deepest aquifers. However, the space and time resolutions of GRACE observations are coarse. GRACE typically resolves water storage changes over regions the size of Nebraska on a monthly basis, while city-scale, daily observations would be more useful for water management, agriculture, and weather prediction. High resolution numerical (computer) hydrology models have been developed, which predict the fates of water and energy after they strike the land surface as precipitation and sunlight. These are similar to weather and climate forecast models, which simulate atmospheric processes. We integrated the GRACE observations into a hydrology model using an advanced technique called data assimilation. The results were new estimates of groundwater, soil moisture, and snow variations, which combined the veracity of GRACE with the high resolution of the model. We tested the technique over the Mississippi River basin, but it will be even more valuable in parts of the world which lack reliable data on water availability.

Zaitchik, Benjamin F.↗