Deriving meteorological variables across Africa for the study and control of vector-borne disease:a comparison of remote sensing and spatial interpolation of climate

This paper presents the results of an investigation into the utility of remote sensing (RS) using meteorological satellites sensors and spatial interpolation (SI) of data from meteorological stations, for the prediction of spatial variation in monthly climate across continental Africa in 1990. Information from the Advanced Very High Resolution Radiometer (AVHRR) of the National Oceanic and Atmospheric Administration's (NOAA) polar-orbiting meteorological satellites was used to estimate land surface temperature (LST) and atmospheric moisture. Cold cloud duration (CCD) data derived from the High Resolution Radiometer (HRR) onboard the European Meteorological Satellite programme's (EUMETSAT) Meteosat satellite series were also used as a RS proxy measurement of rainfall. Temperature, atmospheric moisture and rainfall surfaces were independently derived from SI of measurements from the World Meteorological Organization (WMO) member stations of Africa. These meteorological station data were then used to test the accuracy of each methodology, so that the appropriateness of the two techniques for epidemiological research could be compared. SI was a more accurate predictor of temperature, whereas RS provided a better surrogate for rainfall; both were equally accurate at predicting atmospheric moisture. The implications of these results for mapping short and long-term climate change and hence their potential for the study anti control of disease vectors are considered. Taking into account logistic and analytical problems, there were no clear conclusions regarding the optimality of either technique, but there was considerable potential for synergy.

Geographic variation in sardine population traits: implications for stock assessment

Tese dout., Ciências e Tecnologias das Pescas, Universidade do Algarve, 2007

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].

Spatial covariation of Azotobacter abundance and soil properties: A case study using the wavelet transform

The soil microflora is very heterogeneous in its spatial distribution. The origins of this heterogeneity and its significance for soil function are not well understood. A problem for understanding spatial variation better is the assumption of statistical stationarity that is made in most of the statistical methods used to assess it. These assumptions are made explicit in geostatistical methods that have been increasingly used by soil biologists in recent years. Geostatistical methods are powerful, particularly for local prediction, but they require the assumption that the variability of a property of interest is spatially uniform, which is not always plausible given what is known about the complexity of the soil microflora and the soil environment. We have used the wavelet transform, a relatively new innovation in mathematical analysis, to investigate the spatial variation of abundance of Azotobacter in the soil of a typical agricultural landscape. The wavelet transform entails no assumptions of stationarity and is well suited to the analysis of variables that show intermittent or transient features at different spatial scales. In this study, we computed cross-variograms of Azotobacter abundance with the pH, water content and loss on ignition of the soil. These revealed scale-dependent covariation in all cases. The wavelet transform also showed that the correlation of Azotobacter abundance with all three soil properties depended on spatial scale, the correlation generally increased with spatial scale and was only significantly different from zero at some scales. However, the wavelet analysis also allowed us to show how the correlation changed across the landscape. For example, at one scale Azotobacter abundance was strongly correlated with pH in part of the transect, and not with soil water content, but this was reversed elsewhere on the transect. The results show how scale-dependent variation of potentially limiting environmental factors can induce a complex spatial pattern of abundance in a soil organism. The geostatistical methods that we used here make assumptions that are not consistent with the spatial changes in the covariation of these properties that our wavelet analysis has shown. This suggests that the wavelet transform is a powerful tool for future investigation of the spatial structure and function of soil biota. (c) 2006 Elsevier Ltd. All rights reserved.

Exploring the spatial relations between soil physical properties and apparent electrical conductivity

