Showing 1–37 of 37 results for author: Welsch, B T

Large Photospheric Doppler Shift in Solar Active Region 12673: I. Field-Aligned Flows

Authors: Jiayi Liu, Xudong Sun, Peter W. Schuck, Sarah A. Jaeggli, Brian T. Welsch, Carlos Quintero Noda

Abstract: Delta ($δ$) sunspots sometimes host fast photospheric flows along the central magnetic polarity inversion line (PIL). Here we study the strong Doppler shift signature in the central penumbral light bridge of solar active region NOAA 12673. Observations from the Helioseismic and Magnetic Imager (HMI) indicate highly sheared, strong magnetic fields. Large Doppler shifts up to 3.2 km s$^{-1}$ appeare… ▽ More Delta ($δ$) sunspots sometimes host fast photospheric flows along the central magnetic polarity inversion line (PIL). Here we study the strong Doppler shift signature in the central penumbral light bridge of solar active region NOAA 12673. Observations from the Helioseismic and Magnetic Imager (HMI) indicate highly sheared, strong magnetic fields. Large Doppler shifts up to 3.2 km s$^{-1}$ appeared during the formation of the light bridge and persisted for about 16 hours. A new velocity estimator, called DAVE4VMwDV, reveals fast converging and shearing motion along the PIL from HMI vector magnetograms, and recovers the observed Doppler signal much better than an old version of the algorithm. The inferred velocity vectors are largely (anti-)parallel to the inclined magnetic fields, suggesting that the observed Doppler shift contains significant contribution from the projected, field-aligned flows. High-resolution observations from the Hinode/Spectro-Polarimeter (SP) further exhibit a clear correlation between the Doppler velocity and the cosine of the magnetic inclination, which is in agreement with HMI results and consistent with a field-aligned flow of about 9.6 km s$^{-1}$. The complex Stokes profiles suggest significant gradients of physical variables along the line of sight. We discuss the implications on the $δ$-spot magnetic structure and the flow-driving mechanism. △ Less

Submitted 18 July, 2023; originally announced July 2023.

Comments: 17 pages, 14 figures. Accepted for publication in ApJ

Changes of Magnetic Energy and Helicity in Solar Active Regions from Major Flares

Authors: Yang Liu, Brian T. Welsch, Gherardo Valori, Manolis K. Georgoulis, Yang Guo, Etienne Pariat, Sung-Hong Park, Julia K. Thalmann

Abstract: Magnetic free energy powers solar flares and coronal mass ejections (CMEs), and the buildup of magnetic helicity might play a role in the development of unstable structures that subsequently erupt. To better understand the roles of energy and helicity in large flares and eruptions, we have characterized the evolution of magnetic energy and helicity associated with 21 X-class flares from 2010 to 20… ▽ More Magnetic free energy powers solar flares and coronal mass ejections (CMEs), and the buildup of magnetic helicity might play a role in the development of unstable structures that subsequently erupt. To better understand the roles of energy and helicity in large flares and eruptions, we have characterized the evolution of magnetic energy and helicity associated with 21 X-class flares from 2010 to 2017. Our sample includes both confined and eruptive events, with 6 and 15 in each category, respectively. Using HMI vector magnetic field observations from several hours before to several hours after each event, we employ (a) the Differential Affine Velocity Estimator for Vector Magnetograms (DAVE4VM) to determine the photospheric fluxes of energy and helicity, and (b) non-linear force-free field (NLFFF) extrapolations to estimate the coronal content of energy and helicity in source-region fields. Using Superposed Epoch analysis (SPE), we find, on average: (1) decreases in both magnetic energy and helicity, in both photospheric fluxes and coronal content, that persist for a few hours after eruptions, but no clear changes, notably in relative helicity, for confined events; (2) significant increases in the twist of photospheric fields in eruptive events, with twist uncertainties too large in confined events to constrain twist changes (and lower overall twist in confined events); and (3) on longer time scales (event time +12 hours), replenishment of free magnetic energy and helicity content to near pre-event levels for eruptive events. For eruptive events, magnetic helicity and free energy in coronal models clearly decrease after flares, with the amounts of decrease proportional to each region's pre-flare content. △ Less

Submitted 17 November, 2022; originally announced November 2022.

Comments: 25 pages, 15 figures, 1 table, accepted for publication in the Astrophysical Journal

The Photospheric Imprints of Coronal Electric Currents

Authors: Brian T. Welsch

Abstract: Flares and coronal mass ejections are powered by magnetic energy stored in coronal electric currents. Here, we explore the nature of coronal currents in observed and model active region (ARs) by studying manifestations of these currents in photospheric vector magnetograms. We employ Gauss's separation method, recently introduced to the solar physics literature, to partition the photospheric field… ▽ More Flares and coronal mass ejections are powered by magnetic energy stored in coronal electric currents. Here, we explore the nature of coronal currents in observed and model active region (ARs) by studying manifestations of these currents in photospheric vector magnetograms. We employ Gauss's separation method, recently introduced to the solar physics literature, to partition the photospheric field into three distinct components, each arising from a separate source: (i) currents passing through the photosphere, (ii) currents flowing below it, and (iii) currents flowing above it. We refer to component (iii) as the photospheric imprint of coronal currents. In both AR 10930 and AR 11158, photospheric imprints exhibit large-scale, spatially coherent structures along these regions' central, sheared polarity inversion lines (PILs) that are consistent with coronal currents flowing horizontally above these PILs, similar to recent findings in AR 12673 by Schuck et al. (2022). We find similar photospheric imprints in a simple model of a non-potential AR with known currents. We find that flare-associated changes in photospheric imprints in AR 11158 accord with earlier reports that near-PIL fields become more horizontal, consistent with the "implosion" scenario. We hypothesize that this evolution effectively shortens, in an overall sense, current-carrying coronal fields, leading to decreased inductive energy (DIE) in the coronal field. We further hypothesize that, in the hours prior to flares, parts of the coronal field slowly expand, in a process we deem coronal inflation (CI) -- essentially, the inverse of the implosion process. Both of these hypotheses are testable with non-potential coronal field extrapolations. △ Less

Submitted 8 November, 2022; v1 submitted 3 November, 2022; originally announced November 2022.

Comments: 28 pages, 10 figures, to be submitted to ApJ. Addition of co-authors is expected

Reconstruction of photospheric velocity fields from highly corrupted data

Authors: Erico L. Rempel, Roman Chertovskih, Kamilla R. Davletshina, Suzana S. A. Silva, Brian T. Welsch, Abraham C. -L. Chian

