Spatial variation in adult sex ratio across multiple scales in the invasive golden apple snail, Pomacea canaliculata

Adult sex ratio (ASR) has critical effects on behavior and life history and has implications for population demography, including the invasiveness of introduced species. ASR exhibits immense variation in nature, yet the scale dependence of this variation is rarely analyzed. In this study, using the generalized multilevel models, we investigated the variation in ASR across multiple nested spatial scales and analyzed the underlying causes for an invasive species, the golden apple snail Pomacea canaliculata. We partitioned the variance in ASR to describe the variations at different scales and then included the explanatory variables at the individual and group levels to analyze the potential causes driving the variation in ASR. We firstly determined there is a significant female-biased ASR for this species when accounting for the spatial and temporal autocorrelations of sampling. We found that, counter to nearly equal distributed variation at plot, habitat and region levels, ASR showed little variation at the town level. Temperature and precipitation at the region level were significantly positively associated with ASR, whereas the individual weight, the density characteristic, and sampling time were not significant factors influencing ASR. Our study suggests that offspring sex ratio of this species may shape the general pattern of ASR in the population level while the environmental variables at the region level translate the unbiased offspring sex ratio to the female-biased ASR. Future research should consider the implications of climate warming on the female-biased ASR of this invasive species and thus on invasion pattern.

Maize yield and rainfall on different spatial and temporal scales in Southern Brazil

This study aimed to establish relationships between maize yield and rainfall on different temporal and spatial scales, in order to provide a basis for crop monitoring and modelling. A 16-year series of maize yield and daily rainfall from 11 municipalities and micro-regions of Rio Grande do Sul State was used. Correlation and regression analyses were used to determine associations between crop yield and rainfall for the entire crop cycle, from tasseling to 30 days after, and from 5 days before tasseling to 40 days after. Close relationships between maize yield and rainfall were found, particularly during the reproductive period (45-day period comprising the flowering and grain filling). Relationships were closer on a regional scale than at smaller scales. Implications of the crop-rainfall relationships for crop modelling are discussed.

Spatial variations at different observational scales and the seasonal distributions of stream macroalgae in a Brazilian subtropical region

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

Determining spatial and temporal scales for management: lessons from whaling

Selection of the appropriate management unit is critical to the conservation of animal populations. Defining such units depends upon knowledge of population structure and upon the timescale being considered. Here, we examine the trajectory of eleven subpopulations of five species of baleen whales to investigate temporal and spatial scales in management. These subpopulations were all extirpated by commercial whaling, and no recovery or repopulation has occurred since. In these cases, time elapsed since commercial extinction ranges from four decades to almost four centuries. We propose that these subpopulations did not recover either because cultural memory of the habitat has been lost, because widespread whaling among adjacent stocks eliminated these as sources for repopulation, and/or because segregation following exploitation produced the abandonment of certain areas. Spatial scales associated with the extirpated subpopulations are frequently smaller than those typically employed in management. Overall, the evidence indicates that: (1) the time frame for management should be at most decadal in scope (i.e., <100 yr) and based on both genetic and nongenetic evidence of population substructure, and (2) at least some stocks should be defined on a smaller spatial scale than they currently are.

Identifying the spatial scale of land use that most strongly influences overall river ecosystem health score

Catchment and riparian degradation has resulted in declining ecosystem health of streams worldwide. With restoration a priority in many regions, there is an increasing interest in the scale at which land use influences stream ecosystem health. Our goal was to use a substantial data set collected as part of a monitoring program (the Southeast Queensland, Australia, Ecological Health Monitoring Program data set, collected at 116 sites over six years) to identify the spatial scale of land use, or the combination of spatial scales, that most strongly influences overall ecosystem health. In addition, we aimed to determine whether the most influential scale differed for different aspects of ecosystem health. We used linear-mixed models and a Bayesian model-averaging approach to generate models for the overall aggregated ecosystem health score and for each of the five component indicators (fish, macroinvertebrates, water quality, nutrients, and ecosystem processes) that make up the score. Dense forest close to the survey site, mid-dense forest in the hydrologically active nearstream areas of the catchment, urbanization in the riparian buffer, and tree cover at the reach scale were all significant in explaining ecosystem health, suggesting an overriding influence of forest cover, particularly close to the stream. Season and antecedent rainfall were also important explanatory variables, with some land-use variables showing significant seasonal interactions. There were also differential influences of land use for each of the component indicators. Our approach is useful given that restoring general ecosystem health is the focus of many stream restoration projects; it allowed us to predict the scale and catchment position of restoration that would result in the greatest improvement of ecosystem health in the regions streams and rivers. The models we generated suggested that good ecosystem health can be maintained in catchments where 80% of hydrologically active areas in close proximity to the stream have mid-dense forest cover and moderate health can be obtained with 60% cover.

