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.

Wavelets and the generalization of the variogram

The experimental variogram computed in the usual way by the method of moments and the Haar wavelet transform are similar in that they filter data and yield informative summaries that may be interpreted. The variogram filters out constant values; wavelets can filter variation at several spatial scales and thereby provide a richer repertoire for analysis and demand no assumptions other than that of finite variance. This paper compares the two functions, identifying that part of the Haar wavelet transform that gives it its advantages. It goes on to show that the generalized variogram of order k=1, 2, and 3 filters linear, quadratic, and cubic polynomials from the data, respectively, which correspond with more complex wavelets in Daubechies's family. The additional filter coefficients of the latter can reveal features of the data that are not evident in its usual form. Three examples in which data recorded at regular intervals on transects are analyzed illustrate the extended form of the variogram. The apparent periodicity of gilgais in Australia seems to be accentuated as filter coefficients are added, but otherwise the analysis provides no new insight. Analysis of hyerpsectral data with a strong linear trend showed that the wavelet-based variograms filtered it out. Adding filter coefficients in the analysis of the topsoil across the Jurassic scarplands of England changed the upper bound of the variogram; it then resembled the within-class variogram computed by the method of moments. To elucidate these results, we simulated several series of data to represent a random process with values fluctuating about a mean, data with long-range linear trend, data with local trend, and data with stepped transitions. The results suggest that the wavelet variogram can filter out the effects of long-range trend, but not local trend, and of transitions from one class to another, as across boundaries.

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.

Spatial analysis of the error in a model of soil nitrogen

Models of the dynamics of nitrogen in soil (soil-N) can be used to aid the fertilizer management of a crop. The predictions of soil-N models can be validated by comparison with observed data. Validation generally involves calculating non-spatial statistics of the observations and predictions, such as their means, their mean squared-difference, and their correlation. However, when the model predictions are spatially distributed across a landscape the model requires validation with spatial statistics. There are three reasons for this: (i) the model may be more or less successful at reproducing the variance of the observations at different spatial scales; (ii) the correlation of the predictions with the observations may be different at different spatial scales; (iii) the spatial pattern of model error may be informative. In this study we used a model, parameterized with spatially variable input information about the soil, to predict the mineral-N content of soil in an arable field, and compared the results with observed data. We validated the performance of the N model spatially with a linear mixed model of the observations and model predictions, estimated by residual maximum likelihood. This novel approach allowed us to describe the joint variation of the observations and predictions as: (i) independent random variation that occurred at a fine spatial scale; (ii) correlated random variation that occurred at a coarse spatial scale; (iii) systematic variation associated with a spatial trend. The linear mixed model revealed that, in general, the performance of the N model changed depending on the spatial scale of interest. At the scales associated with random variation, the N model underestimated the variance of the observations, and the predictions were correlated poorly with the observations. At the scale of the trend, the predictions and observations shared a common surface. The spatial pattern of the error of the N model suggested that the observations were affected by the local soil condition, but this was not accounted for by the N model. In summary, the N model would be well-suited to field-scale management of soil nitrogen, but suited poorly to management at finer spatial scales. This information was not apparent with a non-spatial validation. (c),2007 Elsevier B.V. All rights reserved.

Adjoint goal-based error norms for adaptive mesh ocean modelling

Flow in the world's oceans occurs at a wide range of spatial scales, from a fraction of a metre up to many thousands of kilometers. In particular, regions of intense flow are often highly localised, for example, western boundary currents, equatorial jets, overflows and convective plumes. Conventional numerical ocean models generally use static meshes. The use of dynamically-adaptive meshes has many potential advantages but needs to be guided by an error measure reflecting the underlying physics. A method of defining an error measure to guide an adaptive meshing algorithm for unstructured tetrahedral finite elements, utilizing an adjoint or goal-based method, is described here. This method is based upon a functional, encompassing important features of the flow structure. The sensitivity of this functional, with respect to the solution variables, is used as the basis from which an error measure is derived. This error measure acts to predict those areas of the domain where resolution should be changed. A barotropic wind driven gyre problem is used to demonstrate the capabilities of the method. The overall objective of this work is to develop robust error measures for use in an oceanographic context which will ensure areas of fine mesh resolution are used only where and when they are required. (c) 2006 Elsevier Ltd. All rights reserved.

A baseline appraisal of water-dependant ecosystem services, the roles they play within desakota livelihood systems and their potential sensitivity to climate change

