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 109 records · Page 6

Does Ocean Color Data Assimilation Improve Estimates of Global Ocean Inorganic Carbon?

Ocean color data assimilation has been shown to dramatically improve chlorophyll abundances and distributions globally and regionally in the oceans. Chlorophyll is a proxy for phytoplankton biomass (which is explicitly defined in a model), and is related to the inorganic carbon cycle through the interactions of the organic carbon (particulate and dissolved) and through primary production where inorganic carbon is directly taken out of the system. Does ocean color data assimilation, whose effects on estimates of chlorophyll are demonstrable, trickle through the simulated ocean carbon system to produce improved estimates of inorganic carbon? Our emphasis here is dissolved inorganic carbon, pC02, and the air-sea flux. We use a sequential data assimilation method that assimilates chlorophyll directly and indirectly changes nutrient concentrations in a multi-variate approach. The results are decidedly mixed. Dissolved organic carbon estimates from the assimilation model are not meaningfully different from free-run, or unassimilated results, and comparisons with in situ data are similar. pC02 estimates are generally worse after data assimilation, with global estimates diverging 6.4% from in situ data, while free-run estimates are only 4.7% higher. Basin correlations are, however, slightly improved: r increase from 0.78 to 0.79, and slope closer to unity at 0.94 compared to 0.86. In contrast, air-sea flux of C02 is noticeably improved after data assimilation. Global differences decline from -0.635 mol/m2/y (stronger model sink from the atmosphere) to -0.202 mol/m2/y. Basin correlations are slightly improved from r=O.77 to r=0.78, with slope closer to unity (from 0.93 to 0.99). The Equatorial Atlantic appears as a slight sink in the free-run, but is correctly represented as a moderate source in the assimilation model. However, the assimilation model shows the Antarctic to be a source, rather than a modest sink and the North Indian basin is represented incorrectly as a sink rather than the source indicated by the free-run model and data estimates.

Gregg, Watson↗

The Impact of Ocean Observations in Seasonal Climate Prediction

The ocean provides the most significant memory for the climate system. Hence, a critical element in climate forecasting with coupled models is the initialization of the ocean with states from an ocean data assimilation system. Remotely-sensed ocean surface fields (e.g., sea surface topography, SST, winds) are now available for extensive periods and have been used to constrain ocean models to provide a record of climate variations. Since the ocean is virtually opaque to electromagnetic radiation, the assimilation of these satellite data is essential to extracting the maximum information content. More recently, the Argo drifters have provided unprecedented sampling of the subsurface temperature and salinity. Although the duration of this observation set has been too short to provide solid statistical evidence of its impact, there are indications that Argo improves the forecast skill of coupled systems. This presentation will address the impact these different observations have had on seasonal climate predictions with the GMAO's coupled model.

Rienecker, Michele↗

A High-Resolution Merged Wind Dataset for DYNAMO: Progress and Future Plans

In order to support research on optimal data assimilation methods for the Cyclone Global Navigation Satellite System (CYGNSS), launching in 2016, work has been ongoing to produce a high‐resolution merged wind dataset for the Dynamics of the Madden Julian Oscillation (DYNAMO) field campaign, which took place during late 2011/early 2012. The winds are produced by assimilating DYNAMO observations into the Weather Research and Forecasting (WRF) three‐dimensional variational (3DVAR) system. Data sources from the DYNAMO campaign include the upper‐air sounding network, radial velocities from the radar network, vector winds from the Advanced Scatterometer (ASCAT) and Oceansat‐2 Scatterometer (OSCAT) satellite instruments, the NOAA High Resolution Doppler Lidar (HRDL), and several others. In order the prep them for 3DVAR, significant additional quality control work is being done for the currently available TOGA and SMART‐R radar datasets, including automatically dealiasing radial velocities and correcting for intermittent TOGA antenna azimuth angle errors. The assimilated winds are being made available as model output fields from WRF on two separate grids with different horizontal resolutions ‐ a 3‐km grid focusing on the main DYNAMO quadrilateral (i.e., Gan Island, the R/V Revelle, the R/V Mirai, and Diego Garcia), and a 1‐km grid focusing on the Revelle. The wind dataset is focused on three separate approximately 2‐week periods during the Madden Julian Oscillation (MJO) onsets that occurred in October, November, and December 2011. Work is ongoing to convert the 10‐m surface winds from these model fields to simulated CYGNSS observations using the CYGNSS End‐To‐End Simulator (E2ES), and these simulated satellite observations are being compared to radar observations of DYNAMO precipitation systems to document the anticipated ability of CYGNSS to provide information on the relationships between surface winds and oceanic precipitation at the mesoscale level. This research will improve our understanding of the future utility of CYGNSS for documenting key MJO processes.

