On Systems Thinking, Systems Biology, and the in Silico Plant

The recent summary report of a Department of Energy Workshop on Plant Systems Biology (P.V. Minorsky [2003] Plant Physiol 132: 404-409) offered a welcomed advocacy for systems analysis as essential in understanding plant development, growth, and production. The goal of the Workshop was to consider methods for relating the results of molecular research to real-world challenges in plant production for increased food supplies, alternative energy sources, and environmental improvement. The rather surprising feature of this report, however, was that the Workshop largely overlooked the rich history of plant systems analysis extending over nearly 40 years (Sinclair and Seligman, 1996) that has considered exactly those challenges targeted by the Workshop. Past systems research has explored and incorporated biochemical and physiological knowledge into plant simulation models from a number of perspectives. The research has resulted in considerable understanding and insight about how to simulate plant systems and the relative contribution of various factors in influencing plant production. These past activities have contributed directly to research focused on solving the problems of increasing biomass production and crop yields. These modeling approaches are also now providing an avenue to enhance integration of molecular genetic technologies in plant improvement (Hammer et al., 2002).

Modelling Leaf Production and Crop Development in Maize (Zea mays L.) After Tassel Initiation Under Diverse Conditions of Temperature and Photoperiod

Prediction of the initiation, appearance and emergence of leaves is critically important to the success of simulation models of crop canopy development and some aspects of crop ontogeny. Data on leaf number and crop ontogeny were collected on five cultivars of maize differing widely in maturity and genetic background grown under natural and extended photoperiods, and planted on seven sowing dates from October 1993 to March 1994 at Gatton, South-east Queensland. The same temperature coefficients were established for crop ontogeny before silking, and the rates of leaf initiation, leaf tip appearance and full leaf expansion, the base, optimum and maximum temperatures for each being 8, 34 and 40 degrees C. After silking, the base temperature for ontogeny was 0 degrees C, but the optimum and maximum temperatures remained unchanged. The rates of leaf initiation, appearance of leaf tips and full leaf expansion varied in a relatively narrow range across sowing times and photoperiod treatments, with average values of 0.040 leaves (degrees Cd)-1, 0.021 leaves (degrees Cd)-1, and 0.019 leaves (degrees Cd)-1, respectively. The relationships developed in this study provided satisfactory predictions of leaf number and crop ontogeny (tassel initiation to silking, emergence to silking and silking to physiological maturity) when assessed using independent data from Gatton (South eastern Queensland), Katherine and Douglas Daly (Northern Territory), Walkamin (North Queensland) and Kununurra (Western Australia).

Modeling chickpea growth and development: Nitrogen accumulation and use

Quantitative information regarding nitrogen (N) accumulation and its distribution to leaves, stems and grains under varying environmental and growth conditions are limited for chickpea (Cicer arietinum L.). The information is required for the development of crop growth models and also for assessment of the contribution of chickpea to N balances in cropping systems. Accordingly, these processes were quantified in chickpea under different environmental and growth conditions (still without water or N deficit) using four field experiments and 1325 N measurements. N concentration ([N]) in green leaves was 50 mg g-1 up to beginning of seed growth, and then it declined linearly to 30 mg g-1 at the end of seed growth phase. [N] in senesced leaves was 12 mg g-1. Stem [N] decreased from 30 mg g-1 early in the season to 8 mg g-1 in senesced stems at maturity. Pod [N] was constant (35 mg g-1), but grain [N] decreased from 60 mg g-1 early in seed growth to 43 mg g-1 at maturity. Total N accumulation ranged between 9 and 30 g m-2. N accumulation was closely linked to biomass accumulation until maturity. N accumulation efficiency (N accumulation relative to biomass accumulation) was 0.033 g g-1 where total biomass was -2 and during early growth period, but it decreased to 0.0176 g g-1 during the later growth period when total biomass was >218 g m-2. During vegetative growth (up to first-pod), 58% of N was partitioned to leaves and 42% to stems. Depending on growth conditions, 37-72% of leaf N and 12-56% of stem N was remobilized to the grains. The parameter estimates and functions obtained in this study can be used in chickpea simulation models to simulate N accumulation and distribution.

