skip to main content


Title: Semiparametric model selection for identification of environmental covariates related to adult groundfish catches and weights
Abstract

Ecologists and fisheries managers are interested in monitoring economically important marine fish species and using this data to inform management strategies. Determining environmental factors that best predict changes in these populations, particularly under rapid climate change, are a priority. I illustrate the application of the least squares-based spline estimation and group LASSO (LSSGLASSO) procedure for selection of coefficient functions in single index varying coefficient models (SIVCMs) on an ecological data set that includes spatiotemporal environmental covariates suspected to play a role in the catches and weights of six groundfish species. Temporal trends in variable selection were apparent, though the selection of variables was largely unrelated to common North Pacific climate indices. These results indicate that the strength of an environmental variable’s effect on a groundfish population may change over time, and not necessarily in-step with known low-frequency patterns of ocean-climate variability commonly attributable to large-scale regime shifts in the North Pacific. My application of the LSSGLASSO procedure for SIVCMs to deep water species using environmental data from various sources illustrates how variable selection with a flexible model structure can produce informative inference for remote and hard-to-reach animal populations.

 
more » « less
NSF-PAR ID:
10227418
Author(s) / Creator(s):
Publisher / Repository:
Nature Publishing Group
Date Published:
Journal Name:
Scientific Reports
Volume:
11
Issue:
1
ISSN:
2045-2322
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Quaternary climate fluctuations restructured biodiversity across North American high latitudes through repeated episodes of range contraction, population isolation and divergence, and subsequent expansion. Identifying how species responded to changing environmental conditions not only allows us to explore the mode and tempo of evolution in northern taxa, but also provides a basis for forecasting future biotic response across the highly variable topography of western North America. Using a multilocus approach under a Bayesian coalescent framework, we investigated the phylogeography of a wide‐ranging mammal, the long‐tailed vole,Microtus longicaudus. We focused on populations along the North Pacific Coast to refine our understanding of diversification by exploring the potentially compounding roles of multiple glacial refugia and more recent fragmentation of an extensive coastal archipelago. Through a combination of genetic data and species distribution models (SDMs), we found that historical climate variability influenced contemporary genetic structure, with multiple isolated locations of persistence (refugia) producing multiple divergent lineages (Beringian or northern, southeast Alaska or coastal, and southern or continental) during glacial advances. These vole lineages all occur along the North Pacific Coast where the confluence of numerous independent lineages in other species has produced overlapping zones of secondary contact, collectively a suture zone. Finally, we detected high levels of neoendemism due to complex island geography that developed in the last 10,000 years with the rising sea levels of the Holocene.

     
    more » « less
  2. Abstract

    Pacific salmon (Oncorhynchusspp.) are exposed to increased environmental change and multiple human stressors. To anticipate future impacts of global change and to improve sustainable resource management, it is critical to understand how wild salmon populations respond to stressors associated with human‐caused changes such as climate warming and ocean acidification, as well as competition in the ocean, which is intensified by the large‐scale production and release of hatchery reared salmon. Pink salmon (O.gorbuscha) are a keystone species in the North Pacific Ocean and support highly valuable commercial fisheries. We investigated the joint effects of changes in ocean conditions and salmon abundances on the productivity of wild pink salmon. Our analysis focused on Prince William Sound in Alaska, because the region accounts for ~50% of the global production of hatchery pink salmon with local hatcheries releasing 600–700 million pink salmon fry annually. Using 60 years of data on wild pink salmon abundances, hatchery releases, and ecological conditions in the ocean, we find evidence that hatchery pink salmon releases negatively affect wild pink salmon productivity, likely through competition between wild and hatchery juveniles in nearshore marine habitats. We find no evidence for effects of ocean acidification on pink salmon productivity. However, a change in the leading mode of North Pacific climate in 1988–1989 weakened the temperature–productivity relationship and altered the strength of intraspecific density dependence. Therefore, our results suggest non‐stationary (i.e., time varying) and interactive effects of ocean climate and competition on pink salmon productivity. Our findings further highlight the need for salmon management to consider potential adverse effects of large‐scale hatchery production within the context of ocean change.

     
    more » « less
  3. Abstract

    Fisheries management is dominated by the need to forecast catch and abundance of commercially and ecologically important species. The influence of spatial information and environmental factors on forecasting error is not often considered. I propose a forecasting method called spatiotemporally explicit model averaging (STEMA) to combine spatial and temporal information through model averaging. I examine the performance of STEMA against two popular forecasting models and a modern spatial prediction model: the autoregressive integrated moving averages with explanatory variables (ARIMAX) model, the Bayesian hierarchical model, and the varying coefficient model. I focus on applying the methods to four species of Alaskan groundfish for which catch data are available. My method reduces forecasting errors significantly for most of the tested models when compared to ARIMAX, Bayesian, and varying coefficient methods. I also consider the effect of sea surface temperature (SST) on the forecasting of catch, as multiple studies reveal a potential influence of water temperature on the survival and growth of juvenile groundfish. For most of the preferred models, inclusion of SST in the model improved forecasting of catch. It is advisable to consider both spatial information and relevant environmental factors in forecasting models to obtain more accurate projections of population abundance. The STEMA method is capable of accounting for spatial information in forecasting and can be applied to various types of data because of its flexible varying coefficient model structure. It is therefore a suitable forecasting method for application to many fields including ecology, epidemiology, and climatology.

     
    more » « less
  4. Abstract

    Identifying areas of high evolutionary potential is a judicious strategy for developing conservation priorities in the face of environmental change. For wide‐ranging species occupying heterogeneous environments, the evolutionary forces that shape distinct populations can vary spatially. Here, we investigate patterns of genomic variation and genotype–environment associations in the hermit thrush (Catharus guttatus), a North American songbird, at broad (across the breeding range) and narrow spatial scales (at a hybrid zone). We begin by building a genoscape or map of genetic variation across the breeding range and find five distinct genetic clusters within the species, with the greatest variation occurring in the western portion of the range. Genotype–environment association analyses indicate higher allelic turnover in the west than in the east, with measures of temperature surfacing as key predictors of putative adaptive genomic variation rangewide. Since broad patterns detected across a species' range represent the aggregate of many locally adapted populations, we investigate whether our broadscale analysis is consistent with a finer scale analysis. We find that top rangewide temperature‐associated loci vary in their clinal patterns (e.g., steep clines vs. fixed allele frequencies) across a hybrid zone in British Columbia, suggesting that the environmental predictors and the associated candidate loci identified in the rangewide analysis are of variable importance in this particular region. However, two candidate loci exhibit strong concordance with the temperature gradient in British Columbia, suggesting a potential role for temperature‐related barriers to gene flow and/or temperature‐driven ecological selection in maintaining putative local adaptation. This study demonstrates how patterns identified at the broad (macrogeographic) scale can be validated by investigating genotype–environment correlations at the local (microgeographic) scale. Furthermore, our results highlight the importance of considering the spatial distribution of putative adaptive variation when assessing population‐level sensitivity to climate change and other stressors.

     
    more » « less
  5. Abstract

    Genomic data are increasingly being integrated into macroecological forecasting, offering an evolutionary perspective that has been largely missing from global change biogeography. Genomic offset, which quantifies the disruption of genotype–environment associations under environmental change, allows for the incorporation of intraspecific climate‐associated genomic differentiation into forecasts of habitat suitability. Gradient Forest (GF) is a commonly used approach to estimate genomic offset; however, major hurdles in the application of GF‐derived genomic offsets are (1) an inability to interpret their absolute magnitude in an ecologically meaningful way and (2) uncertainty in how their implications compare with those of species‐level approaches like Ecological Niche Models (ENMs). Here, we assess the climate change vulnerability of red spruce (Picea rubens), a cool‐temperate tree species endemic to eastern North America, using both ENMs and GF modeling of genomic variation along climatic gradients. To gain better insights into climate change risks, we derive and apply two new threshold‐based genomic offset metrics—Donor and Recipient Importance—that quantify the transferability of propagules between donor populations and recipient localities while minimizing disruption of genotype–environment associations. We also propose and test a method for scaling genomic offsets relative to contemporary genomic variation across the landscape. In three common gardens, we found a significant negative relationship between (scaled) genomic offsets and red spruce growth and higher explanatory power for scaled offsets than climate transfer distances. However, the garden results also revealed the potential effects of spatial extrapolation and neutral genomic differentiation that can compromise the degree to which genomic offsets represent maladaptation and highlight the necessity of using common garden data to evaluate offset‐based predictions. ENMs and our novel genomic offset metrics forecasted drastic northward range shifts in suitable habitats. Combining inferences from our offset‐based metrics, we show that a northward shift mainly will be required for populations in the central and northern parts of red spruce's current range, whereas southern populations might persist in situ due to climate‐associated variation with less offset under future climate. These new genomic offset metrics thus yield refined, region‐specific prognoses for local persistence and show how management could be improved by considering assisted migration.

     
    more » « less