Moose Habitat Selection and Fitness Consequences During Two Critical Winter Tick Life Stages in Vermont, United States†

The moose (Alces alces) is a charismatic species in decline across much of their southern distribution in North America. In the northeastern United States, much of the reduction has been attributed to winter tick (Dermacentor albipictus) infestations. Winter ticks are fairly immobile throughout all life stages, and therefore their distribution patterns at any given time are shaped largely by the occurrence of moose across the landscape during the peak of two critical time periods: fall questing (when ticks latch onto moose) and spring drop-off (when engorged female ticks detach from moose). We used recent land cover and lidar data within a dynamic occupancy modeling framework to estimate first-order habitat selection (use vs. non-use) of female moose (n = 74) during the tick questing and drop-off periods. Patch extinction and colonization rates between the fall questing and spring drop-off periods were strongly influenced by habitat and elevation, but these effects were diminished during the fall questing period when moose were more active across the landscape. From the fall questing period to the spring drop-off period, patches where colonization was high and extinction was low had higher proportions of young (shrub/forage) mixed forest at higher elevations. Further, we evaluated the fitness consequences of habitat selection by adult females during the fall questing period, when females and their calves acquire ticks. We compared Resource Selection Functions (RSF) for five females that successfully reared a calf to age 1 with five females whose calves perished due to ticks. Adult female moose whose offspring perished selected habitats in the fall that spatially coincided with areas of high occupancy probability during the spring tick drop-off period. In contrast, adult female moose whose offspring survived selected areas where the probability of occupancy during the spring drop-off was low; at present, natural selection may favor female adults who do not select the same habitats in fall as in spring. Our model coefficients and mapped results define “hotspots” that are likely encouraging the deleterious effects of the tick-moose cycle. These findings fill knowledge gaps about moose habitat selection that may improve the effectiveness of management aimed at reversing declining population trends.

The moose (Alces alces) is a charismatic species in decline across much of their southern distribution in North America. In the northeastern United States, much of the reduction has been attributed to winter tick (Dermacentor albipictus) infestations. Winter ticks are fairly immobile throughout all life stages, and therefore their distribution patterns at any given time are shaped largely by the occurrence of moose across the landscape during the peak of two critical time periods: fall questing (when ticks latch onto moose) and spring drop-off (when engorged female ticks detach from moose). We used recent land cover and lidar data within a dynamic occupancy modeling framework to estimate first-order habitat selection (use vs. non-use) of female moose (n = 74) during the tick questing and drop-off periods. Patch extinction and colonization rates between the fall questing and spring drop-off periods were strongly influenced by habitat and elevation, but these effects were diminished during the fall questing period when moose were more active across the landscape. From the fall questing period to the spring drop-off period, patches where colonization was high and extinction was low had higher proportions of young (shrub/forage) mixed forest at higher elevations. Further, we evaluated the fitness consequences of habitat selection by adult females during the fall questing period, when females and their calves acquire ticks. We compared Resource Selection Functions (RSF) for five females that successfully reared a calf to age 1 with five females whose calves perished due to ticks. Adult female moose whose offspring perished selected habitats in the fall that spatially coincided with areas of high occupancy probability during the spring tick drop-off period. In contrast, adult female moose whose offspring survived selected areas where the probability of occupancy during the spring drop-off was low; at present, natural selection may favor female adults who do not select the same habitats

