Future fish distributions constrained by depth in warming seas

European continental shelf seas have experienced intense warming over the past 30 years1. In the North Sea, fish have been comprehensively monitored throughout this period and resulting data provide a unique record of changes in distribution and abundance in response to climate change2, 3. We use these data to demonstrate the remarkable power of generalized additive models (GAMs), trained on data earlier in the time series, to reliably predict trends in distribution and abundance in later years. Then, challenging process-based models that predict substantial and ongoing poleward shifts of cold-water species4, 5, we find that GAMs coupled with climate projections predict future distributions of demersal (bottom-dwelling) fish species over the next 50 years will be strongly constrained by availability of habitat of suitable depth. This will lead to pronounced changes in community structure, species interactions and commercial fisheries, unless individual acclimation or population-level evolutionary adaptations enable fish to tolerate warmer conditions or move to previously uninhabitable locations.

Spatial modeling of data with excessive zeros applied to reindeer pellet-group counts

We analyze a real data set pertaining to reindeer fecal pellet-group counts obtained from a survey conducted in a forest area in northern Sweden. In the data set, over 70% of counts are zeros, and there is high spatial correlation. We use conditionally autoregressive random effects for modeling of spatial correlation in a Poisson generalized linear mixed model (GLMM), quasi-Poisson hierarchical generalized linear model (HGLM), zero-inflated Poisson (ZIP), and hurdle models. The quasi-Poisson HGLM allows for both under- and overdispersion with excessive zeros, while the ZIP and hurdle models allow only for overdispersion. In analyzing the real data set, we see that the quasi-Poisson HGLMs can perform better than the other commonly used models, for example, ordinary Poisson HGLMs, spatial ZIP, and spatial hurdle models, and that the underdispersed Poisson HGLMs with spatial correlation fit the reindeer data best. We develop R codes for fitting these models using a unified algorithm for the HGLMs. Spatial count response with an extremely high proportion of zeros, and underdispersion can be successfully modeled using the quasi-Poisson HGLM with spatial random effects.

Efeitos dos fatores ambientais na comunidade de macrófitas em pequenos cursos lóticos perenes subtropicais

Neste estudo foi investigado como a distribuição das espécies e a produção de biomassa de macrófitas aquáticas são influenciadas pelas condições físico-químicas do ambiente. Também foi avaliado como uma espécie com maior potencial competitivo pode interferir na diversidade de espécies da comunidade macrofítica. Para tanto, em cada um dos três arroios, foram dispostos seis transecções, perpendiculares à margem. Em cada transecção foram demarcadas três unidades amostrais de 1m², nas quais foram registrados os parâmetros fitossociológicos cobertura e frequência relativas e valor de importância. A diversidade de espécies foi estimada pelo índice de Shannon, utilizando os valores de cobertura de espécies. Para determinar a biomassa das macrófitas aquáticas foram usados quadrats de 0,25m², alocados dentro da unidade amostral de 1m² usadas para quantificar os dados fitossociológicos, nos mesmos pontos onde foi feito o levantamento de cobertura da vegetação. Utilizamos como variáveis preditoras a velocidade da corrente, radiação solar incidente, coeficiente de sombreamento, vegetação ripária arbórea adjacente, nitrogênio orgânico dissolvido, carbono orgânico dissolvido e condutividade elétrica. Foram registradas 32 espécies de macrófitas aquáticas, distribuídas em 19 famílias e 28 gêneros. Conforme Análise de Correspondência Canônica (CCA), as espécies com maiores valores de biomassa foram relacionadas a unidades amostrais com alta incidência luminosa. As unidades amostrais com dominância de Pistia stratiotes apresentaram menor diversidade de espécies indicando que esta espécie, quando encontra condições que permitam sua proliferação, pode excluir espécies de menor potencial competitivo. De acordo com GLM (Generalized Linear Model), a ausência de vegetação ripária ou presente em apenas uma das margens e baixas velocidades de corrente configura-se em condições favoráveis ao estabelecimento e desenvolvimento de macrófitas aquáticas, possibilitando produção maiores valores de biomassa.

