An investigation of QBO signals in the east Asian and Indian monsoon in GCM experiments

Observations have shown that the monsoon is a highly variable phenomenon of the tropical troposphere, which exhibits significant variance in the temporal range of two to three years. The reason for this specific interannual variability has not yet been identified unequivocally. Observational analyses have also shown that EI Niño indices or western Pacific SSTs exhibit some power in the two to three year period range and therefore it was suggested that an ocean-atmosphere interaction could excite and support such a cycle. Similar mechanisms include land-surface-atmosphere interaction as a possible driving mechanism. A rather different explanation could be provided by a forcing mechanism based on the quasi-biennial oscillation of the zonal wind in the lower equatorial stratosphere (QBO). The QBO is a phenomenon driven by equatorial waves with periods of some days which are excited in the troposphere. Provided that the monsoon circulation reacts to the modulation of tropopause conditions as forced by the QBO, this could explain monsoon variability in the quasi-biennial window. The possibility of a QBO-driven monsoon variability is investigated in this study in a number of general circulation model experiments where the QBO is assimilated to externally controlled phase states. These experiments show that the boreal summer monsoon is significantly influenced by the QBO. A QBO westerly phase implies less precipitation in the western Pacific, but more in India, in agreement with observations. The austral summer monsoon is exposed to similar but weaker mechanisms and the precipitation does not change significantly.

Simulating the interannual variability of major dust storms on Mars using variable lifting thresholds

The redistribution of a finite amount of martian surface dust during global dust storms and in the intervening periods has been modelled in a dust lifting version of the UK Mars General Circulation Model. When using a constant, uniform threshold in the model’s wind stress lifting parameterisation and assuming an unlimited supply of surface dust, multiannual simulations displayed some variability in dust lifting activity from year to year, arising from internal variability manifested in surface wind stress, but dust storms were limited in size and formed within a relatively short seasonal window. Lifting thresholds were then allowed to vary at each model gridpoint, dependent on the rates of emission or deposition of dust. This enhanced interannual variability in dust storm magnitude and timing, such that model storms covered most of the observed ranges in size and initiation date within a single multiannual simulation. Peak storm magnitude in a given year was primarily determined by the availability of surface dust at a number of key sites in the southern hemisphere. The observed global dust storm (GDS) frequency of roughly one in every 3 years was approximately reproduced, but the model failed to generate these GDSs spontaneously in the southern hemisphere, where they have typically been observed to initiate. After several years of simulation, the surface threshold field—a proxy for net change in surface dust density—showed good qualitative agreement with the observed pattern of martian surface dust cover. The model produced a net northward cross-equatorial dust mass flux, which necessitated the addition of an artificial threshold decrease rate in order to allow the continued generation of dust storms over the course of a multiannual simulation. At standard model resolution, for the southward mass flux due to cross-equatorial flushing storms to offset the northward flux due to GDSs on a timescale of ∼3 years would require an increase in the former by a factor of 3–4. Results at higher model resolution and uncertainties in dust vertical profiles mean that quasi-periodic redistribution of dust on such a timescale nevertheless appears to be a plausible explanation for the observed GDS frequency.

Sparse probability density function estimation using the minimum integrated square error

We develop a new sparse kernel density estimator using a forward constrained regression framework, within which the nonnegative and summing-to-unity constraints of the mixing weights can easily be satisfied. Our main contribution is to derive a recursive algorithm to select significant kernels one at time based on the minimum integrated square error (MISE) criterion for both the selection of kernels and the estimation of mixing weights. The proposed approach is simple to implement and the associated computational cost is very low. Specifically, the complexity of our algorithm is in the order of the number of training data N, which is much lower than the order of N2 offered by the best existing sparse kernel density estimators. Numerical examples are employed to demonstrate that the proposed approach is effective in constructing sparse kernel density estimators with comparable accuracy to those of the classical Parzen window estimate and other existing sparse kernel density estimators.

Near-infrared water vapour self-continuum at close to room temperature

The gaseous absorption of solar radiation within near-infrared atmospheric windows in the Earth's atmosphere is dominated by the water vapour continuum. Recent measurements by Baranov et al. (2011) [17] in 2500 cm−1 (4 μm) window and by Ptashnik et al. (2011) [18] in a few near-infrared windows revealed that the self-continuum absorption is typically an order of magnitude stronger than given by the MT_CKD continuum model prior to version 2.5. Most of these measurements, however, were made at elevated temperatures, which makes their application to atmospheric conditions difficult. Here we report new laboratory measurements of the self-continuum absorption at 289 and 318 K in the near-infrared spectral region 1300–8000 cm−1, using a multipass 30 m base cell with total optical path 612 m. Our results confirm the main conclusions of the previous measurements both within bands and in windows. Of particular note is that we present what we believe to be the first near-room temperature measurement using Fourier Transform Spectrometry of the self-continuum in the 6200 cm−1 (1.6 μm) window, which provides tentative evidence that, at such temperatures, the water vapour continuum absorption may be as strong as it is in 2.1 μm and 4 μm windows and up to 2 orders of magnitude stronger than the MT_CKD-2.5 continuum. We note that alternative methods of measuring the continuum in this window have yielded widely differing assessment of its strength, which emphasises the need for further measurements.

