Variation of grain nitrogen content in relation with grain yield in old and modern Spanish wheats grown under a wide range of agronomic conditions in a Mediterranean region

Wheat yield and grain nitrogen concentration (GNC; mg N/g grain) are frequently negatively correlated. In most growing conditions, this is mainly due to a feedback process between GNC and the number of grains/m2. In Mediterranean conditions, breeders may have produced cultivars with conservative grain set. The present study aimed at clarifying the main physiological determinants of grain nitrogen accumulation (GNA) in Mediterranean wheat and to analyse how breeding has affected them. Five field experiments were carried out in north-eastern Spain in the 2005/06 and 2006/ 07 growing seasons with three cultivars released at different times and an advanced line. Depending on the experiment, source-sink ratios during grain filling were altered by reducing grain number/m2 either through pre-anthesis shading (unshaded control or 0.75 shading only between jointing and anthesis) or by directly trimming the spikes after anthesis and before the onset of the effective grain filling period (un-trimmed control or spikes halved 7–10 days after anthesis). Grain nitrogen content (GN content ; mg N/grain) decreased with the year of release of the genotypes. As the number of grains/m2 was also increased by breeding there was a clear dilution effect on the amount of nitrogen allocated to each grain. However, the increase in GN content in old genotypes did not compensate for the loss in grain nitrogen yield (GNY) due to the lower number of grains/m2. GN content of all genotypes increased (increases ranged from 0.13 to 0.40 mg N/grain, depending on experiment and genotype) in response to the post-anthesis spike trimming or pre-anthesis shading. The degree of source-limitation for GNA increased with the year of release of the genotypes (and thus with increases in grain number/m2) from 0.22 (mean of the four manipulative experiments) in the oldest cultivar to 0.51 (mean of the four manipulative experiments) in the most modern line. It was found that final GN content depended strongly on the source-sink ratio established at anthesis between the number of grains set and the amount of nitrogen absorbed at this stage. Thus, Mediterranean wheat breeding that improved yield through increases in grain number/m2 reduced the GN content by diluting a rather limited source of nitrogen into more grains. This dilution effect produced by breeding was further confirmed by the reversal effect produced by grain number/m2 reductions due to either pre-anthesis shading or post-anthesis spike trimming.

Grain yield, carbon isotope discrimination and mineral content in mature kernels of barley under irrigated and rainfed conditions

The relationship between yield, carbon isotope discrimination and ash content in mature kernels was examined for a set of 13 barley (Hordeum vulgare) cultivars. Plants were grown under rainfed and well-irrigated conditions in a Mediterranean area. Water deficit caused a decrease in both grain yield and carbon isotope discrimination (Δ). The yield was positively related to Δ and negatively related to ash content, across genotypes within each treatment. However, whereas the correlation between yield and Δ was higher for the set of genotypes under well-irrigated (r=0.70, P<0.01) than under rainfed (r=0.42) conditions, the opposite occurred when yield and ash content were related, ie r=-0.38 under well-irrigated and r=-0.73, (P<0.01) under rainfed conditions. Carbon isotope discrimination and ash content together account for almost 60% of the variation in yield, in both conditions. There was no significant relationship (r=-0.15) between carbon isotope discrimination and ash content in well-irrigated plants, whereas in rainfed plants, this relationship, although significant (r=-0.54, P< 0.05), was weakly negative. The concentration of several mineral elements was measured in the same kernels. The mineral that correlated best with ash content, yield and A, was K. For yield and Δ, although the relationship with K followed the same pattern as the relationhip with ash content, the correlation coefficients were lower. Thus, mineral accumulation in mature kernels seems to be independent of transpiration efficiency. In fact, filling of grains takes place through the phloem pathway. The ash content in kernels is proposed as a complementary criterion, in addition to kernel Δ, to assess genotype differences in barley grain yield under rainfed conditions.

Growth and yield models in Spain: historical overview, contemporary examples and perspectives

En el presente trabajo se presenta una revisión sobre los modelos forestales desarrollados en España durante los últimos años, tanto para la producción maderable como no maderable y, para la dinámica de los bosques (regeneración, mortalidad). Se presentan modelos tanto de rodal completo como de clases diamétricas y de árbol individual. Los modelos desarrollados hasta la fecha se han desarrollado a partir de datos procedentes de parcelas permanentes, ensayos y el Inventario Forestal Nacional. En el trabajo se muestran los diferentes submodelos desarrollados hasta la fecha, así como las plataformas informáticas que permiten utilizar dichos modelos. Se incluyen las principales perspectivas de desarrollo de la modelización forestal en España.

Clustering of grape yield maps to delineate site-specific management zones

