ARTS - GARLIC - KOPRA

Transcription

1 ARTS - GARLIC - KOPRA A Lbl Model Intercomparison Franz Schreier, Sebastián Gimeno García, Mathias Milz, Ajil Kottayil.. DLR Remote Sensing Technology Institute, Oberpfaffenhofen, GERMANY. LTU Lulea University of Technology, Kiruna, SWEDEN ARTS Workshop Remote Sensing Technology Institute ARTS - GARLIC - KOPRA

2 Introduction Motivation Verification and validation: one of the most important steps of code development Quality of radiative transfer models critical for atmospheric remote sensing products Dozen(s) of RT model intercomparisons, incl. Lbl models, e.g. AMILDA: KOPRA vs. GARLIC and some more IRTMW: ARTS vs. GARLIC and some others Kristineberg 9: ARTS vs. KOPRA ARTS GARLIC KOPRA?!?!? Remote Sensing Technology Institute ARTS - GARLIC - KOPRA

3 Outline The codes Lbl Intercomparisons Data and Assumptions Results Remote Sensing Technology Institute ARTS - GARLIC - KOPRA

4 ARTS Atmospheric Radiative Transfer Simulator Originally developed for microwave applications Nadir: AMSU-B, MHS, HIRS,... Limb: Odin, SMILES,... Uplooking: MIAWARA, AMSOS,... Continuum: MT-CKD. S. Bühler, P. Eriksson et al., JQSRT 9, 65 9, 5 P. Eriksson, S. Bühler et al., JQSRT, , Remote Sensing Technology Institute ARTS - GARLIC - KOPRA

5 GARLIC Generic Atmospheric Radiation Lbl Infrared Code Line shapes: Voigt, VanVleck Doppler, Lorentz Line data: HITRAN, HITEMP, GEISA, JPL,... Continua: Geometries: CKD. (H O, CO, N, O ), dry air (Liebe) Limb, uplooking, downlooking (refraction optional) Instruments: Spectral response: FTS, Heterodyne, Fabry Perot,... Field-of-view: Box, Gauss, Trapez,... Implementation: FORTRAN 8, all data read from external files Inversion: Extensions: Jacobians by automatic differentation (multiple) scattering infrared radiative transfer: J. Mendrok: SARTRE approx. spherical geo (GRL 7) M. Vasquez: cloudy exo-planet atmospheres (A&A a,b) F. Schreier et al. JQSRT, 7, 9 5, Remote Sensing Technology Institute ARTS - GARLIC - KOPRA 5

6 KOPRA Karlsruhe Optimized & Precise Radiative transfer Algorithm Originally developed for mid infrared Fourier transform spectroscopy Limb: MIPAS-Envisat, MIPAS-balloon Limb+uplooking: MIPAS-STR, GLORIA Uplooking, solar absorption: FTIR Kiruna non-lte and (single) scattering Continuum: H O: CKD., CO :??, N : Lafferty AO 96, O : Tibault 97 G. Stiller, T. von Clarmann, B. Funke, N. Glatthor, F. Hase, M. Höpfner and A. Linden, JQSRT 7, 9 8, Remote Sensing Technology Institute ARTS - GARLIC - KOPRA 6

7 Outline The codes Lbl Intercomparisons Data and Assumptions Results Remote Sensing Technology Institute ARTS - GARLIC - KOPRA 7

8 AMILDA Advanced MIPAS Level Data Analysis Intercomparison of algorithms and data analysis strategies used by five European groups Series of exercises: from simple cell transmission to more realistic radiance Exercise : Limb with km, apodized FTS line shape, finite field-of-view; H O, CO, O, N O, CH ; HITRAN 98 and CKD-continuum T. von Clarmann et al. JQSRT, 78, 8 7, I [ nw/(cm sr cm ) ] Radiance I [ nw/(cm sr cm ) ] BexpOD trapez-path trapez-laguerre,5.5,6,6.5 KOPRA GARLIC,5,5.5,6,6.5,7 Wavenumber ν [cm ] Remote Sensing Technology Institute ARTS - GARLIC - KOPRA 8

