DIFA, University of Basilicata, Italy

Similar documents
Data Assimilation. Alan O Neill Data Assimilation Research Centre University of Reading

Review of Gaussian Quadrature method

Satellite Retrieval Data Assimilation

Duality # Second iteration for HW problem. Recall our LP example problem we have been working on, in equality form, is given below.

The heat budget of the atmosphere and the greenhouse effect

ESCI 343 Atmospheric Dynamics II Lesson 14 Inertial/slantwise Instability

Lesson 8. Thermomechanical Measurements for Energy Systems (MENR) Measurements for Mechanical Systems and Production (MMER)

Classical Mechanics. From Molecular to Con/nuum Physics I WS 11/12 Emiliano Ippoli/ October, 2011

The International Association for the Properties of Water and Steam. Release on the Ionization Constant of H 2 O

Numerical Methods I Orthogonal Polynomials

Solutions of Klein - Gordan equations, using Finite Fourier Sine Transform

Remarks to the H-mode workshop paper

When e = 0 we obtain the case of a circle.

Conservation Law. Chapter Goal. 5.2 Theory

38.2. The Uniform Distribution. Introduction. Prerequisites. Learning Outcomes

SOLUTION OF QUADRATIC NONLINEAR PROBLEMS WITH MULTIPLE SCALES LINDSTEDT-POINCARE METHOD. Mehmet Pakdemirli and Gözde Sarı

PHYSICS 211 MIDTERM I 21 April 2004

Lecture Notes PH 411/511 ECE 598 A. La Rosa Portland State University INTRODUCTION TO QUANTUM MECHANICS

Carmine Serio, Guido Masiello. DIFA, University of Basilicata, Italy

Why symmetry? Symmetry is often argued from the requirement that the strain energy must be positive. (e.g. Generalized 3-D Hooke s law)

We will see what is meant by standard form very shortly

Industrial Electrical Engineering and Automation

Discrete Mathematics and Probability Theory Summer 2014 James Cook Note 17

Incorporating Ensemble Covariance in the Gridpoint Statistical Interpolation Variational Minimization: A Mathematical Framework

Partial Derivatives. Limits. For a single variable function f (x), the limit lim


THERMAL EXPANSION COEFFICIENT OF WATER FOR VOLUMETRIC CALIBRATION

Orthogonal Polynomials

Discrete Mathematics and Probability Theory Spring 2013 Anant Sahai Lecture 17

G. MATEESCU 1 A. MATEESCU 1 C. SAMOILĂ 2 D. URSUŢIU 2

Electromagnetics P5-1. 1) Physical quantities in EM could be scalar (charge, current, energy) or vector (EM fields).

Consequently, the temperature must be the same at each point in the cross section at x. Let:

Spanning tree congestion of some product graphs

M344 - ADVANCED ENGINEERING MATHEMATICS

Physics based RS algorithms for water quality parameters: Status and challenges. ZhongPing Lee

Driving Cycle Construction of City Road for Hybrid Bus Based on Markov Process Deng Pan1, a, Fengchun Sun1,b*, Hongwen He1, c, Jiankun Peng1, d

Week 10: Line Integrals

Can the Phase I problem be unfeasible or unbounded? -No

The Islamic University of Gaza Faculty of Engineering Civil Engineering Department. Numerical Analysis ECIV Chapter 11

13: Diffusion in 2 Energy Groups

The First Fundamental Theorem of Calculus. If f(x) is continuous on [a, b] and F (x) is any antiderivative. f(x) dx = F (b) F (a).

Predict Global Earth Temperature using Linier Regression

Heavy tail and stable distributions

TANDEM QUEUE WITH THREE MULTISERVER UNITS AND BULK SERVICE WITH ACCESSIBLE AND NON ACCESSBLE BATCH IN UNIT III WITH VACATION

Integration of tensor fields

Multi-objective optimization of dielectric layer photonic crystal filter

Week 10: Riemann integral and its properties

3.4 Conic sections. In polar coordinates (r, θ) conics are parameterized as. Next we consider the objects resulting from

CBE 291b - Computation And Optimization For Engineers

