Evaluation and Hydrological Application of TRMM and GPM Precipitation Products in a Tropical Monsoon Basin of Thailand

Watershed runoff is essential for water management. However, runoff materials are lacking in poorly gauged catchments and not always accessible. Microwave remote sensing offers emerging capabilities for hydrological simulation. In this study based on multi-satellite retrievals for Global Precipitation Measurement (IMERG), Tropical Rainfall Measuring Mission (TRMM) products, and World Meteorological Organization (WMO) interpolated precipitation data, we simulated runoff using a variable infiltration capacity (VIC) model and studied the differences among the results. Then, we analyzed the impacts of the runoff on a moderate-resolution imaging spectroradiometer vegetation leaf area index (LAI) during dry seasons. The results showed that (1) IMERG V5 and TRMM products are capable of monitoring the night-day rainfall diurnal cycle and have higher correlations than the WMO daily observation interpolations. However, the WMO shows less overestimation of total precipitation than remote-sensing precipitation; (2) in the downstream, the TRMM shows better runoff simulation accuracy in the tributaries, and the WMO shows better results in the mainstreams. Therefore, at basin outlets in mainstreams, the Nash–Sutcliffe efficiency coefficients of monthly runoff by the WMO are higher than the simulations by the TRMM; (3) for the whole basin during dry seasons, the LAI variation is correlated with the outlet runoff, which is similar to the correlation with threeto six-month accumulated precipitation. TRMM products can be used to depict both precipitation deficit and runoff deficit, which cause vegetation variations. Our research suggests the potential of microwave precipitation products for detailed watershed runoff simulations and water management.


Introduction
In most of the tropics and subtropics outside the equatorial rain belt, vegetation is driven primarily by precipitation [1].Tropical monsoon climate regions feature significant dry/rainy season alteration where rice is widely planted.Water is scarce during dry seasons, so dams and canal systems are constructed to support rice planting and provide water resources.
The vegetation ecosystem has a certain tolerance to drought disturbance, and a persistent dry spell would have a significant impact on the composition, structure, and function of a vegetation ecosystem [2].Several vegetation indices by remote sensing technologies were developed to depict vegetation conditions, such as the normalized difference vegetation index (NDVI), a measure of the greenness or vigor of vegetation.The higher its value, the larger the vegetation density is; a low NDVI indicates stressed or small-leaf vegetation [3].However, the NDVI can often be saturated across densely vegetated regions or during high-growth periods [4], whereas the leaf area index (LAI) has been shown to be more effective for drought monitoring [5,6].Previous studies have shown that the multi-temporal accumulation precipitation index has a strong correlation with vegetation response [7][8][9].The evolution of meteorological hydrological factors affects vegetation growth [10].Both factors-precipitation and runoff-that cause variations in the vegetation LAI, should be examined separately.
Dynamic watershed runoff is essential to the water supply, irrigation management, and accurate flood prediction and control.River discharge at a site is an integrated signal of water cycle processes over the catchment.However, many catchments lack runoff observation.A hydrological model is an important tool to understand water cycle processing and to fill the gaps of runoff information.The simulation accuracy mainly depends on model mechanisms and data input.The three-layer variable infiltration capacity (VIC) model [11] is a hydrologically based land surface scheme.Spatial variations of infiltration, precipitation, and vegetation are partially considered in the model to simulate water and energy budgets on the land surface.The VIC model has been successfully applied and evaluated in many major river basins within the U.S. [12], the Mekong River Basin [13,14], and on a global scale [15].The quality and coverage of precipitation observation greatly influence the accuracy of hydrological modeling [16].Satellite is an alternative that estimates precipitation with spatial and temporal continuity.Demand is increasing for rainfall estimation from satellites.The Tropical Rainfall Measuring Mission (TRMM) and the Global Precipitation Measurement (GPM) mission have been successively launched to provide global coverage of the Earth precipitation on a sub-daily basis.The products have been used to estimate rainfall over India [17], China [18,19], Eastern Asia [20], and Mekong river basin [21].The versions of precipitation products are updated frequently.It is necessary to make a full comparison of the TRMM and continuous versions of GPM products based on the local gauge observations.
In this paper based on multiple remote-sensing precipitation products, we compare runoff simulations using the VIC hydrological model.The Mun-chi River Basin in Thailand is famous for the origin of aromatic rice.The flat topography makes the basin suitable for simulation comparisons.The main objectives are as follows: (1) To compare the accuracy of different precipitation products (TRMM, GPM, and WMO interpolated precipitation) with an independent station dataset; (2) to validate runoff simulations at 29 hydrological stations using the VIC model; (3) to analyze the influence of runoff on Moderate-resolution Imaging Spectroradiometer (MODIS) LAI variations during dry seasons.

