The DYMECS project: a statistical approach for the evaluation of convective storms in high-resolution NWP models

A new frontier in weather forecasting is emerging by operational forecast models now being run at convection-permitting resolutions at many national weather services. However, this is not a panacea; significant systematic errors remain in the character of convective storms and rainfall distributions. The DYMECS project (Dynamical and Microphysical Evolution of Convective Storms) is taking a fundamentally new approach to evaluate and improve such models: rather than relying on a limited number of cases, which may not be representative, we have gathered a large database of 3D storm structures on 40 convective days using the Chilbolton radar in southern England. We have related these structures to storm life-cycles derived by tracking features in the rainfall from the UK radar network, and compared them statistically to storm structures in the Met Office model, which we ran at horizontal grid length between 1.5 km and 100 m, including simulations with different subgrid mixing length. We also evaluated the scale and intensity of convective updrafts using a new radar technique. We find that the horizontal size of simulated convective storms and the updrafts within them is much too large at 1.5-km resolution, such that the convective mass flux of individual updrafts can be too large by an order of magnitude. The scale of precipitation cores and updrafts decreases steadily with decreasing grid lengths, as does the typical storm lifetime. The 200-m grid-length simulation with standard mixing length performs best over all diagnostics, although a greater mixing length improves the representation of deep convective storms.

Modeling plant transpiration under limited soil water: comparison of different plant and soil hydraulic parameterizations and preliminary implications for their use in land surface models

Accurate estimates of how soil water stress affects plant transpiration are crucial for reliable land surface model (LSM) predictions. Current LSMs generally use a water stress factor, β, dependent on soil moisture content, θ, that ranges linearly between β = 1 for unstressed vegetation and β = 0 when wilting point is reached. This paper explores the feasibility of replacing the current approach with equations that use soil water potential as their independent variable, or with a set of equations that involve hydraulic and chemical signaling, thereby ensuring feedbacks between the entire soil–root–xylem–leaf system. A comparison with the original linear θ-based water stress parameterization, and with its improved curvi-linear version, was conducted. Assessment of model suitability was focused on their ability to simulate the correct (as derived from experimental data) curve shape of relative transpiration versus fraction of transpirable soil water. We used model sensitivity analyses under progressive soil drying conditions, employing two commonly used approaches to calculate water retention and hydraulic conductivity curves. Furthermore, for each of these hydraulic parameterizations we used two different parameter sets, for 3 soil texture types; a total of 12 soil hydraulic permutations. Results showed that the resulting transpiration reduction functions (TRFs) varied considerably among the models. The fact that soil hydraulic conductivity played a major role in the model that involved hydraulic and chemical signaling led to unrealistic values of β, and hence TRF, for many soil hydraulic parameter sets. However, this model is much better equipped to simulate the behavior of different plant species. Based on these findings, we only recommend implementation of this approach into LSMs if great care with choice of soil hydraulic parameters is taken

The utility and application of mixed-effects models in second language research

Second language acquisition researchers often face particular challenges when attempting to generalize study findings to the wider learner population. For example, language learners constitute a heterogeneous group, and it is not always clear how a study’s findings may generalize to other individuals who may differ in terms of language background and proficiency, among many other factors. In this paper, we provide an overview of how mixed-effects models can be used to help overcome these and other issues in the field of second language acquisition. We provide an overview of the benefits of mixed-effects models and a practical example of how mixed-effects analyses can be conducted. Mixed-effects models provide second language researchers with a powerful statistical tool in the analysis of a variety of different types of data.

Calibration and evaluation of individual-based models using Approximate Bayesian Computation

