Source author record

R. T. James McAteer

R. T. James McAteer 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

26works
10topics
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

26 published item(s)

preprint2017arXiv

The Slowly Varying Corona: I --- Daily Differential Emission Measure Distributions Derived from EVE Spectra

Daily differential emission measure (DEM) distributions of the solar corona are derived from spectra obtained by the Extreme-ultraviolet Variability Experiment (EVE) over a 4-year period starting in 2010 near solar minimum and continuing through the maximum of solar cycle 24. The DEMs are calculated using six strong emission features dominated by Fe lines of charge states VIII, IX, XI, XII, XIV, and XVI that sample the non-flaring coronal temperature range 0.3--5 MK. A proxy for the non-XVIII emission in the wavelength band around the 93.9 Å line is demonstrated. There is little variability in the cool component of the corona (T $<$ 1.3 MK) over the four years, suggesting that the quiet-Sun corona does not respond strongly to the solar cycle, whereas the hotter component (T $>$ 2.0 MK) varies by more than an order of magnitude. A discontinuity in the behavior of coronal diagnostics in 2011 February--March, around the time of the first X-class flare of cycle 24, suggests fundamentally different behavior in the corona under solar minimum and maximum conditions. This global state transition occurs over a period of several months. The DEMs are used to estimate the thermal energy of the visible solar corona (of order $10^{31}$ erg), its radiative energy loss rate (2.5--8 $\times 10^{27}$ erg s$^{-1}$), and the corresponding energy turnover timescale (about an hour). The uncertainties associated with the DEMs and these derived values are mostly due to the coronal Fe abundance and density and the CHIANTI atomic line database.

preprint2016arXiv

Evolution of a Current Sheet in a Solar Flare

We report multi-wavelength and multi-viewpoint observations of a solar eruptive event which involves loop-loop interactions. During a C2.0 flare, motions associated with inflowing and outflowing plasma provide evidence for ongoing magnetic reconnection. The flare loop top and a rising "concave-up" feature are connected by a current-sheet-like structure (CSLS). The physical properties (thickness, length, temperature, and density) of the CSLS are evaluated. In regions adjacent to the CSLS, the EUV emission (characteristic temperature at 1.6 MK) begins to increase more than ten minutes prior to the onset of the flare, and steeply decreases during the decay phase. The reduction of the emission resembles that expected from coronal dimming. The dynamics of this event imply a magnetic reconnection rate in the range 0.01 -- 0.05.

preprint2016arXiv

On the non-Kolmogorov nature of flare-productive solar active regions

A magnetic power spectral analysis is performed on 53 solar active regions, observed from August 2011 to July 2012. Magnetic field data obtained from the Helioseismic and Magnetic Imager, inverted as Active Region Patches, are used to study the evolution of the magnetic power index as each region rotates across the solar disk. Active regions are classified based on the number, and sizes, of solar flares they produce, in order to study the relationship between flare productivity and the magnetic power index. The choice of window size and inertial range plays a key role in determining the correct magnetic power index. The overall distribution of magnetic power indices has a range of $1.0-2.5$. Flare-quiet regions peak at a value of 1.6, however flare-productive regions peak at a value of 2.2. Overall, the histogram of the distribution of power indices of flare-productive active regions is well separated from flare-quiet active regions. Only 12\% of flare-quiet regions exhibit an index greater than 2, whereas 90\% of flare-productive regions exhibit an index greater than 2. Flare-quiet regions exhibit a high temporal variance (i.e, the index fluctuates between high and low values), whereas flare-productive regions maintain an index greater than 2 for several days. This shows the importance of including the temporal evolution of active regions in flare prediction studies, and highlights the potential of a 2-3 day prediction window for space weather applications.

preprint2015arXiv

25 Years of Self-Organized Criticality: Numerical Detection Methods

