Abstract
Severe rainfall has become increasingly frequent and intense in the Taipei metropolitan area. A complex thunderstorm in the Taipei Basin on 14 June 2015 produced an extreme rain rate (>130 mm h−1), leading to an urban flash flood. This paper presents storms’ microphysical and dynamic features during the organizing and heavy rain stages, mainly based on observed polarimetric variables in a Doppler radar network and ground-based raindrop size distribution. Shallower isolated cells in the early afternoon characterized by big raindrops produced a rain rate > 10 mm h−1, but the rain showers persisted for a short time. The storm’s evolution highlighted the behavior of merged convective cells before the heaviest rainfall (exceeding 60 mm within 20 min). The columnar features of differential reflectivity (Z DR) and specific differential phase (K DP) became more evident in merged cells, which correlated with the broad distribution of upward motion and mixed-phase hydrometeors. The K DP below the environmental 0°C level increased toward the ground associated with the melted graupel and resulted in subsequent intense rain rates, showing the contribution of the ice-phase process. Due to the collision–breakup process, the highest concentrations of almost all drop sizes and smaller mass-weighted mean diameter occurred during the maximum rainfall stage.
Abstract
Severe rainfall has become increasingly frequent and intense in the Taipei metropolitan area. A complex thunderstorm in the Taipei Basin on 14 June 2015 produced an extreme rain rate (>130 mm h−1), leading to an urban flash flood. This paper presents storms’ microphysical and dynamic features during the organizing and heavy rain stages, mainly based on observed polarimetric variables in a Doppler radar network and ground-based raindrop size distribution. Shallower isolated cells in the early afternoon characterized by big raindrops produced a rain rate > 10 mm h−1, but the rain showers persisted for a short time. The storm’s evolution highlighted the behavior of merged convective cells before the heaviest rainfall (exceeding 60 mm within 20 min). The columnar features of differential reflectivity (Z DR) and specific differential phase (K DP) became more evident in merged cells, which correlated with the broad distribution of upward motion and mixed-phase hydrometeors. The K DP below the environmental 0°C level increased toward the ground associated with the melted graupel and resulted in subsequent intense rain rates, showing the contribution of the ice-phase process. Due to the collision–breakup process, the highest concentrations of almost all drop sizes and smaller mass-weighted mean diameter occurred during the maximum rainfall stage.
Abstract
The second part of this series presents results from verifying a precipitation forecast calibration method discussed in part 1, based on quantile mapping (QM), weighting of sorted members, and dressing of the ensemble. NOAA’s Global Ensemble Forecast System, version 12 (GEFSv12) reforecasts were used in this study. The method was validated with pre-operational GEFSv12 forecasts from the between December 2017 and November 2019. The method is proposed as an enhancement for GEFSv12 precipitation postprocessing in NOAA’s National Blend of Models.
Part 1 described adaptations to the methodology to leverage the ~ 20-year GEFSv12 reforecast data. As shown in this part 2, when compared to probabilistic quantitative precipitation forecasts (PQPFs) from the raw ensemble, the adapted method produced downscaled, high-resolution forecasts that were significantly more reliable and skillful than raw ensemble-derived probabilities, especially at shorter lead times (i.e., < 5 days) and for forecasts of events from light precipitation to > 10 mm 6 h−1. Cool-season events in the western US were especially improved when the QM algorithm was applied, providing a statistical downscaling with realistic smaller-scale detail related to terrain features. The method provided less value added for forecasts of longer lead times and for the heaviest precipitation.
Abstract
The second part of this series presents results from verifying a precipitation forecast calibration method discussed in part 1, based on quantile mapping (QM), weighting of sorted members, and dressing of the ensemble. NOAA’s Global Ensemble Forecast System, version 12 (GEFSv12) reforecasts were used in this study. The method was validated with pre-operational GEFSv12 forecasts from the between December 2017 and November 2019. The method is proposed as an enhancement for GEFSv12 precipitation postprocessing in NOAA’s National Blend of Models.
Part 1 described adaptations to the methodology to leverage the ~ 20-year GEFSv12 reforecast data. As shown in this part 2, when compared to probabilistic quantitative precipitation forecasts (PQPFs) from the raw ensemble, the adapted method produced downscaled, high-resolution forecasts that were significantly more reliable and skillful than raw ensemble-derived probabilities, especially at shorter lead times (i.e., < 5 days) and for forecasts of events from light precipitation to > 10 mm 6 h−1. Cool-season events in the western US were especially improved when the QM algorithm was applied, providing a statistical downscaling with realistic smaller-scale detail related to terrain features. The method provided less value added for forecasts of longer lead times and for the heaviest precipitation.
Abstract
Data from a dual-polarized, solid-state X-band radar and an operational C-band weather radar are used for high-resolution analyses of two hailstorms in the Vienna, Austria, region. The combination of both radars provides rapid-update (1 min) polarimetric data paired with wind field data of a dual-Doppler analysis. This is the first time that such an advanced setup is used to examine severe storm dynamics at the eastern Alpine fringe, where the influence of local topography is particularly challenging for thunderstorm prediction. We investigate two storms transitioning from the pre-Alps into the Vienna basin with different characteristics: 1) A rapidly evolving multicell storm producing large hail (5 cm), with observations of an intense Z DR column preceding hail formation and the rapid development of multiple pulses of hail; and 2) a cold pool–driven squall line with small hail, for which we find that the updraft location inhibited the formation of larger hailstones. For both cases, we analyzed the evolution of different Z DR column metrics as well as updraft speed and size and found that (i) the 90th percentile of Z DR within the Z DR column was highest for the cell later producing large hail, (ii) the peak 90th percentile of Z DR preceded large hailfall by 20 min and highest updraft size and speed by 10 min, and (iii) sudden drops of the 90th percentile of ZH within the Z DR column indicated imminent hailfall.
Significance Statement
Thunderstorm evolution on the transition from complex terrain into the Vienna basin in northeastern Austria varies strongly. In some instances, thunderstorm cells intensify once they reach flat terrain, while in most cases there is a weakening tendency. To improve our process understanding and short-term forecasting methods, we analyze two representative cases of hail-bearing storms transitioning into the Vienna basin. We mainly build our study on data from a new, cost-efficient weather radar, complemented by an operational radar, lightning observations, and ground reports. Our results show which radar variables could be well suited for early detection of intensification, and how they relate to thunderstorm updraft speeds and lightning activity.
Abstract
Data from a dual-polarized, solid-state X-band radar and an operational C-band weather radar are used for high-resolution analyses of two hailstorms in the Vienna, Austria, region. The combination of both radars provides rapid-update (1 min) polarimetric data paired with wind field data of a dual-Doppler analysis. This is the first time that such an advanced setup is used to examine severe storm dynamics at the eastern Alpine fringe, where the influence of local topography is particularly challenging for thunderstorm prediction. We investigate two storms transitioning from the pre-Alps into the Vienna basin with different characteristics: 1) A rapidly evolving multicell storm producing large hail (5 cm), with observations of an intense Z DR column preceding hail formation and the rapid development of multiple pulses of hail; and 2) a cold pool–driven squall line with small hail, for which we find that the updraft location inhibited the formation of larger hailstones. For both cases, we analyzed the evolution of different Z DR column metrics as well as updraft speed and size and found that (i) the 90th percentile of Z DR within the Z DR column was highest for the cell later producing large hail, (ii) the peak 90th percentile of Z DR preceded large hailfall by 20 min and highest updraft size and speed by 10 min, and (iii) sudden drops of the 90th percentile of ZH within the Z DR column indicated imminent hailfall.
Significance Statement
Thunderstorm evolution on the transition from complex terrain into the Vienna basin in northeastern Austria varies strongly. In some instances, thunderstorm cells intensify once they reach flat terrain, while in most cases there is a weakening tendency. To improve our process understanding and short-term forecasting methods, we analyze two representative cases of hail-bearing storms transitioning into the Vienna basin. We mainly build our study on data from a new, cost-efficient weather radar, complemented by an operational radar, lightning observations, and ground reports. Our results show which radar variables could be well suited for early detection of intensification, and how they relate to thunderstorm updraft speeds and lightning activity.
Abstract
This study revisits the issue of why tropical cyclones (TCs) develop more rapidly at lower latitudes, using ensemble axisymmetric numerical simulations and energy diagnostics based on the isentropic analysis, with the focus on the relative importance of the outflow-layer and boundary-layer inertial stabilities to TC intensification and energy cycle. Results show that although lowering the outflow-layer Coriolis parameter and thus inertial stability can slightly strengthen the outflow, it does not affect the simulated TC development, whereas lowering the boundary-layer Coriolis parameter largely enhances the secondary circulation and TC intensification as in the experiment with a reduced Coriolis parameter throughout the model atmosphere. This suggests that TC outflow is more likely a passive result of the convergent inflow in the boundary layer and convective updraft in the eyewall.
The boundary-layer inertial stability is found to control the convergent inflow in the boundary layer and depth of convection in the eyewall and thus the temperature of energy sink in the TC heat engine, which determines the efficiency and overall mechanical output of heat engine and thus TC intensification. It is also shown that the hypothesized isothermal and adiabatic compression legs at the downstream end of the outflow in the classical Carnot cycle is not supported in the thermodynamic cycle of the simulated TCs, implying that the assumed TC Carnot cycle is not closed. It is the theoretical maximum work of heat engine, not the energy expenditure following the outflow downstream, that determines the mechanical work used to intensify a TC.
Abstract
This study revisits the issue of why tropical cyclones (TCs) develop more rapidly at lower latitudes, using ensemble axisymmetric numerical simulations and energy diagnostics based on the isentropic analysis, with the focus on the relative importance of the outflow-layer and boundary-layer inertial stabilities to TC intensification and energy cycle. Results show that although lowering the outflow-layer Coriolis parameter and thus inertial stability can slightly strengthen the outflow, it does not affect the simulated TC development, whereas lowering the boundary-layer Coriolis parameter largely enhances the secondary circulation and TC intensification as in the experiment with a reduced Coriolis parameter throughout the model atmosphere. This suggests that TC outflow is more likely a passive result of the convergent inflow in the boundary layer and convective updraft in the eyewall.
The boundary-layer inertial stability is found to control the convergent inflow in the boundary layer and depth of convection in the eyewall and thus the temperature of energy sink in the TC heat engine, which determines the efficiency and overall mechanical output of heat engine and thus TC intensification. It is also shown that the hypothesized isothermal and adiabatic compression legs at the downstream end of the outflow in the classical Carnot cycle is not supported in the thermodynamic cycle of the simulated TCs, implying that the assumed TC Carnot cycle is not closed. It is the theoretical maximum work of heat engine, not the energy expenditure following the outflow downstream, that determines the mechanical work used to intensify a TC.
Abstract
The atmospheric circulation response to global warming is an important problem that is theoretically still not well understood. This is a particular issue since climate model simulations provide uncertain, and at times contradictory, projections of future climate. In particular, it is still unclear how a warmer and moister atmosphere will affect midlatitude eddies and their associated poleward transport of heat and moisture. Here we perform a trend analysis of three main components of the global circulation—the zonal-mean state, eddies, and the net energy input into the atmosphere—and examine how they relate in terms of a moist static energy budget for the JRA-55 reanalysis data. A particular emphasis is made on understanding the contribution of moisture to circulation trends. The observed trends are very different between the hemispheres. In the Southern Hemisphere there is an overall strengthening and during boreal summer, also a poleward shifting, of the jet stream, the eddies, and the meridional diabatic heating gradients. Correspondingly, we find an overall strengthening of the meridional gradients of the net atmospheric energy input. In the Northern Hemisphere, the trend patterns are more complex, with the dominant signal being a clear boreal winter Arctic amplification of positive trends in lower-tropospheric temperature and moisture, as well as a significant weakening of both bandpass and low-pass eddy heat and moisture fluxes. Consistently, surface latent and sensible heat fluxes, upward and downward longwave radiation, and longwave cloud radiative fluxes at high latitudes show significant trends. However, radiative fluxes and eddy fluxes are inconsistent, suggesting data assimilation procedures need to be improved.
Significance Statement
We use a long-term reanalysis dataset to get an overall view of the changes in the global circulation and its role in transporting moist static energy from the equator to the poles. We do this by examining the trends in its three main components—the zonal means, the eddies, and the net energy input into the atmosphere. We find that in the Southern Hemisphere, there is an overall strengthening of the eddies, their poleward energy fluxes, and correspondingly the meridional gradients of the net atmospheric energy input. In the Northern Hemisphere, though the pattern is more complex, there is an overall weakening of the eddies and poleward eddy fluxes, and of the meridional gradients of the net atmospheric energy input, consistent with Arctic warming.
Abstract
The atmospheric circulation response to global warming is an important problem that is theoretically still not well understood. This is a particular issue since climate model simulations provide uncertain, and at times contradictory, projections of future climate. In particular, it is still unclear how a warmer and moister atmosphere will affect midlatitude eddies and their associated poleward transport of heat and moisture. Here we perform a trend analysis of three main components of the global circulation—the zonal-mean state, eddies, and the net energy input into the atmosphere—and examine how they relate in terms of a moist static energy budget for the JRA-55 reanalysis data. A particular emphasis is made on understanding the contribution of moisture to circulation trends. The observed trends are very different between the hemispheres. In the Southern Hemisphere there is an overall strengthening and during boreal summer, also a poleward shifting, of the jet stream, the eddies, and the meridional diabatic heating gradients. Correspondingly, we find an overall strengthening of the meridional gradients of the net atmospheric energy input. In the Northern Hemisphere, the trend patterns are more complex, with the dominant signal being a clear boreal winter Arctic amplification of positive trends in lower-tropospheric temperature and moisture, as well as a significant weakening of both bandpass and low-pass eddy heat and moisture fluxes. Consistently, surface latent and sensible heat fluxes, upward and downward longwave radiation, and longwave cloud radiative fluxes at high latitudes show significant trends. However, radiative fluxes and eddy fluxes are inconsistent, suggesting data assimilation procedures need to be improved.
Significance Statement
We use a long-term reanalysis dataset to get an overall view of the changes in the global circulation and its role in transporting moist static energy from the equator to the poles. We do this by examining the trends in its three main components—the zonal means, the eddies, and the net energy input into the atmosphere. We find that in the Southern Hemisphere, there is an overall strengthening of the eddies, their poleward energy fluxes, and correspondingly the meridional gradients of the net atmospheric energy input. In the Northern Hemisphere, though the pattern is more complex, there is an overall weakening of the eddies and poleward eddy fluxes, and of the meridional gradients of the net atmospheric energy input, consistent with Arctic warming.
Abstract
Reconstructing the history of polar temperature from ice core water isotope (δ 18O) calibration has remained a challenge in paleoclimate research, because of our incomplete understanding of various temperature–δ 18O relationships. This paper resolves this classical problem in a new framework called the unified slope equations (USE), which illustrates the general relations among spatial and temporal δ 18O–surface temperature slopes. The USE is applied to the Antarctica temperature change during the last deglaciation in model simulations and observations. It is shown that the comparable Antarctica-mean spatial slope with deglacial temporal slope in δ 18O–surface temperature reconstruction is caused, accidentally, by the compensation responses between the δ 18O–inversion layer temperature relation and the inversion layer temperature itself. Furthermore, in light of the USE, we propose that the present seasonal slope of δ 18O–inversion layer temperature is an optimal paleothermometer that is more accurate and robust than the spatial slope. This optimal slope suggests the possibility of reconstructing past Antarctic temperature changes using present and future instrumental observations.
Significance Statement
This paper develops a new framework called the unified slope equations (USE) to provide, for the first time, a general relation among various spatial and temporal water isotope–temperature slopes. The application of the USE to Antarctic deglacial temperature change shows that the optimal paleothermometer is the seasonal slope of the inversion layer temperature.
Abstract
Reconstructing the history of polar temperature from ice core water isotope (δ 18O) calibration has remained a challenge in paleoclimate research, because of our incomplete understanding of various temperature–δ 18O relationships. This paper resolves this classical problem in a new framework called the unified slope equations (USE), which illustrates the general relations among spatial and temporal δ 18O–surface temperature slopes. The USE is applied to the Antarctica temperature change during the last deglaciation in model simulations and observations. It is shown that the comparable Antarctica-mean spatial slope with deglacial temporal slope in δ 18O–surface temperature reconstruction is caused, accidentally, by the compensation responses between the δ 18O–inversion layer temperature relation and the inversion layer temperature itself. Furthermore, in light of the USE, we propose that the present seasonal slope of δ 18O–inversion layer temperature is an optimal paleothermometer that is more accurate and robust than the spatial slope. This optimal slope suggests the possibility of reconstructing past Antarctic temperature changes using present and future instrumental observations.
Significance Statement
This paper develops a new framework called the unified slope equations (USE) to provide, for the first time, a general relation among various spatial and temporal water isotope–temperature slopes. The application of the USE to Antarctic deglacial temperature change shows that the optimal paleothermometer is the seasonal slope of the inversion layer temperature.
Abstract
Statistical downscaling (SD) of climate change projections is a key piece for impact and adaptation studies, due to its low computational expense compared to dynamical downscaling, which allows to explore uncertainties through the generation of large ensembles. SD has been extensively evaluated and applied in the extratropics, but few examples exist in tropical regions. In this study several state-of-the-art methods belonging to different families have been evaluated for maximum/minimum daily temperature and daily accumulated precipitation (both from the ERA5 reanalysis at 0.25°) in two regions with very different climates: Spain (Mid-latitudes) and Central America (Tropics). Some key assumptions of SD have been tested: the strength of the predictors/predictand links, the skill of different approaches and the extrapolation capability of each method. It has been found that relevant predictors are different in both regions, as well as the behavior of statistical methods. For temperature, most methods perform significantly better in Spain than in Central America, where Transfer Function methods present important extrapolation problems, probably due to the low variability of the training sample (present climate). In both regions, Model Output Statistics (MOS) methods have achieved the best results for temperature. In Central America Transfer Function (TF) methods have achieved better results than MOS methods in the evaluation in the present climate, but they do not preserve trends in the future. For precipitation, MOS methods and the machine learning method eXtreme Gradient Boost have achieved the best results in both regions. Additionally, it has been found that although the use of humidity indexes as predictors improve results for the downscaling of precipitation, future trends given by statistical methods are very sensitive to the use of one or another index. Three indexes have been compared: relative humidity, specific humidity and dew point depression. The use of the specific humidity has been found to seriously deviate trends given by the downscaled projections from those given by raw Global Climate Models in both regions.
Abstract
Statistical downscaling (SD) of climate change projections is a key piece for impact and adaptation studies, due to its low computational expense compared to dynamical downscaling, which allows to explore uncertainties through the generation of large ensembles. SD has been extensively evaluated and applied in the extratropics, but few examples exist in tropical regions. In this study several state-of-the-art methods belonging to different families have been evaluated for maximum/minimum daily temperature and daily accumulated precipitation (both from the ERA5 reanalysis at 0.25°) in two regions with very different climates: Spain (Mid-latitudes) and Central America (Tropics). Some key assumptions of SD have been tested: the strength of the predictors/predictand links, the skill of different approaches and the extrapolation capability of each method. It has been found that relevant predictors are different in both regions, as well as the behavior of statistical methods. For temperature, most methods perform significantly better in Spain than in Central America, where Transfer Function methods present important extrapolation problems, probably due to the low variability of the training sample (present climate). In both regions, Model Output Statistics (MOS) methods have achieved the best results for temperature. In Central America Transfer Function (TF) methods have achieved better results than MOS methods in the evaluation in the present climate, but they do not preserve trends in the future. For precipitation, MOS methods and the machine learning method eXtreme Gradient Boost have achieved the best results in both regions. Additionally, it has been found that although the use of humidity indexes as predictors improve results for the downscaling of precipitation, future trends given by statistical methods are very sensitive to the use of one or another index. Three indexes have been compared: relative humidity, specific humidity and dew point depression. The use of the specific humidity has been found to seriously deviate trends given by the downscaled projections from those given by raw Global Climate Models in both regions.
Abstract
Based on the principle “learn from past errors to correct current forecasts,” statistical postprocessing consists of optimizing forecasts generated by numerical weather prediction (NWP) models. In this context, machine learning (ML) offers state-of-the-art tools for training statistical models and making predictions based on large datasets. In our study, ML-based solutions are developed to reduce forecast errors of 2-m temperature and 10-m wind speed of the ECMWF’s operational medium-range, high-resolution forecasts produced with the Integrated Forecasting System (IFS). IFS forecasts and other spatiotemporal indicators are used as predictors after careful selection with the help of ML interpretability tools. Different ML approaches are tested: linear regression, random forest decision trees, and neural networks. Statistical models of systematic and random errors are derived sequentially where the random error is defined as the residual error after bias correction. In terms of output, bias correction and forecast uncertainty prediction are made available at any point from locations around the world. All three ML methods show a similar ability to capture situation-dependent biases leading to noteworthy performance improvements (between 10% and 15% improvement in terms of root-mean-square error for all lead times and variables), and a similar ability to provide reliable uncertainty predictions.
Abstract
Based on the principle “learn from past errors to correct current forecasts,” statistical postprocessing consists of optimizing forecasts generated by numerical weather prediction (NWP) models. In this context, machine learning (ML) offers state-of-the-art tools for training statistical models and making predictions based on large datasets. In our study, ML-based solutions are developed to reduce forecast errors of 2-m temperature and 10-m wind speed of the ECMWF’s operational medium-range, high-resolution forecasts produced with the Integrated Forecasting System (IFS). IFS forecasts and other spatiotemporal indicators are used as predictors after careful selection with the help of ML interpretability tools. Different ML approaches are tested: linear regression, random forest decision trees, and neural networks. Statistical models of systematic and random errors are derived sequentially where the random error is defined as the residual error after bias correction. In terms of output, bias correction and forecast uncertainty prediction are made available at any point from locations around the world. All three ML methods show a similar ability to capture situation-dependent biases leading to noteworthy performance improvements (between 10% and 15% improvement in terms of root-mean-square error for all lead times and variables), and a similar ability to provide reliable uncertainty predictions.
Abstract
As a part of the Tropical Cyclone Rapid Intensification Project (TCRI), observations were made of the rapid intensification of Hurricane Sally (2020) as it passed over the Gulf of Mexico. High-altitude dropsondes and radar observations from NOAA’s Gulfstream IV, radar observations from WP-3D aircraft, the WSR-88D ground radar network, satellite images, and satellite-detected lightning strikes are used to apply recently developed theoretical knowledge about tropical cyclone intensification. As observed in many other tropical cyclones, strong, bottom-heavy vertical mass flux profiles are correlated with low (but positive) values of low- to midlevel moist convective instability along with high column relative humidity. Such mass flux profiles produce rapid spinup at low levels and the environmental conditions giving rise to them are associated with an intense midlevel vortex. This low-level spinup underneath the midlevel vortex results in the vertical alignment of the vortex column, which is a key step in the rapid intensification process. In the case of Sally, the spinup of the low-level vortex resulted from vorticity stretching, while the spinup of the midlevel vortex at 6 km resulted from vorticity tilting produced by the interaction of convective ascent with moderate vertical shear.
Significance Statement
The purpose of this study is to investigate the rapid intensification of Hurricane Sally as it was approaching the Florida Panhandle. We do that by analyzing an unprecedented dataset from the NOAA WP-3D and Gulfstream-IV aircraft, together with ground-based radar and satellite data. We find that both the dynamics (vorticity structure and evolution) and thermodynamics (instability index, saturation fraction, heating/mass flux profiles) need to be considered in diagnosing intensification processes. Further field projects with continuous high-altitude dropsondes and research are needed to see if these are applicable to other reformation events as well as genesis.
Abstract
As a part of the Tropical Cyclone Rapid Intensification Project (TCRI), observations were made of the rapid intensification of Hurricane Sally (2020) as it passed over the Gulf of Mexico. High-altitude dropsondes and radar observations from NOAA’s Gulfstream IV, radar observations from WP-3D aircraft, the WSR-88D ground radar network, satellite images, and satellite-detected lightning strikes are used to apply recently developed theoretical knowledge about tropical cyclone intensification. As observed in many other tropical cyclones, strong, bottom-heavy vertical mass flux profiles are correlated with low (but positive) values of low- to midlevel moist convective instability along with high column relative humidity. Such mass flux profiles produce rapid spinup at low levels and the environmental conditions giving rise to them are associated with an intense midlevel vortex. This low-level spinup underneath the midlevel vortex results in the vertical alignment of the vortex column, which is a key step in the rapid intensification process. In the case of Sally, the spinup of the low-level vortex resulted from vorticity stretching, while the spinup of the midlevel vortex at 6 km resulted from vorticity tilting produced by the interaction of convective ascent with moderate vertical shear.
Significance Statement
The purpose of this study is to investigate the rapid intensification of Hurricane Sally as it was approaching the Florida Panhandle. We do that by analyzing an unprecedented dataset from the NOAA WP-3D and Gulfstream-IV aircraft, together with ground-based radar and satellite data. We find that both the dynamics (vorticity structure and evolution) and thermodynamics (instability index, saturation fraction, heating/mass flux profiles) need to be considered in diagnosing intensification processes. Further field projects with continuous high-altitude dropsondes and research are needed to see if these are applicable to other reformation events as well as genesis.
Abstract
Entrainment is a key process that can modulate the intensity of supercells, and a better understanding of its impact could help improve forecasts of thunderstorms and the precipitation they produce. In Part III of this series, the three distinct mechanisms of entrainment identified during the mature stage of idealized supercell thunderstorms in Part I (overturning “ribbons” of horizontal vorticity, “disorganized turbulent eddies,” and the “storm-relative airstream”) are examined as the absolute humidity of the environment is decreased. The existence of these mechanisms in a more realistic simulated storm environment is also established. Entrainment is calculated as fluxes of air across the storm core surface; passive fluid tracers help determine the resulting dilution of the storm updraft. Model microphysical rates are used to examine the direct impacts of entrainment on hydrometeors within the storm updraft as well as precipitation that falls to the ground. Results show that as mixed-layer humidity decreases, the “ribbons” and turbulent eddy mechanisms decrease in intensity, but their effects on precipitation production change little. With decreasing humidity in the 3–4 km AGL layer, the storm-relative airstream entrains less humid low-level air into the storm core, decreasing the vertical mass flux, and therefore the precipitation produced by the storm. When the humidity in the mid- to upper troposphere (4–20 km AGL) is decreased, precipitation is significantly reduced, but not due to the effects of the entrained air. Rather, enhanced evaporation and sublimation of falling precipitation decreases the overall precipitation efficiency of the storm.
Abstract
Entrainment is a key process that can modulate the intensity of supercells, and a better understanding of its impact could help improve forecasts of thunderstorms and the precipitation they produce. In Part III of this series, the three distinct mechanisms of entrainment identified during the mature stage of idealized supercell thunderstorms in Part I (overturning “ribbons” of horizontal vorticity, “disorganized turbulent eddies,” and the “storm-relative airstream”) are examined as the absolute humidity of the environment is decreased. The existence of these mechanisms in a more realistic simulated storm environment is also established. Entrainment is calculated as fluxes of air across the storm core surface; passive fluid tracers help determine the resulting dilution of the storm updraft. Model microphysical rates are used to examine the direct impacts of entrainment on hydrometeors within the storm updraft as well as precipitation that falls to the ground. Results show that as mixed-layer humidity decreases, the “ribbons” and turbulent eddy mechanisms decrease in intensity, but their effects on precipitation production change little. With decreasing humidity in the 3–4 km AGL layer, the storm-relative airstream entrains less humid low-level air into the storm core, decreasing the vertical mass flux, and therefore the precipitation produced by the storm. When the humidity in the mid- to upper troposphere (4–20 km AGL) is decreased, precipitation is significantly reduced, but not due to the effects of the entrained air. Rather, enhanced evaporation and sublimation of falling precipitation decreases the overall precipitation efficiency of the storm.