Assessment of Gpm-imerg and Other Precipitation Products against Gauge Data under Different Topographic and Climatic Conditions in Iran: Preliminary Results

The new generation of weather observatory satellites, namely Global Precipitation Measurement (GPM) constellation satellites, is the lead observatory of the 10 highly advanced earth orbiting weather research satellites. Indeed, GPM is the first satellite that has been designed to measure light rain and snowfall, in addition to heavy tropical rainfall. This work compares the final run of the Integrated Multi-satellitE Retrievals for GPM (IMERG) product, the post real time of TRMM and Multi-satellite Precipitation Analysis (TMPA-3B42) and the Era-Interim product from the European Centre for Medium Range Weather Forecasts (ECMWF) against the Iran Meteorological Organization (IMO) daily precipitation measured by the synoptic rain-gauges over four regions with different topography and climate conditions in Iran. Assessment is implemented for a one-year period from March 2014 to February 2015. Overall, in daily scale the results reveal that all three products lead to underestimation but IMERG performs better than other products and underestimates precipitation slightly in all four regions. Based on monthly and seasonal scale, in Guilan all products, in Bushehr and Kermanshah ERA-Interim and in Tehran IMERG and ERA-Interim tend to underestimate. The correlation coefficient between IMERG and the rain-gauge data in daily scale is far superior to that of Era-Interim and TMPA-3B42. On the basis of daily timescale of bias in comparison with the ground data, the IMERG product far outperforms ERA-Interim and 3B42 products. According to the categorical verification technique in this study, IMERG yields better results for detection of precipitation events on the basis of Probability of Detection (POD), Critical Success Index (CSI) and False Alarm Ratio (FAR) in those areas with stratiform and orographic precipitation, such as Tehran and Kermanshah, compared with other satellite/model data sets. In particular, for heavy precipitation (>15 mm/day), IMERG is superior to the other products in all study areas and could be used in future for meteorological and hydrological models, etc.


Introduction
Water is fundamental to life on Earth.Knowing the locations and extents of rain and snowfall is vital to understanding how weather and climate impact our environment and Earth's water and energy cycles, including effects on agriculture, fresh water availability, and response to natural disasters.Since rainfall and snowfall vary greatly on a small scale and over time, satellites can provide a better picture Researchers and scientists use models for analysis and forecasts of the atmospheric state.The models may have different spatial and temporal resolution.One of the best models, which is used around the world, has been developed by the European Center for Medium Range Weather Forecasts (ECMWF).They have produced a global reanalysis for the last decades of ERA-40 and ERA-Interim, 12-hourly and daily precipitation fields [8].Moreover, ERA-Interim is the latest global atmospheric reanalysis produced by ECMWF.The ERA-Interim project was conducted in part to prepare a new atmospheric reanalysis to replace ERA-40, which will extend the data to the early part of the twentieth century.Another tool to determine precipitation at relatively fine temporal and spatial scale is satellite observation.During the last decade, many researchers have evaluated and used satellite data sets and some have found that the TMPA-3B42 post real-time product performed better than other products compared with the rain-gauges [9][10][11][12][13][14][15].One of the disadvantages of TMPA-3B42 was the limited area of observation; it covered only the tropical and subtropical belts.Another disadvantage was the ability to estimate heavy rainfall while light rainfall and snowfall were not detected properly [6][7][8][9][10][11][12].
As seen in Table 1, IMERG is presently available from mid-March 2014 to the present (with access delay in the order of about three months for the final run version).Based on the preliminary analysis during beta testing by Huffman et al. (2015), IMERG is smoother than 3B43 over oceans and at higher latitudes.This was a goal for IMERG, which provides estimates every half hour, versus the 3-hourly interval for the satellite data contributing to 3B43.Usually, satellite estimation in regions subject to convective precipitation could be difficult whereas they expect that IMERG data sets will be more accurate [7].Contrary to other satellites, such as TRMM, that could not measure light rain and snowfall, GPM-IMERG uses different sensors from different satellites to detect both light and heavy rain and snowfall.Three critical improvements in GPM are that (1) the orbital inclination has been increased from 35 ˝to 65 ˝, affording coverage of important additional climate zones; (2) the radar has been upgraded to two frequencies, adding sensitivity to light precipitation; and (3) high-frequency channels (165.5 and 183.3 GHz) have been added to the passive microwave (PMW) imager, which are expected to facilitate sensing of light and solid precipitation.
In brief, the input precipitation estimates computed from the various satellite passive microwave sensors are inter-calibrated to the GPM Combined Instrument (GCI, using GMI and DPR), because it is presumed to be the best snapshot GPM estimate, then morphed and combined with microwave precipitation-calibrated geosynchronous earth orbit (geo) infrared (IR) fields, and adjusted with monthly surface precipitation gauge analysis data (where available) to provide half-hourly and monthly precipitation estimates.Precipitation phase is diagnosed using analyses of surface temperature, humidity, and pressure.On the other hand, the TMPA combines microwave data from multiple satellites, each inter-calibrated to the TRMM Combined Instrument (TCI), using TRMM Microwave Imager (TMI) and TRMM Precipitation Radar (PR).Coverage gaps in space and time are filled in with calibrated infrared (IR) data (which are generally available with near-global coverage every 3 h); coefficients are derived from co-located IR brightness temperatures and the microwave-based precipitation estimates.The final data products reflect scaling the multi-satellite estimates to rain gauge data on a monthly basis, and ensuring that the 3-hourly averages in 3B42 sum to the monthly totals in 3B43 [7].
Other advantages for using satellite precipitation data in places such as Iran include generally insufficient spatial coverage of in-situ data, long delay in data processing and transfer until they become accessible for the public and scientific use, and absence of data sharing in many trans-boundary basins.As a result, this study aims to assess the accuracy of the new generation of satellite precipitation products.To our knowledge, no report exists yet to study the GPM constellation satellites data over Iran.

