03.01.2015 Views

The SDO Flare Detective - Solar Physics at MSU

The SDO Flare Detective - Solar Physics at MSU

The SDO Flare Detective - Solar Physics at MSU

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

<strong>The</strong> AIA flare detective<br />

Paolo C. Grigis<br />

and the FFT team<br />

Harvard-Smithsonian Center for Astrophysics<br />

1<br />

1


Overview<br />

• <strong>The</strong> Fe<strong>at</strong>ure Finding Team (FFT) is an intern<strong>at</strong>ional<br />

consortium with the task to produce a comprehensive<br />

autom<strong>at</strong>ed fe<strong>at</strong>ure and event recognition system for <strong>SDO</strong>.<br />

• It provides several different modules to look <strong>at</strong> different<br />

fe<strong>at</strong>ures and events.<br />

• <strong>The</strong> flare detective is the module with the<br />

detecting flares from AIA images<br />

task of<br />

Paolo Grigis<br />

Harvard-Smithsonian Center for Astrophysics<br />

2<br />

2


Inner working<br />

• FFT modules works a bit differently then the way one<br />

would build such a module if all d<strong>at</strong>a was available<br />

• <strong>The</strong> flare detective is fed one image <strong>at</strong> a time<br />

• <strong>The</strong>re is no options to retrieve th<strong>at</strong> image <strong>at</strong> a l<strong>at</strong>er time<br />

• All relevant inform<strong>at</strong>ion about past images is stored in an<br />

internal st<strong>at</strong>us of the module<br />

• One instance of the flare detective works with images<br />

from one AIA channel (independently from other<br />

instances of the module th<strong>at</strong> may be running on otehr<br />

channels)<br />

Paolo Grigis<br />

Harvard-Smithsonian Center for Astrophysics<br />

3<br />

3


New AIA<br />

image<br />

4<br />

4


New AIA<br />

image<br />

<strong>Flare</strong> <strong>Detective</strong><br />

processes image<br />

5<br />

5


New AIA<br />

image<br />

<strong>Flare</strong> <strong>Detective</strong><br />

processes image<br />

<strong>Flare</strong>s detected<br />

flares sent<br />

to HEK<br />

6<br />

6


New AIA<br />

image<br />

<strong>Flare</strong> <strong>Detective</strong><br />

processes image<br />

<strong>Flare</strong>s detected<br />

flares sent<br />

to HEK<br />

<strong>Flare</strong> <strong>Detective</strong><br />

upd<strong>at</strong>es its<br />

internal st<strong>at</strong>us<br />

7<br />

7


New AIA<br />

image<br />

<strong>Flare</strong> <strong>Detective</strong><br />

processes image<br />

<strong>Flare</strong>s detected<br />

flares sent<br />

to HEK<br />

<strong>Flare</strong> <strong>Detective</strong><br />

upd<strong>at</strong>es its<br />

internal st<strong>at</strong>us<br />

New AIA<br />

image<br />

<strong>Flare</strong> <strong>Detective</strong><br />

processes image<br />

8<br />

8


New AIA<br />

image<br />

<strong>Flare</strong> <strong>Detective</strong><br />

processes image<br />

<strong>Flare</strong>s detected<br />

flares sent<br />

to HEK<br />

<strong>Flare</strong> <strong>Detective</strong><br />

upd<strong>at</strong>es its<br />

internal st<strong>at</strong>us<br />

New AIA<br />

image<br />

<strong>Flare</strong> <strong>Detective</strong><br />

processes image<br />

<strong>Flare</strong>s detected<br />

flares sent<br />

to HEK<br />

9<br />

9


Inner working (2)<br />

• <strong>The</strong> large AIA images are segmented in a small number<br />

of macropixels<br />

• Lightcurves are computed for each macropixels<br />

• <strong>Flare</strong>s start times are determined when the deriv<strong>at</strong>ive of<br />

the lightcurve reaches a certain threshold<br />

• <strong>Flare</strong>s end times are determined when the flux is lower<br />

then a certain fraction of the peak flux<br />

Paolo Grigis<br />

Harvard-Smithsonian Center for Astrophysics<br />

10<br />

10


193A<br />

Paolo Grigis — Harvard-Smithsonian Center for Astrophysics — Hinode Meeting Palermo 2010<br />

11


Paolo Grigis — Harvard-Smithsonian Center for Astrophysics — Hinode Meeting Palermo 2010<br />