Fundamentals of Analytical Chemistry

SUPPLEMENTARY INFORMATION

Reinforcement learning II

The Algebra (al-jabr) of Matrices

Atmospheric Radiation Fall 2008

Application of Exp-Function Method to. a Huxley Equation with Variable Coefficient *

Math 33A Discussion Example Austin Christian October 23, Example 1. Consider tiling the plane by equilateral triangles, as below.

Mathematics Number: Logarithms

Direct Design of Orthogonal Filter Banks and Wavelets by Sequential Convex Quadratic Programming

14.4. Lengths of curves and surfaces of revolution. Introduction. Prerequisites. Learning Outcomes

Chapter 0. What is the Lebesgue integral about?

New Expansion and Infinite Series

ON PREY-PREDATOR MODEL WITH HOLLING-TYPE II AND LESLIE-GOWER SCHEMES AHMED BUSERI ASHINE

MIXED MODELS (Sections ) I) In the unrestricted model, interactions are treated as in the random effects model:

Department of Physical Pharmacy and Pharmacokinetics Poznań University of Medical Sciences Pharmacokinetics laboratory

Before we can begin Ch. 3 on Radicals, we need to be familiar with perfect squares, cubes, etc. Try and do as many as you can without a calculator!!!

Student Activity 3: Single Factor ANOVA

Jim Lambers MAT 169 Fall Semester Lecture 4 Notes

Calculus and linear algebra for biomedical engineering Week 11: The Riemann integral and its properties

Chapter Direct Method of Interpolation More Examples Electrical Engineering

dt. However, we might also be curious about dy

MATHS NOTES. SUBJECT: Maths LEVEL: Higher TEACHER: Aidan Roantree. The Institute of Education Topics Covered: Powers and Logs

Designing Information Devices and Systems I Spring 2018 Homework 7

Lecture 8 Wrap-up Part1, Matlab

Wavelets. Toh Kim Chuan National University of Singapore

Vadose Zone Hydrology

INTRODUCTION TO LINEAR ALGEBRA

Flow in porous media

Seismic Attributes used for Reservoir Simulation: Application to a Heavy Oil Reservoir in Canada

Terminal Velocity and Raindrop Growth

and that at t = 0 the object is at position 5. Find the position of the object at t = 2.

UNIFORM CONVERGENCE. Contents 1. Uniform Convergence 1 2. Properties of uniform convergence 3

Contents. Outline. Structured Rank Matrices Lecture 2: The theorem Proofs Examples related to structured ranks References. Structure Transport

Operations with Polynomials

MTH 122 Fall 2008 Essex County College Division of Mathematics Handout Version 10 1 October 14, 2008

Jack Simons, Henry Eyring Scientist and Professor Chemistry Department University of Utah

The Predom module. Predom calculates and plots isothermal 1-, 2- and 3-metal predominance area diagrams. Predom accesses only compound databases.

7.6 The Use of Definite Integrals in Physics and Engineering

UNIVERSITY OF KWAZULU-NATAL EXAMINATIONS: JUNE 2007

CMDA 4604: Intermediate Topics in Mathematical Modeling Lecture 19: Interpolation and Quadrature

The Wave Equation I. MA 436 Kurt Bryan

Best Approximation. Chapter The General Case

AQA Further Pure 1. Complex Numbers. Section 1: Introduction to Complex Numbers. The number system

A Matrix Algebra Primer

Lorenz Curve and Gini Coefficient in Right Truncated Pareto s Income Distribution

Application of Exact Discretization for Logistic Differential Equations to the Design of a Discrete-Time State-Observer

MATH34032: Green s Functions, Integral Equations and the Calculus of Variations 1

2.57/2.570 Midterm Exam No. 1 March 31, :00 am -12:30 pm

Lecture 20: Numerical Integration III

Tests for the Ratio of Two Poisson Rates

Worksheet #2 Math 285 Name: 1. Solve the following systems of linear equations. The prove that the solutions forms a subspace of

Transcription:

Guido Msiello nd Crmine Serio DIFA, Uniersity of Bsilict, Itly 1

