Multimodel climate and variability of the stratosphere

The stratospheric climate and variability from simulations of sixteen chemistry‐climate models is evaluated. On average the polar night jet is well reproduced though its variability is less well reproduced with a large spread between models. Polar temperature biases are less than 5 K except in the Southern Hemisphere (SH) lower stratosphere in spring. The accumulated area of low temperatures responsible for polar stratospheric cloud formation is accurately reproduced for the Antarctic but underestimated for the Arctic. The shape and position of the polar vortex is well simulated, as is the tropical upwelling in the lower stratosphere. There is a wide model spread in the frequency of major sudden stratospheric warnings (SSWs), late biases in the breakup of the SH vortex, and a weak annual cycle in the zonal wind in the tropical upper stratosphere. Quantitatively, “metrics” indicate a wide spread in model performance for most diagnostics with systematic biases in many, and poorer performance in the SH than in the Northern Hemisphere (NH). Correlations were found in the SH between errors in the final warming, polar temperatures, the leading mode of variability, and jet strength, and in the NH between errors in polar temperatures, frequency of major SSWs, and jet strength. Models with a stronger QBO have stronger tropical upwelling and a colder NH vortex. Both the qualitative and quantitative analysis indicate a number of common and long‐standing model problems, particularly related to the simulation of the SH and stratospheric variability.

How is open solar magnetic flux lost over the solar cycle?

The Sun's open magnetic field, magnetic flux dragged out into the heliosphere by the solar wind, varies by approximately a factor of 2 over the solar cycle. We consider the evolution of open solar flux in terms of a source and loss term. Open solar flux creation is likely to proceed at a rate dependent on the rate of photospheric flux emergence, which can be roughly parameterized by sunspot number or coronal mass ejection rate, when available. The open solar flux loss term is more difficult to relate to an observable parameter. The supersonic nature of the solar wind means open solar flux can only be removed by near-Sun magnetic reconnection between open solar magnetic field lines, be they open or closed heliospheric field lines. In this study we reconstruct open solar flux over the last three solar cycles and demonstrate that the loss term may be related to the degree to which the heliospheric current sheet (HCS) is warped, i.e., locally tilted from the solar rotation direction. This can account for both the large dip in open solar flux at the time of sunspot maximum as well as the asymmetry in open solar flux during the rising and declining phases of the solar cycle. The observed cycle-to-cycle variability is also well matched. Following Sheeley et al. (2001), we attribute modulation of open solar flux by the degree of warp of the HCS to the rate at which opposite polarity open solar flux is brought together by differential rotation.

Seasonality and interannual variability of the westerly jet in the Tibetan Plateau region

In this study, 40-yr ECMWF Re-Analysis (ERA-40) data are used for the description of the seasonal cycle and the interannual variability of the westerly jet in the Tibetan Plateau region. To complement results based on the analysis of monthly mean horizontal wind speeds, an occurrence-based jet climatology is constructed by identifying the locations of the jet axes at 6-hourly intervals throughout 1958–2001. Thus, a dataset describing the highly transient and localized features of jet variability is obtained. During winter and summer the westerly jet is located, respectively, to the south and north of the Tibetan Plateau. During the spring and autumn seasons there are jet transitions from south to north and vice versa. The median dates for these transitions are 28 April and 12 October. The spring transition is associated with large interannual variations, while the fall transition occurs more reliably within a 3-week period. The strength of the jet exhibits a peculiar seasonal cycle. During northward migration in April/May, the jet intensity weakens and its latitudinal position varies largely. In some springs, there are several transitions and split configurations occur before the jet settles in its northern summer position. In June, a well-defined and unusually strong jet reappears at the northern flanks of the Tibetan Plateau. In autumn, the jet gradually but reliably recedes to the south and is typically more intense than in spring. The jet transitions between the two preferred locations follow the seasonal latitudinal migration of the jet in the Northern Hemisphere. An analysis of interannual variations shows the statistical relationship between the strength of the summer jet, the tropospheric meridional temperature gradient, and the all-India rainfall series. Both this analysis and results from previous studies point to the particular dynamical relevance of the onsetting Indian summer monsoon precipitation and the associated diabatic heating for the formation of the strong summer jet. Finally, an example is provided that illustrates the climatological significance of the jet in terms of the covariation between the jet location and the spatial precipitation distribution in central Asia.

