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.

Landscape structure influences tree density patterns in fragmented woodlands in semi-arid eastern Australia.

Landscape and local-scale influences are important drivers of plant community structure. However, their relative contribution and the degree to which they interact remain unclear. We quantified the extent to which landscape structure, within-patch habitat and their confounding effects determine post-clearing tree densities and composition in agricultural landscapes in eastern subtropical Australia. Landscape structure (incorporating habitat fragmentation and loss) and within-patch (site) features were quantified for 60 remnant patches of Eucalyptus populnea (Myrtaceae) woodland. Tree density and species for three ecological maturity classes (regeneration, early maturity, late maturity) and local site features were assessed in one 100 × 10 m plot per patch. All but one landscape characteristic was determined within a 1.3-km radius of plots; Euclidean nearest neighbour distance was measured inside a 5-km radius. Variation in tree density and composition for each maturity class was partitioned into independent landscape, independent site and joint effects of landscape and site features using redundancy analysis. Independent site effects explained more variation in regeneration density and composition than pure landscape effects; significant predictors were the proportion of early and late maturity trees at a site, rainfall and the associated interaction. Conversely, landscape structure explained greater variation in early and late maturity tree density and composition than site predictors. Area of remnant native vegetation within a landscape and patch characteristics (area, shape, edge contrast) were significant predictors of early maturity tree density. However, 31% of the explained variation in early mature tree differences represented confounding influences of landscape and local variables. We suggest that within-patch characteristics are important in influencing semi-arid woodland tree regeneration. However, independent and confounding effects of landscape structure resulting from previous vegetation clearing may have exerted a greater historical influence on older cohorts and should be accounted for when examining woodland dynamics across a broader range of environments.

Functional dynamics of the nitrogen balance of sorghum. II. Grain filling period.

Maintenance of green leaf area during grain filling can increase grain yield of sorghum grown under terminal water limitation. This 'stay-green' trait has been related to the nitrogen (N) supply-demand balance during grain filling. This study quantifies the N demand of grain and N translocation rates from leaves and stem and explores effects of genotype and N stress on onset and rate of leaf senescence during the grain filling period. Three hybrids differing in potential height were grown at three levels of N supply under well-watered conditions. Vertical profiles of biomass, leaf area, and N% of leaves, stem and grain were measured at regular intervals. Weekly SPAD chlorophyll readings on main shoot leaves were correlated with observed specific leaf nitrogen (SLN) to derive seasonal patterns of leaf N content. For all hybrids, individual grain N demand was sink determined and was initially met through N translocation from the stem and rachis. Only if this was insufficient did leaf N translocation occur. Maximum N translocation rates from leaves and stem were dependent on their N status. However, the supply of N at canopy scale was also related to the amount of leaf area senescing at any one time. This supply-demand framework for N dynamics explained effects of N stress and genotype on the onset and rate of leaf senescence.

Functional dynamics of the nitrogen balance of sorghum: I. N demand of vegetative plant parts.

Stay-green, an important trait for grain yield of sorghum grown under water limitation, has been associated with a high leaf nitrogen content at the start of grain filling. This study quantifies the N demand of leaves and stems and explores effects of N stress on the N balance of vegetative plant parts of three sorghum hybrids differing in potential crop height. The hybrids were grown under well-watered conditions at three levels of N supply. Vertical profiles of biomass and N% of leaves and stems, together with leaf size and number, and specific leaf nitrogen (SLN), were measured at regular intervals. The hybrids had similar minimum but different critical and maximum SLN, associated with differences in leaf size and N partitioning, the latter associated with differences in plant height. N demand of expanding new leaves was represented by critical SLN, and structural stem N demand by minimum stem N%. The fraction of N partitioned to leaf blades increased under N stress. A framework for N dynamics of leaves and stems is developed that captures effects of N stress and genotype on N partitioning and on critical and maximum SLN.

Molecular dynamics simulations of sliding in an Fe-Cu tribopair system

