Spatial and temporal scale issues in determining biomass burning regimes in Bolivia and Peru

ATSR-2 active fire data from 1996 to 2000, TRMM VIRS fire counts from 1998 to 2000 and burn scars derived from SPOT VEGETATION ( the Global Burnt Area 2000 product) were mapped for Peru and Bolivia to analyse the spatial distribution of burning and its intra- and inter-annual variability. The fire season in the region mainly occurs between May and October; though some variation was found between the six broad habitat types analysed: desert, grassland, savanna, dry forest, moist forest and yungas (the forested valleys on the eastern slope of the Andes). Increased levels of burning were generally recorded in ATSR-2 and TRMM VIRS fire data in response to the 1997/1998 El Nino, but in some areas the El Nino effect was masked by the more marked influences of socio-economic change on land use and land cover. There were differences between the three global datasets: ATSR-2 under-recorded fires in ecosystems with low net primary productivities. This was because fires are set during the day in this region and, when fuel loads are low, burn out before the ATSR-2 overpass in the region which is between 02.45 h and 03.30 h. TRMM VIRS was able to detect these fires because its overpasses cover the entire diurnal range on a monthly basis. The GBA2000 product has significant errors of commission (particularly areas of shadow in the well-dissected eastern Andes) and omission (in the agricultural zone around Santa Cruz, Bolivia and in north-west Peru). Particular attention was paid to biomass burning in high-altitude grasslands, where fire is an important pastoral management technique. Fires and burn scars from Landsat Thematic Mapper (TM) and Enhanced Thematic Mapper (ETM) data for a range of years between 1987 and 2000 were mapped for areas around Parque Nacional Rio Abiseo (Peru) and Parque Nacional Carrasco (Bolivia). Burn scars mapped in the grasslands of these two areas indicate far more burning had taken place than either the fires or the burn scars derived from global datasets. Mean scar sizes are smaller and have a smaller range in size between years the in the study area in Peru (6.6-7.1 ha) than Bolivia (16.9-162.5 ha). Trends in biomass burning in the two highland areas can be explained in terms of the changing socio-economic environments and impacts of conservation. The mismatch between the spatial scale of biomass burning in the high-altitude grasslands and the sensors used to derive global fire products means that an entire component of the fire regime in the region studied is omitted, despite its importance in the farming systems on the Andes.

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.

A comparative analysis of projected impacts of climate change on river runoff from global and catchment-scale hydrological models

We present a comparative analysis of projected impacts of climate change on river runoff from two types of distributed hydrological model, a global hydrological model (GHM) and catchment-scale hydrological models (CHM). Analyses are conducted for six catchments that are global in coverage and feature strong contrasts in spatial scale as well as climatic and development conditions. These include the Liard (Canada), Mekong (SE Asia), Okavango (SW Africa), Rio Grande (Brazil), Xiangu (China) and Harper's Brook (UK). A single GHM (Mac-PDM.09) is applied to all catchments whilst different CHMs are applied for each catchment. The CHMs typically simulate water resources impacts based on a more explicit representation of catchment water resources than that available from the GHM, and the CHMs include river routing. Simulations of average annual runoff, mean monthly runoff and high (Q5) and low (Q95) monthly runoff under baseline (1961-1990) and climate change scenarios are presented. We compare the simulated runoff response of each hydrological model to (1) prescribed increases in global mean temperature from the HadCM3 climate model and (2)a prescribed increase in global-mean temperature of 2oC for seven GCMs to explore response to climate model and structural uncertainty. We find that differences in projected changes of mean annual runoff between the two types of hydrological model can be substantial for a given GCM, and they are generally larger for indicators of high and low flow. However, they are relatively small in comparison to the range of projections across the seven GCMs. Hence, for the six catchments and seven GCMs we considered, climate model structural uncertainty is greater than the uncertainty associated with the type of hydrological model applied. Moreover, shifts in the seasonal cycle of runoff with climate change are presented similarly by both hydrological models, although for some catchments the monthly timing of high and low flows differs.This implies that for studies that seek to quantify and assess the role of climate model uncertainty on catchment-scale runoff, it may be equally as feasible to apply a GHM as it is to apply a CHM, especially when climate modelling uncertainty across the range of available GCMs is as large as it currently is. Whilst the GHM is able to represent the broad climate change signal that is represented by the CHMs, we find, however, that for some catchments there are differences between GHMs and CHMs in mean annual runoff due to differences in potential evaporation estimation methods, in the representation of the seasonality of runoff, and in the magnitude of changes in extreme monthly runoff, all of which have implications for future water management issues.

