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.

We examine a nonlinear initial value problem both singularly perturbed in a complex parameter and singular in complex time at the origin. The study undertaken in this paper is the continuation of a joined work with A. Lastra published in 2015. A change of balance between the leading and a critical subdominant term of the problem considered in our previous work is performed. It leads to a phenomenon of coalescing singularities to the origin in the Borel plane w.r.t time for a finite set of holomorphic solutions constructed as Fourier series in space on horizontal complex strips. In comparison to our former study, an enlargement of the Gevrey order of the asymptotic expansion for these solutions relatively to the complex parameter is induced.

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.

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.

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.

We evaluate various sums involving the powers of Fibonacci and Lucas numbers.

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 semi-nonparametric (SNP) modeling of the return distribution has been proved to be a flexible and accurate methodology for portfolio risk management that allows two-step estimation of the dynamic conditional correlation (DCC) matrix. For this SNP-DCC model, we propose a stepwise procedure to compute pairwise conditional correlations under bivariate marginal SNP distributions, overcoming the curse of dimensionality. The procedure is compared to the assumption of Dynamic Equicorrelation (DECO), which is a parsimonious model when correlations among the assets are not significantly different but requires joint estimation of the multivariate SNP model. The risk assessment of both methodologies is tested for a portfolio on cryptocurrencies by implementing backtesting techniques and for different risk measures: Value-at-Risk, Expected Shortfall and Median Shortfall. The results support our proposal showing that the SNP-DCC model has better performance for a smaller confidence level than the SNP-DECO model, although both models perform similarly for higher confidence levels.

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.

New q-series in the spirit of Jacobi have been found in a publication rst published in 1884 written in Russian and translated into English in 1928. This work was found by chance and appears to be almost totally unknown. From these entirely new q-series, fresh lattice sums have been discovered and are presented here.

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.

We study the problem of predicting time series data using the autoregressive integrated moving average (ARIMA) model in an online manner. Existing algorithms require model selection, which is time consuming and inapt for the setting of online learning. Using adaptive online learning techniques, we develop algorithms for fitting ARIMA models with fewest possible hyperparameters. We analyse the regret bound of the proposed algorithms and examine their performance using experiments on both synthetic and real world datasets

The present study analyzed light emitting diodes (LEDs) as an output load and used a Taylor series to describe the characteristic curve based on the exponential characteristic of voltage and current. A prototype circuit of a flyback LED driver system was established to verify whether the theory is consistent with actual results. This study focused on the exponential relationship of LED voltage and current. Conventional simulations usually used linear models to present LED loads. However, the linear model resulted in considerable error between simulation and actual characteristics. Therefore, this study employed a Taylor series to describe the nonlinear characteristic of an LED load. Through precise calculations with Mathcad computation software, the error was effectively reduced. Moreover, the process clarified the influence of temperature on LEDs, which benefited the characteristic analysis of the entire system. Finally, a realized circuit of 120-W flyback LED drivers was established for conducting theory verification, including theoretic analysis and evaluation of the system design process of the flyback converter. The circuit simulation software SIMPLIS was used to demonstrate the system model, which enabled quick understanding of the system framework established in this study. Regarding LEDs, a commercially available aluminum luminaire was used as the output load. The measured results of the actual circuit and the simulation results were remarkably consistent. For the same system at the same temperature, the error between the simulation and actual results was less than 3%, which proved the reliability of the Taylor series simulation.

A method for measuring the mechanical parameters of viscoelastic polymers by nanoindentation technology was proposed and verified. Through the mechanical response of load-displacement curves at different loading rates, then creep compliances and relaxation modulus were fitted. Polyimide thin film was employed in this research and experiments for five different loading rates were conducted. The fitting load-displacement loading curves obtained by the inversion method were identical to the experimental curves at five different loading rates，confirming the validity of the method. Moreover, with the loading rates increased，the fitting curves were more consistent commensurately with the nanoindentation experiment. DMA experiments were tested, and the generalized Kelvin/ Maxwell model were used for fitting experiment data. Results from DMA tests generally agree well with data from nanoindentation method, thereby verifying the feasibility of the method. The Prony series obtained by the two methods were used to simulate the creep experiments, which further verified the method.

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.

In this paper, we investigate an analytical solution of a class of nonlinear fractional integro-differential equation base on a generalized fractional power series expansion. The fractional derivatives are described in the conformable's type. The new approach is a modified form of the well-known Taylor series expansion. The illustrative examples are presented to demonstrate the accuracy and effectiveness of the proposed method

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.