Quantifying the state of populations and effects of disturbances at large spatio-temporal scales: The case of coral populations in the great barrier reef

This project was a step forward in applying statistical methods and models to provide new insights for more informed decision-making at large spatial scales. The model has been designed to address complicated effects of ecological processes that govern the state of populations and uncertainties inherent in large spatio-temporal datasets. Specifically, the thesis contributes to better understanding and management of the Great Barrier Reef.

Reconciling transport models across scales: The role of volume exclusion

Diffusive transport is a universal phenomenon, throughout both biological and physical sciences, and models of diffusion are routinely used to interrogate diffusion-driven processes. However, most models neglect to take into account the role of volume exclusion, which can significantly alter diffusive transport, particularly within biological systems where the diffusing particles might occupy a significant fraction of the available space. In this work we use a random walk approach to provide a means to reconcile models that incorporate crowding effects on different spatial scales. Our work demonstrates that coarse-grained models incorporating simplified descriptions of excluded volume can be used in many circumstances, but that care must be taken in pushing the coarse-graining process too far.

Spatial subdivision among assemblages of Spanish mackerel, Scomberomorus commerson (Pisces: Scombridae) across northern Australia: implications for fisheries management

Aim: This study investigated the use of stable δ13C and δ18O isotopes in the sagittal otolith carbonate of narrow-barred Spanish mackerel, Scomberomorus commerson, as indicators of population structure across Australia. Location: Samples were collected from 25 locations extending from the lower west coast of Western Australia (30°), across northern Australian waters, and to the east coast of Australia (18°) covering a coastline length of approximately 9500 km, including samples from Indonesia. Methods: The stable δ13C and δ18O isotopes in the sagittal otolith carbonate of S. commerson were analysed using standard mass spectrometric techniques. The isotope ratios across northern Australian subregions were subjected to an agglomerative hierarchical cluster analysis to define subregions. Isotope ratios within each of the subregions were compared to assess population structure across Australia. Results: Cluster analysis separated samples into four subregions: central Western Australia, north Western Australia, northern Australia and the Gulf of Carpentaria and eastern Australia. Isotope signatures for fish from a number of sampling sites from across Australia and Indonesia were significantly different, indicating population separation. No significant differences were found in otolith isotope ratios between sampling times (no temporal variation). Main conclusions: Significant differences in the isotopic signatures of S. commerson demonstrate that there is unlikely to be any substantial movement of fish among these spatially discrete adult assemblages. The lack of temporal variation among otolith isotope ratios indicates that S. commerson populations do not undergo longshore spatial shifts in distribution during their life history. The temporal persistence of spatially explicit stable isotopic signatures indicates that, at these spatial scales, the population units sampled comprise functionally distinct management units or separate ‘stocks’ for many of the purposes of fisheries management. The spatial subdivision evident among populations of S. commerson across northern and western Australia indicates that it may be advantageous to consider S. commerson population dynamics and fisheries management from a metapopulation perspective (at least at the regional level).

Spatial ecology of a specialist insect herbivore : the leaf-mining moth Tischeria ekebladella on the pedunculate oak Quercus robur

Herbivorous insects comprise a major part of terrestrial biodiversity, and their interactions with their host plants and natural enemies are of vast ecological importance. A large body of research demonstrates that the ecology and evolution of these insects may be affected by trophic interactions, by abiotic influences, and by intraspecific processes, but so far research on these individual aspects has rarely been combined. This thesis uses the leaf-mining moth Tischeria ekebladella and the pedunculate oak (Quercus robur) as a case study to assess how spatial variation in trophic interactions and the physical distribution of host trees jointly affect the distribution, dynamics and evolution of a host-specific herbivore. With respect to habitat quality, Tischeria ekebladella experiences abundant variation at several spatial scales. Most of this variation occurs at small scales notably among leaves and shoots within individual trees. While hypothetically this could cause moths to evolve an ability to select leaves and shoots of high quality, I did not find any coupling between female preference and offspring performance. Based on my studies on temporal variation in resource quality I therefore propose that unpredictable temporal changes in the relative rankings of individual resource units may render it difficult for females to predict the fate of their developing offspring. With respect to intraspecific processes, my results suggest that limited moth dispersal in relation to the spatial distribution of oak trees plays a key role in determining the regional distribution of Tischeria ekebladella. The distribution of the moth is aggregated at the landscape level, where local leaf miner populations are less likely to be present where oaks are scarce. A modelling exercise based on empirical dispersal estimates revealed that the moth population on Wattkast an island in south-western Finland is spatially structured overall, but that the relative importance of local and regional processes on tree-specific moth dynamics varies drastically across the landscape. To conclude, my work in the oak-Tischeria ekebladella system demonstrates that the local abundance and regional distribution of a herbivore may be more strongly influenced by the spatial location of host trees than by their relative quality. Hence, it reveals the importance of considering spatial context in the study of herbivorous insects, and forms a bridge between the classical fields of plant-insect interactions and spatial ecology.