Zonal management in vineyards requires the prior delineation of stable yield zones within the parcel. Among the different methodologies used for zone delineation, cluster analysis of yield data from several years is one of the possibilities cited in scientific literature. However, there exist reasonable doubts concerning the cluster algorithm to be used and the number of zones that have to be delineated within a field. In this paper two different cluster algorithms have been compared (k-means and fuzzy c-means) using the grape yield data corresponding to three successive years (2002, 2003 and 2004), for a ‘Pinot Noir’ vineyard parcel. Final choice of the most recommendable algorithm has been linked to obtaining a stable pattern of spatial yield distribution and to allowing for the delineation of compact and average sized areas. The general recommendation is to use reclassified maps of two clusters or yield classes (low yield zone and high yield zone) and, consequently, the site-specific vineyard management should be based on the prior delineation of just two different zones or sub-parcels. The two tested algorithms are good options for this purpose. However, the fuzzy c-means algorithm allows for a better zoning of the parcel, forming more compact areas and with more equilibrated zonal differences over time.

A semi-grand canonical Monte Carlo simulation model for ion binding to ionizable surfaces: Proton binding of carboxylated latex particles as a case study

In this paper, we present a computer simulation study of the ion binding process at an ionizable surface using a semi-grand canonical Monte Carlo method that models the surface as a discrete distribution of charged and neutral functional groups in equilibrium with explicit ions modelled in the context of the primitive model. The parameters of the simulation model were tuned and checked by comparison with experimental titrations of carboxylated latex particles in the presence of different ionic strengths of monovalent ions. The titration of these particles was analysed by calculating the degree of dissociation of the latex functional groups vs. pH curves at different background salt concentrations. As the charge of the titrated surface changes during the simulation, a procedure to keep the electroneutrality of the system is required. Here, two approaches are used with the choice depending on the ion selected to maintain electroneutrality: counterion or coion procedures. We compare and discuss the difference between the procedures. The simulations also provided a microscopic description of the electrostatic double layer (EDL) structure as a function of pH and ionic strength. The results allow us to quantify the effect of the size of the background salt ions and of the surface functional groups on the degree of dissociation. The non-homogeneous structure of the EDL was revealed by plotting the counterion density profiles around charged and neutral surface functional groups. © 2011 American Institute of Physics.

Comblike poly(a-alkyl g glutamate)s: computer simulation studies of an intermediate thermal phase.

Monte Carlo (MC) simulations have been used to study the structure of an intermediate thermal phase of poly(R-octadecyl ç,D-glutamate). This is a comblike poly(ç-peptide) able to adopt a biphasic structure that has been described as a layered arrangement of backbone helical rods immersed in a paraffinic pool of polymethylene side chains. Simulations were performed at two different temperatures (348 and 363 K), both of them above the melting point of the paraffinic phase, using the configurational bias MC algorithm. Results indicate that layers are constituted by a side-by-side packing of 17/5 helices. The organization of the interlayer paraffinic region is described in atomistic terms by examining the torsional angles and the end-to-end distances for the octadecyl side chains. Comparison with previously reported comblike poly(â-peptide)s revealed significant differences in the organization of the alkyl side chains.

Cluster glasses of ultrasoft particles

We present molecular dynamics (MD) simulations results for dense fluids of ultrasoft, fully penetrable particles. These are a binary mixture and a polydisperse system of particles interacting via the generalized exponential model, which is known to yield cluster crystal phases for the corresponding monodisperse systems. Because of the dispersity in the particle size, the systems investigated in this work do not crystallize and form disordered cluster phases. The clusteringtransition appears as a smooth crossover to a regime in which particles are mostly located in clusters, isolated particles being infrequent. The analysis of the internal cluster structure reveals microsegregation of the big and small particles, with a strong homo-coordination in the binary mixture. Upon further lowering the temperature below the clusteringtransition, the motion of the clusters" centers-of-mass slows down dramatically, giving way to a cluster glass transition. In the cluster glass, the diffusivities remain finite and display an activated temperature dependence, indicating that relaxation in the cluster glass occurs via particle hopping in a nearly arrested matrix of clusters. Finally we discuss the influence of the microscopic dynamics on the transport properties by comparing the MD results with Monte Carlo simulations.

Disentangling the formation of contrasting tree line physiognomies combining model selection and Bayesian parameterization for simulation models.

Alpine tree-line ecotones are characterized by marked changes at small spatial scales that may result in a variety of physiognomies. A set of alternative individual-based models was tested with data from four contrasting Pinus uncinata ecotones in the central Spanish Pyrenees to reveal the minimal subset of processes required for tree-line formation. A Bayesian approach combined with Markov chain Monte Carlo methods was employed to obtain the posterior distribution of model parameters, allowing the use of model selection procedures. The main features of real tree lines emerged only in models considering nonlinear responses in individual rates of growth or mortality with respect to the altitudinal gradient. Variation in tree-line physiognomy reflected mainly changes in the relative importance of these nonlinear responses, while other processes, such as dispersal limitation and facilitation, played a secondary role. Different nonlinear responses also determined the presence or absence of krummholz, in agreement with recent findings highlighting a different response of diffuse and abrupt or krummholz tree lines to climate change. The method presented here can be widely applied in individual-based simulation models and will turn model selection and evaluation in this type of models into a more transparent, effective, and efficient exercise.

Understanding Optical Trapping Phenomena: a Simulation for Undergraduates