Recent Works
The gauges potentially provide in-situ, high time-resolution precipitation measurements at an area of a few square meters and also they have error due to wind, solid precipitation, etc. while the TMI and PR instruments essentially offer instantaneous area-averaged measurements with relatively poor time sampling.The measurement errors of both the satellite retrievals and the gauges are not fully understood [16,17].The comparison of satellite estimates with rain-gauge measurements is often frustrating, because satellites can at best attempt to measure rain amounts over areas many kilometers in size around a gauge, whereas rain-gauges can only record what falls in an area some tens of centimeters in diameter [18].However, recently released IMERG is going to reduce this sampling issue such that GPM constellation satellites will have improved spatial and temporal resolutions (i.e., 30 min).
Many attempts have been made to develop and improve global and regional gridded precipitation gauge data sets in recent years [19][20][21].However, due to the lack of observational precipitation datasets over land areas for parts of Asia, Yatagai et al. [22,23] generated a high resolution rain-gauge based daily precipitation grid dataset for the East, Middle East and Russia.
Xie et al. [24] assessed the performance of five satellite-based rainfall outputs by comparison against gauge analysis: TRMM (3B42-V7, 3B42-RT) of US National Weather Service (NWS), Climate Prediction Center (CPC), Morphing Techniques (CMORPH), Precipitation Estimation from Remotely Sensed Information using Artificial Neural Networks (PERSIANN), and those of Naval Research Laboratory (NRL).They found that all satellite products performed better in wet climate regions and wet seasons, but showed limited skills in estimating precipitation over central Asian arid and semi-arid regions.
Renzullo et al. [25] constructed daily precipitation surfaces for Australia from the post real-time TRMM 3B42 rainfall product.This product has been used in stream flow and flood modeling studies where real-time gauge data is sparse [26,27].Fleming et al. [28] indicated that there is high correlation between TRMM rainfall data and the Australian Bureau of Meteorology monthly grid values data.Regarding the correlation between datasets as a function of climate, the steppe and temperate climate zones show a higher correlation than the tropical and desert zones.The quality of the analysis depends on the density of gauge reports (sparse data imply larger errors), the complexity of the local terrain, i.e., gauges in mountain areas are generally sited in valleys and, therefore, underreport the true areal average rainfall [29].
A strong increase in tropical land precipitation during the past decade is found in ERA-Interim and Global Precipitation Climatology Project 2. While large improvements have been made in terms of temporal homogeneity compared with ERA-40, especially with regard to precipitation, the homogeneity issue still needs to be carefully addressed when interpreting surface flux data from reanalysis or other sources [30].
The performance of GPM satellite-based precipitation has not been studied over arid and semi-arid areas in Iran.To address this, in our research, we compare precipitation estimates of GPM-IMERG-final run data sets with meteorological synoptic stations' data over Iran to evaluate the GPM-based precipitation outputs.

Study Area
To examine the accuracy of satellite observatory precipitation data in different topographic conditions, we selected four study areas in Iran with diverse topography and within different precipitation zones proposed by Modarres [31] (see Figure 1a,b).
Due to Iran's variable topography and climate conditions, we partitioned the study area to cover provinces such as Kermanshah situated in western Iran.Most of the surface of the province lies within the Zagros Mountains, which forms the western periphery of the Iranian Plateau.Running from southeast to northwest, the nearly parallel broken ridges of the Zagros are highest in the east of the province, and elevation drops progressively towards the west, until the vast plains of Iraq fill the horizon (shown in Figure 1 in zone G5).The climate of the highlands is mild in summer and cold in winter, with heavy snowfall; only the province's western strip belongs to the warm climate.The average temperatures in Kermanshah City are approximately 0 ˝C in January and 26 ˝C in July.The winds blowing from the Mediterranean Sea carry rainclouds, with an annual precipitation of up to 700 mm in the highlands and about 400 mm at Kermanshah City [32].The second province, namely Guilan, belongs to zone G8 and includes the northwestern end of the Alborz Mountains and the western part of the Caspian lowlands of Iran.The mountainous belt is cut through by the deep transversal valley of the Sefidrud.To the northwest, the highlands stretch over a continuous watershed and, except for at their northern end, they are over 2000 m high, with three spots over 3000 m.Their eastern and northeastern side is deeply carved by parallel streams flowing down towards the Caspian Sea.Prevailing north-south atmospheric currents, humidified over the Caspian, are forced to a vigorous ascent by the mighty barrier of Alborz and thus are responsible for abundant rainfall on both the plain and the northwestern slope of the mountains.Mean annual precipitation varies between 1200 mm and 1800 mm along the shoreline, decreases towards a sub-humid area in the southwestern corner of the plain, and reaches again very high amounts in the lower part of the mountain, up to 1500-2400 mm including convective rainfall during summer.Bushehr (shown in Figure 1 in zone G4) is just northwest of the Persian Gulf with an arid and subtropical climate.There are two different regions in Bushehr, which are the plains along the Persian Gulf, and the mountain ranges including the Zagros Mountains.Rainfall in autumn and spring usually occurs in the form of convective precipitation, and during the winter in the form of cyclonic storms.The rest of precipitation is in the form of moderate and light rain.Snowfall was unprecedented in this region.The average annual rainfall in Bushehr is approximately 230 mm and the average temperatures are approximately 24 ˝C.Tehran (shown in Figure 1 in zone G2), the capital, at the southern slopes of Alborz Mountains with a semi-arid climate, represents the remaining study area.Tehran's climate is largely defined by its geographic location, with the towering Alborz Mountains to its north and the central desert to the south.It can be generally described as mild in spring and autumn, hot and dry in summer, and cold in winter.The highest point of the Tehran province is Damavand peak at an altitude of 5678 m and the lower most area of the province being the plains of Varamin, 790 m above sea level and located to the south-east of the province.The hottest months of the year are from mid-July to mid-September when temperatures average around 28-30 ˝C and the coldest months drop to 1 ˝C temperature in December to January, and the average annual precipitation is approximately 400 mm, the maximum being during the winter season (see Figure 2a,b).

