Spatial and temporal patterns of locally-acquired dengue transmission in northern Queensland, Australia, 1993-2012

BACKGROUND: Dengue has been a major public health concern in Australia since it re-emerged in Queensland in 1992-1993. We explored spatio-temporal characteristics of locally-acquired dengue cases in northern tropical Queensland, Australia during the period 1993-2012. METHODS: Locally-acquired notified cases of dengue were collected for northern tropical Queensland from 1993 to 2012. Descriptive spatial and temporal analyses were conducted using geographic information system tools and geostatistical techniques. RESULTS: 2,398 locally-acquired dengue cases were recorded in northern tropical Queensland during the study period. The areas affected by the dengue cases exhibited spatial and temporal variation over the study period. Notified cases of dengue occurred more frequently in autumn. Mapping of dengue by statistical local areas (census units) reveals the presence of substantial spatio-temporal variation over time and place. Statistically significant differences in dengue incidence rates among males and females (with more cases in females) (χ(2) = 15.17, d.f.  = 1, p<0.01). Differences were observed among age groups, but these were not statistically significant. There was a significant positive spatial autocorrelation of dengue incidence for the four sub-periods, with the Moran's I statistic ranging from 0.011 to 0.463 (p<0.01). Semi-variogram analysis and smoothed maps created from interpolation techniques indicate that the pattern of spatial autocorrelation was not homogeneous across the northern Queensland. CONCLUSIONS: Tropical areas are potential high-risk areas for mosquito-borne diseases such as dengue. This study demonstrated that the locally-acquired dengue cases have exhibited a spatial and temporal variation over the past twenty years in northern tropical Queensland, Australia. Therefore, this study provides an impetus for further investigation of clusters and risk factors in these high-risk areas.

Spatial incidence of dengue infections in Queensland, Australia - reply

This study aimed to investigate the spatial clustering and dynamic dispersion of dengue incidence in Queensland, Australia. We used Moran's I statistic to assess the spatial autocorrelation of reported dengue cases. Spatial empirical Bayes smoothing estimates were used to display the spatial distribution of dengue in postal areas throughout Queensland. Local indicators of spatial association (LISA) maps and logistic regression models were used to identify spatial clusters and examine the spatio-temporal patterns of the spread of dengue. The results indicate that the spatial distribution of dengue was clustered during each of the three periods of 1993–1996, 1997–2000 and 2001–2004. The high-incidence clusters of dengue were primarily concentrated in the north of Queensland and low-incidence clusters occurred in the south-east of Queensland. The study concludes that the geographical range of notified dengue cases has significantly expanded in Queensland over recent years.

Flying-Fox Species Density-A Spatial Risk Factor for Hendra Virus Infection in Horses in Eastern Australia

Hendra virus causes sporadic but typically fatal infection in horses and humans in eastern Australia. Fruit-bats of the genus Pteropus (commonly known as flying-foxes) are the natural host of the virus, and the putative source of infection in horses; infected horses are the source of human infection. Effective treatment is lacking in both horses and humans, and notwithstanding the recent availability of a vaccine for horses, exposure risk mitigation remains an important infection control strategy. This study sought to inform risk mitigation by identifying spatial and environmental risk factors for equine infection using multiple analytical approaches to investigate the relationship between plausible variables and reported Hendra virus infection in horses. Spatial autocorrelation (Global Moran’s I) showed significant clustering of equine cases at a distance of 40 km, a distance consistent with the foraging ‘footprint’ of a flying-fox roost, suggesting the latter as a biologically plausible basis for the clustering. Getis-Ord Gi* analysis identified multiple equine infection hot spots along the eastern Australia coast from far north Queensland to central New South Wales, with the largest extending for nearly 300 km from southern Queensland to northern New South Wales. Geographically weighted regression (GWR) showed the density of P. alecto and P. conspicillatus to have the strongest positive correlation with equine case locations, suggesting these species are more likely a source of infection of Hendra virus for horses than P. poliocephalus or P. scapulatus. The density of horses, climate variables and vegetation variables were not found to be a significant risk factors, but the residuals from the GWR suggest that additional unidentified risk factors exist at the property level. Further investigations and comparisons between case and control properties are needed to identify these local risk factors.

