Ensemble representation of uncertainty in Lagrangian satellite rainfall estimates

A new algorithm called Lagrangian Simulation (LSIM) has been developed that enables the interpolation uncertainty present in Lagrangian satellite rainfall algorithms such as the Climate Prediction Center (CPC) morphing technique (CMORPH) to be characterized using an ensemble product. The new algorithm generates ensemble sequences of rainfall fields conditioned on multiplatform multisensor microwave satellite data, demonstrating a conditional simulation approach that overcomes the problem of discontinuous uncertainty fields inherent in this type of product. Each ensemble member is consistent with the information present in the satellite data, while variation between members is indicative of uncertainty in the rainfall retrievals. LSIM is based on the combination of a Markov weather generator, conditioned on both previous and subsequent microwave measurements, and a global optimization procedure that uses simulated annealing to constrain the generated rainfall fields to display appropriate spatial structures. The new algorithm has been validated over a region of the continental United States and has been shown to provide reliable estimates of both point uncertainty distributions and wider spatiotemporal structures

Bayesian Analysis of the Impact of Rainfall Data Product on Simulated Slope Failure for North Carolina Locations

In the past decades, many different approaches have been developed in the literature to quantify the load-carrying capacity and geotechnical stability (or the factor of safety, Fs) of variably saturated hillslopes. Much of this work has focused on a deterministic characterization of hillslope stability. Yet, simulated Fs values are subject to considerable uncertainty due to our inability to characterize accurately the soil mantles properties (hydraulic, geotechnical, and geomorphologic) and spatiotemporal variability of the moisture content of the hillslope interior. This is particularly true at larger spatial scales. Thus, uncertainty-incorporating analyses of physically based models of rain-induced landslides are rare in the literature. Such landslide modeling is typically conducted at the hillslope scale using gauge-based rainfall forcing data with rather poor spatiotemporal coverage. For regional landslide modeling, the specific advantages and/or disadvantages of gauge-only, radar-merged and satellite-based rainfall products are not clearly established. Here, we compare and evaluate the performance of the Transient Rainfall Infiltration and Grid-based Regional Slope-stability analysis (TRIGRS) model for three different rainfall products using 112 observed landslides in the period between 2004 and 2011 from the North Carolina Geological Survey database. Our study includes the Tropical Rainfall Measuring Mission (TRMM) Multi-satellite Precipitation Analysis Version 7 (TMPA V7), the North American Land Data Assimilation System Phase 2 (NLDAS-2) analysis, and the reference truth Stage IV precipitation. TRIGRS model performance was rather inferior with the use of literature values of the geotechnical parameters and soil hydraulic properties from ROSETTA using soil textural and bulk density data from SSURGO (Soil Survey Geographic database). The performance of TRIGRS improved considerably after Bayesian estimation of the parameters with the DiffeRential Evolution Adaptive Metropolis (DREAM) algorithm using Stage IV precipitation data. Hereto, we use a likelihood function that combines binary slope failure information from landslide event and null periods using multivariate frequency distribution-based metrics such as the false discovery and false omission rates. Our results demonstrate that the Stage IV-inferred TRIGRS parameter distributions generalize well to TMPA and NLDAS-2 precipitation data, particularly at sites with considerably larger TMPA and NLDAS-2 rainfall amounts during landslide events than null periods. TRIGRS model performance is then rather similar for all three rainfall products. At higher elevations, however, the TMPA and NLDAS-2 precipitation volumes are insufficient and their performance with the Stage IV-derived parameter distributions indicates their inability to accurately characterize hillslope stability

Summertime evaluation of REFAME over the Unites States for near real-time high resolution precipitation estimation

