How individual movement response to habitat edges affects population persistence and spatial spread

How individual-level movement decisions in response to habitat edges influence population-level patterns of persistence and spread of a species is a major challenge in spatial ecology and conservation biology. Here, we integrate novel insights into edge behavior, based on habitat preference and movement rates, into spatially explicit growth-dispersal models. We demonstrate how crucial ecological quantities (e.g., minimal patch size, spread rate) depend critically on these individual-level decisions. In particular, we find that including edge behavior properly in these models gives qualitatively different and intuitively more reasonable results than those of some previous studies that did not consider this level of detail. Our results highlight the importance of new empirical work on individual movement response to habitat edges. © 2013 by The University of Chicago.

Contributions of demography and dispersal parameters to the spatial spread of a stage-structured insect invasion

Stage-structured models that integrate demography and dispersal can be used to identify points in the life cycle with large effects on rates of population spatial spread, information that is vital in the development of containment strategies for invasive species. Current challenges in the application of these tools include: (1) accounting for large uncertainty in model parameters, which may violate assumptions of ‘‘local’’ perturbation metrics such as sensitivities and elasticities, and (2) forecasting not only asymptotic rates of spatial spread, as is usually done, but also transient spatial dynamics in the early stages of invasion. We developed an invasion model for the Diaprepes root weevil (DRW; Diaprepes abbreviatus [Coleoptera: Curculionidae]), a generalist herbivore that has invaded citrus-growing regions of the United States. We synthesized data on DRW demography and dispersal and generated predictions for asymptotic and transient peak invasion speeds, accounting for parameter uncertainty. We quantified the contributions of each parameter toward invasion speed using a ‘‘global’’ perturbation analysis, and we contrasted parameter contributions during the transient and asymptotic phases. We found that the asymptotic invasion speed was 0.02–0.028 km/week, although the transient peak invasion speed (0.03– 0.045 km/week) was significantly greater. Both asymptotic and transient invasions speeds were most responsive to weevil dispersal distances. However, demographic parameters that had large effects on asymptotic speed (e.g., survival of early-instar larvae) had little effect on transient speed. Comparison of the global analysis with lower-level elasticities indicated that local perturbation analysis would have generated unreliable predictions for the responsiveness of invasion speed to underlying parameters. Observed range expansion in southern Florida (1992–2006) was significantly lower than the invasion speed predicted by the model. Possible causes of this mismatch include overestimation of dispersal distances, demographic rates, and spatiotemporal variation in parameter values. This study demonstrates that, when parameter uncertainty is large, as is often the case, global perturbation analyses are needed to identify which points in the life cycle should be targets of management. Our results also suggest that effective strategies for reducing spread during the asymptotic phase may have little effect during the transient phase. Includes Appendix.

A new method for the characterisation and verification of local spatial predictability for convective scale ensembles

The use of kilometre-scale ensembles in operational forecasting provides new challenges for forecast interpretation and evaluation to account for uncertainty on the convective scale. A new neighbourhood based method is presented for evaluating and characterising the local predictability variations from convective scale ensembles. Spatial scales over which ensemble forecasts agree (agreement scales, S^A) are calculated at each grid point ij, providing a map of the spatial agreement between forecasts. By comparing the average agreement scale obtained from ensemble member pairs (S^A(mm)_ij), with that between members and radar observations (S^A(mo)_ij), this approach allows the location-dependent spatial spread-skill relationship of the ensemble to be assessed. The properties of the agreement scales are demonstrated using an idealised experiment. To demonstrate the methods in an operational context the S^A(mm)_ij and S^A(mo)_ij are calculated for six convective cases run with the Met Office UK Ensemble Prediction System. The S^A(mm)_ij highlight predictability differences between cases, which can be linked to physical processes. Maps of S^A(mm)_ij are found to summarise the spatial predictability in a compact and physically meaningful manner that is useful for forecasting and for model interpretation. Comparison of S^A(mm)_ij and S^A(mo)_ij demonstrates the case-by-case and temporal variability of the spatial spread-skill, which can again be linked to physical processes.

Allee effects and population spread in patchy landscapes

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

Do habitat suitability models reliably predict the recovery areas of threatened species?

