Environmental factors have stronger effects than biotic processes in patterns of intertidal populations along the southeast coast of Brazil

Rocky shore communities are shaped by complex interactions among environmental drivers and a range of biological processes. Here, we investigated the importance of abiotic and biotic drivers on the population structure of key rocky intertidal species at 62 sites, spanning ∼50% of the Brazilian rocky shoreline (i.e., ∼500 km). Large-scale population patterns were generally explained by differences in ocean temperature and wave exposure. For the gastropod species Lottia subrugosa, differences at smaller scales (i.e., 0.1-1 km) were better explained by other abiotic influences such as freshwater discharge and substrate roughness. Based on the general population patterns of intertidal species identified, three main oceanographic groups were observed: a cold-oligotrophic grouping at northern sites (Lakes sub-region), a eutrophic group associated with large estuaries and urban zones (Santos and Guanabara bays); and a transitional warm-water group found between the two more productive areas. Larger individuals of Stramonita brasiliensis, L. subrugosa and Echinolittorina lineolata were generally found in the cold-oligotrophic system (i.e., upwelling region), while small suspension feeders dominate the warm-eutrophic systems. Evidence of bottom-up regulation was not observed, and top-down regulation effects were only observed between the whelk S. brasiliensis and its mussel prey Pernaperna. Environmental drivers as compared to biotic interactions, therefore, play a key role determining the population structure of multiple intertidal species, across a range of spatial scales along the SW Atlantic shores.

A 25-Year Study of the Population Dynamics of a Harvested Population of Sika Deer on Kyushu Island, Japan

Sika deer (Cervus nippon) populations have damaged habitats, agricultural crops, and commercial forests in many parts of the world, including Asia, Europe, northern America, and New Zealand. Population management of sika deer is an important task in those areas. To better understand large-scale management and improve management efficiency, the authors estimated spatio-temporal changes of density distribution and population dynamics of a managed population of sika deer on Kyushu Island (approximately 36,750 km2), Japan. The authors estimated these changes by using fecal pellet count surveys conducted from 1995 to 2019 and results from a vector autoregressive spatio-temporal model. No decreasing trend of populations were observed at the island and prefectural scales, even though the management goal has been to reduce the population by half, and harvesting on the island increased annually until it reached about 110,000 sika deer in 2014. A possible explanation for the stable population dynamics is that the population used to determine the harvest number under the prefectural management plan was originally underestimated. This study highlights not only the difficulties of wide-area management of sika deer but also three important factors for successful management: reducing the risk of management failure, using an adaptive management approach, and appropriate management scale.

Titmice are a better indicator of bird density in Northern European than in Western European forests

Population sizes of many birds are declining alarmingly and methods for estimating fluctuations in species' abundances at a large spatial scale are needed. The possibility to derive indicators from the tendency of specific species to co-occur with others has been overlooked. Here, we tested whether the abundance of resident titmice can act as a general ecological indicator of forest bird density in European forests. Titmice species are easily identifiable and have a wide distribution, which makes them potentially useful ecological indicators. Migratory birds often use information on the density of resident birds, such as titmice, as a cue for habitat selection. Thus, the density of residents may potentially affect community dynamics. We examined spatio-temporal variation in titmouse abundance and total bird abundance, each measured as biomass, by using long-term citizen science data on breeding forest birds in Finland and France. We analyzed the variation in observed forest bird density (excluding titmice) in relation to titmouse abundance. In Finland, forest bird density linearly increased with titmouse abundance. In France, forest bird density nonlinearly increased with titmouse abundance, the association weakening toward high titmouse abundance. We then analyzed whether the abundance (measured as biomass) of random species sets could predict forest bird density better than titmouse abundance. Random species sets outperformed titmice as an indicator of forest bird density only in 4.4% and 24.2% of the random draws, in Finland and France, respectively. Overall, the results suggest that titmice could act as an indicator of bird density in Northern European forest bird communities, encouraging the use of titmice observations by even less-experienced observers in citizen science monitoring of general forest bird density.

Finite-sample properties of estimators for first and second order autoregressive processes

