M. Fredrizzi, et al. 2009 I NTRODUCTION , CONT . 4 August 2010 My - - PowerPoint PPT Presentation

m fredrizzi et al 2009 i ntroduction cont
SMART_READER_LITE
LIVE PREVIEW

M. Fredrizzi, et al. 2009 I NTRODUCTION , CONT . 4 August 2010 My - - PowerPoint PPT Presentation

4 August 2010 M ODELING THE D ENSITY OF THE T HERMOSPHERE Suzanne Smith 1 Mentor: Tomoko Matsuo Site: National Oceanic & Atmospheric Administration, NOAA E ARTH S A TMOSPHERE 4 August 2010 2 I MPORTANCE OF M ODELING THE T HERMOSPHERE


slide-1
SLIDE 1

MODELING THE DENSITY OF

THE THERMOSPHERE

Suzanne Smith Mentor: Tomoko Matsuo Site: National Oceanic & Atmospheric Administration, NOAA

4 August 2010

1

slide-2
SLIDE 2

EARTH’S ATMOSPHERE

4 August 2010

2

slide-3
SLIDE 3

IMPORTANCE OF MODELING THE THERMOSPHERE

 Height of satellites and space shuttles orbit.  The neutral density of the thermosphere effects the

amount of drag present.

 With increased density and drag the shuttles and

satellites are slowed and the orbiting altitude is decreased.

 Having an efficient and accurate model of thermospheric

density is a valuable asset.

4 August 2010

3

slide-4
SLIDE 4

INTRODUCTION

 General Circulation Model, GCM  Previous work  CTIPe model: The Coupled Thermosphere Ionosphere

Plasmasphere Electrodynamics Model, Tim Fuller-Rowell et al. 1996

 Global Thermosphere 80-500km: solves momentum, energy, composition  Ionosphere 80-10,000km: solves continuity, momentum, energy, etc.  Forcing: solar UV and EUV, empirical high latitude electric field and

auroral precipitation models, tidal forcing.

 CHAMP Satellite: Challenging Minisatellite Payload Satellite

 height~ 400km; 90min orbital period; Launched date: July 2000.

 2005 CTIPe 5-min Run, Mariangel Fedrizzi

4 August 2010

4

slide-5
SLIDE 5

4 August 2010

5

  • M. Fredrizzi, et al. 2009
slide-6
SLIDE 6

INTRODUCTION, CONT.

 My Work  Used multi-dimensional GCM (CTIPe) output and

reduced it to a low-dimension model.

 Specifically, conducted Singular Value Decomposition

(SVD) Analysis of CTIPe 5-min model output from 2005, and constructed a model of thermospheric density.

 Density in terms of position and time:  (r, t) = 1(r) 1(t) + 2(r) 2(t) + .. + n(r) n(t)  n(r) = EOF  n(t) = Amplitude

4 August 2010

6

slide-7
SLIDE 7

DRIVERS OF DENSITY CHANGE

 Extreme Ultra Violet(EUV)  Diurnal  Seasonal  Solar wind/Magnetosphere Interactions  Auroral Activity

4 August 2010

7

slide-8
SLIDE 8

YEAR MEAN & EOF AMPLITUDE VARIANCE

4 August 2010

8

slide-9
SLIDE 9

YEARS WORTH OF EMPIRICAL ORTHOGANAL FUNCTIONS (EOFS), 

4 August 2010

9

slide-10
SLIDE 10

EOF AMPLITUDES, 

4 August 2010

10

slide-11
SLIDE 11

MODE #1: DIURNAL EUV

 Caused by the earth’s daily rotation.  The day side’s density increases because of the increased

EUV.

4 August 2010

11

August 2005

slide-12
SLIDE 12

MODE #2: SEASONAL EUV

4 August 2010

12

slide-13
SLIDE 13

MODE #2: SEASONAL EUV CONT.

 Caused by the earth’s yearly revolution around the sun.  In our summer months the northern hemisphere is

pointed towards the sun which results in a greater amount of EUVs.

4 August 2010

13

Summer ‘05 Winter‘05

slide-14
SLIDE 14

MODE #3: AURORAL ACTIVITY

 Cause by high latitude electromagnetic forcing resulted

from the interaction between Solar Wind and the earth’s magnetosphere (i.e., auroral activity).

 Aurora occur both in the Northern and Southern

hemisphere creating a symmetric pattern in the EOF contour plots.

4 August 2010

14

Oct 2005

slide-15
SLIDE 15

RESOURCES: DRIVERS OF DENSITY CHANGE

 Ap Index (Kyoto): A measure of the level of

geomagnetic activity over the globe taken every 3hrs.

 Solar Wind (NASA OMNIWeb): collection of different

data sets that help to display storm conditions.

 Joule Heating (CTIPe Model): integrated over the globe  F 10.7 (Ottawa 10.7cm flux): EUV index

4 August 2010

15

slide-16
SLIDE 16

PROVING MODE #3 IS AURORAL ACTIVITY

4 August 2010

16

F10.7 Ap Ap > 150 EOF 1 0.5163 0.4411 0.3068 EOF 2 0.0410 0.0345 0.2714 EOF 3 0.0388 0.0097 0.2548 EOF 4 0.6221 0.0364 0.3659

 Correlating the different EOFs with EUV Index: F10.7

(daily value), and Geomagnetic Index: Ap (taken every three hours).

 Surprising lack of correlation between Ap and EOF3.

slide-17
SLIDE 17

AUGUST 24TH: SOLAR WIND DATA(OMNI)

4 August 2010

17

Proton Density Dst Index Ap Index Solar Wind Speed B-field

slide-18
SLIDE 18

AUGUST 24TH: AP INDEX & JOULE HEATING

4 August 2010

18

slide-19
SLIDE 19

AUGUST 24TH: THERMOSPHERIC DENSITY RECONSTRUCTION USING EOFS

4 August 2010

19

slide-20
SLIDE 20

FINISHED PRODUCT

4 August 2010

20

slide-21
SLIDE 21

ACKNOWLEDGEMENTS

 Tomoko Matsuo, mentor  Mariangel Fredrizzi, officemate & CTIPe Data  Timothy Fuller-Rowell, CTIPe model & mentoring  Rodney, dark chocolate covered acia berries  Doug Biesecker  Mike Crumly, vouching for me  Russ Henson, technology help  National Oceanic & Atmospheric Administration, NOAA  Space Weather Prediction Center, SWPC  MatLab

4 August 2010

21

slide-22
SLIDE 22

REFERENCE

 NOAA Crest,

http://www.thebradentontimes.com/clientuploads/webpages/n

  • aa-logo.jpg.

 Lycoming Crest,

http://upload.wikimedia.org/wikipedia/en/thumb/1/1d/Lycomi ng_College_logo.png/175px-Lycoming_College_logo.png.

 Earth’s Atmosphere,

http://www.vtaide.com/png/images/atmosphere.jpg.

 CHAMP & CTIPe data plot, Mariangel Fredrizzi, et al.  Ap Index, http://wdc.kugi.kyoto-u.ac.jp/kp/index.html.  Solar Wind Data, NASA OMNIWeb,

http://omniweb.gsfc.nasa.gov/

 Joule Heating, CTIPe Model  F 10.7, Daily F 10.7 index, the Ottawa 10.7cm (2800 MHz)

radio flux

4 August 2010

22

slide-23
SLIDE 23

QUESTIONS?

4 August 2010

23