Branching processes and computational collapse of discretized unimodal mappings

In computer simulations of smooth dynamical systems, the original phase space is replaced by machine arithmetic, which is a finite set. The resulting spatially discretized dynamical systems do not inherit all functional properties of the original systems, such as surjectivity and existence of absolutely continuous invariant measures. This can lead to computational collapse to fixed points or short cycles. The paper studies loss of such properties in spatial discretizations of dynamical systems induced by unimodal mappings of the unit interval. The problem reduces to studying set-valued negative semitrajectories of the discretized system. As the grid is refined, the asymptotic behavior of the cardinality structure of the semitrajectories follows probabilistic laws corresponding to a branching process. The transition probabilities of this process are explicitly calculated. These results are illustrated by the example of the discretized logistic mapping.

Transformation processes, pathways, and possible sources of distinctive polychlorinated dibenzo-p-dioxin signatures in sink environments

In recent years, studies on environmental samples with unusual dibenzo-p-dioxin (PCDD) congener profiles were reported from a range of countries. These profiles, characterized by a dominance of octachlorinated dibenzodioxin (OCDD) and relatively low in dibenzofuran (PCDF) concentrations, could not be attributed to known sources or formation processes. In the present study, the processes that result in these unusual profiles were assessed using the concentrations and isomer signatures of PCDDs from dated estuarine sediment cores in Queensland, Australia. Increases in relative concentrations of lower chlorinated PODS and a relative decrease of OCDD were correlated with time of sediment deposition. Preferred lateral, anaerobic dechlorination of OCDD represents a likely pathway for these changes. In Queensland sediments, these transformations result in a distinct dominance of isomers fully chlorinated in the 1,4,6,9-positions (1,4-patterns), and similar 1,4-patterns were observed in sediments from elsewhere. Consequently, these environmental samples may not reflect the signatures of the original source, and a reevaluation of source inputs was undertaken. Natural formation of PCDDs, which has previously been suggested, is discussed; however, based on the present results and literature comparisons, we propose an alternative scenario. This scenario hypothesizes that an anthropogenic PCDD precursor input (e.g. pentachlorophenol) results in the contamination. These results and hypothesis imply further investigations are warrented into possible anthropogenic sources in areas where natural PCDD formation has been suggested.

Simulation of insect movement with respect to plant architecture and morphogenesis

Developments in computer and three dimensional (3D) digitiser technologies have made it possible to keep track of the broad range of data required to simulate an insect moving around or over the highly heterogeneous habitat of a plant's surface. Properties of plant parts vary within a complex canopy architecture, and insect damage can induce further changes that affect an animal's movements, development and likelihood of survival. Models of plant architectural development based on Lindenmayer systems (L-systems) serve as dynamic platforms for simulation of insect movement, providing ail explicit model of the developing 3D structure of a plant as well as allowing physiological processes associated with plant growth and responses to damage to be described and Simulated. Simple examples of the use of the L-system formalism to model insect movement, operating Lit different spatial scales-from insects foraging on an individual plant to insects flying around plants in a field-are presented. Such models can be used to explore questions about the consequences of changes in environmental architecture and configuration on host finding, exploitation and its population consequences. In effect this model is a 'virtual ecosystem' laboratory to address local as well as landscape-level questions pertinent to plant-insect interactions, taking plant architecture into account. (C) 2002 Elsevier Science B.V. All rights reserved.

Nitrate in soil humic acids revealed by N-14 nuclear magnetic resonance spectroscopy

Ecosystem management such as plant residue retention and prescribed burning can significantly affect soil organic matter (SOM) composition and, thereby, the closely associated carbon (C) and nitrogen (N) cycling processes, which underpin terrestrial ecosystem productivity and sustainability. Humic acid (HA) is an important SOM component and its chemical composition has attracted much attention. Here we report the first application of N-14 nuclear magnetic resonance (NMR) spectroscopy to soil HA study, revealing the surprising existence of nitrate-N and ammonia-N in the HAs. This newly discovered HA nitrate-N, though in a relatively low concentrations, is closely related to soil N availability and responsive to plant residue management regimes in contrasting forest ecosystems. The HA nitrate-N may be a useful and sensitive biochemical indicator of SOM quality in response to different ecosystem management regimes.

