Refine
Year of publication
Document Type
- Article (280)
- Postprint (11)
- Preprint (9)
- Monograph/Edited Volume (8)
- Other (1)
Language
- English (309) (remove)
Keywords
- Complex networks (5)
- Event synchronization (4)
- precipitation (3)
- synchronization (3)
- Amazon rainforest (2)
- Extreme rainfall (2)
- Synchronization (2)
- channel (2)
- classification (2)
- climate networks (2)
Institute
- Institut für Physik und Astronomie (232)
- Interdisziplinäres Zentrum für Dynamik komplexer Systeme (47)
- Institut für Geowissenschaften (21)
- Department Psychologie (16)
- Institut für Biochemie und Biologie (4)
- Extern (3)
- Department Linguistik (2)
- Institut für Informatik und Computational Science (2)
- Mathematisch-Naturwissenschaftliche Fakultät (2)
- Department Sport- und Gesundheitswissenschaften (1)
We investigate the cognitive control in polyrhythmic hand movements as a model paradigm for bimanual coordination. Using a symbolic coding of the recorded time series, we demonstrate the existence of qualitative transitions induced by experimental manipulation of the tempo. A nonlinear model with delayed feedback control is proposed, which accounts for these dynamical transitions in terms of bifurcations resulting from variation of the external control parameter. Furthermore, it is shown that transitions can also be observed due to fluctuations in the timing control level. We conclude that the complexity of coordinated bimanual movements results from interactions between nonlinear control mechanisms with delayed feedback and stochastic timing components.
Recurrence-plot-based recurrence networks are an approach used to analyze time series using a complex networks theory. In both approaches - recurrence plots and recurrence networks -, a threshold to identify recurrent states is required. The selection of the threshold is important in order to avoid bias of the recurrence network results. In this paper, we propose a novel method to choose a recurrence threshold adaptively. We show a comparison between the constant threshold and adaptive threshold cases to study period-chaos and even period-period transitions in the dynamics of a prototypical model system. This novel method is then used to identify climate transitions from a lake sediment record.
We study Hamiltonian chaos generated by the dynamics of passive tracers moving in a two-dimensional fluid flow and describe the complex structure formed in a chaotic layer that separates a vortex region from the shear flow. The stable and unstable manifolds of unstable periodic orbits are computed. It is shown that their intersections in the Poincare map as an invariant set of homoclinic points constitute the backbone of the chaotic layer. Special attention is paid to the finite time properties of the chaotic layer. In particular, finite time Lyapunov exponents are computed and a scaling law of the variance of their distribution is derived. Additionally, the box counting dimension as an effective dimension to characterize the fractal properties of the layer is estimated for different duration times of simulation. Its behavior in the asymptotic time limit is discussed. By computing the Lyapunov exponents and by applying methods of symbolic dynamics, the formation of the layer as a function of the external forcing strength, which in turn represents the perturbation of the originally integrable system, is characterized. In particular, it is shown that the capture of KAM tori by the layer has a remarkable influence on the averaged Lyapunov exponents. (C) 2004 Elsevier Ltd. All rights reserved
Experimental evidences point Out the participation of nonsynaptic mechanisms (e.g., fluctuations in extracellular tons) in epileptiform bursting and spreading depression (SD). During these abnormal oscillatory patterns, it is observed an increase of extracellular potassium concentration [K+](o) and a decrease of extracellular calcium concentration [Ca2+](o) which raises the neuronal excitability. However, whether the high [K+](o) triggers and propagates these abnormal neuronal activities or plays a secondary role into this process is unclear. To better understand the influence of extracellular potassium dynamics in these oscillatory patterns, the experimental conditions of high [K+](o) and zero [Ca2+](o) were replicated in an extended Golomb model where we added important regulatory mechanisms of ion concentration as Na+-K+ pump, ion diffusion and glial buffering. Within these Conditions, simulations of the cell model exhibit seizure-like discharges (ictal bursting). The SD was elicited by the interruption of the Na+- K+ pump activity, mimicking the effect of cellular hypoxia (an experimental protocol to elicit SD, the hypoxia-induced SD). We used the bifurcation theory and the fast-slow method to analyze the interference of K+ dynamics in the cellular excitability. This analysis indicates that the system loses its stability at a high [K+](o), transiting to an elevated state of neuronal excitability. Effects of high [K+](o), are observed in different stages of ictal bursting and SD. In the initial stage, the increase of [K+](o) creates favorable conditions to trigger both oscillatory patterns. During the neuronal activity, a continuous growth of [K+](o) by outward K+ flow depresses K+ Currents in a positive feedback way. At the last stage, due to the depression of K+ currents, the Na+-K+ pump is the main mechanism in the end of neuronal activity. Thus, this work suggests that [K+](o) dynamics may play a fundamental role in these abnormal oscillatory patterns.
In this paper, we analytically study a star motif of Stuart-Landau oscillators, derive the bifurcation diagram and discuss the different forms of synchronization arising in such a system. Despite the parameter mismatch between the central node and the peripheral ones, an analytical approach independent of the number of units in the system has been proposed. The approach allows to calculate the separatrices between the regions with distinct dynamical behavior and to determine the nature of the different transitions to synchronization appearing in the system. The theoretical analysis is supported by numerical results.
We employ a spectral decomposition method to analyze synchronization of a non-identical oscillator network. We study the case that a small parameter mismatch of oscillators is characterized by one parameter and phase synchronization is observed. We derive a linearized equation for each eigenmode of the coupling matrix. The parameter mismatch is reflected on inhomogeneous term in the linearized equation. We find that the oscillation of each mode is essentially characterized only by the eigenvalue of the coupling matrix with a suitable normalization. We refer to this property as spectral universality, because it is observed irrespective of network topology. Numerical results in various network topologies show good agreement with those based on linearized equation. This universality is also observed in a system driven by additive independent Gaussian noise.
Linear methods of dimensionality reduction are useful tools for handling and interpreting high dimensional data. However, the cumulative variance explained by each of the subspaces in which the data space is decomposed may show a slow convergence that makes the selection of a proper minimum number of subspaces for successfully representing the variability of the process ambiguous. The use of nonlinear methods can improve the embedding of multivariate data into lower dimensional manifolds. In this article, a nonlinear method for dimensionality reduction, Isomap, is applied to the sea surface temperature and thermocline data in the tropical Pacific Ocean, where the El Nino-Southern Oscillation (ENSO) phenomenon and the annual cycle phenomena interact. Isomap gives a more accurate description of the manifold dimensionality of the physical system. The knowledge of the minimum number of dimensions is expected to improve the development of low dimensional models for understanding and predicting ENSO
We study the overdamped version of two coupled anharmonic oscillators under the influence of both low- and high-frequency forces respectively and a Gaussian noise term added to one of the two state variables of the system. The dynamics of the system is first studied in the presence of both forces separately without noise. In the presence of only one of the forces, no resonance behaviour is observed, however, hysteresis happens there. Then the influence of the high-frequency force in the presence of a low-frequency, i.e. biharmonic forcing, is studied. Vibrational resonance is found to occur when the amplitude of the high-frequency force is varied. The resonance curve resembles a stochastic resonance-like curve. It is maximum at the value of g at which the orbit lies in one well during one half of the drive cycle of the low-frequency force and in the other for the remaining half cycle. Vibrational resonance is characterized using the response amplitude and mean residence time. We show the occurrence of stochastic resonance behaviour in the overdamped system by replacing the high-frequency force by Gaussian noise. Similarities and differences between both types of resonance are presented. (c) 2006 Elsevier B.V. All rights reserved.
Identifying abrupt transitions is a key question in various disciplines. Existing transition detection methods, however, do not rigorously account for time series uncertainties, often neglecting them altogether or assuming them to be independent and qualitatively similar. Here, we introduce a novel approach suited to handle uncertainties by representing the time series as a time-ordered sequence of probability density functions. We show how to detect abrupt transitions in such a sequence using the community structure of networks representing probabilities of recurrence. Using our approach, we detect transitions in global stock indices related to well-known periods of politico-economic volatility. We further uncover transitions in the El Niño-Southern Oscillation which coincide with periods of phase locking with the Pacific Decadal Oscillation. Finally, we provide for the first time an ‘uncertainty-aware’ framework which validates the hypothesis that ice-rafting events in the North Atlantic during the Holocene were synchronous with a weakened Asian summer monsoon.
Identifying abrupt transitions is a key question in various disciplines. Existing transition detection methods, however, do not rigorously account for time series uncertainties, often neglecting them altogether or assuming them to be independent and qualitatively similar. Here, we introduce a novel approach suited to handle uncertainties by representing the time series as a time-ordered sequence of probability density functions. We show how to detect abrupt transitions in such a sequence using the community structure of networks representing probabilities of recurrence. Using our approach, we detect transitions in global stock indices related to well-known periods of politico-economic volatility. We further uncover transitions in the El Nino-Southern Oscillation which coincide with periods of phase locking with the Pacific Decadal Oscillation. Finally, we provide for the first time an 'uncertainty-aware' framework which validates the hypothesis that ice-rafting events in the North Atlantic during the Holocene were synchronous with a weakened Asian summer monsoon.
Sedimentary proxy records constitute a significant portion of the recorded evidence that allows us to investigate paleoclimatic conditions and variability. However, uncertainties in the dating of proxy archives limit our ability to fix the timing of past events and interpret proxy record intercomparisons. While there are various age-modeling approaches to improve the estimation of the age-depth relations of archives, relatively little focus has been placed on the propagation of the age (and radiocarbon calibration) uncertainties into the final proxy record.
We present a generic Bayesian framework to estimate proxy records along with their associated uncertainty, starting with the radiometric age-depth and proxy-depth measurements, and a radiometric calibration curve if required. We provide analytical expressions for the posterior proxy probability distributions at any given calendar age, from which the expected proxy values and their uncertainty can be estimated. We illustrate our method using two synthetic data sets and then use it to construct the proxy records for groundwater inflow and surface erosion from Lonar lake in central India.
Our analysis reveals interrelations between the uncertainty of the proxy record over time and the variance of proxies along the depth of the archive. For the Lonar lake proxies, we show that, rather than the age uncertainties, it is the proxy variance combined with calibration uncertainty that accounts for most of the final uncertainty. We represent the proxy records as probability distributions on a precise, error-free timescale that makes further time series analyses and intercomparisons of proxies relatively simple and clear. Our approach provides a coherent understanding of age uncertainties within sedimentary proxy records that involve radiometric dating. It can be potentially used within existing age modeling structures to bring forth a reliable and consistent framework for proxy record estimation.
We investigate a network of influences connected to global mean temperature. Considering various climatic factors known to influence global mean temperature, we evaluate not only the impacts of these factors on temperature but also the directed dependencies among the factors themselves. Based on an existing recurrence-based connectivity measure, we propose a new and more general measure that quantifies the level of dependence between two time series based on joint recurrences at a chosen time delay. The measures estimated in the analysis are tested for statistical significance using twin surrogates. We find, in accordance with earlier studies, the major drivers for global mean temperature to be greenhouse gases, ENSO, volcanic activity, and solar irradiance. We further uncover a feedback between temperature and ENSO. Our results demonstrate the need to involve multiple, delayed interactions within the drivers of temperature in order to develop a more thorough picture of global temperature variations.
We propose a RAndom Interacting Network (RAIN) model to study the interactions between a pair of complex networks. The model involves two major steps: (i) the selection of a pair of nodes, one from each network, based on intra-network node-based characteristics, and (ii) the placement of a link between selected nodes based on the similarity of their relative importance in their respective networks. Node selection is based on a selection fitness function and node linkage is based on a linkage probability defined on the linkage scores of nodes. The model allows us to relate within-network characteristics to between-network structure. We apply the model to the interaction between the USA and Schengen airline transportation networks (ATNs). Our results indicate that two mechanisms: degree-based preferential node selection and degree-assortative link placement are necessary to replicate the observed inter-network degree distributions as well as the observed inter-network assortativity. The RAIN model offers the possibility to test multiple hypotheses regarding the mechanisms underlying network interactions. It can also incorporate complex interaction topologies. Furthermore, the framework of the RAIN model is general and can be potentially adapted to various real-world complex systems.
We describe effects of the asymmetry of cycles and non-stationarity in time series on the phase synchronization method which may lead to artifacts. We develop a modified method that overcomes these effects and apply it to study parkinsonian tremor. Our results indicate that there is synchronization between two different hands and provide information about the time delay separating their dynamics. These findings suggest that this method may be useful for detecting and quantifying weak synchronization between two non-stationary signals.
Spatiotemporal dynamics of the Calvin cycle multistationarity and symmetry breaking instabilities
(2011)
The possibility of controlling the Calvin cycle has paramount implications for increasing the production of biomass. Multistationarity, as a dynamical feature of systems, is the first obvious candidate whose control could find biotechnological applications. Here we set out to resolve the debate on the multistationarity of the Calvin cycle. Unlike the existing simulation-based studies, our approach is based on a sound mathematical framework, chemical reaction network theory and algebraic geometry, which results in provable results for the investigated model of the Calvin cycle in which we embed a hierarchy of realistic kinetic laws. Our theoretical findings demonstrate that there is a possibility for multistationarity resulting from two sources, homogeneous and inhomogeneous instabilities, which partially settle the debate on multistability of the Calvin cycle. In addition, our tractable analytical treatment of the bifurcation parameters can be employed in the design of validation experiments.
Contents: 1 Introduction 2 Experiment 3 Data 4 Symbolic dynamics 4.1 Symbolic dynamics as a tool for data analysis 4.2 2-symbols coding 4.3 3-symbols coding 5 Measures of complexity 5.1 Word statistics 5.2 Shannon entropy 6 Testing for stationarity 6.1 Stationarity 6.2 Time series of cycle durations 6.3 Chi-square test 7 Control parameters in the production of rhythms 8 Analysis of relative phases 9 Discussion 10 Outlook
Similar power laws for foreshock and aftershock sequences in a spring block model for earthquakes
(1999)
Seismic quiescence as an indicator for large earthquakes in a system of self-organized criticality
(2000)
Competitive random sequential adsorption of point and fixed-sized particles: analytical results
(2001)
Objectives: Scoring sleep visually based on polysomnography is an important but time-consuming element of sleep medicine. Where-as computer software assists human experts in the assignment of sleep stages to polysomnogram epochs, their performance is usually insufficient. This study evaluates the possibility to fully automatize sleep staging considering the reliability of the sleep stages available from human expert sleep scorers. Methods: We obtain features from EEG, ECG and respiratory signals of polysomnograms from ten healthy subjects. Using the sleep stages provided by three human experts, we evaluate the performance of linear discriminant analysis on the entire polysomnogram and:only on epochs where the three experts agree in their-sleep stage scoring. Results: We show that in polysomnogram intervals, to which all three scorers assign the same sleep stage, our algorithm achieves 90% accuracy. This high rate of agreement with the human experts is accomplished with only a small set of three frequency features from the EEG. We increase-the performance to 93% by including ECG and respiration features. In contrast, on intervals of ambiguous sleep stage, the sleep stage classification obtained from our algorithm, agrees with the human consensus scorer in approximately 61%. Conclusions: These findings suggest that machine classification is highly consistent with human sleep staging and that error in the algorithm's assignments is rather a problem of lack of well-defined criteria for human experts to judge certain polysomnogram epochs than an insufficiency of computational procedures
Background: Inferring regulatory interactions between genes from transcriptomics time-resolved data, yielding reverse engineered gene regulatory networks, is of paramount importance to systems biology and bioinformatics studies. Accurate methods to address this problem can ultimately provide a deeper insight into the complexity, behavior, and functions of the underlying biological systems. However, the large number of interacting genes coupled with short and often noisy time-resolved read-outs of the system renders the reverse engineering a challenging task. Therefore, the development and assessment of methods which are computationally efficient, robust against noise, applicable to short time series data, and preferably capable of reconstructing the directionality of the regulatory interactions remains a pressing research problem with valuable applications.
Results: Here we perform the largest systematic analysis of a set of similarity measures and scoring schemes within the scope of the relevance network approach which are commonly used for gene regulatory network reconstruction from time series data. In addition, we define and analyze several novel measures and schemes which are particularly suitable for short transcriptomics time series. We also compare the considered 21 measures and 6 scoring schemes according to their ability to correctly reconstruct such networks from short time series data by calculating summary statistics based on the corresponding specificity and sensitivity. Our results demonstrate that rank and symbol based measures have the highest performance in inferring regulatory interactions. In addition, the proposed scoring scheme by asymmetric weighting has shown to be valuable in reducing the number of false positive interactions. On the other hand, Granger causality as well as information-theoretic measures, frequently used in inference of regulatory networks, show low performance on the short time series analyzed in this study.
Conclusions: Our study is intended to serve as a guide for choosing a particular combination of similarity measures and scoring schemes suitable for reconstruction of gene regulatory networks from short time series data. We show that further improvement of algorithms for reverse engineering can be obtained if one considers measures that are rooted in the study of symbolic dynamics or ranks, in contrast to the application of common similarity measures which do not consider the temporal character of the employed data. Moreover, we establish that the asymmetric weighting scoring scheme together with symbol based measures (for low noise level) and rank based measures (for high noise level) are the most suitable choices.
Background: Inferring regulatory interactions between genes from transcriptomics time-resolved data, yielding reverse engineered gene regulatory networks, is of paramount importance to systems biology and bioinformatics studies. Accurate methods to address this problem can ultimately provide a deeper insight into the complexity, behavior, and functions of the underlying biological systems. However, the large number of interacting genes coupled with short and often noisy time-resolved read-outs of the system renders the reverse engineering a challenging task. Therefore, the development and assessment of methods which are computationally efficient, robust against noise, applicable to short time series data, and preferably capable of reconstructing the directionality of the regulatory interactions remains a pressing research problem with valuable applications.
Results: Here we perform the largest systematic analysis of a set of similarity measures and scoring schemes within the scope of the relevance network approach which are commonly used for gene regulatory network reconstruction from time series data. In addition, we define and analyze several novel measures and schemes which are particularly suitable for short transcriptomics time series. We also compare the considered 21 measures and 6 scoring schemes according to their ability to correctly reconstruct such networks from short time series data by calculating summary statistics based on the corresponding specificity and sensitivity. Our results demonstrate that rank and symbol based measures have the highest performance in inferring regulatory interactions. In addition, the proposed scoring scheme by asymmetric weighting has shown to be valuable in reducing the number of false positive interactions. On the other hand, Granger causality as well as information-theoretic measures, frequently used in inference of regulatory networks, show low performance on the short time series analyzed in this study.
Conclusions: Our study is intended to serve as a guide for choosing a particular combination of similarity measures and scoring schemes suitable for reconstruction of gene regulatory networks from short time series data. We show that further improvement of algorithms for reverse engineering can be obtained if one considers measures that are rooted in the study of symbolic dynamics or ranks, in contrast to the application of common similarity measures which do not consider the temporal character of the employed data. Moreover, we establish that the asymmetric weighting scoring scheme together with symbol based measures (for low noise level) and rank based measures (for high noise level) are the most suitable choices.
Identifying causal links (couplings) is a fundamental problem that facilitates the understanding of emerging structures in complex networks. We propose and analyze inner composition alignment-a novel, permutation-based asymmetric association measure to detect regulatory links from very short time series, currently applied to gene expression. The measure can be used to infer the direction of couplings, detect indirect (superfluous) links, and account for autoregulation. Applications to the gene regulatory network of E. coli are presented.
In this paper, we study the complete synchronization of a class of time-varying delayed coupled chaotic systems using feedback control. In terms of Linear Matrix Inequalities, a sufficient condition is obtained through using a Lyapunov-Krasovskii functional and differential equation in equalities. The conditions can be easily verified and implemented. We present two simulation examples to illustrate the effectiveness of the proposed method.
In one of the data mining techniques, change-point detection is of importance in evaluating time series measured in real world. For decades this technique has been developed as a nonlinear dynamics. We apply the method for detecting the change points, Singular Spectrum Transformation (SST), to the climate time series. To know where the structures of climate data sets change can reveal a climate background. In this paper we discuss the structures of precipitation data in Kenya and Wrangel Island (Arctic land) by using the SST.
We study the effects of mutual and external chaotic phase synchronization in ensembles of bursting oscillators. These oscillators (used for modeling neuronal dynamics) are essentially multiple time scale systems. We show that a transition to mutual phase synchronization takes place on the bursting time scale of globally coupled oscillators, while on the spiking time scale, they behave asynchronously. We also demonstrate the effect of the onset of external chaotic phase synchronization of the bursting behavior in the studied ensemble by a periodic driving applied to one arbitrarily taken neuron. We also propose an explanation of the mechanism behind this effect. We infer that the demonstrated phenomenon can be used efficiently for controlling bursting activity in neural ensembles
We study phase synchronization effects of chaotic oscillators with a type-I intermittency behavior. The external and mutual locking of the average length of the laminar stage for coupled discrete and continuous in time systems is shown and the mechanism of this synchronization is explained. We demonstrate that this phenomenon can be described by using results of the parametric resonance theory and that this correspondence enables one to predict and derive all zones of synchronization
Being one of the fundamental phenomena in nonlinear science, synchronization of oscillations has permanently remained an object of intensive research. Development of many asymptotic methods and numerical simulations has allowed an understanding and explanation of various phenomena of self-synchronization. But even in the classical case of coupled van der Pol oscillators a full description of all possible dynamical regimes, their mutual transitions and characteristics is still lacking. We present here a study of the phenomenon of mutual synchronization for two non-scalar- coupled non-identical limit-cycle oscillators and analyze phase, frequency and amplitude characteristics of synchronization regimes. A series of bifurcation diagrams that we obtain exhibit various regions of qualitatively different behavior. Among them we find mono-, bi- and multistability regions, beating and "oscillation death" ones; also a region, where one of the oscillators dominates the other one is observed. The frequency characteristics that we obtain reveal three qualitatively different types of synchronization: (i) on the mean frequency (the in-phase synchronization), (ii) with a shift from the mean frequency caused by a conservative coupling term (the anti-phase synchronization), and (iii) on the frequency of one of the oscillators (when one oscillator dominates the other). (C) 2003 Elsevier B.V. All rights reserved
We investigate the influence of noise on synchronization between the spiking activities of neurons with external impulsive forces. We first analyze the dependence of the synchronized firing on the amplitude and the angular frequency of the impulsive force in the noise-free system. Three cases (regular spiking, traveling wave, and chaotic spiking) with low synchronized firing are chosen to study effects due to noise. In each case we find that small noise can be a promoter of synchronization phenomena in neural activities, by choosing an appropriate noise intensity acting on some of the neurons. (C) 2005 American Institute of Physics
Synthetic multicellular oscillatory systems controlling protein dynamics with genetic circuits
(2011)
Synthetic biology is a relatively new research discipline that combines standard biology approaches with the constructive nature of engineering. Thus, recent efforts in the field of synthetic biology have given a perspective to consider cells as 'programmable matter'. Here, we address the possibility of using synthetic circuits to control protein dynamics. In particular, we show how intercellular communication and stochasticity can be used to manipulate the dynamical behavior of a population of coupled synthetic units and, in this manner, finely tune the expression of specific proteins of interest, e.g. in large bioreactors.
Detuning-dependent dominance of oscillation death in globally coupled synthetic genetic oscillators
(2009)
We study dynamical regimes of globally coupled genetic relaxation oscillators in the presence of small detuning. Using bifurcation analysis, we find that under strong coupling via the slow variable, the detuning can eliminate standard oscillatory solutions in a large region of the parameter space, providing the dominance of oscillation death. This result is substantially different from previous results on oscillation quenching, where for homogeneous populations, the coexistence of oscillation death and limit cycle oscillations is always present. We propose further that this effect of detuning-dependent dominance could be a powerful regulator of genetic network's dynamics.
Many cellular processes require decision making mechanisms, which must act reliably even in the unavoidable presence of substantial amounts of noise. However, the multistable genetic switches that underlie most decision-making processes are dominated by fluctuations that can induce random jumps between alternative cellular states. Here we show, via theoretical modeling of a population of noise-driven bistable genetic switches, that reliable timing of decision-making processes can be accomplished for large enough population sizes, as long as cells are globally coupled by chemical means. In the light of these results, we conjecture that cell proliferation, in the presence of cell-cell communication, could provide a mechanism for reliable decision making in the presence of noise, by triggering cellular transitions only when the whole cell population reaches a certain size. In other words , the summation performed by the cell population would average out the noise and reduce its detrimental impact.
A better understanding of precipitation dynamics in the Indian subcontinent is required since India's society depends heavily on reliable monsoon forecasts. We introduce a non-linear, multiscale approach, based on wavelets and event synchronization, for unravelling teleconnection influences on precipitation. We consider those climate patterns with the highest relevance for Indian precipitation. Our results suggest significant influences which are not well captured by only the wavelet coherence analysis, the state-of-the-art method in understanding linkages at multiple timescales. We find substantial variation across India and across timescales. In particular, El Niño–Southern Oscillation (ENSO) and the Indian Ocean Dipole (IOD) mainly influence precipitation in the south-east at interannual and decadal scales, respectively, whereas the North Atlantic Oscillation (NAO) has a strong connection to precipitation, particularly in the northern regions. The effect of the Pacific Decadal Oscillation (PDO) stretches across the whole country, whereas the Atlantic Multidecadal Oscillation (AMO) influences precipitation particularly in the central arid and semi-arid regions. The proposed method provides a powerful approach for capturing the dynamics of precipitation and, hence, helps improve precipitation forecasting.