Experimental evidence suggests that high strain rates, stresses, strains and temperatures are experienced near sliding interfaces. The associated microstructural changes are due to several dynamic an interacting phenomena. 3D non-equilibrium molecular dynamics (MD) simulations of sliding were conducted with the aim of understanding the dynamic processes taking place in crystalline tribopairs, with a focus on plastic deformation and microstructural evolution. Embedded atom potentials were employed for simulating sliding of an Fe-Cu tribopair. Sliding velocity, crystal orientation and presence of lattice defects were some of the variables in these simulations. Extensive plastic deformation involving dislocation and twin activity, dynamic recrystallization, amorphization and/or nanocrystallization, mechanical mixing and material transfer were observed. Mechanical mixing in the vicinity of the sliding interface was observed even in the Fe-Cu system, which would cluster under equilibrium conditions, hinting at the ballistic nature of the process. Flow localization was observed at high velocities implying the possible role of adiabatic heating. The presence of preexisting defects (such as dislocations and interfaces) played a pivotal role in determining friction and microstructural evolution. The study also shed light on the relationship between adhesion and plastic deformation, and friction. Comparisons with experiments suggest that such simulations can indeed provide valuable insights that are difficult to obtain from experiments.

Structure and dynamics of two-dimensional sheared granular flows

The structure and dynamics of the two-dimensional linear shear flow of inelastic disks at high area fractions are analyzed. The event-driven simulation technique is used in the hard-particle limit, where the particles interact through instantaneous collisions. The structure (relative arrangement of particles) is analyzed using the bond-orientational order parameter. It is found that the shear flow reduces the order in the system, and the order parameter in a shear flow is lower than that in a collection of elastic hard disks at equilibrium. The distribution of relative velocities between colliding particles is analyzed. The relative velocity distribution undergoes a transition from a Gaussian distribution for nearly elastic particles, to an exponential distribution at low coefficients of restitution. However, the single-particle distribution function is close to a Gaussian in the dense limit, indicating that correlations between colliding particles have a strong influence on the relative velocity distribution. This results in a much lower dissipation rate than that predicted using the molecular chaos assumption, where the velocities of colliding particles are considered to be uncorrelated.

Sulphate ion dynamics in alpha-alums studied by esr at high pressures

The variation of zero-field splitting and linewidth of Cr3+ ion in KCr and KAI alums with hydrostatic pressure and with temperature is investigated. A model for the apparent phase transition is proposed on the basis of the reorientational motion of the SO2�4 groups.

Reassessing spatial and temporal dynamics of kangaroo populations

