Next Article in Journal
HRP-OG: Online Learning with Generative Feature Replay for Hypertension Risk Prediction in a Nonstationary Environment
Previous Article in Journal
The Influence of Induced Head Acceleration on Lower-Extremity Biomechanics during a Cutting Task
Previous Article in Special Issue
Rainfall Observation Leveraging Raindrop Sounds Acquired Using Waterproof Enclosure: Exploring Optimal Length of Sounds for Frequency Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

High-Resolution Rainfall Estimation Using Ensemble Learning Techniques and Multisensor Data Integration

Department of Physics, FMIPA Universitas Indonesia, Depok 16424, Indonesia
*
Author to whom correspondence should be addressed.
Sensors 2024, 24(15), 5030; https://doi.org/10.3390/s24155030
Submission received: 26 June 2024 / Revised: 26 July 2024 / Accepted: 1 August 2024 / Published: 3 August 2024
(This article belongs to the Special Issue Atmospheric Precipitation Sensors)

Abstract

:
In Indonesia, the monitoring of rainfall requires an estimation system with a high resolution and wide spatial coverage because of the complexities of the rainfall patterns. This study built a rainfall estimation model for Indonesia through the integration of data from various instruments, namely, rain gauges, weather radars, and weather satellites. An ensemble learning technique, specifically, extreme gradient boosting (XGBoost), was applied to overcome the sparse data due to the limited number of rain gauge points, limited weather radar coverage, and imbalanced rain data. The model includes bias correction of the satellite data to increase the estimation accuracy. In addition, the data from several weather radars installed in Indonesia were also combined. This research handled rainfall estimates in various rain patterns in Indonesia, such as seasonal, equatorial, and local patterns, with a high temporal resolution, close to real time. The validation was carried out at six points, namely, Bandar Lampung, Banjarmasin, Pontianak, Deli Serdang, Gorontalo, and Biak. The research results show good estimation accuracy, with respective values of 0.89, 0.91, 0.89, 0.9, 0.92, and 0.9, and root mean square error (RMSE) values of 2.75 mm/h, 2.57 mm/h, 3.08 mm/h, 2.64 mm/h, 1.85 mm/h, and 2.48 mm/h. Our research highlights the potential of this model to accurately capture diverse rainfall patterns in Indonesia at high spatial and temporal scales.

1. Introduction

Indonesia is a tropical country with various complexities in its rainfall patterns [1,2,3,4]. Based on data processing and analysis over 30 years (1991–2020), it has been identified that the rainfall distribution patterns in the Indonesian region consist of 487 monsoonal, 178 equatorial, and 34 local types. This has a direct impact on many aspects of life, such as transportation, agriculture, water resource management, and disasters [5,6,7,8]. Moreover, extreme rainfall can cause significant impacts on the economic sector, infrastructure, and public safety [9,10]. However, at present, the information regarding rainfall in Indonesia has regrettably not reached a high resolution, and the geographical coverage is limited.
Rainfall measurements in Indonesia are generally performed using two methods: direct observation via a rain gauge and indirect estimation via remote sensing [11,12,13]. The measurements using rain gauges are accurate; however, the total surface area of rain gauges in Indonesia covers only 1.29 × 10−11 of the country’s total area. This inadequate gauge density results in inaccuracies when representing the total rainfall across an area [14]. The distribution of rain gauges in Indonesia also highlights a concerning disparity in gauge density. Papua Island, the eastern region of Indonesia, has fewer rain gauges even though the region has local rain patterns with its diverse topography. In contrast, Java Island, the southern region of Indonesia, has a monsoon pattern and non-complex topography, but the rain gauge distribution there is very dense. The optimum rain gauge network has been the subject of research and operation in Indonesia over the years [15,16]. However, the national weather radar network does not cover the entire observation area in Indonesia [17]. Several weather radars also cannot explain a single relationship between radar reflectivity and rainfall [18]. The variability in this relationship can certainly affect the accuracy of rainfall estimates using weather radars. Meanwhile, weather satellite images with global coverage have resolution limitations [19]. In measuring rainfall in Indonesia, the necessity for a wide spatial coverage with a high level of resolution is the main challenge [20,21]. This condition explains why the existence of only one data source is insufficient to fulfill Indonesia’s need for reliable rainfall data.
Indonesia’s national rainfall information product is provided by the Meteorology, Climatology, and Geophysics Agency (BMKG) via its official website at www.bmkg.go.id, accessed on 31 July 2024. The information displayed is limited to rain classification images covering Indonesia with a time resolution of 1 h and an accumulation of 24 h. Apart from that, a research institution in Indonesia, the National Innovation Research Agency (BRIN), has the Sadewa website at the address https://sadewa.brin.go.id/sadewabgr/ (accessed on 24 June 2024), which provides information on rainfall predictions with a time resolution of 1 hours resulting from modeling with a spatial resolution of 1 km. The results of previous studies also have limitations. The generated products are typically not in the form of numerical data but are mostly in the form of classifications [22,23,24], are not yet able to produce comprehensive spatial products, tend to be limited to data per point or region [25,26,27,28,29], and do not provide actual and up-to-date information to meet real-time needs [20,30,31,32,33]. Based on these facts, there is no numerical, spatial, high-resolution rainfall estimation model that uses observational equipment. In fact, information about rainfall in Indonesia is really needed because it has a direct impact on many aspects of life [5,6,7,8], such as agriculture, because rainfall affects agricultural productivity [34]; transportation, because it can be the cause of transportation accidents [35,36,37,38]; and disasters, because its impact is often associated with disasters such as floods and landslides [9,10,39,40].
Currently, the existing rainfall estimation models cannot detect hidden patterns or non-linear trends in rainfall data, which are important features for producing accurate rainfall information products [41]. The dynamics of the air in the atmosphere may be significantly non-linear; however, there are still discernable patterns in its movement. These simple patterns may possibly be manually completed; however, for big data and complex non-linear patterns, the patterns may be generated using machine learning [42]. The implementation of machine learning in the field of meteorology includes the use of algorithms capable of processing extensive data, such as observation network data, weather radars, satellites, and weather models [43,44,45]. Through processing more significant data and more complex analysis, machine learning can be developed into a more accurate rainfall estimation model [46,47]. Several approaches to rainfall estimation have been explored by different researchers, and the importance of a high spatial and temporal resolution for accurate estimation precision has been emphasized [48,49,50,51,52,53,54]. However, using machine learning requires managing several challenging issues, namely, the unbalancing of classes [55], such as a disproportionate number of instances without rain or deficient rainfall compared with instances with high rainfall, a large number of missing attributes [56] during the process in which the sensor collects data, and the need to work incrementally immediately after new data are available. Addressing these challenges is crucial for harnessing the full potential of machine learning, especially in rainfall estimation.
Tree-based machine learning is a popular approach in the world of machine learning. The development of tree-based algorithms in machine learning has experienced a significant evolution. The developments began with decision trees, which are a simple tree structure used for decision making [57]. The drawbacks of decision trees include their tendency to overfit data and their typically lower accuracy [58,59], often necessitating ensemble methods to tackle these issues. The solution to the overfitting problem in decision trees is “ensembling”, which involves assembling many weak trees together into a strong forest, making predictions from each of the trees, and voting to decide on the winning prediction. This is further developed into a random forest, which uses a large number of trees to increase the accuracy and overcome overfitting [60]. With the bagging technique, weak learners/decision trees are arranged in parallel, and so they do not learn from each other; therefore, they have limitations in terms of computation [61]. In 1996, Schapire developed a boosting technique called AdaBoost, an ensemble learning technique that gives more weight to data misclassified by the previous model [62]. In many applications, AdaBoost is sensitive to both noise and data outliers and is slow in processing because it involves all data samples in each iteration [63]. Furthermore, gradient boosting advances development through the introduction of the concept of gradients and iteratively updating the model to reduce prediction errors [64]. Finally, extreme gradient boosting (XGBoost) improves boosting by increasing efficiency and accuracy with regularization, overfitting management, and layered trees [65]. The strength of XGBoost is its scalability in all scenarios. The scalability of XGBoost is due to several important systems and algorithmic optimizations. One is a new tree-learning algorithm for handling sparse data [65]. Moreover, parallel and distributed computing makes the learning process faster [66]. XGBoost is now the preferred method for developing predictive models due to its remarkable accuracy, efficiency, and adaptability [67,68]. Recently, XGBoost has even dominated the applied machine learning domain and has won several Kaggle competitions [69].
Applying XGBoost in rainfall estimation and related meteorological phenomena has demonstrated its efficacy and accuracy in predictive modeling. XGBoost was utilized alongside rain gauges, radars, and satellite data, achieving a high correlation and low RMSE, highlighting its robustness in integrating diverse data sources [70]. Despite its effectiveness, that study focused on one city, which may not be generalized well to different geographic locations with varying climate conditions. Coupling XGBoost with the Bat algorithm to estimate evapotranspiration showed superior performance in arid regions compared with other models [71]. Additionally, utilizing XGBoost for precipitation nowcasting outperformed other methods, proving its superiority in real-time applications [72]. Finally, integrating XGBoost with another method further enhanced precipitation nowcasting, demonstrating significant improvements in detection probability and error reduction [73]. The existing literature and studies show that XGBoost can be computationally demanding, especially when dealing with large datasets and deep trees. Training such models requires significant computational resources, which might not be available on less capable hardware. This complexity extends to hyperparameter tuning, which necessitates extensive experimentation to achieve optimal performance [66,74].
Currently, the application of machine learning to estimate rainfall in Indonesia is increasing and developing [75]. This study aimed to estimate rainfall in Indonesia, which exhibits seasonal, equatorial, and local characteristics. Our proposed approach integrates data from different instrument sources (i.e., rain gauges, weather radars, and weather satellites) to produce high-resolution rainfall estimates. Furthermore, ensemble learning, XGBoost, was applied to anticipate the sparse data that may occur due to the limited number of rain gauge points, limited radar coverage, and rain data imbalance problems. Finally, a spatial rainfall estimation model was built with a temporal resolution close to real time. Therefore, it is hoped that this study can benefit all stakeholders who use rainfall information in Indonesia, especially the transportation sector, which needs actual and high-resolution rainfall information.

2. Materials and Methods

The data sources used in this research consisted of satellite data, weather radars, and rain gauges provided by the BMKG. The weather radar and rain gauge data were obtained directly from the BMKG monitoring system, which operates in various regions of Indonesia. A complex method for rainfall estimation is proposed in this study. This study included a bias correction process for satellite data to ensure data accuracy and weather radar data integration and implemented the XGBoost machine learning ensemble to process big data. The research stage diagram can be seen in Figure 1.

2.1. Study Area

