Source author record

Brian T. Welsch

Brian T. Welsch 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

16works
5topics
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

16 published item(s)

preprint2025arXiv

Multi-height probing of horizontal flows in the solar photosphere

We tested whether simultaneous spectropolarimetric imaging in two magnetically sensitive optical spectral lines, which probe two different layers of the solar atmosphere (the photosphere and the temperature minimum), can help constrain the depth variation of horizontal flows. We first tested the feasibility of our method using Fourier local correlation tracking (FLCT) to track physical quantities at different optical depths ($\logτ_{500}={-1,-2,-3,-4}$) in an atmosphere simulated with the MURaM code. We then inferred the horizontal distribution of the LOS magnetic field component from synthetic spectropolarimetric observations of Fe I 525.0 nm and Mg I b2 spectral lines, applied FLCT to the time sequence of these synthetic magnetograms, and compared our findings with the original height-dependent horizontal velocities. Tracking the LOS magnetic field component (which coincides with the vertical component at the disk center) yields horizontal velocities that, after appropriate temporal and spatial averaging, agree excellently with the horizontal component of the simulated velocities, both calculated at constant $τ_{500}$ surfaces, up to the temperature minimum ($\logτ_{500}=-3$). When tracking the temperature at constant $τ_{500}$ surfaces, this agreement already breaks down completely at the mid photosphere ($\logτ_{500}=-2$). Tracking the vertical component of the magnetic field inferred from synthetic observations of the Fe I 525.0 nm and the Mg I b2 spectral lines yields a satisfactory inference of the horizontal velocities in the mid-photosphere ($\logτ_{500}\approx-1$) and the temperature minimum ($\logτ_{500}\approx-3$), respectively. Our results indicate that high-spatial-resolution spectropolarimetric imaging in solar spectral lines can provide meaningful information about the horizontal plasma velocities over a range of heights.

preprint2022arXiv

Reconstruction of photospheric velocity fields from highly corrupted data

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.

preprint2021arXiv

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

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.

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

Active Region Emergence & Remote Flares

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.

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

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.

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.

preprint2013arXiv

Detection of coherent structures in photospheric turbulent flows

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.

preprint2013arXiv

Solar Magnetic Tracking. IV. The Death of Magnetic Features

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.

preprint2012arXiv

Decorrelation Times of Photospheric Fields and Flows

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

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.