african workshop on gnss and space weather: introduction...

36
African Workshop on GNSS and Space Weather: Introduction to Total Electron Content (TEC) Anthea Coster, MIT Haystack Observatory Outline What is TEC and why is it important? How is TEC Measured? Where to obtain TEC? The Madrigal Database. Examples of how TEC measures Space Weather. Summary

Upload: others

Post on 07-Mar-2021

1 views

Category:

Documents


0 download

TRANSCRIPT

Page 1: African Workshop on GNSS and Space Weather: Introduction ...indico.ictp.it/event/9124/session/22/contribution/92/...B = geomagnet ic field, N = electron density e = electroni c charge,

African Workshop on GNSS and Space Weather: Introduction to Total Electron Content (TEC)

Anthea Coster , MIT Haystack Observatory

OutlineWhat is TEC and why is it important?How is TEC Measured?Where to obtain TEC? The Madrigal Database.Examples of how TEC measures Space Weather.Summary

Page 2: African Workshop on GNSS and Space Weather: Introduction ...indico.ictp.it/event/9124/session/22/contribution/92/...B = geomagnet ic field, N = electron density e = electroni c charge,

DEFINITION OF TEC

From Attila Komjathy, JPL

TEC = Total Electron Content (1016 x el/m2)

Page 3: African Workshop on GNSS and Space Weather: Introduction ...indico.ictp.it/event/9124/session/22/contribution/92/...B = geomagnet ic field, N = electron density e = electroni c charge,

Structure of the Ionosphere

D-region

F-region

Topside

E-region

The Ionosphere:Its regions and day/night transitions

Page 4: African Workshop on GNSS and Space Weather: Introduction ...indico.ictp.it/event/9124/session/22/contribution/92/...B = geomagnet ic field, N = electron density e = electroni c charge,

Map of GPS receivers

Page 5: African Workshop on GNSS and Space Weather: Introduction ...indico.ictp.it/event/9124/session/22/contribution/92/...B = geomagnet ic field, N = electron density e = electroni c charge,

Distributed vs. single point measurements

Wide Area Distribution of 'Raw' Information

Distributed networks of

sensors yield global

physics unattainable

with single-point

measurements

[Coster et al, 2003]

Example :

Global GPS-derived

ionospheric mapping

during geomagnetic

disturbances

Page 6: African Workshop on GNSS and Space Weather: Introduction ...indico.ictp.it/event/9124/session/22/contribution/92/...B = geomagnet ic field, N = electron density e = electroni c charge,

What is TEC and why is it important?How is TEC Measured?Where to obtain TEC? The Madrigal DatabaseExamples of TEC StudiesSummary

OUTLINE

Page 7: African Workshop on GNSS and Space Weather: Introduction ...indico.ictp.it/event/9124/session/22/contribution/92/...B = geomagnet ic field, N = electron density e = electroni c charge,

GPS Background

• At most 32 satellites

• 6 orbital planes

• 4~6 satellites per plane

• 55o inclination angle

• near circular orbit

• ~ 20000 km altitude

• ~12 hours round trip

(11 hour 58 min 2.05 sec)

Page 8: African Workshop on GNSS and Space Weather: Introduction ...indico.ictp.it/event/9124/session/22/contribution/92/...B = geomagnet ic field, N = electron density e = electroni c charge,

GPS Background

Each GPS spacecraft:

• Carries highly accurate clock

• Transmits its clock and position

• Signals are transmitted on 2 (or 3)

frequencies

Page 9: African Workshop on GNSS and Space Weather: Introduction ...indico.ictp.it/event/9124/session/22/contribution/92/...B = geomagnet ic field, N = electron density e = electroni c charge,

GPS Satellite Signal Structure

9

Carrier

• Doppler

• Range

Code modulation:

• Identifies SV

• Spread power

• Range

Navigation data

• SV orbit

• Error correction

• SV health

X

GPS

signal

Page 10: African Workshop on GNSS and Space Weather: Introduction ...indico.ictp.it/event/9124/session/22/contribution/92/...B = geomagnet ic field, N = electron density e = electroni c charge,

Range Measurements

10

Transmitted at SV:

tts tr

