Next Article in Journal
A Study on Climate-Driven Flash Flood Risks in the Boise River Watershed, Idaho
Previous Article in Journal
Analysis of Water Losses and Assessment of Initiatives Aimed at Their Reduction in Selected Water Supply Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Tropical Andes Radar Precipitation Estimates Need High Temporal and Moderate Spatial Resolution

by
Mario Guallpa
1,2,3,*,
Johanna Orellana-Alvear
2,4 and
Jörg Bendix
4
1
Subgerencia de Gestión Ambiental, Empresa Pública Municipal de Telecomunicaciones, Agua Potable, Alcantarillado y Saneamiento de Cuenca, Benigno Malo 7-78 y Mariscal Sucre, Cuenca 010104, Ecuador
2
Departamento de Recursos Hídricos y Ciencias Ambientales, Universidad de Cuenca, Av. 12 de Abril, Cuenca 010207, Ecuador
3
Facultad de Ingeniería, Universidad de Cuenca, Av. 12 de Abril, Cuenca 010203, Ecuador
4
Laboratory for Climatology and Remote Sensing (LCRS), Faculty of Geography, University of Marburg, Marburg 35032,Germany
*
Author to whom correspondence should be addressed.
Water 2019, 11(5), 1038; https://doi.org/10.3390/w11051038
Submission received: 21 March 2019 / Revised: 6 May 2019 / Accepted: 15 May 2019 / Published: 18 May 2019
(This article belongs to the Section Hydrology)

Abstract

:
Weather radar networks are an excellent tool for quantitative precipitation estimation (QPE), due to their high resolution in space and time, particularly in remote mountain areas such as the Tropical Andes. Nevertheless, reduction of the temporal and spatial resolution might severely reduce the quality of QPE. Thus, the main objective of this study was to analyze the impact of spatial and temporal resolutions of radar data on the cumulative QPE. For this, data from the world’s highest X-band weather radar (4450 m a.s.l.), located in the Andes of Ecuador (Paute River basin), and from a rain gauge network were used. Different time resolutions (1, 5, 10, 15, 20, 30, and 60 min) and spatial resolutions (0.5, 0.25, and 0.1 km) were evaluated. An optical flow method was validated for 11 rainfall events (with different features) and applied to enhance the temporal resolution of radar data to 1-min intervals. The results show that 1-min temporal resolution images are able to capture rain event features in detail. The radar–rain gauge correlation decreases considerably when the time resolution increases (r from 0.69 to 0.31, time resolution from 1 to 60 min). No significant difference was found in the rain total volume (3%) calculated with the three spatial resolution data. A spatial resolution of 0.5 km on radar imagery is suitable to quantify rainfall in the Andes Mountains. This study improves knowledge on rainfall spatial distribution in the Ecuadorian Andes, and it will be the basis for future hydrometeorological studies.

1. Introduction

An adequate representation of spatio-temporal rainfall variability is of utmost importance for many meteorological, hydrological, and ecological studies [1,2]. This particularly holds true for the tropical Andes Mountains, which are affected by a high spatio-temporal variability of precipitation but are at the same time characterized by a lack of hydrometeorological monitoring [3,4]. In this regard, rainfall data acquired from different remote sensors, such as weather radar systems, have become more accessible over the last few years [5]. This technology has allowed us to study rainfall with high spatial and temporal resolutions [6], and it can be used in several applications, such as input data in rainfall forecasting studies [7,8,9,10,11], different hydrological models [12,13], early warning systems for floods [14], erosion studies [1,15], research related to atmospheric chemistry [16], and many others. Meanwhile, cost-effective systems derived from ship radar technology are available to establish radar networks with limited funds at hand [17].
The estimation of the rainfall from radar measurements, particularly in the case of non-Doppler or non-polarimetric systems, is affected by different sources of error as discussed in several studies [10,17,18,19]. This particularly holds true for mountain areas [20,21,22,23]. These errors can be grouped into three classes [24,25]. The first class includes ground clutter and beam blocking, which can be removed by adjusting the radar reflectivity maps. The second class is related to the transformation of radar reflectivity to rain rate. The differences between radar and ground rain rate (bias) can be improved using rain gauge data collected over long periods of time. To conduct bias correction of radar data using rain gauge observations, different methodologies have been used successfully in the past, such as radar–rain gauge ratios, as well as Bayesian and geostatistical techniques [23,26,27]. The final class comprises the correction of sampling errors associated with the temporal spacing of radar scans.
Rainfall accumulation maps derived from quantitative precipitation estimation (QPE) with radar data often represent the actual spatial distribution of rainfall incorrectly due to the temporal sampling frequency [28,29]. According to [30,31], the effect of this temporal sampling could be more significant than other error sources combined, particularly when using temporally low-resolution data. Generally, the change rate of rainfall characteristics (e.g., intensity, spatial coverage, rain cell structure) between images increases as the temporal sampling interval increases [32,33]. This is especially critical for highly dynamic convective cells of limited spatial extension, which is the main source of rainfall in the tropics.
In more detail, the influence of specific factors on the assessment of radar rainfall accumulation in comparison to ground data have been found: (i) The sampling frequency has a negative impact depending on the radar spatial resolution. The authors of [28] found a spatial discontinuity in rainfall patterns for hourly accumulation maps by using a sampling frequency of 5 min and a grid resolution of 0.5 km. Similar results were found by [25,31], who also observed an error reduction when using a coarser spatial resolution. Nevertheless, [30] found that for 1 km and 5-min sampling, the error in 10-min rainfall accumulation compared with the high-resolution data (0.2 km and 50 s) is of a magnitude of 20% of the mean rainfall accumulation, increasing to 35% for 2 km and 10-min sampling. Additionally, the temporal sampling interval has greater influence on error than the spatial resolution [25]. For instance, [30] found that for 2 km and 5-min sampling, the error is 50% of the 10-min rainfall accumulation, which is twice that of a 1-min advection interpolation. A solution to reduce the error is to increase the temporal resolution. The authors of [24] increased the sampling frequency to 2 min with 1-km spatial resolution, whereby the mean error was 14% and 8% for scattered and widespread rainfall, respectively. (ii) The storm speed highly influences the cumulative QPE. The authors of [33] determined that the mean error (ME) of the radar–rain gauge correlation decreases with the increase of the storm speed. However, the ME varies throughout the radar range, because the intensity and spatial extent of the precipitation highly influence the error magnitude. Nonetheless, other error metrics increase considerably with storm speed. For instance, the mean absolute error (MAE) between radar and rain gauge ranges from 0.5 to 1 mm h−1 with respect to the events’ average rain intensity (5 mm h−1) [31]. To improve precipitation estimation using radar images, researchers [25,28,34,35] have recommended the use of temporal interpolation methods to add artificial images between the scan times through precipitation advection methods. These methods are proven to reduce the mean absolute error by more than 60%.
Despite the progress in improving cumulative rainfall estimation using radar data, the abovementioned factors, such as radar sample frequency, storm speed, and radar spatial resolution, still challenge the proper evaluation and correction of QPE [28]. Even though the scanning frequency is not a main factor correctly capturing the rain rate at slow storm movements, in particular, the current needed to obtain products with a high spatial resolution (e.g., 0.5 km) increases the difficulty of determining a proper rainfall estimate for fast-tracking storms.
Fast-tracking storms have a strong influence during times of rainfall accumulation of hydrological interest (e.g., 1–3 h) that need more exploration [32]. Furthermore, many hydrological-related studies, such as [5,14,36] have been performed in relatively low altitudes (0–1000 m a.s.l.). As a consequence, the sampling frequency effect of radar data on the cumulative QPE in higher altitude mountainous areas, such as watersheds located in the Tropical Andes, is still unknown. In this context, the aim of the current study is to analyze the impact of spatial and temporal resolutions of radar data on QPE accumulation for different rainfall features (e.g., rain field speed, rainfall intensity, spatial rainfall coverage, and rain total volume precipitated in the sub-basin). The research was conducted using the highest global X-band radar (CAXX), which is part of the first weather radar network installed in the Tropical Andes of Ecuador [17]. In more detail, in this study, a precipitation advection model based on the optical flow method was validated regarding its applicability in the complex terrain of the Andes. The validated scheme was then used to enhance the original temporal resolution of radar data (5 min) to 1-min temporal resolution. Finally, the effect of using different data recording frequencies (1–60 min) and spatial resolutions (0.5, 0.25, and 0.1 km) to derive rainfall accumulation was assessed by means of several error metrics.

