The impacts of climate change on river flow regimes at the global scale

This paper presents an assessment of the impacts of climate change on a series of indicators of hydrological regimes across the global domain, using a global hydrological model run with climate scenarios constructed using pattern-scaling from 21 CMIP3 (Coupled Model Intercomparison Project Phase 3) climate models. Changes are compared with natural variability, with a significant change being defined as greater than the standard deviation of the hydrological indicator in the absence of climate change. Under an SRES (Special Report on Emissions Scenarios) A1b emissions scenario, substantial proportions of the land surface (excluding Greenland and Antarctica) would experience significant changes in hydrological behaviour by 2050; under one climate model scenario (Hadley Centre HadCM3), average annual runoff increases significantly over 47% of the land surface and decreases over 36%; only 17% therefore sees no significant change. There is considerable variability between regions, depending largely on projected changes in precipitation. Uncertainty in projected river flow regimes is dominated by variation in the spatial patterns of climate change between climate models (hydrological model uncertainty is not included). There is, however, a strong degree of consistency in the overall magnitude and direction of change. More than two-thirds of climate models project a significant increase in average annual runoff across almost a quarter of the land surface, and a significant decrease over 14%, with considerably higher degrees of consistency in some regions. Most climate models project increases in runoff in Canada and high-latitude eastern Europe and Siberia, and decreases in runoff in central Europe, around the Mediterranean, the Mashriq, central America and Brasil. There is some evidence that projecte change in runoff at the regional scale is not linear with change in global average temperature change. The effects of uncertainty in the rate of future emissions is relatively small

Multimodel assessment of the factors driving stratospheric ozone evolution over the 21st century

The evolution of stratospheric ozone from 1960 to 2100 is examined in simulations from 14 chemistry‐climate models, driven by prescribed levels of halogens and greenhouse gases. There is general agreement among the models that total column ozone reached a minimum around year 2000 at all latitudes, projected to be followed by an increase over the first half of the 21st century. In the second half of the 21st century, ozone is projected to continue increasing, level off, or even decrease depending on the latitude. Separation into partial columns above and below 20 hPa reveals that these latitudinal differences are almost completely caused by differences in the model projections of ozone in the lower stratosphere. At all latitudes, upper stratospheric ozone increases throughout the 21st century and is projected to return to 1960 levels well before the end of the century, although there is a spread among models in the dates that ozone returns to specific historical values. We find decreasing halogens and declining upper atmospheric temperatures, driven by increasing greenhouse gases, contribute almost equally to increases in upper stratospheric ozone. In the tropical lower stratosphere, an increase in upwelling causes a steady decrease in ozone through the 21st century, and total column ozone does not return to 1960 levels in most of the models. In contrast, lower stratospheric and total column ozone in middle and high latitudes increases during the 21st century, returning to 1960 levels well before the end of the century in most models.

Estimation of optimal gravity wave parameters for climate models using data assimilation

There is a current need to constrain the parameters of gravity wave drag (GWD) schemes in climate models using observational information instead of tuning them subjectively. In this work, an inverse technique is developed using data assimilation principles to estimate gravity wave parameters. Because mostGWDschemes assume instantaneous vertical propagation of gravity waves within a column, observations in a single column can be used to formulate a one-dimensional assimilation problem to estimate the unknown parameters. We define a cost function that measures the differences between the unresolved drag inferred from observations (referred to here as the ‘observed’ GWD) and the GWD calculated with a parametrisation scheme. The geometry of the cost function presents some difficulties, including multiple minima and ill-conditioning because of the non-independence of the gravity wave parameters. To overcome these difficulties we propose a genetic algorithm to minimize the cost function, which provides a robust parameter estimation over a broad range of prescribed ‘true’ parameters. When real experiments using an independent estimate of the ‘observed’ GWD are performed, physically unrealistic values of the parameters can result due to the non-independence of the parameters. However, by constraining one of the parameters to lie within a physically realistic range, this degeneracy is broken and the other parameters are also found to lie within physically realistic ranges. This argues for the essential physical self-consistency of the gravity wave scheme. A much better fit to the observed GWD at high latitudes is obtained when the parameters are allowed to vary with latitude. However, a close fit can be obtained either in the upper or the lower part of the profiles, but not in both at the same time. This result is a consequence of assuming an isotropic launch spectrum. The changes of sign in theGWDfound in the tropical lower stratosphere, which are associated with part of the quasi-biennial oscillation forcing, cannot be captured by the parametrisation with optimal parameters.

