Spatio-temporal analysis of an invasive plant pathogen (Phytophthora ramorum) in England and Wales

Phytophthora ramorum is a damaging invasive plant pathogen and was first discovered in the UK in 2002. Spatial point analyses were applied to the occurrence of this disease in England and Wales during the period of 2003-2006 in order to assess its spatio-temporal spread. Out of the 4301 garden centres and nurseries (GCN) surveyed, there were 164, 105, 123 and 41 sites with P. ramorum in 2003, 2004, 2005 and 2006, respectively. Spatial analysis of the observed point patterns of GCN outbreaks suggested that these sites were significantly clumped within a radius of ca 60 km in 2003, but not in later years. Further analyses were conducted to determine the relationship of GCN outbreak sites over two consecutive years and thus to infer possible disease spread over time. This analysis suggested that disease spread among GCN sites was most likely to have occurred within a distance of 60 km for 2003-2004, but not for the later years. There were 35, 63, 81 and 58 sites with P. ramorum in the semi-natural environment (SNE). Analyses were carried out to assess whether infected GCN sites could act as an inoculum source of infected SNE plants or vice versa. In all years, there was a significant spatial closeness among GCN and SNE outbreak sites within a distance of 1 km. But a significant relationship over a longer distance (within 60 km) was only observed between cases in 2003 and 2004. These analyses suggest that statutory actions taken so far appear to have reduced the extent of long-distance spread of P. ramorum among garden centres and nurseries, but not the disease spread at a shorter distance between GCN and SNE sites.

Modeling data with structural and temporal correlation using lower level and higher level multilevel models

Novel imaging techniques are playing an increasingly important role in drug development, providing insight into the mechanism of action of new chemical entities. The data sets obtained by these methods can be large with complex inter-relationships, but the most appropriate statistical analysis for handling this data is often uncertain - precisely because of the exploratory nature of the way the data are collected. We present an example from a clinical trial using magnetic resonance imaging to assess changes in atherosclerotic plaques following treatment with a tool compound with established clinical benefit. We compared two specific approaches to handle the correlations due to physical location and repeated measurements: two-level and four-level multilevel models. The two methods identified similar structural variables, but higher level multilevel models had the advantage of explaining a greater proportion of variation, and the modeling assumptions appeared to be better satisfied.

Spatio-temporal surface soil heat flux estimates from satellite data: results for the AMMA experiment at the Fakara (Niger) supersite

This paper describes a method that employs Earth Observation (EO) data to calculate spatiotemporal estimates of soil heat flux, G, using a physically-based method (the Analytical Method). The method involves a harmonic analysis of land surface temperature (LST) data. It also requires an estimate of near-surface soil thermal inertia; this property depends on soil textural composition and varies as a function of soil moisture content. The EO data needed to drive the model equations, and the ground-based data required to provide verification of the method, were obtained over the Fakara domain within the African Monsoon Multidisciplinary Analysis (AMMA) program. LST estimates (3 km × 3 km, one image 15 min−1) were derived from MSG-SEVIRI data. Soil moisture estimates were obtained from ENVISAT-ASAR data, while estimates of leaf area index, LAI, (to calculate the effect of the canopy on G, largely due to radiation extinction) were obtained from SPOT-HRV images. The variation of these variables over the Fakara domain, and implications for values of G derived from them, were discussed. Results showed that this method provides reliable large-scale spatiotemporal estimates of G. Variations in G could largely be explained by the variability in the model input variables. Furthermore, it was shown that this method is relatively insensitive to model parameters related to the vegetation or soil texture. However, the strong sensitivity of thermal inertia to soil moisture content at low values of relative saturation (<0.2) means that in arid or semi-arid climates accurate estimates of surface soil moisture content are of utmost importance, if reliable estimates of G are to be obtained. This method has the potential to improve large-scale evaporation estimates, to aid land surface model prediction and to advance research that aims to explain failure in energy balance closure of meteorological field studies.

Spatio-temporal dependencies in functional connectivity in rodent cortical cultures

Models of functional connectivity in cortical cultures on multi-electrodes arrays may aid in understanding how cognitive pathways form and improve techniques that aim to interface with neuronal systems. To enable research on such models, this study uses both data- and model-driven approaches to determine what dependencies are present in and between functional connectivity networks derived from bursts of extracellularly recorded activity. Properties of excitation in bursts were analysed using correlative techniques to assess the degree of linear dependence and then two parallel techniques were used to assess functional connectivity. Three models presenting increasing levels of spatio-temporal dependency were used to capture the dynamics of individual functional connections and their consistencies were verified using surrogate data. By comparing network-wide properties between model generated networks and functional networks from data, complex interdependencies were revealed. This indicates the persistent co-activation of neuronal pathways in spontaneous bursts, as can be found in whole brain structures.