Relations between the apparent electrical conductivity of the soil (ECa) and top- and sub-soil physical properties were examined for two arable fields in southern England (Crowmarsh Battle Farms and the Yattendon Estate). The spatial variation of ECa and the soil properties was explored geostatistically. The variogram ranges showed that ECa varied on a similar spatial scale to many of the soil physical properties in both fields. Several features in the map of kriged predictions of ECa were also evident in maps of the soil properties. In addition, the correlation coefficients showed a strong relation between ECa and several soil properties. A moving correlation analysis enabled differences in the relations between ECa and the soil properties to be examined within the fields. The results indicated that relations were inconsistent; they were stronger in some areas than others. A regression of ECa on the principal component scores of the leading components for both fields showed that the first two components accounted for a large proportion of the variance in ECa, whereas the others accounted for little or none. For Crowmarsh topsoil sand and clay, loss on ignition and volumetric water measured in the autumn had large correlations on the first component, and for Yattendon they were large for topsoil sand and clay, and autumn and spring volumetric water. The cross-variograms suggested strong coregionalization between ECa and several soil physical properties; in particular subsoil sand and silt at Crowmarsh, and subsoil sand and clay at Yattendon. The structural correlations from the linear model of coregionalization confirmed the strength of the relations between ECa and the subsoil properties. Nevertheless, no one property was consistently important for both fields. Although a map of ECa can indicate the general patterns of spatial variation in the soil, it is not a substitute for information on soil properties obtained by sampling and analysing the soil. Nevertheless, it could be used to guide further sampling. (c) 2005 Elsevier B.V. All rights reserved.

Using the variogram to explore imagery of two different spatial resolutions

The resolution of remotely sensed data is becoming increasingly fine, and there are now many sources of data with a pixel size of 1 m x 1 m. This produces huge amounts of data that have to be stored, processed and transmitted. For environmental applications this resolution possibly provides far more data than are needed: data overload. This poses the question: how much is too much? We have explored two resolutions of data-20 in pixel SPOT data and I in pixel Computerized Airborne Multispectral Imaging System (CAMIS) data from Fort A. P. Hill (Virginia, USA), using the variogram of geostatistics. For both we used the normalized difference vegetation index (NDVI). Three scales of spatial variation were identified in both the SPOT and 1 in data: there was some overlap at the intermediate spatial scales of about 150 in and of 500 m-600 in. We subsampled the I in data and scales of variation of about 30 in and of 300 in were identified consistently until the separation between pixel centroids was 15 in (or 1 in 225pixels). At this stage, spatial scales of about 100m and 600m were described, which suggested that only now was there a real difference in the amount of spatial information available from an environmental perspective. These latter were similar spatial scales to those identified from the SPOT image. We have also analysed I in CAMIS data from Fort Story (Virginia, USA) for comparison and the outcome is similar.:From these analyses it seems that a pixel size of 20m is adequate for many environmental applications, and that if more detail is required the higher resolution data could be sub-sampled to a 10m separation between pixel centroids without any serious loss of information. This reduces significantly the amount of data that needs to be stored, transmitted and analysed and has important implications for data compression.

Estimating the spatial scale of herbicide and soil interactions by nested sampling, hierarchical analysis of variance and residual maximum likelihood

An unbalanced nested sampling design was used to investigate the spatial scale of soil and herbicide interactions at the field scale. A hierarchical analysis of variance based on residual maximum likelihood (REML) was used to analyse the data and provide a first estimate of the variogram. Soil samples were taken at 108 locations at a range of separating distances in a 9 ha field to explore small and medium scale spatial variation. Soil organic matter content, pH, particle size distribution, microbial biomass and the degradation and sorption of the herbicide, isoproturon, were determined for each soil sample. A large proportion of the spatial variation in isoproturon degradation and sorption occurred at sampling intervals less than 60 m, however, the sampling design did not resolve the variation present at scales greater than this. A sampling interval of 20-25 m should ensure that the main spatial structures are identified for isoproturon degradation rate and sorption without too great a loss of information in this field.

A geostatistical analysis of soil, vegetation, and image data characterizing land surface variation

The elucidation of spatial variation in the landscape can indicate potential wildlife habitats or breeding sites for vectors, such as ticks or mosquitoes, which cause a range of diseases. Information from remotely sensed data could aid the delineation of vegetation distribution on the ground in areas where local knowledge is limited. The data from digital images are often difficult to interpret because of pixel-to-pixel variation, that is, noise, and complex variation at more than one spatial scale. Landsat Thematic Mapper Plus (ETM+) and Satellite Pour l'Observation de La Terre (SPOT) image data were analyzed for an area close to Douna in Mali, West Africa. The variograms of the normalized difference vegetation index (NDVI) from both types of image data were nested. The parameters of the nested variogram function from the Landsat ETM+ data were used to design the sampling for a ground survey of soil and vegetation data. Variograms of the soil and vegetation data showed that their variation was anisotropic and their scales of variation were similar to those of NDVI from the SPOT data. The short- and long-range components of variation in the SPOT data were filtered out separately by factorial kriging. The map of the short-range component appears to represent the patterns of vegetation and associated shallow slopes and drainage channels of the tiger bush system. The map of the long-range component also appeared to relate to broader patterns in the tiger bush and to gentle undulations in the topography. The results suggest that the types of image data analyzed in this study could be used to identify areas with more moisture in semiarid regions that could support wildlife and also be potential vector breeding sites.

