Impacts of the Madden-Julian oscillation on intraseasonal precipitation over northeast Brazil

Este estudo tem como objetivo investigar os impactos da oscilação de Madden-Julian (OMJ) na precipitação da região Nordeste do Brasil (NEB). Para tanto foram utilizados dados diários de precipitação baseados em 492 pluviômetros distribuídos na região e cobrindo um período de 30 anos (1981 − 2010). As análises através de composições de anomalias de precipitação, radiação de onda longa e fluxo de umidade, foram obtidas com base no índice da OMJ desenvolvido por Jones-Carvalho. Para distinguir o sinal da OMJ de outros padrões de variabilidade climática, todos os dados diários foram filtrados na escala de 20 − 90 dias; portanto somente dias classificados como eventos da OMJ foram considerados nas composições. Uma análise preliminar baseada apenas nos dados de precipitação foi feita para uma pequena área localizada no interior semiárido do NEB, conhecida como Seridó. Essa microrregião é uma das áreas mais secas do NEB e foi reconhecida pela Convenção das Nações Unidas para o Combate à Desertificação e Mitigação dos Efeitos das Secas como particularmente vulnerável à desertificação. Composições de anomalias de precipitação foram feitas para cada uma das oito fases da OMJ durante Fevereiro-Maio (principal período chuvoso da microrregião). Os resultados mostraram a existência de variações significativas nos padrões de precipitação (de precipitação excessiva à deficiente) associados à propagação da OMJ. A combinação dos sinais de precipitação obtidos durantes as fases úmidas e secas da OMJ mostrou que a diferença corresponde cerca de 50 − 150% de modulação das chuvas na microrregião. Em seguida, uma investigação abrangente sobre o papel da OMJ sobre toda a região Nordeste foi feita considerando-se as quatro estações do ano. Os resultados mostraram que os impactos da OMJ na precipitação intrassazonal do NEB apresentam forte sazonalidade. A maior coerência espacial dos sinais de precipitação ocorreram durante o verão austral, quando cerca de 80% das estações pluviométricas apresentaram anomalias positivas de precipitação durante as fases 1 − 2 da OMJ e anomalias negativas de precipitação nas fases 5 − 6 da oscilação. Embora impactos da OMJ na precipitação intrassazonal tenham sido encontrados na maioria das localidades e em todas as estações do ano, eles apresentaram variações na magnitude dos sinais e dependem da fase da oscilação. As anomalias de precipitação do NEB observadas são explicadas através da interação existente entre as ondas de Kelvin-Rossby acopladas convectivamente e as características climáticas predominantes sobre a região em cada estação do ano. O aumento de precipitação observado sobre a maior parte do NEB durante o verão e primavera austrais encontra-se associado com o fluxo de umidade de oeste (regime de oeste), o qual favorece a atividade convectiva em amplas áreas da América do Sul tropical. Por outro lado, as anomalias de precipitação durante o inverno e outono austrais apresentaram uma variabilidade espacial mais complexa. Durante estas estações, as anomalias de precipitação observadas nas estações localizadas na costa leste do NEB dependem da intensidade do anticiclone do Atlântico Sul, o qual é modulado em grande parte por ondas de Rossby. As características topográficas do NEB parecem desempenhar um papel importante na variabilidade observada na precipitação, principalmente nestas áreas costeiras. A intensificação do anticiclone aumenta a convergência dos ventos alísios na costa contribuindo para a ocorrência de precipitação observada à barlavento do planalto da Borborema. Por outro lado, o aumento da subsidência parece ser responsável pelos déficits de precipitação observados à sotavento. Tais condições mostraram-se típicas durante o predomínio do regime de leste sobre a região tropical da América do Sul e o NEB, durante o qual ocorre uma diminuição no fluxo de umidade proveniente da Amazônia.

Observed Teleconnections in Northern Winter: Subseasonal Evolution and Tropical Linkages

