Seasonal drought in North America’s sagebrush biome structures dynamic mesic resources for sage‐grouse

Abstract The North American semi‐arid sagebrush, Artemisia spp., biome exhibits considerable climatic complexity driving dynamic spatiotemporal shifts in primary productivity. Greater and Gunnison sage‐grouse, Centrocercus urophasianus and C. minimus, are adapted to patterns of resource intermittence and rely on stable adult survival supplemented by occasional recruitment pulses when climatic conditions are favorable. Predictions of intensifying water scarcity raise concerns over new demographic bottlenecks impacting sage‐grouse populations in drought‐sensitive landscapes. We estimate biome‐wide mesic resource productivity from 1984 to 2016 using remote sensing to identify patterns of food availability influencing selective pressures on sage‐grouse. We linked productivity to abiotic factors to examine effects of seasonal drought across time, space, and land tenure, with findings partitioned along gradients of ecosystem water balance within Great Basin, Rocky Mountains and Great Plains regions. Precipitation was the driver of mesic resource abundance explaining ≥70% of variance in drought‐limited vegetative productivity. Spatiotemporal shifts in mesic abundance were apparent given biome‐wide climatic trends that reduced precipitation below three‐quarters of normal in 20% of years. Drought sensitivity structured grouse populations wherein landscapes with the greatest uncertainty in mesic abundance and distribution supported the fewest grouse. Privately owned lands encompassed 40% of sage‐grouse range, but contained a disproportional 68% of mesic resources. Regional drought sensitivity identified herein acted as ecological minimums to influence differences in landscape carrying capacity across sage‐grouse range. Our model depictions likely reflect a new normal in water scarcity that could compound impacts of demographic bottlenecks in Great Basin and Great Plains. We conclude that long‐term population maintenance depends on a diversity of drought resistant mesic resources that offset climate driven variability in vegetative productivity. We recommend a holistic public–private lands approach to mesic restoration to offset a deepening risk of water scarcity.


