Time series data are significant and are derived from temporal data, which involve real numbers representing values collected regularly over time. Time series have a great impact on many types of data. However, time series have anomalies. We introduce hybrid algorithm named novel matrix profile (NMP) to solve the all-pairs similarity search problem for time series data. The proposed NMP inherits the features from two state-of-the art algorithms: similarity time-series automatic multivariate prediction (STAMP), and short text online microblogging protocol (STOMP). The proposed algorithm caches the output in an easy-to-access fashion for single- and multidimensional data. The proposed NMP algorithm can be used on large data sets and generates approximate solutions of high quality in a reasonable time. The proposed NMP can also handle several data mining tasks. It is implemented on a Python platform. To determine its effectiveness, it is compared with the state-of-the-art matrix profile algorithms i.e., STAMP and STOMP. The results confirm that the proposed NMP provides higher accuracy than the compared algorithms.

The anomaly detection task is very important in computer science. And there are a lot of anomaly detection methods. Different from some thresholding methods, some unsupervised methods could make us get more accurate and faster result, which is the object of the project. In this paper, I tried to use EWMA and some other methods in two datasets: Webank time consuming indicators dataset and AIOps Challenge dataset. The paper consists nine parts: background of the project, related work, description of algorithms, implementation details, experimental setup and data sets used, experimental results and discussion, future directions, reference and meeting notes.

Divergences have become a very useful tool for measuring similarity (or dissimilarity) between probability distributions. Depending on the field of application a more appropriate measure may be necessary. In this paper we introduce a family of divergences we call gamma-divergences. They are based on the convexity property of the functions that generate them. We demonstrate that these divergences verify all the usually required properties, and we extend them to weighted probability distribution. We investigate their properties in the context of kernel theory. Finally, we apply our findings to the analysis of simulated and real time series.

While increasing empirical evidence suggests that global time series forecasting models can achieve better forecasting performance than local ones, there is a research void regarding when and why the global models fail to provide a good forecast. This paper uses anomaly detection algorithms and Explainable Artificial Intelligence (XAI) to answer when and why a forecast should not be trusted. To address this issue, a dashboard was built to inform the user regarding (i) the relevance of the features for that particular forecast, (ii) which training samples most likely influenced the forecast outcome, (iii) why the forecast is considered an outlier, and (iv) provide a range of counterfactual examples to understand value changes, in the feature vector or the predicted value, can lead to a different outcome. Moreover, a modular architecture and a methodology were developed to iteratively remove noisy data instances from the train set, to enhance the overall global time series forecasting model performance. Finally, to test the effectiveness of the proposed approach, it was validated on two publicly available real-world datasets.

This article presents the results of reviewing the predictive capacity of Google trends for national elections in Chile. The electoral results of the elections between Michelle Bachelet and Sebastián Piñera in 2006, Sebastián Piñera and Eduardo Frei in 2010, Michelle Bachelet and Evelyn Matthei in 2013, Sebastián Piñera and Alejandro Guillier in 2017, and Gabriel Boric and José Antonio Kast in 2021 were reviewed. The time series analysed were organised on the basis of relative searches between the candidacies, assisted by R software, mainly with the gtrendsR and forecast libraries. With the series constructed, forecasts were made using the ARIMA technique to check the weight of one presidential option over the other. The ARIMA analyses were performed on 3 ways of organising the data: the linear series, the series transformed by moving average and the series transformed by Hodrick-Prescott. The result indicates that the method offers optimal pre-dictive ability.

Landslide is a sliding movement of rock mass, debris and soil along the slope under the action of gravity. Small Baseline Subset (SBAS) is an established method for the investigation and monitoring of landslide moments. This study focuses on monitoring the long-temporal displacement of mountainous terrain in Danba County, Sichuan Province via SBAS technique, based on 31 scenes of L-band ALOS/PALSAR data from Feb. 2007 to Oct. 2010.The results show that the largest velocity rates in LOS direction are ±120 mm/yr and maximum accumulated displacement is up to -300, which indicates fast movement of the mountainous terrain during the observation time. These results get good consistency against the results of previous study. This demonstrates the strong potential of SBAS technique for monitoring the landslides geohazard.

Reservoir computing is a machine learning method that uses the response of a dynamical system to a certain input in order to solve a task. As the training scheme only involves optimising the weights of the responses of the dynamical system, this method is particularly suited for hardware implementation. Furthermore, the inherent memory of dynamical systems which are suitable for use as reservoirs mean that this method has the potential to perform well on time series prediction tasks, as well as other tasks with time dependence. However, reservoir computing still requires extensive task dependent parameter optimisation in order to achieve good performance. We demonstrate that by including a time-delayed version of the input for various time series prediction tasks, good performance can be achieved with an unoptimised reservoir. Furthermore, we show that by including the appropriate time-delayed input, one unaltered reservoir can perform well on six different time series prediction tasks at a very low computational expense. Our approach is of particular relevance to hardware implemented reservoirs, as one does not necessarily have access to pertinent optimisation parameters in physical systems but the inclusion of an additional input is generally possible.

Traditional applications of Interferometric Synthetic Aperture Radar (InSAR) data involved inverting an interferogram stack to determine the average displacement velocity. While this approach has useful applications in continuously deforming regions, much information is lost by simply fitting a line through the time series. Thanks to regular acquisitions across most of the the world by the ESA Sentinel-1 satellite constellation, we are now in a position to explore opportunities for near-real time deformation monitoring. In this paper we present a simple statistical approach for detecting offsets and gradient changes in InSAR time series. Our key assumption is that 5 years of Sentinel-1 data is sufficient to calculate the population standard deviation of the detection variables. Our offset detector identifies statistically significant peaks in the first, second and third difference series. The gradient change detector identifies statistically significant movements in the second derivative series. We exploit the high spatial resolution of Sentinel-1 data and the spatial continuity of geophysical deformation signals to filter out false positive detections that arise due to signal noise. When combined with near-real time processing of InSAR data these detectors, particularly the gradient change, could be used to detect incipient ground deformation associated with geophysical phenomena, for example from landslides or volcanic eruptions.

The study of order patterns of three equally spaced values $x_t,x_{t+d},x_{t+2d}$ in a time series is a powerful tool. The lag $d$ is changed in a wide range so that differences of frequencies of order patterns become autocorrelation functions. Similar to a spectrogram in speech analysis, four ordinal autocorrelation functions are used to visualize big data series, as for instance heart and brain activity over many hours. The method applies to real data without preprocessing, outliers and missing data do not matter. On the theoretical side, we study properties of order correlation functions and show that the four autocorrelation functions are orthogonal in a certain sense. An analysis of variance of a modified permutation entropy can be performed with four variance components associated with the functions.

With the great development of intelligent transportation systems (ITS), travel time prediction has attracted the attentions of many researchers and a large number of prediction methods have been developed. However, as an unavoidable topic, the predictability of travel time series is the basic premise for travel time prediction has received less attention than the methodology. Based on the analysis of the complexity of travel time series, this paper defines travel time predictability to express the probability of correct travel time prediction and proposes an entropy-based method to measure the upper bound of travel time predictability. Multiscale entropy is employed to quantify the complexity of travel time series, and the relationships between entropy and the upper bound of travel time predictability are presented. Empirical studies are made with vehicle trajectory data in an express road section. The effectiveness of time scales, tolerance, and series length to entropy and travel time predictability are analysis, and some valuable suggestions about the accuracy of travel time predictability are discussed. Finally, the comparisons between travel time predictability and actual prediction results from two prediction models, ARIMA and BPNN, are conducted. Experimental results demonstrate the validity and reliability of the proposed travel time predictability.

Machine learning methods have widely been applied to detect anomalies in machine and cutting tool behavior during lathe or milling. However, detecting anomalies in the workpiece itself did not get the same attention by researchers. That is why in this article, the authors present a pub-licly available multivariate time series dataset which was recorded during milling of 16MnCr5. Due to artificially introduced, though realistic anomalies in the workpiece the dataset can be ap-plied for anomaly detection. By using a convolutional autoencoder as a first model good results in detecting the location of the anomalies in the workpiece were achieved. Furthermore, milling tools with two different diameters where used which led to a dataset eligible for transfer learn-ing. The objective of this article is to provide researchers with a real-world time series dataset of the milling process which is suitable for modern machine learning research topics like anomaly detection and transfer learning.

The removal of non-climatic biases, so-called inhomogeneities, from long climatic records needs sophistically developed statistical methods. One principle is that usually the differences between a candidate series and its neighbour series are analysed instead of directly the candidate series, in order to neutralize the possible impacts of regionally common natural climate variation on the detection of inhomogeneities. In most homogenization methods, two main kinds of time series comparisons are applied, i.e. composite reference series or pairwise comparisons. In composite reference series the inhomogeneities of neighbour series are attenuated by averaging the individual series, and the accuracy of homogenization can be improved by the iterative improvement of composite reference series. By contrast, pairwise comparisons have the advantage that coincidental inhomogeneities affecting several station series in a similar way can be identified with higher certainty than with composite reference series. In addition, homogenization with pairwise comparisons tends to facilitate the most accurate regional trend estimations. A new time series comparison method is presented here, which combines the use of pairwise comparisons and composite reference series in a way that their advantages are unified. This time series comparison method is embedded into the ACMANT homogenization method, and tested in large, commonly available monthly temperature test datasets.

This study compared the performance of five Family Generalized Auto-Regressive Conditional Heteroscedastic (fGARCH) models (sGARCH, gjrGARCH, iGARCH, TGARCH and NGARCH) in the presence of high positive autocorrelation. To achieve this, financial time series was simulated with autocorrelated coefficients as ρ = (0.8, 0.85, 0.9, 0.95, 0.99), at different time series lengths (as 250, 500, 750, 1000, 1250, 1500) and each trial was repeated 1000 times carried out in R environment using rugarch package. And the performance of the preferred model was judged using Root Mean Squared Error (RMSE) and Mean Absolute Error (MAE). Results from the simulation revealed that these GARCH models performances varies with the different autocorrelation values and at different time series lengths. But in the overall, NGARCH model dominates with 62.5% and 59.3% using RMSE and MAE respectively. We therefore recommended that investors, financial analysts and researchers interested in stock prices and asset return should adapt NGARCH model when there is high positive autocorrelation in the financial time series data.

In this paper, we study the usage of machine learning models for sales time series forecasting. The effect of machine learning generalization has been considered. A stacking approach for building regression ensemble of single models has been studied. The results show that using stacking technics, we can improve the performance of predictive models for sales time series forecasting.

The morphology and the activity of a submarine caldera, Avyssos, at the northern part of Nisyros volcano in the South Aegean Sea (Greece), has been studied by means of remotely operated underwater vehicle dives. The recorded time series of temperature and conductivity over the submarine volcano have been analyzed in terms of the Generalized Moments Method. The findings of the mathematical analysis shed light on the volcanic activity, but also on the morphology (shape) of the submarine volcano. The conductivity time series indicates the volcano is at rest in agreement with other types of observations. On the other hand, temperature fluctuations, which in general describe a multifractal process, show that the submarine caldera operates as an open system that interacts with its surroundings. This type of analysis can be used as an indicator for the state of activity and the morphological structure (closed or open system) of a submarine volcano.

We define a new seasonal forecasting method based on fuzzy transforms. We use the best interpolating polynomial for extracting the trend of the time series and generate the inverse fuzzy transform on each seasonal subset of the universe of discourse for predicting the value of a an assigned output. Like first example, we use the daily weather dataset of the municipality of Naples (Italy) starting from data collected from 2003 till to 2015 making predictions on the following outputs: mean temperature, max temperature and min temperature, all considered daily. Like second example, we use the daily mean temperature measured at the weather station “Chiavari Caperana” in the Liguria Italian Region. We compare the results with our method, the average seasonal variation, ARIMA and the usual fuzzy transforms concluding that the best results are obtained under our approach in both examples.