Abstract: The analysis of the photospheric velocity field is essential for understanding plasma turbulence in the solar surface, which may be responsible for driving processes such as magnetic reconnection, flares, wave propagation, particle acceleration, and coronal heating. Currently, the only available methods to estimate velocities at the solar photosphere transverse to an observer's line of sight infer… ▽ More The analysis of the photospheric velocity field is essential for understanding plasma turbulence in the solar surface, which may be responsible for driving processes such as magnetic reconnection, flares, wave propagation, particle acceleration, and coronal heating. Currently, the only available methods to estimate velocities at the solar photosphere transverse to an observer's line of sight infer flows from differences in image structure in successive observations. Due to data noise, algorithms such as local correlation tracking (LCT) may lead to a vector field with wide gaps where no velocity vectors are provided. In this letter, a novel method for image inpainting of highly corrupted data is proposed and applied to the restoration of horizontal velocity fields in the solar photosphere. The restored velocity field preserves all the vector field components present in the original field. The method shows robustness when applied to both simulated and observational data. △ Less

Submitted 19 May, 2022; originally announced May 2022.

Journal ref: ApJ 933 2 (2022)

Toward Improved Understanding of Magnetic Fields Participating in Solar Flares: Statistical Analysis of Magnetic Field within Flare Ribbons

Authors: Maria D. Kazachenko, Benjamin J. Lynch, Antonia Savcheva, Xudong Sun, Brian T. Welsch

Abstract: Violent solar flares and coronal mass ejections (CMEs) are magnetic phenomena. However, how magnetic fields reconnecting in the flare differ from non-flaring magnetic fields remains unclear owing to the lack of studies of the flare magnetic properties. Here we present a first statistical study of flaring (highlighted by flare-ribbons) vector magnetic fields in the photosphere. Our systematic appro… ▽ More Violent solar flares and coronal mass ejections (CMEs) are magnetic phenomena. However, how magnetic fields reconnecting in the flare differ from non-flaring magnetic fields remains unclear owing to the lack of studies of the flare magnetic properties. Here we present a first statistical study of flaring (highlighted by flare-ribbons) vector magnetic fields in the photosphere. Our systematic approach allows us to describe key physical properties of solar flare magnetism, including distributions of magnetic flux, magnetic shear, vertical current and net current over flaring versus non-flaring parts of the active region, and compare these with flare/CME properties. Our analysis suggests that while flares are guided by the physical properties that scale with AR size, like the total amount of magnetic flux that participates in the reconnection process and the total current (extensive properties), CMEs are guided by mean properties, like the fraction of the AR magnetic flux that participates (intensive property), with little dependence on the amount of shear at polarity inversion line (PIL) or the net current. We find that the non-neutralized current is proportional to the amount of shear at PIL, providing direct evidence that net vertical currents are formed as a result of any mechanism that could generate magnetic shear along PIL. We also find that eruptive events tend to have smaller PIL fluxes and larger magnetic shears than confined events. Our analysis provides a reference for more realistic solar and stellar flare models. The database is available online and can be used for future quantitative studies of flare magnetism. △ Less

Submitted 10 November, 2021; originally announced November 2021.

Comments: 25 pages, 12 figures, accepted to The Astrophysical Journal (10 November 2021); Link to FlareMagDB dataset: http://solarmuri.ssl.berkeley.edu/~kazachenko/FlareMagDB/

Oscillations observed in Umbra, Plage, Quiet-Sun and the Polarity Inversion Line of Active Region 11158 using HMI/SDO Data

Authors: A. A. Norton, R. B. Stutz, B. T. Welsch

Abstract: Using data from the Helioseismic Magnetic Imager, we report on the amplitudes and phase relations of oscillations in quiet-Sun, plage, umbra and the polarity inversion line (PIL) of an active region NOAA$\#$11158. We employ Fourier, wavelet and cross correlation spectra analysis. Waves with 5-minute periods are observed in umbra, PIL and plage with common phase values of $φ(v,I)=\fracπ{2}$,… ▽ More Using data from the Helioseismic Magnetic Imager, we report on the amplitudes and phase relations of oscillations in quiet-Sun, plage, umbra and the polarity inversion line (PIL) of an active region NOAA$\#$11158. We employ Fourier, wavelet and cross correlation spectra analysis. Waves with 5-minute periods are observed in umbra, PIL and plage with common phase values of $φ(v,I)=\fracπ{2}$, $φ(v,B_{los})=-\fracπ{2}$. In addition, $φ(I,B_{los})=π$ in plage are observed. These phase values are consistent with slow standing or fast standing surface sausage wave modes. The line width variations, and their phase relations with intensity and magnetic oscillations, show different values within the plage and PIL regions, which may offer a way to further differentiate wave mode mechanics. Significant Doppler velocity oscillations are present along the PIL, meaning that plasma motion is perpendicular to the magnetic field lines, a signature of Alvènic waves. A time-distance diagram along a section of the PIL shows Eastward propagating Doppler oscillations converting into magnetic oscillations; the propagation speeds range between 2$-$6 km s$^{-1}$. Lastly, a 3-minute wave is observed in select regions of the umbra in the magnetogram data. △ Less

Submitted 4 January, 2021; originally announced January 2021.

Comments: 20 pages, 12 figures

Journal ref: Philisophical Transactions of the Royal Society A, 2021, 379, 2190

Critical Science Plan for the Daniel K. Inouye Solar Telescope (DKIST)

Authors: Mark P. Rast, Nazaret Bello González, Luis Bellot Rubio, Wenda Cao, Gianna Cauzzi, Edward DeLuca, Bart De Pontieu, Lyndsay Fletcher, Sarah E. Gibson, Philip G. Judge, Yukio Katsukawa, Maria D. Kazachenko, Elena Khomenko, Enrico Landi, Valentin Martínez Pillet, Gordon J. D. Petrie, Jiong Qiu, Laurel A. Rachmeler, Matthias Rempel, Wolfgang Schmidt, Eamon Scullion, Xudong Sun, Brian T. Welsch, Vincenzo Andretta, Patrick Antolin , et al. (62 additional authors not shown)

