BME STUDIES OF STOCHASTIC DIFFERENTIAL EQUATIONS REPRESENTING PHYSICAL LAWS -PART II
|
|
- Lindsey Anderson
- 6 years ago
- Views:
Transcription
1 BME SUDIES OF SOCHASIC DIFFERENIAL EQUAIONS REPRESENING PHYSICAL LAWS -PAR II M.L. Serre and G. Christakos Environmental Modelling Program, Department of Environmental Sciences & Engineering School of Public Health, University of North Carolina, Chapel Hill, NC , USA Classical Geostatistics methods have been designed to use mainly statistical knowledge about natural variables and they lack the ability to incorporate important forms of knowledge like physical laws and scientific theories into the mapping process. On the other hand, the powerful and versatile Bayesian maximum entropy (BME) method of Modern Geostatistics can accomplish such a task, rigorously and efficiently. In this work, BME is used to incorporate the Darcy law of subsurface hydrology in the spatial mapping of a hydraulic head field. he hydraulic map thus obtained is physically meaningful as well as numerically more accurate than that obtained using classical methods (e.g., kriging). Moreover, taking advantage of the Darcy law the BME hydraulic head mapping can involve other related soil properties, like hydraulic conductivity. he approach leads to very accurate hydraulic head solutions, and may be also applied to study the inverse problem, in which one seeks to estimate hydraulic conductivity from hydraulic head measurements.. INRODUCION he Bayesian Maximum Entropy (BME) method of Modern Geostatistics is a method which offers a rigorous framework to account for a large class of general and specificatory knowledge bases [-6]. he general knowledge G includes physical laws, scientific theories, empirical relationships, and statistical moments. he specificatory knowledge S includes hard data (exact measurement) and soft data (such as intervals, probability assessments, uncertain observations and fuzzy sets) that are specific to the mapping situation. Previous studies have shown that the BME method is very useful and numerically efficient in a number of practical situations in which various forms of soft data are available [5, 6]. In the first part of this work ([]) we discussed two methods used by BME analysis in order to account for physical knowledge in the form of physical laws. In this work, the BME formalism is used to incorporate the Darcy law of groundwater flow in porous media, which relates changes in hydraulic head with soil properties such as the hydraulic conductivity or resistivity. Darcy law is considered part of the general knowledge G which leads to a prior G -based pdf that describes the joint spatial distribution of hydraulic head and resistivity. hen, the specificatory knowledge S is used to update the prior pdf leading to the posterior pdf of the hydraulic head map. Numerical examples show that, by accounting for the Darcy law BME provides predictions that are physically meaningful and numerically more accurate than classical data analysis that does not account for the physical law.. INCORPORAING DARCY LAW IN BME MAPPING Consider the Darcy law of one-dimensional flow usually expressed as K(s) H(s) = q(s), where the specific discharge q(s) is considered deterministic while the hydraulic head H(s) and the hydraulic conductivity K(s) are random fields. It is convenient to define the hydraulic resistivity R(s) = / K(s) and rewrite Darcy law as dh(s) / ds = q(s) R(s). () Let H be the known value of H(s) at the spatial origin s =, let H map = [H,..., H nh ] be the vector of random variables representing the hydraulic head H(s) at points s,...,s nh, and assume
2 that R data = [R nh +,..., R n ] is the vector of random variables representing the hydraulic resistivity R(s) at points s nh +,...,s n. At the prior stage of the BME analysis we seek to obtain the prior pdf f G ;s...s n ) given general knowledge G (which includes the Darcy law), where H map and R data represent a realization of the random vectors H map and R data, respectively. he general knowledge G is transformed in a suitable set of moment equations as follows dh map dr data g α ) f G ;s...s n ) = g α ( α =,..., N α ), where g α is a set of functions chosen such that their stochastic expectations g α are either known statistical moments or they can be inferred from Darcy law. For example, we assume that the hydraulic resistivity has a known mean m R (s) = R(s), non-centered covariance C R (s,s' ) = R(s)R(s' ), and th-order non-centered moment m R, (s) = R (s). hese moments are incorporated in the general knowledge by letting g i (R i ;s i ) = R i and g i ) = m R ), i J R = {n H +,...,n}; g ij (R i,r j ;s i ) = R i R j and g ij ) = C R ) for i, j J R ; and finally g n+i (R i ;s i ) = R i and g n+i ) = m R, ) for n + i J R. On the other hand the moments describing the hydraulic head are not known a priori, but they may be inferred from Darcy law, Eq. (). Using the Darcy law we may obtain expressions for the mean m H (s) = H(s) and noncentered covariance function C H (s,s' ) = H(s)H(s' ) of the hydraulic head, as well as the noncentered cross-covariance function C RH (s,s' ) = R(s)H(s' ). Clearly, if we take the stochastic expectation of Eq. () and solve with respect to the boundary condition H we obtain m H (s) = H o s duq(u)m R (u), which is readily incorporated into G by choosing g i (H i ;s i ) = H i and g i ) = m H ) for i J H, where J R = {n H +,...,n}. Similarly, multiplying Eq. () by R(s' ), taking the stochastic expectation and solving, we obtain s C HR (s,s' ) = H m R (s' ) duq(u)c R (u,s' ), which we incorporate in G by choosing g ij (H i,r j ;s i ) = H i R j and g ij ) = C HR ) for i J H, j J R. Other moments may be obtained in a similar way. For example, multiplying Eq. () by H(s' ) and taking the stochastic expectation leads C H (s,s' ) = s du s' du' q(u)q(u' )C R (u,u' ) H + H [m H (s) + m H (s' )], which we also incorporate in G by choosing g ij (H i,h j ;s i ) = H i H j and g ij ) = C H ) for i, j J H. he selected g α functions, summarized in able, represent a set of constraints which incorporates the general knowledge consisting of statistical moments and the Darcy law. Using these constraints under the epistemic goal of information maximization leads to the following solution for the prior pdf [] f G ;s...s n ) = Z exp[ µ i H i + µ i R i + µ ij H i H j + µ ij R i R j i J H i J R i J H j J H i J R j J R + µ ij H i R i + µ i+n R i ] () i J H j J R where the Lagrange coefficients µ i and µ ij are calculated by solving the system of equations obtained from able (). In the posterior stage of the BME method, the specificatory knowledge S is used to update the prior pdf and lead to the posterior pdf f K (H k ) which completely describes the hydraulic i J R
3 able : General knowledge constraints α g α g α α = i, i J H g i (H i ;s i ) = H i g i ) = m H ) α = i, i J R g i (R i ;s i ) = R i g i ) = m R ) α = (ij), i J H, j J H g ij (H i,h j ;s i ) = H i H j g ij ) = C H ) α = (ij), i J R, j J R g ij (R i,r j ;s i ) = R i R j ; g ij ) = C R ) α = (ij), i J H, j J R g ij (H i,r j ;s i ) = H i R j ; g ij ) = C HR ) α = n + i, i J R g n+i (R i ;s i ) = R i ; g n+i ) = m R, ) head distribution. Let H hard = [H,..., H mhh ], H soft = [H mhh +,..., H m H ], and H k represent H(s) at the hard data points, soft data points and estimation points, respectively, such that H map = [H hard H soft H k ]. Similarly let R hard = [R,..., R mrh ] and R soft = [R mrh +,..., R m R ] represent R(s) at the hard and soft data points, respectively, such that R data = [R hard R soft ]. For soft data of the interval type, the BME posterior pdf is given by ([,]) f K (H k ) = A u H dh soft l H u R dr soft l R f G ), where l H and u H are the lower and upper bounds of the soft data for the hydraulic head such that l H H soft u H, and l R and u R are the lower and upper bounds of the soft data for the hydraulic resistivity so that l R R soft u R.. SIMULAED CASE SUDY In this case study we consider the R(s) profile shown in Fig a. hese values were generated at grid-points separated by a distance of.5 m, using the exponential covariance model c R (s,s' ) = c R, exp( s' s /a R ) with c R, = (sec/ mm) and a R = 6 m, and a mean of m R = sec/ mm. he histogram of R(s) values (Fig. b) is approximately symmetric with a kurtosis coefficient of m R, / c R,, while the corresponding K(s)=/ R(s) values have an asymmetric distribution (Fig. c). he hydraulic head H(s) values calculated by solving Darcy law, Eq. (), with H = m and q(s) = q =.5 mm / sec are shown in Fig. a. he mean hydraulic head m H (s) = H o q m R s is shown with a dashed line in Fig. a, and the histogram of hydraulic head fluctuation H(s) m H (s) is shown with plain lines (Fig. b) at points corresponding to s = 5 m, s = 5 m and s = 75 m. hese histograms were obtained by using, H(s) realizations, and they show that the variance of H(s) increases as s increases. In Fig. b (dotted line) the prior pdf's for the hydraulic head at points s = 5 m, s = 5 m and s = 75 m. Note that these are in excellent agreement with the experimental histogram calculated from the data. In a first test case, we consider the estimation of H(s) m H (s) using sparse measurements for H(s) and R(s). he selected hard and soft (interval) data are shown with triangles and error bars, respectively, in Fig a for R(s), and in Fig a for H(s). In Fig a the estimated head fluctuation profile was derived from simple kriging (SK) using only hard head data; the actual head fluctuation profile is also shown for comparison. Note the poor SK estimates at unobserved locations. he head fluctuation profile in Fig. b was obtained from BME using hard and soft (interval) head data as well as the Darcy law. Note that the BME estimated profile
4 R(s) (in sec/mm) (a) Frequency.5... (b) Frequency.8.6. (c) R (sec/mm) K (mm/sec) Figure : (a) Simulated profile for the hydraulic resistivity R(s). Histograms of (b) the hydraulic resistivity R(s), and (c) the hydraulic conductivity K(s) = / R(s). (a) Hydraulic head H(s) Mean hydraulic head m H (s) (b) Calculated from data BME prior pdf H(s) (in meter) Probability density function.5.5 s=5 m s=5 m s=75 m H(s) m H (s) (in meters) Figure : (a) Simulated H(s) profile using hard and soft data denoted by triangles and error bars, respectively. (b) H(s) histogram at different spatial locations s. is a substantial improvement over the classical SK method. Indeed, by being able to incorporate the Darcy law, the BME method allowed to account for hard and soft (interval) hydraulic conductivity data, as well. In a second case study we consider the problem of estimating the hydraulic resistivity R(s) using sparse measurements of R(s) and H(s). his is sometimes called the inverse problem. While it is possible to obtain accurate measurements of the hydraulic head H(s), it is difficult to directly measure the hydraulic resistivity. herefore, in case study, we consider that the measurements for H(s) are the hard data shown with triangles in Fig. a, while the information for R(s) consists of the soft interval data shown with error bars in Fig. b. For illustration purposes we calculate BME estimates for different mapping situations. In the first situation we assume that the specificatory knowledge consists of only the soft data for R(s), and we show the corresponding BME profile for R(s) in Fig. b. Note that this type of soft information commonly arises when estimating hydraulic resistivities, and the BME method provides a rigorous framework to account for such soft information. In the second mapping situation we assume that the specificatory knowledge includes both the soft data for R(s), as well as the hard data for H(s) shown with triangles in Fig. a. Using this information we obtain the profile shown in Fig. c, which clearly provides a substantial improvement over that of Fig. b. Note that the profile of Fig. c respects the soft data but also reproduces the important features of the
5 .... Hydraulic Head Fluctuation (in meter) Hydraulic Head Fluctuation (in meter) Figure : Estimated head fluctuation profiles obtained from: (a) SK using hard head data (triangles); (b) BME using hard and soft (interval) head data and the Darcy law. Dotted lines represent the actual head fluctuation profile. (a) 5 (b) H(s) (in meter) (c) 5 R(s) (in sec/mm) 5 5 (d) 5 R(s) (in sec/mm) R(s) (in sec/mm) Figure : (a) Hydraulic head H(s) profile with hard data denoted by triangles. Hydraulic resistivity R(s) obtained from: (b) BME using R(s) data; (c) BME using data for both R(s) and H(s); and (d) BME using only H(s) data. he dotted line represents the true R(s) profile. true resistivity profile shown in dotted line. his unique feature of the BME approach is due to the fact that BME accounts for the Darcy law which relates H(s) and R(s), leading to substantial improvements in the estimation of the hydraulic resistivity. Finally, in Fig. d we show the R(s) profile using only hard H(s) data (triangles in Fig. a). Note that in this case the important features of the hydraulic resistivity profiles are well preserved. ACKNOWLEDGMENS: his work was supported by grants from the National Institute of Environmental Health Sciences (Grant no. P ES598-) the Department of Energy (Grant no. DE-FC9-9SR86), and the Army Research Office (Grant no. DAAG ), and by the Computational Science Graduate Fellowship Program of the Office of Scientific Computing.
6 REFERENCES. CHRISAKOS, G., "A Bayesian/maximum-entropy view to the spatial estimation problem". Mathematical Geology, (7), (99).. CHRISAKOS, G.- Random Field Models in Earth Sciences. Academic Press, San Diego, CA, 99.. CHRISAKOS, G. and X. LI - "Bayesian maximum entropy analysis and mapping: farewell to kriging estimators?" Mathematical Geology, (), 5-6 (998).. CHRISAKOS, G, D.. HRISOPULOS and SERRE, M.L. - BME studies of stochastic differential equation representing physical laws -part I". In Proc. of the 5th Annual Confer. of the Internat. Assoc.for Mathematical Geology, rodheim, Norway (999). 5. SERRE, M. L., P. BOGAER and G. CHRISAKOS. "Latest Computational Results in Spatiotemporal Prediction Using the Bayesian Maximum Entropy Method". In Proc. of the th Annual Confer. of the Internat. Assoc. for Mathematical Geology, Eds. Buccianti A., Nardi G., Potenza R., De Frede Editore, Napoli, 7- (998). 6. SERRE, M. L. and CHRISAKOS, G.- "Modern geostatistics: computational BME analysis in the light of uncertain physical knowledge - the Equus Beds study". Journal of Stochastic Environmental Research and Risk Analysis, (), -6 (999)
BME STUDIES OF STOCHASTIC DIFFERENTIAL EQUATIONS REPRESENTING PHYSICAL LAW
7 VIII. BME STUDIES OF STOCHASTIC DIFFERENTIAL EQUATIONS REPRESENTING PHYSICAL LAW A wide variety of natural processes are described using physical laws. A physical law may be expressed by means of an
More informationENVIRONMENTAL SPATIOTEMPORAL MAPPING AND GROUND WATER FLOW MODELLING USING THE BME AND ST METHODS
ENVIRONMENTAL SPATIOTEMPORAL MAPPING AND GROUND WATER FLOW MODELLING USING THE BME AND ST METHODS Marc Laurent Serre A dissertation submitted to the faculty of the University of North Carolina at Chapel
More informationGLOBAL MAPPING OF OZONE DISTRIBUTION USING SOFT INFORMATION AND DATA FROM TOMS AND SBUV INSTRUMENTS ON THE NIMBUS 7 SATELLITE
GLOBAL MAPPING OF OZONE DISTRIBUTION USING SOFT INFORMATION AND DATA FROM TOMS AND SBUV INSTRUMENTS ON THE NIMBUS 7 SATELLITE George Christakos, Alexander Kolovos, Marc L. Serre, and Chandra Abhishek Center
More informationCOMPUTATIONAL MULTI-POINT BME AND BME CONFIDENCE SETS
VII. COMPUTATIONAL MULTI-POINT BME AND BME CONFIDENCE SETS Until now we have considered the estimation of the value of a S/TRF at one point at the time, independently of the estimated values at other estimation
More informationREFERENCES. Armstrong M Basic linear geostatistics. Springer. Berlin. 153 p.
REFERENCES Armstrong M. 1998. Basic linear geostatistics. Springer. Berlin. 153 p. Armstrong M., Carignan J. 1997. Géostatistique linéaire : Application au Domaine Minier. Les Presses de l'ecole des Mines.
More informationMulti-Point BME Space/Time Mapping of Environmental Variables George Christakos & ^.C 27J , [/..S.A EM ail: unc.edu.
Multi-Point BME Space/Time Mapping of Environmental Variables George Christakos & ^.C 27J99-7400, [/..S.A EM ail: george_christakos@ unc.edu Abstract The multi-point Bayesian Maximum Entropy (BME) approach
More informationSpatiotemporal Analysis of Solar Radiation for Sustainable Research in the Presence of Uncertain Measurements
Spatiotemporal Analysis of Solar Radiation for Sustainable Research in the Presence of Uncertain Measurements Alexander Kolovos SAS Institute, Inc. alexander.kolovos@sas.com Abstract. The study of incoming
More informationA MATLAB Implementation of the. Minimum Relative Entropy Method. for Linear Inverse Problems
A MATLAB Implementation of the Minimum Relative Entropy Method for Linear Inverse Problems Roseanna M. Neupauer 1 and Brian Borchers 2 1 Department of Earth and Environmental Science 2 Department of Mathematics
More informationSpace-time mapping of soil salinity using probabilistic bayesian maximum entropy
Stoch Envir Res and Risk Ass (2004) 18: 219 227 Ó Springer-Verlag 2004 DOI 10.1007/s00477-004-0177-5 ORIGINAL PAPER A. Douaik M. van Meirvenne T. To th M. Serre Space-time mapping of soil salinity using
More informationSpatial Inference of Nitrate Concentrations in Groundwater
Spatial Inference of Nitrate Concentrations in Groundwater Dawn Woodard Operations Research & Information Engineering Cornell University joint work with Robert Wolpert, Duke Univ. Dept. of Statistical
More informationCovariance Matrix Simplification For Efficient Uncertainty Management
PASEO MaxEnt 2007 Covariance Matrix Simplification For Efficient Uncertainty Management André Jalobeanu, Jorge A. Gutiérrez PASEO Research Group LSIIT (CNRS/ Univ. Strasbourg) - Illkirch, France *part
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 informationAssessment of Hydraulic Conductivity Upscaling Techniques and. Associated Uncertainty
CMWRXVI Assessment of Hydraulic Conductivity Upscaling Techniques and Associated Uncertainty FARAG BOTROS,, 4, AHMED HASSAN 3, 4, AND GREG POHLL Division of Hydrologic Sciences, University of Nevada, Reno
More informationCorrecting Variogram Reproduction of P-Field Simulation
Correcting Variogram Reproduction of P-Field Simulation Julián M. Ortiz (jmo1@ualberta.ca) Department of Civil & Environmental Engineering University of Alberta Abstract Probability field simulation is
More informationInverse Modelling for Flow and Transport in Porous Media
Inverse Modelling for Flow and Transport in Porous Media Mauro Giudici 1 Dipartimento di Scienze della Terra, Sezione di Geofisica, Università degli Studi di Milano, Milano, Italy Lecture given at the
More informationAdvanced analysis and modelling tools for spatial environmental data. Case study: indoor radon data in Switzerland
EnviroInfo 2004 (Geneva) Sh@ring EnviroInfo 2004 Advanced analysis and modelling tools for spatial environmental data. Case study: indoor radon data in Switzerland Mikhail Kanevski 1, Michel Maignan 1
More informationInverting hydraulic heads in an alluvial aquifer constrained with ERT data through MPS and PPM: a case study
Inverting hydraulic heads in an alluvial aquifer constrained with ERT data through MPS and PPM: a case study Hermans T. 1, Scheidt C. 2, Caers J. 2, Nguyen F. 1 1 University of Liege, Applied Geophysics
More informationRandom Processes. DS GA 1002 Probability and Statistics for Data Science.
Random Processes DS GA 1002 Probability and Statistics for Data Science http://www.cims.nyu.edu/~cfgranda/pages/dsga1002_fall17 Carlos Fernandez-Granda Aim Modeling quantities that evolve in time (or space)
More information7 Geostatistics. Figure 7.1 Focus of geostatistics
7 Geostatistics 7.1 Introduction Geostatistics is the part of statistics that is concerned with geo-referenced data, i.e. data that are linked to spatial coordinates. To describe the spatial variation
More informationDownscaling Seismic Data to the Meter Scale: Sampling and Marginalization. Subhash Kalla LSU Christopher D. White LSU James S.
Downscaling Seismic Data to the Meter Scale: Sampling and Marginalization Subhash Kalla LSU Christopher D. White LSU James S. Gunning CSIRO Contents Context of this research Background Data integration
More informationSimultaneous use of hydrogeological and geophysical data for groundwater protection zone delineation by co-conditional stochastic simulations
Simultaneous use of hydrogeological and geophysical data for groundwater protection zone delineation by co-conditional stochastic simulations C. Rentier,, A. Dassargues Hydrogeology Group, Departement
More information6. GRID DESIGN AND ACCURACY IN NUMERICAL SIMULATIONS OF VARIABLY SATURATED FLOW IN RANDOM MEDIA: REVIEW AND NUMERICAL ANALYSIS
Harter Dissertation - 1994-132 6. GRID DESIGN AND ACCURACY IN NUMERICAL SIMULATIONS OF VARIABLY SATURATED FLOW IN RANDOM MEDIA: REVIEW AND NUMERICAL ANALYSIS 6.1 Introduction Most numerical stochastic
More informationReview. DS GA 1002 Statistical and Mathematical Models. Carlos Fernandez-Granda
Review DS GA 1002 Statistical and Mathematical Models http://www.cims.nyu.edu/~cfgranda/pages/dsga1002_fall16 Carlos Fernandez-Granda Probability and statistics Probability: Framework for dealing with
More informationSTAT 730 Chapter 4: Estimation
STAT 730 Chapter 4: Estimation Timothy Hanson Department of Statistics, University of South Carolina Stat 730: Multivariate Analysis 1 / 23 The likelihood We have iid data, at least initially. Each datum
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 informationPATTERN RECOGNITION AND MACHINE LEARNING CHAPTER 2: PROBABILITY DISTRIBUTIONS
PATTERN RECOGNITION AND MACHINE LEARNING CHAPTER 2: PROBABILITY DISTRIBUTIONS Parametric Distributions Basic building blocks: Need to determine given Representation: or? Recall Curve Fitting Binary Variables
More informationEnsemble Kalman filter assimilation of transient groundwater flow data via stochastic moment equations
Ensemble Kalman filter assimilation of transient groundwater flow data via stochastic moment equations Alberto Guadagnini (1,), Marco Panzeri (1), Monica Riva (1,), Shlomo P. Neuman () (1) Department of
More informationMACRODISPERSION AND DISPERSIVE TRANSPORT BY UNSTEADY RIVER FLOW UNDER UNCERTAIN CONDITIONS
MACRODISPERSION AND DISPERSIVE TRANSPORT BY UNSTEADY RIVER FLOW UNDER UNCERTAIN CONDITIONS M.L. Kavvas and L.Liang UCD J.Amorocho Hydraulics Laboratory University of California, Davis, CA 95616, USA Uncertainties
More informationDS-GA 1002 Lecture notes 12 Fall Linear regression
DS-GA Lecture notes 1 Fall 16 1 Linear models Linear regression In statistics, regression consists of learning a function relating a certain quantity of interest y, the response or dependent variable,
More informationEstimation Theory. as Θ = (Θ 1,Θ 2,...,Θ m ) T. An estimator
Estimation Theory Estimation theory deals with finding numerical values of interesting parameters from given set of data. We start with formulating a family of models that could describe how the data were
More informationStatistical Rock Physics
Statistical - Introduction Book review 3.1-3.3 Min Sun March. 13, 2009 Outline. What is Statistical. Why we need Statistical. How Statistical works Statistical Rock physics Information theory Statistics
More informationEfficient geostatistical simulation for spatial uncertainty propagation
Efficient geostatistical simulation for spatial uncertainty propagation Stelios Liodakis University of the Aegean University Hill Mytilene, Greece stelioslio@geo.aegean.gr Phaedon Kyriakidis Cyprus University
More informationThe Generalized Likelihood Uncertainty Estimation methodology
CHAPTER 4 The Generalized Likelihood Uncertainty Estimation methodology Calibration and uncertainty estimation based upon a statistical framework is aimed at finding an optimal set of models, parameters
More informationStochastic Spectral Approaches to Bayesian Inference
Stochastic Spectral Approaches to Bayesian Inference Prof. Nathan L. Gibson Department of Mathematics Applied Mathematics and Computation Seminar March 4, 2011 Prof. Gibson (OSU) Spectral Approaches to
More informationDynamic System Identification using HDMR-Bayesian Technique
Dynamic System Identification using HDMR-Bayesian Technique *Shereena O A 1) and Dr. B N Rao 2) 1), 2) Department of Civil Engineering, IIT Madras, Chennai 600036, Tamil Nadu, India 1) ce14d020@smail.iitm.ac.in
More informationBayesian SAE using Complex Survey Data Lecture 4A: Hierarchical Spatial Bayes Modeling
Bayesian SAE using Complex Survey Data Lecture 4A: Hierarchical Spatial Bayes Modeling Jon Wakefield Departments of Statistics and Biostatistics University of Washington 1 / 37 Lecture Content Motivation
More informationSupporting Information: A hybrid approach to estimating national scale spatiotemporal variability of PM2.5 in the contiguous United States.
Supporting Information: A hybrid approach to estimating national scale spatiotemporal variability of PM2.5 in the contiguous United States. Authors: Bernardo S. Beckerman* Michael Jerrett Marc Serre Randall
More informationThomas Bayes versus the wedge model: An example inference using a geostatistical prior function
Thomas Bayes versus the wedge model: An example inference using a geostatistical prior function Jason M. McCrank, Gary F. Margrave, and Don C. Lawton ABSTRACT The Bayesian inference is used to estimate
More informationA kernel indicator variogram and its application to groundwater pollution
Int. Statistical Inst.: Proc. 58th World Statistical Congress, 2011, Dublin (Session IPS101) p.1514 A kernel indicator variogram and its application to groundwater pollution data Menezes, Raquel University
More informationA conditional simulation of non-normal velocity/pressure fields
A conditional simulation of non-normal velocity/ fields K. R. urley, A. Kareem Department of Civil Engineering University of Florida, ainesville, FL, 36 Department of Civil Engineering and eological Sciences
More informationDeterministic solution of stochastic groundwater flow equations by nonlocalfiniteelements A. Guadagnini^ & S.P. Neuman^
Deterministic solution of stochastic groundwater flow equations by nonlocalfiniteelements A. Guadagnini^ & S.P. Neuman^ Milano, Italy; ^Department of. Hydrology & Water Resources, The University ofarizona,
More informationCS839: Probabilistic Graphical Models. Lecture 7: Learning Fully Observed BNs. Theo Rekatsinas
CS839: Probabilistic Graphical Models Lecture 7: Learning Fully Observed BNs Theo Rekatsinas 1 Exponential family: a basic building block For a numeric random variable X p(x ) =h(x)exp T T (x) A( ) = 1
More informationSub-kilometer-scale space-time stochastic rainfall simulation
Picture: Huw Alexander Ogilvie Sub-kilometer-scale space-time stochastic rainfall simulation Lionel Benoit (University of Lausanne) Gregoire Mariethoz (University of Lausanne) Denis Allard (INRA Avignon)
More informationA Short Note on the Proportional Effect and Direct Sequential Simulation
A Short Note on the Proportional Effect and Direct Sequential Simulation Abstract B. Oz (boz@ualberta.ca) and C. V. Deutsch (cdeutsch@ualberta.ca) University of Alberta, Edmonton, Alberta, CANADA Direct
More informationBME representation of particulate matter distributions in the
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 106, NO. D9, PAGES 9717-9731, MAY 16, 2001 BME representation of particulate matter distributions in the state of California on the' basis of uncertain measurements
More information44 CHAPTER 2. BAYESIAN DECISION THEORY
44 CHAPTER 2. BAYESIAN DECISION THEORY Problems Section 2.1 1. In the two-category case, under the Bayes decision rule the conditional error is given by Eq. 7. Even if the posterior densities are continuous,
More informationStructural Econometrics: Dynamic Discrete Choice. Jean-Marc Robin
Structural Econometrics: Dynamic Discrete Choice Jean-Marc Robin 1. Dynamic discrete choice models 2. Application: college and career choice Plan 1 Dynamic discrete choice models See for example the presentation
More informationCOLLOCATED CO-SIMULATION USING PROBABILITY AGGREGATION
COLLOCATED CO-SIMULATION USING PROBABILITY AGGREGATION G. MARIETHOZ, PH. RENARD, R. FROIDEVAUX 2. CHYN, University of Neuchâtel, rue Emile Argand, CH - 2009 Neuchâtel, Switzerland 2 FSS Consultants, 9,
More informationDynamic Data Driven Simulations in Stochastic Environments
Computing 77, 321 333 (26) Digital Object Identifier (DOI) 1.17/s67-6-165-3 Dynamic Data Driven Simulations in Stochastic Environments C. Douglas, Lexington, Y. Efendiev, R. Ewing, College Station, V.
More informationA Stochastic Collocation based. for Data Assimilation
A Stochastic Collocation based Kalman Filter (SCKF) for Data Assimilation Lingzao Zeng and Dongxiao Zhang University of Southern California August 11, 2009 Los Angeles Outline Introduction SCKF Algorithm
More informationOn dealing with spatially correlated residuals in remote sensing and GIS
On dealing with spatially correlated residuals in remote sensing and GIS Nicholas A. S. Hamm 1, Peter M. Atkinson and Edward J. Milton 3 School of Geography University of Southampton Southampton SO17 3AT
More informationMaking rating curves - the Bayesian approach
Making rating curves - the Bayesian approach Rating curves what is wanted? A best estimate of the relationship between stage and discharge at a given place in a river. The relationship should be on the
More informationBivariate Rainfall and Runoff Analysis Using Entropy and Copula Theories
Entropy 2012, 14, 1784-1812; doi:10.3390/e14091784 Article OPEN ACCESS entropy ISSN 1099-4300 www.mdpi.com/journal/entropy Bivariate Rainfall and Runoff Analysis Using Entropy and Copula Theories Lan Zhang
More informationAnalysis of the temporal and spatial variations of the chloride concentration in the Río Mojotoro basin
Analysis of the temporal and spatial variations of the chloride concentration in the Río Mojotoro basin Jacqueline Köhn, Eduardo Kruse & Jerónimo Ainchil Facultad de Ciencias Astronómicas y Geofísicas
More informationCOLLOCATED CO-SIMULATION USING PROBABILITY AGGREGATION
COLLOCATED CO-SIMULATION USING PROBABILITY AGGREGATION G. MARIETHOZ, PH. RENARD, R. FROIDEVAUX 2. CHYN, University of Neuchâtel, rue Emile Argand, CH - 2009 Neuchâtel, Switzerland 2 FSS Consultants, 9,
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 informationMaximum variance formulation
12.1. Principal Component Analysis 561 Figure 12.2 Principal component analysis seeks a space of lower dimensionality, known as the principal subspace and denoted by the magenta line, such that the orthogonal
More informationThe Expectation Maximization Algorithm
The Expectation Maximization Algorithm Frank Dellaert College of Computing, Georgia Institute of Technology Technical Report number GIT-GVU-- February Abstract This note represents my attempt at explaining
More informationGaussian Processes. Le Song. Machine Learning II: Advanced Topics CSE 8803ML, Spring 2012
Gaussian Processes Le Song Machine Learning II: Advanced Topics CSE 8803ML, Spring 01 Pictorial view of embedding distribution Transform the entire distribution to expected features Feature space Feature
More informationThree data analysis problems
Three data analysis problems Andreas Zezas University of Crete CfA Two types of problems: Fitting Source Classification Fitting: complex datasets Fitting: complex datasets Maragoudakis et al. in prep.
More informationFundamentals of Data Assimila1on
014 GSI Community Tutorial NCAR Foothills Campus, Boulder, CO July 14-16, 014 Fundamentals of Data Assimila1on Milija Zupanski Cooperative Institute for Research in the Atmosphere Colorado State University
More information13 : Variational Inference: Loopy Belief Propagation and Mean Field
10-708: Probabilistic Graphical Models 10-708, Spring 2012 13 : Variational Inference: Loopy Belief Propagation and Mean Field Lecturer: Eric P. Xing Scribes: Peter Schulam and William Wang 1 Introduction
More informationModule 8. Lecture 5: Reliability analysis
Lecture 5: Reliability analysis Reliability It is defined as the probability of non-failure, p s, at which the resistance of the system exceeds the load; where P() denotes the probability. The failure
More informationHierarchical sparse Bayesian learning for structural health monitoring. with incomplete modal data
Hierarchical sparse Bayesian learning for structural health monitoring with incomplete modal data Yong Huang and James L. Beck* Division of Engineering and Applied Science, California Institute of Technology,
More informationTools for Parameter Estimation and Propagation of Uncertainty
Tools for Parameter Estimation and Propagation of Uncertainty Brian Borchers Department of Mathematics New Mexico Tech Socorro, NM 87801 borchers@nmt.edu Outline Models, parameters, parameter estimation,
More informationMachine Learning 4771
Machine Learning 4771 Instructor: Tony Jebara Topic 11 Maximum Likelihood as Bayesian Inference Maximum A Posteriori Bayesian Gaussian Estimation Why Maximum Likelihood? So far, assumed max (log) likelihood
More informationGeostatistics in Hydrology: Kriging interpolation
Chapter Geostatistics in Hydrology: Kriging interpolation Hydrologic properties, such as rainfall, aquifer characteristics (porosity, hydraulic conductivity, transmissivity, storage coefficient, etc.),
More informationA full-bayesian approach to the groundwater inverse problem for steady state flow
WATER RESOURCES RESEARCH, VOL. 36, NO. 8, PAGES 2081 2093, AUGUST 2000 A full-bayesian approach to the groundwater inverse problem for steady state flow Allan D. Woodbury Department of Civil Engineering,
More informationModeling and Interpolation of Non-Gaussian Spatial Data: A Comparative Study
Modeling and Interpolation of Non-Gaussian Spatial Data: A Comparative Study Gunter Spöck, Hannes Kazianka, Jürgen Pilz Department of Statistics, University of Klagenfurt, Austria hannes.kazianka@uni-klu.ac.at
More informationBayesian Learning in Undirected Graphical Models
Bayesian Learning in Undirected Graphical Models Zoubin Ghahramani Gatsby Computational Neuroscience Unit University College London, UK http://www.gatsby.ucl.ac.uk/ Work with: Iain Murray and Hyun-Chul
More information11280 Electrical Resistivity Tomography Time-lapse Monitoring of Three-dimensional Synthetic Tracer Test Experiments
11280 Electrical Resistivity Tomography Time-lapse Monitoring of Three-dimensional Synthetic Tracer Test Experiments M. Camporese (University of Padova), G. Cassiani* (University of Padova), R. Deiana
More informationAcceptable Ergodic Fluctuations and Simulation of Skewed Distributions
Acceptable Ergodic Fluctuations and Simulation of Skewed Distributions Oy Leuangthong, Jason McLennan and Clayton V. Deutsch Centre for Computational Geostatistics Department of Civil & Environmental Engineering
More informationBayesian reservoir characterization
Bayesian reservoir characterization Bayesian reservoir characterization Luiz Lucchesi Loures ABSTRACT This research aims to provide a complete solution for reservoir properties determination, including
More informationComputing Maximum Entropy Densities: A Hybrid Approach
Computing Maximum Entropy Densities: A Hybrid Approach Badong Chen Department of Precision Instruments and Mechanology Tsinghua University Beijing, 84, P. R. China Jinchun Hu Department of Precision Instruments
More informationF denotes cumulative density. denotes probability density function; (.)
BAYESIAN ANALYSIS: FOREWORDS Notation. System means the real thing and a model is an assumed mathematical form for the system.. he probability model class M contains the set of the all admissible models
More informationIJMGE Int. J. Min. & Geo-Eng. Vol.49, No.1, June 2015, pp
IJMGE Int. J. Min. & Geo-Eng. Vol.49, No.1, June 2015, pp.131-142 Joint Bayesian Stochastic Inversion of Well Logs and Seismic Data for Volumetric Uncertainty Analysis Moslem Moradi 1, Omid Asghari 1,
More informationCOS513 LECTURE 8 STATISTICAL CONCEPTS
COS513 LECTURE 8 STATISTICAL CONCEPTS NIKOLAI SLAVOV AND ANKUR PARIKH 1. MAKING MEANINGFUL STATEMENTS FROM JOINT PROBABILITY DISTRIBUTIONS. A graphical model (GM) represents a family of probability distributions
More informationChapter 17: Undirected Graphical Models
Chapter 17: Undirected Graphical Models The Elements of Statistical Learning Biaobin Jiang Department of Biological Sciences Purdue University bjiang@purdue.edu October 30, 2014 Biaobin Jiang (Purdue)
More informationUSING GEOSTATISTICS TO DESCRIBE COMPLEX A PRIORI INFORMATION FOR INVERSE PROBLEMS THOMAS M. HANSEN 1,2, KLAUS MOSEGAARD 2 and KNUD S.
USING GEOSTATISTICS TO DESCRIBE COMPLEX A PRIORI INFORMATION FOR INVERSE PROBLEMS THOMAS M. HANSEN 1,2, KLAUS MOSEGAARD 2 and KNUD S. CORDUA 1 1 Institute of Geography & Geology, University of Copenhagen,
More informationGEOINFORMATICS Vol. II - Stochastic Modelling of Spatio-Temporal Phenomena in Earth Sciences - Soares, A.
STOCHASTIC MODELLING OF SPATIOTEMPORAL PHENOMENA IN EARTH SCIENCES Soares, A. CMRP Instituto Superior Técnico, University of Lisbon. Portugal Keywords Spacetime models, geostatistics, stochastic simulation
More informationStatistics and Data Analysis
Statistics and Data Analysis The Crash Course Physics 226, Fall 2013 "There are three kinds of lies: lies, damned lies, and statistics. Mark Twain, allegedly after Benjamin Disraeli Statistics and Data
More informationFREQUENTIST BEHAVIOR OF FORMAL BAYESIAN INFERENCE
FREQUENTIST BEHAVIOR OF FORMAL BAYESIAN INFERENCE Donald A. Pierce Oregon State Univ (Emeritus), RERF Hiroshima (Retired), Oregon Health Sciences Univ (Adjunct) Ruggero Bellio Univ of Udine For Perugia
More informationThe Stochastic Inverse Problem in Groundwater Hydrology: A Kalman Filter Approach
The Stochastic Inverse Problem in Groundwater Hydrology: A Kalman Filter Approach Franz Konecny Institute of Mathematics and Applied Statistics C/muerszZ?/ o/ Agrzcw/Zwm/ ^czences, Fzermo /W J#0 Ore^or
More informationGeophysical Data Analysis: Discrete Inverse Theory
Geophysical Data Analysis: Discrete Inverse Theory MATLAB Edition William Menke Lamont-Doherty Earth Observatory and Department of Earth and Environmental Sciences Columbia University. ' - Palisades, New
More informationBetter Simulation Metamodeling: The Why, What and How of Stochastic Kriging
Better Simulation Metamodeling: The Why, What and How of Stochastic Kriging Jeremy Staum Collaborators: Bruce Ankenman, Barry Nelson Evren Baysal, Ming Liu, Wei Xie supported by the NSF under Grant No.
More informationPreliminary Statistics Lecture 2: Probability Theory (Outline) prelimsoas.webs.com
1 School of Oriental and African Studies September 2015 Department of Economics Preliminary Statistics Lecture 2: Probability Theory (Outline) prelimsoas.webs.com Gujarati D. Basic Econometrics, Appendix
More informationPhysics 403. Segev BenZvi. Choosing Priors and the Principle of Maximum Entropy. Department of Physics and Astronomy University of Rochester
Physics 403 Choosing Priors and the Principle of Maximum Entropy Segev BenZvi Department of Physics and Astronomy University of Rochester Table of Contents 1 Review of Last Class Odds Ratio Occam Factors
More informationHomework Week 8 G = H T S. Given that G = H T S, using the first and second laws we can write,
Statistical Molecular hermodynamics University of Minnesota Homework Week 8 1. By comparing the formal derivative of G with the derivative obtained taking account of the first and second laws, use Maxwell
More informationPhysician Performance Assessment / Spatial Inference of Pollutant Concentrations
Physician Performance Assessment / Spatial Inference of Pollutant Concentrations Dawn Woodard Operations Research & Information Engineering Cornell University Johns Hopkins Dept. of Biostatistics, April
More informationQuantitative Trendspotting. Rex Yuxing Du and Wagner A. Kamakura. Web Appendix A Inferring and Projecting the Latent Dynamic Factors
1 Quantitative Trendspotting Rex Yuxing Du and Wagner A. Kamakura Web Appendix A Inferring and Projecting the Latent Dynamic Factors The procedure for inferring the latent state variables (i.e., [ ] ),
More informationModeling Real Estate Data using Quantile Regression
Modeling Real Estate Data using Semiparametric Quantile Regression Department of Statistics University of Innsbruck September 9th, 2011 Overview 1 Application: 2 3 4 Hedonic regression data for house prices
More informationErgodicity in data assimilation methods
Ergodicity in data assimilation methods David Kelly Andy Majda Xin Tong Courant Institute New York University New York NY www.dtbkelly.com April 15, 2016 ETH Zurich David Kelly (CIMS) Data assimilation
More informationHMM part 1. Dr Philip Jackson
Centre for Vision Speech & Signal Processing University of Surrey, Guildford GU2 7XH. HMM part 1 Dr Philip Jackson Probability fundamentals Markov models State topology diagrams Hidden Markov models -
More informationCovariance function estimation in Gaussian process regression
Covariance function estimation in Gaussian process regression François Bachoc Department of Statistics and Operations Research, University of Vienna WU Research Seminar - May 2015 François Bachoc Gaussian
More informationAdjustments for the Display of Quantized Ion Channel Dwell Times in Histograms with Logarithmic Bins
662 Biophysical Journal Volume 78 February 2000 662 667 Adjustments for the Display of Quantized Ion Channel Dwell Times in Histograms with Logarithmic Bins J. Alex Stark and Stephen B. Hladky National
More informationBAYESIAN KRIGING AND BAYESIAN NETWORK DESIGN
BAYESIAN KRIGING AND BAYESIAN NETWORK DESIGN Richard L. Smith Department of Statistics and Operations Research University of North Carolina Chapel Hill, N.C., U.S.A. J. Stuart Hunter Lecture TIES 2004
More informationModeling Seepage Control in Hydraulic Structures
Number 3 Volume 13 September2006 Journal of Engineering Dr. Rafa H. El-Suhaili Professor Environmental Engineering Department College of Engineering Fawaz M. Aziz Al-Joubori M.Sc Student Irrigation and
More informationSoil Uncertainty and Seismic Ground Motion
Boris Jeremić and Kallol Sett Department of Civil and Environmental Engineering University of California, Davis GEESD IV May 2008 Support by NSF CMMI 0600766 Outline Motivation Soils and Uncertainty Stochastic
More informationLightweight Probabilistic Deep Networks Supplemental Material
Lightweight Probabilistic Deep Networks Supplemental Material Jochen Gast Stefan Roth Department of Computer Science, TU Darmstadt In this supplemental material we derive the recipe to create uncertainty
More informationINTERPOLATION AND UPDATE IN DYNAMIC DATA-DRIVEN APPLICATION SIMULATIONS
INTERPOLATION AND UPDATE IN DYNAMIC DATA-DRIVEN APPLICATION SIMULATIONS CRAIG C. DOUGLAS University of Kentucky, Department of Computer Science, 325 McVey Hall, Lexington, KY 4561 and Yale University,
More information