In this paper, we study the asymptotic behavior of the sequential empirical process and the sequential empirical copula process, both constructed from residuals of multivariate stochastic volatility models. Applications for the detection of structural changes and specification tests of the distribution of innovations are discussed. It is also shown that if the stochastic volatility matrices are diagonal, which is the case if the univariate time series are estimated separately instead of being jointly estimated, then the empirical copula process behaves as if the innovations were observed; a remarkable property. As a by-product, one also obtains the asymptotic behavior of rank-based measures of dependence applied to residuals of these time series models.

The alteration of natural land cover to impervious surfaces during development increases stormwater runoff. Stormwater Control Measures (SCMs) are used to manage water quantity and enhance water quality by restoring the hydrologic cycle altered by development. Often, SCMs have an outflow pipe to handle overflows or to manage the release of water detained when infiltration is not possible. Traditionally, these are static controls (e.g. a small orifice is used to restrict the volume of outflow), however, these systems can be improved by instituting real-time controls (RTC). RTC improve the functionality of SCMs by dynamically controlling outflows to adjust to environmental conditions. A major impediment to the widespread implementation of RTC is the high cost of installation and operation. This study utilized machine learning methods to develop a forecasting approach for the implementation of low-cost RTC that were implemented on a programmable gate of the outlet structure of a multi-stage basin in southeastern Pennsylvania. The goals were to decrease the peak flow exiting the basin during rain events, increase the volume of water detained, decrease the number of overtopping events, maintain healthy vegetation in the basin, and protect the downstream vegetation from erosion. Multiple popular data science algorithms were evaluated including multiple linear regression and long short-term memory. These algorithms were used with a dataset, which consisted of four years of historical sensor data, collected in 5-minute intervals, to train models to predict water levels to optimize operations. The accuracy of 30 models with three different methods of handling missing values were compared. A long short-term memory model configured with a 30-minute lead time produced the best results. Having an approximate same lag time of 30 minutes for the contributing drainage area of the SCM provided a sufficient RTC functioning period to improve the performance of the outlet structure.

The futures market plays a significant role in hedging and speculating by investors. Although various models and instruments are developed for real-time trading, it is difficult to realize profit by processing and trading a vast amount of real-time data. This study proposes a real-time index futures trading strategy that uses the pattern of KOSPI 200 index futures time series data. We construct a pattern matching trading system (PMTS) based on a dynamic time warping algorithm that recognizes patterns of market data movement in the morning and determines the afternoon's clearing strategy. We adopt 13 and 27 representative patterns and conduct simulations with various ranges of parameters to find optimal ones. Our experimental results show that the PMTS provides stable and effective trading strategies with relatively low trading frequencies. Investor communities that have sustained financial markets are able to make more efficient investments by using the PMTS. In this sense, the system developed in this paper is a sustainable investment technique and helps financial markets achieve efficient sustainability.

The principal aim of a spectral observer is twofold: the reconstruction of a signal of time via state estimation and the decomposition of such a signal into the frequencies that make it up. A spectral observer can be catalogued as an online algorithm for time-frequency analysis because is a method that can compute on the fly the Fourier transform (FT) of a signal, without having the entire signal available from the start. In this regard, this paper presents a novel spectral observer with an adjustable constant gain for reconstructing a given signal by means of the recursive identification of the coefficients of a Fourier series. The reconstruction or estimation of a signal in the context of this work means to find the coefficients of a linear combination of sines a cosines that fits a signal such that it can be reproduced. The design procedure of the spectral observer is presented along with the following applications: (1) the reconstruction of a simple periodical signal, (2) the approximation of both a square and a triangular signal, (3) the edge detection in signals by using the Fourier coefficients, (4) the fitting of the historical Bitcoin market data from 2014-12-01 to 2018-01-08 and (5) the estimation of a input force acting upon a Duffing oscillator. To round out this paper, we present a detailed discussion about the results of the applications as well as a comparative analysis of the proposed spectral observer vis-à-vis the Short Time Fourier Transform (STFT), which is a well-known method for time-frequency analysis.

We present a custom platform that integrates data from several sensors measuring synchronously different variables of the beehive and wirelessly transmits all measurements to a cloud server. There is a rich literature on beehive monitoring. The choice of our work is not to use ready platforms such as Arduino and Raspberry Pi and to present a low cost and power solution for long term monitoring. We integrate sensors that are not limited to the typical toolbox of beehive monitoring such as gas, vibrations and bee counters. The synchronous sampling of all sensors every 5 minutes allows us to form a multivariable timeseries that serves two-ways: a) it provides immediate alerting in case a measurement exceeds predefined boundaries that are known to characterize a healthy beehive, and b) based on historical data predict future levels that are correlated with hive’s health. Finally, we demonstrate the benefit of using additional regressors in the prediction of the variables of interest. The database, the code and a video of the vibrational activity of two months are made open to the interested readers.

Increasingly measured data in the context of smart cities can be used to develop new and innovative business models to increase efficiency and the value of life. A time-series classification algorithm can support to automatize many different processes such as forecasting services. In order to ensure data security and privacy, Federated Learning trains a global model collaboratively on multiple clients. Having different data-distributions and data-quantities across participating clients, neural networks suffer from slow convergence and overfitting. Based on different data-distributions, data-quantities and number of clients, we develop and evaluate different data-clustering strategies to update global model weights in comparison to the state of the art. We use public time-series data, generate various synthetic datasets and train a Relational-Regularized Autoencoder for classification purposes. Our results show an improvement of model performance concerning generalization.

Vegetation status assessment is crucial for agricultural monitoring and management. Vegetation indices derived from high resolution image time series can be used to derive key phenological parameters for annual crops. In this work, we propose a procedure for the estimation of these parameters and their associated uncertainties. The approach uses Bayesian inference through Markov Chain Monte Carlo in order to obtain the full joint posterior distribution of the phenological parameters given the satellite observations. The proposed algorithm is quantitatively validated on synthetic data. Its use on real data is presented together with an application to real-time within season estimation allowing for phenology forecasting.

Background: Jerusalem artichoke (Helianthus tuberosus L.) is tolerant to salinity stress and has high economic value. The salt tolerance mechanisms of Jerusalem artichoke are still unclear. Especially in the early stage of Jerusalem artichoke exposure to salt stress, the plant physiology, biochemistry and gene transcription are likely to undergo large changes. Elucidating these changes may be of great significance to understanding the salt tolerance mechanisms of it. Results: We obtained high-quality transcriptome from leaves and roots of Jerusalem artichoke exposed to salinity (300 mM NaCl) for 0 h, 6 h, 12 h, 24 h and 48 h, with 150,129 unigenes and 9023 DEGs (Differentially Expressed Genes). The RNA-seq data were clustered into time-dependent groups (nine clusters each in leaves and roots); gene functions were distributed evenly among the groups convergence. KEGG enrichment analysis showed the genes related to plant hormone signal transduction were enriched in almost all treatment comparisons. Under salt stress, genes belongs to PYL (abscisic acid receptor PYR / PYL family), PP2C (Type 2C protein phosphatases), GH3 (Gretchen Hagen3), ETR (ethylene receptor), EIN2/3 (ethylene-insensitive protein 2/3), JAZ (Genes such as jasmonate ZIM-domain gene) and MYC2 (Transcription factor MYC2) had extremely similar expression patterns. The results of qPCR of 12 randomly selected genes confirmed the accuracy of RNA-seq. Conclusions: Under the impact of high salinity (300mM) environment, Jerusalem artichoke in the seedling stage was difficult to survive for a long time, and the phenotype was severe in the short term. Based on the expression of genes on the time scale, we found that the distribution of gene functions in time is relatively even. Upregulation of the phytohormone signal transduction had a crucial role in the response of Jerusalem artichoke seedlings to salt stress, the genes of abscisic acid, auxin, ethylene, and jasmonic acid had the most obvious change pattern.

It is a common practice to detect outliers in a financial time series in order to avoid the adverse effect of additive outliers. This paper investigated the performance of GARCH family models (sGARCH; gjrGARCH; iGARCH; TGARCH and NGARCH) in the presence of different sizes of outliers (small, medium and large) for different time series lengths (250, 500, 750, 1000, 1250 and 1500) using root mean square error (RMSE) and mean absolute error (MAE) to adjudge the models. In a simulation iteration of 1000 times in R environment using rugarch package, results revealed that for small size of outliers, irrespective of the length of time series, iGARCH dominated, for medium size of outliers, it was sGARCH and gjrGARCH that dominated irrespective of time series length, while for large size of outliers, irrespective of time series length, gjrGARCH dominated. The study further leveled that in the presence of additive outliers on time series analysis, both RMSE and MAE increased as the time series length increased.

A continuous Deep Belief Network (cDBN) with two hidden layers is proposed in this paper, focusing on the problem of weak feature learning ability when dealing with continuous data. In cDBN, the input data is trained in an unsupervised way by using continuous version of transfer functions, the contrastive divergence is designed in hidden layer training process to raise convergence speed, an improved dropout strategy is then implemented in unsupervised training to realize features learning by de-cooperating between the units, and then the network is fine-tuned using back propagation algorithm. Besides, hyper-parameters are analysed through stability analysis to assure the network can find the optimal. Finally, the experiments on Lorenz chaos series, CATS benchmark and other real world like CO2 and waste water parameters forecasting show that cDBN has the advantage of higher accuracy, simpler structure and faster convergence speed than other methods.

Introduction Pharmacovigilance may detect safety issues after marketing of medications, and this can result in regulatory action such as direct healthcare professional communications (DHPC). DHPC can be effective in changing prescribing behaviour, however the extent to which prescribers vary in their response to DHPC is unknown. This study aims to explore changes in prescribing and prescribing variation among GP practices following a DHPC on the safety of mirabegron, a medication to treat overactive bladder (OAB). Methods This is an interrupted time series study of English GP practices from 2014-2017. NHS Digital provided monthly statistics on aggregate practice-level prescribing and practice characteristics (practice staff and registered patient profiles, Quality & Outcomes Framework indicators, and deprivation of the practice area). The primary outcome was monthly mirabegron items as a percentage of all OAB drug items. The exposure was a DHPC issued by the European Medicines Agency in September 2015. Variation between practices in mirabegron prescribing before and after the DHPC was assessed using the systematic component of variation (SCV). Multilevel segmented regression with random effects quantified the change in level and trend of prescribing after the DHPC. Practice characteristics were assessed for their association with a reduction in prescribing following the DHPC. Results This study included 7,408 practices. During September 2015, 88.9% of practices prescribed mirabegron and mirabegron composed a mean of 8.2% (SD 6.8) of OAB items. Variation between practices was classified as very high and the median SCV did not change significantly (p=0.11) in the 6 months after the September 2015 DHPC (12.4) compared to before (11.6). Before the DHPC, there was a monthly trend of 0.294 (95%CI, 0.287, 0.301) percentage points increase in mirabegron percentage. There was no significant change in the month immediately after the DHPC (-0.023, 95% CI -0.105 to 0.058) however there was a significant reduction in trend (-0.036, 95% CI -0.049 to -0.023). Higher numbers of registered patients and patients aged ≥65 years, and practice area deprivation were associated with having a significant decrease in level and slope of mirabegron prescribing post-DHPC. Conclusion Variation in mirabegron prescribing was high over the study period and did not change substantively following the DHPC. There was no immediate prescribing change post-DHPC, although the monthly growth did slow. Knowledge of the degree of variation in and determinants of response to safety communications may allow those that do not change prescribing to be provided with additional supports.

The relationship between economic growth, growth volatility and financial sector development continues to attract attention in the theoretical and empirical literature. Over time, some studies hypothesize that finance has a causal linear relationship with growth. Recently several other authors contradict this claim and argue that the relationship that exists between finance and growth is nonlinear. We investigate these claims for Nigeria for the period between 1970 and 2015, using semi-parametric econometric methods, Hansen sample splitting techniques and threshold estimator. We observed no evidence of ‘Too much finance’ as claimed by many researchers in recent times. We show that the relationship between financial development and economic growth is U-shaped. This is equally true for the relationship between financial development and growth volatility. We also discuss policy implications of our findings and recommend financial innovations and decentralization of stock exchanges to boost access to financial services, in addition, improved regulation to enhance financial market efficiency.

