# Multiscale Modeling

The emergence of organized multiscale patterns resulting from convection is ubiquitous, observed throughout different cloud types. The reproduction of such patterns by general circulation models remains a challenge due to the complex nature of clouds, characterized by processes interacting over a wide range of spatio-temporal scales. The new advances in data-driven modeling techniques have raised a lot of promises to discover dynamical equations from partial observations of complex systems.

This study presents such a discovery from high-resolution satellite datasets of continental cloud fields. The model is made of stochastic differential equations able to simulate with high fidelity the spatio-temporal coherence and variability of the cloud patterns such as the characteristic lifetime of individual clouds or global organizational features governed by convective inertia gravity waves. This feat is achieved through the model's lagged effects associated with convection recirculation times, and hidden variables parameterizing the unobserved processes and variables.

We parameterize sub-grid scale (SGS) fluxes in sinusoidally forced two-dimensional turbulence on the *β*-plane at high Reynolds numbers (Re ∼25,000) using simple 2-layer convolutional neural networks (CNN) having only O(1000) parameters, two orders of magnitude smaller than recent studies employing deeper CNNs with 8–10 layers; we obtain stable, accurate, and long-term online or a posteriori solutions at 16× downscaling factors. Our methodology significantly improves training efficiency and speed of online large eddy simulations runs, while offering insights into the physics of closure in such turbulent flows. Our approach benefits from extensive hyperparameter searching in learning rate and weight decay coefficient space, as well as the use of cyclical learning rate annealing, which leads to more robust and accurate online solutions compared to fixed learning rates. Our CNNs use either the coarse velocity or the vorticity and strain fields as inputs, and output the two components of the deviatoric stress tensor, *S*_{d}. We minimize a loss between the SGS vorticity flux divergence (computed from the high-resolution solver) and that obtained from the CNN-modeled *S*_{d}, without requiring energy or enstrophy preserving constraints. The success of shallow CNNs in accurately parameterizing this class of turbulent flows implies that the SGS stresses have a weak non-local dependence on coarse fields; it also aligns with our physical conception that small-scales are locally controlled by larger scales such as vortices and their strained filaments. Furthermore, 2-layer CNN-parameterizations are more likely to be interpretable.

*modus operandi*, the backward integration time is the key parameter to select per scale/variable to parameterize in order to derive the relevant parameterizations which are doomed to be no longer exact away from instability onset due to the breakdown of slaving typically encountered, e.g., for chaotic regimes. The selection criterion is then made through data-informed minimization of a least-square parameterization defect. It is thus shown through optimization of the backward integration time per scale/variable to parameterize, that skilled OPM reduced systems can be derived for predicting with accuracy higher-order critical transitions or catastrophic tipping phenomena, while training our parameterization formulas for regimes prior to these transitions takes place.

A central challenge in physics is to describe non-equilibrium systems driven by randomness, such as a randomly growing interface, or fluids subject to random fluctuations that account e.g. for local stresses and heat fluxes in the fluid which are not related to the velocity and temperature gradients. For deterministic systems with infinitely many degrees of freedom, normal form and center manifold theory have shown a prodigious efficiency to often completely characterize how the onset of linear instability translates into the emergence of nonlinear patterns, associated with genuine physical regimes. However, in presence of random fluctuations, the underlying reduction principle to the center manifold is seriously challenged due to large excursions caused by the noise, and the approach needs to be revisited.

In this study, we present an alternative framework to cope with these difficulties exploiting the approximation theory of stochastic invariant manifolds, on one hand, and energy estimates measuring the defect of parameterization of the high-modes, on the other. To operate for fluid problems subject to stochastic stirring forces, these error estimates are derived under assumptions regarding dissipation effects brought by the high-modes in order to suitably counterbalance the loss of regularity due to the nonlinear terms. As a result, the approach enables us to analyze, from reduced equations of the stochastic fluid problem, the occurrence in large probability of a stochastic analogue to the pitchfork bifurcation, as long as the noise’s intensity and the eigenvalue’s magnitude of the mildly unstable mode scale accordingly.

In the case of SPDEs forced by a multiplicative noise in the orthogonal subspace of e.g. its mildly unstable mode, our parameterization formulas show that the noise gets transmitted to this mode via non-Markovian coefficients, and that the reduced equation is only stochastically driven by the latter. These coefficients depend explicitly on the noise path's history, and their memory content is self-consistently determined by the intensity of the random force and its interaction through the SPDE's nonlinear terms. Applications to a stochastic Rayleigh-B\'enard problem are detailed, for which conditions for a stochastic pitchfork bifurcation (in large probability) to occur, are clarified.

