An ensemble Kalman filter–variational hybrid data assimilation system based on the gridpoint statistical interpolation (GSI) three-dimensional variational data assimilation (3DVar) system was developed. The performance of the system was investigated using the National Centers for Environmental Prediction (NCEP) Global Forecast System model. Experiments covered a 6-week Northern Hemisphere winter period. Both the control and ensemble forecasts were run at the same, reduced resolution. Operational conventional and satellite observations along with an 80-member ensemble were used. Various configurations of the system including one- or two-way couplings, with zero or nonzero weights on the static covariance, were intercompared and compared with the GSI 3DVar system. It was found that the hybrid system produced more skillful forecasts than the GSI 3DVar system. The inclusion of a static component in the background-error covariance and recentering the analysis ensemble around the variational analysis did not improve the forecast skill beyond the one-way coupled system with zero weights on the static covariance. The one-way coupled system with zero static covariances produced more skillful wind forecasts averaged over the globe than the EnKF at the 1–5-day lead times and more skillful temperature forecasts than the EnKF at the 5-day lead time. Sensitivity tests indicated that the difference may be due to the use of the tangent linear normal mode constraint in the variational system. For the first outer loop, the hybrid system showed a slightly slower (faster) convergence rate at early (later) iterations than the GSI 3DVar system. For the second outer loop, the hybrid system showed a faster convergence.

A GSI-based EnVar data assimilation system is extended to directly assimilate radar reflectivity to initialize convective-scale forecasts. When hydrometeor mixing ratios are used as state variables (method mixing ratio), large differences of the cost function gradients with respect to the small hydrometeor mixing ratios and wind prevent efficient convergence. Using logarithmic mixing ratios as state variables (method logarithm) fixes this problem, but generates spuriously large hydrometeor increments partly due to the transform to and from the logarithmic space. The tangent linear of the reflectivity operators further contributes to spuriously small and large hydrometeor increments in method mixing ratio and method logarithm, respectively. A new method is proposed by directly adding the reflectivity as a state variable (method dBZ). Without the tangent linear and adjoint of the nonlinear operator, the new method therefore avoids the aforementioned problems. The newly proposed method is examined on the analysis and prediction of the 8 May 2003 Oklahoma City tornadic supercell storm. Both the probabilistic forecast of strong low-level vorticity and maintenance of strong updraft and vorticity in method dBZ are more consistent with reality than in method logarithm and method mixing ratio. Detailed diagnostics suggest that a more realistic cold pool due to the better analyzed hydrometeors in method dBZ than in other methods leads to constructive interaction between the surface gust front and the updraft aloft associated with the midlevel mesocyclone. Similar low-level vorticity forecast and maintenance of the storm are produced by the WSM6 and Thompson microphysics schemes in method dBZ. The Thompson scheme matches the reflectivity distribution with the observations better for all lead times, but shows more southeastward track bias compared to the WSM6 scheme.

Forecasts generated by the Center for Analysis and Prediction of Storms with 1- and 4-km grid spacing using the Advanced Research Weather Research and Forecasting Model (ARW-WRF; ARW1 and ARW4, respectively) for the 2009–11 NOAA Hazardous Weather Testbed Spring Experiments are compared and verified. Object-based measures, including average values of object attributes, the object-based threat score (OTS), and the median of maximum interest (MMI) are used for the verification. Verification was first performed against observations at scales resolvable by each forecast model and then performed at scales resolvable by both models by remapping ARW1 to the ARW4 grid (ARW1to4). Thirty-hour forecasts of 1-h accumulated precipitation initialized at 0000 UTC on 22, 36, and 33 days during the spring of 2009, 2010, and 2011, respectively, are evaluated over a domain covering most of the central and eastern United States. ARW1, ARW1to4, and ARW4 all significantly overforecasted the number of objects during diurnal convection maxima. The overforecasts by ARW1 and ARW1to4 were more pronounced than ARW4 during the first convection maximum at 1-h lead time. The average object area and aspect ratio were closer to observations for ARW1 and ARW1to4 than for ARW4. None of the models showed a significant advantage over the others for average orientation angle and centroid location. Increased accuracy for ARW1, compared to ARW4, was statistically significant for the MMI but not the OTS. However, ARW1to4 had similar MMI and OTS as ARW4 at most lead times. These results are consistent with subjective evaluations that the greatest impact of grid spacing is on the smallest resolvable objects.

