Spatio-temporal Bayesian network models with latent variables for revealing trophic dynamics and functional networks in fisheries ecology

Ecosystems consist of complex dynamic interactions among species and the environment, the understanding of which has implications for predicting the environmental response to changes in climate and biodiversity. However, with the recent adoption of more explorative tools, like Bayesian networks, in predictive ecology, few assumptions can be made about the data and complex, spatially varying interactions can be recovered from collected field data. In this study, we compare Bayesian network modelling approaches accounting for latent effects to reveal species dynamics for 7 geographically and temporally varied areas within the North Sea. We also apply structure learning techniques to identify functional relationships such as prey–predator between trophic groups of species that vary across space and time. We examine if the use of a general hidden variable can reflect overall changes in the trophic dynamics of each spatial system and whether the inclusion of a specific hidden variable can model unmeasured group of species. The general hidden variable appears to capture changes in the variance of different groups of species biomass. Models that include both general and specific hidden variables resulted in identifying similarity with the underlying food web dynamics and modelling spatial unmeasured effect. We predict the biomass of the trophic groups and find that predictive accuracy varies with the models' features and across the different spatial areas thus proposing a model that allows for spatial autocorrelation and two hidden variables. Our proposed model was able to produce novel insights on this ecosystem's dynamics and ecological interactions mainly because we account for the heterogeneous nature of the driving factors within each area and their changes over time. Our findings demonstrate that accounting for additional sources of variation, by combining structure learning from data and experts' knowledge in the model architecture, has the potential for gaining deeper insights into the structure and stability of ecosystems. Finally, we were able to discover meaningful functional networks that were spatially and temporally differentiated with the particular mechanisms varying from trophic associations through interactions with climate and commercial fisheries.

The effects of spatial and temporal heterogeneity on the population dynamics of four animal species in a Danish landscape

Background: Variation in carrying capacity and population return rates is generally ignored in traditional studies of population dynamics. Variation is hard to study in the field because of difficulties controlling the environment in order to obtain statistical replicates, and because of the scale and expense of experimenting on populations. There may also be ethical issues. To circumvent these problems we used detailed simulations of the simultaneous behaviours of interacting animals in an accurate facsimile of a real Danish landscape. The models incorporate as much as possible of the behaviour and ecology of skylarks Alauda arvensis, voles Microtus agrestis, a ground beetle Bembidion lampros and a linyphiid spider Erigone atra. This allows us to quantify and evaluate the importance of spatial and temporal heterogeneity on the population dynamics of the four species. Results: Both spatial and temporal heterogeneity affected the relationship between population growth rate and population density in all four species. Spatial heterogeneity accounted for 23–30% of the variance in population growth rate after accounting for the effects of density, reflecting big differences in local carrying capacity associated with the landscape features important to individual species. Temporal heterogeneity accounted for 3–13% of the variance in vole, skylark and spider, but 43% in beetles. The associated temporal variation in carrying capacity would be problematic in traditional analyses of density dependence. Return rates were less than one in all species and essentially invariant in skylarks, spiders and beetles. Return rates varied over the landscape in voles, being slower where there were larger fluctuations in local population sizes. Conclusion: Our analyses estimated the traditional parameters of carrying capacities and return rates, but these are now seen as varying continuously over the landscape depending on habitat quality and the mechanisms of density dependence. The importance of our results lies in our demonstration that the effects of spatial and temporal heterogeneity must be accounted for if we are to have accurate predictive models for use in management and conservation. This is an area which until now has lacked an adequate theoretical framework and methodology.

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.

Insights into the dynamics and sources of atmospheric lead and particulate matter in Sao Paulo, Brazil, from high temporal resolution sampling

