progress (…) on forward modelling l1 and l2 …...rom saf workshop, ecmwf, 16–18 june 2014 1 ian...

28
ROM SAF Workshop, ECMWF, 16–18 June 2014 1 Ian Culverwell, Met Office Sean Healy, ECMWF Progress (…) on forward modelling L1 and L2 bending angles

Upload: others

Post on 07-Jul-2020

4 views

Category:

Documents


0 download

TRANSCRIPT

Page 1: Progress (…) on forward modelling L1 and L2 …...ROM SAF Workshop, ECMWF, 16–18 June 2014 1 Ian Culverwell, Met Office Sean Healy, ECMWF Progress (…) on forward modelling L1

ROM SAF Workshop, ECMWF, 16–18 June 2014 1

Ian Culverwell, Met Office

Sean Healy, ECMWF

Progress (…) on forward modelling L1 and L2 bending

angles

Page 2: Progress (…) on forward modelling L1 and L2 …...ROM SAF Workshop, ECMWF, 16–18 June 2014 1 Ian Culverwell, Met Office Sean Healy, ECMWF Progress (…) on forward modelling L1

ROM SAF Workshop, ECMWF, 16–18 June 2014 2

Motivation

• How to do ionospheric correction when L2 signal drops out?

• Different (better?) tropospheric retrievals if L1 and L2 both used?

• Some information on ionospheric parameters?

• Residual ionospheric error?

Page 3: Progress (…) on forward modelling L1 and L2 …...ROM SAF Workshop, ECMWF, 16–18 June 2014 1 Ian Culverwell, Met Office Sean Healy, ECMWF Progress (…) on forward modelling L1

ROM SAF Workshop, ECMWF, 16–18 June 2014 3

Model ionosphere: electron density

Single Chapman Layer (Chapman 1931)

ne(r) = nemax exp[½(1 – u – e-u)],

where u = (r – r0) / H.

3 parameters:

nemax = TEC / √(2πe) H

r0 = peak height

H = ionospheric scale height

Page 4: Progress (…) on forward modelling L1 and L2 …...ROM SAF Workshop, ECMWF, 16–18 June 2014 1 Ian Culverwell, Met Office Sean Healy, ECMWF Progress (…) on forward modelling L1

ROM SAF Workshop, ECMWF, 16–18 June 2014 4

Model ionosphere: bending angle

αLi(a) = – 2a ∫a∞ dlogn/dx dx / √(x2 – a2), x=nr

≈ (k4 / fLi2) 2a ∫a∞ dne/dx dx / √(x2 – a2), k4 = 40.3 m3s-2

≈ (k4 / fLi2) 2a [2r0 / (r0+a)3/2] ∫a∞ dne/dr dr / √(r – a)

= (k4 / fLi2) ne

max [4er02a2/H(r0+a)3]1/2 · Z(l)

where

Z(l) = ∫-l∞ (e-u – 1) exp[½(1 – u – e-u)] / √(u + l) du

is just a dimensionless, O(1) function of

l = (r0 – a) / H

Page 5: Progress (…) on forward modelling L1 and L2 …...ROM SAF Workshop, ECMWF, 16–18 June 2014 1 Ian Culverwell, Met Office Sean Healy, ECMWF Progress (…) on forward modelling L1

ROM SAF Workshop, ECMWF, 16–18 June 2014 5

Model ionosphere: bending angle Z(l)

Actual ≈

-√(2π) el/2

√(2π/l3)

Page 6: Progress (…) on forward modelling L1 and L2 …...ROM SAF Workshop, ECMWF, 16–18 June 2014 1 Ian Culverwell, Met Office Sean Healy, ECMWF Progress (…) on forward modelling L1

ROM SAF Workshop, ECMWF, 16–18 June 2014 6

Model ionosphere: bending angle Z(l)

Melbourne, 2004,

Gaussian ne

Schreiner et al, 1999,

Chapman layer

Hajj & Romans,

1998, GPS/MET

Page 7: Progress (…) on forward modelling L1 and L2 …...ROM SAF Workshop, ECMWF, 16–18 June 2014 1 Ian Culverwell, Met Office Sean Healy, ECMWF Progress (…) on forward modelling L1

ROM SAF Workshop, ECMWF, 16–18 June 2014 7

Direct modelling of L1 and L2 in ROPP

• The Radio Occultation Processing Package:

• A collection of Fortran 95 code, build and test scripts, data files and documentation designed to aid users wishing to process, quality-control and assimilate radio occultation data into their NWP models.

• Provided by ROM SAF (EUMETSAT).

• The following features will be available in ROPP8.0 (Dec 2014).

Page 8: Progress (…) on forward modelling L1 and L2 …...ROM SAF Workshop, ECMWF, 16–18 June 2014 1 Ian Culverwell, Met Office Sean Healy, ECMWF Progress (…) on forward modelling L1

ROM SAF Workshop, ECMWF, 16–18 June 2014 8

Examples of forward modelled L1 and L2 using ropp_fm