Teleconnections refer to the climate variability links between non-contiguous geographic regions, and tend to be associated with variability in both space and time of the climate’s semi-permanent circulation features. Teleconnections are well-developed in Northern winter, when they influence subseasonal-to-seasonal climate variability, notably, in surface temperature and precipitation. This work is comprised of four independent studies that improve understanding of tropical-extratropical teleconnections and their surface climate responses, subseasonal teleconnection evolution, and the utility of teleconnections in attribution of extreme climate events. After an introduction to teleconnection analysis as well as the major teleconnection patterns and associated climatic footprints manifest during Northern winter, the lagged impact of the Madden-Julian Oscillation (MJO) on subseasonal climate variability is presented. It is found that monitoring of MJO-related velocity potential anomalies is sufficient to predict MJO impacts. These impacts include, for example, the development of significant positive temperature anomalies over the eastern United States one to three weeks following an anomalous convective dipole with enhanced (suppressed) convection centered over the Indian Ocean (western Pacific). Subseasonal teleconnection evolution is assessed with respect to the Pacific-North America (PNA) pattern and the North Atlantic Oscillation (NAO). This evolution is analyzed both in the presence and absence of MJO-related circulation anomalies. It is found that removal of the MJO results only in small shifts in the centers of action of the NAO and PNA, and that in either case there is a small but significant lag in which the NAO leads a PNA pattern of opposite phase. Barotropic vorticity analysis suggests that this relationship may result in part from excitation of Rossby waves by the NAO in the Asian waveguide. An attempt is made to elegantly differentiate between the MJO extratropical response and patterns of variability more internal to the extratropics. Analysis of upper-level streamfunction anomalies is successful in this regard, and it is suggested that this is the preferred method for the real time monitoring of tropical-extratropical teleconnections. The extreme 2013-2014 North American winter is reconstructed using teleconnection analysis, and it is found that the North Pacific Oscillation-West Pacific (NPO/WP) pattern was the leading contributor to climate anomalies over much of North America. Such attribution is cautionary given the propensity to implicate the tropics for all midlatitude climate anomalies based on the El Niño-Southern Oscillation (ENSO) paradigm. A recent hypothesis of such tropical influence is presented and challenged.

The role of the Indian Ocean sector and sea surface salinity for prediction of the coupled Indo-Pacific system

The purpose of this dissertation is to evaluate the potential downstream influence of the Indian Ocean (IO) on El Niño/Southern Oscillation (ENSO) forecasts through the oceanic pathway of the Indonesian Throughflow (ITF), atmospheric teleconnections between the IO and Pacific, and assimilation of IO observations. Also the impact of sea surface salinity (SSS) in the Indo-Pacific region is assessed to try to address known problems with operational coupled model precipitation forecasts. The ITF normally drains warm fresh water from the Pacific reducing the mixed layer depths (MLD). A shallower MLD amplifies large-scale oceanic Kelvin/Rossby waves thus giving ~10% larger response and more realistic ENSO sea surface temperature (SST) variability compared to observed when the ITF is open. In order to isolate the impact of the IO sector atmospheric teleconnections to ENSO, experiments are contrasted that selectively couple/decouple the interannual forcing in the IO. The interannual variability of IO SST forcing is responsible for 3 month lagged widespread downwelling in the Pacific, assisted by off-equatorial curl, leading to warmer NINO3 SST anomaly and improved ENSO validation (significant from 3-9 months). Isolating the impact of observations in the IO sector using regional assimilation identifies large-scale warming in the IO that acts to intensify the easterlies of the Walker circulation and increases pervasive upwelling across the Pacific, cooling the eastern Pacific, and improving ENSO validation (r ~ 0.05, RMS~0.08C). Lastly, the positive impact of more accurate fresh water forcing is demonstrated to address inadequate precipitation forecasts in operational coupled models. Aquarius SSS assimilation improves the mixed layer density and enhances mixing, setting off upwelling that eventually cools the eastern Pacific after 6 months, counteracting the pervasive warming of most coupled models and significantly improving ENSO validation from 5-11 months. In summary, the ITF oceanic pathway, the atmospheric teleconnection, the impact of observations in the IO, and improved Indo-Pacific SSS are all responsible for ENSO forecast improvements, and so each aspect of this study contributes to a better overall understanding of ENSO. Therefore, the upstream influence of the IO should be thought of as integral to the functioning of ENSO phenomenon.

The role of the Rocky Mountains in shaping the atmospheric mean state and the response to tropical forcing in idealised simulations.

