Next Article in Journal
Multiple-Band Electric Field Response to the Geomagnetic Storm on 4 November 2021
Previous Article in Journal
Variation in Vegetation Composition and Structure across Mudflat Areas in the Yellow River Delta, China
Previous Article in Special Issue
Combining Hydrological Models and Remote Sensing to Characterize Snowpack Dynamics in High Mountains
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Cumulative Rainfall Radar Recalibration with Rain Gauge Data Using the Colour Pattern Regression Algorithm QGIS Plugin

by
Pablo Blanco-Gómez
1,*,†,
Pau Estrany-Planas
2,† and
José Luis Jiménez-García
3,†
1
Vielca Ingenieros S.A., Avda. Aragón, 17, 46010 Valencia, Spain
2
Empresa Municipal de Agua y Alcantarillado (EMAYA), Son Pacs, Camí dels Reis, 400, 07010 Palma, Spain
3
Vielca Medio Ambiente S.L., Avda. Aragón, 17, 46010 Valencia, Spain
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Remote Sens. 2024, 16(18), 3496; https://doi.org/10.3390/rs16183496
Submission received: 25 July 2024 / Revised: 6 September 2024 / Accepted: 19 September 2024 / Published: 20 September 2024
(This article belongs to the Special Issue Remote Sensing of Floods: Progress, Challenges and Opportunities)

Abstract

:
Climate change is a major issue in wastewater management at local and regional levels, as it affects the frequency of flooding and therefore the need to update infrastructure and design regulations. To this end, rainfall data are the main input to hydraulic models used for the design of drainage systems and, in advanced contexts, for their real-time monitoring. Field observations are of great interest and water authorities are increasing the number of existing rain gauges, but at present they are scarce and require maintenance, so their number needs to be considered with their O&M costs. Remote sensors, including both the existing satellite rain products (SRPs) and radar imagery (RI), can complete the spatial distribution of rainfall and optimise the cost of observations. While most SRPs are based on re-analysis and have a lag in availability, RI can be obtained in near real time and is becoming increasingly popular in weather forecasting applications. Unfortunately, actual rainfall forecasts from RI observations are not accurate enough for real-time monitoring of drainage systems. In this paper, the Colour Pattern Regression (CPR) algorithm is used to recalibrate the 6 h rainfall values from RI provided by the Agencia Estatal de Meteorología (AEMET) with the observed rain gauge data, using as a case study the metropolitan area of Palma (Spain).

1. Introduction

The specific and local outcomes of climate change are uncertain and lead to impacts on the hydrological cycle by altering the frequency, intensity, spatial extent, or duration of weather and climate extremes, such as heavy precipitation events or droughts [1]. However, in the context of global warming and climate change, it is generally accepted that both temperature and precipitation patterns will change on a planetary scale [2] and, in this sense, the latest report of the Intergovernmental Panel on Climate Change (IPCC) [3] states that the frequency and intensity of heavy precipitation events have increased since the 1950s over most land areas for which observational data are sufficient for trend analysis, and that anthropogenic climate change is likely to be the main driver. An increase in the intensity of rainfall in urban contexts, together with an increase in urbanisation activities, leads to the need to increase the capacity of urban drainage systems.
Rainfall is characterised by a high spatio-temporal variability, which in turn complicates its quantitative description [4]. On one hand, rain gauges can be considered as a very precise measurement estimation of the actual rainfall, but cannot reproduce its spatial distribution [5]; and, on the other hand, meteorological radar presents a good spatial distribution of the rainfall, but a poorer quantitative estimation.
The combination of meteorological radar and distributed hydrological models results in a fundamental tool for flood risk management [4], but also the application of rainfall prediction techniques, or nowcasting, based on weather radar is ideal for this purpose due to its good spatial coverage (around a km2) and temporal resolution (less than 15 min) [6,7,8]. The benefit of using radar-based nowcasts as inputs to hydrological models has also been shown in urban catchments and flood forecasting [8,9,10]. Regardless, urban catchments are usually smaller than rural watersheds, the runoff transformation of the catchments is higher, and their reactive capacity is lower because natural retention is much less, so the immediacy and accuracy of rainfall information is of paramount interest.
However, the use of meteorological radar for operational applications in urban hydrology has been rare so far, either because of the more complex error characterisation of radar data compared to direct observations in rain gauges, or because of the actual expansion and development of the technical infrastructure—i.e., radar networks are not available everywhere. In addition, it is also a feature of radar measurements that the uncertainty of the measurement increases with higher rain intensities, which is of primary interest in urban hydrology, especially in high-intensity contexts [7]. Both measurement devices—i.e., meteorological radar and rain gauges—are complementary and thus the simultaneous use of both can provide the best spatial estimate of rainfall as required for urban hydrological applications.
This paper applies the Colour Pattern Regression (CPR) algorithm [11] to re-calibrate the 6 h rainfall images coming from the meteorological radar provided by [12]. In order to verify the suitability of the methodology, a case study was selected in Palma de Mallorca (Spain), where rainfall patterns and the capacity of the sewerage system affect the number of discharges into the Mediterranean Sea, which in turn affects water quality and the possibility of beach closures, ultimately affecting the tourism sector [13]. Only in the period 2015–2021 is there an average of 45 days per year when more than one discharge body is flowing—see Figure 1. In this context, the present study aims not only to quantify the spatio-temporal distribution of rainfall, coupling the spatial distribution of the meteorological radar with the observed rainfall of the rain gauges, while questioning the quantification of the radar, but also to improve the capabilities of a prospective digital twin of the sewerage system, facilitating the future infrastructure planning and maintenance of the EMAYA system.

2. Materials and Methods

2.1. Case Study

The city of Palma is located in the south of the island of Mallorca and is the capital of the Balearic Islands region of Spain, with an area of around 210 km2 and elevations ranging from 0 m, on the coast, to around 400 m in the municipality, but with heights of over 1000 m in the Tramuntana mountains, around 20 km north of the city centre (Figure 2).
Since 1943, the Empresa Municipal de Agua y Alcantarillado, S.A. (EMAYA) has been responsible for both water and waste management in the city of Palma, which has a population of 420,000 [14]. The tributary basins of the city have a Mediterranean climate, with a high seasonality of rainfall, experiencing bigger registers in autumn, which include records of more than 100 mm/day [15]. The rainwater is collected in either the storm or mixed sewer system, which ends up in the Palma II wastewater treatment plant, and both the treatment discharges and the overflows from the collection system end up in the Mediterranean Sea. Likewise, tourist activities are also affected by these discharges, as they affect the quality of bathing water and cause beach closures, so minimising the occurrence of such events is a priority for the EMAYA. In this sense, the recent update of the national regulation by means of Royal Decree 665/2023 [16] and the proposal for the adaptation of the Urban Wastewater Treatment Directive [17] establish criteria for the design of infrastructures that limit the discharge of wastewater into the sea.
Annual rainfall in this area ranges from 201.3 to 702.2 mm/year, with an annual average of 449 mm over the period 1980–2010 [18]. The rainy season from September to November accounts for 52.2% of the annual rainfall, while the dry season in summer (from June to August) accounts for 10% [19].

2.2. Conceptual Model

The aim of this study is to calibrate the meteorological radar with observed precipitation at rain gauging stations in the context of the Palma city drainage system. The 6 h cumulative rainfall product from the meteorological radar of [12] was used to verify the agreement of the Marshall–Palmer relationship between the observed reflectivity and the calculated rainfall and the observed rainfall in the existing rain gauges. The methodology used in this study is shown in Figure 3. It includes the following steps: (1) data homogenisation and quality control using the RainFA software [13]; (2) identification of 73 daily events with more than 10 mm accumulated rainfall in a single gauging station used for both spatial interpolation (3) and meteorological radar data acquisition (5); (3) spatial interpolation of the rainfall data at the gauging stations using the Inverse Distance Weighting (IDW) method; (4) as a result of the spatial interpolation (3), raster maps of observed precipitation were obtained and used for GIS calculations in QGIS, using the CPR algorithm plugin (7); (5) 6 h accumulated rainfall product from the meteorological radar [12] data acquisition and geo-referencing using the European Space Agency (ESA) SeNtinel Applications Platform (SNAP) software; (6) transforming the geo-referenced meteorological radar maps into a 3-band RGB image using the colour legend provided by [12]; (7) applying the Colour Pattern Regression (CPR) algorithm [11] to obtain simulated raster maps of precipitation, minimising the differences between the spatially distributed IDW observed precipitation from (4) and the colour ramp from the radar images computed in (6)—converting the 1-band product provided by [12] to 3-band RGB images; (7) this procedure is generalised to recalibrate satellite precipitation meteorological radar products for Palma de Mallorca (Spain) in the period 2015–2021.