Summertime stratospheric processes at northern mid-latitudes: comparisons between MANTRA balloon measurements and the Canadian Middle Atmosphere Model

In this paper we report on a study conducted using the Middle Atmospheric Nitrogen TRend Assessment (MANTRA) balloon measurements of stratospheric constituents and temperature and the Canadian Middle Atmosphere Model (CMAM). Three different kinds of data are used to assess the inter-consistency of the combined dataset: single profiles of long-lived species from MANTRA 1998, sparse climatologies from the ozonesonde measurements during the four MANTRA campaigns and from HALOE satellite measurements, and the CMAM climatology. In doing so, we evaluate the ability of the model to reproduce the measured fields and to thereby test our ability to describe mid-latitude summertime stratospheric processes. The MANTRA campaigns were conducted at Vanscoy, Saskatchewan, Canada (52◦ N, 107◦ W)in late August and early September of 1998, 2000, 2002 and 2004. During late summer at mid-latitudes, the stratosphere is close to photochemical control, providing an ideal scenario for the study reported here. From this analysis we find that: (1) reducing the value for the vertical diffusion coefficient in CMAM to a more physically reasonable value results in the model better reproducing the measured profiles of long-lived species; (2) the existence of compact correlations among the constituents, as expected from independent measurements in the literature and from models, confirms the self-consistency of the MANTRA measurements; and (3) the 1998 measurements show structures in the chemical species profiles that can be associated with transport, adding to the growing evidence that the summertime stratosphere can be much more disturbed than anticipated. The mechanisms responsible for such disturbances need to be understood in order to assess the representativeness of the measurements and to isolate longterm trends.

Symmetric stability of compressible zonal flows on a generalized equatorial β plane

Sufficient conditions are derived for the linear stability with respect to zonally symmetric perturbations of a steady zonal solution to the nonhydrostatic compressible Euler equations on an equatorial � plane, including a leading order representation of the Coriolis force terms due to the poleward component of the planetary rotation vector. A version of the energy–Casimir method of stability proof is applied: an invariant functional of the Euler equations linearized about the equilibrium zonal flow is found, and positive definiteness of the functional is shown to imply linear stability of the equilibrium. It is shown that an equilibrium is stable if the potential vorticity has the same sign as latitude and the Rayleigh centrifugal stability condition that absolute angular momentum increase toward the equator on surfaces of constant pressure is satisfied. The result generalizes earlier results for hydrostatic and incompressible systems and for systems that do not account for the nontraditional Coriolis force terms. The stability of particular equilibrium zonal velocity, entropy, and density fields is assessed. A notable case in which the effect of the nontraditional Coriolis force is decisive is the instability of an angular momentum profile that decreases away from the equator but is flatter than quadratic in latitude, despite its satisfying both the centrifugal and convective stability conditions.

Assessment of temperature, trace species, and ozone in chemistry-climate model simulations of the recent past

