Options and extensions for the stochastic shallow convection scheme in ICON

Author(s):  
Maike Ahlgrimm ◽  
Daniel Klocke ◽  
Alberto de Lozar ◽  
Ekaterina Machulskaya ◽  
Mirjana Sakradzija ◽  
...  

<p>The Icosahedral Model (ICON) of the German Weather Service (Deutscher Wetterdienst, DWD) is used for numerical weather prediction at global and regional scales. In the limited area mode, resolution is typically on the order of a few kilometers horizontal grid spacing. Deep convective transport is partially resolved at these scales, but shallow convection remains poorly represented without a parameterization.</p><p>A stochastic shallow convection scheme was developed in collaboration with the Max Planck Institute for Meteorology, and is now being implemented in ICON with a view towards operational use. The scheme is scale-adaptive and renders resolution-dependent tuning of the convection parameterization unnecessary. Mass flux limiters essential for the stable operation of the unaltered convection scheme can be removed when the stochastic perturbations are introduced.</p><p>Alongside the original, explicit stochastic scheme an approximation using stochastic differential equations (SDE) has been developed. The advantage of the SDE version is a lower computational and memory cost, and the ability to save and restart the model‘s stochastic cloud state easily.</p><p>Equivalence of the two versions can be demonstrated by running one version interactively, the other passively (“piggy-backing”). While the SDE approximation is computationally more efficient, the explicit version of the scheme can be easily extended to keep track of additional properties of the shallow cloud ensemble. For example, the convective updraft core fraction can be calculated for use in the diagnostic subgrid cloud scheme. Or knowledge of individual clouds’ depth can be used to derive a more realistic lateral detrainment profile than is currently in use.</p><p>We demonstrate the performance of the scheme and illustrate options and applications in single column mode, case studies and month-long hindcasts.</p>

2018 ◽  
Author(s):  
Clemens Wastl ◽  
Yong Wang ◽  
Aitor Atencia ◽  
Christoph Wittman

Abstract. A modification of the widely used SPPT (Stochastically Perturbed Parametrisation Tendencies) scheme is proposed and tested in a Convection-permitting – Limited Area Ensemble Forecasting system (C-LAEF) developed at ZAMG (Zentralanstalt für Meteorologie und Geodynamik). The tendencies from four physical parametrisation schemes are perturbed: radiation, shallow convection, turbulence and microphysics. Whereas in SPPT the total model tendencies are perturbed, in the present approach (pSPPT hereinafter) the partial tendencies of the physics parametrisation schemes are sequentially perturbed. Thus, in pSPPT an interaction between the uncertainties of the different physics parametrisation schemes is sustained and a more physically consistent relationship between the processes is kept. Two configurations of pSPPT are evaluated over two months (one of summer and another of winter). Both schemes increase the stability of the model and lead to statistically significant improvements in the probabilistic performance compared to the original SPPT. An evaluation of selected test cases shows that the positive effect of stochastic physics is much more pronounced on days with high convective activity. Small discrepancies in the humidity analysis can be dedicated to the use of a very simple supersaturation adjustment. This and other adjustments are discussed to provide some suggestions for future investigations.


2005 ◽  
Vol 133 (8) ◽  
pp. 2422-2435 ◽  
Author(s):  
Sami Niemelä ◽  
Carl Fortelius

Abstract This paper presents a case study of a single cold air outbreak event with widespread convective precipitation over southern Finland on 25 May 2001. The purpose of the study is to investigate the applicability of the convection and condensation scheme of the High-Resolution Limited Area Model (HIRLAM) on meso-γ-scales. The study concentrates on the issue of grid-size-dependent convection parameterization. An explicit approach without the convection scheme is also examined. At the same time, the performance of an experimental nonhydrostatic version of HIRLAM is evaluated. Model simulations are conducted with three different horizontal grid spacings: 11, 5.6, and 2.8 km. Model results are compared to observed radar reflectivity data utilizing a radar simulation model, which calculates radar reflectivities from three-dimensional model output. The best results are obtained using nonhydrostatic dynamics and a grid-size-dependent convection scheme with a 5.6-km grid interval. However, even the best configuration still overestimates the area of strong reflectivity (intense precipitation). All the other combinations produce even stronger reflectivity. The grid-size-dependent convection parameterization is evidently beneficial at smaller grid spacings than 5.6 km. The nonhydrostatic model clearly outperforms its hydrostatic counterpart at the 5.6- and 2.8-km grid spacings, whereas with an 11-km grid interval, both models perform equally well.


