Ensemble clustering in deterministic ensemble Kalman filters

Ensemble clustering (EC) can arise in data assimilation with ensemble square root filters (EnSRFs) using non-linear models: an M-member ensemble splits into a single outlier and a cluster of M−1 members. The stochastic Ensemble Kalman Filter does not present this problem. Modifications to the EnSRFs by a periodic resampling of the ensemble through random rotations have been proposed to address it. We introduce a metric to quantify the presence of EC and present evidence to dispel the notion that EC leads to filter failure. Starting from a univariate model, we show that EC is not a permanent but transient phenomenon; it occurs intermittently in non-linear models. We perform a series of data assimilation experiments using a standard EnSRF and a modified EnSRF by a resampling though random rotations. The modified EnSRF thus alleviates issues associated with EC at the cost of traceability of individual ensemble trajectories and cannot use some of algorithms that enhance performance of standard EnSRF. In the non-linear regimes of low-dimensional models, the analysis root mean square error of the standard EnSRF slowly grows with ensemble size if the size is larger than the dimension of the model state. However, we do not observe this problem in a more complex model that uses an ensemble size much smaller than the dimension of the model state, along with inflation and localisation. Overall, we find that transient EC does not handicap the performance of the standard EnSRF.

EEG activation differences in the pre-motor cortex and supplementary motor area between normal individuals with high and low traits of autism

The human mirror neuron system (hMNS) is believed to provide a basic mechanism for social cognition. Event-related desynchronization (ERD) in alpha (8–12 Hz) and low beta band (12–20 Hz) over sensori-motor cortex has been suggested to index mirror neurons' activity. We tested whether autistic traits revealed by high and low scores on the Autistic Quotient (AQ) in the normal population are linked to variations in the electroencephalogram (EEG) over motor, pre-motor cortex and supplementary motor area (SMA) during action observation. Results revealed that in the low AQ group, the pre-motor cortex and SMA were more active during hand action than static hand observation whereas in the high AQ group the same areas were active both during static and hand action observation. In fact participants with high traits of autism showed greater low beta ERD while observing the static hand than those with low traits and this low beta ERD was not significantly different when they watched hand actions. Over primary motor cortex, the classical alpha and low beta ERD during hand actions relative to static hand observation was found across all participants. These findings suggest that the observation–execution matching system works differently according to the degree of autism traits in the normal population and that this is differentiated in terms of the EEG according to scalp site and bandwidth.

Chronic systemic therapy with low-dose morpholino oligomers ameliorates the pathology and normalizes locomotor behavior in mdx Mice

The administration of antisense oligonucleotides (AOs) to skip one or more exons in mutated forms of the DMD gene and so restore the reading frame of the transcript is one of the most promising approaches to treat Duchenne muscular dystrophy (DMD). At present, preclinical studies demonstrating the efficacy and safety of long-term AO administration have not been conducted. Furthermore, it is essential to determine the minimal effective dose and frequency of administration. In this study, two different low doses (LDs) of phosphorodiamidate morpholino oligomer (PMO) designed to skip the mutated exon 23 in the mdx dystrophic mouse were administered for up to 12 months. Mice treated for 50 weeks showed a substantial dose-related amelioration of the pathology, particularly in the diaphragm. Moreover, the generalized physical activity was profoundly enhanced compared to untreated mdx mice showing that widespread, albeit partial, dystrophin expression restores the normal activity in mdx mice. Our results show for the first time that a chronic long-term administration of LDs of unmodified PMO, equivalent to doses in use in DMD boys, is safe, significantly ameliorates the muscular dystrophic phenotype and improves the activity of dystrophin-deficient mice, thus encouraging the further clinical translation of this approach in humans.

Modelling climate impact on floods with ensemble climate projections