Study Area
Dominated by the Korat Plateau, the Mun-chi River Basin has a flat terrain varying in height between 170 m and 300 m, with boundary mountain ranges between 500 m and 1000 m.The drainage area of the Chi and Mun rivers covers roughly 120,000 km 2 of the Mekong River Basin (Figure 1a).The 30-year runoff at outlets for the Chi and Mun rivers is 291 and 624 m 3 /s (Table 1).The multi-year mean precipitation increases from west to east, ranging from 1100 to 1600 mm (Figure 1b).During the rainy season beginning in May and ending in October, appreciable rainfall comes throughout the whole region.The weather becomes dry and cold from November to February due to the northeast monsoon.The temperature peaks in March and April and drops to its annual lows in November and December.Potential evaporation is constantly high with an annual average of 1600 mm.Soils are mostly sandy with high rates of percolation, low organic matter, and low soil fertility, which constrains on-farm water management for agricultural production and irrigation development.Using high-resolution remote-sensing technology, the planting area of rice, sugarcane, and cassava occupies 45.7%, 13.4%, and 8.2% of the basin, respectively.Rice is planted in the middle and downstream of the Mun-chi River Basin.The main growth period for rice is from June to December, and in irrigated regions, local farmers plant rice from January to April during dry seasons.Cassava is mainly planted in the southwest and north of the basin.The sugarcane is planted in the southwest and north of the basin.Sugarcane is planted in the northwest and north of the basin.Cassava and sugarcane can grow throughout the entire year.