2011 ◽  
Vol 68 (3) ◽  
pp. 515-532 ◽  
Author(s):  
Francesco A. Isotta ◽  
P. Spichtinger ◽  
U. Lohmann ◽  
K. von Salzen

Abstract A transient shallow-convection scheme is implemented into the general circulation model ECHAM5 and the coupled aerosol model HAM, developed at the Max Planck Institute for Meteorology in Hamburg. The shallow-convection scheme is extended to take the ice phase into account. In addition, a detailed double-moment microphysics approach has been added. In this approach, the freezing processes and precipitation formation are dependent on aerosols. Furthermore, in the scheme, tracers are transported and scavenged consistently as in the rest of the model. Results of a single-column model simulation for the Barbados Oceanography and Meteorology Experiment (BOMEX) campaign are compared with previously published large-eddy simulation (LES) results. Compared to the standard version, the global ECHAM5-HAM simulations with the newly implemented scheme show a decreased frequency of shallow convection in better agreement with LES. Less shallow convection is compensated by more stratus and stratocumulus. Deep and especially midlevel convection are markedly affected by those changes, which in turn influence high-level clouds. Generally, a better agreement with the observations can be obtained. For a better understanding of the scheme’s impact and to test different setting parameters, sensitivity analyses are performed. The mixing properties, cloud-base vertical velocity, and launching layer of the test parcel, respectively, are varied. In this context, results from simulations without shallow convection are also presented.


Author(s):  
Magnus Lindskog ◽  
Adam Dybbroe ◽  
Roger Randriamampianina

AbstractMetCoOp is a Nordic collaboration on operational Numerical Weather Prediction based on a common limited-area km-scale ensemble system. The initial states are produced using a 3-dimensional variational data assimilation scheme utilizing a large amount of observations from conventional in-situ measurements, weather radars, global navigation satellite system, advanced scatterometer data and satellite radiances from various satellite platforms. A version of the forecasting system which is aimed for future operations has been prepared for an enhanced assimilation of microwave radiances. This enhanced data assimilation system will use radiances from the Microwave Humidity Sounder, the Advanced Microwave Sounding Unit-A and the Micro-Wave Humidity Sounder-2 instruments on-board the Metop-C and Fengyun-3 C/D polar orbiting satellites. The implementation process includes channel selection, set-up of an adaptive bias correction procedure, and careful monitoring of data usage and quality control of observations. The benefit of the additional microwave observations in terms of data coverage and impact on analyses, as derived using the degree of freedom of signal approach, is demonstrated. A positive impact on forecast quality is shown, and the effect on the precipitation for a case study is examined. Finally, the role of enhanced data assimilation techniques and adaptions towards nowcasting are discussed.


2013 ◽  
Vol 70 (8) ◽  
pp. 2547-2565 ◽  
Author(s):  
Marie-Dominique Leroux ◽  
Matthieu Plu ◽  
David Barbary ◽  
Frank Roux ◽  
Philippe Arbogast

Abstract The rapid intensification of Tropical Cyclone (TC) Dora (2007, southwest Indian Ocean) under upper-level trough forcing is investigated. TC–trough interaction is simulated using a limited-area operational numerical weather prediction model. The interaction between the storm and the trough involves a coupled evolution of vertical wind shear and binary vortex interaction in the horizontal and vertical dimensions. The three-dimensional potential vorticity structure associated with the trough undergoes strong deformation as it approaches the storm. Potential vorticity (PV) is advected toward the tropical cyclone core over a thick layer from 200 to 500 hPa while the TC upper-level flow turns cyclonic from the continuous import of angular momentum. It is found that vortex intensification first occurs inside the eyewall and results from PV superposition in the thick aforementioned layer. The main pathway to further storm intensification is associated with secondary eyewall formation triggered by external forcing. Eddy angular momentum convergence and eddy PV fluxes are responsible for spinning up an outer eyewall over the entire troposphere, while spindown is observed within the primary eyewall. The 8-km-resolution model is able to reproduce the main features of the eyewall replacement cycle observed for TC Dora. The outer eyewall intensifies further through mean vertical advection under dynamically forced upward motion. The processes are illustrated and quantified using various diagnostics.


