Source author record

Peter T. Gallagher

Peter T. Gallagher 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

44works
8topics
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

44 published item(s)

preprint2022arXiv

Interferometric imaging of the type IIIb and U radio bursts observed with LOFAR on 22 August 2017

The Sun is the source of different types of radio bursts that are associated with solar flares, for example. Among the most frequently observed phenomena are type III solar bursts. Their radio images at low frequencies (below 100 MHz) are relatively poorly studied due to the limitations of legacy radio telescopes. We study the general characteristics of types IIIb and U with stria structure solar radio bursts in the frequency range of 20 - 80 MHz, in particular the source size and evolution in different altitudes, as well as the velocity and energy of electron beams responsible for their generation. In this work types IIIb and U with stria structure radio bursts are analyzed using data from the LOFAR telescope including dynamic spectra and imaging observations, as well as data taken in the X-ray range (GOES and RHESSI satellites) and in the extreme ultraviolet (SDO satellite). In this study we determined the source size limited by the actual shape of the contour at particular frequencies of type IIIb and U solar bursts in a relatively wide frequency band from 20 to 80 MHz. Two of the bursts seem to appear at roughly the same place in the studied active region and their source sizes are similar. It is different in the case of another burst, which seems to be related to another part of the magnetic field structure in this active region. The velocities of the electron beams responsible for the generation of the three bursts studied here were also found to be different.

preprint2022arXiv

Searching for stellar flares from low mass stars using ASKAP and TESS

Solar radio emission at low frequencies (<1 GHz) can provide valuable information on processes driving flares and coronal mass ejections (CMEs). Radio emission has been detected from active M dwarf stars, suggestive of much higher levels of activity than previously thought. Observations of active M dwarfs at low frequencies can provide information on the emission mechanism for high energy flares and possible stellar CMEs. Here, we conducted two observations with the Australian Square Kilometre Pathfinder Telescope (ASKAP) totalling 26 hours and scheduled to overlap with the Transiting Exoplanet Survey Satellite (TESS) Sector 36 field, utilising the wide fields of view of both telescopes to search for multiple M dwarfs. We detected variable radio emission in Stokes I centered at 888 MHz from four known active M dwarfs. Two of these sources were also detected with Stokes V circular polarisation. When examining the detected radio emission characteristics, we were not able distinguish between the models for either electron cyclotron maser or gyrosynchrotron emission. These detections add to the growing number of M dwarfs observed with variable low frequency emission.

preprint2021arXiv

LOFAR observations of a jet-driven piston shock in the low solar corona

The Sun produces highly dynamic and eruptive events that can drive shocks through the corona. These shocks can accelerate electrons, which result in plasma emission in the form of a type II radio burst. Despite the large number of type II radio bursts observations, the precise origin of coronal shocks is still subject to investigation. Here we present a well observed solar eruptive event that occurred on 16 October 2015, focusing on a jet observed in the extreme ultraviolet (EUV) by the Atmospheric Imaging Assembly (SDO/AIA), a streamer observed in white-light by the Large Angle and Spectrometric Coronagraph (SOHO/LASCO), and a metric type II radio burst observed by the LOw Frequency Array (LOFAR). LOFAR interferometrically imaged the fundamental and harmonic sources of a type II radio burst and revealed that the sources did not appear to be co-spatial, as would be expected from the plasma emission mechanism. We correct for the separation between the fundamental and harmonic using a model which accounts for scattering of radio waves by electron density fluctuations in a turbulent plasma. This allows us to show the type II radio sources were located $\sim$0.5 R$_\odot$ above the jet and propagated at a speed of $\sim$1000 kms$^{-1}$, which was significantly faster than the jet speed of $\sim$200 kms$^{-1}$. This suggests that the type II burst was generated by a piston shock driven by the jet in the low corona.

preprint2020arXiv

A Comparison of Flare Forecasting Methods. IV. Evaluating Consecutive-Day Forecasting Patterns

A crucial challenge to successful flare prediction is forecasting periods that transition between "flare-quiet" and "flare-active". Building on earlier studies in this series (Barnes et al. 2016; Leka et al. 2019a,b) in which we describe methodology, details, and results of flare forecasting comparison efforts, we focus here on patterns of forecast outcomes (success and failure) over multi-day periods. A novel analysis is developed to evaluate forecasting success in the context of catching the first event of flare-active periods, and conversely, of correctly predicting declining flare activity. We demonstrate these evaluation methods graphically and quantitatively as they provide both quick comparative evaluations and options for detailed analysis. For the testing interval 2016-2017, we determine the relative frequency distribution of two-day dichotomous forecast outcomes for three different event histories (i.e., event/event, no-event/event and event/no-event), and use it to highlight performance differences between forecasting methods. A trend is identified across all forecasting methods that a high/low forecast probability on day-1 remains high/low on day-2 even though flaring activity is transitioning. For M-class and larger flares, we find that explicitly including persistence or prior flare history in computing forecasts helps to improve overall forecast performance. It is also found that using magnetic/modern data leads to improvement in catching the first-event/first-no-event transitions. Finally, 15% of major (i.e., M-class or above) flare days over the testing interval were effectively missed due to a lack of observations from instruments away from the Earth-Sun line.

preprint2020arXiv

A new facility for airborne solar astronomy: NASA's WB-57 at the 2017 total solar eclipse

