ALBERT

Description: Two perturbation generation schemes, the ensemble transformation with rescaling (ETR) and the ensemble Kalman filter (EnKF), are compared for the NCEP operational environment for the Global Ensemble Forecast System (GEFS). Experiments that utilize each of the two schemes are carried out and evaluated for two boreal summer seasons. It is found that these two schemes generally have comparable performance. Experiments utilizing both perturbation methods fail to generate sufficient spread at medium-range lead times beyond day 8. In general, the EnKF-based experiment outperforms the ETR in terms of the continuous ranked probability skill score (CRPSS) in the Northern Hemisphere (NH) for the first week. In the SH, the ensemble mean forecast is more skillful from the ETR perturbations. Additional experiments are performed with the stochastic total tendency perturbation (STTP) scheme, in which the total tendencies of all model variables are perturbed to represent the uncertainty in the forecast model. An improved spread–error relationship is found for the ETR-based experiments, but the STTP increases the ensemble spread for the EnKF-based experiment that is already overdispersive at early lead times, especially in the SH. With STTP employed, an increase in the EnKF-based CRPSS in the NH is reduced with a larger degradation in both the probability and ensemble-mean forecast skills in the SH. The results indicate that a rescaling of the EnKF initial perturbations and/or tuning of the STTP scheme is required when STTP is applied using the EnKF-based perturbations. This study provided guidance for the replacement of ETR with EnKF perturbations as part of the 2015 GEFS implementation.

Description: The capability of all-sky microwave radiance assimilation in the Gridpoint Statistical Interpolation (GSI) analysis system has been developed at the National Centers for Environmental Prediction (NCEP). This development effort required the adaptation of quality control, observation error assignment, bias correction, and background error covariance to all-sky conditions within the ensemble–variational (EnVar) framework. The assimilation of cloudy radiances from the Advanced Microwave Sounding Unit-A (AMSU-A) microwave radiometer for ocean fields of view (FOVs) is the primary emphasis of this study. In the original operational hybrid 3D EnVar Global Forecast System (GFS), the clear-sky approach for radiance data assimilation is applied. Changes to data thinning and quality control have allowed all-sky satellite radiances to be assimilated in the GSI. Along with the symmetric observation error assignment, additional situation-dependent observation error inflation is employed for all-sky conditions. Moreover, in addition to the current radiance bias correction, a new bias correction strategy has been applied to all-sky radiances. In this work, the static background error variance and the ensemble spread of cloud water are examined, and the levels of cloud variability from the ensemble forecast in single- and dual-resolution configurations are discussed. Overall, the all-sky approach provides more realistic simulated brightness temperatures and cloud water analysis increments, and improves analysis off the west coasts of the continents by reducing a known bias in stratus. An approximate 10% increase in the use of AMSU-A channels 1–5 and a 12% increase for channel 15 are also observed. The all-sky AMSU-A radiance assimilation became operational in the 4D EnVar GFS system upgrade of 12 May 2016.

Description: This work describes developments to improve the Doppler radial wind data assimilation scheme used in the National Centers for Environmental Prediction (NCEP) Gridpoint Statistical Interpolation (GSI) data assimilation system with a focus on convection-permitting, 0–18-h forecasts of a heavy precipitation single case study. This work focuses on two aspects: 1) the extension of the radial wind observation operator to include vertical velocity and 2) a refinement of the radial wind super-observation processing. The refinement includes reducing the magnitude of observation smoothing and allowing observations from higher scan angles into the analysis with the intent to improve the assimilation of the radar data for operational, convection-permitting models. The results of this study demonstrate that there is sensitivity to the refinement in super-observation settings. The inclusion of vertical velocity in the observation operator is shown to have a neutral to slightly positive impact on the forecast. Results from this study are suggested to be used as a foundation to prioritize future research into the effective assimilation of radial winds in an operational setting.

