Search Results
You are looking at 11 - 20 of 40 items for
- Author or Editor: S. Richardson x
- Refine by Access: All Content x
Abstract
A dryline and misocyclones have been simulated in a cloud-resolving model by applying specified initial and time-dependent lateral boundary conditions obtained from analyses of the 22 May 2002 International H2O Project (IHOP_2002) dataset. The initial and lateral boundary conditions were obtained from a combination of the time–spaced Lagrangian analyses for temperature and moisture with horizontal velocities from multiple-Doppler wind syntheses. The simulated dryline, horizontal dry-convective rolls (HCRs) and open cells (OCCs), misocyclones, and cumulus clouds are similar to the corresponding observed features. The misocyclones move northward at nearly the mean boundary layer (BL) wind speed, rotate dryline gradients owing to their circulations, and move the local dryline eastward via their passage. Cumuli develop along a secondary dryline, along HCR and OCC segments between the primary and secondary drylines, along HCR and OCC segments that have moved over the dryline, and within the dryline updraft. After the initial shearing instability develops, misocyclogenesis proceeds from tilting and stretching of vorticity by the persistent secondary dryline circulation. The resulting misocyclone evolution is discussed.
Abstract
A dryline and misocyclones have been simulated in a cloud-resolving model by applying specified initial and time-dependent lateral boundary conditions obtained from analyses of the 22 May 2002 International H2O Project (IHOP_2002) dataset. The initial and lateral boundary conditions were obtained from a combination of the time–spaced Lagrangian analyses for temperature and moisture with horizontal velocities from multiple-Doppler wind syntheses. The simulated dryline, horizontal dry-convective rolls (HCRs) and open cells (OCCs), misocyclones, and cumulus clouds are similar to the corresponding observed features. The misocyclones move northward at nearly the mean boundary layer (BL) wind speed, rotate dryline gradients owing to their circulations, and move the local dryline eastward via their passage. Cumuli develop along a secondary dryline, along HCR and OCC segments between the primary and secondary drylines, along HCR and OCC segments that have moved over the dryline, and within the dryline updraft. After the initial shearing instability develops, misocyclogenesis proceeds from tilting and stretching of vorticity by the persistent secondary dryline circulation. The resulting misocyclone evolution is discussed.
Abstract
Precipitation forecasts from five global numerical weather prediction (NWP) models are verified against rain gauge observations using the new stable equitable error in probability space (SEEPS) score. It is based on a 3 × 3 contingency table and measures the ability of a forecast to discriminate between “dry,” “light precipitation,” and “heavy precipitation.” In SEEPS, the threshold defining the boundary between the light and heavy categories varies systematically with precipitation climate. Results obtained for SEEPS are compared to those of more well-known scores, and are broken down with regard to individual contributions from the contingency table. It is found that differences in skill between the models are consistent for different scores, but are small compared to seasonal and geographical variations, which themselves can be largely ascribed to the varying prevalence of deep convection. Differences between the tropics and extratropics are quite pronounced. SEEPS scores at forecast day 1 in the tropics are similar to those at day 6 in the extratropics. It is found that the model ranking is robust with respect to choices in the score computation. The issue of observation representativeness is addressed using a “quasi-perfect model” approach. Results suggest that just under one-half of the current forecast error at day 1 in the extratropics can be attributed to the fact that gridbox values are verified against point observations.
Abstract
Precipitation forecasts from five global numerical weather prediction (NWP) models are verified against rain gauge observations using the new stable equitable error in probability space (SEEPS) score. It is based on a 3 × 3 contingency table and measures the ability of a forecast to discriminate between “dry,” “light precipitation,” and “heavy precipitation.” In SEEPS, the threshold defining the boundary between the light and heavy categories varies systematically with precipitation climate. Results obtained for SEEPS are compared to those of more well-known scores, and are broken down with regard to individual contributions from the contingency table. It is found that differences in skill between the models are consistent for different scores, but are small compared to seasonal and geographical variations, which themselves can be largely ascribed to the varying prevalence of deep convection. Differences between the tropics and extratropics are quite pronounced. SEEPS scores at forecast day 1 in the tropics are similar to those at day 6 in the extratropics. It is found that the model ranking is robust with respect to choices in the score computation. The issue of observation representativeness is addressed using a “quasi-perfect model” approach. Results suggest that just under one-half of the current forecast error at day 1 in the extratropics can be attributed to the fact that gridbox values are verified against point observations.
Abstract
A supercell produced a nearly tornadic vortex during an intercept by the Second Verification of the Origins of Rotation in Tornadoes Experiment on 26 May 2010. Using observations from two mobile radars performing dual-Doppler scans, a five-probe mobile mesonet, and a proximity sounding, factors that prevented this vortex from strengthening into a significant tornado are examined. Mobile mesonet observations indicate that portions of the supercell outflow possessed excessive negative buoyancy, likely owing in part to low boundary layer relative humidity, as indicated by a high environmental lifted condensation level. Comparisons to a tornadic supercell suggest that the Prospect Valley storm had enough far-field circulation to produce a significant tornado, but was unable to converge this circulation to a sufficiently small radius. Trajectories suggest that the weak convergence might be due to the low-level mesocyclone ingesting parcels with considerable crosswise vorticity from the near-storm environment, which has been found to contribute to less steady and weaker low-level updrafts in supercell simulations. Yet another factor that likely contributed to the weak low-level circulation was the inability of parcels rich in streamwise vorticity from the forward-flank precipitation region to reach the low-level mesocyclone, likely owing to an unfavorable pressure gradient force field. In light of these results, we suggest that future research should continue focusing on the role of internal, storm-scale processes in tornadogenesis, especially in marginal environments.
Abstract
A supercell produced a nearly tornadic vortex during an intercept by the Second Verification of the Origins of Rotation in Tornadoes Experiment on 26 May 2010. Using observations from two mobile radars performing dual-Doppler scans, a five-probe mobile mesonet, and a proximity sounding, factors that prevented this vortex from strengthening into a significant tornado are examined. Mobile mesonet observations indicate that portions of the supercell outflow possessed excessive negative buoyancy, likely owing in part to low boundary layer relative humidity, as indicated by a high environmental lifted condensation level. Comparisons to a tornadic supercell suggest that the Prospect Valley storm had enough far-field circulation to produce a significant tornado, but was unable to converge this circulation to a sufficiently small radius. Trajectories suggest that the weak convergence might be due to the low-level mesocyclone ingesting parcels with considerable crosswise vorticity from the near-storm environment, which has been found to contribute to less steady and weaker low-level updrafts in supercell simulations. Yet another factor that likely contributed to the weak low-level circulation was the inability of parcels rich in streamwise vorticity from the forward-flank precipitation region to reach the low-level mesocyclone, likely owing to an unfavorable pressure gradient force field. In light of these results, we suggest that future research should continue focusing on the role of internal, storm-scale processes in tornadogenesis, especially in marginal environments.
Abstract
Spatial variability of precipitation is analyzed to characterize to what extent precipitation observed at a single location is representative of precipitation over a larger area. Characterization of precipitation representativeness is made in probabilistic terms using a parametric approach, namely, by fitting a censored shifted gamma distribution to observation measurements. Parameters are estimated and analyzed for independent precipitation datasets, among which one is based on high-density gauge measurements. The results of this analysis serve as a basis for accounting for representativeness error in an ensemble verification process. Uncertainty associated with the scale mismatch between forecast and observation is accounted for by applying a perturbed-ensemble approach before the computation of scores. Verification results reveal a large impact of representativeness error on precipitation forecast reliability and skill estimates. The parametric model and estimated coefficients presented in this study could be used directly for forecast postprocessing to partly compensate for the limitation of any modeling system in terms of precipitation subgrid-scale variability.
Abstract
Spatial variability of precipitation is analyzed to characterize to what extent precipitation observed at a single location is representative of precipitation over a larger area. Characterization of precipitation representativeness is made in probabilistic terms using a parametric approach, namely, by fitting a censored shifted gamma distribution to observation measurements. Parameters are estimated and analyzed for independent precipitation datasets, among which one is based on high-density gauge measurements. The results of this analysis serve as a basis for accounting for representativeness error in an ensemble verification process. Uncertainty associated with the scale mismatch between forecast and observation is accounted for by applying a perturbed-ensemble approach before the computation of scores. Verification results reveal a large impact of representativeness error on precipitation forecast reliability and skill estimates. The parametric model and estimated coefficients presented in this study could be used directly for forecast postprocessing to partly compensate for the limitation of any modeling system in terms of precipitation subgrid-scale variability.
Abstract
We investigate the run-to-run consistency (jumpiness) of ensemble forecasts of tropical cyclone tracks from three global centers: ECMWF, the Met Office, and NCEP. We use a divergence function to quantify the change in cross-track position between consecutive ensemble forecasts initialized at 12-h intervals. Results for the 2019–21 North Atlantic hurricane season show that the jumpiness varied substantially between cases and centers, with no common cause across the different ensemble systems. Recent upgrades to the Met Office and NCEP ensembles reduced their overall jumpiness to match that of the ECMWF ensemble. The average divergence over the set of cases provides an objective measure of the expected change in cross-track position from one forecast to the next. For example, a user should expect on average that the ensemble mean position will change by around 80–90 km in the cross-track direction between a forecast for 120 h ahead and the updated forecast made 12 h later for the same valid time. This quantitative information can support users’ decision-making, for example, in deciding whether to act now or wait for the next forecast. We did not find any link between jumpiness and skill, indicating that users should not rely on the consistency between successive forecasts as a measure of confidence. Instead, we suggest that users should use ensemble spread and probabilistic information to assess forecast uncertainty, and consider multimodel combinations to reduce the effects of jumpiness.
Significance Statement
Forecasting the tracks of tropical cyclones is essential to mitigate their impacts on society. Numerical weather prediction models provide valuable guidance, but occasionally there is a large jump in the predicted track from one run to the next. This jumpiness complicates the creation and communication of consistent forecast advisories and early warnings. In this work we aim to better understand forecast jumpiness and we provide practical information to forecasters to help them better use the model guidance. We show that the jumpiest cases are different for different modeling centers, that recent model upgrades have reduced forecast jumpiness, and that there is not a strong link between jumpiness and forecast skill.
Abstract
We investigate the run-to-run consistency (jumpiness) of ensemble forecasts of tropical cyclone tracks from three global centers: ECMWF, the Met Office, and NCEP. We use a divergence function to quantify the change in cross-track position between consecutive ensemble forecasts initialized at 12-h intervals. Results for the 2019–21 North Atlantic hurricane season show that the jumpiness varied substantially between cases and centers, with no common cause across the different ensemble systems. Recent upgrades to the Met Office and NCEP ensembles reduced their overall jumpiness to match that of the ECMWF ensemble. The average divergence over the set of cases provides an objective measure of the expected change in cross-track position from one forecast to the next. For example, a user should expect on average that the ensemble mean position will change by around 80–90 km in the cross-track direction between a forecast for 120 h ahead and the updated forecast made 12 h later for the same valid time. This quantitative information can support users’ decision-making, for example, in deciding whether to act now or wait for the next forecast. We did not find any link between jumpiness and skill, indicating that users should not rely on the consistency between successive forecasts as a measure of confidence. Instead, we suggest that users should use ensemble spread and probabilistic information to assess forecast uncertainty, and consider multimodel combinations to reduce the effects of jumpiness.
Significance Statement
Forecasting the tracks of tropical cyclones is essential to mitigate their impacts on society. Numerical weather prediction models provide valuable guidance, but occasionally there is a large jump in the predicted track from one run to the next. This jumpiness complicates the creation and communication of consistent forecast advisories and early warnings. In this work we aim to better understand forecast jumpiness and we provide practical information to forecasters to help them better use the model guidance. We show that the jumpiest cases are different for different modeling centers, that recent model upgrades have reduced forecast jumpiness, and that there is not a strong link between jumpiness and forecast skill.
Abstract
Several studies have documented the sensitivity of convective storm simulations to the microphysics parameterization, but there is less research documenting how these sensitivities change with environmental conditions. In this study, the influence of the lifting condensation level (LCL) on the sensitivity of simulated ordinary convective storm cold pools to the microphysics parameterization is examined. To do this, seven perturbed-microphysics ensembles with nine members each are used, where each ensemble uses a different base state with a surface-based LCL between 500 and 2000 m. A comparison of ensemble standard deviations of cold-pool properties shows a clear trend of increasing sensitivity to the microphysics as the LCL is raised. In physical terms, this trend is the result of lower relative humidities in high-LCL environments that increase low-level rain evaporational cooling rates, which magnifies differences in evaporation already present among the members of a given ensemble owing to the microphysics variations. Omitting supersaturation from the calculation of rain evaporation so that only the raindrop size distribution influences evaporation leads to more evaporation in the low-LCL simulations (owing to more drops), as well as a slightly larger spread in evaporational cooling amounts between members in the low-LCL ensembles. Cold pools in the low-LCL environments are also found to develop earlier and are initially more sensitive to raindrop breakup owing to a larger warm-cloud depth. Altogether, these results suggest that convective storms may be more predictable in low-LCL environments, and forecasts of convection in high-LCL environments may benefit the most from microphysics perturbations within an ensemble forecasting system.
Significance Statement
Computer simulations of thunderstorms can have grid spacings ranging from tens to thousands of meters. Because individual precipitation particles form on scales smaller than these grid spacings, the bulk effects of precipitation processes in models must be approximated. Past studies have found that models are sensitive to these approximations. In this study, we test whether the sensitivity to these approximations changes with the relative humidity in the lowest 1–2 km of the atmosphere. We found that increasing the relative humidity decreases the sensitivity of simulations to the precipitation process approximations. These results can inform meteorologists about the uncertainties surrounding computer-generated thunderstorm forecasts and suggest environmental conditions where using several computer models with different precipitation process approximations may be beneficial.
Abstract
Several studies have documented the sensitivity of convective storm simulations to the microphysics parameterization, but there is less research documenting how these sensitivities change with environmental conditions. In this study, the influence of the lifting condensation level (LCL) on the sensitivity of simulated ordinary convective storm cold pools to the microphysics parameterization is examined. To do this, seven perturbed-microphysics ensembles with nine members each are used, where each ensemble uses a different base state with a surface-based LCL between 500 and 2000 m. A comparison of ensemble standard deviations of cold-pool properties shows a clear trend of increasing sensitivity to the microphysics as the LCL is raised. In physical terms, this trend is the result of lower relative humidities in high-LCL environments that increase low-level rain evaporational cooling rates, which magnifies differences in evaporation already present among the members of a given ensemble owing to the microphysics variations. Omitting supersaturation from the calculation of rain evaporation so that only the raindrop size distribution influences evaporation leads to more evaporation in the low-LCL simulations (owing to more drops), as well as a slightly larger spread in evaporational cooling amounts between members in the low-LCL ensembles. Cold pools in the low-LCL environments are also found to develop earlier and are initially more sensitive to raindrop breakup owing to a larger warm-cloud depth. Altogether, these results suggest that convective storms may be more predictable in low-LCL environments, and forecasts of convection in high-LCL environments may benefit the most from microphysics perturbations within an ensemble forecasting system.
Significance Statement
Computer simulations of thunderstorms can have grid spacings ranging from tens to thousands of meters. Because individual precipitation particles form on scales smaller than these grid spacings, the bulk effects of precipitation processes in models must be approximated. Past studies have found that models are sensitive to these approximations. In this study, we test whether the sensitivity to these approximations changes with the relative humidity in the lowest 1–2 km of the atmosphere. We found that increasing the relative humidity decreases the sensitivity of simulations to the precipitation process approximations. These results can inform meteorologists about the uncertainties surrounding computer-generated thunderstorm forecasts and suggest environmental conditions where using several computer models with different precipitation process approximations may be beneficial.
Interaction Between the Atmosphere and the Oceans
Report of the Joint Panel on Air-Sea Interaction
Abstract
Large-scale cloud features referred to as cloudbands are known to be related to widespread and heavy rain via the transport of tropical heat and moisture to higher latitudes. The Australian northwest cloudband is such a feature that has been identified in simple searches of satellite imagery but with limited investigation of its atmospheric dynamical support. An accurate, long-term climatology of northwest cloudbands is key to robustly assessing these events. A dynamically based search algorithm has been developed that is guided by the presence and orientation of the subtropical jet stream. This jet stream is the large-scale atmospheric feature that determines the development and alignment of a cloudband. Using a new 40-yr dataset of cloudband events compiled by this search algorithm, composite atmospheric and ocean surface conditions over the period 1979–2018 have been assessed. Composite cloudband upper-level flow revealed a tilted low pressure trough embedded in a Rossby wave train. Composites of vertically integrated water vapor transport centered around the jet maximum during northwest cloudband events reveal a distinct atmospheric river supplying tropical moisture for cloudband rainfall. Parcel backtracking indicated multiple regions of moisture support for cloudbands. A thermal wind anomaly orientated with respect to an enhanced sea surface temperature gradient over the Indian Ocean was also a key composite cloudband feature. A total of 300 years of a freely coupled control simulation of the ACCESS-D system was assessed for its ability to simulate northwest cloudbands. Composite analysis of model cloudbands compared reasonably well to reanalysis despite some differences in seasonality and frequency of occurrence.
Abstract
Large-scale cloud features referred to as cloudbands are known to be related to widespread and heavy rain via the transport of tropical heat and moisture to higher latitudes. The Australian northwest cloudband is such a feature that has been identified in simple searches of satellite imagery but with limited investigation of its atmospheric dynamical support. An accurate, long-term climatology of northwest cloudbands is key to robustly assessing these events. A dynamically based search algorithm has been developed that is guided by the presence and orientation of the subtropical jet stream. This jet stream is the large-scale atmospheric feature that determines the development and alignment of a cloudband. Using a new 40-yr dataset of cloudband events compiled by this search algorithm, composite atmospheric and ocean surface conditions over the period 1979–2018 have been assessed. Composite cloudband upper-level flow revealed a tilted low pressure trough embedded in a Rossby wave train. Composites of vertically integrated water vapor transport centered around the jet maximum during northwest cloudband events reveal a distinct atmospheric river supplying tropical moisture for cloudband rainfall. Parcel backtracking indicated multiple regions of moisture support for cloudbands. A thermal wind anomaly orientated with respect to an enhanced sea surface temperature gradient over the Indian Ocean was also a key composite cloudband feature. A total of 300 years of a freely coupled control simulation of the ACCESS-D system was assessed for its ability to simulate northwest cloudbands. Composite analysis of model cloudbands compared reasonably well to reanalysis despite some differences in seasonality and frequency of occurrence.
Abstract
Tropical cyclones are some of the most devastating natural hazards and the “three beasts”—Harvey, Irma, and Maria—during the Atlantic hurricane season 2017 are recent examples. The European Centre for Medium-Range Weather Forecasts (ECMWF) is working on fulfilling its 2016–25 strategy in which early warnings for extreme events will be made possible by a high-resolution Earth system ensemble forecasting system. Several verification reports acknowledge deterministic and probabilistic tropical cyclone tracks from ECMWF as world leading. However, producing reliable intensity forecasts is still a difficult task for the ECMWF global forecasting model, especially regarding maximum wind speed. This article will put the ECMWF strategy into a tropical cyclone perspective and highlight some key research activities, using Harvey, Irma, and Maria as examples. We describe the observation usage around tropical cyclones in data assimilation and give examples of their impact. From a model perspective, we show the impact of running at 5-km resolution and also the impact of applying ocean coupling. Finally, we discuss the future challenges to tackle the errors in intensity forecasts for tropical cyclones.
Abstract
Tropical cyclones are some of the most devastating natural hazards and the “three beasts”—Harvey, Irma, and Maria—during the Atlantic hurricane season 2017 are recent examples. The European Centre for Medium-Range Weather Forecasts (ECMWF) is working on fulfilling its 2016–25 strategy in which early warnings for extreme events will be made possible by a high-resolution Earth system ensemble forecasting system. Several verification reports acknowledge deterministic and probabilistic tropical cyclone tracks from ECMWF as world leading. However, producing reliable intensity forecasts is still a difficult task for the ECMWF global forecasting model, especially regarding maximum wind speed. This article will put the ECMWF strategy into a tropical cyclone perspective and highlight some key research activities, using Harvey, Irma, and Maria as examples. We describe the observation usage around tropical cyclones in data assimilation and give examples of their impact. From a model perspective, we show the impact of running at 5-km resolution and also the impact of applying ocean coupling. Finally, we discuss the future challenges to tackle the errors in intensity forecasts for tropical cyclones.
Abstract
The ability to find and recognize patterns in high-dimensional geophysical data is fundamental to climate science and critical for meaningful interpretation of weather and climate processes. Archetypal analysis (AA) is one technique that has recently gained traction in the geophysical science community for its ability to find patterns based on extreme conditions. While traditional empirical orthogonal function (EOF) analysis can reveal patterns based on data covariance, AA seeks patterns from the points located at the edges of the data distribution. The utility of any objective pattern method depends on the properties of the data to which it is applied and the choices made in implementing the method. Given the relative novelty of the application of AA in geophysics it is important to develop experience in applying the method. We provide an assessment of the method, implementation, sensitivity, and interpretation of AA with respect to geophysical data. As an example for demonstration, we apply AA to a 39-yr sea surface temperature (SST) reanalysis dataset. We show that the decisions made to implement AA can significantly affect the interpretation of results, but also, in the case of SST, that the analysis is exceptionally robust under both spatial and temporal coarse graining.
Significance Statement
Archetypal analysis (AA), when applied to geophysical fields, is a technique designed to find typical configurations or modes in underlying data. This technique is relatively new to the geophysical science community and has been shown to be beneficial to the interpretation of extreme modes of the climate system. The identification of extreme modes of variability and their expression in day-to-day weather or state of the climate at longer time scales may help in elucidating the interplay between major teleconnection drivers and their evolution in a changing climate. The purpose of this work is to bring together a comprehensive report of the AA methodology using an SST reanalysis for demonstration. It is shown that the AA results are significantly affected by each implementation decision, but also can be resilient to spatiotemporal averaging. Any application of AA should provide a clear documentation of the choices made in applying the method.
Abstract
The ability to find and recognize patterns in high-dimensional geophysical data is fundamental to climate science and critical for meaningful interpretation of weather and climate processes. Archetypal analysis (AA) is one technique that has recently gained traction in the geophysical science community for its ability to find patterns based on extreme conditions. While traditional empirical orthogonal function (EOF) analysis can reveal patterns based on data covariance, AA seeks patterns from the points located at the edges of the data distribution. The utility of any objective pattern method depends on the properties of the data to which it is applied and the choices made in implementing the method. Given the relative novelty of the application of AA in geophysics it is important to develop experience in applying the method. We provide an assessment of the method, implementation, sensitivity, and interpretation of AA with respect to geophysical data. As an example for demonstration, we apply AA to a 39-yr sea surface temperature (SST) reanalysis dataset. We show that the decisions made to implement AA can significantly affect the interpretation of results, but also, in the case of SST, that the analysis is exceptionally robust under both spatial and temporal coarse graining.
Significance Statement
Archetypal analysis (AA), when applied to geophysical fields, is a technique designed to find typical configurations or modes in underlying data. This technique is relatively new to the geophysical science community and has been shown to be beneficial to the interpretation of extreme modes of the climate system. The identification of extreme modes of variability and their expression in day-to-day weather or state of the climate at longer time scales may help in elucidating the interplay between major teleconnection drivers and their evolution in a changing climate. The purpose of this work is to bring together a comprehensive report of the AA methodology using an SST reanalysis for demonstration. It is shown that the AA results are significantly affected by each implementation decision, but also can be resilient to spatiotemporal averaging. Any application of AA should provide a clear documentation of the choices made in applying the method.
