Atmospheric processes

Abstract:

The mesoscale atmospheric energy spectrum has puzzled scientists for decades, sitting between classical turbulence and wave theories. Using year-long ECMWF operational analyses of high resolution and a spherical coarse-graining framework (Flowsieve), we present the first consistent global maps of local mesoscale kinetic energy fluxes. At 200~hPa, we identify a striking band of upscale transfer aligned with the ITCZ, while storm tracks and orography leave distinct dynamical imprints at both 200 and 600~hPa. By decomposing divergent and rotational components, we show that divergent energy dominates in the tropics and stratosphere, while rotational energy dominates in the extratropical troposphere. Conditioning spectra on this balance reveals contrasting regimes: a Nastrom–Gage-like spectrum under divergent dominance, and a spectrum reminiscent of the classical dual cascade of textbook two-dimensional turbulence under rotational dominance at 600~hPa. These results demonstrate that mesoscale energy transfer is shaped by a patchwork of mechanisms, reconciling long-standing debates and providing new inspiration for parametrisations and predictability in weather and climate models.

Abstract:

Machine learning hold the promise of unlocking more accurate and realistic parameterizations of atmospheric processes, but brings its own set of challenges and drawbacks. Among top issues are generalization, stability and interpretability. Here we present a parameter-efficient neural network parameterization which aims to address these issues by incorporating physical knowledge to a high degree. By predicting fluxes and microphysical process rates instead of total tendencies, the conservation of water can be hardcoded, which is shown to improve online performance. Furthermore, a physically motivated architecture based on vertically recurrent neural networks enables high computational efficiency and a low number of parameters. The models are trained and evaluated using a superparameterization setup with real orography. The impact of incorporating stochasticity is also discussed.

Abstract:

Representing and quantifying uncertainty in physical parameterisations is a central challenge in weather and climate modelling, and approaches are often developed separately for different timescales. Here, we consider the separation of uncertainty by source using machine learning frameworks for subgrid-scale parameterisations. In this context, aleatoric uncertainty arises from internal variability in the training data, and epistemic uncertainty, arises from poorly constrained parameters during training. Using the Lorenz 1996 system as a testbed for simplified chaotic dynamics, we deal with uncertainties through a unified framework using Bayesian Neural Networks, to explore how the different sources of uncertainty evolve over different prediction timescales.

Abstract:

In the Model Uncertainty-Model Intercomparison Project (MUMIP) we compare parameterisation packages from different modelling centres using their single-column modelling (SCM) frameworks. We will showcase the dataset from an Indian Ocean experiment at a 0.2 degrees grid covering one month, with about 10 million simulations of each model. These parametrised models are compared against a convection-permitting benchmark from DYAMOND under common dynamical constraints. We will show differences and similarities in precipitation patterns and physics tendencies among four models and show how these differences can be generalised. Following earlier works, we find that at coarse grids that do not resolve convection, parameterisation packages tend to produce overconfident tendencies compared to the convection-permitting benchmark. Furthermore, we test several hypotheses on the MUMIP dataset to explain the differences. We use the data to explore the foundations of stochastic physical parametrisations. Would stochastic physics effectively overcome the overconfidence for good reasons? May the stochastic perturbations actually have a physically meaningful quantitative interpretation? Can stochastic physics be used to partially overcome truncation and grid spacing limitations?

Abstract:

The Atlantic Multidecadal Variability (AMV) and the North Atlantic Oscillation (NAO) are the dominant modes of oceanic and atmospheric variability in the North Atlantic, respectively, and are key sources of predictability from seasonal to decadal timescales. However, the physical processes and feedback mechanisms linking the AMV and NAO, and the role of diabatic processes in these feedbacks, remain debated. We present a data-driven dynamical modelling framework which captures coupled decadal variability in AMV, NAO, and North Atlantic precipitation. Applying equation discovery methods to observational data, we identify deterministic low-order dynamical models consisting of three coupled ordinary differential equations. These models reproduce observed North Atlantic decadal variability and show robust out-of-sample predictive skill on multi-annual to decadal lead times. The resulting model dynamics include a distinct quasi-periodic 20-year oscillation consistent with a damped oceanic mode of variability. Notably, precipitation-related terms feature prominently in the low-order models, suggesting an important role for latent heat release and freshwater fluxes in mediating ocean–atmosphere interactions. We propose new feedback mechanisms between North Atlantic sea surface temperature and the NAO, with precipitation acting as a dynamical bridge. By linearising the low-order models and computing finite-time Lyapunov exponents, we find that North Atlantic precipitation is more predictable in a positive AMV phase. We then analyse several decadal prediction ensemble experiments based on initialised hindcasts and find comparable state-dependent predictability of precipitation. Overall, these results illustrate how data-driven equation discovery can provide mechanistic hypotheses and new insight beyond conventional analyses of observations and climate model simulations.