A novel breeding programme for improved growth in barramundi Lates calcarifer (Bloch) using foundation stock from progeny-tested parents

Rapid genetic gains for growth in barramundi ( Lates calcarifer) appear achievable by starting a breeding programme using foundation stock from progeny tested broodstock. The potential gains of this novel breeding design were investigated using biologically feasible scenarios tested with computer simulation models. The design involves the production of a large number of full-sib families using artificial mating which are compared in common growout conditions. The estimated breeding values of their paternal parents are calculated using a binomial probit analysis to assess their suitability as foundation broodstock. The programme can theoretically yield faster rates of genetic gain compared to other breeding programmes for aquaculture species. Assuming a heritability of 0.25 for growth, foundation broodstock evaluated in two years had breeding values for faster growth ranging from 21% to 51% depending on the genetic diversity of stock under evaluation. As a comparison it will take between nine and twenty-two years to identify broodstock with similar breeding values in a contemporary barramundi breeding programme.

Prediction of mortality and conception rates of beef breeding cattle in northern Australia

In current simulation packages for the management of extensive beef-cattle enterprises, the relationships for the key biological rates (namely conception and mortality) are quite rudimentary. To better estimate these relationships, cohort-level data covering 17 100 cow-years from six sites across northern Australia were collated and analysed. Further validation data, from 7200 cow-years, were then used to test these relationships. Analytical problems included incomplete and non-standardised data, considerable levels of correlation among the 'independent' variables, and the close similarity of alternate possible models. In addition to formal statistical analyses of these data, the theoretical equations for predicting mortality and conception rates in the current simulation models were reviewed, and then reparameterised and recalibrated where appropriate. The final models explained up to 80% of the variation in the data. These are now proposed as more accurate and useful models to be used in the prediction of biological rates in simulation studies for northern Australia. © The State of Queensland (through the Department of Agriculture, Fisheries and Forestry) 2012. © CSIRO.

Variation in water extraction with maize plant density and its impact on model application

The use of maize simulation models to determine the optimum plant population for rainfed environments allows the evaluation of plant populations over multiple years and locations at a lower cost than traditional field experimentation. However the APSIM maize model that has been used to conduct some of these 'virtual' experiments assumes that the maximum rate of soil water extraction by the crop root system is constant across plant populations. This untested assumption may cause grain yield to be overestimated in lower plant populations. A field experiment was conducted to determine whether maximum rates of water extraction vary with plant population, and the maximum rate of soil water extraction was estimated for three plant populations (2.4, 3.5 and 5.5 plants m(-2)) under water limited conditions. Maximum soil water extraction rates in the field experiment decreased linearly with plant population, and no difference was detected between plant populations for the crop lower limit of soil water extraction. Re-analysis of previous maize simulation experiments demonstrated that the use of inappropriately high extraction-rate parameters at low plant populations inflated predictions of grain yield, and could cause erroneous recommendations to be made for plant population. The results demonstrate the importance of validating crop simulation models across the range of intended treatments. (C) 2013 Elsevier E.V. All rights reserved.

Quarantine and surveillance strategies for plant pathogen detection and control

