---
abstract: 'We propose to construct cross and joint ordinal pattern transition networks from multivariate time series for two coupled systems, where synchronizations are often present. In particular, we focus on phase synchronization, which is one of prototypical scenarios in dynamical systems. We systematically show that cross and joint ordinal patterns transition networks are sensitive to phase synchronization. Furthermore, we find that some particular missing ordinal patterns play crucial roles in forming the detailed structures in the parameter space whereas the calculations of permutation entropy measures often do not. We conclude that cross and joint ordinal partition transition network approaches provide complementary insights to the traditional symbolic analysis of synchronization transitions.'
address: 'Department of Physics, East China Normal University, Shanghai, 200062, China'
author:
- Heng Guo
- Jiayang Zhang
- Yong Zou
- Shuguang Guan
title: Cross and joint ordinal partition transition networks for multivariate time series analysis
---

nonlinear time series analysis, complex networks, ordinal pattern partition, transition network, phase synchronization

Introduction
============

Complex network theory provides an important paradigm for understanding the structural properties of complex systems composed of different interacting entities. During the last decade, various complex network approaches have been proposed to extract useful insights from time series data [@Donner2011IJBC]. Depending on slightly different definitions of nodes and links of network representations for time series, there are methods based on recurrences [@Donner2010NJP; @Marwan2009], visibility conditions [@Lacasa2008], cycles detections [@zhangPRL2006] and correlation networks [@yangPhA2008]. Some successful applications of these methods include time series from climate, sunspots, oil-water flows, and financial market [@Donner2011IJBC; @Donges2011PNAS; @Zou2012bChaos; @gaoChaos2017; @Elsner2009; @zouNJP2014; @zouNPG2014; @ZhangCNSNS2017; @CzechowskiChaos2016; @ZhangFOP2017a; @JiangFOP2017]. The consideration of a time series as a complex network allows a reinterpretation of many network theoretic measures in terms of characteristic phase space properties of a dynamical system [@DonnerEPJB2011].

Recently, a growing number of works are focused on transforming time series into networks by ordinal partitions of time series [@MichaelChaos2015; @KulpChaos2016; @KulpChaos2016b; @McCulloughChaos2016; @SakellariouChaos2016]. The basic idea of ordinal partition network method can be traced back to identifying ordinal patterns of time series [@BandtPRL2002; @ParlitzSPJST2013; @AmigoPTRSA2014]. Given a one-dimensional time series $\{ x(t)\}_{t=1, \cdots, L}$ comprising of $L$ points from a dynamical system, we first reconstruct its phase space by time delay embedding technique which yields $\vec{x}(t) = [x(t), x(t+\tau), \cdots, x(t+(D_x-1)\tau)]$, where $D_x$ and $\tau$ are embedding dimension and delay respectively [@Takens1981; @Kantz97]. The next step is to compute the rank order of each embedded vector $[x(t), x(t+\tau), \cdots, x(t+(D_x-1)\tau)]$, which is conveniently denoted by a symbol $\pi_x(t)$. When sliding windows from $t=1$ to $N = L - (D_x - 1)\tau$ in the embedded space, a symbolic representation of the trajectory $\pi_x(t)$ is produced. Following the symbolic representation, one traditional approach is to compute permutation entropy $\mathcal{H}_O$ based on the frequency plot of order patterns. Generally speaking, for a time series generated by a stochastic process for $N \to \infty$, it is known that all $D_x!$ patterns almost occur with equal probabilities, which yields the maximal value of $\mathcal{H}_O$. However, for a time series produced by deterministic dynamics, the frequency plot of the $D_x!$ patterns is not uniform, which leads to reduced value of $\mathcal{H}_O$. In some cases, a set of patterns may never occur and these missing patterns are often called forbidden patterns, which provide important information for quantifying determinism in time series data [@AmigoEPL2007; @AmigoEPL2008; @RossoEPJB2012; @RossoPhyA2012; @KulpChaos2014]. Then, the level of determinism of a time series from a particular dynamical system may be suggested by permutation entropy $\mathcal{H}_O$, which consists in very well established statistical measures in nonlinear time series analysis [@BandtPRL2002; @AmigoPTRSA2014; @politiPRL2017]. Some applications include characterizations of the difference between healthy and patients from EEG data [@MichaelChaos2015; @KulpChaos2016]. It is worth noting that complications may arise in real time analysis because missing ordinal patterns might be related to finite time length during the period of observation and correlated stochastic processes, which require some revised methods for the detection of determinism in relatively short noisy data [@AmigoEPL2007; @AmigoEPL2008; @RossoEPJB2012; @RossoPhyA2012; @KulpChaos2014].

The primary idea of ordinal partition transition network takes into account the inhomogeneous evolutionary behavior among the ordinal patterns [@MichaelChaos2015; @KulpChaos2016], which provides complementary information on the standard ordinal symbolic analysis of time series. Most of these works have focused on univariate time series $\{x(t)\}$. Embedding parameters $D_x$ and $\tau$ have crucial impacts on the resulting ordinal partition transition networks, especially for forbidden patterns [@KulpChaos2016b; @McCulloughChaos2016; @SakellariouChaos2016]. In addition, multivariate time series are ubiquitous in nature, ranging from stock markets and climate sciences. In a recent work, we proposed to construct ordinal partition transition networks from multivariate data [@zhangSciRep2017]. The resulting network is a directed and weighted network characterizing the pattern transition properties of time series in its associated velocity space. This novel approach has been successfully applied to capture phase coherence to non-phase coherence transitions and to characterize paths to phase synchronization, showing complementary insight to the traditional symbolic analysis of nonlinear time series analysis.

In this work, we further extend these ideas [@zhangSciRep2017] to construct cross and joint ordinal partition transition networks for two coupled systems. We show that both cross and joint ordinal pattern transition networks are able to capture synchronization transitions, in particular, focusing on the transitions to phase synchronization (PS) which is one of paradigmatic types in synchronization phenomena [@Pikovsky_Kurths_synchr; @Osipov2003; @ZhangFOP2017; @ChenFOP2017; @HuangFOP2016; @YingFOP2017]. The outline of this paper is as follows: first, we illustrate cross and joint network construction approaches and then introduce two entropy measures to quantify the inhomogeneous frequencies of ordinal patterns and their transitions in Sec. \[secCOPTJ\]. We apply these two entropy measures to characterize the synchronization transitions with both unidirectional and bidirectional coupling schemes in Sec. \[secSync\], and then some conclusions are drawn in Sec. \[secCons\].