2.3. Rainfall Stations Data

Precipitation data were collected from 11 different rain gauges operated by AEMET—i.e., B228, B236, and B278—and the non-profit regional meteorological association Balears Meteo [20]—i.e., El Pil·lari, La Bonanova, Pòrtol, Son Rapinya, Pont d’Inca, Puntiró, Secar de la Real, and Son Ferriol. The data sets were collected from 2015 to 2021 and are recorded at different frequencies—i.e., 1, 5, and 10 min. A data quality control software called RainFA [21] was used to detect outliers and temporal homogenisation. First, the various time series were adjusted to a 10 min time step, summing 2 or 10 registers where appropriate, and then four different data quality controls were applied: (1) user-defined thresholds, using the standard procedure of converting to Nan the negative numbers and the records greater than 40 mm in 10 min; (2) box and whisker and bivariate Highest Density Region (HDR) boxplots, using the summary statistics of the distribution to obtain the quartiles and the potential outliers, comparing stations in pairs with the HDR boxplots and looking for the correspondence of common values densities; (3) double-mass curves, comparing the cumulative rainfall of a single station with the average records of all stations and looking for trend anomalies; and (4) trend analysis, looking for non-stationary series with the Mann–Kendall test. As a result, the Son Ferriol station data set was removed from the study as it had less than two years of data and its registers did not agree with the rest of the stations in the bivariate HDR boxplots and double-mass curve analysis. In addition, 4 discrete values were discarded from the rest of the time series.
Following this data regularisation, daily precipitation events were discretised and those with cumulative registers greater than 10 mm were selected for calibration of the radar imagery. Storm dates were used to request meteorological radar data from AEMET and, where available, to accumulate precipitation into four daily 6 h values in accordance with the radar products described in Section 2.4.
The 6 h accumulated rainfall was represented as a point-vector shapefile in QGIS and spatially distributed using the Inverse Distance Weighting (IDW) method to produce raster maps of the observed rainfall, which were compared with the different radar images using the CPR algorithm QGIS plugin.

2.4. Meteorological Radar

The meteorological data come from the 6 h accumulated precipitation product provided by the Spanish Meteorological Service (AEMET, [12]). The instrument is equipped with Doppler capability and operates in the 5.6 GHz band [22]. It currently provides raster images for a circle of 240 km radius (long-range mode) reprojected to EPSG:4326 (World Geodetic System 1984, WGS84) centred on the Illes Balears radar (39.57ºN, 2.65ºE). Each image has 494 × 536 pixels and a spatial resolution of 1 × 1 km2 in the long-range mode.
The 6 h accumulated precipitation product in mm is the sum of the hourly precipitation products of the period, calculated using the Marshall–Palmer precipitation relation [23], shown in Equation (1).
Z = a · R b ,
where Z is the reflectivity factor, R is the rain rate (in mm/h), and a and b are constant variables that depend on the type of precipitation. According to [4,5,7,12], for rainfall precipitation they take the following values: a = 200 and b = 1.6.
However, the same source [12] states that there is no clear relationship between reflectivity and rain water content, nor between reflectivity and rain rate, as the latter is also affected by the velocity at which rain falls. They conclude by suggesting the use of reflectivity raster maps—as detected by radar—instead of rain rates.
Ref. [24] studied the variation of these coefficients with the raindrop size distribution associated with different types of rainfall—e.g., ‘orographic’, ‘thunderstorm’, ‘stratiform’, or ‘showers’. They concluded that (1) for a given rainfall rate, ‘orographic’ rainfall presents smaller raindrop sizes and larger concentrations than ‘thunderstorms’, as would be expected for these types of rainfall, and this explains the smaller prefactors—i.e., a coefficients—and larger exponents—i.e., b coefficients—of the exponential Z-R relationships of the ‘orographic’ rainfall when compared to ‘thunderstorms’; (2) it has been difficult to obtain clear conclusions for the other types of rainfall; and (3) the coefficients adopted in [12] are used in many parts of the world because they are close to the mean power law of different authors [24,25,26] for rain rates between 1 and 50 mm/h.
The 6 h radar data acquired was manually geo-referenced using the tool Ground Control Points (GCP) Manager of the ESA SNAP [27] software and then transformed into a 3-banded RGB raster image by means of a colour classification algorithm using the rasterio [28] Python library.

2.5. CPR Algorithm

The CPR algorithm is an add-on module to the open-source QGIS software implemented as a raster interpolation method to determine the relationship between 3-band colour raster products and raster maps [11]. It was developed to fill the gap of establishing correlations between raster maps and aerial colour patterns, using a linear regression between the observed (O) values of the raster maps and the RGB bands of the aerial images—referred to as simulated, S, values.
After (1) normalizing the pixel sizes and (2) clipping with a mask, both layers can be treated as matrices of size m × n and the objective function can be written as follows—Equation (2):
m i n j = 1 n i = 1 m O i j S i j 2 ,
where the observed values (O) are represented in Equation (3) and the simulated values (S) are represented as a matrix polynomial equation—i.e., Equation (4)—using three parameters—i.e., R, G, and B—which are minimised in the calculation process.
O 11 O 1 n O m 1 O m n ,
S 11 S 1 n S m 1 S m n = R · R 11 R 1 n R m 1 R m n + G · G 11 G 1 n G m 1 G m n + B · B 11 B 1 n B m 1 B m n ,
The QGIS CPR algorithm plugin also implements three goodness-of-fit metrics, (1) the Normalised Nash–Sutcliffe Efficiency Coefficient (NNSE), (2) Kling–Gupta Efficiency (KGE), and (3) the Percent Bias Index (PBIAS), which ultimately establish four categories of performance metrics: Very good, Good, Satisfactory, and Unsatisfactory. In the present case, the pixel size of both products—i.e., the IDW-interpolated rainfall raster map and the 3-band coloured radar image—is the same, so there is no spatially distributed metric, but a single performance number of each date–time evaluation.
Refs. [29,30], among others, defined values of the goodness-of-fit metrics for hydrological modelling of streamflow and suspended sediment for different time steps, but there is no reference of the corresponding values for rainfall data from radar observations at the 6 h time step. However, the QGIS CPR algorithm allows the user to adjust the thresholds to suit the needs of the problem.

2.6. Evaluation Procedure

In line with the goodness-of-fit evaluation measurements of the CPR algorithm, three quantitative performance metrics were used to measure, in this case, the accuracy of the meteorological radar predictions with respect to the observed raster maps. The statistics used are (see Table 1): (1) root mean square error (RMSE); (2) the Nash–Sutcliffe efficiency coefficient (NSE); and (3) percent bias (PBIAS). The RMSE explains the deviation between two sets of values, which in this case correspond to the difference between the inferred value from the colour pattern of the CPR process—i.e., observed data—and the corresponding value that was assigned in the colour legend of the meteorological radar product—i.e., simulated data. NSE indicates how well the correlation of observed and simulated data fits the 1:1 line when plotted together, while PBIAS indicates the degree of over- or under-estimation of the observed value as a percentage.
Where O i and S i are, respectively, observed and simulated data, i is the data index and n is the total number of measurements, and O ¯ is the average value of observed data.

3. Results

