Source author record

Xudong Sun

Xudong Sun 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

30works
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

30 published item(s)

preprint2026arXiv

Convergence of Consensus-Based Particle Methods for Nonconvex Bi-Level Optimization

In this paper, we study a consensus-based optimization method for nonconvex bi-level optimization, where the objective is to minimize an upper-level function over the set of global minimizers of a lower-level problem. The proposed approach is derivative-free, and constructs its consensus point via smooth quantile selection combined with a Gibbs-type Laplace approximation. We establish convergence guarantees for both the associated \textit{mean-field} dynamics and its \textit{finite-particle} approximation. In particular, under suitable assumptions on smooth quantile localization, error bounds, and stability, we show that the mean-field law reaches any arbitrary prescribed Wasserstein neighborhood of the target bi-level solution with an explicit exponential rate up to the hitting time. Numerical experiments on a two-dimensional constrained problem and neural network training further support the theoretical results.

preprint2022arXiv

A New Magnetic Parameter of Active Regions Distinguishing Large Eruptive and Confined Solar Flares

With the aim of investigating how the magnetic field in solar active regions (ARs) controls flare activity, i.e., whether a confined or eruptive flare occurs, we analyze 106 flares of Geostationary Operational Environmental Satellite (GOES) class $\geq$M1.0 during 2010$-$2019. We calculate mean characteristic twist parameters $α$$_{FPIL}$ within the "flaring polarity inversion line" region and $α$$_\mathrm{HFED}$ within the area of high photospheric magnetic free energy density, which both provide measures of the nonpotentiality of AR core region. Magnetic twist is thought to be related to the driving force of electric current-driven instabilities, such as the helical kink instability. We also calculate total unsigned magnetic flux ($Φ$$_\mathrm{AR}$) of ARs producing the flare, which describes the strength of the background field confinement. By considering both the constraining effect of background magnetic fields and the magnetic non-potentiality of ARs, we propose a new parameter $α$/$Φ$$_\mathrm{AR}$ to measure the probability for a large flare to be associated with a coronal mass ejection (CME). We find that in about 90\% of eruptive flares, $α$$_\mathrm{FPIL}$/$Φ$$_\mathrm{AR}$ and $α$$_\mathrm{HFED}$/$Φ$$_\mathrm{AR}$ are beyond critical values (2.2$\times$$10^{-24}$ and 3.2$\times$$10^{-24}$ Mm$^{-1}$ Mx$^{-1}$), whereas they are less than critical values in $\sim$ 80\% of confined flares. This indicates that the new parameter $α$/$Φ$$_\mathrm{AR}$ is well able to distinguish eruptive flares from confined flares. Our investigation suggests that the relative measure of magnetic nonpotentiality within the AR core over the restriction of the background field largely controls the capability of ARs to produce eruptive flares.

preprint2022arXiv

Exploring the Solar Poles: The Last Great Frontier of the Sun

Despite investments in multiple space and ground-based solar observatories by the global community, the Sun's polar regions remain unchartered territory - the last great frontier for solar observations. Breaching this frontier is fundamental to understanding the solar cycle - the ultimate driver of short-to-long term solar activity that encompasses space weather and space climate. Magnetohydrodynamic dynamo models and empirically observed relationships have established that the polar field is the primary determinant of the future solar cycle amplitude. Models of solar surface evolution of tilted active regions indicate that the mid to high latitude surges of magnetic flux govern dynamics leading to the reversal and build-up of polar fields. Our theoretical understanding and numerical models of this high latitude magnetic field dynamics and plasma flows - that are a critical component of the sunspot cycle - lack precise observational constraints. This limitation compromises our ability to observe the enigmatic kilo Gauss polar flux patches and constrain the polar field distribution at high latitudes. The lack of these observations handicap our understanding of how high latitude magnetic fields power polar jets, plumes, and the fast solar wind that extend to the boundaries of the heliosphere and modulate solar open flux and cosmic ray flux within the solar system. Accurate observation of the Sun's polar regions, therefore, is the single most outstanding challenge that confronts Heliophysics. This paper argues the scientific case for novel out of ecliptic observations of the Sun's polar regions, in conjunction with existing, or future multi-vantage point heliospheric observatories. Such a mission concept can revolutionize the field of Heliophysics like no other mission concept has - with relevance that transcends spatial regimes from the solar interior to the heliosphere.

preprint2022arXiv

On the Coordinate System of Space-Weather HMI Active Region Patches (SHARPs): A Technical Note

We describe the coordinate systems of two streams of HMI active region vector data. A distinction is made between (a) the 2D grid on which the field vector is measured (or sampled), and (b) the 3D coordinate established at each grid point, in which the field vector is presented. The HMI data reduction can involve coordinate changes on both, with those performed on the former termed "remapping", the latter "vector transformation". Relevant pipeline procedures are described. Useful examples are given for data analysis.

preprint2022arXiv

On the nature of photospheric horizontal magnetic field increase in major solar flares

Rapid increase of horizontal magnetic field ($B_h$) around the flaring polarity inversion line is the most prominent photospheric field change during flares. It is considered to be caused by the contraction of flare loops, the details behind which is still not fully understood. Here we investigate the $B_h$-increase in 35 major flares using HMI high-cadence vector magnetograms. We find that $B_h$-increase is always accompanied by the increase of field inclination. It usually initiates near the flare ribbons, showing step-like change in between the ribbons. In particular, its evolution in early flare phase shows close spatio-temporal correlation to flare ribbons. We further find that $B_h$-increase tends to have similar intensity in confined and eruptive flares, but larger spatial-extent in eruptive flares in a statistical sense. Its intensity and timescale have inverse and positive correlations to the initial ribbon separations, respectively. The results altogether are well consistent with a recent proposed scenario which suggests that the reconnection-driven contraction of flare loops enhances photospheric $B_h$ according to the ideal induction equation, providing statistical evidence to the reconnection-driven origin for $B_h$-increase for the first time.

preprint2022arXiv

The CGEM Lorentz Force Data from HMI Vector Magnetograms

We describe a new data product from the CGEM (Coronal Global Evolutionary Model) collaboration that estimates the Lorentz force in active regions (ARs) based on HMI vector magnetogram patches. Following Fisher et al. (2012), we compute three components of the integrated Lorentz force over the outer solar atmosphere every 12 minutes throughout an AR's disk passage. These estimates, differenced during solar eruptive events, can provide valuable diagnostics on dynamic processes. We describe the pipeline modules, provide data retrieval examples, and document some systematic uncertainties that users should be aware of. Finally we document the formal uncertainty propagation procedures.

preprint2021arXiv

Benchmarking time series classification -- Functional data vs machine learning approaches

Time series classification problems have drawn increasing attention in the machine learning and statistical community. Closely related is the field of functional data analysis (FDA): it refers to the range of problems that deal with the analysis of data that is continuously indexed over some domain. While often employing different methods, both fields strive to answer similar questions, a common example being classification or regression problems with functional covariates. We study methods from functional data analysis, such as functional generalized additive models, as well as functionality to concatenate (functional-) feature extraction or basis representations with traditional machine learning algorithms like support vector machines or classification trees. In order to assess the methods and implementations, we run a benchmark on a wide variety of representative (time series) data sets, with in-depth analysis of empirical results, and strive to provide a reference ranking for which method(s) to use for non-expert practitioners. Additionally, we provide a software framework in R for functional data analysis for supervised learning, including machine learning and more linear approaches from statistics. This allows convenient access, and in connection with the machine-learning toolbox mlr, those methods can now also be tuned and benchmarked.

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

Electric Current Neutralization in Solar Active Regions and Its Relation to Eruptive Activity

It is well established that magnetic free energy associated with electric currents powers solar flares and coronal mass ejections (CMEs) from solar active regions (ARs). However, the conditions that determine whether an AR will produce an eruption are not well understood. Previous work suggests that the degree to which the driving electric currents, or the sum of all currents within a single magnetic polarity, are neutralized may serve as a good proxy for assessing the ability of ARs to produce eruptions. Here, we investigate the relationship between current neutralization and flare/CME production using a sample of 15 flare-active and 15 flare-quiet ARs. All flare-quiet and 4 flare-active ARs are also CME-quiet. We additionally test the relation of current neutralization to the degree of shear along polarity inversion lines (PILs) in an AR. We find that flare-productive ARs are more likely to exhibit non-neutralized currents, specifically those that also produce a CME. We find that flare/CME-active ARs also exhibit higher degrees of PIL shear than flare/CME-quiet ARs. We additionally observe that currents become more neutralized during magnetic flux emergence in flare-quiet ARs. Our investigation suggests that current neutralization in ARs is indicative of their eruptive potential.

preprint2020arXiv

Maximum Entropy-Regularized Multi-Goal Reinforcement Learning

In Multi-Goal Reinforcement Learning, an agent learns to achieve multiple goals with a goal-conditioned policy. During learning, the agent first collects the trajectories into a replay buffer, and later these trajectories are selected randomly for replay. However, the achieved goals in the replay buffer are often biased towards the behavior policies. From a Bayesian perspective, when there is no prior knowledge about the target goal distribution, the agent should learn uniformly from diverse achieved goals. Therefore, we first propose a novel multi-goal RL objective based on weighted entropy. This objective encourages the agent to maximize the expected return, as well as to achieve more diverse goals. Secondly, we developed a maximum entropy-based prioritization framework to optimize the proposed objective. For evaluation of this framework, we combine it with Deep Deterministic Policy Gradient, both with or without Hindsight Experience Replay. On a set of multi-goal robotic tasks of OpenAI Gym, we compare our method with other baselines and show promising improvements in both performance and sample-efficiency.

preprint2020arXiv

Non-Neutralized Electric Current of Active Regions Explained as a Projection Effect

Active regions (ARs) often possess an observed net electric current in a single magnetic polarity. We show that such "non-neutralized" currents can arise from a geometric projection effect when a twisted flux tube obliquely intersects the photosphere. To this end, we emulate surface maps of an emerging AR by sampling horizontal slices of a semi-torus flux tube at various heights. Although the tube has no net toroidal current, its poloidal current, when projected along the vertical direction, amounts to a significant non-neutralized component on the surface. If the tube emerges only partially as in realistic settings, the non-neutralized current will 1) develop as double ribbons near the sheared polarity inversion line, (2) positively correlate with the twist, and 3) reach its maximum before the magnetic flux. The projection effect may be important to the photospheric current distribution, in particular during the early stages of flux emergence.

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

Comparison of Coronal Extrapolation Methods for Cycle 24 Using HMI Data

Two extrapolation models of the solar coronal magnetic field are compared using magnetogram data from the SDO/HMI instrument. The two models, a horizontal current-current sheet-source surface (HCCSSS) model and a potential field-source surface (PFSS) model differ in their treatment of coronal currents. Each model has its own critical variable, respectively the radius of a cusp surface and a source surface, and it is found that adjusting these heights over the period studied allows better fit between the models and the solar open flux at 1 AU as calculated from the Interplanetary Magnetic Field (IMF). The HCCSSS model provides the better fit for the overall period from 2010 November to 2015 May as well as for two subsets of the period - the minimum/rising part of the solar cycle, and the recently-identified peak in the IMF from mid-2014 to mid-2015 just after solar maximum. It is found that a HCCSSS cusp surface height of 1.7 Rsun provides the best fit to the IMF for the overall period, while 1.7 & 1.9 Rsun give the best fits for the two subsets. The corresponding values for the PFSS source surface height are 2.1, 2.2 and 2.0 Rsun respectively. This means that the HCCSSS cusp surface rises as the solar cycle progresses while the PFSS source surface falls.

preprint2016arXiv

Evolution of Magnetic Field and Energy in A Major Eruptive Active Region Based on SDO/HMI Observation

We report the evolution of magnetic field and its energy in NOAA active region 11158 over 5 days based on a vector magnetogram series from the Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamic Observatory (SDO). Fast flux emergence and strong shearing motion led to a quadrupolar sunspot complex that produced several major eruptions, including the first X-class flare of Solar Cycle 24. Extrapolated non-linear force-free coronal fields show substantial electric current and free energy increase during early flux emergence near a low-lying sigmoidal filament with sheared kilogauss field in the filament channel. The computed magnetic free energy reaches a maximum of ~2.6e32 erg, about 50% of which is stored below 6 Mm. It decreases by ~0.3e32 erg within 1 hour of the X-class flare, which is likely an underestimation of the actual energy loss. During the flare, the photospheric field changed rapidly: horizontal field was enhanced by 28% in the core region, becoming more inclined and more parallel to the polarity inversion line. Such change is consistent with the conjectured coronal field "implosion", and is supported by the coronal loop retraction observed by the Atmospheric Imaging Assembly (AIA). The extrapolated field becomes more "compact" after the flare, with shorter loops in the core region, probably because of reconnection. The coronal field becomes slightly more sheared in the lowest layer, relaxes faster with height, and is overall less energetic.

preprint2015arXiv

Complex Flare Dynamics Initiated by a Filament-Filament Interaction

We report on an eruption involving a relatively rare filament-filament interaction on 2013 June 21, observed by SDO and STEREO-B. The two filaments were separated in height with a 'double-decker' configuration. The eruption of the lower filament began simultaneously with a descent of the upper filament resulting in a convergence and direct interaction of the two filaments. The interaction was accompanied by the heating of surrounding plasma and an apparent crossing of a loop-like structure through the upper filament. The subsequent coalescence of the filaments drove a bright front ahead of the erupting structures. The whole process was associated with a C3.0 flare followed immediately by an M2.9 flare. Shrinking loops and descending dark voids were observed during the M2.9 flare at different locations above a C-shaped flare arcade as part of the energy release, giving us a unique insight into the flare dynamics.

preprint2015arXiv

Evidence of thermal conduction suppression in a solar flaring loop by coronal seismology of slow-mode waves

Analysis of a longitudinal wave event observed by the Atmospheric Imaging Assembly (AIA) onboard the Solar Dynamics Observatory (SDO) is presented. A time sequence of 131 A images reveals that a C-class flare occurred at one footpoint of a large loop and triggered an intensity disturbance (enhancement) propagating along it. The spatial features and temporal evolution suggest that a fundamental standing slow-mode wave could be set up quickly after meeting of two initial disturbances from the opposite footpoints. The oscillations have a period of ~12 min and a decay time of ~9 min. The measured phase speed of 500$\pm$50 km/s matches the sound speed in the heated loop of ~10 MK, confirming that the observed waves are of slow mode. We derive the time-dependent temperature and electron density wave signals from six AIA extreme-ultraviolet (EUV) channels, and find that they are nearly in phase.The measured polytropic index from the temperature and density perturbations is 1.64$\pm$0.08 close to the adiabatic index of 5/3 for an ideal monatomic gas. The interpretation based on a 1D linear MHD model suggests that the thermal conductivity is suppressed by at least a factor of 3 in the hot flare loop at 9 MK and above. The viscosity coefficient is determined by coronal seismology from the observed wave when only considering the compressive viscosity dissipation. We find that to interpret the rapid wave damping, the classical compressive viscosity coefficient needs to be enhanced by a factor of 15 as the upper limit.

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

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

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

preprint2015arXiv

The Role of Erupting Sigmoid in Triggering a Flare with Parallel and Large-Scale Quasi-Circular Ribbons

In this paper, we present observations and analysis of an interesting sigmoid formation, eruption and the associated flare that occurred on 2014 April 18 using multi-wavelength data sets. We discuss the possible role of the sigmoid eruption in triggering the flare, which consists of two different set of ribbons: parallel ribbons as well as a large-scale quasi-circular ribbon. Several observational evidence and nonlinear force-free field extrapolation results show the existence of a large-scale fan-spine type magnetic configuration with a sigmoid lying under a section of the fan dome. The event can be explained with the following two phases. During the pre-flare phase, we observed the formation and appearance of sigmoid via tether-cutting reconnection between the two sets of sheared fields under the fan dome. The second, main flare phase, features the eruption of the sigmoid, the subsequent flare with parallel ribbons, and a quasi-circular ribbon. We propose the following multi-stage successive reconnections scenario for the main flare. First, tether-cutting reconnection is responsible for the formation and the eruption of the sigmoid structure. Second, the reconnection occurred in the wake of the erupting sigmoid produces the parallel flare ribbons on the both sides of the circular polarity inversion line. Third, the null-type reconnection higher in the corona, possibly triggered by the erupting sigmoid, leads to the formation of a large quasi-circular ribbon. For the first time we suggest a mechanism for this type of flare consisting of a double set of ribbons triggered by an erupting sigmoid in a large scale fan-spine type magnetic configuration.

preprint2015arXiv

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

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

preprint2014arXiv

On Polar Magnetic Field Reversal and Surface Flux Transport During Solar Cycle 24

As each solar cycle progresses, remnant magnetic flux from active regions (ARs) migrates poleward to cancel the old-cycle polar field. We describe this polarity reversal process during Cycle 24 using four years (2010.33--2014.33) of line-of-sight magnetic field measurements from the Helioseismic and Magnetic Imager. The total flux associated with ARs reached maximum in the north in 2011, more than two years earlier than the south; the maximum is significantly weaker than Cycle 23. The process of polar field reversal is relatively slow, north-south asymmetric, and episodic. We estimate that the global axial dipole changed sign in October 2013; the northern and southern polar fields (mean above 60$^\circ$ latitude) reversed in November 2012 and March 2014, respectively, about 16 months apart. Notably, the poleward surges of flux in each hemisphere alternated in polarity, giving rise to multiple reversals in the north. We show that the surges of the trailing sunspot polarity tend to correspond to normal mean AR tilt, higher total AR flux, or slower mid-latitude near-surface meridional flow, while exceptions occur during low magnetic activity. In particular, the AR flux and the mid-latitude poleward flow speed exhibit a clear anti-correlation. We discuss how these features can be explained in a surface flux transport process that includes a field-dependent converging flow toward the ARs, a characteristic that may contribute to solar cycle variability.

preprint2014arXiv

The Helioseismic and Magnetic Imager (HMI) Vector Magnetic Field Pipeline: Overview and Performance

The Helioseismic and Magnetic Imager (HMI) began near-continuous full-disk solar measurements on 1 May 2010 from the Solar Dynamics Observatory (SDO). An automated processing pipeline keeps pace with observations to produce observable quantities, including the photospheric vector magnetic field, from sequences of filtergrams. The primary 720s observables were released in mid 2010, including Stokes polarization parameters measured at six wavelengths as well as intensity, Doppler velocity, and the line-of-sight magnetic field. More advanced products, including the full vector magnetic field, are now available. Automatically identified HMI Active Region Patches (HARPs) track the location and shape of magnetic regions throughout their lifetime. The vector field is computed using the Very Fast Inversion of the Stokes Vector (VFISV) code optimized for the HMI pipeline; the remaining 180 degree azimuth ambiguity is resolved with the Minimum Energy (ME0) code. The Milne-Eddington inversion is performed on all full-disk HMI observations. The disambiguation, until recently run only on HARP regions, is now implemented for the full disk. Vector and scalar quantities in the patches are used to derive active region indices potentially useful for forecasting; the data maps and indices are collected in the SHARP data series, hmi.sharp_720s. Patches are provided in both CCD and heliographic coordinates. HMI provides continuous coverage of the vector field, but has modest spatial, spectral, and temporal resolution. Coupled with limitations of the analysis and interpretation techniques, effects of the orbital velocity, and instrument performance, the resulting measurements have a certain dynamic range and sensitivity and are subject to systematic errors and uncertainties that are characterized in this report.

preprint2014arXiv

The Helioseismic and Magnetic Imager (HMI) Vector Magnetic Field Pipeline: SHARPs -- Space-weather HMI Active Region Patches

A new data product from the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO) called Space-weather HMI Active Region Patches (SHARPs) is now available. SDO/HMI is the first space-based instrument to map the full-disk photospheric vector magnetic field with high cadence and continuity. The SHARP data series provide maps in patches that encompass automatically tracked magnetic concentrations for their entire lifetime; map quantities include the photospheric vector magnetic field and its uncertainty, along with Doppler velocity, continuum intensity, and line-of-sight magnetic field. Furthermore, keywords in the SHARP data series provide several parameters that concisely characterize the magnetic-field distribution and its deviation from a potential-field configuration. These indices may be useful for active-region event forecasting and for identifying regions of interest. The indices are calculated per patch and are available on a twelve-minute cadence. Quick-look data are available within approximately three hours of observation; definitive science products are produced approximately five weeks later. SHARP data are available at http://jsoc.stanford.edu and maps are available in either of two different coordinate systems. This article describes the SHARP data products and presents examples of SHARP data and parameters.

preprint2014arXiv

The Magnetic Field of Active Region 11158 During the 2011 February 12-17 Flares : Differences between Photospheric Extrapolation and Coronal Forward-Fitting Methods

We developed a {\sl coronal non-linear force-free field (COR-NLFFF)} forward-fitting code that fits an approximate {\sl non-linear force-free field (NLFFF)} solution to the observed geometry of automatically traced coronal loops. In contrast to photospheric NLFFF codes, which calculate a magnetic field solution from the constraints of the transverse photospheric field, this new code uses coronal constraints instead, and this way provides important information on systematic errors of each magnetic field calculation method, as well as on the non-forcefreeness in the lower chromosphere. In this study we applied the COR-NLFFF code to active region NOAA 11158, during the time interval of 2011 Feb 12 to 17, which includes an X2.2 GOES-class flare plus 35 M and C-class flares. We calcuated the free magnetic energy with a 6-minute cadence over 5 days. We find good agreement between the two types of codes for the total nonpotential $E_N$ and potential energy $E_P$, but find up to a factor of 4 discrepancy in the free energy $E_{free}=E_N-E_P$, and up to a factor of 10 discrepancy in the decrease of the free energy $ΔE_{free}$ during flares. The coronal NLFFF code exhibits a larger time variability, and yields a decrease of free energy during the flare that is sufficient to satisfy the flare energy budget, while the photospheric NLFFF code shows much less time variability and an order of magnitude less free energy decrease during flares. The discrepancy may partly be due to the pre-processing of photospheric vector data, but more likely due to the non-forcefreeness in the lower chromosphere. We conclude that the coronal field cannot be correctly calculated based on photospheric data alone, but requires additional information on coronal loop geometries.

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

Hot Spine Loops and the Nature of a Late-Phase Solar Flare

The fan-spine magnetic topology is believed to be responsible for many curious features in solar explosive events. A spine field line links distinct flux domains, but direct observation of such feature has been rare. Here we report a unique event observed by the Solar Dynamic Observatory where a set of hot coronal loops (over 10 MK) connected to a quasi-circular chromospheric ribbon at one end and a remote brightening at the other. Magnetic field extrapolation suggests these loops are partly tracer of the evolving spine field line. Continuous slipping- and null-point-type reconnections were likely at work, energizing the loop plasma and transferring magnetic flux within and across the fan quasi-separatrix layer. We argue that the initial reconnection is of the "breakout" type, which then transitioned to a more violent flare reconnection with an eruption from the fan dome. Significant magnetic field changes are expected and indeed ensued. This event also features an extreme-ultraviolet (EUV) late phase, i.e. a delayed secondary emission peak in warm EUV lines (about 2-7 MK). We show that this peak comes from the cooling of large post-reconnection loops beside and above the compact fan, a direct product of eruption in such topological settings. The long cooling time of the large arcades contributes to the long delay; additional heating may also be required. Our result demonstrates the critical nature of cross-scale magnetic coupling - topological change in a sub-system may lead to explosions on a much larger scale.

preprint2013arXiv

Reflection and refraction of (magneto-)acoustic waves at the magnetic canopy: further evidences from multi-height seismic data

We recently presented evidences (Rajaguru et al. 2012) that seismic halos around expanding magnetic structures in the lower solar atmosphere are related to the acoustic to magnetoacoustic wave conversions, using multi-height data from Helioseismic and Magnetic Imager (HMI) and Atmospheric Imaging Assembly (AIA) (1700 and 1600 Åchannels) onboard Solar Dynamics Observatory (SDO). Using the same data, we here present and discuss further evidences through analyses of maps of phase-shifts between observables from different heights and their correspondence with oscillation power. The phase shift maps provide more direct signatures of reflection and refraction of (magneto-)acoustic wave modes.

preprint2012arXiv

A Non-radial Eruption in a Quadrupolar Magnetic Configuration with a Coronal Null

We report one of several homologous non-radial eruptions from NOAA active region (AR) 11158 that are strongly modulated by the local magnetic field as observed with the Solar Dynamic Observatory (SDO). A small bipole emerged in the sunspot complex and subsequently created a quadrupolar flux system. Non-linear force-free field (NLFFF) extrapolation from vector magnetograms reveals its energetic nature: the fast-shearing bipole accumulated ~2e31 erg free energy (10% of AR total) over just one day despite its relatively small magnetic flux (5% of AR total). During the eruption, the ejected plasma followed a highly inclined trajectory, over 60 degrees with respect to the radial direction, forming a jet-like, inverted-Y shaped structure in its wake. Field extrapolation suggests complicated magnetic connectivity with a coronal null point, which is favorable of reconnection between different flux components in the quadrupolar system. Indeed, multiple pairs of flare ribbons brightened simultaneously, and coronal reconnection signatures appeared near the inferred null. Part of the magnetic setting resembles that of a blowout-type jet; the observed inverted-Y structure likely outlines the open field lines along the separatrix surface. Owing to the asymmetrical photospheric flux distribution, the confining magnetic pressure decreases much faster horizontally than upward. This special field geometry likely guided the non-radial eruption during its initial stage.

preprint2012arXiv

Emerging dimmings of active regions observed by SDO

With the observations from the Atmospheric Imaging Assembly and the Helioseismic and Magnetic Imager aboard the \emph{Solar Dynamics Observatory}, we statistically investigate the emerging dimmings (EDs) of 24 isolated active regions (IARs) from June 2010 to May 2011. All the IARs present EDs in lower temperature lines (e.g., 171 Å) at their early emerging stages, meanwhile in higher temperature lines (e.g., 211 Å), the ED regions brighten continuously. There are two type of EDs: fan-shaped and halo-shaped. There are 19 fan-shaped EDs and 5 halo-shaped ones. The EDs appear several to more than ten hours delay to the first emergence of the IARs. The shortest delay is 3.6 hr and the longest 19.0 hr. The EDs last from 3.3 hr to 14.2 hr, with a mean duration of 8.3 hr. Before the appearance of the EDs, the emergence rate of the magnetic flux of the IARs is from 1.2 $\times$ 10$^{19}$ Mx hr$^{-1}$ to 1.4 $\times$ 10$^{20}$ Mx hr$^{-1}$. The larger the emergence rate is, the shorter the delay time is. While the dimmings appear, the magnetic flux of the IARs ranges from 8.8 $\times$ 10$^{19}$ Mx to 1.3 $\times$ 10$^{21}$ Mx. These observations imply that the reconfiguration of the coronal magnetic fields due to reconnection between the newly-emerging flux and the surrounding existing fields results in a new thermal distribution which leads to a dimming for the cooler channel (171 Å) and brightnening in the warmer channels.

preprint2012arXiv

Properties of high-frequency wave power halos around active regions: an analysis of multi-height data from HMI and AIA onboard SDO

We study properties of waves of frequencies above the photospheric acoustic cut-off of $\approx$5.3 mHz, around four active regions, through spatial maps of their power estimated using data from Helioseismic and Magnetic Imager (HMI) and Atmospheric Imaging Assembly (AIA) onboard Solar Dynamics Observatory (SDO). The wavelength channels 1600 Å and 1700 Å from AIA are now known to capture clear oscillation signals due to helioseismic p modes as well as waves propagating up through to the chromosphere. Here we study in detail, in comparison with HMI Doppler data, properties of the power maps, especially the so called 'acoustic halos' seen around active regions, as a function of wave frequencies, inclination and strength of magnetic field (derived from the vector field observations by HMI) and observation height. We infer possible signatures of (magneto-)acoustic wave refraction from the observation height dependent changes, and hence due to changing magnetic strength and geometry, in the dependences of power maps on the photospheric magnetic quantities. We discuss the implications for theories of p mode absorption and mode conversions by the magnetic field.