1. Identifying those areas suitable for recolonization by threatened species is essential to support efficient conservation policies. Habitat suitability models (HSM) predict species' potential distributions, but the quality of their predictions should be carefully assessed when the species-environment equilibrium assumption is violated.2. We studied the Eurasian otter Lutra lutra, whose numbers are recovering in southern Italy. To produce widely applicable results, we chose standard HSM procedures and looked for the models' capacities in predicting the suitability of a recolonization area. We used two fieldwork datasets: presence-only data, used in the Ecological Niche Factor Analyses (ENFA), and presence-absence data, used in a Generalized Linear Model (GLM). In addition to cross-validation, we independently evaluated the models with data from a recolonization event, providing presences on a previously unoccupied river.3. Three of the models successfully predicted the suitability of the recolonization area, but the GLM built with data before the recolonization disagreed with these predictions, missing the recolonized river's suitability and badly describing the otter's niche. Our results highlighted three points of relevance to modelling practices: (1) absences may prevent the models from correctly identifying areas suitable for a species spread; (2) the selection of variables may lead to randomness in the predictions; and (3) the Area Under Curve (AUC), a commonly used validation index, was not well suited to the evaluation of model quality, whereas the Boyce Index (CBI), based on presence data only, better highlighted the models' fit to the recolonization observations.4. For species with unstable spatial distributions, presence-only models may work better than presence-absence methods in making reliable predictions of suitable areas for expansion. An iterative modelling process, using new occurrences from each step of the species spread, may also help in progressively reducing errors.5. Synthesis and applications. Conservation plans depend on reliable models of the species' suitable habitats. In non-equilibrium situations, such as the case for threatened or invasive species, models could be affected negatively by the inclusion of absence data when predicting the areas of potential expansion. Presence-only methods will here provide a better basis for productive conservation management practices.

A sodar for profiling in a spatially inhomogeneous urban environment

The urban boundary layer, above the canopy, is still poorly understood. One of the challenges is obtaining data by sampling more than a few meters above the rooftops, given the spatial and temporal inhomogeneities in both horizontal and vertical. Sodars are generally useful tools for ground-based remote sensing of winds and turbulence, but rely on horizontal homogeneity (as do lidars) in building up 3-component wind vectors from sampling three or more spatially separated volumes. The time taken for sound to travel to a typical range of 200 m and back is also a limitation. A sodar of radically different design is investigated, aimed at addressing these problems. It has a single vertical transmitted sound pulse. Doppler shifted signals are received from a number of volumes around the periphery of the transmitted beam with microphones which each having tight angular sensitivity at zenith angles slightly off-vertical. The spatial spread of sampled volumes is therefore smaller. By having more receiver microphones than a conventional sodar, the effect of smaller zenith angle is offset. More rapid profiling is also possible with a single vertical transmitted beam, instead of the usual multiple beams.A prototype design is described, together with initial field measurements. It is found that the beam forming using a single dish antenna and the drift of the sound pulse downwind both give rise to reduced performance compared with expectations. It is concluded that, while the new sodar works in principle, the compromises arising in the design mean that the expected advantages have not been realized

A rede de marcação e recaptura: o caso de borboletas do gênero heliconius

(The Mark and Recapture Network: a Heliconius case study). The current pace of habitat destruction, especially in tropical landscapes, has increased the need for understanding minimum patch requirements and patch distance as tools for conserving species in forest remnants. Mark recapture and tagging studies have been instrumental in providing parameters for functional models. Because of their popularity, ease of manipulation and well known biology, butterflies have become model in studies of spatial structure. Yet, most studies on butterflies movement have focused on temperate species that live in open habitats, in which forest patches are barrier to movement. This study aimed to view and review data from mark-recapture as a network in two species of butterfly (Heliconius erato and Heliconius melpomene). A work of marking and recapture of the species was carried out in an Atlantic forest reserve located about 20km from the city of Natal (RN). Mark recapture studies were conducted in 3 weekly visits during January-February and July-August in 2007 and 2008. Captures were more common in two sections of the dirt road, with minimal collection in the forest trail. The spatial spread of captures was similar in the two species. Yet, distances between recaptures seem to be greater for Heliconius erato than for Heliconius melpomene. In addition, the erato network is more disconnected, suggesting that this specie has shorter traveling patches. Moving on to the network, both species have similar number of links (N) and unweighed vertices (L). However, melpomene has a weighed network 50% more connections than erato. These network metrics suggest that erato has more compartmentalized network and restricted movement than melpomene. Thus, erato has a larger number of disconnected components, nC, in the network, and a smaller network diameter. The frequency distribution of network connectivity for both species was better explained by a Power-law than by a random, Poissom distribution, showing that the Power-law provides a better fit than the Poisson for both species. Moreover, the Powerlaw erato is much better adjusted than in melpomene, which should be linked to the small movements that erato makes in the network