Nonlinear time delay systems produce inherently delay-induced periodic oscillations, which are, however, too idealistic compared to observations. We exhibit a unified stochastic framework to systematically rectify such oscillations into oscillatory patterns with enriched temporal variabilities through generic, nonlinear responses to stochastic perturbations. Two paradigms of noise-driven chaos in high dimension are identified, fundamentally different from chaos triggered by parameter-space noise. Noteworthy is a low-dimensional stretch-and-fold mechanism, leading to stochastic strange attractors exhibiting horseshoe-like structures mirroring turbulent transport of passive tracers. The other is high-dimensional , with noise acting along the critical eigendirection and transmitted to “deeper” stable modes through nonlinearity, leading to stochastic attractors exhibiting swarm-like behaviors with power-law and scale break properties. The theory is applied to cloud delay models to parameterize missing physics such as intermittent rain and Lagrangian turbulent effects. The stochastically rectified model reproduces with fidelity complex temporal variabilities of open-cell oscillations exhibited by high-end cloud simulations.

(stable or unstable) time-periodic solutions and a quasi-periodic solution.

A general, variational approach to derive low-order reduced systems for nonlinear systems subject to an autonomous forcing, is introduced. The approach is based on the concept of optimal parameterizing manifold (PM) that substitutes the more classical notion of slow manifold or invariant manifold when breakdown of slaving occurs. An optimal PM provides the manifold that describes the average motion of the neglected scales as a function of the resolved scales and allows, in principle, for determining the best vector field of the reduced state space that describes e.g. the dynamics' slow motion. The underlying optimal parameterizations are approximated by dynamically-based formulas derived analytically from the original equations. These formulas are contingent upon the determination of only a few (scalar) parameters obtained from minimization of cost functionals, depending on training dataset collected from direct numerical simulation. In practice, a training period of length comparable to a characteristic recurrence or decorrelation time of the dynamics, is sufficient for the efficient derivation of optimized parameterizations. Applications to the closure of low-order models of Atmospheric Primitive Equations and Rayleigh-Bénard convection are then discussed. The approach is finally illustrated --- in the context of the Kuramoto-Sivashinsky turbulence --- as providing efficient closures without slaving for a cutoff scale kc placed within the inertial range and the reduced state space is just spanned by the unstable modes, without inclusion of any stable modes whatsoever. The underlying optimal PMs obtained by our variational approach are far from slaving and allow for remedying the excessive backscatter transfer of energy to the low modes encountered by classical invariant manifold approximations in their standard forms when the latter are used at this cutoff wavelength.

The multiscale variability of the ocean circulation due to its nonlinear dynamics remains a big challenge for theoretical understanding and practical ocean modeling. This paper demonstrates how the data-adaptive harmonic (DAH) decomposition and inverse stochastic modeling techniques introduced in (Chekroun and Kondrashov, (2017), Chaos, 27), allow for reproducing with high fidelity the main statistical properties of multiscale variability in a coarse-grained eddy-resolving ocean flow. This fully-data-driven approach relies on extraction of frequency-ranked time-dependent coefficients describing the evolution of spatio-temporal DAH modes (DAHMs) in the oceanic flow data. In turn, the time series of these coefficients are efficiently modeled by a family of low-order stochastic differential equations (SDEs) stacked per frequency, involving a fixed set of predictor functions and a small number of model coefficients. These SDEs take the form of stochastic oscillators, identified as multilayer Stuart–Landau models (MSLMs), and their use is justified by relying on the theory of Ruelle–Pollicott resonances. The good modeling skills shown by the resulting DAH-MSLM emulators demonstrates the feasibility of using a network of stochastic oscillators for the modeling of geophysical turbulence. In a certain sense, the original quasiperiodic Landau view of turbulence, with the amendment of the inclusion of stochasticity, may be well suited to describe turbulence.

The problem of emergence of fast gravity-wave oscillations in rotating, stratified flow is reconsidered. Fast inertia-gravity oscillations have long been considered an impediment to initialization of weather forecasts, and the concept of a “slow manifold” evolution, with no fast oscillations, has been hypothesized. It is shown on a reduced Primitive Equation model introduced by Lorenz in 1980 that fast oscillations are absent over a finite interval in Rossby number but they can develop brutally once a critical Rossby number is crossed, in contradistinction with fast oscillations emerging according to an exponential smallness scenario such as reported in previous studies, including some others by Lorenz. The consequences of this dynamical transition on the closure problem based on slow variables is also discussed. In that respect, a novel variational perspective on the closure problem exploiting manifolds is introduced. This framework allows for a unification of previous concepts such as the slow manifold or other concepts of “fuzzy” manifold. It allows furthermore for a rigorous identification of an optimal limiting object for the averaging of fast oscillations, namely the optimal parameterizing manifold (PM). It is shown through detailed numerical computations and rigorous error estimates that the manifold underlying the nonlinear Balance Equations provides a very good approximation of this optimal PM even somewhat beyond the emergence of fast and energetic oscillations.