INTRODUCTION
The moose (Alces alces) is a charismatic species that has been in recent decline across much of their southern range in North America (Jones et al., 2019;DeBow et al., in press). In the northeastern United States (US), the reduction has been caused primarily by winter tick (Dermacentor albipictus) infestations (Samuel, 2004;Musante et al., 2007;Jones et al., 2019). Unlike other ungulates such as the white-tailed deer (Odocoileus virginianus), moose have not evolved behavioral strategies (e.g., programmed grooming) to counteract tick parasitism (but see Addison et al., 2019). Consequently, the number of ticks on an individual moose can be staggering, with reports as high as 90,000 ticks (Jones et al., 2019), which can lead to epizootics (>50% calf mortality). Research suggests that heavy tick infestations are a driving force for high mortality and low birth rates, and heavily infested calves just under 1 year old often perish in mid-April due to complications of anemia (Samuel, 2004;Jones et al., 2019;DeBow, 2020;DeBow et al., in press;Rosenblatt et al., in press).
Winter ticks are a native, one-host parasite found on a variety of mammals in the northeastern US, but commonly carry out their life cycle on moose (Samuel, 2004). Winter ticks climb and cluster on early successional vegetation to seek a host from mid-September to the first permanent snowfall, a period known as "questing" (McPherson et al., 2000;Samuel, 2004). The questing period coincides with the moose breeding period (i.e., the mid-September rut) when moose tend to be most active on the landscape. Once on a moose, winter ticks go through two life stages, nymph and adult, taking blood meals and molting between each stage (Samuel, 2004). The blood meals in latefall and early winter are not as significant as the adult blood meal in spring, when an adult female tick will typically consume 1.70-2.55 g of blood (Addison et al., 1998). After engorgement, female ticks drop from the moose into the leaf litter, where they lay their eggs and the cycle continues.
Winter ticks tend to be fairly immobile throughout all life stages, and therefore their distribution patterns at any given time are shaped largely by the patterns of occurrence of moose across the landscape (McPherson et al., 2000). As a highly mobile species, moose have the capacity to traverse great distances (e.g., a radio-collared individual in the study moved approximately 150 km within a few weeks). When studying moose occurrence, the distinction between habitat selection and habitat use is important. Habitat selection is the process an animal takes to actually choose or select a particular habitat given a range of options, while habitat use is the result of that choice, regardless of the available options (Johnson, 1980;Beyer et al., 2010). The selection of resources is thought to occur in a hierarchical manner, from a more coarse-scale examination of where a species selects their geographic range (first-order selection), to the home range of an individual within their geographical range (secondorder), and finally to more fine-scale examination of the habitat components within home ranges (third-order) (Johnson, 1980). Meyer and Thuiller (2006) indicate that first-order selection could be divided further into population and metapopulation levels to account for selection of areas used by populations within the geographic range.
Resource selection functions (RSFs) are a primary tool for understanding habitat selection (Manly et al., 2002;Boyce, 2006). By comparing random locations at a given scale with those actually used by the species, RSFs yield the probability of use of a resource unit (or a result that is proportional to it). Importantly, the scale at which habitat selection is evaluated shapes inferences regarding key habitat components (Boyce, 2006). For example, a study in central Ontario, Canada indicated that at a regional scale (first-order selection), moose selected areas of canopy disturbance (Forbes and Theberge, 1993), while a study in Sweden examined habitat selection at the home range scale (second-order) and found that moose selected clear-cuts and young and medium-aged forests (Cederlund and Okarma, 1988). Finally, within home ranges (third-order), moose in central Finland preferred non-pine dominated habitats and mature forests, and avoided human settlements (Nikula et al., 2004).
Understanding which habitats adult female moose select during the tick fall questing and spring drop-off periods may be important for the successful management of moose, especially during epizootic years. Adult female moose have high value for moose management as they drive recruitment of the population and significantly contribute to habitat use, home ranges, and movements of calves (Cederlund and Okarma, 1988). Previous studies have examined habitat use or selection of female moose across different seasons (e.g., Cederlund et al., 1987;Ball et al., 2001;Terry, 2015;Blouin et al., in press). However, few have examined habitat selection during the two relatively short, critical winter tick life periods. One important study in New Hampshire and Maine, United States, examined third order habitat selection during both the questing and drop-off tick periods and found that moose selected 4-16-year-old forest openings relative to their availability within their home ranges regardless of season (Healy et al., 2018). This suggests that moose key into the same habitats during both questing and drop-off periods, thus these habitats promote the winter-tick life cycle and their association with moose (Samuel, 2004;Healy et al., 2018). Forest openings provide important forage for moose, potentially attracting higher densities of moose and exacerbating the impact of winter ticks on the moose population as a whole. Although Healy et al. (2018) succinctly identified third order habitat selection within individual home ranges during these critical times, it is important to note that third order selection is conditional on both first and second order habitat selection (across individuals), neither of which have been examined in the northeastern US to date. Analysis of habitat selection at these scales may provide important insights into habitats that are unused or avoided during the questing and drop-off winter tick periods, and provide important context for Healy et al. (2018)'s findings.
Moreover, an understanding of how an individual's habitat selection decisions ultimately affect their fitness is lacking. Fitness can be defined as a measure of an individual's ability, relative to others, to produce viable offspring (Stearns, 1992). During epizootic years, the habitat selected by an adult female moose during the fall tick questing period may positively or negatively affect fitness and the fate of offspring (Pekins, 2020). The probability that a calf will survive to their first birthday (May) may depend on areas traversed during their first fall questing period (September). Ideally, adult female moose will select optimal habitats (i.e., those with abundant forage) during the fall questing period while limiting exposure to damaging tick loads.
Understanding moose habitat selection during the critical fall tick questing and spring tick drop-off periods and the resulting fitness consequences is a priority need for wildlife managers in the northeastern US. In recent years, moose have been declining across their southern extent, including Vermont and much of the New England region (Jones et al., 2019;DeBow et al., in press). At this time, direct management through regulated hunter harvest provides one potential means of breaking the tick-moose cycle and allowing recovery (VFWD, 2020), while indirect management through the creation or alternation of habitats may provide additional tools for breaking the cycle. As such, examination of habitat selection during the winter tick questing and drop-off period and how that selection relates to fitness may provide insights for the management of an iconic species in decline.
Our objectives were to: (1) Investigate first order habitat selection of adult female moose during each period using a multi-season occupancy framework (MacKenzie et al., 2003); and (2) Evaluate the fitness consequences of selection by comparing second-order fall questing resource selection for adult female moose that successfully reared a calf to age 1 versus those that did not.