The detection and characterization of self-organized criticality (SOC), in both real and simulated data, has undergone many significant revisions over the past 25 years. The explosive advances in the many numerical methods available for detecting, discriminating, and ultimately testing, SOC have played a critical role in developing our understanding of how systems experience and exhibit SOC. In this article, methods of detecting SOC are reviewed; from correlations to complexity to critical quantities. A description of the basic autocorrelation method leads into a detailed analysis of application-oriented methods developed in the last 25 years. In the second half of this manuscript space-based, time-based and spatial-temporal methods are reviewed and the prevalence of power laws in nature is described, with an emphasis on event detection and characterization. The search for numerical methods to clearly and unambiguously detect SOC in data often leads us outside the comfort zone of our own disciplines - the answers to these questions are often obtained by studying the advances made in other fields of study. In addition, numerical detection methods often provide the optimum link between simulations and experiments in scientific research. We seek to explore this boundary where the rubber meets the road, to review this expanding field of research of numerical detection of SOC systems over the past 25 years, and to iterate forwards so as to provide some foresight and guidance into developing breakthroughs in this subject over the next quarter of a century.

preprint2015arXiv

An automated classification approach to ranking photospheric proxies of magnetic energy build-up

We study the photospheric magnetic field of ~2000 active regions in solar cycle 23 to search for parameters indicative of energy build-up and subsequent release as a solar flare. We extract three sets of parameters: snapshots in space and time- total flux, magnetic gradients, and neutral lines; evolution in time- flux evolution; structures at multiple size scales- wavelet analysis. This combines pattern recognition and classification techniques via a relevance vector machine to determine whether a region will flare. We consider classification performance using all 38 extracted features and several feature subsets. Classification performance is quantified using both the true positive rate and the true negative rate. Additionally, we compute the true skill score which provides an equal weighting to true positive rate and true negative rate and the Heidke skill score to allow comparison to other flare forecasting work. We obtain a true skill score of ~0.5 for any predictive time window in the range 2-24hr, with a TPR of ~0.8 and a TNR of ~0.7. These values do not appear to depend on the time window, although the Heidke skill score (<0.5) does. Features relating to snapshots of the distribution of magnetic gradients show the best predictive ability over all predictive time windows. Other gradient-related features and the instantaneous power at various wavelet scales also feature in the top five ranked features in predictive power. While the photospheric magnetic field governs the coronal non-potentiality (and likelihood of flaring), photospheric magnetic field alone is not sufficient to determine this uniquely. Furthermore we are only measuring proxies of the magnetic energy build up. We still lack observational details on why energy is released at any particular point in time. We may have discovered the natural limit of the accuracy of flare predictions from these large scale studies.

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

Frozen-in Fractals All Around: Inferring the Large Scale Effects of Small-Scale Magnetic Structure

The large-scale structure of the magnetic field in the solar corona provides the energy to power large-scale solar eruptive events. Our physical understanding of this structure, and hence our ability to predict these events, is limited by the type of data currently available. It is shown that the multifractal spectrum is a powerful tool to study this structure, by providing a physical connection between the details of photospheric magnetic gradients and current density at all size scales. This uses concepts associated with geometric measure theory and the theory of weakly differentiable functions to compare Ampère's law to the wavelet-transform modulus maximum method. The Hölder exponent provides a direct measure of the rate of change of current density across spatial size scales. As this measure is independent of many features of the data (pixel resolution, data size, data type, presence of quiet-Sun data), it provides a unique approach to studying magnetic-field complexity and hence a potentially powerful tool for a statistical prediction of solar-flare activity. Three specific predictions are provided to test this theory: the multifractal spectra will not be dependent on the data type or quality; quiet-Sun gradients will not persist with time; structures with large current densities at large size scale will be the source of energy storage for solar eruptive events.

preprint2015arXiv

Prediction of Solar Flare Size and Time-to-Flare Using Support Vector Machine Regression

We study the prediction of solar flare size and time-to-flare using 38 features describing magnetic complexity of the photospheric magnetic field. This work uses support vector regression to formulate a mapping from the 38-dimensional feature space to a continuous-valued label vector representing flare size or time-to-flare. When we consider flaring regions only, we find an average error in estimating flare size of approximately half a \emph{geostationary operational environmental satellite} (\emph{GOES}) class. When we additionally consider non-flaring regions, we find an increased average error of approximately 3/4 a \emph{GOES} class. We also consider thresholding the regressed flare size for the experiment containing both flaring and non-flaring regions and find a true positive rate of 0.69 and a true negative rate of 0.86 for flare prediction. The results for both of these size regression experiments are consistent across a wide range of predictive time windows, indicating that the magnetic complexity features may be persistent in appearance long before flare activity. This is supported by our larger error rates of some 40 hr in the time-to-flare regression problem. The 38 magnetic complexity features considered here appear to have discriminative potential for flare size, but their persistence in time makes them less discriminative for the time-to-flare problem.

