1 / 28

Surface Flows From Magnetograms

Surface Flows From Magnetograms. Brian Welsch, George Fisher, Bill Abbett, & Yan Li Space Sciences Laboratory, UC-Berkeley M.K. Georgoulis Applied Physics Laboratory, Johns Hopkins University P. W. Schuck Naval Research Lab – D. C. D.W. Longcope & B. Ravindra

ceri
Download Presentation

Surface Flows From Magnetograms

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Surface Flows From Magnetograms Brian Welsch, George Fisher, Bill Abbett, & Yan Li Space Sciences Laboratory, UC-Berkeley M.K. Georgoulis Applied Physics Laboratory, Johns Hopkins University P. W. Schuck Naval Research Lab – D. C. D.W. Longcope & B. Ravindra Physics Department, Montana State University Marc DeRosa Lockheed-Martin Advanced Technology Center K. Kusano Earth Simulator Ctr., Japan Agency for Marine-Earth Sci. & Tech.

  2. Main Ideas 1. Why should we study surface flows derived from magnetograms? Poynting/Helicity Fluxes, Data Driving, Flux Transport Models 2. What classes of flows can we derive, in principle, from magnetogram sequences? In theory, flows perpendicular to both B &Bz can be found. 3. Practically, what are the best available methods for deriving flow information? Tests suggest tracking (e.g., LCT) combined with MEF. 4. How should available methods be improved? We must develop the capability to match both ∂ t B hor &∂ t Bz.

  3. Abstract Estimates of velocities from time series of photospheric and/or chromospheric vector magnetograms can be used to determine transport rates for magnetic flux (i.e., emergence and submergence), energy (Poynting flux), and helicity across the magnetogram layer, and to provide time-dependent boundary conditions for data-driven simulations of the solar atmosphere above this layer. Velocity components perpendicular to the magnetic field are necessary both to compute these transport rates and to derive model boundary conditions. Since Doppler shifts also contain contributions from flows parallel to the magnetic field, perpendicular velocities are not generally recoverable from Doppler shifts alone. Consequently, several methods have been developed to estimate the perpendicular velocity from magnetograms, and have recently been validated using synthetic magnetograms from MHD simulations. Combined with data from the next generation of magnetographs (SOLIS, SOT/Hindode, and HMI/SDO), these techniques should provide valuable tools for understanding and modeling solar variability on time scales both short (e.g., flares and CMEs) and long (e.g., the 11-year sunspot cycle).

  4. Flows, v, affect magnetic evolution at both the surface and in the corona. • Since E = -(v x B)/c + R, the fluxes of magnetic energy & helicityacross the surface depend upon v. ∂tU = c ∫ dA (E x B) ∙n / 4π (1) ∂tH = c ∫ dA (E x A) ∙n / 4π (2) • B in the corona is coupled to B at the surface, so the surface v provides an essential boundary condition for data-driven MHD simulations of the coronal B field. • Studying v could improve evolutionary models of photospheric B fields, e.g., flux transport models.

  5. Progress Slide #1

  6. The Induction Equation relates ∂tB to v,∂tB = - c ( x E) =  x [(v x B) + cR] (3) • Resistive effects, R, are model-dependent. • Using, e.g., R =η2B would require magnetograms at two heights in the atmosphere – these are rare! • R=0 (ideal MHD) is usually assumed. • Even ideally, only the normal component of (3) does not require magnetograms from two heights.

  7. Consequently, the goal is to derive v from∂tBn= ∙(vnBhor- vhorBn) (4) • Clearly, (4) can’t determine v, which has three components! • By assumtion, v||, the component of v parallel to B, cannot affect ∂tBn (v x B). In fact, v|| is not useful for any of the reasons I’ve given for wanting to study surface flows! • So ignore v||! With v= v + v||, try to find v , and use v∙ B = 0 (5) as an additional equation to solve for v. • But the system is still not closed!

  8. vLOS vLOS vLOS Aside: Doppler shifts cannot close the system! Generally, Doppler shifts cannot distinguish flows || to B (red), perp. to B (blue), or in an intermediate direction (gray). With v estimated another way & projected onto the LOS, the Doppler shift determines v|| (Georgoulis & LaBonte, 2006) Doppler shifts are only unambiguous along polarity inversion lines, where Bn changes sign (Chae et al. 2004, Lites 2005).

  9. Dopplergrams are sometimes consistent with “siphon flows” moving along the magnetic field. Left: MDI Dopplergram at 19:12 UT on 2003 October 29 superposed with the magnetic neutral line. Right: Evolution of the vertical shear flow speed calculated in the box region of the left panel. The two vertical dashed lines mark the beginning and end of the X10 flare. (From Deng et al. 2006)

  10. Aside #2: In addition to v||, other “null flows,” for which ∂tBn = 0, lie in the null space of (4). • Any flow obeying (vnBhor- vhorBn) =  x fn, for some scalar function f, lies in the null space of (4). • For instance, “contour flows” -- flows along contours of Bn -- do not alter Bn, but correspond to twisting motions. Hence, they inject energy and helicity into the corona! • No method based upon ∂tBn alone will estimate null flows well. • Happily, magnetograms usually have sufficient spatial structure that the null space of (4) is expected to be small.

  11. Progress Slide #2

  12. Several approaches have been used with ∂tBn and v∙ B = 0 to solve for v. • Kusano et al. (2002), Welsch et al. (2004), and Schuck (2006) used tracking methods to close the system. • Longcope (2004) minimized ∫ dA (vn2 + vhor2). • Georgoulis & LaBonte (2006) assumed vn = 0. Tracking methods, which follow the apparent motions of magnetogram features, have the longest heritage.

  13. The apparent motion of magnetic flux in magnetograms is the flux transport velocity, uf. Démoulin & Berger (2003): In addition to horizontal flows, vertical velocities can lead to uf0. In this figure, v= 0, but vn0, souf0. uf is not equivalent to v; rather, ufvhor - (vn/Bn)Bhor • uf is the apparent velocity (2 components) • v is the actual plasma velocity (3 comps) (NB: non-ideal effects can also cause flux transport!)

  14. Aside #3: Conformal mappings should be used to project magnetograms for tracking! • Apparent motions in the plane-of-sky do not accurately represent surface flows. • If spherical data are to be projected onto a plane prior to tracking, conformal (shape-preserving) mappings should be used. • Equal-area (authalic) projections distort shapes, and will therefore distort the directions of estimated flows. • We currently use Mercator projections, which distort length scales, but in a way that is easy to correct.

  15. Does the flux transport velocity agree with white light velocities? Yes & no! Tracking MDI full-disk white light images & mag- netograms gives different results. Data are from AR 8210 on 01 May ’98; Δt = 96 min.

  16. Recently, we conducted quantitative tests & comparisons of several available methods. - We created “synthetic magnetograms” from ANMHD simulations of an emerging flux rope.- In these data, both v & B are known exactly.

  17. Via several methods, we estimated v from N= 7 pairs of magnetograms, with increasing Δt’s. • We verified that the ANMHD data were consistent with ∂tBn= ∙(vnBhor- vhorBn). • Here, I show representative results from just a few of the methods tested: • Fourier LCT (FLCT, Welsch et al. 2004) • Inductive LCT (ILCT, Welsch et al. 2004) • Minimum Energy Fit (MEF, Longcope 2004) • Differential Affine Velocity Estimator (DAVE, Schuck 2006)

  18. We verified that the estimated v’s also obeyed ∂tBn = ∙ (vnBhor– vhorrBn). ILCT FLCT DAVE MEF

  19. We tested Démoulin & Berger’s relation of uf to v. FLCT ILCT Estimated uf’s are highly correlated with ANMHD’s uf. MEF DAVE

  20. Not surprisingly, the methods’ performance worsened as the time between magnetograms increased. % vector errors (direction & magnitude) were at least 50% (!!!). % speed errors (magnitude) were smaller, but biases were seen.

  21. Some methods estimated the direction of v to within ~30º, on average. CVEC and CCS were as defined by Schrijver et al. (2005):

  22. How well do we do at estimating the E in the fluxes of magnetic energy and helicity? FLCT ILCT MEF DAVE

  23. While most methods’ E is highly correlated with ANMHD’s E, only MEF gets the energy flux right. FLCT ILCT MEF DAVE Unfortunately, even MEF exhibits a large amount of scatter.

  24. As with the energy flux, only MEF accurately estimates the helicity flux.

  25. Conclusions from Tests • Tests using MHD magnetograms, in which B is error-free, show that estimated v’s are correlated with true v’s, but are inaccurate. • Only MEF reproduced energy and helicity fluxes well; other methods erred by > 50%. • Further tests (Ravindra & Longcope, 2007, in prep.) show that combining MEF with tracking does enhances performance.

  26. What lies ahead? • All of the methods we tested were developed to find v’s that match tBn. • But we also have observations of tBhor! • The Vector Magnetogram Fitting (VMF) technique developed by McClymont, Jiao, & Mikic (1997) was developed to drive changes ΔBhor. • A combined approach, using tracking and VMF, might allow incorporating tBhor into estimates of v.

  27. Progress Slide #3

More Related