| INTRODUC TI ON
Water input is one of the most dynamic determinants of terrestrial productivity in arid and semi-arid regions and is fundamental to biological processes responsible for ecosystem function. Droughtprone ecosystems account for 40% of terrestrial land surface globally and support upwards of 2 billion people (Gilbert, 2011). These regions are defined by relatively low mean annual precipitation rates and excessive evapotranspirative demands that directly influence primary production (Le Houérou, Bingham, & Skerbek, 1988;Noy-Meir, 1973). Annual and intra-annual variability in precipitation is typical with distinct dry seasons and unpredictable prolonged droughts that lead to frequent periods of water scarcity (Schlesinger et al., 1990). Human population growth, combined with shifts in climatic conditions, is likely to increase pressure on these ecosystems and further strain already limiting water resources (MEAB, 2005).
In arid and semi-arid ecosystems, seasonal drought often limits resource availability, serving as an ecological minimum or bottleneck, thereby reducing the abundance and distribution of associated wildlife populations (Bolger, Patten, & Bostock, 2005). Nonlinear landscape response during these periods can generate distinct spatial patterns of vegetative productivity associated with contrasting soil moisture dynamics and floral communities that structure geographic distribution of resources (Knapp & Smith, 2001). Efficiencies in arid lands conservation rely on an ability to predict landscape response to offset water scarcity, yet our understanding of drought to structure spatiotemporal intermittence of primary production remains limited (Vicente-Serrano et al., 2013). Misalignment of conservation actions can impede ecosystem benefits needed to counteract predicted climatic fluctuations that have the potential to operate as new and powerful ecological constraints to biodiversity in semi-arid landscapes (Maron, McAlpine, & Watson, 2015).
The greater and Gunnison sage-grouse (Centrocercus urophasianus and C. minimus; hereafter "sage-grouse") are gallinaceous birds and indicator species of the semi-arid sagebrush (Artemisia spp.) biome of western North America (Rowland, Wisdom, Suring, & Meinke, 2006;Schroeder et al., 2004). Sage-grouse are emblematic of water scarcity because they exhibit a bottleneck in reproductive cost associated with nutritional stress that aligns with periods of seasonal drought in late summer (Blomberg, Sedinger, Nonne, & Atamian, 2013). Seasonal drying and senescence of herbaceous vegetation induce sage-grouse to seek out few remaining productive sites associated with highvalue foraging habitat for chick growth and survival. Sites include wet meadows, riparian corridors, drought resistant rangelands and irrigated alfalfa, hereafter "mesic resources" (Connelly, Rinkes, & Braun, 2011). Sage-grouse are representative of spatial patterns in droughtinduced ecological minimums, as productive mesic sites in late summer provide an important, but limiting food resource that structure sage-grouse abundance and distribution within broader landscapes (Donnelly, Naugle, Hagen, & Maestas, 2016). Significant temperature increases across western North America in recent decades (Melillo, Richmond, & Yohe, 2014) now threaten availability of mesic resources that have foreshadowed growing concerns of deepening drought and its effect on sage-grouse populations (Gibson, Blomberg, Atamian, & Sedinger, 2017;Guttery et al., 2013).
Restoration and maintenance of mesic resource productivity in sagebrush landscapes must rely on geographically specific conservation strategies adapted to address regional drought sensitivities. The sagebrush biome exhibits considerable complexity in precipitation magnitude, timing, and variation (Rajagopalan & Lall, 1998). Spatially dynamic patterns of resource availability driven by these processes produce disparities in habitat quality and wildlife fitness potential within landscapes (Coates et al., 2016;Pastor, Moen, & Cohen, 1997). Wide-ranging species like sage-grouse may experience regional variability in selective pressures associated with nonlinear landscape response to drought as soil moisture is a primary driver of seasonally important food resources associated with mesic sites (e.g., forbs and macro invertebrates, Wenninger & Inouye, 2008).
Further complicating water scarcity are geographical shifts in proportional public and private land ownership (Donnelly et al., 2016) linked to regionally pervasive land-use practices that can inhibit or promote drought effects.
Changes to large-scale ecosystem dynamics increasingly require more than just field surveys to understand, monitor, and report on their effects (Marvin et al., 2016). To better quantify landscape sensitivity to seasonal drought and its influence in structuring mesic resources for sage-grouse, we leveraged over 15,000 satellite images to produce a spatiotemporal dataset that tracked annual vegetative productivity patterns across the sagebrush biome from 1984 to 2016. Productivity data were combined with land tenure information to evaluate proportional mesic resource abundance by ownership. The new dataset introduces geographic scale and perspective to ecological drought and its relationship to sage-grouse that to this point remain unexplored. Results provide a framework that for the first time links local evidence-based studies to rangewide drought effects influencing demographic constraints in sage-grouse populations (Blomberg, Sedinger, Atamian, & Nonne, 2012;Blomberg, Sedinger, Gibson, Coates, & Casazza, 2014;Gibson et al., 2017;Guttery et al., 2013). Study outcomes deliver new insight to support development of regionally specific conservation strategies necessary to offset drought-induced bottlenecks impacting sage-grouse and other drought sensitive wildlife in sagebrush ecosystems.

K E Y W O R D S
drought, Great Basin, Great Plains, net primary productivity, precipitation, Rocky Mountains, sagebrush biome, sage-grouse 2 | MATERIAL S AND ME THODS