Robust model-order reduction of complex biological processes

This paper addresses robust model-order reduction of a high dimensional nonlinear partial differential equation (PDE) model of a complex biological process. Based on a nonlinear, distributed parameter model of the same process which was validated against experimental data of an existing, pilot-scale BNR activated sludge plant, we developed a state-space model with 154 state variables in this work. A general algorithm for robustly reducing the nonlinear PDE model is presented and based on an investigation of five state-of-the-art model-order reduction techniques, we are able to reduce the original model to a model with only 30 states without incurring pronounced modelling errors. The Singular perturbation approximation balanced truncating technique is found to give the lowest modelling errors in low frequency ranges and hence is deemed most suitable for controller design and other real-time applications. (C) 2002 Elsevier Science Ltd. All rights reserved.

Ecological responses to recent climate change

There is now ample evidence of the ecological impacts of recent climate change, from polar terrestrial to tropical marine environments. The responses of both flora and fauna span an array of ecosystems and organizational hierarchies, from the species to the community levels. Despite continued uncertainty as to community and ecosystem trajectories under global change, our review exposes a coherent pattern of ecological change across systems. Although we are only at an early stage in the projected trends of global warming, ecological responses to recent climate change are already clearly visible.

A note on quasi-stationary distributions of birth-death processes and the SIS logistic epidemic

For Markov processes on the positive integers with the origin as an absorbing state, Ferrari, Kesten, Martinez and Picco studied the existence of quasi-stationary and limiting conditional distributions by characterizing quasi-stationary distributions as fixed points of a transformation Phi on the space of probability distributions on {1, 2,.. }. In the case of a birth-death process, the components of Phi(nu) can be written down explicitly for any given distribution nu. Using this explicit representation, we will show that Phi preserves likelihood ratio ordering between distributions. A conjecture of Kryscio and Lefevre concerning the quasi-stationary distribution of the SIS logistic epidemic follows as a corollary.

Ecological criteria for evaluating candidate sites for marine reserves

Several schemes have been developed to help select the locations of marine reserves. All of them combine social, economic, and biological criteria, and few offer any guidance as to how to prioritize among the criteria identified. This can imply that the relative weights given to different criteria are unimportant. Where two sites are of equal value ecologically; then socioeconomic criteria should dominate the choice of which should be protected. However, in many cases, socioeconomic criteria are given equal or greater weight than ecological considerations in the choice of sites. This can lead to selection of reserves with little biological value that fail to meet many of the desired objectives. To avoid such a possibility, we develop a series of criteria that allow preliminary evaluation of candidate sites according to their relative biological values in advance of the application of socioeconomic criteria. We include criteria that,. while not strictly biological, have a strong influence on the species present or ecological processes. Out scheme enables sites to be assessed according to their biodiversity, the processes which underpin that diversity, and the processes that support fisheries and provide a spectrum of other services important to people. Criteria that capture biodiversity values include biogeographic representation, habitat representation and heterogeneity, and presence of species or populations of special interest (e.g., threatened species). Criteria that capture sustainability of biodiversity and fishery values include the size of reserves necessary to protect viable habitats, presence of exploitable species, vulnerable life stages, connectivity among reserves, links among ecosystems, and provision of ecosystem services to people. Criteria measuring human and natural threats enable candidate sites to be eliminated from consideration if risks are too great, but also help prioritize among sites where threats can be mitigated by protection. While our criteria can be applied to the design of reserve networks, they also enable choice of single reserves to be made in the context of the attributes of existing protected areas. The overall goal of our scheme is to promote the development of reserve networks that will maintain biodiversity and ecosystem functioning at large scales. The values of eco-system goods and services for people ultimately depend on meeting this objective.

Using the flexibility index to compare batch and continuous activated sludge processes

This paper considers the question of which is better: the batch or the continuous activated sludge processes? It is an important question because dissension still exists in the wastewater industry as to the relative merits of each of the processes. A review of perceived differences in the processes from the point of view of two related disciplines, process engineering and biotechnology, is presented together with the results of previous comparative studies. These reviews highlight possible areas where more understanding is required. This is provided in the paper by application of the flexibility index to two case studies. The flexibility index is a useful process design tool that measures the ability of the process to cope with long term changes in operation.

