## A Mechanism of Stochastic Synchronization in the Climate System: an Interpretation of the Boundary Current Synchronization as a Maxwell's Demon

Yuki Yasuda a and Tsubasa Kohyama b a Supercomputing Research Center, Institute of Integrated Research, Institute of Science Tokyo b Department of Information Sciences, Ochanomizu University

ABSTRACT: This study has applied information thermodynamics to a bivariate linear stochastic differential equation (SDE) that describes a synchronization phenomenon of sea surface temperatures (SSTs) between the Gulf Stream and the Kuroshio Current, which is referred to as the boundary current synchronization (BCS). Information thermodynamics divides the entire system fluctuating with stochastic noise into subsystems and describes the interactions between these subsystems from the perspective of information transfer. The SDE coefficients have been estimated through regression analysis using observational and numerical simulation data. In the absence of stochastic noise, the solution of the estimated SDE shows that the SSTs relax toward zero without oscillation. The estimated SDE can be interpreted as a Maxwell's demon system, with the Gulf Stream playing the role of the 'Particle' and the Kuroshio Current playing the role of the 'Demon.' The Gulf Stream forces the SST of the Kuroshio Current to be in phase. By contrast, the Kuroshio Current maintains the phase by interfering with the relaxation of the Gulf Stream SST. In the framework of Maxwell's demon, the Gulf Stream is measured by the Kuroshio Current, whereas the Kuroshio Current performs feedback control on the Gulf Stream. When the Gulf Stream and the Kuroshio Current are coupled in an appropriate parameter regime, synchronization is realized with atmospheric and oceanic fluctuations as the driving source. This new mechanism, 'stochastic synchronization,' suggests that such fluctuations can be converted into directional variations in a subsystem of the climate system through synchronization by utilizing information (i.e., information-to-energy conversion).

## 1. Introduction

Stochastic modeling is an effective approach for various climate phenomena (e.g., Dijkstra 2013; Franzke et al. 2015). This approach is based on the time scale separation between slow and fast variations (Hasselmann 1976). Typically, slow processes include decadal climate variations, while fast processes include weather disturbances. The scale separation suggests that the correlation between slow and fast processes is weak; thus, the fast processes can effectively be represented by stochastic noise. Stochastic modeling describes the integral (or aggregated) responses of slow climate variables driven by stochastic noise. A major approach in stochastic modeling is the use of idealized dynamical systems with few variables (e.g., Dijkstra 2013). Although these models cannot describe phenomena in detail, they are expected to capture the essence of phenomena (e.g., Benzi 2010). In this study, we investigate the coupling of two western boundary currents using a stochastic dynamical system and propose a new mechanism for synchronization, which may be applicable to other climate phenomena.

Western boundary currents (WBCs) are strong upper flows located on the western side of ocean basins, playing an important role in the transport of heat in the climate system (Kwon et al. 2010). Two major examples of WBCs are the Gulf Stream in the North Atlantic and the Kuroshio Current in the North Pacific. These WBCs are primarily driven by atmospheric wind stress (Stommel 1948) and modulated through interactions with the atmosphere (Qiu et al. 2014; Ma et al. 2016). The heat transport by the WBCs influences the entire atmosphere, including storm tracks and annular modes (Hoskins and Valdes 1990; Minobe et al. 2008; Ogawa et al. 2012; Omrani et al. 2019). Furthermore, the WBCs sometimes contribute to explosive cyclogenesis, making them also important for daily weather (Sanders 1986; Kuwano-Yoshida and Minobe 2017).

The WBCs, namely the Gulf Stream and the Kuroshio Current, exhibit interannual to decadal variability (Qiu and Chen 2005; Kelly et al. 2010). Decadal variability due to the interaction of either the Gulf Stream or the Kuroshio Current with the atmosphere has been studied (Latif and Barnett 1994; Gallego and Cessi 2000; Hogg et al. 2006; Kwon et al. 2010; Qiu et al. 2014). Recently, the importance of the coexistence of both the Gulf Stream and the Kuroshio Current has been pointed out for ocean heat content (Kelly and Dong 2004) and atmospheric annular modes (Omrani et al. 2019). However, the roles of both ocean currents in the climate system are not yet fully understood.

Kohyama et al. (2021) discovered the boundary current synchronization (BCS) by using observational data and numerical experiments with coupled atmosphere-ocean models. The sea surface temperatures (SSTs) of the Gulf Stream and the Kuroshio Current show positive correlations on interannual to decadal time scales. An increase (decrease) in the SSTs is accompanied by a northward (southward) shift of the WBCs as well as a northward (southward) shift of the atmospheric baroclinic jet. The BCS is considered to arise from the interaction between the Gulf Stream and the Kuroshio Current via the atmosphere (i.e., the inter-basin coupling) rather than a passive response to atmospheric forcings. In fact, the BCS was not reproduced in low-resolution numerical experiments that could not resolve oceanic mesoscale eddies. These results suggest that the BCS occurs through the transfer of information on WBC variability through the atmosphere. However, the roles of the Gulf Stream and the Kuroshio Current in the BCS have not been well investigated, and further research is required.

Over the past decade, the fusion of information theory and non-equilibrium physics has led to the development of a new theoretical framework called 'information thermodynamics' (Parrondo et al. 2015; Peliti and Pigolotti 2021; Shiraishi 2023). Information thermodynamics divides nonlinear systems fluctuating with stochastic noise into subsystems and describes the interactions between these subsystems from the perspective of information transfer. Information transfer can be expressed as information flow (Allahverdyan et al. 2009) and is incorporated into the second law of thermodynamics (Horowitz and Esposito 2014). Applying this extended second law of thermodynamics (i.e., the second law of information thermodynamics), we can elucidate the dynamics of subsystems and the asymmetry between them. This new theoretical framework is applicable to stochastic systems that include deterministic forcing and dissipation. These points suggest that information thermodynamics is suitable for studying the climate system, which is a nonlinear system with multiple degrees of freedom driven by stochastic noise and deterministic forcings. However, since information thermodynamics is an emerging field, it has scarcely been used in atmospheric and oceanic sciences.

This study applies information thermodynamics to a dynamical system for the BCS and proposes that the coupled system of the Gulf Stream and the Kuroshio Current can be interpreted as a Maxwell's demon. Section 2 provides a review of Maxwell's demon. Section 3 introduces the bivariate linear dynamical system that describes the time evolution of the regional-mean SSTs of

the Gulf Stream and the Kuroshio Current. Section 4 applies information thermodynamics to this system and discusses the roles of the Gulf Stream and the Kuroshio Current. Section 5 presents the conclusions.

## 2. Review of Maxwell's demon

The resolution of Maxwell's demon is a landmark in information thermodynamics (Parrondo et al. 2015). Originally, Maxwell argued that if an intelligent being (a demon) had information about the velocities of particles, it could separate fast, hot particles from slow, cold particles without exerting any work on the system, thus appearing to violate the second law of thermodynamics. Szilard then developed a simple model, the Szilard engine, which captures the essential point of Maxwell's demon: information can be converted to work. Nowadays, this information-to-energy conversion has been realized experimentally (Toyabe et al. 2010; Koski et al. 2014). This type of system is called an information engine or loosely Maxwell's demon (e.g., Horowitz and Esposito 2014; Ito and Sagawa 2015). Here, we refer to the entire system of information conversion as the Maxwell's demon system and review the Szilard engine as such a system, which consists of a demon (i.e., a controller) and a subsystem operated by the demon. See further details in reviews and textbooks (Sagawa and Ueda 2013a; Parrondo et al. 2015; Peliti and Pigolotti 2021; Shiraishi 2023).

## a. Setup of the Szilard engine

Let us consider an isothermal cycle with a single particle. Figure 1 is a schematic of this experiment, where the memory (i.e., the demon) in the upper half is not considered for the moment. First, a partition is inserted in the middle of the box, separating the box into halves (Fig. 1a). The particle moves at random, so the probability of the particle being on the left (or right) side equals 1 2. Second, the position of the particle is measured (Fig. 1b). Third, based on the measurement / outcome, when the particle is on the left (resp. right), the partition is moved quasi-statically to the right (resp. left). (Fig. 1c). Such a measurement-based operation is called feedback control or simply, control. The energy source for this expansion is the surrounding thermal fluctuations, or equivalently, heat transferred through the box. Therefore, the Szilard engine appears to contradict

Fig. 1. Schematic of the Szilard engine (e.g., Parrondo et al. 2015). A box contains a single particle and is surrounded by an isothermal environment, where heat is transferred through the walls of the box. (a) A partition is inserted in the middle of the box. (b) The demon measures the position of the particle and records the result in its memory. (c) Based on the measurement outcome, the demon performs feedback control. Specifically, when the particle is on the left (resp. right), the demon moves the partition to the right (resp. left). The particle system is isothermally expanded and returns to the initial state. This process extracts a positive work from the isothermal environment. (d) The demon's memory is initialized, requiring external work (i.e., a negative work).

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01133v2-figure-34c81068-44ea-4597-b521-821dc21c1bc3.jpg)

the second law of thermodynamics and to be a perpetual motion machine of the second kind that extracts work via isothermal cycles.

To resolve the paradox, it is necessary to consider the exchange of information quantities and the existence of memory (Sagawa and Ueda 2013a). The memory, also called the demon, is a subsystem that performs measurement and feedback control. The other subsystem is the particle within the box and is referred to as the particle. We denote the states of the demon and the particle by 𝑋 and 𝑌 and their realizations by lowercase letters 𝑥 and 𝑦 , respectively. The variables 𝑋 and 𝑌 are discrete, reflecting the particle position being either on the left or right. We first derive the second law of information thermodynamics and then consider the Szilard engine again.

## b. The second law of information thermodynamics

We define Shannon entropy and mutual information (e.g., Cover and Thomas 2005) to derive the second law of information thermodynamics.

$$S ( X ) \coloneqq \sum _ { x } p ( x ) [ - \ln p ( x ) ],$$

$$S ( Y ) \coloneqq \sum _ { y } p ( y ) [ - \ln p ( y ) ],$$

$$M ( X, Y ) \coloneqq \sum _ { x, y } ^ { \cdot } p ( x, y ) \ln \left [ \frac { p ( x, y ) } { p ( x ) p ( y ) } \right ].$$

Shannon entropies 𝑆 ( 𝑋 ) and 𝑆 𝑌 ( ) are indicators of uncertainty about the subsystems 𝑋 and 𝑌 , respectively. Mutual information 𝑀 𝑋,𝑌 ( ) is a non-negative indicator of the correlation between 𝑋 and 𝑌 , where 𝑀 𝑋,𝑌 ( ) is zero only when 𝑋 and 𝑌 are statistically independent.

Mutual information characterizes the non-additivity of Shannon entropy (e.g., Cover and Thomas 2005). For the joint Shannon entropy 𝑆 ( 𝑋,𝑌 ) , the following holds:

$$S ( X, Y ) \coloneqq \sum _ { x, y } p ( x, y ) [ - \ln p ( x, y ) ],$$

$$= S ( X ) + S ( Y ) - M ( X, Y ).$$

Addingtheuncertainties for the subsystems, 𝑆 ( 𝑋 ) + 𝑆 𝑌 ( ) , leads to an overestimation and is generally not equal to the uncertainty for the entire system 𝑆 ( 𝑋,𝑌 ) . It is necessary to subtract 𝑀 𝑋,𝑌 ( ) , which is non-negative and represents the magnitude of correlation between 𝑋 and 𝑌 .

The second law of information thermodynamics is derived using the non-additivity of Shannon entropy (Sagawa and Ueda 2013b). We consider the evolution from time 𝑡 to 𝑡 + Δ 𝑡 and assume the second law of thermodynamics for the entire system (Seifert 2005, 2012):

$$\Delta S ( X _ { t }, Y _ { t } ) + \frac { \Delta Q _ { t } } { T } \geq 0,$$

where changes are denoted by Δ , and time dependence is denoted by a subscript 𝑡 . For example, Δ 𝑆 ( 𝑋𝑡 ,𝑌 𝑡 ) : = 𝑆 ( 𝑋𝑡 + Δ 𝑡 ,𝑌 𝑡 + Δ 𝑡 ) -𝑆 ( 𝑋𝑡 ,𝑌 𝑡 ) . The symbol Δ 𝑄𝑡 represents the amount of heat dissipated from the joint system ( 𝑋,𝑌 ) to the surrounding environment, which corresponds to the entropy

change of the environment. Equation (6) expresses that the total entropy change, including the environment, is non-negative.

The second law of information thermodynamics for the subsystem 𝑋 is derived from Eq. (6) by assuming that only 𝑋 changes in time. Using the entropy decomposition [Eq. (5)], the following is derived:

$$\Delta S ( X _ { t }, Y _ { t } ) + \frac { \Delta Q _ { t } } { T } = \Delta S ( X _ { t } ) + \underbrace { \Delta S ( Y _ { t } ) } _ { = 0 } - \Delta M ( X _ { t }, Y _ { t } ) + \frac { \Delta Q _ { t } } { T } \geq 0,$$

where Δ 𝑆 𝑌 ( ) = 0 because 𝑌𝑡 does not change in time. This assumption also indicates that the entropy change of the environment Δ 𝑄𝑡 / 𝑇 is due to the change in 𝑋𝑡 . We finally obtain the following inequality:

$$\Delta \sigma \left ( X _ { t } \right ) \geq \Delta M \left ( X _ { t }, Y _ { t } \right ),$$

where

$$\Delta \sigma \left ( X _ { t } \right ) \coloneqq \Delta S \left ( X _ { t } \right ) + \frac { \Delta Q _ { t } \left ( X _ { t } \right ) } { T },$$

$$\Delta M \left ( X _ { t }, Y _ { t } \right ) = M ( X _ { t + \Delta t }, Y _ { t } ) - M ( X _ { t }, Y _ { t } ).$$

The total entropy change by 𝑋𝑡 is denoted as Δ 𝜎 𝑋𝑡 ( ) and is lower bounded by the change in mutual information Δ 𝑀 𝑋𝑡,𝑌𝑡 ( ) . Since 𝑌𝑡 does not change, Δ 𝑀 𝑋𝑡,𝑌𝑡 ( ) is attributed only to the change in 𝑋𝑡 . The second law of information thermodynamics [Eq. (8)] indicates that Δ 𝜎 𝑋𝑡 ( ) can be negative when Δ 𝑀 𝑋𝑡,𝑌𝑡 ( ) &lt; 0. Thus, the entropy of the subsystem 𝑋 can be decreased by utilizing the change in mutual information without heat transfer. This conclusion is consistent with thermodynamics because this discussion started from the second law of thermodynamics for the entire system [Eq. (6)]. Note that Δ 𝜎 𝑌𝑡 ( ) ≥ 𝑀 𝑋𝑡,𝑌𝑡 ( + Δ 𝑡 ) -𝑀 𝑋𝑡,𝑌𝑡 ( ) is derived when only 𝑌 changes in time.

## c. Explanation of the Szilard engine

The second law of information thermodynamics is used to explain the Szilard engine. Time 𝑡 is not explicitly stated if there is no confusion. In the initial state, the states of the demon 𝑋 and the particle 𝑌 are uncorrelated, that is, 𝑀 𝑋,𝑌 ( ) is zero. First, the demon inserts a partition (Fig. 1a).

This causes the particle 𝑌 to be either on the right or the left. The mutual information 𝑀 𝑋,𝑌 ( ) remains zero because the demon's internal state (i.e., the memory 𝑋 ) has not changed.

Next, the demon measures the position of the particle and changes its memory 𝑋 to record the measurement outcome, which makes 𝑋 either right or left (Fig. 1b). The states of the demon and the particle become correlated, and the mutual information increases, Δ 𝑀 𝑋,𝑌 ( ) &gt; 0. This increase in correlation is due to the change in the demon 𝑋 . According to the second law of information thermodynamics, the increase in correlation is accompanied by a positive entropy change, Δ 𝜎 𝑋 ( ) ≥ Δ 𝑀 𝑋,𝑌 ( ) ( &gt; 0) [Eq. (8)].

Based on the memory value 𝑋 , the demon quasi-statically moves the partition to extract work (Fig. 1c). Regardless of the memory value, the particle 𝑌 returns to the initial state after this feedback control, and the correlation 𝑀 𝑋,𝑌 ( ) becomes zero again. This decrease in correlation is due to the change in the particle 𝑌 . According to the second law of information thermodynamics, the decrease in correlation can be accompanied by a negative entropy change, Δ 𝜎 𝑌 ( ) ≥ Δ 𝑀 𝑋,𝑌 ( ) ( &lt; 0). This negative Δ 𝜎 𝑌 ( ) corresponds to the positive work extracted by the isothermal expansion.

Finally, work is performed from the outside to initialize the demon's memory 𝑋 (Fig. 1d). There is a trade-off for the memory between initialization and measurement, and it is more appropriate to consider both together (Sagawa and Ueda 2009, 2013b; Sagawa 2019). The positive work extracted by the isothermal expansion and the negative work required for the initialization (and measurement) cancel each other out. No net positive work can be extracted using the Szilard engine. Therefore, the Szilard engine is not a perpetual motion machine of the second kind.