| Study site
We used current sage-grouse range (>74 million ha, Schroeder et al., 2004) to delineate the sagebrush biome and define the study area boundary (Figure 1). Sage-grouse distribution provided a requisite for large and unfragmented sagebrush landscapes as >90% of birds occur in areas of >40% sagebrush land-cover (Knick, Hanser, & Preston, 2013). Sage-grouse range was divided into three regions (Great Basin, Rocky Mountains, and Great Plains) along gradients of ecosystem water balance using climatic differences associated with temporal overlap of seasonal precipitation and evapotranspiration (Sala, Lauenroth, & Golluscio, 1997, Figure 1). Alignment of these processes correlates to differences in soil water availability (Lauenroth, Schlaepfer, & Bradford, 2014). Summer precipitation occurring when temperature and plant growth are high and potential evapotranspiration (PET; Black, 2007) is greatest supports pulse soil water dynamics that minimize deep soil water storage (Sala, Lauenroth, & Parton, 1992). The opposite occurs when overlap is poor and precipitation comes during cold periods when PET is lowest and potential for deep soil water storage is maximized.
We estimated gradients of pulse and deep soil water potential by calculating correlation coefficients of daily precipitation and PET (Schlaepfer, Lauenroth, & Bradford, 2012). Estimates were calculated within water years (from 01 October to 30 September) and fit to a continuous grid (4 × 4 km) using day of year precipitation and day of year PET means for the study period (1984 to 2016) obtained from Gridded Surface Meteorological Dataset (GRIDMET; Abatzoglou, 2013). Results were scaled from negative one to one and interpreted to be more predictive of summer precipitation pulse soil water potential as positive values increased toward one. Negative values decreasing to negative one were considered more predictive of cold season precipitation and deep soil water potential, with values near zero predictive of both pulse and deep soil water characteristics. We averaged coefficient values within North American level III ecoregions (https://www.epa. gov/eco-research/ecoregions) by assigning the mean of intersecting grid cells to each ecoregion. Ecoregions were aggregated geographically using a k-means clustering analysis of their correlation coefficient mean and intersected with current sage-grouse range to form the Great Basin, Rocky Mountains, and Great Plains regions ( Figure 1).
F I G U R E 1 Study area delineated using sage-grouse range as a proxy to define sagebrush biome. The study area was divided into three areas (Great Basin, Rocky Mountains, and Great Plains) along gradients of broader ecosystem water balance as shown Portions of sage-grouse range in Canada (<1% of study area) were excluded from correlation coefficient calculations because continuous estimates of PET were unavailable. Sagebrush steppe landscapes in eastern Washington, United States, were omitted from our analysis due to the area's high rate of fragmentation and agricultural influences that were uncharacteristic of the remainder of sage-grouse range (Shirk, Schroeder, Robb, & Cushman, 2017).

| Modeling drought limited primary production
In semi-arid ecosystems, aboveground NPP is controlled by patterns of soil water availability (Noy-Meir, 1973). We examined this concept by delineating spatiotemporal patterns of mesic resource productivity as a proxy to soil moisture derived from Landsat 4, 5, and 8 satellite imagery. Measurements were based on normalized difference vegetation indices (NDVI) which quantify photosynthetic activity and correlate closely to fluctuations in NPP (Box, Holben, & Kalb, 1989;Pettorelli et al., 2005). Satellite images (n = 5,180) were used to conduct annual monitoring during a 33-year span (1984-2011 and 2013-2016) to account for climate driven variation in landscape condition (Loik, Breshears, Lauenroth, & Belnap, 2004). Poor-quality Landsat 7 imagery prevented monitoring in 2012. All images were calibrated across sensors and corrected for atmospheric effects and illumination/viewing geometry (Masek et al., 2006;Vermote, Justice, Claverie, & Franch, 2016). Pixels containing surface anomalies (i.e., cloud, cloud shadow, snow, and water) were removed using the Landsat CFMask band. Images were acquired 15 July to 30 September to coincide with high evapotranspiration demand in sagebrush ecosystems when mesic food resource availability for sage-grouse is restricted (Connelly, Rinkes, et al., 2011)  To delineate spatial patterns of mesic resources, we generated a raster image (30x30 m pixel) of maximum NDVI values selected from 32 years of overlapping Landsat images (n = 15,180) used in our analysis. Individual pixels were representative of the highest primary productivity measured 15 July to 30 September, 1984 to 2016. Results were representative of maximum landscape productivity potential during seasonal drought. We applied an object-based segmentation algorithm to this image using program eCognition 9.2 (https://www.ecognition.com) that clustered pixels into polygons representative of underlying landscape features. The resulting polygon layer was used in our model as a non-uniform sampling grid (hereafter "mesic grid") to estimate availability and spatial distribution of vegetative productivity. Only polygons containing mean and maximum NDVI pixel values indicative of higher primary productivity (≥0.3) were retained. We considered this NDVI value as the threshold separating non-productive sites from productive and drought resilient sites with high value for foraging sage-grouse (Donnelly et al., 2016;Weier & Herring, 2000). NDVI polygon values under this threshold were considered to have near zero probability of seasonal drought resiliency and were continuously unproductive and dry during late summer months.
Datasets were summarized within the mesic grid and polygons containing individual or combined majority forest, woodland, or agricultural crop removed to eliminate sites uncorrelated to sage-grouse habitats (Connelly, Rinkes, et al., 2011). Polygons containing alfalfa were retained due to sage-grouse food resource values associated with this crop (Connelly, Rinkes, et al., 2011).
To reduce bias from mesic resources with low probability of sage-grouse use, we quantified sagebrush cover proximal to individual mesic polygons by applying a neighborhood analysis using the LANDFIRE dataset. Polygons containing a cumulative sum of <10% sagebrush cover within a 1.0 km radius of their boundary were removed. This step eliminated center portions of large mesic complexes associated with irrigated hay and alfalfa fields. This procedure was not applied to the Great Plains due to rangeland composition containing lower shrub densities. The final mesic grid included approximately 600,000 polygons with probability greater than zero of productivity (i.e., NDVI ≥ 0.3) occurring during seasonal drought between 1984 and 2016.
Annual estimates of NPP were calculated by averaging Landsat images into single multispectral images for each of 32 years. Images were used to calculate seamless NDVI surfaces (30 × 30 m pixels) across the study area. Results provided a continuous measure of mean annual vegetative productivity during seasonal drought. The mesic grid was applied to each surface to calculate NDVI polygon means from intersected pixel values using zonal statistical functions.
Final polygon summaries were representative of annual mesic productivity and patterns of ecological minimums that structure high value food resources for sage-grouse.

| Mesic resource types
Mesic resource polygons were categorized by types as "alfalfa," "rangeland," "riparian," and "wet meadow," to examine drought sensitivity and spatiotemporal patterns. Visual interpretation of high resolution imagery, in combination with decision support from ancillary spatial datasets, was used to classify polygons. Alfalfa sites were representative of irrigated cropland. In semi-arid regions of the western United States, alfalfa is grown as a perennial crop that dominates agricultural production (Lindenmayer, Hansen, Brummer, & Pritchett, 2011). The NASS croplands data layer was used to automate identification of alfalfa.
Riparian and wet meadow areas were linked to shallow groundwater systems in floodplains that supported phreatophytic herbaceous and shrubland plant communities. We differentiated wet meadow and riparian sites within floodplains by confining riparian areas to vegetated corridors bordering stream channels. Wet meadows were visually discernible in aerial imagery due to dominance of herbaceous vegetation and land-use practices associated with livestock pasture and hay production. The "rangeland" type was made up of upland herbaceous and shrubland plant communities.
Rangeland sites inclusive to the mesic polygon grid exhibited lower sensitivity to seasonal drought and higher productivity relative to surrounding upland areas.

| Linking mesic productivity to abiotic factors and land tenure
We estimated percent normal precipitation (PNP) annually  across the study area within a continuous 4 × 4 km grid by summing GRIDMET day of year precipitation data for each year and dividing it by the annual study period mean. All calculations were made within water years beginning 1 October and ending

| Sage-grouse
We examined availability of mesic sites to sage-grouse during seasonal drought by summarizing its abundance proximal to known bird distributions. Sage-grouse populations exhibit distinct clustering patterns that concentrate birds within the sagebrush biome (Doherty, Tack, Evans, & Naugle, 2010). Fitting data summaries to known bird distributions assured us that mesic patterns observed were characteristic of landscapes supporting sage-grouse populations. Distributions were estimated by buffering all sage-grouse lek locations (n = 6,304) by 10 km to account for observed spatial relationships among birds and mesic resources during seasonal drought (Donnelly et al., 2016). Only leks identified as active (male attendance >0) in the last 10 years were used. Lek locations are well documented for sage-grouse and considered a reasonable index of bird distribution (Reese & Bowyer, 2007). Areas of Gunnison sage-grouse were incorporated into estimates without spatial modification due to their limited distribution. Buffered lek areas encompassed 36 of 74 million hectares of current sagegrouse range.
Within lek buffers, we calculated mesic resource density and distance between mesic sites relative to mean annual precipitation observed from 1984 to 2016 using only productive sites (NDVI ≥ 0.3).
Distance between mesic sites was calculated by generating euclidean distance surfaces annually and estimating the mean for each year. Linear regression models were fit to individual regions to compare sensitivity of mesic density and distance to changing PNP.
Proportional abundance of sage-grouse populations was estimated by region (Figure 1) using mean maximum male counts from 2015 lek surveys conducted within the project area. Lek surveys have been widely used by resource agencies to monitor trends in sage-grouse populations and are considered an index of abundance (WAFWA, 2015).

| Data processing
All image processing and raster-based analyses, not otherwise noted, were conducted using Google Earth Engine cloud-based geospatial processing platform (Gorelick et al., 2017). Vector-based processing was completed using QGIS (qgis.osgeo.org). Plotting and statistical analyses were completed using the R Base Package (R Core Team, 2015).

| RE SULTS
Patterns of drought-induced productivity were distinct within Great Basin, Rocky Mountains, and Great Plains. Landscape partitioning was representative of regional gradients in ecosystem water balance.
Negatively correlated precipitation and PET confirmed dominance of cold season precipitation and hydrology influenced by high elevation snowpack in the Great Basin (r = −0.65) compared to summer rainfall-driven events in the Great Plains. A low correlation (r = 0.04) in the Rocky Mountains indicated that precipitation was evenly distributed with a slight bias to cold season winter snowpack. Previous studies using ground-based weather observation data (Lauenroth et al., 2014;Schlaepfer et al., 2012) corroborate spatial patterns derived from GRIDMET in this study.
Mesic productivity was 22% and 34% more sensitive to pre-  (Table 1). These patterns were characteristic of elevational zonation associated with wet meadows and rangelands found in the regions (Figures 4 and 5). Least drought sensitive alfalfa occurred lower in watersheds in the Great Basin and Rocky Mountains and were 501 m and 643 m below rangeland sites (Table 1). Higher drought sensitivity in the Great Basin and the Great Plains dramatically altered the density and spatial distribution of productive mesic sites between above and below average precipitation years (Figures 4 and 6).  Within the sagebrush biome, private lands encompassed 40% of sage-grouse range, but contained a disproportional 68% of mesic resources. Private ownership of mesic sites within regions averaged 59% (SD ± 5%), 68% (SD ± 4%), and 77% (SD ± 3%) in the Great Basin, Rocky Mountains, and Great Plains, respectively.
Nearly 90% of available wet meadow sites in all regions were privately owned; most productive rangelands were publicly owned in the Great Basin (66%) and Rocky Mountains (51%; Table 2).
Riparian ownership was split evenly in the Great Basin with public lands containing a considerable proportion of sites in the Rocky Mountains (38%) and Great Plains (23%). Alfalfa ownership was 95%-97% private (Table 2).

| D ISCUSS I ON
Our analyses are the first to show that sage-grouse populations are structured in part by biome-wide gradients in ecosystem water balance that drive drought sensitivity in mesic resources (Figure 7).
Correspondingly, landscapes with the greatest uncertainty in longterm mesic availability supported the fewest birds. Documented shifts in mesic resource abundance across space and time were apparent (Figures 4-5) given climatic trends that reduced precipitation below three-quarters of normal in 20% of years (1984-2016, Figure 3). Topography and seasonality of precipitation are primary mechanisms underpinning variation in density and distribution of mesic resources despite regional similarities in average annual precipitation (Lauenroth et al., 2014). At broad scales, spatial configuration of vegetative productivity is known to structure sage-grouse distributions, with population abundance positively related to mesic availability (Donnelly et al., 2016). The causal F I G U R E 2 Sensitivity of mesic productivity (NDVI) to percent normal precipitation (PNP). Linear models fit to overall productivity mean and individually by mesic type. Blue isobars depict point cloud density where shorter distance between lines infers higher density of points. Horizontal dashed red line identifies threshold used to delineate productive (NDVI ≥ 0.3) and unproductive mesic sites (NDVI < 0.3) F I G U R E 3 Yearly time series  of productive mesic sites (NDVI ≥ 0.3), binned by type and weighted by hectares. Plots paired with annual percent normal precipitation (PNP) trends for the Great Basin (a), Rocky Mountains (b), and the Great Plains (c) linkage between rangewide bird densities and patterns of mesic resources we observed underscores the selective pressure of seasonal drought as an important determinant of landscape carrying capacity. Our regional findings provide the landscape context to previous studies suggesting weather-driven productivity as a key factor influencing sage-grouse survival (Blomberg et al., 2013;Guttery et al., 2013).
Through our findings we speculate that a diversity of mesic resources (i.e., rangelands, riparian, and wet meadows) rather than any single type is what sustains grouse populations over time. During drought, resilient wet meadows maintain adult survival as opposed to years of above average precipitation when resurgent mesic range-

lands (Figures 4-5) maximize food availability and bird recruitment.
Previous studies describe a slow-paced life history of sage-grouse (Connelly, Hagen, & Schroeder, 2011;Moynahan, Lindberg, & Thomas, 2006), recognizing that populations are reliant on relatively stable adult survival supplemented by occasional pulses in recruitment associated with favorable climatic conditions (Blomberg et al., 2012;Guttery et al., 2013). Regional drought sensitivity identified herein ( Figure 2) may compound these relationships and act as ecological minimums to influence demographic performance differently across the species range. For example, in the Great Basin, distances between mesic sites were double in comparison to other regions ( Figure 8) due to nonlinear patterns of intensifying drought that concentrated available mesic resources in wet meadow valley bottoms and high elevation rangelands (Figure 4), effectively extending the distance that young have to move between productive mesic sites to forage. Increased movements compound drought effects and are a factor known to lower brood survival (Gibson et al., 2017).   (Schroeder, Young, & Braun, 1999). Over the last 33 years, this region received half (52%) its annual precipitation in the nesting season period, which in typical years comes as short duration low intensity events. However, long duration intense rainfall events become more prevalent in wet years (Lauenroth & Bradford, 2009) with high mesic abundance ( Figure 3c). Therefore, the right amount of precipitation in March-June can enhance demographic performance (Blomberg et al., 2012); conversely, too much rain can reduce nest success (Moynahan, Lindberg, Rotella, & Thomas, 2007;Smith, 2016).
Sage-grouse reliance on climatic intermittence explaining ≥70% of variance in mesic resource availability raises concerns over intensifying droughts predicted in the western United States (Dettinger, Udall, & Georgakakos, 2015). Our model depictions likely reflect a new normal in water scarcity that could compound impacts of demographic bottlenecks in Great Basin and Great Plains regions (Figures 4b and   6b). The immediate concern for sage-grouse conservation is increased frequency and severity of droughts that influence ecological norms in semi-arid regions (Trenberth, Dai, Rasmussen, & Parsons, 2003).
Further complicating drought are altered fire regimes in some sagebrush ecosystems that may negate positive effects of precipitation and population growth (Coates et al., 2016). In the long-term, climate projections that make systems more sensitive to drought could trigger mismatches in timing between resource availability and wildlife needs.
Sage-grouse have demonstrated some adaptive capacity to mitigate drought within existing climatic envelopes, but it is unlikely these traits portend projected landscape conditions (Gibson et al., 2017).
Our findings highlight the critical role of private lands in sagegrouse conservation efforts in the western United States. As evidenced by 68% private ownership of mesic resources (Table 2).

Significant conservation emphasis targets publicly managed uplands
where birds mostly breed and nest (Doherty et al., 2010). show relative density of productive mesic sites by elevation during wet (a) and dry (b) years F I G U R E 7 Average density of productive mesic sites relative to annual percent normal precipitation (PNP) for the period of study 1984-2016. Densities estimated within 10 km radius of known sage-grouse lek locations (n = 6,304) partitioned by region F I G U R E 8 Average distance between productive mesic sites relative to annual percent normal precipitation (PNP) for the period of study 1984-2016. Distances estimated within 10 km radius of known sage-grouse lek locations (n = 6,304) partitioned by region a holistic approach that includes drought resilient wet meadows (86%-90% private-see Table 2) important to supporting population maintenance during periods of water scarcity. Our data allow for new public-private conservation strategies that account for regional variation in mesic resource ownership to promote cross-boundary protection and restoration actions to help birds survive drought.
Sustainability of mesic resources hinges on maintenance of ecological processes and associated land-use practices that foster drought resiliency (Gillson, Dawson, Jack, & McGeoch, 2013). Already being implemented are innovative techniques that restore degraded riparian systems by increasing mesic resistance to drought and elevating landscape productivity to benefit wildlife and ranching (Silverman et al., 2018). Another emerging solution is targeted removal of invasive conifer in high elevation sagebrush rangelands (Miller, Naugle, Maestas, Hagen, & Hall, 2017) known to increase snow retention and extend availability of soil water longer into the growing season (Kormos et al., 2017). Ecosystem benefits of invasive tree removal extend beyond sage-grouse to include conservation of non-target sagebrush-obligate avifauna (Donnelly et al., 2017), enhancement of big game forage (Stephens, Johnston, Jonas, & Paschke, 2016), and promotion of ecosystem resilience to resistance to invasive species (Miller et al., 2017). To accelerate these efforts, we provide access to our mesic resource data through interactive web applications that allow users to track mesic productivity trends in

ACK N OWLED G M ENTS
We thank the National Resources Conservation Service-led Sage Grouse Initiative, part of their Working Lands for Wildlife program, for helping foster and coordinate this science. We too thank Bureau of Land Management for their vision in partnering on this work in support of cross-boundary conservation. We are grateful for all colleagues and technicians that helped with data development and analysis. Science support from the Intermountain West Joint Venture also helped make this work possible. Views in this manuscript from United States Fish and Wildlife Service authors are their own and do not necessarily represent the views of the agency.

CO N FLI C T O F I NTE R E S T
None declared.

AUTH O R CO NTR I B UTI O N S
This study was conceived and designed by J.P.