William A. Gallus

Toward Better Understanding of the Contiguous Rain Area (CRA) Method for Spatial Forecast Verification

The contiguous rain area (CRA) method for spatial forecast verification is a features-based approach that evaluates the properties of forecast rain systems, namely, their location, size, intensity, and finescale pattern. It is one of many recently developed spatial verification approaches that are being evaluated as part of a Spatial Forecast Verification Methods Intercomparison Project. To better understand the strengths and weaknesses of the CRA method, it has been tested here on a set of idealized geometric and perturbed forecasts with known errors, as well as nine precipitation forecasts from three high-resolution numerical weather prediction models. The CRA method was able to identify the known errors for the geometric forecasts, but only after a modification was introduced to allow nonoverlapping forecast and observed features to be matched. For the perturbed cases in which a radar rain field was spatially translated and amplified to simulate forecast errors, the CRA method alsoreproducedthe knownerrorsexcept when ahigh-intensity thresholdwas usedto define the CRA (

An East Asian Cold Surge: Case Study

10 mm h 21 ) and a large translation error was imposed (.200 km). The decomposition of total error into displacement, volume, and pattern components reflected the source of the error almost all of the time when a mean squared error formulation was used, but not necessarily when a correlation-based formulation was used. When applied to real forecasts, the CRA method gave similar results when either best-fit criteria, minimization of the mean squared error, or maximization of the correlation coefficient, was chosen for matching forecast and observed features. The diagnosed displacement error was somewhat sensitive to the choice of search distance. Of the many diagnostics produced by this method, the errors in the mean and peak rain rate between the forecast and observed features showed the best correspondence with subjective evaluations of the forecasts, while the spatial correlation coefficient (after matching) did not reflect the subjective judgments.

The Use of a Modified Ebert–McBride Technique to Evaluate Mesoscale Model QPF as a Function of Convective System Morphology during IHOP 2002

Abstract Since a cold surge is a hazardous weather phenomenon in east Asia, the rapid population increase and economic growth over the past two decades require improvement in forecasting cold surges and their related weather events over this region. However, without a better understanding of these events, this task cannot be accomplished. A cold surge with a well-defined cold front passing through Taiwan was selected to illustrate its impact on the east Asian weather system. This case is typical of a large portion of surges occurring in the region. Major findings of this study are as follows. Coupling with the upper ridge–trough structures of the wave train straddling the eastern seaboard of northeast Asia, cold surges occur sequentially. A cold front with a prefront high pressure zone is formed by the new surge outflow interacting with the anticyclone of the aging surge. The warm moist air advected northeastward along the cold front assists the development of the new surges low center, while the prefron...

Influence of Initial Conditions on the WRF–ARW Model QPF Response to Physical Parameterization Changes

Abstract The Ebert–McBride technique (EMT) is an entity-oriented method useful for quantitative precipitation verification. The EMT was modified to optimize its ability to identify contiguous rain areas (CRAs) during the 2002 International H2O Project (IHOP). This technique was then used to identify systematic sources of error as a function of observed convective system morphology in three 12-km model simulations run over the IHOP domain: Eta, the fifth-generation Pennsylvania State University–NCAR Mesoscale Model (MM5), and the Weather Research and Forecasting (WRF). The EMT was fine-tuned to optimize the pattern matching of forecasts to observations for the scales of precipitation systems observed during IHOP. To investigate several error measures provided by the EMT, a detailed morphological analysis of observed systems was performed using radar data for all CRAs identified in the IHOP domain. The modified EMT suggests that the Eta Model produced average rain rates, peak rainfall amounts, and total rai...

Spring and Summer Midwestern Severe Weather Reports in Supercells Compared to Other Morphologies

To assist in optimizing a mixed-physics ensemble for warm season mesoscale convective system rainfall forecasting, the impact of various physical schemes as well as their interactions on rainfall when different initializations were used has been investigated. For this purpose, high-resolution Weather Research and Forecasting (WRF) model simulations of eight International H2O Project events were performed. For each case, three different treatments of convection, three different microphysical schemes, and two different planetary boundary layer (PBL) schemes were used. All cases were initialized with both Local Analyses and Prediction System (LAPS) “hot” start analyses and 40-km Eta Model analyses. To evaluate the impacts of the variation of two different physical schemes and their interaction on the simulated rainfall under the two different initial conditions, the factor separation method was used. The sensitivity to the use of various physical schemes and their interactions was found to be dependent on the initialization dataset. Runs initialized with Eta analyses appeared to be influenced by the use of the Betts–Miller–Janjic scheme in that model’s assimilation system, which tended to reduce the WRF’s sensitivity to changes in the microphysical scheme compared with that present when LAPS analyses were used for initialization. In addition, differences in initialized thermodynamics resulted in changes in sensitivity to PBL and convective schemes. With both initialization datasets, the greatest sensitivity to the simulated rain rate was due to changes in the convective scheme. However, for rain volume, substantial sensitivity was present due to changes in both the physical parameterizations and the initial datasets.

On the Impact of WRF Model Vertical Grid Resolution on Midwest Summer Rainfall Forecasts

