Ross River virus disease clusters and spatial relationship with mosquito biting exposure in Redland Shire, southern Queensland, Australia

The spatial heterogeneity in the risk of Ross River virus (family Togaviridae, genus Alphavirus, RRV) disease, the most common mosquito-borne disease in Australia, was examined in Redland Shire in southern Queensland, Australia. Disease cases, complaints from residents of intense mosquito biting exposure, and human population data were mapped using a geographic information system. Surface maps of RRV disease age-sex standardized morbidity ratios and mosquito biting complaint morbidity ratios were created. To determine whether there was significant spatial variation in disease and complaint patterns, a spatial scan analysis method was used to test whether the number of cases and complaints was distributed according to underlying population at risk. Several noncontiguous areas in proximity to productive saline water habitats of Aedes vigilax (Skuse), a recognized vector of RRV, had higher than expected numbers of RRV disease cases and complaints. Disease rates in human populations in areas which had high numbers of adult Ae. vigilax in carbon dioxide- and octenol-baited light traps were up to 2.9 times those in areas that rarely had high numbers of mosquitoes. It was estimated that targeted control of adult Ae. vigilax in these high-risk areas could potentially reduce the RRV disease incidence by an average of 13.6%. Spatial correlation was found between RRV disease risk and complaints from residents of mosquito biting. Based on historical patterns of RRV transmission throughout Redland Shire and estimated future human population growth in areas with higher than average RRV disease incidence, it was estimated that RRV incidence rates will increase by 8% between 2001 and 2021. The use of arbitrary administrative areas that ranged in size from 4.6 to 318.3 km2, has the potential to mask any small scale heterogeneity in disease patterns. With the availability of georeferenced data sets and high-resolution imagery, it is becoming more feasible to undertake spatial analyses at relatively small scales.

Do predators aggregate in response to pest density in agroecosystems? Assessing within-field spatial patterns

1. The spatial heterogeneity of predator populations is an important component of ecological theories pertaining to predator-prey dynamics. Most studies within agricultural fields show spatial correlation (positive or negative) between mean predator numbers and prey abundance across a whole field over time but generally ignore the within-field spatial dimension. We used explicit spatial mapping to determine if generalist predators aggregated within a soybean field, the size of these aggregations and if predator aggregation was associated with pest aggregation, plant damage and predation rate. 2. The study was conducted at Gatton in the Lockyer Valley, 90 km west of Brisbane, Australia. Intensive sampling grids were used to investigate within-field spatial patterns. The first row of each grid was located in a lucerne field (10 m from interface) and the remaining rows were in an adjacent soybean field. At each point on the grid the abundance of foliage-dwelling and ground-dwelling pests and predators was measured, predation rates [using sentinel Helicoverpa armigera (Hubner) egg cards] and plant damage were estimated. Eight grids were sampled across two summer cropping seasons (2000/01, 2001/02). 3. Predators exhibited strong spatial patterning with regions of high and low abundance and activity within what are considered to be uniform soybean fields. Ground-dwelling and foliage-dwelling predators were often aggregated in patches approximately 40 m across. 4. Lycosidae (wolf spiders) displayed aggregation and were consistently more abundant within the lucerne, with a decreasing trap catch with distance from the lucrene/soybean interface. This trend was consistent between subsequent grids in a single field and between fields. 5. The large amount of spatial variability in within-field arthropod abundance (pests and predators) and activity (egg predation and plant damage) indicates that whole field averages were misleading. This result has serious implications for sampling of arthropod abundance and pest management decision-making based on scouting data. 6. There was a great deal of temporal change in the significant spatial patterns observed within a field at each sampling time point during a single season. Predator and pest aggregations observed in these fields were generally not stable for the entire season. 7. Predator aggregation did not correlate consistently with pest aggregation, plant damage or predation rate. Spatial patterns in predator abundance were not associated consistently with any single parameter measured. The most consistent positive association was between foliage-dwelling predators and pests (significant in four of seven grids). Inferring associations between predators and prey based on an intensive one-off sampling grid is difficult, due to the temporal variability in the abundance of each group. 8. Synthesis and applications. This study demonstrated that generalist predator populations are rarely distributed randomly and field edges and adjacent crops can have an influence on within-field predator abundance. This must be considered when estimating arthropod (pest and predator) abundance from a set of samples taken at random locations within a field.

