skip to main content

This content will become publicly available on November 1, 2022

Title: Learning Forecasts of Rare Stratospheric Transitions from Short Simulations
Abstract Rare events arising in nonlinear atmospheric dynamics remain hard to predict and attribute. We address the problem of forecasting rare events in a prototypical example, sudden stratospheric warmings (SSWs). Approximately once every other winter, the boreal stratospheric polar vortex rapidly breaks down, shifting midlatitude surface weather patterns for months. We focus on two key quantities of interest: the probability of an SSW occurring, and the expected lead time if it does occur, as functions of initial condition. These optimal forecasts concretely measure the event’s progress. Direct numerical simulation can estimate them in principle but is prohibitively expensive in practice: each rare event requires a long integration to observe, and the cost of each integration grows with model complexity. We describe an alternative approach using integrations that are short compared to the time scale of the warming event. We compute the probability and lead time efficiently by solving equations involving the transition operator, which encodes all information about the dynamics. We relate these optimal forecasts to a small number of interpretable physical variables, suggesting optimal measurements for forecasting. We illustrate the methodology on a prototype SSW model developed by Holton and Mass and modified by stochastic forcing. While highly idealized, more » this model captures the essential nonlinear dynamics of SSWs and exhibits the key forecasting challenge: the dramatic separation in time scales between a single event and the return time between successive events. Our methodology is designed to fully exploit high-dimensional data from models and observations, and has the potential to identify detailed predictors of many complex rare events in meteorology. « less
; ; ; ;
Award ID(s):
Publication Date:
Journal Name:
Monthly Weather Review
Page Range or eLocation-ID:
3647 to 3669
Sponsoring Org:
National Science Foundation
More Like this
  1. Although sudden stratospheric warmings (SSWs) can improve subseasonal-to-seasonal forecasts, it is unclear whether the two types of SSW - displacements and splits - have different near- surface effects. To examine the longer-term (i.e., multi-week lead) tropospheric response to displacements and splits, we utilize an intermediate-complexity model and impose wave-1 and wave-2 stratospheric heating perturbations spun-off from a control run. At longer lags, the tropospheric response is found to be insensitive to both the wavenumber and location of the imposed heating, in agreement with freely evolving displacements and splits identified in the control run. At shorter lags, however, large differences aremore »found between displacements and splits in both the control run and the different wavenumber- forced events. In particular, in the control run, the free-running splits have an immediate barotropic response throughout the stratosphere and troposphere whereas displacements take 1–2 weeks before a near-surface response becomes evident. Interestingly, this barotropic response found during CTRL splits is not captured by the barotropically forced wave-2 events, indicating that the zonal-mean tropospheric circulation is somehow coupled with the generation of the wave-2 splits. It is also found that in the control run, displacements yield stronger Polar-Cap temperature anomalies than splits, yet both still yield similar magnitude tropospheric responses. Hence, the strength of the stratospheric warming is not the only governing factor in the surface response. Overall, SSW classification based on vortex morphology may be useful for subseasonal but not seasonal tropospheric prediction.« less
  2. The tropospheric response to midwinter sudden stratospheric warmings (SSWs) is examined using an idealized model. SSW events are triggered by imposing high-latitude stratospheric heating perturbations of varying magnitude for only a few days, spun off from a free-running control integration (CTRL). The evolution of the thermally triggered SSWs is then compared with naturally occurring SSWs identified in CTRL. By applying a heating perturbation, with no modification to the momentum budget, it is possible to isolate the tropospheric response directly attributable to a change in the stratospheric polar vortex, independent of any planetary wave momentum torques involved in the initiation ofmore »an SSW. Zonal-wind anomalies associated with the thermally triggered SSWs first propagate downward to the high-latitude troposphere after ~2 weeks, before migrating equatorward and stalling at midlatitudes, where they straddle the near-surface jet. After ~3 weeks, the circulation and eddy fluxes associated with thermally triggered SSWs evolve very similarly to SSWs in CTRL, despite the lack of initial planetary wave driving. This suggests that at longer lags, the tropospheric response to SSWs is generic and it is found to be linearly governed by the strength of the lower-stratospheric warming, whereas at shorter lags, the initial formation of the SSW potentially plays a large role in the downward coupling. In agreement with previous studies, synoptic waves are found to play a key role in the persistent tropospheric jet shift at long lags. Synoptic waves appear to respond to the enhanced midlatitude baroclinicity associated with the tropospheric jet shift, and preferentially propagate poleward in an apparent positive feedback with changes in the high-latitude refractive index.« less
  3. Abstract The Whole Atmosphere Community Climate Model, version 4 (WACCM4), is used to investigate the influence of stratospheric conditions on the development of sudden stratospheric warmings (SSWs). To this end, targeted experiments are performed on selected modeled SSW events. Specifically, the model is reinitialized three weeks before a given SSW, relaxing the surface fluxes, winds, and temperature below 10 km to the corresponding fields from the free-running simulation. Hence, the tropospheric wave evolution is unaltered across the targeted experiments, but the stratosphere itself can evolve freely. The stratospheric zonal-mean state is then altered 21 days prior to the selected SSWsmore »and rerun with an ensemble of different initial conditions. It is found that a given tropospheric evolution concomitant with the development of an SSW does not uniquely determine the occurrence of an event and that the stratospheric conditions are relevant to the subsequent evolution of the stratospheric flow toward an SSW, even for a fixed tropospheric evolution. It is also shown that interpreting the meridional heat flux at 100 hPa as a proxy of the tropospheric injection of wave activity into the stratosphere should be regarded with caution and that stratospheric dynamics critically influence the heat flux at that altitude.« less
  4. Abstract

    The frequency of stratospheric sudden warming events (SSWs) is an important characteristic of the coupled stratosphere–troposphere system. However, many modern climate models are unable to reproduce the observed SSW frequency. A previous study suggested that one of the reasons could be the momentum damping at the surface. The goal of the present study is to understand what determines the climatological SSW frequency and how the surface damping comes into play. To this end, we conduct a parameter sweep with an idealized model, using a wide range of values for the surface damping. It is found that the SSW frequencymore »is a strong and nonlinear function of the surface damping. Various tropospheric and stratospheric factors are identified to link the surface damping to the SSW frequency. The factors include the magnitude of the surface winds, the meridional and vertical wind shear, the synoptic eddy activity in the troposphere, the transient wave activity flux at the lower stratosphere, and the strength of the stratospheric polar vortex. Mathematical–statistical modeling, informed by the parameter sweep, is used to quantify how the different factors relate to each other. This successfully reproduces the complex variations of the SSW frequency with the surface damping seen in the parameter sweep. The results may help in explaining some of the difficulties that climate models have in simulating the observed SSW frequency.

    « less
  5. Abstract

    Forecasting the El Niño-Southern Oscillation (ENSO) has been a subject of vigorous research due to the important role of the phenomenon in climate dynamics and its worldwide socioeconomic impacts. Over the past decades, numerous models for ENSO prediction have been developed, among which statistical models approximating ENSO evolution by linear dynamics have received significant attention owing to their simplicity and comparable forecast skill to first-principles models at short lead times. Yet, due to highly nonlinear and chaotic dynamics (particularly during ENSO initiation), such models have limited skill for longer-term forecasts beyond half a year. To resolve this limitation, heremore »we employ a new nonparametric statistical approach based on analog forecasting, called kernel analog forecasting (KAF), which avoids assumptions on the underlying dynamics through the use of nonlinear kernel methods for machine learning and dimension reduction of high-dimensional datasets. Through a rigorous connection with Koopman operator theory for dynamical systems, KAF yields statistically optimal predictions of future ENSO states as conditional expectations, given noisy and potentially incomplete data at forecast initialization. Here, using industrial-era Indo-Pacific sea surface temperature (SST) as training data, the method is shown to successfully predict the Niño 3.4 index in a 1998–2017 verification period out to a 10-month lead, which corresponds to an increase of 3–8 months (depending on the decade) over a benchmark linear inverse model (LIM), while significantly improving upon the ENSO predictability “spring barrier”. In particular, KAF successfully predicts the historic 2015/16 El Niño at initialization times as early as June 2015, which is comparable to the skill of current dynamical models. An analysis of a 1300-yr control integration of a comprehensive climate model (CCSM4) further demonstrates that the enhanced predictability afforded by KAF holds over potentially much longer leads, extending to 24 months versus 18 months in the benchmark LIM. Probabilistic forecasts for the occurrence of El Niño/La Niña events are also performed and assessed via information-theoretic metrics, showing an improvement of skill over LIM approaches, thus opening an avenue for environmental risk assessment relevant in a variety of contexts.

    « less