Quantifying multiple breeding vital rates in two declining grassland songbirds

. Many studies of reproductive success in North American songbirds have focused on nesting success, while relatively few have evaluated breeding-season adult survival and post-fledging survival. Grassland songbirds are among North America's most rapidly declining avian groups, and knowledge of factors that influence vital rates is needed to address declines, develop management strategies, and accurately model population limitation. We concurrently monitored nesting success, breeding-season adult survival, and postfledging survival of two grassland obligates, Baird's Sparrow and Grasshopper sparrow, breeding in western North Dakota and northeastern Montana. Nesting success was monitored by locating and visiting nests at regular intervals while adult and post-fledging survival were assessed by daily telemetry tracking of radio-tagged birds. We analyzed the three variables using logistic exposure and modeled climate, temporal, and vegetative covariates to explain variation in rates. Cumulative nesting success, breeding-season adult survival, and post-fledging survival were 37%, 78%, and 25%, respectively, for Baird's Sparrow and 16%, 74%, and 55% for Grasshopper Sparrow. Both nesting success and post-fledging survival in Baird's Sparrow were responsive to environmental covariates including temporal effects and vertical vegetation structure. Conversely, vital rates of Grasshopper Sparrow were largely unresponsive to covariates we modeled, perhaps because of the species' broader habitat niche relative to Baird's Sparrow. Breeding season adult survival in both species showed little annual variation and was high relative to overwintering survival estimates for the same species, while post-fledging survival in Baird's Sparrow was low and may be a management concern. We suggest as a next step the formal comparison of vital rates across life-stages in an integrated population model capable of identifying sources of population limitation throughout the full annual cycle of the species. que la survie des adultes et la survie après l'envol ont été calculées au moyen du suivi télémétrique quotidien d'oiseaux munis d'un émetteur. Nous avons analysé ces trois variables au moyen de l'exposition logistique et modélisé les covariables climatiques, temporelles et végétales pour expliquer la variation des taux. Le succès cumulatif de la nidification, la survie des adultes en saison de reproduction et la survie après l'envol étaient respectivement de 37 %, 78 % et 25 % pour le Bruant de Baird, et de 16 %, 74 % et 55 % pour le Bruant sauterelle. Le succès de nidification et la survie après l'envol du Bruant de Baird étaient sensibles aux covariables environnementales, y compris les effets temporels et la structure verticale de la végétation. À l'inverse, les taux vitaux du Bruant sauterelle ne variaient pas en fonction des covariables que nous avons modélisées, peut-être en raison de la niche d'habitat plus élargie de cette espèce comparativement au Bruant de Baird. La survie des adultes en saison de reproduction chez les deux espèces a montré peu de variation annuelle et était élevée par rapport à leur taux de survie en hiver, tandis que la survie après l'envol chez le Bruant de Baird était faible et devrait peut-être être une préoccupation à prioriser. Comme prochaine étape, nous proposons de comparer de façon formelle les taux vitaux à travers les étapes de vie dans un modèle de population intégré permettant d'identifier les facteurs qui limitent les populations tout au long du cycle annuel complet de ces deux espèces.