The evidence provided by modelled assessments of future climate impact on flooding is fundamental to water resources and flood risk decision making. Impact models usually rely on climate projections from global and regional climate models (GCM/RCMs). However, challenges in representing precipitation events at catchment-scale resolution mean that decisions must be made on how to appropriately pre-process the meteorological variables from GCM/RCMs. Here the impacts on projected high flows of differing ensemble approaches and application of Model Output Statistics to RCM precipitation are evaluated while assessing climate change impact on flood hazard in the Upper Severn catchment in the UK. Various ensemble projections are used together with the HBV hydrological model with direct forcing and also compared to a response surface technique. We consider an ensemble of single-model RCM projections from the current UK Climate Projections (UKCP09); multi-model ensemble RCM projections from the European Union's FP6 ‘ENSEMBLES’ project; and a joint probability distribution of precipitation and temperature from a GCM-based perturbed physics ensemble. The ensemble distribution of results show that flood hazard in the Upper Severn is likely to increase compared to present conditions, but the study highlights the differences between the results from different ensemble methods and the strong assumptions made in using Model Output Statistics to produce the estimates of future river discharge. The results underline the challenges in using the current generation of RCMs for local climate impact studies on flooding. Copyright © 2012 Royal Meteorological Society

Microbiological, chemical and rheological properties of low fat set yoghurt produced with exopolysaccharide (EPS) producing Bifidobacterium strains.

In this work, the microbiological and physicochemical differences of three types of low fat set yoghurts were studied, as well as the changes taking place during storage at 4 °C for 28 days. The first yoghurt was produced with yoghurt starters and exopolysaccharide (EPS) producing Bifidobacterium longum subsp. infantis CCUG 52486 (CCUGY), the second with yoghurt starters and Bifidobacterium infantis NCIMB 702205 (NCIMBY) and the third with just yoghurt starters (control yoghurt). No significant differences were observed in terms of cell concentrations; for all three yoghurts, similar final cell concentrations were obtained for the yoghurt starter cultures (~7.5 log cfu g−1) and the Bifidobacterium strains (~7.8 log cfu g−1). Both Bifidobacterium survived well during storage, as in both cases the cell viability decreased by less than 0.5 log cfu g−1after 28 days of storage. A decrease in pH followed by an increase in lactic acid was observed during storage for all three yoghurts, which was mostly attributed to the activity of the yoghurt starter cultures. The two yoghurts with the EPS producing Bifidobacterium strains exhibited lower syneresis than the control yoghurt. The lowest was shown by CCUGY, which also exhibited the highest storage modulus and firmness, and a well defined porous web-like structure in cryo-SEM. Examination of the micro-structure of the yoghurts using cryo-scanning electron microscopy (cryo-SEM) indicated that the above observations were due to the interaction between the EPS and the milk proteins. Overall, the results indicated that the EPS producing Bifidobacterium longum subsp. infantis CCUG 52486 is the most promising strain, and can be used with yoghurt starter cultures to manufacture low fat set yoghurt with probiotic activities and at the same time enhanced physicochemical and rheological properties.

Climate predictability experiments with a general circulation model

The atmospheric response to the evolution of the global sea surface temperatures from 1979 to 1992 is studied using the Max-Planck-Institut 19 level atmospheric general circulation model, ECHAM3 at T 42 resolution. Five separate 14-year integrations are performed and results are presented for each individual realization and for the ensemble-averaged response. The results are compared to a 30-year control integration using a climate monthly mean state of the sea surface temperatures and to analysis data. It is found that the ECHAM3 model, by and large, does reproduce the observed response pattern to El Nin˜o and La Nin˜a. During the El Nin˜ o events, the subtropical jet streams in both hemispheres are intensified and displaced equatorward, and there is a tendency towards weak upper easterlies over the equator. The Southern Oscillation is a very stable feature of the integrations and is accurately reproduced in all experiments. The inter-annual variability at middle- and high-latitudes, on the other hand, is strongly dominated by chaotic dynamics, and the tropical SST forcing only modulates the atmospheric circulation. The potential predictability of the model is investigated for six different regions. Signal to noise ratio is large in most parts of the tropical belt, of medium strength in the western hemisphere and generally small over the European area. The ENSO signal is most pronounced during the boreal spring. A particularly strong signal in the precipitation field in the extratropics during spring can be found over the southern United States. Western Canada is normally warmer during the warm ENSO phase, while northern Europe is warmer than normal during the ENSO cold phase. The reason is advection of warm air due to a more intense Pacific low than normal during the warm ENSO phase and a more intense Icelandic low than normal during the cold ENSO phase, respectively.

