Optimally coherent sets in geophysical flows: A new approach to delimiting the stratospheric polar vortex

The "edge" of the Antarctic polar vortex is known to behave as a barrier to the meridional (poleward) transport of ozone during the austral winter. This chemical isolation of the polar vortex from the middle and low latitudes produces an ozone minimum in the vortex region, intensifying the ozone hole relative to that which would be produced by photochemical processes alone. Observational determination of the vortex edge remains an active field of research. In this letter, we obtain objective estimates of the structure of the polar vortex by introducing a new technique based on transfer operators that aims to find regions with minimal external transport. Applying this new technique to European Centre for Medium-Range Weather Forecasts (ECMWF) ERA-40 three-dimensional velocity data we produce an improved three-dimensional estimate of the vortex location in the upper stratosphere where the vortex is most pronounced. This novel computational approach has wide potential application in detecting and analysing mixing structures in a variety of atmospheric, oceanographic, and general fluid dynamical settings

The backbone of the climate network

We propose a method to reconstruct and analyze a complex network from data generated by a spatio-temporal dynamical system, relying on the nonlinear mutual information of time series analysis and betweenness centrality of complex network theory. We show, that this approach reveals a rich internal structure in complex climate networks constructed from reanalysis and model surface air temperature data. Our novel method uncovers peculiar wave-like structures of high energy flow, that we relate to global surface ocean currents. This points to a major role of the oceanic surface circulation in coupling and stabilizing the global temperature field in the long term mean (140 years for the model run and 60 years for reanalysis data). We find that these results cannot be obtained using classical linear methods of multivariate data analysis, and have ensured their robustness by intensive significance testing.Comment: 6 pages, 5 figure

Influence of turbulent advection on a phytoplankton ecosystem with nonuniform carrying capacity

In this work we study a plankton ecosystem model in a turbulent flow. The plankton model we consider contains logistic growth with a spatially varying background carrying capacity and the flow dynamics are generated using the two-dimensional (2D) Navier-Stokes equations. We characterize the system in terms of a dimensionless parameter, γ TB / TF, which is the ratio of the ecosystem biological time scales TB and the flow time scales TF. We integrate this system numerically for different values of γ until the mean plankton reaches a statistically stationary state and examine how the steady-state mean and variance of plankton depends on γ. Overall we find that advection in the presence of a nonuniform background carrying capacity can lead to very different plankton distributions depending on the time scale ratio γ. For small γ the plankton distribution is very similar to the background carrying capacity field and has a mean concentration close to the mean carrying capacity. As γ increases the plankton concentration is more influenced by the advection processes. In the largest γ cases there is a homogenization of the plankton concentration and the mean plankton concentration approaches the harmonic mean, 1/K -1. We derive asymptotic approximations for the cases of small and large γ. We also look at the dependence of the power spectra exponent, β, on γ where the power spectrum of plankton is k-β. We find that the power spectra exponent closely obeys β=1+2/γ as predicted by earlier studies using simple models of chaotic advection

Inferring directed climatic interactions with renormalized partial directed coherence and directed partial correlation

ACKNOWLEDGMENTS This work was supported in part by Spanish MINECO/FEDER (FIS2015-66503-C3-2-P) and ITN LINC (FP7 289447). C.M. also acknowledges partial support from ICREA ACADEMIAPeer reviewedPublisher PD

Dispersive stabilization of the inverse cascade for the Kolmogorov flow

It is shown by perturbation techniques and numerical simulations that the inverse cascade of kink-antikink annihilations, characteristic of the Kolmogorov flow in the slightly supercritical Reynolds number regime, is halted by the dispersive action of Rossby waves in the beta-plane approximation. For beta tending to zero, the largest excited scale is proportional to the logarithm of one over beta and differs strongly from what is predicted by standard dimensional phenomenology which ignores depletion of nonlinearity.Comment: 4 pages, LATEX, 3 figures. v3: revised version with minor correction

Generation and Structure of Solitary Rossby Vortices in Rotating Fluids

The formation of zonal flows and vortices in the generalized Charney-Hasegawa-Mima equation is studied. We focus on the regime when the size of structures is comparable to or larger than the deformation (Rossby) radius. Numerical simulations show the formation of anticyclonic vortices in unstable shear flows and ring-like vortices with quiescent cores and vorticity concentrated in a ring. Physical mechanisms that lead to these phenomena and their relevance to turbulence in planetary atmospheres are discussed.Comment: 3 pages in REVTeX, 5 postscript figures separately, submitted to Phys. Rev.

Hetonic quartets in a two-layer quasi-geostrophic flow : V-states and stability

M.A.S. and X.C. were supported by RFBR/CNRS (PRC Grant No. 16-55-150001/1069). M.A.S. was supported also by RFBR (Grant No. 16-05-00121), RSF (Grant No. 14-50-00095, geophysical applications) and MESRF (Grant No. 14.W.03.31.0006, numerical simulation, vortex dynamics).We investigate families of finite core vortex quartets in mutual equilibrium in a two- layer quasi-geostrophic flow. The finite core solutions stem from known solutions for discrete (singular) vortex quartets. Two vortices lie in the top layer and two vortices lie in the bottom layer. Two vortices have a positive potential vorticity anomaly while the two others have negative potential vorticity anomaly. The vortex configurations are therefore related to the baroclinic dipoles known in the literature as hetons. Two main branches of solutions exist depending on the arrangement of the vortices: the translating zigzag-shaped hetonic quartets and the rotating zigzag- shaped hetonic quartets. By addressing their linear stability, we show that while the rotating quartets can be unstable over a large range of the parameter space, most translating quartets are stable. This has implications on the longevity of such vortex equilibria in the oceans.PostprintPeer reviewe