Space-Time Scales of Northward Propagation of Convection during Boreal Summer

This study uses precipitation estimates from the Tropical Rainfall Measuring Mission to quantify the spatial and temporal scales of northward propagation of convection over the Indian monsoon region during boreal summer. Propagating modes of convective systems in the intraseasonal time scales such as the Madden-Julian oscillation can interact with the intertropical convergence zone and bring active and break spells of the Indian summer monsoon. Wavelet analysis was used to quantify the spatial extent (scale) and center of these propagating convective bands, as well as the time period associated with different spatial scales. Results presented here suggest that during a good monsoon year the spatial scale of this oscillation is about 30 degrees centered around 10 degrees N. During weak monsoon years, the scale of propagation decreases and the center shifts farther south closer to the equator. A strong linear relationship is obtained between the center/scale of convective wave bands and intensity of monsoon precipitation over Indian land on the interannual time scale. Moreover, the spatial scale and its center during the break monsoon were found to be similar to an overall weak monsoon year. Based on this analysis, a new index is proposed to quantify the spatial scales associated with propagating convective bands. This automated wavelet-based technique developed here can be used to study meridional propagation of convection in a large volume of datasets from observations and model simulations. The information so obtained can be related to the interannual and intraseasonal variation of Indian monsoon precipitation.

Large scale features and assessment of spatial scale correspondence between TMPA and IMD rainfall datasets over Indian landmass

Daily rainfall datasets of 10 years (1998-2007) of Tropical Rainfall Measuring Mission (TRMM) Multi-satellite Precipitation Analysis (TMPA) version 6 and India Meteorological Department (IMD) gridded rain gauge have been compared over the Indian landmass, both in large and small spatial scales. On the larger spatial scale, the pattern correlation between the two datasets on daily scales during individual years of the study period is ranging from 0.4 to 0.7. The correlation improved significantly (similar to 0.9) when the study was confined to specific wet and dry spells each of about 5-8 days. Wavelet analysis of intraseasonal oscillations (ISO) of the southwest monsoon rainfall show the percentage contribution of the major two modes (30-50 days and 10-20 days), to be ranging respectively between similar to 30-40% and 5-10% for the various years. Analysis of inter-annual variability shows the satellite data to be underestimating seasonal rainfall by similar to 110 mm during southwest monsoon and overestimating by similar to 150 mm during northeast monsoon season. At high spatio-temporal scales, viz., 1 degrees x1 degrees grid, TMPA data do not correspond to ground truth. We have proposed here a new analysis procedure to assess the minimum spatial scale at which the two datasets are compatible with each other. This has been done by studying the contribution to total seasonal rainfall from different rainfall rate windows (at 1 mm intervals) on different spatial scales (at daily time scale). The compatibility spatial scale is seen to be beyond 5 degrees x5 degrees average spatial scale over the Indian landmass. This will help to decide the usability of TMPA products, if averaged at appropriate spatial scales, for specific process studies, e.g., cloud scale, meso scale or synoptic scale.

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.

Mapping reef fish and the seascape: using acoustics and spatial modeling to guide coastal management