The class of autoregressive (AR) processes is extensively used to model temporal dependence in observed time series. Such models are easily available and routinely fitted using freely available statistical software like . A potential problem is that commonly applied estimators for the coefficients of AR processes are severely biased when the time series are short. This paper studies the finite-sample properties of well-known estimators for the coefficients of stationary AR(1) and AR(2) processes and provides bias-corrected versions of these estimators which are quick and easy to apply. The new estimators are constructed by modeling the relationship between the true and originally estimated AR coefficients using weighted orthogonal polynomial regression, taking the sampling distribution of the original estimators into account. The finite-sample distributions of the new bias-corrected estimators are approximated using transformations of skew-normal densities, combined with a Gaussian copula approximation in the AR(2) case. The properties of the new estimators are demonstrated by simulations and in the analysis of a real ecological data set. The estimators are easily available in our accompanying -package for AR(1) and AR(2) processes of length 10–50, both giving bias-corrected coefficient estimates and corresponding confidence intervals.

Age-dependent patterns of spatial autocorrelation in fish populations

The degree of spatial autocorrelation in population fluctuations increases with dispersal and geographical covariation in the environment, and decreases with strength of density dependence. Because the effects of these processes can vary throughout an individual's lifespan, we studied how spatial autocorrelation in abundance changed with age in three marine fish species in the Barents Sea. We found large interspecific differences in age-dependent patterns of spatial autocorrelation in density. Spatial autocorrelation increased with age in cod, the reverse trend was found in beaked redfish, while it remained constant among age classes in haddock. We also accounted for the average effect of local cohort dynamics, i.e. the expected local density of an age class given last year's local density of the cohort, with the goal of disentangling spatial autocorrelation patterns acting on an age class from those formed during younger age classes and being carried over. We found that the spatial autocorrelation pattern of older age classes became increasingly determined by the distribution of the cohort during the previous year. Lastly, we found high degrees of autocorrelation over long distances for the three species, suggesting the presence of far-reaching autocorrelating processes on these populations. We discuss how differences in the species' life history strategies could cause the observed differences in age-specific variation in spatial autocorrelation. As spatial autocorrelation can differ among age classes, our study indicates that fluctuations in age structure can influence the spatio-temporal variation in abundance of marine fish populations.

Harmonizing child mortality data at disparate geographic levels

There is an increasing focus on reducing inequalities in health outcomes in developing countries. Subnational variation is of particular interest, with geographically-indexed data being used to understand the spatial risk of detrimental outcomes and to identify who is at greatest risk. While some health surveys provide observations with associated geographic coordinates (point data), many others provide data that have their locations masked and instead only report the strata (polygon information) within which the data resides (masked data). How to harmonize these data sources for spatial analysis has been previously considered although only ad hoc methods and comparison of methods is lacking. In this paper, we present a new method for analyzing masked survey data, using a method that is consistent with the data-generating process. In addition, we critique two previously proposed approaches to analyzing masked data and illustrate that they are fundamentally flawed methodologically. To validate our method, we compare our approach with previously formulated solutions in several realistic simulation environments in which the underlying structure of the risk field is known. We simulate samples from spatiotemporal fields in a way that mimics the sampling frame implemented in the most common health surveys in low- and middle-income countries, the Demographic and Health Surveys and Multiple Indicator Cluster Surveys. In simulations, the newly proposed approach outperforms previously proposed approaches in terms of minimizing error while increasing the precision of estimates. The approaches are subsequently compared using child mortality data from the Dominican Republic where our findings are reinforced. The ability to accurately increase precision of child mortality estimates, and health outcomes in general, by leveraging various types of data, improves our ability to implement precision public health initiatives and better understand the landscape of geographic health inequalities.

Environmental conditions modulate compensatory effects of site dependence in a food-caching passerine

Although density regulates the abundance of most wild animal populations by influencing vital rates, such as fecundity and survival, the mechanisms responsible for generating negative density dependence are unclear for many species. Site dependence occurs when there is preferential filling of high-quality territories, which results in higher per capita vital rates at low densities because a larger proportion of occupied territories are of high quality. Using 41 yr of territory occupancy and demographic data, we investigated whether site dependence was a mechanism acting to influence fecundity and, by extension, regulate a population of Canada Jays in Algonquin Provincial Park, Ontario, Canada. As predicted by site dependence, the proportion of occupied territories that were of high quality was negatively correlated with population density and periods of vacancy were shorter for high-quality territories than for low-quality territories. We also found evidence that per capita fecundity was positively related to the proportion of occupied territories that were of high quality, but only when environmental conditions, which influence the entire population, were otherwise poor for breeding. Our results suggest that site dependence likely plays a role in regulating this population but that environmental conditions can modulate the strength of density dependence.