Geospatial environmental data modelling applications using remote sensing, GIS and spatial statistics

This thesis presents novel modelling applications for environmental geospatial data using remote sensing, GIS and statistical modelling techniques. The studied themes can be classified into four main themes: (i) to develop advanced geospatial databases. Paper (I) demonstrates the creation of a geospatial database for the Glanville fritillary butterfly (Melitaea cinxia) in the Åland Islands, south-western Finland; (ii) to analyse species diversity and distribution using GIS techniques. Paper (II) presents a diversity and geographical distribution analysis for Scopulini moths at a world-wide scale; (iii) to study spatiotemporal forest cover change. Paper (III) presents a study of exotic and indigenous tree cover change detection in Taita Hills Kenya using airborne imagery and GIS analysis techniques; (iv) to explore predictive modelling techniques using geospatial data. In Paper (IV) human population occurrence and abundance in the Taita Hills highlands was predicted using the generalized additive modelling (GAM) technique. Paper (V) presents techniques to enhance fire prediction and burned area estimation at a regional scale in East Caprivi Namibia. Paper (VI) compares eight state-of-the-art predictive modelling methods to improve fire prediction, burned area estimation and fire risk mapping in East Caprivi Namibia. The results in Paper (I) showed that geospatial data can be managed effectively using advanced relational database management systems. Metapopulation data for Melitaea cinxia butterfly was successfully combined with GPS-delimited habitat patch information and climatic data. Using the geospatial database, spatial analyses were successfully conducted at habitat patch level or at more coarse analysis scales. Moreover, this study showed it appears evident that at a large-scale spatially correlated weather conditions are one of the primary causes of spatially correlated changes in Melitaea cinxia population sizes. In Paper (II) spatiotemporal characteristics of Socupulini moths description, diversity and distribution were analysed at a world-wide scale and for the first time GIS techniques were used for Scopulini moth geographical distribution analysis. This study revealed that Scopulini moths have a cosmopolitan distribution. The majority of the species have been described from the low latitudes, sub-Saharan Africa being the hot spot of species diversity. However, the taxonomical effort has been uneven among biogeographical regions. Paper III showed that forest cover change can be analysed in great detail using modern airborne imagery techniques and historical aerial photographs. However, when spatiotemporal forest cover change is studied care has to be taken in co-registration and image interpretation when historical black and white aerial photography is used. In Paper (IV) human population distribution and abundance could be modelled with fairly good results using geospatial predictors and non-Gaussian predictive modelling techniques. Moreover, land cover layer is not necessary needed as a predictor because first and second-order image texture measurements derived from satellite imagery had more power to explain the variation in dwelling unit occurrence and abundance. Paper V showed that generalized linear model (GLM) is a suitable technique for fire occurrence prediction and for burned area estimation. GLM based burned area estimations were found to be more superior than the existing MODIS burned area product (MCD45A1). However, spatial autocorrelation of fires has to be taken into account when using the GLM technique for fire occurrence prediction. Paper VI showed that novel statistical predictive modelling techniques can be used to improve fire prediction, burned area estimation and fire risk mapping at a regional scale. However, some noticeable variation between different predictive modelling techniques for fire occurrence prediction and burned area estimation existed.

Local spatial structure of forest biomass and its consequences for remote sensing of carbon stocks