NASA's WB-57 High Altitude Research Program provides a deployable, mobile, stratospheric platform for scientific research. Airborne platforms are of particular value for making coronal observations during total solar eclipses because of their ability both to follow the Moon's shadow and to get above most of the atmospheric airmass that can interfere with astronomical observations. We used the 2017 Aug 21 eclipse as a pathfinding mission for high-altitude airborne solar astronomy, using the existing high-speed visible-light and near-/mid-wave infrared imaging suite mounted in the WB-57 nose cone. In this paper, we describe the aircraft, the instrument, and the 2017 mission; operations and data acquisition; and preliminary analysis of data quality from the existing instrument suite. We describe benefits and technical limitations of this platform for solar and other astronomical observations. We present a preliminary analysis of the visible-light data quality and discuss the limiting factors that must be overcome with future instrumentation. We conclude with a discussion of lessons learned from this pathfinding mission and prospects for future research at upcoming eclipses, as well as an evaluation of the capabilities of the WB-57 platform for future solar astronomy and general astronomical observation.

preprint2020arXiv

Ensemble Forecasting of Major Solar Flares: Methods for Combining Models

One essential component of operational space weather forecasting is the prediction of solar flares. With a multitude of flare forecasting methods now available online it is still unclear which of these methods performs best, and none are substantially better than climatological forecasts. Space weather researchers are increasingly looking towards methods used by the terrestrial weather community to improve current forecasting techniques. Ensemble forecasting has been used in numerical weather prediction for many years as a way to combine different predictions in order to obtain a more accurate result. Here we construct ensemble forecasts for major solar flares by linearly combining the full-disk probabilistic forecasts from a group of operational forecasting methods (ASAP, ASSA, MAG4, MOSWOC, NOAA, and MCSTAT). Forecasts from each method are weighted by a factor that accounts for the method's ability to predict previous events, and several performance metrics (both probabilistic and categorical) are considered. It is found that most ensembles achieve a better skill metric (between 5\% and 15\%) than any of the members alone. Moreover, over 90\% of ensembles perform better (as measured by forecast attributes) than a simple equal-weights average. Finally, ensemble uncertainties are highly dependent on the internal metric being optimized and they are estimated to be less than 20\% for probabilities greater than 0.2. This simple multi-model, linear ensemble technique can provide operational space weather centres with the basis for constructing a versatile ensemble forecasting system -- an improved starting point to their forecasts that can be tailored to different end-user needs.

preprint2020arXiv

LOFAR observations of radio burst source sizes and scattering in the solar corona

Low frequency radio wave scattering and refraction can have a dramatic effect on the observed size and position of radio sources in the solar corona. The scattering and refraction is thought to be due to fluctuations in electron density caused by turbulence. Hence, determining the true radio source size can provide information on the turbulence in coronal plasma. However, the lack of high spatial resolution radio interferometric observations at low frequencies, such as with the LOw Frequency ARray (LOFAR), has made it difficult to determine the true radio source size and level of radio wave scattering. Here we directly fit the visibilities of a LOFAR observation of a Type IIIb radio burst with an elliptical Gaussian to determine its source size and position. This circumvents the need to image the source and then de-convolve LOFAR's point spread function, which can introduce spurious effects to the source size and shape. For a burst at 34.76 MHz, we find full width at half maximum (FWHM) heights along the major and minor axes to be $18.8^\prime$ $\pm~0.1^\prime$ and $10.2^\prime$ $\pm~0.1^\prime$, respectively, at a plane of sky heliocentric distance of 1.75 R$_\odot$. Our results suggest that the level of density fluctuations in the solar corona is the main cause of the scattering of radio waves, resulting in large source sizes. However, the magnitude of $\varepsilon$ may be smaller than what has been previously derived in observations of radio wave scattering in tied-array images.

preprint2020arXiv

Statistical Study of GOES X-ray Quasi-Periodic Pulsations in Solar Flares

Small amplitude quasi-periodic pulsations (QPPs) detected in soft X-ray emission are commonplace in many flares. To date, the underpinning processes resulting in the QPPs are unknown. In this paper, we attempt to constrain the prevalence of \textit{stationary} QPPs in the largest statistical study to date, including a study of the relationship of QPP periods to the properties of the flaring active region, flare ribbons, and CME affiliation. We build upon the work of \cite{inglis2016} and use a model comparison test to search for significant power in the Fourier spectra of lightcurves of the GOES 1--8~Å channel. We analyze all X-, M- and C- class flares of the past solar cycle, a total of 5519 flares, and search for periodicity in the 6-300~s timescale range. Approximately 46\% of X-class, 29\% of M-class and 7\% of C-class flares show evidence of stationary QPPs, with periods that follow a log-normal distribution peaked at 20~s. The QPP periods were found to be independent of flare magnitude, however a positive correlation was found between QPP period and flare duration. No dependence of the QPP periods to the global active region properties was identified. A positive correlation was found between QPPs and ribbon properties including unsigned magnetic flux, ribbon area and ribbon separation distance. We found that both flares with and without an associated CME can host QPPs. Furthermore, we demonstrate that for X- and M- class flares, decay phase QPPs have statistically longer periods than impulsive phase QPPs.

preprint2019arXiv

Evolution of the Alfvén Mach number associated with a coronal mass ejection shock