Incursions of plant pests and diseases pose serious threats to food security, agricultural productivity and the natural environment. One of the challenges in confidently delimiting and eradicating incursions is how to choose from an arsenal of surveillance and quarantine approaches in order to best control multiple dispersal pathways. Anthropogenic spread (propagules carried on humans or transported on produce or equipment) can be controlled with quarantine measures, which in turn can vary in intensity. In contrast, environmental spread processes are more difficult to control, but often have a temporal signal (e.g. seasonality) which can introduce both challenges and opportunities for surveillance and control. This leads to complex decisions regarding when, where and how to search. Recent modelling investigations of surveillance performance have optimised the output of simulation models, and found that a risk-weighted randomised search can perform close to optimally. However, exactly how quarantine and surveillance strategies should change to reflect different dispersal modes remains largely unaddressed. Here we develop a spatial simulation model of a plant fungal-pathogen incursion into an agricultural region, and its subsequent surveillance and control. We include structural differences in dispersal via the interplay of biological, environmental and anthropogenic connectivity between host sites (farms). Our objective was to gain broad insights into the relative roles played by different spread modes in propagating an invasion, and how incorporating knowledge of these spread risks may improve approaches to quarantine restrictions and surveillance. We find that broad heuristic rules for quarantine restrictions fail to contain the pathogen due to residual connectivity between sites, but surveillance measures enable early detection and successfully lead to suppression of the pathogen in all farms. Alternative surveillance strategies attain similar levels of performance by incorporating environmental or anthropogenic dispersal risk in the prioritisation of sites. Our model provides the basis to develop essential insights into the effectiveness of different surveillance and quarantine decisions for fungal pathogen control. Parameterised for authentic settings it will aid our understanding of how the extent and resolution of interventions should suitably reflect the spatial structure of dispersal processes.

Ecological communities in variable environments : dynamics and diversity under coloured environmental stochasticity

While environmental variation is an ubiquitous phenomenon in the natural world which has for long been appreciated by the scientific community recent changes in global climatic conditions have begun to raise consciousness about the economical, political and sociological ramifications of global climate change. Climate warming has already resulted in documented changes in ecosystem functioning, with direct repercussions on ecosystem services. While predicting the influence of ecosystem changes on vital ecosystem services can be extremely difficult, knowledge of the organisation of ecological interactions within natural communities can help us better understand climate driven changes in ecosystems. The role of environmental variation as an agent mediating population extinctions is likely to become increasingly important in the future. In previous studies population extinction risk in stochastic environmental conditions has been tied to an interaction between population density dependence and the temporal autocorrelation of environmental fluctuations. When populations interact with each other, forming ecological communities, the response of such species assemblages to environmental stochasticity can depend, e.g., on trophic structure in the food web and the similarity in species-specific responses to environmental conditions. The results presented in this thesis indicate that variation in the correlation structure between species-specific environmental responses (environmental correlation) can have important qualitative and quantitative effects on community persistence and biomass stability in autocorrelated (coloured) environments. In addition, reddened environmental stochasticity and ecological drift processes (such as demographic stochasticity and dispersal limitation) have important implications for patterns in species relative abundances and community dynamics over time and space. Our understanding of patterns in biodiversity at local and global scale can be enhanced by considering the relevance of different drift processes for community organisation and dynamics. Although the results laid out in this thesis are based on mathematical simulation models, they can be valuable in planning effective empirical studies as well as in interpreting existing empirical results. Most of the metrics considered here are directly applicable to empirical data.

Simulating and Forecasting the Demand for New Consumer Durables

A diffusion/replacement model for new consumer durables designed to be used as a long-term forecasting tool is developed. The model simulates new demand as well as replacement demand over time. The model is called DEMSIM and is built upon a counteractive adoption model specifying the basic forces affecting the adoption behaviour of individual consumers. These forces are the promoting forces and the resisting forces. The promoting forces are further divided into internal and external influences. These influences are operationalized within a multi-segmental diffusion model generating the adoption behaviour of the consumers in each segment as an expected value. This diffusion model is combined with a replacement model built upon the same segmental structure as the diffusion model. This model generates, in turn, the expected replacement behaviour in each segment. To be able to use DEMSIM as a forecasting tool in early stages of a diffusion process estimates of the model parameters are needed as soon as possible after product launch. However, traditional statistical techniques are not very helpful in estimating such parameters in early stages of a diffusion process. To enable early parameter calibration an optimization algorithm is developed by which the main parameters of the diffusion model can be estimated on the basis of very few sales observations. The optimization is carried out in iterative simulation runs. Empirical validations using the optimization algorithm reveal that the diffusion model performs well in early long-term sales forecasts, especially as it comes to the timing of future sales peaks.

Analysis of forest soil chemistry and hydrology with a dynamic model ACIDIC.