Abstract: The Daniel K. Inouye Solar Telescope (DKIST) will revolutionize our ability to measure, understand and model the basic physical processes that control the structure and dynamics of the Sun and its atmosphere. The first-light DKIST images, released publicly on 29 January 2020, only hint at the extraordinary capabilities which will accompany full commissioning of the five facility instruments. With… ▽ More The Daniel K. Inouye Solar Telescope (DKIST) will revolutionize our ability to measure, understand and model the basic physical processes that control the structure and dynamics of the Sun and its atmosphere. The first-light DKIST images, released publicly on 29 January 2020, only hint at the extraordinary capabilities which will accompany full commissioning of the five facility instruments. With this Critical Science Plan (CSP) we attempt to anticipate some of what those capabilities will enable, providing a snapshot of some of the scientific pursuits that the Daniel K. Inouye Solar Telescope hopes to engage as start-of-operations nears. The work builds on the combined contributions of the DKIST Science Working Group (SWG) and CSP Community members, who generously shared their experiences, plans, knowledge and dreams. Discussion is primarily focused on those issues to which DKIST will uniquely contribute. △ Less

Submitted 20 August, 2020; v1 submitted 18 August, 2020; originally announced August 2020.

The Coronal Global Evolutionary Model: Using HMI Vector Magnetogram and Doppler Data to Determine Coronal Magnetic Field Evolution

Authors: J. Todd Hoeksema, William P. Abbett, David J. Bercik, Mark C. M. Cheung, Marc L. DeRosa, George H. Fisher, Keiji Hayashi, Maria D. Kazachenko, Yang Liu, Erkka Lumme, Benjamin J. Lynch, Xudong Sun, Brian T. Welsch

Abstract: The Coronal Global Evolutionary Model (CGEM) provides data-driven simulations of the magnetic field in the solar corona to better understand the build-up of magnetic energy that leads to eruptive events. The CGEM project has developed six capabilities. CGEM modules (1) prepare time series of full-disk vector magnetic field observations to (2) derive the changing electric field in the solar photosp… ▽ More The Coronal Global Evolutionary Model (CGEM) provides data-driven simulations of the magnetic field in the solar corona to better understand the build-up of magnetic energy that leads to eruptive events. The CGEM project has developed six capabilities. CGEM modules (1) prepare time series of full-disk vector magnetic field observations to (2) derive the changing electric field in the solar photosphere over active-region scales. This local electric field is (3) incorporated into a surface flux transport model that reconstructs a global electric field that evolves magnetic flux in a consistent way. These electric fields drive a (4) 3D spherical magneto-frictional (SMF) model, either at high-resolution over a restricted range of solid angle or at lower resolution over a global domain, to determine the magnetic field and current density in the low corona. An SMF-generated initial field above an active region and the evolving electric field at the photosphere are used to drive (5) detailed magneto-hydrodynamic (MHD) simulations of active regions in the low corona. SMF or MHD solutions are then used to compute emissivity proxies that can be compared with coronal observations. Finally, a lower-resolution SMF magnetic field is used to initialize (6) a global MHD model that is driven by an SMF electric-field time series to simulate the outer corona and heliosphere, ultimately connecting Sun to Earth. As a demonstration, this report features results of CGEM applied to observations of the evolution of NOAA Active Region 11158 in February 2011. △ Less

Submitted 2 October, 2020; v1 submitted 25 June, 2020; originally announced June 2020.

Comments: 19 pages, 8 figures

Journal ref: Astrophysical Journal Supplement Series 250:28 (15pp), 2020 October

The PDFI_SS Electric Field Inversion Software

Authors: George H. Fisher, Maria D. Kazachenko, Brian T. Welsch, Xudong Sun, Erkka Lumme, David J. Bercik, Marc L. DeRosa, Mark C. M. Cheung

Abstract: 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 Kazache… ▽ More 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. △ Less

Submitted 24 April, 2020; v1 submitted 17 December, 2019; originally announced December 2019.

Comments: 68 pages, 28 figures, 3 tables. This version (arXiv) will be updated when major changes to the software occur

Journal ref: Astrophysical Journal Supplement Series, 248:2 (52pp) 2020 May

Probing the effect of cadence on the estimates of photospheric energy and helicity injections in eruptive active region NOAA AR 11158

Authors: E. Lumme, M. D. Kazachenko, G. H. Fisher, B. T. Welsch, J. Pomoell, E. K. J. Kilpua

Abstract: In this work we study how the input data cadence affects the photospheric energy and helicity injection estimates in eruptive NOAA active region 11158. We sample the novel 2.25-minute vector magnetogram and Dopplergram data from the \emph{Helioseismic and Magnetic Imager} (HMI) instrument onboard the \emph{Solar Dynamics Observatory} (SDO) spacecraft to create input datasets of variable cadences r… ▽ More In this work we study how the input data cadence affects the photospheric energy and helicity injection estimates in eruptive NOAA active region 11158. We sample the novel 2.25-minute vector magnetogram and Dopplergram data from the \emph{Helioseismic and Magnetic Imager} (HMI) instrument onboard the \emph{Solar Dynamics Observatory} (SDO) spacecraft to create input datasets of variable cadences ranging from 2.25 minutes to 24 hours. We employ state-of-the-art data processing, velocity and electric field inversion methods for deriving estimates of the energy and helicity injections from these datasets. We find that the electric field inversion methods that reproduce the observed magnetic field evolution through the use of Faraday's law are more stable against variable cadence: the PDFI (PTD-Doppler-FLCT-Ideal) electric field inversion method produces consistent injection estimates for cadences from 2.25 minutes up to 2 hours, implying that the photospheric processes acting on time scales below 2 hours contribute little to the injections, or that they are below the sensitivity of the input data and the PDFI method. On other hand, the electric field estimate derived from the output of DAVE4VM (Differential Affine Velocity Estimator for Vector Magnetograms), which does not fulfil Faraday's law exactly, produces significant variations in the energy and helicity injection estimates in the 2.25-minute to 2-hour cadence range. We present also a third, novel DAVE4VM-based electric field estimate, which corrects the poor inductivity of the raw DAVE4VM estimate. This method is less sensitive to the changes of cadence, but still faces significant issues for the lowest of considered cadences ($\geq$2 hours). We find several potential problems in both PDFI- and DAVE4VM-based injection estimates and conclude that the quality of both should be surveyed further in controlled environments. △ Less

Submitted 30 June, 2019; originally announced July 2019.

Journal ref: Solar Physics, 294, 84 (2019)

A Database of Flare Ribbon Properties From Solar Dynamics Observatory I: Reconnection Flux

Authors: Maria D. Kazachenko, Benjamin J. Lynch, Brian T. Welsch, Xudong Sun