Bsic methodologicl steps to retriee surfce emissiity Step 1: Represent emissiity with Fourier cosine truncted series to lower it dimensionlity below tht of the IASI spectrum. Step 2: Constrin the retriel with Lbortory mesurements Step 3: Blnce between Atmospheric nd Emissiity Constrints with 2-Dimensionl L-cure criterion 2

Bsic methodologicl steps to retriee surfce emissiity Step 2: Constrin the retriel with Lbortory mesurements Guido Msiello, Crmine Serio, nd Vincenzo Cuomo, "Eploiting Qurtz Spectrl Signture for the Detection of Cloud-Affected Stellite Infrred Obsertions oer Africn Desert Ares," Appl. Opt. 43, 2305-2315 2004 3

4 Bckground Retriel Methodology: the retriel problem is formulted within the contet of optiml estimtion 1 1 T T S K y S K y + ε 1 1 T T S F R S F R + ε o o o o F K F R y = = = = = How we modify the Physics to include surfce emissiity Linerize

5 How we modify the Physics to include surfce emissiity Use the logit trnsform to mp emissiity from the interl [0,1] to the rnge [-, + ] ε i y i = logit ε i = log i = 1,, N ch ; Nch = 1 ε i IASI chnnels which hs the inerse ep y i ε i = i = 1,, Nch i 1+ ep y i

6 To retriee for emissiity, the gien rdince Ri, with i the chnnel, is first linerized lso with respect to the function yi, tht is we consider in the inerse problem lso liner term of the type

Continue. Second, we deelop the function in truncted cosine series 7

Continue. nd Inserting the truncted cosine trnsform within the liner term 3, we get 8

9 The Fourier trnsform llows us to work in terms of spectrl resolution, ectly the wy we del with this concept. For se emissity, 60 Fourier Coefficients re enough

but, if we wnt to resole the Qurtz Resthrlen bnds in desert soil we need more thn 200 Fourier Coefficients 10

11 Bckground Retriel Methodology: the retriel problem is formulted within the contet of optiml estimtion 1 1 T T S K y S K y + ε 1 1 T T S F R S F R + ε o o o o F K F R y = = = = = From the size of bckground constrint we introduce informtion from lbortory mesurements Linerize

12 Lbortory emissiity is used for the bckground, men nd corince mtri ssumed to be digonl ASTER-Slisbury dt bse

The whole corince mtri for tmospheric prmeters nd emissiity is built up in block-digonl mtri 13 Atmosphere S γ 1S 0 Atm 0 = Emis γ 2S Emissiity γ 1 nd γ 2 cn be optimized to blnce between the two terms. Blncing is obtined with n originl nd fully nlyticl implementtion of 2-D L-cure criterion

14 Fully 2-D L-cure method, outline of the mthemtics inoled

Retriel eercise oer desert res 16 Shr Desert Arbin Desert Nmibi Desert Khlri Snn July 22, 2007, 6:45 Arbi 8:35 Other

Results: Nmibi desert 17 ε in LW lower thn in the SW. It mens fine grin of snds

Simultneously retrieed 18 with εσ: T s, T, H 2 O, O 3

Arbin desert 19 Corse piels Western

Shr desert 20

Klhri Snn 21

Conclusions 22 We he deeloped physicl inerse methodology to retriee the emissiity spectrum simultneously with Surfce Temperture nd Atmospheric prmeters: Temperture, wter pour nd ozone. The methodology relies minly on three bsic ides Deelop the emissiity spectrum in truncted Fourier cosine series Constrin the solution with Lbortory mesurements Blnce the optiml estimtion finl product with 2-Dimensionl L-cure criterion A test retriel eercise with IASI obsertions oer desert re shows tht the retrieed emissiity spectrum is cpble to cpture the fine detils of the surfce emission, een with non committl bckground corince mtri for emissiity The methodology will be soon pplied to derie mps t globl scle of the emissiity spectrum.

Set of prmeters retrieed 23 with φ-iasi Simultneously Emissiity spectrum Skin Temperture Temperture profile Wter pour profile Ozone profile Sequentilly, column mount CO CO 2 CH 4 N 2 O