Information thermodynamics clarifies that the positive work extracted during the control process is due to the change in mutual information; that is, the total entropy of a subsystem can decrease by consuming mutual information. In this way, information thermodynamics allows for the analysis of subsystems from the perspective of information.

## d. Autonomous Maxwell's demon systems

Maxwell's demon systems are also realized in an autonomous manner (Horowitz and Esposito 2014; Loos and Klapp 2020). The time evolution of an autonomous system can be expressed by stochastic differential equations without external interference, such as human manipulation. In such a system, one cannot isolate individual steps, such as measurement and control, since all

steps can occur simultaneously and continuously over time. This is a major difference between the Szilard engine and autonomous processes. Thus, the distinction between the 'Demon' and the 'Particle' appears unclear.

The determination of which subsystem is the 'Demon' or the 'Particle' is based on the direction of information flow (Allahverdyan et al. 2009; Horowitz and Esposito 2014; Loos and Klapp 2020).

$$i _ { X \leftarrow Y } \coloneqq \frac { M \left ( X _ { t + d t }, Y _ { t } \right ) - M \left ( X _ { t }, Y _ { t } \right ) } { d t },$$

$$i _ { Y \leftarrow X } \coloneqq \frac { M \left ( X _ { t }, Y _ { t + d t } \overset { \ddots } { \right ) } - M \left ( X _ { t }, Y _ { t } \right ) } { d t },$$

$$\frac { d } { d t } M \left ( X _ { t }, Y _ { t } \right ) = \dot { I } _ { X \leftarrow Y } + \dot { I } _ { Y \leftarrow X }.$$

The information flow / 𝐼 𝑋 ← 𝑌 from 𝑌 to 𝑋 is the change in mutual information 𝑀 𝑋𝑡,𝑌𝑡 ( ) due to the change in 𝑋 , whereas / 𝐼 𝑌 ← 𝑋 is due to the change in 𝑌 . The sum of both information flows is equal to the total rate of change of mutual information [Eq. (13)]. Information flow is regarded as a quantity that represents the flow of information between subsystems.

For autonomous systems, the second law of information thermodynamics becomes as follows (Horowitz and Esposito 2014; Loos and Klapp 2020):

$$\dot { \sigma } ( X ) \geq \dot { I } _ { X \leftarrow Y },$$

$$\dot { \sigma } ( Y ) \geq \dot { I } _ { Y \leftarrow X }.$$

The rate of change of total entropy is lower bounded by the information flow. Hereafter, the rate of change of total entropy is referred to as the production rate of total entropy or simply the entropy production rate. Equations (14) and (15) are derived for general autonomous systems (Appendix A).

We finally summarize the distinction between a 'Demon' and a 'Particle.' As in Section 2c, measurement is a change that increases correlation due to an evolution of the 'Demon', which is accompanied by a positive entropy production. On the other hand, control is a change that decreases correlation due to an evolution of the 'Particle', which is accompanied by a negative entropy production. Therefore, in an autonomous system, the 'Demon' (resp. 'Particle') is the subsystem showing a positive (resp. negative) information flow and entropy production rate

(Horowitz and Esposito 2014; Loos and Klapp 2020). In our study, a system consisting of a 'Demon' and a 'Particle' has been referred to as the Maxwell's demon system. Maxwell's demon systems are no longer thought experiments but have been experimentally realized in both nonautonomous and autonomous settings (Toyabe et al. 2010; Koski et al. 2015). We propose that such systems can also be realized in the atmosphere and ocean, that is, the coupled system of the Gulf Stream and the Kuroshio Current can be interpreted as a Maxwell's demon system.

## 3. Bivariate linear dynamical system for BCS

## a. Governing equations

This section introduces a bivariate dynamical system describing BCS, which is hereafter referred to as the dynamical system. This system consists of the regional-mean SST anomalies of the Gulf Stream and the Kuroshio Current, 𝑇𝐺 and 𝑇𝐾 , respectively. The subscripts 𝐺 and 𝐾 denote the Gulf Stream and the Kuroshio Current, respectively. Following Kohyama et al. (2021), 𝑇𝐺 is the spatially averaged SST in the Gulf Stream region (35 N-45 N, 80 W-50 W), and ◦ ◦ ◦ ◦ 𝑇𝐾 is the spatially averaged SST in the Kuroshio Current region (35 N-45 N, 140 E-170 E). Both ◦ ◦ ◦ ◦ 𝑇𝐺 and 𝑇𝐾 are standardized (i.e., non-dimensionalized) after removing the climatology and the linear trend.

The variables 𝑇𝐺 and 𝑇𝐾 are assumed to obey the following equations:

$$\frac { d T _ { G } } { d t } = - r _ { G } T _ { G } + c _ { G \leftarrow K } T _ { K } + \xi _ { G },$$

$$\frac { d T _ { K } } { d t } = - r _ { K } T _ { K } + c _ { K \leftarrow G } T _ { G } + \xi _ { K },$$

where 𝑟 𝐺 , 𝑐 𝐺 ← 𝐾 , 𝑟 𝐾 , and 𝑐 𝐾 ← 𝐺 are real constants. The unit of time is set to 1 month. The mutually independent white Gaussian noises 𝜉 𝐺 and 𝜉 𝐾 satisfy

$$\langle \xi _ { G } ( t ) \xi _ { G } ( s ) \rangle = 2 D _ { G } \delta ( t - s ),$$

$$\langle \xi _ { K } ( t ) \xi _ { K } ( s ) \rangle = 2 D _ { K } \delta ( t - s ),$$

$$\langle \xi _ { K } ( t ) \xi _ { G } ( s ) \rangle = 0,$$

where 𝐷𝐺 and 𝐷𝐾 are positive constants, 𝛿 𝑡 ( - ) 𝑠 is the Dirac delta function, 𝑡 -𝑠 is the time difference, and ⟨·⟩ represents the ensemble average over the noise realizations.

This type of a system that consists only of oceanic variables (e.g., SSTs) can be derived from a more complicated setting that include both oceanic and atmospheric variables. Indeed, Gallego and Cessi (2001) started with an atmosphere-ocean coupling model and then derived idealized equations that consist only of Atlantic and Pacific variables by representing the atmospheric effects in terms of the oceanic variables. Our dynamical system is interpreted as a simplified model derived through these procedures. The cross terms, 𝑐 𝐺 ← 𝐾𝑇𝐾 and 𝑐 𝐾 ← 𝐺𝑇𝐺 , may include interbasin couplings via the atmosphere.

## b. Regression analysis

The coefficients of the dynamical system were estimated by the regression analysis using two sets of SST time series, both of which are used in Kohyama et al. (2021). One is the observational data, Optimum Interpolation SST (OISST, Reynolds et al. 2007). The other is the numerical output from a coupled atmosphere-ocean model, GFDL-CM4C192 (Held et al. 2019; Eyring et al. 2016). Monthly mean values were used for both datasets. The OISST data range from December 1981 to September 2018, whereas the GFDL-CM4C192 data range from January 1950 to December 2014. Using the longer-period GFDL-CM4C192 data, we expect to enhance the reliability of the analysis results. All of the following analyses were performed independently for the OISST and GFDL-CM4C192 data.

Figure 2 shows the OISST data, GFDL-CM4C192 data, and their regression results, where a firstorder bivariate autoregressive (AR1) model was employed (Brockwell and Davis 1991; Mudelsee 2014). In estimating the AR1 model, the SST at the current time step was regressed onto the SST at the previous time step, and the regression coefficients were obtained. This regression method is known as the ordinary least squares (OLS) and is standard in time series analysis (Hamilton 1994; Mudelsee 2014). The AR1 model reproduces the SST variations well (Fig. 2), suggesting the validity of the AR1 model. The estimated AR1 model can be uniquely converted to the continuous-time dynamical system [Eqs. (16) and (17)]. Specifically, under the assumption of statistical stationarity, we require that the means and variances of the variations in 𝑇𝐺 and 𝑇𝐾 are the same between the AR1 model and the dynamical system (see Subsection a in Appendix B for

Fig. 2. Time series of the standardized SST anomalies for the Gulf Stream (blue) and the Kuroshio Current (orange) from the OISST and GFDL-CM4C192 data. The dashed gray lines represent the regression results for the corresponding time series with the first-order autoregressive (AR1) model.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01133v2-figure-ff3fc2c8-3620-40af-83dd-ca6751a3356b.jpg)

details). In other words, this conversion requires that the statistical properties of the AR1 model match those of the time-integrated values over each discrete time step (i.e., one month) using the dynamical system.

We used the moving block bootstrap method (Mudelsee 2014) to quantify the uncertainty in the estimated coefficients of the dynamical system. This bootstrap method is a non-parametric method that does not assume normality. This method first splits the original time series into blocks. These blocks are resampled with replacement, and the resampled blocks are connected to generate a new time series of the same length as the original time series, which is called the resampled time series. The coefficients of the dynamical system are estimated for each resampled time series. We can compute confidence intervals of the coefficients for the dynamical system by repeating this procedure and creating histograms of the coefficients. The block length was quantitatively determined from the auto-correlation coefficients (Michael Sherman and Speed 1998; Mudelsee 2014): the block lengths were 18 months for the OISST data and 30 months for the GFDLCM4C192 data. The shorter block length for the OISST data may be attributed to its shorter time series length. The number of resamples was set to 2000 following Mudelsee (2014).

Table 1 shows the estimation results of the coefficients of the dynamical system, where the brackets show the 95% confidence intervals by the bootstrap method. The relaxation coefficients

Table 1. Estimated coefficients of the dynamical system [Eqs. (16) and (17)] from the OISST and GFDLCM4C192 data. The values represent the estimates obtained from the regression analysis, and the brackets show the 95% confidence intervals obtained from the moving block bootstrap method.

|         | OISST                   | GFDL-CM4C192           |
|---------|-------------------------|------------------------|
| 𝑟 𝐺     | 0.343 [0.327, 0.555]    | 0.299 [0.289, 0.429]   |
| 𝑟 𝐾     | 0.325 [0.301, 0.538]    | 0.206 [0.196, 0.321]   |
| 𝑐 𝐺 ← 𝐾 | 0.0280 [-0.0634, 0.117] | 0.105 [0.0433, 0.171]  |
| 𝑐 𝐾 ← 𝐺 | 0.110 [0.00387, 0.240]  | 0.0756 [0.0249, 0.144] |
| 𝐷 𝐺     | 0.343 [0.324, 0.545]    | 0.253 [0.252, 0.382]   |
| 𝐷 𝐾     | 0.308 [0.278, 0.519]    | 0.176 [0.173, 0.283]   |

𝑟 𝐺 and 𝑟 𝐾 are positive, and their inverses give the relaxation time scales. The time scale is about 3 months for the Gulf Stream and the Kuroshio Current. The estimates of the interaction coefficients 𝑐 𝐺 ← 𝐾 and 𝑐 𝐾 ← 𝐺 are positive. The confidence intervals suggest that 𝑐 𝐺 ← 𝐾 based on the OISST data may be negative. The impact of this negative value is discussed in Appendix C. The variances of the noise 𝐷𝐺 and 𝐷𝐾 are of similar magnitude. We confirmed that the uncertainty quantification is not strongly dependent on the choice of method, as similar results were obtained using a parametric approach based on normality (not shown).

In the absence of noise, the solution of the dynamical system decays to zero without oscillation. According to linear dynamical system theory (e.g., Strogatz 2018), solutions of a deterministic linear system are categorized using the trace and determinant of its coefficient matrix. In our case, the coefficient matrices from all 2000 resampled time series exhibit decaying solutions, indicating that the dynamical system is stable. Moreover, most of these estimated coefficients are categorized into the decaying solution without oscillation when 𝑐 𝐺 ← 𝐾 &gt; 0 and 𝑐 𝐾 ← 𝐺 &gt; 0: all estimates from the GFDL-CM4C192 data belong to this category, whereas 1624 out of 2000 estimates from the OISST data belong to this category. For the OISST, the remaining 376 estimates show decaying oscillators due to the negative 𝑐 𝐺 ← 𝐾 (Table 1), the case of which is discussed in Appendix C. The influences of noise are analyzed using information thermodynamics in Section 4.

## c. Validation of the assumptions for the dynamical system

We evaluate the validity of the estimated dynamical system. The following assumptions have been made regarding the dynamical system: (i) statistical stationarity and (ii) zero auto- and cross-correlations for the noises 𝜉 𝐺 and 𝜉 𝐾 . We verify these two assumptions.

First, regarding statistical stationarity, we used the Dickey-Fuller test (Hamilton 1994), which tests the null hypothesis of non-stationarity against the alternative hypothesis of a stationary AR1 process. For the OISST and GFDL-CM4C192 data, the 𝑝 -value was at most around 10 -13 . This result supports the validity of regarding each time series as stationary.

Second, we confirm the validity of no correlations, namely Eqs. (18)-(20). The first two equations indicate that the noises do not depend on the past, resulting in zero auto-correlations. Equation (20) indicates that the noises 𝜉 𝐺 and 𝜉 𝐾 are independent of each other, resulting in zero cross-correlations. In regression analysis, the residuals between the estimates and the target variables are considered as noise. Thus, we can verify the assumptions by examining the auto- and cross-correlation coefficients of these residuals.

Figure 3 shows the results of the correlation analysis for these residuals. The gray areas represent the 95% confidence intervals (Brockwell and Davis 1991). For most lags, the auto-correlation coefficients fall within the confidence intervals, suggesting that the auto-correlations can be considered zero. These results support the assumption that the noises do not strongly depend on the past. The cross-correlation coefficients also take values close to zero and fall within the confidence intervals for most time lags. This result supports the assumption that the noises of the Gulf Stream and the Kuroshio Current are statistically independent. Kohyama et al. (2021) showed that the BCS is not a passive phenomenon driven by the atmospheric baroclinic jet. Our results suggest that the SSTs of both currents are not driven by common noise, which is consistent with the result of Kohyama et al. (2021). These noises are attributed to atmospheric and oceanic fluctuations on time scales of a few months or less, likely encompassing the spatially localized effects of the atmospheric baroclinic jet and oceanic mesoscale eddies in the Gulf Stream and Kuroshio Current regions.

We finally clarify the relationship between noise and temperature. In non-equilibrium physics, the variance of stochastic noise is generally proportional to the temperature of an isothermal environment; this noise arises from collisions of small molecules with the larger target particle (e.g.,

Fig. 3. Auto- and cross-correlation coefficients of the residuals between the AR1-model estimates and the target time series for the OISST and GFDL-CM4C192 data. The gray shading indicates the 95% confidence intervals.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01133v2-figure-cabce1eb-d090-4bdc-82e1-c90344b2f46e.jpg)

Sekimoto 2010). In our case, since the noises 𝜉 𝐺 and 𝜉 𝐾 represent atmospheric and oceanic fluctuations, their variances 𝐷𝐺 and 𝐷𝐾 are interpreted as the strengths of these small-scale disturbances. Thus, isothermal environments correspond to the disturbance fields, and 'temperature' in non-equilibrium physics corresponds to 𝐷𝐺 and 𝐷𝐾 . These effective temperatures should be distinguished from the SST anomalies 𝑇𝐺 and 𝑇𝐾 at longer time scales.

## d. Reproduction of the BCS by the dynamical system

Finally, we confirm that the estimated dynamical system reproduces the BCS between 𝑇𝐺 and 𝑇𝐾 . To show the BCS, we examine the lag correlation between 𝑇𝐺 and 𝑇𝐾 following Kohyama et al. (2021), who directly obtained the lag correlation coefficients from the time series. In the present study, we first estimated the coefficients of the dynamical system from the same time series. Then, using these coefficients, we calculated the analytical solutions of the lag correlation coefficients using the theoretical formulas (Subsection b in Appendix B).

Figure 4 shows the resultant lag correlations. Here, the 95% confidence interval for each lag correlation coefficient was estimated using the 2000 sets of resampled time series obtained by the moving block bootstrap method. The solid lines indicate that the correlation coefficient takes significant values larger than zero, which supports the existence of the synchronicity. The longer interval of significant correlations for the GFDL-CM4C192 data is likely due to its longer data range compared to that of the OISST data. The Gulf Stream slightly leads the Kuroshio Current

Fig. 4. Lag correlation coefficients between the Gulf Stream and Kuroshio Current SSTs for the OISST and GFDL-CM4C192 data. The solid lines indicate significant correlations at the 95% confidence level, while the dashed lines represent insignificant correlations.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01133v2-figure-88d607d7-7637-4a98-9803-6589e84e4e2d.jpg)

by 1 and 0.2 months for the OISST and GFDL-CM4C192 data, respectively. The significance of these small leads is difficult to assess given the short data periods. The obtained lag correlation coefficients are very close to the results of Kohyama et al. (2021) (see their Fig. 1D). Therefore, we can conclude that the dynamical system serves as an appropriate model for describing BCS.