Abstract: We present a database of 3137 solar flare ribbon events corresponding to every flare of GOES class C1.0 and greater within 45 degrees from the central meridian, from April 2010 until April 2016, observed by the \emph{Solar Dynamics Observatory}. For every event in the database, we compare the GOES peak X-ray flux with corresponding active-region and flare-ribbon properties. We find that while the… ▽ More We present a database of 3137 solar flare ribbon events corresponding to every flare of GOES class C1.0 and greater within 45 degrees from the central meridian, from April 2010 until April 2016, observed by the \emph{Solar Dynamics Observatory}. For every event in the database, we compare the GOES peak X-ray flux with corresponding active-region and flare-ribbon properties. We find that while the peak X-ray flux is not correlated with the active region unsigned magnetic flux, it is strongly correlated with the flare ribbon reconnection flux, flare ribbon area, and the fraction of active region flux that undergoes reconnection. We find the relationship between the peak X-ray flux and the flare ribbon reconnection flux to be $I_\mathrm{X,peak} \propto Φ_\mathrm{ribbon}^{1.5}$. This scaling law is consistent with earlier hydrodynamic simulations of impulsively heated flare loops. Using the flare reconnection flux as a proxy for the total released flare energy $E$, we find that the occurrence frequency of flare energies follows a power-law dependence: $dN/dE \propto E^{-1.6}$ for $10^{31}<E<10^{33}$ erg, consistent with earlier studies of solar and stellar flares. The database is available online and can be used for future quantitative studies of flares. △ Less

Submitted 6 July, 2017; v1 submitted 17 April, 2017; originally announced April 2017.

Comments: 15 pages, 12 figures, Accepted by the Astrophysical Journal

Flux Accretion and Coronal Mass Ejection Dynamics

Authors: Brian T. Welsch

Abstract: Coronal mass ejections (CMEs) are the primary drivers of severe space weather disturbances in the heliosphere. Models of CME dynamics have been proposed that do not fully include the effects of magnetic reconnection on the forces driving the ejection. Both observations and numerical modeling, however, suggest that reconnection likely plays a major role in most, if not all, fast CMEs. Here, we theo… ▽ More Coronal mass ejections (CMEs) are the primary drivers of severe space weather disturbances in the heliosphere. Models of CME dynamics have been proposed that do not fully include the effects of magnetic reconnection on the forces driving the ejection. Both observations and numerical modeling, however, suggest that reconnection likely plays a major role in most, if not all, fast CMEs. Here, we theoretically investigate the accretion of magnetic flux onto a rising ejection by reconnection involving the ejection's background field. This reconnection alters the magnetic structure of the ejection and its environment, thereby modifying the forces acting upon the ejection, generically increasing its upward acceleration. The modified forces, in turn, can more strongly drive the reconnection. This feedback process acts, effectively, as an instability, which we refer to as a reconnective instability. Our analysis implies that CME models that neglect the effects of reconnection cannot accurately describe observed CME dynamics. Our ultimate aim is to understand changes in CME acceleration in terms of observable properties of magnetic reconnection, such as the amount of reconnected flux. This flux can be estimated from observations of flare ribbons and photospheric magnetic fields. △ Less

Submitted 19 July, 2018; v1 submitted 31 January, 2017; originally announced January 2017.

Comments: 29 pages, 12 figures; fifth revision, accepted for publication in Solar Physics

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

Authors: George H. Fisher, William. P. Abbett, David J. Bercik, Maria D. Kazachenko, Benjamin J. Lynch, Brian T. Welsch, J. Todd Hoeksema, Keiji Hayashi, Yang Liu, Aimee A. Norton, Alberto Sainz Dalda, Xudong Sun, Marc L. DeRosa, Mark C. M. Cheung

Abstract: 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 us… ▽ More 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. △ Less

Submitted 9 July, 2015; v1 submitted 22 May, 2015; originally announced May 2015.

Comments: 7 pages, 2 figures. Feature article in Space Weather 13

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

Authors: Maria D. Kazachenko, George H. Fisher, Brian T. Welsch, Yang Liu, Xudong Sun

Abstract: 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 s… ▽ More 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. △ Less

Submitted 11 August, 2015; v1 submitted 22 May, 2015; originally announced May 2015.

Comments: 51 pages, 10 figures, accepted by ApJ on August 11, 2015

Active Region Emergence & Remote Flares

Authors: Yixing Fu, Brian T. Welsch

Abstract: 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 difference between (i) the summed magnetic energies of their individual potential fields and (ii) the energy of their superposed potential fields. We expect… ▽ More 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 difference 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 NOAA's 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. △ Less

Submitted 8 June, 2016; v1 submitted 24 April, 2015; originally announced April 2015.

Comments: 23 pages, 9 figures, 2 tables; published in Solar Physics

Journal ref: Solar Physics, Volume 291, Issue 2, pp.383-410 (02/2016)

The Roles of Reconnected Flux and Overlying Fields in CME Speeds

Authors: Minda Deng, Brian T. Welsch

Abstract: 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 deca… ▽ More 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. △ Less

Submitted 7 June, 2016; v1 submitted 11 April, 2015; originally announced April 2015.

Comments: 12 pages, 6 figures, 2 tables; under review

Deriving Potential Coronal Magnetic Fields from Vector Magnetograms

Authors: Brian T. Welsch, George H. Fisher

Abstract: 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 photosp… ▽ More 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. △ Less

Submitted 16 June, 2016; v1 submitted 30 March, 2015; originally announced March 2015.

Comments: 25 pages, 11 figures, under review

Sign singularity and flares in solar active region NOAA 11158

Authors: Luca Sorriso-Valvo, Gaetano De Vita, Maria D. Kazachenko, Sam Krucker, Leonardo Primavera, Sergio Servidio, Antonio Vecchio, Brian T. Welsch, George H. Fisher, Fabio Lepreti, Vincenzo Carbone

Abstract: 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 o… ▽ More 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. △ Less

Submitted 18 January, 2015; originally announced January 2015.

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

Authors: Maria D. Kazachenko, George H. Fisher, Brian T. Welsch

Abstract: 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… ▽ More 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. △ Less

Submitted 10 September, 2014; v1 submitted 15 April, 2014; originally announced April 2014.

Comments: 56 pages, 10 figures, ApJ (in press)

Journal ref: 2014, The Astrophysical Journal, Volume 795, Number 1

The coronal energy input from magnetic braiding

Authors: A. R. Yeates, F. Bianchi, B. T. Welsch, P. J. Bushby

