eSDO 1121: Magnetic Field Extrapolation

This document can be viewed as a PDF.
Deliverable eSDO-1121: Magnetic Field Extrapolation
M. Smith, E.Auden, L. van Driel-Gesztelyi
10 August 2005


The Magnetic field extrapolation process will calculate magnetic field lines between Solar active regions using vector and line-of-sight magnetograms from the SDO HMI instrument. The extrapolated fields will be shown as an overlay on the magnetogram and their footpoints displayed in tabular form.

The primary investigation will concentrate on Wiegelmann's [1] improved 'Optimization' method of extrapolation, which uses the nonlinear magnetic model. The nonlinear model provides the most accurate representation of magnetic fields, particularly around active regions. The Optimization technique, originally developed by Wheatland [2], computes the magnetic field over a predefined volume, or box, using full-disk vector magnetogram data for the extrapolation. The magnetogram is also used to determine the bottom boundary conditions of the box, while the lateral and top boundaries are computed from a Potential Field model, which provides a valid approximation of magnetic field activity high in the the Sun's corona.

Wiegelmann's Optimization method improves on other nonlinear methods by introducing the concepts of a boundary layer and weighting function. These diminish the influence of the lateral and top boundary conditions on the calculation of the magnetic field over the region of interest.

The vector magnetogram data will be preprocessed in two stages before it is submitted to the Optimization method for extrapolation. These are (in order):

  • Azimuth Disambiguation
  • Wiegelmann and Sakurai's preprocessing procedure

Azimuth Disambiguation

'Azimuth Disambiguation' (Georgoulis [3]) is a new technique designed to resolve the Pi ambiguity, one of several problems inherent in modern vector magnetograms. The Pi ambiguity describes the case where the azimuth angle of the transverse component of the magnetic field has two equally likely values 180 degrees apart. Azimuth Disambiguation will typically add a 1-2 minute processing time penalty to the overall extrapolation computation on an average desktop computer.

Wiegelmann, Inhester and Sakurai's preprocessing procedure

Another drawback of currently available vector magnetograms is the level of noise that is sometimes apparent in the transverse components. Such noisy data can violate the force-free assumptions used in the Optimization extrapolation and invalidate the results. Wiegelmann, Inhester and Sakurai [4] have recently developed a numerical technique that detects any data lying outside force-free constraints and 'smooths' them into a force-free condition. The amount of additional processing time incurred by this technique is unknown at present.

Should the Optimization method prove computationally expensive, then the 6 nonlinear force-free algorithms discussed at the Lockheed Martin NLFFF modelling meeting in May 2005 will be reviewed; the algorithm which provides the best balance between speed and accuracy will be deployed.

Again, if none of these is acceptably fast, then the simpler (but less accurate) Potential Field model will be used.


  • HMI active region vector magnetogram.
  • HMI line-of-sight magnetogram.


  • FITS file containing extrapolated fields overlayed on the original magnetogram.
  • FITS file tabular extension containing footprints for each extrapolated field. A footprint will be described in terms of solar coordinates and radius.

Test Data

  • SOHO MDI vector magnetograms, migrating to Solar-B vector magnetograms

Tool Interface

  • commandline:
    1. AstroGrid CEA web service: this algorithm will be deployed as a CEA service hosted in the UK that users can call to process datasets on the grid. Note: due to the computational intensity of this algorithm, access to the CEA service may be restricted to recognized solar users registered with AstroGrid.
    2. SolarSoft routine: if computational intensity is not deemed too high, the C module will be wrapped in IDL and distributed through the MSSL SolarSoft gateway. Users will access to a SolarSoft installation can call the routine from the commandline to process locally held data.
    3. JSOC module: the C module will be installed in the JSOC pipeline. Users can access the routine through pipeline execution to operate on data local to the JSOC data centre.

Science Use Case

The user wants to view extrapolated magnetic fields superimposed on a vector magnetogram image, and also a description of each field's footpoint that includes solar coordinates and radius.

  1. The user identifies the image taken during the specified time period.
  2. The user input the full-disk magnetogram HMI image to the Magnetic Field Extrapolation algorithm. (Currently there are no user-specified variables to describe for the loop recognition algorithm to be distributed through SolarSoft or deployed as an AstroGrid CEA service.)
  3. The algorithm is run and returns a FITS file to the user.
  4. The user can view an image within the FITS file displaying extrapolated magnetic fields superimposed in colour over the original HMI magnetogram image.
  5. The user can also view a table of footprints for each of the fields identified. The footprints are described in terms of solar coordinates and radius.

Technical Use Case


  1. The Magnetic Field Extrapolation algorithm receives a full-disk vector magnetogram as input.
  2. Georgoulis's Azimuth Disambiguation technique is applied and used to remove the Pi Ambiguity of the transverse components of the vector magnetogram data.
  3. Wiegelmann, Inhester and Sakurai's preprocessing procedure is applied to check for and 'smooth' magnetogram data that lie outside force-free constraints.
  4. A computational box is defined, i.e. a bounded volume which represents the region in which the magnetic field will be calculated.
  5. A physical domain is defined inside the box which represents the volume in which the nonlinear magnetic field will be calculated using vector magnetogram data.
  6. A boundary layer is defined which stretches from the edge of the physical domain to the computational box boundary.
  7. The measured normal component of the magnetic field, Bz, is used to calculate a potential magnetic field over the whole box using Seehafer's method [5] and is used to define the lateral and top boundary conditions.
  8. The vector magnetogram data is used to determine the bottom boundary (photosphere) conditions of the box.
  9. The nonlinear magnetic field is computed from vector magnetogram data using the Landweber iteration method [6]. A weighting cosine function is applied which is set to unity within the physical domain, decreases within the boundary layer and falls to zero at the lateral and top boundaries of the computational box.

Quicklook Products


Support Information

  1. T.Wiegelmann, 'Optimization code with weighting function for the reconstruction of coronal magnetic fields', Solar Physics, 219, 87-108 (2004)
  2. M.S.Wheatland, P.A. Sturrock and G.Roumeliotis, AstroPhysical Journal 540, 1150-1155.
  3. M. Georgoulis, A New Technique for a Routine Azimuth Disambiguation of Solar Vector Magnetograms, AstroPhysical Journal 629:L69-72.
  4. T. Wiegelmann, B.Inhester, T.Sakurai, 'Preprocessing of vector magnetograph data for a nonlinear force-free magnetic field reconstruction', Solar Physics, In Press (2005).
  5. N.Seehafer, Solar Physics, 58, 215 (1978).
  6. A.K.Louis, 'Inverse und schlecht gestellte Probleme', Teubner Studienbuecher, ISBN 3-519-02085-X (1989) (discusses Landweber's iteration method).

-- MikeSmith - 19 Aug 2005

This topic: SDO > PhaseAMagneticFieldExtrapolation
Topic revision: r8 - 2005-09-30 - ElizabethAuden
This site is powered by the TWiki collaboration platform Powered by PerlCopyright © 2008-2020 by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding TWiki? Send feedback