Lang, Timothy J.↗

Analysis of Terrestrial Water Storage Changes from GRACE and GLDAS

Since March 2002, the Gravity Recovery and Climate Experiment (GRACE) has provided first estimates of land water storage variations by monitoring the time-variable component of Earth's gravity field. Here we characterize spatial-temporal variations in terrestrial water storage changes (TWSC) from GRACE and compare them to those simulated with the Global Land Data Assimilation System (GLDAS). Additionally, we use GLDAS simulations to infer how TWSC is partitioned into snow, canopy water and soil water components, and to understand how variations in the hydrologic fluxes act to enhance or dissipate the stores. Results quantify the range of GRACE-derived storage changes during the studied period and place them in the context of seasonal variations in global climate and hydrologic extremes including drought and flood, by impacting land memory processes. The role of the largest continental river basins as major locations for freshwater redistribution is highlighted. GRACE-based storage changes are in good agreement with those obtained from GLDAS simulations. Analysis of GLDAS-simulated TWSC illustrates several key characteristics of spatial and temporal land water storage variations. Global averages of TWSC were partitioned nearly equally between soil moisture and snow water equivalent, while zonal averages of TWSC revealed the importance of soil moisture storage at low latitudes and snow storage at high latitudes. Evapotranspiration plays a key role in dissipating globally averaged terrestrial water storage. Latitudinal averages showed how precipitation dominates TWSC variations in the tropics, evapotranspiration is most effective in the midlatitudes, and snowmelt runoff is a key dissipating flux at high latitudes. Results have implications for monitoring water storage response to climate variability and change, and for constraining land model hydrology simulations.

Syed, Tajdarul H.↗

Comparative Analysis of Upper Ocean Heat Content Variability from Ensemble Operational Ocean Analyses

Upper ocean heat content (HC) is one of the key indicators of climate variability on many time-scales extending from seasonal to interannual to long-term climate trends. For example, HC in the tropical Pacific provides information on thermocline anomalies that is critical for the longlead forecast skill of ENSO. Since HC variability is also associated with SST variability, a better understanding and monitoring of HC variability can help us understand and forecast SST variability associated with ENSO and other modes such as Indian Ocean Dipole (IOD), Pacific Decadal Oscillation (PDO), Tropical Atlantic Variability (TAV) and Atlantic Multidecadal Oscillation (AMO). An accurate ocean initialization of HC anomalies in coupled climate models could also contribute to skill in decadal climate prediction. Errors, and/or uncertainties, in the estimation of HC variability can be affected by many factors including uncertainties in surface forcings, ocean model biases, and deficiencies in data assimilation schemes. Changes in observing systems can also leave an imprint on the estimated variability. The availability of multiple operational ocean analyses (ORA) that are routinely produced by operational and research centers around the world provides an opportunity to assess uncertainties in HC analyses, to help identify gaps in observing systems as they impact the quality of ORAs and therefore climate model forecasts. A comparison of ORAs also gives an opportunity to identify deficiencies in data assimilation schemes, and can be used as a basis for development of real-time multi-model ensemble HC monitoring products. The OceanObs09 Conference called for an intercomparison of ORAs and use of ORAs for global ocean monitoring. As a follow up, we intercompared HC variations from ten ORAs -- two objective analyses based on in-situ data only and eight model analyses based on ocean data assimilation systems. The mean, annual cycle, interannual variability and longterm trend of HC have been analyzed

Xue, Yan↗

Dynamic Responses of the Earth's Outer Core to Assimilation of Observed Geomagnetic Secular Variation

Assimilation of surface geomagnetic observations and geodynamo models has advanced very quickly in recent years. However, compared to advanced data assimilation systems in meteorology, geomagnetic data assimilation (GDAS) is still in an early stage. Among many challenges ranging from data to models is the disparity between the short observation records and the long time scales of the core dynamics. To better utilize available observational information, we have made an effort in this study to directly assimilate the Gauss coefficients of both the core field and its secular variation (SV) obtained via global geomagnetic field modeling, aiming at understanding the dynamical responses of the core fluid to these additional observational constraints. Our studies show that the SV assimilation helps significantly to shorten the dynamo model spin-up process. The flow beneath the core-mantle boundary (CMB) responds significantly to the observed field and its SV. The strongest responses occur in the relatively small scale flow (of the degrees L is approx. 30 in spherical harmonic expansions). This part of the flow includes the axisymmetric toroidal flow (of order m = 0) and non-axisymmetric poloidal flow with m (is) greater than 5. These responses can be used to better understand the core flow and, in particular, to improve accuracies of predicting geomagnetic variability in future.

secular variation↗

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↗

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↗

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↗

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↗

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

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

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↗

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

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

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↗

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↗

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↗