We investigate the theoretical and numerical computation of rare transitions in simple geophysical turbulent models. We consider the barotropic quasi-geostrophic and two-dimensional Navier-Stokes equations in regimes where bistability between two coexisting large-scale attractors exist. By means of large deviations and instanton theory with the use of an Onsager-Machlup path integral formalism for the transition probability, we show how one can directly compute the most probable transition path between two coexisting attractors analytically in an equilibrium (Langevin) framework and numerically otherwise. We adapt a class of numerical optimization algorithms known as minimum action methods to simple geophysical turbulent models. We show that by numerically minimizing an appropriate action functional in a large deviation limit, one can predict the most likely transition path for a rare transition between two states. By considering examples where theoretical predictions can be made, we show that the minimum action method successfully predicts the most likely transition path. Finally, we discuss the application and extension of such numerical optimization schemes to the computation of rare transitions observed in direct numerical simulations and experiments and to other, more complex, turbulent systems. models. Stochastic resonance, however, remains a very interesting possibility. To address such issues one should study the attractors and the dynamics of the rare transition between attractors in a hierarchy of models from the simplest to more complex ones used by the climate community. For these complex climate models, which genuinely reproduce the turbulent nature of the Earthʼs atmosphere and ocean dynamics, such a task is currently inconceivable and seems unreachable in the foreseeable future using direct numerical simulations. The reasons are the rarity of the transitions and the computational complexity of these models. The main aim of this paper is to make a step in the direction of this challenge by studying bistability and the associated transitions in turbulent dynamics using tools that will allow one to compute transitions in more complex systems in the near future.These rare transitions are essential phenomena because they correspond to drastic changes in complex system behavior. Moreover, they cannot be studied using conventional tools. They contain dynamics occurring on multiple and extremely different timescales, usually with no spectral gap. This prevents the use of classical tools from dynamical system theory. The theoretical understanding of these transitions is an extremely difficult problem due to the complexity, the large number of degrees of freedom, and the non-equilibrium nature of many of these flows. Up to now, there have been an extremely limited number of theoretical results, where analysis has been limited to analogies with models of very few degrees of freedom [27] or to specific classes of systems that can be directly related to equilibrium Langevin dynamics [28]. For this reason, the u...