Habitat selection and population regulation in temporally fluctuating environments

Understanding and predicting the distribution of organisms in heterogeneous environments lies at the heart of ecology, and the theory of density-dependent habitat selection (DDHS) provides ecologists with an inferential framework linking evolution and population dynamics. Current theory does not allow for temporal variation in habitat quality, a serious limitation when confronted with real ecological systems. We develop both a stochastic equivalent of the ideal free distribution to study how spatial patterns of habitat use depend on the magnitude and spatial correlation of environmental stochasticity and also a stochastic habitat selection rule. The emerging patterns are confronted with deterministic predictions based on isodar analysis, an established empirical approach to the analysis of habitat selection patterns. Our simulations highlight some consistent patterns of habitat use, indicating that it is possible to make inferences about the habitat selection process based on observed patterns of habitat use. However, isodar analysis gives results that are contingent on the magnitude and spatial correlation of environmental stochasticity. Hence, DDHS is better revealed by a measure of habitat selectivity than by empirical isodars. The detection of DDHS is but a small component of isodar theory, which remains an important conceptual framework for linking evolutionary strategies in behavior and population dynamics.

Of sheep and rain: large-scale population dynamics of the red kangaroo

1. We analysed time-series data from populations of red kangaroos (Macropus rufus, Desmarest) inhabiting four areas in the pastoral zone of South Australia. We formulated a set of a priori models to disentangle the relative effects of the covariates: rainfall, harvesting, intraspecific competition, and domestic herbivores, on kangaroo population-growth rate. 2. The statistical framework allowed for spatial variation in the growth-rate parameters, response to covariates, and environmental variability, as well as spatially correlated error terms due to shared environment. 3. The most parsimonious model included all covariates but no area-specific parameter values, suggesting that kangaroo densities respond in the same way to the covariates across the areas. 4. The temporal dynamics were spatially correlated, even after taking into account the potentially synchronizing effect of rainfall, harvesting and domestic herbivores. 5. Counter-intuitively, we found a positive rather than negative effect of domestic herbivore density on the population-growth rate of kangaroos. We hypothesize that this effect is caused by sheep and cattle acting as a surrogate for resource availability beyond rainfall. 6. Even though our system is well studied, we must conclude that approximating resources by surrogates such as rainfall is more difficult than previously thought. This is an important message for studies of consumer-resource systems and highlights the need to be explicit about population processes when analysing population patterns.

First-principles quantum dynamics in interacting Bose gases II: stochastic gauges

First principles simulations of the quantum dynamics of interacting Bose gases using the stochastic gauge representation are analysed. In a companion paper, we showed how the positive-P representation can be applied to these problems using stochastic differential equations. That method, however, is limited by increased sampling error as time evolves. Here, we show how the sampling error can be greatly reduced and the simulation time significantly extended using stochastic gauges. In particular, local stochastic gauges (a subset) are investigated. Improvements are confirmed in numerical calculations of single-, double- and multi-mode systems in the weak-mode coupling regime. Convergence issues are investigated, including the recognition of two modes by which stochastic equations produced by phase-space methods in general can diverge: movable singularities and a noise-weight relationship. The example calculated here displays wave-like behaviour in spatial correlation functions propagating in a uniform 1D gas after a sudden change in the coupling constant. This could in principle be tested experimentally using Feshbach resonance methods.

The cost of uniqueness in groundwater model calibration

