Next Article in Journal
A Technical Review of Planet Smallsat Data: Practical Considerations for Processing and Using PlanetScope Imagery
Previous Article in Journal
Modelling the Spectral Index to Detect a Baltic-Type Crude Oil Emulsion Dispersed in the Southern Baltic Sea
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evaluation of Satellite Precipitation Estimates over the South West Pacific Region

1
School of Earth, Atmosphere, and Environment, Monash University, Melbourne 3800, Australia
2
Climate Risk and Early Warning Systems (CREWS), Bureau of Meteorology, Melbourne 3008, Australia
3
SPACE Research Centre, School of Science, Royal Melbourne Institute of Technology (RMIT) University, Melbourne 3000, Australia
*
Author to whom correspondence should be addressed.
Remote Sens. 2021, 13(19), 3929; https://doi.org/10.3390/rs13193929
Submission received: 15 July 2021 / Revised: 24 September 2021 / Accepted: 27 September 2021 / Published: 30 September 2021

Abstract

:
Rainfall estimation over the Pacific region is difficult due to the large distances between rain gauges and the high convection nature of many rainfall events. This study evaluates space-based rainfall observations over the South West Pacific Region from the Japan Aerospace Exploration Agency’s (JAXA) Global Satellite Mapping of Precipitation (GSMaP), the USA National Oceanographic and Atmospheric Administration’s (NOAA) Climate Prediction Center morphing technique (CMORPH), the Climate Hazards group Infrared Precipitation with Stations (CHIRPS), and the National Aeronautics and Space Administration’s (NASA) Integrated Multi-Satellite Retrievals for GPM (IMERG). The technique of collocation analysis (CA) is used to compare the performance of monthly satellite precipitation estimates (SPEs). Multi-Source Weighted-Ensemble Precipitation (MSWEP) was used as a reference dataset to compare with each SPE. European Centre for Medium-range Weather Forecasts’ (ECMWF) ERA5 reanalysis was also combined with Soil Moisture-2-Rain–ASCAT (SM2RAIN–ASCAT) to perform triple CA for the six sub-regions of Fiji, New Caledonia, Papua New Guinea (PNG), the Solomon Islands, Timor, and Vanuatu. It was found that GSMaP performed best over low rain gauge density areas, including mountainous areas of PNG (the cross-correlation, CC = 0.64), and the Solomon Islands (CC = 0.74). CHIRPS had the most consistent performance (high correlations and low errors) across all six sub-regions in the study area. Based on the results, recommendations are made for the use of SPEs over the South West Pacific Region.

Graphical Abstract

1. Introduction

Rainfall is fundamental to sustaining communities, economy, and the natural environment. Numerous climate-sensitive sectors rely upon accurate precipitation measurements, and the demand for greater accuracy will increase as rainfall variability increases globally due to anthropogenic climate change [1]. Accurate rainfall estimation crucially impacts a range of sectors, e.g., agriculture, fire and forestry management, landslide prevention [2], and hydrological modelling of reservoirs [3].
Rain gauges and weather radars have been the traditional means of measuring precipitation, yet satellite remote sensing is increasingly considered as complementary to rain gauges due to its improving accuracy and spatial coverage [4,5]. Rain gauges themselves are prone to errors from evaporative losses, maintenance, and distribution issues (especially in remote areas), and the small representative area is particularly problematic in regions with high topography and convection rainfall [6]. Weather radars can be effective for capturing the spatial extent of precipitation, but their use is limited, particularly in developing and least developed countries, due to large upfront expense and maintenance costs [7].
Satellite precipitation estimates (SPEs) provide data with potentially global coverage that are maintained in the long-term by reputable national space agencies and meteorological services, available in near-real-time, and for the most part available free-of-charge online [6]. On the other hand, SPEs generally have relatively low spatial resolution and are less accurate than rain gauges or weather radars at the local scale [6].
Due to the assumptions made in SPE retrieval algorithms, it is vital to verify their accuracy [8]. Traditionally, rain gauges and weather radars have been used as ground validation to test the accuracy of SPEs [8,9]. Recently, the technique of collocation analysis (CA) has been used to evaluate the accuracy of SPEs in regions with insufficient ground-validation infrastructure [10,11]. In this study, CA was selected as a novel methodology to investigate the accuracy of various SPEs over a rain gauge sparse region—the South West Pacific.
The small island developing states and least developed countries in the Pacific are among the world’s poorest nations, with a low GDP per capita, and a high vulnerability to natural hazards [12]. These island nations are also highly exposed to meteorological hazards, both drought and acute water shortages and, conversely, heavy precipitation leading to flooding [13,14]. Climate change is altering precipitation patterns [1], causing reduced total rainfall in some areas and increasing the risk of severe drought in the region [15,16].
Despite efforts to improve the rain gauge record, the South West Pacific Islands is one of the regions where surface-based precipitation observations are very sparse [15,16]. Therefore, SPEs are an obvious choice for precipitation monitoring to help vulnerable Pacific communities, potentially enabling better management of hydrometeorological hazards such as droughts and floods.
CA was used on a global scale for the first time in 2017 for evaluating SPEs by comparing them to a satellite rainfall product derived from soil moisture data and a climate reanalysis product [10]. Since then, the methodology has been used for comparing the accuracy of SPEs globally and across various regions [17,18,19]. The recent comparison of SPEs by Beck et al. [20] revealed that ERA5 was the most accurate ungauged product, while gauge-based SPEs were more accurate than gauge-corrected ones. Tanim et al. [19] used ERA-Interim and interpolated gauge data in a CA to show that IMERG Final run was more accurate than IR-based SPEs including CMORPH. Global CAs have revealed that for tropical rainforests, such as those common to Pacific Island Nations, SM2RAIN–ASCAT and ERA5 do not perform as well as SPEs [10,21]. However, there have been no CA studies analysing SPE accuracy over the South West Pacific Region, and this research addresses the identified gap in the literature.
This paper is organised as follows. Section 2 describes the study area, datasets, and methods used in the study. Section 3 presents the results, while Section 4 discusses the findings and provides recommendations for SPE usage. Section 5 summarises the major findings.

2. Materials and Methods

2.1. Study Area