Interaction between mercury (Hg), arsenic (As) and selenium (Se) affects the activity of glutathione S-transferase in breast milk: possible relationship with fish and shellfish intake

Breast milk is regarded as an ideal source of nutrients for the growth and development of neonates, but it can also be a potential source of pollutants. Mothers can be exposed to different contaminants as a result of their lifestyle and environmental pollution. Mercury (Hg) and arsenic (As) could adversely affect the development of fetal and neonatal nervous system. Some fish and shellfish are rich in selenium (Se), an essential trace element that forms part of several enzymes related to the detoxification process, including glutathione S-transferase (GST). The goal of this study was to determine the interaction between Hg, As and Se and analyze its effect on the activity of GST in breast milk. Milk samples were collected from women between day 7 and 10 postpartum. The GST activity was determined spectrophotometrically; total Hg, As and Se concentrations were measured by atomic absorption spectrometry. To explain the possible association of Hg, As and Se concentrations with GST activity in breast milk, generalized linear models were constructed. The model explained 44% of the GST activity measured in breast milk. The GLM suggests that GST activity was positively correlated with Hg, As and Se concentrations. The activity of the enzyme was also explained by the frequency of consumption of marine fish and shellfish in the diet of the breastfeeding women.

Modeling height-diameter relationship in pinus pinaster ait. in the forest intervention zone of Lomba, NE-Portugal

In this work, the relationship between diameter at breast height (d) and total height (h) of individual-tree was modeled with the aim to establish provisory height-diameter (h-d) equations for maritime pine (Pinus pinaster Ait.) stands in the Lomba ZIF, Northeast Portugal. Using data collected locally, several local and generalized h-d equations from the literature were tested and adaptations were also considered. Model fitting was conducted by using usual nonlinear least squares (nls) methods. The best local and generalized models selected, were also tested as mixed models applying a first-order conditional expectation (FOCE) approximation procedure and maximum likelihood methods to estimate fixed and random effects. For the calibration of the mixed models and in order to be consistent with the fitting procedure, the FOCE method was also used to test different sampling designs. The results showed that the local h-d equations with two parameters performed better than the analogous models with three parameters. However a unique set of parameter values for the local model can not be used to all maritime pine stands in Lomba ZIF and thus, a generalized model including covariates from the stand, in addition to d, was necessary to obtain an adequate predictive performance. No evident superiority of the generalized mixed model in comparison to the generalized model with nonlinear least squares parameters estimates was observed. On the other hand, in the case of the local model, the predictive performance greatly improved when random effects were included. The results showed that the mixed model based in the local h-d equation selected is a viable alternative for estimating h if variables from the stand are not available. Moreover, it was observed that it is possible to obtain an adequate calibrated response using only 2 to 5 additional h-d measurements in quantile (or random) trees from the distribution of d in the plot (stand). Balancing sampling effort, accuracy and straightforwardness in practical applications, the generalized model from nls fit is recommended. Examples of applications of the selected generalized equation to the forest management are presented, namely how to use it to complete missing information from forest inventory and also showing how such an equation can be incorporated in a stand-level decision support system that aims to optimize the forest management for the maximization of wood volume production in Lomba ZIF maritime pine stands.

Tendências das aves noturnas em Portugal: metodologias e análise de um programa de monitorização baseado em colaboradores voluntários