2018 ◽  
Vol 99 (7) ◽  
pp. 1415-1432 ◽  
Author(s):  
Yong Wang ◽  
Martin Belluš ◽  
Andrea Ehrlich ◽  
Máté Mile ◽  
Neva Pristov ◽  
...  

AbstractThis paper describes 27 years of scientific and operational achievement of Regional Cooperation for Limited Area Modelling in Central Europe (RC LACE), which is supported by the national (hydro-) meteorological services of Austria, Croatia, the Czech Republic, Hungary, Romania, Slovakia, and Slovenia. The principal objectives of RC LACE are to 1) develop and operate the state-of-the-art limited-area model and data assimilation system in the member states and 2) conduct joint scientific and technical research to improve the quality of the forecasts.In the last 27 years, RC LACE has contributed to the limited-area Aire Limitée Adaptation Dynamique Développement International (ALADIN) system in the areas of preprocessing of observations, data assimilation, model dynamics, physical parameterizations, mesoscale and convection-permitting ensemble forecasting, and verification. It has developed strong collaborations with numerical weather prediction (NWP) consortia ALADIN, the High Resolution Limited Area Model (HIRLAM) group, and the European Centre for Medium-Range Weather Forecasts (ECMWF). RC LACE member states exchange their national observations in real time and operate a common system that provides member states with the preprocessed observations for data assimilation and verification. RC LACE runs operationally a common mesoscale ensemble system, ALADIN–Limited Area Ensemble Forecasting (ALADIN-LAEF), over all of Europe for early warning of severe weather.RC LACE has established an extensive regional scientific and technical collaboration in the field of operational NWP for weather research, forecasting, and applications. Its 27 years of experience have demonstrated the value of regional cooperation among small- and medium-sized countries for success in the development of a modern forecasting system, knowledge transfer, and capacity building.


2005 ◽  
Vol 18 (19) ◽  
pp. 4046-4064 ◽  
Author(s):  
Guang J. Zhang ◽  
Mingquan Mu

Abstract This study presents the simulation of the Madden–Julian oscillation (MJO) in the NCAR CCM3 using a modified Zhang–McFarlane convection parameterization scheme. It is shown that, with the modified scheme, the intraseasonal (20–80 day) variability in precipitation, zonal wind, and outgoing longwave radiation (OLR) is enhanced substantially compared to the standard CCM3 simulation. Using a composite technique based on the empirical orthogonal function (EOF) analysis, the paper demonstrates that the simulated MJOs are in better agreement with the observations than the standard model in many important aspects. The amplitudes of the MJOs in 850-mb zonal wind, precipitation, and OLR are comparable to those of the observations, and the MJOs show clearly eastward propagation from the Indian Ocean to the Pacific. In contrast, the simulated MJOs in the standard CCM3 simulation are weak and have a tendency to propagate westward in the Indian Ocean. Nevertheless, there remain several deficiencies that are yet to be addressed. The time period of the MJOs is shorter, about 30 days, compared to the observed time period of 40 days. The spatial scale of the precipitation signal is smaller than observed. Examination of convective heating from both deep and shallow convection and its relationship with moisture anomalies indicates that near the mature phase of the MJO, regions of shallow convection developing ahead of the deep convection coincide with regions of positive moisture anomalies in the lower troposphere. This is consistent with the recent observations and theoretical development that shallow convection helps to precondition the atmosphere for MJO by moistening the lower troposphere. Sensitivity tests are performed on the individual changes in the modified convection scheme. They show that both change of closure and use of a relative humidity threshold for the convection trigger play important roles in improving the MJO simulation. Use of the new closure leads to the eastward propagation of the MJO and increases the intensity of the MJO signal in the wind field, while imposing a relative humidity threshold enhances the MJO variability in precipitation.


Időjárás ◽  
2021 ◽  
Vol 125 (4) ◽  
pp. 571-607
Author(s):  
André Simon ◽  
Martin Belluš ◽  
Katarína Čatlošová ◽  
Mária Derková ◽  
Martin Dian ◽  
...  