This paper investigates the feasibility of using approximate Bayesian computation (ABC) to calibrate and evaluate complex individual-based models (IBMs). As ABC evolves, various versions are emerging, but here we only explore the most accessible version, rejection-ABC. Rejection-ABC involves running models a large number of times, with parameters drawn randomly from their prior distributions, and then retaining the simulations closest to the observations. Although well-established in some fields, whether ABC will work with ecological IBMs is still uncertain. Rejection-ABC was applied to an existing 14-parameter earthworm energy budget IBM for which the available data consist of body mass growth and cocoon production in four experiments. ABC was able to narrow the posterior distributions of seven parameters, estimating credible intervals for each. ABC’s accepted values produced slightly better fits than literature values do. The accuracy of the analysis was assessed using cross-validation and coverage, currently the best available tests. Of the seven unnarrowed parameters, ABC revealed that three were correlated with other parameters, while the remaining four were found to be not estimable given the data available. It is often desirable to compare models to see whether all component modules are necessary. Here we used ABC model selection to compare the full model with a simplified version which removed the earthworm’s movement and much of the energy budget. We are able to show that inclusion of the energy budget is necessary for a good fit to the data. We show how our methodology can inform future modelling cycles, and briefly discuss how more advanced versions of ABC may be applicable to IBMs. We conclude that ABC has the potential to represent uncertainty in model structure, parameters and predictions, and to embed the often complex process of optimizing an IBM’s structure and parameters within an established statistical framework, thereby making the process more transparent and objective.

Toward multifactorial null models of range contraction in terrestrial vertebrates

The contraction of a species’ distribution range, which results from the extirpation of local populations, generally precedes its extinction. Therefore, understanding drivers of range contraction is important for conservation and management. Although there are many processes that can potentially lead to local extirpation and range contraction, three main null models have been proposed: demographic, contagion, and refuge. The first two models postulate that the probability of local extirpation for a given area depends on its relative position within the range; but these models generate distinct spatial predictions because they assume either a ubiquitous (demographic) or a clinal (contagion) distribution of threats. The third model (refuge) postulates that extirpations are determined by the intensity of human impacts, leading to heterogeneous spatial predictions potentially compatible with those made by the other two null models. A few previous studies have explored the generality of some of these null models, but we present here the first comprehensive evaluation of all three models. Using descriptive indices and regression analyses we contrast the predictions made by each of the null models using empirical spatial data describing range contraction in 386 terrestrial vertebrates (mammals, birds, amphibians, and reptiles) distributed across the World. Observed contraction patterns do not consistently conform to the predictions of any of the three models, suggesting that these may not be adequate null models to evaluate range contraction dynamics among terrestrial vertebrates. Instead, our results support alternative null models that account for both relative position and intensity of human impacts. These new models provide a better multifactorial baseline to describe range contraction patterns in vertebrates. This general baseline can be used to explore how additional factors influence contraction, and ultimately extinction for particular areas or species as well as to predict future changes in light of current and new threats.

Datasets related to in-land water for limnology and remote sensing applications: distance-to-land, distance-to-water, water-body identifier and lake-centre co-ordinates

Datasets containing information to locate and identify water bodies have been generated from data locating static-water-bodies with resolution of about 300 m (1/360 deg) recently released by the Land Cover Climate Change Initiative (LC CCI) of the European Space Agency. The LC CCI water-bodies dataset has been obtained from multi-temporal metrics based on time series of the backscattered intensity recorded by ASAR on Envisat between 2005 and 2010. The new derived datasets provide coherently: distance to land, distance to water, water-body identifiers and lake-centre locations. The water-body identifier dataset locates the water bodies assigning the identifiers of the Global Lakes and Wetlands Database (GLWD), and lake centres are defined for in-land waters for which GLWD IDs were determined. The new datasets therefore link recent lake/reservoir/wetlands extent to the GLWD, together with a set of coordinates which locates unambiguously the water bodies in the database. Information on distance-to-land for each water cell and the distance-to-water for each land cell has many potential applications in remote sensing, where the applicability of geophysical retrieval algorithms may be affected by the presence of water or land within a satellite field of view (image pixel). During the generation and validation of the datasets some limitations of the GLWD database and of the LC CCI water-bodies mask have been found. Some examples of the inaccuracies/limitations are presented and discussed. Temporal change in water-body extent is common. Future versions of the LC CCI dataset are planned to represent temporal variation, and this will permit these derived datasets to be updated.

Can climate models represent the precipitation associated with extratropical cyclones?

