Search Results

You are looking at 41 - 50 of 66 items for

  • Author or Editor: Jeffrey Anderson x
  • All content x
Clear All Modify Search
Lili Lei, Jeffrey L. Anderson, and Glen S. Romine

Abstract

For ensemble-based data assimilation, localization is used to limit the impact of observations on physically distant state variables to reduce spurious error correlations caused by limited ensemble size. Traditionally, the localization value applied is spatially homogeneous. Yet there are potentially larger errors and different covariance length scales in precipitation systems, and that may justify the use of different localization functions for precipitating and nonprecipitating regions. Here this is examined using empirical localization functions (ELFs). Using output from an ensemble observing system simulation experiment (OSSE), ELFs provide estimates of horizontal and vertical localization for different observation types in regions with and without precipitation. For temperature and u- and υ-wind observations, the ELFs for precipitating regions are shown to have smaller horizontal localization scales than for nonprecipitating regions. However, the ELFs for precipitating regions generally have larger vertical localization scales than for nonprecipitating regions. The ELFs are smoothed and then applied in three additional OSSEs. Spatially homogeneous ELFs are found to improve performance relative to a commonly used localization function with compact support. When different ELFs are applied in precipitating and nonprecipitating regions, performance is further improved, but varying ELFs by observation type was not found to be as important. Imbalance in initial states caused by use of different localization functions is diagnosed by the domain-averaged surface pressure tendency. Forecasts from analyses with ELFs have smaller surface pressure tendencies than the standard localization, indicating improved initial balance with ELFs.

Full access
Man-Yau Chan, Jeffrey L. Anderson, and Xingchao Chen

Abstract

The introduction of infrared water vapor channel radiance ensemble data assimilation (DA) has improved numerical weather forecasting at operational centers. Further improvements might be possible through extending ensemble data assimilation methods to better assimilate infrared satellite radiances. Here, we will illustrate that ensemble statistics under clear-sky conditions are different from cloudy conditions. This difference suggests that extending the ensemble Kalman filter (EnKF) to handle bi-Gaussian prior distributions may yield better results than the standard EnKF. In this study, we propose a computationally efficient bi-Gaussian ensemble Kalman filter (BGEnKF) to handle bi-Gaussian prior distributions. As a proof-of-concept, we used the 40-variable Lorenz 1996 model as a proxy to examine the impacts of assimilating infrared radiances with the BGEnKF and EnKF. A nonlinear observation operator that constructs radiance-like bimodal ensemble statistics was used to generate and assimilate pseudoradiances. Inflation was required for both methods to effectively assimilate pseudoradiances. In both 800- and 20-member experiments, the BGEnKF generally outperformed the EnKF. The relative performance of the BGEnKF with respect to the EnKF improved when the observation spacing and time between DA cycles (cycling interval) are increased from small values. The relative performance then degraded when observation spacing and cycling interval become sufficiently large. The BGEnKF generated less noise than the EnKF, suggesting that the BGEnKF produces more balanced analysis states than the EnKF. This proof-of-concept study motivates future investigation into using the BGEnKF to assimilate infrared observations into high-order numerical weather models.

Restricted access
Thomas M. Hamill, Jeffrey S. Whitaker, Jeffrey L. Anderson, and Chris Snyder
Full access
Jeffrey Anderson, Huug van den Dool, Anthony Barnston, Wilbur Chen, William Stern, and Jeffrey Ploshay

A statistical model and extended ensemble integrations of two atmospheric general circulation models (GCMs) are used to simulate the extratropical atmospheric response to forcing by observed SSTs for the years 1980 through 1988. The simulations are compared to observations using the anomaly correlation and root-mean-square error of the 700-hPa height field over a region encompassing the extratropical North Pacific Ocean and most of North America. On average, the statistical model is found to produce considerably better simulations than either numerical model, even when simple statistical corrections are used to remove systematic errors from the numerical model simulations. In the mean, the simulation skill is low, but there are some individual seasons for which all three models produce simulations with good skill.