Reef fish distributions are patchy in time and space with some coral reef habitats supporting higher densities (i.e., aggregations) of fish than others. Identifying and quantifying fish aggregations (particularly during spawning events) are often top priorities for coastal managers. However, the rapid mapping of these aggregations using conventional survey methods (e.g., non-technical SCUBA diving and remotely operated cameras) are limited by depth, visibility and time. Acoustic sensors (i.e., splitbeam and multibeam echosounders) are not constrained by these same limitations, and were used to concurrently map and quantify the location, density and size of reef fish along with seafloor structure in two, separate locations in the U.S. Virgin Islands. Reef fish aggregations were documented along the shelf edge, an ecologically important ecotone in the region. Fish were grouped into three classes according to body size, and relationships with the benthic seascape were modeled in one area using Boosted Regression Trees. These models were validated in a second area to test their predictive performance in locations where fish have not been mapped. Models predicting the density of large fish (≥29 cm) performed well (i.e., AUC = 0.77). Water depth and standard deviation of depth were the most influential predictors at two spatial scales (100 and 300 m). Models of small (≤11 cm) and medium (12–28 cm) fish performed poorly (i.e., AUC = 0.49 to 0.68) due to the high prevalence (45–79%) of smaller fish in both locations, and the unequal prevalence of smaller fish in the training and validation areas. Integrating acoustic sensors with spatial modeling offers a new and reliable approach to rapidly identify fish aggregations and to predict the density large fish in un-surveyed locations. This integrative approach will help coastal managers to prioritize sites, and focus their limited resources on areas that may be of higher conservation value.

Small-scale spatial and temporal variability in growth and mortality of fish larvae in the subtropical northcentral Gulf of Mexico: implications for assessing recruitment success

Extensive plankton collections were taken during seven September cruises (1990–93) along the inner continental shelf of the northcentral Gulf of Mexico (GOM). Despite the high productivity and availability of food during these cruises, significant small-scale spatial variability was found in larval growth rates for both Atlantic bumper (Chloroscombrus chrysurus, Carangidae) and vermilion snapper (Rhomboplites aurorubens, Lutjanidae). The observed variability in larval growth rates was not correlated with changes in water temperature or associated with conspicuous hydrographic features and suggested the existence of less-recognizable regions where conditions for growth vary. Cruise estimates of mortality coefficients (Z) for larval Atlantic bumper (n=32,241 larvae from six cruises) and vermilion snapper (n= 2581 larvae from four cruises) ranged from 0.20 to 0.37 and 0.19 to 0.29, respectively. Even in a subtropical climate like the GOM, where larval-stage durations may be as short as two weeks, observed variability in growth rates, particularly when combined with small changes in mortality rates, can cause order-of-magnitude differences in cumulative larval survival. To what extent the observed differences in growth rates at small spatial scales are fine-scale “noise” that ultimately is smoothed by larger-scale processes is not known. Future research is needed to further characterize the small-scale variability in growth rates of larvae, particularly with regard to microzooplankton patchiness and the temporal and spatial pattern of potential predators. Small-scale spatial variability in larval growth rates may in fact be the norm, and understanding the implications of this subtle mosaic may help us to better evaluate our ability to partition the causes of recruitment variability.

The effect of spatial scale on interactions between two weevils and their parasitoid

1. The effect of spatial scale on the interactions between three hymenopteran parasitoids and their weevil hosts was investigated. The parasitoid Mesopolobus incultus (Walker) parasitised Gymnetron pascuorum Gyll.; the parasitoids Entodon sparetus (Walker) and Bracon sp. parasitised Mecinus pyraster Herbst. Both of these weevils develop inside the seedhead of Plantago lanceolata L. but occupy different niches. Seedheads were sampled annually from 162 plants at each of two experimental sites consisting of a series of habitat patches of two distinct sizes. Data were analysed from three site-years. 2. Parasitoid densities at each site-year were closely related to the abundance of their respective weevil hosts. The overall proportion of hosts parasitised was more variable for M. incultus than for E. sparetus and Bracon sp. 3. Changes in spatial scale affected the variability of parasitoid densities. For M. incultus, there was generally a greater degree of additional heterogeneity for all increases of scale; for E. sparetus, this was true only at the largest scales; for Bracon sp., all components of variance were negative. 4. The rate of parasitism was related to host density in different ways at different spatial scales. Mesopolobus incultus exhibited inverse density dependence at the finest (seedhead) scale, direct density dependence at the intermediate (plant) scale, and density independence at the large (habitat area 729 m2) scale. Entodon sparetus showed no response to variation in host density at any spatial scale. Bracon sp. showed direct density dependence only at the intermediate and largest scales. 5. Parasitoids E. sparetus and Bracon sp. seemed able to detect more than one M. pyraster individual in seedheads with multiple host occupancy; a greater incidence of conspecific parasitoids than expected emerged from such seedheads.