The impact of human land use and landscape productivity on population dynamics of red fox in southeastern Norway

In the boreal forest, the red fox (Vulpes vulpes) is a key species due to its many strong food web linkages and its exploitation of niches that form in the wake of human activities. Recent altitudinal range expansion and a perceived population increase have become topics of concern in Scandinavia, primarily due to the potential impacts of red foxes on both prey and competitor species. However, despite it being a common species, there is still surprisingly little knowledge about the temporal and spatial characteristics of its population dynamics. In this study, we synthesized 12 years of snow-track transect data covering a 27,000-km2 study area to identify factors associated with red fox distribution and population dynamics. Using Bayesian hierarchical regression models, we evaluated the relationships of landscape productivity and climate gradients as well as anthropogenic subsidies with an index of red fox population size and growth rates. We found that landscapes with high human settlement density and large amounts of gut piles from moose (Alces alces) hunting were associated with higher red fox abundances. Population dynamics were characterized by direct density-dependent growth, and the structure of density dependence was best explained by the amount of agricultural land in the landscape. Population equilibrium levels increased, and populations were more stable, in areas with a higher presence of agricultural lands, whereas density-dependent population growth was more prominent in areas of low agricultural presence. We conclude that human land use is a dominant driver of red fox population dynamics in the boreal forest. We encourage further research focusing on contrasting effects of anthropogenic subsidization on predator population carrying capacities and temporal stability, and potential impacts on prey dynamics.

Hierarchical modeling strengthens evidence for density dependence in observational time series of population dynamics

The extent to which populations in nature are regulated by density-dependent processes is unresolved. While experiments increasingly find evidence of strong density dependence, unmanipulated population time series yield much more ambiguous evidence of regulation, especially when accounting for effects of observation error. Here, we reexamine the evidence for density dependence in time series of population sizes in nature, by conducting an aggregate analysis of the populations in the Global Population Dynamics Database (GPDD). First, following the conventional approach, we fit a density-dependent and a density-independent variant of the Gompertz state-space model to each time series. Then, we conduct an aggregate analysis of the entire database by considering two random-effects density-dependent models that leverage information across data sets. When individual time series are tested independently, we find very little evidence for density dependence. However, in the aggregate, we find very strong evidence for density dependence, even though, paradoxically, estimated strengths of density dependence for individual time series tend to be weaker than when each individual time series is analyzed independently. Furthermore, a hierarchical model that accounts for taxonomic variation in the strength of density dependence reveals that density dependence is consistently stronger in insects and fish than in birds and mammals. Our findings resolve apparent inconsistencies between observational and experimental studies of density dependence by revealing that the observational record does indeed contain strong support for the hypothesis that density dependence is widespread in nature.

Black swans in space: modeling spatiotemporal processes with extremes

In ecological systems, extremes can happen in time, such as population crashes, or in space, such as rapid range contractions. However, current methods for joint inference about temporal and spatial dynamics (e.g., spatiotemporal modeling with Gaussian random fields) may perform poorly when underlying processes include extreme events. Here we introduce a model that allows for extremes to occur simultaneously in time and space. Our model is a Bayesian predictive-process GLMM (generalized linear mixed-effects model) that uses a multivariate-t distribution to describe spatial random effects. The approach is easily implemented with our flexible R package glmmfields. First, using simulated data, we demonstrate the ability to recapture spatiotemporal extremes, and explore the consequences of fitting models that ignore such extremes. Second, we predict tree mortality from mountain pine beetle (Dendroctonus ponderosae) outbreaks in the U.S. Pacific Northwest over the last 16 yr. We show that our approach provides more accurate and precise predictions compared to traditional spatiotemporal models when extremes are present. Our R package makes these models accessible to a wide range of ecologists and scientists in other disciplines interested in fitting spatiotemporal GLMMs, with and without extremes.

Integrating population dynamics models and distance sampling data: a spatial hierarchical state-space approach