Observações hidrológicas e resultados de modelagem no espalhamento sazonal e espacial da pluma de água Amazônica

Neste estudo são apresentadas as análises espaciais e temporais dos dados sinóticos coletados durante quatro campanhas oceanográficas do programa “REVIZEE/SCORE-N” realizadas em março-abril de 1995, outubro-novembro de 1997, abril-maio de 1999 e agosto-setembro de 2001. Em complemento às análises dos dados hidrológicos, foram utilizados resultados do modelo numérico CLIPPER - Experimento ATL6 (1/6 deg.). A ênfase deste artigo é avaliar o deslocamento sazonal e espacial da pluma de água Amazônica nas regiões localizadas à sudoeste, noroeste e em frente à foz do rio Amazonas. Dados hidrológicos e resultados do modelo demonstram que durante os meses de março-maio a pluma do rio Amazonas move-se para sudoeste, ao longo da costa norte brasileiro, tendendo a diminuir o deslocamento da pluma para sudoeste durante os meses de agosto-setembro. De outubro a novembro, não foi observado ou registrado o deslocamento da pluma para sudoeste. Resultados do modelo revelam anéis gerados pela retroflexão da Corrente Norte do Brasil (CNB) em torno de 7º N, e também a influência de água proveniente do rio Amazonas em volta desses anéis durante o mês de agosto.

Serosurveillance of Schmallenberg virus in Switzerland using bulk tank milk samples.

Infections with Schmallenberg virus (SBV), a novel Orthobunyavirus transmitted by biting midges, can cause abortions and malformations of newborns and severe symptoms in adults of domestic and wild ruminants. Understanding the temporal and spatial distribution of the virus in a certain territory is important for the control and prevention of the disease. In this study, seroprevalence of antibodies against SBV and the spatial spread of the virus was investigated in Swiss dairy cattle applying a milk serology technique on bulk milk samples. The seroprevalence in cattle herds was significantly higher in December 2012 (99.5%) compared to July 2012 (19.7%). This high between-herd seroprevalence in cattle herds was observed shortly after the first detection of viral infections. Milk samples originating from farms with seropositive animals taken in December 2012 (n=209; mean 160%) revealed significantly higher S/P% ratios than samples collected in July 2012 (n=48; mean 103.6%). This finding suggests a high within-herd seroprevalence in infected herds which makes testing of bulk tank milk samples for the identification farms with past exposures to SBV a sensitive method. It suggests also that within-herd transmission followed by seroconversion still occurred between July and December. In July 2012, positive bulk tank milk samples were mainly restricted to the western part of Switzerland whereas in December 2012, all samples except one were positive. A spatial analysis revealed a separation of regions with and without positive farms in July 2012 and no spatial clustering within the regions with positive farms. In contrast to the spatial dispersion of bluetongue virus, a virus that is also transmitted by Culicoides midges, in 2008 in Switzerland, the spread of SBV occurred from the western to the eastern part of the country. The dispersed incursion of SBV took place in the western part of Switzerland and the virus spread rapidly to the remaining territory. This spatial pattern is consistent with the hypothesis that transmission by Culicoides midges was the main way of spreading.

Novel Use of Dual Anti-Inflammatory Therapy to Overcome Drug Resistance and Improve Functional Recovery Following Spinal Cord Injury

