scholarly journals Stochastic Modeling of Rainfall in Peninsular Malaysia Using Bartlett Lewis Rectangular Pulses Models

2011 ◽  
Vol 2011 ◽  
pp. 1-6 ◽  
Author(s):  
Ibrahim Suliman Hanaish ◽  
Kamarulzaman Ibrahim ◽  
Abdul Aziz Jemain

Three versions of Bartlett Lewis rectangular pulse rainfall models, namely, the Original Bartlett Lewis (OBL), Modified Bartlett Lewis (MBL), and 2N-cell-type Bartlett Lewis model (BL2n), are considered. These models are fitted to the hourly rainfall data from 1970 to 2008 obtained from Petaling Jaya rain gauge station, located in Peninsular Malaysia. The generalized method of moments is used to estimate the model parameters. Under this method, minimization of two different objective functions which involve different weight functions, one weight is inversely proportional to the variance and another one is inversely proportional to the mean squared, is carried out using Nelder-Mead optimization technique. For the purpose of comparison of the performance of the three different models, the results found for the months of July and November are used for illustration. This performance is assessed based on the goodness of fit of the models. In addition, the sensitivity of the parameter estimates to the choice of the objective function is also investigated. It is found thatBL2nslightly outperformsOBL. However, the best model is the Modified Bartlett LewisMBL, particularly when the objective function considered involves weight which is inversely proportional to the variance.

2018 ◽  
Vol 80 (6) ◽  
Author(s):  
Siti Mariam Saad ◽  
Abdul Aziz Jemain ◽  
Noriszura Ismail

This study evaluates the utility and suitability of a simple discrete multiplicative random cascade model for temporal rainfall disaggregation. Two of a simple random cascade model, namely log-Poisson and log-Normal  models are applied to simulate hourly rainfall from daily rainfall at seven rain gauge stations in Peninsular Malaysia. The cascade models are evaluated based on the capability to simulate data that preserve three important properties of observed rainfall: rainfall variability, intermittency and extreme events. The results show that both cascade models are able to simulate reasonably well the commonly used statistical measures for rainfall variability (e.g. mean and standard deviation) of hourly rainfall. With respect to rainfall intermittency, even though both models are underestimated, the observed dry proportion, log-Normal  model is likely to simulate number of dry spells better than log-Poisson model. In terms of rainfall extremes, it is demonstrated that log-Poisson and log-Normal  models gave a satisfactory performance for most of the studied stations herein, except for Dungun and Kuala Krai stations, which both located in the east part of Peninsula.


2012 ◽  
Author(s):  
Fadhilah Y. ◽  
Zalina Md. ◽  
Nguyen V–T–V. ◽  
Suhaila S. ◽  
Zulkifli Y.

Dalam mengenal pasti model yang terbaik untuk mewakili taburan jumlah hujan bagi data selang masa satu jam di 12 stesen di Wilayah Persekutuan empat taburan digunakan iaitu Taburan Eksponen, Gamma, Weibull dan Gabungan Eksponen. Parameter–parameter dianggar menggunakan kaedah kebolehjadian maksimum. Model yang terbaik dipilih berdasarkan nilai minimum yang diperolehi daripada ujian–ujian kebagusan penyuaian yang digunakan dalam kajian ini. Ujian ini dipertahankan lagi dengan plot kebarangkalian dilampaui. Taburan Gabungan Eksponen di dapati paling baik untuk mewakili taburan jumlah hujan dalam selang masa satu jam. Daripada anggaran parameter bagi taburan Gabungan Eksponen ini, boleh diterjemah bahawa jumlah hujan tertinggi yang direkodkan diperolehi daripada hujan yang dikategorikan sebagai hujan lebat, walaupun hujan renyai–renyai berlaku lebih kerap. Kata kunci: Jumlah hujan dalam selang masa sejam, ujian kebagusan penyuaian, kebolehjadian maksimum In determining the best–fit model for the hourly rainfall amounts for the twelve stations in the Wilayah Persekutuan, four distributions namely, the Exponential, Gamma, Weibull and Mixed–Exponential were used. Parameters for each distribution were estimated using the maximum likelihood method. The best–fit model was chosen based upon the minimum error produced by the goodness–offit tests used in this study. The tests were justified further by the exceedance probability plot. The Mixed–Exponential was found to be the most appropriate distribution in describing the hourly rainfall amounts. From the parameter estimates for the Mixed–Exponential distribution, it could be implied that most of the hourly rainfall amount recorded were received from the heavy rainfall even though there was a high occurrences of light rainfall. Key words: Hourly rainfall amount, goodness-of-fit test, exceedance probability, maximum likelihood