preprint2014arXiv

Coronal Fourier power spectra: implications for coronal seismology and coronal heating

The dynamics of regions of the solar corona are investigated using Atmospheric Imaging Assembly (AIA) 171Å and 193Å data. The coronal emission from the quiet Sun, coronal loop footprints, coronal moss, and from above a sunspot is studied. It is shown that the mean Fourier power spectra in these regions can be described by a power law at lower frequencies that tails to flat spectrum at higher frequencies, plus a Gaussian-shaped contribution that varies depending on the region studied. This Fourier spectral shape is in contrast to the commonly-held assumption that coronal time-series are well described by the sum of a long time-scale background trend plus Gaussian-distributed noise, with some specific locations also showing an oscillatory signal. The implications of this discovery to the field of coronal seismology and the automated detections of oscillations are discussed. The power law contribution to the shape of the Fourier power spectrum is interpreted as being due to the summation of a distribution of exponentially decaying emission events along the line of sight. This is consistent with the idea that the solar atmosphere is heated everywhere by small energy deposition events.

preprint2014arXiv

Qualities of Sequential Chromospheric Brightenings Observed in Optical and UV Images

Chromospheric flare ribbons observed in H-alpha appear well-organized when first examined: ribbons impulsively brighten, morphologically evolve, and exponentially decay back to pre-flare levels. Upon closer inspection of H-alpha flares, there is often a significant number of compact areas brightening in concert with the flare eruption but are spatially separated from the evolving flare ribbon. One class of these brightenings is known as sequential chromospheric brightenings (SCBs). SCBs are often observed in the intimidate vicinity of erupting flares and are associated with coronal mass ejections. In the past decade there have been several previous investigations of SCBs. These studies have exclusively relied upon H-alphaimages to discover and analyze these ephemeral brightenings. This work employs the automated detection algorithm of Kirk et al. (2011) to extract the physical qualities of SCBs in observations of ground-based H-alpha images and complementary AIA images in HeII, Civ, and 1700 Å. The meta-data produced in this tracking process are then culled using complementary Doppler velocities to isolate three distinguishable types of SCBs. From a statistical analysis, we find that the SCBs at the chromospheric H-alpha layer appear earlier, and last longer than their corresponding signatures measured in AIA. From this multi-layer analysis, we infer that SCBs are spatially constrained to the mid-chromosphere. We also derive an energy budget to explain SCBs in which SCBs have a postulated energy of not more than 0.01% of the total flare energy.

preprint2013arXiv

On the sensitivity of the GOES flare classification to properties of the electron beam in the thick target model

The collisional thick target model, wherein a large number of electrons are accelerated down a flaring loop, can be used to explain many observed properties of solar flares. In this study, we focus on the sensitivity of GOES flare classification to the properties of the thick target model. Using a hydrodynamic model with RHESSI-derived electron beam parameters, we explore the effects of the beam energy flux (or total non-thermal energy), the cut-off energy, and the spectral index of the electron distribution on the soft X-rays (SXRs) observed by GOES. We conclude that (1) the GOES class is proportional to the non-thermal energy $E^α$ for $α\approx 1.7$ in the low energy passband (1-8 Å) and $\approx 1.6$ in the high energy passband (0.5-4 Å); (2) the GOES class is only weakly dependent on the spectral index in both passbands; (3) increases in the cut-off will increase the flux in the 0.5-4 Å passband, but decrease the flux in the 1-8 Å passband, while decreases in the cut-off will cause a decrease in the 0.5-4 Å passband and a slight increase in the 1-8 Å passband.

preprint2013arXiv

Remote Oscillatory responses to a solar flare

The processes governing energy storage and release in the Sun are both related to the solar magnetic field. We demonstrate the existence of a magnetic connection between energy released caused by a flare and increased oscillatory power in the lower solar atmosphere. The oscillatory power in active regions tends to increase in response to explosive events at a different location, but not in the region itself. We carry out timing studies and show that this is probably caused by a large scale magnetic connection between the regions, and not a globally propagating wave. We show that oscillations tend to exist in longer lived wave trains at short periods (P< 200s) at the time of a flare. This may be a mechanism by which flare energy can be redistributed throughout the solar atmosphere.