We present a novel method for interpolating univariate time series data. The proposed method combines multi-point fractional Brownian bridges, a genetic algorithm, and Takens’ theorem for reconstructing a phase space from univariate time series data. The basic idea is to first generate a population of different stochastically interpolated time series data, and secondly, to use a genetic algorithm to find the pieces in the population which generate the smoothest reconstructed phase space trajectory. A smooth trajectory curve is hereby found to have a low variance of second derivatives along the curve. For simplicity, we refer to the developed method as PhaSpaSto-interpolation, which is an abbreviation for phase-space-trajectory-smoothing stochastic interpolation. The proposed approach is tested and validated with a univariate time series of the Lorenz system, five non-model data sets and tested against a cubic spline interpolation and a linear linear interpolation. We find that the criterion for smoothness guarantees low errors on known model and non-model data. Finally, we interpolate the discussed non-model data sets, and show the corresponding improved phase space portraits. The proposed method is useful for interpolating low-sampled time series data sets for, e.g., machine learning, regression analysis, or time series prediction approaches.

This article presents the implementation of a model to estimate the future price of commodities in the Brazilian market from time series of short-term technical evaluation. For this, data from two databases were used, one referring to the foreign market (opening values, maximum, minimum, closing, closing adjustment and volume) and the other, from the Brazilian market (the price of the day), considering commodities, sugar, cotton, corn, soybean and wheat. Subsequently, the technical indicators were calculated from the TA-Lib technical analysis library. Pearson’s correlation coefficient was applied, records with low correlation were removed, and then the database was consolidated. From the pre-processed data, Long Short-Term Memory (LSTM) recurrent neural networks were used to perform data prediction at the one and three day interval. These models were evaluated using the mean square error (MSE), obtaining results between 0.00010 and 0.00037 on test data one day ahead, and from 0.00017 to 0.00042 three days ahead. However, based on the results obtained, it was observed that the developed model obtained a promising forecasting performance for all the commodities evaluated. As a main contribution, there is the consolidation of databases that can be used in future scientific research. Furthermore, based on its interpretation, it can assist in decision making regarding the buying and selling of commodities to increase financial gains.

To improve the monitoring of the electrical power grid it is necessary to evaluate the influence of contamination in relation to leakage current and its progression to a disruptive discharge. In this paper, insulators were tested in a saline chamber to simulate the increase of salt contamination on their surface, and to evaluate the supportability of these components. From the time series forecasting of the leakage current, it is possible to evaluate the development of the fault before a flashover occurs. Choosing which method to use is always a difficult task since some models may have higher computational effort. In this paper, for a complete evaluation, the long short-term memory (LSTM), group method of data handling (GMDH), adaptive neuro-fuzzy inference system (ANFIS), bootstrap aggregation (bagging), sequential learning (boosting), random subspace, and stacked generalization (stacking) ensemble learning models are analyzed. A review and comparison of these well-established methods for time series forecasting is performed. From the results of the best structure of the model, the hyperparameters are evaluated and the Wavelet transform is used to obtain an enhanced model.

Photoplethysmography is a widely used technique to noninvasively assess heart rate, blood pressure, and oxygen saturation. This technique has a large potential for further applications, for example in the field of physiological and mental health monitoring. However, advanced applications of photoplethysmography have been hampered by the lack of accurate and reliable methods to analyze the characteristics of the complex nonlinear dynamics of the photoplethysmogram. Methods of nonlinear time series analysis may be used to estimate the dynamical characteristics of the photoplethysmogram but they are highly influenced by the length of the time series, which is often limited in practical photoplethysmography applications. The aim of this study was to evaluate the error in the estimation of the dynamical characteristics of the photoplethysmogram associated to the limited length of the time series. The dynamical properties were evaluated using recurrence quantification analysis, and the estimation error was computed as a function of the length of the time series. Results demonstrated that properties such as de-terminism and entropy can be estimated with an error lower than 1% even for short photople-thysmogram recordings. Additionally, the lower limit for the time series length to estimate the average prediction time was computed.

This paper describes the successful application of the Transformer model used in the natural language processing and vision tasks as a means of processing the time series of signals from gyroscope and accelerometer sensors for the classification of human activities. The Transformer model is based on deep neural networks with many layers which can generalize well on signals. All measured signals come from a smartphone placed in a waist bag. Activity prediction is sequence-to-sequence, each time step of the signal is assigned a designation of the performed activity. Emphasis is placed on attention mechanisms, which express individual dependencies between signal values within a time series. In comparison with another recent result, the recognition precision was improved from 89.67 percent to 99.2 percent. The transformer model should in the future be included among the top options in machine learning methods for human activity recognition.

This study investigates the impact of Foreign Direct Investment (FDI) on economic growth and examines the causality between FDI and economic growth in Bangladesh during 1972-2013. Gross Domestic Product (GDP), export performance (EXP), Foreign Direct Investment (FDI), and Gross Fixed Capital Formation (GFCF) are considered to capture the objective of the study. The study methodology includes some systematic steps. As the data used in the study is time-series in nature, the author employs unit root tests, and in this case, Augmented Dickey-Fuller (ADF) and Phillips-Perron (PP) tests are used. Then Johansen’s cointegration test, Granger causality test, regression with Newey-West Standard Error and Vector Error Correction Model (VECM) are applied. By using the ADF and PP test the study reveals that the variables of four-time series are integrated of I (1) i.e. they are stationary at first difference. Regression analysis result demonstrates that FDI has a positive effect on economic growth. The Granger Causality test discloses that there is a unidirectional relationship between FDI and economic growth. But the VECM estimation finds that in the long run FDI negatively affects economic growth.

Air sampling using vortex air samplers combined with species specific amplification of pathogen DNA, was carried out over two years in four or five locations in the Salinas Valley of California. The resulting time series data for the abundance of pathogen DNA trapped per day displayed complex dynamics with features of both deterministic (chaotic) and stochastic dynamics. Methods of nonlinear time series analysis developed for the reconstruction of low dimensional attractors provided new insights into the complexity of the pathogen abundance data, but also indicated that practicality may limit the capacity for definitively classifying the dynamics of air borne plant pathogen inoculum. Over the two years of the study five location/year combinations were classified as having stochastic linear dynamics and four were not. Calculation of entropy values for either the number of pathogen DNA copies or for a binary string indicating the pathogen abundance data were increasing or not, revealed (1) some robust differences in the dynamics between seasons that were not obvious in the time series data themselves, and also (2) that the series were almost all at their theoretical maximum entropy value when considered from the simple perspective of whether instantaneous change along the sequence is positive or not.

We discuss methods for calculation of critical indices and amplitudes from the perturbative expansions. Several important examples of the Stokes flow through 2D and 3D channels are brought up. Power series for the permeability derived for small values of amplitude are employed to calculation of various critical exponents in the regime of large amplitudes. Special nonlinear approximations valid for arbitrary values of the wave amplitude are derived from the expansions. The technique developed for critical phenomena is applied then for relaxation phenomena. The concept of time-translation invariance is discussed, its spontaneous violation and restoration considered. Emerging probabilistic patterns correspond to a local breakdown of time-translation invariance. Their evolution leads to the timetranslation symmetry complete (or partial) restoration. We estimate typical time extent, amplitude and direction for such restorative process. The new technique is based on explicit introduction of origin in time. After some transformations we come to the exponential and generalized, exponential-type solution with explicit finite time scale, which was only implicit in initial parametrization with polynomial approximation. The concept of crash as a relaxation phenomenon, consisting of time-translation invariance breaking and restoration, is put forward. COVID-19 related mini-crash in the time series for Shanghai Composite is discussed as an illustration.

Terrestrial laser scanning (TLS) has been adopted as a feasible technique to digitize trees and forest stands, providing accurate information on tree and forest structural attributes. However, there is limited understanding on how a variety of forest structural changes can be quantified using TLS in boreal forest conditions. In this study, we assessed the accuracy and feasibility of TLS in quantifying changes in the structure of boreal forests. We collected TLS data and field reference from 37 sample plots in 2014 (T1) and 2019 (T2). Tree stems typically have planar, vertical, and cylindrical characteristics in a point cloud, and thus we applied surface normal filtering, point cloud clustering, and RANSAC-cylinder filtering to identify these geometries and to characterize trees and forest stands at both time points. The results strengthened the existing knowledge that TLS has the capacity to characterize trees and forest stands in space and showed that TLS could characterize structural changes in time in boreal forest conditions. Root-mean-square-errors (RMSEs) in the estimates for changes in the tree attributes were 0.99-1.22 cm for diameter at breast height (Δdbh), 44.14-55.49 cm2 for basal area (Δg), and 1.91-4.85 m for tree height (Δh). In general, tree attributes were estimated more accurately for Scots pine trees, followed by Norway spruce and broadleaved trees. At the forest stand level, an RMSE of 0.60-1.13 cm was recorded for changes in basal area-weighted mean diameter (ΔDg), 0.81-2.26 m for changes in basal area-weighted mean height (ΔHg), 1.40-2.34 m2/ha for changes in mean basal area (ΔG), and 74-193 n/ha for changes in the number of trees per hectare (ΔTPH). The plot-level accuracy was higher in Scots pine-dominated sample plots than in Norway spruce-dominated and mixed-species sample plots. TLS-derived tree and forest structural attributes at time points T1 and T2 differed significantly from each other (p < 0.05). If there was an increase or decrease in dbh, g, h, height of the crown base, crown ratio, Dg, Hg, or G recorded in the field, a similar outcome was achieved by using TLS. Our results provided new information on the feasibility of TLS for the purposes of forest ecosystem growth monitoring.

Modern energy systems collect high volumes of data that can provide valuable information about energy consumption. Electric companies can now use historical data to make informed decisions on energy production by forecasting the expected demand. Many deep learning models have been proposed to deal with these type of time series forecasting problems. Deep neural networks, such as recurrent or convolutional, can automatically capture complex patterns in time series data and provide accurate predictions. In particular, Temporal Convolutional Networks (TCN) are a specialised architecture that has advantages over recurrent networks for forecasting tasks. TCNs are able to extract long-term patterns using dilated causal convolutions and residual blocks, and can also be more efficient in terms of computation time. In this work, we propose a TCN-based deep learning model to improve the predictive performance in energy demand forecasting. Two energy-related time series with data from Spain have been studied: the national electric demand, and the power demand at charging stations for electric vehicles. An extensive experimental study has been conducted, involving more than 1900 models with different architectures and parametrisations. The TCN proposal outperforms the forecasting accuracy of Long Short-Term Memory (LSTM) recurrent networks, which are considered the state-of-the-art in the field.

Inland lake variations are considered sensitive indicators of global climate change. However, human activity is playing as a more and more important role in inland lake area variations. Therefore, it is critical to identify whether anthropogenic activity or natural event is playing as the dominant factor in inland lake surface area change. In this study, we proposed a Douglas-Peucker simplification algorithm and bend simplification algorithm combined method to locate major lake surface area disturbances; these disturbances were then characterized to extract the time series change features according to documented records; and the disturbances were finally classified into anthropogenic or natural. We took the nine lakes in Yunnan Province as test sites, a 31 years long (from 1987 to 2017) time series Landsat TM/OLI images and HJ-1A/1B used as data sources, the official records was used as references to aid the feature extraction and disturbance identification accuracy. Results of our method for both disturbance location and the disturbance identification could be concluded as follows: 1) The method can accurately locate the main lake changing events based on the time series lake surface area curve. The accuracy of this model for segmenting the lake area time series curves in our study area was 95.24%. 2) Our proposed method achieved an overall accuracy of 91.67%, with F-score of 94.67 for anthropogenic disturbances and F-score of 85.71 for natural disturbances. 3) According to our results, lakes in Yunnan Provence, China, have undergone extensive disturbances, and the human-induced disturbances occurred almost twice as often as natural disturbances, indicating intensified disturbances caused by human activities. This inland lake area disturbance identification method is expected to uncover whether a disturbance to inland lake area is human activity-induced or natural event.