The Sun regularly produces large-scale eruptive events, such as coronal mass ejections (CMEs) that can drive shock waves through the solar corona. Such shocks can result in electron acceleration and subsequent radio emission in the form of a type II radio burst. However, the early-phase evolution of shock properties and its relationship to type II burst evolution is still subject to investigation. Here we study the evolution of a CME-driven shock by comparing three commonly used methods of calculating the Alfvén Mach number ($M_A$), namely: shock geometry, a comparison of CME speed to a model of the coronal Alfvén speed, and the type II band-splitting method. We applied the three methods to the 2017 September 2 event, focusing on the shock wave observed in extreme ultraviolet (EUV) by the Solar Ultraviolet Imager (SUVI) on board GOES-16, in white-light by the Large Angle and Spectrometric Coronagraph (LASCO) on board SOHO, and the type II radio burst observed by the Irish Low Frequency Array (I-LOFAR). We show that the three different methods of estimating shock $M_A$ yield consistent results and provide a means of relating shock property evolution to the type II emission duration. The type II radio emission emerged from near the nose of the CME when $M_A$ was in the range 1.4-2.4 at a heliocentric distance of $\sim$1.6 $R_\odot$. The emission ceased when the CME nose reached $\sim$2.4 $R_\odot$, despite an increasing Alfvén Mach number (up to 4). We suggest the radio emission cessation is due to the lack of quasi-perpendicular geometry at this altitude, which inhibits efficient electron acceleration and subsequent radio emission.

preprint2016arXiv

Flaring Rates and the Evolution of Sunspot Group McIntosh Classifications

Sunspot groups are the main source of solar flares, with the energy to power them being supplied by magnetic-field evolution (e.g. flux emergence or twisting/shearing). To date, few studies have investigated the statistical relation between sunspot-group evolution and flaring, with none considering evolution in the McIntosh classification scheme. Here we present a statistical analysis of sunspot groups from Solar Cycle 22, focusing on 24-hour changes in the three McIntosh classification components. Evolution-dependent >C1.0, >M1.0, and >X1.0 flaring rates are calculated, leading to the following results: (i) flaring rates become increasingly higher for greater degrees of upward evolution through the McIntosh classes, with the opposite found for downward evolution; (ii) the highest flaring rates are found for upward evolution from larger, more complex, classes (e.g. Zurich D- and E-classes evolving upward to F-class produce >C1.0 rates of 2.66 +/- 0.28 and 2.31 +/- 0.09 flares per 24 hours, respectively); (iii) increasingly complex classes give higher rates for all flare magnitudes, even when sunspot groups do not evolve over 24 hours. These results support the hypothesis that injection of magnetic energy by flux emergence (i.e. increasing in Zurich or compactness classes) leads to a higher frequency and magnitude of flaring.

preprint2016arXiv

Geomagnetically Induced Currents in the Irish Power Network during Geomagnetic Storms

Geomagnetically induced currents (GICs) are a well-known terrestrial space weather hazard. They occur in power transmission networks and are known to have adverse effects in both high and mid-latitude countries. Here, we study GICs in the Irish power transmission network (geomagnetic latitude 54.7--58.5$^{\circ}$ N) during five geomagnetic storms (06-07 March 2016, 20-21 December 2015, 17-18 March 2015, 29-31 October 2003 and 13-14 March 1989). We simulate electric fields using a plane wave method together with two ground resistivity models, one of which is derived from magnetotelluric measurements (MT model). We then calculate GICs in the 220, 275 and 400~kV transmission network. During the largest of the storm periods studied, the peak electric field was calculated to be as large as 3.8~V~km\textsuperscript{-1}, with associated GICs of up to 23~A using our MT model. Using our homogenous resistivity model, those peak values were 1.46~V~km\textsuperscript{-1} and 25.8~A. We find that three 400 and 275~kV substations are the most likely locations for the Irish transformers to experience large GICs.

preprint2016arXiv

Quasi-Periodic Pulsations during the Impulsive and Decay phases of an X-class Flare

Quasi-periodic pulsations (QPP) are often observed in X-ray emission from solar flares. To date, it is unclear what their physical origins are. Here, we present a multi-instrument investigation of the nature of QPP during the impulsive and decay phases of the X1.0 flare of 28 October 2013. We focus on the character of the fine structure pulsations evident in the soft X-ray time derivatives and compare this variability with structure across multiple wavelengths including hard X-ray and microwave emission. We find that during the impulsive phase of the flare, high correlations between pulsations in the thermal and non-thermal emissions are seen. A characteristic timescale of ~20s is observed in all channels and a second timescale of ~55s is observed in the non-thermal emissions. Soft X-ray pulsations are seen to persist into the decay phase of this flare, up to 20 minutes after the non-thermal emission has ceased. We find that these decay phase thermal pulsations have very small amplitude and show an increase in characteristic timescale from ~40s up to ~70s. We interpret the bursty nature of the co-existing multi-wavelength QPP during the impulsive phase in terms of episodic particle acceleration and plasma heating. The persistent thermal decay phase QPP are most likely connected with compressive MHD processes in the post-flare loops such as the fast sausage mode or the vertical kink mode.

preprint2016arXiv

Radio Diagnostics of Electron Acceleration Sites During the Eruption of a Flux Rope in the Solar Corona

Electron acceleration in the solar corona is often associated with flares and the eruption of twisted magnetic structures known as flux ropes. However, the locations and mechanisms of such particle acceleration during the flare and eruption are still subject to much investigation. Observing the exact sites of particle acceleration can help confirm how the flare and eruption are initiated and how they evolve. Here we use the Atmospheric Imaging Assembly to analyse a flare and erupting flux rope on 2014-April-18, while observations from the Nancay Radio Astronomy Facility allows us to diagnose the sites of electron acceleration during the eruption. Our analysis shows evidence for a pre-formed flux rope which slowly rises and becomes destabilised at the time of a C-class flare, plasma jet and the escape of >75 keV electrons from rope center into the corona. As the eruption proceeds, continued acceleration of electrons with energies of ~5 keV occurs above the flux rope for a period over 5 minutes. At flare peak, one site of electron acceleration is located close to the flare site while another is driven by the erupting flux rope into the corona at speeds of up to 400 km/s. Energetic electrons then fill the erupting volume, eventually allowing the flux rope legs to be clearly imaged from radio sources at 150-445MHz. Following the analysis of Joshi et al. (2015), we conclude that the sites of energetic electrons are consistent with flux rope eruption via a tether-cutting or flux cancellation scenario inside a magnetic fan-spine structure. In total, our radio observations allow us to better understand the evolution of a flux rope eruption and its associated electron acceleration sites, from eruption initiation to propagation into the corona.

