Impact of mass lumping on gravity and Rossby waves in 2D finite-element shallow-water models

The goal of this study is to evaluate the effect of mass lumping on the dispersion properties of four finite-element velocity/surface-elevation pairs that are used to approximate the linear shallow-water equations. For each pair, the dispersion relation, obtained using the mass lumping technique, is computed and analysed for both gravity and Rossby waves. The dispersion relations are compared with those obtained for the consistent schemes (without lumping) and the continuous case. The P0-P1, RT0 and P-P1 pairs are shown to preserve good dispersive properties when the mass matrix is lumped. Test problems to simulate fast gravity and slow Rossby waves are in good agreement with the analytical results.

Influence of freezable/non-freezable water and sucrose on the viability of Theobroma cacao somatic embryos following desiccation and freezing

Encapsulated cocoa (Theobroma cacao L.) somatic embryos subjected to 0.08-1.25 M sucrose treatments were analyzed for embryo soluble sugar content, non-freezable water content, moisture level after desiccation and viability after desiccation and freezing. Results indicated that the higher the sucrose concentration in the treatment medium, the greater was the extent of sucrose accumulation in the embryos. Sucrose treatment greatly assisted embryo post-desiccation recovery since only 40% of the control embryos survived desiccation, whereas a survival rate of 60-95% was recorded for embryos exposed to 0.5-1.25 M sucrose. The non-freezable water content of the embryos was estimated at between 0.26 and 0.61 g H2O g(-1)dw depending on the sucrose treatment, and no obvious relationship could be found between the endogenous sucrose level and the amount of non-freezable water in the embryos. Cocoa somatic embryos could withstand the loss of a fraction of their non-freezable water without losing viability following desiccation. Nevertheless, the complete removal of potentially freezable water was not sufficient for most embryos to survive freezing.

Temperature dependence of protein dynamics simulated with three different water models

The effect of variation of the water model on the temperature dependence of protein and hydration water dynamics is examined by performing molecular dynamics simulations of myoglobin with the TIP3P, TIP4P, and TIP5P water models and the CHARMM protein force field at temperatures between 20 and 300 K. The atomic mean-square displacements, solvent reorientational relaxation times, pair angular correlations between surface water molecules, and time-averaged structures of the protein are all found to be similar, and the protein dynamical transition is described almost indistinguishably for the three water potentials. The results provide evidence that for some purposes changing the water model in protein simulations without a loss of accuracy may be possible.

Forecasting economic and financial time series with non-linear models

In this paper we discuss the current state-of-the-art in estimating, evaluating, and selecting among non-linear forecasting models for economic and financial time series. We review theoretical and empirical issues, including predictive density, interval and point evaluation and model selection, loss functions, data-mining, and aggregation. In addition, we argue that although the evidence in favor of constructing forecasts using non-linear models is rather sparse, there is reason to be optimistic. However, much remains to be done. Finally, we outline a variety of topics for future research, and discuss a number of areas which have received considerable attention in the recent literature, but where many questions remain.

Molecular dynamics simulations of proteins: can the explicit water model be varied?

In molecular mechanics simulations of biological systems, the solvation water is typically represented by a default water model which is an integral part of the force field. Indeed, protein nonbonding parameters are chosen in order to obtain a balance between water-water and protein-water interactions and hence a reliable description of protein solvation. However, less attention has been paid to the question of whether the water model provides a reliable description of the water properties under the chosen simulation conditions, for which more accurate water models often exist. Here we consider the case of the CHARMM protein force field, which was parametrized for use with a modified TIP3P model. Using quantum mechanical and molecular mechanical calculations, we investigate whether the CHARMM force field can be used with other water models: TIP4P and TIP5P. Solvation properties of N-methylacetamide (NMA), other small solute molecules, and a small protein are examined. The results indicate differences in binding energies and minimum energy geometries, especially for TIP5P, but the overall description of solvation is found to be similar for all models tested. The results provide an indication that molecular mechanics simulations with the CHARMM force field can be performed with water models other than TIP3P, thus enabling an improved description of the solvent water properties.