Data Sources
In this study, we used only the in-situ meteorological synoptic stations within the four study areas that have undergone quality control by the Iran Meteorological Organization (IMO).However, the stations with large data gaps were removed.The GPM mission's Precipitation Processing System (PPS) at NASA's Goddard Space Flight Center released the IMERG data to the public in late February, 2015.The data set includes precipitation rates since mid-March 2014.Current and future data sets are freely available to users from NASA Goddard Earth Sciences Data and Information Services Center (GES DISC) website [33].At present, the IMERG data are available from 12 March 2014 to present, thus providing IMERG-V03D data for March 2014 to February 2015 period.The 3B42-V7 from TRMM and other satellites at GES DISC and ERA-Interim data sets from ECMWF [34] are available to the public and are free of charge.
We should mention that further works are needed to evaluate the interannual variation of IMERG data relying on larger samples.
The satellite products were compared with the measured precipitation data of 43 meteorological synoptic stations (Figure 3) at daily, monthly and seasonal time scales for different climate conditions in Iran.At the time of data analysis, we had access only to this number of stations.
In this research, the data of the nearest grid point in the satellite grid (less than 0.05 degrees off the grid points) is compared with that corresponding to the ground point observation (i.e., the meteorological synoptic station).In cases where the grid points were close to the stations, the comparison was carried out directly between them.However, in cases where the ground station was surrounded by four grid cells but not particularly close to none, an average of the four grid points around the station was used as the basis for comparison.Table 1 indicates the characteristics of data sets and Table 2 shows the number of stations and pixels in each area as used in this study.

Data Sources
In this study, we used only the in-situ meteorological synoptic stations within the four study areas that have undergone quality control by the Iran Meteorological Organization (IMO).However, the stations with large data gaps were removed.The GPM mission's Precipitation Processing System (PPS) at NASA's Goddard Space Flight Center released the IMERG data to the public in late February, 2015.The data set includes precipitation rates since mid-March 2014.Current and future data sets are freely available to users from NASA Goddard Earth Sciences Data and Information Services Center (GES DISC) website [33].At present, the IMERG data are available from 12 March 2014 to present, thus providing IMERG-V03D data for March 2014 to February 2015 period.The 3B42-V7 from TRMM and other satellites at GES DISC and ERA-Interim data sets from ECMWF [34] are available to the public and are free of charge.
We should mention that further works are needed to evaluate the interannual variation of IMERG data relying on larger samples.
The satellite products were compared with the measured precipitation data of 43 meteorological synoptic stations (Figure 3) at daily, monthly and seasonal time scales for different climate conditions in Iran.At the time of data analysis, we had access only to this number of stations.
In this research, the data of the nearest grid point in the satellite grid (less than 0.05 degrees off the grid points) is compared with that corresponding to the ground point observation (i.e., the meteorological synoptic station).In cases where the grid points were close to the stations, the comparison was carried out directly between them.However, in cases where the ground station was surrounded by four grid cells but not particularly close to none, an average of the four grid points around the station was used as the basis for comparison.Table 1 indicates the characteristics of data sets and Table 2 shows the number of stations and pixels in each area as used in this study.

Statistical Analysis
A comprehensive study was performed on the correlation of IMERG, 3B42, and Era-Interim with the station precipitation data over Iran.In the first step, statistical indices such as bias, multiplicative bias (MBias), relative bias (RBias), mean absolute error (MAE), root mean square error (RMSE), and linear correlation coefficient (CC) were determined.The bias is defined as the average difference between in-situ observations and satellite/model precipitation estimates (S/M-PE), and can be either positive or negative.A negative bias indicates underestimation by satellite precipitation while a positive bias indicates overestimation.
The multiplicative bias (MBias) is the ratio of S/M-PE over rain-gauge value such that a perfect estimation would result in an Mbias of unity.Underestimation will lead to values less than unity, and overestimation to values greater than unity.
The relative bias (RBias) describes the systematic bias of satellite-based precipitation and behaves the same as bias.