preprint2012arXiv

Coronal Mass Ejection Mass, Energy, and Force Estimates Using STEREO

Understanding coronal mass ejection (CME) energetics and dynamics has been a long-standing problem, and although previous observational estimates have been made, such studies have been hindered by large uncertainties in CME mass. Here, the two vantage points of the Solar Terrestrial Relations Observatory (STEREO) COR1 and COR2 coronagraphs were used to accurately estimate the mass of the 2008 December 12 CME. Acceleration estimates derived from the position of the CME front in 3-D were combined with the mass estimates to calculate the magnitude of the kinetic energy and driving force at different stages of the CME evolution. The CME asymptotically approaches a mass of 3.4\pm1.0x10^15 g beyond ~10 R_sun. The kinetic energy shows an initial rise towards 6.3\pm3.7x10^29 erg at ~3 R_sun, beyond which it rises steadily to 4.2\pm2.5x10^30 erg at ~18 R_sun. The dynamics are described by an early phase of strong acceleration, dominated by a force of peak magnitude of 3.4\pm2.2x10^14N at ~3 R_sun, after which a force of 3.8\pm5.4x10^13 N takes affect between ~7-18 R_sun. These results are consistent with magnetic (Lorentz) forces acting at heliocentric distances of <7 R_sun, while solar wind drag forces dominate at larger distances (>7 R_sun).

preprint2012arXiv

Properties of Sequential Chromospheric Brightenings and Associated Flare Ribbons

We report on the physical properties of solar sequential chromospheric brightenings (SCBs) observed in conjunction with moderate-sized chromospheric flares with associated CMEs. To characterize these ephemeral events, we developed automated procedures to identify and track subsections (kernels) of solar flares and associated SCBs using high resolution H-alpha images. Following the algorithmic identification and a statistical analysis, we compare and find the following: SCBs are distinctly different from flare kernels in their temporal characteristics of intensity, Doppler structure, duration, and location properties. We demonstrate that flare ribbons are themselves made up of subsections exhibiting differing characteristics. Flare kernels are measured to have a mean propagation speed of 0.2 km/s and a maximum speed of 2.3 km/s over a mean distance of 5 x 10^3 km. Within the studied population of SCBs, different classes of characteristics are observed with coincident negative, positive, or both negative and positive Doppler shifts of a few km/s. The appearance of SCBs precede peak flare intensity by ~12 minutes and decay ~1 hour later. They are also found to propagate laterally away from flare center in clusters at 41 km/s or 89 km/s. Given SCBs distinctive nature compared to flares, we suggest a different physical mechanism relating to their origin than the associated flare. We present a heuristic model of the origin of SCBs.

preprint2012arXiv

Sequential Chomospheric Brightening: An Automated Approach to Extracting Physics from Ephemeral Brightening

We make a comparison between small scale chromospheric brightenings and energy release processes through examining the temporal evolution of sequential chromospheric brightenings (SCBs), derive propagation velocities, and propose a connection of the small-scale features to solar flares. Our automated routine detects and distinguishes three separate types of brightening regularly observed in the chromosphere: plage, flare ribbon, and point brightenings. By studying their distinct dynamics, we separate out the flare-associated bright points commonly known as SCBs and identify a propagating Moreton wave. Superimposing our detections on complementary off-band images, we extract a Doppler velocity measurement beneath the point brightening locations. Using these dynamic measurements, we put forward a connection between point brightenings, the erupting flare, and overarching magnetic loops. A destabilization of the pre-flare loop topology by the erupting flare directly leads to the SCBs observed.

preprint2012arXiv

Toward Reliable Benchmarking of Solar Flare Forecasting Methods