Understanding the natural and forced variability of the atmospheric general circulation and its drivers is one of the grand challenges in climate science. It is of paramount importance to understand to what extent the systematic error of climate models affects the processes driving such variability. This is done by performing a set of simulations (ROCK experiments) with an intermediate complexity atmospheric model (SPEEDY), in which the Rocky Mountains orography is increased or decreased to influence the structure of the North Pacific jet stream. For each of these modified-orography experiments, the climatic response to idealized sea surface temperature anomalies of varying intensity in the El Niño Southern Oscillation (ENSO) region is studied. ROCK experiments are characterized by variations in the Pacific jet stream intensity whose extension encompasses the spread of the systematic error found in Coupled Model Intercomparison Project (CMIP6) models. When forced with ENSO-like idealised anomalies, they exhibit a non-negligible sensitivity in the response pattern over the Pacific North American region, indicating that the model mean state can affect the model response to ENSO. It is found that the classical Rossby wave train response to ENSO is more meridionally oriented when the Pacific jet stream is weaker and more zonally oriented with a stronger jet. Rossby wave linear theory suggests that a stronger jet implies a stronger waveguide, which traps Rossby waves at a lower latitude, favouring a zonal propagation of Rossby waves. The shape of the dynamical response to ENSO affects the ENSO impacts on surface temperature and precipitation over Central and North America. A comparison of the SPEEDY results with CMIP6 models suggests a wider applicability of the results to more resources-demanding climate general circulation models (GCMs), opening up to future works focusing on the relationship between Pacific jet misrepresentation and response to external forcing in fully-fledged GCMs.

Numerical simulation of baroclinic waves with a parameterized boundary layer

A dry three-dimensional baroclinic life cycle model is used to investigate the role of turbulent fluxes of heat and momentum within the boundary layer on mid-latitude cyclones. Simulations are performed of life cycles for two basic states, both with and without turbulent fluxes. The different basic states produce cyclones with contrasting frontal and mesoscale-flow structures. The analysis focuses on the generation of potential-vorticity (PV) in the boundary layer and its subsequent transport into the free troposphere. The dynamic mechanism through which friction mitigates a barotropic vortex is that of Ekman pumping. This has often been assumed to be also the dominant mechanism for baroclinic developments. The PV framework highlights an additional, baroclinic mechanism. Positive PV is generated baroclinically due to friction to the north-east of a surface low and is transported out of the boundary layer by a cyclonic conveyor belt flow. The result is an anomaly of increased static stability in the lower troposphere which restricts the growth of the baroclinic wave. The reduced coupling between lower and upper levels can be sufficient to change the character of the upper-level evolution of the mature wave. The basic features of the baroclinic damping mechanism are robust for different frontal structures, with and without turbulent heat fluxes, and for the range of surface roughness found over the oceans.

An interpretation of baroclinic initial value problems: results for simple basic states with nonzero interior PV gradients

In the Eady model, where the meridional potential vorticity (PV) gradient is zero, perturbation energy growth can be partitioned cleanly into three mechanisms: (i) shear instability, (ii) resonance, and (iii) the Orr mechanism. Shear instability involves two-way interaction between Rossby edge waves on the ground and lid, resonance occurs as interior PV anomalies excite the edge waves, and the Orr mechanism involves only interior PV anomalies. These mechanisms have distinct implications for the structural and temporal linear evolution of perturbations. Here, a new framework is developed in which the same mechanisms can be distinguished for growth on basic states with nonzero interior PV gradients. It is further shown that the evolution from quite general initial conditions can be accurately described (peak error in perturbation total energy typically less than 10%) by a reduced system that involves only three Rossby wave components. Two of these are counterpropagating Rossby waves—that is, generalizations of the Rossby edge waves when the interior PV gradient is nonzero—whereas the other component depends on the structure of the initial condition and its PV is advected passively with the shear flow. In the cases considered, the three-component model outperforms approximate solutions based on truncating a modal or singular vector basis.

Inertia-Gravity Waves Emitted from Balanced Flow: Observations, Properties, and Consequences

This paper describes laboratory observations of inertia–gravity waves emitted from balanced fluid flow. In a rotating two-layer annulus experiment, the wavelength of the inertia–gravity waves is very close to the deformation radius. Their amplitude varies linearly with Rossby number in the range 0.05–0.14, at constant Burger number (or rotational Froude number). This linear scaling challenges the notion, suggested by several dynamical theories, that inertia–gravity waves generated by balanced motion will be exponentially small. It is estimated that the balanced flow leaks roughly 1% of its energy each rotation period into the inertia–gravity waves at the peak of their generation. The findings of this study imply an inevitable emission of inertia–gravity waves at Rossby numbers similar to those of the large-scale atmospheric and oceanic flow. Extrapolation of the results suggests that inertia–gravity waves might make a significant contribution to the energy budgets of the atmosphere and ocean. In particular, emission of inertia–gravity waves from mesoscale eddies may be an important source of energy for deep interior mixing in the ocean.

Equatorial waves in the lower stratosphere. II: Annual and interannual variability