We present the first results of a study investigating the processes that control concentrations and sources of Pb and particulate matter in the atmosphere of Sao Paulo City Brazil Aerosols were collected with high temporal resolution (3 hours) during a four-day period in July 2005 The highest Pb concentrations measured coincided with large fireworks during celebration events and associated to high traffic occurrence Our high-resolution data highlights the impact that a singular transient event can have on air quality even in a megacity Under meteorological conditions non-conducive to pollutant dispersion Pb and particulate matter concentrations accumulated during the night leading to the highest concentrations in aerosols collected early in the morning of the following day The stable isotopes of Pb suggest that emissions from traffic remain an Important source of Pb in Sao Paulo City due to the large traffic fleet despite low Pb concentrations in fuels (C) 2010 Elsevier BV All rights reserved

Disentangling social networks from spatiotemporal dynamics: the temporal structure of a dolphin society

Social networks are static illustrations of dynamic societies, within which social interactions are constantly changing. Fundamental sources of variation include ranging behaviour and temporal demographic changes. Spatiotemporal dynamics can favour or limit opportunities for individuals to interact, and then a network may not essentially represent social processes. We examined whether a social network can embed such nonsocial effects in its topology, whereby emerging modules depict spatially or temporally segregated individuals. To this end, we applied a combination of spatial, temporal and demographic analyses to a long-term study of the association patterns of Guiana dolphins, Sotalia guianensis. We found that association patterns are organized into a modular social network. Space use was unlikely to reflect these modules, since dolphins' ranging behaviour clearly overlapped. However, a temporal demographic turnover, caused by the exit/entrance of individuals (most likely emigration/immigration), defined three modules of associations occurring at different times. Although this factor could mask real social processes, we identified the temporal scale that allowed us to account for these demographic effects. By looking within this turnover period (32 months), we assessed fission-fusion dynamics of the poorly known social organization of Guiana dolphins. We highlight that spatiotemporal dynamics can strongly influence the structure of social networks. Our findings show that hypothetical social units can emerge due to the temporal opportunities for individuals to interact. Therefore, a thorough search for a satisfactory spatiotemporal scale that removes such nonsocial noise is critical when analysing a social system. (C) 2012 The Association for the Study of Animal Behaviour. Published by Elsevier Ltd. All rights reserved.

Assessment of the population structure and temporal changes in spatial dynamics and genetic characteristics of the Atlantic bluefin tuna under a fishery independent framework.

As a large and long-lived species with high economic value, restricted spawning areas and short spawning periods, the Atlantic bluefin tuna (BFT; Thunnus thynnus) is particularly susceptible to over-exploitation. Although BFT have been targeted by fisheries in the Mediterranean Sea for thousands of years, it has only been in these last decades that the exploitation rate has reached far beyond sustainable levels. An understanding of the population structure, spatial dynamics, exploitation rates and the environmental variables that affect BFT is crucial for the conservation of the species. The aims of this PhD project were 1) to assess the accuracy of larval identification methods, 2) determine the genetic structure of modern BFT populations, 3) assess the self-recruitment rate in the Gulf of Mexico and Mediterranean spawning areas, 4) estimate the immigration rate of BFT to feeding aggregations from the various spawning areas, and 5) develop tools capable of investigating the temporal stability of population structuring in the Mediterranean Sea. Several weaknesses in modern morphology-based taxonomy including demographic decline of expert taxonomists, flawed identification keys, reluctance of the taxonomic community to embrace advances in digital communications and a general scarcity of modern user-friendly materials are reviewed. Barcoding of scombrid larvae revealed important differences in the accuracy of the taxonomic identifications carried out by different ichthyoplanktologists following morphology-based methods. Using a Genotyping-by-Sequencing a panel of 95 SNPs was developed and used to characterize the population structuring of BFT and composition of adult feeding aggregations. Using novel molecular techniques, DNA was extracted from bluefin tuna vertebrae excavated from late iron age, ancient roman settlements Byzantine-era Constantinople and a 20th century collection. A second panel of 96 SNPs was developed to genotype historical and modern samples in order to elucidate changes in population structuring and allele frequencies of loci associated with selective traits.

Spatial and Temporal Melt Variability at Helheim Glacier, East Greenland, and Its Effect on Ice Dynamics