Abstract: We estimate the energy input into the solar corona from photospheric footpoint motions, using observations of a plage region by the Hinode Solar Optical Telescope. Assuming a perfectly ideal coronal evolution, two alternative lower bounds for the Poynting flux are computed based on field line footpoint trajectories, without requiring horizontal magnetic field data. When applied to the observed vel… ▽ More We estimate the energy input into the solar corona from photospheric footpoint motions, using observations of a plage region by the Hinode Solar Optical Telescope. Assuming a perfectly ideal coronal evolution, two alternative lower bounds for the Poynting flux are computed based on field line footpoint trajectories, without requiring horizontal magnetic field data. When applied to the observed velocities, a bound based solely on displacements between the two footpoints of each field line is tighter than a bound based on relative twist between field lines. Depending on the assumed length of coronal magnetic field lines, the higher bound is found to be reasonably tight compared with a Poynting flux estimate using an available vector magnetogram. It is also close to the energy input required to explain conductive and radiative losses in the active region corona. Based on similar analysis of a numerical convection simulation, we suggest that observations with higher spatial resolution are likely to bring the bound based on relative twist closer to the first bound, but not to increase the first bound substantially. Finally, we put an approximate upper bound on the magnetic energy by constructing a hypothetical ``unrelaxed'' magnetic field with the correct field line connectivity. △ Less

Submitted 18 March, 2014; originally announced March 2014.

Comments: 10 pages, 11 figures, accepted in Astronomy & Astrophysics

The Photospheric Poynting Flux and Coronal Heating

Authors: B. T. Welsch

Abstract: Some models of coronal heating suppose that convective motions at the photosphere shuffle the footpoints of coronal magnetic fields and thereby inject sufficient magnetic energy upward to account for observed coronal and chromospheric energy losses in active regions. Using high-resolution observations of plage magnetic fields made with the Solar Optical Telescope aboard the Hinode satellite, we in… ▽ More Some models of coronal heating suppose that convective motions at the photosphere shuffle the footpoints of coronal magnetic fields and thereby inject sufficient magnetic energy upward to account for observed coronal and chromospheric energy losses in active regions. Using high-resolution observations of plage magnetic fields made with the Solar Optical Telescope aboard the Hinode satellite, we investigate this idea by estimating the upward transport of magnetic energy --- the vertical Poynting flux, S_z --- across the photosphere in a plage region. To do so, we combine: (i) estimates of photospheric horizontal velocities, v_h, determined by local correlation tracking applied to a sequence of line-of-sight magnetic field maps from the Narrowband Filter Imager, with (ii) a vector magnetic field measurement from the SpectroPolarimeter. Plage fields are ideal observational targets for estimating energy injection by convection, because they are: (i) strong enough to be measured with relatively small uncertainties; (ii) not so strong that convection is heavily suppressed (as within umbrae); and (iii) unipolar, so S_z in plage is not influenced by mixed-polarity processes (e.g., flux emergence) unrelated to heating in stable, active-region fields. In this plage region, we found that the average S_z varied in space, but was positive (upward) and sufficient to explain coronal heating, with values near (5 +/- 1) x 10^7 erg/cm^2/s. We find the energy input per unit magnetic flux to be on the order of 10^5 erg/s/Mx. A comparison of intensity in a Ca II image co-registered with the this plage shows stronger spatial correlations with both total field, B, and unsigned vertical field, |B_z|, than either S_z or horizontal field, B_h. The observed Ca II brightness enhancement, however, probably contains a strong contribution from a near-photosphere hot-wall effect unrelated to atmospheric heating. △ Less

Submitted 4 December, 2014; v1 submitted 19 February, 2014; originally announced February 2014.

Comments: 30 pages, 11 figures, accepted by Pub. Astron. Soc. Japan

Detection of coherent structures in photospheric turbulent flows

Authors: Abraham C. -L. Chian, Erico L. Rempel, Guillaume Aulanier, Brigitte Schmieder, Shawn C. Shadden, Brian T. Welsch, Anthony R. Yeates

Abstract: We study coherent structures in solar photospheric flows in a plage in the vicinity of the active region AR 10930 using the horizontal velocity data derived from Hinode/SOT magnetograms. Eulerian and Lagrangian coherent structures are detected by computing the Q-criterion and the finite-time Lyapunov exponents of the velocity field, respectively. Our analysis indicates that, on average, the deform… ▽ More We study coherent structures in solar photospheric flows in a plage in the vicinity of the active region AR 10930 using the horizontal velocity data derived from Hinode/SOT magnetograms. Eulerian and Lagrangian coherent structures are detected by computing the Q-criterion and the finite-time Lyapunov exponents of the velocity field, respectively. Our analysis indicates that, on average, the deformation Eulerian coherent structures dominate over the vortical Eulerian coherent structures in the plage region. We demonstrate the correspondence of the network of high magnetic flux concentration to the attracting Lagrangian coherent structures (a-LCS) in the photospheric velocity based on both observations and numerical simulations. In addition, the computation of a-LCS provides a measure of the local rate of contraction/expansion of the flow. △ Less

Submitted 9 December, 2013; originally announced December 2013.

Comments: Submitted to ApJ

Solar Magnetic Tracking. IV. The Death of Magnetic Features

Authors: Derek A. Lamb, Timothy A. Howard, Craig E. DeForest, Clare E. Parnell, Brian T. Welsch

Abstract: The removal of magnetic flux from the quiet-sun photosphere is important for maintaining the statistical steady-state of the magnetic field there, for determining the magnetic flux budget of the Sun, and for estimating the rate of energy injected into the upper solar atmosphere. Magnetic feature death is a measurable proxy for the removal of detectable flux. We used the SWAMIS feature tracking cod… ▽ More The removal of magnetic flux from the quiet-sun photosphere is important for maintaining the statistical steady-state of the magnetic field there, for determining the magnetic flux budget of the Sun, and for estimating the rate of energy injected into the upper solar atmosphere. Magnetic feature death is a measurable proxy for the removal of detectable flux. We used the SWAMIS feature tracking code to understand how nearly 20000 detected magnetic features die in an hour-long sequence of Hinode/SOT/NFI magnetograms of a region of quiet Sun. Of the feature deaths that remove visible magnetic flux from the photosphere, the vast majority do so by a process that merely disperses the previously-detected flux so that it is too small and too weak to be detected. The behavior of the ensemble average of these dispersals is not consistent with a model of simple planar diffusion, suggesting that the dispersal is constrained by the evolving photospheric velocity field. We introduce the concept of the partial lifetime of magnetic features, and show that the partial lifetime due to Cancellation of magnetic flux, 22 h, is 3 times slower than previous measurements of the flux turnover time. This indicates that prior feature-based estimates of the flux replacement time may be too short, in contrast with the tendency for this quantity to decrease as resolution and instrumentation have improved. This suggests that dispersal of flux to smaller scales is more important for the replacement of magnetic fields in the quiet Sun than observed bipolar cancellation. We conclude that processes on spatial scales smaller than those visible to Hinode dominate the processes of flux emergence and cancellation, and therefore also the quantity of magnetic flux that threads the photosphere. △ Less