A hybrid ensemble transform Kalman filter (ETKF)–three-dimensional variational data assimilation (3DVAR) system developed for the Weather Research and Forecasting Model (WRF) was studied for the forecasts of the tracks of two major hurricanes, Ike and Gustav, in 2008 over the Gulf of Mexico. The impacts of the flow-dependent ensemble covariance generated by the ETKF were revealed by comparing the forecasts, analyses, and analysis increments generated by the hybrid data assimilation method with those generated by the 3DVAR that used the static background covariance. The root-mean-square errors of the track forecasts by the hybrid data assimilation (DA) method were smaller than those by the 3DVAR for both Ike and Gustav. Experiments showed that such improvements were due to the use of the flow-dependent covariance provided by the ETKF ensemble in the hybrid DA system. Detailed diagnostics further revealed that the increments produced by the hybrid and the 3DVAR were different for both the analyses of the hurricane itself and its environment. In particular, it was found that the hybrid, using the flow-dependent covariance that gave the hurricane-specific error covariance estimates, was able to systematically adjust the position of the hurricane during the assimilation whereas the 3DVAR was not. The study served as a pilot study to explore and understand the potential of the hybrid method for hurricane data assimilation and forecasts. Caution needs to be taken to extrapolate the results to operational forecast settings.

Neighborhood and object-based probabilistic precipitation forecasts from a convection-allowing ensemble are verified and calibrated. Calibration methods include logistic regression, one- and two-parameter reliability-based calibration, and cumulative distribution function (CDF)-based bias adjustment. Newly proposed object-based probabilistic forecasts for the occurrence of a forecast object are derived from the percentage of ensemble members with a matching object. Verification and calibration of single- and multimodel subensembles are performed to explore the effect of using multiple models. The uncalibrated neighborhood-based probabilistic forecasts have skill minima during the afternoon convective maximum. Calibration generally improves the skill, especially during the skill minima, resulting in positive skill. In general all calibration methods perform similarly, with a slight advantage of logistic regression (one-parameter reliability based) calibration for 1-h (6 h) accumulations. The uncalibrated object-based probabilistic forecasts are, in general, less skillful than the uncalibrated neighborhood-based probabilistic forecasts. Object-based calibration also results in positive skill at all lead times. For object-based calibration the skill is significantly different among the calibration methods, with the logistic regression performing the best and CDF-based bias adjustment performing the worst. For both the neighborhood and object-based probabilistic forecasts, the impact of using 10 or 25 days of training data for calibration is generally small and is most significant for the two-parameter reliability-based method. An uncalibrated Advanced Research Weather Research and Forecasting Model (ARW-WRF) subensemble is significantly more skillful than an uncalibrated WRF Nonhydrostatic Mesoscale Model (NMM) subensemble. The difference is reduced by calibration. The multimodel subensemble only shows an advantage for the neighborhood-based forecasts beyond 1-day lead time and shows no advantage for the object-based forecasts.

Twenty-member real-time convection-allowing storm-scale ensemble forecasts with perturbations to model physics, dynamics, initial conditions (IC), and lateral boundary conditions (LBC) during the NOAA Hazardous Weather Testbed Spring Experiment provide a unique opportunity to study the relative impact of different sources of perturbation on convection-allowing ensemble diversity. In Part II of this two-part study, systematic similarity/dissimilarity of hourly precipitation forecasts among ensemble members from the spring season of 2009 are identified using hierarchical cluster analysis (HCA) with a fuzzy object-based threat score (OTS), developed in Part I. In addition to precipitation, HCA is also performed on ensemble forecasts using the traditional Euclidean distance for wind speed at 10 m and 850 hPa, and temperature at 500 hPa. At early lead times (3 h, valid at 0300 UTC) precipitation forecasts cluster primarily by data assimilation and model dynamic core, indicating a dominating impact of models, with secondary clustering by microphysics. There is an increasing impact of the planetary boundary layer (PBL) scheme on clustering relative to the microphysics scheme at later lead times. Forecasts of 10-m wind speed cluster primarily by the PBL scheme at early lead times, with an increasing impact of LBC at later lead times. Forecasts of midtropospheric variables cluster primarily by IC at early lead times and LBC at later lead times. The radar and Mesonet data assimilation (DA) show its impact, with members without DA in a distinct cluster, through the 12-h lead time (valid at 1200 UTC) for both precipitation and nonprecipitation variables. The implication for optimal ensemble design for storm-scale forecasts is also discussed.