Simulations of the stratosphere from thirteen coupled chemistry-climate models (CCMs) are evaluated to provide guidance for the interpretation of ozone predictions made by the same CCMs. The focus of the evaluation is on how well the fields and processes that are important for determining the ozone distribution are represented in the simulations of the recent past. The core period of the evaluation is from 1980 to 1999 but long-term trends are compared for an extended period (1960–2004). Comparisons of polar high-latitude temperatures show that most CCMs have only small biases in the Northern Hemisphere in winter and spring, but still have cold biases in the Southern Hemisphere spring below 10 hPa. Most CCMs display the correct stratospheric response of polar temperatures to wave forcing in the Northern, but not in the Southern Hemisphere. Global long-term stratospheric temperature trends are in reasonable agreement with satellite and radiosonde observations. Comparisons of simulations of methane, mean age of air, and propagation of the annual cycle in water vapor show a wide spread in the results, indicating differences in transport. However, for around half the models there is reasonable agreement with observations. In these models the mean age of air and the water vapor tape recorder signal are generally better than reported in previous model intercomparisons. Comparisons of the water vapor and inorganic chlorine (Cly) fields also show a large intermodel spread. Differences in tropical water vapor mixing ratios in the lower stratosphere are primarily related to biases in the simulated tropical tropopause temperatures and not transport. The spread in Cly, which is largest in the polar lower stratosphere, appears to be primarily related to transport differences. In general the amplitude and phase of the annual cycle in total ozone is well simulated apart from the southern high latitudes. Most CCMs show reasonable agreement with observed total ozone trends and variability on a global scale, but a greater spread in the ozone trends in polar regions in spring, especially in the Arctic. In conclusion, despite the wide range of skills in representing different processes assessed here, there is sufficient agreement between the majority of the CCMs and the observations that some confidence can be placed in their predictions.

Seasonal persistence of midlatitude total ozone anomalies

Temporal autocorrelations of monthly mean total ozone anomalies over the 35–60°S and 35–60°N latitude bands reveal that anomalies established in the wintertime midlatitude ozone buildup persist (with photochemical decay) until the end of the following autumn, and then are rapidly erased once the next winter's buildup begins. The photochemical decay rate is found to be identical between the two hemispheres. High predictability of ozone through late summer exists based on the late-spring values. In the northern hemisphere, extending the 1979–2001 springtime ozone trend to other months through regression based on the seasonal persistence of anomalies captures the seasonality of the ozone trends remarkably well. In the southern hemisphere, the springtime trend only accounts for part of the summertime trends. There is a strong correlation between the ozone anomalies in northern hemisphere spring and those in the subsequent southern hemisphere spring, but not the converse.

Wave-breaking characteristics of Northern Hemisphere winter blocking: a two-dimensional approach

This paper generalises and applies recently developed blocking diagnostics in a two- dimensional latitude-longitude context, which takes into consideration both mid- and high-latitude blocking. These diagnostics identify characteristics of the associated wave-breaking as seen in the potential temperature (θ) on the dynamical tropopause, in particular the cyclonic or anticyclonic Direction of wave-Breaking (DB index), and the Relative Intensity (RI index) of the air masses that contribute to blocking formation. The methodology is extended to a 2-D domain and a cluster technique is deployed to classify mid- and high-latitude blocking according to the wave-breaking characteristics. Mid-latitude blocking is observed over Europe and Asia, where the meridional gradient of θ is generally weak, whereas high-latitude blocking is mainly present over the oceans, to the north of the jet-stream, where the meridional gradient of θ is much stronger. They occur respectively on the equatorward and poleward flank of the jet- stream, where the horizontal shear ∂u/∂y is positive in the first case and negative in the second case. A regional analysis is also conducted. It is found that cold-anticyclonic and cyclonic blocking divert the storm-track respectively to the south and to the north over the East Atlantic and western Europe. Furthermore, warm-cyclonic blocking over the Pacific and cold-anticyclonic blocking over Europe are identified as the most persistent types and are associated with large amplitude anomalies in temperature and precipitation. Finally, the high-latitude, cyclonic events seem to correlate well with low- frequency modes of variability over the Pacific and Atlantic Ocean.

Winter and Summer Northern Hemisphere blocking in CMIP5 models