Submitted 15 July, 2013; originally announced July 2013.

Comments: Accepted by ApJ

Are Decaying Magnetic Fields Above Active Regions Related to Coronal Mass Ejection Onset?

Authors: J. Suzuki, B. T. Welsch, Y. Li

Abstract: Coronal mass ejections (CMEs) are powered by magnetic energy stored in electric currents in coronal magnetic fields, with the pre-CME field in balance between outward magnetic pressure of the proto-ejecta and inward magnetic tension from confining overlying fields. In studies of global, current-free coronal magnetic field models --- Potential-Field Source-Surface (PFSS) models --- it has been repo… ▽ More Coronal mass ejections (CMEs) are powered by magnetic energy stored in electric currents in coronal magnetic fields, with the pre-CME field in balance between outward magnetic pressure of the proto-ejecta and inward magnetic tension from confining overlying fields. In studies of global, current-free coronal magnetic field models --- Potential-Field Source-Surface (PFSS) models --- it has been reported that model field strengths above flare sites tend to be weaker in when CMEs occur than when eruptions fail to occur. This suggests that potential field models might usefully quantify magnetic confinement. An implication of this idea is that a decrease in model field strength overlying a possible eruption site should correspond to diminished confinement, implying an eruption is more likely. We have searched for such an effect by {\em post facto} investigation of the time evolution of model field strengths above a sample of 10 eruption sites, which included both slow and fast CMEs. In most events we study, we find no statistically significant evolution in either: (i) the rate of magnetic field decay with height; (ii) the strength of overlying magnetic fields near 50 Mm; (iii) or the ratio of fluxes at low and high altitudes (below 1.1$R_{\odot}$, and between 1.1--1.5$R_{\odot}$, respectively). Instead, we found that overlying field strengths and overlying flux tend to increase slightly, and their rates of decay with height become slightly more gradual, consistent with increased confinement. Since CMEs occur regardless of whether the parameters we use to quantify confinement are increasing or decreasing, either: (i) these parameters do not accurately characterize confinement in CME source regions; or (ii) systematic evolution in the large-scale magnetic environment of CME source regions is not, by itself, a necessary condition for CMEs to occur; or both. △ Less

Submitted 20 November, 2012; originally announced November 2012.

Journal ref: J. Suzuki et al. 2012 ApJ 758 22

Global Energetics of Thirty-Eight Large Solar Eruptive Events

Authors: A. G. Emslie, B. R. Dennis, A. Y. Shih, P. C. Chamberlin, R. A. Mewaldt, C. S. Moore, G. H. Share, A. Vourlidas, B. T. Welsch

Abstract: We have evaluated the energetics of 38 solar eruptive events observed by a variety of spacecraft instruments between February 2002 and December 2006, as accurately as the observations allow. The measured energetic components include: (1) the radiated energy in the GOES 1 - 8 A band; (2) the total energy radiated from the soft X-ray (SXR) emitting plasma; (3) the peak energy in the SXR-emitting pla… ▽ More We have evaluated the energetics of 38 solar eruptive events observed by a variety of spacecraft instruments between February 2002 and December 2006, as accurately as the observations allow. The measured energetic components include: (1) the radiated energy in the GOES 1 - 8 A band; (2) the total energy radiated from the soft X-ray (SXR) emitting plasma; (3) the peak energy in the SXR-emitting plasma; (4) the bolometric radiated energy over the full duration of the event; (5) the energy in flare-accelerated electrons above 20 keV and in flare-accelerated ions above 1 MeV; (6) the kinetic and potential energies of the coronal mass ejection (CME); (7) the energy in solar energetic particles (SEPs) observed in interplanetary space; and (8) the amount of free (nonpotential) magnetic energy estimated to be available in the pertinent active region. Major conclusions include: (1) the energy radiated by the SXR-emitting plasma exceeds, by about half an order of magnitude, the peak energy content of the thermal plasma that produces this radiation; (2) the energy content in flare-accelerated electrons and ions is sufficient to supply the bolometric energy radiated across all wavelengths throughout the event; (3) the energy contents of flare-accelerated electrons and ions are comparable; (4) the energy in SEPs is typically a few percent of the CME kinetic energy (measured in the rest frame of the solar wind); and (5) the available magnetic energy is sufficient to power the CME, the flare-accelerated particles, and the hot thermal plasma. △ Less

Submitted 12 September, 2012; originally announced September 2012.

A Magnetic Calibration of Photospheric Doppler Velocities

Authors: Brian T. Welsch, George H. Fisher, Xudong Sun

Abstract: 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 u… ▽ More 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. △ Less

Submitted 13 February, 2013; v1 submitted 11 January, 2012; originally announced January 2012.

Comments: 49 pages, 15 figures, accepted by ApJ

Decorrelation Times of Photospheric Fields and Flows

Authors: Brian T. Welsch, Kanya Kusano, Tetsuya T. Yamamoto, K. Muglach

Abstract: 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 Decemb… ▽ More 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 $Δt$ between magnetogram pairs tracked, spatial binning applied to the magnetograms, and windowing parameter $σ$ 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, $τ$. For $Δt > τ$, 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 $Δt$. △ Less

Submitted 18 January, 2012; v1 submitted 27 October, 2011; originally announced October 2011.

Comments: 48 pages, 20 figures, submitted to the Astrophysical Journal; full-resolution images in manuscript (8MB) at http://solarmuri.ssl.berkeley.edu/~welsch/public/manuscripts/flow_lifetimes_v2.pdf

Lagrangian coherent structures in photospheric flows and their implications for coronal magnetic structure

Authors: A. R. Yeates, G. Hornig, B. T. Welsch