2016 ◽  
Vol 10 (5) ◽  
pp. 77 ◽  
Author(s):  
Zalina Mohd Daud ◽  
Siti Musliha Mat Rasid ◽  
Norzaida Abas

<p>High resolution rainfall data is an important input for studies on hydrological systems. Often time synthetic data has to be generated in the absence of historical data.  The stochastic Neyman Scott Rectangular Pulse (NSRP) model has been developed to produce synthetic data of high resolution. To take into account the spatial and temporal nature of rainfall, an extra domain had been added producing the Spatial Temporal NSRP (ST-NSRP) model. However these models require the repeat estimation of all the model parameters for each location. This work develops the Regionalized ST-NSRP model which produces a single model for a region. The regionalization approach is carried out by the station-year method using ten years (2001-2010) of records from sixteen stations. The statistical characteristics used in estimating the model’s parameters are hourly and daily coefficient of variation, hourly and daily lag one auto correlation, hourly cross correlation between sites and hourly skewness of rainfall series. Simulation process is conducted at two independent stations with distinctly different rainfall profile. The Regionalized ST-NSRP model produced statistical characteristics of the ungauged sites which matched those of the observed series fairly well even though there is a tendency to underestimate hourly covariance and lag one autocorrelation. Considering that spatial variability of rainfall is high within the studied region, this model is sufficiently robust in producing reasonable synthetic hourly rainfall series which captured the characteristics and pattern of observed rainfall series. Therefore the model has shown a promising potential in generating high resolution synthetic data for any ungauged site within the region.</p>


2012 ◽  
Vol 9 (12) ◽  
pp. 14205-14230
Author(s):  
A. Chebbi ◽  
Z. K. Bargaoui ◽  
M. da Conceição Cunha

Abstract. Based on rainfall intensity-duration-frequency (IDF) curves, a robust optimization approach is proposed to identify the best locations to install new rain gauges. The advantage of robust optimization is that the resulting design solutions yield networks which behave acceptably under hydrological variability. Robust optimisation can overcome the problem of selecting representative rainfall events when building the optimization process. This paper reports an original approach based on Montana IDF model parameters. The latter are assumed to be geostatistical variables and their spatial interdependence is taken into account through the adoption of cross-variograms in the kriging process. The problem of optimally locating a fixed number of new monitoring stations based on an existing rain gauge network is addressed. The objective function is based on the mean spatial kriging variance and rainfall variogram structure using a variance-reduction method. Hydrological variability was taken into account by considering and implementing several return periods to define the robust objective function. Variance minimization is performed using a simulated annealing algorithm. In addition, knowledge of the time horizon is needed for the computation of the robust objective function. A short and a long term horizon were studied, and optimal networks are identified for each. The method developed is applied to north Tunisia (area = 21 000 km2). Data inputs for the variogram analysis were IDF curves provided by the hydrological bureau and available for 14 tipping bucket type rain gauges. The recording period was from 1962 to 2001, depending on the station. The study concerns an imaginary network augmentation based on the network configuration in 1973, which is a very significant year in Tunisia because there was an exceptional regional flood event in March 1973. This network consisted of 13 stations and did not meet World Meteorological Organization (WMO) recommendations for the minimum spatial density. So, it is proposed to virtually augment it by 25, 50, 100 and 160% which is the rate that would meet WMO requirements. Results suggest that for a given augmentation robust networks remain stable overall for the two time horizons.