This report forms part of a larger research programme on 'Reinterpreting the Urban-Rural Continuum', which conceptualises and investigates current knowledge and research gaps concerning 'the role that ecosystems services play in the livelihoods of the poor in regions undergoing rapid change'. The report aims to conduct a baseline appraisal of water-dependant ecosystem services, the roles they play within desakota livelihood systems and their potential sensitivity to climate change. The appraisal is conducted at three spatial scales: global, regional (four consortia areas), and meso scale (case studies within the four regions). At all three scales of analysis water resources form the interweaving theme because water provides a vital provisioning service for people, supports all other ecosystem processes and because water resources are forecast to be severely affected under climate change scenarios. This report, combined with an Endnote library of over 1100 scientific papers, provides an annotated bibliography of water-dependant ecosystem services, the roles they play within desakota livelihood systems and their potential sensitivity to climate change. After an introductory, section, Section 2 of the report defines water-related ecosystem services and how these are affected by human activities. Current knowledge and research gaps are then explored in relation to global scale climate and related hydrological changes (e.g. floods, droughts, flow regimes) (section 3). The report then discusses the impacts of climate changes on the ESPA regions, emphasising potential responses of biomes to the combined effects of climate change and human activities (particularly land use and management), and how these effects coupled with water store and flow regime manipulation by humans may affect the functioning of catchments and their ecosystem services (section 4). Finally, at the meso-scale, case studies are presented from within the ESPA regions to illustrate the close coupling of human activities and catchment performance in the context of environmental change (section 5). At the end of each section, research needs are identified and justified. These research needs are then amalgamated in section 6.

Detecting local adaptation in widespread grassland species - the importance of scale and local plant community

1 Adaptation of plant populations to local environments has been shown in many species but local adaptation is not always apparent and spatial scales of differentiation are not well known. In a reciprocal transplant experiment we tested whether: (i) three widespread grassland species are locally adapted at a European scale; (ii) detection of local adaptation depends on competition with the local plant community; and (iii) local differentiation between neighbouring populations from contrasting habitats can be stronger than differentiation at a European scale. 2 Seeds of Holcus lanatus, Lotus corniculatus and Plantago lanceolata from a Swiss, Czech and UK population were sown in a reciprocal transplant experiment at fields that exhibit environmental conditions similar to the source sites. Seedling emergence, survival, growth and reproduction were recorded for two consecutive years. 3 The effect of competition was tested by comparing individuals in weeded monocultures with plants sown together with species from the local grassland community. To compare large-scale vs. small-scale differentiation, a neighbouring population from a contrasting habitat (wet-dry contrast) was compared with the 'home' and 'foreign' populations. 4 In P. lanceolata and H. lanatus, a significant home-site advantage was detected in fitness-related traits, thus indicating local adaptation. In L. corniculatus, an overall superiority of one provenance was found. 5 The detection of local adaptation depended on competition with the local plant community. In the absence of competition the home-site advantage was underestimated in P. lanceolata and overestimated in H. lanatus. 6 A significant population differentiation between contrasting local habitats was found. In some traits, this small-scale was greater than large-scale differentiation between countries. 7 Our results indicate that local adaptation in real plant communities cannot necessarily be predicted from plants grown in weeded monocultures and that tests on the relationship between fitness and geographical distance have to account for habitat-dependent small-scale differentiation. Considering the strong small-scale differentiation, a local provenance from a different habitat may not be the best choice in ecological restoration if distant populations from a more similar habitat are available.

Design and optimisation of a large-area process-based model for annual crops

The formulation of a new process-based crop model, the general large-area model (GLAM) for annual crops is presented. The model has been designed to operate on spatial scales commensurate with those of global and regional climate models. It aims to simulate the impact of climate on crop yield. Procedures for model parameter determination and optimisation are described, and demonstrated for the prediction of groundnut (i.e. peanut; Arachis hypogaea L.) yields across India for the period 1966-1989. Optimal parameters (e.g. extinction coefficient, transpiration efficiency, rate of change of harvest index) were stable over space and time, provided the estimate of the yield technology trend was based on the full 24-year period. The model has two location-specific parameters, the planting date, and the yield gap parameter. The latter varies spatially and is determined by calibration. The optimal value varies slightly when different input data are used. The model was tested using a historical data set on a 2.5degrees x 2.5degrees grid to simulate yields. Three sites are examined in detail-grid cells from Gujarat in the west, Andhra Pradesh towards the south, and Uttar Pradesh in the north. Agreement between observed and modelled yield was variable, with correlation coefficients of 0.74, 0.42 and 0, respectively. Skill was highest where the climate signal was greatest, and correlations were comparable to or greater than correlations with seasonal mean rainfall. Yields from all 35 cells were aggregated to simulate all-India yield. The correlation coefficient between observed and simulated yields was 0.76, and the root mean square error was 8.4% of the mean yield. The model can be easily extended to any annual crop for the investigation of the impacts of climate variability (or change) on crop yield over large areas. (C) 2004 Elsevier B.V. All rights reserved.