Description: The impact of assimilating synoptic surveillance dropwindsonde data on the analysis and forecast of the structure and intensity of Tropical Storm Karen (2013) was examined. Data-denial experiments were conducted using the NCEP hybrid 3D ensemble–variational GSI and forecasts were made using the NCEP GFS model. The assimilation of dropwindsonde data resulted in a slightly more tilted tropical cyclone vortex, stronger vertical wind shear, and more upper-tropospheric dry air west of Karen in the initial conditions. These differences grew with time in the GFS forecasts, and resulted in a weaker and more sheared vortex by 24 h in the forecast that included the dropwindsonde data. After 24 h, the cyclone reintensified in the experiment where dropwindsonde data were excluded, likely because of moist processes in a favorable region for synoptic-scale ascent ahead of a baroclinic trough. In contrast, the forecast including the dropwindsonde data kept Karen weak and also did a better job forecasting the structure and track of Karen. These results suggest that differences in the analysis and short-term evolution of Karen and the environment due to the dropwindsonde data played a role in the longer-term structure and intensity of the cyclone, including the distribution and magnitude of associated diabatic heating. These results strongly suggest that a systematic study be undertaken to examine the impact of these data on tropical cyclone structure and intensity, since previous work has focused largely on the impact on track.

Description: At the National Centers for Environmental Prediction (NCEP), a new three-dimensional variational data assimilation (3DVAR) analysis system was implemented into the operational Global Data Assimilation System (GDAS) on 1 May 2007. The new analysis system, the Gridpoint Statistical Interpolation (GSI), replaced the Spectral Statistical Interpolation (SSI) 3DVAR system, which had been operational since 1991. The GSI was developed at the Environmental Modeling Center at NCEP as part of an effort to create a more unified, robust, and efficient analysis scheme. The key aspect of the GSI is that it formulates the analysis in model grid space, which allows for more flexibility in the application of the background error covariances and makes it straightforward for a single analysis system to be used across a broad range of applications, including both global and regional modeling systems and domains. Due to the constraints of working with an operational system, the final GDAS package included many changes other than just a simple replacing of the SSI with the new GSI. The new GDAS package contained an upgrade to the Global Forecast System model, including a new vertical coordinate, as well as new features in the GSI that were never developed for the SSI. Some of these new features included changes to the observation selection, quality control, minimization algorithm, dynamic balance constraint, and assimilation of new observation types. The evaluation of the new system relative to the SSI-based system was performed for nearly an entire year of analyses and forecasts. The objective and subjective evaluations showed that the new package exhibited superior forecast performance relative to the old SSI-based system. The new system has been shown to improve forecast skill in the tropics and substantially reduce the short-term forecast error in the extratropics. This implementation has laid the groundwork for future scientific advancements in data assimilation at NCEP.

Description: The assimilation of official advisory minimum sea level pressure observations has been developed and tested in the National Centers for Environmental Prediction (NCEP) Global Data Assimilation System (GDAS) to address forecaster concerns regarding some tropical systems being far too weak in operational Global Forecast System (GFS) analyses. The assimilation of these observations has been operational in the GFS since December 2009. Using the T574 version of the NCEP GFS model, it is demonstrated that the assimilation of these observations results in a substantial reduction in the initial intensity bias for tropical systems, resulting in improved track and intensity guidance for lead times out to 5 days.

Description: A 36-h adjoint-based forecast sensitivity study of three response functions defined in the lower troposphere—average temperature in an isolated region of the upper Midwest (R1), meridional temperature difference (R2), and average zonal component of the wind (R3)—is conducted with the goal of providing a synoptic and dynamic interpretation of the sensitivity gradient structure and evolution. In addition to calculating and interpreting the sensitivity gradients with respect to basic model variables along the model forecast trajectory, a technique is outlined that allows for the calculation of the sensitivity gradients with respect to variables derivable from the model state vector (including geopotential, relative vorticity, and divergence), and a method for visualizing the sensitivities with respect to the horizontal components of the wind is proposed and demonstrated. The sensitivity of R1 to all model and derived variables revealed that R1 was controlled by nearly adiabatic processes associated with the addition or generation of temperature perturbations upstream of the region in which R1 was defined. For R2, the sensitivity gradients revealed the well-known influence of confluent horizontal flow and vertical tilting of isentropes to increase the north–south temperature gradient over the region within which R2 was defined. The sensitivity of R3 to the components of the horizontal wind reveals that simply adding or generating an upstream zonal wind perturbation is insufficient to change the zonal wind at 36 h as these wind perturbations upstream of the domain within which R3 is defined are torqued by the Coriolis force as they are advected toward the domain. These results suggest adjoint-derived sensitivities of quasi-conserved response functions may be more easily interpretable than sensitivities calculated for nonconserved response functions.