Advances in the availability of multitemporal and global built-/human-settlements datasets as derived from Remote Sensing (RS) can now provide globally consistent definitions of “human-settlement” at unprecedented spatial fineness. Yet, these data only provide a time-series of past extents and urban growth/expansion models have not had parallel advances at high-spatial resolution. We present a flexible modelling framework for producing annual built-settlement extents in the near future past last observed extents as provided by RS-based data. Using a random forest and autoregressive temporal models with short time-series of built-settlement extents and subnational level data, we predict annual 100m resolution binary settlement extents five years beyond the last observations. We applied this framework within varying contexts and predicted annual extents from 2010 to 2015. We found that our model framework preformed consistently across all sample countries and, when compared to time-specific imagery, demonstrated the capacity to capture human-settlement missed by the input time-series and validation extents. When comparing building footprints of small settlements to forecast extents, we saw that the modelling framework had a 12 percent increase in ground-truth accuracy. This framework shows promise for predicting near-future settlement extents, and provides a foundation for forecasts further into the future.

The future of the Colorado River water supply (WS) affects millions of people and the U.S. economy. A recent study suggested a cross-basin correlation between the Colorado River and its neighboring Great Salt Lake (GSL). Following that study, the feasibility of using the previously developed multi-year prediction of the GSL water level to forecast the Colorado River WS was tested. Time-series models were developed to predict the changes in WS out to 10 years. Regressive methods and the GSL water level data were used for the depiction of decadal variability of the Colorado River WS. Various time-series models suggest a decline in the 10-year-averaged WS since 2013 before starting to increase around 2020. Comparison between this WS prediction and the WS projection published in a 2012 government report (derived from climate models) reveals a widened imbalance between supply and demand by 2020. Further research to update similar multi-year prediction of the Colorado River WS is needed. Such information could aid in management decision making in the face of future water shortages.

Land Surface Temperature (LST) is a basic parameter in energy exchange between the land and atmosphere and is frequently used in many sciences such as climatology, hydrology, agriculture, ecology, etc. LST time series data have usually deficient, missing and unacceptable data caused by the presence of clouds in images, presence of dust in atmosphere and sensor failure. In this study, Singular Spectrum Analysis (SSA) algorithm was used to resolve the problem of missing and outlier data caused by cloud cover. The region studied in the present research included an image frame of MODIS with horizontal number 22 and vertical number 05 (h22v05). This image involved a large part of Iran and Turkmenistan and Caspian Sea. In this study, MODIS LST sensor (MOD11A1) was used during 2015 with 1×1 Km spatial resolution and day/night LST data (daily temporal resolution). The results of the data quality showed that cloud cover caused 36.37% of missing data in the studied time series with 730 day/night LST images. Further, the results of SSA algorithm in reconstruction of LST images indicated the Root Mean Square Error (RMSE) of 2.95 K between the original and reconstructed data in LST time series in the study region. In general, the findings showed that SSA algorithm using spatio-temporal interpolation in LST time series can be effectively used to resolve the problem of missing data caused by cloud cover.

We propose a novel tensor-based formalism for inferring causal structures from time series. An information theoretical analysis of transfer entropy (TE), shows that TE results from transmission of information over a set of communication channels. Tensors are the mathematical equivalents of these multi-channel causal channels. A multi-channel causal channel is a generalization of a discrete memoryless channel (DMC). We consider a DMC as a single-channel causal channel. Investigation of a system comprising three variables shows that in our formalism, bivariate analysis suffices to differentiate between direct and indirect relations. For this to be true, we have to combine the output of multi-channel causal channels with the output of single-channel causal channels. We can understand this result when we consider the role of noise. Subsequent transmission of information over noisy channels can never result in less noisy transmission overall. This implies that a Data Processing Inequality (DPI) exists for transfer entropy.

Hand and arm gesture recognition using radio frequency (RF) sensing modality proves valuable in man-machine interface and smart environment. In this paper, we use time series analysis method for accurately measuring the similarity of the micro-Doppler (MD) signatures between the training and test data, thus providing improved gesture classification. We characterize the MD signatures by the maximum instantaneous Doppler frequencies depicted in the spectrograms. In particular, we apply the dynamic time warping (DTW) method and compare its performance with that of the long short-term memory (LSTM) network. Both methods take into account the values as well as the temporal evolution and trends of time series data. It is shown that the DTW method achieves high gesture classification rates and is robust to time misalignment.

InSAR and associated analytic methods enable relative surface deformation measurements from low Earth orbit with a potential accuracy of centimeters to millimeters. However, assessing the actual accuracy can be quite difficult. The analytic methods are complicated enough that naive analytic error propagation is infeasible, and, in many settings, InSAR practitioners lack sufficient ground truth to assess results. Phase noise due to partial decorrelation from changes in the scattering properties of the ground is a prominent source of accuracy loss. In this paper we present a method to assess the loss of precision due to this component of phase noise. The proposed method consists of generating synthetic data stacks whose statistical properties match those of the actual input SAR data stacks, and then using the synthetic data for an ensemble calculation. The spread of the results of the ensemble calculation indicates the loss of precision. We show examples of the ensemble analysis at a mining operation in South Africa, and demonstrate the ability to estimate the precision of two InSAR deformation retrieval methods on a point-by-point and epoch-by-epoch basis.

Background: Digital clinical measures collected via various digital sensing technologies such as smartphones, smartwatches, wearables, ingestible, and implantable are increasingly used by individuals and clinicians to capture health outcomes or behavioral and physiological characteristics of individuals. Time series classification (TSC) is very commonly used in modeling digital clinical measures. While deep learning models for TSC are very common and powerful, there exist some fundamental challenges. This review presents non-deep learning models commonly used for time series classification in biomedical applications that achieve high performance. Objective: We performed a systematic review to characterize the techniques used in time series classification of digital clinical measures throughout all stages of data processing and model building. Methods: We conducted a literature search on PubMed, and the Institute of Electrical and Electronics Engineers (IEEE), Web of Science, and SCOPUS databases using a range of search terms to retrieve peer-reviewed articles reporting academic research on digital clinical measures in the five year period between June 2016 and June 2021. We identified and categorized research studies based on the types of classification algorithms and sensor input types. Results: We found 452 papers in total from four different databases: PubMed, IEEE, Web of Science Database, and SCOPUS. After removing duplicates and irrelevant papers, 135 articles remained for detailed review and data extraction. Among these, engineered features using time series methods that were subsequently fed into widely-used machine learning classifiers was the most commonly used technique and also most frequently achieved the best performance metrics (77 out of 135 articles). Statistical modeling (24 out of 135 articles) algorithms were the second most common and also second best classification technique. Wavelet-based classification models (8 out of 135 articles) were also common. Electroencephalogram (29 out of 135 articles) was the most common data type used as an input. Accuracy was the most commonly reported performance metric, with 67.65% of articles reporting on accuracy. In this review paper, we provide summaries of signal pre-processing methods, feature engineering and selection methods, time series models, as well as model interpretations. Importantly, we found that about 50% of the papers only report one performance metric, which may result in a skewed view of overall performance. Conclusion: While high time series classification performance has been achieved in digital clinical, physiological, or biomedical measures, no standard benchmark datasets, modeling methods, or reporting methodology exist. There is no single widely used method for time series model development or feature interpretation– many different methods have proven successful.

Nighttime lights (NTL) are a popular type of data for evaluating economic performance of regions and economic impacts of various shocks and interventions. Several validation studies use traditional statistics on economic activity like national or regional Gross Domestic Product (GDP) as a benchmark to evaluate the usefulness of NTL data. Many of these studies rely on dated and imprecise Defence Meteorological Satellite Program (DMSP) data and use aggregated units such as nation-states or the first sub-national level. Yet applied researchers who draw support from validation studies to justify their use of NTL data as a proxy for economic activity increasingly focus on smaller and lower level spatial units. This study uses a 2001-19 time-series of GDP for over 3100 US counties as a benchmark to examine the usefulness of the recently released version 2 VIIRS nighttime lights (V.2 VNL) products as proxies for local economic activity. Contrasts are made between cross-sectional predictions for GDP differences between areas and time-series predictions of GDP changes within areas. Disaggregated GDP data for various industries are used to examine what types of economic activity are best proxied by NTL data and comparisons are also made with the predictive performance of earlier NTL data products.

Extensive, but remote oil and gas fields of the United States, Canada, and Russia require the construction and operation of extremely long pipelines. Global warming and local heating effects lead to rising soil temperatures and thus a reduction in the sub-grade capacity of the soils; this causes changes in the spatial positions and forms of the pipelines, consequently increasing the number of accidents. Oil operators are compelled to monitor the soil temperature along the routes of the remoted pipelines in order to be able to perform remedial measures in time. They are therefore seeking methods for the analysis of volumetric diagnostic information. To forecast soil temperatures at the different depths we propose compiling a multidimensional dataset, defining descriptive statistics; selecting uncorrelated time series; generating synthetic features; robust scaling temperature series, tuning the additive regression model to forecast soil temperatures.

Lake Malombe fish stocks have been depleted by chronic overfishing. Various management approaches (co-management, command control, and ecosystem-based management to fisheries) have been used to manage the fishery. However, the lack of an accurate predictive model has hampered their success. Therefore, we developed and tested a time series model for Lake Malombe fishery. The seasonal fish biomass and CPUE trends were first observed and both were non-stationary. The second-order differencing was applied to transform the non-stationary data into stationary. Autocorrelation functions (AC), partial autocorrelation function (PAC), and Akaike information criterion (AIC) were estimated, which led to the identification and construction of autoregressive integrated moving average (ARIMA) models, suitable in explaining the time series and forecasting. The results showed that ARIMA (1,2,1) provided a better prediction than its counterparts. The model satisfactorily predicted that by 2032, both fish biomass and CPUE will decrease to 3204.6 tons and 59.672 respectively, signifying the potential threat to Lake Malombe fishery. The model justified the necessity of taking precautionary measures to avoid the total collapse of the fishery.

Monitoring drought impacts in forest ecosystems is a complex process, because forest ecosystems are composed of different species with heterogeneous structural compositions. Even though forest drought status is a key control on the carbon cycle, very few indices exist to monitor and predict forest drought stress. The Forest Drought Indicator (ForDRI) is a new monitoring tool developed by the National Drought Mitigation Center (NDMC) to identify forest drought stress. ForDRI integrates 12 types of data, including satellite, climate, evaporative demand, ground water, and soil moisture, into a single hybrid index to estimate tree stress. The model uses Principal Component Analysis (PCA) to determine the contribution of each input variable based on its covariance in the historical records (2003–2017). A 15-year time series of 780 ForDRI maps at a weekly interval were produced. The ForDRI values at a 12.5km spatial resolution were compared with normalized weekly Bowen ratio data, a biophysically based indicator of stress, from nine AmeriFlux sites. There were strong and significant correlations between Bowen ratio data and ForDRI at sites that had experienced intense drought. In addition, tree ring annual increment data at eight sites in four eastern U.S. national parks were compared with ForDRI values at the corresponding sites. The correlation between ForDRI and tree ring increments at the selected eight sites during the summer season ranged between 0.46 and 0.75. Generally, the correlation between the ForDRI and normalized Bowen ratio or tree ring increment are reasonably good and indicate the usefulness of the ForDRI model for estimating drought stress and providing decision support on forest drought management.

For more than ten years, Central Chile faces drought conditions, which impact crop production and quality, increasing food security risk. Under this scenario, implementing management practices that allow increasing water use efficiency is urgent. The study was carried out in kiwifruit trees, located in the O’Higgins region, Chile; for season 2018-2019 and 2019-2020. We evaluate nine vegetation indices in the VNIR and SWIR regions derived from Sentinel-2 (A/B) satellites to know how much variability in the canopy water status could explain. Over the study's site were installed sensors that continuously measure the leaf's turgor pressure (Yara Water-Sensor). A strong correlation between turgor pressure and vegetation indices was obtained with the Spearman's rho coefficient ($\rho$). However, the NIR range's indices were influenced by the vegetative development of the crop rather than its water status. Red-edge showed better performance as the vegetative growth did not affect it. It is necessary to expand the study to consider higher variability in kiwifruit's water conditions and incorporate the sensitivity of different wavelengths.

