In this study, we perform experiments with a coupled atmosphere-ocean general circulation model (CGCM) to examine ENSO’s influence on the interannual sea-surface temperature (SST) variability of the tropical Indian Ocean. The control experiment includes both the Indian and Pacific Oceans in the ocean model component of the CGCM (the Indo-Pacific Run). The anomaly experiment excludes ENSO’s influence by including only the Indian Ocean while prescribing monthly-varying climatological SSTs for the Pacific Ocean (the Indian-Ocean Run). In the Indo-Pacific Run, an oscillatory mode of the Indian Ocean SST variability is identified by a multi-channel singular spectral analysis (MSSA). The oscillatory mode comprises two patterns that can be identified with the Indian Ocean Zonal Mode (IOZM) and a basin-wide warming/cooling mode respectively. In the model, the IOZM peaks about 3–5 months after ENSO reaches its maximum intensity. The basin mode peaks 8 months after the IOZM. The timing and associated SST patterns suggests that the IOZM is related to ENSO, and the basin-wide warming/cooling develops as a result of the decay of the IOZM spreading SST anomalies from western Indian Ocean to the eastern Indian Ocean. In contrast, in the Indian-Ocean Run, no oscillatory modes can be identified by the MSSA, even though the Indian Ocean SST variability is characterized by east–west SST contrast patterns similar to the IOZM. In both control and anomaly runs, IOZM-like SST variability appears to be associated with forcings from fluctuations of the Indian monsoon. Our modeling results suggest that the oscillatory feature of the IOZM is primarily forced by ENSO.

Systems of correlated particles appear in many fields of modern science and represent some of the most intractable computational problems in nature. The computational challenge in these systems arises when interactions become comparable to other energy scales, which makes the state of each particle depend on all other particles¹. The lack of general solutions for the three-body problem and acceptable theory for strongly correlated electrons shows that our understanding of correlated systems fades when the particle number or the interaction strength increases. One of the hallmarks of interacting systems is the formation of multiparticle bound states^2-9. Here we develop a high-fidelity parameterizable fSim gate and implement the periodic quantum circuit of the spin-½ XXZ model in a ring of 24 superconducting qubits. We study the propagation of these excitations and observe their bound nature for up to five photons. We devise a phase-sensitive method for constructing the few-body spectrum of the bound states and extract their pseudo-charge by introducing a synthetic flux. By introducing interactions between the ring and additional qubits, we observe an unexpected resilience of the bound states to integrability breaking. This finding goes against the idea that bound states in non-integrable systems are unstable when their energies overlap with the continuum spectrum. Our work provides experimental evidence for bound states of interacting photons and discovers their stability beyond the integrability limit.