Showing 1–50 of 119 results for author: Barnes, G

Coronal Heating as Determined by the Solar Flare Frequency Distribution Obtained by Aggregating Case Studies

Authors: James Paul Mason, Alexandra Werth, Colin G. West, Allison A. Youngblood, Donald L. Woodraska, Courtney Peck, Kevin Lacjak, Florian G. Frick, Moutamen Gabir, Reema A. Alsinan, Thomas Jacobsen, Mohammad Alrubaie, Kayla M. Chizmar, Benjamin P. Lau, Lizbeth Montoya Dominguez, David Price, Dylan R. Butler, Connor J. Biron, Nikita Feoktistov, Kai Dewey, N. E. Loomis, Michal Bodzianowski, Connor Kuybus, Henry Dietrick, Aubrey M. Wolfe , et al. (977 additional authors not shown)

Abstract: Flare frequency distributions represent a key approach to addressing one of the largest problems in solar and stellar physics: determining the mechanism that counter-intuitively heats coronae to temperatures that are orders of magnitude hotter than the corresponding photospheres. It is widely accepted that the magnetic field is responsible for the heating, but there are two competing mechanisms th… ▽ More Flare frequency distributions represent a key approach to addressing one of the largest problems in solar and stellar physics: determining the mechanism that counter-intuitively heats coronae to temperatures that are orders of magnitude hotter than the corresponding photospheres. It is widely accepted that the magnetic field is responsible for the heating, but there are two competing mechanisms that could explain it: nanoflares or Alfvén waves. To date, neither can be directly observed. Nanoflares are, by definition, extremely small, but their aggregate energy release could represent a substantial heating mechanism, presuming they are sufficiently abundant. One way to test this presumption is via the flare frequency distribution, which describes how often flares of various energies occur. If the slope of the power law fitting the flare frequency distribution is above a critical threshold, $α=2$ as established in prior literature, then there should be a sufficient abundance of nanoflares to explain coronal heating. We performed $>$600 case studies of solar flares, made possible by an unprecedented number of data analysts via three semesters of an undergraduate physics laboratory course. This allowed us to include two crucial, but nontrivial, analysis methods: pre-flare baseline subtraction and computation of the flare energy, which requires determining flare start and stop times. We aggregated the results of these analyses into a statistical study to determine that $α= 1.63 \pm 0.03$. This is below the critical threshold, suggesting that Alfvén waves are an important driver of coronal heating. △ Less

Submitted 9 May, 2023; originally announced May 2023.

Comments: 1,002 authors, 14 pages, 4 figures, 3 tables, published by The Astrophysical Journal on 2023-05-09, volume 948, page 71

Implications of Different Solar Photospheric Flux-Transport Models for Global Coronal and Heliospheric Modeling

Authors: Graham Barnes, Marc L. DeRosa, Shaela I. Jones, Charles N. Arge, Carl J. Henney, Mark C. M. Cheung

Abstract: The concept of surface-flux transport (SFT) is commonly used in evolving models of the large-scale solar surface magnetic field. These photospheric models are used to determine the large-scale structure of the overlying coronal magnetic field, as well as to make predictions about the fields and flows that structure the solar wind. We compare predictions from two SFT models for the solar wind, open… ▽ More The concept of surface-flux transport (SFT) is commonly used in evolving models of the large-scale solar surface magnetic field. These photospheric models are used to determine the large-scale structure of the overlying coronal magnetic field, as well as to make predictions about the fields and flows that structure the solar wind. We compare predictions from two SFT models for the solar wind, open magnetic field footpoints, and the presence of coronal magnetic null points throughout various phases of a solar activity cycle, focusing on the months of April in even-numbered years between 2012 and 2020, inclusive. We find that there is a solar cycle dependence to each of the metrics considered, but there is not a single phase of the cycle in which all the metrics indicate good agreement between the models. The metrics also reveal large, transient differences between the models when a new active region is rotating into the assimilation window. The evolution of the surface flux is governed by a combination of large scale flows and comparatively small scale motions associated with convection. Because the latter flows evolve rapidly, there are intervals during which their impact on the surface flux can only be characterized in a statistical sense, thus their impact is modeled by introducing a random evolution that reproduces the typical surface flux evolution. We find that the differences between the predicted properties are dominated by differences in the model assumptions and implementation, rather than selection of a particular realization of the random evolution. △ Less

Submitted 13 February, 2023; originally announced February 2023.

Comments: Accepted for publication in The Astrophysical Journal

Properties of Flare-Imminent versus Flare-Quiet Active Regions from the Chromosphere through the Corona II: NonParametric Discriminant Analysis Results from the NWRA Classification Infrastructure (NCI)

Authors: KD Leka, Karin Dissauer, Graham Barnes, Eric L. Wagner

Abstract: A large sample of active-region-targeted time-series images from the Solar Dynamics Observatory / Atmospheric Imaging Assembly, the AIA Active Region Patch database ("AARPs", Paper I: Dissauer et al 2022) is used to investigate whether parameters describing the coronal, transition region, and chromospheric emission can differentiate a region that will imminently produce a solar flare from one that… ▽ More A large sample of active-region-targeted time-series images from the Solar Dynamics Observatory / Atmospheric Imaging Assembly, the AIA Active Region Patch database ("AARPs", Paper I: Dissauer et al 2022) is used to investigate whether parameters describing the coronal, transition region, and chromospheric emission can differentiate a region that will imminently produce a solar flare from one that will not. Parametrizations based on moment analysis of direct and running-difference images provide for physically-interpretable results from nonparametric discriminant analysis. Across four event definitions including both 24hr and 6hr validity periods, 160 image-based parameters capture the general state of the atmosphere, rapid brightness changes, and longer-term intensity evolution. We find top Brier Skill Scores in the 0.07--0.33 range, True Skill Statistics in the 0.68--0.82 range (both depending on event definition), and Receiver Operating Characteristic Skill Scores above 0.8. Total emission can perform notably as can steeply increasing or decreasing brightness, although mean brightness measures do not, demonstrating the well-known active-region-size/flare-productivity relation. Once a region is flare productive, the active-region coronal plasma appears to stay hot. The 94AA filter data provides the most parameters with discriminating power, with indications that it benefits from sampling multiple physical regimes. In particular, classification success using higher-order moments of running difference images indicate a propensity for flare-imminent regions to display short-lived small-scale brightening events. Parameters describing the evolution of the corona can provide flare-imminent indicators, but at no preference over "static" parameters. Finally, all parameters and NPDA-derived probabilities are available to the community for additional research. △ Less

Submitted 21 December, 2022; originally announced December 2022.

Comments: Astrophysical Journal, in press

Large-Scale Spatial Cross-Calibration of Hinode/SOT-SP and SDO/HMI

Authors: David F. Fouhey, Richard E. L. Higgins, Spiro K. Antiochos, Graham Barnes, Marc L. DeRosa, J. Todd Hoeksema, K. D. Leka, Yang Liu, Peter W. Schuck, Tamas I. Gombosi

Abstract: We investigate the cross-calibration of the Hinode/SOT-SP and SDO/HMI instrument meta-data, specifically the correspondence of the scaling and pointing information. Accurate calibration of these datasets gives the correspondence needed by inter-instrument studies and learning-based magnetogram systems, and is required for physically-meaningful photospheric magnetic field vectors. We approach the p… ▽ More We investigate the cross-calibration of the Hinode/SOT-SP and SDO/HMI instrument meta-data, specifically the correspondence of the scaling and pointing information. Accurate calibration of these datasets gives the correspondence needed by inter-instrument studies and learning-based magnetogram systems, and is required for physically-meaningful photospheric magnetic field vectors. We approach the problem by robustly fitting geometric models on correspondences between images from each instrument's pipeline. This technique is common in computer vision, but several critical details are required when using scanning slit spectrograph data like Hinode/SOT-SP. We apply this technique to data spanning a decade of the Hinode mission. Our results suggest corrections to the published Level 2 Hinode/SOT-SP data. First, an analysis on approximately 2,700 scans suggests that the reported pixel size in Hinode/SOT-SP Level 2 data is incorrect by around 1%. Second, analysis of over 12,000 scans show that the pointing information is often incorrect by dozens of arcseconds with a strong bias. Regression of these corrections indicates that thermal effects have caused secular and cyclic drift in Hinode/SOT-SP pointing data over its mission. We offer two solutions. First, direct co-alignment with SDO/HMI data via our procedure can improve alignments for many Hinode/SOT-SP scans. Second, since the pointing errors are predictable, simple post-hoc corrections can substantially improve the pointing. We conclude by illustrating the impact of this updated calibration on derived physical data products needed for research and interpretation. Among other things, our results suggest that the pointing errors induce a hemispheric bias in estimates of radial current density. △ Less

Submitted 29 September, 2022; originally announced September 2022.

Comments: Under revisions at ApJS

The eta Aquilae System: Radial Velocities and Astrometry in Search of eta Aql B

Authors: G. Fritz Benedict, Thomas G. Barnes III, Nancy R. Evans, William D. Cochran, Richard I. Anderson, Barbara E. McArthur, Thomas E. Harrison

Abstract: The classical Cepheid eta Aql was not included in past Leavitt Law work (Benedict et al. 2007) because of a presumed complicating orbit due to a known B9.8V companion. To determine the orbit of eta Aql B, we analyze a significant number of radial velocity measures (RV) from eight sources. With these we establish the RV variation due to Cepheid pulsation, using a twelve Fourier coefficient model, w… ▽ More The classical Cepheid eta Aql was not included in past Leavitt Law work (Benedict et al. 2007) because of a presumed complicating orbit due to a known B9.8V companion. To determine the orbit of eta Aql B, we analyze a significant number of radial velocity measures (RV) from eight sources. With these we establish the RV variation due to Cepheid pulsation, using a twelve Fourier coefficient model, while solving for velocity offsets required to bring the RV data sets into coincidence. RV residuals provide no evidence of orbital motion, suggesting either nearly face-on orientation or very long period. Reanalysis of Hubble Space Telescope Fine Guidance Sensor astrometry now includes reference star parallax and proper motion priors from Gaia EDR3. As modeling confirmation, we reanalyze zeta Gem in parallel, deriving zeta Gem parallax and proper motion values consistent with Gaia EDR3, and consistent with the Benedict 2007 Leavitt Law. In an effort to further characterize eta Aql B, we hypothesize that eta Aql residuals larger than those of the associated reference stars or a parallax inconsistent with EDR3 and the Benedict 2007 Leavitt Law indicate unmodeled orbital motion. Using the astrometric noise or parallax mismatch with EDR3 we estimate possible periods and mass for eta Aql B. Ascribing photocenter motion to the photometric variation of the Cepheid, eta Aql A, yields a plausible separation, consistent with a long period, explaining the lack of RV variation. None of these approaches yields an unassailable characterization of the eta Aql A-B system △ Less

Submitted 20 April, 2022; originally announced April 2022.

Comments: to appear in The Astronomical Journal

SynthIA: A Synthetic Inversion Approximation for the Stokes Vector Fusing SDO and Hinode into a Virtual Observatory

Authors: Richard E. L. Higgins, David F. Fouhey, Spiro K. Antiochos, Graham Barnes, Mark C. M. Cheung, J. Todd Hoeksema, KD Leka, Yang Liu, Peter W. Schuck, Tamas I. Gombosi

