A Python Package to Calculate the OLR-Based Index of the Madden- Julian-Oscillation (OMI) in Climate Science and Weather Forecasting

The Madden-Julian Oscillation (MJO) is a prominent feature of the intraseasonal variability of the atmosphere. The MJO strongly modulates tropical precipitation and has implications around the globe for weather, climate and basic atmospheric research. The time-dependent state of the MJO is described by MJO indices, which are calculated through sometimes complicated statistical approaches from meteorological variables. One of these indices is the OLR-based MJO Index (OMI; OLR stands for outgoing longwave radiation). The Python package mjoindices, which is described in this paper, provides the first open source implementation of the OMI algorithm, to our knowledge. The package meets state-of-the-art criteria for sustainable research software, like automated tests and a persistent archiving to aid the reproducibility of scientific results. The agreement of the OMI values calculated with this package and the original OMI values is also summarized here. There are several reuse scenarios; the most probable one is MJO-related research based on atmospheric models, since the index values have to be recalculated for each model run

A Stochastic Parameterization of Organized Tropical Convection Using Cellular Automata for Global Forecasts in NOAA's Unified Forecast System

In the atmosphere, convection can organize from smaller scale updrafts into more coherent structures on various scales. In bulk-plume cumulus convection parameterizations, this type of organization has to be represented in terms of how the resolved flow would &ldquo;feel&rdquo; convection if more coherent structures were present on the subgrid. This type of subgrid organization acts as building blocks for larger scale tropical convective organization known to modulate local and remote weather. In this work a parameterization for subgrid (and cross-grid) organization in a bulk-plume convection scheme is proposed using the stochastic, self-organizing, properties of cellular automata (CA). We investigate the effects of using a CA which can interact with three different components of the bulk-plume scheme that modulate convective activity: entrainment, triggering, and closure. The impacts of the revised schemes are studied in terms of the model's ability to organize convectively coupled equatorial waves (CCEWs). The differing impacts of adopting the stochastic CA scheme, as compared to the widely used Stochastically Perturbed Physics Tendency (SPPT) scheme, are also assessed. Results show that with the CA scheme, precipitation is more spatially and temporally organized, and there is a systematic shift in equatorial wave phase speed not seen with SPPT. Previous studies have noted a linear relationship between Gross Moist Stability (GMS) and Kelvin wave phase speed. Analysis of GMS in this study shows an increase in Kelvin wave phase speed and an increase in GMS with the CA scheme, which is tied to a shift from large-scale precipitation to convective precipitation. </div

The Intricacies of Identifying Equatorial Waves

Equatorial waves (EWs) are synoptic- to planetary-scale propagating disturbances at low latitudes with periods from a few days to several weeks. Here, this term includes Kelvin waves, equatorial Rossby waves, mixed Rossby–gravity waves, and inertio-gravity waves, which are well described by linear wave theory, but it also other tropical disturbances such as easterly waves and the intraseasonal Madden–Julian Oscillation with more complex dynamics. EWs can couple with deep convection, leading to a substantial modulation of clouds and rainfall. EWs are amongst the dynamic features of the troposphere with the longest intrinsic predictability, and models are beginning to forecast them with an exploitable level of skill. Most of the methods developed to identify and objectively isolate EWs in observations and model fields rely on (or at least refer to) the adiabatic, frictionless linearized primitive equations on the sphere or the shallow-water system on the equatorial -plane. Common ingredients to these methods are zonal wave-number–frequency filtering (Fourier or wavelet) and/or projections onto predefined empirical or theoretical dynamical patterns. This paper gives an overview of six different methods to isolate EWs and their structures, discusses the underlying assumptions, evaluates the applicability to different problems, and provides a systematic comparison based on a case study (February 20–May 20, 2009) and a climatological analysis (2001–2018). In addition, the influence of different input fields (e.g., winds, geopotential, outgoing long-wave radiation, rainfall) is investigated. Based on the results, we generally recommend employing a combination of wave-number–frequency filtering and spatial-projection methods (and of different input fields) to check for robustness of the identified signal. In cases of disagreement, one needs to carefully investigate which assumptions made for the individual methods are most probably not fulfilled. This will help in choosing an approach optimally suited to a given problem at hand and avoid misinterpretation of the results

The intricacies of identifying equatorial waves

Equatorial waves (EWs) are synoptic- to planetary-scale propagating disturbances at low latitudes with periods from a few days to several weeks. Here, this term includes Kelvin waves, equatorial Rossby waves, mixed Rossby–gravity waves, and inertio-gravity waves, which are well described by linear wave theory, but it also other tropical disturbances such as easterly waves and the intraseasonal Madden–Julian Oscillation with more complex dynamics. EWs can couple with deep convection, leading to a substantial modulation of clouds and rainfall. EWs are amongst the dynamic features of the troposphere with the longest intrinsic predictability, and models are beginning to forecast them with an exploitable level of skill. Most of the methods developed to identify and objectively isolate EWs in observations and model fields rely on (or at least refer to) the adiabatic, frictionless linearized primitive equations on the sphere or the shallow-water system on the equatorial -plane. Common ingredients to these methods are zonal wave-number–frequency filtering (Fourier or wavelet) and/or projections onto predefined empirical or theoretical dynamical patterns. This paper gives an overview of six different methods to isolate EWs and their structures, discusses the underlying assumptions, evaluates the applicability to different problems, and provides a systematic comparison based on a case study (February 20–May 20, 2009) and a climatological analysis (2001–2018). In addition, the influence of different input fields (e.g., winds, geopotential, outgoing long-wave radiation, rainfall) is investigated. Based on the results, we generally recommend employing a combination of wave-number–frequency filtering and spatial-projection methods (and of different input fields) to check for robustness of the identified signal. In cases of disagreement, one needs to carefully investigate which assumptions made for the individual methods are most probably not fulfilled. This will help in choosing an approach optimally suited to a given problem at hand and avoid misinterpretation of the results