scholarly journals Combined use of volume radar observations and high-resolution numerical weather predictions to estimate precipitation at the ground: methodology and proof of concept

2019 ◽  
Vol 12 (10) ◽  
pp. 5669-5684 ◽  
Author(s):  
Tony Le Bastard ◽  
Olivier Caumont ◽  
Nicolas Gaussiat ◽  
Fatima Karbou

Abstract. The extrapolation of the precipitation to the ground from radar reflectivities measured at the beam altitude is one of the most delicate phases of radar data processing for producing quantitative precipitation estimations (QPEs) and remains a major scientific issue. In many operational meteorological services such as Météo-France, a vertical profile of reflectivity (VPR) correction is uniformly applied over a large part or the entire radar domain. This method is computationally efficient, and the overall bias induced by the bright band is most of the time well corrected. However, this way of proceeding is questionable in situations with high spatial and vertical variability of precipitation (during the passage of a cold front or in a complex terrain, for example). This study initiates from two statements: first, radars provide information on precipitation with a high spatio-temporal resolution but still require VPR corrections to extrapolate rain rates at the ground level. Second, the horizontal resolution of some numerical weather prediction (NWP) models is now comparable with the radar one, and their dynamical core and microphysics schemes allow the production of realistic simulations of VPRs. The present paper proposes a new approach to assess surface rainfall from radar reflectivity aloft by exploiting simulated VPRs and rainfall forecasts from the high-resolution NWP model AROME-NWC. To our knowledge, this is the first time that simulated precipitation profiles from an NWP model are used to derive radar QPEs. The implementation of the new method on two stratiform situations provided significant improvements on the hourly and 6 h accumulations compared to the operational QPEs, showing the relevance of this new approach.

2019 ◽  
Author(s):  
Tony Le Bastard ◽  
Olivier Caumont ◽  
Nicolas Gaussiat ◽  
Fatima Karbou

Abstract. The extrapolation of the precipitation to the ground from radar reflectivities measured at the beam altitude is one of the most delicate phases of radar data processing for producing Quantitative Precipitation Estimations (QPEs) and remains a major scientific issue. In many operational meteorological services such as Météo-France, a Vertical Profile of Reflectivity (VPR) correction is uniformly applied over a large part or the entire radar domain. This method is computationally efficient and the overall bias induced by the bright band is most of the time well corrected. However, this way of proceeding is questionable in situations with high spatial and vertical variability of precipitation (during the passage of a cold front or in a complex terrain, for example). This study initiates from two statements: first, radars provide information on precipitation with a high spatio-temporal resolution but still require VPR corrections to extrapolate rain rates at the ground level. Second, the horizontal resolution of some Numerical Weather Prediction (NWP) models is now comparable with the radar one and their dynamical core and microphysics schemes allow to produce realistic simulations of VPRs. The present paper proposes a new approach to assess surface rainfall from radar reflectivity aloft by exploiting simulated VPRs and rainfall forecasts from the high resolution NWP model AROME-NWC. To our knowledge, this is the first time that simulated precipitation profiles from a NWP model are used to derive radar QPEs. The implementation of the new method on two stratiform situations provided significant improvements on the hourly and 6-h accumulations compared to the operational QPEs, showing the relevance of this new approach.


2020 ◽  
Author(s):  
Xavier Lapillonne ◽  
William Sawyer ◽  
Philippe Marti ◽  
Valentin Clement ◽  
Remo Dietlicher ◽  
...  

<p>The ICON modelling framework is a unified numerical weather and climate model used for applications ranging from operational numerical weather prediction to low and high resolution climate projection. In view of further pushing the frontier of possible applications and to make use of the latest evolution in hardware technologies, parts of the model were recently adapted to run on heterogeneous GPU system. This initial GPU port focus on components required for high-resolution climate application, and allow considering multi-years simulations at 2.8 km on the Piz Daint heterogeneous supercomputer. These simulations are planned as part of the QUIBICC project “The Quasi-Biennial Oscillation (QBO) in a changing climate”, which propose to investigate effects of climate change on the dynamics of the QBO.</p><p>Because of the low compute intensity of atmospheric model the cost of data transfer between CPU and GPU at every step of the time integration would be prohibitive if only some components would be ported to the accelerator. We therefore present a full port strategy where all components required for the simulations are running on the GPU. For the dynamics, most of the physical parameterizations and infrastructure code the OpenACC compiler directives are used. For the soil parameterization, a Fortran based domain specific language (DSL) the CLAW-DSL has been considered. We discuss the challenges associated to port a large community code, about 1 million lines of code, as well as to run simulations on large-scale system at 2.8 km horizontal resolution in terms of run time and I/O constraints. We show performance comparison of the full model on CPU and GPU, achieving a speed up factor of approximately 5x, as well as scaling results on up to 2000 GPU nodes. Finally we discuss challenges and planned development regarding performance portability and high level DSL which will be used with the ICON model in the near future.</p>


