Modelling sediment delivery ratio over the Murray Darling Basin

This paper presents a scientific and technical description of the modelling framework and the main results of modelling the long-term average sediment delivery at hillslope to medium-scale catchments over the entire Murray Darling Basin (MDB). A theoretical development that relates long-term averaged sediment delivery to the statistics of rainfall and catchment parameters is presented. The derived flood frequency approach was adapted to investigate the problem of regionalization of the sediment delivery ratio (SDR) across the Basin. SDR, a measure of catchment response to the upland erosion rate, was modeled by two lumped linear stores arranged in series: hillslope transport to the nearest streams and flow routing in the channel network. The theory shows that the ratio of catchment sediment residence time (SRT) to average effective rainfall duration is the most important control in the sediment delivery processes. In this study, catchment SRTs were estimated using travel time for overland flow multiplied by an enlargement factor which is a function of particle size. Rainfall intensity and effective duration statistics were regionalized by using long-term measurements from 195 pluviograph sites within and around the Basin. Finally, the model was implemented across the MDB by using spatially distributed soil, vegetation, topographical and land use properties under Geographic Information System (GIs) environment. The results predict strong variations in SDR from close to 0 in floodplains to 70% in the eastern uplands of the Basin. (c) 2005 Elsevier Ltd. All rights reserved.

Fighting fit: thermal plasticity of metabolic function and fighting success in the crayfish Cherax destructor

1. We examined the effect of thermal acclimation on fighting success and underlying performance traits in the crayfish Cherax destructor. We tested the hypothesis that animals will be more successful when fighting at their acclimation temperature than at a colder or warmer temperature, and that changes in metabolic capacity underlie differences in behavioural performance. 2. Thermal acclimation (to 20 degrees C and to 30 degrees C) had a significant effect on behavioural contests, and the likelihood of winning was significantly greater when individuals fought at their acclimation temperature against an individual from an alternate acclimation temperature. 3. The ratio of ADP stimulated respiration to proton leak (respiratory control ratio) of isolated mitochondria increased significantly in chelae muscle of the cold-acclimated group, and differences in respiratory control ratio between winners and losers were significantly correlated with the outcome of agonistic encounters. However, acclimation did not affect tall muscle mitochondria or the activity of pyruvate kinase in either chelae or tail muscle. 4. The force produced by closing chelae was thermally insensitive within acclimation groups, and there were no significant differences between acclimation treatments. None the less, differences in chelae width between contestants were significantly correlated with the outcome of agonistic encounters, but this perceived resource holding power did not reflect the actual power of force production. 5. Thermal acclimation in C destructor has beneficial consequences for dominance and competitive ability, and the success of cold acclimated animals at the cold temperatures can be at least partly explained by concomitant up-regulation of oxidative ATP production capacity.

Modelling the conversion of Colombian lowland ecosystems since 1940: Drivers, patterns and rates

In biologically mega-diverse countries that are undergoing rapid human landscape transformation, it is important to understand and model the patterns of land cover change. This problem is particularly acute in Colombia, where lowland forests are being rapidly cleared for cropping and ranching. We apply a conceptual model with a nested set of a priori predictions to analyse the spatial and temporal patterns of land cover change for six 50-100 km(2) case study areas in lowland ecosystems of Colombia. Our analysis included soil fertility, a cost-distance function, and neighbourhood of forest and secondary vegetation cover as independent variables. Deforestation and forest regrowth are tested using logistic regression analysis and an information criterion approach to rank the models and predictor variables. The results show that: (a) overall the process of deforestation is better predicted by the full model containing all variables, while for regrowth the model containing only the auto-correlated neighbourhood terms is a better predictor; (b) overall consistent patterns emerge, although there are variations across regions and time; and (c) during the transformation process, both the order of importance and significance of the drivers change. Forest cover follows a consistent logistic decline pattern across regions, with introduced pastures being the major replacement land cover type. Forest stabilizes at 2-10% of the original cover, with an average patch size of 15.4 (+/- 9.2) ha. We discuss the implications of the observed patterns and rates of land cover change for conservation planning in countries with high rates of deforestation. (c) 2005 Elsevier Ltd. All rights reserved.

On the modelling of wave breaking and set-up on coral reefs