The annual and interannual variability of idealized, linear, equatorial waves in the lower stratosphere is investigated using the temperature and velocity fields from the ECMWF 15-year re-analysis dataset. Peak Kelvin wave activity occurs during solstice seasons at 100 hPa, during December-February at 70 hPa and in the easterly to westerly quasi-biennial oscillation (QBO) phase transition at 50 hPa. Peak Rossby-gravity wave activity occurs during equinox seasons at 100 hPa, during June-August/September-November at 70 hPa and in the westerly to easterly QBO phase transition at 50 hPa. Although neglect of wind shear means that the results for inertio-gravity waves are likely to be less accurate, they are still qualitatively reasonable and an annual cycle is observed in these waves at 100 hPa and 70 hPa. Inertio-gravity waves with n = 1 are correlated with the QBO at 50 hPa, but the eastward inertio-gravity n = 0 wave is not, due to its very fast vertical group velocity in all background winds. The relative importance of different wave types in driving the QBO at 50 hPa is also discussed. The strongest acceleration appears to be provided by the Kelvin wave while the acceleration provided by the Rossby-gravity wave is negligible. Of the higher-frequency waves, the westward inertio-gravity n = 1 wave appears able to contribute more to the acceleration of the 50 hPa mean zonal wind than the eastward inertio-gravity n = 1 wave.

Equatorial waves in the lower stratosphere. I: A novel detection method

Wavenumber-frequency spectral analysis and linear wave theory are combined in a novel method to quantitatively estimate equatorial wave activity in the tropical lower stratosphere. The method requires temperature and velocity observations that are regularly spaced in latitude, longitude and time; it is therefore applied to the ECMWF 15-year re-analysis dataset (ERA-15). Signals consistent with idealized Kelvin and Rossby-gravity waves are found at wavenumbers and frequencies in agreement with previous studies. When averaged over 1981-93, the Kelvin wave explains approximately 1 K-2 of temperature variance on the equator at 100 hPa, while the Rossby-gravity wave explains approximately 1 m(2)s(-2) of meridional wind variance. Some inertio-gravity wave and equatorial Rossby wave signals are also found; however the resolution of ERA-15 is not sufficient for the method to provide an accurate climatology of waves with high meridional structure.

Testing the limits of quasi-geostrophic theory: application to observed laboratory flows outside the quasi-geostrophic regime

We compare laboratory observations of equilibrated baroclinic waves in the rotating two-layer annulus, with numerical simulations from a quasi-geostrophic model. The laboratory experiments lie well outside the quasi-geostrophic regime: the Rossby number reaches unity; the depth-to-width aspect ratio is large; and the fluid contains ageostrophic inertia–gravity waves. Despite being formally inapplicable, the quasi-geostrophic model captures the laboratory flows reasonably well. The model displays several systematic biases, which are consequences of its treatment of boundary layers and neglect of interfacial surface tension and which may be explained without invoking the dynamical effects of the moderate Rossby number, large aspect ratio or inertia–gravity waves. We conclude that quasi-geostrophic theory appears to continue to apply well outside its formal bounds.

Equatorial waves in opposite QBO phases

A methodology for identifying equatorial waves is used to analyze the multilevel 40-yr ECMWF Re-Analysis (ERA-40) data for two different years (1992 and 1993) to investigate the behavior of the equatorial waves under opposite phases of the quasi-biennial oscillation (QBO). A comprehensive view of 3D structures and of zonal and vertical propagation of equatorial Kelvin, westward-moving mixed Rossby–gravity (WMRG), and n = 1 Rossby (R1) waves in different QBO phases is presented. Consistent with expectation based on theory, upward-propagating Kelvin waves occur more frequently during the easterly QBO phase than during the westerly QBO phase. However, the westward-moving WMRG and R1 waves show the opposite behavior. The presence of vertically propagating equatorial waves in the stratosphere also depends on the upper tropospheric winds and tropospheric forcing. Typical propagation parameters such as the zonal wavenumber, zonal phase speed, period, vertical wavelength, and vertical group velocity are found. In general, waves in the lower stratosphere have a smaller zonal wavenumber, shorter period, faster phase speed, and shorter vertical wavelength than those in the upper troposphere. All of the waves in the lower stratosphere show an upward group velocity and downward phase speed. When the phase of the QBO is not favorable for waves to propagate, their phase speed in the lower stratosphere is larger and their period is shorter than in the favorable phase, suggesting Doppler shifting by the ambient flow and a filtering of the slow waves. Tropospheric WMRG and R1 waves in the Western Hemisphere also show upward phase speed and downward group velocity, with an indication of their forcing from middle latitudes. Although the waves observed in the lower stratosphere are dominated by “free” waves, there is evidence of some connection with previous tropical convection in the favorable year for the Kelvin waves in the warm water hemisphere and WMRG and R1 waves in the Western Hemisphere, which is suggestive of the importance of convective forcing for the existence of propagating coupled Kelvin waves and midlatitude forcing for the existence of coupled WMRG and R1 waves.

