Search NASA⌕ Search

SEARCH · Search NASA

Results for “Lomb-Scargle periodogram”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

Long-Term X-Ray Variability of Circinus X-1

We present an analysis of long term X-ray monitoring observations of Circinus X-1 (Cir X-1) made with four different instruments: Vela 5B, Ariel V ASM, Ginga ASM, and RXTE ASM, over the course of more than 30 years. We use Lomb-Scargle periodograms to search for the approx. 16.5 day orbital period of Cir X-1 in each of these data sets and from this derive a new orbital ephemeris based solely on X-ray measurements, which we compare to the previous ephemerides obtained from radio observations. We also use the Phase Dispersion Minimization (PDM) technique, as well as FFT analysis, to verify the periods obtained from periodograms. We obtain dynamic periodograms (both Lomb-Scargle and PDM) of Cir X-1 during the RXTE era, showing the period evolution of Cir X-1, and also displaying some unexplained discrete jumps in the location of the peak power.

CIRCINUS X-1↗

DMTN-221: Periodicity Analysis in Alert Production

The baselined timeseries features to be computed in Alert Production include a Lomb-Scargle periodogram ran on two classes of variable systems: RR Lyrae and Eclipsing Binaries. Based on a simulated LSST-like cadence light curves taken from the Extended LSST Astronomical Time-series Classification Challenge (ELAsTiCC) we perform an end-to-end test to characterize the periodicity recovery on the Alert Production multi-band light curves. In both variable classes, we found that a single-band Lomb-Scargle implementation yields to a low fraction of recovered periods, with a significant preference on the simple periodic phenomena such as RR Lyrae. We also investigated the results from a multi-band Lomb-Scargle and found an increased fraction of recovered periodicities above 15% for the eclipsing binaries, and over 80$\%$ for the RR Lyrae stars. Our findings suggest that a multi-band Lomb-Scargle should be implemented for searching periodic phenomena through AP. We also asses the computational and scientific performance of several configurations on simulated alert data and find that our current configuration scales linearly with the number of detections while assuming an heuristic frequency grid.

79 ASTRONOMY AND ASTROPHYSICS↗

Coincident 1.3-year Periodicities in the ap Geomagnetic Index and the Solar Wind

Recent observations show an approximately 1.3-year period in the speed of the solar wind detected by the IMP 8 and Voyager 2 spacecraft. A similar period is also seen in the north-south (GSE) component of the magnetic field observed by IMP 8. Since both parameters are commonly used as input to models of geomagnetic activity, the 'ap' index (a measure of geomagnetic disturbance) is examined to look for this periodicity. The Lomb-Scargle periodogram method is used on the ap, plasma, and magnetic field data during the 1973-1994 time range. A dynamic FFT periodogram method is also used to analyze the ap data during this time, as well as to look for periods present between 1932 and 1972. A clear 1.3-year periodicity is present in the post-1986 data when the same period is observed in the plasma and field data. The V(2)B(zsm) and V(2)B(s) proxies for geomagnetic activity also show this periodicity. However, the southward (GSM) component of the magnetic field does not have a 1.3-year period, and neither do solar wind or ap data from 1973-1985. This demonstrates that the ap geomagnetic index can act as a proxy for solar wind periodicities at this time scale. Historic ap data are examined, and show that a similar periodicity in ap exists around 1942. Since auroral data show a 1.4-year periodicity, all these similar periods may result from a common underlying solar mechanism.

Paularena, K. I.↗

Seasonal Variability of Saturn's Atmosphere

The seasonal variability of Saturn's clouds and weather layer, currently displaying a variety of phenomena (convective storms, planetary waves, giant storms and lightning-induced events, etc.) is not yet fully understood. Variations of Saturn's radiance at 5.2 microns, a spectral region dominated by thermal emission in an atmospheric window containing weak gaseous absorption, contain a strong axisymmetric component as well as large discrete features at low and mid-latitudes that are several degrees colder than the planetary average and uncorrelated with features at shorter wavelengths that are dominated by reflected sunlight (Yanamandra-Fisher et al., 2001. Icarus, Vol. 150). The characterization of several fundamental atmospheric properties and processes, however, remains incomplete, namely: How do seasons affect (a) the global distribution of gaseous constituents and aerosols; and (b) temperatures and the stability against convection and large scale-atmospheric transport? Do 5-micron clouds have counterparts at other altitude levels? What changes occur during the emergence of Great White Storms? Data acquired at the NASA/IRTF and NAOJ/Subaru from 1995 - 2011; since 2004, high-resolution multi-spectral and high-spatial imaging data acquired by the NASA/ESA Cassini mission, represents half a Saturnian year or two seasons. With the addition of detailed multi-spectral data sets acquired by amateur observers, we study these dramatic phenomena to better understand the timeline of the evolution of these events. Seasonal (or temporal) trends in the observables such as albedo of the clouds, thermal fields of the atmosphere as function of altitude, development of clouds, hazes and global abundances of various hydrocarbons in the atmosphere can now be modeled. We will present results of our ongoing investigation for the search and characterization of periodicities over half a Saturnian year, based on a non-biased a priori approach and time series techniques (such as Principal Component Analysis, PCA and Lomb-Scargle periodograms, LSP).