2. Materials

2.1. Study Area

The study area comprises the middle and upper part of the Paute River basin (Figure 1 and Table 1), which is known as the Cuenca sub-basin. The Paute basin has an area of 6148 km², which is located in the depression between the east and west of the Andean escarpment of southern Ecuador [37]. It has a complex topography with elevations ranging from 1840 to 4680 m a.s.l. The basin possesses a high ecological and economic value and provides several ecosystem services, such as water supply and regulation, biodiversity conservation, and storage of carbon [38]. About 40% of the basin is covered by Páramo vegetation (3300–4500 m a.s.l.), which plays a key role in water regulation in the region and thus makes this fragile ecosystem highly important for the local population [4,39]. The middle and upper basins (Cuenca sub-basin, 1605.4 km²) are sources of potable water to more than 580,000 inhabitants and indirectly to more than 2000 users through irrigation, agriculture, cattle raising, and other activities [40].
The climate in the mountains around the Paute basin is influenced by the Pacific Ocean regime and by the Amazon basin air masses originating from the east [41]. The basin is characterized by a high spatio-temporal precipitation variability [4]. Different rain regimes have been identified based on bimodal inter-Andean seasonality in the basin’s central part, with higher rainfall in March–April and October–November and a dry season in July–August [37,42]. The basin’s upper part (2600–4680 m a.s.l.) shows bimodal seasonality but with a less pronounced dry season. The average annual rainfall is between 1100 and 1300 mm [4,43,44]. The temperature in the upper part of the basin has a daily average of 8.7 °C (Tmax = 17.8 °C and Tmin = 1.8 °C) [45]. The average wind speed at 3900 m a.s.l. is 4 m s−1, with maximum peaks of 20 m s−1 in the June–September period.

2.2. Weather Radar Data

The highest weather radar in the world, named CAXX, is part of the new radar network installed in southern Ecuador, RadarNet-Sur [17]. The network consists of three X-band, single-polarized radars. The sensor is a Rainscanner-RS120 Gematronik, located on Paragüillas peak at 4450 m a.s.l. on the northern border of the Cajas National Park near the city of Cuenca (Figure 2).
The CAXX radar covers the Paute River sub-basins and has been operating since April 2015. It has a transmission frequency of 9410 ± 30 MHz, a bandwidth of 2.5 MHz, a pulse length of 1.2–0.5 μs, and a wavelength of 3 cm. The parabolic antenna of the system has a vertical and horizontal beam width of 2°. The radar has a maximum range of 100 km and provides polar images of radar reflectivity every 5 min. The built-in software of the equipment allows for cartesian transformation and produces plan position indicator images (.ppi) at the same temporal frequency with three different spatial resolutions (0.1, 0.25, and 0.5 km).
The raw reflectivity data needed some corrections to obtain calibrated images, as well as later adjustment using rain gauge data. (i) The built-in software allows for the clutter correction of the raw reflectivities. The removed clutter is the reflection signals produced by solid objects, such as buildings and mountains lying in the radar beam. The resulting blank bins were interpolated with cubic splines (a method incorporated in the RainView Analyzer software) to fill them with values estimated from uncontaminated surrounding reflectivities. According to [1], this method minimizes false estimates within the clutter fields and thus performs better than the bilinear interpolation used by [20]. Then, the corrected reflectivity images were projected to cartesian coordinates with three different spatial resolutions. (ii) To convert reflectivity (Z) to rain rate (R), the site-specific Z–R relationship found by [46] Z = 204R1.57 was applied to the radar images. (iii) Afterwards, a bias correction between data of the rain gauge network (G) and the radar (R) CAXX images (RG–RR) was implemented in order to obtain adjusted rainfall products. The difference in the total rainfall between the rain gauge and its corresponding radar pixel was interpolated using the inverse distance weight method. These differences were added only on radar pixels that detected rainfall.
We assumed that bias correction is able to compensate for the radar uncertainty caused by radar calibration and signal attenuation. The radar processing library ωradlib [47,48] was used to obtain the adjusted radar images.

2.3. Rain Gauge Data and Rain Events

The rain gauge network is composed of 38 automatic tipping-bucket gauges with 0.254- and 0.1-mm resolution and 5-min time step records. All stations were regularly maintained and operated by ETAPA EP Staff and researchers at the Universidad de Cuenca.
A total of 11 rainfall events were examined from March 2015 to December 2016 (Table 2 and Figure 3). For the rain event selection, three criteria were considered: (i) events with very high rain intensities, (ii) events taking place in the wet or dry season of the study area, with the objective of analyzing the entire rainfall seasonality, and (iii) events with a large spatial extent, covering at least 20 rain gauges within the study area. For each event, the maximum rain intensity was calculated considering a temporal aggregation of 5 min and 1 h. Details for all events are shown in Table 2.

2.4. Re-Analysis Data

Re-analysis data (FNL) of wind speed and direction were used for the comparison with the calculated streamlines of the rain cells. These data were obtained from the National Centers for Environmental Prediction (NCEP) of the National Oceanic and Atmospheric Administration (NOAA). The data have spatial and temporal resolutions of 1° and 6 h, respectively [49], and were interpolated for the location of the radar at 4450 m a.s.l. (CAXX Radar location height). To validate the re-analysis data, wind speed and direction records of spatially distributed climate stations over the study area were used (Figure 1). According to the re-analysis data, winds move from east to west in 91% of all cases (from the Amazon toward the Andes), while in 9% of all cases, the winds originate from the Pacific Ocean (Figure 4).

3. Methodology

The workflow of the study is depicted in Figure 5. The first step was to enhance the temporal resolution of the radar data. The radar captures rain images every 5 min. However, for this study, 1-min-interval rain images were required. Therefore, rain motion vectors between two consecutive image slots were calculated through a rain advection model and validated for the local conditions. Then, 1-min-interval images were calculated for each rain event by using the motion vector fields. The second step was the estimation of hourly rain maps in the different spatial and temporal resolutions (0.50, 0.25, and 0.10 km; 1, 5, 10, 15, 20, 30, and 60 min), including the artificial subslot 1-min-interval images. Then, a statistical evaluation of hourly radar rainfall maps was conducted with hourly rain gauge totals to show the general accuracy of QPE. Finally, the impact of different temporal and spatial resolutions of radar data on QPE was analyzed for different rainfall characteristics (horizontal rain field speed, rain total volume, spatial coverage percentage, and maximum rain intensity).

3.1. Rain Field Advection

