On the attribution of stratospheric ozone and temperature changes to changes in ozone-depleting substances and well-mixed greenhouse gases

The vertical profile of global-mean stratospheric temperature changes has traditionally represented an important diagnostic for the attribution of the cooling effects of stratospheric ozone depletion and CO2 increases. However, CO2-induced cooling alters ozone abundance by perturbing ozone chemistry, thereby coupling the stratospheric ozone and temperature responses to changes in CO2 and ozone-depleting substances (ODSs). Here we untangle the ozone-temperature coupling and show that the attribution of global-mean stratospheric temperature changes to CO2 and ODS changes (which are the true anthropogenic forcing agents) can be quite different from the traditional attribution to CO2 and ozone changes. The significance of these effects is quantified empirically using simulations from a three-dimensional chemistry-climate model. The results confirm the essential validity of the traditional approach in attributing changes during the past period of rapid ODS increases, although we find that about 10% of the upper stratospheric ozone decrease from ODS increases over the period 1975–1995 was offset by the increase in CO2, and the CO2-induced cooling in the upper stratosphere has been somewhat overestimated. When considering ozone recovery, however, the ozone-temperature coupling is a first-order effect; fully 2/5 of the upper stratospheric ozone increase projected to occur from 2010–2040 is attributable to CO2 increases. Thus, it has now become necessary to base attribution of global-mean stratospheric temperature changes on CO2 and ODS changes rather than on CO2 and ozone changes.

Conditioning model output statistics of regional climate model precipitation on circulation patterns

Dynamical downscaling of Global Climate Models (GCMs) through regional climate models (RCMs) potentially improves the usability of the output for hydrological impact studies. However, a further downscaling or interpolation of precipitation from RCMs is often needed to match the precipitation characteristics at the local scale. This study analysed three Model Output Statistics (MOS) techniques to adjust RCM precipitation; (1) a simple direct method (DM), (2) quantile-quantile mapping (QM) and (3) a distribution-based scaling (DBS) approach. The modelled precipitation was daily means from 16 RCMs driven by ERA40 reanalysis data over the 1961–2000 provided by the ENSEMBLES (ENSEMBLE-based Predictions of Climate Changes and their Impacts) project over a small catchment located in the Midlands, UK. All methods were conditioned on the entire time series, separate months and using an objective classification of Lamb's weather types. The performance of the MOS techniques were assessed regarding temporal and spatial characteristics of the precipitation fields, as well as modelled runoff using the HBV rainfall-runoff model. The results indicate that the DBS conditioned on classification patterns performed better than the other methods, however an ensemble approach in terms of both climate models and downscaling methods is recommended to account for uncertainties in the MOS methods.

Bifurcations and state changes in the human alpha rhythm: Theory and experiment

Despite many decades investigating scalp recordable 8–13-Hz (alpha) electroencephalographic activity, no consensus has yet emerged regarding its physiological origins nor its functional role in cognition. Here we outline a detailed, physiologically meaningful, theory for the genesis of this rhythm that may provide important clues to its functional role. In particular we find that electroencephalographically plausible model dynamics, obtained with physiological admissible parameterisations, reveals a cortex perched on the brink of stability, which when perturbed gives rise to a range of unanticipated complex dynamics that include 40-Hz (gamma) activity. Preliminary experimental evidence, involving the detection of weak nonlinearity in resting EEG using an extension of the well-known surrogate data method, suggests that nonlinear (deterministic) dynamics are more likely to be associated with weakly damped alpha activity. Thus rather than the “alpha rhythm” being an idling rhythm it may be more profitable to conceive it as a readiness rhythm.

Past and recent changes in the North Atlantic oscillation

The North Atlantic oscillation (NAO) is under current climate conditions the leading mode of atmospheric circulation variability over the North Atlantic region. While the pattern is present during the entire year, it is most important during winter, explaining a large part of the variability of the large-scale pressure field, being thus largely determinant for the weather conditions over the North Atlantic basin and over Western Europe. In this study, a review of recent literature on the basic understanding of the NAO, its variability on different time scales and driving physical mechanisms is presented. In particular, the observed NAO variations and long-term trends are put into a long term perspective by considering paleo-proxy evidence. A representative number of recently released NAO reconstructions are discussed. While the reconstructions agree reasonably well with observations during the instrumental overlapping period, there is a rather high uncertainty between the different reconstructions for the pre-instrumental period, which leads to partially incoherent results, that is, periods where the NAO reconstructions do not agree even in sign. Finally, we highlight the future need of a broader definition of the NAO, the assessment of the stability of the teleconnection centers over time, the analysis of the relations to other relevant variables like temperature and precipitation, as well as on the relevant processes involved

Optical characterization of the growth and spatial structure of a substorm onset arc