Abstract: Both NASA's Solar Dynamics Observatory (SDO) and the JAXA/NASA Hinode mission include spectropolarimetric instruments designed to measure the photospheric magnetic field. SDO's Helioseismic and Magnetic Imager (HMI) emphasizes full-disk high-cadence and good spatial resolution data acquisition while Hinode's Solar Optical Telescope Spectro-Polarimeter (SOT-SP) focuses on high spatial resolution an… ▽ More Both NASA's Solar Dynamics Observatory (SDO) and the JAXA/NASA Hinode mission include spectropolarimetric instruments designed to measure the photospheric magnetic field. SDO's Helioseismic and Magnetic Imager (HMI) emphasizes full-disk high-cadence and good spatial resolution data acquisition while Hinode's Solar Optical Telescope Spectro-Polarimeter (SOT-SP) focuses on high spatial resolution and spectral sampling at the cost of a limited field of view and slower temporal cadence. This work introduces a deep-learning system named SynthIA (Synthetic Inversion Approximation), that can enhance both missions by capturing the best of each instrument's characteristics. We use SynthIA to produce a new magnetogram data product, SynodeP (Synthetic Hinode Pipeline), that mimics magnetograms from the higher spectral resolution Hinode/SOT-SP pipeline, but is derived from full-disk, high-cadence, and lower spectral-resolution SDO/HMI Stokes observations. Results on held-out data show that SynodeP has good agreement with the Hinode/SOT-SP pipeline inversions, including magnetic fill fraction, which is not provided by the current SDO/HMI pipeline. SynodeP further shows a reduction in the magnitude of the 24-hour oscillations present in the SDO/HMI data. To demonstrate SynthIA's generality, we show the use of SDO/AIA data and subsets of the HMI data as inputs, which enables trade-offs between fidelity to the Hinode/SOT-SP inversions, number of observations used, and temporal artifacts. We discuss possible generalizations of SynthIA and its implications for space weather modeling. This work is part of the NASA Heliophysics DRIVE Science Center (SOLSTICE) at the University of Michigan under grant NASA 80NSSC20K0600E, and will be open-sourced. △ Less

Submitted 27 August, 2021; originally announced August 2021.

Pulsational instability of pre-main sequence models from accreting protostars. I. Constraining the input physics for accretion with spectroscopic parameters and stellar pulsations

Authors: T. Steindl, K. Zwintz, T. G. Barnes, M. Muellner, E. I. Vorobyov

Abstract: Context: The pre-main sequence evolution is often simplified by choosing classical initial models. These have large initial radii and sufficient uniform contraction to make them fully convective. Contrary to that, real stars are born as small protostellar seeds in collapsing molecular clouds and obtain their final mass by means of accretion. Aims: We aim to constrain the input physics of accretion… ▽ More Context: The pre-main sequence evolution is often simplified by choosing classical initial models. These have large initial radii and sufficient uniform contraction to make them fully convective. Contrary to that, real stars are born as small protostellar seeds in collapsing molecular clouds and obtain their final mass by means of accretion. Aims: We aim to constrain the input physics of accretion on protostellar seeds with observed spectroscopic parameters and stellar pulsations of young stellar objects and pre-main sequence stars. Methods: We conducted a literature search for spectroscopic samples of young stellar objects and pre-main sequence stars including all previously known pulsators. The sample size of pulsating pre-main sequence stars is increased by analysing TESS observations and presenting discoveries in CoRoT data. We employ MESA and GYRE to calculate evolutionary tracks of accreting protostellar seeds in a constant accretion scenario, the subsequent pre-main sequence evolution, and their pulsation properties. The results are then compared with observations to constrain the input physics. Results: We discuss 16 formerly unknown pulsating pre-main sequence stars and candidates that are either of SPB, $δ$ Scuti,$γ$ Doradus or $δ$ Scuti - $γ$ Doradus hybrid type. We find that evolutionary tracks with a mass accretion rate of $5\times10^{-6} M_\odot/{\rm yr}$ and fraction of injected accretion energy of $β=0.1$ provide the best results in enveloping the spectroscopic parameters of pre-main sequence stars in a constant accretion scenario. The calculated instability regions constrain the atmospheric boundary conditions to Eddington Gray atmospheres; we discuss the future potential for additional constraints by instability regions that are dependent on radial order. We present a possible candidate for pulsations in M-type young stellar objects. △ Less

Submitted 15 July, 2021; originally announced July 2021.

Comments: 21 pages plus 12 for the appendix, 17 figures plus 27 for the appendix, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 654, A36 (2021)

A Radial Velocity Search for Binary RR Lyrae Variables

Authors: Thomas G. Barnes III, Elisabeth Guggenberger, Katrien Kolenberg

Abstract: We report 272 radial velocities for 19 RR Lyrae variables. For most of the stars we have radial velocities for the complete pulsation cycle. These data are used to determine robust center--of--mass radial velocities that have been compared to values from the literature in a search for evidence of binary systems. Center--of--mass velocities were determined for each star using Fourier Series and Tem… ▽ More We report 272 radial velocities for 19 RR Lyrae variables. For most of the stars we have radial velocities for the complete pulsation cycle. These data are used to determine robust center--of--mass radial velocities that have been compared to values from the literature in a search for evidence of binary systems. Center--of--mass velocities were determined for each star using Fourier Series and Template fits to the radial velocities. Our center--of--mass velocities have uncertainties from $\pm0.16$ km s$^{-1}$ to $\pm$2.5 km s$^{-1}$, with a mean uncertainty of $\pm$0.92 km s$^{-1}$. We combined our center--of--mass velocities with values from the literature to look for deviations from the mean center--of--mass velocity of each star. Fifteen RR Lyrae show no evidence of binary motion (BK And, CI And, Z CVn, DM Cyg, BK Dra, RR Gem, XX Hya, SZ Leo, BX Leo, TT Lyn, CN Lyr, TU Per, U Tri, RV UMa, and AV Vir). In most cases this conclusion is reached due to the sporadic sampling of the center--of--mass velocities over time. Three RR Lyrae show suspicious variation in the center--of--mass velocities that may indicate binary motion but do not prove it (SS Leo, ST Leo, and AO Peg). TU UMa was observed by us near a predicted periastron passage (at 0.14 in orbital phase) but the absence of additional center--of--mass velocities near periastron make the binary detection, based on radial velocities alone, uncertain. Two stars in our sample show $Hγ$ emission in phases 0.9--1.0: SS Leo and TU UMa. △ Less

Submitted 9 June, 2021; originally announced June 2021.

Comments: 5 tables, 21 figures

Fast and Accurate Emulation of the SDO/HMI Stokes Inversion with Uncertainty Quantification

Authors: Richard E. L. Higgins, David F. Fouhey, Dichang Zhang, Spiro K. Antiochos, Graham Barnes, J. Todd Hoeksema, K. D. Leka, Yang Liu, Peter W. Schuck, Tamas I. Gombosi

Abstract: The Helioseismic and Magnetic Imager (HMI) onboard NASA's Solar Dynamics Observatory (SDO) produces estimates of the photospheric magnetic field which are a critical input to many space weather modelling and forecasting systems. The magnetogram products produced by HMI and its analysis pipeline are the result of a per-pixel optimization that estimates solar atmospheric parameters and minimizes dis… ▽ More The Helioseismic and Magnetic Imager (HMI) onboard NASA's Solar Dynamics Observatory (SDO) produces estimates of the photospheric magnetic field which are a critical input to many space weather modelling and forecasting systems. The magnetogram products produced by HMI and its analysis pipeline are the result of a per-pixel optimization that estimates solar atmospheric parameters and minimizes disagreement between a synthesized and observed Stokes vector. In this paper, we introduce a deep learning-based approach that can emulate the existing HMI pipeline results two orders of magnitude faster than the current pipeline algorithms. Our system is a U-Net trained on input Stokes vectors and their accompanying optimization-based VFISV inversions. We demonstrate that our system, once trained, can produce high-fidelity estimates of the magnetic field and kinematic and thermodynamic parameters while also producing meaningful confidence intervals. We additionally show that despite penalizing only per-pixel loss terms, our system is able to faithfully reproduce known systematic oscillations in full-disk statistics produced by the pipeline. This emulation system could serve as an initialization for the full Stokes inversion or as an ultra-fast proxy inversion. This work is part of the NASA Heliophysics DRIVE Science Center (SOLSTICE) at the University of Michigan, under grant NASA 80NSSC20K0600E, and has been open sourced. △ Less

Submitted 27 August, 2021; v1 submitted 31 March, 2021; originally announced March 2021.

On measuring divergence for magnetic field modeling

Authors: S. A. Gilchrist, K. D. Leka, G. Barnes, M. S. Wheatland, M. L. DeRosa

Abstract: A physical magnetic field has a divergence of zero. Numerical error in constructing a model field and computing the divergence, however, introduces a finite divergence into these calculations. A popular metric for measuring divergence is the average fractional flux $\langle |f_{i}| \rangle$. We show that $\langle |f_{i}| \rangle$ scales with the size of the computational mesh, and may be a poor me… ▽ More A physical magnetic field has a divergence of zero. Numerical error in constructing a model field and computing the divergence, however, introduces a finite divergence into these calculations. A popular metric for measuring divergence is the average fractional flux $\langle |f_{i}| \rangle$. We show that $\langle |f_{i}| \rangle$ scales with the size of the computational mesh, and may be a poor measure of divergence because it becomes arbitrarily small for increasing mesh resolution, without the divergence actually decreasing. We define a modified version of this metric that does not scale with mesh size. We apply the new metric to the results of DeRosa et al. (2015), who measured $\langle |f_{i}| \rangle$ for a series of Nonlinear Force-Free Field (NLFFF) models of the coronal magnetic field based on solar boundary data binned at different spatial resolutions. We compute a number of divergence metrics for the DeRosa et al. (2015) data and analyze the effect of spatial resolution on these metrics using a non-parametric method. We find that some of the trends reported by DeRosa et al. (2015) are due to the intrinsic scaling of $\langle |f_{i}| \rangle$. We also find that different metrics give different results for the same data set and therefore there is value in measuring divergence via several metrics. △ Less

Submitted 20 August, 2020; originally announced August 2020.

Comments: Accepted for publication in ApJ

Journal ref: ApJ 900 136 (2020)

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

Authors: Sung-Hong Park, K. D. Leka, Kanya Kusano, Jesse Andries, Graham Barnes, Suzy Bingham, D. Shaun Bloomfield, Aoife E. McCloskey, Veronique Delouille, David Falconer, Peter T. Gallagher, Manolis K. Georgoulis, Yuki Kubo, Kangjin Lee, Sangwoo Lee, Vasily Lobzin, JunChul Mun, Sophie A. Murray, Tarek A. M. Hamad Nageem, Rami Qahwaji, Michael Sharpe, Rob A. Steenburgh, Graham Steward, Michael Terkildsen

Abstract: 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 per… ▽ More 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. △ Less

Submitted 21 January, 2020; v1 submitted 8 January, 2020; originally announced January 2020.

Comments: 33 pages, 13 figures, accepted for publication in ApJ

A Comparison of Flare Forecasting Methods. III. Systematic Behaviors of Operational Solar Flare Forecasting Systems

Authors: K. D. Leka, Sung-Hong Park, Kanya Kusano, Jesse Andries, Graham Barnes, Suzy Bingham, D. Shaun Bloomfield, Aoife E. McCloskey, Veronique Delouille, David Falconer, Peter T. Gallagher, Manolis K. Georgoulis, Yuki Kubo, Kangjin Lee, Sangwoo Lee, Vasily Lobzin, JunChul Mun, Sophie A. Murray, Tarek A. M. Hamad Nageem, Rami Qahwaji, Michael Sharpe, Rob Steenburgh, Graham Steward, Michael Terkildsen