This article proposes a new approach based on finite-horizon parameterizing manifolds (PMs) for the design of low-dimensional suboptimal controllers to optimal control problems of nonlinear partial differential equations (PDEs) of parabolic type. Given a finite horizon [0,T] and a low-mode truncation of the PDE, a PM provides an approximate parameterization of the uncontrolled high modes by the controlled low ones so that the unexplained high-mode energy is reduced, in an L2-sense, when this parameterization is applied. Analytic formulas of such PMs are derived by application of the method of pullback approximation of the high-modes. These formulas allow for an effective derivation of reduced ODE systems, aimed to model the evolution of the low-mode truncation of the controlled state variable, where the high-mode part is approximated by the PM function applied to the low modes. A priori error estimates between the resulting PM-based low-dimensional suboptimal controller u_R* and the optimal controller u* are derived. These estimates demonstrate that the closeness of u_R* to u*? is mainly conditioned on two factors: (i) the parameterization defect of a given PM, associated respectively with u_R* and u*; and (ii) the energy kept in the high modes of the PDE solution either driven by u_R* or u* itself. The practical performances of such PM-based suboptimal controllers are numerically assessed for various optimal control problems associated with a Burgers-type equation. The numerical results show that a PM-based reduced system allows for the design of suboptimal controllers with good performances provided that the associated parameterization defects and energy kept in the high modes are small enough, in agreement with the rigorous results. The practical performances of such PM-based suboptimal controllers are numerically assessed for optimal control problems associated with a Burgers-type equation; the locally as well as globally distributed cases being both considered. The numerical results show that a PM-based reduced system allows for the design of suboptimal controllers with good performances provided that the associated parameterization defects and energy kept in the high modes are small enough, in agreement with the rigorous results.

This first volume is concerned with the analytic derivation of explicit formulas for the leading-order Taylor approximations of (local) stochastic invariant manifolds associated with a broad class of nonlinear stochastic partial differential equations. These approximations take the form of Lyapunov-Perron integrals, which are further characterized in Volume II as pullback limits associated with some partially coupled backward-forward systems. This pullback characterization provides a useful interpretation of the corresponding approximating manifolds and leads to a simple framework that unifies some other approximation approaches in the literature. A self-contained survey is also included on the existence and attraction of one-parameter families of stochastic invariant manifolds, from the point of view of the theory of random dynamical systems.

In this second volume, a general approach is developed to provide approximate parameterizations of the "small" scales by the "large" ones for a broad class of stochastic partial differential equations (SPDEs). This is accomplished via the concept of parameterizing manifolds (PMs), which are stochastic manifolds that improve, for a given realization of the noise, in mean square error the partial knowledge of the full SPDE solution when compared to its projection onto some resolved modes. Backward-forward systems are designed to give access to such PMs in practice. The key idea consists of representing the modes with high wave numbers as a pullback limit depending on the time-history of the modes with low wave numbers. Non-Markovian stochastic reduced systems are then derived based on such a PM approach. The reduced systems take the form of stochastic differential equations involving random coefficients that convey memory effects. The theory is illustrated on a stochastic Burgers-type equation.

In this article, we present a new approach to averaging in non-Hamiltonian systems with periodic forcing. The results here do not depend on the existence of a small parameter. In fact, we show that our averaging method fits into an appropriate nonlinear equivalence problem, and that this problem can be solved formally by using the Lie transform framework to linearize it. According to this approach, we derive formal coordinate transformations associated with both first-order and higher-order averaging, which result in more manageable formulae than the classical ones.

Using these transformations, it is possible to correct the solution of an averaged system by recovering the oscillatory components of the original non-averaged system. In this framework, the inverse transformations are also defined explicitly by formal series; they allow the estimation of appropriate initial data for each higher-order averaged system, respecting the equivalence relation.

Finally, we show how these methods can be used for identifying and computing periodic solutions for a very large class of nonlinear systems with time-periodic forcing. We test the validity of our approach by analyzing both the first-order and the second-order averaged system for a problem in atmospheric chemistry.