Dynamics of QoI sensitivity in Mycosphaerella fijiensis in Costa Rica during 2000 to 2003

From 1997 onward, the strobilurin fungicide azoxystrobin was widely used in the main banana-production zone in Costa Rica against Mycosphaerella fijiensis var. difformis causing black Sigatoka of banana. By 2000, isolates of M. fijiensis with resistance to the quinolene oxidase inhibitor fungicides were common on some farms in the area. The cause was a single point mutation from glycine to alanine in the fungal target protein, cytochrome b gene. An amplification refractory mutation system Scorpion quantitative polymerase chain reaction assay was developed and used to determine the frequency of G 143A allele in samples of M. fijiensis. Two hierarchical surveys of spatial variability, in 2001 and 2002,found no significant variation in frequency on spatial scales <10 in. This allowed the frequency of G143A alleles on a farm to be estimated efficiently by averaging single samples taken at two fixed locations. The frequency of G 143A allele in bulk samples from I I farms throughout Costa Rica was determined at 2-month intervals. There was no direct relationship between the number of spray applications and the frequency of G143A on individual farms. Instead, the frequency converged toward regional averages, presumably due to the large-scale mixing of ascospores dispersed by wind. Using trap plants in an area remote from the main producing area, immigration of resistant ascospores was detected as far as 6 km away both with and against the prevailing wind.

The population genetic structure of clonal organisms generated by exponentially bounded and fat-tailed dispersal

Long distance dispersal (LDD) plays an important role in many population processes like colonization, range expansion, and epidemics. LDD of small particles like fungal spores is often a result of turbulent wind dispersal and is best described by functions with power-law behavior in the tails ("fat tailed"). The influence of fat-tailed LDD on population genetic structure is reported in this article. In computer simulations, the population structure generated by power-law dispersal with exponents in the range of -2 to -1, in distinct contrast to that generated by exponential dispersal, has a fractal structure. As the power-law exponent becomes smaller, the distribution of individual genotypes becomes more self-similar at different scales. Common statistics like G(ST) are not well suited to summarizing differences between the population genetic structures. Instead, fractal and self-similarity statistics demonstrated differences in structure arising from fat-tailed and exponential dispersal. When dispersal is fat tailed, a log-log plot of the Simpson index against distance between subpopulations has an approximately constant gradient over a large range of spatial scales. The fractal dimension D-2 is linearly inversely related to the power-law exponent, with a slope of similar to -2. In a large simulation arena, fat-tailed LDD allows colonization of the entire space by all genotypes whereas exponentially bounded dispersal eventually confines all descendants of a single clonal lineage to a relatively small area.

Seston capture by Hydropsyche siltalai and the accuracy of capture efficiency estimates

1. Suspension feeding by caseless caddisfly larvae (Trichoptera) constitutes a major pathway for energy flow, and strongly influences productivity, in streams and rivers. 2. Consideration of the impact of these animals on lotic ecosystems has been strongly influenced by a single study investigating the efficiency of particle capture of nets built by one species of hydropsychid caddisfly. 3. Using water sampling techniques at appropriate spatial scales, and taking greater consideration of local hydrodynamics than previously, we examined the size-frequency distribution of particles captured by the nets of Hydropsyche siltalai. Our results confirm that capture nets are selective in terms of particle size, and in addition suggest that this selectivity is for particles likely to provide the most energy. 4. By incorporating estimates of flow diversion around the nets of caseless caddisfly larvae, we show that capture efficiency (CE) is considerably higher than previously estimated, and conclude that more consideration of local hydrodynamics is needed to evaluate the efficiency of particle capture. 5. We use our results to postulate a mechanistic explanation for a recent example of interspecific facilitation, whereby a reduction of near-bed velocities seen in single species monocultures leads to increased capture rates and local depletion of seston within the region of reduced velocity.

The advection of high-resolution tracers by low-resolution winds

