Source author record

George H. Fisher

George H. Fisher appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

15works
2topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

15 published item(s)

preprint2020arXiv

The PDFI_SS Electric Field Inversion Software

We describe the PDFI_SS software library, which is designed to find the electric field at the Sun's photosphere from a sequence of vector magnetogram and Doppler velocity measurements, and estimates of horizontal velocities obtained from local correlation tracking using the recently upgraded FLCT code. The library, a collection of Fortran subroutines, uses the "PDFI" technique described by Kazachenko et al. (2014), but modified for use in spherical, Plate-Carrée geometry on a staggered grid. The domain over which solutions are found is a subset of the global spherical surface, defined by user-specified limits of colatitude and longitude. Our staggered-grid approach, based on that of Yee (1966), is more conservative and self-consistent compared to the centered, Cartesian grid used by Kazachenko et al. (2014). The library can be used to compute an end-to-end solution for electric fields from data taken by the HMI instrument aboard NASA's SDO Mission. This capability has been incorporated into the HMI pipeline processing system operating at SDO's JSOC. The library is written in a general and modular way so that the calculations can be customized to modify or delete electric field contributions, or used with other data sets. Other applications include "nudging" numerical models of the solar atmosphere to facilitate assimilative simulations. The library includes an ability to compute "global" (whole-Sun) electric field solutions. The library also includes an ability to compute Potential Magnetic Field solutions in spherical coordinates. This distribution includes a number of test programs which allow the user to test the software.

preprint2016arXiv

Deriving Potential Coronal Magnetic Fields from Vector Magnetograms