9 IRTMW Intercomparison of 8 µw codes Series of 5 cases, starting with line shape and absorption Full RT incl. FoV and SRF all geometries Case Uplooking RAM-like Temperature TB [K] TB [K] TB [K].5... ARTS dg ARTS 6dg GARLIC dg GARLIC 6dg avint BexpOD Hermite trapez-path trapez-laguerre Wavenumber ν [cm ] avint BexpOD Hermite trapez-path trapez-laguerre Case Downlooking AMSU-like Equivalent Temperature T B [K] T B [K] T B [K] ARTS 8dg 69 ARTS dg GARLIC TL 8dg GARLIC TLdg avint B exp OD Hermite path trapez trapez Laguerre Wavenumber [cm ] C. Melsheimer et al. Radio Science,, RS7, 5 Remote Sensing Technology Institute ARTS - GARLIC - KOPRA 9

10 Kristineberg 9 M. Milz: ARTS IR Application Radiances for one channel Motivation: Nadir sounding with similar geometry as AMSU-B / MHS since early 8 s Conclusion: Good agreement ARTS KOPRA (and LblRTM) Outlook: Update continuum HIRS Channel #7 for all Profiles HIRS Channel #9 for all Profiles HIRS Channel # for all Profiles 8 6 HIRS Channel #8 for all Profiles HIRS Channel # for all Profiles HIRS Channel # for all Profiles Remote Sensing Technology Institute ARTS - GARLIC - KOPRA

11 Outline The codes Lbl Intercomparisons Data and Assumptions Results Remote Sensing Technology Institute ARTS - GARLIC - KOPRA

12 GARAND atmospheres Garand Atmospheres 6 5 Altitude [km] 5 Temperature [K] HO [ppm] 6 8 O [ppm] L. Garand et al. Remote Sensing Technology Institute ARTS - GARLIC - KOPRA

13 HIRS High Resolution Infrared Sounder 9 infrared channel scanning radiometer used for operational atmospheric sounding by NASA/NOAA and EUMETSAT. CO absorption bands for temperature sounding. Also measures water vapour, ozone, N O and cloud and surface temperatures. Radiance [µw / (cm sr cm )] 5 Transmission Response Transmission HO O CO CH CO NOAA7! wavenumber [cm ] Remote Sensing Technology Institute ARTS - GARLIC - KOPRA

14 Further Assumptions HITRAN, main molecules only Continuum: CKD or MT-CKD Voigt line shape (no prefactor) Geometry: nadir (downlooking 8 ), no FoV Surface temperature T surf = T atm ( km) Remote Sensing Technology Institute ARTS - GARLIC - KOPRA

15 Outline The codes Lbl Intercomparisons Data and Assumptions Results Remote Sensing Technology Institute ARTS - GARLIC - KOPRA 5

16 Brightness Temperatures 8 tropical 8 mls 8 mlw arts garlic kopra T B [K] us sas saw T B [K] Remote Sensing Technology Institute Channel # ARTS - GARLIC - KOPRA 6 Fri May :8:

17 Perfect! Thanks for your attention Remote Sensing Technology Institute ARTS - GARLIC - KOPRA 7

18 Perfect! wait a moment whats about the other atmospheres?... and differences etc.??? Remember: IASI etc. have sub-kelvin sensitivity! Remote Sensing Technology Institute ARTS - GARLIC - KOPRA 7

19 First results IRS T [K] Ch # 67 Ch # 5 76 Ch # 8 Ch # Ch # 68 Ch # 6 7 Ch # 66 Ch #6 Ch # 69 Ch # 7 79 Ch # 59 Ch # n Ch # 7 Ch # 8 9 Ch # 87 Ch #8 5 a g g k k a Ch # 9 Ch # Ch #9 66 Remote Sensing Technology Institute ARTS - GARLIC - KOPRA 8 Fri Feb 6:58:

20 : N + O lines & continuum added in GARLIC T [K] Ch # 67 Ch # 5 76 Ch # 8 Ch #5 5 Ch # 68 Ch # 6 7 Ch # 66 Ch #6 Ch # 69 Ch # 7 79 Ch # 59 Ch #7 9 n Ch # 7 Ch # 8 9 Ch # 87 Ch #8 5 a g g k k a Ch # 9 Ch # Ch #9 66 Remote Sensing Technology Institute ARTS - GARLIC - KOPRA 9 Fri Feb 6:7:8

21 Channel : ARTS GARLIC Channel Channel Similar absorption: CO > H O > O > N O Mostly stratospheric emission Continuum & H O irrelevant Q-branch, convolution,...?? Remote Sensing Technology Institute ARTS - GARLIC - KOPRA

22 Channel : ARTS GARLIC Channel Channel Similar absorption: CO > H O > O > N O Mostly stratospheric emission Continuum & H O irrelevant Q-branch, convolution,...?? NOAA5-HIRS for KOPRA!!! Response.. HIRS Normalized Response Functions noaa5/ch noaa5/ch noaa5/ch noaa5/ch noaa5/ch5 noaa5/ch6 noaa7/ch noaa7/ch noaa7/ch noaa7/ch noaa7/ch5 noaa7/ch Wavenumber [cm ] Fri Jun 6 :: Remote Sensing Technology Institute ARTS - GARLIC - KOPRA

23 Channel 5, 6, and 7: GARLIC KOPRA CO H O > O > N O Mostly free troposphere emission (Weighting functions peak at 8 and km for 78 and 75 cm ) Remote Sensing Technology Institute ARTS - GARLIC - KOPRA

24 Channel 8 and 9: ARTS KOPRA Weak absorption transparent down to troposphere Continua and/or H O?? Remote Sensing Technology Institute ARTS - GARLIC - KOPRA

25 Without continuum? ARTS vs. GARLIC ARTS GARLIC: no Continuum Differences mostly <.5 K What is special with atm # Channels: 8, 9,, 7, 8, 9 (9,, 8, 9, 5, 66 cm ) What is special with channel 59 cm? ARTS GARLIC: with Continuum.5 T [K] T [K] Atmosphere # Atmosphere # Thu Jun 5 7:5:9 Thu Jun 5 7::9 Remote Sensing Technology Institute ARTS - GARLIC - KOPRA

26 Without continuum? GARLIC vs. KOPRA Garand Atmospheres Good agreement (.9 K) in channel KOPRA failed for atm #,, 5, 5 z[km] GARLIC KOPRA: no Continuum 6 8 T[K] Thu Jun 5 :: 5 Garand Atmospheres T [K] 5 5 z [km] 6 Atmosphere # Thu Jun 5 7::9 e 6 e 5 e e e HO [ppv] Thu Jun 5 :58:6 Remote Sensing Technology Institute ARTS - GARLIC - KOPRA

27 What s common and what s different? ARTS GARLIC KOPRA S(T ) conversion semi-tips Atmos TIPS Weak lines? all rejection? Line wings 5 cm??? cm, H O: 5 cm Continuum MT-CKD. CKD. CKD.,... Inhomogeneous Path? quadrature B linear τ layers?? Remote Sensing Technology Institute ARTS - GARLIC - KOPRA 5

28 Conclusions and Outlook max T ka =. K in channel 5 cm (channel ignored) No code is the best for all channels and/or atmospheres max Tag = +. K in 66 cm max Tgk = +.5 K in channel 5 cm Remote Sensing Technology Institute ARTS - GARLIC - KOPRA 6

29 Conclusions and Outlook max T ka =. K in channel 5 cm (channel ignored) No code is the best for all channels and/or atmospheres max Tag = +. K in 66 cm max Tgk = +.5 K in channel 5 cm Use another code (e.g. LblRTM) as reference? NO! Include more codes? hmmmm? Use common absorption coefficients Plot differences vs. (min, max, delta,... ) temperature, H O,... Compare monochromatic spectra and/or transmission Open questions Continuum Lines strength T conversion Line wing cut-off (CO at cm ) Weak line suppression Remote Sensing Technology Institute ARTS - GARLIC - KOPRA 6

30 Outline 5 Appendix More T Difference Plots Weighting and Response Functions Continua Remote Sensing Technology Institute ARTS - GARLIC - KOPRA 7