Advances in forest carbon mapping have the potential to greatly reduce uncertainties in the global carbon budget and to facilitate effective emissions mitigation strategies such as REDD+ (Reducing Emissions from Deforestation and Forest Degradation). Though broad-scale mapping is based primarily on remote sensing data, the accuracy of resulting forest carbon stock estimates depends critically on the quality of field measurements and calibration procedures. The mismatch in spatial scales between field inventory plots and larger pixels of current and planned remote sensing products for forest biomass mapping is of particular concern, as it has the potential to introduce errors, especially if forest biomass shows strong local spatial variation. Here, we used 30 large (8-50 ha) globally distributed permanent forest plots to quantify the spatial variability in aboveground biomass density (AGBD in Mgha(-1)) at spatial scales ranging from 5 to 250m (0.025-6.25 ha), and to evaluate the implications of this variability for calibrating remote sensing products using simulated remote sensing footprints. We found that local spatial variability in AGBD is large for standard plot sizes, averaging 46.3% for replicate 0.1 ha subplots within a single large plot, and 16.6% for 1 ha subplots. AGBD showed weak spatial autocorrelation at distances of 20-400 m, with autocorrelation higher in sites with higher topographic variability and statistically significant in half of the sites. We further show that when field calibration plots are smaller than the remote sensing pixels, the high local spatial variability in AGBD leads to a substantial ``dilution'' bias in calibration parameters, a bias that cannot be removed with standard statistical methods. Our results suggest that topography should be explicitly accounted for in future sampling strategies and that much care must be taken in designing calibration schemes if remote sensing of forest carbon is to achieve its promise.

Does spatial clustering help explaining differences in the inequality of income distribution? : evidence from Argentina

Resumen: Este artículo analiza la relación entre la agrupación espacial de la distribución del ingreso y la desigualdad en las provincias de Argentina. El objetivo de este trabajo es usar técnicas espaciales para analizar hasta que punto la agrupación espacial de la distribución del ingreso afecta la desigualdad de la distribución del ingreso en un contexto regional de Argentina. En general, la literatura de desigualdad implícitamente considera a cada región o provincia como una entidad independiente y el potencial para la observación de la interacción a través del espacio a menudo se ha ignorado. Mientras tanto, la autocorrelación espacial ocurre cuando la distribución espacial de la variable de interés exhibe un patrón sistemático. Yo computo tres medidas de autocorrelación espacial global: La I de Moran, c de Geary, y G de Getis y Ord, como grado de CLUSTERING provincial entre 1991 y 2002. La principal conclusión del trabajo es que hay evidencia que provincias con desigualdad relativamente alta (baja) tienden a ser localizadas cerca de otras provincias con alta (baja) desigualdad más a menudo de lo esperado debido al azar. Por ende cada provincia no debería ser vista como una observación independiente, como ha sido supuesto implícitamente en estudios previos sobre la desigualdad de ingresos regional.

Effect of spatial variation on areal evapotranspiration simulation in Haibei, Tibet plateau, China

Quantification of areal evapotranspiration from remote sensing data requires the determination of surface energy balance components with support of field observations. Much attention should be given to spatial resolution sensitivity to the physics of surface heterogeneity. Using the Priestley-Taylor model, we generated evapotranspiration maps at several spatial resolutions for a heterogeneous area at Haibei, and validated the evapotranspiration maps with the flux tower data. The results suggested that the mean values for all evapotranspiration maps were quite similar but their standard deviations decreased with the coarsening of spatial resolution. When the resolution transcended about 480 m, the standard deviations drastically decreased, indicating a loss of spatial structure information of the original resolution evapotranspiration map. The absolute values of relative errors of the points for evapotranspiration maps showed a fluctuant trend as spatial resolution of input parameter data layers coarsening, and the absolute value of relative errors reached minimum when pixel size of map matched up to measuring scale of eddy covariance system. Finally, based on the analyses of the semi-variogram of the original resolution evapotranspiration map and the shapes of spatial autocorrelation indices of Moran and Geary for evapotranspiration maps at different resolutions, an appropriate resolution was suggested for the areal evapotranspiration simulation in this study area.

Using choice experiments to explore the spatial distribution of willingness to pay for rural landscape improvements