Statistical Analysis
A comprehensive study was performed on the correlation of IMERG, 3B42, and Era-Interim with the station precipitation data over Iran.In the first step, statistical indices such as bias, multiplicative bias (MBias), relative bias (RBias), mean absolute error (MAE), root mean square error (RMSE), and linear correlation coefficient (CC) were determined.The bias is defined as the average difference between in-situ observations and satellite/model precipitation estimates (S/M-PE), and can be either positive or negative.A negative bias indicates underestimation by satellite precipitation while a positive bias indicates overestimation.
Bias " The multiplicative bias (MBias) is the ratio of S/M-PE over rain-gauge value such that a perfect estimation would result in an Mbias of unity.Underestimation will lead to values less than unity, and overestimation to values greater than unity.
MBias " The relative bias (RBias) describes the systematic bias of satellite-based precipitation and behaves the same as bias.
RBias " The mean absolute error (MAE) is used to represent the average magnitude of the error.
MAE " The root mean square error (RMSE), which gives a greater weight to the larger errors relative to MAE, is used to measure the average error magnitude.
The correlation coefficient (CC) is used to assess the agreement between satellite/model-based precipitation and rain-gauge observations.The value of CC is such that ´1 < CC < +1.A CC value of +1 indicates a perfect positive fit.If there is no linear correlation or a weak linear correlation, CC is close to zero.

CC "
where P S i is the value of satellite/model precipitation estimates for the ith daily event, P O i is the value of rain-gauge observation for the ith daily event, N is the number of observed days, P S is the average value of satellite/model precipitation estimates for N observed days over each pixel or grid points.

Categorical Technique
Another assessment technique of satellite estimation/model forecast is using a contingency table that reflects the frequency of "Yes" and "No" of the satellite estimation/forecast model (see Table 3).A dichotomous estimates says, "Yes, an event will happen", or "No, the event will not happen".By using this table for daily precipitation, a set of statistical indices are shown as follows: Probability of detection (POD) responds to the question of what fraction of the observed "Yes" events was correctly estimated/forecasted.The perfect score is 1.
False alarm ratio (FAR) deals with the question of what fraction of the estimated/forecasted "Yes" events did not occur.The ideal score is 0.
Critical success index (CSI) or threat score (TS), answers the question of how well the estimated/forecasted "Yes" events corresponded to the observed "Yes" events.The perfect score is 1.
Accuracy (fraction correct) measures the fraction of correct estimates/forecasts and its perfect score is 1.
Bias (frequency bias) answers the question of how the estimated/forecasted frequency of "Yes" events compared to the observed frequency of "Yes" events.Range of values is 0 to 8 with a perfect score of 1.
Probability of false detection (POFD) deals with the question of what fraction of observed "No" events were incorrectly estimated/forecasted as "Yes".The range varies from 0 to 1 and the perfect score is 0.
Success ratio (SR) responds to the question of what fraction of estimated/forecasted "Yes" events was correctly observed.The range is 0 to 1 and the perfect score is 1.
Equitable threat score (ETS) or Gilbert skill score, answers the question of how well the estimated/forecasted "Yes" events corresponded to the observed "Yes" events.The range is ´1/3 to 1 and the perfect score is 1.For rare events, the minimum ETS value is near 0, while the absolute minimum is obtained if the event has a climatological frequency of 0.5, and there are no hits.If the score goes below 0 then the chance forecast is preferred to the actual forecast, and the forecast is said to be unskilled.
ETS " a ´arandom a `b `c ´arandom (14) Odds ratio (OR) deals with the ratio of the odds of "yes" estimates/forecasts being correct over the odds of "Yes" estimates/forecasts being wrong.Odds ratio range is 0 to 8, 0 indicates no skill and the perfect score is 8.
Hanssen and Kuiper discriminant (HK) or true skill statistic (TSS) covers the question of how well the estimates/forecast separated the "Yes" events from the "No" events.The range is ´1 to 1, while 0 indicates no skill and 1 is the perfect score [35][36][37].
ORSS " a ˚d ´c ˚b a ˚d `c ˚b (17) where a represents the number of times that observed rain is correctly detected, b is the number of times that rain is detected but not observed, c is the number of times that observed rain is not detected, d is the number of times that observed and estimated did not occur and total is the sample size and a random " pa `cq ˚pa `bq total (18)

Results and Discussion
This Study for the first time has evaluated daily, monthly and seasonally precipitation estimates (derived from half-hourly precipitation estimates of GPM constellation satellites (IMERG) and 3-hourly precipitation estimates of TRMM and multi satellite precipitation analysis (3B42)) and ERA-Interim precipitation fields from ECMWF in comparison with the corresponding rain-gauge observations according to four different topographic and climatic conditions over parts of Iran.The evaluation results follow.

Daily Evaluation
We start with the evaluation of daily IMERG, TMPA-3B42 and ERA-Interim data against rain-gauge observations.This is based on days with observed precipitation solely.Figures 4-6 and Figure 7a-c show the scatter plots of average daily precipitation values during March 2014 to February 2015 in four regions.In Figure 4 and Table 4, all three satellite/model products indicate underestimation according to Bias, Mbias and Rbias values but IMERG shows slight underestimation which means that all three products generate negative systematic errors.3B42 and ERA-Interim heavily underestimate precipitation in the Guilan region (G8) with respect to Mbias and Rbias.This is more so when the precipitation is above 20 mm according to Figure 4c for ERA-Interim.On the basis of RMSE and MAE, ERA-Interim shows better results in comparison with the IMERG and 3B42.evaluation results follow.