Study Area
The study area was in northeast Vermont, United States (Figure 1, mean latitude = 44.77 • ; mean longitude = -71.74 • ; 1,738 km 2 ). This area contains the highest density of moose in the state and represents an important region that connects populations in New Hampshire and Maine with those in New York and southern Quebec (Pearman-Gillman et al., 2020). The study area encompassed two state Wildlife Management Units (E1 and E2) in Essex County and consisted of extensive bogs and softwood swamps, and young, intermediate, and mature forest stands. Elevation ranged from 200 to 1000 m. In higher elevations (>800 m), species composition was dominated by red spruce (Picea rubens) and balsam fir (Abies balsamea), while intermediate elevations consisted primarily of maple (Acer spp.), birch (Betula spp.), and beech (Fagus grandifolia). Lowland swamp and bog areas were dominated by balsam fir, red spruce, black spruce (P. mariana), poplar (Populus spp.), paper birch (Betula papyrifera), and alder (Alnus spp.). Vermont experiences four distinct seasons, including summer (June-August), fall (September-November), winter (December-February), and spring (March-May). Between 2017 and 2019, winter temperatures averaged -3.81 • C, with an average snowfall ranging between 220 and 250 cm (NCDC, 2019). During summer months, temperatures averaged 18 • C, with average precipitation ranging between 100 and 110 cm (NCDC, 2019).