An approximate upper bound to the simulation skill that could be expected from a GCM ensemble, if the model's response to SST forcing is assumed to be perfect, is computed. This perfect model predictability allows one to make some rough extrapolations about the skill that could be expected if one could greatly improve the mean response of the GCMs without significantly impacting the variance of the ensemble. These perfect model predictability skills are better than the statistical model simulations during the summer, but for the winter, present-day statistical forecasts already have skill that is as high as the upper bound for the GCMs. Simultaneous improvements to the GCM mean response and reduction in the GCM ensemble variance would be required for these GCMs to do significantly better than the statistical model in winter. This does not preclude the possibility that, as is presently the case, a statistical blend of GCM and statistical predictions could produce a simulation better than either alone.

Because of the primitive state of coupled ocean–atmosphere GCMs, the vast majority of seasonal predictions currently produced by GCMs are performed using a two-tiered approach in which SSTs are first predicted and then used to force an atmospheric model; this motivates the examination of the simulation problem. However, it is straightforward to use the statistical model to produce true forecasts by changing its predictors from simultaneous to precursor SSTs. An examination of the decrease in skill of the statistical model when changed from simulation to prediction mode is extrapolated to draw conclusions about the skill to be expected from good coupled GCM predictions.

Full access
Robert Pincus, Robert J. Patrick Hofmann, Jeffrey L. Anderson, Kevin Raeder, Nancy Collins, and Jeffrey S. Whitaker

Abstract

This paper explores the degree to which short-term forecasts with global models might be improved if clouds were fully included in a data assimilation system, so that observations of clouds affected all parts of the model state and cloud variables were adjusted during assimilation. The question is examined using a single ensemble data assimilation system coupled to two present-generation climate models with different treatments of clouds. “Perfect-model” experiments using synthetic observations, taken from a free run of the model used in subsequent assimilations, are used to circumvent complications associated with systematic model errors and observational challenges; these provide a rough upper bound on the utility of cloud observations with these models. A series of experiments is performed in which direct observations of the model’s cloud variables are added to the suite of observations being assimilated. In both models, observations of clouds reduce the 6-h forecast error, with much greater reductions in one model than in the other. Improvements are largest in regions where other observations are sparse. The two cloud schemes differ in their complexity and number of degrees of freedom; the model using the simpler scheme makes better use of the cloud observations because of the stronger correlations between cloud-related and dynamical variables (particularly temperature). This implies that the impact of real cloud observations will depend on both the strength of the instantaneous, linear relationships between clouds and other fields in the natural world, and how well each assimilating model’s cloud scheme represents those relationships.

Full access
Michael K. Tippett, Jeffrey L. Anderson, Craig H. Bishop, Thomas M. Hamill, and Jeffrey S. Whitaker

Abstract

Ensemble data assimilation methods assimilate observations using state-space estimation methods and low-rank representations of forecast and analysis error covariances. A key element of such methods is the transformation of the forecast ensemble into an analysis ensemble with appropriate statistics. This transformation may be performed stochastically by treating observations as random variables, or deterministically by requiring that the updated analysis perturbations satisfy the Kalman filter analysis error covariance equation. Deterministic analysis ensemble updates are implementations of Kalman square root filters. The nonuniqueness of the deterministic transformation used in square root Kalman filters provides a framework to compare three recently proposed ensemble data assimilation methods.

Full access
Nedjeljka Žagar, Jeffrey Anderson, Nancy Collins, Timothy Hoar, Kevin Raeder, Lili Lei, and Joseph Tribbia

Abstract

Global data assimilation systems for numerical weather prediction (NWP) are characterized by significant uncertainties in tropical analysis fields. Furthermore, the largest spread of global ensemble forecasts in the short range on all scales is in the tropics. The presented results suggest that these properties hold even in the perfect-model framework and the ensemble Kalman filter data assimilation with a globally homogeneous network of wind and temperature profiles. The reasons for this are discussed by using the modal analysis, which provides information about the scale dependency of analysis and forecast uncertainties and information about the efficiency of data assimilation to reduce the prior uncertainties in the balanced and inertio-gravity dynamics.

The scale-dependent representation of variance reduction of the prior ensemble by the data assimilation shows that the peak efficiency of data assimilation is on the synoptic scales in the midlatitudes that are associated with quasigeostrophic dynamics. In contrast, the variance associated with the inertia–gravity modes is less successfully reduced on all scales. A smaller information content of observations on planetary scales with respect to the synoptic scales is discussed in relation to the large-scale tropical uncertainties that current data assimilation methodologies do not address successfully. In addition, it is shown that a smaller reduction of the large-scale uncertainties in the prior state for NWP in the tropics than in the midlatitudes is influenced by the applied radius for the covariance localization.