Abstract: A workshop was recently held at Nagoya University (31 October - 02 November 2017), sponsored by the Center for International Collaborative Research, at the Institute for Space-Earth Environmental Research, Nagoya University, Japan, to quantitatively compare the performance of today's operational solar flare forecasting facilities. Building upon Paper I of this series (Barnes et al. 2016), in Paper… ▽ More A workshop was recently held at Nagoya University (31 October - 02 November 2017), sponsored by the Center for International Collaborative Research, at the Institute for Space-Earth Environmental Research, Nagoya University, Japan, to quantitatively compare the performance of today's operational solar flare forecasting facilities. Building upon Paper I of this series (Barnes et al. 2016), in Paper II (Leka et al. 2019) we described the participating methods for this latest comparison effort, the evaluation methodology, and presented quantitative comparisons. In this paper we focus on the behavior and performance of the methods when evaluated in the context of broad implementation differences. Acknowledging the short testing interval available and the small number of methods available, we do find that forecast performance: 1) appears to improve by including persistence or prior flare activity, region evolution, and a human "forecaster in the loop"; 2) is hurt by restricting data to disk-center observations; 3) may benefit from long-term statistics, but mostly when then combined with modern data sources and statistical approaches. These trends are arguably weak and must be viewed with numerous caveats, as discussed both here and in Paper II. Following this present work, we present in Paper IV a novel analysis method to evaluate temporal patterns of forecasting errors of both types (i.e., misses and false alarms; Park et al. 2019). Hence, most importantly, with this series of papers we demonstrate the techniques for facilitating comparisons in the interest of establishing performance-positive methodologies. △ Less

Submitted 5 July, 2019; originally announced July 2019.

Comments: 23 pages, 6 figures, accepted for publication in The Astrophysical Journal

A Comparison of Flare Forecasting Methods. II. Benchmarks, Metrics and Performance Results for Operational Solar Flare Forecasting Systems

Authors: K. D. Leka, Sung-Hong Park, Kanya Kusano, Jesse Andries, Graham Barnes, Suzy Bingham, D. Shaun Bloomfield, Aoife E. McCloskey, Veronique Delouille, David Falconer, Peter T. Gallagher, Manolis K. Georgoulis, Yuki Kubo, Kangjin Lee, Sangwoo Lee, Vasily Lobzin, JunChul Mun, Sophie A. Murray, Tarek A. M. Hamad Nageem, Rami Qahwaji, Michael Sharpe, Rob Steenburgh, Graham Steward, Michael Terkildsen

Abstract: Solar flares are extremely energetic phenomena in our Solar System. Their impulsive, often drastic radiative increases, in particular at short wavelengths, bring immediate impacts that motivate solar physics and space weather research to understand solar flares to the point of being able to forecast them. As data and algorithms improve dramatically, questions must be asked concerning how well the… ▽ More Solar flares are extremely energetic phenomena in our Solar System. Their impulsive, often drastic radiative increases, in particular at short wavelengths, bring immediate impacts that motivate solar physics and space weather research to understand solar flares to the point of being able to forecast them. As data and algorithms improve dramatically, questions must be asked concerning how well the forecasting performs; crucially, we must ask how to rigorously measure performance in order to critically gauge any improvements. Building upon earlier-developed methodology (Barnes et al, 2016, Paper I), international representatives of regional warning centers and research facilities assembled in 2017 at the Institute for Space-Earth Environmental Research, Nagoya University, Japan to - for the first time - directly compare the performance of operational solar flare forecasting methods. Multiple quantitative evaluation metrics are employed, with focus and discussion on evaluation methodologies given the restrictions of operational forecasting. Numerous methods performed consistently above the "no skill" level, although which method scored top marks is decisively a function of flare event definition and the metric used; there was no single winner. Following in this paper series we ask why the performances differ by examining implementation details (Leka et al. 2019, Paper III), and then we present a novel analysis method to evaluate temporal patterns of forecasting errors in (Park et al. 2019, Paper IV). With these works, this team presents a well-defined and robust methodology for evaluating solar flare forecasting methods in both research and operational frameworks, and today's performance benchmarks against which improvements and new methods may be compared. △ Less

Submitted 5 July, 2019; originally announced July 2019.

Comments: 26 pages, 5 figures, accepted for publication in the Astrophysical Journal Supplement Series

The effect of metallicity on Cepheid Period-Luminosity relations from a Baade-Wesselink analysis of Cepheids in the Milky Way and Magellanic Clouds

Authors: W. Gieren, J. Storm, P. Konorski, M. Górski, B. Pilecki, I. Thompson, G. Pietrzyński, D. Graczyk, T. G. Barnes, P. Fouqué, N. Nardetto, A. Gallenne, P. Karczmarek, K. Suchomska, P. Wielgórski, M. Taormina, B. Zgirski

