Search Results
You are looking at 51 - 60 of 78 items for
- Author or Editor: Harold E. Brooks x
- Refine by Access: All Content x
Abstract
Two novel approaches to extending the range of prediction for environments conducive to severe thunderstorm events are described. One approach charts Climate Forecast System, version 2 (CFSv2), run-to-run consistency of the areal extent of severe thunderstorm environments using grid counts of the supercell composite parameter (SCP). Visualization of these environments is charted for each 45-day CFSv2 run initialized at 0000 UTC. CFSv2 ensemble-mean forecast maps of SCP coverage over the contiguous United States are also produced for those forecasts meeting certain criteria for high-impact weather. The applicability of this approach to the severe weather prediction challenge is illustrated using CFSv2 output for a series of severe weather episodes occurring in March and April 2014. Another approach, possibly extending severe weather predictability from CFSv2, utilizes a run-cumulative time-averaging technique of SCP grid counts. This process is described and subjectively verified with severe weather events from early 2014.
Abstract
Two novel approaches to extending the range of prediction for environments conducive to severe thunderstorm events are described. One approach charts Climate Forecast System, version 2 (CFSv2), run-to-run consistency of the areal extent of severe thunderstorm environments using grid counts of the supercell composite parameter (SCP). Visualization of these environments is charted for each 45-day CFSv2 run initialized at 0000 UTC. CFSv2 ensemble-mean forecast maps of SCP coverage over the contiguous United States are also produced for those forecasts meeting certain criteria for high-impact weather. The applicability of this approach to the severe weather prediction challenge is illustrated using CFSv2 output for a series of severe weather episodes occurring in March and April 2014. Another approach, possibly extending severe weather predictability from CFSv2, utilizes a run-cumulative time-averaging technique of SCP grid counts. This process is described and subjectively verified with severe weather events from early 2014.
Abstract
Numerical forecasts from a pilot program on short-range ensemble forecasting at the National Centers for Environmental Prediction are examined. The ensemble consists of 10 forecasts made using the 80-km Eta Model and 5 forecasts from the regional spectral model. Results indicate that the accuracy of the ensemble mean is comparable to that from the 29-km Meso Eta Model for both mandatory level data and the 36-h forecast cyclone position. Calculations of spread indicate that at 36 and 48 h the spread from initial conditions created using the breeding of growing modes technique is larger than the spread from initial conditions created using different analyses. However, the accuracy of the forecast cyclone position from these two initialization techniques is nearly identical. Results further indicate that using two different numerical models assists in increasing the ensemble spread significantly.
There is little correlation between the spread in the ensemble members and the accuracy of the ensemble mean for the prediction of cyclone location. Since information on forecast uncertainty is needed in many applications, and is one of the reasons to use an ensemble approach, the lack of a correlation between spread and forecast uncertainty presents a challenge to the production of short-range ensemble forecasts.
Even though the ensemble dispersion is not found to be an indication of forecast uncertainty, significant spread can occur within the forecasts over a relatively short time period. Examples are shown to illustrate how small uncertainties in the model initial conditions can lead to large differences in numerical forecasts from an identical numerical model.
Abstract
Numerical forecasts from a pilot program on short-range ensemble forecasting at the National Centers for Environmental Prediction are examined. The ensemble consists of 10 forecasts made using the 80-km Eta Model and 5 forecasts from the regional spectral model. Results indicate that the accuracy of the ensemble mean is comparable to that from the 29-km Meso Eta Model for both mandatory level data and the 36-h forecast cyclone position. Calculations of spread indicate that at 36 and 48 h the spread from initial conditions created using the breeding of growing modes technique is larger than the spread from initial conditions created using different analyses. However, the accuracy of the forecast cyclone position from these two initialization techniques is nearly identical. Results further indicate that using two different numerical models assists in increasing the ensemble spread significantly.
There is little correlation between the spread in the ensemble members and the accuracy of the ensemble mean for the prediction of cyclone location. Since information on forecast uncertainty is needed in many applications, and is one of the reasons to use an ensemble approach, the lack of a correlation between spread and forecast uncertainty presents a challenge to the production of short-range ensemble forecasts.
Even though the ensemble dispersion is not found to be an indication of forecast uncertainty, significant spread can occur within the forecasts over a relatively short time period. Examples are shown to illustrate how small uncertainties in the model initial conditions can lead to large differences in numerical forecasts from an identical numerical model.
Abstract
Many tornadoes are unreported because of lack of observers or are underrated in intensity, width, or track length because of lack of damage indicators. These reporting biases substantially degrade estimates of tornado frequency and thereby undermine important endeavors such as studies of climate impacts on tornadoes and cost–benefit analyses of tornado damage mitigation. Building on previous studies, we use a Bayesian hierarchical modeling framework to estimate and correct for tornado reporting biases over the central United States during 1975–2018. The reporting biases are treated as a univariate function of population density. We assess how these biases vary with tornado intensity, width, and track length and over the analysis period. We find that the frequencies of tornadoes of all kinds, but especially stronger or wider tornadoes, have been substantially underestimated. Most strikingly, the Bayesian model estimates that there have been approximately 3 times as many tornadoes capable of (E)F2+ damage as have been recorded as (E)F2+ [(E)F indicates a rating on the (enhanced) Fujita scale]. The model estimates that total tornado frequency changed little over the analysis period. Statistically significant trends in frequency are found for tornadoes within certain ranges of intensity, pathlength, and width, but it is unclear what proportion of these trends arise from changes in damage survey practices. Simple analyses of the tornado database corroborate many of the inferences from the Bayesian model.
Significance Statement
Prior studies have shown that the probabilities of a tornado being reported and of its intensity, track length, and width being accurately estimated are strongly correlated with the local population density. We have developed a sophisticated statistical model that accounts for these population-dependent tornado reporting biases to improve estimates of tornado frequency in the central United States. The bias-corrected tornado frequency estimates differ markedly from the official tornado climatology and have important implications for tornado risk assessment, damage mitigation, and studies of climate change impacts on tornado activity.
Abstract
Many tornadoes are unreported because of lack of observers or are underrated in intensity, width, or track length because of lack of damage indicators. These reporting biases substantially degrade estimates of tornado frequency and thereby undermine important endeavors such as studies of climate impacts on tornadoes and cost–benefit analyses of tornado damage mitigation. Building on previous studies, we use a Bayesian hierarchical modeling framework to estimate and correct for tornado reporting biases over the central United States during 1975–2018. The reporting biases are treated as a univariate function of population density. We assess how these biases vary with tornado intensity, width, and track length and over the analysis period. We find that the frequencies of tornadoes of all kinds, but especially stronger or wider tornadoes, have been substantially underestimated. Most strikingly, the Bayesian model estimates that there have been approximately 3 times as many tornadoes capable of (E)F2+ damage as have been recorded as (E)F2+ [(E)F indicates a rating on the (enhanced) Fujita scale]. The model estimates that total tornado frequency changed little over the analysis period. Statistically significant trends in frequency are found for tornadoes within certain ranges of intensity, pathlength, and width, but it is unclear what proportion of these trends arise from changes in damage survey practices. Simple analyses of the tornado database corroborate many of the inferences from the Bayesian model.
Significance Statement
Prior studies have shown that the probabilities of a tornado being reported and of its intensity, track length, and width being accurately estimated are strongly correlated with the local population density. We have developed a sophisticated statistical model that accounts for these population-dependent tornado reporting biases to improve estimates of tornado frequency in the central United States. The bias-corrected tornado frequency estimates differ markedly from the official tornado climatology and have important implications for tornado risk assessment, damage mitigation, and studies of climate change impacts on tornado activity.
Abstract
Forecasts from the National Centers for Environmental Prediction’s experimental short-range ensemble system are examined and compared with a single run from a higher-resolution model using similar computational resources. The ensemble consists of five members from the Regional Spectral Model and 10 members from the 80-km Eta Model, with both in-house analyses and bred perturbations used as initial conditions. This configuration allows for a comparison of the two models and the two perturbation strategies, as well as a preliminary investigation of the relative merits of mixed-model, mixed-perturbation ensemble systems. The ensemble is also used to estimate the short-range predictability limits of forecasts of precipitation and fields relevant to the forecast of precipitation.
Whereas error growth curves for the ensemble and its subgroups are in relative agreement with previous work for large-scale fields such as 500-mb heights, little or no error growth is found for fields of mesoscale interest, such as convective indices and precipitation. The difference in growth rates among the ensemble subgroups illustrates the role of both initial perturbation strategy and model formulation in creating ensemble dispersion. However, increase spread per se is not necessarily beneficial, as is indicated by the fact that the ensemble subgroup with the greatest spread is less skillful than the subgroup with the least spread.
Further examination into the skill of the ensemble system for forecasts of precipitation shows the advantage gained from a mixed-model strategy, such that even the inclusion of the less skillful Regional Spectral Model members improves ensemble performance. For some aspects of forecast performance, even ensemble configurations with as few as five members are shown to significantly outperform the 29-km Meso-Eta Model.
Abstract
Forecasts from the National Centers for Environmental Prediction’s experimental short-range ensemble system are examined and compared with a single run from a higher-resolution model using similar computational resources. The ensemble consists of five members from the Regional Spectral Model and 10 members from the 80-km Eta Model, with both in-house analyses and bred perturbations used as initial conditions. This configuration allows for a comparison of the two models and the two perturbation strategies, as well as a preliminary investigation of the relative merits of mixed-model, mixed-perturbation ensemble systems. The ensemble is also used to estimate the short-range predictability limits of forecasts of precipitation and fields relevant to the forecast of precipitation.
Whereas error growth curves for the ensemble and its subgroups are in relative agreement with previous work for large-scale fields such as 500-mb heights, little or no error growth is found for fields of mesoscale interest, such as convective indices and precipitation. The difference in growth rates among the ensemble subgroups illustrates the role of both initial perturbation strategy and model formulation in creating ensemble dispersion. However, increase spread per se is not necessarily beneficial, as is indicated by the fact that the ensemble subgroup with the greatest spread is less skillful than the subgroup with the least spread.
Further examination into the skill of the ensemble system for forecasts of precipitation shows the advantage gained from a mixed-model strategy, such that even the inclusion of the less skillful Regional Spectral Model members improves ensemble performance. For some aspects of forecast performance, even ensemble configurations with as few as five members are shown to significantly outperform the 29-km Meso-Eta Model.
Abstract
Over the last 50 yr, the number of tornadoes reported in the United States has doubled from about 600 per year in the 1950s to around 1200 in the 2000s. This doubling is likely not related to meteorological causes alone. To account for this increase a simple least squares linear regression was fitted to the annual number of tornado reports. A “big tornado day” is a single day when numerous tornadoes and/or many tornadoes exceeding a specified intensity threshold were reported anywhere in the country. By defining a big tornado day without considering the spatial distribution of the tornadoes, a big tornado day differs from previous definitions of outbreaks. To address the increase in the number of reports, the number of reports is compared to the expected number of reports in a year based on linear regression. In addition, the F1 and greater Fujita-scale record was used in determining a big tornado day because the F1 and greater series was more stationary over time as opposed to the F2 and greater series. Thresholds were applied to the data to determine the number and intensities of the tornadoes needed to be considered a big tornado day. Possible threshold values included fractions of the annual expected value associated with the linear regression and fixed numbers for the intensity criterion. Threshold values of 1.5% of the expected annual total number of tornadoes and/or at least 8 F1 and greater tornadoes identified about 18.1 big tornado days per year. Higher thresholds such as 2.5% and/or at least 15 F1 and greater tornadoes showed similar characteristics, yet identified approximately 6.2 big tornado days per year. Finally, probability distribution curves generated using kernel density estimation revealed that big tornado days were more likely to occur slightly earlier in the year and have a narrower distribution than any given tornado day.
Abstract
Over the last 50 yr, the number of tornadoes reported in the United States has doubled from about 600 per year in the 1950s to around 1200 in the 2000s. This doubling is likely not related to meteorological causes alone. To account for this increase a simple least squares linear regression was fitted to the annual number of tornado reports. A “big tornado day” is a single day when numerous tornadoes and/or many tornadoes exceeding a specified intensity threshold were reported anywhere in the country. By defining a big tornado day without considering the spatial distribution of the tornadoes, a big tornado day differs from previous definitions of outbreaks. To address the increase in the number of reports, the number of reports is compared to the expected number of reports in a year based on linear regression. In addition, the F1 and greater Fujita-scale record was used in determining a big tornado day because the F1 and greater series was more stationary over time as opposed to the F2 and greater series. Thresholds were applied to the data to determine the number and intensities of the tornadoes needed to be considered a big tornado day. Possible threshold values included fractions of the annual expected value associated with the linear regression and fixed numbers for the intensity criterion. Threshold values of 1.5% of the expected annual total number of tornadoes and/or at least 8 F1 and greater tornadoes identified about 18.1 big tornado days per year. Higher thresholds such as 2.5% and/or at least 15 F1 and greater tornadoes showed similar characteristics, yet identified approximately 6.2 big tornado days per year. Finally, probability distribution curves generated using kernel density estimation revealed that big tornado days were more likely to occur slightly earlier in the year and have a narrower distribution than any given tornado day.
Abstract
The primary objective of this study was to estimate the percentage of U.S. tornadoes that are spawned annually by squall lines and bow echoes, or quasi-linear convective systems (QLCSs). This was achieved by examining radar reflectivity images for every tornado event recorded during 1998–2000 in the contiguous United States. Based on these images, the type of storm associated with each tornado was classified as cell, QLCS, or other.
Of the 3828 tornadoes in the database, 79% were produced by cells, 18% were produced by QLCSs, and the remaining 3% were produced by other storm types, primarily rainbands of landfallen tropical cyclones. Geographically, these percentages as well as those based on tornado days exhibited wide variations. For example, 50% of the tornado days in Indiana were associated with QLCSs.
In an examination of other tornado attributes, statistically more weak (F1) and fewer strong (F2–F3) tornadoes were associated with QLCSs than with cells. QLCS tornadoes were more probable during the winter months than were cells. And finally, QLCS tornadoes displayed a comparatively higher and statistically significant tendency to occur during the late night/early morning hours. Further analysis revealed a disproportional decrease in F0–F1 events during this time of day, which led the authors to propose that many (perhaps as many as 12% of the total) weak QLCSs tornadoes were not reported.
Abstract
The primary objective of this study was to estimate the percentage of U.S. tornadoes that are spawned annually by squall lines and bow echoes, or quasi-linear convective systems (QLCSs). This was achieved by examining radar reflectivity images for every tornado event recorded during 1998–2000 in the contiguous United States. Based on these images, the type of storm associated with each tornado was classified as cell, QLCS, or other.
Of the 3828 tornadoes in the database, 79% were produced by cells, 18% were produced by QLCSs, and the remaining 3% were produced by other storm types, primarily rainbands of landfallen tropical cyclones. Geographically, these percentages as well as those based on tornado days exhibited wide variations. For example, 50% of the tornado days in Indiana were associated with QLCSs.
In an examination of other tornado attributes, statistically more weak (F1) and fewer strong (F2–F3) tornadoes were associated with QLCSs than with cells. QLCS tornadoes were more probable during the winter months than were cells. And finally, QLCS tornadoes displayed a comparatively higher and statistically significant tendency to occur during the late night/early morning hours. Further analysis revealed a disproportional decrease in F0–F1 events during this time of day, which led the authors to propose that many (perhaps as many as 12% of the total) weak QLCSs tornadoes were not reported.
Abstract
We compare over 1 million sounding measurements with ERA-Interim reanalysis for the 38-yr period from 1979 to 2016. The large dataset allows us to provide an improved insight into the spatial and temporal distributions of the prerequisites of deep moist convection across Europe. In addition, ERA-Interim is also evaluated. ERA-Interim estimates parameters describing boundary layer moisture and midtropospheric lapse rates well, with correlation coefficients of 0.94. Mixed-layer CAPE is, on average, underestimated by the reanalysis while the most unstable CAPE is overestimated. Vertical shear parameters in the reanalysis are better correlated with observations than CAPE, but are underestimated by approximately 1–2 m s−1. The underestimation decreases as the depth of the shear layer increases. Compared to radiosonde observations, instability in ERA-Interim is overestimated in southern Europe and underestimated over eastern Europe. High values of instability are observed from May to September, out of phase with the climatological pattern of wind shear, which peaks in winter. From September to April, favorable conditions for thunderstorms occur mainly over southern and western Europe with the peak location and higher frequency shifting to central and eastern Europe from May to August. For southern Europe, the annual cycle peaks in September with high values of inhibition suppressing thunderstorm activity in July and August. The area with the highest annual number of days with environmental conditions favorable for thunderstorms extends from Italy and Austria eastward through the Carpathians and Balkans.
Abstract
We compare over 1 million sounding measurements with ERA-Interim reanalysis for the 38-yr period from 1979 to 2016. The large dataset allows us to provide an improved insight into the spatial and temporal distributions of the prerequisites of deep moist convection across Europe. In addition, ERA-Interim is also evaluated. ERA-Interim estimates parameters describing boundary layer moisture and midtropospheric lapse rates well, with correlation coefficients of 0.94. Mixed-layer CAPE is, on average, underestimated by the reanalysis while the most unstable CAPE is overestimated. Vertical shear parameters in the reanalysis are better correlated with observations than CAPE, but are underestimated by approximately 1–2 m s−1. The underestimation decreases as the depth of the shear layer increases. Compared to radiosonde observations, instability in ERA-Interim is overestimated in southern Europe and underestimated over eastern Europe. High values of instability are observed from May to September, out of phase with the climatological pattern of wind shear, which peaks in winter. From September to April, favorable conditions for thunderstorms occur mainly over southern and western Europe with the peak location and higher frequency shifting to central and eastern Europe from May to August. For southern Europe, the annual cycle peaks in September with high values of inhibition suppressing thunderstorm activity in July and August. The area with the highest annual number of days with environmental conditions favorable for thunderstorms extends from Italy and Austria eastward through the Carpathians and Balkans.
Abstract
The problem of forecasting the maintenance of mesoscale convective systems (MCSs) is investigated through an examination of observed proximity soundings. Furthermore, environmental variables that are statistically different between mature and weakening MCSs are input into a logistic regression procedure to develop probabilistic guidance on MCS maintenance, focusing on warm-season quasi-linear systems that persist for several hours. Between the mature and weakening MCSs, shear vector magnitudes over very deep layers are the best discriminators among hundreds of kinematic and thermodynamic variables. An analysis of the shear profiles reveals that the shear component perpendicular to MCS motion (usually parallel to the leading line) accounts for much of this difference in low levels and the shear component parallel to MCS motion accounts for much of this difference in mid- to upper levels. The lapse rates over a significant portion of the convective cloud layer, the convective available potential energy, and the deep-layer mean wind speed are also very good discriminators and collectively provide a high level of discrimination between the mature and dissipation soundings as revealed by linear discriminant analysis. Probabilistic equations developed from these variables used with short-term numerical model output show utility in forecasting the transition of an MCS with a solid line of 50+ dBZ echoes to a more disorganized system with unsteady changes in structure and propagation. This study shows that empirical forecast tools based on environmental relationships still have the potential to provide forecasters with improved information on the qualitative characteristics of MCS structure and longevity. This is especially important since the current and near-term value added by explicit numerical forecasts of convection is still uncertain.
Abstract
The problem of forecasting the maintenance of mesoscale convective systems (MCSs) is investigated through an examination of observed proximity soundings. Furthermore, environmental variables that are statistically different between mature and weakening MCSs are input into a logistic regression procedure to develop probabilistic guidance on MCS maintenance, focusing on warm-season quasi-linear systems that persist for several hours. Between the mature and weakening MCSs, shear vector magnitudes over very deep layers are the best discriminators among hundreds of kinematic and thermodynamic variables. An analysis of the shear profiles reveals that the shear component perpendicular to MCS motion (usually parallel to the leading line) accounts for much of this difference in low levels and the shear component parallel to MCS motion accounts for much of this difference in mid- to upper levels. The lapse rates over a significant portion of the convective cloud layer, the convective available potential energy, and the deep-layer mean wind speed are also very good discriminators and collectively provide a high level of discrimination between the mature and dissipation soundings as revealed by linear discriminant analysis. Probabilistic equations developed from these variables used with short-term numerical model output show utility in forecasting the transition of an MCS with a solid line of 50+ dBZ echoes to a more disorganized system with unsteady changes in structure and propagation. This study shows that empirical forecast tools based on environmental relationships still have the potential to provide forecasters with improved information on the qualitative characteristics of MCS structure and longevity. This is especially important since the current and near-term value added by explicit numerical forecasts of convection is still uncertain.