24.07.2014 Views

Rotational Raman scattering in the Earth's atmosphere ... - SRON

Rotational Raman scattering in the Earth's atmosphere ... - SRON

Rotational Raman scattering in the Earth's atmosphere ... - SRON

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

88 Chapter 4<br />

4.4 Retrieval of a solar spectrum on a f<strong>in</strong>e spectral grid from <strong>the</strong><br />

solar measurement<br />

In this section we describe <strong>the</strong> GOME measurements of <strong>the</strong> Earth radiance and <strong>the</strong> solar irradiance<br />

spectrum by <strong>the</strong> vectors y ear and y sun , respectively. Both vectors have dimension M and cover <strong>the</strong><br />

same spectral range 390–400 nm. Each element of <strong>the</strong>se vectors describes <strong>the</strong> measurement of one<br />

particular detector pixel and can be simulated us<strong>in</strong>g Eqs. (4.1) and (4.10). Fur<strong>the</strong>rmore, <strong>the</strong> cont<strong>in</strong>uous<br />

solar spectrum F 0 is described by a vector x of dimension N, which belongs to a discretization of <strong>the</strong><br />

solar spectrum on a 1 cm −1 wavenumber grid. The forward simulation of GOME earthsh<strong>in</strong>e and solar<br />

measurements can be described by <strong>the</strong> l<strong>in</strong>ear equations<br />

y ear = K ear x + e ear , (4.11)<br />

y sun = K sun x + e sun , (4.12)<br />

where <strong>the</strong> M ×N-dimensional kernel matrices K ear and K sun are matrix representations of <strong>the</strong> correspond<strong>in</strong>g<br />

<strong>in</strong>tegral kernels k ear and k sun <strong>in</strong> <strong>the</strong> Eqs. (4.1) and (4.10). The M-dimensional vectors e ear<br />

and e sun conta<strong>in</strong> <strong>the</strong> measurement error of <strong>the</strong> Earth radiance and solar irradiance spectrum, respectively.<br />

In this section we discuss a retrieval scheme to determ<strong>in</strong>e <strong>the</strong> solar spectrum x from <strong>the</strong> GOME<br />

solar measurements. The retrieved solar spectrum should not be seen as a GOME data product, but<br />

as an auxiliary product to enhance <strong>the</strong> accuracy of GOME Earth radiance simulations, which will be<br />

demonstrated at <strong>the</strong> end of this section.<br />

4.4.1 Inversion scheme<br />

The retrieval of <strong>the</strong> solar spectrum x from <strong>the</strong> GOME solar irradiance spectrum y sun <strong>in</strong>volves <strong>the</strong><br />

<strong>in</strong>version of Eq. (4.12). This presents an underdeterm<strong>in</strong>ed problem because <strong>the</strong> number of measurements<br />

M is smaller than <strong>the</strong> number of parameters N to be retrieved. To ga<strong>in</strong> <strong>in</strong>sight <strong>in</strong>to <strong>the</strong> <strong>in</strong>version<br />

of Eq. (4.12) we consider <strong>the</strong> s<strong>in</strong>gular value decomposition (SVD) of <strong>the</strong> kernel matrix K, viz.<br />

K = UΣV T . (4.13)<br />

For simplicity, <strong>the</strong> subscript ‘sun’ is omitted <strong>in</strong> this subsection. In Eq. (4.13), U is a M × N-matrix<br />

with orthonormal columns (u 1 ,...,u M ) spann<strong>in</strong>g <strong>the</strong> measurement space, V is a N ×N-matrix with<br />

orthonormal columns (v 1 ,...,v N ) spann<strong>in</strong>g <strong>the</strong> state space and Σ is a diagonal matrix conta<strong>in</strong><strong>in</strong>g <strong>the</strong><br />

s<strong>in</strong>gular values (σ 1 ,...,σ N ). The underdeterm<strong>in</strong>ation of <strong>the</strong> <strong>in</strong>version problem is expressed by <strong>the</strong><br />

fact that <strong>the</strong> s<strong>in</strong>gular values (σ M+1 ,...,σ N ) are zero. Therefore <strong>the</strong> estimate of <strong>the</strong> <strong>in</strong>version needs<br />

to be written as a truncated sum<br />

x est =<br />

M∑<br />

n=1<br />

(u T ny)<br />

σ n<br />

v n . (4.14)<br />

The part of <strong>the</strong> solar spectrum that is described by <strong>the</strong> vectors v M+1 ,...,v N cannot be retrieved. The<br />

vector space spanned by <strong>the</strong>se vectors is referred to as <strong>the</strong> null-space of <strong>the</strong> problem, whereas <strong>the</strong>

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!