Indonesia is an archipelagic country situated in Southeast Asia, straddling the equator. Indonesia is located between 6° N to 11° S and 95° E to 141° E. Geographically, Indonesia consists of more than 17,000 large and small islands, making it the largest archipelagic country in the world. Overall, the land area is 1,993,662,036 km2; meanwhile, the water area reaches 6,653,341,439 km2 [76]. These geographic characteristics provide Indonesia with a diverse landscape, including mountains, tropical rainforests, and long beaches. More than 54,700 km of coastline connects the land with the sea, one of the main factors that influences the rainfall patterns in various regions of Indonesia. In this study, understanding the geography of Indonesia is fundamental to the process of rainfall estimation. Climatologically, Indonesia’s territory consists of seasonal zones and non-seasonal zones [77]. Indonesia’s territory is also divided into three types of rain patterns, namely, monsoonal, local, and equatorial types [78], as shown in Figure 2.

2.2. Data

In this rainfall estimation study, information from two different technologies, weather radar data and data from the Himawari satellite, became a training feature to obtain rainfall estimates. The main target of this estimation was global precipitation measurement (GPM) satellite data, which have been previously corrected using rain gauge data. This target was selected to ensure that the rainfall estimation results followed the results of the field observation measurements. The data period used covered the entire year of 2022. The training period spanned from January to November to maximize the training process. December was used for validation because it represents Indonesia’s peak rainy season, providing a range of conditions, from no rain to maximum rainfall.
Various satellite-based rainfall estimation products with a high spatial and temporal resolution have been developed to meet hydrometeorological data needs [79]. To date, the GPM-integrated multi-satellite retrievals (IMERG) product has shown more consistent performance and can be a good alternative for rainfall estimation [80,81]. The IMERG system operates in near-real time, providing data in two runs, known as early and late IMERG; subsequently, after the monthly gauge analysis is received, the final IMERG data are generated [82]. Early-run products have the potential for real-time applications with shorter delay times [83]. GPM combines the data from active and passive instruments within the GPM constellation to produce rainfall estimates, referred to as IMERG [84]. The GPM dataset is an early-run product produced every 30 min, covering 1 January to 31 December 2022. Spatial data with a pixel resolution of 10 km, covering a wide area of Indonesia, were used in this research.
A weather radar is an instrument that can detect various atmospheric parameters, such as rainfall, cloud movement, and wind speed, as well as the physical properties of rain or ice grains [85]. A weather radar transmits electromagnetic waves into the atmosphere, detects their reflections hitting objects, and measures the energy in radar reflectivity [86]. In this research, the resulting weather radar reflectivity data were used as a feature to estimate rainfall. In total, thirty-five weather radar units of various types were used to support this research. Of this number, six weather radar locations were used for data training, namely, the Lampung, Banjarmasin, Deli Serdang, Pontianak, Gorontalo, and Biak weather radars. The six radar locations used represent the existing rain patterns in Indonesia, which are seasonal, equatorial, and local rain patterns. The data are presented in a spatial resolution of 500 square meters. The temporal resolution of the weather radar data was updated every 10 min. The combination of a high spatial resolution and fast temporal resolution enabled very accurate and near-real-time rainfall estimates.
Another feature we used was data from the Himawari satellite in the form of brightness temperature. The Himawari satellite measures electromagnetic radiation emitted by various surfaces and the atmosphere in various spectral channels. These data produce brightness temperatures in various spectral channels, measured in kelvin (K). Equipped with an advanced Himawari imager (AHI) sensor, the Himawari satellite can produce various types of images. The Himawari satellite has 16 bands consisting of 3 visible bands, 3 near-infrared (NIR) bands, and 10 infrared (IR) or thermal bands. Each channel has a different wavelength, resolution, and use [87,88]. This research used IR band 13 with a wavelength of 10.4 μm. The selection of the IR channel was designed for use for observations at night since it can detect thermal radiation from objects on the Earth’s surface even though there is no direct sunlight available. In addition, the wavelength in band 13 is a clean longwave window with relatively high transparency toward electromagnetic radiation and is not significantly affected by water vapor in the atmosphere. This channel helps detect surfaces and clouds [88,89]. The spatial resolution of the Himawari satellite’s IR sensor is 2 km, indicating the sensor’s ability to distinguish objects at a minimum distance of that size.
The rainfall data were obtained from the Automated Weather Observing System (AWOS) equipment, which is an airport weather observation system that can provide weather information for 10 parameters, namely, wind speed, visibility, weather, sky conditions, temperature, relative humidity, wind chill, the heat index, pressure, and rainfall [90,91]. A rain gauge, which is one of the sensors in the AWOS equipment, is capable of measuring rain to an accuracy of 0.1 mm, so it is susceptible to detecting rainfall. With a time resolution of 1 min, this sensor can provide very detailed data about rainfall patterns in short intervals. With six rain gauge points, their use was advantageous in the GPM data correction process, which involved identifying rain events. Using the high-resolution rain gauges in the AWOS equipment increased the accuracy and precision of the rainfall data used in this research.
The instruments used in this study, including information on time resolution, spatial resolution, working principles, and units of product produced, are presented in Table 1.

2.3. Data Preprocessing

The data and instruments used in this study had different resolutions and units of measurement. The rainfall amount acquired from the rain gauges (mm), which were used to verify the GPM satellite rainfall product, was converted into the rain rate (mm/h) by taking into account the duration of the rainfall that occurred [92]. Rainfall data measured every 1 min by a rain gauge were converted into the rainfall intensity in units of millimeters per hour (mm/h). This conversion process was performed by measuring the total amount of rainfall every 10 min and then calculating the average rainfall intensity. In this way, rain gauge rainfall data measured every 1 min could be converted into the rainfall intensity displayed every 10 min in millimeters per hour.
The spatial resolution between the GPM satellite, Himawari satellite, and weather radars was another element that must be equalized. The weather radar pixels and Himawari satellites were adjusted to the GPM satellite pixels. We used bilinear interpolation techniques for the inter-pixel adjustments. In this technique, the new pixel value was determined based on the average by giving weight to the closest pixels. This technique was used because it can produce continuous data between known data points, making it easier to maintain spatial consistency. The bilinear interpolation process was carried out alternately on one of the vertical or horizontal sides [93]. Interpolation only involves the four nearest neighboring points, so the interpolation process is fast even when using a large dataset. By using the average value of the four nearest neighboring points, bilinear interpolation minimizes the loss of information that may occur [94].
The frequency of re-recording of all research instruments was equalized over a certain period. The temporal resolution must have the same interval [95]. The rain gauge data recorded every 1 min were downscaled to 10 min intervals to match the temporal resolution of the weather radars and Himawari satellites. The method used was to take a 10 min average of 10 consecutive data points. The data from the weather radars and Himawari satellites recorded every 10 min were maintained because they were in accordance with the desired resolution. For the training process, data that overlapped with the GPM data were used, and they were recorded every 30 min. The running process still used data every 10 min. Data from the GPM satellite recorded every 30 min were used for the training process.
Some of the data-filtering techniques applied in this research were resampling techniques for time series and intersection data. The resampling technique changed the time–frequency of data from one interval to another. With different sources and different time–frequencies, resampling helped make the data consistent and allowed for easier comparisons. Meanwhile, intersection was used to obtain the same elements between two or more data sets [96]. Intersection performed the operational function of merging several intersections of data sets that were previously incomplete to become sequentially complete. The research files were processed in the Zarr format, which stores each piece of data separately using data compression techniques to help reduce storage space requirements and speed up data transfer [97].

2.4. Bias Correction Strategy

The GPM product had several uncertainty sources, such as sensor calibration, retrieval errors, and orographic effects [98,99]; therefore, the product needed to be corrected for actual rainfall data from the rain gauges [100,101]. The bias correction method modified the bias correction strategy used in [102]. The first step in this research stage was to identify the rain events. In this context, only ‘hit’ events underwent bias correction; meanwhile, the other conditions were ignored in the analysis. A ‘hit’ event was when the GPM correctly detected rain. For the ‘hit’ events that were identified, corrections were made to reduce the bias between the GPM data and actual rain events that occurred using the linear regression method, which was previously performed by the authors of [103,104,105]. As previously explained in the data subsection, the six rain gauge points that were used for GPM correction are listed in Table 2. The bias correction could be carried out with at least one rain gauge available in each area to be corrected [106].

2.5. Integration of Weather Radar Data

A composite map can show the spatial rainfall distribution over a larger area, which is very important for understanding the current weather conditions and possibly predicting the rainfall cloud movements [17]. The number of weather radars scattered in the research area allowed for the intersection of coverage areas. To generate a composite map from several weather radars, the merging of the CMAX reflectivity values in each intersection coverage area was performed [107]. The complete information on the location, frequency band, and weather radar polarization used in this study is shown in Table 3.

2.6. Ensemble Learning Approach

Ensemble learning techniques have achieved state-of-the-art performance in diverse machine learning applications through the combination of predictions from two or more base models. We used the XGBoost algorithm for this research, more specifically, the XGBoost Python Package. XGBoost is a decision tree-based optimization technique that builds on the gradient descent method. The gradient descent method is used to optimize the loss function; additionally, regularization parameters are employed to prevent overfitting. The fundamental concept underlying the XGBoost algorithm is the minimization of the following objective function, which consists of the loss function and regularization terms [65]. Additionally, XGBoost is efficient in handling large datasets [108]. First, the procedure used in this algorithm involves dividing the original dataset into multiple sub-datasets. Then, each subset is randomly assigned to the base learner for prediction. The algorithm calculates the result of the weak learner based on a certain weight. Finally, the model results can be expressed as the weighted sum of the predicted results of all the decision trees.
Hyperparameters were used to enhance the algorithm’s performance results, significantly affecting various model tests. In this study, the hyperparameter method that was applied was Bayesian optimization. This method uses a probabilistic model to direct the search for the best hyperparameters. The previous search results are used to estimate the combination of hyperparameters that is most likely to provide better performance. The Bayesian optimization process consists of two main components: a surrogate model, which fits each observation to a target function, and an acquisition function, which balances exploitation and exploration to determine the next evaluation point. Bayesian optimization stabilizes exploration and exploitation to find optimal options while avoiding missing out on better configurations in areas that have not been evaluated [109,110,111].
The parameters that were optimized and expected to improve the model’s performance included seven parameters [110,112,113], as shown in Table 4 below.
For several reasons, Bayesian optimization has shown significant advantages in hyperparameter tuning for XGBoost compared with other methods. Primarily, it offers a systematic approach to exploring the hyperparameter space efficiently, utilizing prior information to guide the search process and update beliefs about the best regions to explore as more data are gathered. This probabilistic framework allows for a more informed and directed search, often resulting in faster convergence to optimal hyperparameters. For instance, studies such as that by Qiu et al. (2022) demonstrated that Bayesian optimization could achieve lower prediction errors and higher model performance with fewer iterations than traditional grid or random search methods. Moreover, Bayesian optimization effectively balances exploration and exploitation, enhancing the model’s accuracy and robustness [109,110,111,112,113,114].

2.7. Evaluation of Estimation Results

