"A first goal of the project was to develop rare event algorithms to make them suitable for studying complex systems such as turbulent flows or the climate system. The idea of one of these algorithms, called ""adaptive multilevel splitting"" (AMS), is to bias the statistics of an ensemble of simulations in a controlled manner, by progressively replacing the least performing trajectories (with respect to the distance to the rare event of interest) with mutations of the better performing ones. In the standard formulation of this algorithm, each simulation must be integrated until it hits one of two sets in phase space, corresponding either to typical conditions or to a rare event of interest. This is cumbersome for high-dimensional complex systems. We have given a different formulation of the algorithm, where the trajectories are always integrated for a fixed duration. It amounts at evaluating the probability that a given event occurs within a given time. This is well suited to many problems in climate science. We have also shown that the algorithm allows for evaluating ""return times"" for rare event, i.e. the typical time between two occurrences. This metric is commonly used in applications, for heat waves, floods, etc. As an example of application to complex systems, we have computed return times for extreme drag force acting on an object immersed in a turbulent flow. These results led to a journal article in J. Stat. Mech, published in 2018, and have been presented in conferences and various seminars.
Another limitation is the need to construct the distance to the rare event, called ""score function"", which is the crucial object in the selection step of the AMS algorithm. We have outlined a possible strategy, using machine learning methods to iteratively improve the score function. As a first step, we have compared direct methods and machine learning to estimate ""committor functions"", the mathematical object computed by the AMS. These results have been presented at a conference and shall lead to a journal publication soon.
All our codes have been published in a Python package, and a software metapaper describing it shall be submitted soon.
The second main goal was to study the possibility of abrupt transitions in atmospheric jets. We have identified a particularly interesting candidate for such phenomena: ""equatorial superrotation"", which corresponds to the appearance of a strong eastward jet at the equator. While tropical surface winds are easterly on current Earth, superrotation is observed on many planetary atmospheres, and might also be relevant for climates of the past. It might also be speculated that a transition to superrotation might occur due to anthropogenic climate change, providing another example of tipping point for global climate. While dynamical mechanisms leading to superrotation had been studied by dynamicists, the nature of the transition, abrupt or continuous, remained virtually unexplored. We have shown unambiguously, based on idealized models (an analytical model and a 2D numerical model) that a wave-jet feedback mechanism, sketched a few years ago by other researchers, indeed led to bistability and abrupt transitions. These results have been presented in several conferences and published in J. Atmos. Sci.
We have carried out numerical simulations with a full 3D General Circulation Model of the Atmosphere; while more work is still needed to understand the role of different parameters, this preliminary work seems to indicate that our findings still hold in a more realistic framework."