WEIBULL COMPARISON BASED ON RELIABILITY, AVAILABILITY, MAINTAINABILITY, AND DEPENDABILITY (RAMD) ANALYSIS
Anas Sani Maihulla
Department of Mathematics, Sokoto State University, Sokoto- Nigeria [email protected]
Ibrahim Yusuf
Department of Mathematical sciences, Bayero University, Kano- Nigeria [email protected]
Saminu I. Bala
Department of Mathematical sciences, Bayero University, Kano- Nigeria
sibala. [email protected] .ng
Abstract
As a continuous probability distribution, the Weibull distribution is widely used in the study of reliability, availability and other life data. In this research, we propose the RAMD analysis to estimate the three-parameter Weibull distribution. The estimation of the distribution parameters is an important problem that has received a lot of attention from researchers because of their effects in several measurements. The real data results indicate that our proposed estimation method is significantly consistent in estimation compared to the RAMD analysis method. The numerical values of filtration system reliability and availability were calculated using Maple software. The system of first-order differential equations is formulated using a mnemonic approach and solved recursively. Several scenarios were examined to determine the impact of the models under consideration. The calculations were done with Maple 13 software. Other reliability measures such as mean time to failure (MTTF), mean time to repair (MTTR), and dependability ratio was estimated. The comparative analysis was conducted using a reverse osmosis (RO) filtration system.
Keywords: Availability, life data, Weibull-distribution, Maple, Reliability, Reverse Osmosis
1. Introduction
A reliability study's main goal is to provide information that can be used to make decisions. Performance Analysis of K out of N Reverse Osmosis (RO) system in the treatment of wastewater using RAMD analysis was carried out by Ibrahim et al. [19]. Based on the analyses, Jun et al.[9] Concluded The key to the life model is the estimating parameters; it can accurately predict the life of a product in terms of reliability. Because obtaining the estimated parameters is difficult, the process of estimating 3-parameter Weibull distribution is important. Iterative methods are required for the estimation process using the maximum likelihood function, which takes a long time and effort. RAMD analysis is the proposed method for the estimation of the 3- parameter weibull distribution. In comparison to the maximum likelihood method, the real data results show that our proposed estimation method is significantly more consistent in estimation. Maihulla A. S. et al. [10] investigated the Reliability, Availability, Maintainability, and Dependability Analysis of a Complex Reverse Osmosis Machine System in Water Purification. The swedish physicist Weibull proposed the Weibull distribution (Walooddi weibull, 1939). He used it to determine the tensile strength of various materials. Since then, it's been widely used in reliability and life testing problems, such as determining a component's time to failure or life length, which is measured from a specified time until it fails as stated by [14].
Jukic and Markovic, [8] and Walpole et al., [18] studied the Reliability Estimation of Three Parameters Weibull Distribution based on Particle Swarm Optimization. Basheer et al. [5] demonstrated how to treat the pollutants found in olive mill wastewater (total organic carbon (TOC), dissolved organic carbon (DOC), total phosphorus (TP), total nitrogen (TN), and total polyphenols), a sequential Direct Contact Membrane Distillation (DCMD) and a Reverse Osmosis (RO) hybrid membrane system were used. Similar study was also conducted by Teresa [3] but with recoveries of some economical merits. The influence of permeate flow and pressure on pollutant parameter removals was also included in his study. One of the biggest challenges that humanity must address in the twenty-first century is the scarcity of freshwater. That is why Biniaz et al. [17] investigated potential development of an environmentally friendly, economical, and energy-efficient membrane distillation process. Industrial and household wastes can lessen pollution than to the membrane distillation (MD) method. The Gumbel-Hougaard family copula was used to model the reliability and performance of a solar photovoltaic system was analyzed by Maihulla A. S. et al. [12]. Failure of the purification filter is a serious issue that the water purification industries are facing. Modern technical challenges include removing human-caused contaminants from drinking water. After a brief discussion of the treatment stages that municipal water through before it reaches your tap, you'll learn about the detectable contamination of drinking water caused by anthropogenic (humanmade) toxins that is still present. These were all examined by Maihulla et al. [11]. The issue surrounding the comparison of the RAMD analysis and the 3-parameter Weibull was not covered in any of the aforementioned publications.
Our motivation for using a reverse osmosis filtration system as a test bed for the two methods was to see how well they worked together (3-parameter Weibull and RAMD analysis) stems from a serious problem that the water purification industry is facing due to purification filter failure. As a result, there has been a slow advancement in technological advancement in water purification, as well as its importance in the lives of people all over the world. The filtration industry is working hard to keep up with the increasing complexity of the systems. RAMD analysis was used to examine the filtration system's strength, efficiency, and performance improvement, according to the findings of the paper. If the strength, efficiency, and performance of the filtration system are evaluated, users will be able to save money on medical care due to contaminated water. Protect yourself from pollution in the water. The study is divided into five parts, one of which is the current
Anas sani Maihulla, Ibrahim Yusuf & Saminu Bala
WEIBULL COMPARISON BASED ON RELIABILITY, AVAILABILITY, RT&A, No 1 (72)
MAINTAINABILITY, AND DEPENDABILITY (RAMD)_Volume 18, March 2023
introduction. Modeling of filtration is discussed in the section. Materials and methods are
Covered in Section 2. Section 3 discussed the proposed RAMD analysis method, which was compared to the 3-parameter Weibull, as well as an analytical analysis of the system. Section 4 presented the result. Section 5 included a discussion and explanation of the findings, as well as a conclusion.
1.1. Filtration
One of the most common methods for removing these materials is gravity filtration. This procedure involves passing water containing solid impurities (for example, precipitates after water softening) through a porous material, usually sand and gravel layers. The force of gravity pushes the water through the medium. The gaps between the sand and gravel grains allow small water molecules to flow through. Precipitation-derived solids, on the other hand, become trapped in the pores and thus remain in the porous medium. The solid contaminants have been removed from the water that passes through the bottom of the filter.
Figure 1: RO Filtration System's Block diagram
1.2. RAMD indices for subsystem 1 (sand filter) Sand Filter: George et al. [6] stated that, Sand filters are widely used in water purification and work on a completely different mechanism to remove suspended matter. Rather than passing through small orifices through which particles cannot pass, the water passes through a 750 mm deep bed of filter medium, typically 0.75 mm sand.
1.3. RAMD indices for subsystem 2. (Precision filter) The cylinder shell is usually made of stainless steel, and the internal filter elements are made of PP melt-blown, wire burning, folding, titanium filter, activated carbon filter, and other tubular filter elements, which are chosen based on the different filter media and design process to meet the effluent quality requirements.
1.4. RAMD indices for subsystem 3 (activated carbon filter)
Activated carbon filters are used to purify water without leaving any harmful chemicals behind according to Y. K. Siong et al. [21]. For water treatment, a prototype is being created using activated
carbon and an ultraviolet radiation system. Analysis of surface area and porosity. For comparison of surface morphology, scanning electron microscopy (SEM) was used to obtain magnified images of GAC-A and GAC-B.
2. Methods
With the probability density function, let T be a random variable that represents time to failure. f(t), where f(t) is the 3-parameter Weibull distribution's probability density function, which can be written in (1) below, as analyzed by: Jukic and Markovic, [8] and Jun et al., [9].
fit; a, p, Y) = § {^f)'1 exp [- t> 0 (1)
0 t< 0 Where a > 0 denotes the shape parameter, [3 > 0 denotes the scale parameter, and Y< t denotes the location parameter.
Integrating the probability density function yields the cumulative distribution function for the 3-parameter Weibull distribution in (2).
Fit) = P(T <t)= fg fit) a, P, Y)dt = 1 - e^V) (2)
The three-parameter Weibull distribution's mean and variance are (3) and (4) below, as in (Muraleedharan, [14])
A=/-+JbT(^±) (3)
Where /¿= E(t) = MTTF (Mean time to failure)
(4)
2.1. List of notations and definitions
Represent the system in working state
Pi P2
ui i=i, 2,3
P) j = 1,2,3
PXit)
Represent the system in failed state
Represent the initial state of the system working in full capacity state.
Represent the state in which one parallel unit is failed Represent the state in which two parallel unit is failed Represent the state in which three parallel unit is failed Represent the failure rates subsystems Represent the repair rates subsystems Probability to remain at xth state at time t
—Pxit), x= 0,1,2,3. Represent the derivative with respect to time t
2.2. Reliability
The chance that a device will run without failure for a particular period of time is referred to as reliability under the operational conditions indicated.
R(t) = ¡^ f(x)dx Reliability function (5)
2.3. MTBF
Mean Time between Failures (MTBF): The mean time between failures refers to the average period of good system functioning. The MTBF is the reciprocal of the constant failure rate or the ratio of the test time to the number of failures when the failure rate is reasonably consistent over the operating period. It is given by (6), as in Ashish K. and Monika S. [13].
MTBF = J0 R(t)dt = J0 e~et = & Mean Time between Failure (6)
2.4. MTTR
Mean Time between repairs (MTTR): is the reciprocal of the system repair rate, as in (7) below: MTTR = & Mean Time to Repair (7)
2.5. Availability
Availability: For repairable systems, availability is a performance criterion that considers both the system's dependability and maintainability. It's defined as the probability that the system will work correctly when it's needed. JG Wohl [20]. That is the ratio of life time to total time, as in (8) below.
. , .,., Life time Life time MTTF
Availability = —- = ---= --(8)
J Total time Life time+Repair time MTTF+MTTR v '
2.6. Maintainability
Maintainability: When maintenance is necessary, is a design, installation, and operation feature that is often stated as the possibility that a machine can be kept in, or returned to, a particular operational condition within a specified time frame, according to Wohl, JG [20]. It is given by (9) below.
M(t) = 1 - e^MTTR^ Maintainability function (9)
2.7. Dependability A. Ebeling [4] defined dependability as a measure of a system's availability, reliability, maintainability. The advantage of dependability is that it allows for cost, reliability, and maintainability comparisons. For random variables with exponential distribution, the dependability ratio is as follows:
P = Repair rate , a = Failure rate
^ _ a _ MTBF F MTTR
The importance of maintenance is reflected in the high value of the dependability ratio. The dependability value increases if availability is greater than 0.9 and decreases if availability is less than 0.1, according to C. Li, S. Besarati, and colleagues [2]. The following formula in (10) calculates the minimum value of dependability:
& Ind d In d
Dmin _ 1 - (-¡T^'^ - e~~) (10)
2.8. Functions of reliability and failure rate
R(t) Is the reliability function (also known as the survival function). The probability that the device will operate without failure for a given period of time under the specified operating conditions is known as reliability. This included in the study carried out by Haldar and Mahadevan, [7] and
Walpole et al., [18]. The R(t) in (11) below is in terms of the failure and repair rates of the three parameter weibull distribution, unlike that of (5).
R(t)_P(T<t)_ f0f(f,a,p,Y)dt_l-F(t)_e ( $ ( (11)
h(t) Represents the failure rate (also known as the hazard rate). The failure rate is expressed in terms of failures per unit time. It is computed as the ratio of number of failures of the items undergoing the test time. From T _ t to T _ t + At, given that it survived to time t, this was investigated by Nachlas, [15] and Walpole et al.,[18].
P(t<T<t + AtlT>t)_ p(t<T<t+At) _ F(t*At)-F(t) _ (12)
v 1 J P(T>t) R(#) R(t) v '
We get the failure rate by dividing this ratio by At and finding the limit as At^0: in the three parameter weibull distribution, as below in (13). The computation of availability, MTTF, MTTR, MTBF as well as dependability analyses for RAMD against three parameter weibull distribution were presented in table 2.
_ lim P(t<T<t + M\T>t) _Um Arm _ m (13)
At^O At At->0 AR(t) R(t) v '
And for each correspondent subsystem: by substituting h(t) as failure rate, we got (14) below:
And
1 1
MTTF _ -_ — (Mean time between failures) (14)
Failure rate h(t) v ' v '
MTTR _ -1- (Mean time to repair) (15)
Repair rate 1
And therefore the Availability of the system can be calculated using the relation below as in Monika et al. [16]. Availability, also using three parameter weibull distribution can be computed as (16) below:
Availability _ —MTBF— (16)
J MTBF+MTTR v '
3. The proposed RAMD Analysis method
Figure 2: Transition diagram for subsystem 1
3 P0(t) = -3a&P0 + ftP& 3P&(t) = -(2a& + ft) P& + 3a&P0 + ftP+
P2(t) = -(a& + ft) P2 + 2a&P& + ftP, 3P2(t) = - ftP+ +a&P&
(17)
(18)
(19)
(20)
The steady-state probabilities of the system are obtained by imposing the following restrictions: — ^ 0, as t ^ œ. see Arora and Kumar [1].
Under steady state, equation (17) - (20) reduces to
P = — P H& PiFo
(21)
Substituting (21) into (18)
p = p
r2 - P2 r0
(22)
Substituting (22) into (19)
P = P
Using normalization condition
Po + P& +P+ + P,= 1 Substituting (21) and (22) and (23) into (24) we have:
(23)
(24)
P + = — P + P + P = 1 *0 + Pi 10 + "2 ^0 + Pi ^0 = L
Po =
P0
P0+6aJ+3P2ai+6"ja/
Table 2 contains important device output metrics that have been extracted
(25)
Figure 3: Transition diagram for subsystem 2
3 P0(t) = -2a4P0 + P4P& 3P&(t) = 2a4P0 - P4P&
Under steady state, equation (20) and (21) reduces to:
-2aAP0 + piP1 = 2aAP0 - piP1
And
p = 2W1 p H& = 17 F<0
Using normalization condition
p0 + p&= 1
Substituting (22) into (23) we have:
p + +W± p = 1
^0 + "4 ^0 = 1
P =
Pi
"i+2a1
Table 2 contains important device output metrics that have been extracted.
Figure 4: Transition diagram for subsystem 3
3 P0(t) = -a,P0 + P3P1 3 P&(t) = -(P3 + a,)Pi + a,P0 + P3P2
3-t P2(t)= -P3P2 + a,P&
Under steady state, equation (25)- (27) yield p = P
Substituting (28) into (26) under steady state we have:
(26) (27)
(28)
(29)
(30)
(31)
(32)
(33)
P - -0 P
- "2 F)
(35)
Using normalization condition
P) + P& + P+ -l It follows that:
o2
P) -
P2+«2+p3«3
(36)
(37)
3.1. System Reliability
Rsys(t) = Rsi(t) x Rs2(t) x Rs3(t)
Rsys(t) = e-x e-x e-Rsys(t)= e-(a,+a/+a0)t
3.2. System Availability
Arranged in series, failure of one cause the complete failure of the system.
(38)
(39)
(40)
-^sys — Asl X As2 x
"0
A — (■
sys ("'0+6a'0+3"2a1+6 "ta/' v "4+20.4
2)x( -
"4
■)X(.
^ — , ).)\ - X ,).)_- X , )
sys \).)8))5j V).ll- V).l
F0 + a0 +p3a3
)
)l2l
.)l5B5j
¿sys - 0.87445 X 0.81818 X 0.76341 4sys - 0.18249
(41)
(42)
(43)
(44)
(45)
Table 2 contains important device output metrics that have been extracted.
4. Results
The numerical outcomes from RAMD and Weibull analysis performed using Maple software are presented in this section. The discussion of the comparative analysis in this section is found in the part that follows. The tables' contents are translated into graphs. Table 1 below shows the failure and repair rates for the corresponding subsystems. The comparison findings for the two aforementioned methodologies were shown in Tables 2, 3, and 4.
Table 1: Failure and repair rates for subsystems
Subsystem Failure Rate (a) Repair Rate (¡3)
s& a&= 0.004 P& - 0 . 54
s+ a2= 0.006 P+ - 0.67
s3 a3= 0.008 P, - 0. 90
Table 2: RAMD/Weibull indices for the RO filtration system
RAMD indices of Subsystems Subsystem Si Subsystem Subsystem System
Reliability g-0.009 £ g-0.005£ g-0.014 £ g-0.028 £
Availability (Weibull) 0.96755 0.77853 0.89706 0.67573
Availability (RAMD) 0.84556 0.69453 0.77453 0.45486
MTBF (Weibull) 104.34 198.92 88.06 391.32
MTBF(RAMD) 83.33 21.01 166.67 32.25 62.50 25.56 312.50 78.82
MTTR (Weibull) 4.77 3.21 1.96 9.94
MTTR(RAMD) 1.85 1.49 1.11 4.45
Dependability (Weibull) 0.99302 0.99746 0.99753 0.98805
Dependability (RAMD) 0.97497 0.94301 0.92891 0.85405
Dependability ratio (Weibull) 21.87 61.97 44.93
Dependability ratio (RAMD) 45.04 111.86 56.31
Table 3: Availability of a system
Subsystems 3-Parameter weibull distribution RAMD (Exponential distribution)
1 0.96755 0.84556
2 0.77853 0.69453
3 0.89706 0.77453
Availability analysis
1
s-3 ™ 0,5
(D
< RAMD
0
0 Weibull
50
100
Time
■ Weibull 1 RAMD
Figure 5: Variation of Availability with time
Table 4: Reliability of a system
Time(in days) 3-Parameter weibull distribution RAMD (Exponential distribution)
0 1.00000 0.99999
10 0.79918 0.79384
20 0.71540 0.70030
30 0.66808 0.64285
40 0.63288 0.59827
50 0.60240 0.55945
60 0.57436 0.52410
70 0.54801 0.49140
80 0.52308 0.46100
90 0.49942 0.43268
100 0.47695 0.40628
Weibull vs RAMD Plot
1,2 1
>. 0,8
| 0,6
(u
oe 0,4 0,2 0
0 10 20 30 40 50 60 70 80 90 100
Time
—•—Weibull —•— RAMD
Figure 6: Variation of Reliability with time
5. Discussion
Figures 2, 3, and 4 show the transition diagrams for the three reverse osmosis filtering units/subsystems, namely the sand filters, precision filters, and carbonated filters. A comparison of the efficiency of two methods (RAMD analysis and 3-parameter Weibull) was conducted. Table 2 lists all of the additional RAMD-Weibull metrics. It demonstrated how effective the Weibull is when compared to the RAMD. The variation in the system's availability is 0.22087, which is due to the fact that subsystem 3 has a larger variation than subsystems 1 and 2. When the MTBF of the entire system is considered, it is found that the weibull is more efficient in terms of the average operating time from the time a failed device is restored to the time it becomes failed again than RAMD analysis, with a variation of nearly 78.82. When compared to 1 and 3, the same subsystem 2, contributed significantly. Despite the fact that it had fewer units than 1 and 3, it was still a good choice. Using the RAMD analysis, the system's reliability after 50 days of operation is only 0.55945. The reliability of the same system using a 3-parameter Weibull distribution after the same period is 0.60240.
Anas sani Maihulla, Ibrahim Yusuf & Saminu Bala
WEIBULL COMPARISON BASED ON RELIABILITY, AVAILABILITY, RT&A, No 1 (72)
MAINTAINABILITY, AND DEPENDABILITY (RAMD)_Volume 18, March 2023
Figure 5 depicted the availability of the two methods as a function of time (3-parameter Weibull and RAMD analysis). The efficiency of 3-parameter Weibull over RAMD analysis is clearly observed. According to the numerical analysis in Table 4 and the corresponding values in Figure 6.
6. Conclusion
For desalinating saltwater to produce drinking water, reverse osmosis (RO) filtration system is an important technology. The performance of a desalination system is determined by the failure behavior of its components. Because the RO filtration system was designed to be low-power, the subsystems' dependability must be maintained through proper design and material selection in order for the plant to operate continuously. Furthermore, the estimation of the system's strength in terms of reliability characteristics using the 3-parameter Weibull distribution is more efficient than the RAMD analysis method. As a result, the study investigated the system's availability and reliability, as well as other features such as MTBF, MTTR, and dependability analysis. To address the issue of precision filter redundancy, more research can be done.
7. Acknowledgement
We sincerely appreciate the distinguished Sokoto State University. This is because I was given
a study leave and other benefits to allow me to finish my PhD in mathematics. This research
fits inside that.
References
[1] Arora N, Kumar D (1997) Availability analysis of steam ans power generation system in thermal power plant. Macroelectron Reliab 37(5):795-799
[2] C. Li, S. Besarati, Y. Goswami, E. Stefanakos, H. Chen, (2013) "Reverse osmosis desalination driven by low temperature supercritical organic Rankine cycle". Appl. Energy, Vol. 102 Pp 1071-1080.
[3] D. Teresa Sponza (2021). "Treatment of Olive Mill Effluent with Sequential Direct Contact Membrane Distillation (DCMD)/Reverse Osmosis (RO) Hybrid Process and Recoveries of Some Economical Merits" Res Appl. 2021 vol. 1(1): 001-010. Advanced Journal of Physics Research and Applications. Vol. 1(1): 001-010.
[4] Ebeling A (2000) An introduction to reliability and maintainability engineering. Tata Mcgraw Hill Company Ltd, New Delhi
[5] G. T. Basheer, Zakariya and Y Algamal. (2021). Reliability Estimation of Three Parameters Weibull Distribution based on Particle Swarm Optimization. Pakistan journal of statistics and operational research. DOI: http://dx.doi.org/10.18187/pjsor.v17i1.2354 Pak.j.stat.oper.res. Vol.17 No. 1 2021 pp 35-42.
[6] George Solt CEng, FIChemE, (2002) in Plant Engineer's Reference Book (Second Edition).
[7] Haldar A., Mahadevan S., (2000). Probability, Reliability and Statistical Methods in Engineering Design. John Wiley and Sons, Inc.
[8] Jukic D., Markovic D, (2010). On nonlinear weighted total least squares parameter estimation problem for the three-parameter weibull density. Applied Mathematical Modelling, 34, 1839-1848.
Anas sani Maihulla, Ibrahim Yusuf & Saminu Bala
WEIBULL COMPARISON BASED ON RELIABILITY, AVAILABILITY, RT&A, No 1 (72)
MAINTAINABILITY, AND DEPENDABILITY (RAMD)_Volume 18, March 2023
[9] Jun W.C, Xin L, Xin C., (2017). Parameter Evaluation of 3-parameter weibull Distribution based on Adaptive Genetic Algorithm. 2nd International Conference on Machinery, Electronics and Control Simulation, Advances in Engineering Research, volume 138.
[10] Maihulla, A. S., Yusuf, I., and Bala, S. I. (2021). Performance Evaluation of a Complex Reverse Osmosis Machine System in Water Purification using Reliability, Availability, Maintainability and Dependability Analysis. Reliability: Theory & Applications, 16(3), 115-131.
[11] Maihulla, A. S., & Yusuf, I. (2022). RELIABILITY ANALYSIS OF REVERSE OSMOSIS FILTRATION SYSTEM USING COPULA. Reliability: Theory & Applications, 17(2), 163-177. Retrieved from http://www.gnedenko.net/RTA/index.php/rta/article/view/890.
[12] Maihulla, A.S., Yusuf, I. and Salihu Isa, M. (2022), "Reliability modeling and performance evaluation of solar photovoltaic system using Gumbel-Hougaard family copula", International Journal of Quality & Reliability Management, Vol. 39 No. 8, pp. 20412057. https://doi.org/10.1108/IJQRM-03-2021-0071
[13] Monika S. and Ashish K. (2019). "Performance analysis of evaporation system in sugar industry using RAMD analysis" journal of Brazilian society of mechanical science and engineering Vol. 41. (175). https://doi.org/10.1007/s40430-019-1681-3
[14] Muraleedharan G., (2013). Characteristic and Moment Generating Functions of Three Parameter Weibull Distribution-an Independent Approach. Research Journal of Mathematical and Statistical Sciences, Vol.1(8), 25-27.
[15] Nachlas J.A., (2005). Reliability Engineering Probabilistic Models and Maintenance Methods. Taylor & Francis Group.
[16] Nivedita Gupta, Monika Saini and Ashish Kumar (2020). Behavioral Analysis of Cooling Tower in Steam Turbine Power Plant using Reliability, Availability, Maintainability and Dependability Investigation. Journal of Engineering Science and Technology Review 13 (2) (2020) 191 - 198. doi:10.25103/jestr.132.23
[17] P. Biniaz, N. T. Ardekani, M. A. Makarem and M. R. Rahimpour (2019). "Water and Wastewater Treatment Systems by Novel Integrated Membrane Distillation (MD)" ChemEngineering, vol. 3, no. 8
[18] Walpole R.E., Myers R. H., Myers S.L. and Ye K., (2012). Probability and statistics for Engineers and Scientists. Pearson Education, Inc.
[19] Walpole R.E., Myers R. H., Myers S.L. and Ye K., (2012). Probability and statistics for Engineers and Scientists. Pearson Education, Inc.
[20] Wohl JG (1996) system operational readiness and equipment dependability. IEEE Trans Reliability 15 (1): 1-6.
[21] Y. K. Siong , J. Idris a , M. Mazar Atabaki , (2013) "Performance of Activated Carbon in Water Filters" See discussions, stats, and author profiles for this publication at: https://www.researchgate.net/publication/234060484 Article in Water Resources.