Abstract: The extragalactic distance scale builds on the Cepheid period-luminosity (PL) relation. In this paper, we want to carry out a strictly differential comparison of the absolute PL relations obeyed by classical Cepheids in the Milky Way (MW), LMC and SMC galaxies. Taking advantage of the substantial metallicity difference among the Cepheid populations in these three galaxies, we want to establish a p… ▽ More The extragalactic distance scale builds on the Cepheid period-luminosity (PL) relation. In this paper, we want to carry out a strictly differential comparison of the absolute PL relations obeyed by classical Cepheids in the Milky Way (MW), LMC and SMC galaxies. Taking advantage of the substantial metallicity difference among the Cepheid populations in these three galaxies, we want to establish a possible systematic trend of the PL relation absolute zero point as a function of metallicity, and determine the size of such an effect in optical and near-infrared photometric bands. We are using the IRSB Baade-Wesselink type method as calibrated by Storm et al. to determine individual distances to the Cepheids in our samples in MW, LMC and SMC. For our analysis, we use a greatly enhanced sample of Cepheids in the SMC (31 stars) as compared to the small sample (5 stars) available in our previous work. We use the distances to determine absolute Cepheid PL relations in optical and near-infrared bands in each of the three galaxies.} {Our distance analysis of 31 SMC Cepheids with periods from 4-69 days yields tight PL relations in all studied bands, with slopes consistent with the corresponding LMC and MW relations. Adopting the very accurately determined LMC slopes for the optical and near-infrared bands, we determine the zero point offsets between the corresponding absolute PL relations in the 3 galaxies. We find that in all bands the metal-poor SMC Cepheids are intrinsically fainter than their more metal-rich counterparts in the LMC and MW. In the $K$ band the metallicity effect is $-0.23\pm0.06$~mag/dex while in the $V,(V-I)$ Wesenheit index it is slightly stronger, $-0.34\pm0.06$~mag/dex. We find some evidence that the PL relation zero point-metallicity relation might be nonlinear, becoming steeper for lower metallicities. △ Less

Submitted 11 September, 2018; originally announced September 2018.

Comments: 20 pages, accepted for publication by A&A

Journal ref: A&A 620, A99 (2018)

Immersive Virtual Reality Experiences for All-Sky Data

Authors: C. J. Fluke, D. G. Barnes

Abstract: Spherical coordinate systems, which are ubiquitous in astronomy, cannot be shown without distortion on flat, two-dimensional surfaces. This poses challenges for the two complementary phases of visual exploration -- making discoveries in data by looking for relationships, patterns or anomalies -- and publication -- where the results of an exploration are made available for scientific scrutiny or co… ▽ More Spherical coordinate systems, which are ubiquitous in astronomy, cannot be shown without distortion on flat, two-dimensional surfaces. This poses challenges for the two complementary phases of visual exploration -- making discoveries in data by looking for relationships, patterns or anomalies -- and publication -- where the results of an exploration are made available for scientific scrutiny or communication. This is a long-standing problem, and many practical solutions have been developed. Our allskyVR approach provides a workflow for experimentation with commodity virtual reality head-mounted displays. Using the free, open source S2PLOT programming library, and the A-Frame WebVR browser-based framework, we provide a straightforward way to visualise all-sky catalogues on a user-centred, virtual celestial sphere. The allskyVR distribution contains both a quickstart option, complete with a gaze-based menu system, and a fully customisable mode for those who need more control of the immersive experience. The software is available for download from: https://github.com/cfluke/allskyVR △ Less

Submitted 8 May, 2018; originally announced May 2018.

Comments: 11 pages, 3 figures, accepted for publication in Publications of the Astronomical Society of Australia

The NWRA Classification Infrastructure: Description and Extension to the Discriminant Analysis Flare Forecasting System (DAFFS)

Authors: K. D. Leka, Graham Barnes, Eric L. Wagner

Abstract: A classification infrastructure built upon Discriminant Analysis has been developed at NorthWest Research Associates for examining the statistical differences between samples of two known populations. Originating to examine the physical differences between flare-quiet and flare-imminent solar active regions, we describe herein some details of the infrastructure including: parametrization of large… ▽ More A classification infrastructure built upon Discriminant Analysis has been developed at NorthWest Research Associates for examining the statistical differences between samples of two known populations. Originating to examine the physical differences between flare-quiet and flare-imminent solar active regions, we describe herein some details of the infrastructure including: parametrization of large datasets, schemes for handling "null" and "bad" data in multi-parameter analysis, application of non-parametric multi-dimensional Discriminant Analysis, an extension through Bayes' theorem to probabilistic classification, and methods invoked for evaluating classifier success. The classifier infrastructure is applicable to a wide range of scientific questions in solar physics. We demonstrate its application to the question of distinguishing flare-imminent from flare-quiet solar active regions, updating results from the original publications that were based on different data and much smaller sample sizes. Finally, as a demonstration of "Research to Operations" efforts in the space-weather forecasting context, we present the Discriminant Analysis Flare Forecasting System (DAFFS), a near-real-time operationally-running solar flare forecasting tool that was developed from the research-directed infrastructure. △ Less

Submitted 19 February, 2018; originally announced February 2018.

Comments: J. Space Weather Space Climate: Accepted / in press; access supplementary materials through journal; some figures are less than full resolution for arXiv

Does Nearby Open Flux Affect the Eruptivity of Solar Active Regions?

Authors: Marc L. DeRosa, Graham Barnes

Abstract: The most energetic solar flares are typically associated with the ejection of a cloud of coronal material into the heliosphere in the form of a coronal mass ejection (CME). However, there exist large flares which are not accompanied by a CME. The existence of these non-eruptive flares raises the question of whether such flares suffer from a lack of access to nearby open fields in the vicinity abov… ▽ More The most energetic solar flares are typically associated with the ejection of a cloud of coronal material into the heliosphere in the form of a coronal mass ejection (CME). However, there exist large flares which are not accompanied by a CME. The existence of these non-eruptive flares raises the question of whether such flares suffer from a lack of access to nearby open fields in the vicinity above the flare (reconnection) site. In this study, we use a sample of 56 flares from Sunspot Cycles 23 and 24 to test whether active regions that produce eruptive X-class flares are preferentially located near coronal magnetic field domains that are open to the heliosphere, as inferred from a potential field source surface model. The study shows that X-class flares having access to open fields are eruptive at a higher rate than those for which access is lacking. The significance of this result should be moderated due to the small number of non-eruptive X-class flares in the sample, based on the associated Bayes factor. △ Less

Submitted 30 May, 2018; v1 submitted 4 February, 2018; originally announced February 2018.

Comments: ApJ accepted

A cautionary tale of interpreting O-C diagrams: period instability in a classical RR Lyr Star Z CVn mimicking as a distant companion

Authors: M. Skarka, J. Liška, R. Dřevěný, E. Guggenberger, Á. Sódor, T. G. Barnes, K. Kolenberg

Abstract: We present a comprehensive study of Z CVn, an RR Lyrae star that shows long-term cyclic variations of its pulsation period. A possible explanation suggested from the shape of the O-C diagram is the light travel-time effect, which we thoroughly examine. We used original photometric and spectroscopic measurements and investigated the period evolution using available maximum times spanning more than… ▽ More We present a comprehensive study of Z CVn, an RR Lyrae star that shows long-term cyclic variations of its pulsation period. A possible explanation suggested from the shape of the O-C diagram is the light travel-time effect, which we thoroughly examine. We used original photometric and spectroscopic measurements and investigated the period evolution using available maximum times spanning more than one century. If the binary hypothesis is valid, Z CVn orbits around a black hole with minimal mass of $56.5$ $\mathfrak{M}_{\odot}$ on a very wide ($P_{\rm orbit}=78.3$ years) and eccentric orbit ($e=0.63$). We discuss the probability of a formation of a black hole-RR Lyrae pair and, although we found it possible, there is no observational evidence of the black hole in the direction to Z CVn. However, the main objection against the binary hypothesis is the comparison of the systemic radial velocity curve model and spectroscopic observations that clearly show that Z CVn cannot be bound in such a binary. Therefore, the variations of pulsation period are likely intrinsic to the star. This finding represents a discovery/confirmation of a new type of cyclic period changes in RR Lyrae stars. By the analysis of our photometric data, we found that the Blazhko modulation with period of 22.931 d is strongly dominant in amplitude. The strength of the phase modulation varies and is currently almost undetectable. We also estimated photometric physical parameters of Z CVn and investigated their variations during the Blazhko cycle using the Inverse Baade-Wesselink method. △ Less

Submitted 18 October, 2017; originally announced October 2017.

Comments: 15 pages, 8 tables, 9 figures, accepted for publication in MNRAS

Real-time colouring and filtering with graphics shaders

Authors: Dany Vohl, Christopher J. Fluke, David G. Barnes, Amr H. Hassan

Abstract: Despite the popularity of the Graphics Processing Unit (GPU) for general purpose computing, one should not forget about the practicality of the GPU for fast scientific visualisation. As astronomers have increasing access to three dimensional (3D) data from instruments and facilities like integral field units and radio interferometers, visualisation techniques such as volume rendering offer means t… ▽ More Despite the popularity of the Graphics Processing Unit (GPU) for general purpose computing, one should not forget about the practicality of the GPU for fast scientific visualisation. As astronomers have increasing access to three dimensional (3D) data from instruments and facilities like integral field units and radio interferometers, visualisation techniques such as volume rendering offer means to quickly explore spectral cubes as a whole. As most 3D visualisation techniques have been developed in fields of research like medical imaging and fluid dynamics, many transfer functions are not optimal for astronomical data. We demonstrate how transfer functions and graphics shaders can be exploited to provide new astronomy-specific explorative colouring methods. We present 12 shaders, including four novel transfer functions specifically designed to produce intuitive and informative 3D visualisations of spectral cube data. We compare their utility to classic colour mapping. The remaining shaders highlight how common computation like filtering, smoothing and line ratio algorithms can be integrated as part of the graphics pipeline. We discuss how this can be achieved by utilising the parallelism of modern GPUs along with a shading language, letting astronomers apply these new techniques at interactive frame rates. All shaders investigated in this work are included in the open source software shwirl (Vohl, 2017). △ Less

Submitted 3 July, 2017; originally announced July 2017.

Comments: 24 pages, 19 figures, 14 algorithms, 1 table. Accepted for publication by Monthly Notices of the Royal Astronomical Society

A comprehensive study of young B stars in NGC 2264: I. Space photometry and asteroseismology

Authors: K. Zwintz, E. Moravveji, P. I. Papics, A. Tkachenko, N. Przybilla, M. -F. Nieva, R. Kuschnig, V. Antoci, D. Lorenz, N. Themessl, L. Fossati, T. G. Barnes

Abstract: Space photometric time series of the most massive members of the young open cluster NGC 2264 allow us to study their different sources of variability down to the millimagnitude level and permits a search for Slowly Pulsating B (SPB) type pulsation among objects that are only a few million years old. Our goal is to conduct a homogeneous study of young B type stars in the cluster NGC 2264 using phot… ▽ More Space photometric time series of the most massive members of the young open cluster NGC 2264 allow us to study their different sources of variability down to the millimagnitude level and permits a search for Slowly Pulsating B (SPB) type pulsation among objects that are only a few million years old. Our goal is to conduct a homogeneous study of young B type stars in the cluster NGC 2264 using photometric time series from space in combination with high-resolution spectroscopy and spectropolarimetry obtained from the ground. The latter will be presented in a separate follow-up article. We performed frequency analyses for eleven B stars in the field of the young cluster NGC 2264 using photometric time series from the MOST, CoRoT and Spitzer space telescopes and the routines Period04 and SigSpec. We employ the MESA stellar evolution code in combination with the oscillation code GYRE to identify the pulsation modes for two SPB stars which exhibit short period spacing series. From our analysis we identify four objects that show SPB pulsations, five stars that show rotational modulation of their light curves caused by spots, one star that is identified to be a binary, and one object in the field of the cluster that is found to be a non-member Be star. In two SPB stars we detect a number of regularly spaced pulsation modes that are compatible with being members of a g mode period series. Despite NGC 2264's young age, our analysis illustrates that its B type members have already arrived on the zero-age main sequence (ZAMS). Our asteroseismic analysis yields masses between 4 and 6 Msun and ages between 1 and 6 million years, which agree well to the overall cluster age. △ Less

Submitted 19 March, 2017; originally announced March 2017.

Comments: 19 pages, 20 figures, accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 601, A101 (2017)

Evaluating (and Improving) Estimates of the Solar Radial Magnetic Field Component from Line-of-Sight Magnetograms

Authors: K. D. Leka, G. Barnes, E. L. Wagner

Abstract: Although for many solar physics problems the desirable or meaningful boundary is the radial component of the magnetic field $B_{\rm r}$, the most readily available measurement is the component of the magnetic field along the line-of-sight to the observer, $B_{\rm los}$. As this component is only equal to the radial component where the viewing angle is exactly zero, some approximation is required t… ▽ More Although for many solar physics problems the desirable or meaningful boundary is the radial component of the magnetic field $B_{\rm r}$, the most readily available measurement is the component of the magnetic field along the line-of-sight to the observer, $B_{\rm los}$. As this component is only equal to the radial component where the viewing angle is exactly zero, some approximation is required to estimate $B_{\rm r}$ at all other observed locations. In this study, a common approximation known as the "$μ$-correction", which assumes all photospheric field to be radial, is compared to a method which invokes computing a potential field that matches the observed $B_{\rm los}$, from which the potential field radial component, $B_{\rm r}^{\rm pot}$ is recovered. We demonstrate that in regions that are truly dominated by radially-oriented field at the resolution of the data employed, the $μ$-correction performs acceptably if not better than the potential-field approach. However, it is also shown that for any solar structure which includes horizontal fields, i.e. active regions, the potential-field method better recovers both the strength of the radial field and the location of magnetic neutral line. △ Less

Submitted 17 January, 2017; originally announced January 2017.

Comments: Accepted for publication in Solar Physics

Collaborative visual analytics of radio surveys in the Big Data era

Authors: Dany Vohl, Christopher J. Fluke, Amr H. Hassan, David G. Barnes, Virginia A. Kilborn

Abstract: Radio survey datasets comprise an increasing number of individual observations stored as sets of multidimensional data. In large survey projects, astronomers commonly face limitations regarding: 1) interactive visual analytics of sufficiently large subsets of data; 2) synchronous and asynchronous collaboration; and 3) documentation of the discovery workflow. To support collaborative data inquiry,… ▽ More Radio survey datasets comprise an increasing number of individual observations stored as sets of multidimensional data. In large survey projects, astronomers commonly face limitations regarding: 1) interactive visual analytics of sufficiently large subsets of data; 2) synchronous and asynchronous collaboration; and 3) documentation of the discovery workflow. To support collaborative data inquiry, we present encube, a large-scale comparative visual analytics framework. Encube can utilise advanced visualization environments such as the CAVE2 (a hybrid 2D and 3D virtual reality environment powered with a 100 Tflop/s GPU-based supercomputer and 84 million pixels) for collaborative analysis of large subsets of data from radio surveys. It can also run on standard desktops, providing a capable visual analytics experience across the display ecology. Encube is composed of four primary units enabling compute-intensive processing, advanced visualisation, dynamic interaction, parallel data query, along with data management. Its modularity will make it simple to incorporate astronomical analysis packages and Virtual Observatory capabilities developed within our community. We discuss how encube builds a bridge between high-end display systems (such as CAVE2) and the classical desktop, preserving all traces of the work completed on either platform -- allowing the research process to continue wherever you are. △ Less

Submitted 2 December, 2016; originally announced December 2016.

Comments: Accepted for Proceedings of IAU Symposium 325

An interactive, comparative and quantitative 3D visualization system for large-scale spectral-cube surveys using CAVE2

Authors: Dany Vohl, Christopher J. Fluke, Amr H. Hassan, David G. Barnes

Abstract: As the quantity and resolution of spectral-cubes from optical/infrared and radio surveys increase, desktop-based visualization and analysis solutions must adapt and evolve. Novel immersive 3D environments such as the CAVE2 at Monash University can overcome personal computer's visualization limitations. CAVE2 is part advanced 2D/3D visualization space (80 stereo-capable screens providing a total of… ▽ More As the quantity and resolution of spectral-cubes from optical/infrared and radio surveys increase, desktop-based visualization and analysis solutions must adapt and evolve. Novel immersive 3D environments such as the CAVE2 at Monash University can overcome personal computer's visualization limitations. CAVE2 is part advanced 2D/3D visualization space (80 stereo-capable screens providing a total of 84 million pixels) and part supercomputer ($\sim100$ TFLOPS of integrated GPU-based processing power). We present a novel visualization system enabling simultaneous 3D comparative visualization of $\sim100$ spectral-cubes. With CAVE2 augmented by our newly developed web-based controller interface, astronomers can easily organise spectral-cubes on the different display panels, apply real-time transforms to one or many spectral cubes, and request quantitative information about the displayed data. We also discuss how such a solution can help accelerate the discovery rate in varied research scenarios. △ Less

Submitted 3 October, 2016; originally announced October 2016.

Comments: 4 pages, 2 figures, to appear in Proceedings of ADASS XXV (Sydney, Australia; 2015), edited by N. P. F. Lorente & K. Shortridge (San Francisco: ASP), ASP Conf. Series

Large-scale comparative visualisation of sets of multidimensional data

Authors: Dany Vohl, David G. Barnes, Christopher J. Fluke, Govinda Poudel, Nellie Georgiou-Karistianis, Amr H. Hassan, Yuri Benovitski, Tsz Ho Wong, Owen Kaluza, Toan D. Nguyen, C. Paul Bonnington

Abstract: We present encube $-$ a qualitative, quantitative and comparative visualisation and analysis system, with application to high-resolution, immersive three-dimensional environments and desktop displays. encube extends previous comparative visualisation systems by considering: 1) the integration of comparative visualisation and analysis into a unified system; 2) the documentation of the discovery pro… ▽ More We present encube $-$ a qualitative, quantitative and comparative visualisation and analysis system, with application to high-resolution, immersive three-dimensional environments and desktop displays. encube extends previous comparative visualisation systems by considering: 1) the integration of comparative visualisation and analysis into a unified system; 2) the documentation of the discovery process; and 3) an approach that enables scientists to continue the research process once back at their desktop. Our solution enables tablets, smartphones or laptops to be used as interaction units for manipulating, organising, and querying data. We highlight the modularity of encube, allowing additional functionalities to be included as required. Additionally, our approach supports a high level of collaboration within the physical environment. We show how our implementation of encube operates in a large-scale, hybrid visualisation and supercomputing environment using the CAVE2 at Monash University, and on a local desktop, making it a versatile solution. We discuss how our approach can help accelerate the discovery rate in a variety of research scenarios. △ Less