Measures of observation impact in non-Gaussian data assimilation

ABSTRACT Non-Gaussian/non-linear data assimilation is becoming an increasingly important area of research in the Geosciences as the resolution and non-linearity of models are increased and more and more non-linear observation operators are being used. In this study, we look at the effect of relaxing the assumption of a Gaussian prior on the impact of observations within the data assimilation system. Three different measures of observation impact are studied: the sensitivity of the posterior mean to the observations, mutual information and relative entropy. The sensitivity of the posterior mean is derived analytically when the prior is modelled by a simplified Gaussian mixture and the observation errors are Gaussian. It is found that the sensitivity is a strong function of the value of the observation and proportional to the posterior variance. Similarly, relative entropy is found to be a strong function of the value of the observation. However, the errors in estimating these two measures using a Gaussian approximation to the prior can differ significantly. This hampers conclusions about the effect of the non-Gaussian prior on observation impact. Mutual information does not depend on the value of the observation and is seen to be close to its Gaussian approximation. These findings are illustrated with the particle filter applied to the Lorenz ’63 system. This article is concluded with a discussion of the appropriateness of these measures of observation impact for different situations.

A very rare self-assembled ribbon of fused cyclic water pentamers encapsulated in a Cu-II based metal-organic framework

A metal organic framework of Cu-II, tartarate (tar) and 2,2'-bipyridyl (2,2'-bipy)], {[Cu(tar)(2,2'-bipy)]center dot 5H(2)O}(n)} (1) has been synthesized at the mild ambient condition and characterized by single crystal X-ray crystallography. In the compound, the Cu(2,2'-bipy) entities are bridged by tartarate ions which are coordinated to Cu-II by both hydroxyl and monodentate carboxylate oxygen to form a one-dimensional chain. The non-coordinated water molecules form ID water chains by edge-sharing cyclic water pentamers along with dangling water dimers. It shows reversible water expulsion upon heating. The water chains join the ID coordination polymeric chains to a 31) network through hydrogen-bond interactions.

The role of information on farmers' willingness to use recycled water for irrigation

In this paper are identified several factors which affect a potential user's willingness to use recycled water for agricultural irrigation. This study is based on the results of a survey carried out among farmers in the island of Crete, Greece. It was found that a higher level of income and education are positively correlated with a respondent's willingness to use recycled water. Income and education are also positively correlated with a potential user's sensitivity to information on the advantages of using non-conventional water resources. Overall, extra information on the advantages of recycled water has a statistically significant impact on reported degrees of willingness to use recycled water.

Evaluation of water properties in HEA–HEMA hydrogels swollen in aqueous-PEG solutions using thermoanalytical techniques

Hydrogels are polymeric materials used in many pharmaceutical and biomedical applications due to their ability to form 3D hydrophilic polymeric networks, which can absorb large amounts of water. In the present work, polyethylene glycols (PEG) were introduced into the hydrogel liquid phase in order to improve the mechanical properties of hydrogels composed of 2-hydroxyethylacrylate and 2-hydroxyethylmethacrylate (HEA–HEMA) synthesized with different co-monomer compositions and equilibrated in water or in 20 % water–PEG 400 and 600 solutions. The thermoanalytical techniques [differential scanning calorimetry (DSC) and thermogravimetry (TG)] were used to evaluate the amount and properties of free and bound water in HEA–HEMA hydrogels. The internal structure and the mechanical properties of hydrogels were studied using scanning electron microscopy and friability assay. TG “loss-on-drying” experiments were applied to study the water-retention properties of hydrogels, whereas the combination of TG and DSC allowed estimating the total amount of freezable and non-freezing water in hydrogels. The results show that the addition of viscous co-solvent (PEG) to the liquid medium results in significant improvement of the mechanical properties of HEA–HEMA hydrogels and also slightly retards the water loss from the hydrogels. A redistribution of free and bound water in the hydrogels equilibrated in mixed solutions containing 20 vol% of PEGs takes place.

Evaluating a new LISFLOOD-FP formulation with data from the summer 2007 floods in Tewkesbury, UK

Two-dimensional flood inundation modelling is a widely used tool to aid flood risk management. In urban areas, where asset value and population density are greatest, the model spatial resolution required to represent flows through a typical street network (i.e. < 10m) often results in impractical computational cost at the whole city scale. Explicit diffusive storage cell models become very inefficient at such high resolutions, relative to shallow water models, because the stable time step in such schemes scales as a quadratic of resolution. This paper presents the calibration and evaluation of a recently developed new formulation of the LISFLOOD-FP model, where stability is controlled by the Courant–Freidrichs–Levy condition for the shallow water equations, such that, the stable time step instead scales linearly with resolution. The case study used is based on observations during the summer 2007 floods in Tewkesbury, UK. Aerial photography is available for model evaluation on three separate days from the 24th to the 31st of July. The model covered a 3.6 km by 2 km domain and was calibrated using gauge data from high flows during the previous month. The new formulation was benchmarked against the original version of the model at 20 m and 40 m resolutions, demonstrating equally accurate performance given the available validation data but at 67x faster computation time. The July event was then simulated at the 2 m resolution of the available airborne LiDAR DEM. This resulted in a significantly more accurate simulation of the drying dynamics compared to that simulated by the coarse resolution models, although estimates of peak inundation depth were similar.

Probabilistic wind power forecasts with an inverse power curve transformation and censored regression

Forecasting wind power is an important part of a successful integration of wind power into the power grid. Forecasts with lead times longer than 6 h are generally made by using statistical methods to post-process forecasts from numerical weather prediction systems. Two major problems that complicate this approach are the non-linear relationship between wind speed and power production and the limited range of power production between zero and nominal power of the turbine. In practice, these problems are often tackled by using non-linear non-parametric regression models. However, such an approach ignores valuable and readily available information: the power curve of the turbine's manufacturer. Much of the non-linearity can be directly accounted for by transforming the observed power production into wind speed via the inverse power curve so that simpler linear regression models can be used. Furthermore, the fact that the transformed power production has a limited range can be taken care of by employing censored regression models. In this study, we evaluate quantile forecasts from a range of methods: (i) using parametric and non-parametric models, (ii) with and without the proposed inverse power curve transformation and (iii) with and without censoring. The results show that with our inverse (power-to-wind) transformation, simpler linear regression models with censoring perform equally or better than non-linear models with or without the frequently used wind-to-power transformation.

Towards an improved and more flexible representation of water stress in coupled photosynthesis–stomatal conductance models

Coupled photosynthesis–stomatal conductance (A–gs) models are commonly used in ecosystem models to represent the exchange rate of CO2 and H2O between vegetation and the atmosphere. The ways these models account for water stress differ greatly among modelling schemes. This study provides insight into the impact of contrasting model configurations of water stress on the simulated leaf-level values of net photosynthesis (A), stomatal conductance (gs), the functional relationship among them and their ratio, the intrinsic water use efficiency (A/gs), as soil dries. A simple, yet versatile, normalized soil moisture dependent function was used to account for the effects of water stress on gs, on mesophyll conductance (gm) and on the biochemical capacity. Model output was compared to leaf-level values obtained from the literature. The sensitivity analyses emphasized the necessity to combine both stomatal and non-stomatal limitations of A in coupled A–gs models to accurately capture the observed functional relationships A vs. gs and A/gsvs. gs in response to drought. Accounting for water stress in coupled A–gs models by imposing either stomatal or biochemical limitations of A, as commonly practiced in most ecosystem models, failed to reproduce the observed functional relationship between key leaf gas exchange attributes. A quantitative limitation analysis revealed that the general pattern of C3 photosynthetic response to water stress may be well represented in coupled A–gs models by imposing the highest limitation strength to gm, then to gs and finally to the biochemical capacity.