Validation and comparison of geostatistical and spline models for spatial stream networks

Scientists need appropriate spatial-statistical models to account for the unique features of stream network data. Recent advances provide a growing methodological toolbox for modelling these data, but general-purpose statistical software has only recently emerged, with little information about when to use different approaches. We implemented a simulation study to evaluate and validate geostatistical models that use continuous distances, and penalised spline models that use a finite discrete approximation for stream networks. Data were simulated from the geostatistical model, with performance measured by empirical prediction and fixed effects estimation. We found that both models were comparable in terms of squared error, with a slight advantage for the geostatistical models. Generally, both methods were unbiased and had valid confidence intervals. The most marked differences were found for confidence intervals on fixed-effect parameter estimates, where, for small sample sizes, the spline models underestimated variance. However, the penalised spline models were always more computationally efficient, which may be important for real-time prediction and estimation. Thus, decisions about which method to use must be influenced by the size and format of the data set, in addition to the characteristics of the environmental process and the modelling goals

Accounting for uncertainty in ecological analysis: the strengths and limitations of hierarchical statistical modeling

Analyses of ecological data should account for the uncertainty in the process(es) that generated the data. However, accounting for these uncertainties is a difficult task, since ecology is known for its complexity. Measurement and/or process errors are often the only sources of uncertainty modeled when addressing complex ecological problems, yet analyses should also account for uncertainty in sampling design, in model specification, in parameters governing the specified model, and in initial and boundary conditions. Only then can we be confident in the scientific inferences and forecasts made from an analysis. Probability and statistics provide a framework that accounts for multiple sources of uncertainty. Given the complexities of ecological studies, the hierarchical statistical model is an invaluable tool. This approach is not new in ecology, and there are many examples (both Bayesian and non-Bayesian) in the literature illustrating the benefits of this approach. In this article, we provide a baseline for concepts, notation, and methods, from which discussion on hierarchical statistical modeling in ecology can proceed. We have also planted some seeds for discussion and tried to show where the practical difficulties lie. Our thesis is that hierarchical statistical modeling is a powerful way of approaching ecological analysis in the presence of inevitable but quantifiable uncertainties, even if practical issues sometimes require pragmatic compromises

Relating Spatial Patterns of Stream Metabolism to Distributions of Juveniles Salmonids at the River Network Scale

Understanding the factors that drive spatial patterns in stream ecosystem processes and the distribution of aquatic biota is important to effective management of these systems and the conservation of biota at the network scale. In this study, we conducted field surveys throughout an extensive river network in NE Oregon that supports diminishing populations of wild salmonids. We collected data on physical habitat, nutrient concentrations, biofilm standing stocks, stream metabolism (gross primary production [GPP] and ecosystem respiration [ER]), and ESA‐listed juvenile salmonid density from approximately 50 sites across two sub‐basins. Our goals were to (1) to evaluate network patterns in these metrics, and (2) determine network‐scale linkages among these metrics, thus providing inference of processes driving observed patterns. Ambient nitrate‐N and phosphate‐P concentrations were low across both sub‐basins (\u3c40 μg/L). Nitrate‐N decreased with watershed area in both sub‐basins, but phosphate‐P only decreased in one sub‐basin. These spatial patterns suggest co‐limitation in one sub‐basin but N limitation in the other; experimental results using nutrient diffusing substrates across both sub‐basins supported these predictions. Solar exposure, temperature, GPP, ER, and GPP:ER increased with watershed area, but biofilm Chl a and ash‐free dry mass (AFDM) did not. Spatial statistical network (SSN) models explained between 70% and 75% of the total variation in biofilm Chl a, AFDM, and GPP, but only 21% of the variation in ER. Temperature and nutrient concentrations were the most supported predictors of Chl aand AFDM standing stocks, but these variables explained little of the total variation compared to spatial autocorrelation. In contrast, solar exposure and temperature were the most supported variables explaining GPP, and these variables explained far more variation than autocorrelation. Solar exposure, temperature, and nutrient concentrations explained almost none of the variation in ER. Juvenile salmonids—a key management focus in these sub‐basins—were most abundant in cool stream sections where rates of GPP were low, suggesting temperature constraints on these species restrict their distribution to oligotrophic areas where energy production at the base of the food web may be limited

Haul-Out Behavior of Harbor Seals (Phoca vitulina) in Hood Canal, Washington