Author(s):  
Nils P. Wedi

The steady path of doubling the global horizontal resolution approximately every 8 years in numerical weather prediction (NWP) at the European Centre for Medium Range Weather Forecasts may be substan- tially altered with emerging novel computing architectures. It coincides with the need to appropriately address and determine forecast uncertainty with increasing resolution, in particular, when convective-scale motions start to be resolved. Blunt increases in the model resolution will quickly become unaffordable and may not lead to improved NWP forecasts. Consequently, there is a need to accordingly adjust proven numerical techniques. An informed decision on the modelling strategy for harnessing exascale, massively parallel computing power thus also requires a deeper understanding of the sensitivity to uncertainty—for each part of the model—and ultimately a deeper understanding of multi-scale interactions in the atmosphere and their numerical realization in ultra-high-resolution NWP and climate simulations. This paper explores opportunities for substantial increases in the forecast efficiency by judicious adjustment of the formal accuracy or relative resolution in the spectral and physical space. One path is to reduce the formal accuracy by which the spectral transforms are computed. The other pathway explores the importance of the ratio used for the horizontal resolution in gridpoint space versus wavenumbers in spectral space. This is relevant for both high-resolution simulations as well as ensemble-based uncertainty estimation.


2017 ◽  
Vol 17 (23) ◽  
pp. 14253-14273 ◽  
Author(s):  
Daniel Wolfensberger ◽  
Auguste Gires ◽  
Ioulia Tchiguirinskaia ◽  
Daniel Schertzer ◽  
Alexis Berne

Abstract. The framework of universal multifractals (UM) characterizes the spatio-temporal variability in geophysical data over a wide range of scales with only a limited number of scale-invariant parameters. This work aims to clarify the link between multifractals (MFs) and more conventional weather descriptors and to show how they can be used to perform a multi-scale evaluation of model data. The first part of this work focuses on a MF analysis of the climatology of precipitation intensities simulated by the COSMO numerical weather prediction model. Analysis of the spatial structure of the MF parameters, and their correlations with external meteorological and topographical descriptors, reveals that simulated precipitation tends to be smoother at higher altitudes, and that the mean intermittency is mostly influenced by the latitude. A hierarchical clustering was performed on the external descriptors, yielding three different clusters, which correspond roughly to Alpine/continental, Mediterranean and temperate regions. Distributions of MF parameters within these three clusters are shown to be statistically significantly different, indicating that the MF signature of rain is indeed geographically dependent. The second part of this work is event-based and focuses on the smaller scales. The MF parameters of precipitation intensities at the ground are compared with those obtained from the Swiss radar composite during three events corresponding to typical synoptic conditions over Switzerland. The results of this analysis show that the COSMO simulations exhibit spatial scaling breaks that are not present in the radar data, indicating that the model is not able to simulate the observed variability at all scales. A comparison of the operational one-moment microphysical parameterization scheme of COSMO with a more advanced two-moment scheme reveals that, while no scheme systematically outperforms the other, the two-moment scheme tends to produce larger extreme values and more discontinuous precipitation fields, which agree better with the radar composite.


2021 ◽  
Vol 14 (1) ◽  
pp. 42
Author(s):  
Bojun Zhu ◽  
Zhaoxia Pu ◽  
Agie Wandala Putra ◽  
Zhiqiu Gao