Submitted 3 October, 2016; originally announced October 2016.

Comments: 26 pages, 11 figures, 5 tables. Accepted for publication in PeerJ Computer Science

A Fixed-point Scheme for the Numerical Construction of Magnetohydrostatic Atmospheres in Three Dimensions

Authors: S. A. Gilchrist, D. C. Braun, G. Barnes

Abstract: Magnetohydrostatic models of the solar atmosphere are often based on idealized analytic solutions because the underlying equations are too difficult to solve in full generality. Numerical approaches, too, are often limited in scope and have tended to focus on the two-dimensional problem. In this article we develop a numerical method for solving the nonlinear magnetohydrostatic equations in three d… ▽ More Magnetohydrostatic models of the solar atmosphere are often based on idealized analytic solutions because the underlying equations are too difficult to solve in full generality. Numerical approaches, too, are often limited in scope and have tended to focus on the two-dimensional problem. In this article we develop a numerical method for solving the nonlinear magnetohydrostatic equations in three dimensions. Our method is a fixed-point iteration scheme that extends the method of Grad and Rubin (Proc. 2nd Int. Conf. on Peaceful Uses of Atomic Energy 31, 190, 1958) to include a finite gravity force. We apply the method to a test case to demonstrate the method in general and our implementation in code in particular. △ Less

Submitted 27 October, 2016; v1 submitted 2 September, 2016; originally announced September 2016.

Comments: Accepted for publication in Solar Physics

A Comparison of Flare Forecasting Methods, I: Results from the "All-Clear" Workshop

Authors: G. Barnes, K. D. Leka, C. J. Schrijver, T. Colak, R. Qahwaji, O. W. Ashamari, Y. Yuan, J. Zhang, R. T. J. McAteer, D. S. Bloomfield, P. A. Higgins, P. T. Gallagher, D. A. Falconer, M. K. Georgoulis, M. S. Wheatland, C. Balch, T. Dunn, E. L. Wagner

Abstract: Solar flares produce radiation which can have an almost immediate effect on the near-Earth environment, making it crucial to forecast flares in order to mitigate their negative effects. The number of published approaches to flare forecasting using photospheric magnetic field observations has proliferated, with varying claims about how well each works. Because of the different analysis techniques a… ▽ More Solar flares produce radiation which can have an almost immediate effect on the near-Earth environment, making it crucial to forecast flares in order to mitigate their negative effects. The number of published approaches to flare forecasting using photospheric magnetic field observations has proliferated, with varying claims about how well each works. Because of the different analysis techniques and data sets used, it is essentially impossible to compare the results from the literature. This problem is exacerbated by the low event rates of large solar flares. The challenges of forecasting rare events have long been recognized in the meteorology community, but have yet to be fully acknowledged by the space weather community. During the interagency workshop on "all clear" forecasts held in Boulder, CO in 2009, the performance of a number of existing algorithms was compared on common data sets, specifically line-of-sight magnetic field and continuum intensity images from MDI, with consistent definitions of what constitutes an event. We demonstrate the importance of making such systematic comparisons, and of using standard verification statistics to determine what constitutes a good prediction scheme. When a comparison was made in this fashion, no one method clearly outperformed all others, which may in part be due to the strong correlations among the parameters used by different methods to characterize an active region. For M-class flares and above, the set of methods tends towards a weakly positive skill score (as measured with several distinct metrics), with no participating method proving substantially better than climatological forecasts. △ Less

Submitted 22 August, 2016; originally announced August 2016.

Comments: Accepted for publication in ApJ

The Ultimate Display

Authors: C. J. Fluke, D. G. Barnes

Abstract: Astronomical images and datasets are increasingly high-resolution and multi-dimensional. The vast majority of astronomers perform all of their visualisation and analysis tasks on low-resolution, two-dimensional desktop monitors. If there were no technological barriers to designing the ultimate stereoscopic display for astronomy, what would it look like? What capabilities would we require of our co… ▽ More Astronomical images and datasets are increasingly high-resolution and multi-dimensional. The vast majority of astronomers perform all of their visualisation and analysis tasks on low-resolution, two-dimensional desktop monitors. If there were no technological barriers to designing the ultimate stereoscopic display for astronomy, what would it look like? What capabilities would we require of our compute hardware to drive it? And are existing technologies even close to providing a true 3D experience that is compatible with the depth resolution of human stereoscopic vision? We consider the CAVE2 (an 80 Megapixel, hybrid 2D and 3D virtual reality environment directly integrated with a 100 Tflop/s GPU-powered supercomputer) and the Oculus Rift (a low- cost, head-mounted display) as examples at opposite financial ends of the immersive display spectrum. △ Less

Submitted 13 January, 2016; originally announced January 2016.

Comments: 4 pages, 2 figures, to appear in Proceedings of ADASS XXV (Sydney, Australia; 2015), edited by N. P. F. Lorente & K. Shortridge (San Francisco: ASP), ASP Conf. Series

New systemic radial velocities of suspected RR Lyrae binary stars

Authors: Elisabeth Guggenberger, Thomas G. Barnes, Katrien Kolenberg

Abstract: Among the tens of thousands of known RR Lyrae stars there are only a handful that show indications of possible binarity. The question why this is the case is still unsolved, and has recently sparked several studies dedicated to the search for additional RR Lyraes in binary systems. Such systems are particularly valuable because they might allow to constrain the stellar mass. Most of the recent stu… ▽ More Among the tens of thousands of known RR Lyrae stars there are only a handful that show indications of possible binarity. The question why this is the case is still unsolved, and has recently sparked several studies dedicated to the search for additional RR Lyraes in binary systems. Such systems are particularly valuable because they might allow to constrain the stellar mass. Most of the recent studies, however, are based on photometry by finding a light time effect in the timings of maximum light. This approach is a very promising and successful one, but it has a major drawback: by itself, it cannot serve as a definite proof of binarity, because other phenomena such as the Blazhko effect or intrinsic period changes could lead to similar results. Spectroscopic radial velocity measurements, on the other hand, can serve as definite proof of binarity. We have therefore started a project to study spectroscopically RR Lyrae stars that are suspected to be binaries. We have obtained radial velocity (RV) curves with the 2.1m telescope at McDonald observatory. From these we derive systemic RVs which we will compare to previous measurements in order to find changes induced by orbital motions. We also construct templates of the RV curves that can facilitate future studies. We also observed the most promising RR Lyrae binary candidate, TU UMa, as no recent spectroscopic measurements were available. We present a densely covered pulsational RV curve, which will be used to test the predictions of the orbit models that are based on the O-C variations. △ Less

Submitted 2 December, 2015; originally announced December 2015.

Comments: 4 pages, 1 figure, to appear in the No. 105 issue of the Communications from the Konkoly Observatory of the Hungarian Academy of Sciences

Achieving Consistent Doppler Measurements from SDO/HMI Vector Field Inversions

Authors: P. W. Schuck, Spiro Antiochos, K. D. Leka, Graham Barnes

Abstract: NASA's Solar Dynamics Observatory is delivering vector field observations of the full solar disk with unprecedented temporal and spatial resolution; however, the satellite is in a highly inclined geostationary orbit. The relative spacecraft-Sun velocity varies by $\pm3$~km/s over a day which introduces major orbital artifacts in the Helioseismic Magnetic Imager data. We demonstrate that the orbita… ▽ More NASA's Solar Dynamics Observatory is delivering vector field observations of the full solar disk with unprecedented temporal and spatial resolution; however, the satellite is in a highly inclined geostationary orbit. The relative spacecraft-Sun velocity varies by $\pm3$~km/s over a day which introduces major orbital artifacts in the Helioseismic Magnetic Imager data. We demonstrate that the orbital artifacts contaminate all spatial and temporal scales in the data. We describe a newly-developed three stage procedure for mitigating these artifacts in the Doppler data derived from the Milne-Eddington inversions in the HMI Pipeline. This procedure was applied to full disk images of AR11084 to produce consistent Dopplergrams. The data adjustments reduce the power in the orbital artifacts by 31dB. Furthermore, we analyze in detail the corrected images and show that our procedure greatly improve the temporal and spectral properties of the data without adding any new artifacts. We conclude that this new and easily implemented procedure makes a dramatic improvement in the consistency of the HMI data and in its usefulness for precision scientific studies. △ Less

Submitted 20 November, 2015; originally announced November 2015.

Comments: 58 pages, 19 figures, submitted to ApJ

The Influence of Spatial Resolution on Nonlinear Force-Free Modeling

Authors: M. L. DeRosa, M. S. Wheatland, K. D. Leka, G. Barnes, T. Amari, A. Canou, S. A. Gilchrist, J. K. Thalmann, G. Valori, T. Wiegelmann, C. J. Schrijver, A. Malanushenko, X. Sun, S. Régnier

Abstract: The nonlinear force-free field (NLFFF) model is often used to describe the solar coronal magnetic field, however a series of earlier studies revealed difficulties in the numerical solution of the model in application to photospheric boundary data. We investigate the sensitivity of the modeling to the spatial resolution of the boundary data, by applying multiple codes that numerically solve the NLF… ▽ More The nonlinear force-free field (NLFFF) model is often used to describe the solar coronal magnetic field, however a series of earlier studies revealed difficulties in the numerical solution of the model in application to photospheric boundary data. We investigate the sensitivity of the modeling to the spatial resolution of the boundary data, by applying multiple codes that numerically solve the NLFFF model to a sequence of vector magnetogram data at different resolutions, prepared from a single Hinode/SOT-SP scan of NOAA Active Region 10978 on 2007 December 13. We analyze the resulting energies and relative magnetic helicities, employ a Helmholtz decomposition to characterize divergence errors, and quantify changes made by the codes to the vector magnetogram boundary data in order to be compatible with the force-free model. This study shows that NLFFF modeling results depend quantitatively on the spatial resolution of the input boundary data, and that using more highly resolved boundary data yields more self-consistent results. The free energies of the resulting solutions generally trend higher with increasing resolution, while relative magnetic helicity values vary significantly between resolutions for all methods. All methods require changing the horizontal components, and for some methods also the vertical components, of the vector magnetogram boundary field in excess of nominal uncertainties in the data. The solutions produced by the various methods are significantly different at each resolution level. We continue to recommend verifying agreement between the modeled field lines and corresponding coronal loop images before any NLFFF model is used in a scientific setting. △ Less

Submitted 21 August, 2015; originally announced August 2015.

Comments: Accepted to ApJ; comments/corrections to this article are welcome via e-mail, even after publication

A Digital-Receiver for the Murchison Widefield Array

Authors: Thiagaraj Prabu, K. S. Srivani, D. Anish Roshi, P. A. Kamini, S. Madhavi, David Emrich, Brian Crosse, Andrew J. Williams, Mark Waterson, Avinash A. Deshpande, N. Udaya Shankar, Ravi Subrahmanyan, Frank H. Briggs, Robert F. Goeke, Steven J. Tingay, Melanie Johnston-Hollitt, Gopalakrishna M R, Edward H. Morgan, Joseph Pathikulangara, John D. Bunton, Grant Hampson, Christopher Williams, Stephen M. Ord, Randall B. Wayth, Deepak Kumar , et al. (33 additional authors not shown)