The role of atmosphere feedbacks during ENSO in the CMIP3 models. Part II: using AMIP runs to understand the heat flux feedback mechanisms

Several studies using ocean–atmosphere general circulation models (GCMs) suggest that the atmospheric component plays a dominant role in the modelled El Niño-Southern Oscillation (ENSO). To help elucidate these findings, the two main atmosphere feedbacks relevant to ENSO, the Bjerknes positive feedback (μ) and the heat flux negative feedback (α), are here analysed in nine AMIP runs of the CMIP3 multimodel dataset. We find that these models generally have improved feedbacks compared to the coupled runs which were analysed in part I of this study. The Bjerknes feedback, μ, is increased in most AMIP runs compared to the coupled run counterparts, and exhibits both positive and negative biases with respect to ERA40. As in the coupled runs, the shortwave and latent heat flux feedbacks are the two dominant components of α in the AMIP runs. We investigate the mechanisms behind these two important feedbacks, in particular focusing on the strong 1997–1998 El Niño. Biases in the shortwave flux feedback, α SW, are the main source of model uncertainty in α. Most models do not successfully represent the negative αSW in the East Pacific, primarily due to an overly strong low-cloud positive feedback in the far eastern Pacific. Biases in the cloud response to dynamical changes dominate the modelled α SW biases, though errors in the large-scale circulation response to sea surface temperature (SST) forcing also play a role. Analysis of the cloud radiative forcing in the East Pacific reveals model biases in low cloud amount and optical thickness which may affect α SW. We further show that the negative latent heat flux feedback, α LH, exhibits less diversity than α SW and is primarily driven by variations in the near-surface specific humidity difference. However, biases in both the near-surface wind speed and humidity response to SST forcing can explain the inter-model αLH differences.

An analysis of long-term trends, seasonality and short-term dynamics in water quality data from Plynlimon, Wales

This paper examines two hydrochemical time-series derived from stream samples taken in the Upper Hafren catchment, Plynlimon, Wales. One time-series comprises data collected at 7-hour intervals over 22 months (Neal et al., submitted, this issue), while the other is based on weekly sampling over 20 years. A subset of determinands: aluminium, calcium, chloride, conductivity, dissolved organic carbon, iron, nitrate, pH, silicon and sulphate are examined within a framework of non-stationary time-series analysis to identify determinand trends, seasonality and short-term dynamics. The results demonstrate that both long-term and high-frequency monitoring provide valuable and unique insights into the hydrochemistry of a catchment. The long-term data allowed analysis of long-termtrends, demonstrating continued increases in DOC concentrations accompanied by declining SO4 concentrations within the stream, and provided new insights into the changing amplitude and phase of the seasonality of the determinands such as DOC and Al. Additionally, these data proved invaluable for placing the short-term variability demonstrated within the high-frequency data within context. The 7-hour data highlighted complex diurnal cycles for NO3, Ca and Fe with cycles displaying changes in phase and amplitude on a seasonal basis. The high-frequency data also demonstrated the need to consider the impact that the time of sample collection can have on the summary statistics of the data and also that sampling during the hours of darkness provides additional hydrochemical information for determinands which exhibit pronounced diurnal variability. Moving forward, this research demonstrates the need for both long-term and high-frequency monitoring to facilitate a full and accurate understanding of catchment hydrochemical dynamics.

Implications of non-cylindrical flux ropes for magnetic cloud reconstruction techniques and the interpretation of double flux rope events