12


Rebinned AIA image with lightcurves<br />

Paolo Grigis<br />

Harvard-Smithsonian Center for Astrophysics 13<br />

13


14<br />

14


Lightcurve<br />

Deriv<strong>at</strong>ive<br />

15<br />

15


Lightcurve<br />

Deriv<strong>at</strong>ive<br />

Deriv<strong>at</strong>ive Threshold<br />

16<br />

16


Lightcurve<br />

<strong>Flare</strong> Start Detected<br />

Deriv<strong>at</strong>ive<br />

Deriv<strong>at</strong>ive Threshold<br />

17<br />

17


Lightcurve<br />

<strong>Flare</strong> Start Detected<br />

Deriv<strong>at</strong>ive<br />

Deriv<strong>at</strong>ive Threshold<br />

18<br />

18


Peak Strength<br />

Lightcurve<br />

<strong>Flare</strong> Start Detected<br />

Deriv<strong>at</strong>ive<br />

Deriv<strong>at</strong>ive Threshold<br />

19<br />

19


Peak Strength<br />

Lightcurve<br />

End fraction<br />

<strong>Flare</strong> Start Detected<br />

Deriv<strong>at</strong>ive<br />

Deriv<strong>at</strong>ive Threshold<br />

20<br />

20


Peak Strength<br />

Lightcurve<br />

End fraction<br />

<strong>Flare</strong> Start Detected<br />

<strong>Flare</strong> End Detected<br />

Deriv<strong>at</strong>ive<br />

Deriv<strong>at</strong>ive Threshold<br />

21<br />

21


challenges<br />

• <strong>Flare</strong>s in AIA look very different in different wavelength<br />

channels.<br />

• Lightcurves are simpler and more similar to the GOES<br />

lightcurves in the 94A and 131A channels (those are the<br />

hotter channels). <strong>The</strong> flare detection is most reliable in<br />

these channels.<br />

• Cooler channels show irregular spikes during the<br />

impulsive phase (presumably due to emission from the<br />

flare kernels) and secondary smoother peaks in the l<strong>at</strong>e<br />

phase (presumably due to cooling of hot flare m<strong>at</strong>erial in<br />

post-flare loops).<br />

•<br />

Paolo Grigis<br />

Harvard-Smithsonian Center for Astrophysics<br />

22<br />

22


23<br />

23


Performance<br />

• Code running on AIA images, stable<br />

• Events detected and exported to the Heliophysics<br />

Knowledgebase<br />

• Scientific valid<strong>at</strong>ion currently in progress<br />

• Thresholds are currently being optimized<br />

• Some channels deliver more events then others<br />

• 131A <strong>at</strong> the moment too sensitive<br />

Paolo Grigis<br />

Harvard-Smithsonian Center for Astrophysics<br />

24<br />

24


<strong>Flare</strong>s in 171A during the period 1-FEB-2011 1-MAY-2011<br />

25<br />

25


<strong>Flare</strong>s in 193A during the period 1-FEB-2011 1-MAY-2011<br />

26<br />

26


<strong>Flare</strong>s in 94A during the period 1-FEB-2011 1-MAY-2011<br />

27<br />

27


171A flares<br />

10.7 cm radio flux<br />

28<br />

28


D<strong>at</strong>a retrieval<br />

• Events detected can be retrieved by performing queries<br />

to the HEK d<strong>at</strong>abase<br />

• Example IDl code:<br />

IDL> query=ssw_her_make_query(StartTime,EndTime,/fl, $<br />

search ='frm_name=<strong>Flare</strong> <strong>Detective</strong> - Trigger Module')<br />

IDL> events=ssw_her_query(query,/struct4event)<br />

• <strong>Flare</strong> list in ASCII form<strong>at</strong> will also be delivered<br />

Paolo Grigis<br />

Harvard-Smithsonian Center for Astrophysics<br />

29<br />

29


Conclusion<br />

• <strong>The</strong> flare detective is running on AIA d<strong>at</strong>a and delivering<br />

preliminary events, however some degree of caution<br />

should be used <strong>at</strong> this point<br />

• Scientific valid<strong>at</strong>ion currently in progress<br />

• Look forward for the AIA flare list coming out soon!<br />

• Future development: more detailed characteriz<strong>at</strong>ion of the<br />

events detected<br />

Paolo Grigis<br />

Harvard-Smithsonian Center for Astrophysics<br />

30<br />

30

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!