A joint state and parameter estimation scheme for nonlinear dynamical systems

We present a novel algorithm for concurrent model state and parameter estimation in nonlinear dynamical systems. The new scheme uses ideas from three dimensional variational data assimilation (3D-Var) and the extended Kalman filter (EKF) together with the technique of state augmentation to estimate uncertain model parameters alongside the model state variables in a sequential filtering system. The method is relatively simple to implement and computationally inexpensive to run for large systems with relatively few parameters. We demonstrate the efficacy of the method via a series of identical twin experiments with three simple dynamical system models. The scheme is able to recover the parameter values to a good level of accuracy, even when observational data are noisy. We expect this new technique to be easily transferable to much larger models.

Pension scheme redesign and wealth redistribution between the members and sponsor: the USS rule change in October 2011

The redesign of defined benefit pension schemes usually results in a substantial redistribution of wealth between age cohorts of members, pensioners, and the sponsor. This is the first study to quantify the redistributive effects of a rule change by a real world scheme (the Universities Superannuation Scheme, USS) where the sponsor underwrites the pension promise. In October 2011 USS closed its final salary scheme to new members, opened a career average revalued earnings (CARE) section, and moved to ‘cap and share’ contribution rates. We find that the pre-October 2011 scheme was not viable in the long run, while the post-October 2011 scheme is probably viable in the long run, but faces medium term problems. In October 2011 future members of USS lost 65% of their pension wealth (or roughly £100,000 per head), equivalent to a reduction of roughly 11% in their total compensation, while those aged over 57 years lost almost nothing. The riskiness of the pension wealth of future members increased by a third, while the riskiness of the present value of the sponsor’s future contributions reduced by 10%. Finally, the sponsor’s wealth increased by about £32.5 billion, equivalent to a reduction of 26% in their pension costs.

Surface Urban Energy and Water Balance Scheme (SUEWS): development and evaluation at two UK sites

The Surface Urban Energy and Water Balance Scheme (SUEWS) is evaluated at two locations in the UK: a dense urban site in the centre of London and a residential suburban site in Swindon. Eddy covariance observations of the turbulent fluxes are used to assess model performance over a twoyear period (2011-2013). The distinct characteristics of the sites mean their surface energy exchanges differ considerably. The model suggests the largest differences can be attributed to surface cover (notably the proportion of vegetated versus impervious area) and the additional energy supplied by human activities. SUEWS performs better in summer than winter, and better at the suburban site than the dense urban site. One reason for this is the bias towards suburban summer field campaigns in observational data used to parameterise this (and other) model(s). The suitability of model parameters (such as albedo, energy use and water use) for the UK sites is considered and, where appropriate, alternative values are suggested. An alternative parameterisation for the surface conductance is implemented, which permits greater soil moisture deficits before evaporation is restricted at non-irrigated sites. Accounting for seasonal variation in the estimation of storage heat flux is necessary to obtain realistic wintertime fluxes.

Evaluation of the Plant–Craig stochastic convection scheme (v2.0) in the ensemble forecasting system MOGREPS-R (24 km) based on the Unified Model (v7.3)

The Plant–Craig stochastic convection parameterization (version 2.0) is implemented in the Met Office Regional Ensemble Prediction System (MOGREPS-R) and is assessed in comparison with the standard convection scheme with a simple stochastic scheme only, from random parameter variation. A set of 34 ensemble forecasts, each with 24 members, is considered, over the month of July 2009. Deterministic and probabilistic measures of the precipitation forecasts are assessed. The Plant–Craig parameterization is found to improve probabilistic forecast measures, particularly the results for lower precipitation thresholds. The impact on deterministic forecasts at the grid scale is neutral, although the Plant–Craig scheme does deliver improvements when forecasts are made over larger areas. The improvements found are greater in conditions of relatively weak synoptic forcing, for which convective precipitation is likely to be less predictable.

Mobile sink based fault diagnosis scheme for wireless sensor networks

Network diagnosis in Wireless Sensor Networks (WSNs) is a difficult task due to their improvisational nature, invisibility of internal running status, and particularly since the network structure can frequently change due to link failure. To solve this problem, we propose a Mobile Sink (MS) based distributed fault diagnosis algorithm for WSNs. An MS, or mobile fault detector is usually a mobile robot or vehicle equipped with a wireless transceiver that performs the task of a mobile base station while also diagnosing the hardware and software status of deployed network sensors. Our MS mobile fault detector moves through the network area polling each static sensor node to diagnose the hardware and software status of nearby sensor nodes using only single hop communication. Therefore, the fault detection accuracy and functionality of the network is significantly increased. In order to maintain an excellent Quality of Service (QoS), we employ an optimal fault diagnosis tour planning algorithm. In addition to saving energy and time, the tour planning algorithm excludes faulty sensor nodes from the next diagnosis tour. We demonstrate the effectiveness of the proposed algorithms through simulation and real life experimental results.