Optical trapping is an attractive and multidisciplinary topic that has become the center of attention to a large number of researchers. Moreover, it is a suitable subject for advanced students that requires a knowledge of a wide range of topics. As a result, it has been incorporated into some syllabuses of both undergraduate and graduate programs. In this paper, basic concepts in laser trapping theory are reviewed. To provide a better understanding of the underlying concepts for students, a Java application for simulating the behavior of a dielectric particle trapped in a highly focused beam has been developed. The program illustrates a wide range of theoretical results and features, such as the calculation of the force exerted by a beam in the Mie and Rayleigh regimes or the calibration of the trap stiffness. Some examples that are ready to be used in the classroom or in the computer lab are also supplied.

Flocking behaviour: Agent-based simulation and hierarchical leadership.

We have studied how leaders emerge in a group as a consequence of interactions among its members. We propose that leaders can emerge as a consequence of a self-organized process based on local rules of dyadic interactions among individuals. Flocks are an example of self-organized behaviour in a group and properties similar to those observed in flocks might also explain some of the dynamics and organization of human groups. We developed an agent-based model that generated flocks in a virtual world and implemented it in a multi-agent simulation computer program that computed indices at each time step of the simulation to quantify the degree to which a group moved in a coordinated way (index of flocking behaviour) and the degree to which specific individuals led the group (index of hierarchical leadership). We ran several series of simulations in order to test our model and determine how these indices behaved under specific agent and world conditions. We identified the agent, world property, and model parameters that made stable, compact flocks emerge, and explored possible environmental properties that predicted the probability of becoming a leader.

Molecular dynamics simulation of the spherical electrical double layer of a soft nanoparticle. Effect of the surface and counterion valence

Molecular dynamics simulations were performed to study the ion and water distribution around a spherical charged nanoparticle. A soft nanoparticle model was designed using a set of hydrophobic interaction sites distributed in six concentric spherical layers. In order to simulate the effect of charged functionalyzed groups on the nanoparticle surface, a set of charged sites were distributed in the outer layer. Four charged nanoparticle models, from a surface charge value of −0.035 Cm−2 to − 0.28 Cm−2, were studied in NaCl and CaCl2 salt solutions at 1 M and 0.1 M concentrations to evaluate the effect of the surface charge, counterion valence, and concentration of added salt. We obtain that Na + and Ca2 + ions enter inside the soft nanoparticle. Monovalent ions are more accumulated inside the nanoparticle surface, whereas divalent ions are more accumulated just in the plane of the nanoparticle surface sites. The increasing of the the salt concentration has little effect on the internalization of counterions, but significantly reduces the number of water molecules that enter inside the nanoparticle. The manner of distributing the surface charge in the nanoparticle (uniformly over all surface sites or discretely over a limited set of randomly selected sites) considerably affects the distribution of counterions in the proximities of the nanoparticle surface.

Effect of the surface charge discretization on electric double layers. A Monte Carlo simulation study

The structure of the electric double layer in contact with discrete and continuously charged planar surfaces is studied within the framework of the primitive model through Monte Carlo simulations. Three different discretization models are considered together with the case of uniform distribution. The effect of discreteness is analyzed in terms of charge density profiles. For point surface groups,a complete equivalence with the situation of uniformly distributed charge is found if profiles are exclusively analyzed as a function of the distance to the charged surface. However, some differences are observed moving parallel to the surface. Significant discrepancies with approaches that do not account for discreteness are reported if charge sites of finite size placed on the surface are considered.

Subsurface drip irrigation emitter spacing effects on soil water redistribution, corn yield and water productivity

Emitter spacings of 0.3 to 0.6 m are commonly used for subsurface drip irrigation (SDI) of corn on the deep, silt loam soils of the U.S. Great Plains. Subsurface drip irrigation emitter spacings of 0.3, 0.6, 0.9 and 1.2 m were examined for the resulting differences in soil water redistribution, corn grain yield, yield components, seasonal water use, and water productivity in a 4‐year field study (2005 through 2008) at the Kansas State University Northwest Research‐Extension Center, Colby, Kansas. The results indicate that there is increased preferential water movement along the dripline (parallel) as compared to perpendicular to the dripline and that this phenomenon partially compensates for wider emitter spacings in terms of soil water redistribution. Corn yield and water productivity (WP) were not significantly affected by the emitter spacing with application of a full irrigation regime

A measure of stability as a criterion for the verification and analysis of simulation models

The aim of this study is to define a new statistic, PVL, based on the relative distance between the likelihood associated with the simulation replications and the likelihood of the conceptual model. Our results coming from several simulation experiments of a clinical trial show that the PVL statistic range can be a good measure of stability to establish when a computational model verifies the underlying conceptual model. PVL improves also the analysis of simulation replications because only one statistic is associated with all the simulation replications. As well it presents several verification scenarios, obtained by altering the simulation model, that show the usefulness of PVL. Further simulation experiments suggest that a 0 to 20 % range may define adequate limits for the verification problem, if considered from the viewpoint of an equivalence test.

Simulation of a Virtual Campus at the Open University of Catalonia

Peer-reviewed