Efficient simulation of a space-time Neyman-Scott rainfall model
|
|
- Magdalene Chapman
- 5 years ago
- Views:
Transcription
1 WATER RESOURCES RESEARCH, VOL. 42,, doi: /2006wr004986, 2006 Efficient simulation of a space-time Neyman-Scott rainfall model M. Leonard, 1 A. V. Metcalfe, 2 and M. F. Lambert 1 Received 21 February 2006; revised 6 July 2006; accepted 27 July 2006; published 9 November [1] Existing space-time Neyman-Scott models characterize rainfall as the arrival of rain cells, clustered in time and independently distributed in space. Each cell is described as a cylinder having a random intensity (height) and random radial coverage. With this formulation it is possible to have cells with centers lying outside of a target simulation region yet having radii large enough to cover points within the region. To avoid significant boundary effects, it is necessary to include these points in the simulation. However, this can introduce inefficiency into the algorithm that is computationally restrictive. To overcome this, an efficient method is derived and demonstrated to improve computational performance. Citation: Leonard, M., A. V. Metcalfe, and M. F. Lambert (2006), Efficient simulation of a space-time Neyman-Scott rainfall model, Water Resour. Res., 42,, doi: /2006wr Introduction [2] There is an extensive amount of literature on the development of clustered rainfall models at a single point [Onof et al., 2000]. The Neyman-Scott rectangular pulse model is one type of clustered process and has been demonstrated to successfully reproduce a wide range of rainfall statistics at various locations. The general structure of the model consists of a set of random variables to describe the arrival of storms in time, the displacement of storm cells relative to the storm origin, the lifetime of storm cells and the intensity of storm cells. Where cells overlap in time, the total rainfall intensity is the sum of the rainfall contributions from each cell. Cowpertwait [1995] introduced a spatial extension to the Neyman-Scott point model by specifying random variables for the position and radial coverage of each cell in space. Where cells overlap in space, as for the temporal case, cell depths are aggregated. While this spatial description of rainfall is highly idealized (for example, it does not include advection), it is (1) one of only a few rainfall models that characterizes rainfall continuously in time and space and (2) can be solely calibrated to rain gauge data. Wheater et al. [2005] discuss an alternative space-time rainfall model, having a Bartlett-Lewis temporal process and calibrated using radar data. [3] A simulation proceeds by sampling, from calibrated distributions, the arrival of storms, number of rain cells and cell properties over some target region. However, with this formulation it is possible to have rain cells with centers lying outside of a target region, yet having radii large enough to cover points within the region. This introduces a boundary effect that significantly reduces the simulated rainfall depth at points within the region if these cells are ignored. It is essential to account for this boundary effect 1 School of Civil and Environmental Engineering, University of Adelaide, Adelaide, Australia. 2 School of Mathematical Sciences, University of Adelaide, Adelaide, Australia. Copyright 2006 by the American Geophysical Union /06/2006WR because the size of a rain cell is often of a similar magnitude to the size of the target region, meaning that all points within the region (and not just those near the outer perimeter) will be affected. One approach for avoiding this is to wrap the effects of rain cells across opposing points on the boundary of the region. This approach is efficient and effective with the exception that spurious cross correlations will be observed for points separated by large distances. An alternative approach is to mitigate the boundary effect by implementing a buffer around the target region. This approach is demonstrated to substantially inhibit the computational efficiency of the model, hence an algorithm is proposed that avoids the use of a buffer region. This algorithm directly simulates the number of cells that occur outside of the target region yet are known to intersect it. Whereas the buffer algorithm is approximate due to the finite size of the buffer, the direct algorithm is exact with respect to reproducing rainfall depths within the region. The method is not specific to the Neyman-Scott model as it could be applied to similar cluster-based models having the same spatial structure. 2. Algorithm for Simulation of Cells Using a Buffer Region [4] Consider a target region, of fixed radius, r t,asin Figure 1. For any given point within this region it is desired to model the rainfall process using the space-time Neyman- Scott model. According to this model, cells associated with a particular storm arrive within the region as a Poisson process having a spatial rate 8p r t 2, where 8 is the rate parameter. [5] It is however also possible to have cells arrive outside of the target region, referred to as an outer region, yet having a cell radius, R c, greater than the distance from the circumference of the target region to the center of the cell, R xy. Consider therefore a buffer, having radius, r t + r, where r is defined relative to the target region. It is possible to also simulate cells over this region according to a Poisson process with the same spatial rate parameter, 8. Cells that intersect the target region are accepted, otherwise cells that lie solely outside are rejected. 1of5
2 LEONARD ET AL.: TECHNICAL NOTE Figure 1. Schematic diagram of cells generated inside and outside the target region. [6] The arrival of cells according to a Poisson distribution requires their cells to be uniformly distributed over the region. Efficiently simulating uniformly over a circular region can be achieved using a transformation to polar coordinates. The equations for simulating a point (x,y) uniformly within a circular region of radius r t + r are, pffiffiffiffi x ¼ ðr t þ rþ p U ffiffiffiffi cos Q ð1þ y ¼ ðr t þ rþ U sin Q where Q Uniform[0,2p] and U Uniform[0,1]. [7] To illustrate the impact of a buffer region, the Arno Basin case study given by Cowpertwait et al. [2002] is used as an example. A circular target region of radius r t =65km encompasses all of the rain gauges within this region. Cell radii are distributed exponentially with parameter f. For the two example months of January and July reported values of this parameter for the Arno Basin are f = km 1 and f = km 1 respectively. An exponential distribution of cell radii having f = km 1 is capable of producing cells with large radii, for example, the 0.99 quantile gives a cell radius of 103 km. Even though such cells have a low probability of occurrence, they need to be included in a simulation otherwise the simulated rainfall will be substantially lower than the observed. For this reason, the buffer region needs to be sufficiently large to accommodate cells that occur with large radii up to an equally large distance. [8] Consider different sizes for a buffer region, r, set to quantiles 0.0, 0.5, 0.8, 0.95, 0.99 and of the exponential distribution of cell radii. The 0.0 quantile gives the case of having no buffer. Figure 2 shows the simulated proportion of mean rainfall depth for each buffer size for the two example months, where increasing buffer sizes give increasing proportions. The nondimensional quantity f(r t + r) is used to standardize the cell size with respect to the size of the region. In order to obtain accurate estimates of the simulated proportion, the statistic was estimated from an arbitrarily long simulation of 120,000 years length (limited only because of computational requirements). The point located at the center of the 65 km target region was used for 2of5 the comparison, while points closer to the perimeter will demonstrate an even greater reduction in the simulated proportion. [9] Figure 2 gives a theoretical comparison for a point at the center of the region and shows that the proportion will reduce significantly in the event of either smaller target regions or larger expected cell radii (smaller f). The proportion, denoted m, can be obtained theoretically as the ratio of the number of cells landing within the target and buffer region that overlap the center point to the number of cells that overlap the center point for a region extending infinitely, given as m ¼ Z 2p Z rtþr Z 2p Z 1 ¼ 1 e f ð rþrt 8qe fq dq dq 8qe fq dq dq Þ ð1 þ fr þ fr t Þ ð2þ [10] Figure 2 shows that the size of the buffer region strongly affects the simulated statistics within the region, and that large buffer sizes are necessary to avoid reduction in the simulated statistics. For example, in January, the proportion of rainfall is 0.79 for the case having no buffer, and for a buffer set to the 0.99 quantile (r = 103 km), the proportion rainfall is By comparing the two months, Figure 2 also shows that the proportion of simulated rainfall varies with each month due to the variation in the parameter f. [11] The computational inefficiency introduced by the buffer is proportional to the ratio of the area of the target region to this same region with the additional buffer. Assuming, that the 0.99 quantile is used to set the buffer size, the total simulation region is approximately 6.5 times larger for January and 3 times larger for July. For the Arno Basin case study, an average month has a total simulation region 4 times larger than the target region. This inefficiency is important given that the computational requirements for Figure 2. Proportion reduction in simulated rainfall depth, simulations for January and July compared with theoretical result.
3 LEONARD ET AL.: TECHNICAL NOTE simulating the model are intensive. In order to address this problem, the following section develops an algorithm for directly simulating the number, location and radius of cells that land outside of the target region yet partially cover the region. 3. Algorithm for Direct Simulation of Cells [12] Consider a target region, of fixed radius, r t,asin Figure 1. Cells arrive within this region according to a Poisson process with spatial rate parameter, 8. Consider also the same Poisson process occurring over an outer region that extends infinitely, with the exception that it is desired to retain only those cells that intersect the target region. This occurs when the cell radius, R c, is greater than the distance from the circumference of the target region to the center of the cell, R xy. The motivation then is to derive three distributions for the number, location and radius of cells, where each distribution is conditioned on the event that R c > R xy. Equation (3) defines the discrete distribution of the number of cells in the outer region, N o, that intersect with the target region. Having defined the distribution for the number of cells, it is necessary to define the properties for each given cell. The continuous distribution of cell centers, R xy, is given in equation (4) for distances defined relative to the circumference of the target region. This distribution is conditioned on the event R c > r xy, as it is more likely for cells close to the target region to intersect it. The continuous distribution for cell radii is defined in equation (5). This distribution also depends upon the event R c > r xy, since the cell radius must be larger than the distance r xy. PN o ¼ n o jr c > R xy f Rxy f Rc r xy jr c > r xy r c jr c > r xy [13] All other properties of the cell, the intensity, duration and starting time are independent of the location of the cell within the outer region Number of Outer Cells Intersecting Target [14] The following derivation relies on the Poisson distribution as the limit of a large number of Bernoulli trials. Consider an element, as in Figure 1, having radial increment, dr, elemental angle, dq, and total radius, r + r t, with cells arriving over that element according to a Poisson process having rate, 8. The probability of one cell landing inside the element is proportional to the elemental area, (r + r t )drdq, multiplied by the spatial rate. The probability of more than one cell landing in the element is of the order (drdq) 2 and becomes zero in the limit. If a cell radius is distributed exponentially with parameter f, then the probability that a cell, landing at distance r, intersects the target region is the survivor function, e fr. Thus the probability that a cell lands in the element and extends to the target region is given as ð3þ ð4þ ð5þ pr ðþ/8ðrþr t Þdrdq e fr : ð6þ [15] The proportionality factor required for equation (6) is obtained from the Poisson distribution as e 8(r+rt)drdq. This factor is omitted from the resulting derivation for notational convenience as the probability in equation (6) is linearized with respect to q when taking the limit dq! 0, since this gives e 8(r+rt)drdq! 1. It is clear from equation (6), that for a fixed elemental angle dq, the associated probability will vary with respect to the radius. Alternatively, for an arbitrary constant, p(r) =p o, on rearranging equation (6), the elemental angle can be made to vary for a given radius to ensure that each element maintains this fixed probability, given as p o dqðþ¼ r : ð7þ 8ðr þ r t Þdr efr [16] The set of elements, varying across all radii in the outer region and having elemental angles for a given radius specified by equation (7), define a set of Bernoulli trials. Each trial has a fixed probability, p o, that a rain cell lands in the elemental area at radius r and also intersects the target region. At a given radius the number of Bernoulli trials is given by 2p/dq(r), and the total number of trials, N tot, can be obtained by integrating over all radii in the outer region, Z 1 82pðr þ r t Þ N tot ¼ e fr dr r¼0 p o ¼ 82p r t p o f þ 1 f 2 : [17] A set of Bernoulli trials of size N tot, from which n o will be successful is distributed according to a binomial distribution with PN ð o ¼ n o Þ ¼ N tot n o p no o ð8þ ð1 p o Þ Ntotno : ð9þ [18] The Poisson distribution is derived as the limit to the binomial distribution as N tot!1and is specified with a rate parameter, v = p o N tot. Consequently, the number of cells in the outer region that have radius large enough to reach the target region is Poisson distributed with parameter v, n no PN o ¼ n o jr c > R xy ¼ n o! en ; ð10þ where the parameter v = 82p ( r t f þ 1 Þ is obtained from f equation (8) Cell Center Conditioned on Intersecting Target [19] The probability that a cell intersecting the target comes from an element with radial increment, dr, and elemental angle, dq, is the ratio of the expected number of cells with center in the element that intersect the target, to the expected number of outer cells that intersect the target. This is defined as ¼ Z 1 r xy¼0 PR xy ¼ r xy jr c > R xy 8 r xy þ r t e fr xy dr xy dq Z 2p : ð11þ 8 r xy þ r t e fr xy dr xy dq 3of5
4 LEONARD ET AL.: TECHNICAL NOTE Figure 3. Comparison of direct algorithm and buffer algorithm for proportion of memory usage and simulation run time. Results are reported for a range of sizes of the target region using different parameters across 12 months. [20] Evaluating the denominator for all angles 0 < q 2p, and all radii (0 < r xy 1), and taking the limit drdq!0 this becomes the continuous distribution, f Rxy r xy jr c > r xy ¼ kfe fr xy þ ð1 kþf 2 r xy e frxy ; ð12þ where the resulting expression in equation (12) represents a mixture of an exponential distribution, Exp[f], and a gamma distribution, Gamma[2, 1 f ], and where the mixture ratio k = f r t /(f r t +1) depends upon the radius of the target region and the parameter for cell radii. For larger target regions and smaller expected cell radii, the distribution in equation (12) will tend toward an exponential distribution. The simulation of a point within the outer region requires a random angle, q, which can be independently sampled from a uniform distribution Q Uniform[0, 2p], giving the coordinates of the point as x ¼ R xy þ r t cos Q ð13þ y ¼ R xy þ r t sin Q: 3.3. Cell Radius Conditioned on Location [21] Given that a cell has landed in an element at a distance r xy from the edge of the target region, for the cell to intersect the region, its radius must be greater than the distance r xy. Given that the cell radius is exponentially distributed and that r xy is the constant defined by equation (12), the remaining distance R c r xy is also exponential, based on the standard properties of this distribution. The distribution of the cell radius is therefore given as f Rc r c jr c > r xy ¼ fe fðr cr xy Þ : ð14þ [22] Simulation of the direct algorithm proceeds as previously outlined for the arrival or storms and occurrence of 4of5 cells over the target region. For the cells that land outside the region yet intersect it, the number of cells is sampled using equation (10), their location is sampled using equation (12), and their radius is sampled using equation (14). All other cells properties such as the starting time, lifetime and intensity are independent of the location of the cell. 4. Results and Discussion [23] To demonstrate that the bias of the proposed algorithm is practically negligible, 100 replicates of 1,000,000 years length were simulated for the Arno Basin case study. The resulting median proportion of rainfall was , with a standard deviation of , which arises because of the finite simulation length. This proportion is not statistically significantly different from a value of 1.0 at the 10% level. [24] To illustrate the efficiency of the direct simulation algorithm a comparison was conducted with the buffer algorithm for the Arno Basin case study. Results for the direct algorithm are presented in Figure 3 as a proportion of the requirements for the buffer algorithm for a range of values of the nondimensional quantity fr t. The results were obtained using values of f reported for the 12 months of the year for the Arno Basin case study and for a range of target regions r t = {20, 40, 60, 80, 100} km. For each month, the radius of the buffer region, r, was specified as the constant corresponding to the 0.99 quantile of the cell radii for that month. The results are reported for two variables, the memory usage and simulation time. Because the results of the direct algorithm are reported as ratios of the buffer algorithm, they can loosely be considered independent of the computational platform. The efficiency is largely independent of the number of simulation years, and a length of 300 years was used in order to measure run times with sufficient accuracy. It is important to note that the measured efficiencies depend to some extent on the details of implementation in the computer code, since there is often a tradeoff between memory usage and computational time. Therefore Figure 3 should be regarded as an indication of the order magnitude of the efficiency and as a qualitative indication of the trend in efficiency with respect to various attributes. [25] Figure 3 shows that the proportion of memory usage across the various months and target regions respectively follow two curves, with scatter attributed to random variation of the simulating process and the variability of parameters other than f for each month. To interpret Figure 3, consider a parameter value of f = 0.1. For for a small radius of the target region, r t = 20 km, the buffer region is comparatively large to the target region, hence the direct simulation method is considerably more efficient. The direct method gives a proportion of memory usage at 0.25 and a proportion of duration at approximately 0.6. For a target region having a large radius of r t = 100 km, the buffer region is not as comparatively large, hence the efficiency of the direct method is reduced. At this radius, the proportion of memory usage is 0.4 and the proportion of duration is approximately 0.8. For this latter case, it is important to note that even a small reduction in the duration of the simulation is significant since a simulation will always maximize the available resources, for example by increasing the simulation length or by computing additional replicates. Similar
5 LEONARD ET AL.: TECHNICAL NOTE comparisons could be made from Figure 3 by considering the variation with respect to the parameter f with respect to a fixed target region. [26] It would be unlikely for the space-time Neyman- Scott model to be applied to regions with radius larger than 100 km due to the requirement that the observed rainfall must be homogeneous over this region. Furthermore, it would be expected that the majority of applications of this model would be to urban catchments that are considerably smaller than r t = 100 km. Therefore Figure 3 demonstrates that the direct simulation algorithm significantly reduces the computational requirement of this model. 5. Conclusions [27] The need to simulate cells outside of a region of interest was demonstrated in order to avoid significant boundary effects. An efficient method was derived for directly simulating these cells without the use of a buffer region. The performance of the algorithm varies with respect to the parameter values and the size of the target region, but in all cases was demonstrated to provide considerable improvement. The significance of the direct method is due to the computationally intensive nature of the model and the desire to increase simulation lengths and the number of replicates of the model. [28] Acknowledgments. The authors acknowledge the Australian Research Council for funding this project through a Discovery Project grant. The authors also wish to thank the three anonymous reviewers for their helpful suggestions that have significantly improved this paper. References Cowpertwait, P. (1995), A generalized spatial-temporal model of rainfall based on a clustered point process, Proc. R. Soc. London, Ser. A, 450(1938), Cowpertwait, P., C. G. Kilsby, and P. E. O Connell (2002), A space-time Neyman-Scott model of rainfall: Empirical analysis of extremes, Water Resour. Res., 38(8), 1131, doi: /2001wr Onof, C., R. E. Chandler, A. Kakou, P. Northrop, H. S. Wheater, and V. Isham (2000), Rainfall modelling using Poisson-cluster processes: A review of developments, Stochastic Environ. Res. Risk Assess., 14(6), , doi: /s Wheater, H. S., R. E. Chandler, C. Onof, V. S. Isham, E. Bellone, C. Yang, D. Lekkas, G. Lourmas, and M. L. Segond (2005), Spatial-temporal rainfall modelling for flood risk estimation, Stochastic Environ. Res. Risk Assess., 19(6), , doi: /s M. F. Lambert and M. Leonard, School of Civil and Environmental Engineering, University of Adelaide, Adelaide, SA 5005, Australia. (mleonard@civeng.adelaide.edu.au) A. V. Metcalfe, School of Mathematical Sciences, University of Adelaide, Adelaide, SA 5005, Australia. 5of5
Fitting the Bartlett-Lewis rainfall model using Approximate Bayesian Computation
22nd International Congress on Modelling and Simulation, Hobart, Tasmania, Australia, 3 to 8 December 2017 mssanz.org.au/modsim2017 Fitting the Bartlett-Lewis rainfall model using Approximate Bayesian
More informationA model-based approach to the computation of area probabilities for precipitation exceeding a certain threshold
21st International Congress on Modelling and Simulation, Gold Coast, Australia, 29 Nov to 4 Dec 2015 www.mssanz.org.au/modsim2015 A model-based approach to the computation of area probabilities for precipitation
More informationOccurrence analysis of daily rainfalls through non-homogeneous Poissonian processes
Nat. Hazards Earth Syst. Sci., 11, 1657 1668, 2011 doi:10.5194/nhess-11-1657-2011 Author(s) 2011. CC Attribution 3.0 License. Natural Hazards and Earth System Sciences Occurrence analysis of daily rainfalls
More informationStochastic Generation Of Point Rainfall Data At Sub-Daily Timescales: A Comparison Of DRIP And NSRP
Stochastic Generation Of Point Rainfall Data At Sub-Daily Timescales: A Comparison Of DRIP And NSRP 1 Frost, A. J., 1 R.Srikanthan and 2 P. S. P. Cowpertwait 1 Cooperative Research Centre for Catchment
More informationSpatial and temporal variability of wind speed and energy over Greece
European Geosciences Union General Assembly 2014 Vienna, Austria, 27 April 02 May 2014 Session HS5.4/CL4.8/ERE1.16 : Hydropower and other renewable sources of energy for a sustainable future: modeling
More informationExtending clustered point process-based rainfall models to a non-stationary climate
Extending clustered point process-based rainfall models to a non-stationary climate Jo Kaczmarska 1, 2 Valerie Isham 2 Paul Northrop 2 1 Risk Management Solutions 2 Department of Statistical Science, University
More informationDaily Rainfall Disaggregation Using HYETOS Model for Peninsular Malaysia
Daily Rainfall Disaggregation Using HYETOS Model for Peninsular Malaysia Ibrahim Suliman Hanaish, Kamarulzaman Ibrahim, Abdul Aziz Jemain Abstract In this paper, we have examined the applicability of single
More informationA doubly stochastic rainfall model with exponentially decaying pulses
Stoch Environ Res Risk Assess (2018) 32:1645 1664 https://doi.org/10.1007/s00477-017-1483-z ORIGINAL PAPER A doubly stochastic rainfall model with exponentially decaying pulses N. I. Ramesh 1 A. P. Garthwaite
More informationStochastic disaggregation of spatial-temporal rainfall with limited data
XXVI General Assembly of the European Geophysical Society Nice, France, 25-3 March 2 Session NP5.3/ Scaling, multifractals and nonlinearity in geophysics: Stochastic rainfall modelling: scaling and other
More informationEvaluating Error of Temporal Disaggregation from Daily into Hourly Rainfall using Heytos Model at Sampean Catchments Area
IPTEK, The Journal for Technology and Science, Vol., No., February Evaluating Error of Temporal Disaggregation from Daily into Hourly Rainfall using Heytos Model at Sampean Catchments Area Entin Hidayah,
More informationLJMU Research Online
LJMU Research Online Abdellatif, M, Atherton, W and Alkhaddar, R Application of the stochastic model for temporal rainfall disaggregation for hydrological studies in north western England http://researchonline.ljmu.ac.uk/2344/
More informationMaximum likelihood estimation for the seasonal Neyman-Scott rectangular pulses model for rainfall
Maximum likelihood estimation for the seasonal Neyman-Scott rectangular pulses model for rainfall A. MONTANARI & A. BRATH Dipartimento di Ingegneria delle Strutture, dei Trasporti, delle Acque, del Rilevamento,
More informationStatistical downscaling methods for climate change impact assessment on urban rainfall extremes for cities in tropical developing countries A review
1 Statistical downscaling methods for climate change impact assessment on urban rainfall extremes for cities in tropical developing countries A review International Conference on Flood Resilience: Experiences
More informationSpecialist rainfall scenarios and software package
Building Knowledge for a Changing Climate Specialist rainfall scenarios and software package Chris Kilsby Ahmad Moaven-Hashemi Hayley Fowler Andrew Smith Aidan Burton Michael Murray University of Newcastle
More informationDongkyun Kim Francisco Olivera Huidae Cho. 1 Introduction
Stoch Environ Res Risk Assess (2013) 27:1601 1610 DOI 10.1007/s00477-013-0696-z ORIGINAL PAPER Effect of the inter-annual variability of rainfall statistics on stochastically generated rainfall time series:
More informationA rainfall disaggregation scheme for sub-hourly time scales: Coupling a. Bartlett-Lewis based model with adjusting procedures
A rainfall disaggregation scheme for sub-hourly time scales: Coupling a Bartlett-Lewis based model with adjusting procedures Panagiotis Kossieris a,*, Christos Makropoulos a, Christian Onof b and Demetris
More informationReview of existing statistical methods for flood frequency estimation in Greece
EU COST Action ES0901: European Procedures for Flood Frequency Estimation (FloodFreq) 3 rd Management Committee Meeting, Prague, 28 29 October 2010 WG2: Assessment of statistical methods for flood frequency
More informationIMPACT OF CLIMATE CHANGE ON URBAN DRAINAGE SYSTEM PERFORMANCE
9 th INTERNATIONAL WORKSHOP on PRECIPITATION IN URBAN AREAS IMPACT OF CLIMATE CHANGE ON URBAN DRAINAGE SYSTEM PERFORMANCE by C. Egger (1), M. Honti (1) and M. Maurer (1) (1) Eawag, The Swiss Federal Institute
More informationA space-time Neyman-Scott model of rainfall: Empirical analysis of extremes
WATER RESOURCES RESEARCH, VOL. 38, NO. 8, 1131, 10.1029/2001WR000709, 2002 A space-time Neyman-Scott model of rainfall: Empirical analysis of extremes P. S. P. Cowpertwait Institute of Information and
More information7. Higher Dimensional Poisson Process
1 of 5 7/16/2009 7:09 AM Virtual Laboratories > 14 The Poisson Process > 1 2 3 4 5 6 7 7 Higher Dimensional Poisson Process The Process The Poisson process can be defined in higher dimensions, as a model
More informationLearning Objectives for Stat 225
Learning Objectives for Stat 225 08/20/12 Introduction to Probability: Get some general ideas about probability, and learn how to use sample space to compute the probability of a specific event. Set Theory:
More informationAnonymous Referee #1 Received and published: 21 June 2018
Interactive comment on A Hybrid Stochastic Rainfall Model That Reproduces Rainfall Characteristics at Hourly through Yearly Time Scale by Jeongha Park et al. (RC1) Anonymous Referee #1 Received and published:
More informationRadar-raingauge data combination techniques: A revision and analysis of their suitability for urban hydrology
9th International Conference on Urban Drainage Modelling Belgrade 2012 Radar-raingauge data combination techniques: A revision and analysis of their suitability for urban hydrology L. Wang, S. Ochoa, N.
More informationSimple closed form formulas for predicting groundwater flow model uncertainty in complex, heterogeneous trending media
WATER RESOURCES RESEARCH, VOL. 4,, doi:0.029/2005wr00443, 2005 Simple closed form formulas for predicting groundwater flow model uncertainty in complex, heterogeneous trending media Chuen-Fa Ni and Shu-Guang
More informationModelling and simulation of seasonal rainfall
Modelling and simulation of seasonal rainfall Phil Howlett a, Julia Piantadosi a,1,, Jonathan Borwein b a Scheduling and Control Group, Centre for Industrial and Applied Mathematics (CIAM), University
More informationThe structure of convective rain cells at mid-latitudes
Advances in Geosciences, 7, 31 35, 26 SRef-ID: 168-7359/adgeo/26-7-31 European Geosciences Union 26 Author(s). This work is licensed under a Creative Commons License. Advances in Geosciences The structure
More informationNowcasting for New Zealand
ATMOSPHERIC SCIENCE LETTERS Atmos. Sci. Let. 6: 35 39 (2005) Published online in Wiley InterScience (www.interscience.wiley.com). DOI: 10.1002/asl.88 Nowcasting for New Zealand Warren Gray, 1 * Howard
More informationProbability Method in Civil Engineering Prof. Dr. Rajib Maity Department of Civil Engineering Indian Institute of Technology, Kharagpur
Probability Method in Civil Engineering Prof. Dr. Rajib Maity Department of Civil Engineering Indian Institute of Technology, Kharagpur Lecture No. # 34 Probability Models using Discrete Probability Distributions
More information1 Numbers. exponential functions, such as x 7! a x ; where a; x 2 R; trigonometric functions, such as x 7! sin x; where x 2 R; ffiffi x ; where x 0:
Numbers In this book we study the properties of real functions defined on intervals of the real line (possibly the whole real line) and whose image also lies on the real line. In other words, they map
More informationThe effect of spatial rainfall variability on streamflow prediction for a south-eastern Australian catchment
19th International Congress on Modelling and Simulation, Perth, Australia, 12 16 December 2011 http://mssanz.org.au/modsim2011 The effect of spatial rainfall variability on streamflow prediction for a
More informationA Regionalized Stochastic Rainfall Model for the Generation of High Resolution Data in Peninsular Malaysia
Modern Applied Science; Vol. 10, No. 5; 2016 ISSN 1913-1844 E-ISSN 1913-1852 Published by Canadian Center of Science and Education A Regionalized Stochastic Rainfall Model for the Generation of High Resolution
More informationEvaluation and correction of uncertainty due to Gaussian approximation in radar rain gauge merging using kriging with external drift
Evaluation and correction of uncertainty due to Gaussian approximation in radar rain gauge merging using kriging with external drift F. Cecinati* 1, O. Wani 2,3, M. A. Rico-Ramirez 1 1 University of Bristol,
More informationSTUDY ON THE PRECISION OF 1-MINUTE X-BAND MP RADAR RAINFALL DATA IN A SMALL URBAN WATERSHED
Y. Yonese, et al., Int. J. Sus. Dev. Plann. Vol. 13, No. 4 (2018) 614 625 STUDY ON THE PRECISION OF 1-MINUTE X-BAND MP RADAR RAINFALL DATA IN A SMALL URBAN WATERSHED YOSHITOMO YONESE 1, AKIRA KAWAMURA
More informationEffective unsaturated hydraulic conductivity for one-dimensional structured heterogeneity
WATER RESOURCES RESEARCH, VOL. 41, W09406, doi:10.1029/2005wr003988, 2005 Effective unsaturated hydraulic conductivity for one-dimensional structured heterogeneity A. W. Warrick Department of Soil, Water
More informationProbabilistic forecasting for urban water management: A case study
9th International Conference on Urban Drainage Modelling Belgrade 212 Probabilistic forecasting for urban water management: A case study Jeanne-Rose Rene' 1, Henrik Madsen 2, Ole Mark 3 1 DHI, Denmark,
More informationThe stability of planets in the Alpha Centauri system
The stability of planets in the Alpha Centauri system Paul A. Wiegert 1, Matt Holman 2 1 Department of Astronomy, University of Toronto, Toronto, Canada 2 Canadian Institute for Theoretical Astrophysics,
More informationFractal interpolation of rain rate time series
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 109,, doi:10.1029/2004jd004717, 2004 Fractal interpolation of rain rate time series Kevin S. Paulson Radio Communications Research Unit, Rutherford Appleton Laboratory,
More informationWhat is MuDRain (Multivariate disaggregation of rainfall)? # $ K
What is MuDRain (Multivariate disaggregation of rainfall)? # $ K MuDRain is a methodology for spatial-temporal disaggregation of rainfall. It involves the combination of several univariate and multivariate
More informationStochastic Hydrology. a) Data Mining for Evolution of Association Rules for Droughts and Floods in India using Climate Inputs
Stochastic Hydrology a) Data Mining for Evolution of Association Rules for Droughts and Floods in India using Climate Inputs An accurate prediction of extreme rainfall events can significantly aid in policy
More informationPLANNED UPGRADE OF NIWA S HIGH INTENSITY RAINFALL DESIGN SYSTEM (HIRDS)
PLANNED UPGRADE OF NIWA S HIGH INTENSITY RAINFALL DESIGN SYSTEM (HIRDS) G.A. Horrell, C.P. Pearson National Institute of Water and Atmospheric Research (NIWA), Christchurch, New Zealand ABSTRACT Statistics
More informationDraft Water Resources Management Plan 2019 Annex 3: Supply forecast Appendix B: calibration of the synthetic weather generator
Draft Water Resources Management Plan 2019 Annex 3: Supply forecast Appendix B: calibration of the synthetic weather generator November 30, 2017 Version 1 Introduction This appendix contains calibration
More informationSimulation of 6-hourly rainfall and temperature by two resampling schemes
Journal of Hydrology 273 (2003) 69 80 www.elsevier.com/locate/jhydrol Simulation of 6-hourly rainfall and temperature by two resampling schemes R. Wójcik 1, T.A. Buishand* Royal Netherlands Meteorological
More informationA New Probabilistic Rational Method for design flood estimation in ungauged catchments for the State of New South Wales in Australia
21st International Congress on Modelling and Simulation Gold Coast Australia 29 Nov to 4 Dec 215 www.mssanz.org.au/modsim215 A New Probabilistic Rational Method for design flood estimation in ungauged
More informationThe Union and Intersection for Different Configurations of Two Events Mutually Exclusive vs Independency of Events
Section 1: Introductory Probability Basic Probability Facts Probabilities of Simple Events Overview of Set Language Venn Diagrams Probabilities of Compound Events Choices of Events The Addition Rule Combinations
More informationURBAN DRAINAGE MODELLING
9th International Conference URBAN DRAINAGE MODELLING Evaluating the impact of climate change on urban scale extreme rainfall events: Coupling of multiple global circulation models with a stochastic rainfall
More informationSurface Hydrology Research Group Università degli Studi di Cagliari
Surface Hydrology Research Group Università degli Studi di Cagliari Evaluation of Input Uncertainty in Nested Flood Forecasts: Coupling a Multifractal Precipitation Downscaling Model and a Fully-Distributed
More informationA Comparison of Alternatives for Daily to Sub-Daily Rainfall Disaggregation
18 th World IMACS / MODSIM Congress, Cairns, Australia 13-17 July 2009 http://mssanz.org.au/modsim09 A Comparison of Alternatives for Daily to Sub-Daily Rainfall Disaggregation A.Pui 1, A. Sharma 1 and
More informationProbability and Statistics
Kristel Van Steen, PhD 2 Montefiore Institute - Systems and Modeling GIGA - Bioinformatics ULg kristel.vansteen@ulg.ac.be Chapter 3: Parametric families of univariate distributions CHAPTER 3: PARAMETRIC
More informationOn random sequence spacing distributions
On random sequence spacing distributions C.T.J. Dodson School of Mathematics, Manchester University Manchester M6 1QD, UK ctdodson@manchester.ac.uk Abstract The random model for allocation of an amino
More informationImpact of climate change on Australian flood risk: A review of recent evidence
Impact of climate change on Australian flood risk: A review of recent evidence 30/5/2018 FMA Conference, Gold Coast S Westra 1, B Bennett 1, J Evans 2, F Johnson 3, M Leonard 1, A Sharma 3, C Wasko 4,
More informationPrep for College Algebra
Prep for College Algebra This course covers the topics outlined below. You can customize the scope and sequence of this course to meet your curricular needs. Curriculum (219 topics + 85 additional topics)
More informationPrep for College Algebra with Trigonometry
Prep for College Algebra with Trigonometry This course covers the topics outlined below. You can customize the scope and sequence of this course to meet your curricular needs. Curriculum (246 topics +
More informationA Correlation of. Pearson. Mathematical Ideas. to the. TSI Topics
A Correlation of Pearson 2016 to the A Correlation of 2016 Table of Contents Module M1. Linear Equations, Inequalities, and Systems... 1 Module M2. Algebraic Expressions and Equations (Other Than Linear)...
More informationThe Australian Operational Daily Rain Gauge Analysis
The Australian Operational Daily Rain Gauge Analysis Beth Ebert and Gary Weymouth Bureau of Meteorology Research Centre, Melbourne, Australia e.ebert@bom.gov.au Daily rainfall data and analysis procedure
More informationDesign Flood Estimation in Ungauged Catchments: Quantile Regression Technique And Probabilistic Rational Method Compared
Design Flood Estimation in Ungauged Catchments: Quantile Regression Technique And Probabilistic Rational Method Compared N Rijal and A Rahman School of Engineering and Industrial Design, University of
More informationTrivariate copulas for characterisation of droughts
ANZIAM J. 49 (EMAC2007) pp.c306 C323, 2008 C306 Trivariate copulas for characterisation of droughts G. Wong 1 M. F. Lambert 2 A. V. Metcalfe 3 (Received 3 August 2007; revised 4 January 2008) Abstract
More informationStatistics: Learning models from data
DS-GA 1002 Lecture notes 5 October 19, 2015 Statistics: Learning models from data Learning models from data that are assumed to be generated probabilistically from a certain unknown distribution is a crucial
More informationA Hybrid Stochastic Rainfall Model That Reproduces Rainfall Characteristics at Hourly through Yearly Time Scale
A Hybrid Stochastic Rainfall Model That Reproduces Rainfall Characteristics at Hourly through Yearly Time Scale Jeongha Park 1, Christian Onof 2, and Dongkyun Kim 1 1 Department of Civil Engineering, Hongik
More informationII. Introduction to probability, 2
GEOS 33000/EVOL 33000 5 January 2006 updated January 10, 2006 Page 1 II. Introduction to probability, 2 1 Random Variables 1.1 Definition: A random variable is a function defined on a sample space. In
More informationASM Study Manual for Exam P, Second Edition By Dr. Krzysztof M. Ostaszewski, FSA, CFA, MAAA Errata
ASM Study Manual for Exam P, Second Edition By Dr. Krzysztof M. Ostaszewski, FSA, CFA, MAAA (krzysio@krzysio.net) Errata Effective July 5, 3, only the latest edition of this manual will have its errata
More informationEighteenth International Water Technology Conference, IWTC18 Sharm ElSheikh, March 2015
SPATIOTEMPORAL EVALUATION OF GLOBAL PRECIPITATION MAPPING -GSMAP AT BASIN SCALE IN SAGAMI RIVER, JAPAN S. Takegawa 1, K. Takido 2, ando. Saavedra 3 1 Tokyo Institute of Technology, Tokyo,koizumi.s.ae@m.titech.ac.jp
More informationMonte Carlo Studies. The response in a Monte Carlo study is a random variable.
Monte Carlo Studies The response in a Monte Carlo study is a random variable. The response in a Monte Carlo study has a variance that comes from the variance of the stochastic elements in the data-generating
More informationb. ( ) ( ) ( ) ( ) ( ) 5. Independence: Two events (A & B) are independent if one of the conditions listed below is satisfied; ( ) ( ) ( )
1. Set a. b. 2. Definitions a. Random Experiment: An experiment that can result in different outcomes, even though it is performed under the same conditions and in the same manner. b. Sample Space: This
More informationREQUIREMENTS FOR WEATHER RADAR DATA. Review of the current and likely future hydrological requirements for Weather Radar data
WORLD METEOROLOGICAL ORGANIZATION COMMISSION FOR BASIC SYSTEMS OPEN PROGRAMME AREA GROUP ON INTEGRATED OBSERVING SYSTEMS WORKSHOP ON RADAR DATA EXCHANGE EXETER, UK, 24-26 APRIL 2013 CBS/OPAG-IOS/WxR_EXCHANGE/2.3
More informationArea precipitation probabilities derived from point forecasts for operational weather and warning service applications
Quarterly Journal of the Royal Meteorological Society Q. J. R. Meteorol. Soc. 00: 1 9 (0000) Area precipitation probabilities derived from point forecasts for operational weather and warning service applications
More informationA note on in situ estimates of sorption using push-pull tests
WATER RESOURCES RESEARCH, VOL. 41,, doi:10.1029/2004wr003382, 2005 A note on in situ estimates of sorption using push-pull tests Giorgio Cassiani Dipartimento di Scienze Geologiche e Geotecnologie, Università
More informationAlgebra 2. Curriculum (524 topics additional topics)
Algebra 2 This course covers the topics shown below. Students navigate learning paths based on their level of readiness. Institutional users may customize the scope and sequence to meet curricular needs.
More informationRainfall Disaggregation Methods: Theory and Applications
Workshop on Statistical and Mathematical Methods for Hydrological Analysis Rome, May 23 Universitá degli Studi di Roma La Sapienza CNR - IRPI Istituto di Ricerca per la Protezione Idrogeologica GNDCI Gruppo
More informationSources of uncertainty in estimating suspended sediment load
136 Sediment Budgets 2 (Proceedings of symposium S1 held during the Seventh IAHS Scientific Assembly at Foz do Iguaçu, Brazil, April 2005). IAHS Publ. 292, 2005. Sources of uncertainty in estimating suspended
More informationGenerating projected rainfall time series at sub-hourly time scales using statistical and stochastic downscaling methodologies
Generating projected rainfall time series at sub-hourly time scales using statistical and stochastic downscaling methodologies S. Molavi 1, H. D. Tran 1,2, N. Muttil 1 1 School of Engineering and Science,
More informationRainfall Disaggregation Methods: Theory and Applications
Workshop on Statistical and Mathematical Methods for Hydrological Analysis Rome, May 23 Universitá degli Studi di Roma La Sapienza CNR - IRPI Istituto di Ricerca per la Protezione Idrogeologica GNDCI Gruppo
More informationCh3. Generating Random Variates with Non-Uniform Distributions
ST4231, Semester I, 2003-2004 Ch3. Generating Random Variates with Non-Uniform Distributions This chapter mainly focuses on methods for generating non-uniform random numbers based on the built-in standard
More informationMULTIDIMENSIONAL COVARIATE EFFECTS IN SPATIAL AND JOINT EXTREMES
MULTIDIMENSIONAL COVARIATE EFFECTS IN SPATIAL AND JOINT EXTREMES Philip Jonathan, Kevin Ewans, David Randell, Yanyun Wu philip.jonathan@shell.com www.lancs.ac.uk/ jonathan Wave Hindcasting & Forecasting
More informationMeasures Also Significant Factors of Flood Disaster Reduction
Non-Structual Measures Also Significant Factors of Flood Disaster Reduction Babiaková Gabriela, Leškov ková Danica Slovak Hydrometeorological Institute, Bratislava Hydrological Forecasts and Warning Department
More informationRADAR Rainfall Calibration of Flood Models The Future for Catchment Hydrology? A Case Study of the Stanley River catchment in Moreton Bay, Qld
RADAR Rainfall Calibration of Flood Models The Future for Catchment Hydrology? A Case Study of the Stanley River catchment in Moreton Bay, Qld A 1 Principal Engineer, Water Technology Pty Ltd, Brisbane,
More informationDifference of scattering geometrical optics components and line integrals of currents in modified edge representation
RADIO SCIENCE, VOL. 47,, doi:0.029/20rs004899, 202 Difference of scattering geometrical optics components and line integrals of currents in modified edge representation Pengfei Lu and Makoto Ando Received
More informationInternational Journal of World Research, Vol - 1, Issue - XVI, April 2015 Print ISSN: X
(1) ESTIMATION OF MAXIMUM FLOOD DISCHARGE USING GAMMA AND EXTREME VALUE FAMILY OF PROBABILITY DISTRIBUTIONS N. Vivekanandan Assistant Research Officer Central Water and Power Research Station, Pune, India
More informationLecture 4: Frobenius Series about Regular Singular Points
Introductory lecture notes on Partial Differential Equations - c Anthony Peirce. Not to be copied, used, or revised without explicit written permission from the copyright owner. 1 Lecture 4: Frobenius
More informationUncertainty in merged radar - rain gauge rainfall products
Uncertainty in merged radar - rain gauge rainfall products Francesca Cecinati University of Bristol francesca.cecinati@bristol.ac.uk Supervisor: Miguel A. Rico-Ramirez This project has received funding
More informationExtreme rainstorms: Comparing regional envelope curves to stochastically generated events
WATER RESOURCES RESEARCH, VOL. 48, W01509, doi:10.1029/2011wr010515, 2012 Extreme rainstorms: Comparing regional envelope curves to stochastically generated events A. Viglione, 1 A. Castellarin, 2 M. Rogger,
More informationASM Study Manual for Exam P, First Edition By Dr. Krzysztof M. Ostaszewski, FSA, CFA, MAAA Errata
ASM Study Manual for Exam P, First Edition By Dr. Krzysztof M. Ostaszewski, FSA, CFA, MAAA (krzysio@krzysio.net) Errata Effective July 5, 3, only the latest edition of this manual will have its errata
More informationMixture distributions in Exams MLC/3L and C/4
Making sense of... Mixture distributions in Exams MLC/3L and C/4 James W. Daniel Jim Daniel s Actuarial Seminars www.actuarialseminars.com February 1, 2012 c Copyright 2012 by James W. Daniel; reproduction
More informationCluster Analysis using SaTScan. Patrick DeLuca, M.A. APHEO 2007 Conference, Ottawa October 16 th, 2007
Cluster Analysis using SaTScan Patrick DeLuca, M.A. APHEO 2007 Conference, Ottawa October 16 th, 2007 Outline Clusters & Cluster Detection Spatial Scan Statistic Case Study 28 September 2007 APHEO Conference
More information*Corresponding author address: Charles Barrere, Weather Decision Technologies, 1818 W Lindsey St, Norman, OK
P13R.11 Hydrometeorological Decision Support System for the Lower Colorado River Authority *Charles A. Barrere, Jr. 1, Michael D. Eilts 1, and Beth Clarke 2 1 Weather Decision Technologies, Inc. Norman,
More informationProbability Models. 4. What is the definition of the expectation of a discrete random variable?
1 Probability Models The list of questions below is provided in order to help you to prepare for the test and exam. It reflects only the theoretical part of the course. You should expect the questions
More informationEstimating Probable Maximum Precipitation for Linau River Basin in Sarawak
Estimating Probable Maximum Precipitation for Linau River Basin in Sarawak M. Hussain, S. Nadya and F.J. Chia Hydropower Development Division, Sarawak Energy Berhad, 93050, Kuching, Sarawak, Malaysia Date
More informationCopyright 2018 UC Regents and ALEKS Corporation. ALEKS is a registered trademark of ALEKS Corporation. 2/10
Prep for Calculus This course covers the topics outlined below. You can customize the scope and sequence of this course to meet your curricular needs. Curriculum (281 topics + 125 additional topics) Real
More informationChapter 8: Continuous Probability Distributions
Chapter 8: Continuous Probability Distributions 8.1 Introduction This chapter continued our discussion of probability distributions. It began by describing continuous probability distributions in general,
More informationB r Solved Problems Magnetic Field of a Straight Wire
(4) Equate Iencwith d s to obtain I π r = NI NI = = ni = l π r 9. Solved Problems 9.. Magnetic Field of a Straight Wire Consider a straight wire of length L carrying a current I along the +x-direction,
More informationCorrelation between lightning types
GEOPHYSICAL RESEARCH LETTERS, VOL. 34,, doi:10.1029/2007gl029476, 2007 Correlation between lightning types J. L. Lapp 1 and J. R. Saylor 1 Received 25 January 2007; revised 21 February 2007; accepted 20
More informationBALTRAD tailored ender-user product: Risk assessment map for urban drainage management
BALTRAD tailored ender-user product: Risk assessment map for urban drainage management Authors: Jesper E. Nielsen, Michael R. Rasmussen (AAU) Date: December 2012 BALTRAD Document: BALTRAD+ W4-Risk assessment
More informationAdvanced Structural Analysis EGF Cylinders Under Pressure
Advanced Structural Analysis EGF316 4. Cylinders Under Pressure 4.1 Introduction When a cylinder is subjected to pressure, three mutually perpendicular principal stresses will be set up within the walls
More informationIntroduction to Error Analysis
Introduction to Error Analysis Part 1: the Basics Andrei Gritsan based on lectures by Petar Maksimović February 1, 2010 Overview Definitions Reporting results and rounding Accuracy vs precision systematic
More informationORE Open Research Exeter
ORE Open Research Exeter TITLE Statistical Downscaling Methods for Climate Change Impact Assessment on Urban Rainfall Extremes for Cities in Tropical Developing Countries a Review AUTHORS Mugume, Seith;
More informationStandards for Mathematical Practice. Ratio and Proportional Relationships
North Carolina Standard Course of Study Sixth Grade Mathematics 1. Make sense of problems and persevere in solving them. 2. Reason abstractly and quantitatively. 3. Construct viable arguments and critique
More informationModule 7 SEISMIC HAZARD ANALYSIS (Lectures 33 to 36)
Lecture 35 Topics Module 7 SEISMIC HAZARD ANALYSIS (Lectures 33 to 36) 7.4.4 Predictive Relationships 7.4.5 Temporal Uncertainty 7.4.6 Poisson Model 7.4.7 Other Models 7.4.8 Model Applicability 7.4.9 Probability
More informationSupport Vector Regression with Automatic Accuracy Control B. Scholkopf y, P. Bartlett, A. Smola y,r.williamson FEIT/RSISE, Australian National University, Canberra, Australia y GMD FIRST, Rudower Chaussee
More informationTesting for a Change in Mean under Fractional Integration
J. Time Ser. Econom. 206; aop Fabrizio Iacone*, Stephen J. Leybourne and A.M. Robert Taylor Testing for a Change in Mean under Fractional Integration DOI 0.55/jtse-205-0006 Abstract: We consider testing
More informationACTEX CAS EXAM 3 STUDY GUIDE FOR MATHEMATICAL STATISTICS
ACTEX CAS EXAM 3 STUDY GUIDE FOR MATHEMATICAL STATISTICS TABLE OF CONTENTS INTRODUCTORY NOTE NOTES AND PROBLEM SETS Section 1 - Point Estimation 1 Problem Set 1 15 Section 2 - Confidence Intervals and
More informationCE 3710: Uncertainty Analysis in Engineering
FINAL EXAM Monday, December 14, 10:15 am 12:15 pm, Chem Sci 101 Open book and open notes. Exam will be cumulative, but emphasis will be on material covered since Exam II Learning Expectations for Final
More information