Abstract: An FPGA-based digital-receiver has been developed for a low-frequency imaging radio interferometer, the Murchison Widefield Array (MWA). The MWA, located at the Murchison Radio-astronomy Observatory (MRO) in Western Australia, consists of 128 dual-polarized aperture-array elements (tiles) operating between 80 and 300\,MHz, with a total processed bandwidth of 30.72 MHz for each polarization. Radio-… ▽ More An FPGA-based digital-receiver has been developed for a low-frequency imaging radio interferometer, the Murchison Widefield Array (MWA). The MWA, located at the Murchison Radio-astronomy Observatory (MRO) in Western Australia, consists of 128 dual-polarized aperture-array elements (tiles) operating between 80 and 300\,MHz, with a total processed bandwidth of 30.72 MHz for each polarization. Radio-frequency signals from the tiles are amplified and band limited using analog signal conditioning units; sampled and channelized by digital-receivers. The signals from eight tiles are processed by a single digital-receiver, thus requiring 16 digital-receivers for the MWA. The main function of the digital-receivers is to digitize the broad-band signals from each tile, channelize them to form the sky-band, and transport it through optical fibers to a centrally located correlator for further processing. The digital-receiver firmware also implements functions to measure the signal power, perform power equalization across the band, detect interference-like events, and invoke diagnostic modes. The digital-receiver is controlled by high-level programs running on a single-board-computer. This paper presents the digital-receiver design, implementation, current status, and plans for future enhancements. △ Less

Submitted 13 February, 2015; originally announced February 2015.

Comments: 14 pages, 7 figures

An In-Depth Spectroscopic Analysis of RR Lyr Variations over the Pulsation Cycle

Authors: L. Fossati, K. Kolenberg, D. V. Shulyak, A. Elmasli, V. Tsymbal, T. G. Barnes, E. Guggenberger, O. Kochukhov

Abstract: The stellar parameters of RR Lyrae stars vary considerably over a pulsation cycle, and their determination is crucial for stellar modelling. We present a detailed spectroscopic analysis of the pulsating star RR Lyr, the prototype of its class, over a complete pulsation cycle, based on high-resolution spectra collected at the 2.7-m telescope of McDonald Observatory. We used simultaneous photometry… ▽ More The stellar parameters of RR Lyrae stars vary considerably over a pulsation cycle, and their determination is crucial for stellar modelling. We present a detailed spectroscopic analysis of the pulsating star RR Lyr, the prototype of its class, over a complete pulsation cycle, based on high-resolution spectra collected at the 2.7-m telescope of McDonald Observatory. We used simultaneous photometry to determine the accurate pulsation phase of each spectrum and determined the effective temperature, the shape of the depth-dependent microturbulent velocity, and the abundance of several elements, for each phase. The surface gravity was fixed to 2.4. Element abundances resulting from our analysis are stable over the pulsation cycle. However, a variation in ionisation equilibrium is observed around minimum radius. We attribute this mostly to a dynamical acceleration contributing to the surface gravity. Variable turbulent convection on time scales longer than the pulsation cycle has been proposed as a cause for the Blazhko effect. We test this hypothesis to some extent by using the derived variable depth-dependent microturbulent velocity profiles to estimate their effect on the stellar magnitude. These effects turn out to be wavelength-dependent and much smaller than the observed light variations over the Blazhko cycle: if variations in the turbulent motions are entirely responsible for the Blazhko effect, they must surpass the scales covered by the microturbulent velocity. This work demonstrates the possibility of a self-consistent spectroscopic analysis over an entire pulsation cycle using static atmosphere models, provided one takes into account certain features of a rapidly pulsating atmosphere. △ Less

Submitted 16 October, 2014; originally announced October 2014.

Comments: Accepted for publication on MNRAS

Echography of young stars reveals their evolution

Authors: K. Zwintz, L. Fossati, T. Ryabchikova, D. Guenther, C. Aerts, T. G. Barnes, N. Themessl, D. Lorenz, C. Cameron, R. Kuschnig, S. Pollack-Drs, E. Moravveji, A. Baglin, J. M. Matthews, A. F. J. Moffat, E. Poretti, M. Rainer, S. M. Rucinski, D. Sasselov, W. W. Weiss

Abstract: We demonstrate that a seismic analysis of stars in their earliest evolutionary phases is a powerful method to identify young stars and distinguish their evolutionary states. The early star that is born from the gravitational collapse of a molecular cloud reaches at some point sufficient temperature, mass and luminosity to be detected. Accretion stops and the pre-main sequence star that emerges is… ▽ More We demonstrate that a seismic analysis of stars in their earliest evolutionary phases is a powerful method to identify young stars and distinguish their evolutionary states. The early star that is born from the gravitational collapse of a molecular cloud reaches at some point sufficient temperature, mass and luminosity to be detected. Accretion stops and the pre-main sequence star that emerges is nearly fully convective and chemically homogeneous. It will continue to contract gravitationally until the density and temperature in the core are high enough to start nuclear burning of hydrogen. We show that there is a relationship between detected pulsation properties for a sample of young stars and their evolutionary status illustrating the potential of asteroseismology for the early evolutionary phases. △ Less

Submitted 18 July, 2014; originally announced July 2014.

Comments: 23 pages, 9 figures, includes Supplementary Material. Science, published in Science Express on July 3, 2014

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

Authors: J. Todd Hoeksema, Yang Liu, Keiji Hayashi, Xudong Sun, Jesper Schou, Sebastien Couvidat, Aimee Norton, Monica Bobra, Rebecca Centeno, K. D. Leka, Graham Barnes, Michael J. Turmon

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

Submitted 7 April, 2014; originally announced April 2014.

Comments: 42 pages, 19 figures, accepted to Solar Physics

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

Authors: Monica G. Bobra, Xudong Sun, J. Todd Hoeksema, Michael J. Turmon, Yang Liu, Keiji Hayashi, Graham Barnes, K. D. Leka

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

Submitted 7 April, 2014; originally announced April 2014.

Comments: 27 pages, 7 figures. Accepted to Solar Physics

The Helioseismic and Magnetic Imager (HMI) Vector Magnetic Field Pipeline: Optimization of the Spectral Line Inversion Code

Authors: R. Centeno, J. Schou, K. Hayashi, A. Norton, J. T. Hoeksema, Y. Liu, K. D. Leka, G. Barnes

Abstract: The Very Fast Inversion of the Stokes Vector (VFISV) is a Milne-Eddington spectral line inversion code used to determine the magnetic and thermodynamic parameters of the solar photosphere from observations of the Stokes vector in the 6173 A Fe I line by the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO). We report on the modifications made to the original VFISV… ▽ More The Very Fast Inversion of the Stokes Vector (VFISV) is a Milne-Eddington spectral line inversion code used to determine the magnetic and thermodynamic parameters of the solar photosphere from observations of the Stokes vector in the 6173 A Fe I line by the Helioseismic and Magnetic Imager (HMI) onboard the Solar Dynamics Observatory (SDO). We report on the modifications made to the original VFISV inversion code in order to optimize its operation within the HMI data pipeline and provide the smoothest solution in active regions. The changes either sped up the computation or reduced the frequency with which the algorithm failed to converge to a satisfactory solution. Additionally, coding bugs which were detected and fixed in the original VFISV release, are reported here. △ Less

Submitted 14 March, 2014; originally announced March 2014.

Comments: Accepted for publication in Solar Physics

Two spotted and magnetic early B-type stars in the young open cluster NGC2264 discovered by MOST and ESPaDOnS

Authors: L. Fossati, K. Zwintz, N. Castro, N. Langer, D. Lorenz, F. R. N. Schneider, R. Kuschnig, J. M. Matthews, E. Alecian, G. A. Wade, T. G. Barnes, A. A. Thoul

Abstract: Star clusters are known as superb tools for understanding stellar evolution. In a quest for understanding the physical origin of magnetism and chemical peculiarity in about 7% of the massive main-sequence stars, we analysed two of the ten brightest members of the ~10 Myr old Galactic open cluster NGC 2264, the early B-dwarfs HD47887 and HD47777. We find accurate rotation periods of 1.95 and 2.64 d… ▽ More Star clusters are known as superb tools for understanding stellar evolution. In a quest for understanding the physical origin of magnetism and chemical peculiarity in about 7% of the massive main-sequence stars, we analysed two of the ten brightest members of the ~10 Myr old Galactic open cluster NGC 2264, the early B-dwarfs HD47887 and HD47777. We find accurate rotation periods of 1.95 and 2.64 days, respectively, from MOST photometry. We obtained ESPaDOnS spectropolarimetric observations, through which we determined stellar parameters, detailed chemical surface abundances, projected rotational velocities, and the inclination angles of the rotation axis. Because we found only small (<5 km/s) radial velocity variations, most likely caused by spots, we can rule out that HD47887 and HD47777 are close binaries. Finally, using the least-squares deconvolution technique, we found that both stars possess a large-scale magnetic field with an average longitudinal field strength of about 400 G. From a simultaneous fit of the stellar parameters we determine the evolutionary masses of HD47887 and HD47777 to be 9.4+/-0.7 M0 and 7.6+/-0.5 M0. Interestingly, HD47777 shows a remarkable helium underabundance, typical of helium-weak chemically peculiar stars, while the abundances of HD47887 are normal, which might imply that diffusion is operating in the lower mass star but not in the slightly more massive one. Furthermore, we argue that the rather slow rotation, as well as the lack of nitrogen enrichment in both stars, can be consistent with both the fossil and the binary hypothesis for the origin of the magnetic field. However, the presence of two magnetic and apparently single stars near the top of the cluster mass-function may speak in favour of the latter. △ Less

Submitted 21 January, 2014; originally announced January 2014.

Comments: 13 pages, accepted for publication on A&A

A New 1.4GHz Radio Continuum Map of the Sky South of Declination +25 deg

Authors: Mark R. Calabretta, Lister Staveley-Smith, David G. Barnes

Abstract: Archival data from the HI Parkes All-Sky Survey (HIPASS) and the HI Zone of Avoidance (HIZOA) survey have been carefully reprocessed into a new 1.4GHz continuum map of the sky south of declination +25 deg. The wide sky coverage, high sensitivity of 40mK (limited by confusion), resolution of 14'.4 (compared to 51' for the Haslam et al. 408MHz and 35' for the Reich et al. 1.4GHz surveys), and low le… ▽ More Archival data from the HI Parkes All-Sky Survey (HIPASS) and the HI Zone of Avoidance (HIZOA) survey have been carefully reprocessed into a new 1.4GHz continuum map of the sky south of declination +25 deg. The wide sky coverage, high sensitivity of 40mK (limited by confusion), resolution of 14'.4 (compared to 51' for the Haslam et al. 408MHz and 35' for the Reich et al. 1.4GHz surveys), and low level of artefacts makes this map ideal for numerous studies, including: merging into interferometer maps to complete large-scale structures; decomposition of thermal and non-thermal emission components from Galactic and extragalactic sources; and comparison of emission regions with other frequencies. The new map is available for download. △ Less

Submitted 9 October, 2013; originally announced October 2013.

Comments: 30 pages, 17 figures, 4 tables. Publications of the Astronomical Society of Australia, in press. Data products available from http://www.atnf.csiro.au/research/CHIPASS/

The giant lobes of Centaurus A observed at 118 MHz with the Murchison Widefield Array

Authors: B. McKinley, F. Briggs, B. M. Gaensler, I. J. Feain, G. Bernardi, R. B. Wayth, M. Johnston-Hollitt, A. R. Offringa, W. Arcus, D. G. Barnes, J. D. Bowman, J. D. Bunton, R. J. Cappallo, B. E. Corey, A. Deshpande, L. deSouza, D. Emrich, R. Goeke, L. J. Greenhill, B. J. Hazelton, D. Herne, J. N. Hewitt, D. L. Kaplan, J. C. Kasper, B. B. Kincaid , et al. (28 additional authors not shown)