The estimated rainfall results from the ensemble learning model were evaluated using statistical and detection indicators [115]. The statistical indicators assessed the root-mean-square error (RMSE) values of the estimated rainfall; meanwhile, the detection indicators gauged the accuracy of the rainfall detection. The RMSE was used to measure the average magnitude of the errors between the estimated values and the actual values. The RMSE calculated the square root of the average of the squared differences between the predicted values and the actual values. The RMSE value provides an overall picture of the magnitude of the errors and can be useful in comparing model performance. In this case, accuracy refers to the extent to which the model could correctly classify rain occurrences or non-rain occurrences.

3. Results and Discussion

3.1. Data Correlation

The first task was to display and discuss the relationship between the input and target data that were used to build the rainfall estimation model. Specifically, we examined the relationship between the rainfall data produced by the GPM satellite and the rainfall data measured using the rain gauges.
The rainfall estimates produced by the GPM satellite with rain-gauge-measured rainfall data revealed a low correlation [116,117], as shown in Figure 3. This may be the result of several errors in the data retrieval algorithms and sensors [118]. Other studies have stated that this discrepancy can be caused by cloud characteristics, climate, season, geographic location, and topography [119,120]. However, seeing the consistently low correlation in several locations, we suspect that the differences occurred due to lag time [121], i.e., the time delay between raindrops from clouds and rainfall on the Earth’s surface. This time delay means that when the GPM satellite observed rain from clouds, there was a certain time required for the raindrops to reach the Earth’s surface. Many factors, such as the distance between the cloud and the surface, can cause this delay. Therefore, the rainfall estimation data needed to be evaluated and corrected for bias.
The second data correlation occurred between the rainfall data produced by the GPM satellite and the brightness temperature data measured by the Himawari satellite. The plotting of the two datasets is shown in Figure 4. Plotting the data at six points shows an inverse relationship between the cloud brightness temperature from the Himawari satellite data and the rainfall from the GPM satellite data, with a correlation value between −0.37 and −0.49. The correlation slope, even though insignificant, illustrates that when the cloud brightness temperature was low, the rainfall was high. Conversely, when the cloud brightness temperature was high, the rainfall was low. When the cloud brightness temperature decreases or becomes colder, this indicates that the cloud is thicker and tends to have larger water particles. The correlation value between the GPM satellite and the Himawari satellite was low due to differences in the working principles regarding the measurement targets. Working at a wavelength from 0.46 to 13.3 μm [122], the Himawari satellite can identify cloud droplets with a diameter of < 50 μm [123]. However, the presence of droplets in clouds does not always indicate rain. The droplets can turn into raindrops when small droplets collide with each other and merge through a process called coalition and coalescence [124]. However, not all droplets undergo this process until they reach a size large enough to fall as rain. In this case, the Himawari satellite is too sensitive for the task of identifying raindrops. Meanwhile, the GPM satellite has technical specifications that allow it to detect precipitation with a minimum size of 0.2 mm [125]. Therefore, the GPM satellite tends to be more accurate because it focuses on raindrops as the measurement object. Physically, this difference in measurement principles can result in a low correlation between the two satellites. The GPM satellite essentially measures raindrops, while the Himawari satellite measures cloud droplets, which do not necessarily turn into raindrops.
There was a further correlation between the GPM satellite data and weather radar data. In physical terms, the higher the radar reflectivity (the dBZ value), the greater the possibility of high-intensity rain in that area. However, as presented in Figure 5, the data show that the correlation of the rainfall data produced by the GPM satellite with the radar reflectivity data was low. At a shorter measurement time resolution of 30 min, the correlation of the GPM satellite with the weather radars was lower than the daily or monthly measurement times [126]. This was due to the limitations of the early-run IMERG in capturing rapid changes in rain patterns. Topography was also a general factor that influenced the correlation of the GPM satellite with the weather radars [127]. The fundamental issue was the difference in the working principles of these two instruments. GPM identified precipitation directly via a combination of dual-frequency precipitation radar technology and the GPM microwave imager [128]. Meanwhile, the radars measured the reflectivity that was reflected back from the collection of droplets detected in the radar volume, but not from each droplet individually. High reflectivity can be caused by a few large droplets or many small droplets; therefore, it is difficult to determine the distribution and size categories of droplets based solely on reflectivity values [129,130,131]. In physical terms, the correlation between the GPM satellite data and the weather radars was low because GPM measured raindrops, while the radars measured a collection of droplets that may not necessarily be raindrops.

3.2. Bias Correction Result

Compared with the original product, the bias-corrected GPM product showed a better performance in estimating rainfall values. We display a graph of the probability density function (PDF), which is used to describe the probability distribution of a continuous random variable at a certain point. Additionally, a cumulative distribution function (CDF) graph is presented to illustrate the cumulative PDF of a random variable. This graph shows how often a random variable takes either a certain value or less than a certain value in its distribution; as a result, the graph can help in understanding the distribution, evaluating the probabilities, and identifying the patterns in GPM data before and after bias correction.
The PDF and CDF graphs show a significant gap between the uncorrected GPM graph and the measured rainfall graph. This gap indicates that the initial satellite rainfall estimates had a significant bias; therefore, bias correction needed to be applied to bring the estimated results closer to the actual observation data. The GPM, which was not corrected, appeared to underestimate the measured rainfall. Meanwhile, the corrected GPM graph coincides more closely with the measured rainfall graph. This indicates that the applied bias correction succeeded in accurately approximating the observational data. The data distribution is shown in Figure 6.

3.3. Weather Radar Network

A national weather radar network was generated in this study with a total of 35 weather radars being integrated. This merging employed the spatial averaging of composite data, accommodating differences in frequency bands, such as C and X; differences in polarization, including single and dual; and differences in range, as illustrated in Figure 7.

3.4. Hyperparameter Tuning Results

Hyperparameter tuning optimization produced optimal values for each parameter using the Bayesian optimization method, as shown in Table 5.
Each hyperparameter was optimized to achieve a low objective value, which indicated a better model performance, as seen in Figure 8. The objective function in XGBoost regression is an evaluation metric, such as the mean squared error (MSE). In the learning_rate graph, it can be seen that the optimal value was in the range of around 0.01 to 0.1. A learning rate that is too high (>0.1) tends to produce a higher objective value, indicating overfitting. Meanwhile, a learning rate that is too low (<0.01) also does not provide the best performance because the model may be too slow to converge. For the hyperparameter max_depth, the optimal value was one. A too-deep depth (>6) is not optimal because it makes the model too complex. The graph for n_estimators shows that underestimating the number of estimators gave the best objective value. There is an ideal number of decision trees that should be used to achieve a balance between bias and variance. In the subsample graph, the optimal value was 0.96. For min_child_weight, the optimal value was 0.14. Higher values (>5) did not provide the best performance, indicating that the model requires some complexity to capture variations in the data. The gamma graph shows that the optimal value was 0.08, which shows that proper regularization is necessary to avoid overfitting. Finally, for colsample_bytree, the optimal value was around 0.4. Low or high values did not perform best, indicating that proper column sampling is necessary to reduce the correlation between trees.

3.5. Product and Evaluation

After undergoing preprocessing and several processing steps, including weather radar integration, the GPM satellite bias correction, and the modeling of rainfall estimates using the XGBoost algorithm, the rainfall estimation results for the Indonesian region were obtained. Figure 9 presents the spatial rainfall estimation products for the Indonesian region. This rainfall estimate was generated in near-real time every 10 min.
We evaluated the rainfall estimation product to demonstrate the RMSE value and accuracy. Table 6 shows the results of evaluating the rainfall estimation model in this study. Furthermore, we captured and analyzed several rain events in areas with different rain patterns and variations in rainfall intensity.
Based on Figure 10, the weather radar image shows a high reflectivity (>30 dBZ), and the Himawari image shows thick clouds with a low cloud top temperature, recorded at 190 °K or −83.5 °C. This condition indicates the potential for very heavy rain in the area. At the start of the very heavy rain, the rain gauge measured 24.2 mm/h, and the model estimated rainfall to be 20.5 mm/h, while the GPM Satellite only detected normal rain, namely, 6.7 mm/h. After 1 h, all three detected very heavy rain: the rain gauge measured 25 mm/h, the rainfall model estimate was 21.0 mm/h, and the GPM satellite measurement was 20.1 mm/h. This fact shows that the estimation model is more sensitive to local variations and can quickly capture weather dynamics.
Figure 11 shows the Himawari satellite image features with cloud cover, which indicates the potential for heavy rainfall. The brightness temperature of the Himawari satellite was 210 °K, reflecting the presence of a lot of water vapor in the clouds. Likewise, in the radar measurements, the reflectivity of 30 dBZ indicated the number of droplets that would eventually become rain. The results of the rainfall estimation model tended to be close to the results of the rain gauge measurements; this shows its reliability. In contrast, GPM consistently recorded lower values, which may indicate a delay or underestimation.
The consistency of the model estimation results and GPM underestimation of measured rainfall was still visible in moderate rain events that occurred in Gorontalo. Figure 12 shows the rainfall estimates with images that are not too thick, indicating a relatively low rainfall intensity. The model estimated rainfall from 11:00 to 12:30 with a moderate intensity, proven to be accurate based on field observations from the rain gauge, which recorded rainfall of 1.0 mm at 11:30, increasing to 7.6 mm at 12:00, and decreasing again to 3.8 mm at 12:30. These data show that the rain only lasted a short time, around one hour, with a normal intensity, based on the estimates obtained from the rainfall model. This validation indicates that the rainfall estimation model can provide accurate rainfall projections, especially when identifying short rainfall durations.
The rainfall estimates for light-rain events also showed better results when compared with the GPM rain estimates. As Figure 13 shows, the rainfall estimates based on the input Himawari data with a brightness temperature of 260 °K and radar data with a reflectivity of 15 dBZ produced a rainfall estimate of 3.9 mm/h, close to the rain gauge measurement of 4.2 mm/h. Physically, rain can, indeed, occur when the cloud top temperature approaches 260 °K [132]. However, at this temperature, the number of droplets is not very significant, as measured by the radars. Therefore, light rain occurred in a relatively short time. These estimation results reflect the model’s ability to utilize input data from radars and Himawari satellites to produce rainfall estimates that are close to observed values, including during light-rain events.
Very-light-rain conditions can be seen in Figure 14. The Himawari satellite data showed a brightness temperature that was quite cold, namely, 260 °K, but the radar measurements showed no reflectivity. However, very light rainfall was detected by the rain gauge. This was also detected by the model and GPM. Physically, this condition shows a lot of water vapor in the atmosphere, but very little is condensed into droplets because more is evaporating than condensing, and this condition is not measurable by radars. This could be due to the lack of aerosols in the atmosphere, which act as condensation nuclei [124]. These few droplets still become raindrops through collision and coalescence processes and fall to the Earth’s surface, even though they only become very light rain [133,134].
Based on the rainfall estimation images and the evaluation results of the rainfall estimation model, the XGBoost ensemble learning model in this research can be effectively applied to produce very accurate spatial rainfall estimates with a temporal resolution of every 10 min. Even though the training data correlation was not very good, and the data were not balanced, the XGBoost model was able to produce a low RMSE with reasonable accuracy. This demonstrates that the XGBoost model was proven to be accurate in estimating rainfall, even when it used training data, some of which had zero or no value. This is a result of the splitting process employed by XGBoost when building a decision tree during the training process. XGBoost treats non-valuable data as separate elements in the tree. If a feature has no value, the algorithm separates the data into the following two groups: one group with no value and the other group with a known value. This is the reason XGBoost does not lose information from worthless data. In fact, XGBoost is able to select the split point with the highest gain. In other words, the algorithm chooses the best way to separate data that have no value from those that have value so that the gain or profit (decrease in the loss function) is maximized. This technique is referred to as “sparsity-aware split finding” [65,135].
XGBoost can process data in parallel to handling large-scale spatial data, as well as meeting short-time resolution requirements [136]. This is known as parallel learning because XGBoost relies on an ensemble of decision trees. Each tree is built independently of the others, meaning that the building process of each tree does not depend on the results of the other trees. The data can be divided into smaller subsets, which are then processed separately. Additionally, XGBoost applies independent regularization to each tree. Therefore, each tree can be constructed separately without requiring information from the other trees. The characteristic of the independent decision tree in the XGBoost model is that it is stump-shaped, namely, it is a tree with only one level or one branch. The fact that the tree has a limited depth helps prevent overfitting and makes the tree more general and helpful in understanding the linear relationships between features and targets.