Calibration of a groundwater model requires that hydraulic properties be estimated throughout a model domain. This generally constitutes an underdetermined inverse problem, for which a Solution can only be found when some kind of regularization device is included in the inversion process. Inclusion of regularization in the calibration process can be implicit, for example through the use of zones of constant parameter value, or explicit, for example through solution of a constrained minimization problem in which parameters are made to respect preferred values, or preferred relationships, to the degree necessary for a unique solution to be obtained. The cost of uniqueness is this: no matter which regularization methodology is employed, the inevitable consequence of its use is a loss of detail in the calibrated field. This, ill turn, can lead to erroneous predictions made by a model that is ostensibly well calibrated. Information made available as a by-product of the regularized inversion process allows the reasons for this loss of detail to be better understood. In particular, it is easily demonstrated that the estimated value for an hydraulic property at any point within a model domain is, in fact, a weighted average of the true hydraulic property over a much larger area. This averaging process causes loss of resolution in the estimated field. Where hydraulic conductivity is the hydraulic property being estimated, high averaging weights exist in areas that are strategically disposed with respect to measurement wells, while other areas may contribute very little to the estimated hydraulic conductivity at any point within the model domain, this possibly making the detection of hydraulic conductivity anomalies in these latter areas almost impossible. A study of the post-calibration parameter field covariance matrix allows further insights into the loss of system detail incurred through the calibration process to be gained. A comparison of pre- and post-calibration parameter covariance matrices shows that the latter often possess a much smaller spectral bandwidth than the former. It is also demonstrated that, as all inevitable consequence of the fact that a calibrated model cannot replicate every detail of the true system, model-to-measurement residuals can show a high degree of spatial correlation, a fact which must be taken into account when assessing these residuals either qualitatively, or quantitatively in the exploration of model predictive uncertainty. These principles are demonstrated using a synthetic case in which spatial parameter definition is based oil pilot points, and calibration is Implemented using both zones of piecewise constancy and constrained minimization regularization. (C) 2005 Elsevier Ltd. All rights reserved.

Geochronology of the Pengjiakuang and Rushan gold deposits, Eastern Jiaodong Gold Province, Northeastern China: Implications for regional mineralization and geodynamic setting

The Jiaodong gold province, the largest gold-producing district in China, is located in the jiaodong peninsula at the eastern margin of the North China craton and bounded by the continental scale Tan-Lu fault, 40 kin to the west. Previous geochronological studies suggest that pervasive gold deposition took place in the western part of the province between 122 and 119 Ma. Here we report high-quality Ar-40/Ar-39 ages of the Pengjiakuang and Rushan deposits from the eastern part of the jiaodong gold province, placing additional chronological constraints on the timing of regional mineralization. Seven sericite grains extracted from auriferous alteration assemblages at the Pengiiakuang deposit yielded well-defined plateau ages between 120.9 +/- 0.4 and 119.1 +/- 0.2 Ma (2 sigma). Three separates of igneous biotite from a sample of the Queshan gneissic granite, adjacent to the Pengjiakuang deposit, gave reproducible plateau ages of 124.6 +/- 0.6 to 123.9 +/- 0.4 Ma (2 sigma). Six sericite separates front two samples in the Rushan deposit yielded Ar-40/Ar-39 plateau ages at 109.3 +/- 0.3 to 107.7 +/- 0.5 Ma (2 sigma), whereas biotite from the Kunyushan monzogranite that hosts the Rushan deposit had plateau ages ranging from 129.0 +/- 0.6 to 126.9 +/- 0.6 Ma (3 separates front one sample). The apparent age gap between hydrothermal sericite and magmtic biotite from both deposits, together with the similar argon closure temperatures for these mica minerals, suggest that gold mineralization had no direct relationship to the granitoid magmatism. Instead, gold deposition coincided with the emplacement of mafic to intermediate dikes widespread in the jiaodong gold province, which have been dated at ca. 122 to 119 Ma and, less commonly, at 110 to 102 Ma. The new Ar-40/Ar-39 ages from the eastern jiaodong peninsula, when combined with published data from the western part suggest that gold mineralization was broadly contemporaneous throughout the district. The Early Cretaceous gold mineralization also is widely developed in four other major gold districts along the Tan-Lu fault. The temporal and spatial correlation of these gold deposits with mafic to intermediate dikes commonly found in most mineralized areas, the presence of well-documented metamorphic core complexes and half-graben basins along the Tan-Lu fault, and voluminous basalts therein, suggest that the Early Cretaceous was an important period of lithospheric extension, possibly caused by the late Mesozoic lithospheric thinning beneath the eastern block of the North China craton. Lithospheric thinning and extension could have resulted in abnormally high heat and fluid fluxes necessary for large-scaled gold mineralization.

Stochastic modelling and simulations for solute transport in porous media

A stochastic model for solute transport in aquifers is studied based on the concepts of stochastic velocity and stochastic diffusivity. By applying finite difference techniques to the spatial variables of the stochastic governing equation, a system of stiff stochastic ordinary differential equations is obtained. Both the semi-implicit Euler method and the balanced implicit method are used for solving this stochastic system. Based on the Karhunen-Loeve expansion, stochastic processes in time and space are calculated by means of a spatial correlation matrix. Four types of spatial correlation matrices are presented based on the hydraulic properties of physical parameters. Simulations with two types of correlation matrices are presented.