The frequencies of atmospheric blocking in both winter and summer and the changes in them from the 20th to the 21st centuries as simulated in twelve CMIP5 models is analysed. The RCP 8.5 high emission scenario runs are used to represent the 21st century. The analysis is based on the wave-breaking methodology of Pelly and Hoskins (2003a). It differs from the Tibaldi and Molteni (1990) index in viewing equatorward cut-off lows and poleward blocking highs in equal manner as indicating a disruption to the westerlies. 1-dimensional and 2-dimensional diagnostics are applied to identify blocking of the mid-latitude storm-track and also at higher latitudes. Winter blocking frequency is found to be generally underestimated. The models give a decrease in the European blocking maximum in the 21st century, consistent with the results in other studies. There is a mean 21st century winter poleward shift of high- latitude blocking, but little agreement between the models on the details. In summer, Eurasian blocking is also underestimated in the models, whereas it is now too large over the high-latitude ocean basins. A decrease in European blocking frequency in the 21st century model runs is again found. However in summer there is a clear eastward shift of blocking over Eastern Europe and Western Russia, in a region close to the blocking that dominated the Russian summer of 2010. While summer blocking decreases in general, the poleward shift of the storm track into the region of frequent high latitude blocking may mean that the incidence of storms being obstructed by blocks may actually increase.

Climate variability in the southwest Pacific during the Last Termination (20-10 kyr BP)

The degree to which palaeoclimatic changes in the Southern Hemisphere co-varied with events in the high latitude Northern Hemisphere during the Last Termination is a contentious issue, with conflicting evidence for the degree of ‘teleconnection’ between different regions of the Southern Hemisphere. The available hypotheses are difficult to test robustly, however, because there are few detailed palaeoclimatic records in the Southern Hemisphere. Here we present climatic reconstructions from the southwestern Pacific, a key region in the Southern Hemisphere because of the potentially important role it plays in global climate change. The reconstructions for the period 20–10 kyr BP were obtained from five sites along a transect from southern New Zealand, through Australia to Indonesia, supported by 125 calibrated 14C ages. Two periods of significant climatic change can be identified across the region at around 17 and 14.2 cal kyr BP, most probably associated with the onset of warming in the West Pacific Warm Pool and the collapse of Antarctic ice during Meltwater Pulse-1A, respectively. The severe geochronological constraints that inherently afflict age models based on radiocarbon dating and the lack of quantified climatic parameters make more detailed interpretations problematic, however. There is an urgent need to address the geochronological limitations, and to develop more precise and quantified estimates of the pronounced climate variations that clearly affected this region during the Last Termination.

The role of baroclinic waves in the initiation of tropical cyclones across the southern Indian Ocean

Cases where tropical storms are initiated simultaneously along one latitude are investigated. It is argued that such structure arises as part of a baroclinic wave. A case from February 2008 is examined using European Centre for Medium-Range Forecasts (ECMWF) analyses; the birth of three tropical cyclones in the low-level cyclonic regions to the east of upper-level troughs suggests that the wave was instrumental for initiation. Archived satellite imagery and storm warnings reveal that baroclinic waves over the southern Indian Ocean accompany tropical cyclogenesis twice a season on average, mainly in late summer, when breaking Rossby waves on the subtropical westerly jet are closest to the Intertropical Convergence Zone (ITCZ). Copyright © 2012 Royal Meteorological Society

Are the winters 2010 and 2012 archetypes exhibiting extreme opposite behavior of the North Atlantic jet stream?

