Special Issue on Spatial Moment Techniques for Modelling Biological Processes

Over the last two decades, there has been an increasing awareness of, and interest in, the use of spatial moment techniques to provide insight into a range of biological and ecological processes. Models that incorporate spatial moments can be viewed as extensions of mean-field models. These mean-field models often consist of systems of classical ordinary differential equations and partial differential equations, whose derivation, at some point, hinges on the simplifying assumption that individuals in the underlying stochastic process encounter each other at a rate that is proportional to the average abundance of individuals. This assumption has several implications, the most striking of which is that mean-field models essentially neglect any impact of the spatial structure of individuals in the system. Moment dynamics models extend traditional mean-field descriptions by accounting for the dynamics of pairs, triples and higher n-tuples of individuals. This means that moment dynamics models can, to some extent, account for how the spatial structure affects the dynamics of the system in question.

Modelling the movement of interacting cell populations: A moment dynamics approach

Mathematical models describing the movement of multiple interacting subpopulations are relevant to many biological and ecological processes. Standard mean-field partial differential equation descriptions of these processes suffer from the limitation that they implicitly neglect to incorporate the impact of spatial correlations and clustering. To overcome this, we derive a moment dynamics description of a discrete stochastic process which describes the spreading of distinct interacting subpopulations. In particular, we motivate our model by mimicking the geometry of two typical cell biology experiments. Comparing the performance of the moment dynamics model with a traditional mean-field model confirms that the moment dynamics approach always outperforms the traditional mean-field approach. To provide more general insight we summarise the performance of the moment dynamics model and the traditional mean-field model over a wide range of parameter regimes. These results help distinguish between those situations where spatial correlation effects are sufficiently strong, such that a moment dynamics model is required, from other situations where spatial correlation effects are sufficiently weak, such that a traditional mean-field model is adequate.

Development of a multi-locus sequence typing scheme for avian isolates of Pasteurella multocida.

A total of 63 isolates of Pasteurella multocida from Australian poultry, all associated with fowl cholera outbreaks, and three international reference strains, representing the three subspecies within P. multocida were used to develop a multi-locus sequence typing scheme. Primers were designed for conserved regions of seven house-keeping enzymes - adk, est, gdh, mdh, pgi, pmi and zwf - and internal fragments of 570-784 bp were sequenced for all isolates and strains. The number of alleles at the different loci ranged from 11 to 20 and a total of 29 allelic profiles or sequence types were recognised amongst the 66 strains. There was a strong concordance between the MLST data and the existing multi-locus enzyme electrophoresis and ribotyping data. When used to study a sub-set of isolates with a known detailed epidemiological history, the MLST data matched the results given by restriction endonuclease analysis, pulsed-field gel electrophoresis, ribotyping and REP-PCR. The MLST scheme provides a high level of resolution and is an excellent tool for studying the population structure and epidemiology of P. multocida.

Hip power analysis in individuals with transfemoral amputation: A different strategy different from stabilisation during gait stance