Many clustering approaches succeed in pattern segmentation in many applications. This unsupervised segmentation should be effective to reduce an expert labelling time: i.e, they must be able to detect the number of patterns and identify them in a sequence or map with the right cuts. Several direct and hierarchical clustering approaches are compared for this task. A divisive spectral clustering architecture with a no-cut criteria is also proposed. This new algorithm achieves promise segmentation of spatial UCI databases and marine time series compared to other approaches.

Background: Coronavirus disease 2019 (COVID-19) pandemic are causing significant damages to many nations. For mitigating its risk, Japan’s Prime Minister called on all elementary, junior high and high schools nationwide to close beginning March 1, 2020. However, its effectiveness in decreasing disease burden has not been investigated. Methods: We used daily data on the report of COVID-19 and coronavirus infection incidence in Japan until March 31, 2020. Time series analysis were conducted using Bayesian method. Local linear trend models with interventional effect were constructed for number of newly reported cases of COVID-19, including asymptomatic infections. We considered that the effects of intervention start to appear 9 days after the school closure; i.e., on March 9. Results: The intervention of school closure did not appear to decrease the incidence of coronavirus infection. If the effectiveness of school closure began on March 9, mean coefficient α for effectiveness of the measure was calculated to be 0.08 (95% credible interval -0.36 to 0.65), and the actual reported cases were more than predicted, yet with rather wide credible interval. Sensitivity analyses using different dates also showed similar results. Conclusions: School closure carried out in Japan did not show the effectiveness to mitigate the transmission of novel coronavirus infection.

Mapping forest composition using multiseasonal optical time series is still challenging. Highly contrasted results are reported from one study to another suggesting that drivers of classification errors are still under-explored. We evaluated the performances of single-year Formosat-2 time series to discriminate tree species in temperate forests in France and investigated how predictions vary statistically and spatially across multiple years. Our objective was to better estimate the impact of spatial autocorrelation in the validation data on measurement accuracy and to understand which drivers in the time series are responsible for classification errors. The experiments were based on ten Formosat-2 image time series irregularly acquired during the seasonal vegetation cycle from 2006 to 2014. Due to lot of clouds in the year 2006, an alternative 2006 time series using only cloud-free images has been added. Thirteen tree species were classified in each single-year dataset based on the SVM algorithm. The performances were assessed using a spatial leave-one-out cross validation (SLOO-CV) strategy, thereby guaranteeing full independence of the validation samples, and compared with standard non-spatial leave-one-out cross-validation (LOO-CV). The results show relatively close statistical performances from one year to the next despite the differences between the annual time series. Good agreements between years were observed in monospecific tree plantations of broadleaf species versus high disparity in other forests composed of different species. A strong positive bias in the accuracy assessment (up to 0.4 of Overall Accuracy) was also found when spatial dependence in the validation data was not removed. Using the SLOO-CV approach, the average OA values per year ranged from 0.48 for 2006 to 0.60 for 2013, which satisfactorily represents the spatial instability of species prediction between years.

In Portugal, cork oak (Quercus suber L.) stands cover 737 Mha, being the most predominant species of the montado agroforestry system, contributing for the economic, social and environmental development of the country. Cork oak decline is a known problem since the late years of the 19th century that has recently worsen. The causes of oak decline seem to be a result of slow and cumulative processes, although the role of each environmental factor is not yet established. The availability of Sentinel-2 high spatial and temporal resolution dense time series enables gradual processes monitoring. These processes can be monitored using spectral vegetation indices (VI) once their temporal dynamics are expected to be related with green biomass and photosynthetic efficiency. The Normalized Difference Vegetation Index (NDVI) is sensitive to structural canopy changes, however it tends to saturate at moderate-to-dense canopies. Modified VI have been proposed to incorporate the reflectance in the red-edge spectral region, which is highly sensitive to chlorophyll content while largely unaffected by structural properties. In this research, in-situ data on the location and vitality status of cork oak trees are used to assess the correlation between chlorophyll indices (CI) and NDVI time series trends and cork oak vitality at the tree level. Preliminary results seem to be promising since differences between healthy and unhealthy (diseased/dead) trees were observed.

The complexity of Mixed-Integer Linear Programs (MILPs) increases with the number of nodes in energy system models. An increasing complexity constitutes a high computational load that can limit the scale of the energy system model. Especially in microgrid optimisation problems with multiple buildings and energy systems with a number of rival supply, distribution and storage technologies, methods are sought to reduce this complexity. In this paper, we present a new 2-Level Approach to MILP energy system models that determine the system design through a combination of continuous and discrete decisions. On the first level, data reduction methods are used to determine the discrete design decisions in a simplified solution space. Those decisions are then fixed, and on the second level the full dataset is used to extract the exact scaling of the chosen technologies. The performance of the new 2-Level Approach is evaluated for a case study of an urban energy system with six buildings and an island system based on a high share of renewable energy technologies. The results of the studies show a high accuracy with respect to the total annual costs, chosen system structure, installed capacities and peak load with the 2-Level Approach compared to the results of a single level optimization. The computational load is thereby reduced by more than one order of magnitude, while a significantly higher accuracy is reached in comparison to the common time series aggregation approach.

The goal of VAR or BVAR is the characterization of the dynamics and endogenous relationships among time series. Also the VAR models are known for their applications to forecasting and policy analysis. This paper compare the performance of VAR and Sims-Zha Bayesian VAR models when the multiple time series are jointly influenced by different levels of collinearity and autocorrelation in the short term (T=16, 32, 64 and 128). Five levels (-0.9,-0.5, 0,+0.5,+0.9) of collinearity and autocorrelation were considered and the results from the simulation study revealed that VAR(2) model dominated for no and moderate levels of autocorrelation (-0.5, 0, +0.5) irrespective of the collinearity level except in few cases when T=16. While the BVAR models dominated for high autocorrelation levels (-0.9 and +0.9) irrespective of the collinearity level except in few cases when T=128. The performance of the models varies at different levels of the collinearity and autocorrelated error, and also varies with the short term periods. Furthermore, the values of the RMSE and MAE criteria decrease as a result of increase in the time series length. In conclusion, the performance of the forecasting models depend on the time series data structure and the time series length. It is therefore recommended that the data structure and series length should be considered in using an appropriate model for forecasting.

This paper has three objectives, the first is to present a detailed overview in the form of a tutorial for the developments of several key quantile time series modelling approaches. The second objective is to develop a general framework to represent such quantile models in a unifying manner in order to easily develop extensions and connections between existing models that can then be developed to further extend these models in practice. In this regard, the core theme of the paper is to provide perspectives to a general audience of core components that go into construction of a quantile time series model and then to explore each of these core components in detail. The paper is not addressing the concerns of estimation of these models, as there is existing literature on these aspects in many settings, we provide references to relevant works on these aspects in several classes of model. Instead, the focus is rather to provide a unified framework to construct such models for practitioners, therefore the focus is instead on the properties of the models and links between such models from a constructive perspective. The third objective is to compare and discuss the application of the different quantile time series models on several sets of interesting demographic and mortality based time series data sets of relevance to life insurance analysis. The exploration included detailed mortality, fertility, births and morbidity data in several countries with more detailed analysis of regional data in England, Wales and Scotland.

Anomaly detection plays a significant role in helping gas turbines run reliably and economically. Considering collective anomalous data and both sensitivity and robustness of the anomaly detection model, a sequential symbolic anomaly detection method is proposed and applied to the gas turbine fuel system. A structural Finite State Machine is to evaluate posterior probabilities of observing symbolic sequences and most probable state sequences they may locate. Hence an estimating based model and a decoding based model are used to identify anomalies in two different ways. Experimental results indicates that these two models have both ideal performance overall, and estimating based model has a strong ability in robustness, while decoding based model has a strong ability in accuracy, particularly in a certain range of length of sequence. Therefore, the proposed method can well facilitate existing symbolic dynamic analysis based anomaly detection methods especially in gas turbine domain.

: In this article, the authors propose an algorithmic approach to building a model of the dynamics of economic and, in particular, innovation processes. The approach under consideration is based on a complex algorithm that includes (1) decomposition of the time series into components using singular spectrum analysis; (2) recognition of the optimal component model based on fuzzy rules, and (3) creation of statistical models of individual components with their combination. It is shown that this approach corresponds to the high uncertainty characteristic of the tasks of the dynamics of innovation processes. The proposed algorithm makes it possible to create effective models that can be used both for analysis and for predicting the future states of the processes under study. The advantage of this algorithm is the possibility to expand the base of rules and components used for modeling. This is an important condition for improving the algorithm and its applicability for solving a wide range of problems.

Satellite Image Time Series (SITS) is a data set which includes satellite images across several years with a high acquisition rate. Radiometric normalization is a fundamental and important preprocessing method for remote sensing applications using SITS due to the radiometric distortion caused by noise between images. Normalizing the subject image based on the reference image is a general strategy when using traditional radiometric normalization methods to normalize multi-temporal imagery (usually two scenes or three scenes in different time phases). However, these methods are not suitable for calibrating SITS because they cannot minimize the radiometric distortion between any pair of images in SITS. The existing relative radiometric normalization methods for SITS are based on linear assumptions which cannot effectively reduce nonlinear radiometric distortion caused by continuously changing noise in SITS. To overcome this problem and obtain a more accurate SITS, this study proposed a Nonlinear Radiometric Normalization Model (named NMAG) for SITS based on Artificial Neural Networks (ANN) and Greedy Algroithm (GA) . In this method, GA was used to determine the correction order of SITS and calculate the error between the image to be corrected and normalized images, which avoided the selection of a single reference image. ANN was used to obtain the optimal solution of error function, which minimized the radiometric distortion between different images in SITS. SITS composed of 21 Landsat-8 images in Tianjin City from October 2017 to January 2019 were selected to test the method. We compared NMAG with other two contrasting method (refered as CM1 and CM2), and found that the average of root mean square error $(\mu_{RMSE})$ of NMAG $(497.22)$ was significantly smaller than those of CM1 $(641.39)$ and CM2 $(543.47)$, and the accuracy of normalized SITS obtained by using NMAG has increased by 22.4\% and 8.5\% toward CM1 and CM2, respectively. These experimental results confirmed the effectiveness of NMAG in the reduction of radiometric distortion caused by continuously changing noise between images in SITS.

(1) Introduction: The unprecedented coronavirus disease 2019 (COVID-19) epidemic has caused millions of infections worldwide and represents a significant challenge facing modern health care systems. This study was conducted to investigate the impact of lockdown measures on regional neonate services, which might be used to predict the long-term effects of medical behavior. (2) Methods: Using hospital information system (HIS) Statistics data from 4 January 2019 to 27 August 2020, an interrupted time-series analysis was employed to compare changes and trends in hospital admissions and disease spectrum before and after the lockdown interventions. Furthermore, this study was designed to evaluate whether the pandemic influenced newborns' healthcare behavior. (3) Results: Overall, 13,540 infants were admitted to the NICU during the pre-COVID period (n = 12082) and COVID-impacted period (n =4558). The patients' age at admission were younger than that of the pre-COVID-19 period (median age 5d vs. 6d after birth, p＜0.001). The overall number of neonate visits consistently decreased from the first days of the lockdown measures (24 January 2020). The disease spectrum for respiratory system, infectious diseases, and gastrointestinal disease indicated no declined immediately after intervention(p =0.079, p =0.113, p =0.060, respectively). There was an immediate decline in the volume of Jaundice-related conditions (p<0.001) after lockdown measures. The percentage of patients who suffer from respiratory system and infectious diseases has decreased (p =0.005 and p =0.002). However, a relatively high percentage of patients admitted to the neonatal intensive care unit (NICU) presented with Jaundice-related conditions (p <0.001). (4) Conclusions: In summary, the COVID-19 pandemic profoundly impacted the regional neonate services. However, it is still unclear what might be the effect of long-term effects from pandemic.

