The Astro-WISE optical image pipeline

We have designed and implemented a novel way to process wide-field astronomical data within a distributed environment of hardware resources and humanpower. The system is characterized by integration of archiving, calibration, and post-calibration analysis of data from raw, through intermediate, to final data products. It is a true integration thanks to complete linking of data lineage from the final catalogs back to the raw data. This paper describes the pipeline processing of optical wide-field astronomical data from the WFI (http://www.eso.org/lasilla/instruments/wfi/) and OmegaCAM (http://www.astro-wise.org/~omegacam/) instruments using the Astro-WISE information system (the Astro-WISE Environment or simply AWE). This information system is an environment of hardware resources and humanpower distributed over Europe. AWE is characterized by integration of archiving, data calibration, post-calibration analysis, and archiving of raw, intermediate, and final data products. The true integration enables a complete data processing cycle from the raw data up to the publication of science-ready catalogs. The advantages of this system for very large datasets are in the areas of: survey operations management, quality control, calibration analyses, and massive processing.     

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

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

EIDOSCOPE: particle acceleration at plasma boundaries

We describe the mission concept of how ESA can make a major contribution to the Japanese Canadian multi-spacecraft mission SCOPE by adding one cost-effective spacecraft EIDO (Electron and Ion Dynamics Observatory), which has a comprehensive and optimized plasma payload to address the physics of particle acceleration. The combined mission EIDOSCOPE will distinguish amongst and quantify the governing processes of particle acceleration at several important plasma boundaries and their associated boundary layers: collisionless shocks, plasma jet fronts, thin current sheets and turbulent boundary layers. Particle acceleration and associated cross-scale coupling is one of the key outstanding topics to be addressed in the Plasma Universe. The very important science questions that only the combined EIDOSCOPE mission will be able to tackle are: 1) Quantitatively, what are the processes and efficiencies with which both electrons and ions are selectively injected and subsequently accelerated by collisionless shocks? 2) How does small-scale electron and ion acceleration at jet fronts due to kinetic processes couple simultaneously to large scale acceleration due to fluid (MHD) mechanisms? 3) How does multi-scale coupling govern acceleration mechanisms at electron, ion and fluid scales in thin current sheets? 4) How do particle acceleration processes inside turbulent boundary layers depend on turbulence properties at ion/electron scales? EIDO particle instruments are capable of resolving full 3D particle distribution functions in both thermal and suprathermal regimes and at high enough temporal resolution to resolve the relevant scales even in very dynamic plasma processes. The EIDO spin axis is designed to be sun-pointing, allowing EIDO to carry out the most sensitive electric field measurements ever accomplished in the outer magnetosphere. Combined with a nearby SCOPE Far Daughter satellite, EIDO will form a second pair (in addition to SCOPE Mother-Near Daughter) of closely separated satellites that provides the unique capability to measure the 3D electric field with high accuracy and sensitivity. All EIDO instrumentation are state-of-the-art technology with heritage from many recent missions. The EIDOSCOPE orbit will be close to equatorial with apogee 25-30 RE and perigee 8-10 RE. In the course of one year the orbit will cross all the major plasma boundaries in the outer magnetosphere; bow shock, magnetopause and magnetotail current sheets, jet fronts and turbulent boundary layers. EIDO offers excellent cost/benefits for ESA, as for only a fraction of an M-class mission cost ESA can become an integral part of a major multi-agency L-class level mission that addresses outstanding science questions for the benefit of the European science community.     

Solar Flare Prediction Using Advanced Feature Extraction, Machine Learning, and Feature Selection

Novel machine-learning and feature-selection algorithms have been developed to study: i) the flare-prediction-capability of magnetic feature (MF) properties generated by the recently developed Solar Monitor Active Region Tracker (SMART); ii) SMART’s MF properties that are most significantly related to flare occurrence. Spatiotemporal association algorithms are developed to associate MFs with flares from April 1996 to December 2010 in order to differentiate flaring and non-flaring MFs and enable the application of machine-learning and feature-selection algorithms. A machine-learning algorithm is applied to the associated datasets to determine the flare-prediction-capability of all 21 SMART MF properties. The prediction performance is assessed using standard forecast-verification measures and compared with the prediction measures of one of the standard technologies for flare-prediction that is also based on machine-learning: Automated Solar Activity Prediction (ASAP). The comparison shows that the combination of SMART MFs with machine-learning has the potential to achieve more accurate flare-prediction than ASAP. Feature-selection algorithms are then applied to determine the MF properties that are most related to flare occurrence. It is found that a reduced set of six MF properties can achieve a similar degree of prediction accuracy as the full set of 21 SMART MF properties.     

