Section: New Results
Modeling for Oceanic and Atmospheric flows
Coupling Methods for Oceanic and Atmospheric Models
Participants : Eric Blayo, Florian Lemarié, Sophie Thery.
Coupling methods routinely used in regional and global climate models do not provide the exact solution to the oceanatmosphere problem, but an approximation of one [72]. For the last few years we have been actively working on the analysis of oceanatmosphere coupling both in terms of its continuous and numerical formulation.Our activities over the last few years can be divided into four general topics

Stability and consistency analysis of existing coupling methods: in [72] we showed that the usual methods used in the context of oceanatmosphere coupling are prone to splitting errors because they correspond to only one iteration of an iterative process without reaching convergence. Moreover, those methods have an additional condition for the coupling to be stable even if unconditionally stable time stepping algorithms are used. This last remark was further studied recently in [47] and it turned out to be a major source of instability in atmospheresnow coupling.

Study of physicsdynamics coupling: during the PhDthesis of Charles Pelletier (funded by Inria and defended on Feb. 15, 2018, [2]) the scope was on including the formulation of physical parameterizations in the theoretical analysis of the coupling, in particular the parameterization schemes to compute airsea fluxes [79]. To do so, a metamodel representative of the behavior of the full parameterization but with a continuous form easier to manipulate has been derived thanks to a sensitivity analysis. This metamodel is more adequate to conduct the mathematical analysis of the coupling while being physically satisfactory [80]. In parallel we have contributed to a general review gathering the main international specialists on the topic [64]. More recently we have started to work specifically on the discretization methods for the parameterization of planetary boundary layers in climate models [27] which takes the form of an nonstationary nonlinear parabolic equation. The objective is to derive a discretization for which we could prove nonlinear stability criteria and show robustness to large variations in parabolic Courant number while being consistent with our knowledge of the underlying physical principles (e.g. the MoninObukhov theory in the surface layer).

Design of a coupled single column model: in order to focus on specific problems of oceanatmosphere coupling, a work on simplified equation sets has been started. The aim is to implement a onedimensional (in the vertical direction) coupled model with physical parameterizations representative of those used in realistic models. Thanks to this simplified coupled model the objective is to develop a benchmark suite for coupled models evaluation. Last year the single column oceanic and atmospheric components have been developed and coupled during the PhDthesis of Rémi Pellerej (defended on Mar. 26, 2018) and in the framework of the SIMBAD project [17]. A publication describing this model is currently in preparation for the Geoscientific Model Development journal.

Analysis of airseawave interactions in realistic highresolution realistic simulations: part of our activity has been in collaboration with atmosphericists and physical oceanographers to study the impact on some modeling assumptions (e.g. [73]) in largescale realistic oceanatmosphere coupled simulations [14]. Moreover, within the ALBATROS project, we have contributed to the development of a 2way coupling between an ocean global circulation model (NEMO) with a surface wave model (WW3). Such coupling is not straightforward to implement since it requires modifications of the governing equations, boundary conditions and subgrid scale closures in the oceanic model. A paper is currently under review in Geoscientific Model Development journal on that topic.