Understanding the behavior of large outlet glaciers draining the Greenland Ice Sheet is critical for assessing the impact of climate change on sea level rise. The flow of marine-terminating outlet glaciers is partly governed by calving-related processes taking place at the terminus but is also influenced by the drainage of surface runoff to the bed through moulins, cracks, and other pathways. To investigate the extent of the latter effect, we develop a distributed surface-energy-balance model for Helheim Glacier, East Greenland, to calculate surface melt and thereby estimate runoff. The model is driven by data from an automatic weather station operated on the glacier during the summers of 2007 and 2008, and calibrated with independent measurements of ablation. Modeled melt varies over the deployment period by as much as 68% relative to the mean, with melt rates approximately 77% higher on the lower reaches of the glacier trunk than on the upper glacier. We compare melt variations during the summer season to estimates of surface velocity derived from global positioning system surveys. Near the front of the glacier, there is a significant correlation (on >95% levels) between variations in runoff (estimated from surface melt) and variations in velocity, with a 1 day delay in velocity relative to melt. Although the velocity changes are small compared to accelerations previously observed following some calving events, our findings suggest that the flow speed of Helheim Glacier is sensitive to changes in runoff. The response is most significant in the heavily crevassed, fast-moving region near the calving front. The delay in the peak of the cross-correlation function implies a transit time of 12-36 h for surface runoff to reach the bed.

Revealing regional deforestation dynamics in north-eastern Madagascar - Insights from multi-temporal land cover change analysis

The north-eastern escarpment of Madagascar harbours the island’s last remaining large-scale humid forest massifs surrounded by a small-scale agricultural mosaic. There is high deforestation, commonly thought to be caused by shifting cultivation practiced by local land users to produce upland rice. However, little is known about the dynamics between forest and shifting cultivation systems at a regional level. Our study presents a first attempt to quantify changes in the extent of forest and different agricultural land cover classes, and to identify the main dynamics of land cover change for two intervals, 1995–2005 and 2005–2011. Over the 16-year study period, the speed of forest loss increased, the total area of upland rice production remained almost stable, and the area of irrigated rice fields slightly increased. While our findings seem to confirm a general trend of land use intensification, deforestation through shifting cultivation is still on the rise. Deforestation mostly affects the small forest fragments interspersed in the agricultural mosaic and is slowly leading to a homogenization of the landscape. These findings have important implications for future interventions to slow forest loss in the region, as the processes of agricultural expansion through shifting cultivation versus intensified land use cannot per se be considered mutually exclusive.

Spatio-temporal co-ordination of RhoA, Rac1 and Cdc42 activation during prototypical edge protrusion and retraction dynamics

The three canonical Rho GTPases RhoA, Rac1 and Cdc42 co-ordinate cytoskeletal dynamics. Recent studies indicate that all three Rho GTPases are activated at the leading edge of motile fibroblasts, where their activity fluctuates at subminute time and micrometer length scales. Here, we use a microfluidic chip to acutely manipulate fibroblast edge dynamics by applying pulses of platelet-derived growth factor (PDGF) or the Rho kinase inhibitor Y-27632 (which lowers contractility). This induces acute and robust membrane protrusion and retraction events, that exhibit stereotyped cytoskeletal dynamics, allowing us to fairly compare specific morphodynamic states across experiments. Using a novel Cdc42, as well as previously described, second generation RhoA and Rac1 biosensors, we observe distinct spatio-temporal signaling programs that involve all three Rho GTPases, during protrusion/retraction edge dynamics. Our results suggest that Rac1, Cdc42 and RhoA regulate different cytoskeletal and adhesion processes to fine tune the highly plastic edge protrusion/retraction dynamics that power cell motility.

Complex temporal patterns in molecular dynamics:a direct measure of the phase-space exploration by the trajectory at macroscopic time scales

