Dec 01, 2016 Download SOHO Notes for Mac to organize and share personal & workgroup info. What's new in version 9.3.7. Version 9.3.7 is special build of the software that only runs on macOS 10.12 (Sierra). free download. software for Mac::: SOHO Notes is the premier digital note-taking application for Macintosh Use it to capture organize and share all of your personal and workgroup information It can store rich text PDF's images files movies audio bookmarks and web archives. MacUpdate's software library contains more than 1615 Productivity apps designed for Mac owners. Choose the best app and download it today for free. Premier digital note-taking application. Show details Hide details. Download the latest version of SOHO Notes for Mac for free. Read 170 user reviews and compare with similar apps on MacUpdate. Download SOHO Notes for Mac to organize and share personal & workgroup info. What's new in version 9.3.7. Version 9.3.7 is special build of.
The success or failure of magnetic field modeling of the solar corona depends on both the choice of the theoretical model and on the choice of the used data sets. The simplest method is potential field modeling, which requires only a line-of-sight magnetogram Bz(x, y, zphot), but only few solar active regions match a potential field model. Also linear force-free field (LFFF) models are generally considered as unrealistic, where a single constant value for the force-free parameter α represents the multi-current system of an entire active region in the solar corona. The state-of-the-art is the nonlinear force-free field (NLFFF) model, which can accommodate for an arbitrary configuration of current systems, described by a spatially varying α(x) parameter distribution in an active region. The next strategic decision is the choice of data sets to constrain the theoretical model. NLFFF models generally require vector magnetograph data, B(x, y, zphot) = [Bx(x, y, zphot), By(x, y, zphot), Bz(x, y, zphot)], which are used as a lower boundary constraint and are extrapolated into the corona. However, a fundamental problem that has been identified is that the zone in the photosphere and lower chromosphere is not force free (Metcalf et al. 1995), which spoils the extrapolation into coronal heights and leads to a substantial mismatch between the extrapolated field lines and observed coronal loops, typically amounting to a 3D misalignment angle of μ ≈ 24°–44° (DeRosa et al. 2009; Sandman et al. 2009). Obviously, this fundamental problem can only be circumvented by using additional constraints from coronal data, since the corona above the transition region is generally force free, except when the plasma β parameter is larger than unity (e.g., in filaments) or when the equilibrium breaks down (e.g., during filament eruptions, flares, and coronal mass ejections).
The problem now amounts to how to implement coronal data into a theoretical magnetic field model, such as coronal loops, which are believed to be reliable tracers of the coronal magnetic field. A feasible approach is the method of stereoscopic triangulation, which can be applied by using the solar rotation or stereoscopic observations from dual spacecraft, as provided by the STEREO mission (for reviews of solar stereoscopy, see, e.g., Inhester 2006; Wiegelmann et al. 2009; Aschwanden 2011). Stereoscopically triangulated coronal loop coordinates [x(s), y(s), z(s)] (as a function of the curvilinear abscissa s) have been used to constrain: (1) potential field models in terms of buried unipolar magnetic charges (Aschwanden & Sandman 2010) or buried dipoles (Sandman and Aschwanden 2011), (2) linear force-free fields (Feng et al. 2007; Inhester et al. 2008; Conlon & Gallagher 2010), and (3) nonlinear force-free fields (Aschwanden et al. 2012b). The proof of concept to fit NLFFF codes to prescribed field lines was also demonstrated with artificial (nonsolar) loop data, fitting either 3D field line coordinates [x(s), y(s), z(s)], or 2D projections [x(s), y(s)] (Malanushenko et al. 2009, 2012). The method of Malanushenko et al. (2012) employs a Grad–Rubin-type NLFFF code (Grad & Rubin 1958) that fits an LFFF with a local α-value to each coronal loop and then iteratively relaxes to the closest NLFFF solution, while the method of Aschwanden et al. (2012b) uses an approximative analytical solution of a force-free and divergence-free field, parameterized by a number of buried magnetic charges that have a variable twist around their vertical axis, and is forward fitted to observed loop coordinates. The latter method is numerically quite efficient and achieves a factor of two better agreement in the misalignment angle (μ = 14°–19°) than standard NLFFF codes using magnetic vector data. However, the major limitation of the latter method is the availability of solar stereoscopic data, which restricts the method to the beginning of the STEREO mission (i.e., the year 2007), when STEREO had a small spacecraft separation angle that is suitable for stereoscopy (Aschwanden et al. 2012a).
Hence, the application of NLFFF magnetic modeling using coronal constraints could be considerably enhanced, if the methodical restriction to coronal 3D data, as it can be provided only by true stereoscopic measurements, could be relaxed to 2D data, which could be furnished by any high-resolution EUV imager, such as from the SOHO/EIT, Transition Region And Coronal Explorer (TRACE), and SDO/AIA missions. This generalization is exactly the purpose of the present study. We develop a modified code that requires only a line-of-sight magnetogram and a high-resolution EUV image, where we trace 2D loop coordinates to constrain the forward fitting of the analytical NLFFF code described in Aschwanden (2012a), and compare the results with those obtained from stereoscopically triangulated 3D loop coordinates (described in Aschwanden et al. 2012b). Furthermore, we also test magnetic forward fitting to automatically traced 2D loop data and compare the results with manually traced 2D loop data. The latter effort brings us closer to the ultimate goal of fully automated (NLFFF) magnetic field modeling with widely accessible input data.
The paper is organized as follows. The theory of the analytical NLFFF forward-fitting code is briefly summarized in Section 2, the numerical code is described in Section 3, tests of NLFFF forward fitting to simulated data are presented in Section 4, and to stereoscopic and single-image data in Section 5, while a discussion of the application is given in Section 6, with a summary of the conclusions provided in Section 7.
2. ANALYTICAL THEORY
A nonlinear force-free field (NLFFF) is implicitly defined by Maxwell's force-free and divergence-free conditions,
where α(x) is a scalar function that varies in space, but is constant along a given field line, and the current density j is co-aligned and proportional to the magnetic field B. A general solution of Equations (1) and (2) is not available, but numerical solutions are computed (see review by Wiegelmann & Sakurai 2012) using (1) force-free and divergence-free optimization algorithms (Wheatland et al. 2000; Wiegelmann 2004), (2) evolutionary magneto-frictional methods (Yang et al. 1986; Valori et al. 2007), or Grad–Rubin-style (Grad & Rubin 1958) current-field iteration methods (Amari et al. 1999, 2006; Wheatland 2006; Wheatland & Regnier 2009; Malanushenko et al. 2009). Numerical NLFFF solutions bear two major problems: (1) every method based on extrapolation of force-free magnetic field lines from photospheric boundary conditions suffers from the inconsistency of the photospheric boundary conditions with the force-free assumption (Metcalf et al. 1995) and (2) the calculation of a single NLFFF solution with conventional numerical codes is so computationally intensive that forward fitting to additional constraints (requiring many iteration steps) is unfeasible. Hence, an explicit analytical solution of Equations (1) and (2) would be extremely useful, which could be computed much faster and be forward fitted to coronal loops in force-free domains circumventing the non-force-free (photospheric) boundary condition.
An approximate analytical solution of Equations (1) and (2) was recently calculated (Aschwanden 2012a) that can be expressed by a superposition of an arbitrary number of Nm magnetic field components Bm, j = 1, ..., Nm,
where each magnetic field component Bm can be decomposed into a radial Br and an azimuthal field component Bφ,
where (r, φ, θ) are the spherical coordinates of a magnetic field component system (Bm, xm, ym, zm, αm) with a unipolar magnetic charge Bm that is buried at position (xm, ym, zm), has a depth d = 1 − (x2m + ym2 + z2m)1/2, a vertical twist αm, and r = [(x − xm)2 + (y − ym)2 + (z − zm)2]1/2 is the distance of an arbitrary coronal position (x, y, z) to the subphotospheric location (xm, ym, zm) of the buried magnetic charge. The force-free parameter α can also be expressed in terms of the parameter b (Equation (7)), which quantifies the number Ntwist of full twist turns over a (loop) length L,
This analytical approximation is divergence free and force free to second-order accuracy in the parameter (brsin θ), which is proportional to the force-free parameter α as defined by Equation (7). This approximate NLFFF solution is very appropriate for cases with small vertical twist, but may break down for highly nonpotential cases with large twist or magnetic field domains with strong horizontal twist, such as near-horizontal filaments or a Gold–Hoyle flux rope (Gold & Hoyle 1960; Aschwanden 2012a, Appendix A). In the limit of vanishing vertical twist (α↦0 or b↦0), the azimuthal component vanishes, Bφ↦0, and the radial component degenerates to the potential-field solution of a unipolar magnetic charge, Br↦Bm(d/r)2, which is simply a radial field that points away from the buried charge and decreases with the square of the distance. A numerical code that fits this analytical approximation to a given 3D magnetic field is described and tested in Aschwanden & Malanushenko (2012) using analytical models. Applications to real solar data using stereoscopically triangulated coronal loops from the STEREO mission, which supposedly outline force-free magnetic field lines in the solar corona, are presented in Aschwanden et al. (2012b).
3. NUMERICAL FORWARD FITTING
The numeric code that fits the approximative analytical NLFFF solution (Equations (3)–(7)) to a line-of-sight magnetogram Bz(x, y, zphot) plus coronal loop coordinates [x(s), y(s), z(s)] for an ensemble of stereoscopically triangulated loops in a solar active region is described in detail and tested in Aschwanden & Malanushenko (2012). We are using a Cartesian coordinate system (x, y, z) with the origin in the center of the Sun and the plane of the sky is in the (x, y)-plane, while the z-axis is the line of sight. This allows us to take the curvature of the solar surface into full account, in contrast to some other NLFFF codes that approximate the solar surface with a flat plane.
The forward-fitting part of the code consists of two major parts: (1) the decomposition of buried unipolar magnetic charges from a line-of-sight magnetogram Bz(x, y, zphot) (see Appendix A of Aschwanden et al. 2012b) and (2) iterative optimization of the nonlinear force-free α(x) parameters by minimizing the misalignment angles between the loop data and the fitted NLFFF model. The new approach in this work is the generalization of the forward-fitting code from 3D loop coordinates (using stereoscopic measurements before) to 2D loop coordinates [x(s), y(s)] only, which can simply be provided by any high-resolution EUV image without requiring stereoscopic views.
The two different methods are juxtaposed in Figure 1. In the 3D forward-fitting method, the 3D misalignment angle μ3(i, j), i = 1, ..., ns, j = 1, ..., nL is computed for a number of ns loop segment positions and nL coronal loops, defined by the scalar product that calculates the angle μ3 between the 3D vectors of the loop direction Bobs(x) and the magnetic field direction Btheo(x) at a given loop position x(s) = [x(sij), y(sij), z(sij)],
This is illustrated in Figure 1, where the loop directions (black arrows) and magnetic field directions (red arrows) are depicted at three loop segment positions, for two orthogonal projections (Figure 1, left panels). The root-mean-square value of all misalignment angles for each loop segment (i) and loop (j) is then minimized in the forward-fitting procedure to find the best NLFFF approximation,
In addition to the 3D misalignment angle μ3, we can also define a 2D misalignment angle μ2 with the same Equations (9) and (10), except that the magnetic field vectors Btheo(x) and loop vectors Btheo(x) are a function of two-dimensional space coordinates x = (xij, yij), as they are seen in the 2D projection onto the (x, y)-plane (Figure 1, bottom right panel).
If we have only 2D loop coordinates available (in the case without stereoscopy), we can only forward fit the field lines parameterized with an NLFFF code by minimizing the 2D misalignment angle μ2, because the third space variable (zij) is not available. Our strategy here is to calculate the 2D misalignment angle μ2 in each position for an array of altitudes, ri = 1 + hmax(i/nh), i = 1, ..., nh, that covers a limited altitude range [1 < r < (1 + hmax)] or search volume in which we expect coronal loops to be detectable (Figure 1, top right panel). For EUV images, loops are generally detectable within one density scale height, for which we use a height range of [1 < r < 1.15] solar radii here. This yields multiple (nh) misalignment angles for each 2D loop position (xij, yij), which are shown in the (x–z)-plane (Figure 1, top right panel) and (x–y)-plane (Figure 1, bottom right panel). Our strategy is then to estimate the unknown third zij-coordinate in each 2D position (xij, yij) from that height ri that shows the smallest 2D misalignment angle μ2, and can then proceed with the forward-fitting procedure like in the case of 3D stereoscopic data (xij, yij, zij), which is described in detail in Aschwanden & Malanushenko (2012).
A side effect of the generalized 2D method is that the parameter space is enlarged by an additional dimension, i.e., the unknown zi or ri coordinate of the altitude of each loop position (xij, yij), which in principle increases the computation time by a linear factor of the number nh of altitude levels. However, we optimized the code by vectorization and by organizing the optimization of the altitude variables zij (for each loop segment j and loop i) and the αk-parameter variables (for each magnetic charge m) in an interleaved mode, so that the computation time is reduced by 1–2 orders of magnitude compared with earlier versions of the code (Aschwanden & Malanushenko 2012; Aschwanden et al. 2012b), without loss of accuracy.
When comparing 3D with 2D misalignment angles, we have to be aware that the unknown third dimension zij at every loop position (xij, yij) is handled differently in the two methods. In our new 2D method, we optimize the third coordinate zij by minimizing the 2D misalignment angle independently at every loop segment position (xij, yij), interleaved with the optimization of the nonlinear force-free αm parameter. In the (stereoscopic) 3D method, the third coordinate zij is used as a fixed constraint like the observables (xij, yij). However, we can calculate the median 2D misalignment angle μ2 with both methods, while the 3D misalignment angle μ3 is only defined for the (stereoscopic) 3D-fit method, but not for the (loop-tracing) 2D-fit method.
4. TEST RUNS WITH SIMULATED DATA
In order to test the numerical convergence behavior, the uniqueness of the solutions, and the accuracy of the method in terms of misalignment angles we test our code first with simulated data. We simulate six cases that correspond to the same six nonpotential cases presented in Aschwanden & Malanushenko (2012; cases 7–12), consisting of a unipolar case (N7; Figure 2, top), a dipolar case (N8; Figure 2, middle), a quadrupolar case (N9; Figure 2, bottom), and three decapolar cases with 10 randomly buried magnetic charges each (N10, N11, N12; Figure 3). In each case, we run both the 3D-fitting code (mimicking the availability of stereoscopic loop data), as well as the 2D-fitting code (corresponding to loop tracings from a single EUV image without stereoscopic information). Thus, in the 3D-fitting code the target loops are parameterized with 3D data [xi(s), yi(s), zi(s)], while we ignore the third coordinate in the 2D-fitting code and fit only the 2D coordinates [xi(s), yi(s)] of the target loops. The results are shown in Figures 2 and 3, with the 2D fits in the left-hand panels, and the 3D fits in the right-hand panels.
The convergence of the 2D code can be judged by comparing with the previously tested 3D code (Aschwanden & Malanushenko 2012). We list a summary of the results in Table 1. The 3D fits achieve a mean 3D misalignment angle of μ3 = 36 ± 19 and a 2D misalignment angle of μ2 = 20 ± 11. This is reasonable and slightly better than the results of an earlier version of the 3D code (μ3 = 51 ± 43 in Table 4 of Aschwanden & Malanushenko 2012). In comparison, our new 2D-fitting code achieves an even better agreement with μ2 = 12 ± 05 for the 2D-misalignment angle μ2, while the 3D-misalignment angle μ3 is not defined for the 2D code (due to the lack of line-of-sight coordinates z). The better performance of the 2D code is due to the smaller number of constraints (i.e., 200 constrains for 2D-loop coordinates of 10 loops with 10 segments, compared with 300 constraints for the 3D-fit method). The smaller number of free parameters generally improves the accuracy of the solution. The uniqueness of the solution can be best expressed by the mean misalignment angle μ2. Strictly speaking, an NLFFF solution would only be unique if the number of free parameters (i.e., the number of magnetic charges Nm in our case) matches the number of constraints (i.e., the number of fitting positions (i.e., the product of the number of loops Nloop times the number of fitted segment positions Nsegm, i.e., Nfit = Nloop×Nseg in our case). Moreover, the LOS magnetogram is approximated by a number of magnetic charges Nm that neglects weak magnetic sources, and there are residuals in the decomposition of magnetic charges that contribute to the noise or uncertainty and non-uniqueness of the solutions. Therefore, the uniqueness of an NLFFF solution can best be specified by an uncertainty measure for each field line, which can be quantified either by the misalignment angle μ2 or by a maximum transverse displacement of a field line, i.e., Δx ≈ (L/2)tan μ2 for a particular field line with full length L.
Table 1. Median 2D and 3D Misalignment Angles μ2 and μ3 and Free Energy Ratios ENP/EP for Six Simulated Magnetic Field Configurations (see Figures 2 and 3) Using the 3D-fit Code Suitable for Stereoscopic Data (Marked with 'STEREO') and Using the 2D-fit Code Suitable for 2D Loop Data (Marked with 'Tracing')
|Data Set||CPU (s)||STEREO||STEREO||Tracing||STEREO||Tracing|
|Mean||8.7 ± 7.8||3.6 ± 1.9||2.0 ± 1.1||1.2 ± 0.5|
Note. The computation time (CPU) is listed in seconds.
Download table as:
We show in Table 1 also the ratios ENP/EP of the nonpotential to the potential energies for the 3D and 2D fit methods, which agree within an accuracy of order ≈0.1%–2.5%. Note that the simulated data have 1, 2, 4, and 10 magnetic charges, which corresponds to the number of free parameters in the fit, while we used the double number of magnetic charges in the decomposition of the simulated LOS magnetogram, in order to make the parameterization of the fitted model somewhat different from the target model. Nevertheless, although the 2D solutions have high accuracy (with a mean misalignment of μ2 = 12 ± 05), we have to be aware that the simulated data and the forward-fitting code use the same parameterization of nonlinear α parameters, which warrants a higher accuracy in forward fitting than real data with an unknown parameterization. Hence, we test the code with real solar data in the following section.
5. OBSERVATIONS AND RESULTS
To test the feasibility, fidelity, and accuracy of the new analytical NLFFF forward-fitting code based on 2D tracing of loops (rather than 3D stereoscopy) we are using the same observations for which either stereoscopic 3D reconstruction has been attempted earlier (Aschwanden et al. 2008b, 2008c, 2009, 2012b; Sandman et al. 2009; DeRosa et al. 2009; Aschwanden & Sandman 2010; Sandman & Aschwanden 2011; Aschwanden 2012a), such as for active regions observed with STEREO on 2007 April 30, May 9, May 19, and December 11, or where 2D loop tracing was performed and documented, such as for an active region observed on 1998 May 19 with TRACE (Aschwanden et al. 2008a). The active region numbers, observing times, spacecraft separation angles, number of traced loops, and maximum magnetic field strengths of these observations are listed in Table 2. In all cases, we used line-of-sight magnetograms Bz(x, y, zphot) from the Michelson Doppler Imager (MDI; Scherrer et al. 1995) on board the Solar and Heliospheric Observatory (SOHO), while EUV images were used either from the Transition Region And Coronal Explorer (TRACE; Handy et al. 1999), or from the Extreme Ultraviolet Imager (EUVI; Wülser et al. 2004) on board the STEREO spacecraft A(head) and B(ehind).
Table 2. Data Selection of Four Active Regions Observed with SOHO/MDI and STEREO/EUVI or TRACE
|Region||Date||Time EUV||Time MDI||Separation||of EUVI||Field Strength|
|A||10953 (S05E20)||2007 Apr 30||23:00–23:20||22:24||STEREO 60||200||[−3134,+1425]|
|B||10955 (S09E24)||2007 May 9||20:30–20:50||20:47||STEREO 71||70||[−2396,+1926]|
|C||10953 (N03W03)||2007 May 19||12:40–13:00||12:47||STEREO 86||100||[−2056,+2307]|
|D||10978 (S09E06)||2007 Dec 11||16:30–16:50||14:23||STEREO 427||87||[−2270,+2037]|
|E||8222 (N22W30)||1998 May 19||22:21–22:22||20:48||TRACE (manu)||201||[−1787,+1200]|
|F||8222 (N22W30)||1998 May 19||22:21–22:22||20:48||TRACE (auto)||222||[−1787,+1200]|
Download table as:
We show the results of the NLFFF forward fitting of the six active regions in Figures 4–9, all in the same format, which includes the decomposed line-of-sight magnetogram of SOHO/MDI (gray scale in center of Figures 4–9), the stereoscopically triangulated or visually traced loops (blue curves in Figures 4–9), and the best-fit magnetic field lines (red curves for the segments covered by the observed loops, and in orange for complementary loop parts (although truncated at a height of 0.15 solar radii)). The orthogonal projections of the best-fit magnetic field lines are also shown in the right-hand and top panels of Figures 4–9, as well as the histograms of 2D (μ2) and 3D (μ3) misalignment angles (in the bottom panels of Figures 4–9). The median misalignment angles μ2 and μ3 are also listed in Table 3, for both the previous stereoscopic reconstructions (Aschwanden et al. 2012b; Aschwanden 2012b), marked with 'STEREO' in Table 3, and based on 2D loop tracing in the present study, marked with 'Tracing' in Table 3.
Table 3. Median 2D and 3D Misalignment Angles μ2 and μ3 and Free Energy Ratios ENP/EP for Six Active Regions, Including Four Stereoscopically Triangulated Loops Observed with STEREO (A–D) and Two Cases with Traced Loops Observed with TRACE (E and F)
|Data Set||CPU (s)||STEREO||STEREO||Tracing||STEREO||Tracing|
|(A) 2007 Apr 30||1103||214||46||44||1.006||1.007|
|(B) 2007 May 9||208||179||37||24||1.023||1.009|
|(C) 2007 May 19||415||221||59||36||1.085||1.053|
|(D) 2007 Dec 11||390||147||16||15||1.044||1.026|
|(E) 1998 May 19 (manu)||631||12|
|(F) 1998 May 19 (auto)||915||33|
|(A–F) Mean||610 ± 342||19° ± 3°||40 ± 18||27 ± 13|
Notes. All cases have also been modeled by forward fitting of the NLFFF model to the traced 2D loop coordinates (columns marked with 'Tracing'). The computation time (CPU) is listed in seconds.
Download table as:
The active region A (2007 April 30; Figure 4) shows a lack of stereoscopically triangulated loops in the core of the active region (due to the high level of confusion for loop tracing over the 'mossy' regions), where the highest shear and degree of non-potentiality is expected, and thus deprives us from measuring the largest amount of free magnetic energy, while standard NLFFF codes have stronger constraints in these core regions. In all five active regions we see sunspots with strong magnetic fields, but since we limit our NLFFF solutions to an altitude range of ≲ 0.15 solar radii, we cannot see whether the diverging field lines above the sunspots are open or closed field lines. In principle, we could display our NLFFF solutions to larger altitudes to diagnose where open and closed field regions are, but the accuracy of reconstructed field lines is expected to decrease with height with our method, especially because of the second-order approximation that can represent helical twist well for vertical segments of loops (near the surface), but not so well for horizontal segments of loops (in large altitudes).
5.2. Misalignment Angle Statistics
We see in Table 3 that the 3D misalignment angles μ3 have a mean value of μ3 = 19° ± 3° for the stereoscopy method. If we compare the 2D misalignment angles μ2 between the two methods in Table 3, we see that the STEREO method has a somewhat larger mean, μ = 40 ± 18, than the 2D loop tracing method, with μ = 27 ± 13. This is an effect of the optimization of the 2D misalignment angles in the 2D loop tracing method, where the third coordinate is a free variable and thus has a larger flexibility to find an appropriate model field line with a small 2D misalignment angle, in contrast to the stereoscopic method, where the third coordinate of every observed loop is entirely fixed and leaves less room in the minimization of the misalignment angles between observed loops and theoretical field models. Ideally, if stereoscopy would work perfectly, the third coordinate should be sufficiently accurate so that the best field solution can be found easier with fewer free parameters. However, reality apparently reveals that there is a significant stereoscopic error that hinders optimum field fitting, which is non-existent in the 2D forward-fitting situation. Actually, from the two stereoscopic misalignment angles we can estimate the stereoscopic error σSE, assuming isotropic errors. Thus, defining the 2D misalignment angle as μ22 = σx2 + σ2y, and the 3D misalignment angle as μ23 = σx2 + σ2y + σz2 + σ2SE, with isotropic errors σx = σy = σz, we expect
which yields a stereoscopic error of σSE ≈ 20°. This is somewhat larger than estimated earlier from the parallelism of stereoscopically triangulated loops with close spatial proximity, which amounted to σSE ≈ 7°–12° (Aschwanden & Sandman 2010). Thus, both estimates assess a substantial value to the stereoscopic error that exceeds the accuracy of the best-fit NLFFF solution based on 2D loop tracing (μ2 = 27 ± 13) by far. Nevertheless, both best-fit misalignment angles μ2 are consistent with each other for the two forward-fitting methods. This forward-fitting experiment thus demonstrates that we obtain equally accurate NLFFF fits to coronal data with or without stereoscopy, and thus makes the new method extremely useful.
5.3. Spatial Resolution of Instruments
For the 2D loop tracing method of stereoscopic data (cases A, B, C, and D in Table 1) we just ignored information on the z-coordinate of loops in the 2D forward-fitting algorithm. For cases E and F, we compare visual tracing of loops (case E) with automated loop tracing (case F). Interestingly, we find the most accurate forward fit for case E, which has a 2D misalignment error of only μ2 = 12 (Figure 8), which may have resulted from the higher spatial resolution of loop tracing using TRACE data (with a pixel size of 05) in case E, compared with the three times lower resolution of STEREO (with a pixel size of 16) in cases A, B, C, and D.
5.4. Visual versus Automated Loop Detection
Comparing visually traced (Figure 8) versus automatically traced loops (Figure 9), we note that automated tracing leads to a less accurate NLFFF fit, with μ3 = 33 versus μ3 = 12 for visual tracing. Apparently, the automated loop tracing method (Aschwanden et al. 2008a) can easily be misguided or side-tracked by about Δμ ≈ 2° onto near-cospatial loops with similar coherent large curvature. Comparing individual field lines in Figure 8 with Figure 9, we see also that the automated loop tracing algorithm produces a number of short loop segments with large misalignment angles, which are obviously a weakness of the automated tracing code, resulting into a larger misalignment error for the fitted NLFFF solution (which is a factor of ≈3 larger in the average for this case).
5.5. Magnetic Field Strengths
In Figure 10, we show scatterplots of the magnetic field strengths B3 retrieved at the photospheric level for the stereoscopic 3D method versus the field strength B2 obtained from the 2D loop tracing method. The analytical NLFFF forward-fitting algorithm starts first a decomposition of point-like buried magnetic charges from the line-of-sight component Bz(x, y, zphot), but has no constraints for the transverse components Bx and By, except for the coronal loop coordinates. Thus, the directions of coronal loops near the footpoints will determine the transverse components Bx and By in the photosphere. The scatterplots of B2 versus B3 in Figure 10 show that the smallest scatter between the two methods appears for active region B, which is the one closest to a potential field (although we are not able to trace and triangulate loops in the core of the active region, where supposedly the highest level of non-potentiality occurs). The linear regression fits show a good correspondence in the order of ≈0.1%–1.0% between the two methods, which results from different fitting criteria of the nonpotential field components. A ratio of Bφ/Br ≈ 10% in the azimuthal field (or nonpotential) field component Bφ (Equation (5)) with respect to the potential field component Br (Equation (4)) would result into a change of ≈0.5% in the magnetic field strength.
5.6. Magnetic Energies
The free magnetic energy, which is the difference between the nonpotential and the potential field energy, integrated over the spatial volume of an active region, has been calculated for the stereoscopic forward-fitting method in Aschwanden (2012b, Table 1 therein). Here we also calculate these quantities for the 2D loop tracing method, both methods being juxtaposed in Table 3 and in Figure 11. The absolute values of the potential field energy EP agree within a factor of 1.04 ± 0.02 between the two methods, while the free energies show differences of 1%–2% of the potential energy. If we consider the difference in the free energy between two methods as a measure of a systematic error, we assess an uncertainty of Efree/EP − 1 = 0.016 ± 0.013, or about ±1.5%.
6.1. Assessment of Various NLFFF Methods
Quantitative comparisons of various nonlinear force-free field (NLFFF) calculation methods applied to coronal volumes that encompass an active region included optimizational, magneto-frictional, Grad–Rubin-based, and Green's-function-based models (Schrijver et al. 2006, 2008; DeRosa et al. 2009). A critical assessment of 11 NLFFF methods revealed significant differences in the non-potential magnetic energy and in the degree of misalignment (μ3 ≈ 24°–44°) with respect to stereoscopically triangulated coronal loops. The chief problem responsible for these discrepancies was identified in terms of the non-force freeness of the lower (photospheric) boundary, too small boundary areas, and uncertainties of the boundary data (DeRosa et al. 2009). Obviously, the non-force freeness of the photosphere can only be circumvented by incorporating coronal magnetic field geometries into NLFFF extrapolations, such as the information from stereoscopically triangulated loops.
Implementation of coronal magnetic field data into NLFFF codes is not straightforward, because most of the conventional NLFFF codes are designed to extrapolate from a lower boundary in the upward direction, and thus the volume-filling coronal information cannot be treated as a boundary problem. A natural method to match arbitrary constraints that are not necessarily a boundary problem is a forward-fitting approach, which, however, requires a parameterization of a magnetic field model. Since NLFFF models are implicitly defined by the two differential equations of force freeness and divergence freeness, an explicit parameterization of a magnetic field is not trivial. There are essentially three approaches that have been attempted so far: (1) preprocessing of magnetic boundary data by minimizing the Lorentz force at multiple spatial scales (Wiegelmann 2004; Jing et al. 2010), (2) linear force-free fitting with subsequent relaxation to a nonlinear force-free field (Malanushenko et al. 2009, 2012), and (3) magnetic field parameterization with buried magnetic charges and forward fitting of an approximative analytical NLFFF solution (Aschwanden 2012a, 2012b; Aschwanden & Malanushenko 2012). The third method requires stereoscopically triangulated 3D coordinates of coronal loops, which have been successfully modeled in four different active regions, improving the misalignment angles between the observed loops and a potential field (μ3 = 19°–46°) by a factor of about two, when fitted with the analytical NLFFF solution (μ3 = 14°–19°; Aschwanden et al. 2012b). However, although this method provides a more realistic NLFFF model, the main restriction is the availability of STEREO data, as well as the inferior spatial resolution of STEREO/EUVI compared with other EUV imagers (e.g., TRACE, or SDO/AIA). In order to circumvent this restriction we developed a generalized code in this study that requires only a high-resolution EUV image, from which a sufficiently large number of coronal loops can be traced in two dimensions, plus a line-of-sight magnetogram, which were essentially always available since the launch of the SOHO mission in 1995. A fortunate outcome of this study is that the reduction of 3D information from coronal loops to 2D coordinates does not handicap the accuracy of an NLFFF forward-fitting method, as we demonstrated here. One potential limitation of forward fitting of analytical NLFFF approximations may be the accuracy for strong nonpotential cases (say near flaring times), since our analytical NLFFF approximation is only accurate to the second order of the force-free α-parameter. Another caveat may be the universality of the analytical NLFFF parameterization. Our approximation is designed to model azimuthal (rotational) twist around a vertical axis with respect to the solar surface, which corresponds to electric currents flowing in the vertical direction. This geometry may not be appropriate for horizontally twisted structures, such as horizontally extended filaments. For such cases, a more general NLFFF solution could be desirable. In this respect, more general NLFFF solutions such as those currently being developed by Malanushenko et al. (2009, 2012), may provide new tools, after they are generalized for a spherical solar surface and optimized for computational speed. Nevertheless, our analytical NLFFF forward-fitting algorithm can always be used to find a quick first approximation of an NLFFF solution, which then could be refined by alternative (more time-consuming) NLFFF codes.
6.2. Computational Speed
Note that our forward-fitting code calculates one trial magnetic field configuration in about 0.01 s, while some 104–105 iterations are needed for the convergence of an approximate NLFFF solution, accumulating to a few minutes total computation time. Our forward-fitting code computes a solution to an active region with an average computation time of ≈10 minutes (see CPU times in Table 3) on a Mac OS X, 2 × 3.2 GHz Quad-Core Intel Xeon, 32 GB memory, 800 MHz DDR2 FB-DIMM computer. Simple cases that can be represented with Nm ≲ 20 magnetic charges, require only computation times of order ≲ 10 s (see CPU times in Table 1).
We developed an analytical NLFFF forward-fitting code that requires only the input of a line-of-sight magnetogram Bz(x, y, zphot) and a set of 2D loop tracing coordinates (xi, yi) that can be obtained from any high-resolution EUV image. This is the first attempt to measure the coronal magnetic field based on directly observed 2D images alone, while we needed stereoscopically triangulated 3D loop coordinates (xi, yi, zi) from STEREO in previous studies. There exists only one other study (to our knowledge) that heads into the same direction of forward fitting an NLFFF model to 2D data (Malanushenko et al. 2012), using simulated 2D data based on analytical NLFFF solutions (from Low & Lou 1990) or previous NLFFF modeling (from Schrijver et al. 2008).
We forward-fitted analytical NLFFF approximations to four active regions (cases A, B, C, D) using 2D loop tracings only and compared them with previous NLFFF fits using (stereoscopic) 3D loop coordinates. Furthermore, we forward-fitted analytical NLFFF approximations to another active region using visual 2D loop tracings (case E) or automatically traced 2D loop coordinates (case F). Our findings from these exercises are the following.
Download SOHO Notes For Mac 9.3.7 Pro
Download SOHO Notes For Mac 9.3.7 Torrent
- Our forward-fitting experiment with two different methods demonstrated that a nonlinear force-free magnetic field (NLFFF) solution B(x) can be obtained with equal accuracy with or without stereoscopy. This result relinquishes the necessity of STEREO data for future magnetic modeling of active regions on the solar disk, but the availability of suitable STEREO data was crucial to establish this result.
- The accuracy of a forward-fitted NLFFF approximation that includes vertical currents (with twisted azimuthal magnetic field components) matches coronal loops observed in EUV with a median 2D misalignment angle of μ2 = 27 ± 13, while stereoscopic 3D data exhibit a commensurable 2D misalignment angle, but a substantially larger 3D misalignment angle (μ3 = 19° ± 3°), which implies stereoscopic measurement errors in the order of σSE ≈ 20°. These substantial stereoscopic measurement errors lead to less accurate NLFFF fits than 2D loop tracings with unconstrained line-of-sight positions.
- 2D loop tracings in high-resolution images (i.e., 05 pixels with TRACE) lead to more accurate NLFFF fits (with a misalignment angle of μ2 = 12 in case E) than images with lower spatial resolution (i.e., 16 pixels with STEREO), yielding μ2 = 40 ± 18 for the cases A, B, C, and D.
- Visually (or manually) traced 2D loop coordinates appear to be still superior to the best automated loop tracing algorithms, yielding a misalignment angle of μ2 = 12 (in case E) versus μ2 = 33 (in case F).
- Magnetic field strengths B of best-fit NLFFF approximations are retrieved with an accuracy of ≈0.1%–1.0%, comparing the 2D loop tracing method with the stereoscopic 3D triangulation method.
- Magnetic energies differ by a factor of ≈4% ± 2% between the two methods, while the free energy has a systematic error of ≈ ± 1.5% (of the total magnetic energy) between the two methods, which is about an order of magnitude smaller than found between other (standard) NLFFF codes.
- The computational speed of our NLFFF code allows the computation of a space-filling magnetic field configuration of an active region in about 0.01 s, while forward fitting to a set of coronal loops is feasible with about 104–105 iterations, requiring a total computation time of a few minutes.
These results demonstrate clearly that we can perform accurate NLFFF magnetic modeling based on 2D loop tracings, which will relinquish the need of STEREO data in future, at least for active regions near the solar disk center. The analyzed active regions extended up to 0.65 solar radii away from disk center, so in principle we can perform NLFFF modeling for at least half of the number of active regions observed on the solar disk, especially since the parameterization of our analytical NLFFF approximation takes the sphericity of the solar surface fully into account (which is not the case in most other NLFFF codes). For future developments we expect that other NLFFF codes could implement minimization of the misalignment angle μ2 with coronal loops in parallel to the optimization of force freeness and divergence freeness, rather than by preprocessing of vector boundary data. Further improvements in automated 2D loop tracings could replace visual/manual tracing methods and this way render NLFFF modeling with coronal constraints in a fully automated way.
The author appreciates the constructive comments by an anonymous referee and helpful discussions with Allen Gary, Anna Malanushenko, Marc DeRosa, and Karel Schrijver. Part of the work was supported by NASA contract NNG 04EA00C of the SDO/AIA instrument and the NASA STEREO mission under NRL contract N00173-02-C-2035.
Please read the following tips BEFORE trying to migrate:
- TIP #1: In order to migrate your data from SOHO Notes to Notelife, you MUST be running SOHO Notes 9.3.7 and not an earlier version. If you’re not using SOHO Notes 9.3.7 and you try to migrate, you’ll receive an error message that says, “The managed object model version used to open the persistent store is incompatible with the one that was used to create the persistent store”. You can download 9.3.7 here:
- TIP #2: When you migrate your data, you’ll be asked to locate the folder containing your SOHO Notes database, and not the database file itself. The actual database file wlll be disabled (grayed out). Locating and choosing this folder rather than the database is necessary because Notelife is a sandboxed application and macOS requires that you grant access to the entire folder so that all of the database-related files in the folder can be accessed.
- TIP #3: The migration process requires that all notes in SOHO Notes be unencrypted first since Notelife will have no ability to decrypt them. So if you have encrypted notes, launch SOHO Notes and unencrypt all of them first.
- TIP #4: The following fields will NOT migrate since Notelife doesn’t have an equivalent field: comments, labels, shortcuts, alarms, margin settings, and task/contact attachments.
- TIP #5: Notelife doesn't have the notion of a 'note type' so an 'image note' in SOHO Notes will migrate over as a regular note in Notelife with an image attachment. Similarly, the contents of a 'form' note in SOHO Notes will migrate over as a regular note in Notelife with the contents of the form converted to text and stored in the note body in Notelife.
Onenote For Mac
After you've read the above tips, please do the following to migrate:
Download SOHO Notes For Mac 9.3.7 Full
- Quit SOHO Notes
- Launch Notelife
- Choose File > Migrate SOHO Notes and follow the onscreen steps.