As mentioned in Section 2.3, rainfall data were curated and homogenised to a 10 min time step using the RainFA software [21]. Daily rainfall values were then detected and those with a cumulative rainfall greater than 10 mm and their preceding and following days were used to obtain the meteorological radar products from [12]. Up to four daily meteorological radar images—i.e., 6 h radar images at 00:00, 06:00, 12:00, and 18:00—were obtained from 73 different dates, resulting in the following findings: (1) there are 54 days with a gauged rainfall greater than 10 mm and available radar inputs; (2) there are 19 images (out of 161) with observed rainfall—i.e., records greater than 0.0 mm and up to 17.8 mm—where the available meteorological radar image was not able to detect rainfall in the area of interest; and (3) during these days, a total of 21 (out of 42) meteorological radar images represent rainfall while there are no records of observed rainfall at the gauging stations. This means that in the case study and in the situation of moderate to extreme rainfall [31], the radar products are able to detect 88.2% of the 6 h time slots with rainfall, but the radar also detects rainfall in the 50.0% of cases where there is no register at the gauging stations.
The rainfall records were then accumulated in accordance with the meteorological radar images and spatially distributed using the IDW interpolation method. The monthly distribution of the selected storm events is shown in Figure 4 and is in line with the values referred to in [19] and presented in Section 2.1, indicating that the most frequent period of moderate to extreme rainfall in the 2015–2021 period is September to November, accounting for 49.3% of the total events, while the period from June to August accounts for only 9.6%. Although the percentages are very similar, the reference values—i.e., annual rainfall vs. number of moderate to extreme rainfall events—are completely different, suggesting a direct relationship between monthly rainfall and the distribution of storms throughout the year, which should be studied in detail.
As mentioned in the Section 2.4, the radar images were manually georeferenced using the GCP Manager tool of the ESA SNAP software [27] and then transformed into a three-band RGB raster image using two different colour label dictionaries for dates before and after 1 January 2016. The following Table 2 shows the three-band RGB value dictionary for both periods.
The corresponding date–time pairs of (1) the three-band RGB meteorological radar image and (2) the 6 h cumulative rainfall IDW interpolated raster maps were used to apply the CPR algorithm to assess the validity of the meteorological radar cumulative rainfall values while improving the spatial distribution of the on-site weather stations. The results of the application of the CPR algorithm are presented in Figure 5 and Table A1, including the following: (a–b) date–time of the rainfall event; (c–e) R-G-B calculated parameters; (f–h) goodness-of-fit statistics—i.e., NNSE, KGE, and PBIAS—of the performance of the algorithm; and (i) maximum observed 6 h cumulative rainfall value.
The model parameters of the CPR algorithm take on coherent values with the colour ranges shown in Table 2, generally showing larger values with increasing rainfall and including Green shades for higher amounts of rainfall. In this sense, there is a correspondence between the maximum values of rainfall accumulated in rain gauges within 6 h slots and the colour thresholds identified in the meteorological radar legend mentioned above. Figure 6 shows the distribution values of rainfall associated with only Blue tones and including Green tones: (a) shows the dispersion of the results for both coefficients, showing a direct proportionality between their values and the maximum amount of rainfall registered, and, as expected, greater Green coefficient values than the Blue ones for higher rainfall and (b) highlights the Q1 and Q3 quartiles—which account for 50% of the values—in the graph, resulting in a remarkable similarity with the Blue and Green rainfall thresholds of the meteorological radar legend. However, despite the promising general results, a more detailed study is carried out in Section 4 for the different thresholds.
Regarding the CPR algorithm performance statistics, using the KGE values proposed in [11] as a reference—i.e., values between 0 and 1 indicate Very good performance, between −0.5 and 0 are considered Good, and values greater than −2 are considered Satisfactory—137 out of 142 radar images that were able to detect precipitation have at least a Satisfactory performance that explains the spatial distribution of precipitation according to its colour pattern, 106 have a Good performance, while only 28 have a Very Good performance. In terms of NNSE, the thresholds proposed in [11] cannot be considered as a reference, in line with the final conclusion and future work of the referenced publication. Consistent with the KGE statistic performance ranking, the variation in NNSE values suggests the use of the following thresholds: (1) NNSE > 0.40 for Very good performance, (2) 0.40 > NNSE > 0.30 for Good performance, and (3) 0.30 > NNSE > 0.20 for Satisfactory performance, resulting in NNSE < 0.20 as Unsatisfactory models, see Figure 7a. Finally, PBIAS is positive in all cases, indicating that the IDW spatially distributed gauged values are overestimating the simulated values from the colour patterns of the meteorological radar. In addition, using the same KGE classification, PBIAS results are also in agreement with performance ratings established in [11], concentrating both Very good and Good performance within the range of [11.5; 35.2]% and obtaining an average of 51.4% for Satisfactory models, as presented in Figure 7b.

4. Discussion

As seen in the previous Section 3, the results of the CPR algorithm in terms of the KGE statistic are at least Satisfactory in practically all cases, and Good in most of them—i.e., 96.5% and 74.6%, respectively. This means that the colour patterns of the meteorological radar products are a good indicator of the rainfall observed in the rain gauges. However, this does not mean that there is a unique relationship between radar observations (which ultimately quantify reflectivity [7,12]), referred to as Z in Equation (1), and rain rate, referred to as R. In this sense, Einfalt et al. [7], among others, highlight the difficulties in the transformation between Z and R and the considerable deviations in their results [7,32,33]. In particular, by applying different equations [23,34,35], can be seen that the most important differences exist for higher rainfall intensities [7], but in the present case, data on these ranges are missing and so it is only possible to assess the results for rainfall intensities below 32 mm.
In order to assess the ability of the meteorological radar to match the observations of rain gauges, the following procedure was applied for every date–time: (1) create a mask of the different colour categories of the meteorological radar; (2) apply the masks to the IDW-interpolated rainfall raster map and calculate its basic statistics—i.e., minimum, maximum, and average; (3) construct an observed array with the increasing colour categories’ minimum and maximum values; (4) compare the observed array with the simulated array, generated from the meteorological radar legend as [0.5, 1; 1, 2; 2, 4; 4, 8; 8, 16; 16, 32; 32, 64; 64, 128; 128, 256]; and (5) evaluate their goodness of fit using the performance metrics described in Section 2.6. First, an evaluation was performed using the RMSE, discarding the date–time registers that exceeded the maximum observed rainfall at the rain gauges, and found that 114 out of 142 (i.e., 80.3%) did not exceed this condition. Table A2 shows the corresponding performance statistics for the selected dates.
Using the recommended values in [29] for the NSE statistic, only 28 of the models have a performance rating of Very good and a total of 44 have a performance rating of at least Good, corresponding to 19.7% and 31.0%, respectively. In addition, most of the models (i.e., 39 out of 44) have a positive PBIAS, meaning that the observed values in the rain gauge exceed the simulated values in the meteorological radar legend. Even for at least Satisfactory models, 63 out of 142 (only 44.4%), PBIAS is positive in 55 cases—i.e., 87.3%. Such results invite reconsideration of the Marshall–Palmer relation (Equation (1)) coefficients for the transformation of radar reflectivity into rain rate (Z-R), as it is clearly underestimating the rainfall results.
The minimum and maximum values were then presented in ranges and their evolution studied for the Satisfactory, Good, and Very good performance models, in order to identify appropriate thresholds that take into account the rainfall variability in the context of Palma de Mallorca. Only the ranges with a minimum number of registers greater than five were considered, and so corresponded to [0.5, 1]; [1, 2]; [2, 4]; [4, 8]; [8, 16]; [16, 32]. Figure 8 shows the distribution boxplots of the corresponding combinations.
As a result of the analysis, the following thresholds were defined: [0.5, 2.5]; [2.5, 4.5]; [4.5, 6.5]; [6.5, 12]; [12, 16]; [16, 28]. Comparing them with the original ones used by AEMET [12], it is noticeable that below 16 mm of accumulated rainfall, the simulated values used in the meteorological radar legend underestimate the observed rainfall in the rain gauges. This result is consistent with the PBIAS statistics result mentioned above.
Next, the evaluation procedure to assess the ability of the meteorological radar to match the observations of rain gauges, was adapted to assess its capacity of being matched by the CPR-interpolated raster map. In this case, the procedure was repeated from step (2) and the masks were applied to the CPR-interpolated raster map instead of being applied to the IDW-interpolated rainfall raster map, and, as the map was constructed from the colour patterns, a single value corresponded to every colour instead of a series of them. The observed array was compared with the same simulated array (4) and the performance metrics were calculated (5).
The results were then evaluated for the at least Satisfactory models in NSE of the previous evaluation procedure—i.e., the IDW-interpolated rainfall raster map—and showed that 42 out of 63 of the cases (66.7%) could also be considered as Satisfactory, 30 out of 44 (68.2%) could also be considered as Good, and 17 out of 28 (60.7%) could also be considered as Very good. Using the same distribution evaluation as for the IDW-interpolated raster map—i.e., see Figure 8—the defined thresholds were also appropriate, thus validating both the new proposed ranges and the use of the CPR algorithm to recalibrate the meteorological radar estimates.
Finally, the Marshall–Palmer precipitation relationship (Equation (1)) was recalculated for the new rainfall thresholds, first changing the units of the reflectivity factor from Z (in mm6/m3) to Z (in dBZ) using the following relationship, Z ( d B Z ) = 10 · l o g ( Z ) [4,12], and then minimising the distance between the calculated values of Z (dBZ) and the reference thresholds of the equivalent 6 h meteorological radar product in [12]—i.e., [12, 18]; [18, 24]; [24, 30]; [30, 36]; [36, 42]; [42, 48]. The coefficient a was fixed at 200 after warming up the results and comparing them with the bibliographical references [7,32,33,34,35], and the b calculation resulted in 1.335. Figure 9 shows the resulting Z-R curve compared with three other existing curves: (a) the Marshall–Palmer relationship, Z = 200 · R 1.6 ; (b) the tropical climate relationship, Z = 250 · R 1.2 ; and (c) the convective storm relationship, Z = 300 · R 1.4 .
For the period 2015–2021 and the analysis of moderate to extreme events—i.e., rainfall rates up to 28 mm/h—the Z-R relationship is very similar to that for tropical climates and significantly different from the standard Marshall–Palmer coefficients.