Magnetic clouds (MCs) are a subset of interplanetary coronal mass ejections (ICMEs) which exhibit signatures consistent with a magnetic flux rope structure. Techniques for reconstructing flux rope orientation from single-point in situ observations typically assume the flux rope is locally cylindrical, e.g., minimum variance analysis (MVA) and force-free flux rope (FFFR) fitting. In this study, we outline a non-cylindrical magnetic flux rope model, in which the flux rope radius and axial curvature can both vary along the length of the axis. This model is not necessarily intended to represent the global structure of MCs, but it can be used to quantify the error in MC reconstruction resulting from the cylindrical approximation. When the local flux rope axis is approximately perpendicular to the heliocentric radial direction, which is also the effective spacecraft trajectory through a magnetic cloud, the error in using cylindrical reconstruction methods is relatively small (≈ 10∘). However, as the local axis orientation becomes increasingly aligned with the radial direction, the spacecraft trajectory may pass close to the axis at two separate locations. This results in a magnetic field time series which deviates significantly from encounters with a force-free flux rope, and consequently the error in the axis orientation derived from cylindrical reconstructions can be as much as 90∘. Such two-axis encounters can result in an apparent ‘double flux rope’ signature in the magnetic field time series, sometimes observed in spacecraft data. Analysing each axis encounter independently produces reasonably accurate axis orientations with MVA, but larger errors with FFFR fitting.

Office rent determinants: a hedonic panel analysis

It has been frequently observed that office markets are subject to particularly high fluctuations in rents and vacancy levels, thus exposing real estate investors to considerable risk regarding expected future income streams. This paper analyzes the determinants of office rents and their variability over time and across sub-markets to gain insight into the rent price formation and its stability across space and over time. No support is found for the single-market hypothesis which states that arbitrage opportunities effectively align real estate pricing schemes in various parts of city. Instead, the results suggest that the importance of hedonic pricing factors varies both over time and across submarkets.

Symmetry breaking, mixing, instability, and low frequency variability in a minimal Lorenz-like system

Starting from the classical Saltzman two-dimensional convection equations, we derive via a severe spectral truncation a minimal 10 ODE system which includes the thermal effect of viscous dissipation. Neglecting this process leads to a dynamical system which includes a decoupled generalized Lorenz system. The consideration of this process breaks an important symmetry and couples the dynamics of fast and slow variables, with the ensuing modifications to the structural properties of the attractor and of the spectral features. When the relevant nondimensional number (Eckert number Ec) is different from zero, an additional time scale of O(Ec−1) is introduced in the system, as shown with standard multiscale analysis and made clear by several numerical evidences. Moreover, the system is ergodic and hyperbolic, the slow variables feature long-term memory with 1/f3/2 power spectra, and the fast variables feature amplitude modulation. Increasing the strength of the thermal-viscous feedback has a stabilizing effect, as both the metric entropy and the Kaplan-Yorke attractor dimension decrease monotonically with Ec. The analyzed system features very rich dynamics: it overcomes some of the limitations of the Lorenz system and might have prototypical value in relevant processes in complex systems dynamics, such as the interaction between slow and fast variables, the presence of long-term memory, and the associated extreme value statistics. This analysis shows how neglecting the coupling of slow and fast variables only on the basis of scale analysis can be catastrophic. In fact, this leads to spurious invariances that affect essential dynamical properties (ergodicity, hyperbolicity) and that cause the model losing ability in describing intrinsically multiscale processes.

Quantitative analysis of mannitol polymorphs: FT-Raman spectroscopy

Mannitol is a polymorphic excipient which is usually used in pharmaceutical products as the beta form, although other polymorphs (alpha and delta) are common contaminants. Binary mixtures containing beta and delta mannitol were prepared to quantify the concentration of the beta form using FT-Raman spectroscopy. Spectral regions characteristic of each form were selected and peak intensity ratios of beta peaks to delta peaks were calculated. Using these ratios, a correlation curve was established which was then validated by analysing further samples of known composition. The results indicate that levels down to 2% beta could be quantified using this novel, non-destructive approach. Potential errors associated with quantitative studies using FT-Raman spectroscopy were also researched. The principal source of variability arose from inhomogeneities on mixing of the samples; a significant reduction of these errors was observed by reducing and controlling the particle size range. The results show that FT-Raman spectroscopy can be used to rapidly and accurately quantitate polymorphic mixtures.

