Browse
Abstract
While convective storm mode is explicitly depicted in convection-allowing model (CAM) output, subjectively diagnosing mode in large volumes of CAM forecasts can be burdensome. In this work, four machine learning (ML) models were trained to probabilistically classify CAM storms into one of three modes: supercells, quasi-linear convective systems, and disorganized convection. The four ML models included a dense neural network (DNN), logistic regression (LR), a convolutional neural network (CNN), and semisupervised CNN–Gaussian mixture model (GMM). The DNN, CNN, and LR were trained with a set of hand-labeled CAM storms, while the semisupervised GMM used updraft helicity and storm size to generate clusters, which were then hand labeled. When evaluated using storms withheld from training, the four classifiers had similar ability to discriminate between modes, but the GMM had worse calibration. The DNN and LR had similar objective performance to the CNN, suggesting that CNN-based methods may not be needed for mode classification tasks. The mode classifications from all four classifiers successfully approximated the known climatology of modes in the United States, including a maximum in supercell occurrence in the U.S. Central Plains. Further, the modes also occurred in environments recognized to support the three different storm morphologies. Finally, storm mode provided useful information about hazard type, e.g., storm reports were most likely with supercells, further supporting the efficacy of the classifiers. Future applications, including the use of objective CAM mode classifications as a novel predictor in ML systems, could potentially lead to improved forecasts of convective hazards.
Significance Statement
Whether a thunderstorm produces hazards such as tornadoes, hail, or intense wind gusts is in part determined by whether the storm takes the form of a single cell or a line. Numerical forecasting models can now provide forecasts that depict this structure. We tested several automated algorithms to extract this information from forecast output using machine learning. All of the automated methods were able to distinguish between a set of three convective types, with the simple techniques providing similarly skilled classifications compared to the complex approaches. The automated classifications also successfully discriminated between thunderstorm hazards, potentially leading to new forecast tools and better forecasts of high-impact convective hazards.
Abstract
While convective storm mode is explicitly depicted in convection-allowing model (CAM) output, subjectively diagnosing mode in large volumes of CAM forecasts can be burdensome. In this work, four machine learning (ML) models were trained to probabilistically classify CAM storms into one of three modes: supercells, quasi-linear convective systems, and disorganized convection. The four ML models included a dense neural network (DNN), logistic regression (LR), a convolutional neural network (CNN), and semisupervised CNN–Gaussian mixture model (GMM). The DNN, CNN, and LR were trained with a set of hand-labeled CAM storms, while the semisupervised GMM used updraft helicity and storm size to generate clusters, which were then hand labeled. When evaluated using storms withheld from training, the four classifiers had similar ability to discriminate between modes, but the GMM had worse calibration. The DNN and LR had similar objective performance to the CNN, suggesting that CNN-based methods may not be needed for mode classification tasks. The mode classifications from all four classifiers successfully approximated the known climatology of modes in the United States, including a maximum in supercell occurrence in the U.S. Central Plains. Further, the modes also occurred in environments recognized to support the three different storm morphologies. Finally, storm mode provided useful information about hazard type, e.g., storm reports were most likely with supercells, further supporting the efficacy of the classifiers. Future applications, including the use of objective CAM mode classifications as a novel predictor in ML systems, could potentially lead to improved forecasts of convective hazards.
Significance Statement
Whether a thunderstorm produces hazards such as tornadoes, hail, or intense wind gusts is in part determined by whether the storm takes the form of a single cell or a line. Numerical forecasting models can now provide forecasts that depict this structure. We tested several automated algorithms to extract this information from forecast output using machine learning. All of the automated methods were able to distinguish between a set of three convective types, with the simple techniques providing similarly skilled classifications compared to the complex approaches. The automated classifications also successfully discriminated between thunderstorm hazards, potentially leading to new forecast tools and better forecasts of high-impact convective hazards.
Abstract
The virtual temperature used to model moisture-modified tropospheric dynamics is generalized to include a new thermospheric component. The resulting hybrid virtual potential temperature (HVPT) transitions seamlessly with height, from moist virtual potential temperature (MVPT) in the troposphere, to potential temperature in the stratosphere and mesosphere, to thermospheric virtual potential temperature thereafter. For numerical weather prediction (NWP) models looking to extend into the thermosphere, but still heavily invested in retaining MVPT-based dynamical cores for tropospheric prediction, upgrading to HVPT allows the core to capture critical new aspects of variable composition thermospheric dynamics, while leaving the original MVPT-based tropospheric equations and numerics essentially untouched. In this way, HVPT augmentation can both simplify and streamline extension into the thermosphere at little computational cost beyond the inevitable need for more vertical layers and somewhat smaller time steps. To demonstrate, we upgrade the MVPT-based dynamical core of the Navy global NWP model to HVPT, then test its performance in forecasting analytical globally balanced states containing hot or rapidly heated thermospheres and height-varying gas constants. These tests confirm that HVPT augmentation offers an efficient and effective means of extending MVPT-based NWP models into the thermosphere to accelerate development of future ground-to-space NWP models supporting space weather applications. The related issues of variable gravitational acceleration and shallow-atmosphere approximations are also briefly discussed.
Abstract
The virtual temperature used to model moisture-modified tropospheric dynamics is generalized to include a new thermospheric component. The resulting hybrid virtual potential temperature (HVPT) transitions seamlessly with height, from moist virtual potential temperature (MVPT) in the troposphere, to potential temperature in the stratosphere and mesosphere, to thermospheric virtual potential temperature thereafter. For numerical weather prediction (NWP) models looking to extend into the thermosphere, but still heavily invested in retaining MVPT-based dynamical cores for tropospheric prediction, upgrading to HVPT allows the core to capture critical new aspects of variable composition thermospheric dynamics, while leaving the original MVPT-based tropospheric equations and numerics essentially untouched. In this way, HVPT augmentation can both simplify and streamline extension into the thermosphere at little computational cost beyond the inevitable need for more vertical layers and somewhat smaller time steps. To demonstrate, we upgrade the MVPT-based dynamical core of the Navy global NWP model to HVPT, then test its performance in forecasting analytical globally balanced states containing hot or rapidly heated thermospheres and height-varying gas constants. These tests confirm that HVPT augmentation offers an efficient and effective means of extending MVPT-based NWP models into the thermosphere to accelerate development of future ground-to-space NWP models supporting space weather applications. The related issues of variable gravitational acceleration and shallow-atmosphere approximations are also briefly discussed.
Abstract
The prediction of weather conditions in the Arctic is important to human activities in the Arctic. Arctic cyclones (ACs), which are extratropical cyclones that originate within the Arctic or move into the Arctic from lower latitudes, can be associated with hazardous weather conditions that may adversely affect human activities. The purpose of this study is to increase understanding of processes that influence the forecast skill of the synoptic-scale flow over the Arctic and of ACs. The 11-member NOAA Global Ensemble Forecast System (GEFS) reforecast dataset, version 2, is utilized to identify periods of low and high forecast skill of the synoptic-scale flow over the Arctic, hereinafter referred to as low-skill and high-skill periods, respectively, during the summers of 2007–17, and to evaluate the forecast skill of ACs during these respective periods. The ERA-Interim dataset is used to examine characteristics of the Arctic environment and characteristics of ACs during low-skill and high-skill periods. The Arctic environment tends to be characterized by more vigorous baroclinic processes and latent heating during low-skill periods relative to high-skill periods. ACs occur more frequently over much of the Arctic; tend to be stronger; and tend to be located in regions of larger lower-tropospheric baroclinicity, lower-to-midtropospheric Eady growth rate (EGR), and latent heating during low-skill periods relative to high-skill periods. ACs during low-skill periods that are characterized by low forecast skill of intensity tend to be relatively strong and tend to be located in regions of relatively large lower-tropospheric baroclinicity, lower-to-midtropospheric EGR, and latent heating.
Abstract
The prediction of weather conditions in the Arctic is important to human activities in the Arctic. Arctic cyclones (ACs), which are extratropical cyclones that originate within the Arctic or move into the Arctic from lower latitudes, can be associated with hazardous weather conditions that may adversely affect human activities. The purpose of this study is to increase understanding of processes that influence the forecast skill of the synoptic-scale flow over the Arctic and of ACs. The 11-member NOAA Global Ensemble Forecast System (GEFS) reforecast dataset, version 2, is utilized to identify periods of low and high forecast skill of the synoptic-scale flow over the Arctic, hereinafter referred to as low-skill and high-skill periods, respectively, during the summers of 2007–17, and to evaluate the forecast skill of ACs during these respective periods. The ERA-Interim dataset is used to examine characteristics of the Arctic environment and characteristics of ACs during low-skill and high-skill periods. The Arctic environment tends to be characterized by more vigorous baroclinic processes and latent heating during low-skill periods relative to high-skill periods. ACs occur more frequently over much of the Arctic; tend to be stronger; and tend to be located in regions of larger lower-tropospheric baroclinicity, lower-to-midtropospheric Eady growth rate (EGR), and latent heating during low-skill periods relative to high-skill periods. ACs during low-skill periods that are characterized by low forecast skill of intensity tend to be relatively strong and tend to be located in regions of relatively large lower-tropospheric baroclinicity, lower-to-midtropospheric EGR, and latent heating.
Abstract
This study extends initial work by Sun and Penny and Sun et al. to explore the inclusion of path information from surface drifters using an augmented-state Lagrangian data assimilation based on the local ensemble transform Kalman filter (LETKF-LaDA) with vertical localization to improve analysis of the ocean. The region of interest is the Gulf of Mexico during the passage of Hurricane Isaac in the summer of 2012. Results from experiments with a regional ocean model at eddy-permitting and eddy-resolving model resolutions are used to quantify improvements to the analysis of sea surface velocity, sea surface temperature, and sea surface height in a data assimilation system. The data assimilation system assimilates surface drifter positions, as well as vertical profiles of temperature and salinity. Data were used from drifters deployed as a part of the Grand Lagrangian Deployment beginning 20 July 2012. Comparison of experiment results shows that at both eddy-permitting and eddy-resolving horizontal resolutions Lagrangian assimilation of drifter positions significantly improves analysis of the ocean state responding to hurricane conditions. These results, which should be applicable to other tropical oceans such as the Bay of Bengal, open new avenues for estimating ocean initial conditions to improve tropical cyclone forecasting.
Abstract
This study extends initial work by Sun and Penny and Sun et al. to explore the inclusion of path information from surface drifters using an augmented-state Lagrangian data assimilation based on the local ensemble transform Kalman filter (LETKF-LaDA) with vertical localization to improve analysis of the ocean. The region of interest is the Gulf of Mexico during the passage of Hurricane Isaac in the summer of 2012. Results from experiments with a regional ocean model at eddy-permitting and eddy-resolving model resolutions are used to quantify improvements to the analysis of sea surface velocity, sea surface temperature, and sea surface height in a data assimilation system. The data assimilation system assimilates surface drifter positions, as well as vertical profiles of temperature and salinity. Data were used from drifters deployed as a part of the Grand Lagrangian Deployment beginning 20 July 2012. Comparison of experiment results shows that at both eddy-permitting and eddy-resolving horizontal resolutions Lagrangian assimilation of drifter positions significantly improves analysis of the ocean state responding to hurricane conditions. These results, which should be applicable to other tropical oceans such as the Bay of Bengal, open new avenues for estimating ocean initial conditions to improve tropical cyclone forecasting.
Abstract
Deep convection that penetrates the tropopause, referred to here as overshooting convection, is capable of lifting tropospheric air well into the stratosphere. In addition to water, these overshoots also transport various chemical species, affecting chemistry and radiation in the stratosphere. It is not currently known, however, how much transport is a result of this mechanism. To better understand overshooting convection, this study aims to characterize the durations of overshooting events. To achieve this, radar data from the Next Generation Weather Radar (NEXRAD) network is composited onto a three-dimensional grid at 5-min intervals. Overshoots are identified by comparing echo-top heights with tropopause estimates derived from ERA5 reanalysis data. These overshoots are linked in space from one analysis time to the next to form tracks. This process is performed for 12 four-day sample windows in the months May–August of 2017–19. Track characteristics such as duration, overshoot area, tropopause-relative altitude, and column-maximum reflectivity are investigated. Positive correlations are found between track duration and other track characteristics. Integrated track volume is found as a product of the overshoot area, depth, and duration, and provides a measure of the potential stratospheric impact of each track. Short-lived tracks are observed to contribute the most total integrated volume when considering track duration, while tracks that overshoot by 2–3 km show the largest contribution when considering overshoot depth. A diurnal cycle is observed, with peak track initiation around 1600–1700 local time. Track-mean duration peaks a few hours earlier, while track-mean area and tropopause-relative height peak a few hours later.
Abstract
Deep convection that penetrates the tropopause, referred to here as overshooting convection, is capable of lifting tropospheric air well into the stratosphere. In addition to water, these overshoots also transport various chemical species, affecting chemistry and radiation in the stratosphere. It is not currently known, however, how much transport is a result of this mechanism. To better understand overshooting convection, this study aims to characterize the durations of overshooting events. To achieve this, radar data from the Next Generation Weather Radar (NEXRAD) network is composited onto a three-dimensional grid at 5-min intervals. Overshoots are identified by comparing echo-top heights with tropopause estimates derived from ERA5 reanalysis data. These overshoots are linked in space from one analysis time to the next to form tracks. This process is performed for 12 four-day sample windows in the months May–August of 2017–19. Track characteristics such as duration, overshoot area, tropopause-relative altitude, and column-maximum reflectivity are investigated. Positive correlations are found between track duration and other track characteristics. Integrated track volume is found as a product of the overshoot area, depth, and duration, and provides a measure of the potential stratospheric impact of each track. Short-lived tracks are observed to contribute the most total integrated volume when considering track duration, while tracks that overshoot by 2–3 km show the largest contribution when considering overshoot depth. A diurnal cycle is observed, with peak track initiation around 1600–1700 local time. Track-mean duration peaks a few hours earlier, while track-mean area and tropopause-relative height peak a few hours later.
Abstract
Improving lead time for forecasting floods is important to minimize property damage and ensure the safety of the public and emergency services during flood events. Numerical weather prediction (NWP) models are important components of flood forecasting systems and have been vital in extending forecasting lead time under complex weather and terrain conditions. However, NWP forecasts still have significant uncertainty associated with the precipitation fields that are the main inputs of the hydrologic models and thus the resulting flood forecasts. An issue often overlooked is the importance of correctly representing variability over a range of different temporal scales. To address this gap, here a new wavelet-based method for postprocessing NWP precipitation forecasts is proposed. First, precipitation forecasts are decomposed into the frequency domain using a wavelet transform, providing estimates of the amplitudes and phases of the time series at different frequencies. Quantile mapping is then used to correct bias in the amplitudes of each frequency. Randomized phases are used to generate an ensemble of realizations of the precipitation forecasts. The postprocessed precipitation forecasts are reconstructed by taking the inverse of adjusted time-frequency decompositions with the corrected amplitudes and randomized phases. The proposed method was used to postprocess NWP precipitation forecasts in the Sydney region of Australia. There is a significant improvement in postprocessed precipitation forecasts across multiple time scales in terms of bias and temporal and spatial correlation structures. The postprocessed precipitation fields can be used for the modeling of fully distributed hydrologic systems, improving runoff stimulation, flood depth estimation, and flood early warning.
Significance Statement
A new method accounting for the timing and spatial errors of NWP precipitation forecasts is proposed, and it can improve the skill of forecasts across multiple time scales, especially at short lead times. The proposed method provides a practical and effective way to correct these errors by incorporating spatiotemporal neighborhood information through the frequency domain using sophisticated wavelet transforms. With systematic timing and spatial errors removed, precipitation forecasts will be more skillful, and hydrological modeling using the postprocessed forecasts can provide higher accuracy of streamflow estimation.
Abstract
Improving lead time for forecasting floods is important to minimize property damage and ensure the safety of the public and emergency services during flood events. Numerical weather prediction (NWP) models are important components of flood forecasting systems and have been vital in extending forecasting lead time under complex weather and terrain conditions. However, NWP forecasts still have significant uncertainty associated with the precipitation fields that are the main inputs of the hydrologic models and thus the resulting flood forecasts. An issue often overlooked is the importance of correctly representing variability over a range of different temporal scales. To address this gap, here a new wavelet-based method for postprocessing NWP precipitation forecasts is proposed. First, precipitation forecasts are decomposed into the frequency domain using a wavelet transform, providing estimates of the amplitudes and phases of the time series at different frequencies. Quantile mapping is then used to correct bias in the amplitudes of each frequency. Randomized phases are used to generate an ensemble of realizations of the precipitation forecasts. The postprocessed precipitation forecasts are reconstructed by taking the inverse of adjusted time-frequency decompositions with the corrected amplitudes and randomized phases. The proposed method was used to postprocess NWP precipitation forecasts in the Sydney region of Australia. There is a significant improvement in postprocessed precipitation forecasts across multiple time scales in terms of bias and temporal and spatial correlation structures. The postprocessed precipitation fields can be used for the modeling of fully distributed hydrologic systems, improving runoff stimulation, flood depth estimation, and flood early warning.
Significance Statement
A new method accounting for the timing and spatial errors of NWP precipitation forecasts is proposed, and it can improve the skill of forecasts across multiple time scales, especially at short lead times. The proposed method provides a practical and effective way to correct these errors by incorporating spatiotemporal neighborhood information through the frequency domain using sophisticated wavelet transforms. With systematic timing and spatial errors removed, precipitation forecasts will be more skillful, and hydrological modeling using the postprocessed forecasts can provide higher accuracy of streamflow estimation.
Abstract
Forecasting mesoscale convective systems (MCSs) and precipitation over complex terrain is an ongoing challenge even for convective-permitting numerical models. Here, we show the value of combining mesoscale constraints to improve short-term MCS forecasts for two events during the North American monsoon season in 2013, including the following: 1) the initial specification of moisture, via GPS-precipitable water vapor (PWV) data assimilation (DA); 2) kinematics via modification of cumulus parameterization; and 3) microphysics via modification of cloud microphysics parameterization. A total of five convective-permitting Weather Research and Forecasting (WRF) Model experiments is conducted for each event to elucidate the impact of these constraints. Results show that combining GPS-PWV DA with a modified Kain–Fritsch scheme and double-moment microphysics provides relatively the best forecast of both North American monsoon MCSs and convective precipitation in terms of timing, location, and intensity relative to available precipitation and cloud-top temperature observations. Additional examination on the associated reflectivity, vertical wind field, equivalent potential temperature, and hydrometeor distribution of MCS events show the added value of each individual constraint to forecast performance.
Significance Statement
Forecasting thunderstorm clouds and rain over mountainous regions is challenging because of limitations in having radar and rain gauges and in resolving physical drivers in forecast models. We examine the value of considering all possible constraints by incorporating moisture into these models, and correcting physics in the model treatment of cumulus and cloud microphysics parameterizations. This study demonstrates that assimilating moisture and using modified Kain–Fritsch and double-moment microphysics schemes provides the best thunderstorm cloud and rain forecasts in terms of timing, location, and intensity. Each correction improves key properties of these storms such as vertical wind, along with distribution of water in various phases. We highlight the need to improve our efforts on effectively integrating these constraints into current and future forecasts.
Abstract
Forecasting mesoscale convective systems (MCSs) and precipitation over complex terrain is an ongoing challenge even for convective-permitting numerical models. Here, we show the value of combining mesoscale constraints to improve short-term MCS forecasts for two events during the North American monsoon season in 2013, including the following: 1) the initial specification of moisture, via GPS-precipitable water vapor (PWV) data assimilation (DA); 2) kinematics via modification of cumulus parameterization; and 3) microphysics via modification of cloud microphysics parameterization. A total of five convective-permitting Weather Research and Forecasting (WRF) Model experiments is conducted for each event to elucidate the impact of these constraints. Results show that combining GPS-PWV DA with a modified Kain–Fritsch scheme and double-moment microphysics provides relatively the best forecast of both North American monsoon MCSs and convective precipitation in terms of timing, location, and intensity relative to available precipitation and cloud-top temperature observations. Additional examination on the associated reflectivity, vertical wind field, equivalent potential temperature, and hydrometeor distribution of MCS events show the added value of each individual constraint to forecast performance.
Significance Statement
Forecasting thunderstorm clouds and rain over mountainous regions is challenging because of limitations in having radar and rain gauges and in resolving physical drivers in forecast models. We examine the value of considering all possible constraints by incorporating moisture into these models, and correcting physics in the model treatment of cumulus and cloud microphysics parameterizations. This study demonstrates that assimilating moisture and using modified Kain–Fritsch and double-moment microphysics schemes provides the best thunderstorm cloud and rain forecasts in terms of timing, location, and intensity. Each correction improves key properties of these storms such as vertical wind, along with distribution of water in various phases. We highlight the need to improve our efforts on effectively integrating these constraints into current and future forecasts.
Abstract
This paper examines the meteorological factors that led to the record-breaking heavy precipitation event in Taiwan in the early 2020 mei-yu season (15–31 May). The extreme amount of rainfall (an average of 135.9 mm per station) during the 36-h period around 22 May (hereafter Y20R) also set a record. Compared to climatology, the Pacific subtropical high was stronger and the southwesterly monsoonal flow was more intense during the first half of the 2020 mei-yu season, resulting in a stronger moisture conveyor belt over the northern Indo-China Peninsula. The record-breaking precipitation in Y20R was mainly caused by the eastward movement of a southwest vortex (SWV) generated in southwestern China. When the eastern portion of the SWV touched northern Taiwan, its associated west-southwesterly winds and the large-scale southwesterly monsoonal flow transported moisture toward the Taiwan Strait. The moisture-laden southwesterly flow was lifted by the stationary mei-yu front, leading to the heavy rainfall in southern Taiwan. When the SWV passed through northern Taiwan, it became the dominant weather system that enhanced the west-southwesterly winds and transported moisture from South China to Taiwan. The front moved southward through the Taiwan Strait during this period, with its location greatly determining the pattern of rainfall in southern Taiwan. In summary, the most critical factors leading to heavy rainfall in southern Taiwan are the strong 850-hPa southwesterly winds and moisture fluxes associated with the SWV. The other key factors include, in order of sensitivity to rainfall, the distance of the front, the distance of the SWV, the frontal speed, and the intensity of the SWV.
Abstract
This paper examines the meteorological factors that led to the record-breaking heavy precipitation event in Taiwan in the early 2020 mei-yu season (15–31 May). The extreme amount of rainfall (an average of 135.9 mm per station) during the 36-h period around 22 May (hereafter Y20R) also set a record. Compared to climatology, the Pacific subtropical high was stronger and the southwesterly monsoonal flow was more intense during the first half of the 2020 mei-yu season, resulting in a stronger moisture conveyor belt over the northern Indo-China Peninsula. The record-breaking precipitation in Y20R was mainly caused by the eastward movement of a southwest vortex (SWV) generated in southwestern China. When the eastern portion of the SWV touched northern Taiwan, its associated west-southwesterly winds and the large-scale southwesterly monsoonal flow transported moisture toward the Taiwan Strait. The moisture-laden southwesterly flow was lifted by the stationary mei-yu front, leading to the heavy rainfall in southern Taiwan. When the SWV passed through northern Taiwan, it became the dominant weather system that enhanced the west-southwesterly winds and transported moisture from South China to Taiwan. The front moved southward through the Taiwan Strait during this period, with its location greatly determining the pattern of rainfall in southern Taiwan. In summary, the most critical factors leading to heavy rainfall in southern Taiwan are the strong 850-hPa southwesterly winds and moisture fluxes associated with the SWV. The other key factors include, in order of sensitivity to rainfall, the distance of the front, the distance of the SWV, the frontal speed, and the intensity of the SWV.
Abstract
The national upgrade of the operational weather radar network to include polarimetric capabilities has led to numerous studies focusing on polarimetric radar signatures commonly observed in supercells. One such signature is the horizontal separation of regions of enhanced differential reflectivity (Z DR) and specific differential phase (K DP) values due to hydrometeor size sorting. Recent observational studies have shown that the orientation of this separation tends to be more perpendicular to storm motion in supercells that produce tornadoes. Although this finding has potential operational utility, the physical relationship between this observed radar signature and tornadic potential is not known. This study uses an ensemble of supercell simulations initialized with tornadic and nontornadic environments to investigate this connection. The tendency for tornadic supercells to have a more perpendicular separation orientation was reproduced, although to a lesser degree. This difference in orientation angles was caused by stronger rearward storm-relative flow in the nontornadic supercells, leading to a rearward shift of precipitation and, therefore, the enhanced K DP region within the supercell. Further, this resulted in an unfavorable rearward shift of the negative buoyancy region, which led to an order of magnitude less baroclinic generation of circulation in the nontornadic simulations compared to tornadic simulations.
Abstract
The national upgrade of the operational weather radar network to include polarimetric capabilities has led to numerous studies focusing on polarimetric radar signatures commonly observed in supercells. One such signature is the horizontal separation of regions of enhanced differential reflectivity (Z DR) and specific differential phase (K DP) values due to hydrometeor size sorting. Recent observational studies have shown that the orientation of this separation tends to be more perpendicular to storm motion in supercells that produce tornadoes. Although this finding has potential operational utility, the physical relationship between this observed radar signature and tornadic potential is not known. This study uses an ensemble of supercell simulations initialized with tornadic and nontornadic environments to investigate this connection. The tendency for tornadic supercells to have a more perpendicular separation orientation was reproduced, although to a lesser degree. This difference in orientation angles was caused by stronger rearward storm-relative flow in the nontornadic supercells, leading to a rearward shift of precipitation and, therefore, the enhanced K DP region within the supercell. Further, this resulted in an unfavorable rearward shift of the negative buoyancy region, which led to an order of magnitude less baroclinic generation of circulation in the nontornadic simulations compared to tornadic simulations.
Abstract
This study evaluates the impact of assimilating Global Navigation Satellite System (GNSS) radio occultation (RO) bending angles from Formosa Satellite Mission-7/Constellation Observing System for Meteorology, Ionosphere, and Climate-2 (COSMIC-2) receiver satellites on Hurricane Weather Research and Forecasting (HWRF) Model tropical cyclone (TC) forecasts. Launched in June 2019, the COSMIC-2 mission provides significantly higher tropics data coverage compared to its predecessor COSMIC constellation. GNSS RO measurements yield information about atmospheric pressure, temperature, and water vapor profiles. HWRF is cycled with and without COSMIC-2 bending angle data assimilation for six 2020 Atlantic hurricane cases. COSMIC-2 assimilation has little impact on HWRF track forecasts, consistent with HWRF’s design limiting cycled data assimilation impacts on surrounding large-scale flows; however, COSMIC-2 assimilation results in a statistically significant ∼8%–12% mean absolute forecast error reduction in minimum central sea level pressure for t = 36-, 54-, 60-, and 108–120-h lead times. Forecasts initialized from analyses assimilating COSMIC-2 observations also have a 1%–4% smaller 600–700-hPa specific humidity (SPFH) root-mean-squared deviation compared to radiosondes and dropwindsondes for most lead times. While not all HWRF intensity forecasts benefit from COSMIC-2 assimilation, a few show notable improvement. For example, assimilating two COSMIC-2 profiles within the inner core of developing Hurricane Hanna (2020) increases 800-hPa SPFH by up to 1 g kg−1 locally, helping to correct a dry bias. The forecast initialized from this analysis better captures Hanna’s observed intensification rate, likely because its moister inner core facilitates development of persistent deep convection near the TC center, where diabatic heating is more efficiently converted to cyclonic wind kinetic energy.
Significance Statement
Tropical cyclone (TC) intensification can be strongly sensitive to the lower-to-midtropospheric water vapor distribution near the storm. The COSMIC-2 GNSS radio occultation (RO) receiver satellite mission provides denser spatial coverage of atmospheric water vapor and temperature profiles over the tropics compared to other GNSS RO observation platforms. Herein, using six 2020 Atlantic TC cases, we evaluate the impacts of assimilating COSMIC-2 RO bending angles into a regional forecast model that already assimilates clear-sky satellite radiances. It is shown that COSMIC-2 assimilation yields a modest ∼10% intensity forecast skill improvement for several lead times, although more substantial intensity forecast improvement is found for a few forecasts where the COSMIC-2 observation assimilation helps correct a lower-to-midtropospheric water vapor bias.
Abstract
This study evaluates the impact of assimilating Global Navigation Satellite System (GNSS) radio occultation (RO) bending angles from Formosa Satellite Mission-7/Constellation Observing System for Meteorology, Ionosphere, and Climate-2 (COSMIC-2) receiver satellites on Hurricane Weather Research and Forecasting (HWRF) Model tropical cyclone (TC) forecasts. Launched in June 2019, the COSMIC-2 mission provides significantly higher tropics data coverage compared to its predecessor COSMIC constellation. GNSS RO measurements yield information about atmospheric pressure, temperature, and water vapor profiles. HWRF is cycled with and without COSMIC-2 bending angle data assimilation for six 2020 Atlantic hurricane cases. COSMIC-2 assimilation has little impact on HWRF track forecasts, consistent with HWRF’s design limiting cycled data assimilation impacts on surrounding large-scale flows; however, COSMIC-2 assimilation results in a statistically significant ∼8%–12% mean absolute forecast error reduction in minimum central sea level pressure for t = 36-, 54-, 60-, and 108–120-h lead times. Forecasts initialized from analyses assimilating COSMIC-2 observations also have a 1%–4% smaller 600–700-hPa specific humidity (SPFH) root-mean-squared deviation compared to radiosondes and dropwindsondes for most lead times. While not all HWRF intensity forecasts benefit from COSMIC-2 assimilation, a few show notable improvement. For example, assimilating two COSMIC-2 profiles within the inner core of developing Hurricane Hanna (2020) increases 800-hPa SPFH by up to 1 g kg−1 locally, helping to correct a dry bias. The forecast initialized from this analysis better captures Hanna’s observed intensification rate, likely because its moister inner core facilitates development of persistent deep convection near the TC center, where diabatic heating is more efficiently converted to cyclonic wind kinetic energy.
Significance Statement
Tropical cyclone (TC) intensification can be strongly sensitive to the lower-to-midtropospheric water vapor distribution near the storm. The COSMIC-2 GNSS radio occultation (RO) receiver satellite mission provides denser spatial coverage of atmospheric water vapor and temperature profiles over the tropics compared to other GNSS RO observation platforms. Herein, using six 2020 Atlantic TC cases, we evaluate the impacts of assimilating COSMIC-2 RO bending angles into a regional forecast model that already assimilates clear-sky satellite radiances. It is shown that COSMIC-2 assimilation yields a modest ∼10% intensity forecast skill improvement for several lead times, although more substantial intensity forecast improvement is found for a few forecasts where the COSMIC-2 observation assimilation helps correct a lower-to-midtropospheric water vapor bias.
