Metapopulation viability of an endangered shorebird depends on dispersal and human-created habitats: piping plovers (Charadrius melodus) and prairie rivers
© Catlin et al. 2016
Received: 14 October 2015
Accepted: 15 February 2016
Published: 15 March 2016
Many species are distributed as metapopulations in dynamic landscapes, where habitats change through space and time. Individuals locate habitat through dispersal, and the relationship between a species and landscape characteristics can have profound effects on population persistence. Despite the importance of connectivity in dynamic environments, few empirical studies have examined temporal variability in dispersal or its effect on metapopulation dynamics. In response to this knowledge gap, we studied the dispersal, demography, and viability of a metapopulation of an endangered, disturbance-dependent shorebird. We examined three subpopulations of piping plovers (Charadrius melodus) on the lower Platte and Missouri rivers from 2008–2013. High flow events from an upstream dam on the Missouri River in 2010 and 2011 allowed us to assess the effect of total habitat loss and the subsequent creation of new habitat associated with a large disturbance at one ‘natural’ study location. The other two sites within the metapopulation, which were maintained by anthropogenic activities (e.g., mining, development, habitat restoration), were largely unaffected by this disturbance, resulting in a controlled natural experiment.
High flow events were associated with increased emigration, decreased immigration, and decreased survival in the subpopulation that experienced high flows. Following the high flow event, immigration into that subpopulation increased. Dispersal rates among subpopulations were negatively correlated with distance. The metapopulation had a low probability of extinction over 100 years (0 %) under the current disturbance interval and associated dispersal and survival rates. However, persistence depended on relatively stable, human-created habitats, not the dynamic, natural habitat (47.7 % extinction probability for this subpopulation).
We found that functional connectivity, as measured by the rate of dispersal among subpopulations, increased as a result of the high flow event in our study metapopulation. Plovers also increased reproductive output following this event. Although the study metapopulation had a low overall probability of extinction, metapopulation persistence depended on anthropogenically created habitats that provided a small but stable source of nesting habitat and dispersers through time. However, all subpopulations remained small, even if persistent, making them individually vulnerable to extinction through stochastic events. Given the highly dynamic nature of habitat availability in this system, maintaining several subpopulations within the metapopulation and stable sources of habitat will be critical, and this species will likely remain conservation-reliant.
Habitat patchiness, spatial subdivision, and local extinction/recolonization dynamics are common in ecological systems, and many species are organized into metapopulations (i.e., a group of local subpopulations that inhabit discrete habitat patches but interact through dispersal) as a consequence . Metapopulation theory and empirical studies of species distributed as metapopulations are valuable tools that inform conservation strategies for imperiled species, with applications in reserve design, corridor creation, and extinction risk estimation [2, 3]. As such, there has been a proliferation of theoretical and empirical studies of metapopulations and their properties.
Although basic metapopulation theory deals primarily with static systems, there is growing interest in its application to dynamic environments that more closely mimic natural systems . Populations in dynamic systems experience local extinctions and colonization opportunities through both stochastic and deterministic (e.g., succession) processes [5, 6]. Metapopulation persistence in dynamic landscapes depends on the species’ demographic characteristics as well as underlying spatial (e.g., area, connectivity, patch size, patch quality) and temporal (e.g., turnover rate, extent, intensity) habitat properties [3, 7]. This interplay between species characteristics and the dynamic properties of a metapopulation and its habitat have important consequences for how disturbances and changing landscapes affect a species. Relatively few studies, however, have coupled empirical species data with dynamic landscape modelling to produce actionable management plans [3, 7].
The application of dynamic metapopulation models has clear connections with the management of imperiled species distributed as metapopulations naturally or due to anthropogenic fragmentation [8, 9]. As the study of metapopulations in dynamic landscapes matures, it is important to merge information on imperiled species ecology and landscape dynamism in empirical studies . Such studies should relate manageable features of dynamic systems to species persistence, thereby providing qualitative and quantitative recommendations for the management of dynamic systems . In addition, studies need to focus on less well-studied factors in dynamic systems (e.g., connectivity, variation in patch quality) , such that generalizations across variable landscapes can be made.
Structural and functional connectivity are fundamental characteristics affecting metapopulation persistence in dynamic landscapes because, as the spatial and temporal orientation of habitat changes, animals must locate habitat through dispersal [6, 10]. For instance, in dynamic, disturbed ecosystems, connectivity among subpopulations can allow (i) individuals to disperse away from areas experiencing a disturbance to refugia and (ii) surviving individuals to later recolonize recently disturbed or created habitat [6, 11, 12]. However, temporal changes in connectivity are rarely studied or modeled . To fully understand metapopulation persistence in dynamic, disturbance-dependent systems, considering changes in connectivity due to natural or anthropogenic forces and developing metapopulation models that account for both spatial and temporal changes in connectivity is necessary .
We studied a disturbance-dependent animal that displays high site fidelity despite a capacity for long-distance dispersal, the piping plover (Charadrius melodus; hereafter ‘plover’), within a metapopulation on the Missouri and Platte rivers. Plovers in this area nest on a variety of habitats with a range of turnover and disturbance rates. We conducted this study (i) to determine the effect of disturbance and habitat turnover on dispersal and functional connectivity in the metapopulation; (ii) to use empirical data to compare the effects of human alterations to habitat turnover (i.e., subpopulation-specific habitat management) among dynamic and static metapopulation patches; and (iii) to make suggestions about the management of this system relative to functional connectivity and turnover to promote persistence.
The piping plover is a migratory, precocial shorebird that nests on sparsely vegetated beaches on the Atlantic coast, Great Lakes, and Great Plains in the United States and Canada . Plover adult annual survival averages 76 %, resulting in an approximately 5 years average lifespan . Some birds breed in their first year post-hatch but do so approximately one month later than experienced adults . Across their range, disturbances (e.g., high-water events on rivers, reservoirs, and alkali lakes; coastal storms) are critical for maintaining early successional habitat. Plovers were listed under the U.S. Endangered Species Act in 1986, primarily as a result of habitat loss and low reproductive output .
Plover populations are distributed as metapopulations throughout their geographic range, both naturally and as a result of anthropogenic landcover change. In the southern Great Plains, plover subpopulations on the Missouri River nest on sandbars  and sand and gravel mines and lakeshore housing developments in the lower Platte River floodplain . Nesting areas are separated from other nesting areas by inhospitable expanses of channelized river, impoundments , development , and agricultural lands. In North Dakota, plovers disperse among subpopulations  that inhabit naturally isolated alkali lakes, riverine sandbars, and reservoirs , where flooding can have profound impacts on habitat and demography . A similar pattern has been observed in Saskatchewan, where short- and long-term flooding has affected plover distribution and dispersal . Outside of the Great Plains, plovers use ephemeral barrier island habitat created through storms and degraded through succession, erosion, and redevelopment [21, 22]. Thus, throughout their range, plovers are subject to local extinctions due to natural (e.g., flooding, succession) and anthropogenic (e.g., development) processes, relying on disturbance to refresh habitat and dispersal to recolonize newly available habitat.
Despite the metapopulation structure commonly observed for this species, relatively little is known about plover dispersal among subpopulations and across spatiotemporal scales or the implications of dispersal on metapopulation persistence (review in ), which could limit conservation options (e.g., ) for this Federally threatened species. One study related increased dispersal among nearby subpopulations to reproductive failure and flooding but did not evaluate the effect of disturbance and dispersal on persistence nor the effect of distance on connectivity . An understanding of dispersal rates and connectivity based on distance could allow for improved regional models or forecasts of metapopulation persistence. Plover management has generally been specific to a subpopulation or a managing entity, increasing the difficulty of obtaining the appropriate information needed to understand metapopulation dynamics and managing the species at a metapopulation scale [15, 18].
Observations also support a metapopulation structure for these subpopulations. Dispersal events are commonly observed among our study subpopulations (result herein). However, during our study, no birds banded in our subpopulations were found breeding in populations north of our study area, and only one individual banded at Lake Sakakawea, North Dakota was known to breed in our study area (>780 km; DHC and MBB, unpublished data). We have also received limited reports of banded birds from our subpopulations nesting elsewhere. Most of these birds (8 total) bred in areas immediately adjacent (<5 km) to Lewis and Clark Lake in relatively small blocks of ephemeral habitat (DHC, personal observation). Two birds banded on Lewis and Clark Lake and the Gavins Point Reach were found nesting within another small population with more stable habitat in the upper Niobrara River near Spencer, Nebraska (>50 km). There was also limited movement between our subpopulations and another subpopulation on the central Platte River, near Kearney, Nebraska. Three birds banded on the central Platte were found nesting on the Gavins Reach and Lewis and Clark Lake (>260 km), and five were found nesting in the lower Platte subpopulation (>150 km). In contrast, no birds banded on Lewis and Clark Lake or the Gavins Point Reach were found nesting on the central Platte, and only one bird from the Lower Platte subpopulation was found nesting there (DHC and MBB, unpublished data). Most of these sightings were at locations adjacent to our metapopulation, and movements appeared to increase following high flow events. Therefore, the metapopulation we considered here was open, but the interchange of individuals between our subpopulations and others appeared to follow the same patterns as that within our metapopulation and likely did not have a strong effect on the observed metapopulation dynamics.
Historically, plovers within our metapopulation nested on sandbars within the two rivers, where seasonal water level fluctuations maintained early successional habitat conditions [23, 24]. Before they were altered, the Missouri and Platte rivers experienced peak water levels in March and June due to prairie and mountain snowmelt and precipitation, coinciding with current peak nesting times for plovers . These peak flows submerged existing sandbars and redistributed sediments, creating unvegetated sandbars suitable for plover nesting as water levels receded either within that season or in following seasons [15, 24]. We define a bankfull flow as a high flow event where a river’s water level fills the channel to the top of the river bank and begins to overflow onto the floodplain, while water levels often extend farther into the floodplain during a flood event . Because both bankfull flows and floods are capable of moving sediment and forming or removing sandbars , we hereafter combine these terms and refer to them as “high flow events” or “high flows” rather than simply floods. Historically, these rivers experienced more extensive high flow events, approximately once every 4–6 years . Currently, high flow events are less regular and occur approximately every 20 years (based on recent discharge rates from Gavins Point Dam available at http://www.nwd-mr.usace.army.mil/rcc/projdata/gapt.pdf).
The natural hydrographs for these rivers have been altered due to channelization, bank reinforcement, water diversion, and dam construction. Upstream mainstem dam operation and water diversion in the Platte River drainage has significantly altered flow regimes and channel characteristics [26, 27]. Similarly, roughly 35 % of the Missouri River has been impounded in lake ecosystems, and an additional 32 % of the river has been channelized . Dams and reservoirs reduce high flow volumes needed to scour vegetation from existing sandbars and carry sediments downstream . The decline in sediment volume and decreased high flow frequency has decreased the amount of sandbar habitat on parts of the Missouri River by 96 % from 1892 to 2006 . This lack of habitat is now the primary threat to the persistence of plovers in the Great Plains .
Today, the majority of adults (77 %; MBB, LRD, JGJ, unpublished data) in the lower Platte River subpopulation (hereafter ‘HC,’ for human created) nest off-river on human-created sand and gravel mines and lakeshore housing developments adjacent to the lower Platte River and a small portion (60 km) of the Loup River . Eighty-one percent of the birds nesting on HC were clustered in the southeastern 1/3 of the study area (MBB, LRD, JGJ, unpublished data; Fig. 1).
The Gavins Point Reach subpopulation occurs in one of the last free-flowing portions of the Missouri River and extends 95 km downstream from the Gavins Point Dam (42° 51′N, 97°29′W). The Lewis and Clark Lake subpopulation is within a reservoir impounded by the dam that contains nesting sandbars at the upstream end (42° 51′N, 97°47′W). From 2005 to 2010, the U.S. Army Corps of Engineers (USACE) built sandbars at both locations to provide nesting habitat for plovers and interior least terns (Sternula antillarum athalassos; ). Sandbar nesting habitat used by birds on the Gavins Point Reach and Lewis and Clark Lake was generally composed of low, unvegetated mud and sandflats with higher elevation areas of either barren sand or vegetation dominated by cottonwood (Populus spp.) and willow (Salix spp.) saplings.
In June 2010, unusually large volumes of water were released from Gavins Point Dam, which flooded all active nests and hatched chicks on Gavins Point Reach (hereafter ‘M1F,’ Missouri River site with high flow), reducing reproductive output at this site to 0. In 2011, increased mountain snowpack and spring precipitation resulted in historically high water levels in the Missouri River, and water levels covered all nesting habitat. As a result, few nests were initiated, no nests hatched, and reproductive output at the site was 0 for a second consecutive year.
The two high-water years also created an abundance of nesting habitat for the M1F subpopulation downstream from the dam for the 2012 and 2013 breeding seasons (USACE, unpublished data). The Lewis and Clark Lake subpopulation was largely unaffected by the water level fluctuations (hereafter ‘M2,’ Missouri River site without high flow effects), although some new sandbars were created in 2012 and 2013 from the high velocity of the water passing through the marsh at the upstream end of the lake (USACE, unpublished data). The HC subpopulation was completely unaffected by these water fluctuations (MBB, pers. obs.).
We searched nesting areas at each subpopulation during the plover breeding season (April–August) from 2008 to 2013. At all sites, we visually scanned the area for both banded and unbanded adults and chicks, watched for behaviors suggestive of nesting or breeding, and searched for nests by walking through all potential nesting habitat (i.e., unvegetated and sparsely vegetated wet and dry sand habitat). We captured unbanded, incubating adult (age ≥ 1 yo) plovers using drop door traps placed over their nests, and juvenile (age 0–1 yo) birds were caught by hand as soon after hatching as possible. Once captured, birds were banded using a color band combination that was unique to the individual. We associated adults with nests if the bird was captured on a specific nest or if we observed them incubating eggs or brooding chicks.
Each individual was assigned to one of the subpopulations (HC, M1F, M2) for a given year based on sightings and information about nesting status. If an individual was sighted in multiple subpopulations within the same year, we first assigned birds to a subpopulation based on known breeding locations, or, if that information was unavailable, we assumed that the bird belonged to the subpopulation within which it was most frequently sighted in that year. If there were an equal number of sightings in multiple subpopulations, we ignored the sightings and assigned the bird a value of ‘0’ for ‘unseen in that year’ (only 11 of > 15,000 occasions). All juvenile birds were assigned to the subpopulation in which they hatched.
Survival and dispersal
To determine what factors affected plover survival and dispersal among our subpopulations, we analyzed capture histories for both juvenile and adult birds from 2008 to 2013 using multi-state mark-recapture models in Program MARK . Multi-state models allow for simultaneous estimation of apparent annual survival (φ), resight (p), and transition (ψ) rates for multiple ‘states’ or ‘strata’ [31, 32]. In our model, the states or strata were the subpopulations of our metapopulation, and transition rates represented the probability of dispersing among the subpopulations. We examined the effect of age class (adult vs. juvenile), year, subpopulation, hatch date, age at banding (in days), reproductive success, distance, and high flow on survival, transition, and resight rates (Additional file 1).
We used the median ĉ test in Program MARK to assess the goodness of fit of the general multi-state model. We estimated overdispersion (ĉ) using the most complex model, such that survival, resight, and transition rates all varied by age class, subpopulation (‘sub’), and year (i.e., age class × sub × year). For all models in this analysis, we fixed juvenile survival and emigration from M1F in 2010 and 2011 at 0 because no fledged chicks were produced during those years at M1F due to high water levels. We used ĉ to adjust standard errors and deviance estimates to account for overdispersion. We used Akaike’s Information Criterion (QAICc; corrected for small sample bias and overdispersion) to rank and to interpret our models .
We used a sequential approach for model selection to reduce the number of models under consideration  and to create a sufficiently predictive, general model against which we could test our ecological hypotheses about survival and dispersal among the subpopulations in the final step. At each subsequent step, we retained the model structure for parameters in the best model (lowest QAICc) for use in the next step of model selection (Additional file 1).
In the first step of modeling, we varied the model for resight rate (p) while holding survival (φ) and transition models (ψ) at their most complex (i.e., age class × sub × year). In addition, we added individual covariates for the effects of hatch date (‘hatch’) and the age at banding (in days, ‘age’) to all models for survival. We added these covariates to control for known sources of variation in our data [35–37]. Because of a difference in search frequency at each of the subpopulations, all models for resight rate controlled for differences among the subpopulations. We tested all possible combinations (both additive and multiplicative) of age class, year, and subpopulation (8 models; Additional file 1).
In the second step of the model selection procedure, we varied survival rate (φ) while holding transition rate at its most complex structure (i.e., age class × sub × year), and modeled resight rate as the best model (lowest QAICc) from the first step. We assumed a priori that juvenile survival would be different than adult survival  and controlled for this difference as well as the effects of hatch date and age at banding in all of our models. We tested all possible combinations (both additive and multiplicative) of age class, year, and subpopulation (8 models; Additional file 1).
In the third step of model selection, we varied the model for transition rate while modeling resight and survival rates as the best model from the first and second steps, respectively. We tested all possible combinations (both additive and multiplicative) of age class, year, and subpopulation, including a constant model (15 models; Additional file 1).
In the final step, we tested several hypotheses about survival and transition rates in relation to high flow (habitat availability), distance-mediated dispersal, and reproductive success. Because ‘subpopulation’ and ‘year’ are not necessarily biologically informative parameters, we hypothesized that yearly and subpopulation-specific variation in transition rates could be explained by biologically relevant parameters. For instance, reproductive success can be an important driver of breeding dispersal in many birds [38–40] and has been shown to affect plovers in particular . There is also some evidence that habitat availability (related to flooding of breeding sites) may affect plover dispersal . In addition, dispersal probabilities tend to decrease with increasing distance from a source, although the form of that decrease is often species-specific . We hypothesized that the distance among our sites would be negatively correlated with transition rates. We also hypothesized that site-specific annual reproductive success (average number of chicks fledged/pair at a site where reproductive success in year t affects transition rates between years t and t + 1; [15; MBB, KLH, unpublished data]) and changes in habitat availability would be more parsimonious descriptors of transition probabilities than year-specific and year-by-site specific transition rates.
To model these effects, we used each of the model structures for transition rate from the preceding step within Δ QAICc < 4 “points” of the top-ranked model. For each of these models, we replaced any occurrence of ‘subpopulation’ with the distance between each site. We replaced any occurrence of ‘year’ with reproductive success, variables for flow (see below), or both flow and reproductive success. Given the observed effects of high flows on M1F (i.e., loss of nesting habitat in 2011, increase in habitat in 2012 in the year after; ) and the lack of high flow effects observed for M2 and HC during the study period, we hypothesized that high flows would cause an immediate increase in dispersal out of M1F and a decrease in dispersal into M1F. In addition, we hypothesized the transition rates into M1F from elsewhere would increase for the year(s) following a high flow event as habitat quality and quantity increased.
To test these hypotheses, we created three variables. The first variable for high flow (‘high flow emigration’) affected transition from M1F to the other subpopulations from 2010 to 2011, and the second, high flow variable (‘high flow immigration’) affected transition from M2 or HC to M1F from 2010 to 2011. The variable for ‘post-high flow immigration’ affected any transition from HC or M2 to M1F (where new habitat was created) from 2011 to 2012. All three variables were used simultaneously in modeling transition rates. If a model contained a parameter for age class, we included the interaction between the reproductive success variable and age class under the assumption that juvenile and adult birds would respond differently to reproductive success.
In addition to replacing year and subpopulation for transition rate in this step, we added parameters for reproductive success and high flow to the top survival models to explore the potential effects of these factors on survival. These parameters did not replace subpopulation or year in survival models. We analyzed all model structure combinations for survival and transition rates (32 models; Additional file 1).
All real estimates were model-averaged over all models using QAICc weights . We used model ranking, t-statistics (β/SE), and 95 % confidence limits to interpret the relative size of individual effects (estimated βs) from top-ranked models.
We used the demographic and transition rates observed in this study to parameterize a population viability analysis (PVA) model for the plover metapopulation formed by the HC, M2, and M1F subpopulations. We used this PVA model to investigate long-term metapopulation viability given the current metapopulation structure, plover demographic rates, and local disturbance regime. The PVA model was constructed in Vortex (version 10.0.7.3; ), a widely used, previously validated  program that simulates the effects of both deterministic forces and demographic, environmental, and genetic stochastic events to assess extinction risk.
Estimate of fledged piping plover chicks produced/paira from 2008–2012
Resighting rates of banded birds varied by age class, subpopulation, and year. Rates were highest at M1F (adult: 0.92 ± 0.06, juvenile: 0.74 ± 0.13; mean ± SD) and M2 (adult: 0.88 ± 0.08, juvenile: 0.66 ± 0.14; mean ± SD) relative to those on HC (adult: 0.64 ± 0.14, juvenile: 0.32 ± 0.11; mean ± SD).
Multistate mark recapture model ranking for piping plover survival and transition from 2008–2012
φ age×sub + age×year + juvenile:(band + hatch) + high flow ψ age + distance + high flow b
φ age×sub + age×year + juvenile:(band + hatch) + high flow ψ age×success + distance + high flow
φ age×sub + age×year + juvenile:(band + hatch) + high flow + adult:success ψ age + distance + high flow
φ age×sub + age×year + juvenile:(band + hatch) ψ age + distance + high flow
φ age×sub + age×year + juvenile:(band + hatch) + adult:success ψ age + distance + high flow
φ age×sub + age×year + juvenile:(band + hatch) ψ age×success + distance + high flow
φ age×sub + age×year + juvenile:(band + hatch) + high flow + adult:success ψ age×success + distance + high flow
φ age×sub + age×year + juvenile:(band + hatch) + adult:success ψ age×success + distance + high flow
Beta coefficients of the effect of group- and individual-specific covariates on survival and transition parameters from the top-ranked multi-state mark recapture model
Lower 95 % CL
Upper 95 % CL
Hatch Date (juvenile)
Age at Banding (juvenile)
High flow Emigration
High flow Immigration
Post-high flow Immigration
There was little indication from model ranking that reproductive success affected survival (Table 2). The highest ranked model that included reproductive success had a weight of only 0.105 (Table 2), and the confidence limits for the coefficient included 0, indicating that the effect was not precisely estimated (Table 3). Birds that hatched later in the year had lower survival than those hatching earlier. Age at banding had a positive effect on juvenile survival, but the 95 % CI and magnitude of the β/SE value for age at banding indicated that this variable had only a marginal effect relative to the other covariates in the top-ranked model (Table 3).
Biologically relevant parameters (i.e., age-class, distance between sites, and flooding) were more parsimonious (lower QAICc values) descriptors of plover dispersal than the subpopulation and year parameters (Table 2). Neither subpopulation nor year appeared in the competing (QAICc ≤ 4) transition rate models. The highest ranked models containing subpopulation were > 10 QAICc “units” from the top-ranked model and had essentially 0 weight. The highest ranked model containing year was > 32 QAICc points from the top-ranked model.
High flows increased transitions from M1F (‘high flow emigrate’) and decreased transitions into M1F (‘high flow immigrate’), and the post-high flow environment was associated with increased immigration rates into M1F (‘post-high flow immigrate’; Table 3). The inclusion of reproductive output in the transition rate model did not substantially improve fit (top-ranked model vs. second ranked model; Table 2). The interaction between reproductive success and age class appeared in the second-ranked model but did not decrease the deviance enough to outweigh the penalty (Table 2).
Extinction risk for the piping plover metapopulation
Extinction probability (SE)
Mean years to extinction (SE)
Mean subopulation size (individuals)
0.0 % (0.0)
0.0 % (0.0)
47.7 % (1.6)
0.5 % (0.2)
The overall low risk of metapopulation extinction was due largely to the persistence of HC. Adults in HC comprised 60.6 % of the total metapopulation by year 100, compared to only 11.0 % from M1F despite the fact that this subpopulation had a substantially larger carrying capacity in most years (Table 4; Fig. 6b). Furthermore, M1F had a deterministic finite rate of increase (λdeter) of 0.9, indicating that, based on a life table analysis of the mean mortality and reproductive rates initially used to parameterize the model, the population should decline. However, the stochastic finite rate of population increase (λstoch), which also takes into account stochastic fluctuations and immigration/emigration throughout the course of a simulation, was 1.1 (Table 4). This result indicates that M2 and HC supported the persistence of M1F given the frequency of high flow events that currently occurs along the Missouri River.
The metapopulation dynamics in this study were strongly influenced by landscape dynamics, with functional connectivity (as measured through dispersal rates) changing through time as a result of both natural (high flows) and anthropogenic (isolation by distance) factors. It is clear that extreme events, including natural disturbances, can profoundly influence a species’ evolutionary history and population dynamics [6, 43, 44]. Large floods and bankfull flows historically occurred regularly on the Missouri River and its tributaries and can exert a primary source of selective pressure for adaptation as a cause of mortality in species like piping plovers . Such species can exhibit life history adaptations (e.g., synchronization of life history events in relation to a common flow regime) and/or behavioral adaptations (e.g., adaptations that allow a species to respond directly to high flow events) in response to the magnitude, frequency, seasonal timing, predictability, and/or duration of these events . This study and others have identified several such potential adaptations in plovers in the Great Plains. For instance, before the prevalence of dams, the Missouri River and surrounding tributaries experienced relatively predictable, high water flow peaks in March and late May/early June . Plovers in this region typically begin laying eggs in mid- to late May, which, in addition to avoiding much of the inclement weather in the spring and producing chicks near peak invertebrate abundance in mid- to late-summer, allowed them to produce offspring after waters had receded to expose high quality habitat . Plovers may also exhibit behavioral adaptations that allow them to tolerate high flows, such as their ability to renest up to four times in a breeding season in the event that eggs or young hatchlings are lost to flooding or predation . Such life history adaptations are expected for species in dynamic environments where disturbances are frequent, large, and predictable , characteristic of the historical Missouri River .
Our results also indicate that plovers increase their dispersal rates and reproductive output following high flow events (or large-scale fluctuations in the amount of habitat). Species in naturally disturbed environments generally display compensatory or stabilizing effects to counterbalance higher mortality during a disturbance . In our study, adult mortality increased slightly while hatch-year mortality was 100 % within M1F as a result of high flows, and identical trends have been observed in a piping plover population on the Platte River  as well as in other waterbird species . At M1F, higher flow-related mortality in 2010–2011 was compensated for in 2012 and 2013 by lower mortality and higher reproductive output, likely due to density-dependent increases in chick survival (KLH, unpublished data) and nesting  related to high flows in this system . Higher recruitment following high flow-related mortality has also been demonstrated in crimson finches (Neochmia phaeton; ) and some Australian waterbirds . Therefore, these compensatory mechanisms (i.e., low mortality and high reproductive output following high flows) could be important adaptations in plovers and other riverine species that promote population persistence in response to high flows, which would have historically occurred multiple times within a plover’s lifespan in the Great Plains .
Piping plovers generally exhibit high site fidelity , which could benefit a species in somewhat unpredictable, dynamic environments . Our observations in normal flow years, showing lower dispersal rates between all pairs of populations, further support the propensity for site fidelity in this species. In contrast, we found dispersal rates from M1F increased substantially during the 2010–2011 high flow years as some birds left inundated habitat and moved to the HC and M2 populations, while dispersal rates from M2 and HC into M1F increased significantly after 2011 when high quality habitat was created. High flows and other extreme events appear to be important dispersal cues for many otherwise site faithful species, including several birds and fish [44, 49, 54–57]. Increased dispersal rates into M1F were critical to the persistence and recolonization of this subpopulation. Therefore, this study confirms that natural disturbance is a strong driver of dispersal in this system, and disturbance-related shifts in functional connectivity had a major influence on metapopulation dynamics and the persistence of individual subpopulations.
Although high flows clearly influenced plover dispersal rates, the precise cues used by this species to make dispersal decisions remain unclear. Individuals could have responded to some regional abiotic cue associated with the high flow event itself (). Alternatively, individuals could have reacted directly to the amount of nesting habitat or to the lower population size within M1F after the event; immigration rates in plovers and other birds can be density dependent [14, 39, 58]. Dispersal in plovers was previously linked to reproductive success, exhibiting a positive relationship between emigration rates and nest failure at a given site . The relationship between individual reproductive success and dispersal is well established in birds [38–40], and there is growing evidence that individuals use information on local, conspecific reproductive success to make dispersal decisions [58–61]. Although the evidence that reproductive output within a location was correlated with dispersal in this study was weak, increased emigration out of M1F following total reproductive failure could be related to reproductive success as well as habitat loss. Irrespective of the root cause, high flows along the Missouri River were closely related to increased functional connectivity in this plover metapopulation.
If plovers have adapted certain life history attributes and behaviors to the natural flow regime, relatively recent alterations to this disturbance regime may have far-reaching effects on the species’ population dynamics and viability [3, 45]. Dams along the Missouri River and its tributaries have almost completely eliminated high flows, and only the most catastrophic events (e.g., the 2010–2011 event observed in this study) have the capacity to exceed dam storage capabilities . In the absence of habitat-creating flows, we would expect productivity and dispersal rates to remain at low levels for extended periods of time, consistent with rates observed prior to the 2010–2011 high flow events . Furthermore, nearly 70 % of the Missouri River is either impounded or channelized , which has increased the distance among areas with suitable habitat. Plovers will use reservoir shorelines to nest, but many of the reservoirs do not have shorelines suitable for nesting . If connectivity is depressed in the absence of flow-based cues or because of increased isolation from channelization and impoundments, the benefits of dispersal may also decline for this plover metapopulation. The spatial structure and connectivity patterns within metapopulations can govern critical genetic and demographic processes . Decreased connectivity could reduce the likelihood of the recolonization of an extirpated population  or eliminate demographic rescue effects for a declining population . A reduction in dispersal also could increase the risks of a population experiencing inbreeding depression or genetic drift , particularly given the small sizes of these plover subpopulations. Therefore, a loss of flow-based dispersal cues could have important implications for this, and other, piping plover metapopulations.
The implications of the loss of high flows for this species, which has various adaptations to a natural flow regime, can be seen in the results of the PVA and a related study that investigated multiple management and flow scenarios . In this study, M1F had a relatively high extinction risk (47.7 %) under the current flow regime and was predicted to support only a very small total population (22 adults). Twenty-two adults is exactly the number observed in 1996 and 1997 at M1F, the nadir during current monitoring and before a substantial high water event. These rates translated to a total predicted average metapopulation size of approximately 200 adults, far below the average expected carrying capacity of 700 adults. In contrast, separate simulations predicted metapopulation size would be substantially higher under a more frequent flow regime, reminiscent of a flow regime to which plovers in the Great Plains were originally adapted .
Conservation-reliant species require some level of management to prevent extinction or extirpation . Piping plovers will likely require varying levels of habitat management, beach closure, predator control, monitoring, and captive rearing throughout their range to maintain and to grow populations [15, 68]. In our study, the predicted persistence of the HC and M2 populations is dependent upon the assumption of consistent management and availability of nesting habitat at these sites to maintain the carrying capacity and reproductive output. This result may have been related to the lack of correlation of disturbances among subpopulations; increased autocorrelation in habitat turnover can negatively impact metapopulation persistence .
Human-created sites that support HC and M2 are critical to the survival of the metapopulation, by providing both stable habitat and dispersers. Although habitats used by these populations were relatively stable during our study period and were modeled as such, habitat was created through very different processes. Habitat used by HC is created by sand and gravel mining operations and residential housing development construction. Mining practices involve mechanical disturbance of dredged sand, which maintains sparsely vegetated expanses of sand adjacent to water . The conversion of mines to housing developments requires additional mechanical disturbance to redistribute the sand so that the topography is suitable for housing construction.
Sandbars used by plovers within M2 during this study were created by the USACE in order to comply with requirements of the Endangered Species Act [16, 35]. These human-created sandbars require vegetation management, sand augmentation, and predator control  to maintain their usefulness beyond a few years. Management of human-created HC habitat is based on agreements between private interests (mining companies or real estate developers) and regulatory agencies and is limited mostly to implementing measures that avoid “take”. Thus, while the USACE is incentivized to create and maintain habitat in M2 to meet regulatory obligations, private interests in the areas used by HC are not . Consequently, evolving mining practices and economic conditions could appreciably alter the amount, configuration and distribution of habitat used at HC.
Our analysis of this controlled, natural experiment (albeit with low replication) indicated that despite a low overall probability of metapopulation extinction, the persistence was predicated on temporarily stable, anthropogenically modified habitats as well as dispersal among these subpopulations. Given the highly dynamic nature of habitat availability in this system, maintaining several populations within the metapopulation, stable sources of habitat, and assuring connectivity through dispersal may be critical.
Funding was provided by the Nebraska Environmental Trust, the Nebraska State Wildlife Grant Program, the Nebraska Wildlife Conservation Fund, the U.S. Army Corps of Engineers, the U.S. Fish and Wildlife Service, and Virginia Tech. We thank cooperators, the National Park Service, U.S. Fish and Wildlife Service, South Dakota Department of Game, Fish, and Parks, Nebraska Game and Parks Commission, and the Missouri River Institute for support throughout the project. We thank the sand and gravel mining companies and lakeshore housing development communities for access to property. We acknowledge the tireless efforts of our many technicians from 2005–2013. This work was conducted under Institutional Animal Care and Use Committee permits 877 and 14-003, and U.S. Fish and Wildlife Service Threatened and Endangered Species permits TE070027-0 and TE103272-3. Finally, E. Grant and four anonymous reviewers provided helpful comments that were used to improve this manuscript.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Hanski I. Metapopulation dynamics. Nature. 1998;396:41–9.View ArticleGoogle Scholar
- Akçakaya H, Mills G, Doncaster C. The role of metapopulations in conservation. In: MacDonald D, Service K, editors. Key Topics in Conservation Biology. Malden: Blackwell Publishing; 2007. p. 64–84.Google Scholar
- van Teefelen A, Vos C, Opdam P. Species in a dynamic world: consequences of habitat network dynamics on conservation planning. Biol Conserv. 2012;153:239–53.View ArticleGoogle Scholar
- Moloney K, Levin S. The effects of disturbance architecture on landscape-level population dynamics. Ecology. 1996;77:375–94.View ArticleGoogle Scholar
- Thomas C. Extinction, colonization and metapopulations: environmental tracking by rare species. Conserv Biol. 1994;8:373–8.View ArticleGoogle Scholar
- Stelter C, Reich M, Grimm V, Wissel C. Modelling persistence in dynamic landscapes: lessons from a metapopulation of the grasshopper Bryodema tuberculata. J Anim Ecol. 1997;66(4):508–18.View ArticleGoogle Scholar
- Meulebrouk K, Verheyen K, Brys R, Hermy M. Metapopulation viability of an endangered holoparasitic plant in a dynamic landscape. Ecography. 2009;32:1040–50.View ArticleGoogle Scholar
- Rappaport D, Tambosi L, Metzger J. A landscape triage approach: combining spatial and temporal dynamics to prioritize restoration and conservation. J Appl Ecol. 2015;52(3):590–601.View ArticleGoogle Scholar
- BenDor T, Woodruff S. Moving targets and biodiversity offsets for endangered species habitat: is lesser prairie chicken habitat a stock or flow? Sustainability. 2014;6:1250–9.View ArticleGoogle Scholar
- Zeigler S, Fagan WF. Transient windows for connectivity in a changing world. Movement Ecol. 2014;2:1.View ArticleGoogle Scholar
- Brachet S, Olivieri I, Godelle B, Klein E, Franscaria-Lacoste N, Gouyon P-H. Dispersal and metapopulation viability in a heterogeneous landscape. J Theorretical Biol. 1999;198:479–95.View ArticleGoogle Scholar
- Johst K, Brandl R, Eber S. Metapopulation persistence in dynamic landscapes: the role of dispersal distance. Oikos. 2002;98:263–70.View ArticleGoogle Scholar
- Elliot-Smith E, Haig S. Piping plover (Charadrius melodus). In: Birds of North America Online. Cornell Lab of Ornithology, Ithaca, NY. 2004. http://bna.birds.cornell.edu/bna/species/002. Accessed 10 Oct 2014.
- Catlin D, Fraser J, Felio J. Demographic responses of piping plovers to habitat creation on the Missouri River. Wildl Monogr. 2015;192:1–42.View ArticleGoogle Scholar
- U.S. Fish and Wildlife Service. Piping Plover (Charadrius melodus) 5-Year Review: US Fish and Wildlife Service2009.Google Scholar
- Catlin D, Fraser J, Felio J, Cohen J. Piping plover habitat selection and nest success on natural, managed, and engineered sandbars. J Wildl Manag. 2011;75(2):305–10.View ArticleGoogle Scholar
- Bomberger Brown M, Jorgensen J, Dinan L. Interior least tern and piping plover monitoring, research, management, and outreach report for the Lower Platter River, Nebraska. Joint Report of the Tern and Plover Conservation Partnership and the Nebraska Game and Parks Commission. Lincoln: Tern and Plover Conservation Partnership and Nebraska Game and Parks Commission; 2013.Google Scholar
- McGowan C, Catlin D, Shaffer T, Gratto-Trevor C, Aron C. Establishing endangered species recovery criteria using predictive simulation modeling. Biol Conserv. 2014;177:220–9.View ArticleGoogle Scholar
- Anteau M, Shaffer T, Sherfy M, Sovada M, Stucker J, Wiltermuth M. Nest survival of piping plovers at a dynamic reservoir indicates an ecological trap for a threatened population. Oecologia. 2012;170:1167–79.View ArticleGoogle Scholar
- Roche E, Gratto-Trevor C, Goossen J, White C. Flooding effects dispersal decisions in piping plovers (Charadrius melodus) in prairie Canada. Auk. 2012;129(2):296–306.View ArticleGoogle Scholar
- Cohen J, Houghton L, Fraser J. Nesting density and reproductive success of piping plovers in response to storm- and human-created habitat changes. Wildl Monogr. 2009;173(1):1–24.View ArticleGoogle Scholar
- Wilcox L. A twenty-year banding study of the piping plover. Auk. 1959;76(2):129–52.View ArticleGoogle Scholar
- Catlin D, Jacobson R, Sherfy M, Anteau M, Felio J, Fraser J, et al. Discussion of “Natural hydrograph of the Missouri River near Sioux City and the least tern and piping plover” by Donald Jorgensen. J Hyrdol Eng. 2010;15:1076–8.Google Scholar
- Hesse L, Mestle G. An alternative hydrograph for the Missouri River based on the precontrol condition. N Am J Fish Manag. 1993;13:360–6.View ArticleGoogle Scholar
- Dunne T, Leopold L. Water in Environmental Planning. New York: W.H. Freeman and Company; 1978.Google Scholar
- Joeckel R, Henebry G. Channel and island change in the lower Platte River, eastern Nebraska USA. Geomorphology. 2008;102(3):407–18.View ArticleGoogle Scholar
- Johnson W. Woodland expansions in the Platte River, Nebraska: patterns and causes. Ecol Monogr. 1994;64(1):45–84.View ArticleGoogle Scholar
- Laustrup M, LeValley M, editors. Missouri River Environmental Assessment Program. Missouri River Natural Resouces Committee, US Geological Survey; 1996.Google Scholar
- Dixon M, Johnson W, Scott M, Bowen D, Rabbe L. Dynamics of plains cottonwood (Populus deltoides) forests and historical landscape change along unchannelized segments of the Missouri River. USA Environ Manag. 2012;49:990–1008.View ArticleGoogle Scholar
- White G, Burnham K. Program MARK: survival estimation from populations of marked animals. Bird Study. 1999;46(Supplement):120–38.View ArticleGoogle Scholar
- Brownie C, Hines J, Nichols J, Pollock K, Hestbeck J. Capture-recapture studies for multiple strata including non-Markovian transitions. Biometrics. 1993;49:1173–87.View ArticleGoogle Scholar
- Hestbeck J, Nichols J, Malecki R. Estimates of movement and site fidelity using mark resight data of wintering Canada geese. Ecology. 1991;72:523–33.View ArticleGoogle Scholar
- Burnham K, Anderson D. Model selection and multi-model inference: a practical information-theoretic approach (2nd edition). New York: Springer-Verlag; 2002.Google Scholar
- Lebreton J, Burnham K, Clobert J, Anderson D. Modeling survival and testing biological hypotheses using marked animals - a unified approach with case-studies. Ecol Monogr. 1992;62:67–118.View ArticleGoogle Scholar
- Catlin D. Population dynamics of piping plovers (Charadrius melodus) on the Missouri River. Blacksburg VA: Virginia Tech; 2009.Google Scholar
- Catlin D, Felio J, Fraser J. Effect of great horned owl trapping on chick survival in piping plovers. J Wildl Manag. 2011;75:458–62.View ArticleGoogle Scholar
- Catlin D, Felio J, Fraser J. Effects of water discharge on fledging time, growth, and survival of piping plovers on the Missouri River. J Wildl Manag. 2013;77(3):525–33.View ArticleGoogle Scholar
- Catlin D, Rosenberg D, Haley K. The effects of nesting success and mate fidelity on breeding dispersal in burrowing owls. Can J Zool-Revue Canadienne de Zoologie. 2005;83:1574–80.View ArticleGoogle Scholar
- Greenwood PJ, Harvey PH. The natal and breeding dispersal of birds. Annu Rev Ecol Syst. 1982;13:1–21.View ArticleGoogle Scholar
- Haas C. Effect of prior nesting success on site fidelity and breeding dispersal: an experimental approach. Auk. 1998;115:929–36.View ArticleGoogle Scholar
- Lacy R, Borbat M, Pollack J. Vortex: A Stochastic Simulation of the Extinction Process (version 9.99). Brookfield: Chicago Zoological Society; 2010.Google Scholar
- Brook B. Pessimistic and optimistic bias in population viability analysis. Conserv Biol. 2000;14(2):564–6.View ArticleGoogle Scholar
- Parmesan C, Root T, Willig M. Impacts of extreme weather and climate on terrestrial biota. Bull Am Meteorol Soc. 2000;81(3):443–50.View ArticleGoogle Scholar
- Santoro S, Green A, Figerola J. Environmental instability as a motor for dispersal: a case study from a growing population of glossy ibis. PLoS One. 2013;8(12):e82983.View ArticleGoogle Scholar
- Lytle D, Poff N. Adaptation to natural flow regimes. Trends Ecol Evol. 2004;19(2):94–100.View ArticleGoogle Scholar
- U.S. Fish and Wildlife Service. Piping Plover (Charadrius melodus), Great Plains population, Draft Revised Recovery Plan: U.S. Fish and Wildlife Service2000.Google Scholar
- Lloret F, Escudero A, Iriondo J, Martínez - Vilalta J, Valladares F. Exterme climatic events and vegetation: the role of stabilizing processes. Global Change Biology. 2012; 18(3):797-805.Google Scholar
- Sidle J, Carlson D, Kirsch E, Dinan J. Flooding: mortality and habitat renewal for least terns and piping plovers. Colonial Waterbirds. 1992;15(1):132–6.View ArticleGoogle Scholar
- Poiani A. Effects of floods on distribution and reproduction of aquatic birds. Adv Ecol Res. 2006;39:63–83.View ArticleGoogle Scholar
- Hunt KL, Dinan LR, Friedrich MJ, Bomberger Brown M, Jorgensen JG, Catlin DH, et al. Density dependent double brooding in piping plovers (Charadrius melodus) in the northern Great Plains USA. Waterbirds. 2015;38(4):321–9.View ArticleGoogle Scholar
- Houston W, Black R. Effect of major flooding on breeding and habitat of the crimson finch (Neochmia phaeton): a riparian specialist. River Res Appl. 2014;30:609–16.View ArticleGoogle Scholar
- Friedrich M, Hunt K, Catlin D, Fraser J. The importance of site to mate choice: mate and site fidelity in piping plovers. Auk. 2015;132(1):265–76.View ArticleGoogle Scholar
- Switzer P. Site fidelity in predictable and unpredictable habitats. Evol Ecol. 1993;7:533–55.View ArticleGoogle Scholar
- David B, Closs G. Behavior of a stream-dwelling fish before, during and after high-discharge events. Trans Am Fish Soc. 2012;131:762–71.View ArticleGoogle Scholar
- Dugger K, Ainley D, Lyver P, Barton K, Ballard G. Survival differences and the effect of environmental instability on breeding dispersal in an Adélie penguin metapopulation. Proc Natl Acad Sci. 2010;107(27):12375–80.View ArticleGoogle Scholar
- Koster W, Cook D. Diurnal and nocturnal movements of ricer blackfish (Gadopsis marmoratus) in a south-eastern Australian upland stream. Ecol Freshw Fish. 2008;17:146–54.View ArticleGoogle Scholar
- Nunn A, Copp G, Vilizzi L, Carter M. Seasonal and diel patterns in the migrations of fishes between a river and a floodplain tributary. Ecol Freshw Fish. 2010;19:153–62.View ArticleGoogle Scholar
- Hénaux V, Bregnballe T, Lebreton J. Dispersal and recruitment during population growth in a colonial bird, the great comorant Phalacrocorax carbo sinensis. J Avian Biol. 2007;38:44–57.View ArticleGoogle Scholar
- Brown C, Bomberger Brown M, Danchin E. Breeding habitat selection in cliff swallows: The effect of conspecific reproductive success on colony choice. J Anim Ecol. 2000;69:133–42.View ArticleGoogle Scholar
- Danchin E, Boulinier T, Massot M. Conspecific reproductive success and breeding habitat selection: Implications for the study of coloniality. Ecology. 1998;79(7):2415–28.View ArticleGoogle Scholar
- Rioux S, Amirault-Langlais D, Shaffer F. Piping plover make decisions regarding dispersal based on personal and public information in a variable coastal ecosystem. J Field Ornithol. 2011;82:32–43.View ArticleGoogle Scholar
- U.S. Fish and Wildlife Service. Amendment to the 2000 biological opinion on the operation of the Missouri River main stem reservoir system, operation and maintenance of the Missouri River bank stabilization and navagation project, and operation of the Kansas River reservoir system. Lakewood: U.S. Fish and Wildlife Service; 2003.Google Scholar
- Henderson M, Merriam G, Wegner J. Patchy environments and species survival: chipmunks in an agricultural mosaic. Biol Conserv. 1985;31:95–105.View ArticleGoogle Scholar
- Brown J, Kodric-Brown A. Turnover rates in insular biogeography: effect of immigration on extinction. Ecology. 1977;58:445–9.View ArticleGoogle Scholar
- Richards C. Inbreeding depression and genetic rescue in a plant metapopulation. Am Nat. 2000;155(3):383–94.View ArticleGoogle Scholar
- Zeigler S, Catlin D, Bomberger Brown M, Fraser J, Dinan L, Hunt K et al. A brave new world: metapopulation dynamics of piping plovers following the loss of natural disturbances. Ecological Applications. in revision.Google Scholar
- Scott J, Goble D, Haines A, Wiens J, Neel M. Conservation-reliant species and the future of conservation. Conserv Lett. 2010;3:91–7.View ArticleGoogle Scholar
- Hecht A, Melvin S. Expenditures and effort associated with recovery of breeding Atlantic Coast piping plovers. J Wildl Manag. 2009;73(7):1099–107.View ArticleGoogle Scholar