5. Conclusions

The EMAYA is developing a digital twin of the urban drainage system of the city of Palma in Mallorca (Spain), whose main input data for being operational is rainfall, both for nowcasting and forecasting solutions. On-site observations with rain gauges are of utmost interest due to their estimation precision, but their network needs to be increased in density to improve their spatial resolution, thus increasing the O&M costs. A trade-off between cost and precision is needed and remote sensors, such as satellite observation products, represent an alternative to optimise the cost of the rain gauge network, but their measurements should be calibrated to improve their precision. The present manuscript uses the CPR algorithm to recalibrate the meteorological radar information and validates its results by calculating the new thresholds from the IDW-interpolated rainfall raster map. The meteorological radar has demonstrated its ability to detect rainfall, with up to 88.2% of the 6 h time slots containing rainfall. Results from the application of the CPR algorithm suggest that the colour patterns of the meteorological radar products are a good indicator of the rainfall observed at the rain gauges. New thresholds for the meteorological radar colour patterns were calculated from the IDW-interpolated raster map. When the new ranges were applied to the CPR-calculated raster map, the same values were obtained, validating the use of the CPR algorithm to recalibrate the meteorological radar estimates.
Finally, in terms of methodology, the CPR algorithm should be adapted to the QGIS plugin if it is to be used for this particular purpose, with the following performance rating values for NNSE: Very good > 0.40 > Good > 0.30 > Satisfactory > 0.20 > Unsatisfactory, in accordance with the reference paper [11].

Author Contributions

P.B.-G.: conceptualisation; methodology; software; validation; formal analysis; investigation; resources; data curation; writing—original draft; writing—review and editing; project administration. P.E.-P.: methodology; investigation; resources; data curation; writing—original draft; visualisation; supervision. J.L.J.-G.: software; resources; data curation. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.

Acknowledgments

The authors are grateful to the EMAYA for the provision of precipitation data from both BalearsMeteo and AEMET. In addition, we would like to thank Vicente M. Candela Canales for supporting the R&D investments and programmes within the Vielca companies.

Conflicts of Interest

Author Pablo Blanco-Gómez was employed by the company Vielca Ingenieros S.A., author Pau Estrany-Planas was employed by the company Empresa Municipal de Agua y Alcantarillado (EMAYA) and author José Luis Jiménez-García was employed by the company Vielca Medio Ambiente S.L. All 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.

Appendix A