Precipitation is the key input for hydrometeorological modeling and applications. In many regions of the world, including populated areas, ground-based measurement of precipitation (whether from radar or rain gauge) is either sparse in time and space or nonexistent. Therefore, high-resolution satellite-based precipitation products are recognized as critical data sources, especially for rapidly-evolving hydrometeorological events such as flash floods which primarily occur during summer/warm seasons. As " proof of concept" , a recently proposed algorithm called Rain Estimation using Forward Adjusted-advection of Microwave Estimates (REFAME) and its variation REFAMEgeo are evaluated over the contiguous United States during summers of 2009 and 2011. Both methods are originally designed for near real-time high resolution precipitation estimation from remotely sensed data. High-resolution Q2 (ground radar) precipitation data, in conjunction with two operational near real-time satellite-based precipitation products (PERSIANN, PERSIANN-CCS) are used as evaluation reference and for comparison. The study is performed at half-hour temporal resolution and at a range of spatial resolutions (0.08-, 0.25-, 0.5-, and 1-degree latitude/longitude). The statistical analyses suggest that REFAMEgeo performs favorably among the studied products in terms of capturing both spatial coverage and intensity of precipitation at near real-time with the temporal resolution offered by geostationary satellites. With respect to volume precipitation, REFAMEgeo together with REFAME demonstrates slight overestimation of intense precipitation and underestimation of light precipitation events. Compared to REFAME, It is observed that REFAMEgeo maintains stable performance, even when the amount of accessible microwave (MW) overpasses is limited. Based on the encouraging outcome of this study which was intended as " proof of concept" , further testing for other seasons and data-rich regions is the next logical step. Upon confirmation of the relative reliability of the algorithm, it is reasonable to recommend the use of its precipitation estimates for data-sparse regions of the world. © 2012 Elsevier B.V

Comparison of daily and sub-daily SWAT models for daily streamflow simulation in the Upper Huai River Basin of China

Despite the significant role of precipitation in the hydrological cycle, few studies have been conducted to evaluate the impacts of the temporal resolution of rainfall inputs on the performance of SWAT (soil and water assessment tool) models in large-sized river basins. In this study, both daily and hourly rainfall observations at 28 rainfall stations were used as inputs to SWAT for daily streamflow simulation in the Upper Huai River Basin. Study results have demonstrated that the SWAT model with hourly rainfall inputs performed better than the model with daily rainfall inputs in daily streamflow simulation, primarily due to its better capability of simulating peak flows during the flood season. The sub-daily SWAT model estimated that 58% of streamflow was contributed by baseflow compared to 34 % estimated by the daily model. Using the future daily and three-hour precipitation projections under the RCP (Representative Concentration Pathways) 4.5 scenario as inputs, the sub-daily SWAT model predicted a larger amount of monthly maximum daily flow during the wet years than the daily model. The differences between the daily and sub-daily SWAT model simulation results indicated that temporal rainfall resolution could have much impact on the simulation of hydrological process, streamflow, and consequently pollutant transport by SWAT models. There is an imperative need for more studies to examine the effects of temporal rainfall resolution on the simulation of hydrological and water pollutant transport processes by SWAT in river basins of different environmental conditions

Bias adjustment of satellite-based precipitation estimation using gauge observations: A case study in Chile

Satellite-based precipitation estimates (SPEs) are promising alternative precipitation data for climatic and hydrological applications, especially for regions where ground-based observations are limited. However, existing satellite-based rainfall estimations are subject to systematic biases. This study aims to adjust the biases in the Precipitation Estimation from Remotely Sensed Information using Artificial Neural Networks–Cloud Classification System (PERSIANN-CCS) rainfall data over Chile, using gauge observations as reference. A novel bias adjustment framework, termed QM-GW, is proposed based on the nonparametric quantile mapping approach and a Gaussian weighting interpolation scheme. The PERSIANN-CCS precipitation estimates (daily, 0.04°×0.04°) over Chile are adjusted for the period of 2009–2014. The historical data (satellite and gauge) for 2009–2013 are used to calibrate the methodology; nonparametric cumulative distribution functions of satellite and gauge observations are estimated at every 1°×1° box region. One year (2014) of gauge data was used for validation. The results show that the biases of the PERSIANN-CCS precipitation data are effectively reduced. The spatial patterns of adjusted satellite rainfall show high consistency to the gauge observations, with reduced root-mean-square errors and mean biases. The systematic biases of the PERSIANN-CCS precipitation time series, at both monthly and daily scales, are removed. The extended validation also verifies that the proposed approach can be applied to adjust SPEs into the future, without further need for ground-based measurements. This study serves as a valuable reference for the bias adjustment of existing SPEs using gauge observations worldwide