Description: The gridpoint statistical interpolation (GSI) analysis system is a unified global/regional three-dimensional variational data assimilation (3DVAR) analysis code that has been under development for several years at the National Centers for Environmental Prediction (NCEP)/Environmental Modeling Center. It has recently been implemented into operations at NCEP in both the global and North American data assimilation systems (GDAS and NDAS, respectively). An important aspect of this development has been improving the balance of the analysis produced by GSI. The improved balance between variables has been achieved through the inclusion of a tangent-linear normal-mode constraint (TLNMC). The TLNMC method has proven to be very robust and effective. The TLNMC as part of the global GSI system has resulted in substantial improvement in data assimilation at NCEP.

Description: The 24–25 January 2000 eastern United States snowstorm was noteworthy as operational numerical weather prediction (NWP) guidance was poor for lead times as short as 36 h. Despite improvements in the forecast of the surface cyclone position and intensity at 1200 UTC 25 January 2000 with decreasing lead time, NWP guidance placed the westward extent of the midtropospheric, frontogenetically forced precipitation shield too far to the east. To assess the influence of initial condition uncertainties on the forecast of this event, an adjoint model is used to evaluate forecast sensitivities for 36- and 48-h forecasts valid at 1200 UTC 25 January 2000 using as response functions the energy-weighted forecast error, lower-tropospheric circulation about a box surrounding the surface cyclone, 750-hPa frontogenesis, and vertical motion. The sensitivities with respect to the initial conditions for these response functions are in general very similar: geographically isolated, maximized in the middle and lower troposphere, and possessing an upshear vertical tilt. The sensitivities are maximized in a region of enhanced low-level baroclinicity in the vicinity of the surface cyclone’s precursor upper trough. However, differences in the phase and structure of the gradients for the four response functions are evident, which suggests that perturbations could be constructed to alter one response function but not necessarily the others. Gradients of the forecast error response function with respect to the initial conditions are used in an iterative procedure to construct initial condition perturbations that reduce the forecast error. These initial condition perturbations were small in terms of both spatial scale and magnitude. Those initial condition perturbations that were confined primarily to the midtroposphere grew rapidly into much larger amplitude upper-and-lower tropospheric perturbations. The perturbed forecasts were not only characterized by reduced final time forecast error, but also had a synoptic evolution that more closely followed analyses and observations.

Description: Experimental ensemble predictions of tropical cyclone (TC) tracks from the ensemble Kalman filter (EnKF) using the Global Forecast System (GFS) model were recently validated for the 2009 Northern Hemisphere hurricane season by Hamill et al. A similar suite of tests is described here for the 2010 season. Two major changes were made this season: 1) a reduction in the resolution of the GFS model, from 2009’s T384L64 (~31 km at 25°N) to 2010’s T254L64 (~47 km at 25°N), and some changes in model physics; and 2) the addition of a limited test of deterministic forecasts initialized from a hybrid three-dimensional variational data assimilation (3D-Var)/EnKF method. The GFS/EnKF ensembles continued to produce reduced track errors relative to operational ensemble forecasts created by the National Centers for Environmental Prediction (NCEP), the Met Office (UKMO), and the Canadian Meteorological Centre (CMC). The GFS/EnKF was not uniformly as skillful as the European Centre for Medium-Range Weather Forecasts (ECMWF) ensemble prediction system. GFS/EnKF track forecasts had slightly higher error than ECMWF at longer leads, especially in the western North Pacific, and exhibited poorer calibration between spread and error than in 2009, perhaps in part because of lower model resolution. Deterministic forecasts from the hybrid were competitive with deterministic EnKF ensemble-mean forecasts and superior in track error to those initialized from the operational variational algorithm, the Gridpoint Statistical Interpolation (GSI). Pending further successful testing, the National Oceanic and Atmospheric Administration (NOAA) intends to implement the global hybrid system operationally for data assimilation.