Daily Evaluation
We start with the evaluation of daily IMERG, TMPA-3B42 and ERA-Interim data against raingauge observations.This is based on days with observed precipitation solely.Figures 4-6 and 7a-c show the scatter plots of average daily precipitation values during March 2014 to February 2015 in four regions.In Figure 4 and Table 4, all three satellite/model products indicate underestimation according to Bias, Mbias and Rbias values but IMERG shows slight underestimation which means that all three products generate negative systematic errors.3B42 and ERA-Interim heavily underestimate precipitation in the Guilan region (G8) with respect to Mbias and Rbias.This is more so when the precipitation is above 20 mm according to Figure 4c for ERA-Interim.On the basis of RMSE and MAE, ERA-Interim shows better results in comparison with the IMERG and 3B42.
As seen in Figure 5 and Table 5, IMERG, 3B42 and ERA-Interim tend to underestimate while ERA-Interim underestimates heavily when precipitation is above 10 mm in Bushehr (G4).The IMERG yields Bias, Mbias and Rbias better than other products.On the other hand, ERA-Interim provides better values of RMSE and MAE in comparison with 3B42 and IMERG.In all, IMERG indicates a good agreement with the observed data.
Figure 6 presents the scatter plot of daily precipitation over Kermanshah area (G5) where IMERG, 3B42 and ERA-Interim underestimate the observed values.As seen in Table 6, IMERG has a better value for RMSE, Bias, Mbias, Rbias and CC in comparison with other data sets while IMERG performs better when precipitation is greater than 10 mm.
As seen in Figure 7, for, IMERG, 3B42 and Era-Interim underestimate Tehran (G2) daily precipitation.Furthermore, based on Table 7, 3B42 and IMERG show lower bias and better value of Mbias and Rbias in comparison with the other two datasets whereas ERA-Interim yields better in RMSE and MAE.As seen in Figure 5 and Table 5, IMERG, 3B42 and ERA-Interim tend to underestimate while ERA-Interim underestimates heavily when precipitation is above 10 mm in Bushehr (G4).The IMERG yields Bias, Mbias and Rbias better than other products.On the other hand, ERA-Interim provides better values of RMSE and MAE in comparison with 3B42 and IMERG.In all, IMERG indicates a good agreement with the observed data.As seen in Figure 7, for, IMERG, 3B42 and Era-Interim underestimate Tehran (G2) daily precipitation.Furthermore, based on Table 7, 3B42 and IMERG show lower bias and better value of Mbias and Rbias in comparison with the other two datasets whereas ERA-Interim yields better in RMSE and MAE.

Monthly Evaluation
To evaluate monthly precipitation, we derived IMERG from half-hourly, 3B42 from 3-hourly and ERA-Interim from daily precipitation.As a reference we compared rain-gauge observations with their corresponding grid points on satellites and models for all days of the corresponding months.Figure 8 and Table 8 show underestimation for monthly scale in all three products but very slightly for IMERG with Bias, Rbias and Mbias of ´29.18 mm, 0.98 and ´0.02, respectively, which shows a better result in comparison with the other two products.Also, IMERG yields MAE and RMSE of 52.75 mm and 81.43 mm, respectively, in comparison to 59 mm and 102.76 mm for 3B42 while the same for ERA-Interim are 62.18 mm and 95.41 mm.On the other hand, CC of ERA-Interim indicates better agreement than others in the Guilan region.
As seen in Figure 9 and Table 9, IMERG and 3B42 tend to overestimate while ERA-Interim underestimates in Bushehr.The 3B42 yields for all indices better results than the other two data sets while ERA-Interim shows significant agreement with a correlation coefficient of 0.95.

Monthly Evaluation
To evaluate monthly precipitation, we derived IMERG from half-hourly, 3B42 from 3-hourly and ERA-Interim from daily precipitation.As a reference we compared rain-gauge observations with their corresponding grid points on satellites and models for all days of the corresponding months.Figure 8 and Table 8 show underestimation for monthly scale in all three products but very slightly for IMERG with Bias, Rbias and Mbias of −29.18 mm, 0.98 and −0.02, respectively, which shows a better result in comparison with the other two products.Also, IMERG yields MAE and RMSE of 52.75 mm and 81.43 mm, respectively, in comparison to 59 mm and 102.76 mm for 3B42 while the same for ERA-Interim are 62.18 mm and 95.41 mm.On the other hand, CC of ERA-Interim indicates better agreement than others in the Guilan region.
As seen in Figure 9 and Table 9, IMERG and 3B42 tend to overestimate while ERA-Interim underestimates in Bushehr.The 3B42 yields for all indices better results than the other two data sets while ERA-Interim shows significant agreement with a correlation coefficient of 0.95.
Regarding to monthly precipitation over the Kermanshah area, based on Figure 10 and Table 10, IMERG and 3B42 overestimate and ERA-Interim underestimates the observed values.3B42 presents lower Bias, MAE and RMSE in comparison with the other two data sets.Furthermore, based on Figure 10, ERA-Interim shows a very slight bias during the January, February and March (winter season).
As seen in Figure 11 and Table 11, 3B42 tends to overestimate Tehran monthly precipitation whereas IMERG and ERA-Interim tend to underestimate.IMERG shows strong overestimates during October while ERA-Interim is almost unbiased in this month.Generally, IMERG represents better value of Bias, Mbias and Rbias while Era-Interim yields a better value for RMSE and MAE.As seen in Figure 11 and Table 11, 3B42 tends to overestimate Tehran monthly precipitation whereas IMERG and ERA-Interim tend to underestimate.IMERG shows strong overestimates during October while ERA-Interim is almost unbiased in this month.Generally, IMERG represents better value of Bias, Mbias and Rbias while Era-Interim yields a better value for RMSE and MAE.