In this study we analyze how the ion concentrations in forest soil solution are determined by hydrological and biogeochemical processes. A dynamic model ACIDIC was developed, including processes common to dynamic soil acidification models. The model treats up to eight interacting layers and simulates soil hydrology, transpiration, root water and nutrient uptake, cation exchange, dissolution and reactions of Al hydroxides in solution, and the formation of carbonic acid and its dissociation products. It includes also a possibility to a simultaneous use of preferential and matrix flow paths, enabling the throughfall water to enter the deeper soil layers in macropores without first reacting with the upper layers. Three different combinations of routing the throughfall water via macro- and micropores through the soil profile is presented. The large vertical gradient in the observed total charge was simulated succesfully. According to the simulations, gradient is mostly caused by differences in the intensity of water uptake, sulfate adsorption and organic anion retention at the various depths. The temporal variations in Ca and Mg concentrations were simulated fairly well in all soil layers. For H+, Al and K there were much more variation in the observed than in the simulated concentrations. Flow in macropores is a possible explanation for the apparent disequilibrium of the cation exchange for H+ and K, as the solution H+ and K concentrations have great vertical gradients in soil. The amount of exchangeable H+ increased in the O and E horizons and decreased in the Bs1 and Bs2 horizons, the net change in whole soil profile being a decrease. A large part of the decrease of the exchangeable H+ in the illuvial B horizon was caused by sulfate adsorption. The model produces soil water amounts and solution ion concentrations which are comparable to the measured values, and it can be used in both hydrological and chemical studies of soils.

Queuing strategies for qos management of traffic between lans interconnected by low-speed wan links

We provide a comparative performance evaluation of packet queuing and link admission strategies for low-speed wide area network Links (e.g. 9600 bps, 64 kbps) that interconnect relatively highspeed, connectionless local area networks (e.g. 10 Mbps). In particular, we are concerned with the problem of providing differential quality of service to interLAN remote terminal and file transfer sessions, and throughput fairness between interLAN file transfer sessions. We use analytical and simulation models to study a variety of strategies. Our work also serves to address the performance comparison of connectionless vs. connection-oriented interconnection of CLNS LANS. When provision of priority at the physical transmission level is not feasible, we show, for low-speed WAN links (e.g. 9600 bps), the superiority of connection-oriented interconnection of connectionless LANs, with segregation of traffic streams with different QoS requirements into different window flow controlled connections. Such an implementation can easily be obtained by transporting IP packets over an X.25 WAN. For 64 kbps WAN links, there is a drop in file transfer throughputs, owing to connection overheads, but the other advantages are retained, The same solution also helps to provide throughput fairness between interLAN file transfer sessions. We also provide a corroboration of some of our modelling results with results from an experimental test-bed.

Impacts of modern transitions on thermal comfort in vernacular dwellings in warm-humid climate of Sugganahalli (India)

Sugganahalli, a rural vernacular community in a warm-humid region in South India, is under transition towards adopting modern construction practices. Vernacular local building elements like rubble walls and mud roofs are given way to burnt brick walls and reinforced cement concrete (RCC)/tin roofs. Over 60% of Indian population is rural, and implications of such transitions on thermal comfort and energy in buildings are crucial to understand. Vernacular architecture evolves adopting local resources in response to the local climate adopting passive solar designs. This paper investigates the effectiveness of passive solar elements on the indoor thermal comfort by adopting modern climate-responsive design strategies. Dynamic simulation models validated by measured data have also been adopted to determine the impact of the transition from vernacular to modern material-configurations. Age-old traditional design considerations were found to concur with modern understanding into bio-climatic response and climate-responsiveness. Modern transitions were found to increase the average indoor temperatures in excess of 7 degrees C. Such transformations tend to shift the indoor conditions to a psychrometric zone that is likely to require active air-conditioning. Also, the surveyed thermal sensation votes were found to lie outside the extended thermal comfort boundary for hot developing countries provided by Givoni in the bio-climatic chart.

In-cylinder investigations and analysis of a SI gas engine fuelled with H-2 and CO rich syngas fuel: Sensitivity analysis of combustion descriptors for engine diagnostics and control