## 4. Analysis of the dynamical system based on information thermodynamics

## a. Application of the second law of information thermodynamics

We apply the second law of information thermodynamics (Horowitz and Esposito 2014; Loos and Klapp 2020) to the Gulf Stream 𝑇𝐺 and the Kuroshio Current 𝑇𝐾 to clarify their asymmetric roles in the BCS. We assume a statistical steady state in which the probability distribution of 𝑇𝐺 and 𝑇𝐾 does not vary over time. The time-series analysis supports the validity of the steady-state assumption (Section 3c).

The second law of information thermodynamics gives the following inequalities and equality regarding the information flows and entropy production rates for the Gulf Stream and the Kuroshio

Current (Horowitz and Esposito 2014; Loos and Klapp 2020):

$$\dot { \sigma } _ { G } \geq i _ { G \leftarrow K },$$

$$\dot { \sigma } _ { K } \geq i _ { K \leftarrow G },$$

$$0 = \dot { I } _ { G \leftarrow K } + \dot { I } _ { K \leftarrow G },$$

where / 𝜎𝐺 and / 𝜎𝐾 represent the production rates of total entropy for the Gulf Stream 𝑇𝐺 and the Kuroshio Current 𝑇𝐾 , respectively, and / 𝐼 𝐺 ← 𝐾 and / 𝐼 𝐾 ← 𝐺 represent the information flow to the Gulf Stream and to the Kuroshio Current, respectively. Equations (21), (22), and (23) correspond to Eqs. (14), (15), and (13), respectively. The sum of the two information flows is zero in Eq. (23) because the time derivative of mutual information is zero due to the steady-state assumption (Appendix A). Moreover, under this assumption, the production rates of total entropy correspond to Δ 𝑄 in Eq. (9) because the change of Shannon entropies is zero. The interpretation of Δ 𝑄 for the dynamical system is discussed in Subsection e in Appendix A. Adding Eqs. (21) and (22) yields / 𝜎𝐺 + / 𝜎𝐾 ≥ 0, i.e., the total entropy production for the entire system is non-negative, indicating the consistency with the second law of thermodynamics. Considering each subsystem, / 𝜎𝐺 and / 𝜎𝐾 are lower bounded by / 𝐼 𝐺 ← 𝐾 and / 𝐼 𝐾 ← 𝐺 , respectively, one of which can be negative. By using information thermodynamics, we investigate the asymmetry between the dynamics of 𝑇𝐺 and 𝑇𝐾 from the perspective of entropy production and information transfer.

For the dynamical system [Eqs. (16) and (17)], the analytical solutions for the entropy production rates and information flows are obtained (Loos and Klapp 2020):

$$\dot { \sigma } _ { G } = \frac { c _ { G \leftarrow K } } { D _ { G } } \left \langle T _ { K } \frac { \text{d} T _ { G } } { \text{d} t } \right \rangle,$$

$$\dot { \sigma } _ { K } = \frac { c _ { K \leftarrow G } } { D _ { K } } \, \dot { \left < T _ { G } } \frac { \mathrm d T _ { K } } { \mathrm d t } \right >,$$

$$\dot { I } _ { G \leftarrow K } = \frac { \Sigma _ { G K } } { | \mathbf S | } \left \langle T _ { K } \frac { d T _ { G } } { d t } \right \rangle,$$

$$i _ { K \leftarrow G } = \frac { \overset { \cdot } { \Sigma _ { G K } } } { | \mathbf S | } \left \langle T _ { G } \frac { d T _ { K } } { d t } \right \rangle,$$

where we introduce the covariance matrix Σ between 𝑇𝐺 and 𝑇𝐾 . The off-diagonal component of Σ is represented by Σ 𝐺𝐾 , and the determinant of Σ is represented by | Σ | . Equations (24)-(27) indicate

that the entropy production rates and information flows are linked to the correlations between the SSTs and their rates of change. Thus, the temporal characteristics of the dynamical system can be examined using information thermodynamics. Details on the derivation of Eqs. (24)-(27) are given in Subsection c in Appendix B.

In the steady state, the entropy production rate and information flow of the Kuroshio Current are related to those of the Gulf Stream because of the following relationship (Subsection c in Appendix B):

$$\left \langle T _ { K } \frac { d T _ { G } } { d t } \right \rangle + \left \langle T _ { G } \frac { d T _ { K } } { d t } \right \rangle = 0.$$

Applying Eq. (28) into Eqs. (24)-(27), we obtain the following relationships:

$$\dot { \sigma } _ { K } = - \frac { c _ { K \leftarrow G } D _ { G } } { c _ { G \leftarrow K } D _ { K } } \dot { \sigma } _ { G },$$

$$\dot { I } _ { K \leftarrow G } = - \dot { I } _ { G \leftarrow K }.$$

Equation (30) can also be derived from the conservation of the information flows [Eq. (23)]. We analyze the dynamical system using the theoretical formulas (24)-(30).

## b. Analysis of regime diagrams for the BCS

Weshow that the dynamical system can be regarded as a Maxwell's demon system by examining the dependence of information thermodynamics quantities on the interaction coefficients 𝑐 𝐾 ← 𝐺 and 𝑐 𝐺 ← 𝐾 shown in Figs. 5 and 6. The parameters 𝑐 𝐾 ← 𝐺 and 𝑐 𝐺 ← 𝐾 are the most important because they characterize the interactions between the Gulf Stream and the Kuroshio Current. In Figs. 5 and 6, the white dots represent the estimated values of ( 𝑐 𝐾 ← 𝐺,𝑐𝐺 ← ) 𝐾 listed in Table 1. The error bars indicate the standard deviations obtained by the bootstrap method.

The error bars for the estimated values suggest that the true value is likely to lie in the region where 𝑐 𝐾 ← 𝐺 &gt; 0 and 𝑐 𝐺 ← 𝐾 &gt; 0 (Figs. 5 and 6), though the results based on the OISST data suggest that the true value may possibly lie in the region where 𝑐 𝐺 ← 𝐾 &lt; 0. Since the period of the OISST data is shorter, the uncertainty in the estimation is larger compared to the results of the GFDL-CM4C192 data. Appendix C briefly discusses the case where 𝑐 𝐺 ← 𝐾 &lt; 0. Hereafter, we focus on the region where 𝑐 𝐺 ← 𝐾 &gt; 0.

## OISST

Fig. 5. Dependence of information thermodynamic quantities on the interaction coefficients 𝑐 𝐾 ← 𝐺 and 𝑐 𝐺 ← 𝐾 for the OISST data: (a) / 𝜎 𝐺 , (b) / 𝐼 𝐺 ← 𝐾 , (c) √︃ ⟨ 𝑇 2 𝐺 ⟩ -√︃ ⟨ 𝑇 2 𝐺 ⟩| 𝑐 𝐺 ← 𝐾 = 0 , (d) / 𝜎 𝐾 , (e) / 𝐼 𝐾 ← 𝐺 , and (f) √︃ ⟨ 𝑇 2 𝐾 ⟩ -√︃ ⟨ 𝑇 2 𝐾 ⟩| 𝑐 𝐾 ← 𝐺 = 0 . Sub-figures (c) and (f) show the differences in the standard deviations of the SSTs from those without the interactions, namely, when 𝑐 𝐺 ← 𝐾 = 0 and 𝑐 𝐾 ← 𝐺 = 0, respectively. The white dots represent the estimated values of ( 𝑐 𝐾 ← 𝐺 , 𝑐 𝐺 ← 𝐾 ) , whose values are listed in Table 1. The error bars indicate the standard deviations obtained from the bootstrap method. The gray dashed diagonal lines show the symmetric condition 𝑐 𝐺 ← 𝐾 𝐷 𝐾 = 𝑐 𝐾 ← 𝐺 𝐷 𝐺 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01133v2-figure-b4d1d102-d91a-4371-99de-0007d57a0473.jpg)

The Gulf Stream and the Kuroshio Current are symmetric on the gray dashed diagonal lines in Figs. 5 and 6, where 𝑐 𝐺 ← 𝐾 𝐷 𝐾 = 𝑐 𝐾 ← 𝐺 𝐷 𝐺 (Appendix D). On these lines, the interchange of 𝑇𝐺 and 𝑇𝐾 does not affect the entropy production rates and information flows because these quantities are zero (Loos and Klapp 2020). Below these lines, / 𝜎𝐺 and / 𝐼 𝐺 ← 𝐾 are negative for the Gulf Stream, whereas / 𝜎𝐾 and / 𝐼 𝐾 ← 𝐺 are positive for the Kuroshio Current. Above these lines, all the signs of these quantities are reversed. Thus, the roles of the Gulf Stream and the Kuroshio Current are reversed across the diagonal lines in Figs. 5 and 6.

## GFDL-CM4C192

cG

cG

Fig. 6. Same as Fig. 5 but for the GFDL-CM4C192 data.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01133v2-figure-da8f0be7-d3c8-49c9-aa5d-c82d8d9fc864.jpg)

The Gulf Stream and the Kuroshio Current are likely to be asymmetric due to various factors, such as the land-sea distribution. This consideration suggests that the true value is not strictly on the gray dashed diagonal lines in Figs. 5 and 6, although the estimated values from the GFDL-CM4C192 data are quite close to these diagonal lines (Fig. 6). This result is probably due to the fact that small-scale disturbances regarded as noise are not sufficiently reproduced in the GFDL-CM4C192 model. Indeed, the noise variances 𝐷𝐺 and 𝐷𝐾 are smaller than those estimated from the OISST data (Table 1), and these parameters affect the symmetric condition 𝑐 𝐺 ← 𝐾 𝐷 𝐾 = 𝑐 𝐾 ← 𝐺 𝐷 𝐺 . Hereafter, we consider that the true value exists on the lower side of the diagonal lines, where both estimated values from the OISST and GFDL-CM4C192 data belong.

Figures 5 and 6 indicate that the dynamical system is a Maxwell's demon system, where the Gulf Stream is a 'Particle' and the Kuroshio Current is a 'Demon.' The entropy production rate and information flow of the Gulf Stream, / 𝜎𝐺 and / 𝐼 𝐺 ← 𝐾 , are negative (Figs. 5a, 5b, 6a, and 6b), whereas those of the Kuroshio Current, / 𝜎𝐾 and / 𝐼 𝐾 ← 𝐺 , are positive (Figs. 5d, 5e, 6d, and 6e). According

to the discussions on Maxwell's demon systems (Sections 2c and 2d), the Kuroshio Current as the 'Demon' tends to establish the correlation with the Gulf Stream SST by measuring it, which is reflected in the positive information flow ( / 𝐼 𝐾 ← 𝐺 &gt; 0). Furthermore, the Kuroshio Current utilizes its own SST as a memory and exhibits a positive entropy production rate ( / 𝜎𝐾 &gt; 0), which is attributed to the memory update process associated with the measurement. On the other hand, the Gulf Stream as the 'Particle' tends to decrease the correlation due to the feedback control by the Kuroshio Current, which is reflected in the negative information flow ( / 𝐼 𝐺 ← 𝐾 &lt; 0). In addition, the Gulf Stream SST varies by rectifying noise and shows a negative entropy production rate ( / 𝜎𝐺 &lt; 0). The term 'rectify' means to extract coherent directional output from random noise fluctuations.

We further discuss the rectification of noise. Figures 5c, 5f, 6c, and 6f show that the SST magnitudes √︃ ⟨ 𝑇 𝐺 2 ⟩ and √︃ ⟨ 𝑇 𝐾 2 ⟩ are increased compared to the case without interactions. As discussed in Section 3b, without noise, the dynamical system would decay to zero without oscillation. This result suggests that the driving source is noise. Moreover, in the absence of interactions, namely 𝑐 𝐺 ← 𝐾 = 𝑐 𝐾 ← 𝐺 = 0, the SSTs of the Gulf Stream and the Kuroshio Current would evolve independently, and BCS would not occur. Compared to the SST magnitudes in this case, they are increased due to the interactions. Within the framework of the dynamical system, the presence of both interactions and noises is essential for BCS.

We first investigate the reason for the increase in 𝑇𝐺 . The entropy production rate / 𝜎𝐺 is given as follows, which is negative (Figs. 5a and 6a):

$$0 > \dot { \sigma } _ { G } = \frac { c _ { G \leftarrow K } } { D _ { G } } \left \langle T _ { K } \frac { \mathrm d T _ { G } } { \mathrm d t } \right \rangle = \frac { ( c _ { G \leftarrow K } ) ^ { 2 } } { D _ { G } } \langle T _ { K } T _ { K } \rangle - \frac { r _ { G } c _ { G \leftarrow K } } { D _ { G } } \langle T _ { K } T _ { G } \rangle,$$

where all coefficients are positive. The rightmost side is obtained by substituting Eq. (16) into (24). The relative magnitudes of the two terms on the rightmost side characterize the behavior of the Gulf Stream. When the first term dominates, indicating a stronger influence from the Kuroshio Current, the Gulf Stream shows increased fluctuations and a larger entropy production rate / 𝜎𝐺 . Conversely, when the second term dominates, the relaxation is more effective, reducing the stochasticity of the Gulf Stream (i.e., decreasing / 𝜎𝐺 ). Clearly, the negative / 𝜎𝐺 indicates that the second term is dominant.

The negative / 𝜎𝐺 suggests that the action from the Kuroshio Current on the Gulf Stream, 𝑐 𝐺 ← 𝐾𝑇𝐾 , does not directly contribute to the increase in 𝑇𝐺 . Indeed, this action 𝑐 𝐺 ← 𝐾𝑇𝐾 is negatively

Fig. 7. Schematic of the increase in the Gulf Stream SST due to noise and feedback. Without the action from the Kuroshio Current, namely no feedback (top), the SST increases due to noise and quickly returns to zero through relaxation. With feedback from the Kuroshio Current (bottom), relaxation is interfered with, allowing the possibility of further increase by noise before returning to zero. The Kuroshio Current measures the Gulf Stream SST and inserts a 'partition' at an appropriate location.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01133v2-figure-9a364bde-2dda-4cf5-96d5-1709c6408466.jpg)

correlated with the tendency of the Gulf Stream SST, 𝑐 𝐺 ← ⟨ 𝐾 𝑇𝐾 d 𝑇𝐺 / d 𝑡 ⟩ = 𝐷𝐺𝜎𝐺 &lt; / 0. This result is due to the relaxation. Without relaxation, the governing equation (16) would be d 𝑇𝐺 / d 𝑡 = 𝑐 𝐺 ← 𝐾𝑇𝐾 , and only the first term in Eq. (31) would remain, leading to a positive / 𝜎𝐺 . The mechanism for the increase in 𝑇𝐺 is then described as follows. For example, when 𝑇𝐺 is positive, it tends to decrease to zero due to relaxation, but this decrease is mitigated by 𝑐 𝐺 ← 𝐾𝑇𝐾 . In other words, the action from the Kuroshio Current suppresses the relaxation of the Gulf Stream. Consequently, the remaining term 𝜉 𝐺 in Eq. (16), namely the noise, becomes more effective and increase 𝑇𝐺 . This suppression of relaxation to zero corresponds to the feedback control that inserts a partition (Fig. 7).

We next investigate the reason for the increase in 𝑇𝐾 . The expression for the entropy production rate / 𝜎𝐾 is given as follows, which is positive (Figs. 5d and 6d):

$$0 < \dot { \sigma } _ { K } = \frac { c _ { K \leftarrow G } } { D _ { K } } \left \langle T _ { G } \frac { \mathrm d T _ { K } } { \mathrm d t } \right \rangle = \frac { \left ( c _ { K \leftarrow G } \right ) ^ { 2 } } { D _ { K } } \langle T _ { G } T _ { G } \rangle - \frac { r _ { K } c _ { K \leftarrow G } } { D _ { K } } \langle T _ { G } T _ { K } \rangle,$$

where all coefficients are positive. The rightmost side is obtained by substituting Eq. (17) into (25). Equation (32) is interpreted in the same manner as Eq. (31). When the first term dominates, the Kuroshio Current shows strong fluctuations and a larger entropy production rate / 𝜎𝐾 . Conversely,

when the second term dominates, the stochasticity of the Kuroshio Current is reduced and / 𝜎𝐾 decreases. The result of / 𝜎𝐾 &gt; 0 indicates the dominance of the first term.

The positive / 𝜎𝐾 suggests that the action from the Gulf Stream on the Kuroshio Current, 𝑐 𝐾 ← 𝐺𝑇𝐺 , directly contributes to the increase in 𝑇𝐾 . Similar to the Gulf Stream, there is an effect of suppressing relaxation; however, its contribution is relatively small and / 𝜎𝐾 is positive. The amplification of 𝑇𝐾 occurs by following the changes in 𝑇𝐺 due to the action 𝑐 𝐾 ← 𝐺𝑇𝐺 from the Gulf Stream. According to the framework of Maxwell's demon systems, the Kuroshio Current (i.e., the 'Demon') uses its SST as memory and the SST fluctuations are increased by updating the memory value 𝑇𝐾 based on the measurement outcome of the Gulf Stream SST.