When long-range zero-lag synchronization is feasible in cortical networks

Many studies have reported long-range synchronization of neuronal activity between brain areas, in particular in the beta and gamma bands with frequencies in the range of 14–30 and 40–80 Hz, respectively. Several studies have reported synchrony with zero phase lag, which is remarkable considering the synaptic and conduction delays inherent in the connections between distant brain areas. This result has led to many speculations about the possible functional role of zero-lag synchrony, such as for neuronal communication, attention, memory, and feature binding. However, recent studies using recordings of single-unit activity and local field potentials report that neuronal synchronization may occur with non-zero phase lags. This raises the questions whether zero-lag synchrony can occur in the brain and, if so, under which conditions. We used analytical methods and computer simulations to investigate which connectivity between neuronal populations allows or prohibits zero-lag synchrony. We did so for a model where two oscillators interact via a relay oscillator. Analytical results and computer simulations were obtained for both type I Mirollo–Strogatz neurons and type II Hodgkin–Huxley neurons. We have investigated the dynamics of the model for various types of synaptic coupling and importantly considered the potential impact of Spike-Timing Dependent Plasticity (STDP) and its learning window. We confirm previous results that zero-lag synchrony can be achieved in this configuration. This is much easier to achieve with Hodgkin–Huxley neurons, which have a biphasic phase response curve, than for type I neurons. STDP facilitates zero-lag synchrony as it adjusts the synaptic strengths such that zero-lag synchrony is feasible for a much larger range of parameters than without STDP.

Still feeling it: the time course of emotional recovery from an attentional perspective

Emotional reactivity and the time taken to recover, particularly from negative, stressful, events, are inextricably linked, and both are crucial for maintaining well-being. It is unclear, however, to what extent emotional reactivity during stimulus onset predicts the time course of recovery after stimulus offset. To address this question, 25 participants viewed arousing (negative and positive) and neutral pictures from the International Affective Picture System (IAPS) followed by task-relevant face targets, which were to be gender categorized. Faces were presented early (400–1500 ms) or late (2400–3500 ms) after picture offset to capture the time course of recovery from emotional stimuli. Measures of reaction time (RT), as well as face-locked N170 and P3 components were taken as indicators of the impact of lingering emotion on attentional facilitation or interference. Electrophysiological effects revealed negative and positive images to facilitate face-target processing on the P3 component, regardless of temporal interval. At the individual level, increased reactivity to: (1) negative pictures, quantified as the IAPS picture-locked Late Positive Potential (LPP), predicted larger attentional interference on the face-locked P3 component to faces presented in the late time window after picture offset. (2) Positive pictures, denoted by the LPP, predicted larger facilitation on the face-locked P3 component to faces presented in the earlier time window after picture offset. These results suggest that subsequent processing is still impacted up to 3500 ms after the offset of negative pictures and 1500 ms after the offset of positive pictures for individuals reacting more strongly to these pictures, respectively. Such findings emphasize the importance of individual differences in reactivity when predicting the temporality of emotional recovery. The current experimental model provides a novel basis for future research aiming to identify profiles of adaptive and maladaptive recovery.

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.

Estimation of sea surface temperature from the Spinning Enhanced Visible and Infra Red Imager, improved using numerical weather prediction

Most of the operational Sea Surface Temperature (SST) products derived from satellite infrared radiometry use multi-spectral algorithms. They show, in general, reasonable performances with root mean square (RMS) residuals around 0.5 K when validated against buoy measurements, but have limitations, particularly a component of the retrieval error that relates to such algorithms' limited ability to cope with the full variability of atmospheric absorption and emission. We propose to use forecast atmospheric profiles and a radiative transfer model to simulate the algorithmic errors of multi-spectral algorithms. In the practical case of SST derived from the Spinning Enhanced Visible and Infrared Imager (SEVIRI) onboard Meteosat Second Generation (MSG), we demonstrate that simulated algorithmic errors do explain a significant component of the actual errors observed for the non linear (NL) split window algorithm in operational use at the Centre de Météorologie Spatiale (CMS). The simulated errors, used as correction terms, reduce significantly the regional biases of the NL algorithm as well as the standard deviation of the differences with drifting buoy measurements. The availability of atmospheric profiles associated with observed satellite-buoy differences allows us to analyze the origins of the main algorithmic errors observed in the SEVIRI field of view: a negative bias in the inter-tropical zone, and a mid-latitude positive bias. We demonstrate how these errors are explained by the sensitivity of observed brightness temperatures to the vertical distribution of water vapour, propagated through the SST retrieval algorithm.

Sea surface temperature estimation from the Geostationary Operational Environmental Satellite 12 (GOES-12)

