Browse
Abstract
Severe weather probabilities are derived from the Warn-on-Forecast System (WoFS) run by NOAA’s National Severe Storms Laboratory (NSSL) during spring 2018 using the random forest (RF) machine learning algorithm. Recent work has shown this method generates skillful and reliable forecasts when applied to convection-allowing model ensembles for the “Day 1” time range (i.e., 12–36-h lead times), but it has been tested in only one other study for lead times relevant to WoFS (e.g., 0–6 h). Thus, in this paper, various sets of WoFS predictors, which include both environment and storm-based fields, are input into a RF algorithm and trained using the occurrence of severe weather reports within 39 km of a point to produce severe weather probabilities at 0–3-h lead times. We analyze the skill and reliability of these forecasts, sensitivity to different sets of predictors, and avenues for further improvements. The RF algorithm produced very skillful and reliable severe weather probabilities and significantly outperformed baseline probabilities calculated by finding the best performing updraft helicity (UH) threshold and smoothing parameter. Experiments where different sets of predictors were used to derive RF probabilities revealed 1) storm attribute fields contributed significantly more skill than environmental fields, 2) 2–5 km AGL UH and maximum updraft speed were the best performing storm attribute fields, 3) the most skillful ensemble summary metric was a smoothed mean, and 4) the most skillful forecasts were obtained when smoothed UH from individual ensemble members were used as predictors.
Abstract
Severe weather probabilities are derived from the Warn-on-Forecast System (WoFS) run by NOAA’s National Severe Storms Laboratory (NSSL) during spring 2018 using the random forest (RF) machine learning algorithm. Recent work has shown this method generates skillful and reliable forecasts when applied to convection-allowing model ensembles for the “Day 1” time range (i.e., 12–36-h lead times), but it has been tested in only one other study for lead times relevant to WoFS (e.g., 0–6 h). Thus, in this paper, various sets of WoFS predictors, which include both environment and storm-based fields, are input into a RF algorithm and trained using the occurrence of severe weather reports within 39 km of a point to produce severe weather probabilities at 0–3-h lead times. We analyze the skill and reliability of these forecasts, sensitivity to different sets of predictors, and avenues for further improvements. The RF algorithm produced very skillful and reliable severe weather probabilities and significantly outperformed baseline probabilities calculated by finding the best performing updraft helicity (UH) threshold and smoothing parameter. Experiments where different sets of predictors were used to derive RF probabilities revealed 1) storm attribute fields contributed significantly more skill than environmental fields, 2) 2–5 km AGL UH and maximum updraft speed were the best performing storm attribute fields, 3) the most skillful ensemble summary metric was a smoothed mean, and 4) the most skillful forecasts were obtained when smoothed UH from individual ensemble members were used as predictors.
Abstract
The development of sea surface temperature (SST) anomalies over the northeast Pacific and their impacts on lower-tropospheric air temperatures over the Pacific Northwest are examined. Northeast Pacific SST anomalies are influenced by the synoptic-scale flow, with high pressure and weak surface winds associated with developing warm SST anomalies, while large pressure gradients and strong surface winds result in SST declines. SST over the northeast Pacific correlates significantly with surface air temperatures over the Pacific Northwest, with correlations increasing when high-frequency variability is filtered out. The correlations between unfiltered time series of SST and surface air temperature are largest for a zero-day lag and are strongest near the coast, contrasting with weaker correlations over the Columbia basin east of the Cascade Mountains. SST correlations with minimum surface air temperature are largest during the warm season, and maximum temperature correlations are highest in March; both have low correlations during autumn. Model simulations of periods with warm and cold northeast Pacific SST anomalies possess lower-tropospheric air temperature warming or cooling over the coastal zone, with SST influence weakening east of the Cascade crest. Eastern Pacific SST anomalies influence sea level pressure and lower-tropospheric heights, with warm SST anomalies resulting in simulated lowered pressure near the surface and increased heights aloft. The relationship between northeast Pacific SST and surface air temperature over land evince complex feedbacks: SST temperature anomalies can be advected inland from the Pacific, the SST anomalies can influence the synoptic-scale flow that affects the SST anomalies, and the synoptic-scale anomalies that produce the SST anomalies can directly influence temperatures over land.
Significance Statement
Understanding the connection between northeast Pacific sea surface temperatures and low-level air temperatures over land is valuable for both subseasonal prediction and for examining the fidelity of model physics.
Abstract
The development of sea surface temperature (SST) anomalies over the northeast Pacific and their impacts on lower-tropospheric air temperatures over the Pacific Northwest are examined. Northeast Pacific SST anomalies are influenced by the synoptic-scale flow, with high pressure and weak surface winds associated with developing warm SST anomalies, while large pressure gradients and strong surface winds result in SST declines. SST over the northeast Pacific correlates significantly with surface air temperatures over the Pacific Northwest, with correlations increasing when high-frequency variability is filtered out. The correlations between unfiltered time series of SST and surface air temperature are largest for a zero-day lag and are strongest near the coast, contrasting with weaker correlations over the Columbia basin east of the Cascade Mountains. SST correlations with minimum surface air temperature are largest during the warm season, and maximum temperature correlations are highest in March; both have low correlations during autumn. Model simulations of periods with warm and cold northeast Pacific SST anomalies possess lower-tropospheric air temperature warming or cooling over the coastal zone, with SST influence weakening east of the Cascade crest. Eastern Pacific SST anomalies influence sea level pressure and lower-tropospheric heights, with warm SST anomalies resulting in simulated lowered pressure near the surface and increased heights aloft. The relationship between northeast Pacific SST and surface air temperature over land evince complex feedbacks: SST temperature anomalies can be advected inland from the Pacific, the SST anomalies can influence the synoptic-scale flow that affects the SST anomalies, and the synoptic-scale anomalies that produce the SST anomalies can directly influence temperatures over land.
Significance Statement
Understanding the connection between northeast Pacific sea surface temperatures and low-level air temperatures over land is valuable for both subseasonal prediction and for examining the fidelity of model physics.
Abstract
In this work, we characterized the occurrences and conditions before the initiations of mesoscale convective systems (MCSs) in the central United States, using 15 years of observations and convection-permitting climate model simulations. The variabilities of MCSs in summer were obtained using high-resolution (4 km) observation data [Stage-IV (stIV)] and ECMWF Re-Analysis v5 (ERA5)-forced Weather Research and Forecasting (WRF) Model simulations (E5RUN). MCSs were identified using the object tracking algorithm MODE-time domain (MTD). MTD-determined MCSs were divided into daytime short-lived MCSs (SLM12), daytime long-lived MCSs (LLM12), nighttime short-lived MCSs (SLM00), and nighttime long-lived MCSs (LLM00). E5RUN showed skill to simulate MCSs by obtaining similar statistics in occurrences, areal coverages, and propagation speeds compared to those of stIV. We calculated the 15 parameters using sounding data from E5RUN before an MCS was initiated (−1, −3, −6, and −9 h) at each location of an MCS. The parameters were tested to figure out the significance of predicting the longevities of MCSs. The key findings are 1) LLM12 showed favorable thermodynamic variables compared to that of SLM12 and 2) LLM00 showed significant conditions of vertically rotating winds and sheared environments that affect the longevity of MCSs. Moreover, storm-relative helicity of 0–3 km, precipitable water, and vertical wind shear of 0–6 km are the most significant parameters to determine the longevities of MCSs (both daytime and nighttime MCSs).
Significance Statement
The purpose of this study is to understand the features of mesoscale convective systems (MCSs) in observational data and convection-permitting climate model simulations. We tested long-term simulations using new forcing data (ERA5) to see the benefits and limitations. We designed a novel approach to obtain the distributions of meteorological parameters (instead of obtaining one value for one event of MCS) before initiations of MCSs to understand preconvective conditions (times from −9 to −1 h from initiation). We also divided MCSs into daytime/nighttime and short-/long-lived MCSs to help predict MCSs longevity considering the initiation times. Our results provide hints for the forecasters to predict MCS longevity based on preconvective conditions from parameters discussed in this work.
Abstract
In this work, we characterized the occurrences and conditions before the initiations of mesoscale convective systems (MCSs) in the central United States, using 15 years of observations and convection-permitting climate model simulations. The variabilities of MCSs in summer were obtained using high-resolution (4 km) observation data [Stage-IV (stIV)] and ECMWF Re-Analysis v5 (ERA5)-forced Weather Research and Forecasting (WRF) Model simulations (E5RUN). MCSs were identified using the object tracking algorithm MODE-time domain (MTD). MTD-determined MCSs were divided into daytime short-lived MCSs (SLM12), daytime long-lived MCSs (LLM12), nighttime short-lived MCSs (SLM00), and nighttime long-lived MCSs (LLM00). E5RUN showed skill to simulate MCSs by obtaining similar statistics in occurrences, areal coverages, and propagation speeds compared to those of stIV. We calculated the 15 parameters using sounding data from E5RUN before an MCS was initiated (−1, −3, −6, and −9 h) at each location of an MCS. The parameters were tested to figure out the significance of predicting the longevities of MCSs. The key findings are 1) LLM12 showed favorable thermodynamic variables compared to that of SLM12 and 2) LLM00 showed significant conditions of vertically rotating winds and sheared environments that affect the longevity of MCSs. Moreover, storm-relative helicity of 0–3 km, precipitable water, and vertical wind shear of 0–6 km are the most significant parameters to determine the longevities of MCSs (both daytime and nighttime MCSs).
Significance Statement
The purpose of this study is to understand the features of mesoscale convective systems (MCSs) in observational data and convection-permitting climate model simulations. We tested long-term simulations using new forcing data (ERA5) to see the benefits and limitations. We designed a novel approach to obtain the distributions of meteorological parameters (instead of obtaining one value for one event of MCS) before initiations of MCSs to understand preconvective conditions (times from −9 to −1 h from initiation). We also divided MCSs into daytime/nighttime and short-/long-lived MCSs to help predict MCSs longevity considering the initiation times. Our results provide hints for the forecasters to predict MCS longevity based on preconvective conditions from parameters discussed in this work.
Abstract
A high-resolution, short-term climate prediction system for summer (June–July–August) climate over Southwest China has been developed using the Weather Research and Forecasting (WRF) Model nested with a global climate prediction system (PCCSM4). The system includes 12 ensemble members generated by PCCSM4 with different initial conditions, and the finest horizontal resolution of WRF is 8 km. This study evaluates the ability of the WRF Model to predict summer climate over Southwest China, focusing on the system design, model tuning, and evaluation of baseline model performance. Sensitivity simulations are first conducted to provide the optimal model configuration, and the model performance is evaluated against available observational data using reforecast simulations for 1981–2020. When compared to PCCSM4, the WRF Model shows major improvements in predicting the spatial distribution of major variables such as 2-m temperature, 10-m wind speed, and precipitation. WRF also shows better skill in predicting interannual temperature variability and extreme temperature events, with higher anomaly correlation coefficients. However, large model biases remain in seasonal precipitation anomaly predictions. Overall, this study highlights the potential advantages of using the high-resolution WRF Model to predict summer climate conditions over Southwest China.
Abstract
A high-resolution, short-term climate prediction system for summer (June–July–August) climate over Southwest China has been developed using the Weather Research and Forecasting (WRF) Model nested with a global climate prediction system (PCCSM4). The system includes 12 ensemble members generated by PCCSM4 with different initial conditions, and the finest horizontal resolution of WRF is 8 km. This study evaluates the ability of the WRF Model to predict summer climate over Southwest China, focusing on the system design, model tuning, and evaluation of baseline model performance. Sensitivity simulations are first conducted to provide the optimal model configuration, and the model performance is evaluated against available observational data using reforecast simulations for 1981–2020. When compared to PCCSM4, the WRF Model shows major improvements in predicting the spatial distribution of major variables such as 2-m temperature, 10-m wind speed, and precipitation. WRF also shows better skill in predicting interannual temperature variability and extreme temperature events, with higher anomaly correlation coefficients. However, large model biases remain in seasonal precipitation anomaly predictions. Overall, this study highlights the potential advantages of using the high-resolution WRF Model to predict summer climate conditions over Southwest China.
Abstract
Tropical precipitation and circulation are often coupled and span a vast spectrum of scales from a few to several thousands of kilometers and from hours to weeks. Current operational numerical weather prediction (NWP) models struggle with representing the full range of scales of tropical phenomena. Synoptic to planetary scales are of particular importance because improved skill in the representation of tropical larger-scale features such as convectively coupled equatorial waves (CCEWs) has the potential to reduce forecast error propagation from the tropics to the midlatitudes. Here we introduce diagnostics from a recently developed tropical variability diagnostics toolbox, where we focus on two recent versions of NOAA’s Unified Forecast System (UFS): operational GFSv15 forecasts and experimental GFSv16 forecasts from April to October 2020. The diagnostics include space–time coherence spectra to identify preferred scales of coupling between circulation and precipitation, pattern correlations of Hovmöller diagrams to assess model skill in zonal propagation of precipitating features, CCEW skill assessment, plus a diagnostic aimed at evaluating moisture–convection coupling in the tropics. Results show that the GFSv16 forecasts are slightly more realistic than GFSv15 in their coherence between precipitation and model dynamics at synoptic to planetary scales, with modest improvements in moisture convection coupling. However, this slightly improved performance does not necessarily translate to improvements in traditional precipitation skill scores. The results highlight the utility of these diagnostics in the pursuit of better understanding of NWP model performance in the tropics, while also demonstrating the challenges in translating model advancements into improved skill.
Abstract
Tropical precipitation and circulation are often coupled and span a vast spectrum of scales from a few to several thousands of kilometers and from hours to weeks. Current operational numerical weather prediction (NWP) models struggle with representing the full range of scales of tropical phenomena. Synoptic to planetary scales are of particular importance because improved skill in the representation of tropical larger-scale features such as convectively coupled equatorial waves (CCEWs) has the potential to reduce forecast error propagation from the tropics to the midlatitudes. Here we introduce diagnostics from a recently developed tropical variability diagnostics toolbox, where we focus on two recent versions of NOAA’s Unified Forecast System (UFS): operational GFSv15 forecasts and experimental GFSv16 forecasts from April to October 2020. The diagnostics include space–time coherence spectra to identify preferred scales of coupling between circulation and precipitation, pattern correlations of Hovmöller diagrams to assess model skill in zonal propagation of precipitating features, CCEW skill assessment, plus a diagnostic aimed at evaluating moisture–convection coupling in the tropics. Results show that the GFSv16 forecasts are slightly more realistic than GFSv15 in their coherence between precipitation and model dynamics at synoptic to planetary scales, with modest improvements in moisture convection coupling. However, this slightly improved performance does not necessarily translate to improvements in traditional precipitation skill scores. The results highlight the utility of these diagnostics in the pursuit of better understanding of NWP model performance in the tropics, while also demonstrating the challenges in translating model advancements into improved skill.
Abstract
This study reviews the recent addition of dropwindsonde wind data near the tropical cyclone (TC) center as well as the first-time addition of high-density, flight-level reconnaissance observations (HDOBs) into the National Centers for Environmental Prediction (NCEP) Global Forecast System (GFS). The main finding is that the additional data have profound positive impacts on subsequent TC track forecasts. For TCs in the North Atlantic (NATL) basin, statistically significant improvements in track extend through 4–5 days during reconnaissance periods. Further assessment suggests that greater improvements might also be expected at days 6–7. This study also explores the importance of comprehensively assessing data impact. For example, model or data assimilation changes can affect the so-called “early” and “late” versions of the forecast very differently. It is also important to explore different ways to describe the error statistics. In several instances the impacts of the additional data strongly differ depending on whether one examines the mean or median errors. The results demonstrate the tremendous potential for further improving TC forecasts. The data added here were already operationally transmitted and assimilated by other systems at NCEP, and many further improvements likely await with improved use of these and other reconnaissance observations. This demonstrates the need of not only investing in data assimilation improvements, but also enhancements to observational systems in order to reach next-generation hurricane forecasting goals.
Significance Statement
This study demonstrates that data gathered from reconnaissance missions into tropical cyclones substantially improves tropical cyclone track forecasts.
Abstract
This study reviews the recent addition of dropwindsonde wind data near the tropical cyclone (TC) center as well as the first-time addition of high-density, flight-level reconnaissance observations (HDOBs) into the National Centers for Environmental Prediction (NCEP) Global Forecast System (GFS). The main finding is that the additional data have profound positive impacts on subsequent TC track forecasts. For TCs in the North Atlantic (NATL) basin, statistically significant improvements in track extend through 4–5 days during reconnaissance periods. Further assessment suggests that greater improvements might also be expected at days 6–7. This study also explores the importance of comprehensively assessing data impact. For example, model or data assimilation changes can affect the so-called “early” and “late” versions of the forecast very differently. It is also important to explore different ways to describe the error statistics. In several instances the impacts of the additional data strongly differ depending on whether one examines the mean or median errors. The results demonstrate the tremendous potential for further improving TC forecasts. The data added here were already operationally transmitted and assimilated by other systems at NCEP, and many further improvements likely await with improved use of these and other reconnaissance observations. This demonstrates the need of not only investing in data assimilation improvements, but also enhancements to observational systems in order to reach next-generation hurricane forecasting goals.
Significance Statement
This study demonstrates that data gathered from reconnaissance missions into tropical cyclones substantially improves tropical cyclone track forecasts.
Abstract
The approach of applying different amounts of horizontal localization to different ranges of background-error covariance horizontal scales as proposed by Buehner and Shlyaeva was recently implemented in the four-dimensional ensemble–variational (4DEnVar) data assimilation scheme of the global deterministic prediction system (GDPS) at Environment and Climate Change Canada operations. To maximize the benefits from this approach to reduce the sampling noise in the ensemble-derived background-error covariances, it was necessary to adopt a new weighting between the climatological and flow-dependent covariances that increases significantly the role of the latter. Thus, in December 2021 the GDPS became the first operational global deterministic medium-range weather forecasting system to rely completely on flow-dependent covariances in the troposphere and the lower stratosphere. The experiments that led to the adoption of these two related changes and their impacts on the forecasts up to 7 days for various regions of the globe during the boreal summer of 2019 and winter of 2020 are presented here. It is also illustrated that relying more on ensemble-derived covariances amplifies the positive impacts on the GDPS when the background ensemble generation strategy is improved.
Abstract
The approach of applying different amounts of horizontal localization to different ranges of background-error covariance horizontal scales as proposed by Buehner and Shlyaeva was recently implemented in the four-dimensional ensemble–variational (4DEnVar) data assimilation scheme of the global deterministic prediction system (GDPS) at Environment and Climate Change Canada operations. To maximize the benefits from this approach to reduce the sampling noise in the ensemble-derived background-error covariances, it was necessary to adopt a new weighting between the climatological and flow-dependent covariances that increases significantly the role of the latter. Thus, in December 2021 the GDPS became the first operational global deterministic medium-range weather forecasting system to rely completely on flow-dependent covariances in the troposphere and the lower stratosphere. The experiments that led to the adoption of these two related changes and their impacts on the forecasts up to 7 days for various regions of the globe during the boreal summer of 2019 and winter of 2020 are presented here. It is also illustrated that relying more on ensemble-derived covariances amplifies the positive impacts on the GDPS when the background ensemble generation strategy is improved.
Abstract
This study investigates the impact of assimilating ground-based radar reflectivity and wind data on tropical cyclone (TC) intensity prediction. The effect on a high-impact TC in the western North Pacific region that penetrated the Bashi Channel is examined. A multiscale correction based on the successive covariance localization (SCL) method is adopted to improve the analysis and forecast performance. In addition, GNSS-R wind speed is assimilated jointly in the rapid update assimilation framework to complement the TC boundary layer where radar data are limited. Model experiments are conducted and evaluated using the observing system simulation experiments (OSSEs) framework with a coupled atmosphere–ocean model nature run. Taking the experiment without data assimilation as the baseline, assimilating the radar data with the standard localization and SCL methods reduces the wind speed analysis error by 12% and 44%, respectively. The SCL method dominates the improvement in TC intensity prediction with a lead time longer than 2 days and the TC’s peak intensity forecast is improved by 18 hPa. The additional assimilation of the GNSS-R wind speed observation further reduces the wind speed error in the low-level analysis by 12% and 5% under the standard and SCL radar assimilation framework, respectively. GNSS-R wind assimilation leads to a further 6-hPa improvement in TC’s peak intensity. However, the sampling error introduced by the SCL method restrains the effect of GNSS-R assimilation. Sensitivity experiments with different GNSS-R data arrangements show that better GNSS-R wind coverage and additional wind direction information can further improve the TC analysis.
Significance Statement
Tropical cyclone (TC) intensity prediction over the western North Pacific (WNP) region remains a significant challenge due to limited observations. This study aims to improve the TC intensity prediction in WNP by assimilating the ground-based radar data using a multiscale correction framework and incorporating with the satellite ocean surface wind speed observation. We particularly focus on a high-impact TC like Typhoon Hato (2017), which penetrated the Bashi Channel and later made landfall in China, causing great damage. Our results showed that the assimilation strategy improved the TC intensity prediction for a lead time longer than 2 days. These results demonstrate the great potential of these observations and can provide guidance for future applications in operation centers.
Abstract
This study investigates the impact of assimilating ground-based radar reflectivity and wind data on tropical cyclone (TC) intensity prediction. The effect on a high-impact TC in the western North Pacific region that penetrated the Bashi Channel is examined. A multiscale correction based on the successive covariance localization (SCL) method is adopted to improve the analysis and forecast performance. In addition, GNSS-R wind speed is assimilated jointly in the rapid update assimilation framework to complement the TC boundary layer where radar data are limited. Model experiments are conducted and evaluated using the observing system simulation experiments (OSSEs) framework with a coupled atmosphere–ocean model nature run. Taking the experiment without data assimilation as the baseline, assimilating the radar data with the standard localization and SCL methods reduces the wind speed analysis error by 12% and 44%, respectively. The SCL method dominates the improvement in TC intensity prediction with a lead time longer than 2 days and the TC’s peak intensity forecast is improved by 18 hPa. The additional assimilation of the GNSS-R wind speed observation further reduces the wind speed error in the low-level analysis by 12% and 5% under the standard and SCL radar assimilation framework, respectively. GNSS-R wind assimilation leads to a further 6-hPa improvement in TC’s peak intensity. However, the sampling error introduced by the SCL method restrains the effect of GNSS-R assimilation. Sensitivity experiments with different GNSS-R data arrangements show that better GNSS-R wind coverage and additional wind direction information can further improve the TC analysis.
Significance Statement
Tropical cyclone (TC) intensity prediction over the western North Pacific (WNP) region remains a significant challenge due to limited observations. This study aims to improve the TC intensity prediction in WNP by assimilating the ground-based radar data using a multiscale correction framework and incorporating with the satellite ocean surface wind speed observation. We particularly focus on a high-impact TC like Typhoon Hato (2017), which penetrated the Bashi Channel and later made landfall in China, causing great damage. Our results showed that the assimilation strategy improved the TC intensity prediction for a lead time longer than 2 days. These results demonstrate the great potential of these observations and can provide guidance for future applications in operation centers.
Abstract
Microphysical perturbation experiments were conducted to investigate the sensitivity of convective heavy rain simulation to cloud microphysical parameterization and its feasibility for ensemble forecasts. An ensemble of 20 perturbation members differing in either the microphysics package or process treatments within a single scheme was applied to simulate 10 summer-afternoon heavy-rain convection cases. The simulations revealed substantial disagreements in the location and amplitude of peak rainfall among the microphysics-package and single-scheme members, with an overall spread of 57%–161%, 66%–161%, and 65%–149% of the observed average rainfall, maximum rainfall, and maximum intensity, respectively. The single-scheme members revealed that the simulation of heavy convective precipitation is quite sensitive to factors including ice-particle fall speed parameterization, aerosol type, ice particle shape, and size distribution representation. The microphysical ensemble can derive reasonable probability of occurrence for a location-specific heavy-rain forecast. Spatial-forecast performance indices up to 0.6 were attained by applying an optimal fuzzy radius of about 8 km for the warning-area coverage. The forecasts tend to be more successful for more organized convection. Spectral mapping methods were further applied to provide ensemble forecasts for the 10 heavy rainfall cases. For most cases, realistic spatial patterns were derived with spatial correlation up to 0.8. The quantitative performance in average rainfall, maximum rainfall, and maximum intensity from the ensembles reached correlations of 0.83, 0.84, and 0.51, respectively, with the observed values.
Significance Statement
Heavy rainfall from summer convections is stochastic in terms of intensity and location; therefore, an accurate deterministic forecast is often challenging. We designed perturbation experiments to explore weather forecasting models’ sensitivity to cloud microphysical parameterizations and the feasibility of application to ensemble forecast. Promising results were obtained from simulations of 10 real cases. The cloud microphysical ensemble approach may provide reasonable forecasts of heavy rainfall probability and convincing rainfall spatial distribution, particularly for more organized convection.
Abstract
Microphysical perturbation experiments were conducted to investigate the sensitivity of convective heavy rain simulation to cloud microphysical parameterization and its feasibility for ensemble forecasts. An ensemble of 20 perturbation members differing in either the microphysics package or process treatments within a single scheme was applied to simulate 10 summer-afternoon heavy-rain convection cases. The simulations revealed substantial disagreements in the location and amplitude of peak rainfall among the microphysics-package and single-scheme members, with an overall spread of 57%–161%, 66%–161%, and 65%–149% of the observed average rainfall, maximum rainfall, and maximum intensity, respectively. The single-scheme members revealed that the simulation of heavy convective precipitation is quite sensitive to factors including ice-particle fall speed parameterization, aerosol type, ice particle shape, and size distribution representation. The microphysical ensemble can derive reasonable probability of occurrence for a location-specific heavy-rain forecast. Spatial-forecast performance indices up to 0.6 were attained by applying an optimal fuzzy radius of about 8 km for the warning-area coverage. The forecasts tend to be more successful for more organized convection. Spectral mapping methods were further applied to provide ensemble forecasts for the 10 heavy rainfall cases. For most cases, realistic spatial patterns were derived with spatial correlation up to 0.8. The quantitative performance in average rainfall, maximum rainfall, and maximum intensity from the ensembles reached correlations of 0.83, 0.84, and 0.51, respectively, with the observed values.
Significance Statement
Heavy rainfall from summer convections is stochastic in terms of intensity and location; therefore, an accurate deterministic forecast is often challenging. We designed perturbation experiments to explore weather forecasting models’ sensitivity to cloud microphysical parameterizations and the feasibility of application to ensemble forecast. Promising results were obtained from simulations of 10 real cases. The cloud microphysical ensemble approach may provide reasonable forecasts of heavy rainfall probability and convincing rainfall spatial distribution, particularly for more organized convection.