Seasonal Evaluation
To evaluate the seasonal precipitation rain-gauge observations with their corresponding grid points on satellites and models, all days of the corresponding seasons are compared.Based on seasonal scale, Figure 12 and Table 12 indicate that all three products tend to underestimate in all seasons except spring in the Guilan region.IMERG has a better value of Bias, RMSE, Mbias and Rbias while ERA-Interim represents significant agreement.
With respect to Figure 13 and Table 13, 3B42 shows the best results in Bushehr in comparison with other data sets.IMERG and 3B42 displayed overestimates and ERA-Interim underestimates.
As seen in Figure 14 and Table 14, 3B42 and IMERG tend to overestimate while ERA-Interim underestimates in the Kermanshah region.The IMERG yields a correlation coefficient of 0.95 which is a good agreement between satellite estimation and rain-gauge observation and there is no great difference for MAE and Bias among all data sets.In all, ERA-Interim indicates a strong bias during

Seasonal Evaluation
To evaluate the seasonal precipitation rain-gauge observations with their corresponding grid points on satellites and models, all days of the corresponding seasons are compared.Based on seasonal scale, Figure 12 and Table 12 indicate that all three products tend to underestimate in all seasons except spring in the Guilan region.IMERG has a better value of Bias, RMSE, Mbias and Rbias while ERA-Interim represents significant agreement.
With respect to Figure 13 and Table 13, 3B42 shows the best results in Bushehr in comparison with other data sets.IMERG and 3B42 displayed overestimates and ERA-Interim underestimates.
As seen in Figure 14 and Table 14, 3B42 and IMERG tend to overestimate while ERA-Interim underestimates in the Kermanshah region.The IMERG yields a correlation coefficient of 0.95 which is a good agreement between satellite estimation and rain-gauge observation and there is no great difference for MAE and Bias among all data sets.In all, ERA-Interim indicates a strong bias during autumn and a slight bias during winter.IMERG and ERA-Interim displayed underestimation and 3B42 overestimation over Tehran region.Among these products, IMERG has better results in comparison with the other data sets while ERA-Interim yields a smaller magnitude of error (Figure 15 and Table 15).IMERG and ERA-Interim displayed underestimation and 3B42 overestimation over Tehran region.Among these products, IMERG has better results in comparison with the other data sets while ERA-Interim yields a smaller magnitude of error (Figure 15 and Table 15).

Evaluation of Dichotomous Estimates/Forecasts
With respect to an evaluation through contingency tables, Figure 16a-c and Table 16 reveals

Evaluation of Dichotomous Estimates/Forecasts
With respect to an evaluation through contingency tables, Figure 16a-c and Table 16 reveals higher POD (0.74), lower FAR (0.40) and higher CSI (0.49) by ERA-Interim over the North of Iran, Guilan area, compared to other satellite products.High POD may have been influenced by the dominance of convective storms [38].Lower POD for 3B42 and IMERG may be associated with missed precipitation over this region.In Kermanshah area, IMERG yields a higher POD and CSI and lower FAR in comparison to 3B42.These results could be due to the inability of TRMM Microwave Imager (TMI) to measure snowfall and light rain over the ground and Zagros Mountains while the main feature of Dualfrequency Precipitation Radar (DPR) based on the GPM is to measure snowfall and light rain.Overall, the ERA-Interim product produced fewer robust results when compared to IMERG.
As Figure 16a-c and Table 16 indicate, IMERG product performs better than others with higher POD and CSI and lower FAR in Tehran.We should note that precipitation in Tehran is stratiform and orographic.The missed precipitation may be caused by snow cover on the ground at higher altitudes over the Alborz Mountains and by the inability to catch warm rain processes, short-lived convective storms, or maritime precipitation along the Caspian Sea.POD value of IMERG (0.70) seems to be better than other products over Bushehr area while it is quite low for ERA-Interim (0.39).On the other hand, IMERG indicates a high FAR value (0.59) because of rare precipitation events despite high humidity.Figure 16b shows similar results for 3B42 and ERA Interim in Bushehr area as well.
In Kermanshah area, IMERG yields a higher POD and CSI and lower FAR in comparison to 3B42.These results could be due to the inability of TRMM Microwave Imager (TMI) to measure snowfall and light rain over the ground and Zagros Mountains while the main feature of Dual-frequency Precipitation Radar (DPR) based on the GPM is to measure snowfall and light rain.Overall, the ERA-Interim product produced fewer robust results when compared to IMERG.
As Figure 16a-c and Table 16 indicate, IMERG product performs better than others with higher POD and CSI and lower FAR in Tehran.We should note that precipitation in Tehran is stratiform and orographic.