Over 1.2 million Americans are currently living with a traumatic spinal cord injury (SCI). Despite the need for effective therapies, there are currently no proven effective treatments that can improve recovery of function in SCI patients. Many therapeutic compounds have shown promise in preclinical models of SCI, but all of these have fallen short in clinical trials. P-glycoprotein (Pgp) is an active transporter expressed on capillary endothelial cell membranes at the blood-spinal cord barrier (BSCB). Pgp limits passive diffusion of blood-borne drugs into the CNS, by actively extruding drugs from the endothelial cell membrane. Pgp can become pathologically up-regulated, thus greatly impeding therapeutic drug delivery (‘multidrug resistance’). Importantly, many drugs that have been evaluated for the treatment of SCI are Pgp substrates. We hypothesized that Pgp-mediated drug resistance diminishes the delivery and efficacy of neuroprotective drugs following SCI. We observed a progressive, spatial spread of Pgp overexpression within the injured spinal cord. To assess Pgp function, we examined spinal cord uptake of systemically-delivered riluzole, a drug that is currently being evaluated in clinical trials as an SCI intervention. Blood-to-spinal cord riluzole penetration was reduced following SCI in wild-type but not Pgp-null rats, highlighting a critical role for Pgp in mediating spinal cord drug resistance after injury. Others have shown that pro-inflammatory signaling drives Pgp up-regulation in cancer and epilepsy. We have detected inflammation in both acutely- and chronically-injured spinal cord tissue. We therefore evaluated the ability of the dual COX-/5-LOX inhibitor licofelone to attenuate Pgp-mediated drug resistance following SCI. Licofelone treatment both reduced spinal cord Pgp levels and enhanced spinal cord riluzole bioavailability following SCI. Thus, we propose that licofelone may offer a new combinatorial treatment strategy to enhance spinal cord drug delivery following SCI. Additionally, we assessed the ability of licofelone, riluzole, or both to enhance recovery of locomotor function following SCI. We found that licofelone treatment conferred a significant improvement in hindlimb function that was sustained through the end of the study. In contrast, riluzole did not improve functional outcome. We therefore conclude that licofelone holds promise as a potential neuroprotective intervention for SCI.

Adaptation is not required to explain the long-term response of axons to molecular gradients

It has been suggested that growth cones navigating through the developing nervous system might display adaptation, so that their response to gradient signals is conserved over wide variations in ligand concentration. Recently however, a new chemotaxis assay that allows the effect of gradient parameters on axonal trajectories to be finely varied has revealed a decline in gradient sensitivity on either side of an optimal concentration. We show that this behavior can be quantitatively reproduced with a computational model of axonal chemotaxis that does not employ explicit adaptation. Two crucial components of this model required to reproduce the observed sensitivity are spatial and temporal averaging. These can be interpreted as corresponding, respectively, to the spatial spread of signaling effects downstream from receptor binding, and to the finite time over which these signaling effects decay. For spatial averaging, the model predicts that an effective range of roughly one-third of the extent of the growth cone is optimal for detecting small gradient signals. For temporal decay, a timescale of about 3 minutes is required for the model to reproduce the experimentally observed sensitivity.

Effects of pixel noise and stimulus structure on visual detection performance

This thesis consisted of two major parts, one determining the masking characteristics of pixel noise and the other investigating the properties of the detection filter employed by the visual system. The theoretical cut-off frequency of white pixel noise can be defined from the size of the noise pixel. The empirical cut-off frequency, i.e. the largest size of noise pixels that mimics the effect of white noise in detection, was determined by measuring contrast energy thresholds for grating stimuli in the presence of spatial noise consisting of noise pixels of various sizes and shapes. The critical i.e. minimum number of noise pixels per grating cycle needed to mimic the effect of white noise in detection was found to decrease with the bandwidth of the stimulus. The shape of the noise pixels did not have any effect on the whiteness of pixel noise as long as there was at least the minimum number of noise pixels in all spatial dimensions. Furthermore, the masking power of white pixel noise is best described when the spectral density is calculated by taking into account all the dimensions of noise pixels, i.e. width, height, and duration, even when there is random luminance only in one of these dimensions. The properties of the detection mechanism employed by the visual system were studied by measuring contrast energy thresholds for complex spatial patterns as a function of area in the presence of white pixel noise. Human detection efficiency was obtained by comparing human performance with an ideal detector. The stimuli consisted of band-pass filtered symbols, uniform and patched gratings, and point stimuli with randomised phase spectra. In agreement with the existing literature, the detection performance was found to decline with the increasing amount of detail and contour in the stimulus. A measure of image complexity was developed and successfully applied to the data. The accuracy of the detection mechanism seems to depend on the spatial structure of the stimulus and the spatial spread of contrast energy.

A rede de marcação e recaptura: o caso de borboletas do gênero heliconius