A map of the study area and relevant sub-regions is presented in Figure 1.
The study area extends from 0° N to 30° S, and from 120° E to the International Date Line (further extension eastward was limited by the availability of the Soil Moisture product. Soil Moisture 2 Rain–ASCAT (SM2RAIN–ASCAT). The study area includes the six sub-regions: Timor, Papua New Guinea (PNG), Bougainville Island and Solomon Islands, Vanuatu, New Caledonia, and the main islands of Fiji. The study period is limited to the temporal range of SM2RAIN–ASCAT, from January 2007 to December 2020.

2.2. Datasets

Four satellite precipitation datasets commonly used in evaluation studies were selected for this study: the Japan Aerospace Exploration Agency’s (JAXA) Global Satellite Mapping of Precipitation (GSMaP V6), the USA National Oceanographic and Atmospheric Administration’s (NOAA) Climate Prediction Center morphing technique (CMORPH V1), Climate Hazards group Infrared Precipitation with Stations (CHIRPS V2), and National Aeronautics and Space Administration’s (NASA) Integrated Multi-Satellite Retrievals for GPM (IMERG V6b).
The Global Precipitation Mission (GPM) is a constellation of low earth orbiting satellites which uses PMW imagers and sounders and a dual frequency radar for remotely detecting precipitation; it is maintained by both NASA and JAXA. GSMaP [22] and IMERG [23] provide several SPEs created using different precipitation retrieval algorithms based on data from the GPM constellation.
In addition to utilising MW-based estimates from the GPM constellation, GSMaP also uses infrared (IR) data from geostationary satellites to observe the movement of cloud systems and advection estimates, as well as applying a Kalman filter to form better estimates when data is sparse [24]. A daily global dataset of rain gauges from the CPC is used for calibration via the matching over 24-hour totals, which allows for more accurate calibration of MW readings than the use of monthly datasets [25].
CMORPH uses a similar cloud advection algorithm to GSMaP and the same rain gauge network for daily calibration [26]. CMORPH MW inputs have been calibrated by matching the probability distribution functions of the estimates to those from rain gauge values over land and through the use of corrective ratios formed from matching means to the GPCP dataset over the ocean. In this study, the CMORPH Blended (CMORPH-BLD) product which further incorporates gauge data through optimal interpolation was selected for evaluation. CMORPH-BLD has been shown to be more accurate than the gauge-corrected product CMORPH-CRT, over regions where a dense gauge network exists [27], though performance was similar when gauge density is low such as, for example, over PNG [28].
GSMaP and CMORPH data were made available through the World Meteorological Organization’s Space-based Weather and Climate Extremes Monitoring Demonstration Project [29].
IMERG Final Run monthly product blends the Global Precipitation Climatology Centre (GPCC) V5 Monitoring Product and is retroactively used to gauge-correct the Early and Late runs [30]. IMERG also uses a cloud advection algorithm and a Kalman filter. Morphed MW estimates are used where their availability is sufficient, with an IR-based SPE called PERSIANN-CCS being used otherwise. Although GSMaP and IMERG are similar in selecting the satellite sensors that they use, the difference in rain gauge data used for calibration and in their precipitation retrieval algorithms results in different performance.
The CHIRPS dataset was created as part of the USA Agency for International Development Famine Early Warning Systems Network (FEWS NET) [5]. It forms a climatological field using long-term trends from both gauge and satellite products (including CMORPH and Land Surface Temperature estimates) onto which precipitation estimate from near real time IR Cold Cloud Duration (CCD) are incorporated. Blending also takes place with interpolation from up to five nearby gauges. Its higher resolution (0.05° compared to 0.1° or 0.25° of other datasets) and temporal coverage (beginning from 1981 compared to around the 2000s for other SPEs) may make it applicable for Pacific Island Nations that lack consistent rain gauge history.
All SPEs require calibration to surface rainfall [31] and the choice of ground-truth data impacts the accuracy of the satellite products. CMORPH blends in CPC daily data, CHIRPS blends in nearest gauges, IMERG blends in monthly rainfall data, and GSMaP corrects to daily data. It has been shown that using daily gauge data resulted in GSMaP being more accurate over the Tibetan Plateau than IMERG [25].
The Multi-Source Weighted-Ensemble Precipitation dataset (MSWEP V2.8) is a merged product combining numerous gauge, reanalysis, and satellite products [32] and will serve as a gridded reference dataset for this study. However, as ERA-INT, GSMaP, and CMORPH are inputs to MSWEP, it thus will not be used in collocation analysis performed in this study.
Soil Moisture 2 Rain (SM2RAIN–ASCAT) [33] and ERA5 [34] use precipitation derivation and reanalysis methods, respectively, to estimate rainfall and are often used in CA to compare SPEs [10,18,20].
SM2RAIN–ASCAT considers an inverse water equation where rainfall is estimated from observed soil moisture and modelled soil temperature and evaporation [33]. It can be more effective than other MW based SPEs when an estimation of accumulated rainfall is desired as opposed to the instantaneous estimates generated by SPEs [33] which are subject to sampling bias. The SM2RAIN–ASCAT product is used in this study, as it has the largest temporal range of all SM2RAIN products covering January 2007 until December 2020. The ASCAT product is based upon three MetOp C-Band MW sounders [21].
The European Centre for Medium-Range Weather Forecasts (ECMWF) produces the ERA5 reanalysis product, which uses a wide range of weather observations from multiple sources including satellites, weather balloons, aircrafts, and weather stations, for humidity, air pressure, wind, and temperature, to provide an accurate history of global meteorological observations from 1950 to present [34]. It does not directly ingest rain gauge data, with a second-order inclusion occurring from the use of a gauge-adjusted radar analysis over the USA, in addition to the use of satellite data. A comparison of ERA5 to global rain gauge datasets reveals that ERA5 overestimates rainfall over the oceans, but is in very close agreement for precipitation over land [34].

2.3. Method

CA is a methodology that assumes an unknown uncertainty for a geophysical process (in this case precipitation) and then compares at least three datasets made up of mutually independent measurements to determine their relative accuracy [11,35]. There are two main types of CA which make different assumptions about errors. There is a simple additive error model:
R = a + B   T + ϵ
and a multiplicative error model, which is considered more accurate for rainfall data due to the sporadic nature of rainfall amounts [11]:
ln ( R ) = a   T B   ln ( ϵ )
where R is the observed rainfall amount, T is the unknown true rainfall amount, ϵ is the random error, and a and B denote systematic biases.
CA methods overcome the need for an exact source of truth, and are spatially consistent across a domain [36]. A noteworthy feature of CA is that the methodology relies upon the contributing datasets being independent of each other. More specifically, CA’s underlying assumptions are of (i) stationarity of the statistics, (ii) linearity between at least three estimates (versus the same target) across all timescales, and (iii) existence of uncorrelated error between at least three estimates [10,11].
Monthly averaged daily precipitation data (January 2007–December 2020) was downloaded as NetCDF4 files for all datasets except IMERG (downloaded as HDF5), converted to units of mm per day, and regridded to a common resolution of 0.25°, as this was the native resolution of SM2RAIN–ASCAT, ERA5, and CMORPH datasets. The land-based SPEs of CHIRPS, CMORPH, and SM2RAIN–ASCAT were extrapolated to smooth around coastlines, and a common land–sea mask from CMORPH was applied for all the datasets to remove the variation of masks amongst the datasets. CMORPH’s mask was used as it presented a compromise of coverage over Pacific Island Countries with the other land-based SPEs (CHIRPS, SM2RAIN–ASCAT) as well as the land–sea mask of ERA5.
Due to the performance of each product varying in accuracy throughout the year [17] it is conventional to use monthly anomalies rather than rainfall totals. Using anomalies restricts our analysis to assume an additive error model (Equation (1)), rather than the multiplicative variant (Equation (2)). Despite the benefits of a multiplicative error model, recent studies, such as Chen et al. [17] and Massari et al. [10], used an additive error model which still provided robust results. Therefore, it is reasonable to remove the climatological component of the data from all the datasets and perform an additive error analysis.
MSWEP was used as a reference dataset to compare each SPE. It is considered to be the most accurate near-real-time satellite product available due to its optimised merging of other rainfall products [20]. The Mean Bias, Mean Average Error, Root Mean Square Error, and Pearson-R correlation are used, as they are conventional comparative statistical measures used to understand dataset accuracy in SPE research (formulas can be found in [28]). Assumptions of linearity and normality also were validated through plotting to ensure that Pearson correlation was the correct statistic for this process. Fisher’s z-transformation comparison of correlation coefficients was used to determine significance of the Pearson correlations.
The upper and lower quantiles for each SPE were also compared with MSWEP to record how consistently a value in the lower 20% or upper 80% of the reference dataset was identified by the SPE. The metric of this comparison is percent hit rate, which compares the number of events that lie in the top and bottom quantile of the SPE data to the number expected from MSWEP. While daily data may be more relevant for the accuracy of precipitation estimates for individual tropical cyclone events, the use of monthly data can be relevant to climate extremes such as droughts and accumulated heavy precipitation during the Southern Hemisphere tropical cyclone season, typically from November through to April [37,38].
The time series of one SPE, SM2RAIN–ASCAT and ERA5 were used to create a 3 × 3 covariance matrix (C) for each pixel. This matrix was then used to calculate the cross-correlation (CC) and RMSE between the SPE and the unknown true rainfall for each SPE. As we are using an additive error model, the equations are:
C C = C X , Y   C X , Z C X , X   C Y , Z
R M S E = C X , X C X , Y   C X , Z   C Y , Z
where X represents the SPE being considered, Y is SM2RAIN–ASCAT, and Z is ERA5.
Assumption of stationarity were tested by the Augmented Dickey–Fuller unit root test, and linearity by visual inspection. It is worth noting that the assumptions of these datasets having uncorrelated errors are not completely met, as some sensors from the GPM constellation are included in ERA5 [34]. Although SM2R is best not used with ERA5 in triple CA [21], it is still used in cases where a suitable rain gauge product is not available [10,17]. Similarly, Quadruple Collocation using the four SPEs was explored, but was not utilized due to the overlap of sensors violating the assumptions of independence.

3. Results

3.1. Performance of Datasets Compared to MSWEP

Figure 2 shows the spatial differences in Mean Bias errors (MBE) between each SPE and the reference MSWEP for six sub-regions. It is evident that topography impacts bias in SPEs, with often opposite direct biases between highlands and lowlands. In general, there is an overestimation of precipitation over mountains. Around coastal areas, the SPEs are more likely to underestimate rainfall, which is consistent with issues of SPEs measuring rainfall over oceans [39].
The comparison of each of the SPEs to MSWEP is summarised in Figure 3; see Appendix A for details and comparison, including SM2RAIN–ASCAT. GSMaP and IMERG have stronger yet insignificant Pearson correlations in Solomon Islands, Timor, and PNG; however, CMORPH and CHIRPS have significantly greater performance when over more densely gauged subregions of Fiji, Vanuatu, and New Caledonia (Figure 3b). CMORPH has the largest RMSE over PNG, GSMaP has the largest over Fiji, Vanuatu, New Caledonia, and Solomon Islands, while IMERG has the smallest or near smallest errors over these sub-regions (Figure 3c). SM2RAIN–ASCAT had significantly lower average correlations than the other SPEs for all sub-regions, with the exception of insignificantly lower R values in Timor (Significance P = 0.86) and PNG (P = 0.19).
A comparison of the percent hit rate for the upper and lower quantiles with MSWEP indicated there was no consistent difference between the accuracy of lower and upper quantiles for any of the SPEs when compared to MSWEP. The gauge-based products of CMORPH-BLD and IMERG show greater quantile accuracy in Fiji and New Caledonia, both sub-regions with denser gauge networks. CHIRPS had overall the lowest average hit rate (90.6%), while GSMaP performed well for sub-regions with low gauge density.

3.2. Triple Collocation Analysis

As in the comparison to MSWEP, we found poor performance of the SPEs over PNG compared to the other sub-regions. While most of the datasets show agreement, there are some areas where performance across the datasets diverges.
Cross correlation between each SPE and the unknown truth from triple CA of ERA5 and SM2RAIN–ASCAT are shown in Figure 4. All SPEs demonstrated best performance over Fiji, New Caledonia, and Vanuatu. GSMaP shows far greater CC over PNG, particularly over highlands, than the other datasets. IMERG appears to perform weaker than other SPEs over PNG, Timor, and the Solomon Islands, but has comparable performance over Fiji, Vanuatu, and New Caledonia.
Comparisons of the mean CC and RMSE using triple CA for each SPE are shown in Figure 5a,b, respectively. GSMaP outperformed other SPEs over PNG and the Solomon Islands, while CHIRPS outperformed over Timor. CHIRPS had greater or equal correlations and lower RMSE than IMERG and CMORPH. IMERG did not perform as well as other SPEs over PNG, Timor, the Solomon Islands, and Fiji, whereas its performance was comparable over New Caledonia and Vanuatu.
A comparison of the SPEs with the highest correlation at each pixel for each sub-region is provided in Figure 6; the second highest correlation is also presented.
GSMaP performs best over mountainous areas of PNG, Timor, and the Solomon Islands, while CHIRPS shows consistent performance across the study region, except for poor performance over PNG. Performance of CHIRPS appears best around some coastlines and lower-lying areas of Timor and New Caledonia. CMORPH performs well over Fiji, and IMERG is comparable with CHIRPS over New Caledonia.

4. Discussion

Evidence that SPEs over PNG are the least accurate of the six regions reinforces the known difficulties of estimating rainfall from satellites over areas with complex topography, dense tropical rainforest, and low rain-gauge coverage [28]. Studies on evaluating SPEs in Pacific Island Countries are very limited, with the only previous study over PNG performed by Chua et al. (2020) [28] comparing CMORPH, GSMaP, and ERA5 with rain gauge data. The results obtained in this study showed similar but consistently higher correlations and smaller errors than those presented in Chua et al. [28], likely due to the use of anomaly values rather than raw rainfall totals.
As expected, CHIRPS’s weaker quantile performance was due to the blending technique from multiple nearby gauges rather than the nearest neighbour [5].
Including rain gauge data within the comparison of MSWEP and ERA5 would be a useful future study to benefit users in the Pacific Islands and for the wider field of collocation analysis.
Using triple CA, stronger CCs in less mountainous areas were found, supporting earlier studies that SPEs are more accurate at low elevations [40,41].
CMORPH performed best over Fiji, likely due to a higher density of contributing rain gauges [16]. This supports previous findings that the CMORPH-BLD is dependent on the availability of CPC rain gauge data [27,28]. It appears that IMERG’s strong performance over New Caledonia is also due to the presence of a dense GPCC rain gauge network, the difference from the Fiji result being due to fewer contributing rain gauges in the GPCC dataset than the CPC [42,43]. As this difference in the rain gauge data is most pronounced in the case of Fiji, we can assume that IMERG would have higher correlations than CMORPH over Fiji if the same rain gauge dataset was used for calibration.
Although GSMaP has been shown to have a significant positive bias, particularly over PNG [44], it still demonstrates the strongest CCs. While the SPEs have similar correlations over Vanuatu, Fiji, and New Caledonia, CHIRPS is shown to perform most similarly to ERA and SM2RAIN–ASCAT by having the smallest errors.
The main limitations of using SM2RAIN–ASCAT are that it (i) underestimates peak rainfall events due to soils becoming saturated and incurring surface runoff, (ii) incorrectly records spurious low-intensity rainfall events due to high-frequency soil moisture fluctuations associated with random measurement error, and (iii) is limited to only liquid-phase precipitation over land [21]. The accuracy of SM2RAIN–ASCAT is being improved through the use of constellations of longer-wave MW sensors [45,46] incorporating newer technologies from Sentinel-1 [46] and the Cyclone Global Navigation Satellite System (CYGNSS) constellation [47].
The CC values from triple CA are considered a more important metric for determining accuracy than RMSE [18,36], and the errors are used here to contextualise the results.
It was found that CHIRPS has relatively high CCs compared to the other SPEs for all sub-regions. In sub-regions with low rain gauge density, GSMaP is recommended for use due to its high correlations and low latency. GSMaP performs better over PNG, and CHIRPS over Timor. For Vanuatu, the Solomon Islands, and Bougainville, CHIRPS has slightly better quantile performance and lower errors, while GSMaP has stronger correlations and smaller biases.
As GSMaP uses a similar satellite constellation to IMERG, it is expected that it will be able to provide a more accurate spatiotemporal distribution of tropical cyclone rainfall than CHIRPS [48]. Given the latency for GSMaP is shorter than for CHIRPS (3 days versus 1 month), it is preferred for rapid tropical cyclone impact assessment.
Over New Caledonia, it is recommended to use IMERG due to its dense GPCC gauge network and quantile analysis performance. CHIRPS proves to be a close second, however IMERG is recommended due to its more homogenous distribution of bias and stronger CC from triple CA.
Similarly within Fiji, it can be confidently recommended that CMORPH be used for applications due to its strong CCs. The negative bias around coasts could lead to greater rates of false alarms for drought, so it is advised that if CMORPH historical datasets are used for statistical precipitation forecasting, the inclusion of in situ verification of drought conditions is recommended. In the case of probabilistic drought forecasting, these biases will not impact the seasonal forecasts themselves, but we caution the coupling of CMORPH and forecasts for drought monitoring to similarly include in situ measurements where possible. If other countries had similar amounts of rain gauge data in the CPC dataset as Fiji, then CMORPH could also become the recommended SPE. If similar gauge data was included in the GPCC dataset, then IMERG could become the recommended dataset for use over Fiji due to its strong performance over New Caledonia.
Analysis indicates that SPE accuracy varies over a range of variables and geographic features. Careful consideration of strengths and limitations of SPEs is required before implementing them in operational services. It is recommended that Pacific Island Country National Meteorological Services seeking to operationalise a certain SPE dataset perform correlation analyses on the SPE datasets available to them.

5. Conclusions

This research found that a key factor for SPE accuracy in the study region was topography, with more mountainous sub-regions (PNG, Timor, the Solomon Islands) consistently having weaker correlations than those characterised by less mountainous terrain (Fiji, Vanuatu, and New Caledonia).
A comparison to MSWEP was performed and results were in line with those obtained in a previous analysis by Chua et al. (2020) [28], which showed GSMaP to be more accurate than CMORPH over PNG. The triple CA found that GSMaP performed best in particularly mountainous areas and CHIRPS had the most consistent performance across all six sub-regions in the study region.
These different levels of analysis demonstrated that the choice of dataset can reveal quite different accuracy results. The blended datasets of IMERG and CMORPH performed well for the Fiji and New Caledonia sub-regions where high density rain gauge observations were representative of the small domains. GSMaP appears to perform best in sub-regions of low rain gauge density.
Recommendations for SPE for use during tropical cyclone seasons and drought periods were inferred, however, there remains a variety of other avenues of research that could provide more specific insights for applying SPEs to natural hazard risk reduction. CHIRPS was recommended for drought applications due to its consideration of historical climate data, whereas GSMaP was favoured for tropical cyclone seasons with strong performance over regions with low density or no rain gauges, and lower latency.
Case studies using daily raw data in a multiplicative error model could provide more relevant performance of tropical cyclone rainfall accuracy, and new techniques of soil moisture measurement and SPE merging can provide improved ability for SPEs to contribute to disaster risk reduction efforts.
This research reaffirms the importance of rain gauges for increasing SPE accuracy, and strongly recommends maintaining, and even expanding where possible, a network of surface-based meteorological observation stations. The methodology established in this study for the South West Pacific Region can be replicated for other precipitation data sparse regions, and would be a crucial step in the process of operationalising SPE applications.

Author Contributions

Conceptualization, A.W. and Y.K.; methodology, A.W.; software, A.W. and Z.-W.C.; validation, A.W.; formal analysis, A.W.; investigation, A.W.; writing—original draft preparation, A.W. and Y.K.; writing—review and editing, A.W., Z.-W.C., and Y.K.; supervision, Y.K.; project administration, Y.K.; funding acquisition, Y.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

GSMaP data were provided by JAXA; CMORPH by CPC, NOAA; CHIRPS by FEWS NET; IMERG by NASA; MSWEP by GloH20; ERA5 by ECMWF; and SM2RAIN–ASCAT by ECSD. Contains modified Copernicus Climate Change Service Information (2019). Neither the European Commission nor ECMWF is responsible for any use that may be made of the Copernicus Information or Data it contains. Comments from Andrew B. Watkins, David A. Jones, four anonymous reviewers, and MDPI academic editor helped to improve quality of the manuscript. We are grateful to Caroline Poulsen and other colleagues from the Bureau of Meteorology and Monash University for their helpful advice and guidance.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Tabulated Results of Metrics

Table A1. Results of SPE comparison to MSWEP. Red indicates highest value, blue indicated lowest.
Table A1. Results of SPE comparison to MSWEP. Red indicates highest value, blue indicated lowest.
TimorPNGFijiVanuatuSolomon IslandsNew CaledoniaAverage
MBEMeanCMORPH−0.363−1.704−0.245−1.392−3.052−0.709−1.244
GSMaP−0.552−1.151−1.166−1.417−0.271−0.487−0.841
CHIRPS0.097−0.974−1.019−1.748−1.962−0.590−1.033
SM2RAIN-ASCAT0.5161.052−1.027−0.8790.502−0.207−0.007
IMERG0.089−0.374−0.220−1.864−0.5000.257−0.435
MedianCMORPH−0.421−1.675−0.338−1.365−2.685−0.905−1.232
GSMaP−0.638−1.094−1.029−1.244−0.414−0.425−0.807
CHIRPS0.050−0.911−1.391−2.318−1.275−0.616−1.077
SM2RAIN-ASCAT0.6750.620−0.965−0.6630.353−0.354−0.056
IMERG0.126−0.468−0.159−1.670−0.5230.336−0.393
QuintLower CMORPH93.88187.02893.95494.36994.90993.98493.021
GSMaP95.78887.77290.06094.07095.35690.24192.215
CHIRPS91.17685.55690.90493.06494.16490.52590.898
SM2RAIN-ASCAT91.60781.94788.45389.65192.78690.15789.100
IMERG93.77388.13692.48494.36995.58892.94892.883
UpperCMORPH91.67284.38993.30192.24695.54293.26591.736
GSMaP93.19186.02490.87792.27795.65891.29391.554
CHIRPS90.61682.26290.82291.11494.76792.49790.346
SM2RAIN-ASCAT90.84281.37788.37189.21092.99890.72588.921
IMERG91.71586.16091.06893.23795.96993.93492.014
MAEMeanCMORPH1.0671.9281.3631.6511.6401.1411.465
GSMaP0.9891.8221.7801.8922.1851.3921.677
CHIRPS1.0481.8711.5921.8261.7731.1681.547
SM2RAIN-ASCAT0.9492.1942.1252.3782.3791.4391.911
IMERG0.8601.6111.5471.5741.5311.1061.372
MedianCMORPH1.0281.8281.2771.6051.6181.1441.416
GSMaP0.9691.7951.8121.8422.1731.4331.671
CHIRPS0.9831.7901.5541.8251.7521.1751.513
SM2RAIN-ASCAT0.9202.1202.0842.3612.3461.4911.887
IMERG0.8071.5081.5111.5981.5311.0791.339
PearsonMeanCMORPH0.7350.7030.8860.7910.7500.8590.787
GSMaP0.8020.7400.8510.7740.7670.8210.792
CHIRPS0.7320.6070.8410.7400.6910.8520.744
SM2RAIN-ASCAT0.7230.5090.6680.5310.4270.7320.598
IMERG0.7880.7200.8650.7870.7880.8650.802
MedianCMORPH0.7370.7090.8940.7870.7540.8580.790
GSMaP0.8030.7450.8500.7800.7720.8150.794
CHIRPS0.7420.6110.8310.7400.7050.8490.747
SM2RAIN-ASCAT0.7380.5100.6710.5340.4270.7730.609
IMERG0.7990.7210.8710.7730.8100.8800.809
RMSEMeanCMORPH1.6582.5892.0322.4822.1981.6422.100
GSMaP1.5252.4342.4752.8562.8761.8922.343
CHIRPS1.5692.5112.3322.7262.3761.6322.191
SM2RAIN-ASCAT1.3712.8663.1363.3673.0331.9912.627
IMERG1.2942.1342.1482.4942.0341.5641.945
MedianCMORPH1.6002.4521.9032.5092.1501.6602.046
GSMaP1.5022.4112.5032.7562.8471.9512.328
CHIRPS1.4912.4432.3112.6582.3561.6892.158
SM2RAIN-ASCAT1.3182.6923.1023.3262.9832.1072.588
IMERG1.2252.0062.1332.5332.0571.5471.917

References

  1. IPCC. Climate Change 2014: Synthesis Report. Contribution of Working Groups I, II and III to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; Core Writing Team, Ed.; IPCC: Geneva, Switzerland, 2014; p. 151. [Google Scholar]
  2. Stanley, T.; Kirschbaum, D.B.; Pascale, S.; Kapnick, S. Extreme precipitation in the Himalayan landslide hotspot. In Satellite Precipitation Measurement; Springer: Cham, Switzerland, 2020; pp. 1087–1111. [Google Scholar]
  3. Li, X.; Chen, Y.; Deng, X.; Zhang, Y.; Chen, L. Evaluation and hydrological utility of the GPM IMERG Precipitation products over the Xinfengjiang river reservoir basin, China. Remote Sens. 2021, 13, 866. [Google Scholar] [CrossRef]
  4. Dinku, T. The value of satellite rainfall estimates in agriculture and food security. In Satellite Precipitation Measurement; Springer: Cham, Switzerland, 2020; pp. 1113–1129. [Google Scholar]
  5. Funk, C.; Peterson, P.; Landsfeld, M.; Pedreros, D.; Verdin, J.; Shukla, S.; Husak, G.; Rowland, J.; Harrison, L.; Hoell, A.; et al. The climate hazards infrared precipitation with stations—A new environmental record for monitoring extremes. Sci. Data 2015, 2, 150066. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Tarnavsky, E.; Bonifacio, R. Drought risk management using satellite-based rainfall estimates. In Satellite Precipitation Measurement; Springer: Cham, Switzerland, 2020; pp. 1029–1053. [Google Scholar]
  7. Morin, E.; Marra, F.; Armon, M. Dryland precipitation climatology from satellite observations. In Satellite Precipitation Measurement; Springer: Cham, Switzerland, 2020; pp. 843–859. [Google Scholar]
  8. Kidd, C.; Shige, S.; Vila, D.; Tarnavsky, E.; Yamamoto, M.K.; Maggioni, V.; Maseko, B. The IPWG Satellite Precipitation Validation Effort The IPWG satellite precipitation validation effort. In Satellite Precipitation Measurement: Volume 2; Levizzani, V., Levizzani, V., Kidd, C., Kirschbaum, D.B., Kummerow, C.D., Nakamura, K., Turk, F.J., Eds.; Springer International Publishing: Berlin, Germany, 2020; pp. 453–470. [Google Scholar]
  9. Ebert, E.E.; Janowiak, J.E.; Kidd, C. Comparison of near-real-time precipitation estimates from satellite observations and numerical models. Bull. Am. Meteorol. Soc. 2007, 88, 47–64. [Google Scholar] [CrossRef] [Green Version]
  10. Massari, C.; Crow, W.; Brocca, L. An assessment of the performance of global rainfall estimates without ground-based observations. Hydrol. Earth Syst. Sci. 2017, 21, 4347–4361. [Google Scholar] [CrossRef] [Green Version]
  11. Alemohammad, S.H.; McColl, K.A.; Konings, A.G.; Entekhabi, D.; Stoffelen, A. Characterization of precipitation product errors across the United States using multiplicative triple collocation. Hydrol. Earth Syst. Sci. 2015, 19, 3489–3503. [Google Scholar] [CrossRef] [Green Version]
  12. De Freitas, C.R.; Helbig, M.; Matzarakis, A. Hydroclimatic assessment of water resources of low Pacific islands: Evaluating sensitivity to climatic change and variability. Int. J. Climatol. 2014, 34, 881–892. [Google Scholar] [CrossRef]
  13. CSIRO; Bureau of Meteorology. Climate Change in the Pacific: Scientific Assessment and New Research; Program, P.C.C.S., Ed.; Australian Government: Aspendale, Australia, 2011; p. 272. [Google Scholar]
  14. Kuleshov, Y.; McGree, S.; Jones, D.; Charles, A.; Cottrill, A.; Prakash, B.; Atalifo, T.; Nihmei, S.; Seuseu, F.L.S.K. Extreme weather and climate events and their impacts on island countries in the western Pacific: Cyclones, floods and droughts. Atmos. Clim. Sci. 2014, 4, 803–818. [Google Scholar] [CrossRef] [Green Version]
  15. Tahani, L.; Prakash, B.; Ngari, A.; Malala, H.; Mackenzie, B.; Inape, K.; Finaulahi, S.; Ene, E.; Kuleshov, Y.; Schreider, S.; et al. Recent changes in mean and extreme temperature and precipitation in the western Pacific islands. J. Clim. 2019, 32, 4919–4941. [Google Scholar] [CrossRef]
  16. McGree, S.; Schreider, S.; Kuleshov, Y. Trends and variability in droughts in the Pacific islands and northeast Australia. J. Clim. 2016, 29, 8377–8397. [Google Scholar] [CrossRef]
  17. Chen, F.; Crow, W.T.; Ciabatta, L.; Filippucci, P.; Panegrossi, G.; Marra, A.C.; Puca, S.; Massari, C. Enhanced large-scale validation of satellite-based land rainfall products. J. Hydrometeorol. 2021, 22, 245–257. [Google Scholar] [CrossRef]
  18. Li, C.; Tang, G.; Hong, Y. Cross-evaluation of ground-based, multi-satellite and reanalysis precipitation products: Applicability of the triple collocation method across mainland China. J. Hydrol. 2018, 562, 71–83. [Google Scholar] [CrossRef]
  19. Tanim, A.H.; Mullick, M.R.A.; Sikdar, M.S. Evaluation of spatial rainfall products in sparsely gauged region using copula uncertainty modeling with triple collocation. J. Hydrol. Eng. 2021, 26. [Google Scholar] [CrossRef]
  20. Chen, F.; Crow, W.T.; Ciabatta, L.; Filippucci, P.; Panegrossi, G.; Marra, A.C.; Puca, S.; Massari, C. Daily evaluation of 26 precipitation datasets using Stage-IV gauge-radar data for the CONUS. Hydrol. Earth Syst. Sci. 2019, 23, 207–224. [Google Scholar] [CrossRef] [Green Version]
  21. Brocca, L.; Filippucci, P.; Hahn, S.; Ciabatta, L.; Massari, C.; Camici, S.; Schüller, L.; Bojkov, B.; Wagner, W. SM2RAIN–ASCAT (2007–2018): Global daily satellite rainfall data from ASCAT soil moisture observations. Earth Syst. Sci. Data 2019, 11, 1583–1601. [Google Scholar] [CrossRef] [Green Version]
  22. Tashima, T.; Kubota, T.; Mega, T.; Ushio, T.; Oki, R. Precipitation extremes monitoring using the near-real-time GSMaP product. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens. 2020, 13, 5640–5651. [Google Scholar] [CrossRef]
  23. Huffman, G.J.; Bolvin, D.T.; Braithwaite, D.; Hsu, K.; Joyce, R.; Xie, P. Integrated Multi-Satellite Retrievals for GPM (IMERG), Algorithm Theoretical Basis Document; NASA: Washington, DC, USA, 2018. [Google Scholar]
  24. Kubota, T.; Aonashi, K.; Ushio, T.; Shige, S.; Takayabu, Y.N.; Kachi, M.; Arai, Y.; Tashima, T.; Masaki, T.; Kawamoto, N.; et al. Global satellite mapping of precipitation (GSMaP) products in the GPM Era. Satell. Precip. Meas. 2020, 1, 355–373. [Google Scholar]
  25. Lu, D.; Yong, B. Evaluation and hydrological utility of the latest GPM IMERG V5 and GSMaP V7 precipitation products over the Tibetan plateau. Remote Sens. 2018, 10, 2022. [Google Scholar] [CrossRef] [Green Version]
  26. Xie, P.; Joyce, R.; Wu, S.; Yoo, S.-H.; Yarosh, Y.; Sun, F.; Lin, R. Reprocessed, bias-corrected CMORPH global high-resolution precipitation estimates from 1998. J. Hydrometeorol. 2017, 18, 1617–1641. [Google Scholar] [CrossRef]
  27. Chua, Z.-W.; Kuleshov, Y.; Watkins, A. Evaluation of satellite precipitation estimates over Australia. Remote Sens. 2020, 12, 678. [Google Scholar] [CrossRef] [Green Version]
  28. Chua, Z.-W.; Kuleshov, Y.; Watkins, A.B. Drought detection over Papua New Guinea using satellite-derived products. Remote Sens. 2020, 12, 3589. [Google Scholar] [CrossRef]
  29. Kuleshov, Y.; Kurino, T.; Kubota, T.; Tashima, T.; Xie, P. WMO space-basedweather and climate extremes monitoring demonstration project (SEMDP): First outcomes of regional cooperation on drought and heavy precipitation monitoring for Australia and Southeast Asia. In RAINFALL—Extremes, Distribution and Properties; IntechOpen: London, UK, 2019; pp. 51–57. [Google Scholar]
  30. Huffman, G.J.; Bolvin, D.T.; Braithwaite, D.; Hsu, K.-L.; Joyce, R.J.; Kidd, C.; Nelkin, E.J.; Sorooshian, S.; Stocker, E.F.; Tan, J.; et al. Integrated multi-satellite retrievals for the global precipitation measurement (GPM) mission (IMERG). In Satellite Precipitation Measurement; Springer: Cham, Switzerland, 2020; pp. 343–353. [Google Scholar]
  31. Massari, C.; Maggioni, V. Error and uncertainty characterization. In Satellite Precipitation Measurement: Volume 2; Levizzani, V., Kidd, C., Kirschbaum, D.B., Kummerow, C.D., Nakamura, K., Turk, F.J., Eds.; Springer International Publishing: Cham, Switzerland, 2020; pp. 515–532. [Google Scholar]
  32. Beck, H.E.; van Dijk, A.I.J.M.; Levizzani, V.; Schellekens, J.; Miralles, D.G.; Martens, B.; de Roo, A. MSWEP: 3-hourly 0.25° global gridded precipitation (1979–2015) by merging gauge, satellite, and reanalysis data. Hydrol. Earth Syst. Sci. 2017, 21, 589–615. [Google Scholar] [CrossRef] [Green Version]
  33. Brocca, L.; Moramarco, T.; Melone, F.; Wagner, W. A new method for rainfall estimation through soil moisture observations. Geophys. Res. Lett. 2013, 40, 853–858. [Google Scholar] [CrossRef]
  34. Hersbach, H.; Bell, B.; Berrisford, P.; Hirahara, S.; Horányi, A.; Muñoz-Sabater, J.; Nicolas, J.; Peubey, C.; Radu, R.; Schepers, D.; et al. The ERA5 global reanalysis. Q. J. R. Meteorol. Soc. 2020, 146, 1999–2049. [Google Scholar] [CrossRef]
  35. Hersbach, H.; Bell, B.; Berrisford, P.; Hirahara, S.; Horányi, A.; Muñoz-Sabater, J.; Nicolas, J.; Peubey, C.; Radu, R.; Schepers, D.; et al. Estimating error cross-correlations in soil moisture data sets using extended collocation analysis. J. Geophys. Res. Atmos. 2016, 121, 1208–1219. [Google Scholar] [CrossRef] [Green Version]
  36. Hersbach, H.; Bell, B.; Berrisford, P.; Hirahara, S.; Horányi, A.; Muñoz-Sabater, J.; Nicolas, J.; Peubey, C.; Radu, R.; Schepers, D.; et al. Extended triple collocation: Estimating errors and correlation coefficients with respect to an unknown target. Geophys. Res. Lett. 2014, 41, 6229–6236. [Google Scholar] [CrossRef] [Green Version]
  37. Hersbach, H.; Bell, B.; Berrisford, P.; Hirahara, S.; Horányi, A.; Muñoz-Sabater, J.; Nicolas, J.; Peubey, C.; Radu, R.; Schepers, D.; et al. Tropical cyclone genesis in the Southern Hemisphere and its relationship with the ENSO. In Annales Geophysicae; Copernicus GmbH: Gottingen, Germany, 2009. [Google Scholar]
  38. Kuleshov, Y. Climate change and southern hemisphere tropical cyclones international initiative: Twenty years of successful regional cooperation. In Climate Change, Hazards and Adaptation Options; Springer: Dodrecht, The Netherlands, 2020; pp. 411–439. [Google Scholar]
  39. Khan, S.; Maggioni, V. Assessment of level-3 gridded global precipitation mission (GPM) products over oceans. Remote Sens. 2019, 11, 255. [Google Scholar] [CrossRef] [Green Version]
  40. Beck, H.E.; Vergopolan, N.; Pan, M.; Levizzani, V.; van Dijk, A.I.J.M.; Weedon, G.P.; Brocca, L.; Pappenberger, F.; Huffman, G.J.; Wood, E.F. Global-scale evaluation of 22 precipitation datasets using gauge observations and hydrological modeling. Hydrol. Earth Syst. Sci. 2017, 21, 6201–6217. [Google Scholar] [CrossRef] [Green Version]
  41. 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] [Green Version]
  42. Yang, S.; Chen, M.; Xie, P.; Yatagai, A.; Hayasaka, T.; Fukushima, Y.; Liu, C. A gauge-based analysis of daily precipitation over East Asia. J. Hydrometeorol. 2007, 8, 607–626. [Google Scholar] [CrossRef]
  43. Yang, S.; Chen, M.; Xie, P.; Yatagai, A.; Hayasaka, T.; Fukushima, Y.; Liu, C. So, how much of the Earth’s surface is covered by rain gauges? Bull. Am. Meteorol. Soc. 2017, 98, 69–78. [Google Scholar] [CrossRef]
  44. Masunaga, H.; Schröder, M.; Furuzawa, F.A.; Kummerow, C.; Rustemeier, E.; Schneider, U. Inter-product biases in global precipitation extremes. Environ. Res. Lett. 2019, 14. [Google Scholar] [CrossRef]
  45. Tarpanelli, A.; Massari, C.; Ciabatta, L.; Filippucci, P.; Amarnath, G.; Brocca, L. Exploiting a constellation of satellite soil moisture sensors for accurate rainfall estimation. Adv. Water Resour. 2017, 108, 249–255. [Google Scholar] [CrossRef]
  46. Tarpanelli, A.; Massari, C.; Ciabatta, L.; Filippucci, P.; Amarnath, G.; Brocca, L. Soil moisture and precipitation: The SM2RAIN algorithm for rainfall retrieval from satellite soil moisture. In Satellite Precipitation Measurement; Springer: Cham, Switzerland, 2020; pp. 1013–1027. [Google Scholar]
  47. Chew, C.C.; Small, E.E. Soil moisture sensing using spaceborne GNSS reflections: Comparison of CYGNSS reflectivity to SMAP soil moisture. Geophys. Res. Lett. 2018, 45, 4049–4057. [Google Scholar] [CrossRef] [Green Version]
  48. Chen, S.; Hong, Y.; Cao, Q.; Kirstetter, P.-E.; Gourley, J.J.; Qi, Y.; Zhang, J.; Howard, K.; Hu, J.; Wang, J. Performance evaluation of radar and satellite rainfalls for Typhoon Morakot over Taiwan: Are remote-sensing products ready for gauge denial scenario of extreme events? J. Hydrol. 2013, 506, 4–13. [Google Scholar] [CrossRef]