Stochastic versions of Gompertz, Ricker, and various other dynamics models play a fundamental role in quantifying strength of density dependence and studying long-term dynamics of wildlife populations. These models are frequently estimated using time series of abundance estimates that are inevitably subject to observation error and missing data. This issue can be addressed with a state-space modeling framework that jointly estimates the observed data model and the underlying stochastic population dynamics (SPD) model. In cases where abundance data are from multiple locations with a smaller spatial resolution (e.g., from mark-recapture and distance sampling studies), models are conventionally fitted to spatially pooled estimates of yearly abundances. Here, we demonstrate that a spatial version of SPD models can be directly estimated from short time series of spatially referenced distance sampling data in a unified hierarchical state-space modeling framework that also allows for spatial variance (covariance) in population growth. We also show that a full range of likelihood based inference, including estimability diagnostics and model selection, is feasible in this class of models using a data cloning algorithm. We further show through simulation experiments that the hierarchical state-space framework introduced herein efficiently captures the underlying dynamical parameters and spatial abundance distribution. We apply our methodology by analyzing a time series of line-transect distance sampling data for fin whales (Balaenoptera physalus) off the U.S. west coast. Although there were only seven surveys conducted during the study time frame, 1991-2014, our analysis detected presence of strong density regulation and provided reliable estimates of fin whale densities. In summary, we show that the integrative framework developed herein allows ecologists to better infer key population characteristics such as presence of density regulation and spatial variability in a population's intrinsic growth potential.

A geostatistical state-space model of animal densities for stream networks

Population dynamics are often correlated in space and time due to correlations in environmental drivers as well as synchrony induced by individual dispersal. Many statistical analyses of populations ignore potential autocorrelations and assume that survey methods (distance and time between samples) eliminate these correlations, allowing samples to be treated independently. If these assumptions are incorrect, results and therefore inference may be biased and uncertainty underestimated. We developed a novel statistical method to account for spatiotemporal correlations within dendritic stream networks, while accounting for imperfect detection in the surveys. Through simulations, we found this model decreased predictive error relative to standard statistical methods when data were spatially correlated based on stream distance and performed similarly when data were not correlated. We found that increasing the number of years surveyed substantially improved the model accuracy when estimating spatial and temporal correlation coefficients, especially from 10 to 15 yr. Increasing the number of survey sites within the network improved the performance of the nonspatial model but only marginally improved the density estimates in the spatiotemporal model. We applied this model to brook trout data from the West Susquehanna Watershed in Pennsylvania collected over 34 yr from 1981 to 2014. We found the model including temporal and spatiotemporal autocorrelation best described young of the year (YOY) and adult density patterns. YOY densities were positively related to forest cover and negatively related to spring temperatures with low temporal autocorrelation and moderately high spatiotemporal correlation. Adult densities were less strongly affected by climatic conditions and less temporally variable than YOY but with similar spatiotemporal correlation and higher temporal autocorrelation.

Spatial separation of catches in highly mixed fisheries

Mixed fisheries are the dominant type of fishery worldwide. Overexploitation in mixed fisheries occurs when catches continue for available quota species while low quota species are discarded. As EU fisheries management moves to count all fish caught against quota (the “landing obligation”), the challenge is to catch available quota within new constraints, else lose productivity. A mechanism for decoupling exploitation of species caught together is spatial targeting, which remains challenging due to complex fishery and population dynamics. How far spatial targeting can go to practically separate species is often unknown and anecdotal. We develop a dimension-reduction framework based on joint dynamic species distribution modelling to understand how spatial community and fishery dynamics interact to determine species and size composition. In application to the highly mixed fisheries of the Celtic Sea, clear common spatial patterns emerge for three distinct assemblages. While distribution varies interannually, the same species are consistently found in higher densities together, with more subtle differences within assemblages, where spatial separation may not be practically possible. We highlight the importance of dimension reduction techniques to focus management discussion on axes of maximal separation and identify spatiotemporal modelling as a scientific necessity to address the challenges of managing mixed fisheries.

Fine-scale population dynamics in a marine fish species inferred from dynamic state-space models