The usefulness of any simulation of atmospheric tracers using low-resolution winds relies on both the dominance of large spatial scales in the strain and time dependence that results in a cascade in tracer scales. Here, a quantitative study on the accuracy of such tracer studies is made using the contour advection technique. It is shown that, although contour stretching rates are very insensitive to the spatial truncation of the wind field, the displacement errors in filament position are sensitive. A knowledge of displacement characteristics is essential if Lagrangian simulations are to be used for the inference of airmass origin. A quantitative lower estimate is obtained for the tracer scale factor (TSF): the ratio of the smallest resolved scale in the advecting wind field to the smallest “trustworthy” scale in the tracer field. For a baroclinic wave life cycle the TSF = 6.1 ± 0.3 while for the Northern Hemisphere wintertime lower stratosphere the TSF = 5.5 ± 0.5, when using the most stringent definition of the trustworthy scale. The similarity in the TSF for the two flows is striking and an explanation is discussed in terms of the activity of potential vorticity (PV) filaments. Uncertainty in contour initialization is investigated for the stratospheric case. The effect of smoothing initial contours is to introduce a spinup time, after which wind field truncation errors take over from initialization errors (2–3 days). It is also shown that false detail from the proliferation of finescale filaments limits the useful lifetime of such contour advection simulations to 3σ−1 days, where σ is the filament thinning rate, unless filaments narrower than the trustworthy scale are removed by contour surgery. In addition, PV analysis error and diabatic effects are so strong that only PV filaments wider than 50 km are at all believable, even for very high-resolution winds. The minimum wind field resolution required to accurately simulate filaments down to the erosion scale in the stratosphere (given an initial contour) is estimated and the implications for the modeling of atmospheric chemistry are briefly discussed.

The observed relationships between wheat and climate in China

Recent changes in climate have had a measurable impact on crop yield in China. The objective of this study is to investigate how climate variability affects wheat yield in China at different spatial scales. First the response of wheat yield to the climate at the provincial level from 1978 to 1995 for China was analysed. Wheat yield variability was only correlated with climate variability in some regions of China. At the provincial level, the variability of precipitation had a negative impact on wheat yield in parts of southeast China, but the seasonal mean temperature had a negative impact on wheat yield in only a few provinces, where significant variability in precipitation explained about 23–60% of yield variability, and temperature variability accounted for 37–41% of yield variability from 1978 to 1995. The correlation between wheat yield and climate for the whole of China from 1985 to 2000 was investigated at five spatial scales using climate data. The Climate Research Unit (CRU) and National Centers for Environmental Prediction (NCEP) proportions of the grid cells with a significant yield–precipitation correlation declined progressively from 14.6% at 0.5° to 0% at 5° scale. In contrast, the proportion of grid cells significant for the yield–temperature correlation increased progressively from 1.9% at 0.5° scale to 16% at 5° scale. This indicates that the variability of precipitation has a higher association with wheat yield at small scales (0.5°, 2°/2.5°) than at larger scales (4°/5.0°); but wheat yield has a good association with temperature at all levels of aggregation. The precipitation variable at the smaller scales (0.5°, 2°/2.5°) is a dominant factor in determining inter-annual wheat yield variability more so than at the larger scales (4°/5°). We conclude that in the current climate the relationship between wheat yield and each of precipitation and temperature becomes weaker and stronger, respectively, with an increase in spatial scale.

Evaluation of satellite-based and model re-analysis rainfall estimates for Uganda

The dependence of much of Africa on rain fed agriculture leads to a high vulnerability to fluctuations in rainfall amount. Hence, accurate monitoring of near-real time rainfall is particularly useful, for example in forewarning possible crop shortfalls in drought-prone areas. Unfortunately, ground based observations are often inadequate. Rainfall estimates from satellite-based algorithms and numerical model outputs can fill this data gap, however rigorous assessment of such estimates is required. In this case, three satellite based products (NOAA-RFE 2.0, GPCP-1DD and TAMSAT) and two numerical model outputs (ERA-40 and ERA-Interim) have been evaluated for Uganda in East Africa using a network of 27 rain gauges. The study focuses on the years 2001 to 2005 and considers the main rainy season (February to June). All data sets were converted to the same temporal and spatial scales. Kriging was used for the spatial interpolation of the gauge data. All three satellite products showed similar characteristics and had a high level of skill that exceeded both model outputs. ERA-Interim had a tendency to overestimate whilst ERA-40 consistently underestimated the Ugandan rainfall.

Restructuring and rescaling water governance in mining contexts: the co-production of waterscapes in Peru

The governance of water resources is prominent in both water policy agendas and academic scholarship. Political ecologists have made important advances in reconceptualising the relationship between water and society. Yet, while they have stressed both the scalar dimensions, and the politicised nature, of water governance, analyses of its scalar politics are relatively nascent. In this paper, we consider how the increased demand for water resources by the growing mining industry in Peru reconfigures and rescales water governance. In Peru, the mining industry’s thirst for water draws in, and reshapes, social relations, technologies, institutions and discourses that operate over varying spatial and temporal scales. We develop the concept of waterscape to examine these multiple ways in water is co-produced through mining, and become embedded in changing modes and structures of water governance, often beyond the watershed scale. We argue that an examination of waterscapes avoids the limitations of thinking about water in purely material terms, structuring analysis of water issues according to traditional spatial scales and institutional hierarchies, and taking these scales and structures for granted.