Spatio-temporal variation in predation by urban domestic cats (Felis catus) and the acceptability of possible management actions in the UK

Urban domestic cat (Felis catus) populations can attain exceedingly high densities and are not limited by natural prey availability. This has generated concerns that they may negatively affect prey populations, leading to calls for management. We enlisted cat-owners to record prey returned home to estimate patterns of predation by free-roaming pets in different localities within the town of Reading, UK and questionnaire surveys were used to quantify attitudes to different possible management strategies. Prey return rates were highly variable: only 20% of cats returned ≥4 dead prey annually. Consequently, approximately 65% of owners received no prey in a given season, but this declined to 22% after eight seasons. The estimated mean predation rate was 18.3 prey cat−1 year−1 but this varied markedly both spatially and temporally: per capita predation rates declined with increasing cat density. Comparisons with estimates of the density of six common bird prey species indicated that cats killed numbers equivalent to adult density on c. 39% of occasions. Population modeling studies suggest that such predation rates could significantly reduce the size of local bird populations for common urban species. Conversely, most urban residents did not consider cat predation to be a significant problem. Collar-mounted anti-predation devices were the only management action acceptable to the majority of urban residents (65%), but were less acceptable to cat-owners because of perceived risks to their pets; only 24% of cats were fitted with such devices. Overall, cat predation did appear to be of sufficient magnitude to affect some prey populations, although further investigation of some key aspects of cat predation is warranted. Management of the predation behavior of urban cat populations in the UK is likely to be challenging and achieving this would require considerable engagement with cat owners.

A spatio-temporal structure-based approach to drought characterisation

Drought characterisation is an intrinsically spatio-temporal problem. A limitation of previous approaches to characterisation is that they discard much of the spatio-temporal information by reducing events to a lower-order subspace. To address this, an explicit 3-dimensional (longitude, latitude, time) structure-based method is described in which drought events are defined by a spatially and temporarily coherent set of points displaying standardised precipitation below a given threshold. Geometric methods can then be used to measure similarity between individual drought structures. Groupings of these similarities provide an alternative to traditional methods for extracting recurrent space-time signals from geophysical data. The explicit consideration of structure encourages the construction of summary statistics which relate to the event geometry. Example measures considered are the event volume, centroid, and aspect ratio. The utility of a 3-dimensional approach is demonstrated by application to the analysis of European droughts (15 °W to 35°E, and 35 °N to 70°N) for the period 1901–2006. Large-scale structure is found to be abundant with 75 events identified lasting for more than 3 months and spanning at least 0.5 × 106 km2. Near-complete dissimilarity is seen between the individual drought structures, and little or no regularity is found in the time evolution of even the most spatially similar drought events. The spatial distribution of the event centroids and the time evolution of the geographic cross-sectional areas strongly suggest that large area, sustained droughts result from the combination of multiple small area (∼106 km2) short duration (∼3 months) events. The small events are not found to occur independently in space. This leads to the hypothesis that local water feedbacks play an important role in the aggregation process.

In-situ spatio-temporal measurements of the detailed azimuthal substructure of the substorm current wedge

The substorm current wedge (SCW) is a fundamental component of geomagnetic substorms. Models tend to describe the SCW as a simple line current flowing into the ionosphere towards dawn and out of the ionosphere towards dusk, linked by a westward electrojet. We use multi-spacecraft observations from perigee passes of the Cluster 1 and 4 spacecraft during a substorm on 15 Jan 2010, in conjunction with ground-based observations, to examine the spatial structuring and temporal variability of the SCW. At this time, the spacecraft travelled east-west azimuthally above the auroral region. We show that the SCW has significant azimuthal sub-structure on scales of 100~km at altitudes of 4,000-7,000~km. We identify 26 individual current sheets in the Cluster 4 data and 34 individual current sheets in the Cluster 1 data, with Cluster 1 passing through the SCW 120-240~s after Cluster 4 at 1,300-2,000~km higher altitude. Both spacecraft observed large-scale regions of net upward and downward field-aligned current, consistent with the large-scale characteristics of the SCW, although sheets of oppositely directed currents were observed within both regions. We show that the majority of these current sheets were closely aligned to a north-south direction, in contrast to the expected east-west orientation of the pre-onset aurora. Comparing our results with observations of the field-aligned current associated with bursty bulk flows (BBFs) we conclude that significant questions remain for the explanation of SCW structuring by BBF driven ``wedgelets". Our results therefore represent constraints on future modelling and theoretical frameworks on the generation of the SCW.