These analyses of 𝑇𝐺 and 𝑇𝐾 suggest that information-to-energy conversion (Toyabe et al. 2010) occurs in a subsystem (i.e., the Gulf Stream), even though the entire system is dissipative, as in other Maxwell demon systems. Generally, stochastic noise fluctuates a system back and forth, so the net energy or work is difficult to extract. However, using the information from measurements, fluctuations can be converted into directional movements in a subsystem (Toyabe et al. 2010), corresponding to the negative / 𝜎𝐺 or increased 𝑇𝐺 amplitude in our system. Such conversion cannot occur for the entire dissipative system (i.e., no net positive conversion), as reflected in the positive value of / 𝜎𝐺 + / 𝜎𝐾 .

Synchronicity is essential for this conversion. For instance, in Fig. 7, the partition needs to be inserted just behind the 'Particle.' This insertion is based on the memory value of the 'Demon,' as in the Szilard engine (Fig. 1), where the memory value becomes equal to the particle position after the measurement. Although it is difficult to directly compare our dynamical system to the Szilard engine because our system is autonomous (Section 2d), the memory value 𝑋 and particle position 𝑌 correspond to 𝑇𝐾 and 𝑇𝐺 , respectively. Both subsystems 𝑇𝐾 and 𝑇𝐺 are influenced by the small-scale disturbance fields (Section 3c), corresponding to isothermal environments in the Szilard engine. The 'Demon' ( 𝑇𝐾 ) follows a change in the 'Particle' ( 𝑇𝐺 ), thereby enabling feedback control to convert fluctuations. The 𝑇𝐾 amplitude increases due to this following (i.e., its continual updates based on 𝑇𝐺 ), although 𝑇𝐾 is dissipative (i.e., / 𝜎𝐾 &gt; 0). In the Szilard engine, this dissipation corresponds to the work required to initialize the memory.

The dependence of the lag correlation is consistent with the description of the Kuroshio Current following the Gulf Stream. Figure 8 shows the dependence of the maximum lag correlation

Fig. 8. Maximum lag correlation coefficients and the corresponding lags between the Gulf Stream and Kuroshio Current SSTs as functions of the interaction coefficients 𝑐 𝐾 ← 𝐺 and 𝑐 𝐺 ← 𝐾 . The results are shown for the OISST and GFDL-CM4C192 data. The gray dashed diagonal lines show the symmetric condition 𝑐 𝐺 ← 𝐾 𝐷 𝐾 = 𝑐 𝐾 ← 𝐺 𝐷 𝐺 . The label 'GS leads' means that the Gulf Stream leads, whereas the label 'KC leads' means that the Kuroshio Current leads.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01133v2-figure-d828e8dd-ae98-46f0-9b26-a4d657d0226c.jpg)

coefficient and the lag at that time over the same parameter space as in Figs. 5 and 6. Note that these lag correlation coefficients are theoretically computed as in Fig. 4. As confirmed in Section 3d, the lag correlation is positive in Fig. 8, suggesting that the BCS occurs. Above the gray dashed diagonal lines in Fig. 8, the SST of the Kuroshio Current leads by 1-2 months, while below the lines, the SST of the Gulf Stream leads. Across the diagonal lines, the roles of the 'Particle' and the 'Demon' are reversed because of the reversal of all signs of the entropy production rates and information flows. In the region below the diagonal lines, which we have been considering so far, the Gulf Stream (i.e., the 'Particle') leads. The Kuroshio Current (i.e., the 'Demon') measures the Gulf Stream SST and follows the changes in 𝑇𝐺 .

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01133v2-figure-1276c8b5-16a1-4b6f-9c47-087614e537ca.jpg)

/gid00020/gid00020/gid00021

Fig. 9. Schematic of the interpretation of the BCS as a Maxwell's demon system (i.e., the stochastic synchronization). The Kuroshio Current plays the role of the 'Demon' and the Gulf Stream plays the role of the 'Particle.' See Section 4c for details.

## c. Interpretation of BCS as a Maxwell's demon system: stochastic synchronization

From the above results, we obtain an interpretation of the BCS as a Maxwell's demon system (Fig. 9), which implies the asymmetric roles of the Gulf Stream and the Kuroshio Current. First, the Gulf Stream forces the SST of the Kuroshio Current to be in phase. This effect of making the Kuroshio Current follow is represented by the term 𝑐 𝐾 ← 𝐺𝑇𝐺 in the dynamical system, Eq. (17). In the framework of Maxwell's demon, the Gulf Stream is interpreted as being measured by the Kuroshio Current. By contrast, the Kuroshio Current locks the phase of the Gulf Stream SST by interfering with its relaxation toward the climatology. This effect of interference with the relaxation is represented by the term 𝑐 𝐺 ← 𝐾𝑇𝐾 in the dynamical system, Eq. (16). In the framework of Maxwell's demon, the Kuroshio Current is interpreted as performing feedback control on the Gulf Stream. On one hand, for the dynamical system without noise, the SSTs would relax to zero without oscillating. On the other hand, without interactions between the Gulf Stream and the Kuroshio Current, their SSTs would fluctuate independently. When both currents are coupled in an appropriate parameter regime, synchronization is realized with atmospheric and oceanic noise as the driving source. The physical origins of this driving source can be attributed to the effects of atmospheric jet streams and oceanic mesoscale eddies on time scales of a few months or less.

In this regard, Yamagami et al. (2024) recently examined the asymmetric roles between the Gulf Stream and Kuroshio Current using a coupled atmosphere-ocean general circulation model through 'pacemaker' experiments. In these experiments, SST variations only in a western boundary current region are strongly relaxed toward (i.e., keep 'pace' with) a control run, so that the responses of the

other climate systems to the western boundary current variability are isolated. They showed that, when the Gulf Stream acts as a pacemaker, BCS is reproduced. The heat release from the Gulf Stream modulates the atmospheric annular mode, exciting oceanic Rossby waves in the Pacific and influencing the SST in the Kuroshio Current region. This mechanism serves as a plausible physical process that realizes the role of the 'Particle' in the context of the Maxwell's demon system, i.e., the Gulf Stream forces the SST of the Kuroshio Current to be in phase. Conversely, they also showed that, when the Kuroshio Current acts as a pacemaker, the BCS is not reproduced. This result is also consistent with our interpretation that, unlike the Gulf Stream, the Kuroshio Current does not force the SST of the Gulf Stream by itself. Rather, by responding to the Gulf Stream's forcing, the Kuroshio Current returns a feedback to interfere with the Gulf Stream SST relaxation toward its climatology. The absence of BCS in the Kuroshio pacemaker experiment can be understood that, because its SST varies independently from the state of the Gulf Stream, the Kuroshio Current cannot reproduce an optimal feedback in an appropriate timing. Our future work is to reveal what kind of physical processes regarding the Kuroshio Current realizes the role of the 'Demon' in the real world.

Our proposed mechanism can be realized in other stochastic dissipative systems with positive feedback. The dynamical system investigated here, Eqs. (16) and (17), is known as the Langevin system, one of the simplest and most common models with dissipation and stochastic forcing (e.g., Sekimoto 2010). Indeed, Langevin systems have been employed as simple models for climate phenomena (e.g., Dijkstra 2013). The key ingredient in our study is the positive feedback; that is, both interaction coefficients 𝑐 𝐾 ← 𝐺 and 𝑐 𝐺 ← 𝐾 are positive. Due to this positive feedback, atmospheric or oceanic fluctuations are converted into directional variations in a subsystem via synchronization (i.e., information-to-energy conversion; Toyabe et al. (2010)), although such conversion cannot be realized for the entire dissipative system. The wide applicability of Langevin systems suggests that the proposed mechanism of synchronization (Fig. 9), 'stochastic synchronization,' may describe other climate phenomena. Such applications remain an important topic for future work.

## 5. Conclusions

This study has introduced a bivariate linear dynamical system that describes the BCS and analyzed it using a theory of information thermodynamics (Horowitz and Esposito 2014; Loos and Klapp

2020). This system can be interpreted as a Maxwell's demon system, with the Gulf Stream playing the role of the 'Particle' and the Kuroshio Current playing the role of the 'Demon.' Information thermodynamics clarifies the asymmetric roles of the Gulf Stream and the Kuroshio Current from the perspective of entropy production and information transfer (Section 4c). Our results suggest that random noise fluctuations can be autonomously converted into coherent, directional variations (i.e., information-to-energy conversion), through synchronization between the Gulf Stream and the Kuroshio Current. The proposed mechanism, named 'stochastic synchronization,' may describe other climate phenomena because our analysis was based on a widely applicable Langevin system. For future research, we propose applying information thermodynamics to more realistic models that describe the BCS (e.g., Gallego and Cessi 2001; Kohyama et al. 2021). While the current study demonstrates that the Gulf Stream can be interpreted as the 'Particle' and the Kuroshio Current as the 'Demon,' what determines these roles remains unclear, and these roles might be reversed under certain conditions, such as due to temporal changes in interaction coefficients. Moreover, it is unclear how the interactions, such as the interference of the relaxation, are realized in the real climate system. These points will be revealed by applying the theory to realistic models.

Information thermodynamics is applicable to other climate phenomena, such as the El Ni˜ noSouthern Oscillation, and may provide new physical insights. However, the direct application may not be straightforward because information thermodynamics requires the estimation of probability, and such estimation is generally difficult for systems with large degrees of freedom. Recently developed data-driven methods, such as dimensionality reduction (Reddy et al. 2020; Hou and Behdinan 2022), may be effective for the probability estimation and expand the applicability of information thermodynamics.

Acknowledgments. During the preparation of this work, the authors used Claude 3 only for English editing. After using this service, the authors reviewed and edited the content as needed. The authors take full responsibility for the content of the publication. The second author is supported by Japan Society for the Promotion of Science (JSPS) Kakenhi (22H04487, 23H01241, and 23K13169), and the MEXT program for the advanced studies of climate change projection (SENTAN) Grant Number JPMXD0722680395. The authors thank Masaru Inatsu for initiating our connection, leading to this collaboration.