INTRODUCTION
Accurate demographic rates are fundamental to our understanding of species ecology (Pulliam 1988, Murdoch 1994 particularly when reversing population decline is a management goal (Anders and Marshall 2005). Avian population dynamics are influenced by multiple vital rates and life-history traits (Clark and Martin 2007). These include not only those that are directly tied to fecundity, such as nesting success and clutch-size, but also rates that affect recruitment to death ratios and lifetime reproductive success, such as adult and juvenile survival (Stahl andOli 2006, Clark andMartin 2007). Vital rates may display both spatial and temporal variation, changing across environmental conditions and throughout the annual cycle of species (Rushing et al. 2017). Therefore, addressing species declines often requires a complete understanding of vital rates that influence population dynamics (Fletcher et al. 2006), and in some cases, failure to assess multiple rates can result in the overlooking of life-history phases that are critical for management (Crouse et al. 1987). Although species abundance data can be effectively used to track population changes over time (e.g., Rosenberg et al. 2019), field-collected data describing vital rate performance under different habitat conditions, and at different life-stages, may be needed to understand population changes mechanistically (Donovan et al. 1995, McCoy et al. 1999, Eng et al. 2011).
Among North American songbirds, breeding-season reproductive studies often focus on nesting success, ignoring adult and postfledging survival (Streby and Andersen 2011). This oversight is likely due in part to the difficulty and cost associated with tracking individual birds within and across seasons (Kershner et al. 2004, Suedkamp Wells et al. 2007, Rush and Stutchbury 2008, Cox et al. 2014). Yet, in populations of avian species, both adult and post-fledging survival can influence population dynamics (Pulliam et al. 1992, Anders and Marshall 2005, Fletcher et al. 2006, Bonnot et al. 2011, Cox et al. 2014) and the importance of these vital rates often vary among species (Stahl and Oli 2006). Further, habitat requirements can differ among life-history stages, thus managing for a single vital rate may be inadequate. For example, good foraging habitat for adults may differ from highquality nest sites (Steele 1993) while fledglings may select habitat cover that differs from adult preferences or changes with age Bock 2005, Small et al. 2015). If habitat preference is linked to fitness in any respect (Chalfoun and Schmidt 2012), it is reasonable to expect that response to environmental conditions among vital rates may differ accordingly. Therefore, there is a need to examine the effect of environmental covariates on multiple vital rates within populations both to guide management (Perlut et al. 2008a, Young et al. 2019 and to provide baseline data for the construction of accurate population models (Streby and Andersen 2011).
North America's grassland obligate songbirds provide a relevant model system to explore the importance of monitoring multiple vital rates, as many species within this group are both steeply declining and have been poorly studied with respect to adult and juvenile survival. Despite increased attention in recent decades, grassland songbirds and their habitats remain in crisis (Askins et al. 2007, North American Bird Conservation Initiative 2009, Green et al. 2018). This assemblage of species has experienced steep declines documented from the advent of the Breeding Bird Survey (Rosenberg et al. 2019) and faces a wide variety of threats on both the breeding and wintering grounds. These threats include habitat loss through conversion to agriculture, climate change, fragmentation and disturbance associated with energy development, shrub encroachment, non-native plant species, and disruption of historic fire and grazing regimes (see Brennan and Kuvlesky 2005, Askins et al. 2007, North American Bird Conservation Initiative 2009. Although declines among grassland birds are broadly understood to be driven by agricultural land conversion (Murphy 2003, Pool et al. 2014, Hill et al. 2014), more information on vital rates is required to develop effective, species-level, management strategies for implementation in remaining grassland habitats (e.g., Davis 2003, Fletcher et al. 2006, Perlut et al. 2008a).
The mixed-grass prairie of the Northern Great Plains (NGP) is one of the largest and most intact areas of grassland in North America, though today only 50% of its historic 600,000 km 2 extent remains (Comer et al. 2018). Still, this eco-region is critical breeding habitat for an assemblage of 26 grassland obligate bird species (Askins et al. 2007). We focused our demographic research on two species belonging to this group: Baird's Sparrow (Centronyx bairdii) and Grasshopper Sparrow (Ammodramus savannarum). Both species have suffered population losses of more than 60% since 1970 ) and both are species of management concern for multiple states and provinces in the Great Plains region (Green et al. 2018). The two species share similar habitat preferences and reproductive ecology in the NGP (Jones et al. 2010, Lipsey andNaugle 2017) and both are shortdistance migrants that winter in the desert grasslands of the southwestern United States and Northern Mexico (Vickery 1996, Green et al. 2002. However, despite apparent ecological similarities, the two species also present an interesting contrast with respect to degree of habitat specialism (Correll et al. 2019); Baird's Sparrow is narrowly range-restricted to the NGP (however see Youngberg et al. 2020) and highly specialized to mixed-grass prairie ecosystems during the breeding season (see Green et al. 2002) while Grasshopper Sparrow is continently distributed and found in a greater variety of surrogate grasslands and opencountry habitats (Vickery 1996).
Although the nesting ecologies of these species have been wellstudied (e.g., Davis and Sealy 1998, Davis 2003, Jones et al. 2010, Davis et al. 2016), estimates of nesting success in grassland species can vary with climate conditions (George et al. 1992, Skagen and Adams 2012, Conrey et al. 2016, Zuckerberg et al. 2018) as well as habitat and management treatments (e.g., Davis 2005, Lloyd and Martin 2005, Hovick et al. 2012, C. A. Davis et al. 2016, Pipher et al. 2016. Replicate studies, and those that examine both annual and within season variation, therefore remain valuable. Moreover, estimates of both adult and post-fledging survival are few, or lacking entirely, for many grassland birds. There are currently no published estimates of within breeding-season adult survival for Baird's Sparrow or Grasshopper Sparrow. Several studies of annual adult survival in Grasshopper Sparrow have been conducted in eastern populations (Perkins andVickery 2001, Balent andNorment 2003), but these rely on mark-recapture techniques and thus survival estimates are considered apparent because mortality cannot be distinguished from emigration (Lebreton et al. 1992). Research on post-fledging survival is also limited, with no estimates existing for Baird's Sparrow and only a single study of post-fledging survival in Grasshopper Sparrow (Hovick et al. 2011). Only one study to our knowledge has attempted to address the nesting, adult, and post-fledging lifestages simultaneously in a grassland bird species (van Vliet et al. 2020).
To address these knowledge gaps, we conducted a multi-year study to jointly examine the rates and drivers of nesting success (4 years), breeding-season adult survival (3 years), and post-fledging survival (3 years) in single populations of Baird's and Grasshopper Sparrows. We defined adult breeding-season survival as the probability of surviving for 90 days on the breeding grounds and post-fledging survival as the probability of surviving for 20 days, the approximate age of independence. We conducted nest monitoring and radio-telemetry tracking of adult and fledgling birds at two sites in the NGP, one in western North Dakota (2015Dakota ( -2018 and the other in northeastern Montana (2016Montana ( -2018. Broadly, we hypothesized that there would be variation in rates and drivers among life-history phases, as there are potential differences in predation pressure, vegetation structural requirements, and sensitivity to climate exposure across these stages (Jones and Bock 2005, Low et al. 2010, McCauley et al. 2017, Zuckerberg et al. 2018. We predicted that, overall, greater within-season precipitation would result in more productive range conditions (Barnett and Facey 2016), and thus higher survival across all stages (Conrey et al. 2016), despite some negative effects of extreme precipitation events (Carver et al. 2017). We also predicted that nests of both species would benefit from taller and more dense vegetation cover, given the species' habitat-preferences in the region (Lipsey and Naugle 2017) and potential benefits of concealment. Similarly, we predicted that fledglings of both species would benefit from taller and denser vegetation (e.g., Davis 2011, Hovick et al. 2011) which can provide more cover from predators and weather exposure in the vulnerable early postfledging phase (Jones and Bock 2005). Finally, we predicted that fledglings of both species would be sensitive to within-season temperature and precipitation, as fledgling grassland birds can be affected by climate conditions (e.g., Adams et al. 2006) and likely have a limited ability to thermoregulate relative to adults. Our study is the first to concurrently examine all three vital rates in breeding populations of these species and provides insight into the effects of habitat conditions across life-history stages.

Study area
We conducted our study at two mixed-grass prairie sites ( Fig. 1 . We conducted research activities on two pastures at each study site (4 total pastures); pastures ranged in size from 128-177 ha (x̅ = 150.5, SD = 17.6). Private lands support 85% of remaining grassland habitat in North America (North American Bird Conservation Initiative 2013) and much of this land is directly tied to ranching in mixed-grass prairie regions (Lipsey and Naugle 2017). Therefore, in choosing site where grazing occurs, we sought to use plots that were reflective of grasslands in the Northern Great Plains (NGP). We selected pastures based on abundance of singing male focal species, and feasibility of land access. Both pastures at our North Dakota study site were located on the Little Missouri National Grassland, managed by the U.S. Forest Service, while pastures at our Montana site were located both on private and U. S. Bureau of Land Management properties. Pastures at both sites consisted of flat to moderately rolling hills, small seasonal wetlands, and sparse to patchy shrub-cover (e.g., Symphoricarpos occidentalis and Artemisia sp.). Vegetation at sites was predominantly a mixture of non-native, cool-season grasses and native, warm-season, mixed-grass prairie species as well as a diversity of primarily native forbs. Cool-season grasses included Kentucky Bluegrass (Poa pratensis), Crested Wheatgrass (Agropyron cristatum) and Smooth Brome (Bromus inermis); these three species comprised the majority of non-native plant cover at all sites. Native grasses included Blue Grama (Bouteloua gracilis), Western Wheatgrass (Agropyron smithii), Needle-and-thread Grass (Hesperostipa comata), Prairie June Grass (Koeleria macrantha), and Green Needlegrass (Nassella viridula). Our North Dakota and Montana sites differed substantially in average cover of non-native species (x̅ ND = 49%, SD = 26; x̅ MT = 10%, SD = 19). In addition to our avian focal species, other species of grassland songbird commonly breeding at sites included, Chestnut-collared Longspur (Calcarius ornatus), Sprague's Pipit (Anthus spragueii), Western Meadowlark (Sturnella neglecta), Savannah Sparrow (Passerculus sandwichensis), and Bobolink (Dolichonyx oryzivorus). Domestic cattle (Bos taurus) grazed intermittently at variable stocking rates on all pastures throughout the duration of the study and although we did not record grazing intensity yearly, we systematically surveyed all pastures twice a season in 2016-2018 to characterize the structure and cover conditions at our sites (TableA1.1; see Data collection). In 2018, one of our Montana study pastures had to be partially shifted because of a severe unscheduled burn the previous fall.

Nesting success
We located and monitored nests (Table 1) of Baird's and Grasshopper Sparrows from May 24th to August 7th, 2015-2018. To maximize focal species sample-size, we located nests using rope-dragging and systematic walking (Winter et al. 2003), behavioral observation (Martin and Geupel 1993), and opportunistic discovery while traversing plots during other research activities (e.g., telemetry, vegetation surveys). Because birds are most active early in the morning, and there is evidence that grassland birds are most likely to be on the nest from sunrise to 0900 (see Holmes 2012, Kirkham andDavis 2013), we conducted our nest searching efforts primarily during this period. We did not conduct rope-drag surveys when temperatures were less than 10° C, when grass was excessively wet (enough to soak boots and clothing), or during active precipitation. Upon locating nests, we marked the location with two pin-flags, one 5 m to the north, and the other 5 m to the south; we did not place any markers immediately near nest openings to avoid attracting predators directly to nests. We visited nests every three days, with longer intervals occasionally (e.g., five days) in cases of poor weather. We used Global Positioning System (GPS; etrex 10 and etrex 20; Garmin Ltd., Olathe, KS, USA) units and handdrawn microhabitat maps with compass bearings to relocate nest entrances during checks. We also took care to avoid creating trails by varying our approaches to the nest (Major 1990). To increase our ability to correctly assign nest fates, we visited nests more frequently (1-2 days) when they were near fledging age. Human visitation frequency does not appear to affect the nesting success of ground-nesting grassland birds (O'Grady et al. 1996, Pietz et al. 2012, Border et al. 2018), thus we felt this approach could improve the accuracy of our nesting success data with little risk. At each visit, we recorded and photographed nest contents, examined nests for evidence of predators or brood parasitism by Brown-headed Cowbirds (Molothrus ater), and aged nestlings based on physiological development (Jongsomjit et al. 2007, Ruth andKitting 2018). To address nest fate uncertainty (Manolis et al. 2000), we implemented 15-minute fledgling searches during the final visit to every potentially successful nest (mean fledging age for the species; Davis 2003 andJones et al. 2010). During these observation periods, we searched the nest area for fledglings and observed adults for feeding activity. These data helped us to confirm whether nests had fledged or been depredated during the late nestling stage when exact age at termination could not be known (Manolis et al. 2000). We found that the movements of both Baird's and Grasshopper Sparrow fledglings at our sites was extremely limited during this period ( Fig. A2.1-2) reducing the possibility of conflating fledglings from neighboring territories (Streby and Andersen 2013a). We considered nests successful if they fledged at least one host young, failed if they did not, and unknown in cases where nestlings were potentially old enough to fledge, but we found no clear evidence of fledging or depredation (see Manolis et al. 2000).

Adult survival.
We captured and tracked adult male Baird's and Grasshopper Sparrows from May 14th to August 10th, 2015-2017 (Table 1). We chose to focus our capture efforts on males to avoid disrupting nesting females, as pilot data indicated that capture of nesting females triggered nest abandonment. To capture birds, we identified territorial individuals, set up 1-2 mobile mist-nets where we observed males perching and singing, and played audio-lure at the net location to draw birds in for capture. We conducted target-netting activities from sunrise to 1000 MST, but we did not target-net in temperatures less than 10° C, during precipitation, or winds greater than ~24 km/h. Following capture, we outfitted individuals with very high frequency (VHF) radio-transmitters (PicoPip Ag379; 0.42 g, ~30-40-day battery-life; Lotek Wireless, Seattle, WA, USA) using an elastic leg-loop harness for attachment (Rappole and Tipton 1991). We weighed all birds prior to transmitter attachment (x̅ = 18.1 g, SD = 1.1) to ensure that the transmitter weighed no more than 4% of the individual's total body mass. We also fitted all captured birds with a United States Geological Survey (USGS) federal aluminum band, and one or more plastic color bands, and collected standard morphometrics. We briefly observed birds for signs of adverse reaction to harness fit following release (e.g., hampered flight ability).
We tracked all tagged birds daily (sunrise to 1500 MST) using a receiver (Biotrack VHF Receiver; Lotek Wireless) and 3-or 5element folding Yagi antennae (Advanced Telemetry Systems, Isanti, MN, USA). During tracking, we took care to vary the time of day we tracked a given individual bird, and to circle a suspected bird location before flushing the bird, to avoid chasing it through the grass and influencing the recorded location. We recorded a single location per individual each day. We did not track birds in extreme weather conditions (e.g., lightening, heavy rain, excessive winds). After locating birds, we recorded their status (alive or dead) and their location using hand-held GPS units. We tracked birds until transmitters failed (failure was usually characterized by weak or intermittent signal prior to disappearance) or until signal suddenly disappeared (before the expected life of the transmitter). In the latter case, we searched the perimeter of the plot using a 5-element antenna for extended range, but we were rarely successful in relocating birds once we had lost the signal. It is likely these individuals emigrated, given the nomadic lifehistory of grassland birds. For example, breeding Grasshopper Sparrows are known disperse up to 9 km within season (Williams and Boyle 2018). However, we considered the fates of missing individuals unknown, as emigration could not be confirmed. We considered birds dead only when we found physical evidence of mortality (e.g., dead bird, blood on harness, transmitter with chew marks, pile of feathers, transmitter found buried, etc.). In cases where we found an undamaged, clean transmitter and harness, we assumed it had fallen off due to loose fit, or that the bird had been able to remove it by pulling at the elastic with its bill. In some instances, individuals previously fitted with a transmitter were recaptured confirming that they were alive and had shed their harnesses. We were not able to assess the presence of negative transmitter effects in our study (Barron et al. 2010), as we did not have a control group for comparison. However, we seldom observed unusual behavior in tagged birds during our study. In most cases, individuals displayed good flight ability and resumed territorial behavior after release (e.g., perching and singing, interacting with other males).

Post-fledging survival
We monitored post-fledging survival in Baird's and Grasshopper Sparrows from June 15th to August 7th, 2016-2018 (Table 1). We tagged 1-2 nestlings per nest with VHF transmitters (PicoPip Ag337; 0.29 g, ~20-30-day battery-life; Lotek Wireless), depending on the brood-size and individual weights, when nestlings were 7-8 days old. We chose this age to help minimize risk of force fledging (Berkeley et al. 2007), as our focal species typically fledge at 9-10 days old in the NGP (Davis 2003, Jones et al. 2010. We used the same leg-loop harness attachment as for adults, but with a slightly looser harness fit to allow for growth. We also weighed and applied federal metal bands to all nestlings prior to transmitter attachment (x̅ = 13.6 g, SD = 1.3). We ensured that transmitters weighed no more than 3.0% of the individual's total body. In cases where greater than two individuals weighed enough to receive transmitters in a nest, we selected individuals to tag by placing all nestlings of adequate weight in a bird bag and drawing at random. After the birds fledged, we tracked tagged birds daily using the same protocols described for adults but with additional caution to avoid trampling young fledglings.
We considered fledglings dead if they were less than 10 days old when the signal was lost following Hovick et al. (2011), as our movement data ( Fig. A2.1-2) suggested that fledglings were not capable of traveling a distance beyond our antenna range (~300-400 m) from their previous location in a single day at this age. This pattern of limited movement by young fledglings is widely reported in grassland songbirds (Davis and Fisher 2009, Hovick et al. 2011, Small et al. 2015, Young et al. 2019. We therefore assumed they had been carried off by a predator or that the transmitter was damaged during a depredation event. We were not able to compare survival between tagged and untagged birds because of the difficulty involved in tracking untagged fledglings in a grassland environment. However, research on Henslow's Sparrow (Centronyx henslowii) fledglings suggests the negative effects of transmitters are negligible (Young et al. 2019). We also acknowledge that in tracking only individuals of sufficient weight to allow tagging under our permits, our survival estimates may have been biased. Finally, we could not confidently partition causes of mortality in fledglings, as individuals that died of exposure or starvation could also show signs of predation from scavenging activities.

Vegetation surveys and weather data
We characterized the vegetation composition at each nest location using a Daubenmire frame (Daubenmire 1959) centered around the nest to estimate percent cover of live grass, dead grass, forbs, bare ground, and litter. We also estimated the mean height of vegetation contained within the frame. Finally, we estimated visual obstruction (VOR) at the nest site from four cardinal directions using a Robel pole at a distance of 4 m and a height of 1 m (Robel et al. 1970). We selected this suite of variables because they are the most commonly assessed and influential vegetation covariates in studies of grassland birds (Fisher and Davis 2010). We conducted vegetation surveys at nest sites within three days of nest termination to best represent conditions at the nest while it was active. We also collected vegetation data at fledgling bird telemetry locations in 2017-2018. Because estimated bird locations were sometimes less precise than nest locations, we estimated percent cover and mean vegetation height over a 5-m radius area around the bird location. The cover types we estimated at these locations included the same categories recorded at nest sites, in addition to an estimate of non-native species cover.

Nesting success
To analyze nesting success, we estimated daily survival (DSR) of nests using logistic exposure (Shaffer 2004); we performed the analysis in Program R (R Version 3.6.1; R Core Team 2020) using the lme4 package (Bates et al. 2015). Logistic exposure is a modified form of logistic regression that accounts for length of exposure in survival analyses while also providing a flexible modeling approach in which multiple continuous and categorical variables can be included (Shaffer 2004). We included nests with unknown fates in analysis, but truncated intervals to the date of last known activity (Manolis et al. 2000). We combined failures from depredation, abandonment, cowbird parasitism, weather, and unexplained causes in analysis. In cases of nest abandonment, we assigned failure either to the interval following the last known date of nest activity, or the interval in which an event suspected of causing the abandonment occurred (e.g., partial predation, extreme weather event). We excluded nests we thought had been abandoned because of research activity, such as capture of an adult, from analysis (n = 34); this occurred primarily in 2015 and 2016 when we experimented with tracking adult females. We extrapolated DSR estimates produced with logistic exposure over a period approximating the length of a complete nesting cycle in each species (21 days for Baird's Sparrow and Grasshopper Sparrow; Davis 2003, Jones et al. 2010 to calculate the cumulative probability of nesting success. We modeled the influence of environmental variables on nesting success with covariates for year (2015-2018), date (days from May 1 st ), daily precipitation, average precipitation over the previous week to observation date (hereafter weekly precipitation), daily minimum temperature, daily maximum temperature, and daily average temperature. We intended that weekly precipitation would capture effects of drought when present, while daily precipitation could reflect severe weather events as well. We averaged continuous, temporally explicit variables over intervals to best represent conditions during the exposure period being evaluated (see Shaffer 2004). We did not include lagged bio-year climate effects (i.e., effect of previous year's climate conditions) because we expected that any strong annual trends would be subsumed by year effects. Finally, we modeled a full suite of nest-site vegetation covariates including cover of live grass, dead grass, forbs, litter, bare ground, total grass, total vegetation, mean vegetation height and VOR.
We centered and standardized all continuous variables by subtracting the mean and dividing by the standard deviation to make parameter estimates directly comparable (Schielzeth 2010). We tested for correlation among all continuous variables using Pearson's correlation coefficient and did not include variables in the same model if r > 0.4. We selected this threshold because it has been demonstrated that even low levels of collinearity can bias parameter estimates in multiple regression (Graham 2003). We constructed our models in a two-step process. First, to select among correlated variables to include in full analysis, and to compare linear and polynomial terms for continuous variables, we compared global models (consisting of non-correlated variables) with interchanged terms of interest (e.g., one global model with only linear date, one with linear and quadratic date). We tested all continuous predictors for both linear and quadratic effects. We used an information theoretic approach to model selection, selecting the model with the lowest AIC c value (Burnham and Anderson 2002). Because site and year were conflated in our study (e.g., in 2015, data were only collected at the North Dakota site, and only in one pasture), we could not include both variables in global models, thus we chose to include year because it was of more interest to our study questions. However, we compared a univariate model with site to a constant survival model and found that, in both species, site models were indistinguishable from the null model by absolute distance (Baird's Sparrow: ΔAIC c = 2.00; Grasshopper Sparrow: ΔAIC c = 1.02). Following variable selection, we compared all model subsets (including a null constant survival model) evaluating the strength of each model using ΔAIC c and AIC c weight (w i ). We used a conservative threshold of ΔAIC c < 2 to determine top models in our comparisons (Burnham et al. 2011). To address modelselection uncertainty, we used the MuMin package (Barton 2020) to average predictions from a full suite of models (with 95 % confidence intervals). However, because model-averaging has been criticized as a valid means of parameter estimation (Cade 2015), we report parameter estimates (β) with standard errors (SE) and 85% confidence intervals (see Arnold 2010) from only the top model that each variable appeared in. To avoid the use of variable weights and model-averaged confidence intervals in detecting uninformative parameters (Cade 2015), we considered all variables appearing in the top model set to be informative, but only if 85% confidence intervals did not overlap zero in any top model. This procedure allowed us to filter out uninformative parameters in competitive models (Arnold 2010). We only discuss results and report estimates for variables that fit these criteria. Additionally, if the constant survival model appeared in the top model set, we took this as a lack of evidence for any effect of variables modeled for that species. Because the unit of analysis is the number of intervals in logistic exposure, not the number of nests (Shaffer 2004), we report interval sample-size in all modeling tables.

Adult and post-fledging survival
We used the same logistic exposure and information theoretic modeling approach described for nesting success to evaluate adult and post-fledging survival. Logistic exposure can be used to estimate survival similarly to nesting success (see Streby and Andersen 2013b), and in our case, it had several advantages. First, because the unit of analyses is the exposure interval and not the individual, we were able to include data from individuals with unknown fates in analysis by truncating encounter histories. Second, for the same reason, inability to disentangle permanent emigration from mortality in individuals did not directly influence our estimates (Lebreton et al. 1992). Third, we were not restricted to fixed re-sighting intervals and therefore could accommodate cases where individuals were tagged or re-sighted asynchronously. For adults, we extrapolated DSR estimates produced with logistic exposure over a 90-day period to estimate the cumulative probability of an individual surviving the breeding season under given conditions. For fledglings, we used a 20-day period; the exact age of independence is not known for either of our focal species, but a study of a similar species, Henslow's Sparrow, reported independence between 19-21 days post-fledge (Young et al. 2019). Additionally, a 20-day monitoring period was consistent both with the battery life of fledgling transmitters, and our observations of individuals in the field by this age (e.g., capable of strong flight, no longer observed being fed by parents).
We modeled the same temporal and climate variables described for nesting success for both adults and fledglings. For fledglings, we also modeled the effect of days post-fledge, and vegetation cover. Because vegetation data were only collected at fledgling bird locations in 2017-2018, we conducted two separate analyses of post-fledging survival to make use of the full sample, but also examine effects of vegetation. The first analysis included the full dataset and included all variables described above. The second analysis included only the 2017-2018 data and modeled all these variables in addition to a suite of vegetation variables including the previously described cover types (see Vegetation surveys and weather data), extent of non-native vegetation, and vegetation height. Because we did not collect vegetation data at adult bird locations in all years, we were not able to include these variables in the adult analysis. As with nesting success, we could not include variables for year and site in the same global models. However, we evaluated the effect of site in preliminary analysis and found that the site model was not

Nesting success
We located and monitored nests of 152 Baird's Sparrows and 209 Grasshopper Sparrows. Combining data across sites and years, peak nest-initiation (back calculated from hatch date and nestling ages) for Baird's Sparrow (Fig. A3.1A

Baird's Sparrow
Modeled DSR and cumulative survival were 95% (CI: 89 to 98) and 37% (CI: 14 to 61) respectively under average conditions (covariates set to mean values). Success did not vary statistically among years ( Fig. 2A). Covariates appearing in the top model set included year, nest stage, date, daily precipitation, visual obstruction reading (VOR), live grass cover, forb cover, bare ground cover, and litter cover (Table 2). However, only date (β = -0.43, SE = 0.15, CI: -0.66 to -0.20) and VOR (β = 0.28, SE = 0.16, CI: 0.05 to 0.52) had confidence intervals that never overlapped zero across top models. DSR for Baird's Sparrow declined as the season progressed (Fig.3A) and date had the strongest influence on nesting success in the species, appearing in all top models (Table 2) and having the largest effect size. Additional variation was explained by VOR, which had a moderate positive effect on DSR (Fig. 3B).

Grasshopper Sparrow
Modeled DSR and cumulative nesting success were 91% (CI: 87 to 95) and 16% (CI: 5 to 33), respectively, and showed little annual variation ( Fig. 2A). Nesting success in Grasshopper Sparrow was not influenced by any of the variables we modeled, as the constant survival model (ΔAIC c = 1.27) was in the top model set (Table 3) and 85% confidence intervals overlapped zero for all variables modeled.

Adult survival
Over the course of the study, we tagged and monitored 167 adult male Baird's Sparrows and 149 adult male Grasshopper Sparrows. Of the individuals tagged, 49% of Baird's Sparrows and 27% of Grasshopper Sparrows had known fates, while the remaining individuals either left the study area, had a transmitter fail, or appear to have shed their transmitter. Although we did not formally quantify mortality sources, we observed individuals that appeared to have died from weather exposure or illness with no visible injuries (n = 6), those that had lacerations or had been plucked and torn apart by predators (n = 12), and in several instances we noted that transmitters were found near or inside ground squirrel (Spermophilina sp.) burrows (n = 3). Additionally, in two cases, transmitters were found in and near what were likely Short-eared Owl (Asio flameus) pellets.

Modelled survival estimates
Logistic exposure analysis of survival produced DSR estimates of 99% (CI: 99 to 99) and 99% (CI: 98 to 99) for male Baird's and Grasshopper Sparrows under average conditions. Cumulative breeding season survival estimates for the two species (Fig. 2B) were 78% (CI: 51 to 91) and 74% (CI: 36 to 92) respectively. Survival was invariant among years, and none of the climate and time of season covariates we modeled explained variation in survival; the constant survival model was the top model for both species (Table 4).

Post-fledging survival
We tagged and monitored the survival of 94 fledgling Baird's Sparrows and 62 fledgling Grasshopper Sparrows over the duration of our study. Of these individuals, 95% of Baird's Sparrows and 71% of Grasshopper Sparrows had known fates. The remaining birds were greater than 10 days old when their signals were lost and potentially had the flight capability to leave study pastures, but because they could also have been carried off by predators, their fates were unknown. As with adults, we did not formally attempt to determine the cause of mortality for fledglings. However, we observed fledglings of both species that had likely been depredated or scavenged (n = 40). In these cases, bodies were found injured and dismembered, or transmitters were recovered with blood, feathers, or chew marks. We also recorded instances where transmitters were found buried or in ground squirrel burrows (n = 8), and in two cases, transmitters were tracked after being consumed by Plains Garter Snakes (Thamnophis radix). We also noted many cases where no sign of predation could be found, and mortality was likely associated with starvation or exposure (n = 23).

Grasshopper Sparrow
DSR and cumulative survival for Grasshopper Sparrow fledglings under average conditions were and 97% (CI: 94 to 99) and 55% (CI: 27 to 76) respectively, and survival was consistent among years (Fig. 2C). As in Baird's Sparrow, variables for year, precipitation, and temperature appeared in top models (Table 5), but only age was influential (Fig. 4B) and survival increased with days post-fledge (β = 2.74, SE = 0.65, CI: 1.89 to 3.77). Analysis of vegetation cover and survival (Table 6)

DISCUSSION
Examining multiple vital rates in populations of management interest can provide a broader view of species demographics and more accurately inform management activities. Our analysis of nesting success, breeding-season adult survival, and post-fledging survival in Baird's and Grasshopper Sparrows in the Northern Great Plains (NGP) revealed differences in response to environmental conditions among vital rates, as well as between species, highlighting the value of multi-stage monitoring in understanding species-habitat relationships. We found that mixed-grass prairie specialist Baird's Sparrow was generally more responsive to habitat conditions than the more generalist Grasshopper Sparrow. Baird's Sparrow responded to several environmental covariates in nesting success (time of season, vertical vegetation structure) and post-fledging survival (vegetation height, non-native vegetation cover) while Grasshopper Sparrow was only responsive to one environmental covariate (dead grass cover) across all rates estimated. Adult survival of both species was consistent and high relative to existing overwintering survival estimates for the two species, suggesting that adult survival during other parts of the annual cycle may be of greater importance. Finally, cumulative post-fledging survival of Baird's Sparrow was low relative to theoretical demographic thresholds (Cox et al. 2014) indicating that this vulnerable lifehistory stage is potentially important for the species.

Drivers of vital rates
We found that relationships between survival and environmental conditions varied by vital rate and species. Time of season was the strongest predictor of nesting success for Baird's Sparrow, which declined as the season progressed, an effect previously reported in the species (Davis 2005, Lusk and Koper 2013, Davis et al. 2016b) and in grassland birds generally (Zimmerman 1984, Grant et al. 2005, Davis et al. 2006, Grant and Shaffer 2012. This pattern may correspond to increases in predator abundance Age = days post-fledge, Dead = % dead vegetation, Exotic = % exotic cover, Forb = % forb cover , Grass = % grass cover, Hght = vegetation height, Liter = % litter cover, Live = % live vegetation, Prcp = daily precipitation, Temp(mn) = min temp., Year = study year following early season reproduction (Grant et al. 2005) or increased activity of exothermic predators coupled with warming conditions and longer days through the progression of the breeding season (Burhans et al. 2002). Nests of grassland songbirds in the NGP are subject to predation pressure from a diverse array of nest predators including small mammals, large mammals, snakes, and avian predators (Pietz et al. 2012, Bernath-Plaisted and. Thus, earlier nest initiation on the breeding grounds may be adaptive to increase nesting success by avoiding predation pressure later in the season (Grant et al. 2005). Alternatively, this effect could have been driven by higher quality of early-arriving individuals (Wheelwright andSchultz 1994, Verhulst andNilsson 2008). Nesting success of Baird's Sparrow was also partially driven by a positive effect of VOR (visual obstruction reading) matching the species' preference for both habitat and nest-site selection (Dieni and Jones 2003, Davis 2005, Lipsey and Naugle 2017 and further supporting the hypothesis that increased cover can provide concealment from nest predators (Winter 1999, Fondell and Ball 2004, Davis 2005, Klug et al. 2010. Survival in fledglings of both species was strongly driven by age, with most mortalities concentrated within the first five days postfledge. This pattern is common in fledgling songbirds (Cox et al. 2014) and specifically in grassland birds (Kershner et al. 2004, Fisher and Davis 2011, Hovick et al. 2011, Young et al. 2019. A likely explanation for this pattern is increased mobility with age; fledgling birds may be better able to evade predators as they develop sustained flight. Age may also act as a proxy for body condition, which can also be an important predictor of survival in fledgling grassland birds (Adams et al. 2006, Suedkamp Wells et al. 2007, Jones et al. 2017. Thus, because ethical and permitting limitations prevented us from placing transmitters on individuals below a certain weight threshold, it is possible that our rates overestimate post-fledging survival to some extent. Conversely, it is also possible that transmitters themselves reduced fledgling survival creating a negative bias (Barron et al. 2010) though there is some evidence that transmitters have little effect on fledgling grassland birds (Young et al. 2019).
Post-fledging survival in our study was also influenced by vegetation structure and composition. In Baird's Sparrow, survival increased with greater vegetation height. This finding is consistent with other studies of post-fledging survival in grassland bird species and supports the hypothesis that vegetation cover may provide a refuge from predators during a critical lifehistory period of high-mortality (Jones and Bock 2005, Berkeley et al. 2007, Small et al. 2015, Jones et al. 2017. Predation is typically identified as the largest source of mortality in studies of post-fledging survival in grassland songbirds (Adams et al. 2001, Kershner et al. 2004, Suedkamp Wells et al. 2007, Hovick et al. 2011. Though shading and shelter from severe weather are also plausible benefits of taller vegetation, these mechanisms seem less likely in our study as we found no effect of climate conditions on the post-fledging survival of either species. Post-fledging survival in Grasshopper Sparrow responded negatively to increased cover of dead grass, though the mechanism underlying this effect is less clear. It is possible that stiff dead grass impeded movement of fledglings, or that it provided fewer food resources relative to live vegetation. Finally, we also found that survival of fledgling Baird's Sparrow was negatively affected by cover of non-native plant species. Non-native plants have been linked to reduced food availability in grassland systems (Flanders et al. 2006, Hickman et al. 2006) and although fledglings were still being fed by parents during the monitoring period in our study, a reduction in local food resources could have hampered provisioning efficiency of adults. Negative effects of non-native vegetation on vital rates have been reported for several grassland-obligate songbirds in the Great Plains (Lloyd and Martin 2005, Fisher and Davis 2011, Ludlow et al. 2015, Davis et al. 2016b, but more mechanistic studies are needed to understand and mitigate these apparent deleterious effects.

Species differences
Generally, we found many more relationships between vital rates and the environmental covariates tested for Baird's Sparrow than for Grasshopper Sparrow. The only influential covariates we found for Grasshopper Sparrow affected fledgling survival, where survival declined with dead grass cover and increased with fledgling age. There was no effect of time of season on the nesting success of Grasshopper Sparrow (as there was for Baird's), a surprising finding as nests of both Baird's and Grasshopper Sparrows are likely preyed upon by a similar predator community. Given the similarity in nest structure and habitat preference between the two species, it is possible that heightened predation at Grasshopper Sparrow nests was the result of behavioral differences between the two species. We observed that Grasshopper Sparrows were often more conspicuous at the nest site, perching and chipping when observers approached the nest area, and thus perhaps more liklely to attract predators. Conversly, Baird's Sparrows were often quiet or not visible during nest checks, particularly early in the nesting cycle. Once again, unlike Baird's Sparrow, success of Grasshopper Sparrow nests was also unaffected by vegetation height. One explanation for this discrepancy could be that, in the NGP, although general habitat preferences of Baird's Sparrow and Grasshopper Sparrow are similar (Lipsey and Naugle 2017), Baird's Sparrow has been shown to select more strongly for increased vegetation height at the nest site than Grasshopper Sparrow (Dieni and Jones 2003). We observed this difference in nest-site selection in our study as well, though the difference was marginal (Guido 2020).
Interestingly, lack of response to vegetation characteristics in Grasshopper Sparrow relative to Baird's Sparrow was also apparent in the post-fledging phase. We did not observe a negative effect of non-native cover in survival of fledgling Grasshopper Sparrows, as we did in Baird's Sparrow. It is possible that feeding niches of the two species are subtly different, affecting their ability to exploit non-native food resources. For example, a study of the overwintering diet in Baird's Sparrow and Grasshopper Sparrow found that the two species had different seed preferences and food handling times as a result of differing bill morphologies (Titulaer et al. 2018). More broadly, inconsistent response of vital rates to environmental conditions between the two species could be explained in part by divergent life-history strategies. Baird's Sparrow exhibits more of the common specialist characteristics of the two species, occupying a smaller geographic range ( Fig. 1; BirdLife International 2016), less varied habitat-use (Green et al. 2002, Vickery et al. 1996, and lower population numbers (Partners in Flight 2020). It may be that Baird's Sparrows are simply more sensitive to environmental conditions than Grasshopper Sparrows as a consequence of inherent life history differences. Specialists species are highly successful when operating within ideal habitat conditions, but they have a limited ability adapt to sub-optimal or marginal environments (Correll et al. 2019). For example, it is possible that Baird's sparrow was less successful at foraging in habitat that is partially degraded by non-native plant species than Grasshopper Sparrow, resulting in higher fledging mortality for Baird's sparrow. Conversely, Baird's sparrow could be better adapted to select highly-concealed nest sites in mixed-grass prairie than Grasshopper Sparrow, which uses a shorter vegetation structure in many parts of its range (Vickery 1996), perhaps explaining the lack of positive response to vertical structure in Grasshopper sparrow. Regardless of mechanism, inconsistent effects of vegetation structure on vital rates are frequently reported in mixed-grass prairie songbird communities (Davis 2005, Koper and Schmiegelow 2007, Kerns et al. 2010, Lusk and Koper 2013. This variation in reponse to structure among species reflects the unique microhabitat needs of grassland songbirds and serves as a reminder that grassland songbird species can be poor management surrogates for one another (Davis 2005, Derner et al. 2009, Lipsey and Naugle 2017.

Vital rates across life-stages
Nesting success, the most familiar and commonly measured vital rate during the breeding season, fell within established ranges in the Great Plains for Baird's Sparrow (17-43%; Davis and Sealy 1998, Davis 2003, Jones et al. 2010, Lusk and Koper 2013, Davis et al. 2016b, and Grasshopper Sparrow (14-53%; Berthelsen and Smith 1995, Jones et al. 2010, Hovick et al. 2012, Davis et al. 2016b. It is worth noting, however, that our estimate of Grasshopper Sparrow nesting success (16%) fell on the extreme low end of reported estimates for the species, and yet success in that species was unresponsive to all environmental covariates we analyzed. Thus, additional exploration into the causes of low nesting success for this species in the NGP may be of value. In contrast to nesting success, few studies have isolated adult survival on the breeding grounds in North American songbirds, making our results unique but also difficult to contextualize. Several mark-recapture studies examining annual survival in grassland songbirds species such as Bobolink, Savannah Sparrow, and Florida Grasshopper Sparrow report annual apparent survival typically between 40-60% (Perkins and Vickery 2001, Fletcher et al. 2006, Perlut et al. 2008a, 2008b. However, our results are not directly comparable to these estimates, as our rates are confined to the breeding season period. Additionally, telemetry studies are better able to disentangle mortality from detectability than mark-recapture studies, which often suffer low recapture rates driven by poor site fidelity in grassland birds (Balent and Norment 2003, Fletcher et al. 2006. It is also important to note we were only able to track male birds, which may have biased our estimates. Male survivorship is sometime higher than female in grassland birds, and in songbirds generally (Perlut et al. 2008a, Low et al. 2010. This difference is likely a consequence of mortality incurred by females while on the nest. Nonetheless, our adult survival estimates for the breeding season were notably higher and more consistent than adult survival rates estimated for these species during other parts of their annual cycle (range 4-32%; Macías-Duarte et al. 2017 where winter survival for both species was influenced by environmental factors including precipitation and vegetative cover. Increasingly, there is evidence that substantial adult mortality may occur on the wintering grounds or during migration for some species (Hostetler et al. 2015). Studies of other migratory songbirds breeding in North America have shown that adult mortality is often highest during the migration period Holmes 2002, Rushing et al. 2017). Adult survival during the non-breeding periods of the annual cycle may therefore be of equal or greater importance to overall population growth than breeding season survival for our focal species.
Post-fledging survival has also been little studied in our focal species, but a growing number of studies examining survival of grassland songbird species during this critical life-history phase have been conducted over the past two decades. Our estimate of post-fledging survival in Baird's Sparrow (25%) was similar to rates reported for closely related species such as Grasshopper Sparrow (21%; Hovick et al. 2011), Henslow's Sparrow (25%; Young et al. 2019), and Savannah Sparrow (21-35%; van Vliet et al. 2020), though our own estimate of post-fledging survival in Grasshopper Sparrow was much higher (55%). However, survival of fledgling Baird's Sparrow in our study was also lower than existing estimates for many other grassland obligates such as Dickcissel (Spiza americana; 56%; Suedkamp Wells et al. 2007), Sprague's Pipit (29%; Fisher and Davis 2011), Lark Bunting (27-37%; Adams et al. 2006), and Eastern meadowlark (Sturnella magna; 63-69%; Kershner et al. 2004, Suedkamp Wells et al. 2007). Further, survival estimates for Baird's Sparrow fell well below the 40% threshold theoretically necessary to maintain populations without unrealistically high survival during other demographic stages (Cox et al. 2014). Therefore, we suggest that post-fledging survival may be a management consideration for Baird's Sparrow and should be examined more frequently in conjunction with nesting success.
Monitoring multiple vital rates in avian populations of interest is critical not only for management purposes (Fletcher et al. 2006, Perlut et al. 2008b, van Vliet et al. 2020 but also for the creation of accurate population models (Streby and Andersen 2011). While informally comparing vital rates across life-history stages may be the first step in assessing overall population limitation for species, evaluating the relative impact of these seasonal vital rates on population trajectory through full annual cycle population models is necessary to fully understand limitation (Hostetler et al. 2015). Increasingly, avian researchers have shown the importance of assessing demographics across the full spatial and temporal life-cycles of species, and the potential for carryover and interaction between these stages (Latta et al. 2016, Rushing et al. 2017. Integrated population models (IPMs) are a powerful tool for identifying limiting life-stages and geographies, and making accurate population predictions with respect to broad changes in land use and climate (Ahrestani et al. 2017, Zhao et al. 2019. Importantly, for habitat limited species like grassland birds, such models can also be used to identify critical habitat areas (Grand et al. 2019). Although one appealing aspect of IPM's is the ability to estimate latent parameters for which data are lacking, such as survival during the migration period (Ahrestani et al. 2017), such models require large amounts of data from various life-stages. Therefore, accurate data describing fundamental vital rates across multiple life-history phases of species are a critical component of both management oriented and analytical conservation efforts.

CONCLUSION
While demographic monitoring of songbirds has traditionally focused on single vital rates (e.g., nesting success), more comprehensive approaches may yield greater insight into population dynamics (Rushing et al. 2017, Wilson et al. 2018). Our study is the first to simultaneously monitor nesting success, breeding-season adult survival, and post-fledging survival for two declining grassland songbird species in the NGP. We found that Baird's Sparrow demographics were more responsive to environmental conditions than in Grasshopper Sparrow, for which habitat covariates had little impact on vital rates. Both nesting success and post-fledging survival in Baird's Sparrow increased with vertical vegetation structure, and post-fledging survival in Baird's Sparrow was also negatively influenced by nonnative plant cover. Adult survival on the breeding grounds for both species was high and invariant relative to overwintering survival in the same species, indicating that breeding season adult survival in these populations may be of less conservation significance. By contrast, post-fledging survival in Baird's Sparrow was low, suggesting that management of juveniles may be a priority for further research and monitoring efforts. We suggest the combination of these data with other datasets from the nonbreeding grounds in full annual cycle models to formally compare seasonal vital rates with the ultimate goal of creating conservation goals inclusive of the entire life cycle for these declining species.
Responses to this article can be read online at: https://www.ace-eco.org/issues/responses.php/1875