Abstract
Accurate and timely prediction of short-term rainfall is crucial for reducing the damages caused by heavy rainfall events. Therefore, various precipitation nowcasting models have been proposed. However, the performance of these models still remains limited. In particular, the current operational precipitation nowcasting method, which is based on radar echo tracking, such as the McGill Algorithm for Precipitation Nowcasting by Lagrangian Extrapolation (MAPLE), has a critical drawback when predicting newly developed or decayed precipitation fields. Recently proposed deep learning models, such as the U-Net convolutional neural network outperform the models based on radar echo tracking. However, these models are unsuitable for operational precipitation nowcasting due to their blurry predictions over longer lead times. To address these blurry effects and enhance the performance of U-Net-based rainfall prediction, we propose a blended model that combines a partial differential equation (PDE) model based on fluid dynamics with the U-Net model. The evaluation of the forecast skill, based on both qualitative and quantitative methods for 0–3-h lead times, demonstrates that the blended model provides less blurry and more accurate rainfall predictions compared with the U-Net and partial differential equation models. This indicates the potential to enhance the field of very short-term rainfall prediction. Additionally, we also evaluated the monthly-averaged forecast skills for different seasons, and confirmed the operational feasibility of the blended model, which contributes to the performance enhancement of operational nowcasting.
1 Introduction
Accurate and timely predictions of short-term rainfall are crucial for minimizing damages caused by heavy rainfall events. These damages include economic and human losses owing to flooding, as well as disruptions in daily life, such as transportation, including air travel. For example, during the summer of 2020, South Korea experienced destructive floods, causing economic losses amounting to billions of dollars (). To mitigate such damages, precise short-term rainfall predictions within 0–6-h lead times are necessary (e.g., ; ).
Radar echo tracking methods utilizing weather radar data, such as the McGill Algorithm for Precipitation nowcasting by Lagrangian Extrapolation (MAPLE), have been employed in operational short-term rainfall prediction models (e.g., ; ; ; ; ), and Python open-source libraries for precipitation nowcasting based on such methods have been also proposed, such as rainymotion v1.0 () and PySTEPS (). Additionally, more sophisticated models have been proposed by considering advection-diffusion equations and non-stationary velocity fields, and such models outperform MAPLE, which solely handles advection (). However, even with the ability of these models to predict the movement of precipitation patterns, limitations exist in their predictive performance owing to the difficulty in considering dynamical processes, such as precipitation generation and dissipation (e.g., ). In particular, a significant portion of rainfall exhibits a stationary front, thereby hindering rainfall prediction using the radar echo tracking methods.
Numerical weather prediction (NWP) models can forecast the dynamical processes, such as growth and decay of precipitation fields, by incorporating atmospheric physics (e.g., ; ; ; ; ). Compared with the radar echo tracking methods, NWP can model the dynamical processes, including the growth and decay of precipitation fields. However, their spin-up time has limitations in achieving accurate performance within lead times of a few hours (e.g., ).
Recently, deep learning-based models have been studied to predict the non-linear evolution of precipitation patterns (e.g., ; ; ; ; ; ; ; ; ; ; ; ). These models are data-driven and physics-free, as they can be trained using high-resolution radar data that contain spatiotemporal information on instantaneous rainfall. S-band radar has a spatiotemporal resolution of approximately 10-min intervals and 1-km scale (), enabling sufficiently high-resolution predictions within 0–6-h lead times. In particular, models based on convolutional long short-term memory (Conv-LSTM) and U-Net have previously been proposed, wherein multiple past radar images are used as input data to train the model and predict future timestep images (; ; ; ; ; ; ; ; ). These models learn the complex evolution of rainfall patterns over time, resulting in better predictive performance than baseline models, such as NWP, radar echo tracking, and Eulerian persistence. Such convolutional models, however, generate blurry prediction images with increasing prediction time; thus, adopting such models for operational precipitation nowcasting is unlikely. To resolve the issue of blurry images, Generative Adversarial Network (GAN)-based generative models have been proposed in previous studies to obtain predictions without blurry effects (e.g., ; ; ). In particular, proposed a novel idea to produce spatiotemporally consistent prediction by employing spatial and temporal discriminators, which discourage blurry and temporally inconsistent predictions, respectively. Nevertheless, these predictive performances are not sustained over long periods (i.e., the prediction timescale is currently limited to around 0–2-h lead times) and exhibit decreased accuracy in predicting heavy rainfall events (e.g., ).
Studies have been conducted to improve the performance of convolutional models, such as U-Net. Various types of loss functions have been applied, leading to significant performance improvements (e.g., ). For instance, the utilization of the logcosh loss has positively impacted the predictive performance compared with that of simple MSE loss ( and the references therein). In addition, some of recent studies have also employed the loss functions based on Critical Success Index (CSI), which is used to measure the agreement between predicted and observed precipitation (; ). As demonstrated by , the application of CSI-based loss is a substantial factor aiding in the development of a model that can outperform conventional convolutional models and generative models in terms of the rainfall forecasting accuracy. Along with the loss function, examined the dependence on the training method by employing pre-training and fine-tuning steps. According to their findings, the training method has a significant impact on improving the prediction accuracy. Furthermore, in recent studies, both radar and satellite data or surface observations have been used as training data for U-Net models. The incorporation of multiple observational data has shown promising results in improving model performance (e.g., ; ; ). These improvements suggest that if the blurry effects of U-Net models can be addressed, U-Net-based precipitation prediction models can be competitive.
In this study, we propose an approach to resolve the blurry predictions produced by U-Net-based model. Assuming that the evolution of precipitation patterns is primarily driven by certain processes, such as advection and temporal variations in rainfall intensity, U-Net-based models suffer from difficulties in tracking the advection of precipitation owing to blurry effects, whereas models based on the advection-diffusion equation struggle to predict the complex dynamical processes, such as the growth and/or decay of precipitation. To overcome these limitations, we first employed a model that utilizes detailed solutions of partial differential equations (PDEs) associated with fluid dynamics to capture the motion of precipitation fields. We then developed the model by combining the PDE-based and U-Net-based predictions, resulting in less blurry and more precise forecasts. The performance of the constructed blended model was quantitatively and qualitatively evaluated using precipitation cases on the Korean Peninsula. In particular, previous studies regarding the development of deep learning models for precipitation nowcasting have focused on constructing and evaluating deep learning models at the research level. However, in this study, we aim to propose practical approaches that maintain the prediction tendencies of PDE-based models already used in short-term operational forecasting, such as MAPLE, and enhance their performance by incorporating deep learning models into the current approaches of operational nowcasting. To explore the potential of this blended model for operational nowcasting, we also compared them with recently proposed deep learning-based models through baseline models, such as Eulerian persistence and U-Net.
The paper is organized as follows. Section 2 describes the dataset and model that were used in this study. Section 3 presents the performance of the model combining U-Net architecture and forecast based on PDEs. Section 4 provides a brief summary of the study.
2 Methods
2.1 Radar dataset
Korean Meteorological Administration (KMA) has provided radar data based on Hybrid Surface Rainfall (HSR) method since 2018. The HSR method is characterized by the synthesis of reflectivity at the hybrid surface that is unaffected by the ground clutter, beam blockage, non-meteorological echoes, and bright band (e.g., ; ; ). The size of data is (or, ) and the spatial and temporal resolutions of the HSR data are 0.5 km and 5 min, respectively. From the radar reflectivity factor (), precipitation () can be estimated through the Z-R relationship of KMA, .
For training the model described in Section 3, the HSR data from 2018, 2019, and 2021 were used and the prediction results were validated using the HSR data from 2020. We then tested and evaluated the performance of the model using the HSR data from 2022. The preprocessing step is summarized as follows: 1) max-pooling to reduce the spatial resolution of the original data from 0.5 km to 1km; 2) producing a cropped image with a final size of ; 3) applying the transformation, (), to address the data imbalance in precipitation; 4) Generating input data by stacking time-sequential images. When stacking, three images were used with a temporal resolution of 10 min.
The types of heavy rainfall events occurring in the Korean Peninsula and their characteristics have been extensively examined using self-organization map (SOM) and K-means clustering (; ). In previous studies, the classification of heavy rainfall events through the SOM method is explained based on the diverse synoptic weather patterns of the rainy season of Korea (i.e., typically from June to September), and a substantial fraction of such heavy rainfall events do not occur throughout the Korean Peninsula but rather in localized regions. Figure 1 shows two representative examples of heavy rainfall events typically occurring in South Korea. The first example shows a quasi-stationary monsoon front that occurred in the central region of the Korean Peninsula (Central case, hereafter). These types of rainfall typically occur due to a clash between the air masses around the Korean peninsula (e.g., ). The second example shows the heavy rainfall in the southern coastal area of Korean peninsula (Southern case, hereafter). Owing to the presence of a southwesterly low-level jet, which contains a large amount of moisture, such precipitation events are typically induced (e.g., ).
FIGURE 1
2.2 Structure of blended model
In this section, we present the idea of blended model, wherein the forecast based on the fluid dynamics by solving the time-dependent PDEs is combined with that based on U-Net convolutional neural network. Figure 2 shows a structure of the blended model and the details of each steps are described in subsections 2.2.1–2.2.3.
FIGURE 2
2.2.1 Forecasts based on partial differential equations
We adopted the forecast models based on PDEs (PDE model, hereafter) proposed by . Under the assumption that the motion field of precipitation follows a fluid equation, the model provides nowcasting by solving the advection-diffusion equation and Burgers’ equation in a time-dependent manner. The steps for obtaining nowcasts are as follows: 1) calculating the initial motion field, and 2) solving time-dependent PDEs to obtain nowcasts.
The initial motion field in space at time is primarily calculated by the optical flow algorithm. In the optical flow method, the motion field is extracted from two adjacent images in the time domain. Under the assumption that the pixel intensities of an object remain constant, the rainfall intensity can be expressed using a Taylor series:when ignoring high-order terms. Dividing the above equation by , we obtainwhere and represent the basic elements of the vector field. As two unknown variables (, ) cannot be solved using only one equation, relevant regularization force is required to solve the equation. The total variation-L1 algorithm (TV-L1; ) is used to determine the optical flow by minimizing a regularization force, , defined as follows:where and are the rainfall intensities of two different timesteps. The first term is the optical flow constraint, which assumes brightness constancy during motion, and the second term represents the smoothness of displacement fields. TV-L1 is available from the Python public library OpenCV (https://opencv.org).
The motion field obtained by Eqs 1−3 is then updated through the Burgers’ equation, which is expressed as follows,where is the kinematic viscosity. The fluid system with a larger viscosity has a smoother motion field. Using the updated motion field, nowcasting is obtained by solving the advection-diffusion equation,where the first term is responsible for advection and the second term is used to describe the diffusion effect with a diffusion coefficient . Notably, and are variables varying on the spatial and temporal domains, and were assumed as constants. The values and were used, which were empirically determined by . Additionally, to numerically calculate the derivative, we used the centered difference formula (), which is expressed as follows:
The Eq. 7 and Eq. 8 were used for solving the partial differential Eqs 4−6.
2.2.2 Forecasts based on the U-Net convolutional neural network
A deep-learning model called RainNet v1.0 was recently proposed by for precipitation nowcasting, wherein the U-Net architecture is adopted, as introduced by . RainNet v1.0 contains four downsampling and four upsampling steps, as shown in Figure 3. During the downsampling steps, the resolution of the input data is reduced by merging every four adjacent pixels, whereras the number of latent features per pixel, representing the pixel’s vector dimension, is doubled. Conversely, in the upsampling steps, the resolution of the input data is doubled, whereas the number of latent features per pixel is halved. To ensure coherence between the downsampling and upsampling steps, skip connections are employed, enabling the intermediate output of each downsampling step (reflecting the vector representation of each pixel) that contributes to the determination of the output of the corresponding upsampling step. The convolutional layers with up to 1,024 filters, kernel sizes of and , and rectified linear unit (ReLU) activation functions were used. Additionally, the dropout layers with the dropout rate = 0.5 were implemented to prevent overfitting. Note that the model parameters were tuned by testing within the following parameter ranges: 1) the number of downsampling and upsampling steps: 3–5; 2) the size of the kernel: , , and ; 3) Dropout rates: 0–0.5. For a more comprehensive understanding of the architecture, refer to .
FIGURE 3
By employing RainNet v1.0, we obtained the nowcasting outputs, as follows:where and denote the radar ground truth and prediction, respectively. As a loss function, we adopted the logcosh loss function, , as follows:where , and denote the precipitation fields of ground truth, prediction, and the number of pixels, respectively. To predict the precipitation at through Eqs 9, 10, four latest images with interval were used as the input. Notably the spatial resolutions of the input and output are identical, as we used a 1-km interval and of 10 min. We optimized the hyperparameters using 30 epochs with ranges of batch sizes and learning rates, 8–32 and , respectively. As a set of hyperparameters, we employed the Adam optimizer, with a batch size of 8 and learning rate of .
2.2.3 Blending the PDE and U-Net models
By combining the forecast from the PDE model () and that from the U-Net model (), the final nowcasting at time , , can be generated. In this study, the blended model is simply defined as . The implication of the characteristics of the blended model is provided below. The time evolution of rainfall intensity, , is expressed as follows:
Here, we first assumed that the typical evolution of precipitation is expressed as the advection of the precipitation field and variation of intensity changes owing to the dynamic processes. Under such an assumption, we examined three cases: 1) advection dominant regime: ; 2) dynamical processes dominant regime: ; and 3) intermediate regime: .
In case that the precipitation field obtained by the PDE model, , dominates over , expressed by the Eq. 11 can be approximately expressed as follows,
Equation 12 indicates the temporal features of U-Net forecasts, and is simultaneously involved when the evolution of precipitation is mainly governed by the evolution of . In particular, the growth and decay of precipitation fields can be considered, that are hardly predicted through the PDE model.
When the rainfall intensity obtained by the U-Net model, , is more crucial, can be simplified as
We interpret the Eq. 13 that time variation of can be considered when the U-Net model, is a dominant forecaster. The main advantage is that the advection features of and precipitation evolution occurring at the small scale could be predicted, which are poorly predicted through the U-Net model.
When the contribution of the two models is comparable, the precipitation field is evolved through both and ,
The form of shown in Eq. 14 demonstrates that the precipitation field is evolved by the average evolution of two models.
Notably, during real precipitation evolution, three regimes are simultaneously included. Thus, the blended model could substantially enhance the performance of precipitation nowcasting.
2.2.4 Reference model: Eulerian persistence
As a reference model, we employed Eulerian persistence. Eulerian persistence model (Persistence, hereafter) ensures that the most recent observation is frozen:where and are the observed and predicted precipitation fields at times and , respectively. Although this model (Eqs 15, 16) is simple and does not consider the advection or generation and dissipation of the precipitation field, it is powerful for very short lead times. Thus, it has been employed as a baseline model in the recent deep learning-based studies for nowcasting (e.g., ; ; ; ).
The purpose of using Persistence can be summarized as follows. The comparison between the PDE model and Persistence enables the evaluation of the accuracy of precipitation nowcasting, considering advection in the precipitation field. By comparing Persistence with the U-Net and the blended models, we can observe the impact of the prediction of the change in precipitation intensity on the accuracy of precipitation nowcasting. Additionally, the comparison between this study and recent studies of deep learning-based precipitation nowcasting can be conducted using Persistence.
3 Results and discussion
This section describes the performance comparison between PDE, U-Net, and blended models. In particular, we report the prediction performance for heavy rainfall events with an accumulated rainfall of within 1 hour, as such events could lead to destructive floods. The evaluation focused on heavy rainfall events during 2022. Both qualitative and quantitative methods were employed to assess the model performance in predicting the motion of precipitation fields and changes in their intensities. Subsequently, we discuss the operational potential of the blended model based on the monthly averaged performance of the blended model.
Figures 4, 5 show the prediction results or precipitation in the central and southern regions of the Korean Peninsula. As demonstrated in previous studies using U-Net (e.g., ), the U-Net model tends to exhibit significant blurry effects and slightly underestimate the rainfall intensities. Such drawbacks of U-Net predictions are particularly important in certain cases, such as in the example shown in Figure 4, where the stationary rainfall front persists strongly. Conversely, the PDE model overestimates both the spatial extent and rainfall intensities as it fails to accurately simulate the processes such as precipitation dissipation. In particular, PDE-based predictions could provide inaccurate forecasts when encountering precipitation cases with evolving rainfall that dissipates, as shown in the example of Figure 5 (i.e., false alarm rate could increase). In comparison, the blended model performs better than the U-Net and PDE models in tracking the motion of precipitation patterns and accurately predicting dynamical processes, including the dissipation of precipitation and expansion of precipitation areas, in both the central and southern regions.
FIGURE 4
FIGURE 5
We also assessed the robustness of the blended model’s performance by predicting heavy rainfall events with characteristics similar to those in Figures 4, 5. The two right panels in Figure 6 demonstrate that the blended model effectively predicts an elongated structure, including convective and stratiform rains, in the Central area. The two left panels in Figure 6 show the performance in predicting the development of the precipitation field in the Southern coastal area. Notably, the blended model accurately forecasts the injection of precipitation from the south-western area, as illustrated in the case at 1200 UTC on 25 April 2022.
FIGURE 6
The blurry effects of the precipitation prediction field were quantitatively measured using power spectral density (PSD) calculations during the 0–3-h lead times. The radially averaged PSD of rainfall intensity is computed as follows:where indicates the radial wavenumber, is the number of pixels and represents rainfall intensity at individual pixels.
Figure 7 illustrates the results of the PSD calculations using Eq. 17 during the 0–3-h lead times. The PDE model can provide predictions without the loss of spatial information as it considers precipitation advection as a key factor in precipitation pattern evolution. The slight underestimation of the PSD at scales of 10–100 km in the PDE model is attributed to its inability to accurately predict the complex evolution processes observed in real precipitation patterns. For example, the assumptions of constant kinematic viscosity or diffusion coefficient may not adequately represent the intricate details of the actual environment. Additionally, although turbulent motions play a crucial role in precipitation pattern evolution, the motion tracked based on radar data in the two-dimensional space cannot account for the vertical component, which is a limitation. The U-Net model exhibits the blurriest prediction results among the models, whereas the blended model demonstrates less blurriness compared with the U-Net model. Moreover, the power spectrum evolution of the U-Net model shows a decrease in power density at small scales over time, whereas the power spectrum of the blended model maintains consistent power density at small scales, regardless of time. This indicates that the blended model can provide consistent forecasts from the perspective of spatial resolution during 0–3-h lead times.
FIGURE 7
To investigate the correlation between the predicted () and observed () values, the correlation coefficient, , was estimated using the equation described in :
Here, is the number of pixels; and for calculating , rainfall intensities and two motion components were considered (i.e., and were calculated). The closer this coefficient is to 1, the higher is the correlation between the predicted and ground truth values, when this value is close to 0, there is almost no correlation.
The evolution of the correlation coefficients of rainfall intensities, , obtained by Eq. 18, at 0–3-h lead times is shown in Figure 8A. Up to a 1.5-h lead time, U-Net exhibits the highest correlation; however subsequently, the correlation coefficient sharply decreases. This decline indicates that the U-Net model produces overly blurry predictions, resulting in poor correlation with the distribution of ground truth values. Notably, of the U-Net model measured at a 3-h lead time is almost comparable with that of Persistence. Conversely, the blended model shows decent performance at 0–3-h lead times, demonstrating the best performance after a 1.5-h lead time. This is attributed to the preservation of spatial resolution and influence of the PDE model, which accounts for the movement of precipitation patterns. Additionally, the blended model captures complex variations in precipitation intensity owing to the U-Net model, leading to a better correlation than the PDE model.
FIGURE 8
The correlation coefficients of two motion components, and , are displayed in Figures 8B, C. We suggest that compared with those of the U-Net model, the motion fields obtained by the blended model exhibit a better correlation with those obtained by the ground truth. As evident from the PSD analysis shown in Figure 7, this is because the blended model produces less blurry predictions compared with the U-Net model.
To evaluate the accuracy of precipitation prediction, performance was evaluated by measuring critical success index (CSI), probability of detection (POD), and false alarm ratio (FAR). These metrics are defined based on the values of hits (), misses (), and false alarms () as follows:
The CSI is used to evaluate the accuracy of precipitation prediction (Eq. 19). POD represents the ratio of correctly predicted events among the actual occurrences of precipitation (Eq. 20). A low POD value indicates a higher number of missed events. Higher values (closer to 1) for both CSI and POD indicate more accurate predictions. FAR is used to measure the ratio of incorrectly predicted events out of the total predicted events (Eq. 21). A higher FAR value indicates a larger number of false predictions, often associated with overestimating precipitation areas. A lower FAR value, closer to 0, indicates better predictions.
We evaluated the performance of precipitation events with intensities of or more, or more, and or more. Figure 9 shows the performance evaluation at 0–3-h lead times. In terms of the CSI, the PDE model generally underperforms compared with the other models. This is because it lacks the ability to predict the dynamic processes involved in precipitation evolution. Consistently, the POD values of the PDE model are lower compared with those of the other models (i.e., the number of missed events is higher than those of other models), and larger FAR values indicate a tendency to overestimate rainfall intensities and their areas. The blended model outperforms the other models up to a 3-h lead time, and this improvement is particularly significant when predicting moderate rainfall events (). The performance difference between the U-Net and blended models arises from the inability of the U-Net model to predict precipitation advection accurately owing to its blurry effects. These results are consistent with the results of qualitative analysis of Figures 4, 5.
FIGURE 9
The performance of the blended model was also assessed separately for each precipitation type (i.e., Central and Southern cases) to examine the characteristics of nowcasting outputs produced by the blended model (see Figure 10). Overall, the prediction accuracy for the Central case is higher than that for the Southern case. According to FAR scores, the models generally overestimate when predicting heavy rainfall in the Southern area. Additionally, lower POD values for the Southern case at 1–3-h lead times indicate a substantial number of missing events. We interpret that such performance difference is mainly related to the failure to detect precipitation growth in the southern coastal area. While the blended model accurately predicts the eastward movement of precipitation fields in both the central and southern regions of the Korean Peninsula, dynamic processes such as the growth and decay of precipitation are more prominent in the southern coastal regions. To predict such convective rainfall accurately, additional information such as satellite images would be required, which is beyond the scope of this paper.
FIGURE 10
Comparison between recent studies regarding precipitation nowcasting is presented in Table 1. Notably, direct comparisons cannot be conducted owing to the different detailed settings adopted in each study. Despite such limitations, indirect comparison could be conducted using baselines, such as Persistence and U-Net. In predictions of both moderate () and heavy rainfall (), the blended model demonstrates better performance than the other models proposed in previous studies at 0–3-h lead times. Although most previous studies have improved the prediction accuracy through the design and enhancement of deep learning models, we highlight the potential for achieving predictions at higher accuracies by incorporating efficient data processing methods from existing models.
TABLE 1
| Model | Area | Rain rates | Baseline | Lead time | Improvement over the baseline |
|---|---|---|---|---|---|
| U-Net | Germany | Persistence | 1 h | ||
| 1 h | |||||
| DGM | UK, US | U-Net | 1 h | ||
| DeepRaNE | South Korea | Persistence | 1 h | ||
| 3 h | |||||
| Persistence | 1 h | ||||
| 3 h | |||||
| This work (Blended Model) | South Korea | Persistence | 1 h | ||
| 3 h | |||||
| U-Neta | 1 h | ||||
| 3 h | |||||
| Persistence | 1 h | ||||
| 3 h | |||||
| U-Neta | 1 h | ||||
| 3 h |
Comparison between recent studies of precipitation nowcasting. The CSI scores were used for calculating the improvement over the baseline.
U-Net indicates the U-Net model employed in this study.
To test the operational feasibility of the blended model, we examined the performance of blended model in predicting precipitation in various seasons, each characterized by distinct weather patterns. We evaluated the monthly-averaged forecast skills based on the CSI scores for different seasons. Figure 11 shows the monthly-averaged performances of the blended and U-Net models during January (i.e., winter season), March (i.e., spring season), and September (i.e., late summer season) 2022, respectively. Notably, the prediction performance during September is the best, whereas it is the worst during January. This is mainly because during summer, more precipitation cases with sustained strong rainfall intensity substantially evolve by advection compared with those in winter or spring. The performances of both the blended and U-Net models follow such seasonal variability. Additionally, the blended model outperforms U-Net in predicting both moderate () and heavy () precipitation events. Hence, we interpret that the performance enhancement of the blended approach is robust for different seasons.
FIGURE 11
4 Summary
In this study, we proposed an approach to improve the performance of existing U-Net models by combining a PDE model that considers fluid dynamics with the U-Net architecture for precipitation prediction. The main objective was to address the blurriness in the predictions from traditional U-Net models. We conducted qualitative and quantitative performance analyses using precipitation cases in the Korean Peninsula. Based on the qualitative and quantitative evaluations, we found that improving the blurriness and enhancing the movement of precipitation patterns through fluid dynamics contribute to the improvement of the prediction accuracy at 0–3-h lead times. In particular, the performance enhancement measured by CSI, POD, and FAR remains robust across different thresholds of rainfall intensity and seasons. Moreover, the blended model demonstrates higher performance compared with PDE-based models that are currently used in operational short-term precipitation forecasting, including KMA’s nowcasting model based on MAPLE. Therefore, the blended model holds potential for future application in operational short-term forecasting. Additionally, the blended model can be easily applied as an extension of existing operational PDE models, providing further advantages.
While this study focused on the strong rainfall events last more than a few hours in the central and southern regions of Korean Peninsula, the prediction accuracy for predicting the isolated precipitation patterns, where rainfall generates or dissipates rapidly in localized areas, is also crucial. Previous deep learning-based studies have also demonstrated relatively poor performance in predicting such isolated precipitation compared with that in predicting precipitation in central and southern cases, which encompass relatively broader regions (; ). To further improve the performance of nowcasting models, enhancing the prediction capabilities for isolated precipitation patterns is crucial. As indicated by previous studies (e.g., ; ), training the data from satellites and surface observations associated with the physical processes of clouds and precipitation would significantly enhance the performance. Such an approach is beyond the scope of this paper, and will be part of future studies.
Statements
Data availability statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author. The radar data are freely available at the released website of the Korea Meteorological Administration (KMA) (https://data.kma.go.kr/cmmn/main.do).
Author contributions
J-HH: Conceptualization, Investigation, Methodology, Software, Supervision, Validation, Visualization, Writing–original draft. JP: Methodology, Writing–review and editing.
Funding
The author(s) declare financial support was received for the research, authorship, and/or publication of this article. This work was funded by the KMA Research and Development program “Developing AI technology for weather forecasting” under Grant (KMA 2021-00121).
Acknowledgments
The authors thank reviewers for their helpful and constructive comments on the manuscript.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1
AgrawalS.BarringtonL.BrombergC.BurgeJ.GazenC.HickeyJ. (2019). Machine learning for precipitation nowcasting from radar images. Arxiv preprint. Available at: https://arxiv.org/abs/1912.12132.
2
AyzelG.HeistermannM.WinterrathT. (2019). Optical flow models as an open benchmark for radar-based precipitation nowcasting (rainymotion v0.1). Geosci. Mod. Dev.12, 1387–1402. 10.5194/gmd-12-1387-2019
3
AyzelG.SchefferT.HeistermannM. (2020). RainNet v1.0: a convolutional neural network for radar-based precipitation nowcasting. Geosci. Mod. Dev.13, 2631–2644. 10.5194/gmd-13-2631-2020
4
BakkayM. C.SerrurierM.BurdáV. K.DupuyR.Cabrera-GutiérrezN. C.ZamoM.et al (2022). Precipitation nowcasting using deep neural network. Arxiv preprint. Available at: https://arxiv.org/abs/2203.13263.
5
BenjaminS. G.DéVényiD.WeygandtS. S.BrundageK. J.BrownJ. M.GrellG. A.et al (2004). An hourly assimilation-forecast cycle: the RUC. Mon. Wea. Rev.132, 495–518. 10.1175/1520-0493(2004)132<0495:AHACTR>2.0.CO;2
6
ChenL.CaoY.MaL.ZhangJ. (2020). A deep learning-based methodology for precipitation nowcasting with radar. Earth Space Sci.7, e2019EA00812. 10.1029/2019EA000812
7
ChoiS.KimY. (2022). Rad-cGAN v1.0: radar-based precipitation nowcasting model with conditional generative adversarial networks for multiple dam domains. Geosci. Model Dev.15, 5967–5985. 10.5194/gmd-15-5967-2022
8
ChoiY.ChaK.BackM.ChoiH.JeonT. (2021). RAIN-F+: the data-driven precipitation prediction model for integrated weather observations. Remote Sens.13, 3627. 10.3390/rs13183627
9
GermannU.ZawadzkiI. (2002). Scale-dependence of the predictability of precipitation from continental radar images. Part I: description of the methodology. Mon. Wea. Rev.130, 2859–2873. 10.1175/1520-0493(2002)130<2859:SDOTPO>2.0.CO;2
10
GermannU.ZawadzkiI. (2004). Scale-dependence of the predictability of precipitation from continental radar images. Part II: probability forecasts. J. Appl. Meteor.43, 74–89. 10.1175/1520-0450(2004)043<0074:SDOTPO>2.0.CO;2
11
GermannU.ZawadzkiI.TurnerB. (2006). Predictability of precipitation from continental radar images. Part IV: limits to prediction. J. Atmos. Sci.63, 2092–2108. 10.1175/JAS3735.1
12
JoE.ParkC.SonS. W.RohJ. W.LeeG. W.LeeY. H. (2020). Classification of localized heavy rainfall events in South Korea. Asia-Pacific J. Atmos. Sci.56, 77–88. 10.1007/s13143-019-00128-7
13
KimY.HongS. (2022a). Very short-term prediction of weather radar-based rainfall distribution and intensity over the Korean peninsula using convolutional long short-term memory network. Asia-Pac. J. Atmos. Sci.58, 489–506. 10.1007/s13143-022-00269-2
14
KimY.HongS. (2022b). Very short-term rainfall prediction using ground radar observations and conditional generative adversarial networks. IEEE Trans. Geoscience Remote Sens.60, 1–8. 10.1109/TGRS.2021.3108812
15
KoJ.LeeK.HwangH.OhS. G.SonS. W.ShinK. (2022). Effective training strategies for deep-learning-based precipitation nowcasting and estimation. Comput. Geosci.161, 105072. 10.1016/j.cageo.2022.105072
16
Korea Meteorological Administration (KMA) (2021). Abnormal climate report 2020 (in Korean). Available at: http://www.climate.go.kr/home/bbs/view.php?code=93&bname=abnormal&vcode=6494.
17
KoxT.GerholdL.UlbrichU. (2015). Perception and use of uncertainty in severe weather warnings by emergency services in Germany. Atmos. Res.158-159, 292–301. 10.1016/j.atmosres.2014.02.024
18
KwonS.JungS. H.LeeG. (2015). Inter-comparison of radar rainfall rate using Constant Altitude Plan Position Indicator and hybrid surface rainfall maps. J. Hydrology531, 234–247. 10.1016/j.jhydrol.2015.08.063
19
LebedevV.IvashkinV.RudenkoI.GanshinA.MolchanovA.OvcharenkoS.et al (2019). “Precipitation nowcasting with satellite imagery,” in Proceedings of the 25th ACM SIGKDD International Conference on Knowledge Discovery & Data Mining, 2680–2688. 10.1145/3292500.3330762
20
LeeD.-K.HaJ.-C.KimJ. (2008). Application of the Sawyer-Eliassen equation to the interpretation of the synoptic-scale dynamics of a heavy rainfall case over East Asia. Asia-Pac. J. Atmos. Sci.44, 49–68.
21
LeeH. C.LeeY. H.HaJ.-C.ChangD.-E.BellonA.ZawadzkiI.et al (2010). McGill algorithm for precipitation nowcasting by Lagrangian extrapolation (MAPLE) applied to the southSouth Korean radar network. Part II: real-time verification for the summer season. Asia-Pacific J. Atmos. Sci.46, 383–391. 10.1007/s13143-010-1009-9
22
LyuG.JungS. H.Namk.Y.KwonS.LeeC. R.LeeG. (2015). Improvement of radar rainfall estimation using radar reflectivity data from the hybrid lowest elevation angles. J. Korean Earth Sci. Soc.36, 109–124. 10.5467/JKESS.2015.36.1.109
23
LyuG.JungS. H.OhY.ParkH. M.LeeG. (2017). Accuracy evaluation of composite hybrid surface rainfall (HSR) using KMA weather radar network. J. Korean Earth Sci. Soc.38, 496–510. 10.5467/JKESS.2017.38.7.496
24
MiaoK.WangW.HuR.ZhangL.ZhangY.WangX.et al (2020). Multimodal semisupervised deep graph learning for automatic precipitation nowcasting. Math. Probl. Eng.2020, 1–9. 10.1155/2020/4018042
25
OhS. G.ParkC.SonS. W.KoJ.ShinK.KimS.et al (2023). Evaluation of deep-learning-based very short-term rainfall forecasts in South Korea. Asia-Pacific J. Atmos. Sci.59, 239–255. 10.1007/s13143-022-00310-4
26
ParkC.SonS. W.KimJ.ChangE. C.KimJ. H.JoE.et al (2021). Diverse synoptic weather patterns of warm-season heavy rainfall events in South Korea. Mon. Wea. Rev.149, 3875–3893. 10.1175/MWR-D-20-0388.1
27
PulkkinenS.NeriniD.Perez HortalA. A.Velasco-ForeroC.SeedA.GermannU.et al (2019). Pysteps: an open-source Python library for probabilistic precipitation nowcasting (v1.0). Geosci. Mod. Dev.12, 4185–4219. 10.5194/gmd-12-4185-2019
28
RavuriS.LencK.WillsonM.KanginK.LamR.MirowskiP.et al (2021). Skilful precipitation nowcasting using deep generative models of radar. Nature597, 672–677. 10.1038/s41586-021-03854-z
29
RonnebergerO.FischerP.BroxT. (2015). “U-net: convolutional networks for biomedical image segmentation,” in International Conference on Medical Image Computing and Computer-Assisted Intervention, 234–241. 10.1007/978-3-319-24574-4_28
30
RyuS.LyuG.DoY.LeeG. (2020). Improved rainfall nowcasting using Burgers’ equation. J. Hydrology581, 124140. 10.1016/j.jhydrol.2019.124140
31
SeoK.-H.SonJ.-H.LeeJ.-Y.ParkH.-S. (2015). Northern East Asian monsoon precipitation revealed by airmass variability and its prediction. J. Clim.28, 6221–6233. 10.1175/JCLI-D-14-00526.1
32
ShiX.ChenZ.WangH.YeungD.-Y.WongW.-K.WooW.-C. (2015). Convolutional lstm network: a machine learning approach for precipitation nowcasting. Adv. Neural Inf. Process. Syst.28.
33
ShiX.GaoZ.LausenL.WangH.YeungD.-Y. (2017). Deep learning for precipitation nowcasting: a benchmark and A new model. Adv. Neural Inf. Process. Syst.30.
34
ShresthaD.RobertsonD.WangQ.PaganoT.HapuarachchiH. (2013). Evaluation of numerical weather prediction model precipitation forecasts for short-term streamflow forecasting purpose. Hydrology Earth Syst. Sci.17, 1913–1931. 10.5194/hess-17-1913-2013
35
SivleA. D.AgerstenS.SchmidR.SimonA. (2022). Use and perception of weather forecast information across Europe. Meteorol. Appl.29, e2053. 10.1002/met.2053
36
SønderbyC. K.EspeholtL.HeekJ.DehghaniM.OliverA.SalimansT.et al (2020). Metnet: a neural weather model for precipitation forecasting. arXiv preprint. 10.48550/arXiv.2003.12140
37
SunJ.XueM.WilsonJ. W.ZawadzkiI.BallardS. P.Onvlee-HooimeyerJ.et al (2014). Use of NWP for nowcasting convective precipitation: recent progress and challenges. Bull. Am. Meteorological Soc.95, 409–426. 10.1175/BAMS-D-11-00263.1
38
ThorndablS.EinfaltT.WillemsP.NielsenJ. E.VeldhuisM.-C.Arnbjerg-NielsenK.et al (2017). Weather radar rainfall data in urban hydrology. Hydrol. Earth Syst. Sci.21, 1359–1380. 10.5194/hess-21-1359-2017
39
TurnerB. J.ZawadzkiI.GermannU. (2004). Predictability of precipitation from continental radar images. Part III: operational nowcasting implementation (MAPLE). J. Appl. Meteor.5, 231–248. 10.1175/1520-0450(2004)043<0231:POPFCR>2.0.CO;2
40
WangG.YangJ.WangD.LiuL. (2016). A quantitative comparison of precipitation forecasts between the storm-scale numerical weather prediction model and auto-nowcast system in Jiangsu, China. Atmos. Res.181, 1–11. 10.1016/j.atmosres.2016.06.004
41
WedelA.PockT.ZachC.BischofH.CremersD. (2009). An improved algorithm for TV-L1 optical flow. Stat. geometrical approaches Vis. motion analysis, 23–45. 10.1007/978-3-642-03061-1_2
42
WendroffB. (1968). Difference methods for initial-value problems (Robert D. Richtmyer and K. W. Morton). SIAM Rev.10, 381–383. 10.1137/1010073
43
YuW.YoonS.ChoiM.JungK. (2017). Performance comparison of rainfall and flood forecasts using short-term numerical weather prediction data from Korea and Japan. J. Korea Water Resour. Assoc.50, 537–549. 10.3741/JKWRA.2017.50.8.537
Summary
Keywords
advection-diffusion equation, deep learning, convolutional neural network, weather radar, precipitation nowcasting
Citation
Ha J-H and Park J (2024) Advancing very short-term rainfall prediction with blended U-Net and partial differential approaches. Front. Earth Sci. 11:1301523. doi: 10.3389/feart.2023.1301523
Received
25 September 2023
Accepted
27 December 2023
Published
11 January 2024
Volume
11 - 2023
Edited by
Yuh-Lang Lin, North Carolina Agricultural and Technical State University, United States
Reviewed by
Guy Jean-Pierre Schumann, University of Bristol, United Kingdom
Xia Wan, China Meteorological Administration, China
Updates
Copyright
© 2024 Ha and Park.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Ji-Hoon Ha, jhha223@korea.kr
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.