Evaluation of Dichotomous Estimates/Forecasts for Precipitation below 15 mm
As can be seen in Figure 17a-c     In contrast, for precipitation above 15 mm/day, IMERG indicates better results for POD and CSI in all case studies while ERA-Interim indicates a very weak value (Figure 18a-c and Table 18) that could be due to KuPR (Ku-band (13.6 GHz) precipitation radar) instrument that can detect heavy precipitation.
Notice that almost in all study areas the value of FAR for precipitation above 15 mm/day from all products was rather high which could be due to observed liquid water in the atmosphere profiles which did not fall as precipitation due to their small size or evaporation.precipitation on the daily scale.With respect to contingency table metrics, IMERG outperformed ERA-Interim on the basis of POD, FAR, RMSE, Bias, Mbias, Rbias and CC (correlation coefficient) while ERA-Interim performed best in terms of MAE values on the daily scale.3B42 outperformed other products on a monthly scale.Moreover, all three products showed rather the same behavior in seasonal scale but IMERG indicates a better CC.

3.
Just northwest of the Persian Gulf, Bushehr region is subject to warm and subtropical arid climate in low latitudes.In this area, on the daily scale, IMERG was more accurate in terms of Bias, Mbias, Rbias and CC while 3B42 was better in RMSE and MAE.Also, 3B42 outperformed other products on monthly and seasonal scales.On the other hand, IMERG showed reasonable results with a POD of 0.70 that could be due to the dual-frequency sensor based on the GPM which can detect light rain.4.
Tehran is located in a semi-arid climate with the towering Alborz Mountains to its North and the central desert to the South.In this region, according to the daily scale precipitation, IMERG outperformed other precipitation products with POD, FAR and CSI. 5.
ERA-Interim yields weak results of POD, FAR and CSI for precipitation above 15 mm/day over Iran while IMERG is far superior to the other products in all study areas.ERA-Interim in Guilan region shows a significant value of POD for precipitation below 15 mm/day.6.
Generally, in Guilan, Kermanshah, Tehran and Bushehr regions, all three products (IMERG, 3B42 and ERA-Interim) underestimate precipitation based on daily scale, which might be due to an inadequate number of gauges which are provided by the Global Precipitation Climatology Centre (GPCC) and used for bias correction in satellite products or/and the inability to measure available water in the air profiles.Moreover, in semi-arid and hot climates, rain drops may evaporate before reaching the ground [40].7.
Based on monthly and seasonally scale, in Guilan, all products tend to underestimate while in Bushehr and Kermanshah, IMERG and 3B42, and in Tehran, 3B42 indicated overestimate.
In summary, the overestimation or underestimation over mountainous regions indicates that accurate estimation by satellite-based precipitation products remains a challenge.Such inaccuracy may be rooted in the inadequate number of gauges, provided by the Global Precipitation Climatology Centre (GPCC) and used for bias correction in satellite products or the non-uniform beam filling (NUBF) problem for remote sensing instruments.Another cause might be the spatial resolution of the satellite product, since precipitation within a region may occur in smaller scales than the pixel size of satellites.Overall, this preliminary accuracy assessment highlights that the IMERG product can adequately substitute 3B42 products despite its limited historic record.As more IMERG data become available, more detailed studies of GPM-IMERG applications in water, weather, and climate studies are possible in the near future.We expect that the regional analysis of GPM constellation satellites-based precipitation estimates reported here can give the satellite precipitation users a better understanding of the features associated with currently available IMERG precipitation estimates from a broader perspective.
eastern and northeastern side is deeply carved by parallel streams flowing down towards the Caspian Sea.Prevailing north-south atmospheric currents, humidified over the Caspian, are forced to a vigorous ascent by the mighty barrier of Alborz and thus are responsible for abundant rainfall on both the plain and the northwestern slope of the mountains.Mean annual precipitation varies between 1200 mm and 1800 mm along the shoreline, decreases towards a sub-humid area in the southwestern corner of the plain, and reaches again very high amounts in the lower part of the mountain, up to 1500-2400 mm including convective rainfall during summer.Bushehr (shown in Figure1in zone G4) is just northwest of the Persian Gulf with an arid and subtropical climate.There are two different regions in Bushehr, which are the plains along the Persian Gulf, and the mountain ranges including the Zagros Mountains.Rainfall in autumn and spring usually occurs in the form of convective precipitation, and during the winter in the form of cyclonic storms.The rest of precipitation is in the form of moderate and light rain.Snowfall was unprecedented in this region.The average annual rainfall in Bushehr is approximately 230 mm and the average temperatures are approximately 24 °C.Tehran (shown in Figure1in zone G2), the capital, at the southern slopes of Alborz Mountains with a semi-arid climate, represents the remaining study area.Tehran's climate is largely defined by its geographic location, with the towering Alborz Mountains to its north and the central desert to the south.It can be generally described as mild in spring and autumn, hot and dry in summer, and cold in winter.The highest point of the Tehran province is Damavand peak at an altitude of 5678 m and the lower most area of the province being the plains of Varamin, 790 m above sea level and located to the south-east of the province.The hottest months of the year are from mid-July to mid-September when temperatures average around 28-30 °C and the coldest months drop to 1 °C temperature in December to January, and the average annual precipitation is approximately 400 mm, the maximum being during the winter season (see Figure2a,b).

