characterizing photospheric flows david hathaway (nasa/msfc) with john beck & rick bogart...
TRANSCRIPT
Characterizing Photospheric Flows
David Hathaway (NASA/MSFC)
with
John Beck & Rick Bogart (Stanford/CSSA)
Kurt Bachmann, Gaurav Khatri, & Joshua Petitto (Birmingham Southern College)
Sam Han & Joseph Raymond (Tennessee Technological University)
Åke Nordlund (Astronomical Observatory, Denmark)
and of course
Phil Scherrer and the MDI Team
OUTLINE
•Why bother?•Simulating Photospheric Flow Data•MDI Data•Data Preparation: p-mode filtering•Axisymmetric Flow Analysis•The Convection Spectrum (Poloidal Component) •Radial Flow Component•Toroidal Flow Component•Angular Momentum Flux by Cellular Flows•Rotation Rates•Lifetimes•Conclusion
THREE BIG QUESTIONS IN SOLAR PHYSICS
1) What causes solar flares, prominence eruptions, and CMEs?
2) What causes the 11-year cycle ofSunspots and solar activity?
3) How is the corona heated and thesolar wind accelerated?
Data Simulation•Calculate vector velocities using input spectrum of complex spectral coefficients, R, S, and T (Chandrasekhar, 1961).
•Project vector velocities into the line-of-sight and integrate over pixels to get Doppler velocity signal
m
mm
mm
m
mm
mm
m
m
mr
YT
YSV
YT
YSV
YRV
,
,
,
sin
1,
sin
1,
,,
sincos,
coscoscossinsin,
cossincoscossin,,
0
00
00
BV
BBV
BBVV rlos
Data Simulation - ContinuedThe resulting Doppler velocity signal is convolved with a point-spread-function representing the MDI optical system in Full Disk mode to simulate MDI Doppler velocity data.
The final Doppler image is analyzed using the same processes used on the MDI data and results are compared to determine the complex spectral coefficients and their dependence upon ℓ and time.
MDI Doppler Velocity Data
p-mode Filter Transmission
This gets mapped onto heliographic coordinates and projected onto spherical harmonics to give spherical harmonic spectral coefficients.
2D Spectrum from Single Image
Power Spectrum from the 61-day SequenceThe peak at ℓ~120 represents supergranules with
typical size λ ~ 2πR/ℓ ~ 36 Mm
Velocity Spectrum from the 61-day SequenceThe peak at ℓ~140 represents supergranules with
typical size λ ~ 2πR/ℓ ~ 31 Mm andtypical Doppler velocity ~ 100 m/s
Data Simulation #1
This simulation only includes horizontal, poloidal flows with:
mim
m
m
eSS
T
R
)(
0
0
and random phases, m.
Close-up of Region near Disk Center
Observed Power Spectrum from Real and Simulated Data
Observed Velocity Spectrum from Real and Simulated Data
Input Power Spectrum
Input Velocity Spectrum
Low Resolution
High Resolution
1D Power Spectrum from Hi-Resolution Data
Convection Spectrum from Giant Cells to Granules
Photospheric Convection Spectrum Conclusions
•There is a continuous spectrum of convective flows from giant cells to granules
•There are only two distinct modes of convection in the photosphere – granulation and supergranulation
Radial Flow Component Study
The line-of-sight velocity at a point, (x,y), on the disk isgiven by:
cos,sin,, 1 yxVyxVyxV rhlos
where ρ is the heliocentric angle from disk center with:
Ryx /sin2/122
and a second horizontal component, Vh2, is transverse to theline-of-sight.
22221
2
1
12222
12
sin
0
cossin2cossin
rrhlos
rh
rhrhlos
VVVV
VV
VVVVV
Radial Flow Component Study
If we consider the average mean-squared velocity at an angle ρ from disk center we find:
Radial Flow Component Study
222222 sinsin rrhlos VVVV
Exclude Magnetic Elements(|B| > 25 G)
Spectrum Filters for Isolating Cells of Different Sizes
Radial Flow Component Study Results
m
mm
m
mm
m
m
mr
YSV
YSV
YRV
,
,
,
sin
1,
,
,,
Radial Flow Component Study Results
We must include a radial component to the spectrum of the cellular flows with Rℓ=(0.05 + 0.07 ℓ/1000)Sℓ.
Radial Flow Study Conclusions
•Radial flow is ~10% of the horizontal flow for supergranules (increases to ~15% for cells 4000 km across)
•Doppler velocities under-estimate the horizontal flow speed by ~1.4
Toroidal Flow Study
TangentialRadial
Solenoidal flows have strong radial gradients and weak tangential gradients.
Toroidal flows have weak radial gradients and strong tangential gradients.
Foreshortening near the limb weakens the radial gradients more than the tangential gradients.
Solenoidal and toroidal flows will have different center-to-limb behavior.
Toroidal Flow Study (Data Simulation w/o MTF)
Solenoidal only – solid line, Toroidal only – dashed line
Toroidal Flow Study Comparison with MDI Data
MDI – solid line, Data Simulation (30% Toroidal) – dashed line
Toroidal Flow Study Conclusions
•Toroidal flow ~30% of the horizontal flow is consistent with data (but effects of MTF must be determined)
•Phase relationship between Solenoidal and Toroidal components needs to be determined
Angular Momentum Transport
•Axisymmetric flow (the meridional circulation) transports angular momentum toward the poles
•Non-axisymmetric flows (cellular flows) can transport angular momentum toward the equator if prograde velocities (u') are correlated with equatorward velocities (v')
onAccelerati Equatorial0 vu
A Doppler Velocity Indicator ofAngular Momentum Transport
Gilman (A&A 1977) showed that the mean squared Doppler velocity signal should be stronger east of the central meridian for an equatorward transport of angular momentum.
Coriolis force on flows in north-south cells gives equatorward transport and stronger Doppler signal east of the CM.
Coriolis force on flows in east-west cells gives poleward transport and stronger Doppler signal west of the CM..
Angular Momentum Flux Study
Angular Momentum Flux Study
Test for Instrumental EffectsUsing p-mode Signal
Cell Shape Analysis
Compare observed spectral amplitudes as a function of m/ℓ for MDI data and simulated data. Simulated data has no preferred cell shape.
Angular Momentum Study Conclusions
•Angular momentum is transported toward the poles in supergranules (instrumental problems – focus variations across disk and astigmatism – may severely impact this result)
Rotation Rate and Lifetime Study
Cross-Correlation TechniqueCross-correlating strips at the same latitude gives the shift (rotation rate) required to give maximum correlation (lifetime).
New Spectral Filters for Isolating Cells of Different Sizes
Rotation Rate vs. Feature Size
As the size of the features increases the rotation rate increases.
Data Simulation Rotation Rate vs. Feature Size
As the size of the features increases the rotation rate increases.
Rotation Rate vs. Time Difference
As the time lag between images increases the rotation rate increases.
Data Simulation Rotation Rate vs. Time Difference
As the time lag between images increases the rotation rate increases.
)1(
)(
)0()(;)0()(
sin
1,
,
,,
44
220
)()(
,
,
,
L
m
eStSeRtR
YSV
YSV
YRV
Lm
Lm
tmimmmtmimmm
m
mm
m
mm
m
m
mr
Rotation Rate Study Results
We include a systematic change in phase of the complex spectral coefficients.
Rotation Rate Study Conclusions
•Large cells appear to rotate more rapidly than small cells
•Larger lags between images give larger rotation rates
•Both effects are reproduced in the data simulation with a rotation rate that does not vary with time or wavenumber
Lifetime Study (in progress)
As the time lag between images increases the correlation coefficient decreases.
N
itmimmm
itmimmm
m
mm
m
mm
m
m
mr
eeStS
eeRtR
YSV
YSV
YRV
)(
)(
,
,
,
)0()(
)0()(
sin
1,
,
,,
Lifetime Study Approach
We include random changes in the phase of the complex spectral coefficients and add random noise to the image data
Lifetime Study Conclusions
•Large cells live longer than small cells (duh!)
•Lifetimes can be quantified using the data simulation
Conclusion
• Data simulations and forward modeling can be extremely valuable in interpreting observations