4. Conclusions

This study produced a model and a rainfall estimation product that use various data from different instruments, including rain gauges, weather radars, and weather satellites. In the process of model development, data preprocessing was conducted, which involved techniques such as resampling, intersection, and bilinear interpolation. Additionally, bias correction techniques were applied to the GPM satellite data, using rain gauge data as a reference. Furthermore, several weather radars were combined into an integrated weather radar network. The incorporation of the XGBoost model played a crucial role in ensuring the accuracy of rainfall estimates. Based on the discussion and analysis of the results, the rainfall estimates in this research can be applied throughout Indonesia, including areas with monsoonal, local, and equatorial rainfall patterns with near-real-time resolution, providing the latest information every 10 min. This product has the potential to accurately capture diverse rainfall patterns in Indonesia at high spatial and temporal scales. Future work will focus on estimation techniques capable of strengthening confidence levels. It not only produces a single estimation value but can also provide an overview of the uncertainty of the estimate.

Author Contributions

Conceptualization, M.P. and M.S.R.; methodology, M.P.; software, M.P.; validation, M.S.R. and D.H.; formal analysis, M.P., M.S.R. and D.H.; investigation, M.P., M.S.R. and D.H.; resources, M.P.; data curation, M.P.; writing—original draft preparation, M.P.; writing—review and editing, M.P., M.S.R. and D.H.; visualization, M.P.; supervision, M.S.R. and D.H. All authors have read and agreed to the published version of the manuscript.

Funding

The third author (DH) is funded by the Faculty of Mathematics and Natural Sciences, Universitas Indonesia under Publication Grant scheme 2024.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data are available upon request due to restrictions.

Acknowledgments