Cross and joint ordinal partition transition networks {#secCOPTJ}
=====================================================

Network constructions {#secRungeK}
---------------------

The ordinal partition transition networks are illustrated by the following coupled Rössler systems [@zhengPRE2000; @mamenPLA2004], which shows transition scenarios to phase synchronization for various coupling schemes. In particular, PS has been observed for both unidirectional and bidirectional couplings. Note that the network construction method is not restricted by the following model. The ordinary differential equations (ODEs) of the coupled system read $$\begin{aligned}
\label{eq_roessler}
\begin{cases}
\dot{x}_1 &= -(1.0 - \Delta) y_1 - z_1 + \kappa_1 (x_2 - x_1), \\
\dot{y}_1 &= (1.0 - \Delta)x_1 + 0.15 y_1, \\
\dot{z}_1 &= 0.2 + z_1 (x_1 - 10.0),
\end{cases}\\
\begin{cases}
\dot{x}_2 &= -(1.0 + \Delta) y_2 - z_2 + \kappa_2 (x_1 - x_2), \\
\dot{y}_2 &= (1.0 + \Delta) x_2 + 0.15 y_2, \\
\dot{z}_2 &= 0.2 + z_2 (x_2 - 10.0),
\end{cases}\end{aligned}$$ where $\omega_{1} = 1.0 - \Delta, \omega_2 = 1.0 + \Delta$ are frequencies for the systems, and $\kappa_{1,2}$ are coupling strength. In this model, we distinguish a unidirectional coupling case from a bidirectional coupling one because the threshold values to PS have been reported in the literature [@zhengPRE2000; @mamenPLA2004]. More specifically the unidirectional coupling can be achieved by $\kappa_1 = 0$, namely, a drive-response scheme [@zhengPRE2000] and bidirectional coupling is often chosen as $\kappa_1 = \kappa_2 = \kappa$ [@mamenPLA2004]. We numerically integrate the ODEs using the fourth-order Runge-Kutta method with random initial conditions, and the integration step $h = 0.01$. The first $10000$ transient data points are discarded and time series consisting of $N=800000$ data points are analyzed.

[*Ordinal Pattern Transition network.* ]{} We start the idea by constructing the ordinal pattern transition network for a single system, using the first system (Eqs. , $\omega_1 = 1.0, \Delta = 0, \kappa_{1,2} = 0$) as an example [@zhangSciRep2017]. Given time series $(x(t), y(t), z(t))$ in the corresponding three dimensional phase space ($n = 3$), the ordinal pattern transition network is reconstructed based on the signs of the increments of each each variablee $(\Delta x(t), \Delta y(t), \Delta z(t))$, where $\Delta x(t) = x(t+1) - x(t)$, $\Delta y(t) = y(t+1) - y(t)$, and $\Delta z(t) = z(t+1) - z(t)$. In particular, the definitions of the ordinal patterns based on the increment series capture the variations of the trajectory in its associated velocity space. The definition of patterns $\Pi(t) \in (\pi_1, \cdots, \pi_i), i = 1, \cdots, 8$ are enumerated in Tab. \[tab:3D\]. As time evolves, the transition behavior between patterns are illustrated in Fig. \[fig:rosCOPT\](a). The deterministic transitions are explained by the ordinal partitions of phase space by null-clines. A transition between two patterns means the trajectory crosses a null-cline, which leads to a local maximum or minimum. Because of the continuity of the system in phase space, we only observe the transition route $\pi_1 \to \pi_5 \to \pi_6 \to \pi_8 \to \pi_4 \to \pi_3 \to \pi_1$, yielding two missing patterns $\pi_2$ and $\pi_7$ as shown in Fig. \[fig:rosCOPT\](a). These details have been well illustrated in [@zhangSciRep2017].

     $\Pi$      $\pi_1$   $\pi_2$   $\pi_3$   $\pi_4$   $\pi_5$   $\pi_6$   $\pi_7$   $\pi_8$
  ------------ --------- --------- --------- --------- --------- --------- --------- ---------
   $\Delta x$    $ + $     $ + $     $  +$     $ +$     $  - $     $ - $     $ -$      $ -$
   $\Delta y$    $ + $     $ + $      $-$      $ -$      $ + $     $ + $     $ -$      $ -$
   $\Delta z$    $ + $     $ - $     $ +$      $ -$      $ + $     $ - $     $ +$      $ -$

  : Definitions of ordinal patterns of the three dimensional time series $(x(t), y(t), z(t))$, where $\Delta x = x(t+1) - x(t)$, $\Delta y = y(t+1) - y(t)$, and $\Delta z = z(t+1) - z(t)$ [@zhangSciRep2017]. \[tab:3D\]

![ \[fig:rosCOPT\]](rossler.pdf){width="\columnwidth"}

Generalizing the above ideas from a single system to two coupled systems $(\kappa_{1,2} \neq 0)$, we propose two different ways to construct ordinal pattern transition networks. Suppose we have time series $(x_{1}(t), y_{1}(t), z_{1}(t))$ of one system and $(x_{2}(t), y_{2}(t), z_{2}(t))$ of the other, we first compute the respective increments of the two systems as $(\Delta x_{1}(t), \Delta y_{1}(t), \Delta z_{1}(t))$ and $(\Delta x_{2}(t), \Delta y_{2}(t), \Delta z_{2}(t))$. Note that the increment series capture the dynamic properties of time series in the difference space and the signs of the each variable reflects either the increasing ($+$) or the decreasing ($-$) trend. Next, we compare the two systems in the following:

[ *Cross ordinal pattern transition network (COPT).* ]{} A COPT compares the relative speeds between two systems by the signs of $(\Delta x_1(t) - \Delta x_2(t)), (\Delta y_1(t) - \Delta y_2(t))$ and $(\Delta z_1(t) - \Delta z_2(t))$. The pattern definitions of a COPT are shown in Tab. \[tab:3DCOPT\]. An example of COPT is shown in Fig. \[fig:rosCOPT\](b). In this work, we only consider the signs of the same variable from two coupled systems. A further generalization of the pattern definitions of a COPT is to consider the signs of cross-variables, for instance, $(\Delta x_1(t) - \Delta y_2(t)) $ and $(\Delta x_1(t) - \Delta z_2(t))$, respectively.

             $\Pi$             $\pi_1$   $\pi_2$   $\pi_3$   $\pi_4$   $\pi_5$   $\pi_6$   $\pi_7$   $\pi_8$
  --------------------------- --------- --------- --------- --------- --------- --------- --------- ---------
   $\Delta x_1 - \Delta x_2$    $+ $      $+ $       $+$       $+$      $ - $     $ - $      $-$      $ - $
   $\Delta y_1 - \Delta y_2$    $ + $     $ + $     $ -$      $ -$      $ + $     $ + $     $ -$      $ -$
   $\Delta z_1 - \Delta z_2$    $ + $     $ - $     $ +$      $ -$      $ + $     $ - $      $+$      $ -$

  : Pattern definitions of a COPT. Note that $``+"$ means a positive value while $``-"$ is for a negative value. \[tab:3DCOPT\]

Considering the effects of the different magnitudes of the three variables, we also compute an [*alternative*]{} COPT by replacing $\Delta x_1(t) - \Delta x_2(t)$ by $\Delta x_1(t) / x_1(t) - \Delta x_2(t) / x_2(t)$, respectively, $\Delta y_1(t) - \Delta y_2(t)$ by $\Delta y_1(t) / y_1(t) - \Delta y_2(t) / y_2(t)$, and $\Delta z_1(t) - \Delta z_2(t)$ by $\Delta z_1(t) / z_1(t) - \Delta z_2(t) / z_2(t)$. An example of the alternative COPT is shown in Fig. \[fig:rosCOPT\](c). Comparing Fig. \[fig:rosCOPT\](b) to \[fig:rosCOPT\](c), the alternative COPT reflects better the non-coherent transitions between ordinal patterns since the coupling strength is in the non-synchronization regime ($\kappa_{1} = 0$ and $\kappa_{2} = 0.01$). In the following, we compute the alternative COPTs without distinguishing these two slightly different versions.

[*Joint ordinal pattern transition networks (JOPT).* ]{} A JOPT compares the relative speeds between two systems by the signs of $\Delta x_1(t) \cdot \Delta x_2(t), \Delta y_1(t) \cdot \Delta y_2(t)$ and $\Delta z_1(t) \cdot \Delta z_2(t)$ and the pattern definitions of a JOPT are summarized in Tab. \[tab:3DJOPT\]. An example of JOPT is shown in Fig. \[fig:rosCOPT\](d).

               $\Pi$               $\pi_1$   $\pi_2$   $\pi_3$   $\pi_4$   $\pi_5$   $\pi_6$   $\pi_7$   $\pi_8$
  ------------------------------- --------- --------- --------- --------- --------- --------- --------- ---------
   $\Delta x_1 \cdot \Delta x_2$    $+ $      $+ $       $+$       $+$      $ - $     $ - $      $-$      $ - $
   $\Delta y_1 \cdot \Delta y_2$    $ + $     $ + $     $ -$      $ -$      $ + $     $ + $     $ -$      $ -$
   $\Delta z_1 \cdot \Delta z_2$    $ + $     $ - $     $ +$      $ -$      $ + $     $ - $      $+$      $ -$

  : Pattern definitions of a JOPT. Note that $``+"$ means a positive value while $``-"$ is for a negative value. \[tab:3DJOPT\]

In contrast to cross ordinal patterns, we notice that the joint ordinal patterns represent whether the respective variables of two systems show the same trend of changes or not, regardless of the magnitudes of the respective variables.

A direct quantitative comparison between a COPT and a JOPT seems not possible since we have different definitions for patterns. Anyway, we present qualitative similarities when we show a COPT and a JOPT for two interacting stochastic processes in Fig. \[fig:stochasticCOPT\]. In particular, we consider two stationary processes $X_t$ and $Y_t$, each of which admits an autoregressive representation $$\begin{aligned}
X_t &= \sum_{j=1}^{p} a_{2j} X_{t-j} + \sum_{j=1}^{p} b_{2j} Y_{t-j} + \eta_{1t}, \\
Y_t &= \sum_{j=1}^{p} c_{2j} X_{t-j} + \sum_{j=1}^{p} d_{2j} Y_{t-j} + \eta_{2t},\end{aligned}$$ where the noise terms $\eta_{1t}, \eta_{2t}$ are uncorrelated, $p$ is the order, $a_{2j}, b_{2j}, c_{2j}$ and $d_{2j}$ are model coefficients. The terms of $b_{2j}$ and $c_{2j}$ reflect the interacting strength between $X_t$ and $Y_t$ and therefore if $X_t$ and $Y_t$ are independent, $b_{2j}$ and $c_{2j}$ are uniformly zero. Here we consider a simple case of $p = 2$, $a_{2j} = 0.5$, $b_{2j}=0.1$, $c_{2j} = 0.5$, and $d_{2j}=0$, when there is only unidirectional interactions from $Y_t$ to $X_t$.

In order to keep the same number of ordinal patterns as defined in Tabs. \[tab:3DCOPT\] and \[tab:3DJOPT\], we reconstruct a three dimensional phase space for $X_t$ (respectively $Y_t$) using time delay embedding techniques ($\tau = 1$). Figure \[fig:stochasticCOPT\]a shows the ordinal pattern transition network for a single stochastic process where we find rather random transition patterns in the resulting network. These random transitions between patterns have been observed in the COPT and JOPT (Figs. \[fig:stochasticCOPT\]b, c and d). When increasing the interacting strength term $b_{2j}$, one would expect some reduced level of random pattern transitions in the $X_t$ process because of the unidirectional interactions from the $Y_t$ process. The dependence on the interaction strength requires further investigations, which is beyond the topic of the current work.

![ \[fig:stochasticCOPT\]](grangerStoch.pdf){width="\columnwidth"}

Permutation entropy
-------------------

In the next step, we quantify the heterogeneous properties of the resulting networks. For a deterministic system, the frequencies of ordinal patterns are generally different from each other and the existence of forbidden patterns is simply a special case in this regard. Traditionally, permutation entropy $\mathcal{H}$ is introduced to characterize the inhomogeneous appearance of ordinal patterns as following $$\mathcal{H}_O = - \sum_{i=1}^{2^n} p(\pi_i) \log_2 p(\pi_i),$$ where the sum runs over all $D = 2^{n}$ permutations and $n$ is the dimension of one system and $p(\pi_i)$ is the probability of order pattern $\pi_i$. We use the observation frequency $\mathcal{F}(\pi_i)$ to estimate $p(\pi_i)$ and furthermore we use $\log_2$ and hence the units of $\mathcal{H}_O$ are bits. For a $n$-dimensional independent identical distributed stochastic process, one obtains the largest entropy $\mathcal{H}_O = n$ since each of $D = 2^{n}$ ordinal patterns is expected to have the same frequency.

The computation of $\mathcal{H}_O$ characterizes the different frequencies of order patterns and it has been well demonstrated that the transition behavior between ordinal patterns is not fully captured by $\mathcal{H}_O$ [@zhangSciRep2017]. To this end, we first indicate each directed link representing the order pattern transitions in the resulting network by its transition frequency $w_{ij} = p({\pi_i \to \pi_j})$, following the time iterations of the series. In order to emphasize the importance of non-self transitions between ordinal patterns, self-loops have been removed as suggested [@zhangSciRep2017]. Finally, we obtain a weighted directed network characterized by a weighted adjacency matrix $W = \{ w_{ij} \}, i, j \in [1, 2^{n}]$. The matrix $W$ fulfils the normalization $\sum_{i,j}^{2^n} w_{ij} = 1$. Here, based on $W$, the regularity of the order pattern transition properties is quantified by the Shannon entropy $\mathcal{H}_T$, which is $$\label{eq:HtS}
\mathcal{H}_T = - \sum_{i,j=1}^{2^{n}} w_{ij} \log_2 w_{ij},$$ where the sum runs over all possible $2^{2n}$ transitions. In a full analogy to $\mathcal{H}_O$, for a $n$-dimensional independent identical distributed stochastic process, one obtains the largest entropy $\mathcal{H}_T = 2 n$.

Both measures have been demonstrated to show the capabilities to capture the different bifurcation transition scenarios. In the examples of this work, we show that $\mathcal{H}_O$ and $\mathcal{H}_T$ are sensitive to PS.

Detecting transitions to PS {#secSync}
===========================

Preliminaries on PS
-------------------

In this section, traditional measures characterizing PS are briefly reviewed and more historical details can be found in [@PecoraChaos2015; @BoccalettiPR2002; @Pikovsky_Kurths_synchr]. PS is characterized by the phase locking $|m \phi_1 - n \phi_2 | < C$, where $m$, $n$ and $C$ are constants and $\phi_1(t) $ and $\phi_2(t) $ are phases of the two oscillators. In the case of two coupled Rössler oscillators, we easily compute phases by $\phi_1(t) = \arctan {y_1(t)}/{x_1(t)}$ and $\phi_2(t) = \arctan{y_2(t)}/{x_2(t)}$ since the systems are in phase coherent regimes. In addition, $m$ and $n$ are often chosen as $1$ because $1:1$ phase synchronization is more often observed. Equivalently, the phase locking is characterized as the average frequency mismatch drops to zero, namely, $\Delta \Omega = 0$, where $\Delta \Omega = \Omega_1 - \Omega_2$ and $\Omega_{1,2} = \frac{1}{2\pi} \left< \frac{d \phi_{1,2}(t)}{d t} \right>$, where $\left< \cdot \right>$ is an average over time $T$. PS has been related to the spectrum of Lyapunov exponents [@RosenblumPRL1997] and it is characterized by the transition of the Lyapunov exponent from zero to negative values [@Pikovsky_Kurths_synchr].

Unidirectional coupling schemes
-------------------------------

In the coupled Rössler systems (Eqs. \[eq\_roessler\]), PS can be achieved by applying the unidirectional coupling scheme $\kappa_1=0$ and relative small frequency mismatch $\Delta$ between $\omega_1$ and $\omega_2$ [@zhengPRE2000]. More specifically, we choose $\Delta = 0.02$ and $(\omega_1$, $\omega_2) = (0.98, 1.02)$. Note that in this model, generalized synchronization is obtained for relative large values of $\Delta = 0.2$ which leads to $(\omega_1, \omega_2) = (0.8, 1.2)$. The interrelationship between PS and generalized synchronization has been systematically investigated in [@zhengPRE2000].

![ \[fig:twoRosslerW1COPTJOPT\]](Single_Cross.pdf){width="\columnwidth"}

We first shown in Fig. \[fig:twoRosslerW1COPTJOPT\] that the variations of network measures of COPTs depending on the coupling parameters $\kappa_2$ ($\kappa_1 = 0$). When increasing the coupling $\kappa_2$, frequencies of ordinal patterns $\mathcal{F}(\pi_i)$ experience relative large variations (Figs. \[fig:twoRosslerW1COPTJOPT\](a)). In consequence, both $\mathcal{H}_{T}$ and $\mathcal{H}_{O}$ show fast decays to small values when the coupling threshold $\kappa_2$ passing the transition values (as highlighted by vertical dashed lines in Fig. \[fig:twoRosslerW1COPTJOPT\](b)). The transition point to PS has been validated by the average rotation frequencies of $\Omega_1$ and $\Omega_2$ since they are locked to the same value at $\kappa_2 = 0.07$. We choose three representative coupling values $\kappa_2=0.04, 0.08$, and $0.12$ to show the structural variations of the corresponding COPTs in Fig. \[fig:networkCOPJOP\](a-c). In the non-synchrony regime ($\kappa_2 = 0.04$), the transitions between ordinal patterns and their frequencies are rather random (Fig. \[fig:networkCOPJOP\](a)). As the coupling increases to just above the critical value $\kappa_2 = 0.08$, a dominant (more deterministic) transition route emerges (Fig. \[fig:networkCOPJOP\](b)). When PS is achieved for large $\kappa_2 = 0.12$, more missing patterns have been observed (Fig. \[fig:networkCOPJOP\](c)). We note that the different pattern definitions are used between the ordinal transition network for a single chaotic Rössler system as shown in Fig. \[fig:rosCOPT\](a) and the COPT for two coupled systems (Fig. \[fig:networkCOPJOP\](c)). The unidirectional coupling leads to the entrainment of the response to the drive system, which yields the same ordinal partitions for phase space as the drive. Due to the continuity of phase space trajectory, we observe a similar deterministic transition route between phase space partitions as for a single system.

Concerning network measures from JOPTs, we have obtained rather similar results as shown in Fig. \[fig:twoRosslerW2COPTJOPT\]. Three representative JOPTs on the route to PS are illustrated in Fig. \[fig:networkCOPJOP\](d-f). As the coupling strength increases, the random transitions between patterns become more deterministic yielding more missing patterns.

![ \[fig:twoRosslerW2COPTJOPT\]](Single_Joint.pdf){width="\columnwidth"}

![image](netGraphs.pdf){width="\textwidth"}

Based on the results of Figs. \[fig:twoRosslerW1COPTJOPT\], \[fig:twoRosslerW2COPTJOPT\] and \[fig:networkCOPJOP\], we conclude that $\mathcal{H}_O$ and $\mathcal{H}_T$ can effectively capture the transitions to PS. However, the detailed variations of particular patterns remain unclear during the transitions to PS. In other words, we have to look into the details of the appearance or disappearance of each ordinal patterns (as shown in Figs. \[fig:twoRosslerW1COPTJOPT\](a) and \[fig:twoRosslerW2COPTJOPT\](a)) because some of the patterns become suppressed, which provide further evidence in characterizing the synchronization dynamics.

It is easy to understand the occurrences of missing patterns on the transition routes to synchronization. As the coupling strength is increased during the synchronization transitions, the high dimensional coupled systems are constrained to a lower dimensional synchronization manifold which means that the determinism of the system is increased. In consequence, both the frequencies of order patterns and the transitions between different patterns become more inhomogeneous. Therefore, we compute the missing probability (frequency) of each pattern while varying the coupling strengths.

Bidirectional coupling schemes
------------------------------

In this section, we implement bidirectional coupling scheme when choosing $\kappa_1 = \kappa_2 = \kappa$. In addition, we show the transition to synchronization in a two dimensional parameter space of coupling strength $\kappa \in [0.0, 1.2]$ and natural frequency mismatch $\Delta \in [-0.04, 0.04]$. We choose this particular range of parameters because there are different routes of chaos-chaos, chaos-period-chaos transitions to synchronizations [@mamenPLA2004]. This space $(\Delta, \kappa)$ is further divided into $800 \times 1200$ grid points by equal step size. For each parameter combination, we integrate the ODEs with the same strategy as described in Sec. \[secRungeK\] to obtain time series for the coupled system. Then, we construct both COPTs and JOPTs.

First, we show the average frequency mismatches $\Delta \Omega$ in the parameter space $(\Delta, \kappa)$ and PS is characterized by the well-known Arnold tongue [@Pikovsky_Kurths_synchr], as shown in Fig. \[fig:twoparaCOPTOne\](a). PS are observed inside this Arnold tongue, while no PS outside this region. We further validate these results by computing the Lyapunov spectrum of the whole system based on Eqs. \[eq\_roessler\] (Fig. \[fig:twoparaCOPTOne\](b)). We note that it may not be possible to distinguish the tip of the Arnold tongue only by considering the sum of positive Lyapunov exponents because of the non-hyperbolic property of the coupled $6$-dimensional system. One interesting structure inside the Arnold tongue is the periodic region close to $0.03 < \kappa < 0.04$ (like two eyes), which has been reported in [@mamenPLA2004]. It is often claimed that all transitions between different types of synchronization are related to the changes in the Lyapunov spectrum [@RosenblumPRL1997]. In particular, for low values of the coupling strength, one has the following configuration for the coupled $6$-dimension system: $\{ \lambda_1 > 0, \lambda_2 > 0, \lambda_3 \sim 0, \lambda_4 \sim 0, \lambda_5 < 0, \lambda_6 < 0 \}$. Increasing the coupling strength, PS is achieved when $\lambda_4$ becomes negative. If the coupling strength is further increased, generalized synchronization is obtained as $\lambda_2 \sim 0$ and $\lambda_3 < 0$. However, for intermediate coupling strengths and relative large frequency mismatches, the spectrum of Lyapunov exponents has limited power in explaining the weak correlations outside the Arnold tongue, although the phases are not locked [@mamenPLA2004].

Figures \[fig:twoparaCOPTOne\](c-f) show the parameter space color coded by the entropy values $\mathcal{H}_O$ and $\mathcal{H}_T$, which are calculated from the respective COPTs and JOPTs. The outer boarders of the Arnold tongue have been successfully captured by $\mathcal{H}_{O}$ and $\mathcal{H}_{T}$, comparing to the average frequency mismatch plot (Fig. \[fig:twoparaCOPTOne\](a)). More importantly, all these entropy values are able to capture the tip of the Arnold tongue, however, the spectrum of Lyapunov exponents does not (Fig. \[fig:twoparaCOPTOne\](b)). Therefore, Figs. \[fig:twoparaCOPTOne\](c-f) provide important complementary information to the sum of the positive Lyapunov exponents.

Inside the Arnold tongue, however, these entropy values drop to small values, presenting more intricate gradient structures. Anyway, entropy values do not show convincingly the periodic region close to $0.03 < \kappa < 0.04$ except for the case of $\mathcal{H}_{T}$ computed from COPTs. This is because $\mathcal{H}_O$ and $\mathcal{H}_T$ are averaged over all patterns.

![image](double_HtHoLamH.pdf){width="\textwidth"}

In order to illustrate more clearly the gradient structures inside and on the edge of the Arnold tongue, we plot the missing (failure) probabilities of individual ordinal patterns in the space $(\Delta, \kappa)$. For one combination of parameter pair $(\Delta, \kappa)$, we run $n=200$ simulations of random initial conditions. The COPT and JOPT have been reconstructed for each realization of $N=800000$ time points and the frequency plots of ordinal patterns have been obtained. Furthermore, this ensemble of independent realizations results in a failure with probability $p_{\pi_i} = 1$ if pattern $\pi_i$ is missing and a success with probability $p_{\pi_i} = 0$ if $\pi_i$ is observed, respectively.

In the case of COPTs, Fig. \[fig:twoparaCOPTTwoC\] shows that $\pi_1, \pi_3, \pi_4, \pi_5, \pi_6$ and $\pi_8$ are observed patterns for the entire range of parameters we considered. In addition, we find that pattern $\pi_2$ is completely not observed in the Arnold tongue and the disappearance of $\pi_2$ determines sharply the outer boarder. In addition, pattern $\pi_7$ plays a key role in forming the periodic eye structures inside the Arnold tongue. In the case of JOPTs, Fig. \[fig:twoparaCOPTTwoJ\] shows that patterns $\pi_1, \pi_2, \pi_3$ and $\pi_5$ do appear in the majority area of the parameter space. Again, the disappearance of patterns $\pi_7$ and $\pi_8$ determine sharply the outer boarders of the Arnold tongue, while the inside structures are determined by $\pi_4$ and $\pi_6$.

![image](eightPatterns_doubleCross.pdf){width="\textwidth"}

![image](eightPatterns_doubleJoint.pdf){width="\textwidth"}

Conclusions {#secCons}
===========

In summary, we propose to construct cross and joint ordinal pattern transition network from multivariate time series, which has been particularly applied to analyze synchronization transitions. Note that a COPT and JOPT are two slightly different ways to construct networks from time series, providing complementary information. The ordinal patterns of a COPT are defined by considering the signs of the difference of $\Delta \vec{x}_1 - \Delta \vec{x}_2$ between two subsystems. In contrast, the ordinal patterns of a JOPT are defined by the signs of the product of $\Delta \vec{x}_1 \cdot \Delta \vec{x}_2$. It is certain that the amplitudes of oscillations of different variables influence directly the definition of a COPT. However amplitudes become not important for a JOPT because only the signs of the product are considered. In addition, it is straightforward to generalize the ideas of JOPTs from two to three (or even $n$) coupled subsystems with an extended number of pattern definitions. We plan to further generalize the ideas of JOPTs to multilayer or multiplex networks for time series analysis. However, it remains to be a big challenge for constructing a COPT for three coupled subsystems.

Based on the cross and joint ordinal pattern transition networks, we propose two entropy measures to characterize the resulting networks. Namely, $\mathcal{H}_O$ is calculated from the frequencies of each patterns and $\mathcal{H}_T$ is obtained from the transition frequencies between any pair of patterns. Our results show that both $\mathcal{H}_O$ and $\mathcal{H}_T$ track successfully the critical coupling threshold to phase synchronization. The applications of our method to generalization synchronization analysis is a more challenging task and will be a subject for future work [@PecoraChaos2015].

In the two parameter space of $(\Delta, \kappa)$, both entropy measures capture the tip of the Arnold tongue successfully, providing complementary information to the traditional measure of Lyapunov exponents. In order to show the intricate structures inside the Arnold tongue, our results suggest that we should study the missing probability of each pattern separately, instead of relying on the global measures of $\mathcal{H}_O$ and $\mathcal{H}_T$. This is because individual pattern shows different sensitivity to dynamic transitions. In particular, we have observed particular missing patterns which correspond to the outer boarders and the inner structures in the parameter space.

In addition to synchronization analysis based on real time series, it remains to be an interesting topic to identify the driver-response relationship, especially to identify indirect from direct coupling directions [@rosenblum_pre2001; @RomanoPRE2007; @Nawrath_prl_2010; @zouIJBC2011]. From the viewpoint of ordinal pattern perspective, it is possible to combine ordinal recurrence plots [@grothPRE2005] and cross and joint ordinal partition transition network approaches to tackle this problem. More importantly, we plan to address the statistical significance of the coupling directions for real time series data.

Acknowledgement
===============

This work is in part financially sponsored by Natural Science Foundation of Shanghai (Grants No. 17ZR1444800 and 18ZR1411800).

[10]{} url \#1[`#1`]{}urlprefixhref \#1\#2[\#2]{} \#1[\#1]{}

R. V. Donner, M. Small, J. F. Donges, N. Marwan, Y. Zou, R. Xiang, J. Kurths, Recurrence-based time series analysis by means of complex network methods, Int. J. Bifurcation Chaos 21 (4) (2011) 1019–1046. [](http://dx.doi.org/10.1142/S0218127411029021).

R. V. Donner, Y. Zou, J. F. Donges, N. Marwan, J. Kurths, Recurrence networks – [A]{} novel paradigm for nonlinear time series analysis, New J. Phys. 12 (3) (2010) 033025. [](http://dx.doi.org/10.1088/1367-2630/12/3/033025).

N. Marwan, J. F. Donges, Y. Zou, R. V. Donner, J. Kurths, Complex network approach for recurrence analysis of time series, Phys. Lett. A 373 (46) (2009) 4246–4254. [](http://dx.doi.org/10.1016/j.physleta.2009.09.042).

L. Lacasa, B. Luque, F. Ballesteros, J. Luque, J. C. Nuno, From time series to complex networks: The visibility graph, Proc. Natl. Acad. Sci. USA 105 (13) (2008) 4972–4975. [](http://dx.doi.org/10.1073/pnas.0709247105).

J. Zhang, M. Small, Complex network from pseudoperiodic time series: Topology versus dynamics, Phys. Rev. Lett. 96 (2006) 238701. [](http://dx.doi.org/10.1103/PhysRevLett.96.238701).

Y. Yang, H. Yang, Complex network-based time series analysis, Physica A 387 (5) (2008) 1381 – 1386. [](http://dx.doi.org/10.1016/j.physa.2007.10.055).

J. F. Donges, R. V. Donner, M. H. Trauth, N. Marwan, H. J. Schellnhuber, J. Kurths, Nonlinear detection of paleoclimate-variability transitions possibly related to human evolution, Proc. Natl. Acad. Sci. USA 108 (2011) 20422–20427. [](http://dx.doi.org/10.1073/pnas.1117052108).

Y. Zou, R. V. Donner, M. Wickramasinghe, I. Z. Kiss, M. Small, J. Kurths, Phase coherence and attractor geometry of chaotic electrochemical oscillators, Chaos 22 (3) (2012) 033130. [](http://dx.doi.org/10.1063/1.4747707).

Z.-K. Gao, W.-D. Dang, Y.-X. Yang, Q. Cai, Multiplex multivariate recurrence network from multi-channel signals for revealing oil-water spatial flow behavior, Chaos 27 (3) (2017) 035809. [](http://dx.doi.org/10.1063/1.4977950).

J. B. Elsner, T. H. Jagger, E. A. Fogarty, Visibility network of united states hurricanes, Geophys. Res. Lett. 36 (16) (2009) L16702. [](http://dx.doi.org/10.1029/2009GL039129).

Y. Zou, M. Small, Z. Liu, J. Kurths, Complex network approach to characterize the statistical features of the sunspot series, New J. Phys. 16 (1) (2014) 013051. [](http://dx.doi.org/10.1088/1367-2630/16/1/013051).

Y. Zou, R. Donner, N. Marwan, M. Small, J. Kurths, Long-term changes in the north-south asymmetry of solar activity: a nonlinear dynamics characterization using visibility graphs, Nonlin. Proc. Geophys. 21 (6) (2014) 1113–1126. [](http://dx.doi.org/10.5194/npg-21-1113-2014).

R. Zhang, Y. Zou, J. Zhou, Z.-K. Gao, S. Guan, Visibility graph analysis for re-sampled time series from auto-regressive stochastic processes, Commun. Nonlinear Sci. Numer. Simulat. 42 (2017) 396 – 403. [](http://dx.doi.org/10.1016/j.cnsns.2016.04.031).

Z. Czechowski, M. Lovallo, L. Telesca, Multifractal analysis of visibility graph-based ito-related connectivity time series, Chaos 26 (2) (2016) 023118. [](http://dx.doi.org/10.1063/1.4942582).

C. Zhang, Y. Chen, G. Hu, Network reconstructions with partially available data, Frontiers of Physics 12 (3) (2017) 128906. [](http://dx.doi.org/10.1007/s11467-017-0664-z).

Z.-Q. Jiang, Y.-H. Yang, G.-J. Wang, W.-X. Zhou, Joint multifractal analysis based on wavelet leaders, Frontiers of Physics 12 (6) (2017) 128907. [](http://dx.doi.org/10.1007/s11467-017-0674-x).

R. V. Donner, J. Heitzig, J. F. Donges, Y. Zou, N. Marwan, J. Kurths, The geometry of chaotic dynamics — a complex network perspective, Eur. Phys. J. B 84 (4) (2011) 653–672. [](http://dx.doi.org/10.1140/epjb/e2011-10899-1).

M. McCullough, M. Small, T. Stemler, H. H.-C. Iu, Time lagged ordinal partition networks for capturing dynamics of continuous dynamical systems, Chaos 25 (5) (2015) 053101. [](http://dx.doi.org/10.1063/1.4919075).

C. W. Kulp, J. M. Chobot, H. R. Freitas, G. D. Sprechini, Using ordinal partition transition networks to analyze ecg data, Chaos 26 (7) (2016) 073114. [](http://dx.doi.org/10.1063/1.4959537).

C. W. Kulp, J. M. Chobot, B. J. Niskala, C. J. Needhammer, Using forbidden ordinal patterns to detect determinism in irregularly sampled time series, Chaos 26 (2) (2016) 023107. [](http://dx.doi.org/10.1063/1.4941674).

M. McCullough, K. Sakellariou, T. Stemler, M. Small, Counting forbidden patterns in irregularly sampled time series. i. the effects of under-sampling, random depletion, and timing jitter, Chaos 26 (12) (2016) 123103. [](http://dx.doi.org/10.1063/1.4968551).

K. Sakellariou, M. McCullough, T. Stemler, M. Small, Counting forbidden patterns in irregularly sampled time series. ii. reliability in the presence of highly irregular sampling, Chaos 26 (12) (2016) 123104. [](http://dx.doi.org/10.1063/1.4970483).

C. [Bandt]{}, B. [Pompe]{}, Permutation entropy: A natural complexity measure for time series, Phys. Rev. Lett. 88 (17) (2002) 174102. [](http://dx.doi.org/10.1103/PhysRevLett.88.174102).

U. Parlitz, H. Suetani, S. Luther, Identification of equivalent dynamics using ordinal pattern distributions, Eur. Phys. J. S.T. 222 (2) (2013) 553–568. [](http://dx.doi.org/10.1140/epjst/e2013-01859-2).

J. M. Amig[ó]{}, K. Keller, V. A. Unakafova, Ordinal symbolic analysis and its application to biomedical recordings, Phil. Trans. R. Soc. A 373 (2034) (2014) 20140091. [](http://dx.doi.org/10.1098/rsta.2014.0091).

F. Takens, [Detecting strange attractors in turbulence]{}, in: D. Rand, L.-S. Young (Eds.), Dynamical Systems and Turbulence, Warwick 1980, Vol. 898 of Lecture Notes in Mathematics, Springer, New York, 1981, pp. 366–381. [](http://dx.doi.org/10.1007/BFb0091924).

H. Kantz, T. Schreiber, Nonlinear Time Series Analysis, 2nd Edition, Cambridge University Press, Cambridge, 2004.

J. M. Amigó, S. Zambrano, M. A. F. Sanjuán, True and false forbidden patterns in deterministic and random dynamics, Europhys. Lett. 79 (5) (2007) 50001. [](http://dx.doi.org/10.1209/0295-5075/79/50001).

J. M. Amigó, S. Zambrano, M. A. F. Sanjuán, Combinatorial detection of determinism in noisy time series, Europhys. Lett. 83 (6) (2008) 60005. [](http://dx.doi.org/10.1209/0295-5075/83/60005).

O. A. Rosso, L. C. Carpi, P. M. Saco, M. G. Ravetti, H. A. Larrondo, A. Plastino, The [A]{}mig[ó]{} paradigm of forbidden/missing patterns: a detailed analysis, Eur. Phys. J. B 85 (12) (2012) 419. [](http://dx.doi.org/10.1140/epjb/e2012-30307-8).

O. A. Rosso, L. C. Carpi, P. M. Saco, M. G. Ravetti, A. Plastino, H. A. Larrondo, Causality and the entropy-complexity plane: [R]{}obustness and missing ordinal patterns, Physica A 391 (1Ð2) (2012) 42 – 55. [](http://dx.doi.org/10.1016/j.physa.2011.07.030).

C. W. Kulp, L. Zunino, Discriminating chaotic and stochastic dynamics through the permutation spectrum test, Chaos 24 (3) (2014) 033116. [](http://dx.doi.org/10.1063/1.4891179).

A. Politi, Quantifying the dynamical complexity of chaotic time series, Phys. Rev. Lett. 118 (2017) 144101. [](http://dx.doi.org/10.1103/PhysRevLett.118.144101).

J. Zhang, J. Zhou, M. Tang, H. Guo, M. Small, Y. Zou, Constructing ordinal partition transition networks from multivariate time series, Scientific Reports 7 (2017) 7795. [](http://dx.doi.org/10.1038/s41598-017-08245-x).

A. Pikovsky, M. Rosenblum, J. Kurths, Synchronization – A Universal Concept in Nonlinear Sciences, Cambridge University Press, 2001.

G. V. Osipov, B. Hu, C. Zhou, M. V. Ivanchenko, J. Kurths, Three types of transitions to phase synchronization in coupled chaotic oscillators, Phys. Rev. Lett. 91 (2) (2003) 024101. [](http://dx.doi.org/10.1103/PhysRevLett.91.024101).

J. Zhang, Y.-Z. Yu, X.-G. Wang, Synchronization of coupled metronomes on two layers, Front. Phys. 12 (6) (2017) 120508. [](http://dx.doi.org/10.1007/s11467-017-0675-9).

H.-B. Chen, Y.-T. Sun, J. Gao, C. Xu, Z.-G. Zheng, Order parameter analysis of synchronization transitions on star networks, Front. Phys. 12 (6) (2017) 120504. [](http://dx.doi.org/10.1007/s11467-017-0651-4).

X. Huang, J. Gao, Y.-T. Sun, Z.-G. Zheng, C. Xu, Effects of frustration on explosive synchronization, Front. Phys. 11 (6) (2016) 110504. [](http://dx.doi.org/10.1007/s11467-016-0597-y).

L.-M. Ying, J. Zhou, M. Tang, S.-G. Guan, Y. Zou, Mean-field approximations of fixation time distributions of evolutionary game dynamics on graphs, Frontiers of Physics 13 (1) (2017) 130201. [](http://dx.doi.org/10.1007/s11467-017-0698-2).

Z. Zheng, G. Hu, Generalized synchronization versus phase synchronization, Phys. Rev. E 62 (2000) 7882–7885. [](http://dx.doi.org/10.1103/PhysRevE.62.7882).

M. C. Romano, M. Thiel, J. Kurths, W. von Bloh, Multivariate recurrence plots, Phys. Lett. A 330 (3) (2004) 214 – 223. [](http://dx.doi.org/10.1016/j.physleta.2004.07.066).

L. M. Pecora, T. L. Carroll, Synchronization of chaotic systems, Chaos 25 (9) (2015) 097611. [](http://dx.doi.org/10.1063/1.4917383).

S. Boccaletti, J. Kurths, G. Osipov, D. Valladares, C. Zhou, The synchronization of chaotic systems, Physics Reports 366 (1) (2002) 1 – 101. [](http://dx.doi.org/10.1016/S0370-1573(02)00137-0).

M. G. Rosenblum, A. S. Pikovsky, J. Kurths, From phase to lag synchronization in coupled chaotic oscillators, Phys. Rev. Lett. 78 (1997) 4193–4196. [](http://dx.doi.org/10.1103/PhysRevLett.78.4193).

M. G. Rosenblum, A. S. Pikovsky, Detecting direction of coupling in interacting oscillators, Phys. Rev. E 64 (4) (2001) 045202. [](http://dx.doi.org/10.1103/PhysRevE.64.045202).

M. C. Romano, M. Thiel, J. Kurths, C. Grebogi, Estimation of the direction of the coupling by conditional probabilities of recurrence, Phys. Rev. E 76 (2007) 036211. [](http://dx.doi.org/10.1103/PhysRevE.76.036211).

J. Nawrath, M. C. Romano, M. Thiel, I. Z. Kiss, M. Wickramasinghe, J. Timmer, J. Kurths, B. Schelter, Distinguishing direct from indirect interactions in oscillatory networks with multiple time scales, Phys. Rev. Lett. 104 (3) (2010) 038701. [](http://dx.doi.org/10.1103/PhysRevLett.104.038701).

Y. Zou, M. C. Romano, M. Thiel, N. Marwan, J. Kurths, Inferring indirect coupling by means of recurrences, Int. J. Bifurcation Chaos 21 (04) (2011) 1099–1111. [](http://dx.doi.org/10.1142/S0218127411029033).

A. Groth, Visualization of coupling in time series by order recurrence plots, Phys. Rev. E 72 (4) (2005) 046220. [](http://dx.doi.org/10.1103/PhysRevE.72.046220).