Table A1. CPR algorithm results including, for each date and time evaluated, the calculated R, G, and B correlation coefficients; the goodness-of-fit performance metrics NNSE, KGE, and PBIAS; and the maximum 6 h cumulative rainfall value observed at any of the existing gauging stations.
Table A1. CPR algorithm results including, for each date and time evaluated, the calculated R, G, and B correlation coefficients; the goodness-of-fit performance metrics NNSE, KGE, and PBIAS; and the maximum 6 h cumulative rainfall value observed at any of the existing gauging stations.
Date (dd.mm.yyyy)TimeRGBNNSEKGEPBIASMax. 6 h Rainfall (mm)
15.08.201500:000.000000.000000.00000 7.2
06:000.000000.072150.045430.2356−0.095517.3416.0
12:000.000000.130400.101930.2358−0.423917.3727.3
18:000.000000.000000.00000 2.4
17.08.201500:000.000000.000000.00000 0.0
12:000.000000.066950.063690.1681−0.443511.0110.0
18:000.000000.000000.020480.1376−0.465821.996.4
03.09.201500:000.000000.000000.00000 0.0
06:000.000000.078270.059040.0821−1.019024.1816.8
18:000.000000.000000.00000 0.0
04.09.201500:000.000000.000000.00000 0.0
06:000.000000.000000.00000 0.0
12:000.000000.101800.158630.3607−0.218820.9159.4
30.09.201506:000.000000.174890.162320.0182−4.917810.5520.8
12:000.000000.094690.123090.0220−4.326914.9116.0
18:000.000000.083300.037430.23510.024622.8911.0
14.10.201500:000.000000.000000.00000 2.6
06:000.000000.084020.087000.2648−0.29948.5218.9
12:000.000000.083060.046680.2149−0.174912.4815.0
18:000.000000.000000.00000 6.6
10.03.201600:000.000000.000000.003360.2374−0.585599.451.4
06:000.001250.051750.080240.2333−0.044015.0514.2
12:000.000000.000000.002390.3860−0.186925.442.0
18:000.000000.022990.005950.3177−0.304817.783.8
01.04.201600:000.000000.000000.00000 0.0
06:000.000000.000000.000470.40220.241046.830.2
18:000.000000.109860.075590.2154−0.294714.3524.2
08.04.201600:000.000000.000000.00000 0.0
06:000.000000.000000.00000 0.0
12:000.000000.038000.035750.2859−0.07519.968.4
18:000.000000.107700.064890.2355−0.112210.8515.6
20.09.201600:000.000000.000000.00000 0.0
06:000.000000.000000.00000 0.0
18:000.000000.000000.014050.4421−0.05899.0410.4
23.09.201600:000.000000.000000.00000 0.0
06:000.000000.000000.00000 0.0
12:000.000000.000220.000570.4691−0.006761.730.4
18:000.000000.049470.041010.3805−0.129627.3727.6
24.09.201606:000.000000.000000.001790.2736−0.323891.180.8
12:000.000000.000000.00000 0.0
18:000.000000.045080.073790.0828−0.745393.8536.8
06.10.201600:000.000000.000000.00000 0.0
06:000.000000.000000.00000 0.0
18:000.000000.025680.022180.46510.165028.7416.8
20.10.201600:000.000000.390390.247410.44440.10878.1193.4
06:000.000000.000000.00000 0.4
12:000.000000.000000.008930.2981−0.052530.255.6
18:000.000000.000000.00000 1.0
22.10.201600:000.000000.000000.00000 0.0
06:000.000000.083690.061240.1820−0.329311.5113.2
12:000.000000.000000.009550.1035−0.369158.084.2
18:000.000000.005530.006080.42210.068565.813.2
14.11.201600:000.000000.000000.000190.2941−0.608996.530.2
06:000.000000.125950.104950.1063−0.845910.0717.4
12:000.000000.074720.074350.0793−1.221313.6812.2
18:000.000000.047100.019780.1239−0.572020.786.8
05.12.201600:000.000000.000000.00000 0.8
06:000.000000.142270.093370.2800−0.077410.9123.6
12:000.000000.000000.004070.1001−0.186780.881.8
18:000.000000.000000.00000 0.4
20.01.201700:000.000000.000280.000390.3807−0.101524.460.2
12:000.000000.118050.219760.0113−6.607619.8928.4
21.01.201700:000.000000.042550.014850.43280.049410.257.8
06:000.000000.092880.067630.37440.34015.8213.2
18:000.000000.108360.094540.2368−0.025911.3418.0
27.01.201700:000.000000.000110.000250.4884−0.185714.780.2
06:000.000000.000000.00000 0.0
12:000.000000.000000.00000 0.0
18:000.000000.078220.064950.32070.099721.4424.5
24.03.201700:000.000000.009290.004740.4543−0.568622.935.7
06:000.000000.000000.003160.2479−0.422068.343.8
12:000.000000.069800.054000.2533−0.34009.6812.0
18:000.000000.177290.146790.0423−2.67386.4721.4
05.06.201700:000.000000.000000.00000 0.2
06:000.000000.193730.131720.3190−0.219111.4038.4
12:000.000000.152690.243400.1852−0.255017.5151.0
18:000.000000.000000.009700.2031−0.666699.716.4
15.09.201700:000.000000.000000.00000 0.0
06:000.000000.000000.00000 0.0
12:000.000000.032760.008160.47570.212236.017.4
18:000.000000.238950.223100.1300−0.543018.2153.4
19.10.201700:000.000000.000000.00000 3.0
06:000.000000.131180.081030.0879−1.105811.5517.4
12:000.000000.105230.085160.3308−0.31319.5720.1
18:000.000000.000000.012240.1166−0.681599.799.2
06.02.201800:000.000000.000000.005920.3361−0.193324.242.6
06:000.000000.082400.056430.2170−0.204510.2510.2
12:000.000000.023530.012870.1926−0.388216.713.4
18:000.000000.039150.048330.2194−0.297016.4310.9
24.03.201800:000.000000.000000.00000 0.0
06:000.000000.000000.016170.3350−0.537615.758.6
12:000.000000.112710.087590.1488−0.462110.9617.8
18:000.000000.079260.019150.42870.200625.1017.4
14.04.201800:000.000000.000000.000490.4229−0.505786.380.8
06:000.000000.000000.003190.1126−0.518751.591.6
12:000.000000.000000.062760.1274−0.719812.3811.0
18:000.000000.000000.002330.29300.124267.581.2
01.05.201800:000.000000.000000.00000 0.0
06:000.000000.059530.030500.3111−0.146020.2013.6
12:000.000000.016650.006200.4207−0.236139.4711.3
18:000.000000.000000.018880.3354−0.297146.0720.0
13.05.201800:000.000000.000000.00000 0.8
06:000.000000.099600.068070.2157−0.173610.6519.2
12:000.000000.000000.000870.3859−0.739399.813.8
18:000.000000.000000.00000 0.0
10.09.201800:000.000000.000000.00000 0.0
06:000.000000.041520.027720.41840.012817.3916.4
12:000.000000.000000.071430.1713−0.248887.1033.0
18:000.000000.000000.00000 0.0
18.10.201800:000.000000.011730.005560.4698−0.135012.384.6
06:000.000000.000000.010460.3471−0.028058.976.0
12:000.000000.000000.008560.39370.080522.626.8
18:000.000000.042830.025740.3943−0.351316.3014.7
27.10.201800:000.000000.000000.00000 0.0
06:000.000000.000000.00000 0.2
12:000.000000.037700.046360.3741−0.123217.0113.2
18:000.000000.120690.100400.42760.04978.5242.0
28.10.201800:000.000000.106800.072880.2799−0.057510.1619.0
06:000.000000.000000.00000 17.8
12:000.000000.141570.049130.3249−0.270726.3424.8
18:000.000000.004650.009620.41300.337450.924.6
20.11.201800:000.000000.000000.00000 0.0
06:000.000000.000000.000640.2869−0.471996.260.4
12:000.000000.077630.012080.49940.367816.1112.0
18:000.000000.000000.000570.2448−0.107078.670.2
27.08.201900:000.000000.000000.00000 0.0
06:000.000000.000000.00000 0.0
12:000.000000.055190.034570.3435−0.217512.2114.5
18:000.000000.000000.030860.3576−0.016525.8126.0
21.10.201900:000.000000.024660.031560.4677−0.190211.8215.8
06:000.000000.000000.049100.27710.03329.7414.4
12:000.000000.000000.017860.2474−0.213210.995.1
18:000.000000.032620.013470.2389−0.001215.154.4
22.10.201900:000.000000.104690.228990.1012−0.799520.7934.8
06:000.000000.234230.205880.2094−0.23077.1534.2
12:000.000000.000000.001850.3723−0.743599.9411.2
18:000.000000.000000.000770.32720.016843.170.6
08.11.201900:000.000000.057820.024050.2880−0.206824.3512.4
06:000.000000.000000.005160.3753−0.235617.004.6
12:000.000000.000000.000290.3174−0.610498.560.2
18:000.000000.000000.000450.4147−0.725899.452.0
22.11.201900:000.000000.000000.00000 0.0
06:000.000000.000000.003390.22010.033245.391.4
12:000.000000.079730.042860.2222−0.192925.2216.2
18:000.000000.000000.002730.48510.176020.683.4
02.12.201900:000.000000.000000.00000 0.0
06:000.000000.062220.037720.3141−0.369224.4716.8
12:000.000000.000000.001070.27200.165168.150.4
18:000.000000.000000.00000 0.2
04.12.201900:000.000000.000000.00000 14.6
12:000.000000.104600.158490.2332−0.117917.4440.8
18:000.000000.000000.013960.2982−0.100912.565.2
19.01.202000:000.000000.000000.00000 0.0
06:000.000000.000000.003100.1425−0.367844.531.6
18:000.000000.158130.110340.35670.25415.8921.0
22.03.202000:000.000000.000930.000360.3250−0.903888.161.0
06:000.000000.000000.00000 0.0
12:000.000000.065160.053960.3103−0.400723.3721.7
18:000.000000.004390.004700.3867−0.308179.322.4
13.04.202000:000.000000.000000.00000 4.6
06:000.000000.000000.033920.2056−0.539713.4211.1
12:000.000000.000000.070100.2372−0.03976.4016.0
18:000.000000.000000.00000 0.0
20.04.202000:000.000000.000000.008670.44310.17176.804.3
06:000.000000.133230.044100.25660.011430.8017.2
12:000.000000.000000.010090.3816−0.337831.0913.7
18:000.000000.236450.138890.3362−0.129911.6831.4
29.08.20200:000.000000.000000.00000 0.0
06:000.000000.000000.00000 0.4
12:000.000000.045370.034120.3928−0.283721.9722.2
18:000.000000.000000.040430.0846−0.568799.0518.0
30.08.20210:000.000000.000000.00000 0.0
06:000.000000.000000.00000 0.0
12:000.000000.061660.079890.2982−0.446514.5422.2
18:000.000000.000000.00000 3.8
01.10.20210:000.000000.000000.00000 0.2
06:000.000000.000000.00000 0.0
12:000.000000.080980.022290.3358−0.675423.5114.6
18:000.000000.000000.017000.1630−0.303586.1210.0
30.10.20210:000.000000.000000.00000 0.0
06:000.000000.000000.011680.4004−0.408914.257.3
12:000.000000.148040.077110.2142−0.10489.7521.4
18:000.000000.000000.000920.2177−0.620395.641.0
10.11.20210:000.000000.000000.056400.0326−2.857826.8717.4
06:000.000000.129930.086960.1337−0.59159.6817.6
12:000.000000.108630.043350.2869−0.359715.7217.2
18:000.000000.059210.025280.1956−0.530015.947.4
11.11.20210:000.000000.019030.007400.3762−0.139214.803.6
06:000.000000.104920.026160.3295−0.430718.8417.0
12:000.000000.117180.034240.2605−0.626822.0216.3
18:000.000000.078820.029750.2020−0.614020.5510.4
23.11.20210:000.000000.006650.002850.2659−0.151734.971.5
06:000.000000.243260.082130.39290.029818.7836.5
12:000.000000.094040.087810.1373−0.559410.6916.4
18:000.000000.158010.096980.27430.00399.0024.0
25.11.20210:000.000000.000000.00000 0.0
06:000.000000.000000.000230.4554−0.407291.030.2
12:000.000000.084430.038640.2621−0.362515.0611.5
18:000.000000.000000.00000 0.0
28.11.20210:000.000000.000000.003640.2483−0.020843.752.6
06:000.000000.000000.005400.28600.064660.553.0
12:000.000000.000000.005580.47260.243244.825.2
18:000.000000.000000.00000 0.0