Modeling study of the aspect ratio influence on urban canopy energy fluxes with a modified wall-canyon energy budget scheme

The influence of the aspect ratio (building height/street canyon width) and the mean building height of cities on local energy fluxes and temperatures is studied by means of an Urban Canopy Model (UCM) coupled with a one-dimensional second-order turbulence closure model. The UCM presented is similar to the Town Energy Balance (TEB) model in most of its features but differs in a few important aspects. In particular, the street canyon walls are treated separately which leads to a different budget of radiation within the street canyon walls. The UCM has been calibrated using observations of incoming global and diffuse solar radiation, incoming long-wave radiation and air temperature at a site in So Paulo, Brazil. Sensitivity studies with various aspect ratios have been performed to assess their impact on urban temperatures and energy fluxes at the top of the canopy layer. In these simulations, it is assumed that the anthropogenic heat flux and latent heat fluxes are negligible. Results show that the simulated net radiation and sensible heat fluxes at the top of the canopy decrease and the stored heat increases as the aspect ratio increases. The simulated air temperature follows the behavior of the sensible heat flux. (C) 2010 Elsevier Ltd. All rights reserved.

Formulation of a tropical town energy budget (t-TEB) scheme

This work describes the tropical town energy budget (t-TEB) scheme addressed to simulate the diurnal occurrence of the urban heat island (UHI) as observed in the Metropolitan Area of Rio de Janeiro (MARJ; -22A degrees S; -44A degrees W) in Brazil. Reasoning about the tropical urban climate have guided the scheme implementation, starting from the original equations from Masson (Bound-Lay Meteorol 94:357-397, 2000). The modifications include (a) local scaling approaches for obtaining flux-gradient relationships in the roughness sub-layer, (b) the Monin-Obukhov similarity framework in the inertial sub-layer, (c) increasing aerodynamic conductance toward more unstable conditions, and (d) a modified urban subsurface drainage system to transfer the intercepted rainwater by roofs to the roads. Simulations along 2007 for the MARJ are obtained and compared with the climatology. The t-TEB simulation is consistent with the observations, suggesting that the timing and dynamics of the UHI in tropical cities could vary significantly from the familiar patterns observed in mid-latitude cities-with the peak heat island intensity occurring in the morning than at night. The simulations are suggesting that the thermal phase shift of this tropical diurnal UHI is a response of the surface energy budget to the large amount of solar radiation, intense evapotranspiration, and thermal response of the vegetated surfaces over a very humid soil layer.

Habitat Split as a Cause of Local Population Declines of Amphibians with Aquatic Larvae

Most amphibian species have biphasic life histories and undergo an ontogenetic shift from aquatic to terrestrial habitats. In deforested landscapes, streams and forest fragments are frequently disjunct, jeopardizing the life cycle of forest-associated amphibians with aquatic larvae. We tested the impact of habitat split-defined as human-induced disconnection between habitats used by different life-history stages of a species-on four forest-associated amphibian species in a severely fragmented landscape of the Brazilian Atlantic Forest. We surveyed amphibians in forest fragments with and without streams (referred to as wet and dry fragments, respectively), including the adjacent grass-field matrix. Our comparison of capture rates in dry fragments and nearby streams in the matrix allowed us to evaluate the number of individuals that engaged in high-risk migrations through nonforested habitats. Adult amphibians moved from dry fragments to matrix streams at the beginning of the rainy season, reproduced, and returned at the end of the breeding period. Juveniles of the year moved to dry fragments along with adults. These risky reproductive migrations through nonforested habitats that expose individuals to dehydration, predation, and other hazards may cause population declines in dry fragments. Indeed, capture rates were significantly lower in dry fragments compared with wet fragments. Declining amphibians would strongly benefit from investments in the conservation and restoration of riparian vegetation and corridors linking breeding and nonbreeding areas.

A novel locus for split-hand/foot malformation associated with tibial hemimelia (SHFLD syndrome) maps to chromosome region 17p13.1-17p13.3