Inter-hemispheric mesospheric coupling in a comprehensive middle atmosphere model

Observations of noctilucent clouds have revealed a surprising coupling between the winter stratosphere and the summer polar mesopause region. In spite of the great distance involved, this inter-hemispheric link has been suggested to be the principal reason for both the year-to-year variability and the hemispheric differences in the frequency of occurrence of these high-altitude clouds. In this study, we investigate the dynamical influence of the winter stratosphere on the summer mesosphere using simulations from the vertically extended version of the Canadian Middle Atmosphere Model (CMAM). We find that for both Northern and Southern Hemispheres, variability in the summer polar mesopause region from one year to another can be traced back to the planetary-wave flux entering the winter stratosphere. The teleconnection pattern is the same for both positive and negative wave-flux anomalies. Using a composite analysis to isolate the events, it is argued that the mechanism for interhemispheric coupling is a feedback between summer mesosphere gravity-wave drag (GWD) and zonal wind, which is induced by an anomaly in mesospheric cross-equatorial flow, the latter arising from the anomaly in winter hemisphere GWD induced by the anomaly in stratospheric conditions.

Climate variability-observations, reconstructions, and model simulations for the Atlantic-European and Alpine region from 1500-2100 AD

A detailed analysis is undertaken of the Atlantic-European climate using data from 500-year-long proxy-based climate reconstructions, a long climate simulation with perpetual 1990 forcing, as well as two global and one regional climate change scenarios. The observed and simulated interannual variability and teleconnectivity are compared and interpreted in order to improve the understanding of natural climate variability on interannual to decadal time scales for the late Holocene. The focus is set on the Atlantic-European and Alpine regions during the winter and summer seasons, using temperature, precipitation, and 500 hPa geopotential height fields. The climate reconstruction shows pronounced interdecadal variations that appear to “lock” the atmospheric circulation in quasi-steady long-term patterns over multi-decadal periods controlling at least part of the temperature and precipitation variability. Different circulation patterns are persistent over several decades for the period 1500 to 1900. The 500-year-long simulation with perpetual 1990 forcing shows some substantial differences, with a more unsteady teleconnectivity behaviour. Two global scenario simulations indicate a transition towards more stable teleconnectivity for the next 100 years. Time series of reconstructed and simulated temperature and precipitation over the Alpine region show comparatively small changes in interannual variability within the time frame considered, with the exception of the summer season, where a substantial increase in interannual variability is simulated by regional climate models.

The effect of climate change on the variability of the Northern Hemisphere stratospheric polar vortex

With extreme variability of the Arctic polar vortex being a key link for stratosphere–troposphere influences, its evolution into the twenty-first century is important for projections of changing surface climate in response to greenhouse gases. Variability of the stratospheric vortex is examined using a state-of-the-art climate model and a suite of specifically developed vortex diagnostics. The model has a fully coupled ocean and a fully resolved stratosphere. Analysis of the standard stratospheric zonal mean wind diagnostic shows no significant increase over the twenty-first century in the number of major sudden stratospheric warmings (SSWs) from its historical value of 0.7 events per decade, although the monthly distribution of SSWs does vary, with events becoming more evenly dispersed throughout the winter. However, further analyses using geometric-based vortex diagnostics show that the vortex mean state becomes weaker, and the vortex centroid is climatologically more equatorward by up to 2.5°, especially during early winter. The results using these diagnostics not only characterize the vortex structure and evolution but also emphasize the need for vortex-centric diagnostics over zonally averaged measures. Finally, vortex variability is subdivided into wave-1 (displaced) and -2 (split) components, and it is implied that vortex displacement events increase in frequency under climate change, whereas little change is observed in splitting events.

The role of atmosphere feedbacks during ENSO in the CMIP3 models. Part III: the shortwave flux feedback