Appendix B

Table A2. Goodness-of-fit RMSE, NSE, and PBIAS performance statistics of the observed—i.e., IDW-interpolated rainfall raster map—and simulated—i.e., radar legend—thresholds taking into account the different colour patterns of the meteorological radar legend.
Table A2. Goodness-of-fit RMSE, NSE, and PBIAS performance statistics of the observed—i.e., IDW-interpolated rainfall raster map—and simulated—i.e., radar legend—thresholds taking into account the different colour patterns of the meteorological radar legend.
Date (dd.mm.yyyy)TimeMax. 6 h Rainfall (mm)RMSENSEPBIAS
15.08.201506:0016.09.23−0.59135.34
12:0027.39.450.772231.60
17.08.201512:0010.09.99−0.7789−81.19
18:006.42.850.662129.64
03.09.201506:0016.89.58−0.069536.11
04.09.201512:0059.441.40−2.2250−22.79
30.09.201506:0020.89.240.871523.05
12:0016.09.270.608814.63
18:0011.04.740.579737.66
14.10.201506:0018.95.150.846017.87
12:0015.05.860.434439.38
10.03.201600:001.40.270.987510.74
01.04.201600:0024.28.960.643653.91
08.04.201612:008.45.470.1734−86.58
18:0015.66.930.454052.20
20.09.201618:0010.44.340.305519.50
23.09.201618:0027.616.98−0.9877−5.11
24.09.201606:000.80.750.4110−123.65
18:0036.819.09−6.0217−15.04
20.10.201600:0093.442.350.359786.39
12:005.62.500.645615.59
22.10.201606:0013.25.000.722729.52
12:004.21.380.884936.53
14.11.201606:0017.410.54−0.097142.08
12:0012.28.640.2936−24.03
18:006.84.16−0.0462−10.10
05.12.201606:0023.68.990.736222.95
12:001.81.44−0.1816−93.84
20.01.201712:0028.417.520.57185.63
21.01.201700:007.83.790.0680−38.67
06:0013.22.780.94109.52
18:0018.07.560.70165.54
27.01.201718:0024.514.35−1.5837−15.46
24.03.201706:003.80.820.8335−11.65
12:0012.03.770.86520.76
18:0021.411.170.799760.22
05.06.201706:0038.415.990.730964.99
12:0051.019.950.792652.22
18:006.41.690.939069.07
15.09.201712:007.44.65−1.6078−80.89
18:0053.426.200.309771.38
19.10.201706:0017.49.330.400263.78
12:0020.16.630.829134.45
18:009.22.320.927975.49
06.02.201806:0010.23.720.87091.13
18:0010.99.97−1.6244−46.31
24.03.201806:008.62.210.768920.83
12:0017.88.860.636911.37
18:0017.46.790.1904−7.45
14.04.201806:001.61.54−0.6957−120.78
12:0011.05.740.787859.80
18:001.20.710.6839−73.59
01.05.201806:0013.67.94−0.5120−4.61
12:0011.310.04−5.6762−145.24
18:0020.07.220.464260.94
13.05.201806:0019.27.640.677248.25
12:003.80.580.0968−244.89
10.09.201806:0016.48.14−0.1452−6.57
12:0033.018.900.692393.54
18.10.201806:006.02.730.4035−6.27
12:006.83.010.726751.87
18:0014.75.600.368623.13
27.10.201812:0013.210.21−0.6154−107.91
18:0042.016.670.174970.63
28.10.201800:0019.07.850.781636.96
12:0024.87.290.587130.00
20.11.201812:0012.03.170.6958−2.35
27.08.201912:0014.54.710.5455−10.51
18:0026.011.160.479480.78
21.10.201900:0015.812.23−0.8259−85.84
06:0014.46.620.561566.65
12:005.12.010.767920.71
22.10.201900:0034.819.790.46402.46
06:0034.216.790.867762.81
12:0011.21.23−0.3087−295.49
08.11.201900:0012.45.570.451429.06
06:004.62.050.3953−48.81
18:002.00.68−3.5196−569.22
22.11.201912:0016.29.18−0.261817.01
18:003.41.220.833121.61
02.12.201906:0016.85.730.522529.25
04.12.201912:0040.818.910.25742.67
18:005.22.79−0.0720−41.81
19.01.202006:001.61.110.5882−49.47
18:0021.05.900.949233.90
22.03.202012:0021.710.240.211932.67
13.04.202006:0011.12.600.893035.14
12:0016.09.360.801579.60
20.04.202000:004.31.760.748217.41
06:0017.29.040.423760.34
12:0013.73.780.415615.91
18:0031.414.650.566267.41
29.08.202012:0022.27.960.283832.47
18:0018.06.880.817878.23
30.08.202112:0022.27.620.722921.23
01.10.202112:0014.64.430.590725.26
18:0010.05.100.776679.70
30.10.202106:007.32.190.5441−23.15
12:0021.48.080.534161.97
18:001.00.440.7635−143.45
10.11.202100:0017.411.390.612186.23
06:0017.67.030.730051.50
12:0017.26.150.543347.45
18:007.44.090.29677.21
11.11.202106:0017.06.050.553344.87
12:0016.34.920.655642.23
18:0010.44.250.548926.02
23.11.202106:0036.512.250.397454.68
12:0016.48.720.4762−10.23
18:0024.010.900.797457.30
25.11.202112:0011.53.120.725613.63
28.11.202100:002.61.610.1114−72.87
06:003.00.870.89471.95
12:005.21.890.6879−10.62