Computer simulated trajectories of bulk water molecules form complex spatiotemporal structures at the picosecond time scale. This intrinsic complexity, which underlies the formation of molecular structures at longer time scales, has been quantified using a measure of statistical complexity. The method estimates the information contained in the molecular trajectory by detecting and quantifying temporal patterns present in the simulated data (velocity time series). Two types of temporal patterns are found. The first, defined by the short-time correlations corresponding to the velocity autocorrelation decay times (â‰0.1â€ps), remains asymptotically stable for time intervals longer than several tens of nanoseconds. The second is caused by previously unknown longer-time correlations (found at longer than the nanoseconds time scales) leading to a value of statistical complexity that slowly increases with time. A direct measure based on the notion of statistical complexity that describes how the trajectory explores the phase space and independent from the particular molecular signal used as the observed time series is introduced. © 2008 The American Physical Society.

Neuronal network dynamics during epileptogenesis in the medial temporal lobe

Epilepsy is one of the most common neurological disorders, a large fraction of which is resistant to pharmacotherapy. In this light, understanding the mechanisms of epilepsy and its intractable forms in particular could create new targets for pharmacotherapeutic intervention. The current project explores the dynamic changes in neuronal network function in the chronic temporal lobe epilepsy (TLE) in rat and human brain in vitro. I focused on the process of establishment of epilepsy (epileptogenesis) in the temporal lobe. Rhythmic behaviour of the hippocampal neuronal networks in healthy animals was explored using spontaneous oscillations in the gamma frequency band (SγO). The use of an improved brain slice preparation technique resulted in the natural occurence (in the absence of pharmacological stimulation) of rhythmic activity, which was then pharmacologically characterised and compared to other models of gamma oscillations (KA- and CCh-induced oscillations) using local field potential recording technique. The results showed that SγO differed from pharmacologically driven models, suggesting higher physiological relevance of SγO. Network activity was also explored in the medial entorhinal cortex (mEC), where spontaneous slow wave oscillations (SWO) were detected. To investigate the course of chronic TLE establishment, a refined Li-pilocarpine-based model of epilepsy (RISE) was developed. The model significantly reduced animal mortality and demonstrated reduced intensity, yet high morbidy with almost 70% mean success rate of developing spontaneous recurrent seizures. We used SγO to characterize changes in the hippocampal neuronal networks throughout the epileptogenesis. The results showed that the network remained largely intact, demonstrating the subtle nature of the RISE model. Despite this, a reduction in network activity was detected during the so-called latent (no seizure) period, which was hypothesized to occur due to network fragmentation and an abnormal function of kainate receptors (KAr). We therefore explored the function of KAr by challenging SγO with kainic acid (KA). The results demonstrated a remarkable decrease in KAr response during the latent period, suggesting KAr dysfunction or altered expression, which will be further investigated using a variety of electrophysiological and immunocytochemical methods. The entorhinal cortex, together with the hippocampus, is known to play an important role in the TLE. Considering this, we investigated neuronal network function of the mEC during epileptogenesis using SWO. The results demonstrated a striking difference in AMPAr function, with possible receptor upregulation or abnormal composition in the early development of epilepsy. Alterations in receptor function inevitably lead to changes in the network function, which may play an important role in the development of epilepsy. Preliminary investigations were made using slices of human brain tissue taken following surgery for intratctable epilepsy. Initial results showed that oscillogenesis could be induced in human brain slices and that such network activity was pharmacologically similar to that observed in rodent brain. Overall, our findings suggest that excitatory glutamatergic transmission is heavily involved in the process of epileptogenesis. Together with other types of receptors, KAr and AMPAr contribute to epilepsy establishment and may be the key to uncovering its mechanism.

Divisible quantum dynamics satisfies temporal Tsirelson's bound

We show that divisibility of qubit quantum processes implies temporal Tsirelson's bound. We also prove that the classical bound of the temporal Bell's inequality holds for dynamics that can be described by entanglement-breaking channels---a more general class of dynamics than that allowed by classical physics.