No Arabic abstract
The expert system for time series analysis of irregularly spaced signals is reviewed. It consists of a number of complementary algorithms and programs, which may be effective for different types of variability. Obviously, for a pure sine signal, all the methods should produce the same results. However, for irregularly spaced signals with a complicated structure, e.g. a sum of different components, different methods may produce significantly different results. The basic approach is based on classical method of the least squares (1994OAP.....7...49A). However, contrary to common step-by-step methods of removal important components (e.g. mean, trend (detrending), sine wave (prewhitening), where covariations between different components are ignored, i.e. erroneously assumed to be zero, we use complete mathematical models. Some of the methods are illustrated on the observations of the semi-regular pulsating variable RY UMa. The star shows a drastic cyclic change of semi-amplitude of pulsations between 0.01 to 0.37mag, which is interpreted as a bias between the waves with close periods and a beat period of 4000d (11yr). The dominating period has changed from 307.35(8)d before 1993 to 285.26(6)d after 1993. The initial epoch of the maximum brightness for the recent interval is 2454008.8(5). It is suggested that the apparent period switch is due to variability of amplitudes of these two waves and an occasional swap of the dominating wave.
The exact period determination of a multi-periodic variable star based on its luminosity time series data is believed a task requiring skill and experience. Thus the majority of available time series analysis techniques require human intervention to some extent. The present work is dedicated to establish an automated method of period (or frequency) determination from the time series database of variable stars. Relying on the SigSpec method (Reegen 2007), the technique established here employs a statistically unbiased treatment of frequency-domain noise and avoids spurious (i. e. noise induced) and alias peaks to the highest possible extent. Several add-ons were incorporated to tailor SigSpec to our requirements. We present tests on 386 stars taken from ASAS2 project database. From the output file produced by SigSpec, the frequency with maximum spectral significance is chosen as the genuine frequency. Out of 386 variable stars available in the ASAS2 database, our results contain 243 periods recovered exactly and also 88 half periods, 42 different periods etc. SigSpec has the potential to be effectively used for fully automated period detection from variable stars time series database. The exact detection of periods helps us to identify the type of variability and classify the variable stars, which provides a crucial information on the physical processes effective in stellar atmospheres.
The software MAVKA is described, which was elaborated for statistically optimal determination of the characteristics of the extrema of 1000+ variable stars of different types, mainly eclipsing and pulsating. The approximations are phenomenological, but not physical. As often, the discovery of a new variable star is made on time series of a single-filter (single-channel) data, and there is no possibility to determine parameters needed for physical modelling (e.g. temperature, radial velocities, mass ratio of binaries). Besides classical polynomial approximation AP (we limited the degree of the polynomial from 2 to 9), there are realized symmetrical approximations (symmetrical polynomials SP, wall-supported horizontal line WSL and parabola WSP, restricted polynomials of non-integer order based on approximations of the functions proposed by Andronov (2012) and Mikulasek (2015) and generally asymmetric functions (asymptotic parabola AP, parabolic spline PS, generalized hyperbolic secant function SECH and log-normal-like BSK). This software is a successor of the Observation Obscurer with some features for the variable star research, including a block for running parabola RP scalegram and approximation. Whereas the RP is oriented on approximation of the complete data set. MAVKA is pointed to parts of the light curve close to extrema (including total eclipses and transits of stars and exoplanets). The functions for wider intervals, covering the eclipse totally, were discussed in 2017Ap.....60...57A . Global and local approximations are reviewed in 2020kdbd.book..191A . The software is available at http://uavso.org.ua/mavka and https://katerynaandrych.wixsite.com/mavka. We have analyzed the data from own observations, as well as from monitoring obtained at ground-based and space (currently, mainly, TESS) observatories. It may be used for signals of any nature.
With the advent of the Large Synoptic Survey Telescope (LSST), time-domain astronomy will be faced with an unprecedented volume and rate of data. Real-time processing of variables and transients detected by such large-scale surveys is critical to identifying the more unusual events and allocating scarce follow-up resources efficiently. We develop an algorithm to identify these novel events within a given population of variable sources. We determine the distributions of magnitude changes (dm) over time intervals (dt) for a given passband f, pf(dm|dt), and use these distributions to compute the likelihood of a test source being consistent with the population, or an outlier. We demonstrate our algorithm by applying it to the DECam multi-band time-series data of more than 2000 variable stars identified by Saha et al. (2019) in the Galactic Bulge that are largely dominated by long-period variables and pulsating stars. Our algorithm discovers 18 outlier sources in the sample, including a microlensing event, a dwarf nova, and two chromospherically active RS CVn stars, as well as sources in the Blue Horizontal Branch region of the color-magnitude diagram without any known counterparts. We compare the performance of our algorithm for novelty detection with multivariate KDE and Isolation Forest on the simulated PLAsTiCC dataset. We find that our algorithm yields comparable results despite its simplicity. Our method provides an efficient way for flagging the most unusual events in a real-time alert-broker system.
As extensive ground-based observations and characterisation of different variable stars are of the utmost importance in preparing optimal input catalogues for space missions, our aim was to search for new variable stars in selected fields of the northern sky. We obtained 24470 CCD images and analysed photometric time series of stars using the DAOPHOT based package Muniwin as the first step, and the Period04 package was used to further analyse the suspected new variable stars. The light curves and other observational results are presented for 3598 stars online. We found 81 new variable stars, among them is an eclipsing binary with a variable component and possibly eccentric orbits TYC4038-693-1 which we also observed spectroscopically, four {delta} Scuti candidates, six other variable stars with periods falling into the interval of 35 minutes to 20 days. Furthermore, we identified 70 slowly varying stars with so far undefined periodicity. Additional photometric and spectral observations were carried out for TYC 2764-1997-1, and its previous candidacy for eclipsing binaries was approved.
During the last decades there is a continuing international endeavor in developing realistic space weather prediction tools aiming to forecast the conditions on the Sun and in the interplanetary environment. These efforts have led to the need of developing appropriate metrics in order to assess the performance of those tools. Metrics are necessary for validating models, comparing different models and monitoring adjustments or improvements of a certain model over time. In this work, we introduce the Dynamic Time Warping (DTW) as an alternative way to validate models and, in particular, to quantify differences between observed and synthetic (modeled) time series for space weather purposes. We present the advantages and drawbacks of this method as well as applications on WIND observations and EUHFORIA modeled output at L1. We show that DTW is a useful tool that permits the evaluation of both the fast and slow solar wind. Its distinctive characteristic is that it warps sequences in time, aiming to align them with the minimum cost by using dynamic programming. It can be applied in two different ways for the evaluation of modeled solar wind time series. The first way calculates the so-called sequence similarity factor (SSF), a number that provides a quantification of how good the forecast is, compared to a best and a worst case prediction scenarios. The second way quantifies the time and amplitude differences between the points that are best matched between the two sequences. As a result, it can serve as a hybrid metric between continuous measurements (such as, e.g., the correlation coefficient) and point-by-point comparisons. We conclude that DTW is a promising technique for the assessment of solar wind profiles offering functions that other metrics do not, so that it can give at once the most complete evaluation profile of a model.