An extended refraction-diffraction equation [Massel, S.R., 1993. Extended refraction-diffraction equation for surface waves. Coastal Eng. 19, 97-126] has been applied to predict wave transformation and breaking as well as wave-induced set-up on two-dimensional reef profiles of various shapes. A free empirical coefficient alpha in a formula for the average rate of energy dissipation [epsilon(b)] = (alpha rho g omega/8 pi)(root gh/C)(H-3/h) in the modified periodic bore model was found to be a function of the dimensionless parameter F-c0 = (g(1.25)H(0)(0.5)T(2.5))/h(r)(1.75), proposed by Gourlay [Gourlayl M.R., 1994. Wave transformation on a coral reef. Coastal Eng. 23, 17-42]. The applicability of the developed model has been demonstrated for reefs of various shapes subjected to various incident wave conditions. Assuming proposed relationships of the coefficient alpha and F-c0, the model provides results on wave height attenuation and set-up elevation which compare well with experimental data. (C) 2000 Elsevier Science B.V. All rights reserved.

Trait physiology and crop modelling to link phenotypic complexity to underlying genetic systems

New tools derived from advances in molecular biology have not been widely adopted in plant breeding because of the inability to connect information at gene level to the phenotype in a manner that is useful for selection. We explore whether a crop growth and development modelling framework can link phenotype complexity to underlying genetic systems in a way that strengthens molecular breeding strategies. We use gene-to-phenotype simulation studies on sorghum to consider the value to marker-assisted selection of intrinsically stable QTLs that might be generated by physiological dissection of complex traits. The consequences on grain yield of genetic variation in four key adaptive traits – phenology, osmotic adjustment, transpiration efficiency, and staygreen – were simulated for a diverse set of environments by placing the known extent of genetic variation in the context of the physiological determinants framework of a crop growth and development model. It was assumed that the three to five genes associated with each trait, had two alleles per locus acting in an additive manner. The effects on average simulated yield, generated by differing combinations of positive alleles for the traits incorporated, varied with environment type. The full matrix of simulated phenotypes, which consisted of 547 location-season combinations and 4235 genotypic expression states, was analysed for genetic and environmental effects. The analysis was conducted in stages with gradually increased understanding of gene-to-phenotype relationships, which would arise from physiological dissection and modelling. It was found that environmental characterisation and physiological knowledge helped to explain and unravel gene and environment context dependencies. We simulated a marker-assisted selection (MAS) breeding strategy based on the analyses of gene effects. When marker scores were allocated based on the contribution of gene effects to yield in a single environment, there was a wide divergence in rate of yield gain over all environments with breeding cycle depending on the environment chosen for the QTL analysis. It was suggested that knowledge resulting from trait physiology and modelling would overcome this dependency by identifying stable QTLs. The improved predictive power would increase the utility of the QTLs in MAS. Developing and implementing this gene-to-phenotype capability in crop improvement requires enhanced attention to phenotyping, ecophysiological modelling, and validation studies to test the stability of candidate QTLs.

Multilevel modelling for inference of genetic regulatory networks

Time-course experiments with microarrays are often used to study dynamic biological systems and genetic regulatory networks (GRNs) that model how genes influence each other in cell-level development of organisms. The inference for GRNs provides important insights into the fundamental biological processes such as growth and is useful in disease diagnosis and genomic drug design. Due to the experimental design, multilevel data hierarchies are often present in time-course gene expression data. Most existing methods, however, ignore the dependency of the expression measurements over time and the correlation among gene expression profiles. Such independence assumptions violate regulatory interactions and can result in overlooking certain important subject effects and lead to spurious inference for regulatory networks or mechanisms. In this paper, a multilevel mixed-effects model is adopted to incorporate data hierarchies in the analysis of time-course data, where temporal and subject effects are both assumed to be random. The method starts with the clustering of genes by fitting the mixture model within the multilevel random-effects model framework using the expectation-maximization (EM) algorithm. The network of regulatory interactions is then determined by searching for regulatory control elements (activators and inhibitors) shared by the clusters of co-expressed genes, based on a time-lagged correlation coefficients measurement. The method is applied to two real time-course datasets from the budding yeast (Saccharomyces cerevisiae) genome. It is shown that the proposed method provides clusters of cell-cycle regulated genes that are supported by existing gene function annotations, and hence enables inference on regulatory interactions for the genetic network.

Optimum switching point assessment of parallel braking resistor and serial current limiter for power system stability enhancement

A new control algorithm using parallel braking resistor (BR) and serial fault current limiter (FCL) for power system transient stability enhancement is presented in this paper. The proposed control algorithm can prevent transient instability during first swing by immediately taking away the transient energy gained in faulted period. It can also reduce generator oscillation time and efficiently make system back to the post-fault equilibrium. The algorithm is based on a new system energy function based method to choose optimal switching point. The parallel BR and serial FCL resistor can be switched at the calculated optimal point to get the best control result. This method allows optimum dissipation of the transient energy caused by disturbance so to make system back to equilibrium in minimum time. Case studies are given to verify the efficiency and effectiveness of this new control algorithm.