Radio-Collaring
A helicopter was used to capture and radio-collar 126 moose between 2017 and 2019 (see DeBow et al., in press). Collaring occurred in January of each year under snow conditions that increased visibility and limited potential injury to individual moose. In 2017, 60 moose (30 ∼ 8-month-old calves, 30 adult female moose) were fixed with Survey Globalstar V7.1 GPS/VHF collars (VECTRONIC Aerospace GmbH, Berlin, Germany) and ear tags. In 2018 and 2019, 36 moose (30 calves and 6 adult female moose) and 30 moose (30 calves) were similarly collared, respectively. Moose were captured throughout the study area, although more moose were fixed with radio-collars in the northern portion. Collars were programmed to send a Global Positioning System (GPS) coordinate every 13-h. We used GPS Plus X software to manage collars and receive location data. All capture, handling, and radio-collaring procedures were reviewed and approved by the University of Vermont Institutional Animal Care and Use Committee (IACUC protocol #17-035).

Objective 1: First Order Habitat Selection of Adult Female Moose With Multi-Season Occupancy Analysis
We used dynamic occupancy models to assess first-order habitat selection of female moose during the winter tick questing and winter tick drop-off periods (Design I; Manly et al., 2002). These models allow inference about the occurrence of moose across a collection of spatially identified "sites" or patches, and how patterns of occurrence change through time (MacKenzie et al., 2003;Martin et al., 2009). Dynamic occupancy models can be used to assess first-order habitat selection when patches (resource units) are not considered independent (MacKenzie, 2006). In the original formulation (MacKenzie et al., 2003), the initial probability that a species occupies a sample unit during a defined sampling period defines occupancy (psi, ψ). Thereafter, two additional parameters drive changes in occupancy probability over time: colonization (gamma, γ) and extinction (epsilon, ε) (MacKenzie et al., 2003). A fourth parameter (p) models the detection probability of a species at a site, given presence. Importantly, the parameters of dynamic occupancy models can be modeled as a function of covariates; in our case, those included habitat and landscape conditions related to moose and winter ticks (MacKenzie et al., 2017). FIGURE 1 | The location of the study area in northeastern Vermont, United States (1,738 km 2 ), encompassing Essex County (Wildlife Management Unit E). GPS radio-collars were attached to moose (Alces alces) (n = 126) in the area and monitored from 2017 to 2019. The study area was bounded by the Canada-U.S. border to the north, New Hampshire to the east, VT-Route 2 to the south, and VT-Route 114 to the west. Red triangles indicate the capture location for female moose (n = 74), while circles indicate where male moose were captured (n = 52). Blue triangles show the capture location for adult female moose (n = 10) included in the habitat selection and fitness consequences analysis (i.e., objective 2). The map on the right shows the study area in relation to other northeastern states (ME, Maine; NH, New Hampshire; and MA, Massachusetts).
Global Positioning System collar data were filtered into the two primary periods for all female moose captured during the study. For each year, the questing period included the peak 50 days of questing (26 September -15 November), while the drop-off period spanned from (16 March -5 May) (Healy et al., 2018). Within each period, data were further subset into two survey periods of equal duration, which allowed us to address detection probability. After considering only female calves and adults within our defined periods, analyses of occupancy were based on 12,210 GPS locations of 74 moose (Figure 1).
We imposed a grid of 1693 patches, each 1 km 2 , over our study area using the aggregate() function in the R package, raster (Hijmans and Van Etten, 2012). Encounter histories based on collared moose locations were created for each patch over seven total primary periods (drop-off 2017, questing 2017, dropoff 2018, questing 2018, drop-off 2019, questing 2019, drop-off 2020), each with two secondary survey periods (14 total survey periods). We used the rasterize() function in the raster package to count the number of GPS coordinates that occurred within a grid cell patch in each survey period. The counts were then converted to a binary encounter history (i.e., if the count within a patch was greater than 0, it was assigned a 1, otherwise a 0). For example, a patch with an encounter history of 11 00 10 10 10 11 01 indicated the patch was initially occupied by moose during the first primary period (11, drop-off period of 2017), where moose were detected in both the first and second survey periods. In the second primary period (00, the questing period of 2017), the patch either went locally extinct, or failed to go extinct but moose were undetected in both survey periods. The patch was occupied in the third primary period (10, tick drop-off period of 2018), due to one of two conditions: if the patch was extinct in the previous period, it was re-colonized. If the patch was occupied in the previous primary period but simply undetected, it again failed to go locally extinct. In either case, the species was detected in the first survey but missed in the second survey. The remaining primary periods (10 10 11 01) suggest that the patch remained occupied, but moose were undetected on some surveys.
In the context of first order habitat selection of adult female moose (use of resource units), the assumptions of the multiseason occupancy model were: (1) patches were closed with respect to changes in occupancy within each primary period, and (2) there were no false positive detections (MacKenzie et al., 2004). When the first assumption was not met, we assumed changes in occupancy were random and we re-defined occupancy as "site-use." In the context of occupancy modeling utilizing GPS points across individuals, two additional assumptions were made. First, by including both adult and calf female GPS points, we assumed that a patch occupied by a collared female calf was also occupied by an un-collared, nearby female adult. Further, we assumed autocorrelation among patches was negligible due to the number of sites, the number of moose collared, and their ability to traverse great distances.
We identified habitat variables of potential importance to moose occupancy (Table 1). These covariates included 30 m 2 National Land Cover Data (USGS, 2016a) spatial layers describing land cover types (deciduous forest, coniferous forest, mixed forest, and wetland), 30 m 2 layers describing terrain characteristics (elevation and slope) (VCGI, 2002), and 10 m 2 lidar (light detection and ranging) variables (USGS, 2016b) characterizing forest age structure (shrub, forage, cover, and canopy). We used the aggregate() function in the raster package to rescale each raster to match the 1 km 2 patch resolution for occupancy analysis, where each patch value provided the average of the underlying cells (e.g., Figure 2).
We then performed a principal component analysis (PCA) using the princomp() function in the R stats package (R Core Team, 2018) to account for the complex and often correlated relationships between habitat composition rasters (mean evergreen, deciduous, mixed, and wetland), forest structure rasters (mean canopy, cover, forage, shrub), and mean elevation. Each principal component was influenced differently by the habitat variables; the top three principal components were included as habitat variables in the dynamic occupancy model set.
We used the package, RPresence (Hines, 2006) to evaluate changes in occupancy patterns across patches and through time [(occMod()) function, type = "do.1"]. We created a model set of 24 models (including the intercept model), representing alternative hypotheses that explain changing occupancy patterns through time ( Table 2). Each model consisted of four submodels that specified the initial occupancy, detection, colonization, and extinction parameters and related them to resource attributes at a spatial scale of 1 km 2 . Across all models, the formula predicting initial occupancy (ψ) remained constant (∼PC1 + PC2 + PC3), as did the formula predicting detection (p) across the study area (∼PC1 + PC2 + northing + slope + season). The PCA and slope variables in the detection submodel were informed by a GPS collar bias study in our study area (Blouin et al., in press), where GPS transmission rates were found to be a function of habitat variables and slope. The northing variable was included in the detection submodel because more animals were collared in the northern portion of our study area than in the south. The complexity of the gamma (γ) and epsilon (ε) formulas varied across the model set, from one variable examining the influence of year, to more complex formulas examining the interaction between habitat variables and the season (winter tick questing vs. drop-off), which was of specific interest to test hypotheses that colonization and extinction rates depended on critical winter tick periods. Year was included in most models to account for new animals that were collared each January, potentially affecting both extinction and colonization rates ( Table 2).
Models were compared using Akaike information criterion (AIC) (Burnham and Anderson, 2002) to determine a single, best model to predict multi-season occupancy. We used the RPresence (Hines, 2006) function, createAicTable() to rank models and compute model likelihoods. All inferences were made based on the top performing model in the model set.

Objective 2: Fitness Implications of Second-Order Resource Selection Models (RSFs) for Adult Female Moose That Successfully Reared a Calf to Age 1 Versus Those That Did Not
We developed both logistic regression and negative binomial RSFs to assess second-order habitat selection of adult female moose during the fall questing period, and whether these decisions influenced the survival of calves the following spring (Design II from Manly et al., 2002; see also Thomas and Taylor, 2006). As all individuals were collared in January, this required that we targeted uncollared calves with collared adult female moose, allowing inference of where that calf was with its mother during the fall questing period, when acquiring winter ticks. We were successful in collaring 10 calves, providing a small dataset of adult female-calf pairs in which the fate of the calf was known (0 = calf died prior to first birthday, 1 = calf survived to first birthday) as well as the GPS locations of their mothers during the fall questing period. Two RSFs were estimated: one for adult females that successfully reared offspring to age 1, and one for those whose calves perished prior to age 1. We used the same base rasters and approach as for Objective 1 ( Table 1) but added an additional 30 m 2 raster that provided the probability of GPS transmission rate, given each pixel's habitat, elevation, and slope (Blouin et al., in press). Unlike the first objective where analysis was at the 1 km 2 resolution, we gridded the rasters at 200 m 2 resolution for a more fine-scale assessment of selection using the aggregate() function in the raster package. This scale was selected because it was fine enough to discriminate habitat conditions among patches, yet coarse enough to accommodate multiple GPS points within a patch (Boyce, 2006). To account for all habitat variables while avoiding issues associated with modeling highly correlated variables, we performed a Principal Component Analysis with the princomp() function in the R package, stats (R Core Team, 2018). Each component represented a different combination of the nine input variables (mean evergreen, deciduous, mixed, wetland, canopy, cover, forage, shrub, and elevation). The four top components accounting for the majority of the total variation of data were selected to represent habitat composition, structure and terrain variables.
For each adult female, we counted the number of GPS points that fell within each 200 m 2 raster patch during the fall questing period (September 26 -November 15); unused patches were assigned a 0. For each adult female, we conducted 1000 bootstrap trials wherein each trial's dataset consisted of 10,000 randomly sampled unused patches and all patches that included any GPS locations. For the logistic regression RSFs, data were collapsed into used/unused binary classifications. The logistic regression model for each female was "used ∼ PC1 + PC2 + PC3 + PC4" with a logistic link function, while the negative binomial models had the form "counts ∼ PC1 + PC2 + PC3 + PC4" with a log link function. Models were run with the glm() and glm.nb() functions, respectively (R Core Team, 2018). In both models, data points were weighted by the inverse of the patch's GPS transmission rate, which provided the probability that a GPS signal would be detected given the patch's habitat conditions, slope, and elevation (Blouin et al., in press). For each model and trial, model coefficients were stored, and population-level RSFs were based on aggregated coefficients across trials.

Objective 1: First-Order Habitat Selection of Adult Female Moose With Multi-Season Occupancy
The 1693 (1 km 2 ) patches covering our study area exhibited differences in forest composition, forest structure, and elevation (Figure 2). The encounter histories for each patch identified the Models include combinations of variables defined by three PCA variables (PC1, PC2, PC3) that accounted for 72% of total variance. Each formula describes how colonization probability and extinction probability were modeled as a function of the different resources within each PCA variable. Seasons related to tick behavior and included 'questing' in the fall and 'drop-off' in the spring. The study occurred across a three-year period (2017-2019). An "x" marks the combination of each variable included in a model.
locations of collared female moose during the primary questing and drop-off periods (Figure 3). More patches were occupied during the tick questing periods than drop-off periods: mean patches occupied by moose during the questing period was 468 (naïve occupancy = 0.28), while mean patches occupied during the drop-off period was 409 (naïve occupancy = 0.24). This pattern likely reflects moose movement behavior during the fall breeding season. The PCA resulted in three principal components describing 72% of the total variance of the habitat variables (Table 3). PC1 was positively influenced by forage/shrub/cover structure and wetland/evergreen composition (i.e., as PC1 increases, forage/shrub/cover and wetland/evergreen forest increase), and negatively influenced by deciduous forest, and to a lesser extent, the canopy structure and elevation. PC2 was positively influenced by mixed/evergreen forests that was more mature in structure and negatively influenced by young (forage/shrub) deciduous forests and wetlands. Finally, PC3 was positively influenced by canopy structure and evergreen/wetland forests and negatively influenced by mixed forest and elevation ( Table 3 and Figure 4). The spatial variation of the PC scores and their meanings can be traced back to the original variables (Figure 2).
For female moose, the top model for predicting shifting occupancy patterns across winter tick questing and drop-off periods was model 23 ( Table 4). As with all candidate models, initial occupancy for the top model was a function of PC1, PC2, and PC3 and detection was a function of PC1, PC2, northing, slope, and season. Extinction and colonization probabilities for the top model were a function of the interaction between PC1 and season, PC2 and season, PC3 and season, and year. The top model carried virtually all the weight in the model set; no other models were competitive with AIC scores exceeding 18.0 ( Table 4). The intercept model, in which parameters were modeled as constants without covariates, had a AIC score of 223.27 and a model weight of 0 ( Table 4).
The top model estimated 28 parameters ( Table 5). For initial occupancy, PC2 had a strong positive effect and PC1 had a slight positive effect, meaning that during the initial primary period (2017 fall drop-off period), the probability of occupancy increased as the proportion of structurally diverse (canopy, cover, forage, and shrub), mixed and evergreen forest increased (Figures 4, 5A). For detection, moose had a higher probability of being detected in the northern part of our study area and where the slope angle was relatively low (Table 4 and Figure 5B). Probability of detection was also higher during the questing period and in patches comprised of more young deciduous forest and less mature evergreen or mixed forest (i.e., patches with lower PC2 scores, Figure 5B).
Extinction rates of any given season describe the probability that patches that were occupied in the previous season become locally extinct and varied between questing and drop-off seasons (Figure 5). From the fall questing period to the spring drop-off period, patch extinction rates were negatively related to PC2 and positively related to PC3 (Figure 5C, top panel); patches that were occupied during fall questing but went extinct during the spring drop-off period were characterized by higher proportions of mature (canopy) evergreen forests and wetland habitats at lower elevations (Figures 4, 6A). Patches on the landscape that failed to go extinct (i.e., remained occupied) between the fall questing and spring drop-off periods were characterized by higher proportions of young (shrub/forage) mixed or deciduous forests at higher elevations (Figures 4, 6A). In contrast, from the spring drop-off to the fall questing period, patch extinction rates were weakly driven by PC3 only, and overall extinction probabilities were low ( Figure 5C, bottom panel, Figure 6B). The few patches that were occupied in spring but went locally extinct in the fall were characterized by higher proportions of mature (canopy) evergreen forests and wetland habitats at lower elevations ( Figure 6B). Patches on the landscape that failed to go extinct (or remained occupied) during the fall questing period were characterized by higher proportions of young (shrub/forage) deciduous or mixed forests across the elevation gradient.
Colonization rates of any given season describe the probability that patches that were unoccupied in the previous season became locally colonized and varied between questing and drop-off seasons (Figure 5). From the fall questing period to the spring drop-off period, patch colonization rates were negatively related to both PC1 and PC3 (Figure 5D, top panel); patches that were unoccupied during fall questing but were colonized during the spring drop-off period had higher proportions of young (shrub/forage) mixed forest at greater elevations ( Figure 6C). Patches that failed to be colonized (remained unoccupied) between the fall questing and spring drop-off periods were characterized by higher The cumulative proportion indicates that the first three components are explaining 72% of the total variance of the data for objective 1 and 78% for the first four components in objective 2. proportions of mature evergreen forests and wetlands at lower elevations. In contrast, from the spring drop-off to the fall questing period patch colonization rates were largely driven by PC3 (Figure 5D, lower panel). Patches that had higher probabilities of being colonized from the spring drop-off to the fall questing period had higher proportions of mature (canopy) deciduous and evergreen forests and wetland habitats at lower elevations. Lastly, patches that failed to be colonized (remained unoccupied) from the spring drop-off to the fall questing period were characterized by higher proportions early Model names indicate structure for colonization(gam) and extinction(eps). Structure for occupancy and detection was "psi(PC1 + PC2 + PC3)" and "p(PC1 + PC2 + northing + slope + season)." In model names, "S" is season, "Y" is year, and "YT" is year-trend. " AIC" is the delta AIC score, "NLL" is -2 * loglikelihood, "K" is the number of parameters, and "Weight" is the model weight.
To summarize the complex dynamic occupancy model results, patch extinction and colonization rates between the fall questing and spring drop-off periods were driven by at least 2 principal components with strong effect sizes, while a single component with smaller effect sizes influenced rates from the spring drop-off to the fall questing period. Further, patches where colonization was high and extinction low from the fall questing period to the spring drop-off period (highest probability of female moose occupancy) had higher proportions of young (shrub/forage) mixed forest at greater elevations (Figures 6A,C). In contrast, patches where colonization was high and extinction low from the spring drop-off to the fall questing period were characterized by higher proportions of both mature (canopy) and young (shrub/forage) deciduous or mixed forests and wetland habitats at low elevations (Figures 6B,D).

Objective 2: Fitness Implications of Second-Order Resource Selection Models (RSFs) for Adult Female Moose That Successfully Reared a Calf to Age 1 Versus Those That Did Not
The 200 m 2 pixel resolution resulted in 42,691 patches for RSF analysis. The PCA resulted in four principal components describing 78% of the total variance of the habitat variables (Table 3). PC1 was positively influenced by younger (shrub/forage) wetland habitats (i.e., as PC1 increases, forage/shrub and wetland forest increase), and negatively influenced by the proportion of mature (canopy lidar classifications) deciduous forest. PC2 was positively influenced by mixed and evergreen forests that were more mature in structure and negatively influenced by early succession (shrub/forage) deciduous forests. PC3 was positively influenced by mature (canopy) evergreen forests and wetland habitats and negatively influenced by mixed forests and elevation. PC4 was positively influenced by higher elevation deciduous and evergreen forests of all age classes and negatively influenced the mixed forest composition (Table 3 and Figure 7).
Both the logistic regression RSFs and negative binomial RSFs indicated strong differences in habitat selection between female moose whose calves survived to their first birthday (∼18 May; n = 5) vs. died (n = 5) ( Table 6 and Figure 8). In general, female moose with calves that survived exhibited much less variation in how they selected PC2 than those with calves that died (Figure 8). Specifically, adult females with surviving calves had a consistent, slightly negative response to PC2, while females with non-surviving calves showed variable responses to PC2. Most noticeably, adult females with surviving calves did not respond or had slightly positive responses to PC3 and had negative responses to PC4, in contrast to adults with non-surviving calves who showed negative responses to PC3 and no response to PC4. These patterns are consistent across both RSF models (logistic and negative binomial) but are more apparent in the negative binomial RSF (Figure 8).
Component loadings indicate that adult female moose whose offspring survived to their first birthday selected habitats during the fall questing period with higher proportions of young (shrub/forage) deciduous forests, but also selected habitats with higher proportions of mature (canopy) evergreen or mixed forests and wetlands at lower elevations (Figures 7, 9). Adult female moose whose offspring perished prior to their first birthday were less selective of young (shrub/forage) deciduous forests, instead selecting habitats with higher proportions of young mixed forests at higher elevations during the fall tick questing period (Figures 7, 9).

DISCUSSION
Winter tick infestations have significantly impacted moose populations in Vermont and the surrounding region (Jones et al., 2019;DeBow et al., in press). Observations of declining populations spurred research on mortality and productivity of moose across the region. Since 2014, >500 adult females and calves have been fixed with radio-collars regionally (Vermont, New Hampshire, Maine, United States) (Jones et al., 2019;Ellingwood et al., 2020;Pekins, 2020;DeBow et al., in press). These studies showed low observed birth rates and reoccurring epizootic mortality due to winter tick infestations (5 of 6 years between 2014 and 2019) (Jones et al., 2019;Ellingwood et al., 2020;Pekins, 2020;DeBow et al., in press). Research suggests there are three main factors typically influencing the severity of winter tick infestations -the density of the host (moose) (Pekins, 2020), favorable conditions for winter tick survival (i.e., shorter winters) (Samuel, 2007), and the overlap of moose habitat use during two critically important winter tick life periods (i.e., questing and drop-off periods) (Healy et al., 2018;Pekins, 2020).
To better understand the winter tick-moose cycle and the influence of seasonal overlap, Healy et al. (2018) examined thirdorder habitat selection of moose during the critical tick fall questing and spring drop-off periods. Results showed that moose selected optimal 4-16-year-old forest openings relative to their availability within their home ranges across both the questing and the drop-off periods. This illustrates that once home ranges are established, optimal habitats are often re-used between the tick questing and tick-drop off periods, promoting the life cycle of the winter tick within established home ranges.
Our study of first and second-order habitat selection, along with the fitness consequences of selection, builds on the findings of Healy et al. (2018), while addressing other factors that may also be influencing the cycle and impacting moose fitness. Firstorder habitat selection analysis at the landscape scale revealed that moose were more widespread in their occupancy during the fall questing period versus the early spring drop-off period, and patches with higher proportions of young (shrub/forage) mixed forest at greater elevations had the highest probability of occupancy during the spring tick drop-off period. As moose evolved to lose considerable body mass during the winter months (particularly more severe ones), it follows that habitat use would be more constrained during the late-winter/early-spring dropoff period compared to fall (Schwartz et al., 1988). Furthermore, significant protein loss due to heavy tick infestations would likely restrict habitat use during the late winter months (Jones et al., 2019;Pekins, 2020;Rosenblatt et al., in press). Due to the poor condition of moose by the drop-off period, moose may be seeking swollen buds of young vegetation during the early spring greenup to regain the energy lost over the winter months (Schwartz and Renecker, 2007;Pekins, 2020).
Second-order habitat selection analyses during the fall questing period showed a clear pattern in habitat selection by female moose whose offspring survived versus those that died. For calves of collared adults that perished, the cause of death was determined to be serous atrophy of fat, which is a common consequence of substantial tick infestations. Adult female moose whose offspring perished selected patches during the questing period that were characterized by higher proportions of young (shrub/forage) mixed forest at higher elevations (areas with the highest probability of occupancy during the spring drop-off period; Figures 6A,C). That is, their questing habitat selection matched the first-order models of selection for the population during the spring drop-off period, similar to Healy et al.'s (2018)  conclusions. In contrast, adult female moose whose offspring survived showed selection patterns characterized by young (shrub/forage) deciduous habitats, while also selecting habitats with higher proportions of mature (canopy) evergreen forests and wetlands at lower elevations, i.e., their second-order habitat selection patterns deviated from the overall first-order patterns described by the multi-season occupancy analysis.
Why are adult female moose with calves that survived selecting areas differently than adult female moose with calves that perished? Habitat selection theory suggests that organisms choose habitats that maximize fitness, and that fitness often declines with increasing density (Fretwell, 1969;Rosenzweig, 1981). Moose may trade patches of more optimal habitats on the landscape for more suboptimal patches based on several factors, including energy, time, risk of injury due to competition, and predation risk (Dussault et al., 2005;Doligez et al., 2008). Alternatively, habitat selection may be driven primarily by cues identifying suitable habitats acquired during early years as a calf (imprinting) or later (learned), that ultimately prove beneficial to their fitness (Cederlund et al., 1987;Doligez et al., 2008). Whatever the reason, natural selection is favoring female adults who selected habitats with higher proportions of mature (canopy) evergreen forests and wetlands at lower elevations during the fall questing period.
Locations where ticks drop off in spring determine where they will quest in the fall. The abundance of questing winter ticks across space, in turn, is a function of moose density but also winter tick survival. Research supports the density hypothesis (Samuel, 2007;Pekins, 2020). With more hosts, the potential for proliferation of parasites increases. As moose share similar needs, certain social cues and/or habitat quality may be driving their habitat selection decisions, which may be creating higher densities in particular areas on the landscape (i.e., hotspots) (Doligez et al., 2008;Blouin et al., in press). The abundance of questing winter ticks may also depend on the ground-level conditions between drop-off and questing (i.e., winter tick survival). Studies of winter ticks suggests that heavier tick loads in the fall may be influenced by the proportion of canopy cover (Samuel, 2007). Winter ticks have a higher chance of survival in more open forested habitats than canopied or closed habitats, due to sunlight and warmer temperatures (Drew and Samuel, 1986;Addison et al., 2016). Thus, moose that select habitats with both canopy structure and adequate forage structure (i.e., adult female moose with calves that survived; Figure 9), while unknowingly avoiding significant tick infestations, may be helping to identify habitats across the study area optimal for surviving ongoing tick epizootics.   Wildlife managers have considered a suite of actions that aim to reduce winter tick parasitism on moose. First, direct management strategies, such as hunter harvest, have been suggested to reduce the density of moose in certain areas where productivity and mortality are severely impacted by winter ticks (VFWD, 2020). Although moose have been in decline, further reduction of high local densities may be a proactive way to lessen the prolonged suffering and impacts from winter tick infestations (Mercer and McLaren, 2002;Ellingwood et al., 2020). Second, more indirect management strategies, such as the conservation, modification, or formation of certain habitat types may interrupt the tick-moose cycle. Forest management may change or manipulate the availability and quality of habitats, ultimately influencing these patterns of use and selection (Schrempp et al., 2019). Lastly, direct approaches to winter tick management are being explored to lessen the impact of ticks on moose fitness. The naturally occurring Metarhizium anisopliae fungus found in the soil has shown detrimental effects on winter tick survival in the laboratory setting (Sullivan et al., 2020).
Our analyses of first and second-order habitat selection, combined with Healy et al. (2018)'s third-order analysis, may inform all three of these potential management approaches. First, an understanding of where the probability of moose occupancy is highest may focus efforts of more direct management strategies, such as hunter harvest to reduce moose and tick densities. For instance, if the preferred method of management is to reduce local moose densities, reductions might be targeted on the areas or regions that have been encouraging the proliferation of ticks (Figures 6A,C, 9). In terms of habitat management strategies, our model coefficients and mapped results define "hotspots" for both moose and winter ticks as shrubby mixed/deciduous forest in higher elevations (Figures 6A,C,  9). Knowledge about these conditions may influence habitat management decisions, including the spatial distribution of such hotspots across the landscape. Further, timber management that encourages understory regeneration while maintaining some canopy structure (i.e., uneven-aged silviculture; habitats selected by adult female moose whose offspring survived), may be beneficial for moose fitness and detrimental to tick survival. In terms of direct winter tick management, information is lacking. Little is known about the ecology, abundance, and management of winter ticks in their natural setting. Managers may use the potential "hotspots" results to focus field research (e.g., Sullivan et al., 2020) that could ultimately help break the negative impacts of the tick-moose cycle.
Several caveats and questions remain. All 10 adult female moose-calf pairs were evaluated within the 2017 or 2018 questing period. During this time, Vermont experienced epizootic or near epizootic-level mortality (>50%; 10-month old calves). Thus, our models of habitat selection investigated a snapshot in time of a host-parasite cycle; model results may differ in non-epizootic years. Moreover, it would be beneficial to know if an adult female's habitat selection patterns vary by year. We were unable to monitor adult female habitat selection and trace their offspring recruitment for more than a single year. Considering habitat selection of radio-collared individuals across time may lead to conclusions about whether decisions leading to decreased fitness by individuals is a recurring event or not. Further, we did not evaluate the composition of forest species that were selected by adult female moose during these time periods, the nutritional landscape (Schrempp et al., 2019), or the vegetation that is selected by winter ticks for questing in the fall. It is known that tree and shrub species vary in their nutritional value and availability for moose throughout the year (Timmermann and McNicol, 1988). Thus, future consideration of the availability and/or quality of forage during the fall questing or spring dropoff periods may be important for the management of both moose and winter ticks.
Additionally, our analyses focused on females and calves, yet males also act as wide-roaming vehicles of winter tick spread on the landscape. Throughout the growth season, male moose have much larger home ranges than females (Blouin et al., in press). More specifically, during the breeding season, mature male moose become territorial and may push younger males to more suboptimal habitats (Schwartz and Renecker, 2007). Additionally, male moose use scent urination to attract female moose and induce ovulation during the breeding period (Miquelle, 1991;Schwartz and Renecker, 2007). Research suggests that adult female moose may aggressively compete for access to bull urine (Miquelle, 1991), which may have important implications on habitat selection of female and their offspring during the questing period when moose are picking up ticks.
Finally, we did not examine winter tick success in terms of habitat selection of moose. One of the major influences on the severity of winter tick infestations is how favorable the conditions are for tick survival. This includes habitat structure and composition, but is also influenced by climatic FIGURE 9 | A comparison of a logistic Resource Selection Function (RSF) and negative binomial Resource Selection Function (RSF) by female moose (Alces alces; n = 10) during the fall winter tick (Dermacentor albipictus) questing period in northeastern Vermont, United States. Patches (200 m 2 ) across the study area show areas female moose selected (green color) or did not select (white) during this time period for moose whose offspring lived versus those that died. factors (i.e., onset of winter during the questing period or summer drought) (Samuel, 2007;Dunfey-Ball, 2017). Factors of a changing climate may ultimately be driving the winter tick-moose cycle by benefiting the survivability and success of ticks, thus consideration of weather patterns in relation to the acquisition of ticks and fate of moose may be beneficial.
Ultimately, the relationship between winter ticks and moose is complex. Perhaps equally complex is the appropriate and proactive management of moose populations in the face of changing climate patterns, encroaching human development, colliding public opinions, and ongoing consequences from increased parasitism. Structured decision making (SDM) approaches (e.g., Robinson et al., 2016;Franklin et al., 2020) can help managers weigh options and calculate trade-offs. Our aim was to provide science-based knowledge that may shed light on the habitats of importance to moose health that may be incorporated into SDM frameworks; in the end, this may lead to a better understanding of the multifaceted winter tick-moose cycle.