We report findings from a choice experiment survey designed to estimate the economic benefits of policy measures to improve the rural landscape in the Republic of Ireland. Using a panel mixed logit specification to account for unobserved taste heterogeneity we derived individual-specific willingness-to-pay (WTP) estimates for each respondent in the sample. We subsequently investigated the spatial dependence of these estimates. Results suggest the existence of positive spatial autocorrelation for all rural landscape attributes. As a means of benefit transfer, kriging methods were employed to interpolate WTP estimates across the whole of the Republic of Ireland. The kriged WTP surfaces confirm the existence of spatial dependence and illustrate the implied spatial variation and regional disparities in WTP for all the rural landscape improvements investigated.

Spatial and temporal ecology of Scots pine ectomycorrhizas

Spatial analysis was used to explore the distribution of individual species in an ectomycorrhizal (ECM) fungal community to address: whether mycorrhizas of individual ECM fungal species were patchily distributed, and at what scale; and what the causes of this patchiness might be. Ectomycorrhizas were extracted from spatially explicit samples of the surface organic horizons of a pine plantation. The number of mycorrhizas of each ECM fungal species was recorded using morphotyping combined with internal transcribed spacer (ITS) sequencing. Semivariograms, kriging and cluster analyses were used to determine both the extent and scale of spatial autocorrelation in species abundances, potential interactions between species, and change over time. The mycorrhizas of some, but not all, ECM fungal species were patchily distributed and the size of patches differed between species. The relative abundance of individual ECM fungal species and the position of patches of ectomycorrhizas changed between years. Spatial and temporal analysis revealed a dynamic ECM fungal community with many interspecific interactions taking place, despite the homogeneity of the host community. The spatial pattern of mycorrhizas was influenced by the underlying distribution of fine roots, but local root density was in turn influenced by the presence of specific fungal species.

Plant community assembly at small scales: spatial versus environmental factors in a central European grassland

Dispersal limitation and environmental conditions are crucial drivers of plant species distribution and establishment. As these factors operate at different spatial scales, we asked: Do the environmental factors known to determine community assembly at broad scales operate at fine scales (few meters)? How much do these factors account for community variation at fine scales? In which way do biotic and abiotic interactions drive changes in species composition? We surveyed the plant community within a dry grassland along a very steep gradient of soil characteristics like pH and nutrients. We used a spatially explicit sampling design, based on three replicated macroplots of 15x15, 12x12 and 12x12 meters in extent. Soil samples were taken to quantify several soil properties (carbon, nitrogen, plant available phosphorus, pH, water content and dehydrogenase activity as a proxy for overall microbial activity). We performed variance partitioning to assess the effect of these variables on plant composition and statistically controlled for spatial autocorrelation via eigenvector mapping. We also applied null model analysis to test for non-random patterns in species co-occurrence using randomization schemes that account for patterns expected under species interactions. At a fine spatial scale, environmental factors explained 18% of variation when controlling for spatial autocorrelation in the distribution of plant species, whereas purely spatial processes accounted for 14% variation. Null model analysis showed that species spatially segregated in a non-random way and these spatial patterns could be due to a combination of environmental filtering and biotic interactions. Our grassland study suggests that environmental factors found to be directly relevant in broad scale studies are present also at small scales, but are supplemented by spatial processes and more direct interactions like competition.

The geographical structure of British bird distributions:diversity, spatial turnover and scale

1. Using data on the spatial distribution of the British avifauna, we address three basic questions about the spatial structure of assemblages: (i) Is there a relationship between species richness (alpha diversity) and spatial turnover of species (beta diversity)? (ii) Do high richness locations have fewer species in common with neighbouring areas than low richness locations?, and (iii) Are any such relationships contingent on spatial scale (resolution or quadrat area), and do they reflect the operation of a particular kind of species-area relationship (SAR)?

2. For all measures of spatial turnover, we found a negative relationship with species richness. This held across all scales, with the exception of turnover measured as beta (sim).

3. Higher richness areas were found to have more species in common with neighbouring areas.

4. The logarithmic SAR fitted better than the power SAR overall, and fitted significantly better in areas with low richness and high turnover.

5. Spatial patterns of both turnover and richness vary with scale. The finest scale richness pattern (10 km) and the coarse scale richness pattern (90 km) are statistically unrelated. The same is true of the turnover patterns.

