The dependence of precipitation growth and decay on the orientation of orographic features, mesoscale flow and freezing‐level height is quantified using a 10‐year archive of composite weather radar images over the Swiss Alpine region. The mesoscale flow is described by the motion of radar precipitation echoes, computed through variational echo tracking, while the freezing‐level height is extracted from the analyses of the numerical weather prediction model COSMO. On the northern side of the Alps, the areas of growth are generally observed on the upwind slopes of the Alpine chain, while on the southern side their location also depends on factors other than the motion of precipitation, such as the convergence of low‐level flows. On the other hand, the decay of precipitation is generally found in the inner Alpine valleys and on the downwind slopes of the Alpine chain. Compared to situations characterized by low freezing level, the areas of precipitation growth penetrate more into the mountain range with high‐freezing‐level conditions. When considering a time lag of 1 h and specific flow conditions, the systematic growth and decay of precipitation can explain up to 30–40% of its total variability over the orography. The relative contribution to the total variability is lower with high freezing level because the non‐systematic variability is larger. The statistical analysis of precipitation growth and decay using large archives of composite radar images highlights the mesoscale flow conditions and geographical locations most prone to orographic precipitation enhancement, but also has potential to improve existing precipitation nowcasting systems in mountainous regions and to provide a basis for the verification of systematic biases of numerical weather prediction models.
Machine learning algorithms are trained on a 10-yr archive of composite weather radar images in the Swiss Alps to nowcast precipitation growth and decay in the next few hours in moving coordinates (Lagrangian frame). The hypothesis of this study is that growth and decay is more predictable in mountainous regions, which represent a potential source of practical predictability by machine learning methods. In this paper, artificial neural networks (ANN) are employed to learn the complex nonlinear dependence relating the growth and decay to the input predictors, which are geographical location, mesoscale motion vectors, freezing level height, and time of the day. The average long-term growth and decay patterns are effectively reproduced by the ANN, which allows exploring their climatology for any combination of predictors. Due to the low intrinsic predictability of growth and decay, its prediction in real time is more challenging, but is substantially improved when adding persistence information to the predictors, more precisely the growth and decay and precipitation intensity in the immediate past. The improvement is considerable in mountainous regions, where, depending on flow direction, the root-mean-square error of ANN predictions can be 20%–30% lower compared with persistence. Because large uncertainty is associated with precipitation forecasting, deterministic machine learning predictions should be coupled with a model for the predictive uncertainty. Therefore, we consider a probabilistic perspective by estimating prediction intervals based on a combination of quantile decision trees and ANNs. The probabilistic framework is an attempt to address the problem of conditional bias, which often characterizes deterministic machine learning predictions obtained by error minimization.
In this study, we develop and compare satellite rainfall retrievals based on generalized linear models and artificial neural networks. Both approaches are used in classification mode in a first step to identify the precipitating areas (precipitation detection) and in regression mode in a second step to estimate the rainfall intensity at the ground (rain rate). The input predictors are geostationary satellite infrared (IR) brightness temperatures and Satellite Application Facility (SAF) nowcasting products which consist of cloud properties, such as cloud top height and cloud type. Additionally, a set of auxiliary location-describing input variables is employed. The output predictand is the ground-based instantaneous rain rate provided by the European-scale radar composite OPERA, that was additionally quality-controlled. We compare our results to a precipitation product which uses a single infrared (IR) channel for the rainfall retrieval. Specifically, we choose the operational PR-OBS-3 hydrology SAF product as a representative example for this type of approach. With generalized linear models, we show that we are able to substantially improve in terms of hits by considering more IR channels and cloud property predictors. Furthermore, we demonstrate the added value of using artificial neural networks to further improve prediction skill by additionally reducing false alarms. In the rain rate estimation, the indirect relationship between surface rain rates and the cloud properties measurable with geostationary satellites limit the skill of all models, which leads to smooth predictions close to the mean rainfall intensity. Probability matching is explored as a tool to recover higher order statistics to obtain a more realistic rain rate distribution.
Archives of composite weather radar images represent an invaluable resource to study the predictability of precipitation. In this paper, we compare two distinct approaches to construct empirical low-dimensional attractors from radar precipitation fields. In the first approach, the phase space variables of the attractor are defined using the domain-scale statistics of precipitation fields, such as the mean precipitation, fraction of rain, and spatial and temporal correlations. The second type of attractor considers the spatial distribution of precipitation and is built by principal component analysis (PCA). For both attractors, we investigate the density of trajectories in phase space, growth of errors from analogue states, and fractal properties. To represent different scales and climatic and orographic conditions, the analyses are done using multi-year radar archives over the continental United States ( km2, 21 years) and the Swiss Alpine region ( km2, 6 years).
Machine learning algorithms are trained on a 10-yr archive of composite weather radar images in the Swiss Alps to nowcast precipitation growth and decay in the next few hours in moving coordinates (Lagrangian frame). The hypothesis of this study is that growth and decay is more predictable in mountainous regions, which represent a potential source of practical predictability by machine learning methods. In this paper, artificial neural networks (ANN) are employed to learn the complex nonlinear dependence relating the growth and decay to the input predictors, which are geographical location, mesoscale motion vectors, freezing level height, and time of the day. The average long-term growth and decay patterns are effectively reproduced by the ANN, which allows exploring their climatology for any combination of predictors. Due to the low intrinsic predictability of growth and decay, its prediction in real time is more challenging, but is substantially improved when adding persistence information to the predictors, more precisely the growth and decay and precipitation intensity in the immediate past. The improvement is considerable in mountainous regions, where, depending on flow direction, the root-mean-square error of ANN predictions can be 20%–30% lower compared with persistence. Because large uncertainty is associated with precipitation forecasting, deterministic machine learning predictions should be coupled with a model for the predictive uncertainty. Therefore, we consider a probabilistic perspective by estimating prediction intervals based on a combination of quantile decision trees and ANNs. The probabilistic framework is an attempt to address the problem of conditional bias, which often characterizes deterministic machine learning predictions obtained by error minimization.
Archives of composite weather radar images represent an invaluable resource to study the predictability of precipitation. In this paper, we compare two distinct approaches to construct empirical low-dimensional attractors from radar precipitation fields. In the first approach, the phase space dimensions of the attractor are defined using the domain-scale statistics of precipitation fields, such as the mean precipitation, fraction of rain, spatial and temporal correlations. The second type of attractor considers the spatial distribution of precipitation and is built by principal component analysis (PCA). For both attractors, we investigate the density of trajectories in phase space, growth of errors from analogue states, and fractal properties. To represent different scales, climatic and orographic conditions, the analyses are done using multi-year radar archives over the continental United States (≈4000 x 4000 km2, 21 years) and the Swiss Alpine region (≈500 x 500 km2, 6 years).