A deterministic optimization approach for solving the rainfall disaggregation problem
Journal: Bulletin of Computational Applied Mathematics (Bull CompAMa) (Vol.3, No. 2)Publication Date: 2015-12-31
Authors : Debora Cores; Lelys Guenni; Lisbeth Torres;
Page : 7-29
Keywords : Rainfall disaggregation; spectral projected gradient; simulated annealing;
Abstract
One of the main problems in hydrology is the time scale of the historical rainfall data, available from many meteorological data bases. Most of the rainfall data is given at a time scale coarser than the one needed for many applications in hydrology and environmental sciences, as the estimation of spatially continuous rainfall at finer time scales, for drainage systems design and extreme rainfall analysis. A method to disaggregate monthly rainfall to daily or finer temporal scale is very important in many applications. Many authors have addressed this problem by using some stochastic methods including several stochastic rainfall models. The lowering resolution methods must be low-cost and low-storage since the amount of rainfall data is large. The purpose of this work is to formulate this problem as a constrained optimization problem and solve it with a low-cost and low-storage deterministic optimization method. We modify the objective function proposed by Guenni and Bárdossy for solving the disaggregation rainfall problem and we use the low-cost spectral projected gradient (SPG) method. In contrast with the stochastic method, a deterministic approach will take into account important information, as for example the gradient of the objective function. The proposed method was applied to a data set from a rainfall network of the central plains of Venezuela, in which rainfall is highly seasonal and data availability at a daily time scale or even higher temporal resolution is very limited. The numerical results show that the SPG method for solving the disaggregation rainfall problem avoids daily precipitations outliers that might occur as an artifact of the simulation procedure and accurately reproduces the probability distribution. Also, the proposed model and methodology outperforms the one proposed by Guenni and Bárdossy (2002) in the sense that it reduces the absolute error value for the statistical properties from the observed data.
Other Latest Articles
- Simulation of the oxygen distribution in a tumor tissue using residual algorithms
- Motion planning algorithms, topological properties and affine approximation
- Matrix completion via a low rank factorization model and an Augmented Lagrangean Succesive Overrelaxation Algorithm
- L(2,1)-Labeling for Subdivisions of Cycle Dominated Graphs
- Uniform Stability and Boundedness of a Kind of Third Order Delay Differential Equations
Last modified: 2018-08-05 10:23:18