Data availability statement. The data and source code that support the findings of this study are preserved at the Zenodo repository ( https://doi.org/10.5281/zenodo.13085327 ) and developed openly at the GitHub repository ( https://github.com/YukiYasuda2718/bcs\_ maxwell\_demon ).

## APPENDIX A

## Review of information thermodynamics

We re-derive the second law of information thermodynamics and the conservation law of information flows, both of which are used in this study. The derivation is mainly based on Seifert (2012), Horowitz and Esposito (2014), and Loos and Klapp (2020). We show the derivation in detail because information thermodynamics is rarely applied to atmospheric and oceanic sciences, and comprehensive reviews for continuous-state systems are not readily available. Indeed, recent textbooks primarily focus on discrete-state Markov jump processes (Peliti and Pigolotti 2021; Shiraishi 2023), which are not directly applicable to the continuous-state climate system. Moreover, Ito (2016b) explains the concept of information thermodynamics clearly despite its concise description. Although not straightforward, according to Ito (2016a), the theory in Ito (2016b) is equivalent to the theory used in our study (Horowitz and Esposito 2014; Loos and Klapp 2020).

## a. Multivariate nonlinear SDEs

We derive the second law of information thermodynamics for 𝑁 -variable nonlinear stochastic differential equations (SDEs, e.g., Gardiner 2009).

$$d x _ { i } = - a _ { i } ( { \mathbf x } ) d t + \sqrt { 2 D _ { i } } \mathrm d W _ { i }, \quad ( i = 1, \dots, N ),$$

where x : = ( 𝑥 , . . . , 𝑥 1 𝑁 ) T , 𝑎𝑖 ( x ) is a nonlinear function of x , 𝐷𝑖 is a positive constant, d 𝑊𝑖 is the differential of mutually independent Wiener processes for each 𝑖 . Equation (A1) assumes that the coefficient matrix D for the Wiener processes is diagonal: D = diag ( √ 2 𝐷 ,..., 1 √ 2 𝐷𝑁 ) . This assumption corresponds to the bipartite condition (Horowitz and Esposito 2014) and is essential for the following derivation. For example, this condition is necessary for the transformation from Eq. (A3) to (A4), and Eq. (A4) is indispensable for the final inequality (A32). If D is non-diagonal, a noise d 𝑊𝑖 may change not only 𝑥 𝑖 but also other state variables. By taking the matrix D to be diagonal, it is guaranteed that d 𝑊𝑖 contributes only to the fluctuation of 𝑥 𝑖 . Moreover, if D is diagonal, each component can be taken positive without loss of generality. We expressed it as √ 2 𝐷𝑖 to simplify later equations.

The probability distribution 𝑝 tot ( x , 𝑡 ) of the ensemble of solutions to these SDEs follows the Fokker-Planck equation (FPE, e.g., Gardiner 2009). The subscript 'tot' emphasizes that 𝑝 tot is the joint probability of all variables:

$$\frac { \partial p _ { \text{tot} } ( \mathbf x, t ) } { \partial t } = \sum _ { i = 1 } ^ { N } \frac { \partial } { \partial x _ { i } } \left [ a _ { i } ( \mathbf x ) p _ { \text{tot} } ( \mathbf x, t ) \right ] + \sum _ { i = 1 } ^ { N } D _ { i } \frac { \partial ^ { 2 } } { \partial x _ { i } ^ { 2 } } p _ { \text{tot} } ( \mathbf x, t ) = - \sum _ { i = 1 } ^ { N } \frac { \partial } { \partial x _ { i } } J _ { i } ( \mathbf x, t ). \quad ( A 2 )$$

In the rightmost side, the probability flux J is introduced, explicitly expressing the conservation of probability. The 𝑖 -th component of J is defined as follows:

$$J _ { i } ( { \mathbf x }, t ) \coloneqq - a _ { i } ( { \mathbf x } ) p _ { \text{tot} } ( { \mathbf x }, t ) - D _ { i } \frac { \partial } { \partial x _ { i } } p _ { \text{tot} } ( { \mathbf x }, t ).$$

Using 𝐽 𝑖 , the derivative 𝜕𝑝 tot / 𝜕𝑥𝑖 is expressed as:

$$\frac { \partial p _ { \text{tot} } } { \partial x _ { i } } = - \frac { J _ { i } } { D _ { i } } - \frac { p _ { \text{tot} } a _ { i } } { D _ { i } }.$$

For boundary conditions for the FPE, we assume 𝑝 tot ( x , 𝑡 ) = 0 and J x ( , 𝑡 ) = 0 at sufficiently far away from the origin.

## b. Second law of information thermodynamics for subsystems

Wederive the second law of information thermodynamics by transforming the time derivative of the Shannon entropy 𝑆𝜁 for a subsystem 𝑧 (: = 𝑥 𝜁 ) (Horowitz and Esposito 2014; Loos and Klapp

2020):

$$S _ { \zeta } \coloneqq - \int d z \, p _ { \zeta } ( z, t ) \ln p _ { \zeta } ( z, t ),$$

where 𝜁 represents any index from 1 to 𝑁 . We introduced 𝜁 to distinguish it from other subscripts, such as 𝑖 and 𝑗 , and to emphasize that 𝜁 corresponds to 𝑧 ( = 𝑥 𝜁 ) . The marginal probability 𝑝𝜁 is defined by the integral of the joint probability 𝑝 tot :

$$p _ { \zeta } ( z, t ) \coloneqq & \int \left ( \prod _ { i = 1, i \neq \zeta } ^ { N } d x _ { i } \right ) \, p _ { \text{tot} } ( \mathbf x, t ),$$

«

‹

The time evolution of this marginal probability 𝑝𝜁 is given by the integral of the convergence of the probability flux:

$$\frac { \partial } { \partial t } p _ { \zeta } ( z, t ) = \frac { \partial } { \partial t } \int \left ( \prod _ { i = 1, i \neq \zeta } ^ { N } \mathrm d x _ { i } \right ) \, p _ { \text{tot} } ( \mathbf x, t ),$$

$$\overset { \cdot } { \cdot } \quad & \overset { \cdot } { \cdot } \quad \overset { \cdot } { \cdot } \quad \underset { \ell } { \vdots } = - \int \left ( \prod _ { i = 1, i \neq \zeta } ^ { N } \mathrm d x _ { i } \right ) \sum _ { j = 1 } ^ { N } \frac { \partial } { \partial x _ { j } } J _ { j } ( \mathbf x, t ), \\ & \overset { \cdot } { \cdot } \int \overset { N } { \underline { N } } \quad \underset { \ell } { \lambda }$$

$$= - \int \left ( \prod _ { i = 1, i \neq \zeta } ^ { N } d x _ { i } \right ) \, \frac { \partial } { \partial z } J _ { \zeta } ( { \mathbf x }, t ).$$

«

‹

In the last equation, we integrate the derivative 𝜕𝐽𝑖 / 𝜕𝑥𝑖 ( 𝑖 ≠ 𝜁 ) and use the boundary conditions, leading it to be zero. The point of the following derivation is to use the FPE via the joint probability 𝑝 tot ( x , 𝑡 ) to deal with the marginal probability 𝑝𝜁 ( 𝑧, 𝑡 ) .

We transform the time derivative of the Shannon entropy:

$$\frac { d } { d t } S _ { \zeta } = - \int \frac { \partial p _ { \zeta } } { \partial t } \ln p _ { \zeta } \, \mathrm d z - \underbrace { \int \frac { \partial p _ { \zeta } } { \partial t } \mathrm d z } _ { = 0 }, \\ \int _ { 0 } \,. \quad.$$

$$= \int d z \int \left ( \prod _ { i = 1, i \neq \zeta } ^ { N } d x _ { i } \right ) \left [ \frac { \partial } { \partial z } J _ { \zeta } ( { \mathbf x }, t ) \right ] \ln p _ { \zeta } ( z, t ),$$

$$& \vec { \nu } \quad \sum _ { i = 1 } ^ { \nu } \lambda \lambda ^ { \lambda } \lambda ^ { \lambda } \lambda ^ { \lambda } \\ = & \int \underbrace { \left ( \prod _ { i = 1 } ^ { N } d x _ { i } \right ) } _ { = d x } \left [ \frac { \partial J _ { \zeta } ( \mathbf x, t ) } { \partial z } \right ] \left [ \ln \frac { p _ { \zeta } ( z, t ) } { p _ { \text{tot} } ( \mathbf x, t ) } + \ln p _ { \text{tot} } ( \mathbf x, t ) \right ],$$

=

d

x

$$= i _ { \zeta \leftarrow } + \int \mathrm d { \mathbf x } \, \frac { \partial J _ { \zeta } ( { \mathbf x }, t ) } { \partial z } \ln p _ { \text{tot} } ( { \mathbf x }, t ).$$

In Eq. (A10), the conservation of probability was used. In Eq. (A12), the definition of 𝑧 is used: 𝑧 : = 𝑥 𝜁 . We define the information flow / 𝐼 𝑖 ← as follows (Allahverdyan et al. 2009; Loos and Klapp 2020):

$$\dot { I } _ { i \leftarrow } \coloneqq \int \mathrm d { \mathbf x } \, \frac { \partial J _ { i } ( { \mathbf x }, t ) } { \partial x _ { i } } \ln \frac { p _ { i } ( x _ { i }, t ) } { p _ { \text{tot} } ( { \mathbf x }, t ) },$$

where the marginal distribution for 𝑥 𝑖 is

$$p _ { i } ( x _ { i }, t ) \coloneqq \int \left ( \prod _ { j = 1, j \neq i } ^ { N } \mathrm d x _ { j } \right ) p _ { \text{tot} } ( \mathbf x, t ).$$

«

‹

The information flow / 𝐼 𝑖 ← represents changes in the mutual information by the subsystem 𝑥 𝑖 . We further discuss information flow in Subsection c in Appendix A. In Eq. (A13), the information flow / 𝐼 𝜁 ← (i.e., 𝑖 = 𝜁 ) was substituted.

We further transform the second term in Eq. (A13):

$$\int \mathrm d { \mathbf x } \, \frac { \partial J _ { \zeta } } { \partial z } \ln p _ { \text{tot} } = - \int \mathrm d { \mathbf x } \, J _ { \zeta } \frac { \partial \left ( \ln p _ { \text{tot} } \right ) } { \partial z } + \int \left ( \prod _ { i = 1, i \neq \zeta } ^ { N } \mathrm d { x } _ { i } \right ) \, \underbrace { \left [ J _ { \zeta } \ln p _ { \text{tot} } \right ] _ { z = - \infty } ^ { z = + \infty } } _ { = 0 },$$

$$= - \int \mathrm d { \mathrm x } \, \frac { J _ { \zeta } } { p _ { \text{tot} } } \frac { \partial p _ { \text{tot} } } { \partial z },$$

$$= \left \langle - \frac { 1 } { p _ { \text{tot} } } \frac { \dot { \partial } p _ { \text{tot} } } { \partial z } \circ \frac { d z } { \text{d} t } \right \rangle,$$

$$= \dot { \left \langle \frac { 1 } { p _ { \text{tot} } } } \left ( \frac { J _ { \zeta } } { D _ { \zeta } } + \frac { p _ { \text{tot} } \overset { \cdot } { a _ { \zeta } } } { D _ { \zeta } } \right ) \circ \frac { d z } { d t } \right \rangle,$$

$$= \stackrel {. } { \left \langle \frac { J _ { \zeta } } { p _ { \text{tot} } D _ { \zeta } } } \circ \frac {. } { d t } \right \rangle + \left \langle \frac { a _ { \zeta } } { D _ { \zeta } } \circ \frac { d z } { d t } \right \rangle,$$

where ◦ represents the Stratonovich product (e.g., Gardiner 2009) and Eq. (A4) with 𝑖 = 𝜁 is substituted in Eq. (A18) to replace 𝜕𝑝 tot / 𝜕𝑧 . We further explain the Stratonovich product and the conversion formula between 𝐽 𝜁 and d 𝑧 / d 𝑡 used in transforming Eq. (A17) to (A18).

First, the Stratonovich product is defined as follows (e.g., Gardiner 2009):

$$f ( { \mathbf x } ( t ) ) \circ d W \coloneqq \lim _ { \Delta t \to 0 } \frac { f ( { \mathbf x } ( t + \Delta t ) ) + f ( { \mathbf x } ( t ) ) } { 2 } \left ( W _ { t + \Delta t } - W _ { t } \right ),$$

where 𝑓 represents any smooth function. Using the Stratonovich product, the chain rules hold mathematically (Gardiner 2009), and the energy conservation law (i.e., the first law of thermodynamics) holds physically (Sekimoto 2010). The introduction of the Stratonovich product is also discussed in Shiraishi (2023). Because of these properties, the Stratonovich product is commonly used in information thermodynamics. The Stratonovich product is necessary only for terms where the differential d 𝑊 appears. For example, using Eq. (A21), the product of 𝑓 ( x ( 𝑡 )) and d 𝑧 / d 𝑡 becomes:

$$f ( \mathbf x ( t ) ) \circ \frac { \mathrm d z } { \mathrm d t } = f ( \mathbf x ( t ) ) \circ \left ( - a _ { \zeta } + \sqrt { 2 D _ { \zeta } } \frac { \mathrm d W _ { \zeta } } { \mathrm d t } \right ),$$

$$= - f ( \mathbf x ( t ) ) a _ { \zeta } + \sqrt { 2 D _ { \zeta } } f ( x ( t ) ) \circ \frac { \mathrm d W _ { \zeta } } { \mathrm d t },$$

where d 𝑊𝜁 / d cannot be defined mathematically because of the non-smoothness of 𝑡 𝑊𝜁 ( 𝑡 ) (Gardiner 2009), but in physics, d 𝑊𝜁 / d 𝑡 is often used for computational convenience as white Gaussian noise.

In the transformation from Eq. (A17) to (A18), the following formula was used (Seifert 2012):

$$\left \langle f ( { \mathbf x } ( t ) ) \circ \frac { d x _ { i } } { d t } \right \rangle = \int d { \mathbf x } \, f ( { \mathbf x } ) J _ { i } ( { \mathbf x }, t ).$$

The statistical average of the rate of change of the state variable d 𝑥 𝑖 / d 𝑡 becomes the integral of the 𝑖 -th component of the probability flux 𝐽 𝑖 ( x , 𝑡 ) in the phase space. To show this formula, we use the definition of the Stratonovich product [Eq. (A21)] and perform a Taylor expansion. To simplify the derivation, we ignore terms of order smaller than d . 𝑡

$$d t \left \langle f ( \mathbf x ( t ) ) \circ \frac { d x _ { i } } { d t } \right \rangle = \left \langle f ( \mathbf x ( t ) ) \circ d x _ { i } \right \rangle,$$

$$\L ^ { \prime } = \left \langle \left ( f + \sum _ { k = 1 } ^ { N } \frac { 1 } { 2 } \frac { \partial f } { \partial x _ { k } } \mathrm d x _ { k } \right ) \mathrm d x _ { i } \right \rangle,$$

$$= \left \langle f ( x ) \mathrm d x _ { i } + \frac { 1 } { 2 } \frac { \partial f } { \partial x _ { i } } \left ( \mathrm d x _ { i } \right ) ^ { 2 } \right \rangle,$$

$$= \overset { \cdot } { d t } \left \langle - a _ { i } f + D _ { i } \frac { \partial f } { \partial x _ { i } } \right \rangle,$$

$$= d t \, \int \mathrm d { \mathbf x } \, J _ { i } ( { \mathbf x }, t ) f ( { \mathbf x } ).$$

From Eq. (A25) to (A26), the definition of the Stratonovich product was used, and a Taylor expansion was performed. From Eq. (A26) to (A28), the Ito rule was used, where ⟨ d 𝑥 𝑘 d 𝑥 𝑖 ⟩ = 0 for 𝑘 ≠ 𝑖 due to the diagonal matrix D in the SDEs [Eq. (A1)]. In Eq. (A29), after integration by parts, the definition of the probability flux [Eq. (A3)] was substituted. Finally, Eq. (A24) is derived by comparing Eqs. (A25) and (A29).

To rewrite the numerator of the second term in Eq. (A20), we define the following quantity:

$$\dot { Q } _ { i } \coloneqq - \left \langle a _ { i } ( { \mathbf x } ) \circ \frac { d x _ { i } } { d t } \right \rangle,$$

$$= \left \langle \left ( \frac { \overset { \cdot } { d x _ { i } } } { d t } - \xi _ { i } \right ) \circ \frac { \overset { \cdot } { d x _ { i } } } { d t } \right \rangle,$$

where 𝜉 𝑖 : = √ 2 𝐷𝑖 d 𝑊𝑖 / d 𝑡 and the SDE for 𝑖 [Eq. (A1)] is used to obtain the second line. The quantity / 𝑄𝑖 can be regarded as the amount of heat dissipated by the subsystem 𝑥 𝑖 to the surrounding environment per unit time (Sekimoto 2010; Loos and Klapp 2020).

Combining Eqs. (A13), (A20), (A24), and (A30), we obtain

$$\underbrace { \frac { d S _ { \zeta } } { d t } + \frac { \dot { Q } _ { \zeta } } { D _ { \zeta } } } _ { = ; \dot { \sigma } _ { \zeta } } - \dot { I } _ { \zeta \leftarrow } = \left \langle \frac { J _ { \zeta } } { p _ { \text{tot} } D _ { \zeta } } \circ \frac { d z } { d t } \right \rangle = \int \mathrm d { \mathbf x } \, \frac { J _ { \zeta } ^ { 2 } } { p _ { \text{tot} } D _ { \zeta } } \geq 0,$$

where the equality holds when 𝐽 𝜁 = 0, that is, in the case of the detailed balance being satisfied for 𝑧 (Loos and Klapp 2020; Gardiner 2009). Since the subscript 𝜁 is arbitrary, the following inequality is derived:

$$\dot { \sigma } _ { i } = \frac { d S _ { i } } { d t } + \frac { \dot { Q } _ { i } } { D _ { i } } \geq \dot { I } _ { i \leftarrow }.$$

We defined the production rate of total entropy / 𝜎𝑖 as the sum of the rate of change of the Shannon entropy d 𝑆𝑖 / d 𝑡 of the subsystem 𝑥 𝑖 and the rate of change of the heat / 𝑄𝑖 / 𝐷𝑖 dissipated by the subsystem 𝑥 𝑖 to the surroundings. The production rate of total entropy / 𝜎𝑖 is lower bounded by the information flow / 𝐼 𝑖 ← . If / 𝐼 𝑖 ← is negative, / 𝜎𝑖 can be negative, that is, the total entropy can decrease. In the usual thermodynamic setup, 𝐷𝑖 is equal to the temperature of the environment (Sekimoto 2010; Loos and Klapp 2020). Thus, Eq. (A33) is considered a generalization of the Clausius inequality. The equality holds when 𝐽 𝑖 = 0. In this case, the detailed balance is satisfied mathematically (Gardiner 2009), and the subsystem 𝑥 𝑖 becomes an equilibrium state physically (Loos and Klapp 2020).

In the steady state, the time derivative of the Shannon entropy becomes zero, and the following holds:

$$\dot { \sigma } _ { i } = \frac { \dot { Q } _ { i } } { D _ { i } } \geq \dot { I } _ { i \leftarrow }.$$

The steady state is defined as a state in which joint probabilities of state variables between any time points depend only on the time differences (i.e., lags) and not on the absolute time itself. As a result, in the steady state, the marginal probability 𝑝𝑖 ( 𝑥 𝑖 , 𝑡 ) does not change over time; thus, the time derivative of the Shannon entropy 𝑆𝑖 becomes zero.

## c. Conservation law of information flows

We derive the conservation law of information flows (Allahverdyan et al. 2009; Loos and Klapp 2020). We first define the mutual information 𝑀 (e.g., Cover and Thomas 2005):

$$M \coloneqq \int \mathrm d { \mathbf x } \, p _ { \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \boldsymbol \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \cdot } } \ln \left [ \frac { p _ { \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \mathrm \boldsymbol } } } { \prod _ { i = 1 } ^ { N } p _ { i } ( x _ { i }, t ) } \right ]$$

Mutual information represents the magnitude of nonlinear correlations between subsystems. The value of mutual information is non-negative, and it becomes zero only when the subsystems are independent of each other.

We transform the time derivative of mutual information.

$$\frac { d M } { d t } = \int \mathrm d { \mathbf x } \, \frac { \partial p _ { \text{tot} } } { \partial t } \ln \left ( \frac { p _ { \text{tot} } } { \prod _ { i = 1 } ^ { N } p _ { i } } \right ) + \int \mathrm d { \mathbf x } \, p _ { \text{tot} } \left ( \frac { \partial } { \partial t } \ln p _ { \text{tot} } - \sum _ { i = 1 } ^ { N } \frac { \partial } { \partial t } \ln p _ { i } \right ),$$

$$= \sum _ { j = 1 } ^ { N } \left [ \int \mathrm d { \mathbf x } \, \frac { \partial J _ { j } } { \partial x _ { j } } \ln p _ { j } - \int \mathrm d { \mathbf x } \, \frac { \partial J _ { j } } { \partial x _ { j } } \ln p _ { \mathrm t o t } \right ],$$

$$\frac { d M } { d t } & = \int \mathrm d x \, \frac { \partial p _ { \text{totd} } } { \partial t } \ln \left ( \frac { p _ { \text{totd} } } { \prod _ { i = 1 } ^ { N } p _ { i } } \right ) + \int \mathrm d x \, p _ { \text{totd} } \left ( \frac { \partial } { \partial t } \ln p _ { \text{totd} } - \sum _ { i = 1 } ^ { N } \frac { \partial } { \partial t } \ln p _ { i } \right ), \\ & = \int \mathrm d x \, \frac { \partial p _ { \text{totd} } } { \partial t } \ln \left ( \frac { p _ { \text{totd} } } { \prod _ { i = 1 } ^ { N } p _ { i } } \right ) + \underbrace { \int \mathrm d x \, \frac { \partial p _ { \text{totd} } } { \partial t } } _ { = 0 } - \sum _ { i = 1 } ^ { N } \underbrace { \int \mathrm d x _ { i } } \frac { \partial p _ { i } } { \partial t }, \\ & = - \int \mathrm d x \left ( \sum _ { j = 1 } ^ { N } \frac { \partial J _ { j } ( \mathbf x, t ) } { \partial x _ { j } } \right ) \ln \left [ \frac { p _ { \text{totd} } ( \mathbf x, t ) } { \prod _ { i = 1 } ^ { N } p _ { i } ( x _ { i }, t ) } \right ], \\ & = \sum _ { j = 1 } ^ { N } \left [ \int \mathrm d x \, \frac { \partial J _ { j } } { \partial x _ { j } } \left ( \sum _ { i = 1 } ^ { N } \ln p _ { i } \right ) - \int \mathrm d x \, \frac { \partial J _ { j } } { \partial x _ { j } } \ln p _ { \text{totd} } \right ], \\ & = \sum _ { j = 1 } ^ { N } \left [ \int \mathrm d x \, \frac { \partial J _ { j } } { \partial x _ { j } } \ln p _ { j } - \int \mathrm d x \, \frac { \partial J _ { j } } { \partial x _ { j } } \ln p _ { \text{totd} } \right ],$$

In Eq. (A37), the conservation of probability and the definition of marginal probability were used. From Eq. (A39) to (A40), for all 𝑖 such that 𝑖 ≠ 𝑗 , the integral of 𝐽 𝑗 is executable and the result becomes zero due to the boundary conditions. Consequently, only the terms with 𝑖 = 𝑗 remain in Eq. (A40).

The conservation law of information flows is obtained by substituting Eq. (A14) into (A40):

$$\frac { \mathrm d M } { \mathrm d t } = \sum _ { i = 1 } ^ { N } \dot { I } _ { i \leftarrow }.$$