Abstract: We present new wide-field observations of Centaurus A (Cen A) and the surrounding region at 118 MHz with the Murchison Widefield Array (MWA) 32-tile prototype, with which we investigate the spectral-index distribution of Cen A's giant radio lobes. We compare our images to 1.4 GHz maps of Cen A and compute spectral indices using temperature-temperature plots and spectral tomography. We find that th… ▽ More We present new wide-field observations of Centaurus A (Cen A) and the surrounding region at 118 MHz with the Murchison Widefield Array (MWA) 32-tile prototype, with which we investigate the spectral-index distribution of Cen A's giant radio lobes. We compare our images to 1.4 GHz maps of Cen A and compute spectral indices using temperature-temperature plots and spectral tomography. We find that the morphologies at 118 MHz and 1.4 GHz match very closely apart from an extra peak in the southern lobe at 118 MHz, which provides tentative evidence for the existence of a southern counterpart to the northern middle lobe of Cen A. Our spatially-averaged spectral indices for both the northern and southern lobes are consistent with previous analyses, however we find significant spatial variation of the spectra across the extent of each lobe. Both the spectral-index distribution and the morphology at low radio frequencies support a scenario of multiple outbursts of activity from the central engine. Our results are consistent with inverse-Compton modelling of radio and gamma-ray data that supports a value for the lobe age of between 10 and 80 Myr. △ Less

Submitted 4 September, 2013; originally announced September 2013.

Comments: 18 pages, 11 figures, accepted for publication in MNRAS

On the detection and tracking of space debris using the Murchison Widefield Array. I. Simulations and test observations demonstrate feasibility

Authors: S. J. Tingay, D. L. Kaplan, B. McKinley, F. Briggs, R. B. Wayth, N. Hurley-Walker, J. Kennewell, C. Smith, K. Zhang, W. Arcus, R. Bhat, D. Emrich, D. Herne, N. Kudryavtseva, M. Lynch, S. M. Ord, M. Waterson, D. G. Barnes, M. Bell, B. M. Gaensler, E. Lenc, G. Bernardi, L. J. Greenhill, J. C. Kasper, J. D. Bowman , et al. (38 additional authors not shown)

Abstract: The Murchison Widefield Array (MWA) is a new low frequency interferomeric radio telescope. The MWA is the low frequency precursor to the Square Kilometre Array (SKA) and is the first of three SKA precursors to be operational, supporting a varied science mission ranging from the attempted detection of the Epoch of Reionisation to the monitoring of solar flares and space weather. We explore the poss… ▽ More The Murchison Widefield Array (MWA) is a new low frequency interferomeric radio telescope. The MWA is the low frequency precursor to the Square Kilometre Array (SKA) and is the first of three SKA precursors to be operational, supporting a varied science mission ranging from the attempted detection of the Epoch of Reionisation to the monitoring of solar flares and space weather. We explore the possibility that the MWA can be used for the purposes of Space Situational Awareness (SSA). In particular we propose that the MWA can be used as an element of a passive radar facility operating in the frequency range 87.5 - 108 MHz (the commercial FM broadcast band). In this scenario the MWA can be considered the receiving element in a bi-static radar configuration, with FM broadcast stations serving as non-cooperative transmitters. The FM broadcasts propagate into space, are reflected off debris in Earth orbit, and are received at the MWA. The imaging capabilities of the MWA can be used to simultaneously detect multiple pieces of space debris, image their positions on the sky as a function of time, and provide tracking data that can be used to determine orbital parameters. Such a capability would be a valuable addition to Australian and global SSA assets, in terms of southern and eastern hemispheric coverage. We provide a feasibility assessment of this proposal, based on simple calculations and electromagnetic simulations that shows the detection of sub-metre size debris should be possible (debris radius of >0.5 m to ~1000 km altitude). We also present a proof-of-concept set of observations that demonstrate the feasibility of the proposal, based on the detection and tracking of the International Space Station via reflected FM broadcast signals originating in south-west Western Australia. These observations broadly validate our calculations and simulations. △ Less

Submitted 12 August, 2013; originally announced August 2013.

Comments: 24 pages, 5 figures, accepted by The Astronomical Journal. Abstract abridged here due to character number limits

Helioseismology of Pre-Emerging Active Regions III: Statistical Analysis

Authors: G. Barnes, A. C. Birch, K. D. Leka, D. C. Braun