This study compares severe weather reports associated with the nine convective system morphologies used in a recent study by Gallus et al. to an additional morphology, supercell storms. As in that previous study, all convective systems occurring in a 10-state region covering parts of the Midwestern United States and central plains were classified according to their dominant morphology, and severe weather reports associated with each morphology were then analyzed. Unlike the previous study, which examined systems from 2002, the time period over which the climatology was performed was shifted to 2007 to allow access to radar algorithm information needed to classify a storm as a supercell. Archived radar imagery was used to classify systems as nonlinear convective events, isolated cells, clusters of cells, broken lines of cells, squall lines with no stratiform precipitation, trailing stratiform precipitation, parallel stratiform precipitation, and leading stratiform precipitation, and bow echoes. In addition, the three cellular classifications were subdivided to allow an analysis of severe weather reports for events in which supercells were present and those in which they were not. As in the earlier study, all morphologies were found to pose some risk of severe weather, and differences in the two datasets were generally small. The 2007 climatology confirmed the theory that supercellular systems produce severe weather more frequently than other morphologies, and also produce more intense severe weather. Supercell systems were especially prolific producers of tornadoes and hail relative to all other morphologies, but also produced severe wind and flooding much more often than nonsupercell cellular morphologies. These results suggest that it is important to differentiate between cellular morphologies containing rotation and those that do not when associating severe weather reports with convective morphology.

A WRF Ensemble for Improved Wind Speed Forecasts at Turbine Height

Weather Research and Forecast (WRF) model exploratory sensitivity simulations were performed to determine the impact of vertical grid resolution (VGR) on the forecast skill of Midwest summer rainfall. Varying the VGR indicated that a refined VGR, while adopting the widely used North America Regional Reanalysis (NARR) for initial and lateral boundary conditions, does not necessarily result in a consistent improvement in quantitative precipitation forecasts (QPFs). When averaged over a variety of microphysical schemes in an illustrative case, equitable threat score (ETS) and bias values actually worsened with a greater overpredicted rainfall for half of the rainfall thresholds when the VGR was refined. Averaged over strongly forced cases, ETS values worsened for all rainfall thresholds while biases mostly increased, indicating a further overprediction of rainfall when the number of levels was increased. Skill improved, however, for all rainfall thresholds when the resolution above the melting level was increased. Skill also improved for most rainfall thresholds when the resolution in the surface layer was increased, which is attributed to better-resolved surface turbulent momentum and thermal fluxes. Likewise, a refined VGR resulted in improvements in weakly forced cases, which are governed mostly by thermodynamic forcing and are sensitive to vertical profiles of temperature and moisture. Application of the factor separation method suggested that the refined VGR more frequently had a negative impact on skill through the interaction between lower-atmospheric processes and microphysical processes above the melting level.

Role of stratospheric air in a severe weather event: Analysis of potential vorticity and total ozone

The Weather Research and Forecasting Model (WRF) with 10-km horizontal grid spacing was used to explore improvements in wind speed forecasts at a typical wind turbine hub height (80 m). An ensemble consisting of WRF model simulations with different planetary boundary layer (PBL) schemes showed little spread among the individual ensemble members for forecasting wind speed. A second configuration using three random perturbations of the Global Forecast System model produced more spread in the wind speed forecasts, but the ensemble mean possessed a higher mean absolute error (MAE). A third ensemble of different initialization times showed larger model spread, but model MAE was not compromised. In addition, postprocessing techniques such as training of the model for the day 2 forecast based on day 1 results and bias correction based on observed wind direction are examined. Ramp event forecasting was also explored. An event was considered to be a ramp event if the change in wind power was 50% or more of total capacity in either 4 or 2 h or less. This was approximated using a typical wind turbine power curve such that any wind speed increase or decrease of more than 3 m s 21 within the 6‐12 m s 21 window (where power production variesgreatly)in4 horlesswouldbeconsideredaramp.ModelMAE,climatologyoframpevents,andcauses were examined. All PBL schemes examined predicted fewer ramp events compared to the observations, and model forecasts for ramps in general were poor.

The 2015 Plains Elevated Convection at Night Field Project

The role of dry stratospheric air descending to low and middle tropospheric levels in a severe weather outbreak in the midwestern United States is examined using NCEP Eta model output, Rapid Update Cycle (RUC) analyses, and Earth probe Total Ozone Mapping Spectrometer (EP/TOMS) total ozone data. While stratospheric air was not found to play a direct role in the convection, backward trajectories show stratospheric air descended to 800 hPa just west of the convection. Damaging surface winds not associated with thunderstorms also occurred in the region of greatest stratospheric descent. Small-scale features in the high-resolution total ozone data compare favorably with geopotential heights and potential vorticity fields, supporting the notion that stratospheric air descended to near the surface. A detailed vertical structure in the potential vorticity appears to be captured by small-scale total ozone variations. The capability of the total ozone to identify mesoscale features assists model verification. The total ozone data suggest biases in the RUC analysis and Eta forecast of this event. The total ozone is also useful in determining whether potential vorticity is of stratospheric origin or is diabatically generated in the troposphere.

The Impact of Large-Scale Forcing on Skill of Simulated Convective Initiation and Upscale Evolution with Convection-Allowing Grid Spacings in the WRF*

AbstractThe central Great Plains region in North America has a nocturnal maximum in warm-season precipitation. Much of this precipitation comes from organized mesoscale convective systems (MCSs). This nocturnal maximum is counterintuitive in the sense that convective activity over the Great Plains is out of phase with the local generation of CAPE by solar heating of the surface. The lower troposphere in this nocturnal environment is typically characterized by a low-level jet (LLJ) just above a stable boundary layer (SBL), and convective available potential energy (CAPE) values that peak above the SBL, resulting in convection that may be elevated, with source air decoupled from the surface. Nocturnal MCS-induced cold pools often trigger undular bores and solitary waves within the SBL. A full understanding of the nocturnal precipitation maximum remains elusive, although it appears that bore-induced lifting and the LLJ may be instrumental to convection initiation and the maintenance of MCSs at night.To gain ...