According to the fact that the low convergence level on the complete elliptic integral of the first kind for the modulus which having values approach to one. In this paper we propose novelty of the complete elliptic integral consists of the new infinite series. We apply the scheme of iteration by substituting the common modulus with own modulus function into the new infinite series. We obtained so many new exact formulas of the complete elliptic integral derived from this method correspond to the number of iteration order. On the other hand, it has been also obtained a lot of new modulus functions rather than common used previously. The calculation results show that the number of significant figures of the new infinite series of the complete elliptic integral of the first kind is increased more and more. It means more fastly convergent would be obtained comparison values with the previous infinite series.

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.

We evaluate various binomial sums involving the powers of Fibonacci and Lucas numbers.

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.

In the coming years, the Anthropocene Working Group (AWG) will submit its proposal on the ‘Anthropocene’ as a new geological epoch to the International Commission on Stratigraphy (ICS) for approval. If approved, the proposal will be send to the International Union of Geological Sciences (IUGS) for ratification. If the proposal is approved and ratified, the ‘Anthropocene’ will be formalised and the Holocene Series/Epoch will be officially terminated. Currently, the ‘Anthropocene’ is a broadly used term and concept in a wide range of scientific and non-scientific situations and, for many, the official acceptance of this term is only a matter of time. However, the AWG proposal, in its present state, seems to not fully meet the ICS requirements for a new geological epoch. This paper asks what could happen if the current ‘Anthropocene’ proposal is not formalised by the ICS/IUGS. The possible stratigraphic alternatives are evaluated on the basis of the more recent literature and the personal opinions of distinguished AWG and ICS members. The eventual impact on environmental sciences and on non-scientific sectors, where the ‘Anthropocene’ seems already firmly rooted and de facto accepted as a new geological epoch, are also discussed.

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 report on the fabrication and electrical characteristics of zinc-oxide (ZnO) based metal-insulator-semiconductor (MIS) type Schottky barrier diodes (SBHs). ZnO thin layer on the p-type silicon substrate was fabricated by atomic layer deposition (ALD). The structure and surface properties of the thin film were characterized by X-ray diffraction (XRD), atomic force microscope (AFM) and secondary ion mass spectrometer (SIMS). The current-voltage (I-V) characteristics of Al/ALD-grown ZnO/p-Si diodes were measured under dark at room temperature. The electrical parameters such as ideality factor (n), series resistance (R_s) and barrier height (ϕ_b) of the diodes were analyzed using standard thermionic emission (TE) theory, Norde and Cheung method. The barrier height value obtained from I-V and Cheung method was found to be 0.73 eV and 0.76 eV, respectively. The interface state density (Dit) of the diodes was determined from the I-V characteristics. The nonideal behavior of measured parameters suggested the presence of interface states. The obtained results showed that the prepared diode can be used for NIR Schottky photodetector applications.

In this paper we consider type II bivariate generalized power series Poisson distribution as a compound Poisson distribution with bivariate generalized power series compounding distribution. We obtain some properties, p.m.f. and conditional distributions. In addition we also give a brief discussion about the multivariate extension of this case. Then we introduce type II bivariate generalized power series Poisson process and consider a bivariate risk model with type II bivariate generalized power series Poisson model as the counting process. For this model we derive distribution of the time to ruin and bounds for the probability of ruin. We obtain partial integro-differential equation for the ruin probabilities and express its bivariate transform through two univariate boundary transforms,where one of the initial capitals is fixed at zero.

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.

With a longer-term goal of addressing the comparative behavior of the aqueous halides F$^-$, Cl$^-$, Br$^-$, and I$^-$ on the basis of quasi-chemical theory (QCT), here we study structures and free energies of hydration clusters for those anions. We confirm that energetically optimal $(\mathrm{H_2O})_n\mathrm{X}$ clusters, with X = Cl$^-$, Br$^-$, and I$^-$, exhibit \emph{surface} hydration structures. Computed free energies based on optimized surface hydration structures utilizing a harmonic approximation, typically (but not always) disagree with experimental free energies. To remedy the harmonic approximation, we utilize single-point electronic structure calculations on cluster geometries sampled from an AIMD (\emph{ab initio} molecular dynamics) simulation stream. This \emph{rough-landscape} procedure is broadly satisfactory and suggests unfavorable ligand crowding as the physical effect corrected. Nevertheless, this procedure can break down when $n \gtrsim 4$, with the characteristic discrepancy resulting from a relaxed definition of clustering in the identification of $(\mathrm{H_2O})_n\mathrm{X}$ clusters, including ramified structures natural in \emph{physical cluster theories.} With ramified structures, the central equation for the present rough-landscape approach can acquire some inconsistency. Extension of these physical cluster theories in the direction of QCT should remedy that issue, and should be the next step in this research direction.

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.

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.

