Search

As of November 2020, the number of COVID-19 cases is increasing rapidly in many countries. In Europe, the virus spread slowed considerably in the late spring due to strict lockdown, but a second wave of the pandemic grew throughout the fall. In this study, we first reconstruct the time evolution of the effective reproduction numbers ${\cal R}(t)$ for each country by integrating the equations of the classic SIR model. We cluster countries based on the estimated ${\cal R}(t)$ through a suitable time series dissimilarity. The result suggests that simple dynamical mechanisms determine how countries respond to changes in COVID-19 case counts. Inspired by these results, we extend the SIR model to include a social response to explain the number $X(t)$ of new confirmed daily cases. As a first-order model, we assume that the social response is on the form $d_t{\cal R}=-\nu (X-X^*)$, where $X^*$ is a threshold for response. The response rate $\nu$ depends on whether $X^*$ is below or above this threshold, on three parameters $\nu_1,\;\nu_2,\,\nu_3,$, and on $t$. When $XX^*$, $\nu=\nu_2\exp{(-\nu_3t)}$, models the impact of interventions when incidence rate is high. The parameter $\nu_3$ represents the fatigue, i.e., the reduced effect of intervention as time passes. The proposed model reproduces typical evolving patterns of COVID-19 epidemic waves observed in many countries. Estimating the parameters $\nu_1,\,\nu_2,\,\nu_3$ and initial conditions, such as ${\cal R}_0$, for different countries helps to identify important dynamics in their social responses. One conclusion is that the leading cause of the strong second wave in Europe in the fall of 2020 was not the relaxation of interventions during the summer, but rather the general fatigue to interventions developing in the fall., Comment: 19 pages, 6 figures

In a given country, the cumulative death toll of the first wave of the COVID-19 epidemic follows a sigmoid curve as a function of time. In most cases, the curve is well described by the Gompertz function, which is characterized by two essential parameters, the initial growth rate and the decay rate as the first epidemic wave subsides. These parameters are determined by socioeconomic factors and the countermeasures to halt the epidemic. The Gompertz model implies that the total death toll depends exponentially, and hence very sensitively, on the ratio between these rates. The remarkably different epidemic curves for the first epidemic wave in Sweden and Norway and many other countries are classified and discussed in this framework, and their usefulness for the planning of mitigation strategies is discussed., Comment: 27 pages, 14 figures

The temporal fluctuations in global mean surface temperature is an example of a geophysical quantity which can be described using the notions of long-range persistence and scale invariance/scaling, but this description has suffered from lack of a generally accepted physical explanation. Processes with these statistical signatures can arise from non-linear effects, for instance through cascade-like energy transfer in turbulent fluids, but they can also be produced by linear models with scale-invariant impulse-response functions. This paper demonstrates that on time scales from months to centuries, the scale-invariant impulse-response function of global surface temperature can be explained from simple linear multi-box energy balance models. This explanation describes both the scale invariance of the internal variability and the lack of a characteristic time scale of the response to external forcings. With parameters estimated from observational data, the climate response is approximately scaling in these models, even if the response function is not chosen to be scaling a priori. It is also demonstrated that the differences in scaling exponents for temperatures over land and for sea-surface temperatures can be reproduced by a version of the multi-box energy balance model with two distinct surface boxes.

Filtered Poisson processes are often used as reference models for intermittent fluc- tuations in physical systems. Such a process is here extended by adding a noise term, either as a purely additive term to the process or as a dynamical term in a stochastic differential equation. The lowest order moments, probability density function, auto-correlation function and power spectral density are derived and used to identify and compare the effects of the two different noise terms. Monte-Carlo studies of synthetic time series are used to investigate the accuracy of model pa- rameter estimation and to identify methods for distinguishing the noise types. It is shown that the probability density function and the three lowest order moments provide accurate estimations of the parameters, but are unable to separate the noise types. The auto-correlation function and the power spectral density also provide methods for estimating the model parameters, as well as being capable of identifying the noise type. The number of times the signal crosses a prescribed threshold level in the positive direction also promises to be able to differentiate the noise type., Comment: 34 pages, 25 figures