Uma avaliação das metodologias de análise e recolha de dados aplicadas pelo Programa NOCTUAPortugal é de extrema importância para se apurar se estas são as mais indicadas em estudos de citizen science. Comparou-se os resultados de diferentes metodologias analíticas de estimação das tendências populacionais das espécies de aves noturnas durante o período de realização do Programa NOCTUA-Portugal (análise gráfica simples, modelos lineares generalizados (GLM-Poisson e GLMM), modelos aditivos generalizados (GAM-LOESS e GAM-mgcv) e software TRIM). Analisou-se a metodologia de censo de modo a avaliar o número de registos face à duração dos pontos de escuta, comparar a eficiência do ponto de deteção com outros estudos, variação das respostas ao longo da noite e efeito da época do ano, vento, nebulosidade e luminosidade da lua. Os resultados mostraram que a metodologia analítica mais indicada era o GLMM e que não era necessário realizar nenhum ajuste em particular na metodologia de censo; Trends in nocturnal birds in Portugal Methods and analysis of a volunteer-based monitoring program ABSTRACT: An evaluation of the methodologies of analysis and data collection applied by NOCTUA-Portugal Program is extremely important to determine whether these are the most suitable in citizen science studies. We compared the results of different analytical methodologies to estimate population trends of the species of nocturnal birds during the period of the NOCTUA-Portugal Program (simple graphical analysis, generalized linear models (GLM-Poisson and GLMM), generalized additive models (GAM-LOESS and GAMmgcv) and software TRIM). We analyzed the field methodology to assess the effect of point duration on the number of records, compared the point count efficiency with other sources, the variation of responses throughout the night, the effect of time of year, wind, cloud cover and moon luminosity. The results showed that the most suitable analytical methodology was the GLMM and it was not necessary to make any particular adjustment in the field methodology.

Outcome Prediction of Consciousness Disorders in the Acute Stage Based on a Complementary Motor Behavioural Tool.

INTRODUCTION: Attaining an accurate diagnosis in the acute phase for severely brain-damaged patients presenting Disorders of Consciousness (DOC) is crucial for prognostic validity; such a diagnosis determines further medical management, in terms of therapeutic choices and end-of-life decisions. However, DOC evaluation based on validated scales, such as the Revised Coma Recovery Scale (CRS-R), can lead to an underestimation of consciousness and to frequent misdiagnoses particularly in cases of cognitive motor dissociation due to other aetiologies. The purpose of this study is to determine the clinical signs that lead to a more accurate consciousness assessment allowing more reliable outcome prediction. METHODS: From the Unit of Acute Neurorehabilitation (University Hospital, Lausanne, Switzerland) between 2011 and 2014, we enrolled 33 DOC patients with a DOC diagnosis according to the CRS-R that had been established within 28 days of brain damage. The first CRS-R assessment established the initial diagnosis of Unresponsive Wakefulness Syndrome (UWS) in 20 patients and a Minimally Consciousness State (MCS) in the remaining13 patients. We clinically evaluated the patients over time using the CRS-R scale and concurrently from the beginning with complementary clinical items of a new observational Motor Behaviour Tool (MBT). Primary endpoint was outcome at unit discharge distinguishing two main classes of patients (DOC patients having emerged from DOC and those remaining in DOC) and 6 subclasses detailing the outcome of UWS and MCS patients, respectively. Based on CRS-R and MBT scores assessed separately and jointly, statistical testing was performed in the acute phase using a non-parametric Mann-Whitney U test; longitudinal CRS-R data were modelled with a Generalized Linear Model. RESULTS: Fifty-five per cent of the UWS patients and 77% of the MCS patients had emerged from DOC. First, statistical prediction of the first CRS-R scores did not permit outcome differentiation between classes; longitudinal regression modelling of the CRS-R data identified distinct outcome evolution, but not earlier than 19 days. Second, the MBT yielded a significant outcome predictability in the acute phase (p<0.02, sensitivity>0.81). Third, a statistical comparison of the CRS-R subscales weighted by MBT became significantly predictive for DOC outcome (p<0.02). DISCUSSION: The association of MBT and CRS-R scoring improves significantly the evaluation of consciousness and the predictability of outcome in the acute phase. Subtle motor behaviour assessment provides accurate insight into the amount and the content of consciousness even in the case of cognitive motor dissociation.

Genetic diversity and connectivity remain high in Holothuria polii (Delle Chiaje 1823) across a coastal lagoon-open sea environmental gradient