The use of haptic devices in rehabilitation of impaired limbs has become rather popular, given the proven effectiveness in promoting recovery. In a standard framework, such devices are used in rehabilitation centers, where patients interact with virtual tasks, presented on a screen. To track their sessions, kinematic/dynamic parameters or performance scores are recorded. However, as Internet access is now available at almost every home, and in order to reduce the hospitalization time of the patient, the idea of doing rehabilitation at home is gaining wide consent. Medical care programs can be synchronized with the home rehabilitation device; patient data can be sent to the central server that could redirect to the therapist laptop (tele-healthcare). The controversial issue is that the recorded data do not actually represent the clinical conditions of the patients according to the medical assessment scales, forcing them to frequently undergo clinical tests at the hospital. To respond to this demand, we propose the use of a bilateral master/slave haptic system that could allow the clinician, who interacts with the master, to assess remotely and in real time the clinical conditions of the patient that uses the home rehabilitation device as the slave. In this paper, we describe a proof of concept to highlight the main issues of such an application, limited to one degree of freedom, and to the measure of the stiffness and range of motion of the hand.

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.

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

We examine a family of nonlinear q-difference-differential Cauchy problems obtained as a coupling of linear Cauchy problems containing dilation q-difference operators, recently investigated by the author, and quasi-linear Kowalevski type problems that involve contraction q-difference operators. We build up local holomorphic solutions to these problems. Two aspects of these solutions are explored. One facet deals with asymptotic expansions in the complex time variable for which a mixed type Gevrey and q-Gevrey structure is exhibited. The other feature concerns the problem of confluence of these solutions as q tends to 1.

We evaluate some binomial cubic Fibonacci sums.

A signal analysis of the genoma sequenced of coronavirus variants: B.1.1.7, B.1.135, B.1.429-B.1.427, B.1.525 and P1 is presented. We deal with a certain type of finite alternating sum series having independently distributed terms associated with binary (0,1) indicators for the nucleotide bases A,C,G,T. This method provides additional information to conventional Similarity comparisons and Power Spectrum approaches. It leads to uncover distinctive patterns regarding the intrinsic data organization of genomic sequences according to its progression along the nucleotide bases position. Hence, the method could be useful for survelliance of genoma variants.

A linear Cauchy problem with polynomial coefficients wich combines q-difference operators for q>1 and differential operators of irregular type is examined. A finite set of sectorial holomorphic solutions w.r.t the complex time is constructed by means of classical Laplace transforms. These functions share a common asymptotic expansion in the time variable which turns out to carry a double layers structure which couples q-Gevrey and Gevrey bounds. In the last part of the work, the problem of confluence of these solutions as q tends to 1 is investigated.

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

Care coordination (CC) for children with special healthcare needs and medical complexity (CSHCN-CMC) is challenging, and medication management is especially difficult for providers, parents/caregivers, and patients alike. While numerous strategies for CC have been suggested and implemented, barriers to medication optimization remain. The report describes the creation of a pediatric clinical pharmacotherapy practice, related standard operating procedures to assure consistent application of screening tools and care provision through comprehensive medication management (CMM), and establishment of a collaborative practice agreement (CPA) to guide drug therapy delegation, monitoring, and modification. The methodology of a prospective case series is also presented to highlight drug therapy problems and their resolution in CSHCN-CMC. Future opportunities to expand the practice for engagement in population health management as well as prior authorization activities on behalf of physicians will be discussed.

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.

We derive new infinite series involving Fibonacci numbers and Riemann zeta numbers. The calculations are facilitated by evaluating linear combinations of polygamma functions of the same order at certain arguments.

Through the half-unit circle area computation using the integration of the corresponding curve power series representation, we deduce a slow converging positive infinite series to $\pi$. However, by studying the remainder of that series we establish sufficiently close inequalities with equivalent lower and upper bound terms allowing us to estimate, more precisely, how the series approaches $\pi$. We use the obtained inequalities to compute up to four-digit denominator, what are in this sense, the best rational numbers that can replace $\pi$. It turns out that the well-known convergents of the continued fraction of $\pi$, $22/7$ and $355/113$ called, respectively, Yuel\"{u} and Mil\"{u} in China are the only ones found. Thus we apply a similar process to find rational estimations to $\pi+e$ where $e$ is taken as the power series of the exponential function evaluated at $1$. For rational numbers with denominators less than $2000$, the convergent $920/157$ of the continued fraction of $\pi+e$ turns out to be the only rational number of this type.

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.