We present a detailed case study of the characteristics of auroral forms that constitute the first ionospheric signatures of substorm expansion phase onset. Analysis of the optical frequency and along-arc (azimuthal) wave number spectra provides the strongest constraint to date on the potential mechanisms and instabilities in the near-Earth magnetosphere that accompany auroral onset and which precede poleward arc expansion and auroral breakup. We evaluate the frequency and growth rates of the auroral forms as a function of azimuthal wave number to determine whether these wave characteristics are consistent with current models of the substorm onset mechanism. We find that the frequency, spatial scales, and growth rates of the auroral forms are most consistent with the cross-field current instability or a ballooning instability, most likely triggered close to the inner edge of the ion plasma sheet. This result is supportive of a near-Earth plasma sheet initiation of the substorm expansion phase. We also present evidence that the frequency and phase characteristics of the auroral undulations may be generated via resonant processes operating along the geomagnetic field. Our observations provide the most powerful constraint to date on the ionospheric manifestation of the physical processes operating during the first few minutes around auroral substorm onset.

Transport timescales and tracer properties in the extratropical UTLS

A comprehensive evaluation of seasonal backward trajectories initialized in the northern hemisphere lowermost stratosphere (LMS) has been performed to investigate the factors that determine the temporal and spatial structure of troposphere-to-stratosphere-transport (TST) and it's impact on the LMS. In particular we explain the fundamental role of the transit time since last TST (tTST) for the chemical composition of the LMS. According to our results the structure of the LMS can be characterized by a layer with tTST<40 days forming a narrow band around the local tropopause. This layer extends about 30 K above the local dynamical tropopause, corresponding to the extratropical tropopause transition layer (ExTL) as identified by CO. The LMS beyond this layer shows a relatively well defined separation as marked by an aprupt transition to longer tTST indicating less frequent mixing and a smaller fraction of tropospheric air. Thus the LMS constitutes a region of two well defined regimes of tropospheric influence. These can be characterized mainly by different transport times from the troposphere and different fractions of tropospheric air. Carbon monoxide (CO) mirrors this structure of tTST due to it's finite lifetime on the order of three months. Water vapour isopleths, on the other hand, do not uniquely indicate TST and are independent of tTST, but are determined by the Lagrangian Cold Point (LCP) of air parcels. Most of the backward trajectories from the LMS experienced their LCP in the tropics and sub-tropics, and TST often occurs 20 days after trajectories have encountered their LCP. Therefore, ExTL properties deduced from CO and H2O provide totally different informations on transport and particular TST for the LMS.

Complex spatiotemporal haemodynamic response following sensory stimulation in the awake rat

Detailed understanding of the haemodynamic changes that underlie non-invasive neuroimaging techniques such as blood oxygen level dependent functional magnetic resonance imaging is essential if we are to continue to extend the use of these methods for understanding brain function and dysfunction. The use of animal and in particular rodent research models has been central to these endeavours as they allow in-vivo experimental techniques that provide measurements of the haemodynamic response function at high temporal and spatial resolution. A limitation of most of this research is the use of anaesthetic agents which may disrupt or mask important features of neurovascular coupling or the haemodynamic response function. In this study we therefore measured spatiotemporal cortical haemodynamic responses to somatosensory stimulation in awake rats using optical imaging spectroscopy. Trained, restrained animals received non-noxious stimulation of the whisker pad via chronically implanted stimulating microwires whilst optical recordings were made from the contralateral somatosensory cortex through a thin cranial window. The responses we measure from un-anaesthetised animals are substantially different from those reported in previous studies which have used anaesthetised animals. These differences include biphasic response regions (initial increases in blood volume and oxygenation followed by subsequent decreases) as well as oscillations in the response time series of awake animals. These haemodynamic response features do not reflect concomitant changes in the underlying neuronal activity and therefore reflect neurovascular or cerebrovascular processes. These hitherto unreported hyperemic response dynamics may have important implications for the use of anaesthetised animal models for research into the haemodynamic response function.

Model estimation of cerebral hemodynamics between blood flow and volume changes: A data-based modeling approach

It is well known that there is a dynamic relationship between cerebral blood flow (CBF) and cerebral blood volume (CBV). With increasing applications of functional MRI, where the blood oxygen-level-dependent signals are recorded, the understanding and accurate modeling of the hemodynamic relationship between CBF and CBV becomes increasingly important. This study presents an empirical and data-based modeling framework for model identification from CBF and CBV experimental data. It is shown that the relationship between the changes in CBF and CBV can be described using a parsimonious autoregressive with exogenous input model structure. It is observed that neither the ordinary least-squares (LS) method nor the classical total least-squares (TLS) method can produce accurate estimates from the original noisy CBF and CBV data. A regularized total least-squares (RTLS) method is thus introduced and extended to solve such an error-in-the-variables problem. Quantitative results show that the RTLS method works very well on the noisy CBF and CBV data. Finally, a combination of RTLS with a filtering method can lead to a parsimonious but very effective model that can characterize the relationship between the changes in CBF and CBV.