Long-running datasets from aerial surveys of kangaroos (Macropus giganteus, Macropus [uliginosus, Macropus robustus and Macropus rufus) across Queensland, New South Wales and South Australia have been analysed, seeking better predictors of rates of increase which would allow aerial surveys to be undertaken less frequently than annually. Early models of changes in kangaroo numbers in response to rainfall had shown great promise, but much variability. We used normalised difference vegetation index (NDVI) instead, reasoning that changes in pasture condition would provide a better predictor than rainfall. However, except at a fine scale, NDVI proved no better; although two linked periods of rainfall proved useful predictors of rates of increase, this was only in some areas for some species. The good correlations reported in earlier studies were a consequence of data dominated by large droughtinduced adult mortality, whereas over a longer time frame and where changes between years are less dramatic, juvenile survival has the strongest influence on dynamics. Further, harvesting, density dependence and competition with domestic stock are additional and important influences and it is now clear that kangaroo movement has a greater influence on population dynamics than had been assumed. Accordingly, previous conclusions about kangaroo populations as simple systems driven by rainfall need to be reassessed. Examination of this large dataset has permitted descriptions of shifts in distribution of three species across eastern Australia, changes in dispersion in response to rainfall, and an evaluation of using harvest statistics as an index of density and harvest rate. These results have been combined into a risk assessment and decision theory framework to identify optimal monitoring strategies.

Impacts of grazing management options on pasture and animal productivity in a Heteropogon contortus (black speargrass) pasture in central Queensland. 2. Population dynamics of Heteropogon contortus and Stylosanthes scabra cv. Seca.

The dynamics of Heteropogon contortus and Stylosanthes scabra cv. Seca populations were studied in a subset of treatments in an extensive grazing study conducted in central Queensland between 1988 and 2001. These treatments were 4 stocking rates in native pasture and 2 of these stocking rates in legume oversown and supplement/spring burning treatments. For the 1999-2000 summer, population data for H. contortus in 5 of these native pasture and supplement/burning treatments were compared with those for an additional burnt treatment. Seasonal rainfall throughout this study was below the long-term mean and mean annual pasture utilisation ranged from 24 to 61%. Increasing stocking rate from 5 to 2 ha/steer in native pasture reduced H. contortus plant density. Increasing stocking rate reduced seedling recruitment as a result of its effect on soil seedbanks. Seedling recruitment was the major determinant of change in plant density, although some individual H. contortus plants did survive throughout the study. Burning in spring 1999, particularly at light stocking rate, promoted seedling recruitment above that in both unburnt native and legume oversown pasture and resulted in increased H. contortus plant density. In the legume oversown treatments, S. scabra cv. Seca density increased rapidly from 15 plants/m2 in 1988 to 140 plants/m2 in 2001 following a lag phase between 1988 and 1993. This increased S. scabra density was associated with an eventual decline in H. contortus plant density through reduced seedling recruitment. It was concluded that H. contortus population density is sustainable at stocking rates of 4 and 5 ha/steer (30% pasture utilisation) and that spring burning at light stocking rate can promote H. contortus populations. Increasing densities of S. scabra need to be managed to prevent its dominance.

"Monesti on tullut ajateltua, että ihana, kamala kakara!" : Huomioidaanko päiväkodissa kansainvälisesti adoptoituun lapseen ja adoptioon liittyviä erityispiirteitä?

The purpose of this study was to examine the conceptions adoptive parents and kindergarten teachers have of adopted children in day care. Earlier Finnish research concerning international adoption has discussed adoption, adopted children and adoptive parents from the point of view of identity, racism and school. There has hardly been any research in Finland concerning adoption from the point of view of day care. The study focuses on how the special characteristics of adopted children are taken into consideration in day care. The ecological theory, in which a child's growth and learning environments are studied as part of the broader environment, provided the theoretical framework for the study. The study was based also on a second theory: intercultural competence or the capacity to interact with different people taking into consideration their individual characteristics and situations so that one understands and is understood. The first research problem consisted of describing the conceptions adoptive parents and day care teachers have of the special characteristics of adopted children. The second research problem was to examine adopted parents' expectations of kindergarten and of kindergarten teachers. The third research problem was to find out what kind of support day care teachers think they need when working with adopted children. The study focused on 20 adoptive parents and 15 kindergarten teachers. Research material was collected in the form of essay writings amounting to a total of 73 pages. The respondents wrote about their experiences with the help of a supporting questionnaire. The research method was qualitative and narrative. The results of the study show that adoptive parents have similar conceptions of adopted children's special characteristics concerning their character, social competence, adaptation and appearance. Adoptive parents had a better understanding of the dynamics of attachment than kindergarten teachers, in this study. When describing the special characters of adopted children, the teachers regarded adopted children in particular as needing special support. According to the results of the study, adoptive parents feel adoption matters and adoptive parenthood are something quite new and bewildering to day care teachers. Adoptive parents did not feel they got any support from the kindergarten in their parenthood. Adoptive parents regarded as particularly important the start of day care and the familiarisation process during which adopted children need time and special support due to earlier experiences of separation and abandonment. The findings of the study show that kindergarten teachers have little knowledge of adoption and the significance of adoption for a child's growth and development. Kindergarten teachers felt at a loss in dealing with adoption issues and clearly need additional support in the matter.

Net radiation dynamics: performance of 20 daily net radiation models as related to model structure and intricacy in two climates.

We compared daily net radiation (Rn) estimates from 19 methods with the ASCE-EWRI Rn estimates in two climates: Clay Center, Nebraska (sub-humid) and Davis, California (semi-arid) for the calendar year. The performances of all 20 methods, including the ASCE-EWRI Rn method, were then evaluated against Rn data measured over a non-stressed maize canopy during two growing seasons in 2005 and 2006 at Clay Center. Methods differ in terms of inputs, structure, and equation intricacy. Most methods differ in estimating the cloudiness factor, emissivity (e), and calculating net longwave radiation (Rnl). All methods use albedo (a) of 0.23 for a reference grass/alfalfa surface. When comparing the performance of all 20 Rn methods with measured Rn, we hypothesized that the a values for grass/alfalfa and non-stressed maize canopy were similar enough to only cause minor differences in Rn and grass- and alfalfa-reference evapotranspiration (ETo and ETr) estimates. The measured seasonal average a for the maize canopy was 0.19 in both years. Using a = 0.19 instead of a = 0.23 resulted in 6% overestimation of Rn. Using a = 0.19 instead of a = 0.23 for ETo and ETr estimations, the 6% difference in Rn translated to only 4% and 3% differences in ETo and ETr, respectively, supporting the validity of our hypothesis. Most methods had good correlations with the ASCE-EWRI Rn (r2 > 0.95). The root mean square difference (RMSD) was less than 2 MJ m-2 d-1 between 12 methods and the ASCE-EWRI Rn at Clay Center and between 14 methods and the ASCE-EWRI Rn at Davis. The performance of some methods showed variations between the two climates. In general, r2 values were higher for the semi-arid climate than for the sub-humid climate. Methods that use dynamic e as a function of mean air temperature performed better in both climates than those that calculate e using actual vapor pressure. The ASCE-EWRI-estimated Rn values had one of the best agreements with the measured Rn (r2 = 0.93, RMSD = 1.44 MJ m-2 d-1), and estimates were within 7% of the measured Rn. The Rn estimates from six methods, including the ASCE-EWRI, were not significantly different from measured Rn. Most methods underestimated measured Rn by 6% to 23%. Some of the differences between measured and estimated Rn were attributed to the poor estimation of Rnl. We conducted sensitivity analyses to evaluate the effect of Rnl on Rn, ETo, and ETr. The Rnl effect on Rn was linear and strong, but its effect on ETo and ETr was subsidiary. Results suggest that the Rn data measured over green vegetation (e.g., irrigated maize canopy) can be an alternative Rn data source for ET estimations when measured Rn data over the reference surface are not available. In the absence of measured Rn, another alternative would be using one of the Rn models that we analyzed when all the input variables are not available to solve the ASCE-EWRI Rn equation. Our results can be used to provide practical information on which method to select based on data availability for reliable estimates of daily Rn in climates similar to Clay Center and Davis.

Collaborative learning approach to introduce computational fluid dynamics

This paper presents the design, implementation and evaluation of a collaborative learning activity designed to replace traditional face-to-face lectures in a large classroom. This activity aims to better engage the students with their learning and improve the students’ experience and outcomes. This project is implemented in the Fluid Mechanics unit of the Mechanical Engineering degree at the Queensland University of Technology to introduce students with the concept, terminology and process of Computational Fluid Dynamics (CFD). The approach integrates a constructive collaborative assignment which is a key element in the overall quality of teaching and learning, and an integral component of the students’ experience. A detailed survey, given to the students, showed an overall high level of satisfaction. However, the results also highlighted the gap between students’ expectations both for contents and assignment and teacher expectations. Discussions to address this issue are presented in the paper based on a critical reflection.

A molecular dynamics investigation of the structure of PbO---PbF2 glasses*1

The structure of PbO---PbF2 glasses has been studied using molecular dynamics (MD). The existence of [OPb4] structural units is observed over the entire glass-forming range, in conformity with a model proposed earlier based on various structural investigations of this system. Various other features of the structural model are also supported by the MD calculations.

A conformational approach to the study of the dynamics of enzyme inhibition: studies on thermolysin

The preferred conformations of β-phenylpropionyl-Image -phenylalanine (β-PPP) and N-carbobenzoxy-L-phenylalanine (Cbz-Phe), two inhibitors of thermolysin, have been determined by computing potential energy using empirial potential energy functions. Of the 15 to 20 conformations that are favoured for each of these inhibitors only a few have the right conformation to reach the active site of the enzyme. The conformer of β-PPP that initiates binding with the enzyme is different from the bound one, while for Cbz-Phe the bound and initiating conformers are quite similar. Thus, β-PPP favours the ‘induced fit’ model while Cbz-Phe follows the ‘lock and key’ model of binding. The inhibitors differ in their alignment at the active site.