The sum of changes in the mutual information by the subsystem 𝑥 𝑖 (i.e., the information flow / 𝐼 𝑖 ← ) over all 𝑖 becomes the temporal change of the mutual information of the entire system. The information flow / 𝐼 𝑖 ← is interpreted as the sum of the inflow and outflow of information from all the other subsystems to 𝑥 𝑖 .

In the steady state, since the joint probability 𝑝 tot does not change over time, the time derivative of mutual information becomes zero. As a result, the following holds:

$$0 = \sum _ { i = 1 } ^ { N } \dot { I } _ { i \leftarrow \cdot \cdot }.$$

The source of information flow is uniquely determined when the entire system consists of only two variables (Allahverdyan et al. 2009). Here, the two random variables are denoted by uppercase letters 𝑋𝑡 and 𝑌𝑡 , and their realizations are denoted by lowercase letters 𝑥 𝑡 and 𝑦 𝑡 , respectively, where the time dependence is explicitly denoted by the subscript 𝑡 . In this case, the mutual information [Eq. (A35)] is expressed as follows:

$$M ( X _ { t }, Y _ { t } ) \coloneqq \int \mathrm d x _ { t } \mathrm d y _ { t } \ p ( x _ { t }, y _ { t } ) \ln \left [ \frac { p ( x _ { t }, y _ { t } ) } { p ( x _ { t } ) p ( y _ { t } ) } \right ].$$

The dependence on the variables 𝑋𝑡 and 𝑌𝑡 is explicitly written inside the parentheses as 𝑀 𝑋𝑡,𝑌𝑡 ( ) , where 𝑀 𝑋𝑡,𝑌𝑡 ( ) is symmetric with respect to 𝑋𝑡 and 𝑌𝑡 . In the bivariate case, the information flows are expressed as follows (Allahverdyan et al. 2009):

$$\dot { I } _ { X \leftarrow Y } = \lim _ { \Delta t \rightarrow 0 } \frac { M ( X _ { t + \Delta t }, Y _ { t } ) - M ( X _ { t }, Y _ { t } ) } { \dots \dots \dots \dots \dots },$$

$$i _ { Y \leftarrow X } = \lim _ { \Delta t \rightarrow 0 } \frac { M ( X _ { t }, Y _ { t + \Delta t } ) - M ( X _ { t }, Y _ { t } ) } { \Delta t }.$$

Based on the Leibniz rule for differentiation, the sum of these information flows satisfies the conservation law (A41). The variation of mutual information due to the time change of the random variable 𝑋 (or 𝑌 ) is defined as the information flow to 𝑋 (or 𝑌 ). Equation (A44) can be transformed into / 𝐼 𝑋 ← defined by Eq. (A14), and Eq. (A45) is also transformed into / 𝐼 𝑌 ← in the same manner (Allahverdyan et al. 2009). From the conservation law of information flows in the steady state, namely / 𝐼 𝑋 ← + / 𝑌 𝐼 𝑌 ← 𝑋 = 0, we see that information flows can be negative or positive.

## d. Second law of thermodynamics for the entire system

We derive the second law of thermodynamics for the entire system and confirm that focusing on subsystems is essential in information thermodynamics (Sagawa and Ueda 2013b). As a preparation, we define the Shannon entropy of the entire system:

$$S _ { \text{tot} } \coloneqq - \int \text{d} \mathbf x \, p _ { \text{tot} } ( \mathbf x, t ) \ln p _ { \text{tot} } ( \mathbf x, t ).$$

In general, additivity does not hold for Shannon entropy. The degree to which additivity does not hold is evaluated by mutual information.

$$\left ( \sum _ { i = 1 } ^ { N } S _ { i } \right ) - M & = - \left [ \sum _ { i = 1 } ^ { N } \int d x _ { i } \ p _ { i } ( x _ { i }, t ) \ln p _ { i } ( x _ { i }, t ) \right ] - \int d x \ p _ { t o t } ( \mathbf x, t ) \ln \left [ \frac { p _ { t o t } ( \mathbf x, t ) } { \prod _ { i = 1 } ^ { N } p _ { i } ( x _ { i }, t ) } \right ], \quad ( A 4 7 ) \\ & = - \left ( \sum _ { i = 1 } ^ { N } \int d \mathbf x \ p _ { t o t } \ln p _ { i } \right ) - \int d \mathbf x p _ { t o t } \ln \left ( \frac { p _ { t o t } } { \prod _ { i = 1 } ^ { N } p _ { i } } \right ), \\ & = - \int d \mathbf x \ p _ { t o t } \ln p _ { t o t }, \\ & = S _ { t o t },$$

where we use the definitions of 𝑆𝑖 , 𝑀 , and 𝑆 tot in Eqs. (A5), (A35), and (A46), respectively. From Eq. (A47) to (A48), the definition of marginal probability [Eq. (A15)] was applied to the first term. The Shannon entropies of the subsystems are additive only when the mutual information 𝑀 is zero.

We derive the second law for the entire system. Since Eq. (A33) holds for any 𝑖 , it can be summed over all 𝑖 . Applying Eqs. (A41) and (A50), we obtain the following:

$$\dot { \sigma } _ { \text{tot} } \coloneqq \frac { \text{d} S _ { \text{tot} } } { \text{d} t } + \sum _ { i = 1 } ^ { N } \frac { \dot { Q } _ { i } } { D _ { i } } \geq 0.$$

The first term in the middle represents the rate of change of the Shannon entropy of the entire system, and the second term represents the sum of the rates of change of the heat dissipated to the environment. In Eq. (A51), the production rate of total entropy for the entire system is lower bounded by zero, whereas in Eq. (A33), the production rate of total entropy for a subsystem is

lower bounded by the information flow, which can be negative. This information flow appears due to the non-additivity of Shannon entropy (Sagawa and Ueda 2013b). Historically, the second law of thermodynamics [Eq. (A51)] was derived first (Seifert 2005), and the second law of information thermodynamics [Eq. (A33)] was derived later (Horowitz and Esposito 2014; Loos and Klapp 2020).

## e. Interpretation of heat

The quantity / 𝑄𝑖 in Eq. (A30) is regarded as the heat dissipated by a subsystem 𝑥 𝑖 to the environment per unit time (Sekimoto 2010). When the environment is always in a thermal equilibrium state, this heat is equal to the entropy production rate of the environment (Peliti and Pigolotti 2021; Shiraishi 2023). However, this interpretation is not always possible.

Generally, in systems described by stochastic differential equations, noise is not always due to thermal fluctuations. For example, in the dynamical system considered in the current research [Eqs. (16) and (17)], noise represents the influence of small-scale disturbances in the atmosphere and ocean (Section 3c). Moreover, it is unclear whether such small-scale disturbance fields are always in an equilibrium state, although this point may be investigated using equilibrium statistical mechanics for geophysical fluids (Majda and Wang 2006; Bouchet and Venaille 2012; Campa et al. 2014). Indeed, the meanders of the Gulf Stream and the Kuroshio Current are consistent with analysis using equilibrium statistical mechanics (Venaille and Bouchet 2011), which implies that small-scale disturbance fields in both current regions can also be analyzed using equilibrium theory. Currently, however, such an analysis is not sufficient, and the interpretation of / 𝑄𝑖 needs to be considered for each application. In particular, / 𝜎𝑖 = / 𝑄𝑖 / 𝐷𝑖 holds in the steady state [Eq. (A34)], making this interpretation essential.

The present study has analyzed the dynamical system [Eqs. (16) and (17)] using the second law of information thermodynamics as a mathematical tool (Section 4). Although the second law is always derived for SDEs as shown above, mathematically deriving the second law is different from obtaining a meaningful physical insight from it. We have regarded / 𝜎𝐺 ( = / 𝑄𝐺 𝐷𝐺 / ) and / 𝜎𝐾 ( = / 𝑄𝐾 / 𝐷𝐾 ) as indicators of the contribution of interactions to the temporal tendencies of the SSTs. This interpretation does not rely on the interpretation of / 𝑄𝑖 as heat. As a similar example, Ito and Sagawa (2015) considered / 𝑄𝑖 as the robustness of signal transmission in a cell (i.e., not as

heat) and analyzed the processing efficiency of cell signaling using information thermodynamics. The application of information thermodynamics in this manner might yield interesting results for other atmospheric and oceanic systems. Investigating further applications is an important topic for future research in atmospheric and oceanic sciences.

## APPENDIX B

## Theoretical formulas used in Sections 3 and 4

## a. Correspondence between linear dynamical system and AR1 model

We convert the AR1 model to the linear dynamical system by comparing the means and covariances. The obtained theoretical formulas were used in Section 3b.

The linear dynamical system is described in the vector form as follows (Gardiner 2009):

$$d x & = - \mathbf A \, x d t + \mathbf B d W,$$

where x is an 𝑁 -dimensional real vector variable, time 𝑡 is continuous, and W represents 𝑁 Wiener processes that are independent of each other. We consider here a general case of 𝑁 ≥ 1, although the case of 𝑁 = 2 was used in Section 3. The constant real matrix A is assumed to have only eigenvalues with positive real parts, and B is a constant real matrix. Under these conditions, x has a stationary solution (Gardiner 2009). Integrating Eq. (B1) over a time step Δ 𝑡 , the following solution is obtained (Gardiner 2009):

$$\mathbf x ( t + \Delta t ) = e ^ { - \mathbf A \Delta t } \mathbf x ( t ) + \int _ { 0 } ^ { \Delta t } e ^ { - \mathbf A ( \Delta t - s ) } \mathbf B \, d \mathbf W ( s ).$$

The multivariate AR1 model is given as follow (Hamilton 1994):

$$y _ { t + \Delta t } = \mathbf C y _ { t } + \varepsilon _ { t },$$

where y is an 𝑁 -dimensional real vector variable and 𝜀 represents vector white noise. The absolute values of all eigenvalues of the real matrix C are assumed to be smaller than 1, which is the condition for stationarity (Hamilton 1994). Time 𝑡 is discrete and takes steps at a constant interval

## Δ 𝑡 . The covariance matrix V of the vector white noise is generally non-diagonal and is expressed as

$$\langle \varepsilon _ { t } \varepsilon _ { s } ^ { T } \rangle = \begin{cases} \mathbf v & ( t = s ), \\ 0 & ( \text{otherwise} ), \end{cases}$$



where ⟨·⟩ represents the statistical average over noise.

We obtain the condition for the mean values of x ( 𝑡 + Δ 𝑡 ) and y 𝑡 + Δ 𝑡 being equal:

$$\mathbf C = e ^ { - \mathbf A \Delta t },$$

where x ( 𝑡 ) = y 𝑡 is assumed. We also obtain the condition for the covariances of x ( 𝑡 + Δ 𝑡 ) and y 𝑡 + Δ 𝑡 being equal:

$$\mathbf V = \int _ { 0 } ^ { \Delta t } d s \, e ^ { - \mathbf A ( \Delta t - s ) } \mathbf B \mathbf B ^ { T } e ^ { - \mathbf A ^ { T } ( \Delta t - s ) }.$$

Under stationarity, Eqs. (B5) and (B6) become independent of 𝑡 .

The regression analysis with the AR1 model determines the matrices C and V , where the matrix V is estimated from the residuals of regression. Using Eqs. (B5) and (B6), A and B are then obtained. This method requires that the AR1 model and the dynamical system share the same statistical properties over the time interval Δ 𝑡 . Since the linear SDE [Eq. (B1)] describes a Gaussian process, its parameters are uniquely determined from the mean and covariance. Table 1 shows the coefficients of the dynamical system [Eqs. (16) and (17)], which were obtained by applying Eqs. (B5) and (B6) to the regression result with AR1.

## b. Theoretical expressions of lag correlation coefficients

We show here the theoretical expressions of lag correlation coefficients used in Sections 3d and 4b. The dynamical system [Eqs. (16) and (17)] describes a Gaussian process, so the statistical properties of the steady state are completely determined by the lag covariance matrix (Gardiner 2009). By normalizing the lag covariance, we obtain the lag correlation coefficient, which is an important indicator for the synchronization between the Gulf Stream and the Kuroshio Current (Kohyama et al. 2021).

We first show the instantaneous covariance matrix:

$$\mathfrak { z } = \begin{pmatrix} \Sigma _ { G G } & \Sigma _ { G K } \\ \Sigma _ { K G } & \Sigma _ { K K } \end{pmatrix} = \begin{pmatrix} \langle T _ { G } T _ { G } \rangle & \langle T _ { G } T _ { K } \rangle \\ \langle T _ { K } T _ { G } \rangle & \langle T _ { K } T _ { K } \rangle \end{pmatrix},$$

«

‹

«

‹

where the off-diagonal components are equal, Σ 𝐺𝐾 = Σ 𝐾𝐺 . Each component is given as follows:

$$\Sigma _ { G G } = \frac { ( r _ { K } ^ { 2 } + r _ { G } r _ { K } - c _ { G \leftarrow K } c _ { K \leftarrow G } ) D _ { G } + c _ { G \leftarrow K } ^ { 2 } D _ { K } } { ( r _ { G } + r _ { K } ) ( r _ { G } r _ { K } - c _ { G \leftarrow K } c _ { K \leftarrow G } ) },$$

$$\Sigma _ { K K } = \frac { ( r _ { G } ^ { 2 } + r _ { G } r _ { K } - c _ { G \leftarrow K } c _ { K \leftarrow G } ) D _ { K } + c _ { K \leftarrow G } ^ { 2 } D _ { G } } { ( r _ { G } + r _ { K } ) ( r _ { G } r _ { K } - c _ { G \leftarrow K } c _ { K \leftarrow G } ) },$$

$$\Sigma _ { G K } = \Sigma _ { K G } = \frac { \lambda \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \colon \\ ( r _ { G } + r _ { K } ) ( r _ { G } r _ { K } - c _ { G \leftarrow K } c _ { K \leftarrow G } ) } { ( r _ { G } r _ { K } - c _ { G \leftarrow K } c _ { K \leftarrow G } ) }.$$

Details of the derivation for the instantaneous covariance can be found in Gardiner (2009), directly leading to Eqs. (B8)-(B10).

Using the instantaneous covariance Σ , the lag covariance matrix is given as follows (Gardiner 2009):

$$\langle T _ { i } ( t _ { 1 } ) T _ { j } ( t _ { 2 } ) \rangle = \begin{cases} \left ( e ^ { - \mathbf A ( t _ { 1 } - t _ { 2 } ) } \mathbf S \right ) _ { i j } & ( t _ { 1 } \geq t _ { 2 } ), \\ \left ( \mathbf S e ^ { - \mathbf A ^ { T } ( t _ { 2 } - t _ { 1 } ) } \right ) _ { i j } & ( t _ { 1 } < t _ { 2 } ), \end{cases}$$

where where Δ 𝑡 represents the lag.

$$\mathbf A = \begin{pmatrix} r _ { G } & - c _ { G \leftarrow K } \\ - c _ { K \leftarrow G } & r _ { K } \end{pmatrix}.$$

«

‹

The indices 𝐺 and 𝐾 are subscripts for the Gulf Stream and the Kuroshio Current, respectively. In Section 3d, Fig. 4 shows the lag correlation coefficients computed by normalizing the lag covariances with the instantaneous variances:

$$\rho ( \Delta t ) = \frac { ( e ^ { - \mathbf A \Delta t } \mathbf \Sigma ) _ { G K } } { \sqrt { \Sigma _ { G G } \Sigma _ { K K } } } \quad ( \Delta t \geq 0 ),$$

$$\rho ( \Delta t ) = \frac { \stackrel { \dots, - \infty \infty - \infty } { \left ( \Sigma e ^ { - \mathbf A ^ { T } | \Delta t | } \right ) } _ { G K } } { \sqrt { \Sigma _ { G G } \Sigma _ { K K } } } \quad ( \Delta t < 0 ),$$

## c. Theoretical expressions of information thermodynamic quantities

This section derives the expressions for information thermodynamic quantities, such as / 𝜎𝐺 , used in Section 4. In other words, the equations shown below are the specific forms of the theoretical equations in Appendix A applied to the dynamical system, Eqs. (16) and (17). Details of the derivation are discussed in Loos and Klapp (2020).

As a preparation, we first examine the expressions for the variances and covariances. In the steady state, the probability distribution does not change in time, so the following formulas hold:

$$\frac { 1 } { 2 } \frac { d } { d t } \langle T _ { G } T _ { G } \rangle = \left \langle T _ { G } \circ \frac { d T _ { G } } { d t } \right \rangle = 0,$$

$$\frac { 1 } { 2 } \frac { \mathrm d } { \mathrm d t } \langle T _ { K } T _ { K } \rangle = \left \langle T _ { K } \circ \frac { \mathrm d T _ { K } } { \mathrm d t } \right \rangle = 0,$$

$$\frac { \mathrm d } { \mathrm d t } \langle T _ { G } T _ { K } \rangle = \overset { \cdot } { \left \langle } T _ { G } \circ \frac { \mathrm d T _ { K } } { \mathrm d t } \right \rangle + \left \langle T _ { K } \circ \frac { \mathrm d T _ { G } } { \mathrm d t } \right \rangle = 0.$$