Identifying the spatial scale of population structuring is critical for the conservation of natural populations and for drawing accurate ecological inferences. However, population studies often use spatially aggregated data to draw inferences about population trends and drivers, potentially masking ecologically relevant population sub-structure and dynamics. The goals of this study were to investigate how population dynamics models with and without spatial structure affect inferences on population trends and the identification of intrinsic drivers of population dynamics (e.g. density dependence). Specifically, we developed dynamic, age-structured, state-space models to test different hypotheses regarding the spatial structure of a population complex of coastal Atlantic cod (Gadus morhua). Data were from a 93-year survey of juvenile (age 0 and 1) cod sampled along >200 km of the Norwegian Skagerrak coast. We compared two models: one which assumes all sampled cod belong to one larger population, and a second which assumes that each fjord contains a unique population with locally determined dynamics. Using the best supported model, we then reconstructed the historical spatial and temporal dynamics of Skagerrak coastal cod. Cross-validation showed that the spatially structured model with local dynamics had better predictive ability. Furthermore, posterior predictive checks showed that a model which assumes one homogeneous population failed to capture the spatial correlation pattern present in the survey data. The spatially structured model indicated that population trends differed markedly among fjords, as did estimates of population parameters including density-dependent survival. Recent biomass was estimated to be at a near-record low all along the coast, but the finer scale model indicated that the decline occurred at different times in different regions. Warm temperatures were associated with poor recruitment, but local changes in habitat and fishing pressure may have played a role in driving local dynamics. More generally, we demonstrated how state-space models can be used to test evidence for population spatial structure based on survey time-series data. Our study shows the importance of considering spatially structured dynamics, as the inferences from such an approach can lead to a different ecological understanding of the drivers of population declines, and fundamentally different management actions to restore populations.

Estimating partial regulation in spatiotemporal models of community dynamics

Niche-based approaches to community analysis often involve estimating a matrix of pairwise interactions among species (the "community matrix"), but this task becomes infeasible using observational data as the number of modeled species increases. As an alternative, neutral theories achieve parsimony by assuming that species within a trophic level are exchangeable, but generally cannot incorporate stabilizing interactions even when they are evident in field data. Finally, both regulated (niche) and unregulated (neutral) approaches have rarely been fitted directly to survey data using spatiotemporal statistical methods. We therefore propose a spatiotemporal and model-based approach to estimate community dynamics that are partially regulated. Specifically, we start with a neutral spatiotemporal model where all species follow ecological drift, which precludes estimating pairwise interactions. We then add regulatory relations until model selection favors stopping, where the "rank" of the interaction matrix may range from zero to the number of species. A simulation experiment shows that model selection can accurately identify the rank of the interaction matrix, and that the identified spatiotemporal model can estimate the magnitude of species interactions. A 40-yr case study for the Gulf of St. Lawrence marine community shows that recovering grey seals have an unregulated and negative relationship with demersal fishes. We therefore conclude that partial regulation is a plausible approximation to community dynamics using field data and hypothesize that estimating partial regulation will be expedient in future analyses of spatiotemporal community dynamics given limited field data. We conclude by recommending ongoing research to add explicit models for movement, so that meta-community theory can be confronted with data in a spatiotemporal statistical framework.

Interactive effects of shelter and conspecific density shape mortality, growth, and condition in juvenile reef fish

How landscape context influences density-dependent processes is important, as environmental heterogeneity can confound estimates of density dependence in demographic parameters. Here we evaluate 19 populations in a shoaling temperate reef fish (Trachinops caudimaculatus) metapopulation within a heterogeneous seascape (Port Phillip Bay, Australia) to show empirically that shelter availability and population density interact to influence juvenile mortality, growth and condition. Although heterogeneity in shelter availability obscured the underlying patterns of density dependence in different ways, the combination of habitat and its interaction with density were two to six times more important than density alone in explaining variation in demographic parameters for juveniles. These findings contradict many small-scale studies and highlight the need for landscape-scale observations of how density dependence interacts with resource availability and competition to better understand how demographic parameters influence the dynamics of metapopulations in heterogeneous environments.

Precision of systematic and random sampling in clustered populations: habitat patches and aggregating organisms

Natural populations of plants and animals spatially cluster because (1) suitable habitat is patchy, and (2) within suitable habitat, individuals aggregate further into clusters of higher density. We compare the precision of random and systematic field sampling survey designs under these two processes of species clustering. Second, we evaluate the performance of 13 estimators for the variance of the sample mean from a systematic survey. Replicated simulated surveys, as counts from 100 transects, allocated either randomly or systematically within the study region, were used to estimate population density in six spatial point populations including habitat patches and Matérn circular clustered aggregations of organisms, together and in combination. The standard one-start aligned systematic survey design, a uniform 10 x 10 grid of transects, was much more precise. Variances of the 10 000 replicated systematic survey mean densities were one-third to one-fifth of those from randomly allocated transects, implying transect sample sizes giving equivalent precision by random survey would need to be three to five times larger. Organisms being restricted to patches of habitat was alone sufficient to yield this precision advantage for the systematic design. But this improved precision for systematic sampling in clustered populations is underestimated by standard variance estimators used to compute confidence intervals. True variance for the survey sample mean was computed from the variance of 10 000 simulated survey mean estimates. Testing 10 published and three newly proposed variance estimators, the two variance estimators (v) that corrected for inter-transect correlation (ν₈ and ν(W)) were the most accurate and also the most precise in clustered populations. These greatly outperformed the two "post-stratification" variance estimators (ν₂ and ν₃) that are now more commonly applied in systematic surveys. Similar variance estimator performance rankings were found with a second differently generated set of spatial point populations, ν₈ and ν(W) again being the best performers in the longer-range autocorrelated populations. However, no systematic variance estimators tested were free from bias. On balance, systematic designs bring more narrow confidence intervals in clustered populations, while random designs permit unbiased estimates of (often wider) confidence interval. The search continues for better estimators of sampling variance for the systematic survey mean.