preprint2016arXiv

Understanding the Physical Nature of Coronal "EIT Waves"

For almost 20 years the physical nature of globally propagating waves in the solar corona (commonly called "EIT waves") has been controversial and subject to debate. Additional theories have been proposed over the years to explain observations that did not fit with the originally proposed fast-mode wave interpretation. However, the incompatibility of observations made using the Extreme-ultraviolet Imaging Telescope (EIT) onboard the Solar and Heliospheric Observatory with the fast-mode wave interpretation was challenged by differing viewpoints from the twin Solar Terrestrial Relations Observatory spacecraft and higher spatial/temporal resolution data from the Solar Dynamics Observatory. In this article, we reexamine the theories proposed to explain "EIT waves" to identify measurable properties and behaviours that can be compared to current and future observations. Most of us conclude that "EIT waves" are best described as fast-mode large-amplitude waves/shocks that are initially driven by the impulsive expansion of an erupting coronal mass ejection in the low corona.

preprint2015arXiv

Automatic Detection of Magnetic delta in Sunspot Groups

Large and magnetically complex sunspot groups are known to be associated with flares. To date, the Mount Wilson scheme has been used to classify sunspot groups based on their morphological and magnetic properties. The most flare prolific class, the delta sunspot-group, is characterised by opposite polarity umbrae within a common penumbra, separated by less than 2 degrees. In this article, we present a new system, called the Solar Monitor Active Region Tracker - Delta Finder (SMART-DF), that can be used to automatically detect and classify magnetic deltas in near-realtime. Using continuum images and magnetograms from the Helioseismic and Magnetic Imager (HMI) onboard NASA's Solar Dynamics Observatory (SDO), we first estimate distances between opposite polarity umbrae. Opposite polarity pairs having distances of less that 2 degrees are then identified, and if these pairs are found to share a common penumbra, they are identified as a magnetic delta configuration. The algorithm was compared to manual delta detections reported by the Space Weather Prediction Center (SWPC), operated by the National Oceanic and Atmospheric Administration (NOAA). SMART-DF detected 21 out of 23 active regions (ARs) that were marked as delta spots by NOAA during 2011 - 2012 (within +/- 60 degrees longitude). SMART-DF in addition detected five ARs which were not announced as delta spots by NOAA. The near-relatime operation of SMART-DF resulted in many deltas being identified in advance of NOAA's daily notification. SMART-DF will be integrated with SolarMonitor (www.solarmonitor.org) and the near-realtime information will be available to the public.

preprint2015arXiv

Low frequency radio observations of bi-directional electron beams in the solar corona