The temporal derivatives d 𝑇𝐺 / d 𝑡 or d 𝑇𝐾 / d 𝑡 include white Gaussian noise 𝜉 𝐺 or 𝜉 𝐾 , respectively. Thus, the Stratonovich product ◦ needs to be used (Appendix A). The symbol ◦ is omitted in Eqs. (24)-(28), (31), and (32) because the Stratonovich product is not explained in Section 4. By substituting the governing equations (16) and (17), each term in Eq. (B17) is given by

$$\left \langle T _ { G } \circ \frac { d T _ { K } } { \mathrm d t } \right \rangle = - r _ { K } \langle T _ { K } T _ { G } \rangle + c _ { K \leftarrow G } \langle T _ { G } T _ { G } \rangle,$$

$$\left \langle T _ { K } \circ \frac { d T _ { G } } { \mathrm d t } \right \rangle = - r _ { G } \langle T _ { G } T _ { K } \rangle + c _ { G \leftarrow K } \langle T _ { K } T _ { K } \rangle,$$

where the variances and covariances on the right-hand sides are given by Eqs. (B8)-(B10).

We next provide the expressions for the entropy production rates. From Eqs. (16), (17), (A30), and (A34), the following equations are derived for the steady state:

$$\dot { \sigma } _ { G } = \frac { \dot { Q } _ { G } } { D _ { G } } = - \frac { r _ { G } } { D _ { G } } \left \langle T _ { G } \circ \frac { d T _ { G } } { d t } \right \rangle + \frac { c _ { G } \leftarrow K } { D _ { G } } \left \langle T _ { K } \circ \frac { d T _ { G } } { d t } \right \rangle,$$

$$= \frac { c _ { G \leftarrow K } } { D _ { G } } \left \langle T _ { K } \circ \frac { d T _ { G } } { d t } \right \rangle,$$

$$\dot { \sigma } _ { K } = \frac { \dot { Q } _ { K } } { D _ { K } } = - \frac { \lambda } { D _ { K } } \Big < \binom { \cdot } { T _ { K } } \circ \frac { d T _ { K } } { d t } \Big > + \frac { c _ { K \leftarrow G } } { D _ { K } } \Big < T _ { G } \circ \frac { d T _ { K } } { d t } \Big >,$$

$$= \frac { c _ { K \leftarrow G } } { D _ { K } } \left \langle T _ { G } \circ \frac { \mathrm d T _ { K } } { \mathrm d t } \right \rangle,$$

where Eqs. (B15) and (B16) are substituted. Equations (B21) and (B23) are the same as Eqs. (24) and (25), respectively. Equation (29) is obtained using Eqs. (B17), (B21), and (B23).

Theexpressions for the information flows are easily derived in the steady state. Using the property that the probability distribution does not change in time, the time derivative of the Shannon entropy is zero. Thus, from Eq. (A13), the following holds in general:

$$i _ { i \leftarrow } = - \int \mathfrak { d } \mathbf x \, \frac { \partial J _ { i } } { \partial x _ { i } } \ln p _ { \text{tot} },$$

$$& = \int d { \mathbf x } \, J _ { i } \frac { \partial \ln p _ { \text{tot} } } { \partial x _ { i } } + \int \left ( \prod _ { j = 1, j \neq i } ^ { N } d { \mathbf x } _ { j } \right ) \underbrace { [ J _ { i } p _ { \text{tot} } ] _ { x _ { i } = - \infty } ^ { x _ { i } = \infty } } _ { = 0 },$$

$$= \left \langle \frac { \partial \ln p _ { \text{tot} } } { \partial x _ { i } } \circ \frac { d x _ { i } } { d t } \right \rangle,$$

where the boundary conditions and Eq. (A24) are used here. We only need to evaluate the derivatives of ln 𝑝 tot in Eq. (B26). For a linear SDE, the stationary distribution is a Gaussian (Gardiner 2009), and the derivative of its logarithm can be expressed using the variances and covariances. The information flows are then derived as

$$i _ { G \leftarrow K } = - ( \mathbf S ^ { - 1 } ) _ { K G } \left \langle T _ { K } \circ \frac { d T _ { G } } { d t } \right \rangle,$$

$$i _ { K \leftarrow G } = - ( \mathbf S ^ { - 1 } ) _ { G K } \overset {. } { \left < T _ { G } \circ \frac { d T _ { K } } { d t } \right > },$$

where Eqs. (B15) and (B16) are substituted, and the inverse of the covariance matrix Σ is denoted as Σ -1 . The off-diagonal components of Σ -1 are ( Σ -1 ) 𝐾𝐺 and ( Σ -1 ) 𝐺𝐾 , which are equal because Σ is symmetric. Using the specific expression of this off-diagonal component, we obtain Eqs. (26) and (27).

The information flows and entropy production rates become zero when 𝑐 𝐺 ← 𝐾 𝐷 𝐾 = 𝑐 𝐾 ← 𝐺 𝐷 𝐺 (Loos and Klapp 2020). This fact can be directly verified by substituting Eqs. (B8)-(B10), (B18), and (B19) into Eqs. (24)-(27). A simpler method to show this is through the application of the second law of information thermodynamics (see Appendix D). In the case of 𝑐 𝐺 ← 𝐾 𝐷 𝐾 = 𝑐 𝐾 ← 𝐺 𝐷 𝐺 , the two subsystems, namely the Gulf Stream and the Kuroshio Current, become symmetric. Mathematically, the detailed balance is satisfied and the probability flux J is zero (Loos and Klapp 2020).

## APPENDIX C

## Histograms of estimated quantities by the moving block bootstrap method

We quantified the uncertainty of the estimated coefficients using the moving block bootstrap method (Mudelsee 2014) in Section 3b, where the number of resampled time series was set to 2000. We discuss here the histograms of the entropy production rates and information flows obtained from the coefficients estimated by the bootstrap method.

Figure C1 shows the histograms of all estimated quantities, including the coefficients of the dynamical system, where the number of samples in each histogram is 2000. The 95% confidence intervals in Table 1 were obtained from these histograms. We focus on the entropy production rates and information flows, which are not listed in Table 1. The conservation of information flows, / 𝐼 𝐺 ← 𝐾 + / 𝐼 𝐾 ← 𝐺 = 0 [Eq. (23)], always holds. That is, for each resampled time series, / 𝐼 𝐾 ← 𝐺 = - / 𝐼 𝐺 ← 𝐾 . The histograms suggest that / 𝐼 𝐺 ← 𝐾 &lt; 0 and / 𝐼 𝐾 ← 𝐺 &gt; 0 for the OISST data, whereas these information flows take both signs for the GFDL-CM4C192 data. These results are also found in Figs. 5b, 5e, 6b, and 6e.

Wenext discuss the entropy production rates. Equation (29) always holds. Thus, when 𝑐 𝐺 ← 𝐾 &gt; 0 and 𝑐 𝐾 ← 𝐺 &gt; 0, / 𝜎𝐾 and / 𝜎𝐺 have opposite signs to each other. For the GFDL-CM4C192 data, this opposite-sign relationship was confirmed in 1996 out of 2000 samples. On the other hand, for the OISST data, this relationship was confirmed in 1413 samples. The other 587 samples are in the

Fig. C1. Histograms of all estimated quantities obtained from the moving block bootstrap method with 2000 resampled time series for the OISST and GFDL-CM4C192 data. The quantities include all coefficients of the dynamical system, the entropy production rates, and the information flows. The vertical lines are at zero.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.01133v2-figure-0c8905b9-a876-4250-a19b-d336eac80032.jpg)

region where 𝑐 𝐺 ← 𝐾 &lt; 0 and 𝑐 𝐾 ← 𝐺 &gt; 0. In this region, both / 𝜎𝐺 and / 𝜎𝐾 are positive and have the same signs. This result is also found in Fig. 5, where the error bars cross the lines of 𝑐 𝐺 ← 𝐾 = 0.

In the region 𝑐 𝐺 ← 𝐾 &lt; 0, the interpretation of Maxwell's demon systems is not possible because of the same signs of / 𝜎𝐺 and / 𝜎𝐾 , but an interesting phenomenon still occurs in terms of information thermodynamics, which is called feedback cooling (Horowitz and Sandberg 2014; Loos and Klapp 2020). For the interpretation of Maxwell's demon systems, the signs of the entropy production rates need to be opposite to each other (Sections 2c and 2d). As for the Kuroshio Current, the entropy production rate and information flow are positive, as in the case of Maxwell's demon. Thus, even in the scenario of feedback cooling, there is no change in the interpretation for the Kuroshio Current.

The difference arises in the Gulf Stream, which has a positive entropy production rate ( / 𝜎𝐺 &gt; 0). The Gulf Stream still exhibits a negative information flow, which indicates that the Gulf Stream is interpreted as being controlled by the Kuroshio Current, like in Maxwell's demon. However, this control is done by negative feedback and reduces the effect of noise (Horowitz and Sandberg 2014; Loos and Klapp 2020). As a result, the SST magnitude of the Gulf Stream is reduced (Fig. 5c). This reduction is considered as being achieved by discarding heat to the environment, leading to the positive entropy production rate ( / 𝜎𝐺 &gt; 0). Even in this parameter region, the correlation between 𝑇𝐺 and 𝑇𝐾 is positive (Fig. 8), suggesting that the BCS occurs. However, the interpretation is slightly different. The synchronization is achieved by the Kuroshio Current, which reduces the SST fluctuations of the Gulf Stream caused by atmospheric and oceanic noise.

## APPENDIX D

## Constraint based on the second law of information thermodynamics

The analysis in Section 4 is based only on the signs of entropy production rates and information flows, without explicitly using the inequalities between them (i.e., the second law of information thermodynamics). Here, we derive an inequality for the instantaneous correlation coefficient using the second law of information thermodynamics, Eqs. (21) and (22). The derived inequality is a constraint that the dynamical system [Eqs. (16) and (17)] must satisfy.

The second law of information thermodynamics is expressed as the following two inequalities:

$$\frac { c _ { G \leftarrow K } } { D _ { G } } \left \langle T _ { K } \circ \frac { d T _ { G } } { d t } \right \rangle \geq \frac { \Sigma _ { G K } } { | \mathbf z | } \left \langle T _ { K } \circ \frac { d T _ { G } } { d t } \right \rangle,$$

$$\frac { c _ { K \leftarrow G } } { D _ { K } } \left \langle T _ { G } \circ \frac { d T _ { K } } { d t } \right \rangle \geq \frac { \overset {. } { \Sigma } _ { G K } } { | \mathbf z | } \left \langle T _ { G } \circ \frac { d T _ { K } } { d t } \right \rangle.$$

Equation (D1) is obtained by substituting Eqs. (24) and (26) into (21); likewise, Eq. (D2) is obtained by substituting Eqs. (25) and (27) into (22). From the steady-state condition [Eq. (B17)], either ⟨ 𝑇𝐾 ◦ d 𝑇𝐺 / d 𝑡 ⟩ or ⟨ 𝑇𝐺 ◦ d 𝑇𝐾 / d 𝑡 ⟩ is positive, while the other is negative. Here, we assume ⟨ 𝑇𝐾 ◦ d 𝑇𝐺 / d 𝑡 ⟩ &lt; 0. The above two inequalities are then combined into one:

$$\frac { c _ { G \leftarrow K } } { D _ { G } } \leq \frac { \Sigma _ { G K } } { \Sigma _ { K K } \Sigma _ { G G } - ( \Sigma _ { G K } ) ^ { 2 } } \leq \frac { c _ { K \leftarrow G } } { D _ { K } },$$

where we use | Σ | = Σ 𝐾𝐾 Σ 𝐺𝐺 -( Σ 𝐺𝐾 ) 2 . Including the other case where ⟨ 𝑇𝐾 ◦ d 𝑇𝐺 / d 𝑡 ⟩ &gt; 0, the final inequality becomes

$$\min \left \{ \frac { c _ { K \leftarrow G } } { D _ { K } }, \frac { c _ { G \leftarrow K } } { D _ { G } } \right \} \leq \frac { \Sigma _ { G K } } { \Sigma _ { K K } \Sigma _ { G G } - ( \Sigma _ { G K } ) ^ { 2 } } \leq \max \left \{ \frac { c _ { K \leftarrow G } } { D _ { K } }, \frac { c _ { G \leftarrow K } } { D _ { G } } \right \}.$$

When 𝑐 𝐾 ← / 𝐺 𝐷𝐾 and 𝑐 𝐺 ← / 𝐾 𝐷𝐺 are equal, it corresponds to the diagonal lines in Figs 5, 6, and 8. Crossing these lines reverses the magnitude relation between 𝑐 𝐾 ← / 𝐺 𝐷𝐾 and 𝑐 𝐺 ← / 𝐾 𝐷𝐺 . The signs of ⟨ 𝑇𝐾 ◦ d 𝑇𝐺 / d 𝑡 ⟩ and ⟨ 𝑇𝐺 ◦ d 𝑇𝐾 / d 𝑡 ⟩ are also reversed across these lines.

The obtained inequality is transformed into an inequality for the instantaneous correlation coefficient. We define the instantaneous correlation coefficient as follows:

$$\rho _ { 0 } \coloneqq \frac { \Sigma _ { G K } } { \sqrt { \Sigma _ { G G } \Sigma _ { K K } } }.$$

Assuming the time series are standardized, the variances Σ 𝐺𝐺 and Σ 𝐾𝐾 are equal to 1. The following inequality is then obtained:

$$\min \left \{ \frac { c _ { K \leftarrow G } } { D _ { K } }, \frac { c _ { G \leftarrow K } } { D _ { G } } \right \} \leq \frac { \rho _ { 0 } } { 1 - ( \rho _ { 0 } ) ^ { 2 } } \leq \max \left \{ \frac { c _ { K \leftarrow G } } { D _ { K } }, \frac { c _ { G \leftarrow K } } { D _ { G } } \right \},$$

where the function 𝑓 ( 𝑥 ) = 𝑥 /( 1 -𝑥 2 ) is monotonically increasing in the interval (-1 1 , that is, , ) 𝑓 ( 𝑥 ) is invertible. The inequality (D6) means that 𝑓 ( 𝜌 0 ) is bounded by 𝑐 𝐾 ← / 𝐺 𝐷𝐾 and 𝑐 𝐺 ← / 𝐾 𝐷𝐺 , suggesting that 𝜌 0 increases as the interaction coefficients become larger or the noise amplitudes become smaller.

The inequality (D6) is a condition that the correlation coefficient 𝜌 0 must satisfy; thus, this inequality can be regarded as a constraint on the dynamical system. It is possible to obtain the explicit form of 𝜌 0 . However, this explicit form requires the use of the complex formulas (B8), (B9), and (B10). The second law of information thermodynamics provides the inequality (D6) without relying on these complex formulas. In other words, the inequality implies the dependence of 𝜌 0 on the parameters, such as 𝑐 𝐺 ← 𝐾 and 𝐷𝐺 , without knowing the explicit form of 𝜌 0 .

## References

Allahverdyan, A. E., D. Janzing, and G. Mahler, 2009: Thermodynamic efficiency of information and heat flow. Journal of Statistical Mechanics: Theory and Experiment , 2009 (09) , P09011, https://doi.org/10.1088/1742-5468/2009/09/P09011, URL https://dx.doi.org/10.1088/ 1742-5468/2009/09/P09011.

Benzi, R., 2010: Stochastic resonance: from climate to biology. Nonlinear Processes in Geophysics , 17 (5) , 431-441, https://doi.org/10.5194/npg-17-431-2010, URL https://npg. copernicus.org/articles/17/431/2010/.

Bouchet, F., and A. Venaille, 2012: Statistical mechanics of two-dimensional and geophysical flows. Physics Reports , 515 (5) , 227-295, https://doi.org/https://doi.org/10.1016/j.physrep. 2012.02.001, URL https://www.sciencedirect.com/science/article/pii/S0370157312000518, statistical mechanics of two-dimensional and geophysical flows.

Brockwell, P. J., and R. A. Davis, 1991: Time series: theory and methods . Springer science &amp; business media.

Campa, A., T. Dauxois, D. Fanelli, and S. Ruffo, 2014: Physics of Long-Range Interacting Systems . Oxford University Press.

Cover, T. M., and J. A. Thomas, 2005: Elements of Information Theory . Second edition ed., John Wiley &amp; Sons, Ltd.

Dijkstra, H. A., 2013: Nonlinear climate dynamics . Cambridge University Press.

Eyring, V., S. Bony, G. A. Meehl, C. A. Senior, B. Stevens, R. J. Stouffer, and K. E. Taylor, 2016: Overview of the coupled model intercomparison project phase 6 (cmip6) experimental design and organization. Geoscientific Model Development , 9 (5) , 1937-1958, https://doi.org/ 10.5194/gmd-9-1937-2016, URL https://gmd.copernicus.org/articles/9/1937/2016/.

