1 / 37

Observational Test of Coronal Magnetic Field Models I. Comparison with Potential Field Model

Observational Test of Coronal Magnetic Field Models I. Comparison with Potential Field Model. Hao-Sheng Lin & Yu Liu Institute for Astronomy University of Hawaii. ‘Vector’ Coronal Magnetogram of AR 10581 and AR 10582, 2004. Transverse field orientation. Longitudinal Field Strength.

marvel
Télécharger la présentation

Observational Test of Coronal Magnetic Field Models I. Comparison with Potential Field Model

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. Observational Test of Coronal Magnetic Field Models I. Comparison with Potential Field Model Hao-Sheng Lin & Yu Liu Institute for Astronomy University of Hawaii SHINE 2008, June 23-27, Utah

  2. ‘Vector’ Coronal Magnetogramof AR 10581 and AR 10582, 2004 Transverse field orientation Longitudinal Field Strength Contour plot of the line-of-sight magnetogram over-plotted on the EIT Fe XVI 284 A image. The contours are 5G, 3G, and 1G. SHINE 2008, June 23-27, Utah

  3. Current Status of IR Coronal Magnetometry Instrumentation • We can (almost) make routine coronal magnetic field measurements… • SOLARC is working… • CoMP is been relocated to Haleakala… Capabilities • Linear polarization can be obtained easily…even during solar minimum. • Neither has the sensitivity for longitudinal magnetic field measurement during solar minimum. SHINE 2008, June 23-27, Utah

  4. SOLARC: Off-Axis Mirror Coronagraph Secondary mirror Prime focus inverse occulter/field stop Re-imaging lens LCVR Polarimeter Input array of fiber optics bundle Echelle Grating Camera Lens Collimator NICMOS3 IR camera Fiber Bundle Primary mirror Optical Configuration of SOLARC and OFIS SOLARC on the summit of Haleakala, Maui. 50 cm aperture SHINE 2008, June 23-27, Utah

  5. Coronal Multi-channel Polarimeter (CoMP) – S. Tomczyk, HAO Calibration Polarizer Stage Light Trap Occulting Disk Collimating Lens Pre-Filter Wheel Filter/Polarimeter Re-Imaging Lens Detector 20-cm coronagraph, National Solar Observatory, Sac Peak SHINE 2008, June 23-27, Utah

  6. March 9, 2004 SHINE 2008, June 23-27, Utah

  7. Coronal Magnetometry 101 Polarization Mechanism • The forbidden coronal emission lines in the visible and IR are polarized by the Saturated Hanle Effect. Information Content • Linear Polarization (easy to measure): • The direction of the linear polarization yields the direction of the magnetic field projected in the plane of the sky containing sun center. • Linear polarization does not yield information about the strength of the magnetic fields. • The magnetic field direction intepretation is subjected to a 90 degree ambiguity--- the van Vleck effect. • Circular Polarization (difficult to measure): • Similar to the photospheric Zeeman effect, yields line-of-sight magnetic field strength, • with an alignment effect correction. SHINE 2008, June 23-27, Utah

  8. What Can we Do with these Coronal Magnetic Field Measurements? 1. Build a Coronal Magnetic Field Model. SHINE 2008, June 23-27, Utah

  9. Can we invert the polarization measurements to derive the 3d magnetic field structure of the corona? SHINE 2008, June 23-27, Utah

  10. The Full Inversion Problem • The coronal atmosphere is optically thin, and the observed coronal polarization signals may not originate from a single localized source along the line of sight. • There are many independent parameters in the model, but only a few observables… The inversion problem is severely under constrained! • There are currently no tested, reliable inversion methods for reconstruction of the 3D magnetic field structure of corona using polarization measurements… If inversion is not possible, what can we do? SHINE 2008, June 23-27, Utah

  11. What Can we Do with these Coronal Magnetic Field Measurements? • Build a Coronal Magnetic Field Model • Check extrapolation/simulation Models Given a coronal magnetic field model derived from extrapolation or MHD simulation of observed photospheric magnetic fields, can we verify the validity of this model? SHINE 2008, June 23-27, Utah

  12. Forward Modeling… Yes! In principle… If we know the 3-dimensional • magnetic field, • density, and • temperature structure of the corona, then we can calculate what the emergent linear and circular polarization signals should look like, and compare them with the observed polarization signals… SHINE 2008, June 23-27, Utah

  13. In reality… • Extrapolations only yield the magnetic field configuration. • There are no information about n and T. • n and T has to be derived, inferred, assumed, or guessed by other means… • The photospheric and coronal magnetic field observations are not co-temporal… • Uncertainties due to evolution of the active region. • Potential and force-free assumption may not be valid. SHINE 2008, June 23-27, Utah

  14. Limitations of Forward Modeling… • If the observed and synthesized polarization signals match…The model may be good, • But we don’t know if this is the best model… • If the observed and synthesized polarization signals don’t match… • The magnetic field model may be good, but • The density and temperature model is not good, • The corona may have changed… • Can be anything… SHINE 2008, June 23-27, Utah

  15. Testing Potential Field Extrapolation of AR 10582 Can we reproduces these Fe XIII 1075 nm polarization measurements from a coronal magnetic field model derived from potential field extrapolation of observed photospheric magnetogram? SHINE 2008, June 23-27, Utah

  16. About AR 10581 AND 10582… • Flaring activities in AR10582 ceased about 5 days before our coronal Bobservation… • Potential field extrapolation may be OK? Time of coronal magnetic field observation… SHINE 2008, June 23-27, Utah

  17. Potential Field Model of AR 10581 and 10582 SHINE 2008, June 23-27, Utah

  18. Density and Temperature… Make some educated guesses… • The coronal is bright in the active regions, • The density falls off exponentially. • Assume a uniform temperature through out the corona… • For density: assume four empirical models: • Uniform density, ne = constant. • Gravitationally stratified, ne ~ e-h/0.15, Sn ~ ne2. • Snweighted by B: S ~ e-2h/0.15 B. • Sn weighted by B2: S ~ e-2h/0.15 B2. SHINE 2008, June 23-27, Utah

  19. Line-of-Sight direction SHINE 2008, June 23-27, Utah

  20. Photospheric Magnetic flux LOS direction (Lower Panels) SHINE 2008, June 23-27, Utah

  21. None of the empirical models have produced synthesized linear polarization maps that agree with the observed one. However… • The source function weighted by B2, with the most narrow width gave the best result. • Should we try source function with even narrower width (along the LOS direction)? SHINE 2008, June 23-27, Utah

  22. TRACE images have demonstrated long ago that the coronal intensity structures have characteristics scale much smaller than the 200 ~ 300 km width we used… • Try using the loop width of ~ 56 km (7 times Fe XVI 284 A characteristic loop width, Aschwanden et al., 2000) as the width of the source function. SHINE 2008, June 23-27, Utah

  23. Since the thickness of the new source function is small, we computed the synthesized LP map as a function of position along the Line of Sight… SHINE 2008, June 23-27, Utah

  24. SHINE 2008, June 23-27, Utah

  25. p: rms error in degree of linear polarization, • : rms error of azimuth angle of LP, • LP: combined rms error of the degree of polarization and the azimuth angle of the LP. Judge’s code Lin’s synthesis code Source function of the best-fit layer LOS direction SHINE 2008, June 23-27, Utah

  26. Comparison with Judge’s Synthesis Code • Phil Judge’s code includes collisional depolarization effect… • We really can’t tell which code is better from this comparison. But Judge’s code includes more physics… SHINE 2008, June 23-27, Utah

  27. Circular Polarization • The sensitivity of this dataset is not sufficient for a point-to-point comparison within the FOV. • We only compared the averaged LOS magnetic field strength as a function of height. SHINE 2008, June 23-27, Utah

  28. The longitudinal B as a function of height h, as well as the height of reversal of B(h) are calculated for each layer. • The height of B reversal agrees with the observed value at two layers. • B(h) at layer 130 fits the observed one better. SHINE 2008, June 23-27, Utah

  29. SHINE 2008, June 23-27, Utah

  30. Is this a coincidences? • Two independent parameters • the degree of linear polarization, and • azimuth angle of linear polarization have minimum at about the same location… • Two other independent parameter/function • the height of Stokes V reversal, and • B(h) match the observed circular polarization signals at the same location… SHINE 2008, June 23-27, Utah

  31. Conclusions • Potential field extrapolation of AR 10582 has reproduced the observed coronal linear and circular polarization maps. • The inferred source functions of the linear and circular polarization are fairly localized! • Single-source inversion (Judge 2007) might be possible… • The LP and CP source functions are close to the location of the sunspot of the active region. • The locations of the source function of the linear and circular polarization are not the same… • This potential field model may be OK! Liu & Lin, 2008, ApJ, 680, 1496 SHINE 2008, June 23-27, Utah

  32. What’s Next? More observations (if the Sun cooperates) and more comparisons with models… • Is potential-field extrapolation really OK? • Does force-free extrapolations provide better model? • MHD models should come with information about n and T… • Direct comparison can be performed without guessing where the source is located. SHINE 2008, June 23-27, Utah

  33. Please Help! • A future SHINE Session? Observational Test of Coronal Magnetic Field Models Construct the coronal magnetic field model of AR 10582 using your favorite extrapolation/MHD code, and we can start testing these models vigorously… SHINE 2008, June 23-27, Utah

  34. SHINE 2008, June 23-27, Utah

  35. SHINE 2008, June 23-27, Utah

  36. SHINE 2008, June 23-27, Utah

  37. Testing Coronal B Model by Forward Modeling • Construct a coronal magnetic field model, • Potential field extrapolation, • Force-free field extrapolation, • MHD simulation, from photospheric magnetic field observation several days before or after the coronal observation… • Construct a thermodynamic model, • Density n • Temperature T we don’t really have reliable measurement of n and T! • Calculate the emergent Stokes vector, • Compare the observed and synthesized polarization signals. SHINE 2008, June 23-27, Utah

More Related