Introduction and Objectives Joint moments and joint powers during gait are widely used to determine the effects of rehabilitation programs as well as prosthetic fitting. Following the definition of power (dot product of joint moment and joint angular velocity) it has been previously proposed to analyse the 3D angle between both vectors, αMw. Basically, joint power is maximised when both vectors are parallel and cancelled when both vectors are orthogonal. In other words, αMw < 60° reveals a propulsion configuration (more than 50% of the moment contribute to positive power) while αMw > 120° reveals a resistance configuration (more than 50% of the moment contribute to negative power). A stabilisation configuration (less than 50% of the moment contribute to power) corresponds to 60° < αMw < 120°. Previous studies demonstrated that hip joints of able-bodied adults (AB) are mainly in a stabilisation configuration (αMw about 90°) during the stance phase of gait. [1, 2] Individuals with transfemoral amputation (TFA) need to maximise joint power at the hip while controlling the prosthetic knee during stance. Therefore, we tested the hypothesis that TFAs should adopt a strategy that is different from a continuous stabilisation. The objective of this study was to compute joint power and αMw for TFA and to compare them with AB. Methods Three trials of walking at self-selected speed were analysed for 8 TFAs (7 males and 1 female, 46±10 years old, 1.78±0.08 m 82±13 kg) and 8 ABs (males, 25±3 years old, 1.75±0.04, m 67±6 kg). The joint moments are computed from a motion analysis system (Qualisys, Goteborg, Sweden) and a multi-axial transducer (JR3, Woodland, USA) mounted above the prosthetic knee for TFAs and from a motion analysis system (Motion Analysis, Santa Rosa, USA) and force plates (Bertec, Columbus, USA) for ABs. The TFAs were fitted with an OPRA (Integrum, AB, Gothengurg, Sweden) osseointegrated implant system and their prosthetic designs include pneumatic, hydraulic and microprocessor knees. Previous studies showed that the inverse dynamics computed from the multi-axial transducer is the proper method considering the absorption at the foot and resistance at the knee. Results The peak of positive power at loading response (H1) was earlier and lower for TFA compared to AB. Although the joint power is lower, the 3D angle between joint moment and joint angular velocity, αMw, reveals an obvious propulsion configuration (mean αMw about 20°) for TFA compared to a stabilisation configuration (mean αMw about 70°) for AB. The peaks of negative power at midstance (H2) and of positive power at preswing / initial swing (H3) occurred later, lower and longer for TFA compared to AB. Again, the joint powers are lower for TFA but, in this case, αMw is almost comparable (with a time lag), demonstrating a stabilisation (almost a resistance for TFA, mean αMw about 120°) and a propulsion configuration, respectively. The swing phase is not analysed in the present study. Conclusion The analysis of hip joint power may indicate that TFAs demonstrated less propulsion and resistance than ABs during the stance phase of gait. This is true from a quantitative point of view. On the contrary, the 3D angle between joint moment and joint angular velocity, αMw, reveals that TFAs have a remarkable propulsion strategy at loading response and almost a resistance strategy at midstance while ABs adopted a stabilisation strategy. The propulsion configuration, with αMw close to 0°, seems to aim at maximising the positive joint power. The configuration close to resistance, with αMw far from 180°, might aim at unlocking the prosthetic knee before swing while minimising the negative power. This analysis of both joint power and 3D angle between the joint moment and the joint angular velocity provides complementary insights into the gait strategies of TFA that can be used to support evidence-based rehabilitation and fitting of prosthetic components.

A parallel-serial binary arbitration scheme for collision-free multi-access techniques

Multi-access techniques are widely used in computer networking and distributed multiprocessor systems. On-the-fly arbitration schemes permit one of the many contenders to access the medium without collisions. Serial arbitration is cost effective but is slow and hence unsuitable for high-speed multiprocessor environments supporting very high data transfer rates. A fully parallel arbitration scheme takes less time but is not practically realisable for large numbers of contenders. In this paper, a generalised parallel-serial scheme is proposed which significantly reduces the arbitration time and is practically realisable.

Distributed Scheduling for Multi-Hop Wireless Networks

We study the performance of greedy scheduling in multihop wireless networks where the objective is aggregate utility maximization. Following standard approaches, we consider the dual of the original optimization problem. Optimal scheduling requires selecting independent sets of maximum aggregate price, but this problem is known to be NP-hard. We propose and evaluate a simple greedy heuristic. We suggest how the greedy heuristic can be implemented in a distributed manner. We evaluate an analytical bound in detail, for the special case of a line graph and also provide a loose bound on the greedy heuristic for the case of an arbitrary graph.

MT10049 Pheromone component of a multi target approach to fruitspotting bug management

A pheromone-based trapping system will be developed for both A. lutescens and A. nitida to improve insecticide timing and to rationalise use.

Multi-time scale analysis applied to long time simulation of power systems

The results are presented of applying multi-time scale analysis using the singular perturbation technique for long time simulation of power system problems. A linear system represented in state-space form can be decoupled into slow and fast subsystems. These subsystems can be simulated with different time steps and then recombined to obtain the system response. Simulation results with a two-time scale analysis of a power system show a large saving in computational costs.

Robustness of the partitions in grouping multi-item inventories

A recent work obtained closed-form solutions to the.problem of optimally grouping a multi-item inventory into subgroups with a common order cycle per group, when the distribution by value of the inventory could be described by a Pareto function. This paper studies the sensitivity of the optimal subgroup boundaries so obtained. Closed-form expressions have been developed to find intervals for the subgroup boundaries for any given level of suboptimality. Graphs have been provided to aid the user in selecting a cost-effective level of aggregation and choosing appropriate subgroup boundaries for a whole range of inventory distributions. The results of sensitivity analyses demonstrate the availability of flexibility in the partition boundaries and the cost-effectiveness of any stock control system through three groups, and thus also provide a theoretical support to the intuitive ABC system of classifying the items.

A multi-field bio-economic model of irrigated grain-cotton farming systems.