Franzke, C. L. E., T. J. O'Kane, J. Berner, P. D. Williams, and V. Lucarini, 2015: Stochastic climate theory and modeling. WIREs Climate Change , 6 (1) , 63-78, https://doi.org/https:// doi.org/10.1002/wcc.318, URL https://wires.onlinelibrary.wiley.com/doi/abs/10.1002/wcc.318, https://wires.onlinelibrary.wiley.com/doi/pdf/10.1002/wcc.318.

Gallego, B., and P. Cessi, 2000: Exchange of heat and momentum between the atmosphere and the ocean: A minimal model of decadal oscillations. Climate Dynamics , 16 , 479-489.

Gallego, B., and P. Cessi, 2001: Decadal variability of two oceans and an atmosphere. Journal of Climate , 14 (13) , 2815 -2832, https://doi.org/10.1175/1520-0442(2001)014 2815: ⟨ DVOTOA 2.0.CO;2, ⟩ URL https://journals.ametsoc.org/view/journals/clim/14/13/1520-0442 2001 014 2815 dvotoa 2.0.co 2.xml.

Gardiner, C., 2009: Stochastic Methods . Springer Series in Synergetics, Springer Berlin, Heidelberg.

Hamilton, J. D., 1994: Time series analysis . Princeton university press.

Hasselmann, K., 1976: Stochastic climate models part i. theory. Tellus , 28 (6) , 473-485, https://doi.org/https://doi.org/10.1111/j.2153-3490.1976.tb00696.x, URL https://onlinelibrary. wiley.com/doi/abs/10.1111/j.2153-3490.1976.tb00696.x, https://onlinelibrary.wiley.com/doi/ pdf/10.1111/j.2153-3490.1976.tb00696.x.

Held, I. M., and Coauthors, 2019: Structure and performance of gfdl's cm4.0 climate model. Journal of Advances in Modeling Earth Systems , 11 (11) , 3691-3727, https://doi.org/https://doi. org/10.1029/2019MS001829, URL https://agupubs.onlinelibrary.wiley.com/doi/abs/10.1029/ 2019MS001829, https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2019MS001829.

Hogg, A. M. C., W. K. Dewar, P. D. Killworth, and J. R. Blundell, 2006: Decadal variability of the midlatitude climate system driven by the ocean circulation. Journal of Climate , 19 (7) , 1149 - 1166, https://doi.org/10.1175/JCLI3651.1, URL https://journals.ametsoc.org/view/journals/ clim/19/7/jcli3651.1.xml.

Horowitz, J. M., and M. Esposito, 2014: Thermodynamics with continuous information flow. Phys. Rev. X , 4 , 031 015, https://doi.org/10.1103/PhysRevX.4.031015, URL https://link.aps.org/doi/ 10.1103/PhysRevX.4.031015.

Horowitz, J. M., and H. Sandberg, 2014: Second-law-like inequalities with information and their interpretations. New Journal of Physics , 16 (12) , 125 007, https://doi.org/10.1088/1367-2630/ 16/12/125007, URL https://dx.doi.org/10.1088/1367-2630/16/12/125007.

Hoskins, B. J., and P. J. Valdes, 1990: On the existence of storm-tracks. Journal of Atmospheric Sciences , 47 (15) , 1854 - 1864, https://doi.org/10.1175/1520-0469(1990)047 1854:OTEOST 2.0. ⟨ ⟩ CO;2, URL https://journals.ametsoc.org/view/journals/atsc/47/15/1520-0469 1990 047 1854 oteost 2 0 co 2.xml.

Hou, C. K. J., and K. Behdinan, 2022: Dimensionality reduction in surrogate modeling: A review of combined methods. Data Science and Engineering , 7 (4) , 402-427.

Ito, S., 2016a: Backward transfer entropy: Informational measure for detecting hidden markov models and its interpretations in thermodynamics, gambling and causality. Scientific reports , 6 (1) , 36 831.

Ito, S., 2016b: Information thermodynamics on causal networks and its application to biochemical signal transduction . Springer.

Ito, S., and T. Sagawa, 2015: Maxwell's demon in biochemical signal transduction with feedback loop. Nature communications , 6 (1) , 1-6.

Kelly, K. A., and S. Dong, 2004: The relationship of western boundary current heat transport and storage to midlatitude ocean-atmosphere interaction. Earth's Climate , American Geophysical Union (AGU), 347-363, https://doi.org/https://doi.org/10.1029/ 147GM19, URL https://agupubs.onlinelibrary.wiley.com/doi/abs/10.1029/147GM19, https:// agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/147GM19.

Kelly, K. A., R. J. Small, R. M. Samelson, B. Qiu, T. M. Joyce, Y.-O. Kwon, and M. F. Cronin, 2010: Western boundary currents and frontal air-sea interaction: Gulf stream and kuroshio extension. Journal of Climate , 23 (21) , 5644 - 5667, https://doi.org/10.1175/2010JCLI3346.1, URL https://journals.ametsoc.org/view/journals/clim/23/21/2010jcli3346.1.xml.

Kohyama, T., Y. Yamagami, H. Miura, S. Kido, H. Tatebe, and M. Watanabe, 2021: The gulf stream and kuroshio current are synchronized. Science , 374 (6565) , 341-346, https://doi.org/10. 1126/science.abh3295, URL https://www.science.org/doi/abs/10.1126/science.abh3295, https: //www.science.org/doi/pdf/10.1126/science.abh3295.

Koski, J. V., A. Kutvonen, I. M. Khaymovich, T. Ala-Nissila, and J. P. Pekola, 2015: On-chip maxwell's demon as an information-powered refrigerator. Phys. Rev. Lett. , 115 ,

260 602, https://doi.org/10.1103/PhysRevLett.115.260602, URL https://link.aps.org/doi/10. 1103/PhysRevLett.115.260602.

Koski, J. V., V. F. Maisi, J. P. Pekola, and D. V. Averin, 2014: Experimental realization of a szilard engine with a single electron. Proceedings of the National Academy of Sciences , 111 (38) , 13 786-13 789, https://doi.org/10.1073/pnas.1406966111, URL https://www.pnas.org/doi/abs/ 10.1073/pnas.1406966111, https://www.pnas.org/doi/pdf/10.1073/pnas.1406966111.

Kuwano-Yoshida, A., and S. Minobe, 2017: Storm-track response to sst fronts in the northwestern pacific region in an agcm. Journal of Climate , 30 (3) , 1081 - 1102, https://doi.org/10.1175/ JCLI-D-16-0331.1, URL https://journals.ametsoc.org/view/journals/clim/30/3/jcli-d-16-0331. 1.xml.

Kwon, Y.-O., M. A. Alexander, N. A. Bond, C. Frankignoul, H. Nakamura, B. Qiu, and L. A. Thompson, 2010: Role of the gulf stream and kuroshio-oyashio systems in largescale atmosphere-ocean interaction: A review. Journal of Climate , 23 (12) , 3249 - 3281, https://doi.org/10.1175/2010JCLI3343.1, URL https://journals.ametsoc.org/view/journals/clim/ 23/12/2010jcli3343.1.xml.

Latif, M., and T. P. Barnett, 1994: Causes of decadal climate variability over the north pacific and north america. Science , 266 (5185) , 634-637, https://doi.org/10.1126/science.266.5185.634, URL https://www.science.org/doi/abs/10.1126/science.266.5185.634, https://www.science.org/ doi/pdf/10.1126/science.266.5185.634.

Loos, S. A. M., and S. H. L. Klapp, 2020: Irreversibility, heat and information flows induced by non-reciprocal interactions. New Journal of Physics , 22 (12) , 123 051, https://doi.org/10.1088/ 1367-2630/abcc1e, URL https://dx.doi.org/10.1088/1367-2630/abcc1e.

Ma, X., and Coauthors, 2016: Western boundary currents regulated by interaction between ocean eddies and the atmosphere. Nature , 535 (7613) , 533-537.

Majda, A., and X. Wang, 2006: Nonlinear Dynamics and Statistical Theories for Basic Geophysical Flows . Cambridge University Press.

Michael Sherman, F. M. S. J., and F. M. Speed, 1998: Analysis of tidal data via the blockwise bootstrap. Journal of Applied Statistics , 25 (3) , 333-340, https://doi.org/10.1080/02664769823061, URL https://doi.org/10.1080/02664769823061, https://doi.org/10.1080/02664769823061.

Minobe, S., A. Kuwano-Yoshida, N. Komori, S.-P. Xie, and R. J. Small, 2008: Influence of the gulf stream on the troposphere. Nature , 452 (7184) , 206-209.

Mudelsee, M., 2014: Climate Time Series Analysis . Springer Cham.

Ogawa, F., H. Nakamura, K. Nishii, T. Miyasaka, and A. Kuwano-Yoshida, 2012: Dependence of the climatological axial latitudes of the tropospheric westerlies and storm tracks on the latitude of an extratropical oceanic front. Geophysical Research Letters , 39 (5) , https://doi.org/https://doi. org/10.1029/2011GL049922, URL https://agupubs.onlinelibrary.wiley.com/doi/abs/10.1029/ 2011GL049922, https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/2011GL049922.

Omrani, N.-E., F. Ogawa, H. Nakamura, N. Keenlyside, S. W. Lubis, and K. Matthes, 2019: Key role of the ocean western boundary currents in shaping the northern hemisphere climate. Scientific reports , 9 (1) , 3014.

Parrondo, J. M., J. M. Horowitz, and T. Sagawa, 2015: Thermodynamics of information. Nature physics , 11 (2) , 131-139.

Peliti, L., and S. Pigolotti, 2021: Stochastic Thermodynamics: An Introduction . Princeton University Press.

Qiu, B., and S. Chen, 2005: Variability of the kuroshio extension jet, recirculation gyre, and mesoscale eddies on decadal time scales. Journal of Physical Oceanography , 35 (11) , 2090 - 2103, https://doi.org/10.1175/JPO2807.1, URL https://journals.ametsoc.org/view/journals/ phoc/35/11/jpo2807.1.xml.

Qiu, B., S. Chen, N. Schneider, and B. Taguchi, 2014: A coupled decadal prediction of the dynamic state of the kuroshio extension system. Journal of Climate , 27 (4) , 1751 - 1764, https://doi.org/10.1175/JCLI-D-13-00318.1, URL https://journals.ametsoc.org/view/journals/ clim/27/4/jcli-d-13-00318.1.xml.

Reddy, G. T., M. P. K. Reddy, K. Lakshmanna, R. Kaluri, D. S. Rajput, G. Srivastava, and T. Baker, 2020: Analysis of dimensionality reduction techniques on big data. IEEE Access , 8 , 54 776-54 788, https://doi.org/10.1109/ACCESS.2020.2980942.

Reynolds, R. W., T. M. Smith, C. Liu, D. B. Chelton, K. S. Casey, and M. G. Schlax, 2007: Daily high-resolution-blended analyses for sea surface temperature. Journal of Climate , 20 (22) , 5473 - 5496, https://doi.org/10.1175/2007JCLI1824.1, URL https://journals.ametsoc.org/view/ journals/clim/20/22/2007jcli1824.1.xml.

Sagawa, T., 2019: Second law, entropy production, and reversibility in thermodynamics of information. Energy Limits in Computation: A Review of Landauer's Principle, Theory and Experiments , Springer International Publishing, Cham, 101-139, https://doi.org/ 10.1007/978-3-319-93458-7 3, URL https://doi.org/10.1007/978-3-319-93458-7 3.

Sagawa, T., and M. Ueda, 2009: Minimal energy cost for thermodynamic information processing: Measurement and information erasure. Phys. Rev. Lett. , 102 , 250 602, https://doi.org/10.1103/ PhysRevLett.102.250602, URL https://link.aps.org/doi/10.1103/PhysRevLett.102.250602, erratum: Phys. Rev. Lett. 106, 189901.

Sagawa, T., and M. Ueda, 2013a: Information thermodynamics: Maxwell's demon in nonequilibrium dynamics. Nonequilibrium Statistical Physics of Small Systems , John Wiley &amp; Sons, Ltd, chap. 6, 181-211, https://doi.org/https://doi.org/10.1002/9783527658701.ch6, URL https:// onlinelibrary.wiley.com/doi/abs/10.1002/9783527658701.ch6, https://onlinelibrary.wiley.com/ doi/pdf/10.1002/9783527658701.ch6.

Sagawa, T., and M. Ueda, 2013b: Role of mutual information in entropy production under information exchanges. New Journal of Physics , 15 (12) , 125 012, https://doi.org/10.1088/1367-2630/ 15/12/125012, URL https://dx.doi.org/10.1088/1367-2630/15/12/125012.

Sanders, F., 1986: Explosive cyclogenesis in the west-central north atlantic ocean, 1981-84. part i: Composite structure and mean behavior. Monthly Weather Review , 114 (10) , 1781 - 1794, https://doi.org/10.1175/1520-0493(1986)114 1781:ECITWC 2.0.CO;2, ⟨ ⟩ URL https://journals. ametsoc.org/view/journals/mwre/114/10/1520-0493 1986 114 1781 ecitwc 2 0 co 2.xml.

Seifert, U., 2005: Entropy production along a stochastic trajectory and an integral fluctuation theorem. Phys. Rev. Lett. , 95 , 040 602, https://doi.org/10.1103/PhysRevLett.95.040602, URL https://link.aps.org/doi/10.1103/PhysRevLett.95.040602.

Seifert, U., 2012: Stochastic thermodynamics, fluctuation theorems and molecular machines. Reports on Progress in Physics , 75 (12) , 126 001, https://doi.org/10.1088/0034-4885/75/12/ 126001, URL https://dx.doi.org/10.1088/0034-4885/75/12/126001.

Sekimoto, K., 2010: Stochastic Energetics . Lecture Notes in Physics, Springer Berlin, Heidelberg.

Shiraishi, N., 2023: An Introduction to Stochastic Thermodynamics: From Basic to Advanced , Vol. 212. Springer Nature.

Stommel, H., 1948: The westward intensification of wind-driven ocean currents. Eos, Transactions American Geophysical Union , 29 (2) , 202-206, https://doi.org/ https://doi.org/10.1029/TR029i002p00202, URL https://agupubs.onlinelibrary.wiley.com/ doi/abs/10.1029/TR029i002p00202, https://agupubs.onlinelibrary.wiley.com/doi/pdf/10.1029/ TR029i002p00202.

Strogatz, S. H., 2018: Nonlinear dynamics and chaos: with applications to physics, biology, chemistry, and engineering . Second edition ed., CRC press.

Toyabe, S., T. Sagawa, M. Ueda, E. Muneyuki, and M. Sano, 2010: Experimental demonstration of information-to-energy conversion and validation of the generalized jarzynski equality. Nature physics , 6 (12) , 988-992.

Venaille, A., and F. Bouchet, 2011: Oceanic rings and jets as statistical equilibrium states. Journal of Physical Oceanography , 41 (10) , 1860 - 1873, https://doi.org/10.1175/2011JPO4583.1, URL https://journals.ametsoc.org/view/journals/phoc/41/10/2011jpo4583.1.xml.

Yamagami, Y., H. Tatebe, T. Kohyama, S. Kido, and S. Okajima, 2024: Causality Between the Kuroshio and Gulf Stream SST Variability. Asia Oceania Geosciences Society (AOGS) 21st Annual Meeting , PyeongChang, Korea, Asia Oceania Geosci. Soc., AS35-A013.
Yuki Yasuda Tsubasa Kohyama
A Mechanism of Stochastic Synchronization in the Climate System: an Interpretation of the Boundary Current Synchronization as a Maxwell's Demon
This study has applied information thermodynamics to a bivariate linear
stochastic differential equation (SDE) that describes a synchronization
phenomenon of sea surface temperatures (SSTs) between the Gulf Stream and the
Kuroshio Current, which is referred to as the boundary current synchronization
(BCS). Information thermodynamics divides the entire system fluctuating with
stochastic noise into subsystems and describes the interactions between these
subsystems from the perspective of information transfer. The SDE coefficients
have been estimated through regression analysis using observational and
numerical simulation data. In the absence of stochastic noise, the solution of
the estimated SDE shows that the SSTs relax toward zero without oscillation.
The estimated SDE can be interpreted as a Maxwell's demon system, with the Gulf
Stream playing the role of the "Particle" and the Kuroshio Current playing the
role of the "Demon." The Gulf Stream forces the SST of the Kuroshio Current to
be in phase. By contrast, the Kuroshio Current maintains the phase by
interfering with the relaxation of the Gulf Stream SST. In the framework of
Maxwell's demon, the Gulf Stream is measured by the Kuroshio Current, whereas
the Kuroshio Current performs feedback control on the Gulf Stream. When the
Gulf Stream and the Kuroshio Current are coupled in an appropriate parameter
regime, synchronization is realized with atmospheric and oceanic fluctuations
as the driving source. This new mechanism, "stochastic synchronization,"
suggests that such fluctuations can be converted into directional variations in
a subsystem of the climate system through synchronization by utilizing
information (i.e., information-to-energy conversion).