The new mineral species oyonite, ideally Ag3Mn2Pb4Sb7As4S24, has been discovered in the Uchucchacua polymetallic deposit, Oyon district, Catajambo, Lima Department, Peru, as very rare black metallic subhedral to anhedral crystals, up to 100 μm in length, associated with orpiment, tennantite/tetrahedrite, menchettiite, and other unnamed minerals of the system Pb-Ag-Sb-Mn-As-S, in calcite matrix. Its Vickers hardness (VHN100) is 137 kg/mm2 (range 132–147). In reflected light, oyonite is weakly to moderately bireflectant and weakly pleochroic from dark grey to a dark green. Internal reflections are absent. Reflectance values for the four COM wavelengths (Rmin, Rmax (%) (λ in nm)) are: 33.9, 40.2 (471.1); 32.5, 38.9 (548.3), 31.6, 38.0 (586.6); and 29.8, 36.5 (652.3). Electron microprobe analysis gave (in wt %, average of 5 spot analyses): Cu 0.76 (2), Ag 8.39 (10), Mn 3.02 (7), Pb 24.70 (25), As 9.54 (12), Sb 28.87 (21), S 24.30 (18), total 99.58 (23). On the basis of 20 cations per formula unit, the chemical formula of oyonite is Cu0.38Ag2.48Mn1.75Pb3.79Sb7.55As4.05S24.12. The main diffraction lines are (d in Å, hkl and relative intensity): 3.34 (-312; 40), 3.29 (-520; 100), 2.920 (-132; 40), 2.821 (-232; 70), 2.045 (004; 50). The crystal structure study revealed oyonite to be monoclinic, space group P21/n, with unit-cell parameters a = 19.1806 (18), b = 12.7755 (14), c = 8.1789 (10) Å, β = 90.471 (11)°, V = 2004.1 (4) Å3, Z = 2. The crystal structure was refined to a final R1 = 0.032 for 6272 independent reflections. Oyonite belongs to the Sb-rich members of the andorite homeotypic sub-series within the lillianite homologous series. The name oyonite is after the Oyon district, Lima Department, Peru, the district where the type locality (Uchucchacua mine) is located.

For a given set of plasma parameters, along a single series (Lyman, Balmer etc) the lines with higher principal quantum number(n) lines get progressively wider, closer to each other, and start merging for a certain critical n. In the present work four different codes(with further options) are used to calculate the entire Balmer series for moderate and high electron densities. Particular attention is paid to the relevant physics, such as the cutoff criteria, strong and penetrating electron collisions.

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.

Lumbar support exoskeletons with active and passive actuators are currently the cutting-edge technology for preventing back injuries in workers while lifting heavy objects. However, many challenges still exist in both types of exoskeletons, including rigid actuators, risks of human-robot interaction, high battery consumption, bulky design, and limited assistance. In this paper, the design of a compact, lightweight energy storage device combined with rotary series elastic (ES-RSEA) is proposed for use in a lumbar support exoskeleton to increase the level of assistance and exploit the human bioenergy during the two stages of the lifting task. ES takes the responsibility to store and release passive mechanical energy while RSEA provides excellent compliance and prevents injury from the human body's undesired movement. The experimental tests on the spiral spring showed excellent linear characteristics (above 99%) with an actual spring stiffness of 9.96 Nm/rad. The results demonstrate that ES-RSEA can provide maximum torque assistance in the ascent phase with 66.6 Nm while generating nearly 21 Nm of spring torque during descent without turning on the DC motor. Ultimately, the proposed design can maximize the energy storage of human energy, exploit the biomechanics of lifting tasks, and reduce the burden on human effort to perform lifting tasks.

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.

The design flood concept (DF) provides for an essential tool in designing the hydraulic works, in defining the reservoir operation programs and for a reliable flood hazard identification. Under a simplified approach, the maximum discharge and the floods volume are statistically processed to reasonably define the DF. Yet, the integral hydrograph provides additional key temporal and quantitative details of important significance for flood management and particularly for the res-ervoirs operation and associated risks of failures. The procedure presented in this paper (as applied on a set of compatibly shaped hydrographs) involves the following key stages: (i) normalize the floods, (ii) define similar flood shape classes and (iii) evaluate the average dimensionless flood (ADF) for each class. The ADFs are finally transformed into a set of (DF)s. Many statistical distributions approximate acceptably the frequent values of the maximum discharges or the flood volumes, yet displaying a significant spread for medium or rare probabilities of exceedance (PE). This scattering, which can be explained by the epistemic uncertainty, defines an area of uncertainty both for measured and extrapolated values. In considering upper and lower values of the uncertainty in-tervals as limits for maximum discharges and flood volumes, then by combining them compatibly, a set of DFs - as completely defined hydrographs, with different shapes - results for each PE. The herein proposed procedure defines both one peak DF and multi-peaks DF. Subsequently, such DFs do assist water managers in examining and establishing tailored approaches for a variety of input hydrographs. Among the DFs that would correspond to a same PE, the most compact floods arise a special interest, for they are basic in defining the set of safe operation rules for hydraulic structures.

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.