TomograPy: A Fast, Instrument-Independent, Solar Tomography Software

Solar tomography has progressed rapidly in recent years thanks to the development of robust algorithms and the availability of more powerful computers. It can today provide crucial insights in solving issues related to the line-of-sight integration present in the data of solar imagers and coronagraphs. However, there remain challenges such as the increase of the available volume of data, the handling of the temporal evolution of the observed structures, and the heterogeneity of the data in multi-spacecraft studies. We present a generic software package that can perform fast tomographic inversions that scales linearly with the number of measurements, linearly with the length of the reconstruction cube (and not the number of voxels), and linearly with the number of cores and can use data from different sources and with a variety of physical models: TomograPy ( http://nbarbey.github.com/TomograPy/ ), an open-source software freely available on the Python Package Index. For performance, TomograPy uses a parallelized-projection algorithm. It relies on the World Coordinate System standard to manage various data sources. A variety of inversion algorithms are provided to perform the tomographic-map estimation. A test suite is provided along with the code to ensure software quality. Since it makes use of the Siddon algorithm it is restricted to rectangular parallelepiped voxels but the spherical geometry of the corona can be handled through proper use of priors. We describe the main features of the code and show three practical examples of multi-spacecraft tomographic inversions using STEREO/EUVI and STEREO/COR1 data. Static and smoothly varying temporal evolution models are presented.     

Evaluation of a Revised Interplanetary Shock Prediction Model: 1D CESE-HD-2 Solar-Wind Model

We modified the one-dimensional conservation element and solution element (CESE) hydrodynamic (HD) model into a new version [1D CESE-HD-2], by considering the direction of the shock propagation. The real-time performance of the 1D CESE-HD-2 model during Solar Cycle 23 (February 1997?–?December 2006) is investigated and compared with those of the Shock Time of Arrival Model (STOA), the Interplanetary-Shock-Propagation Model (ISPM), and the Hakamada–Akasofu–Fry version 2 (HAFv.2). Of the total of 584 flare events, 173 occurred during the rising phase, 166 events during the maximum phase, and 245 events during the declining phase. The statistical results show that the success rates of the predictions by the 1D CESE-HD-2 model for the rising, maximum, declining, and composite periods are 64 %, 62 %, 57 %, and 61 %, respectively, with a hit window of ±?24 hours. The results demonstrate that the 1D CESE-HD-2 model shows the highest success rates when the background solar-wind speed is relatively fast. Thus, when the background solar-wind speed at the time of shock initiation is enhanced, the forecasts will provide potential values to the customers. A high value (27.08) of χ ² and low p-value (<?0.0001) for the 1D CESE-HD-2 model give considerable confidence for real-time forecasts by using this new model. Furthermore, the effects of various shock characteristics (initial speed, shock duration, background solar wind, longitude, etc.) and background solar wind on the forecast are also investigated statistically.     

Modified Homogeneous Data Set of Coronal Intensities

The Astronomical Institute of the Slovak Academy of Sciences has published the intensities, recalibrated with respect to a common intensity scale, of the 530.3 nm (Fe xiv) green coronal line observed at ground-based stations up to the year 2008. The name of this publication is Homogeneous Data Set (HDS). We have developed a method that allows one to successfully substitute the ground-based observations by satellite observations and, thus, continue with the publication of the HDS. For this purpose, the observations of the Extreme-ultraviolet Imaging Telescope (EIT), onboard the Solar and Heliospheric Observatory (SOHO) satellite, were exploited. Among other data the EIT instrument provides almost daily 28.4 nm (Fe xv) emission-line snapshots of the corona. The Fe xiv and Fe xv data (4051 observation days) taken in the period 1996?–?2008 have been compared and good agreement was found. The method to obtain the individual data for the HDS follows from the correlation analysis described in this article. The resulting data, now under the name of Modified Homogeneous Data Set (MHDS), are identical up to 1996 to those in the HDS. The MHDS can be used further for studies of the coronal solar activity and its cycle. These data are available at http://www.suh.sk .     

Modified Coronal Index of the Solar Activity

The original coronal index of the solar activity (CI) has been constructed on the basis of ground-based measurements of the intensities of the coronal line of 530.3 nm (Rybanský in Bull. Astron. Inst. Czechoslov., 28, 367, 1975; Rybanský et al. in J. Geophys. Res., 110, A08106, 2005). In this paper, CI is compared with the EUV measurements on the CELIAS/SEM equipment based on the same idea as the original idea of the coronal index. The correlation is very good for the period 1996?–?2005 (r=0.94 for daily values). The principal result of this paper is the introduction of the modified coronal index (MCI) which in all uses and contexts can replace the existing CI index. Daily MCI values extend over a time period of six solar activity cycles. Future MCI measurements will be derived from more reliable measurements made by space-based observatories that are not influenced by the weather. MCI measurements are and will continue to be archived at the web site of the Slovak Central Observatory in Hurbanovo ( http://www.suh.sk/obs/vysl/MCI.htm ).     

The unification of the inflation with late-time acceleration in Born-Infeld-f(R) gravity

We study accelerating dynamics from Born-Infeld-f(R) gravity in a simplified conformal approach without matter. In Makarenko et al. (arXiv:1404.2850 [gr-qc], 2011b) it was derived eventually any Dark Energy cosmology from above theory. In this Letter we apply the technique of Makarenko et al. (arXiv:1404.2850 [gr-qc], 2011b) to show that Born-Infeld-f(R) gravity may describe very realistic universe admitting the unification of early-time inflation with late-time acceleration. Specifically, the evolution with periodic as well as non-periodic behavior is considered with possibility to cross the phantom-divide at early or late-times.     

A Multi-wavelength Analysis of Active Regions and Sunspots by Comparison of Automatic Detection Algorithms

Since the Solar Dynamics Observatory (SDO) began recording ≈?1 TB of data per day, there has been an increased need to automatically extract features and events for further analysis. Here we compare the overall detection performance, correlations between extracted properties, and usability for feature tracking of four solar feature-detection algorithms: the Solar Monitor Active Region Tracker (SMART) detects active regions in line-of-sight magnetograms; the Automated Solar Activity Prediction code (ASAP) detects sunspots and pores in white-light continuum images; the Sunspot Tracking And Recognition Algorithm (STARA) detects sunspots in white-light continuum images; the Spatial Possibilistic Clustering Algorithm (SPoCA) automatically segments solar EUV images into active regions (AR), coronal holes (CH), and quiet Sun (QS). One month of data from the Solar and Heliospheric Observatory (SOHO)/Michelson Doppler Imager (MDI) and SOHO/Extreme Ultraviolet Imaging Telescope (EIT) instruments during 12 May?–?23 June 2003 is analysed. The overall detection performance of each algorithm is benchmarked against National Oceanic and Atmospheric Administration (NOAA) and Solar Influences Data Analysis Center (SIDC) catalogues using various feature properties such as total sunspot area, which shows good agreement, and the number of features detected, which shows poor agreement. Principal Component Analysis indicates a clear distinction between photospheric properties, which are highly correlated to the first component and account for 52.86% of variability in the data set, and coronal properties, which are moderately correlated to both the first and second principal components. Finally, case studies of NOAA 10377 and 10365 are conducted to determine algorithm stability for tracking the evolution of individual features. We find that magnetic flux and total sunspot area are the best indicators of active-region emergence. Additionally, for NOAA 10365, it is shown that the onset of flaring occurs during both periods of magnetic-flux emergence and complexity development.     

Interpreting the Helioseismic and Magnetic Imager (HMI) Multi-Height Velocity Measurements

The Solar Dynamics Observatory/Helioseismic and Magnetic Imager (SDO/HMI) filtergrams, taken at six wavelengths around the Fe i 6173.3 Å line, contain information about the line-of-sight velocity over a range of heights in the solar atmosphere. Multi-height velocity inferences from these observations can be exploited to study wave motions and energy transport in the atmosphere. Using realistic convection-simulation datasets provided by the STAGGER and MURaM codes, we generate synthetic filtergrams and explore several methods for estimating Dopplergrams. We investigate at which height each synthetic Dopplergram correlates most strongly with the vertical velocity in the model atmospheres. On the basis of the investigation, we propose two Dopplergrams other than the standard HMI-algorithm Dopplergram produced from HMI filtergrams: a line-center Dopplergram and an average-wing Dopplergram. These two Dopplergrams correlate most strongly with vertical velocities at the heights of 30?–?40 km above (line center) and 30?–?40 km below (average wing) the effective height of the HMI-algorithm Dopplergram. Therefore, we can obtain velocity information from two layers separated by about a half of a scale height in the atmosphere, at best. The phase shifts between these multi-height Dopplergrams from observational data as well as those from the simulated data are also consistent with the height-difference estimates in the frequency range above the photospheric acoustic-cutoff frequency.     

Sources and Databases for the Study of the Dynamics of Moons of Asteroids

A new database of all astrometric observations of moons of asteroids is offered. It has a simple structure and is accessible through the Internet. Regular database updating is provided when new observational results occur in publications. The database is located at the sites of the Natural Satellites Data Center created as a result of collaboration between the Sternberg Astronomical Institute of Moscow State University and the Institut de Mécanique Céleste et de Calcul des Ephémérides (IMCCE). The database addresses in the Internet are https://doi.org/www.sai.msu.ru/neb/nss/indexr.htm and https://doi.org/nsdb.imcce.fr/obsposrespectively. On entering these sites, we need to select one of three languages, the Observation item, and the Astrometric positions of asteroids with moons item. The observational data are provided here with explanations and hyperlinks to the publications in the bibliographic database SAO/NASA Astrophysics Data System (ADS) Abstract Service.     

EXOTIME: searching for planets around pulsating subdwarf B stars

In 2007, a companion with planetary mass was found around the pulsating subdwarf B star V391 Pegasi with the timing method, indicating that a previously undiscovered population of substellar companions to apparently single subdwarf B stars might exist. Following this serendipitous discovery, the EXOTIME (http://www.na.astro.it/~silvotti/exotime/) monitoring program has been set up to follow the pulsations of a number of selected rapidly pulsating subdwarf B stars on time scales of several years with two immediate observational goals:

1. determine $\dot{P}$ of the pulsational periods P
2. search for signatures of substellar companions in O–C residuals due to periodic light travel time variations, which would be tracking the central star’s companion-induced wobble around the centre of mass

These sets of data should therefore, at the same time, on the one hand be useful to provide extra constraints for classical asteroseismological exercises from the $\dot{P}$ (comparison with “local” evolutionary models), and on the other hand allow one to investigate the preceding evolution of a target in terms of possible “binary” evolution by extending the otherwise unsuccessful search for companions to potentially very low masses. While timing pulsations may be an observationally expensive method to search for companions, it samples a different range of orbital parameters, inaccessible through orbital photometric effects or the radial velocity method: the latter favours massive close-in companions, whereas the timing method becomes increasingly more sensitive toward wider separations. In this paper we report on the status of the on-going observations and coherence analysis for two of the currently five targets, revealing very well-behaved pulsational characteristics in HS?0444+0458, while showing HS?0702+6043 to be more complex than previously thought.     

IAU MDC Photographic Meteor Orbits Database: Version 2013

A new 2013 version of the IAU MDC photographic meteor orbits database which is an upgrade of the current 2003 version (Lindblad et al. 2003, EMP 93:249–260) is presented. To the 2003 version additional 292 orbits are added, thus the new version of the database consists of 4,873 meteors with their geophysical and orbital parameters compiled in 41 catalogues. For storing the data, a new format enabling a more simple treatment with the parameters, including the errors of their determination is applied. The data can be downloaded from the IAU MDC web site: http://www.astro.sk/IAUMDC/Ph2013/     

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

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

Validation of the 3D AMR SIP–CESE Solar Wind Model for Four Carrington Rotations

We carry out the adaptive mesh refinement (AMR) implementation of our solar–interplanetary space-time conservation element and solution element (CESE) magnetohydrodynamic model (SIP–CESE MHD model) using a six-component grid system (Feng, Zhou, and Wu, Astrophys. J. 655, 1110, 2007; Feng et al., Astrophys. J. 723, 300, 2010). By transforming the governing MHD equations from the physical space (x,y,z) to the computational space (ξ,η,ζ) while retaining the form of conservation (Jiang et al., Solar Phys. 267, 463, 2010), the SIP–AMR–CESE MHD model is implemented in the reference coordinates with the aid of the parallel AMR package PARAMESH available at http://sourceforge.net/projects/paramesh/ . Meanwhile, the volumetric heating source terms derived from the topology of the magnetic-field expansion factor and the minimum angular separation (at the photosphere) between an open-field foot point and its nearest coronal-hole boundary are also included. We show the preliminary results of applying the SIP–AMR–CESE MHD model for simulating the solar-wind background of different solar-activity phases by comparison with SOHO observations and other spacecraft data from OMNI. Our numerical results show overall good agreements in the solar corona and in interplanetary space with these multiple-spacecraft observations.     

The SP_PREP Data Preparation Package for the Hinode Spectro-Polarimeter

The Hinode/Spectro-Polarimeter (SP) is the first space-borne precision spectro-polarimeter for the study of solar phenomena. It is primarily intended for measuring the solar photospheric vector magnetic field at high spatial and spectral resolution. This objective requires that the data are calibrated and conditioned to a high degree of precision. We describe how the calibration package SP_PREP for the SP operates.