The sensitivity of combustion phasing and combustion descriptors to ignition timing, load and mixture quality on fuelling a multi-cylinder natural gas engine with bio-derived H-2 and CO rich syngas is addressed. While the descriptors for conventional fuels are well established and are in use for closed loop engine control, presence of H-2 in syngas potentially alters the mixture properties and hence combustion phasing, necessitating the current study. The ability of the descriptors to predict abnormal combustion, hitherto missing in the literature, is also addressed. Results from experiments using multi-cylinder engines and numerical studies using zero dimensional Wiebe function based simulation models are reported. For syngas with 20% H-2 and CO and 2% CH4 (producer gas), an ignition retard of 5 +/- 1 degrees was required compared to natural gas ignition timing to achieve peak load of 72.8 kWe. It is found that, for syngas, whose flammability limits are 0.42-1.93, the optimal engine operation was at an equivalence ratio of 1.12. The same methodology is extended to a two cylinder engine towards addressing the influence of syngas composition, especially H-2 fraction (varying from 13% to 37%), on the combustion phasing. The study confirms the utility of pressure trace derived combustion descriptors, except for the pressure trace first derivative, in describing the MBT operating condition of the engine when fuelled with an alternative fuel. Both experiments and analysis suggest most of the combustion descriptors to be independent of the engine load and mixture quality. A near linear relationship with ignition angle is observed. The general trend(s) of the combustion descriptors for syngas fuelled operation are similar to those of conventional fuels; the differences in sensitivity of the descriptors for syngas fuelled engine operation requires re-calibration of control logic for MBT conditions. Copyright (C) 2014, Hydrogen Energy Publications, LLC. Published by Elsevier Ltd. All rights reserved.

Implication of scaling hierarchy associated with nonequilibrium: field and particulate

The advent of nanotechnology has necessitated a better understanding of how material microstructure changes at the atomic level would affect the macroscopic properties that control the performance. Such a challenge has uncovered many phenomena that were not previously understood and taken for granted. Among them are the basic foundation of dislocation theories which are now known to be inadequate. Simplifying assumptions invoked at the macroscale may not be applicable at the micro- and/or nanoscale. There are implications of scaling hierrachy associated with in-homegeneity and nonequilibrium. of physical systems. What is taken to be homogeneous and equilibrium at the macroscale may not be so when the physical size of the material is reduced to microns. These fundamental issues cannot be dispensed at will for the sake of convenience because they could alter the outcome of predictions. Even more unsatisfying is the lack of consistency in modeling physical systems. This could translate to the inability for identifying the relevant manufacturing parameters and rendering the end product unpractical because of high cost. Advanced composite and ceramic materials are cases in point. Discussed are potential pitfalls for applying models at both the atomic and continuum levels. No encouragement is made to unravel the truth of nature. Let it be partiuclates, a smooth continuum or a combination of both. The present trend of development in scaling tends to seek for different characteristic lengths of material microstructures with or without the influence of time effects. Much will be learned from atomistic simulation models to show how results could differ as boundary conditions and scales are changed. Quantum mechanics, continuum and cosmological models provide evidence that no general approach is in sight. Of immediate interest is perhaps the establishment of greater precision in terminology so as to better communicate results involving multiscale physical events.

Modelling the dynamics of phytoplankton with the needs of the end user in mind

How to regulate phytoplankton growth in water supply reservoirs has continued to occupy managers and strategists for some fifty years or so, now, and mathematical models have always featured in their design and operational constraints. In recent years, rather more sophisticated simulation models have begun to be available and these, ideally, purport to provide the manager with improved forecasting of plankton blooms, the likely species and the sort of decision support that might permit management choices to be selected with increased confidence. This account describes the adaptation and application of one such model, PROTECH (Phytoplankton RespOnses To Environmental CHange) to the problems of plankton growth in reservoirs. This article supposes no background knowledge of the main algal types; neither does it attempt to catalogue the problems that their abundance may cause in lakes and reservoirs.