The radio signature of a shock travelling through the solar corona is known as a type II solar radio burst. In rare cases these bursts can exhibit a fine structure known as `herringbones', which are a direct indicator of particle acceleration occurring at the shock front. However, few studies have been performed on herringbones and the details of the underlying particle acceleration processes are unknown. Here, we use an image processing technique known as the Hough transform to statistically analyse the herringbone fine structure in a radio burst at $\sim$20-90 MHz observed from the Rosse Solar-Terrestrial Observatory on 2011 September 22. We identify 188 individual bursts which are signatures of bi-directional electron beams continuously accelerated to speeds of 0.16$_{-0.10}^{+0.11} c$. This occurs at a shock acceleration site initially at a constant altitude of $\sim$0.6 R$_{\odot}$ in the corona, followed by a shift to $\sim$0.5 R$_{\odot}$. The anti-sunward beams travel a distance of 170$_{-97}^{+174}$ Mm (and possibly further) away from the acceleration site, while those travelling toward the sun come to a stop sooner, reaching a smaller distance of 112$_{-76}^{+84}$ Mm. We show that the stopping distance for the sunward beams may depend on the total number density and the velocity of the beam. Our study concludes that a detailed statistical analysis of herringbone fine structure can provide information on the physical properties of the corona which lead to these relatively rare radio bursts.

preprint2015arXiv

Performance of Major Flare Watches from the Max Millennium Program (2001-2010)

The physical processes that trigger solar flares are not well understood and significant debate remains around processes governing particle acceleration, energy partition, and particle and energy transport. Observations at high resolution in energy, time, and space are required in multiple energy ranges over the whole course of many flares in order to build an understanding of these processes. Obtaining high-quality, co-temporal data from ground- and space- based instruments is crucial to achieving this goal and was the primary motivation for starting the Max Millennium program and Major Flare Watch (MFW) alerts, aimed at coordinating observations of all flares >X1 GOES X-ray classification (including those partially occulted by the limb). We present a review of the performance of MFWs from 1 February 2001 to 31 May 2010, inclusive, that finds: (1) 220 MFWs were issued in 3,407 days considered (6.5% duty cycle), with these occurring in 32 uninterrupted periods that typically last 2-8 days; (2) 56% of flares >X1 were caught, occurring in 19% of MFW days; (3) MFW periods ended at suitable times, but substantial gain could have been achieved in percentage of flares caught if periods had started 24 h earlier; (4) MFWs successfully forecast X-class flares with a true skill statistic (TSS) verification metric score of 0.500, that is comparable to a categorical flare/no-flare interpretation of the NOAA Space Weather Prediction Centre probabilistic forecasts (TSS = 0.488).

preprint2015arXiv

Solar and Heliospheric Physics with the Square Kilometre Array

The fields of solar radiophysics and solar system radio physics, or radio heliophysics, will benefit immensely from an instrument with the capabilities projected for SKA. Potential applications include interplanetary scintillation (IPS), radio-burst tracking, and solar spectral radio imaging with a superior sensitivity. These will provide breakthrough new insights and results in topics of fundamental importance, such as the physics of impulsive energy releases, magnetohydrodynamic oscillations and turbulence, the dynamics of post-eruptive processes, energetic particle acceleration, the structure of the solar wind and the development and evolution of solar wind transients at distances up to and beyond the orbit of the Earth. The combination of the high spectral, time and spatial resolution and the unprecedented sensitivity of the SKA will radically advance our understanding of basic physical processes operating in solar and heliospheric plasmas and provide a solid foundation for the forecasting of space weather events.

preprint2014arXiv

CorPITA: An Automated Algorithm for the Identification and Analysis of Coronal "EIT Waves"

The continuous stream of data available from the Atmospheric Imaging Assembly (AIA) telescopes onboard the Solar Dynamics Observatory (SDO) spacecraft has allowed a deeper understanding of the Sun. However, the sheer volume of data has necessitated the development of automated techniques to identify and analyse various phenomena. In this article, we describe the Coronal Pulse Identification and Tracking Algorithm (CorPITA) for the identification and analysis of coronal "EIT waves". CorPITA uses an intensity-profile technique to identify the propagating pulse, tracking it throughout its evolution before returning estimates of its kinematics. The algorithm is applied here to a data-set from February 2011, allowing its capabilities to be examined and critiqued. This algorithm forms part of the SDO Feature Finding Team initiative and will be implemented as part of the Heliophysics Event Knowledgebase (HEK). This is the first fully automated algorithm to identify and track the propagating "EIT wave" rather than any associated phenomena and will allow a deeper understanding of this controversial phenomenon.

preprint2014arXiv

Quasiperiodic acceleration of electrons by a plasmoid-driven shock in the solar atmosphere

Cosmic rays and solar energetic particles may be accelerated to relativistic energies by shock waves in astrophysical plasmas. On the Sun, shocks and particle acceleration are often associated with the eruption of magnetized plasmoids, called coronal mass ejections (CMEs). However, the physical relationship between CMEs and shock particle acceleration is not well understood. Here, we use extreme ultraviolet, radio and white-light imaging of a solar eruptive event on 22 September 2011 to show that a CME-induced shock (Alfvén Mach number 2.4$^{+0.7}_{-0.8}$) was coincident with a coronal wave and an intense metric radio burst generated by intermittent acceleration of electrons to kinetic energies of 2-46 keV (0.1-0.4 c). Our observations show that plasmoid-driven quasi-perpendicular shocks are capable of producing quasi-periodic acceleration of electrons, an effect consistent with a turbulent or rippled plasma shock surface.

preprint2014arXiv

The Compatibility of Flare Temperatures Observed with AIA, GOES, and RHESSI

We test the compatibility and biases of multi-thermal flare DEM (differential emission measure) peak temperatures determined with AIA with those determined by GOES and RHESSI using the isothermal assumption. In a set of 149 M- and X-class flares observed during the first two years of the SDO mission, AIA finds DEM peak temperatures at the time of the peak GOES 1-8 A flux to have an average of Tp = 12.0+/-2.9 MK and Gaussian DEM widths of log10(sigma_T) = 0.50+/-0.13. From GOES observations of the same 149 events, a mean temperature of Tp = 15.6+/-2.4 MK is inferred, which is systematically higher by a factor of TGOES/TAIA = 1.4+/-0.4. We demonstrate that this discrepancy results from the isothermal assumption in the inversion of the GOES filter ratio. From isothermal fits to photon spectra at energies of E ~ 6-12 keV of 61 of these events, RHESSI finds the temperature to be higher still by a factor of TRHESSI/TAIA = 1.9+/-1.0. We find that this is partly a consequence of the isothermal assumption. However, RHESSI is not sensitive to the low-temperature range of the DEM peak, and thus RHESSI samples only the high-temperature tail of the DEM function. This can also contribute to the discrepancy between AIA and RHESSI temperatures. The higher flare temperatures found by GOES and RHESSI imply correspondingly lower emission measures. We conclude that self-consistent flare DEM temperatures and emission measures require simultaneous fitting of EUV (AIA) and soft X-ray (GOES and RHESSI) fluxes.

preprint2014arXiv

The energetics of a global shock wave in the low solar corona

As the most energetic eruptions in the solar system, coronal mass ejections (CMEs) can produce shock waves at both their front and flanks as they erupt from the Sun into the heliosphere. However, the amount of energy produced in these eruptions, and the proportion of their energy required to produce the waves, is not well characterised. Here we use observations of a solar eruption from 2014 February 25 to estimate the energy budget of an erupting CME and the globally-propagating "EIT wave" produced by the rapid expansion of the CME flanks in the low solar corona. The "EIT wave" is shown using a combination of radio spectra and extreme ultraviolet images to be a shock front with a Mach number greater than one. Its initial energy is then calculated using the Sedov-Taylor blast-wave approximation, which provides an approximation for a shock front propagating through a region of variable density. This approach provides an initial energy estimate of $\approx$2.8$\times$10$^{31}$ ergs to produce the "EIT wave", which is approximately 10% the kinetic energy of the associated CME (shown to be $\approx$2.5$\times$10$^{32}$ ergs). These results indicate that the energy of the "EIT wave" may be significant and must be considered when estimating the total energy budget of solar eruptions.

preprint2014arXiv

The formation heights of coronal shocks from 2D density and Alfvén speed maps

Super-Alfvénic shock waves associated with coronal mass ejections (CMEs) can produce radio emission known as Type II bursts. In the absence of direct imaging, accurate estimates of coronal electron densities, magnetic field strengths and Alfvén speeds are required in order to calculate the kinematics of shocks. To date, 1D radial models have been used, but these are not appropriate for shocks propagating in non-radial directions. Here, we study a coronal shock wave associated with a CME and Type II radio burst using 2D electron density and Alfvén speed maps to determine the locations that shocks are excited as the CME expands through the corona. Coronal density maps were obtained from emission measures derived from the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamic Observatory (SDO) and polarized brightness measurements from the Large Angle and Spectrometric Coronagraph (LASCO) on board the Solar and Heliospheric Observatory (SOHO). Alfvén speed maps were calculated using these density maps and magnetic field extrapolations from the Helioseismic and Magnetic Imager (SDO/HMI). The computed density and Alfvén speed maps were then used to calculate the shock kinematics in non-radial directions. Using the kinematics of the Type II burst and associated shock, we find our observations to be consistent with the formation of a shock located at the CME flanks where the Alfvén speed has a local minimum. 1D density models are not appropriate for shocks that propagate non-radially along the flanks of a CME. Rather, the 2D density, magnetic field and Alfvén speed maps described here give a more accurate method for determining the fundamental properties of shocks and their relation to CMEs.

preprint2014arXiv

Understanding CME and associated shock in the solar corona by merging multi wavelengths observation

Using multi-wavelength imaging observations, in EUV, white light and radio, and radio spectral data over a large frequency range, we analyzed the triggering and development of a complex eruptive event. This one includes two components, an eruptive jet and a CME which interact during more than 30 min, and can be considered as physically linked. This was an unusual event. The jet is generated above a typical complex magnetic configuration which has been investigated in many former studies related to the build-up of eruptive jets; this configuration includes fan-field lines originating from a corona null point above a parasitic polarity, which is embedded in one polarity region of large Active Region (AR). The initiation and development of the CME, observed first in EUV, does not show usual signatures. In this case, the eruptive jet is the main actor of this event. The CME appears first as a simple loop system which becomes destabilized by magnetic reconnection between the outer part of the jet and the ambient medium. The progression of the CME is closely associated with the occurrence of two successive types II bursts from distinct origin. An important part of this study is the first radio type II burst for which the joint spectral and imaging observations allowed: i) to follow, step by step, the evolution of the spectrum and of the trajectory of the radio burst, in relationship with the CME evolution; ii) to obtain, without introducing an electronic density model, the B-field and the Alfven speed.

preprint2013arXiv

The Thermal Properties of Solar Flares Over Three Solar Cycles Using GOES X-ray Observations

Solar flare X-ray emission results from rapidly increasing temperatures and emission measures in flaring active region loops. To date, observations from the X-Ray Sensor (XRS) onboard the Geostationary Operational Environmental Satellite (GOES) have been used to derive these properties, but have been limited by a number of factors, including the lack of a consistent background subtraction method capable of being automatically applied to large numbers of flares. In this paper, we describe an automated temperature and emission measure-based background subtraction method (TEBBS), which builds on the methods of Bornmann (1990). Our algorithm ensures that the derived temperature is always greater than the instrumental limit and the pre-flare background temperature, and that the temperature and emission measure are increasing during the flare rise phase. Additionally, TEBBS utilizes the improved estimates of GOES temperatures and emission measures from White et al. (2005). TEBBS was successfully applied to over 50,000 solar flares occurring over nearly three solar cycles (1980-2007), and used to create an extensive catalog of the solar flare thermal properties. We confirm that the peak emission measure and total radiative losses scale with background subtracted GOES X-ray flux as power-laws, while the peak temperature scales logarithmically. As expected, the peak emission measure shows an increasing trend with peak temperature, although the total radiative losses do not. While these results are comparable to previous studies, we find that flares of a given GOES class have lower peak temperatures and higher peak emission measures than previously reported. The resulting TEBBS database of thermal flare plasma properties is publicly available on Solar Monitor (www.solarmonitor.org/TEBBS/) and will be available on Heliophysics Integrated Observatory (www.helio-vo.eu).

preprint2012arXiv

Automatic Detection and Tracking of CMEs II: Multiscale Filtering of Coronagraph Data

Studying CMEs in coronagraph data can be challenging due to their diffuse structure and transient nature, and user-specific biases may be introduced through visual inspection of the images. The large amount of data available from the SOHO, STEREO, and future coronagraph missions, also makes manual cataloguing of CMEs tedious, and so a robust method of detection and analysis is required. This has led to the development of automated CME detection and cata- loguing packages such as CACTus, SEEDS and ARTEMIS. Here we present the development of a new CORIMP (coronal image processing) CME detection and tracking technique that overcomes many of the drawbacks of current catalogues. It works by first employing the dynamic CME separation technique outlined in a companion paper, and then characterising CME structure via a multiscale edge-detection algorithm. The detections are chained through time to determine the CME kinematics and morphological changes as it propagates across the plane-of-sky. The effectiveness of the method is demonstrated by its application to a selection of SOHO/LASCO and STEREO/SECCHI images, as well as to synthetic coronagraph images created from a model corona with a variety of CMEs. The algorithms described in this article are being applied to the whole LASCO and SECCHI datasets, and a catalogue of results will soon be available to the public.

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

Evidence for Partial Taylor Relaxation from Changes in Magnetic Geometry and Energy during a Solar Flare

Solar flares are powered by energy stored in the coronal magnetic field, a portion of which is released when the field reconfigures into a lower energy state. Investigation of sunspot magnetic field topology during flare activity is useful to improve our understanding of flaring processes. Here we investigate the deviation of the non-linear field configuration from that of the linear and potential configurations, and study the free energy available leading up to and after a flare. The evolution of the magnetic field in NOAA region 10953 was examined using data from Hinode/SOT-SP, over a period of 12 hours leading up to and after a GOES B1.0 flare. Previous work on this region found pre- and post-flare changes in photospheric vector magnetic field parameters of flux elements outside the primary sunspot. 3D geometry was thus investigated using potential, linear force-free, and non-linear force-free field extrapolations in order to fully understand the evolution of the field lines. Traced field line geometrical and footpoint orientation differences show that the field does not completely relax to a fully potential or linear force-free state after the flare. Magnetic and free magnetic energies increase significantly ~ 6.5-2.5 hours before the flare by ~ 10^31 erg. After the flare, the non-linear force-free magnetic energy and free magnetic energies decrease but do not return to pre-flare 'quiet' values. The post-flare non-linear force-free field configuration is closer (but not equal) to that of the linear force-free field configuration than a potential one. However, the small degree of similarity suggests that partial Taylor relaxation has occurred over a time scale of ~ 3-4 hours.

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

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.

preprint2011arXiv

LEMUR: Large European Module for solar Ultraviolet Research. European contribution to JAXA's Solar-C mission

Understanding the solar outer atmosphere requires concerted, simultaneous solar observations from the visible to the vacuum ultraviolet (VUV) and soft X-rays, at high spatial resolution (between 0.1" and 0.3"), at high temporal resolution (on the order of 10 s, i.e., the time scale of chromospheric dynamics), with a wide temperature coverage (0.01 MK to 20 MK, from the chromosphere to the flaring corona), and the capability of measuring magnetic fields through spectropolarimetry at visible and near-infrared wavelengths. Simultaneous spectroscopic measurements sampling the entire temperature range are particularly important. These requirements are fulfilled by the Japanese Solar-C mission (Plan B), composed of a spacecraft in a geosynchronous orbit with a payload providing a significant improvement of imaging and spectropolarimetric capabilities in the UV, visible, and near-infrared with respect to what is available today and foreseen in the near future. The Large European Module for solar Ultraviolet Research (LEMUR), described in this paper, is a large VUV telescope feeding a scientific payload of high-resolution imaging spectrographs and cameras. LEMUR consists of two major components: a VUV solar telescope with a 30 cm diameter mirror and a focal length of 3.6 m, and a focal-plane package composed of VUV spectrometers covering six carefully chosen wavelength ranges between 17 and 127 nm. The LEMUR slit covers 280" on the Sun with 0.14" per pixel sampling. In addition, LEMUR is capable of measuring mass flows velocities (line shifts) down to 2 km/s or better. LEMUR has been proposed to ESA as the European contribution to the Solar C mission.

preprint2011arXiv

Short-term evolution of coronal hole boundaries

The interaction of open and closed field lines at coronal hole boundaries is widely accepted to be due to interchange magnetic reconnection. To date, it is unclear how the boundaries vary on short timescales and at what velocity this occurs. Here, we describe an automated boundary tracking method used to determine coronal hole boundary displacements on short timescales. The bound- ary displacements were found to be isotropic and to have typical expansion/contraction speeds of \leq2 km s^-1, which indicate magnetic reconnection rates of \leq 3 \times 10^-3. The observed displacements were used in conjunction with the interchange reconnection model to derive typical diffusion coeffi- cients of \leq 3 \times 10^13 cm^2 s^-1. These results are consistent with an interchange reconnection process in the low corona driven by the random granular motion of open and closed fields in the photosphere.

preprint2011arXiv

STEREO Direct Imaging of a CME-driven Shock to 0.5AU

Fast coronal mass ejections (CMEs) generate standing or bow shocks as they propagate through the corona and solar wind. Although CME shocks have previously been detected indirectly via their emission at radio frequencies, direct imaging has remained elusive due to their low contrast at optical wavelengths. Here we report the first images of a CME-driven shock as it propagates through interplanetary space from 8R_Sun to 120R_Sun (0.5AU), using observations from the STEREO Heliospheric Imager (HI). The CME was measured to have a velocity of ~1000 km/s and a Mach number of 4.1\pm1.2, while the shock front stand-off distance (Delta) was found to increase linearly to ~20 R_Sun at 0.5 AU. The normalised standoff distance (Delta/DO) showed reasonable agreement with semi-empirical relations, where D_O is the CME radius. However, when normalised using the radius of curvature, Delta/R_O did not agree well with theory, implying that R_O was under-estimated by a factor of ~3-8. This is most likely due to the difficulty in estimating the larger radius of curvature along the CME axis from the observations, which provide only a cross-sectional view of the CME.

preprint2011arXiv

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.

preprint2011arXiv

The Wave Properties of Coronal Bright Fronts Observed Using SDO/AIA

Coronal bright fronts (CBFs) are large scale wavefronts that propagate though the solar corona at hundreds of kilometers per second. While their kinematics have been studied in detail, many questions remain regarding the temporal evolution of their amplitude and pulse width. Here, contemporaneous high cadence, multi-thermal observations of the solar corona from the Solar Dynamic Observatory (SDO) and Solar TErrestrial RElations Observatory (STEREO) spacecraft are used to determine the kinematics and expansion rate of a CBF wavefront observed on 2010 August 14. The CBF was found to have a lower initial velocity with weaker deceleration in STEREO observations compared to SDO (~340 km/s and -72 m/s/s as opposed to ~410 km/s and -279 m/s/s). The CBF kinematics from SDO were found to be highly passband-dependent, with an initial velocity ranging from 379+/-12 km/s to 460+/-28 km/s and acceleration ranging from -128+/-28 m/s/s to -431+/-86 m/s/s in the 335A and 304A passbands respectively. These kinematics were used to estimate a quiet coronal magnetic field strength range of ~1-2 G. Significant pulse broadening was also observed, with expansion rates of ~130 km/s (STEREO) and ~220 km/s (SDO). By treating the CBF as a linear superposition of sinusoidal waves within a Gaussian envelope, the resulting dispersion rate of the pulse was found to be ~8-13 Mm^2 s^-1. These results are indicative of a fast-mode magnetoacoustic wave pulse propagating through an inhomogeneous medium.

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

Constraining 3D Magnetic Field Extrapolations Using The Twin Perspectives of STEREO

The 3D magnetic topology of a solar active region (NOAA 10956) was reconstructed using a linear force-free field extrapolation constrained using the twin perspectives of \emph{STEREO}. A set of coronal field configurations was initially generated from extrapolations of the photospheric magnetic field observed by the Michelson Doppler Imager (MDI) on \emph{SOHO}. Using an EUV intensity-based cost function, the extrapolated field lines that were most consistent with 171Å passband images from the Extreme UltraViolet Imager (EUVI) on \emph{STEREO} were identified. This facilitated quantitative constraints to be placed on the twist ($α$) of the extrapolated field lines, where $\nabla \times {\bf B} = α{\bf B}$. Using the constrained values of $α$, the evolution in time of twist, connectivity, and magnetic energy were then studied. A flux emergence event was found to result in significant changes in the magnetic topology and total magnetic energy of the region.

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

Investigating the driving mechanisms of coronal mass ejections

The objective of this investigation was to first examine the kinematics of coronal mass ejections (CMEs) using EUV and coronagraph images, and then to make a comparison with theoretical models in the hope to identify the driving mechanisms of the CMEs. We have studied two CMEs which occurred on 2006 Dec. 17 (CME06) and 2007 Dec. 31 (CME07). The models studied in this work were catastrophe, breakout, and toroidal instability models. We found that after the eruption, the accelerations of both events exhibited a drop before increasing again. Our comparisons with the theories suggested that CME06 can be best described by a hybrid of the catastrophe and breakout models while CME07 is most consistent with the breakout model.

preprint2010arXiv

Large-scale Bright Fronts in the Solar Corona: A Review of "EIT waves"

``EIT waves" are large-scale coronal bright fronts (CBFs) that were first observed in 195 Å images obtained using the Extreme-ultraviolet Imaging Telescope (EIT) onboard the \emph{Solar and Heliospheric Observatory (SOHO)}. Commonly called ``EIT waves", CBFs typically appear as diffuse fronts that propagate pseudo-radially across the solar disk at velocities of 100--700 km s$^{-1}$ with front widths of 50-100 Mm. As their speed is greater than the quiet coronal sound speed ($c_s\leq$200 km s$^{-1}$) and comparable to the local Alfvén speed ($v_A\leq$1000 km s$^{-1}$), they were initially interpreted as fast-mode magnetoacoustic waves ($v_{f}=(c_s^2 + v_A^2)^{1/2}$). Their propagation is now known to be modified by regions where the magnetosonic sound speed varies, such as active regions and coronal holes, but there is also evidence for stationary CBFs at coronal hole boundaries. The latter has led to the suggestion that they may be a manifestation of a processes such as Joule heating or magnetic reconnection, rather than a wave-related phenomena. While the general morphological and kinematic properties of CBFs and their association with coronal mass ejections have now been well described, there are many questions regarding their excitation and propagation. In particular, the theoretical interpretation of these enigmatic events as magnetohydrodynamic waves or due to changes in magnetic topology remains the topic of much debate.

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

preprint2010arXiv

Solar Wind Drag and the Kinematics of Interplanetary Coronal Mass Ejections

Coronal mass ejections (CMEs) are large-scale ejections of plasma and magnetic field from the solar corona, which propagate through interplanetary space at velocities of $\sim$100--2500~km~s$^{-1}$. Although plane-of-sky coronagraph measurements have provided some insight into their kinematics near the Sun ($<$32~R$_\odot$), it is still unclear what forces govern their evolution during both their early acceleration and later propagation. Here, we use the dual perspectives of the Solar TErrestrial RElations Observatory (STEREO) spacecrafts to derive the three-dimensional kinematics of CMEs over a range of heliocentric distances ($\sim$2--250\,R$_{\odot}$). We find evidence for solar wind (SW) drag-forces acting in interplanetary space, with a fast CME decelerated and a slow CME accelerated towards typical SW velocities. We also find that the fast CME showed linear ($δ=1$) dependence on the velocity difference between the CME and the SW, while the slow CME showed a quadratic ($δ=2$) dependence. The differing forms of drag for the two CMEs indicate the forces and thus mechanism responsible for there acceleration may be different.

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.