Efficient coupling methods: we have been developing coupling approaches for several years, based on socalled Schwarz algorithms. In particular, we addressed the development of efficient interface conditions for multiphysics problems representative of airsea coupling [28] (paper in preparation). This work is done in the framework of S. Théry PhD (started in fall 2017).
These topics are addressed through strong collaborations between the applied mathematicians and the climate community (MeteoFrance, Ifremer, LMD, and LOCEAN). Indeed, Our work on oceanatmosphere coupling has steadily matured over the last few years and has reached a point where it triggered interest from the climate community. Through the funding of the COCOA ANR project (started in January 2017, PI: E. Blayo), Airsea team members play a major role in the structuration of a multidisciplinary scientific community working on oceanatmosphere coupling spanning a broad range from mathematical theory to practical implementations in climate models. An expected outcome of this project should be the design of a benchmark suite of idealized coupled test cases representative of known issues in coupled models. Such idealized test cases should motivate further collaborations at an international level.
Numerical Schemes for Ocean Modelling
Participants : Eric Blayo, Matthieu Brachet, Laurent Debreu, Emilie Duval, Nicholas Kevlahan, Florian Lemarié, Christopher Eldred, Farshid Nazari.
The increase of model resolution naturally leads to the representation of a wider energy spectrum. As a result, in recent years, the understanding of oceanic submesoscale dynamics has significantly improved. However, dissipation in submesoscale models remains dominated by numerical constraints rather than physical ones. Effective resolution is limited by the numerical dissipation range, which is a function of the model numerical filters (assuming that dispersive numerical modes are efficiently removed). A review paper on coastal ocean models has been written with German colleagues and has been published in Ocean Modelling ([11]).
F. Lemarié and L. Debreu (with H. Burchard, K. Klingbeil and J. SainteMarie) have organized the international COMMODORE workshop on numerical methods for oceanic models (Paris, Sept. 1719, 2018). https://commodore2018.sciencesconf.org/, see [12] for a summary of the scientific discussions
With the increase of resolution, the hydrostatic assumption becomes less valid and the AIRSEA group also works on the development of nonhydrostatic ocean models. The treatment of nonhydrostatic incompressible flows leads to a 3D elliptic system for pressure that can be ill conditioned in particular with non geopotential vertical coordinates. That is why we flavor the use of the nonhydrostatic compressible equations that removes the need for a 3D resolution at the price of reincluding acoustic waves [24].
In addition, Emilie Duval started her PhD in September 2018 on the coupling between the hydrostatic incompressible and nonhydrostatic compressible equations.
The team is involved in the HEAT (Highly Efficient ATmospheric Modelling) ANR project. This project aims at developing a new atmospheric dynamical core (DYNAMICO) discretized on an icosahedral grid. This project is in collaboration with Ecole Polytechnique, MeteoFrance, LMD, LSCE and CERFACS. This year we worked on dispersion analysis of compatible Galerkin schemes for shallow water model both in 1D ([5]) and 2D ([39]). In addition, we worked on the discrete formulation of the thermal rotating shallow water equations. This formulation, based on quasiHamiltonian discretizations methods, allows for the first time total mass, buoyancy and energy conservation to machine precision ([4]).
Data assimilation for coupled models
In the context of operational meteorology and oceanography, forecast skills heavily rely on proper combination of model prediction and available observations via data assimilation techniques. Historically, numerical weather prediction is made separately for the ocean and the atmosphere in an uncoupled way. However, in recent years, fully coupled oceanatmosphere models are increasingly used in operational centers to improve the reliability of seasonal forecasts and tropical cyclones predictions. For coupled problems, the use of separated data assimilation schemes in each medium is not satisfactory since the result of such assimilation process is generally inconsistent across the interface, thus leading to unacceptable artefacts. Hence, there is a strong need for adapting existing data assimilation techniques to the coupled framework. As part of our ERACLIM2 contribution, R. Pellerej started a PhD on that topic late 2014 and defended it early 2018 [1]. Three general data assimilation algorithms, based on variational data assimilation techniques, have been developed and applied to a single column coupled model. The dynamical equations of the considered problem are coupled using an iterative Schwarz domain decomposition method. The aim is to properly take into account the coupling in the assimilation process in order to obtain a coupled solution close to the observations while satisfying the physical conditions across the airsea interface. Results shows significant improvement compared to the usual approach on this simple system. The aforementioned system has been coded within the OOPS framework (Object Oriented Prediction System) in order to ease the transfer to more complex/realistic models.
Finally, CASIS, a new collaborative project with Mercator Océan has started late 2017 in order to extend developments to iterative Kalman smoother data assimilation scheme, in the framework of a coupled oceanatmospheric boundary layer context.
Optimal control of grids and schemes for ocean model.
Participants : Laurent Debreu, Eugene Kazantsev.
In [33], variational data assimilation technique is applied to a simple bidimensional wave equation that simulates propagation of internal gravity waves in the ocean in order to control grids and numerical schemes. Grid steps of the vertical grid, BruntVaisala frequency and approximation of the horizontal derivative were used as control parameters either separately or in the joint control. Obtained results show that optimized parameters may partially compensate errors committed by numerical scheme due to insufficient grid resolution.
Optimal vertical grid steps and coefficients in horizontal derivative approximation found in the variational control procedure allow us to get the model solution that is rather close to the solution of the reference model. The error in the wave velocity on the coarse grid is mostly compensated in experiments with joint control of parameters while the error in the wave amplitude occurs to be more difficult to correct.
However, optimal grid steps and discretization schemes may be in a disagreement with requirements of other model physics and additional analysis of obtained optimized parameters from the point of view of they agreement with the model is necessary.
Nonhydrostatic Modeling
Participants : Eric Blayo, Laurent Debreu, Emilie Duval.
In the context of the French initiative CROCO (Coastal and Regional Ocean COmmunity model, https://www.crocoocean.org) for the development of a new oceanic modeling system, Emilie Duval started a PhD (Oct. 2018) focused on the design of methods to couple local nonhydrostatic models to larger scale hydrostatic ones. Such a coupling is quite delicate from a mathematical point of view, due to the different nature of hydrostatic and nonhydrostatic equations (where the vertical velocity is either a diagnostic or a prognostic variable).