Abstract: Aims. We show how the build-up of magnetic gradients in the Sun's corona may be inferred directly from photospheric velocity data. This enables computation of magnetic connectivity measures such as the squashing factor without recourse to magnetic field extrapolation. Methods.Assuming an ideal evolution in the corona, and an initially uniform magnetic field, the subsequent field line mapping is… ▽ More Aims. We show how the build-up of magnetic gradients in the Sun's corona may be inferred directly from photospheric velocity data. This enables computation of magnetic connectivity measures such as the squashing factor without recourse to magnetic field extrapolation. Methods.Assuming an ideal evolution in the corona, and an initially uniform magnetic field, the subsequent field line mapping is computed by integrating trajectories of the (time-dependent) horizontal photospheric velocity field. The method is applied to a 12 hour high-resolution sequence of photospheric flows derived from Hinode/SOT magnetograms. Results. We find the generation of a network of quasi-separatrix layers in the magnetic field, which correspond to Lagrangian coherent structures in the photospheric velocity. The visual pattern of these structures arises primarily from the diverging part of the photospheric flow, hiding the effect of the rotational flow component: this is demonstrated by a simple analytical model of photospheric convection. We separate the diverging and rotational components from the observed flow and show qualitative agreement with purely diverging and rotational models respectively. Increasing the flow speeds in the model suggests that our observational results are likely to give a lower bound for the rate at which magnetic gradients are built up by real photospheric flows. Finally, we construct a hypothetical magnetic field with the inferred topology, that can be used for future investigations of reconnection and energy release. △ Less

Submitted 16 January, 2012; v1 submitted 18 October, 2011; originally announced October 2011.

Comments: 9 pages, 10 figures, accepted for publication in A&A

Photospheric Magnetic Evolution in the WHI Active Regions

Authors: B. T. Welsch, S. Christe, J. M. McTiernan

Abstract: Sequences of line-of-sight (LOS) magnetograms recorded by the Michelson-Doppler Imager are used to quantitatively characterize photospheric magnetic structure and evolution in three active regions that rotated across the Sun's disk during the Whole Heliosphere Interval (WHI), in an attempt to relate the photospheric magnetic properties of these active regions to flares and coronal mass ejections (… ▽ More Sequences of line-of-sight (LOS) magnetograms recorded by the Michelson-Doppler Imager are used to quantitatively characterize photospheric magnetic structure and evolution in three active regions that rotated across the Sun's disk during the Whole Heliosphere Interval (WHI), in an attempt to relate the photospheric magnetic properties of these active regions to flares and coronal mass ejections (CMEs). Several approaches are used in our analysis, on scales ranging from whole active regions, to magnetic features, to supergranular scales, and, finally, to individual pixels. We calculated several parameterizations of magnetic structure and evolution that have previously been associated with flare and CME activity, including total unsigned magnetic flux, magnetic flux near polarity inversion lines, amount of cancelled flux, the "proxy Poynting flux," and helicity flux. To catalog flare events, we used flare lists derived from both GOES and RHESSI observations. By most such measures, AR 10988 should have been the most flare- and CME-productive active region, and AR 10989 the least. Observations, however, were not consistent with this expectation: ARs 10988 and 10989 produced similar numbers of flares, and AR 10989 also produced a few CMEs. These results highlight present limitations of statistics-based flare and CME forecasting tools that rely upon line-of-sight photospheric magnetic data alone. △ Less

Submitted 24 March, 2011; v1 submitted 11 March, 2011; originally announced March 2011.

Comments: 30 pages, 10 figures, 4 tables; accepted for publication by Solar Physics, in WHI topical issue

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

Authors: George H. Fisher, Brian T. Welsch, William P. Abbett

Abstract: 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… ▽ More 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. △ Less

Submitted 16 March, 2011; v1 submitted 21 January, 2011; originally announced January 2011.

Comments: Solar Physics, in press. 14 pages, 3 figures

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

Authors: George H. Fisher, David J. Bercik, Brian T. Welsch, Hugh S. Hudson

Abstract: 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 interio… ▽ More 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. △ Less

Submitted 4 December, 2011; v1 submitted 27 June, 2010; originally announced June 2010.

Comments: 23 pages, 6 figures. Accepted by Solar Physics, Nov. 11, 2011

Estimating Electric Fields from Vector Magnetogram Sequences

Authors: G. H. Fisher, B. T. Welsch, W. P. Abbett, D. J. Bercik

Abstract: Determining the electric field (E-field) distribution on the Sun's photosphere is essential for quantitative studies of how energy flows from the Sun's photosphere, through the corona, and into the heliosphere. This E-field also provides valuable input for data-driven models of the solar atmosphere and the Sun-Earth system. We show how Faraday's Law can be used with observed vector magnetogram t… ▽ More Determining the electric field (E-field) distribution on the Sun's photosphere is essential for quantitative studies of how energy flows from the Sun's photosphere, through the corona, and into the heliosphere. This E-field also provides valuable input for data-driven models of the solar atmosphere and the Sun-Earth system. We show how Faraday's Law can be used with observed vector magnetogram time series to estimate the photospheric E-field, an ill-posed inversion problem. Our method uses a "poloidal-toroidal decomposition" (PTD) of the time derivative of the vector magnetic field. The PTD solutions are not unique; the gradient of a scalar potential can be added to the PTD E-field without affecting consistency with Faraday's Law. We present an iterative technique to determine a potential function consistent with ideal MHD evolution; but this E-field is also not a unique solution to Faraday's Law. Finally, we explore a variational approach that minimizes an energy functional to determine a unique E-field, similar to Longcope's "Minimum Energy Fit". The PTD technique, the iterative technique, and the variational technique are used to estimate E-fields from a pair of synthetic vector magnetograms taken from an MHD simulation; and these E-fields are compared with the simulation's known electric fields. These three techniques are then applied to a pair of vector magnetograms of solar active region NOAA AR8210, to demonstrate the methods with real data. △ Less

Submitted 6 May, 2010; v1 submitted 24 December, 2009; originally announced December 2009.

Comments: 41 pages, 10 figures

Journal ref: Astrophys.J.715:242-259,2010

What is the relationship between photospheric flow fields and solar flares?

Authors: B. T. Welsch, Y. Li, P. W. Schuck, G. H. Fisher

Abstract: We estimated photospheric velocities by separately applying the Fourier Local Correlation Tracking (FLCT) and Differential Affine Velocity Estimator (DAVE) methods to 2708 co-registered pairs of SOHO/MDI magnetograms, with nominal 96-minute cadence and ~2" pixels, from 46 active regions (ARs) from 1996-1998 over the time interval t45 when each AR was within 45^o of disk center. For each magnetog… ▽ More We estimated photospheric velocities by separately applying the Fourier Local Correlation Tracking (FLCT) and Differential Affine Velocity Estimator (DAVE) methods to 2708 co-registered pairs of SOHO/MDI magnetograms, with nominal 96-minute cadence and ~2" pixels, from 46 active regions (ARs) from 1996-1998 over the time interval t45 when each AR was within 45^o of disk center. For each magnetogram pair, we computed the average estimated radial magnetic field, B; and each tracking method produced an independently estimated flow field, u. We then quantitatively characterized these magnetic and flow fields by computing several extensive and intensive properties of each; extensive properties scale with AR size, while intensive properties do not depend directly on AR size. Intensive flow properties included moments of speeds, horizontal divergences, and radial curls; extensive flow properties included sums of these properties over each AR, and a crude proxy for the ideal Poynting flux, the total |u| B^2. Several magnetic quantities were also computed, including: total unsigned flux; a measure of the amount of unsigned flux near strong-field polarity inversion lines, R; and the total B^2. Next, using correlation and discriminant analysis, we investigated the associations between these properties and flares from the GOES flare catalog, when averaged over both t45 and shorter time windows, of 6 and 24 hours. We found R and total |u| B^2 to be most strongly associated with flares; no intensive flow properties were strongly associated with flares. △ Less