We present a participatory modelling framework that integrates information from interviews and discussions with farmers and consultants, with dynamic bio-economic models to answer complex questions on the allocation of limited resources at the farm business level. Interviews and discussions with farmers were used to: describe the farm business; identify relevant research questions; identify potential solutions; and discuss and learn from the whole-farm simulations. The simulations are done using a whole-farm, multi-field configuration of APSIM (APSFarm). APSFarm results were validated against farmers' experience. Once the model was accepted by the participating farmers as a fair representation of their farm business, the model was used to explore changes in the tactical or strategic management of the farm and results were then discussed to identify feasible options for improvement. Here we describe the modelling framework and present an example of the application of integrative whole farm system tools to answer relevant questions from an irrigated farm business case study near Dalby (151.27E - 27.17S), Queensland, Australia. Results indicated that even though cotton crops generates more farm income per hectare a more diversified rotation with less cotton would be relatively more profitable, with no increase in risk, as a more cotton dominated traditional rotation. Results are discussed in terms of the benefits and constraints from developing and applying more integrative approaches to represent farm businesses and their management in participatory research projects with the aim of designing more profitable and sustainable irrigated farming systems.

A multi-faceted approach for quantifying the estuarine-nearshore transition in the life cycle of the bull shark, Carcharhinus leucas.

Understanding the ontogenetic habitat linkages of sharks is important for conservation and managing human interactions. We used acoustic telemetry, catch data, elemental and stable isotope signatures and dietary analyses to investigate ontogenetic habitat use in south-east Queensland, Australia, by the bull shark Carcharhinus leucas, a IUCN 'near-threatened' species that is implicated in many shark attacks on humans in urban estuaries. Sequential analyses for delta(15)N and delta(13)C of vertebrae from five adult C. leucas and laser ablation inductively coupled plasma mass spectrometry (LA-ICPMS) for elemental composition from 23 C. leucas, including a pregnant female, were also used to trace ontogenetic habitat dependence. Acoustic telemetry indicated large juvenile and subadult C. leucas remained in estuarine habitats. delta(15)N values across shark vertebrae showed an ontogenetic shift in diet with total length (TL), confirmed by stomach contents. LA-ICPMS data reflected the ontogenetic movements of C. leucas from natal habitats. Differences among adults were gender related. Shifts in habitat use by subadults were correlated with a sigmoidal delta(13)C relationship with TL. C. leucas have a multipartite, stage-specific dependency in their transition between habitats along the freshwater-estuarine-marine continuum, making them particularly susceptible to the habitat alteration that is occurring globally.

Non-Destructive Method to Estimate the Moisture Content in Bread using Multi-Channel Electrical Impedance Spectroscopy

Bread undergoes several physicochemical changes during storage that results in a rapid loss of freshness. These changes depend on moisture content present in bread product. An instrument based on electrical impedance spectroscopy technique is developed to estimate moisture content of bread at different zones using designed multi-channel ring electrodes. A dedicated AT89S52 microcontroller and associated peripherals are employed for hardware. A constant current is applied across bread loaf through central pair of electrodes and developed potential across different zones of bread loaf are measured using remaining four ring electrode pairs. These measured values of voltage and current are used to measure the impedance at each zone. Electrical impedance behavior of the bread loaf at crust and crumb is investigated during storage. A linear relationship is observed between the measured impedance and moisture content present in crust and crumb of bread loaf during storage of 120 hours.

Management of multi-scale forest resource data over time

During the last decades there has been a global shift in forest management from a focus solely on timber management to ecosystem management that endorses all aspects of forest functions: ecological, economic and social. This has resulted in a shift in paradigm from sustained yield to sustained diversity of values, goods and benefits obtained at the same time, introducing new temporal and spatial scales into forest resource management. The purpose of the present dissertation was to develop methods that would enable spatial and temporal scales to be introduced into the storage, processing, access and utilization of forest resource data. The methods developed are based on a conceptual view of a forest as a hierarchically nested collection of objects that can have a dynamically changing set of attributes. The temporal aspect of the methods consists of lifetime management for the objects and their attributes and of a temporal succession linking the objects together. Development of the forest resource data processing method concentrated on the extensibility and configurability of the data content and model calculations, allowing for a diverse set of processing operations to be executed using the same framework. The contribution of this dissertation to the utilisation of multi-scale forest resource data lies in the development of a reference data generation method to support forest inventory methods in approaching single-tree resolution.

I/O Device Virtualization in the multi-core era, a QoS perspective

In this paper, we propose an extension to the I/O device architecture, as recommended in the PCI-SIG IOV specification, for virtualizing network I/O devices. The aim is to enable fine-grained controls to a virtual machine on the I/O path of a shared device. The architecture allows native access of I/O devices to virtual machines and provides device level QoS hooks for controlling VM specific device usage. For evaluating the architecture we use layered queuing network (LQN) models. We implement the architecture and evaluate it using simulation techniques, on the LQN model, to demonstrate the benefits. With the architecture, the benefit for network I/O is 60% more than what can be expected on the existing architecture. Also, the proposed architecture improves scalability in terms of the number of virtual machines intending to share the I/O device.