Van Allen Probes データの解析

平成25年度国立極地研究所研究集会「太陽−地球大気の地上多点観測データ総合解析ワークショップ」
2013年8月21日

於 国立極地研究所

桂華邦裕
名古屋大学太陽地球環境研究所

1Outline• Brief overview of spacecraft and data
• Status and data access for each instrument• Data handling with TDAS
✓ EMFISIS wave spectra & background B✓ ECT/MagEIS

2

The Van Allen Probes
Launch: 04:05 EDT on Aug. 30, 2012
Commissioning phase: through Oct. 27, 2012
Science operation: after Oct. 27, 2012
• Orbit: ~600 km x 5.8 Re• Inclination: ~10 deg., |MLAT| < 20˚• Lapping rate: 4 – 5 laps/yearMore detail at– http://nasa.gov/rbsp, http://vanallenprobes.jhuapl.edu/index.php, YouTube– Agendas and presentations (incl. focused science
topics & inter-mission collaborations) at SWG
meetings available at the JHU/APL website.
– Instrument papers are available for FREE in Space
Science Reviews.
3

4

RBSP instruments[Mauk 2009 AGU talk]•Measurements:– Ions (20 eV to 75 MeV, withcomposition 20 eV to 300 keV)
– Electrons (20 eV to >10 MeV +50 to 2000 MeV)
– DC electric (to 10 Hz) andmagnetic (to 30 Hz) fields– Wave E and B fields (10 Hz to 12kHz, 1-D E to 400 kHz)5

B vector in UVW

B vector inGSE, GSM, SM6

EMFISIS: Status• Attitude solution issues

http://emfisis.physics.uiowa.edu/✓ nutation after maneuvers (280-300s0.003Hz)✓ eclipse ephemeris errors- only in geophysical coordinates +despun✓ 1-sec steps in ephemeris (fixed bysmoothing)• Problems that are ours to own✓ 14.5 Hz Heater signal on RBSP-A- well below calibration goals should not bother anyone otherthan visually.
✓ calibration (spin tones - fixed) about
to see a big re-processing- We will reprocess ALL data backto launch.
7

EMFISIS: Background B field(TDAS bleeding edge)• ; - - - EMFISIS LEVEL 3 (cadence=“4sec”, “1sec”, and “hires”(up to 64Hz)) - - IDL> timespan, '2013-03-17', 1, /days
IDL> cdf_leap_second_init; Proper inclusion of leap seconds
E.g., IDL> rbsp_load_emfisis, level='l3', cadence='1sec', coord='sm', probe='b'E.g., IDL> rbsp_load_emfisis, level='l3', cadence='4sec', coord='gsm', probe='b'
IDL> options, 'rbspb_emfisis_l3_1sec_sm_Mag', 'colors', [2,4,6]
IDL> options, 'rbspb_emfisis_l3_1sec_sm_Mag', 'labels', ['Bx','By','Bz']• OR After downloading L2 data manually,
IDL> fname_emfisis = 'rbsp-b_magnetometer_4sec-sm_emfisisL3_20121123_v1.2.2'
IDL> cdf_leap_second_init; Proper inclusion of leap seconds is important.
IDL> cdf2tplot, file=fname_emfisis
8

EMFISIS: Data access (Wave spectra CDF files)

PNG Walk (Java)with links toautoplot

9

EMFISIS: Wave power spectra HFR(TDAS bleeding edge)#Download CDF files manually.
fname_emfisis = ‘YOUR_DIRECTORY/rbsp-a_HFR-spectra_emfisis-L2_20130317_v1.2.4.cdf”
cdf_leap_second_init
cdf2tplot, files=fname_emfisis, /get_support_data
get_data, 'HFR_Spectra', data=data, dlim=dlim
store_data, 'HFR_Spectra', data={x:data.x,y:data.y,v:reform(data.v[0,*])}
options, 'HFR_Spectra', 'ylog', 1
options, 'HFR_Spectra', 'zlog', 1options, 'HFR_Spectra', 'ytitle', 'Freq [Hz]'options, 'HFR_Spectra', 'ysubtitle', ''options, 'HFR_Spectra', 'ztitle', dlim.ysubtitle
tplot, 'HFR_Spectra', title='Van Allen Probes A: HFR Spectra'
10