Spatial variability of liquid cloud and rain: observations and microphysical effects

Spatial variability of liquid cloud water content and rainwater content is analysed from three different observational platforms: in situ measurements from research aircraft, land-based remote sensing techniques using radar and lidar, and spaceborne remote sensing from CloudSat. The variance is found to increase with spatial scale, but also depends strongly on the cloud or rain fraction regime, with overcast regions containing less variability than broken cloud fields. This variability is shown to lead to large biases, up to a factor of 4, in both the autoconversion and accretion rates estimated at a model grid scale of ≈40 km by a typical microphysical parametrization using in-cloud mean values. A parametrization for the subgrid variability of liquid cloud and rainwater content is developed, based on the observations, which varies with both the grid scale and cloud or rain fraction, and is applicable for all model grid scales. It is then shown that if this parametrization of the variability is analytically incorporated into the autoconversion and accretion rate calculations, the bias is significantly reduced.

Influence of multi-scale landscape structure on the occurrence of carnivorous mammals in a human-modified savanna, Brazil

So Paulo is the most developed state in Brazil and contains few fragments of native ecosystems, generally surrounded by intensive agriculture lands. Despite this, some areas still shelter large native animals. We aimed at understanding how medium and large carnivores use a mosaic landscape of forest/savanna and agroecosystems, and how the species respond to different landscape parameters (percentage of landcover and edge density), in a multi-scale perspective. The response variables were: species richness, carnivore frequency and frequency for the three most recorded species (Puma concolor, Chrysocyon brachyurus and Leopardus pardalis). We compared 11 competing models using Akaike`s information criterion (AIC) and assessed model support using weight of AIC. Concurrent models were combinations of landcover types (native vegetation, ""cerrado"" formations, ""cerrado"" and eucalypt plantation), landscape feature (percentage of landcover and edge density) and spatial scale. Herein, spatial scale refers to the radius around a sampling point defining a circular landscape. The scales analyzed were 250 (fine), 1,000 (medium) and 2,000 m (coarse). The shape of curves for response variables (linear, exponential and power) was also assessed. Our results indicate that species with high mobility, P. concolor and C. brachyurus, were best explained by edge density of the native vegetation at a coarse scale (2,000 m). The relationship between P. concolor and C. brachyurus frequency had a negative power-shaped response to explanatory variables. This general trend was also observed for species richness and carnivore frequency. Species richness and P. concolor frequency were also well explained by a second concurrent model: edge density of cerrado at the fine (250 m) scale. A different response was recorded for L. pardalis, as the frequency was best explained for the amount of cerrado at the fine (250 m) scale. The curve of response was linearly positive. The contrasting results (P. concolor and C. brachyurus vs L. pardalis) may be due to the much higher mobility of the two first species, in comparison with the third. Still, L. pardalis requires habitat with higher quality when compared with other two species. This study highlights the importance of considering multiple spatial scales when evaluating species responses to different habitats. An important and new finding was the prevalence of edge density over the habitat extension to explain overall carnivore distribution, a key information for planning and management of protected areas.

Spatial and temporal population interactions between the parasitoids Cotesia flavipes and Tachinidae flies: considerations on the adverse effects of biological control practice

Biological control of Diatraea saccharalis is regarded as one of the best examples of successful classical biological control in Brazil. Since the introduction of the exotic parasitoid Cotesia flavipes, the decrease of D. saccharalis infestation in sugarcane fields has been attributed to the effectiveness of this agent. Recently, the native tachinid fly parasitoids (Lydella minense and Paratheresia claripalpis) have also been implicated in the success. Here, we investigated the spatial and temporal population interactions between C. flavipes and the tachinid flies, and provide a critical analysis of the biological control practice, focusing on the undesirable effects of introductions of exotic natural enemies. To investigate these questions, a large data set comprising information from two sugarcane mills located in the state of São Paulo, Brazil (Barra and Sao Joao Mills), was analysed. Analysis of the correlation between C. flavipes and tachinid fly population densities through time revealed that such populations were inversely correlated in the Sao Joao Mill and not correlated in the Barra Mill. Logistic regressions were computed to investigate the proportion of sites occupied by the parasitoid species at both mills as a function of time. An increasing trend in the proportion of sites occupied by C. flavipes was observed, with a concomitant decrease of the sites occupied by tachinid flies. This effect was more intense in the Sao Joao Mill. Thus, there is a convincing possibility that constant releases of C. flavipes decreased the tachinid fly populations, resulting in an undesirable effect of biological control practice.

Spatial pattern of parasitism in Diatraea saccharalis Fab. (Lep., Crambidae) populations at two different spatial scales in sugarcane fields in Brazil

The biological control of Diatraea saccharalis is regarded as one of the best examples of successful classical biological control in Brazil. Since the introduction of the exotic parasitoid, Cotesia flavipes, the decrease in D. saccharalis infestation in sugarcane fields has been attributed to the effectiveness of this agent. Native Tachinidae fly parasitoids (Lydella minense and Paratheresia claripalpis) have also been implicated in the success. Quantitative data confirming the actual contribution of these agents to the control of D. saccharalis are, however, rather scant. The purpose of this study was to investigate the spatial pattern of parasitism of these parasitoids in D. saccharalis populations at two large spatial scales (fields and zones). To investigate this subject, a large data set comprising information collected from a sugarcane mill located in the state of São Paulo, Brazil (São João sugarcane mill) was analysed. When regressions between the proportion parasitism against host density were computed, the percentage of significant regressions with either a positive or a negative slope was very small at both spatial scales for both parasitoid species. Regressing the densities of tachinid-parasitized hosts against host densities per field showed that these parasitoids presented a 'moderate aggregative' response to host densities, as 53.33% of the regressions were positively significant. Cotesia flavipes was 'weakly aggregated' on host densities at the field level, because only 33.33% of the regressions were positively significant. At the zone level, neither aggregative nor spatial proportion parasitism responses were evident for either parasitoid species due to the small percentage of significant regressions computed.

Broad-scale spatial patterns of canopy cover and pond morphology affect the structure of a Neotropical amphibian metacommunity

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Metacommunity Dynamics Of Benthic Macroinvertebrates In A Large River

Metacommunity ecology focuses on the interaction between local communities and is inherently linked to dispersal as a result. Within this framework, communities are structured by a combination of in-site responses to the immediate environment (species sorting), stochasticity (patch dynamics), and connections to other communities via distance between communities and dispersal (neutrality), and source-sink dynamics (mass effects; see Chapter 1 for a detailed description of metacommunity theory, the study site, and macroinvertebrate communities found). In Chapter 2 I describe spatial scale of study and dispersal ability as both have the ability to influence the degree to which communities interact. However, little is known about how these factors influence the importance of all metacommunity dynamics. I compared dispersal mode of immature aquatic insects and dispersal ability of winged adults across multiple spatial scales in a large river. The strongest drivers of river communities were patch dynamics, followed by species sorting, then neutrality. Active dispersers during aquatic lifestages on average exhibited lower patch dynamics, higher species sorting, and significant mass effects compared to passive dispersers. Active and strong dispersers also had a scale-independent influence of neutrality, while neutrality was stronger at broader spatial scale for passive and weak dispersers. These results indicate as dispersal ability increases patch dynamics decreases, species sorting increases, and neutrality should decrease. The perceived influence of neutrality may also be dependent on spatial scale and dispersal ability. In Chapter 3 I describe how river benthic macroinvertebrate communities may influence tributary invertebrate communities via adult flight and tributaries may influence mainstem communities via immature drift. This relationship may also depend on relative mainstem and tributary size, as well as abiotic tributary influence on mainstem habitat. To investigate the interaction between a larger river and tributary I sampled mainstem benthic invertebrate communities and quantified habitat of a 7th order river (West Branch Susquehanna River) above and below a 5th order tributary confluence, as well as 0.95-3.2 km upstream in the tributary. Non-metric multidimensional scaling showed similar patterns of clustering between sampling locations for both habitat characteristics and invertebrate communities. In addition, mainstem river communities and habitat directly downstream of the tributary confluence cluster tightly together, intermediate between tributary and mid-channel river samples. In Bray-Curtis dissimilarity comparisons between tributary and mainstem river communities the furthest upstream tributary communities were least similar to river communities. Middle tributary samples were also closest by Euclidean distance to the upstream mainstem riffle and exhibited higher similarity to mid-channel samples than the furthest downstream tributary communities. My results indicate river and tributary benthic invertebrate communities may interact and likely result in direct and indirect mass effects of a tributary on the downstream mainstem community by invertebrate drift and habitat restructuring via material delivery from the tributary. I also showed likely direct effects of adult dispersal from the river and oviposition in proximal tributary locations where Euclidian, rather than river, distance may be more important in determining river-tributary interactions.

Mach Bands: multi-scale spatial filtering and co-operative coding of edges and bars

Perception of Mach bands may be explained by spatial filtering ('lateral inhibition') that can be approximated by 2nd derivative computation, and several alternative models have been proposed. To distinguish between them, we used a novel set of ‘generalised Gaussian’ images, in which the sharp ramp-plateau junction of the Mach ramp was replaced by smoother transitions. The images ranged from a slightly blurred Mach ramp to a Gaussian edge and beyond, and also included a sine-wave edge. The probability of seeing Mach Bands increased with the (relative) sharpness of the junction, but was largely independent of absolute spatial scale. These data did not fit the predictions of MIRAGE, nor 2nd derivative computation at a single fine scale. In experiment 2, observers used a cursor to mark features on the same set of images. Data on perceived position of Mach bands did not support the local energy model. Perceived width of Mach bands was poorly explained by a single-scale edge detection model, despite its previous success with Mach edges (Wallis & Georgeson, 2009, Vision Research, 49, 1886-1893). A more successful model used separate (odd and even) scale-space filtering for edges and bars, local peak detection to find candidate features, and the MAX operator to compare odd- and even-filter response maps (Georgeson, VSS 2006, Journal of Vision 6(6), 191a). Mach bands are seen when there is a local peak in the even-filter (bar) response map, AND that peak value exceeds corresponding responses in the odd-filter (edge) maps.

Application of excitation emission matrix fluorescence monitoring in the assessment of spatial and seasonal drivers of dissolved organic matter composition: Sources and physical disturbance controls

The environmental dynamics of dissolved organic matter (DOM) were characterized for a shallow, subtropical, seagrass-dominated estuarine bay, namely Florida Bay, USA. Large spatial and seasonal variations in DOM quantity and quality were assessed using dissolved organic C (DOC) measurements and spectrophotometric properties including excitation emission matrix (EEM) fluorescence with parallel factor analysis (PARAFAC). Surface water samples were collected monthly for 2 years across the bay. DOM characteristics were statistically different across the bay, and the bay was spatially characterized into four basins based on chemical characteristics of DOM as determined by EEM-PARAFAC. Differences between zones were explained based on hydrology, geomorphology, and primary productivity of the local seagrass community. In addition, potential disturbance effects from a very active hurricane season were identified. Although the overall seasonal patterns of DOM variations were not significantly affected on a bay-wide scale by this disturbance, enhanced freshwater delivery and associated P and DOM inputs (both quantity and quality) were suggested as potential drivers for the appearance of algal blooms in high impact areas. The application of EEM-PARAFAC proved to be ideally suited for studies requiring high sample throughput methods to assess spatial and temporal ecological drivers and to determine disturbance-induced impacts in aquatic ecosystems.

Linking environmental variables with regional- scale variability in ecological structure and standing stock of carbon within UK kelp forests

Kelp forests represent some of the most productive and diverse habitats on Earth. Understanding drivers of ecological patterns at large spatial scales is critical for effective management and conservation of marine habitats. We surveyed kelp forests dominated by Laminaria hyperborea (Gunnerus) Foslie 1884 across 9° latitude and >1000 km of coastline and measured a number of physical parameters at multiple scales to link ecological structure and standing stock of carbon with environmental variables. Kelp density, biomass, morphology and age were generally greater in exposed sites within regions, highlighting the importance of wave exposure in structuring L. hyperborea populations. At the regional scale, wave-exposed kelp canopies in the cooler regions (the north and west of Scotland) were greater in biomass, height and age than in warmer regions (southwest Wales and England). The range and maximal values of estimated standing stock of carbon contained within kelp forests was greater than in historical studies, suggesting that this ecosystem property may have been previously undervalued. Kelp canopy density was positively correlated with large-scale wave fetch and fine-scale water motion, whereas kelp canopy biomass and the standing stock of carbon were positively correlated with large-scale wave fetch and light levels and negatively correlated with temperature. As light availability and summer temperature were important drivers of kelp forest biomass, effective management of human activities that may affect coastal water quality is necessary to maintain ecosystem functioning, while increased temperatures related to anthropogenic climate change may impact the structure of kelp forests and the ecosystem services they provide.

Linking environmental variables with regional- scale variability in ecological structure and standing stock of carbon within UK kelp forests

Kelp forests represent some of the most productive and diverse habitats on Earth. Understanding drivers of ecological patterns at large spatial scales is critical for effective management and conservation of marine habitats. We surveyed kelp forests dominated by Laminaria hyperborea (Gunnerus) Foslie 1884 across 9° latitude and >1000 km of coastline and measured a number of physical parameters at multiple scales to link ecological structure and standing stock of carbon with environmental variables. Kelp density, biomass, morphology and age were generally greater in exposed sites within regions, highlighting the importance of wave exposure in structuring L. hyperborea populations. At the regional scale, wave-exposed kelp canopies in the cooler regions (the north and west of Scotland) were greater in biomass, height and age than in warmer regions (southwest Wales and England). The range and maximal values of estimated standing stock of carbon contained within kelp forests was greater than in historical studies, suggesting that this ecosystem property may have been previously undervalued. Kelp canopy density was positively correlated with large-scale wave fetch and fine-scale water motion, whereas kelp canopy biomass and the standing stock of carbon were positively correlated with large-scale wave fetch and light levels and negatively correlated with temperature. As light availability and summer temperature were important drivers of kelp forest biomass, effective management of human activities that may affect coastal water quality is necessary to maintain ecosystem functioning, while increased temperatures related to anthropogenic climate change may impact the structure of kelp forests and the ecosystem services they provide.

Modeling long-term variability and change of soil moisture and groundwater level - from catchment to global scale

The water stored in and flowing through the subsurface is fundamental for sustaining human activities and needs, feeding water and its constituents to surface water bodies and supporting the functioning of their ecosystems. Quantifying the changes that affect the subsurface water is crucial for our understanding of its dynamics and changes driven by climate change and other changes in the landscape, such as in land-use and water-use. It is inherently difficult to directly measure soil moisture and groundwater levels over large spatial scales and long times. Models are therefore needed to capture the soil moisture and groundwater level dynamics over such large spatiotemporal scales. This thesis develops a modeling framework that allows for long-term catchment-scale screening of soil moisture and groundwater level changes. The novelty in this development resides in an explicit link drawn between catchment-scale hydroclimatic and soil hydraulics conditions, using observed runoff data as an approximation of soil water flux and accounting for the effects of snow storage-melting dynamics on that flux. Both past and future relative changes can be assessed by use of this modeling framework, with future change projections based on common climate model outputs. By direct model-observation comparison, the thesis shows that the developed modeling framework can reproduce the temporal variability of large-scale changes in soil water storage, as obtained from the GRACE satellite product, for most of 25 large study catchments around the world. Also compared with locally measured soil water content and groundwater level in 10 U.S. catchments, the modeling approach can reasonably well reproduce relative seasonal fluctuations around long-term average values. The developed modeling framework is further used to project soil moisture changes due to expected future climate change for 81 catchments around the world. The future soil moisture changes depend on the considered radiative forcing scenario (RCP) but are overall large for the occurrence frequency of dry and wet events and the inter-annual variability of seasonal soil moisture. These changes tend to be higher for the dry events and the dry season, respectively, than for the corresponding wet quantities, indicating increased drought risk for some parts of the world.

A multitracer approach to assess the spatial contamination pattern of hake (Merluccius merluccius) in the French Mediterranean

Chemical contamination levels and stable isotope ratios provide integrated information about contaminant exposure, trophic position and also biological and environmental influences on marine organisms. By combining these approaches with otolith shape analyses, the aim of the present study was to document the spatial variability of Hg and PCB contamination of the European hake (Merluccius merluccius) in the French Mediterranean, hypothesizing that local contaminant sources, environmental conditions and biological specificities lead to site-specific contamination patterns. High Hg concentrations discriminated Corsica (average: 1.36 ± 0.80 μg g− 1 dm) from the Gulf of Lions (average values < 0.5 μg g− 1 dm), where Rhône River input caused high PCB burdens. CB 153 average concentrations ranged between 4.00 ± 0.64 and 18.39 ± 12.38 ng g− 1 dm in the Gulf of Lions, whatever the sex of the individuals, whereas the highest values in Corsica were 6.75 ± 4.22 ng g− 1 dm. Otolith shape discriminated juveniles and adults, due to their different habitats. The use of combined ecotracers was revealed as a powerful tool to discriminate between fish populations at large and small spatial scale, and to enable understanding of the environmental and biological influences on contamination patterns.