Show simple item record

dc.contributor.authorVillagran-Hernandez, Alejandro
dc.date.accessioned2009-07-09T21:29:03Z
dc.date.available2009-07-09T21:29:03Z
dc.date.issued2009-07-09T21:29:03Z
dc.date.submittedMay 2009
dc.identifier.urihttp://hdl.handle.net/1928/9356
dc.description.abstractIntensive computational methods have been used by Earth scientists in a wide range of problems in data inversion and uncertainty quantification such as earthquake epicenter location and climate projections. To quantify the uncertainties resulting from a range of plausible model configurations it is necessary to estimate a multidimensional probability distribution. The computational cost of estimating these distributions for geoscience applications is impractical using traditional methods such as Metropolis/Gibbs algorithms as simulation costs limit the number of experiments that can be obtained reasonably. Several alternate sampling strategies have been proposed that could improve on the sampling efficiency including Multiple Very Fast Simulated Annealing (MVFSA) and Adaptive Metropolis algorithms. As a goal of this research, the performance of these proposed sampling strategies are evaluated with a surrogate climate model that is able to approximate the noise and response behavior of a realistic atmospheric general circulation model (AGCM). The surrogate model is fast enough that its evaluation can be embedded in these Monte Carlo algorithms. The goal of this thesis is to show that adaptive methods can be superior to MVFSA to approximate the known posterior distribution with fewer forward evaluations. However, the adaptive methods can also be limited by inadequate sample mixing. The Single Component and Delayed Rejection Adaptive Metropolis algorithms were found to resolve these limitations, although challenges remain to approximating multi-modal distributions. The results show that these advanced methods of statistical inference can provide practical solutions to the climate model calibration problem and challenges in quantifying climate projection uncertainties. The computational methods would also be useful to problems outside climate prediction, particularly those where sampling is limited by availability of computational resources.en_US
dc.description.sponsorshipNational Science Foundation, grant OCE-0415251. Consejo Nacional de Ciencia y Tecnologia, Mexico, grant 159764.en_US
dc.language.isoen_USen_US
dc.subjectMonte Carloen_US
dc.subjectClimate Modelsen_US
dc.subjectAdaptive Metropolisen_US
dc.subjectSimulated Annealingen_US
dc.subjectModel Calibrationen_US
dc.subjectBayesian Inferenceen_US
dc.subject.lcshClimatology--Statistical methods.
dc.subject.lcshGlobal warming--Mathematical models.
dc.subject.lcshMonte Carlo method.
dc.subject.lcshSimulated annealing (Mathematics)
dc.titleMonte Carlo strategies for calibration in climate modelsen_US
dc.typeDissertationen_US
dc.description.degreeDoctor of Philosophy Statisticsen_US
dc.description.levelDoctoralen_US
dc.description.departmentUniversity of New Mexico. Dept. of Mathematics and Statisticsen_US
dc.description.advisorHuerta, Gabriel
dc.description.committee-memberBedrick, Edward
dc.description.committee-memberGuindani, Michele
dc.description.committee-memberGalewsky, Joseph


Files in this item

Thumbnail

This item appears in the following Collection(s)

Show simple item record