Search

This paper presents efficient and practical implementations of sequential data assimilation methods for the Simplified Parameterizations, primitive‐Equation DYnamics (SPEEDY) Model, a well‐known numerical model, into the data assimilation community for climate prediction. In the SPEEDY model, the time evolution of dynamics is performed via the second‐order Leapfrog integration scheme; this time integrator relies on two steps: the position and the velocity. The computational implementation of SPEEDY blends the time integrator and the spatial discretization of dynamics to accelerate algebraic computations. Thus, there is no access to the right‐hand side function of the ordinary differential equations governing the time evolution of model dynamics. Consequently, the SPEEDY model is often treated as a black box wherein positions and velocities work as inputs and outputs. Since observations in operational data assimilation only match position states, we can exploit augmented vector states to propagate analysis innovations from positions to velocities. For this purpose, we formulate three variants of ensemble‐based filters and perform numerical experiments to assess their accuracies. We consider two scenarios for the experiments: an ideal case wherein positions and velocities can be observed and a more realistic one wherein measurements are only accessible for position states. Besides, we discuss the effects of the ensemble size on the accuracies of our formulations and, even more, the typical case in which velocities are not updated across assimilation steps. The results reveal that all filter formulations' accuracies remain the same in terms of Root‐Mean‐Square‐Error by neglecting observations from velocities (a realistic scenario) even for cases wherein the number of measurements decreases to 6% of model components. Furthermore, for all discussed filter implementations, the propagation of analysis increments from position to velocities improves up to 100% the performance of filter implementations wherein velocities are not updated, a typical operational scenario. [ABSTRACT FROM AUTHOR]

Quasi‐biennial oscillations (QBOs) in thirteen atmospheric general circulation models forced with both observed and annually repeating sea surface temperatures (SSTs) are evaluated. In most models the QBO period is close to, but shorter than, the observed period of 28 months. Amplitudes are within ±20% of the observed QBO amplitude at 10 hPa, but typically about half of that observed at lower altitudes (50 and 70 hPa). For almost all models, the oscillation's amplitude profile shows an overall upward shift compared to reanalysis and its meridional extent is too narrow. Asymmetry in the duration of eastward and westward phases is reasonably well captured, though not all models replicate the observed slowing of the descending westward shear. Westward phases are generally too weak, and most models have an eastward time mean wind bias throughout the depth of the QBO. The intercycle period variability is realistic and in some models is enhanced in the experiment with observed SSTs compared to the experiment with repeated annual cycle SSTs. Mean periods are also sensitive to this difference between SSTs, but only when parametrized non‐orographic gravity wave (NOGW) sources are coupled to tropospheric parameters and not prescribed with a fixed value. Overall, however, modelled QBOs are very similar whether or not the prescribed SSTs vary interannually. A portrait of the overall ensemble performance is provided by a normalized grading of QBO metrics. To simulate a QBO, all but one model used parametrized NOGWs, which provided the majority of the total wave forcing at altitudes above 70 hPa in most models. Hence the representation of NOGWs either explicitly or through parametrization is still a major uncertainty underlying QBO simulation in these present‐day experiments. [ABSTRACT FROM AUTHOR]

We compare the response of the Quasi‐Biennial Oscillation (QBO) to a warming climate in eleven atmosphere general circulation models that performed time‐slice simulations for present‐day, doubled, and quadrupled CO2 climates. No consistency was found among the models for the QBO period response, with the period decreasing by 8 months in some models and lengthening by up to 13 months in others in the doubled CO2 simulations. In the quadrupled CO2 simulations, a reduction in QBO period of 14 months was found in some models, whereas in several others the tropical oscillation no longer resembled the present‐day QBO, although it could still be identified in the deseasonalized zonal mean zonal wind timeseries. In contrast, all the models projected a decrease in the QBO amplitude in a warmer climate with the largest relative decrease near 60 hPa. In simulations with doubled and quadrupled CO2, the multi‐model mean QBO amplitudes decreased by 36 and 51%, respectively. Across the models the differences in the QBO period response were most strongly related to how the gravity wave momentum flux entering the stratosphere and tropical vertical residual velocity responded to the increases in CO2 amounts. Likewise it was found that the robust decrease in QBO amplitudes was correlated across the models to changes in vertical residual velocity, parametrized gravity wave momentum fluxes, and to some degree the resolved upward wave flux. We argue that uncertainty in the representation of the parameterized gravity waves is the most likely cause of the spread among the eleven models in the QBO's response to climate change. [ABSTRACT FROM AUTHOR]

Convectively generated gravity waves (CGGWs) are important for numerical weather prediction due to their effect on the quasi‐biennial oscillation (QBO) in the stratosphere. Using global ECMWF IFS simulations at TCo7999 (or 1.25 km), TCo2559 (or 3.9 km) and TCo1279 (or 7.8 km) horizontal resolutions, sensitivity of resolved CGGWs to the horizontal resolution and to the explicit versus parametrized representation of deep convection is elucidated during the westerly shear phase of the QBO. Parametrized deep convection is found to inhibit CGGWs, resulting in a twofold reduction in CGGW forcing. When deep convection is explicitly resolved, the total CGGW forcing is almost unchanged across the horizontal resolutions. However, the contribution of long and mesoscale CGGWs (with horizontal wavelengths 100 km ≤λh<1,900 km) to the total CGGW forcing decreases and the contribution of smaller‐scale CGGWs (with λh<100 km) increases as the horizontal resolution increases. At the maximum CGGW forcing altitude, at TCo7999 resolution 43% of the total CGGW forcing is due to long and mesoscale waves, whereas at TCo2559 and TCo1279 resolutions their contribution is 70% and 90%, respectively. While CGGW forcing from long and mesoscale waves is similar at TCo7999 resolution with explicit deep convection and at TCo1279 resolution with parametrized deep convection, CGGW forcing from these waves is artificially enhanced at TCo1279 and TCo2559 resolutions with explicit deep convection. This is due to the explicit deep convection being too strong and having too much variance for 100 km ≤λh<1,900 km. Therefore, parametrizations of deep convection and CGGWs (to account for forcing from waves with λh<100 km) are required even at TCo2559 resolution. Additionally, resolved CGGW forcing at TCo7999 resolution is examined for the easterly shear phase of the QBO; similar to the westerly shear phase, the smaller‐scale waves contribute >55% to the total CGGW forcing at the maximum CGGW forcing altitude. [ABSTRACT FROM AUTHOR]

An ensemble of idealized experiments with the simplified general circulation model PUMA is used to analyze the response to reduced surface friction, that is a strengthening of the eddy‐driven jet, a weakening of the Eulerian mean overturning, and a suppression of baroclinic instability. The suppression of baroclinic instability is caused by an effect called the barotropic governor by which increased horizontal shear restricts the ability of baroclinic disturbances to convert available potential energy into kinetic energy. This governor effect ensures that the residual circulation and Eliassen–Palm flux (EP flux) divergence are largely invariant to the surface friction parameter despite the connection between surface friction, the Eulerian mean overturning, and the eddy‐momentum flux. The suppression of instability leads to an increase in persistence measured by the period of peak variance on synoptic time‐scales and a strengthened signal‐to‐noise ratio on seasonal time‐scales. These findings suggest that the signal‐to‐noise paradox seen in the context of seasonal prediction can be caused by excess mechanical damping in atmospheric prediction systems inhibiting the barotropic governor effect. [ABSTRACT FROM AUTHOR]