Spatio-temporal variations of ozone and nitrogen dioxide concentrations under urban trees and in a nearby open area

The evergreen Quercus ilex L. is one of the most common trees in Italian urban environments and is considered effective in the uptake of particulate and gaseous atmospheric pollutants. However, the few available estimates on O3 and NO2 removal by urban Q. ilex originate from model-based studies (which indicate NO2/O3 removal capacity of Q. ilex) and not from direct measurements of air pollutant concentrations. Thus, in the urban area of Siena (central Italy) we began long-term monitoring of O3/NO2 concentrations using passive samplers at a distance of 1, 5, 10 m from a busy road, under the canopies of Q. ilex and in a nearby open-field. Measurements performed in the period June 2011-October 2013 showed always a greater decrease of NO2 concentrations under the Q. ilex canopy than in the open-field transect. Conversely, a decrease of average O3 concentrations under the tree canopy was found only in autumn after the typical Mediterranean post-summer rainfalls. Our results indicate that interactions between O3/NO2 concentrations and trees in Mediterranean urban ecosystems are affected by temporal variations in climatic conditions. We argue therefore that the direct measurement of atmospheric pollutant concentrations should be chosen to describe local changes of aerial pollution.

Spatio-temporal variability of polynya dynamics and ice production in the Laptev Sea between the winters of 1979/80 and 2007/08

Polynyas in the Laptev Sea are examined with respect to recurrence and interannual wintertime ice production.We use a polynya classification method based on passive microwave satellite data to derive daily polynya area from long-term sea-ice concentrations. This provides insight into the spatial and temporal variability of open-water and thin-ice regions on the Laptev Sea Shelf. Using thermal infrared satellite data to derive an empirical thin-ice distribution within the thickness range from 0 to 20 cm, we calculate daily average surface heat loss and the resulting wintertime ice formation within the Laptev Sea polynyas between 1979 and 2008 using reanalysis data supplied by the National Centers for Environmental Prediction, USA, as atmospheric forcing. Results indicate that previous studies significantly overestimate the contribution of polynyas to the ice production in the Laptev Sea. Average wintertime ice production in polynyas amounts to approximately 55 km39 27% and is mostly determined by the polynya area, wind speed and associated large-scale circulation patterns. No trend in ice production could be detected in the period from 1979/80 to 2007/08.

A psychophysical and computational analysis of the spatio-temporal mechanisms underlying the flash-lag effect

Several accounts put forth to explain the flash-lag effect (FLE) rely mainly on either spatial or temporal mechanisms. Here we investigated the relationship between these mechanisms by psychophysical and theoretical approaches. In a first experiment we assessed the magnitudes of the FLE and temporal-order judgments performed under identical visual stimulation. The results were interpreted by means of simulations of an artificial neural network, that wits also employed to make predictions concerning the F LE. The model predicted that a spatio-temporal mislocalisation would emerge from two, continuous and abrupt-onset, moving stimuli. Additionally, a straightforward prediction of the model revealed that the magnitude of this mislocalisation should be task-dependent, increasing when the use of the abrupt-onset moving stimulus switches from a temporal marker only to both temporal and spatial markers. Our findings confirmed the model`s predictions and point to an indissoluble interplay between spatial facilitation and processing delays in the FLE.

Spatio-temporal selectivity of loss of colour and luminance contrast sensitivity with mutiple sclerosis and optic neuritis

Colour and luminance-contrast thresholds were measured in the presence of dynamic Random Luminance-contrast Masking (RLM) in individuals who had had past diagnoses of optic neuritis (ON) some of whom have progressed to a diagnosis of multiple sclerosis (MS). To explore the spatio-temporal selectivity of chromatic and luminance losses in MS/ON, thresholds were measured using three different sizes and modulation rates of the RLM displays: small checks modulating slowly, medium-sized checks with moderate modulation and large checks modulating rapidly. The colour of the chromatic stimuli used were specified in a cone-excitation space to measure relative impairments in red–green and blue–yellow chromatic channels. These observers showed chromatic thresholds along the L/(L + M) axis that were higher than those along the S-cone axis for all display sizes/modulation rates and both red-green and blue-yellow colour thresholds were higher than luminance-contrast thresholds. The principal change in thresholds with spatio-temporal changes in the display was a reduction in thresholds for L/(L + M) and S-cones with increasing check size and modulation rate. However, luminance contrast thresholds did not change with display size/rate. These results are consistent with MS/ON selectively affecting processing in colour pathways rather than in the magnocellular pathway, and that within the colour pathways neurones with opposed L- and M-cone inputs are more damaged than colour-opponent neurons with input from S-cones.