Full access
Jeffrey Anderson, Tim Hoar, Kevin Raeder, Hui Liu, Nancy Collins, Ryan Torn, and Avelino Avellano

The Data Assimilation Research Testbed (DART) is an open-source community facility for data assimilation education, research, and development. DART's ensemble data assimilation algorithms, careful software engineering, and diagnostic tools allow atmospheric scientists, oceanographers, hydrologists, chemists, and other geophysicists to build state-of-the-art data assimilation systems with unprecedented ease. For global numerical weather prediction, DART produces ensemble-mean analyses comparable to analyses from major centers while also providing initial conditions for ensemble predictions. In addition, DART supports more novel assimilation applications like parameter estimation, sensitivity analysis, observing system design, and smoothing. Implementing basic systems for large models requires only a few person-weeks; comprehensive systems have been built in a few months. Incorporating new observation types is also straightforward, requiring only a forward operator mapping between a model's state and an observation's expected value. Forward operators for standard, in situ observations and novel types, like GPS radio occultation soundings, are available. DART algorithms scale well on a variety of parallel architectures, allowing large data assimilation problems to be studied. DART also includes many low-order models and an ensemble assimilation tutorial appropriate for undergraduate and graduate instruction.

Full access
Ting-Chi Wu, Christopher S. Velden, Sharanya J. Majumdar, Hui Liu, and Jeffrey L. Anderson

Abstract

Recent studies have shown that assimilating enhanced satellite-derived atmospheric motion vectors (AMVs) has improved mesoscale forecast of tropical cyclones (TC) track and intensity. The authors conduct data-denial experiments to understand where the TC analyses and forecasts benefit the most from the enhanced AMV information using an ensemble Kalman filter and the Weather Research and Forecasting Model. The Cooperative Institute for Meteorological Satellite Studies at the University of Wisconsin provides enhanced AMV datasets with higher density and temporal resolution using shorter-interval image triplets for the duration of Typhoon Sinlaku and Hurricane Ike (both 2008). These AMV datasets are then spatially and vertically subsetted to create six parallel cycled assimilation-forecast experiments for each TC: all AMVs; AMVs withheld between 100 and 350 hPa (upper layer), between 350 and 700 hPa (middle layer), and between 700 and 999 hPa (lower layer); and only AMVs within (interior) and outside (exterior) 1000-km radius of the TC center. All AMV subsets are found to be useful in some capacity. The interior and upper-layer AMVs are particularly crucial for improving initial TC position, intensity, and the three-dimensional wind structure along with their forecasts. Compared with denying interior or exterior AMVs, withholding AMVs in different tropospheric layers had less impact on TC intensity and size forecasts. The ensemble forecast is less certain (larger spread) in providing accurate TC track, intensity, and size when upper-layer AMVs or interior AMVs are withheld. This information could be useful to potential targeting scenarios, such as activating and focusing satellite rapid-scan operations, and decisions regarding observing system assessments and deployments.

Full access
Hui Liu, Jeffrey Anderson, Ying-Hwa Kuo, Chris Snyder, and Alain Caya

Abstract

A nonlocal quasi-phase radio occultation (RO) observation operator is evaluated in the assimilation of Challenging Minisatellite Payload (CHAMP) radio occultation refractivity using a Weather Research and Forecasting (WRF) ensemble data assimilation system at 50-km resolution. The nonlocal operator calculates the quasi phase through integration of the model refractivity along the observed ray paths. As a comparison, a local refractivity operator that calculates the model refractivity at the observed ray perigee points is also evaluated. The assimilation is done over North America during January 2003 in two different situations: in conjunction with dense, high-quality radiosonde observations and with only satellite cloud drift wind observations. Analyses of temperature and water vapor with the RO refractivity assimilated using the local and nonlocal operator are verified against nearby withheld radiosonde observations. The bias and RMS errors of the analyses of water vapor and temperature using the nonlocal operator are significantly reduced compared with those using the local operator in the troposphere when the only additional observations are satellite cloud drift winds. The reduction of the bias and RMS errors is reduced when radiosonde observations are assimilated.

Full access