Towards more precise sugar beet management based on geostatistical analysis of spatial variability within fields

Within-field variation in sugar beet yield and quality was investigated in three commercial sugar beet fields in the east of England to identify the main associated variables and to examine the possibility of predicting yield early in the season with a view to spatially variable management of sugar beet crops. Irregular grid sampling with some purposively-located nested samples was applied. It revealed the spatial variability in each sugar beet field efficiently. In geostatistical analyses, most variograms were isotropic with moderate to strong spatial dependency indicating a significant spatial variation in sugar beet yield and associated growth and environmental variables in all directions within each field. The Kriged maps showed spatial patterns of yield variability within each field and visual association with the maps of other variables. This was confirmed by redundancy analyses and Pearson correlation coefficients. The main variables associated with yield variability were soil type, organic matter, soil moisture, weed density and canopy temperature. Kriged maps of final yield variability were strongly related to that in crop canopy cover, LAI and intercepted solar radiation early in the growing season, and the yield maps of previous crops. Therefore, yield maps of previous crops together with early assessment of sugar beet growth may make an early prediction of within-field variability in sugar beet yield possible. The Broom’s Barn sugar beet model failed to account for the spatial variability in sugar yield, but the simulation was greatly improved when corrected for early canopy development cover and when the simulated yield was adjusted for weeds and plant population. Further research to optimize inputs to maximise sugar yield should target the irrigation and fertilizing of areas within fields with low canopy cover early in the season.

Spatial modeling of electrical conductivity with neural network

Many policy decisions for agricultural management in the coastal region closely depend on the extent of intrusion of sea water. In this study, Artificial Neural Network (ANN) is used to model the spatial variation of Electrical Conductivity to determine the extent of sea water intrusion in the coastal area of Brisbane, Australia. Quarterly EC data obtained from the observation (monitoring) wells located along the coast is used for training ANN architecture. The study demonstrates that ANN is able to model the spatial variation of EC with very good accuracy (even with very less training records) when some spatial information is used as one of the inputs in the network training. The results considerable improvement when compared with the network trained without the distance information.

Stable isotopes reveal inter-annual and inter-individual variation in the diet of female Australian fur seals

Understanding the temporal and spatial variation of foraging habits of apex predators is central to understanding their role in marine ecosystems and how their populations may respond to environmental variability. In the present study, stable isotope analysis (C and N) of blood was used to investigate inter-individual and inter-annual differences in the diet of adult female Australian fur seals Arctocephalus pusillus doriferus. Positive correlations were observed between red cell and plasma values for δ¹³C and δ¹⁵N (r² = 0.47 and r² = 0.66, respectively, p < 0.001 in both cases), suggesting relatively consistent individual prey choices over 3 or 4 foraging trips. Mean δ¹⁵N values (12.8 to 17.5%) confirm the species occupies the highest marine trophic niche in the region. A significant decrease in plasma δ¹⁵N values, corresponding to two-thirds of a trophic level (ca. 2%), was observed between the 1998 to 2000 and 2003 to 2005 sampling periods. This was associated with a significant decrease in adult female body condition and is consistent with a decline, previously documented by faecal analysis, of the proportion of red cod Pseudophysis bachus, barracouta Thyrsites atun and Gould's squid Nototodarus gouldi in the diet and an increase in redbait Emmelichthys nitidus. While substantial variation in δ¹⁵N was observed within each age cohort, a significant decrease was observed with age, suggesting individual specialisation for particular prey types is evident early in adulthood, but that its composition changes as females age. In addition, generalized linear models indicated body mass had a negative influence on δ¹⁵N, which may reflect larger total body oxygen stores, facilitating individuals hunting cryptic prey of lower trophic level (e.g. octopus) on the sea floor.