Is maximizing resilience compatible with established ecological goal functions?

Cropp and Gabric [Ecosystem adaptation: do ecosystems maximise resilience? Ecology. In press] used a simple phytoplanktonzooplankton-nutrient model and a genetic algorithm to determine the parameter values that would maximize the value of certain goal functions. These goal functions were to maximize biomass, maximize flux, maximize flux to biomass ratio, and maximize resilience. It was found that maximizing goal functions maximized resilience. The objective of this study was to investigate whether the Cropp and Gabric [Ecosystem adaptation: do ecosystems maximise resilience? Ecology. In press] result was indicative of a general ecosystem principle, or peculiar to the model and parameter ranges used. This study successfully replicated the Cropp and Gabric [Ecosystem adaptation: do ecosystems maximise resilience? Ecology. In press] experiment for a number of different model types, however, a different interpretation of the results is made. A new metric, concordance, was devised to describe the agreement between goal functions. It was found that resilience has the highest concordance of all goal functions trialled. for most model types. This implies that resilience offers a compromise between the established ecological goal functions. The parameter value range used is found to affect the parameter versus goal function relationships. Local maxima and minima affected the relationship between parameters and goal functions, and between goal functions. (C) 2003 Elsevier B.V. All rights reserved.

The Lake Tekapo experiment (LTEX): An investigation of atmospheric boundary layer processes in complex terrain

A research program on atmospheric boundary layer processes and local wind regimes in complex terrain was conducted in the vicinity of Lake Tekapo in the southern Alps of New Zealand, during two 1-month field campaigns in 1997 and 1999. The effects of the interaction of thermal and dynamic forcing were of specific interest, with a particular focus on the interaction of thermal forcing of differing scales. The rationale and objectives of the field and modeling program are described, along with the methodology used to achieve them. Specific research aims include improved knowledge of the role of surface forcing associated with varying energy balances across heterogeneous terrain, thermal influences on boundary layer and local wind development, and dynamic influences of the terrain through channeling effects. Data were collected using a network of surface meteorological and energy balance stations, radiosonde and pilot balloon soundings, tethered balloon and kite-based systems, sodar, and an instrumented light aircraft. These data are being used to investigate the energetics of surface heat fluxes, the effects of localized heating/cooling and advective processes on atmospheric boundary layer development, and dynamic channeling. A complementary program of numerical modeling includes application of the Regional Atmospheric Modeling System (RAMS) to case studies characterizing typical boundary layer structures and airflow patterns observed around Lake Tekapo. Some initial results derived from the special observation periods are used to illustrate progress made to date. In spite of the difficulties involved in obtaining good data and undertaking modeling experiments in such complex terrain, initial results show that surface thermal heterogeneity has a significant influence on local atmospheric structure and wind fields in the vicinity of the lake. This influence occurs particularly in the morning. However, dynamic channeling effects and the larger-scale thermal effect of the mountain region frequently override these more local features later in the day.

Online titrimetric and off-gas analysis for examining nitrification processes in wastewater treatment

The two steps of nitrification, namely the oxidation of ammonia to nitrite and nitrite to nitrate, often need to be considered separately in process studies. For a detailed examination, it is desirable to monitor the two-step sequence using online measurements. In this paper, the use of online titrimetric and off-gas analysis (TOGA) methods for the examination of the process is presented. Using the known reaction stoichiometry, combination of the measured signals (rates of hydrogen ion production, oxygen uptake and carbon dioxide transfer) allows the determination of the three key process rates, namely the ammonia consumption rate, the nitrite accumulation rate and the nitrate production rate. Individual reaction rates determined with the TOGA sensor under a number of operation conditions are presented. The rates calculated directly from the measured signals are compared with those obtained from offline liquid sample analysis. Statistical analysis confirms that the results from the two approaches match well. This result could not have been guaranteed using alternative online methods. As a case study, the influences of pH and dissolved oxygen (DO) on nitrite accumulation are tested using the proposed method. It is shown that nitrite accumulation decreased with increasing DO and pH. Possible reasons for these observations are discussed. (C) 2003 Elsevier Science Ltd. All rights reserved.