Yanamandra-Fisher, Padma A.↗

A procedure of multiple period searching in unequally spaced time-series with the Lomb-Scargle method

Periodogram analysis of unequally spaced time-series, as part of many biological rhythm investigations, is complicated. The mathematical framework is scattered over the literature, and the interpretation of results is often debatable. In this paper, we show that the Lomb-Scargle method is the appropriate tool for periodogram analysis of unequally spaced data. A unique procedure of multiple period searching is derived, facilitating the assessment of the various rhythms that may be present in a time-series. All relevant mathematical and statistical aspects are considered in detail, and much attention is given to the correct interpretation of results. The use of the procedure is illustrated by examples, and problems that may be encountered are discussed. It is argued that, when following the procedure of multiple period searching, we can even benefit from the unequal spacing of a time-series in biological rhythm research.

NASA Discipline Space Human Factors↗

Red Noise–based False Alarm Thresholds for Astrophysical Periodograms via Whittle’s Approximation to the Likelihood

Astronomers who search for periodic signals using Lomb–Scargle periodograms rely on false alarm level (FAL) estimates to identify statistically significant peaks. Although FALs are often calculated from white noise models, many astronomical time series suffer from red noise. Prewhitening is a statistical technique in which a continuum model is subtracted from the log power spectrum estimate, after which the observer can proceed with a white-noise treatment. Here we present a prewhitening-based method of calculating frequency-dependent FALs. We fit power laws and autoregressive models of order 1 to each Lomb–Scargle periodogram by minimizing the Whittle approximation to the negative log-likelihood (NLL), then calculate FALs based on the best-fit model power spectrum. Our technique is a novel extension of the Whittle NLL to datasets with uneven time sampling. We demonstrate FAL calculations using observations of α Cen B, GJ 581, HD 192310, synthetic data from the radial velocity (RV) fitting challenge, and Kepler observations of a differential rotator. The Kepler data analysis shows that only true rotation signals are detected by red noise FALs, while white noise FALs suggest all spurious peaks in the low-frequency range are significant. A high-frequency sinusoid injected into α Cen B logR$'$ HK observations exceeds the 1% red noise FAL despite having only 8.9% of the power of the dominant rotation signal. In a periodogram of HD 192310 RVs, peaks associated with differential rotation and planets are detected against the 5% red noise FAL without iterative model fitting or subtraction. The software for calculating red noise–based FALs is available on GitHub.

Astrostatistics (1882)↗

Multitaper Magnitude‐Squared Coherence for Time Series With Missing Data: Understanding Oscillatory Processes Traced by Multiple Observables

To explore the hypothesis of a common source of variability in two time series, observers may estimate the magnitude-squared coherence (MSC), which is a frequency-domain view of the cross correlation. For time series that do not have uniform observing cadence, MSC can be estimated using Welch's overlapping segment averaging. However, multitaper has superior statistical properties to Welch's method in terms of the tradeoff between bias, variance, and bandwidth. The classical multitaper technique has recently been extended to accommodate time series with underlying uniform observing cadence from which some observations are missing. This situation is common for solar and geomagnetic data sets, which may have gaps due to breaks in satellite coverage, instrument downtime, or poor observing conditions. We demonstrate the scientific use of missing-data multitaper magnitude-squared coherence by detecting known solar mid-term oscillations in simultaneous, missing-data time series of solar Lyman α flux and geomagnetic Disturbance Storm Time index. Due to their superior statistical properties, we recommend that multitaper methods be used for all heliospheric time series with underlying uniform observing cadence.

Astro-statistics techniques (1886)↗

Missing-Data Nonparametric Coherency Estimation

Chave recently proposed an estimator for multitaper spectral density where the time series contains missing values. In this article we generalize this technique to a multitaper estimator of coherence and phase and show that one can also obtain bootstrapped confidence intervals. Additionally, we give two examples. The first is a toy example in which the true coherence is known. In the second example we show that the multitaper missing-data coherence estimator computed on real data with a single gap comprising 11% of the data outperforms the Daniell-smoothed coherence estimator where there are no gaps. The case where the two time series have different missing indices is also discussed.

42 ENGINEERING↗