A low-order model investigation of the analysis of gravity waves in the ensemble Kalman filter.

The behavior of the ensemble Kalman filter (EnKF) is examined in the context of a model that exhibits a nonlinear chaotic (slow) vortical mode coupled to a linear (fast) gravity wave of a given amplitude and frequency. It is shown that accurate recovery of both modes is enhanced when covariances between fast and slow normal-mode variables (which reflect the slaving relations inherent in balanced dynamics) are modeled correctly. More ensemble members are needed to recover the fast, linear gravity wave than the slow, vortical motion. Although the EnKF tends to diverge in the analysis of the gravity wave, the filter divergence is stable and does not lead to a great loss of accuracy. Consequently, provided the ensemble is large enough and observations are made that reflect both time scales, the EnKF is able to recover both time scales more accurately than optimal interpolation (OI), which uses a static error covariance matrix. For OI it is also found to be problematic to observe the state at a frequency that is a subharmonic of the gravity wave frequency, a problem that is in part overcome by the EnKF.However, error in themodeled gravity wave parameters can be detrimental to the performance of the EnKF and remove its implied advantages, suggesting that a modified algorithm or a method for accounting for model error is needed.

Ensemble projections for wine production in the Douro Valley of Portugal

Wine production is largely governed by atmospheric conditions, such as air temperature and precipitation, together with soil management and viticultural/enological practices. Therefore, anthropogenic climate change is likely to have important impacts on the winemaking sector worldwide. An important winemaking region is the Portuguese Douro Valley, which is known by its world-famous Port Wine. The identification of robust relationships between atmospheric factors and wine parameters is of great relevance for the region. A multivariate linear regression analysis of a long wine production series (1932–2010) reveals that high rainfall and cool temperatures during budburst, shoot and inflorescence development (February-March) and warm temperatures during flowering and berry development (May) are generally favourable to high production. The probabilities of occurrence of three production categories (low, normal and high) are also modelled using multinomial logistic regression. Results show that both statistical models are valuable tools for predicting the production in a given year with a lead time of 3–4 months prior to harvest. These statistical models are applied to an ensemble of 16 regional climate model experiments following the SRES A1B scenario to estimate possible future changes. Wine production is projected to increase by about 10 % by the end of the 21st century, while the occurrence of high production years is expected to increase from 25 % to over 60 %. Nevertheless, further model development will be needed to include other aspects that may shape production in the future. In particular, the rising heat stress and/or changes in ripening conditions could limit the projected production increase in future decades.

Changes in storm track and cyclone activity in three SRES ensemble experiments with the ECHAM5/MPI-OM1 GCM

Synoptic activity over the Northern Hemisphere is evaluated in ensembles of ECHAM5/MPI-OM1 simulations for recent climate conditions (20C) and for three climate scenarios (following SRES A1B, A2, B1). A close agreement is found between the simulations for present day climate and the respective results from reanalysis. Significant changes in the winter mid-tropospheric storm tracks are detected in all three scenario simulations. Ensemble mean climate signals are rather similar, with particularly large activity increases downstream of the Atlantic storm track over Western Europe. The magnitude of this signal is largely dependent on the imposed change in forcing. However, differences between individual ensemble members may be large. With respect to the surface cyclones, the scenario runs produce a reduction in cyclonic track density over the mid-latitudes, even in the areas with increasing mid-tropospheric activity. The largest decrease in track densities occurs at subtropical latitudes, e.g., over the Mediterranean Basin. An increase of cyclone intensities is detected for limited areas (e.g., near Great Britain and Aleutian Isles) for the A1B and A2 experiments. The changes in synoptic activity are associated with alterations of the Northern Hemisphere circulation and background conditions (blocking frequencies, jet stream). The North Atlantic Oscillation index also shows increased values with enhanced forcing. With respect to the effects of changing synoptic activity, the regional change in cyclone intensities is accompanied by alterations of the extreme surface winds, with increasing values over Great Britain, North and Baltic Seas, as well as the areas with vanishing sea ice, and decreases over much of the subtropics.