The advection field describes the horizontal motion vectors of rain cells over a certain time period. This movement reflects the rain cells’ position change between two consecutive image slots at time t and time t + Δt.
Methods such as the tracking radar echoes by correlation (TREC), which uses the block matching technique [50] and the optical flow using intensity gradient only [51] have already been used to calculate advection. However, these methods showed some inconsistencies in the magnitude of the movement vector [52]. In contrast, the optical flow method using intensity gradient information on several image scales (PyrLK) is a more robust method for an advection model using consecutive images [8,52,53]. The algorithm for this method was introduced by [54] and is based on the Lucas Kanade algorithm [55]. The author of [54] uses an information gradient on different image scales and looks for rain cell speed vectors (magnitude and direction) of two consecutive images. The algorithm has three assumptions: (i) constant radar reflectivities between the pixels of an identified object (rain cells) of the two consecutive images, (ii) small movements between the image reflectivity cells, and (iii) spatial coherence.
Given two images where F(x) and G(x) are respective reflectivity values at a location x (where x is a two-dimensional (2D) vector), the algorithm looks for a displacement h where G (x) = F (x + h) [54,55]. In the case that h is small (small rain cell motion between two consecutive images), it can be calculated iteratively using Equation (1).
h 0 = 0
h k + 1 = h k + x w ( x ) F ( x + h k ) [ G ( x ) F ( x + h k ) x w ( x ) F ( x + h k ) 2
where k is the interaction number until the best h is found by using the Newton–Raphson method. F’(x) and G’(x) are the F(x) and G(x) derivatives calculated by finite differences.
The weight w(x) is defined as
w ( x ) = 1 | G ( x ) F ( x ) |
The method does not calculate the advection vector for each image pixel directly. It only takes the pixels for which the algorithm is most likely to converge. To obtain the entire field information, the calculated motion vector of the selected pixels must be interpolated using the inverse distance weighting method [52]. For the implementation of the optical flow method, the Open Computer Vision (OpenCV) library [56] was used.
To validate PyrLK, two criteria were used: (i) the rain motion vectors were compared with the wind vectors (u, v) determined by re-analysis data using the Pearson correlation coefficient (r), and (ii) by using the rain displacement vectors of each pixel, the radar image was projected for the next 5 min. Then, this predicted radar image was compared with the one observed by the radar at the same time interval. The probability of detection (POD), false alarm (FAR), and frequency of bias (FBI) indices were calculated to quantify the coincidences between the projected and recorded images every 5 min. [8,18,57,58].
Equations (2)–(4) define the formula of these indicators, where tp is the number of true positives (i.e., the pixel number for which the projected and observed images detect precipitation), fp is the number of false positives (i.e., the pixel number for which the projected image detects precipitation while the radar does not), and fn is the number of false negatives (i.e., the pixel number for which the projected image does not detect precipitation while the radar does).
POD = tp tp + fn
FAR = fp tp + fp
FBI = tp + fp tp + fn
Thus, the POD is the precipitation fraction that was correctly detected (POD = 1 is a perfect detection). The FAR measures the precipitation detection fraction that was a false alarm, in other words, the precipitation that was mistakenly detected by the projected image (FAR = 0 indicates perfect detection). The FBI is the radius between the precipitation detected by the projected images in regards to the observations detected by the radar (FBI = 1 indicates perfect detection).

3.2. Derivation of 1-Min-Interval Images

For each rainfall event, the derivation of 1-min-interval images between two 5-min-interval consecutive radar images is performed in two steps according to [26]: (i) the first step is to determine the horizontal rain motion vectors in the X and Y direction for each pixel of the rainfall images (e.g., estimated by the optical flow method explained in Section 3.1). (ii) The motion vectors are used to calculate the artificial 1-min sub-slot images through a linear interpolation of the motion vectors. This procedure was performed for each spatial resolution of the radar images (0.50, 0.25, and 0.1 km). We should point out that the use of a linear interpolation scheme is supported, because rain patterns are preserved between two consecutive images even at high wind speeds [32].
The rain field linear projection at time t + Δti, between time t and t + Δt, is given by Equation (5), where Rr(x,y,t) denotes the two-dimensional map at time t, x and y represent the coordinates for each pixel, and a represents the displacement between two 5-min-interval consecutive images (b = before, f = forward) [28].
R r ( x , y , t + Δ t i ) = w 1 R r ( x + a x , f ,   y + a y , f , t ) + w 2 R r ( x + a x , b , y + a y , b , t + Δ t )
The displacement a (Equations (6) to (9)) is defined as the number of pixels based on the product resolution (0.1, 0.25, and 0.5 km), and w is a weight factor (Equations (10) and (11)) that is calculated by the distance proportion equivalent to the time contribution within the period interval (5 min).
a x , f = Δ t i Δ t × Δ x
a y , f = Δ t i Δ t × Δ y
a x , b = ( 1 Δ t i Δ t ) × Δ x
a y , b = ( 1 Δ t i Δ t ) × Δ y
w 1 = 1 Δ t i Δ t
w 2 = Δ t i Δ t

3.3. Impact of Spatial and Temporal Resolutions of the Radar Data on Cumulative QPE

The hourly-accumulated rainfall observations (G) at rain gauge locations are compared with the corresponding radar pixel values (R) from the hourly accumulation maps. For the evaluation, the Pearson correlation coefficient (r), the bias index (BIAS), and the root-mean-square error (RMSE) were determined between G and R [18,58]. In addition, the POD, FAR, and FBI indices were calculated to quantify the coincidences between data detected by radar and by rain gauges.

3.4. Rainfall Features

Four rainfall features were analyzed: rainfall field motion, rain total volume, spatial coverage percentage, and maximum intensity. The rainfall events were split into three ranges according to the average motion of the rainfall events (1–2.5, 2.5–5, and >5 m s−1). Then, we analyzed whether the radar–rain gauge relations show similar behavior under these three speed ranges. The POD and BIAS indexes were used to analyze the R–G relation for the three spatial and seven temporal resolutions.
To the second and third rainfall feature, the rain total volume (m3) and spatial coverage percentage were calculated for each event over the Cuenca sub-basin (1604.5 km2). The two features were calculated using the entire spatio-temporal configuration of the radar data. Next, with the aim to illustrate the results of the radar data, the same steps for total rain and coverage were estimated using the rain gauge network data only. Despite their use in the recent literature, Kriging methods could not be applied in this study for two reasons: (i) the available rain data lack a normal distribution and (ii) the rain data could not be adjusted to any theoretical semivariograms. Therefore, in spite of their limitations, the Thiessen and inverse distance weight (IDW) methods were used to generate the interpolated rain maps.
Finally, the maximum rainfall intensities (mm h−1) of the events were analyzed. The analysis and discussion focused on the influence that maximum rainfall intensities have on the speed of the rain fields, as well as the total volume and the spatial coverage of the rainfall.

4. Results

4.1. Validation of Optical Flow Method

The average speed of most of the events did not exceed 6 m s−1 (Figure 6), which is in agreement with the wind records measured by climate stations and reported by [44,45] for the study area. Most of the rain fields moved from east to west according to the predominant wind direction.
Figure 7 shows that there is a strong correlation between the events’ average motion and the wind speed. Although the correlation coefficient r was high at 0.91, the wind speeds were greater than the advection velocity of the rain cells. For instance, the maximum wind speed for event 8 was 19 m s−1, while the average speed of the event reached 11.25 m s−1. The other events showed the same behavior. This is because the precipitating clouds usually move at a lower speed than the wind that pushes them [37,59].
The PyrLK model performance was evaluated by comparing the observed and estimated radar images created by the model during the same period (Table 3). According to the advection model, the rainfall event with the maximum average speed of movement (11.25 m s−1) occurred in June 2016. This speed agrees with the austral summer period and the usual presence of strong winds for the study area in this season [37,44] (Figure 4). For this event, r and the POD for the projected and observed images were the lowest, 0.62 and 0.82, respectively. In contrast, for event 10, which had the lowest speed (1.92 m s−1), r and the POD (0.79 and 0.86) were above the average of these indices (0.75 and 0.85). This indicates that the advection model is more accurate for events with slow motion of rain fields. However, for similar intermediate speeds (2.5 m s−1), lower correlation coefficients (0.68–0.75) were found. This is because PyrLK is best adjusted when the rain event characteristics meet the model assumptions between consecutive images (brightness/constant intensity and small movements in the rain cell shapes) [28,54].
As expected, the model detected fewer false alarms in events with slower motions (events 1, 9, and 10) as shown by the FAR. On the other hand, the PyrLK model generally overestimated the rainfall detection; this can be observed in the FBI indicators for all the events. This tendency was more marked for events with fast motions (events 7 and 8). Nonetheless, the average r and POD for the 11 events was 0.75 and 0.85, respectively, which, according to [58,60], statistically allows the PyrLK model to be appropriate for the study site, and consequently, the model is suited for use in the atmospheric conditions of the Tropical Andes. In addition, when using the optical flow method to estimate the rain cell motion vectors, even storms with speeds greater than 12 m s−1 can be captured by consecutive images [28].

4.2. Impact of Spatial and Temporal Resolution of the Radar Data on QPE

Table 4 summarizes the different statistical and detection indexes between hourly accumulated data for the radar–rain gauge relation of the 11 events. It can be seen that r slightly improved (from 0.67 to 0.69) when using 1-min-interval data instead of 5-min-interval data for the three spatial resolutions. The authors of [28] found similar results in their study, which analyzed the image resolution effect on the radar–rain gauge relation.
The BIAS decreased 0.1 mm on average for the three spatial resolutions when 1-min-interval data was used. This is equivalent to enhancing the estimate by 40% by using 1-min-interval data instead of 5-min-interval data, which seems to be consistent with the results (30%–60%) found by [31]. Moreover, in agreement with the previous results, the detection index improved significantly (10%) when using 1-min-interval interpolated images.
On the other hand, when the different indexes from Table 4 are compared for the three spatial resolutions, it was found that the differences are negligible. There is no considerable improvement (<5%) in the radar–rain gauge relation if the spatial resolution increases from 0.50 to 0.10 km. Consequently, 0.50-km pixel images can be used for radar–rain gauge analysis in the study area, which will result in lower computational costs for information processing.
The statistical and detection indexes improve as the radar data temporal resolution increases (Table 4 and Figure 8). The correlation coefficient decreased from 0.69 to 0.31 (45%) when the sampling time increased from 1 to 60 min. The POD showed similar behavior (decrease of 49%). The BIAS increased by 0.35 mm (73%) when the radar temporal resolution was reduced to 30 min. Similar results (20%–60%) were found by [19] with 30-min-interval radar data. The FBI index indicated that the radar data underestimation increased to 52% when the radar recorded images every 60 min in comparison to 1-min sampling, and the FAR detected by the radar increased to 50% with images recorded every 60 min.
Similar results were observed with all three pixel resolutions. However, a comparison of radar-based QPE pixels and gauge-based point observations always implies sources of uncertainty, which prevent perfect validation results: (i) The radar–rain gauge relationship worsens with lower temporal resolution, because radar images only show an integral area of the pixel, but the local (subpixel) rainfall might be higher over the gauge than integrated over the entire radar pixel [25,30,31]. Furthermore, the rain gauge network is unable to record the high spatio-temporal precipitation variability of the Tropical Andes and might miss pixel-integrated rainfall if the subpixel location of the pixel remains rain-free [4,44].
(ii) Another source of uncertainty is the horizontal displacement of falling rain observed by the radar through wind drift on its way to the rain gauge. The vertical wind profile, between the radar signal and the surface, modifies the raindrops’ precipitation trajectory [61]. If the low temporal resolution of the radar is added to this effect, the deviation between QPE and the rain gauge total increases [62,63,64]. The altitude differences between the radar beam (around 4450 m a.s.l.) and the rain gauge network are between 550 and 2000 m. The wind profile effect on the QPE–rain gauge relation has not yet been studied for the Tropical Andes’ atmospheric conditions. For this, the height and depth of the melting layer needs to be known (for instance using a vertical profile radar). Current studies using a micro rain radar (MRR) close to the study area are under development and could be used in the future.
Despite the fact that results recorded in up to 10-min-intervals are likely the most relevant, an investigation of higher integration times (10–60 min) was performed because of the possible importance for remote areas of the Tropical Andes. In this area, most rainfall monitoring networks just use manual automatic weather stations (i.e., collect rainfall data three times per day only). Only a few institutions have installed a limited number of automatic weather stations (AWS) so far. Thus, this study offers these institutions evidence of the limitations of temporally low-resolution data collection and provides insights into future AWS network planning regarding the frequency needed to properly capture the rainfall variability.

4.3. Influence of Rain Event Features

Table 5 shows the three groups of rainfall events based on rainfall motion speed.

4.3.1. Rain Field Motion

Figure 9 shows the POD and BIAS discretization of the radar–rain gauge relation as a function of the average rainfall motion (RM). There is a strong relationship between the POD and BIAS with the RM. The POD values were higher in the event of higher RM (RM > 5 m s−1). The lowest POD values occurred for RM between 1–2.5 m s−1. Although the POD differences for 1-min-interval recorded images were less than 6% between the three speed ranges, this percentage increased to 11% after 5 min. However, the difference again decreased for 60 min (5%).
The BIAS has a similar behavior to the POD regarding the RM (Figure 9). The errors were greater for lower RM (1–2.5 m s−1) and decreased for the other two ranges. The BIAS differences for events between 2.5 and 5 m s−1 and >5 m s−1 were low (11%) compared with the lower RM range (41%). These results are in agreement with those obtained by [31], who mentioned that the error in the radar–rain gauge relation decreases with the increase in the storm speed. In the same way, [28] found lower r and higher RMSE for events with lower RM.

4.3.2. Rain Total Volume and Spatial Coverage Percentage

Figure 10 shows a clear tendency of decreasing rain volumes when the radar sampling time increases. On average for all events, the rain total volume improved by 7% when using 1-min-interval images in comparison with 5-min-interval images and by 12% in comparison with 10-min-interval images. On the other hand, we calculated that the error of rain total volume using the three pixel resolutions was less than 3%. This result confirms that radar images with a 0.5-km resolution can be used for hydrological studies in this region.
There is a remarkable relationship between the spatial rainfall coverage over the basin and the rain total volume (Figure 11). This relationship is not linear, because it depends on other event characteristics, such as intensity and duration [28]. In this case an exponential relation was found (r = 0.79). Nonetheless, the coverage percentages estimated with the rain gauge network differed from those estimated by the radar data. For instance, for the event on 04/17/2015, using only the rain gauge network data, the rainfall precipitated over 80% of the sub-basin’s total area calculated with the Thiessen method and 100% of the sub-basin’s total area calculated with the IDW method (Figure 12). Nonetheless, when using the higher temporal resolution, 1-min-interval radar data, with 0.5-km pixel size, it was identified that the event covered only 58% of the basin. A possible explanation for this might be that the IDW method interpolates the rain for the whole study area, even where no rainfall is recorded. For this reason, the estimated rain volumes are greater than those calculated using the Thiessen method and those observed by the CAXX radar.
The coverage–rain volume relation significantly influences the accumulated rainfall estimation during a certain period. In this case, as in most cases of Tropical Andes studies where rain gauges are used to estimate rainfall [37,57,65,66], there is a sparse rain–gauge network, especially in the sub-basin’s headwaters. The insufficient monitoring added to the high precipitation variability present in the study area [1] can cause two scenarios: (i) rainfall underestimation when a part of the precipitation event is not captured by the rain gauge network and (ii) rainfall overestimation when the tributary area for each rain gauge is considerably high (e.g., the Thiessen method).

4.3.3. Maximum Rainfall Intensities

Figure 13 shows the relation among the maximum rainfall intensities (MRI), the rainfall motion (RM), rainfall volume, and spatial coverage percentage of the 11 events. Events with an RM greater than 5 m s−1 presented relatively low MRI for hourly data (13.46–17.10 mm h−1). The highest MRIs estimated with 5-min-interval data were recorded for the events between the RM ranges of 1–2.5 and 2.5–5 m s−1, with values of 110.3 and 103.20 mm h−1, respectively. These results can be attributed to the different types of precipitation that are generated in these ecosystems [44,46]. The heavy rain structure consists mainly of drops with mean volume diameters (Dm) greater than 1 mm [46], which in turn makes it difficult for the wind to move them easily, therefore obtaining a low average rainfall motion (1–2.5 m s−1).
In contrast, moderate (0.5 < Dm (mm) ≤ 1.0) to light rains (0.1 < Dm (mm) ≤ 0.5) can be pushed relatively easily by the wind, thus generating rain motion speeds greater than 5 m s−1. Finally, when light rain occurs in months with higher wind speeds (e.g., the event on 06/12/2016), the RM can reach values up to 15 m s−1 (Figure 6, event 8).
When the RM (2.5–5 m s−1) of event 3 (9/14/2016) of group 2 was analyzed (Figure 13), it was identified that the rainfall precipitated over the entire basin as calculated by the radar data and rain gauge method. The rainfall accumulated by the radar was higher than that estimated by the rain gauge network. This can be explained by the fact that the rain cells with heavy intensity were not captured by the rain gauges. In contrast, in event 1 of the same group (03/08/2016), where the second highest intensity was recorded (103.20 mm h−1), the rain gauge network properly captured the cell with the heaviest intensity.
On the other hand, the event with the highest MRI (110.3 mm h−1, event 3 of the RM range 1–2.5 m s−1, Figure 13) registered on 04/20/2015, presented a similar basin coverage percentage for the rain gauge network as well as the radar (using 1-min-interval data). However, the rain gauge network registered 42% more precipitation than the radar. A possible explanation for this might be that the greatest amount of rainfall was precipitated under the radar beam. The CAXX radar is installed at 4450 m a.s.l., the average height of the basin is 3454 m a.s.l., and Cuenca city, where some heavy rain events also occur, is located at 2560 m a.s.l. This suggests that the radar beam does not record rainfall events below a certain range.

5. Conclusions

In this study, the impact of spatial and temporal resolutions on the QPE accumulation under different rainfall features was analyzed. The optical flow method (PyrLK) was validated and used afterwards to determine the field vectors of the rainfall motion. Rainfall field advection was applied to enhance the radar data temporal resolution to 1 min through a liner interpolation between two 5-min-interval consecutive images. We used three different spatial resolutions and seven temporal resolutions to assess the sampling time impact on QPE under four rainfall features: horizontal rainfall motion, spatial precipitation coverage, rain total volume precipitated into the sub-basin, and rainfall intensity. The following conclusions were obtained from the study:
(a)
The PyrLK model can be used in the atmospheric conditions in the Tropical Andes. The average correlation coefficient and probability of detection index between the images recorded by the radar and estimated by the model for the 11 analyzed events were 0.75 and 0.85, respectively. However, it was observed that the PyrLK is better adjusted for events with slow rain motions (r = 0.79, POD = 0.86) than for fast motions (r = 0.62, POD = 0.82). Thus, the slower the horizontal rain motion, the higher the model performance.
(b)
The radar–rain gauge relation was enhanced when using 1-min-interval radar data instead of 5-min-interval data (r from 0.67 to 0.69) for the three pixel resolutions. The BIAS between the radar and rain gauge decreased by 40% and the POD significantly improved (11%) when using 1-min-interval interpolated images.
(c)
There was a slight improvement (<5%) in the radar–rain gauge relation when the pixel resolution is increased from 0.50 to 0.10 km. Therefore, for radar–rain gauge analysis in the study area, images of 0.50 km of spatial resolution are a good trade-off, which will reduce the computational cost of data processing.
(d)
The hourly radar–rain gauge relation decreased considerably with the increase of the data sampling time (r from 0.69 to 0.31 for recording times of 1 to 60 min). These tendencies may be explained by the fact that by decreasing the temporal resolution, the radar does not capture the rapidly evolving echoes or rain cells. In addition, the rain gauge network does not capture the high spatial variability of the Tropical Andes precipitation.
(e)
In temporal resolution terms, using 1-min-interval radar data enhanced the rain total volume estimation by 7% compared with 5-min-interval radar data and 12% compared with 10-min-interval radar data. This improvement occurred, because the 1-min-interval radar images were able to capture more details of the rain event characteristics.
(f)
The difference between the spatial rainfall coverage captured by the radar and the rain gauge network significantly influenced the calculated rain volume. This difference occurred mainly because (i) the rain gauge network is neither dense nor homogeneous, and it does not capture all of the rain cells and may (strongly) underestimate or overestimate the rainfall between stations through interpolation processes; and (ii) some precipitations may be generated below the height of the beam radar. Efforts are currently being carried out to fill this knowledge gap by using a combination of the CAXX horizontal radar and a vertically-pointed MRR to determine the vertical profile of rain and identify the height at which the precipitation events are generated.
The study gives important insights into the proper development of hydro-meteorological networks in the high Andes, which is urgently needed due to the great importance of ecosystem ecohydrological services. On the other hand, the results lead to an improvement of rainfall estimations in mountain areas in order to better determine the entire water balance in the region. In addition, the results will form the basis for implementing a precipitation nowcasting model, which, in turn, will help to enhance the rainfield flow forecasting in mountainous areas, e.g., in rainfall hazard surveillance systems.

Author Contributions

M.G. Conceptualization, Methodology, Software, Data curation, Formal Analysis, Writing—Original draft. J.O.-A. Conceptualization, Software, Writing—Review & Editing, Supervision. J.B. Writing—Review & Editing, Visualization.

Funding

The current study was developed in the context of three projects. (i) The knowledge transfer project “RadarNet-Sur” (http://www.radarnetsur.gob.ec/), administered by the University of Marburg and the Provincial Government of Loja (GPL) within the “Platform for Biodiversity and Ecosystem Monitoring and Research in South Ecuador”. This was funded by the German Research Foundation (Deutsche Forschungsgemeinschaft—DFG; BE1780/31-1 and BE1780/38-1). (ii) The project “Identificación de los procesos hidrometeorológicos que desencadenan crecidas a partir de la información suministrada por un radar de precipitación” and (iii) the project “Desarrollo de modelos para pronóstico hidrológico a partir de datos de radar meteorológico en cuencas de montaña ” These were funded by the Research Office of the University of Cuenca (DIUC) and Empresa Pública Municipal de Telecomunicaciones, Agua Potable, Alcantarillado y Saneamiento de Cuenca (ETAPA-EP).

Acknowledgments

The authors thank Rolando Célleri for helping to greatly improve the manuscript and Ing. José Medina for practical support of the computational implementation.

Conflicts of Interest

The authors declare no conflicts of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Fries, A.; Rollenbeck, R.; Bayer, F.; Gonzalez, V.; Oñate-Valivieso, F.; Peters, T.; Bendix, J. Catchment precipitation processes in the San Francisco valley in southern Ecuador: Combined approach using high-resolution radar images and in situ observations. Meteorol. Atmos. Phys. 2014, 13–29. [Google Scholar] [CrossRef]
  2. Lepioufle, J.M.; Leblois, E.; Creutin, J.D. Variography of rainfall accumulation in presence of advection. J. Hydrol. 2012, 464–465, 494–504. [Google Scholar] [CrossRef]
  3. Célleri, R.; Feyen, J. The Hydrology of Tropical Andean Ecosystems: Importance, Knowledge Status, and Perspectives. Mt. Res. Dev. 2009, 29, 350–355. [Google Scholar] [CrossRef] [Green Version]
  4. Célleri, R.; Willems, P.; Buytaert, W.; Feyen, J. Space-time rainfall variability in the Paute Basin, Ecuadorian Andes. Hydrol. Process. 2007, 21, 3316–3327. [Google Scholar] [CrossRef]
  5. Morin, E.; Goodrich, D.C.; Maddox, R.A.; Gao, X.; Gupta, H.V.; Sorooshian, S. Rainfall modeling for integrating radar information into hydrological model. Atmos. Sci. Lett. 2005, 6, 23–30. [Google Scholar] [CrossRef] [Green Version]
  6. Yoon, S.-S.; Phuong, A.T.; Bae, D.-H. Quantitative Comparison of the Spatial Distribution of Radar and Gauge Rainfall Data. J. Hydrometeorol. 2012, 13, 1939–1953. [Google Scholar] [CrossRef]
  7. Kusiak, A.; Wei, X.; Verma, A.P.; Roz, E. Modeling and prediction of rainfall using radar reflectivity data: A data-mining approach. IEEE Trans. Geosci. Remote Sens. 2013, 51, 2337–2342. [Google Scholar] [CrossRef]
  8. Li, L.; He, Z.; Chen, S.; Mai, X.; Zhang, A.; Hu, B.; Li, Z.; Tong, X. Subpixel-Based Precipitation Nowcasting with the Pyramid Lucas-Kanade Optical Flow Technique. Atmosphere 2018, 9, 260. [Google Scholar] [CrossRef]
  9. Montopoli, M.; Marzano, F.S.; Picciotti, E.; Vulpiani, G. Spatially-adaptive advection radar technique for precipitation mosaic nowcasting. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2012, 5, 874–884. [Google Scholar] [CrossRef]
  10. Páez, S. Quantitative Precipitation Forecast Using Weather Radar. Master’s Thesis, University Stuttgart, Stuttgart, Germany, 2012. [Google Scholar]
  11. Versace, P.; Sirangelo, B.; De Luca, D.L. A space-time generator for rainfall nowcasting: The PRAISEST model. Hydrol. Earth Syst. Sci. Discuss. 2008, 5, 749–783. [Google Scholar] [CrossRef]
  12. Berne, A.; Krajewski, W.F. Radar for hydrology: Unfulfilled promise or unrecognized potential? Adv. Water Resour. 2013, 51, 51–357. [Google Scholar] [CrossRef]
  13. Méndez-Antonio, B.; Caetano, E.; Soto-Cortés, G.; Rivera-Trejo, F.G.; Carvajal Rodríguez, R.A.; Watts, C. Weather Radar Data and Distributed Hydrological Modelling: An Application for Mexico Valley. Open J. Mod. Hydrol. 2013, 3, 79–88. [Google Scholar] [CrossRef]
  14. Rozalis, S.; Morin, E.; Yair, Y.; Price, C. Flash flood prediction using an uncalibrated hydrological model and radar rainfall data in a Mediterranean watershed under changing hydrological conditions. J. Hydrol. 2010, 394, 245–255. [Google Scholar] [CrossRef]
  15. Rollenbeck, R.; Bendix, J.; Fabian, P.; Boy, J.; Wilcke, W.; Dalitz, H.; Oesker, M.; Emck, P. Comparison of Different Techniques for the Measurement of Precipitation in Tropical Montane Rain Forest Regions. J. Atmos. Ocean. Technol. 2007, 24, 156–168. [Google Scholar] [CrossRef]
  16. Kanda, M. Progress in Urban Meteorology: A Review. J. Meteorol. Soc. Jpn. 2007, 85B, 363–383. [Google Scholar] [CrossRef]
  17. Bendix, J.; Fries, A.; Zárate, J.; Trachte, K.; Rollenbeck, R.; Pucha-Cofrep, F.; Paladines, R.; Palacios, I.; Orellana, J.; Oñate-Valdivieso, F.; et al. RadarNet-Sur first weather radar network in Tropical high mountains. Bull. Am. Meteorol. Soc. 2017, 98, 1235–1254. [Google Scholar] [CrossRef]
  18. Harrison, D.L.; Driscoll, S.J.; Kitchen, M. Improving precipitation estimates from weather radar using quality control and correction techniques. Meteorol. Appl. 2000, 7, 135–144. [Google Scholar] [CrossRef] [Green Version]
  19. Wilson, J.W.; Brandes, E.A. Radar Measurement of Rainfall—A Summary. Bull. Am. Meteorol. Soc. 1979, 60, 1048–1058. [Google Scholar] [CrossRef] [Green Version]
  20. Delrieu, G.; Serrar, S.; Guardo, E.; Creutin, J.D. Rain Measurement in Hilly Terrain with X-Band Weather Radar Systems: Accuracy of Path-Integrated Attenuation Estimates Derived from Mountain Returns. J. Atmos. Ocean. Technol. 1999, 16, 405–416. [Google Scholar] [CrossRef]
  21. Germann, U.; Galli, G.; Boscacci, M.; Bolliger, M. Radar precipitation measurement in a mountainous region. Q. J. R. Meteorol. Soc. 2006, 132, 1669–1692. [Google Scholar] [CrossRef] [Green Version]
  22. Rollenbeck, R.; Bendix, J. Experimental calibration of a cost-effective X-band weather radar for climate ecological studies in southern Ecuador. Atmos. Res. 2006, 79, 296–316. [Google Scholar] [CrossRef]
  23. Kim, T.J.; Kwon, H.H.; Lima, C. A Bayesian partial pooling approach to mean field bias correction of weather radar rainfall estimates: Application to Osungsan weather radar in South Korea. J. Hydrol. 2018, 565, 14–26. [Google Scholar] [CrossRef]
  24. Jordan, P.; Seed, A.; Austin, G. Sampling errors in radar estimates of rainfall. J. Geophys. Res. 2000, 105, 2247–2257. [Google Scholar] [CrossRef] [Green Version]
  25. Piccolo, F.; Chirico, G.B. Sampling errors in rainfall measurements by weather radar. Adv. Geosci. 2005, 2, 151–155. [Google Scholar] [CrossRef] [Green Version]
  26. Todini, E. A Bayesian technique for conditioning radar precipitation estimates to rain-gauge measurements. Hydrol. Earth Syst. Sci. 2010, 5, 187–199. [Google Scholar] [CrossRef]
  27. Seck, I.; Van Baelen, J. Geostatistical Merging of a Single-Polarized X-Band Weather Radar and a Sparse Rain Gauge Network over an Urban Catchment. Atmosphere 2018, 9, 496. [Google Scholar] [CrossRef]
  28. Seo, B.C.; Krajewski, W.F. Correcting temporal sampling error in radar-rainfall: Effect of advection parameters and rain storm characteristics on the correction accuracy. J. Hydrol. 2015, 531, 272–283. [Google Scholar] [CrossRef]
  29. Steiner, M.; Bell, T.L.; Zhang, Y.; Wood, E.F. Comparison of two methods for estimating the sampling-related uncertainty of satellite rainfall averages based on a large radar dataset. J. Clim. 2003, 16, 3759–3778. [Google Scholar] [CrossRef]
  30. Shucksmith, P.E.; Sutherland-Stacey, L.; Austin, G.L. The spatial and temporal sampling errors inherent in low resolution radar estimates of rainfall. Meteorol. Appl. 2011, 18, 354–360. [Google Scholar] [CrossRef]
  31. Fabry, F.; Bellon, A.; Duncan, M.R.; Austin, G.L. High resolution rainfall measurements by radar for very small basins: The sampling problem reexamined. J. Hydrol. 1994, 161, 415–428. [Google Scholar] [CrossRef]
  32. Villarini, G.; Mandapaka, P.V.; Krajewski, W.F.; Moore, R.J. Rainfall and sampling uncertainties: A rain gauge perspective. J. Geophys. Res. Atmos. 2008, 113, 1–12. [Google Scholar] [CrossRef]
  33. Rodriguez-Iturbe, I. Ecohydrology: A hydrologic perspective of climate-soil-vegetation dynamies. Water Resour. Res. 2000, 36, 3–9. [Google Scholar] [CrossRef] [Green Version]
  34. Liu, C.; Krajewski, W.F. A comparison of methods for calculation of Radar-Rainfall Hourly Accumulations. J. Am. Water Resour. Assoc. 1996, 32. [Google Scholar] [CrossRef]
  35. Villarini, G.; Krajewski, W.F. Review of the different sources of uncertainty in single polarization radar-based estimates of rainfall. Surv. Geophys. 2010, 31, 107–129. [Google Scholar] [CrossRef]
  36. Rossa, A.M.; Laudanna Del Guerra, F.; Borga, M.; Zanon, F.; Settin, T.; Leuenberger, D. Radar-driven high-resolution hydro-meteorological forecasts of the 26 September 2007 Venice flash flood. J. Hydrol. 2010, 394, 230–244. [Google Scholar] [CrossRef]
  37. Campozano, L.; Célleri, R.; Trachte, K.; Bendix, J.; Samaniego, E. Rainfall and Cloud Dynamics in the Andes: A Southern Ecuador Case Study. Adv. Meteorol. 2016, 2016. [Google Scholar] [CrossRef]
  38. Farley, K.A.; Anderson, W.G.; Bremer, L.L.; Harden, C.P. Compensation for ecosystem services: An evaluation of efforts to achieve conservation and development in Ecuadorian páramo grasslands. Environ. Conserv. 2011, 38, 393–405. [Google Scholar] [CrossRef]
  39. Mosquera, G.M.; Lazo, P.X.; Célleri, R.; Wilcox, B.P.; Crespo, P. Runoff from Tropical alpine grasslands increases with areal extent of wetlands. Catena 2015, 125, 120–128. [Google Scholar] [CrossRef]
  40. Pesántez, J.D.; Martinez, A. Propuesta de Modelo de Gestión de la Subcuenca del río Tomebamba, Como Herramienta de Manejo Integrado y Conservación. Master’s Thesis, Universidad del Azuay, Cuenca, Ecuador, 2015; pp. 1–183. [Google Scholar]
  41. Vuille, M.; Bradley, R.S.; Keimig, F. Climate Variability in the Andes of Ecuador and Its Relation to Tropical Pacific and Atlantic Sea Surface Temperature Anomalies. J. Clim. 2000, 13, 2520–2535. [Google Scholar] [CrossRef] [Green Version]
  42. Bendix, J. Heavy rainfall episodes in Ecuador during El Niño events and associated regional atmospheric circulation and SST patterns. Adv. Geosci. 2006, 6, 43–49. [Google Scholar] [CrossRef] [Green Version]
  43. Tapia, J.D.; Célleri, R. Estructura de la Precipitación y su Variabilidad Espacio-Temporal a lo Largo de una Gradiente Altitudinal. Bachelor’s Thesis, Universidad de Cuenca, Cuenca, Ecuador, 2016; pp. 1–71. [Google Scholar]
  44. Padrón, R.S.; Wilcox, B.P.; Crespo, P.; Célleri, R. Rainfall in the Andean Páramo: New Insights from High-Resolution Monitoring in Southern Ecuador. J. Hydrometeorol. 2015, 16, 985–996. [Google Scholar] [CrossRef]
  45. Carrillo-Rojas, G.; Silva, B.; Córdova, M.; Célleri, R.; Bendix, J. Dynamic Mapping of Evapotranspiration Using an Energy Balance-Based Model over an Andean Páramo Catchment of Southern Ecuador. Remote Sens. 2016, 8, 160. [Google Scholar] [CrossRef]
  46. 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]
  47. Heistermann, M.; Jacobi, S.; Pfaff, T. Technical Note: An open source library for processing weather radar data (wradlib). Hydrol. Earth Syst. Sci. 2013, 17, 863–871. [Google Scholar] [CrossRef]
  48. Pfaff, T.; Heistermann, M.; Jacobi, S. wradlib—An Open Source Library for Weather Radar Data Processing. In Proceedings of the The Seventh European Conference on Radar in Meteorology and Hydrology, Toulouse, France, 24–29 June 2012. [Google Scholar]
  49. NOAA. Re-analysis data of wind speed and direction. In National Centers for Environmental Prediction NCEP-FNL; NOAA: Silver Spring, MA, USA, 2000. [Google Scholar]
  50. Rinehart, R.E.; Garvey, E.T. Three-dimensional storm motion detection by conventional weather radar. Nature 1978, 273, 287–289. [Google Scholar] [CrossRef]
  51. Horn, B.; Schunck, B. Determining optical flow. Artif. Intell. 1981, 17, 185–203. [Google Scholar] [CrossRef] [Green Version]
  52. Pfaff, T. Processing and Analysis of Weather Radar Data for Use in Hydrology. Ph.D. Thesis, University of Stuttgart, Stuttgart, Germany, 2013. [Google Scholar]
  53. Woo, W.; Wong, W. Operational Application of Optical Flow Techniques to Radar-Based Rainfall Nowcasting. Atmosphere 2017, 8, 48. [Google Scholar] [CrossRef]
  54. Bouguet, J.-Y. Pyramidal Implementation of the Affine Lucas Kanade Feature Tracker Description of the algorithm. Intel Corp. Microprocess. Res. Labs 2001, 5, 4. [Google Scholar]
  55. Lucas, B.D.; Kanade, T. An Iterative Image Registration Technique with an Application to Stereo Vision. Imaging 1981, 130, 674–679. [Google Scholar]
  56. Bradski, G.; Kaehler, A. Learning OpenCV; Loukides, M., Ed.; O’Reilly Media, Inc.: Sebastopol, CA, USA, 2008; ISBN 9780596516130. [Google Scholar]
  57. Ballari, D.; Castro, E.; Campozano, L. Validation of Satellite Precipitation (Trmm 3B43) in Ecuadorian Coastal Plains, Andean Highlands and Amazonian Rainforest. ISPRS Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci. 2016, XLI-B8, 305–311. [Google Scholar] [CrossRef]
  58. Jatho, N.; Pluntke, T.; Kurbjuhn, C.; Bernhofer, C. An approach to combine radar and gauge based rainfall data under consideration of their qualities in low mountain ranges of Saxony. Nat. Hazards Earth Syst. Sci. 2010, 10, 429–446. [Google Scholar] [CrossRef] [Green Version]
  59. Nugent, A.D.; Smith, R.B.; Minder, J.R. Wind Speed Control of Tropical Orographic Convection. J. Atmos. Sci. 2014, 71, 2695–2712. [Google Scholar] [CrossRef]
  60. Nerini, D.; Zulkafli, Z.; Wang, L.-P.; Onof, C.; Buytaert, W.; Lavado, W.; Guyot, J.-L. A comparative analysis of TRMM-rain gauge data merging techniques at the daily time scale for distributed rainfall-runoff modelling applications. J. Hydrometeorol. 2015. [Google Scholar] [CrossRef]
  61. Sandford, C. Correcting for wind drift in high resolution radar rainfall products: A feasibility study. J. Hydrol. 2015, 531, 284–295. [Google Scholar] [CrossRef]
  62. Collier, C.G. The Impact of Wind Drift on the Utility of Very High Spatial Resolution Radar Data Over Urban Areas. Phys. Chem. Earth. 1999, 24, 889–893. [Google Scholar] [CrossRef]
  63. Harrold, T.W.; English, E.J.; Nicholass, C.A. The accuracy of radar-derived rainfall measurements in hilly terrain. Q. J. R. Meteorol. Soc. 1974, 100, 331–350. [Google Scholar] [CrossRef]
  64. Lauri, T.; Koistinen, J.; Moisseev, D. Advection-Based Adjustment of Radar Measurements. Mon. Weather Rev. 2012, 140, 1014–1022. [Google Scholar] [CrossRef]
  65. Ochoa-Tocachi, B.F.; Buytaert, W.; De Bièvre, B.; Célleri, R.; Crespo, P.; Villacís, M.; Llerena, C.A.; Acosta, L.; Villazón, M.; Guallpa, M.; et al. Impacts of land use on the hydrological response of Tropical Andean catchments. Hydrol. Process. 2016, 30. [Google Scholar] [CrossRef]
  66. Zulkafli, Z.; Buytaert, W.; Onof, C.; Manz, B.; Tarnavsky, E.; Lavado, W.; Guyot, J.-L. A Comparative Performance Analysis of TRMM 3B42 (TMPA) Versions 6 and 7 for Hydrological Applications over Andean-Amazon River Basins. J. Hydrometeorol. 2014, 15, 581–592. [Google Scholar] [CrossRef]
Figure 1. Study Area: Paute River basin, CAXX 100-km radar range, rain gauges, and climate stations. Lower left inset: radar network RadarNet-Sur.
Figure 1. Study Area: Paute River basin, CAXX 100-km radar range, rain gauges, and climate stations. Lower left inset: radar network RadarNet-Sur.
Water 11 01038 g001
Figure 2. X-band radar station: Radar CAXX Rainscanner-RS120 on Paragüillas peak (4450 m a.s.l).
Figure 2. X-band radar station: Radar CAXX Rainscanner-RS120 on Paragüillas peak (4450 m a.s.l).
Water 11 01038 g002
Figure 3. Rain distribution and total accumulation (mm per event) for the time period shown in Table 2 (duration) within the Cuenca sub-basin for the 11 events. Black dots represent the rain gauge network. The maps were created using the adjusted radar data.
Figure 3. Rain distribution and total accumulation (mm per event) for the time period shown in Table 2 (duration) within the Cuenca sub-basin for the 11 events. Black dots represent the rain gauge network. The maps were created using the adjusted radar data.
Water 11 01038 g003
Figure 4. Wind speed and direction at 4450 m a.s.l. estimated with the National Centers for Environmental Prediction (NCEP)-FNL data. (a) Relationship of the wind speed on the X-axis (positive values towards the east and negative values towards the west) with the wind direction (90° indicates the direction from the east to the west). (b) Monthly average distribution of wind speed and direction.
Figure 4. Wind speed and direction at 4450 m a.s.l. estimated with the National Centers for Environmental Prediction (NCEP)-FNL data. (a) Relationship of the wind speed on the X-axis (positive values towards the east and negative values towards the west) with the wind direction (90° indicates the direction from the east to the west). (b) Monthly average distribution of wind speed and direction.
Water 11 01038 g004
Figure 5. Workflow of data processing.
Figure 5. Workflow of data processing.
Water 11 01038 g005
Figure 6. Distribution of (a) the speed and (b) the direction of the 11 rain events, calculated with the optical flow method using intensity gradient information on several image scales (PyrLK) advection model. Event direction means the angle of the origin of the rain cell (e.g., 90° = origin from the east).
Figure 6. Distribution of (a) the speed and (b) the direction of the 11 rain events, calculated with the optical flow method using intensity gradient information on several image scales (PyrLK) advection model. Event direction means the angle of the origin of the rain cell (e.g., 90° = origin from the east).
Water 11 01038 g006
Figure 7. (a) Distribution of the wind speed and direction of the re-analysis data in contrast to the rainfall motion generated by the optical flow method for all precipitation events (negative event speed indicates that precipitation is moving west). (b) Hourly speed correlation generated by re-analysis data (X-axis) and the advection model (Y-axis) for the 11 events.
Figure 7. (a) Distribution of the wind speed and direction of the re-analysis data in contrast to the rainfall motion generated by the optical flow method for all precipitation events (negative event speed indicates that precipitation is moving west). (b) Hourly speed correlation generated by re-analysis data (X-axis) and the advection model (Y-axis) for the 11 events.
Water 11 01038 g007
Figure 8. Variation of the Pearson correlation coefficient (a) and the probability of detection index (b) according to the sampling time and pixel resolution.
Figure 8. Variation of the Pearson correlation coefficient (a) and the probability of detection index (b) according to the sampling time and pixel resolution.
Water 11 01038 g008
Figure 9. Probability of detection and BIAS indexes for precipitation events according to average rainfall motion. Subplots for different temporal and spatial resolutions of radar data are shown.
Figure 9. Probability of detection and BIAS indexes for precipitation events according to average rainfall motion. Subplots for different temporal and spatial resolutions of radar data are shown.
Water 11 01038 g009
Figure 10. Accumulated rainfall volume (m3) for the 11 events according to the different radar temporal and spatial resolutions.
Figure 10. Accumulated rainfall volume (m3) for the 11 events according to the different radar temporal and spatial resolutions.
Water 11 01038 g010
Figure 11. Relation of water volume and its percentage of coverage over the study basin area.
Figure 11. Relation of water volume and its percentage of coverage over the study basin area.
Water 11 01038 g011
Figure 12. Accumulated rainfall (mm) for the event on 04/17/2015. Rain spatial distribution for the basin using only rain gauge network data (upper left) and radar data for different sampling times with 0.5-km spatial resolution.
Figure 12. Accumulated rainfall (mm) for the event on 04/17/2015. Rain spatial distribution for the basin using only rain gauge network data (upper left) and radar data for different sampling times with 0.5-km spatial resolution.
Water 11 01038 g012
Figure 13. Event discretization according to the horizontal rainfall motion for the seven temporal resolutions and for 0.5-km spatial resolution (the same results were found when using the two other spatial resolutions). The precipitated water volume, the area percentage covered by rainfall, and the maximum intensities recorded for 5 min and 60 min are shown.
Figure 13. Event discretization according to the horizontal rainfall motion for the seven temporal resolutions and for 0.5-km spatial resolution (the same results were found when using the two other spatial resolutions). The precipitated water volume, the area percentage covered by rainfall, and the maximum intensities recorded for 5 min and 60 min are shown.
Water 11 01038 g013
Table 1. Characteristics of the middle and upper parts of the Paute River basin.
Table 1. Characteristics of the middle and upper parts of the Paute River basin.
AreaPerimeterVegetal CoverShapeAltitude MinAltitude Max
km2km%-m a.s.l.m a.s.l.
1604.5271PR (71.3); AMF (10.6)
P (7.79); UR (5.7) BQ (2.5)
CO25404680
RainRunoffEToQ. min. dailyQ. máx. dailyR.C.
mm year−1mm year−1mm year−1mmmm-
12216185800.251.070.51
CO: Circular to Oval; R.C.: Runoff Coefficient; PR: Páramo; AMF: Andean Montane Forest; QF: Quinoa Forest; P: Pasture; UR: Urban.
Table 2. List of 11 selected rain events that occurred between 2015 and 2016.
Table 2. List of 11 selected rain events that occurred between 2015 and 2016.
ID EventInitial Date Final DateDurationMin. Volume Rain GaugeMax. Volume Rain GaugeMax. Intensity (5 min)Max. Intensity (h)
mm/dd/yy h:mmmm/dd/yy h:mmhmmmmmm h−1mm h−1
104/17/2015 00:0004/17/2015 03:003:001.1032.1062.4017.10
204/20/2015 16:0004/20/2015 20:004:001.2747.00110.4032.20
304/21/2015 14:0004/21/2015 18:004:001.6027.2040.8017.10
405/03/2015 15:0005/04/2015 01:0010:006.3050.2976.2024.89
503/08/2016 14:0003/08/2016 17:003:001.5028.96103.2023.00
603/10/2016 00:0003/10/2016 06:006:001.9030.4834.8013.20
704/27/2016 17:0004/27/2016 22:005:001.2022.7063.6014.50
806/12/2016 12:0006/12/2016 16:004:001.2020.0739.6213.46
909/14/2016 16:0009/14/2016 22:006:002.0338.3074.4029.10
1010/24/2016 15:0010/24/2016 18:003:001.6025.3057.5920.20
1111/10/2016 17:0011/10/2016 00:007:002.9034.3073.2025.30
Table 3. Correlation coefficient r and average detection index for each rain event between images estimated by the advection model and observed by radar every 5 min.
Table 3. Correlation coefficient r and average detection index for each rain event between images estimated by the advection model and observed by radar every 5 min.
ID Eventr *Detection Index
PODFARFBI
10.830.900.191.11
20.700.820.351.38
30.700.840.291.19
40.740.880.271.24
50.680.780.421.45
60.780.860.261.20
70.750.870.401.61
80.620.820.311.22
90.790.910.181.14
100.790.860.201.10
110.750.890.251.25
* The p-values of r were lower than 0.05 for the 11 rainfall events.
Table 4. Statistical error metrics and detection indexes between radar data and rain gauges for different spatial and temporal resolutions of radar images.
Table 4. Statistical error metrics and detection indexes between radar data and rain gauges for different spatial and temporal resolutions of radar images.
Sample Time
(min)
0.50 km
r *BIAS
(mm)
RMSE
(mm)
PODFARFBI
10.69−0.152.460.760.130.98
50.68−0.212.460.710.150.89
100.65−0.282.420.680.160.84
150.6−0.342.450.660.160.78
200.53−0.42.50.620.190.73
300.44−0.472.540.560.20.66
600.34−0.722.850.420.230.48
0.25 km
10.67−0.112.490.790.121.03
50.67−0.232.460.710.150.87
100.65−0.292.420.680.150.83
150.56−0.362.450.650.160.78
200.5−0.392.490.620.180.72
300.4−0.472.530.560.190.65
600.3−0.752.890.410.240.47
0.10 km
10.68−0.132.480.80.121.08
50.67−0.262.460.710.150.86
100.65−0.32.410.680.150.83
150.58−0.362.450.650.160.77
200.51−0.42.480.610.180.72
300.4−0.492.570.550.180.64
600.31−0.742.840.410.240.47
* The p-values for r were lower than 0.05.
Table 5. Precipitation event groups according to rainfall motion.
Table 5. Precipitation event groups according to rainfall motion.
Rain Motion Range (m s−1)DateRain Motion (m s−1)
StartEndVxVyV
1–2.510/24/2016 15:0010/24/2016 18:000.90−0.641.92
04/20/2015 16:0004/20/2015 20:00−0.64−1.772.41
04/17/2015 00:0004/17/2015 03:00−1.88−0.592.43
03/10/2016 00:0003/10/2016 06:00−1.48−1.512.45
05/03/2015 15:0005/04/2015 01:00−2.010.362.46
2.5–509/14/2016 16:0009/14/2016 22:00−1.30−0.702.64
03/08/2016 14:0003/08/2016 17:002.221.453.54
11/10/2016 17:0011/10/2016 00:00−2.991.713.71
>504/27/2016 17:0004/27/2016 22:00−4.97−0.125.18
04/21/2015 14:0004/21/2015 18:00−5.110.665.42
06/12/2016 12:0006/12/2016 16:00−11.040.9711.24

Share and Cite

MDPI and ACS Style

Guallpa, M.; Orellana-Alvear, J.; Bendix, J. Tropical Andes Radar Precipitation Estimates Need High Temporal and Moderate Spatial Resolution. Water 2019, 11, 1038. https://doi.org/10.3390/w11051038

AMA Style

Guallpa M, Orellana-Alvear J, Bendix J. Tropical Andes Radar Precipitation Estimates Need High Temporal and Moderate Spatial Resolution. Water. 2019; 11(5):1038. https://doi.org/10.3390/w11051038

Chicago/Turabian Style

Guallpa, Mario, Johanna Orellana-Alvear, and Jörg Bendix. 2019. "Tropical Andes Radar Precipitation Estimates Need High Temporal and Moderate Spatial Resolution" Water 11, no. 5: 1038. https://doi.org/10.3390/w11051038

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