Accurate high-resolution precipitation forecasts are critical yet challenging for weather prediction under complex topography or severe synoptic forcing. Data fusion and assimilation aimed at improving model forecasts, as one possible approach, has gained increasing attention in past decades. This study investigates the influence of the observations from a C-band Doppler radar over the west coast of Sumatra on high-resolution numerical simulations of precipitation around its vicinity under the Madden–Julian oscillation (MJO) in January and February 2018. Cases during various MJO phases were selected for simulations with an advanced research version of the weather research and forecasting (WRF) model at a cloud-permitting scale (~3 km). A 3-dimensional variational (3DVAR) data assimilation method and a hybrid three-dimensional ensemble–variational data assimilation (3DEnVAR) method, based on the NCEP Gridpoint Statistical Interpolation (GSI) assimilation system, were used to assimilate the radar reflectivity and the radial velocity data. The WRF-simulated precipitation was validated with the Integrated Multi-satellitE Retrievals for GPM (IMERG) precipitation data, and the fractions skill score (FSS) was calculated in order to evaluate the radar data impacts objectively. The results show improvements in the simulated precipitation with hourly radar data assimilation 6 h prior to the simulations. The modifications with assimilation were validated through the observation departure and moist convection. It was found that forecast improvements are relatively significant when precipitation is more related to local-scale convection but rather small when the background westerly wind is strong under the MJO active phase. The additional simulation experiments, under a 1- or 2-day assimilation cycle, indicate better improvements in the precipitation simulation with 3DEnVAR radar assimilation than those with the 3DVAR method.


2021 ◽  
Vol 38 (5) ◽  
pp. 737-754
Author(s):  
Guifu Zhang ◽  
Jidong Gao ◽  
Muyun Du

AbstractMany weather radar networks in the world have now provided polarimetric radar data (PRD) that have the potential to improve our understanding of cloud and precipitation microphysics, and numerical weather prediction (NWP). To realize this potential, an accurate and efficient set of polarimetric observation operators are needed to simulate and assimilate the PRD with an NWP model for an accurate analysis of the model state variables. For this purpose, a set of parameterized observation operators are developed to simulate and assimilate polarimetric radar data from NWP model-predicted hydrometeor mixing ratios and number concentrations of rain, snow, hail, and graupel. The polarimetric radar variables are calculated based on the T-matrix calculation of wave scattering and integrations of the scattering weighted by the particle size distribution. The calculated polarimetric variables are then fitted to simple functions of water content and volume-weighted mean diameter of the hydrometeor particle size distribution. The parameterized PRD operators are applied to an ideal case and a real case predicted by the Weather Research and Forecasting (WRF) model to have simulated PRD, which are compared with existing operators and real observations to show their validity and applicability. The new PRD operators use less than one percent of the computing time of the old operators to complete the same simulations, making it efficient in PRD simulation and assimilation usage.


2008 ◽  
Vol 136 (3) ◽  
pp. 945-963 ◽  
Author(s):  
Jidong Gao ◽  
Ming Xue

Abstract A new efficient dual-resolution (DR) data assimilation algorithm is developed based on the ensemble Kalman filter (EnKF) method and tested using simulated radar radial velocity data for a supercell storm. Radar observations are assimilated on both high-resolution and lower-resolution grids using the EnKF algorithm with flow-dependent background error covariances estimated from the lower-resolution ensemble. It is shown that the flow-dependent and dynamically evolved background error covariances thus estimated are effective in producing quality analyses on the high-resolution grid. The DR method has the advantage of being able to significantly reduce the computational cost of the EnKF analysis. In the system, the lower-resolution ensemble provides the flow-dependent background error covariance, while the single-high-resolution forecast and analysis provides the benefit of higher resolution, which is important for resolving the internal structures of thunderstorms. The relative smoothness of the covariance obtained from the lower 4-km-resolution ensemble does not appear to significantly degrade the quality of analysis. This is because the cross covariance among different variables is of first-order importance for “retrieving” unobserved variables from the radar radial velocity data. For the DR analysis, an ensemble size of 40 appears to be a reasonable choice with the use of a 4-km horizontal resolution in the ensemble and a 1-km resolution in the high-resolution analysis. Several sensitivity tests show that the DR EnKF system is quite robust to different observation errors. A 4-km thinned data resolution is a compromise that is acceptable under the constraint of real-time applications. A data density of 8 km leads to a significant degradation in the analysis.


2021 ◽  
Vol 94 (2) ◽  
pp. 237-249
Author(s):  
Martin Novák

The article includes a summary of basic information about the Universal Thermal Climate Index (UTCI) calculation by the numerical weather prediction (NWP) model ALADIN of the Czech Hydrometeorological Institute (CHMI). Examples of operational outputs for weather forecasters in the CHMI are shown in the first part of this work. The second part includes results of a comparison of computed UTCI values by ALADIN for selected place with UTCI values computed from real measured meteorological data from the same place.


Sign in / Sign up

Export Citation Format

Share Document