Ocean heat transport and water vapor greenhouse in a warm equable climate: a new look at the low gradient paradox

The authors study the role of ocean heat transport (OHT) in the maintenance of a warm, equable, ice-free climate. An ensemble of idealized aquaplanet GCM calculations is used to assess the equilibrium sensitivity of global mean surface temperature and its equator-to-pole gradient (ΔT) to variations in OHT, prescribed through a simple analytical formula representing export out of the tropics and poleward convergence. Low-latitude OHT warms the mid- to high latitudes without cooling the tropics; increases by 1°C and ΔT decreases by 2.6°C for every 0.5-PW increase in OHT across 30° latitude. This warming is relatively insensitive to the detailed meridional structure of OHT. It occurs in spite of near-perfect atmospheric compensation of large imposed variations in OHT: the total poleward heat transport is nearly fixed. The warming results from a convective adjustment of the extratropical troposphere. Increased OHT drives a shift from large-scale to convective precipitation in the midlatitude storm tracks. Warming arises primarily from enhanced greenhouse trapping associated with convective moistening of the upper troposphere. Warming extends to the poles by atmospheric processes even in the absence of high-latitude OHT. A new conceptual model for equable climates is proposed, in which OHT plays a key role by driving enhanced deep convection in the midlatitude storm tracks. In this view, the climatic impact of OHT depends on its effects on the greenhouse properties of the atmosphere, rather than its ability to increase the total poleward energy transport.

Low cost, patterning of human hNT brain cells on parylene-C with UV & IR laser machining

This paper describes the use of 800nm femtosecond infrared (IR) and 248nm nanosecond ultraviolet (UV) laser radiation in performing ablative micromachining of parylene-C on SiO2 substrates for the patterning of human hNT astrocytes. Results are presented that support the validity of using IR laser ablative micromachining for patterning human hNT astrocytes cells while UV laser radiation produces photo-oxidation of the parylene-C and destroys cell patterning. The findings demonstrate how IR laser ablative micromachining of parylene-C on SiO2 substrates can offer a low cost, accessible alternative for rapid prototyping, high yield cell patterning.

On the lack of stratospheric dynamical variability in low-top versions of the CMIP5 models

We describe the main differences in simulations of stratospheric climate and variability by models within the fifth Coupled Model Intercomparison Project (CMIP5) that have a model top above the stratopause and relatively fine stratospheric vertical resolution (high-top), and those that have a model top below the stratopause (low-top). Although the simulation of mean stratospheric climate by the two model ensembles is similar, the low-top model ensemble has very weak stratospheric variability on daily and interannual time scales. The frequency of major sudden stratospheric warming events is strongly underestimated by the low-top models with less than half the frequency of events observed in the reanalysis data and high-top models. The lack of stratospheric variability in the low-top models affects their stratosphere-troposphere coupling, resulting in short-lived anomalies in the Northern Annular Mode, which do not produce long-lasting tropospheric impacts, as seen in observations. The lack of stratospheric variability, however, does not appear to have any impact on the ability of the low-top models to reproduce past stratospheric temperature trends. We find little improvement in the simulation of decadal variability for the high-top models compared to the low-top, which is likely related to the fact that neither ensemble produces a realistic dynamical response to volcanic eruptions.