Abstract: The subsurface properties of active regions prior to their appearance at the solar surface may shed light on the process of active region formation. Helioseismic holography has been applied to samples taken from two populations of regions on the Sun (pre-emergence and without emergence), each sample having over 100 members, that were selected to minimize systematic bias, as described in Paper I (L… ▽ More The subsurface properties of active regions prior to their appearance at the solar surface may shed light on the process of active region formation. Helioseismic holography has been applied to samples taken from two populations of regions on the Sun (pre-emergence and without emergence), each sample having over 100 members, that were selected to minimize systematic bias, as described in Paper I (Leka et al., 2012). Paper II (Birch et al., 2012) showed that there are statistically significant signatures in the average helioseismic properties that precede the formation of an active region. This paper describes a more detailed analysis of the samples of pre-emergence regions and regions without emergence, based on discriminant analysis. The property that is best able to distinguish the populations is found to be the surface magnetic field, even a day before the emergence time. However, after accounting for the correlations between the surface field and the quantities derived from helioseismology, there is still evidence of a helioseismic precursor to active region emergence that is present for at least a day prior to emergence. △ Less

Submitted 7 July, 2013; originally announced July 2013.

Comments: Submitted to ApJ

Journal ref: ApJ, 786, 19 (2014)

A 189 MHz, 2400 square degree polarization survey with the Murchison Widefield Array 32-element prototype

Authors: G. Bernardi, L. J. Greenhill, D. A. Mitchell, S. M. Ord, B. J. Hazelton, B. M. Gaensler, A. de Oliveira-Costa, M. F. Morales, N. Udaya Shankar, R. Subrahmanyan, R. B. Wayth, E. Lenc, C. L. Williams, W. Arcus, S. B. Arora, D. G. Barnes, J. D. Bowman, F. H. Briggs, J. D. Bunton, R. J. Cappallo, B. E. Corey, A. Deshpande, L. deSouza, D. Emrich, R. Goeke , et al. (28 additional authors not shown)

Abstract: We present a Stokes I, Q and U survey at 189 MHz with the Murchison Widefield Array 32-element prototype covering 2400 square degrees. The survey has a 15.6 arcmin angular resolution and achieves a noise level of 15 mJy/beam. We demonstrate a novel interferometric data analysis that involves calibration of drift scan data, integration through the co-addition of warped snapshot images and deconvolu… ▽ More We present a Stokes I, Q and U survey at 189 MHz with the Murchison Widefield Array 32-element prototype covering 2400 square degrees. The survey has a 15.6 arcmin angular resolution and achieves a noise level of 15 mJy/beam. We demonstrate a novel interferometric data analysis that involves calibration of drift scan data, integration through the co-addition of warped snapshot images and deconvolution of the point spread function through forward modeling. We present a point source catalogue down to a flux limit of 4 Jy. We detect polarization from only one of the sources, PMN J0351-2744, at a level of 1.8 \pm 0.4%, whereas the remaining sources have a polarization fraction below 2%. Compared to a reported average value of 7% at 1.4 GHz, the polarization fraction of compact sources significantly decreases at low frequencies. We find a wealth of diffuse polarized emission across a large area of the survey with a maximum peak of ~13 K, primarily with positive rotation measure values smaller than +10 rad/m^2. The small values observed indicate that the emission is likely to have a local origin (closer than a few hundred parsecs). There is a large sky area at 2^h30^m where the diffuse polarized emission rms is fainter than 1 K. Within this area of low Galactic polarization we characterize the foreground properties in a cold sky patch at $(α,δ) = (4^h,-27^\circ.6)$ in terms of three dimensional power spectra △ Less

Submitted 5 June, 2013; v1 submitted 26 May, 2013; originally announced May 2013.

Comments: 18 pages, 14 figures. Accepted for publication in ApJ. A version with high resolution figures is available at http://www.astro.rug.nl/~bernardi/32T_survey/32T_survey.pdf

Helioseismology of Pre-Emerging Active Regions I: Overview, Data, and Target Selection Criteria

Authors: K. D. Leka, G. Barnes, A. C. Birch, I. Gonzalez-Hernandez, T. Dunn, B. Javornik, D. C. Braun

Abstract: This first paper in a series describes the design of a study testing whether pre-appearance signatures of solar magnetic active regions were detectable using various tools of local helioseismology. The ultimate goal is to understand flux-emergence mechanisms by setting observational constraints on pre-appearance subsurface changes, for comparison with results from simulation efforts. This first pa… ▽ More This first paper in a series describes the design of a study testing whether pre-appearance signatures of solar magnetic active regions were detectable using various tools of local helioseismology. The ultimate goal is to understand flux-emergence mechanisms by setting observational constraints on pre-appearance subsurface changes, for comparison with results from simulation efforts. This first paper provides details of the data selection and preparation of the samples, each containing over 100 members, of two populations: regions on the Sun that produced a numbered NOAA active region, and a "control" sample of areas that did not. The seismology is performed on data from the GONG network; accompanying magnetic data from SOHO/MDI are used for co-temporal analysis of the surface magnetic field. Samples are drawn from 2001 -- 2007, and each target is analyzed for 27.7 hr prior to an objectively determined time of emergence. The results of two analysis approaches are published separately: one based on averages of the seismology- and magnetic-derived signals over the samples, another based on Discriminant Analysis of these signals, for a statistical test of detectable differences between the two populations. We include here descriptions of a new potential-field calculation approach and the algorithm for matching sample distributions over multiple variables. We describe known sources of bias and the approaches used to mitigate them. We also describe unexpected bias sources uncovered during the course of the study and include a discussion of refinements that should be included in future work on this topic. △ Less

Submitted 6 March, 2013; originally announced March 2013.

Comments: 24 pages, 11 figures, 2 tables. Published in ApJ

Journal ref: The Astrophysical Journal, Volume 762, Issue 2, article id. 130, 15 pp. (2013)

Helioseismology of Pre-Emerging Active Regions II: Average Emergence Properties

Authors: A. C. Birch, D. C. Braun, K. D. Leka, G. Barnes, B. Javornik

Abstract: We report on average subsurface properties of pre-emerging active regions as compared to areas where no active region emergence was detected. Helioseismic holography is applied to samples of the two populations (pre-emergence and without emergence), each sample having over 100 members, which were selected to minimize systematic bias, as described in Leka et al. We find that there are statistically… ▽ More We report on average subsurface properties of pre-emerging active regions as compared to areas where no active region emergence was detected. Helioseismic holography is applied to samples of the two populations (pre-emergence and without emergence), each sample having over 100 members, which were selected to minimize systematic bias, as described in Leka et al. We find that there are statistically significant signatures (i.e., difference in the means of more than a few standard errors) in the average subsurface flows and the apparent wave speed that precede the formation of an active region. The measurements here rule out spatially extended flows of more than about 15 m/s in the top 20 Mm below the photosphere over the course of the day preceding the start of visible emergence. These measurements place strong constraints on models of active region formation. △ Less

Submitted 6 March, 2013; originally announced March 2013.

Comments: 15 pages, 10 figures, ApJ (published)

Journal ref: The Astrophysical Journal, Volume 762, Issue 2, article id. 131, 12 pp. (2013)

Regular frequency patterns in the young delta Scuti star HD 261711 observed by the CoRoT and MOST satellites

Authors: K. Zwintz, L. Fossati, D. B. Guenther, T. Ryabchikova, A. Baglin, N. Themessl, T. G. Barnes, J. M. Matthews, M. Auvergne, D. Bohlender, S. Chaintreuil, R. Kuschnig, A. F. J. Moffat, J. F. Rowe, S. M. Rucinski, D. Sasselov, W. W. Weiss

Abstract: We concentrate on an asteroseismological study of HD 261711, a rather hot delta Scuti type pulsating member of the young open cluster NGC 2264 located at the blue border of the instability region. HD 261711 was discovered to be a pre-main sequence delta Scuti star using the time series photometry obtained by the MOST satellite in 2006. High-precision, time-series photometry of HD 261711 was obtain… ▽ More We concentrate on an asteroseismological study of HD 261711, a rather hot delta Scuti type pulsating member of the young open cluster NGC 2264 located at the blue border of the instability region. HD 261711 was discovered to be a pre-main sequence delta Scuti star using the time series photometry obtained by the MOST satellite in 2006. High-precision, time-series photometry of HD 261711 was obtained by the MOST and CoRoT satellites in 4 separate new observing runs that are put into context with the star's fundamental atmospheric parameters obtained from spectroscopy. With the new MOST data set from 2011/12 and the two CoRoT light curves from 2008 and 2011/12, the delta Scuti variability was confirmed and regular groups of frequencies were discovered. The two pulsation frequencies identified in the data from the first MOST observing run in 2006 are confirmed and 23 new delta Scuti-type frequencies were discovered using the CoRoT data. Weighted average frequencies for each group are related to l=0 and l=1 p-modes. Evidence for amplitude modulation of the frequencies in two groups is seen. The effective temperature was derived to be 8600$\pm$200 K, log g is 4.1$\pm$0.2, and the projected rotational velocity is 53$\pm$1km/s. Using our Teff value and the radius of 1.8$\pm$0.5 Rsun derived from SED fitting, we get a log L/Lsun of 1.20$\pm$0.14 which agrees well to the seismologically determined values of 1.65 Rsun and, hence, a log L/Lsun of 1.13. The radial velocity of 14$\pm$2 km/s we derived for HD 261711, confirms the star's membership to NGC 2264. Our asteroseismic models suggest that HD 261711 is a delta Scuti-type star close to the zero-age main sequence (ZAMS) with a mass of 1.8 to 1.9Msun. HD 261711 is either a young ZAMS star or a late PMS star just before the onset of hydrogen-core burning. △ Less

Submitted 14 February, 2013; originally announced February 2013.

Comments: 12 pages, 13 figures, A&A accepted

gamma Doradus pulsation in two pre-main sequence stars discovered by CoRoT

Authors: K. Zwintz, L. Fossati, T. Ryabchikova, A. Kaiser, M. Gruberbauer, T. G. Barnes, A. Baglin, S. Chaintreuil

Abstract: Pulsations in pre-main sequence stars have been discovered several times within the last years. But nearly all of these pulsators are of delta Scuti-type. gamma Doradus-type pulsation in young stars has been predicted by theory, but lack observational evidence. We present the investigation of variability caused by rotation and (gammaDoradus-type) pulsation in two pre-main sequence members of the y… ▽ More Pulsations in pre-main sequence stars have been discovered several times within the last years. But nearly all of these pulsators are of delta Scuti-type. gamma Doradus-type pulsation in young stars has been predicted by theory, but lack observational evidence. We present the investigation of variability caused by rotation and (gammaDoradus-type) pulsation in two pre-main sequence members of the young open cluster NGC2264 using high-precision time series photometry from the CoRoT satellite and dedicated high-resolution spectroscopy. Time series photometry of NGC2264VAS20 and NGC 2264VAS87 was obtained by the CoRoT satellite during the dedicated short run SRa01 in March 2008. NGC2264VAS87 was re-observed by CoRoT during the short run SRa05 in December 2011 and January 2012. Frequency analysis was conducted using Period04 and SigSpec. The spectral analysis was performed using equivalent widths and spectral synthesis. The frequency analysis yielded 10 and 14 intrinsic frequencies for NGC2264VAS20 and NGC2264VAS 87, respectively, in the range from 0 to 1.5c/d which are attributed to be caused by a combination of rotation and pulsation. The effective temperatures were derived to be 6380$\pm$150K for NGC2264VAS20 and 6220$\pm$150K for NGC2264VAS87. Membership of the two stars to the cluster is confirmed independently using X-ray fluxes, radial velocity measurements and proper motions available in the literature. The derived Li abundances of log n(Li)=3.34 and 3.54 for NGC2264VAS20 and NGC2264VAS87, respectively, are in agreement with the Li abundance for other stars in NGC2264 of similar Teff reported in the literature. We conclude that the two objects are members of NGC2264 and therefore are in their pre-main sequence evolutionary stage. Assuming that part of their variability is caused by pulsation, these two stars might be the first pre-main sequence gamma Doradus candidates. △ Less

Submitted 6 January, 2013; originally announced January 2013.

Comments: 11 pages, 10 figures, A&A accepted

Science with the Murchison Widefield Array

Authors: Judd D. Bowman, Iver Cairns, David L. Kaplan, Tara Murphy, Divya Oberoi, Lister Staveley-Smith, Wayne Arcus, David G. Barnes, Gianni Bernardi, Frank H. Briggs, Shea Brown, John D. Bunton, Adam J. Burgasser, Roger J. Cappallo, Shami Chatterjee, Brian E. Corey, Anthea Coster, Avinash Deshpande, Ludi deSouza, David Emrich, Philip Erickson, Robert F. Goeke, B. M. Gaensler, Lincoln J. Greenhill, Lisa Harvey-Smith , et al. (36 additional authors not shown)

Abstract: Significant new opportunities for astrophysics and cosmology have been identified at low radio frequencies. The Murchison Widefield Array is the first telescope in the Southern Hemisphere designed specifically to explore the low-frequency astronomical sky between 80 and 300 MHz with arcminute angular resolution and high survey efficiency. The telescope will enable new advances along four key scien… ▽ More Significant new opportunities for astrophysics and cosmology have been identified at low radio frequencies. The Murchison Widefield Array is the first telescope in the Southern Hemisphere designed specifically to explore the low-frequency astronomical sky between 80 and 300 MHz with arcminute angular resolution and high survey efficiency. The telescope will enable new advances along four key science themes, including searching for redshifted 21 cm emission from the epoch of reionisation in the early Universe; Galactic and extragalactic all-sky southern hemisphere surveys; time-domain astrophysics; and solar, heliospheric, and ionospheric science and space weather. The Murchison Widefield Array is located in Western Australia at the site of the planned Square Kilometre Array (SKA) low-band telescope and is the only low-frequency SKA precursor facility. In this paper, we review the performance properties of the Murchison Widefield Array and describe its primary scientific objectives. △ Less

Submitted 21 February, 2013; v1 submitted 20 December, 2012; originally announced December 2012.

Comments: Updated with revisions. 32 pages including figures and references. Submitted to PASA

Tera-scale Astronomical Data Analysis and Visualization

Authors: A. H. Hassan, C. J. Fluke, D. G. Barnes, V. A. Kilborn

Abstract: We present a high-performance, graphics processing unit (GPU)-based framework for the efficient analysis and visualization of (nearly) terabyte (TB)-sized 3-dimensional images. Using a cluster of 96 GPUs, we demonstrate for a 0.5 TB image: (1) volume rendering using an arbitrary transfer function at 7--10 frames per second; (2) computation of basic global image statistics such as the mean intensit… ▽ More We present a high-performance, graphics processing unit (GPU)-based framework for the efficient analysis and visualization of (nearly) terabyte (TB)-sized 3-dimensional images. Using a cluster of 96 GPUs, we demonstrate for a 0.5 TB image: (1) volume rendering using an arbitrary transfer function at 7--10 frames per second; (2) computation of basic global image statistics such as the mean intensity and standard deviation in 1.7 s; (3) evaluation of the image histogram in 4 s; and (4) evaluation of the global image median intensity in just 45 s. Our measured results correspond to a raw computational throughput approaching one teravoxel per second, and are 10--100 times faster than the best possible performance with traditional single-node, multi-core CPU implementations. A scalability analysis shows the framework will scale well to images sized 1 TB and beyond. Other parallel data analysis algorithms can be added to the framework with relative ease, and accordingly, we present our framework as a possible solution to the image analysis and visualization requirements of next-generation telescopes, including the forthcoming Square Kilometre Array pathfinder radiotelescopes. △ Less

Submitted 20 November, 2012; originally announced November 2012.

Comments: 16 pages, 14 Figures, accepted for publication in Monthly Notices of the Royal Astronomical Society

Cepheid distances from the Baade-Wesselink method

Authors: Wolfgang Gieren, Jesper Storm, Nicolas Nardetto, Alexandre Gallenne, Grzegorz Pietrzyński, Pascal Fouqué, Thomas G. Barnes, Daniel Majaess

Abstract: Recent progress on Baade-Wesselink (BW)-type techniques to determine the distances to classical Cepheids is reviewed. Particular emphasis is placed on the near-infrared surface-brightness (IRSB) version of the BW method. Its most recent calibration is described and shown to be capable of yielding individual Cepheid distances accurate to 6%, including systematic uncertainties. Cepheid distances fro… ▽ More Recent progress on Baade-Wesselink (BW)-type techniques to determine the distances to classical Cepheids is reviewed. Particular emphasis is placed on the near-infrared surface-brightness (IRSB) version of the BW method. Its most recent calibration is described and shown to be capable of yielding individual Cepheid distances accurate to 6%, including systematic uncertainties. Cepheid distances from the IRSB method are compared to those determined from open cluster zero-age main-sequence fitting for Cepheids located in Galactic open clusters, yielding excellent agreement between the IRSB and cluster Cepheid distance scales. Results for the Cepheid period-luminosity (PL) relation in near-infrared and optical bands based on IRSB distances and the question of the universality of the Cepheid PL relation are discussed. Results from other implementations of the BW method are compared to the IRSB distance scale and possible reasons for discrepancies are identified. △ Less

Submitted 26 October, 2012; originally announced October 2012.

Comments: Advancing the Physics of Cosmic Distances, Proceedings IAU Symposium No. 289, 2012

A Distributed GPU-based Framework for real-time 3D Volume Rendering of Large Astronomical Data Cubes

Authors: A. H. Hassan, C. J. Fluke, D. G. Barnes

Abstract: We present a framework to interactively volume-render three-dimensional data cubes using distributed ray-casting and volume bricking over a cluster of workstations powered by one or more graphics processing units (GPUs) and a multi-core CPU. The main design target for this framework is to provide an in-core visualization solution able to provide three-dimensional interactive views of terabyte-size… ▽ More We present a framework to interactively volume-render three-dimensional data cubes using distributed ray-casting and volume bricking over a cluster of workstations powered by one or more graphics processing units (GPUs) and a multi-core CPU. The main design target for this framework is to provide an in-core visualization solution able to provide three-dimensional interactive views of terabyte-sized data cubes. We tested the presented framework using a computing cluster comprising 64 nodes with a total of 128 GPUs. The framework proved to be scalable to render a 204 GB data cube with an average of 30 frames per second. Our performance analyses also compare between using NVIDIA Tesla 1060 and 2050 GPU architectures and the effect of increasing the visualization output resolution on the rendering performance. Although our initial focus, and the examples presented in this work, is volume rendering of spectral data cubes from radio astronomy, we contend that our approach has applicability to other disciplines where close to real-time volume rendering of terabyte-order 3D data sets is a requirement. △ Less

Submitted 1 May, 2012; originally announced May 2012.

Comments: 13 Pages, 7 figures, has been accepted for publication in Publications of the Astronomical Society of Australia