Meteorological forcing dataset for hydrological model
The daily rainfall, maximum temperature, minimum temperature, and average wind speed of 12 weather stations within the Mun-chi River Basin (2002-2017) (Figure 1b) were provided by the National Oceanic and Atmospheric Administration (NOAA) Global Daily Climatology Network [22] (ftp://ftp.ncdc.noaa.gov/pub/data/gsod/).The Thailand data are from the National Center for Atmospheric Research, by the Global Telecommunication System of the WMO.The rainfall dataset summarizes rain and/or melted snow reports during the day in inches.We converted the units to millimeters.There are many interpolation methods to obtain values at positions between the stations.The inverse distance weighting (IDW) method [23] is simple and widely used [24,25] and assumes that closer values are more related than further values.The topography of the Mun-chi river basin is flat, and the weather stations are relatively homogeneously distributed.Therefore, we employed the IDW method to interpolate the daily rainfall, maximum temperature, minimum temperature, and average wind speed data into a 0.1° grid dataset covering the whole basin (Figure 1b).

Hydrological rainfall dataset
We obtained hourly precipitation records (2014-2017) at 16 rain gauges around the upstream Ubolratana dam (UB) and 15 rain gauges around the downstream Sirindhorn dam (SI) from the rainfall telemetry systems of the Electricity Generating Authority of Thailand (EGAT) http://watertele.egat.co.th/.

• Meteorological forcing dataset for hydrological model
The daily rainfall, maximum temperature, minimum temperature, and av 12 weather stations within the Mun-chi River Basin (2002-2017) (Figure 1b) w National Oceanic and Atmospheric Administration (NOAA) Global Daily Clim (ftp://ftp.ncdc.noaa.gov/pub/data/gsod/).The Thailand data are from the Atmospheric Research, by the Global Telecommunication System of the WMO summarizes rain and/or melted snow reports during the day in inches.We c millimeters.There are many interpolation methods to obtain values at positions The inverse distance weighting (IDW) method [23] is simple and widely used that closer values are more related than further values.The topography of the flat, and the weather stations are relatively homogeneously distributed.Theref IDW method to interpolate the daily rainfall, maximum temperature, minim average wind speed data into a 0.1° grid dataset covering the whole basin (Figu

Meteorological forcing dataset for hydrological model
The daily rainfall, maximum temperature, minimum temperature, and average wind speed of 12 weather stations within the Mun-chi River Basin (2002-2017) (Figure 1b) were provided by the National Oceanic and Atmospheric Administration (NOAA) Global Daily Climatology Network [22] (ftp: //ftp.ncdc.noaa.gov/pub/data/gsod/).The Thailand data are from the National Center for Atmospheric Research, by the Global Telecommunication System of the WMO.The rainfall dataset summarizes rain and/or melted snow reports during the day in inches.We converted the units to millimeters.There are many interpolation methods to obtain values at positions between the stations.The inverse distance weighting (IDW) method [23] is simple and widely used [24,25] and assumes that closer values are more related than further values.The topography of the Mun-chi river basin is flat, and the weather stations are relatively homogeneously distributed.Therefore, we employed the IDW method to interpolate the daily rainfall, maximum temperature, minimum temperature, and average wind speed data into a 0.1 • grid dataset covering the whole basin (Figure 1b).The daily rainfall, maximum temperature, minimum temperature, and av 12 weather stations within the Mun-chi River Basin (2002-2017) (Figure 1b) w National Oceanic and Atmospheric Administration (NOAA) Global Daily Clim (ftp://ftp.ncdc.noaa.gov/pub/data/gsod/).The Thailand data are from the Atmospheric Research, by the Global Telecommunication System of the WMO summarizes rain and/or melted snow reports during the day in inches.We c millimeters.There are many interpolation methods to obtain values at positions The inverse distance weighting (IDW) method [23] is simple and widely used that closer values are more related than further values.The topography of the flat, and the weather stations are relatively homogeneously distributed.Theref IDW method to interpolate the daily rainfall, maximum temperature, minim average wind speed data into a 0.1° grid dataset covering the whole basin (Figu

Hydrological rainfall dataset
We obtained hourly precipitation records (2014-2017) at 16 rain gauges around the upstream Ubolratana dam (UB) and 15 rain gauges around the downstream Sirindhorn dam (SI) from the rainfall telemetry systems of the Electricity Generating Authority of Thailand (EGAT) http://watertele.egat.co.th/.

Meteorological Materials •
Meteorological forcing dataset for hydrological model The daily rainfall, maximum temperature, minimum temperature, and av 12 weather stations within the Mun-chi River Basin (2002-2017) (Figure 1b) w National Oceanic and Atmospheric Administration (NOAA) Global Daily Clim (ftp://ftp.ncdc.noaa.gov/pub/data/gsod/).The Thailand data are from the Atmospheric Research, by the Global Telecommunication System of the WMO summarizes rain and/or melted snow reports during the day in inches.We c millimeters.There are many interpolation methods to obtain values at positions The inverse distance weighting (IDW) method [23] is simple and widely used that closer values are more related than further values.The topography of the flat, and the weather stations are relatively homogeneously distributed.Theref IDW method to interpolate the daily rainfall, maximum temperature, minim average wind speed data into a 0.1° grid dataset covering the whole basin (Figu • Hydrological rainfall dataset We obtained hourly precipitation records (2014-2017) at 16 rain gauges Ubolratana dam (UB) and 15 rain gauges around the downstream Sirindho rainfall telemetry systems of the Electricity Generating Authority o http://watertele.egat.co.th/.
Remote sensing precipitation dataset (a) TRMM precipitation product TRMM precipitation product provides global precipitation estimates from a wide variety of meteorological satellites.The Core Observatory satellite, TRMM, a type of low earth orbit (LEO) satellite, was launched in 1997 through cooperation between the National Aeronautics and Space Administration (NASA) and the Japan Aerospace Exploration Agency (JAXA) to analyze and understand the characteristics of precipitation.The altitude of TRMM was 350 km and later changed to 402.5 km in August 2001.TRMM precipitation radar (PR) and TRMM Microwave Imager (TMI) were employed to calibrate precipitation estimates derived from available low earth orbit (LEO) microwave radiometers.TRMM 3B42 V7 products were estimated by combining geostationary-infrared (Geo-IR) and the various microwave data (TMI, TRMM Combined Instrument (TCI), Special Sensor Microwave Imager (SSMI), Advanced Microwave Scanning Radiometer-Earth Observing System (AMSR-E), Advanced Microwave Sounding Unit (AMSU), Special Sensor Microwave Imager/Sounder (SSMI/S), Microwave Humidity Sounder (MHS)) following four steps: (1) Estimate the optimal values from the microwave observations using the Goddard Profiling (GPROF) 2010 algorithm; (2) create Geo-IR precipitation estimates using the optimal values from the previous step; (3) combine the microwave observations and Geo-IR estimates inter-compatibly; (4) rescale and calibrate the merged values using the gauge observations [26].The TRMM Multi-satellite Precipitation 3-hour 0.25 • TRMM_3B42 version 07 (1 January 2002-31 December 2017) was downloaded from https://mirador.gsfc.nasa.gov.We resampled the 0.25 • grid products into 0.1 • latitude-longitude grids, which are in accordance with other precipitation products.We interpolated and calculated the hourly, daily, and monthly datasets based on the TRMM 3-hour products.

(b) GPM IMERG precipitation product
The GPM mission, which was launched on 27 Februay 2014, by NASA and JAXA, is an international project to provide global precipitation observation.The main purpose of the GPM mission is to establish the structure and magnitude of variations in precipitation in order to elucidate the water and energy cycle.GPM was also a constellation of multiple satellites.The microwave imager and dual-frequency precipitation radar of GPM Core Observatory satellite serve as a reference standard to combine precipitation measurements from all satellites that fly within a particular constellation.Additionally, those data measurements can be combined with other data to improve accuracy and reliability.Among the various GPM products, the GPM-3IMERGHH "Final" product represented a half-hour span and 0.1 • grid spatial resolution.The algorithm for estimating GPM-3IMERGHH unites various satellite retrievals from TMPA of NASA Goddard Space Flight Center (GSFC), CMORPH at National Oceanic and Atmospheric Administration (NOAA), and PERSIANN at the University of California, Irvine [27].First, all input datasets were processed using GPROF 2014, an improved multi-channel physical algorithm for retrieving satellite-based rainfall and vertical structure information from various satellite-based microwave observations (GMI, SSMIS, AMSR-2, and MHS).Second, Climate Prediction Center Morphing-Kalman Filter (CMORPH-KF) with the Lagrangian time interpolation technique and the PERSIANN-Cloud Classification System (PERSIANN-CCS) re-calibration scheme provide products with finer temporal and spatial resolution.Finally, monthly GPCC data are applied for bias correction to improve the accuracy of GPM-3IMERGHH products.PrecipitationCal estimates provide daily accumulated precipitation estimate combined microwave-IR [28].In this study, we downloaded PrecipitationCal estimates of GPM IMERG half-hourly final products (GPM_3IMERGHHL) version 05B (

•
Meteorological forcing dataset for hydrological model The daily rainfall, maximum temperature, minimum temperature, and av 12 weather stations within the Mun-chi River Basin (2002-2017) (Figure 1b) w National Oceanic and Atmospheric Administration (NOAA) Global Daily Clim (ftp://ftp.ncdc.noaa.gov/pub/data/gsod/).The Thailand data are from the Atmospheric Research, by the Global Telecommunication System of the WMO summarizes rain and/or melted snow reports during the day in inches.We c millimeters.There are many interpolation methods to obtain values at positions The inverse distance weighting (IDW) method [23] is simple and widely used that closer values are more related than further values.The topography of the flat, and the weather stations are relatively homogeneously distributed.Theref IDW method to interpolate the daily rainfall, maximum temperature, minim average wind speed data into a 0.1° grid dataset covering the whole basin (Figu

Land cover and soil dataset
The MODIS Land Cover Type product (MCD12Q1.051)provides a 500-m grid dataset characterizing University of Maryland (UMD) global land cover classification systems.We downloaded the dataset in 2013 from https://e4ftl01.cr.usgs.gov/MOTA/MCD12Q1.051,maintained by the NASA EOSDIS Land Processes Distributed Active Archive Center (LP DAAC) at the USGS Earth Resources Observation and Science (EROS) Center, Sioux Falls, South Dakota.We also obtained the Regridded Harmonized World Soil Database v1.2 [29] from The Oak Ridge National Laboratory Distributed Active Archive Center.Using the sand, clay, gravel, organic carbon, bulk density, and electrical conductivity data, we defined 0-10 cm, 10-30 cm, and 30-70 cm three-soil layers and prepared parameters based on the Saxton-Rawls equation solutions for soil water characteristics (https://data.nal.usda.gov/).VIC [11] is a macroscale hydrologic model that solves full water and energy balances, originally developed by Xu Liang at the University of Washington.We referred to the study case of the VIC model in neighboring regions [30,31] and assigned the four soil parameter values as follows: The Ds max is the maximum baseflow that can occur from the lowest soil layer, which is estimated based on the calculated saturated hydraulic conductivity at each grid.b infilt defines the shape of the VIC curve, to which we assigned the value 0.2.D s is the fraction of Ds max where non-linear (rapidly increasing) baseflow begins, to which we assigned the value 0.4, and W s is the fraction of the maximum soil moisture (of the lowest soil layer), where non-linear baseflow occurs, to which we assigned the value 0.9.The values of b infilt , D s and W s are applied to each grid within the entire region of the Mun-chi River Basin.

Meteorological forcing dataset for hydrological model
The daily rainfall, maximum temperature, minimum temperature, and av 12 weather stations within the Mun-chi River Basin (2002-2017) (Figure 1b) w National Oceanic and Atmospheric Administration (NOAA) Global Daily Clim (ftp://ftp.ncdc.noaa.gov/pub/data/gsod/).The Thailand data are from the Atmospheric Research, by the Global Telecommunication System of the WMO summarizes rain and/or melted snow reports during the day in inches.We c millimeters.There are many interpolation methods to obtain values at positions The inverse distance weighting (IDW) method [23] is simple and widely used that closer values are more related than further values.The topography of the flat, and the weather stations are relatively homogeneously distributed.Theref IDW method to interpolate the daily rainfall, maximum temperature, minim average wind speed data into a 0.1° grid dataset covering the whole basin (Figu •

Hydrological rainfall dataset
We obtained hourly precipitation records (2014-2017) at 16 rain gauges Ubolratana dam (UB) and 15 rain gauges around the downstream Sirindho rainfall telemetry systems of the Electricity Generating Authority o http://watertele.egat.co.th/.The algorithm of the MODIS LAI product MCD15A2H.006chooses the best pixel available from all the acquisitions of both MODIS sensors located on NASA Terra and Aqua satellites from within the eight-day period.We assumed the daily LAI value was constant during the eight days.

Runoff Data
From the Royal Irrigation Department, we obtained monthly runoff data from 29 hydrological stations (Figure 1b

Evaluation of Precipitation Products Accuracy
We compared the GPM IMERG versions 3-5, TRMM version 7, and the WMO interpolated grid data against daily and monthly hydrological rainfall observation data.We employed the two-dimensional space of the Taylor diagram [32], which represents three statistics simultaneously: The centered root-mean-square difference, the correlation coefficient, and the standard deviation.The Taylor diagram evaluates products from the aspect of precipitation variation.Meanwhile, we compared the total amount of precipitation among WMO, TRMM, and IMERG V5.We finally compared the four-year average hourly hydrological rainfall observations with the TRMM and IMERG V5 products during 2014 to 2017 as the necessary supplements.

Validation of VIC Model
The VIC model version 4.2 was employed for this study.One of the new features of this model is that monthly climatological vegetation parameters were replaced with daily time series in the vegetation forcing file [33].We selected computing water balance only option and simulated daily moisture and energy fluxes and runoff from August 2002 to December 2017 at a 0.1-degree grid, based on the three precipitation datasets:

Vegetation and Drought Indices
The previous drought studies [9, 34,35] showed that the remote-sensing-based NDVI anomaly is a useful index to indicate vegetation drought.Similarly, we used the LAI anomaly (VA) to depict the vegetation variation in equation (1).We used the runoff anomaly (RA) to depict the vegetation water stress in Equation ( 2).In these equations, i and j are the specific year and month of the data, respectively.In this study, we employed a cumulative precipitation anomaly [36] to depict precipitation drought.Cumulative precipitation anomaly has advantage in robustness, dimensionality, and transparency, compared with the classic Palmer drought severity index [37].The standardized precipitation index observed cumulative precipitation deviates from the average.However, it requires a long-term time series climatological average of precipitation accumulations over the desired time scale to estimate an appropriate probability density function [38].The time series of current remote sensing products in this study were no longer enough to make accurate function fitness.P k,n is cumulative precipitation at nmonth time scale.k is ascending time serial number of cumulative precipitation.The ranges of n equal 1, 3, and 6, and we calculated the one-, three-, and six-month cumulative precipitation based on Equation (3).
VA i,j = LAI i,j − LAI ave,j LAI ave,j PA i,j = P i,j − P ave,j P ave,j

Accuracy Evaluation of Precipitation Estimations
We compared the rainfall products according to the correlation coefficients, standard deviation, and multi-year average value.From the Taylor diagram, the correlation coefficients of the daily and monthly precipitation products during the rainy season are higher than those during the dry season.Among the correlation coefficients of different products, the correlation of daily IMERG V3-V5 and TRMM products, which are from 0.55 to 0.65, were higher than the WMO estimation of 0.4 (Figure 2a).The IMERG V5 daily product showed the strongest correlation of 0.626, and the standard deviation was close to 1. Similarly, the correlation coefficients of four monthly satellite products (>0.8) were higher than the WMO interpolation of 0.777 (Figure 2b).The TRMM monthly products show the strongest correlation of 0.843.We further compared the precipitation products around two dams.Figure 2c shows that the correlation coefficients of all precipitation products around the SI dam in the downstream Mun-chi river are higher than those around the UB dam in the upstream area.We examined the average precipitation per day with each gauge.The bias of average precipitation of all gauges based on TRMM, IMERG V3, IMERG V4, IMERG V5, and WMO was +1.21 mm/day, +1.20 mm/day, +1.23 mm/day, +1.21 mm/day, and +0.85 mm/day, respectively.TRMM and IMERG showed higher overestimations of the observed precipitation in the Mun-chi River Basin than the WMO interpolation.From the comparison of scatter charts among TRMM, IMERG V5, and WMO (Figure 2d), the dots of TRMM and IMERG V5 precipitation almost paralleled the 1:1 line but overestimated the gauge-based precipitation.In addition, we examined the hourly precipitation of IMERG V5, TRMM, and the gauge-based observations.The results showed similar precipitation diurnal cycles in both the upstream UB dam and downstream SI dam (Figure 3).The times with less precipitation in the UB dam by TRMM were 06:00 to 13:00, by IMERG V5 was 10:00, and by the gauge-based observations was 0 8:00; in the SI dam, the times by TRMM were 09:00 to 13:00, by IMERG V5 was 11:00, and by the gauge-based In addition, we examined the hourly precipitation of IMERG V5, TRMM, and the gauge-based observations.The results showed similar precipitation diurnal cycles in both the upstream UB dam and downstream SI dam (Figure 3).The times with less precipitation in the UB dam by TRMM were 06:00 to 13:00, by IMERG V5 was 10:00, and by the gauge-based observations was 08:00; in the SI dam, the times by TRMM were 09:00 to 13:00, by IMERG V5 was 11:00, and by the gauge-based observations was 09:00.We compared the time with average hourly precipitation above 40 mm.In the UB dam, the period was from 00:00 to 02:00 by gauged observation, 00:00 to 03:00 by TRMM, and 00:00 to 06:00 by IMERG V5.In the SI dam, the period was from 13:00 to 08:00 by gauged observation, 13:00 to 06:00 by TRMM, and 13:00 to 09:00 by IMERG V5.The period by IMERG V5 was longer than by the gauged observations.IMERG V5 and TRMM products both overestimated the gauged-based rainfall observations at the UB and SI dams, especially in the afternoon.In all, the IMERG V5 and TRMM satellite products overestimated the rainfall, possibly due to systematic deviations.However, from the aspect of correlation and temporal variation, the satellite precipitation products have their advantages.
Water 2019, 11, x FOR PEER REVIEW 8 of 14 observations was 09:00.We compared the time with average hourly precipitation above 40 mm.In the UB dam, the period was from 00:00 to 02:00 by gauged observation, 00:00 to 03:00 by TRMM, and 00:00 to 06:00 by IMERG V5.In the SI dam, the period was from 13:00 to 08:00 by gauged observation, 13:00 to 06:00 by TRMM, and 13:00 to 09:00 by IMERG V5.The period by IMERG V5 was longer than by the gauged observations.IMERG V5 and TRMM products both overestimated the gauged-based rainfall observations at the UB and SI dams, especially in the afternoon.In all, the IMERG V5 and TRMM satellite products overestimated the rainfall, possibly due to systematic deviations.However, from the aspect of correlation and temporal variation, the satellite precipitation products have their advantages.

Validation of Runoff Simulation by VIC Model
The correlation (R 2 ) map of monthly runoff with observations from 2003 to March 2008 show that the average correlation (R 2 ) during rainy seasons (Figure 4a,b) 0.493-0.530are higher than those during dry seasons 0.452-0.478(Figure 4c,d).The average R 2 at 22 stations in tributaries for Figure 4 are 0.495, 0.540, 0.430, and 0.435, respectively.TRMM greatly improves the runoff accuracy in tributaries.The average R 2 at seven mainstream stations are 0.488, 0.497, 0.629, 0.507.WMO is better than TRMM during dry season (Figure 5a).E.21 and E.9 upstream stations with similar drainage areas have similar correlations (R 2 ).E.18, M.5, and E.20A with similar drainage areas have similar correlations (R 2 ).The correlation (R 2 ) for the E.21. and E.9 upstream stations are higher than for M.7, E.18, M.5, and E.20A.The R 2 downstream values show a slight increase with the expansion of the drainage area.However, the simulated runoff at M.6A has a very weak correlation with the observations.At the M.7 hydrological station, the Nash-Sutcliffe efficiency coefficients (NSE) of monthly simulated runoff from 2003 to March 2008 by the TRMM and WMO (Figure 6a) were 0.509 and 0.628, respectively.From April 2014 to March 2017, the NSE of monthly simulated runoff by the TRMM, IMERG V5, and WMO (Figure 6b) were 0.041, −0.166, and 0.508, respectively.The WMO demonstrated a relatively higher simulation accuracy than the TRMM.There is a time delay of the flood peak of measured daily runoff compared to the simulations (Figure 6c).The NSE of daily simulated runoff by the TRMM, IMERG V5, and WMO were −0.354, −0.984, and 0.066, respectively.

The Influence of Runoff on Vegetation Variation
We analyzed the influences of runoff and accumulated precipitation on vegetation.The differences between the WMO and TRMM precipitation products had little impact on our evaluation (Figure 7).From January 2003 to March 2008, the correlations of the LAI with runoff and precipitation show that the measured runoff has the highest correlation (R 2 = 0.529) with vegetation (Figure 7a).The six-month precipitation has similarly higher correlations (R 2 = 0.546 and 0.563) with vegetation (Figure 7h,i).On the contrary, there are poorer correlations between the one-month precipitation and the LAI (Figure

Validation of Runoff Simulation by VIC Model
The correlation (R 2 ) map of monthly runoff with observations from 2003 to March 2008 show that the average correlation (R 2 ) during rainy seasons (Figure 4a,b) 0.493-0.530are higher than those during dry seasons 0.452-0.478(Figure 4c,d).The average R 2 at 22 stations in tributaries for Figure 4 are 0.495, 0.540, 0.430, and 0.435, respectively.TRMM greatly improves the runoff accuracy in tributaries.The average R 2 at seven mainstream stations are 0.488, 0.497, 0.629, 0.507.WMO is better than TRMM during dry season (Figure 5a).E.21 and E.9 upstream stations with similar drainage areas have similar correlations (R 2 ).E.18, M.5, and E.20A with similar drainage areas have similar correlations (R 2 ).The correlation (R 2 ) for the E.21. and E.9 upstream stations are higher than for M.7, E.18, M.5, and E.20A.The R 2 downstream values show a slight increase with the expansion of the drainage area.However, the simulated runoff at M.6A has a very weak correlation with the observations.At the M.7 hydrological station, the Nash-Sutcliffe efficiency coefficients (NSE) of monthly simulated runoff from 2003 to March 2008 by the TRMM and WMO (Figure 6a) were 0.509 and 0.628, respectively.From April 2014 to March 2017, the NSE of monthly simulated runoff by the TRMM, IMERG V5, and WMO (Figure 6b) were 0.041, −0.166, and 0.508, respectively.The WMO demonstrated a relatively higher simulation accuracy than the TRMM.There is a time delay of the flood peak of measured daily runoff compared to the simulations (Figure 6c).The NSE of daily simulated runoff by the TRMM, IMERG V5, and WMO were −0.354, −0.984, and 0.066, respectively.the WMO and TRMM are higher than the measured runoff (Figure 7j-l).The correlations are similar to the results by six-month precipitation anomaly (Figure 7q,r).The one-month precipitation also has poor correlations with LAI variations.

The Comparison of Precipitation Products
Compared with the estimations in downstream tributaries, the upstream monitoring of precipitation is not highly accurate, which may be due to complex, hilly topography.It also suggests that the dataset from the hydrological rain gauge reflects the heterogenicity of precipitation and is suitable to validate precipitation products.We found that daily IMERG V5 has the highest accuracy; this is consistent with the evaluation studies [17,19], which show that IMERG estimations exhibit notable improvements over TRMM to detect heavy storms and precipitation at sub-daily and daily resolution.A study in the Nan river basin of Thailand [39] revealed that mean bias of TRMM and GPM was about +4.5 mm/day and +1.54 mm/day, respectively.Our study reveals higher accuracies of TRMM and GPM products in northeastern Thailand.Our study shows that the bias of GPM is higher than TRMM.TRMM and GPM are global precipitation products; a validation study across different climate zones [40] showed that the TRMM satellite product overestimates heavy and extreme daily precipitation for some areas of humid regions but underestimates light and moderate precipitation over arid regions.IMERG performs better than the TRMM product at the mid-and high-latitudes, as well as in relatively dry climate regions [19].Northeastern Thailand is located in a humid climate zone.Our results are comparable with these previous studies.Using the TRMM precipitation radar, previous research [41] found that early morning rainfall significantly contributes to the climatological rainfall pattern.Future regional products integrated with a high density of ground-based rainfall data are expected to correct the systematic deviations.Our study in northeastern Thailand using hourly GPM IMERG products and rain gauge observations validated the precipitation diurnal pattern, revealing the potential of remote-sensing precipitation products for timely flash flood monitoring.

The Influence of Runoff on Vegetation Variation
We analyzed the influences of runoff and accumulated precipitation on vegetation.The differences between the WMO and TRMM precipitation products had little impact on our evaluation (Figure 7).From January 2003 to March 2008, the correlations of the LAI with runoff and precipitation show that the measured runoff has the highest correlation (R 2 = 0.529) with vegetation (Figure 7a).The six-month precipitation has similarly higher correlations (R 2 = 0.546 and 0.563) with vegetation (Figure 7h,i).On the contrary, there are poorer correlations between the one-month precipitation and the LAI (Figure 7d,e).From April 2014 to March 2017, the overall correlations decrease.The correlations of runoff by the WMO and TRMM are higher than the measured runoff (Figure 7j-l).The correlations are similar to the results by six-month precipitation anomaly (Figure 7q,r).The one-month precipitation also has poor correlations with LAI variations.

The Evaluation of Runoff Simulations
A previous study [13] showed that the ratio of reservoir storage capacity to the mean annual streamflow is low in many Asian river basins, and the shape of the hydrographs does not change much.Our study shows that in the Mun-chi River Basin, the simulated natural runoff is close to the measured data for the monthly scale; both studies are comparable.WMO performs better than TRMM and IMERG to simulate runoff.In a previous study [42] using the GPM IMERG data for hydrological application in the Ganjiang basin, rain gauge-interpolated products exhibited better performance than the newly available GPM IMERG product.We also found that TRMM can improve the runoff simulation in downstream tributaries but not the downstream mainstream runoff.One possible reason is the lower accuracy of WMO interpolation in downstream tributaries.
Another reason is that the runoff of downstream mainstream is influenced by runoff variation from nearby and upstream tributaries.It seems that the influences of precipitation fluctuations on runoff were smoothed by expanding the drainage area.The simulated runoff at the M.6A stations had a poor correlation with the observations.According to the FAO dam material [43], the total capacity of 11 dams in the Mun-chi River Basin is 4643.86 million m 3 .There are five large reservoirs

The Comparison of Products
Compared with the estimations in downstream tributaries, the upstream monitoring of precipitation is not highly accurate, which may be due to complex, hilly topography.It also suggests that the dataset from the hydrological rain gauge reflects the heterogenicity of precipitation and is suitable to validate precipitation products.We found that daily IMERG V5 has the highest accuracy; this is consistent with the evaluation studies [17,19], which show that IMERG estimations exhibit notable improvements over TRMM to detect heavy storms and precipitation at sub-daily and daily resolution.A study in the Nan river basin of Thailand [39] revealed that mean bias of TRMM and GPM was about +4.5 mm/day and +1.54 mm/day, respectively.Our study reveals higher accuracies of TRMM and GPM products in northeastern Thailand.Our study shows that the bias of GPM is higher than TRMM.TRMM and GPM are global precipitation products; a validation study across different climate zones [40] showed that the TRMM satellite product overestimates heavy and extreme daily precipitation for some areas of humid regions but underestimates light and moderate precipitation over arid regions.IMERG performs better than the TRMM product at the mid-and high-latitudes, as well as in relatively dry climate regions [19].Northeastern Thailand is located in a humid climate zone.Our results are comparable with these previous studies.Using the TRMM precipitation radar, previous research [41] found that early morning rainfall significantly contributes to the climatological rainfall pattern.Future regional products integrated with a high density of ground-based rainfall data are expected to correct the systematic deviations.Our study in northeastern Thailand using hourly GPM IMERG products and rain gauge observations validated the precipitation diurnal pattern, revealing the potential of remote-sensing precipitation products for timely flash flood monitoring.

The Evaluation of Runoff Simulations
A previous study [13] showed that the ratio of reservoir storage capacity to the mean annual streamflow is low in many Asian river basins, and the shape of the hydrographs does not change much.Our study shows that in the Mun-chi River Basin, the simulated natural runoff is close to the measured data for the monthly scale; both studies are comparable.WMO performs better than TRMM and IMERG to simulate runoff.In a previous study [42] using the GPM IMERG data for hydrological application in the Ganjiang basin, rain gauge-interpolated products exhibited better performance

Figure 1 .
Figure 1.Mun-chi river basin (a) elevation and hydrological stations and (b) rainfall gauges and landuse types.

Figure 1 .
Figure 1.Mun-chi river basin (a) elevation and hydrological stations and (b) rainfall gauges and land-use types.

Figure 1 .
Figure 1.Mun-chi river basin (a) elevation and hydrological stations and (b) rainfa use types.

•
Meteorological forcing dataset for hydrological model

Figure 1 .
Figure 1.Mun-chi river basin (a) elevation and hydrological stations and (b) rainfa use types.

Figure 1 .
Figure 1.Mun-chi river basin (a) elevation and hydrological stations and (b) rainfa use types.

Figure 1 .
Figure 1.Mun-chi river basin (a) elevation and hydrological stations and (b) rainfa use types.
Shortwave albedo and Vegetation LAI datasetThe daily black sky albedo product (MCD43A3.006)(2001-2017) and the eight-day LAI product (MCD15A2H.006)(July 2002 to 2017) were retrieved from the online AppEEARS, courtesy of the NASA EOSDIS Land Processes Distributed Active Archive Center, USGS/Earth Resources Observation and Science Center, Sioux Falls, South Dakota, USA and downloaded from https://lpdaacsvc.cr.usgs.gov/appeears.Due to cloud contamination, there are low-quality pixels in the images.The images of black sky albedo data on 31 May 2001, 24-26 June 2001, 20 June 2014, and 24 August 2014 are missing.
) from April 2002 to March 2008 and daily runoff data at the basin outlet M7 station (104.86 • E 15.22 • N) from April 2014 to March 2017.

( 1 )
WMO interpolated precipitation; (2) WMO + GPM IMERG V5 precipitation (we employed daily GPM IMERG V5 to replace the WMO precipitation after 12 March 2014); and (3) daily TRMM products from 2002 to 2017.The period from August 2002 to December 2002 was used for model spin-up.The period from January 2003 to March 2008 was used for calibration, and the period from April 2014 to March 2017 was used for model validation.

Figure 2 .
Figure 2. Taylor diagram of daily and monthly statistics among the multi-satellite retrievals for Global Precipitation Measurement (IMERG) V03, IMERG V04, IMERG V05, Tropical Rainfall Measuring Mission (TRMM), and World Meteorological Organization (WMO) interpolation grids from 12 March 2014 to 31 January 2016: (a) Daily precipitation during rainy season, dry season, and for the whole period.(b) Monthly precipitation during rainy season, dry season, and for the whole period; (c) monthly precipitation around two dams; (d) total precipitation at each gauge.

Figure 2 .
Figure 2. Taylor diagram of daily and monthly statistics among the multi-satellite retrievals for Global Precipitation Measurement (IMERG) V03, IMERG V04, IMERG V05, Tropical Rainfall Measuring Mission (TRMM), and World Meteorological Organization (WMO) interpolation grids from 12 March 2014 to 31 January 2016: (a) Daily precipitation during rainy season, dry season, and for the whole period.(b) Monthly precipitation during rainy season, dry season, and for the whole period; (c) monthly precipitation around two dams; (d) total precipitation at each gauge.
11,  x FOR PEER REVIEW 9 of 14 the WMO and TRMM are higher than the measured runoff (Figure7j-l).The correlations are similar to the results by six-month precipitation anomaly (Figure7q,r).The one-month precipitation also has poor correlations with LAI variations.

Figure 4 .
Figure 4.The correlation map of runoff simulation based on (a) WMO during rainy season; (b) TRMM during rainy season; (c) WMO during dry season; (d) TRMM during dry season from 2003 to 2008.

Figure 5 .
Figure 5.The monthly runoff at seven hydrological stations in the mainstream of Mun-chi river from 2003 to March 2008: (a) The correlation of runoff simulations based on WMO and TRMM precipitation river; (b) the comparison between of measured runoff at M.6A and M.5 in 2005 hydrological year.

Figure 4 .
Figure 4.The correlation map of runoff simulation based on (a) WMO during rainy season; (b) TRMM during rainy season; (c) WMO during dry season; (d) TRMM during dry season from 2003 to 2008.

Figure 4 .
Figure 4.The correlation map of runoff simulation based on (a) WMO during rainy season; (b) TRMM during rainy season; (c) WMO during dry season; (d) TRMM during dry season from 2003 to 2008.

Figure 5 .
Figure 5.The monthly runoff at seven hydrological stations in the mainstream of Mun-chi river from 2003 to March 2008: (a) The correlation of runoff simulations based on WMO and TRMM precipitation river; (b) the comparison between of measured runoff at M.6A and M.5 in 2005 hydrological year.

Figure 5 .
Figure 5.The monthly runoff at seven hydrological stations in the mainstream of Mun-chi river from 2003 to March 2008: (a) The correlation of runoff simulations based on WMO and TRMM precipitation river; (b) the comparison between of measured runoff at M.6A and M.5 in 2005 hydrological year.

Figure 6 .
Figure 6.The comparison among simulated and measured runoff: (a) Monthly runoff from 2003 to 2008 for calibration; (b) monthly runoff from 2014 to 2017 for validation and, (c) daily runoff from 2014 to 2017 at M.7 hydrological station of Mun-chi river.

Figure 6 .
Figure 6.The comparison among simulated and measured runoff: (a) Monthly runoff from 2003 to 2008 for calibration; (b) monthly runoff from 2014 to 2017 for validation and, (c) daily runoff from 2014 to 2017 at M.7 hydrological station of Mun-chi river.

Table 1 .
The hydrological stations in the mainstream of Mun-chi river Basin.

Table 1 .
The hydrological stations in the mainstream of Mun-chi river Basin.