A mathematical modelling technique for the analysis of the dynamics of a simple continuous EDA

This paper presents some initial attempts to mathematically model the dynamics of a continuous estimation of distribution algorithm (EDA) based on a Gaussian distribution and truncation selection. Case studies are conducted on both unimodal and multimodal problems to highlight the effectiveness of the proposed technique and explore some important properties of the EDA. With some general assumptions, we show that, for ID unimodal problems and with the (mu, lambda) scheme: (1). The behaviour of the EDA is dependent only on the general shape of the test function, rather than its specific form; (2). When initialized far from the global optimum, the EDA has a tendency to converge prematurely; (3). Given a certain selection pressure, there is a unique value for the proposed amplification parameter that could help the EDA achieve desirable performance; for ID multimodal problems: (1). The EDA could get stuck with the (mu, lambda) scheme; (2). The EDA will never get stuck with the (mu, lambda) scheme.

Dartboard arrangements with a concave penalty function

This paper investigates combinatorial arrangements of the dartboard to maximize a penalty function derived from the differences of adjacent sectors. The particular penalty function is constructed by summing the absolute differences of neighbouring sectors raised to a power between zero and one. The arrangement to give the maximum penalty is found

Modelling the dynamics of genetic algorithms using statistical mechanics

A formalism for modelling the dynamics of Genetic Algorithms (GAs) using methods from statistical mechanics, originally due to Prugel-Bennett and Shapiro, is reviewed, generalized and improved upon. This formalism can be used to predict the averaged trajectory of macroscopic statistics describing the GA's population. These macroscopics are chosen to average well between runs, so that fluctuations from mean behaviour can often be neglected. Where necessary, non-trivial terms are determined by assuming maximum entropy with constraints on known macroscopics. Problems of realistic size are described in compact form and finite population effects are included, often proving to be of fundamental importance. The macroscopics used here are cumulants of an appropriate quantity within the population and the mean correlation (Hamming distance) within the population. Including the correlation as an explicit macroscopic provides a significant improvement over the original formulation. The formalism is applied to a number of simple optimization problems in order to determine its predictive power and to gain insight into GA dynamics. Problems which are most amenable to analysis come from the class where alleles within the genotype contribute additively to the phenotype. This class can be treated with some generality, including problems with inhomogeneous contributions from each site, non-linear or noisy fitness measures, simple diploid representations and temporally varying fitness. The results can also be applied to a simple learning problem, generalization in a binary perceptron, and a limit is identified for which the optimal training batch size can be determined for this problem. The theory is compared to averaged results from a real GA in each case, showing excellent agreement if the maximum entropy principle holds. Some situations where this approximation brakes down are identified. In order to fully test the formalism, an attempt is made on the strong sc np-hard problem of storing random patterns in a binary perceptron. Here, the relationship between the genotype and phenotype (training error) is strongly non-linear. Mutation is modelled under the assumption that perceptron configurations are typical of perceptrons with a given training error. Unfortunately, this assumption does not provide a good approximation in general. It is conjectured that perceptron configurations would have to be constrained by other statistics in order to accurately model mutation for this problem. Issues arising from this study are discussed in conclusion and some possible areas of further research are outlined.

Modelling Apoptotic cell clearance within the atherosclerotic plaque

The aged population have an increased susceptibility to infection, therefore function of the innate immune system may be impaired as we age. Macrophages, and their precursors monocytes, play an important role in host defence in the form of phagocytosis, and also link the innate and adaptive immune system via antigen presentation. Classically-activated 'M1' macrophages are pro-inflammatory, which can be induced by encountering pathogenic material or pro-inflammatory mediators. Alternatively activated 'M2' macrophages have a largely reparative role, including clearance of apoptotic bodies and debris from tissues. Despite some innate immune receptors being implicated in the clearance of apoptotic cells, the process has been observed to have a dominant anti-inflammatory phenotype with cytokines such as IL-10 and TGF-ß being implicated. The atherosclerotic plaque contains recruited monocytes and macrophages, and is a highly inflammatory environment despite high levels of apoptosis. At these sites, monocytes differentiate into macrophages and gorge on lipoproteins, resulting in formation of 'foam cells' which then undergo apoptosis, recruiting further monocytes. This project seeks to understand why, given high levels of apoptosis, the plaque is a pro-inflammatory environment. This phenomenon may be the result of the aged environment or an inability of foam cells to elicit an anti-inflammatory effect in response to dying cells. Here we demonstrate that lipoprotein treatment of macrophages in culture results in reduced capacity to clear apoptotic cells. The effect of lipoprotein treatment on apoptotic cell-mediated immune modulation of macrophage function is currently under study.