Association of low adiponectin levels with the metabolic syndrome--the Chennai Urban Rural Epidemiology Study (CURES-4)

The aim of the study was to assess the relation of adiponectin levels with the metabolic syndrome in Asian Indians, a high-risk group for diabetes and premature coronary artery disease. The study was conducted on 100 (50 men and 50 women) type 2 diabetic subjects and 100 age and sex matched subjects with normal glucose tolerance selected from the Chennai Urban Rural Epidemiology Study, an ongoing population study in Chennai in southern India. Metabolic syndrome was defined using modified Adult Treatment Panel III (ATPIII) guidelines. Adiponectin values were significantly lower in diabetic subjects (men: 5.2 vs 8.3 microg/mL, P=.00l; women: 7.6 vs 11.1 microg/mL, P<.00l) and those with the metabolic syndrome (men: 5.0 vs 6.8 microg/mL, P=.01; women: 6.5 vs 9.9 microg/mL, P=.001) compared with those without. Linear regression analysis revealed adiponectin to be associated with body mass index (P<.05), waist circumference (P<.01), fasting plasma glucose (P=.001), glycated hemoglobin (P<.001), triglycerides (P<.00l), high-density lipoprotein (HDL) cholesterol (P<.001), cholesterol/HDL ratio (P<.00l), and insulin resistance measured by homeostasis assessment model (P<.00l). Factor analysis identified 2 factors: factor 1, negatively loaded with adiponectin and HDL cholesterol and positively loaded with triglycerides, waist circumference, and insulin resistance measured by homeostasis assessment model; and factor 2, with a positive loading of waist circumference and systolic and diastolic blood pressure. Logistic regression analysis revealed adiponectin to be negatively associated with metabolic syndrome (odds ratio [OR], 0.365; P<.001) even after adjusting for age (OR, 0.344; P<.00l), sex (OR, 0.293; P<.001), and body mass index (OR, 0.292; P<.00l). Lower adiponectin levels are associated with the metabolic syndrome per se and several of its components, particularly, diabetes, insulin resistance, and dyslipidemia in this urban south Indian population.

A generalised Bayesian instrumental variable approach under student t-distributed errors with application

Bayesian analysis is given of an instrumental variable model that allows for heteroscedasticity in both the structural equation and the instrument equation. Specifically, the approach for dealing with heteroscedastic errors in Geweke (1993) is extended to the Bayesian instrumental variable estimator outlined in Rossi et al. (2005). Heteroscedasticity is treated by modelling the variance for each error using a hierarchical prior that is Gamma distributed. The computation is carried out by using a Markov chain Monte Carlo sampling algorithm with an augmented draw for the heteroscedastic case. An example using real data illustrates the approach and shows that ignoring heteroscedasticity in the instrument equation when it exists may lead to biased estimates.

Coupling ensemble weather predictions based on TIGGE database with Grid-Xinanjiang model for flood forecast

The incorporation of numerical weather predictions (NWP) into a flood forecasting system can increase forecast lead times from a few hours to a few days. A single NWP forecast from a single forecast centre, however, is insufficient as it involves considerable non-predictable uncertainties and lead to a high number of false alarms. The availability of global ensemble numerical weather prediction systems through the THORPEX Interactive Grand Global Ensemble' (TIGGE) offers a new opportunity for flood forecast. The Grid-Xinanjiang distributed hydrological model, which is based on the Xinanjiang model theory and the topographical information of each grid cell extracted from the Digital Elevation Model (DEM), is coupled with ensemble weather predictions based on the TIGGE database (CMC, CMA, ECWMF, UKMO, NCEP) for flood forecast. This paper presents a case study using the coupled flood forecasting model on the Xixian catchment (a drainage area of 8826 km2) located in Henan province, China. A probabilistic discharge is provided as the end product of flood forecast. Results show that the association of the Grid-Xinanjiang model and the TIGGE database gives a promising tool for an early warning of flood events several days ahead.