The authors are grateful to the Meteorology, Climatology, and Geophysics Agency (BMKG) for providing the data and facilities for this research.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Belgaman, H.A.; Ichiyanagi, K.; Suwarman, R.; Tanoue, M.; Aldrian, E.; Utami, A.I.; Kusumaningtyas, S.D. Characteristics of seasonal precipitation isotope variability in Indonesia. Hydrol. Res. Lett. 2017, 11, 92–98. [Google Scholar] [CrossRef]
  2. Hendon, H.H. Indonesian Rainfall Variability: Impacts of ENSO and Local Air-Sea Interaction. Am. Meteorol. Soc. 2003, 16, 1775–1790. [Google Scholar] [CrossRef]
  3. Pramuwardani, I.; Hartono; Sunarto; Sopaheluwakan, A. Indonesian rainfall variability during Western North Pacific and Australian monsoon phase related to convectively coupled equatorial waves. Arab. J. Geosci. 2018, 11, 673. [Google Scholar] [CrossRef]
  4. Marzuki, M.; Hashiguchi, H.; Yamamoto, M.K.; Mori, S.; Yamanaka, M.D. Regional variability of raindrop size distribution over Indonesia. Ann. Geophys. 2013, 31, 1941–1948. [Google Scholar] [CrossRef]
  5. Narulita, I.; Ningrum, W. Extreme flood event analysis in Indonesia based on rainfall intensity and recharge capacity. IOP Conf. Ser. Earth Environ. Sci. 2018, 118, 012045. [Google Scholar] [CrossRef]
  6. Ikhwali, M.F.; Nur, S.; Darmansyah, D.; Hamdan, A.M.; Ersa, N.S.; Aida, N.; Yusra, A.; Satria, A. A review of climate change studies on paddy agriculture in Indonesia. IOP Conf. Ser. Earth Environ. Sci. 2022, 1116, 012052. [Google Scholar] [CrossRef]
  7. Kirono, D.G.C.; Tapper, N.J. Enso rainfall variability and impacts on crop production in Indonesia. Phys. Geogr. 1999, 20, 508–519. [Google Scholar] [CrossRef]
  8. Pambudi, A.S. Overview and evaluation of the Indonesia’s water resources management policies for food security. Indones. J. Appl. Environ. Stud. 2021, 2, 84–93. [Google Scholar] [CrossRef]
  9. Adibroto, T.A.; Wijayanti, P.; Adhi, R.P.; Nugroho, R. Preliminary study on socio-economic aspect towards Jakarta climate resilient (case study: Cengkareng District, West Jakarta and Penjaringan District, North Jakarta). IOP Conf. Ser. Earth Environ. Sci. 2019, 314, 012035. [Google Scholar] [CrossRef]
  10. Chayyani, N.R.; Gravitiani, E.; Suryanto. Socio-economic vulnerability and losses of flood in Lampung, Indonesia. IOP Conf. Ser. Earth Environ. Sci. 2020, 423, 012054. [Google Scholar] [CrossRef]
  11. Tondang, Y.M.; Situmorang, M.; Ikhsan, T.; Darmawan, Y. Accuracy of Weather Radar Products for Rainfall Estimation in North Sumatra Region. J. Pengkaj. Ilmu Pembelajaran Mat. IPA IKIP Mataram 2023, 11, 351–357. [Google Scholar] [CrossRef]
  12. Rahmawati, N. Space-time variogram for daily rainfall estimates using rain gauges and satellite data in mountainous tropical Island of Bali, Indonesia (Preliminary Study). J. Hydrol. 2020, 590, 125177. [Google Scholar] [CrossRef]
  13. Ramadhan, R.; Marzuki, M.; Yusnaini, H.; Muharsyah, R.; Suryanto, W.; Sholihun, S.; Vonnisa, M.; Battaglia, A.; Hashiguchi, H. Capability of GPM IMERG Products for Extreme Precipitation Analysis over the Indonesian Maritime Continent. Remote Sens. 2022, 14, 412. [Google Scholar] [CrossRef]
  14. Mishra, A.K. Effect of Rain Gauge Density over the Accuracy of Rainfall: A Case Study over Bangalore, India. SpringerPlus 2013, 2, 311. [Google Scholar] [CrossRef] [PubMed]
  15. Savitri, E.; Wahyuningrum, N.; Nugroho, H.Y.S.H.; Sudiana, N.; Indrajaya, Y.; Christanto, N. Identification of the rain gauge stations for the participatory flood and landslide mitigation in the Serayu river basin, Central Java. IOP Conf. Ser. Earth Environ. Sci. 2022, 1109, 012007. [Google Scholar] [CrossRef]
  16. Sunusi, N. Comparison of some schemes for determining the optimal number of rain gauges in a specific area: A case study in an urban area of South Sulawesi, Indonesia. AIMS Environ. Sci. 2022, 9, 244–260. [Google Scholar] [CrossRef]
  17. Permana, D.S.; Hutapea, T.D.; Praja, A.S.; Paski, J.A.I.; Makmur, E.E.S.; Haryoko, U.; Umam, I.H.; Saepudin, M.; Adriyanto, R. The Indonesia In-House Radar Integration System (InaRAISE) of Indonesian Agency for Meteorology Climatology and Geophysics (BMKG): Development, Constraint, and Progress. IOP Conf. Ser. Earth Environ. Sci. 2019, 303, 012051. [Google Scholar] [CrossRef]
  18. Hashiguchi, H.; Vonnisa, M.; Nugroho, S.; Yoseva, M. Z-R Relationships for Weather Radar in Indonesia from the Particle Size and Velocity (Parsivel) Optical Disdrometer. In Proceedings of the 2018 Progress in Electromagnetics Research Symposium (PIERS-Toyama), Toyama, Japan, 1–4 August 2018; pp. 37–41. [Google Scholar] [CrossRef]
  19. Maier, M.W.; Gallagher, F.W.; Germain, K.S.; Anthes, R.; Zuffada, C.; Menzies, R.; Piepmeier, J.; Di Pietro, D.; Coakley, M.M.; Adams, E. Architecting the Future of Weather Satellites. Bull. Am. Meteorol. Soc. 2021, 102, E589–E610. [Google Scholar] [CrossRef]
  20. Fatkhuroyan, F.; Wati, T.; Sukmana, A.; Kurniawan, R. Validation of Satellite Daily Rainfall Estimates Over Indonesia. Forum Geogr. 2018, 32, 170–180. [Google Scholar] [CrossRef]
  21. Arno, G.; Muflihah, M.; Mujahidin, M. Uncertainty of Optimal Rain Gauge Coastal Region: Case Study Makassar. J. Civ. Eng. Forum 2020, 7, 121–132. [Google Scholar] [CrossRef]
  22. Gunadi, I.G.A.; Gunawan, I.M.A.O.; Candana, P.E.W.H.; Arnawa, I.A.W.; Putra, K.A.E.K. Klasifikasi curah hujan harian menggunakan learning vector quantization. J. Ilmu Komput. Indones. (JIK) 2022, 7. Available online: http://www.ogimet.com (accessed on 21 June 2024).
  23. Yasper, A.; Handoko, D.; Putra, M.; Aliwarga, H.K.; Rosid, M.S.R. Hyperparameters Optimization in XGBoost Model for Rainfall Estimation: A Case Study in Pontianak City. J. Penelit. Pendidik. IPA 2023, 9, 7113–7121. [Google Scholar] [CrossRef]
  24. Abdilah, M.N.; Ruhiat, Y.; Guntara, Y. Rainfall Classification Analysis Using Naïve Bayes Classifier Based on Air And Wind Temperatures in Serang City. SPEKTRA J. Fis. Dan Apl. 2024, 9, 39–47. [Google Scholar] [CrossRef]
  25. Ananda, N.; Mukhlish, F.; Wicaksana, H.S.; Budiawan, I. Spatial Evaluation of Estimation Rainfall on Weather Radar Using Marshall-Palmer Z-R Method in West Java. J. Otomasi Kontrol dan Instrumentasi 2024, 16, 35–43. [Google Scholar] [CrossRef]
  26. Hutapea, T.D.F.; Permana, D.S.; Praja, A.S.; Muzayanah, L.F. Modification of Z-R Relationship Constants in Surabaya Radar for Improving The Accuracy of Rainfall Estimates. J. Meteorol. Geofis. 2021, 21, 91–97. [Google Scholar] [CrossRef]
  27. Mulsandi, A.; Mamenun, M.; Fitriano, L.; Hidayat, R. Perbaikan Estimasi Curah Hujan Berbasis Data Dengan Memperhitungkan Faktor Pertumbuhan Awan. J. Sains Teknol. Modif. Cuaca 2019, 20, 67–78. [Google Scholar] [CrossRef]
  28. Nurasniyati, M.; Adriat, R. Estimasi Curah Hujan Di Kota Pontianak Berdasarkan Suhu, Ketebalan dan Tekanan Puncak Awan. Prisma Fisika 2018, 6, 184–189. [Google Scholar]
  29. Oktaviani, R. Estimation of Rainfall Rate Cumulative Distribution in Indonesia Using Global Satellite Mapping of Precipitation Data. KnE Eng. 2019, 1, 259–265. [Google Scholar] [CrossRef]
  30. Supari, S.; Cho, J. Performa Produk Curah Hujan Harian Dari Data Satelit Dan Reanalisis Di Indonesia. 2021. Available online: https://www.researchgate.net/publication/353446008 (accessed on 21 June 2024).
  31. Ayasha, N. A Comparison of Rainfall Estimation Using Himawari-8 Satellite Data In Different Indonesian Topographies. Int. J. Remote Sens. Earth Sci. (IJReSES) 2021, 17, 189–200. [Google Scholar] [CrossRef]
  32. Alfuadi, N.; Wandala, A. Comparative Test of Several Rainfall Estimation Methods Using Himawari-8 Data. Int. J. Remote Sens. Earth Sci. (IJReSES) 2016, 13, 95–104. [Google Scholar] [CrossRef]
  33. Setiyoko, A.; Osawa, T.; Nuarsa, W. Evaluation of GSMaP Precipitation Estimates Over Indonesia. Int. J. Environ. Geosci. 2019, 3, 26–43. [Google Scholar]
  34. Torres, M.; Howitt, R.; Rodrigues, L. Analyzing Rainfall Effects on Agricultural Income: Why Timing Matters. EconomiA 2019, 20, 1–14. [Google Scholar] [CrossRef]
  35. Cao, Y.; Wu, Z.; Xu, Z. Effects of Rainfall on Aircraft Aerodynamics. Prog. Aerosp. Sci. 2014, 71, 85–127. [Google Scholar] [CrossRef]
  36. Perera, L.P.; Soares, C.G. Weather Routing and Safe Ship Handling in the Future of Shipping. Ocean Eng. 2017, 130, 684–695. [Google Scholar] [CrossRef]
  37. Sangkharat, K.; Thornes, J.E.; Wachiradilok, P.; Pope, F.D. Determination of the Impact of Rainfall on Road Accidents in Thailand. Heliyon 2021, 7, e06061. [Google Scholar] [CrossRef] [PubMed]
  38. Zhou, X.; Cheng, L.; Li, M. Assessing and Mapping Maritime Transportation Risk Based on Spatial Fuzzy Multi-Criteria Decision Making: A Case Study in the South China Sea. Ocean Eng. 2020, 208, 107403. [Google Scholar] [CrossRef]
  39. Dalagnol, R.; Gramcianinov, C.B.; Crespo, N.M.; Luiz, R.; Chiquetto, J.B.; Marques, M.T.A.; Neto, G.D.; de Abreu, R.C.; Li, S.; Lott, F.C.; et al. Extreme rainfall and its impacts in the Brazilian Minas Gerais state in January 2020: Can we blame climate change? Clim. Resil. Sustain. 2022, 1, e15. [Google Scholar] [CrossRef]
  40. Ying, L.K. Rainfall Characteristics and Its Effect on Road Infrastructure Health. Int. J. Integr. Eng. 2019, 11, 234–246. [Google Scholar]
  41. Kundu, S.; Biswas, S.K.; Tripathi, D.; Karmakar, R.; Majumdar, S.; Mandal, S. A review on rainfall forecasting using ensemble learning techniques. e-Prime 2023, 6, 100296. [Google Scholar] [CrossRef]
  42. Chen, H.; Chandrasekar, V.; Tan, H.; Cifelli, R. Rainfall Estimation From Ground Radar and TRMM Precipitation Radar Using Hybrid Deep Neural Networks. Geophys. Res. Lett. 2019, 46, 10669–10678. [Google Scholar] [CrossRef]
  43. Yu, J.; Li, X.F.; Lewis, E.; Blenkinsop, S.; Fowler, H.J. UKGrsHP: A UK high-resolution gauge–radar–satellite merged hourly precipitation analysis dataset. Clim. Dyn. 2020, 54, 2919–2940. [Google Scholar] [CrossRef] [PubMed]
  44. Folino, G.; Guarascio, M.; Chiaravalloti, F.; Gabriele, S. A Deep Learning based architecture for rainfall estimation integrating heterogeneous data sources. In Proceedings of the 2019 International Joint Conference on Neural Networks (IJCNN), Budapest, Hungary, 14–19 July 2019; pp. 1–8. [Google Scholar]
  45. Guarascio, M.; Folino, G.; Chiaravalloti, F.; Gabriele, S.; Procopio, A.; Sabatino, P. A Machine Learning Approach for Rainfall Estimation Integrating Heterogeneous Data Sources. IEEE Trans. Geosci. Remote Sens. 2022, 60, 1–11. [Google Scholar] [CrossRef]
  46. Latif, S.D.; Hazrin, N.A.B.; Koo, C.H.; Ng, J.L.; Chaplot, B.; Huang, Y.F.; El-Shafie, A.; Ahmed, A.N. Assessing rainfall prediction models: Exploring the advantages of machine learning and remote sensing approaches. Alex. Eng. J. 2023, 82, 16–25. [Google Scholar] [CrossRef]
  47. Sayed, B.T.; Al-Mohair, H.K.; Alkhayyat, A.; Ramírez-Coronel, A.A.; Elsahabi, M. Comparing machine-learning-based black box techniques and white box models to predict rainfall-runoff in a northern area of Iraq, the Little Khabur River. Water Sci. Technol. 2023, 87, 812–822. [Google Scholar] [CrossRef] [PubMed]
  48. Tian, W.; Yi, L.; Liu, W.; Huang, W.; Ma, G.; Zhang, Y. Ground radar precipitation estimation with deep learning approaches in meteorological private cloud. J. Cloud Comput. 2020, 9, 22. [Google Scholar] [CrossRef]
  49. Khan, R.S.; Bhuiyan, A.E. Artificial Intelligence-Based Techniques for Rainfall Estimation Integrating Multisource Precipitation Datasets. Atmosphere 2021, 12, 1239. [Google Scholar] [CrossRef]
  50. Liu, Z.; Yang, Q.; Shao, J.; Wang, G.; Liu, H.; Tang, X.; Xue, Y.; Bai, L. Improving daily precipitation estimation in the data scarce area by merging rain gauge and TRMM data with a transfer learning framework. J. Hydrol. 2022, 613, 128455. [Google Scholar] [CrossRef]
  51. Lei, H.; Zhao, H.; Ao, T. A two-step merging strategy for incorporating multi-source precipitation products and gauge observations using machine learning classification and regression over China. Hydrol. Earth Syst. Sci. Discuss. 2022, 26, 2969–2995. [Google Scholar] [CrossRef]
  52. Mohia, Y.; Absi, R.; Lazri, M.; Labadi, K.; Ouallouche, F.; Ameur, S. Quantitative Estimation of Rainfall from Remote Sensing Data Using Machine Learning Regression Models. Hydrology 2023, 10, 52. [Google Scholar] [CrossRef]
  53. Pan, Y.; Yuan, Q.; Ma, J.; Wang, L. Improved Daily Spatial Precipitation Estimation by Merging Multi-Source Precipitation Data Based on the Geographically Weighted Regression Method: A Case Study of Taihu Lake Basin, China. Int. J. Environ. Res. Public Health 2022, 19, 13866. [Google Scholar] [CrossRef] [PubMed]
  54. Wehbe, Y.; Temimi, M.; Adler, R.F. Enhancing precipitation estimates through the fusion of weather radar, satellite retrievals, and surface parameters. Remote Sens. 2020, 12, 1342. [Google Scholar] [CrossRef]
  55. Narwane, S.; Sawarkar, S. Machine Learning and Class Imbalance: A Literature Survey. Ind. Eng. J. 2019, 12, 1–8. [Google Scholar] [CrossRef]
  56. Emmanuel, T.; Maupong, T.; Mpoeleng, D.; Semong, T.; Mphago, B.; Tabona, O. A survey on missing data in machine learning. J. Big Data 2021, 8, 140. [Google Scholar] [CrossRef] [PubMed]
  57. Quinlan, J.R. Induction of Decision Trees. Mach. Learn. 1986, 1, 81–106. [Google Scholar] [CrossRef]
  58. Prajwala, T.R. A Comparative Study on Decision Tree and Random Forest Using R Tool. Int. J. Adv. Res. Comput. Commun. Eng. 2015, 4, 196–199. [Google Scholar] [CrossRef]
  59. Leiva, R.G.; Anta, A.F.; Mancuso, V.; Casari, P. A Novel Hyperparameter-free Approach to Decision Tree Construction that Avoids Overfitting by Design. IEEE Access 2019, 7, 99978–99987. [Google Scholar] [CrossRef]
  60. Breiman, L. Random Forests. Mach. Learn. 2001, 45, 5–32. [Google Scholar] [CrossRef]
  61. Manzali, Y.; Elfar, M. Random Forest Pruning Techniques: A Recent Review. Oper. Res. Forum 2023, 4, 43. [Google Scholar] [CrossRef]
  62. Schapire, R.E. Explaining AdaBoost. In Empirical Inference; Springer: Berlin/Heidelberg, Germany, 1995; pp. 81–106. [Google Scholar]
  63. Cao, J.; Kwong, S.; Wang, R. A noise-detection based AdaBoost algorithm for mislabeled data. Pattern Recognit. 2012, 45, 4451–4465. [Google Scholar] [CrossRef]
  64. Friedman, J.H. Greedy Function Approximation: A Gradient Boosting Machine. Ann. Stat. 2001, 29, 1189–1232. [Google Scholar] [CrossRef]
  65. Chen, T.; Guestrin, C. XGBoost: A Scalable Tree Boosting System. In Proceedings of the KDD ’16: 22nd ACM SIGKDD International Conference on Knowledge Discovery and Data Mining, San Francisco, CA, USA, 13–17 August 2016. [Google Scholar] [CrossRef]
  66. Bentéjac, C.; Csörgő, A.; Martínez-Muñoz, G. A comparative analysis of gradient boosting algorithms. Artif. Intell. Rev. 2020, 54, 1937–1967. [Google Scholar] [CrossRef]
  67. Kapoor, S.; Perrone, V. A Simple and Fast Baseline for Tuning Large XGBoost Models. arXiv 2021, arXiv:2111.06924. Available online: http://arxiv.org/abs/2111.06924 (accessed on 21 June 2024).
  68. Tarwidi, D.; Pudjaprasetya, S.R.; Adytia, D.; Apri, M. An optimized XGBoost-based machine learning method for predicting wave run-up on a sloping beach. MethodsX 2023, 10, 102119. [Google Scholar] [CrossRef] [PubMed]
  69. Mienye, I.D.; Sun, Y. A Survey of Ensemble Learning: Concepts, Algorithms, Applications, and Prospects. IEEE Access 2022, 10, 99129–99149. [Google Scholar] [CrossRef]
  70. Putra, M.; Rosid, M.S.; Handoko, D. Rainfall Estimation Using Machine Learning Approaches with Raingauge, Radar, and Satellite Data. In Proceedings of the International Conference on Electrical Engineering and Informatics (ICELTICs), Banda Aceh, Indonesia, 27–28 September 2022. [Google Scholar]
  71. Han, Y.; Wu, J.; Zhai, B.; Pan, Y.; Huang, G.; Wu, L.; Zeng, W. Coupling a Bat Algorithm with XGBoost to Estimate Reference Evapotranspiration in the Arid and Semiarid Regions of China. Adv. Meteorol. 2019, 2019, 9575782. [Google Scholar] [CrossRef]
  72. Mai, X.; Zhong, H.; Li, L. Research on Rain or Shine Weather Forecast in Precipitation Nowcasting Based on XGBoost. In Lecture Notes on Data Engineering and Communications Technologies; Springer Science and Business Media Deutschland GmbH: Berlin/Heidelberg, Germany, 2021; Volume 88, pp. 1313–1319. [Google Scholar] [CrossRef]
  73. Mai, X.; Zhong, H.; Li, L. Combination of XGBoost and PPLK method for improving the precipitation nowcasting. MATEC Web Conf. 2022, 355, 03039. [Google Scholar] [CrossRef]
  74. Bartz, E.; Bartz-Beielstein, T.; Zaefferer, M.; Mersmann, O. Hyperparameter Tuning for Machine and Deep Learning with R: A Practical Guide; Springer: Berlin/Heidelberg, Germany, 2022. [Google Scholar]
  75. Puspasari, R.L.; Yoon, D.; Kim, H.; Kim, K.-W. Machine Learning for Flood Prediction in Indonesia: Providing Online Access for Disaster Management Control. Econ. Environ. Geol. 2023, 56, 65–73. [Google Scholar] [CrossRef]
  76. Ramdhan, M.; Arifin, T. Application of Geographic Information System for Assessment of Indonesia Marine Proportion. J. Ilm. Geomatika 2013, 19, 141–146. [Google Scholar]
  77. Ardhitama, A.; Sholihah, R. Kajian Penentuan Awal Musim di Daerah Non ZOM 14 Riau dengan Menggunakan Data Curah Hujan dan Hari Hujan. J. Sains Teknol. Modif. Cuaca 2014, 15, 65–73. [Google Scholar] [CrossRef]
  78. Aldrian, E. Meteorologi Laut Indonesia. Badan Meteorologi dan Geofisika 2008. [Google Scholar]
  79. Dembélé, M.; Zwart, S.J. Evaluation and comparison of satellite-based rainfall products in Burkina Faso, West Africa. Int. J. Remote Sens. 2016, 37, 3995–4014. [Google Scholar] [CrossRef]
  80. Nepal, B.; Shrestha, D.; Sharma, S.; Shrestha, M.S.; Aryal, D.; Shrestha, N. Assessment of GPM-Era Satellite Products’ (IMERG and GSMaP) Ability to Detect Precipitation Extremes over Mountainous Country Nepal. Atmosphere 2021, 12, 254. [Google Scholar] [CrossRef]
  81. Ji, H.; Peng, D.; Gu, Y.; Liang, Y.; Luo, X. Evaluation of multiple satellite precipitation products and their potential utilities in the Yarlung Zangbo River Basin. Sci. Rep. 2022, 12, 13334. [Google Scholar] [CrossRef] [PubMed]
  82. Foelsche, U.; Kirchengast, G.; Fuchsberger, J.; Tan, J.; Petersen, W.A. Evaluation of GPM IMERG Early, Late, and Final rainfall estimates with WegenerNet gauge data in southeast Austria. Hydrol. Earth Syst. Sci. 2017, 21, 6559–6572. [Google Scholar] [CrossRef]
  83. Nan, L.; Yang, M.; Wang, H.; Xiang, Z.; Hao, S. Comprehensive evaluation of global precipitation measurement mission (GPM) IMERG precipitation products over mainland China. Water 2021, 13, 3381. [Google Scholar] [CrossRef]
  84. Gilewski, P.; Nawalany, M. Inter-comparison of Rain-Gauge, Radar, and Satellite (IMERG GPM) precipitation estimates performance for rainfall-runoff modeling in a mountainous catchment in Poland. Water 2018, 10, 1665. [Google Scholar] [CrossRef]
  85. Binetti, M.S.; Campanale, C.; Massarelli, C.; Uricchio, V.F. The Use of Weather Radar Data: Possibilities, Challenges and Advanced Applications. Earth 2022, 3, 157–171. [Google Scholar] [CrossRef]
  86. Sokol, Z.; Szturc, J.; Orellana-Alvear, J.; Popová, J.; Jurczyk, A.; Célleri, R. The Role of Weather Radar in Rainfall Estimation and Its Application in Meteorological and Hydrological Modelling—A Review. Remote Sens. 2021, 13, 351. [Google Scholar] [CrossRef]
  87. Bessho, K.; Date, K.; Hayashi, M.; Ikeda, A.; Imai, T.; Inoue, H.; Kumagai, Y.; Miyakawa, T.; Murata, H.; Ohno, T.; et al. An introduction to Himawari-8/9—Japan’s new-generation geostationary meteorological satellites. J. Meteorol. Soc. Jpn. Ser. II 2016, 94, 151–183. [Google Scholar] [CrossRef]
  88. Purbantoro, B.; Aminuddin, J.; Manago, N.; Toyoshima, K.; Lagrosas, N.; Sumantyo, J.T.S.; Kuze, H. Comparison of Cloud Type Classification with Split Window Algorithm Based on Different Infrared Band Combinations of Himawari-8 Satellite. Adv. Remote Sens. 2018, 7, 218–234. [Google Scholar] [CrossRef]
  89. Kim, B.-Y.; Lee, K.-T. Using the Himawari-8 AHI multi-channel to improve the calculation accuracy of outgoing longwave radiation at the top of the atmosphere. Remote Sens. 2019, 11, 589. [Google Scholar] [CrossRef]
  90. Yang, X.; You, Z.; Hiller, J.; Watkins, D. Updating and augmenting weather data for pavement mechanistic-empirical design using ASOS/AWOS database in Michigan. Int. J. Pavement Eng. 2018, 19, 1025–1033. [Google Scholar] [CrossRef]
  91. Sulistya, W.; Nugraha, H.A.; Dharmawan, G.B.; Putra, M.; Furqon, A.; Sugiarto, S.; Pramagusta, A.P. Development of Automated Weather Observing System Based on Realtime Web Display. In Proceedings of the 2019 International Electronics Symposium (IES), Surabaya, Indonesia, 27–28 September 2019; pp. 577–581. [Google Scholar]
  92. Otuaga, P.M. Rainfall Intensity-Duration Frequency Curves for Maxima Duration Series for Ondo Town, Ondo State, Nigeria. Am. Sci. Res. J. Eng. Technol. Sci. (ASRJETS) 2015, 12, 66–74. [Google Scholar]
  93. Yang, W.; Liu, J.; Li, M.; Guo, Z. Isophote-Constrained Autoregressive Model with Adaptive Window Extension for Image Interpolation. IEEE Trans. Circuits Syst. Video Technol. 2018, 28, 1071–1086. [Google Scholar] [CrossRef]
  94. Han, D. Comparison of Commonly Used Image Interpolation Methods. In Proceedings of the 2nd International Conference on Computer Science and Electronics Engineering (ICCSEE 2013), Hangzhou, China, 22–23 March 2013; pp. 1556–1559. [Google Scholar]
  95. Weerasekera, R.; Sridharan, M.; Ranjitkar, P. Implications of Spatiotemporal Data Aggregation on Short-Term Traffic Prediction Using Machine Learning Algorithms. J. Adv. Transp. 2020, 2020, 7057519. [Google Scholar] [CrossRef]
  96. Shenk, J.; Byttner, W.; Nambusubramaniyan, S.; Zoeller, A. Traja: A Python toolbox for animal trajectory analysis. J. Open Source Softw. 2021, 6, 3202. [Google Scholar] [CrossRef]
  97. Kang, D.; Rubel, O.; Byna, S.; Blanas, S. Predicting and Comparing the Performance of Array Management Libraries. In Proceedings of the 2020 IEEE 34th International Parallel and Distributed Processing Symposium (IPDPS), New Orleans, LA, USA, 18–22 May 2020; pp. 906–915. [Google Scholar] [CrossRef]
  98. Khan, S.; Maggioni, V.; Porcacchia, L. Uncertainties Associated with the IMERG Multi-Satellite Precipitation Product. In Proceedings of the IGARSS 2016—2016 IEEE International Geoscience and Remote Sensing Symposium, Beijing, China, 10–15 July 2016; pp. 2127–2130. [Google Scholar]
  99. Hartke, S.H.; Wright, D.B.; Quintero, F.; Falck, A.S. Incorporating IMERG satellite precipitation uncertainty into seasonal and peak streamflow predictions using the Hillslope Link hydrological model. J. Hydrol. X 2023, 18, 100148. [Google Scholar] [CrossRef]
  100. Omranian, E.; Sharif, H.O. Evaluation of the Global Precipitation Measurement (GPM) Satellite Rainfall Products over the Lower Colorado River Basin, Texas. J. Am. Water Resour. Assoc. 2018, 54, 882–898. [Google Scholar] [CrossRef]
  101. Liu, S.; Hu, K.; Zhang, S.; Zeng, Y. Comprehensive Evaluation of Satellite-Based Precipitation at Sub-Daily Time Scales Over a High-Profile Watershed with Complex Terrain. Earth Space Sci. 2019, 6, 2347–2361. [Google Scholar] [CrossRef]
  102. Liu, Z.; Hou, H.; Zhang, L.; Hu, B. Event-Based Bias Correction of the GPM IMERG V06 Product by Random Forest Method over Mainland China. Remote Sens. 2022, 14, 3859. [Google Scholar] [CrossRef]
  103. Partarini, N.M.C.; Sujono, J.; Pratiwi, E.P.A. Koreksi dan Validasi Data Curah Hujan Satelite GPM-IMERG dan CHIRPS di DAS Selorejo, Kabupaten Malang. In Proceedings of the Civil Engineering, Environmental, Disaster & Risk Management Symposium (CEEDRiMS), Online, 30 June–30 July 2021; pp. 149–156. [Google Scholar]
  104. Misnawati, M.; Boer, R.; June, T.; Faqih, A. Perbandingan Metodologi Koreksi Bias Data Curah Hujan CHIRPS. Limnotek 2018, 25, 18–29. [Google Scholar]
  105. Mamenun, M.; Pawitan, H.; Sopaheluwakan, A. Validasi dan koreksi data satelit trmm pada tiga pola hujan di Indonesia. J. Meteorol. dan Geofis. 2014, 15. [Google Scholar] [CrossRef]
  106. Tapiador, F.J.; Navarro, A.; García-Ortega, E.; Merino, A.; Sánchez, J.L.; Marcos, C.; Kummerow, C. The contribution of rain gauges in the calibration of the IMERG product: Results from the first validation over Spain. J. Hydrometeorol. 2020, 21, 161–182. [Google Scholar] [CrossRef]
  107. Barbieri, S.; Di Fabio, S.; Lidori, R.; Rossi, F.L.; Marzano, F.S.; Picciotti, E. Mosaicking Weather Radar Retrievals from an Operational Heterogeneous Network at C and X Band for Precipitation Monitoring in Italian Central Apennines. Remote Sens. 2022, 14, 248. [Google Scholar] [CrossRef]
  108. Piraei, R.; Afzali, S.H.; Niazkar, M. Assessment of XGBoost to Estimate Total Sediment Loads in Rivers. Water Resour. Manag. 2023, 37, 5289–5306. [Google Scholar] [CrossRef]
  109. Kavzoglu, T.; Teke, A. Advanced hyperparameter optimization for improved spatial prediction of shallow landslides using extreme gradient boosting (XGBoost). Bull. Eng. Geol. Environ. 2022, 81, 201. [Google Scholar] [CrossRef]
  110. Yin, J.; Li, N. Ensemble learning models with a Bayesian optimization algorithm for mineral prospectivity mapping. Ore Geol. Rev. 2022, 145, 104916. [Google Scholar] [CrossRef]
  111. Al-Taai, S.R.; Azize, N.M.; Thoeny, Z.A.; Imran, H.; Bernardo, L.F.A.; Al-Khafaji, Z. XGBoost Prediction Model Optimized with Bayesian for the Compressive Strength of Eco-Friendly Concrete Containing Ground Granulated Blast Furnace Slag and Recycled Coarse Aggregate. Appl. Sci. 2023, 13, 8889. [Google Scholar] [CrossRef]
  112. Wu, B.; Wei, M.; Chen, P. Bayesian optimization-based XGBoost for performance Prediction of Carbon Nanotube Membranes. Preprint 2024. [Google Scholar] [CrossRef]
  113. Zhang, B.; Zhang, B.; Gao, Y.; Gao, Y.; Wang, L.; Wang, L.; Li, T.; Li, T. BO-XGBoost-based voltage/var optimization for distribution network considering the LCOE of PV system. IET Renew. Power Gener. 2023, 18, 502–514. [Google Scholar] [CrossRef]
  114. Parthasarathy, S.; Lakshminarayanan, A.R.; Khan, A.A.A.; Sathick, K.J.; Jayaraman, V. Detection of Health Insurance Fraud using Bayesian Optimized XGBoost. Int. J. Saf. Secur. Eng. 2023, 13, 853–861. [Google Scholar] [CrossRef]
  115. Pan, X.; Wu, H.; Chen, S.; Nanding, N.; Huang, Z.; Chen, W.; Li, C.; Li, X. Evaluation and Applicability Analysis of GPM Satellite Precipitation over Mainland China. Remote Sens. 2023, 15, 2866. [Google Scholar] [CrossRef]
  116. Misnawati, F. Koreksi bias estimasi curah hujan pada satelit gpm-imerg di pulau Jawa. Megasains 2022, 13, 26–33. [Google Scholar] [CrossRef]
  117. Iqbal, Z.; Ahmed, K.; Wang, X.; Gabriel, H.F. Bias Correction Method of High-Resolution Satellite-Based Precipitation Product for Peninsular Malaysia. Preprint 2021. [Google Scholar] [CrossRef]
  118. Sorooshian, S.; AghaKouchak, A.; Arkin, P.; Eylander, J.; Foufoula-Georgiou, E.; Harmon, R.; Hendrickx, J.M.H.; Imam, B.; Kuligowski, R.; Skahill, B.; et al. Advanced concepts on remote sensing of precipitation at multiple scales. Bull. Am. Meteorol. Soc. 2011, 92, 1353–1357. [Google Scholar] [CrossRef]
  119. Serrat-Capdevila, A.; Merino, M.; Valdes, J.B.; Durcik, M. Evaluation of the performance of three satellite precipitation products over Africa. Remote Sens. 2016, 8, 836. [Google Scholar] [CrossRef]
  120. Sun, Q.; Miao, C.; Duan, Q.; Ashouri, H.; Sorooshian, S.; Hsu, K.-L. A Review of Global Precipitation Data Sets: Data Sources, Estimation, and Intercomparisons. Rev. Geophys. 2018, 56, 79–107. [Google Scholar] [CrossRef]
  121. Lasser, M.; O, S.; Foelsche, U. Evaluation of GPM-DPR precipitation estimates with WegenerNet gauge data. Atmos. Meas. Tech. Discuss. 2018, 1–30. [Google Scholar] [CrossRef]
  122. Zhang, W.; Xu, H.; Zheng, F. Aerosol optical depth retrieval over East Asia using Himawari-8/AHI data. Remote Sens. 2018, 10, 137. [Google Scholar] [CrossRef]
  123. Miles, N.L.; Verlinde, J.; Clothiaux, E.E. Cloud Droplet Size Distributions in Low-Level Stratiform Clouds. J. Atmos. Sci. 2000, 57, 295–311. [Google Scholar] [CrossRef]
  124. Morbidelli, R. Rainfall: Modeling, Measurement and Applications; Elsevier: Amsterdam, The Netherlands, 2022. [Google Scholar]
  125. Kojima, M.; Miura, T.; Furukawa, K.; Hyakusoku, Y.; Ishikiri, T.; Kai, H.; Iguchi, T.; Hanado, H.; Nakagawa, K. Dual-frequency precipitation radar (DPR) development on the global precipitation measurement (GPM) core observatory. In Earth Observing Missions and Sensors: Development, Implementation, and Characterization II; SPIE: Bellingham, WA, USA, 2012; p. 85281A. [Google Scholar] [CrossRef]
  126. Gaona, M.F.R.; Overeem, A.; Leijnse, H.; Uijlenhoet, R. First-year evaluation of GPM rainfall over the Netherlands: IMERG day 1 final run (V03D). J. Hydrometeorol. 2016, 17, 2799–2814. [Google Scholar] [CrossRef]
  127. Pejcic, V.; Garfias, P.S.; Mühlbauer, K.; Trömel, S.; Simmer, C. Comparison between precipitation estimates of ground-based weather radar composites and GPM’s DPR rainfall product over Germany. Meteorol. Z. 2020, 29, 451–466. [Google Scholar] [CrossRef]
  128. Global Precipitation Measurement (GPM) Mission Algorithm Theoretical Basis Document GPROF2017 Version 1 and Version 2 (Used in GPM V5 Processing) Passive Microwave Algorithm Team Facility. 2018. Available online: https://gpm.nasa.gov/sites/default/files/2020-05/ATBD_GPM_V5B_April15_2018.pdf (accessed on 31 July 2024).
  129. Nanding, N.; Rico-Ramirez, M.A. Precipitation Measurement with Weather Radars. In Handbook of Environmental Chemistry; Springer Science and Business Media Deutschland GmbH: Berlin/Heidelberg, Germany, 2021; Volume 102, pp. 235–258. [Google Scholar] [CrossRef]
  130. Kumar, A.; Singh, D. A Review on ‘Weather Surveillance Radar’. Int. J. Adv. Eng. 2015, 1, 239310. [Google Scholar]
  131. Nikahd, A.; Hashim, M.; Nazemosadat, M.J. A Review of Uncertainty Sources on Weather Ground-Based Radar for Rainfall Estimation. Appl. Mech. Mater. 2016, 818, 254–271. [Google Scholar] [CrossRef]
  132. Rauber, R.M.; Olthoff, L.S.; Ramamurthy, M.K.; Kunkel, K.E. The Relative Importance of Warm Rain and Melting Processes in Freezing Precipitation Events. J. Appl. Meteorol. 2000, 39, 1185–1195. [Google Scholar] [CrossRef]
  133. Curtius, J. Nucleation of atmospheric particles. EPJ Web Conf. 2009, 1, 199–209. [Google Scholar] [CrossRef]
  134. Testik, F.Y.; Barros, A.P.; Bliven, L.F. Toward a physical characterization of raindrop collision outcome regimes. J. Atmos. Sci. 2011, 68, 1097–1113. [Google Scholar] [CrossRef]
  135. Zhu, J.; Guo, B.; Luo, Y.; Pan, J. A feasible method for training classified data with sparsity. J. Phys. Conf. Ser. 2021, 1978, 012067. [Google Scholar] [CrossRef]
  136. Ali, Z.A.; Abduljabbar, Z.H.; Tahir, H.A.; Sallow, A.B.; Almufti, S.M. eXtreme Gradient Boosting Algorithm with Machine Learning: A Review. Acad. J. Nawroz Univ. 2023, 12, 320–334. [Google Scholar] [CrossRef]