Various interpretations of the notion of a trend in the context of global warming are discussed, contrasting the difference between viewing a trend as the deterministic response to an external forcing and viewing it as a slow variation which can be separated from the background spectral continuum of long-range persistent climate noise. The emphasis in this paper is on the latter notion, and a general scheme is presented for testing a multi-parameter trend model against a null hypothesis which models the observed climate record as an autocorrelated noise. The scheme is employed to the instrumental global sea-surface temperature record and the global land-temperature record. A trend model comprising a linear plus an oscillatory trend with period of approximately 60 yr, and the statistical significance of the trends, are tested against three different null models: first-order autoregressive process, fractional Gaussian noise, and fractional Brownian motion. The linear trend is significant in all cases, but the oscillatory trend is insignificant for ocean data and barely significant for land data. By means of a Bayesian iteration, however, using the significance of the linear trend to formulate a sharper null hypothesis, the oscillatory trend in the land record appears to be statistically significant. The results suggest that the global land record may be better suited for detection of the global warming signal than the ocean record., Comment: 7 pages, 7 figures

A linearized energy-balance model for global temperature is formulated, featuring a scale-free long-range memory (LRM) response and stochastic forcing representing the influence on the ocean heat reservoir from atmospheric weather systems. The model is parametrized by an effective response strength, the stochastic forcing strength, and the memory exponent. The instrumental global surface temperature record and the deterministic component of the forcing are used to estimate these parameters by means of the maximum-likelihood method. The residual obtained by subtracting the deterministic solution from the observed record is analyzed as a noise process and shown to be consistent with a long-memory time-series model and inconsistent with a short-memory model. By decomposing the forcing record in contributions from solar, volcanic, and anthropogenic activity one can estimate the contribution of each to 20'th century global warming. The LRM model is applied with a reconstruction of the forcing for the last millennium to predict the large-scale features of northern hemisphere temperature reconstructions, and the analysis of the residual also clearly favors the LRM model on millennium time scale. The decomposition of the forcing shows that volcanic aerosols give a considerably greater contribution to the cooling during the Little Ice Age than the reduction in solar irradiance associated with the Maunder minimum in solar activity. The LRM model implies a transient climate response in agreement with IPCC AR4 projections, but the stronger response on longer time scales suggests to replace the notion of equilibrium climate sensitivity by a time-scale dependent sensitivity.

We introduce tools for inference in the multifractal random walk introduced by Bacry et al. (2001). These tools include formulas for smoothing, filtering and volatility forecasting. In addition, we present methods for computing conditional densities for one- and multi-step returns. The inference techniques presented in this paper, including maximum likelihood estimation, are applied to data from the Oslo Stock Exchange, and it is observed that the volatility forecasts based on the multifractal random walk have a much richer structure than the forecasts obtained from a basic stochastic volatility model., Comment: 8 pages, 3 figures

Maximum likelihood estimation applied to high-frequency data allows us to quantify intermittency in the fluctu- ations of asset prices. From time records as short as one month these methods permit extraction of a meaningful intermittency parameter {\lambda} characterising the degree of volatility clustering of asset prices. We can therefore study the time evolution of volatility clustering and test the statistical significance of this variability. By analysing data from the Oslo Stock Exchange, and comparing the results with the investment grade spread, we find that the estimates of {\lambda} are lower at times of high market uncertainty., Comment: 8 pages, 4 figures, 2 tables

We discuss stochastic modeling of volatility persistence and anti-correlations in electricity spot prices, and for this purpose we present two mean-reverting versions of the multifractal random walk (MRW). In the first model the anti-correlations are modeled in the same way as in an Ornstein-Uhlenbeck process, i.e. via a drift (damping) term, and in the second model the anti-correlations are included by letting the innovations in the MRW model be fractional Gaussian noise with H < 1/2. For both models we present approximate maximum likelihood methods, and we apply these methods to estimate the parameters for the spot prices in the Nordic electricity market. The maximum likelihood estimates show that electricity spot prices are characterized by scaling exponents that are significantly different from the corresponding exponents in stock markets, confirming the exceptional nature of the electricity market. In order to compare the damped MRW model with the fractional MRW model we use ensemble simulations and wavelet-based variograms, and we observe that certain features of the spot prices are better described by the damped MRW model. The characteristic correlation time is estimated to approximately half a year., Comment: 13 pages, 4 figures, 2 tables