Life-history and spatial determinants of somatic growth dynamics in Komodo dragon populations

Somatic growth patterns represent a major component of organismal fitness and may vary among sexes and populations due to genetic and environmental processes leading to profound differences in life-history and demography. This study considered the ontogenic, sex-specific and spatial dynamics of somatic growth patterns in ten populations of the world's largest lizard the Komodo dragon (Varanus komodoensis). The growth of 400 individual Komodo dragons was measured in a capture-mark-recapture study at ten sites on four islands in eastern Indonesia, from 2002 to 2010. Generalized Additive Mixed Models (GAMMs) and information-theoretic methods were used to examine how growth rates varied with size, age and sex, and across and within islands in relation to site-specific prey availability, lizard population density and inbreeding coefficients. Growth trajectories differed significantly with size and between sexes, indicating different energy allocation tactics and overall costs associated with reproduction. This leads to disparities in maximum body sizes and longevity. Spatial variation in growth was strongly supported by a curvilinear density-dependent growth model with highest growth rates occurring at intermediate population densities. Sex-specific trade-offs in growth underpin key differences in Komodo dragon life-history including evidence for high costs of reproduction in females. Further, inverse density-dependent growth may have profound effects on individual and population level processes that influence the demography of this species.

Semiparametric modeling of the spatial distribution of occupational accident risk in the casual labor market, Piracicaba, southeast Brazil

The objective of this study was to estimate the spatial distribution of work accident risk in the informal work market in the urban zone of an industrialized city in southeast Brazil and to examine concomitant effects of age, gender, and type of occupation after controlling for spatial risk variation. The basic methodology adopted was that of a population-based case-control study with particular interest focused on the spatial location of work. Cases were all casual workers in the city suffering work accidents during a one-year period; controls were selected from the source population of casual laborers by systematic random sampling of urban homes. The spatial distribution of work accidents was estimated via a semiparametric generalized additive model with a nonparametric bidimensional spline of the geographical coordinates of cases and controls as the nonlinear spatial component, and including age, gender, and occupation as linear predictive variables in the parametric component. We analyzed 1,918 cases and 2,245 controls between 1/11/2003 and 31/10/2004 in Piracicaba, Brazil. Areas of significantly high and low accident risk were identified in relation to mean risk in the study region (p < 0.01). Work accident risk for informal workers varied significantly in the study area. Significant age, gender, and occupational group effects on accident risk were identified after correcting for this spatial variation. A good understanding of high-risk groups and high-risk regions underpins the formulation of hypotheses concerning accident causality and the development of effective public accident prevention policies.

Possible causes of spatial variability in Diatraea saccharalis Fabr. (Lep., Crambidae) and its parasitoid populations in sugarcane fields in Brazil

Based on the large data set stored in some sugarcane mills regarding the biological control programme of Diatraea saccharalis using the parasitoids Cotesia flavipes and tachinid flies, the aim of the present study was to determine whether sugarcane varieties, number of cuts and plant ages are major determinants in promoting spatial variability in D. saccharalis, C. flavipes and tachinid fly populations in sugarcane fields in the state of São Paulo (Brazil). The data set used to support this study was provided by the Sao Joao and Barra sugarcane mills. Coefficient of variation was adopted as a measure of the spatial variability of population density. Spatial variability was estimated for the total density of D. saecharalis (parasitized + unparasitized larvae), and also for D. saccharalis larvae parasitized by tachinids and by C. flavipes. Statistical analysis revealed that the spatial variation in D. saccharalis populations was influenced by the number of cuts and by plant ages in the Barra and Sao Joao Mills, respectively. Similar results were obtained for D. saccharalis larvae parasitized by tachinids; however, in the Sao Joao Mill, the spatial variability of these populations was also influenced by the different number of varieties. Finally, considering D. saccharalis larvae parasitized by C. flavipes, no significant regressions with any of the three categories were observed for both sugarcane mills. The multitrophic implications for this agricultural system are also discussed.

The role of spatial and temporal descriptors for neotropical tadpole communities in southern Brazil

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)