2016 ◽  
Author(s):  
David N. Dralle ◽  
Nathaniel J. Karst ◽  
Kyriakos Charalampous ◽  
Sally E. Thompson

Abstract. The study of single streamflow recession events is receiving increasing attention following the presentation of novel theoretical explanations for the emergence of power-law forms of the recession relationship, and drivers of its variability. Individually characterizing streamflow recessions often involves describing the similarities and differences between model parameters fitted to each recession time series. Significant methodological sensitivity has been identified in the fitting and parameterization of models that describe populations of many recessions, but the dependence of estimated model parameters on methodological choices has not been evaluated for event-by-event forms of analysis. Here, we use daily streamflow data from 16 catchments in northern California and southern Oregon to investigate how combinations of commonly used streamflow recession definitions and fitting techniques impact parameter estimates of a widely-used power-law recession model. We show that: (i) methodological decisions, including ones that have received little attention in the literature, can impact parameter value estimates and model goodness-of-fit; (ii) the central tendencies of event-scale recession parameter probability distributions are largely robust to methodological choices, in the sense that differing methods rank catchments similarly according to the medians of these distributions; (iii) recession parameter distributions are method-dependent, but roughly catchment-independent, such that changing the choices made about a particular method affects a given parameter in similar ways across most catchments; and (iv) the observed correlative relationship between the power law recession scale parameter and catchment antecedent wetness varies depending on recession definition and fitting choices.


1985 ◽  
Vol 248 (3) ◽  
pp. R378-R386 ◽  
Author(s):  
M. H. Nathanson ◽  
G. M. Saidel

Optimal experimental design is used to predict the experimental conditions that will allow the "best" estimates of model parameters. A variety of criteria must be considered before an optimal design is chosen. Maximizing the determinant of the information matrix (D optimality), which tends to produce the most precise simultaneous estimates of all parameters, is commonly considered as the primary criterion. To complement this criterion, we present another whose effect is to reduce the interaction among the parameter estimates so that changes in any one parameter can be more distinct. This new criterion consists of maximizing the determinant of an appropriately scaled information matrix (M optimality). These criteria are applied jointly in a multiple-objective function. To illustrate the use of these concepts, we develop an optimal experimental design of blood sampling schedules using a detailed ferrokinetic model.


2018 ◽  
Vol 140 (7) ◽  
Author(s):  
Ahmed Ramadan ◽  
Connor Boss ◽  
Jongeun Choi ◽  
N. Peter Reeves ◽  
Jacek Cholewicki ◽  
...  

Estimating many parameters of biomechanical systems with limited data may achieve good fit but may also increase 95% confidence intervals in parameter estimates. This results in poor identifiability in the estimation problem. Therefore, we propose a novel method to select sensitive biomechanical model parameters that should be estimated, while fixing the remaining parameters to values obtained from preliminary estimation. Our method relies on identifying the parameters to which the measurement output is most sensitive. The proposed method is based on the Fisher information matrix (FIM). It was compared against the nonlinear least absolute shrinkage and selection operator (LASSO) method to guide modelers on the pros and cons of our FIM method. We present an application identifying a biomechanical parametric model of a head position-tracking task for ten human subjects. Using measured data, our method (1) reduced model complexity by only requiring five out of twelve parameters to be estimated, (2) significantly reduced parameter 95% confidence intervals by up to 89% of the original confidence interval, (3) maintained goodness of fit measured by variance accounted for (VAF) at 82%, (4) reduced computation time, where our FIM method was 164 times faster than the LASSO method, and (5) selected similar sensitive parameters to the LASSO method, where three out of five selected sensitive parameters were shared by FIM and LASSO methods.


2013 ◽  
Vol 17 (10) ◽  
pp. 4259-4268 ◽  
Author(s):  
A. Chebbi ◽  
Z. K. Bargaoui ◽  
M. da Conceição Cunha