The atmospheric circulation over the North Atlantic-European sector experienced exceptional but highly contrasting conditions in the recent 2010 and 2012 winters (November-March; with the year dated by the relevant January). Evidence is given for the remarkably different locations of the eddy-driven westerly jet over the North Atlantic. In the 2010 winter the maximum of the jet stream was systematically between 30ºN and 40ºN (in the ‘south jet regime’), while in the 2012 winter it was predominantly located around 55ºN (north jet regime). These jet features underline the occurrence of either weak flow (2010) or strong and persistent ridges throughout the troposphere (2012). This is confirmed by the very different occurrence of blocking systems over the North Atlantic, associated with episodes of strong cyclonic (anticyclonic) Rossby wave breaking in 2010 (2012) winters. These dynamical features underlie strong precipitation and temperature anomalies over parts of Europe, with detrimental impacts on many socioeconomic sectors. Despite the highly contrasting atmospheric states, mid and high-latitude boundary conditions do not reveal strong differences in these two winters. The two winters were associated with opposite ENSO phases, but there is no causal evidence of a remote forcing from the Pacific sea surface temperatures. Finally, the exceptionality of the two winters is demonstrated in relation to the last 140 years. It is suggested that these winters may be seen as archetypes of North Atlantic jet variability under current climate conditions.

Are greenhouse gas signals of Northern Hemisphere winter extra-tropical cyclone activity dependent on the identification and tracking algorithm?

For Northern Hemisphere extra-tropical cyclone activity, the dependency of a potential anthropogenic climate change signal on the identification method applied is analysed. This study investigates the impact of the used algorithm on the changing signal, not the robustness of the climate change signal itself. Using one single transient AOGCM simulation as standard input for eleven state-of-the-art identification methods, the patterns of model simulated present day climatologies are found to be close to those computed from re-analysis, independent of the method applied. Although differences in the total number of cyclones identified exist, the climate change signals (IPCC SRES A1B) in the model run considered are largely similar between methods for all cyclones. Taking into account all tracks, decreasing numbers are found in the Mediterranean, the Arctic in the Barents and Greenland Seas, the mid-latitude Pacific and North America. Changing patterns are even more similar, if only the most severe systems are considered: the methods reveal a coherent statistically significant increase in frequency over the eastern North Atlantic and North Pacific. We found that the differences between the methods considered are largely due to the different role of weaker systems in the specific methods.

IMILAST: a community effort to intercompare extratropical cyclone detection and tracking algorithms

The variability of results from different automated methods of detection and tracking of extratropical cyclones is assessed in order to identify uncertainties related to the choice of method. Fifteen international teams applied their own algorithms to the same dataset—the period 1989–2009 of interim European Centre for Medium-Range Weather Forecasts (ECMWF) Re-Analysis (ERAInterim) data. This experiment is part of the community project Intercomparison of Mid Latitude Storm Diagnostics (IMILAST; see www.proclim.ch/imilast/index.html). The spread of results for cyclone frequency, intensity, life cycle, and track location is presented to illustrate the impact of using different methods. Globally, methods agree well for geographical distribution in large oceanic regions, interannual variability of cyclone numbers, geographical patterns of strong trends, and distribution shape for many life cycle characteristics. In contrast, the largest disparities exist for the total numbers of cyclones, the detection of weak cyclones, and distribution in some densely populated regions. Consistency between methods is better for strong cyclones than for shallow ones. Two case studies of relatively large, intense cyclones reveal that the identification of the most intense part of the life cycle of these events is robust between methods, but considerable differences exist during the development and the dissolution phases.

A combined statistical and dynamical approach for downscaling large-scale footprints of European windstorms

The occurrence of mid-latitude windstorms is related to strong socio-economic effects. For detailed and reliable regional impact studies, large datasets of high-resolution wind fields are required. In this study, a statistical downscaling approach in combination with dynamical downscaling is introduced to derive storm related gust speeds on a high-resolution grid over Europe. Multiple linear regression models are trained using reanalysis data and wind gusts from regional climate model simulations for a sample of 100 top ranking windstorm events. The method is computationally inexpensive and reproduces individual windstorm footprints adequately. Compared to observations, the results for Germany are at least as good as pure dynamical downscaling. This new tool can be easily applied to large ensembles of general circulation model simulations and thus contribute to a better understanding of the regional impact of windstorms based on decadal and climate change projections.