Solar flares occur in complex sunspot groups, but it remains unclear how the probability of producing a flare of a given magnitude relates to the characteristics of the sunspot group. Here, we use Geostationary Operational Environmental Satellite X-ray flares and McIntosh group classifications from solar cycles 21 and 22 to calculate average flare rates for each McIntosh class and use these to determine Poisson probabilities for different flare magnitudes. Forecast verification measures are studied to find optimum thresholds to convert Poisson flare probabilities into yes/no predictions of cycle 23 flares. A case is presented to adopt the true skill statistic (TSS) as a standard for forecast comparison over the commonly used Heidke skill score (HSS). In predicting flares over 24 hr, the maximum values of TSS achieved are 0.44 (C-class), 0.53 (M-class), 0.74 (X-class), 0.54 (>=M1.0), and 0.46 (>=C1.0). The maximum values of HSS are 0.38 (C-class), 0.27 (M-class), 0.14 (X-class), 0.28 (>=M1.0), and 0.41 (>=C1.0). These show that Poisson probabilities perform comparably to some more complex prediction systems, but the overall inaccuracy highlights the problem with using average values to represent flaring rate distributions.

preprint2011arXiv

Automated Solar Feature Detection for Space Weather Applications

The solar surface and atmosphere are highly dynamic plasma environments, which evolve over a wide range of temporal and spatial scales. Large-scale eruptions, such as coronal mass ejections, can be accelerated to millions of kilometres per hour in a matter of minutes, making their automated detection and characterisation challenging. Additionally, there are numerous faint solar features, such as coronal holes and coronal dimmings, which are important for space weather monitoring and forecasting, but their low intensity and sometimes transient nature makes them problematic to detect using traditional image processing techniques. These difficulties are compounded by advances in ground- and space- based instrumentation, which have increased the volume of data that solar physicists are confronted with on a minute-by-minute basis; NASA's Solar Dynamics Observatory for example is returning many thousands of images per hour (~1.5 TB/day). This chapter reviews recent advances in the application of images processing techniques to the automated detection of active regions, coronal holes, filaments, CMEs, and coronal dimmings for the purposes of space weather monitoring and prediction.

preprint2011arXiv

Deceleration and Dispersion of Large-scale Coronal Bright Fronts

One of the most dramatic manifestations of solar activity are large-scale coronal bright fronts (CBFs) observed in extreme ultraviolet (EUV) images of the solar atmosphere. To date, the energetics and kinematics of CBFs remain poorly understood, due to the low image cadence and sensitivity of previous EUV imagers and the limited methods used to extract the features. In this paper, the trajectory and morphology of CBFs was determined in order to investigate the varying properties of a sample of CBFs, including their kinematics and pulse shape, dispersion, and dissipation. We have developed a semi-automatic intensity profiling technique to extract the morphology and accurate positions of CBFs in 2.5-10 min cadence images from STEREO/EUVI. The technique was applied to sequences of 171A and 195A images from STEREO/EUVI in order to measure the wave properties of four separate CBF events. Following launch at velocities of ~240-450kms^{-1} each of the four events studied showed significant negative acceleration ranging from ~ -290 to -60ms^{-2}. The CBF spatial and temporal widths were found to increase from ~50 Mm to ~200 Mm and ~100 s to ~1500 s respectively, suggesting that they are dispersive in nature. The variation in position-angle averaged pulse-integrated intensity with propagation shows no clear trend across the four events studied. These results are most consistent with CBFs being dispersive magnetoacoustic waves.

preprint2010arXiv

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.

preprint2010arXiv

Characterising Complexity in Solar Magnetogram Data using a Wavelet-based Segmentation Method

The multifractal nature of solar photospheric magnetic structures are studied using the 2D wavelet transform modulus maxima (WTMM) method. This relies on computing partition functions from the wavelet transform skeleton defined by the WTMM method. This skeleton provides an adaptive space-scale partition of the fractal distribution under study, from which one can extract the multifractal singularity spectrum. We describe the implementation of a multiscale image processing segmentation procedure based on the partitioning of the WT skeleton which allows the disentangling of the information concerning the multifractal properties of active regions from the surrounding quiet-Sun field. The quiet Sun exhibits a average Hölder exponent $\sim -0.75$, with observed multifractal properties due to the supergranular structure. On the other hand, active region multifractal spectra exhibit an average Hölder exponent $\sim 0.38$ similar to those found when studying experimental data from turbulent flows.

preprint2010arXiv

Coronal Mass Ejection Detection using Wavelets, Curvelets and Ridgelets: Applications for Space Weather Monitoring

