An exploration of the equivalent weights particle filter

Particle filters are fully non-linear data assimilation techniques that aim to represent the probability distribution of the model state given the observations (the posterior) by a number of particles. In high-dimensional geophysical applications the number of particles required by the sequential importance resampling (SIR) particle filter in order to capture the high probability region of the posterior, is too large to make them usable. However particle filters can be formulated using proposal densities, which gives greater freedom in how particles are sampled and allows for a much smaller number of particles. Here a particle filter is presented which uses the proposal density to ensure that all particles end up in the high probability region of the posterior probability density function. This gives rise to the possibility of non-linear data assimilation in large dimensional systems. The particle filter formulation is compared to the optimal proposal density particle filter and the implicit particle filter, both of which also utilise a proposal density. We show that when observations are available every time step, both schemes will be degenerate when the number of independent observations is large, unlike the new scheme. The sensitivity of the new scheme to its parameter values is explored theoretically and demonstrated using the Lorenz (1963) model.

Quantitative ventilation assessments of idealized urban canopy layers with various urban layouts and the same building packing density

This paper investigates urban canopy layers (UCL) ventilation under neutral atmospheric condition with the same building area density (λp=0.25) and frontal area density (λf=0.25) but various urban sizes, building height variations, overall urban forms and wind directions. Turbulent airflows are first predicted by CFD simulations with standard k-ε model evaluated by wind tunnel data. Then air change rates per hour (ACH) and canopy purging flow rate (PFR) are numerically analyzed to quantify the rate of air exchange and the net ventilation capacity induced by mean flows and turbulence. With a parallel approaching wind (θ=0o), the velocity ratio first decreases in the adjustment region, followed by the fully-developed region where the flow reaches a balance. Although the flow quantities macroscopically keep constant, however ACH decreases and overall UCL ventilation becomes worse if urban size rises from 390m to 5km. Theoretically if urban size is infinite, ACH may reach a minimum value depending on local roof ventilation, and it rises from 1.7 to 7.5 if the standard deviation of building height variations increases (0% to 83.3%). Overall UCL ventilation capacity (PFR) with a square overall urban form (Lx=Ly=390m) is better as θ=0o than oblique winds (θ=15o, 30o, 45o), and it exceeds that of a staggered urban form under all wind directions (θ=0o to 45o), but is less than that of a rectangular urban form (Lx=570m, Ly=270m) under most wind directions (θ=30o to 90o). Further investigations are still required to quantify the net ventilation efficiency induced by mean flows and turbulence.

Drug permeability in 16HBE14o- airway cell layers correlates with absorption from the isolated perfused rat lung

The permeability of the lung is critical in determining the disposition of inhaled drugs and the respiratory epithelium provides the main physical barrier to drug absorption. The 16HBE14o- human bronchial epithelial cell line has been developed recently as a model of the airway epithelium. In this study, the transport of 10 low molecular weight compounds was measured in the 16HBE14o- cell layers, with apical to basolateral (absorptive) apparent permeability coefficients (P(app)) ranging from 0.4 x 10(-6)cms(-1) for Tyr-D-Arg-Phe-Phe-NH(2) to 25.2x10(-6)cms(-1) for metoprolol. Permeability in 16HBE14o- cells was found to correlate with previously reported P(app) in Caco-2 cells and absorption rates in the isolated perfused rat lung (k(a,lung)) and the rat lung in vivo (k(a,in vivo)). Log linear relationships were established between P(app) in 16HBE14o- cells and P(app) in Caco-2 cells (r(2)=0.82), k(a,lung) (r(2)=0.78) and k(a,in vivo) (r(2)=0.68). The findings suggest that permeability in 16HBE14o- cells may be useful to predict the permeability of compounds in the lung, although no advantage of using the organ-specific cell line 16HBE14o- compared to Caco-2 cells was found in this study.

Use of simulated intestinal fluids with Caco-2 cells and rat ileum