Figure 1. Architecture of the ensemble learning-based rainfall estimation model using multiple instruments: rain gauges, integrated weather radars, and weather satellites.
Figure 1. Architecture of the ensemble learning-based rainfall estimation model using multiple instruments: rain gauges, integrated weather radars, and weather satellites.
Sensors 24 05030 g001
Figure 2. A map of the study area in Indonesia with three different rain patterns and a map of the weather radar network used in this research.
Figure 2. A map of the study area in Indonesia with three different rain patterns and a map of the weather radar network used in this research.
Sensors 24 05030 g002
Figure 3. Scatter plots and correlation of GPM and rain gauge data in (a) Lampung, (b) Banjarmasin, (c) Deli Serdang, (d) Pontianak, (e) Gorontalo, and (f) Biak.
Figure 3. Scatter plots and correlation of GPM and rain gauge data in (a) Lampung, (b) Banjarmasin, (c) Deli Serdang, (d) Pontianak, (e) Gorontalo, and (f) Biak.
Sensors 24 05030 g003
Figure 4. Scatter plots and correlations between GPM and Himawari satellite data in (a) Lampung, (b) Banjarmasin, (c) Deli Serdang, (d) Pontianak, (e) Gorontalo, and (f) Biak.
Figure 4. Scatter plots and correlations between GPM and Himawari satellite data in (a) Lampung, (b) Banjarmasin, (c) Deli Serdang, (d) Pontianak, (e) Gorontalo, and (f) Biak.
Sensors 24 05030 g004
Figure 5. Scatter plots and correlations between GPM and weather radar data in (a) Lampung, (b) Banjarmasin, (c) Deli Serdang, (d) Pontianak, (e) Gorontalo, and (f) Biak.
Figure 5. Scatter plots and correlations between GPM and weather radar data in (a) Lampung, (b) Banjarmasin, (c) Deli Serdang, (d) Pontianak, (e) Gorontalo, and (f) Biak.
Sensors 24 05030 g005
Figure 6. Comparison of (a) PDF and (b) CDF graphs between measured rainfall, uncorrected GPM, and corrected GPM. The x-axis represents the possible amount of rainfall. The y-axis represents the probability value.
Figure 6. Comparison of (a) PDF and (b) CDF graphs between measured rainfall, uncorrected GPM, and corrected GPM. The x-axis represents the possible amount of rainfall. The y-axis represents the probability value.
Sensors 24 05030 g006
Figure 7. Data fusion technique for multiple weather radars with frequency, polarization, and coverage range differences. There are several coverage overlaps of the weather radars with the data from each radar, which are combined into a single data output. The dot is the weather radar location, and the circle is the weather radar range.
Figure 7. Data fusion technique for multiple weather radars with frequency, polarization, and coverage range differences. There are several coverage overlaps of the weather radars with the data from each radar, which are combined into a single data output. The dot is the weather radar location, and the circle is the weather radar range.
Sensors 24 05030 g007
Figure 8. XGBoost hyperparameter optimization graph using Bayesian optimization.
Figure 8. XGBoost hyperparameter optimization graph using Bayesian optimization.
Sensors 24 05030 g008
Figure 9. Rainfall estimation products in mm/h using ensemble learning techniques and multisensor data integration produced in this study.
Figure 9. Rainfall estimation products in mm/h using ensemble learning techniques and multisensor data integration produced in this study.
Sensors 24 05030 g009
Figure 10. Image capture of GPM, radar, satellite, and estimation model results as well as comparison graphs of rain between rain gauge, GPM, and rainfall estimations during ongoing rain that occurred during very heavy rain (>20 mm/h) in Bandar Lampung on 16 December 2022 between 13.00 and 14.30.
Figure 10. Image capture of GPM, radar, satellite, and estimation model results as well as comparison graphs of rain between rain gauge, GPM, and rainfall estimations during ongoing rain that occurred during very heavy rain (>20 mm/h) in Bandar Lampung on 16 December 2022 between 13.00 and 14.30.
Sensors 24 05030 g010
Figure 11. Image capture of GPM, radar, satellite, and estimation model results as well as comparison graphs of rain between rain gauge, GPM, and rainfall estimation during ongoing rain which occurred during heavy rain (10–20 mm/h) in Biak on 26 December 2022 between 16.30 and 18.00.
Figure 11. Image capture of GPM, radar, satellite, and estimation model results as well as comparison graphs of rain between rain gauge, GPM, and rainfall estimation during ongoing rain which occurred during heavy rain (10–20 mm/h) in Biak on 26 December 2022 between 16.30 and 18.00.
Sensors 24 05030 g011
Figure 12. Image capture of GPM, radar, satellite, and estimation model results as well as comparison graphs of rain between rain gauge, GPM, and rainfall estimation during ongoing rain that occurred during moderate rain (10–20 mm/h) in Gorontalo on 17 December 2022 between 11.30 and 12.30.
Figure 12. Image capture of GPM, radar, satellite, and estimation model results as well as comparison graphs of rain between rain gauge, GPM, and rainfall estimation during ongoing rain that occurred during moderate rain (10–20 mm/h) in Gorontalo on 17 December 2022 between 11.30 and 12.30.
Sensors 24 05030 g012
Figure 13. Image capture of GPM, radar, satellite, and estimation model results as well as comparison graphs of rain between rain gauge, GPM, and rainfall estimation during ongoing rain which occurred during light rain (5–10 mm/h) in Banjarmasin on 5 December 2022 between 05.30 and 06.30.
Figure 13. Image capture of GPM, radar, satellite, and estimation model results as well as comparison graphs of rain between rain gauge, GPM, and rainfall estimation during ongoing rain which occurred during light rain (5–10 mm/h) in Banjarmasin on 5 December 2022 between 05.30 and 06.30.
Sensors 24 05030 g013
Figure 14. Image capture of GPM, radar, satellite, and estimation model results as well as comparison graphs of rain between rain gauge, GPM, and rainfall estimation during ongoing rain which occurred during very light rain (<1 mm/h) in Pontianak on 23 December 2022 between 08.00 and 09.00.
Figure 14. Image capture of GPM, radar, satellite, and estimation model results as well as comparison graphs of rain between rain gauge, GPM, and rainfall estimation during ongoing rain which occurred during very light rain (<1 mm/h) in Pontianak on 23 December 2022 between 08.00 and 09.00.
Sensors 24 05030 g014
Table 1. The research data sources from rain gauge, weather radar, and weather satellite equipment.
Table 1. The research data sources from rain gauge, weather radar, and weather satellite equipment.
InstrumentSize of DatasetProductTime ResolutionSpatial ResolutionUnit of Measurement
GPM satellite
(NASA, USA)
9.3 GiBRainfall30 min10 × 10 kmmm/h
Himawari satellite
(JMA, Japan)
1 TiBBrightness temperature10 min2 × 2 kmK
Weather radars
(EEC, USA)
50.7 TiBReflectivity10 min0.5 × 0.5 kmdBZ
Rain gauges
(All Weather, Inc, USA)
13.3 MiBRainfall1 minPointmm
Table 2. The rain gauge locations in AWOS equipment used for the bias correction of GPM rainfall products.
Table 2. The rain gauge locations in AWOS equipment used for the bias correction of GPM rainfall products.
LocationLatitudeLongitudeElevation (masl)
Bandar Lampung105.174−5.23983
Banjarmasin114.767−3.43920
Pontianak109.402−0.1422
Deli Serdang98.8843.6457
Gorontalo122.8520.63832
Biak136.104−1.1912
Table 3. The weather radars operating in Indonesia integrated into this study. These consisted of 35 weather radars located in various locations with diverse frequency bands, polarizations, and peak power.
Table 3. The weather radars operating in Indonesia integrated into this study. These consisted of 35 weather radars located in various locations with diverse frequency bands, polarizations, and peak power.
No.LocationLatitudeLongitudeElevation (masl)Frequency BandPolarizationPeak Power
1Banda Aceh5.5395.49446CSingle250 kW
2Nias1.1697.06CSingle350 kW
3Medan3.5398.6361CSingle250 kW
4Padang0.78100.324CSingle250 kW
5Pekanbaru0.45101.4631CSingle250 kW
6Bengkulu−3.85102.3415CSingle400 kW
7Jambi−1.63103.6444CSingle400 kW
8Palembang−2.91104.712CSingle250 kW
9Pangkalpinang−2.16106.1430CSingle350 kW
10Lampung−5.2105.17106CSingle250 kW
11Cengkareng−6.17106.6425CSingle250 kW
12Pontianak−0.08109.3926CSingle250 kW
13Sintang−0.04111.4528CDual400 kW
14Pangkalanbun−2.73111.6431CSingle400 kW
15Banjarmasin−3.46114.8481CSingle250 kW
16Balikpapan−1.25116.8950CSingle250 kW
17Tarakan3.31117.5845CSingle250 kW
18Yogyakarta−7.73110.35182CSingle350 kW
19Surabaya−7.41112.763CSingle250 kW
20Denpasar−8.73115.1728CSingle250 kW
21Lombok−8.75116.2494CSingle400 kW
22Bima−8.54118.6845CSingle250 kW
23Maumere−8.61122.0836CSingle400 kW
24Kupang−10.21123.62326CDual400 kW
25Majene−3.55118.9830XSingle2 × 500 kW
26Makassar−4.99119.5711CSingle250 kW
27Masamba−2.55120.3266XSingle2 × 500 kW
28Gorontalo0.63123.0190CSingle250 kW
29Ternate0.85127.34105CSingle400 kW
30Manado1.5129.9116CSingle250 kW
31Ambon−3.71128.099CSingle250 kW
32Biak−1.16136.0872CSingle250 kW
33Sorong−0.89131.2822CSingle250 kW
34Timika−4.52136.8954CSingle250 kW
35Merauke−8.49131.2888CSingle250 kW
Table 4. Hyperparameter tuning for XGBoost.
Table 4. Hyperparameter tuning for XGBoost.
HyperparameterRangeDefinition
learning_rate 0.01–1The step size when updating model weights to minimize errors, which affects the speed and convergence of the training process.
max_depth0–12The maximum depth of a tree, which controls the complexity of the model by limiting the number of levels of splitting in each tree.
n_estimators100–1000The total number of decision trees to be built and used in an ensemble model, which directly affects the performance and complexity of the model.
subsample0.1–1The proportion of training data samples used to build each tree, introducing variation in the training process.
min_child_weight0.1–2The sum of instance weights required at a leaf node, which ensures that nodes will not split if they do not meet this minimum weight threshold.
gamma0–1The minimal reduction in loss required to split nodes, which helps control tree growth by preventing insignificant splits.
colsample_bytree0.1–1The proportion of features (columns) randomly selected to build each tree, which helps prevent overfitting by reducing the correlation between trees.
Table 5. Hyperparameter and optimal values obtained by each Bayesian optimization strategy.
Table 5. Hyperparameter and optimal values obtained by each Bayesian optimization strategy.
HyperparameterOptimal Value
learning_rate 0.04
max_depth1
n_estimators886
subsample0.96
min_child_weight0.14
gamma0.08
colsample_bytree0.45
Table 6. The results of the rainfall estimation model evaluation.
Table 6. The results of the rainfall estimation model evaluation.
LocationRMSEAccuracy
Bandar Lampung2.750.89
Banjarmasin2.570.91
Pontianak3.080.89
Deli Serdang2.640.9
Gorontalo1.850.92
Biak2.480.9
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

Putra, M.; Rosid, M.S.; Handoko, D. High-Resolution Rainfall Estimation Using Ensemble Learning Techniques and Multisensor Data Integration. Sensors 2024, 24, 5030. https://doi.org/10.3390/s24155030

AMA Style

Putra M, Rosid MS, Handoko D. High-Resolution Rainfall Estimation Using Ensemble Learning Techniques and Multisensor Data Integration. Sensors. 2024; 24(15):5030. https://doi.org/10.3390/s24155030

Chicago/Turabian Style

Putra, Maulana, Mohammad Syamsu Rosid, and Djati Handoko. 2024. "High-Resolution Rainfall Estimation Using Ensemble Learning Techniques and Multisensor Data Integration" Sensors 24, no. 15: 5030. https://doi.org/10.3390/s24155030

APA Style

Putra, M., Rosid, M. S., & Handoko, D. (2024). High-Resolution Rainfall Estimation Using Ensemble Learning Techniques and Multisensor Data Integration. Sensors, 24(15), 5030. https://doi.org/10.3390/s24155030

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