Inverse Problems, Information Content and Stratospheric Aerosols
|
|
- Anastasia Brown
- 5 years ago
- Views:
Transcription
1 Inverse Problems, Information Content and Stratospheric Aerosols May 22, 214
2 Table of contents 1 Inverse Problems 2 Aerosols 3 OSIRIS 4 Information Content 5 Non-Linear Techniques 6 Results
3 Inverse Problems - Overview Given the output of a system, determine the system s state Often more difficult than the forward problem First described mathematically by Victor Ambartsumian in 1929 but didn t take off until the 5 s Most people, if you describe a train of events to them, will tell you what you the result would be. They can put those events together in their minds, and argue from them that something will come to pass. There are few people, however, who, if you told them the result would be able to evolve from their own inner consciousness what the steps were which led up to that result. This power is what I mean when I talk of reasoning backwards, or analytically. -Sherlock Holmes 1887
4 Inverse Problems - Overview Show up all the time in geophysics, but pretty much every field has them Often ill-posed or underdetermined, although can be over-determined as well. I ll talk about some methods I ve used for investigating and solving the ill-posed/underdetermined variety First, a look at the inverse problem I ve been investigating
5 Aerosols Small droplets (occasionally crystals) of sulphuric acid and water in the stratosphere Reflect solar radiation, decreasing warming (.5 C surface cooling after Pinatubo eruption) Highly variable due to volcanic events, stratospheric circulation patterns IPCC points to aerosols as one of the biggest unknowns in climate forecasting
6 Aerosols To fully describe aerosol we would need to know the particle size distribution and total amount of aerosol. This is pretty much impossible with OSIRIS so we settle for a few simpler parameters Integral properties of the distribution such as total number density, volume, or extinction Assume a distribution shape and try to retrieve parameters of that shape Lognormal is common and well supported by other data ( (ln r ln r g ) 2 ) dn(r) dr = n aer r ln(σ g ) 2π exp 2 ln(σ g ) 2 Leaves only 3 aerosol parameters - mode radius, mode width and number density
7 Aerosols A few typical lognormal distributions
8 OSIRIS OSIRIS is a Canadian instrument launched in 21 Measures radiance along a single line of sight Satellite nods up and down to produce vertical profiles Used to retrieve vertical profiles of ozone, NO 2, BrO and aerosols Solar Illumination Line of Sight Satellite Earth Atmosphere
9 OSIRIS I meas = s1 s J(s) = k sc(s) k(s) k(s)j(s)e τ(s,s 1) ds 4π I diff (s, Ω)p(s, Θ)dΩ Solar Illumination Line of Sight Attenuate J(s) back to satellite Satellite Integrate I diff and rescatter to get J(s) Earth Atmosphere
10 OSIRIS An example of some OSIRIS measurements Ozone visible around 3 and 6 nm NO 2 around 44 nm Broadband aerosol signature >5 nm Spectral Radiance (photons/m 2 /s/nm/st) Wavelength (nm)
11 OSIRIS Goal is to use this broadband aerosol signature to determine size and quantity of aerosol particles A few problems that make this tricky: Finite measurements of continuous distributions (underdetermined system) Aerosol signal is broadband, limiting useful wavelengths (non-orthogonal measurements) Measurements at each height contain information from all heights (non-orthogonal measurements) This altitude coupling also makes the problem non-linear To get an idea of what s possible first take a look at the information content of the measurements
12 Information Content Twomey explored a very nice way of determining information content of linear systems by looking at the eigenfunctions/vectors of the measurements. Say we have some linear system: y = K(s )x(s )ds y = Kx s Then the Jacobian (or kernel), K, is the sensitivity of the measurement, y, to changes in the state, x. K i j = dk i dx j
13 Information Content If we can predict a measurement before its taken, then that measurement provides no new information. This is possible if K is rank deficient, or, K i = a 1 i a j K j (s). Since if this is true we can also reconstruct the measurement y i = a 1 i a j K j x(s )ds = a 1 k a l y l j i s j i l k
14 Information Content For noisy measurements this doesn t need to be strictly true, and if we can reconstruct a measurement within error, then it is redundant. The covariance matrix C ij = K i (s )K j (s )ds s tells us the orthogonality of the measurements. The number of measurements that contain unique information is then given by number of eigenvalues of C that exceed the relative error Λ i Λ 1 > ɛ 2
15 Linearizing the OSIRIS Measurements To apply this to OSIRIS we first have to make the problem linear. Assume changes in aerosol don t change the diffuse field I diff Assume most of the information is coming from the tangent point (lowest altitude we re looking through) Attenuation back to the satellite is negligible This lets us write the radiance in terms of some aerosol properties, and the measurement y as an integral over a particle size distribution. I aer n aer σ aer p aer (Θ) y i = K i (r) dv (r) dr dr
16 Linearizing the OSIRIS Measurements Take 4 measurements, at 2 different wavelengths and 2 different satellite geometries - how much info? K4 K3 K2 K λ = 75 nm Θ = 6 λ = 75 nm Θ = 118 λ = 1.53 µm Θ = 6 λ = 1.53 µm Θ = Radius (µm) Multiple Scatter Single Scatter 4
17 Linearizing the OSIRIS Measurements This gives normalized eigenvalues of 1.,.8,.4,.1 Λi Λ 1 = 1.,.28,.21,.12 Meaning only 1% of the 4 th measurement is new information, resulting in a 1 error magnification. Implying that it is probably not a good idea to extract 4 pieces of information from these 4 measurements. (Unless the measurements were very accurate)
18 Using the Information What can we do with these 3ish pieces of information? It might be possible to define some particle size distribution with 2 parameters and retrieve that, however, just because we have three pieces of information doesn t mean we can retrieve any 3 pieces. The work by Twomey (Also see Backus and Gilbert) provides a method of determining which quantities we re most sensitive to, and a maximum possible error.
19 An Orthogonal Framework Dealing with non-orthogonal functions is hard, so first we make our measurements orthogonal using Principal Component Analysis K T φ T = K T QΛ 1/2 Where Q and Λ are the eigenvectors and values of C = KK T respectively. Arbitrary Units st Component 2 nd Component 3 rd Component 4 th Component Radius (µm) 4
20 An Orthogonal Framework With an orthogonal set it s easy to reconstruct an arbitrary function (fourier decomposition, spherical harmonics, etc...), so we can test our ability to reconstruct various parameters. So, lets say we want to retrieve a property, P, that depends on our desired quantity dv (r)/dr P = Ψ(r) dv (r) dr dr
21 An Orthogonal Framework If we can write Ψ(r) in terms of our kernels, then it can be obtained directly from measurements P = Ψ(r) dv (r) dr dr dv (r) a i K i (r) dr = dr i i We can decompose P using our orthogonal functions P = dv (r) (Ψ φ i )φ i (r) dr dr i a i y i The accuracy of our retrieved quantity P, is then dictated by how well we can reconstruct Ψ(r) using our kernels
22 An Orthogonal Framework Applying this to the OSIRIS data we can test how well we can retrieve a couple aerosol properties, total number of particles, n aer, and total amount of light scattered by aerosols, k aer. Ψ naer = 3 4πr 3 Ψ kaer = 3 4πr 3 σ aer (r)
23 An Orthogonal Framework Number density is very poorly reconstructed, however extinction is actually quite good, implying extinction is a much more robust quantity for retrieval. However, errors may still be large if the majority of particles are small. Ψnaer(r) (1/µm 3 ) 4 2 Ψnaer(r) = 3 4πr 3 i (Ψ φ i)φi(r) Ψkaer (r) = 3 4πr3 σaer(λ, r) i (Ψ φ i)φi(r) 2 1 Ψkaer (r) (cm2 /µm 3 ).1 1 Radius (µm) 4
24 A Linear Inversion We can also use these orthogonal function to perform an inversion. First, write the retrieved quantity as a linear combination of our kernels (plus whatever we can t reconstruct) dv (r) dr = i a i φ i (r) + O N Then our measurements would be y i = K i (r) dv (r) dr = dr K i (a i φ i )dr which yields dv (r) = K T QΛ 1 Q T y dr Note the Λ 1 term, that may cause large error magnification for small eigenvalues i
25 A Linear Inversion As an example, I attempted to retrieve a few different particle size distributions using OSIRIS measurements dv (r)/dr (1 5 µm 1 ) rg=.8µm V =-1.8e-7 σg=1.3 rg=.15µm V =1.7e-6 σg=1.3 rg=.35µm V =5.4e-6 σg= rg=.8µm V =1.3e-6 σg=1.6 rg=.15µm V =3.9e-6 σg=1.6 rg=.35µm V =5.2e-6 σg= Radius (µm) rg=.8µm V =3.4e-6 σg=1.9 rg=.15µm V =4.7e-6 σg=1.9 rg=.35µm V =4.1e-6 σg=
26 A Non-Linear Inversion Linear analysis is nice for getting a grasp on the problem, but too inaccurate for actual retrievals - need a non-linear technique. System Output, y meas A priori System State, x System Model y mod = F(x) Retrieved System State Update State x (n) x (n+1)
27 A Non-Linear Inversion Provided you have a good forward model, the key here is updating the state vector. Lots of ways to do this Jacobian based methods: Newton, Gradient descent, Levenberg-Marquardt Quick to converge, slow to iterate Multiplicative methods: Chahine relaxation, Algebraic reconstruction Quick to iterate slow to converge Statistical methods: Optimal Estimation Requires detailed statistical knowledge of the system More complicated methods Genetic algorithms, simulated annealing Difficult implementation
28 A Non-Linear Inversion Most OSIRIS products use the Multiplicative Algebraic Reconstruction Technique (MART) since the measurements are not too non-linear and Jacobians are expensive to calculate. x (n+1) i = x (n) i j i y measij y modij (x (n) ) W kij However, aerosol measurements are more non-linear as a function of particle size and have multiple dimensions (number density and mode radius), so this method provides very slow convergence.
29 A Non-Linear Inversion Instead the aerosol retrieval uses the Levenberg-Marquardt algorithm (K T K + γdiag(k T K))δ = K T (y meas y mod ) x (n+1) = x (n) + δ This is basically a combination of Gauss-Newton s method and gradient descent, but generally provides faster convergence and improved robustness Downside is that it requires computation of the Jacobian, K, which takes time.
30 A Non-Linear Inversion The final retrieval algorithm 1 Assume aerosol consists of a lognormal distribution with variable mode radius and number density 2 Use two OSIRIS measurements at each altitude, one at 75 nm and one at 153 nm 3 Simulate measurements using forward model and update radius and density using LM 4 Compare simulation results with measurements and if they don t agree repeat step 2
31 Results
Inverse Theory. COST WaVaCS Winterschool Venice, February Stefan Buehler Luleå University of Technology Kiruna
Inverse Theory COST WaVaCS Winterschool Venice, February 2011 Stefan Buehler Luleå University of Technology Kiruna Overview Inversion 1 The Inverse Problem 2 Simple Minded Approach (Matrix Inversion) 3
More informationCS 542G: Robustifying Newton, Constraints, Nonlinear Least Squares
CS 542G: Robustifying Newton, Constraints, Nonlinear Least Squares Robert Bridson October 29, 2008 1 Hessian Problems in Newton Last time we fixed one of plain Newton s problems by introducing line search
More informationISTINA - : Investigation of Sensitivity Tendencies and Inverse Numerical Algorithm advances in aerosol remote sensing
STNA - : nvestigation of Sensitivity Tendencies and nverse Numerical Algorithm advances in aerosol remote sensing B. Torres, O. Dubovik, D. Fuertes, and P. Litvinov GRASP- SAS, LOA, Universite Lille-1,
More informationGetting Started with Communications Engineering
1 Linear algebra is the algebra of linear equations: the term linear being used in the same sense as in linear functions, such as: which is the equation of a straight line. y ax c (0.1) Of course, if we
More informationStratospheric Aerosol Particle Size Retrieval
Stratospheric Aerosol Particle Size Retrieval A Thesis Submitted to the College of Graduate Studies and Research in Partial Fulfillment of the Requirements for the degree of Master of Science in the Department
More informationIntroduction to gradient descent
6-1: Introduction to gradient descent Prof. J.C. Kao, UCLA Introduction to gradient descent Derivation and intuitions Hessian 6-2: Introduction to gradient descent Prof. J.C. Kao, UCLA Introduction Our
More informationRadiation in the atmosphere
Radiation in the atmosphere Flux and intensity Blackbody radiation in a nutshell Solar constant Interaction of radiation with matter Absorption of solar radiation Scattering Radiative transfer Irradiance
More informationSpaced-Based Measurements of Stratospheric Aerosols
Spaced-Based Measurements of Stratospheric Aerosols Larry W. Thomason NASA Langley Research Center Hampton, Virginia USA 6/17/2003 L. Thomason 1 Measurement by Extinction of Solar Radiation Stratospheric
More informationMATH 308 COURSE SUMMARY
MATH 308 COURSE SUMMARY Approximately a third of the exam cover the material from the first two midterms, that is, chapter 6 and the first six sections of chapter 7. The rest of the exam will cover the
More informationAbsorption and scattering
Absorption and scattering When a beam of radiation goes through the atmosphere, it encounters gas molecules, aerosols, cloud droplets, and ice crystals. These objects perturb the radiation field. Part
More informationPOLI 8501 Introduction to Maximum Likelihood Estimation
POLI 8501 Introduction to Maximum Likelihood Estimation Maximum Likelihood Intuition Consider a model that looks like this: Y i N(µ, σ 2 ) So: E(Y ) = µ V ar(y ) = σ 2 Suppose you have some data on Y,
More informationAlgorithm Document HEIDOSCILI
lgorithm Document for the retrieval of OClO, BrO and NO 2 vertical profiles from SCIMCHY limb measurements by HEIDOSCILI (Heidelberg DOS of SCIMCHY Limb measurements) uthors: Sven Kühl, Janis Pukite, Thomas
More informationInverse problems and uncertainty quantification in remote sensing
1 / 38 Inverse problems and uncertainty quantification in remote sensing Johanna Tamminen Finnish Meterological Institute johanna.tamminen@fmi.fi ESA Earth Observation Summer School on Earth System Monitoring
More information5 Linear Algebra and Inverse Problem
5 Linear Algebra and Inverse Problem 5.1 Introduction Direct problem ( Forward problem) is to find field quantities satisfying Governing equations, Boundary conditions, Initial conditions. The direct problem
More informationNumerical Methods for Inverse Kinematics
Numerical Methods for Inverse Kinematics Niels Joubert, UC Berkeley, CS184 2008-11-25 Inverse Kinematics is used to pose models by specifying endpoints of segments rather than individual joint angles.
More informationComparative Study Among Lease Square Method, Steepest Descent Method, and Conjugate Gradient Method for Atmopsheric Sounder Data Analysis
Comparative Study Among Lease Square Method, Steepest Descent Method, and Conjugate Gradient Method for Atmopsheric Sounder Data Analysis Kohei Arai 1 Graduate School of Science and Engineering Saga University
More informationFinal Review Sheet. B = (1, 1 + 3x, 1 + x 2 ) then 2 + 3x + 6x 2
Final Review Sheet The final will cover Sections Chapters 1,2,3 and 4, as well as sections 5.1-5.4, 6.1-6.2 and 7.1-7.3 from chapters 5,6 and 7. This is essentially all material covered this term. Watch
More informationWhat are Aerosols? Suspension of very small solid particles or liquid droplets Radii typically in the range of 10nm to
What are Aerosols? Suspension of very small solid particles or liquid droplets Radii typically in the range of 10nm to 10µm Concentrations decrease exponentially with height N(z) = N(0)exp(-z/H) Long-lived
More informationεx 2 + x 1 = 0. (2) Suppose we try a regular perturbation expansion on it. Setting ε = 0 gives x 1 = 0,
4 Rescaling In this section we ll look at one of the reasons that our ε = 0 system might not have enough solutions, and introduce a tool that is fundamental to all perturbation systems. We ll start with
More informationLaura Revell, Andrea Stenke, Beiping Luo, Stefanie Kremser, Eugene Rozanov, Timofei Sukhodolov and Thomas Peter
Impacts of Mt. Pinatubo volcanic aerosol on the tropical stratosphere in chemistry-climate model simulations using CCMI and CMIP6 stratospheric aerosol data Laura Revell, Andrea Stenke, Beiping Luo, Stefanie
More information7. Aerosols and Climate
7. Aerosols and Climate I. Scattering 1. When radiation impinges on a medium of small particles, scattering of some of the radiation occurs in all directions. The portion scattered backward is called the
More informationNumerical Methods - Numerical Linear Algebra
Numerical Methods - Numerical Linear Algebra Y. K. Goh Universiti Tunku Abdul Rahman 2013 Y. K. Goh (UTAR) Numerical Methods - Numerical Linear Algebra I 2013 1 / 62 Outline 1 Motivation 2 Solving Linear
More informationOutline. December 14, Applications Scattering. Chemical components. Forward model Radiometry Data retrieval. Applications in remote sensing
in in December 4, 27 Outline in 2 : RTE Consider plane parallel Propagation of a signal with intensity (radiance) I ν from the top of the to a receiver on Earth Take a layer of thickness dz Layer will
More informationThe Odin/OSIRIS time series from 2001 to now
The Odin/OSIRIS time series from 21 to now SPARC/IOC/WMO-IGACO workshop on Past Changes in the Vertical Distribution of Ozone Geneva, January 25-27 211 The Atmosphere as Seen from Odin Bright Dim.5 º The
More informationPrinciples of active remote sensing: Lidars and lidar sensing of aerosols, gases and clouds.
Lecture 14 Principles of active remote sensing: Lidars and lidar sensing of aerosols, gases and clouds. Objectives: 1. Optical interactions of relevance to lasers. 2. General principles of lidars. 3. Lidar
More informationLinear Algebra Review
Chapter 1 Linear Algebra Review It is assumed that you have had a course in linear algebra, and are familiar with matrix multiplication, eigenvectors, etc. I will review some of these terms here, but quite
More informationCSC321 Lecture 7: Optimization
CSC321 Lecture 7: Optimization Roger Grosse Roger Grosse CSC321 Lecture 7: Optimization 1 / 25 Overview We ve talked a lot about how to compute gradients. What do we actually do with them? Today s lecture:
More informationThe Spectral Radiative Effects of Inhomogeneous Clouds and Aerosols
The Spectral Radiative Effects of Inhomogeneous Clouds and Aerosols S. Schmidt, B. Kindel, & P. Pilewskie Laboratory for Atmospheric and Space Physics University of Colorado SORCE Science Meeting, 13-16
More informationMATH 310, REVIEW SHEET 2
MATH 310, REVIEW SHEET 2 These notes are a very short summary of the key topics in the book (and follow the book pretty closely). You should be familiar with everything on here, but it s not comprehensive,
More informationLinear smoother. ŷ = S y. where s ij = s ij (x) e.g. s ij = diag(l i (x))
Linear smoother ŷ = S y where s ij = s ij (x) e.g. s ij = diag(l i (x)) 2 Online Learning: LMS and Perceptrons Partially adapted from slides by Ryan Gabbard and Mitch Marcus (and lots original slides by
More informationCSC321 Lecture 8: Optimization
CSC321 Lecture 8: Optimization Roger Grosse Roger Grosse CSC321 Lecture 8: Optimization 1 / 26 Overview We ve talked a lot about how to compute gradients. What do we actually do with them? Today s lecture:
More informationComputational Methods. Least Squares Approximation/Optimization
Computational Methods Least Squares Approximation/Optimization Manfred Huber 2011 1 Least Squares Least squares methods are aimed at finding approximate solutions when no precise solution exists Find the
More informationApplied Mathematics 205. Unit I: Data Fitting. Lecturer: Dr. David Knezevic
Applied Mathematics 205 Unit I: Data Fitting Lecturer: Dr. David Knezevic Unit I: Data Fitting Chapter I.4: Nonlinear Least Squares 2 / 25 Nonlinear Least Squares So far we have looked at finding a best
More informationStratospheric aerosol profile retrieval from SCIAMACHY limb observations
Stratospheric aerosol profile retrieval from SCIAMACHY limb observations Yang Jingmei Zong Xuemei Key Laboratory of Middle Atmosphere and Global Environment Observation (LAGEO), Institute of Atmospheric
More informationAircraft Validation of Infrared Emissivity derived from Advanced InfraRed Sounder Satellite Observations
Aircraft Validation of Infrared Emissivity derived from Advanced InfraRed Sounder Satellite Observations Robert Knuteson, Fred Best, Steve Dutcher, Ray Garcia, Chris Moeller, Szu Chia Moeller, Henry Revercomb,
More informationCHAPTER 8. AEROSOLS 8.1 SOURCES AND SINKS OF AEROSOLS
1 CHAPTER 8 AEROSOLS Aerosols in the atmosphere have several important environmental effects They are a respiratory health hazard at the high concentrations found in urban environments They scatter and
More informationContents. 1 Introduction. 1.1 History of Optimization ALG-ML SEMINAR LISSA: LINEAR TIME SECOND-ORDER STOCHASTIC ALGORITHM FEBRUARY 23, 2016
ALG-ML SEMINAR LISSA: LINEAR TIME SECOND-ORDER STOCHASTIC ALGORITHM FEBRUARY 23, 2016 LECTURERS: NAMAN AGARWAL AND BRIAN BULLINS SCRIBE: KIRAN VODRAHALLI Contents 1 Introduction 1 1.1 History of Optimization.....................................
More informationMonday, Oct. 2: Clear-sky radiation; solar attenuation, Thermal. nomenclature
Monday, Oct. 2: Clear-sky radiation; solar attenuation, Thermal nomenclature Sun Earth Y-axis: Spectral radiance, aka monochromatic intensity units: watts/(m^2*ster*wavelength) Blackbody curves provide
More informationIn this method, one defines
1 Feautrier s Method for Radiative Transfer This project is to use Feautrier s method to solve the Gray Atmosphere problem for a plane-parallel atmosphere Although we have developed alternate techniques
More informationAuthors response to the reviewers comments
Manuscript No.: amtd-3-c1225-2010 Authors response to the reviewers comments Title: Satellite remote sensing of Asian aerosols: A case study of clean, polluted, and Asian dust storm days General comments:
More informationIterative Methods. Splitting Methods
Iterative Methods Splitting Methods 1 Direct Methods Solving Ax = b using direct methods. Gaussian elimination (using LU decomposition) Variants of LU, including Crout and Doolittle Other decomposition
More informationAlgorithm document for SCIAMACHY Stratozone limb ozone profile retrievals
Algorithm document for SCIAMACHY Stratozone limb ozone profile retrievals Date of origin: February 6, 2006 Author: Christian von Savigny Institute of Environmental Physics University of Bremen Otto-Hahn-Allee
More informationSingular Value Decomposition. 1 Singular Value Decomposition and the Four Fundamental Subspaces
Singular Value Decomposition This handout is a review of some basic concepts in linear algebra For a detailed introduction, consult a linear algebra text Linear lgebra and its pplications by Gilbert Strang
More informationStatistical Machine Learning
Statistical Machine Learning Christoph Lampert Spring Semester 2015/2016 // Lecture 12 1 / 36 Unsupervised Learning Dimensionality Reduction 2 / 36 Dimensionality Reduction Given: data X = {x 1,..., x
More information1 Characteristics of the Atmosphere
CHAPTER 1 1 Characteristics of the Atmosphere SECTION The Atmosphere BEFORE YOU READ After you read this section, you should be able to answer these questions: What is Earth s atmosphere made of? How do
More informationmin f(x). (2.1) Objectives consisting of a smooth convex term plus a nonconvex regularization term;
Chapter 2 Gradient Methods The gradient method forms the foundation of all of the schemes studied in this book. We will provide several complementary perspectives on this algorithm that highlight the many
More informationGroup Representations
Group Representations Alex Alemi November 5, 2012 Group Theory You ve been using it this whole time. Things I hope to cover And Introduction to Groups Representation theory Crystallagraphic Groups Continuous
More informationStratospheric aerosol retrieval with optical spectrograph and infrared imaging system limb scatter measurements
JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 112,, doi:10.1029/2006jd008079, 2007 Stratospheric aerosol retrieval with optical spectrograph and infrared imaging system limb scatter measurements A. E. Bourassa,
More informationCOMBINED OZONE RETRIEVAL USING THE MICHELSON INTERFEROMETER FOR PASSIVE ATMOSPHERIC SOUNDING (MIPAS) AND THE TROPOSPHERIC EMISSION SPECTROMETER (TES)
COMBINED OZONE RETRIEVAL USING THE MICHELSON INTERFEROMETER FOR PASSIVE ATMOSPHERIC SOUNDING (MIPAS) AND THE TROPOSPHERIC EMISSION SPECTROMETER (TES) Claire Waymark, Anu Dudhia, John Barnett, Fred Taylor,
More informationWhy integrate? 1 Introduction. 2 Integrating constant functions
Why integrate? 1 Introduction Teaching physics my colleagues and I often find that some of you find it difficult to see when you need to use integration. This session will explain the need for integration
More informationCourse outline, objectives, workload, projects, expectations
Course outline, objectives, workload, projects, expectations Introductions Remote Sensing Overview Elements of a remote sensing observing system 1. platform (satellite, surface, etc) 2. experimental design
More informationThe Radiative Transfer Equation
The Radiative Transfer Equation R. Wordsworth April 11, 215 1 Objectives Derive the general RTE equation Derive the atmospheric 1D horizontally homogenous RTE equation Look at heating/cooling rates in
More informationLecture 14. Principles of active remote sensing: Lidars. Lidar sensing of gases, aerosols, and clouds.
Lecture 14. Principles of active remote sensing: Lidars. Lidar sensing of gases, aerosols, and clouds. 1. Optical interactions of relevance to lasers. 2. General principles of lidars. 3. Lidar equation.
More informationPrincipal Component Analysis (PCA) for Sparse High-Dimensional Data
AB Principal Component Analysis (PCA) for Sparse High-Dimensional Data Tapani Raiko, Alexander Ilin, and Juha Karhunen Helsinki University of Technology, Finland Adaptive Informatics Research Center Principal
More informationPhysics 202 Laboratory 5. Linear Algebra 1. Laboratory 5. Physics 202 Laboratory
Physics 202 Laboratory 5 Linear Algebra Laboratory 5 Physics 202 Laboratory We close our whirlwind tour of numerical methods by advertising some elements of (numerical) linear algebra. There are three
More informationPose Tracking II! Gordon Wetzstein! Stanford University! EE 267 Virtual Reality! Lecture 12! stanford.edu/class/ee267/!
Pose Tracking II! Gordon Wetzstein! Stanford University! EE 267 Virtual Reality! Lecture 12! stanford.edu/class/ee267/!! WARNING! this class will be dense! will learn how to use nonlinear optimization
More informationECEN 615 Methods of Electric Power Systems Analysis Lecture 18: Least Squares, State Estimation
ECEN 615 Methods of Electric Power Systems Analysis Lecture 18: Least Squares, State Estimation Prof. om Overbye Dept. of Electrical and Computer Engineering exas A&M University overbye@tamu.edu Announcements
More informationDot Products, Transposes, and Orthogonal Projections
Dot Products, Transposes, and Orthogonal Projections David Jekel November 13, 2015 Properties of Dot Products Recall that the dot product or standard inner product on R n is given by x y = x 1 y 1 + +
More informationUKCA_RADAER Aerosol-radiation interactions
UKCA_RADAER Aerosol-radiation interactions Nicolas Bellouin UKCA Training Workshop, Cambridge, 8 January 2015 University of Reading 2014 n.bellouin@reading.ac.uk Lecture summary Why care about aerosol-radiation
More informationNon-polynomial Least-squares fitting
Applied Math 205 Last time: piecewise polynomial interpolation, least-squares fitting Today: underdetermined least squares, nonlinear least squares Homework 1 (and subsequent homeworks) have several parts
More informationRetrieval problems in Remote Sensing
Maarten Krol August 2006 m.c.krol@phys.uu.nl 1. Goal To show that remote sensing is an indirect measurement, i.e. that we do not measure directly, but by means of (for instance) electromagnetic radiation
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 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 informationATOC 3500/CHEM 3152 Week 9, March 8, 2016
ATOC 3500/CHEM 3152 Week 9, March 8, 2016 Hand back Midterm Exams (average = 84) Interaction of atmospheric constituents with light Haze and Visibility Aerosol formation processes (more detail) Haze and
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 informationUnconstrained optimization
Chapter 4 Unconstrained optimization An unconstrained optimization problem takes the form min x Rnf(x) (4.1) for a target functional (also called objective function) f : R n R. In this chapter and throughout
More informationChapter 7 Iterative Techniques in Matrix Algebra
Chapter 7 Iterative Techniques in Matrix Algebra Per-Olof Persson persson@berkeley.edu Department of Mathematics University of California, Berkeley Math 128B Numerical Analysis Vector Norms Definition
More informationLINEAR ALGEBRA: NUMERICAL METHODS. Version: August 12,
LINEAR ALGEBRA: NUMERICAL METHODS. Version: August 12, 2000 74 6 Summary Here we summarize the most important information about theoretical and numerical linear algebra. MORALS OF THE STORY: I. Theoretically
More informationLecture 6: Lies, Inner Product Spaces, and Symmetric Matrices
Math 108B Professor: Padraic Bartlett Lecture 6: Lies, Inner Product Spaces, and Symmetric Matrices Week 6 UCSB 2014 1 Lies Fun fact: I have deceived 1 you somewhat with these last few lectures! Let me
More informationInversion of Phase Data for a Phase Velocity Map 101. Summary for CIDER12 Non-Seismologists
Inversion of Phase Data for a Phase Velocity Map 101 Summary for CIDER12 Non-Seismologists 1. Setting up a Linear System of Equations This is a quick-and-dirty, not-peer reviewed summary of how a dataset
More informationSatellite Retrieval Assimilation (a modified observation operator)
Satellite Retrieval Assimilation (a modified observation operator) David D. Kuhl With Istvan Szunyogh and Brad Pierce Sept. 21 2009 Weather Chaos Group Meeting Overview Introduction Derivation of retrieval
More informationSummary of various integrals
ummary of various integrals Here s an arbitrary compilation of information about integrals Moisés made on a cold ecember night. 1 General things o not mix scalars and vectors! In particular ome integrals
More informationMEASURING TRACE GAS PROFILES FROM SPACE
MEASURING TRACE GAS PROFILES FROM SPACE Caroline Nowlan Atomic and Molecular Physics Division Harvard-Smithsonian Center for Astrophysics Collaborators: Kelly Chance, Xiong Liu, Gonzalo Gonzalez Abad,
More informationLecture 10: Powers of Matrices, Difference Equations
Lecture 10: Powers of Matrices, Difference Equations Difference Equations A difference equation, also sometimes called a recurrence equation is an equation that defines a sequence recursively, i.e. each
More informationLeast squares problems
Least squares problems How to state and solve them, then evaluate their solutions Stéphane Mottelet Université de Technologie de Compiègne 30 septembre 2016 Stéphane Mottelet (UTC) Least squares 1 / 55
More informationAPPLICATIONS WITH METEOROLOGICAL SATELLITES. W. Paul Menzel. Office of Research and Applications NOAA/NESDIS University of Wisconsin Madison, WI
APPLICATIONS WITH METEOROLOGICAL SATELLITES by W. Paul Menzel Office of Research and Applications NOAA/NESDIS University of Wisconsin Madison, WI July 2004 Unpublished Work Copyright Pending TABLE OF CONTENTS
More information7.2 Steepest Descent and Preconditioning
7.2 Steepest Descent and Preconditioning Descent methods are a broad class of iterative methods for finding solutions of the linear system Ax = b for symmetric positive definite matrix A R n n. Consider
More informationBindel, Fall 2016 Matrix Computations (CS 6210) Notes for
1 Iteration basics Notes for 2016-11-07 An iterative solver for Ax = b is produces a sequence of approximations x (k) x. We always stop after finitely many steps, based on some convergence criterion, e.g.
More informationRadiation in climate models.
Lecture. Radiation in climate models. Objectives:. A hierarchy of the climate models.. Radiative and radiative-convective equilibrium.. Examples of simple energy balance models.. Radiation in the atmospheric
More informationFrequency, Vibration, and Fourier
Lecture 22: Frequency, Vibration, and Fourier Computer Graphics CMU 15-462/15-662, Fall 2015 Last time: Numerical Linear Algebra Graphics via linear systems of equations Why linear? Have to solve BIG problems
More informationj=1 r 1 x 1 x n. r m r j (x) r j r j (x) r j (x). r j x k
Maria Cameron Nonlinear Least Squares Problem The nonlinear least squares problem arises when one needs to find optimal set of parameters for a nonlinear model given a large set of data The variables x,,
More informationNumerical Linear Algebra Primer. Ryan Tibshirani Convex Optimization /36-725
Numerical Linear Algebra Primer Ryan Tibshirani Convex Optimization 10-725/36-725 Last time: proximal gradient descent Consider the problem min g(x) + h(x) with g, h convex, g differentiable, and h simple
More informationCourse Description for Real Analysis, Math 156
Course Description for Real Analysis, Math 156 In this class, we will study elliptic PDE, Fourier analysis, and dispersive PDE. Here is a quick summary of the topics will study study. They re described
More informationInstructor (Brad Osgood)
TheFourierTransformAndItsApplications-Lecture26 Instructor (Brad Osgood): Relax, but no, no, no, the TV is on. It's time to hit the road. Time to rock and roll. We're going to now turn to our last topic
More information[Disclaimer: This is not a complete list of everything you need to know, just some of the topics that gave people difficulty.]
Math 43 Review Notes [Disclaimer: This is not a complete list of everything you need to know, just some of the topics that gave people difficulty Dot Product If v (v, v, v 3 and w (w, w, w 3, then the
More informationPartial Fractions. June 27, In this section, we will learn to integrate another class of functions: the rational functions.
Partial Fractions June 7, 04 In this section, we will learn to integrate another class of functions: the rational functions. Definition. A rational function is a fraction of two polynomials. For example,
More informationAnalysis Methods in Atmospheric and Oceanic Science
Analysis Methods in Atmospheric and Oceanic Science Copyright 2016 University of Maryland. 1 AOSC 652 Root Finding & Function Minimization Week 7, Day 2 12 Oct 2016 2 Method #3: Newton-Raphson Root is
More informationJOINT RETRIEVAL OF CO AND VIBRATIONAL TEMPERATURE FROM MIPAS-ENVISAT
JOINT RETRIEVAL OF CO AND VIBRATIONAL TEMPERATURE FROM MIPAS-ENVISAT Joanne Walker and Anu Dudhia Atmospheric, Oceanic and Planetary Physics, Oxford Universtity, UK ABSTRACT MIPAS is a limb viewing fourier
More informationDRAGON ADVANCED TRAINING COURSE IN ATMOSPHERE REMOTE SENSING. Inversion basics. Erkki Kyrölä Finnish Meteorological Institute
Inversion basics y = Kx + ε x ˆ = (K T K) 1 K T y Erkki Kyrölä Finnish Meteorological Institute Day 3 Lecture 1 Retrieval techniques - Erkki Kyrölä 1 Contents 1. Introduction: Measurements, models, inversion
More informationComparison of AERONET inverted size distributions to measured distributions from the Aerodyne Aerosol Mass Spectrometer
Comparison of inverted size distributions to measured distributions from the Aerodyne Aerosol Mass Spectrometer Peter DeCarlo Remote Sensing Project April 28, 23 Introduction The comparison of direct in-situ
More informationIsomorphisms and Well-definedness
Isomorphisms and Well-definedness Jonathan Love October 30, 2016 Suppose you want to show that two groups G and H are isomorphic. There are a couple of ways to go about doing this depending on the situation,
More informationRetrieval of tropospheric methane from MOPITT measurements: algorithm description and simulations
Retrieval of tropospheric methane from MOPITT measurements: algorithm description and simulations Merritt N. Deeter*, Jinxue Wang, John C. Gille, and Paul L. Bailey National Center for Atmospheric Research,
More informationLecture 22. r i+1 = b Ax i+1 = b A(x i + α i r i ) =(b Ax i ) α i Ar i = r i α i Ar i
8.409 An Algorithmist s oolkit December, 009 Lecturer: Jonathan Kelner Lecture Last time Last time, we reduced solving sparse systems of linear equations Ax = b where A is symmetric and positive definite
More informationJoint Temperature, Humidity, and Sea Surface Temperature Retrieval from IASI Sensor Data
Joint Temperature, Humidity, and Sea Surface Temperature Retrieval from IASI Sensor Data Marc Schwaerz and Gottfried Kirchengast Institute for Geophysics, Astrophysics, and Meteorology (IGAM), University
More informationNon-Convex Optimization. CS6787 Lecture 7 Fall 2017
Non-Convex Optimization CS6787 Lecture 7 Fall 2017 First some words about grading I sent out a bunch of grades on the course management system Everyone should have all their grades in Not including paper
More informationLab 5 The Green House Effect Pre-lab
Lab 5 The Green House Effect Pre-lab NAME Partner Given information: Radius of the Sun: 6.957 X 10 8 m Radius of the Earth: 6.371 X 10 6 m Emissivity of Sun: ~1.0 Surface temperature of the Sun: 5,778
More informationInt= surface density change (cm/yr water equivalent)
Earth s changing mass distribution: Making the most of GRACE Frederik J Simons F. A. Dahlen, Jessica C. Hawthorne, Dong V. Wang Princeton University Mark A. Wieczorek IPG Paris Lei Wang Ohio State University
More informationMarkov chain Monte Carlo methods in atmospheric remote sensing
1 / 45 Markov chain Monte Carlo methods in atmospheric remote sensing Johanna Tamminen johanna.tamminen@fmi.fi ESA Summer School on Earth System Monitoring and Modeling July 3 Aug 11, 212, Frascati July,
More informationLAB 8: INTEGRATION. Figure 1. Approximating volume: the left by cubes, the right by cylinders
LAB 8: INTGRATION The purpose of this lab is to give intuition about integration. It will hopefully complement the, rather-dry, section of the lab manual and the, rather-too-rigorous-and-unreadable, section
More informationOdin-OSIRIS: A Summary of the Results from the Past Eleven Years
Odin-OSIRIS: A Summary of the Results from the Past Eleven Years ESA ATMOS 12 June 18, 12 Bruges, Brussels : Year Eleven of a Two Year Mission ESA ATMOS 12 June 18, 12 Bruges, Brussels OSIRIS Designed
More information