CMG Research: Boundary Inverse Problems in Glaciology Final Report
|
|
- Patricia Hubbard
- 5 years ago
- Views:
Transcription
1 CMG Research: Boundary Inverse Problems in Glaciology Final Report M. Truffer, D. Maxwell, S. Avdonin October 19, Introduction The primary purpose of this proposal was to develop inverse methods to deduce basal boundary conditions from surface velocity observations of glaciers. In previous work we had developed methods for a simple one-dimensional model that could be formulated as a linear inverse problem (Truffer, 24). In this work we proposed to extend these ideas to non-linear problems that arise when more dimensions are considered. The stated goals of the proposal were to derive robust inverse methods for finding basal velocity fields demonstrate the limits of such inverse methods (e.g. resolution) apply these methods to existing glaciological data sets It is the nature and intend of CMG proposals to not only make substantial contributions in the field of geophysics, but to also advance mathematical research. In the sections below, we separate our completed and some ongoing work along the lines of existing and planned publications. The material is organized into sections that range from almost purely geophysical to almost purely mathematical. A separate section addresses educational contributions. 2 Seasonal evolution of the basal boundary We used surface measurements from Taku Glacier, Alaska, to show that the observed seasonal evolution in surface velocity gradients can be explained by changes in the basal hydraulic gradient. This work was not based on any inverse methods per-se, but rather on exploring a greatly reduced parameter space (till friction angle and hydraulic gradient). However, it helped us develop tools for modeling full force balance ice flow (in 2D) with a plastic boundary condition, and showed that this could explain observed flow features (Fig. 1). This work was published in the Journal of Glaciology (Truffer and others, 29). 1
2 Figure 1: Modelling of a longitudinal section at Taku Glacier. The upper panel shows winter conditions with no water drainage and the lower panel shows summer conditions with elevated basal water pressure. The color bar shows velocities in ma 1. 3 Iterative methods for glacier basal boundary conditions Large non-linear problems are often most efficiently addressed via iterative methods. In Maxwell and others (28b) we introduced an ad-hoc accelerated version of an inverse method known as Kozlov-Maz ya iteration (Kozlov and Maz ya, 199). It operates by making an assumption about the basal boundary condition and then iterating between Neumann and Dirichlet conditions at the surface and the base. Both conditions are known at the surface, and the iteration has been shown to converge in the linear case. We have found it to be too slow in its originally proposed form, but introduced an accelerated version that worked well on synthetic examples, as well as applications to glacier. Figure 2 shows the application to Athabasca Glacier, where measurements of both surface velocities and internal deformation could be used to validate the method. We also applied the model to Glaciar Perito Moreno, where suitable ice thickness and surface velocities are known. The iterative method of Maxwell and others (28b) was applied and extended by others in Arthern and Gudmundsson (21). The method of Maxwell and others (28b) is fast, but suffers from spurious oscillations in grounded or near grounded regions of the glacier and hence has difficulty resolving sharp velocity transitions. These deficiencies make it difficult to infer so-called slipperiness coefficients that relate basal velocities and stresses. We have investigated two methods based on sequential quadratic optimization (SQO) that have the potential to do a better job of resolving these features. The first method infers the basal velocities directly (a constrained Dirichlet problem), while the second infers slipperiness coefficients (a constrained Robin problem). We presented a comparison of these methods with our original Kozlov-Maz ya scheme in Maxwell and others (28a) (Fig. 3). The SQO methods (and especially the constrained Robin problem) do a superior job of resolving basal boundary conditions, particularly when the basal velocities are governed by Coulomb friction. These methods are also significantly more expensive than Kozlov-Maz ya iteration. 2
3 Figure 2: Athabasca Glacier: (a) modeled velocity contour lines ( ma 1 ); (b) contour lines derived from measurements (Raymond, 1971); (c) measured (squares) and modeled (dashed curve) surface velocities and measurement-derived (asterisks) and modeled basal velocities (solid curve); and (d) modeled basal shear stress. 4 Stopping criteria for inverting surface velocities An integral part of any robust inverse method is the recognition that all data (and models) are subject to error. The fact that we are solving ill-posed problems means that it is not only unnecessary to fit data exactly, it is actually undesirable. The reason is that fitting data too well results in unrealistic features in the solutions. The general approach is, therefore, to optimize a certain property of the solution to the inverse problem subject to the condition of obtaining a desired fit to input data (tolerance). It can be shown that our iterative methods produce minimally featured (smooth) solutions. We applied these ideas to the problem of inverting velocity fields for basal stickiness, using the Shallow Shelf Approximation, which is a low-order theory for ice flow over a medium that offers little or no resistance (e.g. soft till or water). These equations have been inverted in previous work before, using a steepest descent method, but not much attention has been paid to the level of desired data fit. We used a series of synthetic ice stream models to demonstrate the effect of over- or underfitting data and showed how to find optimal inversion results, by fitting data to within a given tolerance. In this work we discovered that the steepest descent method converges very slowly and will often not reach optimal resolution in reasonable time. To alleviate this problem, we introduced a new, rapidly converging method, which we labeled the Incomplete Gauss-Newton method. It is based on using a series of regular Gauss-Newton minimizations of the linear problem, and has substantially better performance than either steepest descent or nonlinear conjugate gradient methods (Fig. 4). This work was written up by our graduate student, Marijke Habermann, and is now under review at the Journal of Glaciology. 3
4 Velocity.4.2 Velocity.4.2 Velocity a..2 b..2 c Stress.2 Stress.2 Stress.2 d. e. f. Figure 3: Comparison of inverse method performances for a test case involving a plastic boundary condition. The top row shows results for basal velocity and the bottom for basal stress. The control is shown with dashed lines. Kozlov-Maz ya iteration (a,d) introduces spurious oscillations, which are avoided by SQO methods. b and e shows constrained Dirichlet minimization and c and f constrained Robin minimization Figure 4: Comparison of steepest descent, nonlinear conjugate gradient and incomplete Gauss- Newton methods. The incomplete Gauss-Newton method performs much faster than the others, and the steepest descent method does not reach the desired level of misfit T, even after a large number of iterations. 4
5 Figure 5: Finding a derivative from noisy data using the siple library. 5 A public inverse methods library The work described above led to a variety of tools that are applicable to any (linear or nonlinear) inverse problem, and are thus of interest to a broader community. The tools were all developed in freely available Python software. The toolbox Small Inverse Problems (siple) together with an extensive tutorial and quick refresher on inverse problems is available at The tutorial works through the example of finding a derivative of a (inherently) noisy data series, which is also an ill-posed problems. It is used to demonstrate that these tools are not unique to the geophysical problem being solved, but can be applied quite generally (Fig. 5). Incidentally, the problem of finding derivatives from data has a direct applicability when finding velocities from position measurements, such as the high resolution GPS time series. We plan to explore this in future work. 6 Linking to ice sheet models We are currently working on incorporating the inverse tools developed in this grant to UAF s Parallel Ice Sheet Model (PISM, see This work was started under this NSF grant and is now being continued under a follow-up NASA modeling grant to E. Bueler. It is our goal to include parameter estimation as an integral part of the ice sheet model. The model is open source, and has consequently become one of the most used ice sheet models in the world. 7 Contributions to mathematics Several aspects of our work was published or is prepared for publication in the mathematics literature. It is in keeping with the spirit of the CMG program that contributions span from advances in geophysics to mathematics. While developing iterative inverse tools we have discovered a previously overlooked fact that Kozlov-Maz ya iteration can be posed as a more classical inverse method (Landweber iteration) when a suitable perspective is taken. This observation has lead to two insights: 5
6 The body of existing work on Landweber iteration can be applied to Kozlov-Maz ya iteration. In particular, our acceleration method for Kozlov-Maz ya iteration can be understood in terms of the conjugate gradient method for accelerating Landweber iteration. There are three other inverse schemes (and accelerated variants) that are closely related to Kozlov-Maz ya iteration. A literature search shows that two of the related schemes have been discovered, the third is not known, and the four methods are not already known to be related to each other. We have a manuscript in progress (Maxwell et al., in prep.) for the mathematics literature that develops these ideas. This is an interesting mathematical result with direct implications for geophysics. In particular showing the relation between these different methods clears up some of the confusion for the applied geophysicist to decide which inverse method is the correct one. We also made some advances on the purely mathematical side. The main geophysical goal of our project is solving boundary inverse problems for real glaciers. Even a simple mathematical model of such a problem involves a nonlinear partial differential equation (PDE). More sophisticated models are described by systems of several nonlinear PDEs. Therefore, our project gives rise to many difficult mathematical problems, and its mathematical part is closely related to several areas of mathematics such as the theory of inverse problems for partial differential equations, control theory, approximation theory, and signal processing (sampling and interpolation theory). In Avdonin and others (29b); Maxwell and others (28b) we were able to adjust a very efficient Kozlov Maz ya iteration method (which was originally proposed for systems of linear PDEs) to our nonlinear equations. There are other efficient methods for solving inverse problems for linear PDEs which could potentially serve for our purposes. First of all, we consider the Boundary Control (BC) method in inverse theory. This method is based on controllability of the corresponding partial differential equations. Avdonin and others (29a); Avdonin and Mikhaylov (28); Avdonin (28) contain new results on controllability of partial differential equations. In particular, these are results on graphs; graphs are reasonably considered as a good approximation of 2-dimensional and 3-dimensional problems. On the other hand, for numerical realization of solutions to our problems, it important to understand how to represent a solution to a PDE in a form of a series of its samples. Avdonin and Ivanov (28); Avdonin and others (28) present some new results in this direction. 8 Educational activities As a direct result of this grant, we taught a one-semester graduate short course on Inverse Methods in Glaciology at UAF. The course material was developed during a sabbatical stay and was taught at ETH Zürich as well. The material together with Matlab code can be found at truffer/inverse. This same material was also used to teach a module and supervise a project on inverse theory at the McCarthy Summer School of Glaciology in 21, which will be repeated in 212 ( While the siple library (see above) is well-suited to teaching inverse methods (from theory to application), it has not (yet) been used in that fashion. It is publicly available though and perfectly suited for self-study. This grant has also supported a PhD student in Geophysics (Marijke Habermann) and aspects of PhD mathematics students Victor Mikhaylov and Anna Bulanova. 6
7 9 Dissemination of results In addition to the peer-reviewed literature mentioned in the previous sections we presented results at a variety of meetings and seminars. Habermann has presented aspects of her work at several meetings (Northwest Glaciologists 21, 211; AGU 21; SIRG New Zealand, 21; Karthaus, 29) and at seminars at the New Zealand Antarctic Center and the University of Washington. Maxwell has presented a poster (29) and an invited talk (21) at AGU. Truffer presented at a WAIS meeting, as well as at the 28 IGS meeting in Limerick. We also helped organize a first-of-its-kind special session at the 28 AGU meeting, which attracted 16 abstracts. References Arthern, Robert J and G. H. Gudmundsson, 21. Initialization of ice-sheet forecasts viewed as an inverse Robin problem, J. Glaciol., 56(197), Avdonin, S., 28. Control problems on quantum graphs, Analysis on Graphs and Its Applications, Proceedings of Symposia in Pure Mathematics, vol. 77, Avdonin, Sergei, B.P. Belinskiy and S.A. Ivanov, 29a. Exact controllability of an elastic ring, Applied Math. Optim, 6(1), Avdonin, S., A. Bulanova and D. Ovsyannikov, 28. Optimal cubature formulae related to solutions of initial boundary value problems, Vestnik St. Petersburg Univ., (2), Avdonin, S. and S. Ivanov, 28. Sampling and interpolation problems for vector valued signals in the Paley Wiener spaces, IEEE Trans. Signal Proc., 56(11), Avdonin, S., V. Kozlov, D. Maxwell and M. Truffer, 29b. Iterative methods for solving a nonlinear boundary inverse problem in glaciology, J. Inverse and Ill-Posed Problems. Avdonin, S. and V. Mikhaylov, 28. Controllability of partial differential equations on graphs, Appl. Math., 35, Kozlov, V A and V G Maz ya, 199. On iterative procedures for solving ill-posed boundary value problems that preserve differential equations, Lenningr. Math. J., 1, Maxwell, David, Martin Truffer and Sergei Avdonin, 28a. Inverse Methods for Reconstructing Basal Boundary Data, American Geophysical Union 28 Fall Meeting, (poster). Maxwell, David., Martin Truffer, Sergei Avdonin and Martin Stuefer, 28b. Determining glacier velocities and stresses with inverse methods: an iterative scheme, J. Glaciol., 54, Raymond, C. F., Flow in a transverse section of Athabasca Glacier, Alberta, Canada, J. Glaciol., 1, Truffer, M., 24. The basal speed of valley glaciers: an inverse approach, J. Glaciol., 5(169), Truffer, M., R. J. Motyka, Michael Hekkers, I. M. Howat and M. A. King, 29. Terminus dynamics at an advancing glacier: Taku Glacier, Alaska, J. Glaciol., 55(194),
Modeled and observed fast flow in the Greenland ice sheet
Modeled and observed fast flow in the Greenland ice sheet Ed Bueler 1 Constantine Khroulev 1 Andy Aschwanden 2 Ian Joughin 3 1 Dept of Mathematics and Statistics, University of Alaska Fairbanks 2 Arctic
More informationITERATIVE METHODS FOR SOLVING A NONLINEAR BOUNDARY INVERSE PROBLEM IN GLACIOLOGY S. AVDONIN, V. KOZLOV, D. MAXWELL, AND M. TRUFFER
ITERATIVE METHODS FOR SOLVING A NONLINEAR BOUNDARY INVERSE PROBLEM IN GLACIOLOGY S. AVDONIN, V. KOZLOV, D. MAXWELL, AND M. TRUFFER Abstract. We address a Cauchy problem for a nonlinear elliptic PDE arising
More informationSUPPLEMENTARY INFORMATION
SUPPLEMENTARY INFORMATION DOI: 1.138/NGEO1218 Supplementary information Ice speed of a calving glacier modulated by small fluctuations in basal water pressure Shin Sugiyama 1, Pedro Skvarca 2, Nozomu Naito
More informationInverse methods in glaciology
Inverse methods in glaciology Martin Truer University of Alaska Fairbanks McCarthy, Summer School, 2010 1 / 20 General Problem Setting Examples of inverse problems Solution methods 2 / 20 Outline General
More informationAn iterative scheme for determining glacier velocities and stresses
888 Journal of Glaciology, Vol. 54, No. 188, 2008 An iterative scheme for determining glacier velocities and stresses David MAXWELL, 1 Martin TRUFFER, 2 ergei AVDONIN, 1 Martin TUEFER 2 1 Department of
More informationGlacier Hydrology II: Theory and Modeling
Glacier Hydrology II: Theory and Modeling McCarthy Summer School 2018 Matt Hoffman Gwenn Flowers, Simon Fraser Operated by Los Alamos National Security, LLC for the U.S. Department of Energy's NNSA Observations
More informationPISM, a Parallel Ice Sheet Model: Current status of our Antarctic ice sheet simulation
PISM, a Parallel Ice Sheet Model: Current status of our Antarctic ice sheet simulation Craig Lingle, 1 Ed Bueler, 2 Jed Brown, 1 and David Covey, 1 1 Geophysical Institute, Univ. of Alaska, Fairbanks 2
More informationConjugate Directions for Stochastic Gradient Descent
Conjugate Directions for Stochastic Gradient Descent Nicol N Schraudolph Thore Graepel Institute of Computational Science ETH Zürich, Switzerland {schraudo,graepel}@infethzch Abstract The method of conjugate
More informationPISM, a Parallel Ice Sheet Model
PISM, a Parallel Ice Sheet Model Ed Bueler 1, Craig Lingle 2, and Jed Brown 3 1 Department of Mathematics and Statistics, Univ. of Alaska, Fairbanks 2 Geophysical Institute, Univ. of Alaska, Fairbanks
More informationSCIENTIFIC REPORT NERC GEF
SCIENTIFIC REPORT NERC GEF Loan 927 Measuring changes in the dynamics of Pine Island Glacier, Antarctica A.M. Smith & E.C. King, British Antarctic Survey (BAS) pp J.B.T. Scott ABSTRACT A brief period of
More informationIterative Methods for Smooth Objective Functions
Optimization Iterative Methods for Smooth Objective Functions Quadratic Objective Functions Stationary Iterative Methods (first/second order) Steepest Descent Method Landweber/Projected Landweber Methods
More informationMath 411 Preliminaries
Math 411 Preliminaries Provide a list of preliminary vocabulary and concepts Preliminary Basic Netwon s method, Taylor series expansion (for single and multiple variables), Eigenvalue, Eigenvector, Vector
More informationverification, validation, and basal strength in models for the present state of ice sheets
verification, validation, and basal strength in models for the present state of ice sheets Ed Bueler 1 Constantine Khroulev 2 Andy Aschwanden 3 1 Dept of Mathematics and Statistics, University of Alaska
More informationRegistration-guided least-squares waveform inversion
Registration-guided least-squares waveform inversion Hyoungsu Baek 1, Henri Calandra, Laurent Demanet 1 1 MIT Mathematics department, TOTAL S.A. January 15 013 Abstract Full waveform inversion with frequency
More informationGreenland subglacial drainage evolution regulated by weakly-connected regions of the bed
Greenland subglacial drainage evolution regulated by weakly-connected regions of the bed Matthew Hoffman Stephen Price Lauren Andrews Ginny Catania Weakly-connected Drainage Distributed Drainage Channelized
More informationA hybrid Marquardt-Simulated Annealing method for solving the groundwater inverse problem
Calibration and Reliability in Groundwater Modelling (Proceedings of the ModelCARE 99 Conference held at Zurich, Switzerland, September 1999). IAHS Publ. no. 265, 2000. 157 A hybrid Marquardt-Simulated
More informationMulti-Modal Flow in a Thermocoupled Model of the Antarctic Ice Sheet, with Verification
Multi-Modal Flow in a Thermocoupled Model of the Antarctic Ice Sheet, with Verification Craig Lingle 1 Jed Brown 2 Ed Bueler 2 1 Geophysical Institute University of Alaska Fairbanks, USA 2 Department of
More informationIPAM Summer School Optimization methods for machine learning. Jorge Nocedal
IPAM Summer School 2012 Tutorial on Optimization methods for machine learning Jorge Nocedal Northwestern University Overview 1. We discuss some characteristics of optimization problems arising in deep
More informationRESEARCH ARTICLE. Fast reconstruction of harmonic functions from Cauchy data using the Dirichlet-to-Neumann map and integral equations
Inverse Problems in Science and Engineering Vol. 00, No. 00, September 2010, 1 10 RESEARCH ARTICLE Fast reconstruction of harmonic functions from Cauchy data using the Dirichlet-to-Neumann map and integral
More informationA posteriori estimator for the accuracy of the shallow shelf approximation
Math Geosci manuscript No. will be inserted by the editor A posteriori estimator for the accuracy of the shallow shelf approximation Guillaume Jouvet Marco Picasso Received: date / Accepted: date Abstract
More informationSeismic tomography with co-located soft data
Seismic tomography with co-located soft data Mohammad Maysami and Robert G. Clapp ABSTRACT There is a wide range of uncertainties present in seismic data. Limited subsurface illumination is also common,
More informationMEANDER CURVE (MODIFIED FOR ADEED)
MEANDER CURVE (MODIFIED FOR ADEED) Overview: Friction between water and stream banks causes water to move in a corkscrew fashion. This helical flow is called the water spiral. Gravity and the water spiral
More informationA multigrid method for large scale inverse problems
A multigrid method for large scale inverse problems Eldad Haber Dept. of Computer Science, Dept. of Earth and Ocean Science University of British Columbia haber@cs.ubc.ca July 4, 2003 E.Haber: Multigrid
More informationNumerical Methods I Solving Nonlinear Equations
Numerical Methods I Solving Nonlinear Equations Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 MATH-GA 2011.003 / CSCI-GA 2945.003, Fall 2014 October 16th, 2014 A. Donev (Courant Institute)
More information3 Erosion and Deposition by Ice
CHAPTER 12 3 Erosion and Deposition by Ice SECTION Agents of Erosion and Deposition BEFORE YOU READ After you read this section, you should be able to answer these questions: What are glaciers? How do
More informationAppendix A. Supplementary Material
73 74 75 76 77 78 79 71 711 712 713 714 715 716 717 718 719 72 721 722 723 724 725 726 727 728 Appendix A. Supplementary Material Antarctic regions We initially localize over all of Antarctica, and analyze
More informationInverse Problems and Optimal Design in Electricity and Magnetism
Inverse Problems and Optimal Design in Electricity and Magnetism P. Neittaanmäki Department of Mathematics, University of Jyväskylä M. Rudnicki Institute of Electrical Engineering, Warsaw and A. Savini
More informationwarwickphysics Physics Courses
warwickphysics Physics Courses Entry October 2015 2 The University of Warwick Department of Physics 3 Introduction The physics course is designed as a broad and flexible education. The department offers
More informationLinear Inverse Problems. A MATLAB Tutorial Presented by Johnny Samuels
Linear Inverse Problems A MATLAB Tutorial Presented by Johnny Samuels What do we want to do? We want to develop a method to determine the best fit to a set of data: e.g. The Plan Review pertinent linear
More informationCurriculum Vitae of Andy Aschwanden, PhD
Curriculum Vitae of Andy Aschwanden, PhD Current work address: Geophysical Institute, 903 Koyukuk Dr., Fairbanks AK 99775. Email: aaschwanden@alaska.edu Research Interests dynamics and thermodynamics of
More information5. The topographic map below shows a lake and two rivers.
Mapping A B1 1. The diagram below shows latitude measurements every 10 degrees and longitude measurements every 15 degrees. What is the latitude and longitude of point X? 5. The topographic map below shows
More informationGenerating Equidistributed Meshes in 2D via Domain Decomposition
Generating Equidistributed Meshes in 2D via Domain Decomposition Ronald D. Haynes and Alexander J. M. Howse 2 Introduction There are many occasions when the use of a uniform spatial grid would be prohibitively
More informationSolving PDEs with Multigrid Methods p.1
Solving PDEs with Multigrid Methods Scott MacLachlan maclachl@colorado.edu Department of Applied Mathematics, University of Colorado at Boulder Solving PDEs with Multigrid Methods p.1 Support and Collaboration
More informationCool Off, Will Ya! Investigating Effect of Temperature Differences between Water and Environment on Cooling Rate of Water
Ding 1 Cool Off, Will Ya! Investigating Effect of Temperature Differences between Water and Environment on Cooling Rate of Water Chunyang Ding 000844-0029 Physics HL Ms. Dossett 10 February 2014 Ding 2
More informationModelling of surface to basal hydrology across the Russell Glacier Catchment
Modelling of surface to basal hydrology across the Russell Glacier Catchment Sam GAP Modelling Workshop, Toronto November 2010 Collaborators Alun Hubbard Centre for Glaciology Institute of Geography and
More informationB.A. (Physics) May 1999 Whitman College, Walla Walla, Washington
Julie Elliott Department of Earth and Atmospheric Sciences, Cornell University, Ithaca, NY 14853 julie.elliott@cornell.edu www.geo.cornell.edu/research_staff/jle84 Education Ph.D. (Geophysics) August 2011
More informationYouzuo Lin and Lianjie Huang
PROCEEDINGS, Thirty-Ninth Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, February 24-26, 2014 SGP-TR-202 Building Subsurface Velocity Models with Sharp Interfaces
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 6 Optimization Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002. Reproduction permitted
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 6 Optimization Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002. Reproduction permitted
More informationSUPPLEMENTARY INFORMATION
SUPPLEMENTARY INFORMATION DOI: 10.1038/NGEO1887 Diverse calving patterns linked to glacier geometry J. N. Bassis and S. Jacobs 1. Supplementary Figures (a) (b) (c) Supplementary Figure S1 Schematic of
More informationOptimization. Escuela de Ingeniería Informática de Oviedo. (Dpto. de Matemáticas-UniOvi) Numerical Computation Optimization 1 / 30
Optimization Escuela de Ingeniería Informática de Oviedo (Dpto. de Matemáticas-UniOvi) Numerical Computation Optimization 1 / 30 Unconstrained optimization Outline 1 Unconstrained optimization 2 Constrained
More informationPreconditioning. Noisy, Ill-Conditioned Linear Systems
Preconditioning Noisy, Ill-Conditioned Linear Systems James G. Nagy Emory University Atlanta, GA Outline 1. The Basic Problem 2. Regularization / Iterative Methods 3. Preconditioning 4. Example: Image
More informationGEOPHYSICAL INVERSE THEORY AND REGULARIZATION PROBLEMS
Methods in Geochemistry and Geophysics, 36 GEOPHYSICAL INVERSE THEORY AND REGULARIZATION PROBLEMS Michael S. ZHDANOV University of Utah Salt Lake City UTAH, U.S.A. 2OO2 ELSEVIER Amsterdam - Boston - London
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 5 Nonlinear Equations Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002. Reproduction
More informationSummer School in Glaciology, Fairbanks/McCarthy, Exercises: Glacial geology
Bob Anderson Summer School in Glaciology, Fairbanks/McCarthy, 2010 Exercises: Glacial geology 1. Glacier thickness. We wish to estimate the local thickness of a glacier given only a topographic map of
More informationMerrily we roll along
Merrily we roll along Name Period Date Lab partners Overview Measuring motion of freely falling objects is difficult because they acclerate so fast. The speed increases by 9.8 m/s every second, so Galileo
More informationUnconstrained Multivariate Optimization
Unconstrained Multivariate Optimization Multivariate optimization means optimization of a scalar function of a several variables: and has the general form: y = () min ( ) where () is a nonlinear scalar-valued
More informationGlaciology Exchange (Glacio-Ex) Norwegian/Canadian/US Partnership Program
Glaciology Exchange (Glacio-Ex) Norwegian/Canadian/US Partnership Program Luke Copland University of Ottawa, Canada Jon Ove Hagen University of Oslo, Norway Kronebreeen, Svalbard. Photo: Monica Sund The
More informationSamira Ardani. Academic Interests
Samira Ardani Ocean Engineering Department 200 Seawolf Parkway, Galveston, TX 77554 Office: PMEC, Room 137 Email: arda12@ tamu.edu Academic Interests - Nearshore Wave Modeling: Transformation of Waves
More informationLecture 35 Minimization and maximization of functions. Powell s method in multidimensions Conjugate gradient method. Annealing methods.
Lecture 35 Minimization and maximization of functions Powell s method in multidimensions Conjugate gradient method. Annealing methods. We know how to minimize functions in one dimension. If we start at
More informationMANY methods have been developed so far for solving
IEEE TRANSACTIONS ON NEURAL NETWORKS, VOL. 9, NO. 5, SEPTEMBER 1998 987 Artificial Neural Networks for Solving Ordinary Partial Differential Equations Isaac Elias Lagaris, Aristidis Likas, Member, IEEE,
More informationComparison of an Analytical Method and Matlab to Model Electromagnetic Distribution in a Trough
Comparison of an Analytical Method and Matlab to Model Electromagnetic Distribution in a Trough JJ Bruyns University of Johannesburg, South Africa Jacob@uj.ac.za JC Greeff Tshwane University of Technology
More informationDELFT UNIVERSITY OF TECHNOLOGY
DELFT UNIVERSITY OF TECHNOLOGY REPORT 16-02 The Induced Dimension Reduction method applied to convection-diffusion-reaction problems R. Astudillo and M. B. van Gijzen ISSN 1389-6520 Reports of the Delft
More informationCollaborative Research: Norwegian-United States IPY Scientific Traverse: Climate Variability and Glaciology in East Antarctica
The University of Maine DigitalCommons@UMaine University of Maine Office of Research and Sponsored Programs: Grant Reports Special Collections 6-27-2012 Collaborative Research: Norwegian-United States
More informationGlacier Thermodynamics: Ice Temperature and Heat Transfer Processes
Glacier Thermodynamics: Ice Temperature and Heat Transfer Processes ESS431: Principles of Glaciology ESS505: The Cryosphere Wednesday, 10/24 Ben Hills Today s Objectives: Why do we care about ice temperature?
More informationGEO-SLOPE International Ltd, Calgary, Alberta, Canada Wick Drain
1 Introduction Wick Drain This example is about modeling the behavior of a wick drain. The primary purpose here is to illustrate how interface elements can conveniently be used to include the effects of
More informationAlternative numerical method in continuum mechanics COMPUTATIONAL MULTISCALE. University of Liège Aerospace & Mechanical Engineering
University of Liège Aerospace & Mechanical Engineering Alternative numerical method in continuum mechanics COMPUTATIONAL MULTISCALE Van Dung NGUYEN Innocent NIYONZIMA Aerospace & Mechanical engineering
More informationBureau of Economic Geology, The University of Texas at Austin, Austin, TX. Research Associate, Geophysics from 11/2004
Paul C. Sava Bureau of Economic Geology, The University of Texas at Austin, University Station, Box X, Austin, TX 78758, (512) 471-0327 paul.sava@beg.utexas.edu, http://sepwww.stanford.edu/sep/paul Research
More informationDynamics of Glaciers
Dynamics of Glaciers McCarthy Summer School 01 Andy Aschwanden Arctic Region Supercomputing Center University of Alaska Fairbanks, USA June 01 Note: This script is largely based on the Physics of Glaciers
More informationComparison between least-squares reverse time migration and full-waveform inversion
Comparison between least-squares reverse time migration and full-waveform inversion Lei Yang, Daniel O. Trad and Wenyong Pan Summary The inverse problem in exploration geophysics usually consists of two
More informationA finite element solver for ice sheet dynamics to be integrated with MPAS
A finite element solver for ice sheet dynamics to be integrated with MPAS Mauro Perego in collaboration with FSU, ORNL, LANL, Sandia February 6, CESM LIWG Meeting, Boulder (CO), 202 Outline Introduction
More informationGeochemical and Mineralogical Dispersion Models in Till: Physical Process Constraints and Impacts on Geochemical Exploration Interpretation
Geochemical and Mineralogical Dispersion Models in Till: Physical Process Constraints and Impacts on Geochemical Exploration Interpretation Cliff Stanley Dept. of Geology Acadia University Wolfville,,
More informationPreconditioning. Noisy, Ill-Conditioned Linear Systems
Preconditioning Noisy, Ill-Conditioned Linear Systems James G. Nagy Emory University Atlanta, GA Outline 1. The Basic Problem 2. Regularization / Iterative Methods 3. Preconditioning 4. Example: Image
More informationEducation. B.A. (Physics) May 1999 Whitman College, Walla Walla, Washington. Professional Experience
Julie Elliott Earth, Atmospheric, and Planetary Sciences, Purdue University West Lafayette, IN 47906 765-494-0484 (work) 509-961-3568 (cell) julieelliott@purdue.edu Education Ph.D. (Geophysics) August
More informationLecture 5: Gradient Descent. 5.1 Unconstrained minimization problems and Gradient descent
10-725/36-725: Convex Optimization Spring 2015 Lecturer: Ryan Tibshirani Lecture 5: Gradient Descent Scribes: Loc Do,2,3 Disclaimer: These notes have not been subjected to the usual scrutiny reserved for
More informationOptimization. Totally not complete this is...don't use it yet...
Optimization Totally not complete this is...don't use it yet... Bisection? Doing a root method is akin to doing a optimization method, but bi-section would not be an effective method - can detect sign
More informationOptimization with COMSOL Multiphysics COMSOL Tokyo Conference Walter Frei, PhD Applications Engineer
Optimization with COMSOL Multiphysics COMSOL Tokyo Conference 2014 Walter Frei, PhD Applications Engineer Product Suite COMSOL 5.0 Agenda An introduction to optimization A lot of concepts, and a little
More informationNotes for CS542G (Iterative Solvers for Linear Systems)
Notes for CS542G (Iterative Solvers for Linear Systems) Robert Bridson November 20, 2007 1 The Basics We re now looking at efficient ways to solve the linear system of equations Ax = b where in this course,
More informationLecture Notes to Accompany. Scientific Computing An Introductory Survey. by Michael T. Heath. Chapter 5. Nonlinear Equations
Lecture Notes to Accompany Scientific Computing An Introductory Survey Second Edition by Michael T Heath Chapter 5 Nonlinear Equations Copyright c 2001 Reproduction permitted only for noncommercial, educational
More informationT. Perron Glaciers 1. Glaciers
T. Perron 12.001 Glaciers 1 Glaciers I. Why study glaciers? [PPT: Perito Moreno glacier, Argentina] Role in freshwater budget o Fraction of earth s water that is fresh (non-saline): 3% o Fraction of earth
More informationSEISMIC RESPONSE OF SINGLE DEGREE OF FREEDOM STRUCTURAL FUSE SYSTEMS
3 th World Conference on Earthquake Engineering Vancouver, B.C., Canada August -6, 4 Paper No. 377 SEISMIC RESPONSE OF SINGLE DEGREE OF FREEDOM STRUCTURAL FUSE SYSTEMS Ramiro VARGAS and Michel BRUNEAU
More informationIterative Methods for Solving A x = b
Iterative Methods for Solving A x = b A good (free) online source for iterative methods for solving A x = b is given in the description of a set of iterative solvers called templates found at netlib: http
More informationTHE SIMPLE PENDULUM: A CLASSIC CASE STUDY Originally developed by Dr. Larry Silverstein
THE SIMPLE PENDULUM: A CLASSIC CASE STUDY Originally developed by Dr. Larry Silverstein This Experiment/Activity draws together, and exemplifies, many threads and activities of science, among them: 1.
More information"Ice Sheets and Sea Level Rise: How Should IPCC Handle Deep Uncertainty?" Michael Oppenheimer For Inside the IPCC Princeton University 1 April 2008
"Ice Sheets and Sea Level Rise: How Should IPCC Handle Deep Uncertainty?" Michael Oppenheimer For Inside the IPCC Princeton University 1 April 2008 This Talk is about: IPCCs (controversial) assessment
More informationESS 431 Principles of Glaciology ESS 505 The Cryosphere
MID-TERM November 9, 2015 ESS 431 Principles of Glaciology ESS 505 The Cryosphere Instructions: Please answer the following 5 questions. [The actual 5 questions will be selected from these 12 questions
More informationSynoptic Meteorology I: Other Force Balances
Synoptic Meteorology I: Other Force Balances For Further Reading Section.1.3 of Mid-Latitude Atmospheric Dynamics by J. Martin provides a discussion of the frictional force and considerations related to
More informationHow Do You Teach Vibrations to Technology Students?
Paper ID #9337 How Do You Teach Vibrations to Technology Students? Dr. M. Austin Creasy, Purdue University (Statewide Technology) Assistant Professor Mechanical Engineering Technology Purdue University
More informationSPARSE signal representations have gained popularity in recent
6958 IEEE TRANSACTIONS ON INFORMATION THEORY, VOL. 57, NO. 10, OCTOBER 2011 Blind Compressed Sensing Sivan Gleichman and Yonina C. Eldar, Senior Member, IEEE Abstract The fundamental principle underlying
More informationArctic Science & Engineering. Martin Jeffries. Ron Liston Seminar, 17 October PhD. MSc. Calgary ( ) UK ( )
Arctic Science & Calgary (1981-1985) Engineering PhD Martin Jeffries MSc Ron Liston Seminar, 17 October 2018. UK (1979-1981) Outline Career Background Research Background The State of the Arctic Arctic
More informationBasal Processes (i.e. boundary conditions)
Basal Processes (i.e. boundary conditions) Alan Rempel, University of Oregon Motivation: a) the sliding law b) weathering & erosion Outline: - temperate sliding - drainage elements - intemperate effects
More informationNumerical Optimization Prof. Shirish K. Shevade Department of Computer Science and Automation Indian Institute of Science, Bangalore
Numerical Optimization Prof. Shirish K. Shevade Department of Computer Science and Automation Indian Institute of Science, Bangalore Lecture - 13 Steepest Descent Method Hello, welcome back to this series
More informationNumerical solutions of nonlinear systems of equations
Numerical solutions of nonlinear systems of equations Tsung-Ming Huang Department of Mathematics National Taiwan Normal University, Taiwan E-mail: min@math.ntnu.edu.tw August 28, 2011 Outline 1 Fixed points
More informationStable Adaptive Momentum for Rapid Online Learning in Nonlinear Systems
Stable Adaptive Momentum for Rapid Online Learning in Nonlinear Systems Thore Graepel and Nicol N. Schraudolph Institute of Computational Science ETH Zürich, Switzerland {graepel,schraudo}@inf.ethz.ch
More informationTopic 8c Multi Variable Optimization
Course Instructor Dr. Raymond C. Rumpf Office: A 337 Phone: (915) 747 6958 E Mail: rcrumpf@utep.edu Topic 8c Multi Variable Optimization EE 4386/5301 Computational Methods in EE Outline Mathematical Preliminaries
More informationSPS8. STUDENTS WILL DETERMINE RELATIONSHIPS AMONG FORCE, MASS, AND MOTION.
MOTION & FORCES SPS8. STUDENTS WILL DETERMINE RELATIONSHIPS AMONG FORCE, MASS, AND MOTION. A. CALCULATE VELOCITY AND ACCELERATION. B. APPLY NEWTON S THREE LAWS TO EVERYDAY SITUATIONS BY EXPLAINING THE
More informationPart 4: IIR Filters Optimization Approach. Tutorial ISCAS 2007
Part 4: IIR Filters Optimization Approach Tutorial ISCAS 2007 Copyright 2007 Andreas Antoniou Victoria, BC, Canada Email: aantoniou@ieee.org July 24, 2007 Frame # 1 Slide # 1 A. Antoniou Part4: IIR Filters
More informationDiscontinuous Galerkin methods for nonlinear elasticity
Discontinuous Galerkin methods for nonlinear elasticity Preprint submitted to lsevier Science 8 January 2008 The goal of this paper is to introduce Discontinuous Galerkin (DG) methods for nonlinear elasticity
More informationScientific Computing: Optimization
Scientific Computing: Optimization Aleksandar Donev Courant Institute, NYU 1 donev@courant.nyu.edu 1 Course MATH-GA.2043 or CSCI-GA.2112, Spring 2012 March 8th, 2011 A. Donev (Courant Institute) Lecture
More informationHybrid particle swarm algorithm for solving nonlinear constraint. optimization problem [5].
Hybrid particle swarm algorithm for solving nonlinear constraint optimization problems BINGQIN QIAO, XIAOMING CHANG Computers and Software College Taiyuan University of Technology Department of Economic
More informationNonlinear Optimization for Optimal Control
Nonlinear Optimization for Optimal Control Pieter Abbeel UC Berkeley EECS Many slides and figures adapted from Stephen Boyd [optional] Boyd and Vandenberghe, Convex Optimization, Chapters 9 11 [optional]
More informationVlasov simulations of electron holes driven by particle distributions from PIC reconnection simulations with a guide field
GEOPHYSICAL RESEARCH LETTERS, VOL. 35, L22109, doi:10.1029/2008gl035608, 2008 Vlasov simulations of electron holes driven by particle distributions from PIC reconnection simulations with a guide field
More informationDesign and Optimization of Energy Systems Prof. C. Balaji Department of Mechanical Engineering Indian Institute of Technology, Madras
Design and Optimization of Energy Systems Prof. C. Balaji Department of Mechanical Engineering Indian Institute of Technology, Madras Lecture - 09 Newton-Raphson Method Contd We will continue with our
More informationUnit 15 LESSON 1 WHAT ARE FORCES?
Unit 15 LESSON 1 WHAT ARE FORCES? Pushing and Pulling A force is a push or pull. Forces can cause an object at rest to move, speed up, slow down, change direction, or stop. Forces can change the shape
More informationChapter 3 Numerical Methods
Chapter 3 Numerical Methods Part 2 3.2 Systems of Equations 3.3 Nonlinear and Constrained Optimization 1 Outline 3.2 Systems of Equations 3.3 Nonlinear and Constrained Optimization Summary 2 Outline 3.2
More informationIsaac Newton was a British scientist whose accomplishments
E8 Newton s Laws of Motion R EA D I N G Isaac Newton was a British scientist whose accomplishments included important discoveries about light, motion, and gravity. You may have heard the legend about how
More informationCOURSE DESCRIPTIONS. 1 of 5 8/21/2008 3:15 PM. (S) = Spring and (F) = Fall. All courses are 3 semester hours, unless otherwise noted.
1 of 5 8/21/2008 3:15 PM COURSE DESCRIPTIONS (S) = Spring and (F) = Fall All courses are 3 semester hours, unless otherwise noted. INTRODUCTORY COURSES: CAAM 210 (BOTH) INTRODUCTION TO ENGINEERING COMPUTATION
More informationCOMP 551 Applied Machine Learning Lecture 3: Linear regression (cont d)
COMP 551 Applied Machine Learning Lecture 3: Linear regression (cont d) Instructor: Herke van Hoof (herke.vanhoof@mail.mcgill.ca) Slides mostly by: Class web page: www.cs.mcgill.ca/~hvanho2/comp551 Unless
More informationThe speed of Shor s R-algorithm
IMA Journal of Numerical Analysis 2008) 28, 711 720 doi:10.1093/imanum/drn008 Advance Access publication on September 12, 2008 The speed of Shor s R-algorithm J. V. BURKE Department of Mathematics, University
More informationNotes on Some Methods for Solving Linear Systems
Notes on Some Methods for Solving Linear Systems Dianne P. O Leary, 1983 and 1999 and 2007 September 25, 2007 When the matrix A is symmetric and positive definite, we have a whole new class of algorithms
More informationChapter 2 Finite Element Formulations
Chapter 2 Finite Element Formulations The governing equations for problems solved by the finite element method are typically formulated by partial differential equations in their original form. These are
More information