Generation of inertia-gravity waves in the rotating thermal annulus by a localised boundary layer instability

Waves with periods shorter than the inertial period exist in the atmosphere (as inertia-gravity waves) and in the oceans (as Poincaré and internal gravity waves). Such waves owe their origin to various mechanisms, but of particular interest are those arising either from local secondary instabilities or spontaneous emission due to loss of balance. These phenomena have been studied in the laboratory, both in the mechanically-forced and the thermally-forced rotating annulus. Their generation mechanisms, especially in the latter system, have not yet been fully understood, however. Here we examine short period waves in a numerical model of the rotating thermal annulus, and show how the results are consistent with those from earlier laboratory experiments. We then show how these waves are consistent with being inertia-gravity waves generated by a localised instability within the thermal boundary layer, the location of which is determined by regions of strong shear and downwelling at certain points within a large-scale baroclinic wave flow. The resulting instability launches small-scale inertia-gravity waves into the geostrophic interior of the flow. Their behaviour is captured in fully nonlinear numerical simulations in a finite-difference, 3D Boussinesq Navier-Stokes model. Such a mechanism has many similarities with those responsible for launching small- and meso-scale inertia-gravity waves in the atmosphere from fronts and local convection.

The influence of the QBO on the propagation of equatorial waves into the stratosphere

The variation of stratospheric equatorial wave characteristics with the phase of the quasi-biennial oscillation (QBO) is investigated using ECMWF Re-Analysis and NOAA outgoing longwave radiation (OLR) data. The impact of the QBO phases on the upward propagation of equatorial waves is found to be consistent and significant. In the easterly phase, there is larger Kelvin wave amplitude but smaller westward-moving mixed Rossby–gravity (WMRG) and n = 1 Rossby (R1) wave amplitude due to reduced propagation from the upper troposphere into the lower stratosphere, compared with the westerly phase. Differences in the wave amplitude exist in a deeper layer in summer than in winter, consistent with the seasonality of ambient zonal winds. There is a strong evidence of Kelvin wave amplitude peaking just below the descending westerly phase, suggesting that Kelvin waves act to bring the westerly phase downward. However, the corresponding evidence for WMRG and R1 waves is less clear. In the lower stratosphere there is zonal variation in equatorial waves. This reflects the zonal asymmetry of wave amplitudes in the upper troposphere, the source for the lower-stratospheric waves. In easterly winters the upper-tropospheric WMRG and R1 waves over the eastern Pacific region appear to be somewhat stronger compared to climatology, perhaps because of the accumulation of waves that are unable to propagate upward into the lower stratosphere. Vertical propagation features of these waves are generally consistent with theory and suggest a mixture of Doppler shifting by ambient flows and filtering. Some lower-stratosphere equatorial waves have a connection with preceding tropical convection, especially for Kelvin and R1 waves in winter.

A robust mechanism for strengthening of the Brewer–Dobson circulation in response to climate change: critical-layer control of subtropical wave breaking

Climate models consistently predict a strengthened Brewer–Dobson circulation in response to greenhouse gas (GHG)-induced climate change. Although the predicted circulation changes are clearly the result of changes in stratospheric wave drag, the mechanism behind the wave-drag changes remains unclear. Here, simulations from a chemistry–climate model are analyzed to show that the changes in resolved wave drag are largely explainable in terms of a simple and robust dynamical mechanism, namely changes in the location of critical layers within the subtropical lower stratosphere, which are known from observations to control the spatial distribution of Rossby wave breaking. In particular, the strengthening of the upper flanks of the subtropical jets that is robustly expected from GHG-induced tropospheric warming pushes the critical layers (and the associated regions of wave drag) upward, allowing more wave activity to penetrate into the subtropical lower stratosphere. Because the subtropics represent the critical region for wave driving of the Brewer–Dobson circulation, the circulation is thereby strengthened. Transient planetary-scale waves and synoptic-scale waves generated by baroclinic instability are both found to play a crucial role in this process. Changes in stationary planetary wave drag are not so important because they largely occur away from subtropical latitudes.