This paper describes the techniques used to obtain sea surface temperature (SST) retrievals from the Geostationary Operational Environmental Satellite 12 (GOES-12) at the National Oceanic and Atmospheric Administration’s Office of Satellite Data Processing and Distribution. Previous SST retrieval techniques relying on channels at 11 and 12 μm are not applicable because GOES-12 lacks the latter channel. Cloud detection is performed using a Bayesian method exploiting fast-forward modeling of prior clear-sky radiances using numerical weather predictions. The basic retrieval algorithm used at nighttime is based on a linear combination of brightness temperatures at 3.9 and 11 μm. In comparison with traditional split window SSTs (using 11- and 12-μm channels), simulations show that this combination has maximum scatter when observing drier colder scenes, with a comparable overall performance. For daytime retrieval, the same algorithm is applied after estimating and removing the contribution to brightness temperature in the 3.9-μm channel from solar irradiance. The correction is based on radiative transfer simulations and comprises a parameterization for atmospheric scattering and a calculation of ocean surface reflected radiance. Potential use of the 13-μm channel for SST is shown in a simulation study: in conjunction with the 3.9-μm channel, it can reduce the retrieval error by 30%. Some validation results are shown while a companion paper by Maturi et al. shows a detailed analysis of the validation results for the operational algorithms described in this present article.

Grouping of 24 apple cultivars on the basis of starch degradation rate and their fruit pattern

The ripening processes of 24 apple cultivars were examined in the United Kingdom National Fruit Collection in 2010. Basically the starch content, and additionally ground colour, water-soluble solids content and flesh firmness were studied during ripening. The degradation of the starch content was evaluated using a 0–10 scale. A starch degradation value of 50% was taken to be the optimum harvest date, with harvest beginning at a value of 40% and finishing at 60%. Depending on the cultivar, this represented a harvest window of 9 to 21 days. Later ripening cultivars matured more slowly, leading to a longer harvesting period, with the exception of cv. Feuillemorte. Pronounced differences were observed among the cultivars on the basis of the starch degradation pattern, allowing them to be divided into four groups. Separate charts were elaborated for each group that are recommended for use in practice.

Radiative forcing and global warming potentials of 11 halogenated compounds

Radiative forcing values have been calculated for 11 halogenated compounds which are in current use or which have been suggested as possible replacements for the chlorofluorocarbons. Absorption cross-sections measured over a range of atmospheric temperature and pressure conditions as part of a multi-laboratory programme have been used together with a narrow band radiative transfer model. We provide a “best estimate” radiative forcing taking into account the likely vertical profile of the gas in each case. The Global Warming Potential over a variety of time horizons has also been calculated where the lifetime is available. We present the first such information for 1,2-dichloroethane. For chloroform our radiative forcing is 5 times higher than the value used in previous assessments, possibly because these ignored the effect of absorption outside the 800–1200 cm−1 “window”. For several of the other compounds considered here, our forcing is between 10 and 30% lower than previous assessments. The perfluorocarbons have been found to have large global warming potentials, many times that of CFC-11, due to both strong absorption and long lifetimes. The importance of absorption features at wavenumbers below 800 cm−1 and the effect of temperature variations in absorption cross-section on the radiative forcing are also investigated.

Optimizing performance of the lightshelf by modifying ceiling geometry in highly luminous climates

The impact of ceiling geometries on the performance of lightshelves was investigated using physical model experiments and radiance simulations. Illuminance level and distribution uniformity were assessed for a working plane in a large space located in sub-tropical climate regions where innovative systems for daylighting and shading are required. It was found that the performance of the lightshelf can be improved by changing the ceiling geometry; the illuminance level increased in the rear of the room and decreased in the front near the window compared to rooms having conventional horizontal ceilings. Moreover, greater uniformity was achieved throughout the room as a result of reducing the difference in the illuminance level between the front and rear of the room. Radiance simulation results were found to be in good agreement with physical model data obtained under a clear sky and high solar radiation. The best ceiling shape was found to be one that is curved in the front and rear of the room.

PDFOS: PDF estimation based over-sampling for imbalanced two-class problems

This contribution proposes a novel probability density function (PDF) estimation based over-sampling (PDFOS) approach for two-class imbalanced classification problems. The classical Parzen-window kernel function is adopted to estimate the PDF of the positive class. Then according to the estimated PDF, synthetic instances are generated as the additional training data. The essential concept is to re-balance the class distribution of the original imbalanced data set under the principle that synthetic data sample follows the same statistical properties. Based on the over-sampled training data, the radial basis function (RBF) classifier is constructed by applying the orthogonal forward selection procedure, in which the classifier’s structure and the parameters of RBF kernels are determined using a particle swarm optimisation algorithm based on the criterion of minimising the leave-one-out misclassification rate. The effectiveness of the proposed PDFOS approach is demonstrated by the empirical study on several imbalanced data sets.

A method to assess individual research outputs

We present here a straightforward method which can be used to obtain a quantitative indication of an individual research output for an academic. Different versions, selections and options are presented to enable a user to easily calculate values both for stand-alone papers and overall for the collection of outputs for a person. The procedure is particularly useful as a metric to give a quantitative indication of the research output of a person over a time window. Examples are included to show how the method works in practice and how it compares to alternative techniques.