Split-hand/foot malformation (SHFM) associated with aplasia of long bones, SHFLD syndrome or Tibial hemimelia-ectrodactyly syndrome is a rare condition with autosomal dominant inheritance, reduced penetrance and an incidence estimated to be about 1 in 1,000,000 liveborns. To date, three chromosomal regions have been reported as strong candidates for harboring SHFLD syndrome genes: 1q42.2-q43, 6q14.1 and 2q14.2. We characterized the phenotype of nine affected individuals from a large family with the aim of mapping the causative gene. Among the nine affected patients, four had only SHFM of the hands and no tibial defects, three had both defects and two had only unilateral tibial hemimelia. In keeping with previous publications of this and other families, there was clear evidence of both variable expression and incomplete penetrance, the latter bearing hallmarks of anticipation. Segregation analysis and multipoint Lod scores calculations (maximum Lod score of 5.03 using the LINKMAP software) using all potentially informative family members, both affected and unaffected, identified the chromosomal region 17p13.1-17p13.3 as the best and only candidate for harboring a novel mutated gene responsible for the syndrome in this family. The candidate gene CRK located within this region was sequenced but no pathogenic mutation was detected.

A geometric mass-preserving redistancing scheme for the level set function

In this paper we describe and evaluate a geometric mass-preserving redistancing procedure for the level set function on general structured grids. The proposed algorithm is adapted from a recent finite element-based method and preserves the mass by means of a localized mass correction. A salient feature of the scheme is the absence of adjustable parameters. The algorithm is tested in two and three spatial dimensions and compared with the widely used partial differential equation (PDE)-based redistancing method using structured Cartesian grids. Through the use of quantitative error measures of interest in level set methods, we show that the overall performance of the proposed geometric procedure is better than PDE-based reinitialization schemes, since it is more robust with comparable accuracy. We also show that the algorithm is well-suited for the highly stretched curvilinear grids used in CFD simulations. Copyright (C) 2010 John Wiley & Sons, Ltd.

Evaluation of a bounded high order upwind scheme for 3D incompressible free surface flow computations

fit the context of normalized variable formulation (NVF) of Leonard and total variation diminishing (TVD) constraints of Harten. this paper presents an extension of it previous work by the authors for solving unsteady incompressible flow problems. The main contributions of the paper are threefold. First, it presents the results of the development and implementation of a bounded high order upwind adaptative QUICKEST scheme in the 3D robust code (Freeflow), for the numerical solution of the full incompressible Navier-Stokes equations. Second, it reports numerical simulation results for 1D hock tube problem, 2D impinging jet and 2D/3D broken clam flows. Furthermore, these results are compared with existing analytical and experimental data. And third, it presents the application of the numerical method for solving 3D free surface flow problems. (C) 2007 IMACS. Published by Elsevier B.V. All rights reserved,

Direct Observation of Tetragonal Distortion in Epitaxial Structures through Secondary Peak Split in a Synchrotron Radiation Renninger Scan

This paper reports a direct observation of an interesting split of the (022)(022) four-beam secondary peak into two (022) and (022) three-beam peaks, in a synchrotron radiation Renninger scan (phi-scan), as an evidence of the layer tetragonal distortion in two InGaP/GaAs (001) epitaxial structures with different thicknesses. The thickness, composition, (a perpendicular to) perpendicular lattice parameter, and (01) in-plane lattice parameter of the two epitaxial ternary layers were obtained from rocking curves (omega-scan) as well as from the simulation of the (022)(022) split, and then, it allowed for the determination of the perpendicular and parallel (in-plane) strains. Furthermore, (022)(022) omega:phi mappings were measured in order to exhibit the multiple diffraction condition of this four-beam case with their split measurement.

Near threshold vibrational excitation of molecules by positron impact: A projection operator approach

We report vibrational excitation (v(i) = 0 -> v(f) = 1) cross-sections for positron scattering by H(2) and model calculations for the (v(i) = 0 -> v(f) = 1) excitation of the C-C symmetric stretch mode of C(2)H(2). The Feshbach projection operator formalism was employed to vibrationally resolve the fixed-nuclei phase shifts obtained with the Schwinger multichannel method. The near threshold behavior of H(2) and C(2)H(2) significantly differ in the sense that no low lying singularity (either virtual or bound state) was found for the former, while a e(+)-acetylene virtual state was found at the equilibrium geometry (this virtual state becomes a bound state upon stretching the molecule). For C(2)H(2), we also performed model calculations comparing excitation cross-sections arising from virtual (-i kappa(0)) and bound (+i kappa(0)) states symmetrically located around the origin of the complex momentum plane (i.e. having the same kappa(0)). The virtual state is seen to significantly couple to vibrations, and similar cross-sections were obtained for shallow bound and virtual states. (c) 2007 Elsevier B.V. All rights reserved.

A Direct Numerov Sixth-order Numerical Scheme to Accurately Solve the Unidimensional Poisson Equation with Dirichlet Boundary Conditions

In this article, we present an analytical direct method, based on a Numerov three-point scheme, which is sixth order accurate and has a linear execution time on the grid dimension, to solve the discrete one-dimensional Poisson equation with Dirichlet boundary conditions. Our results should improve numerical codes used mainly in self-consistent calculations in solid state physics.