Figure 1 .
Figure 1.(a) Spatial distribution of precipitation zones of Iran [31]; and (b) Provincial map of Iran and the four selected study regions.

Figure 1 .
Figure 1.(a) Spatial distribution of precipitation zones of Iran [31]; and (b) Provincial map of Iran and the four selected study regions.

Table 4 .
Statistical indices of daily precipitation events in Guilan area (mid-March 2014 to February 2015).

Figure 16 .
Figure 16.(a-c) POD, FAR and CSI in the four study regions.

Figure 16 .
Figure 16.(a-c) POD, FAR and CSI in the four study regions.

Figure 17 .
Figure 17.(a-c) POD, FAR and CSI in four study region for precipitation below 15 mm/day.

24 Figure 18 .
Figure 18.(a-c) POD, FAR and CSI in four study regions for precipitation above 15 mm/day.

Figure 18 .
Figure 18.(a-c) POD, FAR and CSI in four study regions for precipitation above 15 mm/day.

Table 2 .
Number of stations and pixels in each area.

Table 2 .
Number of stations and pixels in each area.

Table 3 .
Contingency table to evaluate precipitation occurrence by satellite products.

Table 5 .
Statistical indices of daily precipitation events in Bushehr area (mid-March 2014 to February 2015).

Table 6 .
Statistical indices of daily precipitation events in Kermanshah area (mid-March 2014 to February 2015).

Table 7 .
Statistical indices of daily precipitation events in Tehran area (mid-March 2014 to February 2015).

Table 4 .
Statistical indices of daily precipitation events in Guilan area (mid-March 2014 to February 2015).

Table 5 .
Statistical indices of daily precipitation events in Bushehr area (mid-March 2014 to February 2015).Interim underestimate the observed values.As seen in Table6, IMERG has a better value for RMSE, Bias, Mbias, Rbias and CC in comparison with other data sets while IMERG performs better when precipitation is greater than 10 mm.
Figure6presents the scatter plot of daily precipitation over Kermanshah area (G5) where IMERG, 3B42 and ERA-

Table 6 .
Statistical indices of daily precipitation events in Kermanshah area (mid-March 2014 to February 2015).

Table 7 .
Statistical indices of daily precipitation events in Tehran area (mid-March 2014 to February 2015).

Table 8 .
Statistical indices of monthly precipitation in Guilan area (mid-March 2014 to February 2015).

Table 9 .
Statistical indices of monthly precipitation in Bushehr area (mid-March 2014 to February 2015).

Table 8 .
Statistical indices of monthly precipitation in Guilan area (mid-March 2014 to February 2015).

Table 9 .
Statistical indices of monthly precipitation in Bushehr area (mid-March 2014 to February 2015).

Table 10 .
Statistical indices of monthly precipitation in Kermanshah area (mid-March 2014 to February 2015).

Table 10 .
Statistical indices of monthly precipitation in Kermanshah area (mid-March 2014 to February 2015).

Table 10 .
Statistical indices of monthly precipitation in Kermanshah area (mid-March 2014 to February 2015).

Table 11 .
Statistical indices of monthly precipitation in Tehran area (mid-March 2014 to February 2015).

Table 12 .
Statistical indices of seasonal precipitation in Guilan area.underestimates in the Kermanshah region.The IMERG yields a correlation coefficient of 0.95 which is a good agreement between satellite estimation and rain-gauge observation and there is no great difference for MAE and Bias among all data sets.In all, ERA-Interim indicates a strong bias during autumn and a slight bias during winter.

Table 12 .
Statistical indices of seasonal precipitation in Guilan area.

Table 13 .
Statistical indices of seasonal precipitation in Bushehr area.

Table 13 .
Statistical indices of seasonal precipitation in Bushehr area.

Table 14 .
Statistical indices of seasonal precipitation in Kermanshah area.

Table 14 .
Statistical indices of seasonal precipitation in Kermanshah area.

Table 14 .
Statistical indices of seasonal precipitation in Kermanshah area.

Table 15 .
Statistical indices of seasonal precipitation in Tehran area.

Table 15 .
Statistical indices of seasonal precipitation in Tehran area.

Table 16 .
Spatially averaged POD, FAR and CSI metrics.
and Table17, in the Guilan area, ERA-Interim performs better than IMERG and 3B42 for precipitation below 15 mm/day which might be due to high amount of moisture in the atmosphere in this area that the satellites observed, although precipitation did not occur.Whereas, IMERG yields better results in Tehran and Bushehr areas in comparison to the other products.Moreover, in Kermanshah, IMERG and ERA-Interim indicated rather similar behavior, while in all study areas, 3B42 indicated low values of indices for precipitation below 15 mm/day.

Table 17 .
Spatially averaged POD, FAR and CSI metrics for precipitation below 15 mm/day.

Comparison of CSI for precipitation below 15 mm/day Figure
17. (a-c) POD, FAR and CSI in four study region for precipitation below 15 mm/day.

Table 17 .
Spatially averaged POD, FAR and CSI metrics for precipitation below 15 mm/day.Evaluation of Dichotomous Estimates/Forecasts for Precipitation above 15 mm

Table 18 .
Spatially averaged POD, FAR and CSI metrics for precipitation above 15 mm/day.