Dissociating between object affordances and spatial compatibility effects using early response components

Perception and action are tightly linked: objects may be perceived not only in terms of visual features, but also in terms of possibilities for action. Previous studies showed that when a centrally located object has a salient graspable feature (e.g., a handle), it facilitates motor responses corresponding with the feature's position. However, such so-called affordance effects have been criticized as resulting from spatial compatibility effects, due to the visual asymmetry created by the graspable feature, irrespective of any affordances. In order to dissociate between affordance and spatial compatibility effects, we asked participants to perform a simple reaction-time task to typically graspable and non-graspable objects with similar visual features (e.g., lollipop and stop sign). Responses were measured using either electromyography (EMG) on proximal arm muscles during reaching-like movements, or with finger key-presses. In both EMG and button press measurements, participants responded faster when the object was either presented in the same location as the responding hand, or was affordable, resulting in significant and independent spatial compatibility and affordance effects, but no interaction. Furthermore, while the spatial compatibility effect was present from the earliest stages of movement preparation and throughout the different stages of movement execution, the affordance effect was restricted to the early stages of movement execution. Finally, we tested a small group of unilateral arm amputees using EMG, and found residual spatial compatibility but no affordance, suggesting that spatial compatibility effects do not necessarily rely on individuals’ available affordances. Our results show dissociation between affordance and spatial compatibility effects, and suggest that rather than evoking the specific motor action most suitable for interaction with the viewed object, graspable objects prompt the motor system in a general, body-part independent fashion

High-latitude climate sensitivity to ice-sheet forcing over the last 120kyr

Interpretation of ice-core records is currently limited by paucity of modelling at adequate temporal and spatial resolutions. Several key questions relate to mechanisms of polar amplification and inter-hemispheric coupling on glacial/interglacial timescales. Here, we present the first results from a large set of global ocean–atmosphere climate model ‘snap-shot’ simulations covering the last 120 000 years using the Hadley Centre climate model (HadCM3) at up to 1 kyr temporal resolution. Two sets of simulations were performed in order to examine the roles of orbit and greenhouse gases versus ice-sheet forcing of orbital-scale climate change. A series of idealised Heinrich events were also simulated, but no changes to aerosols or vegetation were prescribed. This paper focuses on high latitudes and inter-hemispheric linkages. The simulations reproduce polar temperature trends well compared to ice-core reconstructions, although the magnitude is underestimated. Polar amplification varies with obliquity, but this variability is dampened by including variations in land ice coverage, while the overall amplification factor increases. The relatively constant amplification of Antarctic temperatures (with ice-sheet forcing included) suggests it is possible to use Antarctic temperature reconstructions to estimate global changes (which are roughly half the magnitude). Atlantic Ocean overturning circulation varies considerably only with the introduction of Northern Hemisphere ice sheets, but only weakens in the North Atlantic in the deep glacial, when ocean–sea-ice feedbacks result in the movement of the region of deep convection to lower latitudes and with the introduction of freshwater to the surface North Atlantic in order to simulate Heinrich events.

Late Quaternary fire regimes of Australasia

We have compiled 223 sedimentary charcoal records from Australasia in order to examine the temporal and spatial variability of fire regimes during the Late Quaternary. While some of these records cover more than a full glacial cycle, here we focus on the last 70,000 years when the number of individual records in the compilation allows more robust conclusions. On orbital time scales, fire in Australasia predominantly reflects climate, with colder periods characterized by less and warmer intervals by more biomass burning. The composite record for the region also shows considerable millennial-scale variability during the last glacial interval (73.5–14.7 ka). Within the limits of the dating uncertainties of individual records, the variability shown by the composite charcoal record is more similar to the form, number and timing of Dansgaard–Oeschger cycles as observed in Greenland ice cores than to the variability expressed in the Antarctic ice-core record. The composite charcoal record suggests increased biomass burning in the Australasian region during Greenland Interstadials and reduced burning during Greenland Stadials. Millennial-scale variability is characteristic of the composite record of the sub-tropical high pressure belt during the past 21 ka, but the tropics show a somewhat simpler pattern of variability with major peaks in biomass burning around 15 ka and 8 ka. There is no distinct change in fire regime corresponding to the arrival of humans in Australia at 50 ± 10 ka and no correlation between archaeological evidence of increased human activity during the past 40 ka and the history of biomass burning. However, changes in biomass burning in the last 200 years may have been exacerbated or influenced by humans.

The residual circulation of the Southern Ocean: which spatio-temporal scales are needed?