MSWEP : 3-hourly 0.25° global gridded precipitation (1979-2015) by merging gauge, satellite, and reanalysis data

Current global precipitation (P) datasets do not take full advantage of the complementary nature of satellite and reanalysis data. Here, we present Multi-Source Weighted-Ensemble Precipitation (MSWEP) version 1.1, a global P dataset for the period 1979-2015 with a 3hourly temporal and 0.25 degrees ffi spatial resolution, specifically designed for hydrological modeling. The design philosophy of MSWEP was to optimally merge the highest quality P data sources available as a function of timescale and location. The long-term mean of MSWEP was based on the CHPclim dataset but replaced with more accurate regional datasets where available. A correction for gauge under-catch and orographic effects was introduced by inferring catchment-average P from streamflow (Q) observations at 13 762 stations across the globe. The temporal variability of MSWEP was determined by weighted averaging of P anomalies from seven datasets; two based solely on interpolation of gauge observations (CPC Unified and GPCC), three on satellite remote sensing (CMORPH, GSMaP-MVK, and TMPA 3B42RT), and two on atmospheric model reanalysis (ERA-Interim and JRA-55). For each grid cell, the weight assigned to the gauge-based estimates was calculated from the gauge network density, while the weights assigned to the satellite-and reanalysis-based estimates were calculated from their comparative performance at the surrounding gauges. The quality of MSWEP was compared against four state-of-the-art gauge-adjusted P datasets (WFDEI-CRU, GPCP-1DD, TMPA 3B42, and CPC Unified) using independent P data from 125 FLUXNET tower stations around the globe. MSWEP obtained the highest daily correlation coefficient (R) among the five P datasets for 60.0% of the stations and a median R of 0.67 vs. 0.44-0.59 for the other datasets. We further evaluated the performance of MSWEP using hydrological modeling for 9011 catchments (< 50 000 km(2)) across the globe. Specifically, we calibrated the simple conceptual hydrological model HBV (Hydrologiska Byrans Vattenbalansavdelning) against daily Q observations with P from each of the different datasets. For the 1058 sparsely gauged catchments, representative of 83.9% of the global land surface (excluding Antarctica), MSWEP obtained a median calibration NSE of 0.52 vs. 0.29-0.39 for the other P datasets. MSWEP is available via http://www.gloh2o.org

LMODEL: A satellite precipitation methodology using cloud development modeling. Part I: Algorithm construction and calibration

The Lagrangian Model (LMODEL) is a new multisensor satellite rainfall monitoring methodology based on the use of a conceptual cloud-development model that is driven by geostationary satellite imagery and is locally updated using microwave-based rainfall measurements from low earth-orbiting platforms. This paper describes the cloud development model and updating procedures; the companion paper presents model validation results. The model uses single-band thermal infrared geostationary satellite imagery to characterize cloud motion, growth, and dispersal at high spatial resolution (similar to 4 km). These inputs drive a simple, linear, semi-Lagrangian, conceptual cloud mass balance model, incorporating separate representations of convective and stratiform processes. The model is locally updated against microwave satellite data using a two-stage process that scales precipitable water fluxes into the model and then updates model states using a Kalman filter. Model calibration and updating employ an empirical rainfall collocation methodology designed to compensate for the effects of measurement time difference, geolocation error, cloud parallax, and rainfall shear