6. With coarsening scale, locations of the most species-rich quadrats move north. This observed sensitivity of richness 'hotspot' location to spatial scale has implications for conservation biology, e.g. the location of a reserve selected on the basis of maximum richness may change considerably with reserve size or scale of analysis.

7. Average turnover measured using indices declined with coarsening scale, but the average number of species gained or lost between neighbouring quadrats was essentially scale invariant at 10-13 species, despite mean richness rising from 80 to 146 species (across an 81-fold area increase). We show that this kind of scale invariance is consistent with the logarithmic SAR.

APPLICATION OF GEOMATICS TECHNOLOGIES TO CHARACTERIZE SPATIAL VARIABILITY AT STRATUS VINEYARDS

Vineyards vary over space and time, making geomatics technologies ideally suited to study terroir. This study applied geomatics technologies - GPS, remote sensing and GIS - to characterize the spatial variability at Stratus Vineyards in the Niagara Region. The concept of spatial terroir was used to visualize, monitor and analyze the spatial and temporal variability of variables that influence grape quality. Spatial interpolation and spatial autocorrelation were used to measure the pattern demonstrated by soil moisture, leaf water potential, vine vigour, soil composition and grape composition on two Cabernet Franc blocks and one Chardonnay block. All variables demonstrated some spatial variability within and between the vineyard block and over time. Soil moisture exhibited the most significant spatial clustering and was temporally stable. Geomatics technologies provided valuable spatial information related to the natural spatial variability at Stratus Vineyards and can be used to inform and influence vineyard management decisions.

Spatial Aggregation of Forest Songbird Territories and Possible Implications for Area Sensitivity

Habitat area requirements of forest songbirds vary greatly among species, but the causes of this variation are not well understood. Large area requirements could result from advantages for certain species when settling their territories near those of conspecifics. This phenomenon would result in spatial aggregations much larger than single territories. Species that aggregate their territories could show reduced population viability in highly fragmented forests, since remnant patches may remain unoccupied if they are too small to accommodate several territories. The objectives of this study were twofold: (1) to seek evidence of territory clusters of forest birds at various spatial scales, lags of 250-550 m, before and after controlling for habitat spatial patterns; and (2) to measure the relationship between spatial autocorrelation and apparent landscape sensitivity for these species. In analyses that ignored spatial variation of vegetation within remnant forest patches, nine of the 17 species studied significantly aggregated their territories within patches. After controlling for forest vegetation, the locations of eight out of 17 species remained significantly clustered. The aggregative pattern that we observed may, thus, be indicative of a widespread phenomenon in songbird populations. Furthermore, there was a tendency for species associated with higher forest cover to be more spatially aggregated [ERRATUM].

Non-testability of equal weights spatial dependence

We show that any invariant test for spatial autocorrelation in a spatial error or spatial lag model with equal weights matrix has power equal to size. This result holds under the assumption of an elliptical distribution. Under Gaussianity, we also show that any test whose power is larger than its size for at least one point in the parameter space must be biased.

Evolution and population structure of africanized honey bees in Brazil: Evidence from spatial analysis of morphometric data

In recent years, studies based on isoenzymatic patterns of geographic variation have revealed that what is usually called the Africanized honey bee does not constitute a single population. Instead, several local populations exist with various degrees of admixture with European honey bees. In this paper, we evaluated new data on morphometric patterns of Africanized honey bees collected at 42 localities in Brazil, using univariate and multivariate (canonical) trend surface and spatial autocorrelation analyses. The clinal patterns of variation found for genetically independent characters (wing size characters and some wing venation angles) are concordant with previous studies of malate dehydrogenase (MDH) allelic frequencies and support the hypothesis that larger honey bees in southern and southeastern Brazil originated by racial admixture in the initial phases of African honey bee colonization. Geographic variation patterns of Africanized honey bee populations reflect a demic diffusion process in which European genes were gradually lost because of the higher fitness of the African gene pool in Neotropical environmental conditions.