Extratropical cyclones produce the majority of precipitation in many regions of the extratropics. This study evaluates the ability of a climate model, HiGEM, to reproduce the precipitation associated with extratropical cyclones. The model is evaluated using the ERA-Interim reanalysis and GPCP dataset. The analysis employs a cyclone centred compositing technique, evaluates composites across a range of geographical areas and cyclone intensities and also investigates the ability of the model to reproduce the climatological distribution of cyclone associated precipitation across the Northern Hemisphere. Using this phenomena centred approach provides an ability to identify the processes which are responsible for climatological biases in the model. Composite precipitation intensities are found to be comparable when all cyclones across the Northern Hemisphere are included. When the cyclones are filtered by region or intensity, differences are found, in particular, HiGEM produces too much precipitation in its most intense cyclones relative to ERA-Interim and GPCP. Biases in the climatological distribution of cyclone associated precipitation are also found, with biases around the storm track regions associated with both the number of cyclones in HiGEM and also their average precipitation intensity. These results have implications for the reliability of future projections of extratropical precipitation from the model.

A cautionary note on the use of Ornstein Uhlenbeck models in macroevolutionary studies

Phylogenetic comparative methods are increasingly used to give new insights into the dynamics of trait evolution in deep time. For continuous traits the core of these methods is a suite of models that attempt to capture evolutionary patterns by extending the Brownian constant variance model. However, the properties of these models are often poorly understood, which can lead to the misinterpretation of results. Here we focus on one of these models – the Ornstein Uhlenbeck (OU) model. We show that the OU model is frequently incorrectly favoured over simpler models when using Likelihood ratio tests, and that many studies fitting this model use datasets that are small and prone to this problem. We also show that very small amounts of error in datasets can have profound effects on the inferences derived from OU models. Our results suggest that simulating fitted models and comparing with empirical results is critical when fitting OU and other extensions of the Brownian model. We conclude by making recommendations for best practice in fitting OU models in phylogenetic comparative analyses, and for interpreting the parameters of the OU model.

Fast and slow responses of Southern Ocean sea surface temperature to SAM in coupled climate models

We investigate how sea surface temperatures (SSTs) around Antarctica respond to the Southern An- nular Mode (SAM) on multiple timescales. To that end we examine the relationship between SAM and SST within unperturbed preindustrial control simulations of coupled general circulation models (GCMs) included in the Climate Modeling Intercomparison Project phase 5 (CMIP5). We develop a technique to extract the re- sponse of the Southern Ocean SST (55◦S−70◦S) to a hypothetical step increase in the SAM index. We demonstrate that in many GCMs, the expected SST step re- sponse function is nonmonotonic in time. Following a shift to a positive SAM anomaly, an initial cooling regime can transition into surface warming around Antarctica. However, there are large differences across the CMIP5 ensemble. In some models the step response function never changes sign and cooling persists, while in other GCMs the SST anomaly crosses over from negative to positive values only three years after a step increase in the SAM. This intermodel diversity can be related to differences in the models’ climatological thermal ocean stratification in the region of seasonal sea ice around Antarctica. Exploiting this relationship, we use obser- vational data for the time-mean meridional and vertical temperature gradients to constrain the real Southern Ocean response to SAM on fast and slow timescales.

CHEMICAL EVOLUTION OF THE GALACTIC BULGE: SINGLE AND DOUBLE INFALL MODELS

Recent work, has produced a wealth of data concerning the chemical evolution of the Galactic bulge, both for stars and nebulae. Present theoretical models generally adopt it limited range of such constraints, frequenfly using it single chemical element (usually iron), which is not enough to describe it unambiguously. In this work, we take into account contraints involving,9 Many chemical elements as possible, basically obtained from bulge nebulae and stars. Our main goal is to show that different scenarios can describe, at least partially the abundance distribution and several dishuice-independent correlations for these objects . Three classes of models were developed. The first is it one-zone, single-infall model, the. Second is it one-zone, double-infall model and the third is a multizone, double-infall model. We show that a one-zone model with it single infall episode is able to reproduce some of the observational data, but the best results tire achieved using it multizone, double-infall model.

Chirikov diffusion in the asteroidal three-body resonance (5,-2,-2)

The theory of diffusion in many-dimensional Hamiltonian system is applied to asteroidal dynamics. The general formulation developed by Chirikov is applied to the NesvornA1/2-Morbidelli analytic model of three-body (three-orbit) mean-motion resonances (Jupiter-Saturn-asteroid). In particular, we investigate the diffusion along and across the separatrices of the (5, -2, -2) resonance of the (490) Veritas asteroidal family and their relationship to diffusion in semi-major axis and eccentricity. The estimations of diffusion were obtained using the Melnikov integral, a Hadjidemetriou-type sympletic map and numerical integrations for times up to 10(8) years.