Natural circulation loop (NCL) systems are buoyancy-driven heat exchangers that are used in various industrial applications. The concept of passive heat exchange in NCL systems is attractive, because there is no need for an externally driven equipment (e.g., a pump) to maintain the fluid circulation. However, relying on buoyancy as the sole driving force may lead to several potential difficulties, one of which is generation of (possibly) time-varying nonlinearities in the dynamical system, where a difference in the time scales of heat transfer and fluid flow causes the flow to change from a steady-state regime to either an oscillatory regime or a flow-reversal regime, both of which are undesirable. In this paper, tools of symbolic time series analysis (e.g., probabilistic finite state automata (PFSA)) are proposed to identify selected regimes of operation in NCL systems, where the underlying principle is built upon the concept of pattern classification from measurements of fluid-flow dynamics. The proposed method is shown to be capable of identifying the current regime of operation from the initial time response under a given set of operational parameters. The efficacy of regime classification is demonstrated by testing on two data sets, generated from numerical simulation of a MATLAB SimuLink model that has previously been validated with experimental data. The results of the proposed PFSA-based classification are compared with those of a hidden Markov model (HMM) that serves as the baseline.

An application of deep convolutional neural network and recurrence plot for financial market movement prediction is presented. Though it is challenging and subjective to interpret its information, the pattern formed by a recurrence plot provide a useful insight into the dy- namical system. We used a recurrence plot of seven financial time series to train a deep neural network for financial market movement predic- tion. Our approach is tested on our dataset and achieved an average of 53.25% classification accuracy. The result suggests that a well trained deep convolutional neural network can learn a recurrence plot and pre- dict a financial market direction.

The crucial stages in the geochemical evolution of the Earth’s crust, ocean, and atmosphere could be explained by the assumed low-energy nuclear reactions (LENR) that are triggered by seismic activity. LENR result in the fission of medium-weight elements accompanied by neutron emissions, involving Fe and Ni as starting elements, and C, N, O as resultants. Geochemical data and experimental evidences support the LENR hypothesis. A spectral analysis of the period 1955-2013 shows common cycles between interannual changes in atmospheric CO2 growth rate and global seismic-moment release, whereas the trending behavior of the atmospheric CO2 was in response to the anthropogenic emissions. Assuming a correlation between such seismic and atmospheric fluctuations, the latter could be explained by cycles of worldwide seismicity, which would trigger massively LENR in the Earth’s Crust. In this framework, LENR from active faults could be considered as a relevant cause of carbon formation and degassing of freshly-formed CO2 during seismic activity. However, further studies are necessary to validate the present hypothesis which, at the present time, mainly aims to stimulate debate on the models which regulates atmospheric CO2.

While the COVID-19 outbreak was reported to first originate from Wuhan, China, it has been declared as a Public Health Emergency of International Concern (PHEIC) on 30 January 2020 by WHO, and it has spread to over 180 countries by the time of this paper was being composed. As the disease spreads around the globe, it has evolved into a worldwide pandemic, endangering the state of global public health and becoming a serious threat to the global community. To combat and prevent the spread of the disease, all individuals should be well-informed of the rapidly changing state of COVID-19. In the endeavor of accomplishing this objective, a COVID-19 real-time analytical tracker has been built to provide the latest status of the disease and relevant analytical insights. The real-time tracker is designed to cater to the general audience without advanced statistical aptitude. It aims to communicate insights through various straightforward and concise data visualizations that are supported by sound statistical foundations and reliable data sources. This paper aims to discuss the major methodologies which are utilized to generate the insights displayed on the real-time tracker, which include real-time data retrieval, normalization techniques, ARIMA time-series forecasting, and logistic regression models. In addition to introducing the details and motivations of the utilized methodologies, the paper additionally features some key discoveries that have been derived in regard to COVID-19 using the methodologies.

A major task in Functional Time Series Analysis is measuring the dependence within and between processes, for which lagged covariance and cross-covariance operators have proven to be a practical tool in well-established spaces. This article deduces estimators and asymptotic upper bounds of the estimation errors for lagged covariance and cross-covariance operators of processes in Cartesian products of abstract Hilbert spaces for fixed and increasing lag and Cartesian powers. We allow the processes to be non-centered, and to have values in different spaces when investigating the dependence between processes. Also, we discuss features of estimators for the principle components of our covariance operators.

This study addresses the causal links between external factors and the main hydro-climatic variables. There is a gap in the literature on the description of a complete chain in addressing the structures of direct causal links of solar activity on the terrestrial variables. This is why, the present study uses the extensive facilities of the application of information theory in view of recent advances in different fields. Also, by other methods (e.g. neural networks) first are tested the existence non-linear links of solar-terrestrial influences on hydro-climate system. The results are promising related to the solar impact on terrestrial phenomena which is discriminant in space-time domain. The implications prove robust for determining the causal measure of climate variables under direct solar impact which makes it easier to consider solar activity in climate models, by appropriate parametrizations. This study found that hydro-climatic variables are sensitive to solar impact only for certain frequencies (periods) and these have a coherence with the Solar-Flux only for some lags of the Solar-Flux (in advance).

Since December 2019, the coronavirus disease (COVID-19) has rapidly spread worldwide. The Mexican government has implemented public safety measures to minimize the spread of the virus. In this paper, the authors use statistical models in two stages to estimate the total number of coronavirus (COVID-19) cases per day at the state and national level in Mexico. Two types of models are proposed: first, a polynomial model of the growth for the first part of the outbreak until the inflection point of the pandemic curve and then a second nonlinear growth model is used to estimate the middle and the end of the outbreak. Model selection will be performed using Vuong’s test. The proposed models show overall fit similar to predictive models (e.g. time series, and machine learning); however, the interpretation of parameters is less complex for decision-makers and the residuals follow the expected distribution when fitting the models without autocorrelation being an issue.

This article explores the required amount of time series points from a high-speed traffic network to accurately estimate the Hurst exponent. The methodology consists in designing an experiment using estimators that are applied to time series, followed by addressing the minimum amount of points required to obtain accurate estimates of the Hurst exponent in real-time. The methodology addresses the exhaustive analysis of the Hurst exponent considering bias behavior, standard deviation, mean square error, and convergence using fractional Gaussian noise signals with stationary increases. Our results show that the Whittle estimator successfully estimates the Hurst exponent in series with few points. Based on the results obtained, a minimum length for the time series is empirically proposed. Finally, to validate the results, the methodology is applied to real traffic captures in a high-speed network based on the IEEE 802.3ab standard.

Accurate inundation maps for flooded wetlands and rivers are a critical resource for their management and conservation. In this paper we automate a method (thresholding of the short-wave infrared band) for classifying inundation, using Landsat imagery and Google Earth Engine. We demonstrate the method in the Okavango Delta, northern Botswana, a complex case study due to the spectral overlap between inundated areas covered with aquatic vegetation and dryland vegetation classes on satellite imagery. Inundation classifications in the Okavango Delta have predominately been implemented on broad spatial resolution images. We present the longest time series to date (1990-2019) of inundation maps at high spatial resolution (30m) for the Okavango Delta. We validated the maps using image-based and in situ data accuracy assessments, with accuracy ranging from 91.5 - 98.1%. Use of Landsat imagery resulted in consistently lower estimates of inundation extent than previous studies, likely due to the increased number of mixed pixels that occur when using broad spatial resolution imagery, which can lead to overestimations of the size of inundated areas. We provide the inundation maps and Google Earth Engine code for public use.

The variability of meteorological parameters such as temperature and precipitation, and climatic conditions such as intense droughts, are known to impact vegetation health over southern Africa. Thus, understanding large-scale ocean–atmospheric phenomena like the El Niño/Southern Oscillation (ENSO) and Indian Ocean Dipole/Dipole Mode Index (DMI) is important as these factors drive the variability of temperature and precipitation. In this study, 16 years (2002–2017) of Moderate Resolution Imaging Spectroradiometer (MODIS) Terra/Aqua 16-day normalized difference vegetation index (NDVI), extracted and processed using JavaScript code editor in the Google Earth Engine (GEE) platform in order to analyze the response pattern of the oldest proclaimed nature reserve in Africa, the Hluhluwe-iMfolozi Park (HiP), during the study period. The MODIS-enhanced vegetation index and burned area index were also analyzed for this period. The area-averaged Modern Retrospective Analysis for Research Application (MERRA) model maximum temperature and precipitation were also extracted using the JavaScript code editor in the GEE platform. This procedure demonstrated a strong reversal of both the NDVI and Enhanced Vegetation Index (EVI), leading to signs of a sudden increase of burned areas (strong BAI) during the strongest El Niño period. Both the Theilsen method and the Mann–Kendall test showed no significant greening or browning trends over the whole time series, although the annual Mann–Kendall test, in 2003 and 2014–2015, indicated significant browning trends due to the most recent strongest El Niño. Moreover, a multi-linear regression model seems to indicate a significant influence of both ENSO activity and precipitation. Our results indicate that the recent 2014–2016 drought altered the vegetation condition in the HiP. We conclude that it is vital to exploit freely available GEE resources to develop drought monitoring vegetation systems, and to integrate climate information for analyzing its influence on protected areas, especially in data-poor counties.

This paper is aimed at studying the temporal evolution of the surface displacements occurred over the past few years in the ocean-reclaimed platforms of the Shanghai megacity (China), which are mainly ascribable to consolidation processes of large dredger fills and alluvial deposits. With respect to previous analyses carried out over the same area, this work provides a joint multi-platform differential interferometry synthetic aperture radar (DInSAR) analysis, based on the application of the advanced Small BAseline Subset (SBAS) algorithm. This led us to retrieve long-term deformation time-series that are helpful for a better understanding of the on-going deformation phenomena. To this aim, we have exploited two sequences of SAR data collected by the ASAR/ENVISAT and by the COSMO-SkyMed (CSK) sensors, respectively, spanning the whole time period from 2007 to 2016. Unfortunately, the large time gap (of about three years) existing between the available ASAR/ENVISAT and CSK datasets gave rise to additional difficulties for their combination. Nevertheless, this problem has been faced by benefiting from the knowledge of a time-dependent model describing the temporal evolution of the expected deformations affecting the Shanghai ocean-reclaimed platforms.

The electric power distribution utility is responsible for providing energy to consumers in a continuous and stable way, failures in the electrical power system reduce the reliability indexes of the grid, directly harming its performance. For this reason, there is a need for failure prediction to reestablish power in the shortest possible time. Considering an evaluation of the number of failures over time, this paper proposes to perform a failure prediction during the first year of the pandemic in Brazil (2020) to verify the feasibility of using time series forecasting models for fault prediction. The Long Short-Term Memory (LSTM) model will be evaluated to obtain a forecast result that can be used by the electric power utility to organize the maintenance teams. The Wavelet transform shows to be promising in improving the predictive ability of the LSTM, making the Wavelet LSTM model suitable for the study at hand. The results show that the proposed approach has better results regarding the evaluation of the error in prediction and has robustness when a statistical analysis is performed.

Snow is a vital environmental parameter and dynamically responsive to climate change, particularly in mountainous regions. Snow cover can be monitored at variable spatial scales using Earth Observation (EO) data. Long-lasting remote sensing missions enable the generation of multi-decadal time series and thus the detection of long-term trends. However, there have been few attempts to use these to model future snow cover dynamics. In this study, we therefore explore the potential of such time series to forecast the Snow Line Elevation (SLE) in the European Alps. We generate monthly SLE time series from the entire Landsat archive (1985-2021) in 43 Alpine catchments. Positive long-term SLE change rates are detected, with the highest rates (5-8 m/y) in the Western and Central Alps. We utilize this SLE dataset to implement and evaluate seven uni-variate time series modeling and forecasting approaches. The best results were achieved by Random Forests, with a Nash-Sutcliffe efficiency (NSE) of 0.79 and a Mean Absolut Error (MAE) of 258 m, Telescope (0.76, 268 m), and seasonal ARIMA (0.75, 270 m). Since the model performance varies strongly with the input data, we developed a Combined forecast based on the best performing methods in each catchment. This approach was then used to forecast the SLE for the years 2022-2029. In the majority of the catchments the shift of the forecast median SLE level retained the sign of the long-term trend. In cases where a deviating SLE dynamic is forecast a discussion based on the unique properties of the catchment and past SLE dynamics is required. In the future, we expect major improvements in our SLE forecasting efforts by including external predictor variables in a multi-variate modeling approach.

