How much electromagnetic energy crosses the photosphere in evolving solar active regions? With the advent of high-cadence vector magnetic field observations, addressing this fundamental question has become tractable. In this paper, we apply the PTD-D
oppler-FLCT-Ideal (PDFI) electric field inversion technique of Kazachenko et al. (2014) to a 6-day HMI/SDO vector magnetogram and Doppler velocity sequence, to find the electric field and Poynting flux evolution in active region NOAA 11158, which produced an X2.2 flare early on 2011 February 15. We find photospheric electric fields ranging up to $2$ V/cm. The Poynting fluxes range from $[-0.6$ to $2.3]times10^{10}$ ergs$cdot$cm$^{-2}$s$^{-1}$, mostly positive, with the largest contribution to the energy budget in the range of $[10^9$-$10^{10}]$ ergs$cdot$cm$^{-2}$s$^{-1}$. Integrating the instantaneous energy flux over space and time, we find that the total magnetic energy accumulated above the photosphere from the initial emergence to the moment before the X2.2 flare to be $E=10.6times10^{32}$ ergs, which is partitioned as $2.0$ and $8.6times10^{32}$ ergs, respectively, between free and potential energies. Those estimates are consistent with estimates from preflare non-linear force-free field (NLFFF) extrapolations and the Minimum Current Corona estimates (MCC), in spite of our very different approach. This study of photospheric electric fields demonstrates the potential of the PDFI approach for estimating Poynting fluxes and opens the door to more quantitative studies of the solar photosphere and more realistic data-driven simulations of coronal magnetic field evolution.
We study the effect of newly emerged solar active regions (ARs) on the large-scale magnetic environment of pre-existing ARs (PEARs). We first present a theoretical approach to quantify the interaction energy between new ARs and PEARs as the differenc
e between (i) the summed magnetic energies of their individual potential fields and (ii) the energy of their superposed potential fields. We expect that this interaction energy can, depending upon the relative arrangements of newly emerged and PEAR magnetic flux, indicate the existence of topological free magnetic energy in the global coronal field that is independent of any internal free magnetic energy due to coronal electric currents flowing within the newly emerged and PEAR flux systems. We then examine the interaction energy in two well-studied cases of flux emergence, but find that the predicted energetic perturbation is relatively small compared to energies released in large solar flares. Next, we present an observational study on the influence of the emergence of new ARs on flare statistics in PEARs, using NOAAs Solar Region Summary and GOES flare databases. As part of an effort to precisely determine the emergence time of ARs in a large event sample, we find that emergence in about half of these regions exhibits a two-stage behavior, with an initial gradual phase followed by a more rapid phase. Regarding flaring, we find that the emergence of new ARs is associated with a significant increase in the occurrence rate of X- and M-class flares in PEARs. This effect tends to be more significant when PEARs and new emerging ARs are closer. Given the relative weakness of the interaction energy, this effect suggests that perturbations in the large-scale magnetic field, such as topology changes invoked in the breakout model of coronal mass ejections, might play a significant role in the occurrence of some flares.
Researchers have reported (i) correlations of coronal mass ejection (CME) speeds and the total photospheric magnetic flux swept out by flare ribbons in flare-associated eruptive events, and, separately, (ii) correlations of CME speeds and more rapid
decay, with height, of magnetic fields in potential coronal models above eruption sites. Here, we compare the roles of both ribbon fluxes and the decay rates of overlying fields in a set of 16 eruptive events. We confirm previous results that higher CME speeds are associated with both larger ribbon fluxes and more rapidly decaying overlying fields. We find the association with ribbon fluxes to be weaker than a previous report, but stronger than the dependence on the decay rate of overlying fields. Since the photospheric ribbon flux is thought to approximate the amount of coronal magnetic flux reconnected during the event, the correlation of speeds with ribbon fluxes suggests that reconnection plays some role in accelerating CMEs. One possibility is that reconnected fields that wrap around the rising ejection produce an increased outward hoop force, thereby increasing CME acceleration. The correlation of CME speeds with more rapidly decaying overlying fields might be caused by greater downward magnetic tension in stronger overlying fields, which could act as a source of drag on rising ejections.
The minimum-energy configuration for the magnetic field above the solar photosphere is curl-free (hence, by Amperes law, also current-free), so can be represented as the gradient of a scalar potential. Since magnetic fields are divergence free, this
scalar potential obeys Laplaces equation, given an appropriate boundary condition (BC). With measurements of the full magnetic vector at the photosphere, it is possible to employ either Neumann or Dirichlet BCs there. Historically, the Neumann BC was used with available line-of-sight magnetic field measurements, which approximate the radial field needed for the Neumann BC. Since each BC fully determines the 3D vector magnetic field, either choice will, in general, be inconsistent with some aspect of the observed field on the boundary, due to the presence of both currents and noise in the observed field. We present a method to combine solutions from both Dirichlet and Neumann BCs to determine a hybrid, least-squares potential field, which minimizes the integrated square of the residual between the potential and actual fields. This has advantages in both not overfitting the radial field used for the Neumann BC, and maximizing consistency with the observations. We demonstrate our methods with SDO/HMI vector magnetic field observations of AR 11158, and find that residual discrepancies between the observed and potential fields are significant, and are consistent with nonzero horizontal photospheric currents. We also analyze potential fields for two other active regions observed with two different vector magnetographs, and find that hybrid potential fields have significantly less energy than the Neumann fields in every case --- by more than 10^(32) erg in some cases. This has major implications for estimates of free magnetic energy in coronal field models, e.g., non-linear force-free field extrapolations.
Photospheric electric fields, estimated from sequences of vector magnetic field and Doppler measurements, can be used to estimate the flux of magnetic energy (the Poynting flux) into the corona and as time-dependent boundary conditions for dynamic mo
dels of the coronal magnetic field. We have modified and extended an existing method to estimate photospheric electric fields that combines a poloidal-toroidal (PTD) decomposition of the evolving magnetic field vector with Doppler and horizontal plasma velocities. Our current, more comprehensive method, which we dub the {bf P}TD-{bf D}oppler-{bf F}LCT {bf I}deal (PDFI) technique, can now incorporate Doppler velocities from non-normal viewing angles. It uses the texttt{FISHPACK} software package to solve several two-dimensional Poisson equations, a faster and more robust approach than our previous implementations. Here, we describe systematic, quantitative tests of the accuracy and robustness of the PDFI technique using synthetic data from anelastic MHD (texttt{ANMHD}) simulations, which have been used in similar tests in the past. We find that the PDFI method has less than $1%$ error in the total Poynting flux and a $10%$ error in the helicity flux rate at a normal viewing angle $(theta=0$) and less than $25%$ and $10%$ errors respectively at large viewing angles ($theta<60^circ$). We compare our results with other inversion methods at zero viewing angle, and find that our methods estimates of the fluxes of magnetic energy and helicity are comparable to or more accurate than other methods. We also discuss the limitations of the PDFI method and its uncertainties.
The zero point of measured photospheric Doppler shifts is uncertain for at least two reasons: instrumental variations (from, e.g., thermal drifts), and the convective blueshift, a known correlation between intensity and upflows. Calibrated Doppler ve
locities would be useful for (i) improving estimates of the Poynting flux of magnetic energy across the photosphere, and (ii) constraining processes underlying flux cancellation, the mutual apparent loss of magnetic flux in closely spaced, opposite-polarity magnetogram features. We present a method to absolutely calibrate line-of-sight (LOS) velocities in solar active regions (ARs) near disk center using three successive vector magnetograms and one Dopplergram coincident with the central magnetogram. It exploits the fact that Doppler shifts measured along polarity inversion lines (PILs) of the LOS magnetic field determine one component of the velocity perpendicular to the magnetic field, and optimizes consistency between changes in LOS flux near PILs and the transport of transverse magnetic flux by LOS velocities, assuming ideal electric fields govern the magnetic evolution. We apply our method to vector magnetograms of AR 11158, observed by the Helioseismic and Magnetic Imager (HMI) aboard the Solar Dynamics Observatory, and find clear evidence of offsets in the Doppler zero point, in the range of 50 -- 550 m s$^{-1}$. In addition, we note that a simpler calibration can be determined from an LOS magnetogram and Dopplergram pair from the median Doppler velocity among all near-disk-center PIL pixels. We briefly discuss shortcomings in our initial implementation, and suggest ways to address these. In addition, as a step in our data reduction, we discuss use of temporal continuity in the transverse magnetic field direction to correct apparently spurious fluctuations in resolution of the 180$^circ$ ambiguity.
We use autocorrelation to investigate evolution in flow fields inferred by applying Fourier Local Correlation Tracking (FLCT) to a sequence of high-resolution (0.3 arcsec), high-cadence ($simeq 2$ min) line-of-sight magnetograms of NOAA active region
(AR) 10930 recorded by the Narrowband Filter Imager (NFI) of the Solar Optical Telescope (SOT) aboard the {em Hinode} satellite over 12--13 December 2006. To baseline the timescales of flow evolution, we also autocorrelated the magnetograms, at several spatial binnings, to characterize the lifetimes of active region magnetic structures versus spatial scale. Autocorrelation of flow maps can be used to optimize tracking parameters, to understand tracking algorithms susceptibility to noise, and to estimate flow lifetimes. Tracking parameters varied include: time interval $Delta t$ between magnetogram pairs tracked, spatial binning applied to the magnetograms, and windowing parameter $sigma$ used in FLCT. Flow structures vary over a range of spatial and temporal scales (including unresolved scales), so tracked flows represent a local average of the flow over a particular range of space and time. We define flow lifetime to be the flow decorrelation time, $tau$. For $Delta t > tau$, tracking results represent the average velocity over one or more flow lifetimes. We analyze lifetimes of flow components, divergences, and curls as functions of magnetic field strength and spatial scale. We find a significant trend of increasing lifetimes of flow components, divergences, and curls with field strength, consistent with Lorentz forces partially governing flows in the active photosphere, as well as strong trends of increasing flow lifetime and decreasing magnitudes with increases in both spatial scale and $Delta t$.
The availability of vector magnetogram sequences with sufficient accuracy and cadence to estimate the time derivative of the magnetic field allows us to use Faradays law to find an approximate solution for the electric field in the photosphere, using
a Poloidal-Toroidal Decomposition (PTD) of the magnetic field and its partial time derivative. Without additional information, however, the electric field found from this technique is under-determined -- Faradays law provides no information about the electric field that can be derived the gradient of a scalar potential. Here, we show how additional information in the form of line-of-sight Doppler flow measurements, and motions transverse to the line-of-sight determined with ad-hoc methods such as local correlation tracking, can be combined with the PTD solutions to provide much more accurate solutions for the solar electric field, and therefore the Poynting flux of electromagnetic energy in the solar photosphere. Reliable, accurate maps of the Poynting flux are essential for quantitative studies of the buildup of magnetic energy before flares and coronal mass ejections.
We compute the change in the Lorentz force integrated over the outer solar atmosphere implied by observed changes in vector magnetograms that occur during large, eruptive solar flares. This force perturbation should be balanced by an equal and opposi
te force perturbation acting on the solar photosphere and solar interior. The resulting expression for the estimated force change in the solar interior generalizes the earlier expression presented by Hudson, Fisher and Welsch (CS-383, ASP, 221, 2008), providing horizontal as well as vertical force components, and provides a more accurate result for the vertical component of the perturbed force. We show that magnetic eruptions should result in the magnetic field at the photosphere becoming more horizontal, and hence should result in a downward (towards the solar interior) force change acting on the photosphere and solar interior, as recently argued from an analysis of magnetogram data by Wang and Liu (Astrophys. J. Lett. 716, L195, 2010). We suggest the existence of an observational relationship between the force change computed from changes in the vector magnetograms, the outward momentum carried by the ejecta from the flare, and the properties of the helioseismic disturbance driven by the downward force change. We use the impulse driven by the Lorentz-force change in the outer solar atmosphere to derive an upper limit to the mass of erupting plasma that can escape from the Sun. Finally, we compare the expected Lorentz-force change at the photosphere with simple estimates from flare-driven gasdynamic disturbances and from an estimate of the perturbed pressure from radiative backwarming of the photosphere in flaring conditions.