Received at RX:

)( sr ttc −=

Pseudorange:

Nsr +−= )(

2

Carrier phase:

sr

Page 11: African Workshop on GNSS and Space Weather: Introduction ...indico.ictp.it/event/9124/session/22/contribution/92/...B = geomagnet ic field, N = electron density e = electroni c charge,

GPS Positioning

Page 12: African Workshop on GNSS and Space Weather: Introduction ...indico.ictp.it/event/9124/session/22/contribution/92/...B = geomagnet ic field, N = electron density e = electroni c charge,

GPS Signal Multipath

GPS ANTENNA

Page 13: African Workshop on GNSS and Space Weather: Introduction ...indico.ictp.it/event/9124/session/22/contribution/92/...B = geomagnet ic field, N = electron density e = electroni c charge,

Illustration of Atmospheric Effects

Elevation Refraction Range Delay

Page 14: African Workshop on GNSS and Space Weather: Introduction ...indico.ictp.it/event/9124/session/22/contribution/92/...B = geomagnet ic field, N = electron density e = electroni c charge,

n = 2 11

11

2

1

412 4 2 2

12

-X X

X Y Y X YT T L

( )

( ) ( )

− − + −

Where:

XN

=

2

2 YH

=

N

e

Ne

m=

2

0

12

H

e

e B

m=

= the angular frequency of the radar wave,

Y = Y cos , Y = Ysin ,

= angle between the wave vector k and B

k = wave vector of propagating radiation,

B = geomagnetic field, N = electron density

e = electronic charge, m = electron mass,

and permittivity constant.

L T

e

0

,

=

Index of Refraction in the Ionosphere

Page 15: African Workshop on GNSS and Space Weather: Introduction ...indico.ictp.it/event/9124/session/22/contribution/92/...B = geomagnet ic field, N = electron density e = electroni c charge,

Ionospheric Range Correction

22

2

2

1

2

2

12

1) - 1( f

ANn eNN −−

Page 16: African Workshop on GNSS and Space Weather: Introduction ...indico.ictp.it/event/9124/session/22/contribution/92/...B = geomagnet ic field, N = electron density e = electroni c charge,

GPS Positioning

Page 17: African Workshop on GNSS and Space Weather: Introduction ...indico.ictp.it/event/9124/session/22/contribution/92/...B = geomagnet ic field, N = electron density e = electroni c charge,

TEC from GPS is measured from the difference of the GPS pseudo-range measurement at two frequencies

17

Where P1 and P2 are the pseudo-ranges measured by GPS at the two different frequencies, f1 and f2.

Page 18: African Workshop on GNSS and Space Weather: Introduction ...indico.ictp.it/event/9124/session/22/contribution/92/...B = geomagnet ic field, N = electron density e = electroni c charge,

Total Electron Content (TEC) EstimationDual-Frequency Measurements

Page 19: African Workshop on GNSS and Space Weather: Introduction ...indico.ictp.it/event/9124/session/22/contribution/92/...B = geomagnet ic field, N = electron density e = electroni c charge,

Ionospheric ParametersGPS can be used to measure

Ground-Based Receivers

• Total Electron Content

• Scintillation Parameters: S4 and σΦ

Space-Based Receivers

• Electron Density Profiles

• Scintillation Parameters: S4 and σΦ

Page 20: African Workshop on GNSS and Space Weather: Introduction ...indico.ictp.it/event/9124/session/22/contribution/92/...B = geomagnet ic field, N = electron density e = electroni c charge,

Problem 1: Mapping Function

• Mapping function used to map TEC (line of sight) to Vertical TEC

• Assigns V-TEC to pierce point

• Function of Elevation Range: 1 < Z < 3

Page 21: African Workshop on GNSS and Space Weather: Introduction ...indico.ictp.it/event/9124/session/22/contribution/92/...B = geomagnet ic field, N = electron density e = electroni c charge,

Problem 2: GPS Biases

• GPS delay difference between two frequencies provides TEC

• Delay differences are also introduced by the satellite and receiver

• Satellite biases are determined by IGS community and are fairly stable

• Receiver biases are determined by individual user and most users estimate one bias over a 24 - hour period.

Page 22: African Workshop on GNSS and Space Weather: Introduction ...indico.ictp.it/event/9124/session/22/contribution/92/...B = geomagnet ic field, N = electron density e = electroni c charge,

What is TEC and why is it important?How is TEC Measured?Where to obtain TEC? The Madrigal Database.Examples of TEC StudiesSummary

OUTLINE

Page 23: African Workshop on GNSS and Space Weather: Introduction ...indico.ictp.it/event/9124/session/22/contribution/92/...B = geomagnet ic field, N = electron density e = electroni c charge,

What is now available through CEDAR Madrigal Webhttp://cedar.openmadrigal.org/

Page 24: African Workshop on GNSS and Space Weather: Introduction ...indico.ictp.it/event/9124/session/22/contribution/92/...B = geomagnet ic field, N = electron density e = electroni c charge,

Standard TEC Data in Madrigal available since 2000

1. Provided in 1 degree by 1 degree bins2. Provided every 5 minutes3. Vertical TEC data estimates and Errors on these estimates4. Geographic Lat and Long5. Only provides data where observations are available. Does not attempt to model TEC where data is not available. Uses all GNSS data available.6. New TEC products are on the horizon …i.e. GLONASS observations

Page 25: African Workshop on GNSS and Space Weather: Introduction ...indico.ictp.it/event/9124/session/22/contribution/92/...B = geomagnet ic field, N = electron density e = electroni c charge,

25

Page 26: African Workshop on GNSS and Space Weather: Introduction ...indico.ictp.it/event/9124/session/22/contribution/92/...B = geomagnet ic field, N = electron density e = electroni c charge,

New Line of Site TEC Data in Madrigalavailable now for ~ 3 years

1. Provided for every receiver2. Provided every 20 second 3. Satellite and Receiver ID4. Geographic Lat and Long of Receiver5. Pierce Point: Altitude, Lat and Long6. Azimuth and Elevation to Satellite7. Files are LARGE8. HDF5 format

Page 27: African Workshop on GNSS and Space Weather: Introduction ...indico.ictp.it/event/9124/session/22/contribution/92/...B = geomagnet ic field, N = electron density e = electroni c charge,

27

Page 28: African Workshop on GNSS and Space Weather: Introduction ...indico.ictp.it/event/9124/session/22/contribution/92/...B = geomagnet ic field, N = electron density e = electroni c charge,

28

Instructions for using Madrigal provided at school website.

Page 29: African Workshop on GNSS and Space Weather: Introduction ...indico.ictp.it/event/9124/session/22/contribution/92/...B = geomagnet ic field, N = electron density e = electroni c charge,

What is TEC and why is it important?How is it MeasuredOne Place to obtain it: Madrigal DatabaseExamples of TEC StudiesSummary

OUTLINE

Page 30: African Workshop on GNSS and Space Weather: Introduction ...indico.ictp.it/event/9124/session/22/contribution/92/...B = geomagnet ic field, N = electron density e = electroni c charge,

Astronomy Picture of the Day 17 September 2012

At the end of last month, a long standing solar filament suddenly erupted into space

producing an energetic Coronal Mass Ejection (CME)

Coronal Mass Ejection

Page 31: African Workshop on GNSS and Space Weather: Introduction ...indico.ictp.it/event/9124/session/22/contribution/92/...B = geomagnet ic field, N = electron density e = electroni c charge,

Day 90, 2001 Day 101, 2001

Page 32: African Workshop on GNSS and Space Weather: Introduction ...indico.ictp.it/event/9124/session/22/contribution/92/...B = geomagnet ic field, N = electron density e = electroni c charge,

Venetie, Alaska 1 March 2017

32

Page 33: African Workshop on GNSS and Space Weather: Introduction ...indico.ictp.it/event/9124/session/22/contribution/92/...B = geomagnet ic field, N = electron density e = electroni c charge,

33

10:03 UT

Page 34: African Workshop on GNSS and Space Weather: Introduction ...indico.ictp.it/event/9124/session/22/contribution/92/...B = geomagnet ic field, N = electron density e = electroni c charge,

Example of Differential TEC showing response to Solar Flare X9.3 on 6 Sep 2017

Page 35: African Workshop on GNSS and Space Weather: Introduction ...indico.ictp.it/event/9124/session/22/contribution/92/...B = geomagnet ic field, N = electron density e = electroni c charge,

16 YR Median Diff of Daytime GPS and ISR TEC at Millstone from 2006-2017. Median Diff. = 8.97 STD = 4.6

Page 36: African Workshop on GNSS and Space Weather: Introduction ...indico.ictp.it/event/9124/session/22/contribution/92/...B = geomagnet ic field, N = electron density e = electroni c charge,

36

SUMMARY

• Global TEC maps can be derived from multiple GNSS (GPS) receivers

• This data is available through Madrigal database

• TEC is a powerful tool to monitor space physics events: geomagnetic storms, solar flares, auroras