Point vortices on the sphere: a case with opposite vorticities

We study systems formed of 2N point vortices on a sphere with N vortices of strength +1 and N vortices of strength -1. In this case, the Hamiltonian is conserved by the symmetry which exchanges the positive vortices with the negative vortices. We prove the existence of some fixed and relative equilibria, and then study their stability with the ``Energy Momentum Method''. Most of the results obtained are nonlinear stability results. To end, some bifurcations are described.Comment: 35 pages, 9 figure

Deformation of geometry and bifurcation of vortex rings

We construct a smooth family of Hamiltonian systems, together with a family of group symmetries and momentum maps, for the dynamics of point vortices on surfaces parametrized by the curvature of the surface. Equivariant bifurcations in this family are characterized, whence the stability of the Thomson heptagon is deduced without recourse to the Birkhoff normal form, which has hitherto been a necessary tool.Comment: 26 page

Generation and Structure of Solitary Rossby Vortices in Rotating Fluids

The formation of zonal flows and vortices in the generalized Charney-Hasegawa-Mima equation is studied. We focus on the regime when the size of structures is comparable to or larger than the deformation (Rossby) radius. Numerical simulations show the formation of anticyclonic vortices in unstable shear flows and ring-like vortices with quiescent cores and vorticity concentrated in a ring. Physical mechanisms that lead to these phenomena and their relevance to turbulence in planetary atmospheres are discussed.Comment: 3 pages in REVTeX, 5 postscript figures separately, submitted to Phys. Rev.

Isolating the Roles of Different Forcing Agents in Global Stratospheric Temperature Changes Using Model Integrations with Incrementally Added Single Forcings

Satellite instruments show a cooling of global stratospheric temperatures over the whole data record (1979-2014). This cooling is not linear and includes two descending steps in the early 1980s and mid-1990s. The 1979-1995 period is characterized by increasing concentrations of ozone depleting substances (ODS) and by the two major volcanic eruptions of El Chichon (1982) and Mount Pinatubo (1991). The 1995-present period is characterized by decreasing ODS concentrations and by the absence of major volcanic eruptions. Greenhouse gas (GHG) concentrations increase over the whole time period. In order to isolate the roles of different forcing agents in the global stratospheric temperature changes, we performed a set of AMIP-style simulations using the NASA Goddard Earth Observing System Chemistry-Climate Model (GEOSCCM). We find that in our model simulations the cooling of the stratosphere from 1979 to present is mostly driven by changes in GHG concentrations in the middle and upper stratosphere and by GHG and ODS changes in the lower stratosphere. While the cooling trend caused by increasing GHGs is roughly constant over the satellite era, changing ODS concentrations cause a significant stratospheric cooling only up to the mid-1990s, when they start to decrease because of the implementation of the Montreal Protocol. Sporadic volcanic events and the solar cycle have a distinct signature in the time series of stratospheric temperature anomalies but do not play a statistically significant role in the long-term trends from 1979 to 2014. Several factors combine to produce the step-like behavior in the stratospheric temperatures: in the lower stratosphere, the flattening starting in the mid-1990s is due to the decrease in ozone-depleting substances; Mount Pinatubo and the solar cycle cause the abrupt steps through the aerosol-associated warming and the volcanically induced ozone depletion. In the middle and upper stratosphere, changes in solar irradiance are largely responsible for the step-like behavior of global temperature anomalies, together with volcanically induced ozone depletion and water vapor increases in the post-Pinatubo years

Controls on boundary layer ventilation: Boundary layer processes and large-scale dynamics

Midlatitude cyclones are important contributors to boundary layer ventilation. However, it is uncertain how efficient such systems are at transporting pollutants out of the boundary layer, and variations between cyclones are unexplained. In this study 15 idealized baroclinic life cycles, with a passive tracer included, are simulated to identify the relative importance of two transport processes: horizontal divergence and convergence within the boundary layer and large-scale advection by the warm conveyor belt. Results show that the amount of ventilation is insensitive to surface drag over a realistic range of values. This indicates that although boundary layer processes are necessary for ventilation they do not control the magnitude of ventilation. A diagnostic for the mass flux out of the boundary layer has been developed to identify the synoptic-scale variables controlling the strength of ascent in the warm conveyor belt. A very high level of correlation (R-2 values exceeding 0.98) is found between the diagnostic and the actual mass flux computed from the simulations. This demonstrates that the large-scale dynamics control the amount of ventilation, and the efficiency of midlatitude cyclones to ventilate the boundary layer can be estimated using the new mass flux diagnostic. We conclude that meteorological analyses, such as ERA-40, are sufficient to quantify boundary layer ventilation by the large-scale dynamics

Can Natural Variability Explain Observed Antarctic Sea Ice Trends? New Modeling Evidence from CMIP5

