Browse
Abstract
We investigate the practical predictability limits of deep convection in a state-of-the-art, high-resolution, limited-area ensemble prediction system. A combination of sophisticated predictability measures, namely, believable and decorrelation scale, are applied to determine the predictable scales of short-term forecasts in a hierarchy of model configurations. First, we consider an idealized perfect model setup that includes both small-scale and synoptic-scale perturbations. We find increased predictability in the presence of orography and a strongly beneficial impact of radar data assimilation, which extends the forecast horizon by up to 6 h. Second, we examine realistic COSMO-KENDA simulations, including assimilation of radar and conventional data and a representation of model errors, for a convectively active two-week summer period over Germany. The results confirm increased predictability in orographic regions. We find that both latent heat nudging and ensemble Kalman filter assimilation of radar data lead to increased forecast skill, but the impact is smaller than in the idealized experiments. This highlights the need to assimilate spatially and temporally dense data, but also indicates room for further improvement. Finally, the examination of operational COSMO-DE-EPS ensemble forecasts for three summer periods confirms the beneficial impact of orography in a statistical sense and also reveals increased predictability in weather regimes controlled by synoptic forcing, as defined by the convective adjustment time scale.
Abstract
We investigate the practical predictability limits of deep convection in a state-of-the-art, high-resolution, limited-area ensemble prediction system. A combination of sophisticated predictability measures, namely, believable and decorrelation scale, are applied to determine the predictable scales of short-term forecasts in a hierarchy of model configurations. First, we consider an idealized perfect model setup that includes both small-scale and synoptic-scale perturbations. We find increased predictability in the presence of orography and a strongly beneficial impact of radar data assimilation, which extends the forecast horizon by up to 6 h. Second, we examine realistic COSMO-KENDA simulations, including assimilation of radar and conventional data and a representation of model errors, for a convectively active two-week summer period over Germany. The results confirm increased predictability in orographic regions. We find that both latent heat nudging and ensemble Kalman filter assimilation of radar data lead to increased forecast skill, but the impact is smaller than in the idealized experiments. This highlights the need to assimilate spatially and temporally dense data, but also indicates room for further improvement. Finally, the examination of operational COSMO-DE-EPS ensemble forecasts for three summer periods confirms the beneficial impact of orography in a statistical sense and also reveals increased predictability in weather regimes controlled by synoptic forcing, as defined by the convective adjustment time scale.
Abstract
Numerical weather, climate, or Earth system models involve the coupling of components. At a broad level, these components can be classified as the resolved fluid dynamics, unresolved fluid dynamical aspects (i.e., those represented by physical parameterizations such as subgrid-scale mixing), and nonfluid dynamical aspects such as radiation and microphysical processes. Typically, each component is developed, at least initially, independently. Once development is mature, the components are coupled to deliver a model of the required complexity. The implementation of the coupling can have a significant impact on the model. As the error associated with each component decreases, the errors introduced by the coupling will eventually dominate. Hence, any improvement in one of the components is unlikely to improve the performance of the overall system. The challenges associated with combining the components to create a coherent model are here termed physics–dynamics coupling. The issue goes beyond the coupling between the parameterizations and the resolved fluid dynamics. This paper highlights recent progress and some of the current challenges. It focuses on three objectives: to illustrate the phenomenology of the coupling problem with references to examples in the literature, to show how the problem can be analyzed, and to create awareness of the issue across the disciplines and specializations. The topics addressed are different ways of advancing full models in time, approaches to understanding the role of the coupling and evaluation of approaches, coupling ocean and atmosphere models, thermodynamic compatibility between model components, and emerging issues such as those that arise as model resolutions increase and/or models use variable resolutions.
Abstract
Numerical weather, climate, or Earth system models involve the coupling of components. At a broad level, these components can be classified as the resolved fluid dynamics, unresolved fluid dynamical aspects (i.e., those represented by physical parameterizations such as subgrid-scale mixing), and nonfluid dynamical aspects such as radiation and microphysical processes. Typically, each component is developed, at least initially, independently. Once development is mature, the components are coupled to deliver a model of the required complexity. The implementation of the coupling can have a significant impact on the model. As the error associated with each component decreases, the errors introduced by the coupling will eventually dominate. Hence, any improvement in one of the components is unlikely to improve the performance of the overall system. The challenges associated with combining the components to create a coherent model are here termed physics–dynamics coupling. The issue goes beyond the coupling between the parameterizations and the resolved fluid dynamics. This paper highlights recent progress and some of the current challenges. It focuses on three objectives: to illustrate the phenomenology of the coupling problem with references to examples in the literature, to show how the problem can be analyzed, and to create awareness of the issue across the disciplines and specializations. The topics addressed are different ways of advancing full models in time, approaches to understanding the role of the coupling and evaluation of approaches, coupling ocean and atmosphere models, thermodynamic compatibility between model components, and emerging issues such as those that arise as model resolutions increase and/or models use variable resolutions.
Abstract
Deep moist convection is an inherently multiscale phenomenon with organization processes coupling convective elements to larger-scale structures. A realistic representation of the tropical dynamics demands a simulation framework that is capable of representing physical processes across a wide range of scales. Therefore, storm-resolving numerical simulations at 2.4 km have been performed covering the tropical Atlantic and neighboring parts for 2 months. The simulated cloud fields are combined with infrared geostationary satellite observations, and their realism is assessed with the help of object-based evaluation methods. It is shown that the simulations are able to develop a well-defined intertropical convergence zone. However, marine convective activity measured by the cold cloud coverage is considerably underestimated, especially for the winter season and the western Atlantic. The spatial coupling across the resolved scales leads to simulated cloud number size distributions that follow power laws similar to the observations, with slopes steeper in winter than summer and slopes steeper over ocean than over land. The simulated slopes are, however, too steep, indicating too many small and too few large tropical cloud cells. It is also discussed that the number of larger cells is less influenced by multiday variability of environmental conditions. Despite the identified deficits, the analyzed simulations highlight the great potential of this modeling framework for process-based studies of tropical deep convection.
Abstract
Deep moist convection is an inherently multiscale phenomenon with organization processes coupling convective elements to larger-scale structures. A realistic representation of the tropical dynamics demands a simulation framework that is capable of representing physical processes across a wide range of scales. Therefore, storm-resolving numerical simulations at 2.4 km have been performed covering the tropical Atlantic and neighboring parts for 2 months. The simulated cloud fields are combined with infrared geostationary satellite observations, and their realism is assessed with the help of object-based evaluation methods. It is shown that the simulations are able to develop a well-defined intertropical convergence zone. However, marine convective activity measured by the cold cloud coverage is considerably underestimated, especially for the winter season and the western Atlantic. The spatial coupling across the resolved scales leads to simulated cloud number size distributions that follow power laws similar to the observations, with slopes steeper in winter than summer and slopes steeper over ocean than over land. The simulated slopes are, however, too steep, indicating too many small and too few large tropical cloud cells. It is also discussed that the number of larger cells is less influenced by multiday variability of environmental conditions. Despite the identified deficits, the analyzed simulations highlight the great potential of this modeling framework for process-based studies of tropical deep convection.
Abstract
Aircraft observations of wind and temperature collected by airport surveillance radars [Mode-S Enhanced Surveillance (Mode-S EHS)] were assimilated in the Consortium for Small-Scale Modeling Kilometre-scale Ensemble Data Assimilation (COSMO-KENDA), which couples an ensemble Kalman filter to a 40-member ensemble of the convection permitting COSMO-DE model. The number of observing aircrafts in Mode-S EHS was about 15 times larger than in the AMDAR system. In the comparison of both aircraft observation systems, a similar observation error standard deviation was diagnosed for wind. For temperature, a larger error was diagnosed for Mode-S EHS. With the high density of Mode-S EHS observations, a reduction of temperature and wind error in forecasts of 1 and 3 hours was found mainly in the flight level and less near the surface. The amount of Mode-S EHS data was reduced by random thinning to test the effect of a varying observation density. With the current data assimilation setup, a saturation of the forecast error reduction was apparent when more than 50% of the Mode-S EHS data were assimilated. Forecast kinetic energy spectra indicated that the reduction in error is related to analysis updates on all scales resolved by COSMO-DE.
Abstract
Aircraft observations of wind and temperature collected by airport surveillance radars [Mode-S Enhanced Surveillance (Mode-S EHS)] were assimilated in the Consortium for Small-Scale Modeling Kilometre-scale Ensemble Data Assimilation (COSMO-KENDA), which couples an ensemble Kalman filter to a 40-member ensemble of the convection permitting COSMO-DE model. The number of observing aircrafts in Mode-S EHS was about 15 times larger than in the AMDAR system. In the comparison of both aircraft observation systems, a similar observation error standard deviation was diagnosed for wind. For temperature, a larger error was diagnosed for Mode-S EHS. With the high density of Mode-S EHS observations, a reduction of temperature and wind error in forecasts of 1 and 3 hours was found mainly in the flight level and less near the surface. The amount of Mode-S EHS data was reduced by random thinning to test the effect of a varying observation density. With the current data assimilation setup, a saturation of the forecast error reduction was apparent when more than 50% of the Mode-S EHS data were assimilated. Forecast kinetic energy spectra indicated that the reduction in error is related to analysis updates on all scales resolved by COSMO-DE.
Abstract
A kilometer-scale ensemble data assimilation system (KENDA) based on a local ensemble transform Kalman filter (LETKF) has been developed for the Consortium for Small-Scale Modeling (COSMO) limited-area model. The data assimilation system provides an analysis ensemble that can be used to initialize ensemble forecasts at a horizontal grid resolution of 2.8 km. Convective-scale ensemble forecasts over Germany using ensemble initial conditions derived by the KENDA system are evaluated and compared to operational forecasts with downscaled initial conditions for a short summer period during June 2012.
The choice of the inflation method applied in the LETKF significantly affects the ensemble analysis and forecast. Using a multiplicative background covariance inflation does not produce enough spread in the analysis ensemble leading to a degradation of the ensemble forecasts. Inflating the analysis ensemble instead by either multiplicative analysis covariance inflation or relaxation inflation methods enhances the analysis spread and is able to provide initial conditions that produce more consistent ensemble forecasts. The forecast quality for short forecast lead times up to 3 h is improved, and 21-h forecasts also benefit from the increased spread.
Doubling the ensemble size has not only a clear positive impact on the analysis but also on the short-term ensemble forecasts, while a simple representation of model error perturbing parameters of the model physics has only a small impact. Precipitation and surface wind speed ensemble forecasts using the high-resolution KENDA-derived initial conditions are competitive compared to the operationally used downscaled initial conditions.
Abstract
A kilometer-scale ensemble data assimilation system (KENDA) based on a local ensemble transform Kalman filter (LETKF) has been developed for the Consortium for Small-Scale Modeling (COSMO) limited-area model. The data assimilation system provides an analysis ensemble that can be used to initialize ensemble forecasts at a horizontal grid resolution of 2.8 km. Convective-scale ensemble forecasts over Germany using ensemble initial conditions derived by the KENDA system are evaluated and compared to operational forecasts with downscaled initial conditions for a short summer period during June 2012.
The choice of the inflation method applied in the LETKF significantly affects the ensemble analysis and forecast. Using a multiplicative background covariance inflation does not produce enough spread in the analysis ensemble leading to a degradation of the ensemble forecasts. Inflating the analysis ensemble instead by either multiplicative analysis covariance inflation or relaxation inflation methods enhances the analysis spread and is able to provide initial conditions that produce more consistent ensemble forecasts. The forecast quality for short forecast lead times up to 3 h is improved, and 21-h forecasts also benefit from the increased spread.
Doubling the ensemble size has not only a clear positive impact on the analysis but also on the short-term ensemble forecasts, while a simple representation of model error perturbing parameters of the model physics has only a small impact. Precipitation and surface wind speed ensemble forecasts using the high-resolution KENDA-derived initial conditions are competitive compared to the operationally used downscaled initial conditions.
Abstract
An idealized convective test bed for the local ensemble transform Kalman filter (LETKF) is set up to perform storm-scale data assimilation of simulated Doppler radar observations. Convective systems with lifetimes exceeding 6 h are triggered in a doubly periodic domain. Perfect-model experiments are used to investigate the limited predictability in precipitation forecasts by comparing analysis schemes that resolve different length scales. Starting from a high-resolution reference scheme with 8-km covariance localization and observations with 2-km resolution on a 5-min cycle, an experimental hierarchy is set up by successively choosing a larger covariance localization radius of 32 km, observations that are horizontally averaged by a factor of 4, a coarser resolution in the calculation of the analysis weights, and a cycling interval of 20 min. After 3 h of assimilation, the high-resolution analysis scheme is clearly superior to the configurations with coarser scales in terms of RMS error and field-oriented measures. The difference is associated with the observation resolution and a larger localization radius required for filter convergence with coarse observations. The high-resolution analysis leads to better forecasts for the first hour, but after 3 hours, the forecast quality of the schemes is indistinguishable. The more rapid error growth in forecasts from the high-resolution analysis appears to be associated with a limited predictability of the small scales, but also with gravity wave noise and spurious convective cells. The latter suggests that the field is in some sense less balanced, or less consistent with the model dynamics, than in the coarser-resolution analysis.
Abstract
An idealized convective test bed for the local ensemble transform Kalman filter (LETKF) is set up to perform storm-scale data assimilation of simulated Doppler radar observations. Convective systems with lifetimes exceeding 6 h are triggered in a doubly periodic domain. Perfect-model experiments are used to investigate the limited predictability in precipitation forecasts by comparing analysis schemes that resolve different length scales. Starting from a high-resolution reference scheme with 8-km covariance localization and observations with 2-km resolution on a 5-min cycle, an experimental hierarchy is set up by successively choosing a larger covariance localization radius of 32 km, observations that are horizontally averaged by a factor of 4, a coarser resolution in the calculation of the analysis weights, and a cycling interval of 20 min. After 3 h of assimilation, the high-resolution analysis scheme is clearly superior to the configurations with coarser scales in terms of RMS error and field-oriented measures. The difference is associated with the observation resolution and a larger localization radius required for filter convergence with coarse observations. The high-resolution analysis leads to better forecasts for the first hour, but after 3 hours, the forecast quality of the schemes is indistinguishable. The more rapid error growth in forecasts from the high-resolution analysis appears to be associated with a limited predictability of the small scales, but also with gravity wave noise and spurious convective cells. The latter suggests that the field is in some sense less balanced, or less consistent with the model dynamics, than in the coarser-resolution analysis.
Abstract
This study analyzes the effects of land surface heterogeneities at various horizontal scales on the transition from shallow to deep convection and on the cloud size distribution. An idealized case of midlatitude summertime convection is simulated by means of large-eddy simulations coupled to an interactive land surface. The transition is accelerated over heterogeneous surfaces. The simulation with an intermediate patch size of 12.8 km exhibits the fastest transition with a transition time two-thirds that over a homogeneous surface. A similar timing is observed for the precipitation onset whereas the total accumulated rainfall tends to increase with patch size. The cloud size distribution can be approximated by a power law with a scale break. The exponent of the power law is independent of the heterogeneity scale, implying a similar cloud cover between the simulations. In contrast, the scale break varies with patch size. The size of the largest clouds does not scale with the boundary layer height, although their maximum size scales with the patch size. Finally, the idea that larger clouds grow faster, known from homogeneous surface conditions, is not fully valid over heterogeneous surfaces. These various aspects can be understood from the complex interplay between the characteristics of the triggered mesoscale circulations and a cloud development acting in response to the diurnal cycle in surface heating. The results also call for adequate representation of such effects in convective parameterizations.
Abstract
This study analyzes the effects of land surface heterogeneities at various horizontal scales on the transition from shallow to deep convection and on the cloud size distribution. An idealized case of midlatitude summertime convection is simulated by means of large-eddy simulations coupled to an interactive land surface. The transition is accelerated over heterogeneous surfaces. The simulation with an intermediate patch size of 12.8 km exhibits the fastest transition with a transition time two-thirds that over a homogeneous surface. A similar timing is observed for the precipitation onset whereas the total accumulated rainfall tends to increase with patch size. The cloud size distribution can be approximated by a power law with a scale break. The exponent of the power law is independent of the heterogeneity scale, implying a similar cloud cover between the simulations. In contrast, the scale break varies with patch size. The size of the largest clouds does not scale with the boundary layer height, although their maximum size scales with the patch size. Finally, the idea that larger clouds grow faster, known from homogeneous surface conditions, is not fully valid over heterogeneous surfaces. These various aspects can be understood from the complex interplay between the characteristics of the triggered mesoscale circulations and a cloud development acting in response to the diurnal cycle in surface heating. The results also call for adequate representation of such effects in convective parameterizations.
Abstract
This paper considers the incorporation of constraints to enforce physically based conservation laws in the ensemble Kalman filter. In particular, constraints are used to ensure that the ensemble members and the ensemble mean conserve mass and remain nonnegative through measurement updates. In certain situations filtering algorithms such as the ensemble Kalman filter (EnKF) and ensemble transform Kalman filter (ETKF) yield updated ensembles that conserve mass but are negative, even though the actual states must be nonnegative. In such situations if negative values are set to zero, or a log transform is introduced, the total mass will not be conserved. In this study, mass and positivity are both preserved by formulating the filter update as a set of quadratic programming problems that incorporate nonnegativity constraints. Simple numerical experiments indicate that this approach can have a significant positive impact on the posterior ensemble distribution, giving results that are more physically plausible both for individual ensemble members and for the ensemble mean. In two examples, an update that includes a nonnegativity constraint is able to properly describe the transport of a sharp feature (e.g., a triangle or cone). A number of implementation questions still need to be addressed, particularly the need to develop a computationally efficient quadratic programming update for large ensemble.
Abstract
This paper considers the incorporation of constraints to enforce physically based conservation laws in the ensemble Kalman filter. In particular, constraints are used to ensure that the ensemble members and the ensemble mean conserve mass and remain nonnegative through measurement updates. In certain situations filtering algorithms such as the ensemble Kalman filter (EnKF) and ensemble transform Kalman filter (ETKF) yield updated ensembles that conserve mass but are negative, even though the actual states must be nonnegative. In such situations if negative values are set to zero, or a log transform is introduced, the total mass will not be conserved. In this study, mass and positivity are both preserved by formulating the filter update as a set of quadratic programming problems that incorporate nonnegativity constraints. Simple numerical experiments indicate that this approach can have a significant positive impact on the posterior ensemble distribution, giving results that are more physically plausible both for individual ensemble members and for the ensemble mean. In two examples, an update that includes a nonnegativity constraint is able to properly describe the transport of a sharp feature (e.g., a triangle or cone). A number of implementation questions still need to be addressed, particularly the need to develop a computationally efficient quadratic programming update for large ensemble.