31 Max. Temperature Delta in Layer T [K] Ch # 67 Ch # 5 76 Ch # 8 Ch # Ch # 68 Ch # 6 7 Ch # 66 Ch # Ch # 69 Ch # 7 79 Ch # 59 Ch # T max [K] Ch # 7 Ch # 8 9 Ch # 87.5 Ch # a g g k k a Ch # 9 Ch # Ch # Remote Sensing Technology Institute ARTS - GARLIC - KOPRA 8 Tue Feb 6 7::56

32 T [K] Brightness Temperature Difference vs Minimum Temperature Ch # 67. Ch # 5 76 Ch # 8 Ch # Ch # 68 Ch # 6 7 Ch # 66 Ch # Ch # 69 Ch # 7 79 Ch # 59 Ch #7 9.5 Ch # 87 Ch # Ch #8 5 Ch # T min [K].5.5 Ch # 7 Ch # a g g k k a Ch # 9 Remote Sensing Technology Institute ARTS - GARLIC - KOPRA 9 Mon Feb 5 :5:

33 T [K] Brightness Temperature Differences vs. Surface Temperature Ch # 67 Ch # 5 76 Ch # 8 Ch # Ch # 68 Ch # 6 7 Ch # 66 Ch # Ch # 69 Ch # 7 79 Ch # 59 Ch # T surf [K].5.5 Ch # 7 Ch # Ch # 87 Ch #8 5.5 a g g k k a Ch # 9 Ch # Ch #9 66 Remote Sensing Technology Institute ARTS - GARLIC - KOPRA Mon Feb 5 5:8:

34 Brightness Temperature Differences vs Temperature Range T [K] Ch # 67 Ch # 5 76 Ch # 8 Ch # Ch # 68 Ch # 6 7 Ch # 66 Ch # Ch # 69 Ch # 7 79 Ch # 59 Ch # T max T min [K].5.5 Ch # 7 Ch # Ch # 87 Ch #8 5.5 a g g k k a Ch # 9 Ch # Ch #9 66 Remote Sensing Technology Institute ARTS - GARLIC - KOPRA Mon Feb 5 5:5:

35 T [K] Ch # 67. Ch # 5 76 Ch # 8 Ch #5 5 6 Brightness Temperature Difference vs Water Column... Ch # 68 Ch # 6 7 Ch # 66 Ch # Ch # 69 Ch # 7 79 Ch # 59 Ch #7 9 6 H O [kg/m ] Ch # 87 Ch # Ch #8 5 Ch #9 66 Ch # 7 Ch # a g g k k a Ch # 9 6 Remote Sensing Technology Institute ARTS - GARLIC - KOPRA Mon Feb 5 ::6

36 Weighting Functions s I(ν, s) = B(ν, T ) T (ν, s ) s ds Remote Sensing Technology Institute ARTS - GARLIC - KOPRA

37 Response Functions HIRS Normalized Response Functions HIRS Normalized Response Functions.6 noaa5/ch7 noaa5/ch8 noaa5/ch9 noaa5/ch noaa7/ch7 noaa7/ch8 noaa7/ch9 noaa7/ch. noaa5/ch noaa5/ch noaa7/ch noaa7/ch Response. Response Wavenumber [cm ] HIRS Normalized Response Functions Fri Jun 6 :5:5 5 6 Wavenumbers [cm ] HIRS Normalized Response Functions Fri Jun 6 :7:. noaa5/ch noaa5/ch noaa5/ch5 noaa5/ch6 noaa7/ch noaa7/ch noaa7/ch5 noaa7/ch6.. noaa5/ch7 noaa5/ch8 noaa5/ch9 noaa7/ch7 noaa7/ch8 noaa7/ch9. Response. Response... 5 Wavenumber [cm ] Fri Jun 6 :: Wavenumber [cm ] Fri Jun 6 ::9 Remote Sensing Technology Institute ARTS - GARLIC - KOPRA

38 CO Continua Remote Sensing Technology Institute ARTS - GARLIC - KOPRA 5