The recent observed positive trends in total Antarctic sea ice extent are at odds with the expectation of melting sea ice in a warming world. More problematic yet, climate models indicate that sea ice should decrease around Antarctica in response to both increasing greenhouse gases and stratospheric ozone depletion. The resolution of this puzzle, we suggest, may lie in the large natural variability of the coupled atmosphere‒ocean‒sea‒ice system. Contrasting forced and control integrations from four state‒of‒the‒art Coupled Model Intercomparison Project Phase 5 (CMIP5) models, we show that the observed Antarctic sea ice trend falls well within the distribution of trends arising naturally in the system, and that the forced response in the models is small compared to the natural variability. From this, we conclude that it may prove difficult to attribute the observed trends in total Antarctic sea ice to anthropogenic forcings, although some regional features might be easier to explain

The surface impacts of Arctic stratospheric ozone anomalies

In the Arctic stratosphere, total column ozone in the spring can vary, from year to year, by as much as 30%. This large interannual variability, however, is absent from many present-generation climate models, in which the prescribed seasonal cycle of stratospheric ozone includes, at best, smooth multi-decadal trends. We here investigate the extent to which interannual variability in Arctic stratospheric ozone is able to affect the surface climate of the Northern Hemisphere extratropics. We do this by contrasting pairs of model integrations with positive and negative springtime ozone anomalies, using a simple yet widely used climate model. For ozone anomaly amplitudes somewhat larger than the recent observed variability, we find a significant influence on the tropospheric circulation, and the surface temperatures and precipitation patterns. More interestingly, these impacts have very clear regional patterns—they are largest over the North Atlantic sector—even though the prescribed ozone anomalies are zonally symmetric. However, confirming other studies, for ozone anomaly amplitudes within the observed range of the last three decades, our model experiments do not show statistically significant impacts at the surface

Stability analysis of a class of two-dimensional multipolar vortex equilibria

Published versio

An Efficient Spectral Dynamical Core for Distributed Memory Computers

The practical question of whether the classical spectral transform method, widely used in atmospheric modeling, can be efficiently implemented on inexpensive commodity clusters is addressed. Typically, such clusters have limited cache and memory sizes. To demonstrate that these limitations can be overcome, the authors have built a spherical general circulation model dynamical core, called BOB (“Built on Beowulf”), which can solve either the shallow water equations or the atmospheric primitive equations in pressure coordinates. That BOB is targeted for computing at high resolution on modestly sized and priced commodity clusters is reflected in four areas of its design. First, the associated Legendre polynomials (ALPs) are computed “on the fly” using a stable and accurate recursion relation. Second, an identity is employed that eliminates the storage of the derivatives of the ALPs. Both of these algorithmic choices reduce the memory footprint and memory bandwidth requirements of the spectral transform. Third, a cache-blocked and unrolled Legendre transform achieves a high performance level that resists deterioration as resolution is increased. Finally, the parallel implementation of BOB is transposition-based, employing load-balanced, one-dimensional decompositions in both latitude and wavenumber. A number of standard tests is used to compare BOB's performance to two well-known codes—the Parallel Spectral Transform Shallow Water Model (PSTSWM) and the dynamical core of NCAR's Community Climate Model CCM3. Compared to PSTSWM, BOB shows better timing results, particularly at the higher resolutions where cache effects become important. BOB also shows better performance in its comparison with CCM3's dynamical core. With 16 processors, at a triangular spectral truncation of T85, it is roughly five times faster when computing the solution to the standard Held–Suarez test case, which involves 18 levels in the vertical. BOB also shows a significantly smaller memory footprint in these comparison tests

The Antarctic Atmospheric Energy Budget. Part II: The Effect of Ozone Depletion and its Projected Recovery

In this study we continue our investigation of the atmospheric energy budget of the Antarctic polar cap (the region poleward of 70°S) using integrations of the Whole Atmosphere Community Climate Model from the year 1960 to 2065. In agreement with observational data, we find that the climatological mean net top-of-atmosphere (TOA) radiative flux is primarily balanced by the horizontal energy flux convergence over the polar cap. On interannual timescales, changes in the net TOA radiative flux are also primarily balanced by changes in the energy flux convergence, with the variability in both terms significantly correlated with the Southern Annular Mode (SAM); positive and negative correlations, respectively. On multidecadal timescales, twentieth century stratospheric ozone depletion produces a negative trend in the net TOA radiative flux due to a decrease in the absorbed solar radiation within the atmosphere-surface column. The negative trend in the net TOA radiative flux is balanced by a positive trend in energy flux convergence, primarily in austral summer. This negative (positive) trend in the net TOA radiation (energy flux convergence) occurs despite a positive trend in the SAM, suggesting that the effects of the SAM on the energy budget are overwhelmed by the direct radiative effects of ozone depletion. In the twenty-first century, ozone recovery is expected to reverse the negative trend in the net TOA radiative flux, which would then, again, be balanced by a decrease in the energy flux convergence. Therefore, over the next several decades, ozone recovery will, in all likelihood, mask the effect of GHG warming on the Antarctic energy budget