Previous studies using coupled general circulation models (GCMs) suggest that the atmosphere model plays a dominant role in the modeled El Nin ̃ o–Southern Oscillation (ENSO), and that intermodel differences in the thermodynamical damping of sea surface temperatures (SSTs) are a dominant contributor to the ENSO amplitude diversity. This study presents a detailed analysis of the shortwave flux feedback (aSW) in 12 Coupled Model Intercomparison Project phase 3 (CMIP3) simulations, motivated by findings that aSW is the primary contributor to model thermodynamical damping errors. A ‘‘feedback decomposition method,’’ developed to elucidate the aSW biases, shows that all models un- derestimate the dynamical atmospheric response to SSTs in the eastern equatorial Pacific, leading to un- derestimated aSW values. Biases in the cloud response to dynamics and the shortwave interception by clouds also contribute to errors in aSW. Changes in the aSW feedback between the coupled and corresponding atmosphere-only simulations are related to changes in the mean dynamics. A large nonlinearity is found in the observed and modeled SW flux feedback, hidden when linearly cal- culating aSW. In the observations, two physical mechanisms are proposed to explain this nonlinearity: 1) a weaker subsidence response to cold SST anomalies than the ascent response to warm SST anomalies and 2) a nonlinear high-level cloud cover response to SST. The shortwave flux feedback nonlinearity tends to be underestimated by the models, linked to an underestimated nonlinearity in the dynamical response to SST. The process-based methodology presented in this study may help to correct model ENSO atmospheric biases, ultimately leading to an improved simulation of ENSO in GCMs.

On the drivers of inter-annual and decadal rainfall variability in Queensland, Australia

Queensland experiences considerable inter-annual and decadal rainfall variability, which impacts water-resource management, agriculture and infrastructure. To understand the mechanisms by which large-scale atmospheric and coupled air–sea processes drive these variations, empirical orthogonal teleconnection (EOT) analysis is applied to 1900–2010 seasonal Queensland rainfall. Fields from observations and the 20th Century Reanalysis are regressed onto the EOT timeseries to associate the EOTs with large-scale drivers. In winter, spring and summer the leading, state-wide EOTs are highly correlated with the El Nino–Southern Oscillation (ENSO); the Inter-decadal Pacific Oscillation modulates the summer ENSO teleconnection. In autumn, the leading EOT is associated with locally driven, late-season monsoon variations, while ENSO affects only tropical northern Queensland. Examining EOTs beyond the first, southeastern Queensland and the Cape York peninsula emerge as regions of coherent rainfall variability. In the southeast, rainfall anomalies respond to the strength and moisture content of onshore easterlies, controlled by Tasman Sea blocking. The summer EOT associated with onshore flow and blocking has been negative since 1970, consistent with the observed decline in rainfall along the heavily populated coast. The southeastern Queensland EOTs show considerable multi-decadal variability, which is independent of large-scale drivers. Summer rainfall in Cape York is associated with tropical-cyclone activity.

Decadal to multidecadal variability and the climate change background

Three prominent quasi-global patterns of variability and change are observed using the Met Office's sea surface temperature (SST) analysis and almost independent night marine air temperature analysis. The first is a global warming signal that is very highly correlated with global mean SST. The second is a decadal to multidecadal fluctuation with some geographical similarity to the El Niño–Southern Oscillation (ENSO). It is associated with the Pacific Decadal Oscillation (PDO), and its Pacific-wide manifestation has been termed the Interdecadal Pacific Oscillation (IPO). We present model investigations of the relationship between the IPO and ENSO. The third mode is an interhemispheric variation on multidecadal timescales which, in view of climate model experiments, is likely to be at least partly due to natural variations in the thermohaline circulation. Observed climatic impacts of this mode also appear in model simulations. Smaller-scale, regional atmospheric phenomena also affect climate on decadal to interdecadal timescales. We concentrate on one such mode, the winter North Atlantic Oscillation (NAO). This shows strong decadal to interdecadal variability and a correspondingly strong influence on surface climate variability which is largely additional to the effects of recent regional anthropogenic climate change. The winter NAO is likely influenced by both SST forcing and stratospheric variability. A full understanding of decadal changes in the NAO and European winter climate may require a detailed representation of the stratosphere that is hitherto missing in the major climate models used to study climate change.