We developed a novel approach in the field of spatiotemporal modelling, based on the spatialisation of time: the Timescape algorithm. It is especially aimed at sparsely distributed datasets in ecological research, whose spatial and temporal variability is strongly entangled. The algorithm is based on the definition of a spatiotemporal distance that incorporates a causality constraint and that is capable of accommodating the seasonal behaviour of the modelled variable as well. The actual modelling is conducted exploiting any established spatial interpolation technique, substituting the ordinary spatial distance with our Timescape distance, thus sorting, from the same input set of observations, those causally related to each estimated value at a given site and time. The notion of causality is expressed topologically and it has to be tuned for each particular case. The Timescape algorithm originates from the field of stable isotopes spatial modelling (isoscapes), but in principle it can be used to model any real scalar random field distribution.

An earthquake of Mw6.4 hit the coastal zone of Albania on 26 November 2019, at 02:54:11 UTC. It was intensively felt at about 34km far away, in Tirana City, where damages and lives lost occurred. To emphasize a geomagnetic signature before the onset of this earthquake, the data collected on the interval 15 October–30 November 2019, at the Panagjurishte (PAG)-Bulgaria and Surlari (SUA)-Romania observatories are analyzed by using both the polarization parameter (BPOL)-time invariant in non-seismic conditions, becoming unstable before this seismic event, and the strain effect for geomagnetic signal identification. Consequently, BPOL time series and its standard deviations are performed for the both sites using ULF-FFT band-pass filtering. A statistical analysis, based on a standardized random variable equation, was applied to emphasize on the BPOL*(PAG) and ABS BPOL*(PAG) time series the anomalous signal’s singularity and, to differentiate the transient local anomalies due to Mw6.4earthquake, from the internal and external parts of the geomagnetic field, taken PAG observatory as reference. Finally, the ABS BPOL*(PAG-SUA) time series are obtained on the interval 1-30 November, 2019, where a geomagnetic signature greater than 2.0, was detected on 23 November and the lead time was 3 days before the onset of Mw6.4earthquake.

In the steel industries, workers are exposed to heat and ambient thermal stresses on a daily basis, leading to discomfort and limited performance. In this study, the main purpose is to investigate the effect of climate heat stress on the rate of accidents in the workplace for workers for 5 consecutive years. The data of this study were received without any sampling through the HSE Center for Steel Industry and meteorological data from 2015 to 2019 from Isfahan Meteorological station. The daily number of casualties among workers in the steel industry during 2015-2019 by adjusting seasonal patterns, months, effects of the day of the week and other meteorological factors on the average daily temperature using the studied model has a decreasing effect. Eviews software (version 8) was used to model and investigate the relationship between events and meteorological variables. The mean temperature was at least 40.2-2 and at most 70.34 ° C, respectively. In the time-series study for the main model, the number of accidents shows a direct relationship with the average temperature and wind speed. Climatic indices of humidity and rainfall have the least impact on accidents compared to temperature and wind speed. A strong correlation was shown between the increase in average ambient temperature and the rate of accidents over the past 5 years. Given the fundamental differences in studies of environmental exposure and wind speed over heat stress, further analysis in workers should be considered.

A timely inventory of agricultural areas and crop types is an essential requirement for ensuring global food security. Satellite remote sensing has proven to be an increasingly more reliable tool to identify crop types. With the Copernicus program and its Sentinel satellites, a growing source of satellite remote sensing data is publicly available at no charge. Here we use joint Sentinel-1 radar and Sentinel-2 optical imagery to create a crop map for Belgium. To ensure homogenous radar and optical input across the country, Sentinel-1 12-day backscatter composites were created after incidence angle normalization, and Sentinel-2 NDVI images were smoothed to yield dekadal cloud-free composites. An optimized random forest classifier predicted the 8 crop types with a maximum accuracy of 82% and a kappa coefficient of 0.77. We found that a combination of radar and optical imagery always outperformed a classification based on single-sensor inputs, and that classification performance increased throughout the season until July, when differences between crop types are largest. Furthermore we showed that the concept of classification confidence derived from the random forest classifier provided insight in the reliability of the predicted class for each pixel, clearly showing that parcel borders have a lower classification confidence. We concluded that the synergistic use of radar and optical data for crop classification led to richer information increasing classification accuracies compared to optical-only classification. Further work should focus on object-level classification and crop monitoring to exploit the rich potential of combined radar and optical observations.

Smart grids require flexible data driven forecasting methods. We propose clustering tools for bottom-up short-term load forecasting. We focus on individual consumption data analysis which plays a major role for energy management and electricity load forecasting. The two first sections are dedicated to the industrial context and a review of individual electrical data analysis. We are interested in hierarchical time-series for bottom-up forecasting. The idea is to disaggregate the signal in such a way that the sum of disaggregated forecasts improves the direct prediction. The 3-steps strategy defines numerous super-consumers by curve clustering, builds a hierarchy of partitions and selects the best one minimizing a forecast criterion. Using a nonparametric model to handle forecasting, and wavelets to define various notions of similarity between load curves, this disaggregation strategy applied to French individual consumers leads to a gain of 16\% in forecast accuracy. We then explore the upscaling capacity of this strategy facing massive data and implement proposals using R, the free software environment for statistical computing. The proposed solutions to make the algorithm scalable combines data storage, parallel computing and double clustering step to define the super-consumers.

Most of existing fuzzy forecasting models partition historical training time series into fuzzy time series and build fuzzy-trend logical relationship groups to generate forecasting rules. The determination process of intervals is complex and uncertainty. In this paper, we present a novel fuzzy forecasting model based on high-order fuzzy-fluctuation trends and the fuzzy-fluctuation logical relationships of the training time series. Firstly, we compare each data with the data of its previous day in historical training time series to generate a new fluctuation trend time series(FTTS). Then, fuzzify the FTTS into fuzzy-fluctuation time series(FFTS) according to the up, equal or down range and orientation of the fluctuations. Since the relationship between historical FFTS and the fluctuation trend of future is nonlinear, Particle Swarm Optimization (PSO) algorithm is employed to estimate the required parameters. Finally, use the acquired parameters to forecast the future fluctuations. In order to compare the performance of the proposed model with that of the other models, we apply the proposed method to forecast the Taiwan Stock Exchange Capitalization Weighted Stock Index (TAIEX) time series datasets. The experimental results and the comparison results show that the proposed method can be successfully applied in stock market forecasting or such kinds of time series. We also apply the proposed method to forecast Shanghai Stock Exchange Composite Index (SHSECI) to verify its effectiveness and universality.

Satellite Crop Detection technologies are focused on detection of different types of crops on the field in the early stage before harvesting. Crop detection is usually done on a time series of satellite data by classification of the desired fields. Currently, data obtained from Remote Sensing (RS) are used to solve tasks related to the identification of the type of agricultural crops, also modern technologies using AI methods are desired in the postprocessing part. In this challenge Sentinel-1 and Sentinel-2 time series data were used due to their periodic availability. Our focus was to develop methodology for classification of time series of Sentinel 2 and Sentinel 1 data and compare how accuracy of classification can be increased, but also how to guarantee availability of data. We analyse phenology of single crops and on the basis of this analysis we started to provide crop classification. Original crop classifications were made from Enhanced Vegetation Index (EVI) layers made from Sentinel-2 time-series data and then we added also . To increase accuracy we also integrate into the process parcel borders and provide classification of fields..

In 2019, Chinese government implemented the first round of National Centralized Drug Procurement (NCDP) pilot (so-called "4+7" policy) in mainland China, achieved a prominent price reduction of 52% on average for 25 bidding winning products. Under cross-price elasticity theory, the price behavior of pharmaceutical enterprises for policy-related drugs might change. This study used drug purchasing data from the Centralized Drug Procurement Survey in Shenzhen 2019, and applied single-group Interruption Time Series (ITS) design to examine the impact of "4+7" policy on the drug price index (DPI) of policy-related drugs. The ITS analysis showed that the DPI of winning (-0.183 per month, p&lt;0.0001) and non-winning (-0.034 per month, p=0.046) products significantly decreased after the implementation of "4+7" policy. No significant difference was found for the immediate change of DPI for alternative drugs (p=0.537), while a significant decrease in change trend was detected in the post-"4+7" policy period (-0.003 per month, p=0.014). The DPI of the overall policy-related drugs significantly decreased (-0.261 per month, p&lt;0.0001) after "4+7" policy. These findings indicate that the price behavior of pharmaceutical enterprises changed under NCDP policy, while the price linkage effect is still limited. It is necessary to further expand the scope of centralized purchased drugs and strengthen the monitoring of related drugs regarding price change and consumption structure.

The spatiotemporal variability of vegetation in the Middle East was investigated for the period 2001–2019 using the Moderate Resolution Imaging Spectroradiometer (MODIS) 16-day/500 m composites of the Normalized Difference Vegetation Index (NDVI; MOD13A1). The results reveal a strong increase in the NDVI coverage in the Middle East during the study period (R = 0.75, p-value = 0.05). In Egypt, the annual coverage exhibits the strongest positive trend (R = 0.99, p-value = 0.05). In Turkey, both the vegetation coverage and density increased from 2001 to 2019, which can be attributed to the construction of some of the biggest dams in the Middle East, such as the Atatürk and Ilisu dams. Significant increases in the annual coverage and maximum and average NDVI in Saudi Arabia are due to farming in the northern part of the country for which groundwater and desalinated seawater are used. The results of this study suggest that the main factors affecting the vegetation coverage in the Middle East are governmental policies. These policies can have a positive effect on the vegetation coverage in some countries such as Egypt, Saudi Arabia, Qatar, Kuwait, Iran, and Turkey.

A vast majority of the countries is under the economic and health crises due to the current epidemic of coronavirus disease 2019 (COVID-19). The present study analyzes the COVID-19 using time series, which is an essential gizmo for knowing the enlargement of infection and its changing behavior, especially the trending model. We have considered an autoregressive model with a non-linear time trend component that approximately converted into the linear trend using the spline function. The spline function split the COVID-19 series into different piecewise segments between respective knots and fitted the linear time trend. First, we obtain the number of knots with its locations in the COVID-19 series and then the estimation of the best-fitted model parameters are determined under Bayesian setup. The results advocate that the proposed model/methodology is a useful procedure to convert the non-linear time trend into a linear pattern of newly coronavirus case for various countries in the pandemic situation of COVID-19.

Advances in DNA sequencing technology are facilitating genomic analyses of unprecedented scope and scale, widening the gap between our abilities to generate and fully exploit biological sequence data. Comparable analytical challenges are encountered in other data-intensive fields involving sequential data, such as signal processing, in which dimensionality reduction (i.e., compression) methods are routinely used to lessen the computational burden of analyses. In this work we explore the application of dimensionality reduction methods to numerically represent high-throughput sequence data for three important biological applications of virus sequence data: reference-based mapping, short sequence classification and de novo assembly. Despite using highly compressed sequence transformations to accelerate the processes, our sequence processing approach yielded comparable accuracy to existing approaches, and are ideally suited for sequences originating from highly diverse virus populations. We demonstrate the application of our methodology to both synthetic and real viral pathogen sequence data. Our results show that the use of highly compressed sequence approximations can provide accurate results and that useful analytical performance can be retained and even enhanced through appropriate dimensionality reduction of sequence data.

This article shortly introduces Gaussian processes as a new approach for modelling time series in the field of blazar physics. In the second part of the paper, recent results from an application of GP modelling to the multi-wavelength light curves of the blazar PKS 1502+106 are discussed.