E-layer F-layer D-layer F1-layer

Page 9: Progress (…) on forward modelling L1 and L2 …...ROM SAF Workshop, ECMWF, 16–18 June 2014 1 Ian Culverwell, Met Office Sean Healy, ECMWF Progress (…) on forward modelling L1

ROM SAF Workshop, ECMWF, 16–18 June 2014 9

Direct modelling of L1 and L2 in ropp_1dvar

• Minimise 2J(x) = (x – b)T BBBB-1 (x – b) + (H(x) – o)T RRRR-1 (H(x) – o)

• x = {T, q, p*, nemax, r0, H}; o = {αL1, αL2}

• RRRR

• σ(αL1) = max(αn, 10 µrad); σ(αL2) = max(αn, 30 µrad)

• Assume αL1 and αL2 errors to be uncorrelated

• Needs some experimentation

• BBBB

• σ(nemax) = 2e11 m-3; σ(r0) = 150 km; σ(H) = 25 km

• Assume {nemax, r0, H} errors to be uncorrelated from each other and

from those of {T, q, p*}

• Needs some experimentation

Page 10: Progress (…) on forward modelling L1 and L2 …...ROM SAF Workshop, ECMWF, 16–18 June 2014 1 Ian Culverwell, Met Office Sean Healy, ECMWF Progress (…) on forward modelling L1

ROM SAF Workshop, ECMWF, 16–18 June 2014 10

Example 1: ropp_1dvar retrieval based on neutral bending angle

Page 11: Progress (…) on forward modelling L1 and L2 …...ROM SAF Workshop, ECMWF, 16–18 June 2014 1 Ian Culverwell, Met Office Sean Healy, ECMWF Progress (…) on forward modelling L1

ROM SAF Workshop, ECMWF, 16–18 June 2014 11

Example 1: ropp_1dvar retrieval based on L1 & L2 bending angles

Page 12: Progress (…) on forward modelling L1 and L2 …...ROM SAF Workshop, ECMWF, 16–18 June 2014 1 Ian Culverwell, Met Office Sean Healy, ECMWF Progress (…) on forward modelling L1

ROM SAF Workshop, ECMWF, 16–18 June 2014 12

Example 2: ropp_1dvar retrieval based on neutral bending angle

Page 13: Progress (…) on forward modelling L1 and L2 …...ROM SAF Workshop, ECMWF, 16–18 June 2014 1 Ian Culverwell, Met Office Sean Healy, ECMWF Progress (…) on forward modelling L1

ROM SAF Workshop, ECMWF, 16–18 June 2014 13

Example 2: ropp_1dvar retrieval based on L1 & L2 bending angles

Page 14: Progress (…) on forward modelling L1 and L2 …...ROM SAF Workshop, ECMWF, 16–18 June 2014 1 Ian Culverwell, Met Office Sean Healy, ECMWF Progress (…) on forward modelling L1

ROM SAF Workshop, ECMWF, 16–18 June 2014 14

Example 3: ropp_1dvar retrieval based on neutral bending angle

Page 15: Progress (…) on forward modelling L1 and L2 …...ROM SAF Workshop, ECMWF, 16–18 June 2014 1 Ian Culverwell, Met Office Sean Healy, ECMWF Progress (…) on forward modelling L1

ROM SAF Workshop, ECMWF, 16–18 June 2014 15

Example 3: ropp_1dvar retrieval based on L1 & L2 bending angles

Page 16: Progress (…) on forward modelling L1 and L2 …...ROM SAF Workshop, ECMWF, 16–18 June 2014 1 Ian Culverwell, Met Office Sean Healy, ECMWF Progress (…) on forward modelling L1

ROM SAF Workshop, ECMWF, 16–18 June 2014 16

Example 4: ropp_1dvar retrieval based on neutral bending angle

Page 17: Progress (…) on forward modelling L1 and L2 …...ROM SAF Workshop, ECMWF, 16–18 June 2014 1 Ian Culverwell, Met Office Sean Healy, ECMWF Progress (…) on forward modelling L1

ROM SAF Workshop, ECMWF, 16–18 June 2014 17

Example 4: ropp_1dvar retrieval based on L1 & L2 bending angles

Page 18: Progress (…) on forward modelling L1 and L2 …...ROM SAF Workshop, ECMWF, 16–18 June 2014 1 Ian Culverwell, Met Office Sean Healy, ECMWF Progress (…) on forward modelling L1

ROM SAF Workshop, ECMWF, 16–18 June 2014 18

Direct modelling of L1 and L2 in ROPP: αn stats

Page 19: Progress (…) on forward modelling L1 and L2 …...ROM SAF Workshop, ECMWF, 16–18 June 2014 1 Ian Culverwell, Met Office Sean Healy, ECMWF Progress (…) on forward modelling L1

ROM SAF Workshop, ECMWF, 16–18 June 2014 19

Direct modelling of L1 and L2 in ROPP: αLi stats