References

  1. IPCC. Managing the Risks of Extreme Events and Disasters to Advance Climate Change Adaptation. A Special Report of Working Groups I and II of the Intergovernmental Panel on Climate Change; Field, C.B., Barros, V., Stocker, T.F., Qin, D., Dokken, D.J., Ebi, K.L., Mastrandrea, M.D., Mach, K.J., Plattner, G.-K., Allen, S.K., et al., Eds.; University Press: Cambridge, UK; New York, NY, USA, 2012; p. 582. [Google Scholar]
  2. Verma, S.; Bhattarai, R.; Bosch, N.S.; Cooke, R.C.; Kalita, P.K.; Markus, M. Climate Change Impacts on Flow, Sediment and Nutrient Export in a Great Lakes Watershed Using SWAT. CLEAN-Soil Air Water 2015, 43, 1464–1474. [Google Scholar] [CrossRef]
  3. IPCC. Climate Change 2023: Synthesis Report (Full Volume) Contribution of Working Groups I, II and III to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change; Core Writing Team, Lee, H., Romero, J., Eds.; IPCC: Geneva, Switzerland, 2023; p. 184. [Google Scholar] [CrossRef]
  4. Ortiz, E.; Todini, E. Combinación Bayesiana de datos pluviométricos e imágenes del Radar Meteorológico mediante Block Kriging y Filtro de Kalman. Aplicación en España. In Proceedings of the XXXI Jornadas Científicas de la AME y el 11º Encuentro Hispano-Luso de Meteorología, Sevilla, Spain, 1–3 March 2010; p. 6. [Google Scholar]
  5. Erena, M.; López, J.; García, P.; Caro, M.; Belda, F.; Palenzuela, J.; Toledano, F.; Torralba, P.; González-Barbera, G.; García-Pintado, J. Estimación de precipitación combinada radar-pluviómetros y publicación mediante servicios OGC. In Proceedings of the XV Congreso Nacional de Tecnologías de la Información Geográfica, AGE-CSIC, Madrid, Spain, 19–21 September 2012; pp. 451–458. [Google Scholar]
  6. Pulkkinen, S.; Chandrasekar, V.; Harri, A.M. Nowcasting of Precipitation in the High-Resolution Dallas–Fort Worth (DFW) Urban Radar Remote Sensing Network. IEEE J. Sel. Top. Appl. Earth Obs. Remote. Sens. 2018, 11, 2773–2787. [Google Scholar] [CrossRef]
  7. Einfalt, T.; Arnbjerg-Nielsen, K.; Golz, C.; Jensen, N.E.; Quirmbach, M.; Vaes, G.; Vieux, B. Towards a roadmap for use of radar rainfall data in urban drainage. J. Hydrol. 2004, 299, 186–202. [Google Scholar] [CrossRef]
  8. Thorndahl, S.; Einfalt, T.; Willems, P.; Nielsen, J.E.; ten Veldhuis, M.C.; Arnbjerg-Nielsen, K.; Rasmussen, M.R.; Molnar, P. Weather radar rainfall data in urban hydrology. Hydrol. Earth Syst. Sci. 2017, 21, 1359–1380. [Google Scholar] [CrossRef]
  9. Liguori, S.; Rico-Ramirez, M.; Schellart, A.; Saul, A. Using probabilistic radar rainfall nowcasts and NWP forecasts for flow prediction in urban catchments. Atmos. Res. 2012, 103, 80–95. [Google Scholar] [CrossRef]
  10. Heuvelink, D.; Berenguer, M.; Brauer, C.C.; Uijlenhoet, R. Hydrological application of radar rainfall nowcasting in the Netherlands. Environ. Int. 2020, 136, 105431. [Google Scholar] [CrossRef] [PubMed]
  11. Blanco-Gómez, P.; Amurrio-Garcia, C.; Jiménez-García, J.L.; Cecilia, J.M. CPR Algorithm—A new interpolation methodology and QGIS plugin for Colour Pattern Regression between aerial images and raster maps. SoftwareX 2023, 22, 101356. [Google Scholar] [CrossRef]
  12. Agencia Estatal de Meteorología. Interpretación: Radar. 2024. Available online: https://www.aemet.es/es/eltiempo/observacion/radar/ayuda (accessed on 6 September 2024).
  13. Estrany-Planas, P.; Blanco-Gómez, P.; Vilarrasa, V. Aplicación del análisis regional de frecuencias a las series diarias de precipitación usando la metodología de los L-moments en la ciudad de Palma para la caracterización de la lluvia de diseño. In Proceedings of the VII Jornadas de Ingeniería del Agua—La resiliencia de las Infraestructuras Hidráulicas ante el Cambio Climático, Cartagena, Spain, 18–19 October 2023; pp. 359–372. [Google Scholar] [CrossRef]
  14. Empresa Municipal de Agua y Alcantarillado, S.A. (EMAYA). La Història i les Històries d’EMAYA; Secció de Comunicació EMAYA: Palma de Mallorca, Spain, 2008; p. 51. Available online: https://www.emaya.es/media/2924/la-histo-ria-i-les-histo-ries-demaya-web.pdf (accessed on 6 September 2024).
  15. Grimalt, M.; Rosello, J. Inundaciones en la Ciudad de Palma de Mallorca: Distribución de los Episodios y Tipos de Tiempo Asociados. In El Clima: Aire, Agua, Tierra y Fuego; Montavez, J.P., Ed.; Asociación Española de Climatología: Cartagena, Spain, 2018; pp. 245–258. [Google Scholar]
  16. Ministerio de la Presidencia, Relaciones con las Cortes y Memoria Democrática. Real Decreto 665/2023, de 18 de Julio, por el que se Modifica el Reglamento del Dominio Público Hidráulico, Aprobado por Real Decreto 849/1986, de 11 de Abril; el Reglamento de la Administración Pública del Agua, Aprobado por Real Decreto 927/1988, de 29 de Julio; y el Real Decreto 9/2005, de 14 de Enero, por el que se Establece la Relación de Actividades Potencialmente Contaminantes del Suelo y los Criterios y estáNdares para la Declaración de Suelos Contaminados. 2023. Available online: https://www.boe.es/diario_boe/txt.php?id=BOE-A-2023-18806 (accessed on 6 September 2024).
  17. Council of European Communities. Council Directive 91/271/EEC of 21 May 1991 Concerning Urban Waste-Water Treatment. 1991. Available online: http://data.europa.eu/eli/dir/1991/271/oj (accessed on 6 September 2024).
  18. Torrens Calleja, J.M.; Rosselló Geli, J.; Grimalt Gelabert, M. Recopilación de Información Vinculada a Temporales de Viento, Precipitaciones Torrenciales e Inundaciones en la Ciudad de Palma de Mallorca Entre los Años 2000 y 2015; RUA: Alicante, Spain, 2016. [Google Scholar]
  19. YACU. Estudio de Caracterización del Régimen Extremo de Precipitaciones en la Isla de Mallorca. Informe Final; Direcció General de Recursos Hídrics: Conselleria de Medi Ambient, Govern Balear: Illas, Spain, 2002; p. 621. [Google Scholar]
  20. BALEARSMETEO. Balears Meteo, Network of Amateur Weather Stations in the Balearic Islands. 2024. Available online: http://balearsmeteo.com/ (accessed on 6 September 2024).
  21. Blanco-Gómez, P.; Ortiz-Vallespi, J.; Estray-Plana, P.; Orihuela-Martínez, J. Rainfall Frequency Analysis (RainFA) Package. 2023. Available online: https://github.com/vielca/RainFA (accessed on 6 September 2024).
  22. Cánovas-García, F.; García-Galiano, S.; Alonso-Sarría, F. Assessment of Satellite and Radar Quantitative Precipitation Estimates for Real Time Monitoring of Meteorological Extremes Over the Southeast of the Iberian Peninsula. Remote Sens. 2018, 10, 1023. [Google Scholar] [CrossRef]
  23. Marshall, J.S.; Palmer, W.M.K. The distribution of raindrops with size. J. Atmos. Sci. 1948, 5, 165–166. [Google Scholar] [CrossRef]
  24. Uijlenhoet, R. Raindrop size distributions and radar reflectivity-rain rate relationships for radar hydrology. Hydrol. Earth Syst. Sci. 2001, 5, 615–627. [Google Scholar] [CrossRef]
  25. Battan, L. Radar Observation of the Atmosphere; The University of Chicago Press: Chicago, IL, USA, 1973; p. 324. [Google Scholar]
  26. List, R. A Linear Radar Reflectivity–Rainrate Relationship for Steady Tropical Rain. J. Atmos. Sci. 1988, 45, 3564–3572. [Google Scholar] [CrossRef]
  27. SNAP. ESA Sentinel Application Platform v9.0.0. 2022. Available online: http://step.esa.int (accessed on 6 September 2024).
  28. Gillies, S. Rasterio: Geospatial Raster I/O for Python Programmers. 2013. Available online: https://github.com/rasterio/rasterio (accessed on 6 September 2024).
  29. Moriasi, N.D.; Arnold, J.G.; Van Liew, W.M.; Bingner, R.L.; Harmel, R.D.; Veith, T.L. Model Evaluation Guidelines for Systematic Quantification of Accuracy in Watershed Simulations. Trans. ASABE 2007, 50, 885–900. [Google Scholar] [CrossRef]
  30. Brighenti, T.M.; Bonumá, N.B.; Grison, F.; de Almeida Mota, A.; Kobiyama, M.; Chaffe, P.L.B. Two calibration methods for modeling streamflow and suspended sediment with the swat model. Ecol. Eng. 2019, 127, 103–113. [Google Scholar] [CrossRef]
  31. Zambrano-Bigiarini, M.; Nauditt, A.; Birkel, C.; Verbist, K.; Ribbe, L. Temporal and spatial evaluation of satellite-based rainfall estimates across the complex topographical and climatic gradients of Chile. Hydrol. Earth Syst. Sci. 2017, 21, 1295–1320. [Google Scholar] [CrossRef]
  32. Vieux, B.E.; Bedient, P.B. Estimation of Rainfall for Flood Prediction from WSR-88D Reflectivity: A Case Study, 17–18 October 1994. Weather Forecast. 1998, 13, 407–415. [Google Scholar] [CrossRef]
  33. Aziding, K.; Badron, K.; Basri, A.B.; Ismail, A.F.; Ahmad, Y.A. Enhanced Technique for Prediction of Z-R Relationship in Tropical Region. J. Physics Conf. Ser. 2023, 2559, 012009. [Google Scholar] [CrossRef]
  34. Orellana-Alvear, J.; Célleri, R.; Rollenbeck, R.; Bendix, J. Analysis of Rain Types and Their Z–R Relationships at Different Locations in the High Andes of Southern Ecuador. J. Appl. Meteorol. Climatol. 2017, 56, 3065–3080. [Google Scholar] [CrossRef]
  35. Yeo, J.; Lee, Y.; Ong, J. Radar measured rain attenuation with proposed Z–R relationship at a tropical location. AEU-Int. J. Electron. Commun. 2015, 69, 458–461. [Google Scholar] [CrossRef]
