1. Introduction
Streamflow is one of the most important components of the hydrological cycle. Many efforts have been made to develop different models to emulate the hydrological cycle and simulate streamflow. Examples are statistical data-driven (e.g., Kim and Barros 2001; Sahoo et al. 2006; Piotrowski et al. 2006) or physically based (Estupina-Borrell et al. 2006; Sirdas and Sen 2007; Beven 2011) models in the forms of lumped [e.g., Sacramento Soil Moisture Accounting (SAC-SMA; Burnash et al. 1973) and Hydrologiska Byråns Vattenbalansavdelning (HBV; Bergstrom 1995)], semilumped [e.g., Variable Infiltration Capacity model (VIC); Liang et al. 1994], and distributed [e.g., Hydrology Laboratory Research Distributed Hydrologic Model (HL-RDHM); Koren et al. 2003, 2004, 2007]. One of the most, if not the most, important criteria in all of these types of hydrological modeling schemes, but particularly in the distributed format, is the availability of high-quality data with desirable spatial and temporal coverages. Satellite products with their global coverage are very well suited for this purpose. With the advancement in remote sensing science and technology, high-resolution data and information about the earth’s surface characteristics (e.g., topography, soil types, and land use) and hydrometeorological forcings (e.g., precipitation, temperature, and evapotranspiration) have been made available globally. Particularly, remote sensing of precipitation—one of the key hydrometeorological variables in generating floods—has gained significant attention in the recent past. Numerous efforts have been made to produce satellite-based precipitation estimates at high spatiotemporal resolution in global scale. Examples are the CPC morphing technique (CMORPH; Joyce et al. 2004), Precipitation Estimation from Remotely Sensed Information Using Artificial Neural Networks (PERSIANN; Hsu et al. 1997, 1999; Sorooshian et al. 2000), PERSIANN–Climate Data Record (PERSIANN-CDR; Ashouri et al. 2015), Tropical Rainfall Measuring Mission (TRMM) Multisatellite Precipitation Analysis (TMPA; Huffman et al. 2007), the Naval Research Laboratory (NRL)-Blend satellite rainfall estimates (Turk et al. 2010), and the Integrated Multisatellite Retrievals for Global Precipitation Measurement (IMERG; Huffman et al. 2015). Such products are valuable sources of information and data for flood modeling in a distributed format.
Previous efforts have been made in evaluating the accuracy of different satellite-based precipitation products against gauge observations (e.g., Sorooshian et al. 2000; Hong et al. 2006; Ebert et al. 2007; Miao et al. 2015; Ashouri et al. 2016) and utilizing such data in different applications, especially hydrological modeling (see Behrangi et al. 2011; Bajracharya et al. 2015; Maggioni et al. 2013; Nguyen et al. 2014, 2015, 2016; Seyyedi et al. 2015). Feasibility of using satellite-based precipitation as input for hydrologic simulation has been demonstrated in the Mediterranean (Ciabatta et al. 2016) and for simulating high flows in Africa (Thiemig et al. 2013). However, certain challenges and limitations with using satellite-based precipitation for rainfall–runoff modeling have been identified. Of reoccurring concern is bias in satellite-based precipitation estimates that carries over to hydrological simulations when used as model input (Guetter et al. 1996; Stisen and Sandholt 2010; Thiemig et al. 2013). Harris et al. (2007) also point to spatial resolution of coarse satellite-based precipitation products as a concern for use in hydrologic models, cautioning their use in an operational setting, and recommending adjustments to the precipitation estimates even beyond simple bias adjustments.
In this study, we implement a newly developed precipitation climate data record, PERSIANN-CDR, into a long-term hydrological modeling framework to simulate historical streamflow and flood events. In other words, we seek to evaluate the performance of PERSIANN-CDR in a distributed rainfall–runoff modeling application and compare its performance with other high-resolution precipitation products. For this purpose, the National Weather Service (NWS) HL-RDHM is used. This study is conducted in two phases. In the first phase, a period of time when PERSIANN-CDR, TMPA, and stage IV gauge-adjusted radar data products are available is chosen and the hydrological modeling is performed for this period. The study period for this phase is selected from 2003 to 2010 when all the aforementioned precipitation products are available. The results of this phase reveal how PERSIANN-CDR performs compared with the well-established TMPA data product as far as generating streamflow in the study basin. In the second phase, the entire record of the 0.25° daily PERSIANN-CDR data is used to extend the modeling process back to 1983, allowing for long-term streamflow simulations conducted at a higher resolution than previously possible.
The paper is structured as follows: section 2 describes the data used in this study, section 3 provides detailed information about the methodology and modeling structure, section 4 presents the results, and section 5 concludes with a summary of the main findings.
2. Data
a. TMPA
TRMM, a joint mission between NASA and the Japan Aerospace Exploration Agency (JAXA), was launched in November 1997 with a design lifetime of 3 years. TRMM, however, produced more than 17 years of data to study tropical rainfall for weather and climate research. This mission officially came to an end on 15 April 2015. As one of the TRMM products, TMPA (Huffman et al. 2007, 2010) contains near-global (50°S–50°N) precipitation data at 3-hourly temporal resolution and 0.25° × 0.25° grid cells. In this study, version 7 of this product is used. TMPA has an established record in precipitation and hydrological modeling studies (for a complete list of TMPA citations, refer to
b. PERSIANN-CDR
PERSIANN-CDR (Ashouri et al. 2015) is a newly developed and released satellite-based precipitation product that covers more than three decades (from 1 January 1983 to present) of daily precipitation estimations at 0.25° resolution for the 60°S–60°N latitude band. PERSIANN-CDR uses the archive of infrared brightness temperature CDR from Gridded Satellite (GridSat)-B1 (Knapp 2008a,b; Knapp et al. 2011) in the International Satellite Cloud Climatology Project (ISCCP; Rossow and Schiffer 1991; Rossow and Garder 1993) as the input to the trained neural network of the PERSIANN model. The resulting rainfall estimates are then corrected for bias using the monthly Global Precipitation Climatology Project (GPCP), version 2.2, 2.5° product. The dataset has been released and made available for public access through NOAA/National Centers for Environmental Information (NCEI; http://www1.ncdc.noaa.gov/pub/data/sds/cdr/CDRs/PERSIANN/Overview.pdf). PERSIANN-CDR has already shown its usefulness for a relatively wide range of different applications (Guo et al. 2015; Solmon et al. 2015; Ceccherini et al. 2015; Yong 2015; Yang et al. 2016). Using different extreme precipitation indices, Miao et al. (2015) evaluated the performance of PERSIANN-CDR in capturing the behavior of historical extreme precipitation events over China. Their results showed the capability of PERSIANN-CDR in reproducing similar spatial and temporal patterns of daily precipitation extremes as those depicted by the East Asia (EA) ground-based gridded daily precipitation dataset. In another study by Hagos et al. (2016) that investigated changes in the frequency of landfalling atmospheric river and extreme precipitation in the simulation of the Community Earth System Model (CESM), PERSIANN-CDR was used as the “observation” precipitation. Luchetti et al. (2016) used PERSIANN-CDR in a NOAA–NASA collaborative project for updating the ENSO-based rainfall climatology for regions in Hawaii and U.S.-affiliated Pacific islands. With respect to the applicability of PERSIANN-CDR, the paper concludes that their results “solidified the ability of the high resolution PERSIANN-CDR to be more than adequate for use in long-term precipitation climatology studies.” With respect to hydrological application, Casse and Gosset (2015) used PERSIANN-CDR to study hydrological changes and flood increases in the Niger River and the city of Niamey (Niger) over the period of 1983–2013. The results showed that PERSIANN-CDR produces annual rainfall amounts comparable with those from gauge-adjusted satellite rainfall estimates and gauge data. The paper also concludes that “the PERSIANN-CDR based hydrological simulation presents a realistic inter-annual variability, and detects flooded years, but not the exact flooded period day by day” (Casse and Gosset 2015, p. 122).
c. Stage IV gauge-adjusted radar data
National Centers for Environmental Prediction (NCEP) Environmental Modeling Center (EMC) provides the stage IV gauge-adjusted precipitation product (Fulton et al. 1998) from high-resolution Doppler Next Generation Weather Radar (NEXRAD) network and hourly rain gauge data over the contiguous United States. Stage IV radar data are available at 1-, 6-, and 24-hourly scales at 4-km spatial resolution at Hydrologic Rainfall Analysis Project (HRAP) national grid system. Stage IV radar data are manually quality controlled at NWS River Forecast Centers (RFCs). (More information about stage IV data can be obtained from www.emc.ncep.noaa.gov/mmb/ylin/pcpanl/stage4/.) Stage IV precipitation data have been used in different studies (e.g., Ebert et al. 2007; Zeweldi and Gebremichael 2009; Anagnostou et al. 2010; Ashouri et al. 2015). In this study, stage IV data are first used as the reference data for evaluating satellite-based precipitation products. They are then used as a forcing data into the hydrological model to generate streamflow simulations.
In this study, all data products are first scaled to 0.25° and daily spatiotemporal resolution before use.
3. Methodology
As introduced earlier, the main goal of this study is to investigate the performance of the newly developed precipitation climate data record, PERSIANN-CDR, in a hydrological rainfall–runoff modeling application and compare its performance with other precipitation products. For this purpose, the NOAA/NWS/Office of Hydrologic Development’s HL-RDHM (Koren et al. 2003, 2004, 2007) is used as the hydrological model to simulate the streamflow using the precipitation data products. HL-RDHM has been widely used for hydrologic studies (e.g., Smith et al. 2004; Reed et al. 2007; Tang et al. 2007; Yilmaz et al. 2008; Wagener et al. 2009; Khakbaz et al. 2012; Smith et al. 2012a,b). The conceptually based SAC-SMA provides the foundation of HL-RDHM. A schematic of SAC-SMA is shown in Fig. 1. SAC-SMA features two conceptual layers, upper- and lower-zone storage, both of which have two basic components, tension water and free water. An enhanced version of SAC-SMA known as the Sacramento Soil Moisture Accounting Heat Transfer (SAC-HT) includes the use of Noah land surface model–based physics to estimate a physically meaningful soil moisture profile, allowing for the estimation of heat transfer within the column (Koren et al. 2007). The second HL-RDHM module that this study uses is a routing scheme known as Rutpix9. This scheme consists of a hillslope component in which fast (overland flow) runoff is routed over a uniform conceptual hillslope and is combined with a slow (subsurface flow) component. Following hillslope runoff generation, a channel-routing process moves water downstream using a topographically based, cell-to-cell method. In Rutpix9, the relationship between discharge and channel cross section is based on the rating curve method (NWS 2011).

Schematic of the SAC-SMA (Burnash et al. 1973; image source: http://ldas.gsfc.nasa.gov/nldas/images/SAC_schematic.jpg).
Citation: Journal of Hydrometeorology 17, 7; 10.1175/JHM-D-15-0192.1

Schematic of the SAC-SMA (Burnash et al. 1973; image source: http://ldas.gsfc.nasa.gov/nldas/images/SAC_schematic.jpg).
Citation: Journal of Hydrometeorology 17, 7; 10.1175/JHM-D-15-0192.1
Schematic of the SAC-SMA (Burnash et al. 1973; image source: http://ldas.gsfc.nasa.gov/nldas/images/SAC_schematic.jpg).
Citation: Journal of Hydrometeorology 17, 7; 10.1175/JHM-D-15-0192.1
With respect to the calibration of the hydrological model, many studies have made a case for input-specific model calibration in order to compensate for possible pitfalls associated with a specific satellite precipitation product (e.g., Ciabatta et al. 2016; Qi et al. 2016; Stisen and Sandholt 2010; Thiemig et al. 2013). While this technique allows for multiple satellite products to yield satisfactory hydrologic simulation performance despite discrepancies in precipitation estimates, it hinders efforts to improve the precipitation product as well as the model by allowing the two to compensate for one another’s deficiencies. For this reason, the work presented here relies on calibration of the model by experts who used precipitation estimates they deemed as the most dependable (based on radar in this case) for the particular model/study region. These parameters are held constant for evaluation across all precipitation products rather than performing product-specific calibration. This allows us to focus our evaluation solely on the performance of the precipitation products without their performance being altered by the potential improvements that product-specific calibration would bring. For HL-RDHM, the model has been expertly calibrated by the NWS for the Distributed Hydrologic Model Intercomparison Project, phase 2 (DMIP2), basins. A priori parameter sets have been derived from soil and land-use data for HL-RDHM. Scalar multipliers of these parameter sets are used for calibration, assuming that the relationship of the individual pixels to one another is adequately characterized in the a priori sets. In this study, we rely on these expertly calibrated parameters by NWS experts for the calibration of the hydrological model. The calibrated parameters for SAC-HT and Rutpix9 provided by NWS are summarized as basin-average values in Table 1. Detailed information about HL-RDHM can be found in its user manual (NWS 2011; Koren and Barrett 1994).
With respect to the study area, three basins from the Oklahoma test basins from DMIP2 are chosen (Fig. 2). The two smaller study basins with outlets at Osage Creek near Elm Springs, Arkansas (ELMSP), and the Illinois River at Savoy, Arkansas (SAVOY), have drainage areas of 337 and 433 km2 and feed into the third study basin, which has an outlet located on the Illinois River south of Siloam Springs, Arkansas (SLOA4), and a drainage area of 1489 km2 (Smith et al. 2012). The Illinois River basin encompasses the three study basins and is characterized by an annual rainfall of ~1200 mm (Smith et al. 2004). It is noteworthy that the selection of these three test basins is based on the availability of the a priori calibrated parameter sets from NWS, which were only available for these three basins at the time of this study.

The three study basins (SAVOY, ELMSP, and SLOA4) [modified from Smith et al. (2012)].
Citation: Journal of Hydrometeorology 17, 7; 10.1175/JHM-D-15-0192.1

The three study basins (SAVOY, ELMSP, and SLOA4) [modified from Smith et al. (2012)].
Citation: Journal of Hydrometeorology 17, 7; 10.1175/JHM-D-15-0192.1
The three study basins (SAVOY, ELMSP, and SLOA4) [modified from Smith et al. (2012)].
Citation: Journal of Hydrometeorology 17, 7; 10.1175/JHM-D-15-0192.1
Prior to setting up the hydrological modeling scheme, a preliminary evaluation on the accuracy of the utilized precipitation products against stage IV gauge-adjusted radar data is conducted. After having the model structure in place, HL-RDHM is run in two phases. In the first phase, the time period of 2003–10 when all products (PERSIANN, PERSIANN-CDR, TMPA, and stage IV) are available is selected. The resulting hydrographs from HL-RDHM forced with the three precipitation products at the three study basins are compared with the USGS streamflow observations. The results of this phase depict how PERSIANN-CDR-simulated streamflow compares with other currently available precipitation products, that is, TMPA and stage IV data. Having proven the concept, in the second phase we extend the simulation process back to 1983, where only PERSIANN-CDR data are available.







4. Results
Before conducting the hydrological modeling process, it is necessary to evaluate the accuracy of the precipitation products over the study region. A comparison of PERSIANN, PERSIANN-CDR, and TMPA precipitation products against stage IV gauge-adjusted radar data as our reference dataset is conducted. Basic precipitation evaluation statistics, including correlation coefficient, standard deviation, and root-mean-square deviation, are presented on Taylor diagrams (Fig. 3). The results show that, in the three study basins, TMPA and PERSIANN-CDR have close performances, with slightly higher correlation coefficient for TMPA (~0.8 vs 0.75 for PERSIANN-CDR) and similar RMSD (~6 mm day−1) for both products. TMPA shows a higher standard deviation (~10 mm day−1) than PERSIANN-CDR (~8 mm day−1). TMPA and PERSIANN-CDR both outperform PERSIANN, mainly because, unlike PERSIANN, TMPA and PERSIANN-CDR are gauge-adjusted precipitation products.

Precipitation comparison plots for (top) SAVOY, (middle) ELMSP, and (bottom) SLOA4 basins for 2003–10.
Citation: Journal of Hydrometeorology 17, 7; 10.1175/JHM-D-15-0192.1

Precipitation comparison plots for (top) SAVOY, (middle) ELMSP, and (bottom) SLOA4 basins for 2003–10.
Citation: Journal of Hydrometeorology 17, 7; 10.1175/JHM-D-15-0192.1
Precipitation comparison plots for (top) SAVOY, (middle) ELMSP, and (bottom) SLOA4 basins for 2003–10.
Citation: Journal of Hydrometeorology 17, 7; 10.1175/JHM-D-15-0192.1



Simulated and observed streamflow hydrographs and respective scatterplots at the outlet of SAVOY basin using (from top to bottom) stage IV, TMPA, PERSIANN, and PERSIANN-CDR precipitation products. The solid black line shows the USGS observations.
Citation: Journal of Hydrometeorology 17, 7; 10.1175/JHM-D-15-0192.1

Simulated and observed streamflow hydrographs and respective scatterplots at the outlet of SAVOY basin using (from top to bottom) stage IV, TMPA, PERSIANN, and PERSIANN-CDR precipitation products. The solid black line shows the USGS observations.
Citation: Journal of Hydrometeorology 17, 7; 10.1175/JHM-D-15-0192.1
Simulated and observed streamflow hydrographs and respective scatterplots at the outlet of SAVOY basin using (from top to bottom) stage IV, TMPA, PERSIANN, and PERSIANN-CDR precipitation products. The solid black line shows the USGS observations.
Citation: Journal of Hydrometeorology 17, 7; 10.1175/JHM-D-15-0192.1

As in Fig. 4, but for ELMSP basin.
Citation: Journal of Hydrometeorology 17, 7; 10.1175/JHM-D-15-0192.1

As in Fig. 4, but for ELMSP basin.
Citation: Journal of Hydrometeorology 17, 7; 10.1175/JHM-D-15-0192.1
As in Fig. 4, but for ELMSP basin.
Citation: Journal of Hydrometeorology 17, 7; 10.1175/JHM-D-15-0192.1

As in Fig. 4, but for SLOA4 basin.
Citation: Journal of Hydrometeorology 17, 7; 10.1175/JHM-D-15-0192.1

As in Fig. 4, but for SLOA4 basin.
Citation: Journal of Hydrometeorology 17, 7; 10.1175/JHM-D-15-0192.1
As in Fig. 4, but for SLOA4 basin.
Citation: Journal of Hydrometeorology 17, 7; 10.1175/JHM-D-15-0192.1

Taylor diagram showing the correlation coefficient, standard deviation, and root-mean-square deviation from the PERSIANN-, PERSIANN-CDR-, TMPA-, and stage IV–derived hydrographs for (top) SAVOY, (middle) ELMSP, and (bottom) SLOA4 basins for 2003–10.
Citation: Journal of Hydrometeorology 17, 7; 10.1175/JHM-D-15-0192.1

Taylor diagram showing the correlation coefficient, standard deviation, and root-mean-square deviation from the PERSIANN-, PERSIANN-CDR-, TMPA-, and stage IV–derived hydrographs for (top) SAVOY, (middle) ELMSP, and (bottom) SLOA4 basins for 2003–10.
Citation: Journal of Hydrometeorology 17, 7; 10.1175/JHM-D-15-0192.1
Taylor diagram showing the correlation coefficient, standard deviation, and root-mean-square deviation from the PERSIANN-, PERSIANN-CDR-, TMPA-, and stage IV–derived hydrographs for (top) SAVOY, (middle) ELMSP, and (bottom) SLOA4 basins for 2003–10.
Citation: Journal of Hydrometeorology 17, 7; 10.1175/JHM-D-15-0192.1
Mean, standard deviation (std dev), centered RMSE, CORR, ln E, BIAS, and d for simulated streamflow (2003–10) from PERSIANN, PERSIANN-CDR, stage IV, and TMPA compared to USGS streamflow observations.


In addition to the quantitative analysis over the entire record, an analysis on the long-term (2003–10) annual cycle from daily USGS streamflow observation, as well as the stage IV–, TMPA-, PERSIANN-, and PERSIANN-CDR-derived hydrographs for the three study basins, is conducted. As shown in Fig. 8, PERSIANN-CDR and TMPA depict a relatively similar performance for all seasons. PERSIANN’s performance degrades in the spring season. The improvements from PERSIANN to PERSIANN-CDR are also evident in Fig. 8. Respective daily based statistical measures from day-of-year long-term (2003–10) annual cycle analyses are calculated against gauge observations and shown in Table 3. The reduction in the percentage volume bias from PERSIANN to PERSIANN-CDR is evident. At the SAVOY basin, TMPA shows the highest index of agreement, where this measure is the highest for PERSIANN-CDR. Stage IV shows the lowest RMSE in all the three catchments. PERSIANN-CDR shows a higher RMSE than TMPA at SAVOY but a lower RMSE at ELMSP basin. RMSEs from PERSIANN-CDR and TMPA are very close at the SLOA4 basin. With respect to the CORR, stage IV radar data show the best performance at SAVOY and ELMSP basins. PERSIANN-CDR shows a slightly higher CORR than TMPA at the ELMSP basin, but TMPA outperforms PERSIANN-CDR in the SAVOY and SLOA4 basins in this regard. The large negative bias in stage IV radar data in ELMSP and SLOA4 is again observed in this analysis. The seasonal analysis (figures not presented) shows that the largest negative bias in stage IV–derived streamflow happens in the summer and fall seasons.

Long-term (2003–10) annual cycle from USGS observation and stage IV–, TMPA-, PERSIANN-, and PERSIANN-CDR-derived hydrographs for (top) SAVOY, (middle) ELMSP, and (bottom) SLOA4 basins.
Citation: Journal of Hydrometeorology 17, 7; 10.1175/JHM-D-15-0192.1

Long-term (2003–10) annual cycle from USGS observation and stage IV–, TMPA-, PERSIANN-, and PERSIANN-CDR-derived hydrographs for (top) SAVOY, (middle) ELMSP, and (bottom) SLOA4 basins.
Citation: Journal of Hydrometeorology 17, 7; 10.1175/JHM-D-15-0192.1
Long-term (2003–10) annual cycle from USGS observation and stage IV–, TMPA-, PERSIANN-, and PERSIANN-CDR-derived hydrographs for (top) SAVOY, (middle) ELMSP, and (bottom) SLOA4 basins.
Citation: Journal of Hydrometeorology 17, 7; 10.1175/JHM-D-15-0192.1
Results of phase 1 analyses for the period when other high-resolution, satellite-derived precipitation products are available show that the performance of PERSIANN-CDR has been close to other precipitation products. The findings of this phase reveal the high potential for PERSIANN-CDR data use in rainfall–runoff modeling applications, particularly long-term simulations of more than three decades, given the fact that PERSIANN-CDR rainfall data spans from 1 January 1983 to present (delayed) time. Therefore, in phase 2, HL-RDHM is forced with daily PERSIANN-CDR rainfall estimation for 1983–2012 to reconstruct historical streamflow at the three study basins. It is important to note that, similar to phase 1, we rely on the NWS a priori parameter sets for the calibration of the model. The resulting hydrographs for SAVOY, ELMSP, and SLOA4 are shown in Fig. 9, with the black line being the USGS observation and the blue line being the PERSIANN-CDR-derived hydrograph. As shown, an immediate result is that PERSIANN-CDR could reconstruct the full record of historical streamflow, especially for years prior to 1996 (except 1986 for SAVOY basin) when generally the USGS streamflow observations are not available for the three study basins. The scatterplots comparing the PERSIANN-CDR-simulated streamflow against USGS observations for 1983–2012 for the three study basins are shown in Fig. 9. The summary of quantitative comparisons is presented in Table 4. The results depict high CORR (~0.67–0.73), relatively low BIAS (~5%–12%), and high index of agreement criterion (~0.68–0.83) between PERSIANN-CDR-simulated daily streamflow and USGS daily observations, demonstrating reasonable performance of PERSIANN-CDR despite the calibration of the hydrological model by a precipitation product other than PERSIANN-CDR. It is also noteworthy that the efficiency index of the logarithmic Nash–Sutcliffe coefficient derived in the second phase of the study is comparable, and even slightly better, than the same coefficient in the first phase of the study for other precipitation products.

Long-term (1983–2012) simulated streamflow from PERSIANN-CDR daily precipitation data (blue) vs USGS streamflow observations (black) for (top) SAVOY, (middle) ELMSP, and (bottom) SLOA4 basins, plus the respective scatterplots (right).
Citation: Journal of Hydrometeorology 17, 7; 10.1175/JHM-D-15-0192.1

Long-term (1983–2012) simulated streamflow from PERSIANN-CDR daily precipitation data (blue) vs USGS streamflow observations (black) for (top) SAVOY, (middle) ELMSP, and (bottom) SLOA4 basins, plus the respective scatterplots (right).
Citation: Journal of Hydrometeorology 17, 7; 10.1175/JHM-D-15-0192.1
Long-term (1983–2012) simulated streamflow from PERSIANN-CDR daily precipitation data (blue) vs USGS streamflow observations (black) for (top) SAVOY, (middle) ELMSP, and (bottom) SLOA4 basins, plus the respective scatterplots (right).
Citation: Journal of Hydrometeorology 17, 7; 10.1175/JHM-D-15-0192.1
BIAS, CORR, and RMSE statistics for simulated streamflow from PERSIANN-CDR against USGS observed streamflow for 1983–2012.


5. Discussion
There are few points what we would like to clarify via the explanations below. First, we would like to point out that the main reason for including the real-time PERSIANN product in our analysis was to investigate the progress that the Center for Hydrometeorology and Remote Sensing at the University of California, Irvine, has made over time in improving its precipitation products. Regarding this, the improvements are evident from both the precipitation and simulated streamflow points of view. That said, there are certain facts that should be considered here since a blind comparison of real-time precipitation products [e.g., PERSIANN or TMPA–real time (3B42RT)] with gauge-adjusted products (e.g., TMPA and PERSIANN-CDR) may not be a fair comparison in general. In the case of PERSIANN and PERSIANN-CDR, the latter is indirectly gauge corrected based on the gauge information from the Global Precipitation Climatology Centre (GPCC; Schneider et al. 2008) used in the GPCP monthly precipitation product (see Huffman and Bolvin 2013), whereas PERSIANN is a real-time product without any gauge-adjusting component. Real-time PERSIANN data are useful for real-time applications such as global flood monitoring, while the long-term (+30 years), daily PERSIANN-CDR rainfall data make it useful for long-term hydrological and climatological studies and applications. As is the case with other satellite-based products, the choice of PERSIANN or PERSIANN-CDR depends on the type of application, period of the study, and the constraints of the problem.
The second point refers to the relatively high percent bias observed in the stage IV streamflow simulations for ELMSP and SLOA4 basins compared to PERSIANN-CDR and TMPA. There could be a couple of reasons for that. First, the validation period in this study is not the same as the calibration period for which the a priori parameters are estimated, so there is no guarantee that stage IV is always better than the other products in the discharge simulations. More importantly, calibration performed by the NWS (see Kuzmin et al. 2008) is based on a multi-time-scale form of RMSE and does not include bias in the objective function. In Table 2, stage IV simulation has the lowest RMSE that is consistent with the calibration methods. It is noteworthy that the correlations of stage IV simulations in all cases are better than the others.
The third point is about the low Nash–Sutcliff efficiency values obtained from the simulated streamflow time series. The key point here is the calibration of the hydrological model. As explained in the introduction and methodology sections, instead of conducting product-specific calibration, we relied on the NWS expertly calibrated a priori parameter sets for the calibration of the HL-RDHM for the study basins. These parameters are kept fixed for all the precipitation products, allowing us to focus our evaluation solely on the performance of the precipitation products, rather than mixing it with the improvements that product-specific calibrations of the hydrological model can, and do, introduce in the final simulations. It is obvious that by conducting product-specific calibration, one would achieve higher-efficiency results.
The fourth point is about the performance of PERSIANN-CDR in reproducing the observed peak flows during the 1983–2012 period (Fig. 9), when observations are available. Part of the mismatch is due to the fact that our hydrological model was not specifically calibrated with PERSIANN-CDR data. In addition to that, it could be partly due to the fault of the forcing precipitation, the fault of the model, or even the fault of the precipitation product used for calibration (or a combination of all).
6. Conclusions
The main goal of this study was to evaluate the performance of the newly developed precipitation climate data record, PERSIANN-CDR, in a rainfall–runoff modeling scheme and compare its performance with other high-resolution precipitation products. In examining the accuracy of the precipitation products, PERSIANN-CDR and TMPA showed close performances compared to the stage IV gauge-adjusted radar data product. Focusing only on the PERSIANN products, it is found that PERSIANN-CDR outperforms the PERSIANN real-time product, depicting better statistical measures. This is mainly because PERSIANN-CDR is gauge corrected, whereas PERSIANN is a real-time product with no gauge-correction component.
For the purpose of evaluating PERSIANN-CDR’s application in hydrological modeling, two phases of study were designed. In the first phase, the main goal was to test how PERSIANN-CDR’s performance compares with the performance of other precipitation products. To have all the products available, the time period of 2003–10 was selected. The HL-RDHM was run separately for each of these precipitation data products to simulate streamflow hydrographs at the outlets of three DMIP2 study basins where NWS a priori parameters are available. The simulations at SAVOY, ELMSP, and SLOA4 basins were compared with USGS observations. The results show that PERSIANN-CDR- and TMPA-derived simulations have close performances with higher correlation coefficients and better RMSEs for TMPA and lower biases for PERSIANN-CDR. Annual cycle analysis of simulated hydrographs also depicts close performance between TMPA and PERSIANN-CDR. Phase 1 of this study serves as the proof of concept regarding the applicability of PERSIANN-CDR in rainfall–runoff modeling. Given this result and the fact that PERSIANN-CDR precipitation data span from 1983 to the present, we could extend the simulation process back to 1983 to reconstruct the historical record of streamflow. This is particularly important when even USGS streamflow observations are not available prior to the year 1996 for the three study basins. In this phase of the study, only PERSIANN-CDR precipitation data were available as the forcing to the model. The resulting PERSIANN-CDR-derived hydrographs were compared with USGS observations depicting high CORR, relatively low BIAS, and high index of agreement criterion.
It is noteworthy that for the three DMIP2 study basins there were periods of time, mostly before 1996, when USGS daily data were not available. Using PERSIANN-CDR and HL-RDHM, we could simulate the streamflow for those periods and fill the gaps. This is particularly important for long-term trend studies where full data coverage over a long period of time, at least 30 years according to a World Meteorological Organization (WMO) report (Burroughs 2003), is needed. To conclude, PERSIANN-CDR could prove its usefulness for long-term hydrological rainfall–runoff modeling and streamflow simulation. It can be particularly helpful for simulating streamflow in ungauged basins.
Acknowledgments
Partial financial support for this research was provided from the National Oceanic and Atmospheric Administration’s (NOAA) Cooperative Institute for Climate and Satellites (CICS) and NOAA’s National Centers for Environmental Information (NCEI) Climate Data Record program (Prime Award NA14NES4320003 and NCSU CICS Subaward 2014-2918-03), and fellowship programs from the NASA Earth and Space Science Fellowship (NESSF, NNX12AO11H) and Department of Defense (DoD) National Defense Science and Engineering Graduate Fellowship (NDSEG).
REFERENCES
Anagnostou, E. N., Maggioni V. , Nikolopoulos E. I. , Meskele T. , Hossain F. , and Papadopoulos A. , 2010: Benchmarking high-resolution global satellite rainfall products to radar and rain-gauge rainfall estimates. IEEE Trans. Geosci. Remote Sens., 48, 1667–1683, doi:10.1109/TGRS.2009.2034736.
Ashouri, H., Hsu K.-L. , Sorooshian S. , Braithwaite D. K. , Knapp K. R. , Cecil L. D. , Nelson B. R. , and Prat O. P. , 2015: PERSIANN-CDR: Daily precipitation climate data record from multisatellite observations for hydrological and climate studies. Bull. Amer. Meteor. Soc., 96, 69–83, doi:10.1175/BAMS-D-13-00068.1.
Ashouri, H., Sorooshian S. , Hsu K. , Bosilovich M. G. , Lee J. , and Wehner M. F. , 2016: Evaluation of NASA’s MERRA precipitation product in reproducing the observed trend and distribution of extreme precipitation events in the United States. J. Hydrometeor., 17, 693–711, doi:10.1175/JHM-D-15-0097.1.
Bajracharya, S. R., Palash W. , Shrestha M. S. , Khadgi V. R. , Duo C. , Das P. J. , and Dorji C. , 2015: Systematic evaluation of satellite-based rainfall products over the Brahmaputra basin for hydrological applications. Adv. Meteor., 2015, 398687, doi:10.1155/2015/398687.
Behrangi, A., Khakbaz B. , Jaw T. C. , AghaKouchak A. , Hsu K. , and Sorooshian S. , 2011: Hydrologic evaluation of satellite precipitation products over a mid-size basin. J. Hydrol., 397, 225–237, doi:10.1016/j.jhydrol.2010.11.043.
Bergstrom, S., 1995: The HBV model. Computer Models of Watershed Hydrology, V. P. Singh, Ed., Water Resources Publications, 443–476.
Beven, K. J., 2011: Rainfall–Runoff Modelling: The Primer. Wiley, 488 pp.
Burnash, R. J. C., Ferral R. L. , and McGuire R. A. , 1973: A generalized streamflow simulation system: Conceptual modeling for digital computers. NWS/California Department of Water Resources Tech. Rep., 204 pp.
Burroughs, W., 2003: Climate into the 21st Century. Cambridge University Press, 240 pp.
Casse, C., and Gosset M. , 2015: Analysis of hydrological changes and flood increase in Niamey based on the PERSIANN-CDR satellite rainfall estimate and hydrological simulations over the 1983–2013 period. Proc. IAHS, 370, 117–123, doi:10.5194/piahs-370-117-2015.
Ceccherini, G., Ameztoy I. , Hernández C. P. R. , and Moreno C. C. , 2015: High-resolution precipitation datasets in South America and West Africa based on satellite-derived rainfall, enhanced vegetation index and digital elevation model. Remote Sens., 7, 6454–6488, doi:10.3390/rs70506454.
Ciabatta, L., Brocca L. , Massari C. , Moramarco T. , Gabellani S. , Puca S. , and Wagner W. , 2016: Rainfall–runoff modelling by using SM2RAIN-derived and state-of-the-art satellite rainfall products over Italy. Int. J. Appl. Earth Obs. Geoinf., 48, 163–173, doi:10.1016/j.jag.2015.10.004.
Collischonn, B., Collischonn W. , and Morelli Tucci C. E. , 2008: Daily hydrological modeling in the Amazon basin using TRMM rainfall estimators. J. Hydrol., 360, 207–216, doi:10.1016/j.jhydrol.2008.07.032.
Ebert, E. E., Janowiak J. E. , and Kidd C. , 2007: Comparison of near-real-time precipitation estimates from satellite observations and numerical models. Bull. Amer. Meteor. Soc., 88, 47–64, doi:10.1175/BAMS-88-1-47.
Estupina-Borrell, V., Dartus D. , and Ababou R. , 2006: Flash flood modeling with the MARINE hydrological distributed model. Hydrol. Earth System Sci. Discuss., 3, 3397–3438, doi:10.5194/hessd-3-3397-2006.
Fulton, R. A., Breidenbach J. P. , Seo D. J. , Miller D. A. , and O’Bannon T. , 1998: The WSR-88D rainfall algorithm. Wea. Forecasting, 13, 377–395, doi:10.1175/1520-0434(1998)013<0377:TWRA>2.0.CO;2.
Guetter, A. K., Georgakakos K. P. , and Tsonis A. A. , 1996: Hydrologic applications of satellite data: 2. Flow simulation and soil water estimates. J. Geophys. Res., 101, 26 527–26 538, doi:10.1029/96JD01655.
Guo, H., Chen S. , Bao A. , Hu J. , Gebregiorgis A. S. , Xue X. , and Zhang X. , 2015: Inter-comparison of high-resolution satellite precipitation products over central Asia. Remote Sens., 7, 7181–7211, doi:10.3390/rs70607181.
Hagos, S. M., Leung L. R. , Yoon J. H. , Lu J. , and Gao Y. , 2016: A projection of changes in landfalling atmospheric river frequency and extreme precipitation over western North America from the large ensemble CESM simulations. Geophys. Res. Lett., 43, 1357–1363, doi:10.1002/2015GL067392.
Harris, A., Rahman S. , Hossain F. , Yarborough L. , Bagtzoglou A. C. , and Easson G. , 2007: Satellite-based flood modeling using TRMM-based rainfall products. Sensors, 7, 3416–3427, doi:10.3390/s7123416.
Hogue, T. S., Sorooshian S. , Gupta H. , Holz A. , and Braatz D. , 2000: A multistep automatic calibration scheme for river forecasting models. J. Hydrometeor., 1, 524–542, doi:10.1175/1525-7541(2000)001<0524:AMACSF>2.0.CO;2.
Hong, Y., Adler R. , and Huffman G. , 2006: Evaluation of the potential of NASA multi-satellite precipitation analysis in global landslide hazard assessment. Geophys. Res. Lett., 33, L22402, doi:10.1029/2006GL028010.
Hsu, K., Gao X. , Sorooshian S. , and Gupta H. V. , 1997: Precipitation Estimation From Remotely Sensed Information Using Artificial Neural Networks. J. Appl. Meteor. Climatol., 36, 1176–1190, doi:10.1175/1520-0450(1997)036<1176:PEFRSI>2.0.CO;2.
Hsu, K., Gupta H. V. , Gao X. , and Sorooshian S. , 1999: Estimation of physical variables from multichannel remotely sensed imagery using a neural network: Application to rainfall estimation. Water Resour. Res., 35, 1605–1618, doi:10.1029/1999WR900032.
Huffman, G. J., and Bolvin D. T. , 2013: GPCP version 2.2 SG combined precipitation data set documentation. NASA GSFC Doc., 46 pp. [Available online at ftp://precip.gsfc.nasa.gov/pub/gpcp-v2.2/doc/V2.2_doc.pdf.]
Huffman, G. J., and Coauthors, 2007: The TRMM Multisatellite Precipitation Analysis (TMPA): Quasi-global, multiyear, combined-sensor precipitation estimates at fine scales. J. Hydrometeor., 8, 38–55, doi:10.1175/JHM560.1.
Huffman, G. J., Adler R. F. , Bolvin D. T. , and Nelkin E. J. , 2010: The TRMM Multi-Satellite Precipitation Analysis (TMPA). Satellite Rainfall Applications for Surface Hydrology, M. Gebremichael and F. Hossain, Eds., Springer, 3–22, doi:10.1007/978-90-481-2915-7_1.
Huffman, G. J., Bolvin D. T. , Braithwaite D. , Hsu K. , Joyce R. , Kidd C. , Nelkin E. J. , and Xie P. , 2015: NASA Global Precipitation Measurement Integrated Multi-satellitE Retrievals for GPM (IMERG). Algorithm Theoretical Basis Doc., version 4.5, 30 pp. [Available online at http://pmm.nasa.gov/sites/default/files/document_files/IMERG_ATBD_V4.5.pdf.]
Joyce, R. J., Janowiak J. E. , Arkin P. A. , and Xie P. , 2004: CMORPH: A method that produces global precipitation estimates from passive microwave and infrared data at high spatial and temporal resolution. J. Hydrometeor., 5, 487–503, doi:10.1175/1525-7541(2004)005<0487:CAMTPG>2.0.CO;2.
Khakbaz, B., Imam B. , Hsu K. , and Sorooshian S. , 2012: From lumped to distributed via semi-distributed: Calibration strategies for semi-distributed hydrologic models. J. Hydrol., 418–419, 61–77, doi:10.1016/j.jhydrol.2009.02.021.
Kim, G., and Barros A. P. , 2001: Quantitative flood forecasting using multisensor data and neural networks. J. Hydrol., 246, 45–62, doi:10.1016/S0022-1694(01)00353-5.
Knapp, K. R., 2008a: Scientific data stewardship of International Satellite Cloud Climatology Project B1 global geostationary observations. J. Appl. Remote Sens., 2, 023548, doi:10.1117/1.3043461.
Knapp, K. R., 2008b: Calibration assessment of ISCCP geostationary infrared observations using HIRS. J. Atmos. Oceanic Technol., 25, 183–195, doi:10.1175/2007JTECHA910.1.
Knapp, K. R., and Coauthors, 2011: Globally gridded satellite observations for climate studies. Bull. Amer. Meteor. Soc., 92, 893–907, doi:10.1175/2011BAMS3039.1.
Koren, V., and Barrett C. B. , 1994: Satellite based distributed monitoring, forecasting and simulation (MFS) system for the Nile River. Proc. Int. Workshop on Application of Remote Sensing in Hydrology, NHRI Symp. 14, Saskatoon, Canada, NHRI, 187–200.
Koren, V., Smith M. , and Duan Q. , 2003: Use of a priori parameter estimates in the derivation of spatially consistent parameter sets of rainfall–runoff models. Calibration of Watershed Models, Q. Duan et al., Eds., Water Science and Application Series, Vol. 6, Amer. Geophys. Union, 239–254.
Koren, V., Reed S. , Smith M. , Zhang Z. , and Seo D.-J. , 2004: Hydrology laboratory research modeling system (HL-RMS) of the US National Weather Service. J. Hydrol., 291, 297–318, doi:10.1016/j.jhydrol.2003.12.039.
Koren, V., Smith M. , Cui Z. , and Cosgrove B. , 2007: Physically-based modifications to the Sacramento Soil Moisture Accounting model: Modeling the effects of frozen ground on the rainfall–runoff process. NOAA Tech. Rep. NWS 52, 43 pp. [Available online at www.nws.noaa.gov/oh/hrl/hsmb/docs/hydrology/PBE_SAC-SMA/NOAA_Technical_Report_NWS_52.pdf.]
Kuzmin, V., Seo D.-J. , and Koren V. , 2008: Fast and efficient optimization of hydrologic model parameters using a priori estimates and stepwise line search. J. Hydrol., 353, 109–128, doi:10.1016/j.jhydrol.2008.02.001.
Legates, D. R., and McCabe G. J. Jr., 1999: Evaluating the use of “goodness-of-fit” measures in hydrologic and hydroclimatic model validation. Water Resour. Res., 35, 233–241, doi:10.1029/1998WR900018.
Liang, X., Lettenmaier D. P. , Wood E. F. , and Burges S. J. , 1994: A simple hydrologically based model of land surface water and energy fluxes for GSMs. J. Geophys. Res., 99, 14 415–14 428, doi:10.1029/94JD00483.
Luchetti, N. T., Sutton J. R. , Wright E. E. , Kruk M. C. , and Marra J. J. , 2016: When El Niño rages: How satellite data can help water-stressed islands. Bull. Amer. Meteor. Soc., doi:10.1175/BAMS-D-15-00219.1, in press.
Maggioni, V., Vergara H. J. , Anagnostou E. N. , Gourley J. J. , Hong Y. , and Stampoulis D. , 2013: Investigating the applicability of error correction ensembles of satellite rainfall products in river flow simulations. J. Hydrometeor., 14, 1194–1211, doi:10.1175/JHM-D-12-074.1.
Miao, C., Ashouri H. , Hsu K.-L. , Sorooshian S. , and Duan Q. , 2015: Evaluation of the PERSIANN-CDR daily rainfall estimates in capturing the behavior of extreme precipitation events over China. J. Hydrometeor., 16, 1387–1396, doi:10.1175/JHM-D-14-0174.1.
Nash, J. E., and Sutcliffe J. V. , 1970: River flow forecasting through conceptual models. Part I—A discussion of principles. J. Hydrol., 10, 282–290, doi:10.1016/0022-1694(70)90255-6.
Nguyen, P., Sellars S. , Thorstensen A. , Tao Y. , Ashouri H. , Braithwaite D. , Hsu K. , and Sorooshian S. , 2014: Satellites track precipitation of Super Typhoon Haiyan. Eos, Trans. Amer. Geophys. Union, 95, 133–135, doi:10.1002/2014EO160002.
Nguyen, P., Thorstensen A. , Sorooshian S. , Hsu K. , and AghaKouchak A. , 2015: Flood forecasting and inundation mapping using HiResFlood-UCI and near-real-time satellite precipitation data: The 2008 Iowa flood. J. Hydrometeor., 16, 1171–1183, doi:10.1175/JHM-D-14-0212.1.
Nguyen, P., and Coauthors, 2016: A high resolution coupled hydrologic–hydraulic model (HiResFlood-UCI) for flash flood modeling. J. Hydrol., doi:10.1016/j.jhydrol.2015.10.047, in press.
NWS, 2011: Hydrology Laboratory-Research Distributed Hydrologic Model (HL-RDHM) user manual V. 3.2.0. National Weather Service Rep., 79 pp.
Piotrowski, A., Napiorkowski J. J. , and Rowinski P. M. , 2006: Flash-flood forecasting by means of neural networks and nearest neighbour approach—A comparative study. Nonlinear Processes Geophys., 13, 443–448, doi:10.5194/npg-13-443-2006.
Qi, W., Zhang C. , Fu G. , Sweetapple C. , and Zhou H. , 2016: Evaluation of global fine-resolution precipitation products and their uncertainty quantification in ensemble discharge simulations. Hydrol. Earth Syst. Sci., 20, 903–920, doi:10.5194/hess-20-903-2016.
Reed, S., Schaake J. , and Zhang Z. , 2007: A distributed hydrologic model and threshold frequency-based method for flash flood forecasting at ungauged locations. J. Hydrol., 337, 402–420, doi:10.1016/j.jhydrol.2007.02.015.
Rossow, W. B., and Schiffer R. A. , 1991: ISCCP cloud data products. Bull. Amer. Meteor. Soc., 72, 2–20, doi:10.1175/1520-0477(1991)072<0002:ICDP>2.0.CO;2.
Rossow, W. B., and Garder L. C. , 1993: Cloud detection using satellite measurements of infrared and visible radiances for ISCCP. J. Climate, 6, 2341–2369, doi:10.1175/1520-0442(1993)006<2341:CDUSMO>2.0.CO;2.
Sahoo, G. B., Ray C. , and De Carlo E. H. , 2006: Calibration and validation of a physically distributed hydrological model, MIKE SHE, to predict streamflow at high frequency in a flashy mountainous Hawaii stream. J. Hydrol., 327, 94–109, doi:10.1016/j.jhydrol.2005.11.012.
Schneider, U., Fuchs T. , Meyer-Christoffer A. , and Rudolf B. , 2008: Global precipitation analysis products of the GPCC. Global Precipitation Climatology Centre Doc., DWD, 14 pp. [Available online at ftp://ftp.dwd.de/pub/data/gpcc/PDF/GPCC_intro_products_2008.pdf.]
Seyyedi, H., Anagnostou E. N. , Beighley E. , and McCollum J. , 2015: Hydrologic evaluation of satellite and reanalysis precipitation datasets over a mid-latitude basin. Atmos. Res., 164–165, 37–48, doi:10.1016/j.atmosres.2015.03.019.
Sirdas, S., and Sen Z. , 2007: Determination of flash floods in western Arabian Peninsula. J. Hydrol. Eng., 12, 676–681, doi:10.1061/(ASCE)1084-0699(2007)12:6(676).
Smith, M. B., Seo D.-J. , Koren V. I. , Reed S. M. , Zhang Z. , Duan Q. , Moreda F. , and Cong S. , 2004: The distributed model intercomparison project (DMIP): Motivation and experiment design. J. Hydrol., 298, 4–26, doi:10.1016/j.jhydrol.2004.03.040.
Smith, M. B., and Coauthors, 2012a: The distributed model intercomparison project—Phase 2: Motivation and design of the Oklahoma experiments. J. Hydrol., 418–419, 3–16, doi:10.1016/j.jhydrol.2011.08.055.
Smith, M. B., and Coauthors, 2012b: Results of the DMIP 2 Oklahoma experiments. J. Hydrol., 418–419, 17–48, doi:10.1016/j.jhydrol.2011.08.056.
Solmon, F., Nair V. S. , and Mallet M. , 2015: Increasing Arabian dust activity and the Indian summer monsoon. Atmos. Chem. Phys., 15, 8051–8064, doi:10.5194/acp-15-8051-2015.
Sorooshian, S., Hsu K.-L. , Gao X. , Gupta H. V. , Imam B. , and Braithwaite D. , 2000: Evaluation of PERSIANN system satellite-based estimates of tropical rainfall. Bull. Amer. Meteor. Soc., 81, 2035–2046, doi:10.1175/1520-0477(2000)081<2035:EOPSSE>2.3.CO;2.
Stisen, S., and Sandholt I. , 2010: Evaluation of remote-sensing-based rainfall products through predictive capability in hydrological runoff modeling. Hydrol. Processes, 24, 879–891, doi:10.1002/hyp.7529.
Su, F. G., Hong Y. , and Lettenmaier D. P. , 2008: Evaluation of TRMM Multisatellite Precipitation Analysis (TMPA) and its utility in hydrologic prediction in La Plata basin. J. Hydrometeor., 9, 622–640, doi:10.1175/2007JHM944.1.
Tang, Y., Reed P. , Van Werkhoven K. , and Wagener T. , 2007: Advancing the identification and evaluation of distributed rainfall–runoff models using global sensitivity analysis. Water Resour. Res., 43, W06415, doi:10.1029/2006WR005813.
Thiemig, V., Rojas R. , Zambrano-Bigiarini M. , and Roo A. D. , 2013: Hydrological evaluation of satellite-based rainfall estimates over the Volta and Baro–Akobo basin. J. Hydrol., 499, 324–338, doi:10.1016/j.jhydrol.2013.07.012.
Turk, J. T., Mostovoy G. V. , and Anantharaj V. , 2010: The NRL-Blend high resolution precipitation product and its application to land surface hydrology. Satellite Rainfall Applications for Surface Hydrology, M. Gebremichael and F. Hossain, Eds., Springer, 85–104, doi:10.1007/978-90-481-2915-7_6.
Wagener, T., van Werkhoven K. , Reed P. , and Tang Y. , 2009: Multiobjective sensitivity analysis to understand the information content in streamflow observations for distributed watershed modeling. Water Resour. Res., 45, W02501, doi:10.1029/2008WR007347.
Willmot, C. J., 1984: On the evaluation of model performance in physical geography. Spatial Statistics and Models, G. L. Gaile and C. J. Willmot, Eds., D. Reidel, 443–460.
Yang, X., Yong B. , Yong H. , Chen S. , and Zhang X. , 2016: Error analysis of multi-satellite precipitation estimates with an independent raingauge observation network over a medium-sized humid basin. Hydrol. Sci. J., 61, 1813–1830, doi:10.1080/02626667.2015.1040020.
Yilmaz, K. K., Hogue T. S. , Hsu K. L. , Sorooshian S. , Gupta H. V. , and Wagener T. , 2005: Intercomparison of rain gauge, radar, and satellite-based precipitation estimates with emphasis on hydrologic forecasting. J. Hydrometeor., 6, 497–517, doi:10.1175/JHM431.1.
Yilmaz, K. K., Gupta H. V. , and Wagener T. , 2008: A process-based diagnostic approach to model evaluation: Application to the NWS distributed hydrologic model. Water Resour. Res., 44, W09417, doi:10.1029/2007WR006716.
Yong, B., 2015: Comments on “Error analysis of satellite precipitation products in mountainous basins.” J. Hydrometeor., 16, 1443–1444, doi:10.1175/JHM-D-14-0202.1.
Zeweldi, D. A., and Gebremichael M. , 2009: Evaluation of CMORPH precipitation products at fine space–time scales. J. Hydrometeor., 10, 300–307, doi:10.1175/2008JHM1041.1.