We present an approximated maximum likelihood method for the multifractal random walk processes of [E. Bacry et al., Phys. Rev. E 64, 026103 (2001)]. The likelihood is computed using a Laplace approximation and a truncation in the dependency structure for the latent volatility. The procedure is implemented as a package in the R computer language. Its performance is tested on synthetic data and compared to an inference approach based on the generalized method of moments. The method is applied to estimate parameters for various financial stock indices., Comment: 8 pages, 3 figures, 2 tables

The sunspot number (SSN), the total solar irradiance (TSI), a TSI reconstruction, and the solar flare index (SFI), are analyzed for long-range persistence (LRP). Standard Hurst analysis yields $H \approx 0.9$, which suggests strong LRP. However, solar activity time series are non-stationary due to the almost periodic 11 year smooth component, and the analysis does not give the correct $H$ for the stochastic component. Better estimates are obtained by detrended fluctuations analysis (DFA), but estimates are biased and errors are large due to the short time records. These time series can be modeled as a stochastic process of the form $x(t)=y(t)+\sigma \sqrt{y(t)}\, w_H(t)$, where $y(t)$ is the smooth component, and $w_H(t) $ is a stationary fractional noise with Hurst exponent $H$. From ensembles of numerical solutions to the stochastic model, and application of Bayes' theorem, we can obtain bias and error bars on $H$ and also a test of the hypothesis that a process is uncorrelated ($H=1/2$). The conclusions from the present data sets are that SSN, TSI and TSI reconstruction almost certainly are long-range persistent, but with most probable value $H\approx 0.7$. The SFI process, however, is either very weakly persistent ($H<0.6$) or completely uncorrelated. Some differences between stochastic properties of the TSI and its reconstruction indicate some error in the reconstruction scheme., Comment: 23 pages, 12 figures

We demonstrate how to model the toppling activity in avalanching systems by stochastic differential equations (SDEs). The theory is developed as a generalization of the classical mean field approach to sandpile dynamics by formulating it as a generalization of Itoh's SDE. This equation contains a fractional Gaussian noise term representing the branching of an avalanche into small active clusters, and a drift term reflecting the tendency for small avalanches to grow and large avalanches to be constricted by the finite system size. If one defines avalanching to take place when the toppling activity exceeds a certain threshold the stochastic model allows us to compute the avalanche exponents in the continum limit as functions of the Hurst exponent of the noise. The results are found to agree well with numerical simulations in the Bak-Tang-Wiesenfeld and Zhang sandpile models. The stochastic model also provides a method for computing the probability density functions of the fluctuations in the toppling activity itself. We show that the sandpiles do not belong to the class of phenomena giving rise to universal non-Gaussian probability density functions for the global activity. Moreover, we demonstrate essential differences between the fluctuations of total kinetic energy in a two-dimensional turbulence simulation and the toppling activity in sandpiles., Comment: 14 pages, 11 figures

A stochastic theory for the toppling activity in sandpile models is developed, based on a simple mean-field assumption about the toppling process. The theory describes the process as an anti-persistent Gaussian walk, where the diffusion coefficient is proportional to the activity. It is formulated as a generalization of the It\^{o} stochastic differential equation with an anti-persistent fractional Gaussian noise source. An essential element of the theory is re-scaling to obtain a proper thermodynamic limit, and it captures all temporal features of the toppling process obtained by numerical simulation of the Bak-Tang-Wiesenfeld sandpile in this limit., Comment: 9 pages, 4 figures

The topological entropy of piecewise affine maps is studied. It is shown that singularities may contribute to the entropy only if there is angular expansion and we bound the entropy via the expansion rates of the map. As a corollary we deduce that non-expanding conformal piecewise affine maps have zero topological entropy. We estimate the entropy of piecewise affine skew-products. Examples of abnormal entropy growth are provided., Comment: Examples refined

We construct the simplest chaotic system with a two-point attractor., Comment: 2 pages

We give a detailed study of dynamical properties of the Zhang model, including evaluation of topological entropy and estimates for the Lyapunov exponents and the dimension of the attractor. In the thermodynamic limit the entropy goes to zero and the Lyapunov spectrum collapses., Comment: 60 pages, 17 figures