Coronal mass ejections (CMEs) are large-scale eruptions of plasma and magnetic feld that can produce adverse space weather at Earth and other locations in the Heliosphere. Due to the intrinsic multiscale nature of features in coronagraph images, wavelet and multiscale image processing techniques are well suited to enhancing the visibility of CMEs and supressing noise. However, wavelets are better suited to identifying point-like features, such as noise or background stars, than to enhancing the visibility of the curved form of a typical CME front. Higher order multiscale techniques, such as ridgelets and curvelets, were therefore explored to characterise the morphology (width, curvature) and kinematics (position, velocity, acceleration) of CMEs. Curvelets in particular were found to be well suited to characterising CME properties in a self-consistent manner. Curvelets are thus likely to be of benefit to autonomous monitoring of CME properties for space weather applications.

preprint2010arXiv

Evidence of a Plasmoid-Looptop Interaction and Magnetic Inflows During a Solar Flare/CME Eruptive Event

Observational evidence is presented for the merging of a downward-propagating plasmoid with a looptop kernel during an occulted limb event on 2007 January 25. RHESSI lightcurves in the 9-18 keV energy range, as well as that of the 245 MHz channel of the Learmonth Solar Observatory, show enhanced nonthermal emission in the corona at the time of the merging suggesting that additional particle acceleration took place. This was attributed to a secondary episode of reconnection in the current sheet that formed between the two merging sources. RHESSI images were used to establish a mean downward velocity of the plasmoid of 12 km/s. Complementary observations from the SECCHI suite of instruments onboard STEREO-Behind showed that this process occurred during the acceleration phase of the associated CME. From wavelet-enhanced EUVI, images evidence of inflowing magnetic field lines prior to the CME eruption is also presented. The derived inflow velocity was found to be 1.5 km/s. This combination of observations supports a recent numerical simulation of plasmoid formation, propagation and subsequent particle acceleration due to the tearing mode instability during current sheet formation.

preprint2010arXiv

Propagation of an Earth-directed coronal mass ejection in three dimensions

Solar coronal mass ejections (CMEs) are the most significant drivers of adverse space weather at Earth, but the physics governing their propagation through the heliosphere is not well understood. While stereoscopic imaging of CMEs with the Solar Terrestrial Relations Observatory (STEREO) has provided some insight into their three-dimensional (3D) propagation, the mechanisms governing their evolution remain unclear due to difficulties in reconstructing their true 3D structure. Here we use a new elliptical tie-pointing technique to reconstruct a full CME front in 3D, enabling us to quantify its deflected trajectory from high latitudes along the ecliptic, and measure its increasing angular width and propagation from 2-46 solar radii (approximately 0.2 AU). Beyond 7 solar radii, we show that its motion is determined by an aerodynamic drag in the solar wind and, using our reconstruction as input for a 3D magnetohydrodynamic simulation, we determine an accurate arrival time at the Lagrangian L1 point near Earth.

preprint2010arXiv

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

preprint2009arXiv

Turbulence, Complexity, and Solar Flares

The issue of predicting solar flares is one of the most fundamental in physics, addressing issues of plasma physics, high-energy physics, and modelling of complex systems. It also poses societal consequences, with our ever-increasing need for accurate space weather forecasts. Solar flares arise naturally as a competition between an input (flux emergence and rearrangement) in the photosphere and an output (electrical current build up and resistive dissipation) in the corona. Although initially localised, this redistribution affects neighbouring regions and an avalanche occurs resulting in large scale eruptions of plasma, particles, and magnetic field. As flares are powered from the stressed field rooted in the photosphere, a study of the photospheric magnetic complexity can be used to both predict activity and understand the physics of the magnetic field. The magnetic energy spectrum and multifractal spectrum are highlighted as two possible approaches to this.

preprint2006arXiv

The Influence of Magnetic Field on Oscillations in the Solar Chromosphere

Two sequences of solar images obtained by the Transition Region and Coronal Explorer in three UV passbands are studied using wavelet and Fourier analysis and compared to the photospheric magnetic flux measured by the Michelson Doppler Interferometer on the Solar Heliospheric Observatory to study wave behaviour in differing magnetic environments. Wavelet periods show deviations from the theoretical cutoff value and are interpreted in terms of inclined fields. The variation of wave speeds indicates that a transition from dominant fast-magnetoacoustic waves to slow modes is observed when moving from network into plage and umbrae. This implies preferential transmission of slow modes into the upper atmosphere, where they may lead to heating or be detected in coronal loops and plumes.