Figure 1. Example plot of precipitation anomaly for March 2007 depicting the study area, from 120° E to 180° E, 0° N to 30° S. The red boxes indicate the six studied sub-regions.
Figure 1. Example plot of precipitation anomaly for March 2007 depicting the study area, from 120° E to 180° E, 0° N to 30° S. The red boxes indicate the six studied sub-regions.
Remotesensing 13 03929 g001
Figure 2. Mean bias error for each SPE over each sub−region.
Figure 2. Mean bias error for each SPE over each sub−region.
Remotesensing 13 03929 g002aRemotesensing 13 03929 g002b
Figure 3. Comparison of (a) the mean MBE, (b) mean Pearson (R) Correlation, and (c) mean RMSE for each sub-region for the four SPEs.
Figure 3. Comparison of (a) the mean MBE, (b) mean Pearson (R) Correlation, and (c) mean RMSE for each sub-region for the four SPEs.
Remotesensing 13 03929 g003aRemotesensing 13 03929 g003b
Figure 4. Cross Correlations between each SPE and the unknown truth from triple CA of ERA5 and SM2RAIN−ASCAT. Smaller correlations can be seen around the mountainous regions, most evident for the sub-regions of PNG and the Solomon Islands.
Figure 4. Cross Correlations between each SPE and the unknown truth from triple CA of ERA5 and SM2RAIN−ASCAT. Smaller correlations can be seen around the mountainous regions, most evident for the sub-regions of PNG and the Solomon Islands.
Remotesensing 13 03929 g004
Figure 5. Comparison of (a) the mean CC and (b) RMSE using triple CA for each SPE.
Figure 5. Comparison of (a) the mean CC and (b) RMSE using triple CA for each SPE.
Remotesensing 13 03929 g005
Figure 6. Comparison of highest and second highest CC at each pixel for each sub-region.
Figure 6. Comparison of highest and second highest CC at each pixel for each sub-region.
Remotesensing 13 03929 g006
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wild, A.; Chua, Z.-W.; Kuleshov, Y. Evaluation of Satellite Precipitation Estimates over the South West Pacific Region. Remote Sens. 2021, 13, 3929. https://doi.org/10.3390/rs13193929

AMA Style

Wild A, Chua Z-W, Kuleshov Y. Evaluation of Satellite Precipitation Estimates over the South West Pacific Region. Remote Sensing. 2021; 13(19):3929. https://doi.org/10.3390/rs13193929

Chicago/Turabian Style

Wild, Ashley, Zhi-Weng Chua, and Yuriy Kuleshov. 2021. "Evaluation of Satellite Precipitation Estimates over the South West Pacific Region" Remote Sensing 13, no. 19: 3929. https://doi.org/10.3390/rs13193929

APA Style

Wild, A., Chua, Z. -W., & Kuleshov, Y. (2021). Evaluation of Satellite Precipitation Estimates over the South West Pacific Region. Remote Sensing, 13(19), 3929. https://doi.org/10.3390/rs13193929

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