Drinking water quality is a major concern, especially in African countries. This manuscript aims to analyze the chemical composition of Lioua’s groundwater in order to determine the geological processes influencing the chemical elements' composition and origin. Therefore, chemometrics techniques such as multivariate statistical analysis (MSA) and time series methods (TSM) are used. Indeed, MSA includes a component analysis (PCA) and a cluster analysis (CA), while au-tocorrelation analysis (AA) supplemented by simple spectral density analysis (SDA) is used for TMS. PCA displays three main factors explaining a total variance (TV) of 85.01 %. Factors 1, 2, and 3 are 68.72%, 11.96%, and 8.89 % of TV, respectively. In the CA, three groups were controlled by TDS and EC. G1 reveals a close association between SO42−, K+, Ca2+, and TDS; G2 reveals a close association between Na+, Cl−, Mg2+, and EC; G3 shows the dissociation of bicarbonates HCO3− and NO3− from other chemical elements. AA shows a linear interrelationship of EC, Mg2+, Na+, K+, Cl−, and SO42−. However, NO3− and HCO3− indicate uncorrelated characteristics with other parameters. For SDA, the correlograms of Mg2+, Na+, K+, Cl−, and SO42− have a similar trend with EC. None-theless, pH, Ca2+, HCO3− and NO3− exhibit multiple peaks related to the presence of several dis-tinct cyclic mechanisms. The methods enabled the authors to conclude that the geochemical processes influencing the chemical composition are: (i) dissolution of evaporated mineral depos-its, (ii) water-rock interaction, and (iii) evaporation process. In addition, Groundwater exhibits two bipolar characteristics, one recorded with negative and positive charges on pH and Ca+ and another recorded only with negative charges on HCO3− and NO3−. On the other hand, SO42-, K+, Ca2+, and TDS are the major predominant elements in the groundwater’s chemical composition. The major participation of salts and chlorides is in the electrical conductivity of water. The dominance of the lithological factor in the overall mineralization of the Plio-Quaternary surface aquifer waters. The origins of HCO3− and NO3− are different. Indeed, carbonated for HCO3- has a carbonate origin, whereas NO3– has an anthropogenic origin. The salinity was affected by Mg2+, SO42-, Cl-, Na+, K+, and EC. Ca2+, HCO3− and NO3− are resulted from human activity fertilizers, the carbonate facies outcrops, and domestic sewage.

For more than half a century, Manfred Deistler has been contributing to the construction of the rigorous theoretical foundations of the statistical analysis of time series and more general stochastic processes. Half a century of unremitting activity is not easily summarized in a few pages. In this short note, we chose to concentrate on a relatively little-known aspect of Manfred's contribution which nevertheless had quite an impact on the development of one of the most powerful tools of contemporary time series and econometrics: dynamic factor models.

The environmental issues we are currently facing require long-term prospective efforts for sustainable growth. Renewable energy sources seem to be one of the most practical and efficient alternatives in this regard. Understanding a nation's pattern of energy use and renewable energy production is crucial for developing strategic plans. No previous study has been performed to explore the dynamics of power consumption with the change in renewable energy production on a country-wide scale. In contrast, a number of deep learning algorithms demonstrated acceptable performance while handling sequential data in the era of data-driven predictions. In this study, we developed a scheme to investigate and predict total power consumption and renewable energy production time series for eleven years of data using a Recurrent Neural Network (RNN). The dynamics of the interaction between the total annual power consumption and renewable energy production are investigated through extensive Exploratory Data Analysis (EDA) and a feature engineering framework. The performance of the model is found satisfactory through the comparison of the predicted data with the observed data, visualization of the distribution of the errors and Root Mean Squared Error (RMSE) value of 0.084. Higher performance is achieved through the increase in the number of epochs and hyperparameter tuning. The proposed framework can be used and transferred to investigate the trend of renewable energy production and power consumption and predict the future scenarios for different communities. Incorporation of the cloud-based platform into the proposed pipeline may lead to real-time forecasting.

Feature selection is a crucial step to overcome the curse of dimensionality problem in data mining. This work proposes Recursive k-means Silhouette Elimination (RkSE) as a new unsupervised feature selection algorithm to reduce dimensionality in univariate and multivariate time-series datasets. Where k-means clustering is applied recursively to select the cluster representative features, following a unique application of silhouette measure for each cluster and a user-defined threshold as the feature selection or elimination criteria. The proposed method is evaluated on a hydraulic test rig, multi sensor readings in two different fashions: (1) Reduce the dimensionality in a multivariate classification problem using various classifiers of different functionalities. (2) Classification of univariate data in a sliding window scenario, where RkSE is used as a window compression method, to reduce the window dimensionality by selecting the best time points in a sliding window. Moreover, the results are validated using 10-fold cross validation technique. As well as, compared to the results when the classification is pulled directly with no feature selection applied. Additionally, a new taxonomy for k-means based feature selection methods is proposed. The experimental results and observations in the two comprehensive experiments demonstrated in this work reveal the capabilities and accuracy of the proposed method.

We present and benchmark FilterNet, a flexible deep learning architecture for time series classification tasks, such as activity recognition via multichannel sensor data. It adapts popular CNN and CNN-LSTM motifs which have excelled in activity recognition benchmarks, implementing them in a many-to-many architecture to markedly improve frame-by-frame accuracy, event segmentation accuracy, model size, and computational efficiency. We propose several model variants, evaluate them alongside other published models using the Opportunity benchmark dataset, demonstrate the effect of model ensembling and of altering key parameters, and quantify the quality of the models’ segmentation of discrete events. We also offer recommendations for use and suggest potential model extensions. FilterNet advances the state of the art in all measured accuracy and speed metrics on the benchmarked dataset, and it can be extensively customized for other applications.

Risk modelling along with multi-objective optimization problems have been at theepicenter of attention for supply chain managers. In this paper, we introduce a datasetfor risk modelling in sophisticated supply chain networks based on formal mathematical models. We have discussed the methodology and simulation tools used to synthesize the dataset. Additionally, the underlying mathematical models are discussed in granular details along with providing directions to conducting statistical analyses or neural machine learning models. The simulation is performed using MATLAB ™Simulink and the models are illustrated as well.

Today massive amounts of sequenced metagenomic and -transcriptomic data from different ecological niches and environmental locations are available. Scientific progress depends critically on methods that allow extracting useful information from the various types of sequence data. Here, we will first discuss types of information contained in the various flavours of biological sequence data, and how this information can be interpreted to increase our scientific knowledge and understanding. We argue that a mechanistic understanding is required to consistently interpret experimental observations, and that this understanding is greatly facilitated by the generation and analysis of dynamic mathematical models. We conclude that, in order to construct mathematical models and to test mechanistic hypotheses, time-series data is of critical importance. We review diverse techniques to analyse time-series data and discuss various approaches by which time-series of biological sequence data was successfully used to derive and test mechanistic hypotheses. Analysing the bottlenecks of current strategies in the extraction of knowledge and understanding from data, we conclude that combined experimental and theoretical efforts should be implemented as early as possible during the planning phase of individual experiments and scientific research projects.

Time recordings of impulse-type oscillation responses are short and highly transient. These characteristics may complicate the usage of classical spectral signal processing techniques for a) describing the dynamics and b) deriving discriminative features from the data. However, common model identification and validation techniques mostly rely on steady-state recordings, characteristic spectral properties and non-transient behavior. In this work, a recent method, which allows reconstructing differential equations from time series data, is extended for higher degrees of automation. With special focus on short and strongly damped oscillations, an optimization procedure is proposed that fine-tunes the reconstructed dynamical models with respect to model simplicity and error reduction. This framework is analyzed with particular focus on the amount of information available to the reconstruction, noise contamination and non-linearities contained in the time series input. Using the example of a mechanical oscillator, we illustrate how the optimized reconstruction method can be used to identify a suitable model and to extract features from uni-variate and multivariate time series recordings in an engineering-compliant environment. Moreover, the determined minimal models allow for identifying the qualitative nature of the underlying dynamical systems as well as testing for the degree and strength of non-linearity. The reconstructed differential equations would then be potentially available for classical numerical studies, such as bifurcation analysis. These results represent a physically interpretable enhancement of data-driven modeling approaches in structural dynamics.

Context: Anomaly detection in a data center is a challenging task, having to consider different services on various resources. Current literature shows the application of artificial intelligence techniques to either log files or monitoring data: the former created by services at run time, while the latter produced by specific sensors directly on the physical or virtual machine. Objectives: We propose a model that exploits information both in log files and monitoring data to identify patterns and detect anomalies over time. Methods: The key idea is to use on one side natural language processing solutions to detect problems at service level, extracting words that represent anomalies. Clustering and topic modeling techniques have been used to identify patterns and group them with respect to topics. On the other side time series anomaly detection technique has been applied to sensors data in order to combine problems found in the log files with problems stored in the monitoring data. Results: We have tested our approach on a real data center equipped with log files and monitoring data that can characterize the behaviour of physical and virtual resources in production. We have observed a correspondence between anomalies in log files and monitoring data, e.g. an increase in memory usage or machine load. The results are extremely promising. Conclusion: Our model requires to integrate site administrators' expertise in order to consider all critical scenario in the data center and understand results properly.

We investigate the power of time series analysis based on a variety of information-theoretic approaches from statistics (AIC, BIC) and machine learning (Minimum Message Length) - and we then compare their efficacy with traditional time series model and with hybrids involving deep learning. More specifically, we develop AIC, BIC and Minimum Message Length (MML) ARMA (autoregressive moving average) time series models - with this Bayesian information-theoretic MML ARMA modelling already being new work. We then study deep learning based algorithms in time series forecasting, using Long Short Term Memory (LSTM), and we then combine this with the ARMA modelling to produce a hybrid ARMA-LSTM prediction. Part of the purpose of the use of LSTM is to seek capture any hidden information in the residuals left from the traditional ARMA model. We show that MML not only outperforms earlier statistical approaches to ARMA modelling, but we further show that the hybrid MML ARMA-LSTM models outperform both ARMA models and LSTM models.

An accurate outbreak prediction of COVID-19 can successfully help to get insight into the spread and consequences of infectious diseases. Recently, machine learning (ML) based prediction models have been successfully employed for the prediction of the disease outbreak. The present study aimed to engage an artificial neural network-integrated by grey wolf optimizer for COVID-19 outbreak predictions by employing the Global dataset. Training and testing processes have been performed by time-series data related to January 22 to September 15, 2020 and validation has been performed by time-series data related to September 16 to October 15, 2020. Results have been evaluated by employing mean absolute percentage error (MAPE) and correlation coefficient (r) values. ANN-GWO provided a MAPE of 6.23, 13.15 and 11.4% for training, testing and validating phases, respectively. According to the results, the developed model could successfully cope with the prediction task.

Forest fires are part of a set of natural disasters that have always affected regions of the world typically characterized by a tropical climate with long periods of drought. However, due to climate change in recent years, other regions of our planet have also been affected by this phenomenon, never seen before. One of them is certainly the Italian peninsula, and especially the regions of southern Italy. For this reason, the scientific community, as well as remote sensing one, is highly concerned in developing reliable techniques to provide useful support to the competent authorities. In particular, three specific tasks have been carried out in this work: (i) fire risk prevention, (ii) active fire detection, and (iii) post-fire area assessment. To accomplish these analyses, the capability of a set of spectral indices, derived from spaceborne remote sensing (RS) data, is assessed to monitor the forest fires. The spectral indices are obtained from Sentinel-2 multispectral images of the European Space Agency (ESA), which are free of charge and openly accessible. Moreover, the twin Sentinel-2 sensors allow to overcome some restrictions on time delivery and observation repeat time. The performance of the proposed analyses were assessed experimentally to monitor the forest fires occurred in two specific study areas during the summer of 2017: the volcano Vesuvius, near Naples, and the Lattari mountains, near Sorrento (both in Campania, Italy).

IoT system becomes hot topic nowadays for smart home, IoT helps devices to communicate together without human intervention inside home, so it is offering many challenges. A new smart home IoT platform powered using electromagnetic energy harvesting is proposed in this paper. It contains a high gain transmitted antenna array and efficient circular polarized array rectenna system to harvest enough power from any direction to increase life time of the batteries used in IoT system. Optimized energy consumption, the software with adopting the Zigbee protocol of the sensor node and low power microcontroller are used to operate in lower power modes. The proposed system has an 84.6 days lifetime which is approximately 10 times the lifetime for similar system. On the other hand, the proposed power management circuit operated at 0.3 V DC to boost the voltage to ~3.7V from radio frequency energy harvesting and manage battery level to increase the batteries lifetime. A predictive indoor environment monitoring system is designed based on a novel hybrid system to provide a non-static plan, approve energy consumption and avoiding failure of sensor nodes in smart home.