Coastal lagoons represent habitats with widely heterogeneous environmental conditions, particularly as regards salinity and temperature,which fluctuate in both space and time. These characteristics suggest that physical and ecological factors could contribute to the genetic divergence among populations occurring in coastal lagoon and opencoast environments. This study investigates the genetic structure of Holothuria polii at a micro-geographic scale across theMar Menor coastal lagoon and nearbymarine areas, estimating the mitochondrial DNA variation in two gene fragments, cytochrome oxidase I (COI) and 16S rRNA (16S). Dataset of mitochondrial sequences was also used to test the influence of environmental differences between coastal lagoon andmarine waters on population genetic structure. All sampled locations exhibited high levels of haplotype diversity and low values of nucleotide diversity. Both genes showed contrasting signals of genetic differentiation (non-significant differences using COI and slight differences using 16S, which could due to different mutation rates or to differential number of exclusive haplotypes. We detected an excess of recent mutations and exclusive haplotypes, which can be generated as a result of population growth. However, selective processes can be also acting on the gene markers used; highly significant generalized additive models have been obtained considering genetic data from16S gene and independent variables such as temperature and salinity.

Landscape Patterns of Permafrost Disturbance and Degradation in the Canadian High Arctic

Intensification of permafrost disturbances such as active layer detachments (ALDs) and retrogressive thaw slumps (RTS) have been observed across the circumpolar Arctic. These features are indicators of unstable conditions stemming from recent climate warming and permafrost degradation. In order to understand the processes interacting to give rise to these features, a multidisciplinary approach is required; i.e., interactions between geomorphology, hydrology, vegetation and ground thermal conditions. The goal of this research is to detect and map permafrost disturbance, predict landscape controls over disturbance and determine approaches for monitoring disturbance, all with the goal of contributing to the mitigation of permafrost hazards. Permafrost disturbance inventories were created by applying semi-automatic change detection techniques to IKONOS satellite imagery collected at the Cape Bounty Arctic Watershed Observatory (CBAWO). These methods provide a means to estimate the spatial distribution of permafrost disturbances for a given area for use as an input in susceptibility modelling. Permafrost disturbance susceptibility models were then developed using generalized additive and generalized linear models (GAM, GLM) fitted to disturbed and undisturbed locations and relevant GIS-derived predictor variables (slope, potential solar radiation, elevation). These models successfully delineated areas across the landscape that were susceptible to disturbances locally and regionally when transferred to an independent validation location. Permafrost disturbance susceptibility models are a first-order assessment of landscape susceptibility and are promising for designing land management strategies for remote permafrost regions. Additionally, geomorphic patterns associated with higher susceptibility provide important knowledge about processes associated with the initiation of disturbances. Permafrost degradation was analyzed at the CBAWO using differential interferometric synthetic aperture radar (DInSAR). Active-layer dynamics were interpreted using inter-seasonal and intra-seasonal displacement measurements and highlight the importance of hydroclimatic factors on active layer change. Collectively, these research approaches contribute to permafrost monitoring and the assessment of landscape-scale vulnerability in order to develop permafrost disturbance mitigation strategies.

Transferability of environmental favourability models in geographic space: The case of the Iberian desman (Galemys pyrenaicus) in Portugal and Spain

Transferring distribution models between different geographical areas may be problematic, as the performance of models outside their original scope is hard to predict. A modelling procedure is needed that gets the gist of the environmental descriptors of a distribution area, without either overfitting to the training data or overestimating the species’ distribution potential.We tested the transferability power of the favourability function, a generalized linear model, on the distribution of the Iberian desman (Galemys pyrenaicus) in the Iberian territories of Portugal and Spain.We also tested the effects of two of the main potential constraints on model transferability: the analysed ranges of the predictor variables, and the completeness of the species distribution data. We modelled 10 km×10km presence/absence data from Portugal and Spain separately, extrapolated each model to the other country, and compared predictions with observations. The Spanish model, despite arguably containing more false absences, showed good predictive ability in Portugal. The Portuguese model, whose predictors ranged between only a subset of the values observed in Spain, overestimated desman distribution when transferred.We discuss possible reasons for this differential model behaviour, and highlight the importance of this kind of models for prediction and conservation applications

Between and within-year effects of haying on grassland bird

In recent years, haying has extended to Iberian Mediterranean dry grasslands potentially threatening grassland birds. We evaluate the between and within-year effects of haying on grassland birds in Alentejo region, Portugal. Our main goals were: (1) to investigate variations on bird abundance and species richness in the fields hayed, with respect to past haying events occurred in a field and its surroundings and (2) to investigate the shifts in bird abundance, species richness and spatial dynamics resulting from haying a field and its surrounding area in a given year. We conducted grassland bird censuses during the breeding season through point counts from 2012 to 2015. The relationship between bird abundance/richness and past haying events was investigated using Generalized Linear Models whereas within-year effects of haying were analysed using Generalized Additive Models. Bird abundance in a field was positively related with the surface hayed in the vicinity of that field in the previous year. However, contrasting yearly effects were found for non passerines. Also, some species prefer fields with less haying events or surface hayed, whereas others occur mostly in fields frequently managed for haying. Haying a field leads, in the short term, to its abandonment by birds, and thus to a decrease in bird abundance and, for some species, to spatial concentration in surrounding fields offering suitable habitat. We conclude that within-year effects of haying have higher impact on grassland birds than between-year effects. Maintaining haying at low levels by rotating haying yearly through the different fields in each farm and using partial haying may be an adequate way to ensure an effective management of grassland bird populations.

Generalized neurofuzzy network modeling algorithms using Bezier-Bernstein polynomial functions and additive decomposition

This paper introduces a new neurofuzzy model construction algorithm for nonlinear dynamic systems based upon basis functions that are Bezier-Bernstein polynomial functions. This paper is generalized in that it copes with n-dimensional inputs by utilising an additive decomposition construction to overcome the curse of dimensionality associated with high n. This new construction algorithm also introduces univariate Bezier-Bernstein polynomial functions for the completeness of the generalized procedure. Like the B-spline expansion based neurofuzzy systems, Bezier-Bernstein polynomial function based neurofuzzy networks hold desirable properties such as nonnegativity of the basis functions, unity of support, and interpretability of basis function as fuzzy membership functions, moreover with the additional advantages of structural parsimony and Delaunay input space partition, essentially overcoming the curse of dimensionality associated with conventional fuzzy and RBF networks. This new modeling network is based on additive decomposition approach together with two separate basis function formation approaches for both univariate and bivariate Bezier-Bernstein polynomial functions used in model construction. The overall network weights are then learnt using conventional least squares methods. Numerical examples are included to demonstrate the effectiveness of this new data based modeling approach.

Thermodynamics of Ising model with infinite-range interactions by generalized canonical ensemble

In this work we present the idea of how generalized ensembles can be used to simplify the operational study of non-additive physical systems. As alternative of the usual methods of direct integration or mean-field theory, we show how the solution of the Ising model with infinite-range interactions is obtained by using a generalized canonical ensemble. We describe how the thermodynamical properties of this model in the presence of an external magnetic field are founded by simple parametric equations. Without impairing the usual interpretation, we obtain an identical critical behaviour as observed in traditional approaches.

Working covariance model selection for generalized estimating equations

We investigate methods for data-based selection of working covariance models in the analysis of correlated data with generalized estimating equations. We study two selection criteria: Gaussian pseudolikelihood and a geodesic distance based on discrepancy between model-sensitive and model-robust regression parameter covariance estimators. The Gaussian pseudolikelihood is found in simulation to be reasonably sensitive for several response distributions and noncanonical mean-variance relations for longitudinal data. Application is also made to a clinical dataset. Assessment of adequacy of both correlation and variance models for longitudinal data should be routine in applications, and we describe open-source software supporting this practice.