Figure 1. Daily discharge calendar for the period 2015–2021. Green colour intensity indicates the number of discharges affected. Source: EMAYA.
Figure 1. Daily discharge calendar for the period 2015–2021. Green colour intensity indicates the number of discharges affected. Source: EMAYA.
Remotesensing 16 03496 g001
Figure 2. Location map: location of the city of Palma in the island of Mallorca (Spain); hydro-meteorological and radar stations of the present study and sewage overflow discharges (in orange) into the Mediterranean Sea.
Figure 2. Location map: location of the city of Palma in the island of Mallorca (Spain); hydro-meteorological and radar stations of the present study and sewage overflow discharges (in orange) into the Mediterranean Sea.
Remotesensing 16 03496 g002
Figure 3. Flowchart for the methodology adopted in the present study.
Figure 3. Flowchart for the methodology adopted in the present study.
Remotesensing 16 03496 g003
Figure 4. Monthly distribution of moderate to extreme rainfall storm events in Palma de Mallorca (Spain) for the period 2015–2021.
Figure 4. Monthly distribution of moderate to extreme rainfall storm events in Palma de Mallorca (Spain) for the period 2015–2021.
Remotesensing 16 03496 g004
Figure 5. CPR algorithm results including, for each date and time evaluated, the calculated R, G, and B correlation coefficients and the maximum 6 h cumulative rainfall value observed at any of the existing gauging stations.
Figure 5. CPR algorithm results including, for each date and time evaluated, the calculated R, G, and B correlation coefficients and the maximum 6 h cumulative rainfall value observed at any of the existing gauging stations.
Remotesensing 16 03496 g005
Figure 6. Maximum 6 h observed rainfall for meteorological radar images that include only Blue and also Green values. (a) Shows the dispersion of both CPR algorithm parameters, and (b) presents the box and whisker results.
Figure 6. Maximum 6 h observed rainfall for meteorological radar images that include only Blue and also Green values. (a) Shows the dispersion of both CPR algorithm parameters, and (b) presents the box and whisker results.
Remotesensing 16 03496 g006
Figure 7. Performance rating result analysis for (a) NNSE and (b) PBIAS statistics.
Figure 7. Performance rating result analysis for (a) NNSE and (b) PBIAS statistics.
Remotesensing 16 03496 g007
Figure 8. Minimum and maximum values of the 6 h observed rainfall for meteorological radar image colour ranges. Corresponding to (a) Satisfactory, (b) Good, and (c) Very good performance models.
Figure 8. Minimum and maximum values of the 6 h observed rainfall for meteorological radar image colour ranges. Corresponding to (a) Satisfactory, (b) Good, and (c) Very good performance models.
Remotesensing 16 03496 g008
Figure 9. Comparison of the resulting curve with several Z-R relationships adopted from [7,32].
Figure 9. Comparison of the resulting curve with several Z-R relationships adopted from [7,32].
Remotesensing 16 03496 g009
Table 1. Evaluation performance metrics.
Table 1. Evaluation performance metrics.
StatisticEquationValue Range
RMSE i = 1 n ( O i S i ) 2 n [0, + )
NSE 1 i = 1 n ( O i S i ) 2 i = 1 n ( O i O ¯ ) 2 ( ,1]
PBIAS i = 1 n ( O i S i ) · 100 i = 1 n O i ( , + )
Table 2. Colour coding of cumulative rainfall values from meteorological radar for storm events before and after 1 January 2016.
Table 2. Colour coding of cumulative rainfall values from meteorological radar for storm events before and after 1 January 2016.
6 h Rainfall (mm)Before 20162016–Onwards *
0.5–1Various[0,0,250]
1–2[254,254,254][0,0,205]
2–4[254,254,128][0,0,155]
4–8[254,128,128][0,0,105]
8–16[0,0,254][0,100,0]
16–32[0,128,254][0,250,0]
32–64[128,254,254][250,250,0]
64–128[192,0,0][250,165,0]
128–256[155,187,89][250,0,0]
* Note that blueish colours explain cumulative values less than 8 mm, greenish colours cover the range of values between 8 and 32 mm, and a combination of red and green colours is used for values between 32 and 128 mm. The resulting colours are graphically represented as follows: Remotesensing 16 03496 i001
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Blanco-Gómez, P.; Estrany-Planas, P.; Jiménez-García, J.L. Cumulative Rainfall Radar Recalibration with Rain Gauge Data Using the Colour Pattern Regression Algorithm QGIS Plugin. Remote Sens. 2024, 16, 3496. https://doi.org/10.3390/rs16183496

AMA Style

Blanco-Gómez P, Estrany-Planas P, Jiménez-García JL. Cumulative Rainfall Radar Recalibration with Rain Gauge Data Using the Colour Pattern Regression Algorithm QGIS Plugin. Remote Sensing. 2024; 16(18):3496. https://doi.org/10.3390/rs16183496

Chicago/Turabian Style

Blanco-Gómez, Pablo, Pau Estrany-Planas, and José Luis Jiménez-García. 2024. "Cumulative Rainfall Radar Recalibration with Rain Gauge Data Using the Colour Pattern Regression Algorithm QGIS Plugin" Remote Sensing 16, no. 18: 3496. https://doi.org/10.3390/rs16183496

APA Style

Blanco-Gómez, P., Estrany-Planas, P., & Jiménez-García, J. L. (2024). Cumulative Rainfall Radar Recalibration with Rain Gauge Data Using the Colour Pattern Regression Algorithm QGIS Plugin. Remote Sensing, 16(18), 3496. https://doi.org/10.3390/rs16183496

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop