Solar Magnetic Feature Detection and Tracking for Space Weather Monitoring

We present an automated system for detecting, tracking, and cataloging emerging active regions throughout their evolution and decay using SOHO Michelson Doppler Interferometer (MDI) magnetograms. The SolarMonitor Active Region Tracking (SMART) algorithm relies on consecutive image differencing to remove both quiet-Sun and transient magnetic features, and region-growing techniques to group flux concentrations into classifiable features. We determine magnetic properties such as region size, total flux, flux imbalance, flux emergence rate, Schrijver's R-value, R* (a modified version of R), and Falconer's measurement of non-potentiality. A persistence algorithm is used to associate developed active regions with emerging flux regions in previous measurements, and to track regions beyond the limb through multiple solar rotations. We find that the total number and area of magnetic regions on disk vary with the sunspot cycle. While sunspot numbers are a proxy to the solar magnetic field, SMART offers a direct diagnostic of the surface magnetic field and its variation over timescale of hours to years. SMART will form the basis of the active region extraction and tracking algorithm for the Heliophysics Integrated Observatory (HELIO)

A Nonlinear Force-Free Magnetic Field Approximation Suitable for Fast Forward-Fitting to Coronal Loops. I. Theory

We derive an analytical approximation of nonlinear force-free magnetic field solutions (NLFFF) that can efficiently be used for fast forward-fitting to solar magnetic data, constrained either by observed line-of-sight magnetograms and stereoscopically triangulated coronal loops, or by 3D vector-magnetograph data. The derived NLFFF solutions provide the magnetic field components $B_x({\bf x})$, $B_y({\bf x})$, $B_z({\bf x})$, the force-free parameter $\alpha({\bf x})$, the electric current density ${\bf j}({\bf x})$, and are accurate to second-order (of the nonlinear force-free $\alpha$-parameter). The explicit expressions of a force-free field can easily be applied to modeling or forward-fitting of many coronal phenomena.Comment: Solar Physics (in press), 26 pages, 11 figure

Random matrix ensembles of time correlation matrices to analyze visual lifelogs

Visual lifelogging is the process of automatically recording images and other sensor data for the purpose of aiding memory recall. Such lifelogs are usually created using wearable cameras. Given the vast amount of images that are maintained in a visual lifelog, it is a significant challenge for users to deconstruct a sizeable collection of images into meaningful events. In this paper, random matrix theory (RMT) is applied to a cross-correlation matrix C, constructed using SenseCam lifelog data streams to identify such events. The analysis reveals a number of eigenvalues that deviate from the spectrum suggested by RMT. The components of the deviating eigenvectors are found to correspond to “distinct significant events” in the visual lifelogs. Finally, the cross-correlation matrix C is cleaned by separating the noisy part from the non-noisy part. Overall, the RMT technique is shown to be useful to detect major events in SenseCam images

The Evolution of Sunspot Magnetic Fields Associated with a Solar Flare

Solar flares occur due to the sudden release of energy stored in active-region magnetic fields. To date, the pre-cursors to flaring are still not fully understood, although there is evidence that flaring is related to changes in the topology or complexity of an active region's magnetic field. Here, the evolution of the magnetic field in active region NOAA 10953 was examined using Hinode/SOT-SP data, over a period of 12 hours leading up to and after a GOES B1.0 flare. A number of magnetic-field properties and low-order aspects of magnetic-field topology were extracted from two flux regions that exhibited increased Ca II H emission during the flare. Pre-flare increases in vertical field strength, vertical current density, and inclination angle of ~ 8degrees towards the vertical were observed in flux elements surrounding the primary sunspot. The vertical field strength and current density subsequently decreased in the post-flare state, with the inclination becoming more horizontal by ~7degrees. This behaviour of the field vector may provide a physical basis for future flare forecasting efforts.Comment: Accepted for Publication in Solar Physics. 16 pages, 4 figure

Automated Detection of Coronal Loops using a Wavelet Transform Modulus Maxima Method

We propose and test a wavelet transform modulus maxima method for the au- tomated detection and extraction of coronal loops in extreme ultraviolet images of the solar corona. This method decomposes an image into a number of size scales and tracks enhanced power along each ridge corresponding to a coronal loop at each scale. We compare the results across scales and suggest the optimum set of parameters to maximise completeness while minimising detection of noise. For a test coronal image, we compare the global statistics (e.g., number of loops at each length) to previous automated coronal-loop detection algorithms

CHD4 and the NuRD complex directly control cardiac sarcomere formation

Cardiac development relies on proper cardiomyocyte differentiation, including expression and assembly of cell-type-specific actomyosin subunits into a functional cardiac sarcomere. Control of this process involves not only promoting expression of cardiac sarcomere subunits but also repressing expression of noncardiac myofibril paralogs. This level of transcriptional control requires broadly expressed multiprotein machines that modify and remodel the chromatin landscape to restrict transcription machinery access. Prominent among these is the nucleosome remodeling and deacetylase (NuRD) complex, which includes the catalytic core subunit CHD4. Here, we demonstrate that direct CHD4-mediated repression of skeletal and smooth muscle myofibril isoforms is required for normal cardiac sarcomere formation, function, and embryonic survival early in gestation. Through transcriptomic and genome-wide analyses of CHD4 localization, we identified unique CHD4 binding sites in smooth muscle myosin heavy chain, fast skeletal α-actin, and the fast skeletal troponin complex genes. We further demonstrate that in the absence of CHD4, cardiomyocytes in the developing heart form a hybrid muscle cell that contains cardiac, skeletal, and smooth muscle myofibril components. These misexpressed paralogs intercalate into the nascent cardiac sarcomere to disrupt sarcomere formation and cause impaired cardiac function in utero. These results demonstrate the genomic and physiological requirements for CHD4 in mammalian cardiac development

Multiresolution analysis of active region magnetic structure and its correlation with the Mt. Wilson classification and flaring activity

Two different multi-resolution analyses are used to decompose the structure of active region magnetic flux into concentrations of different size scales. Lines separating these opposite polarity regions of flux at each size scale are found. These lines are used as a mask on a map of the magnetic field gradient to sample the local gradient between opposite polarity regions of given scale sizes. It is shown that the maximum, average and standard deviation of the magnetic flux gradient for alpha, beta, beta-gamma and beta-gamma-delta active regions increase in the order listed, and that the order is maintained over all length-scales. This study demonstrates that, on average, the Mt. Wilson classification encodes the notion of activity over all length-scales in the active region, and not just those length-scales at which the strongest flux gradients are found. Further, it is also shown that the average gradients in the field, and the average length-scale at which they occur, also increase in the same order. Finally, there are significant differences in the gradient distribution, between flaring and non-flaring active regions, which are maintained over all length-scales. It is also shown that the average gradient content of active regions that have large flares (GOES class 'M' and above) is larger than that for active regions containing flares of all flare sizes; this difference is also maintained at all length-scales.Comment: Accepted for publication in Solar Physic

How to optimize nonlinear force-free coronal magnetic field extrapolations from SDO/HMI vector magnetograms?

The SDO/HMI instruments provide photospheric vector magnetograms with a high spatial and temporal resolution. Our intention is to model the coronal magnetic field above active regions with the help of a nonlinear force-free extrapolation code. Our code is based on an optimization principle and has been tested extensively with semi-analytic and numeric equilibria and been applied before to vector magnetograms from Hinode and ground based observations. Recently we implemented a new version which takes measurement errors in photospheric vector magnetograms into account. Photospheric field measurements are often due to measurement errors and finite nonmagnetic forces inconsistent as a boundary for a force-free field in the corona. In order to deal with these uncertainties, we developed two improvements: 1.) Preprocessing of the surface measurements in order to make them compatible with a force-free field 2.) The new code keeps a balance between the force-free constraint and deviation from the photospheric field measurements. Both methods contain free parameters, which have to be optimized for use with data from SDO/HMI. Within this work we describe the corresponding analysis method and evaluate the force-free equilibria by means of how well force-freeness and solenoidal conditions are fulfilled, the angle between magnetic field and electric current and by comparing projections of magnetic field lines with coronal images from SDO/AIA. We also compute the available free magnetic energy and discuss the potential influence of control parameters.Comment: 17 Pages, 6 Figures, Sol. Phys., accepte

Are Solar Active Regions with Major Flares More Fractal, Multifractal, or Turbulent than Others?

Multiple recent investigations of solar magnetic field measurements have raised claims that the scale-free (fractal) or multiscale (multifractal) parameters inferred from the studied magnetograms may help assess the eruptive potential of solar active regions, or may even help predict major flaring activity stemming from these regions. We investigate these claims here, by testing three widely used scale-free and multiscale parameters, namely, the fractal dimension, the multifractal structure function and its inertial-range exponent, and the turbulent power spectrum and its power-law index, on a comprehensive data set of 370 timeseries of active-region magnetograms (17,733 magnetograms in total) observed by SOHO's Michelson Doppler Imager (MDI) over the entire Solar Cycle 23. We find that both flaring and non-flaring active regions exhibit significant fractality, multifractality, and non-Kolmogorov turbulence but none of the three tested parameters manages to distinguish active regions with major flares from flare-quiet ones. We also find that the multiscale parameters, but not the scale-free fractal dimension, depend sensitively on the spatial resolution and perhaps the observational characteristics of the studied magnetograms. Extending previous works, we attribute the flare-forecasting inability of fractal and multifractal parameters to i) a widespread multiscale complexity caused by a possible underlying self-organization in turbulent solar magnetic structures, flaring and non-flaring alike, and ii) a lack of correlation between the fractal properties of the photosphere and overlying layers, where solar eruptions occur. However useful for understanding solar magnetism, therefore, scale-free and multiscale measures may not be optimal tools for active-region characterization in terms of eruptive ability or, ultimately,for major solar-flare prediction.Comment: 25 pages, 7 figures, 2 tables, Solar Phys., in pres

Genome-wide association study of eosinophilic granulomatosis with polyangiitis reveals genomic loci stratified by ANCA status

Eosinophilic granulomatosis with polyangiitis (EGPA) is a rare inflammatory disease of unknown cause. 30% of patients have anti-neutrophil cytoplasmic antibodies (ANCA) specific for myeloperoxidase (MPO). Here, we describe a genome-wide association study in 676 EGPA cases and 6,809 controls, that identifies 4 EGPA-associated loci through conventional case-control analysis, and 4 additional associations through a conditional false discovery rate approach. Many variants are also associated with asthma and six are associated with eosinophil count in the general population. Through Mendelian randomisation, we show that a primary tendency to eosinophilia contributes to EGPA susceptibility. Stratification by ANCA reveals that EGPA comprises two genetically and clinically distinct syndromes. MPO+ ANCA EGPA is an eosinophilic autoimmune disease sharing certain clinical features and an HLA-DQ association with MPO+ ANCA-associated vasculitis, while ANCA-negative EGPA may instead have a mucosal/barrier dysfunction origin. Four candidate genes are targets of therapies in development, supporting their exploration in EGPA