Using spatiotemporal species distribution models to identify temporally evolving hotspots of species co-occurrence

Identifying spatiotemporal hotspots is important for understanding basic ecological processes, but is particularly important for species at risk. A number of terrestrial and aquatic species are indirectly affected by anthropogenic impacts, simply because they tend to be associated with species that are targeted for removals. Using newly developed statistical models that allow for the inclusion of time-varying spatial effects, we examine how the co-occurrence of a targeted and nontargeted species can be modeled as a function of environmental covariates (temperature, depth) and interannual variability. The nontarget species in our case study (eulachon) is listed under the U.S. Endangered Species Act, and is encountered by fisheries off the U.S. West Coast that target pink shrimp. Results from our spatiotemporal model indicated that eulachon bycatch risk decreases with depth and has a convex relationship with sea surface temperature. Additionally, we found that over the 2007-2012 period, there was support for an increase in eulachon density from both a fishery data set (+40%) and a fishery-independent data set (+55%). Eulachon bycatch has increased in recent years, but the agreement between these two data sets implies that increases in bycatch are not due to an increase in incidental targeting of eulachon by fishing vessels, but because of an increasing population size of eulachon. Based on our results, the application of spatiotemporal models to species that are of conservation concern appears promising in identifying the spatial distribution of environmental and anthropogenic risks to the population.

On Extrapolating Past the Range of Observed Data When Making Statistical Predictions in Ecology

Ecologists are increasingly using statistical models to predict animal abundance and occurrence in unsampled locations. The reliability of such predictions depends on a number of factors, including sample size, how far prediction locations are from the observed data, and similarity of predictive covariates in locations where data are gathered to locations where predictions are desired. In this paper, we propose extending Cook's notion of an independent variable hull (IVH), developed originally for application with linear regression models, to generalized regression models as a way to help assess the potential reliability of predictions in unsampled areas. Predictions occurring inside the generalized independent variable hull (gIVH) can be regarded as interpolations, while predictions occurring outside the gIVH can be regarded as extrapolations worthy of additional investigation or skepticism. We conduct a simulation study to demonstrate the usefulness of this metric for limiting the scope of spatial inference when conducting model-based abundance estimation from survey counts. In this case, limiting inference to the gIVH substantially reduces bias, especially when survey designs are spatially imbalanced. We also demonstrate the utility of the gIVH in diagnosing problematic extrapolations when estimating the relative abundance of ribbon seals in the Bering Sea as a function of predictive covariates. We suggest that ecologists routinely use diagnostics such as the gIVH to help gauge the reliability of predictions from statistical models (such as generalized linear, generalized additive, and spatio-temporal regression models).

Using citizen-science data to identify local hotspots of seabird occurrence

Seabirds have been identified and used as indicators of ecosystem processes such as climate change and human activity in nearshore ecosystems around the globe. Temporal and spatial trends have been documented at large spatial scales, but few studies have examined more localized patterns of spatiotemporal variation, by species or functional group. In this paper, we apply spatial occupancy models to assess the spatial patchiness and interannual trends of 18 seabird species in the Puget Sound region (Washington State, USA). Our dataset, the Puget Sound Seabird Survey of the Seattle Audubon Society, is unique in that it represents a seven-year study, collected with a focus on winter months (October-April). Despite historic declines of seabirds in the region over the last 50 years, results from our study are optimistic, suggesting increases in probabilities of occurrence for 14 of the 18 species included. We found support for declines in occurrence for white-winged scoters, brants, and 2 species of grebes. The decline of Western grebes in particular is troubling, but in agreement with other recent studies that have shown support for a range shift south in recent years, to the southern end of California Current.