Browse
Abstract
In operational weather forecasting, it is effective to aggregate information on all members of an ensemble forecast through the cluster analysis. The temporal coherence of ensemble members in each cluster is an important piece of information whether each cluster is well-divided. This information is especially important for forecasts where the target area is city or prefecture, i.e., Eulerian framework, because the members that compose the clusters can change over time due to the small target region. This study provided the temporal coherence of members in clusters by performing principal component analysis and cluster analysis on 3-hourly 500-hPa geopotential height forecasts and linking the clustering results in the time direction. The new method provided a consistently well-divided forecast scenario throughout the forecast period for Eulerian frame forecasts, as well as information on the temporal coherency of the members in the clusters, which was demonstrated effective through the experiment to the pre-select cluster with small errors. The application of the new technique to improve precipitation forecasts was also discussed.
Abstract
In operational weather forecasting, it is effective to aggregate information on all members of an ensemble forecast through the cluster analysis. The temporal coherence of ensemble members in each cluster is an important piece of information whether each cluster is well-divided. This information is especially important for forecasts where the target area is city or prefecture, i.e., Eulerian framework, because the members that compose the clusters can change over time due to the small target region. This study provided the temporal coherence of members in clusters by performing principal component analysis and cluster analysis on 3-hourly 500-hPa geopotential height forecasts and linking the clustering results in the time direction. The new method provided a consistently well-divided forecast scenario throughout the forecast period for Eulerian frame forecasts, as well as information on the temporal coherency of the members in the clusters, which was demonstrated effective through the experiment to the pre-select cluster with small errors. The application of the new technique to improve precipitation forecasts was also discussed.
Abstract
A hybrid three-dimensional ensemble–variational (En3D-Var) data assimilation system has been developed to explore incorporating information from an 11-member regional ensemble prediction system, which is dynamically downscaled from a global ensemble system, into a 3-hourly cycling convective-scale data assimilation system over the western Maritime Continent. From the ensemble, there exists small-scale ensemble perturbation structures associated with positional differences of tropical convection, but these structures are well represented only after the downscaled ensemble forecast has evolved for at least 6 h due to spinup. There was also a robust moderate negative correlation between total specific humidity and potential temperature background errors, presumably because of incorrect vertical motion in the presence of clouds. Time shifting of the ensemble perturbations, by using those available from adjacent cycles, helped to ameliorate the sampling error prevalent in their raw autocovariances. Monthlong hybrid En3D-Var trials were conducted using different weights assigned to the ensemble-derived and climatological background error covariances. The forecast fits to radiosonde relative humidity and wind observations were generally improved with hybrid En3D-Var, but in all experiments, the fits to surface observations were degraded compared to the baseline 3D-Var configuration. Over the Singapore radar domain, there was a general improvement in the precipitation forecasts, especially when the weighting toward the climatological background error covariance was larger, and with the additional application of time-shifted ensemble perturbations. Future work involves consolidating the ensemble prediction and deterministic system, by centering the ensemble prediction system on the hybrid analysis, to better represent the analysis and forecast uncertainties.
Abstract
A hybrid three-dimensional ensemble–variational (En3D-Var) data assimilation system has been developed to explore incorporating information from an 11-member regional ensemble prediction system, which is dynamically downscaled from a global ensemble system, into a 3-hourly cycling convective-scale data assimilation system over the western Maritime Continent. From the ensemble, there exists small-scale ensemble perturbation structures associated with positional differences of tropical convection, but these structures are well represented only after the downscaled ensemble forecast has evolved for at least 6 h due to spinup. There was also a robust moderate negative correlation between total specific humidity and potential temperature background errors, presumably because of incorrect vertical motion in the presence of clouds. Time shifting of the ensemble perturbations, by using those available from adjacent cycles, helped to ameliorate the sampling error prevalent in their raw autocovariances. Monthlong hybrid En3D-Var trials were conducted using different weights assigned to the ensemble-derived and climatological background error covariances. The forecast fits to radiosonde relative humidity and wind observations were generally improved with hybrid En3D-Var, but in all experiments, the fits to surface observations were degraded compared to the baseline 3D-Var configuration. Over the Singapore radar domain, there was a general improvement in the precipitation forecasts, especially when the weighting toward the climatological background error covariance was larger, and with the additional application of time-shifted ensemble perturbations. Future work involves consolidating the ensemble prediction and deterministic system, by centering the ensemble prediction system on the hybrid analysis, to better represent the analysis and forecast uncertainties.
Abstract
This study compares aerosol direct radiative effects on numerical weather forecasts made by the NCEP Global Forecast Systems (GFS) with two different aerosol datasets, the OPAC and MERRA2 aerosol climatologies. The underestimation aerosol optical depth (AOD) by OPAC over northwest Africa, central to east Africa, Arabian Peninsula, Southeast Asia, and the Indo-Gangetic Plain, and overestimation in the storm track regions in both hemispheres are reduced by MERRA2. Surface downward short-wave (SW) and long-wave (LW) fluxes and the top of the atmosphere SW and outgoing LW fluxes from model forecasts are compared with CERES satellite observations. Forecasts made with OPAC aerosols have large radiative flux biases, especially in northwest Africa and the storm track regions. These biases are also reduced in the forecasts made with MERRA2 aerosols. The improvements from MERRA2 are most noticeable in the surface downward SW fluxes. GFS medium-range weather forecasts made with the MERRA2 aerosols demonstrated slightly improved forecast accuracy of sea level pressure and precipitation over the India and East Asian summer monsoon region. A stronger Africa easterly jet is produced, associated with a low pressure over the east Atlantic Ocean and west of north-west Africa. Impacts on large-scale skill scores such as 500hPa geopotential height anomaly correlation are generally positive in the Northern Hemisphere and the Pacific and North American regions in both the winter and summer seasons.
Abstract
This study compares aerosol direct radiative effects on numerical weather forecasts made by the NCEP Global Forecast Systems (GFS) with two different aerosol datasets, the OPAC and MERRA2 aerosol climatologies. The underestimation aerosol optical depth (AOD) by OPAC over northwest Africa, central to east Africa, Arabian Peninsula, Southeast Asia, and the Indo-Gangetic Plain, and overestimation in the storm track regions in both hemispheres are reduced by MERRA2. Surface downward short-wave (SW) and long-wave (LW) fluxes and the top of the atmosphere SW and outgoing LW fluxes from model forecasts are compared with CERES satellite observations. Forecasts made with OPAC aerosols have large radiative flux biases, especially in northwest Africa and the storm track regions. These biases are also reduced in the forecasts made with MERRA2 aerosols. The improvements from MERRA2 are most noticeable in the surface downward SW fluxes. GFS medium-range weather forecasts made with the MERRA2 aerosols demonstrated slightly improved forecast accuracy of sea level pressure and precipitation over the India and East Asian summer monsoon region. A stronger Africa easterly jet is produced, associated with a low pressure over the east Atlantic Ocean and west of north-west Africa. Impacts on large-scale skill scores such as 500hPa geopotential height anomaly correlation are generally positive in the Northern Hemisphere and the Pacific and North American regions in both the winter and summer seasons.
Abstract
Tropical cyclone tornadoes (TCTORs) are a hazard to life and property during landfalling tropical cyclones (TCs). The threat is often spread over a wide area within the TC envelope, and must be continually evaluated as the TC moves inland and dissipates. To anticipate the risk of TCTORs, forecasters may use high-resolution, rapidly updating model analyses and short-range forecasts such as the Rapid Refresh (RAP) and High-resolution Rapid Refresh (HRRR), and an ingredients-based approach similar to that used for forecasting continental midlatitude tornadoes. Though RAP and HRRR errors have been identified in typical midlatitude convective environments, this study evaluates the performance of the RAP and the HRRR within the TC envelope, with particular attention given to sounding-derived parameters previously identified as useful for TCTOR forecasting.
A sample of 1,730 observed upper-air soundings is sourced from 13 TCs that made landfall along the US coastline between 2017–2019. The observed soundings are paired with their corresponding model gridpoint soundings from the RAP analysis, RAP 12-hour forecast, and HRRR 12-hour forecast. Model errors are calculated for both the raw sounding variables of temperature, dewpoint, and wind speed, as well as for the selected sounding-derived parameters. Results show a moist bias that worsens with height across all model runs. There are also statistically significant underpredictions in stability-related parameters such as convective available potential energy (CAPE) and kinematic parameters such as vertical wind shear.
Abstract
Tropical cyclone tornadoes (TCTORs) are a hazard to life and property during landfalling tropical cyclones (TCs). The threat is often spread over a wide area within the TC envelope, and must be continually evaluated as the TC moves inland and dissipates. To anticipate the risk of TCTORs, forecasters may use high-resolution, rapidly updating model analyses and short-range forecasts such as the Rapid Refresh (RAP) and High-resolution Rapid Refresh (HRRR), and an ingredients-based approach similar to that used for forecasting continental midlatitude tornadoes. Though RAP and HRRR errors have been identified in typical midlatitude convective environments, this study evaluates the performance of the RAP and the HRRR within the TC envelope, with particular attention given to sounding-derived parameters previously identified as useful for TCTOR forecasting.
A sample of 1,730 observed upper-air soundings is sourced from 13 TCs that made landfall along the US coastline between 2017–2019. The observed soundings are paired with their corresponding model gridpoint soundings from the RAP analysis, RAP 12-hour forecast, and HRRR 12-hour forecast. Model errors are calculated for both the raw sounding variables of temperature, dewpoint, and wind speed, as well as for the selected sounding-derived parameters. Results show a moist bias that worsens with height across all model runs. There are also statistically significant underpredictions in stability-related parameters such as convective available potential energy (CAPE) and kinematic parameters such as vertical wind shear.
Abstract
Herein, 14 severe quasi-linear convective systems (QLCS) covering a wide range of geographical locations and environmental conditions are simulated for both 1- and 3-km horizontal grid resolutions, to further clarify their comparative capabilities in representing convective system features associated with severe weather production. Emphasis is placed on validating the simulated reflectivity structures, cold pool strength, mesoscale vortex characteristics, and surface wind strength. As to the overall reflectivity characteristics, the basic leading-line trailing stratiform structure was often better defined at 1 versus 3 km, but both resolutions were capable of producing bow echo and line echo wave pattern type features. Cold pool characteristics for both the 1- and 3-km simulations were also well replicated for the differing environments, with the 1-km cold pools slightly colder and often a bit larger. Both resolutions captured the larger mesoscale vortices, such as line-end or bookend vortices, but smaller, leading-line mesoscale updraft vortices, that often promote QLCS tornadogenesis, were largely absent in the 3-km simulations. Finally, while maximum surface winds were only marginally well predicted for both resolutions, the simulations were able to reasonably differentiate the relative contributions of the cold pool versus mesoscale vortices. The present results suggest that while many QLCS characteristics can be reasonably represented at a grid scale of 3 km, some of the more detailed structures, such as overall reflectivity characteristics and the smaller leading-line mesoscale vortices would likely benefit from the finer 1-km grid spacing.
Significance Statement
High-resolution model forecasts using 3-km grid spacing have proven to offer significant forecast guidance enhancements for severe convective weather. However, it is unclear whether additional enhancements can be obtained by decreasing grid spacings further to 1 km. Herein, we compare forecasts of severe quasi-linear convective systems (QLCS) simulated using 1- versus 3-km grids to document the potential value added of such increases in grid resolutions. It is shown that some significant improvements can be obtained in the representation of many QLCS features, especially as regards reflectivity structure and in the development of small, leading-line mesoscale vortices that can contribute to both severe surface wind and tornado production.
Abstract
Herein, 14 severe quasi-linear convective systems (QLCS) covering a wide range of geographical locations and environmental conditions are simulated for both 1- and 3-km horizontal grid resolutions, to further clarify their comparative capabilities in representing convective system features associated with severe weather production. Emphasis is placed on validating the simulated reflectivity structures, cold pool strength, mesoscale vortex characteristics, and surface wind strength. As to the overall reflectivity characteristics, the basic leading-line trailing stratiform structure was often better defined at 1 versus 3 km, but both resolutions were capable of producing bow echo and line echo wave pattern type features. Cold pool characteristics for both the 1- and 3-km simulations were also well replicated for the differing environments, with the 1-km cold pools slightly colder and often a bit larger. Both resolutions captured the larger mesoscale vortices, such as line-end or bookend vortices, but smaller, leading-line mesoscale updraft vortices, that often promote QLCS tornadogenesis, were largely absent in the 3-km simulations. Finally, while maximum surface winds were only marginally well predicted for both resolutions, the simulations were able to reasonably differentiate the relative contributions of the cold pool versus mesoscale vortices. The present results suggest that while many QLCS characteristics can be reasonably represented at a grid scale of 3 km, some of the more detailed structures, such as overall reflectivity characteristics and the smaller leading-line mesoscale vortices would likely benefit from the finer 1-km grid spacing.
Significance Statement
High-resolution model forecasts using 3-km grid spacing have proven to offer significant forecast guidance enhancements for severe convective weather. However, it is unclear whether additional enhancements can be obtained by decreasing grid spacings further to 1 km. Herein, we compare forecasts of severe quasi-linear convective systems (QLCS) simulated using 1- versus 3-km grids to document the potential value added of such increases in grid resolutions. It is shown that some significant improvements can be obtained in the representation of many QLCS features, especially as regards reflectivity structure and in the development of small, leading-line mesoscale vortices that can contribute to both severe surface wind and tornado production.
Abstract
This study employs a long time series (1997–2017) of reforecasts based on a version of the ECMWF Integrated Forecast System to evaluate the dependence of medium-range (i.e., 3–15 days) precipitation forecast skill over California on the state of the large-scale atmospheric flow. As a basis for this evaluation, four recurrent large-scale flow regimes over the North Pacific and western North America associated with precipitation in a domain encompassing northern and central California were objectively identified in ECMWF ERA5 reanalysis data for November–March 1981–2017. Two of the regimes are characterized by zonal upper-level flow across the North Pacific, and the other two are characterized by wavy, blocked flow. Forecast verification statistics conditioned on regime occurrence indicate considerably lower medium-range precipitation skill over California in blocking regimes than in zonal regimes. Moreover, forecasts of blocking regimes tend to exhibit larger errors and uncertainty in the synoptic-scale flow over the eastern North Pacific and western North America compared with forecasts of zonal regimes. Composite analyses for blocking forecasts reveal a tendency for errors to develop in conjunction with the amplification of a ridge over the western and central North Pacific. The errors in the ridge tend to be communicated through the large-scale Rossby wave pattern, resulting in misforecasting of downstream trough amplification and, thereby, moisture flux and precipitation over California. The composites additionally indicate that error growth in the blocking ridge can be linked to misrepresentation of baroclinic development as well as upper-level divergent outflow associated with latent heat release.
Abstract
This study employs a long time series (1997–2017) of reforecasts based on a version of the ECMWF Integrated Forecast System to evaluate the dependence of medium-range (i.e., 3–15 days) precipitation forecast skill over California on the state of the large-scale atmospheric flow. As a basis for this evaluation, four recurrent large-scale flow regimes over the North Pacific and western North America associated with precipitation in a domain encompassing northern and central California were objectively identified in ECMWF ERA5 reanalysis data for November–March 1981–2017. Two of the regimes are characterized by zonal upper-level flow across the North Pacific, and the other two are characterized by wavy, blocked flow. Forecast verification statistics conditioned on regime occurrence indicate considerably lower medium-range precipitation skill over California in blocking regimes than in zonal regimes. Moreover, forecasts of blocking regimes tend to exhibit larger errors and uncertainty in the synoptic-scale flow over the eastern North Pacific and western North America compared with forecasts of zonal regimes. Composite analyses for blocking forecasts reveal a tendency for errors to develop in conjunction with the amplification of a ridge over the western and central North Pacific. The errors in the ridge tend to be communicated through the large-scale Rossby wave pattern, resulting in misforecasting of downstream trough amplification and, thereby, moisture flux and precipitation over California. The composites additionally indicate that error growth in the blocking ridge can be linked to misrepresentation of baroclinic development as well as upper-level divergent outflow associated with latent heat release.
Abstract
A time–space shift method is developed for relocating model-predicted tornado vortices to radar-observed locations to improve the model initial conditions and subsequent predictions of tornadoes. The method consists of the following three steps. (i) Use the vortex center location estimated from radar observations to sample the best ensemble member from tornado-resolving ensemble predictions. Here, the best member is defined in terms of the predicted vortex center track that has a closest point, say at the time of t = t*, to the estimated vortex center at the initial time t0 (when the tornado vortex signature is first detected in radar observations). (ii) Create a time-shifted field from the best ensemble member in which the field within a circular area of about 10-km radius around the vortex center is taken from t = t*, while the field outside this circular area is transformed smoothly via temporal interpolation to the best ensemble member at t0. (iii) Create a time–space-shifted field in which the above time-shifted circular area is further shifted horizontally to co-center with the estimated vortex center at t0, while the field outside this circular area is transformed smoothly via spatial interpolation to the non-shifted field at t0 from the best ensemble member. The method is applied to the 20 May 2013 Oklahoma Newcastle–Moore tornado case, and is shown to be very effective in improving the tornado track and intensity predictions.
Significance Statement
The time–space shift method developed in this paper can smoothly relocate tornado vortices in model-predicted fields to match radar-observed locations. The method is found to be very effective in improving not only model initial condition but also the subsequent tornado track and intensity predictions. The method is also not sensitive to small errors in radar-estimated vortex center location at the initial time. The method should be useful for future real-time or even operational applications although further tests and improvements are needed (and are planned).
Abstract
A time–space shift method is developed for relocating model-predicted tornado vortices to radar-observed locations to improve the model initial conditions and subsequent predictions of tornadoes. The method consists of the following three steps. (i) Use the vortex center location estimated from radar observations to sample the best ensemble member from tornado-resolving ensemble predictions. Here, the best member is defined in terms of the predicted vortex center track that has a closest point, say at the time of t = t*, to the estimated vortex center at the initial time t0 (when the tornado vortex signature is first detected in radar observations). (ii) Create a time-shifted field from the best ensemble member in which the field within a circular area of about 10-km radius around the vortex center is taken from t = t*, while the field outside this circular area is transformed smoothly via temporal interpolation to the best ensemble member at t0. (iii) Create a time–space-shifted field in which the above time-shifted circular area is further shifted horizontally to co-center with the estimated vortex center at t0, while the field outside this circular area is transformed smoothly via spatial interpolation to the non-shifted field at t0 from the best ensemble member. The method is applied to the 20 May 2013 Oklahoma Newcastle–Moore tornado case, and is shown to be very effective in improving the tornado track and intensity predictions.
Significance Statement
The time–space shift method developed in this paper can smoothly relocate tornado vortices in model-predicted fields to match radar-observed locations. The method is found to be very effective in improving not only model initial condition but also the subsequent tornado track and intensity predictions. The method is also not sensitive to small errors in radar-estimated vortex center location at the initial time. The method should be useful for future real-time or even operational applications although further tests and improvements are needed (and are planned).
Abstract
Hail forecasts produced by the CAM-HAILCAST pseudo-Lagrangian hail size forecasting model were evaluated during the 2019, 2020, and 2021 NOAA Hazardous Weather Testbed (HWT) Spring Forecasting Experiments (SFEs). As part of this evaluation, HWT SFE participants were polled about their definition of a “good” hail forecast. Participants were presented with two different verification methods conducted over three different spatiotemporal scales, and were then asked to subjectively evaluate the hail forecast as well as the different verification methods themselves. Results recommended use of multiple verification methods tailored to the type of forecast expected by the end-user interpreting and applying the forecast. The hail forecasts evaluated during this period included an implementation of CAM-HAILCAST in the Limited Area Model of the Unified Forecast System with the Finite Volume 3 (FV3) dynamical core. Evaluation of FV3-HAILCAST over both 1- and 24-h periods found continued improvement from 2019 to 2021. The improvement was largely a result of wide intervariability among FV3 ensemble members with different microphysics parameterizations in 2019 lessening significantly during 2020 and 2021. Overprediction throughout the diurnal cycle also lessened by 2021. A combination of both upscaling neighborhood verification and an object-based technique that only retained matched convective objects was necessary to understand the improvement, agreeing with the HWT SFE participants’ recommendations for multiple verification methods.
Significance Statement
“Good” forecasts of hail can be determined in multiple ways and must depend on both the performance of the guidance and the perspective of the end-user. This work looks at different verification strategies to capture the performance of the CAM-HAILCAST hail forecasting model across three years of the Spring Forecasting Experiment (SFE) in different parent models. Verification strategies were informed by SFE participant input via a survey. Skill variability among models decreased in SFE 2021 relative to prior SFEs. The FV3 model in 2021, compared to 2019, provided improved forecasts of both convective distribution and 38-mm (1.5 in.) hail size, as well as less overforecasting of convection from 1900 to 2300 UTC.
Abstract
Hail forecasts produced by the CAM-HAILCAST pseudo-Lagrangian hail size forecasting model were evaluated during the 2019, 2020, and 2021 NOAA Hazardous Weather Testbed (HWT) Spring Forecasting Experiments (SFEs). As part of this evaluation, HWT SFE participants were polled about their definition of a “good” hail forecast. Participants were presented with two different verification methods conducted over three different spatiotemporal scales, and were then asked to subjectively evaluate the hail forecast as well as the different verification methods themselves. Results recommended use of multiple verification methods tailored to the type of forecast expected by the end-user interpreting and applying the forecast. The hail forecasts evaluated during this period included an implementation of CAM-HAILCAST in the Limited Area Model of the Unified Forecast System with the Finite Volume 3 (FV3) dynamical core. Evaluation of FV3-HAILCAST over both 1- and 24-h periods found continued improvement from 2019 to 2021. The improvement was largely a result of wide intervariability among FV3 ensemble members with different microphysics parameterizations in 2019 lessening significantly during 2020 and 2021. Overprediction throughout the diurnal cycle also lessened by 2021. A combination of both upscaling neighborhood verification and an object-based technique that only retained matched convective objects was necessary to understand the improvement, agreeing with the HWT SFE participants’ recommendations for multiple verification methods.
Significance Statement
“Good” forecasts of hail can be determined in multiple ways and must depend on both the performance of the guidance and the perspective of the end-user. This work looks at different verification strategies to capture the performance of the CAM-HAILCAST hail forecasting model across three years of the Spring Forecasting Experiment (SFE) in different parent models. Verification strategies were informed by SFE participant input via a survey. Skill variability among models decreased in SFE 2021 relative to prior SFEs. The FV3 model in 2021, compared to 2019, provided improved forecasts of both convective distribution and 38-mm (1.5 in.) hail size, as well as less overforecasting of convection from 1900 to 2300 UTC.