Fine-scale spatial genetic correlation analyses reveal strong female philopatry within a brush-tailed rock-wallaby colony in southeast Queensland

We combine spatial data on home ranges of individuals and microsatellite markers to examine patterns of fine-scale spatial genetic structure and dispersal within a brush-tailed rock-wallaby (Petrogale penicillata) colony at Hurdle Creek Valley, Queensland. Brush-tailed rock-wallabies were once abundant and widespread throughout the rocky terrain of southeastern Australia; however, populations are nearly extinct in the south of their range and in decline elsewhere. We use pairwise relatedness measures and a recent multilocus spatial autocorrelation analysis to test the hypotheses that in this species, within-colony dispersal is male-biased and that female philopatry results in spatial clusters of related females within the colony. We provide clear evidence for strong female philopatry and male-biased dispersal within this rock-wallaby colony. There was a strong, significant negative correlation between pairwise relatedness and geographical distance of individual females along only 800 m of cliff line. Spatial genetic autocorrelation analyses showed significant positive correlation for females in close proximity to each other and revealed a genetic neighbourhood size of only 600 m for females. Our study is the first to report on the fine-scale spatial genetic structure within a rock-wallaby colony and we provide the first robust evidence for strong female philopatry and spatial clustering of related females within this taxon. We discuss the ecological and conservation implications of our findings for rock-wallabies, as well as the importance of fine-scale spatial genetic patterns in studies of dispersal behaviour.

Optimal sampling strategy for estimation of spatial genetic structure in tree populations

Fine-scale spatial genetic structure (SGS) in natural tree populations is largely a result of restricted pollen and seed dispersal. Understanding the link between limitations to dispersal in gene vectors and SGS is of key interest to biologists and the availability of highly variable molecular markers has facilitated fine-scale analysis of populations. However, estimation of SGS may depend strongly on the type of genetic marker and sampling strategy (of both loci and individuals). To explore sampling limits, we created a model population with simulated distributions of dominant and codominant alleles, resulting from natural regeneration with restricted gene flow. SGS estimates from subsamples (simulating collection and analysis with amplified fragment length polymorphism (AFLP) and microsatellite markers) were correlated with the 'real' estimate (from the full model population). For both marker types, sampling ranges were evident, with lower limits below which estimation was poorly correlated and upper limits above which sampling became inefficient. Lower limits (correlation of 0.9) were 100 individuals, 10 loci for microsatellites and 150 individuals, 100 loci for AFLPs. Upper limits were 200 individuals, five loci for microsatellites and 200 individuals, 100 loci for AFLPs. The limits indicated by simulation were compared with data sets from real species. Instances where sampling effort had been either insufficient or inefficient were identified. The model results should form practical boundaries for studies aiming to detect SGS. However, greater sample sizes will be required in cases where SGS is weaker than for our simulated population, for example, in species with effective pollen/seed dispersal mechanisms.

Linking riparian vegetation spatial structure in Australian tropical savannas to ecosystem health indicators: semi-variogram analysis of high spatial resolution satellite imagery

Government agencies responsible for riparian environments are assessing the combined utility of field survey and remote sensing for mapping and monitoring indicators of riparian zone health. The objective of this work was to determine if the structural attributes of savanna riparian zones in northern Australia can be detected from commercially available remotely sensed image data. Two QuickBird images and coincident field data covering sections of the Daly River and the South Alligator River - Barramundie Creek in the Northern Territory were used. Semi-variograms were calculated to determine the characteristic spatial scales of riparian zone features, both vegetative and landform. Interpretation of semi-variograms showed that structural dimensions of riparian environments could be detected and estimated from the QuickBird image data. The results also show that selecting the correct spatial resolution and spectral bands is essential to maximize the accuracy of mapping spatial characteristics of savanna riparian features. The distribution of foliage projective cover of riparian vegetation affected spectral reflectance variations in individual spectral bands differently. Pan-sharpened image data enabled small-scale information extraction (< 6 m) on riparian zone structural parameters. The semi-variogram analysis results provide the basis for an inversion approach using high spatial resolution satellite image data to map indicators of savanna riparian zone health.