The Southern Ocean circulation consists of a complicated mixture of processes and phenomena that arise at different time and spatial scales which need to be parametrized in the state-of-the-art climate models. The temporal and spatial scales that give rise to the present-day residual mean circulation are here investigated by calculating the Meridional Overturning Circulation (MOC) in density coordinates from an eddy-permitting global model. The region sensitive to the temporal decomposition is located between 38°S and 63°S, associated with the eddy-induced transport. The ‘‘Bolus’’ component of the residual circulation corresponds to the eddy-induced transport. It is dominated by timescales between 1 month and 1 year. The temporal behavior of the transient eddies is examined in splitting the ‘‘Bolus’’ component into a ‘‘Seasonal’’, an ‘‘Eddy’’ and an ‘‘Inter-monthly’’ component, respectively representing the correlation between density and velocity fluctuations due to the average seasonal cycle, due to mesoscale eddies and due to large-scale motion on timescales longer than one month that is not due to the seasonal cycle. The ‘‘Seasonal’’ bolus cell is important at all latitudes near the surface. The ‘‘Eddy’’ bolus cell is dominant in the thermocline between 50°S and 35°S and over the whole ocean depth at the latitude of the Drake Passage. The ‘‘Inter-monthly’’ bolus cell is important in all density classes and is maximal in the Brazil–Malvinas Confluence and the Agulhas Return Current. The spatial decomposition indicates that a large part of the Eulerian mean circulation is recovered for spatial scales larger than 11.25°, implying that small-scale meanders in the Antarctic Circumpolar Current (ACC), near the Subantarctic and Polar Fronts, and near the Subtropical Front are important in the compensation of the Eulerian mean flow.

Present state of global wetland extent and wetland methane modelling: methodology of a model inter-comparison project (WETCHIMP)

The Wetland and Wetland CH4 Intercomparison of Models Project (WETCHIMP) was created to evaluate our present ability to simulate large-scale wetland characteristics and corresponding methane (CH4) emissions. A multi-model comparison is essential to evaluate the key uncertainties in the mechanisms and parameters leading to methane emissions. Ten modelling groups joined WETCHIMP to run eight global and two regional models with a common experimental protocol using the same climate and atmospheric carbon dioxide (CO2) forcing datasets. We reported the main conclusions from the intercomparison effort in a companion paper (Melton et al., 2013). Here we provide technical details for the six experiments, which included an equilibrium, a transient, and an optimized run plus three sensitivity experiments (temperature, precipitation, and atmospheric CO2 concentration). The diversity of approaches used by the models is summarized through a series of conceptual figures, and is used to evaluate the wide range of wetland extent and CH4 fluxes predicted by the models in the equilibrium run. We discuss relationships among the various approaches and patterns in consistencies of these model predictions. Within this group of models, there are three broad classes of methods used to estimate wetland extent: prescribed based on wetland distribution maps, prognostic relationships between hydrological states based on satellite observations, and explicit hydrological mass balances. A larger variety of approaches was used to estimate the net CH4 fluxes from wetland systems. Even though modelling of wetland extent and CH4 emissions has progressed significantly over recent decades, large uncertainties still exist when estimating CH4 emissions: there is little consensus on model structure or complexity due to knowledge gaps, different aims of the models, and the range of temporal and spatial resolutions of the models.

Combined use of satellite estimates and rain gauge observations to generate high-quality historical rainfall time series over Ethiopia

Climate data are used in a number of applications including climate risk management and adaptation to climate change. However, the availability of climate data, particularly throughout rural Africa, is very limited. Available weather stations are unevenly distributed and mainly located along main roads in cities and towns. This imposes severe limitations to the availability of climate information and services for the rural community where, arguably, these services are needed most. Weather station data also suffer from gaps in the time series. Satellite proxies, particularly satellite rainfall estimate, have been used as alternatives because of their availability even over remote parts of the world. However, satellite rainfall estimates also suffer from a number of critical shortcomings that include heterogeneous time series, short time period of observation, and poor accuracy particularly at higher temporal and spatial resolutions. An attempt is made here to alleviate these problems by combining station measurements with the complete spatial coverage of satellite rainfall estimates. Rain gauge observations are merged with a locally calibrated version of the TAMSAT satellite rainfall estimates to produce over 30-years (1983-todate) of rainfall estimates over Ethiopia at a spatial resolution of 10 km and a ten-daily time scale. This involves quality control of rain gauge data, generating locally calibrated version of the TAMSAT rainfall estimates, and combining these with rain gauge observations from national station network. The infrared-only satellite rainfall estimates produced using a relatively simple TAMSAT algorithm performed as good as or even better than other satellite rainfall products that use passive microwave inputs and more sophisticated algorithms. There is no substantial difference between the gridded-gauge and combined gauge-satellite products over the test area in Ethiopia having a dense station network; however, the combined product exhibits better quality over parts of the country where stations are sparsely distributed.