Submitted 14 October, 2009; v1 submitted 5 May, 2009; originally announced May 2009.

Comments: 57 pages, 13 figures; revised content; added URL to manuscript with higher-quality images

Journal ref: Astrophys.J.705:821-843,2009

FLCT: A Fast, Efficient Method for Performing Local Correlation Tracking

Authors: G. H. Fisher, B. T. Welsch

Abstract: We describe the computational techniques employed in the recently updated Fourier local correlation tracking (FLCT) method. The FLCT code is then evaluated using a series of simple, 2D, known flow patterns that test its accuracy and characterize its errors. We describe the computational techniques employed in the recently updated Fourier local correlation tracking (FLCT) method. The FLCT code is then evaluated using a series of simple, 2D, known flow patterns that test its accuracy and characterize its errors. △ Less

Submitted 28 December, 2007; originally announced December 2007.

Comments: 9 pages, 3 figures

On the Origin of Strong-Field Polarity Inversion Lines

Authors: B. T. Welsch, Y. Li

Abstract: Several studies have correlated observations of impulsive solar activity -- flares and coronal mass ejections (CMEs) -- with the amount of magnetic flux near strong-field polarity inversion lines (PILs) in active regions' photospheric magnetic fields, as measured in line-of-sight (LOS) magnetograms. Practically, this empirical correlation holds promise as a space weather forecasting tool. Scient… ▽ More Several studies have correlated observations of impulsive solar activity -- flares and coronal mass ejections (CMEs) -- with the amount of magnetic flux near strong-field polarity inversion lines (PILs) in active regions' photospheric magnetic fields, as measured in line-of-sight (LOS) magnetograms. Practically, this empirical correlation holds promise as a space weather forecasting tool. Scientifically, however, the mechanisms that generate strong gradients in photospheric magnetic fields remain unknown. Hypotheses include: the (1) emergence of highly twisted or kinked flux ropes, which possess strong, opposite-polarity fields in close proximity; (2) emergence of new flux in close proximity to old flux; and (3) flux cancellation driven by photospheric flows acting fields that have already emerged. If such concentrations of flux near strong gradients are formed by emergence, then increases in unsigned flux near strong gradients should be correlated with increases in total unsigned magnetic flux -- a signature of emergence. Here, we analyze time series of MDI line-of-sight (LOS) magnetograms from several dozen active regions, and conclude that increases in unsigned flux near strong gradients tend to occur during emergence, though strong gradients can arise without flux emergence. We acknowledge support from NSF-ATM 04-51438. △ Less

Submitted 2 October, 2007; originally announced October 2007.

Comments: 8 pages, 3 b/w figures

Surface Flows From Magnetograms

Authors: B. T. Welsch, G. H. Fisher

Abstract: Estimates of velocities from time series of photospheric and/or chromospheric vector magnetograms can be used to determine fluxes of magnetic energy (the Poynting flux) and helicity across the magnetogram layer, and to provide time-dependent boundary conditions for data-driven simulations of the solar atmosphere above this layer. Velocity components perpendicular to the magnetic field are necess… ▽ More Estimates of velocities from time series of photospheric and/or chromospheric vector magnetograms can be used to determine fluxes of magnetic energy (the Poynting flux) and helicity across the magnetogram layer, and to provide time-dependent boundary conditions for data-driven simulations of the solar atmosphere above this layer. Velocity components perpendicular to the magnetic field are necessary both to compute these transport rates and to derive model boundary conditions. Here, we discuss some possible approaches to estimating perpendicular flows from magnetograms. Since Doppler shifts contain contributions from flows parallel to the magnetic field, perpendicular velocities are not generally recoverable from Doppler shifts alone. The induction equation's vertical component relates evolution in $B_z$ to the perpendicular flow field, but has a finite null space, meaning some ``null'' flows, e.g., motions along contours of normal field, do not affect $B_z$. Consequently, additional information is required to accurately specify the perpendicular flow field. Tracking methods, which analyze $\partial_t B_z$ in a neighborhood, have a long heritage, but other approaches have recently been developed. In a recent paper, several such techniques were tested using synthetic magnetograms from MHD simulations. Here, we use the same test data to characterize: 1) the ability of the induction equation's normal component, by itself, to estimate flows; and 2) a tracking method's ability to recover flow components that are perpendicular to $\mathbf{B}$ and parallel to contours of $B_z$. This work has been supported by NASA Heliophysics Theory grant NNG05G144G. △ Less

Submitted 2 October, 2007; originally announced October 2007.

Comments: 12 pages, 5 b/w figures

Solar Magnetic Tracking. I. Software Comparison and Recommended Practices

Authors: C. E. DeForest, H. J. Hagenaar, D. A. Lamb, C. E. Parnell, B. T. Welsch

Abstract: Feature tracking and recognition are increasingly common tools for data analysis, but are typically implemented on an ad-hoc basis by individual research groups, limiting the usefulness of derived results when selection effects and algorithmic differences are not controlled. Specific results that are affected include the solar magnetic turnover time, the distributions of sizes, strengths, and li… ▽ More Feature tracking and recognition are increasingly common tools for data analysis, but are typically implemented on an ad-hoc basis by individual research groups, limiting the usefulness of derived results when selection effects and algorithmic differences are not controlled. Specific results that are affected include the solar magnetic turnover time, the distributions of sizes, strengths, and lifetimes of magnetic features, and the physics of both small scale flux emergence and the small-scale dynamo. In this paper, we present the results of a detailed comparison between four tracking codes applied to a single set of data from SOHO/MDI, describe the interplay between desired tracking behavior and parameterization of tracking algorithms, and make recommendations for feature selection and tracking practice in future work. △ Less

Submitted 23 April, 2007; originally announced April 2007.

Comments: In press for Astrophys. J. 2007

Journal ref: Astrophys.J.666:576-587,2007