Seed Cell Wall Storage Polysaccharides: Models to Understand Cell Wall Biosynthesis and Degradation

Cell wall storage polysaccharides (CWSPs) are found as the principal storage compounds in seeds of many taxonomically important groups of plants. These groups developed extremely efficient biochemical mechanisms to disassemble cell walls and use the products of hydrolysis for growth. To accumulate these storage polymers, developing seeds also contain relatively high activities of noncellulosic polysaccharide synthases and thus are interesting models to seek the discovery of genes and enzymes related to polysaccharide biosynthesis. CWSP systems offer opportunities to understand phenomena ranging from polysaccharide deposition during seed maturation to the control of source-sink relationship in developing seedlings. By studying polysaccharide biosynthesis and degradation and the consequences for cell and physiological behavior, we can use these models to develop future biotechnological applications.

Mitochondrial sequence data and Dipsacales phylogeny: Mixed models, partitioned Bayesian analyses, and model selection

Phylogenetic analyses of chloroplast DNA sequences, morphology, and combined data have provided consistent support for many of the major branches within the angiosperm, clade Dipsacales. Here we use sequences from three mitochondrial loci to test the existing broad scale phylogeny and in an attempt to resolve several relationships that have remained uncertain. Parsimony, maximum likelihood, and Bayesian analyses of a combined mitochondrial data set recover trees broadly consistent with previous studies, although resolution and support are lower than in the largest chloroplast analyses. Combining chloroplast and mitochondrial data results in a generally well-resolved and very strongly supported topology but the previously recognized problem areas remain. To investigate why these relationships have been difficult to resolve we conducted a series of experiments using different data partitions and heterogeneous substitution models. Usually more complex modeling schemes are favored regardless of the partitions recognized but model choice had little effect on topology or support values. In contrast there are consistent but weakly supported differences in the topologies recovered from coding and non-coding matrices. These conflicts directly correspond to relationships that were poorly resolved in analyses of the full combined chloroplast-mitochondrial data set. We suggest incongruent signal has contributed to our inability to confidently resolve these problem areas. (c) 2007 Elsevier Inc. All rights reserved.

A hands-on approach for fitting long-term survival models under the GAMLSS framework

In many data sets from clinical studies there are patients insusceptible to the occurrence of the event of interest. Survival models which ignore this fact are generally inadequate. The main goal of this paper is to describe an application of the generalized additive models for location, scale, and shape (GAMLSS) framework to the fitting of long-term survival models. in this work the number of competing causes of the event of interest follows the negative binomial distribution. In this way, some well known models found in the literature are characterized as particular cases of our proposal. The model is conveniently parameterized in terms of the cured fraction, which is then linked to covariates. We explore the use of the gamlss package in R as a powerful tool for inference in long-term survival models. The procedure is illustrated with a numerical example. (C) 2009 Elsevier Ireland Ltd. All rights reserved.

Chitosan as a removing agent of beta-lactoglobulin from membrane models

Many chitosan biological activities depend on the interaction with biomembranes, but so far it has not been possible to obtain molecular-level evidence of chitosan action. In this article, we employ Langmuir phospholipid monolayers as cell membrane models and show that chitosan is able to remove beta-lactoglobulin (BLG) from negatively charged dimyristoyl phosphatidic acid (DMPA) and dipalmitoyl phosphatidyl glycerol (DPPG). This was shown with surface pressure isotherms and elasticity and PM-IRRAS measurements in the Langmuir monolayers, in addition to quartz crystal microbalance and fluorescence spectroscopy measurements for Langmuir-Blodgett (LB) films transferred onto solid substrates. Some specificity was noted in the removal action because chitosan was unable to remove BLG incorporated into neutral dipalmitoyl phosphatidyl choline (DPPC) and cholesterol monolayers and had no effect on horseradish peroxidase and urease interacting with DMPA. An obvious biological implication of these findings is to offer reasons that chitosan can remove BLG from lipophilic environments, as reported in the recent literature.