The minimum-energy configuration for the magnetic field above the solar photosphere is curl-free (hence, by Ampere's 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 Laplace's 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.

preprint2015arXiv

Emergence of Anchored Flux Tubes Through the Convection Zone

We model the evolution of buoyant magnetic flux tubes in the Sun's convection zone. A flux tube is assumed to lie initially near the top of the stably stratified radiative core below the convection zone, but a segment of it is perturbed into the convection zone by gradual heating and convective overshoot motions. The ends ("footpoints") of the segment remain anchored at the base of the convection zone, and if the segment is sufficiently long, it may be buoyantly unstable, rising through the convection zone in a short time. The length of the flux tube determines the ratio of buoyancy to magnetic tension: short loops of flux are arrested before reaching the top of the convection zone, while longer loops emerge to erupt through the photosphere. Using Spruit's convection zone model, we compute the minimum footpoint separation $L_c$ required for erupting flux tubes. We explore the dependence of $L_c$ on the initial thermal state of the perturbed flux tube segment and on its initial magnetic field strength. Following an investigation of thermal diffusion time scales and the dynamic rise times of unstable flux tube segments, we conclude that the most likely origin for magnetic flux which erupts to the surface is from short length scale perturbations ($L < L_c$) which are initially stable, but which are subsequently destabilized either by diffusion of heat into the tube or by stretching of the anchor points until $L$ just exceeds $L_c$. In either case, the separation of the anchor points of the emergent tube should lie between the critical distance for a tube in mechanical equilibrium and one in thermal equilibrium. Finally, after comparing the dispersion of dynamic rise times with the much shorter observed active region formation time scales, we conclude that active regions form from the emergence of a single flux tube segment.

preprint2015arXiv

Photospheric Electric Fields and Energy Fluxes in the Eruptive Active Region NOAA 11158

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-Doppler-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.6\times10^{32}$ ergs, which is partitioned as $2.0$ and $8.6\times10^{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.

preprint2015arXiv

Sign singularity and flares in solar active region NOAA 11158

Solar Active Region NOAA 11158 has hosted a number of strong flares, including one X2.2 event. The complexity of current density and current helicity are studied through cancellation analysis of their sign-singular measure, which features power-law scaling. Spectral analysis is also performed, revealing the presence of two separate scaling ranges with different spectral index. The time evolution of parameters is discussed. Sudden changes of the cancellation exponents at the time of large flares, and the presence of correlation with EUV and X-ray flux, suggest that eruption of large flares can be linked to the small scale properties of the current structures.

preprint2015arXiv

The Coronal Global Evolutionary Model: Using HMI Vector Magnetogram and Doppler Data to Model the Buildup of Free Magnetic Energy in the Solar Corona

The most violent space weather events (eruptive solar flares and coronal mass ejections) are driven by the release of free magnetic energy stored in the solar corona. Energy can build up on timescales of hours to days, and then may be suddenly released in the form of a magnetic eruption, which then propagates through interplanetary space, possibly impacting the Earth's space environment. Can we use the observed evolution of the magnetic and velocity fields in the solar photosphere to model the evolution of the overlying solar coronal field, including the storage and release of magnetic energy in such eruptions? The objective of CGEM, the Coronal Global Evolutionary Model, funded by the NASA/NSF Space Weather Modeling program, is to develop and evaluate such a model for the evolution of the coronal magnetic field. The evolving coronal magnetic field can then be used as a starting point for magnetohydrodynamic (MHD) models of the corona, which can then be used to drive models of heliospheric evolution and predictions of magnetic field and plasma density conditions at 1AU.

preprint2015arXiv

Why Is the Great Solar Active Region 12192 Flare-Rich But CME-Poor?

Solar active region (AR) 12192 of October 2014 hosts the largest sunspot group in 24 years. It is the most prolific flaring site of Cycle 24, but surprisingly produced no coronal mass ejection (CME) from the core region during its disk passage. Here, we study the magnetic conditions that prevented eruption and the consequences that ensued. We find AR 12192 to be "big but mild"; its core region exhibits weaker non-potentiality, stronger overlying field, and smaller flare-related field changes compared to two other major flare-CME-productive ARs (11429 and 11158). These differences are present in the intensive-type indices (e.g., means) but generally not the extensive ones (e.g., totals). AR 12192's large amount of magnetic free energy does not translate into CME productivity. The unexpected behavior suggests that AR eruptiveness is limited by some relative measure of magnetic non-potentiality over the restriction of background field, and that confined flares may leave weaker photospheric and coronal imprints compared to their eruptive counterparts.

preprint2014arXiv

A Comprehensive Method of Estimating Electric Fields from Vector Magnetic Field and Doppler Measurements

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 models 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 $(θ=0$) and less than $25%$ and $10%$ errors respectively at large viewing angles ($θ<60^\circ$). We compare our results with other inversion methods at zero viewing angle, and find that our method's 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.

preprint2013arXiv

A Magnetic Calibration of Photospheric Doppler Velocities

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 velocities 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.

preprint2011arXiv

Can we Determine Electric Fields and Poynting Fluxes from Vector Magnetograms and Doppler Measurements?

The availability of vector magnetogram sequences with sufficient accuracy and cadence to estimate the time derivative of the magnetic field allows us to use Faraday's 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 -- Faraday's 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.

preprint2011arXiv

Global Forces in Eruptive Solar Flares: The Lorentz force Acting on the Solar Atmosphere and the Solar Interior

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 opposite 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.

preprint2011arXiv

Radiative Cooling in MHD Models of the Quiet Sun Convection Zone and Corona

We present a series of numerical simulations of the quiet Sun plasma threaded by magnetic fields that extend from the upper convection zone into the low corona. We discuss an efficient, simplified approximation to the physics of optically thick radiative transport through the surface layers, and investigate the effects of convective turbulence on the magnetic structure of the Sun's atmosphere in an initially unipolar (open field) region. We find that the net Poynting flux below the surface is on average directed toward the interior, while in the photosphere and chromosphere the net flow of electromagnetic energy is outward into the solar corona. Overturning convective motions between these layers driven by rapid radiative cooling appears to be the source of energy for the oppositely directed fluxes of electromagnetic energy.

preprint2010arXiv

Improving Large-scale Convection Zone-to-Corona Models

We introduce two new methods that are designed to improve the realism and utility of large, active region-scale 3D MHD models of the solar atmosphere. We apply these methods to RADMHD, a code capable of modeling the Sun's upper convection zone, photosphere, chromosphere, transition region, and corona within a single computational volume. We first present a way to approximate the physics of optically-thick radiative transfer without having to take the computationally expensive step of solving the radiative transfer equation in detail. We then briefly describe a rudimentary assimilative technique that allows a time series of vector magnetograms to be directly incorporated into the MHD system.

preprint2010arXiv

On the Mechanical Energy Available to Drive Solar Flares

Where does solar flare energy come from? More specifically, assuming that the ultimate source of flare energy is mechanical energy in the convection zone, how is this translated into energy dissipated or stored in the corona? This question appears to have been given relatively little thought, as attention has been focussed predominantly on mechanisms for the rapid dissipation of coronal magnetic energy by way of MHD instabilities and plasma micro instabilities. We consider three types of flare theory: the steady state "photospheric dynamo" model in which flare power represents coronal dissipation of currents generated simultaneously by sub-photospheric flows; the "magnetic energy storage" model where sub-photospheric flows again induce coronal currents but which in this case are built up over a longer period before being released suddenly; and "emerging flux" models, in which new magnetic flux rising to the photosphere already contains free energy, and does not require subsequent stressing by photospheric motions. We conclude that photospheric dynamos can power only very minor flares; that coronal energy storage can in principle meet the requirements of a major flare, although perhaps not the very largest flares, but that difficulties in coupling efficiently to the energy source may limit this mechanism to moderate sized flares; and that emerging magnetic flux tubes, generated in the solar interior, can carry sufficient free energy to power even the largest flares ever observed.

preprint2010arXiv

Radiation Hydrodynamics in Solar Flares

Solar flares are currently understood as the explosive release of energy stored in the form of stressed magnetic fields. In many cases, the released energy seems to take the form of large numbers of electrons accelerated to high energies or alternatively plasma heated to very high temperatures. The transport of this energy into the remaining portion of the atmosphere results in violent mass motion and strong emission across the electromagnetic spectrum. One important phenomenon observed during flares is the appearance in coronal magnetic loops of large amounts of upflowing, soft X-ray emitting plasma. It is believed that this is due to chromospheric evaporation, the process of heating cool chromospheric material beyond its ability to radiate. The pressure increase in the evaporated plasma leads to a number of interesting phenomena in the flare chromosphere. The sudden pressure increase initiates a downward moving "chromospheric condensation", an overdense region which gradually decelerates as it accretes material and propagates into the gravitationally stratified chromosphere. Solutions to an equation of motion for this condensation shows that its motion decays after about one minute of propagation into the chromosphere. When the front of this downflowing region is supersonic relative to the atmosphere ahead of it, a radiating shock will form. If the downflow is rapid enough, the shock strength should be sufficient to excite UV radiation normally associated with the transition region, and furthermore, the radiating shock will be brighter than the transition region. These results lead to a number of observationally testable relationships between the optical and ultraviolet spectra from the condensation and radiating shock.