The goal of this study was to model haul-out behavior of harbor seals (Phoca vitulina) in the Hood Canal region of Washington State with respect to changes in physiological, environmental, and temporal covariates. Previous research has provided a solid understanding of seal haul-out behavior. Here, we expand on that work using a generalized linear mixed model (GLMM) with temporal autocorrelation and a large dataset. Our dataset included behavioral haul-out records from archival and VHF radio tag deployments on 25 individual seals representing 61,430 seal hours. A novel application for increased computational efficiency allowed us to examine this large dataset with a GLMM that appropriately accounts for temporal autocorellation. We found significant relationships with the covariates hour of day, day of year, minutes from high tide and year. Additionally, there was a significant effect of the interaction term hour of day : day of year. This interaction term demonstrated that seals are more likely to haul out during nighttime hours in August and September, but then switch to predominantly daylight haul-out patterns in October and November. We attribute this change in behavior to an effect of human disturbance levels. This study also examined a unique ecological event to determine the role of increased killer whale (Orcinus orca) predation on haul-out behavior. In 2003 and 2005 these harbor seals were exposed to unprecedented levels of killer whale predation and results show an overall increase in haul-out probability after exposure to killer whales. The outcome of this study will be integral to understanding any changes in population abundance as a result of increased killer whale predation

Combining geostatistical and biotic interaction modelling to predict amphibian refuges under crayfish invasion across dendritic stream networks

Biodiversity ResearchAim: Biological invasions are pervasive in freshwater ecosystems, often causing native species to contract into areas that remain largely free from invasive species impacts. Predicting the location of such ecological refuges is challenging, because they are shaped by the habitat requirements of native and invasive species, their biotic interactions, and the spatial and temporal invasion patterns. Here, we investigated the spatial distribution and environmental drivers of refuges from invasion in river systems, by considering biotic interactions in geostatistical models accounting for stream network topology. We focused on Mediterranean amphibians negatively impacted by the invasive crayfishes Procambarus clarkii and Pacifastacus leniusculus. Location: River Sabor, NE Portugal. Methods: We surveyed amphibians at 168 200-m stream stretches in 2015. Geostatistical models were used to relate the probabilities of occurrence of each species to environmental and biotic variables, while controlling for linear (Euclidean) and hydrologic spatial dependencies. Biotic interactions were specified using crayfish probabilities of occurrence extracted from previously developed geostatistical models. Models were used to map the distribution of potential refuges for the most common amphibian species, under current conditions and future scenarios of crayfish expansion. Results: Geostatistical models were produced for eight out of 10 species detected, of which five species were associated with lower stream orders and only one species with higher stream orders. Six species showed negative responses to one or both crayfish species, even after accounting for environmental effects and spatial dependencies. Most amphibian species were found to retain large expanses of potential habitat in stream headwaters, but current refuges will likely contract under plausible scenarios of crayfish expansion. Main conclusions: Incorporating biotic interactions in geostatistical modelling provides a practical and relatively simple approach to predict present and future distributions of refuges from biological invasion in stream networks. Using this approach, our study shows that stream headwaters are key amphibian refuges under invasion by alien crayfishinfo:eu-repo/semantics/publishedVersio

Bayesian Multimodel Inference for Geostatistical Regression Models

The problem of simultaneous covariate selection and parameter inference for spatial regression models is considered. Previous research has shown that failure to take spatial correlation into account can influence the outcome of standard model selection methods. A Markov chain Monte Carlo (MCMC) method is investigated for the calculation of parameter estimates and posterior model probabilities for spatial regression models. The method can accommodate normal and non-normal response data and a large number of covariates. Thus the method is very flexible and can be used to fit spatial linear models, spatial linear mixed models, and spatial generalized linear mixed models (GLMMs). The Bayesian MCMC method also allows a priori unequal weighting of covariates, which is not possible with many model selection methods such as Akaike's information criterion (AIC). The proposed method is demonstrated on two data sets. The first is the whiptail lizard data set which has been previously analyzed by other researchers investigating model selection methods. Our results confirmed the previous analysis suggesting that sandy soil and ant abundance were strongly associated with lizard abundance. The second data set concerned pollution tolerant fish abundance in relation to several environmental factors. Results indicate that abundance is positively related to Strahler stream order and a habitat quality index. Abundance is negatively related to percent watershed disturbance

Fishery Discards: Factors Affecting Their Variability within a Demersal Trawl Fishery

Discards represent one of the most important issues within current commercial fishing. It occurs for a range of reasons and is influenced by an even more complex array of factors. We address this issue by examining the data collected within the Danish discard observer program and describe the factors that influence discarding within the Danish Kattegat demersal fleet over the period 1997 to 2008. Generalised additive models were used to assess how discards of the 3 main target species, Norway lobster, cod and plaice, and their subcomponents (under and over minimum landings size) are influenced by important factors and their potential relevance to management. Our results show that discards are influenced by a range of different factors that are different for each species and portion of discards. We argue that knowledge about the factors influential to discarding and their use in relation to potential mitigation measures are essential for future fisheries management strategies