The paper presented is dedicated to the evaluation of the influence of various improvements to the numerical weather prediction (NWP) systems exploited at the Slovak Hydrometeorological Institute (SHMÚ). The impact was illustrated in a case study with multicell thunderstorms and the results were confronted with the reference analyses from the INCA nowcasting system, regional radar reflectivity data, and METEOSAT satellite imagery. The convective cells evolution was diagnosed in non-hydrostatic dynamics experiments to study weak mesoscale vortices and updrafts. The growth of simulated clouds and evolution of the temperature at their top were compared with the brightness temperature analyzed from satellite imagery. The results obtained indicated the potential for modeling and diagnostics of small-scale structures within the convective cloudiness, which could be related to severe weather. Furthermore, the non-hydrostatic dynamics experiments related to the stability and performance improvement of the time scheme led to the formulation of a new approach to linear operator definition for semi-implicit scheme (in text referred as NHHY). We demonstrate that the execution efficiency has improved by more than 20%. The exploitation of several high resolution measurement types in data assimilation contributed to more precise position of predicted patterns and precipitation representation in the case study. The non-hydrostatic dynamics provided more detailed structures. On the other hand, the potential of a single deterministic forecast of prefrontal heavy precipitation was not as high as provided by the ensemble system. The prediction of a regional ensemble system A-LAEF (ALARO Limited Area Ensemble Forecast) enhanced the localization of precipitation patterns. Though, this was rather due to the simulation of uncertainty in the initial conditions and also because of the stochastic perturbation of physics tendencies. The various physical parameterization setups of A-LAEF members did not exhibit a systematic effect on precipitation forecast in the evaluated case. Moreover, the ensemble system allowed an estimation of uncertainty in a rapidly developing severe weather case, which was high even at very short range.


2007 ◽  
Vol 135 (4) ◽  
pp. 1424-1438 ◽  
Author(s):  
Andrew R. Lawrence ◽  
James A. Hansen

Abstract An ensemble-based data assimilation approach is used to transform old ensemble forecast perturbations with more recent observations for the purpose of inexpensively increasing ensemble size. The impact of the transformations are propagated forward in time over the ensemble’s forecast period without rerunning any models, and these transformed ensemble forecast perturbations can be combined with the most recent ensemble forecast to sensibly increase forecast ensemble sizes. Because the transform takes place in perturbation space, the transformed perturbations must be centered on the ensemble mean from the most recent forecasts. Thus, the benefit of the approach is in terms of improved ensemble statistics rather than improvements in the mean. Larger ensemble forecasts can be used for numerous purposes, including probabilistic forecasting, targeted observations, and to provide boundary conditions to limited-area models. This transformed lagged ensemble forecasting approach is explored and is shown to give positive results in the context of a simple chaotic model. By incorporating a suitable perturbation inflation factor, the technique was found to generate forecast ensembles whose skill were statistically comparable to those produced by adding nonlinear model integrations. Implications for ensemble forecasts generated by numerical weather prediction models are briefly discussed, including multimodel ensemble forecasting.


2012 ◽  
Vol 140 (10) ◽  
pp. 3149-3162 ◽  
Author(s):  
Daan Degrauwe ◽  
Steven Caluwaerts ◽  
Fabrice Voitus ◽  
Rafiq Hamdi ◽  
Piet Termonia

Abstract Spectral limited-area models face a particular challenge at their lateral boundaries: the fields need to be made periodic. Boyd proposed a windowing-based method to improve the periodization and relaxation. In a companion paper, the implementation of this windowing method in the operational semi-implicit semi-Lagrangian spectral HARMONIE system was described and some first reproducibility tests, comparing this method to the old existing one, were presented. The present paper provides an in-depth study of the impact of this method for different configurations of the implementation. This is carried out in three steps in well-controlled experimental setups of increasing complexity. First, different aspects of Boyd’s method are analyzed in an idealized perfect-model test using a representative 1D shallow-water model. Second, the implementation is tested in an adiabatic 3D numerical weather prediction (NWP) model with perfect-model experiments. Finally, the impact of using Boyd’s method in a more operational-like NWP context is investigated as well. The presented tests show that, while the implementation of Boyd’s method is neutral in terms of scores, it is superior to the existing spline method in the case of strong dynamical forcings at the lateral boundaries.


Sign in / Sign up

Export Citation Format

Share Document