Page 20: Progress (…) on forward modelling L1 and L2 …...ROM SAF Workshop, ECMWF, 16–18 June 2014 1 Ian Culverwell, Met Office Sean Healy, ECMWF Progress (…) on forward modelling L1

ROM SAF Workshop, ECMWF, 16–18 June 2014 20

“Retrieved ionospheric parameters”

Fit “normal” L1-L2 to Z(l) for every GRAS occultation for 5 days in June 2012

Page 21: Progress (…) on forward modelling L1 and L2 …...ROM SAF Workshop, ECMWF, 16–18 June 2014 1 Ian Culverwell, Met Office Sean Healy, ECMWF Progress (…) on forward modelling L1

ROM SAF Workshop, ECMWF, 16–18 June 2014 21

Direct modelling of L1 and L2 in ROPP: summary

• Can forward model αL1 and αL2 by assuming simple Chapman layer ionosphere

• Allows extrapolation of αL1 and αL2 based on whole profiles of data

• Retrievals using αL1 and αL2 are possible in principle BUT

• Bigger impacts on q and tropospheric T?

• Not always capturing curvature of L1 and L2 at top?

• Require nemax < 0 to get decent fit in some cases!

• Backgrounds/errors/covariances of nemax, r0 and H

need some experimentation.

• Errors/covariances of αL1 and αL2 need some experimentation.

Page 22: Progress (…) on forward modelling L1 and L2 …...ROM SAF Workshop, ECMWF, 16–18 June 2014 1 Ian Culverwell, Met Office Sean Healy, ECMWF Progress (…) on forward modelling L1

ROM SAF Workshop, ECMWF, 16–18 June 2014 22

Residual ionospheric correction

• Working on a model for the residual error that remains after standard ionospheric correction (Vorob’ev and Krasil’nikova,1994)

• α(a) = – 2a ∫a∞dlogn/dx dx /√(x2-a2) ≈

– 2a ∫a∞N′ dr /√(r2 – a2) Basic bending angle, ∝∝∝∝ TEC/f2

+ 2a ∫a∞ NN′ dr (2r2 – a2)/(r2 – a2)3/2 1st order correction, ∝∝∝∝ TEC2/f4

• Implies that residual after standard ionospheric correction

∝∝∝∝ (αL1 – αL2)2 X some slowly varying function of height,

which we can estimate/calculate analytically for a Chapman layer.

• We are working with Julia Danzer (GFZ/ROM SAF VS) to test this.

• Could have implications for RO-based climatological mean temperatures.

Page 23: Progress (…) on forward modelling L1 and L2 …...ROM SAF Workshop, ECMWF, 16–18 June 2014 1 Ian Culverwell, Met Office Sean Healy, ECMWF Progress (…) on forward modelling L1

ROM SAF Workshop, ECMWF, 16–18 June 2014 23

ResidualαL1

αL1 correction

(cf climate

change trend

~ few µrad

per decade)

Page 24: Progress (…) on forward modelling L1 and L2 …...ROM SAF Workshop, ECMWF, 16–18 June 2014 1 Ian Culverwell, Met Office Sean Healy, ECMWF Progress (…) on forward modelling L1

ROM SAF Workshop, ECMWF, 16–18 June 2014 24

Conclusions

• Calculated bending angle produced by a spherically symmetric Chapman layer ionosphere.

• Incorporated in ROPP forward model and 1dvar system.

• Expressed residual bending angle, after standard ionospheric correction, in terms of (αL1 – αL2)

2.

Page 25: Progress (…) on forward modelling L1 and L2 …...ROM SAF Workshop, ECMWF, 16–18 June 2014 1 Ian Culverwell, Met Office Sean Healy, ECMWF Progress (…) on forward modelling L1

ROM SAF Workshop, ECMWF, 16–18 June 2014 25

Background error covariance matrix

Page 26: Progress (…) on forward modelling L1 and L2 …...ROM SAF Workshop, ECMWF, 16–18 June 2014 1 Ian Culverwell, Met Office Sean Healy, ECMWF Progress (…) on forward modelling L1

ROM SAF Workshop, ECMWF, 16–18 June 2014 26

Bending angle error covariance matrix

cf σnmin = 6 µrad

σ2min = 30 µrad

σ1min = 10 µrad

Page 27: Progress (…) on forward modelling L1 and L2 …...ROM SAF Workshop, ECMWF, 16–18 June 2014 1 Ian Culverwell, Met Office Sean Healy, ECMWF Progress (…) on forward modelling L1

ROM SAF Workshop, ECMWF, 16–18 June 2014 27

Fit “normal” L1-L2 to Z(l) for every GRAS occultation for

5 days in June 2012

Page 28: Progress (…) on forward modelling L1 and L2 …...ROM SAF Workshop, ECMWF, 16–18 June 2014 1 Ian Culverwell, Met Office Sean Healy, ECMWF Progress (…) on forward modelling L1

ROM SAF Workshop, ECMWF, 16–18 June 2014 28

Fit “normal” L1-L2 to Z(l) for every GRAS occultation for

5 days in June 2012