(The Mark and Recapture Network: a Heliconius case study). The current pace of habitat destruction, especially in tropical landscapes, has increased the need for understanding minimum patch requirements and patch distance as tools for conserving species in forest remnants. Mark recapture and tagging studies have been instrumental in providing parameters for functional models. Because of their popularity, ease of manipulation and well known biology, butterflies have become model in studies of spatial structure. Yet, most studies on butterflies movement have focused on temperate species that live in open habitats, in which forest patches are barrier to movement. This study aimed to view and review data from mark-recapture as a network in two species of butterfly (Heliconius erato and Heliconius melpomene). A work of marking and recapture of the species was carried out in an Atlantic forest reserve located about 20km from the city of Natal (RN). Mark recapture studies were conducted in 3 weekly visits during January-February and July-August in 2007 and 2008. Captures were more common in two sections of the dirt road, with minimal collection in the forest trail. The spatial spread of captures was similar in the two species. Yet, distances between recaptures seem to be greater for Heliconius erato than for Heliconius melpomene. In addition, the erato network is more disconnected, suggesting that this specie has shorter traveling patches. Moving on to the network, both species have similar number of links (N) and unweighed vertices (L). However, melpomene has a weighed network 50% more connections than erato. These network metrics suggest that erato has more compartmentalized network and restricted movement than melpomene. Thus, erato has a larger number of disconnected components, nC, in the network, and a smaller network diameter. The frequency distribution of network connectivity for both species was better explained by a Power-law than by a random, Poissom distribution, showing that the Power-law provides a better fit than the Poisson for both species. Moreover, the Powerlaw erato is much better adjusted than in melpomene, which should be linked to the small movements that erato makes in the network

A spatial view of ensemble spread in convection permitting ensembles

With movement toward kilometer-scale ensembles, new techniques are needed for their characterization. A new methodology is presented for detailed spatial ensemble characterization using the fractions skill score (FSS). To evaluate spatial forecast differences, the average and standard deviation are taken of the FSS calculated over all ensemble member–member pairs at different scales and lead times. These methods were found to give important information about the ensemble behavior allowing the identification of useful spatial scales, spinup times for the model, and upscale growth of errors and forecast differences. The ensemble spread was found to be highly dependent on the spatial scales considered and the threshold applied to the field. High thresholds picked out localized and intense values that gave large temporal variability in ensemble spread: local processes and undersampling dominate for these thresholds. For lower thresholds the ensemble spread increases with time as differences between the ensemble members upscale. Two convective cases were investigated based on the Met Office United Model run at 2.2-km resolution. Different ensemble types were considered: ensembles produced using the Met Office Global and Regional Ensemble Prediction System (MOGREPS) and an ensemble produced using different model physics configurations. Comparison of the MOGREPS and multiphysics ensembles demonstrated the utility of spatial ensemble evaluation techniques for assessing the impact of different perturbation strategies and the need for assessing spread at different, believable, spatial scales.

Spatial and Temporal Spread of Bovine Tuberculosis in Wild White-Tailed Deer in Michigan

In 1975, a wild white-tailed deer infected with bovine tuberculosis was shot in the northeastern Lower Peninsula, Michigan. The shooting of a second infected deer in the same area in 1994 triggered ongoing disease surveillance in the region. By 2002, bovine tuberculosis had been confirmed in 12 Michigan counties: from 449 deer; two elk; 41 non-cervid wildlife; one captive cervid facility and 28 cattle herds. We analyzed geographic spread of disease since the surveillance began and investigated factors influencing the prevalence of disease within the infected area. These analyses reveal that 78 percent of tuberculous deer came from within a 1560 km2 'core' area, within which the prevalence of apparent disease averaged 2.5 percent. Prevalence declined dramatically outside of the core and was an order of magnitude lower 30 km from its boundary. This prevalence gradient was highly significant (P<0.0001) and did not alter over the 6 year surveillance period (P= 0.98). Within the core, deer density and supplemental feeding by hunters were positively and independently correlated with tuberculosis prevalence in deer. Together, these two factors explained 55 percent of the variation in prevalence. We conclude that bovine tuberculosis was already well established in the deer population in 1994, that the infected area has not expanded significantly since that time, and that deer over-abundance and food supplementation have both contributed to ongoing transmission of disease. Managers are currently enforcing prohibitions on deer feeding in the core and are working to lower deer numbers there through increased hunting pressure.