Abstract. Based on rainfall intensity-duration-frequency (IDF) curves, fitted in several locations of a given area, a robust optimization approach is proposed to identify the best locations to install new rain gauges. The advantage of robust optimization is that the resulting design solutions yield networks which behave acceptably under hydrological variability. Robust optimization can overcome the problem of selecting representative rainfall events when building the optimization process. This paper reports an original approach based on Montana IDF model parameters. The latter are assumed to be geostatistical variables, and their spatial interdependence is taken into account through the adoption of cross-variograms in the kriging process. The problem of optimally locating a fixed number of new monitoring stations based on an existing rain gauge network is addressed. The objective function is based on the mean spatial kriging variance and rainfall variogram structure using a variance-reduction method. Hydrological variability was taken into account by considering and implementing several return periods to define the robust objective function. Variance minimization is performed using a simulated annealing algorithm. In addition, knowledge of the time horizon is needed for the computation of the robust objective function. A short- and a long-term horizon were studied, and optimal networks are identified for each. The method developed is applied to north Tunisia (area = 21 000 km2). Data inputs for the variogram analysis were IDF curves provided by the hydrological bureau and available for 14 tipping bucket type rain gauges. The recording period was from 1962 to 2001, depending on the station. The study concerns an imaginary network augmentation based on the network configuration in 1973, which is a very significant year in Tunisia because there was an exceptional regional flood event in March 1973. This network consisted of 13 stations and did not meet World Meteorological Organization (WMO) recommendations for the minimum spatial density. Therefore, it is proposed to augment it by 25, 50, 100 and 160% virtually, which is the rate that would meet WMO requirements. Results suggest that for a given augmentation robust networks remain stable overall for the two time horizons.


2017 ◽  
Vol 65 (4) ◽  
pp. 479-488 ◽  
Author(s):  
A. Boboń ◽  
A. Nocoń ◽  
S. Paszek ◽  
P. Pruski

AbstractThe paper presents a method for determining electromagnetic parameters of different synchronous generator models based on dynamic waveforms measured at power rejection. Such a test can be performed safely under normal operating conditions of a generator working in a power plant. A generator model was investigated, expressed by reactances and time constants of steady, transient, and subtransient state in the d and q axes, as well as the circuit models (type (3,3) and (2,2)) expressed by resistances and inductances of stator, excitation, and equivalent rotor damping circuits windings. All these models approximately take into account the influence of magnetic core saturation. The least squares method was used for parameter estimation. There was minimized the objective function defined as the mean square error between the measured waveforms and the waveforms calculated based on the mathematical models. A method of determining the initial values of those state variables which also depend on the searched parameters is presented. To minimize the objective function, a gradient optimization algorithm finding local minima for a selected starting point was used. To get closer to the global minimum, calculations were repeated many times, taking into account the inequality constraints for the searched parameters. The paper presents the parameter estimation results and a comparison of the waveforms measured and calculated based on the final parameters for 200 MW and 50 MW turbogenerators.


2021 ◽  
Vol 13 (14) ◽  
pp. 7911
Author(s):  
Ibrahim Alsaidan ◽  
Mohamed A. M. Shaheen ◽  
Hany M. Hasanien ◽  
Muhannad Alaraj ◽  
Abrar S. Alnafisah

For the precise simulation performance, the accuracy of fuel cell modeling is important. Therefore, this paper presents a developed optimization method called Chaos Game Optimization Algorithm (CGO). The developed method provides the ability to accurately model the proton exchange membrane fuel cell (PEMFC). The accuracy of the model is tested by comparing the simulation results with the practical measurements of several standard PEMFCs such as Ballard Mark V, AVISTA SR-12.5 kW, and 6 kW of the Nedstack PS6 stacks. The complexity of the studied problem stems from the nonlinearity of the PEMFC polarization curve that leads to a nonlinear optimization problem, which must be solved to determine the seven PEMFC design variables. The objective function is formulated mathematically as the total error squared between the laboratory measured terminal voltage of PEMFC and the estimated terminal voltage yields from the simulation results using the developed model. The CGO is used to find the best way to fulfill the preset requirements of the objective function. The results of the simulation are tested under different temperature and pressure conditions. Moreover, the results of the proposed CGO simulations are compared with alternative optimization methods showing higher accuracy.


Sign in / Sign up

Export Citation Format

Share Document