In most in vitro studies of oral drug permeability, little attempt is made to reproduce the gastrointestinal lumenal environment. The aim of this study was to evaluate the compatibility of simulated intestinal fluid (SIF) solutions with Caco-2 cell monolayers and Ussing chamber-mounted rat ileum under standard permeability experiment protocols. In preliminary experiments, fasted-state simulated intestinal fluid (FaSSIF) and fed-state simulated intestinal fluid (FeSSIF) solutions based on the dissolution medium formulae of Dressman and co-workers (1998) were modified for compatibility with Caco-2 cells to produce FaS-SIF and FeSSIF "transport" solutions for use with in vitro permeability models. For Caco-2 cells exposed to FaSSIF and FESSIF transport solutions, the transepithelial electrical resistance was maintained for over 4 h and mannitol permeability was equivalent to that in control (Hank's Balanced Salt Solution-treated) cell layers. Scanning electron microscopy revealed that microvilli generally maintained a normal distribution, although some shortening of microvilli and occasional small areas of denudation were observed. For rat ileum in the Ussing chambers, the potential difference (PD) collapsed to zero over 120 min when exposed to the FaSSIF transport solution and an even faster collapse of the PD was observed when the FeSSIF transport solution was used. Electron micrographs revealed erosion of the villi tips and substantial denudation of the microvilli after exposure of ileal tissue to FaSSIF and FeSSIF solutions, and permeability to mannitol was increased by almost two-fold. This study indicated that FaSSIF and FeSSIF transport solutions can be used with Caco-2 monolayers to evaluate drug permeability, but rat ileum in Ussing chambers is adversely affected by these solutions. Metoprolol permeability in Caco-2 experiments was reduced by 33% using the FaSSIF and 75% using the FeSSIF compared to permeability measured using HBSS. This illustrates that using physiological solutions can influence permeability measurements.

Serum protein layers on parylene-C and silicon oxide: effect on cell adhesion

Among the range of materials used in bioengineering, parylene-C has been used in combination with silicon oxide and in presence of the serum proteins, in cell patterning. However, the structural properties of adsorbed serum proteins on these substrates still remain elusive. In this study, we use an optical biosensing technique to decipher the properties of fibronectin (Fn) and serum albumin adsorbed on parylene-C and silicon oxide substrates. Our results show the formation of layers with distinct structural and adhesive properties. Thin, dense layers are formed on parylene-C, whereas thicker, more diffuse layers are formed on silicon oxide. These results suggest that Fn acquires a compact structure on parylene-C and a more extended structure on silicon oxide. Nonetheless, parylene-C and silicon oxide substrates coated with Fn host cell populations that exhibit focal adhesion complexes and good cell attachment. Albumin adopts a deformed structure on parylene-C and a globular structure on silicon oxide, and does not support significant cell attachment on either surface. Interestingly, the co-incubation of Fn and albumin at the ratio found in serum, results in the preferential adsorption of albumin on parylene-C and Fn on silicon oxide. This finding is supported by the exclusive formation of focal adhesion complexes in differentiated mouse embryonic stem cells (CGR8), cultured on Fn/albumin coated silicon oxide, but not on parylene-C. The detailed information provided in this study on the distinct properties of layers of serum proteins on substrates such as parylene-C and silicon oxide is highly significant in developing methods for cell patterning.

The equivalent-weights particle filter in a high-dimensional system

In general, particle filters need large numbers of model runs in order to avoid filter degeneracy in high-dimensional systems. The recently proposed, fully nonlinear equivalent-weights particle filter overcomes this requirement by replacing the standard model transition density with two different proposal transition densities. The first proposal density is used to relax all particles towards the high-probability regions of state space as defined by the observations. The crucial second proposal density is then used to ensure that the majority of particles have equivalent weights at observation time. Here, the performance of the scheme in a high, 65 500 dimensional, simplified ocean model is explored. The success of the equivalent-weights particle filter in matching the true model state is shown using the mean of just 32 particles in twin experiments. It is of particular significance that this remains true even as the number and spatial variability of the observations are changed. The results from rank histograms are less easy to interpret and can be influenced considerably by the parameter values used. This article also explores the sensitivity of the performance of the scheme to the chosen parameter values and the effect of using different model error parameters in the truth compared with the ensemble model runs.

The effect of the equivalent-weights particle filter on dynamical balance in a primitive equation model

The disadvantage of the majority of data assimilation schemes is the assumption that the conditional probability density function of the state of the system given the observations [posterior probability density function (PDF)] is distributed either locally or globally as a Gaussian. The advantage, however, is that through various different mechanisms they ensure initial conditions that are predominantly in linear balance and therefore spurious gravity wave generation is suppressed. The equivalent-weights particle filter is a data assimilation scheme that allows for a representation of a potentially multimodal posterior PDF. It does this via proposal densities that lead to extra terms being added to the model equations and means the advantage of the traditional data assimilation schemes, in generating predominantly balanced initial conditions, is no longer guaranteed. This paper looks in detail at the impact the equivalent-weights particle filter has on dynamical balance and gravity wave generation in a primitive equation model. The primary conclusions are that (i) provided the model error covariance matrix imposes geostrophic balance, then each additional term required by the equivalent-weights particle filter is also geostrophically balanced; (ii) the relaxation term required to ensure the particles are in the locality of the observations has little effect on gravity waves and actually induces a reduction in gravity wave energy if sufficiently large; and (iii) the equivalent-weights term, which leads to the particles having equivalent significance in the posterior PDF, produces a change in gravity wave energy comparable to the stochastic model error. Thus, the scheme does not produce significant spurious gravity wave energy and so has potential for application in real high-dimensional geophysical applications.

Working through the layers: curriculum implementation in language education

This chapter presents findings on English Language instruction at the lower primary level in the context of policies for curricular innovation at national, school and classroom levels. The focus is on policies which connect national and school levels, and on how they might be interpreted when implemented in multiple schools within Singapore’s educational system. Referring to case studies in two schools and to individual lesson observations in 10 schools, we found much agreement with national policies in terms of curriculum (i.e. lesson content and activity selection),leading to great uniformity in the lessons taught by different teachers in different schools. In addition, we found that schools had an important mediating influence on implementation of national policies. However, adoptions and adaptations of policy innovations at the classroom level were somewhat superficial as they were more related to changes in educational facilities and procedures than in philosophies.

New self-assembling multifunctional templates for the biofabrication and controlled self-release of cultured tissue

The need to source live human tissues for research and clinical applications has been a major driving force for the development of new biomaterials. Ideally, these should elicit the formation of scaffold-free tissues with native-like structure and composition. In this study, we describe a biologically interactive coating that combines the fabrication and subsequent self-release of live purposeful tissues using template–cell–environment feedback. This smart coating was formed from a self-assembling peptide amphiphile comprising a proteasecleavable sequence contiguous with a cell attachment and signaling motif. This multifunctional material was subsequently used not only to instruct human corneal or skin fibroblasts to adhere and deposit discreet multiple layers of native extracellular matrix but also to govern their own self-directed release from the template solely through the action of endogenous metalloproteases. Tissues recovered through this physiologically relevant process were carrier-free and structurally and phenotypically equivalent to their natural counterparts. This technology contributes to a new paradigm in regenerative medicine, whereby materials are able to actively direct and respond to cell behavior. The novel application of such materials as a coating capable of directing the formation and detachment of complex tissues solely under physiological conditions can have broad use for fundamental research and in future cell and tissue therapies.

Late Barremian-Early Aptian palaeoenvironmental change: the Cassis-La Bedoule section (SE France)

The well-dated section of Cassis-La Bédoule in the South Provencal Basin (southern France) allows for a detailed reconstruction of palaeoenvironmental change during the latest Barremian and Early Aptian. For this study, phosphorus (P) and clay-mineral contents, stable-isotope ratios on carbonate (δ13Ccarb) and organic matter (δ13Corg), and redox-sensitive trace elements (RSTE: V, U, As, Co, and Mo) have been measured in this historical stratotype. The base of the section consists of rudist limestone, which is attributed to the Urgonian platform. The presence of low P and RSTE content, and content of up to 30% kaolinite indicate deposition under oligotrophic and oxic conditions, and the presence of warm, humid climatic conditions on the adjacent continent. The top of the Urgonian succession is marked by a hardground with encrusted brachiopods and bivalves, which is interpreted as a drowning surface. The section continues with a succession of limestone and marl containing the first occurrence of planktonic foraminifera. This interval includes several laminated, organic-rich layers recording RSTE enrichments and high Corg:Ptot ratios. The deposition of these organic-rich layers was associated with oxygen-depleted conditions and a large positive excursion in δ13Corg. During this interval, a negative peak in the δ13Ccarb record is observed, which dates as latest Barremian. This excursion is coeval with negative excursions elsewhere in Tethyan platform and basin settings and is explained by the increased input of light dissolved inorganic carbon by rivers and/or volcanic activity. In this interval, an increase in P content, owing to reworking of nearshore sediments during the transgression, is coupled with a decrease in kaolinite content, which tends to be deposited in more proximal areas. The overlying hemipelagic sediments of the Early Aptian Deshayesites oglanlensis and D. weissi zones indicate rather stable palaeoenvironmental conditions with low P content and stable δ13C records. A change towards marl-dominated beds occurs close to the end of the D. weissi zone. These beds display a long decrease in their δ13Ccarb and δ13Corg records, which lasted until the end of the Deshayesites deshayesi subzone (corresponding to C3 in Menegatti et al., 1998). This is followed by a positive shift during the Roloboceras hambrovi and Deshayesites grandis subzones, which corresponds in time to oceanic anoxic event (OAE) 1a interval. This positive shift is coeval with two increases in the P content. The marly interval equivalent to OAE 1a lacks organic-rich deposits and RSTE enrichments indicating that oxic conditions prevailed in this particular part of the Tethys ocean. The clay mineralogy is dominated by smectite, which is interpreted to reflect trapping of kaolinite on the surrounding platforms rather than indicating a drier climate.

Twin experiments with the equivalent weights particle filter and HadCM3

This paper investigates the use of a particle filter for data assimilation with a full scale coupled ocean–atmosphere general circulation model. Synthetic twin experiments are performed to assess the performance of the equivalent weights filter in such a high-dimensional system. Artificial 2-dimensional sea surface temperature fields are used as observational data every day. Results are presented for different values of the free parameters in the method. Measures of the performance of the filter are root mean square errors, trajectories of individual variables in the model and rank histograms. Filter degeneracy is not observed and the performance of the filter is shown to depend on the ability to keep maximum spread in the ensemble.

Electronic structure of Pd multi-layers on Re(0001): the role of charge transfer

Understanding the origin of the properties of metal-supported metal thin films is important for the rational design of bimetallic catalysts and other applications, but it is generally difficult to separate effects related to strain from those arising from interface interactions. Here we use density functional (DFT) theory to examine the structure and electronic behavior of few-layer palladium films on the rhenium (0001) surface, where there is negligible interfacial strain and therefore other effects can be isolated. Our DFT calculations predict stacking sequences and interlayer separations in excellent agreement with quantitative low-energy electron diffraction experiments. By theoretically simulating the Pd core-level X-ray photoemission spectra (XPS) of the films, we are able to interpret and assign the basic features of both low-resolution and high-resolution XPS measurements. The core levels at the interface shift to more negative energies, rigidly following the shifts in the same direction of the valence d-band center. We demonstrate that the valence band shift at the interface is caused by charge transfer from Re to Pd, which occurs mainly to valence states of hybridized s-p character rather than to the Pd d-band. Since the d-band filling is roughly constant, there is a correlation between the d-band center shift and its bandwidth. The resulting effect of this charge transfer on the valence d-band is thus analogous to the application of a lateral compressive strain on the adlayers. Our analysis suggests that charge transfer should be considered when describing the origin of core and valence band shifts in other metal / metal adlayer systems.

Inferences about pressures and vertical extension of cloud layers from POLDER3/PARASOL measurements in the oxygen A-band

We present new inferences about cloud vertical structures from multidirectionnal measurements in the oxygen A-band. The analysis of collocated data provided by instruments onboard satellite platforms within the A-Train, as well as simulations have shown that for monolayered clouds, the cloud oxygen pressure PO2PO2 derived from the POLDER3 instrument was sensitive to the cloud vertical structure in two ways: First, PO2PO2 is actually close to the pressure of the geometrical middle of cloud and we propose a method to correct it to get the cloud top pressure (CTP), and then to obtain the cloud geometrical extent. Second, for the liquid water clouds, the angular standard deviation σPO2σPO2 of PO2PO2 is correlated with the geometrical extent of cloud layers, which makes possible a second estimation of the cloud geometrical thickness. The determination of the vertical location of cloud layers from passive measurements, eventually completed from other observations, would be useful in many applications for which cloud macrophysical properties are needed