EMFISIS: Wave power spectra HFR
2013/07/09-11 storm#Download CDF files manually.
fname_emfisis = 'YOUR_DIRECTORY/rbsp-'+probe+'_HFR-'+datatype+'_emfisis-L2_201307'+['09','10','11']+'_v1.3.'+['3','2','4']+'.cdf'
cdf_leap_second_init
cdf2tplot, files=fname_emfisis, /get_support_data
get_data, 'HFR_Spectra', data=data, dlim=dlim
store_data, 'HFR_Spectra', data={x:data.x,y:data.y,v:reform(data.v[0,*])}
options, 'HFR_Spectra', 'ylog', 1
options, 'HFR_Spectra', 'zlog', 1
options, 'HFR_Spectra', 'ytitle', 'Freq [Hz]'options, 'HFR_Spectra', 'ysubtitle', ''options, 'HFR_Spectra', 'ztitle', dlim.ysubtitle
timespan, ‘2013-07-09’, 3
tplot, 'HFR_Spectra', title='Van Allen Probes A: HFR Spectra'12

:w

EMFISIS: Wave power spectra WFR
2013/07/09-11 storm#Download CDF files manually.
fname_emfisis = 'YOUR_DIRECTORY/rbsp-'+probe+'_WFR-spectral-matrix-diagonal_emfisis-L2_201307' $+ ['09','10','11']+'_v1.3.'+['3','2','4']+'.cdf'
cdf_leap_second_init
cdf2tplot, files=fname_emfisis, /get_support_data
tvar = ‘BuBu’
get_data, tvar, data=data, dlim=dlim
store_data, 'WFR_'+tvar, data={x:data.x,y:data.y,v:reform(data.v[0,*])}
options, 'WFR_'+tvar, 'spec', 1options, 'WFR_'+tvar, 'ylog', 1options, 'WFR_'+tvar, 'zlog', 1
options, 'WFR_'+tvar, 'ytitle', tvar
options, 'WFR_'+tvar, 'ysubtitle', 'Freq [Hz]'options, 'WFR_'+tvar, 'ztitle', dlim.ysubtitle
timespan, ‘2013-07-09’, 3
tplot, 'WFR_B?B?', title='Van Allen Probes A: HFR Spectra'
14

ECT: Status
Operating well

• HOPE-Lev2:✓ x6.7 code corrected, reprocessing waiting on metadata fixes.

• HOPE-Lev3:✓ PA code being improved to allow for longer integration
times and hope moments' code

• REPT-Lev3:✓ PA files validated.

• MagEIS-Lev3:✓ delivered for implementation in processing chain

16

ECT: Data access (QL & autoplot)

17

ECT: Data access(autoplot)

18

ECT: Data access (TDAS)• The PI team does not officially plan todevelop TDAS routines for ECT data.

• but existing TDAS routines (e.g.,cdf2tplot.pro) do work well on ECT CDF
files.

• In this workshop, we will look at ECT/MagEIS
electron and proton data.

19

- FESA > Spin-Averaged Differential Electron Flux
- FPSA > Spin-Averaged Differential Proton Flux
- FEDU > Undirectional Differential Electron Flux
- FPDU > Undirectional Differential Proton Flux•••

timespan, '2013-03-17', 1•
•
•
•
•
•
•
•
•
•
•
•
•
•
•

prefix='rbsp_'+probe+'_mageis'+'_'

ECT: Data accessprobe=’a’fname_mageis = YOUR_DIRECTORY + 'rbsp'+probe+'_ect-mageisL2_20130317_v2.2.0.cdf'
cdf2tplot, files=fname_mageis, /get_support_data, prefix=prefix; - - - OPTIONS - - tvar = prefix + ['FESA','FPSA']
options, tvar, 'ztitle', 'Flux [cm!U-2!Ns!U-1!Nsr!U-1!NkeV!U-1!N]'
options, tvar, 'zlog', 1
options, tvar, 'ylog', 1
options, tvar, 'ytitle', 'Energy [keV]'
options, tvar, 'ysubtitle', ''options, prefix+'FESA', 'yrange', [30,5000]
options, prefix+'FESA', 'ystyle', 1
options, prefix+'FPSA', 'yrange', [50,2000]
options, prefix+'FPSA', 'ystyle', 1; - - - TPLOT - - tplot, '*FESA *FPSA'