Science Platforms for Heliophysics Data Analysis
We recommend that NASA maintain and fund science platforms that enable interactive and scalable data analysis in order to maximize the scientific return of data collected from space-based instruments.
Source author record
Jack Ireland
Jack Ireland 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
10works
9topics
4close collaborators
Research graph
See the researcher in context
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
10 published item(s)
preprint2020arXiv
A Survey of Computational Tools in Solar Physics
The SunPy Project developed a 13-question survey to understand the software and hardware usage of the solar physics community. 364 members of the solar physics community, across 35 countries, responded to our survey. We found that 99$\pm$0.5% of respondents use software in their research and 66% use the Python scientific software stack. Students are twice as likely as faculty, staff scientists, and researchers to use Python rather than Interactive Data Language (IDL). In this respect, the astrophysics and solar physics communities differ widely: 78% of solar physics faculty, staff scientists, and researchers in our sample uses IDL, compared with 44% of astrophysics faculty and scientists sampled by Momcheva and Tollerud (2015). 63$\pm$4% of respondents have not taken any computer-science courses at an undergraduate or graduate level. We also found that most respondents utilize consumer hardware to run software for solar-physics research. Although 82% of respondents work with data from space-based or ground-based missions, some of which (e.g. the Solar Dynamics Observatory and Daniel K. Inouye Solar Telescope) produce terabytes of data a day, 14% use a regional or national cluster, 5% use a commercial cloud provider, and 29% use exclusively a laptop or desktop. Finally, we found that 73$\pm$4% of respondents cite scientific software in their research, although only 42$\pm$3% do so routinely.
0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv
Machine Learning in Heliophysics and Space Weather Forecasting: A White Paper of Findings and Recommendations
The authors of this white paper met on 16-17 January 2020 at the New Jersey Institute of Technology, Newark, NJ, for a 2-day workshop that brought together a group of heliophysicists, data providers, expert modelers, and computer/data scientists. Their objective was to discuss critical developments and prospects of the application of machine and/or deep learning techniques for data analysis, modeling and forecasting in Heliophysics, and to shape a strategy for further developments in the field. The workshop combined a set of plenary sessions featuring invited introductory talks interleaved with a set of open discussion sessions. The outcome of the discussion is encapsulated in this white paper that also features a top-level list of recommendations agreed by participants.
0 citations0 reviewsNo code linkedOpen access
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.
0 citations0 reviewsNo code linkedOpen access
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.
0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv
SunPy - Python for Solar Physics
This paper presents SunPy (version 0.5), a community-developed Python package for solar physics. Python, a free, cross-platform, general-purpose, high-level programming language, has seen widespread adoption among the scientific community, resulting in the availability of a large number of software packages, from numerical computation (NumPy, SciPy) and machine learning (scikit-learn) to visualisation and plotting (matplotlib). SunPy is a data-analysis environment specialising in providing the software necessary to analyse solar and heliospheric data in Python. SunPy is open-source software (BSD licence) and has an open and transparent development workflow that anyone can contribute to. SunPy provides access to solar data through integration with the Virtual Solar Observatory (VSO), the Heliophysics Event Knowledgebase (HEK), and the HELiophysics Integrated Observatory (HELIO) webservices. It currently supports image data from major solar missions (e.g., SDO, SOHO, STEREO, and IRIS), time-series data from missions such as GOES, SDO/EVE, and PROBA2/LYRA, and radio spectra from e-Callisto and STEREO/SWAVES. We describe SunPy's functionality, provide examples of solar data analysis in SunPy, and show how Python-based solar data-analysis can leverage the many existing tools already available in Python. We discuss the future goals of the project and encourage interested users to become involved in the planning and development of SunPy.
0 citations0 reviewsNo code linkedOpen access
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.
0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv
The Radiated Energy Budget of Chromospheric Plasma in a Major Solar Flare Deduced From Multi-Wavelength Observations
This paper presents measurements of the energy radiated by the lower solar atmosphere, at optical, UV, and EUV wavelengths, during an X-class solar flare (SOL2011-02-15T01:56) in response to an injection of energy assumed to be in the form of nonthermal electrons. Hard X-ray observations from RHESSI were used to track the evolution of the parameters of the nonthermal electron distribution to reveal the total power contained in flare accelerated electrons. By integrating over the duration of the impulsive phase, the total energy contained in the nonthermal electrons was found to be $>2\times10^{31}$ erg. The response of the lower solar atmosphere was measured in the free-bound EUV continua of H I (Lyman), He I, and He II, plus the emission lines of He II at 304Å and H I (Ly$α$) at 1216Å by SDO/EVE, the UV continua at 1600Å and 1700Å by SDO/AIA, and the WL continuum at 4504Å, 5550Å, and 6684Å, along with the Ca II H line at 3968Å using Hinode/SOT. The summed energy detected by these instruments amounted to $\sim3\times10^{30}$ erg; about 15% of the total nonthermal energy. The Ly$α$ line was found to dominate the measured radiative losses. Parameters of both the driving electron distribution and the resulting chromospheric response are presented in detail to encourage the numerical modelling of flare heating for this event, to determine the depth of the solar atmosphere at which these line and continuum processes originate, and the mechanism(s) responsible for their generation.
0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv
Estimating the properties of hard X-ray solar flares by constraining model parameters
We compare four different methods of calculating uncertainty estimates in fitting parameterized models to RHESSI X-ray spectra, considering only statistical sources of error. Three of the four methods are based on estimating the scale-size of the minimum in a hypersurface formed by the weighted sum of the squares of the differences between the model fit and the data as a function of the fit parameters, and are implemented as commonly practiced. The fourth method uses Bayesian data analysis and Markov chain Monte Carlo (MCMC) techniques to calculate an uncertainty estimate. Two flare spectra are modeled: one from the GOES X1.3 class flare of 19 January 2005, and the other from the X4.8 flare of 23 July 2002. The four methods give approximately the same uncertainty estimates for the 19 January 2005 spectral fit parameters, but lead to very different uncertainty estimates for the 23 July 2002 spectral fit. This is because each method implements different analyses of the hypersurface, yielding method-dependent results that differ greatly depending on the shape of the hypersurface. For the 23 July 2002 flare data, there is a 95% probability that the low energy cutoff lies below approximately 40 keV, and a 68% probability that it lies in the range 7-36 keV. The low-energy cutoff for the 19 January 2005 flare is more tightly constrained to 107+/-4 keV (68% probability). Using the Bayesian/MCMC approach, we also estimate for the first time probability density functions (PDFs) for the total number of flare accelerated electrons and the energy they carry. For the 23 July 2002 event, these PDFs are asymmetric with long tails orders of magnitude higher than the most probable value, caused by the poorly constrained value of the low-energy cutoff. For the 19 January 2005 flare spectrum, the PDFs for the total number of flare accelerated electrons and their energy are much more symmetric and narrow.
0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv
Automated Detection of Oscillating Regions in the Solar Atmosphere
Recently observed oscillations in the solar atmosphere have been interpreted and modeled as magnetohydrodynamic wave modes. This has allowed the estimation of parameters that are otherwise hard to derive, such as the coronal magnetic-field strength. This work crucially relies on the initial detection of the oscillations, which is commonly done manually. The volume of Solar Dynamics Observatory (SDO) data will make manual detection inefficient for detecting all of the oscillating regions. An algorithm is presented which automates the detection of areas of the solar atmosphere that support spatially extended oscillations. The algorithm identifies areas in the solar atmosphere whose oscillation content is described by a single, dominant oscillation within a user-defined frequency range. The method is based on Bayesian spectral analysis of time-series and image filtering. A Bayesian approach sidesteps the need for an a-priori noise estimate to calculate rejection criteria for the observed signal, and it also provides estimates of oscillation frequency, amplitude and noise, and the error in all these quantities, in a self-consistent way. The algorithm also introduces the notion of quality measures to those regions for which a positive detection is claimed, allowing simple post-detection discrimination by the user. The algorithm is demonstrated on two Transition Region and Coronal Explorer (TRACE) datasets, and comments regarding its suitability for oscillation detection in SDO are made.
0 citations0 reviewsNo code linkedOpen access